Yield Of Double-strand Breaks Research Articles

A new open-source GPU-based microscopic Monte Carlo simulation tool for the calculations of DNA damages caused by ionizing radiation - Part II: sensitivity and uncertainty analysis.

Calculations of deoxyribonucleic acid (DNA) damages involve many parameters in the computation process. As these parameters are often subject to uncertainties, it is of central importance to comprehensively quantify their impacts on DNA single-strand break (SSB) and double-strand break (DSB) yields. This has been a challenging task due to the required large number of simulations and the relatively low computational efficiency using CPU-based MC packages. In this study, we present comprehensive evaluations on sensitivities and uncertainties of DNA SSB and DSB yields on 12 parameters using our GPU-based MC tool, gMicroMC. We sampled one electron at a time in a water sphere containing a human lymphocyte nucleus and transport the electrons and generated radicals until 2Gy dose was accumulated in the nucleus. We computed DNA damages caused by electron energy deposition events in the physical stage and the hydroxyl radicals at the end of the chemical stage. We repeated the computations by varying 12 parameters: (a) physics cross section, (b) cutoff energy for electron transport, (c)-(e) three branching ratios of hydroxyl radicals in the de-excitation of excited water molecules, (f) temporal length of the chemical stage, (g)-(h) reaction radii for direct and indirect damages, (i) threshold energy defining the threshold damage model to generate a physics damage, (j)-(k) minimum and maximum energy values defining the linear-probability damage model to generate a physics damage, and (l) probability to generate a damage by a radical. We quantified sensitivity of SSB and DSB yields with respect to these parameters for cases with 1.0 and 4.5keV electrons. We further estimated uncertainty of SSB and DSB yields caused by uncertainties of these parameters. Using a threshold of 10% uncertainty as a criterion, threshold energy in the threshold damage model, maximum energy in the linear-probability damage model, and probability for a radical to generate a damage were found to cause large uncertainties in both SSB and DSB yields. The scaling factor of the cross section, cutoff energy, physics reaction radius, and minimum energy in the linear-probability damage model were found to generate large uncertainties in DSB yields. We identified parameters that can generate large uncertainties in the calculations of SSB and DSB yields. Our study could serve as a guidance to reduce uncertainties of parameters and hence uncertainties of the simulation results.

Microdosimetric calculations of the direct DNA damage induced by low energy electrons using the Geant4-DNA Monte Carlo code

To calculate the yield of direct DNA damage induced by low energy electrons using Monte Carlo generated microdosimetric spectra at the nanometer scale and examine the influence of various simulation inputs. The potential of classical microdosimetry to offer a viable and simpler alternative to more elaborate mechanistic approaches for practical applications is discussed. Track-structure simulations with the Geant4-DNA low-energy extension of the Geant4 Monte Carlo toolkit were used for calculating lineal energy spectra in spherical volumes with dimensions relevant to double-strand-break (DSB) induction. The microdosimetric spectra were then used to calculate the yield of simple and clustered DSB based on literature values of the threshold energy of DNA damage. The influence of the different implementations of the dielectric function of liquid water available in Geant4-DNA (Option 2 and Option 4 constructors), as well as the effect of particle tracking cutoff energy and target size are examined.Frequency- and dose-mean lineal energies in liquid-water spheres of 2, 2.3, 2.6, and 3.4 nm diameter, as well as, number of simple and clustered DSB/Gy/cell are presented for electrons over the 100 eV to 100 keV energy range. Results are presented for both the ‘default’ (Option 2) and ‘Ioannina’ (Option 4) physics models of Geant4-DNA applying several commonly used tracking cutoff energies (10, 20, 50, 100 eV). Overall, the choice of the physics model and target diameter has a moderate effect (up to ~10%–30%) on the DSB yield whereas the effect of the tracking cutoff energy may be significant (>100%). Importantly, the yield of both simple and clustered DSB was found to vary significantly (by a factor of 2 or more) with electron energy over the examined range.The yields of electron-induced simple and clustered DSB exhibit a strong energy dependence over the 100 eV–100 keV range with implications to radiation quality issues. It is shown that a classical microdosimetry approach for the calculation of DNA damage based on lineal energy spectra in nanometer-size targets predicts comparable results to computationally intensive mechanistic approaches which use detailed atomistic DNA geometries, thus, offering a relatively simple and robust alternative for some practical applications.

An investigation of the physics mechanism based on DNA damage produced by protons and alpha particles in a realistic DNA model

The direct effect of monoenergetic protons and alpha particles on DNA molecules, as the biological endpoint, were investigated using the Geant4-DNA extensions of the Geant4 Toolkit. two different geometries of B-DNA were simulated to confirm the importance of DNA geometry, namely an atomic model of the Protein Data Bank (PDB) and a simple water cylinder. A direct damage is produced when incident radiations deposit energy in the DNA molecule. Further, the yields of DNA Direct Single-Strand Breaks (SSB), Double-Strand Breaks (DSB), and Total-Strand Breaks (TSB) were calculated for protons and alpha particles with different incident energies. The relative biological effectiveness for the induction of direct DNA double-strand breaks (RBEDSB) was evaluated in the human cell. The SSBevent ratio was determined for each incident particle. The results from the models indicated that TSB, SSB, and DSB yields are higher for simple geometry, compared to the atomic model, and a strong dependence was found by selecting the DNA geometrical model. The SSB yield incorporated the weak dependence of SSB on the type and energy of the incident particles. In addition, based on the results, the TSB yield was approximately independent of the linear energy transfer (LET) of the incident energies while the DSB yield increased as a function of the LET. Further, the results demonstrated that the DSB yield and RBEDSB of the protons were higher for iso-LET ions, compared to the alpha particles. Finally, the contribution of secondary electrons to the creation of DSBs was equal to 83 and 74% for protons and alpha particles at the same LET value, respectively.

PO-1068 Very high yield of double strand breaks found at the distal end of the proton Bragg peak

PO-1068 Very high yield of double strand breaks found at the distal end of the proton Bragg peak

A mechanistic approach towards determining double strand breaks and Relative Biological Effectiveness variation along proton tracks

Relative Biological Effectiveness (RBE) is an important but complex quantity which is used in Proton Therapy to compare the efficacy of cell-killing of protons relative to x-ray. It is known to depend on many factors (such as proton’s energy, dose, and cell type) and varies along the track of a particle. The exact RBE value is also known to be subjected to a large amount of uncertainties. In light of these complexities and uncertainties, this paper presents the calculation of the biologically and clinically relevant parameters - RBE (for double strand breaks) and double strand breaks yield, using a mechanistic and radiobiology-inspired approach. This paper is divided into two main parts. The first features the debut of our in-house radiobiology simulation software known as Integrated Radiobiology Simulation with Geant4 and System Biology (IRSGS). The software and the methodology employed will be covered in details including chromosome geometry construction and processing of direct and indirect effect. The second part focuses on the derivation of the spatial distribution of strand breaks (especially double strand breaks) and RBE for protons.

HOW DETECTION OF PLASMID DNA FRAGMENTATION AFFECTS RADIATION STRAND BREAK YIELDS.

A compromised detection of radiation-induced plasmid DNA fragments results in underestimation of calculated damage yields. Electrophoretic methods are easy and cheap, but they can only detect a part of the fragments, neglecting the shortest ones. These can be detected with atomic force microscopy, but at the expense of time and price. Both methods were used to investigate their capabilities to detect the DNA fragments induced by high-energetic heavy ions. The results were taken into account in calculations of radiation-induced yields of single and double strand breaks. It was estimated that the double strand break yield is twice as high when the fragments are at least partially detected with the agarose electrophoresis, compared to when they were completely omitted. Further increase by 13% was observed when the measured fragments were corrected for the fraction of the shortest fragments up to 300 base pairs, as detected with the atomic force microscopy. The effect of fragment detection on the single strand break yield was diminished.

SIMULATION OF EARLY RADIATION-INDUCED DNA DAMAGE ON DIFFERENT TYPES OF CELL NUCLEI.

This work presents a comparison of simulated early radiation-induced DNA damage represented by yields of double-strand breaks (DSB) in three different human cell nuclei geometries representing fibroblasts, lymphocytes and endothelial cells for protons and alpha particles of different energies and for different irradiation configurations. Each cell nucleus model includes a multi-scale description of the DNA target from the molecular level to the whole human genome representation (6 Gbp) in the G0/G1 phase of the cell cycle and was generated with the DnaFabric software. The three nuclei differ in shape, volume, and therefore DNA density. A calculation chain based on Geant4-DNA that takes into account the physical, physico-chemical and chemical stages was used to simulate the irradiation of the different cell nuclei. Results show an increase of DSB/primary/μm with an increase of DNA density and an increase of DSB/Gy/Gbp with an increase of the cell nucleus volume which indicates that the cell nucleus shape and size have an impact on early DNA damage, which may play a role in latter effects.

Complex interplay of lesion-specific DNA repair enzyme on bistranded clustered DNA damage harboring Tg:G mismatch in nucleosome core particles

5,6-Dihydroxy-5,6-dihydrothymine (thymine glycol) and 7,8-dihydro-8-oxo-20-deoxyguanosine (8-oxodG) are major DNA damage lesions produced by endogenous oxidative stress, as well as inflicted by carcinogens and ionizing radiation. The processing of Tg:G mismatch and 8-oxodG in close proximity of each other in a bistranded clustered environment in DNA oligomer duplexes as well as in a nucleosome core particle (NCP) model are reported here. The processing of the lesions was evaluated by purified enzyme cocktails of hNTH1 and hOGG1 as well as with a HeLa cell extract. Interestingly, the yield of double-strand breaks (DSBs) resulting from the processing of the bistranded lesions are appreciably lower when the DNA is treated with the HeLa cell extract compared with the relevant purified enzyme cocktail in both models. Clustered bistranded lesions become more repair refractive when reconstituted as an NCP. This indicates a complex interplay between the repair enzymes that influence the processing of the bistranded cluster damage positively to avoid the formation of DSBs under cellular conditions. In addition to position and orientation of the lesions, the type of the lesions in the cluster environment in DNA along with the relative abundance of the lesion-specific enzymes in the cells strongly prevents the processing of the oxidized nucleobases.

Targeted alpha therapy with 212Pb or 225Ac: Change in RBE from daughter migration

Targeted α-therapy (TAT) could be delivered early to patients who are at a high-risk for developing brain metastases, targeting the areas of the vasculature where tumor cells are penetrating into the brain. We have utilized a Monte Carlo model representing brain vasculature to calculate physical dose and DNA damage from the α-emitters 225Ac and 212Pb. The micron-scale dose distributions from all radioactive decay products were modeled in Geant4, including the eV-scale interactions using the Geant4-DNA models. These interactions were then superimposed on an atomic-scale DNA model to estimate strand break yields. In addition to 225Ac having a higher dose per decay than 212Pb, it also has a double strand break yield per decay that is 4.7 ± 0.5 times that of 212Pb. However, the efficacy of both nuclides depends on retaining the daughter nuclei at the target location in the brain vasculature. The relative biological effectiveness (RBE) of 225Ac and 212Pb are similar when the entire decay chains are included, with maxima of 2.7 ± 0.6 and 2.5 ± 0.5 (respectively), and RBE values of about 2 to a depth of 80 μm. If the initial daughter is lost, the RBE of 212Pb is completely reduced to 1 or lower and the RBE of 225Ac is approximately 2 only for the first 40 μm.

The Moss Physcomitrella patens Is Hyperresistant to DNA Double-Strand Breaks Induced by γ-Irradiation.

The purpose of this study was to investigate whether the moss Physcomitrella patens cells are more resistant to ionizing radiation than animal cells. Protoplasts derived from P. patens protonemata were irradiated with γ-rays of 50–1000 gray (Gy). Clonogenicity of the protoplasts decreased in a γ-ray dose-dependent manner. The dose that decreased clonogenicity by half (LD50) was 277 Gy, which indicated that the moss protoplasts were 200-times more radioresistant than human cells. To investigate the mechanism of radioresistance in P. patens, we irradiated protoplasts on ice and initial double-strand break (DSB) yields were measured using the pulsed-field gel electrophoresis assay. Induced DSBs linearly increased dependent on the γ-ray dose and the DSB yield per Gb DNA per Gy was 2.2. The DSB yield in P. patens was half to one-third of those reported in mammals and yeasts, indicating that DSBs are difficult to induce in P. patens. The DSB yield per cell per LD50 dose in P. patens was 311, which is three- to six-times higher than those in mammals and yeasts, implying that P. patens is hyperresistant to DSBs. Physcomitrella patens is indicated to possess unique mechanisms to inhibit DSB induction and provide resistance to high numbers of DSBs.

TOPAS-nBio: An Extension to the TOPAS Simulation Toolkit for Cellular and Sub-cellular Radiobiology.

The TOPAS Monte Carlo (MC) system is used in radiation therapy and medical imaging research, having played a significant role in making Monte Carlo simulations widely available for proton therapy related research. While TOPAS provides detailed simulations of patient scale properties, the fundamental unit of the biological response to radiation is a cell. Thus, our goal was to develop TOPAS-nBio, an extension of TOPAS dedicated to advance understanding of radiobiological effects at the (sub-)cellular, (i.e., the cellular and sub-cellular) scale. TOPAS-nBio was designed as a set of open source classes that extends TOPAS to model radiobiological experiments. TOPAS-nBio is based on and extends Geant4-DNA, which extends the Geant4 toolkit, the basis of TOPAS, to include very low-energy interactions of particles down to vibrational energies, explicitly simulates every particle interaction (i.e., without using condensed histories) and propagates radiolysis products. To further facilitate the use of TOPAS-nBio, a graphical user interface was developed. TOPAS-nBio offers full track-structure Monte Carlo simulations, integration of chemical reactions within the first millisecond, an extensive catalogue of specialized cell geometries as well as sub-cellular structures such as DNA and mitochondria, and interfaces to mechanistic models of DNA repair kinetics. We compared TOPAS-nBio simulations to measured and published data of energy deposition patterns and chemical reaction rates (G values). Our simulations agreed well within the experimental uncertainties. Additionally, we expanded the chemical reactions and species provided in Geant4-DNA and developed a new method based on independent reaction times (IRT), including a total of 72 reactions classified into 6 types between neutral and charged species. Chemical stage simulations using IRT were a factor of 145 faster than with step-by-step tracking. Finally, we applied the geometric/chemical modeling to obtain initial yields of double-strand breaks (DSBs) in DNA fibers for proton irradiations of 3 and 50 MeV and compared the effect of including chemical reactions on the number and complexity of DSB induction. Over half of the DSBs were found to include chemical reactions with approximately 5% of DSBs caused only by chemical reactions. In conclusion, the TOPAS-nBio extension to the TOPAS MC application offers access to accurate and detailed multiscale simulations, from a macroscopic description of the radiation field to microscopic description of biological outcome for selected cells. TOPAS-nBio offers detailed physics and chemistry simulations of radiobiological experiments on cells simulating the initially induced damage and links to models of DNA repair kinetics.

Dosimetric evaluation of radionuclides for VCAM-1-targeted radionuclide therapy of early brain metastases.

Brain metastases develop frequently in patients with breast cancer, and present a pressing therapeutic challenge. Expression of vascular cell adhesion molecule 1 (VCAM-1) is upregulated on brain endothelial cells during the early stages of metastasis and provides a target for the detection and treatment of early brain metastases. The aim of this study was to use a model of early brain metastasis to evaluate the efficacy of α-emitting radionuclides, 149Tb, 211At, 212Pb, 213Bi and 225Ac; β-emitting radionuclides, 90Y, 161Tb and 177Lu; and Auger electron (AE)-emitters 67Ga, 89Zr, 111In and 124I, for targeted radionuclide therapy (TRT).METHODS: Histologic sections and two photon microscopy of mouse brain parenchyma were used to inform a cylindrical vessel geometry using the Geant4 general purpose Monte Carlo (MC) toolkit with the Geant4-DNA low energy physics models. Energy deposition was evaluated as a radial function and the resulting phase spaces were superimposed on a DNA model to estimate double-strand break (DSB) yields for representative β- and α-emitters, 177Lu and 212Pb. Relative biological effectiveness (RBE) values were determined by only evaluating DNA damage due to physical interactions.RESULTS: 177Lu produced 2.69 ± 0.08 DSB per GbpGy, without significant variation from the lumen of the vessel to a radius of 100 µm. The DSB yield of 212Pb included two local maxima produced by the 6.1 MeV and 8.8 MeV α-emissions from decay products, 212Bi and 212Po, with yields of 7.64 ± 0.12 and 9.15 ± 0.24 per GbpGy, respectively. Given its higher DSB yield 212Pb may be more effective for short range targeting of early micrometastatic lesions than 177Lu.CONCLUSION: MC simulation of a model of early brain metastases provides invaluable insight into the potential efficacy of α-, β- and AE-emitting radionuclides for TRT. 212Pb, which has the attributes of a theranostic radionuclide since it can be used for SPECT imaging, showed a favorable dose profile and RBE.

Modeling Radiation Effects of Ultrasoft X Rays on the Basis of Amorphous Track Structure.

There is experimental evidence that ultrasoft X rays (0.1-5 keV) show a higher biological effectiveness than high-energy photons. Similar to high-LET radiation, this is attributed to a rather localized dose distribution associated with a considerably smaller range of secondary electrons, which results in an increasing yield of double-strand breaks (DSBs) and potentially more complex lesions. We previously reported on the application of the Giant LOop Binary LEsion (GLOBLE) model to ultrasoft X rays, in which experimental values of the relative biological effectiveness (RBE) for DSB induction were used to show that this increasing DSB yield was sufficient to explain the enhanced effectiveness in the cell inactivation potential of ultrasoft X rays. Complementary to GLOBLE, we report here on a modeling approach to predict the increased DSB yield of ultrasoft X rays on the basis of amorphous track structure formed by secondary electrons, which was derived from Monte Carlo track structure simulations. This procedure is associated with increased production of single-strand break (SSB) clusters, which are caused by the highly localized energy deposition pattern induced by low-energy photons. From this, the RBE of ultrasoft X rays can be determined and compared to experimental data, showing that the inhomogeneity of the energy deposition pattern represents the key variable to describe the increased biological effectiveness of ultrasoft X rays. Thus, this work demonstrates an extended applicability of the amorphous track structure concept and tests its limits with respect to its predictive power. The employed model mechanism offers a possible explanation for how the cellular response to ultrasoft X rays is directly linked to the energy deposition properties on the nanometric scale.

Comprehensive track-structure based evaluation of DNA damage by light ions from radiotherapy-relevant energies down to stopping

Track structures and resulting DNA damage in human cells have been simulated for hydrogen, helium, carbon, nitrogen, oxygen and neon ions with 0.25–256 MeV/u energy. The needed ion interaction cross sections have been scaled from those of hydrogen; Barkas scaling formula has been refined, extending its applicability down to about 10 keV/u, and validated against established stopping power data. Linear energy transfer (LET) has been scored from energy deposits in a cell nucleus; for very low-energy ions, it has been defined locally within thin slabs. The simulations show that protons and helium ions induce more DNA damage than heavier ions do at the same LET. With increasing LET, less DNA strand breaks are formed per unit dose, but due to their clustering the yields of double-strand breaks (DSB) increase, up to saturation around 300 keV/μm. Also individual DSB tend to cluster; DSB clusters peak around 500 keV/μm, while DSB multiplicities per cluster steadily increase with LET. Remarkably similar to patterns known from cell survival studies, LET-dependencies with pronounced maxima around 100–200 keV/μm occur on nanometre scale for sites that contain one or more DSB, and on micrometre scale for megabasepair-sized DNA fragments.

Efficiency of radiation-induced base lesion excision and the order of enzymatic treatment

Purpose: To clarify whether initial base excision repair processes at clustered DNA damage sites comprising multiple base lesions affect subsequent excision processes via the formation of additional strand breaks by glycosylase and apurinic/apyrimidinic (AP) endonuclease base excision enzymes.Materials and methods: Plasmid DNA (pUC18) as a model DNA molecule was exposed to high-linear-energy-transfer (LET) ionizing radiation (He2+ or C6+ ions) or low-LET ionizing radiation (X-rays) under various conditions to produce varied radical-scavenging effects. pUC18 was then treated sequentially or simultaneously with two bacterial base excision enzymes (glycosylases), namely, endonuclease III and formamidopyrimidine-DNA glycosylase, which convert pyrimidine (or abasic [AP] site) and purine (or AP site) lesions to single-strand breaks (SSB), respectively. Yields of additional SSB or double-strand breaks (DSB) as digestion products were examined after changing the order of enzymatic treatment.Results: There were few differences among the enzymatic treatments, indicating that treatment order did not affect the final yields of additional SSB or DSB formed by glycosylase activity. This suggests that of the total damage, the fraction of clustered damage sites with a persistent base lesion dependent on the order of glycosylase treatment was insignificant if present.Conclusion: Base lesion clusters induced by high- or low-LET radiation appear three or more base pairs apart, and are promptly converted to a DSB by glycosylase, regardless of the order of enzymatic treatment.

Proton and light ion RBE for the induction of direct DNA double strand breaks.

To present and characterize a Monte Carlo (MC) tool for the simulation of the relative biological effectiveness for the induction of direct DNA double strand breaks (RBEDSB (direct)) for protons and light ions. The MC tool uses a pregenerated event-by-event tracks library of protons and light ions that are overlaid on a cell nucleus model. The cell nucleus model is a cylindrical arrangement of nucleosome structures consisting of 198 DNA base pairs. An algorithm relying on k-dimensional trees and cylindrical symmetries is used to search coincidences of energy deposition sites with volumes corresponding to the sugar-phosphate backbone of the DNA molecule. Strand breaks (SBs) are scored when energy higher than a threshold is reached in these volumes. Based on the number of affected strands, they are categorized into either single strand break (SSB) or double strand break (DSB) lesions. The number of SBs composing each lesion (i.e., its size) is also recorded. RBEDSB (direct) is obtained by taking the ratio of DSB yields of a given radiation field to a (60)Co field. The MC tool was used to obtain SSB yields, DSB yields, and RBEDSB (direct) as a function of linear energy transfer (LET) for protons ((1)H(+)), (4)He(2+), (7)Li(3+), and (12)C(6+) ions. For protons, the SSB yields decreased and the DSB yields increased with LET. At ≈24.5 keV μm(-1), protons generated 15% more DSBs than (12)C(6+) ions. The RBEDSB (direct) varied between 1.24 and 1.77 for proton fields between 8.5 and 30.2 keV μm(-1), and it was higher for iso-LET ions with lowest atomic number. The SSB and DSB lesion sizes showed significant differences for all radiation fields. Generally, the yields of SSB lesions of sizes ≥2 and the yields of DSB lesions of sizes ≥3 increased with LET and increased for iso-LET ions of lower atomic number. On the other hand, the ratios of SSB to DSB lesions of sizes 2-4 did not show variability with LET nor projectile atomic number, suggesting that these metrics are independent of the radiation quality. Finally, a variance of up to 8% in the DSB yields was observed as a function of the particle incidence angle on the cell nucleus. This simulation effect is due to the preferential alignment of ion tracks with the DNA nucleosomes at specific angles. The MC tool can predict SSB and DSB yields for light ions of various LET and estimate RBEDSB (direct). In addition, it can calculate the frequencies of different DNA lesion sizes, which is of interest in the context of biologically relevant absolute dosimetry of particle beams.

Microdosimetry of DNA conformations: relation between direct effect of (60)Co gamma rays and topology of DNA geometrical models in the calculation of A-, B- and Z-DNA radiation-induced damage yields.

In order to obtain the energy deposition pattern of ionizing radiation in the nanometric scale of genetic material and to investigate the different sensitivities of the DNA conformations, direct effects of (60)Co gamma rays on the three A, B and Z conformations of DNA have been studied. For this purpose, single-strand breaks (SSB), double-strand breaks (DSB), base damage (BD), hit probabilities and three microdosimetry quantities (imparted energy, mean chord length and lineal energy) in the mentioned DNA conformations have been calculated and compared by using GEometry ANd Tracking 4 (Geant4) toolkit. The results show that A-, B- and Z-DNA conformations have the highest yields of DSB (1.2 Gy(-1) Gbp(-1)), SSB (25.2 Gy(-1) Gbp(-1)) and BD (4.81 Gy(-1) Gbp(-1)), respectively. Based on the investigation of direct effects of radiation, it can be concluded that the DSB yield is largely correlated to the topological characteristics of DNA models, although the SSB yield is not. Moreover, according to the comparative results of the present study, a reliable candidate parameter for describing the relationship between DNA damage yields and geometry of DNA models in the theoretical radiation biology research studies would be the mean chord length (4 V/S) of the models.

Monte Carlo simulation of ionizing radiation induced DNA strand breaks utilizing coarse grained high-order chromatin structures

Ionizing radiation threatens genome integrity by causing DNA damage. Monte Carlo simulation of the interaction of a radiation track structure with DNA provides a powerful tool for investigating the mechanisms of the biological effects. However, the more or less oversimplification of the indirect effect and the inadequate consideration of high-order chromatin structures in current models usually results in discrepancies between simulations and experiments, which undermine the predictive role of the models. Here we present a biophysical model taking into consideration factors that influence indirect effect to simulate radiation-induced DNA strand breaks in eukaryotic cells with high-order chromatin structures. The calculated yields of single-strand breaks and double-strand breaks (DSBs) for photons are in good agreement with the experimental measurements. The calculated yields of DSB for protons and α particles are consistent with simulations by the PARTRAC code, whereas an overestimation is seen compared with the experimental results. The simulated fragment size distributions for 60Co γ irradiation and α particle irradiation are compared with the measurements accordingly. The excellent agreement with 60Co irradiation validates our model in simulating photon irradiation. The general agreement found in α particle irradiation encourages model applicability in the high linear energy transfer range. Moreover, we demonstrate the importance of chromatin high-order structures in shaping the spectrum of initial damage.

Mechanisms and Consequences of Double-Strand DNA Break Formation in Chromatin.

All organisms suffer double-strand breaks (DSBs) in their DNA as a result of exposure to ionizing radiation. DSBs can also form when replication forks encounter DNA lesions or repair intermediates. The processing and repair of DSBs can lead to mutations, loss of heterozygosity, and chromosome rearrangements that result in cell death or cancer. The most common pathway used to repair DSBs in metazoans (non-homologous DNA end joining) is more commonly mutagenic than the alternative pathway (homologous recombination mediated repair). Thus, factors that influence the choice of pathways used DSB repair can affect an individual's mutation burden and risk of cancer. This review describes radiological, chemical, and biological mechanisms that generate DSBs, and discusses the impact of such variables as DSB etiology, cell type, cell cycle, and chromatin structure on the yield, distribution, and processing of DSBs. The final section focuses on nucleosome-specific mechanisms that influence DSB production, and the possible relationship between higher order chromosome coiling and chromosome shattering (chromothripsis).

Sub-micrometer 20 MeV protons or 45 MeV lithium spot irradiation enhances yields of dicentric chromosomes due to clustering of DNA double-strand breaks

In conventional experiments on biological effects of radiation types of diverse quality, micrometer-scale double-strand break (DSB) clustering is inherently interlinked with clustering of energy deposition events on nanometer scale relevant for DSB induction. Due to this limitation, the role of the micrometer and nanometer scales in diverse biological endpoints cannot be fully separated. To address this issue, hybrid human-hamster AL cells have been irradiated with 45MeV (60keV/μm) lithium ions or 20MeV (2.6keV/μm) protons quasi-homogeneously distributed or focused to 0.5×1μm2 spots on regular matrix patterns (point distances up to 10.6×10.6μm), with pre-defined particle numbers per spot to provide the same mean dose of 1.7Gy. The yields of dicentrics and their distribution among cells have been scored. In parallel, track-structure based simulations of DSB induction and chromosome aberration formation with PARTRAC have been performed. The results show that the sub-micrometer beam focusing does not enhance DSB yields, but significantly affects the DSB distribution within the nucleus and increases the chance to form DSB pairs in close proximity, which may lead to increased yields of chromosome aberrations. Indeed, the experiments show that focusing 20 lithium ions or 451 protons per spot on a 10.6μm grid induces two or three times more dicentrics, respectively, than a quasi-homogenous irradiation. The simulations reproduce the data in part, but in part suggest more complex behavior such as saturation or overkill not seen in the experiments. The direct experimental demonstration that sub-micrometer clustering of DSB plays a critical role in the induction of dicentrics improves the knowledge on the mechanisms by which these lethal lesions arise, and indicates how the assumptions of the biophysical model could be improved. It also provides a better understanding of the increased biological effectiveness of high-LET radiation.