RITCARD: Radiation-Induced Tracks, Chromosome Aberrations, Repair and Damage

Modelling radiobiology

Radiotherapy has played an essential role in cancer treatment for over a century, and remains one of the best-studied methods of cancer treatment. Because of its close links with the physical sciences, it has been the subject of extensive quantitative mathematical modelling, but a complete understanding of the mechanisms of radiotherapy has remained elusive. In part this is because of the complexity and range of scales involved in radiotherapy-from physical radiation interactions occurring over nanometres to evolution of patient responses over months and years. This review presents the current status and ongoing research in modelling radiotherapy responses across these scales, including basic physical mechanisms of DNA damage, the immediate biological responses this triggers, and genetic- and patient-level determinants of response. Finally, some of the major challenges in this field and potential avenues for future improvements are also discussed.

"dsbandrepair" - An updated Geant4-DNA simulation tool for evaluating the radiation-induced DNA damage and its repair

PURPOSE

Interdisciplinary scientific communities have shown large interest to achieve a mechanistic description of radiation-induced biological damage, aiming to predict biological results produced by different radiation quality exposures. Monte Carlo track-structure simulations are suitable and reliable for the study of early DNA damage induction used as input for assessing DNA damage. This study presents the most recent improvements of a Geant4-DNA simulation tool named "dsbandrepair".

METHODS

"dsbandrepair" is a Monte Carlo simulation tool based on a previous code (FullSim) that estimates the induction of early DNA single-strand breaks (SSBs) and double-strand breaks (DSBs). It uses DNA geometries generated by the DNAFabric computational tool for simulating the induction of early single-strand breaks (SSBs) and double-strand breaks (DSBs). Moreover, the new tool includes some published radiobiological models for survival fraction and un-rejoined DSB. Its application for a human fibroblast cell and human umbilical vein endothelial cell containing both heterochromatin and euchromatin was conducted. In addition, this new version offers the possibility of using the new IRT-syn method for computing the chemical stage.

RESULTS

The direct and indirect strand breaks, SSBs, DSBs, and damage complexity obtained in this work are equivalent to those obtained with the previously published simulation tool when using the same configuration in the physical and chemical stages. Simulation results on survival fraction and un-rejoined DSB are in reasonable agreement with experimental data.

CONCLUSIONS

"dsbandrepair" is a tool for simulating DNA damage and repair, benchmarked against experimental data. It has been released as an advanced example in Geant4.11.2.

Modelling Cellular Response to Ionizing Radiation: Mechanistic, Semi-Mechanistic, and Phenomenological Approaches - A Historical Perspective

Radiation research is a multidisciplinary field, and among its many branches, mathematical and computational modelers have played a significant role in advancing boundaries of knowledge. A fundamental contribution is modelling cellular response to ionizing radiation as that is the key to not only understanding how radiation can kill cancer cells, but also cause cancer and other health issues. The invention of microdosimetry in the 1950s by Harold Rossi paved the way for brilliant scientists to study the mechanism of radiation at cellular and sub-cellular scales. This paper reviews some snippets of ingenious mathematical and computational models published in microdosimetry symposium proceedings and publications of the radiation research community. Among these are simulations of radiation tracks at atomic and molecular levels using Monte Carlo methods, models of cell survival, quantification of the amount of energy required to create a single strand break, and models of DNA-damage-repair. These models can broadly be categorized into mechanistic, semi-mechanistic, and phenomenological approaches, and this review seeks to provide historical context of their development. We salute pioneers of the field and great teachers who supported and educated the younger members of the community and showed them how to build upon their work.

Galactic cosmic ray simulation at the NASA space radiation laboratory - Progress, challenges and recommendations on mixed-field effects

For missions beyond low Earth orbit to the moon or Mars, space explorers will encounter a complex radiation field composed of various ion species with a broad range of energies. Such missions pose significant radiation protection challenges that need to be solved in order to minimize exposures and associated health risks. An innovative galactic cosmic ray simulator (GCRsim) was recently developed at the NASA Space Radiation Laboratory (NSRL) at Brookhaven National Laboratory (BNL). The GCRsim technology is intended to represent major components of the space radiation environment in a ground analog laboratory setting where it can be used to improve understanding of biological risks and serve as a testbed for countermeasure development and validation. The current GCRsim consists of 33 energetic ion beams that collectively simulate the primary and secondary GCR field encountered by humans in space over the broad range of particle types, energies, and linear energy transfer (LET) of interest to health effects. A virtual workshop was held in December 2020 to assess the status of the NASA baseline GCRsim. Workshop attendees examined various aspects of simulator design, with a particular emphasis on beam selection strategies. Experimental results, modeling approaches, areas of consensus, and questions of concern were also discussed in detail. This report includes a summary of the GCRsim workshop and a description of the current status of the GCRsim. This information is important for future advancements and applications in space radiobiology.

Towards the ionizing radiation induced bond dissociation mechanism in oxygen, water, guanine and DNA fragmentation: a density functional theory simulation

The radiation-induced damages in bio-molecules are ubiquitous processes in radiotherapy and radio-biology, and critical to space projects. In this study, we present a precise quantification of the fragmentation mechanisms of deoxyribonucleic acid (DNA) and the molecules surrounding DNA such as oxygen and water under non-equilibrium conditions using the first-principle calculations based on density functional theory (DFT). Our results reveal the structural stability of DNA bases and backbone that withstand up to a combined threshold of charge and hydrogen abstraction owing to simultaneously direct and indirect ionization processes. We show the hydrogen contents of the molecules significantly control the stability in the presence of radiation. This study provides comprehensive information on the impact of the direct and indirect induced bond dissociations and DNA damage and introduces a systematic methodology for fine-tuning the input parameters necessary for the large-scale Monte Carlo simulations of radio-biological responses and mitigation of detrimental effects of ionizing radiation.

Geometrical Properties of the Nucleus and Chromosome Intermingling Are Possible Major Parameters of Chromosome Aberration Formation

Ionizing radiation causes chromosome aberrations, which are possible biomarkers to assess space radiation cancer risks. Using the Monte Carlo codes Relativistic Ion Tracks (RITRACKS) and Radiation-Induced Tracks, Chromosome Aberrations, Repair and Damage (RITCARD), we investigated how geometrical properties of the cell nucleus, irradiated with ion beams of linear energy transfer (LET) ranging from 0.22 keV/μm to 195 keV/μm, influence the yield of simple and complex exchanges. We focused on the effect of (1) nuclear volume by considering spherical nuclei of varying radii; (2) nuclear shape by considering ellipsoidal nuclei of varying thicknesses; (3) beam orientation; and (4) chromosome intermingling by constraining or not constraining chromosomes in non-overlapping domains. In general, small nuclear volumes yield a higher number of complex exchanges, as compared to larger nuclear volumes, and a higher number of simple exchanges for LET < 40 keV/μm. Nuclear flattening reduces complex exchanges for high-LET beams when irradiated along the flattened axis. The beam orientation also affects yields for ellipsoidal nuclei. Reducing chromosome intermingling decreases both simple and complex exchanges. Our results suggest that the beam orientation, the geometry of the cell nucleus, and the organization of the chromosomes within are important parameters for the formation of aberrations that must be considered to model and translate in vitro results to in vivo risks.

Nanoscale Calculation of Proton-Induced DNA Damage Using a Chromatin Geometry Model with Geant4-DNA

Monte Carlo simulations can quantify various types of DNA damage to evaluate the biological effects of ionizing radiation at the nanometer scale. This work presents a study simulating the DNA target response after proton irradiation. A chromatin fiber model and new physics constructors with the ELastic Scattering of Electrons and Positrons by neutral Atoms (ELSEPA) model were used to describe the DNA geometry and the physical stage of water radiolysis with the Geant4-DNA toolkit, respectively. Three key parameters (the energy threshold model for strand breaks, the physics model and the maximum distance to distinguish DSB clusters) of scoring DNA damage were studied to investigate the impact on the uncertainties of DNA damage. On the basis of comparison of our results with experimental data and published findings, we were able to accurately predict the yield of various types of DNA damage. Our results indicated that the difference in physics constructor can cause up to 56.4% in the DNA double-strand break (DSB) yields. The DSB yields were quite sensitive to the energy threshold for strand breaks (SB) and the maximum distance to classify the DSB clusters, which were even more than 100 times and four times than the default configurations, respectively.

The effect of non-ionizing excitations on the diffusion of ion species and inter-track correlations in FLASH ultra-high dose rate radiotherapy

Purpose. We present a microscopic mechanism that accounts for the outward burst of 'cold' ion species (IS) in a high-energy particle track due to coupling with 'hot' non-ion species (NIS). IS refers to radiolysis products of ionized molecules, whereas NIS refers to non-ionized excitations of molecules in a medium. The interaction is mediated by a quantized field of acoustic phonons, a channel that allows conversion of thermal energy of NIS to kinetic energy of IS, a flow of heat from the outer to the inner core of the track structure.Methods. We perform step-by-step Monte Carlo (MC) simulations of ionizing radiation track structures in water to score the spatial coordinates and energy depositions that form IS and NIS at atto-second time scales. We subsequently calculate the resulting temperature profiles of the tracks with MC track structure simulations and verify the results analytically using the Rutherford scattering formulation. These temperature profiles are then used as boundary conditions in a series of multi-scale atomistic molecular dynamic (MD) simulations that describe the sudden expansion and enhanced diffusive broadening of tracks initiated by the non-equilibrium spectrum of high-energy IS. We derive a stochastic coarse-grained Langevin equation of motion for IS from first-principle MD to describe the irreversible femto-second flow of thermal energy pumping from NIS to IS, mediated by quantized fields of acoustic phonons. A pair-wise Lennard-Jones potential implemented in a classical MD is then employed to validate the results calculated from the Langevin equation.Results. We demonstrate the coexistence of 'hot' NIS with 'cold' IS in the radiation track structures right after their generation. NIS, concentrated within nano-scale volumes wrapping around IS, are the main source of intensive heat-waves and the outward burst of IS due to femto-second time scale IS-NIS coupling. By comparing the transport of IS coupled to NIS with identical configurations of non-interacting IS in thermal equilibrium at room temperature, we demonstrate that the energy gain of IS due to the surrounding hot nanoscopic volumes of NIS significantly increases their effective diffusion constants. Comparing the average track separation and the time scale calculated for a deposited dose of 10 Gy and a dose rate of 40 Gy s^-1, typical values used in FLASH ultra high dose rate (UHDR) experiments, we find that the sudden expansion of tracks and ballistic transport proposed in this work strengthens the hypothesis of inter-track correlations recently introduced to interpret mitigation of the biological responses at the FLASH-UHDR (Abolfathet al2020Med. Phys.47, 6551-6561).Conclusions. The much higher diffusion constants predicted in the present model suggest higher inter-track chemical reaction rates at FLASH-UHDR, as well as lower intra-track reaction rates. This study explains why research groups relying on the current Monte Carlo frameworks have reported negligible inter-track overlaps, simply because of underestimation of the diffusion constants. We recommend incorporation of the IS-NIS coupling and heat exchange in all MC codes to enable these tool-kits to appropriately model reaction-diffusion rates at FLASH-UHDR.Novelty. To introduce a hypothetical pathway of outward burst of radiolysis products driven by highly localized thermal spikes wrapping around them and to investigate the interplay of the non-equilibrium spatio-temporal distribution of the chemical activities of diffusive high-energy particle tracks on inter-track correlations at FLASH-UHDR.

Impact of Radiation Quality on Microdosimetry and Chromosome Aberrations for High-Energy (>250 MeV/n) Ions

Studying energy deposition by space radiation at the cellular scale provides insights on health risks to astronauts. Using the Monte Carlo track structure code RITRACKS, and the chromosome aberrations code RITCARD, we performed a modeling study of single-ion energy deposition spectra and chromosome aberrations for high-energy (>250 MeV/n) ion beams with linear energy transfer (LET) varying from 0.22 to 149.2 keV/µm. The calculations were performed using cells irradiated directly by mono-energetic ion beams, and by poly-energetic beams after particle transport in a digital mouse model, representing the radiation exposure of a cell in a tissue. To discriminate events from ion tracks directly traversing the nucleus, to events from δ-electrons emitted by distant ion tracks, we categorized ion contributions to microdosimetry or chromosome aberrations into direct and indirect contributions, respectively. The ions were either ions of the mono-energetic beam or secondary ions created in the digital mouse due to interaction of the beam with tissues. For microdosimetry, the indirect contribution is largely independent of the beam LET and minimally impacted by the beam interactions in mice. In contrast, the direct contribution is strongly dependent on the beam LET and shows increased probabilities of having low and high-energy deposition events when considering beam transport. Regarding chromosome aberrations, the indirect contribution induces a small number of simple exchanges, and a negligible number of complex exchanges. The direct contribution is responsible for most simple and complex exchanges. The complex exchanges are significantly increased for some low-LET ion beams when considering beam transport.

Comparison between the results of a recently-developed biological weighting function (V79-RBE10BWF) and thein vitroclonogenic survival RBE10of other repair-competent asynchronized normoxic mammalian cell lines and ions not used for the development of the model

728 simulated microdosimetric lineal energy spectra (26 different ions between 1H and 238U, 28 energy points from 1 to 1000 MeV/n) were used in combination with a recently-developed biological weighting function (Parisi et al., 2020) and 571 published in vitro clonogenic survival curves in order to: 1) assess prediction intervals for the in silico results by deriving an empirical indication of the experimental uncertainty from the dispersion in the in vitro hamster lung fibroblast (V79) data used for the development of the biophysical model; 2) explore the possibility of modeling the relative biological effectiveness (RBE) of the 10% clonogenic survival of asynchronized normoxic repair-competent mammalian cell lines other than the one used for the development of the model (V79); 3) investigate the predictive power of the model through a comparison between in silico results and in vitro data for 10 ions not used for the development of the model. At first, different strategies for the assessment of the in silico prediction intervals were compared. The possible sources of uncertainty responsible for the dispersion in the in vitro data were also shortly reviewed. Secondly, also because of the relevant scatter in the in vitro data, no statistically-relevant differences were found between the RBE10 of the investigated different asynchronized normoxic repair-competent mammalian cell lines. The only exception (Chinese Hamster peritoneal fibroblasts, B14FAF28), is likely due to the limited dataset (all in vitro ion data were extracted from a single publication), systematic differences in the linear energy transfer (LET) calculations for the employed very-heavy ions, and the use of reference photon survival curves extracted from a different publication. Finally, the in silico predictions for the 10 ions not used for the model development were in good agreement with the corresponding in vitro data.

Track Structure Components: Characterizing Energy Deposited in Spherical Cells from Direct and Peripheral HZE Ion Hits

To understand the biological effects of radiation, it is important to determine how ionizing radiation deposits energy in micrometric targets. The energy deposited in a target located in an irradiated tissue is a function of several factors such as the radiation type and the irradiated volume size. We simulated the energy deposited by energetic ions in spherical targets of 1, 2, 4, and 8 µm radii encompassed in irradiated parallelepiped volumes of various sizes using the stochastic radiation track structure code Relativistic Ion Tracks (RITRACKS). Because cells are usually part of a tissue when they are irradiated, electrons originating from radiation tracks in neighboring volumes also contribute to energy deposition in the target. To account for this contribution, we used periodic boundary conditions in the simulations. We found that the single-ion spectra of energy deposition in targets comprises two components: the direct ion hits to the targets, which is identical in all irradiation conditions, and the contribution of hits from electrons from neighboring volumes, which depends on the irradiated volume. We also calculated an analytical expression of the indirect hit contributions using the local effect model, which showed results similar to those obtained with RITRACKS.

A Geant4-DNA Evaluation of Radiation-Induced DNA Damage on a Human Fibroblast

Simple Summary

DNA damage caused by ionizing radiation in a human fibroblast cell evaluated by the Geant4-DNA Monte Carlo toolkit is presented. A validation study using a computational geometric human DNA model was then carried out, and the calculated DNA damage as a function of particle type and energy is presented. The results of this work showed a significant improvement on past work and were consistent with recent radiobiological experimental data, such as damage yields. This work and the developed methodology could impact a broad number of research fields in which the understanding of radiation effects is crucial, such as cancer radiotherapy, space science, and medical physics.

Abstract

Accurately modeling the radiobiological mechanisms responsible for the induction of DNA damage remains a major scientific challenge, particularly for understanding the effects of low doses of ionizing radiation on living beings, such as the induction of carcinogenesis. A computational approach based on the Monte Carlo technique to simulate track structures in a biological medium is currently the most reliable method for calculating the early effects induced by ionizing radiation on DNA, the primary cellular target of such effects. The Geant4-DNA Monte Carlo toolkit can simulate not only the physical, but also the physico-chemical and chemical stages of water radiolysis. These stages can be combined with simplified geometric models of biological targets, such as DNA, to assess direct and indirect early DNA damage. In this study, DNA damage induced in a human fibroblast cell was evaluated using Geant4-DNA as a function of incident particle type (gammas, protons, and alphas) and energy. The resulting double-strand break yields as a function of linear energy transfer closely reproduced recent experimental data. Other quantities, such as fragment length distribution, scavengeable damage fraction, and time evolution of damage within an analytical repair model also supported the plausibility of predicting DNA damage using Geant4-DNA.The complete simulation chain application “molecularDNA”, an example for users of Geant4-DNA, will soon be distributed through Geant4.

Dose-effect relationships of 12C6+ ions-induced dicentric plus ring chromosomes, micronucleus and nucleoplasmic bridges in human lymphocytes in vitro

PURPOSE

The objective of this research was to explore the dose-effect relationships of dicentric plus ring (dic + r), micronucleus (MN) and nucleoplasmic bridges (NPB) induced by carbon ions in human lymphocytes.

MATERIALS AND METHODS

Venous blood samples were collected from three healthy donors. 12C6+ ions beam was used to irradiate the blood samples at the energy of 330 MeV and linear energy transfer (LET) of 50 keV/μm with a dose rate of 1 Gy/min in the spread-out Bragg peak. The irradiated doses were 0 (sham irradiation), 1, 2, 3, 4, 5 and 6 Gy. Dic + r chromosomes aberrations were scored in metaphases. The cytokinesis-block micronucleus cytome (CBMN) was conducted to analyze MN and NPB. The maximum low-dose relative biological effectiveness (RBEM) values of the induction of dic + r, MN and NPB in human lymphocytes for 12C6+ ions irradiation was calculated relative to 60Co γ-rays.

RESULTS

The frequencies of dic + r, MN and NPB showed significantly increases in a dose-depended manner after exposure to 12C6+ ions. The distributions of dic + r and MN exhibited overdispersion, while the distribution of NPB agreed with Poisson distribution at all doses. Linear-quadratic equations were established based on the frequencies of dic + r and MN. The dose-response curves of NPB frequencies followed a linear model. The derived RBEM values for dic + r, MN and NPB in human lymphocytes irradiated with 12C6+ ions were 8.07 ± 2.73, 2.69 ± 0.20 and 4.00 ± 2.69 in comparison with 60Co γ-rays.

CONCLUSION

The dose-response curves of carbon ions-induced dic + r, MN and NPB were constructed. These results could be helpful to improve radiation risk assessment and dose estimation after exposed to carbon ions irradiation.

A Mathematical Radiobiological Model (MRM) to Predict Complex DNA Damage and Cell Survival for Ionizing Particle Radiations of Varying Quality

Predicting radiobiological effects is important in different areas of basic or clinical applications using ionizing radiation (IR); for example, towards optimizing radiation protection or radiation therapy protocols. In this case, we utilized as a basis the ‘MultiScale Approach (MSA)’ model and developed an integrated mathematical radiobiological model (MRM) with several modifications and improvements. Based on this new adaptation of the MSA model, we have predicted cell-specific levels of initial complex DNA damage and cell survival for irradiation with ¹¹Β, ¹²C, ¹⁴Ν, ¹⁶Ο, ²⁰Νe, ⁴⁰Αr, ²⁸Si and ⁵⁶Fe ions by using only three input parameters (particle’s LET and two cell-specific parameters: the cross sectional area of each cell nucleus and its genome size). The model-predicted survival curves are in good agreement with the experimental ones. The particle Relative Biological Effectiveness (RBE) and Oxygen Enhancement Ratio (OER) are also calculated in a very satisfactory way. The proposed integrated MRM model (within current limitations) can be a useful tool for the assessment of radiation biological damage for ions used in hadron-beam radiation therapy or radiation protection purposes.

A review of simulation codes and approaches for radiation chemistry

Historically, the field of radiation chemistry began shortly after the discovery of radioactivity, and its development has been closely related to discoveries in other related fields such as radiation and nuclear physics. Radiolysis of water and radiation chemistry have been very important in elucidating how radiation affects living matter and how it induces DNA damage. Nowadays, we recognize the importance of chemistry to understanding the effects of radiation on cells; however, it took several decades to obtain this insight, and much is still unknown. The radiolysis of water and aqueous solutions have been the subject of much experimental and theoretical research for many decades. One important concept closely related to radiation chemistry is radiation track structure. Track structure results from early physical and physicochemical events that lead to a highly non-homogenous distribution of radiolytic species. Because ionizing radiation creates unstable species that are distributed non-homogenously, the use of conventional reaction kinetics methods does not describe this chemistry well. In recent years, several methods have been developed for simulating radiation chemistry. In this review, we give a brief history of the field and the development of the simulation codes. We review the current methods used to simulate radiolysis of water and radiation chemistry, and we describe several radiation chemistry codes and their applications.

Determination of Chromosome Aberrations in Human Fibroblasts Irradiated by Mixed Fields Generated with Shielding

To better study biological effects of space radiation using ground-based facilities, the NASA Space Radiation Laboratory (NSRL) at the Brookhaven National Laboratory has been upgraded to rapidly switch ions and energies. This has allowed investigators to design irradiation protocols comprising a mixture of ions and energies more indicative of the galactic cosmic ray (GCR) environment. Despite these advancements, beam selection and delivery schemes should be optimized against facility and experimental constraints and validated to ensure such irradiations are a suitable representation of the space environment. Importantly, since experiments are time consuming and expensive, models capable of predicting biological outcomes over a range of irradiation conditions (single ion, sequential multi ion or mixed fields) are needed to support such efforts. In this work, human fibroblasts were placed behind 20 g/cm2 aluminum and 10.345 g/cm2 polyethylene and irradiated separately by 344 MeV hydrogen, 344 MeV/n helium, 450 MeV/n oxygen and 950 MeV/n iron ions at various doses. The fluorescence in situ hybridization (FISH) whole chromosome painting technique was then used to assess the cells for chromosome aberrations (CAs), notably simple exchanges. A multi-scale modeling approach was also developed to predict the formation of chromosome aberrations in these experiments. The Geant4 simulation toolkit was used to determine the spectra of particles and energies produced by interactions between the incident beams and shielding. The simulated mixed field generated by shielding was then transferred into the track structure code, RITRACKS (relativistic ion tracks), to generate three-dimensional (3D) voxelized dose maps at the nanometer scale. Finally, these voxel dose maps were input into the new damage and repair model, RITCARD (radiation-induced tracks, chromosome aberrations, repair and damage), to predict the formation of various CAs. The multi-scale model described herein is a significant advancement for the computational tools used to predict biological outcomes in cells exposed to highly complex, mixed ion fields related to the GCR environment. Results show that the simulation and experimental data are in good agreement for the complex radiation fields generated by all ions incident on shielding for most data points. The differences between model predictions and measurements are discussed. Although improvements are needed, the model extends current capabilities for evaluating beam selection and delivery schemes at the NSRL ground-based GCR simulator and for informing NASA risk projection models in the future.

Simultaneous Exposure of Cultured Human Lymphoblastic Cells to Simulated Microgravity and Radiation Increases Chromosome Aberrations

During space travel, humans are continuously exposed to two major environmental stresses, microgravity (μG) and space radiation. One of the fundamental questions is whether the two stressors are interactive. For over half a century, many studies were carried out in space, as well as using devices that simulated μG on the ground to investigate gravity effects on cells and organisms, and we have gained insights into how living organisms respond to μG. However, our knowledge on how to assess and manage human health risks in long-term mission to the Moon or Mars is drastically limited. For example, little information is available on how cells respond to simultaneous exposure to space radiation and μG. In this study, we analyzed the frequencies of chromosome aberrations (CA) in cultured human lymphoblastic TK6 cells exposed to X-ray or carbon ion under the simulated μG conditions. A higher frequency of both simple and complex types of CA were observed in cells exposed to radiation and μG simultaneously compared to CA frequency in cells exposed to radiation only. Our study shows that the dose response data on space radiation obtained at the 1G condition could lead to the underestimation of astronauts' potential risk for health deterioration, including cancer. This study also emphasizes the importance of obtaining data on the molecular and cellular responses to irradiation under μG conditions.

Ionizing Radiation and Complex DNA Damage: Quantifying the Radiobiological Damage Using Monte Carlo Simulations

Ionizing radiation is a common tool in medical procedures. Monte Carlo (MC) techniques are widely used when dosimetry is the matter of investigation. The scientific community has invested, over the last 20 years, a lot of effort into improving the knowledge of radiation biology. The present article aims to summarize the understanding of the field of DNA damage response (DDR) to ionizing radiation by providing an overview on MC simulation studies that try to explain several aspects of radiation biology. The need for accurate techniques for the quantification of DNA damage is crucial, as it becomes a clinical need to evaluate the outcome of various applications including both low- and high-energy radiation medical procedures. Understanding DNA repair processes would improve radiation therapy procedures. Monte Carlo simulations are a promising tool in radiobiology studies, as there are clear prospects for more advanced tools that could be used in multidisciplinary studies, in the fields of physics, medicine, biology and chemistry. Still, lot of effort is needed to evolve MC simulation tools and apply them in multiscale studies starting from small DNA segments and reaching a population of cells.

A Bi-Exponential Repair Algorithm for Radiation-Induced Double-Strand Breaks: Application to Simulation of Chromosome Aberrations

Background: Radiation induces DNA double-strand breaks (DSBs), and chromosome aberrations (CA) form during the DSBs repair process. Several methods have been used to model the repair kinetics of DSBs including the bi-exponential model, i.e., N(t) = N₁exp(−t/τ₁) + N₂exp(−t/τ₂), where N(t) is the number of breaks at time t, and N₁, N₂, τ₁ and τ₂ are parameters. This bi-exponential fit for DSB decay suggests that some breaks are repaired rapidly and other, more complex breaks, take longer to repair. Methods: The bi-exponential repair kinetics model is implemented into a recent simulation code called RITCARD (Radiation Induced Tracks, Chromosome Aberrations, Repair, and Damage). RITCARD simulates the geometric configuration of human chromosomes, radiation-induced breaks, their repair, and the creation of various categories of CAs. The bi-exponential repair relies on a computational algorithm that is shown to be mathematically exact. To categorize breaks as complex or simple, a threshold for the local (voxel) dose was used. Results: The main findings are: i) the curves for the kinetics of restitution of DSBs are mostly independent of dose; ii) the fraction of unrepaired breaks increases with the linear energy transfer (LET) of the incident radiation; iii) the simulated dose–response curves for simple reciprocal chromosome exchanges that are linear-quadratic; iv) the alpha coefficient of the dose–response curve peaks at about 100 keV/µm.