DFT Study of the Hydroxyl Radical Addition to 2'-Deoxyguanosine and the Guanine Base in Four Double-Stranded B-Form Dimers

Density functional theory (DFT) calculations of reactions between 2'-deoxyguanosine (dR-Gua) and hydroxyl radical (HO^•) with water molecules (H₂O)_n, n = 0, 1, and 2, were carried out. The HO^• addition to three carbon sites, C(4), C(5), and C(8), and the subsequent ring cleavage of the three HO adducts were investigated. The addition to C(5) is of the smallest activation energy according to the largest lobe of the dR-Gua highest occupied molecular orbital (HOMO) at C(5). However, its adduct has small stability, and the C(8) adduct has the largest one. The C(8) adduct and the ring-opened amide have similar stability, which would lead to the apparent small yield of the former. Calculations were also performed on HO^• additions to the C(4) and C(8) sites of the guanine moiety of four dimer sequence models of B-form DNA with nucleotide moieties (a) 5'-GA-3', (b) 5'-GG-3', (c) 5'-GT-3', and (d) 5'-GC-3'. For instance, the (a) 5'-GA-3' model has a molecular formula C₃₉H₅₀N₁₅Na₂O₂₁P₂. The HO^• attack to C(4) is ruled out owing to the reinforced deformation of the parallel stacking of base pairs. The clear selectivity that the (b) 5'-GG-3' sequence is most reactive was found with the inclusion of the water dimer.

Exposing the G-quadruplex to electric fields: the role played by telomeres in the propagation of DNA errors

We use quantum calculations to assess the impact of external electric fields on the stability of G-quadruplex, a key structure in telomere functionality.

Mutagenic effects induced by the attack of NO2 radical to the guanine-cytosine base pair

We investigate the attack of the nitrogen dioxide radical (NO^•₂) to the guanine—cytosine (GC) base pair and the subsequent tautomeric reactions able to induce mutations, by means of density functional theory (DFT) calculations. The conducted simulations allow us to identify the most reactive sites of the GC base pair. Indeed, the computed relative energies demonstrate that the addition of the NO^•₂ radical to the C8 position of the guanine base forms to the most stable adduct. Although the initial adducts might evolve to non-canonical structures via inter-base hydrogen bonds rearrangements, the probability for the proton exchange to occur lies in the same range as that observed for undamaged DNA. As a result, tautomeric errors in NO₂-attacked DNA arises at the same rate as in canonical DNA, with no macroscopic impact on the overall stability of DNA. The potential mutagenic effects of the GC–NO^•₂ radical adducts likely involve side reactions, e.g., the GC deprotonation to the solvent, rather than proton exchange between guanine and cytosine basis.

Assessing the importance of proton transfer reactions in DNA

Although engineered by millions of years of evolution, the cellular machinery is not flawless, and errors regularly appear during DNA replication. The subsequent alteration of the stored genetic message results in a mutation and might be the starting point of important health disorders. The question therefore is what causes DNA mutations? All living organisms are constantly exposed to a number of external agents such as free radicals and to radiation, which may lead to induced mutations. There are also mutations happening without invoking the action of any exogenous element, the so-called spontaneous mutations. The former can be partially controlled by avoiding exposure to high-risk environments, while the latter are more intriguing because their origin is unclear and difficult to determine. As noted by Watson and Crick when they first discovered the DNA structure, the correct replication of DNA rests on the assumption that the base pairs remain in their most stable, canonical form. However, protons along the interbase hydrogen-bond network are not static entities. They can in fact interchange their positions in DNA bases through proton transfer (PT) reactions before strands unwind, giving rise to noncanonical structures defined as rare tautomers. The importance of these rare tautomers was also cleverly anticipated by Watson and Crick and some years later claimed by Löwdin to be a source of spontaneous mutations. In Watson and Crick's words: "It would be of interest to know the precise difference in free energy between the various tautomeric forms under physiological conditions." Unfortunately, rare tautomeric forms are very difficult to detect, so no direct and accurate free energy measure has been discerned. In contrast, theoretical chemistry is making good progress toward the quantification of PT reactions in DNA and their biological consequences. This Account touches upon the theoretical studies devoted to appraising the importance of rare tautomers as promoters of spontaneous mutations. We focus in particular on the crucial role played by the biological environment on DNA stability. It has now been demonstrated that valuable macroscopic predictions require not only highly accurate theories but also refined chemical models. Hybrid quantum mechanics/molecular mechanics (QM/MM) simulations performed on short but complete DNA sequence fragments emerge in this context as the most adequate tools. In addition, these methods can be used to quantify the effect of different external agents on the PT tautomeric equilibria and, eventually, to conveniently handle them. This is the case for the possible alteration of the naturally observed mutation rate by exposure to intense electric fields. Theoretical predictions envision in this respect promising applications of ultrashort electric pulses in medicine to selectively modify the mutated/canonical ratio in DNA.

Attacking mechanism of hydroxyl radical to DNA base-pair: density functional study in vacuum and in water

Recently, the influence of radiation on human body has been recognized as a serious problem. In particular, highly reactive hydroxyl radicals *OH produced by the radiation react with DNA, resulting in a great damage on its structure and electronic properties. It is thus important to investigate the reaction mechanism of *OH to DNA for elucidating the initial damage in DNA induced by the radiation. In the present study, we search for transition states (TS) of the reaction between G-C/A-T base-pair and [Formula: see text] in vacuum and in water, by the density functional theory (DFT) calculations. At first, we obtain the stable structures for the dehydrogenated G-C and A-T, in which the hydrogen atom of NH2 group of G or A base is abstracted by [Formula: see text]. From the structures of the dehydrogenated as well as the natural base-pairs, the TS between these structures is searched for and the activation free energy (AFE) is estimated for the reaction. In vacuum, AFEs for the G-C and A-T are almost the same each other, while the stabilization energy by the reaction for G-C is about 4.9 kcal/mol larger than that for A-T, indicating that the population of the dehydrogenated G-C is remarkably larger than that of the dehydrogenated A-T in vacuum. On the other hand, in water approximated by the continuum solvation model, the AFE for A-T is 2.6 kcal/mol smaller than that for G-C, indicating that the reaction dehydrogenated by [Formula: see text] occurs more frequently for the solvated A-T base-pair than G-C.

Sequence and conformation effects on ionization potential and charge distribution of homo-nucleobase stacks using M06-2X hybrid density functional theory calculations

DNA is subject to oxidative damage due to radiation or by-products of cellular metabolism, thereby creating electron holes that migrate along the DNA stacks. A systematic computational analysis of the dependence of the electronic properties of nucleobase stacks on sequence and conformation was performed here, on the basis of single- and double-stranded homo-nucleobase stacks of 1-10 bases or 1-8 base pairs in standard A-, B-, and Z-conformation. First, several levels of theory were tested for calculating the vertical ionization potentials of individual nucleobases; the M06-2X/6-31G* hybrid density functional theory method was selected by comparison with experimental data. Next, the vertical ionization potential, and the Mulliken charge and spin density distributions were calculated and considered on all nucleobase stacks. We found that (1) the ionization potential decreases with the number of bases, the lowest being reached by Gua≡Cyt tracts; (2) the association of two single strands into a double-stranded tract lowers the ionization potential significantly (3) differences in ionization potential due to sequence variation are roughly three times larger than those due to conformational modifications. The charge and spin density distributions were found (1) to be located toward the 5'-end for single-stranded Gua-stacks and toward the 3'-end for Cyt-stacks and basically delocalized over all bases for Ade- and Thy-stacks; (2) the association into double-stranded tracts empties the Cyt- and Thy-strands of most of the charge and all the spin density and concentrates them on the Gua- and Ade-strands. The possible biological implications of these results for transcription are discussed.