Theoretical and Computational Chemistry

Alchemical Absolute Protein-Ligand Binding Free Energies for Drug Design

Authors

Abstract

The recent advances in relative protein-ligand binding free energy calculations have shown the value of alchemical methods in drug discovery. Accurately assessing absolute binding free energies, although highly desired, remains a challenging endeavour, mostly limited to small model cases. Here, we demonstrate accurate first principles based absolute binding free energy estimates for 128 pharmaceutically relevant targets. We use a novel rigorous method to generate protein-ligand ensembles for the ligand in its decoupled state. Not only do the calculations deliver accurate protein-ligand binding affinity estimates, but they also provide detailed physical insight into the structural determinants of binding. We identify subtle rotamer rearrangements between apo and holo states of a protein that are crucial for binding. When compared to relative binding free energy calculations, obtaining absolute binding free energies is considerably more challenging in large part due to the need to explicitly account for the protein in its apo state. In this work we present several approaches to obtain apo state ensembles for accurate absolute ΔG calculations, thus outlining protocols for prospective application of the methods for drug discovery.

Version notes

Addressed reviewer comments, added discussion of effects of using crystallographic water positions to SI.

Content

Thumbnail image of revised_main.pdf

Supplementary material

Thumbnail image of revised_si.pdf
Supplementary Information
Outline of the NEQ method, relative binding free energies recalculated from absolute ones broken down by system, absolute binding free energies calculated using apo and holo structures for the decoupled state broken down by system, p38 loop motion and its effects on absolute binding free energy, effects of extended equilibration on tyk2 absolute binding free energies, comparison of uncertainties, descriptions of restraint generation and decorrelation methods, validation of decorelation approaches, details of parametrization and simulation conditions, numerical values for computed and experimental absolute binding free energies.