Estimation of Absolute Free Energies of Hydration Using Continuum Methods: Accuracy of Partial Charge Models and Optimization of Nonpolar Contributions

Absolute free energies of hydration (ΔGhyd) for more than 500 neutral and charged compounds have been computed, using Poisson-Boltzmann (PB) and Generalized Born (GB) continuum methods plus a solvent-accessible surface area (SA) term, to evaluate the accuracy of eight simple point-charge models used in molecular modeling. The goal is to develop improved procedures and protocols for protein-ligand binding calculations and virtual screening (docking). The best overall PBSA and GBSA results, in comparison with experimental ΔGhyd values for small molecules, were obtained using MSK, RESP, or ChelpG charges obtained from ab initio calculations using 6-31G* wave functions. Correlations using semiempirical (AM1BCC, AM1CM2, and PM3CM2) or empirical (Gasteiger-Marsili and MMFF94) methods yielded mixed results, particularly for charged compounds. For neutral compounds, the AM1BCC method yielded the best agreement with experimental results. In all cases, the PBSA and GBSA results are highly correlated (overall r(2) = 0.94), which highlights the fact that various partial charge models influence the final results much more than which continuum method is used to compute hydration free energies. Overall improved agreement with experimental results was demonstrated using atom-based constants in place of a single surface area term. Sets of optimized SA constants, suitable for use with a given charge model, were derived by fitting to the difference in experimental free energies and polar continuum results. The use of optimized atom-based SA constants for the computation of ΔGhyd can fine-tune already reasonable agreement with experimental results, ameliorate gross deficiencies in any particular charge model, account for nonoptimal radii, or correct for systematic errors.

[1]  Ronald M. Levy,et al.  AGBNP: An analytic implicit solvent model suitable for molecular dynamics simulations and high‐resolution modeling , 2004, J. Comput. Chem..

[2]  Irwin D Kuntz,et al.  Investigation of neuraminidase-substrate recognition using molecular dynamics and free energy calculations. , 2003, Journal of medicinal chemistry.

[3]  Christopher I. Bayly,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: II. Parameterization and validation , 2002, J. Comput. Chem..

[4]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[5]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[6]  Gregory D. Hawkins,et al.  Pairwise solute descreening of solute charges from a dielectric medium , 1995 .

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

[8]  Claudio N. Cavasotto,et al.  Accurate Transferable Model for Water, n-Octanol, and n-Hexadecane Solvation Free Energies , 2002 .

[9]  F. Lombardo,et al.  Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings , 1997 .

[10]  P A Kollman,et al.  Calculation and prediction of binding free energies for the matrix metalloproteinases. , 2000, Journal of medicinal chemistry.

[11]  Barry Honig,et al.  Extending the Applicability of the Nonlinear Poisson−Boltzmann Equation: Multiple Dielectric Constants and Multivalent Ions† , 2001 .

[12]  Paul R. Gerber,et al.  Charge distribution from a simple molecular orbital type calculation and non-bonding interaction terms in the force field MAB , 1998, J. Comput. Aided Mol. Des..

[13]  Irwin D Kuntz,et al.  A molecular basis for the selectivity of thiadiazole urea inhibitors with stromelysin-1 and gelatinase-A from generalized born molecular dynamics simulations. , 2004, Journal of medicinal chemistry.

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

[15]  J. Åqvist,et al.  Ion-water interaction potentials derived from free energy perturbation simulations , 1990 .

[16]  Gregory D. Hawkins,et al.  Parametrized Models of Aqueous Free Energies of Solvation Based on Pairwise Descreening of Solute Atomic Charges from a Dielectric Medium , 1996 .

[17]  Michael H. Abraham,et al.  Thermodynamics of solute transfer from water to hexadecane , 1990 .

[18]  Akihiko Konagaya,et al.  Molecular dynamics, free energy, and SPR analyses of the interactions between the SH2 domain of Grb2 and ErbB phosphotyrosyl peptides. , 2003, Biochemistry.

[19]  P. Kollman,et al.  Binding of a diverse set of ligands to avidin and streptavidin: an accurate quantitative prediction of their relative affinities by a combination of molecular mechanics and continuum solvent models. , 2000, Journal of medicinal chemistry.

[20]  Yizhak Marcus,et al.  Thermodynamics of solvation of ions. Part 5.—Gibbs free energy of hydration at 298.15 K , 1991 .

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

[22]  P. Kollman,et al.  Application of RESP charges to calculate conformational energies, hydrogen bond energies, and free energies of solvation , 1993 .

[23]  C. Breneman,et al.  Determining atom‐centered monopoles from molecular electrostatic potentials. The need for high sampling density in formamide conformational analysis , 1990 .

[24]  Donald G. Truhlar,et al.  New Class IV Charge Model for Extracting Accurate Partial Charges from Wave Functions , 1998 .

[25]  W. L. Jorgensen.  Free energy calculations: a breakthrough for modeling organic chemistry in solution , 1989 .

[26]  Yizhak Marcus,et al.  A simple empirical model describing the thermodynamics of hydration of ions of widely varying charges, sizes, and shapes , 1994 .

[27]  Carmay Lim,et al.  Theory of Ionic Hydration: Insights from Molecular Dynamics Simulations and Experiment , 1999 .

[28]  A. D. McLachlan,et al.  Solvation energy in protein folding and binding , 1986, Nature.

[29]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[30]  H. Scheraga,et al.  Accessible surface areas as a measure of the thermodynamic parameters of hydration of peptides. , 1987, Proceedings of the National Academy of Sciences of the United States of America.

[31]  P. Kollman,et al.  Continuum Solvent Studies of the Stability of DNA, RNA, and Phosphoramidate−DNA Helices , 1998 .

[32]  B. Honig,et al.  New Model for Calculation of Solvation Free Energies: Correction of Self-Consistent Reaction Field Continuum Dielectric Theory for Short-Range Hydrogen-Bonding Effects , 1996 .

[33]  A. Bondi.  van der Waals Volumes and Radii , 1964 .

[34]  P. Kollman,et al.  Use of MM-PBSA in reproducing the binding free energies to HIV-1 RT of TIBO derivatives and predicting the binding mode to HIV-1 RT of efavirenz by docking and MM-PBSA. , 2001, Journal of the American Chemical Society.

[35]  David A. Case,et al.  Effective Born radii in the generalized Born approximation: The importance of being perfect , 2002, J. Comput. Chem..

[36]  Piotr Cieplak,et al.  Molecular dynamics and free energy analyses of cathepsin D-inhibitor interactions: insight into structure-based ligand design. , 2002, Journal of medicinal chemistry.

[37]  Michael H. Abraham,et al.  Hydrogen bonding. Part 34. The factors that influence the solubility of gases and vapours in water at 298 K, and a new method for its determination , 1994 .

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

[39]  D. Case,et al.  Molecular Dynamics Simulations of Nucleic Acids with a Generalized Born Solvation Model , 2000 .

[40]  William L. Jorgensen,et al.  Free Energies of Hydration from a Generalized Born Model and an All-Atom Force Field , 2004 .

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

[42]  Charles L. Brooks,et al.  Performance comparison of generalized born and Poisson methods in the calculation of electrostatic solvation energies for protein structures , 2004, J. Comput. Chem..

[43]  D. Beglov,et al.  Atomic Radii for Continuum Electrostatics Calculations Based on Molecular Dynamics Free Energy Simulations , 1997 .

[44]  P A Kollman,et al.  An analysis of the interactions between the Sem-5 SH3 domain and its ligands using molecular dynamics, free energy calculations, and sequence analysis. , 2001, Journal of the American Chemical Society.

[45]  J. Gasteiger,et al.  ITERATIVE PARTIAL EQUALIZATION OF ORBITAL ELECTRONEGATIVITY – A RAPID ACCESS TO ATOMIC CHARGES , 1980 .

[46]  Donald G. Truhlar,et al.  MODEL FOR AQUEOUS SOLVATION BASED ON CLASS IV ATOMIC CHARGES AND FIRST SOLVATION SHELL EFFECTS , 1996 .

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

[48]  William L. Jorgensen,et al.  Accuracy of free energies of hydration using CM1 and CM3 atomic charges , 2004, J. Comput. Chem..

[49]  A. Leo,et al.  Partition coefficients and their uses , 1971 .

[50]  B. Roux,et al.  Atomic Radii for Continuum Electrostatics Calculations on Nucleic Acids , 2002 .

[51]  Araz Jakalian,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: I. Method , 2000 .

[52]  Kathryn N. Rankin,et al.  On the transferability of hydration‐parametrized continuum electrostatics models to solvated binding calculations , 2003, J. Comput. Chem..

[53]  T. Halgren.  Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94 , 1996, J. Comput. Chem..

[54]  Emil Alexov,et al.  Rapid grid‐based construction of the molecular surface and the use of induced surface charge to calculate reaction field energies: Applications to the molecular systems and geometric objects , 2002, J. Comput. Chem..

[55]  D. Case,et al.  Theory and applications of the generalized born solvation model in macromolecular simulations , 2000, Biopolymers.

[56]  K. Sharp,et al.  Accurate Calculation of Hydration Free Energies Using Macroscopic Solvent Models , 1994 .

[57]  P. Kollman,et al.  Combined molecular mechanical and continuum solvent approach (MM-PBSA/GBSA) to predict ligand binding , 2000 .

[58]  Thomas A. Halgren.  Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94 , 1996, J. Comput. Chem..

[59]  Gregory D. Hawkins,et al.  Extension of the platform of applicability of the SM5.42R universal solvation model , 1999 .