Catalytic mechanism of the inverting N-acetylglucosaminyltransferase I: DFT quantum mechanical model of the reaction pathway and determination of the transition state structure

The complex N-glycan structures on glycoproteins play important roles in cell adhesion and recognition events in metazoan organisms. A critical step in the biosynthetic pathway leading from high mannose to these complex structures includes the transfer of N-acetylglucosamine (GlcNAc) to a mannose residue by the inverting N-acetylglucosaminyltransferase I (GnT-I). The catalytic mechanism of this enzymatic reaction is explored herein using DFT quantum chemical methods. The computational model used to follow the reaction is based on the X-ray crystallographic structure of GnT-I and contains 127 atoms that represent fragments of residues critical for the substrate binding and catalysis. The mechanism of the catalytic reaction was monitored by means of a 2D potential energy map calculated as a function of predefined reaction coordinates at the B3LYP/6-31G** level. This potential energy surface revealed one transition state associated with a reaction pathway following a concerted mechanism. The reaction barrier was estimated, and the structure of the transition state was characterized at the B3LYP/6-311++G**// B3LYP/6-31G** level.     

Structural Analysis of Thermus thermophilus HB27 Mannosyl-3-phosphoglycerate Synthase Provides Evidence for a Second Catalytic Metal Ion and New Insight into the Retaining Mechanism of Glycosyltransferases

Mannosyl-3-phosphoglycerate synthase is a glycosyltransferase involved in the two-step synthetic pathway of mannosylglycerate, a compatible solute that accumulates in response to salt and/or heat stresses in many microorganisms thriving in hot environments. The three-dimensional structure of mannosyl-3-phosphoglycerate synthase from Thermus thermophilus HB27 in its binary complex form, with GDP-α-d-mannose and Mg²⁺, shows a second metal binding site, about 6 Å away from the mannose moiety. Kinetic and mutagenesis studies have shown that this metal site plays a role in catalysis. Additionally, Asp¹⁶⁷ in the DXD motif is found within van der Waals contact distance of the C1′ atom in the mannopyranose ring, suggesting its action as a catalytic nucleophile, either in the formation of a glycosyl-enzyme intermediate according to the double-displacement S_N2 reaction mechanism or in the stabilization of the oxocarbenium ion-like intermediate according to the D_N*A_Nss (S_Ni-like) reaction mechanism. We propose that either mechanism may occur in retaining glycosyltransferases with a GT-A fold, and, based on the gathered structural information, we identified an extended structural signature toward a common scaffold between the inverting and retaining glycosyltransferases.     

A DFT study on the mechanism of the organocatalytic synthesis of a benzoxazine-substituted indolizine derivative

A series of theoretical computations were conducted via density functional theory at the B3LYP/6–31 + G(d,p) level to study the mechanism of the organocatalytic synthesis of a benzoxazine-substituted indolizine derivative. Four possible alternative pathways were considered in this work. The calculated results show that the formation of an N-ylide precursor from 4-dimethylaminopyridine (DMAP) is a key step as it provides the necessary nucleophilic centre for the subsequent H-migration and H-elimination processes. The precursor N-ylide and Schiff base isomers with the most favourable activities in the preliminary work were identified theoretically by analysing the reaction mechanism. The synthetic mechanism to obtain the indolizine derivative was found to be a two-step reaction, with the rate-determining step being the first H migration to form a transition state with a four-membered ring. The catalytic activity of DMAP in the first H-migration step in the overall synthetic process greatly reduces the reaction barrier height. The chiral selectivity of the synthesis is dominated by the spatial geometry of the Schiff base functional group.     

Mutagenesis and mechanistic study of a glycoside hydrolase family 54 alpha-L-arabinofuranosidase from Trichoderma koningii

A GH (glycoside hydrolase) family 54 alpha-L-arabinofuranosidase from Trichoderma koningii G-39 (termed Abf) was successfully expressed in Pichia pastoris and purified to near homogeneity by cation-exchange chromatography. To determine the amino acid residues essential for the catalytic activity of Abf, extensive mutagenesis of 24 conserved glutamate and aspartate residues was performed. Among the mutants, D221N, E223Q and D299N were found to decrease catalytic activity significantly. The kcat values of the D221N and D299N mutants were 7000- and 1300-fold lower respectively, than that of the wild-type Abf. E223Q was nearly inactive. These results are consistent with observations obtained from the Aspergillus kawachii alpha-L-arabinofuranosidase three-dimensional structure. This structure indicates that Asp221 of T. koningii Abf is significant for substrate binding and that Glu223 as well as Asp299 function as a nucleophile and a general acid/base catalyst for the enzymatic reaction respectively. The catalytic mechanism of wild-type Abf was further investigated by NMR spectroscopy and kinetic analysis. The results showed that Abf is a retaining enzyme. It catalyses the hydrolysis of various substrates via the formation of a common intermediate that is probably an arabinosyl-enzyme intermediate. A two-step, double-displacement mechanism involving first the formation, and then the breakdown, of an arabinosyl-enzyme intermediate was proposed. Based on the kcat values of a series of aryl-alpha-L-arabinofuranosides catalytically hydrolysed by wild-type Abf, a relatively small Br?nsted constant, beta(lg)=-0.18, was obtained, suggesting that the rate-limiting step of the enzymatic reaction is the dearabinosylation step. Further kinetic studies with the D299G mutant revealed that the catalytic activity of this mutant depended largely on the pK(a) values (>6) of leaving phenols, with beta(lg)=-1.3, indicating that the rate-limiting step of the reaction becomes the arabinosylation step. This kinetic outcome supports the idea that Asp299 is the general acid/base residue. The pH activity profile of D299N provided further evidence strengthening this suggestion.     

Computational study of the Sonogashira cross-coupling reaction in the gas phase and in dichloromethane solution

The Sonogashira cross-couplig reaction, consisting of oxidative addition, cis-trans isomerization, transmetalation, and reductive elimination, was computationally modeled using the DFT B3LYP/cc-pVDZ method for reaction between bromobenzene and phenylacetylene. Palladium diphosphane was used as a catalyst, copper(I) bromide as a co-catalyst and trimethylamine as a base. The reaction mechanism was studied both in the gas phase and in dichloromethane solution using PCM method. The complete catalytic cycle is thermodynamically strongly shifted toward products (diphenylacetylene and regenerated palladium catalyst) and is exothermic being in accordance with experimental data. The rate-determining step is the oxidative addition, since the highest point on the Gibbs energy graph of the complete reaction is the transition state of this step. This conclusion is also supported by recent experimental data. The computed energy profile suggests that the transmetalation step is initiated by the dissociation of neutral ligand, while the activation Gibbs energy of this step is 0.1 kcal mol(-1) in the gas phase.     

Theoretical modeling and experimental studies on N-n-Decyl-2-oxo-5-nitro-1-benzylidene-methylamine

The Schiff base compound, N-n-Decyl-2-oxo-5-nitro-1-benzylidene-methylamine, has been -synthesized and characterized by IR, electronic spectroscopy, and X-ray single-crystal determination. Molecular geometry from X-ray experiment of the title compound in the ground state have been compared using the Hartree-Fock (HF) and density functional method (B3LYP) with 6-31G(d) basis set. Calculated results show that density functional theory (DFT) at B3LYP/6-31G(d) level can well reproduce the structure of the title compound. To investigate the solvent effect for the atomic charge distributions of the title compound, self-consistent reaction field theory with Onsager reaction field model was used. In addition, DFT calculations of the title compound, molecular electrostatic potential and thermodynamic properties were performed at B3LYP/6-31G(d) level of theory.     

Structural evidence of a passive base-flipping mechanism for AGT, an unusual GT-B glycosyltransferase

The Escherichia coli T4 bacteriophage uses two glycosyltransferases to glucosylate and thus protect its DNA: the retaining alpha-glucosyltransferase (AGT) and the inverting beta-glucosyltransferase (BGT). They glucosylate 5-hydroxymethyl cytosine (5-HMC) bases of duplex DNA using UDP-glucose as the sugar donor to form an alpha-glucosidic linkage and a beta-glucosidic linkage, respectively. Five structures of AGT have been determined: a binary complex with the UDP product and four ternary complexes with UDP or UDP-glucose and oligonucleotides containing an A:G, HMU:G (hydroxymethyl uracyl) or AP:G (apurinic/apyrimidinic) mismatch at the target base-pair. AGT adopts the GT-B fold, one of the two folds known for GTs. However, while the sugar donor binding mode is classical for a GT-B enzyme, the sugar acceptor binding mode is unexpected and breaks the established consensus: AGT is the first GT-B enzyme that predominantly binds both the sugar donor and acceptor to the C-terminal domain. Its active site pocket is highly similar to four retaining GT-B glycosyltransferases (trehalose-6-phosphate synthase, glycogen synthase, glycogen and maltodextrin phosphorylases) strongly suggesting a common evolutionary origin and catalytic mechanism for these enzymes. Structure-guided mutagenesis and kinetic analysis do not permit identification of a nucleophile residue responsible for a glycosyl-enzyme intermediate for the classical double displacement mechanism. Interestingly, the DNA structures reveal partially flipped-out bases. They provide evidence for a passive role of AGT in the base-flipping mechanism and for its specific recognition of the acceptor base.     

Molecular modeling insights into the catalytic mechanism of the retaining galactosyltransferase LgtC

The bacterial enzyme lipopolysaccharyl alpha-galactosyltransferase C (EC 2.4.1.x, LgtC) is involved in the synthesis of lipooligosaccharides displayed on the cell surfaces of Neisseria meningitidis. LgtC catalyzes the transfer of a galactosyl residue from UDP-Gal to the terminal galactose residue of glycoconjugates with an overall retention of stereochemistry at the anomeric center. Several hypothetical catalytic mechanisms of the LgtC enzyme were examined herein using DFT quantum chemical methods up to the B3LYP/6-311++G**//B3LYP/6-31G* level. The computational model used to follow the reaction is based on the crystallographic structure of LgtC in complex with both the nucleotide-galactose donor and the oligosaccharide-acceptor analogues. The 136 atoms included in this model represent fragments of residues critical for the substrate binding and catalysis. From our calculations, the preferred pathway is predicted to be a one step mechanism with the nucleophilic attack of the acceptor oxygen onto the anomeric carbon and the proton transfer to a phosphate oxygen occurring simultaneously. This mechanism has an A(N)D(N)A(H)D(H) character, with the unique transition state structure in which the attacking galactose group is more closely bound to the anomeric carbon than to the UDP leaving group and where the hydrogen bond between the nucleophile and the leaving group oxygens facilitates the attack of the acceptor O4(') from the same side of the transferred galactose.     

Potential transition-state analogs for glycosyltransferases. Design and DFT calculations of conformational behavior

The structure of a previously calculated transition state (TS) was used to design the [tetrahydro-2-(methylthio)furan-2-yl]methyl phosphate dianion (1) as a new scaffold for transition-state analogs of reactions catalyzed by the inverting glycosyltransferases. This scaffold contains relevant features of the donor and acceptor and represents a new type of potential inhibitors for these enzymes. Available conformational space of 1 was explored using DFT quantum chemical methods by means of two-dimensional potential-energy maps calculated as a function of Phi, Psi, and omega dihedral angles at the B3LYP/6-31+G* level. The calculated potential energy surfaces revealed the existence of several low-energy domains. Structures from these regions were refined at the 6-311++G** level and led to 14 conformers. The stability of conformers is influenced by their environment, and in aqueous solution two conformers dominate the equilibrium. A superposition of calculated conformers with the predicted TS structure revealed that the preferred conformers in solution nicely mimic structural features of the TS. These results imply that 1 has structural properties required to mimic the TS and therefore can be used as a scaffold for further development of TS-analog inhibitors for retaining glycosyltransferases.     

A density functional theory study of the decomposition mechanism of nitroglycerin

The detailed decomposition mechanism of nitroglycerin (NG) in the gas phase was studied by examining reaction pathways using density functional theory (DFT) and canonical variational transition state theory combined with a small-curvature tunneling correction (CVT/SCT). The mechanism of NG autocatalytic decomposition was investigated at the B3LYP/6-31G(d,p) level of theory. Five possible decomposition pathways involving NG were identified and the rate constants for the pathways at temperatures ranging from 200 to 1000 K were calculated using CVT/SCT. There was found to be a lower energy barrier to the β-H abstraction reaction than to the α-H abstraction reaction during the initial step in the autocatalytic decomposition of NG. The decomposition pathways for CHOCOCHONO₂ (a product obtained following the abstraction of three H atoms from NG by NO₂) include O–NO₂ cleavage or isomer production, meaning that the autocatalytic decomposition of NG has two reaction pathways, both of which are exothermic. The rate constants for these two reaction pathways are greater than the rate constants for the three pathways corresponding to unimolecular NG decomposition. The overall process of NG decomposition can be divided into two stages based on the NO₂ concentration, which affects the decomposition products and reactions. In the first stage, the reaction pathway corresponding to O–NO₂ cleavage is the main pathway, but the rates of the two autocatalytic decomposition pathways increase with increasing NO₂ concentration. However, when a threshold NO₂ concentration is reached, the NG decomposition process enters its second stage, with the two pathways for NG autocatalytic decomposition becoming the main and secondary reaction pathways.     

A Comparative Semiempirical, Ab initio and DFT Study of Decarbonylation of Cyclic Ketones Part 1: Thermal Fragmentation of Cyclopropanone

Calculations using different quantum mechanical methods including semiempirical (MNDO,AM1 and PM3), ab initio (RHF and MP2 calculations using the 6-311G and 6-311++G^** basis sets), and density functional theory (LSDA, BP, MIXBP and B3LYP, i.e., B3LYP/6-311+G**//B3LYP/6-31G*) have been performed on the thermal fragmentation of cyclopropanone to ethylene and carbon monoxide. All RHF calculations predict a concerted single step mechanism for this conversion. The estimated activation energies vary from 34.4 to 54.6 kcal·mol^-1, mainly localized around 37±2 kcal·mol^-1, depending on the method. Whereas the calculated RHF reaction energies also varied from 14.5 to -33.3 kcal·mol^-1, the B3LYP/6-311+G**//B3LYP/6-31G* method predicts the experimental value (-17.7 kcal·mol^-1) within experimental uncertainties. Remarkably, semiempirical AM1 and PM3 methods and simple DFT calculations, LSDA, predict comparable results to the more advanced methods. UHF ab initio calculations predict the same single step mechanism, whereas a multistep biradical mechanism with an unrealistically low activation energy is favored by the semiempirical methods. Structures of the activated complex of the single step mechanism, estimated by different methods, are very similar and consistent with a nonlinear cheletropic [_2s + _2a] reaction, as predicted by the orbital symmetry rules and earlier EHT calculations.Electronic Supplementary Material available.     

A theoretical study on the activation of Ser70 in the acylation mechanism of cephalosporin antibiotics

A computational study using some molecular modeling and quantum mechanical methods has been performed for determining the most favor activation process for Ser70 in the acylation mechanism for the cephalosporin antibiotics among the three proposed ones given in the literature. The computation is based on an X-ray structure of the B chain of the Bacillus licheniformis BS3 beta-lactamase-cefoxitin complex. The position of a catalytic water involved in one of the reaction mechanism is defined using the Grid20 and InsightII programs, while that of the truncated ligand is defined using the InsightII and FirstDiscovery programs. The geometry of structures of each reaction scheme is optimized at the HF/3-21 G level of theory, and then the single point energy for each reactive species in each reaction scheme is computed at the levels of HF/6-31 + G (3df, 2p) and B3LYP/6-31 + G (3df, 2p). The effect of solvent on each reaction scheme is also studied by comparing the calculation results for each reaction scheme either in gas phase or in solution using the HF/6-31 + G (3df, 2p) level of theory. A computation using the B3LYP/6-31 + G (3df, 2p) level of theory and the Polarized Continuum Model (PCM) and by treating water as a solvent is also conducted for each activation process. It is found that, energetically, the most favor activation process for Ser70 in the acylation mechanism is the one where a proton transfer is mediated by the catalytic water and the catalytic residues Glu166 and Ser70. This agrees with those observed in an ultrahigh resolution X-ray structure and a QM/MM theoretical study published recently on the same acylation process.     

Probing the catalytically essential residues of the alpha-L-fucosidase from the hyperthermophilic archaeon Sulfolobus solfataricus

Retaining glycosidases promote the hydrolysis of the substrate by following a double-displacement mechanism involving a covalent intermediate. The catalytic residues are a general acid/base catalyst and the nucleophile. Experimental identification of these residues in a specific glycosidase allows for the assigning of the corresponding residues in all of the other enzymes belonging to the same family. By means of sequence alignment, mutagenesis, and detailed kinetic studies of the alpha-fucosidase from Sulfolobus solfataricus (Ssalpha-fuc) (family 29), we show here that the residues, invariant in this family, have the function inferred from the analysis of the 3D structure of the enzyme from Thermotoga maritima (Tmalpha-fuc). These include in Ssalpha-fuc the substrate-binding residues His46 and His123 and the nucleophile of the reaction, previously described. The acid/base catalyst could be assigned less easily. The k(cat) of the Ssalpha-fucGlu292Gly mutant, corresponding to the acid/base catalyst of Tmalpha-fuc, is reduced by 154-fold but could not be chemically rescued. Instead, the Ssalpha-fucGlu58Gly mutant revealed a 4000-fold reduction of k(cat)/K(M) if compared to the wild-type and showed the rescue of the k(cat) by sodium azide at wild-type levels. Thus, our data suggest that a catalytic triad, namely, Glu58, Glu292, and Asp242, is involved in catalysis. Glu58 and Glu292 cooperate in the role of acid/base catalyst, while Asp242 is the nucleophile of the reaction. Our data suggest that in glycosidase family 29 alpha-fucosidases promoting the retaining mechanism with slightly different catalytic machineries coexist.     

Conserved domains of glycosyltransferases.

Glycosyltransferases catalyze the synthesis of glycoconjugates by transferring a properly activated sugar residue to an appropriate acceptor molecule or aglycone for chain initiation and elongation. The acceptor can be a lipid, a protein, a heterocyclic compound, or another carbohydrate residue. A catalytic reaction is believed to involve the recognition of both the donor and acceptor by suitable domains, as well as the catalytic site of the enzyme. To elucidate the structural requirements for substrate recognition and catalytic reactions of glycosyltransferases, we have searched the databases for homologous sequences, identified conserved amino acid residues, and proposed potential domain motifs for these enzymes. Depending on the configuration of the anomeric functional group of the glycosyl donor molecule and of the resulting glycoconjugate, all known glycosyltransferases can be divided into two major types: retaining glycosyltransferases, which transfer sugar residue with the retention of anomeric configuration, and inverting glycosyltransferases, which transfer sugar residue with the inversion of anomeric configuration. One conserved domain of the inverting glycosyltransferases identified in the database is responsible for the recognition of a pyrimidine nucleotide, which is either the UDP or the TDP portion of a donor sugar-nucleotide molecule. This domain is termed "Nucleotide Recognition Domain 1 beta," or NRD1 beta, since the type of nucleotide is the only common structure among the sugar donors and acceptors. NRD1 beta is present in 140 glycosyltransferases. The central portion of the NRD1 beta domain is very similar to the domain that is present in one family of retaining glycosyltransferases. This family is termed NRD1 alpha to designate the similarity and stereochemistry of sugar transfer, and it consists of 77 glycosyltransferases identified thus far. In the central portion there is a homologous region for these two families and this region probably has a catalytic function. A third conserved domain is found exclusively in membrane-bound glycosyltransferases and is termed NRD2; this domain is present in 98 glycosyltransferases. All three identified NRDs are present in archaebacterial, eubacterial, viral, and eukaryotic glycosyltransferases. The present article presents the alignment of conserved NRD domains and also presents a brief overview of the analyzed glycosyltransferases which comprise about 65% of all known sugar-nucleotide dependent (Leloir-type) and putative glycosyltransferases in different databases. A potential mechanism for the catalytic reaction is also proposed. This proposed mechanism should facilitate the design of experiments to elucidate the regulatory mechanisms of glycosylation reactions. Amino acid sequence information within the conserved domain may be utilized to design degenerate primers for identifying DNA encoding new glycosyltransferases.     

Site-directed mutagenesis establishes aspartic acids-227 and -342 as essential for enzyme activity in an isomalto-dextranase from Arthrobacter globiformis

Isomalto-dextranase, from Arthrobacter globiformis T6, is a member of the glycoside hydrolase family 27. However, the alignments of the whole amino acid sequence are distinct from other members of this family. The enzymes cleave the glycosidic bond of the substrate in two different manners: either retaining or inverting the anomeric configuration. We believe that a retaining enzyme is involved in a two-step, double-displacement mechanism utilizing active site carboxylic acids as the nucleophile and general acid/base catalysts in the hydrolytic reaction. The critical amino acid residues at the isomalto-dextranase active site that catalyzes the hydrolysis reaction of dextran have been identified and the roles of nine amino acid residues (D107, D163, D227, D295, D340, D342, D373, D396, and E420) in the isomalto-dextranase from A. globiformis analyzed by site-directed mutagenesis. Of 15 mutant enzymes that were prepared, eight had reduced activities for dextran hydrolysis. Aspartic acids-227 and -342, which are part of the apparent catalytic dyad, were essential for hydrolase activity toward dextran.     

Density functional computations of enantioselective alkynylation of aldehyde catalyzed by chiral zinc(II)-complexes

The enantioselective alkynylation of aldehyde catalyzed by chiral zinc(II)-complexes was studied by means of the density functional theory (DFT). All the structures were optimized completely at the B3LYP/6-31G(d,p) level. To obtain more exact energies, single-point energy calculations at B3LYP/6-31+G(d,p) level were carried out on the B3LYP/6-31G(d,p) geometries. As shown, this enantioselective alkynylation was endothermic. The chirality-determining step for the alkynylation was the formation of the catalyst–ethanol complexes and the transition states for this step involved a six-membered ring. The dominant products predicted theoretically were of (R)-chirality, in good agreement with experiment.      

Computational studies on polynitropurines as potential high energy density materials

As part of a search for high energy density materials (HEDMs), a series of purine derivatives with nitro groups were designed computationally. The relationship between the structures and the performances of these polynitropurines was studied. Density functional theory (DFT) at the B3LYP/6-311G** level was employed to evaluate the heats of formation (HOFs) of the polynitropurines by designing an isodesmic reaction method. Results indicated that the HOFs were influenced by the number and positions of substituent groups. Detonation properties were evaluated using the Kamlet–Jacobs equations, based on the theoretical densities and heats of formation of the polynitropurines. The relative stabilities of the polynitropurines were studied via the pyrolysis mechanism and the UB3LYP/6-311G** method. Homolysis of the ring–NO₂ bond is predicted to be the initial step in the thermal decomposition of these purine derivatives. Considering their detonation properties and relative stabilities, the tetranitropurine (D1) derivatives may be regarded as potential candidates for practical HEDCs. These results may provide useful information for further investigations.     

Site-directed mutagenesis of glutamate 317 of bovine alpha-1,3Galactosyltransferase and its effect on enzyme activity: implications for reaction mechanism

Bovine alpha1,3galactosyltransferase (alpha1,3GalT) transfers galactose from UDP-alpha-galactose to terminal beta-linked galactosyl residues with retention of configuration of the incoming galactose residue. The epitope synthesized has been shown to be critical for xenotransplantation. According to a proposed double-displacement reaction mechanism, glutamate-317 (E317) is thought to be the catalytic nucleophile. The proposed catalytic role of E317 involves an initial nucleophilic attack with inversion of configuration and formation of a covalent sugar-enzyme intermediate between E317 and galactose from the donor substrate, followed by a second nucleophilic attack performed by the acceptor substrate with a second inversion of configuration. To determine whether E317 of alpha1,3GalT is critical for enzyme activity, site-directed mutagenesis was used to substitute alanine, aspartic acid, cysteine and histidine for E317. If the proposed reaction mechanism for the alpha1,3GalT enzyme is correct, E317D and E317H would produce active enzymes since they can act as nucleophiles. The non-conservative mutation E317A and conservative mutation E317C are predicted to produce inactive or very low activity enzymes since the E317A mutant cannot engage in a nucleophilic attack, and the E317C mutant would trap the galactose residue. The results obtained demonstrate that E317D and E317H mutants retained activity, albeit significantly less than the wild-type enzyme. Additionally, both E317A and E317C mutant also retained enzyme activity, suggesting that E317 is not the catalytic nucleophile proposed in the double-displacement mechanism. Therefore, either a different amino acid may act as the catalytic nucleophile or the reaction must proceed by a different mechanism.     

Theoretical study of the decomposition mechanisms and kinetics of the ingredients RDX in composition B

RDX as a component in composition B (TNT + RDX) was first studied by us on its mechanism and kinetics of decomposition reactions in this paper. We have pointed out three possible pathways and found a new low-energy process of its decomposition. The N-N bond cleavage in composition B has higher dissociation energies than the monomer, but it is also the initial step. The optimized structures and the frequencies of all the stationary points were calculated at the B3LYP/6-31G(d) level. The minimum-energy paths were obtained by using the intrinsic reaction coordinate (IRC) theory, and the reaction potential energy curve was corrected with zero-point energy. Finally, the rate constants were calculated in a wide temperature region from 200 to 2500 K using TST, TST/Eckart theories. The obtained results also indicate that the tunneling effects are remarkable at low temperature (200 K 相似文献