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

Evaluation of the Yield of DNA Double-Strand Breaks for Carbon Ions Using Monte Carlo Simulation and DNA Fragment Distribution

PURPOSE

The aim of this work was to provide a method to evaluate the yield of DNA double-strand breaks (DSBs) for carbon ions, overcoming the bias in existing methods due to the nonrandom distribution of DSBs.

METHODS AND MATERIALS

A previously established biophysical program based on the radiation track structure and a multilevel chromosome model was used to simulate DNA damage induced by x-rays and carbon ions. The fraction of activity retained (FAR) as a function of absorbed dose or particle fluence was obtained by counting the fraction of DNA fragments larger than 6 Mbp. Simulated FAR curves for the 250 kV x-rays and carbon ions at various energies were compared with measurements using constant-field gel electrophoresis. The doses or fluences at the FAR of 0.7 based on linear interpolation were used to estimate the simulation error for the production of DSBs.

RESULTS

The relative difference of doses at the FAR of 0.7 between simulation and experiment was -8.5% for the 250 kV x-rays. The relative differences of fluences at the FAR of 0.7 between simulations and experiments were -17.5%, -42.2%, -18.2%, -3.1%, 10.8%, and -14.5% for the 34, 65, 130, 217, 2232, and 3132 MeV carbon ions, respectively. In comparison, the measurement uncertainty was about 20%. Carbon ions produced remarkably more DSBs and DSB clusters per unit dose than x-rays. The yield of DSBs for carbon ions, ranging from 10 to 16 Gbp-1Gy-1, increased with linear energy transfer (LET) but plateaued in the high-LET end. The yield of DSB clusters first increased and then decreased with LET. This pattern was similar to the relative biological effectiveness for cell survival for heavy ions.

CONCLUSIONS

The estimated yields of DSBs for carbon ions increased from 10 Gbp-1Gy-1 in the low-LET end to 16 Gbp-1Gy-1 in the high-LET end with 20% uncertainty.

Fast neutron energy based modelling of biological effectiveness with implications for proton and ion beams

A practical neutron energy dependent RBE model has been developed, based on the relationship between a mono-energetic neutron energy and its likely recoil proton energy. Essentially, the linear energy transfer (LET) values of the most appropriate recoil proton energies are then used to modify the linear quadratic model radiosensitivities (α and β) from their reference LET radiation values to provide the RBE estimates. Experimental neutron studies published by Hall (including some mono-energetic beams ranging from 0.2 to 15 MeV), Broerse, Berry, and data from the Clatterbridge and Detroit clinical neutron beams, which all contain some data from a spectrum of neutron energies, are used to derive single effective neutron energies (NE_eff) for each spectral beam. These energies yield a recoil proton spectrum, but with an effective mean proton energy (being around 50% of NE_eff). The fractional increase in LET is given by the recoil proton LET divided by the proton (LET_U) value which provides the highest RBE. This ratio is then used to determine the change in the linear-quadratic model α and β parameters, from those of the reference radiation, to estimate the RBE. The predicted proton recoil RBE is then reasonably close to the experimental neutron RBE values found when taking into account the variation inherent in biological experiments. The work has some important consequences. The data of Hall et al (1975 Radiat. Res. 64 245-55) shows that the highest RBE values are found with neutron energies around 0.3-0.4 MeV, but this energy cannot possibly generate recoil proton energies which are higher, as necessary for a 0.68 MeV proton with a 30.5 keV μm^-1 LET_U (the LET value which provides the maximum obtainable RBE for a specified ion). For 0.4 MeV neutrons with proton recoil energies of around 0.2 MeV, the latter have a LET of around 62.88 keV μm^-1. This could have an impact on proton beam RBE modelling. However, this is compensated by finding that the maximum radiosensitivity for mono-energetic neutrons was around 1.7 times larger than previously suggested from experimental ion beam studies, probably due to the necessary spreading out of Bragg peaks for ion beam experimental purposes, sampling errors and particle range considerations. This semi-empirical model can be used with minimal computer support and could have applications in ionic beams and in radioprotection.

Iron Ladies - How Desiccated Asexual Rotifer Adineta vaga Deal With X-Rays and Heavy Ions?

Space exposure experiments from the last 15 years have unexpectedly shown that several terrestrial organisms, including some multi-cellular species, are able to survive in open space without protection. The robustness of bdelloid rotifers suggests that these tiny creatures can possibly be added to the still restricted list of animals that can deal with the exposure to harsh condition of space. Bdelloids are one of the smallest animals on Earth. Living all over the world, mostly in semi-terrestrial environments, they appear to be extremely stress tolerant. Their desiccation tolerance at any stage of their life cycle is known to confer tolerance to a variety of stresses including high doses of radiation and freezing. In addition, they constitute a major scandal in evolutionary biology due to the putative absence of sexual reproduction for at least 60 million years. Adineta vaga, with its unique characteristics and a draft genome available, was selected by ESA (European Space Agency) as a model system to study extreme resistance of organisms exposed to space environment. In this manuscript, we documented the resistance of desiccated A. vaga individuals exposed to increasing doses of X-ray, protons and Fe ions. Consequences of exposure to different sources of radiation were investigated in regard to the cellular type including somatic (survival assay) and germinal cells (fertility assay). Then, the capacity of A. vaga individuals to repair DNA DSB induced by different source of radiation was investigated. Bdelloid rotifers represent a promising model in order to investigate damage induced by high or low LET radiation. The possibility of exposure both on hydrated or desiccated specimens may help to decipher contribution of direct and indirect radiation damage on biological processes. Results achieved through this study consolidate our knowledge about the radioresistance of A. vaga and improve our capacity to compare extreme resistance against radiation among living organisms including metazoan.

Predicting DNA damage foci and their experimental readout with 2D microscopy: a unified approach applied to photon and neutron exposures

The consideration of how a given technique affects results of experimental measurements is a must to achieve correct data interpretation. This might be challenging when it comes to measurements on biological systems, where it is unrealistic to have full control (e.g. through a software replica) of all steps in the measurement chain. In this work we address how the effectiveness of different radiation qualities in inducing biological damage can be assessed measuring DNA damage foci yields, only provided that artefacts related to the scoring technique are adequately considered. To this aim, we developed a unified stochastic modelling approach that, starting from radiation tracks, predicts both the induction, spatial distribution and complexity of DNA damage, and the experimental readout of foci when immunocytochemistry coupled to 2D fluorescence microscopy is used. The approach is used to interpret γ-H2AX data for photon and neutron exposures. When foci are reconstructed in the whole cell nucleus, we obtain information on damage characteristics "behind" experimental observations, as the average damage content of a focus. We reproduce how the detection technique affects experimental findings, e.g. contributing to the saturation of foci yields scored at 30 minutes after exposure with increasing dose and to the lack of dose dependence for yields at 24 hours.

WHAT ROLES FOR TRACK-STRUCTURE AND MICRODOSIMETRY IN THE ERA OF -omics AND SYSTEMS BIOLOGY?

Ionizing radiation is a peculiar perturbation when it comes to damage to biological systems: it proceeds through discrete energy depositions, over a short temporal scale and a spatial scale critical for subcellular targets as DNA, whose damage complexity determines the outcome of the exposure. This lies at the basis of the success of track structure (and nanodosimetry) and microdosimetry in radiation biology. However, such reductionist approaches cannot account for the complex network of interactions regulating the overall response of the system to radiation, particularly when effects are manifest at the supracellular level and involve long times. Systems radiation biology is increasingly gaining ground, but the gap between reductionist and holistic approaches is becoming larger. This paper presents considerations on what roles track structure and microdosimetry can have in the attempt to fill this gap, and on how they can be further exploited to interpret radiobiological data and inform systemic approaches.

A mechanistic relative biological effectiveness model-based biological dose optimization for charged particle radiobiology studies

In charged particle therapy, the objective is to exploit both the physical and radiobiological advantages of charged particles to improve the therapeutic index. Use of the beam scanning technique provides the flexibility to implement biological dose optimized intensity-modulated ion therapy (IMIT). An easy-to-implement algorithm was developed in the current study to rapidly generate a uniform biological dose distribution, namely the product of physical dose and the relative biological effectiveness (RBE), within the target volume using scanned ion beams for charged particle radiobiological studies. Protons, helium ions and carbon ions were selected to demonstrate the feasibility and flexibility of our method. The general-purpose Monte Carlo simulation toolkit Geant4 was used for particle tracking and generation of physical and radiobiological data needed for later dose optimizations. The dose optimization algorithm was developed using the Python (version 3) programming language. A constant RBE-weighted dose (RWD) spread-out Bragg peak (SOBP) in a water phantom was selected as the desired target dose distribution to demonstrate the applicability of the optimization algorithm. The mechanistic repair-misrepair-fixation (RMF) model was incorporated into the Monte Carlo particle tracking to generate radiobiological parameters and was used to predict the RBE of cell survival in the iterative process of the biological dose optimization for the three selected ions. The post-optimization generated beam delivery strategy can be used in radiation biology experiments to obtain radiobiological data to further validate and improve the accuracy of the RBE model. This biological dose optimization algorithm developed for radiobiology studies could potentially be extended to implement biologically optimized IMIT plans for patients.

DNA double strand break (DSB) induction and cell survival in iodine-enhanced computed tomography (CT)

A multi-scale Monte Carlo model is proposed to assess the dosimetric and biological impact of iodine-based contrast agents commonly used in computed tomography. As presented, the model integrates the general purpose MCNP6 code system for larger-scale radiation transport and dose assessment with the Monte Carlo damage simulation to determine the sub-cellular characteristics and spatial distribution of initial DNA damage. The repair-misrepair-fixation model is then used to relate DNA double strand break (DSB) induction to reproductive cell death. Comparisons of measured and modeled changes in reproductive cell survival for ultrasoft characteristic k-shell x-rays (0.25-4.55 keV) up to orthovoltage (200-500 kVp) x-rays indicate that the relative biological effectiveness (RBE) for DSB induction is within a few percent of the RBE for cell survival. Because of the very short range of secondary electrons produced by low energy x-ray interactions with contrast agents, the concentration and subcellular distribution of iodine within and near cellular targets have a significant impact on the estimated absorbed dose and number of DSB produced in the cell nucleus. For some plausible models of the cell-level distribution of contrast agent, the model predicts an increase in RBE-weighted dose (RWD) for the endpoint of DSB induction of 1.22-1.40 for a 5-10 mg ml^-1 iodine concentration in blood compared to an RWD increase of 1.07  ±  0.19 from a recent clinical trial. The modeled RWD of 2.58  ±  0.03 is also in good agreement with the measured RWD of 2.3  ±  0.5 for an iodine concentration of 50 mg ml^-1 relative to no iodine. The good agreement between modeled and measured DSB and cell survival estimates provides some confidence that the presented model can be used to accurately assess biological dose for other concentrations of the same or different contrast agents.

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.

Rapid MCNP simulation of DNA double strand break (DSB) relative biological effectiveness (RBE) for photons, neutrons, and light ions

To account for particle interactions in the extracellular (physical) environment, information from the cell-level Monte Carlo damage simulation (MCDS) for DNA double strand break (DSB) induction has been integrated into the general purpose Monte Carlo N-particle (MCNP) radiation transport code system. The effort to integrate these models is motivated by the need for a computationally efficient model to accurately predict particle relative biological effectiveness (RBE) in cell cultures and in vivo. To illustrate the approach and highlight the impact of the larger scale physical environment (e.g. establishing charged particle equilibrium), we examined the RBE for DSB induction (RBEDSB) of x-rays, (137)Cs γ-rays, neutrons and light ions relative to γ-rays from (60)Co in monolayer cell cultures at various depths in water. Under normoxic conditions, we found that (137)Cs γ-rays are about 1.7% more effective at creating DSB than γ-rays from (60)Co (RBEDSB  =  1.017) whereas 60-250 kV x-rays are 1.1 to 1.25 times more efficient at creating DSB than (60)Co. Under anoxic conditions, kV x-rays may have an RBEDSB up to 1.51 times as large as (60)Co γ-rays. Fission neutrons passing through monolayer cell cultures have an RBEDSB that ranges from 2.6 to 3.0 in normoxic cells, but may be as large as 9.93 for anoxic cells. For proton pencil beams, Monte Carlo simulations suggest an RBEDSB of about 1.2 at the tip of the Bragg peak and up to 1.6 a few mm beyond the Bragg peak. Bragg peak RBEDSB increases with decreasing oxygen concentration, which may create opportunities to apply proton dose painting to help address tumor hypoxia. Modeling of the particle RBE for DSB induction across multiple physical and biological scales has the potential to aid in the interpretation of laboratory experiments and provide useful information to advance the safety and effectiveness of hadron therapy in the treatment of cancer.

Stimulation of intercellular induction of apoptosis in transformed cells at very low doses of ionising radiation: spatial and temporal features

The ultimate response of a cell or tissue to radiation is dependent in part on intercellular signalling. This becomes increasingly important at low doses, or at low dose rates, associated with typical human exposures. In order to help characterise the underlying mechanism of intercellular signalling, and how they are perturbed following exposure to ionising radiation, a previously well-defined model system of intercellular induction of apoptosis (IIA) (Portess et al. 2007, Cancer Res. 67, 1246-1253) was adopted. The aim of the present work is to evaluate the signalling mechanisms underpinning this process through exploring the variables that can affect the IIA, i.e. dose, time and space.

Energy dependence of the complexity of DNA damage induced by carbon ions

To assess the complexity of DNA damage induced by carbon ions as a function of their energy and LET, 2-Gy irradiations by 100 keV u(-1)-400 MeV u(-1) carbon ions were investigated using the PARTRAC code. The total number of fragments and the yield of fragments of <30 bp were calculated. The authors found a particularly important contribution of DNA fragmentation in the range of <1 kbp for specific energies of <6 MeV u(-1). They also considered the effect of different specific energies with the same LET, i.e. before and after the Bragg peak. As a first step towards a full characterisation of secondary particle production from carbon ions interacting with tissue, a comparison between DNA-damage induction by primary carbon ions and alpha particles resulting from carbon break-up is presented, for specific energies of >1 MeV u(-1).

Investigation of radiation-induced multilayered signalling response of the inflammatory pathway

Ionising radiation exposure of cells might induce the perturbation of cell functions and, in particular, the activation or inhibition of several important pathways. This perturbation can cause the deregulation of both intra- and extra-cellular signalling cascades (such as the inflammatory pathway) and alter not only the behaviour of directly exposed cells but also the neighbouring non-irradiated ones, through the so-called bystander effect. The aim of the present work was to investigate the complex nonlinear interactions between the inflammatory pathway and other strictly interlaced signalling pathways, such as Erk1/2 and Akt/PKB, focusing on the radiation-induced perturbation of such pathways in the dose range of 0-2 Gy. The results show how radiation affects these interconnected pathways and how confounding factors, such as the change of culture medium, can hide radiation-induced perturbations.

Modelling proton bunches focussed to submicrometre scales: low-LET radiation damage in high-LET-like spatial structure

Microbeam experiments approximating high-LET tracks by bunches of lower-LET particles focussed to submicrometre scales (Schmid et al. 2012, Phys. Med. Biol. 57, 5889) provide an unprecedented benchmark for models of biological effects of radiation. PARTRAC track structure-based Monte Carlo simulations have verified that focussed 20 MeV proton bunches resemble the radial dose distributions of single 55 MeV carbon ions as used in the experiments. However, the predicted yields of double-strand break and short (<1 kbp) DNA fragments by focussed protons correspond to homogeneous proton irradiation and are much smaller than for carbon tracks. The calculated yields of dicentrics overestimate the effect of focussing but reproduce the fourfold difference between carbon ions and homogeneously distributed protons. The extent to which focussed low-LET particles approximate high-LET radiation is limited by the achievable focussing: submicrometre focussing of proton bunches cannot reproduce local nanometre clustering, i.e. DNA damage complexity characteristic of high-LET radiation.

Sub-micrometer 20MeV protons or 45MeV 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μm(2) 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.

Use of the γ-H2AX assay to investigate DNA repair dynamics following multiple radiation exposures

Radiation therapy is one of the most common and effective strategies used to treat cancer. The irradiation is usually performed with a fractionated scheme, where the dose required to kill tumour cells is given in several sessions, spaced by specific time intervals, to allow healthy tissue recovery. In this work, we examined the DNA repair dynamics of cells exposed to radiation delivered in fractions, by assessing the response of histone-2AX (H2AX) phosphorylation (γ-H2AX), a marker of DNA double strand breaks. γ-H2AX foci induction and disappearance were monitored following split dose irradiation experiments in which time interval between exposure and dose were varied. Experimental data have been coupled to an analytical theoretical model, in order to quantify key parameters involved in the foci induction process. Induction of γ-H2AX foci was found to be affected by the initial radiation exposure with a smaller number of foci induced by subsequent exposures. This was compared to chromatin relaxation and cell survival. The time needed for full recovery of γ-H2AX foci induction was quantified (12 hours) and the 1:1 relationship between radiation induced DNA double strand breaks and foci numbers was critically assessed in the multiple irradiation scenarios.

Use of the γ-H2AX assay to investigate DNA repair dynamics following multiple radiation exposures

Radiation therapy is one of the most common and effective strategies used to treat cancer. The irradiation is usually performed with a fractionated scheme, where the dose required to kill tumour cells is given in several sessions, spaced by specific time intervals, to allow healthy tissue recovery. In this work, we examined the DNA repair dynamics of cells exposed to radiation delivered in fractions, by assessing the response of histone-2AX (H2AX) phosphorylation (γ-H2AX), a marker of DNA double strand breaks. γ-H2AX foci induction and disappearance were monitored following split dose irradiation experiments in which time interval between exposure and dose were varied. Experimental data have been coupled to an analytical theoretical model, in order to quantify key parameters involved in the foci induction process. Induction of γ-H2AX foci was found to be affected by the initial radiation exposure with a smaller number of foci induced by subsequent exposures. This was compared to chromatin relaxation and cell survival. The time needed for full recovery of γ-H2AX foci induction was quantified (12 hours) and the 1:1 relationship between radiation induced DNA double strand breaks and foci numbers was critically assessed in the multiple irradiation scenarios.

Ionizing radiation and genetic risks. XVII. Formation mechanisms underlying naturally occurring DNA deletions in the human genome and their potential relevance for bridging the gap between induced DNA double-strand breaks and deletions in irradiated germ cells

While much is known about radiation-induced DNA double-strand breaks (DSBs) and their repair, the question of how deletions of different sizes arise as a result of the processing of DSBs by the cell's repair systems has not been fully answered. In order to bridge this gap between DSBs and deletions, we critically reviewed published data on mechanisms pertaining to: (a) repair of DNA DSBs (from basic studies in this area); (b) formation of naturally occurring structural variation (SV) - especially of deletions - in the human genome (from genomic studies) and (c) radiation-induced mutations and structural chromosomal aberrations in mammalian somatic cells (from radiation mutagenesis and radiation cytogenetic studies). The specific aim was to assess the relative importance of the postulated mechanisms in generating deletions in the human genome and examine whether empirical data on radiation-induced deletions in mouse germ cells are consistent with predictions of these mechanisms. The mechanisms include (a) NHEJ, a DSB repair process that does not require any homology and which functions in all stages of the cell cycle (and is of particular relevance in G0/G1); (b) MMEJ, also a DSB repair process but which requires microhomology and which presumably functions in all cell cycle stages; (c) NAHR, a recombination-based DSB repair mechanism which operates in prophase I of meiosis in germ cells; (d) MMBIR, a microhomology-mediated, replication-based mechanism which operates in the S phase of the cell cycle, and (e) strand slippage during replication (involved in the origin of small insertions and deletions (INDELs). Our analysis permits the inference that, between them, these five mechanisms can explain nearly all naturally occurring deletions of different sizes identified in the human genome, NAHR and MMBIR being potentially more versatile in this regard. With respect to radiation-induced deletions, the basic studies suggest that those arising as a result of the operation of NHEJ/MMEJ processes, as currently formulated, are expected to be relatively small. However, data on induced mutations in mouse spermatogonial stem cells (irradiation in G0/G1 phase of the cell cycle and DSB repair presumed to be via NHEJ predominantly) show that most are associated with deletions of different sizes, some in the megabase range. There is thus a 'discrepancy' between what the basic studies suggest and the empirical observations in mutagenesis studies. This discrepancy, however, is only an apparent but not a real one. It can be resolved by considering the issue of deletions in the broader context of and in conjunction with the organization of chromatin in chromosomes and nuclear architecture, the conceptual framework for which already exists in studies carried out during the past fifteen years or so. In this paper, we specifically hypothesize that repair of DSBs induced in chromatin loops may offer a basis to explain the induction of deletions of different sizes and suggest an approach to test the hypothesis. We emphasize that the bridging of the gap between induced DSB and resulting deletions of different sizes is critical for current efforts in computational modeling of genetic risks.