Unraveling the stereoisomeric toxicity of V-series nerve agent VX on human acetylcholinesterase: a well-tempered metadynamics study

Abstract

Highly toxic organophosphorus nerve agents exert their severe toxicity by inhibiting AChE irreversibly through the phosphonylation of Ser203 within the catalytic triad. The inhibition of AChE varies with stereoisomeric nerve agents. The inhibition mechanism of AChE by stereoselective nerve agents remains elusive. This study explores the mechanism of inhibition of hAChE by stereoselective nerve agents using well-tempered metadynamics (WT-MtD) simulations. The stereoisomeric nerve agent VX (in P_S and P_R forms) has been examined because it is known from the literature that the P_S form of VX inhibits human AChE faster than the P_R form. The 1D free energy surface profiles from the WT-MtD simulations reveal that the P_S_VX nerve agent binds more strongly (∼4.0 kcal mol⁻¹) with hAChE than the P_R_VX nerve agent. The WT-MtD simulation results showed that the catalytic triad residues Ser203, His447, and Glu334 play an important role in governing the inhibition process of hAChE by stereoisomeric VX nerve agents. The catalytic triad Ser203, His447, and Glu334 of hAChE aligns in a very similar fashion to that observed for the hydrolysis of acetylcholine (ACh) by P_S_VX. The aligned triad residues assist in generating the nucleophilic center at Ser203 for attack on the electrophilic phosphorus center of P_S_VX more effectively than P_R_VX. The WT-MtD snapshots taken from basin I for both stereoisomers of the nerve agents, used for superimposition onto the apo crystal structure of hAChE, show that the catalytic triad of P_R_VX bound hAChE is not aligned with the catalytic triad of the crystal structure. The role of water molecules in the slower inhibition of P_R_VX has also been investigated.