2-羟基咪唑在气相和水中的异构平衡和质子迁移的从头算计算和Monte Carlo模拟

采用从头算方法在MP2/6-31＋G^*水平上研究了2-羟基咪唑分子在孤立分子和一水合物的异构体的相对稳定性和可能的质子迁移反应, 分析了一个水分子的参与对2-羟基咪唑分子异构体的相对稳定性和质子迁移速率的影响, 采用Monte Carlo模拟方法研究了反应体系在水溶液中反应的溶剂化效应. 结果表明: 2-羟基咪唑分子的孤立分子和一水合物的最稳定异构体相同, 都为酮式. 直接质子迁移反应在水溶液中活化能垒有轻微增加, 但产物能量得到降低; 水助催化质子迁移反应在水溶液中的活化能垒和产物能量都得到明显降低. 综合气相和水相的计算结果, 2-羟基咪唑水助催化的质子迁移反应较易进行, 且在水溶液中进行容易, 可以很容易被实验观察到.     

2',3'-二脱氧-2',3'-二去氢鸟嘌呤核苷构象稳定性的理论研究

采用密度泛函方法在B3LYP/6-31＋G**水平上研究了2',3'-二脱氧-2',3'-二去氢鸟嘌呤核苷分子(D4G)的构象. 分别研究在气相中的孤立分子和一水合物异构体的相对稳定性和异构体之间的相互转变过程, 分析了水分子的参与对D4G异构体的相对稳定性和几何结构参数以及自然电荷的影响. 结果表明, 孤立的D4G分子在气相中存在8种稳定构象, 其中构象d4g-2是所有构象中最稳定的, 气相中D4G主要以d4g-2存在. 气相中各构象的相对稳定性为: d4g-2＞d4g-1＞d4g-5＞d4g-3＞d4g-6＞d4g-4＞d4g-8＞d4g-7. 计算得到的各构象键长和键角数据与实验值接近. 一个水分子的加入对D4G分子的构型参数有所影响, 基本不改变D4G分子各构象的稳定性顺序, 但构象转变的能垒有所提高. 氢键在分子构象中发挥了重要作用.     

2',3'-二脱氧-2',3'-二去氢鸟嘌呤核苷构象稳定性的理论研究

采用密度泛函方法在B3LYP/6-31＋G**水平上研究了2',3'-二脱氧-2',3'-二去氢鸟嘌呤核苷分子(D4G)的构象. 分别研究在气相中的孤立分子和一水合物异构体的相对稳定性和异构体之间的相互转变过程, 分析了水分子的参与对D4G异构体的相对稳定性和几何结构参数以及自然电荷的影响. 结果表明, 孤立的D4G分子在气相中存在8种稳定构象, 其中构象d4g-2是所有构象中最稳定的, 气相中D4G主要以d4g-2存在. 气相中各构象的相对稳定性为: d4g-2＞d4g-1＞d4g-5＞d4g-3＞d4g-6＞d4g-4＞d4g-8＞d4g-7. 计算得到的各构象键长和键角数据与实验值接近. 一个水分子的加入对D4G分子的构型参数有所影响, 基本不改变D4G分子各构象的稳定性顺序, 但构象转变的能垒有所提高. 氢键在分子构象中发挥了重要作用.     

水合甘氨酸体系水桥式质子迁移的理论研究

水分子与甘氨酸作用会导致甘氨酸羧基上的质子迁移到氨基上,质子可以通过水分子链进行迁移。采用密度泛函理论的B3LYP/6-31++g**方法研究了水分子链的逐渐增长（1~5个水分子）对质子迁移的影响,发现水分子数少于5时,质子迁移一步完成;水分子数为5时,质子迁移经由一个中间体,需两步完成;水分子链的增长使质子迁移反应的自由能越来越低,但是反应的活化能越来越高,即在热力学上有利于质子迁移反应,在动力学上不利于质子迁移反应。     

4(3H)-嘧啶酮及其类似物互变异构的理论研究

采用密度泛函理论,在B3LYP/6-311G**基组水平上,计算并考察了4(3H)-嘧啶酮及其类似物(5-氟-4(3H)-嘧啶酮、4-巯基嘧啶和5-氟-4-巯基嘧啶)醇式结构和酮式结构进行结构互变质子迁移过程中的2种可能途径:(a)分子内质子迁移;(b)水助质子迁移.计算结果表明,途经b所需要的活化能较小.研究还表明,氢键在降低反应活化能方面起着重要的作用.     

2,6-硫代黄嘌呤质子互变异构反应机理的理论研究

在密度泛函B3LYP/6-311G~(**)理论水平上,对气相和水相中2,6-硫代黄嘌呤各烯醇式与酮式水助质子互变异构体及其过渡态进行几何构型全自由度优化,获得它们在气相和水相中的几何结构和电子结构,PCM反应场溶剂模型用于水相计算.结果显示在气相和水相中,水参与反应降低了互变异构质子迁移的反应活化能,对互变异构质子迁移的反应起到催化作用,但是没有改变各异构体的稳定性顺序,其顺序为W1＞W3＞W2.进一步研究了2,6-硫代黄嘌呤各烯醇式与酮式水助质子互变异构的反应机理,提出了2,6-硫代黄嘌呤各烯醇式与酮式互变异构质子迁移的反应为平面六元环的过渡态结构.探讨了溶剂化效应对互变异构体的几何结构、能量、电荷分布以及互变异构反应活化能的影响等.     

基态丙酮酸单分子反应的机理

运用MP2和B3LYP方法,对基态丙酮酸的8个可能的单分子反应的反应机理进行了详细的理论计算.研究发现,氢键相互作用及两个羰基的相对空间取向在丙酮酸异构体的稳定性中起着重要的作用.并且由最稳定的异构体Tc旋转为其它异构体的能垒都比较高,旋转异构反应在室温下是难于进行的.经由质子迁移过程可实现Ct和Tt异构体之间的转化,但反应能垒比内旋转过程高得多,说明质子迁移无法与内旋转过程相竞争.对基态丙酮酸脱羧和脱羰反应途径的研究均表明,脱羧反应是经由五中心过渡态首先生成甲基羟基卡宾中间体,然后中间体再进一步异构得到产物.四中心协同反应机理是不利的.     

微波辐射下4-硫胸腺嘧啶核苷的合成研究

提出了一种高效合成3’,5’-O-二乙酰基-4-硫胸腺嘧啶核苷的方法,该方法具有反应时间短、操作简便等优点,目标产物产率可达93.9%以上。以二氧六环做溶剂,反应温度102℃,乙酰化胸腺嘧啶核苷与P2S5的摩尔比为1∶1.2,微波功率700W辐射8分钟,即可高收率得到3’,5’-O-二乙酰基-4-硫胸腺嘧啶核苷。     

嘧啶的研究:氯和硫醇嘧啶类的作用:4-甲基-5-正丙基-6-氨基-2-氧嘧啶的合成

1. 2-乙硫基-4-甲基-5-正丙基-6-氯代嘧啶在甲醇溶液与甲醇钠相互作用,和在乙醇溶液与乙醇钠互相作用,分别得到相应的嘧啶-甲醚及嘧啶-乙醚。2. 氯对2-乙硫基-4-甲基-5-正丙基-6-氯代嘧啶,2-乙硫基-4-甲基-5-正丙基-6-甲氧基嘧啶及2-乙硫基-4-甲基-5-正丙基-6-乙氧基嘧啶起作用,得到相应的嘧啶砜,用氯在水溶液中氧化硫醇嘧啶化合物(Ⅰ),此反应是特殊的,对嘧啶环中的双键没有影响,祇有硫醇基团氧化而得到相应的稳定的嘧啶砜(Ⅱ)。3. 将醇氨与2-乙磺醯基-4-甲基-5-正丙基-6-氯代嘧啶作用时,对于2-位上的磺醯基不起反应,取代于6-位的氯被置换为氨基而得到2-乙磺醯基-4-甲基-5-正丙基-6-氨基嘧啶。4. 硷对2-乙磺醯基-4-甲基-5-正丙基-6-乙氧基嘧啶及2-乙磺醯基-4-甲基-5-正丙基-6-甲氧基嘧啶作用,分别得到2-氧代-4-甲基-5-正丙基-6-乙氧基嘧啶及2-氧代-4-甲基-5一正丙基-6-甲氧基嘧啶,因此,嘧啶砜中在2-位的磺醯基与硷反应时被置换为羟基,而取代于6-位的乙氧基和甲氧基则不起反应。5. 叙述了一种合成2-氧代-4-甲基-5-正丙基-6-氨基嘧啶(或1,2-二氢化-2-酮-4-氨基-5-正丙基-6-甲基嘧啶)的新方法。     

Cu2+对腺嘧啶-胸腺嘧啶碱基对阴离子质子转移的影响

采用B3LYP/DZP++方法研究了腺嘌呤-胸腺嘧啶(A-T)碱基对阴离子(AT)-的单质子转移机理以及金属离子Cu2+对(AT)-碱基对质子转移的影响.(AT)-碱基对的单质子转移路径是由胸腺嘧啶N25位上的质子H26沿分子间的氢键N25-H26…N10转移到腺嘌呤的N10位.金属Cu2+可通过络合作用分别吸附在(AT)-碱基对O24、O28、N4、N13上,从而影响(AT)-碱基对中质子转移过程.Cu2+络合作用在胸腺嘧啶(T)的O24、O28上时,发生了从胸腺嘧啶到腺嘌呤方向上的单质子转移反应;而作用在腺嘌呤(A)的N4、N13上时,得到了双质子转移的稳定产物.     

DFT and MP2 study of isomery scheme in Formazan and intermolecular and intramolecular proton transfer between its tautomers

DFT and Moller Plesset (MP2) calculations were applied to study of isomery scheme in formazan. Formazan that can be presented by three tautomers and eight isomers has various applications in dyes, complexes, and biochemistry. The structures of its isomers and related transition states were optimized, and important molecular parameters, IR frequencies, and energetic results were extracted. The relative stabilities of formazan isomers in the gas phase were found to be as 1EZ >1ZZ >1EE >1ZE >2EE > 3 >2EZ >2ZZ >TS1 >TS3 >TS2. Thermodynamic data confirms that tautomer 1 is major tautomer, and all possible tautomerism interconversions have small rates at room temperature. Then, relative stabilities were calculated in different solvents (chloroform, tetrahydrofuran, acetone, and water). The relative stabilities and thermodynamic data in solvents are nearly similar to those in the gas phase, but the rate constants are slightly more than that in the gas phase. Moreover, relative stabilities of formazan isomers and intermolecular proton transfer in presence of one to three molecules of water have been studied. The results showed that activation barriers in water‐assisted tautomerisms are in general lower than those in the gas phase, but the relative energies of isomers do not change importantly by water clusters. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

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.     

Kinetic and thermodynamic study of inter‐ and intramolecular proton transfer in N′‐acetyl formohydrazide tautomers

Nine tautomers and eleven possible tautomeric interconversions of N′‐acetyl formohydrazide have been studied at B3LYP/6‐311++G** level of theory. From these calculations, optimized geometries, molecular parameters, IR frequencies, NMR chemical shifts, and energetic results are obtained. In all tautomers except tautomers 4, E isomer is more stable than Z isomer. Energetic data were used to calculate the energy barriers of tautomeric interconversions and very high energy barriers were obtained for all tautomeric interconversions. Moreover, study of solvent effects on relative stabilities of tautomers and transition states showed that they are similar to those in the gas phase. In addition, intermolecular proton transfer with the assistance of one to three water molecules has been studied and the results showed that activation barriers in water‐assisted tautomerism are in general lower than those in the gas phase. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

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.     

DFT and MP2 Calculations on Tautomers and Water-Assisted Proton Transfer on 1,2,5-Oxadiazol-4,3-diamine

In this work a detailed quantum-chemical comparison of the relative stability of six tautomers of 1,2,5-oxadiazol-4,3-diamine studied in the gas phase and solution. Theoretical calculations are carried out by the density functional theory (DFT/B3LYP) and MP2 methods using the standard 311++G(d,p) basis set. The results indicate that A is the most stable form in the gas phase and also is the predominant tautomer in solution at the DFT and MP2 methods. The transition states of proton transfer reaction are calculated. The variation of dipole moments and charges on atoms are studied in various solvent. Specific solvent effects with addition of one water molecule near the electrophilic centers of tautomer investigated. Also the transition state of proton transfer assisted by a water molecule was investigated.     

Systematic study of the tautomerism of uracil induced by proton transfer. Exploration of water stabilization and mutagenicity

To systematically investigate all the possible tautomerisms from uracil (U) and its enol form (U) induced by proton transfer, we describe a study of structural tautomer interconversion in the gas phase, in a continuum solvent, and in a microhydrated environment with 1 or 2 explicit water molecules, using density functional theory (DFT) calculations by means of the B3LYP exchange and correlation functions. A total of 62 geometries including 25 transition states were optimized, and the geometrical parameters have been discussed. Some rules of the configuration variation in tautomerization were summarized. The relative stabilities of all the tautomers were established. When a proton transfers from the di-keto form to the keto-enol form, water molecules in different regions show absolutely opposite effects: some assist, whereas others hinder the tautomerization. However, when a proton transfers from the keto-enol form to the di-enol form, water molecules in different regions show similar effects: the Gibbs free energy always increases and the activation energy always decreases. Additionally, some important factors that obviously affect the activation energy and Gibbs free energy were found and discussed in detail. The reasons that water molecules can assist or prevent the proton transfer were given. Furthermore, on the basis of our calculated results, we explain why it is hard to detect the di-enol form of uracil in general experiments.     

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.可见,水分子的加入极大地降低了质子转移的活化能垒.另外发现,水分子参与形成协同的双质子转移机理比水助单质子转移机理更利于降低活化能垒.     

Molecular view of the anomalous acidities of Sn2+, Pb2+, and Hg2+

Experimental results taken from both the condensed and gaseous phase show that, when associated with water, the three dications Sn(2+), Pb(2+), and Hg(2+) exhibit a facile proton-transfer reaction. In the gas phase, no stable [M.(H(2)O)(n)](2+) ions are observed; but instead the cations appear to undergo rapid hydrolysis to give ions of the form M(+)OH(H(2)O)(n-1). A series of ab initio calculations have been undertaken on the structures and proton-transfer reaction profiles associated with the complexes [M.(H(2)O)(2,4)](2+), where M is one of Sn, Pb, Hg, and Ca. The latter has been used as a reference point both in terms of comparisons with previous calculations, and the fact that Ca(2+) is a very weak acid. The calculations show that for Sn(2+), Pb(2+), and Hg(2+), the only barriers to proton transfer are those associated with the movement of water molecules. In the gas phase, these barriers could be overcome through energy gained during ion formation, and in the condensed phase the thermal motion of water molecules would be sufficient. In contrast, the calculations show that for Ca(2+) it is the proton-transfer step that provides the most significant reaction barrier. Proton transfer in Sn(2+) and Pb(2+) is further assisted by distortions in the geometries of [M.(H(2)O)(2,4)](2+) complexes due to voids created by the 5s(2) (6s(2)) inert lone pair. For Hg(2+), ease of proton transfer is derived partly from the high degree of covalent bonding found in both the reactants and products.     

Relative binding affinities of molecular capsules investigated by ESI-mass spectrometry

ESI-MS(/MS) has been used as a method which allows the fast, unambiguous and sensitive simultaneous detection and relative stability approximation of supramolecular assemblies in mixtures. In spite of the obvious fundamental differences between solution and gas phase, ESI-MS in the case of self-assembled molecular capsules has been shown to produce very similar results to single binding experiments monitored by NMR titrations as well as conformational searches performed by Monte-Carlo simulations. MS/MS experiments reveal the same relative order of gas phase stabilities as previously found in solution. Moreover, proton transfer reactions which lead to new molecular capsules, are not detectable in the time-averaged NMR spectrum. However, the newly produced species are found in the complex mixtures by ESI-MS and can be conveniently characterized by subsequent MS/MS experiments: in a collision-induced dissociation the single half-spheres are easily discovered and structurally assigned. Thus, ESI-MS has worked as a valuable tool for the rapid screening of complex supramolecular mixtures and in combination with MS/MS experiments elucidated both the path of unexpected side reactions as well as the thermodynamic gas-phase stabilities of all components in the mixture.