1-(4-硝基苯基)-3-(5,6-二甲基-1,2,4三氮唑)-三氮烯质子迁移的理论研究

采用密度泛函B3LYP在6-31G*基组下,对有机显色剂1-(4-硝基苯基)-3-(5,6-二甲基-1,2,4三氮唑)-三氮烯(NPDMTT)的各种可能结构进行质子迁移的3种可能途径:(a)分子内质子迁移,(b)水助质子迁移,(c)甲醇助质子转移的计算,得到了各种途径异构体的相对能,获得了它们的互变异构过程的活化能、活化吉布斯自由能和质子转移反应的速率常数等性质。计算结果表明,分子内质子转移形成的各种异构体相对能量较大,当水分子或甲醇分子参与反应时,异构体的相对能量明显减小,但无论是孤立分子、一水合物还是一甲醇合物,其最稳定的异构体都相同,均为A2。溶剂化效应对异构化能垒的影响较大。最稳定的异构体分子内质子转移在N11和N13间转移的速控步骤的活化能为130.9 kJ.mol-1,反应速度常数为2.172×10-11s-1;当水分子参与反应以双质子转移机理异构化时,活化能显著降低,有利于三氮异构化,其中异构体质子在N11和N13间转移的速控步骤的活化能为22.55 kJ.mol-1,反应速度常数为3.617×107s-1;当醇分子参与反应以双质子转移机理异构时,活化能减小得更多,其中异构体质子在N11和N13间转移的速控步骤的活化能为2.384 kJ.mol-1,反应速度常数为9.032×1011s-1。计算结果还表明,氢键作用在增大NPDMTT一水合物和NPDMTT一甲醇合物相对稳定性、降低质子转移异构化反应活化能等方面起着重要作用。     

HClO对贻贝粘附单元DOPA氧化的理论研究

采用量子化学密度泛函理论(DFT)研究贻贝粘附单元DOPA(3,4-二羟基苯丙氨酸)的结构与性质,得分子的几何构型、原子电荷分布、反应活性及热力学等参数,表明:DOPA苯环易与HClO(次氯酸)发生亲电取代反应(1),生成3-氯4,5-二羟基苯丙氨酸,阻碍生成贻贝内超强粘附单元DOPA二联体,降低粘附蛋白间粘性;DOPA侧链易与HClO发生亲电亲核反应(2),促使DOPA侧链的断裂,降低粘附蛋白内粘性;在相同温度下,反应(1)和反应(2)的△G<0,且△G(1)<△G(2),反应(1)较易发生.     

分子筛催化甲苯歧化SE1反应机理的分子模拟研究

应用分子模拟半经验量子力学Mopac 6．0-AM1近似计算方法分析了分子筛催化甲苯歧化反应的S_E~1反应历程,确定了反应历程中的反应态、过渡态和产物态,得到了反应活化能和反应热等相关信息,对内禀反应坐标的计算进一步验证了反应过程中的能量变化。计算结果表明,分子筛催化的甲苯歧化反应沿S_E~1反应历程可通过两步基元反应完成;质子由分子筛向甲苯分子转移的过程为反应的快速步骤,其活化能达到428．54 kJ·mol~(-1),需要在高温下进行;甲苯歧化总反应的热效应很小,与实验数据相吻合。     

分子筛合成过程中三聚铝硅酸盐生成机理的DFT研究

采用密度泛函理论研究了在分子筛合成的碱性环境中铝硅酸盐三聚体的生成机理.在B3LYP/6-31 G(d,p)计算水平上对反应物、过渡态和产物分别进行了几何结构的全优化和频率计算,通过内禀反应坐标的方法验证了反应路径,并计算了反应的活化能.采用COSMO-RS模型考虑了溶剂效应.结果表明,二聚硅酸分子SiOSi(OH)6和单体铝酸根离子Al(OH)**的缩聚反应可以按照协同方式进行,SiO-Al桥键的形成与水分子的脱除同时发生,最终生成铝硅酸盐的三聚体.计算得到的缩聚反应活化能为80.1 kJ/mol.     

Combined quantum mechanical and molecular mechanical reaction pathway calculation for aromatic hydroxylation by p-hydroxybenzoate-3-hydroxylase

The reaction pathway for the aromatic 3-hydroxylation of p-hydroxybenzoate by the reactive C4a-hydroperoxyflavin cofactor intermediate in p-hydroxybenzoate hydroxylase (PHBH) has been investigated by a combined quantum mechanical and molecular mechanical (QM/MM) method. A structural model for the C4a-hydroperoxyflavin intermediate in the PHBH reaction cycle was built on the basis of the crystal structure coordinates of the enzyme-substrate complex. A reaction pathway for the subsequent hydroxylation step was calculated by imposing a reaction coordinate that involves cleavage of the peroxide oxygen-oxygen bond and formation of the carbon-oxygen bond between the C3 atom of the substrate and the distal oxygen of the peroxide moiety of the cofactor. The geometric changes and the Mulliken charge distributions along the calculated reaction pathway are in line with an electrophilic aromatic substitution type of mechanism. The energy barrier of the calculated reaction is considerably lower when the substrate hydroxyl moiety is deprotonated, in comparison with the barrier found with a protonated hydroxyl moiety. This effect of the protonation state of the substrate on the calculated energy barrier supports experimental observations that deprotonation is required for hydroxylation of the substrate. A notable event in the calculated reaction pathway is a lengthening of the peroxide oxygen-oxygen bond at an intermediate stage. Further analysis of the reaction pathway indicates that this oxygen-oxygen bond elongation is accompanied by an increase in electrophilic reactivity on the distal oxygen of the peroxide moiety, which may assist the C-O bond formation in the reaction of the C4a-hydroperoxyflavin intermediate with the substrate. Analysis of the effect of individual active site residues on the reaction reveals a specific transition state stabilization by the backbone carbonyl moiety of Pro293. The crystal water 717 appears to drive the hydroxylation step through a stabilizing hydrogen bond interaction to the proximal oxygen of the C4a-hydroperoxyflavin intermediate, which increases in strength as the hydroperoxyflavin cofactor converts to the anionic (deprotonated) hydroxyflavin.     

Theoretical Study on Free Fatty Acid Elimination Mechanism for Waste Cooking Oils to Biodiesel over Acid Catalyst

A theoretical investigation on the esterification mechanism of free fatty acid (FFA) in waste cooking oils (WCOs) has been carried out using DMol³ module based on the density functional theory (DFT). Three potential pathways of FFA esterification reaction are designed to achieve the formation of fatty acid methyl ester (FAME), and calculated results show that the energy barrier can be efficiently reduced from 88.597 kcal/mol to 15.318 kcal/mol by acid catalyst. The molar enthalpy changes (Δ_rH_m°) of designed pathways are negative, indicating that FFA esterification reaction is an exothermic process. The obtained favorable energy pathway is: H⁺ firstly activates FFA, then the intermediate combines with methanol to form a tetrahedral structure, and finally, producing FAME after removing a water molecule. The rate-determining step is the combination of the activated FFA with methanol, and the activation energy is about 11.513 kcal/mol at 298.15 K. Our results should provide basic and reliable theoretical data for further understanding the elimination mechanism of FFA over acid catalyst in the conversion of WCOs to biodiesel products.     

6-巯基嘌呤质子转移异构化的量子化学计算研究

采用密度泛函B3LYP方法,在6-311G(d,p)基组水平上对6-巯基嘌呤质子转移引起的硫酮式与硫醇式互变异构反应机理进行了计算研究,获得了互变异构过程的反应焓、活化能、活化吉布斯自由能和质子转移反应的速率常数等参数。计算结果表明,6-巯基嘌呤无论是孤立分子还是一水合物,其硫酮式TP2是最稳定的异构体。计算结果同已有实验结果相符。由硫酮式通过分子内质子转移向硫醇式异构化找到4条反应通道(P1,P2,P3,P4),各通道的活化能分别为114.0,133.9,128.0,95.1 kJ·mol~(-1)。当水分子参与反应以双质子转移机理异构化时,活化能垒显著降低,各通道的活化能依次降为51.2,63.0,70.5,42.8 kJ·mol~(-1),可见,水助催化有利于硫酮式向硫醇式转变。计算结果还表明,氢键的强弱对TP2一水合物的稳定性会有一定的影响。     

多氯代9,10-菲醌的热力学性质及稳定性的理论研究

采用密度泛函理论(DFT),在B3LYP/6-31G*水平上对135个多氯9,10-菲醌(PCPQ)系列化合物进行全优化和振动分析,得到各分子在298.15K,101.3kPa状态下的热力学参数。设计等键反应,计算了PCPQ系列化合物的标准生成焓(ΔfHθ)和标准生成自由能(ΔfGθ)。用SPSS13.0对这些热力学参数与氯原子的取代位置及取代数目(NPCS)进行多元线性回归,得出相关方程。结果表明:PCPQ系列化合物的Sθ、ΔfHθ和ΔfGθ与NPCS之间有很强的相关性。根据异构体标准生成自由能的相对大小,从理论上求得异构体的相对稳定性。同时,发现氯原子取代模式对扭角(C-C(=O)-C(=O)-C)有很大的影响。     

Understanding the molecular mechanism and regioselectivity in the synthesis of celecoxib via a domino reaction: A DFT study

The molecular mechanism and energetic of the domino reaction involved in the synthesis of celecoxib, a well-known anti-inflammatory drug, were theoretically studied at the DFT-B3LYP/6-31G^* level. The first reaction in this domino process, which is also the rate-determining step, is a complete regioselective [3 + 2] cycloaddition (32CA) reaction associated with the nucleophilic attack of C5 carbon atom of enamine 7 on the C3 carbon atom of nitrile imine 6, leading to cycloadduct 8. The second reaction is a rapid acid/base catalysed stepwise elimination reaction of the morpholine 9 from cycloadduct 8 affording celecoxibe 3. The results also show that neither molecular mechanism of reaction nor activation barriers are considerably affected by the inclusion of solvent. The calculated relative Gibbs free energies as well as local reactivity indices obtained using the calculated Parr functions explain the complete regioselective fashion provided by the 32CA reaction under consideration in excellent agreement with the experimental findings.     

Ab initio calculations of the effects of H+ and NH4+ on the initial decomposition of HMX

In this work, the effects of H(+) and NH(4)(+) on the initial decomposition of HMX were investigated on the basis of the B3P86/6-31G** and B3LYP/6-31G* calculations. Three initial decomposition pathways including the N-NO(2) bond fission, HONO elimination and C-N bond dissociation were considered for the complexes formed by HMX with H(+) (PHMX1 and PHMX2) or with NH(4)(+) (AHMX). We found that H(+) and NH(4)(+) did not evidently induce the HMX to trigger the N-NO(2) heterolysis because the energy barrier of N-NO(2) heterolysis was found to be higher than the bond dissociation energy of N-NO(2) homolytic cleavage. Meanwhile, the transition state barriers of the HONO elimination from the complexes were found to be similar to that from the isolated HMX, which means that the HONO elimination reaction of HMX was not affected by the H(+) and NH(4)(+). As for the ring-opening reaction of HMX due to the C-N bond dissociation, the calculated potential energy profile showed that the energy of the complex (AHMX) went uphill along the C-N bond length and no transition state existed on the curve. However, the transition state energy barriers of C-N bond dissociation were calculated to be only 5.0 kcal/mol and 5.5 kcal/mol for the PHMX1 and PHMX2 complexes, respectively, which were much lower than the C-N bond dissociation energy of isolated HMX. Moreover, among the three initial decomposition reactions, the C-N bond dissociation was also the most energetically favorable pathway for the PHMX1 and PHMX2. Our calculation results showed that the H(+) can significantly promote the initial thermal decomposition of C-N bond of HMX, which, however, is influenced by NH(4)(+) slightly.     

C6H6和SO3相互作用的理论研究（英文）

采用B3LYP和MP2方法在6-31G*、6-31+G*和6-311+G**基组下对C6H6…SO3复合物体系的4种可能结构进行自由优化,得3种。在考虑基组重叠误差校正基础上,得结合能,并用自然键轨道分析方法讨论其相互作用。结果表明,用B3LYP/6-31G*计算3种复合物的结合能分别为-17.75, -18.33, -18.80 kJ/mol,且C6H6和SO3结合时电子从苯环向SO3转移,形成电荷转移复合物,它们之间的作用包含π-p作用方式。     

亚硝酰氢与羟基自由基反应的密度泛函理论研究

用密度泛函理论研究了HNO OH反应机理。在(U)B3LYP/aug-cc-pVTZ水平上,优化了反应通道上各驻点(反应物、中间体、过渡态和产物)的几何构型,获得了零点能校正后的反应势能曲线。研究表明:根据进攻方式的不同,有3个反应通道,反应通道不同则产物不同。反应的主要产物是NO H_2O,次要产物是NH_2 O_2;主要产物和次要产物与反应物的总能量之差(经零点能校正后)分别为-133.42 kJ/mol和150.44 kJ/mol。生成主要产物时主要反应通道的活化能为102.11 kJ/mol。     

丙烯酸正丁酯聚合的Diels-Alder反应机理理论研究

本文从理论上对丙烯酸正丁酯(nBA)自引发聚合的Diels-Alder反应机理进行了研究.利用DFT方法在UB3LYP/6-31G*水平上对反应的最低能量路径进行了计算,各驻点能量分别采用MP2/6-311G**和B3LYP/6-311 G(3df;2p)进一步精确计算.结构表明:此Diels-Alder反应仅包括一种途径,即路径(I),另一条途径在热力学不支持.     

Elucidating the catalytic mechanism of β-secretase (BACE1): A quantum mechanics/molecular mechanics (QM/MM) approach

In this quantum mechanics/molecular mechanics (QM/MM) study, the mechanisms of the hydrolytic cleavage of the Met2-Asp3 and Leu2-Asp3 peptide bonds of the amyloid precursor protein (WT-substrate) and its Swedish mutant (SW) respectively catalyzed by β-secretase (BACE1) have been investigated by explicitly including the electrostatic and steric effects of the protein environment in the calculations. BACE1 catalyzes the rate-determining step in the generation of Alzheimer amyloid beta peptides and is widely acknowledged as a promising therapeutic target. The general acid-base mechanism followed by the enzyme proceeds through the following two steps: (1) formation of the gem-diol intermediate and (2) cleavage of the peptide bond. The formation of the gem-diol intermediate occurs with the barriers of 19.6 and 16.1 kcal/mol for the WT- and SW-substrate respectively. The QM/MM energetics predict that with the barriers of 21.9 and 17.2 kcal/mol for the WT- and SW-substrate respectively the cleavage of the peptide bond occurs in the rate-determining step. The computed barriers are in excellent agreement with the measured barrier of ∼18.0 kcal/mol for the SW-substrate and in line with the experimental observation that the cleavage of this substrate is sixty times more efficient than the WT-substrate.     

Ab initio QM/MM modelling of acetyl-CoA deprotonation in the enzyme citrate synthase

The first step of the reaction catalysed by the enzyme citrate synthase is studied here with high level combined quantum mechanical/molecular mechanical (QM/MM) methods (up to the MP2/6-31+G(d)//6-31G(d)/CHARMM level). In the first step of the reaction, acetyl-CoA is deprotonated by Asp375, producing an intermediate, which is the nucleophile for attack on the second substrate, oxaloacetate, prior to hydrolysis of the thioester bond of acetyl-CoA and release of the products. A central question has been whether the nucleophilic intermediate is the enolate of acetyl-CoA, the enol, or an 'enolic' intermediate stabilized by a 'low-barrier' hydrogen bond with His274 at the active site. The imidazole sidechain of His274 is neutral, and donates a hydrogen bond to the carbonyl oxygen of acetyl-CoA in substrate complexes. We have investigated the identity of the nucleophilic intermediate by QM/MM calculations on the substrate (keto), enolate, enol and enolic forms of acetyl-CoA at the active site of citrate synthase. The transition states for proton abstraction from acetyl-CoA by Asp375, and for transfer of the hydrogen bonded proton between His274 and acetyl-CoA have been modelled approximately. The effects of electron correlation are included by MP2/6-31G(d) and MP2/6-31+G(d) calculations on active site geometries produced by QM/MM energy minimization. The results do not support the hypothesis that a low-barrier hydrogen bond is involved in catalysis in citrate synthase, in agreement with earlier calculations. The acetyl-CoA enolate is identified as the only intermediate consistent with the experimental barrier for condensation, stabilized by conventional hydrogen bonds from His274 and a water molecule.     

A theoretical study on reaction pathways to carbanions

Different substituents (NO2, C6H5, NH2, NH-CH=CH-CHO) to a methylene group were taken into account to investigate under which conditions the mechanism of formation of carbanions by proton transfer to a base (methylamine) can be favorable, as a preliminary study of the reaction catalyzed by semicarbazide-sensitive amine oxidases. Three different approaching paths of methylamine to C(alpha) in NO2-C(alpha)H2-NO2, and the relevant potential energy surfaces, were examined at the SCF/3-21G and 6-31G* levels. The proton transfer along the first two paths occurred with a similar barrier, which became fairly consistent after including the MP2 correlation correction, with either basis set, while the last approaching path was abandoned. For the other model systems the minimum was searched only at the 3-21G level in the vicinity of the first reaction path. The substitution of a nitro group with a phenyl group sharply raised the barrier for the proton transfer to methylamine. Also by substituting the second nitro group with either -NH2 or -NH-CH=CH-CHO, a steep uphill pathway was found. A more realistic model of the substrate-cofactor complex, namely the Schiff base between benzylamine and pyridoxal, again produced a barrier, almost matching that obtained for C6H5-C(alpha)H2-NO2. In both cases, the energy profiles for the rotation about the CC(alpha)NC dihedral and the proton shift tautomers were also considered at the 3-21G and 6-31G* levels. A preliminary scan of the effect of methyl (or methylphosphate) substitutions to the pyridoxal ring was performed and the stability of the Schiff bases involving other cofactors was also considered.     

MD and QM/MM study on catalytic mechanism of a FAD-dependent enzyme ORF36: For nitro sugar biosynthesis

The catalytic mechanism of a FAD-dependent nitrososynthase (ORF36) was studied with molecular dynamics (MD) and quantum mechanical/molecular mechanical (QM/MM) methods. Residues Leu160 and Phe374 play an important role during the FAD binding with ORF36. Similar phenylalanine/leucine pair was found in the other two enzymes of this family. For the second oxidation step of ORF36 toward thymidine diphosphate-l-epi-vancosamine, three elementary catalytic steps were found: a hydroxylation step, a hydrogen back-transfer step and a hydroxyl group elimination step. The hydroxylation step is found to be the rate-determining step with an energy barrier of 26.3 kcal/mol under the B3LYP/cc-pVTZ//CHARMM22 level. Two possible pathways for the second oxidation step are carefully investigated. Our simulations indicate that an oxygen atom from the coenzyme FADHOOH is inserted into the product. In addition, the electrostatic influence of 17 individual residues and five neighboring water molecules on the rate-determining step was estimated. The results indicate that groups Gly132/Ala133/Leu134, Met375/Gln376 and a water fence play a key role in facilitating the rate-determining step. On the other hand, residues Leu160, Val161 and Ser162 are found to be critical to suppress the rate-determining step. Our results lead to further understanding of the detailed catalytic pathways for nitro sugar biosynthesis.     

A quantum chemical study of the mechanism of action of Vitamin K epoxide reductase (VKOR) II. Transition states

A reaction path including transition states is generated for the Silverman mechanism [R.B. Silverman, Chemical model studies for the mechanism of Vitamin K epoxide reductase, J. Am. Chem. Soc. 103 (1981) 5939-5941] of action for Vitamin K epoxide reductase (VKOR) using quantum mechanical methods (B3LYP/6-311G**). VKOR, an essential enzyme in mammalian systems, acts to convert Vitamin K epoxide, formed by Vitamin K carboxylase, to its (initial) quinone form for cellular reuse. This study elaborates on a prior work that focused on the thermodynamics of VKOR [D.W. Deerfield II, C.H. Davis, T. Wymore, D.W. Stafford, L.G. Pedersen, Int. J. Quant. Chem. 106 (2006) 2944-2952]. The geometries of proposed model intermediates and transition states in the mechanism are energy optimized. We find that once a key disulfide bond is broken, the reaction proceeds largely downhill. An important step in the conversion of the epoxide back to the quinone form involves initial protonation of the epoxide oxygen. We find that the source of this proton is likely a free mercapto group rather than a water molecule. The results are consistent with the current view that the widely used drug Warfarin likely acts by blocking binding of Vitamin K at the VKOR active site and thereby effectively blocking the initiating step. These results will be useful for designing more complete QM/MM studies of the enzymatic pathway once three-dimensional structural data is determined and available for VKOR.