Abstract:Implicit solvent methods for classical molecular modeling are frequently used to provide fast, physics-based hydration free energies of macromolecules. Less commonly considered is the transferability of these methods to other solvents. The Statistical Assessment of Modeling of Proteins and Ligands 5 (SAMPL5) distribution coefficient dataset and the accompanying explicit solvent partition coefficient reference calculations provide a direct test of solvent model transferability. Here we use the 3D reference interaction site model (3D-RISM) statistical-mechanical solvation theory, with a well tested water model and a new united atom cyclohexane model, to calculate partition coefficients for the SAMPL5 dataset. The cyclohexane model performed well in training and testing (R=0.98\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$R=0.98$$\end{document} for amino acid neutral side chain analogues) but only if a parameterized solvation free energy correction was used. In contrast, the same protocol, using single solute conformations, performed poorly on the SAMPL5 dataset, obtaining R=0.73\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$R=0.73$$\end{document} compared to the reference partition coefficients, likely due to the much larger solute sizes. Including solute conformational sampling through molecular dynamics coupled with 3D-RISM (MD/3D-RISM) improved agreement with the reference calculation to R=0.93\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$R=0.93$$\end{document}. Since our initial calculations only considered partition coefficients and not distribution coefficients, solute sampling provided little benefit comparing against experiment, where ionized and tautomer states are more important. Applying a simple pKa\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\hbox {p}K_{\text {a}}$$\end{document} correction improved agreement with experiment from R=0.54\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$R=0.54$$\end{document} to R=0.66\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$R=0.66$$\end{document}, despite a small number of outliers. Better agreement is possible by accounting for tautomers and improving the ionization correction.

The partition of organic compounds from water into the methyl isobutyl ketone extraction solvent with updated Abraham model equation

Updated Abraham model correlations for describing solute transfer into both 2-Pentanol and 3-Methyl-1-butanol based on much larger data sets

Solubility of Isophthalic Acid in Propyl Acetate and Partition Coefficient Between Propyl Acetate and Water

Abraham model correlations for describing solute transfer into 4-methyl-2-pentanol from both water and the gas phase

Development of Abraham model correlations for solute transfer into 2-ethyl-1-hexanol from both water and the gas phase based on measured solubility ratios

Abraham General Solvation Parameter Model: Predictive Expressions for Solute Transfer into Isobutyl Acetate

Abraham solvation parameter model: determination of experiment-based solute descriptor values for 3,5-dimethoxybenzoic acid based on measured solubility data

Abraham model description of the solubilising properties of the isopropyl acetate organic mono-solvent

Experiment-based Abraham model solute descriptors for 2,3-dimethoxybenzoic acid from measured solubility data

Updated Abraham model correlations for enthalpies of solvation of organic solutes dissolved in benzene and acetonitrile

Evaluation of the Goss-modified solvation parameter model for the characterization of biphasic systems and descriptor assignments

Solubility of 4-methyl-3-nitrobenzoic acid in organic mono-solvents: calculation of Abraham model solute descriptors

Development of Abraham model correlations for predicting the solubility of crystalline organic compounds and inorganic gases in ethyl formate

Experiment-based Abraham model solute descriptors for 2‐[4-(dibutylamino)-2-hydroxybenzoyl]benzoic acid

Collection of Partition Coefficients in Hexadecyltrimethylammonium Bromide, Sodium Cholate, and Lithium Perfluorooctanesulfonate Micellar Solutions: Experimental Determination and Computational Predictions

SAMPL5: 3D-RISM partition coefficient calculations with partial molar volume corrections and solute conformational sampling

Calculating the Partition Coefficients of Organic Solvents in Octanol/Water and Octanol/Air

Calculation of Relative Solubility of Semipolar Solvents by Abraham Solvation Parameter Model for Extractables and Leachables Analysis in Chemical Characterization of Medical Devices

Determination of Abraham model solute descriptors for o-acetoacetanisidide based on experimental solubility data in organic mono-solvents

A New Method for Calculating Partition Coefficients of Organic Compounds

Gas-to-ionic Liquid Partition: QSPR Modeling and Mechanistic Interpretation.