Catalytic Reaction Mechanism of Acetylcholinesterase Determined by Born−Oppenheimer Ab Initio QM/MM Molecular Dynamics Simulations
Authors:
- 15 June 2010
- journal article
- research article
- Published by American Chemical Society (ACS) in The Journal of Physical Chemistry B
- Vol. 114 (26) , 8817-8825
Abstract
Acetylcholinesterase (AChE) is a remarkably efficient serine hydrolase responsible for the termination of impulse signaling at cholinergic synapses. By employing Born−Oppenheimer molecular dynamics simulations with a B3LYP/6-31G(d) QM/MM potential and the umbrella sampling method, we have characterized its complete catalytic reaction mechanism for hydrolyzing neurotransmitter acetylcholine (ACh) and determined its multistep free-energy reaction profiles for the first time. In both acylation and deacylation reaction stages, the first step involves the nucleophilic attack on the carbonyl carbon, with the triad His447 serving as the general base, and leads to a tetrahedral covalent intermediate stabilized by the oxyanion hole. From the intermediate to the product, the orientation of the His447 ring needs to be adjusted very slightly, and then, the proton transfers from His447 to the product, and the break of the scissile bond happens spontaneously. For the three-pronged oxyanion hole, it only makes two hydrogen bonds with the carbonyl oxygen at either the initial reactant or the final product state, but the third hydrogen bond is formed and stable at all transition and intermediate states during the catalytic process. At the intermediate state of the acylation reaction, a short and low-barrier hydrogen bond (LBHB) is found to be formed between two catalytic triad residues His447 and Glu334, and the spontaneous proton transfer between two residues has been observed. However, it is only about 1−2 kcal/mol stronger than the normal hydrogen bond. In comparison with previous theoretical investigations of the AChE catalytic mechanism, our current study clearly demonstrates the power and advantages of employing Born−Oppenheimer ab initio QM/MM MD simulations in characterizing enzyme reaction mechanisms.Keywords
This publication has 97 references indexed in Scilit:
- Flexibility of Catalytic Zinc Coordination in Thermolysin and HDAC8: A Born−Oppenheimer ab Initio QM/MM Molecular Dynamics StudyJournal of Chemical Theory and Computation, 2009
- Born−Oppenheimer ab Initio QM/MM Molecular Dynamics Simulations of the Hydrolysis Reaction Catalyzed by Protein Arginine Deiminase 4The Journal of Physical Chemistry B, 2009
- Active Site Cysteine Is Protonated in the PAD4 Michaelis Complex: Evidence from Born−Oppenheimer Ab Initio QM/MM Molecular Dynamics SimulationsThe Journal of Physical Chemistry B, 2009
- Increasing the time step with mass scaling in Born‐Oppenheimer ab initio QM/MM molecular dynamics simulationsJournal of Computational Chemistry, 2009
- Toward Accurate Screening in Computer-Aided Enzyme DesignBiochemistry, 2009
- Highly Dissociative and Concerted Mechanism for the Nicotinamide Cleavage Reaction in Sir2Tm Enzyme Suggested by Ab Initio QM/MM Molecular Dynamics SimulationsJournal of the American Chemical Society, 2008
- Most Efficient Cocaine Hydrolase Designed by Virtual Screening of Transition StatesJournal of the American Chemical Society, 2008
- How Do SET-Domain Protein Lysine Methyltransferases Achieve the Methylation State Specificity? Revisited by Ab Initio QM/MM Molecular Dynamics SimulationsJournal of the American Chemical Society, 2008
- Ab Initio Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulation of Enzyme Catalysis: The Case of Histone Lysine Methyltransferase SET7/9The Journal of Physical Chemistry B, 2007
- Comparison of multiple Amber force fields and development of improved protein backbone parametersProteins: Structure, Function, and Bioinformatics, 2006