Abstract:A principal goal of drug discovery project is to design molecules that can tightly and selectively bind to the target protein receptor. Accurate prediction of protein-ligand binding free energies is therefore of central importance in computational chemistry and computer aided drug design. Multiple recent improvements in computing power, classical force field accuracy, enhanced sampling methods, and simulation setup have enabled accurate and reliable calculations of protein-ligands binding free energies, and position free energy calculations to play a guiding role in small molecule drug discovery. In this Account, we outline the relevant methodological advances, including the REST2 (Replica Exchange with Solute Temperting) enhanced sampling, the incorporation of REST2 sampling with convential FEP (Free Energy Perturbation) through FEP/REST, the OPLS3 force field, and the advanced simulation setup that constitute our FEP+ approach, followed by the presentation of extensive comparisons with experiment, demonstrating sufficient accuracy in potency prediction (better than 1 kcal/mol) to substantially impact lead optimization campaigns. The limitations of the current FEP+ implementation and best practices in drug discovery applications are also discussed followed by the future methodology development plans to address those limitations. We then report results from a recent drug discovery project, in which several thousand FEP+ calculations were successfully deployed to simultaneously optimize potency, selectivity, and solubility, illustrating the power of the approach to solve challenging drug design problems. The capabilities of free energy calculations to accurately predict potency and selectivity have led to the advance of ongoing drug discovery projects, in challenging situations where alternative approaches would have great difficulties. The ability to effectively carry out projects evaluating tens of thousands, or hundreds of thousands, of proposed drug candidates, is potentially transformative in enabling hard to drug targets to be attacked, and in facilitating the development of superior compounds, in various dimensions, for a wide range of targets. More effective integration of FEP+ calculations into the drug discovery process will ensure that the results are deployed in an optimal fashion for yielding the best possible compounds entering the clinic; this is where the greatest payoff is in the exploitation of computer driven design capabilities. A key conclusion from the work described is the surprisingly robust and accurate results that are attainable within the conventional classical simulation, fixed charge paradigm. No doubt there are individual cases that would benefit from a more sophisticated energy model or dynamical treatment, and properties other than protein-ligand binding energies may be more sensitive to these approximations. We conclude that an inflection point in the ability of MD simulations to impact drug discovery has now been attained, due to the confluence of hardware and software development along with the formulation of "good enough" theoretical methods and models.

Nonbonded Force Field Parameters from MBIS Partitioning of the Molecular Electron Density Improve Binding Affinity Predictions of the T4-Lysozyme Double Mutant

Non-bonded force field parameters from MBIS partitioning of the molecular electron density improve binding affinity predictions of the T4-lysozyme double mutant

Biomolecular conformational changes and ligand binding: from kinetics to thermodynamics.

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

Pbsa_e: A Pbsa-Based Free Energy Estimator for Protein-Ligand Binding Affinity

Estimating the Roles of Protonation and Electronic Polarization in Absolute Binding Affinity Simulations

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

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

Ligand affinities predicted with the MM/PBSA method: dependence on the simulation method and the force field

Accurate Calculation of Relative Binding Free Energies Between Ligands with Different Net Charges

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

Accuracy Of Alternate Nonpolarizable Force Fields For The Determination Of Protein-Ligand Binding Affinities Dominated By Cation-Pi Interactions

How Well Does the Extended Linear Interaction Energy Method Perform in Accurate Binding Free Energy Calculations?

Protein-ligand free energies of binding from full-protein DFT calculations: convergence and choice of exchange-correlation functional

Calculation of Protein–ligand Binding Affinities Based on a Fragment Quantum Mechanical Method

Improving the Performance of MM/PBSA in Protein-Protein Interactions Via the Screening Electrostatic Energy.

An Accurate Free Energy Estimator: Based on MM/PBSA Combined with Interaction Entropy for Protein-Ligand Binding Affinity.

VAD-MM/GBSA: A Variable Atomic Dielectric MM/GBSA Model for Improved Accuracy in Protein–Ligand Binding Free Energy Calculations

Accurate and Efficient Calculation of Protein-Protein Binding Free Energy-Interaction Entropy with Residue Type-Specific Dielectric Constants.

Improving the accuracy of predicting protein-ligand binding-free energy with semiempirical quantum chemistry charge.

Advancing Drug Discovery through Enhanced Free Energy Calculations