Collective Reaction Coordinate for Hybrid Quantum and Molecular Mechanics Simulations: A Case Study of the Hydride Transfer in Dihydrofolate Reductase.

The optimal description of the reaction coordinate in chemical systems is of great importance in simulating condensed phase reactions. In the current work, we present a collective reaction coordinate which is composed of several geometric coordinates which represent structural progress during the course of a hydride transfer reaction: the antisymmetric reactive stretch coordinate, the donor-acceptor distance (DAD) coordinate, and an orbital rehybridization coordinate. In this approach, the former coordinate serves as a distinguished reaction coordinate, while the latter two serve as environmental, Marcus-type inner-sphere reorganization coordinates. The classical free energy surface is obtained from multidimensional quantum mechanics-molecular mechanics (QM/MM) potential of mean force (PMF) simulations in conjunction with a general and efficient multidimensional weighted histogram method implementation. The minimum free energy path, or the collective reaction coordinate, connecting the dividing hypersurface to reactants and products, is obtained using an iterative scheme. In this approach, the string method is used to find the minimum free energy path. This path guides the multidimensional sampling, while the path is adaptively refined until convergence is achieved. As a model system, we choose the hydride transfer reaction in Escherichia coli dihydrofolate reductase (ecDHFR) using a recently developed accurate semiempirical potential energy surface. To estimate the advantages of the collective reaction coordinate, we perform activated dynamics simulations to obtain the reaction transmission coefficient. The results show that the combination of a distinguished reaction coordinate and an inner-sphere reorganization coordinate considerably reduces the dividing surface recrossing.

[1]  Darrin M York,et al.  An Efficient Linear-Scaling Ewald Method for Long-Range Electrostatic Interactions in Combined QM/MM Calculations. , 2005, Journal of chemical theory and computation.

[2]  S. Benkovic,et al.  Chemical basis for enzyme catalysis. , 2000, Biochemistry.

[3]  N. Goodey,et al.  The role of enzyme dynamics and tunnelling in catalysing hydride transfer: studies of distal mutants of dihydrofolate reductase , 2006, Philosophical Transactions of the Royal Society B: Biological Sciences.

[4]  Jiali Gao,et al.  Toward a Molecular Orbital Derived Empirical Potential for Liquid Simulations , 1997 .

[5]  M. Karplus,et al.  How Enzymes Work: Analysis by Modern Rate Theory and Computer Simulations , 2004, Science.

[6]  Alexander D. MacKerell,et al.  Development and current status of the CHARMM force field for nucleic acids , 2000, Biopolymers.

[7]  S. Benkovic,et al.  Effects of the donor-acceptor distance and dynamics on hydride tunneling in the dihydrofolate reductase catalyzed reaction. , 2012, Journal of the American Chemical Society.

[8]  G. Ciccotti,et al.  String method in collective variables: minimum free energy paths and isocommittor surfaces. , 2006, The Journal of chemical physics.

[9]  J. Kraut,et al.  Loop and subdomain movements in the mechanism of Escherichia coli dihydrofolate reductase: crystallographic evidence. , 1997, Biochemistry.

[10]  Dan Thomas Major,et al.  Momentum Distribution as a Fingerprint of Quantum Delocalization in Enzymatic Reactions: Open-Chain Path-Integral Simulations of Model Systems and the Hydride Transfer in Dihydrofolate Reductase. , 2012, Journal of chemical theory and computation.

[11]  D. Truhlar,et al.  The incorporation of quantum effects in enzyme kinetics modeling. , 2002, Accounts of chemical research.

[12]  J. Kästner.  Umbrella sampling , 2011 .

[13]  D. Truhlar,et al.  Mechanisms and free energies of enzymatic reactions. , 2006, Chemical reviews.

[14]  D. Truhlar,et al.  Small temperature dependence of the kinetic isotope effect for the hydride transfer reaction catalyzed by Escherichia coli dihydrofolate reductase. , 2005, The journal of physical chemistry. B.

[15]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[16]  David Chandler,et al.  Statistical mechanics of isomerization dynamics in liquids and the transition state approximation , 1978 .

[17]  D. York,et al.  Solvent polarization and kinetic isotope effects in nitroethane deprotonation and implications to the nitroalkane oxidase reaction. , 2005, Journal of the American Chemical Society.

[18]  Donald G Truhlar,et al.  Multidimensional tunneling, recrossing, and the transmission coefficient for enzymatic reactions. , 2006, Chemical reviews.

[19]  John E. Straub,et al.  Classical and modern methods in reaction rate theory , 1988 .

[20]  J. Kraut,et al.  Determination by Raman spectroscopy of the pKa of N5 of dihydrofolate bound to dihydrofolate reductase: mechanistic implications. , 1994, Biochemistry.

[21]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[22]  D. Chandler.  Roles of classical dynamics and quantum dynamics on activated processes occurring in liquids , 1986 .

[23]  Kent R. Wilson,et al.  Nonequilibrium solvation effects on reaction rates for model SN2 reactions in water , 1989 .

[24]  D. Truhlar,et al.  Nonperfect synchronization of reaction center rehybridization in the transition state of the hydride transfer catalyzed by dihydrofolate reductase. , 2005, Journal of the American Chemical Society.

[25]  Sharon Hammes-Schiffer,et al.  Nuclear Quantum Effects and Enzyme Dynamics in Dihydrofolate Reductase Catalysis , 2002 .

[26]  Jiali Gao,et al.  Hybrid Quantum and Molecular Mechanical Simulations: An Alternative Avenue to Solvent Effects in Organic Chemistry , 1996 .

[27]  Judith P Klinman,et al.  A 21st century revisionist's view at a turning point in enzymology. , 2009, Nature chemical biology.

[28]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[29]  R. Haddon,et al.  Hybridization as a metric for the reaction coordinate of the chemical reaction. Concert in chemical reactions , 1999 .

[30]  Kevin J. Naidoo,et al.  Implementation of an adaptive umbrella sampling method for the calculation of multidimensional potential of mean force of chemical reactions in solution , 2003, J. Comput. Chem..

[31]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[32]  K. Müller,et al.  Location of saddle points and minimum energy paths by a constrained simplex optimization procedure , 1979 .

[33]  Wei Wu,et al.  Density embedded VB/MM: a hybrid ab initio VB/MM with electrostatic embedding. , 2008, The journal of physical chemistry. A.

[34]  Henry Eyring,et al.  The Absolute Rate of Reactions in Condensed Phases , 1935 .

[35]  R. Callender,et al.  Structure of Dihydrofolate When Bound to Dihydrofolate Reductase , 1998 .

[36]  R. Marcus,et al.  Electron transfers in chemistry and biology , 1985 .

[37]  J. Kraut,et al.  pH-dependent conformational changes in Escherichia coli dihydrofolate reductase revealed by Raman difference spectroscopy. , 1997, Biophysical journal.

[38]  Y. Mo,et al.  An ab initio molecular orbital-valence bond (MOVB) method for simulating chemical reactions in solution , 2000 .

[39]  Arieh Warshel,et al.  The catalytic effect of dihydrofolate reductase and its mutants is determined by reorganization energies. , 2007, Biochemistry.

[40]  E. Vanden-Eijnden,et al.  Solvent coarse-graining and the string method applied to the hydrophobic collapse of a hydrated chain , 2007, Proceedings of the National Academy of Sciences.

[41]  H. C. Andersen.  Molecular dynamics simulations at constant pressure and/or temperature , 1980 .

[42]  Hao Hu,et al.  Elucidating solvent contributions to solution reactions with ab initio QM/MM methods. , 2010, The journal of physical chemistry. B.

[43]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

[44]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[45]  Bernard R. Brooks,et al.  Artificial reaction coordinate “tunneling” in free‐energy calculations: The catalytic reaction of RNase H , 2009, J. Comput. Chem..

[46]  Walter Thiel,et al.  QM/MM methods for biomolecular systems. , 2009, Angewandte Chemie.

[47]  R. C. Haddon,et al.  Hybridization as a Metric for the Reaction Coordinate of Chemical Reactions , 1998 .

[48]  S. Hammes‐Schiffer.  Quantum-classical simulation methods for hydrogen transfer in enzymes: a case study of dihydrofolate reductase. , 2004, Current opinion in structural biology.

[49]  H. Eyring.  The Activated Complex in Chemical Reactions , 1935 .

[50]  S. Benkovic,et al.  Effects of a distal mutation on active site chemistry. , 2006, Biochemistry.

[51]  P. Agarwal,et al.  Network of coupled promoting motions in enzyme catalysis , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[52]  G. Henkelman,et al.  A climbing image nudged elastic band method for finding saddle points and minimum energy paths , 2000 .

[53]  S. Benkovic,et al.  Construction and evaluation of the kinetic scheme associated with dihydrofolate reductase from Escherichia coli. , 1987, Biochemistry.

[54]  J. Klinman.  An integrated model for enzyme catalysis emerges from studies of hydrogen tunneling. , 2009, Chemical physics letters.

[55]  S. Benkovic,et al.  Tunneling and coupled motion in the Escherichia coli dihydrofolate reductase catalysis. , 2004, Journal of the American Chemical Society.

[56]  Eamonn F. Healy,et al.  Development and use of quantum mechanical molecular models. 76. AM1: a new general purpose quantum mechanical molecular model , 1985 .

[57]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[58]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[59]  J. Hynes,et al.  Theoretical aspects of tunneling proton transfer reactions in a polar environment , 2010 .

[60]  S. Benkovic,et al.  The effect of active-site isoleucine to alanine mutation on the DHFR catalyzed hydride-transfer. , 2010, Chemical communications.

[61]  H. Eyring,et al.  The Application of the Theory of Absolute Reacton Rates to Proteins. , 1939 .

[62]  Romelia Salomón-Ferrer,et al.  Dynamics and dissipation in enzyme catalysis , 2011, Proceedings of the National Academy of Sciences.

[63]  M Karplus,et al.  Dynamical theory of activated processes in globular proteins. , 1982, Proceedings of the National Academy of Sciences of the United States of America.

[64]  D. Truhlar,et al.  Reaction-path energetics and kinetics of the hydride transfer reaction catalyzed by dihydrofolate reductase. , 2003, Biochemistry.

[65]  Jiali Gao,et al.  Differential quantum tunneling contributions in nitroalkane oxidase catalyzed and the uncatalyzed proton transfer reaction , 2009, Proceedings of the National Academy of Sciences.

[66]  A. Warshel,et al.  Origin of the temperature dependence of isotope effects in enzymatic reactions: the case of dihydrofolate reductase. , 2007, The journal of physical chemistry. B.

[67]  G. Henkelman,et al.  Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points , 2000 .

[68]  Bruce J. Berne,et al.  Isomerization dynamics in liquids by molecular dynamics , 1980 .

[69]  Christoph Dellago,et al.  On the calculation of reaction rate constants in the transition path ensemble , 1999 .

[70]  C. Brooks,et al.  Barriers to Hydride Transfer in Wild Type and Mutant Dihydrofolate Reductase from E. coli , 2003 .

[71]  Martin Karplus,et al.  A POSITION DEPENDENT FRICTION MODEL FOR SOLUTION REACTIONS IN THE HIGH FRICTION REGIME : PROTON TRANSFER IN TRIOSEPHOSPHATE ISOMERASE (TIM) , 1996 .

[72]  M. Levitt,et al.  Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. , 1976, Journal of molecular biology.

[73]  E. Neria,et al.  Molecular dynamics of an enzyme reaction: proton transfer in TIM , 1997 .

[74]  Rudolph A. Marcus,et al.  On the Theory of Oxidation‐Reduction Reactions Involving Electron Transfer. I , 1956 .

[75]  Y. Mo,et al.  Ab initio QM/MM simulations with a molecular orbital‐valence bond (MOVB) method: application to an SN2 reaction in water , 2000 .

[76]  K. Hinsen,et al.  Potential of mean force and reaction rates for proton transfer in acetylacetone , 1997 .

[77]  Weiqing Ren,et al.  Higher Order String Method for Finding Minimum Energy Paths , 2003 .

[78]  Alexander D. MacKerell,et al.  Extending the treatment of backbone energetics in protein force fields: Limitations of gas‐phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations , 2004, J. Comput. Chem..

[79]  I. Benjamin,et al.  Effect of a Phase Transfer Catalyst on the Dynamics of an SN2 Reaction. A Molecular Dynamics Study , 2011 .

[80]  N. Goodey,et al.  Coordinated effects of distal mutations on environmentally coupled tunneling in dihydrofolate reductase , 2006, Proceedings of the National Academy of Sciences.

[81]  Ian F. Thorpe,et al.  Conformational substates modulate hydride transfer in dihydrofolate reductase. , 2005, Journal of the American Chemical Society.

[82]  A. Kohen,et al.  Elusive transition state of alcohol dehydrogenase unveiled , 2010, Proceedings of the National Academy of Sciences.

[83]  W. Thiel,et al.  Hybrid Quantum and Classical Simulations of the Dihydrofolate Reductase Catalyzed Hydride Transfer Reaction on an Accurate Semi-Empirical Potential Energy Surface. , 2011, Journal of chemical theory and computation.

[84]  John A. Montgomery,et al.  Trajectory analysis of a kinetic theory for isomerization dynamics in condensed phases , 1979 .

[85]  Arieh Warshel,et al.  Energetics and Dynamics of Enzymatic Reactions , 2001 .

[86]  B. C. Garrett,et al.  The role of collective solvent coordinates and nonequilibrium solvation in charge-transfer reactions , 2001 .

[87]  S. Benkovic,et al.  Coupling interactions of distal residues enhance dihydrofolate reductase catalysis: mutational effects on hydride transfer rates. , 2002, Biochemistry.

[88]  I. Tuñón,et al.  A Novel Strategy to Study Electrostatic Effects in Chemical Reactions: Differences between the Role of Solvent and the Active Site of Chalcone Isomerase in a Michael Addition. , 2012, Journal of chemical theory and computation.

[89]  Jiali Gao,et al.  Dynamics of an enzymatic substitution reaction in haloalkane dehalogenase. , 2004, Journal of the American Chemical Society.

[90]  Benoît Roux,et al.  Extension to the weighted histogram analysis method: combining umbrella sampling with free energy calculations , 2001 .

[91]  A. Warshel,et al.  Transition state theory can be used in studies of enzyme catalysis: lessons from simulations of tunnelling and dynamical effects in lipoxygenase and other systems , 2006, Philosophical Transactions of the Royal Society B: Biological Sciences.

[92]  Michael Hirsch,et al.  Bifurcation of reaction pathways: the set of valley ridge inflection points of a simple three-dimensional potential energy surface , 1998 .

[93]  S. Benkovic,et al.  A Perspective on Enzyme Catalysis , 2003, Science.

[94]  Scott F. Smith,et al.  SN2 reaction profiles in the gas phase and aqueous solution , 1984 .

[95]  Alessandro Laio,et al.  The energy gap as a universal reaction coordinate for the simulation of chemical reactions. , 2009, The journal of physical chemistry. B.

[96]  A. Chakraborty,et al.  A growing string method for determining transition states: comparison to the nudged elastic band and string methods. , 2004, The Journal of chemical physics.

[97]  Eric Vanden-Eijnden,et al.  Transition-path theory and path-finding algorithms for the study of rare events. , 2010, Annual review of physical chemistry.

[98]  Eric Vanden-Eijnden,et al.  Simplified and improved string method for computing the minimum energy paths in barrier-crossing events. , 2007, The Journal of chemical physics.

[99]  E. Vanden-Eijnden,et al.  String method for the study of rare events , 2002, cond-mat/0205527.

[100]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .