2,6-二巯基嘌呤质子转移异构化的密度泛函理论研究

采用密度泛函B3LYP方法, 在6-311＋G(d,p)基组水平上对2,6-二巯基嘌呤质子转移引起的硫醇式与硫酮式互变异构反应进行了计算研究, 获得了互变异构过程的反应焓、活化能、活化吉布斯自由能和质子转移反应的速率常数等性质. 计算结果表明, 2,6-二巯基嘌呤无论是孤立分子还是一水合物, 其二硫酮式R是最稳定异构体. 由二硫酮式通过分子内质子转移向二硫醇式异构化共有6条反应通道, 其主通道(1)速控步骤的活化能为139.1 kJ•mol^－1, 速率常数为2.16×10^－12 s^－1; 当水分子参与反应以双质子转移机理异构化时, 活化能显著降低, 有利于硫酮式向硫醇式转变, 其主通道(7)速控步骤的活化能为61.3 kJ•mol^－1, 速率常数为1.33×10 s^－1. 计算结果还表明, 氢键作用在增大2,6-二巯基嘌呤氢键一水合物稳定性、降低质子转移异构化反应活化能等方面起着重要的作用.     

5-氯尿嘧啶质子转移异构化的密度泛函理论研究

采用密度泛函B3LYP/6-311+G**方法,对5-氯尿嘧啶分子内质子转移及水助催化质子转移引起的互变异构反应机理进行了计算研究,获得了互变异构过程的反应焓、活化能、活化吉布斯自由能和质子转移反应的速率常数等参数。计算结果表明,5-氯尿嘧啶无论是孤立分子还是一水合物,其双酮式CU1是最稳定异构体,由双酮式向烯醇式异构化找到3条通道(P1,P2,P3),各通道速控步骤的活化能分别为177.85、177.05和197.58kJ/mol。当水分子参与反应以双质子转移机理异构化时,活化能显著降低,各通道速控步骤的活化能依次降为66.24、69.36和77.85kJ/mol,有利于双酮式向烯醇式或酮醇式转变。计算结果还表明,氢键作用在增大5-氯尿嘧啶一水复合物稳定性、降低质子转移异构化反应活化能等方面起着重要作用。     

5-氟胞嘧啶气相及水助质子转移异构化的理论研究

采用密度泛函B3LYP/6-311G**方法,对6种5-氟胞嘧啶异构体孤立分子的稳定性及质子转移引起的酮式-烯醇式、氨基式-亚胺式互变异构反应机理进行了计算研究,获得了零点能、吉布斯自由能及质子转移过程的反应焓、活化能、活化吉布斯自由能和速率常数等参数.计算结果表明,气相中烯醇-氨基式FC4是最稳定的异构体.分子内质子转移设计了FC1→FC2和FC1→FC6两条通道,分别标记为P(1)和P(2),各通道速控步骤的活化能和速率常数分别为155.9 kJ·mol-1,4.70×10-15 s-1和173.1 kJ·mol-1,1.41×10-18 s-1.水助催化时,相应通道P(3) 和P(4) 速控步骤的活化能和速率常数分别为51.0 kJ·mol-1,1.41×103 s-1和88.2 kJ·mol-1,4.53×10-3 s-1.可见,水分子的加入极大地降低了质子转移的活化能垒.另外发现,水分子参与形成协同的双质子转移机理比水助单质子转移机理更利于降低活化能垒.     

二苯甲酰甲烷互变异构反应机理溶剂效应的理论研究

采用密度泛函B3LYP方法,在6-31+G(d,p)基组水平上对二苯甲酰甲烷质子转移引起的酮式-烯醇式互变异构反应机理进行了计算研究,获得了零点能、总能量、吉布斯自由能及质子转移过程的反应焓、活化能、活化吉布斯自由能和速率常数等参数.3种非质子溶剂中的优化和频率计算采用Onsager模型进行计算.计算结果表明,不论在气相还是3种溶剂中,二苯甲酰甲烷的烯醇式较酮式稳定,烯醇式向酮式气相转变需要较高的活化能垒,在不同极性的溶剂中,随着溶剂介电常数的增大,异构化反应活化能垒减小,反应速率常数增大.     

烯丙基类不对称醚异构化反应机理的理论研究

采用MP2和密度泛函M06-2X方法,在6-31++G(d,p)基组水平上对烯丙基类不对称醚异构化反应机理进行了计算研究. 揭示了其可能的反应途径,预测了互变异构吉布斯自由能,活化能等性质. 计算结果表明,在没有金催化剂的条件下,尽管有醇溶剂时异构化活化能垒有所降低,异构化反应依然不容易进行. 相反,存在金催化剂并且有醇溶剂情况下,烯丙基类不对称醚异构化反应活化自由能大大降低,仅为7.5 kcal/mol. 通过比较有无醇溶剂和金催化剂对异构化的影响,揭示了金烯烃络合和醇分子参与反应以质子转移的异构化反应机理,很好解释了实验中观察的现象. 计算结果还表明：醇分子不仅参与反应提供质子转移,它还能与醚竞争金催化剂络合,因此在高浓度醇条件下会抑制异构化反应进行.     

p-羟基苯丙酮酸互变异构反应机理溶剂效应的理论研究

采用密度泛函B3LYP方法,在6-31+G**基组水平上对p-羟基苯丙酮酸质子转移引起的烯醇式-酮式互变异构反应机理进行了计算研究,获得了互变异构过程的反应焓、活化能、活化吉布斯自由能和质子转移反应的速率常数等参数。Onsager反应场溶剂模型用于水相和甲醇相的计算,计算结果表明,p-羟基苯丙酮酸无论在气相还是水相、甲醇相中,其酮式异构体都是最稳定的存在形式,烯醇式异构体也会以一定量共存。气相中孤立分子内质子转移几何构型改变较显著,异构化反应需要较大的活化能,不利于发生质子转移,Onsager模型计算方法对质子转移反应的影响较小。     

胞嘧啶一水复合物质子转移异构化反应机理的理论研究

采用密度泛函B3LYP方法,在6-311＋G＊＊基组水平上对胞嘧啶一水复合物质子转移引起的氨-酮式、氨-烯醇式与亚胺-酮式互变异构反应机理进行了计算研究,获得了互变异构过程的反应焓、活化能、活化吉布斯自由能和质子转移反应速率常数等参数。将Onsager反应场溶剂模型用于水相的计算,结果表明,胞嘧啶一水复合物无论在气相中...     

3-羟基-3-甲基-2-丁酮和苯甲酰甲酸甲酯热分解反应的 理论研究

用从头算的方法在6-31G水平上研究了3-羟基-3-甲基-2-丁酮(1)和苯甲酰甲酸甲酯(2)热分解反应的机理。结果是:前一反应是经历五元环过渡态到达氢键中间体,它接着直接分解成乙醛的异构体和丙酮,最后乙醛的异构体异构化成乙醛;后一反应经历六元环过渡态形成中间体1(INT1),中间体1(INT1)直接分解成中间体2(INT2)和甲醛,中间体2(INT2)经过第二个过渡态分解成苯甲醛的异构体和一氧化碳,最后苯甲醛异构体异构化成苯甲醛。其中氢迁过程是反应的速控步骤。在MP~2/6-31G//HF/6-31G+ZPE水平上,对应于这两个反应速控步骤的活化位垒分别是251.42kJ/moL和247.94kJ/mol。采用传统过渡态理论计算了两反应的热反应速率常数,理论的计算结果与实验值吻合较好。     

2-(2-羟基苯亚甲基胺)-4,6-二羟基-嘧啶的质子转移异构化反应

为了探索2-(2-羟基苯亚甲基胺)-4,6-二羟基-嘧啶(M1)分子醇式和酮式结构互变异构化的反应机理,利用密度泛函理论(DFT)方法,在B3LYP/6-311+G(d,p)基组水平上,对M1化合物异构化反应的势能面进行了研究,在探讨各种可能的反应途径中,发现单体至少有8种异构体和10种过渡态.结果表明:2-(2-羟基苯亚甲基胺)-6-羟基-4(3H)嘧啶酮(M6)不论是单体、与水形成的配合物,还是二聚体,比其相对应的异构体能量低,表明在通常情况下是以M6形式稳定存在的;在考察的可能反应途径中,直接进行的分子内质子转移过程需要的活化自由能为143.8 kJ· mol-1,水助催化时,反应的活化自由能为38.9 kJ· mol-1,二聚体双质子转移的活化自由能为0.6 kJ·mol-1,二聚体双质子转移所需活化自由能最低,在室温下就可以进行,由此可见氢键在降低反应活化能方面起着重要的作用.     

三羟甲基氨基甲烷水杨醛类席夫碱质子转移异构化的理论研究

采用密度泛函(DFT)中的B3LYP方法,在6-311+G(d,p)基组水平上对三羟甲基氨基甲烷水杨醛席夫碱气相、水溶剂及甲醇溶剂中的分子内质子转移和衍生席夫碱的互变异构反应机理进行了计算研究,获得了反应焓、活化能、活化吉布斯自由能和质子转移反应的速率常数等参数.液相计算采用Onsager模型.结果表明,不论在气相、水溶剂还是甲醇溶剂中,三羟甲基氨基甲烷水杨醛席夫碱(L3=H,L5=H)烯醇亚胺式异构体R1和醌型的酮烯胺异构体P1可以共存,但以苯环型的烯醇亚胺式R1为主要形式.当由苯环型的烯醇亚胺R1向醌型的酮烯胺P1分子内质子转移时活化能较低,室温常压下反应容易进行.水和甲醇溶剂对异构化反应影响较小.当—NO2,—OMe取代生成衍生席夫碱时(L3=H,L5=NO2;L3=OMe,L5=H),结果表明,苯环型的烯醇亚胺式和醌型的酮烯胺式异构体都能共存,质子转移异构化反应的活化能垒也较低.     

喹唑啉酮互变异构的理论计算

采用密度泛函理论B3LYP/6-311+G(d,p)方法,计算并考察了喹唑啉酮进行结构互变的质子迁移过程的两种可能途径:(a)分子内质子迁移,(b)水助质子迁移.结果表明,途经b所需要的能垒小,氢键在降低反应能垒方面起重要作用.     

别嘌醇质子迁移过程的理论研究

别嘌醇(Allopurinol)是次黄嘌呤的位置异构体,是唯一在临床上应用的黄嘌呤氧化酶抑制剂.     

3-卤代吡唑质子转移过程的理论研究

采用密度泛函B3LYP/6-311G**方法,对3-卤(-F、-Cl、-Br)代吡唑几何构型进行了全自由度优化,获得了它们的几何结构和电子结构。计算结果显示,N1-H型的稳定性大于N2-H型。计算并考察了3-卤代吡唑进行结构互变的质子转移过程的四种可能途径:(a)分子内质子转移;(b)水助质子转移;(c)同种二聚体双质子转移;(d)异种二聚体双质子转移。计算结果表明(以3-氟代吡唑为例),途径d所需要的活化能最小(54.89 kJ/mol),而途径a所需要的活化能最大(198.83kJ/mol),途径b和c的活化能居中间分别为(104.05 kJ/mol和69.05 kJ/mol)。研究还表明氢键在降低活化能方面起着重要的作用,卤素(-F、-Cl、-Br)对活化能的影响不大。     

Ab initio study of 2,4-substituted azolidines. II. Amino-imino tautomerism of 2-aminothiazolidine-4-one and 4-aminothiazolidine-2-one in water solution

The relative stabilities of the tautomers of 2-aminothiazolidine-4-one and 4-aminothiazolidine-2-one were calculated at the MP2/6-31+G(d,p) level by considering their mono- and trihydrated complexes. Single-point calculations at the MP4/6-31+G(d,p)//MP2/6-31+G(d,p) level of theory were performed to obtain more accurate energies. The values of proton transfer barriers in the isolated, mono- and trihydrated tautomers of 2-aminothiazolidine-4-one (2AT) and 4-aminothiazolidine-2-one (4AT) were calculated for two different mechanisms of tautomerisation. In the absence of water, the process of proton transfer should not occur. Addition of water molecules decreases the barrier making the process faster, as the participation of two water molecules in a proton transfer reaction is more favorable than the participation of only one water molecule. To estimate the effect of the medium (water) on the relative stabilities of the tautomers of the studied compounds we applied the polarizable continuum model (PCM). (13)C NMR chemical shieldings were calculated using the GIAO approach at MP2/6-31+G(d,p) optimized geometries. HF and the DFT B3LYP functional with 6-31+G(d,p) basis set were employed. The quantum chemical results for the chemical shifts in gas phase and in polar solvents (water and DMSO) were compared with experimental data. TD DFT B3LYP/aug-cc-pVTZ calculations were performed to predict the absorption maxima of tautomers A and B of 2AT and 4AT.     

2-巯基吡啶质子迁移过程的理论研究

采用密度泛函理论,在B3LYP/6-311G**基组水平上,计算并考察了2-巯基吡啶分子硫醇式结构和硫酮式结构进行结构互变质子迁移过程中的4种可能途径:(a)分子内质子迁移,(b)水助质子迁移,(c)同种二聚体双质子迁移和(d)异种二聚体间双质子迁移.计算结果表明,途经c所需要的活化能最小(9.73 kJ.mol-1,逆反应则为55.28 kJ.mol-1),而过程a所需要的活化能最大(106.02 kJ.mol-1),途径b和d的活化能居中间(分别为32.05和15.91 kJ.mol-1).研究还表明,氢键在降低反应活化能方面起着重要的作用.     

Excited state proton transfer in guanine in the gas phase and in water solution: a theoretical study

Theoretical investigations were performed to study the phenomena of ground and electronic excited state proton transfer in the isolated and monohydrated forms of guanine. Ground and transition state geometries were optimized at both the B3LYP/6-311++G(d,p) and HF/6-311G(d,p) levels. The geometries of tautomers including those of transition states corresponding to the proton transfer from the keto to the enol form of guanine were also optimized in the lowest singlet pipi* excited state using the configuration interaction singles (CIS) method and the 6-311G(d,p) basis set. The time-dependent density function theory method augmented with the B3LYP functional (TD-B3LYP) and the 6-311++G(d,p) basis set was used to compute vertical transition energies using the B3LYP/6-311++G(d,p) geometries. The TD-B3LYP/6-311++G(d,p) calculations were also performed using the CIS/6-311G(d,p) geometries to predict the adiabatic transition energies of different tautomers and the excited state proton transfer barrier heights of guanine tautomerization. The effect of the bulk aqueous environment was considered using the polarizable continuum model (PCM). The harmonic vibrational frequency calculations were performed to ascertain the nature of potential energy surfaces. The excited state geometries including that of transition states were found to be largely nonplanar. The nonplanar fragment was mostly localized in the six-membered ring. Geometries of the hydrated transition states in the ground and lowest singlet pipi* excited states were found to be zwitterionic in which the water molecule is in the form of hydronium cation (H3O(+)) and guanine is in the anionic form, except for the N9H form in the excited state where water molecule is in the hydroxyl anionic form (OH(-)) and the guanine is in the cationic form. It was found that proton transfer is characterized by a high barrier height both in the gas phase and in the bulk water solution. The explicit inclusion of a water molecule in the proton transfer reaction path reduces the barrier height drastically. The excited state barrier height was generally found to be increased as compared to that in the ground state. On the basis of the current theoretical calculation it appears that the singlet electronic excitation of guanine may not facilitate the excited state proton transfer corresponding to the tautomerization of the keto to the enol form.     

Oxo-hydroxy tautomerism of 5-fluorouracil: water-assisted proton transfer

Post-Hartree-Fock ab initio quantum chemical calculations were performed for 5-fluorouracil in the gas phase and in a three-water cluster. Full geometry optimizations of the 5-fluorouracil-water complexes were carried out at the MP2/6-31+G(d,p) level of theory. MP4/6-31+G(d,p)//MP2/6-31+G(d,p) and MP4/6-31++G(d,p)//MP2/6-31+G(d,p) single-point calculations were performed to obtain more accurate energies. In water solution, 5-fluorouracil exists mainly in the 2,4-dioxo form (A). We propose that the populations of the 2-hydroxy-4-oxo (B) and 4-hydroxy-2-oxo (D) tautomers are 1 x 10(-4)% and 3.9 x 10(-8)%, respectively, on the basis of the relative stabilities of the tautomers calculated at the MP4/6-31++G(d,p)//MP2/6-31+G(d,p) level of theory. A profound difference between isolated and hydrated 5-fluorouracil is noted for the height of the tautomerization barrier. In the absence of water, the process of proton transfer is very slow. The addition of water molecules decreases the barrier by 2.3 times, making the process much faster. The minimum energy path (MP2/6-31+G(d,p)) for water-assisted proton transfer in trihydrated 5-fluorouracil was followed. CNDO/S-CI calculations predict singlet pi-pi(*) electron transitions at 312 nm for B and at 318 nm for D. The fluorescence spectrum of 5-fluorouracil in water confirms the presence of the hydroxy tautomer.     

2-羟基吡啶质子转移过程的理论研究

采用量子化学中的密度泛函理论,在B3LYP/6-31G(d)基组水平上,计算并考察了2-羟基吡啶分子醇式结构和酮式结构进行结构互变的质子转移过程中的4种可能途径:(a)分子内质子转移,(b)水助催化质子转移,(c)同种二聚体双质子转移和(d)异种二聚体间双质子转移.计算结果表明,途经c所需要的活化能最小(2.6 kJ•mol－1,逆反应则为27.1 kJ•mol－1),而过程a所需要的活化能最大(137.2 kJ•mol－1),途径b和d的活化能居中间(分别为38.7和17.3 kJ•mol－1).研究还表明,氢键在降低反应活化能方面起着重要的作用.