A QM/MM-based computational investigation on the catalytic mechanism of saccharopine reductase

Saccharopine reductase from Magnaporthe grisea, an NADPH-containing enzyme in the α-aminoadipate pathway, catalyses the formation of saccharopine, a precursor to L-lysine, from the substrates glutamate and α-aminoadipate-δ-semialdehyde. Its catalytic mechanism has been investigated using quantum mechanics/molecular mechanics (QM/MM) ONIOM-based approaches. In particular, the overall catalytic pathway has been elucidated and the effects of electron correlation and the anisotropic polar protein environment have been examined via the use of the ONIOM(HF/6-31G(d):AMBER94) and ONIOM(MP2/6-31G(d)//HF/6-31G(d):AMBER94) methods within the mechanical embedding formulism and ONIOM(MP2/6-31G(d)//HF/6-31G(d):AMBER94) and ONIOM(MP2/6-311G(d,p)//HF/6-31G(d):AMBER94) within the electronic embedding formulism. The results of the present study suggest that saccharopine reductase utilises a substrate-assisted catalytic pathway in which acid/base groups within the cosubstrates themselves facilitate the mechanistically required proton transfers. Thus, the enzyme appears to act most likely by binding the three required reactant molecules glutamate, α-aminoadipate-δ-semialdehyde and NADPH in a manner and polar environment conducive to reaction.     

Host-guest inclusion complex between β-cyclodextrin and paeonol: A theoretical approach

Host-guest interactions of β-cyclodextrin (β-CD) with paeonol (PAE) were simulated using semi-empirical PM3 and both ONIOM2 [(B3LYP/6-31G*:PM3), (HF/6-31G*:PM3)] methods. The results obtained with PM3 method clearly indicate that the complexes formed are energetically favored with or without solvent, the model 1 (PAE entering into the cavity of β-CD from its wide side by OCH3 group) is found more favored than the model 2 (PAE entering into the cavity of β-CD from its wide side by COCH3 group). Finally, natural bonding orbital (NBO) analysis was performed based on ONIOM2 optimized complexes to quantify the donor–acceptor interactions between PAE and β-CD.     

A computational study of intramolecular proton transfer in gaseous protonated glycine

The minimum energy paths for intramolecular proton transfer between the amino nitrogen and carbonyl oxygen atoms in gaseous protonated glycine were estimated at the Hartree-Fock (HF) and second-order M?ller-Plesset Perturbation (MP2) levels of theory. Potential energy profiles and their associated reactant, transition state, and product species calculated at the MP2/6-31G* level were shown to differ significantly from those obtained at the HF/6-31G* level. Effects of electron correlation and basis functions on the calculated geometries and energies of relevant species were examined at the HF, MP2, MP4, CCSD, and B3LYP levels using the 6-31G*, 6-31G**, 6-31+G**, 6-311+G**, 6-31+G(2d,2p), 6-311+G(3df,2p), cc-pVDZ, aug-cc-pVDZ, and cc-pVTZ basis sets. The HF and MP2 optimized levels with the 6-31G*, 6-31G**, 6-31+G**, and 6-311+G** bases were used to calculate the thermodynamic and kinetic properties of the proton transfer reaction at 298.15 K and 1 atm, which include enthalpy, entropy, Gibbs free energy, equilibrium constant, potential energy barriers, tunneling transmission coefficients, and rate constants. Results indicate that the proton in a carbonyl O-protonated glycine undergoes a rapid migration to the amino nitrogen atom, while the reverse process is extremely unfavorable. The objective of this work is to develop practical theoretical procedures for studying proton transfer reactions in amino acids and peptides and to assemble physical data from these model calculations for future references.     

Calculated structures and relative stabilities of furoxan,some 1,2-dinitrosoethylenes and other isomers

The structures and relative stabilities of furoxan and some of its isomers, e.g., the 1,2-dinitrosoethylenes, have been determined by means of ab initio Hartee–Fock and Møller–Plesset calculations. Geometries were optimized at the HF/3-21G, HF/6-31G* and MP2/6-31G* levels, and subsequently used for computing MP2/6-31G*, MP3/6-31G*, and MP4/6-31G* energies. The results are markedly affected by the inclusion of electronic correlation, which renders three of the isomers unstable. It also emphasizes the importance of a zwitterionic contribution to the structure of furoxan, which promotes ring-opening through a cis 1,2-dinitrosoethylene intermediate/transition state that has an MP4/6-31G*//MP2/6-31G* energy that is 31.6 kcal/mol above furoxan.     

异硫氰酸与甲亚胺环加成反应机理的理论研究

对异硫氰酸与甲亚胺形成四元或六元环产物的环加成反应进行了理论研究,结果表明,一分子异硫氰酸与一分子甲亚胺形成四元环(1:1)产物P1的反应(1)为经过一个两性离子中间体的分步反应,其中第二步为速控步骤,其活化热垒为107.86kJ/mol.此外,反应(1)的中间体还可与另一甲亚胺或异硫氰酸分子继续反应形成两个不同的六元环(1:2或2:1)产物P2或P3;这两个反应均为协同反应,其活化势垒分别为15.88和21.82kJ/mol.这些结果与当异硫氰酸酯与亚胺发生环加成反应时,只有类似于P2和P3的两种六元环产物生成的实验事实一致。     

双自由基CF2与O3的反应机理

主要用作致冷剂和发泡剂的氯氟烃（CFCs）是破坏臭氧层的主要物质之一．对氯氟烃类化合物及其降解产物（包括光解、光氧化、化学反应产物等）在大气中行为问题的研究是大气化学研究的重要内容．前人[1-3]从理论和实验两方面研究了自由基与臭氧的反应机制，但是氯氟烃光解过程中     

Comparative theoretical study of transition structures,barrier heights,and reaction energies for the intramolecular tautomerization in acetaldehyde/vinyl alcohol and acetaldimine/vinylamine systems

The transition structures associated with the possible intramolecular tautomerization for acetaldehyde/vinyl alcohol and acetaldimine/vinylamine systems as models of keto/enol and imine/enamine interconversion processes, respectively, were characterized. The relative stabilities of the tautomers and the associated barrier heights were calculated. Ab initio analytical gradients and second derivatives at the HF level of theory and 3-21G, 6-31G, 6-31G**, 6-31++G**, and 6-311++G** basis-set, DFT (BP86/6-311++G** and BLYP/6-311++G**), and semiempirical (AM1 and PM3) procedures were used to identify the stationary points. Correlation effects were estimated using the perturbational approach at MP2/6-31G**, MP2/6-311++G**, and MP2/6-311++G (3df,2p) levels. The geometry, electronic structure, harmonic vibrational frequencies, and transition vector associated with the transition structures as well as the relative stabilities of different isomers and barrier heights were analyzed. The dependence of these properties upon theoretical methods is analyzed and discussed. The transition structures are four-membered rings and the corresponding transition vectors are associated to collective fluctuations. The 1,3 intramolecular hydrogen migration is much more advanced than are the hybridization changes on donor and acceptor centers at the transition structure. The corresponding barrier heights can be related to the change of bond orders and acid/base properties of these centers. A comparison of the results obtained with different methods renders that the nature of the transition structure seems to be a rather robust entity. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 9–24, 1998     

Quantum chemical calculations based on ONIOM and the DFT methods in the inclusion complex: doxycycline/2-O-Me-β-cyclodextrin

A computational study of inclusion complexes of 2-methyl-βCD with Doxycycline tautomeric (enol and keto form) has been performed with several combinations of ONIOM hybrid calculations. The reliability of the ONIOM2 calculations at the integrated level, ONIOM2 (M05-2X/6-31G(d): M05-2X/3-21G*), ONIOM2 (M05-2X/6-31G(d):HF/3-21G*), ONIOM2 (B3LYP/6-31G(d):HF/3-21G*), ONIOM2 (B3LYP/6-31G(d):B3LYP/3-21G*) and ONIOM2 (B3PW91/6-31G(d):B3PW91/3-21G*) was examined. Their complexation, binding, deformation and stabilization energies, and geometrical data were compared with those of the target geometry structure optimized at the M05-2X/6-31G(d) level of theory. Mixed combinations ONIOM2 (M05-2X 6-31G(d):HF 3-21G*) and ONIOM2 (B3LYP 6-31G(d):HF 3-21G*) reproduces nearly the target geometry structure and provides realistic energetic results at a relatively low computational cost.     

适用于液晶大分子的量子化学计算方法以及六烷氧基苯并菲取代的1,3,5-三酰胺苯液晶的研究

盘状液晶分子之间的相互作用决定液晶的性质. 为了选择合适的计算方法, 以便用量子化学研究液晶大分子, 设计了对位取代苯和三酰胺苯作为模型分子, 用高水平的ONIOM [MP2/6-31G*(0.25):HF/6-31G(d,p)]计算提供了与晶体结构资料相符合的较准确的二体相互作用结果. 然后与各种较低级别的计算进行比较, 说明ONIOM (HF/STO-3G: AM1:UFF)水平比较合适. 盘重叠部分的相互作用主要是色散作用, 用UFF力场处理, 氢键主要是静电作用, 用HF/STO-3G处理, 其余部分用AM1处理. 通过对苯取代的1,3,5-三酰胺苯的双分子构型优化, 并与晶体结构数据进行比较, 进一步说明ONIOM (HF/STO-3G:AM1:UFF)水平计算的适用性. 在此基础上, 对六烷氧基苯并菲取代的1,3,5-三酰胺苯的双分子构型进行了优化, 为解释它所形成的液晶具有较高的电荷传输能力提供了有用的信息.     

Inclusion complexes of ortho-anisidine and β-cyclodextrin: A quantum mechanical calculation

The structural aspects for the complexation of ortho-anisidine (O-AN)/β-cyclodextrin were explored by using PM6, density function theory B3LYP/6-31G*, M05-2X/6-31G*, B3PW91/6-31G*, MPW1PW91/6-31G*, HF/6-31G* methods and several combinations of ONIOM2 hybrid calculations. Calculations were performed upon the inclusion complexation of β-cyclodextrin (β-CD) with neutral (O-AN1) and cationic (O-AN2) species of ortho-anisidine. The obtained results with PM6 method clearly indicate that the formed complexes are energetically favored, the complex of O-AN2/β-CD in B orientation is significantly more favorable than the others energetically. The structures show the presence of several intermolecular hydrogen bond interactions that were studied on the basis of natural bonding orbital (NBO) analysis, employed to quantify the donor–acceptor interactions between ortho-anisidine and β-CD.     

Rearrangement of azirine intermediates to nitriles: Theoretical study of cleavage of 3,4-dihydro-1aH-azirine[2,3-c]pyrrol-2-one to cyanoketene–formaldimine complex

A detailed investigation of the reaction path for the thermal rearrangement of 3,4-dihydro-1aH-azirine[2,3-c]pyrrol-2-one to yield a cyanoketene–formaldimine complex is carried out at the MP2/6-31G* and B3LYP/6-31G* levels of theory. The ring opening of the five-membered pyrrolinone ring and the formation of the nitrile group takes place in a concerted manner, presenting a significant strain energy release and allowing for an electronic stabilization by coarctate conjugation of the transition structure (TS). These two factors make possible a moderate energy barrier. Although the structural features B3LYP/6-31G* theoretical levels, it is found that the MP2 energy barrier (28.8) CCSD(T)/6-31G*//MP2/6-31G* value (17.1 kcal/mol). The complex electronic rearrangement can be rationalized using the theory of coarctate transition structures developed by Herges as the evolution of an azirine structure without referring to a hypothetical vinyl nitrene intermediate. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 912–922, 1998     

Comparison of levels of electronic structure theory in direct dynamics simulations of C2H5F --> HF + C2H4 product energy partitioning

Direct dynamics simulations at the MP2/6-311++G** level of theory were performed to study C(2)H(5)F --> HF + C(2)H(4) product energy partitioning. The simulation results are compared with experiment and a previous MP2/6-31G* simulation. The current simulation with the larger basis set releases more energy to HF vibration and less to HF + C(2)H(4) relative translation as compared to the previous simulation with the 6-31G* basis set. The HF rotation and vibration energy distributions determined from the current simulation are in overall very good agreement with previous experimental studies of C(2)H(5)F dissociation by chemical activation and IRMPA. A comparison of the simulations with experiments suggests there may be important mass effects for energy partitioning in HX elimination from haloalkanes. The transition state (TS) structures and energies calculated with MP2 and the 6-31G* and 6-311++G** basis sets are compared with those calculated using CCD, CCSD, CCSD(T), and the 6-311++G** basis set.     

乙烯醇钠的从头算研究(英文)

在ＲＨＦ／３－２１Ｇ和ＲＨＦ／６－３１＋Ｇ水平上优化了乙烯醇钠的的平衡构型及过渡态，比较了不同构型的稳定性。在ＲＨＦ／６－３１＋Ｇ水平上探讨了气相中由乙醛和氢氧化钠生成乙烯醇钠的机理。并在ＭＰ２（ｆｕｌ）／６－３１＋Ｇ／／ＭＰ２（ｆｕｌ）／６－３１＋Ｇ水平下计算了反应热     

Investigation of ethylene dimerization over faujasite zeolite by the ONIOM method.

Ethylene dimerization was investigated by using an 84T cluster of faujasite zeolite modeled by the ONIOM3(MP2/6-311++G(d,p):HF/6-31G(d):UFF) method. Concerted and stepwise mechanisms were evaluated. In the stepwise mechanism, the reaction proceeds by protonation of ethylene to form the surface ethoxide and then C--C bond formation between the ethoxide and the second ethylene molecule to give the butoxide product. The first step is rate-determining and has an activation barrier of 30.06 kcal mol(-1). The ethoxide intermediate is rather reactive and readily reacts with another ethylene molecule with a smaller activation energy of 28.87 kcal mol(-1). In the concerted mechanism, the reaction occurs in one step of simultaneous protonation and C--C bond formation. The activation barrier is calculated to be 38.08 kcal mol(-1). Therefore, the stepwise mechanism should dominate in ethylene dimerization.     

Characteristic negative ion fragmentations of deprotonated peptides containing post-translational modifications: mono-phosphorylated Ser, Thr and Tyr. A joint experimental and theoretical study

Peptides and proteins may contain post-translationally modified phosphorylated amino acid residues, in particular phosphorylated serine (pSer), threonine (pThr) and tyrosine (pTyr). Following earlier work by Lehmann et al., the [M-H]- anions of peptides containing pSer and pThr functionality show loss of the elements of H3PO4. This process, illustrated for Ser (and using a model system), is CH3CONH-C(CH2OPO3H2)CONHCH(3) --> [CH3CONHC(==CH2)CONHCH3 (-OPO3H2)] (a) --> [CH3CONHC(==CH2)CONHCH3-H]- + H3PO4, a process endothermic by 83 kJ mol(-1) at the MP2/6-31++G(d,p)//HF/6-31++G(d,p) level of theory. In addition, intermediate (a) may decompose to yield CH3CONHC(==CH2)CONHCH3 + H2PO4 - in a process exothermic by 3 kJ mol(-1). The barrier to the transition state for these two processes is 49 kJ mol(-1). Characteristic cleavages of pSer and pThr are more energetically favourable than the negative ion backbone cleavages of peptides described previously. In contrast, loss of HPO3 from [M-H]- is characteristic of pTyr. The cleavage [NH2CH(CH2-C6H4-OPO3H-)CO2H] --> [NH2C(CH2-C6H4-O-)CO2H (HPO3)] (b) --> NH2CH(CH2-C6H4-O-)CO2H + HPO3 is endothermic by 318 kJ mol(-1) at the HF/6-31+G(d)//AM1 level of theory. In addition, intermediate (b) also yields NH2CH(CH2-C6H4-OH)CO2H + PO3 - (reaction endothermic by 137 kJ mol(-1)). The two negative ion cleavages of pTyr have a barrier to the transition state of 198 kJ mol(-1) (at the HF/6-31+G(d)//AM1 level of theory) comparable with those already reported for negative ion backbone cleavages.     

Benzooxirene. Ab initio calculations

Ab initio calculations have been performed on benzooxirene, the corresponding oxo carbene (“ketocarbene”), and the transition state linking the two. At the highest level used, QCISD(T)/6-31G^*//MP2(FULL)/6-1G^* with MP2(FULL)/ 6-31G^* zero point energy corrections, the relative energies of the oxirene, the transition state and the carbene are 0, 24.6, and −17.8 kJ mol⁻¹. Correlation energy effects are very important in this system: at the QCISD(T) level the oxirene lies above the carbene, as at the MP4 and HF levels, but at the MP2 level the ordering is reversed. Benzooxirene is probably slightly nonplanar: the HF/6-31G^* geometry is C_2v but the MP2(Fermi contact)/6-31G^* geometry is C_s with a 6-/3-ring coplanarity deviation of about 6.9 °, although in the MP2(FULL)/6-31G^* geometry this is reduced to about 3.1 °.     

Conformational study of the structure of 12-crown-4-alkali metal cation complexes

A conformational search was performed for the 12-crown-4 (12c4)-alkali metal cation complexes using two different methods, one of them is the CONFLEX method, whereby eight conformations were predicted. Computations were performed for the eight predicted conformations at the HF/6-31+G*, MP2/6-31+G*//HF/6-31+G*, B3LYP/6-31+G*, MP2/6-31+G*//B3LYP/6-31+G*, and MP2/6-31+G* levels. The calculated energies predict a C4 conformation for the 12c4-Na+, -K+, -Rb+, and -Cs+ complexes and a C(s) conformation for the 12c4-Li+ complex to be the lowest energy conformations. For most of the conformations considered, the relative energies, with respect to the C4 conformation, at the MP2/6-31+G*//B3LYP/6-31+G* are overestimated, compared to those at the MP2/6-31+G* level, the highest level of theory considerd in this report, by 0.2 kcal/mol. Larger relative energy differences are attributed to larger differences between the B3LYP and MP2 optimized geomtries. Binding enthalpies (BEs) were calculated at the above-mentioned levels for the eight conformations. The agreement between the calculated and experimental BEs is discussed.     

A theoretical study of the effect of basis sets on stationary structures for the addition of carbon dioxide to methylamine: A relation among geometries,energy status,and electronic structures

The carbon dioxide (CO₂) to methylamine (CH₃NH₂) addition reaction is investigated using an ab initio MO method with the aim of obtaining an insight into the relation between transition-state structure and reactivity. The complete potential energy hypersurface was explored and the stationary points representing reactants, transition state, and product were localized. To establish the relevance of different basis sets, calculations were carried out at different levels of theory: STO -3G , 3-21G , 4-21G , 4-31G , 6-31G , 6-311G , 6-31G *, and MP 2/6-31G *//HF /6-31G *. A reaction analysis done by correlation of bond-order indices using the More O'Ferrall-Jencks diagram reveals the reaction mechanism as an asynchronous process. The transition state can be described as a four-membered ring. The lone nitrogen electron pair of the amine appears as the driving force of the reaction, where a positive hydrogen rather than a radical is more likely to be transferred. © 1993 John Wiley & Sons, Inc.     

低硝酸盐离子N2O22-和亚硝酸HONO的异构体及其异构化反应的量化研究

采用Gaussian-98程序进行，在HF／6—31G(d)，B3LYP／6—31G(d)和MP2／6—31G(d)水平下优化分子结构并寻找过渡态，对于MP2／6—31G(d)结果在QCISD(T，E4T)，MP4／6—311 G(d，p)，MP4／6—311 G(2df，p)水平下重新计算能级．并用内禀反应坐标(IRC)法研究了N2O2^2-和亚硝酸HONO的异构化反应机理。     

Theoretical studies on the conformations and properties of phenylazonaphthalenes

The conformations of phenylazo‐2‐naphthalene (I) and phenylazo‐1‐naphthalene (II) have been studied using ab initio methods and density functional theory. The rotational potential energy surfaces about the C N bonds were calculated for both the trans and the cis forms at the PM3 and HF/STO‐3G levels. The PM3 method was found to be incapable for locating the energy minima correctly. The geometries of rotamers obtained from the HF/STO‐3G surface were fully optimized at the HF/6‐31G* and B3LYP/6‐31G* levels. Single‐point MP2/6‐31G* calculations have also been carried out at the HF/6‐31G* geometries. The vibrational spectra, standard thermodynamic functions, heats of formation, and equilibrium molar fractions have been obtained. According to the calculated results, I is comparatively more stable than II. The trans forms have much lower energies than the cis forms, implying that I and II mainly exist in the trans planar forms. The energy gap (ΔE) between the frontier molecular orbitals (MOs) of the cis forms are about 0.7 eV higher than those of the trans ones, indicating that change in geometry from trans to cis form will cause a hypsochromic effect. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 79: 25–33, 2000