Abstract:Accurate in silico prediction of protein-ligand binding affinities has been a primary objective of structure-based drug design for decades due to the putative value it would bring to the drug discovery process. However, computational methods have historically failed to deliver value in real-world drug discovery applications due to a variety of scientific, technical, and practical challenges. Recently, a family of approaches commonly referred to as relative binding free energy (RBFE) calculations, which rely on physics-based molecular simulations and statistical mechanics, have shown promise in reliably generating accurate predictions in the context of drug discovery projects. This advance arises from accumulating developments in the underlying scientific methods (decades of research on force fields and sampling algorithms) coupled with vast increases in computational resources (graphics processing units and cloud infrastructures). Mounting evidence from retrospective validation studies, blind challenge predictions, and prospective applications suggests that RBFE simulations can now predict the affinity differences for congeneric ligands with sufficient accuracy and throughput to deliver considerable value in hit-to-lead and lead optimization efforts. Here, we present an overview of current RBFE implementations, highlighting recent advances and remaining challenges, along with examples that emphasize practical considerations for obtaining reliable RBFE results. We focus specifically on relative binding free energies because the calculations are less computationally intensive than absolute binding free energy (ABFE) calculations and map directly onto the hit-to-lead and lead optimization processes, where the prediction of relative binding energies between a reference molecule and new ideas (virtual molecules) can be used to prioritize molecules for synthesis. We describe the critical aspects of running RBFE calculations, from both theoretical and applied perspectives, using a combination of retrospective literature examples and prospective studies from drug discovery projects. This work is intended to provide a contemporary overview of the scientific, technical, and practical issues associated with running relative binding free energy simulations, with a focus on real-world drug discovery applications. We offer guidelines for improving the accuracy of RBFE simulations, especially for challenging cases, and emphasize unresolved issues that could be improved by further research in the field.

Hit Identification and Binding Mode Predictions by Rigorous Free Energy Simulations

Predicting Fragment Binding Modes Using Customized Lennard-Jones Potentials in Short Molecular Dynamics Simulations

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

Relative Binding Free Energy Calculations in Drug Discovery: Recent Advances and Practical Considerations

Prediction of protein–ligand binding affinity by free energy simulations: assumptions, pitfalls and expectations

Binding Free Energy Calculation Based on the Fragment Molecular Orbital Method and Its Application in Designing Novel SHP-2 Allosteric Inhibitors

Rapid and Reliable Binding Affinity Prediction of Bromodomain Inhibitors: A Computational Study

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

Rigorous Free Energy Perturbation Approach to Estimating Relative Binding Affinities Between Ligands with Multiple Protonation and Tautomeric States.

Rigorous Free Energy Simulations in Virtual Screening

PopShift: A Thermodynamically Sound Approach to Estimate Binding Free Energies by Accounting for Ligand-Induced Population Shifts from a Ligand-Free Markov State Model

Accuracy and precision of binding free energy prediction for a tacrine related lead inhibitor of acetylcholinesterase with an arsenal of supercomputerized molecular modelling methods: a comparative study

Implicit Ligand Theory: Rigorous Binding Free Energies and Thermodynamic Expectations from Molecular Docking

Fast and accurate predictions of binding free energies using MM-PBSA and MM-GBSA

Prediction of ligand binding mode among multiple cross-docking poses by molecular dynamics simulations

Quantum Mechanics-Based Fast and Reliable Prediction of Binding Pose Structures

Can molecular dynamics simulations improve predictions of protein-ligand binding affinity with machine learning?

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

Modeling Local Structural Rearrangements Using FEP/REST: Application to Relative Binding Affinity Predictions of CDK2 Inhibitors

True Accuracy of Fast Scoring Functions to Predict High-Throughput Screening Data from Docking Poses: The Simpler the Better

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