Increasing the time step with mass scaling in Born‐Oppenheimer ab initio QM/MM molecular dynamics simulations

Born‐Oppenheimer ab initio QM/MM molecular dynamics simulation with umbrella sampling is a state‐of‐the‐art approach to calculate free energy profiles of chemical reactions in complex systems. To further improve its computational efficiency, a mass‐scaling method with the increased time step in MD simulations has been explored and tested. It is found that by increasing the hydrogen mass to 10 amu, a time step of 3 fs can be employed in ab initio QM/MM MD simulations. In all our three test cases, including two solution reactions and one enzyme reaction, the resulted reaction free energy profiles with 3 fs time step and mass scaling are found to be in excellent agreement with the corresponding simulation results using 1 fs time step and the normal mass. These results indicate that for Born‐Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, the mass‐scaling method can significantly reduce its computational cost while has little effect on the calculated free energy profiles. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009

[1]  Alan M. Ferrenberg,et al.  New Monte Carlo technique for studying phase transitions. , 1988, Physical review letters.

[2]  Charles H. Bennett,et al.  Mass tensor molecular dynamics , 1975 .

[3]  A. Warshel,et al.  Structure/function correlations of proteins using MM, QM/MM, and related approaches: methods, concepts, pitfalls, and current progress. , 2003, Advances in protein chemistry.

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

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

[6]  Richard Wolfenden,et al.  Migration of Methyl Groups between Aliphatic Amines in Water [J. Am. Chem. Soc. 2003, 125, 310-311]. , 2003 .

[7]  Yingkai Zhang,et al.  How do SET-domain protein lysine methyltransferases achieve the methylation state specificity? Revisited by Ab initio QM/MM molecular dynamics simulations. , 2008, Journal of the American Chemical Society.

[8]  H. Senn,et al.  QM/MM Methods for Biological Systems , 2006 .

[9]  Yingkai Zhang,et al.  Pseudobond ab initio QM/MM approach and its applications to enzyme reactions , 2006 .

[10]  Tai-Sung Lee,et al.  A pseudobond approach to combining quantum mechanical and molecular mechanical methods , 1999 .

[11]  C. Brooks,et al.  Constant-temperature free energy surfaces for physical and chemical processes , 1993 .

[12]  Ian R. Gould,et al.  Electron-Transfer Reactions in the Marcus Inverted Region: Differences in Solvation and Electronic Coupling between Excited Charge-Transfer Complexes and Geminate Radical Ion Pairs , 1988 .

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

[14]  Adrian J Mulholland,et al.  Modelling enzyme reaction mechanisms, specificity and catalysis. , 2005, Drug discovery today.

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

[16]  W. Hase,et al.  Simulations of Gas-Phase Chemical Reactions: Applications to SN2 Nucleophilic Substitution , 1994, Science.

[17]  Dj McLennan,et al.  Semi-empirical calculation of rates of SN2 Finkelstein reactions in solution by a quasi-thermodynamic cycle , 1978 .

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

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

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

[21]  Urs Stocker,et al.  Increasing the Time Step and Efficiency of Molecular Dynamics Simulations: Optimal Solutions for Equilibrium Simulations or Structure Refinement of Large Biomolecules , 2003 .

[22]  J. Valleau,et al.  A Monte Carlo method for obtaining the interionic potential of mean force in ionic solution , 1975 .

[23]  D. Truhlar,et al.  Quantum mechanical methods for enzyme kinetics. , 2003, Annual review of physical chemistry.

[24]  James H. Hurley,et al.  Structure and Catalytic Mechanism of a SET Domain Protein Methyltransferase , 2002, Cell.

[25]  B. Bogdanov,et al.  Thermochemistry and structures of solvated SN2 complexes and transition states in the gas phase: experiment and theory , 2005 .

[26]  Frank Jensen,et al.  Steric Effects in SN2 Reactions. The Influence of Microsolvation , 2001 .

[27]  Thanh N. Truong,et al.  Hydration effects on reaction profiles: an ab initio dielectric continuum study of the SN2 Cl- + CH3Cl reaction , 1995 .

[28]  Berk Hess,et al.  Improving efficiency of large time‐scale molecular dynamics simulations of hydrogen‐rich systems , 1999, Journal of computational chemistry.

[29]  Craig,et al.  Gas-phase ionic reactions: dynamics and mechanism of nucleophilic displacements , 1998, Science.

[30]  Yingkai Zhang,et al.  Improved pseudobonds for combined ab initio quantum mechanical/molecular mechanical methods. , 2005, The Journal of chemical physics.

[31]  B. Roux.  The calculation of the potential of mean force using computer simulations , 1995 .

[32]  R. Friesner,et al.  Ab initio quantum chemical and mixed quantum mechanics/molecular mechanics (QM/MM) methods for studying enzymatic catalysis. , 2005, Annual review of physical chemistry.

[33]  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.

[34]  Shankar Kumar,et al.  Multidimensional free‐energy calculations using the weighted histogram analysis method , 1995, J. Comput. Chem..

[35]  Thanh N Truong,et al.  A coupled reference interaction site model/molecular dynamics study of the potential of mean force curve of the SN2 Cl- + CH3Cl reaction in water. , 2005, The journal of physical chemistry. B.

[36]  M. Karplus,et al.  A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations , 1990 .

[37]  Sheng-Yong Yang,et al.  Free Energy Profiles for the Identity SN2 Reactions Cl- + CH3Cl and NH3 + H3BNH3: A Constraint Ab Initio Molecular Dynamics Study , 2004 .

[38]  Martin Karplus,et al.  Catalysis and specificity in enzymes: a study of triosephosphate isomerase and comparison with methyl glyoxal synthase. , 2003, Advances in protein chemistry.

[39]  Yingkai Zhang,et al.  Ab initio quantum mechanical/molecular mechanical molecular dynamics simulation of enzyme catalysis: the case of histone lysine methyltransferase SET7/9. , 2007, The journal of physical chemistry. B.

[40]  A. Patrícia Bento,et al.  Ab initio and DFT benchmark study for nucleophilic substitution at carbon (SN2@C) and silicon (SN2@Si) , 2005, J. Comput. Chem..

[41]  Yingkai Zhang,et al.  Highly dissociative and concerted mechanism for the nicotinamide cleavage reaction in Sir2Tm enzyme suggested by ab initio QM/MM molecular dynamics simulations. , 2008, Journal of the American Chemical Society.

[42]  Yingkai Zhang,et al.  Determination of free energy profiles by repository based adaptive umbrella sampling: bridging nonequilibrium and quasiequilibrium simulations. , 2008, The Journal of chemical physics.

[43]  J. Andrew McCammon,et al.  Mass and step length optimization for the calculation of equilibrium properties by molecular dynamics simulation , 1990 .

[44]  Alan M. Ferrenberg,et al.  Optimized Monte Carlo data analysis. , 1989, Physical Review Letters.

[45]  U. Singh,et al.  A combined ab initio quantum mechanical and molecular mechanical method for carrying out simulations on complex molecular systems: Applications to the CH3Cl + Cl− exchange reaction and gas phase protonation of polyethers , 1986 .

[46]  C L Brooks,et al.  Calculations on folding of segment B1 of streptococcal protein G. , 1998, Journal of molecular biology.

[47]  Keiji Morokuma,et al.  POTENTIAL ENERGY SURFACE OF THE SN2 REACTION IN HYDRATED CLUSTERS , 1982 .

[48]  Hao Hu,et al.  Free energies of chemical reactions in solution and in enzymes with ab initio quantum mechanics/molecular mechanics methods. , 2008, Annual review of physical chemistry.