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.

New evaluation of reconstructed spatial distribution function from radial distribution functions.

A new method to reconstruct three-dimensional spatial distribution function from radial distribution function in solvation structure.

Recent developments and applications of reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED): A hybrid model of quantum chemistry and integral equation theory of molecular liquids

N-representable one-electron reduced density matrix reconstruction with frozen core electrons

An independent, general method for checking consistency between diffraction data and partial radial distribution functions derived from them: the example of liquid water

Scaled Alternating Steepest Descent Algorithm Applied for Protein Structure Determination from Nuclear Magnetic Resonance Data.

Simple and accurate expressions for radial distribution functions of hard disk and hard sphere fluids

A new method on reconstructing protein structure from NOESY distances

Accelerating the 3D-RISM theory of molecular solvation with treecode summation and cut-offs

3DRISM-HI-D2MSA: an improved analytic theory to compute solvent structure around hydrophobic solutes with proper treatment of solute–solvent electrostatic interactions

Solvent Distribution Effects on Quantum Chemical Calculations with Quantum Computers

Spatial distribution of reduced density of hard spheres near a hard-sphere dimer: Results from three-dimensional Ornstein-Zernike equations coupled with several different closures and from grand canonical Monte Carlo simulation

Modeling solvation effects in real-space and real-time within Density Functional Approaches

A comparative study of aqueous DMSO mixtures by computer simulations and integral equation theories

Evaluation of Reverse Monte Carlo Models based on Molecular Dynamics Simulations: A Case Study of Ion Conducting Network Glasses

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

Efficient structural relaxation based on the random phase approximation: Applications to the water clusters

Integral equation models for solvent in macromolecular crystals

A solution for an inverse problem in liquid AFM: calculation of three-dimensional solvation structure on a sample surface

Vibrational Self-Consistent Field (VSCF) and Post-VSCF Method Calculations Combined with the Reference Interaction Site Model Self-Consistent Field Method Coupled with the Constrained Spatial Electron Density Distribution: Applications to NaHCOO in Aqueous Phase