Modelling DNA damage-repair and beyond

The paper presents a review of mechanistic modelling studies of DNA damage and DNA repair, and consequences to follow in mammalian cell nucleus. We hypothesize DNA deletions are consequences of repair of double strand breaks leading to the modifications of genome that play crucial role in long term development of genetic inheritance and diseases. The aim of the paper is to review 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 damaged human genome from endogenous and exogenous events. The model of the cell nucleus presented enables simulation of DNA damage at molecular level identifying the spectrum of damage induced in all chromosomal territories and loops. Our mechanistic modelling of DNA repair for double stand breaks (DSB), single strand breaks (SSB) and base damage (BD), shows the complexity of DNA damage is responsible for the longer repair times and the reason for the biphasic feature of mammalian cells repair curves. In the absence of experimentally determined data, the mechanistic model of repair predicts the in vivo rate constants for the proteins involved in the repair of DSB, SSB, and of BD.

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.

Distinct characteristics of the DNA damage response in mammalian oocytes

DNA damage is a critical threat that poses significant challenges to all cells. To address this issue, cells have evolved a sophisticated molecular and cellular process known as the DNA damage response (DDR). Among the various cell types, mammalian oocytes, which remain dormant in the ovary for extended periods, are particularly susceptible to DNA damage. The occurrence of DNA damage in oocytes can result in genetic abnormalities, potentially leading to infertility, birth defects, and even abortion. Therefore, understanding how oocytes detect and repair DNA damage is of paramount importance in maintaining oocyte quality and preserving fertility. Although the fundamental concept of the DDR is conserved across various cell types, an emerging body of evidence reveals striking distinctions in the DDR between mammalian oocytes and somatic cells. In this review, we highlight the distinctive characteristics of the DDR in oocytes and discuss the clinical implications of DNA damage in oocytes.

Modeling of DNA Damage Repair and Cell Response in Relation to p53 System Exposed to Ionizing Radiation

Repair of DNA damage induced by ionizing radiation plays an important role in the cell response to ionizing radiation. Radiation-induced DNA damage also activates the p53 system, which determines the fate of cells. The kinetics of repair, which is affected by the cell itself and the complexity of DNA damage, influences the cell response and fate via affecting the p53 system. To mechanistically study the influences of the cell response to different LET radiations, we introduce a new repair module and a p53 system model with NASIC, a Monte Carlo track structure code. The factors determining the kinetics of the double-strand break (DSB) repair are modeled, including the chromosome environment and complexity of DSB. The kinetics of DSB repair is modeled considering the resection-dependent and resection-independent compartments. The p53 system is modeled by simulating the interactions among genes and proteins. With this model, the cell responses to low- and high-LET irradiation are simulated, respectively. It is found that the kinetics of DSB repair greatly affects the cell fate and later biological effects. A large number of DSBs and a slow repair process lead to severe biological consequences. High-LET radiation induces more complex DSBs, which can be repaired by slow processes, subsequently resulting in a longer cycle arrest and, furthermore, apoptosis and more secreting of TGFβ. The Monte Carlo track structure simulation with a more realistic repair module and the p53 system model developed in this study can expand the functions of the NASIC code in simulating mechanical radiobiological effects.

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.

A Mechanistic DNA Repair and Survival Model (Medras): Applications to Intrinsic Radiosensitivity, Relative Biological Effectiveness and Dose-Rate

Variations in the intrinsic radiosensitivity of different cells to ionizing radiation is now widely believed to be a significant driver in differences in response to radiotherapy. While the mechanisms of radiosensitivity have been extensively studied in the laboratory, there are a lack of models which integrate this knowledge into a predictive framework. This paper presents an overview of the Medras model, which has been developed to provide a mechanistic framework in which different radiation responses can be modelled and individual responses predicted. This model simulates the repair of radiation-induced DNA damage, incorporating the overall kinetics of repair and its fidelity, to predict a range of biological endpoints including residual DNA damage, mutation, chromosome aberration, and cell death. Validation of this model against a range of exposure types is presented, including considerations of varying radiation qualities and dose-rates. This approach has the potential to inform new tools to deliver mechanistic predictions of radiation sensitivity, and support future developments in treatment personalization.

Multi-scaled Monte Carlo calculation for radon-induced cellular damage in the bronchial airway epithelium

Radon is a leading cause of lung cancer in indoor public and mining workers. Inhaled radon progeny releases alpha particles, which can damage cells in the airway epithelium. The extent and complexity of cellular damage vary depending on the alpha particle's kinetic energy and cell characteristics. We developed a framework to quantitate the cellular damage on the nanometer and micrometer scales at different intensities of exposure to radon progenies Po-218 and Po-214. Energy depositions along the tracks of alpha particles that were slowing down were simulated on a nanometer scale using the Monte Carlo code Geant4-DNA. The nano-scaled track histories in a 5 μm radius and 1 μm-thick cylindrical volume were integrated into the tracking scheme of alpha trajectories in a micron-scale bronchial epithelium segment in the user-written SNU-CDS program. Damage distribution in cellular DNA was estimated for six cell types in the epithelium. Deep-sited cell nuclei in the epithelium would have less chance of being hit, but DNA damage from a single hit would be more serious, because low-energy alpha particles of high LET would hit the nuclei. The greater damage in deep-sited nuclei was due to the 7.69 MeV alpha particles emitted from Po-214. From daily work under 1 WL of radon concentration, basal cells would respond with the highest portion of complex DSBs among the suspected progenitor cells in the most exposed regions of the lung epithelium.

Microdosimetric investigation of the radiation quality of low-medium energy electrons using Geant4-DNA

The increasing clinical use of low-energy photon and electron sources (below few tens of keV) has raised concerns on the adequacy of the existing approximation of an energy-independent radiobiological effectiveness. In this work, the variation of the quality factor (Q) and relative biological effectiveness (RBE) of electrons over the low-medium energy range (0.1 keV-1 MeV) is examined using several microdosimetry-based Monte Carlo methodologies with input data obtained from Geant4-DNA track-structure simulations. The sensitivity of the results to the different methodologies, Geant4-DNA physics models, and target sizes is examined. Calculations of Q and RBE are based on the ICRU Report 40 recommendations, the Kellerer-Hahn approximation, the site version of the theory of dual radiation action (TDRA), the microdosimetric kinetic model (MKM) of cell survival, and the calculated yield of DNA double strand breaks (DSB). The stochastic energy deposition spectra needed as input in the above approaches have been calculated for nanometer spherical volumes using the different electron physics models of Geant4-DNA. Results are normalized at 100 keV electrons which is here considered the reference radiation. It is shown that in the energy range ~50 keV-1 MeV, the calculated Q and RBE are approximately unity (to within 1-2%) irrespective of the methodology, Geant4-DNA physics model, and target size. At lower energies, Q and RBE become energy-dependent reaching a maximum value of ~1.5-2.5 between ~200 and 700 eV. The detailed variation of Q and RBE at low energies depends mostly upon the adopted methodology and target size, and less so upon the Geant4-DNA physics model. Overall, the DSB yield predicts the highest RBE values (with RBE_max≈2.5) whereas the MKM the lowest RBE values (with RBE_max≈1.5). The ICRU Report 40, Kellerer-Hahn, and TDRA methods are in excellent agreement (to within 1-2%) over the whole energy range predicting a Q_max≈2. In conclusion, the approximation Q=RBE=1 was found to be valid only above ~50 keV whereas at lower energies both Q and RBE become strongly energy-dependent. It is envisioned that the present work will contribute towards establishing robust methodologies to determine theoretically the energy-dependence of radiation quality of individual electrons which may then be used in subsequent calculations involving practical electron and photon radiation sources.

A Stochastic Model of DNA Double-Strand Breaks Repair Throughout the Cell Cycle

Cell cycle phase is a decisive factor in determining the repair pathway of DNA double-strand breaks (DSBs) by non-homologous end joining (NHEJ) or homologous recombination (HR). Recent experimental studies revealed that 53BP1 and BRCA1 are the key mediators of the DNA damage response (DDR) with antagonizing roles in choosing the appropriate DSB repair pathway in G1, S, and G2 phases. Here, we present a stochastic model of biochemical kinetics involved in detecting and repairing DNA DSBs induced by ionizing radiation during the cell cycle progression. A three-dimensional stochastic process is defined to monitor the cell cycle phase and DSBs repair at times after irradiation. To estimate the model parameters, a Metropolis Monte Carlo method is applied to perform maximum likelihood estimation utilizing the kinetics of γ-H2AX and RAD51 foci formation in G1, S, and G2 phases. The recruitment of DSB repair proteins is verified by comparing our model predictions with the corresponding experimental data on human cells after exposure to X and γ-radiation. Furthermore, the interaction between 53BP1 and BRCA1 is simulated for G1 and S/G2 phases determining the competition between NHEJ and HR pathways in repairing induced DSBs throughout the cell cycle. In accordance with recent biological data, the numerical results demonstrate that the maximum proportion of HR occurs in S phase cells and the high level of NHEJ takes place in G1 and G2 phases. Moreover, the stochastic realizations of the total yield of simple and complex DSBs ligation are compared for G1 and S/G2 damaged cells. Finally, the proposed stochastic model is validated when DSBs induced by different particle radiation such as iron, silicon, oxygen, proton, and carbon.

Modelling variable proton relative biological effectiveness for treatment planning

Dose in proton radiotherapy is generally prescribed by scaling the physical proton dose by a constant value of 1.1. Relative biological effectiveness (RBE) is defined as the ratio of doses required by two radiation modalities to cause the same level of biological effect. The adoption of an RBE of 1.1. assumes that the biological efficacy of protons is similar to photons, allowing decades of clinical dose prescriptions from photon treatments and protocols to be utilized in proton therapy. There is, however, emerging experimental evidence that indicates that proton RBE varies based on technical, tissue and patient factors. The notion that a single scaling factor may be used to equate the effects of photons and protons across all biological endpoints and doses is too simplistic and raises concern for treatment planning decisions. Here, we review the models that have been developed to better predict RBE variations in tissue based on experimental data as well as using a mechanistic approach.

Clustered DNA Damages induced by 0.5 to 30 eV Electrons

Low-energy electrons (LEEs) of energies ≤30 eV are generated in large quantities by ionizing radiation. These electrons can damage DNA; particularly, they can induce the more detrimental clustered lesions in cells. This type of lesions, which are responsible for a large portion of the genotoxic stress generated by ionizing radiation, is described in the Introduction. The reactions initiated by the collisions of 0.5-30 eV electrons with oligonucleotides, duplex DNA, and DNA bound to chemotherapeutic platinum drugs are explained and reviewed in the subsequent sections. The experimental methods of LEE irradiation and DNA damage analysis are described with an emphasis on the detection of cluster lesions, which are considerably enhanced in DNA-Pt-drug complexes. Based on the energy dependence of damage yields and cross-sections, a mechanism responsible for the clustered lesions can be attributed to the capture of a single electron by the electron affinity of an excited state of a base, leading to the formation of transient anions at 6 and 10 eV. The initial capture is followed by electronic excitation of the base and dissociative attachment-at other DNA sites-of the electron reemitted from the temporary base anion. The mechanism is expected to be universal in the cellular environment and plays an important role in the formation of clustered lesions.

MODELLING DSB REPAIR KINETICS FOR DNA DAMAGE INDUCED BY PROTON AND CARBON IONS

Proton and carbon therapy are the main choices of particle therapy for cancer treatment. Particle dose distribution is superior to conventional photon therapy dose distribution due to Bragg peak. However, the basic biology of cellular damage and cell death is not well understood. The aim of this work is to present a mechanistic model of double strand break (DSB) repair that predicts the repair kinetics of damage induced by particles employed in cancer therapy. Monte Carlo Track Damage Simulation (MCDS) was employed to model DNA damage. The frequency of DSB and SSB was computed for proton and carbon ions. DSBs were subjected to repair model to calculate the repair kinetics. Two distinct DSB repair models dependent on the cell cycle were proposed. The DSB repair model contains non-homologous end joining (NHEJ), homologous recombination (HR) and back up non-homologous end joining (B-NHEJ) repair processes. The DSB complexity results in the switch in the repair pathway from NHEJ to a slower process that starts with DSB end resection. DSB end resection in early S and G1 phases of the cell cycle enhances the B-NHEJ repair pathway, while in late S and G2 phases of the cell cycle promotes HR repair pathway. The repair model was transformed to a set of nonlinear differential equations. The model calculates the overall repair kinetics and protein temporal repair activity at the site of damage. The damage and repair model provides a detailed mechanistic understanding of all processes that are involved in the damage induction and repair. The number of DSB and their complexity increase as the particle energy decreases due to the proximity of particle interactions in water. The repair kinetics show a biphasic behaviour that is due to the NHEJ fast repair of simple type DSB and HR slow repair of complex type DSB.

Mechanistic Modelling of Radiation Responses

Radiobiological modelling has been a key part of radiation biology and therapy for many decades, and many aspects of clinical practice are guided by tools such as the linear-quadratic model. However, most of the models in regular clinical use are abstract and empirical, and do not provide significant scope for mechanistic interpretation or making predictions in novel cell lines or therapies. In this review, we will discuss the key areas of ongoing mechanistic research in radiation biology, including physical, chemical, and biological steps, and review a range of mechanistic modelling approaches which are being applied in each area, highlighting the possible opportunities and challenges presented by these techniques.

Modeling the interplay between DNA-PK, Artemis, and ATM in non-homologous end-joining repair in G1 phase of the cell cycle

Modeling a biological process equips us with more comprehensive insight into the process and a more advantageous experimental design. Non-homologous end joining (NHEJ) is a major double-strand break (DSB) repair pathway that occurs throughout the cell cycle. The objective of the current work is to model the fast and slow phases of NHEJ in G1 phase of the cell cycle following exposure to ionizing radiation (IR). The fast phase contains the major components of NHEJ; Ku70/80 complex, DNA-dependent protein kinase catalytic subunit (DNA-PKcs), and XLF/XRCC4/ligase IV complex (XXL). The slow phase in G1 phase of the cell cycle is associated with more complex lesions and involves ATM and Artemis proteins in addition to the major components. Parameters are mainly obtained from experimental data. The model is successful in predicting the kinetics of DSB foci in 13 normal, ATM-deficient, and Artemis-deficient mammalian fibroblast cell lines in G1 phase of the cell cycle after exposure to low doses of IR. The involvement of ATM provides the model with the potency to be connected to different signaling pathways. Ku70/80 concentration and DNA-binding rate as well as XXL concentration and enzymatic activity are introduced as the best targets for affecting NHEJ DSB repair process. On the basis of the current model, decreasing concentration and DNA binding rate of DNA-PKcs is more effective than inhibiting its activity towards the Artemis protein.

A New Standard DNA Damage (SDD) Data Format

Our understanding of radiation-induced cellular damage has greatly improved over the past few decades. Despite this progress, there are still many obstacles to fully understand how radiation interacts with biologically relevant cellular components, such as DNA, to cause observable end points such as cell killing. Damage in DNA is identified as a major route of cell killing. One hurdle when modeling biological effects is the difficulty in directly comparing results generated by members of different research groups. Multiple Monte Carlo codes have been developed to simulate damage induction at the DNA scale, while at the same time various groups have developed models that describe DNA repair processes with varying levels of detail. These repair models are intrinsically linked to the damage model employed in their development, making it difficult to disentangle systematic effects in either part of the modeling chain. These modeling chains typically consist of track-structure Monte Carlo simulations of the physical interactions creating direct damages to DNA, followed by simulations of the production and initial reactions of chemical species causing so-called "indirect" damages. After the induction of DNA damage, DNA repair models combine the simulated damage patterns with biological models to determine the biological consequences of the damage. To date, the effect of the environment, such as molecular oxygen (normoxic vs. hypoxic), has been poorly considered. We propose a new standard DNA damage (SDD) data format to unify the interface between the simulation of damage induction in DNA and the biological modeling of DNA repair processes, and introduce the effect of the environment (molecular oxygen or other compounds) as a flexible parameter. Such a standard greatly facilitates inter-model comparisons, providing an ideal environment to tease out model assumptions and identify persistent, underlying mechanisms. Through inter-model comparisons, this unified standard has the potential to greatly advance our understanding of the underlying mechanisms of radiation-induced DNA damage and the resulting observable biological effects when radiation parameters and/or environmental conditions change.

Track to the future: historical perspective on the importance of radiation track structure and DNA as a radiobiological target

PURPOSE

Understanding the mechanisms behind induced biological response following exposure to ionizing radiation is not only important in assessing the risk associated with human exposure, but potentially can help identify ways of improving the efficacy of radiotherapy. Over the decades, there has been much discussion on what is the key biological target for radiation action and its associated size. It was already known in the 1930s that microscopic features of radiation significantly influenced biological outcomes. This resulted in the development of classic target theory, leading to field of microdosimetry and subsequently nanodosimetry, studying the inhomogeneity and stochastics of interactions, along with the identification of DNA as a key target.

CONCLUSIONS

Ultimately, the biological response has been found to be dependent on the radiation track structure (spatial and temporal distribution of ionization and excitation events). Clustering of energy deposition on the nanometer scale has been shown to play a critical role in determining biological response, producing not just simple isolated DNA lesions but also complex clustered lesions that are more difficult to repair. The frequency and complexity of these clustered damage sites are typically found to increase with increasing LET. However in order to fully understand the consequences, it is important to look at the relative distribution of these lesions over larger dimensions along the radiation track, up to the micrometer scale. Correlation of energy deposition events and resulting sites of DNA damage can ultimately result in complex gene mutations and complex chromosome rearrangements following repair, with the frequency and spectrum of the resulting rearrangements critically dependent on the spatial and temporal distribution of these sites and therefore the radiation track. Due to limitations in the techniques used to identify these rearrangements it is likely that the full complexity of the genetic rearrangements that occur has yet to be revealed. This paper discusses these issues from a historical perspective, with many of these historical studies still having relevance today. These can not only cast light on current studies but guide future studies, especially with the increasing range of biological techniques available. So, let us build on past knowledge to effectively explore the future.

Geant4-DNA example applications for track structure simulations in liquid water: A report from the Geant4-DNA Project

This Special Report presents a description of Geant4-DNA user applications dedicated to the simulation of track structures (TS) in liquid water and associated physical quantities (e.g., range, stopping power, mean free path…). These example applications are included in the Geant4 Monte Carlo toolkit and are available in open access. Each application is described and comparisons to recent international recommendations are shown (e.g., ICRU, MIRD), when available. The influence of physics models available in Geant4-DNA for the simulation of electron interactions in liquid water is discussed. Thanks to these applications, the authors show that the most recent sets of physics models available in Geant4-DNA (the so-called "option4" and "option 6" sets) enable more accurate simulation of stopping powers, dose point kernels, and W-values in liquid water, than the default set of models ("option 2") initially provided in Geant4-DNA. They also serve as reference applications for Geant4-DNA users interested in TS simulations.

Small de novo CNVs as biomarkers of parental exposure to low doses of ionizing radiation of caesium-137

The radiological accident in Goiania in 1987 caused a trail of human contamination, animal, plant and environmental by a radionuclide. Exposure to ionizing radiation results in different types of DNA lesions. The mutagenic effects of ionizing radiation on the germline are special concern because they can endures for several generations, leading to an increase in the rate of mutations in children of irradiated parents. Thus, to evaluate the biological mechanisms of ionizing radiation in somatic and germline cells, with consequent determination of the rate mutations, is extremely important for the estimation of genetic risks. Recently it was established that Chromosomal Microarray Analysis is an important tool for detecting wide spectra of gains or losses in the human genome. Here we present the results of the effect of accidental exposure to low doses of ionizing radiation on the formation of CNVs in the progeny of a human population accidentally exposed to Caesium-137 during the radiological accident in Goiânia, Brazil.

Computational analysis of androgen receptor dependent radiosensitivity in prostate cancer

In this study, we quantitatively analyze the mechanism by which androgen deprivation therapy (ADT) is enhancing radiosensitivity in prostate cancer (PCa) patients. It has been shown in laboratory experiments, as well as in patient data in the literature, that the androgen receptor (AR) reduces the effectiveness of ionizing radiation treatment by enhancing the non-homologous end joining (NHEJ) repair of radiation damage. The suppression of AR by ADT suppresses the activity of NHEJ that leads to radiosensitivity in PCa patients. In this paper, we have studied this positive interaction between AR and NHEJ using mathematical models of the NHEJ that we have developed using both the experimental and clinical data for PCa. Our results show that the biological observation of suppression of AR by ADT leading to down-regulation of the first NHEJ protein Ku and NHEJ is a plausible biological mechanism that explains both the experimental and clinical observations in the literature. The presented analysis is the first step in quantitatively analyzing possible treatment scenarios to find the optimal treatment strategies for PCa using the combination treatment with ADT, NHEJ inhibitors, and IR.

In vitro effects of tetraiodothyroacetic acid combined with X-irradiation on basal cell carcinoma cells

We investigated radiosensitization in an untreated basal cell carcinoma (TE.354.T) cell line and post-pretreatment with tetraiodothyroacetic acid (tetrac) X 1 h at 37°C, 0.2 and 2.0 µM tetrac. Radioresistant TE.354.T cells were grown in modified medium containing fibroblast growth factor-2, stem cell factor-1 and a reduced calcium level. We also added reproductively inactivated (30 Gy) “feeder cells” to the medium. The in vitro doubling time was 34.1 h, and the colony forming efficiency was 5.09 percent. These results were therefore suitable for clonogenic radiation survival assessment. The 250 kVp X-ray survival curve of control TE.354.T cells showed linear-quadratic survival parameters of α_X-ray = 0.201 Gy⁻¹ and β_X-ray = 0.125 Gy⁻². Tetrac concentrations of either 0.2 or 2.0 µM produced α_X-ray and β_X-ray parameters of 2.010 and 0.282 Gy⁻¹ and 2.050 and 0.837 Gy⁻², respectively. The surviving fraction at 2 Gy (SF₂) for control cells was 0.581, while values for 0.2 and 2.0 µM tetrac were 0.281 and 0.024. The SF₂ data show that tetrac concentrations of 0.2 and 2.0 µM sensitize otherwise radioresistant TE.354.T cells by factors of 2.1 and 24.0, respectively. Thus, radioresistant basal cell carcinoma cells may be radiosensitized pharmacologically by exposure to tetrac.

Radiation track, DNA damage and response-a review

The purpose of this paper has been to review the current status and progress of the field of radiation biophysics, and draw attention to the fact that physics, in general, and radiation physics in particular, with the aid of mathematical modeling, can help elucidate biological mechanisms and cancer therapies. We hypothesize that concepts of condensed-matter physics along with the new genomic knowledge and technologies and mechanistic mathematical modeling in conjunction with advances in experimental DNA (Deoxyrinonucleic acid molecule) repair and cell signaling have now provided us with unprecedented opportunities in radiation biophysics to address problems in targeted cancer therapy, and genetic risk estimation in humans. Obviously, one is not dealing with 'low-hanging fruit', but it will be a major scientific achievement if it becomes possible to state, in another decade or so, that we can link mechanistically the stages between the initial radiation-induced DNA damage; in particular, at doses of radiation less than 2 Gy and with structural changes in genomic DNA as a precursor to cell inactivation and/or mutations leading to genetic diseases. The paper presents recent development in the physics of radiation track structure contained in the computer code system KURBUC, in particular for low-energy electrons in the condensed phase of water for which we provide a comprehensive discussion of the dielectric response function approach. The state-of-the-art in the simulation of proton and carbon ion tracks in the Bragg peak region is also presented. The paper presents a critical discussion of the models used for elastic scattering, and the validity of the trajectory approach in low-electron transport. Brief discussions of mechanistic and quantitative aspects of microdosimetry, DNA damage and DNA repair are also included as developed by the authors' work.

Mechanistic Modelling and Bayesian Inference Elucidates the Variable Dynamics of Double-Strand Break Repair

DNA double-strand breaks are lesions that form during metabolism, DNA replication and exposure to mutagens. When a double-strand break occurs one of a number of repair mechanisms is recruited, all of which have differing propensities for mutational events. Despite DNA repair being of crucial importance, the relative contribution of these mechanisms and their regulatory interactions remain to be fully elucidated. Understanding these mutational processes will have a profound impact on our knowledge of genomic instability, with implications across health, disease and evolution. Here we present a new method to model the combined activation of non-homologous end joining, single strand annealing and alternative end joining, following exposure to ionising radiation. We use Bayesian statistics to integrate eight biological data sets of double-strand break repair curves under varying genetic knockouts and confirm that our model is predictive by re-simulating and comparing to additional data. Analysis of the model suggests that there are at least three disjoint modes of repair, which we assign as fast, slow and intermediate. Our results show that when multiple data sets are combined, the rate for intermediate repair is variable amongst genetic knockouts. Further analysis suggests that the ratio between slow and intermediate repair depends on the presence or absence of DNA-PKcs and Ku70, which implies that non-homologous end joining and alternative end joining are not independent. Finally, we consider the proportion of double-strand breaks within each mechanism as a time series and predict activity as a function of repair rate. We outline how our insights can be directly tested using imaging and sequencing techniques and conclude that there is evidence of variable dynamics in alternative repair pathways. Our approach is an important step towards providing a unifying theoretical framework for the dynamics of DNA repair processes.

The quest to exploit the Auger effect in cancer radiotherapy - a reflective review

To identify the emergence of the recognition of the potential of the Auger effect for clinical application, and after tracing the salient milestones towards that goal, to evaluate the status quo and future prospects. It was not until 40 years after the discovery of Auger electrons, that the availability of radioactive DNA precursors enabled the biological power, and the clinical potential, of the Auger effect to be appreciated. Important milestones on the path to clinical translation have been identified and reached, but hurdles remain. Nevertheless the potential is still evident, and there is reasonable optimism that the goal of clinical translation is achievable.

DNA double-strand break repair: a theoretical framework and its application

DNA double-strand breaks (DSBs) are formed as a result of genotoxic insults, such as exogenous ionizing radiation, and are among the most serious types of DNA damage. One of the earliest molecular responses following DSB formation is the phosphorylation of the histone H2AX, giving rise to γH2AX. Many copies of γH2AX are generated at DSBs and can be detected in vitro as foci using well-established immuno-histochemical methods. It has previously been shown that anti-γH2AX antibodies, modified by the addition of the cell-penetrating peptide TAT and a fluorescent or radionuclide label, can be used to visualize and quantify DSBs in vivo. Moreover, when labelled with a high amount of the short-range, Auger electron-emitting radioisotope, (111)In, the amount of DNA damage within a cell can be increased, leading to cell death. In this report, we develop a mathematical model that describes how molecular processes at individual sites of DNA damage give rise to quantifiable foci. Equations that describe stochastic mean behaviours at individual DSB sites are derived and parametrized using population-scale, time-series measurements from two different cancer cell lines. The model is used to examine two case studies in which the introduction of an antibody (anti-γH2AX-TAT) that targets a key component in the DSB repair pathway influences system behaviour. We investigate: (i) how the interaction between anti-γH2AX-TAT and γH2AX effects the kinetics of H2AX phosphorylation and DSB repair and (ii) model behaviour when the anti-γH2AX antibody is labelled with Auger electron-emitting (111)In and can thus instigate additional DNA damage. This work supports the conclusion that DSB kinetics are largely unaffected by the introduction of the anti-γH2AX antibody, a result that has been validated experimentally, and hence the hypothesis that the use of anti-γH2AX antibody to quantify DSBs does not violate the image tracer principle. Moreover, it provides a novel model of DNA damage accumulation in the presence of Auger electron-emitting (111)In that is supported qualitatively by the available experimental data.

Track structure modeling in liquid water: A review of the Geant4-DNA very low energy extension of the Geant4 Monte Carlo simulation toolkit

Understanding the fundamental mechanisms involved in the induction of biological damage by ionizing radiation remains a major challenge of today's radiobiology research. The Monte Carlo simulation of physical, physicochemical and chemical processes involved may provide a powerful tool for the simulation of early damage induction. The Geant4-DNA extension of the general purpose Monte Carlo Geant4 simulation toolkit aims to provide the scientific community with an open source access platform for the mechanistic simulation of such early damage. This paper presents the most recent review of the Geant4-DNA extension, as available to Geant4 users since June 2015 (release 10.2 Beta). In particular, the review includes the description of new physical models for the description of electron elastic and inelastic interactions in liquid water, as well as new examples dedicated to the simulation of physicochemical and chemical stages of water radiolysis. Several implementations of geometrical models of biological targets are presented as well, and the list of Geant4-DNA examples is described.

Spectrum of Radiation-Induced Clustered Non-DSB Damage - A Monte Carlo Track Structure Modeling and Calculations

The aim of this report is to present the spectrum of initial radiation-induced cellular DNA damage [with particular focus on non-double-strand break (DSB) damage] generated by computer simulations. The radiation types modeled in this study were monoenergetic electrons (100 eV-1.5 keV), ultrasoft X-ray photons Ck, AlK and TiK, as well as some selected ions including 3.2 MeV/u proton; 0.74 and 2.4 MeV/u helium ions; 29 MeV/u nitrogen ions and 950 MeV/u iron ions. Monte Carlo track structure methods were used to simulate damage induction by these radiation types in a cell-mimetic condition from a single-track action. The simulations took into account the action of direct energy deposition events and the reaction of hydroxyl radicals on atomistic linear B-DNA segments of a few helical turns including the water of hydration. Our results permitted the following conclusions: a. The absolute levels of different types of damage [base damage, simple and complex single-strand breaks (SSBs) and DSBs] vary depending on the radiation type; b. Within each damage class, the relative proportions of simple and complex damage vary with radiation type, the latter being higher with high-LET radiations; c. Overall, for both low- and high-LET radiations, the ratios of the yields of base damage to SSBs are similar, being about 3.0 ± 0.2; d. Base damage contributes more to the complexity of both SSBs and DSBs, than additional SSB damage and this is true for both low- and high-LET radiations; and e. The average SSB/DSB ratio for low-LET radiations is about 18, which is about 5 times higher than that for high-LET radiations. The hypothesis that clustered DNA damage is more difficult for cells to repair has gained currency among radiobiologists. However, as yet, there is no direct in vivo experimental method to validate the dependence of kinetics of DNA repair on DNA damage complexity (both DSB and non-DSB types). The data on the detailed spectrum of DNA damage presented here, in particular the non-DSB type, provide a good basis for testing mechanistic models of DNA repair kinetics such as base excision repair.

DSB repair model for mammalian cells in early S and G1 phases of the cell cycle: application to damage induced by ionizing radiation of different quality

The purpose of this work is to test the hypothesis that kinetics of double strand breaks (DSB) repair is governed by complexity of DSB. To test the hypothesis we used our recent published mechanistic mathematical model of DSB repair for DSB induced by selected protons, deuterons, and helium ions of different energies representing radiations of different qualities. In light of recent advances in experimental and computational techniques, the most appropriate method to study cellular responses in radiation therapy, and exposures to low doses of ionizing radiations is using mechanistic approaches. To this end, we proposed a 'bottom-up' approach to study cellular response that starts with the DNA damage. Monte Carlo track structure method was employed to simulate initial damage induced in the genomic DNA by direct and indirect effects. Among the different types of DNA damage, DSB are known to be induced in simple and complex forms. The DSB repair model in G1 and early S phases of the cell cycle was employed to calculate the repair kinetics. The model considers the repair of simple and complex DSB, and the DSB produced in the heterochromatin. The inverse sampling method was used to calculate the repair kinetics for each individual DSB. The overall repair kinetics for 500 DSB induced by single tracks of the radiation under test were compared with experimental results. The results show that the model is capable of predicting the repair kinetics for the DSB induced by radiations of different qualities within an accepted range of uncertainty.

Genome-based, mechanism-driven computational modeling of risks of ionizing radiation: The next frontier in genetic risk estimation?

Research activity in the field of estimation of genetic risks of ionizing radiation to human populations started in the late 1940s and now appears to be passing through a plateau phase. This paper provides a background to the concepts, findings and methods of risk estimation that guided the field through the period of its growth to the beginning of the 21st century. It draws attention to several key facts: (a) thus far, genetic risk estimates have been made indirectly using mutation data collected in mouse radiation studies; (b) important uncertainties and unsolved problems remain, one notable example being that we still do not know the sensitivity of human female germ cells to radiation-induced mutations; and (c) the concept that dominated the field thus far, namely, that radiation exposures to germ cells can result in single gene diseases in the descendants of those exposed has been replaced by the concept that radiation exposure can cause DNA deletions, often involving more than one gene. Genetic risk estimation now encompasses work devoted to studies on DNA deletions induced in human germ cells, their expected frequencies, and phenotypes and associated clinical consequences in the progeny. We argue that the time is ripe to embark on a human genome-based, mechanism-driven, computational modeling of genetic risks of ionizing radiation, and we present a provisional framework for catalyzing research in the field in the 21st century.

Modeling dose deposition and DNA damage due to low-energy β(-) emitters

One of the main issues of low-energy internal emitters concerns the very short ranges of the beta particles, versus the dimensions of the biological targets. Depending on the chemical form, the radionuclide may be more concentrated either in the cytoplasm or in the nucleus of the target cell. Consequently, since in most cases conventional dosimetry neglects this issue it may overestimate or underestimate the dose to the nucleus and hence the biological effects. To assess the magnitude of these deviations and to provide a realistic evaluation of the localized energy deposition by low-energy internal emitters, the biophysical track-structure code PARTRAC was used to calculate nuclear doses, DNA damage yields and fragmentation patterns for different localizations of radionuclides in human interphase fibroblasts. The nuclides considered in the simulations were tritium and nickel-63, which emit electrons with average energies of 5.7 (range in water of 0.42 μm) and 17 keV (range of 5 μm), respectively, covering both very short and medium ranges of beta-decay products. The simulation results showed that the largest deviations from the conventional dosimetry occur for inhomogeneously distributed short-range emitters. For uniformly distributed radionuclides selectively in the cytoplasm but excluded from the cell nucleus, the dose in the nucleus is 15% of the average dose in the cell in the case of tritium but 64% for nickel-63. Also, the numbers of double-strand breaks (DSBs) and the distributions of DNA fragments depend on subcellular localization of the radionuclides. In the low- and medium-dose regions investigated here, DSB numbers are proportional to the nuclear dose, with about 50 DSB/Gy for both studied nuclides. In addition, DSB numbers on specific chromosomes depend on the radionuclide localization in the cell as well, with chromosomes located more peripherally in the cell nucleus being more damaged by short-ranged emitters in cytoplasm compared with chromosomes located more centrally. These results illustrate the potential for over- or underestimating the risk associated with low-energy emitters, particularly for tritium intake, when their distribution at subcellular levels is not appropriately considered.

Radiation induced base excision repair (BER): a mechanistic mathematical approach

This paper presents a mechanistic model of base excision repair (BER) pathway for the repair of single-stand breaks (SSBs) and oxidized base lesions produced by ionizing radiation (IR). The model is based on law of mass action kinetics to translate the biochemical processes involved, step-by-step, in the BER pathway to translate into mathematical equations. The BER is divided into two subpathways, short-patch repair (SPR) and long-patch repair (LPR). SPR involves in replacement of single nucleotide via Pol β and ligation of the ends via XRCC1 and Ligase III, while LPR involves in replacement of multiple nucleotides via PCNA, Pol δ/ɛ and FEN 1, and ligation via Ligase I. A hallmark of IR is the production of closely spaced lesions within a turn of DNA helix (named complex lesions), which have been attributed to a slower repair process. The model presented considers fast and slow component of BER kinetics by assigning SPR for simple lesions and LPR for complex lesions. In the absence of in vivo reaction rate constants for the BER proteins, we have deduced a set of rate constants based on different published experimental measurements including accumulation kinetics obtained from UVA irradiation, overall SSB repair kinetic experiments, and overall BER kinetics from live-cell imaging experiments. The model was further used to calculate the repair kinetics of complex base lesions via the LPR subpathway and compared to foci kinetic experiments for cells irradiated with γ rays, Si, and Fe ions. The model calculation show good agreement with experimental measurements for both overall repair and repair of complex lesions. Furthermore, using the model we explored different mechanisms responsible for inhibition of repair when higher LET and HZE particles are used and concluded that increasing the damage complexity can inhibit initiation of LPR after the AP site removal step in BER.

The BIANCA model/code of radiation-induced cell death: application to human cells exposed to different radiation types

This paper presents a biophysical model of radiation-induced cell death, implemented as a Monte Carlo code called BIophysical ANalysis of Cell death and chromosome Aberrations (BIANCA), based on the assumption that some chromosome aberrations (dicentrics, rings, and large deletions, called ‘‘lethal aberrations’’) lead to clonogenic inactivation. In turn, chromosome aberrations are assumed to derive from clustered, and thus severe, DNA lesions (called ‘‘cluster lesions,’’ or CL) interacting at the micrometer scale; the CL yield and the threshold distance governing CL interaction are the only model parameters. After a pilot study on V79 hamster cells exposed to protons and carbon ions, in the present work the model was extended and applied to AG1522 human cells exposed to photons, He ions, and heavier ions including carbon and neon. The agreement with experimental survival data taken from the literature supported the assumptions. In particular, the inactivation of AG1522 cells was explained by lethal aberrations not only for X-rays, as already reported by others, but also for the aforementioned radiation types. Furthermore, the results are consistent with the hypothesis that the critical initial lesions leading to cell death are DNA cluster lesions having yields in the order of *2 CL Gy-1 cell-1 at low LET and*20 CL Gy-1 cell-1 at high LET, and that the processing of these lesions is modulated by proximity effects at the micrometer scale related to interphase chromatin organization. The model was then applied to calculate the fraction of inactivated cells, as well as the yields of lethal aberrations and cluster lesions, as a function of LET; the results showed a maximum around 130 keV/lm, and such maximum was much higher for cluster lesions and lethal aberrations than for cell inactivation.

Modeling damage complexity-dependent non-homologous end-joining repair pathway

Non-homologous end joining (NHEJ) is the dominant DNA double strand break (DSB) repair pathway and involves several repair proteins such as Ku, DNA-PKcs, and XRCC4. It has been experimentally shown that the choice of NHEJ proteins is determined by the complexity of DSB. In this paper, we built a mathematical model, based on published data, to study how NHEJ depends on the damage complexity. Under an appropriate set of parameters obtained by minimization technique, we can simulate the kinetics of foci track formation in fluorescently tagged mammalian cells, Ku80-EGFP and DNA-PKcs-YFP for simple and complex DSB repair, respectively, in good agreement with the published experimental data, supporting the notion that simple DSB undergo fast repair in a Ku-dependent, DNA-PKcs-independent manner, while complex DSB repair requires additional DNA-PKcs for end processing, resulting in its slow repair, additionally resulting in slower release rate of Ku and the joining rate of complex DNA ends. Based on the numerous experimental descriptions, we investigated several models to describe the kinetics for complex DSB repair. An important prediction of our model is that the rejoining of complex DSBs is through a process of synapsis formation, similar to a second order reaction between ends, rather than first order break filling/joining. The synapsis formation (SF) model allows for diffusion of ends before the synapsis formation, which is precluded in the first order model by the rapid coupling of ends. Therefore, the SF model also predicts the higher number of chromosomal aberrations observed with high linear energy transfer (LET) radiation due to the higher proportion of complex DSBs compared to low LET radiation, and an increased probability of misrejoin following diffusion before the synapsis is formed, while the first order model does not provide a mechanism for the increased effectiveness in chromosomal aberrations observed.

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.