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)对活化能的影响不大。     

5-醛基-1,2,3-三唑互变异构的理论研究

用密度泛函B3LYP/6-311++G**方法,对气相和水相中的5-醛基-1,2,3-三唑互变异构体进行了几何构型全自由度优化,获得了它们在气相和水相中的几何结构和电子结构.讨论了不同的取代基团和溶剂化效应对互变异构体的几何结构、能量以及互变异构反应活化能的影响.并进一步研究了N1-H、N2-H和N3-H型三唑之间的互变异构机理:(a)分子内质子转移;(b)水助质子转移.计算结果表明,途径(b)所需要的活化能较小,平均为101.31kJ/mol,途径(a)为211.70kJ/mol.     

5-卤代1,2,3,-三唑互变异构的密度泛函理论研究

用密度泛函B3LYP/6-311 G**方法,对气相和水相中的1,2,3三,-唑及5-卤(-F、C l和-B r)代1,2,3,-三唑互变异构体进行了几何构型全自由度优化,获得了它们在气相和水相中的几何结构和电子结构。计算结果显示,在气相和水相中1,2,3三-唑和5-卤代1,2,3三-唑的N2-H型要比对应的N1-H型和N3-H型稳定。讨论了不同的取代基团和溶剂化效应对互变异构体的几何结构、能量和电荷分布以及互变异构反应活化能的影响带。并进一步研究了N1-H、N2-H和N3-H型三唑之间的互变异构机理:(a)分子内质子转移;(b)水助质子转移。计算结果表明,途径(b)所需要的活化能较小,为120.06KJ/mol,途径(a)为204.12KJ/mol。     

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).研究还表明,氢键在降低反应活化能方面起着重要的作用.     

6-甲基-4-羟基嘧啶单体及二聚体质子转移过程的理论研究

应用密度泛函理论的B3LYP/6-311+G(d)方法研究了6-甲基-4-羟基嘧啶单体及二聚体质子转移的异构化反应.对反应势能面的研究发现,该化含物可能存在9种单体异构体,对其最稳定的单体构型进行分析.各单体间异构化反应的过渡态共有9种,反应的活化能最小为22.06 kJ/mol,最大为356.55 kJ/mol,最可能的反应路径在室温下即可进行. 研究了2种二聚体及其异构化反应的过渡态,发现二聚体均比其对应的单体稳定,而且质子转移所需要的活化能仅为20.13 kJ/mol,比单体低很多. 氢键在这种变化中起了主要作用,由单体和二聚体的总能量计算了氢键的键能.     

甲基重氮烷阳离子（CH3N2^＋）与DNA碱基的甲基化反应

在B3LYP／6．31G^＊＊和MP2／6．311G^＊＊水平上,计算了终致癌物甲基重氮烷阳离子（CH3N2^＋）在气相和水溶剂中与四种DNA碱基上10个亲核位发生甲基化过程的反应机理．结果表明：在气相和溶剂中,终致癌物CH3N2^＋对碱基上相应亲核位的甲基化反应过程,活化能都比较小（〈33．5kJ／mol）,是容易进行的放热反应．所以,二甲基亚硝胺一旦形成终致癌物,则很容易与碱基形成癌变物．     

尿酸质子转移反应的理论研究I——分子内质子转移研究

使用量子化学中的Hartree-Fock方法和密度泛函理论中的B3LYP方法，分别在3－21G^*和6－31G（d)水平上，计算了尿酸分子从三羰基异构体向三羟基异构体的转化。结果表明，转化过程经历了单羟基和双羟基异构体2种中间物和3种过渡态时的分子内质子转移（IPT），转移中的H原邻近的N，O和C原子形成了具有四元环结构的过渡态。随着IPT的进行，N－H键逐渐被削弱和断裂，O－H键则逐渐生成。3个反应的活化能分别为190.3kJ/mol,181.4kJ/mol和249.9kJ/mol(B3LYP/6-31G(d))。较高的活化能表明在室温下，无催化剂的IPT难以进行。     

3,6-二羟基哒嗪质子转移过程的理论研究

为了探索3,6-二羟基哒嗪分子醇式和酮式结构互变异构化的反应机理,本研究组采用DFT B3LYP/6-311+G(d)方法对标题化合物异构化反应的势能面进行了研究,在探讨各种可能的反应途径中,发现至少有34种异构体和43种过渡态.结果表明,6-羟基-3(2H)-哒嗪酮不论是单体,与水形成配合物,还是二聚体,比其相对应的异构体能量低,表明在通常情况下是以6-羟基-3(2H)-哒嗪酮形式稳定存在的,这与前人通过实验数据对互变异构体的比率进行预测的结果是一致的;在考察的可能反应途径中,直接进行的分子内质子转移过程需要的活化能为142.2 kJ·mol^-1,水助催化时,反应活化能为41.3 kJ·mol^-1,考虑溶剂效应后,其活化能为59.2 kJ·mol^-1,二聚体双质子转移的活化能为16.8 kJ·mol^-1,二聚体双质子转移所需活化能最低,在室温下就可以进行.由此可见氢键在降低反应活化能方面起着重要的作用.     

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).研究还表明,氢键在降低反应活化能方面起着重要的作用.     

苯硼酸衍生物的晶体工程：用能量较高的syn,syn-构象设计共晶

研苯硼酸衍生物有3种不同构象（syn,syn-;syn,anti-和anti,anti-构象）,这3种构象能量非常不同.苯硼酸衍生物通常是采用能量最低的syn,anti-构象存在.然而在配合物中,苯硼酸衍生物展现了硼酸功能团的构象多样性.在本文中,我们采用syn,syn-构象的硼酸衍生物作为结构单元,与含氮给体化合物反应,合成了一系列共晶配合物.采用的含氮给体化合物包括含1,2-重氮片段的化合物（阿普唑仑,1H-四唑,乙酰唑胺和苯并三唑）、邻菲罗啉和2,2′-联吡啶.结果表明,只有1,2-重氮片段不一定能与syn,syn-构象的硼酸及衍生物生成共晶.     

Effect of the H-bonding on aromaticity of purine tautomers

Four tautomers of purine (1-H, 3-H, 7-H, and 9-H) and their equilibrium H-bonded complexes with F(-) and HF for acidic and basic centers, respectively, were optimized by means of the B3LYP/6-311++G(d,p) level of theory. Purine tautomer stability increases in the following series: 1-H < 3-H < 7-H < 9-H, consistent with increasing aromaticity. Furthermore, the presence of a hydrogen bond with HF does not change this order. For neutral H-bonded complexes, the strongest and the weakest intermolecular interactions occur (-14.12 and -10.49 kcal/mol) for less stable purine tautomers when the proton acceptor is located in the five- and six-membered rings, respectively. For 9-H and 7-H tautomers the order is reversed. The H-bond energy for the imidazole complex with HF amounts to -14.03 kcal/mol; hence, in the latter case, the fusion of imidazole to pyrimidine decreases its basicity. The ionic H-bonds of N(-)···HF type are stronger by ~10 kcal/mol than the neutral N···HF intermolecular interactions. The hydrogen bond N(-)···HF energies in pyrrole and imidazole are -32.28 and -30.03 kcal/mol, respectively, and are substantially stronger than those observed in purine complexes. The aromaticity of each individual ring and of the whole molecule for all tautomers in ionic complexes is very similar to that observed for the anion of purine. This is not the case for neutral complexes and purine as a reference. The N···HF bonds perturb much more the π-electron structure of five-membered rings than that of the six-membered ones. The H-bonding complexes for 7-H and 9-H tautomers are characterized by higher aromaticity and a much lower range of HOMA variability.     

The low energy tautomers and conformers of the dipeptides HisGly and GlyHis and of their sodium ion complexes in the gas phase

The low-lying conformers of the dipeptides HisGly and GlyHis, and of their sodium cation complexes, have been studied with a combination of Monte Carlo search with the Amber force field and local geometry optimization at the ab initio HF/6-31G(d) level, completed with MP2(full)/6-311+G(2d,2p) energetics at the HF/6-31G(d) geometries. For each dipeptide, both the N(delta)-H and N(epsilon)-H tautomers of the imidazole side chain of His were considered. For each of the four isomeric dipeptides, 20-30 conformers were fully characterized at the ab initio level. All low energy structures are found to involve H-bonding at the N(delta) position of imidazole, either as a N-H donor or a N acceptor, depending upon the tautomer. In three out of the four species, the most stable conformer involves a C-terminus carboxylic acid in its less favorable trans conformation, in order to maximize intramolecular H bonding. It turns out that the lowest energy tautomer of HisGly is N(epsilon)-H, while that of GlyHis is N(delta)-H. This result argues in favor of the diversity of His tautomeric states in peptides and proteins. The sodium cation complexes of both GlyHis and HisGly have been studied as well, again considering both tautomers in each case. In three out of the four species, the most stable structure involves chelation of sodium by the two carbonyl oxygens and the imidazole ring. On the contrary, the sodium complex of the N(delta)-H tautomer of HisGly favors chelation to the peptidic carbonyl oxygen, the imidazole ring and the amino terminus. In the N(epsilon)-H tautomers of both peptides, the most favorable binding site of imidazole is the N(delta) nitrogen, while in the N(delta)-H tautomers, it is the pi cloud which provides side chain interaction. As a result, both GlyHisNa+ and HisGlyNa+ favor the N(epsilon)-H tautomer of His, in contrast to what was found for the free peptides.     

黄嘌呤异构体质子转移异构化反应机理的理论研究

采用量子化学密度泛函B3LYP/6-31G(d,p)和M06-2X/6-311++G(d,p)方法对黄嘌呤两种酮式异构体X(1,3,7)与X(1,3,9)间质子转移引起的互变异构反应机理进行了计算研究,获得了异构化反应过程的反应焓﹑活化吉布斯自由能和质子转移反应的速率常数等参数。水相计算采用极化连续(PCM)模型。结果表明,由于可能的氢迁移顺序差异,分子内由X(1,3,7)向X(1,3,9)异构化可能共有16条反应通道,涉及11个中间体和20个过渡态,其主反应通道速控步骤的活化吉布斯自由能为183.10k J/mol,速率常数为5.17×10-20s-1,其余各通道速控步骤活化吉布斯自由能均较高,而且整体水溶剂效应不利于质子转移的发生。     

吲唑、3-卤代吲唑互变异构反应机理的理论研究

The molecular structures of indazole and 3-halogeno-indazole tautomers were calculated by the B3LYP method at the 6-311G** level, both in the gaseous and aqueous phases, with full geometry optimization.The geometry and electronic structure of the tautomers of indazole, 3-halogeno-indazole and their transition states were obtained. The Onsager solvate theory model was employed for the aqueous solution calculations.The results of the calculation indicated that the N1-H form of the studied molecule is more stable than that of the N2-H form. The influences of the different 3-halogeno and solvent effects on the geometry, energy,charge and activation energy were discussed. The reaction mechanism of the tautomerization of indazole and 3-halogeno-indazole was also studied and a three-membered cyclic transition state of the tautomer reaction has been obtained.     

A microwave and quantum chemical study of the conformational properties and intramolecular hydrogen bonding of 1-fluorocyclopropanecarboxylic acid

The structural and conformational properties of 1-fluorocyclopropanecarboxylic acid have been explored by microwave spectroscopy and a series of ab initio (MP2/6-311++G(d,p) level), density functional theory (B3LYP/aug-cc-pVTZ level), and G3 quantum chemical calculations. Four "stable" conformers, denoted conformers I-IV, were found in the quantum chemical calculations, three of which (conformers I -III) were predicted to be low-energy forms. Conformer I was in all the quantum chemical calculations predicted to have the lowest energy, conformer III to have the second lowest energy, and conformer II to have the third lowest energy. Conformers II and III were calculated to have relatively large dipole moments, while conformer I was predicted to have a small dipole moment. The microwave spectrum was investigated in the 18-62 GHz spectral range. The microwave spectra of conformers II and III were assigned. Conformer I was not assigned presumably because its dipole moment is comparatively small. Conformer II is stabilized by an intramolecular hydrogen bond formed between the fluorine atom and the hydrogen atom of the carboxylic acid group. Conformer III has a synperiplanar orientation for the F-C-C=O and H-O-C=O chains of atoms. Its dipole moment is: mua = 3.4(10), mub = 10.1(13), and muc = 0.0 (assumed) and mu(tot) = 10.6(14) x 10(-30) C m [3.2(4) D]. Several vibrationally excited states of the lowest torsional mode of each of II and III were also assigned. The hydrogen-bonded conformer II was found to be 2.7(2) kJ/mol less stable than III by relative intensity measurements. Absolute intensity measurements were used to show that the unassigned conformer I is the most abundant form present at a concentration of roughly 65% at room temperature. Conformer I was estimated to be ca. 5.0 kJ/mol more stable than the hydrogen-bonded rotamer (conformer II) and ca. 2.3 kJ/mol more stable than conformer III. The best agreement with the theoretical calculations is found in the MP2 calculations, which predict conformer I to be 5.1 kJ/mol more stable than III and 1.7 kJ/mol more stable than II.