Effect of the integration method on the accuracy and computational efficiency of free energy calculations using thermodynamic integration

Although calculations of free energy using molecular dynamics simulations have gained significant importance in the chemical and biochemical fields, they still remain quite computationally intensive. Furthermore, when using thermodynamic integration, numerical evaluation of the integral of the Hamiltonian with respect to the coupling parameter may introduce unwanted errors in the free energy. In this paper, we compare the performance of two numerical integration techniques—the trapezoidal and Simpson’s rules—and propose a new method, based on the analytic integration of physically based fitting functions that are able to accurately describe the behavior of the data. We develop and test our methodology by performing detailed studies on two prototype systems, hydrated methane and hydrated methanol, and treat Lennard-Jones and electrostatic contributions separately. We conclude that the widely used trapezoidal rule may introduce systematic errors in the calculation, but these errors are reduced if Simpson’s ...

[1]  William C. Swope,et al.  A molecular dynamics method for calculating the solubility of gases in liquids and the hydrophobic hydration of inert-gas atoms in aqueous solution , 1984 .

[2]  Sture Nordholm,et al.  In silico prediction of drug solubility: 1. Free energy of hydration. , 2007, The journal of physical chemistry. B.

[3]  J. Guthrie,et al.  A blind challenge for computational solvation free energies: introduction and overview. , 2009, The journal of physical chemistry. B.

[4]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997, J. Comput. Chem..

[5]  Conrad Shyu,et al.  Reducing the bias and uncertainty of free energy estimates by using regression to fit thermodynamic integration data , 2009, J. Comput. Chem..

[6]  Gerhard Hummer,et al.  Free Energy of Ionic Hydration , 1996 .

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

[8]  David L Mobley,et al.  Small molecule hydration free energies in explicit solvent: An extensive test of fixed-charge atomistic simulations. , 2009, Journal of chemical theory and computation.

[9]  Howard Reiss,et al.  Statistical Mechanics of Rigid Spheres , 1959 .

[10]  Herman J. C. Berendsen,et al.  ALGORITHMS FOR BROWNIAN DYNAMICS , 1982 .

[11]  Vincenzo Mollica,et al.  Group contributions to the thermodynamic properties of non-ionic organic solutes in dilute aqueous solution , 1981 .

[12]  A. Pohorille,et al.  Free energy calculations : theory and applications in chemistry and biology , 2007 .

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

[14]  H. Berendsen,et al.  A LEAP-FROG ALGORITHM FOR STOCHASTIC DYNAMICS , 1988 .

[15]  M. Gilson,et al.  Calculation of protein-ligand binding affinities. , 2007, Annual review of biophysics and biomolecular structure.

[16]  Ioannis G. Economou,et al.  Molecular simulation of absolute hydration Gibbs energies of polar compounds , 2010 .

[17]  Ioannis G. Economou,et al.  Molecular simulation of the hydration Gibbs energy of barbiturates , 2010 .

[18]  Arieh Warshel,et al.  A local reaction field method for fast evaluation of long‐range electrostatic interactions in molecular simulations , 1992 .

[19]  Shunzhou Wan,et al.  Calculation of the aqueous solvation energy and entropy, as well as free energy, of simple polar solutes. , 2004, The Journal of chemical physics.

[20]  H. G. Petersen,et al.  Error estimates on averages of correlated data , 1989 .

[21]  Frank H. Stillinger,et al.  Structure in aqueous solutions of nonpolar solutes from the standpoint of scaled-particle theory , 1973 .

[22]  Michael R. Shirts,et al.  Solvation free energies of amino acid side chain analogs for common molecular mechanics water models. , 2005, The Journal of chemical physics.

[23]  T. Straatsma,et al.  Free energy of ionic hydration: Analysis of a thermodynamic integration technique to evaluate free energy differences by molecular dynamics simulations , 1988 .

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

[25]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[26]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[27]  W. L. Jorgensen.  The Many Roles of Computation in Drug Discovery , 2004, Science.

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

[29]  Michael R. Shirts,et al.  Extremely precise free energy calculations of amino acid side chain analogs: Comparison of common molecular mechanics force fields for proteins , 2003 .

[30]  Jorge Nocedal,et al.  On the limited memory BFGS method for large scale optimization , 1989, Math. Program..

[31]  David L Mobley,et al.  Comparison of charge models for fixed-charge force fields: small-molecule hydration free energies in explicit solvent. , 2007, The journal of physical chemistry. B.

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

[33]  R. Pierotti,et al.  A scaled particle theory of aqueous and nonaqueous solutions , 1976 .

[34]  A D J Haymet,et al.  Free energy of solvation of simple ions: molecular-dynamics study of solvation of Cl- and Na+ in the ice/water interface. , 2005, The Journal of chemical physics.

[35]  Magdalena Bacilieri,et al.  Prediction of the aqueous solvation free energy of organic compounds by using autocorrelation of molecular electrostatic potential surface properties combined with response surface analysis. , 2008, Bioorganic & medicinal chemistry.

[36]  Chris Oostenbrink,et al.  A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS force‐field parameter sets 53A5 and 53A6 , 2004, J. Comput. Chem..