Core Flipping in Lead Optimization: Rank Ordering using Multisite λ Dynamics

10 October 2024, Version 2
This content is a preprint and has not undergone peer review at the time of posting.

Abstract

In structure-based drug discovery, reliable structural models of ligands bound to their target receptor are critical for establishing the structure activity relationship of congeneric series. In such a series, substitutions on a common scaffold core might lead to different binding modes, ranging from slight changes of orientations to flipping or inversion of the core structure. Moreover, molecular docking might lead to alternative orientations within the top ranked poses without being able to discriminate which is the most likely. To determine the relative binding affinities between two alternative ligand poses, we propose a methodology based on relative binding free energy calculations using the multisite λ dynamics method. We use a dual-topology approach with distance restraining schemes. We introduced a novel strategy using one-step perturbation to calculate the contributions of the applied restraints. While using FEP/MBAR instead for that purpose led to smaller uncertainties it suffered from convergence issues. We test the validity and predictive power of our approach using two pharmaceutically relevant targets and eight compounds from experimentally characterized congeneric series. For each target, our approach correctly ranks the known X-ray poses to be more favorable than alternative flipped poses. The proposed methodology can be easily extended to rank more than two poses in a single simulation and should also be applicable for the prediction of conformational binding cliffs, ranking fragments, core hopping cases and evaluating different scaffold binding to a common target.

Keywords

Core flipping
multisite λ dynamics
relative binding free energy
one-step perturbation

Supplementary materials

Title
Description
Actions
Title
Supporting Information text
Description
Supporting information: Details for small molecule parameters, molecular docking, windowing approach for removing restraints; results from docking, ligand RMSD from equilibrium MD simulations, ALF convergence plots, interaction analysis, windowing approach and ligand RMSF from the OSP methodology.
Actions

Comments

Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.