Free Energy Reconstruction from Logarithmic Mean-Force Dynamics Using Multiple Nonequilibrium Trajectories.

Efficient and reliable estimation of the mean force (MF), the derivatives of the free energy with respect to a set of collective variables (CVs), has been a challenging problem because free energy differences are often computed by integrating the MF. Among various methods for computing free energy differences, logarithmic mean-force dynamics (LogMFD) [ Morishita et al., Phys. Rev. E 2012 , 85 , 066702 ] invokes the conservation law in classical mechanics to integrate the MF, which allows us to estimate the free energy profile along the CVs on-the-fly. Here, we present a method called parallel dynamics, which improves the estimation of the MF by employing multiple replicas of the system and is straightforwardly incorporated in LogMFD or a related method. In the parallel dynamics, the MF is evaluated by a nonequilibrium path-ensemble using the multiple replicas based on the Crooks-Jarzynski nonequilibrium work relation. Thanks to the Crooks relation, realizing full-equilibrium states is no longer mandatory for estimating the MF. Additionally, sampling in the hidden subspace orthogonal to the CV space is highly improved with appropriate weights for each metastable state (if any), which is hardly achievable by typical free energy computational methods. We illustrate how to implement parallel dynamics by combining it with LogMFD, which we call logarithmic parallel dynamics (LogPD). Biosystems of alanine dipeptide and adenylate kinase in explicit water are employed as benchmark systems to which LogPD is applied to demonstrate the effect of multiple replicas on the accuracy and efficiency in estimating the free energy profiles using parallel dynamics.

[1]  R. Chelli,et al.  Path-breaking schemes for nonequilibrium free energy calculations. , 2013, The Journal of chemical physics.

[2]  G. Schulz,et al.  Adenylate kinase motions during catalysis: an energetic counterweight balancing substrate binding. , 1996, Structure.

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

[4]  Masuhiro Mikami,et al.  On‐the‐fly reconstruction of free‐energy profiles using logarithmic mean‐force dynamics , 2013, J. Comput. Chem..

[5]  Alessandro Laio,et al.  Efficient exploration of reactive potential energy surfaces using Car-Parrinello molecular dynamics. , 2003, Physical review letters.

[6]  Christophe Chipot,et al.  Good practices in free-energy calculations. , 2010, The journal of physical chemistry. B.

[7]  M. Tuckerman,et al.  On the use of the adiabatic molecular dynamics technique in the calculation of free energy profiles , 2002 .

[8]  Masuhiro Mikami,et al.  Efficient free energy calculation of water across lipid membranes , 2008, J. Comput. Chem..

[9]  Yasushige Yonezawa,et al.  A method for predicting protein conformational pathways by using molecular dynamics simulations guided by difference distance matrices , 2016, J. Comput. Chem..

[10]  A. Laio,et al.  Metadynamics: a method to simulate rare events and reconstruct the free energy in biophysics, chemistry and material science , 2008 .

[11]  R. Zwanzig.  High‐Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases , 1954 .

[12]  H. Okumura,et al.  Free-energy calculation via mean-force dynamics using a logarithmic energy landscape. , 2012, Physical review. E, Statistical, nonlinear, and soft matter physics.

[13]  G. Ciccotti,et al.  Temperature‐accelerated molecular dynamics gives insights into globular conformations sampled in the free state of the AC catalytic domain , 2014, Proteins.

[14]  Y. Fukunishi,et al.  The Filling Potential Method: A Method for Estimating the Free Energy Surface for Protein−Ligand Docking , 2003 .

[15]  M. Klein,et al.  Constant pressure molecular dynamics algorithms , 1994 .

[16]  A. Laio,et al.  Efficient reconstruction of complex free energy landscapes by multiple walkers metadynamics. , 2006, The journal of physical chemistry. B.

[17]  A. Laio,et al.  Flexible docking in solution using metadynamics. , 2005, Journal of the American Chemical Society.

[18]  J. Kirkwood.  Statistical Mechanics of Fluid Mixtures , 1935 .

[19]  Rigoberto Hernandez,et al.  Adaptive Steered Molecular Dynamics of the Long-Distance Unfolding of Neuropeptide Y. , 2010, Journal of chemical theory and computation.

[20]  C. Mann,et al.  A Practical Treatise on Diseases of the Skin , 1889, Atlanta Medical and Surgical Journal (1884).

[21]  Bharat V. Adkar,et al.  Dynamic coupling between the LID and NMP domain motions in the catalytic conversion of ATP and AMP to ADP by adenylate kinase. , 2010, The Journal of chemical physics.

[22]  Makoto Nakamura,et al.  An ab initio approach to free-energy reconstruction using logarithmic mean force dynamics. , 2014, The Journal of chemical physics.

[23]  Chris Oostenbrink,et al.  Free‐energy differences between states with different conformational ensembles , 2013, J. Comput. Chem..

[24]  S. Nosé.  A molecular dynamics method for simulations in the canonical ensemble , 1984 .

[25]  Cristian Micheletti,et al.  Small- and large-scale conformational changes of adenylate kinase: a molecular dynamics study of the subdomain motion and mechanics. , 2008, Biophysical journal.

[26]  Tetsuya Morishita,et al.  A New Surface and Structure for Silicene: Polygonal Silicene Formation on the Al(111) Surface , 2013 .

[27]  C. W. Gear,et al.  'Coarse' integration/bifurcation analysis via microscopic simulators: Micro-Galerkin methods , 2002 .

[28]  Hiroshi Fujisaki,et al.  Extended Phase-Space Methods for Enhanced Sampling in Molecular Simulations: A Review , 2015, Front. Bioeng. Biotechnol..

[29]  E. Meirovitch,et al.  Backbone dynamics of escherichia coli adenylate kinase at the extreme stages of the catalytic cycle studied by (15)N NMR relaxation. , 2000, Biochemistry.

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

[31]  Rafael C. Bernardi,et al.  Enhanced sampling techniques in molecular dynamics simulations of biological systems. , 2015, Biochimica et biophysica acta.

[32]  Takeshi Yamamoto.  Comment on "comment on 'simple reversible molecular dynamics algorithms for Nosé-Hoover chain dynamics' " [J. Chem. Phys. 110, 3623 (1999)]. , 2006, The Journal of chemical physics.

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

[34]  C. Jarzynski.  Nonequilibrium Equality for Free Energy Differences , 1996, cond-mat/9610209.

[35]  G. Ciccotti,et al.  Constrained reaction coordinate dynamics for the simulation of rare events , 1989 .

[36]  Mark E Tuckerman,et al.  Temperature-accelerated method for exploring polymorphism in molecular crystals based on free energy. , 2011, Physical review letters.

[37]  E. Vanden-Eijnden,et al.  A temperature accelerated method for sampling free energy and determining reaction pathways in rare events simulations , 2006 .

[38]  Evans,et al.  Steady states, invariant measures, and response theory. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[39]  M. Tuckerman,et al.  Generalized Gaussian moment thermostatting: A new continuous dynamical approach to the canonical ensemble , 2000 .

[40]  Eric Vanden-Eijnden,et al.  Some recent techniques for free energy calculations , 2009, J. Comput. Chem..

[41]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[42]  Shuichi Nosé,et al.  Constant Temperature Molecular Dynamics Methods , 1991 .

[43]  William G. Hoover,et al.  Nonequilibrium molecular dynamics via Gauss's principle of least constraint , 1983 .

[44]  Jackson,et al.  Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlation. , 1992, Physical review. B, Condensed matter.

[45]  M. Tuckerman,et al.  Free Energy Reconstruction from Metadynamics or Adiabatic Free Energy Dynamics Simulations. , 2014, Journal of chemical theory and computation.

[46]  M. Tuckerman,et al.  Efficient and direct generation of multidimensional free energy surfaces via adiabatic dynamics without coordinate transformations. , 2008, The journal of physical chemistry. B.

[47]  Tetsuya Morishita,et al.  From Nosé–Hoover chain to Nosé–Hoover network: design of non-Hamiltonian equations of motion for molecular-dynamics with multiple thermostats , 2010 .

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

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

[50]  I. Kevrekidis,et al.  Coarse molecular dynamics of a peptide fragment: Free energy, kinetics, and long-time dynamics computations , 2002, physics/0212108.

[51]  Hiroshi Fujisaki,et al.  Minimum Free Energy Path of Ligand-Induced Transition in Adenylate Kinase , 2012, PLoS Comput. Biol..

[52]  M. Klein,et al.  Nosé-Hoover chains : the canonical ensemble via continuous dynamics , 1992 .

[53]  C. Jarzynski,et al.  Escorted free energy simulations: improving convergence by reducing dissipation. , 2008, Physical review letters.

[54]  T. Keyes,et al.  Multiple branched adaptive steered molecular dynamics. , 2014, The Journal of chemical physics.

[55]  Alessandro Laio,et al.  A recipe for the computation of the free energy barrier and the lowest free energy path of concerted reactions. , 2005, The journal of physical chemistry. B.

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

[57]  G. Crooks.  Path-ensemble averages in systems driven far from equilibrium , 1999, cond-mat/9908420.