Locations of radiation-produced DNA double strand breaks along chromosomes: a stochastic cluster process formalism

Nuclear dynamics of radiation-induced foci in euchromatin and heterochromatin

Repair of double strand breaks (DSBs) is essential for cell survival and genome integrity. While much is known about the molecular mechanisms involved in DSB repair and checkpoint activation, the roles of nuclear dynamics of radiation-induced foci (RIF) in DNA repair are just beginning to emerge. Here, we summarize results from recent studies that point to distinct features of these dynamics in two different chromatin environments: heterochromatin and euchromatin. We also discuss how nuclear architecture and chromatin components might control these dynamics, and the need of novel quantification methods for a better description and interpretation of these phenomena. These studies are expected to provide new biomarkers for radiation risk and new strategies for cancer detection and treatment.

Integration of Monte Carlo simulations with PFGE experimental data yields constant RBE of 2.3 for DNA double-strand break induction by nitrogen ions between 125 and 225 keV/μm LET

The number of small radiation-induced DNA fragments can be heavily underestimated when determined from measurements of DNA mass fractions by gel electrophoresis, leading to a consequent underestimation of the initial DNA damage induction. In this study we reanalyzed the experimental results for DNA fragmentation and DNA double-strand break (DSB) yields in human fibroblasts irradiated with γ rays and nitrogen ion beams with linear energy transfer (LET) equal to 80, 125, 175 and 225 keV/μm, originally measured by Höglund et al. (Radiat Res 155, 818-825, 2001 and Int J Radiat Biol 76, 539-547, 2000). In that study the authors converted the measured distributions of fragment masses into DNA fragment distributions using mid-range values of the measured fragment length intervals, in particular they assumed fragments with lengths in the interval of 0-48 kbp had the mid-range value of 24 kbp. However, our recent detailed simulations with the Monte Carlo code PARTRAC, while reasonably in agreement with the mass distributions, indicate significantly increased yields of very short fragments by high-LET radiation, so that the actual average fragment lengths, in the interval 0-48 kbp, 2.4 kbp for 225 keV/μm nitrogen ions were much shorter than the assumed mid-range value of 24 kbp. When the measured distributions of fragment masses are converted into fragment distributions using the average fragment lengths calculated by PARTRAC, significantly higher yields of DSB related to short fragments were obtained and resulted in a constant relative biological effectiveness (RBE) for DSB induction yield of 2.3 for nitrogen ions at 125-225 keV/μm LET. The previously reported downward trend of the RBE values over this LET range for DSB induction appears to be an artifact of an inadequate average fragment length in the smallest interval.

Mathematical models of the generation of radiation-induced DNA double-strand breaks

The double-strand break (dsb) is one of the most critical lesions leading to a variety of radiobiological effects. In this paper, we reconsider the previously constructed and generally accepted mathematical models for dsb generation, and give a concrete mathematical basis for the generation of dsbs and the calculation of the number of induced dsbs, under the assumption of randomness in the break location in DNA and in the number of breaks. Using these models based on the Poisson distribution and the binomial distribution, we calculate the dose dependence of dsb generation. We deduced from our models that the dose dependence of the number of dsbs is described approximately as a quadratic form in both distribution models where dsb generation is accounted for by two ssbs. Previously reported experimental data on the dsb generation in phage DNA was found to be in good agreement with our models. Though the widely used model, the linear quadratic (LQ) model or the molecular theory of dsb formation based on the Poisson distribution, also gives the quadratic term, in spite of rough estimates or some mathematical incompleteness, a marked feature of our formulation is the absence of a parameter like the [Formula: see text] in the quadratic term that requires experimental data to determine. Thus in this study we provide mathematical validity to the generally accepted models of the number of dsb.

Evidence for formation of DNA repair centers and dose-response nonlinearity in human cells

The concept of DNA "repair centers" and the meaning of radiation-induced foci (RIF) in human cells have remained controversial. RIFs are characterized by the local recruitment of DNA damage sensing proteins such as p53 binding protein (53BP1). Here, we provide strong evidence for the existence of repair centers. We used live imaging and mathematical fitting of RIF kinetics to show that RIF induction rate increases with increasing radiation dose, whereas the rate at which RIFs disappear decreases. We show that multiple DNA double-strand breaks (DSBs) 1 to 2 μm apart can rapidly cluster into repair centers. Correcting mathematically for the dose dependence of induction/resolution rates, we observe an absolute RIF yield that is surprisingly much smaller at higher doses: 15 RIF/Gy after 2 Gy exposure compared to approximately 64 RIF/Gy after 0.1 Gy. Cumulative RIF counts from time lapse of 53BP1-GFP in human breast cells confirmed these results. The standard model currently in use applies a linear scale, extrapolating cancer risk from high doses to low doses of ionizing radiation. However, our discovery of DSB clustering over such large distances casts considerable doubts on the general assumption that risk to ionizing radiation is proportional to dose, and instead provides a mechanism that could more accurately address risk dose dependency of ionizing radiation.

Track structures, DNA targets and radiation effects in the biophysical Monte Carlo simulation code PARTRAC

This review describes the PARTRAC suite of comprehensive Monte Carlo simulation tools for calculations of track structures of a variety of ionizing radiation qualities and their biological effects. A multi-scale target model characterizes essential structures of the whole genomic DNA within human fibroblasts and lymphocytes in atomic resolution. Calculation methods and essential results are recapitulated regarding the physical, physico-chemical and chemical stage of track structure development of radiation damage induction. Recent model extension towards DNA repair processes extends the time dimension by about 12 orders of magnitude and paves the way for superior predictions of radiation risks.

Separating DNA with different topologies by atomic force microscopy in comparison with gel electrophoresis

Atomic force microscopy, which is normally used for DNA imaging to gain qualitative results, can also be used for quantitative DNA research, at a single-molecular level. Here, we evaluate the performance of AFM imaging specifically for quantifying supercoiled and relaxed plasmid DNA fractions within a mixture, and compare the results with the bulk material analysis method, gel electrophoresis. The advantages and shortcomings of both methods are discussed in detail. Gel electrophoresis is a quick and well-established quantification method. However, it requires a large amount of DNA, and needs to be carefully calibrated for even slightly different experimental conditions for accurate quantification. AFM imaging is accurate, in that single DNA molecules in different conformations can be seen and counted. When used carefully with necessary correction, both methods provide consistent results. Thus, AFM imaging can be used for DNA quantification, as an alternative to gel electrophoresis.

DNA fragmentation in V79 cells irradiated with light ions as measured by pulsed‐field gel electrophoresis. II. Simulation with a generalized broken stick model

PURPOSE

To characterize the differences among the experimental DNA fragmentation spectra induced in Chinese hamster V79 cells by gamma-rays, low-energy protons and alpha-particles through the use of a phenomenological model.

MATERIALS AND METHODS

A model of DNA fragmentation was developed as a generalization of the broken-stick model, in which the double-strand breaks induced by radiation were considered randomly placed, but in which the manifestly non-random fragmentation of the control sample was fully taken into account and considered as the initial fragment distribution. Further, an analytical method was introduced that allowed an evaluation of the deviation from randomness of the fragmentation induced by radiation.

RESULTS

The analysis of the experimental distribution of DNA fragments showed that there was a progressive departure from randomness in radiation-induced fragmentation going from gamma-rays to protons and then to alpha-particles. This deviation was characterized by an enhanced induction of fragments, and therefore by a larger correlation of double-strand breaks, in the experimental range of lower molecular weights.

CONCLUSION

The analysis shows that low-energy light ions induce DNA fragmentation, at the loop level of the chromatin organization, that can be significantly non-random. The same analysis can readily be applied at different length scales, and thus it could offer a basis for the study of the link between DNA damage, correlated at various spatial scales and biological end-points.

UVA generates pyrimidine dimers in DNA directly

There is increasing evidence that UVA radiation, which makes up approximately 95% of the solar UV light reaching the Earth's surface and is also commonly used for cosmetic purposes, is genotoxic. However, in contrast to UVC and UVB, the mechanisms by which UVA produces various DNA lesions are still unclear. In addition, the relative amounts of various types of UVA lesions and their mutagenic significance are also a subject of debate. Here, we exploit atomic force microscopy (AFM) imaging of individual DNA molecules, alone and in complexes with a suite of DNA repair enzymes and antibodies, to directly quantify UVA damage and reexamine its basic mechanisms at a single-molecule level. By combining the activity of endonuclease IV and T4 endonuclease V on highly purified and UVA-irradiated pUC18 plasmids, we show by direct AFM imaging that UVA produces a significant amount of abasic sites and cyclobutane pyrimidine dimers (CPDs). However, we find that only approximately 60% of the T4 endonuclease V-sensitive sites, which are commonly counted as CPDs, are true CPDs; the other 40% are abasic sites. Most importantly, our results obtained by AFM imaging of highly purified native and synthetic DNA using T4 endonuclease V, photolyase, and anti-CPD antibodies strongly suggest that CPDs are produced by UVA directly. Thus, our observations contradict the predominant view that as-yet-unidentified photosensitizers are required to transfer the energy of UVA to DNA to produce CPDs. Our results may help to resolve the long-standing controversy about the origin of UVA-produced CPDs in DNA.

Cluster identification in AA5754 aluminium sheets using mathematical morphology analysis

Quantitative image analysis of particle distribution in the microstructure of continuous cast (CC) and direct chill cast (DC) AA5754 aluminium alloy sheets have been conducted. This information can be used as an input for modelling mechanical deformation and instability in these materials. The quantitative analysis reveals that there are significant differences in the microstructure of the two materials even though the total content of second-phase particles is statistically similar. Qualitative observation shows the second-phase particles to be arranged in the form of streaks parallel to the rolling direction in the CC sheets and in a uniform random manner in the DC sheets. The main difference in the geometric microstructure of the CC and DC material is the spatial arrangement of the second-phase particles. A new mathematical technique called proximity analysis is developed to identify clusters and group of particles belonging to a cluster. Quantification through proximity analysis reveals that the particle clusters in CC sheet are in the form of long clusters (streaks) parallel to the rolling direction and are significantly longer than those in DC sheets (with the largest cluster in CC being four times larger than DC), and also have anisotropic angular orientation parallel to the rolling direction. The lower value of fracture strain observed in the CC sheets compared to DC sheets is attributed to a combination of large sizes of clusters and their preferential alignment along the rolling direction in the CC microstructure.

Nanoscale detection of ionizing radiation damage to DNA by atomic force microscopy

The detection and quantification of ionizing radiation damage to DNA at a single-molecule level by atomic force microscopy (AFM) is reported. The DNA damage-detection technique combining supercoiled plasmid relaxation assay with AFM imaging is a direct and quantitative approach to detect gamma-ray-induced single- and double-strand breaks in DNA, and its accuracy and reliability are validated through a comparison with traditional agarose gel electrophoresis. In addition, the dependence of radiation-induced single-strand breaks on plasmid size and concentration at a single-molecule level in a low-dose (1 Gy) and low-concentration range (0.01 ng microL(-1)-10 ng microL(-1)) is investigated using the AFM-based damage-detection assay. The results clearly show that the number of single-strand breaks per DNA molecule is linearly proportional to the plasmid size and inversely correlated to the DNA concentration. This assay can also efficiently detect DNA damage in highly dilute samples (0.01 ng microL(-1)), which is beyond the capability of traditional techniques. AFM imaging can uniquely supplement traditional techniques for sensitive measurements of damage to DNA by ionizing radiation.

Image-based modeling reveals dynamic redistribution of DNA damage into nuclear sub-domains

Several proteins involved in the response to DNA double strand breaks (DSB) form microscopically visible nuclear domains, or foci, after exposure to ionizing radiation. Radiation-induced foci (RIF) are believed to be located where DNA damage occurs. To test this assumption, we analyzed the spatial distribution of 53BP1, phosphorylated ATM, and γH2AX RIF in cells irradiated with high linear energy transfer (LET) radiation and low LET. Since energy is randomly deposited along high-LET particle paths, RIF along these paths should also be randomly distributed. The probability to induce DSB can be derived from DNA fragment data measured experimentally by pulsed-field gel electrophoresis. We used this probability in Monte Carlo simulations to predict DSB locations in synthetic nuclei geometrically described by a complete set of human chromosomes, taking into account microscope optics from real experiments. As expected, simulations produced DNA-weighted random (Poisson) distributions. In contrast, the distributions of RIF obtained as early as 5 min after exposure to high LET (1 GeV/amu Fe) were non-random. This deviation from the expected DNA-weighted random pattern can be further characterized by “relative DNA image measurements.” This novel imaging approach shows that RIF were located preferentially at the interface between high and low DNA density regions, and were more frequent than predicted in regions with lower DNA density. The same preferential nuclear location was also measured for RIF induced by 1 Gy of low-LET radiation. This deviation from random behavior was evident only 5 min after irradiation for phosphorylated ATM RIF, while γH2AX and 53BP1 RIF showed pronounced deviations up to 30 min after exposure. These data suggest that DNA damage–induced foci are restricted to certain regions of the nucleus of human epithelial cells. It is possible that DNA lesions are collected in these nuclear sub-domains for more efficient repair.

DNA damages are daily cellular events. If such events are left unchecked in an organism, they can lead to DNA mutations and possibly cancer over a long period of time. Consequently, cells have very efficient DNA repair machinery. Many studies have focused on the different molecular factors involved in the repair machinery, neglecting to consider the spatial context where damage occurs. Therefore, little is known about the role the nuclear architecture might have in the DNA damage response. In this study, we introduce computer modeling and image analysis tools in order to relate the position of DNA damage markers to morphologically distinct regions of the nucleus. Using these tools, we show that radiation-induced damages locate preferentially in non-condensed DNA regions or at the boundary of regions with condensed DNA. These results contradict the current dogma that the molecular response to randomly generated DNA damages is independent of their nuclear locations. Instead, this suggests the existence of repair centers in the nucleus. Overall, our approach shows that nuclear architecture plays a role in the DNA damage response, reminding us that the nucleus is not simply a soup of DNA and proteins.

Detecting ultraviolet damage in single DNA molecules by atomic force microscopy

We report detection and quantification of ultraviolet (UV) damage in DNA at a single molecule level by atomic force microscopy (AFM). By combining the supercoiled plasmid relaxation assay with AFM imaging, we find that high doses of medium wave ultraviolet (UVB) and short wave ultraviolet (UVC) light not only produce cyclobutane pyrimidine dimers (CPDs) as reported but also cause significant DNA degradation. Specifically, 12.5 kJ/m(2) of UVC and 165 kJ/m(2) of UVB directly relax 95% and 78% of pUC18 supercoiled plasmids, respectively. We also use a novel combination of the supercoiled plasmid assay with T4 Endonuclease V treatment of irradiated plasmids and AFM imaging of their relaxation to detect damage caused by low UVB doses, which on average produced approximately 0.5 CPD per single plasmid. We find that at very low UVB doses, the relationship between the number of CPDs and UVB dose is almost linear, with 4.4 CPDs produced per Mbp per J/m(2) of UVB radiation. We verified these AFM results by agarose gel electrophoresis separation of UV-irradiated and T4 Endonuclease V treated plasmids. Our AFM and gel electrophoresis results are consistent with the previous result obtained using other traditional DNA damage detection methods. We also show that damage detection assay sensitivity increases with plasmid size. In addition, we used photolyase to mark the sites of UV lesions in supercoiled plasmids for detection and quantification by AFM, and these results were found to be consistent with the results obtained by the plasmid relaxation assay. Our results suggest that AFM can supplement traditional methods for high resolution measurements of UV damage to DNA.

A Robust Procedure for Removing Background Damage in Assays of Radiation-Induced DNA Fragment Distributions

The non-random distribution of DNA breakage in PFGE (pulsed-field gel electrophoresis) experiments poses a problem of proper subtraction of the background DNA damage to obtain a fragment-size distribution due to radiation only. A naive bin-to-bin subtraction of the background signal will not result in the right DNA mass distribution histogram. This problem could become more pronounced for high-LET (linear energy transfer) radiation, because the fragment-size distribution manifests a higher frequency of smaller fragments. Previous systematic subtraction methods have been based on random breakage, appropriate for low-LET radiation. Moreover, an investigation is needed to determine whether the background breakage is itself random or non-random. We consider two limiting cases: (1) the background damage is present in all cells, and (2) it is present in only a small subset of cells, while other cells are not contributing to the background DNA fragmentation. We give a generalized formalism based on stochastic processes for the subtraction of the background damage in PFGE experiments for any LET and apply it to two sets of PFGE data for iron ions.

Clusters of DNA Double-Strand Breaks Induced by Different Doses of Nitrogen Ions for Various LETs: Experimental Measurements and Theoretical Analyses

The yields and clustering of DNA double-strand breaks (DSBs) were investigated in normal human skin fibroblasts exposed to gamma rays or to a wide range of doses of nitrogen ions with various linear energy transfers (LETs). Data obtained by pulsed-field gel electrophoresis on the dose and LET dependence of DNA fragmentation were analyzed with the randomly located clusters (RLC) formalism. The formalism considers stochastic clustering of DSBs along a chromosome due to chromatin structure, particle track structure, and multitrack action. The relative biological effectiveness (RBE) for the total DSB yield did not depend strongly on LET, but particles with higher LET produced higher fractions of small DNA fragments, corresponding in the formalism to an increase in the average number of DSBs per DSB cluster. The results are consistent with the idea that DSB clustering along chromosomes is what leads to large RBEs of high-LET radiations for major biological end points. At a given dose, large fragments are less affected by the variability in LET than small fragments, suggesting that the two free ends in large fragments are often produced by two different tracks. The formalism successfully described an extra increase in small DNA fragments as dose increases and a related decrease in large fragments, mainly due to interlacing of DSB clusters produced along a chromosome by different tracks, since interlacing cuts larger DNA fragments into smaller ones.

Modelling the probability distribution of the number of DNA double-strand breaks due to sporadic alkylation of nucleotide bases

Metabolites and certain chemical agents (for example methyl methanesulfonate) can induce nucleotide bases on chromosomal strands to become alkylated. These alkylated sites have the potential to become single-strand chromosomal breaks, a form of DNA damage, if they are exposed to a sufficient temperature in vitro. It has been proposed that a single-strand break (SSB) sufficiently close to another SSB on the opposite chromosomal strand will form a double-strand break (DSB). DNA repair mechanisms are less able to repair DSBs compared to SSBs. Because of the complex three-dimensional structure of DNA, some chromosomal regions are more susceptible to alkylation than others. A question of interest is therefore whether these alkylated bases are randomly distributed or tend to be clustered. Pulsed-field gel electrophoresis allows the number of DNA fragments (and hence the number of DSBs) to be observed directly. The randomness of alkylation events can therefore be tested using the standard statistical hypothesis-testing framework. Under the null hypothesis, that the SSBs are randomly distributed on each of the strands, we can calculate the probability of observing a number of DSBs at least as large as that observed and hence the associated p-value. Previously, the probability distribution of the number of DSBs has been determined by Monte Carlo simulations; when considering the whole genome this can be very time consuming. In this paper, we theoretically derive an approximation to the distribution enabling appropriate probabilities to be calculated quickly. Based on previous findings we assume that the number of breaks on each strand is small compared to the number of nucleotide bases. We show that our method can give the correct probability distribution when alkylation events are relatively rare, discuss how rare these events have to be and suggest potential extensions to the model when a greater proportion of bases are alkylated.

Chromatin loops are responsible for higher counts of small DNA fragments induced by high-LET radiation, while chromosomal domains do not affect the fragment sizes

PURPOSE

To apply a polymer model of DNA damage induced by high-LET (linear energy transfer) radiation and determine the influence of chromosomal domains and loops on fragment length distribution.

MATERIALS AND METHODS

The yields of DSB (double-strand breaks) induced by high-LET radiation were calculated using a track structure model along with a polymer model of DNA packed in the cell nucleus. The cell nucleus was constructed to include the chromosomal domains and chromatin loops. The latter were generated by the random walk method.

RESULTS AND CONCLUSIONS

We present data for DSB yields per track per cell, DNA fragment sizes, the radial distribution of DSB with respect to the track center, and the distribution of 0, 1, 2, and more DSB from a single particle. Calculations were carried out for a range of particles including He (40 keV/microm), N (225 keV/microm), and Fe ions (150 keV/mum). Situations relevant to PFGE (pulsed-field gel electrophoresis) and microbeam experiments with direct irradiation of the cell nucleus were simulated to demonstrate the applicability of the model. Data show that chromosomal domains do not have a significant influence on fragment-size distribution, while the presence of DNA loops increases the frequencies of smaller fragments by nearly 30% for fragment sizes in the range from 2 kbp (bp = base pair) to 20 kbp.

A Monte Carlo model of DNA double-strand break clustering and rejoining kinetics for the analysis of pulsed-field gel electrophoresis data

In studies of radiation-induced DNA fragmentation and repair, analytical models may provide rapid and easy-to-use methods to test simple hypotheses regarding the breakage and rejoining mechanisms involved. The random breakage model, according to which lesions are distributed uniformly and independently of each other along the DNA, has been the model most used to describe spatial distribution of radiation-induced DNA damage. Recently several mechanistic approaches have been proposed that model clustered damage to DNA. In general, such approaches focus on the study of initial radiation-induced DNA damage and repair, without considering the effects of additional (unwanted and unavoidable) fragmentation that may take place during the experimental procedures. While most approaches, including measurement of total DNA mass below a specified value, allow for the occurrence of background experimental damage by means of simple subtractive procedures, a more detailed analysis of DNA fragmentation necessitates a more accurate treatment. We have developed a new, relatively simple model of DNA breakage and the resulting rejoining kinetics of broken fragments. Initial radiation-induced DNA damage is simulated using a clustered breakage approach, with three free parameters: the number of independently located clusters, each containing several DNA double-strand breaks (DSBs), the average number of DSBs within a cluster (multiplicity of the cluster), and the maximum allowed radius within which DSBs belonging to the same cluster are distributed. Random breakage is simulated as a special case of the DSB clustering procedure. When the model is applied to the analysis of DNA fragmentation as measured with pulsed-field gel electrophoresis (PFGE), the hypothesis that DSBs in proximity rejoin at a different rate from that of sparse isolated breaks can be tested, since the kinetics of rejoining of fragments of varying size may be followed by means of computer simulations. The problem of how to account for background damage from experimental handling is also carefully considered. We have shown that the conventional procedure of subtracting the background damage from the experimental data may lead to erroneous conclusions during the analysis of both initial fragmentation and DSB rejoining. Despite its relative simplicity, the method presented allows both the quantitative and qualitative description of radiation-induced DNA fragmentation and subsequent rejoining of double-stranded DNA fragments.

Extrapolation of the dna fragment-size distribution after high-dose irradiation to predict effects at low doses

The patterns of DSBs induced in the genome are different for sparsely and densely ionizing radiations: In the former case, the patterns are well described by a random-breakage model; in the latter, a more sophisticated tool is needed. We used a Monte Carlo algorithm with a random-walk geometry of chromatin, and a track structure defined by the radial distribution of energy deposition from an incident ion, to fit the PFGE data for fragment-size distribution after high-dose irradiation. These fits determined the unknown parameters of the model, enabling the extrapolation of data for high-dose irradiation to the low doses that are relevant for NASA space radiation research. The randomly-located-clusters formalism was used to speed the simulations. It was shown that only one adjustable parameter, Q, the track efficiency parameter, was necessary to predict DNA fragment sizes for wide ranges of doses. This parameter was determined for a variety of radiations and LETs and was used to predict the DSB patterns at the HPRT locus of the human X chromosome after low-dose irradiation. It was found that high-LET radiation would be more likely than low-LET radiation to induce additional DSBs within the HPRT gene if this gene already contained one DSB.

A review of studies of ionizing radiation-induced double-strand break clustering

Underpinning current models of the mechanisms of the action of radiation is a central role for DNA damage and in particular double-strand breaks (DSBs). For radiations of different LET, there is a need to know the exact yields and distributions of DSBs in human cells. Most measurements of DSB yields within cells now rely on pulsed-field gel electrophoresis as the technique of choice. Previous measurements of DSB yields have suggested that the yields are remarkably similar for different types of radiation with RBE values < or = 1.0. More recent studies in mammalian cells, however, have suggested that both the yield and the spatial distribution of DSBs are influenced by radiation quality. RBE values for DSBs induced by high-LET radiations are greater than 1.0, and the distributions are nonrandom. Underlying this is the interaction of particle tracks with the higher-order chromosomal structures within cell nuclei. Further studies are needed to relate nonrandom distributions of DSBs to their rejoining kinetics. At the molecular level, we need to determine the involvement of clustering of damaged bases with strand breakage, and the relationship between higher-order clustering over sizes of kilobase pairs and above to localized clustering at the DNA level. Overall, these studies will allow us to elucidate whether the nonrandom distributions of breaks produced by high-LET particle tracks have any consequences for their repair and biological effectiveness.