QM/MM minimum free-energy path: Methodology and application to triosephosphate isomerase

Journal of Chemical Theory and Computation

Volumn 3, Issue 2, 2007, Pages 390-406

  Hu, Hao ^a   Lu, Zhenyu ^a   Yang, Weitao ^a  

^a DUKE UNIVERSITY   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

EID: 34248224411     PISSN: 15499618     EISSN: None     Source Type: Journal    
DOI: 10.1021/ct600240y     Document Type: Article

References (92)

1. Zhang, Y.; Liu, H.; Yang, W. Free energy calculation on enzyme reactions with an efficient iterative procedure to determine minimum energy paths on a combined ab initio QM/MM potential energy surface. J. Chem. Phys. 2000, 112, 3483-3492.
2. Lu, Z.; Yang, W. Reaction path potential for complex systems derived from combined ab initio quantum mechanical and molecular mechanical calculations. J. Chem. Phys. 2004, 121, 89-100.
3. Warshel, A. Computer modeling of chemical reactions in enzymes and solutions; John Wiley & Sons: New York, 1991.
4. Fersht, A. Structure and mechanism in protein science; W. H. Freeman and Company: New York, 1998.
5. Garcia-Viloca, M.; Gao, J.; Karplus, M.; Truhlar, D. G. How enzymes work: Analysis by modern rate theory and computer simulations. Science 2004, 303, 186-195.
6. Shurki, A.; Warshel, A. Why does the Ras switch "break" by oncogenic mutations? Proteins: Struct., Funct., Genet. 2004, 55, 1-10.
7. Cui, Q.; Karplus, M. Triosephosphate isomerase: A theoretical comparison of alternative pathways. J. Am. Chem. Soc. 2001, 123, 2284-2290.
8. Cui, Q.; Elstner, M.; Karplus, M. A Theoretical Analysis of the Proton and Hydride Transfer in Liver Alcohol Dehydrogenase (LADH). J. Phys. Chem. B 2002, 106, 2721-2740.
9. Liu, H. Y.; Zhang, Y. K.; Yang, W. T. How is the active site of enolase organized to catalyze two different reaction steps? J. Am. Chem. Soc. 2000, 122, 6560-6570.
10. Cisneros, G. A.; Liu, H. Y.; Lu, Z. Y.; Yang, W. T. Reaction path determination for quantum mechanical/molecular mechanical modeling of enzyme reactions by combining first order and second order "chain-of- replicas" methods. J. Chem. Phys. 2005, 122.
11. Cisneros, G. A.; Liu, H. Y.; Zhang, Y. K.; Yang, W. T. Ab initio QM/MM study shows there is no general acid in the reaction catalyzed by 4-oxalocrotonate tautornerase. J. Am. Chem. Soc. 2003, 125, 10384-10393.
12. Cisneros, G. A.; Wang, M.; Silinski, P.; Fitzgerald, M. C.; Yang, W. The protein backbone makes important contributions to 4-oxalocrotonate tautomerase enzyme catalysis: Understanding from theory and experiment. Biochemistry 2004, 43, 6885-6892.
13. Cisneros, G. A.; Wang, M.; Silinski, P.; Fitzgerald, M. C.; Yang, W. T. Theoretical and experimental determination on two substrates turned over by 4-oxalocrotonate tautomerase. J. Phys. Chem. A 2006, 110, 700-708.
14. Warshel, A.; Levitt, M. Theoretic studies of enzymic reactions: Dielectric electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. J. Mol. Biol. 1976, 103, 227-249.
15. Field, M. J.; Bash, P. A.; Karplus, M. A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations. J. Comput. Chem. 1990, 11, 700-733.
16. Dewar, M. J. S.; Thiel, W. Ground states of molecules. 38. The MNDO method. Approximations and parameters. J. Am. Chem. Soc. 1977, 99, 4899-4507.
17. Dewar, M. J. S.; Zoebisch, E. G.; Healy, E. F.; Stewart, J. J. P. Development and use of quantum mechanical molecular models. 76. AM1: A new general purpose quantum mechanical molecular model. J. Am. Chem. Soc. 1985, 107, 3902-3909.
18. Stewart, J. J. P. Optimization of parameters for semiempirical methods I. Method. J. Comput. Chem. 1989, 10, 209-220.
19. Stewart, J. J. P. Optimization of parameters for semiempirical methods II. Applications. J. Comput. Chem. 1989, 10, 221-264.
20. Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, T.; Suhai, S.; Seifert, G. Self-consistent-charge density-functional tight-binding method for simulations of complex materials properties. Phys. Rev. B 1998, 58, 7260-7268.
21. Elstner, M.; Frauenheim, T.; Kaxiras, E.; Seifert, G.; Suhai, S. A self-consistent charge density-functional based tight-binding scheme for large biomolecules. Phys. Status Solidi B 2000, 217, 357-376.
22. Parr, R. G.; Yang, W. Density-Functional Theory of Atoms and Molecules; Oxford University Press: New York, 1994.
23. Hohenberg, P.; Kohn, W. Inhomogeneous electron gas. Phys. Rev. 1964, 136, B864-B871.
24. Kohn, W.; Sham, L. J. Self-consistent equations including exchange and correlation effects. Phys. Rev. 1965, 140, A1133-A1138.
25. Shurki, A.; Warshel, A. Structure/function correlations of proteins using MM, QMcharacterMM, and related approaches: Methods, concepts, pitfalls, and current progress. Adv. Protein Chem. 2003, 66, 249-313.
26. Warshel, A. Computer simulations of enzyme catalysis: Methods, progress, and insights. Annu. Rev. Biophys. Biomol. Struct. 2003, 32, 425-443.
27. Gao, J.; Truhlar, D. G. Quantum mechanical methods for enzyme kinetics. Annu. Rev. Phys. Chem. 2002, 53, 467-505.
28. Riccardi, D.; Schaefer, P.; Yang, Y.; Yu, H.; Ghosh, N.; Prat-Resina, X.; Konig, P.; Li, G.; Xu, D.; Guo, H.; Elstner, M.; Cui, Q. Development of effective quantum mechanical/molecular mechanical (QM/MM) methods for complex biological processes. J. Phys. Chem. B 2006, 110, 6458-6469.
29. Torrie, G. M.; Valleau, J. P. Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling. J. Comput. Phys. 1977, 23, 187-199.
30. N2 reaction involving chloride ion and methyl chloride in the gas phase and aqueous solution. J. Am. Chem. Soc. 1985, 107, 154-163.
31. Jorgensen, W. L. Free energy calculations: A breakthrough for modeling organic chemistry in solution. Acc. Chem. Res. 1989, 22.
32. Kuhn, B.; Kollman, P. A. QM-FE and molecular dynamics calculations on catechol O-methyltransferase: Free energy of activation in the enzyme and in aqueous solution and regioselectivity of the enzyme-catalyzed reaction. J. Am. Chem. Soc. 2000, 122, 2586-2596.
33. - exchange reaction and gas phase protonation of polyethers. J. Comput. Chem. 1986, 7, 718-730.
34. Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Montgomery, J. A., Jr.; Vreven, T.; Kudin, K. N.; Burant, J. C.; Millam, J. M.; Iyengar, S. S.; Tomasi, J.; Barone, V.; Mennucci, B.; Cossi, M.; Scalmani, G.; Rega, N.; Petersson, G. A.; Nakatsuji, H.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Klene, M.; Li, X.; Knox, J. E.; Hratchian, H. P.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Ayala, P. Y.; Morokuma, K.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Zakrzewski, V. G.; Dapprich, S.; Daniels, A. D.; Strain, M. C.; Farkas, O.; Malick, D. K.; Rabuck, A. D.; Raghavachari, K.; Foresman, J. B.; Ortiz, J. V.; Cui, Q.; Baboul, A. G.; Clifford, S.; Cioslowski, J.; Stefanov, B. B.; Liu, G.; Liashenko, A.; Piskorz, P.; Komaromi, I.; Martin, R. L.; Fox, D. J.; Keith, T.; Al-Laham, M. A.; Peng, C. Y.; Nanayakkara, A.; Challacombe, M.; Gill, P. M. W.; Johnson, B.; Chen, W.; Wong, M. W.; Gonzalez, C.; Pople, J. A. Gaussian 03, revision C.02; Gaussian, Inc.: Wallingford, CT, 2004.
35. Jorgensen, W. L.; Tirado-Rives, J. Molecular modeling of organic and biomolecular systems using BOSS and MCPRO. J. Comput. Chem. 2005, 26, 1689-1700.
36. Wesolowski, T.; Warshel, A. Ab initio free energy perturbation calculations of solvation free energy using the frozen density functional approach. J. Phys. Chem. 1994, 98, 5183-5187.
37. Wesolowski, T.; Muller, R. P.; Warshel, A. Ab initio frozen density functional calculations of proton transfer reactions in solution. J. Phys. Chem. 1996, 100, 15444-15449.
38. Muller, R. P.; Warshel, A. Ab initio calculations of free energy barriers for chemical reactions in solution. J. Phys. Chem. 1995, 99, 17516-17524.
39. Bentzien, J.; Muller, R. P.; Florián, J.; Warshel, A. Hybrid ab initio quantum mechanics/molecular mechanics calculations of free energy surfaces for enzymatic reactions: The nucleophilic attack in subtilisin. J. Phys. Chem. B 1998, 102, 2293-2301.
40. Strajbl, M.; Hong, G.; Warshel, A. Ab initio QM/MM simulation with proper sampling: "First principle" calculations of the free energy of the autodissociation of water in aqueous solution. J. Phys. Chem. B 2002, 106, 13333-13343.
41. Rod, T. H.; Ryde, U. Accurate QM/MM free energy calculations of enzyme reactions: Methylation by catechol O-methyltransferase. J. Chem. Theory. Comput. 2005, 1, 1240-1251.
42. Zhang, Y.; Lee, T.-S.; Yang, W. A pseudobond approach to combining quantum mechanical and molecular mechanical methods. J. Chem. Phys. 1999, 110, 46-54.
43. Liu, H.; Lu, Z.; Cisneros, G. A.; Yang, W. Parallel iterative reaction path optimization in ab initio quantum mechanical/molecular mechanical modeling of enzyme reactions. J. Chem. Phys. 2004, 121, 697-706.
44. Wang, M.; Lu, Z.; Yang, W. Transmission coefficient calculation for proton transfer in triosephosphate isomerase based on the reaction path potential method. J. Chem. Phys. 2004, 121, 101-107.
45. Wang, M.; Lu, Z.; Yang, W. Nuclear quantum effects on an enzyme-catalyzed reaction with reaction path potential: Proton transfer in triosephosphate isomerase. J. Chem. Phys. 2006, 124, 124516.
46. Kästner, J.; Senn, H. M.; Thiel, S.; Otte, N.; Thiel, W. QM/MM free-energy perturbation compared to thermodynamic integration and umbrella sampling: Application to an enzymatic reaction. J. Chem. Theory Comput. 2006, 2, 452-461.
47. Klahn, M.; Braun-Sand, S.; Rosta, E.; Warshel, A. On possible pitfalls in ab initio quantum mechanics/molecular mechanics minimization approaches for studies of enzymatic reactions. J. Phys. Chem. B 2005, 109, 15645-15650.
48. Schenter, G. K.; Garrett, B. C.; Truhlar, D. G. Generalized transition state theory in terms of the potential of mean force. J. Chem. Phys. 2003, 119, 5828-5833.
49. Fleurat-Lessard, P.; Ziegler, T. Tracing the minimum-energy path on the free-energy surface. J. Chem. Phys. 2005, 123, 084101.
50. Yang, S.-Y.; Hristov, I.; Fleurat-Lessard, P.; Ziegler, T. Optimizing the structures of minimum and transition state on the free energy surface. J. Phys. Chem. A 2005, 109, 197-204.
51. Hirao, H.; Nagae, Y.; Nagaoka, M. Transition-state optimization by the free energy gradient method: Application to aqueous-phase Menshutkin reaction between ammonia and methyl chloride. Chem. Phys. Lett. 2001, 348, 350-356.
52. Okuyama-Yoshida, N.; Nagaoka, M.; Yamabe, T. Potential energy function for intramolecular proton transfer reaction of glycine in aqueous solution. J. Phys. Chem. A 1998, 102, 285-292.
53. Okuyama-Yoshida, N.; Kataoka, K.; Nagaoka, M.; Yamabe, T. Structure optimization via free energy gradient method: Application to glycine zwitterion in aqueous solution. J. Chem. Phys. 2000, 113, 3519-3524.
54. Nagaoka, M.; Okuyama-Yoshida, N.; Yamabe, T. Origin of the transition state on the free energy surface: Intramolecular proton transfer reaction of glycine in aqueous solution. J. Phys. Chem. A 1998, 102, 8202-8208.
55. Okuyama-Yoshida, N.; Nagaoka, M.; Yamabe, T. Transition-state optimization on free energy surface: Toward solution chemical reaction ergodography. Int. J. Quantum Chem. 1998, 70, 95-103.
56. Jónsson, H.; Mills, G.; Jacobsen, K. W. Nudged elastic band method for finding minimum energy paths of transitions; World Scientific: Singapore, 1998.
57. Ayala, P. Y.; Schlegel, H. B. A combined method for determining reaction paths, minima, and transition state geometries. J. Chem. Phys. 1997, 107, 375-384.
58. Xie, L.; Liu, H.; Yang, W. Adapting the nudged elastic band method for determining minimum-energy paths of chemical reactions in enzymes. J. Chem. Phys. 2004, 120, 8039-8052.
59. Burger, S. K.; Yang, W. Quadratic string method for determining the minimum-energy path based on multiobjective optimization. J. Chem. Phys. 2006, 124, 054109.
60. Singh, U. C.; Kollman, P. A. An approach to computing electrostatic charges for molecules. J. Comput. Chem. 1984, 5, 129-145.
61. Besler, B. H.; Merz, K. M., Jr.; Kollman, P. A. Atomic charges derived from semiempirical methods. J. Comput. Chem. 1990, 11, 431-439.
62. Bolhuis, P. G.; Dellago, C.; Chandler, D. Reaction coordinates of biomolecular isomerization. Proc. Natl. Acad. Sci. U.S.A. 2000, 97, 5877-5882.
63. Becke, A. D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648-5652.
64. Lee, C.; Yang, W.; Parr, R. G. Development of the Colle-Salvetti correlation energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785.
65. Miller, W. H.; Handy, N. C.; Adams, J. E. Reaction path Hamiltonian for polyatomic molecules. J. Chem. Phys. 1980, 72, 99-112.
66. Morita, A.; Kato, S. Ab initio molecular orbital theory on intramolecular charge polarization: Effect of hydrogen abstraction on the charge sensitivity of aromatic and non-aromatic species. J. Am. Chem. Soc. 1997, 119, 4021-4032.
67. Cornell, W. D.; Cieplak, P.; Bayly, C.; Gould, I. R.; Merz, K. M. J.; Ferguson, D. M.; Spellmeyer, D. C.; Fox, T.; Caldwell, J. W.; Kollman, P. A. A second generation force field for the simulation of proteins and nucleic acids. J. Am. Chem. Soc. 1995, 117, 5179-5197.
68. Ponder, J. W. TINKER, Software Tools for Molecular Design, version 3.8; Washington University: St. Louis, MO, 2000.
69. Tuckerman, M. E.; Berne, B. J.; Martyna, G. J. Reversible multiple time scale molecular dynamics. J. Chem. Phys. 1992, 97, 1990-2001.
70. Schlick, T.; Skeel, R. D.; Brünger, A. T.; Kalé, L. V.; Board, J. A.; Hermans, J.; Schulten, K. Algorithmic challenges in computational molecular biophysics. J. Comput. Phys. 1999, 151, 9-48.
71. Berendsen, H. J. C.; Postma, J. P. M.; van Gunsteren, W. F.; DiNola, A.; Haak, J. R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684-3690.
72. Mann, G.; Yun, R. H.; Nyland, L.; Prins, J.; Board, J.; Hermans, J. The Sigma MD program and a generic interface applicable to multi-functional programs with complex, hierarchical command structure. In Computational Methods for Macromolecules: Challenges and Applications; Proceedings of the 3rd International Workshop on Algorithms for Macromolecular Modelling; Schlick, T., Gan, H. H., Eds.; Springer-Verlag: Berlin, 2002; pp 129-145.
73. Hu, H.; Elstner, M.; Hermans, J. Comparison of a QM/MM force field and molecular mechanics force fields in simulations of alanine and glycine dipeptides (Ace-Ala-Nme and Ace-Gly-Nme) in water in relation to the problem of modeling the unfolded peptide backbone in solution. Proteins: Struct., Funct., Genet. 2003, 50, 451-463.
74. Hu, H.; Yang, W. Dual-topology/dual-coordinate free-energy simulation using QM/MM force field. J. Chem. Phys. 2005, 123, 041102.
75. Ryckaert, J. P.; Ciccotti, G.; Berendsen, H. J. C. Numerical integration of the Cartesian equations of motion of a system with constraints: Molecular dynamics of n-alkanes. J. Comput. Phys. 1977, 23, 327-341.
76. Beveridge, D. L.; DiCapua, F. M. Free energy via molecular simulatin: A primer. In Computer Simulation of Biomolecular Systems; van Gunsteren, W. F., Weiner, P. K., Eds. ESCOM: Leiden, The Netherlands, 1989; Vol. 1, pp 1-26.
77. Hermans, J. A simple analysis of noise and hysteresis in free energy simulations. J. Phys. Chem. 1991, 95, 9029-9032.
78. Hu, H.; Yun, R. H.; Hermans, J. Reversibility of free energy simulations: Slow growth may have a unique advantage. (With a note on use of Ewald summation). Mol. Simul. 2002, 28, 67-80.
79. Li, G.; Cui, Q. Direct determination of reaction paths and stationary points on potential of mean force surfaces. J. Mol. Graphics Modell. 2005, 24, 82-93.
80. Chandler, D. Introduction to Modern Statistical Mechanics; Oxford University Press: New York, 1986.
81. Zwanzig, R. W. High-temperature equation of state by a perturbation method. I. Nonpolar gases. J. Chem. Phys. 1954, 22, 1420-1426.
82. Carter, E. A.; Ciccotti, G.; Hynes, J. T.; Kapral, R. Constrained reaction coordinate dynamics for the simulation of rare events. Chem. Phys. Lett. 1989, 156, 472-477.
83. Ryckaert, J. P.; Ciccotti, G. Introduction of Andersen's demon in the molecular dynamics of systems with constraints. J. Chem. Phys. 1983, 78, 7368-7374.
84. Sprik, M.; Ciccotti, G. Free energy from constrained molecular dynamics. J. Chem. Phys. 1998, 109, 7737-7744.
85. Berendsen, H. J. C.; Postma, J. P. M.; van Gunsteren, W. F. Statistical mechanics and molecular dynamics: The calculation of free energy. In Molecular Dynamics and Protein Structure; Hermans, J., Ed. Polycrystal Book Service: Western Springs, IL, 1985; pp 43-46.
86. Straatsma, T. P.; Zacharias, M.; McCammon, J. A. Holonomic constraint contributions to energy differences from thermodynamic integration molecular dynamics simulations. Chem. Phys. Lett. 1992, 196, 297-302.
87. de Otter, W. K.; Briels, W. J. The calculation of free-energy differences by constrained molecular-dynamics simulations. J. Chem. Phys. 1998, 109, 4139-4146.
88. Schlitter, J.; Klähn, M. A new concise expression for the free energy of a reaction coordinate. J. Chem. Phys. 2003, 118, 2057-2060.
89. Antikainen, N. M.; Smiley, R. D.; Benkovic, S. J.; Hammes, G. G. Conformation coupled enzyme catalysis: Single-molecule and transient kinetics investigation of dihydrofolate reductase. Biochemistry 2005, 44, 16835-16843.
90. Eisenmesser, E. Z.; Bosco, D. A.; Akke, M.; Kern, D. Enzyme dynamics during catalysis. Science 2002, 295, 1520-1523.
91. Amadei, A.; de Groot, B. L.; Ceruso, M. A.; Paci, M.; Di Nola, A.; Berendsen, H. J. C. A kinetic model for the internal motions of proteins: Diffusion between multiple harmonic wells. Proteins: Struct., Funct., Genet. 1999, 35, 283-292.
92. Kitao, A.; Hayward, S.; Go, N. Energy landscape of a native protein: Jumping-among-minima model. Proteins: Struct., Funct., Genet. 1999, 33, 496-517.