Prediction of SAMPL-1 yydration free energies using a continuum electrostatics-dispersion model

  1. (PDF, 473 KB)
  2. Get@NRC: Prediction of SAMPL-1 yydration free energies using a continuum electrostatics-dispersion model (Opens in a new window)
DOIResolve DOI:
AuthorSearch for: ; Search for: ; Search for: ; Search for:
Journal titleThe journal of physical chemistry. B
Pages45114520; # of pages: 10
Subjectanalysis; Biotechnology; Molecular Weight; pha; Prospective Studies; Water
AbstractThe SAMPL-1 hydration free energy blind prediction challenge data set includes 63 compounds that are more chemically diverse, polyfunctional, drug-like, and with examples of transfer free energies and molecular weights larger than ever before seen in previously tabulated data sets of neutral compounds. For the prospective SAMPL-1 study, we employed a continuum model including a boundary element solution of the Poisson equation to describe electrostatic solvation, a molecular surface area-based cost of cavity formation in water, and a continuum Lennard-Jones potential to account for dispersion-repulsion solute-solvent effects. For the latter contribution, continuum van der Waals atom-type coefficients were calibrated and validated on previously available hydration data sets. In the prospective study, this continuum hydration model yielded SAMPL-1 predictions highly correlated with experimental data, albeit with a slope of slightly above 0.5, suggesting a valid model but with a systematic error. Analysis of the major outliers, all overestimating the experimental hydration data, highlights a common structural theme as a possible cause of the prediction errors: densely polar and hydrogen-bond-capable structures, featuring primarily substituted (sulfon)amide groups, often in conjugated systems. By examining analog pairs within the SAMPL-1 data set, it was also noted that certain solvation trends are captured neither by chemical sense nor by our hydration model, which seem too additive. A retrospective analysis of model transferability between hydration data sets as a function of its parameters and complexity indicates that the electrostatic component of the model is fairly transferrable across data sets, but the nonelectrostatic terms are less so. For the chemical space covered in SAMPL-1, absolute prediction errors indicate that the simpler transferrable electrostatics-only model outperforms the more complex model including cavity and continuum dispersion terms. Possible directions to further improve this continuum hydration model are proposed.
Publication date
AffiliationNational Research Council Canada; NRC Biotechnology Research Institute
Peer reviewedNo
NRC number49593
NPARC number12919052
Export citationExport as RIS
Report a correctionReport a correction
Record identifiered74687e-60de-48b6-96ca-41f2f55bf0cf
Record created2009-11-10
Record modified2016-05-09
Bookmark and share
  • Share this page with Facebook (Opens in a new window)
  • Share this page with Twitter (Opens in a new window)
  • Share this page with Google+ (Opens in a new window)
  • Share this page with Delicious (Opens in a new window)