Next Article in Journal
Zoledronic Acid Modulation of TRPV1 Channel Currents in Osteoblast Cell Line and Native Rat and Mouse Bone Marrow-Derived Osteoblasts: Cell Proliferation and Mineralization Effect
Next Article in Special Issue
Targeted Radionuclide Therapy: New Advances for Improvement of Patient Management and Response
Previous Article in Journal
Non-Smoking-Associated Lung Cancer: A distinct Entity in Terms of Tumor Biology, Patient Characteristics and Impact of Hereditary Cancer Predisposition
Previous Article in Special Issue
Integrating Small Animal Irradiators with Functional Imaging for Advanced Preclinical Radiotherapy Research
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Mechanistic Modelling of Radiation Responses

Centre for Cancer Research & Cell Biology, Queen’s University Belfast, Belfast BT9 7AE, Northern Ireland, UK
*
Author to whom correspondence should be addressed.
Cancers 2019, 11(2), 205; https://doi.org/10.3390/cancers11020205
Submission received: 20 November 2018 / Revised: 4 February 2019 / Accepted: 6 February 2019 / Published: 10 February 2019
(This article belongs to the Special Issue New Developments in Radiotherapy)

Abstract

:
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.

1. Introduction

In contrast to many other areas of biology, mathematical modelling has been an essential part of the fields of radiation biology and radiation therapy since their inception. Close links between the disciplines of physics, biology, and medicine have underpinned this connection, fostering a robust field of physical and biological modelling of radiation responses. On the side of physics, advanced mathematical methods underpin the optimization of modern radiotherapy techniques, enabling precise conformation of dose to target tumors. However, mathematical modelling can also contribute to advancing our understanding of fundamental radiobiology, and has done so since the very earliest studies.
Mathematical modelling has been applied in radiobiology since the 1920s, in studies of proliferation in bacteria and chick embryos [1,2]. In both cases, an exponential dependence of cell viability on radiation was obtained, leading to the development of a “target theory” hypothesis of radiation response. This suggested that cells contain a sensitive region (or “target”) which could be inactivated by a damaging event (or “hit”), and that the yield of these hits is proportional to the dose delivered. If radiation-induced is induced randomly among cells, these the number of hits per cell will follow a Poisson distribution, with a mean number per cell of n = D / D 0 , where D is the dose delivered and D 0 is the dose which causes one “hit” on average. If only un-hit cells survive, then the survival probability is the same as the probability of a cell having 0 hits, that is: S = P h i t s ( 0 , D ) = e D D 0 , where P h i t s ( 0 , D ) is the probability of 0 hits following exposure to a dose D . An exponential dependence results, as the probability of seeing 0 hits in a Poisson distribution with mean n is given by e n .
These models proved highly effective in studies on bacteria and other simple systems, but evidence grew that not all response curves showed this simple exponential dependence [3], with some systems showing evidence of a “shouldered” response curve, with initially lower sensitivity to radiation. Interest in this type of model expanded rapidly after the development of robust techniques to assess the in vitro radiation sensitivity of mammalian cells [4,5], which consistently showed shouldered response curves.
Multi-target models sought to address this by suggesting that rather than a single target, cells contained multiple targets, all of which must be inactivated to cause cell death. If each target still had the same e D D 0 hit probability, then the probability of a cell surviving is S = 1 ( 1 e D D 0 ) n , where n is the number of targets in the cell. This gave a survival curve with an initial slope of 0, which gradually increased to 1 / D 0 at high doses (on a log scale). While the multi-target model was able to successfully fit many observed cell-survival datasets, questions remained about its wider applicability. In particular, a growing body of data indicated that many cell lines did not have zero initial slope, and experiments at high doses often showed increasing curvature, up to the limit of detectability [6,7].
One alternative survival model, independently proposed by a number of authors, was the linear-quadratic model S = e α D β D 2 , which incorporated both a linear initial slope (governed by α ) and a continually increasing quadratic curvature (governed by β ). The key features of this Linear-Quadratic (LQ) model are illustrated in Figure 1, showing the contributions of the linear and quadratic response components. The LQ has become the dominant mathematical model in modelling the survival of cells in preclinical studies, in part motivated by its close links to clinical observations. Studies of fractionation effects in tissues and tumors demonstrated that not only did the LQ model effectively reproduce in vitro survival, but also the effect of clinical fractionation, providing a consistent way to interpret these effects [8,9,10].
One important feature of the LQ is that despite its relative simplicity, it provides an intuitive explanation for a range of biological effects. For example, the tissue-sparing benefits of fractionation can be understood by noting that if cells are allowed to repair between two fractions of radiation, they repair “sub-lethal” damage and repeat the shouldered portion of the curve, seeing less killing than if the dose is delivered in a single fraction. The degree of sparing depends on the degree of curvature, typically quantified in terms of the α / β ratio. This has units of Gy, and corresponds to the dose at which the linear and quadratic components have equal contributions to cell killing. Cells with a low α / β ratio see superior sparing compared to those with a high α / β ratio, as illustrated in Figure 1. As a result, the LQ remains the foundation for determining the equivalence of different treatment fractionation schedules in the clinic [11].
However, the LQ is not without its challenges, and there remain outstanding questions about its interpretation. In addition to its applications in the analysis of clonogenic survival and fractionation, quadratic-type responses had also previously been fit to other endpoints, such as the formation of chromosome aberrations and mutations [12,13]. In each case, the general features driving the relationship are broadly similar; for example, that the linear component relates to “single-hit” damage resulting from individual incident particles, and the quadratic component rates to the strength of interaction of accumulated radiation effects from multiple incident particles. However, the exact mechanisms underlying these trends were the subject of some debate.
As the LQ gained increased prominence as a model of survival in the 60s and 70s, a range of quantitative models were developed to attempt to understand these trends. In 1972, Kellerer and Rossi published proposed that the yield of “elementary lesions” followed a quadratic dependence with, dose being formed either directly by a single event, or as “dual lesions” due to the interaction of two sub-lesion events, in the “theory of dual radiation interaction” [14]. In this model, the nature of the “lesions” were deliberately left unspecified. However, the following year Chadwick and Leenhouts published their “molecular theory of cell survival”, which took a similar approach to predicting the sensitivity of cell survival, but explicitly identified the critical lesions as DNA double strand breaks, and sub-lesions as single strand breaks [15]. Although the two models had somewhat different underlying assumptions, in both cases they could reproduce LQ trends at low doses, in good correspondence with experimental observations. However, this was not the only approach taken.
An alternative family of models considered a linear yield of initial lesions, and modeled the quadratic term of radiation response as being due to lesion interaction. Among the first comprehensive model was the Repair-Misrepair (RMR) model of Tobias et al. [16,17]. In this model lesions are induced linearly with dose, and can be repaired with either linear kinetics, or quadratic kinetics proportional to the density of lesions. The former process can be interpreted as correct end joining, and the latter as binary misrepair. Each of these types of joining can be associated with a probability of lethality, giving rise to response curves which can replicate a number of different response models, including target theory models and the LQ.
A number of other models have been developed using similar concepts of linear and quadratic repair processes, with some minor variants on approach. Another extensively studied approach is the Lethal-Potentially Lethal (LPL) model of Curtis [18]. In this model, rather than a single class of lesion, there are “lethal” lesions, which are always fatal to the cell, and “potentially lethal” lesions, which are repaired with linear and quadratic kinetics, as in the RMR. The linear component is treated as correct rejoining, while the quadratic component is taken as incorrect rejoining and also considered to be lethal.
One final class of model, rather than explicitly postulating direct lesion interaction as the driver of misrepair, suggested that cells have a finite repair capacity. One such Saturable Repair model was proposed by Goodhead [19], which suggested that as the number of lesions increased, the probability of successful repair fell, due to a reduction in the cell’s ability to repair all of the induced damage, for example, due to enzyme or energy depletion. As this means additional dose causes increasing lethality, a shouldered response curve can be formed, potentially closely mimicking the LQ model for suitable parameters.
All of these models have several desirable features—they provide a mechanistic link between the LQ and fundamental damage and repair processes, and can potentially be applied to predict yields of other endpoints, such as mutations, as well as the cell-protecting effects of fractionation or the Relative Biological Effectiveness (RBE) of charged particles (discussed in more detail below).
Most of these models are sufficiently general to be applied in a range of scenarios. However, in the majority of cases, the limitations of experimental clonogenic survival data mean that clinical results cannot adequately distinguish between these models, as they can all predict broadly similar dose response curves with suitable parameter selection (reviewed in [20]). It has thus proven challenging to robustly demonstrate which, if any, of these models truly describe the mechanisms underpinning the LQ.
This presents difficulties for the translation of the LQ into modern radiotherapy. The benefits of fractionation are typically explained in terms of the 4 R’s of radiotherapy—repair, reoxygenation, redistribution, and repopulation [21]—and while some of those (most notably repair) are naturally incorporated in the LQ model, others factors are more difficult to incorporate without the inclusion of further empirical fitting parameters. Moreover, there is a growing interest in a range of other factors which alter radiation responses, including radiations of different qualities [22], radiation-response modifying drugs [23], microenvironmental features of the tumor [24], and inter-patient variations in radiation sensitivity [25].
In addition, as the LQ model is applied as an empirical fit to survival, the parameters obtained for a given cell or tissue do not provide any information about other potentially biologically relevant endpoints of radiation exposure, such as the formation of potentially carcinogenic mutations. These limitations, coupled with the burgeoning knowledge of the mechanisms underlying radiation response, have driven considerable interest in the development of more mechanistic models of radiation response and cellular fate.
However, as illustrated in Figure 2, radiation damage is initiated by physical interactions on the timescale of femtoseconds, triggering a cascade of processes that can result in alterations to the cell’s fate hours, days, or even years following the initial exposure, mediated by both chemical and biological processes. While there are many models of both individual processes and combinations relevant to specific endpoints, because of the range of processes and scales involved, there is not as yet a single integrated model of all of these processes. In this review, we will discuss the modelling efforts on these different scales, which attempt to better understand cellular responses to ionizing radiation, and briefly discuss how this may build towards improved individualized predictions.

2. Physical DNA Damage

The genetic content of the cell was putatively identified as the primary target for radiation action some years before the structure of DNA was determined [13], a prediction which has been extensively validated in a range of subsequent studies [26,27,28]. In particular, it has been shown that DNA double strand breaks (DSBs)—events which sever both strands of the DNA in close proximity—are strongly correlated with lethality. As discussed in more detail in Section 3, this is because although cells have a number of pathways which can effectively repair DSBs, these are more error-prone than those which repair simpler events, such as base damages or single-strand breaks [29].
On the macroscopic scale, radiation dose is frequently quoted as a spatially averaged quantity across relatively large volumes—whole cell populations or tumors—but this does not fully reflect the complexity of radiation damage. DNA has structure on the nanometre scale, with a double-helix having a diameter of approximately 2 nm, and DSBs involving breaks separated by approximately two turns or fewer, corresponding to on the order of 10 to 20 base pairs, or 4–8 nm in length. As a result, physical differences in energy deposition on this scale lead to differences in biological effect. While “sparsely ionizing” radiation, such as high energy X-rays, distributes its energy effectively randomly around the DNA, high LET radiation (such as heavy charged particles) deposit their energy much more densely. This leads to damage with greater complexity, potentially in greater amounts per unit dose, which the cell finds more difficult to repair (schematically illustrated in Figure 3).
This elevation of damage means that densely ionizing particles exhibit a high Relative Biological Effectiveness (RBE). This is defined as the ratio between the dose of a sparsely ionizing reference radiation which gives rise to a particular biological effect, and the dose of a densely ionizing radiation which has the same biological effect. The exact magnitude of RBE and its dependence on parameters such as dose and LET is of particular concern in the delivery of radiotherapy with charged particles, such as protons and ions, making this a topic of considerable research interest in radiobiology.
While many techniques are available to model energy depositions on the nanoscale, the “gold standard” is provided by Monte Carlo track structure calculations [30]. Such simulations stochastically model radiation interactions on an event-by-event basis, tracking each scattering, excitation, or ionization event of both the primary and any secondary particles which are produced (giving “tracks” as illustrated in Figure 3). Such models have been in use for more than three decades. Early research focused on microdosimetric energy distributions, assessing the stochastic distribution of energy on the sub-cellular scale [31,32]. This work demonstrated the high degree of heterogeneity which can be seen in the distribution of energy even within cells which are exposed to the same average dose, highlighting the need for more detailed simulations to accurately predict radiation damage following ionizing radiation exposure.
In the early 1990s, significant developments were made advancing this approach from simple microdosimetric calculations to models incorporating key features of DNA structure. This enabled calculations of energy depositions in various sub-regions of DNA, and predictions of the damage to the DNA structure which would result [33,34]. By fitting the sensitivity of various parts of the DNA to ionizing events, such as the amount of energy which must be deposited in a sugar-phosphate backbone to cause a strand break, predicted yields of measurable endpoints (such as SSB and DSB) could be tuned to match experimentally observed values, and the models used to interrogate underlying mechanisms and endpoints which were not so readily measurable, such as the yield of complex, multiply-damaged sites [26].
Historically, these Monte Carlo codes were typically developed ad-hoc by individual groups for their particular research purposes, and so were frequently restricted in terms of their energy range and simulated particles (reviewed in [35]). A wide range of codes have been applied in simulations of dosimetry on both the patient and cellular scale, including MCNP [36], EGS [37], FLUKA [38], KURBUC [39], MC4 [40], RITRACK [41], TRAX [42], PHITS [43], PARTRAC [44], and Geant4 [45]. Many Monte Carlo codes make use of a “condensed history” physics model, where some charged particle interactions are condensed into single steps to improve performance (including MCNP, EGS, FLUKA, Geant4). While this accurately reproduces results on the patient scale, it is insufficient to fully explain interactions on the cellular scale. Instead, track-structure models which simulate each interaction individually are needed to fully describe the interaction of radiation on the cellular scale [46,47].
As a result, a number of dedicated codes have been developed for these purposes (including KURBUC, MC4, RITRACK, TRAX, PHITS, and PARTRAC), and general-purpose Monte Carlo codes have been extended to calculate track structure at lower energies, making them applicable to cellular scales. Probably the most widely used code in this area is Geant4, a general-purpose Monte Carlo toolkit originally designed for high-energy charged particle transport, which has been extended to encompass ionizing events at energies relevant to cellular nanodosimetry (down to a few eV per particle) through the Geant4-DNA project [48,49]. Geant4 also underpins the TOPAS and TOPAS-nBio toolkits [50,51], which seek to make these Monte Carlo approaches more generally accessible.
These different codes have seen a range of levels of uptake in different fields and applications, but although all of these codes are capable of simulating nanoscale track structures, there remain a number of outstanding research areas where further development is needed to more robustly predict DNA damage from first principles. These include cross section models, descriptions of radiation chemistry, DNA structures and geometry, and understanding of how these different events combine to give rise to DSBs.
Interaction cross-sections are one of the key pieces of input data needed for Monte Carlo models, describing the probability of various scattering and ionization events as a function of energy and target material. However, the majority of nanoscale Monte Carlo calculations simply treat all organic material as water-equivalent. While this is a good approximation at high energies, at low energies these assumptions may break down [52,53]. As a result, there is ongoing research to attempt to measure and model relevant cross-sections in organic molecules, including DNA, so that they can be more robustly incorporated in the Monte Carlo models to provide better predictions of the damage ionizing radiation may have on DNA.
Radiation chemistry is also an area of substantial development. It is well-known that direct interactions between the ionizing particles and DNA are responsible for only approximately 30% of initial damage for photon irradiation, and a somewhat greater proportion for heavier charged particles. Instead, so-called “indirect” interactions play a dominant role. These are events where radiation interacts with another molecule (predominantly water), creating ions or radical species which proceed to react with the DNA, damaging its structure and leading to strand breaks. While early modelling simulated these radical species in small test volumes [33], simulating such reactions in realistic cell volumes can be extremely time-consuming, because of the large number of species which can be generated. Codes are under development, which provide this functionality more generally, including as part of PARTRAC [44], Geant4-DNA [54], TOPAS-nBio [55], and TRAX-CHEM [56]. However, in most cases these codes are still limited to simulations in liquid water, due to the lack of information on reaction rates in more complex chemical compositions.
Because of the above approximations, actual DNA structures are also frequently neglected. Instead, homogeneous water volumes are simulated, with DNA structures superimposed on the simulation results to determine the distribution of damaging events. While this avoids complications from the lack of cross-sections for relevant targets, the appropriate DNA structure to use remains an outstanding research question. Early work [34] made use of small sections of DNA (equivalent to single DNA strands of tens of base pairs in length), which is sufficient to calculate rates of local damage in an “isolated” DNA section. However, DNA has a number of higher-degree levels of organization, with individual strands wrapped around histones to form nucleosomes, which are in turn packed into larger fibers. This organization means that DNA damage is not randomly distributed, but may depend on local structure around a radiation track. A range of studies have modelled different DNA structures [44,57,58,59], and shown that the choice of a particular DNA structure can have substantial impacts on the predicted yields of different types of damage, and how it is distributed through the genome.
Despite the above uncertainties, numerous models have been generated which predict the yield of radiation-induced DSBs in a variety of systems, based on models of either direct damage alone, or combinations of direct and indirect damage. These range from multi-scale models which simulate damage in detail on small DNA structures and extrapolate to averaged or cell-level responses [58,60,61,62]. There are also a range of whole-cell models incorporating varying levels of detail on DNA damage and chromosome structure, summarized in Table 1. Although these models have differences in their design, underlying physics, whether or not they include chemical effects, and models of biological structure and DNA damage, many of them produce broadly similar predictions.
In part, this is the result of the above-mentioned significant uncertainty in the nature of the physico-chemical processes, which give rise to double-strand breaks at the interface between physical and biological models. Due to the lack of distinction between different modelling assumptions, there is very large uncertainty in the underlying biological parameters, and even quite different models can be brought into a reasonable agreement through careful selection of different parameters. These include the amount of energy which needs to be deposited in a single strand to cause a break (with threshold energies estimated from a few eV up to 40 eV, or in some cases an energy-dependent probability [34,61,67]), the likelihood of radical interaction causing damage, and how different damages in close proximity interact. This problem is compounded by the often significant uncertainties present in experimental measurements of DNA damage, which make determining absolute yields of many types of damage challenging.
Thus, while a range of models exist to predict the yield and distribution of DNA damage for different types of irradiation, there is a significant need to better characterize many of the fundamental parameters governing these models, to refine their underlying structure, and enable more robust comparisons between models. A recent collaboration has proposed a Standard for DNA Damage reporting, enabling the outputs of various models to be efficiently compared, to help facilitate such developments [68].

3. DNA Repair

Although DNA damage is the primary driver of radiation-induced effects, it represents only the initiating event in a cascade of processes. Living systems have evolved a wide range of systems enabling the efficient repair of DNA damage, so only a very small fraction of even double-strand breaks lead to lethality. While radiation induces 30–40 DSBs per cell per Gy, this typically corresponds to less than one lethal event per cell [69]. Failures in these repair processes lead to lethality, either through the failure to repair a subset of DSBs [70,71] or misrepair events, which cause significant chromosomal abnormalities [72]. This is particularly apparent in cell lines with defects in key DNA repair genes, which can be more than an order of magnitude more sensitive than repair-competent cells.
Due to their lethality, the repair of DSBs has been the subject of the majority of the modelling efforts in this area. A number of high-level models of DNA damage were developed to attempt to understand how damage may combine into different types of damage. A number of similar approaches considered proximity effects as the major driver of misrepair, with spatially dense clusters of damage having a high probability of incorrect end joining, leading to potentially deleterious misrepair [73,74,75,76,77,78]. Such models often proved effective at predicting yields of events, such as mutations [76] and chromosome aberrations, which can in turn be linked to cell death in some cases [76,79]. However, most such models do not incorporate detailed models of the underlying repair processes, which limits their capacity to incorporate our growing mechanistic understanding of DNA repair.
As schematically illustrated in Figure 4, there are three primary mechanisms by which DSBs can be repaired: Nonhomologous End Joining (NHEJ), Homologous Recombination (HR), or Microhomology-Mediated End Joining (MMEJ). NHEJ is active throughout the cell cycle, and rapidly repairs most DNA DSBs, although is somewhat error-prone and often leads to small modifications around the site of repair [80]. By contrast, HR makes use of DNA from a matching chromosome to precisely repair DNA without introducing new sequence errors. However, because of the dependence on a template strand, it is limited to activity in late S and G2 phases, where DNA has been replicated [81]. Finally, as a backup to these processes, MMEJ can act to rejoin DNA. In this process, DNA strands are matched based on small regions of homology, typically involving significant resection at break ends, and is very error-prone and almost always leads to at least moderate sequence alteration [82]. A number of potential pathways have been implicated in MMEJ, including a resection-dependent subset of NHEJ [83], as well as a wholly distinct pathway known as alternative-NHEJ, which is independent of most major NHEJ proteins [84].
The kinetics of these different repair pathways has been the subject of a range of experimental and modelling studies. The most common empirical approach describes DNA repair as a single- or multi-exponential process, with different sub-populations of breaks being repaired with different kinetics. However, other explanations for the slowing rate of repair at later times have also been suggested, such as a second-order exponential dependence [85] or more complex repair interplays with multiple kinetic parameters or variable repair-times (summarized in [86]). Some of these models, such as the multi-exponential approach, offer a seemingly natural link between groups of breaks repaired with different kinetics and the various available pathways, although there is evidence that other features of breaks, such as their complexity and the local structure of DNA, may also impact on repair kinetics.
Some models of DNA repair fidelity and survival have sought to incorporate detail on these pathways, reflecting them as different repair processes which still act in a largely probabilistic fashion [76]. While still relatively abstract, such models still enable predictions of differential responses in cells with different genetic backgrounds, building towards individualization of radiation response predictions.
Much more detailed studies seek to fully model the biochemistry of these pathways, tracking the migration and interaction of DNA ends and proteins, either using population biokinetics and probabilistic models [87,88,89] or by modelling individual agents in a Monte Carlo fashion [90,91]. These models have primarily focused on NHEJ, and have demonstrated the ability to reproduce observed trends in DNA repair kinetics and the impact of knocking out some key genes [92], potentially providing a framework to better understand these processes on a quantitative level.
By contrast with NHEJ, no quantitative mechanistic models of the other major DNA repair pathways, HR and MMEJ, have been published in combination with radiation. For HR, this is in part due to its complexity, as it involves template DNA strands and a range of potential mechanisms for the resolution of the resulting D-loop (or Holliday Junction) formed between the repaired and template chromatids, the relative activity of which has yet to be fully determined [93]. In the case of alternative-NHEJ, the pathway has only been recently identified as distinct from NHEJ, and while many of the key genes involved have been identified [94], there has been relatively limited modelling of their function.
To fully model the consequences of ionizing radiation, robust descriptions of these pathways are needed to accurately predict the consequences of repair. This includes not only the activity and fidelity of individual pathways, but also how they interact with one another to determine which pathway repairs a particular break, and how this depends on the physical complexity of damage and the potential presence of other damage within the nucleus.
At present, there are a number of models, based on both simple proximity-based or more sophisticated molecular kinetics, which make predictions that can potentially be experimentally tested, relating to yields of different types of misrepair under different irradiation conditions. However, as with initial physical damage, the extensive parameter space underpinning these responses make robustly fitting and falsifying different models challenging—even relatively simple kinetic models of NHEJ can involve between 10 and 20 fitting parameters [88,92].

4. Cell Fate

Once initial DNA damage distributions and their repair have been determined, there is still the question about the cell’s eventual fate. A wide range of approaches has been taken to incorporate our knowledge of these underlying processes into radiation response models.
As noted above, many early radiation response models involved simplified descriptions of how radiation response processes were believed to proceed. Approaches like the dual action model of Rossi and Kellerer [14], the Repair-Misrepair model of Tobias [16], the lethal-potentially lethal model of Curtis [18], or the molecular theory of Chadwick and Leenhouts [15] represented some of the first efforts of incorporating knowledge of radiation biology into mechanistic response models. Such approaches proved successful at capturing the basic features of the radiation dose response curve, and there are a number of modern developments which have sought to place these on a more biologically robust footing by incorporating additional detail about the underlying biology [89,95,96,97]. However, a range of other approaches also seek to build on these more general assumptions to incorporate further mechanistic details about the underlying radiation physics and biology.
For example, predicting the RBE of high LET radiations has been a topic of considerable interest, due to both the interest in the underlying biology as well as the implications for clinical practice in facilities which seek to use heavy charged particles for the treatment of cancer. Many empirical models have been developed which predict protons RBE-LET dependence by applying simple modifications to LQ parameters, typically of the form α p = p 1 α x + p 2 L E T α x / β x , where α x is the α LQ parameter for protons or X-rays respectively, p 1 and p 2 are empirical fitting parameters, and α x / β x is the ratio for X-ray responses [98,99,100]. Similar modifications are sometimes also proposed for the β term. While broadly successful for proton therapy, such approaches are typically ineffective at the higher LETs used in therapy with heavier ions, such as carbon.
A number of mechanistic explanations of cell inactivation were explored to attempt to better understand the mechanisms underpinning these effects. This includes early work, such as geometric radial energy distribution models developed by Butts and Katz [101] and applications of the dual action model [14], as well as more recent modelling approaches, including the Repair-Misrepair-Fixation model (RMF) [89] and many extensions of the DNA damage models, outlined in Section 3.
However, clinical RBE modelling is dominated by two approaches, which build on microdosimetric simulations and simple assumptions about cellular fate. These are the Local Effect Model (LEM) and the Microdosimetric Kinetic Model (MKM). In both cases, the models seek to take advantage of a combination of biological response parameters determined for X-ray exposure and microdosimetric calculations to predict the effects of high LET irradiation.
The LEM defines radiation response in terms of “lethal events”, such that survival is given by S = e N , where N is the number of lethal events. For X-rays, it can be seen that these lesions are formed at a rate N = α D β D 2 , and the LEM assumes that the same dependence holds in sub-volumes of the nucleus of cells exposed to heterogeneous radiation exposures. Thus, the total yield of lethal events across a cell is given by N = α D r β D r 2 V   d V , where D r is the dose delivered to a point r within the nucleus, and V is the total nuclear volume. It can be seen that because of the D r 2 term, the total yield of lethal lesions will be greater in a heterogeneous exposure than a uniform exposure, which gives rise to elevated RBEs in such systems [102]. More recent versions of the LEM have moved to explicitly consider break repair and misrepair probabilities similar to some other models described above, but these approaches have yet to make the transition to clinical practice [103].
The MKM takes a conceptually similar approach. In it, the nucleus is sub-divided into a large number of small “domains”, each of which can contain either directly lethal events, or potentially lethal events. Potentially lethal events within the same sub-domain then have the chance of mis-rejoining to cause lethal events or successfully repairing, giving rise to a linear-quadratic response for lethal lesions within each micro-dosimetric domain. As in the LEM, rates for the formation of lethal and sub-lethal damages can be determined from the approximately uniform damage caused by X-ray irradiations, although the MKM uses the dose-weighted specific energy per track passing through a domain to calculate the induced damage. Because high LET exposure leads to heterogeneous cellular exposure and elevated risks of multiple sub-lethal events in a single domain, an increased probability of misrepair and thus cell death results [104].
In both cases, by combining these relatively simple biological models with microdosimetric models, good predictions can be obtained for RBEs in a range of systems. As a result, these models have successfully made the transition into clinical use, with the LEM and MKM supporting treatment planning optimization at carbon therapy institutes in Europe and Japan, respectively [105,106]. This success shows the potential of more mechanistic models to inform clinical practice and improve patient outcomes. Moreover, their mechanistic nature has enabled them to be applied and validated outside their initial area of application, with the MKM model being applied in areas such as dose rate effects in brachytherapy [107] and other modalities, such as boron neutron capture therapy [108], while the LEM has been applied to gold nanoparticle enhanced X-ray therapies [109].
However, a major challenge which is presented in this area is what events lead to cell death. For example, while the PARTRAC code is capable of simulating detailed models of initial damage and how the cell resolves it, significant uncertainties in the absolute yield of damage remain, and survival predictions still require assumptions to be made about how cells tolerate this damage [64]. Similar limitations are present in other whole-cell DNA damage codes [66].
One common approach is to assume that some subset of types of repair failure can lead to cell death. These can include combination of DSBs or other lesions into “complex” damages, which are deemed irreparable, those which are determined to be unrepaired at a particular timepoint, or those leading to significant misrepair, such as chromosome aberrations. A number of models have taken this type of approach. For example, the BIANCA [79] model simulates the interaction of “cluster lesions” in DNA to predict yields of chromosome aberrations (and later cell death) as a function of LET. The GLOBLE model [110] applies a similar concept in low LET scenarios, modelling isolated and clustered DSBs to better understand the LQ model and its temporal dependence. The “Multiscale” approach [111] also models high- and low-LET effects by simulating the formation of complex clustered lesions as the event which drives radiation-induced cell death. Finally, a combination of MCDS and RMF models [112] combines predictions of DSB yield with frequency-mean specific energy values to predict the RBE of different ions.
While this represents only a selection of the models making use of this approach, it can be seen that this concept of grouping damages spatially into lesions which either misrepair or are deemed irreparable is very common across a wide range of models, despite substantial differences in modelled endpoint and scenarios of interest. This is likely in part motivated by the good evidence that chromosome aberrations that cause significant genetic rearrangement or loss are associated with a loss of clonal viability [72], and the good correlations which have been observed between unrepaired DSBs and cell death in some lines [71,113,114]. However, while this conceptual approach effectively captures overall trends in radiobiological responses, there remains significant disagreement between models on the exact form of this dependence, and these models have been validated to different degrees against heterogeneous datasets.
Moreover, cell death is a complex process, involving both “passive” and “active” death pathways, which can depend to varying degrees on the underlying cell biology, making it unclear by which process a given cell will actually die [115]. These pathways include apoptosis [116,117], necrosis [118], mitotic catastrophe [119], autophagy [120], and senescence [121]. Because many of these pathways require the active involvement of different genes, they can vary substantially based on both the tissue type of origin of a given cell, and any genetic alterations in the many genes involved in these processes. This can significantly impact on the overall sensitivity of different cell lines, for example, with many lymphoid cells showing dramatically higher levels of apoptosis and radiation sensitivity than cells derived from epithelial tissue or fibroblasts [122,123].
These dependencies have been extensively studied in the literature, and some published response models incorporate a high-level approximation of these pathways to reflect differences between such cell lines. However, there is again relatively limited modelling of the mechanisms of these pathways in determining radiation response. While some models do exist, they frequently focus on the pathways more generally, including seeking to understand intrinsic and signalling-driven apoptosis via chemotherapy, rather than that driven by radiation-induced damage [124,125,126,127,128].
Despite the significant complexity underlying these different pathways, improving our understanding of them represents an area of significant unexplored potential in our attempts to better understand radiation responses.

5. Tissue-Level Responses

Although this review has primarily focused on the response of individual cells, it is important to note that clinical outcomes are defined at the multi-cellular level, in terms of responses in tumors and normal tissues. While the intrinsic sensitivity of individual cells is believed to play a significant role in radiosensitivity in both preclinical and clinical models [129,130,131], other effects also potentially become significant.
While cell killing is associated with toxicity and tumor control, the precise amount of killing and the distribution throughout the target which gives rise to different endpoints is highly variable. In particular, while many tumor control probability (TCP) models assume that control depends on killing all (clonogenically viable) tumor cells, normal tissue complication probability (NTCP) can depend on the distribution of dose throughout the organ in a more complex fashion. Some organs respond to irradiation with a severity which approximately follows the average dose across the whole organ (“parallel” organs such as the lung or liver), while other organs can see significant toxicity following irradiation of only a small sub-region (“serial” organs such as the spinal cord or bowel).
A variety of modelling approaches have been proposed to analyze this problem, beginning from Poisson-based probabilistic models [132], which have been generalized using abstract approaches, such as the generalized Equivalent Uniform Dose (EUD) [133,134], which seek to incorporate these volume dependencies for different tissue types. For NTCP, two models are commonly used. The first is the Lyman-Kutcher-Burman (LKB) model, which applies a tissue-specific empirical volume scaling to correct for irradiation of different sub-volumes to different doses, reflecting the difference between serial and parallel organs [135,136]. An alternative process, which seeks to incorporate knowledge of the structure of different tissue types, is the Relative Seriality model [137]. This model described organs in terms of their “seriality”, with highly serial organs sensitive to toxicity arising due to damage in very small volumes. However, as both of these approaches tend to abstract out significant amounts of mechanistic information, there is also interest in more mechanistic approaches, such as tumor models, which incorporate more details of proliferation, heterogeneity, and hypoxia [138,139], and normal tissue models, which more explicitly define “functional sub-units” of the organs and how they interrelate to incorporate these structural effects [140].
A final tissue-level endpoint of interest is carcinogenesis. While the majority of predictions of cancer risk following radiation exposure are statistical, based on extrapolation from populations such as the atomic bomb cohorts, the significant uncertainty in cancer risk at low doses means there is considerable interest in using more mechanistic models to mitigate this uncertainty. Some of these models remain stochastic, while incorporating more information on the stages of carcinogenesis, while others seek to develop more comprehensive models from the cell level, simulating DNA repair failure and mutation, together with clonal expansion driving cancer formation [141,142,143,144].
However, in both tumor and normal tissue models, there remains significant uncertainty in model parameters, and there is typically limited incorporation of information from the single-cell level beyond overall radiosensitivity, suggesting a further challenge in translating from cell-level models to clinical responses.
In addition, irradiation in tissues opens the possibility of effects resulting from intercellular communication. The two major pathways of this effect are the radiation induced bystander effect (RIBE), and immune effects. In the RIBE, signals from neighboring cells can cause unirradiated cells to suffer genomic stress and potentially die [145]. This is a departure from traditional radiobiology, where radiation responses are predicted based on the dose delivered to the irradiated cell alone. A number of RIBE models have been published, but are typically not incorporated in clinical planning [146,147,148]. The involvement of the immune system in radiotherapy is an area of rapidly growing interest, with substantial evidence that changes in regulation of immune signals can play a significant role in determining tumor control following radiotherapy [149,150]. As with the RIBE, a number of models of radio-immunotherapy combinations have been published [151], but as this is a new and rapidly developing field, once again they have seen little translation to clinical practice.
The step between cell-based models and tissue-level responses remains a major challenge in radiobiological modelling. Tissue structure introduces the possibility for complex volume effects which are not apparent on the single-cell level, and this significant increase in complexity is combined with limited experimental data, presenting a substantial challenge for modelers to develop useful and robustly validated tools to link our extensive preclinical knowledge to clinically-meaningful predictions.

6. Potential Impacts of Modelling Advances

A recurring theme in many of the above discussions is the substantial biological complexity that underlies many aspects of radiation response, which is obscured by the relatively simple response models that can successfully characterize both pre-clinical and clinical responses. In many cases, it is reasonable to question how much benefit will be derived from the development of these much more sophisticated models, or if relatively simple, well-established approaches may be sufficient to describe all aspects of radiation response.
However, there are a number of major outstanding clinical questions which these simpler models are not suited to address. One of the most significant of these is that there is now extensive clinical evidence that even in a group of patients with tumors of the same type, there may be a broad range of radiation sensitivities among their individual tumors [132], driven in large part by the cancers’ spectrum of underlying mutations [25]. Estimates of this variation from clinical response curves suggest it could represent large differences in radiation sensitivity even for cancers which would be treated identically in the clinic (e.g., variations in 25% or more in the α component [132]). Consequently, a large fraction of patients are almost certainly being under- and over-dosed by these population-based models, which may significantly impact on clinical outcomes. Empirical population-level models offer no way to address these issues and predict the sensitivity of an individual cancer, as their model parameters are typically at least somewhat abstract and cannot be linked directly to particular genetic pathways or tumor characteristics. This presents a significant challenge in radiation therapy, preventing it from becoming a fully personalized therapy.
Similar issues arise in a number of other areas which can impact on tumor response—the effects of factors such as the tumor microenvironment, hypoxia, and radiation quality are not easily incorporated in these models, at best being addressed by empirical fitting parameters such as OER and RBE effects. These represent significant additional avenues by which treatments could be personalized, if they could be more robustly incorporated in biological models of tissue and tumor response.
As a result, there has been significant interest in the development of clinical signatures of radiation sensitivity. It has been reported that the in vitro clonogenic radiation sensitivity of cells cultured from tumors correlates well with tumor sensitivity in both animal models [130,152] and patients [153,154,155], along with some other intermediate endpoints, such as residual DNA damage in ex vivo irradiated samples. However, such assays are often labor-intensive and difficult to incorporate into a clinical workflow [156]. More recently, there has been a dramatic increase in interest in the use of genetic signatures of radiation response. These techniques seek to associate alterations in the mutation status or expression of genes with radiation sensitivity, typically using purely statistical association. A number of signatures have been reported, either focusing on radiation response in general [157,158], or specific modulators of response, such as signatures of hypoxia [159]. While a promising approach, many of these signatures have faced challenges around reproducibility and robustness [160], as well as questions about if they are truly predictive of radiation response or merely more generally prognostic [161]. In addition, the volume of data required to fully parameterize such empirical models is often prohibitive—as an example, there is still little agreement on the true magnitude of the RBE of protons, despite several hundred studies reporting results in this area [162].
The incorporation of more detailed mechanistic models has the potential to help address these challenges. Models that enable the integration of statistical, population-level data with our growing mechanistic knowledge of the underlying drivers of radiation response can enable the generation of more robust and translatable predictions. Some models have demonstrated the basic potential for predicting aspects of radiation response based on phenotypic or genotypic characteristics, although this remains a significant area of ongoing research.
Such models can potentially be useful in other areas also. For example, radiation is rarely delivered as a single agent in modern clinical practice, instead being combined with chemotherapy or new targeted therapies to improve outcome. While these are typically prescribed based on their activity as single agents, there is growing interest in specifically seeking to develop radiation sensitizing agents [23]. In this approach, agents are usually identified which seek to target specific parts of the radiation response pathway, such as genes involved in DNA repair [163,164,165]. While this approach has identified some radiation sensitizers, translation has been limited. More robust mechanistic models may enable better predictions of the impact of modifying these aspects of DNA repair without the need for extensive preclinical testing.
However, in all cases, it is important to ensure that model development is carried out in a rational way. As noted in the discussion of alternative clonogenic survival models in Section 1 and DNA damage models in Section 2, as model complexity increases, it can become increasingly challenging to distinguish between different models due to a proliferation of adjustable parameters, and added complexity may actually impact negatively on model predictive power [166]. Thus, it is important to ensure that that any additional modelling complexity adds value—either by enabling explanation of novel phenomena, or by enabling additional data sources to be exploited to refine and constrain model predictions—and ideally, be feasible to test and validate in a clinical setting.

7. Conclusions

Mathematical modelling has been a key part of radiation biology and therapy for almost a century, and continues to play an important role in research developments. While a number of relatively simple and largely empirical models dominate current clinical and pre-clinical practice, significant advances are being made in our understanding of all aspects of radiation response, including the initial physics, resulting chemistry, and subsequent biological consequences. Although this is a very complex and challenging area of research, applied correctly it has significant potential to advance our understanding of the drivers of radiation response, and potentially translate this into useful clinical predictive tools in the future, if strong collaboration can be maintained between the diverse disciplines in radiation research.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Curie, M. Sur l’étude des courbes de probabilité relatives à l’action des rayons X sur les bacilles. Comptes Rendus l’Académie des Sci. 1929, 188, 202–204. (In French) [Google Scholar]
  2. Crowther, J.A. Some Considerations Relative to the Action of X-Rays on Tissue Cells. Proc. R. Soc. B Biol. Sci. 1924, 96, 207–211. [Google Scholar] [CrossRef]
  3. Alper, T.; Gillies, N.E.; Elkind, M.M. The sigmoid survival curve in radiobiology. Nature 1960, 186, 1062–1063. [Google Scholar] [CrossRef] [PubMed]
  4. Puck, T.T.; Marcus, P.L. Action of X-rays on Mammalian Cells. J. Exp. Med. 1956, 103, 653–666. [Google Scholar] [CrossRef] [PubMed]
  5. Puck, T.T.; Morkovin, D.; Marcus, P.I.; Cieciura, S.J. Action of x-rays on mammalian cells. II. Survival curves of cells from normal human tissues. J. Exp. Med. 1957, 106, 485–500. [Google Scholar] [CrossRef] [PubMed]
  6. Barendsen, G.; Beusker, T.; Vergroesen, A.; Budke, L. Effects of different ionizing radiations on human cells in tissue culture: II. biological experiments. Radiat. Res. 1960, 13, 841–849. [Google Scholar] [CrossRef]
  7. Sinclair, W.K. The shape of radiation survival curves of mammalian cells cultured in vitro. In Biophysical Aspects of Radiation Quality; International Atomic Energy Agency: Vienna, Austria, 1966. [Google Scholar]
  8. Douglas, B.G.; Fowler, J.F. The effect of multiple small doses of X rays on skin reactions in the mouse and a basic interpretation. Radiat. Res. 1976, 66, 401–426. [Google Scholar] [CrossRef]
  9. Barendsen, G.W. Dose fractionation, dose rate and iso-effect relationships for normal tissue responses. Int. J. Radiat. Oncol. Biol. Phys. 1982, 8, 1981–1997. [Google Scholar] [CrossRef]
  10. Thames, H.D.; Rodney Withers, H.; Peters, L.J.; Fletcher, G.H. Changes in early and late radiation responses with altered dose fractionation: Implications for dose-survival relationships. Int. J. Radiat. Oncol. Biol. Phys. 1982, 8, 219–226. [Google Scholar] [CrossRef]
  11. Fowler, J.F. 21 Years of biologically effective dose. Br. J. Radiol. 2010, 83, 554–568. [Google Scholar] [CrossRef]
  12. Catcheside, D.G.; Lea, D.E.; Thoday, J.M. The production of chromosome structural changes in Tradescantia microspores in relation to dosage, intensity and temperature. J. Genet. 1946, 47, 137–149. [Google Scholar] [CrossRef] [PubMed]
  13. Lea, D.E. Actions of Radiations on Living Cells, 1st ed.; Cambridge University Press: London, UK, 1946. [Google Scholar]
  14. Kellerer, A.M.; Rossi, H.H. The theory of dual radiation action. Curr. Top. Radiat. Res. 1972, 8, 85–158. [Google Scholar]
  15. Chadwick, K.H.; Leenhouts, H.P. A molecular theory of cell survival. Phys. Med. Biol. 1973, 18, 78–87. [Google Scholar] [CrossRef] [PubMed]
  16. Tobias, C.A.; Blakelly, E.A.; Ngo, F.Q.H.; Yang, T.C.H. The Repair-Misrepair Model of Cell Survival. In Radiation Biology and Cancer Research; Meyn, A., Withers, R., Eds.; Raven Press: New York, NY, USA, 1980; pp. 195–230. ISBN 978-3-319-46447-3. [Google Scholar]
  17. Tobias, C.A. The repair-misrepair model in radiobiology: Comparison to other models. Radiat. Res. Suppl. 1985, 8, S77–S95. [Google Scholar] [CrossRef] [PubMed]
  18. Curtis, S.B. Lethal and potentially lethal lesions induced by radiation--a unified repair model. Radiat. Res. 1986, 106, 252–270. [Google Scholar] [CrossRef]
  19. Goodhead, D.T. Saturable repair models of radiation action in mammalian cells. Radiat. Res. Suppl. 1985, 8, S58–S67. [Google Scholar] [CrossRef]
  20. Bodgi, L.; Canet, A.; Pujo-Menjouet, L.; Lesne, A.; Victor, J.M.; Foray, N. Mathematical models of radiation action on living cells: From the target theory to the modern approaches. A historical and critical review. J. Theor. Biol. 2016, 394, 93–101. [Google Scholar] [CrossRef]
  21. Withers, R.H. The Four R’s of Radiotherapy; Academic Press, Inc.: Cambridge, MA, USA, 1975; Volume 5, ISBN 9780120354054. [Google Scholar]
  22. Loeffler, J.S.; Durante, M. Charged particle therapy—optimization, challenges and future directions. Nat. Rev. Clin. Oncol. 2013, 10, 411–424. [Google Scholar] [CrossRef]
  23. Sharma, R.A.; Plummer, R.; Stock, J.K.; Greenhalgh, T.A.; Ataman, O.; Kelly, S.; Clay, R.; Adams, R.A.; Baird, R.D.; Billingham, L.; et al. Clinical development of new drug-radiotherapy combinations. Nat. Rev. Clin. Oncol. 2016, 13, 627–642. [Google Scholar] [CrossRef]
  24. Barker, H.E.; Paget, J.T.E.; Khan, A.A.; Harrington, K.J. The tumour microenvironment after radiotherapy: Mechanisms of resistance and recurrence. Nat. Rev. Cancer 2015, 15, 409–425. [Google Scholar] [CrossRef]
  25. Yard, B.D.; Adams, D.J.; Chie, E.E.; Tamayo, P.; Battaglia, J.S.; Gopal, P.; Rogacki, K.; Pearson, B.E.; Phillips, J.; Raymond, D.P.; et al. A genetic basis for the variation in the vulnerability of cancer to DNA Damage. Nat. Commun. 2016, 7, 11428. [Google Scholar] [CrossRef]
  26. Goodhead, D.T. Initial Events in the Cellular Effects of Ionizing-Radiations—Clustered Damage in DNA. Int. J. Radiat. Biol. 1994, 65, 7–17. [Google Scholar] [CrossRef] [PubMed]
  27. Warters, R.; Hofer, K.; Harris, C.; Smith, J. Radionuclide toxicity in cultured mammalian cells: Elucidation of the primary site of radiation damage. Curr. Top. Radiat. Res. Q. 1978, 12, 389–407. [Google Scholar] [PubMed]
  28. Munro, T.R. The Relative Radiosensitivity of the Nucleus and Cytoplasm of Chinese Hamster Fibroblasts. Radiat. Res. 1970, 42, 451–470. [Google Scholar] [CrossRef] [PubMed]
  29. Ward, J.F.; Blakely, W.F.; Joner, E.I. Mammalian cells are not killed by DNA single-strand breaks caused by hydroxyl radicals from hydrogen peroxide. Radiat. Res. 1985, 103, 383–392. [Google Scholar] [CrossRef] [PubMed]
  30. Nikjoo, H.; Emfietzoglou, D.; Watanabe, R.; Uehara, S. Can Monte Carlo track structure codes reveal reaction mechanism in DNA damage and improve radiation therapy? Radiat. Phys. Chem. 2008, 77, 1270–1279. [Google Scholar] [CrossRef]
  31. Zaider, M.; Brenner, D.J.; Wilson, W.E. The applications of track calculations to radiobiology I. Monte Carlo simulation of proton tracks. Radiat. Res. 1983, 95, 231–247. [Google Scholar] [CrossRef]
  32. Nikjoo, H.; Goodhead, D.T.; Charlton, D.E.; Paretzke, H.G. Energy deposition in small cylindrical targets by ultrasoft X-rays. Phys. Med. Biol. 1989, 34, 691–705. [Google Scholar] [CrossRef]
  33. Nikjoo, P.; O’Neill, D.T.; Goodhead, H. Computational modelling of low-energy electron-induced DNA damage by early physical and chemical events. Int. J. Radiat. Biol. 1997, 71, 467–483. [Google Scholar] [CrossRef]
  34. Nikjoo, H.; O’Neill, P.; Wilson, W.E.; Goodhead, D.T. Computational Approach for Determining the Spectrum of DNA Damage Induced by Ionizing Radiation. Radiat. Res. 2001, 156, 577–583. [Google Scholar] [CrossRef]
  35. Nikjoo, H.; Uehara, S.; Emfietzoglou, D.; Cucinotta, F.A. Track-structure codes in radiation research. Radiat. Meas. 2006, 41, 1052–1074. [Google Scholar] [CrossRef]
  36. Goorley, T.; James, M.; Booth, T.; Brown, F.; Bull, J.; Cox, L.J.; Durkee, J.; Elson, J.; Fensin, M.; Forster, R.A.; et al. Initial MCNP6 Release Overview. Nucl. Technol. 2012, 180, 298–315. [Google Scholar] [CrossRef]
  37. Kawrakow, I.; Rogers, D.W.O. The EGSnrc Code System: Monte Carlo Simulation of Electron and Photon Transport; National Research Council Canada: Ottawa, ON, Canada, January 2000. [Google Scholar]
  38. Battistoni, G.; Bauer, J.; Boehlen, T.T.; Cerutti, F.; Chin, M.P.W.; Dos Santos Augusto, R.; Ferrari, A.; Ortega, P.G.; Kozłowska, W.; Magro, G.; et al. The FLUKA Code: An Accurate Simulation Tool for Particle Therapy. Front. Oncol. 2016, 6. [Google Scholar] [CrossRef] [PubMed]
  39. Nikjoo, H.; Emfietzoglou, D.; Liamsuwan, T.; Taleei, R.; Liljequist, D.; Uehara, S. Radiation track, DNA damage and response—A review. Reports Prog. Phys. 2016, 79. [Google Scholar] [CrossRef] [PubMed]
  40. Emfietzoglou, D.; Papamichael, G.; Nikjoo, H. Monte Carlo Electron Track Structure Calculations in Liquid Water Using a New Model Dielectric Response Function. Radiat. Res. 2017, 188, 355–368. [Google Scholar] [CrossRef] [PubMed]
  41. Plante, I.; Ponomarev, A.L.; Cucinotta, F.A. Calculation of the energy deposition in nanovolumes by protons and HZE particles: Geometric patterns of initial distributions of DNA repair foci. Phys. Med. Biol. 2013, 58, 6393–6405. [Google Scholar] [CrossRef] [PubMed]
  42. Wälzlein, C.; Krämer, M.; Scifoni, E.; Durante, M. Low-energy electron transport in non-uniform media. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. Atoms 2014, 320, 75–82. [Google Scholar]
  43. Sato, T.; Iwamoto, Y.; Hashimoto, S.; Ogawa, T.; Furuta, T.; Abe, S.I.; Kai, T.; Tsai, P.E.; Matsuda, N.; Iwase, H.; et al. Features of Particle and Heavy Ion Transport code System (PHITS) version 3.02. J. Nucl. Sci. Technol. 2018, 55, 684–690. [Google Scholar] [CrossRef]
  44. Friedland, W.; Dingfelder, M.; Kundrat, P.; Jacob, P. Track structures, DNA targets and radiation effects in the biophysical Monte Carlo simulation code PARTRAC. Mutat. Res.—Fundam. Mol. Mech. Mutagen. 2011, 711, 28–40. [Google Scholar] [CrossRef]
  45. Agostinelli, S.; Allison, J.; Amako, K.; Apostolakis, J.; Araujo, H.; Arce, P.; Asai, M.; Axen, D.; Banerjee, S.; Barrand, G.; et al. GEANT4—A simulation toolkit. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 2003, 506, 250–303. [Google Scholar] [CrossRef]
  46. Kyriakou, I.; Emfietzoglou, D.; Ivanchenko, V.; Bordage, M.C.; Guatelli, S.; Lazarakis, P.; Tran, H.N.; Incerti, S. Microdosimetry of electrons in liquid water using the low-energy models of Geant4. J. Appl. Phys. 2017, 122. [Google Scholar] [CrossRef]
  47. Lazarakis, P.; Incerti, S.; Ivanchenko, V.; Kyriakou, I.; Emfietzoglou, D.; Corde, S.; Rosenfeld, A.B.; Lerch, M.; Tehei, M.; Guatelli, S. Investigation of Track Structure and Condensed History physics models for applications in radiation dosimetry on a micro and nano scale in Geant4. Biomed. Phys. Eng. Express 2018, 4, 024001. [Google Scholar] [CrossRef]
  48. Incerti, S.; Baldacchino, G.; Bernal, M.; Capra, R.; Champion, C.; Francis, Z.; Guèye, P.; Mantero, A.; Mascialino, B.; Moretto, P.; et al. The geant4-dna project. Int. J. Model. Simul. Sci. Comput. 2010, 1, 157–178. [Google Scholar] [CrossRef]
  49. Bernal, M.A.A.; Bordage, M.C.C.; Brown, J.M.C.M.C.; Davídková, M.; Delage, E.; El Bitar, Z.; Enger, S.A.A.; Francis, Z.; Guatelli, S.; Ivanchenko, V.N.N.; et al. Track structure modeling in liquid water: A review of the Geant4-DNA very low energy extension of the Geant4 Monte Carlo simulation toolkit. Phys. Med. 2015, 31, 1–14. [Google Scholar] [CrossRef]
  50. Perl, J.; Shin, J.; Schümann, J.; Faddegon, B.; Paganetti, H. TOPAS: An innovative proton Monte Carlo platform for research and clinical applications. Med. Phys. 2012, 39, 6818–6837. [Google Scholar] [CrossRef] [PubMed]
  51. McNamara, A.; Geng, C.; Turner, R.; Mendez, J.R.; Perl, J.; Held, K.; Faddegon, B.; Paganetti, H.; Schuemann, J. Validation of the radiobiology toolkit TOPAS-nBio in simple DNA geometries. Phys. Med. 2017, 33, 207–215. [Google Scholar] [CrossRef] [PubMed]
  52. Bug, M.U.; Yong Baek, W.; Rabus, H.; Villagrasa, C.; Meylan, S.; Rosenfeld, A.B. An electron-impact cross section data set (10 eV–1 keV) of DNA constituents based on consistent experimental data: A requisite for Monte Carlo simulations. Radiat. Phys. Chem. 2017, 130, 459–479. [Google Scholar] [CrossRef]
  53. Francis, Z.; El Bitar, Z.; Incerti, S.; Bernal, M.A.; Karamitros, M.; Tran, H.N. Calculation of lineal energies for water and DNA bases using the Rudd model cross sections integrated within the Geant4-DNA processes. J. Appl. Phys. 2017, 122. [Google Scholar] [CrossRef]
  54. Karamitros, M.; Luan, S.; Bernal, M.A.; Allison, J.; Baldacchino, G.; Davidkova, M.; Francis, Z.; Friedland, W.; Ivantchenko, V.; Ivantchenko, A.; et al. Diffusion-controlled reactions modeling in Geant4-DNA. J. Comput. Phys. 2014, 274, 841–882. [Google Scholar] [CrossRef]
  55. Ramos-Méndez, J.; Perl, J.; Schuemann, J.; McNamara, A.; Paganetti, H.; Faddegon, B. Monte Carlo simulation of chemistry following radiolysis with TOPAS-nBio. Phys. Med. Biol. 2018, 63. [Google Scholar] [CrossRef]
  56. Boscolo, D.; Krämer, M.; Durante, M.; Fuss, M.C.; Scifoni, E. TRAX-CHEM: A pre-chemical and chemical stage extension of the particle track structure code TRAX in water targets. Chem. Phys. Lett. 2018, 698, 11–18. [Google Scholar] [CrossRef]
  57. McNamara, A.L.; Ramos-Méndez, J.; Perl, J.; Held, K.; Dominguez, N.; Moreno, E.; Henthorn, N.T.; Kirkby, K.J.; Meylan, S.; Villagrasa, C.; et al. Geometrical structures for radiation biology research as implemented in the TOPAS-nBio toolkit. Phys. Med. Biol. 2018, 63. [Google Scholar] [CrossRef]
  58. Henthorn, N.T.; Warmenhoven, J.W.; Sotiropoulos, M.; Mackay, R.I.; Kirkby, K.J.; Merchant, M.J. Nanodosimetric Simulation of Direct Ion-Induced DNA Damage Using Different Chromatin Geometry Models. Radiat. Res. 2017, 188, 770–783. [Google Scholar] [CrossRef]
  59. Delage, E.; Pham, Q.T.; Karamitros, M.; Payno, H.; Stepan, V.; Incerti, S.; Maigne, L.; Perrot, Y. PDB4DNA: Implementation of DNA geometry from the Protein Data Bank (PDB) description for Geant4-DNA Monte-Carlo simulations. Comput. Phys. Commun. 2015, 192, 282–288. [Google Scholar] [CrossRef]
  60. Semenenko, V.A.; Stewart, R.D. A Fast Monte Carlo Algorithm to Simulate the Spectrum of DNA Damages Formed by Ionizing Radiation. Radiat. Res. 2004, 161, 451–457. [Google Scholar] [CrossRef]
  61. Friedland, W.; Jacob, P.; Bernhardt, P.; Paretzke, H.G.; Dingfelder, M. Simulation of DNA damage after proton irradiation. Radiat Res 2003, 159, 401–410. [Google Scholar] [CrossRef]
  62. Lampe, N.; Karamitros, M.; Breton, V.; Brown, J.M.C.; Kyriakou, I.; Sakata, D.; Sarramia, D.; Incerti, S. Mechanistic DNA damage simulations in Geant4-DNA part 1: A parameter study in a simplified geometry. Phys. Med. 2018, 48, 135–145. [Google Scholar] [CrossRef]
  63. Nikjoo, H.; Girard, P. A model of the cell nucleus for DNA damage calculations. Int. J. Radiat. Biol. 2012, 88, 87–97. [Google Scholar] [CrossRef]
  64. Friedland, W.; Schmitt, E.; Kundrat, P.; Dingfelder, M.; Baiocco, G.; Barbieri, S.; Ottolenghi, A. Comprehensive track-structure based evaluation of DNA damage by light ions from radiotherapy-relevant energies down to stopping. Sci. Rep. 2017, 7, 1–15. [Google Scholar] [CrossRef]
  65. Bernal, M.A.; Sikansi, D.; Cavalcante, F.; Incerti, S.; Champion, C.; Ivanchenko, V.; Francis, Z. An atomistic geometrical model of the B-DNA configuration for DNA-radiation interaction simulations. Comput. Phys. Commun. 2013, 184, 2840–2847. [Google Scholar] [CrossRef]
  66. Meylan, S.; Incerti, S.; Karamitros, M.; Tang, N.; Bueno, M.; Clairand, I.; Villagrasa, C. Simulation of early DNA damage after the irradiation of a fibroblast cell nucleus using Geant4-DNA. Sci. Rep. 2017, 7, 11923. [Google Scholar] [CrossRef]
  67. Prise, K.M.; Folkard, M.; Michael, B.D.; Vojnovic, B.; Brocklehurst, B.; Hopkirk, A.; Munro, I.H. Critical energies for SSB and DSB induction in plasmid DNA by low-energy photons: Action spectra for strand-break induction in plasmid DNA irradiated in vacuum. Int. J. Radiat. Biol. 2000, 76, 881–890. [Google Scholar]
  68. Schuemann, J.; McNamara, A.L.; Warmenhoven, J.W.; Henthorn, N.T.; Kirkby, K.; Merchant, M.J.; Ingram, S.; Paganetti, H.; Held, K.; Ramos-Mendez, J.; et al. A new Standard DNA Damage (SDD) data format. Radiat. Res. 2019, 191, 76–92. [Google Scholar] [CrossRef]
  69. Rothkamm, K.; Löbrich, M. Evidence for a lack of DNA double-strand break repair in human cells exposed to very low X-ray doses. Proc. Natl. Acad. Sci. USA 2003, 100, 5057–5062. [Google Scholar] [CrossRef]
  70. Siddiqui, M.S.; François, M.; Fenech, M.F.; Leifert, W.R. Persistent γH2AX: A promising molecular marker of DNA damage and aging. Mutat. Res.—Rev. Mutat. Res. 2015, 766, 1–19. [Google Scholar] [CrossRef]
  71. Noda, A.; Hirai, Y.; Hamasaki, K.; Mitani, H.; Nakamura, N.; Kodama, Y. Unrepairable DNA double-strand breaks that are generated by ionising radiation determine the fate of normal human cells. J. Cell Sci. 2012, 125, 5280–5287. [Google Scholar] [CrossRef]
  72. Cornforth, M.; Bedford, J. A quantitative comparison of potentially lethal damage repair and the rejoining of interphase chromosome breaks in low passage normal human fibroblasts. Radiat. Res. 1987, 111, 385–405. [Google Scholar] [CrossRef]
  73. Chen, A.M.; Lucas, J.N.; Simpson, P.J.; Griffin, C.S.; Savage, J.R.K.; Brenner, D.J.; Hlatky, L.R.; Sachs, R.K. Computer Simulation of Data on Chromosome Aberrations Produced by X Rays or Alpha Particles and Detected by Fluorescence In Situ Hybridization. Radiat. Res. 1997, 148, S93. [Google Scholar] [CrossRef]
  74. Sachs, R.; Chen, A.; Brenner, D. Review: Proximity effects in the production of chromosome aberrations by ionizing radiation. Int. J. Radiat. Biol. 1997, 71, 1–19. [Google Scholar]
  75. Tello Cajiao, J.J.; Carante, M.P.; Bernal Rodriguez, M.A.; Ballarini, F. Proximity effects in chromosome aberration induction by low-LET ionizing radiation. DNA Repair 2017, 58, 38–46. [Google Scholar] [CrossRef]
  76. McMahon, S.J.; Schuemann, J.; Paganetti, H.; Prise, K.M. Mechanistic Modelling of DNA Repair and Cellular Survival Following Radiation-Induced DNA Damage. Sci. Rep. 2016, 6, 33290. [Google Scholar] [CrossRef]
  77. Edwards, A.A.; Moiseenko, V.V.; Nikjoo, H. On the mechanism of the formation of chromosomal aberrations by ionising radiation. Radiat. Environ. Biophys. 1996, 35, 25–30. [Google Scholar] [CrossRef]
  78. McMahon, S.J.; McNamara, A.L.; Schuemann, J.; Paganetti, H.; Prise, K.M. A general mechanistic model enables predictions of the biological effectiveness of different qualities of radiation. Sci. Rep. 2017, 7, 10790. [Google Scholar] [CrossRef]
  79. Ballarini, F.; Carante, M.P. Chromosome aberrations and cell death by ionizing radiation: Evolution of a biophysical model. Radiat. Phys. Chem. 2016, 128, 18–25. [Google Scholar] [CrossRef]
  80. Waters, C.A.; Strande, N.T.; Wyatt, D.W.; Pryor, J.M.; Ramsden, D.A. Nonhomologous end joining: A good solution for bad ends. DNA Repair 2014, 17, 39–51. [Google Scholar] [CrossRef]
  81. Krejci, L.; Altmannova, V.; Spirek, M.; Zhao, X. Homologous recombination and its regulation. Nucleic Acids Res. 2012, 40, 5795–5818. [Google Scholar] [CrossRef]
  82. Deriano, L.; Roth, D.B. Modernizing the Nonhomologous End-Joining Repertoire: Alternative and Classical NHEJ Share the Stage. Annu. Rev. Genet. 2013, 47, 433–455. [Google Scholar] [CrossRef]
  83. Löbrich, M.; Jeggo, P. A Process of Resection-Dependent Nonhomologous End Joining Involving the Goddess Artemis. Trends Biochem. Sci. 2017, 42, 690–701. [Google Scholar] [CrossRef]
  84. Iliakis, G. Backup pathways of NHEJ in cells of higher eukaryotes: Cell cycle dependence. Radiother. Oncol. 2009, 92, 310–315. [Google Scholar] [CrossRef]
  85. Fowler, J. Is repair of DNA strand break damage from ionizing radiation second-order rather than first-order? A simpler explanation of apparently multiexponential repair. Radiat. Res. 1999, 136, 124–136. [Google Scholar] [CrossRef]
  86. Foray, N.; Charvet, A.M.; Duchemin, D.; Favaudon, V.; Lavalette, D. The repair rate of radiation-induced DNA damage: A stochastic interpretation based on the Gamma function. J. Theor. Biol. 2005, 236, 448–458. [Google Scholar] [CrossRef]
  87. Taleei, R.; Girard, P.M.; Sankaranarayanan, K.; Nikjoo, H. The Non-homologous End-Joining (NHEJ) Mathematical Model for the Repair of Double-Strand Breaks: II. Application to Damage Induced by Ultrasoft X Rays and Low-Energy Electrons. Radiat. Res. 2013, 179, 540–548. [Google Scholar] [CrossRef]
  88. Cucinotta, F.A.; Pluth, J.M.; Anderson, J.A.; Harper, J.V; O’Neill, P. Biochemical kinetics model of DSB repair and induction of gamma-H2AX foci by non-homologous end joining. Radiat. Res. 2008, 169, 214–222. [Google Scholar] [CrossRef]
  89. Carlson, D.J.; Stewart, R.D.; Semenenko, V.A.; Sandison, G.A. Combined Use of Monte Carlo DNA Damage Simulations and Deterministic Repair Models to Examine Putative Mechanisms of Cell Killing. Radiat. Res. 2008, 169, 447–459. [Google Scholar] [CrossRef]
  90. Friedland, W.; Jacob, P.; Kundrát, P. Stochastic Simulation of DNA Double-Strand Break Repair by Non-homologous End Joining Based on Track Structure Calculations. Radiat. Res. 2010, 173, 677–688. [Google Scholar] [CrossRef]
  91. Henthorn, N.T.; Warmenhoven, J.W.; Sotiropoulos, M.; Mackay, R.I.; Kirkby, N.F.; Kirkby, K.J.; Merchant, M.J. In Silico Non-Homologous End Joining Following Ion Induced DNA Double Strand Breaks Predicts That Repair Fidelity Depends on Break Density. Sci. Rep. 2018, 8, 2654. [Google Scholar] [CrossRef]
  92. Taleei, R.; Nikjoo, H. Biochemical DSB-repair model for mammalian cells in G1 and early S phases of the cell cycle. Mutat. Res. 2013, 756, 206–212. [Google Scholar] [CrossRef]
  93. Matos, J.; West, S.C. Holliday junction resolution: Regulation in space and time. DNA Repair 2014, 19, 176–181. [Google Scholar] [CrossRef]
  94. Seol, J.H.; Shim, E.Y.; Lee, S.E. Microhomology-mediated end joining: Good, bad and ugly. Mutat. Res.—Fundam. Mol. Mech. Mutagen. 2018, 809, 81–87. [Google Scholar] [CrossRef]
  95. Guerrero, M.; Carlone, M. Mechanistic Formulation of a Lineal-Quadratic-Linear (LQL) Model: Split-Dose Experiments and Exponentially Decaying Sources. Med. Phys. 2009, 36, 2636. [Google Scholar] [CrossRef]
  96. Zhao, L.; Mi, D.; Hu, B.; Sun, Y. A generalized target theory and its applications. Sci. Rep. 2015, 5, 1–11. [Google Scholar] [CrossRef]
  97. Bodgi, L.; Foray, N. The nucleo-shuttling of the ATM protein as a basis for a novel theory of radiation response: Resolution of the linear-quadratic model. Int. J. Radiat. Biol. 2016, 92, 117–131. [Google Scholar] [CrossRef]
  98. Carabe-Fernandez, A.; Dale, R.G.; Jones, B. The incorporation of the concept of minimum RBE (RbEmin) into the linear-quadratic model and the potential for improved radiobiological analysis of high-LET treatments. Int. J. Radiat. Biol. 2007, 83, 27–39. [Google Scholar] [CrossRef]
  99. Wedenberg, M.; Lind, B.K.; Hårdemark, B. A model for the relative biological effectiveness of protons: The tissue specific parameter α/β of photons is a predictor for the sensitivity to LET changes. Acta Oncol. 2013, 52, 580–588. [Google Scholar] [CrossRef]
  100. McNamara, A.L.; Schuemann, J.; Paganetti, H. A phenomenological relative biological effectiveness (RBE) model for proton therapy based on all published in vitro cell survival data. Phys. Med. Biol. 2015, 60, 8399–8416. [Google Scholar] [CrossRef]
  101. Butts, J.J.; Katz, R. Theory of RBE for Heavy Ion Bombardment of Dry Enzymes and Viruses. Radiat. Res. 1967, 30, 855. [Google Scholar] [CrossRef]
  102. Elsässer, T.; Krämer, M.; Scholz, M. Accuracy of the local effect model for the prediction of biologic effects of carbon ion beams in vitro and in vivo. Int. J. Radiat. Oncol. Biol. Phys. 2008, 71, 866–872. [Google Scholar] [CrossRef]
  103. Friedrich, T.; Scholz, U.; Elssser, T.; Durante, M.; Scholz, M. Calculation of the biological effects of ion beams based on the microscopic spatial damage distribution pattern. Int. J. Radiat. Biol. 2012, 88, 103–107. [Google Scholar] [CrossRef]
  104. Hawkins, R.B. A microdosimetric-kinetic model of cell death from exposure to ionizing radiation of any LET, with experimental and clinical applications. Int. J. Radiat. Biol. 1996, 69, 739–755. [Google Scholar] [CrossRef]
  105. Krämer, M.; Scholz, M. Treatment planning for heavy-ion radiotherapy: Calculation and optimization of biologically effective dose. Phys. Med. Biol. 2000, 45, 3319–3330. [Google Scholar] [CrossRef]
  106. Inaniwa, T.; Furukawa, T.; Kase, Y.; Matsufuji, N.; Toshito, T.; Matsumoto, Y.; Furusawa, Y.; Noda, K. Treatment planning for a scanned carbon beam with a modified microdosimetric kinetic model. Phys. Med. Biol. 2010, 55, 6721–6737. [Google Scholar] [CrossRef] [PubMed]
  107. Hawkins, R.B.; Inaniwa, T. A Microdosimetric-Kinetic Model for Cell Killing by Protracted Continuous Irradiation II: Brachytherapy and Biologic Effective Dose. Radiat. Res. 2014, 182, 72–82. [Google Scholar] [CrossRef]
  108. Sato, T.; Masunaga, S.I.; Kumada, H.; Hamada, N. Microdosimetric Modeling of Biological Effectiveness for Boron Neutron Capture Therapy Considering Intra- and Intercellular Heterogeneity in10B Distribution. Sci. Rep. 2018, 8, 2–11. [Google Scholar] [CrossRef]
  109. McMahon, S.J.; Hyland, W.B.; Muir, M.F.; Coulter, J.A.; Jain, S.; Butterworth, K.T.; Schettino, G.; Dickson, G.R.; Hounsell, A.R.; O’Sullivan, J.M.; et al. Biological consequences of nanoscale energy deposition near irradiated heavy atom nanoparticles. Sci. Rep. 2011, 1, 18. [Google Scholar] [CrossRef]
  110. Friedrich, T.; Durante, M.; Scholz, M. Modeling Cell Survival after Photon Irradiation Based on Double-Strand Break Clustering in Megabase Pair Chromatin Loops. Radiat. Res. 2012, 178, 385–394. [Google Scholar] [CrossRef] [PubMed]
  111. Verkhovtsev, A.; Surdutovich, E.; Solov’yov, A.V. Multiscale approach predictions for biological outcomes in ion-beam cancer therapy. Sci. Rep. 2016, 6, 1–10. [Google Scholar] [CrossRef]
  112. Frese, M.C.; Yu, V.K.; Stewart, R.D.; Carlson, D.J. A Mechanism-Based Approach to Predict the Relative Biological Effectiveness of Protons and Carbon Ions in Radiation Therapy. Int. J. Radiat. Oncol. Biol. Phys. 2011, 83, 442–450. [Google Scholar] [CrossRef]
  113. Olive, P.L. Retention of γH2AX foci as an indication of lethal DNA damage. Radiother. Oncol. 2011, 101, 18–23. [Google Scholar] [CrossRef] [PubMed]
  114. Banáth, J.P.; Klokov, D.; MacPhail, S.H.; Banuelos, C.A.; Olive, P.L. Residual γH2AX foci as an indication of lethal DNA lesions. BMC Cancer 2010, 10, 4. [Google Scholar] [CrossRef]
  115. Cohen–Jonathan, E.; Bernhard, E.J.; McKenna, W.G. How does radiation kill cells? Curr. Opin. Chem. Biol. 1999, 3, 77–83. [Google Scholar] [CrossRef]
  116. Watters, D. Molecular mechanisms of ionizing radiation-induced apoptosis. Immunol. Cell Biol. 1999, 77, 263–271. [Google Scholar] [CrossRef] [PubMed]
  117. Roos, W.P.; Kaina, B. DNA damage-induced cell death: From specific DNA lesions to the DNA damage response and apoptosis. Cancer Lett. 2013, 332, 237–248. [Google Scholar] [CrossRef] [PubMed]
  118. Zong, W.; Thompson, C.B. Necrotic death as a cell fate. Genes Dev. 2006, 20, 1–15. [Google Scholar] [CrossRef] [PubMed]
  119. Vakifahmetoglu, H.; Olsson, M.; Zhivotovsky, B. Death through a tragedy: Mitotic catastrophe. Cell Death Differ. 2008, 15, 1153–1162. [Google Scholar] [CrossRef] [PubMed]
  120. Chaurasia, M.; Bhatt, A.N.; Das, A.; Dwarakanath, B.S.; Sharma, K. Radiation-induced autophagy: Mechanisms and consequences. Free Radic. Res. 2016, 50, 273–290. [Google Scholar] [CrossRef] [PubMed]
  121. Ewald, J.A.; Desotelle, J.A.; Wilding, G.; Jarrard, D.F. Therapy-induced senescence in cancer. J. Natl. Cancer Inst. 2010, 102, 1536–1546. [Google Scholar] [CrossRef] [PubMed]
  122. Radford, I.R. Mouse lymphoma cells that undergo interphase death show markedly increased sensitivity to radiation-induced DNA double-Strand breakage as compared with cells that undergo mitotic death. Int. J. Radiat. Biol. 1991, 59, 1353–1369. [Google Scholar] [CrossRef]
  123. Langley, R.E.; Bump, E.A.; Quartuccio, S.G.; Medeiros, D.; Braunhut, S.J. Radiation-induced apoptosis in microvascular endothelial cells. Br. J. Cancer 1997, 75, 666–672. [Google Scholar] [CrossRef]
  124. Lindner, A.U.; Concannon, C.G.; Boukes, G.J.; Cannon, M.D.; Llambi, F.; Ryan, D.; Boland, K.; Kehoe, J.; McNamara, D.A.; Murray, F.; et al. Systems analysis of BCL2 protein family interactions establishes a model to predict responses to chemotherapy. Cancer Res. 2013, 73, 519–528. [Google Scholar] [CrossRef]
  125. Liu, C.; Li, C.-Y.; Yuan, F. Mathematical modeling of the phoenix rising pathway. PLoS Comput. Biol. 2014, 10, e1003461. [Google Scholar] [CrossRef]
  126. Dolan, D.W.P.; Zupanic, A.; Nelson, G.; Hall, P.; Miwa, S.; Kirkwood, T.B.L.; Shanley, D.P. Integrated Stochastic Model of DNA Damage Repair by Non-homologous End Joining and p53/p21- Mediated Early Senescence Signalling. PLoS Comput. Biol. 2015, 11, 1–19. [Google Scholar] [CrossRef] [PubMed]
  127. Zhao, Y.; Lou, I.C.; Conolly, R.B. Computational modeling of signaling pathways mediating cell cycle checkpoint control and apoptotic responses to ionizing radiation-induced DNA damage. Dose-Response 2012, 10, 251–273. [Google Scholar] [CrossRef] [PubMed]
  128. Würstle, M.L.; Zink, E.; Prehn, J.H.M.; Rehm, M. From computational modelling of the intrinsic apoptosis pathway to a systems-based analysis of chemotherapy resistance: Achievements, perspectives and challenges in systems medicine. Cell Death Dis. 2014, 5, 1–11. [Google Scholar] [CrossRef] [PubMed]
  129. West, C.M.; Davidson, S.E.; Roberts, S.A.; Hunter, R.D. The independence of intrinsic radiosensitivity as a prognostic factor for patient response to radiotherapy of carcinoma of the cervix. Br. J. Cancer 1997, 76, 1184–1190. [Google Scholar] [CrossRef]
  130. Gerweck, L.E.; Zaidi, S.T.; Zietman, A. Multivariate determinants of radiocurability 1: Prediction of single fraction tumor control doses. Int. J. Radiat. Oncol. Biol. Phys. 1994, 29, 57–66. [Google Scholar] [CrossRef]
  131. Gerweck, L.E.; Vijayappa, S.; Kurimasa, A.; Ogawa, K.; Chen, D.J. Tumor cell radiosensitivity is a major determinant of tumor response to radiation. Cancer Res. 2006, 66, 8352–8355. [Google Scholar] [CrossRef] [PubMed]
  132. Webb, S.; Nahum, A.E. A model for calculating tumour control probability in radiotherapy including the effects of inhomogeneous distributions of dose and clonogenic cell density. Phys. Med. Biol. 1993, 38, 653–666. [Google Scholar] [CrossRef] [PubMed]
  133. Gay, H.A.; Niemierko, A. A free program for calculating EUD-based NTCP and TCP in external beam radiotherapy. Phys. Med. 2007, 23, 115–125. [Google Scholar] [CrossRef] [PubMed]
  134. Niemierko, A. Reporting and analyzing dose distributions: A concept of equivalent uniform dose. Med. Phys. 1997, 24, 103–110. [Google Scholar] [CrossRef] [PubMed]
  135. Lyman, J.T. Complication probability as assessed from dose-volume histograms. Radiat. Res. Suppl. 1985, 8, S13–S19. [Google Scholar] [CrossRef] [PubMed]
  136. Burman, C.; Kutcher, G.J.; Emami, B.; Goitein, M. Fitting of normal tissue tolerance data to an analytic function. Int. J. Radiat. Oncol. Biol. Phys. 1991, 21, 123–135. [Google Scholar] [CrossRef]
  137. Källman, P.; Ågren, A.; Brahme, A. Tumour and normal tissue responses to fractionated non-uniform dose delivery. Int. J. Radiat. Biol. 1992, 62, 249–262. [Google Scholar] [CrossRef] [PubMed]
  138. Zaider, M.; Zelefsky, M.J.; Hanin, G.; Li, X.A.; Wang, J.Z.; Stewart, R.D. Tumour control probability: A formulation applicable to any temporal protocol of dose delivery. Phys. Med. Biol. 2000, 45, 279. [Google Scholar] [CrossRef] [PubMed]
  139. Jeong, J.; Shoghi, K.I.; Deasy, J.O. Modelling the interplay between hypoxia and proliferation in radiotherapy tumour response. Phys. Med. Biol. 2013, 58, 4897–4919. [Google Scholar] [CrossRef] [PubMed]
  140. Rutkowska, E.; Baker, C.; Nahum, A. Mechanistic simulation of normal-tissue damage in radiotherapy—Implications for dose-volume analyses. Phys. Med. Biol. 2010, 55, 2121–2136. [Google Scholar] [CrossRef] [PubMed]
  141. Little, M.P.; Heidenreich, W.F.; Moolgavkar, S.H.; Schöllnberger, H.; Thomas, D.C. Systems biological and mechanistic modelling of radiation-induced cancer. Radiat. Environ. Biophys. 2008, 47, 39–47. [Google Scholar] [CrossRef] [PubMed]
  142. Shuryak, I.; Brenner, D.J.; Ullrich, R.L. Radiation-Induced carcinogenesis: Mechanistically based differences between Gamma-Rays and neutrons, and interactions with DMBA. PLoS ONE 2011, 6. [Google Scholar] [CrossRef]
  143. Rühm, W.; Eidemüller, M.; Kaiser, J. Biologically-based mechanistic models of radiation-related carcinogenesis applied to epidemiological data. Int. J. Radiat. Biol. 2017, 93, 1093–1117. [Google Scholar] [CrossRef]
  144. Schneider, U. Mechanistic model of radiation-induced cancer after fractionated radiotherapy using the linear-quadratic formula. Med. Phys. 2009, 36, 1138–1143. [Google Scholar] [CrossRef]
  145. Prise, K.M.; O’Sullivan, J.M. Radiation-induced bystander signalling in cancer therapy. Nat. Rev. Cancer 2009, 9, 351–360. [Google Scholar] [CrossRef]
  146. Ebert, M.; Suchowerska, N.; Jackson, M.; McKenzie, D.R. A mathematical framework for separating the direct and bystander components of cellular radiation response. Acta Oncol. 2010, 49, 1334–1343. [Google Scholar] [CrossRef] [PubMed]
  147. McMahon, S.J.; Butterworth, K.T.; McGarry, C.K.; Trainor, C.; O’Sullivan, J.M.; Hounsell, A.R.; Prise, K.M. A Computational Model of Cellular Response to Modulated Radiation Fields. Int. J. Radiat. Oncol. Biol. Phys. 2012, 84, 250–256. [Google Scholar] [CrossRef] [PubMed]
  148. McMahon, S.J.; Butterworth, K.T.; Trainor, C.; McGarry, C.K.; O’Sullivan, J.M.; Schettino, G.; Hounsell, A.R.; Prise, K.M. A Kinetic-Based Model of Radiation-Induced Intercellular Signalling. PLoS ONE 2013, 8, e54526. [Google Scholar] [CrossRef] [PubMed]
  149. Demaria, S.; Formenti, S.C. Radiation as an immunological adjuvant: Current evidence on dose and fractionation. Front. Oncol. 2012, 2, 1–7. [Google Scholar] [CrossRef] [PubMed]
  150. Gandhi, S.J.; Minn, A.J.; Vonderheide, R.H.; Wherry, E.J.; Hahn, S.M.; Maity, A. Awakening the immune system with radiation: Optimal dose and fractionation. Cancer Lett. 2015, 368, 185–190. [Google Scholar] [CrossRef] [PubMed]
  151. Serre, R.; Benzekry, S.; Padovani, L.; Meille, C.; Andre, N.; Ciccolini, J.; Barlesi, F.; Muracciole, X.; Barbolosi, D. Mathematical modeling of cancer immunotherapy and its synergy with radiotherapy. Cancer Res. 2016, 76, 4931–4940. [Google Scholar] [CrossRef] [PubMed]
  152. Bristow, R.G.; Hill, R.P. Comparison between in vitro radiosensitivity and in vivo radioresponse in murine tumor cell lines II: In vivo radioresponse following fractionated treatment and in vitro/in vivo correlations. Int. J. Radiat. Oncol. Biol. Phys. 1990, 18, 331–345. [Google Scholar] [CrossRef]
  153. West, C.M.; Davidson, S.E.; Roberts, S.A.; Hunter, R.D. Intrinsic radiosensitivity and prediction of patient response to radiotherapy for carcinoma of the cervix. Br. J. Cancer 1993, 68, 819–823. [Google Scholar] [CrossRef]
  154. Fertil, B.; Malaise, E.P. Inherent cellular radiosensitivity as a basic concept for human tumor radiotherapy. Int. J. Radiat. Oncol. Biol. Phys. 1981, 7, 621–629. [Google Scholar] [CrossRef]
  155. Kelland, L.R.; Tonkin, K.S.; Steel, G.G. A comparison of the in vivo and in vitro radiation response of three human cervix carcinomas. Radiother. Oncol. 1989, 16, 55–63. [Google Scholar] [CrossRef]
  156. Torres-Roca, J.F.; Stevens, C.W. Predicting Response to Clinical Radiotherapy: Past, Present, and Future Directions. Cancer Control 2008, 15, 151–156. [Google Scholar] [CrossRef] [PubMed]
  157. Scott, J.G.; Berglund, A.; Schell, M.J.; Mihaylov, I.; Fulp, W.J.; Yue, B.; Welsh, E.; Caudell, J.J.; Ahmed, K.; Strom, T.S.; et al. A genome-based model for adjusting radiotherapy dose (GARD): A retrospective, cohort-based study. Lancet Oncol. 2017, 2045, 202–211. [Google Scholar] [CrossRef]
  158. Torres-Roca, J.F.; Eschrich, S.; Zhao, H.; Bloom, G.; Sung, J.; McCarthy, S.; Cantor, A.B.; Scuto, A.; Li, C.; Zhang, S.; et al. Prediction of Radiation Sensitivity Using a Gene Expression Classifier. Cancer Res. 2005, 65, 7169–7176. [Google Scholar] [CrossRef] [PubMed]
  159. Buffa, F.M.; Harris, A.L.; West, C.M.; Miller, C.J. Large meta-analysis of multiple cancers reveals a common, compact and highly prognostic hypoxia metagene. Br. J. Cancer 2010, 102, 428–435. [Google Scholar] [CrossRef] [PubMed]
  160. Hall, J.S.; Iype, R.; Senra, J.; Taylor, J.; Armenoult, L.; Oguejiofor, K.; Li, Y.; Stratford, I.; Stern, P.L.; O’Connor, M.J.; et al. Investigation of Radiosensitivity Gene Signatures in Cancer Cell Lines. PLoS ONE 2014, 9, e86329. [Google Scholar] [CrossRef] [PubMed]
  161. Forker, L.J.; Choudhury, A.; Kiltie, A.E. Biomarkers of Tumour Radiosensitivity and Predicting Benefit from Radiotherapy. Clin. Oncol. 2015, 27, 561–569. [Google Scholar] [CrossRef]
  162. Paganetti, H. Relative biological effectiveness (RBE) values for proton beam therapy. Variations as a function of biological endpoint, dose, and linear energy transfer. Phys. Med. Biol. 2014, 59, R419–R472. [Google Scholar] [CrossRef]
  163. Dungey, F.A.; Caldecott, K.W.; Chalmers, A.J. Enhanced radiosensitisation of human glioma cells by combining inhibition of PARP with inhibition of Hsp90. Mol. Cancer Ther. 2010, 8, 2243–2254. [Google Scholar] [CrossRef]
  164. Morgan, M.A.; Parsels, L.A.; Zhao, L.; Parsels, J.D.; Davis, M.A.; Hassan, M.C.; Arumugarajah, S.; Hylander-Gans, L.; Morosini, D.; Simeone, D.M.; et al. Mechanism of radiosensitization by the Chk1/2 inhibitor AZD7762 involves abrogation of the G2checkpoint and inhibition of homologous recombinational DNA repair. Cancer Res. 2010, 70, 4972–4981. [Google Scholar] [CrossRef]
  165. Dunne, V.; Ghita, M.; Small, D.M.; Coffey, C.B.M.; Weldon, S.; Taggart, C.C.; Osman, S.O.; McGarry, C.K.; Prise, K.M.; Hanna, G.G.; et al. Inhibition of ataxia telangiectasia related-3 (ATR) improves therapeutic index in preclinical models of non-small cell lung cancer (NSCLC) radiotherapy. Radiother. Oncol. 2017, 124, 475–481. [Google Scholar] [CrossRef]
  166. van der Schaaf, A.; Langendijk, J.A.; Fiorino, C.; Rancati, T. Embracing Phenomenological Approaches to Normal Tissue Complication Probability Modeling: A Question of Method. Int. J. Radiat. Oncol. 2015, 91, 468–471. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Left: Illustration of linear-quadratic dose response model. This postulates that cell survival has a “one-hit” linear component, and a “two-hit” quadratic component, which contribute varying amounts with varying doses. While simple, this model has been shown to encapsulate many aspects of radiation responses. Right: As an illustration, it provides an intuitive explanation for the benefits of fractionation, where in low α / β tissues the “shoulder” is repeated, giving rise to significant sparing, while in high α / β tissues, much less sparing is seen. However, this model has only limited scope for mechanistic applications.
Figure 1. Left: Illustration of linear-quadratic dose response model. This postulates that cell survival has a “one-hit” linear component, and a “two-hit” quadratic component, which contribute varying amounts with varying doses. While simple, this model has been shown to encapsulate many aspects of radiation responses. Right: As an illustration, it provides an intuitive explanation for the benefits of fractionation, where in low α / β tissues the “shoulder” is repeated, giving rise to significant sparing, while in high α / β tissues, much less sparing is seen. However, this model has only limited scope for mechanistic applications.
Cancers 11 00205 g001
Figure 2. Illustration of stages of a comprehensive mechanistic response model, and approximate associated timescale. Early responses are dominated by physical interactions, models of which generate distributions of reactive chemical species, and direct DSBs generated in typically less than a nanosecond. These chemical species may also be included in a chemical interaction model, which generates predictions of indirect DSBs, over a scale of 1 ms or less. Responses at the cell level then take place over hours to days, as cells repair this damage and potentially respond by arresting, dying, or seeing other changes in function. Finally, tissue-level responses integrate such responses over a scale of hours to years, giving rise to clinically observable effects. Mechanistic modelling activity in radiotherapy has primarily focused early in this timescale, on physico-chemical and early cellular responses, with relatively little modelling of whole-tissue responses. This is in part due to the significant cumulative complexity in describing all aspects of these responses in a single system.
Figure 2. Illustration of stages of a comprehensive mechanistic response model, and approximate associated timescale. Early responses are dominated by physical interactions, models of which generate distributions of reactive chemical species, and direct DSBs generated in typically less than a nanosecond. These chemical species may also be included in a chemical interaction model, which generates predictions of indirect DSBs, over a scale of 1 ms or less. Responses at the cell level then take place over hours to days, as cells repair this damage and potentially respond by arresting, dying, or seeing other changes in function. Finally, tissue-level responses integrate such responses over a scale of hours to years, giving rise to clinically observable effects. Mechanistic modelling activity in radiotherapy has primarily focused early in this timescale, on physico-chemical and early cellular responses, with relatively little modelling of whole-tissue responses. This is in part due to the significant cumulative complexity in describing all aspects of these responses in a single system.
Cancers 11 00205 g002
Figure 3. Illustration of track structure and energy deposition distributions for a 10 MeV proton (~5 keV/µm LET, blue, top) and 200 MeV carbon ion (~100 keV/µm LET, red, bottom), together with a schematic representation of DNA for scale. For both particles, the primary particle track (lines) deposits only a small fraction of the total energy, with large amounts of energy being deposited in the vicinity of the track core by low energy secondary electrons, with individual energy depositions here shown as points. It can be seen that while the proton deposits energy relatively sparsely, the carbon ion deposits energy extremely densely along its track, potentially causing multiple closely spaced damages within any DNA it traverses.
Figure 3. Illustration of track structure and energy deposition distributions for a 10 MeV proton (~5 keV/µm LET, blue, top) and 200 MeV carbon ion (~100 keV/µm LET, red, bottom), together with a schematic representation of DNA for scale. For both particles, the primary particle track (lines) deposits only a small fraction of the total energy, with large amounts of energy being deposited in the vicinity of the track core by low energy secondary electrons, with individual energy depositions here shown as points. It can be seen that while the proton deposits energy relatively sparsely, the carbon ion deposits energy extremely densely along its track, potentially causing multiple closely spaced damages within any DNA it traverses.
Cancers 11 00205 g003
Figure 4. Illustration of double-strand break repair pathways, with DNA sequence progression schematically illustrated by color scale, to highlight sequence alterations. Ionizing radiation induced double strand breaks (top) can be repaired by one of three pathways—Non-Homologous End Joining (NHEJ, left), Homologous Recombination (HR, center), or Microhomology-Mediated End Joining (MMEJ, right). These pathways differ in their efficiency, fidelity, and activity throughout the cell cycle. NHEJ has been the subject of extensive mechanistic modelling in radiotherapy, while relatively few studies have been published of the other pathways.
Figure 4. Illustration of double-strand break repair pathways, with DNA sequence progression schematically illustrated by color scale, to highlight sequence alterations. Ionizing radiation induced double strand breaks (top) can be repaired by one of three pathways—Non-Homologous End Joining (NHEJ, left), Homologous Recombination (HR, center), or Microhomology-Mediated End Joining (MMEJ, right). These pathways differ in their efficiency, fidelity, and activity throughout the cell cycle. NHEJ has been the subject of extensive mechanistic modelling in radiotherapy, while relatively few studies have been published of the other pathways.
Cancers 11 00205 g004
Table 1. List of Monte Carlo models of DNA damage, which combine various track structure physics codes with a variety of DNA models to generate predictions of different classes of DNA damage. A broad range of underlying codes and models have been used to simulate damage from a range of radiation types.
Table 1. List of Monte Carlo models of DNA damage, which combine various track structure physics codes with a variety of DNA models to generate predictions of different classes of DNA damage. A broad range of underlying codes and models have been used to simulate damage from a range of radiation types.
PaperPhysics CodeDNA ModelEndpoints
Nikjoo et al. [39,63]KURBUC/PITSWhole nucleus containing chromatin fibers arranged in hierarchy of spherical volumesDSBs, SSBs, base damages
Friedland et al. [44,64]PARTRACWhole nucleus containing model chromatin fiber random walkDSBs, SSBs, DNA fragment sizes
Bernal et al. [65]Geant4-DNAAtomistic DNA segment model in cubeDSBs, SSBs
Plante et al. [41]RITRACKSFlexible polymer chain chromosome model in nucleusDSBs
Meylan et al. [66]Geant4-DNADNAFabric-based nucleus modelDSBs, SSBs
McNamara et al. [51,57]Geant4-DNA/TOPAS-nBioMultiple DNA structures, fractal walk nucleusDSBs, SSBs

Share and Cite

MDPI and ACS Style

McMahon, S.J.; Prise, K.M. Mechanistic Modelling of Radiation Responses. Cancers 2019, 11, 205. https://doi.org/10.3390/cancers11020205

AMA Style

McMahon SJ, Prise KM. Mechanistic Modelling of Radiation Responses. Cancers. 2019; 11(2):205. https://doi.org/10.3390/cancers11020205

Chicago/Turabian Style

McMahon, Stephen J., and Kevin M. Prise. 2019. "Mechanistic Modelling of Radiation Responses" Cancers 11, no. 2: 205. https://doi.org/10.3390/cancers11020205

APA Style

McMahon, S. J., & Prise, K. M. (2019). Mechanistic Modelling of Radiation Responses. Cancers, 11(2), 205. https://doi.org/10.3390/cancers11020205

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop