Atomic-Level Characterization of the Activation Mechanism of SERCA by Calcium

We have performed molecular dynamics (MD) simulations to elucidate, in atomic detail, the mechanism by which the sarcoplasmic reticulum Ca2+-ATPase (SERCA) is activated by Ca2+. Crystal structures suggest that activation of SERCA occurs when the cytoplasmic head-piece, in an open (E1) conformation stabilized by Ca2+, undergoes a large-scale open-to-closed (E1 to E2) transition that is induced by ATP binding. However, spectroscopic measurements in solution suggest that these structural states (E1 and E2) are not tightly coupled to biochemical states (defined by bound ligands); the closed E2 state predominates even in the absence of ATP, in both the presence and absence of Ca2+. How is this loose coupling consistent with the high efficiency of energy transduction in the Ca2+-ATPase? To provide insight into this question, we performed long (500 ns) all-atom MD simulations starting from the open crystal structure, including a lipid bilayer and water. In both the presence and absence of Ca2+, we observed a large-scale open-to-closed conformational transition within 400 ns, supporting the weak coupling between structural and biochemical states. However, upon closer inspection, it is clear that Ca2+ is necessary and sufficient for SERCA to reach the precise geometrical arrangement necessary for activation of ATP hydrolysis. Contrary to suggestions from crystal structures, but in agreement with solution spectroscopy, the presence of ATP is not required for this activating transition. Principal component analysis showed that Ca2+ reshapes the free energy landscape of SERCA to create a path between the open conformation and the activated closed conformation. Thus the malleability of the free energy landscape is essential for SERCA efficiency, ensuring that ATP hydrolysis is tightly coupled to Ca2+ transport. These results demonstrate the importance of real-time dynamics in the formation of catalytically competent conformations of SERCA, with broad implications for understanding enzymatic catalysis in atomic detail

Modeling Signal Propagation Mechanisms and Ligand-Based Conformational Dynamics of the Hsp90 Molecular Chaperone Full-Length Dimer

Hsp90 is a molecular chaperone essential for protein folding and activation in normal homeostasis and stress response. ATP binding and hydrolysis facilitate Hsp90 conformational changes required for client activation. Hsp90 plays an important role in disease states, particularly in cancer, where chaperoning of the mutated and overexpressed oncoproteins is important for function. Recent studies have illuminated mechanisms related to the chaperone function. However, an atomic resolution view of Hsp90 conformational dynamics, determined by the presence of different binding partners, is critical to define communication pathways between remote residues in different domains intimately affecting the chaperone cycle. Here, we present a computational analysis of signal propagation and long-range communication pathways in Hsp90. We carried out molecular dynamics simulations of the full-length Hsp90 dimer, combined with essential dynamics, correlation analysis, and a signal propagation model. All-atom MD simulations with timescales of 70 ns have been performed for complexes with the natural substrates ATP and ADP and for the unliganded dimer. We elucidate the mechanisms of signal propagation and determine “hot spots” involved in interdomain communication pathways from the nucleotide-binding site to the C-terminal domain interface. A comprehensive computational analysis of the Hsp90 communication pathways and dynamics at atomic resolution has revealed the role of the nucleotide in effecting conformational changes, elucidating the mechanisms of signal propagation. Functionally important residues and secondary structure elements emerge as effective mediators of communication between the nucleotide-binding site and the C-terminal interface. Furthermore, we show that specific interdomain signal propagation pathways may be activated as a function of the ligand. Our results support a “conformational selection model” of the Hsp90 mechanism, whereby the protein may exist in a dynamic equilibrium between different conformational states available on the energy landscape and binding of a specific partner can bias the equilibrium toward functionally relevant complexes

A Coevolutionary Residue Network at the Site of a Functionally Important Conformational Change in a Phosphohexomutase Enzyme Family

Coevolution analyses identify residues that co-vary with each other during evolution, revealing sequence relationships unobservable from traditional multiple sequence alignments. Here we describe a coevolutionary analysis of phosphomannomutase/phosphoglucomutase (PMM/PGM), a widespread and diverse enzyme family involved in carbohydrate biosynthesis. Mutual information and graph theory were utilized to identify a network of highly connected residues with high significance. An examination of the most tightly connected regions of the coevolutionary network reveals that most of the involved residues are localized near an interdomain interface of this enzyme, known to be the site of a functionally important conformational change. The roles of four interface residues found in this network were examined via site-directed mutagenesis and kinetic characterization. For three of these residues, mutation to alanine reduces enzyme specificity to ∼10% or less of wild-type, while the other has ∼45% activity of wild-type enzyme. An additional mutant of an interface residue that is not densely connected in the coevolutionary network was also characterized, and shows no change in activity relative to wild-type enzyme. The results of these studies are interpreted in the context of structural and functional data on PMM/PGM. Together, they demonstrate that a network of coevolving residues links the highly conserved active site with the interdomain conformational change necessary for the multi-step catalytic reaction. This work adds to our understanding of the functional roles of coevolving residue networks, and has implications for the definition of catalytically important residues

The Mechanism of Ubiquitination in the Cullin-RING E3 Ligase Machinery: Conformational Control of Substrate Orientation

In cullin-RING E3 ubiquitin ligases, substrate binding proteins, such as VHL-box, SOCS-box or the F-box proteins, recruit substrates for ubiquitination, accurately positioning and orienting the substrates for ubiquitin transfer. Yet, how the E3 machinery precisely positions the substrate is unknown. Here, we simulated nine substrate binding proteins: Skp2, Fbw7, β-TrCP1, Cdc4, Fbs1, TIR1, pVHL, SOCS2, and SOCS4, in the unbound form and bound to Skp1, ASK1 or Elongin C. All nine proteins have two domains: one binds to the substrate; the other to E3 ligase modules Skp1/ASK1/Elongin C. We discovered that in all cases the flexible inter-domain linker serves as a hinge, rotating the substrate binding domain, optimally and accurately positioning it for ubiquitin transfer. We observed a conserved proline in the linker of all nine proteins. In all cases, the prolines pucker substantially and the pucker is associated with the backbone rotation toward the E2/ubiquitin. We further observed that the linker flexibility could be regulated allosterically by binding events associated with either domain. We conclude that the flexible linker in the substrate binding proteins orients the substrate for the ubiquitin transfer. Our findings provide a mechanism for ubiquitination and polyubiquitination, illustrating that these processes are under conformational control

Linkage between fitness of yeast cells and adenylate kinase catalysis

Enzymes have evolved with highly specific values of their catalytic parameters kcat and KM. This poses fundamental biological questions about the selection pressures responsible for evolutionary tuning of these parameters. Here we are address these questions for the enzyme adenylate kinase (Adk) in eukaryotic yeast cells. A plasmid shuffling system was developed to allow quantification of relative fitness (calculated from growth rates) of yeast in response to perturbations of Adk activity introduced through mutations. Biophysical characterization verified that all variants studied were properly folded and that the mutations did not cause any substantial differences to thermal stability. We found that cytosolic Adk is essential for yeast viability in our strain background and that viability could not be restored with a catalytically dead, although properly folded Adk variant. There exist a massive overcapacity of Adk catalytic activity and only 12% of the wild type kcat is required for optimal growth at the stress condition 20°C. In summary, the approach developed here has provided new insights into the evolutionary tuning of kcat for Adk in a eukaryotic organism. The developed methodology may also become useful for uncovering new aspects of active site dynamics and also in enzyme design since a large library of enzyme variants can be screened rapidly by identifying viable colonies

Changes in Dynamics upon Oligomerization Regulate Substrate Binding and Allostery in Amino Acid Kinase Family Members

Oligomerization is a functional requirement for many proteins. The interfacial interactions and the overall packing geometry of the individual monomers are viewed as important determinants of the thermodynamic stability and allosteric regulation of oligomers. The present study focuses on the role of the interfacial interactions and overall contact topology in the dynamic features acquired in the oligomeric state. To this aim, the collective dynamics of enzymes belonging to the amino acid kinase family both in dimeric and hexameric forms are examined by means of an elastic network model, and the softest collective motions (i.e., lowest frequency or global modes of motions) favored by the overall architecture are analyzed. Notably, the lowest-frequency modes accessible to the individual subunits in the absence of multimerization are conserved to a large extent in the oligomer, suggesting that the oligomer takes advantage of the intrinsic dynamics of the individual monomers. At the same time, oligomerization stiffens the interfacial regions of the monomers and confers new cooperative modes that exploit the rigid-body translational and rotational degrees of freedom of the intact monomers. The present study sheds light on the mechanism of cooperative inhibition of hexameric N-acetyl-L-glutamate kinase by arginine and on the allosteric regulation of UMP kinases. It also highlights the significance of the particular quaternary design in selectively determining the oligomer dynamics congruent with required ligand-binding and allosteric activities

Rapid Sampling of Molecular Motions with Prior Information Constraints

Proteins are active, flexible machines that perform a range of different functions. Innovative experimental approaches may now provide limited partial information about conformational changes along motion pathways of proteins. There is therefore a need for computational approaches that can efficiently incorporate prior information into motion prediction schemes. In this paper, we present PathRover, a general setup designed for the integration of prior information into the motion planning algorithm of rapidly exploring random trees (RRT). Each suggested motion pathway comprises a sequence of low-energy clash-free conformations that satisfy an arbitrary number of prior information constraints. These constraints can be derived from experimental data or from expert intuition about the motion. The incorporation of prior information is very straightforward and significantly narrows down the vast search in the typically high-dimensional conformational space, leading to dramatic reduction in running time. To allow the use of state-of-the-art energy functions and conformational sampling, we have integrated this framework into Rosetta, an accurate protocol for diverse types of structural modeling. The suggested framework can serve as an effective complementary tool for molecular dynamics, Normal Mode Analysis, and other prevalent techniques for predicting motion in proteins. We applied our framework to three different model systems. We show that a limited set of experimentally motivated constraints may effectively bias the simulations toward diverse predicates in an outright fashion, from distance constraints to enforcement of loop closure. In particular, our analysis sheds light on mechanisms of protein domain swapping and on the role of different residues in the motion

Structure and Dynamics of the G121V Dihydrofolate Reductase Mutant: Lessons from a Transition-State Inhibitor Complex

It is well known that enzyme flexibility is critical for function. This is due to the observation that the rates of intramolecular enzyme motions are often matched to the rates of intermolecular events such as substrate binding and product release. Beyond this role in progression through the reaction cycle, it has been suggested that enzyme dynamics may also promote the chemical step itself. Dihydrofolate reductase (DHFR) is a model enzyme for which dynamics have been proposed to aid in both substrate flux and catalysis. The G121V mutant of DHFR is a well studied form that exhibits a severe reduction in the rate of hydride transfer yet there remains dispute as to whether this defect is caused by altered structure, dynamics, or both. Here we address this by presenting an NMR study of the G121V mutant bound to reduced cofactor and the transition state inhibitor, methotrexate. NMR chemical shift markers demonstrate that this form predominantly adopts the closed conformation thereby allowing us to provide the first glimpse into the dynamics of a catalytically relevant complex. Based on 15N and 2H NMR spin relaxation, we find that the mutant complex has modest changes in ps-ns flexibility with most affected residues residing in the distal adenosine binding domain rather than the active site. Thus, aberrant ps-ns dynamics are likely not the main contributor to the decreased catalytic rate. The most dramatic effect of the mutation involves changes in µs-ms dynamics of the F-G and Met20 loops. Whereas loop motion is quenched in the wild type transition state inhibitor complex, the F-G and Met20 loops undergo excursions from the closed conformation in the mutant complex. These excursions serve to decrease the population of conformers having the correct active site configuration, thus providing an explanation for the G121V catalytic defect

Evidence of an Antimicrobial-Immunomodulatory Role of Atlantic Salmon Cathelicidins during Infection with Yersinia ruckeri

Cathelicidins are a family of antimicrobial peptides that act as effector molecules of the innate immune system with broad-spectrum antimicrobial properties. These evolutionary conserved cationic host-defence peptides are integral components of the immune response of fish, which are generally believed to rely heavily on innate immune defences to invading pathogens. In this study we showed that Atlantic salmon cathelicidin 1 and 2 (asCATH1 and asCATH2) stimulated peripheral blood leukocytes increasing the transcription of the chemokine interleukin-8. Further, functional differences were identified between the two cathelicidins. In the presence of serum, asCATH1 displayed greatly diminished host haemolytic activity, while the constitutively expressed asCATH2 had no haemolytic activity with or without serum. These findings support our hypothesis that fish cathelicidins exert their primary antimicrobial action at the site of pathogen invasion such as epithelial surfaces. Further, we hypothesise that like their mammalian counterparts in the presence of serum they act as mediators of the innate and adaptive immune response via the release of cytokines thus indirectly protecting against a variety of pathogens. We highlight the importance of this immunomodulatory role from the involvement of asCATHs during an infection with the fish pathogen Yersinia ruckeri. While we were able to demonstrate in vitro that asCATH1 and 2, possessed direct microbicidal activity against the fish pathogen, Vibrio anguillarum, and a common gram negative bacterium, Escherichia coli, little or no bactericidal activity was found against Y. ruckeri. The contribution of either asCATH in the immune response or as a potential virulence factor during yersiniosis is highlighted from the increased expression of asCATH1 and 2 mRNA during an in vivo challenge with Y. ruckeri . We propose that Atlantic salmon cathelicidins participate in the interplay between the innate and adaptive immune systems via the release of cytokines enabling a more effective response to invading pathogens

Role of Active Site Rigidity in Activity: MD Simulation and Fluorescence Study on a Lipase Mutant

Relationship between stability and activity of enzymes is maintained by underlying conformational flexibility. In thermophilic enzymes, a decrease in flexibility causes low enzyme activity while in less stable proteins such as mesophiles and psychrophiles, an increase in flexibility is associated with enhanced enzyme activity. Recently, we identified a mutant of a lipase whose stability and activity were enhanced simultaneously. In this work, we probed the conformational dynamics of the mutant and the wild type lipase, particularly flexibility of their active site using molecular dynamic simulations and time-resolved fluorescence techniques. In contrast to the earlier observations, our data show that active site of the mutant is more rigid than wild type enzyme. Further investigation suggests that this lipase needs minimal reorganization/flexibility of active site residues during its catalytic cycle. Molecular dynamic simulations suggest that catalytically competent active site geometry of the mutant is relatively more preserved than wild type lipase, which might have led to its higher enzyme activity. Our study implies that widely accepted positive correlation between conformation flexibility and enzyme activity need not be stringent and draws attention to the possibility that high enzyme activity can still be accomplished in a rigid active site and stable protein structures. This finding has a significant implication towards better understanding of involvement of dynamic motions in enzyme catalysis and enzyme engineering through mutations in active site