Irregular variations in radiation sensitivity when the linear energy transfer is increased

Ionization Detail Parameters for DNA Damage Evaluation in Charged Particle Radiotherapy: Simulation Study Based on Cell Survival Database

Details of excitation and ionization acts hide a description of the biological effects of charged particle traversal through living tissue. Nanodosimetry enables the introduction of novel quantities that characterize and quantify the particle track structure while also serving as a foundation for assessing biological effects based on this quantification. This presents an opportunity to enhance the planning of charged particle radiotherapy by taking into account the ionization detail. This work uses Monte Carlo simulations with Geant4-DNA code for a wide variety of charged particles and their radiation qualities to analyze the distribution of ionization cluster sizes within nanometer-scale volumes, similar to DNA diameter. By correlating these results with biological parameters extracted from the PIDE database for the V79 cell line, a novel parameter R2 based on ionization details is proposed for the evaluation of radiation quality in terms of biological consequences, i.e., radiobiological cross section for inactivation. By incorporating the probability p of sub-lethal damage caused by a single ionization, we address limitations associated with the usually proposed nanodosimetric parameter Fk for characterizing the biological effects of radiation. We show that the new parameter R2 correlates well with radiobiological data and can be used to predict biological outcomes.

Ra-223 induces clustered DNA damage and inhibits cell survival in several prostate cancer cell lines

The bone-seeking radiopharmaceutical Xofigo (Radium-223 dichloride) has demonstrated both extended survival and palliative effects in treatment of bone metastases in prostate cancer. The alpha-particle emitter Ra-223, targets regions undergoing active bone remodeling and strongly binds to bone hydroxyapatite (HAp). However, the toxicity mechanism and properties of Ra-223 binding to hydroxyapatite are not fully understood. By exposing 2D and 3D (spheroid) prostate cancer cell models to free and HAp-bound Ra-223 we here studied cell toxicity, apoptosis and formation and repair of DNA double-strand breaks (DSBs). The rapid binding with a high affinity of Ra-223 to bone-like HAp structures was evident (KD= 19.2 × 10-18 M) and almost no dissociation was detected within 24 h. Importantly, there was no significant uptake of Ra-223 in cells. The Ra-223 alpha-particle decay produced track-like distributions of the DNA damage response proteins 53BP1 and ɣH2AX induced high amounts of clustered DSBs in prostate cancer cells and activated DSB repair through non-homologous end-joining (NHEJ). Ra-223 inhibited growth of prostate cancer cells, independent of cell type, and induced high levels of apoptosis. In summary, we suggest the high cell killing efficacy of the Ra-223 was attributed to the clustered DNA damaged sites induced by α-particles.

Development of a new microdosimetric biological weighting function for the RBE10 assessment in case of the V79 cell line exposed to ions from 1H to 238U

An improved biological weighting function (IBWF) is proposed to phenomenologically relate microdosimetric lineal energy probability density distributions with the relative biological effectiveness (RBE) for the in vitro clonogenic cell survival (surviving fraction = 10%) of the most commonly used mammalian cell line, i.e. the Chinese hamster lung fibroblasts (V79). The IBWF, intended as a simple and robust tool for a fast RBE assessment to compare different exposure conditions in particle therapy beams, was determined through an iterative global-fitting process aimed to minimize the average relative deviation between RBE calculations and literature in vitro data in case of exposure to various types of ions from ¹H to ²³⁸U. By using a single particle- and energy- independent function, it was possible to establish an univocal correlation between lineal energy and clonogenic cell survival for particles spanning over an unrestricted linear energy transfer range of almost five orders of magnitude (0.2 keV µm^-1 to 15 000 keV µm^-1 in liquid water). The average deviation between IBWF-derived RBE values and the published in vitro data was ∼14%. The IBWF results were also compared with corresponding calculations (in vitro RBE₁₀ for the V79 cell line) performed using the modified microdosimetric kinetic model (modified MKM). Furthermore, RBE values computed with the reference biological weighting function (BWF) for the in vivo early intestine tolerance in mice were included for comparison and to further explore potential correlations between the BWF results and the in vitro RBE as reported in previous studies. The results suggest that the modified MKM possess limitations in reproducing the experimental in vitro RBE₁₀ for the V79 cell line in case of ions heavier than ²⁰Ne. Furthermore, due to the different modelled endpoint, marked deviations were found between the RBE values assessed using the reference BWF and the IBWF for ions heavier than ²H. Finally, the IBWF was unchangingly applied to calculate RBE values by processing lineal energy density distributions experimentally measured with eight different microdosimeters in 19 ¹H and ¹²C beams at ten different facilities (eight clinical and two research ones). Despite the differences between the detectors, irradiation facilities, beam profiles (pristine or spread out Bragg peak), maximum beam energy, beam delivery (passive or active scanning), energy degradation system (water, PMMA, polyamide or low-density polyethylene), the obtained IBWF-based RBE trends were found to be in good agreement with the corresponding ones in case of computer-simulated microdosimetric spectra (average relative deviation equal to 0.8% and 5.7% for ¹H and ¹²C ions respectively).

Charged particle radiobiology beamline using tandem accelerator-based MeV protons and carbon ions: a pilot study on the track-end radiation quality, variable biological effectiveness and Bayesian beam dosimetry

For in vitro cell irradiation using tandem accelerator-based MeV protons and carbon ions, by TOPAS simulation, a pilot study of performance evaluation is presented on a collimation beamline for 3 MeV protons and 10 MeV carbon ions from a 2  ×  3 MV tandem accelerator. Based on the elements and source parameters, a collimated beam of 2.8 MeV protons or 2.5 MeV carbon ions, with 5.175 mm or 5.166 mm full width tenth maximum (FWTM), respectively, can be delivered to the target cell dish. TOPAS simulations and/or deterministic algorithms present a Bragg curve of linear energy transfer (LET) (10-70 keV μm^-1) along a 138 μm range of the proton beam, and a declining LET of the carbon beam (900-100 keV μm^-1) within 4 μm range. Based on the biophysical models for relative biological effectiveness (RBE) of protons, TOPAS RBE scorers presents a set of depth-variation curves of the proton RBE (for V79 and DU145 cells), linearly related to the Bragg curve of the proton LET. Based on the microdosimetric-kinetic (MK) theory, in the 4 μm range for a monolayer cell thickness, the mean RBEα (V79 cells) of the carbon ion beam is estimated as 3.612 (late S phase) and 1.737 (G ₁/S phase) for the mean LET of 492 keV μm^-1. For practical irradiations, a tunable proton RBE can be acquired by changing the thickness of the cell dish. For the low-energy high-fluence (rate) beams, indirect beam measurements are proposed to detect the proton-beam induced scattering/recoil protons from a beam-intercepting Mylar film, and the carbon-beam induced backscattered electrons from a gold-deposited Havar-foil beam window. Statistical dosimetry for the indirect measurement is established, using a Bayesian model based on the preset number of detection counts, by which the mean value of the whole-dish dose can be prescribed and the uncertainty introduced in the survival data can be corrected.

Modelling of Cellular Survival Following Radiation-Induced DNA Double-Strand Breaks

A mechanistic model of cellular survival following radiation-induced DNA double-strand breaks (DSBs) was proposed in this study. DSBs were assumed as the initial lesions in the DNA of the cell nucleus induced by ionizing radiation. The non-homologous end-joining (NHEJ) pathway was considered as the domain pathway of DSB repair in mammalian cells. The model was proposed to predict the relationship between radiation-induced DSBs in nucleus and probability of cell survival, which was quantitatively described by two input parameters and six fitting parameters. One input parameter was the average number of primary particles which caused DSB, the other input parameter was the average number of DSBs yielded by each primary particle that caused DSB. The fitting parameters were used to describe the biological characteristics of the irradiated cells. By determining the fitting parameters of the model with experimental data, the model is able to estimate surviving fractions for the same type of cells exposed to particles with different physical parameters. The model further revealed the mechanism of cell death induced by the DSB effect. Relative biological effectiveness (RBE) of charged particles at different survival could be calculated with the model, which would provide reference for clinical treatment.

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.

Effects of Charged Particles on Human Tumor Cells

The use of charged particle therapy in cancer treatment is growing rapidly, in large part because the exquisite dose localization of charged particles allows for higher radiation doses to be given to tumor tissue while normal tissues are exposed to lower doses and decreased volumes of normal tissues are irradiated. In addition, charged particles heavier than protons have substantial potential clinical advantages because of their additional biological effects, including greater cell killing effectiveness, decreased radiation resistance of hypoxic cells in tumors, and reduced cell cycle dependence of radiation response. These biological advantages depend on many factors, such as endpoint, cell or tissue type, dose, dose rate or fractionation, charged particle type and energy, and oxygen concentration. This review summarizes the unique biological advantages of charged particle therapy and highlights recent research and areas of particular research needs, such as quantification of relative biological effectiveness (RBE) for various tumor types and radiation qualities, role of genetic background of tumor cells in determining response to charged particles, sensitivity of cancer stem-like cells to charged particles, role of charged particles in tumors with hypoxic fractions, and importance of fractionation, including use of hypofractionation, with charged particles.

Data-driven RBE parameterization for helium ion beams

Helium ion beams are expected to be available again in the near future for clinical use. A suitable formalism to obtain relative biological effectiveness (RBE) values for treatment planning (TP) studies is needed. In this work we developed a data-driven RBE parameterization based on published in vitro experimental values. The RBE parameterization has been developed within the framework of the linear-quadratic (LQ) model as a function of the helium linear energy transfer (LET), dose and the tissue specific parameter (α/β)ph of the LQ model for the reference radiation. Analytic expressions are provided, derived from the collected database, describing the RBEα = αHe/αph and Rβ = βHe/βph ratios as a function of LET. Calculated RBE values at 2 Gy photon dose and at 10% survival (RBE10) are compared with the experimental ones. Pearson's correlation coefficients were, respectively, 0.85 and 0.84 confirming the soundness of the introduced approach. Moreover, due to the lack of experimental data at low LET, clonogenic experiments have been performed irradiating A549 cell line with (α/β)ph = 5.4 Gy at the entrance of a 56.4 MeV u(-1)He beam at the Heidelberg Ion Beam Therapy Center. The proposed parameterization reproduces the measured cell survival within the experimental uncertainties. A RBE formula, which depends only on dose, LET and (α/β)ph as input parameters is proposed, allowing a straightforward implementation in a TP system.

Track-event theory of cell survival with second-order repair

When fractionation schemes for hypofractionation and stereotactic body radiotherapy are considered, a reliable cell survival model at high dose is needed for calculating doses of similar biological effectiveness. In this work, a simple model for cell survival which is valid also at high dose is developed from Poisson statistics. It is assumed that a cell is killed by an event that is defined by two double-strand breaks on the same or different chromosomes. Two different mechanisms can produce events. A one-track event is always represented by two simultaneous double-strand breaks. A two-track event results in one double-strand break. Therefore, at least two two-track events on the same or different chromosomes are necessary to produce an event. It is assumed that two double-strand breaks can be repaired with a certain repair probability. Both the one-track events and the two-track events are statistically independent. From the stochastic nature of cell killing which is described by the Poisson distribution, the cell survival probability was derived. The model was fitted to experimental data. It was shown that a solution based on Poisson statistics exists for cell survival. It exhibits exponential cell survival at high dose and a finite gradient of cell survival at vanishing dose, which is in agreement with experimental cell studies. The model fits the experimental data as well as the LQ model and is based on two free parameters. It was shown that cell survival can be described with a simple analytical formula on the basis of Poisson statistics. This solution represents in the limit of large dose the typical exponential behavior and predicts cell survival as well as the LQ model.

Systematic analysis of RBE and related quantities using a database of cell survival experiments with ion beam irradiation

For tumor therapy with light ions and for experimental aspects in particle radiobiology the relative biological effectiveness (RBE) is an important quantity to describe the increased effectiveness of particle radiation. By establishing and analysing a database of ion and photon cell survival data, some remarkable properties of RBE-related quantities were observed. The database consists of 855 in vitro cell survival experiments after ion and photon irradiation. The experiments comprise curves obtained in different labs, using different ion species, different irradiation modalities, the whole range of accessible energies and linear energy transfers (LETs) and various cell types. Each survival curve has been parameterized using the linear-quadratic (LQ) model. The photon parameters, α and β, appear to be slightly anti-correlated, which might point toward an underlying biological mechanism. The RBE values derived from the survival curves support the known dependence of RBE on LET, on particle species and dose. A positive correlation of RBE with the ratio α/β of the photon LQ parameters is found at low doses, which unexpectedly changes to a negative correlation at high doses. Furthermore, we investigated the course of the β coefficient of the LQ model with increasing LET, finding typically a slight initial increase and a final falloff to zero. The observed fluctuations in RBE values of comparable experiments resemble overall RBE uncertainties, which is of relevance for treatment planning. The database can also be used for extensive testing of RBE models. We thus compare simulations with the local effect model to achieve this goal.

In vitro RBE-LET dependence for multiple particle types

BACKGROUND

In vitro RBE values for various high LET radiation types have been determined for many different cell types. Occasionally it is criticized that RBE for a given endpoint cannot be single-value dependent on LET alone, but also on particle species, due to the different dose deposition profiles on microscopic scale. Hence LET is not sufficient as a predictor of RBE, and this is one of the motivations for development of radiobiological models which explicitly depend on the detailed particle energy spectrum of the applied radiation field. The aim of the present study is to summarize the available data in the literature regarding the dependency of RBE on LET for different particles.

METHOD

As RBE is highly dependent on cell type and endpoint, we discriminated the RBE-LET relationship for the three investigated cell lines and at the same endpoint (10% survival in colony formation). Data points were collected from 20, four and four publications for V79, CHO and T1, respectively, in total covering 228 RBE values from a broad range of particle species.

RESULTS AND DISCUSSION

All RBE-LET data points demonstrate surprising agreement within the general error band formed by the numerous data points, and display the expected RBE peak at around 100-200 keV/μm. For all three cell lines, the influence of varying the particle type on the RBE was far from obvious, compared to the general experimental noise. Therefore, a dependence of particle type cannot be concluded, and LET alone in fact does seem to be an adequate parameter for describing RBE at 10% survival.

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

PURPOSE

The formulation of relative biological effectiveness (RBE) for high linear energy transfer (high-LET) radiation treatments is revisited. The effects of changed production of sub-lethal damage with varying LET is now considered via the RBEmin concept, where RBEmin represents the lower limit to which RBE tends at high doses per fraction.

MATERIALS AND METHODS

An existing linear-quadratic formulation for calculating RBE variations with fractional dose for high-LET radiations is modified to incorporate the twin concepts of RBEmax (which represents the value of RBE at an effective dose-per-fraction of 0 Gy) and RBEmin.

RESULTS

Fits of the model to data showed RBEmin values in the range of 0.1- 2.27. In all cases the raw data was a better statistical fit to the model which included RBEmin, although this was only very highly significant in one case. In the case of the mouse oesophagus it is shown that, if change in the beta-radiosensitivity coefficient with LET is considered as trivial, an underestimation > 5% in RBE can be expected at X-ray doses of 2 Gy/fraction if RBEmin is not considered. To ensure that the results were not biased by the statistical method used to obtain the parameter values relevant to this analysis (i.e., using fraction-size effect or Fe-plots), an alternative method was used which provided very similar correlation with the data.

CONCLUSIONS

If the production of sublethal damage is considered independent of LET, there will be a risk that non-corrected evaluation of RBE will lead to an over- or under-estimate of RBE at low doses per fractions (the clinically relevant region).

Differences in radiosensitivity between three HER2 overexpressing cell lines

PURPOSE

HER2 is a potential target for radionuclide therapy, especially when HER2 overexpressing breast cancer cells are resistant to Herceptin(R) treatment. Therefore, it is of interest to analyse whether HER2 overexpressing tumour cells have different inherent radiosensitivity.

METHODS

The radiosensitivity of three often used HER2 overexpressing cell lines, SKOV-3, SKBR-3 and BT-474, was analysed. The cells were exposed to conventional photon irradiation, low linear energy transfer (LET), to characterise their inherent radiosensitivity. The analysis was made with clonogenic survival and growth extrapolation assays. The cells were also exposed to alpha particles, high LET, from (211)At decays using the HER2-binding affibody molecule (211)At-(Z(HER2:4))(2) as targeting agent. Assays for studies of internalisation of the affibody molecule were applied.

RESULTS

SKOV-3 cells were most radioresistant, SKBR-3 cells were intermediate and BT-474 cells were most sensitive as measured with the clonogenic and growth extrapolation assays after photon irradiation. The HER2 dependent cellular uptake of (211)At was qualitatively similar for all three cell lines. However, the sensitivity to the alpha particles from (211)At differed; SKOV-3 was most resistant, SKBR-3 intermediate and BT-474 most sensitive. These differences were unexpected because it is assumed that all types of cells should have similar sensitivity to high-LET radiation. The sensitivity to alpha particle exposure correlated with internalisation of the affibody molecule and with size of the cell nucleus.

CONCLUSION

There can be differences in radiosensitivity, which, if they also exist between patient breast cancer cells, are important to consider for both conventional radiotherapy and for HER2-targeted radionuclide therapy.

Cellular parameters and RBE-LET dependences for modelling heavy-ion radiotherapy

Sets of four parameters (m, E0, sigma0 and kappa) of the cellular track structure model of Katz have been fitted to recently published data concerning human melanoma (AA) and mammalian (V79) cells exposed to a variety of lighter ions and to mixed ion-Co60 and ion-ion irradiation. Using these parameters, model predictions of V79 survival were verified against experimental data. RBE-LET dependences were calculated and compared with experimental data obtained for V79 cells after exposure to 3He, 12C and 20Ne ion beams. The presented track-segment approach used in track structure calculations, while satisfactory for heavier ions, may be of limited value for predicting the RBE-LET dependence of proton and helium radiotherapy beams in regions close to the distal range of these particles. We discuss the predictive capability of this model and propose standards in reporting cellular radiobiology data for application in modelling heavy ion beam radiotherapy.

The potential of proton and light ion beams in radiotherapy

In an investigation by the Swedish Cancer Society, the present status, critical issues and future aspects and potentials were described by an expert group for each of nine major areas of radiation therapy research. The present report deals with the potential of proton and light ion beams in radiotherapy.

Cellular retention of radioactivity and increased radiation dose. Model experiments with EGF-dextran

Targeting of tumor cells with radiolabeled biomolecules is a possible approach to inactivate disseminated tumor cells. However, rapid degradation of the biomolecules after cellular internalization and subsequent excretion of the radioactivity is a problem. We studied the possibility of using dextran as a carrier of radionuclides to improve the intracellular retention. An EGF-dextran conjugate, aimed for targeting of tumor cells overexpressing the EGF-receptor, was used as model. Retention tests were performed with (125)I on different parts: [(125)I]-EGF-dextran-[(125)I], [(125)I]-EGF-dextran and EGF-dextran-[(125)I]. Comparisons were made with [(125)I]-EGF. The radiolabeled compounds were incubated with cultured glioma cells for different times. The cellular retention of radioactivity was then measured for up to 24 h. Expected radiation doses at the cellular level were calculated assuming that (131)I, instead of (125)I, was coupled to EGF and EGF-dextran. The results indicated that the EGF-part of the conjugate was degraded and the EGF-attached radioactivity was rapidly excreted, whereas radioactivity on dextran was retained intracellularly to a high degree, i.e. 70-80% of the radioactivity bound to dextran was still cell-associated after 24 h. The retention after 24 h was significantly higher (p < 0.001) when the radioactivity was on the dextran instead of the EGF-part. The radiolabeled EGF-dextran had a notably high specific radioactivity; up to 11 MBq/microg. There was potential for at least hundred times increased radiation dose per receptor interaction when the radioactivity was on the dextran part. The advantage with radioactivity on the dextran part was the high cellular retention and the high specific radioactivity (higher than previously reported for other residualizing labels) without severe loss of receptor specific binding. Thus, dextran seems suitable as a carrier of radionuclides aimed for therapy and gives potential for a highly increased radiation dose.

Cell death induced by high-linear-energy transfer carbon beams in human glioblastoma cell lines

The cytotoxic effect of high-linear-energy transfer (LET) carbon beams on two human glioblastoma cell lines (A172 and TK1) was analyzed, especially concerning cell death, including apoptosis. Gamma-ray radiation was used for comparison. The results of standard colony formation assay showed that the survival fraction of each cell line decreased in an LET-dependent manner. The results of other direct cytotoxic assays, dye exclusion test, and lactate dehydrogenase (LDH) release assay, also displayed a similar relationship between the cytotoxic effect of carbon beams and LET. The maximum values of the cell death index (CDI) were 50.2% in A172 and 37.5% in TK1, both obtained on day 7 after exposure to carbon beams of 80 keV/microm. Apoptosis was observed only on days 4 and 7 after carbon beam irradiation, with maximum values of 7% in A172 and 4.5% in TK1, and the induction of apoptosis after high-LET radiation could be p53-independent. This indicated that a combination of multiple assays to detect cell death was important in evaluating the radiosensitivity of tumor cells, because this approach could more precisely reflect the clinical effectiveness of radiotherapy.

The assessment of RBE effects using the concept of biologically effective dose

PURPOSE

To modify existing linear-quadratic (LQ) equations in order to take account of relative biological effectiveness (RBE) using the concept of biologically effective dose (BED).

METHODS AND MATERIALS

Clinically useful forms of the LQ model have been modified to incorporate RBE effects in such a way as to allow comparison between high- and low-LET (linear energy transfer) radiations in terms of similar biological dose units. The new parameter in the formulation is RBEM, the intrinsic (or maximum) RBE at zero dose. The principal assumption (following Kellerer and Rossi; ref. 1) is that high-LET radiation modifies the alpha-coefficient of damage while leaving the beta-coefficient unaltered.

RESULTS

The equations allow a quantitative estimation of how the apparent RBE will change with changes in dose/fraction or dose-rate and of how the magnitude and rate of change is governed by the low-LET alpha/beta ratio of the irradiated tissue. The modifications are applicable to all types of radiotherapy (fractionated, continuous low dose-rate, therapy with decaying sources, etc.). In cases where the normal tissue RBEM is greater than that for the tumor, the revised formulation helps explain why there will be situations where therapeutic index will be adversely affected by use of high-LET radiation. Such clinical advantages as have been observed are more likely to result from favorable geometrical sparing of critical normal tissues and/or the fact that slowly growing tumors may have alpha/beta values more typical of late-responding normal tissues.

CONCLUSIONS

The incorporation of RBE into existing LQ methodology allows quantitative assessment of clinical applications of high-LET radiations via an examination of the associated BEDs. On the basis of such assessments high-LET radiations are shown to confer few advantages.