QM/MM investigation on the catalytic mechanism of Bacteroides thetaiotaomicron α-glucosidase BtGH97a

Bacteroides thetaiotaomicron α-glucosidase BtGH97a is an inverting enzyme. In this paper, the hydrolysis mechanism of p-nitro-phenyl α-d-glucopyranoside (pNP-Glc) catalyzed by BtGH97a was firstly studied by using quantum mechanical/molecular mechanical (QM/MM) approach. Two possible reaction pathways were considered. In the first pathway, a water molecule deprotonated by a nucleophilic base (here E439 or E508) attacks firstly on the anomeric carbon of pNP-Glc, then a proton from an acid residue (E532) attacks on the glycosidic oxygen to finish the hydrolysis reaction (named as nucleophilic attack-first pathway). In the second pathway, the proton from E532 attacks firstly on the glycosidic oxygen, then the water deprotonated by the nucleophilic base attacks on the anomeric carbon of pNP-Glc (named as proton attack-first pathway). Our calculation results indicate that the nucleophilic attack-first pathway is favorable in energy, in which the nucleophilic attack process is the rate-determining step with an energy barrier of 15.4kcal/mol in the case of residue E508 as nucleophilic base. In this rate-determining step, the deprotonation of water and the attack on the anomeric carbon are concerted. In the proton attack-first pathway, the proton attack on the glycosidic oxygen is the rate-determining step, and the energy barrier is 24.1kcal/mol. We conclude that the hydrolysis mechanism would follow nucleophilic attack-first pathway.     

QM/MM investigation on the catalytic mechanism of Bacteroides thetaiotaomicron α-glucosidase BtGH97a

Bacteroides thetaiotaomicron α-glucosidase BtGH97a is an inverting enzyme. In this paper, the hydrolysis mechanism of p-nitro-phenyl α-d-glucopyranoside (pNP-Glc) catalyzed by BtGH97a was firstly studied by using quantum mechanical/molecular mechanical (QM/MM) approach. Two possible reaction pathways were considered. In the first pathway, a water molecule deprotonated by a nucleophilic base (here E439 or E508) attacks firstly on the anomeric carbon of pNP-Glc, then a proton from an acid residue (E532) attacks on the glycosidic oxygen to finish the hydrolysis reaction (named as nucleophilic attack-first pathway). In the second pathway, the proton from E532 attacks firstly on the glycosidic oxygen, then the water deprotonated by the nucleophilic base attacks on the anomeric carbon of pNP-Glc (named as proton attack-first pathway). Our calculation results indicate that the nucleophilic attack-first pathway is favorable in energy, in which the nucleophilic attack process is the rate-determining step with an energy barrier of 15.4 kcal/mol in the case of residue E508 as nucleophilic base. In this rate-determining step, the deprotonation of water and the attack on the anomeric carbon are concerted. In the proton attack-first pathway, the proton attack on the glycosidic oxygen is the rate-determining step, and the energy barrier is 24.1 kcal/mol. We conclude that the hydrolysis mechanism would follow nucleophilic attack-first pathway.     

A computational study of the deacylation mechanism of human butyrylcholinesterase

To investigate the mechanism of the deacylation reaction in the active site of human butyrylcholinesterase (BuChE), we carried out quantum mechanical (QM) calculations on cluster models of the active site built from a crystallographic structure. The models consisted of the substrate butyrate moiety, the catalytic triad of residues (Ser198, Glu325, and His438), the "oxy-anion hole" (Gly116, Gly117, and Ala199), the side chain of Glu197, four water molecules, the side chain of Ser225, and the peptide linkage between Val321 and Asn322. Analyses of the equilibrium geometries, electronic properties, and energies of the QM models gave insights into the catalytic mechanism. In addition, the QM calculations provided the data required to build a molecular mechanics representation of the reactive BuChE region that was employed in molecular dynamics simulations followed by molecular-mechanics-Poisson-Boltzmann (MM-PB) calculations. Subsequently, we combined the QM energies with average MM-PB energies to estimate the free energy of the reactive structures in the enzyme. The rate-determining step corresponds to the formation of a tetrahedral intermediate with a free-energy barrier of approximately 14.0 kcal/mol. The modulation of the BuChE activity, exerted by either neutral molecules (glycerol, GOL) or a second butyrylcholine (CHO) molecule bound to the cation-pi site, does not involve any significant allosteric effect. Interestingly, the presence of GOL or CHO stabilizes a product complex formed between a butyric acid molecule and BuChE. These results are in consonance with the crystallographic structure of BuChE, in which the catalytic Ser198 interacts with a butyric fragment, while the cation-pi site is occupied by one GOL molecule.     

Mechanism of carboxypeptidase-Y-catalyzed reaction deduced from a pressure-dependence study

The activation volumes for kcat of the carboxypeptidase-Y-catalyzed hydrolysis of ester substrates were slightly negative (-1 to -4 ml/mol), while those for peptide and depsi-peptide analog were highly positive (+10 to +27 ml/mol). These values and the contrasting pH dependences of these two groups of the substrates are explained by a mechanism involving three ionic states of the enzyme and the second stable intermediate (acyl-enzyme). Esters are mostly rate-controlled by the deacylation step and peptides are controlled by both the acylation and the deacylation steps. Pressure increase induced a partial shift of the rate-determining step. Reaction volumes for Km-1 of peptide and depsi-peptide analog showed large and positive values (+16 to +29 ml/mol) which reflects the electrostatic interaction in the substrate recognition by this enzyme.     

QM/MM study of the active site of free papain and of the NMA-papain complex.

Hybrid quantum mechanical/molecular mechanical (QM/MM) calculations using restricted and unrestricted Hartree-Fock and B3LYP ab initio (QM) and Amber force field (MM), respectively, have been applied to study the catalytic site of papain in both free and substrate bonded forms. Ab initio geometry optimizations have been performed for the active site of papain and the N-methyl-acetamide (NMA)-papain complex within the molecular mechanical treatment of the protein environment. A covalent tetrahedral intermediate structure could be obtained only when the amide N atom of the substrate molecule was protonated through a proton transfer from the His-159 in the catalytic site. Our results support the previous assumption that a proton transfer from His-159 to the amide N atom of the substrate occurs prior to or concerted with the nucleophilic attack of the Cys-25 sulfur atom to the carbonyl group of the substrate. The electron correlation effect will reduce the proton transfer barrier. Therefore, this proton transfer can be easily observed in the B3LYP/6-31G* calculations. The HF/6-31G* method overestimates the reaction barrier against this proton transfer. The sulfur atom of Cys-25 and the imidazole ring of His-159 are found to be coplanar in the free form of the enzyme. However, the rotation of the imidazole ring of His-159 was observed during the formation of the tetrahedral intermediate. Without the papain environment, the coplanar thiolate-imidazolium ion pair RS-...ImH+ is much less stable than the neutral form of RSH....Im. Within the protein environment, however, the thiolate-imidazolium ion pair becomes more stable than its neutral form by 4.1 and 0.4 kcal/mol in HF/6-31G* and B3LYP/6-31G* calculations, respectively. The barrier of proton transfer from S-H group of Cys-25 to the imidazole ring of His-159 was reduced from 22.0 kcal/mol to 15.2 kcal/mol by the protein environment in HF/6-31G* calculations. This barrier is found to be much smaller (2.5 kcal/mol) in B3LYP/6-31G* calculations.     

Modeling evolution of hydrogen bonding and stabilization of transition states in the process of cocaine hydrolysis catalyzed by human butyrylcholinesterase

Molecular dynamics (MD) simulations and quantum mechanical/molecular mechanical (QM/MM) calculations were performed on the prereactive enzyme-substrate complex, transition states, intermediates, and product involved in the process of human butyrylcholinesterase (BChE)-catalyzed hydrolysis of (-)-cocaine. The computational results consistently reveal a unique role of the oxyanion hole (consisting of G116, G117, and A199) in BChE-catalyzed hydrolysis of cocaine, compared to acetylcholinesterase (AChE)-catalyzed hydrolysis of acetylcholine. During BChE-catalyzed hydrolysis of cocaine, only G117 has a hydrogen bond with the carbonyl oxygen (O31) of the cocaine benzoyl ester in the prereactive BChE-cocaine complex, and the NH groups of G117 and A199 are hydrogen-bonded with O31 of cocaine in all of the transition states and intermediates. Surprisingly, the NH hydrogen of G116 forms an unexpected hydrogen bond with the carboxyl group of E197 side chain and, therefore, is not available to form a hydrogen bond with O31 of cocaine in the acylation. The NH hydrogen of G116 is only partially available to form a weak hydrogen bond with O31 of cocaine in some structures involved in the deacylation. The change of the estimated hydrogen-bonding energy between the oxyanion hole and O31 of cocaine during the reaction process demonstrates how the protein environment can affect the energy barrier for each step of the BChE-catalyzed hydrolysis of cocaine. These insights concerning the effects of the oxyanion hole on the energy barriers provide valuable clues on how to rationally design BChE mutants with a higher catalytic activity for the hydrolysis of (-)-cocaine.     

Steady-state kinetic analysis of human cholinesterases over wide concentration ranges of competing substrates

Substrate competition for human acetylcholinesterase (AChE) and human butyrylcholinesterase (BChE) was studies under steady-state conditions using wide range of substrate concentrations. Competing couples of substates were acetyl-(thio)esters. Phenyl acetate (PhA) was the reporter substrate and competitor were either acetylcholine (ACh) or acetylthiocholine (ATC). The common point between investigated substrates is that the acyl moiety is acetate, i.e. same deacylation rate constant for reporter and competitor substrate.Steady-state kinetics of cholinesterase-catalyzed hydrolysis of PhA in the presence of ACh or ATC revealed 3 phases of inhibition as concentration of competitor increased: a) competitive inhibition, b) partially mixed inhibition, c) partially uncompetitive inhibition for AChE and partially uncompetitive activation for BChE. This sequence reflects binding of competitor in the active centrer at low concentration and on the peripheral anionic site (PAS) at high concentration. In particular, it showed that binding of a competing ligand on PAS may affect the catalytic behavior of AChE and BChE in an opposite way, i.e. inhibition of AChE and activation of BChE, regardless the nature of the reporter substrate.For both enzymes, progress curves for hydrolysis of PhA at very low concentration (?K_m) in the presence of increasing concentration of ATC showed that: a) the competing substrate and the reporter substrate are hydrolyzed at the same time, b) complete hydrolysis of PhA cannot be reached above 1 mM competing substrate. This likely results from accumulation of hydrolysis products (P) of competing substrate and/or accumulation of acetylated enzyme·P complex that inhibit hydrolysis of the reporter substrate.     

QM/MM study of the mechanism of enzymatic limonene 1,2-epoxide hydrolysis

Limonene 1,2-epoxide hydrolase (LEH) is completely different from those of classic epoxide hydrolases (EHs) which catalyze the hydrolysis of epoxides to vicinal diols. A novel concerted general acid catalysis step involving the Asp101-Arg99-Asp132 triad is proposed to play an important role in the mechanism. Combined quantum-mechanical/molecular-mechanical (QM/MM) calculations gave activation barriers of 16.9 and 25.1 kcal/mol at the B3LYP/6-31G(d,p)//CHARMM level for nucleophilic attack on the more and less substituted epoxide carbons, respectively. Furthermore, the important roles of residues Arg99, Tyr53 and Asn55 on mutated LEH were evaluated by QM/MM-scanned energy mapping. These results may provide an explanation for site-directed mutagenesis.     

Characterization of a complete cycle of acetylcholinesterase catalysis by <Emphasis Type="Italic">ab initio</Emphasis> QM/MM modeling

The reaction mechanism of acetylcholine hydrolysis by acetylcholinesterase, including both acylation and deacylation stages from the enzyme-substrate (ES) to the enzyme-product (EP) molecular complexes, is examined by using an ab initio type quantum mechanical – molecular mechanical (QM/MM) approach. The density functional theory PBE0/aug-6–31+G* method for a fairly large quantum part trapped inside the native protein environment, and the AMBER force field parameters in the molecular mechanical part are employed in computations. All reaction steps, including the formation of the first tetrahedral intermediate (TI1), the acylenzyme (EA) complex, the second tetrahedral intermediate (TI2), and the EP complex, are modeled at the same theoretical level. In agreement with the experimental rate constants, the estimated activation energy barrier of the deacylation stage is slightly higher than that for the acylation phase. The critical role of the non-triad Glu202 amino acid residue in orienting lytic water molecule and in stabilizing the second tetrahedral intermediate at the deacylation stage of the enzymatic process is demonstrated. Figure The computed energy diagram for the reaction path from the enzyme – substrate complex (ES) to the enzyme-product complex (EP).     

Catalytic mechanism of cellulose degradation by a cellobiohydrolase, CelS

The hydrolysis of cellulose is the bottleneck in cellulosic ethanol production. The cellobiohydrolase CelS from Clostridium thermocellum catalyzes the hydrolysis of cello-oligosaccharides via inversion of the anomeric carbon. Here, to examine key features of the CelS-catalyzed reaction, QM/MM (SCCDFTB/MM) simulations are performed. The calculated free energy profile for the reaction possesses a 19 kcal/mol barrier. The results confirm the role of active site residue Glu87 as the general acid catalyst in the cleavage reaction and show that Asp255 may act as the general base. A feasible position in the reactant state of the water molecule responsible for nucleophilic attack is identified. Sugar ring distortion as the reaction progresses is quantified. The results provide a computational approach that may complement the experimental design of more efficient enzymes for biofuel production.     

Exploring the RING-Catalyzed Ubiquitin Transfer Mechanism by MD and QM/MM Calculations

Ubiquitylation is a universal mechanism for controlling cellular functions. A large family of ubiquitin E3 ligases (E3) mediates Ubiquitin (Ub) modification. To facilitate Ub transfer, RING E3 ligases bind both the substrate and ubiquitin E2 conjugating enzyme (E2) linked to Ub via a thioester bond to form a catalytic complex. The mechanism of Ub transfer catalyzed by RING E3 remains elusive. By employing a combined computational approach including molecular modeling, molecular dynamics (MD) simulations, and quantum mechanics/molecular mechanics (QM/MM) calculations, we characterized this catalytic mechanism in detail. The three-dimensional model of dimeric RING E3 ligase RNF4 RING, E2 ligase UbcH5A, Ub and the substrate SUMO2 shows close contact between the substrate and Ub transfer catalytic center. Deprotonation of the substrate lysine by D117 on UbcH5A occurs with almost no energy barrier as calculated by MD and QM/MM calculations. Then, the side chain of the activated lysine gets close to the thioester bond via a conformation change. The Ub transfer pathway begins with a nucleophilic addition that forms an oxyanion intermediate of a 4.23 kcal/mol energy barrier followed by nucleophilic elimination, resulting in a Ub modified substrate by a 5.65 kcal/mol energy barrier. These results provide insight into the mechanism of RING-catalyzed Ub transfer guiding the discovery of Ub system inhibitors.     

Crystallographic evidence for active-site dynamics in the hydrolytic aldehyde dehydrogenases. Implications for the deacylation step of the catalyzed reaction

The overall chemical mechanism of the reaction catalyzed by the hydrolytic aldehyde dehydrogenases (ALDHs) involves three main steps: (1) nucleophilic attack of the thiol group of the catalytic cysteine on the carbonyl carbon of the aldehyde substrate; (2) hydride transfer from the tetrahedral thiohemiacetal intermediate to the pyridine ring of NAD(P)(+); and (3) hydrolysis of the resulting thioester intermediate (deacylation). Crystal structures of different ALDHs from several organisms-determined in the absence and presence of bound NAD(P)(+), NAD(P)H, aldehydes, or acid products-showed specific details at the atomic level about the catalytic residues involved in each of the catalytic steps. These structures also showed the conformational flexibility of the nicotinamide half of the cofactor, and of the catalytic cysteinyl and glutamyl residues, the latter being the general base that activates the hydrolytic water molecule in the deacylation step. The architecture of the ALDH active site allows for this conformational flexibility, which, undoubtedly, is crucial for catalysis in these enzymes. Focusing in the deacylation step of the ALDH-catalyzed reaction, here we review and systematize the crystallographic evidence of the structural features responsible for the conformational flexibility of the catalytic glutamyl residue, and for the positioning of the hydrolytic water molecule inside the ALDH active site. Based on the analysis of the available crystallographic data and of energy-minimized models of the thioester reaction intermediate, as well as on the results of theoretical calculations of the pK(a) of the carboxyl group of the catalytic glutamic acid in its three different conformations, we discuss the role that the conformational flexibility of this residue plays in the activation of the hydrolytic water. We also propose a critical participation in the water activation process of the peptide bond to which the catalytic glutamic acid in the intermediate conformation is hydrogen bonded.     

Density functional theory analysis of carbonyl sulfide hydrolysis: effect of solvation and nucleophile variation

The detailed mechanisms of the hydrolysis of carbonyl sulfide (OCS) by nucleophilic water and hydroxide ion in both the gas phase and bulk water solvent have been investigated using density functional theory. Various reaction channels on the potential surface have been identified. The thermodynamic results demonstrate that the hydrolysis of OCS by nucleophilic water and hydroxide ion should proceed more favorably at low temperature. The hydrolysis of OCS by the hydroxide ion is the main reaction channel from thermodynamic and kinetic perspectives, and the bulk solvent can influence the rate-determining step in this channel. However, the solvent barely modifies the activation energy of the rate-determining step. For the hydrolysis of OCS by nucleophilic water, the solvent does not modify the rate-determining step, and the corresponding activation energy of the rate-determining step barely changes. This bulk solvent effect suggests that most of the contribution of the solvent is accounted for by considering one water molecule and a hydroxide ion.     

Theoretical insights into the protonation states of active site cysteine and citrullination mechanism of Porphyromonas gingivalis peptidylarginine deiminase

Porphyromonas gingivalis peptidylarginine deiminase (PPAD) catalyzes the citrullination of peptidylarginine, which plays a critical role in the rheumatoid arthritis (RA) and gene regulation. For a better understanding of citrullination mechanism of PPAD, it is required to establish the protonation states of active site cysteine, which is still a controversial issue for the members of guanidino‐group‐modifying enzyme superfamily. In this work, we first explored the transformation between the two states: State N (both C351 and H236 are neutral) and State I (both residues exist as a thiolate–imidazolium ion pair), and then investigated the citrullination reaction of peptidylarginine, using a combined QM/MM approach. State N is calculated to be more stable than State I by 8.46 kcal/mol, and State N can transform to State I via two steps of substrate‐assisted proton transfer. Citrullination of the peptidylarginine contains deamination and hydrolysis. Starting from State N, the deamination reaction corresponds to an energy barrier of 18.82 kcal/mol. The deprotonated C351 initiates the nucleophilic attack to the substrate, which is the key step for deamination reaction. The hydrolysis reaction contains two chemical steps. Both the deprotonated D238 and H236 can act as the bases to activate the hydrolytic water, which correspond to similar energy barriers (～17 kcal/mol). On the basis of our calculations, C351, D238, and H236 constitute a catalytic triad, and their protonation states are critical for both the deamination and hydrolysis processes. In view of the sequence similarity, these findings may be shared with human PAD1–PAD4 and other guanidino‐group‐modifying enzymes. Proteins 2017; 85:1518–1528. © 2017 Wiley Periodicals, Inc.     

QM/MM modeling the Ras-GAP catalyzed hydrolysis of guanosine triphosphate

The mechanism of the hydrolysis reaction of guanosine triphosphate (GTP) by the protein complex Ras-GAP (p21(ras) - p120(GAP)) has been modeled by the quantum mechanical-molecular mechanical (QM/MM) and ab initio quantum calculations. Initial geometry configurations have been prompted by atomic coordinates of a structural analog (PDBID:1WQ1). It is shown that the minimum energy reaction path is consistent with an assumption of two-step chemical transformations. At the first stage, a unified motion of Arg789 of GAP, Gln61, Thr35 of Ras, and the lytic water molecule results in a substantial spatial separation of the gamma-phosphate group of GTP from the rest of the molecule (GDP). This phase of hydrolysis process proceeds through the low-barrier transition state TS1. At the second stage, Gln61 abstracts and releases protons within the subsystem including Gln61, the lytic water molecule and the gamma-phosphate group of GTP through the corresponding transition state TS2. Direct quantum calculations show that, in this particular environment, the reaction GTP + H(2)O --> GDP + H(2)PO(4) (-) can proceed with reasonable activation barriers of less than 15 kcal/mol at every stage. This conclusion leads to a better understanding of the anticatalytic effect of cancer-causing mutations of Ras, which has been debated in recent years.     

Partial Purification and Characterization of Soluble Isoform of Butyrylcholinesterase from Rat Intestine

Butyrylcholinesterase (BChE; E.C. 3.1.1.8.) was 260-fold purified from soluble fraction of rat intestine. The enzyme was composed of tetrameric globular form by nonreducing electrophoresis. Optimum pH value was determined as 7.2 after zero buffer extrapolation. Optimum temperature was examined as 37 degrees C after zero time extrapolation. The enzyme showed marked substrate activation with positively charged, acyl-choline substrates. As a measure of catalytic efficiency, kcat/Km values were determined as 16,210, 25,650, and 46,150 for acetylthiocholine (ATCh), propionylthiocholine (PTCh), and butyrylthiocholine (BTCh), respectively. When the catalytic efficiencies are compared, soluble isoform of rat intestinal BChE became increasingly efficient as the size of the acyl portion of the substrate increases; BTCh > PTCh > ATCh. Differently, the enzyme showed substrate inhibition with benzoylcholine (BzCh) and a kcat/Km value of 21,190 was found. Triton X-100 inhibited more efficiently the rat intestinal BChE soluble isoform than it did the human serum BChE.     

Carbinolamine formation and dehydration in a DNA repair enzyme active site

In order to suggest detailed mechanistic hypotheses for the formation and dehydration of a key carbinolamine intermediate in the T4 pyrimidine dimer glycosylase (T4PDG) reaction, we have investigated these reactions using steered molecular dynamics with a coupled quantum mechanics-molecular mechanics potential (QM/MM). We carried out simulations of DNA abasic site carbinolamine formation with and without a water molecule restrained to remain within the active site quantum region. We recovered potentials of mean force (PMF) from thirty replicate reaction trajectories using Jarzynski averaging. We demonstrated feasible pathways involving water, as well as those independent of water participation. The water-independent enzyme-catalyzed reaction had a bias-corrected Jarzynski-average barrier height of approximately (6.5 kcal mol(-1) (27.2 kJ mol(-1)) for the carbinolamine formation reaction and 44.5 kcal mol(-1) (186 kJ mol(-1)) for the reverse reaction at this level of representation. When the proton transfer was facilitated with an intrinsic quantum water, the barrier height was approximately 15 kcal mol(-1) (62.8 kJ mol(-1)) in the forward (formation) reaction and 19 kcal mol(-1) (79.5 kJ mol(-1)) for the reverse. In addition, two modes of unsteered (free dynamics) carbinolamine dehydration were observed: in one, the quantum water participated as an intermediate proton transfer species, and in the other, the active site protonated glutamate hydrogen was directly transferred to the carbinolamine oxygen. Water-independent unforced proton transfer from the protonated active site glutamate carboxyl to the unprotonated N-terminal amine was also observed. In summary, complex proton transfer events, some involving water intermediates, were studied in QM/MM simulations of T4PDG bound to a DNA abasic site. Imine carbinolamine formation was characterized using steered QM/MM molecular dynamics. Dehydration of the carbinolamine intermediate to form the final imine product was observed in free, unsteered, QM/MM dynamics simulations, as was unforced acid-base transfer between the active site carboxylate and the N-terminal amine.     

次氯酸与不饱和脂肪酸氧化机制和转化产物的理论研究

摘要 目的：揭示次氯酸与不饱和脂肪酸的氧化反应机制及转化产物。方法：运用Gaussian 16软件包,采用密度泛函方法M06-2X(D3),结合6-31+G(d)基组,在SMD液相水模型水平下进行计算。结果：次氯酸与单不饱和脂肪酸油酸的氧化反应是先形成氯鎓离子中间体,氯鎓离子再与水分子反应生成氯醇,第一步氯鎓离子的形成是控速步骤,其反应活化自由能~8 kcal/mol。环氧化合物和短链的醛是两种转化产物,前者由氯醇脱氯化氢而来,而后者由环氧化合物和氯醇通过系列与次氯酸根的反应而得到,生成它们的控速步骤的反应活化自由能分别为23 和24 kcal/mol。选取两个乙基为取代基的乙烯为油酸模型,其与次氯酸反应的活化自由能仅比油酸高1 kcal/mol。计算得到次氯酸与亚油酸、顺-9,反-11 亚油酸、梓树酸和花生四烯酸模型氧化反应生成氯醇的活化自由能分别是~10、13、16和14 kcal/mol。结论：氯鎓离子中间体机制是次氯酸与不饱和脂肪酸氧化反应的主要机制,反应的活化自由能通常低于15 kcal/mol,意味着此氧化反应动力学上容易发生。氧化产物氯醇能转化为环氧化合物和短链的醛,但活化自由能较高,约23和24 kcal/mol。选取距离双键3个碳以内的结构为不饱和脂肪酸模型,它能够很好地反映不饱和脂肪酸的反应活性。     

Catalytic mechanism and energy barriers for butyrylcholinesterase-catalyzed hydrolysis of cocaine

The geometries of the transition states, intermediates, and prereactive enzyme-substrate complex and the corresponding energy barriers have been determined by performing hybrid quantum mechanical/molecular mechanical (QM/MM) calculations on butyrylcholinesterase (BChE)-catalyzed hydrolysis of (-)- and (+)-cocaine. The energy barriers were evaluated by performing QM/MM calculations with the QM method at the MP2/6-31+G* level and the MM method using the AMBER force field. These calculations allow us to account for the protein environmental effects on the transition states and energy barriers of these enzymatic reactions, showing remarkable effects of the protein environment on intermolecular hydrogen bonding (with an oxyanion hole), which is crucial for the transition state stabilization and, therefore, on the energy barriers. The calculated energy barriers are consistent with available experimental kinetic data. The highest barrier calculated for BChE-catalyzed hydrolysis of (-)- and (+)-cocaine is associated with the third reaction step, but the energy barrier calculated for the first step is close to the highest and is so sensitive to the protein environment that the first reaction step can be rate determining for (-)-cocaine hydrolysis catalyzed by a BChE mutant. The computational results provide valuable insights into future design of BChE mutants with a higher catalytic activity for (-)-cocaine.