From DNA radiation damage to cell death: theoretical approaches

WHIRLY1-deficient chloroplasts display enhanced formation of cyclobutane pyrimidine dimers during exposure to UV-B radiation

The single-stranded DNA/RNA binding protein WHIRLY1 is a major chloroplast nucleoid-associated protein required for the compactness of nucleoids. Most nucleoids in chloroplasts of WHIRLY1-knockdown barley plants are less compact compared to nucleoids in wild-type plants. The reduced compaction leads to an enhanced optical cross-section, which may cause the plastid DNA to be a better target for damaging UV-B radiation. To investigate this hypothesis, primary foliage leaves, chloroplasts, and nuclei from wild-type and WHIRLY1-knockdown plants were exposed to experimental UV-B radiation. Thereafter, total, genomic and plastid DNA were isolated, respectively, and analyzed for the occurrence of cyclobutane pyrimidine dimers (CPDs), which is a parameter for genome stability. The results of this study revealed that WHIRLY1-deficient chloroplasts had strongly enhanced DNA damages, whereas isolated nuclei from the same plant line were not more sensitive than nuclei from the wild-type, indicating that WHIRLY1 has different functions in chloroplasts and nucleus. This supports the hypothesis that the compaction of nucleoids may provide protection against UV-B radiation.

Formulation of Time-Dependent Cell Survival with Saturable Repairability of Radiation Damage

This study aims to provide a model that compounds historically proposed ideas regarding cell survival irradiated with X rays or particles. The parameters used in this model have simple meanings and are closely related to cell death-related phenomena. The model is adaptable to a wide range of doses and dose rates and thus can consistently explain previously published cell survival data. The formulas of the model were derived by using five basic ideas: 1. "Poisson's law"; 2. "DNA affected damage"; 3. "repair"; 4. "clustered affected damage"; and 5. "saturation of reparability". The concept of affected damage is close to but not the same as the effect caused by the double-strand break (DSB). The parameters used in the formula are related to seven phenomena: 1. "linear coefficient of radiation dose"; 2. "probability of making affected damage"; 3. "cell-specific repairability", 4. "irreparable damage by adjacent affected damage"; 5. "recovery of temporally changed repairability"; 6. "recovery of simple damage which will make the affected damage"; 7. "cell division". By using the second parameter, this model includes cases where a single hit results in repairable-lethal and double-hit results in repairable-lethal. The fitting performance of the model for the experimental data was evaluated based on the Akaike information criterion, and practical results were obtained for the published experiments irradiated with a wide range of doses (up to several 10 Gy) and dose rates (0.17 Gy/h to 55.8 Gy/h). The direct association of parameters with cell death-related phenomena has made it possible to systematically fit survival data of different cell types and different radiation types by using crossover parameters.

Basic and translational research on carbon-ion radiobiology

Carbon-ion beam irradiation (IR) has evident advantages over the conventional photon beams in treating tumors. It releases enormous amount of energy in a well-defined range with insignificant scatter in surrounding tissues based on well-localized energy deposition. Over the past 28 years, more than 14,000 patients with various types of cancer have been treated by carbon ion radiotherapy (CIRT) with promising results at QST. I have provided an overview of the basic and translational research on carbon-ion radiobiology including mechanisms underlying high linear energy transfer (LET) carbon-ion IR-induced cell death (apoptosis, autophagy, senescence, mitotic catastrophe etc.) and high radiocurability produced by carbon-ion beams in combination with DNA damaging drugs or with molecular-targeted drugs, micro-RNA therapeutics and immunotherapy. Additionally, I have focused on the application of these treatment in human cancer cells, especially cancer stem cells (CSCs). Finally, I have summarized the current studies on the application of basic carbon-ion beam IR according to the cancer types and clinical outcomes.

Time-dose reciprocity mechanism for the inactivation of Escherichia coli explained by a stochastic process with two inactivation effects

There is a great demand for developing and demonstrating novel disinfection technologies for protection against various pathogenic viruses and bacteria. In this context, ultraviolet (UV) irradiation offers an effective and convenient method for the inactivation of pathogenic microorganisms. The quantitative evaluation of the efficacy of UV sterilization relies on the simple time-dose reciprocity law proposed by Bunsen-Roscoe. However, the inactivation rate constants reported in the literature vary widely, even at the same dose and wavelength of irradiation. Thus, it is likely that the physical mechanism of UV inactivation cannot be described by the simple time-dose reciprocity law but requires a secondary inactivation process, which must be identified to clarify the scientific basis. In this paper, we conducted a UV inactivation experiment with Escherichia coli at the same dose but with different irradiances and irradiation durations, varying the irradiance by two to three orders of magnitude. We showed that the efficacy of inactivation obtained by UV-light emitting diode irradiation differs significantly by one order of magnitude at the same dose but different irradiances at a fixed wavelength. To explain this, we constructed a stochastic model introducing a second inactivation rate, such as that due to reactive oxygen species (ROS) that contribute to DNA and/or protein damage, together with the fluence-based UV inactivation rate. By solving the differential equations based on this model, the efficacy of inactivation as a function of the irradiance and irradiation duration under the same UV dose conditions was clearly elucidated. The proposed model clearly shows that at least two inactivation rates are involved in UV inactivation, where the generally used UV inactivation rate does not depend on the irradiance, but the inactivation rate due to ROS does depend on the irradiance. We conclude that the UV inactivation results obtained to date were simply fitted by one inactivation rate that superimposed these two inactivation rates. The effectiveness of long-term UV irradiation at a low irradiance but the same dose provides useful information for future disinfection technologies such as the disinfection of large spaces, for example, hospital rooms using UV light, because it can reduce the radiation dose and its risk to the human body.

A computational approach to quantifying miscounting of radiation-induced double-strand break immunofluorescent foci

Immunofluorescent tagging of DNA double-strand break (DSB) markers, such as γ-H2AX and other DSB repair proteins, are powerful tools in understanding biological consequences following irradiation. However, whilst the technique is widespread, there are many uncertainties related to its ability to resolve and reliably deduce the number of foci when counting using microscopy. We present a new tool for simulating radiation-induced foci in order to evaluate microscope performance within in silico immunofluorescent images. Simulations of the DSB distributions were generated using Monte Carlo track-structure simulation. For each DSB distribution, a corresponding DNA repair process was modelled and the un-repaired DSBs were recorded at several time points. Corresponding microscopy images for both a DSB and (γ-H2AX) fluorescent marker were generated and compared for different microscopes, radiation types and doses. Statistically significant differences in miscounting were found across most of the tested scenarios. These inconsistencies were propagated through to repair kinetics where there was a perceived change between radiation-types. These changes did not reflect the underlying repair rate and were caused by inconsistencies in foci counting. We conclude that these underlying uncertainties must be considered when analysing images of DNA damage markers to ensure differences observed are real and are not caused by non-systematic miscounting.

PyFoci is a tool that simulates distributions of fluorescently labeled DNA double-strand break marker protein foci and allows the estimation of miscounting under different radiation types, doses and microscopy settings.

Hi-C implementation of genome structure for in silico models of radiation-induced DNA damage

Developments in the genome organisation field has resulted in the recent methodology to infer spatial conformations of the genome directly from experimentally measured genome contacts (Hi-C data). This provides a detailed description of both intra- and inter-chromosomal arrangements. Chromosomal intermingling is an important driver for radiation-induced DNA mis-repair. Which is a key biological endpoint of relevance to the fields of cancer therapy (radiotherapy), public health (biodosimetry) and space travel. For the first time, we leverage these methods of inferring genome organisation and couple them to nano-dosimetric radiation track structure modelling to predict quantities and distribution of DNA damage within cell-type specific geometries. These nano-dosimetric simulations are highly dependent on geometry and are benefited from the inclusion of experimentally driven chromosome conformations. We show how the changes in Hi-C contract maps impact the inferred geometries resulting in significant differences in chromosomal intermingling. We demonstrate how these differences propagate through to significant changes in the distribution of DNA damage throughout the cell nucleus, suggesting implications for DNA repair fidelity and subsequent cell fate. We suggest that differences in the geometric clustering for the chromosomes between the cell-types are a plausible factor leading to changes in cellular radiosensitivity. Furthermore, we investigate changes in cell shape, such as flattening, and show that this greatly impacts the distribution of DNA damage. This should be considered when comparing in vitro results to in vivo systems. The effect may be especially important when attempting to translate radiosensitivity measurements at the experimental in vitro level to the patient or human level.

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.

Fitting the Generalized Target Model to Cell Survival Data of Proton Radiation Reveals Dose-Dependent RBE and Inspires an Alternative Method to Estimate RBE in High-Dose Regions

The imprecise estimation of the relative biological effectiveness (RBE) of proton radiation has been one of the main challenges for further calculating the biologically effective dose in proton therapy. Since dose levels can greatly influence the proton RBE, the relationship between the two should be clarified first. In addition, since the dose-response curves are usually too complex to readily assess RBE in high-dose regions, a reliable and simple method is needed to predict the RBE of proton radiation accurately in clinically relevant doses. The standard linear-quadratic (LQ) model is widely used to determine the RBE of particles for clinical applications. However, there has been some debate over its use when modeling the cell survival curves in high-dose regions, since those survival curves usually show linear behavior in the semilogarithmic plot. By considering both cellular repair effects and indirect effects of radiation, we have proposed a generalized target model with linear-quadratic linear (LQL) characteristics. For the more accurate evaluation of proton RBE in radiotherapy, here we used this generalized target model to fit the cell survival data in V79 and C3H 10T1/2 cells exposed to proton radiation with different LETs. The fitting results show that the generalized target model works as well as the LQ model in general. Based on the fitting parameters of the generalized target model, the RBE of six given doses D_T (RBE_T) could be calculated in the corresponding cell lines with different LETs. The results show that the RBE_T gradually decreases with increased dose in both cell types. In addition, inspired by the calculation method of the maximum values of RBE (RBE_M) in the low-dose region, a novel method was proposed for estimating the RBE in the high-dose region (RBE_H) based on the slope ratio of the dose-response curves in this region. Linear regression analysis indicated a significant linear correlation between the proposed RBE_H and the RBE_T in high-dose regions, which suggests that the current method can be used as an alternative tool, which is both simple and robust, to estimate RBE in high-dose regions.

Modeling the effect of intratumoral heterogeneity of radiosensitivity on tumor response over the course of fractionated radiation therapy

Background

Standard radiobiology theory of radiation response assumes a uniform innate radiosensitivity of tumors. However, experimental data show that there is significant intratumoral heterogeneity of radiosensitivity. Therefore, a model with heterogeneity was developed and tested using existing experimental data to show the potential effects from the presence of an intratumoral distribution of radiosensitivity on radiation therapy response over a protracted radiation therapy treatment course.

Methods

The standard radiation response curve was modified to account for a distribution of radiosensitivity, and for variations in the repopulation rates of the tumor cell subpopulations. Experimental data from the literature were incorporated to determine the boundaries of the model. The proposed model was then used to show the changes in radiosensitivity of the tumor during treatment, and the effects of fraction size, α/β ratio and variation of the repopulation rates of tumor cells.

Results

In the presence of an intratumoral distribution of radiosensitivity, there is rapid selection of radiation-resistant cells over a course of fractionated radiation therapy. Standard treatment fractionation regimes result in the near-complete replacement of the initial population of sensitive cells with a population of more resistant cells. Further, as treatment progresses, the tumor becomes more resistant to further radiation treatment, making each fractional dose less efficacious. A wider initial distribution induces increased radiation resistance. Hypofractionation is more efficient in a heterogeneous tumor, with increased cell kill for biologically equivalent doses, while inducing less resistance. The model also shows that a higher growth rate in resistant cells can account for the accelerated repopulation that is seen during the clinical treatment of patients.

Conclusions

Modeling of tumor cell survival with radiosensitivity heterogeneity alters the predicted tumor response, and explains the induction of radiation resistance by radiation treatment, the development of accelerated repopulation, and the potential beneficial effects of hypofractionation. Tumor response to treatment may be better predicted by assaying for the distribution of radiosensitivity, or the extreme of the radiosensitivity, rather than measuring the initial, general radiation sensitivity of the untreated tumor.

Analysis of Radiation-Induced Chromosomal Aberrations on a Cell-by-Cell Basis after Alpha-Particle Microbeam Irradiation: Experimental Data and Simulations

There is a continued need for further clarification of various aspects of radiation-induced chromosomal aberration, including its correlation with radiation track structure. As part of the EMRP joint research project, Biologically Weighted Quantities in Radiotherapy (BioQuaRT), we performed experimental and theoretical analyses on chromosomal aberrations in Chinese hamster ovary cells (CHO-K1) exposed to α particles with final energies of 5.5 and 17.8 MeV (absorbed doses: ∼2.3 Gy and ∼1.9 Gy, respectively), which were generated by the microbeam at the Physikalisch-Technische Bundesanstalt (PTB) in Braunschweig, Germany. In line with the differences in linear energy transfer (approximately 85 keV/μm for 5.5 MeV and 36 keV/μm for 17.8 MeV α particles), the 5.5 MeV α particles were more effective than the 17.8 MeV α particles, both in terms of the percentage of aberrant cells (57% vs. 33%) and aberration frequency. The yield of total aberrations increased by a factor of ∼2, although the increase in dicentrics plus centric rings was less pronounced than in acentric fragments. The experimental data were compared with Monte Carlo simulations based on the BIophysical ANalysis of Cell death and chromosomal Aberrations model (BIANCA). This comparison allowed interpretation of the results in terms of critical DNA damage [cluster lesions (CLs)]. More specifically, the higher aberration yields observed for the 5.5 MeV α particles were explained by taking into account that, although the nucleus was traversed by fewer particles (nominally, 11 vs. 25), each particle was much more effective (by a factor of ∼3) at inducing CLs. This led to an increased yield of CLs per cell (by a factor of ∼1.4), consistent with the increased yield of total aberrations observed in the experiments.

BIANCA, a biophysical model of cell survival and chromosome damage by protons, C-ions and He-ions at energies and doses used in hadrontherapy

An upgraded version of the BIANCA II biophysical model, which describes more realistically interphase chromosome organization and the link between chromosome aberrations and cell death, was applied to V79 and AG01522 cells exposed to protons, C-ions and He-ions over a wide LET interval (0.6-502 keV µm^-1), as well as proton-irradiated U87 cells. The model assumes that (i) ionizing radiation induces DNA 'cluster lesions' (CLs), where by definition each CL produces two independent chromosome fragments; (ii) fragment (distance-dependent) mis-rejoining, or un-rejoining, produces chromosome aberrations; (iii) some aberrations lead to cell death. The CL yield, which mainly depends on radiation quality but is also modulated by the target cell, is an adjustable parameter. The fragment un-rejoining probability, f, is the second, and last, parameter. The value of f, which is assumed to depend on the cell type but not on radiation quality, was taken from previous studies, and only the CL yield was adjusted in the present work. Good agreement between simulations and experimental data was obtained, suggesting that BIANCA II is suitable for calculating the biological effectiveness of hadrontherapy beams. For both V79 and AG01522 cells, the mean number of CLs per micrometer was found to increase with LET in a linear-quadratic fashion before the over-killing region, where a less rapid increase, with a tendency to saturation, was observed. Although the over-killing region deserves further investigation, the possibility of fitting the CL yields is an important feature for hadrontherapy, because it allows performing predictions also at LET values where experimental data are not available. Finally, an approach was proposed to predict the ion-response of the cell line(s) of interest from the ion-response of a reference cell line and the photon response of both. A pilot study on proton-irradiated AG01522 and U87 cells, taking V79 cells as a reference, showed encouraging results.

Proximity effects in chromosome aberration induction: Dependence on radiation quality, cell type and dose

It is widely accepted that, in chromosome-aberration induction, the (mis-)rejoining probability of two chromosome fragments depends on their initial distance, r. However, several aspects of these "proximity effects" need to be clarified, also considering that they can vary with radiation quality, cell type and dose. A previous work performed by the BIANCA (BIophysical ANalysis of Cell death and chromosome Aberrations) biophysical model has suggested that, in human lymphocytes and fibroblasts exposed to low-LET radiation, an exponential function of the form exp(-r/r₀), which is consistent with free-end (confined) diffusion, describes proximity effects better than a Gaussian function. Herein, the investigation was extended to intermediate- and high-LET. Since the r₀ values (0.8 μm for lymphocytes and 0.7 μm for fibroblasts) were taken from the low-LET study, the results were obtained by adjusting only one model parameter, i.e. the yield of "Cluster Lesions" (CLs), where a CL was defined as a critical DNA damage producing two independent chromosome fragments. In lymphocytes, the exponential model allowed reproducing both dose-response curves for different aberrations (dicentrics, centric rings and excess acentrics), and values of F-ratio (dicentrics to centric rings) and G-ratio (interstitial deletions to centric rings). In fibroblasts, a good correspondence was found with the dose-response curves, whereas the G-ratio (and, to a lesser extent, the F-ratio) was underestimated. With increasing LET, F decreased and G increased in both cell types, supporting their role as "fingerprints" of high-LET exposure. A dose-dependence was also found at high LET, where F increased with dose and G decreased, possibly due to inter-track effects. We therefore conclude that, independent of radiation quality, in lymphocytes an exponential function can describe proximity effects at both inter- and intra-chromosomal level; on the contrary, in fibroblasts further studies (experimental and theoretical) are needed to explain the strong bias for intra-arm relative to inter-arm exchanges.

Proximity effects in chromosome aberration induction by low-LET ionizing radiation

Although chromosome aberrations are known to derive from distance-dependent mis-rejoining of chromosome fragments, evaluating whether a certain model describes such "proximity effects" better than another one is complicated by the fact that different approaches have often been tested under different conditions. Herein, a biophysical model ("BIANCA", i.e. BIophysical ANalysis of Cell death and chromosome Aberrations) was upgraded, implementing explicit chromosome-arm domains and two new models for the dependence of the rejoining probability on the fragment initial distance, r. Such probability was described either by an exponential function like exp(-r/r₀), or by a Gaussian function like exp(-r²/2σ²), where r₀ and σ were adjustable parameters. The second, and last, parameters was the yield of "Cluster Lesions" (CL), where "Cluster Lesion" defines a critical DNA damage producing two independent chromosome fragments. The model was applied to low-LET-irradiated lymphocytes (doses: 1-4Gy) and fibroblasts (1-6.1Gy). Good agreement with experimental yields of dicentrics and centric rings, and thus their ratio ("F-ratio"), was found by both the exponential model (with r₀=0.8μm for lymphocytes and 0.7μm for fibroblasts) and the Gaussian model (with σ=1.1μm for lymphocytes and 1.3μm for fibroblasts). While the former also allowed reproducing dose-responses for excess acentric fragments, the latter substantially underestimated the experimental curves. Both models provided G-ratios (ratio of acentric to centric rings) higher than those expected from randomness, although the values calculated by the Gaussian model were lower than those calculated by the exponential one. For lymphocytes the calculated G-ratios were in good agreement with the experimental ones, whereas for fibroblasts both models substantially underestimated the experimental results, which deserves further investigation. This work suggested that, although both models performed better than a step model (which previously allowed reproducing the F-ratio but underestimated the G-ratio), an exponential function describes proximity effects better than a Gaussian one.

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.

Radiosensitivity and relative biological effectiveness based on a generalized target model

By considering both cellular repair effects and indirect effects of radiation, we have generalized the traditional target model, and made it have a linear-quadratic-linear characteristic. To assess the repair capacity-dependent radiosensitivity and relative biological effectiveness (RBE), the generalized target model was used to fit the survival of human normal embryonic lung fibroblast MRC-5 cells in the G0 and G1 phases after various types of radiations. The fitting results indicate that the generalized target model works well in the dose ranges considered. The resulting calculations qualitatively show that the parameter ratio (a/V) in the model could represent the cellular repair capacity. In particular, the significant linear correlations between radiosensitivity/RBE and cellular repair capacity are observed for different slopes of the linear regression curves. These results show that the radiosensitivity and RBE depend on the cellular repair capacity and can be regulated by linear energy transfer. These analyses suggest that the ratio a/V in the generalized target model can also be used for radiation damage assessment in radiotherapy.

First proof of bismuth oxide nanoparticles as efficient radiosensitisers on highly radioresistant cancer cells

This study provides the first proof of the novel application of bismuth oxide as a radiosensitiser. It was shown that on the highly radioresistant 9L gliosarcoma cell line, bismuth oxide nanoparticles sensitise to both kilovoltage (kVp) or megavoltage (MV) X-rays radiation. 9L cells were exposed to a concentration of 50μg.mL^-1 of nanoparticle before irradiation at 125kVp and 10MV. Sensitisation enhancement ratios of 1.48 and 1.25 for 125kVp and 10MV were obtained in vitro, respectively. The radiation enhancement of the nanoparticles is postulated to be a combination of the high Z nature of the bismuth (Z=83), and the surface chemistry. Monte Carlo simulations were performed to elucidate the physical interactions between the incident radiation and the nanoparticle. The results of this work show that Bi₂O₃ nanoparticles increase the radiosensitivity of 9L gliosarcoma tumour cells for both kVp and MV energies. Monte Carlo simulations demonstrate the advantage of a platelet morphology.

Calculating Variations in Biological Effectiveness for a 62 MeV Proton Beam

A biophysical model of radiation-induced cell death and chromosome aberrations [called BIophysical ANalysis of Cell death and chromosome Aberrations (BIANCA)] was further developed and applied to therapeutic protons. The model assumes a pivotal role of DNA cluster damage, which can lead to clonogenic cell death following three main steps: (i) a DNA “cluster lesion” (CL) produces two independent chromosome fragments; (ii) fragment mis-rejoining within a threshold distance d gives rise to chromosome aberrations; (iii) certain aberration types (dicentrics, rings, and large deletions) lead to clonogenic inactivation. The yield of CLs and the probability, f, that a chromosome fragment remains un-rejoined even if other fragment(s) are present within d, were adjustable parameters. The model, implemented as a MC code providing simulated dose–responses directly comparable with experimental data, was applied to pristine and modulated Bragg peaks of the proton beam used to treat eye melanoma at INFN-LNS in Catania, Italy. Experimental survival curves for AG01522 cells exposed to the Catania beam were reproduced, supporting the model assumptions. Furthermore, cell death and chromosome aberrations at different depths along a spread-out Bragg peak (SOBP) dose profile were predicted. Both endpoints showed an increase along the plateau, and high levels of damage were found also beyond the distal dose fall-off, due to low-energy protons. Cell death and chromosome aberrations were also predicted for V79 cells, in the same irradiation scenario as that used for AG01522 cells. In line with other studies, this work indicated that assuming a constant relative biological effectiveness (RBE) along a proton SOBP may be sub-optimal. Furthermore, it provided qualitative and quantitative evaluations of the dependence of the beam effectiveness on the considered endpoint and dose. More generally, this work represents an example of therapeutic beam characterization avoiding the use of experimental RBE values, which can be source of uncertainties.

Experimental Adaptation of Rotaviruses to Tumor Cell Lines

A number of viruses show a naturally extended tropism for tumor cells whereas other viruses have been genetically modified or adapted to infect tumor cells. Oncolytic viruses have become a promising tool for treating some cancers by inducing cell lysis or immune response to tumor cells. In the present work, rotavirus strains TRF-41 (G5) (porcine), RRV (G3) (simian), UK (G6-P5) (bovine), Ym (G11-P9) (porcine), ECwt (murine), Wa (G1-P8), Wi61 (G9) and M69 (G8) (human), and five wild-type human rotavirus isolates were passaged multiple times in different human tumor cell lines and then combined in five different ways before additional multiple passages in tumor cell lines. Cell death caused by the tumor cell-adapted isolates was characterized using Hoechst, propidium iodide, 7-AAD, Annexin V, TUNEL, and anti-poly-(ADP ribose) polymerase (PARP) and -phospho-histone H2A.X antibodies. Multiple passages of the combined rotaviruses in tumor cell lines led to a successful infection of these cells, suggesting a gain-of-function by the acquisition of greater infectious capacity as compared with that of the parental rotaviruses. The electropherotype profiles suggest that unique tumor cell-adapted isolates were derived from reassortment of parental rotaviruses. Infection produced by such rotavirus isolates induced chromatin modifications compatible with apoptotic cell death.

A generalized target theory and its applications

Different radiobiological models have been proposed to estimate the cell-killing effects, which are very important in radiotherapy and radiation risk assessment. However, most applied models have their own scopes of application. In this work, by generalizing the relationship between "hit" and "survival" in traditional target theory with Yager negation operator in Fuzzy mathematics, we propose a generalized target model of radiation-induced cell inactivation that takes into account both cellular repair effects and indirect effects of radiation. The simulation results of the model and the rethinking of "the number of targets in a cell" and "the number of hits per target" suggest that it is only necessary to investigate the generalized single-hit single-target (GSHST) in the present theoretical frame. Analysis shows that the GSHST model can be reduced to the linear quadratic model and multitarget model in the low-dose and high-dose regions, respectively. The fitting results show that the GSHST model agrees well with the usual experimental observations. In addition, the present model can be used to effectively predict cellular repair capacity, radiosensitivity, target size, especially the biologically effective dose for the treatment planning in clinical applications.

The role of DNA cluster damage and chromosome aberrations in radiation-induced cell killing: a theoretical approach

The role played by DNA cluster damage and chromosome aberrations in radiation-induced cell killing was investigated, assuming that certain chromosome aberrations (dicentrics, rings and large deletions, or 'lethal aberrations') lead to clonogenic inactivation and that chromosome aberrations are due to micrometre-scale rejoining of chromosome fragments derived from DNA cluster lesions (CLs). The CL yield and the threshold distance governing fragment rejoining were left as model parameters. The model, implemented as a Monte Carlo code called BIANCA (BIophysical ANalysis of Cell death and chromosome Aberrations), provided simulated survival curves that were compared with survival data on AG1522 and V79 cells exposed to different radiation types, including heavy ions. The agreement between simulation outcomes and experimental data suggests that lethal aberrations are likely to play an important role in cell killing not only for AG1522 cells exposed to X rays, as already reported by others, but also for other radiation types and other cells. Furthermore, the results are consistent with the hypothesis that the critical DNA lesions leading to cell death and chromosome aberrations are double-strand break clusters (possibly involving the ∼1000-10 000 bp scale) and that the effects of such clusters are modulated by micrometre-scale proximity effects during DNA damage processing.

Modeling radiation-induced cell death: role of different levels of DNA damage clustering

Some open questions on the mechanisms underlying radiation-induced cell death were addressed by a biophysical model, focusing on DNA damage clustering and its consequences. DNA "cluster lesions" (CLs) were assumed to produce independent chromosome fragments that, if created within a micrometer-scale threshold distance (d), can lead to chromosome aberrations following mis-rejoining; in turn, certain aberrations (dicentrics, rings and large deletions) were assumed to lead to clonogenic cell death. The CL yield and d were the only adjustable parameters. The model, implemented as a Monte Carlo code called BIophysical ANalysis of Cell death and chromosome Aberrations (BIANCA), provided simulated survival curves that were directly compared with experimental data on human and hamster cells exposed to photons, protons, α-particles and heavier ions including carbon and iron. d = 5 μm, independent of radiation quality, and CL yields in the range ~2-20 CLs Gy(-1) cell(-1), depending on particle type and energy, led to good agreement between simulations and data. This supports the hypothesis of a pivotal role of DNA cluster damage at sub-micrometric scale, modulated by chromosome fragment mis-rejoining at micrometric scale. To investigate the features of such critical damage, the CL yields were compared with experimental or theoretical yields of DNA fragments of different sizes, focusing on the base-pair scale (related to the so-called local clustering), the kbp scale ("regional clustering") and the Mbp scale, corresponding to chromatin loops. Interestingly, the CL yields showed better agreement with kbp fragments rather than bp fragments or Mbp fragments; this suggests that also regional clustering, in addition to other clustering levels, may play an important role, possibly due to its relationship with nucleosome organization in the chromatin fiber.

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

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

Toward computer simulation of high-LET in vitro survival curves

We developed a Monte Carlo based computer program called MCSC (Monte Carlo Survival Curve) able to predict the survival fraction of cells irradiated in vitro with a broad beam of high linear energy transfer particles. Three types of cell responses are studied: the usual high dose response, the bystander effect and the low-dose hypersensitivity (HRS). The program models the broad beam irradiation and double strand break distribution following Poisson statistics. The progression of cells through the cell cycle is taken into account while the repair takes place. Input parameters are experimentally determined for A549 lung carcinoma cells irradiated with 10 and 20 keV µm(-1) protons, 115 keV µm(-1) alpha particles and for EAhy926 endothelial cells exposed to 115 keV µm(-1) alpha particles. Results of simulations are presented and compared with experimental survival curves obtained for A549 and EAhy296 cells. Results are in good agreement with experimental data for both cell lines and all irradiation protocols. The benefits of MCSC are several: the gain of time that would have been spent performing time-consuming clonogenic assays, the capacity to estimate survival fraction of cell lines not forming colonies and possibly the evaluation of radiosensitivity parameters of given individuals.

Direct action of radiation on mummified cells: modeling of computed tomography by Monte Carlo algorithms

X-ray imaging is a nondestructive and preferred method in paleopathology to reconstruct the history of ancient diseases. Sophisticated imaging technologies such as computed tomography (CT) have become common for the investigation of skeletal disorders in human remains. Researchers have investigated the impact of ionizing radiation on living cells, but never on ancient cells in dry tissue. The effects of CT exposure on ancient cells have not been examined in the past and may be important for subsequent genetic analysis. To remedy this shortcoming, we developed different Monte Carlo models to simulate X-ray irradiation on ancient cells. Effects of mummification were considered by using two sizes of cells and three different phantom tissues, which enclosed the investigated cell cluster. This cluster was positioned at the isocenter of a CT scanner model, where the cell hit probabilities P(0,1,…, n) were calculated according to the Poisson distribution. To study the impact of the dominant physics process, CT scans for X-ray spectra of 80 and 120 kVp were simulated. Comparison between normal and dry tissue phantoms revealed that the probability of unaffected cells increased by 21 % following cell shrinkage for 80 kVp, while for 120 kVp, a further increase of unaffected cells of 23 % was observed. Consequently, cell shrinkage caused by dehydration decreased the impact of X-ray radiation on mummified cells significantly. Moreover, backscattered electrons in cortical bone protected deeper-lying ancient cells from radiation damage at 80 kVp X-rays.

Modelling of cell killing due to sparsely ionizing radiation in normoxic and hypoxic conditions and an extension to high LET radiation

PURPOSE

An approach for describing cell killing with sparsely ionizing radiation in normoxic and hypoxic conditions based on the initial number of randomly distributed DNA double-strand breaks (DSB) is proposed. An extension of the model to high linear energy transfer (LET) radiation is also presented.

MATERIALS AND METHODS

The model is based on the probabilities that a given DNA giant loop has one DSB or at least two DSB. A linear combination of these two classes of damage gives the mean number of lethal lesions. When coupled with a proper modelling of the spatial distribution of DSB from ion tracks, the formalism can be used to predict cell response to high LET radiation in aerobic conditions.

RESULTS

Survival data for sparsely ionizing radiation of cell lines in normoxic/hypoxic conditions were satisfactorily fitted with the proposed parametrization. It is shown that for dose ranges up to about 10 Gy, the model describes tested experimental survival data as good as the linear-quadratic model does. The high LET extension yields a reasonable agreement with data in aerobic conditions.

CONCLUSIONS

A new survival model has been introduced that is able to describe the most relevant features of cellular dose-response postulating two damage classes.

Basics of particle therapy II: relative biological effectiveness

In the previous review, the physical aspect of heavy particles, with a focus on the carbon beam was introduced. Particle beam therapy has many potential advantages for cancer treatment without increasing severe side effects in normal tissue, these kinds of radiation have different biologic characteristics and have advantages over using conventional photon beam radiation during treatment. The relative biological effectiveness (RBE) is used for many biological, clinical endpoints among different radiation types and is the only convenient way to transfer the clinical experience in radiotherapy with photons to another type of radiation therapy. However, the RBE varies dependent on the energy of the beam, the fractionation, cell types, oxygenation status, and the biological endpoint studied. Thus this review describes the concerns about RBE related to particle beam to increase interests of the Korean radiation oncologists' society.

Chromatin dynamics during cell cycle mediate conversion of DNA damage into chromatid breaks and affect formation of chromosomal aberrations: biological and clinical significance

The formation of diverse chromosomal aberrations following irradiation and the variability in radiosensitivity at different cell-cycle stages remain a long standing controversy, probably because most of the studies have focused on elucidating the enzymatic mechanisms involved using simple DNA substrates. Yet, recognition, processing and repair of DNA damage occur within the nucleoprotein complex of chromatin which is dynamic in nature, capable of rapid unfolding, disassembling, assembling and refolding. The present work reviews experimental work designed to investigate the impact of chromatin dynamics and chromosome conformation changes during cell-cycle in the formation of chromosomal aberrations. Using conventional cytogenetics and premature chromosome condensation to visualize interphase chromatin, the data presented support the hypothesis that chromatin dynamic changes during cell-cycle are important determinants in the conversion of sub-microscopic DNA lesions into chromatid breaks. Consequently, the type and yield of radiation-induced chromosomal aberrations at a given cell-cycle-stage depends on the combined effect of DNA repair processes and chromatin dynamics, which is cell-cycle-regulated and subject to up- or down-regulation following radiation exposure or genetic alterations. This new hypothesis is used to explain the variability in radiosensitivity observed at various cell-cycle-stages, among mutant cells and cells of different origin, or among different individuals, and to revisit unresolved issues and unanswered questions. In addition, it is used to better understand hypersensitivity of AT cells and to provide an improved predictive G2-assay for evaluating radiosensitivity at individual level. Finally, experimental data at single cell level obtained using hybrid cells suggest that the proposed hypothesis applies only to the irradiated component of the hybrid.