Theoretical studies on the reaction of pentafulvenone with water in gas phase and aqueous solvent

In gas phase, the hydrations of pentafulvenone to generate three types of cyclopentadienyl carboxylic acids are studied theoretically at the MP2/6-311+G**//B3LYP/6-311+G** level. A water molecule attacking the C=O double bond of pentafulvenone can yield cyclopentadienyl carboxylic acids via the formation of fulvenediols, and attacking the C=C double bond of pentafulvenone can directly yield cyclopentadienyl carboxylic acid. The barriers of rate-determining transition states are 42.2 and 30.4 kcal mol⁻¹, respectively. The barriers of rate-determining transition states for two water molecules system are 20.2 and 19.6 kcal mol⁻¹, respectively. The products can isomerize to each other. In aqueous solvent, the hydrations of pentafulvenone are investigated using PCM-UAHF model at the MP2 (PCM)/6-311+G**// B3LYP (PCM)/6-311+G** and MP2 (PCM)/6-311+G**// B3LYP/6-311+G** levels. The barriers of all rate-determining transition states are decreased. The added water molecule acts as catalyst in both gas phase and aqueous solvent. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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

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

FTIR, Raman spectra and ab initio calculations of 2-mercaptobenzothiazole

FTIR and Raman spectra of a rubber vulcanization accelerator, 2-mercaptobenzothiazole (MBT), were recorded in the solid phase. The harmonic vibrational wavenumbers, for both the toutomeric forms of MBT, as well as for its dimeric complex, have been calculated, using ab initio RHF and density functional B3LYP methods invoking different basis sets upto RHF/6-31G** and B3LYP/6-31G** and the results were compared with the experimental values. Conformational studies have been also carried out regarding its toutomeric monomer forms and its dimer form. With all the basis sets the thione form of MBT (II) is predicted to be more stable than thiol form (I) and dimeric conformation (III) is predicted to be more stable with monomeric conformations (I) and (II). Vibrational assignments have been made, and it has been found that the calculated normal mode frequencies of dimeric conformation (III) are required for the analysis of IR and Raman bands of the MBT. The predicted shift in NH- stretching vibration towards the lower wave number side with the B3LYP/6-31G** calculations for the most stable dimer form (III), is in better agreement with experimental results. The intermolecular sulfur-nitrogen distance in N-H...S hydrogen bond was found to be 3.35 angstroms from these calculations, is also in agreement to the experimental value.     

Conformational behavior of some hydroxamic acids

The conformational preferences of a few hydroxamic acids are investigated by the density functional B3LYP/6-311++G**//B3LYP16-31G* and semiempirical AM1 and PM3 methods in this work. It is found that both semiempirical methods give satisfactory results in comparison with sophisticated DFT and ab initio calculations, except for the activation barriers, which are overestimated. Of the two semiempirical methods, while the PM3 method gives better results for relative stabilities, AM1 geometries are in slightly better agreement with the experiments. The keto forms are found to be most stable and the reaction pathways for the interconversion between the keto and enol forms have been deduced. The effect of solvation on the reaction has also been investigated, as has the effect of methyl substitution at the carbon and nitrogen atoms. All the investigated acids exhibit N-acid behavior.     

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

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

Experimental and theoretical structural studies on 4-(2-phenylethyl)-5-(2-furyl)-2,4-dihydro-3H-1,2,4-triazole-3-thione

The structural and conformational features of 4-(2-phenylethyl)-5-(2-furyl)-2, 4-dihydro-3H-1,2,4-triazole-3-thione (1a), which can be related to the biological activity, have been investigated by X-ray diffraction and molecular modeling techniques. Ab initio method (RHF/6-31G) and density functional theory (B3LYP/6-31G(D)) have been used to calculate structural parameters, conformations, and relative energy of two tautomeric specious (1a and 1b) of the title compound. The geometry and the conformation of the thione form, 1a, is well reproduced by the DFT (B3LYP/6-31G(D)) method as compared with X-ray structure in which this form is found. The thione form is also predicted to be 14.42 kcal/mol more stable than the thiol form in the gas-phase by the DFT method.     

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.     

A theoretical study of medium effects on the structure of the glycine analogue aminomethylphosphonic acid

alpha-Aminophosphonic acids are analogues of natural alpha-aminoacids and very promising agents for use in various pharmaceutical applications. However, in contrast to the numerous theoretical investigations on the structure of natural alpha-aminoacids, only very few studies on alpha-aminophosphonic acids have been performed. In the present work, we report a detailed investigation of the simplest compound, the glycine analogue aminomethylphosphonic acid (AMPA), by means of quantum mechanical calculations at the B3LYP/6-311++G(3df,2p)//B3LYP/6-31+G(d,p) and MP2/6-311++G(3df,2p)//B3LYP/6-31+G(d,p) levels. We focus on the structure of the neutral species looking at the evolution of non-ionized and ionized forms from gas phase to non-polar solvents and aqueous media. Continuum and discrete-continuum solvent models have been employed to account for the effects of the environment. The discussion is centered on: (1) the geometry and relative stability of possible conformers in gas phase and aqueous solution, (2) the free energy of tautomerization in different media, (3) the role of hydrogen bonds in liquid water, and (4) the free energy of transfer from water to a hydrophobic solvent such as cyclohexane. Systematic comparison between AMPA and Gly is performed. Though both systems exhibit many similarities, some important differences have also been found that may be explained, at least in part, by the higher acidity of phosphonic acids compared to carboxylic acids. In particular, in solvents lacking hydrogen-bond formation capability, Gly derivatives should mainly exist as non-ionized molecules while the equivalent AMPA derivatives should adopt a zwitterionic structure in media with dielectric constant above 10. This might have significant environmental or biological consequences that will need to be addressed.     

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

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

亚硝基化合物与甲醛的反应机理和溶剂效应的理论研究

采用密度泛函方法B3LYP/6-311++G(d,p)研究了亚硝基苯C₆H₅-NO和2-甲基-2-亚硝基丙烷(CH₃)₃C-NO与甲醛分别在气相和溶剂中的反应机理. 在气相中均找到两条反应通道, 即协同机理和分步机理, 均生成实验产物氧肟酸, 而且分步机理均为优势通道; 除2-甲基-2-亚硝基丙烷的反应没有协同途径外, 在溶剂中反应机理与气相中的类似. 采用导电极化连续介质模型分别研究了在乙腈与水溶液中反应的溶剂化效应, 发现这些溶剂可降低反应的活化能, 但降低的程度比较小, 反应速率变化不大.     

Hydroxylation Reaction Mechanism for Nitrosodimethylamine by Oxygen Atom

The hydroxylation reaction mechanism of nitrosodimethylamine(NDMA) by oxygen atom was theoretically investigated at the B3LYP/6-31G** level. It has been found that the path of the oxydation of the C―H bond is easier than the path involving a Singlet/Triplet crossing. The study of the potential surface shows that both solvent effect at B3LYP/6-31G** level and different method at more credible MP2/6-311G** level in the gas phase have no effect on the hydroxylation reaction mechanism. The oxidation hydroxylation process of NDMA by O is exothermic reaction and easy to occur.     

Electronic structure and reactivity of guanylthiourea: A quantum chemical study

Electronic structure analysis of guanylthiourea (GTU) and its isomers has been carried out using quantum chemical methods. Two major tautomeric classes (thione and thiol) have been identified on the potential energy (PE) surface. In both the cases conjugation of pi‐electrons and intramolecular H‐bonds have been found to play a stabilizing role. Various isomers of GTU on its PE surface have been analyzed in two different groups (thione and thiol). The interconversion from the most stable thione conformer ( GTU‐1 ) to the most stable thiol conformer ( GTU‐t1 ) was found to take place via bimolecular process which involves protonation at sulfur atom of GTU‐1 followed by subsequent C? N bond rotation and deprotonation. The detailed analysis of the protonation has been carried out in gas phase and aqueous phase (using CPMC model). Sulfur atom (S1) was found to be the preferred protonation site (over N4) in GTU‐1 in gas phase whereas N4 was found to be the preferred site of protonation in aqueous medium. The mechanism of S‐alkylation reaction in GTU has also been studied. The formation of alkylated analogs of thiol isomers (alkylated guanylthiourea) is believed to take place via bimolecular process which involves alkyl cation attack at S atom followed by C? N bond rotation and deprotonation. The reactive intermediate RS(NH₂)C? N? C(NH₂)₂⁺ belongs to the newly identified ^⊕N(←L)₂ class of species and provides the necessary dynamism for easy conversion of thione to thiol. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

可乐定分子构象异构和互变异构的理论研究

采用HF/6-31G*方法, 对氨基型和亚胺型可乐定中性分子以及质子化的可乐定分子进行了构象分析. 之后采用B3LYP/6-31G**方法对处于势能面上的能量极小构象异构体进行全自由度几何优化和频率分析. Onsager反应场溶剂模型用于水相的计算. 结果表明, 在气相和水相中可乐定中性分子主要以亚胺型存在, 这同已有实验结果一致. 进一步, 寻找构象异构化过渡态和氨基型亚胺型互变异构化过渡态, 探讨质子化和水溶剂化效应对异构化过程的几何结构和能量的影响. 为了考察氯苯的共轭效应对可乐定互变异构体稳定性的影响, 在B3LYP/6-31G**水平上, 研究了2-氨基-2-咪唑啉的互变异构化反应机理.     

An ab initio and density functional theory study of keto-enol equilibria of hydroxycyclopropenone in gas and aqueous solution phase

Keto-enol tautomerism in hydroxycyclopropenone (2-hydroxy-2-cyclopropen-1-one) has been studied using ab initio methods, the B3LYP functional of density functional theory, as well as complete basis set (CBS-QB3 and CBS-APNO) and G3 methods. Absolute and relative energies were calculated with each of the methods, whereas computations of geometries and harmonic frequencies for hydroxycyclopropenone and 1,2-cyclopropanedione were computed in the gas phase but were limited to HF, MP2 and CCSD levels of theory, and the B3LYP functional, in combination with the 6-31++G** basis set. Using the MP2/6-31++G** gas phase optimized structure, each species was then optimized fully in aqueous solution by employing the polarizable continuum model (PCM) self-consistent reaction field approach, in which HF, MP2 and B3LYP levels of theory were utilized, with the same 6-31++G** basis set. In both gas and aqueous solution phases, the keto form is higher in energy for all of the model chemistries considered. The presence of the solvent, however, is found to have very little effect on the bond lengths, angles and harmonic frequencies. From the B3LYP/6-31++G** Gibbs free energy, the keto-enol tautomeric equilibrium constant for 2-hydroxy-2-cyclopropen-1-one <==> 1,2-cyclopropanedione is computed to be K(T)(gas) = 2.35 x 10(-6), K(T)(aq) = 5.61 x 10(-14). It is concluded that the enol form is overwhelmingly predominant in both environments, with the effect of the solvent shifting the direction of equilibrium even more strongly in the favor of hydroxycyclopropenone. The almost exclusive nature of this species is attributed to stabilization resulting from aromaticity. Confirmation is provided by comparison of the simulated vibrational spectra of hydroxycyclopropenone with the measured infrared spectrum in an argon matrix.     

Tautomeric equilibria of 2‐ and 4‐thiouracil in gas phase and in solvent: A density functional study

The relative stabilities of the five favored tautomers of 2‐ and 4‐thiouracil in gas phase and in water solution were determined by density functional theory employing the Becke, Lee, Yang, and Parr (B3LYP) exchange–correlation potential and the three 6‐31G(d,p), 6‐311++G(d,p), and triple‐zeta valence (TZVP) basis sets. Zero‐point vibrational corrections were also computed. Bulk solvent effects were studied in the framework of the self‐consistent reaction field approach by the polarizable continuum model. All calculations indicate that the most stable tautomer for both species, in the gas phase as well as in solution, has the oxo‐thione form, in full agreement with the previous ab initio and experimental studies. The tautomeric stability orders obtained in the aqueous solution are sensibly different from that in the gas phase. At B3LYP/6‐311++G(d,p) level in the gas phase, the following orders of stability for 2‐ and 4‐thiouracil tautomers were observed, respectively: S2U1>S2U2>S2U4>S2U5>S2U3 and S4U1>S4U2>S4U3>S4U4>S4U5. The corresponding trends in the aqueous phase are S2U1>S2U3>S2U2>S2U5>S2U4 and S4U1>S4U2>S4U3>S4U5>S4U4. On the basis of the computed energy differences we can hypothesize that only the oxo‐thione forms of 2‐ and 4‐thiouracil should exist in the gas phase and in water solution. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 82: 44–52, 2001     

Berberine alkaloid: Quantum chemical study of different forms by the DFT and MP2 methods

The stable structures and electronic properties for the berberine cation as well as possible ammonium, carbinol and amino-aldehyde forms of protoberberine salts in the presence of hydroxyl ions were investigated by the B3LYP/6-31G(d,p) and MP2/6-31++G(d,p) methods. The geometry optimizations by both methods lead to the nonplanar propeller-twisted and buckled structure for the all forms. The obtained bond lengths and bond angles agree with the experimental values. The comparison of total energies elucidates that the amino-aldehyde form is the most preferable tautomer in gas phase, while the carbinol form is less stable. The least stable tautomer is the ammonium form.     

Low-energy electron attachment to 5'-thymidine monophosphate: modeling single strand breaks through dissociative electron attachment

Mechanisms of low-energy electron (LEE) attachment and subsequent single-strand break (SSB) formation are investigated by density functional theory treatment of a simple model for DNA, i.e., the nucleotide, 5'-thymidine monophosphate (5'-dTMPH). In the present study, the C5'-O5' bond dissociation due to LEE attachment has been followed along the adiabatic as well as on the vertical (electron attached to the optimized geometry of the neutral molecule) anionic surfaces using B3LYP functional and 6-31G* and 6-31++G** basis sets. Surprisingly, it is found that the PES of C5'-O5' bond dissociation in the anion radicals have approximately the same barrier for both adiabatic and vertical pathways. These results provide support for the hypothesis that transiently bound electrons (shape resonances) to the virtual molecular orbitals of the neutral molecule likely play a key role in the cleavage of the sugar-phosphate C5'-O5' bond in DNA resulting in the direct formation of single strand breaks without significant molecular relaxation. To take into account the solvation effects, we considered the neutral and anion radical of 5'-dTMP surrounded by 5 or 11 water molecules with Na+ as a counterion. These structures were optimized using the B3LYP/6-31G** level of theory. We find the barrier height for adiabatic C5'-O5' bond dissociation of 5'-dTMP anion radical in aqueous environment is so substantially higher than in the gas phase that the adiabatic route will not contribute to DNA strand cleavage in aqueous systems. This result is in agreement with experiment.     

2-巯基苯并咪唑及其类似物互变异构的理论研究

采用B3LYP/6-311G**方法, 计算了2-巯基苯并咪唑及其类似物(2-巯基苯并噁唑、2-巯基苯并噻唑、2-羟基苯并咪唑、2-羟基苯并噁唑、2-羟基苯并噻唑以及2-巯基咪唑、2-巯基噁唑、2-巯基噻唑、2-羟基咪唑、2-羟基噁唑、2-羟基噻唑)的(硫)醇式与(硫)酮式结构进行质子迁移的3种可能途径: (a)分子内质子迁移; (b)水助质子迁移; (c)甲醇助质子迁移．结果表明, 途经b和c所需要的活化能较小, 氢键在降低反应活化能方面起重要作用．采用PCM方法研究了反应体系的溶剂化效应．结果表明孤立分子、一水合物和一甲醇合物的最稳定异构体相同, 都为(硫)酮式, 与气相结论一致．溶剂化效应对异构化能垒的影响较小.     

2-Pyridinethiol/2-pyridinethione tautomeric equilibrium. A comparative experimental and computational study

The gas phase and solvent dependent preference of the tautomerization between 2-pyridinethiol (2SH) and 2-pyridinethione (2S) has been assessed using variable temperature Fourier transform infrared (FTIR) experiments, as well as ab initio and density functional theory computations. No spectroscopic evidence (nu(S)(-)(H) stretch) for 2SH was observed in toluene, C(6)D(6), heptane, or methylene chloride solutions. Although, C(s)() 2SH is 2.61 kcal/mol more stable than C(s)() 2S (CCSD(T)/cc-pVTZ//B3LYP/6-311+G(3df,2p)+ZPE), cyclohexane solvent-field relative energies (IPCM-MP2/6-311+G(3df,2p)) favor 2S by 1.96 kcal/mol. This is in accord with the FTIR observations and in quantitative agreement with the -2.6 kcal/mol solution (toluene or C(6)D(6)) calorimetric enthalpy for the 2S/2SH tautomerization favoring the thione. As the intramolecular transition state for the 2S, 2SH tautomerization (2TS) lies 25 (CBS-Q) to 30 kcal/mol (CCSD/cc-pVTZ) higher in energy than either tautomer, tautomerization probably occurs in the hydrogen bonded dimer. The B3LYP/6-311+G(3df,2p) optimized C(2) 2SH dimer is 10.23 kcal/mol + ZPE higher in energy than the C(2)(h)() 2S dimer and is only 2.95 kcal/mol + ZPE lower in energy than the C(2) 2TS dimer transition state. Dimerization equilibrium measurements (FTIR, C(6)D(6)) over the temperature range 22-63 degrees C agree: K(eq)(298) = 165 +/- 40 M(-)(1), DeltaH = -7.0 +/- 0.7 kcal/mol, and DeltaS = -13.4 +/- 3.0 cal/(mol deg). The difference between experimental and B3LYP/6-311+G(3df,2p) [-34.62 cal/(mol deg)] entropy changes is due to solvent effects. The B3LYP/6-311+G(3df,2p) nucleus independent chemical shifts (NICS) are -8.8 and -3.5 ppm 1 A above the 2SH and 2S ring centers, respectively, and the thiol is aromatic. Although the thione is not aromatic, it is stabilized by the thioamide resonance. In solvent, the large 2S dipole, 2-3 times greater than 2SH, favors the thione tautomer and, in conclusion, 2S is thermodynamically more stable than 2SH in solution.     

密度泛函研究细梗胡枝子提取物中新型抗氧化剂的分子结构和活性氧清除机理

The molecular structure and radical scavenging activity of three novel antioxidants from Lespedeza Virgata, lespedezavirgatol, lespedezavirgatal, and lespedezacoumestan, have been studied using density functional theory with the B3LYP and BhandHLYP methods. The optimized geometries of neutral, radical cation, radical and anion forms were obtained at the B3LYP/6-31G(d) level, in which it was found that all the most stable conformations contain intramolecular hydrogen bonds. The same results were obtained from the MP2 method. The homolytic O-H bond dissociation enthalpy and the adiabatic ionization potential of neutral and anion forms for the three new antioxidants and adiabatic electron affinity and H-atom affinity for hydroxyl radical, superoxide anion radical, and hydrogen peroxide radical were determined both in gas phase and in aqueous solution using IEF-PCM and CPCM model with UAHF or Bondi cavity. The antioxidant activities and reactive oxygen species scavenging mechanisms were then discussed, and the results obtained from different methods are consistent. Furthermore, the antioxidant activities are consistent with the experimental findings of the compounds under investigation.