Calculating zeros: Non-equilibrium free energy calculations

Free energy calculations on three model processes with theoretically known free energy changes have been performed using short simulation times. A comparison between equilibrium (thermodynamic integration) and non-equilibrium (fast growth) methods has been made in order to assess the accuracy and precision of these methods. The three processes have been chosen to represent processes often observed in biomolecular free energy calculations. They involve a redistribution of charges, the creation and annihilation of neutral particles and conformational changes. At very short overall simulation times, the thermodynamic integration approach using discrete steps is most accurate. More importantly, reasonable accuracy can be obtained using this method which seems independent of the overall simulation time. In cases where slow conformational changes play a role, fast growth simulations might have an advantage over discrete thermodynamic integration where sufficient sampling needs to be obtained at every λ-point, but only if the initial conformations do properly represent an equilibrium ensemble. From these three test cases practical lessons can be learned that will be applicable to biomolecular free energy calculations.

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

[2]  Jan Hermans,et al.  Reversibility of Free Energy Simulations∶ Slow Growth May Have a Unique Advantage. (With a Note on Use of Ewald Summation) , 2002 .

[3]  D. Zuckerman,et al.  Efficient use of nonequilibrium measurement to estimate free energy differences for molecular systems , 2004, Journal of computational chemistry.

[4]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[5]  Mark A. Miller,et al.  Why is it so difficult to simulate entropies, free energies, and their differences? , 2001, Accounts of chemical research.

[6]  Peter A. Kollman,et al.  The lag between the Hamiltonian and the system configuration in free energy perturbation calculations , 1989 .

[7]  C. Jarzynski,et al.  A “fast growth” method of computing free energy differences , 2001 .

[8]  Andrew E. Torda,et al.  The GROMOS biomolecular simulation program package , 1999 .

[9]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[10]  C. Jarzynski Equilibrium free-energy differences from nonequilibrium measurements: A master-equation approach , 1997, cond-mat/9707325.

[11]  K. Schulten,et al.  Calculating potentials of mean force from steered molecular dynamics simulations. , 2004, The Journal of chemical physics.

[12]  W. V. van Gunsteren,et al.  Estimating entropies from molecular dynamics simulations. , 2004, The Journal of chemical physics.

[13]  Wilfred F. van Gunsteren,et al.  Predictions of free energy differences from a single simulation of the initial state , 1994 .

[14]  M Mezei,et al.  Free Energy Simulations a , 1986, Annals of the New York Academy of Sciences.

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

[16]  E. Cohen,et al.  A note on the Jarzynski equality , 2004 .

[17]  Alan E. Mark,et al.  Estimating the Relative Free Energy of Different Molecular States with Respect to a Single Reference State , 1996 .

[18]  Eric Darve,et al.  Assessing the efficiency of free energy calculation methods. , 2004, The Journal of chemical physics.

[19]  G. Hummer Fast-growth thermodynamic integration: Error and efficiency analysis , 2001 .

[20]  Tim N. Heinz,et al.  Comparison of four methods to compute the dielectric permittivity of liquids from molecular dynamics simulations , 2001 .

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

[22]  P. Schleyer Encyclopedia of computational chemistry , 1998 .

[23]  F. Ritort,et al.  Bias and error in estimates of equilibrium free-energy differences from nonequilibrium measurements , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[24]  Christophe Chipot,et al.  Free Energy Calculations. The Long and Winding Gilded Road , 2002 .

[25]  G. Hummer,et al.  Free energy reconstruction from nonequilibrium single-molecule pulling experiments , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[26]  U. Seifert,et al.  Probing molecular free energy landscapes by periodic loading. , 2004, Physical review letters.

[27]  Jed W. Pitera,et al.  A Comparison of Non-Bonded Scaling Approaches for Free Energy Calculations , 2002 .

[28]  A. Mark,et al.  Avoiding singularities and numerical instabilities in free energy calculations based on molecular simulations , 1994 .

[29]  J. Aqvist,et al.  A new method for predicting binding affinity in computer-aided drug design. , 1994, Protein engineering.

[30]  Robert H. Wood,et al.  Systematic errors in free energy perturbation calculations due to a finite sample of configuration space: sample-size hysteresis , 1991 .

[31]  Irwin Oppenheim,et al.  Statistical Mechanical Theory of Transport Processes. VII. The Coefficient of Thermal Conductivity of Monatomic Liquids , 1954 .

[32]  Wilfred F. van Gunsteren,et al.  An improved GROMOS96 force field for aliphatic hydrocarbons in the condensed phase , 2001, J. Comput. Chem..

[33]  Martin Almlöf,et al.  Free energy calculations and ligand binding. , 2003, Advances in protein chemistry.

[34]  I. Tinoco,et al.  Equilibrium Information from Nonequilibrium Measurements in an Experimental Test of Jarzynski's Equality , 2002, Science.

[35]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[36]  Wilfred F. van Gunsteren,et al.  A generalized reaction field method for molecular dynamics simulations , 1995 .

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

[38]  G. Hummer Fast-growth Thermodynamic Integration∶ Results for Sodium Ion Hydration , 2002 .

[39]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .