Unveiling the Full Dynamical and Reactivity Profiles of Acetylcholinesterase: A Comprehensive All-Atom Investigation

17 November 2023, Version 1
This content is a preprint and has not undergone peer review at the time of posting.


Acetylcholinesterase is one of the most significant known serine hydrolases, govern- ing the mammalian nervous system. Its high rate speed, operating at the diffusion limit, combined with its buried active site feature, has made it a subject of extensive research over the last decades. Despite several studies focused on atomistic details of its differ- ent steps, a comprehensive theoretical investigation of the entire catalytic cycle has not yet been reported. In this work, we present an intuitive workflow aiming at describing the full dynamical and reactive profiles of AChE by coupling extensive Steered Molec- ular Dynamics (SMD) simulations for ligand diffusion and hybrid Quantum Mechan- ics/Molecular Mechanics (QM/MM) computations to decipher the complete reactivity of the substrate within the enzyme. This comprehensive approach provides a broader view of the interconnections between each step that would not be readily accessible if studied independently. Our simulations reveal that although individual steps do not indicate any strong limiting step, a solvent water molecule reorganization between the acylation and deacylation processes through the reactivity results in an energy cost of 20 kcal/mol. The observed barrier surpasses all others and discloses insights on a strong polarization effect acting on water molecules near the active site. An AMOEBA polar- izable Molecular Dynamics simulation tends to confirm this assumption by capturing a substantial dipole moment (3.10 D) on the water molecule closest to the reaction site. These results shed light on the crucial correlation between this high–energy water re- organization and the polarization of confined water molecules. Consequently, carefully considering and modeling buried (polarizable) water molecules is of paramount impor- tance when modeling full enzymatic activity. Therefore, this work will also provide valuable insights for future research on related enzymes with buried active sites.


molecular dynamics

Supplementary materials

Supplementary Informations
Various supplementary figures and computational details


Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.