BioMed Research International

Volume 2013 (2013), Article ID 628536, 7 pages

http://dx.doi.org/10.1155/2013/628536

Research Article

Molecular Dynamics Studies on the Conformational Transitions of Adenylate Kinase: A Computational Evidence for the Conformational Selection Mechanism

Jie Ping,¹ Pei Hao,1,2 Yi-Xue Li,2,3 and Jing-Fang Wang2,4

¹Pathogen Diagnostic Center, Institut Pasteur of Shanghai Chinese Academy of Sciences, Shanghai 200025, China

²Shanghai Center for Bioinformation Technology, 100 Qinzhou Road, Shanghai 200235, China

³Bioinformatics Center, Key Laboratory of Systems Biology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai 200031, China

⁴Key Laboratory of Systems Biomedicine (Ministry of Education), Shanghai Center for Systems Biomedicine, Shanghai Jiao Tong University, Shanghai 200240, China

Received 5 April 2013; Accepted 13 June 2013

Academic Editor: Yudong Cai

Copyright © 2013 Jie Ping et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Escherichia coli adenylate kinase (ADK) is a monomeric phosphotransferase enzyme that catalyzes reversible transfer of phosphoryl group from ATP to AMP with a large-scale domain motion. The detailed mechanism for this conformational transition remains unknown. In the current study, we performed long time-scale molecular dynamics simulations on both open and closed states of ADK. Based on the structural analyses of the simulation trajectories, we detected over 20 times conformational transitions between the open and closed states of ADK and identified two novel conformations as intermediate states in the catalytic processes. With these findings, we proposed a possible mechanism for the large-scale domain motion of Escherichia coli ADK and its catalytic process: (1) the substrate free ADK adopted an open conformation; (2) ATP bound with LID domain closure; (3) AMP bound with NMP domain closure; (4) phosphoryl transfer occurred with ATP, and AMP converted into two ADPs, and no conformational transition was detected in the enzyme; (5) LID domain opened with one ADP released; (6) another ADP released with NMP domain open. As both open and closed states sampled a wide range of conformation transitions, our simulation strongly supported the conformational selection mechanism for Escherichia coli ADK.

1. Introduction

Escherichia coli adenylate kinase (ADK) is a monomeric phosphotransferase enzyme regulating energy homeostasis in cells by catalyzing reversible transfer of a phosphoryl group from ATP to AMP. The structure of ADK is well studied in the past few years, and by now more than 20 crystal structures of ADK from Escherichia coli and other organisms in the absence and presence of substrates have been released in the protein structure databases. According to these crystal studies, the structure of ADK is mainly composed of three major components (Figure 1): a core domain (residues 1–29, 68–117, and 161–214), an AMP binding lid domain (also called NMP domain, residues 30–67), as well as an ATP binding lid domain (also called LID domain, residues 118–167). These basic components are shared by many other kinases and ATPases. As characterized by structural [1–4], biophysical [5, 6], and computational studies [7–10], ADK is believed to adopt an open conformation in the absence of substrates (Figure 1(a)), and with ATP or AMP binding the LID and NMP domains of this enzyme undergo large conformational transitions, leading to a closed conformation (Figure 1(b)). However, during the large conformational transitions, the core domain is comparatively rigid with the structural motions mainly located on the LID and NMP domains. This large conformational transition, especially opening of the nucleotide binding lids required for the catalysis and product release, occur on the microsecond-millisecond times scale [11].

(a)

(b)

Figure 1: The three-dimensional structures of Escherichia coli adenylate kinase in the open (a) and closed (b) conformations. The crystal structures 4ake.pdb and 1ake.pdb are selected as the open and closed state of ADK, respectively. NMP (residues 30–67), CORE (residues 1–29, 68–117, and 161–214), and LID (residues 118–167) domains are colored in red, yellow, and blue, respectively.

Large-scale conformational alterations are thought to mediate allosteric regulation, which are related to the protein function in signal transduction, immune response, and enzymatic activity. Thus, a fundamental problem is to understand the mechanism for the large-scale conformational transitions. To answer this question, two models are proposed: the induced-fit model and the population-shift model. In the induced-fit model, substrate binding is believed to induce conformational alterations in the active site to cause a new conformation for the entire enzyme. Contrarily, in the population-shift model, the enzyme is thought to adopt a conformational equilibrium among many native states, and substrates are able to selectively bind to a suitable native conformation, shifting the equilibrium toward the binding conformation.

Several computational works have been done in recent years to study the conformational transitions of ADK. All atomic molecular dynamics simulations have been successfully used to reveal the structural trajectories of ADK during the large conformational transitions [7–9, 12–21]. However, these studies still remain on the nanosecond times scale much smaller than the times scale on which the large conformational transitions of ADK occur. Additionally, coarse-grained simulations [22–26], normal mode analysis [27–29], and plastic/elastic network approaches [30–33] are also employed to generate transition pathways of ADK. However, these theoretical studies are controversial for adopting the simple harmonic potential approximation instead of all-atom models, which is less convincing due to the intermediate structures far away from the native states. To this end, we performed long time-scale molecular dynamics simulations on both open and closed state of ADK to study the large conformational transitions.

2. Materials and Methods

2.1. Initial Structures

Owing to playing an important role in cellular energy homeostasis, many good attempts have been made to study the three-dimensional (3D) structure of ADK. By now, 29 crystal structures of ADK and its homolog, in the absence and presence of substrates, are available in the protein structure databases [34–45]. In this study, the crystal structures (PDB ID: 4AKE [2] and 1AKE [46]) were selected as the initial structures for the further molecular dynamics simulations. The former was released in 1996 with 2.2 Å resolutions and treated as an open state of ADK in the current case. The latter was obtained in the presence of Ap₅A with 1.9 Å resolutions and considered as a closed state of ADK in this study. Except for the polar hydrogen and heavy atoms of the enzyme and Ap₅A, all other atoms including nonpolar hydrogen in both crystal structures were removed. The pKa values for each residue in ADK were calculated by Delphi [47, 48] as a Poisson-Boltzmann solver with a dielectric constant of 4. Hydrogen atoms were subsequently added to the enzyme with t-Leap procedure of Amber 11 package [49] based on the computational pKa values in the last step to give a total charge of −4. Then, the enzyme (together with Ap₅A in the closed structure) was solvated in a simulation box with explicit TIP3P water molecules. To neutralize the system, 4 sodium ions were added to random place 4 water molecules in the simulation box. The atoms of the enzyme were parameterized by Amber force field parameters [50], while Ap₅A was done by Antechamber module [51] in Amber 11 package.

2.2. Molecular Dynamics Simulations

Two simulation systems were involved in this study, with the crystal structures (PDB ID: 4AKE and 1AKE) as initial structures for the open and closed conformations, respectively. Both systems were subjected to the steepest descent energy minimization (~5000 steps) followed by conjugate gradient energy minimization for the next 5000 steps and subsequently equilibrated with the enzyme atoms (or ligand in the closed system) fixed for a short-time molecular dynamics simulations at 300 K to reduce the van der Waals conflicts. Finally, 10-nanosecond (ns) molecular dynamics simulations were performed for both systems under a constant temperature (300 K) by Amber 11 package [49] with periodic boundary conditions and NPT ensemble. Tweenty five frames were randomly selected from the first 10 ns simulation trajectories for both systems, and 10 independent molecular dynamics simulations (10 ns) were launched for each frame. For the simulations that were detected to have large scale conformational transitions, 200 ns molecular dynamics simulations were added using a GPU-accelerated approach. SHAKE algorithm with a tolerance of 10⁻⁶ was applied to constrain all bonds in both simulation systems [52, 53], and atom velocities for start-up runs were obtained according to the Maxwell distribution at 300 K [54, 55]. The isothermal compressibility was set to 4.5 × 10⁻⁵/bar for solvent simulations [56, 57]. The electrostatic interactions were treated by the particle mesh Ewald (PME) algorithm with interpolation order of 4 and a grid spacing of 0.12 nanometers (nm) [58, 59]. The van der Waals interactions were calculated by using a cut-off of 12 Å. All the molecular dynamics simulations were performed with a time step of 2 femtoseconds (fs), and coordinates for both systems were saved every 1 picosecond (ps).

3. Results and Discussion

3.1. Conformational Transitions

As ADK has been well studied, a lot of kinetic and thermodynamic experimental data were released [60–64], based on which the large-scale conformational transitions for the NMP and LID domain opening are considered to be rate limiting for the catalytic reactions of the enzyme. However, it is really a challenging task to capture large-scale and long-time conformational transitions of proteins or enzymes for both experimental and computational methods. To check whether our simulations involved the large-scale conformational transitions of ADK, we calculated the RMS deviations of the Cα atoms from both open (4ake.pdb) and closed (1ake.pdb) conformations of ADK along all the molecular dynamics simulation trajectories. As shown in Figure 2, the RMS deviation employed a ribbon distribution, indicating that the large-scale conformational transitions occurred in our molecular dynamics simulations. If no conformational transitions occurred, the 2D-RMS deviations should be divided into two major parts. Additionally, the 2D-RMS deviations located in the middle of the ribbon distribution might stand for the intermediate states. As reported by Shapiro and Meirovitch [64], the rate for ADK domain motions was about 52 ns. While we employed a series of simulations with a time scale of 200 ns, almost 4 times larger than the measured rates for ADK domain motions, thus it was expected that the conformational transitions should be detected in our molecular dynamics simulations.

Figure 2: The RMS deviations of Cα atoms from both open and closed conformation along the molecular dynamics trajectories. The open (4ake.pdb) and closed (1ake.pdb) conformations of ADK were specially labeled. The ribbon RMS deviation distribution indicated that large-scale conformational transitions were detected during our molecular dynamics simulations.

To confirm the conformational transitions for Escherichia coli ADK, we also measured the domain motion of ADK using the geometric center distances of NMP and LID domains with CORE domain along the molecular dynamics simulation trajectories. In the crystal structures, the geometric center distances of NMP and CORE domains in the open (4ake.pdb) and closed (1ake.pdb) conformations were 62.7 Å and 18.4 Å, while those for LID and CORE domain were 70.1 Å (open state) and 21.0 Å (closed state), respectively. As shown in Figure 3, the domain distances were clustered into four major parts. The ones colored in yellow and red, respectively, represented the open and closed conformations defined by the crystal structures, whereas the ones colored in blue and green were identified as novel configurations from the open and closed states. The blue ones in Figure 3 employed an open LID domain with a closed NMP domain, while the green ones had an open NMP domain with a closed LID domain. We believed that the blue and green parts in Figure 3 might be the intermediate states of ADK domain motions.

Figure 3: The mass center distances of NMP and LID domain with LID domain. The geometric center distances of NMP and CORE domains in the open and closed conformations were 62.7 and 18.4 Å, while those of LID and CORE domain for the open and closed conformations were 70.1 and 21.0 Å, respectively.

3.2. Conformational Transition Pathway

The large-scale domain motions were detected in the catalytic conversion of Mg²⁺–ATP + AMP → Mg²⁺–ADP + ADP by Escherichia coli ADK [65, 66]. Thus, the possible conformational transition pathway was thought to be associated with the catalytic mechanism for Escherichia coli ADK. According to the previous theoretical studies, LID domain motion was believed to precede NMP domain motion [6, 26, 33, 66]. One possible reason was that there was a stable salt bridge D118-K136 between LID and CORE, which could connect both domains by the strong contributions to the total enthalpic interactions. This salt bridge, also detected in several adenylate kinase structures of different species, was thought to have stabilizing function for the open conformation [9]. Another possible reason was that ATP binding in LID domain was found to have ability to assist NMP domain motion [66]. Based on these points, we proposed a possible conformational transition pathway and catalytic mechanism for Escherichia coli ADK (Figure 4). The substrate free enzyme adopted an open state. After ATP binding, ADK would close its LID domain first to form a novel configuration (model b in Figure 4 and intermediate state 2 in Figure 3). The closure motion of LID domain allowed AMP to bind in NMP domain binding site, resulting in the NMP domain closure (model c in Figure 4). With ATP and AMP binding, the enzyme would adopt a complete closed conformation. After phosphoryl transfer, ATP and AMP were converted into two ADP molecules. In this step, no conformational transition occurred, and thus the enzyme remained in its closed state (model d in Figure 4). When an ADP molecule was released from the ATP binding site, LID domain altered its closed configuration into an open state (model e in Figure 4 and intermediate state 1 in Figure 3). Subsequently, NMP domain opened via the bending of α2 helix toward α4 helix of CORE domain by nearly 15 degrees with respect to α3 helix, and the enzyme would adopt a complete open conformation (model a in Figure 4). This opening process of NMP domain (also AMP binding cleft) was considered to be involved in facilitating an efficient release of the formed product after catalysis.

Figure 4: Conformational transition pathway and proposed catalytic mechanism of ADK. Model a, substrate free ADK with an open conformation. Model b, ATP bound form of ADK with a closed LID domain. Model c, ATP and AMP bound form of ADK with a closed conformation. Model d, two ADP bound forms of ADK with a closed conformation. Model e, one ADP bound form of ADK with a closed NMP domain.

Our proposed mechanism was also supported by some experimental results. Firstly, the complete closed conformation of ADK was found to be the only product-forming state [24]. In other words, to achieve its biological functions, ADK had to convert the open state (substrate free state) into the closed state. Secondly, The ATP binding site in LID domain was able to accommodate ATP, ADP, as well as AMP. ATP employed the highest binding affinity, whereas AMP had the lowest binding affinity [26]. Thirdly, the AMP binding site in NMP domain could only accommodate ADP and AMP [26]. Finally, the AMP binding site could only accommodate AMP in the intermediate state 2 with a closed LID domain and an open NMP domain [18]. Along the aforementioned pathway, a secondary structure analysis of the enzyme structure was performed using DSSP package to check the structural stability during the large-scale conformational transitions. As expected, no significant alterations in secondary structures were detected in CORE and NMP domain during the large-scale conformational transitions. For LID domain, some residues were found to bend and turn with few alterations in secondary structures. Thus, it was believed that ADK maintained its integrate structure with minute changes in its secondary structures, indicating that this enzyme behaved as a rigid body with flexible domains, and pathway detected in our study was reasonable.

Furthermore, our proposed mechanism was also in good agreement with previous theoretical studies. In 2008, Lu and Wang developed a coarse grained two-well model to study the conformational dynamics of ADK in microscopic detail [19]. They identified the LID-closing and NMP-closing pathways, providing theoretical evidence to our model (especially the transitions between model a and b and between model b and c). In 2008, Kubitzki and de Groot used TEE-REX molecular dynamics to study the conformational transitions of ADK and detected a pathway for the open-to-closed transitions [9]. Although a complete transition was not observed, the pathway they found was in accord with the transitions from model a to model c. In 2010, Prof. Head-Gordon and her co-workers used normal modes analyses to study the transition pathway [27]. Although a complete transition pathway was identified, they did not provide any explanation for the catalytic mechanism of ADK. Besides the aforementioned studies, many other works were involved in the free energy calculations during the conformational transitions [19–22], giving an indication that there were intermediate states during the conformational transitions.

3.3. Conformational Transition Mechanism

For the conformational transitions in proteins or enzymes, two different substrate binding mechanisms, induced-fit model, and conformational selection model can be proposed. The former is a model for enzyme-substrate interactions to describe that only the proper substrate is capable of inducing the conformational changes in the binding site, allowing the enzyme to perform its catalytic function [67]. The latter believes that the proteins or enzymes exist in the multiple conformations in the vicinity of its native state, and the ligand selectively binds to an active conformation to shift the equilibrium toward the binding conformation [68, 69]. In the current case, long time-scale molecular dynamics simulations were performed on Escherichia coli ADK, which can sample a large set of conformations between the open and closed states. Our computational results strongly supported the conformational selection model proposed for ADK [5, 70–74].

4. Conclusion

In conclusion, to study the large-scale domain motions in Escherichia coli ADK, we performed long time-scale molecular dynamics on both open and closed states of ADK. The two-dimensional RMS deviations of the Cα atoms from both open and closed conformations along simulation trajectories confirmed the fact that conformational transitions between the open and closed conformation occurred during our simulations. Additionally, two significant intermediate states were identified by monitoring the domain distances between LID/NMP and CORE domains, one of which adopted an open LID domain with a closed NMP domain, and the other one employed an open NMP domain with a closed LID domain. Based on these computational results, we proposed a possible mechanism for the large scale conformational transitions and the catalytic function (Figure 4). The proposed mechanism was in good accordance with the previous experimental and theoretical studies, providing strong support to the conformational selection mechanism for Escherichia coli ADK.

Acknowledgments

This work was supported by Grants from the National Basic Research Program of China (973 Program, no. 2012CB517900), National Key Basic Research Program (no. 30800210), National Natural Science Foundation of China (nos. 31200547 and 90913009), Shanghai Pujiang Scholarship Program (no. 10PJ1408000), Doctoral Program Foundation of Institutions of Higher Education of China (no. 20110073120078), and China Postdoctoral Science Foundation (no. 2012M520949). The authors gratefully acknowledge the support of SA-SIBS scholarship Program.

References

C. Vonrhein, G. J. Schlauderer, and G. E. Schulz, “Movie of the structural changes during a catalytic cycle of nucleoside monophosphate kinases,” Structure, vol. 3, no. 5, pp. 483–490, 1995.

C. W. Müller, G. J. Schlauderer, J. Reinstein, and G. E. Schulz, “Adenylate kinase motions during catalysis: an energetic counterweight balancing substrate binding,” Structure, vol. 4, no. 2, pp. 147–156, 1996.

G. J. Schlauderer, K. Proba, and G. E. Schulz, “Structure of a mutant adenylate kinase ligated with an ATP-analogue showing domain closure over ATP,” Journal of Molecular Biology, vol. 256, no. 2, pp. 223–227, 1996.

G. J. Schlauderer and G. E. Schulz, “The structure of bovine mitochondrial adenylate kinase: comparison with isoenzymes in other compartments,” Protein Science, vol. 5, no. 3, pp. 434–441, 1996.

K. A. Henzler-Wildman, V. Thai, M. Lei et al., “Intrinsic motions along an enzymatic reaction trajectory,” Nature, vol. 450, no. 7171, pp. 838–844, 2007.

J. A. Hanson, K. Duderstadt, L. P. Watkins et al., “Illuminating the mechanistic roles of enzyme conformational dynamics,” Proceedings of the National Academy of Sciences of the United States of America, vol. 104, no. 46, pp. 18055–18060, 2007.

K. Arora and C. L. Brooks, “Large-scale allosteric conformational transitions of adenylate kinase appear to involve a population-shift mechanism,” Proceedings of the National Academy of Sciences of the United States of America, vol. 104, no. 47, pp. 18496–18501, 2007.

C. Snow, G. Y. Qi, and S. Hayward, “Essential dynamics sampling study of adenylate kinase: comparison to citrate synthase and implication for the hinge and shear mechanisms of domain motions,” Proteins, vol. 67, no. 2, pp. 325–337, 2007.

M. B. Kubitzki and B. L. de Groot, “The atomistic mechanism of conformational transition in adenylate kinase: a TEE-REX molecular dynamics study,” Structure, vol. 16, no. 8, pp. 1175–1182, 2008.

K. A. Henzler-Wildman, M. Lei, V. Thai, S. J. Kerns, M. Karplus, and D. Kern, “A hierarchy of timescales in protein dynamics is linked to enzyme catalysis,” Nature, vol. 450, no. 7171, pp. 913–916, 2007.

M. Wolf-Watz, V. Thai, K. Henzler-Wildman, G. Hadjipavlou, E. Z. Eisenmesser, and D. Kern, “Linkage between dynamics and catalysis in a thermophilic-mesophilic enzyme pair,” Nature Structural and Molecular Biology, vol. 11, no. 10, pp. 945–949, 2004.

B. Jana, B. V. Adkar, R. Biswas, and B. Bagchi, “Dynamic coupling between the LID and NMP domain motions in the catalytic conversion of ATP and AMP to ADP by adenylate kinase,” Journal of Chemical Physics, vol. 134, no. 3, Article ID 035101, 10 pages, 2011.

J. B. Brokaw and J. W. Chu, “On the roles of substrate binding and hinge unfolding in conformational changes of adenylate kinase,” Biophysical Journal, vol. 99, no. 10, pp. 3420–3429, 2010.

F. Pontiggia, A. Zen, and C. Micheletti, “Small- and large-scale conformational changes of adenylate kinase: a molecular dynamics study of the subdomain motion and mechanics,” Biophysical Journal, vol. 95, no. 12, pp. 5901–5912, 2008.

H. Dong, S. Qin, and H. X. Zhou, “Effects of macromolecular crowding on protein conformational changes,” PLoS Computational Biology, vol. 6, Article ID e1000833, 2010.

H. Lou and R. I. Cukier, “Molecular dynamics of apo-adenylate kinase: a principal component analysis,” Journal of Physical Chemistry B, vol. 110, no. 25, pp. 12796–12808, 2006.

E. Bae and G. N. Phillips Jr., “Identifying and engineering ion pairs in adenylate kinases: Insights from molecular dynamics simulations of thermophilic and mesophilic homologues,” The Journal of Biological Chemistry, vol. 280, no. 35, pp. 30943–30948, 2005.

H. Krishnamurthy, H. Lou, A. Kimple, C. Vieille, and R. I. Cukier, “Associative mechanism for phosphoryl transfer: a molecular dynamics simulation of Escherichia coli adenylate kinase complexed with its substrates,” Proteins, vol. 58, no. 1, pp. 88–100, 2005.

Q. Lu and J. Wang, “Single molecule conformational dynamics of adenylate kinase: energy landscape, structural correlations, and transition state ensembles,” Journal of the American Chemical Society, vol. 130, no. 14, pp. 4772–4783, 2008.

O. Beckstein, E. J. Denning, J. R. Perilla, and T. B. Woolf, “Zipping and unzipping of adenylate kinase: atomistic insights into the ensemble of open ↔ closed transitions,” Journal of Molecular Biology, vol. 394, no. 1, pp. 160–176, 2009.

M. D. Daily, G. N. Phillips Jr., and Q. Cui, “Many local motions cooperate to produce the adenylate kinase conformational transition,” Journal of Molecular Biology, vol. 400, no. 3, pp. 618–631, 2010.

R. Potestio, F. Pontiggia, and C. Micheletti, “Coarse-grained description of protein internal dynamics: an optimal strategy for decomposing proteins in rigid subunits,” Biophysical Journal, vol. 96, no. 12, pp. 4993–5002, 2009.

R. D. Hills Jr., L. Lu, and G. A. Voth, “Multiscale coarse-graining of the protein energy landscape,” PLoS Computational Biology, vol. 6, no. 6, Article ID e1000827, 2010.

O. Miyashita, P. G. Wolynes, and J. N. Onuchic, “Simple energy landscape model for the kinetics of functional transitions in proteins,” Journal of Physical Chemistry B, vol. 109, no. 5, pp. 1959–1969, 2005.

Q. Lu and J. Wang, “Kinetics and statistical distributions of single-molecule conformational dynamics,” Journal of Physical Chemistry B, vol. 113, no. 5, pp. 1517–1521, 2009.

P. C. Whitford, O. Miyashita, Y. Levy, and J. N. Onuchic, “Conformational transitions of adenylate kinase: switching by cracking,” Journal of Molecular Biology, vol. 366, no. 5, pp. 1661–1671, 2007.

C. Peng, L. Zhang, and T. Head-Gordon, “Instantaneous normal modes as an unforced reaction coordinate for protein conformational transitions,” Biophysical Journal, vol. 98, no. 10, pp. 2356–2364, 2010.

A. Ahmed, F. Rippmann, G. Barnickel, and H. Gohlke, “A normal mode-based geometric simulation approach for exploring biologically relevant conformational transitions in proteins,” Journal of Chemical Information and Modeling, vol. 51, no. 7, pp. 1604–1622, 2011.

S. Kirillova, J. Cortés, A. Stefaniu, and T. Siméon, “An NMA-guided path planning approach for computing large-amplitude conformational changes in proteins,” Proteins, vol. 70, no. 1, pp. 131–143, 2008.

N. A. Temiz, E. Meirovitch, and I. Bahar, “Escherichia coli adenylate kinase dynamics: comparison of elastic network model modes with mode-coupling ¹⁵N-NMR relaxation data,” Proteins, vol. 57, no. 3, pp. 468–480, 2004.

J. N. Stember and W. Wriggers, “Bend-twist-stretch model for coarse elastic network simulation of biomolecular motion,” Journal of Chemical Physics, vol. 131, no. 7, Article ID 074112, 9 pages, 2009.

W. J. Zheng, B. R. Brooks, and G. Hummer, “Protein conformational transitions explored by mixed elastic network models,” Proteins, vol. 69, no. 1, pp. 43–57, 2007.

P. Maragakis and M. Karplus, “Large amplitude conformational change in proteins explored with a plastic network model: adenylate kinase,” Journal of Molecular Biology, vol. 352, no. 4, pp. 807–822, 2005.

M. B. Berry and G. N. J. Philips, “Crystal structures of Bacillus stearothermophilus adenylate kinase with bound Ap₅A, Mg²⁺ Ap₅A, and Mn²⁺ Ap₅A reveal an intermediate lid position and six coordinate octahedral geometry for bound Mg²⁺ and Mn²⁺,” Proteins, vol. 32, no. 3, pp. 276–288, 1998.

C. W. Muller and G. E. Schulz, “Crystal structures of two mutants of adenylate kinase from Escherichia coli that modify the Gly-loop,” Proteins, vol. 15, no. 1, pp. 42–49, 1993.

M. B. Berry, E. Bae, T. R. Bilderback, M. Glaser, and G. N. Phillips Jr., “Crystal structure of ADP/AMP complex of Escherichia coli adenylate kinase,” Proteins, vol. 62, no. 2, pp. 555–556, 2006.

M. B. Berry, B. Meador, T. Bilderback, P. Liang, M. Glaser, and G. N. Phillips Jr., “The closed conformation of a highly flexible protein: the structure of E. coli adenylate kinase with bound AMP and AMPPNP,” Proteins, vol. 19, no. 3, pp. 183–198, 1994.

E. Bae and G. N. Phillips Jr., “Structures and analysis of highly homologous psychrophilic, mesophilic, and thermophilic adenylate kinases,” The Journal of Biological Chemistry, vol. 279, no. 27, pp. 28202–28208, 2004.

R. Couñago, S. Chen, and Y. Shamoo, “In vivo molecular evolution reveals biophysical origins of organismal fitness,” Molecular Cell, vol. 22, no. 4, pp. 441–449, 2006.

U. Abele and G. E. Schulz, “High-resolution structures of adenylate kinase from yeast ligated with inhibitor Ap₅A, showing the pathway of phosphoryl transfer,” Protein Science, vol. 4, no. 7, pp. 1262–1271, 1995.

P. Spuergin, U. Abele, and G. E. Schulz, “Stability, activity and structure of adenylate kinase mutants,” European Journal of Biochemistry, vol. 231, no. 2, pp. 405–413, 1995.

H. Gu, H. F. Chen, D. Q. Wei, and J. F. Wang, “Molecular dynamics simulations exploring drug resistance in HIV-1 proteases,” Chinese Science Bulletin, vol. 55, no. 24, pp. 2677–2683, 2010.

K. Wild, R. Grafmüller, E. Wagner, and G. E. Schulz, “Structure, catalytsis and supramolecular assembly of adenylate kinase from maize,” European Journal of Biochemistry, vol. 250, no. 2, pp. 326–331, 1997.

G. J. Schlauderer and G. E. Schulz, “The structure of bovine mitochondrial adenylate kinase: comparison with isoenzymes in other compartments,” Protein Science, vol. 5, no. 3, pp. 434–441, 1996.

K. Diederichs and G. E. Schulz, “The refined structure of the complex between adenylate kinase from beef heart mitochondrial matrix and its substrate AMP at 1.85 Å resolution,” Journal of Molecular Biology, vol. 217, no. 3, pp. 541–549, 1991.

C. W. Müller and G. E. Schulz, “Structure of the complex between adenylate kinase from Escherichia coli and the inhibitor Ap₅A refined at 1.9 Å resolution. A model for a catalytic transition state,” Journal of Molecular Biology, vol. 224, no. 1, pp. 159–177, 1992.

R. E. Georgescu, E. G. Alexov, and M. R. Gunner, “Combining conformational flexibility and continuum electrostatics for calculating pKas in proteins,” Biophysical Journal, vol. 83, no. 4, pp. 1731–1748, 2002.

W. Rocchia, E. Alexov, and B. Honig, “Extending the applicability of the nonlinear Poisson-Boltzmann equation: multiple dielectric constants and multivalent ions,” Journal of Physical Chemistry B, vol. 105, no. 28, pp. 6507–6514, 2001.

D. A. Case, T. E. Cheatham III, T. Darden et al., “The Amber biomolecular simulation programs,” Journal of Computational Chemistry, vol. 26, no. 16, pp. 1668–1688, 2005.

J. W. Ponder and D. A. Case, “Force fields for protein simulations,” Advances in Protein Chemistry, vol. 66, pp. 27–85, 2003.

J. Wang, W. Wang, P. A. Kollman, and D. A. Case, “Antechamber, an accessory software package for molecular mechanical calculations,” Journal of Computational Chemistry, vol. 25, pp. 1157–1174, 2005.

J. F. Wang and K. C. Chou, “Insights from modeling the 3D structure of New Delhi metallo-β-lactamse and its binding interactions with antibiotic drugs,” PLoS ONE, vol. 6, no. 4, Article ID e18414, 2011.

P. Lian, D. Q. Wei, J. F. Wang, and K. C. Chou, “An allosteric mechanism inferred from molecular dynamics simulations on phospholamban pentamer in lipid membranes,” PLoS ONE, vol. 6, no. 4, Article ID e18587, 2011.

Y. Wang, D. Q. Wei, and J. F. Wang, “Molecular dynamics studies on T1 lipase: insight into a double-flap mechanism,” Journal of Chemical Information and Modeling, vol. 50, no. 5, pp. 875–878, 2010.

J. F. Wang, P. Hao, Y. X. Li, J. L. Dai, and X. Li, “Exploration of conformational transition in the aryl-binding site of human FXa using molecular dynamics simulations,” Journal of Molecular Modeling, vol. 18, no. 6, pp. 2717–2725, 2012.

J. Li, D. Q. Wei, J. F. Wang, and Y. X. Li, “A negative cooperativity mechanism of human CYP2E1 inferred from molecular dynamics simulations and free energy calculations,” Journal of Chemical Information and Modeling, vol. 51, no. 12, pp. 3217–3225, 2011.

J. F. Wang and K. C. Chou, “Insights into the mutation-induced HHH syndrome from modeling human mitochondrial ornithine transporter-1,” PLoS ONE, vol. 7, no. 1, Article ID e31048, 2012.

J. Li, D. Q. Wei, J. F. Wang, Z. T. Yu, and K. C. Chou, “Molecular dynamics simulations of CYP2E1,” Medicinal Chemistry, vol. 8, no. 2, pp. 208–221, 2012.

J. He, D. Q. Wei, J. F. Wang, and K. C. Chou, “Predicting protein-ligand binding sites based on an improved geometric algorithm,” Protein and Peptide Letters, vol. 18, no. 10, pp. 997–1001, 2011.

T. Dahnke, Z. Shi, H. Yan, R. T. Jiang, and M. D. Tsai, “Mechanism of adenylate kinase. Structural and functional roles of the conserved arginine-97 and arginine-132,” Biochemistry, vol. 31, no. 27, pp. 6318–6328, 1992.

M. A. Sinev, E. V. Sineva, V. Ittah, and E. Haas, “Domain closure in adenylate kinase,” Biochemistry, vol. 35, no. 20, pp. 6425–6437, 1996.

Y. E. Shapiro, M. A. Sinev, E. V. Sineva, V. Tugarinov, and E. Meirovitch, “Backbone dynamics of Escherichia coli adenylate kinase at the extreme stages of the catalytic cycle studied by ¹⁵N NMR relaxation,” Biochemistry, vol. 39, no. 22, pp. 6634–6644, 2000.

M. Wolf-Watz, V. Thai, K. Henzler-Wildman, G. Hadjipavlou, E. Z. Eisenmesser, and D. Kern, “Linkage between dynamics and catalysis in a thermophilic-mesophilic enzyme pair,” Nature Structural and Molecular Biology, vol. 11, no. 10, pp. 945–949, 2004.

Y. E. Shapiro and E. Meirovitch, “Activation energy of catalysis-related domain motion in E. coli adenylate kinase,” Journal of Physical Chemistry B, vol. 110, no. 23, pp. 11519–11524, 2006.

O. Miyashita, J. N. Onuchic, and P. G. Wolynes, “Nonlinear elasticity, proteinquakes, and the energy landscapes of functional transitions in proteins,” Proceedings of the National Academy of Sciences of the United States of America, vol. 100, no. 22, pp. 12570–12575, 2003.

P. C. Whitford, S. Gosavi, and J. N. Onuchic, “Conformational transitions in adenylate kinase: allosteric communication reduces misligation,” The Journal of Biological Chemistry, vol. 283, no. 4, pp. 2042–2048, 2008.

D. E. Koshland Jr., “Application of a theory of enzyme specificity to protein synthesis,” Proceedings of the National Academy of Sciences of the United States of America, vol. 44, no. 2, pp. 98–104, 1958.

B. Ma, S. Kumar, C. Tsai, and R. Nussinov, “Folding funnels and binding mechanisms,” Protein Engineering, vol. 12, no. 9, pp. 713–720, 1999.

C. Tsai, S. Kumar, B. Ma, and R. Nussinov, “Folding funnels, binding funnels, and protein function,” Protein Science, vol. 8, no. 6, pp. 1181–1190, 1999.

H. J. Zhang, X. R. Sheng, W. D. Niu, X. M. Pan, and J. M. Zhou, “Evidence for at least two native forms of rabbit muscle adenylate kinase in equilibrium in aqueous solution,” The Journal of Biological Chemistry, vol. 273, no. 13, pp. 7448–7456, 1998.

Y. Han, X. Li, and X. M. Pan, “Native states of adenylate kinase are two active sub-ensembles,” FEBS Letters, vol. 528, no. 1–3, pp. 161–165, 2002.

E. Z. Eisenmesser, O. Millet, W. Labeikovsky et al., “Intrinsic dynamics of an enzyme underlies catalysis,” Nature, vol. 438, no. 7064, pp. 117–121, 2005.

D. Tobi and I. Bahar, “Structural changes involved in protein binding correlate with intrinsic motions of proteins in the unbound state,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 52, pp. 18908–18913, 2005.

P. C. Whitford, O. Miyashita, Y. Levy, and J. N. Onuchic, “Conformational transitions of adenylate kinase: switching by cracking,” Journal of Molecular Biology, vol. 366, no. 5, pp. 1661–1671, 2007.