Adiabatically bound valence anions of guanine

It is believed that guanine, a basic component of DNA and RNA, has the smallest affinity to an excess electron among all nucleic acid bases. Our experimental and computational findings indicate, however, that many so far neglected tautomers of guanine support adiabatically bound anionic states in the gas phase. The computed values of electron vertical detachment energy for the most stable anionic tautomers are within a broad range of the dominant feature of the photoelectron spectrum. We suggest that guanine might be the strongest excess electron acceptor among nucleic acid bases. Thus it might be critical to radiobiological damage of DNA and it might contribute to those chemical transformations of DNA that proceed through bound anionic states.     

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.     

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.     

Electron-rich three-center bonding: role of s,p interactions across the p-block

This paper analyzes the importance of s,p mixing-a necessary addition to the simplest Rundle-Pimentel picture-and periodic and group trends in electron-rich three-center bonding. Our analysis proceeds through a detailed quantum chemical study of the stability of electron-rich three-center bonding in triatomic 22-valence electron anions. To provide interpretations, a perturbational molecular orbital (MO) analysis of s,p mixing is carried out. This analysis of the orbitals and the overlap populations is then tested by density functional calculations for a number of linear trihalides, trichalcogenides, and tripnictides. The most important effect of s,p mixing on the in-line bonding is in destabilization of the 3sigma(g) orbital and is determined by the overlap between the s orbital of the central atom and the p orbital of the terminal atom. Further destabilization arises from the repulsion of p(pi) lone pairs. Both of these antibonding effects increase with increasing negative charge of the system. The stability of isoelectronic X(3) systems thus decreases when moving from right to left in the periodic table. Interesting group trends are discerned; for instance, for the electron-rich tripnictides, the ability to accommodate a hypervalent electron count is the largest in the middle rather at the end of the group. Particularly strong s,p mixing can reverse the bonding/antibonding character of MOs: thus MO 2sigma(u) that is responsible for bonding for trihalides and trichalcogenides is actually antibonding in N(3)(7)(-).     

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).     

顺铂化合物与鸟嘌呤异构体相互作用的理论研究

The influence of binding of cisplatin adducts on tautomeric equilibrium of guanine was investigated using quantum chemical method. The monoaqua adduct [Pt（NH3）2Cl（H2O）]^＋ and the diaqua adduct [Pt（NHa）2（H2O）2]^2＋ were chosen for coordination to the N（7） site of guanine tautomers. The results demonstrate that the platinum adducts influence moderately on tautomeric equilibrium, but do not change the relative stability of tautomers whether in gas phase or in aqueous solution. The keto form having H atom at N（1） and N（9） was always the predominant structure when cisplatin adducts were bound to guanine. However, other forms could coexist in water. Meanwhile, our calculations suggest that the tautomeric equilibrium should be via the same intermediate.     

An ab initio SCF study on the tautomerisation of the 8-oxo-guanine and xanthine

All possible H₉-tautomers of 8-oxo-guanine and xanthine were studied by means of PM3 semiempirical and DFT (density functional theory) quantum chemistry methods. Additionally, the five most stable tautomers of both guanine derivatives were estimated on 3-21G, 6-31G, 6-31G^** and MP2 (6-31G^**) ab initio levels. The impact of the environment polarity on the tautomeric equilibrium was also taken into account. Among the variety of tautomeric isomers most probable are diketo forms of both studied derivatives in non-polar and polar surroundings.

The tautomeric equilibrium was unchanged after connection of the sugar backbone. The most preferred diketo forms of 8-oxo-guanosine and xanthidine are in syn conformations both in polar and non-polar environments. The increase of the syn conformations over anti ones may have the source in the formation of the internal hydrogen bonds between H′5 and N₃ atoms. The calculated values of the pseudorotation phase angle were between 144 and 180° in all cases. This corresponds to C′2-endo conformations of all optimised structures.

The N-glycosidic bond stability of most stable tautomers was compared to standard guanosine. Most tautomers of 8-oxo-guanosine and xanthidine are characterised by more stable C₁′-N₉ bond. This indicates that both these derivatives are hardly susceptible to spontaneous depurination and its removal from the DNA will depend mostly on the activity of DNA repair enzymes.     

Theoretical studies on interactions between low energy electrons and protein-DNA fragments: valence anions of AT-amino acids side chain complexes

Electron attachment to trimeric complexes that mimic most frequent hydrogen bonding interactions between an amino acid side chain (AASC) and the Watson-Crick (WC) 9-methyladenine-1-methylthymine (MAMT) base pair has been studied at the B3LYP/6-31++G(d,p) level of theory. Although the neutral trimers will not occur in the gas phase due to unfavorable free energy of stabilization (G(stab)) they should form a protein-DNA complex where entropy changes related to formation of such a complex will more than balance its disadvantageous G(stab). The most stable neutrals possess an identical pattern of hydrogen bonds (HBs). In addition, the proton-acceptor (N7) and proton-donor (N10) atoms of adenine involved in those HBs are located in the main groove of DNA. All neutral structures support the adiabatically stable valence anions in which the excess electron is localized on a π* orbital of thymine. The vertical detachment energies (VDEs) of anions corresponding to the most stable neutrals are substantially smaller than that of the isolated WC MAMT base pair. Hence, electron transfer from the anionic thymine to the phosphate group and as a consequence formation of a single strand break (SSB) should proceed more efficiently in a protein-dsDNA complex than in the naked dsDNA as far as electron attachment to thymine is concerned.     

Can anion interaction accelerate transformation of cytosine tautomers? Detailed view form QTAIM analysis

The relative stabilities and noncovalent interactions of the six low-lying energy tautomers of cytosine nucleobase with some biological anions (such as F^?, Cl^?, and CN^?) have been investigated in gas phase by density functional theory (DFT) method in conjunction with 6-311++G (d,p) atomic basis set. Furthermore, to systematically investigate all possible tautomerisms from cytosine induced by proton transfer, we describe a study of structural tautomer interconversion in the gas phase and in a continuum solvent using DFT calculation. We carried out geometrical optimizations with the integral equation formalism of polarizable continuum (IEF-PCM) model to account for the solvent effect, and the results were compared with those in the gas phase. The result of calculation revealed that anions bind mostly in a bidentate manner via hydrogen bond, and stabilization energies of these complexes are larger than those in the case when anions bind in a monodentate manner. The quantum theory of atoms in molecules (QTAIM), natural bonding orbital (NBO) and energy decomposition analysis (EDA) have also been applied to understand the nature of hydrogen bond interactions in these complexes. NBO analysis reveals that the interaction patterns between the anions and the tautomers are ??-type interaction between lone pairs and $ \sigma_{{_{{{\text{N}}--{\text{H}}}} }}^{*} $ , $ \sigma_{{_{{{\text{O}}--{\text{H}}}} }}^{*} $ and $ \sigma_{{_{{{\text{H}}--{\text{F}}}} }}^{*} $ antibonding orbitals. Also, according to these theories, the interactions are found to be partially electrostatic and partially covalent. EDA results identify that these bonds have less than 35% covalent character and more than 65% electrostatic, and the covalent character increases in different anions in the order F^??>?CN^??>?Cl^?. On the other hand, orbital interaction energies of complexes of F^? anion are more than those of Cl^? and CN^? complexes. The lower orbital interaction energies in complexes of Cl^? and CN^?anions imply less charge transfer and stronger ionic bond character. Furthermore, relationship between the orbital interaction energy (E ²) with hydrogen bonding energy (E _H···X) and the electron density (??(r)) with hydrogen bonding energy of F^?, Cl^? and CN^? complexes have also been investigated.     

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.     

BX4- and AlX4- superhalogen anions (X = F, Cl, Br): an ab initio study

The vertical electron detachment energies (VDEs) of 30 MX 4 (-) (M = B, Al; X = F, Cl, Br) anions were calculated at the OVGF level with the 6-311+G(3df) basis sets. The largest vertical electron binding energy was found for the AlF 4 (-) system (9.789 eV). The strong VDE dependence on the symmetry of the species, ligand type, ligand-central atom distance, and bonding/nonbonding/antibonding character of the highest occupied molecular orbital was observed and discussed.     

Photoelectron spectrum of valence anions of uracil and first-principles calculations of excess electron binding energies

The photoelectron spectrum (PES) of the uracil anion is reported and discussed from the perspective of quantum chemical calculations of the vertical detachment energies (VDEs) of the anions of various tautomers of uracil. The PES peak maximum is found at an electron binding energy of 2.4 eV, and the width of the main feature suggests that the parent anions are in a valence rather than a dipole-bound state. The canonical tautomer as well as four tautomers that result from proton transfer from an NH group to a C atom were investigated computationally. At the Hartree-Fock and second-order Moller-Plesset perturbation theory levels, the adiabatic electron affinity (AEA) and the VDE have been converged to the limit of a complete basis set to within +/-1 meV. Post-MP2 electron-correlation effects have been determined at the coupled-cluster level of theory including single, double, and noniterative triple excitations. The quantum chemical calculations suggest that the most stable valence anion of uracil is the anion of a tautomer that results from a proton transfer from N1H to C5. It is characterized by an AEA of 135 meV and a VDE of 1.38 eV. The peak maximum is as much as 1 eV larger, however, and the photoelectron intensity is only very weak at 1.38 eV. The PES does not lend support either to the valence anion of the canonical tautomer, which is the second most stable anion, and whose VDE is computed at about 0.60 eV. Agreement between the peak maximum and the computed VDE is only found for the third most stable tautomer, which shows an AEA of approximately -0.1 eV and a VDE of 2.58 eV. This tautomer results from a proton transfer from N3H to C5. The results illustrate that the characteristics of biomolecular anions are highly dependent on their tautomeric form. If indeed the third most stable anion is observed in the experiment, then it remains an open question why and how this species is formed under the given conditions.     

Application of natural bond orbital analysis to delocalization and aromaticity in C-substituted tetrazoles.

Energies of two tautomeric forms of 10 tetrazole derivatives substituted at C5 were established by DFT/B3LYP calculations carried out at the 6-311++G level. In each case the calculated energy of the 2H-tautomer was lower than that of the 1H. Furthermore, three geometric aromaticity indices of both forms were calculated, as were the values of nuclear independent nuclear shift and aromatic stabilization energy. The electronic properties were evaluated with the help of the natural bonding orbital theory. Following this a new pi-delocalization parameter, the root-mean square of pi-electron density localized on the atoms of the five-membered tetrazole ring, SDn, was introduced. It was concluded that the electronic delocalization can be described equally well by three different parameters: SDn, the extent of the transfer of electron density from the p(z) orbital of one nitrogen to the rest of the pi electron system, and population of two antibonding pi orbitals. Arguably, the information provided by the electronic parameters is similar to that contained in the geometric (structural) aromaticity indices except for tetrazole substituted by -BH(2).     

Quantum-Chemical Studies of the Structure and Reactivity of Pyrazol-5-ones and Their Thio and Seleno Analogs: V. Effect of Electron Correlation on Tautomerism and Acidity of 1-Methylheteropyrazolones

The relative stability of tautomeric forms of 1-methyl-substituted heteropyrazolones (O, S, Se) and their gas-phase acidity were estimated by DFT calculations with various basis sets and methods of geometry optimization. The electron correlation effects make an appreciable contribution to the Gibbs free energies of their tautomers and anions, especially those containing the heavy atoms. The qualitative pattern of tautomerism in pyrazolones is essentially similar to that obtained by semiempirical and nonempirical RHF calculations: The most stable is the CH form. For hetero analogs, consideration of electron correlations effects increases the relative stability of SH (SeH) forms. The series of relative acidity of the compounds depending on the heteroatom is preserved (Se O).     

Evaluation of the origin of conformational and tautomeric preferences in N‐formylformamide – A quantum chemical study

Quantum chemical study of N‐formylformamide (NFF) was carried out at various theoretical levels and the determinate equilibrium conformations were recomputed at the high level ab initio methods such as G2MP2, G2, G3, and complete basis set (CBS)‐QB3. The computational results reveal that the amide resonance and intramolecular hydrogen bonding are two superior factors in determining the most stable conformation of diamide (DA) and amide–imidic (AI) acid tautomers, respectively. The evaluation of hydrogen bond energies predicts that the hydrogen bond (HB( strength of NFF is weaker than the malonaldehyde (MA). But the results of atoms in molecules (AIM(, natural bond orbital (NBO), and geometrical parameters are given a different order, E_HB(NFF) > E_HB(MA). Although the bond average energies of tautomerization process emphasized on more stability of AI tautomer, but our theoretical calculations reveal that the DA conformers are more stable than the AI ones. The population analyses of equilibrium conformations by NBO method also predict that the origin of tautomeric preference is mainly because of the electron delocalization of amide functional group, especially LP(N)→ π*_C?O charge transfer. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Quantum chemical studies on prototautomerism of 1H‐imidazo[4,5‐c]pyridine

The structural features of the 1H‐imidazo[4,5‐c]pyridine (ICPY) tautomers and homodimers of the most stable tautomers have been studied by quantum chemical methods. FTIR and Raman spectra of the ICPY were recorded in the range of 4000–60 cm^?1 and 3500–5 cm^?1. The predominant tautomer among four possible isomers of ICPY were determined. The optimized geometries and vibrational frequencies of possible ICPY tautomers and dimers were computed by B3LYP/DFT method with 6‐311++G(d,p) and 6‐31G(d) basis sets. All vibrational frequencies assigned in detail with the help of total energy distribution (TED) and isotopic shifts. ICPY dimeric forms were also characterized according to their hydrogen bonding interactions, and it has been found that the most stable ICPY homodimer establishes moderate strong N ? H …N type hydrogen bond. ¹H NMR, ¹³C NMR, and ¹⁵N NMR properties have been calculated for all tautomeric forms using the gauge independent atomic orbital (GIAO) method. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Integration of electron density and molecular orbital techniques to reveal questionable bonds: the test case of the direct Fe-Fe bond in Fe2(CO)9

The article illustrates the advantages of partitioning the total electron density rho(rb), its Laplacian (inverted Delta)2 rho(rb), and the energy density H(rb) in terms of orbital components. By calculating the contributions of the mathematically constructed molecular orbitals to the measurable electron density, it is possible to quantify the bonding or antibonding character of each MO. This strategy is exploited to review the controversial existence of direct Fe-Fe bonding in the triply bridged Fe2(CO)9 system. Although the bond is predicted by electron counting rules, the interaction between the two pseudo-octahedral metal centers can be repulsive because of their fully occupied t(2g) sets. Moreover, previous atoms in molecules (AIM) studies failed to show a Fe-Fe bond critical point (bcp). The present electron density orbital partitioning (EDOP) analysis shows that one sigma bonding combination of the t(2g) levels is not totally overcome by the corresponding sigma* MO, which is partially delocalized over the bridging carbonyls. This suggests the existence of some, albeit weak, direct Fe-Fe bonding.     

Tautomeric forms of azolide anions: vertical electron detachment energies and Dyson orbitals

Several decouplings of the electron propagator, including the relatively new P3+ approximation for the self-energy, have been used to calculate vertical electron detachment energies of tautomeric forms of closed-shell, pentagonal, aromatic anions in which ring carbons without bonds to hydrogens appear. This study extends previous work in which the most stable forms of anionic, five-member rings with one to five nitrogens were considered. Whereas the lowest electron detachment energies sometimes are assigned by Koopmans's theorem results to pi orbital vacancies, electron propagator calculations always obtain sigma orbital vacancies for the ground states of the doublet radicals. Higher electron detachment energies that correspond to excited doublets with pi vacancies also are presented. The predicted transition energies are in good agreement with low-intensity peaks in recent anion photoelectron spectra that have been assigned to less stable, tautomeric forms of these anions.     

Spodium Bonds: Noncovalent Interactions Involving Group 12 Elements

The term spodium (Sp) bond is proposed to refer to a net attractive interaction between any element of Group 12 and electron‐rich atoms (Lewis bases or anions). These noncovalent interactions are markedly different from coordination bonds (antibonding Sp–ligand orbital involved). Evidence is provided for the existence of this interaction by calculations at the RI‐MP2/aug‐cc‐pVTZ level of theory, atoms‐in‐molecules, and natural bond orbital analyses and by examining solid‐state structures in the Cambridge Structure Database.     

Estimating local bonding/antibonding character of canonical molecular orbitals from their energy derivatives. The case of coordinating lone pair orbitals

According to Koopmans theorem, the derivative of the energy of a canonical molecular orbital (MO) with respect to nuclear coordinates quantifies its bonding/antibonding character. This quantity allows predictions of bond length variation on ionisation in a panel of 19 diatomic species. In polyatomic molecules, the derivative of a MO energy with respect to a given bond length reveals the nature and the degree of the bonding/antibonding contribution of this MO with respect to this bond. Accordingly, the HOMO “lone pairs” of CO and CN^? and the HOMO‐2 of CH₃CN are found to be antibonding with respect to the C? X bond (X = N, O), whereas the HOMO of N₂ is found to be bonding. With the same approach, the variation of the bonding character in the MOs of CO and CH₃CN on interaction with an electron acceptor (modeled through the approach of a proton) or by applying an electric field was studied. © 2016 Wiley Periodicals, Inc.