Accurate and reliable calculation of protein-ligand binding affinities remains a hotbed of computer-aided drug design research. Despite the potentially large impact FEP (free energy perturbation) may have in drug design projects, practical applications of FEP in industrial contexts have been limited. In this work, we use a recently developed method, FEP/REST (free energy perturbation/replica exchange with solute tempering), to calculate the relative binding affinities for a set of congeneric ligands binding to the CDK2 receptor. We compare the FEP/REST results with traditional FEP/MD (molecular dynamics) results and MM/GBSA (molecular mechanics/Generalized Born Surface Area model) results and examine why FEP/REST performed notably better than these other methods, as well as why certain ligand mutations lead to large increases of the binding affinity while others do not. We also introduce a mathematical framework for assessing the consistency and reliability of the calculations using cycle closures in FEP mutation paths.
Alchemical free energy calculations hold increasing promise as an aid to drug discovery efforts. However, applications of these techniques in discovery projects have been relatively few, partly because of the difficulty of planning and setting up calculations. Here, we introduce Lead Optimization Mapper, LOMAP, an automated algorithm to plan efficient relative free energy calculations between potential ligands within a substantial library of perhaps hundreds of compounds. In this approach, ligands are first grouped by structural similarity primarily based on the size of a (loosely defined) maximal common substructure, and then calculations are planned within and between sets of structurally related compounds. An emphasis is placed on ensuring that relative free energies can be obtained between any pair of compounds without combining the results of too many different relative free energy calculations (to avoid accumulation of error) and by providing some redundancy to allow for the possibility of error and consistency checking and provide some insight into when results can be expected to be unreliable. The algorithm is discussed in detail and a Python implementation, based on both Schrödinger's and OpenEye's APIs, has been made available freely under the BSD license.
Predicting protein-ligand binding free energies is a central aim of computational structure-based drug design (SBDD)--improved accuracy in binding free energy predictions could significantly reduce costs and accelerate project timelines in lead discovery and optimization. The recent development and validation of advanced free energy calculation methods represents a major step toward this goal. Accurately predicting the relative binding free energy changes of modifications to ligands is especially valuable in the field of fragment-based drug design, since fragment screens tend to deliver initial hits of low binding affinity that require multiple rounds of synthesis to gain the requisite potency for a project. In this study, we show that a free energy perturbation protocol, FEP+, which was previously validated on drug-like lead compounds, is suitable for the calculation of relative binding strengths of fragment-sized compounds as well. We study several pharmaceutically relevant targets with a total of more than 90 fragments and find that the FEP+ methodology, which uses explicit solvent molecular dynamics and physics-based scoring with no parameters adjusted, can accurately predict relative fragment binding affinities. The calculations afford R(2)-values on average greater than 0.5 compared to experimental data and RMS errors of ca. 1.1 kcal/mol overall, demonstrating significant improvements over the docking and MM-GBSA methods tested in this work and indicating that FEP+ has the requisite predictive power to impact fragment-based affinity optimization projects.
Recent
advances in improved force fields and sampling methods have made it
possible for the accurate calculation of protein–ligand binding
free energies. Alchemical free energy perturbation (FEP) using an
explicit solvent model is one of the most rigorous methods to calculate
relative binding free energies. However, for cases where there are
high energy barriers separating the relevant conformations that are
important for ligand binding, the calculated free energy may depend
on the initial conformation used in the simulation due to the lack
of complete sampling of all the important regions in phase space.
This is particularly true for ligands with multiple possible binding
modes separated by high energy barriers, making it difficult to sample
all relevant binding modes even with modern enhanced sampling methods.
In this paper, we apply a previously developed method that provides
a corrected binding free energy for ligands with multiple binding
modes by combining the free energy results from multiple alchemical
FEP calculations starting from all enumerated poses, and the results
are compared with Glide docking and MM-GBSA calculations. From these
calculations, the dominant ligand binding mode can also be predicted.
We apply this method to a series of ligands that bind to c-Jun N-terminal
kinase-1 (JNK1) and obtain improved free energy results. The dominant
ligand binding modes predicted by this method agree with the available
crystallography, while both Glide docking and MM-GBSA calculations
incorrectly predict the binding modes for some ligands. The method
also helps separate the force field error from the ligand sampling
error, such that deviations in the predicted binding free energy from
the experimental values likely indicate possible inaccuracies in the
force field. An error in the force field for a subset of the ligands
studied was identified using this method, and improved free energy
results were obtained by correcting the partial charges assigned to
the ligands. This improved the root-mean-square error (RMSE) for the
predicted binding free energy from 1.9 kcal/mol with the original
partial charges to 1.3 kcal/mol with the corrected partial charges.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.