Original Research ARTICLE
Mutagenic effects induced by the attack of NO2 radical to the guanine-cytosine base pair
- 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 (NO•2) 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•2 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–NO•2 radical adducts likely involve side reactions, e.g., the GC deprotonation to the solvent, rather than proton exchange between guanine and cytosine basis.
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 NO•2 form in physiological conditions (Augusto et al., 2002). Like any free radical, NO•2 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 NO•2 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 NO•2 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. 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 NO•2 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 NO•2 radical contributes to the mutagenic processes by shifting PT mechanism toward rare-tautomeric GC structures.
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. 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. 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. 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.
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
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. 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. 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. 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.
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.
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
Agnihotri, N., and Mishra, P. C. (2011). Reactivities of radicals of adenineguanine towards reactive oxygen speciesreactive nitrogen oxide species: OH•NO•2. Chem. Phys. Lett. 503, 305–309. doi: 10.1016/j.cplett.2011.01.042
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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 firstname.lastname@example.org