Benchmarking Affinity-Based Docking for Irreversible Acetylcholinesterase Inhibitors: A Physics-Informed Reactivity Framework for Stereoselective Potency Prioritization
- 1Seoul National University , , ,
- 2Chemical, Biological, Radiological and Nuclear (CBRN) Defense Research Institute , ,
Abstract
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 SP > RP 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.
Related
- In Pursuit of Enduring Inactivation: A Perspective on the Development and Study of Irreversible-Covalent Inhibitors.
- Heteroaryl-Substituted Secondary Amines as Multifunctional Antimicrobial, Antibiofilm, Acetylcholinesterase, and α-Glucosidase Inhibitors
- Old Drugs, New Targets: Multi-parameter computational prioritization of FDA‑approved, brain‑penetrant drugs as isoform‑selective GSK‑3β inhibitors for Alzheimer’s disease
- Design and Benchmarking of CK2α′-Preferring Inhibitors for Huntington’s Disease: The Limits of Structure-Based Isoform Selectivity
- Acetylcholinesterase Activity Modulation by Proteinogenic Amino Acid Ester Ferulates salts with Enhanced Solubility and Antioxidant Properties
BibTeX
@article{kim2026benchmarking,
title = {{Benchmarking Affinity-Based Docking for Irreversible Acetylcholinesterase Inhibitors: A Physics-Informed Reactivity Framework for Stereoselective Potency Prioritization}},
author = {Jeongyun Kim and Jin Yoo and Doo-Hee Lee and Ku Kang},
journal = {ACS omega},
year = {2026},
month = sep,
doi = {10.1021/acsomega.6c08355},
url = {https://europepmc.org/article/PMC/PMC13670796},
note = {CC BY-NC-ND},
}
RIS
TY - JOUR
TI - Benchmarking Affinity-Based Docking for Irreversible Acetylcholinesterase Inhibitors: A Physics-Informed Reactivity Framework for Stereoselective Potency Prioritization
AU - Jeongyun Kim
AU - Jin Yoo
AU - Doo-Hee Lee
AU - Ku Kang
PY - 2026
DA - 2026/09/30
JO - ACS omega
DO - 10.1021/acsomega.6c08355
UR - https://europepmc.org/article/PMC/PMC13670796
AB - 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 SP > RP 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.
SN - 2470-1343
ER -