Efficient computation of absolute free energies of binding by computer simulations. Application to the methane dimer in water

An efficient procedure is noted for computing absolute free energies of binding for complexes in solution. Two series of computer simulations are required in which the substrate is annihilated in the solvent by itself and in the solvated complex. For illustration, the free energy of binding for two methane‐like particles at their contact separation of 4 A has been computed in TIP4P water. Though several alternatives are possible, in this case, Monte Carlo simulations were employed with statistical perturbation theory in the NPT ensemble at 25 °C and 1 atm. The results for the free energy of binding as well as for the potential of mean force are consistent with prior findings from the integral equation theory of Pratt and Chandler.

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

[2]  William L. Jorgensen,et al.  Hydration and energetics for tert-butyl chloride ion pairs in aqueous solution , 1987 .

[3]  Bruce J. Berne,et al.  A Monte Carlo simulation of the hydrophobic interaction , 1979 .

[4]  Solute-solute potential of mean force for non-polar molecules in water , 1984 .

[5]  P. A. Bash,et al.  Calculation of the relative change in binding free energy of a protein-inhibitor complex. , 1987, Science.

[6]  Mihaly Mezei,et al.  Monte Carlo computer simulation study of the hydrophobic effect. Potential of mean force for [(CH4)2]aq at 25 and 50 °C , 1982 .

[7]  Arieh Ben-Naim,et al.  Solvation thermodynamics of nonionic solutes , 1984 .

[8]  B. Berne,et al.  HYDROPHOBIC INTERACTION BETWEEN A METHANE MOLECULE AND A PARAFFIN WALL IN LIQUID WATER , 1988 .

[9]  D. Chandler,et al.  Hydrophobic solvation of nonspherical solutes , 1980 .

[10]  Peter A. Kollman,et al.  Calculation of the free energy of association of nucleic acid bases in vacuo and water solution , 1988 .

[11]  James Andrew McCammon,et al.  Ligand-receptor interactions , 1984, Comput. Chem..

[12]  U. Singh,et al.  Free energy perturbation calculations on binding and catalysis after mutating Asn 155 in subtilisin , 1987, Nature.

[13]  W. L. Jorgensen,et al.  Monte Carlo simulation of differences in free energies of hydration , 1985 .

[14]  William L. Jorgensen,et al.  Optimized intermolecular potential functions for liquid hydrocarbons , 1984 .

[15]  William L. Jorgensen,et al.  Theoretical studies of medium effects on conformational equilibria , 1983 .

[16]  J A McCammon,et al.  Theoretical calculation of relative binding affinity in host-guest systems. , 1986, Proceedings of the National Academy of Sciences of the United States of America.

[17]  D. Chandler,et al.  Effects of solute--solvent attractive forces on hydrophobic correlations , 1980 .

[18]  Robert B. Hermann,et al.  Theory of hydrophobic bonding. II. Correlation of hydrocarbon solubility in water with solvent cavity surface area , 1972 .

[19]  Kyoko Watanabe,et al.  Molecular dynamics study of the hydrophobic interaction in an aqueous solution of krypton , 1986 .

[20]  J. Andrew McCammon,et al.  Dynamics and design of enzymes and inhibitors , 1986 .

[21]  William L. Jorgensen,et al.  Monte Carlo simulations of alkanes in water: hydration numbers and the hydrophobic effect , 1985 .

[22]  D. Chandler,et al.  Theory of the hydrophobic effect , 1977 .