Acceleration of Ab Initio QM/MM Calculations under Periodic Boundary Conditions by Multiscale and Multiple Time Step Approaches.

Development of multiscale ab initio quantum mechanical and molecular mechanical (AI-QM/MM) method for periodic boundary molecular dynamics (MD) simulations and their acceleration by multiple time step approach are described. The developed method achieves accuracy and efficiency by integrating the AI-QM/MM level of theory and the previously developed semiempirical (SE) QM/MM-Ewald sum method [J. Chem. Theory Comput. 2005, 1, 2] extended to the smooth particle-mesh Ewald (PME) summation method. In the developed methods, the total energy of the simulated system is evaluated at the SE-QM/MM-PME level of theory to include long-range QM/MM electrostatic interactions, which is then corrected on the fly using the AI-QM/MM level of theory within the real space cutoff. The resulting energy expression enables decomposition of total forces applied to each atom into forces determined at the low-level SE-QM/MM method and correction forces at the AI-QM/MM level, to integrate the system using the reversible reference system propagator algorithm. The resulting method achieves a substantial speed-up of the entire calculation by minimizing the number of time-consuming energy and gradient evaluations at the AI-QM/MM level. Test calculations show that the developed multiple time step AI-QM/MM method yields MD trajectories and potential of mean force profiles comparable to single time step QM/MM results. The developed method, together with message passing interface (MPI) parallelization, accelerates the present AI-QM/MM MD simulations about 30-fold relative to the speed of single-core AI-QM/MM simulations for the molecular systems tested in the present work, making the method less than one order slower than the SE-QM/MM methods under periodic boundary conditions.

[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]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[3]  György G. Ferenczy,et al.  Semiempirical AM1 electrostatic potentials and AM1 electrostatic potential derived charges: A comparison with ab initio values , 1989 .

[4]  M. Head‐Gordon,et al.  Curvy-steps approach to constraint-free extended-Lagrangian ab initio molecular dynamics, using atom-centered basis functions: convergence toward Born-Oppenheimer trajectories. , 2004, The Journal of chemical physics.

[5]  Weitao Yang,et al.  QM/MM Minimum Free Energy Path: Methodology and Application to Triosephosphate Isomerase. , 2007, Journal of chemical theory and computation.

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

[7]  Ian J. Bush,et al.  The GAMESS-UK electronic structure package: algorithms, developments and applications , 2005 .

[8]  Thom Vreven,et al.  The ONIOM (our own N-layered integrated molecular orbital + molecular mechanics) method for the first singlet excited (S1) state photoisomerization path of a retinal protonated Schiff base , 2000 .

[9]  W. Thiel,et al.  Semiempirical treatment of electrostatic potentials and partial charges in combined quantum mechanical and molecular mechanical approaches , 1996, J. Comput. Chem..

[10]  T. Darden,et al.  Molecular dynamics simulations of biomolecules: long-range electrostatic effects. , 1999, Annual review of biophysics and biomolecular structure.

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

[12]  Jiali Gao,et al.  Methods and Applications of Combined Quantum Mechanical and Molecular Mechanical Potentials , 2007 .

[13]  J. Perram,et al.  Simulation of electrostatic systems in periodic boundary conditions. I. Lattice sums and dielectric constants , 1980, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[14]  K. Ranaghan,et al.  Protein dynamics and enzyme catalysis: insights from simulations. , 2011, Biochimica et biophysica acta.

[15]  Sergio Martí,et al.  Computing kinetic isotope effects for chorismate mutase with high accuracy. A new DFT/MM strategy. , 2005, The journal of physical chemistry. B.

[16]  A. Warshel,et al.  Electrostatic basis for enzyme catalysis. , 2006, Chemical reviews.

[17]  Sándor Suhai,et al.  Self-consistent-charge density-functional tight-binding method for simulations of complex materials properties , 1998 .

[18]  G. Scuseria,et al.  Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals , 2001 .

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

[20]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

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

[22]  Jerry M. Parks,et al.  Quantum mechanics/molecular mechanics minimum free-energy path for accurate reaction energetics in solution and enzymes: sequential sampling and optimization on the potential of mean force surface. , 2008, The Journal of chemical physics.

[23]  Sándor Suhai,et al.  A Self‐Consistent Charge Density‐Functional Based Tight‐Binding Method for Predictive Materials Simulations in Physics, Chemistry and Biology , 2000 .

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

[25]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[26]  John M Herbert,et al.  Accelerated, energy-conserving Born-Oppenheimer molecular dynamics via Fock matrix extrapolation. , 2005, Physical chemistry chemical physics : PCCP.

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

[28]  D. Rinaldi,et al.  Long-range electrostatic interactions in hybrid quantum and molecular mechanical dynamics using a lattice summation approach. , 2005, The Journal of chemical physics.

[29]  C. Breneman,et al.  Determining atom‐centered monopoles from molecular electrostatic potentials. The need for high sampling density in formamide conformational analysis , 1990 .

[30]  Mark E. Tuckerman,et al.  Explicit reversible integrators for extended systems dynamics , 1996 .

[31]  Martin Karplus,et al.  Lagrangian formulation with dissipation of Born-Oppenheimer molecular dynamics using the density-functional tight-binding method. , 2011, The Journal of chemical physics.

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

[33]  T. Darden,et al.  The effect of long‐range electrostatic interactions in simulations of macromolecular crystals: A comparison of the Ewald and truncated list methods , 1993 .

[34]  J. Herbert,et al.  Periodic boundary conditions for QM/MM calculations: Ewald summation for extended Gaussian basis sets. , 2013, The Journal of chemical physics.

[35]  R. S. Mulliken.  Electronic Population Analysis on LCAO–MO Molecular Wave Functions. I , 1955 .

[36]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[37]  Ryan P Steele,et al.  Communication: Multiple-timestep ab initio molecular dynamics with electron correlation. , 2013, The Journal of chemical physics.

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

[39]  Efthimios Kaxiras,et al.  A QM/MM Implementation of the Self-Consistent Charge Density Functional Tight Binding (SCC-DFTB) Method , 2001 .

[40]  Peter Pulay,et al.  Fock matrix dynamics , 2004 .

[41]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[42]  S. Broyde,et al.  DNA cytosine methylation: structural and thermodynamic characterization of the epigenetic marking mechanism. , 2013, Biochemistry.

[43]  J. Herbert,et al.  An efficient, fragment-based electronic structure method for molecular systems: self-consistent polarization with perturbative two-body exchange and dispersion. , 2011, The Journal of chemical physics.

[44]  Kwangho Nam,et al.  Acceleration of Semiempirical Quantum Mechanical Calculations by Extended Lagrangian Molecular Dynamics Approach. , 2013, Journal of chemical theory and computation.

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

[46]  J. Andrew McCammon,et al.  Influence of Structural Fluctuation on Enzyme Reaction Energy Barriers in Combined Quantum Mechanical/Molecular Mechanical Studies , 2003 .

[47]  Ross C. Walker,et al.  The implementation of a fast and accurate QM/MM potential method in Amber , 2008, J. Comput. Chem..

[48]  Jiali Gao,et al.  A hybrid semiempirical quantum mechanical and lattice-sum method for electrostatic interactions in fluid simulations , 1997 .

[49]  Q. Cui,et al.  pKa calculations in solution and proteins with QM/MM free energy perturbation simulations: a quantitative test of QM/MM protocols. , 2005, The journal of physical chemistry. B.

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

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

[52]  Thomas E Markland,et al.  Multiple time step integrators in ab initio molecular dynamics. , 2014, The Journal of chemical physics.

[53]  J. Stewart.  Optimization of parameters for semiempirical methods I. Method , 1989 .

[54]  Tom Welton,et al.  Solvents and Solvent Effects in Organic Chemistry: REICHARDT:SOLV.EFF. 4ED O-BK , 2010 .

[55]  Elena F. Sheka,et al.  NANOPACK: Parallel codes for semiempirical quantum chemical calculations of large systems in the sp‐ and spd‐basis , 2002 .

[56]  R. Giernoth.  Solvents and Solvent Effects in Organic Chemistry. 4th Ed. By Christian Reichardt and Thomas Welton. , 2011 .

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