Spatial distribution and yield of DNA double-strand breaks induced by 3-7 MeV helium ions in human fibroblasts

Development of modified microdosimetric kinetic model for relative biological effectiveness in proton therapy

To predict the biological effects of ionising radiation, the quantity of biological dose is introduced instead of the physical absorbed dose. In proton therapy, a constant relative biological effectiveness (RBE) of 1.1 is usually applied clinically as recommended by the International Commission of Radiation Units and Measurements. This study presents a new model, based on the modified microdosimetric kinetic model (MMKM), for calculating variable RBE values based on experimental data on the induction of DNA double-strand breaks (DSBs) within cells. The MMKM was proposed based on experimental data for the yield of DSBs in mammalian cells, which allows modification of the yield of primary lesions in the MKM. In this approach, a unique function named f(LET), which describes the relation between RBE and linear energy transfer (LET), was considered for charged particles. In the presented model (DMMKM), the MMKM approach was developed further by considering different f(LET)s for different relevant ions involved in energy deposition events in proton therapy. Although experimental data represent the dependence of the yield of primary lesions on the ion species, the DSB yield (assumed as the main primary lesion) is assumed independent of the ion species in the MMKM. In the DMMKM, by considering the yield of primary lesions as a function of the ion species, the α and β values are in better agreement with the experimental data as compared to those of the MKM and MMKM approaches. The biological dose in the DMMKM is predicted to be lower than that in the MMKM. Further, in the proposed model, the variation of the β parameter is higher than the constant value assumed in the MKM, at the distal end of the spread-out Bragg peak (SOBP). Moreover, the level of cell death estimated by the MMKM at the SOBP region is higher than that obtained based on the DMMKM. It is concluded that considering modified f(LET)s in the model developed here is more consistent with experimental results than when MMKM and MKM approaches are considered. The DMMKM examines the biological effects with full detail and will, therefore, be effective in improving proton therapy.

Chromatin and the Cellular Response to Particle Radiation-Induced Oxidative and Clustered DNA Damage

Exposure to environmental ionizing radiation is prevalent, with greatest lifetime doses typically from high Linear Energy Transfer (high-LET) alpha particles via the radioactive decay of radon gas in indoor air. Particle radiation is highly genotoxic, inducing DNA damage including oxidative base lesions and DNA double strand breaks. Due to the ionization density of high-LET radiation, the consequent damage is highly clustered wherein ≥2 distinct DNA lesions occur within 1–2 helical turns of one another. These multiply-damaged sites are difficult for eukaryotic cells to resolve either quickly or accurately, resulting in the persistence of DNA damage and/or the accumulation of mutations at a greater rate per absorbed dose, relative to lower LET radiation types. The proximity of the same and different types of DNA lesions to one another is challenging for DNA repair processes, with diverse pathways often confounding or interplaying with one another in complex ways. In this context, understanding the state of the higher order chromatin compaction and arrangements is essential, as it influences the density of damage produced by high-LET radiation and regulates the recruitment and activity of DNA repair factors. This review will summarize the latest research exploring the processes by which clustered DNA damage sites are induced, detected, and repaired in the context of chromatin.

A forty-year journey from "classical" biophysics and radiobiology to hadrontherapy, space radiation and low dose rate underground radiobiology

PURPOSE

As a biologist who, since the beginning of her involvement in science, has collaborated closely with physicists, I want to share my forty years of experience describing the events that introduced me to the world of charged particle radiation biology as well as that of low doses/dose rates, with related implications in medicine and radiation protection.

CONCLUSION

The main features of my experience can be summarized in the development of an interdisciplinary culture and in the interest in technological advances for the study of biological responses to radiation in different scenarios, relevant for public health. Mine was a journey that began by chance, but which led me to a world that proved to be of great interest to me. With the current advances in science, the new generations of scientists have new opportunities that I wish them to face with the same interest and enthusiasm that I felt for such an interdisciplinary field as that of radiation biology.

Challenges in the quantification approach to a radiation relevant adverse outcome pathway for lung cancer

PURPOSE

Adverse outcome pathways (AOPs) provide a modular framework for describing sequences of biological key events (KEs) and key event relationships (KERs) across levels of biological organization. Empirical evidence across KERs can support construction of quantified AOPs (qAOPs). Using an example AOP of energy deposition from ionizing radiation onto DNA leading to lung cancer incidence, we investigate the feasibility of quantifying data from KERs supported by all types of stressors. The merits and challenges of this process in the context of AOP construction are discussed.

MATERIALS AND METHODS

Empirical evidence across studies of dose-response from four KERs of the AOP were compiled independently for quantification. Three upstream KERs comprised of evidence from various radiation types in line with AOP guidelines. For these three KERs, a focused analysis of data from alpha-particle studies was undertaken to better characterize the process to the adverse outcome (AO) for a radon gas stressor. Numerical information was extracted from tables and graphs to plot and tabulate the response of KEs. To complement areas of the AOP quantification process, Monte Carlo (MC) simulations in TOPAS-nBio were performed to model exposure conditions relevant to the AO for an example bronchial compartment of the lung with secretory cell nuclei targets.

RESULTS

Quantification of AOP KERs highlighted the relevance of radiation types under the stressor-agnostic intent of AOP design, motivating a focus on specific types. For a given type, significant differences of KE response indicate meaningful data to derive linkages from the MIE to the AO is lacking and that better response-response focused studies are required. The MC study estimates the linear energy transfer (LET) of alpha-particles emitted by radon-222 and its progeny in the secretory cell nuclei of the example lung compartment to range from 94 - 5 + 5 to 192 - 18 + 15 keV/µm.

CONCLUSION

Quantifying AOP components provides a means to assemble empirical evidence across different studies. This highlights challenges in the context of studies examining similar endpoints using different radiation types. Data linking KERs to a MIE of 'deposition of energy' is shown to be non-compatible with the stressor-agnostic principles of AOP design. Limiting data to that describing response-response relationships between adjacent KERs may better delineate studies relevant to the damage that drives a pathway to the next KE and still support an 'all hazards' approach. Such data remains limited and future investigations in the radiation field may consider this approach when designing experiments and reporting their results and outcomes.

An Assessment of Radiation Doses From Radon Exposures Using a Mouse Model System

BACKGROUND

Radon and its progenies contribute significantly to the natural background radiation and cause several thousands of lung cancer cases per year worldwide. Moreover, patients with chronic inflammatory joint diseases are treated in radon galleries. Due to the complex nature of radon exposure, the doses associated with radon exposures are difficult to assess. Hence, there is a clear need to directly measure dose depositions from radon exposures to provide reliable risk estimates for radiation protection guidelines.

OBJECTIVES

We aimed to assess tissue-specific radiation doses associated with radon activity concentrations, that deposit similar dose levels as the annual natural radon exposure or radon gallery visits.

METHODS

We exposed mice to defined radon concentrations, quantified the number of 53BP1 foci as a measure of induced DNA damage, and compared it with the number of foci induced by known doses of reference-type radiations. An image-based analysis of the 3-dimensional foci pattern provided information about the radiation type inflicting the DNA damage.

RESULTS

A 1-hour exposure to 440 kBq/m³ radon-induced DNA damage corresponding to a dose of ∼10 mGy in the lung and ∼3.3 mGy in the kidney, heart, and liver. A 1-hour exposure to 44 kBq/m³ provided values consistent with a linear relationship between dose and radon concentration. Two-thirds of the dose in the lung was caused by α-particles. The dose in the kidney, heart, and liver and one-third of the dose in the lung likely resulted from β- and γ-rays.

DISCUSSION

We found that radon exposures mainly lead to α-particle-induced DNA damage in the lung, consistent with the lung cancer risk obtained in epidemiologic studies. Our presented biodosimetric approach can be used to benchmark risk model calculations for radiation protection guidelines and can help to understand the therapeutic success of radon gallery treatments.

Adverse outcome pathways for ionizing radiation and breast cancer involve direct and indirect DNA damage, oxidative stress, inflammation, genomic instability, and interaction with hormonal regulation of the breast

Knowledge about established breast carcinogens can support improved and modernized toxicological testing methods by identifying key mechanistic events. Ionizing radiation (IR) increases the risk of breast cancer, especially for women and for exposure at younger ages, and evidence overall supports a linear dose-response relationship. We used the Adverse Outcome Pathway (AOP) framework to outline and evaluate the evidence linking ionizing radiation with breast cancer from molecular initiating events to the adverse outcome through intermediate key events, creating a qualitative AOP. We identified key events based on review articles, searched PubMed for recent literature on key events and IR, and identified additional papers using references. We manually curated publications and evaluated data quality. Ionizing radiation directly and indirectly causes DNA damage and increases production of reactive oxygen and nitrogen species (RONS). RONS lead to DNA damage and epigenetic changes leading to mutations and genomic instability (GI). Proliferation amplifies the effects of DNA damage and mutations leading to the AO of breast cancer. Separately, RONS and DNA damage also increase inflammation. Inflammation contributes to direct and indirect effects (effects in cells not directly reached by IR) via positive feedback to RONS and DNA damage, and separately increases proliferation and breast cancer through pro-carcinogenic effects on cells and tissue. For example, gene expression changes alter inflammatory mediators, resulting in improved survival and growth of cancer cells and a more hospitable tissue environment. All of these events overlap at multiple points with events characteristic of "background" induction of breast carcinogenesis, including hormone-responsive proliferation, oxidative activity, and DNA damage. These overlaps make the breast particularly susceptible to ionizing radiation and reinforce that these biological activities are important characteristics of carcinogens. Agents that increase these biological processes should be considered potential breast carcinogens, and predictive methods are needed to identify chemicals that increase these processes. Techniques are available to measure RONS, DNA damage and mutation, cell proliferation, and some inflammatory proteins or processes. Improved assays are needed to measure GI and chronic inflammation, as well as the interaction with hormonally driven development and proliferation. Several methods measure diverse epigenetic changes, but it is not clear which changes are relevant to breast cancer. In addition, most toxicological assays are not conducted in mammary tissue, and so it is a priority to evaluate if results from other tissues are generalizable to breast, or to conduct assays in breast tissue. Developing and applying these assays to identify exposures of concern will facilitate efforts to reduce subsequent breast cancer risk.

Necessities in the Processing of DNA Double Strand Breaks and Their Effects on Genomic Instability and Cancer

Double strand breaks (DSBs) are induced in the DNA following exposure of cells to ionizing radiation (IR) and are highly consequential for genome integrity, requiring highly specialized modes of processing. Erroneous processing of DSBs is a cause of cell death or its transformation to a cancer cell. Four mechanistically distinct pathways have evolved in cells of higher eukaryotes to process DSBs, providing thus multiple options for the damaged cells. The homologous recombination repair (HRR) dependent subway of gene conversion (GC) removes IR-induced DSBs from the genome in an error-free manner. Classical non-homologous end joining (c-NHEJ) removes DSBs with very high speed but is unable to restore the sequence at the generated junction and can catalyze the formation of translocations. Alternative end-joining (alt-EJ) operates on similar principles as c-NHEJ but is slower and more error-prone regarding both sequence preservation and translocation formation. Finally, single strand annealing (SSA) is associated with large deletions and may also form translocations. Thus, the four pathways available for the processing of DSBs are not alternative options producing equivalent outcomes. We discuss the rationale for the evolution of pathways with such divergent properties and fidelities and outline the logic and necessities that govern their engagement. We reason that cells are not free to choose one specific pathway for the processing of a DSB but rather that they engage a pathway by applying the logic of highest fidelity selection, adapted to necessities imposed by the character of the DSB being processed. We introduce DSB clusters as a particularly consequential form of chromatin breakage and review findings suggesting that this form of damage underpins the increased efficacy of high linear energy transfer (LET) radiation modalities. The concepts developed have implications for the protection of humans from radon-induced cancer, as well as the treatment of cancer with radiations of high LET.

Accounting for radiation-induced indirect damage on DNA with the Geant 4-DNA code

The use of Monte Carlo (MC) simulations remains a powerful tool to study the biological effects induced by ionizing radiation on living beings. Several MC codes are commonly used in research fields such as nanodosimetry, radiotherapy, radiation protection, and space radiation. This work presents an enhancement of an existing model [1] for radiobiological purposes, to account for the indirect DNA damage induced by ionizing particles. The Geant4-DNA simulation toolkit was used to simulate the physical, pre-chemical, and chemical stages of early DNA damage induced by protons and α-particles. Liquid water was used as the medium for simulations. Two phase-space files were generated, one containing the energy deposition events and another with the position of chemical species produced by water radiolysis from 0.1 ps up to 1 ns. These files were used as input in the radiobiological code that contains the genetic material model with atomic resolution, consisting of several copies of 30 nm chromatin fibers. The B-DNA configuration was used. This work focused on the indirect damage produced by the hydroxyl radical (OH) attack on the sugar-phosphate group. The approach followed to account for the indirect DNA damage was the same as those used by other radiobiological codes [2,3]. The critical parameter considered here was the reaction radius, which was calculated from the Smoluchowski's diffusion equation. Single, double, and total strand break yields produced by direct, indirect, and mixed mechanisms are reported. The obtained results are consistent with experimental and calculation data sets published in the literature.

Proximity effects in chromosome aberration induction: Dependence on radiation quality, cell type and dose

It is widely accepted that, in chromosome-aberration induction, the (mis-)rejoining probability of two chromosome fragments depends on their initial distance, r. However, several aspects of these "proximity effects" need to be clarified, also considering that they can vary with radiation quality, cell type and dose. A previous work performed by the BIANCA (BIophysical ANalysis of Cell death and chromosome Aberrations) biophysical model has suggested that, in human lymphocytes and fibroblasts exposed to low-LET radiation, an exponential function of the form exp(-r/r₀), which is consistent with free-end (confined) diffusion, describes proximity effects better than a Gaussian function. Herein, the investigation was extended to intermediate- and high-LET. Since the r₀ values (0.8 μm for lymphocytes and 0.7 μm for fibroblasts) were taken from the low-LET study, the results were obtained by adjusting only one model parameter, i.e. the yield of "Cluster Lesions" (CLs), where a CL was defined as a critical DNA damage producing two independent chromosome fragments. In lymphocytes, the exponential model allowed reproducing both dose-response curves for different aberrations (dicentrics, centric rings and excess acentrics), and values of F-ratio (dicentrics to centric rings) and G-ratio (interstitial deletions to centric rings). In fibroblasts, a good correspondence was found with the dose-response curves, whereas the G-ratio (and, to a lesser extent, the F-ratio) was underestimated. With increasing LET, F decreased and G increased in both cell types, supporting their role as "fingerprints" of high-LET exposure. A dose-dependence was also found at high LET, where F increased with dose and G decreased, possibly due to inter-track effects. We therefore conclude that, independent of radiation quality, in lymphocytes an exponential function can describe proximity effects at both inter- and intra-chromosomal level; on the contrary, in fibroblasts further studies (experimental and theoretical) are needed to explain the strong bias for intra-arm relative to inter-arm exchanges.

Proximity effects in chromosome aberration induction by low-LET ionizing radiation

Although chromosome aberrations are known to derive from distance-dependent mis-rejoining of chromosome fragments, evaluating whether a certain model describes such "proximity effects" better than another one is complicated by the fact that different approaches have often been tested under different conditions. Herein, a biophysical model ("BIANCA", i.e. BIophysical ANalysis of Cell death and chromosome Aberrations) was upgraded, implementing explicit chromosome-arm domains and two new models for the dependence of the rejoining probability on the fragment initial distance, r. Such probability was described either by an exponential function like exp(-r/r₀), or by a Gaussian function like exp(-r²/2σ²), where r₀ and σ were adjustable parameters. The second, and last, parameters was the yield of "Cluster Lesions" (CL), where "Cluster Lesion" defines a critical DNA damage producing two independent chromosome fragments. The model was applied to low-LET-irradiated lymphocytes (doses: 1-4Gy) and fibroblasts (1-6.1Gy). Good agreement with experimental yields of dicentrics and centric rings, and thus their ratio ("F-ratio"), was found by both the exponential model (with r₀=0.8μm for lymphocytes and 0.7μm for fibroblasts) and the Gaussian model (with σ=1.1μm for lymphocytes and 1.3μm for fibroblasts). While the former also allowed reproducing dose-responses for excess acentric fragments, the latter substantially underestimated the experimental curves. Both models provided G-ratios (ratio of acentric to centric rings) higher than those expected from randomness, although the values calculated by the Gaussian model were lower than those calculated by the exponential one. For lymphocytes the calculated G-ratios were in good agreement with the experimental ones, whereas for fibroblasts both models substantially underestimated the experimental results, which deserves further investigation. This work suggested that, although both models performed better than a step model (which previously allowed reproducing the F-ratio but underestimated the G-ratio), an exponential function describes proximity effects better than a Gaussian one.

Modeling radiation-induced cell death: role of different levels of DNA damage clustering

Some open questions on the mechanisms underlying radiation-induced cell death were addressed by a biophysical model, focusing on DNA damage clustering and its consequences. DNA "cluster lesions" (CLs) were assumed to produce independent chromosome fragments that, if created within a micrometer-scale threshold distance (d), can lead to chromosome aberrations following mis-rejoining; in turn, certain aberrations (dicentrics, rings and large deletions) were assumed to lead to clonogenic cell death. The CL yield and d were the only adjustable parameters. The model, implemented as a Monte Carlo code called BIophysical ANalysis of Cell death and chromosome Aberrations (BIANCA), provided simulated survival curves that were directly compared with experimental data on human and hamster cells exposed to photons, protons, α-particles and heavier ions including carbon and iron. d = 5 μm, independent of radiation quality, and CL yields in the range ~2-20 CLs Gy(-1) cell(-1), depending on particle type and energy, led to good agreement between simulations and data. This supports the hypothesis of a pivotal role of DNA cluster damage at sub-micrometric scale, modulated by chromosome fragment mis-rejoining at micrometric scale. To investigate the features of such critical damage, the CL yields were compared with experimental or theoretical yields of DNA fragments of different sizes, focusing on the base-pair scale (related to the so-called local clustering), the kbp scale ("regional clustering") and the Mbp scale, corresponding to chromatin loops. Interestingly, the CL yields showed better agreement with kbp fragments rather than bp fragments or Mbp fragments; this suggests that also regional clustering, in addition to other clustering levels, may play an important role, possibly due to its relationship with nucleosome organization in the chromatin fiber.

Monte Carlo simulations of the relative biological effectiveness for DNA double strand breaks from 300 MeV u(-1) carbon-ion beams

Monte Carlo simulations are used to calculate the relative biological effectiveness (RBE) of 300 MeV u(-1) carbon-ion beams at different depths in a cylindrical water phantom of 10 cm radius and 30 cm long. RBE values for the induction of DNA double strand breaks (DSB), a biological endpoint closely related to cell inactivation, are estimated for monoenergetic and energy-modulated carbon ion beams. Individual contributions to the RBE from primary ions and secondary nuclear fragments are simulated separately. These simulations are based on a multi-scale modelling approach by first applying the FLUKA (version 2011.2.17) transport code to estimate the absorbed doses and fluence energy spectra, then using the MCDS (version 3.10A) damage code for DSB yields. The approach is efficient since it separates the non-stochastic dosimetry problem from the stochastic DNA damage problem. The MCDS code predicts the major trends of the DSB yields from detailed track structure simulations. It is found that, as depth is increasing, RBE values increase slowly from the entrance depth to the plateau region and change substantially in the Bragg peak region. RBE values reach their maxima at the distal edge of the Bragg peak. Beyond this edge, contributions to RBE are entirely from nuclear fragments. Maximum RBE values at the distal edges of the Bragg peak and the spread-out Bragg peak are, respectively, 3.0 and 2.8. The present approach has the flexibility to weight RBE contributions from different DSB classes, i.e. DSB0, DSB+ and DSB++.

Effects of indirect actions and oxygen on relative biological effectiveness: estimate of DSB induction and conversion induced by gamma rays and helium ions

Clustered DNA damage other than double-strand breaks (DSBs) can be detrimental to cells and can lead to mutagenesis or cell death. In addition to DSBs induced by ionizing radiation, misrepair of non-DSB clustered damage contributes extra DSBs converted from DNA misrepair via pathways for base excision repair and nucleotide excision repair. This study aimed to quantify the relative biological effectiveness (RBE) when DSB induction and conversion from non-DSB clustered damage misrepair were used as biological endpoints. The results showed that both linear energy transfer (LET) and indirect action had a strong impact on the yields for DSB induction and conversion. RBE values for DSB induction and maximum DSB conversion of helium ions (LET = 120 keV/μm) to (60)Co gamma rays were 3.0 and 3.2, respectively. These RBE values increased to 5.8 and 5.6 in the absence of interference of indirect action initiated by addition of 2-M dimethylsulfoxide. DSB conversion was ∼1-4% of the total non-DSB damage due to gamma rays, which was lower than the 10% estimate by experimental measurement. Five to twenty percent of total non-DSB damage due to helium ions was converted into DSBs. Hence, it may be possible to increase the yields of DSBs in cancerous cells through DNA repair pathways, ultimately enhancing cell killing.

Health, growth and reproductive success of mice exposed to environmentally relevant levels of Ra-226 via drinking water over multiple generations

PURPOSE

To assess health, growth and reproductive success of mammals exposed for multiple generations to levels of radium-226 known to occur in environments surrounding uranium mines and mills in Canada.

METHODS

The study consisted of a control group and four treatment groups each containing 40 mice (20 males and 20 females) of the CBA/CaJ strain that were continuously exposed to a range of radium-226 levels via drinking water. Breeding was at 8-10 weeks of age and the study was concluded after three breeding cycles.

RESULTS

When compared to control mice, constant consumption of drinking water containing 0.012, 0.076, 0.78 and 8.0 Bq/l of radium-226 over four generations of mice did not demonstrably affect physical condition, weight, pregnancy rate, number of pups per litter, sex ratio and bodyweight gain of pups. Between generations, the observed differences in pregnancy rates that were noted in all groups, including controls, seemed to directly correlate with the weight and age of the females at breeding.

CONCLUSIONS

Based on the endpoints measured on four generations of mice, there is no indication that the consumption of radium-226 via drinking water (at activity concentrations up to 8.0 Bq/l) affects health, growth and reproductive fitness.

Nanoscale analysis of clustered DNA damage after high-LET irradiation by quantitative electron microscopy--the heavy burden to repair

Low- and high-linear energy transfer (LET) ionising radiation are effective cancer therapies, but produce structurally different forms of DNA damage. Isolated DNA damage is repaired efficiently; however, clustered lesions may be more difficult to repair, and are considered as significant biological endpoints. We investigated the formation and repair of DNA double-strand breaks (DSBs) and clustered lesions in human fibroblasts after exposure to sparsely (low-LET; delivered by photons) and densely (high-LET; delivered by carbon ions) ionising radiation. DNA repair factors (pKu70, 53BP1, γH2AX, and pXRCC1) were detected using immunogold-labelling and electron microscopy, and spatiotemporal DNA damage patterns were analysed within the nuclear ultrastructure at the nanoscale level. By labelling activated Ku-heterodimers (pKu70) the number of DSBs was determined in electron-lucent euchromatin and electron-dense heterochromatin. Directly after low-LET exposure (5 min post-irradiation), single pKu70 dimers, which reflect isolated DSBs, were randomly distributed throughout the entire nucleus with a linear dose correlation up to 30 Gy. Most euchromatic DSBs were sensed and repaired within 40 min, whereas heterochromatic DSBs were processed with slower kinetics. Essentially all DNA lesions induced by low-LET irradiation were efficiently rejoined within 24h post-irradiation. High-LET irradiation caused localised energy deposition within the particle tracks, and generated highly clustered DNA lesions with multiple DSBs in close proximity. The dimensions of these clustered lesions along the particle trajectories depended on the chromatin packing density, with huge DSB clusters predominantly localised in condensed heterochromatin. High-LET irradiation-induced clearly higher DSB yields than low-LET irradiation, with up to ∼ 500 DSBs per μm(3) track volume, and large fractions of these heterochromatic DSBs remained unrepaired. Hence, the spacing and quantity of DSBs in clustered lesions influence DNA repair efficiency, and may determine the radiobiological outcome.

Mutation induction by inhaled radon progeny modeled at the tissue level

The observable responses of living systems to ionizing radiation depend on the level of biological organization studied. Understanding the relationships between the responses characteristic of the different levels of organization is of crucial importance. The main objective of the present study is to investigate how some cellular effects of radiation manifest at the tissue level by modeling mutation induction due to chronic exposure to inhaled radon progeny. For this purpose, a mathematical model of the bronchial epithelium was elaborated to quantify cell nucleus hits and cell doses. Mutagenesis was modeled considering endogenous as well as radiation-induced DNA damages and cell cycle shortening due to cell inactivation. The model parameters describing the cellular effects of radiation are obtained from experimental data. Cell nucleus hits, cell doses, and mutation induction were computed for the activity hot spots of the large bronchi at different exposures. Results demonstrate that the mutagenic effect of densely ionizing radiation is dominated by cell cycle shortening due to cell inactivation and not by DNA damages. This suggests that radiation burdens of non-progenitor cells play a significant role in mutagenesis in case of protracted exposures to densely ionizing radiation. Mutation rate as a function of dose rate exhibits a convex shape below a threshold. This threshold indicates the exhaustion of the tissue regeneration capacity of local progenitor cells. It is suggested that progenitor cell hyperplasia occurs beyond the threshold dose rate, giving a possible explanation of the inverse dose-rate effect observed in the epidemiology of lung cancer among uranium miners.

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.

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.

Fast Monte Carlo simulation of DNA damage formed by electrons and light ions

The passage of ionizing radiation through living organisms initiates physical and chemical processes that create clusters of damaged nucleotides within one or two turns of the DNA. These clusters are widely considered an important initiating event for the induction of other biological endpoints, including cell killing and neoplastic transformation. Monte Carlo simulations of the DNA damage formation process are a useful adjunct to experiments because they provide additional information about the spatial configuration of damage within a cluster. In this paper, the fast Monte Carlo damage simulation (MCDS) algorithm is re-parameterized so that yields of double-strand breaks, single-strand breaks and sites of multiple base damage can be simulated for electrons, protons and alpha particles with kinetic energies on the order of GeV. The MCDS algorithm provides a useful, quasi-phenomenological scheme to interpolate damage yields from computationally expensive, but more detailed, track-structure simulations. The predicted characteristics of various classes of damage produced by electrons, protons and alpha particles, such as average number of lesions per DNA damage cluster and cluster length in base pairs, are presented. A study examining the effects on damage complexity of an extrinsic free radical scavenger, dimethyl sulfoxide, is also presented. The reported studies provide new information that will aid efforts to characterize the relative biological effectiveness of high-energy protons and other light ions, which are sometimes used in particle therapy for the treatment of cancer.

Simulation of light ion induced DNA damage patterns

The biophysical simulation code PARTRAC was extended by a module to handle ions heavier than alpha particles. Cross sections for ion-electron interactions were taken from He(++) ions of the same velocity and scaled by Z(eff(2))/4. Calculated linear energy transfer values, radial dose distributions and secondary electron spectra were found in agreement with experimental results. DNA damage due to irradiation of human fibroblast cells by several light ions from H to S was calculated for various energies complemented by 220 kV(p) X rays as reference radiation. With increasing linear energy transfer, the calculated total yield of double-strand breaks per dose showed saturation behaviour at about twice the value for reference radiation. When data analysis methods for experimental double-strand break yield determination were applied to the simulated DNA damage patterns, the two data sets were found in accord. The calculated patterns of DNA damage clusters were analysed on local and regional scale finding regional clusters in closer correlation to experimental cell inactivation data.

Incorporation of microdosimetric concepts into a biologically-based model of radiation carcinogenesis

The generalised state-vector model of radiation carcinogenesis (SVM) simulates radiation induced biological effects by expressing the transition rates between the various initiation and promotion stages in terms of dose rate for low and high linear energy transfer (LET) particles. In the present work, the SVM has been reformulated to incorporate single track characteristics of particles with varying LET. Transition rates of the initiation phase were expressed as functions of LET by describing the complexity and clustering of DNA double strand breaks (DSBs) and its effect on repair kinetics, while the promotion phase was reformulated based on a multi-target single-hit hypothesis. Such an approach allows the consideration of hit frequencies and the variability of the specific energy and LET spectra of radon progeny alpha particles in bronchial target cells for different exposure conditions.

DNA Double-Strand Break Misrejoining after Exposure of Primary Human Fibroblasts to CKCharacteristic X Rays, 29 kVp X Rays and60Co γ Rays

The efficiency of ionizing photon radiation for inducing mutations, chromosome aberrations, neoplastic cell transformation, and cell killing depends on the photon energy. We investigated the induction and rejoining of DNA double-strand breaks (DSBs) as possible contributors for the varying efficiencies of different photon energies. A specialized pulsed-field gel electrophoresis assay based on Southern hybridization of single Mbp genomic restriction fragments was employed to assess DSB induction and rejoining by quantifying the restriction fragment band. Unrejoined and misrejoined DSBs were determined in dose fractionation protocols using doses per fraction of 2.2 and 4.4 Gy for CK characteristic X rays, 4 and 8 Gy for 29 kVp X rays, and 5, 10 and 20 Gy for 60Co gamma rays. DSB induction by CK characteristic X rays was about twofold higher than for 60Co gamma rays, whereas 29 kVp X rays showed only marginally elevated levels of induced DSBs compared with 60Co gamma rays (a factor of 1.15). Compared with these modest variations in DSB induction, the variations in the levels of unrejoined and misrejoined DSBs were more significant. Our results suggest that differences in the fidelity of DSB rejoining together with the different efficiencies for induction of DSBs can explain the varying biological effectiveness of different photon energies.

Artemis deficiency confers a DNA double-strand break repair defect and Artemis phosphorylation status is altered by DNA damage and cell cycle progression

Mutations in the Artemis gene are causative in a subset of human severe combined immunodeficiencies (SCIDs) and Artemis-deficient cells exhibit radiation sensitivity and defective V(D)J recombination, implicating Artemis function in non-homologous end joining (NHEJ). Here we show that Artemis-deficient cells from Athabascan-speaking Native American SCID patients (SCIDA) display significantly elevated sensitivity to ionizing radiation (IR) but only a very subtle defect in DNA double-strand (DSB) break repair in contrast to the severe DSB repair defect of NHEJ-deficient cells. Primary human SCIDA fibroblasts accumulate and exhibit persistent arrest at both the G1/S and G2/M boundaries in response to IR, consistent with the presence of persistent DNA damage. Artemis protein is phosphorylated in a PI3-like kinase-dependent manner after either IR or a number of other DNA damaging treatments including etoposide, but SCIDA cells are not hypersensitive to treatment with etoposide. Inhibitor studies with various DNA damaging agents establish multiple phosphorylation states and suggest multiple kinases function in Artemis phosphorylation. We observe that Artemis phosphorylation occurs rapidly after irradiation like that of histone H2AX. However, unlike H2AX, Artemis de-phosphorylation is uncoupled from overall DNA repair and correlates instead with cell cycle progression to or through mitosis. Our results implicate a direct and non-redundant function of Artemis in the repair of a small subset of DNA double-strand breaks, possibly those with hairpin termini, which may account for the pronounced radiation sensitivity observed in Artemis-deficient cells.

Dose-Dependent Misrejoining of Radiation-Induced DNA Double-Strand Breaks in Human Fibroblasts: Experimental and Theoretical Study for High- and Low-LET Radiation

Misrejoining of DNA double-strand breaks (DSBs) was measured in human primary fibroblasts after exposure to X rays and high-LET particles (helium, nitrogen and iron) in the dose range 10-80 Gy. To measure joining of wrong DNA ends, the integrity of a 3.2-Mbp restriction fragment was analyzed directly after exposure and after 16 h of repair incubation. It was found that the misrejoining frequency for X rays was nonlinearly related to dose, with less probability of misrejoining at low doses than at high doses. The dose dependence for the high-LET particles, on the other hand, was closer to being linear, with misrejoining frequencies higher than for X rays, particularly at the lower doses. These experimental results were simulated with a Monte Carlo approach that includes a cell nucleus model with all 46 chromosomes present, combined with realistic track structure simulations to calculate the geometrical positions of all DSBs induced for each dose. The model assumes that the main determinant for misrejoining probability is the distance between two simultaneously present DSBs. With a Gaussian interaction probability function with distance, it was found that the data for both low- and high-LET radiation could be fitted with an interaction distance (sigma of the Gaussian curve) of 0.25 microm. This is half the distance previously found to best fit chromosomal aberration data in human lymphocytes using the same methods (Holley et al., Radiat. Res. 158, 568-580, 2002). The discrepancy may indicate inadequacies in the chromosome model, for example insufficient chromosomal overlap, but may also be partly due to differences between fibroblasts and lymphocytes.

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.

Distribution of breakpoints and fragment sizes in human chromosome 5 after heavy-ion bombardment

PURPOSE

To measure the location of heavy ion-induced residual breakpoints in human chromosome 5 and the size distribution of chromosome fragments involved in inter- or intrachromosomal exchanges.

MATERIALS AND METHODS

Human peripheral blood lymphocytes were exposed to 4 Gy accelerated 56Fe (iron) ions (500 MeV per nucleon (MeV n(-1)), linear energy transfer=200 keV microm(-1)). Cells were then stimulated to grow in vitro for 48 h, and chromosomes were prematurely condensed by calyculin A. Chromosome 5 was painted using high-resolution multicolour banding. The location of the observed residual breakpoints and the size of all chromosome 5 fragments involved in structural aberrations were measured using dedicated image analysis software.

RESULTS

Mapping of 283 breakpoints revealed a slight deviation from randomness, with an excess of breakpoints clustered in two small bands and an under representation of breaks at the telomeric end in the q-arm. Breakpoints per unit length were similar in p- and q-arms. The distribution of chromosome fragments has a maximum for very small fragments (< 10% of the chromosome size), indicating a severe fragmentation of chromosome 5 after heavy-ion bombardment. Only fragments < 40% of the chromosome size were involved in intrachromosomal exchanges (interstitial deletions or inversions), whereas fragments up to 75% of the whole chromosome 5 were found in interchromosomal exchanges.

CONCLUSIONS

Residual breakpoints after exposure to high-energy iron ions were not distributed randomly along chromosome 5, although the p- and q-arms displayed similar radiosensitivity. Large fragments are either restituted or misrejoined to other chromosome ends, whereas small intrachromosomal fragments can produce either inter- or intrachromosomal exchanges.

A multi-port low-fluence alpha-particle irradiator: fabrication, testing and benchmark radiobiological studies

A new multi-port irradiator, designed to facilitate the study of the effects of low fluences of alpha particles on monolayer cultures, has been developed. The irradiator consists of four individual planar (241)Am alpha-particle sources that are housed inside a helium-filled Lucite chamber. Three of the radioactive sources consist of 20 MBq of (241)Am dioxide foil. The fourth source, used to produce higher dose rates, has an activity of 500 MBq. The four sources are mounted on rotating turntables parallel to their respective 1.5-microm-thick Mylar exit windows. A stainless steel honeycomb collimator is placed between the four sources and their exit windows by a cantilever attachment to the platform of an orbital shaker that moves its table in an orbit of 2 cm. Each exit window is equipped with a beam delimiter to optimize the uniformity of the beam and with a high-precision electronic shutter. Opening and closing of the shutters is controlled with a high-precision timer. Custom-designed stainless steel Mylar-bottomed culture dishes are placed on an adapter on the shutter. The alpha particles that strike the cells have a mean energy of 2.9 MeV. The corresponding LET distribution of the particles has a mean value of 132 keV/microm. Clonogenic cell survival experiments with AG1522 human fibroblasts indicate that the RBE of the alpha particles compared to (137)Cs gamma rays is about 7.6 for this biological end point.

A model for interphase chromosomes and evaluation of radiation-induced aberrations

We have developed a theoretical model for evaluating radiation-induced chromosomal exchanges by explicitly taking into account interphase (G(0)/G(1)) chromosome structure, nuclear organization of chromosomes, the production of double-strand breaks (DSBs), and the subsequent rejoinings in a faithful or unfaithful manner. Each of the 46 chromosomes for human lymphocytes (40 chromosomes for mouse lymphocytes) is modeled as a random polymer inside a spherical volume. The chromosome spheres are packed randomly inside a spherical nucleus with an allowed overlap controlled by a parameter Omega. The rejoining of DSBs is determined by a Monte Carlo procedure using a Gaussian proximity function with an interaction range parameter sigma. Values of Omega and sigma have been found which yield calculated results of interchromosomal aberration frequencies that agree with a wide range of experimental data. Our preferred solution is one with an interaction range of 0.5 microm coupled with a relatively small overlap parameter of 0.675 microm, which more or less confirms previous estimates. We have used our model with these parameter values and with resolution or detectability limits to calculate yields of translocations and dicentrics for human lymphocytes exposed to low-LET radiation that agree with experiments in the dose range 0.09 to 4 Gy. Five different experimental data sets have been compared with the theoretical results. Essentially all of the experimental data fall between theoretical curves corresponding to resolution limits of 1 Mbp and 20 Mbp, which may reflect the fact that different investigators use different limits for sensitivity or detectability. Translocation yields for mouse lymphocytes have also been calculated and are in good agreement with experimental data from 1 cGy to 10 cGy. There is also good agreement with recent data on complex aberrations. Our model is expected to be applicable to both low- and high-LET radiation, and we include a sample prediction of the yield of interchromosomal rejoining in the dose range 0.22 Gy to 2 Gy of 1000 MeV/nucleon iron particles. This dose range corresponds to average particle traversals per nucleus ranging from 1.0 to 9.12.