Generalised canonical–isokinetic ensemble: speeding up multiscale molecular dynamics and coupling with 3D molecular theory of solvation

We have proposed a new canonical–isokinetic ensemble for efficient sampling of conformational space in molecular dynamics (MD) simulations which leads to the optimised isokinetic Nosé–Hoover (OIN) chain algorithm for atomic and molecular systems. We applied OIN to multiple time step (MTS) MD simulations of the rigid and flexible models of water to demonstrate its advantage over the standard canonical, isokinetic and canonical–isokinetic ensembles. With the stabilising effect of OIN thermostatting in MTS-MD, gigantic outer time steps up to picoseconds can be employed to accurately calculate equilibrium and conformational properties. Furthermore, we developed the atomic version of OIN for MTS-MD of a biomolecule in a solvent potential of mean force obtained at sequential MD steps by using the molecular theory of solvation, aka three-dimensional reference interaction site model with the Kovalenko–Hirata closure (3D-RISM-KH). The solvation forces are obtained analytically by converging the 3D-RISM-KH integral equations once per several OIN outer time steps, and are calculated in between by using solvation force-coordinate extrapolation (SFCE) in the subspace of previous successive solutions to 3D-RISM-KH. For illustration, we applied the multiscale OIN/SFCE/3D-RISM-KH algorithm to a fully flexible model of alanine dipeptide in aqueous solution. Although the computational rate of solvent sampling in OIN/SFCE/3D-RISM-KH is already 20 times faster than standard MD with explicit solvent, further substantial acceleration of sampling stems from making solute evolution steps in a statistically averaged potential of mean force obtained from 3D-RISM-KH. The latter efficiently samples the phase space for essential events with rare statistics such as exchange and localisation of solvent and ligand molecules in confined spaces, pockets and at binding sites of the solute macromolecule, as distinct from MD with explicit solvent which requires enormous computational time and number of steps in such cases.

[1]  I. Omelyan.  Processed splitting algorithms for rigid-body molecular dynamics simulations. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[2]  Michael Feig,et al.  Is Alanine Dipeptide a Good Model for Representing the Torsional Preferences of Protein Backbones? , 2008, Journal of chemical theory and computation.

[3]  Andriy Kovalenko,et al.  Calculation of local water densities in biological systems: a comparison of molecular dynamics simulations and the 3D-RISM-KH molecular theory of solvation. , 2011, The journal of physical chemistry. B.

[4]  Giovanni Ciccotti,et al.  Computer Simulations in Condensed Matter: Systems: From Materials to Chemical Biology. Volume 1 , 2006 .

[5]  Ruhong Zhou,et al.  Molecular Dynamics with Multiple Time Scales: How to Avoid Pitfalls. , 2010, Journal of chemical theory and computation.

[6]  Jesús A. Izaguirre,et al.  Verlet-I/R-RESPA/Impulse is Limited by Nonlinear Instabilities , 2003, SIAM J. Sci. Comput..

[7]  A. Kovalenko,et al.  Evaluation of the SCF Combination of KS-DFT and 3D-RISM-KH; Solvation Effect on Conformational Equilibria, Tautomerization Energies, and Activation Barriers. , 2007, Journal of chemical theory and computation.

[8]  R Folk,et al.  Optimized Verlet-like algorithms for molecular dynamics simulations. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[9]  J. Berg,et al.  Molecular dynamics simulations of biomolecules , 2002, Nature Structural Biology.

[10]  P. Pulay.  Convergence acceleration of iterative sequences. the case of scf iteration , 1980 .

[11]  A. Kovalenko,et al.  Self-consistent combination of the three-dimensional RISM theory of molecular solvation with analytical gradients and the Amsterdam density functional package , 2005, 2005 International Conference on MEMS,NANO and Smart Systems.

[12]  T. Schlick,et al.  Overcoming stability limitations in biomolecular dynamics. I. Combining force splitting via extrapolation with Langevin dynamics in LN , 1998 .

[13]  H. Fenniri,et al.  Helical rosette nanotubes with tunable stability and hierarchy. , 2005, Journal of the American Chemical Society.

[14]  A. Kovalenko,et al.  Multiple time scale molecular dynamics for fluids with orientational degrees of freedom. I. Microcanonical ensemble. , 2011, The Journal of chemical physics.

[15]  Charles L. Brooks,et al.  Conformational equilibrium in the alanine dipeptide in the gas phase and aqueous solution : a comparison of theoretical results , 1992 .

[16]  Minhaeng Cho,et al.  Classical and quantum mechanical/molecular mechanical molecular dynamics simulations of alanine dipeptide in water: comparisons with IR and vibrational circular dichroism spectra. , 2008, The Journal of chemical physics.

[17]  A. Kovalenko,et al.  Microtubule stability studied by three-dimensional molecular theory of solvation. , 2007, Biophysical journal.

[18]  Berend Smit,et al.  Understanding molecular simulation: from algorithms to applications , 1996 .

[19]  Andriy Kovalenko,et al.  3D-RISM-KH approach for biomolecular modelling at nanoscale: thermodynamics of fibril formation and beyond , 2010 .

[20]  T Schlick,et al.  Biomolecular dynamics at long timesteps: bridging the timescale gap between simulation and experimentation. , 1997, Annual review of biophysics and biomolecular structure.

[21]  M. Tuckerman.  Statistical Mechanics: Theory and Molecular Simulation , 2010 .

[22]  F. Hirata,et al.  Theoretical study for volume changes associated with the helix-coil transition of peptides. , 2001, Biopolymers.

[23]  Carlos Simmerling,et al.  Three-dimensional molecular theory of solvation coupled with molecular dynamics in Amber. , 2010, Journal of chemical theory and computation.

[24]  H. Fenniri,et al.  Molecular basis for water-promoted supramolecular chirality inversion in helical rosette nanotubes. , 2007, Journal of the American Chemical Society.

[25]  Fumio Hirata,et al.  Combination of molecular dynamics method and 3D‐RISM theory for conformational sampling of large flexible molecules in solution , 2008, J. Comput. Chem..

[26]  Jesús A. Izaguirre,et al.  An impulse integrator for Langevin dynamics , 2002 .

[27]  R. Skeel,et al.  Langevin stabilization of molecular dynamics , 2001 .

[28]  S. Gusarov,et al.  Multiscale Modeling of the Adsorption Interaction between Model Bitumen Compounds and Zeolite Nanoparticles in Gas and Liquid Phases , 2011 .

[29]  B. Montgomery Pettitt,et al.  A site-site theory for finite concentration saline solutions , 1992 .

[30]  I. P. Omelyan,et al.  Advanced multiple time scale molecular dynamics , 2009 .

[31]  David A. Dixon,et al.  Conformational Study of the Alanine Dipeptide at the MP2 and DFT Levels , 2002 .

[32]  Benoît Roux,et al.  An Integral Equation To Describe the Solvation of Polar Molecules in Liquid Water , 1997 .

[33]  Mark E. Tuckerman,et al.  Algorithms and novel applications based on the isokinetic ensemble. I. Biophysical and path integral molecular dynamics , 2003 .

[34]  Jesús A. Izaguirre,et al.  Targeted Mollified Impulse: A Multiscale Stochastic Integrator for Long Molecular Dynamics Simulations , 2003, Multiscale Model. Simul..

[35]  Masaru Yoshida,et al.  Molecular theory of solvation for supramolecules and soft matter structures: application to ligand binding, ion channels, and oligomeric polyelectrolyte gelators , 2012 .

[36]  Fumio Hirata,et al.  Self-consistent description of a metal–water interface by the Kohn–Sham density functional theory and the three-dimensional reference interaction site model , 1999 .

[37]  Fumio Hirata,et al.  Solution of three‐dimensional reference interaction site model and hypernetted chain equations for simple point charge water by modified method of direct inversion in iterative subspace , 1999 .

[38]  Andriy Kovalenko,et al.  Multiscale methods for nanochemistry and biophysics in solution , 2011 .

[39]  Giovanni Ciccotti,et al.  Computer Simulations in Condensed Matter Systems: From Materials to Chemical Biology. Volume 2 (Lecture Notes in Physics) , 2006 .

[40]  Ioannis G Kevrekidis,et al.  Integrating diffusion maps with umbrella sampling: application to alanine dipeptide. , 2011, The Journal of chemical physics.

[41]  Igor Omelyan,et al.  Interpretation of atomic motion in flexible molecules: accelerating molecular dynamics simulations. , 2012, Physical review. E, Statistical, nonlinear, and soft matter physics.

[42]  David Chandler,et al.  Density functional theory of nonuniform polyatomic systems. I. General formulation , 1986 .

[43]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

[44]  B. Brooks,et al.  Langevin dynamics of peptides: The frictional dependence of isomerization rates of N‐acetylalanyl‐N′‐methylamide , 1992, Biopolymers.

[45]  Michael P Eastwood,et al.  Equipartition and the Calculation of Temperature in Biomolecular Simulations. , 2010, Journal of chemical theory and computation.

[46]  Application of a Multiple Time Step Algorithm to Biomolecular Systems , 1999 .

[47]  J. Izaguirre.  Longer Time Steps for Molecular Dynamics , 1999 .

[48]  Yaohang Li,et al.  Nonequilibrium, multiple‐timescale simulations of ligand‐receptor interactions in structured protein systems , 2003, Proteins.

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

[50]  Masuhiro Mikami,et al.  Rigid‐body dynamics in the isothermal‐isobaric ensemble: A test on the accuracy and computational efficiency , 2003, J. Comput. Chem..

[51]  M. Tuckerman,et al.  Long time molecular dynamics for enhanced conformational sampling in biomolecular systems. , 2004, Physical review letters.

[52]  Noriyuki Minezawa,et al.  Efficient implementation of three-dimensional reference interaction site model self-consistent-field method: application to solvatochromic shift calculations. , 2007, The Journal of chemical physics.

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

[54]  Burkhard Dünweg,et al.  Multiple time step integrators and momentum conservation , 1997 .

[55]  T. Schlick,et al.  Extrapolation versus impulse in multiple-timestepping schemes. II. Linear analysis and applications to Newtonian and Langevin dynamics , 1998 .

[56]  M. Tuckerman,et al.  Equilibrium Statistical Mechanics, Non-Hamiltonian Molecular Dynamics, and Novel Applications from Resonance-Free Timesteps to Adiabatic Free Energy Dynamics , 2006 .

[57]  Yuichi Harano,et al.  Theoretical study for partial molar volume of amino acids and polypeptides by the three-dimensional reference interaction site model , 2001 .

[58]  Gillian C. Lynch,et al.  Protein solvation from theory and simulation: Exact treatment of Coulomb interactions in three-dimensional theories. , 2010, The Journal of chemical physics.

[59]  M. Karplus,et al.  SIMULATIONS OF MACROMOLECULES BY MULTIPLE TIME-STEP METHODS , 1995 .

[60]  Andriy Kovalenko,et al.  Multiple time scale molecular dynamics for fluids with orientational degrees of freedom. II. Canonical and isokinetic ensembles. , 2011, The Journal of chemical physics.

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

[62]  D. LeMaster,et al.  Millisecond time scale conformational flexibility in a hyperthermophile protein at ambient temperature. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[63]  Benoît Roux,et al.  NUMERICAL SOLUTION OF THE HYPERNETTED CHAIN EQUATION FOR A SOLUTE OF ARBITRARY GEOMETRY IN THREE DIMENSIONS , 1995 .

[64]  Fumio Hirata,et al.  Potentials of mean force of simple ions in ambient aqueous solution. II. Solvation structure from the three-dimensional reference interaction site model approach, and comparison with simulations , 2000 .

[65]  David Chandler,et al.  Density functional theory of nonuniform polyatomic systems. II: Rational closures for integral equations , 1986 .

[66]  Shirley Dex,et al.  JR 旅客販売総合システム（マルス）における運用及び管理について , 1991 .

[67]  S. Melchionna.  Design of quasisymplectic propagators for Langevin dynamics. , 2007, The Journal of chemical physics.

[68]  Andriy Kovalenko,et al.  Three-dimensional Rism Theory for Molecular Liquids and Solid-Liquid Interfaces , 2004 .

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

[70]  Y. Komeiji.  Ewald summation and multiple time step methods for molecular dynamics simulation of biological molecules , 2000 .

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

[72]  Andriy Kovalenko,et al.  Hydration effects on the HET-s prion and amyloid-beta fibrillous aggregates, studied with three-dimensional molecular theory of solvation. , 2008, Biophysical journal.

[73]  F. Hirata,et al.  Three-dimensional density profiles of water in contact with a solute of arbitrary shape: a RISM approach , 1998 .

[74]  Ruhong Zhou,et al.  Molecular dynamics with multiple time scales: The selection of efficient reference system propagators , 1996 .

[75]  Andriy Kovalenko,et al.  Solution of three-dimensional reference interaction site model and hypernetted chain equations for simple point charge water by modified method of direct inversion in iterative subspace , 1999, J. Comput. Chem..

[76]  Andriy Kovalenko,et al.  An MM/3D-RISM approach for ligand binding affinities. , 2010, The journal of physical chemistry. B.

[77]  R. Skeel,et al.  Nonlinear Resonance Artifacts in Molecular Dynamics Simulations , 1998 .

[78]  E. Hairer,et al.  Simulating Hamiltonian dynamics , 2006, Math. Comput..

[79]  I. Omelyan.  A New Leapfrog Integrator of Rotational Motion. The Revised Angular-Momentum Approach , 1999, physics/9901025.

[80]  H. Fenniri,et al.  Structural water drives self-assembly of organic rosette nanotubes and holds host atoms in the channel. , 2010, Chemphyschem : a European journal of chemical physics and physical chemistry.

[81]  Longitudinal wavevector-and frequency-dependent dielectric constant of the TIP 4 P water model , 1998 .

[82]  Giovanni Ciccotti,et al.  Molecular dynamics of rigid systems in cartesian coordinates: A general formulation , 1982 .

[83]  Robert D. Skeel,et al.  Long-Time-Step Methods for Oscillatory Differential Equations , 1998, SIAM J. Sci. Comput..

[84]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[85]  Andriy Kovalenko,et al.  Association thermodynamics and conformational stability of beta-sheet amyloid beta(17-42) oligomers: effects of E22Q (Dutch) mutation and charge neutralization. , 2010, Biophysical journal.

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

[87]  T. Schlick,et al.  Resonance in the dynamics of chemical systems simulated by the implicit midpoint scheme , 1995 .

[88]  Andriy Kovalenko,et al.  Efficient treatment of solvation shells in 3D molecular theory of solvation , 2012, J. Comput. Chem..

[89]  T. Schlick,et al.  Efficient multiple-time-step integrators with distance-based force splitting for particle-mesh-Ewald molecular dynamics simulations , 2002 .

[90]  Bernard R. Brooks,et al.  New spherical‐cutoff methods for long‐range forces in macromolecular simulation , 1994, J. Comput. Chem..

[91]  Ronald M. Levy,et al.  Long-time conformational transitions of alanine dipeptide in aqueous solution: Continuous and discrete-state kinetic models , 2004 .

[92]  Yuko Okamoto,et al.  Effects of the fixed end in single-molecule imaging techniques: A replica-exchange molecular dynamics simulation study , 2005 .

[93]  C. Brooks.  Computer simulation of liquids , 1989 .

[94]  B. Pettitt,et al.  An Integral Equation Study of the Hydrophobic Interaction between Graphene Plates. , 2008, Journal of chemical theory and computation.

[95]  V. Pande,et al.  Multiplexed-replica exchange molecular dynamics method for protein folding simulation. , 2003, Biophysical journal.

[96]  Andriy Kovalenko,et al.  Overcoming the Barrier on Time Step Size in Multiscale Molecular Dynamics Simulation of Molecular Liquids. , 2012, Journal of chemical theory and computation.

[97]  Fumio Hirata,et al.  Potentials of mean force of simple ions in ambient aqueous solution. I. Three-dimensional reference interaction site model approach , 2000 .

[98]  Yuefan Deng,et al.  Error and timing analysis of multiple time-step integration methods for molecular dynamics , 2007, Comput. Phys. Commun..

[99]  Fumio Hirata,et al.  Potential of Mean Force between Two Molecular Ions in a Polar Molecular Solvent: A Study by the Three-Dimensional Reference Interaction Site Model , 1999 .

[100]  Andriy Kovalenko,et al.  Modeling solvatochromic shifts using the orbital-free embedding potential at statistically mechanically averaged solvent density. , 2010, The journal of physical chemistry. A.

[101]  E. M.,et al.  Statistical Mechanics , 2021, Manual for Theoretical Chemistry.

[102]  Siu A Chin.  Dynamical multiple-time stepping methods for overcoming resonance instabilities. , 2004, The Journal of chemical physics.

[103]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[104]  I. Fdez. Galván,et al.  Study on the conformational equilibrium of the alanine dipeptide in water solution by using the averaged solvent electrostatic potential from molecular dynamics methodology. , 2011, The Journal of chemical physics.

[105]  Edward D Harder,et al.  Efficient multiple time step method for use with Ewald and particle mesh Ewald for large biomolecular systems , 2001 .

[106]  R. Mazo,et al.  On the Statistical Mechanical Theory of Solutions , 1958 .