New insights into the mechanism of the Schiff base formation catalyzed by type I dehydroquinate dehydratase from <Emphasis Type="Italic">S. enterica</Emphasis>

The Schiff base formation catalyzed by type I dehydroquinate dehydratase (DHQD) from Salmonella enterica has been studied by molecular docking, molecular dynamics simulation, and quantum chemical calculations. The substrate locates stably a similar position as the Schiff base intermediate observed in the crystal structure and forms strong hydrogen bonds with several active site residues. This binding mode is different from that of several other Schiff base enzymes. Then, the quantum chemical model has been constructed and the fundamental reaction pathways have been explored by performing quantum chemical calculation. The energy barrier of the previously proposed reaction pathway is calculated to be 30.7 kcal/mol, which is much higher than the experimental value of 14.3 kcal/mol of the whole dehydration reaction by type I DHQD from S. enterica. It means that this pathway is not favorable in energy. Therefore, a new and unexpected reaction pathway has been investigated with the favorable and reasonable energy barrier of 12.1 kcal/mol. The complicated role of catalytic His143 residue has also been elucidated that it mediates two proton transfers to facilitate the reaction. Moreover, the similarity and the difference between these two reaction pathways have been analyzed in detail. The new structural and mechanistic insights may direct the design of the inhibitors of type I dehydroquinate dehydratase as non-toxic antimicrobials, antifungals, and herbicides.     

The reaction mechanism for dehydration process catalyzed by type I dehydroquinate dehydratase from Gram-negative Salmonella enterica

The fundamental reaction mechanism for the dehydration process catalyzed by type I dehydroquinate dehydratase from Gram-negative Salmonella enterica has been studied by density functional theory calculations. The results indicate that the dehydration process undergoes a two-step cis-elimination mechanism, which is different from the previously proposed one. The catalytic roles of both the highly conserved residue His143 and the Schiff base formed between the substrate and Lys170 have also been elucidated. The structural and mechanistic insight presented here may direct the design of type I dehydroquinate dehydratase enzyme inhibitors as non-toxic antimicrobials, anti-fungals, and herbicides.     

Catalytic roles of active-site amino acid residues of coenzyme B12-dependent diol dehydratase: protonation state of histidine and pull effect of glutamate

The hydrogen abstraction and the OH migration processes catalyzed by diol dehydratase are discussed by means of a quantum mechanical/molecular mechanical method. To evaluate the push effect of His143 and the pull effect of Glu170, we considered three kinds of whole-enzyme model, the protonated and two unprotonated His143 models. A calculated activation energy for the hydrogen abstraction by the adenosyl radical is 15.6 (13.6) kcal/mol in the protonated (unprotonated) His143 model. QM/MM calculational results show that the mechanism of the OH migration is significantly changed by the protonation of His143. In the protonated His143 model, the OH group migration triggered by the full proton donation from the imidazolium to the migrating OH group occurs by a stepwise OH abstraction/re-addition process in which the water production reduces the barrier for the C-O bond cleavage. On the other hand, the OH migration in the unprotonated His143 model proceeds in a concerted manner, as we previously proposed using a simple model including only K+ ion and substrate. The latter mechanism seems to be kinetically more favorable from the calculated energy profiles and is consistent with experimental results. The activation barrier of the OH group migration step is only 1.6 kcal/mol reduced by the hydrogen-bonding interaction between the O2 of the substrate and unprotonated His143. Thus, it is predicted that His143 is not protonated, and therefore the main active-site amino acid residue that lowers the energy of the transition state for the OH group migration is determined to be Glu170.     

Computational mutation analysis of hydrogen abstraction and radical rearrangement steps in the catalysis of coenzyme B12-dependent diol dehydratase

A mutation analysis of the catalytic functions of active-site residues of coenzyme B(12)-dependent diol dehydratase in the conversion of 1,2-propanediol to 1,1-propanediol has been carried out by using QM/MM computations. Mutants His143Ala, Glu170Gln, Glu170Ala, and Glu170Ala/Glu221Ala were considered to estimate the impact of the mutations of His143 and Glu170. In the His143Ala mutant the activation energy for OH migration increased to 16.4 from 11.5 kcal mol(-1) in the wild-type enzyme. The highest activation energy, 19.6 kcal mol(-1), was measured for hydrogen back-abstraction in this reaction. The transition state for OH migration is not sufficiently stabilized by the hydrogen-bonding interaction formed between the spectator OH group and Gln170 in the Glu170Gln mutant, which demonstrates that a strong proton acceptor is required to promote OH migration. In the Glu170Ala mutant, a new strong hydrogen bond is formed between the spectator OH group and Glu221. A computed activation energy of 13.6 kcal mol(-1) for OH migration in the Glu170Ala mutant is only 2.1 kcal mol(-1) higher than the corresponding barrier in the wild-type enzyme. Despite the low activation barrier, the Glu170Ala mutant is inactive because the subsequent hydrogen back-abstraction is energetically demanding in this mutant. OH migration is not feasible in the Glu170Ala/Glu221Ala mutant because the activation barrier for OH migration is greatly increased by the loss of COO(-) groups near the spectator OH group. This result indicates that the effect of partial deprotonation of the spectator OH group is the most important factor in reducing the activation barrier for OH migration in the conversion of 1,2-propanediol to 1,1-propanediol catalyzed by diol dehydratase.     

Molecular dynamics simulations of the catalytic pathway of a cysteine protease: a combined QM/MM study of human cathepsin K

Molecular dynamics simulations using a combined QM/MM potential have been performed to study the catalytic mechanism of human cathepsin K, a member of the papain family of cysteine proteases. We have determined the two-dimensional free energy surfaces of both acylation and deacylation steps to characterize the reaction mechanism. These free energy profiles show that the acylation step is rate limiting with a barrier height of 19.8 kcal/mol in human cathepsin K and of 29.3 kcal/mol in aqueous solution. The free energy of activation for the deacylation step is 16.7 kcal/mol in cathepsin K and 17.8 kcal/mol in aqueous solution. The reduction of free energy barrier is achieved by stabilization of the oxyanion in the transition state. Interestingly, although the "oxyanion hole" has been formed in the Michaelis complex, the amide units do not donate hydrogen bonds directly to the carbonyl oxygen of the substrate, but they stabilize the thiolate anion nucleophile. Hydrogen-bonding interactions are induced as the substrate amide group approaches the nucleophile, moving more than 2 A and placing the oxyanion in contact with Gln19 and the backbone amide of Cys25. The hydrolysis of peptide substrate shares a common mechanism both for the catalyzed reaction in human cathepsin K and for the uncatalyzed reaction in water. Overall, the nucleophilic attack by Cys25 thiolate and the proton-transfer reaction from His162 to the amide nitrogen are highly coupled, whereas a tetrahedral intermediate is formed along the nucleophilic reaction pathway.     

Modeling the catalysis of anti-cocaine catalytic antibody: competing reaction pathways and free energy barriers

The competing reaction pathways and the corresponding free energy barriers for cocaine hydrolysis catalyzed by an anti-cocaine catalytic antibody, mAb15A10, were studied by using a novel computational strategy based on the binding free energy calculations on the antibody binding with cocaine and transition states. The calculated binding free energies were used to evaluate the free energy barrier shift from the cocaine hydrolysis in water to the antibody-catalyzed cocaine hydrolysis for each reaction pathway. The free energy barriers for the antibody-catalyzed cocaine hydrolysis were predicted to be the corresponding free energy barriers for the cocaine hydrolysis in water plus the calculated free energy barrier shifts. The calculated free energy barrier shift of -6.87 kcal/mol from the dominant reaction pathway of the cocaine benzoyl ester hydrolysis in water to the dominant reaction pathway of the antibody-catalyzed cocaine hydrolysis is in good agreement with the experimentally derived free energy barrier shift of -5.93 kcal/mol. The calculated mutation-caused shifts of the free energy barrier are also reasonably close to the available experimental activity data. The good agreement suggests that the protocol for calculating the free energy barrier shift from the cocaine hydrolysis in water to the antibody-catalyzed cocaine hydrolysis may be used in future rational design of possible high-activity mutants of the antibody as anti-cocaine therapeutics. The general strategy of the free energy barrier shift calculation may also be valuable in studying a variety of chemical reactions catalyzed by other antibodies or proteins through noncovalent bonding interactions with the substrates.     

A QM/MM simulation study of transamination reaction at the active site of aspartate aminotransferase: Free energy landscape and proton transfer pathways

Transaminase is a key enzyme for amino acid metabolism, which reversibly catalyzes the transamination reaction with the help of PLP (pyridoxal 5^' -phosphate) as its cofactor. Here we have investigated the mechanism and free energy landscape of the transamination reaction involving the aspartate transaminase (AspTase) enzyme and aspartate-PLP (Asp-PLP) complex using QM/MM simulation and metadynamics methods. The reaction is found to follow a stepwise mechanism where the active site residue Lys258 acts as a base to shuttle a proton from α -carbon (CA) to imine carbon (C4A) of the PLP-Asp Schiff base. In the first step, the Lys258 abstracts the CA proton of the substrate leading to the formation of a carbanionic intermediate which is followed by the reprotonation of the Asp-PLP Schiff base at C4A atom by Lys258. It is found that the free energy barrier for the proton abstraction by Lys258 and that for the reprotonation are 17.85 and 3.57 kcal/mol, respectively. The carbanionic intermediate is 7.14 kcal/mol higher in energy than the reactant. Hence, the first step acts as the rate limiting step. The present calculations also show that the Lys258 residue undergoes a conformational change after the first step of transamination reaction and becomes proximal to C4A atom of the Asp-PLP Schiff base to favor the second step. The active site residues Tyr70* and Gly38 anchor the Lys258 in proper position and orientation during the first step of the reaction and stabilize the positive charge over Lys258 generated at the intermediate step.     

Catalytic mechanism of matrix metalloproteinases: two-layered ONIOM study

The two-layered ONIOM(B3LYP:MNDO) method has been used to investigate the hydrolytical mechanism of matrix metalloproteinases (MMPs), a large family of zinc-dependent endopeptidases capable of degrading a wide range of macromolecules of the extracellular matrix. Human stromelysin-1 (MMP-3) was chosen as a physiologically important member of the MMP family. As a structural reference, X-ray data on the stromelysin-1 catalytic domain (SCD) complexed to the transition state analogue diphenyl piperidine sulfonamide inhibitor was used. The backbone spacer of 11 residues (201-211) was included in the final model, spanning the catalytic Glu202 residue and the three structural His201,205,211 zinc ligands. The polypeptide framework incorporated, partly accounting for the protein rigidity, reduces the activation free energy slightly by 1.6 kcal/mol. Essentially a single-step catalytic mechanism was obtained, generally following a classical proposal for MMPs. Glu202 here acts as a base, abstracting a proton from the metal-bound reactant water and delivering this proton to the peptide nitrogen. An auxiliary water molecule is suggested to be of crucial importance acting as an electrophilic agent to the carbonyl oxygen of the substrate. The direct inclusion of the auxiliary water molecule decreases the activation free energy by about 5 kcal/mol via donation of a strong hydrogen bond. The calculated activation barrier of 13.1 kcal/mol agrees well with experimental rates.     

Acetylcholinesterase: mechanisms of covalent inhibition of wild-type and H447I mutant determined by computational analyses

The reaction mechanisms of two inhibitors TFK+ and TFK0 binding to both the wild-type and H447I mutant mouse acetylcholinesterase (mAChE) have been investigated by using a combined ab initio quantum mechanical/molecular mechanical (QM/MM) approach and classical molecular dynamics (MD) simulations. In the wild-type mAChE, the binding reactions of TFK+ and TFK0 are both spontaneous processes, which proceed through the nucleophilic addition of the Ser203-Ogamma to the carbonyl-C of TFK+ or TFK0, accompanied with a simultaneous proton transfer from Ser203 to His447. No barrier is found along the reaction paths, consistent with the experimental reaction rates approaching the diffusion-controlled limit. By contrast, TFK+ binding to the H447I mutant may proceed with a different reaction mechanism. A water molecule takes over the role of His447 and participates in the bond breaking and forming as a "charge relayer". Unlike in the wild-type mAChE case, Glu334, a conserved residue from the catalytic triad, acts as a catalytic base in the reaction. The calculated energy barrier for this reaction is about 8 kcal/mol. These predictions await experimental verification. In the case of the neutral ligand TFK0, however, multiple MD simulations on the TFK0/H447I complex reveal that none of the water molecules can be retained in the active site as a "catalytic" water. Furthermore, our alchemical free energy calculation also suggests that the binding of TFK0 to H447I is much weaker than that of TFK+. Taken together, our computational studies confirm that TFK0 is almost inactive in the H447I mutant and also provide detailed mechanistic insights into the experimental observations.     

Probing O-dealkylation and deamination aging processes in tabun-conjugated AChE: a computational study

The mechanisms of the aging process of tabun-conjugated acetylcholinesterase were explored using density functional theory calculations. The free energy surfaces were calculated for O-dealkylation (C–O bond breaking) and deamination (P–N bond breaking) pathways for the aging process of tabun-conjugated acetylcholinesterase as suggested by mass and crystallographic studies. Initially, the calculations were performed using tabun-conjugated serine (SUN) molecule. O-dealkylation mechanism proceeds via one-step S_N2 type process, whereas the deamination process proceeds via two steps addition–elimination reaction at the phosphorus center of SUN molecule. The recent proposal of another deamination mechanism using human butyrylcholinesterase (hBChE) conjugated with N-mono methyl analogue of tabun (TA4) has also been explored (Nachon et al. in Chem Biol Interact 187:44–48, 2010). The rate-determining activation barrier calculated for this deamination mechanism (26.3 kcal/mol) was comparable with O-dealkylation process (26.9 kcal/mol) with B3LYP/6-31+G* level of theory. To examine the influence of catalytic residue His447, additional calculations were performed with imidazole group of His447 residue. The incorporation of imidazole group of catalytic residue His447 showed marked decrease in the free energies of activation for all the studied aging processes of tabun-inhibited serine. The aging mechanisms have been explored with TA4-inhibited serine, and calculated results showed that the deamination with the rearrangement process is markedly preferred in this case, which supports the Nachon et al. proposal based on the crystallographic studies.     

Ab initio and density functional theory studies of the catalytic mechanism for ester hydrolysis in serine hydrolases

We present results from ab initio and density functional theory studies of the mechanism for serine hydrolase catalyzed ester hydrolysis. A model system containing both the catalytic triad and the oxyanion hole was studied. The catalytic triad was represented by formate anion, imidazole, and methanol. The oxyanion hole was represented by two water molecules. Methyl formate was used as the substrate. In the acylation step, our computations show that the cooperation of the Asp group and oxyanion hydrogen bonds is capable of lowering the activation barrier by about 15 kcal/mol. The transition state leading to the first tetrahedral intermediate in the acylation step is rate limiting with an activation barrier (ΔE₀) of 13.4 kcal/mol. The activation barrier in the deacylation step is smaller. The double-proton-transfer mechanism is energetically unfavorable by about 2 kcal/mol. The bonds between the Asp group and the His group, and the hydrogen bonds in the oxyanion hole, increase in strength going from the Michaelis complex toward the transition state and the tetrahedral intermediate. In the acylation step, the tetrahedral intermediate is a very shallow minimum on the energy surface and is not viable when molecular vibrations are included. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 69: 89–103, 1998     

A semiempirical study of the optimized ground and excited state potential energy surfaces of retinal and its protonated Schiff base.

The electronic ground and first excited states of retinal and its Schiff base are optimized for the first time using the semiempirical AM1 Hamiltonian. The barrier for rotation about the C(11)-C(12) double bond is characterized by variation of both the twist angle delta(C(10)-C(11)-C(12)-C(13)) and the bond length d(C(11)-C(12)). The potential energy surface is obtained by varying these two parameters. The calculated ground state rotational barrier is equal to 15.6 kcal/mol for retinal and 20.5 kcal/mol for its Schiff base. The all-trans conformation is more stable by 3.7 kcal/mol than the 11-cis geometry. For the first excited state, S(1,) the 90 degrees twisted geometry represents a saddle point for retinal with the rotational barrier of 14.6 kcal/mol. In contrast, this conformation is an energy minimum for the Schiff base. It can be easily reached at room temperature from the planar minima since it is separated from them by a barrier of only 0.6 kcal/mol. The 90 degrees minimum conformation is more stable than the all-trans by 8.6 kcal/mol. We are thus able to present a reaction path on the S(1) surface of the retinal Schiff base with an almost barrier-less geometrical relaxation into a twisted minimum geometry, as observed experimentally. The character of the ground and first excited singlet states underscores the need for the inclusion of double excitations in the calculations.     

Density functional calculations of ATP systems. 2. ATP hydrolysis at the active site of actin

The hydrolysis of adenosine 5'-triphosphate (ATP) at the active site of actin has been studied using density functional calculations. The active site is modeled by the triphosphate tail of ATP, an Mg cation, surrounding water molecules, and the nearby protein residues. Four reaction paths have been followed by constraining coordinates that allow phosphate stretching, nucleophilic attack of the catalytic water, and OH(-) formation via water deprotonation. The lowest-energy barrier (21.0 kcal/mol) is obtained for a dissociative reaction where the terminal phosphate breaks on approaching the catalytic water, followed by proton release via a proton wire mechanism. A higher barrier (39.6 kcal/mol) results for an associative reaction path where OH(-) is formed first, with a pentacoordinated phosphorus atom (P-O distances 2.1 A). Stretching the terminal bridging P-O bond results in bond rupture at 2.8 A with an energy barrier of 28.8 kcal/mol. The residues Gln137 and His161 are not important in the reactions, but insight into their roles in vivo has been obtained. The favored coordination of the end products H(2)PO(4)(-) and ADP(3-) includes a hydrogen bond and an O-Mg-O bridge between the phosphates as well as a hydrogen bond between H(2)PO(4)(-) and the Ser14 side chain. The total energy is 2.1 kcal/mol lower than in the initial reactants. Classical simulations of ATP- and ADP.P(i)-actin show few hydrolysis-induced differences in the protein structure, indicating that phosphate migration is necessary for a change in conformation.     

人工设计逆醛缩酶RA95.5-8F催化β-羟基酮C—C裂解的理论研究

采用量子力学与分子力学组合(QM/MM)方法对人工设计逆醛缩酶RA95.5-8F催化β-羟基酮化合物裂解反应的机理进行了研究.结果表明,裂解反应主要包括赖氨酸Lys1083对底物的亲核进攻、Schiff碱形成、烯胺水解及C—N断裂等过程, C—N键裂解生成丙酮为整个反应的决速步骤,能垒为106.27 kJ/mol;活性中心的赖氨酸Lys1083、酪氨酸Tyr1051、天冬酰胺Asn1110和酪氨酸Tyr1180构成一个催化四联体, Lys1083通过与底物形成席夫碱对底物进行活化, Tyr1051作为催化酸碱参与质子转移过程,催化四联体的氢键网络有利于反应过渡态的稳定并使R-构型的底物更容易结合在活性位点,导致RA95.5-8F对R构型底物具有高的选择性和催化活性.     

Peptide hydrolysis catalyzed by matrix metalloproteinase 2: a computational study

The MMP-2 reaction mechanism is investigated by using different computational methodologies. First, quantum mechanical (QM) calculations are carried out on a cluster model of the active site bound to an Ace-Gly approximately Ile-Nme peptide. Along the QM reaction path, a Zn-bound water molecule attacks the Gly carbonyl group to give a tetrahedral intermediate. The breaking of the C-N bond is completed thanks to the Glu 404 residue that shuttles a proton from the water molecule to Ile-N atom. The gas-phase QM energy barrier is quite low ( approximately 14 kcal/mol), thus suggesting that the essential catalytic machinery is included in the cluster model. A similar reaction path occurs in the MMP-2 catalytic domain bound to an octapeptide substrate according to hybrid QM and molecular mechanical (QM/MM) geometry optimizations. However, the rupture of the Gly( P 1) approximately Ile( P 1') amide bond is destabilized in the static QM/MM calculations, owing to the positioning of the Ile( P 1') side chain inside the MMP-2 S 1' pocket and to the inability of simple energy miminization methodologies to properly relax complex systems. Molecular dynamics simulations show that these steric limitations are overcome easily through structural fluctuations. The energetic effect of structural fluctuations is taken into account by combining QM energies with average MM Poisson-Boltzmann free energies, resulting in a total free energy barrier of 14.8 kcal/mol in good agreement with experimental data. The rate-determining event in the MMP-2 mechanism corresponds to a H-bond rearrangement involving the Glu 404 residue and/or the Glu 404-COOH --> N-Ile( P 1') proton transfer. Overall, the present computational results and previous experimental data complement each other well in order to provide a detailed view of the MMPs catalytic mechanism.     

Computational study of the phosphoryl transfer catalyzed by a cyclin-dependent kinase

A cyclin-dependent kinase, Cdk2, catalyzes the transfer of the gamma-phosphate from ATP to a threonine or serine residue of its polypeptide substrates. Here, we investigate aspects of the reaction mechanism of Cdk2 by gas-phase density functional calculations, classical molecular dynamics, and Car-Parrinello QM/MM simulations. We focus on the role of the conserved Asp127 and on the nature of the phosphoryl transfer reaction mechanism catalyzed by Cdk2. Our findings suggest that Asp127 is active in its deprotonated form by assisting the formation of the near-attack orientation of the substrate serine or threonine. Therefore, the residue does not act as a general base during the catalysis. The mechanism for the phosphoryl transfer is a single SN2-like concerted step, which shows a phosphorane-like transition state geometry. Although the resulting reaction mechanism is in agreement with a previous density functional study of the same catalytic reaction mechanism (Cavalli et al., Chem. Comm. 2003, 1308-1309), the reaction barrier is considerably lower when QM/MM calculations are performed, as in this study ( approximately 42 kcal mol(-1) QM vs. approximately 24 kcal mol(-1) QM/MM); this indicates that important roles for the catalysis are played by the protein environment and solvent waters. Because of the high amino acid sequence conservation among the whole family of cyclin-dependent kinases (CDKs), these results could be general for the CDK family.     

A Mechanistic Study of the Spontaneous Hydrolysis of Glycylserine as the Simplest Model for Protein Self‐Cleavage

A common feature of several classes of intrinsically reactive proteins with diverse biological functions is that they undergo self‐catalyzed reactions initiated by an N→O or N→S acyl shift of a peptide bond adjacent to a serine, threonine, or cysteine residue. In this study, we examine the N→O acyl shift initiated peptide‐bond hydrolysis at the serine residue on a model compound, glycylserine (GlySer), by means of DFT and ab initio methods. In the most favorable rate‐determining transition state, the serine ?COO^? group acts as a general base to accept a proton from the attacking ?OH function, which results in oxyoxazolidine ring closure. The calculated activation energy (29.4 kcal mol^?1) is in excellent agreement with the experimental value, 29.4 kcal mol^?1, determined by ¹H NMR measurements. A reaction mechanism for the entire process of GlySer dipeptide hydrolysis is also proposed. In the case of proteins, we found that when no other groups that may act as a general base are available, the N→O acyl shift mechanism might instead involve a water‐assisted proton transfer from the attacking serine ?OH group to the amide oxygen. However, the calculated energy barrier for this process is relatively high (33.6 kcal mol^?1), thus indicating that in absence of catalytic factors the peptide bond adjacent to serine is no longer a weak point in the protein backbone. An analogous rearrangement involving the amide N‐protonated form, rather than the principle zwitterion form of GlySer, was also considered as a model for the previously proposed mechanism of sea‐urchin sperm protein, enterokinase, and agrin (SEA) domain autoproteolysis. The calculated activation energy (14.3 kcal mol^?1) is significantly lower than the experimental value reported for SEA (≈21 kcal mol^?1), but is still in better agreement as compared to earlier theoretical attempts.     

(H_2O)_n(n=0～2)催化HS和HOCl反应机理的理论研究

在CCSD(T)/6-311+G(3df,2p)//M06-2X/6-311+G(3df,2p)水平上研究了(H_2O)n(n=0~2)催化HS和HOCl的反应机理.结果表明,HS与HOCl反应中HS夺取HOCl上的H原子形成产物H_2S和ClO.在无水催化时,该反应存在2种不同的路径(分别经过过渡态TS1和TS2,二者互为顺反结构),对应的能垒分别为100.28和100.91kJ/mol,到达产物(H_2S+ClO)需吸收18.99kJ/mol能量,反应不易发生;在单个水分子参与时,水分子可通过形成弱相互作用或者作为H原子转移桥梁影响反应机理,获得了4种水催化路径,能垒(间于53.97~92.39kJ/mol之间)均低于无水催化过程.同时发现,在反应到达产物前,水分子可以与产物形成中间体IM,IM相对能仅为0.46kJ/mol,有利于产物形成;有2个水分子参与反应时,找到了3条催化路径,最优反应路径过渡态TS7的能垒为45.05kJ/mol,低于无水催化过程,相比单个水分子最优路径能垒(53.97kJ/mol)并无显著降低.