Rational determination of charge distributions for free energy calculations

Point charges derived from RHF/6‐31G* electrostatic potentials are attractive because they tend to exaggerate the polarity of solvated molecules, thereby compensating in an average fashion missing induction effects. In the context of free energy calculations, wherein the molecule is transferred from a polar environment to a nonpolar one, we propose a more rational approach based on a self‐consistent reaction field computation at a higher level of theory, supplemented by an estimation of the corresponding distortion energy to account for the change of polarity of the surroundings. Application of this method to the test cases acetamide, acetic acid, methyl acetate and phenol, using multinanosecond molecular dynamics/“umbrella sampling” simulations, yields consistent hydration free energies in reasonably good agreement with experiment. © 2003 Wiley Periodicals, Inc. J Comput Chem 24: 409–415, 2003

[1]  Tony F. Heinz,et al.  Studies of Liquid Surfaces by Second Harmonic Generation , 1986 .

[2]  Jean-Louis Rivail,et al.  A quantum chemical approach to dielectric solvent effects in molecular liquids , 1976 .

[3]  A. Pohorille,et al.  Molecular dynamics of phenol at the liquid-vapor interface of water. , 1991, The Journal of chemical physics.

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

[5]  J. Board,et al.  Ewald summation techniques in perspective: a survey , 1996 .

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

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

[8]  Donald E. Williams,et al.  Representation of the molecular electrostatic potential by a net atomic charge model , 1981 .

[9]  Christophe Chipot,et al.  Rational determination of transfer free energies of small drugs across the water-oil interface. , 2002, Journal of medicinal chemistry.

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

[11]  C. Brooks.  Computer simulation of liquids , 1989 .

[12]  János G. Ángyán,et al.  Common theoretical framework for quantum chemical solvent effect theories , 1992 .

[13]  Modesto Orozco,et al.  New strategies to incorporate the solvent polarization in self‐consistent reaction field and free‐energy perturbation simulations , 1995 .

[14]  Heather A. Carlson,et al.  Accuracy of free energies of hydration for organic molecules from 6‐31g*‐derived partial charges , 1993, J. Comput. Chem..

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

[16]  F. J. Luque,et al.  Theoretical Methods for the Description of the Solvent Effect in Biomolecular Systems. , 2000, Chemical reviews.

[17]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[18]  P. Kollman,et al.  Atomic charges derived from semiempirical methods , 1990 .

[19]  K. Merz,et al.  Analysis of a large data base of electrostatic potential derived atomic charges , 1992 .

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

[21]  F. J. Luque,et al.  Generalized linear response approximation in discrete methods , 1997 .

[22]  A. Pohorille,et al.  Viewpoint 9--molecular structure of aqueous interfaces. , 1993, Journal of molecular structure.

[23]  William L. Jorgensen,et al.  Monte Carlo simulations of the hydration of substituted benzenes with OPLS potential functions , 1993, J. Comput. Chem..

[24]  J. Hine,et al.  Structural effects on rates and equilibriums. XIX. Intrinsic hydrophilic character of organic compounds. Correlations in terms of structural contributions , 1975 .

[25]  Zi-Min Lu,et al.  Phase space structure of triatomic molecules , 1997 .

[26]  Richard Wolfenden,et al.  Cooperativity and anticooperativity in solvation by water: imidazoles, quinones, nitrophenols, nitrophenolate, and nitrothiophenolate ions , 1987 .

[27]  János G. Ángyán.  Comment on “New strategies to incorporate the solvent polarization in self-consistent reaction field and free-energy perturbation simulations” [J. Chem. Phys. 103, 10183 (1995)] , 1997 .

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

[29]  L. Onsager.  Electric Moments of Molecules in Liquids , 1936 .

[30]  Jean-Louis Rivail,et al.  LIQUID-STATE QUANTUM CHEMISTRY : AN IMPROVED CAVITY MODEL , 1994 .

[31]  Modesto Orozco,et al.  Semiclassical-Continuum Approach to the Electrostatic Free Energy of Solvation , 1997 .

[32]  Daniel Rinaldi,et al.  Lamé's functions and ellipsoidal harmonics for use in chemical physics , 1982, Comput. Chem..

[33]  D. N. Card,et al.  Monte Carlo Estimation of the Free Energy by Multistage Sampling , 1972 .

[34]  G. Torrie,et al.  Monte Carlo study of a phase‐separating liquid mixture by umbrella sampling , 1977 .

[35]  C. Cramer,et al.  Implicit Solvation Models: Equilibria, Structure, Spectra, and Dynamics. , 1999, Chemical reviews.

[36]  G. Torrie,et al.  Monte Carlo free energy estimates using non-Boltzmann sampling: Application to the sub-critical Lennard-Jones fluid , 1974 .

[37]  Jean-Louis Rivail,et al.  Reaction field factors for a multipole distribution in a cavity surrounded by a continuum , 1993 .

[38]  P. Kollman.  Advances and Continuing Challenges in Achieving Realistic and Predictive Simulations of the Properties of Organic and Biological Molecules , 1996 .