Abstract:Acetylcholinesterase is one of the most significant known serine hydrolases that governs 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 the different 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 molecular dynamics simulations for ligand diffusion and hybrid quantum mechanics/molecular mechanics 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 the two steps were 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 into a strong polarization effect acting on water molecules near the active site. An AMOEBA polarizable 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 reorganization and the polarization of confined water molecules. Consequently, carefully considering and modeling buried (polarizable) water molecules are of paramount importance when modeling full enzymatic activity. Therefore, this work will also provide valuable insights for future research on related enzymes with buried active sites.

Correlating Predicted Reactivities with Experimental Inhibition Data of Covalent ChlaDUB1 Inhibitors

Activity prediction of substrates in NADH-dependent carbonyl reductase by docking requires catalytic constraints and charge parameterization of catalytic zinc environment

Discovery and Computational Studies of Potent Covalent Kinase Inhibitors with -Substituent Electrophiles Targeting Cysteine

Discovery and Computational Studies of Potent Covalent Kinase Inhibitors with alpha-Substituent Electrophiles Targeting Cysteine

Correlation of Experimental and Calculated Inhibition Constants of Protease Inhibitor Complexes

Machine Learning Identifies Chemical Characteristics That Promote Enzyme Catalysis

Reactivities of acrylamide warheads toward cysteine targets: a QM/ML approach to covalent inhibitor design

Machine Learning Prediction of Hydrogen Atom Transfer Reactivity in Photoredox-Mediated C–H Functionalization

Chemical Reactivity Prediction: Current Methods and Different Application Areas

ReactZyme: A Benchmark for Enzyme-Reaction Prediction

Using Protein-Confined Proximity To Determine Chemical Reactivity

Control of catalytic efficiency by a co-evolving network of catalytic and non-catalytic residues

Application of AlphaFold models in evaluating ligandable cysteines across E3 ligases

Rapid and Reliable Binding Affinity Prediction of Bromodomain Inhibitors: A Computational Study

Exploring PROTAC Cooperativity with Coarse-Grained Alchemical Methods

Combining de novo molecular design with semiempirical protein-ligand binding free energy calculation

Global Reactivity Profiling of the Catalytic Lysine in Human Kinome for Covalent Inhibitor Development

Statistical Analysis and Prediction of Covalent Ligand Targeted Cysteine Residues

Computational Study of Protein-Ligand Unbinding for Enzyme Engineering

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

Structure-Based Design and Synthesis of Covalent Inhibitors for Deubiquitinase and Acetyltransferase ChlaDUB1 of Chlamydia trachomatis