Formulation of the multi-hit model with a non-Poisson distribution of hits

On the Emergence of the Deviation from a Poisson Law in Stochastic Mathematical Models for Radiation-Induced DNA Damage: A System Size Expansion

In this paper, we study the system size expansion of a stochastic model for radiation-induced DNA damage kinetics and repair. In particular, we characterize both the macroscopic deterministic limit and the fluctuation around it. We further show that such fluctuations are Gaussian-distributed. In deriving such results, we provide further insights into the relationship between stochastic and deterministic mathematical models for radiation-induced DNA damage repair. Specifically, we demonstrate how the governing deterministic equations commonly employed in the field arise naturally within the stochastic framework as a macroscopic limit. Additionally, by examining the fluctuations around this macroscopic limit, we uncover deviations from a Poissonian behavior driven by interactions and clustering among DNA damages. Although such behaviors have been empirically observed, our derived results represent the first rigorous derivation that incorporates these deviations from a Poissonian distribution within a mathematical model, eliminating the need for specific ad hoc corrections.

Accumulation of sublethal radiation damage and its effect on cell survival

Objective.Determine the extent of sublethal radiation damage (SRD) in a cell population that received a given dose of radiation and the impact of this damage on cell survival.Approach.We developed a novel formalism to account for accumulation of SRD with increasing dose. It is based on a very general formula for cell survival that correctly predicts the basic properties of cell survival curves, such as the transition from the linear-quadratic to a linear dependence at high doses. Using this formalism we analyzed extensive experimental data for photons, protons and heavy ions to evaluate model parameters, quantify the extent of SRD and its impact on cell survival.Main results.Significant accumulation of SRD begins at doses below 1 Gy. As dose increases, so does the number of damaged cells and the amount of SRD in individual cells. SRD buildup in a cell increases the likelihood of complex irrepairable damage. For this reason, during a dose fraction delivery, each dose increment makes cells more radiosensitive. This gradual radosensitization is evidenced by the increasing slope of survival curves observed experimentally. It continues until the fraction is delivered, unless radiosensitivity reaches its maximum first. The maximum radiosensitivity is achieved when SRD accumulated in most cells is the maximum damage they can repair. After this maximum is reached, the slope of a survival curve, logarithm of survival versus dose, becomes constant, dose independent. The survival curve becomes a straight line, as experimental data at high doses show. These processes are random. They cause large cell-to-cell variability in the extent of damage and radiosensitivity of individual cells.Significance.SRD is in effect a radiosensitizer and its accumulation is a significant factor affecting cell survival, especially at high doses. We developed a novel formalism to study this phenomena and reported pertinent data for several particle types.

Multiple levels of stochasticity accounted for in different radiation biophysical models: from physics to biology

PURPOSE

In the present paper we investigate how some stochastic effects are included in a class of radiobiological models with particular emphasis on how such randomnesses reflect into the predicted cell survival curve.

MATERIALS AND METHODS

We consider four different models, namely the Generalized Stochastic Microdosimetric Model GSM², in its original full form, the Dirac GSM² the Poisson GSM² and the Repair-Misrepair Model (RMR). While GSM² and the RMR models are known in literature, the Dirac and the Poisson GSM²  have been newly introduced in this work. We further numerically investigate via Monte Carlo simulation of four different particle beams, how the proposed stochastic approximations reflect into the predicted survival curves. To achieve these results, we consider different ion species at energies of interest for therapeutic applications, also including a mixed field scenario.

RESULTS

We show how the Dirac GSM², the Poisson GSM² and the RMR can be obtained from the GSM² under suitable approximations on the stochasticity considered. We analytically derive the cell survival curve predicted by the four models, characterizing rigorously the high and low dose limits. We further study how the theoretical findings emerge also using Monte Carlo numerical simulations.

CONCLUSIONS

We show how different models include different levels of stochasticity in the description of cellular response to radiation. This translates into different cell survival predictions depending on the radiation quality.

Cell Survival Computation via the Generalized Stochastic Microdosimetric Model (GSM2); Part I: The Theoretical Framework

The current article presents the first application of the Generalized Stochastic Microdosimetric Model (GSM2) for computing explicitly a cell survival curve. GSM2 is a general probabilistic model that predicts the kinetic evolution of DNA damages taking full advantage of a microdosimetric description of a radiation energy deposition. We show that, despite the high generality and flexibility of GSM2, an explicit form for the survival fraction curve predicted by the GSM2 is achievable. We illustrate how several correction terms typically added a posteriori in existing radiobiological models to improve the prediction accuracy, are naturally included into GSM2. Among the most relevant features of the survival curve derived from GSM2 and presented in this article, is the linear-quadratic behavior at low doses and a purely linear trend for high doses. The study also identifies and discusses the connections between GSM2 and existing cell survival models, such as the Microdosimetric Kinetic Model (MKM) and the Multi-hit model. Several approximations to predict cell survival in different irradiation regimes are also introduced to include intercellular non-Poissonian behaviors.

On calculation of the average linear energy transfer for radiobiological modelling

Applying the concept of linear energy transfer (LET) to modeling of biological effects of charged particles usually involves calculation of the average LET. To calculate this, the energy distribution of particles is characterized by either the source spectrum or fluence spectrum. Also, the average can be frequency-or dose-weighted. This makes four methods of calculating the average LET, each producing a different number. The purpose of this note is to describe which of these four methods is best suited for radiobiological modelling. We focused on data for photons (x-rays and gamma radiation) because in this case differences in the four averaging methods are most pronounced. However, our conclusions are equally applicable to photons and hadrons. We based our arguments on recently emerged Monte Carlo data that fully account for transport of electrons down to very low energies comparable to the ionization potential of water. We concluded that the frequency average LET calculated using the fluence spectrum has better predictive power than does that calculated using any of the other three options. This optimal method is not new but is different from those currently dominating research in this area.

Holistic View on Cell Survival and DNA Damage: How Model-Based Data Analysis Supports Exploration of Dynamics in Biological Systems

In this work, a method is established to calibrate a model that describes the basic dynamics of DNA damage and repair. The model can be used to extend planning for radiotherapy and hyperthermia in order to include the biological effects. In contrast to “syntactic” models (e.g., describing molecular kinetics), the model used here describes radiobiological semantics, resulting in a more powerful model but also in a far more challenging calibration. Model calibration is attempted from clonogenic assay data (doses of 0–6 Gy) and from time-resolved comet assay data obtained within 6 h after irradiation with 6 Gy. It is demonstrated that either of those two sources of information alone is insufficient for successful model calibration, and that both sources of information combined in a holistic approach are necessary to find viable model parameters. Approximate Bayesian computation (ABC) with simulated annealing is used for parameter search, revealing two aspects that are beneficial to resolving the calibration problem: (1) assessing posterior parameter distributions instead of point-estimates and (2) combining calibration runs from different assays by joining posterior distributions instead of running a single calibration run with a combined, computationally very expensive objective function.

A simple model for calculating relative biological effectiveness of X-rays and gamma radiation in cell survival

OBJECTIVES

The relative biological effectiveness (RBE) of X-rays and γ radiation increases substantially with decreasing beam energy. This trend affects the efficacy of medical applications of this type of radiation. This study was designed to develop a model based on a survey of experimental data that can reliably predict this trend.

METHODS

In our model, parameters α and β of a cell survival curve are simple functions of the frequency-average linear energy transfer (L_F) of delta electrons. The choice of these functions was guided by a microdosimetry-based model. We calculated L_F by using an innovative algorithm in which L_F is associated with only those electrons that reach a sensitive-to-radiation volume (SV) within the cell. We determined model parameters by fitting the model to 139 measured (α,β) pairs.

RESULTS

We tested nine versions of the model. The best agreement was achieved with [Formula: see text] and β being linear functions of [Formula: see text] .The estimated SV diameter was 0.1-1 µm. We also found that α, β, and the α/β ratio increased with increasing [Formula: see text] .

CONCLUSIONS

By combining an innovative method for calculating [Formula: see text] with a microdosimetric model, we developed a model that is consistent with extensive experimental data involving photon energies from 0.27 keV to 1.25 MeV.

ADVANCES IN KNOWLEDGE

We have developed a photon RBE model applicable to an energy range from ultra-soft X-rays to megaelectron volt γ radiation, including high-dose levels where the RBE cannot be calculated as the ratio of α values. In this model, the ionization density represented by [Formula: see text] determines the RBE for a given photon spectrum.

Generalized Multi-Hit Model of Radiation-Induced Cell Survival with a Closed-Form Solution: An Alternative Method for Determining Isoeffect Doses in Practical Radiotherapy

The standard linear-quadratic (LQ) model is currently the preferred model for describing the ionizing radiation-induced cell survival curves and tissue responses. And the LQ model is also widely used to calculate isoeffect doses for comparing different fractionated schemes in clinical radiotherapy. Despite its ubiquity, because the actual dose-response curve may appear linear at high doses in the semilogarithmic plot, the application of the LQ model is greatly challenged in the high-dose region, while the dose employed in stereotactic body radiotherapy (SBRT) is often in this area. Alternatively, the biophysical models of radiation-induced effects with a linear-quadratic-linear (LQL) characteristic can well fit the dose-survival curve of cells in vitro. However, most of these LQL models are phenomenological and have not fully considered the biophysical mechanism of radiation-induced damage and repair, and the fitting quality decreases in some high-dose ranges. In this work, to provide an alternative model to describe the cell survival curves in high-dose ranges and predict the biologically effective dose (BED) for SBRT, we propose a novel generalized multi-hit model with a closed-form solution by considering an upper bound on the number of lethal damages induced by radiation that can be repaired in a cell. This model has a clear biophysical basis and a simple expression, and also has the LQL characteristic under low- and high-dose approximate conditions. The experimental data fitting indicated that compared to the standard LQ model and our previously generalized target model, the current model can better fit the radiation-induced cell survival curves in the high-dose ranges (P < 0.05). The current model parameters and parameter ratios were determined from the fits in different kinds of cell lines irradiated with various dose rates and linear energy transfer (LET), which indicates that the model parameters significantly depend on the dose rate and LET. Based on the current model, we derived two equivalence formulae for the BED calculations in the low- and high-dose ranges, and then calculated the BED for the clinical data of SBRT from 17 selected studies. The correlation analysis showed that there were significant linear correlations between the BED at isocenter and planning target volume (PTV) edge calculated by this model and the LQ model (R > 0.86, P < 0.001). In conclusion, the generalized multi-hit model proposed in this work can be used as an alternative tool to handle in vitro radiation-induced cell survival curves in high-dose ranges, and calculate the in vivo BED for comparing the dose fractionation schemes in clinical radiotherapy.

Systematic microdosimetric data for protons of therapeutic energies calculated with Geant4-DNA

The purpose of this study was to generate physical data needed for microdosimetry-based models of proton RBE. Our focus was on the frequency and dose average lineal energies, y  F and y  D . We report data for proton energies from 0.1 to 100 MeV, for spherical volumes 2-103 nm in diameter. These data were calculated using Geant4-DNA Monte Carlo software. The physics implemented in Geant4-DNA has been extensively tested for this type of calculations but data on y  F and y  D for protons generated with this code have been very limited. An innovative aspect of our study is that we introduced a straightforward procedure for calculation of y  F and y  D for polyenergetic beams and presented the data in a format that simplifies these calculations. We compared our data with previous studies that used different Monte Carlo codes and with experimental data.

Target Size Variation in Microdosimetric Distributions and its Impact on the Linear-Quadratic Parameterization of Cell Survival

The linear-quadratic (LQ) parameterization of survival fraction [SF( D)] inherently assumes that all cells in a population receive the same dose ( D), albeit the distribution of specific energy z over the individual cells f( z, D) can be very wide. From these microdosimetric distributions, which are target size dependent, we estimate the size of the cellular sensitive volume by analyzing its influence on the LQ parameterization of cell survival. A Monte Carlo track structure code was used to simulate detailed tracks from a ⁶⁰Co source as well as proton and carbon ions of various energies. From these tracks, f( z, D) distributions were calculated for spherical targets with diameters ranging from 10 nm to 12 μm. A cell survival function based on f( z, D) was fitted to experimental LQ α values, revealing an intrinsic limitation that target size imposes on the usage of f( z, D) to describe the linear term of the LQ parameterization. The results indicate that such threshold volume arises naturally from the relationship between the particle's probability of no-hit and the probability of cell survival. Further analysis led to the proposal of a radiobiological property [Formula: see text], defined as the mean lineal energy corresponding to the target size that allows equivalence between the mean inactivation dose (MID) and the mean specific energy [Formula: see text]. The fact that [Formula: see text] is an increasing continuous function of target size within the range of biological targets of interest in radiobiology, ensures the uniqueness of [Formula: see text] for any radiation quality, thus, its potential usefulness in modeling. In conclusion, an accurate estimation of such threshold volumes may be useful for improving modeling of cell survival curves.

A new formalism for modelling parameters α and β of the linear-quadratic model of cell survival for hadron therapy

We propose a new formalism for calculating parameters α and β of the linear-quadratic model of cell survival. This formalism, primarily intended for calculating relative biological effectiveness (RBE) for treatment planning in hadron therapy, is based on a recently proposed microdosimetric revision of the single-target multi-hit model. The main advantage of our formalism is that it reliably produces α and β that have correct general properties with respect to their dependence on physical properties of the beam, including the asymptotic behavior for very low and high linear energy transfer (LET) beams. For example, in the case of monoenergetic beams, our formalism predicts that, as a function of LET, (a) α has a maximum and (b) the α/β ratio increases monotonically with increasing LET. No prior models reviewed in this study predict both properties (a) and (b) correctly, and therefore, these prior models are valid only within a limited LET range. We first present our formalism in a general form, for polyenergetic beams. A significant new result in this general case is that parameter β is represented as an average over the joint distribution of energies E 1 and E 2 of two particles in the beam. This result is consistent with the role of the quadratic term in the linear-quadratic model. It accounts for the two-track mechanism of cell kill, in which two particles, one after another, damage the same site in the cell nucleus. We then present simplified versions of the formalism, and discuss predicted properties of α and β. Finally, to demonstrate consistency of our formalism with experimental data, we apply it to fit two sets of experimental data: (1) α for heavy ions, covering a broad range of LETs, and (2) β for protons. In both cases, good agreement is achieved.

A novel multitarget model of radiation-induced cell killing based on the Gaussian distribution

The multitarget version of the traditional target theory based on the Poisson distribution is still used to describe the dose-survival curves of cells after ionizing radiation in radiobiology and radiotherapy. However, noting that the usual ionizing radiation damage is the result of two sequential stochastic processes, the probability distribution of the damage number per cell should follow a compound Poisson distribution, like e.g. Neyman's distribution of type A (N. A.). In consideration of that the Gaussian distribution can be considered as the approximation of the N. A. in the case of high flux, a multitarget model based on the Gaussian distribution is proposed to describe the cell inactivation effects in low linear energy transfer (LET) radiation with high dose-rate. Theoretical analysis and experimental data fitting indicate that the present theory is superior to the traditional multitarget model and similar to the Linear - Quadratic (LQ) model in describing the biological effects of low-LET radiation with high dose-rate, and the parameter ratio in the present model can be used as an alternative indicator to reflect the radiation damage and radiosensitivity of the cells.

Transgenerational latent early-life associated regulation unites environment and genetics across generations

The origin of idiopathic diseases is still poorly understood. The latent early-life associated regulation (LEARn) model unites environmental exposures and gene expression while providing a mechanistic underpinning for later-occurring disorders. We propose that this process can occur across generations via transgenerational LEARn (tLEARn). In tLEARn, each person is a 'unit' accumulating preclinical or subclinical 'hits' as in the original LEARn model. These changes can then be epigenomically passed along to offspring. Transgenerational accumulation of 'hits' determines a sporadic disease state. Few significant transgenerational hits would accompany conception or gestation of most people, but these may suffice to 'prime' someone to respond to later-life hits. Hits need not produce symptoms or microphenotypes to have a transgenerational effect. Testing tLEARn requires longitudinal approaches. A recently proposed longitudinal epigenome/envirome-wide association study would unite genetic sequence, epigenomic markers, environmental exposures, patient personal history taken at multiple time points and family history.

Thermoradiotherapy planning: Integration in routine clinical practice

Planning of combined radiotherapy and hyperthermia treatments should be performed taking the synergistic action between the two modalities into account. This work evaluates the available experimental data on cytotoxicity of combined radiotherapy and hyperthermia treatment and the requirements for integration of hyperthermia and radiotherapy treatment planning into a single planning platform. The underlying synergistic mechanisms of hyperthermia include inhibiting DNA repair, selective killing of radioresistant hypoxic tumour tissue and increased radiosensitivity by enhanced tissue perfusion. Each of these mechanisms displays different dose-effect relations, different optimal time intervals and different optimal sequences between radiotherapy and hyperthermia. Radiosensitisation can be modelled using the linear-quadratic (LQ) model to account for DNA repair inhibition by hyperthermia. In a recent study, an LQ model-based thermoradiotherapy planning (TRTP) system was used to demonstrate that dose escalation by hyperthermia is equivalent to ∼10 Gy for prostate cancer patients treated with radiotherapy. The first step for more reliable TRTP is further expansion of the data set of LQ parameters for normally oxygenated normal and tumour tissue valid over the temperature range used clinically and for the relevant time intervals between radiotherapy and hyperthermia. The next step is to model the effect of hyperthermia in hypoxic tumour cells including the physiological response to hyperthermia and the resulting reoxygenation. Thermoradiotherapy planning is feasible and a necessity for an optimal clinical application of hyperthermia combined with radiotherapy in individual patients.

Using state variables to model the response of tumour cells to radiation and heat: a novel multi-hit-repair approach

In order to overcome the limitations of the linear-quadratic model and include synergistic effects of heat and radiation, a novel radiobiological model is proposed. The model is based on a chain of cell populations which are characterized by the number of radiation induced damages (hits). Cells can shift downward along the chain by collecting hits and upward by a repair process. The repair process is governed by a repair probability which depends upon state variables used for a simplistic description of the impact of heat and radiation upon repair proteins. Based on the parameters used, populations up to 4-5 hits are relevant for the calculation of the survival. The model describes intuitively the mathematical behaviour of apoptotic and nonapoptotic cell death. Linear-quadratic-linear behaviour of the logarithmic cell survival, fractionation, and (with one exception) the dose rate dependencies are described correctly. The model covers the time gap dependence of the synergistic cell killing due to combined application of heat and radiation, but further validation of the proposed approach based on experimental data is needed. However, the model offers a work bench for testing different biological concepts of damage induction, repair, and statistical approaches for calculating the variables of state.