A density functional theory study on pelargonidin

A complete conformational analysis on the isolated and polarizable continuum model (PCM) modeled aqueous solution cation, quinonoidal, and anion forms of pelargonidin, comprising the diverse tautomers of the latter forms, was carried out at the B3LYP/6-31++G(d,p) level. The results indicate that the most stable conformer of cationic and quinonoidal forms of pelargonidin are completely planar in the gas phase, whereas that of the anionic form is not planar. In contrast, PCM calculations show that the plane of the B ring is slightly rotated with regard to the AC bicycle in the most stable conformer of the cation and quinonoidal form. The most stable conformers of the cation, both in gas phase and aqueous solution, display anti and syn orientations for, respectively, C2-C3-O-H and C6-C5-O-H dihedral angles, whereas syn and anti orientation of hydroxyls at 7 and 4' positions are nearly isoenergetic. The most stable tautomer of quinonoidal pelargonidin is obtained by deprotonating hydroxyl at C5 in gas phase but at C7 according to PCM. Also, the most stable tautomer of the anion is different in gas phase (hydrogens are abstracted from hydroxyls at C5 and C4') and PCM simulation (C3 and C5). Tautomeric equilibria affect substantially the geometries of the AC-B backbone providing bond length variations that basically agree with the predictions of the resonance model. Most of the conformers obtained display an intramolecular hydrogen bond between O3 and H6'. Nevertheless, this interaction is not present in the most stable anions. Ionization potentials and O-H bond dissociation energies computed for the most stable conformers of cation, quinonoidal, and anion forms are consistent with an important antioxidant activity.     

On the unusual stability of valence anions of thymine based on very rare tautomers: A computational study

We characterized anionic states of thymine using various electronic structure methods, with the most accurate results obtained at the CCSD(T)/aug-cc-pVDZ level of theory followed by extrapolations to complete basis set limits. We found that the most stable anion in the gas phase is related to an imino-oxo tautomer, in which the N1H proton is transferred to the C5 atom. This valence anion, aT(c5)(nl), is characterized by an electron vertical detachment energy (VDE) of 1251 meV and it is adiabatically stable with respect to the canonical neutral nT(can) by 2.4 kcal/mol. It is also more stable than the dipole-bound (aT(dbs)(can)), and valence anion aT(val)(can) of the canonical tautomer. The VDE values for aT(dbs)(can)and T(val)(can) are 55 and 457 meV, respectively. Another, anionic, low-lying imino-oxo tautomer with a VDE of 2458 meV has a proton transferred from N3H to C5 aT(c5)(n3). It is less stable than aT(val)(can) by 3.3 kcal/mol. The mechanism of formation of anionic tautomers with the carbons C5 or C6 protonated may involve intermolecular proton transfer or dissociative electron attachment to the canonical neutral tautomer followed by a barrier-free attachment of a hydrogen atom to C5. The six-member ring structure of the anionic tautomers with carbon atoms protonated is unstable upon an excess electron detachment. Within the PCM hydration model, the low-lying valence anions become adiabatically bound with respect to the canonical neutral; becomes the most stable, being followed by aT(c5)(nl), aT(c5)(n3), aT(can), and aT(c5)(nl).     

Stabilization of very rare tautomers of 1-methylcytosine by an excess electron

We characterized valence anionic states of 1-methylcytosine using various electronic structure methods. We found that the most stable valence anion is related to neither the canonical amino-oxo nor a rare imino-oxo tautomer, in which a proton is transferred from the N4 to N3 atom. Instead, it is related to an imino-oxo tautomer, in which the C5 atom is protonated. This anion is characterized by an electron vertical detachment energy (VDE) of 2.12 eV and it is more stable than the anion based on the canonical tautomer by 1.0 kcal/mol. The latter is characterized by a VDE of 0.31 eV. Another unusual low-lying imino-oxo tautomer with a VDE of 3.60 eV has the C6 atom protonated and is 3.6 kcal/mol less stable than the anion of the canonical tautomer. All these anionic states are adiabatically unbound with respect to the canonical amino-oxo neutral, with the instability of 5.8 kcal/mol for the most stable valence anion. The mechanism of formation of anionic tautomers with carbon atoms protonated may involve intermolecular proton transfer or dissociative electron attachment to the canonical neutral tautomer followed by a barrier-free attachment of a hydrogen atom to the C5 or C6 atom. The six-member ring structure of anionic tautomers with carbon atoms protonated is unstable upon an excess electron detachment. Indeed the neutral systems collapse without a barrier to a linear or a bicyclo structure, which might be viewed as lesions to DNA or RNA. Within the PCM hydration model, the anions become adiabatically bound with respect to the corresponding neutrals, and the two most stable tautomers have a carbon atom protonated.     

Computational Study of Proton Transfer in Tautomers of 3‐ and 5‐Hydroxypyrazole Assisted by Water

The tautomerism of 3‐ and 5‐hydroxypyrazole is studied at the B3LYP, CCSD and G3B3 computational levels, including the gas phase, PCM–water effects, and proton transfer assisted by water molecules. To understand the propensity of tautomerization, hydrogen‐bond acidity and basicity of neutral species is approached by means of correlations between donor/acceptor ability and H‐bond interaction energies. Tautomerism processes are highly dependent on the solvent environment, and a significant reduction of the transition barriers upon solvation is seen. In addition, the inclusion of a single water molecule to assist proton transfer decreases the barriers between tautomers. Although the second water molecule further reduces those barriers, its effect is less appreciable than the first one. Neutral species present more stable minima than anionic and cationic species, but relatively similar transition barriers to anionic tautomers.     

Valence and dipole-bound anions of the most stable tautomers of guanine

Anionic states of guanine, which is the only nucleic acid base of which the anions have not yet been studied in either photoelectron spectroscopic (PES) or Rydberg electron transfer (RET) experiments, have been characterized for the four most stable tautomers of neutral guanine using a broad spectrum of electronic structure methods from the density functional theory, with the B3LYP exchange-correlation functional, to the coupled-cluster method, with single, double, and perturbative triple excitations. Both valence and dipole-bound anionic states were addressed. We identified some of the difficulties facing future PES or RET experiments on the anion of guanine. Even if guanine is successfully transferred to the gas phase without thermal decomposition, it is critical to have the canonical amino-oxo (G) and both amino-hydroxy (GH and GHN7H) tautomers in the beam, not only the most stable, a noncanonical, amino-oxo tautomer (GN7H), as the latter does not support an adiabatically bound anionic state. We also suggested a scheme for enrichment of gas-phase guanine with the canonical tautomer, which is not the most stable in the gas phase, but which is of main interest due to its biological relevance. The tautomers G, GN7H, and GHN7H support vertically bound valence anionic states with the CCSD(T) value of vertical detachment energy of +0.58, +0.21, and +0.39 eV, respectively. These anionic states are, however, adiabatically unbound and thus metastable. The vertical electronic stability of these valence anionic states is accompanied by serious "buckling" of the molecular skeleton. The G and GHN7H tautomers support dipole-bound states with the CCSD(T) values of adiabatic electron affinity of 65 and 36 meV, respectively. A contribution from higher-than-second-order correlation terms represents, respectively, 48 and 68% of the total vertical electron detachment energy determined at the CCSD(T) level.     

Barrier-free proton transfer in the valence anion of 2'-deoxyadenosine-5'-monophosphate. II. A computational study

The propensity of four representative conformations of 2(')-deoxyadenosine-5(')-monophosphate (5(')-dAMPH) to bind an excess electron has been studied at the B3LYP6-31++G(d,p) level. While isolated canonical adenine does not support stable valence anions in the gas phase, all considered neutral conformations of 5(')-dAMPH form adiabatically stable anions. The type of an anionic 5(')-dAMPH state, i.e., the valence, dipole bound, or mixed (valence/dipole bound), depends on the internal hydrogen bond(s) pattern exhibited by a particular tautomer. The most stable anion results from an electron attachment to the neutral syn-south conformer. The formation of this anion is associated with a barrier-free proton transfer triggered by electron attachment and the internal rotation around the C4(')-C5(') bond. The adiabatic electron affinity of the a_south-syn anion is 1.19 eV, while its vertical detachment energy is 1.89 eV. Our results are compared with the photoelectron spectrum (PES) of 5(')-dAMPH(-) measured recently by Stokes et al., [J. Chem. Phys. 128, 044314 (2008)]. The computational VDE obtained for the most stable anionic structure matches well with the experimental electron binding energy region of maximum intensity. A further understanding of DNA damage might require experimental and computational studies on the systems in which purine nucleotides are engaged in hydrogen bonding.     

Solvent effects on glycine II. Water-assisted tautomerization

The water-assisted tautomerization of glycine has been investigated at the B3LYP/6-31+G** level using supermolecules containing up to six water molecules as well as considering a 1:1 glycine-water complex embedded in a continuum. The conformations of the tautomers in this mechanism do not display an intramolecular H bond, instead the functional groups are bridged by a water molecule. The replacement of the intramolecular H bond by the bridging water reduces the polarity of the N-H bond in the zwitterion and increases that of the O-H bond in the neutral, stabilizing the zwitterion. Both the charge transfer effects and electrostatic interactions stabilize the nonintramolecularly H-bonded zwitterion conformer over the intramolecularly hydrogen bonded one. The nonintramolecularly H-bonded neutral is favored only by charge transfer effects. Although there is no strong evidence whether the intramolecularly hydrogen bonded or non hydrogen bonded structures are favored in the bulk solution represented as a dielectric continuum, it is likely that the latter species are more stable. The free energy of activation of the water-assisted mechanism is higher than the intramolecular proton transfer channel. However, when the presumably higher conformational energy of the zwitterion reacting in the intramolecular mechanism is taken into account, both mechanisms are observed to compete. The various conformers of the neutral glycine may form via multiple proton transfer reactions through several water molecules instead of a conformational rearrangement.     

Interaction with glycine increases stability of a mutagenic tautomer of uracil. A density functional theory study

The most stable structures for the gas-phase complexes of minor tautomers of uracil (U) with glycine (G) were characterized at the density functional B3LYP/6-31++G level of theory. These are cyclic structures stabilized by two hydrogen bonds. The relative stability of isolated tautomers of uracil was rationalized by using thermodynamic and structural arguments. The stabilization energies for complexes between the tautomers of U and G result from interplay between the stabilizing two-body interaction energies and destabilizing one-body terms. The latter are related to the energies of (i) tautomerization of the unperturbed moieties and (ii) distortions of the resulting rare tautomers in the complex. The two-body term describes the interaction energy between distorted tautomers. The two-body interaction energy term correlates with perturbations of length of the proton-donor bonds as well as with deprotonation enthalpies and proton affinities of the appropriate monomer sites. It was demonstrated that the relative instability of rare tautomers of uracil is diminished due to their interactions with glycine. In particular, the instability of the third most stable tautomer (U(III)) is decreased from 11.9 kcal/mol for non-interacting uracil to 6.7 kcal/mol for uracil in a complex with the zwitterionic tautomer of glycine. A decrease of instability by 5.2 kcal/mol could result in an increase of concentration of U(III) by almost 5 orders of magnitude. This is the tautomer with proton donor and acceptor sites matching guanine rather than adenine. Moreover, kinetic characteristics obtained for the glycine-assisted conversion of the most stable tautomer of uracil (U(I)) to U(III) indicate that the U(I)<-->U(III) thermodynamic equilibrium could be easily attained at room temperature. The resulting concentration of this tautomer falls in a mutationally significant range.     

Computational study on the bond dissociation enthalpies in the enolic and ketonic forms of beta-diketones: their influence on metal-ligand bond enthalpies

A computational study on the thermodynamic properties of 13 beta-diketones is presented. The B3LYP//6-311+G(2d,2p)//B3LYP/6-31G(d) theoretical approach was employed to compute the O-H and C-H bond dissociation enthalpies and enthalpy of tautomerization and to estimate standard gas-phase enthalpies of formation for the radicals and for the parent molecules. The gas-phase enthalpies of formation for the neutral molecules are in excellent agreement with available experimental data, supporting the estimates made for the radicals. The latter are very important for the clarification of the thermochemistry of many beta-diketonato metal complexes previously reported in the literature. Importantly, when substituents R = -CHR' are attached to the beta-diketone's scaffold, C-H homolytic bond cleavage is always favored with respect to O-H bond scission.     

Gas phase ion chemistry of charged silver(I) adenine ions via multistage mass spectrometry experiments and DFT calculations

Electrospray ionization (ESI) of solutions containing adenine and AgNO(3) yields polymeric [Ad(x)+ Ag(y)-zH]((y-z)+) species. Density functional theory (DFT) calculations have been used to examine potential structures for several of the smaller ions while multistage mass spectrometry experiments have been used to probe their unimolecular reactivity (via collision-induced dissociation (CID)) and bimolecular reactivity (via ion-molecule reactions with the neutral reagents acetonitrile, methanol, butylamine and pyridine). DFT calculations of neutral adenine tautomers and their silver ion adducts provide insights into the binding modes of adenine. We find that the most stable [Ad + Ag](+) ion does not correspond to the most stable neutral adenine tautomer, consistent with previous studies that have shown that transition metal ions can stabilize rare tautomeric forms of nucleobases. Both the charge and the stoichiometry of the [Ad(x)+ Ag(y)-zH]((y-z)+) complexes play pivotal roles in directing the types of fragmentation and ion-molecule reactions observed. Thus, [Ad(2)+ Ag(2)](2+) is observed to dissociate to [Ad + Ag](+) and to react with butylamine via proton transfer, while [Ad(2)+ Ag(2)- H](+) fragments via loss of neutral adenine to form the [Ad + Ag(2)- H](+) ion and does not undergo proton transfer to butylamine. DFT calculations on several isomeric [Ad(2)+ Ag(2)](2+) ions suggest that planar centrosymmetric cations, in which two adjacent silver atoms are bridged by two N7H adenine tautomers via N(3),N(9)-bidentate interactions, are the most stable. The [Ad + Ag(2)-H](+) ion adds two neutral reagents in ion-molecule reactions, consistent with the presence of two vacant coordination sites. It undergoes a silver atom loss to form the [Ad + Ag - H](+) radical cation, which in turn fragments quite differently to the even electron [Ad + Ag](+) ion. Several other pairs of radical cation/even electron adenine-silver complexes were also found to undergo different fragmentation reactions.     

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

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.     

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-二巯基嘌呤氢键一水合物稳定性、降低质子转移异构化反应活化能等方面起着重要的作用.     

尿酸分子互变异构体平面构象的理论研究

使用半经验量子化学中的AM1方法、从头计算Hartree-Fock理论(在3-21G*水平)和密度泛函理论中的B3LYP方法(使用6-31G(d)基组),研究了尿酸分子的所有35种互变异构体。计算结果表明,三羰基互变异构体是所有异构体中能量最低的,其次为单羟基异构体和双羟基异构体,而含有三羟基的互变异构体相对能量最高。随着羟基数的增加, C-N键的平均键长从1.395逐渐缩短到1.351,而CC键的平均键长基本保持不变(1.400~1.406)。     

Molecular energetics of cytosine revisited: a joint computational and experimental study

A static bomb calorimeter has been used to measure the standard molar energy of combustion, in oxygen, at T = 298.15 K, of a commercial sample of cytosine. From this energy, the standard (p degrees = 0.1 MPa) molar enthalpy of formation in the crystalline state was derived as -(221.9 +/- 1.7) kJ.mol(-1). This value confirms one experimental value already published in the literature but differs from another literature value by 13.5 kJ.mol(-1). Using the present standard molar enthalpy of formation in the condensed phase and the enthalpy of sublimation due to Burkinshaw and Mortimer [J. Chem. Soc., Dalton Trans. 1984, 75], (155.0 +/- 3.0) kJ.mol(-1), results in a value for the gas-phase standard molar enthalpy of formation for cytosine of -66.9 kJ.mol(-1). A similar value, -65.1 kJ.mol(-1), has been estimated after G3MP2B3 calculations combined with the reaction of atomization on three different tautomers of cytosine. In agreement with experimental evidence, the hydroxy-amino tautomer is the most stable form of cytosine in the gas phase. The enthalpies of formation of the other two tautomers were also estimated as -60.7 kJ.mol(-1) and -57.2 kJ.mol(-1) for the oxo-amino and oxo-imino tautomers, respectively. The same composite approach was also used to compute other thermochemical data, which is difficult to be measured experimentally, such as C-H, N-H, and O-H bond dissociation enthalpies, gas-phase acidities, and ionization enthalpies.     

Steady‐State Spectroscopy of the 2‐(N‐methylacetimidoyl)‐1‐naphthol Molecule

The steady‐state spectroscopy of 2‐(N‐methylacetimidoyl)‐1‐naphthol (MAN) reveals composite absorption and emission spectra from 298 to 193 K in hexane. The ground electronic state (So) absorption can be assigned to the sum of three molecular structures: the OH normal tautomer, and two NH proton transfer tautomers. The NH‐structures are the most stable ones in equilibrium with the OH tautomer for the S₀ state. On photoexcitation of the OH tautomer the excited state intramolecular proton transfer is undergone, and the corresponding NH emission is monitored at 470 nm. On photoexcitation of the NH tautomers the previous emission is monitored in addition to another emission at 600 nm, which is ascribed to intramolecular hydrogen‐bonded (IHB) nonplanar NH structures generated from the IHB planar NH tautomers. A Jab?oński diagram is introduced which gathers all the experimental evidence as well as the theoretical calculations executed at the DFT‐B3LYP and TD‐DFT levels. The MAN molecule is compared with other analogs such as 1‐hydroxy‐2‐acetonaphthone (HAN), 2‐(1?‐hydroxy‐2?‐naphthyl)benzimidazole and methyl 1‐hydroxy‐2‐naphthoate to validate the theoretical calculations. Photoexcitation of MAN generates two emission bands at longer wavelengths than that of the emission band of HAN. The MAN molecule exhibits a great photostability in hydrocarbon solution which depends on the photophysics of the NH tautomers (keto forms).     

Effect of cooperative hydrogen bonding in azo-hydrazone tautomerism of azo dyes

Azo-hydrazone tautomerism in azo dyes has been modeled by using density functional theory (DFT) at the B3LYP/6-31+G(d,p) level of theory. The most stable tautomer was determined both for model compounds and for azo dyes Acid Orange 7 and Solvent Yellow 14. The effects of the sulfonate group substitution and the replacement of the phenyl group with naphthyl on the tautomer stability and on the behavior in solvent have been discussed. Intramolecular hydrogen bond energies have been estimated for the azo and hydrazone tautomers to derive a relationship between the tautomer stability and the hydrogen bond strength. The transition structures for proton transfer displayed resonance assisted strong hydrogen bonding properties within the framework of the electrostatic-covalent hydrogen bond model (ECHBM). Evolution of the intramolecular hydrogen bond with changing structural and environmental factors during the tautomeric conversion process has been studied extensively by means of the atoms-in-molecules (AIM) analysis of the electron density. The bulk solvent effect was examined using the self-consistent reaction field model. Special solute-solvent interactions were further investigated by means of quantum mechanical calculations after defining the first-solvation shell by molecular dynamics simulations. The effect of cooperative hydrogen bonding with solvent molecules on the tautomer stability has been discussed.     

Spectral and acid-base features of 3,7-dihydroxy-2,8-diphenyl-4H,6H-pyrano[3,2-g]chromene-4,6-dione (diflavonol)--a potential probe for monitoring the properties of liquid phases

Diflavonol is a molecule that can exist in neutral or anionic form and in several tautomeric forms in ground and excited states. Absorption and emission spectroscopy combined with theoretical calculations have shown that only one tautomer of neutral diflavonol exists in the ground state, but two exist in the excited state. In the latter case, one is the tautomer originating from the ground state tautomer, which exists in strongly protic solvents, the other is the phototautomer occurring in weakly protic or aprotic solvents as a result of the intramolecular transfer of one proton. The OH groups present in diflavonol and involved in weak intramolecular hydrogen bonds exhibit a proton-donating ability reflected by the experimental values of acidity constants or theoretical enthalpies and free energies of proton detachment. The electronically excited molecule is a relatively strong acid when it loses one proton. With increasing basicity of the medium, monoanionic and dianionic forms occur which exhibit spectral characteristics and an emission ability different from those of neutral diflavonol. These interesting features of diflavonol open up possibilities for the analytical use of the compound and its application as a spectral probe sensitive to the properties of liquid phases.     

Tautomer selective photochemistry in 1-(tetrazol-5-yl)ethanol

A combined matrix isolation FTIR and theoretical DFT/B3LYP/6-311++G(d,p) study of the molecular structure and photochemistry of 1-(tetrazol-5-yl)ethanol [1-TE] was performed. The potential energy surface landscapes of the 1H and 2H tautomers of the compound were investigated and the theoretical results were used to help characterize the conformational mixture existing in equilibrium in the gas phase prior to deposition of the matrices, as well as the conformers trapped in the latter. In the gas phase, at room temperature, the compound exists as a mixture of 12 conformers (five of the 1H tautomer and seven of the 2H tautomer). Upon deposition of the compound in an argon matrix at 10 K, only three main forms survive, because the low barriers for conformational isomerization allow extensive conformational cooling during deposition. Deposition of the matrix at 30 K led to further simplification of the conformational mixture with only one conformer of each tautomer of 1-TE surviving. These conformers correspond to the most stable forms of each tautomer, which bear different types of intramolecular H-bonds: 1H-I has an NH···O hydrogen bond, whereas 2H-I has an OH···N hydrogen bond. Upon irradiating with UV light (λ > 200 nm), a matrix containing both 1H-I and 2H-I forms, an unprecedented tautomer selective photochemistry was observed, with the 2H tautomeric form undergoing unimolecular decomposition to azide + hydroxypropanenitrile and the 1H-tautomer being photostable.     

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.