# Electron holes in G-quadruplexes: The role of adenine ending groups

Evangelos Balanikas, Lara Martinez-Fernandez, Gérard Baldacchino, Dimitra Markovitsi

► To cite this version:

Evangelos Balanikas, Lara Martinez-Fernandez, Gérard Baldacchino, Dimitra Markovitsi. Electron holes in G-quadruplexes: The role of adenine ending groups. International Journal of Molecular Sciences, 2021, 22 (24), pp.13436. 10.3390/ijms222413436. cea-04115159

HAL Id: cea-04115159

https://cea.hal.science/cea-04115159

Submitted on 2 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Electron Holes in G-Quadruplexes: The Role of Adenine Ending Groups

Evangelos Balanikas<fnref n="1" />, Lara Martinez-Fernandez <fnref n="2" />,*, Gérard Baldacchino <fnref n="1" /> ID and Dimitra Markovitsi <fnref n="1" />,*,† ID

**Abstract:** The study deals with four-stranded DNA structures (G-Quadruplexes), known to undergo ionization upon direct absorption of low-energy UV photons. Combining quantum chemistry calculations and time-resolved absorption spectroscopy with 266 nm excitation, it focuses on the electron holes generated in tetramolecular systems with adenine groups at the ends. Our computations show that the electron hole is placed in a single guanine site, whose location depends on the position of the adenines at the 3' or 5' ends. This position also affects significantly the electronic absorption spectrum of (G⁺)• radical cations. Their decay is highly anisotropic, composed of a fast process (<2 µs), followed by a slower one occurring in ~20 µs. On the one hand, they undergo deprotonation to (G-H2)• radicals and, on the other, they give rise to a reaction product absorbing in the 300–500 nm spectral domain.

**Keywords:** G-Quadruplexes; electron holes; photoionization; guanine radicals; oxidative damage; time-resolved spectroscopy; quantum chemistry

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Balanikas, E.; Martinez-Fernandez, L.; Baldacchino, G.; Markovitsi, D. Electron Holes in G-Quadruplexes: The Role of Adenine Ending Groups. *Int. J. Mol. Sci.* **2021**, *22*, 13436. https://doi.org/10.3390/ijms222413436

Academic Editors: Janusz Adam Rak, Anil Kumar and Magdalena Zdrowowicz

Received: 31 October 2021
Accepted: 1 December 2021
Published: 14 December 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

Electron holes in DNA bases (radical cations) are generated through electron abstraction reactions with external oxidants or under the effect of ionizing radiation. Regardless of their initial location, electron holes are eventually trapped by guanine sites following a charge migration process [1–9]. Guanine sites constitute hot spots for oxidative damage because their oxidation potential is lower compared to that of other nucleobases [10], and it is further decreased upon stacking [4–6]. It is thus understandable that G-Quadruplexes, which are four-stranded structures characterized by vertical stacking of guanine tetrads (G-tetrads, Figure 1), are prone to such type of damage. Oxidative lesions may perturb the multiple biological functions in which these non-canonical structures play a key role, such as replication, transcription or genome maintenance [11]. Furthermore, G-Quadruplexes are targets for anticancer therapy, which may involve redox reactions [12]. Finally, G-Quadruplexes are extensively studied in view of applications in the field of nanotechnology [13], exploiting, in particular, their behavior as nanowires based on hole transport [14].

In the above-described context, the characterization of electron holes in G-Quadruplexes, (G⁺)•, is particularly important. On the one hand, it is expected to shed light on the primary steps leading to oxidative damage. On the other, it may contribute to the optimization of nanodevices based on charge transport. The first article reporting trapping of electron holes by G-Quadruplexes appeared in 2003 [15]. Later, femtosecond transient absorption experiments studied the efficiency of hole transport through such structures [16,17],

<footnote n="1">¹ LIDYL, CEA, CNRS, Université Paris-Saclay, 91191 Gif-sur-Yvette, France; vangelis.balanikas@cea.fr (E.B.); gerard.baldacchino@cea.fr (G.B.)</footnote>

<footnote n="2">² Departamento de Química, Modúlo 13, Facultad de Ciencias and IADCHEM (Institute for Advanced Research in Chemistry), Campus de Excelencia UAM-CSIC, Universidad Autónoma de Madrid, Cantoblanco, 28049 Madrid, Spain</footnote>

<footnote>* Correspondence: lara.martinez@uam.es (L.M.-F.); dimitra.markovitsi@cea.fr or dimitra.markovitsi@universite-paris-saclay.fr (D.M.)</footnote>

<footnote>† Current affiliation: Institut de Chimie Physique, UMR8000, Université Paris-Saclay, CNRS, 91405 Orsay, France.</footnote>
while direct conductivity measurements were performed by integrating them in electric circuits [14,18,19]. In parallel, theoretical studies examined the role the G-Quadruplex topology on the conductivity and computed the electronic coupling which is the driving force for charge transport [14,20–25].

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Schematic representation of a G-tetrad (a) and the (AG₄A)₄/K⁺ G-Quadruplex (b). Transfer of the red or green proton (a) to the aqueous solvent, concerning the guanine radical cation, gives rise to deprotonated radicals (G-H1)• or (G-H2)•, respectively. In (b), guanines are depicted in gray (tetrads labelled as G1, G2, G3, G4), adenines in violet (at 5', A1) and green (at 3', A2); the backbone is shown in light grey and ions are represented by pink spheres.</figcaption></figure>

A different type of studies focused on the electron holes per se. Their absorption spectra were determined combining quantum chemistry calculations and nanosecond transient absorption spectroscopy [26,27]. Such experiments became possible thanks to the development of specific protocols allowing to study photo-activated processes in room temperature aqueous solutions of DNA multimers. Using low-intensity laser pulses (≤6 × 10⁶ W cm⁻²) at 266 nm, neither the backbone nor the solvent are photo-reactive and, consequently, the transient absorption (ΔA) signals recorded on the microsecond and millisecond time-scales stem, mainly, from nucleobase radicals [28–30] and/or photodimers [31,32]. However, an additional condition is indispensable in order to characterize spectrally the electron holes in G-Quadruplexes: the pH of the solution has to be adjusted to 3 to avoid deprotonation. As a matter of fact, (G⁺)• are unstable in neutral aqueous solution losing a proton to the bulk water [33] (Figure 1a). However, measurements at pH 7 brought precious information regarding the deprotonation dynamics [34,35].

Among the G-Quadruplexes studied by nanosecond transient absorption, (AG₄A)₄/K⁺, formed by association of four AGGGGA strands in the presence of K⁺ ions, has the largest propensity to undergo low-energy photoionization [36]; the term “low-energy” is used in comparison with the ionization potential of the DNA constituents [37]. The quantum yield φ determined for one-photon ionization at 266 nm is 1.5 × 10⁻², which is one order of magnitude higher than that of genomic DNA [29]. Interestingly, the location of the adenine groups at the 5' or the 3' ends was found to play a role in the photoionization process [36]. Although G-Quadruplexes occurring in vivo are formed by folding of a single strand, the tetramolecular systems presented here, in association with other analogs studied previously, constitute convenient models, in particular for time-resolved studies. They allow exploring the effect of peripheral nucleobases on processes triggered by UV radiation. The conclusions drawn from such studies regarding the role of terminal groups could be extrapolated, at a certain extent, to peripheral nucleobases directly attached to the guanine core, as part of loops. Moreover, they may be useful for the development of devices in the field of nanotechnology.

The objective of the present work is, then, to study the electron holes in (AG₄A)₄/K⁺. The first part deals with the spectral features of an electron hole in this system. As the absorp-
tion spectrum determined at pH 3 differs from that reported previously for tetramolecular G-Quadruplexes formed by the oligonucleotide TGGGGT [27], we present additional results obtained for the asymmetric structures (AG₄)₄/K⁺ and (G₄A)₄/K⁺, in which four adenines are located respectively, at the 5' and the 3' ends. Their electronic structure, as well as the electronic absorption spectra, are rationalized by Quantum Mechanical (QM)/Molecular Mechanics (MM) calculations based on the Density Functional Theory (DFT) and its time-dependent version (TD-DFT). In addition to the guanine core, four stacked adenines are included in the QM part; being the first time that all the peripheral groups at 5' and 3' ends are taken into account for the study of electron holes in tetramolecular G-Quadruplexes. For comparison, the dinucleotides 5'-AG-3' and 5'-GA-3' with stacked bases, are also computed. In the second part, we follow the evolution of the time-resolved absorption spectra obtained by photoionization of (AG₄A)₄/K⁺ at pH 7; we discuss the (G⁺)• deprotonation dynamics and present evidence for the formation of a reaction product.

**2. Results and Discussion**

**2.1. *Spectral Properties***

As the intensity of transient signals arising from guanine radicals is weak, spectra
can be recorded only at times longer than 2 µs [29]. The spectrum obtained at 3 µs for
(AG₄A)₄/K⁺ at pH 3 is compared with that reported for the (G⁺)• in deoxyguanosine
monophosphate (dGMP) in Figure 2a [33]. We note that this comparison is limited in the
visible spectral domain because (G⁺)• are very photoactive [38] in the UV. As a result,
the probing light destroys the studied system due to the large number of accumulations
required to obtain the G-Quadruplex spectrum. The monomer spectrum is represented with
its molar absorption coefficient ε, while, in the case of G-Quadruplex, ΔA, determined over
an optical path-length of 1 cm, was divided by the initial concentration of hydrated ejected
electrons [e<sub>hyd</sub><sup>-</sup>]₀, which is equal to that of generated electron holes. The methodology
followed in order to determine [e<sub>hyd</sub><sup>-</sup>]₀ is described in reference [29].

<figure><img src="image_6.png" /></figure>

**Figure 2.** Experimental absorption spectra (circles) determined at 3 µs by photoionization of (a) the symmetric structure (AG₄A)₄/K⁺ (red) and (b) the two asymmetric analogues (G₄A)₄/K⁺ (green) and (AG₄)₄/K⁺ (blue) at pH 3; ΔA was divided by the initial concentration of hydrated ejected electrons. The black line in (a) corresponds to the absorption spectrum reported for the radical cation of dGMP [39], plotted with its molar absorption coefficient.
We observe in Figure 2a that the intensity of the G-Quadruplex spectrum at 400 nm is significantly lower compared to that of the monomer. A similar trend was reported for the spectra of electron holes in monomolecular telomeric [26] and tetramolecular TGGGGT [27] G-Quadruplexes. However, for these two systems, the absorption band in the visible was found to practically overlap that of the monomeric (G⁺)⁺. This is not obviously the case for (AG₄A)₄/K⁺: although its band exhibits, within the experimental error (±5%), the same intensity as that of the monomeric (G⁺)⁺, its maximum is located at ~35 nm longer wavelengths (Figure 2a). Yet, the guanine cores of the tetramolecular structures formed by association of TGGGGT [40] and AGGGGA [36] strands have the same parallel topology; the cation spectra of these two systems are compared in Figure S1. We, therefore, conclude that it is the presence of adenines that alters the spectral properties of the electron hole.

The role of adenines on the spectrum of the electron hole was further checked by performing similar experiments on the asymmetric structures (AG₄)₄/K⁺ and (G₄A)₄/K⁺. The spectra of these two photoionized systems at pH 3 are clearly different (Figure 2b), that of (AG₄)₄/K⁺ being somewhat closer to the monomer spectrum. The spectrum of (AG₄A)₄/K⁺ appears, thus, to be an average of spectra of the two asymmetric G-Quadruplexes.

In order to understand the origin of the adenine effect on the above presented spectra of electron holes, we performed theoretical calculations on (G₄A)₄/K⁺ and (AG₄)₄/K⁺ containing an electron hole. As a first step, we optimized their geometry. In this respect, we recall that two-dimensional NOESY NMR measurements on (AG₄A)₄/K⁺ did not detect any hydrogen bonding between adenines, meaning that no stable adenine tetrads are formed, but adenines could be intermittently stacked on the guanines [36]. Consequently, we considered the limiting case in which adenines are stacked on the guanine core since such an arrangement would maximize the electronic interactions between the two types of bases. Our computations showed that the electron holes are located at guanine sites for both the model dinucleotides (Figure S3) and the G-Quadruplexes (Figure 3). Although this is expected in view of the oxidation potentials of the bases, it is shown for the first time that this tendency is preserved even within G-Quadruplexes with stacked adenine nucleobase. Moreover, the cation is almost perfectly localized on a single guanine. This is surprising because theoretical studies reported hole delocalization for other parallel G-Quadruplexes [41], or even for simple parallel guanine stacks [42]. The structure of the optimized guanine cation within the systems studied here has the geometrical features reported for the monomeric species [42–45] and also found for previous G-Quadruplexes [26,27,41,46,47]. Despite the important shortening of the C2-NH2 bond in the guanine cation, the spin density has contributions from almost the entire ring of the nucleobase (Figure 3).

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Spin densities for the optimized radical cation in the systems (a) (G₄A)₄/K⁺ and (b) (AG₄)₄/K⁺. MM guanines, backbone and ions are depicted in gray, QM guanines in blue, adenines in violet (at 5') and green (at 3').</figcaption></figure>
An important finding is that the electron hole in each asymmetric structure is not located on the same G-tetrad (Figure 3). In $(G_4A)_4/K^+$, it is located within the G4 tetrad that is stacked with the 3' adenines (A2), whereas in $(AG_4)_4/K^+$, it is located in G2, corresponding to an inner tetrad, far from the 5'-adenines (A1). However, the relative energy difference between the two species is low, <1 kcal mol⁻¹, suggesting that both are likely in the symmetric $(AG_4A)_4/K^+$ structure with perfectly stacked adenines.

The existence of a radical cation in the guanine core of a G-Quadruplex, as part of an excited charge transfer (CT) state, is considered to be the primary step leading to low-energy photoionization [29,30,36]. We can rationalize our present findings in analogy with those previously reported regarding the nature and stability of CT states. It is nowadays well assessed, not only for G-Quadruplexes, but also for other DNA systems, that CT states are more stable when the electron moves along the 5' → 3' channel [36,48–50]. In the first place, this agrees with the “G-tetrad” in which the hole is located. If the system has a stacked 5'-GA-3' sequence, the cation is localized close to the adenine connected with an initial stable 5'-G⁺ → A⁻-3' CT state/minimum. However, adenines at the 5' end would involve the formation of less stable 5'-A⁺ → G⁻-3' (3'-G⁺ → A⁻-5) or 5'-G⁺ → G⁻-3' CT states; thus, the cation tends to be localized at an inner tetrad. This matches our previous results about excited states in stacked AG/GA dinucleotides, which revealed that the coupling between the bright and CT states is larger in 5'-GA-3' compared to 5'-AG-3' [36]. The loss of stability when the non-guanine bases are positioned at the 5' end is even more dramatic when adenine is replaced by thymine. For instance, non-CT minimum was located for the 5'-TG-3' dinucleotide [36]. This could be related with the tendency towards some extent of delocalization between guanines for holes in tetramolecular G-Quadruplexes with thymine terminal groups [47].

In a subsequent step, the absorption spectra of the two "cationic" G-Quadruplexes were modeled by computing the energies and intensities of the one hundred lowest-energy excited states obtained by TD-DFT (Figure 4). At wavelengths longer than 600 nm, the intensity found for the spectrum of the $(G_4A)_4/K^+$ "cation" is higher compared to the $(AG_4)_4/K^+$ analogue. A close inspection of the lowest-energy excited states reveals the potential background for this result. For the $(AG_4)_4/K^+$ "cation", the energies of these excited states are blue-shifted up to 0.2 eV and the global oscillator strength (f) is half the one summed up compared to $(G_4A)_4/K^+$ "cation". As discussed above, the most intense absorption of the $(G_4A)_4/K^+$ "cation" in the red part of the spectrum stems from the directionality of the CT excited states: larger stability (i.e., lower energy) for the CT states in which the electron moves towards the 3' end. The broad bands located in the 400–600 nm region follow similar trend in terms of energies, i.e., slightly red-shifted for the $(G_4A)_4/K^+$ "cation". However, the behavior of f is reversed compared to the red region, the $(G_4A)_4/K^+$ "cation" spectrum being slightly less intense in the 500–550 nm region; although several excited states are present in this energy range, their computed f is at most 0.04. In contrast, the $(AG_4)_4/K^+$ "cation" has two states peaking at these energies, one of them being particularly bright (f = 0.08).

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Absorption spectra computed for the optimized radical cation in $(G_4A)_4/K^+$ (green) and $(AG_4)_4/K^+$ (cyan) at the PCM/TD-DFT(M052X)/6-31G(d)//MM level of theory.</figcaption></figure>
Similar results were also obtained for the model dinucleotides (Figure S4): the lowest in energy broad band is red-shifted and less intense for 5'-GA-3' compared to 5'-AG-3' but these differences are smaller compared to those observed for G-Quadruplexes. This is related to the further stabilization (beyond directionality) of lower in energy CT excited states due to the presence of metal inner cations coordinated within G-Quadruplex systems [41,51–53].

## 2.2. Dynamical Properties

The transient absorption spectra determined for (AG₄A)₄/K⁺ at pH 7 at selected times, ranging from 3 µs to 10 ms, are shown in Figure 5. The main evolution observed at early times is a shift of the band peaking in the visible toward longer wavelengths (Figure 5a). A similar trend was found for all the G-Quadruplexes studied so far [29,35]. It is attributed to the deprotonation of (G⁺)•, losing the amino proton at position 2, which is not engaged in hydrogen bonding [34] (see Figure 1a). As a result, the spectral feature of the deprotonated radical (G-H2)•, absorbing at around 600 nm (Figure 5b), becomes dominant.

<figure><img src="image_9.png" /></figure>

**Figure 5.** Transient absorption spectra determined for aerated solutions of (AG₄A)₄/K⁺ at pH 7 (circles) at 3 µs (red), 50 µs (pink), 0.2 ms (green), 0.5 ms (blue), 2 ms (cyan), 6 ms (dark red) and 10 ms (yellow). In (a), the experimentally measured ΔA is shown while in (b), ΔA at 605 nm has been normalized to 1; black dots correspond to the (G-H2)• spectrum reported for monomeric guanosine [54].

As already described for other four-stranded systems [35], the (G⁺)• decay in (AG₄A)₄/K⁺ is highly anisotropic, spanning from at least 30 ns to a few tens of µs. This is reflected on the rapid rise of the transient absorption signal observed on the sub-microsecond time-scale in the red part of the spectrum (Figure 6a), where the absorbance of (G-H2)• radicals is higher than that of (G⁺)•. These experiments were performed by saturating the solution with N₂O, an efficient scavenger of hydrated electrons [55], whose strong absorption blurs that of radicals [30]. However, as the sensitivity of the experimental setup is much lower at this time-domain, the transient signals are very noisy, despite the large number of accumulations. The rise in Figure 6a can be approximated by a single exponential function with a time-constant of 340 ± 40 ns.
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Transient absorption signals recorded for (AG₄A)₄/K⁺ at pH 7: (a) 605 nm (pink); N₂O saturated solutions; (b) 290 nm (violet), 380 nm (blue), 530 nm (green) and 620 nm (dark red); aerated solutions. Yellow lines correspond to fits with mono-exponential functions: {A₁ (1 - exp(t/τ₁)) + A₀ (rise)} and {A₁ exp(-t/τ₂) + A₀} (decays); τ₁ = 340 ± 40 ns and τ₂ = 4.6 ± 0.4 µs.</figcaption></figure>

The (G⁺)• population continues to decrease on the microsecond time-scale, as indicated by the signals recorded at wavelengths where its absorbance is preponderant: around its maximum (Figure 2a) and also below 300 nm, as expected in analogy with the spectrum of the monomeric cation (inset in Figure 2a). In this case, the signal has been approximated by a mono-exponential function (Figure 6b). We stress that the parameters derived from such fits do not correspond to reaction rate constants, as explained in detail in reference [35], but are simply used for a phenomenological description of the signals. In the red part of the spectrum, dominated by the absorbance of (G-H2)•, the transient signals appear to be flat. This could be explained by the superposition of a rise due to (G-H2)• formation, and a decay due to the disappearance of (G⁺)•, also absorbing at this region. Similar quasi-flat signals are also observed around 350–400 nm.

In Figure 5b, we compare the evolution of the spectral shape recorded on the millisecond time-scale by normalizing the transient absorption around 600 nm. We observe that between 0.2 ms and 6 ms, the spectral shape in the 500–700 nm domain remains practically the same. In contrast, the relative intensity at shorter wavelengths, featuring a peak around 335 nm and a shoulder at ~380 nm, increases. We, therefore, conclude that, in addition to (G-H2)•, another species, noted as X, is formed. This is also attested by the transient absorption signals recorded on the millisecond time-scale: at all wavelengths above 500 nm, they practically overlap, while they follow a much slower decay pattern at shorter wavelengths (Figure 7). Moreover, we remark in Figure 7, that the decay of (G-H2)• in (AG₄A)₄/K⁺ is similar to that observed for (TG₄T)₄/K⁺ [46], showing that their reaction rate does not depend on the nature of the terminal groups.

For several G-Quadruplexes studied previously, (G-H2)• → (G-H1)• tautomerization was observed on the millisecond time-scale [29]. This transformation is particularly clear in the case of monomolecular TEL21/Na⁺ and tetramolecular (TG₄T)₄/Na⁺ [26,27]. However, comparison of the spectrum obtained for (AG₄A)₄/K⁺ at 10 ms with that of monomeric (G-H1)• in Figure 8, shows that tautomerization of deprotonated radicals does not constitute a major reaction pathway. At this point, we conclude that X is a reaction product or a reaction intermediate formed at earlier times and not correlated with (G-H2)• radicals. The transient absorption experiments performed in the frame of this study dedicated to electron holes cannot unveil the nature of X. However, they provide some hints which may be useful for future works.
<figure><img src="image_11.png" /><figcaption>**Figure 7.** Normalized transient absorption decays determined for aerated solutions of (AG₄A)₄/K⁺ at 385 nm (violet), 500 nm (blue) and 685 nm (red). For comparison, the decay obtained for (TG₄T)₄/K⁺ at 605 nm [46] is shown in black.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 8.** Comparison of the transient absorption spectrum determined for (AG₄A)₄/K⁺ at 10 ms (yellow circles) with that reported for monomeric (G-H1)• (black dots) [39].</figcaption></figure>

So far, the fingerprints of two types of UV-induced reactions in DNA have been detected in the 300–500 nm spectral region by nanosecond transient absorption spectroscopy. On the one hand, they are correlated with dimerization of nucleobases [31,32] occurring on a singlet excited state, i.e., at times much shorter than those explored by such experiments (see Figure 6). The initially formed photoproducts may also undergo further dark reactions on the millisecond time-scale. In all cases, the intensity of the transient signals is proportional to the concentration of the photons absorbed by the system [hv]. On the other hand, we may have reaction products or reaction intermediates resulting from photoionization. In this case, the intensity of their transient absorption is proportional to the concentration of electron holes generated within the system. Moreover, the latter is equal to the concentration of hydrated ejected electrons [e<sub>hyd</sub><sup>-</sup>]₀ which can be determined by independent experiments [29]. In the range of the excitation intensities used in the studies of low-energy DNA photoionization, e<sub>hyd</sub><sup>-</sup> results from both one- and two-photon processes. Thus, the behavior of the transient absorption observed upon changing the ratio [e<sub>hyd</sub><sup>-</sup>]₀/[hv], achieved by varying the energy of the incident laser pulses, may help disentangling the two types of reaction.

The transient absorption spectrum recorded for (AG₄A)₄/K⁺ at 30 ms exhibits a striking similarity with that detected under the same conditions for TEL21/K⁺ [56] (Figure 9a). For the latter system, it was shown that it corresponds to a reaction intermediate decaying on hundreds of ms and stemming from photoionization because the pattern of transient absorption at 365 nm, where both X and (G-H2)• absorb, is not altered upon varying [e<sub>hyd</sub><sup>-</sup>]₀/[hv] by 50%. In the present study, the variation was limited to 20%, due to the lower contribution of the two-photon process to the (AG₄A)₄/K⁺ photoionization, as attested by the respective ionization curves [36,56]. Unfortunately, much larger variations of
the laser pulse energy either lead to very noisy signals or provoke two-photons ionization of the aqueous solvent damaging the DNA [57]. Despite this limitation, we estimate that it would be possible to detect changes in the decay pattern if X had a different origin than (G-H2)•. However, this is not the case as shown in Figure 9b. Therefore, by analogy with the TEL21/K⁺ behavior, X could be a reaction intermediate resulting from electron holes.

<figure><img src="image_13.png" /><figcaption>**Figure 9.** (a) Comparison of the transient absorption spectra determined for (AG₄A)₄/K⁺ (yellow circles) and TEL21/K⁺ (dark blue circles) at 30 ms; (b) Transient signals recorded for (AG₄A)₄/K⁺ at 365 nm with incident laser pulses: 3 mJ (green) and 6 mJ (pink).</figcaption></figure>

Several types of DNA lesions are known to result from (G⁺)• following reaction paths other than deprotonation; these include adducts with non-DNA species, as well as 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodG) [58]. Numerous publications discuss 8-oxodG formed in G-Quadruplexes through electron abstraction reactions with external oxidants (see, for example, references [59,60]). The formation of 8-oxodG was also explored for G-Quadruplexes that have been simply irradiated at 266 nm, in the absence of external oxidants. In the case of (TG₄T)₄/Na⁺ [27], which has the same topology as (AG₄A)₄/K⁺, no 8-oxodG was found using high performance liquid chromatography combined to mass spectrometry. In contrast, the same technique did detect 8-oxodG in photoionized TEL21/Na⁺ [26]. However, the time-resolved spectra recorded for latter system [26] do not exhibit the spectral features observed for (AG₄A)₄/K⁺ and TEL21/K⁺ (Figure 9), suggesting that a different reaction takes place in these systems.

**3. Materials and Methods**

**3.1. Materials**

Oligonucleotides were purchased from Eurogentec Europe. As the purity is an important requirement in radical studies, they were purified by reversed phase HPLC and tested by MALDI-TOF. The phosphate buffer was composed of an equimolar mixture of KH₂PO₄/K₂HPO₄ in concentrations 0.15 or 0.015 mol·L⁻¹ each. Solutions were prepared using ultra-pure water, delivered by a Milli-Q Integral system. The pH of the buffer, measured by a HANNA Instr. Apparatus (pH 210), was adjusted to 7 by addition of a concentrated KOH solution. For measurements in pH 3, a concentrated solution of phosphoric acid (1 mol·L⁻¹) was added in the solution. The purity of the buffer salts was 99.99%. For the formation of G-Quadruplexes, an appropriate mother solution was heated to 96 °C during 5 min, then slowly cooled to 4 °C, where it was incubated overnight. Prior to time-resolved experiments, the absorbance at 266 nm was adjusted to 2.0 ± 0.1 over 1 cm, heated to 50 °C, so that to destroy higher-order aggregates possibly formed at high concentration, and then cooled slowly to room temperature; the G-Quadruplex concentration was ~10⁻⁵ mol·L⁻¹. G-Quadruplex formation was checked par circular dichroism (see SI in reference [36]).

During the experiments, 2 mL of solution, contained in a 1 cm × 1 cm quartz cell, were mildly stirred while the temperature was maintained at 23 °C. To avoid re-excitation of damaged G-Quadruplexes, solutions were replaced frequently.
## 3.2. *Transient Absorption Spectroscopy*

The transient absorption setup used as excitation source was the fourth harmonic of a nanosecond Nd:YAG laser (Quanta Ray, Spectra-Physics, Santa Clara, CA, USA). The excited area at the surface of the sample was 0.6 × 1.0 cm². The analyzing beam, orthogonal to the exciting beam, was provided by a 150 W Xe-arc lamp (OSRAM XBO, Munich, Germany); its optical path length through the sample was 1 cm while its thickness was limited to 0.1 cm in order to use the most homogeneous part of the light. The probed volume was located in the very first 0.1 cm part of the cell along the propagation of the exciting laser beam delimited by four dedicated slits. Then, the analyzing light was dispersed in a SPEX 270M monochromator (Horiba Jobin-Yvon, Longjumeau, France), detected by a Hamamatsu R928 photomultiplier (Hamamatsu Photonics France, Massy, France) and recorded by a Waverunner 6084 oscilloscope (Lecroy, Teledyn Lecroy, Courtaboeuf, France). For measurements on the sub-microsecond scale, the Xe-arc lamp was intensified via an electric discharge (Applied Photophysics, Leatherhead, UK). Transient absorption spectra were recorded using a wavelength-by-wavelength approach. The excitation rate was 0.2 Hz. The incident pulse energy at the surface of the sample was measured using a NIST traceable pyroelectric sensor (Nova2/PE25, Ophir Spiricon Europe GmbH, Darmstadt, Germany).

## 3.3. *QM/MM Calculations*

Our QM/MM ONIOM [61] models for the (G₄A)₄ and (AG₄)₄ systems are depicted in Figure 3. We included, in the QM part, the involved adenine tetrad (at 3' or 5' end) plus the two closest G-tetrads and inner ions. The rest of the system, i.e., the two remaining G tetrads, the backbones and the outer counterions are considered just at the MM level of theory. This allows for a balanced description of the asymmetric species with no dramatic effect on the selected geometry. The geometry of the "cationic" G-Quadruplexes was optimized at the QM/MM level, selecting density functional theory (DFT) and, in particular, the M052X [62,63] functional and the 6-31G(d) basis set for the QM region and amber force field parm96.dat [64] for the MM part. The entire system is considered to be embedded in water by using an implicit polarization method (Polarizable Continuum Model, PCM) [65,66]. For modeling the absorption spectra, we resort to the time-dependent version of DFT with the same settings above specified (TD-DFT) to compute the vertical absorption energies and oscillator strength. Then, these energies are uniformly shifted by -0.6 eV and convoluted with a Gaussian function with a hwm of 0.3 eV. This computational approach has already been shown to satisfactorily reproduce the experimental absorption and circular dichroism spectra for other G-Quadruplex systems (see, for instance, reference [26,36,46,47,53]). All the calculations were performed using Gaussian16 program [67].

# 4. Main Conclusions

One important outcome derived from the computational part of our study is that the guanine site in which the electron hole is localized within the examined tetramolecular G-Quadruplexes is modulated by the stacking of adenines at the 3' or 5' ends. Although the energy difference determined for the two asymmetric models used in the computations is small, the situation may change in real systems under the effect of conformational motions and the intermittent stacking (and probably not symmetric) of the terminal groups. In a more general way, this is the first demonstration that peripheral bases may affect the position of the electron hole within the G-Quadruplex core. It would be interesting to pursue such investigations on monomolecular systems, more relevant for biological processes than the models examined here.

Peripheral adenine groups affect, in a different way, the absorption spectra of the electron hole, depending on their position in the 5' or 3' ends. Their knowledge could be useful for the characterization of excited CT by femtosecond spectroscopy [68], allowing to distinguish between G⁺ → G⁻ and G⁺ → A⁻ CT states. Such a disentanglement is important in order to better understand the mechanism responsible for low-energy photoionization of DNA [29,36]. Furthermore, it would also be interesting to explore by
pulse radiolysis if similar effects are encountered for the corresponding "anionic" species formed by electron attachment [57,69].

As for the previously studied G-Quadruplexes, the decay of (G⁺)• in (AG₄A)₄/K⁺ spans from at least 30 ns to tens of µs. It would not be surprising that such a highly anisotropic reaction rate arises, at least partly, from a variety of positions of electron holes within the G-Quadruplex core. (G⁺)• deprotonation in (AG₄A)₄/K⁺ gives rise to (G-H2)• deprotonated radicals, whose reaction products have not been identified so far. Our experimental results indicate that the peripheral groups are unlikely to be involved in the (G-H2)• reactions because their decay pattern is the same for (AG₄A)₄/K⁺ and (TG₄T)₄/K⁺. Moreover, for both systems, (G-H2)• → (G-H1)• tautomerization does not constitute a noticeable reaction path.

Finally, our time-resolved experiments brought evidence for the formation of a re-action product, similar to that observed for the telomeric G-Quadruplex TEL21/K⁺, and probably stemming from (G⁺)•. Given that its spectral fingerprint decays, although more slowly than that of (G-H2)•, it should not correspond to a final lesion but rather to a reaction intermediate. However, the spectral and dynamical features presented here are not sufficient for its identification. We hope that our findings will motivate experts of analytical chemistry to focus on this problem.

**Supplementary Materials:** Supplementary Materials can be found at https://www.mdpi.com/article/10.3390/ijms222413436/s1.

**Author Contributions:** Conceptualization, D.M. and L.M.-F.; methodology, D.M., E.B. and L.M.-F.; validation, E.B., G.B. and D.M.; formal analysis, E.B.; investigation, E.B. and L.M.-F.; writing—original draft preparation, D.M. and L.M.-F.; review and editing, E.B. and L.M.-F., G.B. and D.M.; supervision, D.M. and G.B.; project administration and funding acquisition, D.M. and L.M.-F. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work has received funding from the European Union's Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement number 765266 (LightDy-NAmics).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

**Acknowledgments:** L.M.F thanks the Proy. PID2019-110091GB-I00 funded by MCIN/AEI/10.13039/501100011033 and the Centro de Computación Científica UAM (CCC-UAM) for computing time. We thank Peter Podbevsek for helpful discussions.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Hall, D.B.; Holmlin, R.E.; Barton, J.K. Oxidative DNA damage through long-range electron transfer. *Nature* **1996**, *382*, 731–735. [CrossRef]

2. Douki, T.; Angelov, D.; Cadet, J. UV laser photolysis of DNA: Effect of duplex stability on charge-transfer efficiency. *J. Am. Chem. Soc.* **2001**, *123*, 11360–11366. [CrossRef]

3. Stemp, E.D.A.; Arkin, M.R.; Barton, J.K. Oxidation of guanine in DNA by Ru(phen)(2)(dppz)(3+) using the flash-quench technique. *J. Am. Chem. Soc.* **1997**, *119*, 2921–2925. [CrossRef]

4. Meggers, E.; Michel-Beyerle, M.E.; Giese, B. Sequence dependent long range hole transport in DNA. *J. Am. Chem. Soc.* **1998**, *120*, 12950–12955. [CrossRef]

5. Saito, I.; Nakamura, T.; Nakatani, K.; Yoshioka, Y.; Yamaguchi, K.; Sugiyama, H. Mapping of the hot spots for DNA damage by one-electron oxidation: Efficacy of GG doublets and GGG triplets as a trap in long-range hole migration. *J. Am. Chem. Soc.* **1998**, *120*, 12686–12687. [CrossRef]

6. Kanvah, S.; Joseph, J.; Schuster, G.B.; Barnett, R.N.; Cleveland, C.L.; Landman, U. Oxidation of DNA: Damage to nucleobases. *Acc. Chem. Res.* **2010**, *43*, 280–287. [CrossRef] [PubMed]

7. Lewis, F.D.; Letsinger, R.L.; Wasielewski, M.R. Dynamics of photoinduced charge transfer and hole transport in synthetic DNA hairpins. *Acc. Chem. Res.* **2001**, *34*, 159–170. [CrossRef] [PubMed]

8. Kawai, K.; Majima, T. Hole transfer kinetics of DNA. *Acc. Chem. Res.* **2013**, *46*, 2616–2625. [CrossRef]
9. Ma, J.; Marignier, J.L.; Pernot, P.; Houee-Levin, C.; Kumar, A.; Sevilla, M.D.; Adhikary, A.; Mostafavi, M. Direct observation of the oxidation of DNA bases by phosphate radicals formed under radiation: A model of the backbone-to-base hole transfer. *Phys. Chem. Chem. Phys.* **2018**, *20*, 14927–14937. [CrossRef]

10. Palecek, E.; Bartosik, M. Electrochemistry of nucleic acids. *Chem. Rev.* **2012**, *112*, 3427–3481. [CrossRef] [PubMed]

11. Varshney, D.; Spiegel, J.; Zyner, K.; Tannahill, D.; Balasubramanian, S. The regulation and functions of DNA and RNA G-quadruplexes. *Nat. Rev. Mol. Cell Biol.* **2020**, *21*, 459–474. [CrossRef] [PubMed]

12. Kizhuveetil, U.; Omer, S.; Karunagaran, D.; Suraishkumar, G.K. Improved redox anti-cancer treatment efficacy through reactive species rhythm manipulation. *Sci. Rep.* **2020**, *10*, 1588. [CrossRef]

13. Mergny, J.L.; Sen, D. DNA Quadruple helices in nanotechnology. *Chem. Rev.* **2019**, *119*, 6290–6325. [CrossRef] [PubMed]

14. Livshits, G.I.; Stern, A.; Rotem, D.; Borovok, N.; Eidelshtein, G.; Migliore, A.; Penzo, E.; Wind, S.J.; Di Felice, R.; Skourtis, S.S.; et al. Long-range charge transport in single G-quadruplex DNA molecules. *Nat. Nanotechnol.* **2014**, *9*, 1040–1046. [CrossRef] [PubMed]

15. Delaney, S.; Barton, J.K. Charge transport in DNA duplex/quadruplex conjugates. *Biochemistry* **2003**, *42*, 14159–14165. [CrossRef] [PubMed]

16. Choi, J.; Park, J.; Tanaka, A.; Park, M.J.; Jang, Y.J.; Fujitsuka, M.; Kim, S.K.; Majima, T. Hole trapping of G-quartets in a G-quadruplex. *Angew. Chem. Int. Ed.* **2013**, *52*, 1134–1138. [CrossRef]

17. Thazhathveetil, A.K.; Harris, M.A.; Young, R.M.; Wasielewski, M.R.; Lewis, F.D. Efficient charge transport via DNA G-quadruplexes. *J. Am. Chem. Soc.* **2017**, *139*, 1730–1733. [CrossRef]

18. Liu, S.P.; Weisbrod, S.H.; Tang, Z.; Marx, A.; Scheer, E.; Erbe, A. Direct measurement of electrical transport through G-quadruplex DNA with mechanically controllable break junction electrodes. *Angew. Chem. Int. Ed.* **2010**, *49*, 3313–3316. [CrossRef]

19. Sha, R.J.; Xiang, L.M.; Liu, C.R.; Balaeff, A.; Zhang, Y.Q.; Zhang, P.; Li, Y.Q.; Beratan, D.N.; Tao, N.J.; Seeman, N.C. Charge splitters and charge transport junctions based on guanine quadruplexes. *Nat. Nanotechnol.* **2018**, *13*, 316–321. [CrossRef] [PubMed]

20. Wang, X.; Liang, S.D. Topological effects of charge transfer in telomere G-quadruplex: Mechanism on telomerase activation and inhibition. *Inter. J. Mod. Phys. B* **2013**, *27*, 1350001. [CrossRef]

21. Calzolari, A.; Di Felice, R.; Molinari, E.; Garbesi, A. Electron channels in biomolecular nanowires. *J. Phys. Chem. B* **2004**, *108*, 2509–2515. [CrossRef]

22. Lech, C.J.; Anh Tuan, P.; Michel-Beyerle, M.-E.; Voityuk, A.A. Electron-hole transfer in G-quadruplexes with different tetrad stacking geometries: A combined QM and MD study. *J. Phys. Chem. B* **2013**, *117*, 9851–9856. [CrossRef]

23. Kang, D.W.; Sun, M.L.; Zuo, Z.W.; Wang, H.X.; Lv, S.J.; Li, X.Z.; Li, L.B. Charge transport and magnetoresistance of G4-DNA molecular device modulated by counter ions and dephasing effect. *Phys. Lett. A* **2016**, *380*, 977–982. [CrossRef]

24. Sun, W.M.; Varsano, D.; Di Felice, R. Effects of G-quadruplex topology on electronic transfer integrals. *Nanomaterials* **2016**, *6*, 6100184. [CrossRef]

25. Ravindranath, R.; Mondal, P.; Gillet, N. Radical cation transfer in a guanine pair: An insight to the G-quadruplex structure role using constrained DFT/MM. *Theor. Chem. Acc.* **2021**, *140*, 89. [CrossRef]

26. Banyasz, A.; Martinez-Fernandez, L.; Balty, C.; Perron, M.; Douki, T.; Improta, R.; Markovitsi, D. Absorption of low-energy UV Radiation by human telomere g-quadruplexes generates long-lived guanine radical cations. *J. Am. Chem. Soc.* **2017**, *139*, 10561–10568. [CrossRef] [PubMed]

27. Banyasz, A.; Balanikas, E.; Martinez-Fernandez, L.; Baldacchino, G.; Douki, T.; Improta, R.; Markovitsi, D. Radicals generated in tetramolecular guanine quadruplexes by photo-ionization: Spectral and dynamical features. *J. Phys. Chem. B* **2019**, *123*, 4950–4957. [CrossRef] [PubMed]

28. Marguet, S.; Markovitsi, D.; Talbot, F. One and two photon ionization of DNA single and double helices studied by laser flash photolysis at 266 nm. *J. Phys. Chem. B* **2006**, *110*, 11037–11039. [CrossRef]

29. Balanikas, E.; Banyasz, A.; Douki, T.; Baldacchino, G.; Markovitsi, D. Guanine radicals induced in DNA by low-energy photoionization. *Acc. Chem. Res.* **2020**, *53*, 1511–1519. [CrossRef] [PubMed]

30. Balanikas, E.; Markovitsi, D. DNA photoionization: From high to low energies. In *DNA Photodamage*; Improta, R., Douki, T., Eds.; RSC: Cambridge, UK, 2021; pp. 37–54.

31. Marguet, S.; Markovitsi, D. Time-resolved study of thymine dimer formation. *J. Am. Chem. Soc.* **2005**, *127*, 5780–5781. [CrossRef]

32. Banyasz, A.; Martinez-Fernandez, L.; Ketola, T.; Muñoz-Losa, A.; Esposito, L.; Markovitsi, D.; Improta, R. Excited state pathways leading to formation of adenine dimers. *J. Phys. Chem. Lett.* **2016**, *7*, 2020–2023. [CrossRef]

33. Candeias, L.P.; Steenken, S. Ionization of purine nucleosides and nucleotides and their components by 193-nm laser photolysis in aqueous solution: Model studies for oxidative damage of DNA. *J. Am. Chem. Soc.* **1992**, *114*, 699–704. [CrossRef]

34. Wu, L.D.; Liu, K.H.; Jie, J.L.; Song, D.; Su, H.M. Direct observation of guanine radical cation deprotonation in G-quadruplex DNA. *J. Am. Chem. Soc.* **2015**, *137*, 259–266. [CrossRef]

35. Balanikas, E.; Banyasz, A.; Baldacchino, G.; Markovitsi, D. Deprotonation dynamics of guanine radical cations. *Photochem. Photobiol.* **2021**, *13*, 540. [CrossRef]

36. Balanikas, E.; Martinez-Fernadez, L.; Improta, R.; Podbevsek, P.; Baldacchino, G.; Markovitsi, D. The structural duality of nucleobases in guanine quadruplexes controls their low-energy photoionization. *J. Phys. Chem. Lett.* **2021**, *12*, 8309–8313. [CrossRef] [PubMed]

37. Pluharova, E.; Slavicek, P.; Jungwirth, P. Modeling photoionization of aqueous DNA and its components. *Acc. Chem. Res.* **2015**, *48*, 1209–1217. [CrossRef]
38. Khanduri, D.; Adhikary, A.; Sevilla, M.D. Highly oxidizing excited states of one-electron-oxidized guanine in DNA: Wavelength and pH dependence. *J. Am. Chem. Soc.* **2011**, *133*, 4527–4537. [CrossRef] [PubMed]

39. Candeias, L.P.; Steenken, S. Structure and acid-base properties of one-electron-oxidized deoxyguanosine, guanosine, and 1-methylguanosine. *J. Am. Chem. Soc.* **1989**, *111*, 1094–1099. [CrossRef]

40. Sket, P.; Plavec, J. Tetramolecular DNA quadruplexes in solution: Insights into structural diversity and cation movement. *J. Am. Chem. Soc.* **2010**, *132*, 12724–12732. [CrossRef] [PubMed]

41. Martinez-Fernandez, L.; Banyasz, A.; Markovitsi, D.; Improta, I. Topology controls the electronic absorption delocalization of electron hole in guanine quadruplexes. *Chem. Eur. J.* **2018**, *24*, 15185–15189. [CrossRef] [PubMed]

42. Kumar, A.; Sevilla, M.D. Density functional theory studies of the extent of hole delocalization in one-electron oxidized adenine and guanine base stacks. *J. Phys. Chem. B* **2011**, *115*, 4990–5000. [CrossRef]

43. Adhikary, A.; Malkhasian, A.Y.S.; Collins, S.; Koppen, J.; Becker, D.; Sevilla, M.D. UVA-visible photo-excitation of guanine radical cations produces sugar radicals in DNA and model structures. *Nucl. Ac. Res.* **2005**, *33*, 5553–5564. [CrossRef]

44. Adhikary, A.; Kumar, A.; Becker, D.; Sevilla, M.D. The guanine cation radical: Investigation of deprotonation states by ESR and DFT. *J. Phys. Chem. B* **2006**, *110*, 24171–24180. [CrossRef] [PubMed]

45. Cheng, P.; Li, Y.; Li, S.; Zhang, M.; Zhou, Z. Collision-induced dissociation (CID) of guanine radical cation in the gas phase: An experimental and computational study. *Phys. Chem. Chem. Phys.* **2010**, *12*, 4667–4677. [CrossRef] [PubMed]

46. Behmand, B.; Balanikas, E.; Martinez-Fernandez, L.; Improta, R.; Banyasz, A.; Baldacchino, G.; Markovitsi, D. Potassium ions enhance guanine radical generation upon absorption of low-energy photons by G-quadruplexes and modify their reactivity. *J. Phys. Chem. Lett.* **2020**, *11*, 1305–1309. [CrossRef] [PubMed]

47. Martinez-Fernandez, L.; Esposito, L.; Improta, R. Studying the excited electronic states of guanine rich DNA quadruplexes by quantum mechanical methods: Main achievements and perspectives. *Photochem. Photobiol. Sci.* **2020**, *19*, 436–444. [CrossRef]

48. Santoro, F.; Barone, V.; Lami, A.; Improta, R. The excited electronic states of adenine-guanine stacked dimers in aqueous solution: A PCM/TD-DFT study. *Phys. Chem. Chem. Phys.* **2010**, *12*, 4934–4948. [CrossRef] [PubMed]

49. Aquino, A.J.A.; Nachtigallova, D.; Hobza, P.; Truhlar, D.G.; Hattig, C.; Lischka, H. The charge-transfer states in a stacked nucleobase dimer complex: A benchmark study. *J. Comput. Chem.* **2011**, *32*, 1217–1227. [CrossRef]

50. Martinez-Fernandez, L.; Improta, R. Novel adenine/thymine photodimerization channels mapped by PCM/TD-DFT calculations on dApT and TpdA dinucleotides. *Photochem. Photobiol. Sci* **2017**, *16*, 1277–1283. [CrossRef] [PubMed]

51. Hua, Y.; Changenet-Barret, P.; Improta, R.; Vayá, I.; Gustavsson, T.; Kotlyar, A.B.; Zikich, D.; Šket, P.; Plavec, J.; Markovitsi, D. Cation effect on the electronic excited states of guanine nanostructures studied by time-resolved fluorescence spectroscopy. *J. Phys. Chem. C* **2012**, *116*, 14682–14689. [CrossRef]

52. Improta, R. Quantum mechanical calculations unveil the structure and properties of the absorbing and emitting excited electronic states of guanine quadruplex. *Chem. Eur. J.* **2014**, *20*, 8106–8115. [CrossRef] [PubMed]

53. Martinez-Fernandez, L.; Changenet, P.; Banyasz, A.; Gustavsson, T.; Markovitsi, D.; Improta, R. A Comprehensive study of guanine excited state relaxation and photoreactivity in G-quadruplexes. *J. Phys. Chem. Lett.* **2019**, *10*, 6873–6877. [CrossRef]

54. Chatgilialoglu, C. The two faces of the guanyl radical: Molecular context and behavior. *Molecules* **2021**, *26*, 3511. [CrossRef]

55. Von Sonntag, C. *Free-Radical-Induced DNA Damage and Its Repair*; Springer: Berlin/Heidelberg, Germany, 2006; p. 523.

56. Balanikas, E.; Banyasz, A.; Baldacchino, G.; Markovitsi, D. Guanine radicals generated in telomeric G-quadruplexes by direct absorption of low-energy UV photons: Effect of potassium ions. *Molecules* **2020**, *25*, 2094. [CrossRef] [PubMed]

57. Ma, J.; Kumar, A.; Muroya, Y.; Yamashita, S.; Sakurai, T.; Denisov, S.A.; Sevilla, M.D.; Adhikary, A.; Seki, S.; Mostafavi, M. Observation of dissociative quasi-free electron attachment to nucleoside via excited anion radical in solution. *Nat. Commun.* **2019**, *10*, 102. [CrossRef] [PubMed]

58. Cadet, J.; Douki, T.; Ravanat, J.L. Oxidatively generated damage to the guanine moiety of DNA: Mechanistic aspects and formation in cells. *Acc. Chem. Res.* **2008**, *41*, 1075–1083. [CrossRef] [PubMed]

59. Fleming, A.M.; Burrows, C.J. G-quadruplex folds of the human telomere sequence alter the site reactivity and reaction pathway of guanine oxidation compared to duplex DNA. *Chem. Res. Toxicol.* **2013**, *26*, 593–607. [CrossRef]

60. Bielskute, S.; Plavec, J.; Podbevsek, P. Impact of oxidative lesions on the human telomeric G-quadruplex. *J. Am. Chem. Soc.* **2019**, *141*, 2594–2603. [CrossRef]

61. Dapprich, S.; Komaromi, I.; Byun, K.S.; Morokuma, K.; Frisch, M.J. A new ONIOM implementation in Gaussian98. Part, I. The calculation of energies, gradients, vibrational frequencies and electric field derivatives. *J. Mol. Struct.-Theochem* **1999**, *461*, 1–21. [CrossRef]

62. Zhao, Y.; Schultz, N.E.; Truhlar, D.G. Design of density functionals by combining the method of constraint satisfaction with parametrization for thermochemistry, thermochemical kinetics, and noncovalent interactions. *J. Chem. Theory Comput.* **2006**, *2*, 364–382. [CrossRef]

63. Zhao, Y.; Truhlar, D.G. Density functionals with broad applicability in chemistry. *Acc. Chem. Res.* **2008**, *41*, 157–167. [CrossRef] [PubMed]

64. Cornell, W.D.; Cieplak, P.; Bayly, C.I.; Gould, I.R.; Merz, K.M.; Ferguson, D.M.; Spellmeyer, D.C.; Fox, T.; Caldwell, J.W.; Kollman, P.A. A 2nd generation force-field for the simulation of proteins, nucleic acids and organic molecules. *J. Am. Chem. Soc.* **1995**, *117*, 5179–5197. [CrossRef]
65. Miertus, S.; Scrocco, E.; Tomasi, J. Electrostatic interaction of a solute with a continuum-A direct utilization of abinitio molecular potentials for the prevision of solvent effects. *Chem. Phys.* **1981**, *55*, 117–129.

66. Tomasi, J.; Mennucci, B.; Cammi, R. Quantum mechanical continuum solvation models. *Chem. Rev.* **2005**, *105*, 2999–3093.
[CrossRef]

67. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. *Gaussian 16 Rev. C.01*; Gaussian, Inc.: Wallingford, CT, USA, 2016.

68. Borrego-Varillas, R.; Cerullo, G.; Markovitsi, D. Exciton trapping dynamics in DNA multimers. *J. Phys. Chem. Lett.* **2019**, *10*, 1639–1643.
[CrossRef]

69. Ma, J.; Wang, F.; Denisov, S.A.; Adhikary, A.; Mostafavi, M. Reactivity of prehydrated electrons toward nucleobases and nucleotides in aqueous solution. *Sci. Adv.* **2017**, *3*, e1701669.
[CrossRef]