Fragmenstein: predicting protein-ligand structures of compounds derived from known crystallographic fragment hits using a strict conserved-binding–based methodology

26 February 2024, Version 1
This content is a preprint and has not undergone peer review at the time of posting.

Abstract

Current strategies centred on either merging or linking initial hits from fragment-based drug design (FBDD) crystallographic screens ignore 3D structural information. We show that an algorithmic approach (Fragmenstein) that ‘stitches’ the ligand atoms from this structural information together can provide more accurate and reliable predictions for protein-ligand complex conformation than existing methods such as pharmacophore-constrained docking. This approach works under the assumption of conserved binding: when a larger molecule is designed containing the initial fragment hit, the common substructure between the two will adopt the same binding mode. Fragmenstein either takes the coordinates of ligands from a experimental fragment screen and stitches the atoms together to produce a novel merged compound, or uses them to predict the complex for a provided compound. The compound is then energy minimised under strong constraints to obtain a structurally plausible compound. This method is successful in showing the importance of using the coordinates of known binders when predicting the conformation of derivative compounds through a retrospective analysis of the COVID Moonshot data. It has also had a real-world application in hit-to-lead screening, yielding a sub-micromolar merger from parent hits in a single round.

Keywords

Fragment-based drug design
crystallography
physics-based
Open source
Fragment hit
Elaboration
Merging
Linking
Template-based
conformer
RDKit
PyRosetta
Screening

Supplementary materials

Title
Description
Actions
Title
Supplementary figures
Description
Supplementary figure 1. Detailed rules employed in determined atomic overlap between two compounds. Atoms . Supplementary figure 2. Detailed mapping schema used in the placement operation Supplementary figure 3. Distribution of ligand efficiency (left) and of number of interactions per heavy atom for the different merger performed on the Mac1 poised dataset, namely Fragmenstein, Fragmenstein modified to be constrained to a single hit, MCS merger (void of positional information) and BRICS decomposition and building Supplementary figure 4. Placement of NU442 Supplementary figure 5. Example of legitimate merger from Mac1, wherein the acenaphthylene core is chemically sound, but for which no analogues are present in make-on-demand space.
Actions

Supplementary weblinks

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.