Abstract
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.
Supplementary materials
Title
Supplementary Informations
Description
Various supplementary figures and computational details
Actions