Benchmarking Affinity-Based Docking for Irreversible Acetylcholinesterase Inhibitors: A Physics-Informed Reactivity...

Irreversible organophosphorus inhibitors of acetylcholinesterase (AChE) act by covalent phosphylation of the catalytic serine, yet their computational prioritization is often based on affinity scores from noncovalent docking. We benchmark that assumption against a curated set of measured human-AChE inhibition rate constants (n = 14; G-, V-, and organophosphate-pesticide classes, including a configuration-resolved VX pair). AutoDock Vina affinity (Spearman ρ = +0.12) and GNINA CNN scores (ρ = -0.40) do not provide statistically reliable potency rankings; explicit covalent docking also remains weak (ρ = +0.23, n = 11), with bootstrap confidence intervals crossing zero. Selecting covalent-compatible (near-attack) poses rather than top-ranked poses does not change this outcome. A quantum-chemical cluster scan supports a mechanistic interpretation: once the activated serine alkoxide/reactive complex is formed, P-O bond formation is downhill within the model, so the discriminating information is expected to reside in electrophilicity, leaving-group chemistry, and transport into reactive geometry rather than in a single docked pose. Consistently, mechanism-oriented molecular descriptors recover a modest but reproducible signal: TPSA alone correlates with potency (ρ = -0.62), and a Gaussian-process model used primarily for calibrated uncertainty yields leave-one-out ρ = +0.55 with permutation support (p ≈ 0.01) and a bootstrap CI whose lower bound lies near zero. A leave-both-VX-enantiomers-out test preserves the experimental S _P > R _P ordering. Given the small panel (n = 14), the descriptor correlations and stereochemical results are reported as exploratory and uncertainty-bounded rather than predictive. The resulting framework is presented as an uncertainty-bounded prioritization tool for data-scarce covalent AChE inhibitors, not as a substitute for kinetic measurement.