Efficient parallel implementations of QM/MM-REMD (quantum mechanical/molecular mechanics-replica-exchange MD) and umbrella sampling: isomerization of H2O2 in aqueous solution.

An efficient parallel implementation of QM/MM-based replica-exchange molecular dynamics (REMD) as well as umbrella samplings techniques was proposed by adopting the generalized distributed data interface (GDDI). Parallelization speed-up of 40.5 on 48 cores was achieved, making our QM/MM-MD engine a robust tool for studying complex chemical dynamics in solution. They were comparatively used to study the torsional isomerization of hydrogen peroxide in aqueous solution. All results by QM/MM-REMD and QM/MM umbrella sampling techniques yielded nearly identical potentials of mean force (PMFs) regardless of the particular QM theories for solute, showing that the overall dynamics are mainly determined by solvation. Although the entropic penalty of solvent rearrangements exists in cisoid conformers, it was found that both strong intermolecular hydrogen bonding and dipole-dipole interactions preferentially stabilize them in solution, reducing the torsional free-energy barrier at 0° by about 3 kcal/mol as compared to that in gas phase.

[1]  M. Gordon,et al.  The effective fragment potential: small clusters and radial distribution functions. , 2004, The Journal of chemical physics.

[2]  K. Hukushima,et al.  Exchange Monte Carlo Method and Application to Spin Glass Simulations , 1995, cond-mat/9512035.

[3]  Heming Xiao,et al.  Theoretical study on intermolecular interactions and thermodynamic properties of water–hydrogen peroxide clusters , 2003 .

[4]  L. Slipchenko,et al.  Solvent effects on the electronic transitions of p-nitroaniline: a QM/EFP study. , 2011, The journal of physical chemistry. A.

[5]  Jiali Gao,et al.  A molecular-orbital derived polarization potential for liquid water , 1998 .

[6]  Yuichi Inadomi,et al.  Fragment molecular orbital method: application to molecular dynamics simulation, ‘ab initio FMO-MD’ , 2003 .

[7]  J. Stone,et al.  Hydrogen peroxide: a signaling messenger. , 2006, Antioxidants & redox signaling.

[8]  Donald G Truhlar,et al.  X-Pol Potential: An Electronic Structure-Based Force Field for Molecular Dynamics Simulation of a Solvated Protein in Water. , 2009, Journal of chemical theory and computation.

[9]  D. Truhlar,et al.  QM/MM: what have we learned, where are we, and where do we go from here? , 2007 .

[10]  M. Ruiz‐López,et al.  Molecular dynamics of hydrogen peroxide in liquid water using a combined quantum/classical force field , 2007 .

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

[12]  Spencer R Pruitt,et al.  Fragmentation methods: a route to accurate calculations on large systems. , 2012, Chemical reviews.

[13]  Mark S Gordon,et al.  Fast fragments: The development of a parallel effective fragment potential method , 2004, J. Comput. Chem..

[14]  Michele Parrinello,et al.  Quickstep: Fast and accurate density functional calculations using a mixed Gaussian and plane waves approach , 2005, Comput. Phys. Commun..

[15]  N. Handy,et al.  Potential Energy Surface and Vibrational−Rotational Energy Levels of Hydrogen Peroxide , 1998 .

[16]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[17]  Michael D. Hands,et al.  Intermolecular interactions in complex liquids: effective fragment potential investigation of water-tert-butanol mixtures. , 2012, The journal of physical chemistry. B.

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

[19]  Wang,et al.  Replica Monte Carlo simulation of spin glasses. , 1986, Physical review letters.

[20]  Masato Kobayashi,et al.  Two‐level hierarchical parallelization of second‐order Møller–plesset perturbation calculations in divide‐and‐conquer method , 2011, J. Comput. Chem..

[21]  L. Bartolotti,et al.  Structures, energetics, and vibrational spectra of H2O2...(H2O)n, n = 1-6 clusters: Ab initio quantum chemical investigations. , 2005, The journal of physical chemistry. A.

[22]  Car,et al.  Unified approach for molecular dynamics and density-functional theory. , 1985, Physical review letters.

[23]  David van der Spoel,et al.  Protein folding kinetics and thermodynamics from atomistic simulations. , 2006, Physical review letters.

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

[25]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[26]  Ivan S Ufimtsev,et al.  Quantum Chemistry on Graphical Processing Units. 3. Analytical Energy Gradients, Geometry Optimization, and First Principles Molecular Dynamics. , 2009, Journal of chemical theory and computation.

[27]  M. Cho,et al.  Direct simulations of anharmonic infrared spectra using quantum mechanical/effective fragment potential molecular dynamics (QM/EFP-MD): methanol in water. , 2012, The journal of physical chemistry. A.

[28]  Quantum mechanical/effective fragment potential molecular dynamics (QM/EFP-MD) study on intra-molecular proton transfer of glycine in water , 2012 .

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

[30]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[31]  E. Glendening,et al.  An intrinsic reaction coordinate calculation of the torsion-internal rotation potential of hydrogen peroxide and its isotopomers. , 2004, The Journal of chemical physics.

[32]  Y. Sugita,et al.  Interionic hydration structures of NaCl in aqueous solution: a combined study of quantum mechanical cluster calculations and QM/EFP-MD simulations. , 2013, The journal of physical chemistry. B.

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

[34]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[35]  Kazuo Kitaura,et al.  Exploring chemistry with the fragment molecular orbital method. , 2012, Physical chemistry chemical physics : PCCP.

[36]  Yirong Mo,et al.  Origins of Rotational Barriers in Hydrogen Peroxide and Hydrazine. , 2005, Journal of chemical theory and computation.

[37]  N. Uddin,et al.  Hydrophobic and hydrophilic associations of a methanol pair in aqueous solution. , 2012, The journal of physical chemistry. B.

[38]  Mark S. Gordon,et al.  A new hierarchical parallelization scheme: Generalized distributed data interface (GDDI), and an application to the fragment molecular orbital method (FMO) , 2004, J. Comput. Chem..

[39]  J. W. C. Johns,et al.  The far infrared spectrum of H2O2. First observation of the staggering of the levels and determination of the cis barrier , 1989 .

[40]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[41]  K. Kitaura,et al.  Fragment molecular orbital method: an approximate computational method for large molecules , 1999 .

[42]  Sanford Weisberg,et al.  Computing science and statistics : proceedings of the 30th Symposium on the Interface, Minneapolis, Minnesota, May 13-16, 1998 : dimension reduction, computational complexity and information , 1998 .

[43]  Michael W. Mahoney,et al.  A five-site model for liquid water and the reproduction of the density anomaly by rigid, nonpolarizable potential functions , 2000 .