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.

Calculation of Solute Partition Coefficient Using the A-P Scheme

Prediction of Liquid-Liquid Equilibria of Polymersalt Aqueous Two-Phase Systems by a Modified Pitzer's Virial Equation

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

An Approach to Computing Solvent Reorganization Energy.

Predicting Octanol/water Partition Coefficient Using Solvation Free Energy and Solvent-Accessible Surface Area.

Estimation of Molecular Linear Free Energy Relationship Descriptors by a Group Contribution Approach. 2. Prediction of Partition Coefficients

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

A New Method for Calculating Partition Coefficients of Organic Compounds

Estimating protein-ligand binding free energy: atomic solvation parameters for partition coefficient and solvation free energy calculation.

Empirical Aqueous Solvation Models Based on Accessible Surface Areas with Implicit Electrostatics

Guidelines for descriptor assignments for the solvation parameter model by separation techniques

Model for Aqueous Solvation Based on Class IV Atomic Charges and First Solvation Shell Effects

A New Atom-Additive Method for Calculating Partition Coefficients

Assessment of liquid-liquid partition for the assignment of descriptors for the solvation parameter model

Solvation Entropy Made Simple

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

GAFF/IPolQ-Mod+LJ-Fit: Optimized force field parameters for solvation free energy predictions

Comparison of Implicit and Explicit Solvent Models for the Calculation of Solvation Free Energy in Organic Solvents.

Molecular Partition Coefficient from Machine Learning with Polarization and Entropy Embedded Atom- Centered Symmetry Functions

An Extended Aqueous Solvation Model Based on Atom-Weighted Solvent Accessible Surface Areas: SAWSA V2.0 Model

Ion solvation free energy calculations based on first-principles molecular dynamics thermodynamic integration