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

Residue	% of covalent drugs	Reactivity	Notes
Cysteine	~98%	High (nucleophile thiol)	Most accessible; preferred
Lysine	~1%	Moderate (amine)	Less reactive; selective for sulfonyl fluoride
Serine	<1%	Low (alcohol, requires activation)	β-lactam, boronate
Threonine	very rare	Low	Boronate, aldehyde
Tyrosine	very rare	Moderate (phenol)	Sulfonyl fluoride, fluorosulfate
Aspartate/Glutamate	very rare	Low (carboxylate)	Aldehyde Schiff base

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

Warhead	SMARTS pattern	Reactivity	Reversibility	Cys-selective
Acrylamide	C(=O)C=C	Moderate (Michael acceptor)	Irreversible	Yes
Chloroacetamide	C(=O)CCl	High (SN2)	Irreversible	Yes
α-haloketone	[CX3](=O)C[F,Cl,Br]	Very high	Irreversible	Yes (but reactive)
Vinyl sulfone	S(=O)(=O)C=C	Moderate (Michael)	Irreversible	Yes
Sulfonyl fluoride	S(=O)(=O)F	Moderate	Irreversible	Lys/Tyr/Ser
Fluorosulfate (SuFEx)	OS(=O)(=O)F	Moderate	Irreversible	Tyr/Lys
Aldehyde	C(=O)[H]	Variable	Reversible (covalent equilibrium)	Cys/Lys/Ser
Boronate (B-OH or B(OH)2)	B(O)O	Moderate	Reversible	Ser/Thr
Nitrile	C#N	Low	Reversible (Cys-S adduct)	Cys
Epoxide	C1OC1	High	Irreversible	Cys/Lys/Asp
α,β-unsaturated ketone	[CX3](=O)C=C	Moderate (Michael)	Irreversible	Cys
Isothiocyanate	N=C=S	High	Irreversible	Cys/Lys
Maleimide	C(=O)N(C(=O))C=C	Very high	Irreversible	Cys
Cysteine-selective heterocycle	various	Moderate	Variable	Yes (designed)

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

Goal	Warhead choice	Reactivity tier
Cysteine TCI, drug candidate	Acrylamide	Moderate (~kinact/Ki ~10^3-10^5 M^-1 s^-1)
Cysteine probe (chemical biology)	Chloroacetamide	High (kinact/Ki ~10^4-10^6 M^-1 s^-1)
Lysine TCI (uncommon)	Sulfonyl fluoride	Moderate
Tyrosine TCI	Fluorosulfate (SuFEx)	Moderate
Reversible covalent (KRAS G12C-like)	Acrylamide with α-substitution	Moderate reversibility
Activity-based protein profiling (ABPP)	Iodoacetamide / chloroacetamide	Very high
Boronic acid inhibitor (proteasome)	Boronate	Reversible
Aldehyde inhibitor (calpain)	Aldehyde	Reversible covalent

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

kinact/Ki range (M^-1 s^-1)	Interpretation	Reference
> 10^5	Highly efficient covalent inhibitor	Chloroacetamide probes, fragment-warhead TCIs
10^3 - 10^5	Standard for TCI; clinical candidate	Sotorasib (AMG510) KRAS G12C ~2x10^4 (Hallin 2020); adagrasib ~5x10^3
10^2 - 10^3	Moderate; clinical possible with high target dwell time	Ibrutinib BTK ~5x10^3 (Pan 2007)
< 10^2	Weak; needs warhead optimization	Reversibility likely dominates
<= 10	Probably not covalent (or wrong residue)	Background rate vs GSH

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

Warhead	GSH t1/2 at 10 mM	Risk
Chloroacetamide	minutes	High (reacts with off-target Cys)
Acrylamide	hours	Moderate
Substituted acrylamide (alpha-Me)	days	Low (drug-like)
Nitrile	days	Low
Sulfonyl fluoride	hours-days	Variable

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

Covalent Docking Tools

Tool	Approach	Strength	Fails when
DOCKovalent (London et al 2014 Nat Chem Biol 10:1066)	Constraint-based DOCK	Free, well-validated	Browser-based; small library
GOLD covalent (CCDC)	GOLD with covalent constraint	Commercial; selectivity	License cost
AutoDock 4 covalent	AD4 with covalent bond	Open source	Slower than Vina
CovDock (Schrödinger)	Glide-based + covalent	Commercial best	License cost
MOE covalent	Triposite Discovery	Commercial	License cost
HCovDock (Wu Q, Huang S-Y 2023 Briefings Bioinform 24:bbac559)	Hierarchical fragment + covalent	Open; supports many residues	Newer, less validated
ICM-Pro covalent	Active site grid + covalent	Commercial; metal centers	License cost

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:

Descriptor	Use case
LUMO energy (DFT)	Michael acceptor reactivity (lower LUMO = more reactive)
Electrophile partial charge	SN2 reactivity
RDKit rdMolDescriptors.CalcLabuteASA	Steric accessibility
AlphaFold3 / Boltz-2 binding pose	Geometric fit to reactive Cys

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

Aspect	Irreversible	Reversible covalent
Examples	KRAS G12C (acrylamide), BTK (ibrutinib)	Boronate (bortezomib), aldehyde (calpain inhibitors)
Toxicity profile	Off-target Cys labeling potential	Off-target equilibrium
Resistance mechanism	Mutation of reactive Cys	Mutation reduces affinity
Patent / IP	Stronger (specific bond)	Standard
When to choose	If Cys is hot-spot, conserved, druggable	If reversibility critical (e.g., proteasome)

Common Errors

Symptom	Cause	Fix
Warhead not matching SMARTS	Different stereochemistry or charged	Use canonicalized + neutral SMARTS
DOCKovalent rejects ligand	No suitable Cys in pocket	Re-check residue accessibility
GSH adduct dominates	Warhead too reactive	Use less reactive warhead; or alpha-substitute
Off-target labeling in cells	Promiscuous warhead	Iterate warhead reactivity vs selectivity
Docking pose but no labeling	Geometric mismatch	Distance check; rotamer search
Reversible inhibitor not acting irreversibly	Wrong warhead class	Re-check reaction mechanism
HCovDock fails on PROTAC	Tool optimized for monomer covalent	Use specialized tools for bivalent

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