Abstract:Abstract Alchemical free energy methods with molecular mechanics (MM) force fields are now widely used in the prioritization of small molecules for synthesis in structure-enabled drug discovery projects because of their ability to deliver 1–2 kcal mol −1 accuracy in well-behaved protein-ligand systems. Surpassing this accuracy limit would significantly reduce the number of compounds that must be synthesized to achieve desired potencies and selectivities in drug design campaigns. However, MM force fields pose a challenge to achieving higher accuracy due to their inability to capture the intricate atomic interactions of the physical systems they model. A major limitation is the accuracy with which ligand intramolecular energetics—especially torsions—can be modeled, as poor modeling of torsional profiles and coupling with other valence degrees of freedom can have a significant impact on binding free energies. Here, we demonstrate how a new generation of hybrid machine learning / molecular mechanics (ML/MM) potentials can deliver significant accuracy improvements in modeling protein-ligand binding affinities. Using a nonequilibrium perturbation approach, we can correct a standard, GPU-accelerated MM alchemical free energy calculation in a simple post-processing step to efficiently recover ML/MM free energies and deliver a significant accuracy improvement with small additional computational effort. To demonstrate the utility of ML/MM free energy calculations, we apply this approach to a benchmark system for predicting kinase:inhibitor binding affinities—a congeneric ligand series for non-receptor tyrosine kinase TYK2 (Tyk2)—wherein state-of-the-art MM free energy calculations (with OPLS2.1) achieve inaccuracies of 0.93±0.12 kcal mol −1 in predicting absolute binding free energies. Applying an ML/MM hybrid potential based on the ANI2x ML model and AMBER14SB/TIP3P with the OpenFF 1.0.0 (“Parsley”) small molecule force field as an MM model, we show that it is possible to significantly reduce the error in absolute binding free energies from 0.97 [95% CI: 0.68, 1.21] kcal mol −1 (MM) to 0.47 [95% CI: 0.31, 0.63] kcal mol −1 (ML/MM).

Microkinetic modeling with machine learning predicted binding energies of reaction intermediates of ethanol steam reforming: The limitations

Quantifying the error propagation in microkinetic modeling of catalytic reactions with model-predicted binding energies

Modeling of Transition State by Molecular Dynamics. Prediction of Catalytic Efficiency of the Mutants of Mandelate Racemase

The importance of reaction energy in predicting chemical reaction barriers with machine learning models

Assessment of machine learning models trained by molecular dynamics simulations results for inferring ethanol adsorption on an aluminium surface

Interpretable Machine Learning Model for Activation Energy Prediction Based on Biomass Properties

Using machine learning to go beyond potential energy surface benchmarking for chemical reactivity

Prediction of Hydration energies of Adsorbates at Pt(111) and Liquid Water Interfaces using Machine Learning

Predicting Chemical Reaction Barriers with a Machine Learning Model

Prediction of energies for reaction intermediates and transition states on catalyst surfaces using graph-based machine learning models

Beyond potential energy surface benchmarking: a complete application of machine learning to chemical reactivity

Barrier Energies ( in eV ) for the Reactions in the Microkinetic Model a

Machine learning applications for thermochemical and kinetic property prediction

Understanding the Source of Error in First-Principles-Based Micro-Kinetic Modeling: Density Functional Theory Calculations Versus the Mean-Field Approximation

Modeling Chemical Processes in Explicit Solvents with Machine Learning Potentials

Towards chemical accuracy for alchemical free energy calculations with hybrid physics-based machine learning / molecular mechanics potentials

Forward and reverse uncertainty analyses for RRKM/master equation based kinetic predictions: A case study of ethyl with oxygen

Ultrastructural study of the action of aminopiperazines on transformed cell lines and on nervous system tissue culture.

Machine-Learning-Assisted Free Energy Simulation of Solution-Phase and Enzyme Reactions

A machine learning model to predict the pyrolytic kinetics of different types of feedstocks

Dynamic hybrid modeling of fuel ethanol fermentation process by integrating biomass concentration XGBoost model and kinetic parameter artificial neural network model into mechanism model