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.

LiProS: FAIR simulation workflow to Predict Accurate Lipophilicity Profiles for Small Molecules

LiProS: FAIR simulation workflow to Predict Accurate Lipophilicity Profiles for Small Molecules

Developing the QSPR Model for Predicting the Storage Lipid/water Distribution Coefficient of Organic Compounds

Filling the gap in LogP and pK_a evaluation for saturated fluorine-containing derivatives with machine learning

A comparison of molecular representations for lipophilicity quantitative structure–property relationships with results from the SAMPL6 logP Prediction Challenge

ChemFlow─From 2D Chemical Libraries to Protein-Ligand Binding Free Energies

Recommendations for Accurate Lipid Annotation and Semi-absolute Quantification from LC-MS/MS Datasets

ALipSol: An Attention-Driven Mixture-of-Experts Model for Lipophilicity and Solubility Prediction

Modeling of small molecule's affinity to phospholipids using IAM-HPLC and QSRR approach enhanced by similarity-based machine algorithms

A Validated Method for Predicting Small Molecule Ionization Sites using Gibb's Free Energies

Rapid proteome-wide prediction of lipid-interacting proteins through ligand-guided structural genomics

FaissMolLib: An efficient and easy deployable tool for ligand-based virtual screening

SMPLIP-Score: predicting ligand binding affinity from simple and interpretable on-the-fly interaction fingerprint pattern descriptors

ProLIF: a library to encode molecular interactions as fingerprints

Deriving Accurate Lipid Classification based on Molecular Formula

Machine Learning for Fast, Quantum Mechanics-Based Approximation of Drug Lipophilicity

Computational Prodrug Design Methodology for Liposome Formulability Enhancement of Small-Molecule APIs

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

FATSLiM: a fast and robust software to analyze MD simulations of membranes

Lipid Annotator: Towards Accurate Annotation in Non-Targeted Liquid Chromatography High-Resolution Tandem Mass Spectrometry (LC-HRMS/MS) Lipidomics Using a Rapid and User-Friendly Software

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