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.

Computation of Rate Coefficients in Solutions Based on Transition State Theory Combined with a Heuristically Corrected Polarizable Continuum Model: Intermolecular Diels-Alder Reactions as Case Studies

A simple heuristic approach to estimate the thermochemistry of condensed-phase molecules based on the polarizable continuum model

Comparison of DFT, MP2/CBS, and CCSD(T)/CBS methods for a dual-level QM/MM Monte Carlo simulation approach calculating the free energy of activation of reactions in solution and “on water”: a case study

Transition State Theory within Constrained Nuclear-Electronic Orbital Framework: Accurate Hydrogen Atom Transfer Reaction Rates with Direct Incorporation of Quantum Nuclear Delocalization Effects

Correlation in Quantum Chemical Calculation and Its Effect on the Uncertainty of Theoretically Predicted Rate Coefficients and Branching Ratios

Theoretical Study on the Decomposition of Dmp and Dep in Condensed Phase : Solvent Effect and Tunneling Effect

Kinetics and Thermodynamics of Reactions Involving Criegee Intermediates: An Assessment of Density Functional Theory and Ab Initio Methods Through Comparison with CCSDT(Q)/CBS Data

Rate Coefficients of the Cl + CH3C(O)OCH3 → HCl + CH3C(O)OCH2 Reaction at Different Temperatures Calculated by Transition-State Theory with Ab Initio and Density Functional Theory Reaction Paths

Rate Coefficients of the CF3CHFCF3 + H → CF3CFCF3 + H2 Reaction at Different Temperatures Calculated by Transition State Theory with Ab Initio and DFT Reaction Paths

Quantum-classical Rate Coefficient Datasets of Vibrational Energy Transfer in Carbon Monoxide Based on Highly Accurate Potential Energy Surface.

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

A Method Of Calculating The Kamlet-Abboud-Taft Solvatochromic Parameters Using Cosmo-Rs

Heavy Atom Tunneling in Organic Reactions at Coupled Cluster Potential Accuracy with a Parallel Implementation of Anharmonic Constant Calculations and Semiclassical Transition State Theory

Multistate, polarizable QM/MM embedding scheme based on the direct reaction field method: Solvatochromic shifts, analytical gradients and optimization of conical intersections in solution

Constrained Broyden Minimization Combined with the Dimer Method for Locating Transition State of Complex Reactions

Quantum Mechanical Calculation of Noncovalent Interactions: A Large-Scale Evaluation of PMx, DFT, and SAPT Approaches

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

Elements of Transition-State Theory in Relation to the Thermal Dissociation of Selected Solid Compounds

Constrained Nuclear-Electronic Orbital Transition State Theory Using Energy Surfaces with Nuclear Quantum Effects

Validation of electronic structure methods for isomerization reactions of large organic molecules

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