Abstract:Methods for calculating the relative binding free energy (RBFE) between ligands to a target protein are gaining importance in the structure-based drug discovery domain, especially as methodological advances and automation improve accuracy and ease of use. In an RBFE calculation, the difference between the binding affinities of two ligands to a protein is calculated by transforming one ligand into another, in the protein-ligand complex, and in solvent. Alchemical binding free energy calculations are often used for such ligand transformations. Such calculations are not without challenges, however; for example, it can be challenging to handle interfacial waters when these play a crucial role in mediating protein-ligand binding. In some cases, the exchange of the interfacial waters with solvent water might be very infrequent in the course of typical molecular simulations, and such interfacial waters can be considered trapped on the simulation timescale. In these cases, RBFE calculation between two ligands, where one ligand binds with a trapped water while the other ligand displaces it, can result in inaccuracies if the surrounding water structure is not sampled adequately for both ligands. So far, a popular choice for treating the trapped waters in RBFE calculations is to combine free energy calculations with enhanced sampling methods that insert/delete waters in the binding site. Despite recent developments in the enhanced sampling methods, they can result in hysteresis in the RBFE estimate, depending on whether the simulations were started with or without the trapped waters. In this study, we introduce an alternative method, separation of states, to calculate the RBFE between ligand pairs where the ligands bind to the protein with different numbers/positions of trapped waters. The separation of states approach treats the sampling of the trapped waters separately from the free energy calculation of the ligand transformation. In our method, a trapped water in protein's binding site is decoupled from the system first, and the cavity created by its decoupling is stabilized. We then grow a larger ligand into this cavity-- a ligand that is known to displace the trapped water. In this study, we show that our method results in precise and accurate estimates of RBFEs for ligand pairs involving the rearrangement of trapped water via RBFE calculations for five such ligand pairs. We have optimized our simulation protocol to be suited for large distributed computational resources and have automated our RBFE calculation workflow.

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

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

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

What to Make of Zero: Resolving the Statistical Noise from Conformational Reorganization in Alchemical Binding Free Energy Estimates with Metadynamics Sampling

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

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

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

Performance and Analysis of the Alchemical Transfer Method for Binding Free Energy Predictions of Diverse Ligands

Leveraging a separation of states method for relative binding free energy calculations in systems with trapped waters

Understanding Aldose Reductase-Inhibitors interactions with free energy simulation

Studying the Collective Functional Response of a Receptor in Alchemical Ligand Binding Free Energy Simulations with Accelerated Solvation Layer Dynamics

Equilibrium and Nonequilibrium Ensemble Methods for Accurate, Precise and Reproducible Absolute Binding Free Energy Calculations

Improvements in Precision of Relative Binding Free Energy Calculations Afforded by the Alchemical Enhanced Sampling (ACES) Approach

Accurate Estimation Of Protein-Ligand Binding Free Energies Based On Geometric Restraints

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

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

BFEE2: Automated, Streamlined, and Accurate Absolute Binding Free-Energy Calculations

BFEE: A User-Friendly Graphical Interface Facilitating Absolute Binding Free-Energy Calculations.

Scalable Indirect Free Energy Method Applied to Divalent Cation-Metalloprotein Binding

Alchemical Binding Free Energy Calculations in AMBER20: Advances and Best Practices for Drug Discovery