6-硫代黄嘌呤互变异构体的密度泛函理论计算

在密度泛函B3LYP/6-311G**水平下,对14种气相和水相中可能存在的6-硫代黄嘌呤异构体进行了几何构型的全自由度优化,并计算出它们的总能量、焓、熵、吉布斯自由能。Onsager反应场溶剂模型用于水相的计算.计算结果表明,6-硫代黄嘌呤在气相中和水相中主要以硫酮的形式存在.在气相和水相中,硫酮-N7(H)均比硫酮-N9(H)更稳定.计算结果同已有实验结果一致.6-硫代黄嘌呤异构化的熵效应小,对互变异构平衡几乎没有显著的影响,而焓变对互变异构产生了主要的影响.较详细地讨论了水溶剂化作用对异构体的能量、几何结构、电荷分布和偶极矩的影响.     

2-硫代黄嘌呤互变异构体的密度泛函理论计算

在密度泛函B3LYP/6-311G**水平上,对14种气相和水相中可能存在的2-硫代黄嘌呤互变异构体进行了几何构型全自由度优化,并计算出它们的总能量、焓、熵、吉布斯自由能.Onsager反应场溶剂模型用于水相的计算.计算结果表明,2-硫代黄嘌呤在气相和水相中主要以酮式结构形式存在,与已有实验结果一致.在气相和水相中,酮式结构—N(7)(H)均比酮式结构—N(9)(H)更稳定.2-硫代黄嘌呤互变异构的熵效应小,对互变异构平衡没有显著的影响,而焓变对互变异构却产生了主要的影响.水溶剂化自由能与异构体的气相偶极矩存在相关性.另外,较详细地考察了2-硫代黄嘌呤与6-硫代黄嘌呤的相对稳定性.     

6-巯基嘌呤互变异构体的密度泛函理论计算

在密度泛函B3LYP/6-311G水平下,对8种气相和水相中可能存在的6-巯基嘌呤异构体进行了几何构型的全自由度优化,并计算出它们的总能量、焓、熵、吉布斯自由能．Onsager反应场溶剂模型用于水相的计算．计算结果表明,6-巯基嘌呤在气相和水相中主要以硫酮形式存在．在气相中,硫酮．N(7)(H)要比硫酮-N(9)(H)更稳定,而在水相中,则硫酮-N(9)(H)要比硫酮-N(7)(H)更稳定．计算结果同已有实验结果一致．6-巯基嘌呤的熵效应小,对互变异构平衡几乎没有显著的影响,而焓变对互变异构产生了主要的影响．较详细地讨论了水溶剂化作用对异构体的能量、几何结构、电荷分布和偶极矩的影响,溶剂化吉布斯自由能与异构体的气相偶极矩存在相关性．     

7-硫代鸟嘌呤互变异构体的密度泛函理论计算

在密度泛函B3LYP/6-311G^**水平下,对14种气相和水相中可能存在的7-硫代鸟嘌呤异构体进行了几何构型全自由度优化,并计算出它们的总能量、焓、熵、吉布斯自由能.Onsager反应场溶剂模型用于水相的计算.计算结果表明,7-硫代鸟嘌呤在气相中和水相中主要以硫酮形式存在.在气相中,硫酮-N(6)(H)要比硫酮-N(4)(H)更稳定,而在水相中,则硫酮-N(4)(H)要比硫酮-N(6)(H)更稳定.计算结果同已有实验结果一致.7-硫代鸟嘌呤互变异构的熵效应小,对互变异构平衡没有显著的影响,而焓变对互变异构产生了主要的影响.较详细地讨论了水溶剂化作用对异构体的能量、几何结构、电荷分布和偶极矩的影响.溶剂化自由能与异构体的偶极矩存在相关性.     

5-氟胞嘧啶互变异构的密度泛函理论计算

采用BH-HLYP/6-311＋G**方法对10种气相和水相中可能存在的5-氟胞嘧啶互变异构体进行了几何全优化, 并计算出它们的总能量和吉布斯自由能. Onsager反应场溶剂模型用于水相的计算. 计算结果表明, 5-氟胞嘧啶在气相中主要以烯醇式-氨基式形式存在, 在水相中主要以酮式-氨基式形式存在. 溶剂化自由能与异构体的气相偶极矩存在相关性．进一步求得互变异构化以及构象异构化和顺反异构化的过渡态, 探讨异构化过程中的几何结构和能量的变化.     

黄嘌呤及其互变异构体的密度泛函理论研究

嘌呤碱及其衍生物在生物系统中起重要作用。对人具有兴奋和利尿作用的茶碱和咖啡碱就是黄嘌呤的甲基衍生物。黄嘌呤存在多种互变异构体，从理论计算的角度研究这些互变异构体的几何结构、电子结构及相对稳定性，进一步研究溶剂对其结构和性质的影响是有意义旧。密度泛函理论方法既考虑了电子相关，又较其它CI（组态相互作用）或MPn（n级微扰）方法节省机时。其计算结果较好，被广泛应用于研究各种化合物。本文采用密度泛函B3LYP／6—311G方法对14种黄嘌呤可能的互变异构体（见图1），分别在气相和水相中进行几何构型全自由度优化和能量计算，讨论了异构体的相对稳定性，水的溶剂化作用对异构体的能量、几何构型、电荷分布和偶极矩的影响，探讨了溶剂极性对异构体的能量和偶极矩的影响。     

2,6-二巯基嘌呤互变异构体热力学稳定性的密度泛函理论研究

利用密度泛函(DFT)B3LYP/6-311G(d,p)方法,水相计算采用自洽反应场(SCRF)中的Onsager模型,对气相和水相中可能存在的13种2,6-二巯基嘌呤互变异构体进行了全优化,并计算了各异构体的热力学参数、偶极矩及原子净电荷。计算结果表明,不论是气相还是水相,二硫酮DTP(1,3,7)是最稳定的异构体。溶剂化效应使各异构体的稳定性均增强,偶极矩大者其稳定性显著增大。溶剂化吉布斯自由能与异构体在两相中偶极矩之差存在相关性。二硫酮DTP(1,3,7)在水相中与致癌物BPDE进行亲核取代反应时,二硫酮DTP(1,3,7)中的S10原子优先进攻亲电试剂BPDE.     

2,6-二巯基吡啶互变异构平衡体系溶剂效应的理论研究

在气相及甲苯、氯仿、乙腈和水等溶剂中对2,6-二巯基吡啶及其硫酮式互变异构体进行了HF/6-31G^**水平上的优化,其中溶液中的计算采用Onsager自洽反应场(SCRF)模型.探讨了溶剂对体系几何结构和能量的影响.结果表明:溶剂的存在与极性的增加有利于平衡体系中硫酮式异构体的存在.     

5-氟尿嘧啶和5-氯尿嘧啶及其互变异构体的理论计算研究

采用HF/3-21G方法, 对6种气相和水相中可能存在的5-氟尿嘧啶(和5-氯尿嘧啶)互变异构体进行了构象分析.采用B3LYP/6-311＋G**方法对处于优势构象时的各互变异构体进行了几何全优化, 并计算出它们的总能量、焓、熵、吉布斯自由能. Onsager反应场溶剂模型用于水相的计算. 计算结果表明, 5-氟尿嘧啶和5-氯尿嘧啶在气相中和水相中主要以双酮形式存在. 5-氟尿嘧啶和5-氯尿嘧啶的熵效应小, 对互变异构平衡没有显著的影响, 而焓变对互变异构产生了主要的影响. 讨论了水溶剂化作用对异构体的能量、电荷分布和偶极矩的影响. 溶剂化自由能与异构体的气相偶极矩存在相关性. 另外, 详细地将5-氟尿嘧啶和5-氯尿嘧啶与尿嘧啶进行了对比, 获得三者最稳定异构体间电子结构异同的有用信息.     

Theoretical investigation of solvent effects on tautomeric equilibrium of 2‐diazo‐4,6‐dinitrophenol

The density functional theory has been used to study the tautomeric equilibrium of 2‐diazo‐4,6‐dinitrophenol(DDNP) in the gas phase and in 14 solvents at the B3LYP/6‐31G* level. The solvent effects on the tautomeric equilibria were investigated by the self‐consistent reaction field theory (SCRF) based on conductor polarized continuum model (CPCM) in apolar and polar solvents and by the hybrid continuum‐discrete model in protic solvent, respectively. Solvent effects on the computed molecular properties, such as molecular geometries, dipole moments, E_LUMO, E_HOMO, total energies for DDNP tautomers and transition state, tautomerization energies and solvation energies have been found to be evident. The tautomeric equilibrium of DDNP is solvent‐dependent to a certain extent. The tautomer I (cyclic azoxy form) is preferred in the gas phase, while in nonpolar solvents tautomer I and II (quinold form) exist in comparable amounts, and in highly polar solvents, the tautomeric equilibrium is shifted in favor of the more polar tautomer II . © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Reaction mechanism and tautomeric equilibrium of 2-mercaptopyrimidine in the gas phase and in aqueous solution: a combined Monte Carlo and quantum mechanics study

A combined Monte Carlo and quantum mechanical study was carried out to analyze the tautomeric equilibrium of 2-mercaptopyrimidine in the gas phase and in aqueous solution. Second- and fourth-order M?ller-Plesset perturbation theory calculations indicate that in the gas phase thiol (Pym-SH) is more stable than the thione (Pym-NH) by ca. 8 kcal/mol. In aqueous solution, thermodynamic perturbation theory implemented on a Monte Carlo NpT simulation indicates that both the differential enthalpy and Gibbs free energy favor the thione form. The calculated differential enthalpy is DeltaH(SH)(-->)(NH)(solv) = -1.7 kcal/mol and the differential Gibbs free energy is DeltaG(SH)(-->)(NH)(solv) = -1.9 kcal/mol. Analysis is made of the contribution of the solute-solvent hydrogen bonds and it is noted that the SH group in the thiol and NH group in the thione tautomers act exclusively as a hydrogen bond donor in aqueous solution. The proton transfer reaction between the tautomeric forms was also investigated in the gas phase and in aqueous solution. Two distinct mechanisms were considered: a direct intramolecular transfer and a water-assisted mechanism. In the gas phase, the intramolecular transfer leads to a large energy barrier of 34.4 kcal/mol, passing through a three-center transition state. The proton transfer with the assistance of one water molecule decreases the energy barrier to 17.2 kcal/mol. In solution, these calculated activation barriers are, respectively, 32.0 and 14.8 kcal/mol. The solvent effect is found to be sizable but it is considerably more important as a participant in the water-assisted mechanism than the solvent field of the solute-solvent interaction. Finally, the calculated total Gibbs free energy is used to estimate the equilibrium constant.     

Theoretical investigation of the tautomerism of isoorotic acid in gaseous and aqueous phases

In the present investigation, the tautomeric and conformational equilibrium of isoorotic acid have been studied using Møller–Plesset second‐order (MP2) and density functional theory (DFT) methods in the gas phase and aqueous solution (ε = 78.5) using the IPCM model. The relative energies of these tautomers have been calculated at the two levels of theory using 6‐311++G** basis set. Energetics and relative stabilities of the tautomers were compared and analyzed in both the gaseous and aqueous phases. The results indicate that the diketo tautomer (iso) is the most stable form in the gas phase and water. The carboxylic substitution in the uracil ring does not alter its relative stability order of the tautomers. The proton affinity of the oxygen atoms and the deprotonation enthalpy of the NH bonds of isoorotic acid have been compared with recent data of uracil. The relative stability of both syn‐ and anti‐conformations was investigated and the syn form was found to be more stable by 17.65 kcal/mol. It was determined in ab initio calculations that an electron can attach to isoorotic acid, forming a stable anion better than uracil. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Comprehensive core-level study of the effects of isomerism, halogenation, and methylation on the tautomeric equilibrium of cytosine

Core-level X-ray photoemission and near-edge X-ray absorption fine structure spectra of 5-methylcytosine, 5-fluorocytosine, and isocytosine are presented and discussed with the aid of high-level ab initio calculations. The effects of the methylation, halogenation, and isomerization on the relative stabilities of cytosine tautomers are clearly identified spectroscopically. The hydroxy-oxo tautomeric forms of these molecules have been identified, and their quantitative populations at the experimental temperature are calculated and compared with the experimental results and with previous calculations. The calculated values of Gibbs free energy and Boltzmann population ratios are in good agreement with the experimental results characterizing tautomer equilibrium.     

Tautomerism and infrared spectra of 2-thiopurine: an experimental matrix isolation and theoretical ab initio and density functional theory study

Infrared spectra of 2-thiopurine (2-mercaptopurine, 2-purinethiol ) isolated in low-temperature Ar and N2 matrixes are reported. These spectra indicate that the compound adopts exclusively the thiol N9H tautomeric form. The theoretical calculations of relative energies of 2-thiopurine tautomers have been carried out at the MP4(SDTQ)//HF level using the 6-31G(d,p) basis set. The thiol N9H tautomer was predicted to be the most stable of all isomers of 2-thiopurine. The infrared spectra of the tautomers of 2-thiopurine have been calculated at the DFT(B3LYP)/6-31G(d,p) level. Good agreement between the experimental spectra and the spectra calculated for thiol N9H tautomer supported the identification of the dominant tautomer. It has also allowed for the reliable assignment of the bands observed in the experimental IR spectrum.     

Benzoderivatives of nucleic acid bases as modified DNA building blocks

The tautomeric properties of benzoderivatives of the canonical nucleic acid bases have been studied by using different computational approaches. Attention has been paid to the impact of the benzene group in altering the tautomeric preferences of the canonical bases both in the gas phase and in aqueous solution. To this end, relative solvation free energies of the tautomers determined from Self-Consistent Reaction Field continuum calculations and Monte Carlo-Free Energy Perturbation are combined with gas-phase tautomerization free energies determined from quantum mechanical calculations. The results provide a detailed picture of the tautomeric preferences of the benzoderivatives of nucleic acid bases. This information is used to examine the recognition properties of the preferred tautomers of the benzo-fused derivatives, paying particular attention to the ability to form Watson-Crick hydrogen-bonding and stacking interactions as well as to the hydrophobic nature of the modified bases. The implications of present results on the potential use of benzo-fused bases as potential building blocks in modified DNA duplexes are examined.     

Determination of thermodynamic parameters of tautomerization in gas phase by mass spectrometry and DFT calculations: Keto-enol versus nitrile-ketenimine equilibria

The study of tautomerics equilibria is really important because the reactivity of each compound with tautomeric capacity can be determined from the proportion of each tautomer. In the present work the tautomeric equilibria in some γ,δ-unsaturated β-hydroxynitriles and γ,δ-unsaturated β-ketonitriles were studied. The first family of compounds presents two possible theoretical tautomers, nitrile and ketenimine, while the second one presents four possible theoretical tautomers, keto-nitrile, enol (E and Z)-nitrile and keto-ketenimine. The equilibrium in gas phase was studied by gas chromatography-mass spectrometry (GC-MS). Tautomerization enthalpies were calculated by this methodology, and results were compared with those obtained by density functional theory (DFT) calculations, observing a good agreement between them. Nitrile tautomers were favored within the first family of compounds, while keto-nitrile tautomers were favored in the second family.     

Fitting quantum machine learning potentials to experimental free energy data: predicting tautomer ratios in solution

The computation of tautomer ratios of druglike molecules is enormously important in computer-aided drug discovery, as over a quarter of all approved drugs can populate multiple tautomeric species in solution. Unfortunately, accurate calculations of aqueous tautomer ratios—the degree to which these species must be penalized in order to correctly account for tautomers in modeling binding for computer-aided drug discovery—is surprisingly difficult. While quantum chemical approaches to computing aqueous tautomer ratios using continuum solvent models and rigid-rotor harmonic-oscillator thermochemistry are currently state of the art, these methods are still surprisingly inaccurate despite their enormous computational expense. Here, we show that a major source of this inaccuracy lies in the breakdown of the standard approach to accounting for quantum chemical thermochemistry using rigid rotor harmonic oscillator (RRHO) approximations, which are frustrated by the complex conformational landscape introduced by the migration of double bonds, creation of stereocenters, and introduction of multiple conformations separated by low energetic barriers induced by migration of a single proton. Using quantum machine learning (QML) methods that allow us to compute potential energies with quantum chemical accuracy at a fraction of the cost, we show how rigorous relative alchemical free energy calculations can be used to compute tautomer ratios in vacuum free from the limitations introduced by RRHO approximations. Furthermore, since the parameters of QML methods are tunable, we show how we can train these models to correct limitations in the underlying learned quantum chemical potential energy surface using free energies, enabling these methods to learn to generalize tautomer free energies across a broader range of predictions.

We show how alchemical free energies can be calculated with QML potentials to identify deficiencies in RRHO approximations for computing tautomeric free energies, and how these potentials can be learned from experiment to improve prediction accuracy.     

Role of the tautomerism of 2-azaadenine and 2-azahypoxanthine in substrate recognition by xanthine oxidase

The tautomerism of 2-azaadenine and 2-hypoxanthine has been examined in the gas phase and in aqueous solution. The tautomerism in the gas phase has been studied by means of semiempirical and ab initio quantum-mechanical computations, as well as density-functional calculations. The influence of the aqueous solvent on the relative stability between tautomers has been estimated from self-consistent reaction field calculations performed with different high-level continuum models. The results provide a detailed picture of the tautomeric preference for these purine bases. The importance of tautomerism in the substrate recognition by xanthine oxidase is discussed. Finally, the rate of oxidation of 2-azaadenine and 2- hypoxanthine by xanthine oxidase is discussed in terms of the recognition model at the enzyme active site.     

Tautomeric equilibria in 8-oxopurines: Implications for mutagenicity

Free-radical-induced DNA damage by ionizing radiation leads to a number of oxidized purines, of which 7H-8-oxoguanine (8OG) and 7H-8-oxoadenine (8OA) are predominant and known to cause an appreciable amount of cellular damage. A detailed quantum mechanical study at various levels of theory in both the gas phase and in an aqueous solution has been carried out in order to assess the tautomeric preferences of the bases. The calculated energies of various plausible tautomers suggest that at higher levels of ab initio theory with inclusion of electron correlation, the 8-keto-6-enolic form of 8-oxoguanine (8OG2) would predominate over the 6,8-diketo form (8OG1) in the gas phase whereas the 6-amino-8-keto form (8OA1) predominates over the other possible tautomers of 8-oxoadenine. Aqueous solvation, however, changes the gas-phase order for 8-oxoguanine, 8OG1 turning out to be the major tautomeric species in an aqueous medium. The estimated free energies of hydration by polarized continuum models are indicative that the mutagenically significant amounts of minor tautomeric forms of 8-oxoguanine and 8-oxoadenine exist in the aqueous phase and might be held responsible for inducing transversional as well as transitional mutations.     

Standard enthalpies of formation in the gas phase and relative stability of tautomers of C-nitro-1,2,4-triazole and isomers of N-alkyl-C-nitro-1,2,4-triazole: quantum-chemical studies

Quantum-chemical calculations of the standard enthalpies of formation in the gas phase of C-nitro-1,2,4-triazole and isomers of N-alkyl-C-nitro-1,2,4-triazoles (Alk = Me, Et, i-Pr, t-Bu) were carried out by the B3LYP method using equations for the isodesmic reactions and isomerization reactions. The relative Gibbs free energies of tautomers and isomers in aqueous solution were calculated. For the tautomers of C-nitro-1,2,4-triazole the structural indexes of aromaticity were estimated and the electron population density of the Natural Bond Orbital was analyzed. The relative stabilities of the described tautomers and isomers in the gas phase and in solutions were discussed on the basis of the results of the calculations. Dedicated to Academician B. A. Trofimov in his 70^th jubilee. Translated from Khimiya Geterosiklicheskikh Soedinenii, No. 1, 83-94, January, 2009.