Abstract:The “embedded cluster reference interaction site model” (EC-RISM) integral equation theory is applied to the problem of predicting aqueous pKa values for drug-like molecules based on an ensemble of tautomers. EC-RISM is based on self-consistent calculations of a solute’s electronic structure and the distribution function of surrounding water. Following-up on the workflow developed after the SAMPL5 challenge on cyclohexane-water distribution coefficients we extended and improved the methodology by taking into account exact electrostatic solute–solvent interactions taken from the wave function in solution. As before, the model is calibrated against Gibbs energies of hydration from the “Minnesota Solvation Database” and a public dataset of acidity constants of organic acids and bases by adjusting in total 4 parameters, among which only 3 are relevant for predicting pKa values. While the best-performing training model yields a root-mean-square error (RMSE) of 1 pK unit, the corresponding test set prediction on the full SAMPL6 dataset of macroscopic pKa values using the same level of theory exhibits slightly larger error (1.7 pK units) than the best test set model submitted (1.7 pK units for corresponding training set vs. test set performance of 1.6). Post-submission analysis revealed a number of physical optimization options regarding the numerical treatment of electrostatic interactions and conformational sampling. While the experimental test set data revealed after submission was not used for reparametrizing the methodology, the best physically optimized models consequentially result in RMSEs of 1.5 if only improved electrostatic interactions are considered and of 1.1 if, in addition, conformational sampling accounts for quantum-chemically derived rankings. We conclude that these numbers are probably near the ultimate accuracy achievable with the simple 3-parameter model using a single or the two best-ranking conformations per tautomer or microstate. Finally, relations of the present macrostate approach to microstate pKa results are discussed and some illustrative results for microstate populations are presented.

Prediction of hydration free energies for aliphatic and aromatic chloro derivatives using molecular dynamics simulations with the OPLS-AA force field

Prediction of hydration free energies for the SAMPL4 diverse set of compounds using molecular dynamics simulations with the OPLS-AA force field

Alchemical prediction of hydration free energies for SAMPL

Prediction of hydration free energies for the SAMPL4 data set with the AMOEBA polarizable force field

Computing hydration free energies of small molecules with first principles accuracy

Extended solvent-contact model approach to SAMPL4 blind prediction challenge for hydration free energies

Exploring the Limits of the Generalized CHARMM and AMBER Force Fields through Predictions of Hydration Free Energy of Small Molecules

Predicting hydration energies for multivalent ions

Machine Learning-Based Enhancements of Empirical Energy Functions: Structure, Dynamics and Spectroscopy of Modified Benzenes

Enhancing Accuracy and Feature Insights in Hydration Free Energy Predictions for Small Molecules with Machine Learning

A Multiscale Procedure for Predicting the Hydration Free Energies of Polycyclic Aromatic Hydrocarbons

Physics-based Machine Learning to Predict Hydration Free Energies for Small Molecules with a minimal number of descriptors: Interpretable and Accurate

Expanded ensemble predictions of toluene--water partition coefficients in the SAMPL9 LogP challenge

Predicting hydration free energies of the FreeSolv database of druglike molecules with molecular density functional theory

Assessment of Halogen Off-Center Point-Charge Models Using Explicit Solvent Simulations

Fast prediction of hydration free energies for SAMPL4 blind test from a classical density functional theory

MolE8: Finding DFT Potential Energy Surface Minima Values from Force-Field Optimised Organic Molecules with New Machine Learning Representations

Machine Learning Prediction of Hydration Free Energy with Physically Inspired Descriptors.

Comparative Assessment of Computational Methods for Free Energy Calculations of Ionic Hydration

The SAMPL6 challenge on predicting aqueous pKa values from EC-RISM theory

A Tiny Charge‐scaling in the OPLS‐AA + L‐OPLS Force Field Delivers the Realistic Dynamics and Structure of Liquid Primary Alcohols