Relative free-energy calculations for scaffold hopping-type transformations with an automated RE-EDS sampling procedure
Abstract
The calculation of relative free-energy differences between different compounds plays an important role in drug design to identify potent binders for a given protein target. Most rigorous methods based on molecular dynamics simulations estimate the free-energy difference between pairs of ligands. Thus, the comparison of multiple ligands requires the construction of a "state graph", in which the compounds are connected by alchemical transformations. The computational cost can be optimized by reducing the state graph to a minimal set of transformations. However, this may require individual adaptation of the sampling strategy if a transformation process does not converge in a given simulation time. In contrast, path-free methods like replica-exchange enveloping distribution sampling (RE-EDS) allow the sampling of multiple states within a single simulation without the pre-definition of alchemical transition paths. To optimize sampling and convergence, a set of RE-EDS parameters needs to be estimated in a pre-processing step. Here, we present an automated procedure for this step that determines all required parameters, improving the robustness and ease of use of the methodology. To illustrate the performance, the relative binding free energies are calculated for a series of checkpoint kinase 1 inhibitors containing challenging transformations in ring size, opening/closing, and extension, which reflect changes observed in scaffold hopping. The simulation of such transformations with RE-EDS can be conducted with conventional force fields and, in particular, without soft bond-stretching terms. Show more
Permanent link
https://doi.org/10.3929/ethz-b-000524065Publication status
publishedExternal links
Journal / series
Journal of Computer-Aided Molecular DesignVolume
Pages / Article No.
Publisher
SpringerSubject
Molecular dynamics; Free energy calculation; Protein-ligand binding; Replica exchange; Enveloping distribution samplingOrganisational unit
09458 - Riniker, Sereina Z. / Riniker, Sereina Z.
Funding
178762 - Passive Membrane-Permeability Prediction for Peptides and Peptidomimetics Using Computational Methods (SNF)
Related publications and datasets
Is source of: https://doi.org/10.3929/ethz-b-000541613
More
Show all metadata