Bio Covalent Design

Version Compatibility

Reference examples tested with: RDKit 2024.09+, OpenEye / AutoDock Vina 1.2+ (for covalent extensions), GOLD (commercial), DOCKovalent (web service), HCovDock 1.0+.

Before using code patterns, verify installed versions match. If versions differ:

• Python: pip show rdkit then help(rdkit.Chem) to check signatures
• CLI: check version output of each docking tool

If code throws ImportError, AttributeError, or TypeError, introspect the installed package and adapt the example to match the actual API rather than retrying.

Covalent Inhibitor Design

Design molecules that form covalent bonds with target protein residues. The "covalent revolution" (Lonsdale & Ward 2018) made TCIs (Targeted Covalent Inhibitors) clinically validated: KRAS G12C inhibitors (sotorasib, adagrasib), BTK inhibitors (ibrutinib), and EGFR inhibitors (osimertinib) are recent successes. Postdoc-grade covalent design requires balancing intrinsic reactivity (must form bond) vs selectivity (only the intended residue), reversibility (irreversible vs reversible covalent), and drug-likeness (warheads can hurt PK).

For warhead substructure filtering (in non-covalent contexts), see chemoinformatics/substructure-search. For non-covalent docking, see chemoinformatics/virtual-screening. For pose validation, see chemoinformatics/pose-validation.

Reactive Residue Taxonomy

Cysteine is the dominant target because:

• Soft nucleophile (matches soft electrophiles)
• Low background reactivity (rare in proteins, ~1.7%)
• Distinguishable from common nucleophiles (GSH, off-target Cys)

Warhead Chemistry

Practical hierarchy: Acrylamide is the modern default for cysteine-selective TCIs (KRAS G12C, EGFR, BTK). Chloroacetamide is more reactive (faster) but less selective.

Decision Tree by Scenario

Kinetics: kinact / Ki

Covalent inhibition kinetics:

• Ki: reversible binding affinity (initial, like non-covalent IC50)
• kinact: rate of covalent bond formation (sec^-1)
• kinact/Ki: second-order rate constant, "covalent efficiency" (M^-1 s^-1)

Modern best practice: report kinact/Ki, not just IC50. Two compounds with same IC50 can have very different kinact/Ki:

• Low Ki, low kinact: tight binding, slow covalent bond
• High Ki, high kinact: loose binding, fast covalent bond

Intrinsic Reactivity Assays

Before committing to a warhead, measure intrinsic reactivity (off-target risk):

# Generic GSH stability assay readout - measure half-life of warhead with 10 mM GSH
# kinact_GSH from time-course of warhead disappearance

The "GSH-stable" warhead (t1/2 > 4 hours) is the modern target for druglike TCIs.

Covalent Docking Tools

For open-source covalent docking, HCovDock (2023) is the modern alternative; DOCKovalent is the longstanding standard.

Example: KRAS G12C Inhibitor Design Workflow

Goal: Decorate a co-crystal scaffold with a cysteine-targeting warhead and rank candidates by covalent efficiency.

Approach: Load scaffold SMILES, enumerate acrylamide-bearing analogs, filter by reactivity selectivity, dock under covalent constraint, and rank by kinact/Ki surrogates.

from rdkit import Chem

# Step 1: scaffold from co-crystal (4LRW or AMG510)
scaffold_smi = 'c1ccc(C(=O)NC2=Nc3c(...)cnc23)cc1'
scaffold = Chem.MolFromSmiles(scaffold_smi)

# Step 2: enumerate analogs with acrylamide warhead
def add_acrylamide(scaffold, attachment_atom_idx):
    """Append acrylamide (-NC(=O)C=C) at a hydrogen position"""
    warhead = Chem.MolFromSmiles('NC(=O)C=C')
    # ... combine via Chem.RWMol or fragment combination
    pass

# Step 3: filter for reactive group selectivity
# Step 4: dock with DOCKovalent / GOLD covalent / HCovDock
# Step 5: rank by kinact/Ki surrogate (compute reactive Michael acceptor reactivity)

Reactivity Surrogates (computed without experiment)

For ranking warheads without wet-lab data:

Goal: Approximate Michael-acceptor reactivity from 2D structure without running DFT.

Approach: Parse the SMILES, locate the acrylamide substructure, count alpha-carbon substituents (more substitution lowers LUMO and slows reactivity), and return a negative count as a relative reactivity proxy.

def acceptor_lumo_surrogate(smi):
    # Crude: count alpha-substituents to acrylamide; more substituents = lower reactivity
    mol = Chem.MolFromSmiles(smi)
    if mol is None:
        return None
    acryl_pat = Chem.MolFromSmarts('[CX3](=O)[CX3]=[CX3]')
    matches = mol.GetSubstructMatches(acryl_pat)
    if not matches:
        return None
    # Count substituents on alpha-C
    alpha_c = mol.GetAtomWithIdx(matches[0][1])
    n_subs = len([n for n in alpha_c.GetNeighbors() if n.GetIdx() not in matches[0]])
    return -n_subs  # crude proxy (lower = more reactive)

For real reactivity prediction, DFT calculations (LUMO energy, HOMO-LUMO gap) are needed.

Per-Tool Failure Modes

Wrong warhead for residue

Trigger: Acrylamide warhead targeted at Tyr.

Mechanism: Acrylamides are Cys-selective; do not form bonds with Tyr at physiological pH.

Symptom: No covalent adduct observed despite docking pose.

Fix: Match warhead to residue: acrylamide/chloroacetamide for Cys; sulfonyl fluoride for Lys/Tyr/Ser.

Excessive reactivity (off-target)

Trigger: Chloroacetamide in drug-candidate context.

Mechanism: Too reactive; forms adducts with off-target Cys (e.g., GSH t1/2 < 30 min).

Symptom: Toxicity in cell-based assays; non-specific binding signal.

Fix: Replace with acrylamide (more selective); add alpha-substitution to acrylamide for tunable reactivity.

Geometric mismatch

Trigger: Warhead positioned but Cys-Cβ distance > 6 Å.

Mechanism: Even with reactive warhead, geometric reach matters; Cys side chain has limited reach.

Symptom: No covalent labeling in mass spec despite predicted docking.

Fix: Validate by measuring distance from warhead to Cys-Cβ; redock with constrained covalent bond.

Reversibility unintended

Trigger: Designed irreversible TCI but warhead is reversible.

Mechanism: Nitrile, aldehyde, boronate are reversible; equilibrium with non-covalent.

Symptom: Activity wanes after substrate washout in cellular assays.

Fix: Use truly irreversible warhead (acrylamide, chloroacetamide); or design for reversible covalent intentionally.

kinact/Ki conflation

Trigger: Optimizing for IC50 instead of kinact/Ki.

Mechanism: Compounds with same IC50 differ in covalent efficiency.

Symptom: Apparently identical compounds have different in vivo PK.

Fix: Always measure kinact/Ki (kinetic assay); rank by covalent efficiency.

DOCKovalent over-prediction

Trigger: Default DOCKovalent run.

Mechanism: Covalent constraint forces docking; many ligands "succeed" but are unrealistic.

Symptom: Many compounds pass docking; few label in vitro.

Fix: Post-filter by reactivity (chemoinformatics), geometric fit (Cys-Cβ distance), and PoseBusters.

Reconciliation: Irreversible vs Reversible Covalent

Common Errors

References

• Lonsdale & Ward, Chem. Soc. Rev. 47:3816 (2018) -- covalent revolution review.
• Singh et al., Nat. Rev. Drug Discov. 10:307 (2011) -- TCI design principles.
• London et al., Nat Chem Biol 10:1066-1072 (2014) -- DOCKovalent.
• Wu Q, Huang S-Y et al., Briefings Bioinform. 24:bbac559 (2023) -- HCovDock.
• Cai et al., J. Cheminformatics 14:39 (2022) -- GOLD covalent toolkit.
• Backus et al., Nat. Chem. 8:530 (2016) -- proteome-wide covalent ABPP.
• Pettinger et al., Angew. Chem. Int. Ed. 56:15200 (2017) -- reactive warhead reactivity quantification.
• Schwartz et al., Nat. Chem. Biol. 10:1006 (2014) -- KRAS G12C disulfide-tethered fragments.

Related Skills

• chemoinformatics/molecular-io - Parse warhead SMILES
• chemoinformatics/substructure-search - Warhead SMARTS detection
• chemoinformatics/virtual-screening - Pre-dock candidate non-covalent fit
• chemoinformatics/pose-validation - Validate covalent docking
• chemoinformatics/conformer-generation - Warhead conformer ensembles
• chemoinformatics/admet-prediction - ADMET of covalent leads
• chemoinformatics/molecular-descriptors - Reactivity surrogate descriptors