I-SOLV: a new surface-based empirical model for computing solvation free energies.

We have developed a new empirical model, I-SOLV, for computing solvation free energies of neutral organic molecules. It computes the solvation free energy of a solute molecule by summing up the contributions from its component atoms. The contribution from a certain atom is determined by the solvent-accessible surface area as well as the surface tension of this atom. A total of 49 atom types are implemented in our model for classifying C, N, O, S, P, F, Cl, Br and I in common organic molecules. Their surface tensions are parameterized by using a data set of 532 neutral organic molecules with experimentally measured solvation free energies. A head-to-head comparison of our model with several other solvation models was performed on a test set of 82 molecules. Our model outperformed other solvation models, including widely used PB/SA and GB/SA models, with a mean unsigned error as low as 0.39 kcal/mol. Our study has demonstrated again that well-developed empirical solvation models are not necessarily less accurate than more sophisticated theoretical models. Empirical models may serve as appealing alternatives due to their simplicity and accuracy.

[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]  A. Bondi.  van der Waals Volumes and Radii , 1964 .

[3]  Ronald M. Levy,et al.  The SGB/NP hydration free energy model based on the surface generalized born solvent reaction field and novel nonpolar hydration free energy estimators , 2002, J. Comput. Chem..

[4]  L. Lai,et al.  Calculating partition coefficient by atom-additive method , 2000 .

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

[6]  M. Sanner,et al.  Reduced surface: an efficient way to compute molecular surfaces. , 1996, Biopolymers.

[7]  Luhua Lai,et al.  A New Atom-Additive Method for Calculating Partition Coefficients , 1997, J. Chem. Inf. Comput. Sci..

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

[9]  Donald G Truhlar,et al.  SM6:  A Density Functional Theory Continuum Solvation Model for Calculating Aqueous Solvation Free Energies of Neutrals, Ions, and Solute-Water Clusters. , 2005, Journal of chemical theory and computation.

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

[11]  J. Tomasi,et al.  Quantum mechanical continuum solvation models. , 2005, Chemical reviews.

[12]  Wei Zhang,et al.  An extended aqueous solvation model based on atom-weighted solvent accessible surface areas: SAWSA v2.0 model , 2005, Journal of molecular modeling.

[13]  Ruben Abagyan,et al.  Accurate Transferable Model for Water, n-Octanol, and n-Hexadecane Solvation Free Energies , 2002 .

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

[15]  D. Case,et al.  Generalized born models of macromolecular solvation effects. , 2000, Annual review of physical chemistry.

[16]  Sathesh Bhat,et al.  Molecular surface generation using a variable‐radius solvent probe , 2005, Proteins.

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

[18]  Donald G. Truhlar,et al.  omnisol: Fast Prediction of Free Energies of Solvation and Partition Coefficients , 1998 .

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

[20]  Irwin D Kuntz,et al.  Estimation of Absolute Free Energies of Hydration Using Continuum Methods: Accuracy of Partial Charge Models and Optimization of Nonpolar Contributions , 2022 .

[21]  P. Kollman,et al.  Solvation Model Based on Weighted Solvent Accessible Surface Area , 2001 .

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

[23]  W. C. Still,et al.  The GB/SA Continuum Model for Solvation. A Fast Analytical Method for the Calculation of Approximate Born Radii , 1997 .

[24]  Xiaojie Xu,et al.  Empirical Aqueous Solvation Models Based on Accessible Surface Areas with Implicit Electrostatics , 2002 .

[25]  B. Lee,et al.  The interpretation of protein structures: estimation of static accessibility. , 1971, Journal of molecular biology.

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

[27]  Donald G. Truhlar,et al.  PARAMETRIZED MODEL FOR AQUEOUS FREE ENERGIES OF SOLVATION USING GEOMETRY-DEPENDENT ATOMIC SURFACE TENSIONS WITH IMPLICIT ELECTROSTATICS , 1997 .

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

[29]  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..

[30]  M. L. Connolly.  Solvent-accessible surfaces of proteins and nucleic acids. , 1983, Science.

[31]  Vladimir B. Sulimov,et al.  Computation of hydration free energies of organic solutes with an implicit water model , 2006, J. Comput. Chem..

[32]  I. Kuntz,et al.  Inclusion of Solvation in Ligand Binding Free Energy Calculations Using the Generalized-Born Model , 1999 .

[33]  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..

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

[35]  C. Cramer,et al.  New Universal Solvation Model and Comparison of the Accuracy of the SM5.42R, SM5.43R, C-PCM, D-PCM, and IEF-PCM Continuum Solvation Models for Aqueous and Organic Solvation Free Energies and for Vapor Pressures , 2004 .

[36]  Renxiao Wang,et al.  A computational analysis of the binding affinities of FKBP12 inhibitors using the MM‐PB/SA method , 2006, Proteins.

[37]  B. Dominy,et al.  Development of a generalized Born model parameterization for proteins and nucleic acids , 1999 .