Abstract:Abstract Alchemical free energy calculations are an important tool in the computational chemistry toolbox, enabling the efficient calculation of quantities critical for drug discovery such as ligand binding affinities, selectivities, and partition coefficients. However, modern alchemical free energy calculations suffer from three significant limitations: (1) modern molecular mechanics force fields are limited in their ability to model complex molecular interactions, (2) classical force fields are unable to treat phenomena that involve rearrangements of chemical bonds, and (3) these calculations are unable to easily learn to improve their performance if readily-available experimental data is available. Here, we show how all three limitations can be overcome through the use of quantum machine learning (QML) potentials capable of accurately modeling quantum chemical energetics even when chemical bonds are made and broken. Because these potentials are based on mathematically convenient deep learning architectures instead of traditional quantum chemical formulations, QML simulations can be run at a fraction of the cost of quantum chemical simulations using modern graphics processing units (GPUs) and machine learning frameworks. We demonstrate that alchemical free energy calculations in explicit solvent are especially simple to implement using QML potentials because these potentials lack singularities and other pathologies typical of molecular mechanics potentials, and that alchemical free energy calculations are highly effective even when bonds are broken or made. Finally, we show how a limited number of experimental free energy measurements can be used to significantly improve the accuracy of computed free energies for unrelated compounds with no significant generalization gap. We illustrate these concepts on the prediction of aqueous tautomer free energies (related to tautomer ratios), which are highly relevant to drug discovery in that more than a quarter of all approved drugs exist as a mixture of tautomers.

Addressing Suboptimal Poses in Nonequilibrium Alchemical Calculations

Calculation of Relative Binding Free Energy for Mutations in Protein Complexes: The Alchemical Path

Perturbation Potentials to Overcome Order/Disorder Transitions in Alchemical Binding Free Energy Calculations

Alchemical Transfer Approach to Absolute Binding Free Energy Estimation

I. Dissociation free energies of drug-receptor systems via non-equilibrium alchemical simulations: a theoretical framework.

A Method for Treating Significant Conformational Changes in Alchemical Free Energy Simulations of Protein-Ligand Binding

Validation of the Alchemical Transfer Method for the Estimation of Relative Binding Affinities of Molecular Series

Automated Adaptive Absolute Binding Free Energy Calculations

Multistate Method to Efficiently Account for Tautomerism and Protonation in Alchemical Free-Energy Calculations

Equilibrium and Non-equilibrium Ensemble Methods for Accurate, Precise and Reproducible Absolute Binding Free Energy Calculations

Automated Adaptive Absolute Binding Free Energy Calculations

Teaching free energy calculations to learn from experimental data

Relative Binding Free Energy Calculations for Ligands with Diverse Scaffolds with the Alchemical Transfer Method

Lambda-ABF: Simplified, Portable, Accurate and Cost-effective Alchemical Free Energy Computations

Automated relative binding free energy calculations from SMILES to ΔΔG

On achieving high accuracy and reliability in the calculation of relative protein–ligand binding affinities

Combining Alchemical Transformation with Physical Pathway to Accurately Compute Absolute Binding Free Energy

Large scale relative protein ligand binding affinities using non-equilibrium alchemy

On analytical corrections for restraints in absolute binding free energy calculations

Alchemical Response Parameters from an Analytical Model of Molecular Binding