Estimating relative free energies from a single ensemble: Hydration free energies

The ability to determine the free energy of solvation for a number of small organic molecules with varying sizes and properties from the coordinate trajectory of a single simulation of a given reference state was investigated. The relative free energies were estimated from a single step perturbation using the perturbation formula. The reference state consisted of a cavity surrounded by solvent. To enhance sampling a soft‐core interaction was used for the cavity. The effect of the size of the cavity, the effective core height, and the length of simulation on the ability to reproduce results obtained from thermodynamic integration calculations was considered. The results using a single step perturbation from an appropriately chosen initial state were comparable to results from thermodynamic integration calculations for a wide range of compounds. Using a large number of compounds the computational efficiency was potentially increased by 2–3 orders of magnitude over traditional free energy approaches. Factors determining the efficiency of the approach are discussed. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1604–1617, 1999

[1]  Robert H. Wood,et al.  Estimation of errors in free energy calculations due to the lag between the hamiltonian and the system configuration , 1991 .

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

[3]  Peter A. Kollman,et al.  Free energy calculation methods: A theoretical and empirical comparison of numerical errors and a new method qualitative estimates of free energy changes , 1997, J. Comput. Chem..

[4]  B. Widom,et al.  Some Topics in the Theory of Fluids , 1963 .

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

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

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

[8]  D. Matyushov,et al.  A thermodynamic analysis of solvation in dipolar liquids , 1996 .

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

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

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

[12]  Ronald M. Levy,et al.  Gaussian fluctuation formula for electrostatic free‐energy changes in solution , 1991 .

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

[14]  Gerhard Hummer,et al.  Calculation of free‐energy differences from computer simulations of initial and final states , 1996 .

[15]  I. R. Mcdonald,et al.  Theory of simple liquids , 1998 .

[16]  P. M. Dean,et al.  New Perspectives in Drug Design , 1995 .

[17]  T. Hansson,et al.  On the Validity of Electrostatic Linear Response in Polar Solvents , 1996 .

[18]  W. Hatfield,et al.  An EPR study of Cr3+:K2SbF5 , 1973 .

[19]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[20]  W. L. Jorgensen,et al.  AN EXTENDED LINEAR RESPONSE METHOD FOR DETERMINING FREE ENERGIES OF HYDRATION , 1995 .

[21]  G. King,et al.  Calculation of electrostatic free energy differences with a time-saving approximate method , 1993 .

[22]  Wilfred F. van Gunsteren,et al.  Computer Simulation of Biomolecular Systems: Theoretical and Experimental Applications , 1989 .

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

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