Protein conformational transitions explored by mixed elastic network models

We develop a mixed elastic network model (MENM) to study large-scale conformational transitions of proteins between two (or more) known structures. Elastic network potentials for the beginning and end states of a transition are combined, in effect, by adding their respective partition functions. The resulting effective MENM energy function smoothly interpolates between the original surfaces, and retains the beginning and end structures as local minima. Saddle points, transition paths, potentials of mean force, and partition functions can be found efficiently by largely analytic methods. To characterize the protein motions during a conformational transition, we follow "transition paths" on the MENM surface that connect the beginning and end structures and are invariant to parameterizations of the model and the mathematical form of the mixing scheme. As illustrations of the general formalism, we study large-scale conformation changes of the motor proteins KIF1A kinesin and myosin II. We generate possible transition paths for these two proteins that reveal details of their conformational motions. The MENM formalism is computationally efficient and generally applicable even for large protein systems that undergo highly collective structural changes.     

Fine structure of conformational ensembles in adenylate kinase

Adenylate kinase (ADK) catalyzes the reversible Mg²⁺‐dependent phosphoryl transfer reaction Mg²⁺+2ADP ?Mg²⁺+ATP + AMP in essential cellular systems. This reaction is a major player in cellular energy homeostasis and the isoform network of ADK plays an important role in AMP metabolic signaling circuits. ADK has 3 domains, the LID, NMP, and CORE domains, that undergo large conformational rearrangements during ADK's catalytic cycle. In spite of extensive experimental and computational studies, details of the conformational pathway from open to closed forms remain uncertain. In this paper we explore this pathway using coarse‐grained molecular dynamics (MD) trajectories of ADK calculated by GROMACS using a SMOG model and classify the conformations within the resultant trajectories by K‐means clustering. ADK conformations segregate naturally into open; intermediate; and closed forms with long‐term residence in the intermediate state. Structural clustering divides the intermediate conformation into 3 sub‐states that are distinguished from one another on the basis of differences in both structure and dynamics. These distinctions are defined on the basis of a number of different metrics including radius of gyration, dihedral angle fluctuation, and fluctuations of interatomic pair distances. Furthermore, differences in the sub‐states appear to correspond to the distinct ways each sub‐state contributes to the molecular mechanism of catalysis: One sub‐state acts as a gate‐way to the open conformation; one sub‐state a gate‐way to the closed conformation. A third intermediate sub‐state appears to represent a metastable off‐pathway structure that is nevertheless frequently visited during the passage from open to closed state.     

Activation and conformational changes of adenylate kinase in urea solution

Ａｄｅｎｙｌａｔｅｋｉｎａｓｅ(ＥＣ2.7.4.3)ｃａｔａｌｙｚｅｓｔｈｅｉｎｔｅｒｃｏｎｖｅｒｓｉｏｎｏｆａｄｅｎｉｎｅｎｕｃｌｅｏｔｉｄｅｓａｃｃｏｒｄｉｎｇｔｏ:ＡＤＰ+ＭｇＡＤＰＡＭＰ+ＭｇＡＴＰ.Ｉｔｉｓｕｂｉｑｕｉｔｏｕｓａｎｄｐａｒｔｉｃｕｌａｒｌｙａｂｕｎｄａｎｔｉｎｔｉｓｓｕｅｓｗｉｔｈｈｉｇｈｅｎｅｒｇｙｔｕｒｎｏｖｅｒ.Ｒａｂｂｉｔｍｕｓｃｌｅａｄｅｎｙｌａｔｅｋｉｎａｓｅｉｓａｍｏｎｏｍｅｒｅｎｚｙｍｅｗｉｔｈｔｗｏｓｔｒｕｃｔｕｒｅｄｏｍａｉｎｓｃｏｍｐｏｓｅｄｏｆ194ａｍｉｎｏａ…     

Predicting order of conformational changes during protein conformational transitions using an interpolated elastic network model

The decryption of sequence of structural events during protein conformational transitions is essential to a detailed understanding of molecular functions ofvarious biological nanomachines. Coarse‐grained models have proven useful by allowing highly efficient simulations of protein conformational dynamics. By combining two coarse‐grained elastic network models constructed based on the beginning and end conformations of a transition, we have developed an interpolated elastic network model to generate a transition pathway between the two protein conformations. For validation, we have predicted the order of local and global conformational changes during key ATP‐driven transitions in three important biological nanomachines (myosin, F₁ ATPase and chaperonin GroEL). We have found that the local conformational change associated with the closing of active site precedes the global conformational change leading to mechanical motions. Our finding is in good agreement with the distribution of intermediate experimental structures, and it supports the importance of local motions at active site to drive or gate various conformational transitions underlying the workings of a diverse range of biological nanomachines. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Conformational transitions of adenylate kinase: switching by cracking

Conformational heterogeneity in proteins is known to often be the key to their function. We present a coarse grained model to explore the interplay between protein structure, folding and function which is applicable to allosteric or non-allosteric proteins. We employ the model to study the detailed mechanism of the reversible conformational transition of Adenylate Kinase (AKE) between the open to the closed conformation, a reaction that is crucial to the protein's catalytic function. We directly observe high strain energy which appears to be correlated with localized unfolding during the functional transition. This work also demonstrates that competing native interactions from the open and closed form can account for the large conformational transitions in AKE. We further characterize the conformational transitions with a new measure Phi(Func), and demonstrate that local unfolding may be due, in part, to competing intra-protein interactions.     

Coupling between ATP hydrolysis and protein conformational change in maltose transporter

As the intracellular part of maltose transporter, MalK dimer utilizes the energy of ATP hydrolysis to drive protein conformational change, which then facilitates substrate transport. Free energy evaluation of the complete conformational change before and after ATP hydrolysis is helpful to elucidate the mechanism of chemical‐to‐mechanical energy conversion in MalK dimer, but is lacking in previous studies. In this work, we used molecular dynamics simulations to investigate the structural transition of MalK dimer among closed, semi‐open and open states. We observed spontaneous structural transition from closed to open state in the ADP‐bound system and partial closure of MalK dimer from the semi‐open state in the ATP‐bound system. Subsequently, we calculated the reaction pathways connecting the closed and open states for the ATP‐ and ADP‐bound systems and evaluated the free energy profiles along the paths. Our results suggested that the closed state is stable in the presence of ATP but is markedly destabilized when ATP is hydrolyzed to ADP, which thus explains the coupling between ATP hydrolysis and protein conformational change of MalK dimer in thermodynamics. Proteins 2017; 85:207–220. © 2016 Wiley Periodicals, Inc.     

Escherichia coli adenylate kinase dynamics: comparison of elastic network model modes with mode-coupling (15)N-NMR relaxation data

The dynamics of adenylate kinase of Escherichia coli (AKeco) and its complex with the inhibitor AP(5)A, are characterized by correlating the theoretical results obtained with the Gaussian Network Model (GNM) and the anisotropic network model (ANM) with the order parameters and correlation times obtained with Slowly Relaxing Local Structure (SRLS) analysis of (15)N-NMR relaxation data. The AMPbd and LID domains of AKeco execute in solution large amplitude motions associated with the catalytic reaction Mg(+2)*ATP + AMP --> Mg(+2)*ADP + ADP. Two sets of correlation times and order parameters were determined by NMR/SRLS for AKeco, attributed to slow (nanoseconds) motions with correlation time tau( perpendicular) and low order parameters, and fast (picoseconds) motions with correlation time tau( parallel) and high order parameters. The structural connotation of these patterns is examined herein by subjecting AKeco and AKeco*AP(5)A to GNM analysis, which yields the dynamic spectrum in terms of slow and fast modes. The low/high NMR order parameters correlate with the slow/fast modes of the backbone elucidated with GNM. Likewise, tau( parallel) and tau( perpendicular) are associated with fast and slow GNM modes, respectively. Catalysis-related domain motion of AMPbd and LID in AKeco, occurring per NMR with correlation time tau( perpendicular), is associated with the first and second collective slow (global) GNM modes. The ANM-predicted deformations of the unliganded enzyme conform to the functional reconfiguration induced by ligand-binding, indicating the structural disposition (or potential) of the enzyme to bind its substrates. It is shown that NMR/SRLS and GNM/ANM analyses can be advantageously synthesized to provide insights into the molecular mechanisms that control biological function.     

Iterative cluster-NMA: A tool for generating conformational transitions in proteins

Computational models provide insight into the structure-function relationship in proteins. These approaches, especially those based on normal mode analysis, can identify the accessible motion space around a given equilibrium structure. The large magnitude, collective motions identified by these methods are often well aligned with the general direction of the expected conformational transitions. However, these motions cannot realistically be extrapolated beyond the local neighborhood of the starting conformation. In this article, the iterative cluster-NMA (icNMA) method is presented for traversing the energy landscape from a starting conformation to a desired goal conformation. This is accomplished by allowing the evolving geometry of the intermediate structures to define the local accessible motion space, and thus produce an appropriate displacement. Following the derivation of the icNMA method, a set of sample simulations are performed to probe the robustness of the model. A detailed analysis of beta1,4-galactosyltransferase-T1 is also given, to highlight many of the capabilities of icNMA. Remarkably, during the transition, a helix is seen to be extended by an additional turn, emphasizing a new unknown role for secondary structures to absorb slack during transitions. The transition pathway for adenylate kinase, which has been frequently studied in the literature, is also discussed.     

Associative mechanism for phosphoryl transfer: a molecular dynamics simulation of Escherichia coli adenylate kinase complexed with its substrates

The ternary complex of Escherichia coli adenylate kinase (ECAK) with its substrates adenosine monophosphate (AMP) and Mg-ATP, which catalyzes the reversible transfer of a phosphoryl group between adenosine triphosphate (ATP) and AMP, was studied using molecular dynamics. The starting structure for the simulation was assembled from the crystal structures of ECAK complexed with the bisubstrate analog diadenosine pentaphosphate (AP(5)A) and of Bacillus stearothermophilus adenylate kinase complexed with AP(5)A, Mg(2+), and 4 coordinated water molecules, and by deleting 1 phosphate group from AP(5)A. The interactions of ECAK residues with the various moieties of ATP and AMP were compared to those inferred from NMR, X-ray crystallography, site-directed mutagenesis, and enzyme kinetic studies. The simulation supports the hypothesis that hydrogen bonds between AMP's adenine and the protein are at the origin of the high nucleoside monophosphate (NMP) specificity of AK. The ATP adenine and ribose moieties are only loosely bound to the protein, while the ATP phosphates are strongly bound to surrounding residues. The coordination sphere of Mg(2+), consisting of 4 waters and oxygens of the ATP beta- and gamma-phosphates, stays approximately octahedral during the simulation. The important role of the conserved Lys13 in the P loop in stabilizing the active site by bridging the ATP and AMP phosphates is evident. The influence of Mg(2+), of its coordination waters, and of surrounding charged residues in maintaining the geometry and distances of the AMP alpha-phosphate and ATP beta- and gamma-phosphates is sufficient to support an associative reaction mechanism for phosphoryl transfer.     

Substrate-binding loop interactions with pseudouridine trigger conformational changes that promote catalytic efficiency of pseudouridine kinase PUKI

Pseudouridine, one major RNA modification, is catabolized into uracil and ribose-5′-phosphate by two sequential enzymatic reactions. In the first step, pseudouridine kinase (PUKI) phosphorylates pseudouridine to pseudouridine 5′-monophosphate. High-fidelity catalysis of pseudouridine by PUKI prevents possible disturbance of in vivo pyrimidine homeostasis. However, the molecular basis of how PUKI selectively phosphorylates pseudouridine over uridine with >100-fold greater efficiency despite minor differences in their K_m values has not been elucidated. To investigate this selectivity, in this study we determined the structures of PUKI from Escherichia coli strain B (EcPUKI) in various ligation states. The structure of EcPUKI was determined to be similar to PUKI from Arabidopsis thaliana, including an α/β core domain and β-stranded small domain, with dimerization occurring via the β-stranded small domain. In a binary complex, we show that Ser30 in the substrate-binding loop of the small domain mediates interactions with the hallmark N1 atom of pseudouridine nucleobase, causing conformational changes in its quaternary structure. Kinetic and fluorescence spectroscopic analyses also showed that the Ser30-mediated interaction is a prerequisite for conformational changes and subsequent catalysis by EcPUKI. Furthermore, S30A mutation or EcPUKI complexed with other nucleosides homologous to pseudouridine but lacking the pseudouridine-specific N1 atom did not induce such conformational changes, demonstrating the catalytic significance of the proposed Ser30-mediated interaction. These analyses provide structural and functional evidence for a pseudouridine-dependent conformational change of EcPUKI and its functional linkage to catalysis.     

Analysis of conformational motions and related key residue interactions responsible for a specific function of proteins with elastic network model

Protein collective motions play a critical role in many biochemical processes. How to predict the functional motions and the related key residue interactions in proteins is important for our understanding in the mechanism of the biochemical processes. Normal mode analysis (NMA) of the elastic network model (ENM) is one of the effective approaches to investigate the structure-encoded motions in proteins. However, the motion modes revealed by the conventional NMA approach do not necessarily correspond to a specific function of protein. In the present work, a new analysis method was proposed to identify the motion modes responsible for a specific function of proteins and then predict the key residue interactions involved in the functional motions by using a perturbation approach. In our method, an internal coordinate that accounts for the specific function was introduced, and the Cartesian coordinate space was transformed into the internal/Cartesian space by using linear approximation, where the introduced internal coordinate serves as one of the axes of the coordinate space. NMA of ENM in this internal/Cartesian space was performed and the function-relevant motion modes were identified according to their contributions to the specific function of proteins. Then the key residue interactions important for the functional motions of the protein were predicted as the interactions whose perturbation largely influences the fluctuation along the internal coordinate. Using our proposed methods, the maltose transporter (MalFGK₂) from E. Coli was studied. The functional motions and the key residue interactions that are related to the channel-gating function of this protein were successfully identified.     

Effects of zinc on creatine kinase: activity changes, conformational changes, and aggregation

The effects of zinc on creatine kinase (CK) are very distinctive compared with other bivalent metal ions. Zinc up to 0.1 mM induced increases in CK activity, accompanied by significant hydrophobic surface exposure and increase in a-helix content of CK. Zinc over 0.1 mM denatured and inactived CK. In the presence of 0.1 mM zinc, the CK activity was very close to that of the native CK, but its conformation changed greatly. The kinetic courses of CK inactivation and conformational change in the presence of 1 mM zinc were measured to determine apparent rate constants of inactivation and conformational change. Zinc over 0.05 mM induced CK aggregation at 37°C, and the aggregation was dependent on zinc concentration, CK concentration, and temperature. The inactivation and aggregation can be reversed by EDTA. An explanation for CK aggregation induced by zinc is proposed, as well as a mechanism for CK abnormality in Alzheimer's disease.To whom correspondence should be addressed.     

hCINAP is an atypical mammalian nuclear adenylate kinase with an ATPase motif: structural and functional studies

Human coilin interacting nuclear ATPase protein (hCINAP) directly interacts with coilin, a marker protein of Cajal Bodies (CBs), nuclear organelles involved in the maturation of small nuclear ribonucleoproteins UsnRNPs and snoRNPs. hCINAP has previously been designated as an adenylate kinase (AK6), but is very atypical as it exhibits unusually broad substrate specificity, structural features characteristic of ATPase/GTPase proteins (Walker motifs A and B) and also intrinsic ATPase activity. Despite its intriguing structure, unique properties and cellular localization, the enzymatic mechanism and biological function of hCINAP have remained poorly characterized. Here, we offer the first high-resolution structure of hCINAP in complex with the substrate ADP (and dADP), the structure of hCINAP with a sulfate ion bound at the AMP binding site, and the structure of the ternary complex hCINAP-Mg(2+) ADP-Pi. Induced fit docking calculations are used to predict the structure of the hCINAP-Mg(2+) ATP-AMP ternary complex. Structural analysis suggested a functional role for His79 in the Walker B motif. Kinetic analysis of mutant hCINAP-H79G indicates that His79 affects both AK and ATPase catalytic efficiency and induces homodimer formation. Finally, we show that in vivo expression of hCINAP-H79G in human cells is toxic and drastically deregulates the number and appearance of CBs in the cell nucleus. Our findings suggest that hCINAP may not simply regulate nucleotide homeostasis, but may have broader functionality, including control of CB assembly and disassembly in the nucleus of human cells.     

Study of global motions in proteins by weighted masses molecular dynamics: Adenylate kinase as a test case

The weighted masses molecular dynamics (WMMD) technique is applied to the protein adenylate kinase. A novel set of restraints has been developed to allow the use of this technique with proteins. The WMMD simulation is successful in predicting the flexibility of the two mobile domains of the protein. The end product of the simulation is similar to the known open and AMP bound forms of the enzyme. The biological relevance of the restraints used and potential methods of improving the technique are discussed. © 1996 Wiley-Liss, Inc.