Characterization of a complete cycle of acetylcholinesterase catalysis by ab initio QM/MM modeling

Selective modulation of alkali metal ions on acetylcholinesterase

Acetylcholinesterase (AChE) is an important hydrolase in cholinergic synapses and a candidate target in the treatment of Alzheimer's disease. The lithium treatment widely used in neurological disorders can alter the AChE activity, yet the underlying mechanism of how the ion species regulate the enzymatic activity remains unclear. In this work, we performed combined quantum mechanics/molecular mechanics (QM/MM) and molecular dynamics (MD) simulations and well-tempered metadynamics to understand the modulation of human AChE (hAChE) activity using three alkali metal ions (Li+, Na+, and K+) in different concentrations. Our simulations show that the binding affinity and catalytic activity are affected by different ion species through allosteric ion coordination geometries on the hAChE complex and distant electrostatic screening effect. A Li+ cluster involving D330, E393, and D397 residues and three Li+ ions was found to be highly conserved and can be critical to the enzyme activity. Binding energy calculations indicate that the electrostatic screening from allosterically bound cations can affect the key residues at the catalytic site and active-site gorge, including E199. Furthermore, an increase in ion concentration can lead to lower reactivity, especially for Li+ ions, which exhibit more cation-hAChE contacts than Na+ and K+. The selective ion binding and their preferred modulation on hAChE are highly related to ion species. This work provides a molecular perspective on selective modulation by different ion species of the enzyme catalytic processes.

Novel Uracil-Based Inhibitors of Acetylcholinesterase with Potency for Treating Memory Impairment in an Animal Model of Alzheimer's Disease

Novel derivatives based on 6-methyluracil and condensed uracil, 2,4-quinazoline-2,4-dione, were synthesized with terminal meta- and para-benzoate moieties in polymethylene chains at the N atoms of the pyrimidine ring. In the synthesized compounds, the polymethylene chains were varied from having tris- to hexamethylene chains and quaternary ammonium groups; varying substituents (ester, salt, acid) at benzene ring were introduced into the chains and benzoate moieties. In vivo biological experiments demonstrated the potency of these compounds in decreasing the number of β-amyloid plaques and their suitability for the treatment of memory impairment in a transgenic model of Alzheimer's disease.

Interaction between glyphosate pesticide and amphiphilic peptides for colorimetric analysis

The large-scale use of glyphosate pesticides in food production has attracted attention due to environmental damage and toxicity risks. Several regulatory authorities have established safe limits or concentrations of these pesticides in water and various food products consumed daily. The irreversible inhibition of acetylcholinesterase (AChE) activity is one of the strategies used for pesticide detection. Herein, we found that lipopeptide sequences can act as biomimetic microenvironments of AChE, showing higher catalytic activities than natural enzymes in an aqueous solution, based on IC₅₀ values. These biomolecules contain in the hydrophilic part the amino acids l-proline (P), l-arginine (R), l-tryptophan (W), and l-glycine (G), covalently linked to a hydrophobic part formed by one or two long aliphatic chains. The obtained materials are referred to as compounds 1 and 2, respectively. According to fluorescence assays, 2 is more hydrophobic than 1. The circular dichroism (CD) data present a significant difference in the molar ellipticity values, likely related to distinct conformations assumed by the proline residue in the lipopeptide supramolecular structure in solution. The morphological aspect was further characterized using small-angle X-ray scattering (SAXS) and cryogenic transmission electron microscopy (cryo-TEM), which showed that compounds 1 and 2 self-assembly into cylindrical and planar core-shell structures, respectively. The mimetic AchE behaviour of lipopeptides was confirmed by Ellman's hydrolysis reaction, where the proline residue in the peptides act as a nucleophilic scavenger of organophosphate pesticides. Moreover, the isothermal titration calorimetry (ITC) experiments revealed that host-guest interactions in both systems were dominated by enthalpically-driven thermodynamics. UV-vis kinetic experiments were performed to assess the inhibition of the lipopeptide catalytic activity and the IC₅₀ values were obtained, and we found that the detection limit correlated with the increase in hydrophobicity of the lipopeptides, implying the micellization process is more favorable.

Mechanism and biomass association of glucuronoyl esterase: an α/β hydrolase with potential in biomass conversion

Glucuronoyl esterases (GEs) are α/β serine hydrolases and a relatively new addition in the toolbox to reduce the recalcitrance of lignocellulose, the biggest obstacle in cost-effective utilization of this important renewable resource. While biochemical and structural characterization of GEs have progressed greatly recently, there have yet been no mechanistic studies shedding light onto the rate-limiting steps relevant for biomass conversion. The bacterial GE OtCE15A possesses a classical yet distinctive catalytic machinery, with easily identifiable catalytic Ser/His completed by two acidic residues (Glu and Asp) rather than one as in the classical triad, and an Arg side chain participating in the oxyanion hole. By QM/MM calculations, we identified deacylation as the decisive step in catalysis, and quantified the role of Asp, Glu and Arg, showing the latter to be particularly important. The results agree well with experimental and structural data. We further calculated the free-energy barrier of post-catalysis dissociation from a complex natural substrate, suggesting that in industrial settings non-catalytic processes may constitute the rate-limiting step, and pointing to future directions for enzyme engineering in biomass utilization.

Zong and coworkers combine computational and experimental methods to decipher in detail the mechanism of action of glucuronoyl esterases, enzymes with significant biotechnological potential for decoupling lignin from polysaccharides in biomass.

Chemical profile of Lippia thymoides, evaluation of the acetylcholinesterase inhibitory activity of its essential oil, and molecular docking and molecular dynamics simulations

The essential oils of the fresh and dry flowers, leaves, branches, and roots of Lippia thymoides were obtained by hydrodistillation and analyzed using gas chromatography (GC) and GC-mass spectrometry (MS). The acetylcholinesterase inhibitory activity of the essential oil of fresh leaves was investigated on silica gel plates. The interactions of the key compounds with acetylcholinesterase were simulated by molecular docking and molecular dynamics studies. In total, 75 compounds were identified, and oxygenated monoterpenes were the dominant components of all the plant parts, ranging from 19.48% to 84.99%. In the roots, the main compounds were saturated and unsaturated fatty acids, having contents varying from 39.5% to 32.17%, respectively. In the evaluation of the anticholinesterase activity, the essential oils (detection limit (DL) = 0.1 ng/spot) were found to be about ten times less active than that of physostigmine (DL = 0.01ng/spot), whereas thymol and thymol acetate presented DL values each of 0.01 ng/spot, equivalent to that of the positive control. Based on the docking and molecular dynamics studies, thymol and thymol acetate interact with the catalytic residues Ser203 and His447 of the active site of acetylcholinesterase. The binding free energies (ΔGbind) for these ligands were -18.49 and -26.88 kcal/mol, demonstrating that the ligands are able to interact with the protein and inhibit their catalytic activity.

On the Influence of the Protonation States of Active Site Residues on AChE Reactivation: A QM/MM Approach

Acetylcholinesterase (AChE), an enzyme of the serine hydrolase superfamily, is a mediator of signal transmission at cholinergic synapses by catalyzing acetylcholine cleavage into acetate and choline. This enzyme is vulnerable to covalent inhibition by organophosphate compounds (like VX). Covalent inhibition of AChE does not revert spontaneously. Known reactivator compounds have limited action in restoring catalytic activity. QM/MM simulations of VX-inhibited AChE reactivation by pralidoxime (2-PAM), a classical reactivator, were performed. These afforded a broad view of the effect of protonation states of active-site residues, and provide evidence for the role of Glu202, which needs to be protonated for reactivation to occur. In situ deprotonation of 2-PAM for both protonation states of Glu202 showed that His447 is able to deprotonate 2-PAM with the assistance of Glu202. Because the active site of serine hydrolases is highly conserved, this work provides new insights on the interplay between the catalytic triad residues and this glutamate, newly identified as protonatable.

Modeling reactivation of the phosphorylated human butyrylcholinesterase by QM(DFTB)/MM calculations

Human butyrylcholinesterase (BChE) is a bioscavenger that protects the enzyme which is critical for the central nerve system, acetylcholinesterase, from poisoning by organophosphorus agents. Elucidating the details of the hydrolysis reaction mechanism is important to understand how the phosphorylated BChE can be reactivated. Application of the QM(DFTB)/MM(AMBER) method to construct the minimum energy pathways for the hydrolysis reaction of the diethylphosphorylated BChE allowed us to suggest a mechanism of reactivation of the wild-type and the G117H mutated enzyme. Unlike previous approaches assuming that either His438 or His117 serves as a general base in the catalysis, in our proposal the Glu197 residue is responsible for activation of the nucleophilic water molecule (Wat) leading to the chemical transformations that restore the catalytic Ser198 residue in BChE. In agreement with the experimental data, it is shown that the G117H mutation facilitates the reactivation of the inhibited enzyme.

How is acetylcholinesterase phosphonylated by soman? An ab initio QM/MM molecular dynamics study

Acetylcholinesterase (AChE) is a crucial enzyme in the cholinergic nerve system that hydrolyzes acetylcholine (ACh) and terminates synaptic signals by reducing the effective concentration of ACh in the synaptic clefts. Organophosphate compounds irreversibly inhibit AChEs, leading to irreparable damage to nerve cells. By employing Born–Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, a state-of-the-art approach to simulate enzyme reactions, we have characterized the covalent inhibition mechanism between AChE and the nerve toxin soman and determined its free energy profile for the first time. Our results indicate that phosphonylation of the catalytic serine by soman employs an addition–elimination mechanism, which is highly associative and stepwise: in the initial addition step, which is also rate-limiting, His440 acts as a general base to facilitate the nucleophilic attack of Ser200 on the soman’s phosphorus atom to form a trigonal bipyrimidal pentacovalent intermediate; in the subsequent elimination step, Try121 of the catalytic gorge stabilizes the leaving fluorine atom prior to its dissociation from the active site. Together with our previous characterization of the aging mechanism of soman inhibited AChE, our simulations have revealed detailed molecular mechanistic insights into the damaging function of the nerve agent soman.

Molecular modeling evidence for His438 flip in the mechanism of butyrylcholinesterase hysteretic behavior

Cholinesterases display a hysteretic behavior with certain substrates and irreversible inhibitors. For years, this behavior has remained puzzling. However, several lines of evidence indicated that it is caused by perturbation of the catalytic triad and its water environment. In the present study, using molecular dynamics simulations of Ala328Cys BuChE mutant and wild-type BuChE in the absence and presence of a co-solvent (sucrose, glycerol), we provide evidence that hysteresis originates in a flip of the catalytic triad histidine (His438). This event is controlled by water molecules that interact with active site residues. The physiological significance of this phenomenon is still an issue.

Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology

Computational enzymology is a rapidly maturing field that is increasingly integral to understanding mechanisms of enzyme-catalyzed reactions and their practical applications. Combined quantum mechanics/molecular mechanics (QM/MM) methods are important in this field. By treating the reacting species with a quantum mechanical method (i.e., a method that calculates the electronic structure of the active site) and including the enzyme environment with simpler molecular mechanical methods, enzyme reactions can be modeled. Here, we review QM/MM methods and their application to enzyme-catalyzed reactions to investigate fundamental and practical problems in enzymology. A range of QM/MM methods is available, from cheaper and more approximate methods, which can be used for molecular dynamics simulations, to highly accurate electronic structure methods. We discuss how modeling of reactions using such methods can provide detailed insight into enzyme mechanisms and illustrate this by reviewing some recent applications. We outline some practical considerations for such simulations. Further, we highlight applications that show how QM/MM methods can contribute to the practical development and application of enzymology, e.g., in the interpretation and prediction of the effects of mutagenesis and in drug and catalyst design.

Reaction pathway and free energy profiles for butyrylcholinesterase-catalyzed hydrolysis of acetylthiocholine

The catalytic mechanism for butyrylcholineserase (BChE)-catalyzed hydrolysis of acetylthiocholine (ATCh) has been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical-free energy (QM/MM-FE) calculations on both acylation and deacylation of BChE. Additional quantum mechanical (QM) calculations have been carried out, along with the QM/MM-FE calculations, to understand the known substrate activation effect on the enzymatic hydrolysis of ATCh. It has been shown that the acylation of BChE with ATCh consists of two reaction steps including the nucleophilic attack on the carbonyl carbon of ATCh and the dissociation of thiocholine ester. The deacylation stage includes nucleophilic attack of a water molecule on the carboxyl carbon of substrate and dissociation between the carboxyl carbon of substrate and hydroxyl oxygen of Ser198 side chain. QM/MM-FE calculation results reveal that the acylation of BChE is rate-determining. It has also been demonstrated that an additional substrate molecule binding to the peripheral anionic site (PAS) of BChE is responsible for the substrate activation effect. In the presence of this additional substrate molecule at PAS, the calculated free energy barrier for the acylation stage (rate-determining step) is decreased by ~1.7 kcal/mol. All of our computational predictions are consistent with available experimental kinetic data. The overall free energy barriers calculated for BChE-catalyzed hydrolysis of ATCh at regular hydrolysis phase and substrate activation phase are ~13.6 and ~11.9 kcal/mol, respectively, which are in reasonable agreement with the corresponding experimentally derived activation free energies of 14.0 kcal/mol (for regular hydrolysis phase) and 13.5 kcal/mol (for substrate activation phase).

Reaction pathway and free energy profile for butyrylcholinesterase-catalyzed hydrolysis of acetylcholine

A catalytic mechanism for the butyrylcholinesterase (BChE)-catalyzed hydrolysis of acetylcholine (ACh) has been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical-free energy calculations on both acylation and deacylation of BChE. It has been shown that the acylation with ACh includes two reaction steps, including nucleophilic attack on the carbonyl carbon of ACh and dissociation of choline ester. The deacylation stage includes nucleophilic attack of a water molecule on the carboxyl carbon of the substrate and dissociation between the carboxyl carbon of the substrate and the hydroxyl oxygen of the Ser198 side chain. Notably, despite the fact that acetylcholinesterase (AChE) and BChE are very similar enzymes, the acylation of BChE with ACh is rate-determining, which is remarkably different from the AChE-catalyzed hydrolysis of ACh, in which the deacylation is rate-determining. The computational prediction is consistent with available experimental kinetic data. The overall free energy barrier calculated for BChE-catalyzed hydrolysis of ACh is 13.8 kcal/mol, which is in good agreement with the experimentally derived activation free energy of 13.3 kcal/mol.

Catalytic reaction mechanism of acetylcholinesterase determined by Born-Oppenheimer ab initio QM/MM molecular dynamics simulations

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.

Fixation of the two Tabun isomers in acetylcholinesterase: a QM/MM study

Dysfunction of acetylcholinesterase (AChE) due to inhibition by organophosphorus (OP) compounds is a major threat since AChE is a key enzyme in neurotransmission. To more rigorously design reactivation agents, it is of prime importance to understand the mechanism of inhibition of AChE by OP compounds. Tabun is one of the more potent nerve agents. It is produced as a mixture of two enantiomers, one of them (the levorotatory isomer) being 6.3 times more potent. Could it be that the inhibition mechanism is different for the two enantiomers? To address this critical issue, we used a hybrid quantum mechanics/molecular mechanics (QM/MM) methodology. Calculations were performed using BP86 functional and TZVP basis set. Single points were also done with B3LYP and PBE0 functionals. We studied the four possible attacks of tabun on the oxygen of Ser203 using two crystallographic structures (PDB codes 2C0P and 3DL7): (S) tabun with the cyano group syn to the oxygen of Ser203 and (R) tabun with the cyano group anti, corresponding to the experimental X-ray structure; (S) tabun with the cyano group anti to the oxygen of Ser203 and (R) tabun with the cyano group syn, leading to a different isomer than was experimentally seen. We found that the most active enantiomer is (S) tabun with the cyano group syn to the oxygen of Ser203. Thus it seems that the cyano group does not leave anti to the oxygen of Ser203 due to repulsive polar interactions between cyanide and aromatic residues in the active site.