Sequence dependence of electron-induced DNA strand breakage revealed by DNA nanoarrays

The electronic structure of DNA is determined by its nucleotide sequence, which is for instance exploited in molecular electronics. Here we demonstrate that also the DNA strand breakage induced by low-energy electrons (18 eV) depends on the nucleotide sequence. To determine the absolute cross sections for electron induced single strand breaks in specific 13 mer oligonucleotides we used atomic force microscopy analysis of DNA origami based DNA nanoarrays. We investigated the DNA sequences 5′-TT(XYX)3TT with X = A, G, C and Y = T, BrU 5-bromouracil and found absolute strand break cross sections between 2.66 · 10−14 cm2 and 7.06 · 10−14 cm2. The highest cross section was found for 5′-TT(ATA)3TT and 5′-TT(ABrUA)3TT, respectively. BrU is a radiosensitizer, which was discussed to be used in cancer radiation therapy. The replacement of T by BrU into the investigated DNA sequences leads to a slight increase of the absolute strand break cross sections resulting in sequence-dependent enhancement factors between 1.14 and 1.66. Nevertheless, the variation of strand break cross sections due to the specific nucleotide sequence is considerably higher. Thus, the present results suggest the development of targeted radiosensitizers for cancer radiation therapy.

D NA exhibits nucleotide sequence dependent electronic properties, which affect its charge transport properties, UV stability, and sensitivity toward radiation [1][2][3] . These properties manifest themselves in quantities such as the ionization potential 4,5 , in the dynamics of the electronic states (excitonic coupling, lifetime of excited states) 1,6 , but also in the reactivity for instance towards low-energy electrons 3,7,8 .
A large number of secondary electrons is formed along the radiation track of high-energy radiation, which is routinely applied in radiation therapy to kill tumor tissue. These low-energy electrons (LEEs) have a most probable energy around 10 eV 9 and are able to directly induce DNA single and double strand breaks (SSBs and DSBs) via dissociative electron attachment through the formation of negative ion resonances [10][11][12] . Although it is experimentally extremely challenging to quantify the LEE induced DNA strand break (SB) yield of oligonucleotides of specific nucleotide sequence, a number of studies suggested a sequence dependence of electroninduced DNA strand breakage. In electron transmission through self-assembled DNA monolayers it was demonstrated that the number of electrons trapped in the DNA film depends strongly on the guanine (G) content 3 . These experiments did not yield any information on the damage of the DNA film, but the SSB yield upon irradiation with 1 eV electrons was later quantified using microarrays of DNA SAMs and fluorescence detection of the hybridisation efficiency 7 . It was found that the SSB yield increases linearly with the number of G bases present in the oligonucleotide. This behaviour is explained by the relative instability of G compared to the other DNA bases 13 , which is also reflected in the low ionization potential (IP) of G 14 . However, in addition to the basedependent damage, it was also demonstrated that the specific base sequence has a strong effect on the DNA damage. In electron transmission experiments the human telomeric repeat TTAGGG turned out to be particularly prone to LEE capture 15 . The particular role of the telomere sequence was confirmed by ab initio calculations showing that the IP of a TTAGGG sequence is smaller than the IP of a TTGGGG sequence, although isolated G has a lower IP than A 16 . Despite the particular role of G it was recently demonstrated by using HPLC analysis of oligonucleotide trimers irradiated with 10 eV electrons, that the total DNA damage was largest for TTT and smallest for TGT 8 . Nevertheless, the ratio of strand breaks (C-O cleavage) to nucleobase loss (C-N cleavage) was found to be highest for TGT 8 .
The investigation of LEE induced strand breakage of specific oligonucleotides longer than a few nucleobases is very challenging due to the small penetration depth of LEEs resulting in small amounts of damaged material 8 . Recently, we demonstrated the detection of LEEinduced bond cleavage in DNA origami based DNA nanoarrays on the single molecule level using atomic force microscopy (AFM) [17][18][19] . The visualization of DNA strand breakage by AFM using DNA nanoarrays has several advantages: (i) Due to the detection of DNA strand breaks at a single-molecule level only miniscule amounts of material are required to establish sub-monolayer surface coverage. (ii) Two or more different oligonucleotide sequences can be irradiated within a single experiment to efficiently compare a number of different DNA structures. (iii) Absolute strand break cross sections (s SSB ) are readily accessible (vide infra) thus providing benchmark values for further experimental and theoretical studies. (iv) The DNA nanoarray technique is not limited to single strands, but can be extended to quantify double strand breaks and to investigate higher-order DNA structures. Here, we compare the absolute strand break cross sections of different 13 mer oligonucleotide sequences. We compare the sequences 59-TT(XTX) 3 TT with X 5 A, C or G to evaluate the role of the different DNA nucleobases for DNA strand breakage. In a next step, we study the sensitizing effect of the incorporation of 5-bromouracil (BrU) by comparing the absolute strand break cross sections of the sequences 59-TT(XBrUX) 3 TT with X 5 A, C or G. Fig. 1a shows a scheme of a triangular DNA origami platform carrying six target sequences, which represents the DNA nanoarray. Three biotinylated target sequences are situated in the center of the trapezoids, and three are located on the sides of the trapezoids. Basically, the nucleotide sequence of each of the six target oligonucleotides can be freely chosen. In the present study the three central target strands (green in Fig. 1a), and the three side positions (black in Fig. 1a) cannot be distinguished in the AFM images and are hence chosen to have the same sequence. Thus, two different target sequences are studied within one irradiation experiment, and the specific sequences and their positions are indicated in Fig. 1a. In Fig. 1b typical AFM images from DNA origami samples after incubation with streptavidin (SAv) which binds to biotin (Bt) are shown. The left image shows a control sample that was not irradiated while the image on the right was obtained from a sample irradiated with 18 eV electrons at a fluence of 5.0 3 10 12 cm 22 . The energy of 18 eV was chosen due to its relevance for the damage induced by secondary electrons originating from the ionization track of high-energy radiation. For secondary electrons the damage probability has a global maximum around 18 eV, i.e. the damage induced mainly by ionization and electronic excitation weighted by the LEE distribution in aqueous samples irradiated with high-energy radiation 20 . In the AFM image on the right-hand side of Fig. 1 the number of specifically bound SAv is reduced compared to the non-irradiated control sample indicating that a number of target sequences have been damaged by electron-induced strand breakage. To determine the absolute cross section for strand breakage (see Methods section), the relative number of SBs (N SB ) was recorded as a function of the electron fluence. The fluence dependence of N SB for the target sequences TT(XTX) 3 TT with X 5 A, C, G is displayed in Fig. 2a. From the linear fit in the lowfluence regime s SSB is determined, which is shown in Fig. 2b. The oligonucleotide TT(GTG) 3 TT shows the lowest response to 18 eV electrons. To ensure for an accurate linear fit, a smaller fluence increment and thus more data points were chosen for TT(GTG) 3 TT.

Results and Discussion
In general, there are three basic mechanisms that could account for the electron-induced DNA strand cleavage at 18 eV (in the following the general form TT(XYX) 3 3 TT and TT(GBrUG) 3 TT oligonucleotides is shown. On the right, a typical AFM image after irradiation with 18 eV electrons is shown. The number of specifically bound SAv is reduced due to strand breaks in the protruding sequences.
www.nature.com/scientificreports SCIENTIFIC REPORTS | 4 : 7391 | DOI: 10.1038/srep07391 ment at 18 eV is very unlikely, but the initial incoming electron might undergo inelastic scattering in the surrounding of the target sequence and might attach at lower energies with higher cross sections. DEA cross sections for instance for the loss of hydrogen from the thymine anion at 1.0 eV have been previously determined to be 7.8 ? 10 217 cm 2,23 . Alternatively, low energy secondary electrons can be generated either from the substrate, or by electron impact ionization according to (i), which can induce a SB by DEA. iii. Neutral dissociation: TT(XYX) 3 TT 1 e 2 (18 eV) R TT(XYX) 3 TT** 1 e 2 (,18 eV) R SB. The transfer of electronic energy from the incoming electron to the oligonucleotide might be associated with a dissociative state resulting in a SB. The inelastic electron scattering cross section is comparable with the cross section for electron impact ionization and is e.g. for pyrimidine of the order of 10 216 cm 2,24 . Bond breaking by a catalytic electron is a similar mechanism that involves transient negative ions and neutral fragmentation products 25 .
Judging from the magnitude of cross sections it is most likely that the initial step in strand breakage is either ionization or electronic excitation of the target strands. For a DNA strand break to occur a bond within the phosphate-sugar backbone needs to be cleaved. A direct SB without the involvement of the nucleobases as it was suggested previously is thus feasible 26,27 . The IP of the phosphate-sugar backbone is with about 11.7 eV considerably higher than the IPs of the nucleobases, nevertheless electron impact ionization cross sections of the DNA backbone are calculated to be higher than that of the nucleobases due to the higher number of electrons in the sugarphosphate backbone 22 .
However, the absolute strand break cross sections determined here (Fig. 2) show a strong dependence on the specific nucleotide sequence with the TT(GTG) 3 TT sequence having the lowest s SSB and TT(ATA) 3 TT having the highest s SSB . This indicates that the nucleobases have a strong influence on the strand breakage and a damage mechanism involving only the DNA backbone is not likely.
The absolute cross sections for strand breakage vary from (2.21 6 0.87) ? 10 214 cm 2 to (6.00 6 0.86) ? 10 214 cm 2 depending on the specific nucleotide sequence (see table 1 for details). The cross sections are comparable with the cross sections found for electron induced strand breakage in plasmid DNA. Boudaiffa et al. found an effective cross section for SSBs at 10 eV electron energy of 2.6 ? 10 215 cm 2 using the plasmid pGEM 3Zf(2) with 3199 base pairs 28 . Later, Panajotovic et al. reported effective SSB cross sections of 10.8 ? 10 215 cm 2 (at 10 eV) and 24.8 ? 10 215 cm 2 (at 1 eV) using the same plasmid 29 . Very recently the absolute cross section for loss of supercoiled plasmid DNA was determined at 10 eV to be 3.0 ? 10 214 cm 2,30 , which is in accordance with the values found in the present study. In the previous studies agarose gel electrophoresis has been used to detect SSBs in plasmid DNA. This method and the present DNA nanoarray based method both detect the first SSB that occurs upon electron irradiation. Subsequent SSBs that occur due to the impact of additional electrons are no more detected. Since the geometrical cross sections of DNA with more than a few nucleotides (1 nm 2 5 10 214 cm 2 corresponds to about three nucleotides) are larger than the SSB cross sections, the reported cross sections should be independent of the size of the system (plasmid DNA vs. oligonucleotides).
The trend of s SSB with respect to the nucleotide sequence is surprising since G containing sequences are generally assumed to be particularly fragile due to the small ionization potential of G. Both, the vertical and adiabatic IPs of the isolated nucleobases increase in the order of G , A , C , T 14 . A similar order is found for the total electron impact ionization cross sections at low energy (20 eV): G < A . T < C 21 . The IP can be considered an important quantity in this context since at an energy of 18 eV the electron ionization cross section is presumably orders of magnitude higher than the DEA cross section (vide supra). However, the absolute strand break cross sections found in the present study for the different sequences cannot be explained with the electron impact ionization properties of the individual nucleobases.
The adjacent nucleobases interact via stacking interactions between the aromatic systems, and it was previously found that the strength of p-p interaction depends strongly on the type of nucleobases. Stacking interactions are strongest between the purine bases A and G, and in the case of A nucleobases the stacking interactions further increase by a propeller twist of the A bases. This effect plays also a role for AT base stacks, but is not relevant for GG and GT interactions 31 . In DFT calculations it was also found that the IP of G stacks decreases more with the number of bases than in the case of A stacks 5 . Nevertheless it was found that the difference between vertical and adiabatic IP, i.e. the reorganization energy, is smaller for A stacks. A smaller reorganization energy facilitates hole delocalization. On the other hand the large nuclear reorganization energy  results in hole localization in G stacks 5 . This is consistent with the observation that G stacks act as hole traps in DNA.
In order to find out whether the nucleotide sequences studied here give rise to specific modifications of their electronic structure we have computed the IPs of stacked oligonucleotide trimers. The IPs obtained with two different theoretical methods (MP2 and GW) are presented in Table 2. According to the results obtained by MP2 calculations the ATA stack has an IP of 8.01 eV, which is indeed lower than the IP of the GTG stack (8.19 eV) and thus does not follow the same trend as the IPs of the isolated nucleobases 14 . The vertical IP of isolated A is about 0.4 eV higher than that of isolated G 14 . Nevertheless the IP of the ATA stack is 0.18 eV lower than the IP of GTG according to the MP2 calculations, which is a considerable shift compared to what is expected from the isolated nucleobases. According to the GW method, however, the trend of IPs reflects the trend observed for the isolated nucleobases, i.e. the GTG stack has the lowest IP (7.96 eV). This difference between MP2 and GW method might be explained when considering also the excited states of the system by means of CASPT2. A detailed comparison of the MP2 and GW methods is presented in the supporting information. As a consequence, a correlation of the observed absolute strand break cross sections and the computed IPs for the stacked trimers is not straightforward. Here it should be noted that in the calculations only a single conformation of the stacked trimers is considered and environmental effects are not included. Furthermore, the experimental SB cross sections are determined for oligonucleotide 13 mers, but due to limited computing resources we could only consider the trimers listed in table 2. Consequently, further detailed investigations have to be performed on complex systems to elucidate possible connections between the sequence-dependence of strand breakage and the respective IPs.
A particular role of the electronic states of A containing oligonucleotides was found already in previous studies. It was demonstrated that due to electronic coupling between A molecules in an A homopolymer excitons in the VUV (190 nm) spectral region are extended over up to eight nucleobases 6 . However, the electronic coupling can be eliminated by a single T spacer 32 . On the other hand, by using HPLC analysis of oligonucleotide trimers irradiated with 10 eV electrons the damage due to strand breaks was found to increase in the row of sequences TCT < TGT , TAT, which is consistent with the trend of s SSB values determined here 33 .
5-Bromouracil (BrU) is a well-known radiosensitizer, whose incorporation into DNA by replacing T increases the SSB and DSB yields upon irradiation with electrons and photons [34][35][36] . We have studied the effect of BrU incorporation on the strand break cross sections by using the sequences discussed above and replacing the three central T bases by BrU. The absolute strand break cross sections determined for 18 eV electron irradiation are shown in Fig. 3a. Again, the TT(ABrUA) 3 TT sequence exhibits the highest s SSB , whereas the TT(CBrUC) 3 TT sequence has the lowest s SSB . There is only a small difference between TT(CBrUC) 3 TT and TT(GBrUG) 3 TT, but the cross section of TT(ABrUA) 3 TT is about twice as high as that of TT(CBrUC) 3 TT. The exact values of s SSB of the different sequences are summarized in Table 1.
The similarity of s SSB for oligonucleotides with and without BrU indicates that s SSB is much more sensitive to the nucleotide sequence than to the presence of BrU. The influence of BrU on s SSB might be higher at lower electron energies, since DEA to gas phase BrU shows the highest cross sections for Br 2 formation close to zero eV 37 .
In Table 2 the computed IPs of the BrU containing stacked trimers are listed. With both methods of calculation, MP2 and GW, the same trend of IPs is observed, i.e. the IP of GBrUG is the lowest and CBrUC is the highest. There are basically no differences between the IPs of stacked trimers with and without BrU. The only exception is the IP of AYA calculated by MP2, which is 8.01 eV for Y 5 T, and 8.52 eV for Y 5 BrU. Fig. 3b shows the enhancement factors for BrU incorporation (EF 5 s SSB (XTX)/s SSB (XBrUX)), which are found to be the highest for TT(GYG) 3 TT (1.66) and the lowest for TT(CYC) 3 TT (1.14), see Table 1 for details. Thus, the effect of BrU incorporation is highest when BrU is directly adjacent to G. In a recent study using femtosecond laser spectroscopy, the observation of anionic transients suggested that the effect of BrU should be the strongest in close proximity to A. A similar effect but slightly weaker was inferred to G, since G was assumed to be the major damaging site 38 . However, in this particular study, only mixtures of BrdU and dAMP/dGMP were investigated instead of oligonucleotides so that no information about the occurrence of strand breaks could be obtained. In contrast, we directly probe s SSB of different oligonucleotide sequences in our experiments.
In a recent study using HPLC analysis of trimers irradiated with 10 eV electrons it was found that the total damage of TBrUT was about 50% higher than that of TTT. However, most of the damage   was associated with formation of TUT 39 . In another study using the same technique it was found that the amount of fragments associated with strand breaks is approximately the same for TBrUT and TUT, i.e. no increase of strand breaks was observed with TBrUT 33 . Therefore, when comparing our results with previous studies, the different sequences, lengths of oligonucleotides and the electron energy must be taken into account. The present study represents a starting point for a global assessment of the sequence dependence of LEE induced DNA damage, for which a large number of DNA sequences at a range of electron energies must be studied in the future. The DNA nanoarray technique is suitable for that since different oligonucleotides can be probed in a single irradiation experiment and absolute strand break cross sections serving as benchmark values are obtained.

Conclusions
We have determined the absolute cross sections for DNA single strand breakage induced by 18 eV electrons using DNA origami based DNA nanoarrays. Since the analysis of irradiated samples is performed by atomic force microscopy at a single molecule level this novel method represents a simple way to access absolute strand break cross sections. The absolute single strand break cross sections depend strongly on the nucleotide sequence and we find values varying between (2.21 6 0.87) ? 10 214 cm 2 and (6.00 6 0.86) ? 10 214 cm 2 for the oligonucleotide sequences TT(XTX) 3 TT with X 5 A, C, G. Furthermore, we find that exchange of the central T bases by 5bromouracil increases the strand break cross section in a sequencedependent manner by a factor of 1.14 to 1.66. The observed trend in the absolute strand break cross sections agrees qualitatively with previous HPLC studies investigating the fragmentation of oligonucleotide trimers of the sequence TXT with X 5 A, C, G irradiated with 10 eV electrons 33 . As in the here presented results, the authors observed almost identical yields of strand break fragments for the C-and G-containing trimers while the yield for the A-containing trimer was about twice as high. In addition, the absolute strand break cross sections measured here are comparable in magnitude with cross sections for strand breakage in different plasmid DNAs induced by 1-10 eV electrons as determined by agarose gel electrophoresis 29,30 . The DNA nanoarray technique thus bridges the gap between genomic dsDNA several kbp in size and very short oligonucleotides only few nt long, and enables the detailed investigation of sequence-dependent processes in DNA radiation damage. The observed sequence-specificity most likely results from the modification of electronic states by electronic coupling of the individual nucleobases. The sensitivity of the A containing sequences might be associated with the strong stacking interactions and pronounced hole delocalization in adjacent A bases. Further experimental and theoretical studies will be carried out covering a broad range of electron energies and DNA sequences to elucidate the most relevant damage mechanisms.
The present results suggest that radiosensitizers applied in tumor radiation therapy could operate more efficiently if they targeted specific nucleotide sequences that have the highest damage cross sections. Further extended experiments have to be performed to explore the electron energy dependence of radiosensitization and thus the physico-chemical mode of action of established and potential radiosensitizers 40,41 .

Methods
Preparation of DNA nanoarrays. Triangular DNA origami nanostructures have been prepared from the circular single-stranded viral DNA scaffold strand M13mp18 and 208 short artificial staple strands according to the original design and procedure by Rothemund 42 . Selected staple strands have been extended with a specific target sequence and a Bt modification at the 59-end. The extended staple strands form a nanoscale array, and the individual target sequences can be visualized with AFM after incubation with SAv, which binds to the Bt modifications of the intact protruding strands. For each irradiation experiment two different target sequences have been selected and in total 6 staple strands per DNA origami structure have been modified (see Fig. 1).
The DNA origami structures are assembled by annealing of the scaffold strand with a 30-fold excess of staple stands from 65uC to 4u within approximately 2 hours (1 3 TAE buffer, 10 mM MgCl 2 ). The excess staple strands are separated from the assembled DNA origami structures by spinfiltering twice using Amicon Ultra centrifuge filters (100.000 Da MWCO).
Electron irradiation and AFM analysis. The detailed procedure of immobilization and LEE irradiation of the DNA nanoarrays is given in ref. 17. In brief, the DNA origami structures are bound electrostatically to Si/SiO 2 by incubation in 10 3 TAE and 100 mM MgCl 2 for approximately one hour. Afterwards, the dry samples are transferred into ultrahigh vacuum (UHV), and irradiated with LEEs of defined fluence at a current of 1-10 nA. After irradiation the samples are removed from the UHV chamber and rinsed to remove fragmentation products. To identify the intact remaining target sequences the samples are exposed to a 50 nM solution of SAv for about 2-10 minutes. Then the solution is rinsed again and the dry samples are analysed by AFM. From the AFM images the relative number of DNA strand breaks of a given sequence N SB can be determined from the number of intact oligonucleotides after electron irradiation compared to the initial number of oligonucleotides prior to irradiation (i.e. three oligonucleotides per DNA origami triangle). From N SB the sequence-specific absolute cross section for strand breakage can be determined.
The cross section for DNA damage can be described by the following exposureresponse relation 29 : with N 0 being the initial number of DNA oligonucleotides protruding from the DNA origami platform, s the cross section for DNA damage, J the electron flux and t the irradiation time. The number of strand breaks N SB (i.e. the relative number of damaged oligonucleotides: 1 -N(t)/N 0 ) can be approximated for short irradiation times by using a Taylor series: with F 5 Jt being the electron fluence. Thus, the cross section s can be determined from the slope of N SB (F) in the low-fluence regime. The determined cross section can be conceived as an absolute cross section since it is based on single molecule measurements, i.e. it is not effective for a specific DNA density or film thickness. The obtained absolute strand break cross sections are corrected for electron-induced damage to the biotin label, which was previously determined to be (1.1 6 0.2) ? 10 214 cm 2 18 . To support the presented data the absolute strand break cross section for brominated sequences was determined also in a second, different design of DNA nanoarrays (see supporting information), and the obtained s SSB values agree with the ones shown here.
Calculation of IPs. The ground-state geometries of all the investigated structures (i.e. XTX and XBrUX with X 5 C, A and G) were first optimized in the gas phase at the density functional theory (DFT) level of theory using the B97-D functional 43 and a 6-311G(d,p) polarized basis set. Frequency calculations for each derivative have been performed with the same level of theory, verifying that all structures correspond to true minima of the potential energy surface. IPs were calculated using spin-restricted open-shell second-order Møller-Plesset perturbation theory (MP2) 44 and a 6-31111G(d,p) basis set. Vertical IPs were obtained for all structures from the difference in total energy between the neutral species and the radical cation and anion, respectively, evaluated at the optimized geometry of the neutral species. All MP2 calculations have been achieved with the Gaussian 09 program suite. In addition, we employed many-body perturbation theory calculations in the GW approximation 45,46 . GW is considered as one of the most accurate electronic structure method for the calculation of charged excitations (such as IPs) feasible for the systems studied in this work. All GW calculations were performed in FHIaims 47 using the consistent starting point scheme 48 and a converged basis set (tier 4). More details on the calculations are provided in the SI.