Skip to main content

ORIGINAL RESEARCH article

Front. Chem., 06 March 2015
Sec. Chemical Biology
This article is part of the Research Topic Radiation-induced and oxidative DNA damages View all 10 articles

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

\r\nJos P. Cern-Carrasco*José P. Cerón-Carrasco1*Alberto RequenaAlberto Requena1Jos ZúigaJosé Zúñiga1Denis Jacquemin,Denis Jacquemin2,3
  • 1Departamento de Química Física, Universidad de Murcia, Murcia, Spain
  • 2Chimie et Interdisciplinarité, Synthèse, Analyse, Modélisation, UMR Centre National de la Recherche Scientifique, Université de Nantes, Nantes, France
  • 3Institut Universitaire de France, Paris, France

We investigate the attack of the nitrogen dioxide radical (NO2) 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 NO2 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 NO2-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–NO2 radical adducts likely involve side reactions, e.g., the GC deprotonation to the solvent, rather than proton exchange between guanine and cytosine basis.

1. Introduction

Free radicals are naturally present in the biological medium as intermediates of the cellular metabolism. Although necessary to keep the normal biochemical activity, these highly reactive species tend to react, when present in excess, with a wide panel of biomolecules including DNA (Cadet et al., 2010). Accordingly, free-radical-induced DNA damage (mutations) might eventually initiate degenerative diseases, cardiovascular problems and cancers (von Sonntag, 2009). The reactive nitrogen oxide species (RNOS) family is a representative example of the double role played by free radicals (Lonkar and Dedon, 2011). Indeed, it has been established that mammalians cells endogenously produce RNOS as critical biological mediators (Ignarro, 1990). However, intracellular amounts of RNOS can overpass the tolerable limit by exposure to external chemical agents, such as cigarette smoke or air pollution (Pacher et al., 2007), and physical agents like high-energy radiations (Douki and Cadet, 2008).

One of the most important RNOS is the nitrogen dioxide radical, which is mainly present under its monomer NO2 form in physiological conditions (Augusto et al., 2002). Like any free radical, NO2 might react by electron transfer, hydrogen atom abstraction and/or radical addition to unsaturated bonds (Galano, 2007; Cerón-Carrasco et al., 2010). As for a reaction with DNA is concerned, the latter mechanism involves the addition of the radical to one of the unsaturated bonds of the guanine–cytosine (GC) base pair. There are, in particular, five possible reactive positions in the GC pair: C4, C5, and C8 in the guanine [C4(G), C5(G), and C8(G)] and C5 and C6 in the cytosine [C5(C) and C6(C)]. The atomic numbering is shown in Figure 1. Previous experimental and theoretical evidences indicated that C8(G) is the most reactive site for OH radical attack (Fortini et al., 2003; Shukla et al., 2004; Jena and Mishra, 2005; Shukla and Mishra, 2007; Zhang and Eriksson, 2007; Bergeron et al., 2010; Cerón-Carrasco and Jacquemin, 2012). Agnihotri and Mishra have also shown that NO2 quickly reacts with the guanine radical cation (G•+) at this same position (Agnihotri and Mishra, 2009, 2010, 2011), which might lead to DNA damage through the formation of an intermediate 8-nitroguanine structure (Misiaszek et al., 2005). The reaction of NO2 radical with the undamaged (neutral and closed-shell) GC pair has been, however, much less explored, so additional work is required to understand the biological action of such radical once inside the cellular medium.

FIGURE 1
www.frontiersin.org

Figure 1. Chemical structure and atomic numbering of the canonical GC base pair. Both N9 in G and N1 in C are the connections to the lateral backbone in DNA. In an isolated GC base pair, these positions are occupied by two hydrogen atoms. Arrows show the possible displacements of the protons H1, H4, and H4′ during the tautomeric reactions.

In this work, we use an state-of-the-art theoretical approach within the density functional theory (DFT) framework to investigate the possible mutagenic effects induced by the NO2 radical. In particular, we study the formation of non-canonical Watson and Crick structures arising from activation of the tautomeric equilibria in the GC base pair. These equilibria involve proton transfer (PT) between the guanine and cytosine basis yielding the so-called rare tautomeric forms that have been proposed as a plausible source of genetic errors (Löwdin, 1963; Florián and Leszczyński, 1996; Dannenberg and Tomasz, 2000; Gorb et al., 2004; Kumar and Sevilla, 2010; Villani, 2010a; Jacquemin et al., 2014). In undamaged DNA, GC base pairs naturally exchange the position of the H1, H4 and H4′ protons (see Figure 1), eventually disrupting the standard hydrogen-bond pattern that maintain the double helix bounded and thus inducing genetic errors during the cell replication (Villani, 2010b). Consequently, the tautomeric equilibria can be altered by external agents facilitating PT reactions (Khanduri et al., 2011). The main goal of this contribution is therefore to determine whether NO2 radical contributes to the mutagenic processes by shifting PT mechanism toward rare-tautomeric GC structures.

2. Methods

Since all the structures under study correspond to radical forms, the bullet superscript is hereafter omitted in the text in order to lighten the notation. We employ two different chemical models. In an early stage, a single GC base pair is used to determine the relative stabilities of all possible GC–radical adducts. More specifically, the NO2 radical is added to all possible reactive sites of the GC base pair through the two possible binding modes, e.g., ONO–CX and O2N–CX, where CX stands for the attacked carbons: C4(G), C5(G), C8(G), C5(C), and C6(C). The resulting adducts are next fully solvated with eleven water molecules, as illustrated in Figure 2, in order to mimic the solvent effects of a single solvated GC pair (Kumar et al., 2008). All these structures are optimized without symmetry constraints at the M06-2X/6-311G(d,p) level (Zhao and Truhlar, 2008a,b). The vibrational frequencies are computed at the same level of theory to confirm the absence of imaginary frequencies so that real minima are obtained. The total energies of all stationary points are then improved using the same exchange-correlation functional with the more extended atomic basis set, 6-311++G(d,p). Also, the well-known polarisable continuum model (PCM) is additionally used to account for solvent effects beyond the explicitly-treated first hydration shell (Tomasi et al., 1999, 2005).

FIGURE 2
www.frontiersin.org

Figure 2. Chemical structure for two of the single GC–radical adduct models: O2N–C8(G) and ONO–C6(C).

The isolated GC–radical adducts might form indeed over-distorted structures compared to real DNA, since the stacking effect is neglected. Consequently, the single base pair model is next refined by sandwiching the most stable adduct into at double-stranded B-form trimer d(5′-GGG-3′)d(3′-CCC-5′), as shown in Figure 3. This model includes both nucleotides moieties and sugar-phosphate backbone chains (Chen et al., 2009, 2011; Cerón-Carrasco et al., 2011). Since phosphates are negatively charged at physiological pH, the neutrality of the systems can be ensured by either adding counterions, e.g., Na+ or K+, or protonating the phosphate groups. Since the backbone does not significantly affect PT reactions in DNA (Close and Øhman, 2008; Chen et al., 2009), we go for the latter procedure to have a neutral system and consequently facilitate the convergence of the calculations. The GC-radical adduct is therefore located between two GC base pairs, which correctly mimic the π-stacking interactions in real DNA (Chen et al., 2009). We should underline that the stacked DNA model includes five water molecules in the vicinity of N3(G), O6(G), N7(G), O2(C), and N4(C) atoms (see atomic numbering in 1), and not eleven as in the single GC–radical adduct. These five sites correspond to the solvent exposed heteroatoms in DNA, so five water molecules are enough to provide a suitable representation of the hydrogen-bonds network formed between the base pair and the first hydration shell of the stacked structure (Schneider and Berman, 1995; Auffinger and Westhof, 2000; Makarov et al., 2002).

FIGURE 3
www.frontiersin.org

Figure 3. The ONIOM model for the DNA-embedded GC-radical adduct. The high layer for the ONO-C8(G) adduct is represented as ball and sticks whereas the confining base pair of the medium layer are shown as tube, and the lateral back-bone placed in the low layer are displayed as wireframe.

Due to the size of the DNA fragment (ca. 200 atoms), we use the hybrid ONIOM approach (Dapprich et al., 1999), which allows us to split the system into different regions, called layers. In this system, the layer of interest (defined as high layer) is the central GC–NO2 radical adduct plus the five explicit water molecules. This layer is treated at the same level that the isolated model, i.e., by using M06-2X/6-311G(d,p) for the optimization and PCM-M06-2X/6-31++G(d,p) for single-point energy calculations (Vreven et al., 2001; Mo et al., 2004). The two border GC base pairs that confine the central GC-radical adduct form the medium layer, which is described at the less demanding M06-2X/6-31G(d) level, while the lateral sugar-phosphate backbone (low layer) is simulated with the semiempirical PM6 method (Stewart, 2007). In this computation protocol, only the high layer is optimized whereas the rest of the system is frozen in the space. Such partial optimization procedure retains the characteristic double helix form of DNA while allowing relaxation of the central base pairs during the proton transfer reaction. Analytic calculations of the vibrational frequencies confirm the nature of the obtained structures to be minima (no imaginary frequencies) or transition states (a single imaginary frequency corresponding to the stretching of the transferred protons). All calculations are carried out with Gaussian09 (Frisch et al., 2009).

3. Results and Discussion

Although previous experimental and theoretical evidences suggest that the C8(G) is the preferred site for adding radicals to the GC base pair (Kumar et al., 2011), we have calculated the relative stabilities of all possible GC–radical adducts to provide an accurate picture of the NO2 radical reactivity. The relative energies obtained for the located adducts are listed in Table 1, with the reference value being taken as the most stable GC–radical adduct: ONO–C8(G). As observed, there are two main radical-attacking sites: C8(G) and C6(C). It wasn't possible to optimize GC-radical adduct at other positions, with the only exception of ONO–C5(C) at a significantly higher energy (11.76 kcal.mol−1). This finding is consistent with previous results obtained by Zhang and Eriksson using DFT calculations to determine the reactivity of the OH radical with the neutral closed-shell GC base pair (Zhang and Eriksson, 2007). These authors concluded that the two most stable binding sites for the GC–OH mutation are C8(G) and C6(C) (Zhang and Eriksson, 2007). In contrast, there are two possible binding modes (O2N– and ONO–) in the present case. Inspection of Table 1 reveals that the attack of the NO2 radical could proceed equally through the nitrogen or the oxygen atom. This is a remarkable result: the binding mode does not significantly affect the final stability of the generated adduct since the reported difference is negligible (<0.40 kcal.mol−1).

TABLE 1
www.frontiersin.org

Table 1. Relative energy (ΔE/kcal.mol−1), theoretical interbase bond distances (in Å) and Mulliken atomic charges (in |e|) for the single GC–radical adducts.

Aiming to characterize better the structural changes induced by the radical, we also include in Table 1 the interbase distances and the partial atomic charges of the protons involved in the intramolecular hydrogen bonds. Although there are slight differences (ca. 0.1 Å) in the interbase separation, adding the NO2 radical does significantly alter the GC geometry. The structures of the two most stable GC–radical adducts, namely O2N–C8(G) and ONO–C8(G), are shown in Figure 4. The side views show a clear bent structure as a result of the radial addition. Such departure from planarity can be quantified through the dihedral angle formed by C6(G)C2(G)C2(C)C4(C) (see atomic numbering in Figure 1), which has a value of around 1° for undamaged GC, but close to 10° in both GC–radical adducts. This is a logical consequence of the selected model, which mimics the reactivity of a single GC base pair in solution but does not reproduce the DNA environment, and in particular the π-stacking constrain.

FIGURE 4
www.frontiersin.org

Figure 4. Top and side views of the optimized structure for the O2N–C8(G) and ONO–C8(G) adducts.

More intriguing are the results shown in Table 1 for the atomic charges. As observed, the O2N– biding mode results in larger (more positive) atomic charges for the transferred proton during the tautomeric reaction, e.g., H4 and H1, compared to its ONO– counterpart. Although limited, these changes might affect the stability of the induced mutation since the atomic charges are directly related to the proton acidity and, consequently, to the tendency of forming tautomers.

The single GC base pair model provides the first clues for the attack of the NO2 radical to DNA, but the more refined model is required for biological analysis because (i) bend structures imply energetic penalties once inside DNA due to the π-stacking interactions; and (ii) the location of the radical could be modified by the chemical environment of DNA, e.g., by the interaction with other base pairs and/or later sugar-phosphate backbone. Consequently, the two most stable adducts, O2N–C8(G) and ONO–C8(G), are “docked” into the DNA fragment shown in Figure 3. In contrast with the single GC model in which both binding modes are practically isoenergetic, the ONO–C8(G) adduct is now 3.23 kcal.mol−1 more stable than the O2N–C8(G) adduct in DNA, which confirms the importance of accounting for all key interactions.

In addition to the optimal geometry of the canonical GC–radical adducts, we have computed the energetic profiles corresponding both to the transfer of protons H1 and H4, directly exchanged between the basis, and to the tautomeric equilibria arising from the movement of the H4′ (see Figure 1). The latter equilibrium implies a water-assisted mechanism in which the surrounding water molecules catalyze the process by accepting and donating protons (see Ref. Jacquemin et al., 2014 for details). The optimized rare tautomeric structures are shown in Figure 5, and their relative energies and equilibrium constant are given in Table 2. We should note that all the attempts tried to optimize products arising from a single-PT reaction, that is, to optimize rare tautomers where only one proton (H1, H4 or H4′) is shifted from its position, failed as the geometry optimization quickly restored the original canonical structure. However, we were able to localize the products for H1+H4 and H1+H4′ double-PT with both O2N–C8(G) and ONO–C8(G) adducts, which structures are shown in Figure 5. As discussed above, the O2N–C8(G) binding mode results in more acidic H4 and H1 protons, compared to the ONO–C8(G) adduct, eventually favouring the tautomeric equilibria by ca. 3 kcal.mol−1. The associated equilibrium constants (Keq) can be estimated from the relative energies obtained using the expression

Keq=eΔE/RT    (1)

where R is the ideal gas constant and T is the temperature (298.15 K). Keq provides a measure of the mutation frequency, which can be directly compared to the computed value for the radical-free DNA. Taking the undamaged-DNA as a reference for the impact of H1+H4 and H1+H4′ reactions, one notices that the latter is unlikely to occur due to the low Keq values, which are of the order of 10−16 and 10−13 for the O2N–C8(G) and ONO–C8(G) adducts, respectively. However, the energies listed in Table 2 hint that the O2N–C8(G) adduct stabilizes the H1+H4 mutation by ca. 2 kcal.mol−1, slightly shifting the equilibrium toward the rare tautomeric form (Keq = 7.27 × 10−08). Despite of such a difference, all rare tautomers are in the range of observed frequency for the spontaneous mutation, experimentally measured as ranging from 10−8 to 10−10 (Topal and Fresco, 1976). Interestingly, there is a certain dissimilarity between NO2 and OH. The latter clearly shifts the PT reactions toward the rare tautomeric forms of the GC base pair, and its degenerative effects are therefore governed by the transition states along the PT reactions, which in turn determines their lifetimes (Zhang and Eriksson, 2007; Cerón-Carrasco and Jacquemin, 2012). In contrast, NO2 does not affect the PT reactions, so it is not necessary to compute the transition states to conclude that they make no impact in the global mutation when present in DNA during cell replication.

FIGURE 5
www.frontiersin.org

Figure 5. Chemical structures of the canonical rare tautomeric form optimized for the O2N–C8(G) and ONO–C8(G) adducts: H1+H4 and H1+H4′. For the sake of clarity only the central GC–radical adduct and the surrounding water molecules are displayed as balls and sticks, while the rest of the system is shown has wireframe.

TABLE 2
www.frontiersin.org

Table 2. Relative energy (ΔE/kcal.mol−1) and equilibrium constants (Keq) at 298 K calculated for the H1+H4 and H1+H4′ tautomeric equilibria of undamaged-DNA, O2N–C8(G) and ONO–C8(G) adducts.

Since the canonical O2N-C8(G) and ONO-C8(G) adducts are shown to be the most stable structures at an early stage of the radical attack, we decided to compute the spin densities surfaces for both forms to determine the localisation of the unpaired electron. Such calculations were performed by single-point full QM calculations at the M06-2X/6-311G(d,p) level. As illustrated in Figure 6, the spin density is distributed in the guanine as well as in the radical, in agreement with previous results obtained for the OH radical (Cerón-Carrasco and Jacquemin, 2012). This finding indicates that the subsequent mutagenic steps will be take place at this region rather than at further sugar-phosphate backbones moieties or at other DNA basis. It becomes therefore necessary to perform additional dynamics calculations with larger models (Loos et al., 2009; Ambrosek et al., 2010; Cauët et al., 2010; Dupont et al., 2011; Garrec et al., 2012; Cerón-Carrasco et al., 2013; Dupont et al., 2013; Patel et al., 2013) for the O2N–C8(G) and ONO–C8(G) adducts to further investigate the damage mechanism initiated by GC-NO2 structures. Our results indicate that such adducts should evolve by side mechanisms different than PT between guanine an cytosine, e.g., deprotonation to the solvent (Agnihotri and Mishra, 2011).

FIGURE 6
www.frontiersin.org

Figure 6. Spin density distributions of the canonical radical adducts (isovalue = 0.002 a.u.) computed at M06-2X/6-311G(d,p) level for the three base pairs DNA fragment.

4. Concluding Remarks

We have performed theoretical calculations to determine the reaction mechanism between the NO2 radical and the neutral closed-shell guanine–cytosine base pair. The reported data demonstrate that the carbon C8 of the guanine moiety is the preferential site to form the GC-radical adduct. The NO2 radical equally reacts with the single GC base pair through either its nitrogen or its oxygen atom, with the latter yielding the most stable adduct when a more realistic DNA-embedded model is used. This finding confirms the importance of using a well-balanced chemical model to obtain valuable conclusions. For the radical addition to DNA, the π-stacking interaction not only constrain the planarity of the attacked base pair, but also accommodates better the ONO– binding mode compared to its O2N– counterpart. We have assessed the possible evolution of these adducts through proton transfer reactions between guanine and cytosine basis, which is one of the sources of genetic errors in DNA. The associated equilibrium constants lie in the same range as the one observed for spontaneous mutation. We conclude that the initial ONO–C8(G) does not promote rare tautomeric forms by proton exchange, at least not at a rate larger than that in undamaged DNA. The NO2-induced genetic damage is expectedly to be initiated by other side reactions.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

JC acknowledges the support from the FP7 EU Marie Curie Actions through the Campus Mare Nostrum 37/38 CMN UMU Incoming Mobility Programme ACTion (U-IMPACT), and the support provided by Université de Nantes through its invited professor programme. DJ acknowledges the European Research Council (ERC) and the Région des Pays de la Loire for financial support in the framework of a Starting Grant (Marches - 278845) and a recrutement sur poste stratégique, respectively. This research work was partially supported by the Spanish Ministerio de Ciencia e Innovación under Projects CTQ2011-25872 and CONSOLIDER CSD2009-00038, and used resources of the CEN1 in Spain and the CINES and the CCIPL in France.

References

Agnihotri, N., and Mishra, P. C. (2009). Mutagenic product formation due to reaction of guanine radical cation with nitrogen dioxide. J. Phys. Chem. B 113, 3129–3138. doi: 10.1021/jp805942y

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Agnihotri, N., and Mishra, P. C. (2010). Formation of 8-nitroguanine due to reaction between guanyl radicalnitrogen dioxide: catalytic role of hydration. J. Phys. Chem. B 114, 7391–7404. doi: 10.1021/jp9122437

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Agnihotri, N., and Mishra, P. C. (2011). Reactivities of radicals of adenineguanine towards reactive oxygen speciesreactive nitrogen oxide species: OHNO2. Chem. Phys. Lett. 503, 305–309. doi: 10.1016/j.cplett.2011.01.042

CrossRef Full Text | Google Scholar

Ambrosek, D., Loos, P. F., Assfeld, X., and Daniel, C. (2010). A theoretical study of Ru(II) polypyridyl DNA intercalators structureelectronic absorption spectroscopy of [Ru(phen)2(dppz)]2+[Ru(tap)2(dppz)]2+ complexes intercalated in guanine−cytosine base pairs. J. Inorg. Biochem. 104, 893–901. doi: 10.1016/j.jinorgbio.2010.04.002

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Auffinger, P., and Westhof, E. (2000). Waterion binding around RNADNA (c,g) oligomers. J. Mol. Biol. 300, 1113–1131. doi: 10.1006/jmbi.2000.3894

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Augusto, O., Bonini, M. G., Amanso, A. M., Linares, E., Santos, C. C., and De Menezes, S. L. (2002). Nitrogen dioxidecarbonate radical anion: Two emerging radicals in biology. Free Radic. Biol. Med. 32, 841–859. doi: 10.1016/S0891-5849(02)00786-4

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Bergeron, F., Auvré, F., Radicella, J. P., and Ravanat, J. L. (2010). HO radicals induce an unexpected high proportion of tandem base lesions refractory to repair by DNA glycosylases. Proc. Natl. Acad. Sci. U.S.A. 107, 5528–5533. doi: 10.1073/pnas.1000193107

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Cadet, J., Douki, T., and Ravanat, J. L. (2010). Oxidatively generated base damage to cellular DNA. Free Radic. Biol. Med. 49, 9–21. doi: 10.1016/j.freeradbiomed.2010.03.025

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Cauët, E., Valiev, M., and Weare, J. H. (2010). Vertical ionization potentials of nucleobases in a fully solvated DNA environment. J. Phys. Chem. B 114, 5886–5894. doi: 10.1021/jp9120723

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Cerón-Carrasco, J. P., and Jacquemin, D. (2012). Interplay between hydroxyl radical attackh–bond stability in guanine–cytosine. RSC Adv. 2, 11867–11875. doi: 10.1039/c2ra22389a

CrossRef Full Text | Google Scholar

Cerón-Carrasco, J. P., Bastida, A., Zúñiga, J., Requena, A., and Miguel, B. (2010). A theoretical study of the reaction of β-carotene with the nitrogen dioxide radical in solution. J. Phys. Chem. B 114, 4366–4372. doi: 10.1021/jp911846h

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Cerón-Carrasco, J. P., Zúñiga, J., Requena, A., Perpète, E. A., Michaux, C., and Jacquemin, D. (2011). Combined effect of stackingsolvation on the spontaneous mutation in DNA. Phys. Chem. Chem. Phys. 13, 14584–14589. doi: 10.1039/c1cp20946a

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Cerón-Carrasco, J. P., Jacquemin, D., and Dumont, E. (2013). Impact of DNA environment on the intrastrand cross-link lesions: Hydrogen atom release as the last step of formation of G[8-5m]T. J. Phys. Chem. B 117, 16397–16404. doi: 10.1021/jp408947u

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Chen, H. Y., Kao, C. L., and Hsu, S. C. N. (2009). Proton transfer in guanine-cytosine radical anion embedded in B-form DNA. J. Am. Chem. Soc. 131, 15930–15938. doi: 10.1021/ja906899p

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Chen, H. Y., Yeh, S. W., Hsu, S. C. N., Kao, C. L., and Dong, T. Y. (2011). Effect of nucleobase sequence on the proton transfer reactionstability of the guanine–cytosine base pair radical anion. Phys. Chem. Chem. Phys. 13, 2674–2681. doi: 10.1039/c0cp01419b

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Close, D. M., and Øhman, K. T. (2008). Ionization energies of the nucleotides. J. Phys. Chem. A 112, 11207–11212. doi: 10.1021/jp805308p

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Dannenberg, J. J., and Tomasz, M. (2000). Hydrogen-bond acid/base catalysis: a density functional theory study of protonated guanine-(substituted) cytosine base pairs as models for nucleophilic attack on mitomycin in DNA. J. Am. Chem. Soc. 122, 2062–2068. doi: 10.1021/ja993145i

CrossRef Full Text | Google Scholar

Dapprich, S, Komáromi, I., Byun, K. S., Morokuma, K., and Frisch, M. J. (1999). A new ONIOM implementation in Gaussian 98. part 1. The calculation of energies, gradientsvibrational frequencieselectric field derivatives. J. Mol. Struct. 462, 1–21. doi: 10.1016/S0166-1280(98)00475-8

CrossRef Full Text | Google Scholar

Douki, T., and Cadet, J. (2008). Radiation Chemistry: From Basics to Applications in Materiallife Sciences. Les Ulis: edp Sciences.

Dupont, C., Patel, C., and Dumont, E. (2011). Improved DFT description of intrastrand cross-link formation by inclusion of London dispersion corrections. J. Phys. Chem. B 115, 15138–15144. doi: 10.1021/jp209074q

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Dupont, C., Patel, C., Ravanat, J. L., and Dumont, E. (2013). Addressing the competitive formation of tandem DNA lesions by a nucleobase peroxyl radical: a DFT-D screening. Org. Biomol. Chem. 11, 3038–3045. doi: 10.1039/c3ob40280k

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Florián, J., and Leszczyński, J. (1996). Spontaneous DNA mutations induced by proton transfer in the guanine·cytosine base pairs: An energetic perspective. J. Am. Chem. Soc. 118, 3010–3017. doi: 10.1021/ja951983g

CrossRef Full Text | Google Scholar

Fortini, P., Pascucci, B., Parlanti, E., D'Errico, M., Simonelli, V., and Dogliotti, E. (2003). 8-oxoguanine DNA damage: At the crossroad of alternative repair pathways. Mutat. Res. 531, 127–139. doi: 10.1016/j.mrfmmm.2003.07.004

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Frisch, M., J., Trucks, G. W., Schlegel, H. B., Scuseria, G. E., Robb, M. A., Cheeseman, J. R., et al. (2009). Gaussian 09 Revision D.01. Wallingford CT: Gaussian Inc.

Galano, A. (2007). Relative antioxidant efficiency of a large series of carotenoids in terms of one electron transfer reactions. J. Phys. Chem. B 111, 12898–12908. doi: 10.1021/jp074358u

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Garrec, J., Patel, C., Rothlisberger, U., and Dumont, E. (2012). Insights into intrastrand cross-link lesions of DNA from QM/MM molecular dynamics simulations. J. Am. Chem. Soc. 134, 2111–2119. doi: 10.1021/ja2084042

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Gorb, L., Podolyan, Y., Dziekonski, P., Sokalski. W. A., and Leszczyński, J. (2004). Double-proton transfer adenine-thymineguanine-cytosine base pairs. A post-hartree-fock ab initio study. J. Am. Chem. Soc. 126, 10119–10129. doi: 10.1021/ja049155n

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Ignarro, L. J. (1990). Biosynthesismetabolism of endothelium-derived nitric oxide. Annu. Rev. Pharmacol. Toxicol. 30, 535–560. doi: 10.1146/annurev.pa.30.040190.002535

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Jacquemin, D., Zúñiga, J., Requena, A., and Cerón-Carrasco, J. P. (2014). Assessing the importance of proton transfer reactions in DNA. Acc. Chem. Res. 47, 2467–2474. doi: 10.1021/ar500148c

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Jena, N. R., and Mishra, P. C. (2005). Mechanisms of formation of 8-oxoguanine due to reactions of onetwo OH radicalsthe H2O2 molecule with guanine: a quantum computational study. J. Phys. Chem. B 109, 1405–14218. doi: 10.1021/jp050646j

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Khanduri, D., Adhikary, A., and Sevilla, M. D. (2011). Highly oxidizing excited states of one-electron-oxidized guanine in DNA: wavelengthpH dependence. J. Am. Chem. Soc. 133, 4527–4537. doi: 10.1021/ja110499a

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Kumar, A., and Sevilla, M. D. (2010). Proton-coupled electron transfer in DNA on formation of radiation-produced ion radicals. Chem. Rev. 110, 7002–7023. doi: 10.1021/cr100023g

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Kumar, A., and Sevilla, M., and Suhai, S. (2008). Microhydration of the guanine-cytosine (GC) base pair in the neutralanionic radical states: a density functional study. J. Phys. Chem. B 112, 5189–5198. doi: 10.1021/jp710957p

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Kumar, A., Pottiboyina, V., and Sevilla, M. D. (2011). Hydroxyl radical (OH) reaction with guanine in an aqueous environment: a DFT study. J. Phys. Chem. B 115, 15129–15137. doi: 10.1021/jp208841q

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Löwdin, P. O. (1963). Proton tunneling in DNAits biological implications. Rev. Mod. Phys. 35, 724–732. doi: 10.1103/RevModPhys.35.724

CrossRef Full Text | Google Scholar

Lonkar, P., and Dedon, P. C. (2011). Reactive speciesDNA damage in chronic inflammation: reconciling chemical mechanismsbiological fates. Int. J. Cancer 128, 1999–2009. doi: 10.1002/ijc.25815

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Loos, P. F., Dumont, E., Laurent, A. D., and Assfeld, X. (2009). Important effects of neighbouring nucleotides on electron induced DNA single-strand breaks. Chem. Phys. Lett. 475, 120–123. doi: 10.1016/j.cplett.2009.05.041

CrossRef Full Text | Google Scholar

Makarov, V., Pettitt, B., and Feig, M. (2002). Solvationhydration of proteinsnucleic acids: a theoretical view of simulationexperiment. Acc. Chem. Res. 35, 376–384. doi: 10.1021/ar0100273

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Misiaszek, R., Crean, C., Geacintov, N. E., and Shafirovich, V. (2005). Combination of nitrogen dioxide radicals with 8-oxo-7,8-dihydroguanineguanine radicals in DNA: a oxidation and nitration end-products. J. Am. Chem. Soc. 127, 2191–2200. doi: 10.1021/ja044390r

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Mo, S. J., Vreven, M., Mennucci, B., Morokuma, K., and Tomasi, J. (2004). Theoretical study of the sn2 reaction of Cl (H2O)+CH3Cl using our own n-layered integrated molecular orbitalmolecular mechanics polarizable continuum model method (ONIOM–PCM). Theor. Chem. Acc. 111, 154–161. doi: 10.1007/s00214-003-0519-2

CrossRef Full Text | Google Scholar

Pacher, P., Beckman, J. S., and Liaudet, L. (2007). Nitric oxideperoxynitrite in healthdisease. Physiol. Rev. 87, 315–424. doi: 10.1152/physrev.00029.2006

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Patel, C., Garrec, J., Dupont, C., and Dumont, E. (2013). What singles out the G[8-5]C intrastrand DNA cross-link? Mechanisticstructural insights from quantum mechanics/molecular mechanics simulations. Biochemistry 52, 425–431. doi: 10.1021/bi301198h

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Schneider, B., and Berman, H. M. (1995). Hydration of the DNAA bases is local. Biophys. J. 69, 2661–2669. doi: 10.1016/S0006-3495(95)80136-0

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Shukla, L. I., Adhikary, A., Pazdro, R., Becker, D., and Sevilla, M. D. (2004). Formation of 8-oxo-7,8-dihydroguanine-radicals in γ-irradiated DNA by multiple one-electron oxidations. Nucleic Acids Res. 32, 6565–6574. doi: 10.1093/nar/gkh989

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Shukla, P. K., and Mishra, P. C. (2007). H2O3 as a reactive oxygen species: formation of 8-oxoguanine from its reaction with guanine. J. Phys. Chem. B 111, 4603–4615. doi: 10.1021/jp070399e

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Stewart, J. J. P. (2007). Optimization of parameters for semiempirical methods v: modification of NDDO approximationsapplication to 70 elements. J. Mol. Model. 13, 1173–1213. doi: 10.1007/s00894-007-0233-4

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Tomasi, J., Mennucci, M., and Cances, E. (1999). The IEF version of the PCM solvation method: An overview of a new method addressed to study molecular solutes at the QM ab initio level. J. Mol. Struct. 464, 211–226. doi: 10.1016/S0166-1280(98)00553-3

CrossRef Full Text | Google Scholar

Tomasi, J., Mennucci, B., and Cammi, R. (2005). Quantum mechanical continuum solvation models. Chem. Rev. 105, 2999–3093. doi: 10.1021/cr9904009

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Topal, M. D., and Fresco, J. R. (1976). Complementary base pairingthe origin of substitution mutations. Nature 263, 285–289. doi: 10.1038/263285a0

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Villani, G. (2010a). Theoretical investigation of hydrogen atom transfer in the adenine−thymine base pairits coupling with the electronic rearrangement. Concerted vs. stepwise mechanism. Phys. Chem. Chem. Phys. 12, 2664–2669. doi: 10.1039/b917672a

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Villani G. (2010b). Theoretical investigation of hydrogen atom transfer in the cytosine-guanine base pairits coupling with electronic rearrangement. concerted vs. stepwise mechanism. J. Phys. Chem. B 114, 9653–9662. doi: 10.1021/jp102457s

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

von Sonntag, C. (2009). Free-Radical-Induced DNA DamageIts Repair. Berlin; Heidelberg: Springer-Verlag.

Vreven, T., Mennucci, B., da Silva, C. O., Morokuma, K., and Tomasi, J. (2001). The ONIOM-PCM method: combining the hybrid molecular orbital methodthe polarizable continuum model for solvation. Application to the geometryproperties of a merocyanine in solution. J. Chem. Phys. 115, 62–72. doi: 10.1063/1.1376127

CrossRef Full Text | Google Scholar

Zhang, R. B., and Eriksson, L. A. (2007). Effects of oh radical addition on proton transfer in the guanine–cytosine base pair. J. Phys. Chem. B 111, 6571–6576. doi: 10.1021/jp071772l

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Zhao, Y., and Truhlar, D. G. (2008a). Density functionals with broad applicability in chemistry. Acc. Chem. Res. 41, 157–167. doi: 10.1021/ar700111a

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Zhao, Y., and Truhlar, D. G. (2008b). The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, transition elements: Two new functionalssystematic testing of four mM6-class functionals12 other functionals. Theor. Chem. Acc. 120, 215–241. doi: 10.1007/s00214-007-0310-x

CrossRef Full Text | Google Scholar

Keywords: guanine–cytosine, NO2 radical, induced mutation, proton transfer reaction, tautomeric equilibria, rare tautomers, computational chemistry, density functional theory

Citation: Cerón-Carrasco JP, Requena A, Zúñiga J and Jacquemin D (2015) Mutagenic effects induced by the attack of NO2 radical to the guanine-cytosine base pair. Front. Chem. 3:13. doi: 10.3389/fchem.2015.00013

Received: 08 January 2015; Accepted: 16 February 2015;
Published: 06 March 2015.

Edited by:

Antonio Monari, Université de Lorraine, France

Reviewed by:

Nino Russo, Università della Calabria, Italy
Leif Axel Eriksson, University of Gothenburg, Sweden

Copyright © 2015 Cerón-Carrasco, Requena, Zúñiga and Jacquemin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: José P. Cerón-Carrasco, Departamento de Química Física, Universidad de Murcia, 30100 Murcia, Spain jpceron@um.es

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.