Response of Chinese Hamster V79 Multicellular Spheroids Exposed to High-Energy Carbon Ions

Irradiation with Carbon Ions Effectively Counteracts Hypoxia-related Radioresistance in a Rat Prostate Carcinoma

PURPOSE

Hypoxia in tumors is associated with increased malignancy and resistance to conventional photon radiation therapy. This study investigated the potential of particle therapy to counteract radioresistance in syngeneic rat prostate carcinoma.

METHODS AND MATERIALS

Subcutaneously transplanted R3327-HI tumors were irradiated with photons or carbon ions under acute hypoxic conditions, induced by clamping the tumor-supplying artery 10 min before and during irradiation. Dose-response curves were determined for the endpoint "local tumor control within 300 days" and compared with previously published data acquired under oxic conditions. Doses at 50% tumor control probability (TCD50) were used to quantify hypoxia-induced radioresistance relative to that under oxic conditions and also to quantify the increased effectiveness of carbon ions under oxic and hypoxic conditions relative to photons.

RESULTS

Compared with those under oxic conditions, TCD50 values under hypoxic conditions increased by a factor of 1.53 ± 0.08 for photons and by a factor of 1.28 ± 0.06 for carbon ions (oxygen enhancement ratio). Compared with those for photons, TCD50 values for carbon ions decreased by a factor of 2.08 ± 0.13 under oxic conditions and by a factor of 2.49 ± 0.08 under hypoxic conditions (relative biological effectiveness). While the slope of the photon dose-response curves increased when changing from oxic to hypoxic conditions, the slopes were steeper and remained unchanged for carbon ions.

CONCLUSIONS

The reduced oxygen enhancement ratio of carbon ions indicated that the required dose increase in hypoxic tumors was 17% lower for carbon ions than for photons. Additionally, carbon ions reduced the effect of intertumor heterogeneity on the radiation response. Therefore, carbon ions may confer a significant advantage for the treatment of hypoxic tumors that are highly resistant to conventional photon radiation therapy.

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.

Proton Irradiation Increases the Necessity for Homologous Recombination Repair Along with the Indispensability of Non-Homologous End Joining

Technical improvements in clinical radiotherapy for maximizing cytotoxicity to the tumor while limiting negative impact on co-irradiated healthy tissues include the increasing use of particle therapy (e.g., proton therapy) worldwide. Yet potential differences in the biology of DNA damage induction and repair between irradiation with X-ray photons and protons remain elusive. We compared the differences in DNA double strand break (DSB) repair and survival of cells compromised in non-homologous end joining (NHEJ), homologous recombination repair (HRR) or both, after irradiation with an equal dose of X-ray photons, entrance plateau (EP) protons, and mid spread-out Bragg peak (SOBP) protons. We used super-resolution microscopy to investigate potential differences in spatial distribution of DNA damage foci upon irradiation. While DNA damage foci were equally distributed throughout the nucleus after X-ray photon irradiation, we observed more clustered DNA damage foci upon proton irradiation. Furthermore, deficiency in essential NHEJ proteins delayed DNA repair kinetics and sensitized cells to both, X-ray photon and proton irradiation, whereas deficiency in HRR proteins sensitized cells only to proton irradiation. We assume that NHEJ is indispensable for processing DNA DSB independent of the irradiation source, whereas the importance of HRR rises with increasing energy of applied irradiation.

Full Monte Carlo-Based Biologic Treatment Plan Optimization System for Intensity Modulated Carbon Ion Therapy on Graphics Processing Unit

PURPOSE

One of the major benefits of carbon ion therapy is enhanced biological effectiveness at the Bragg peak region. For intensity modulated carbon ion therapy (IMCT), it is desirable to use Monte Carlo (MC) methods to compute the properties of each pencil beam spot for treatment planning, because of their accuracy in modeling physics processes and estimating biological effects. We previously developed goCMC, a graphics processing unit (GPU)-oriented MC engine for carbon ion therapy. The purpose of the present study was to build a biological treatment plan optimization system using goCMC.

METHODS AND MATERIALS

The repair-misrepair-fixation model was implemented to compute the spatial distribution of linear-quadratic model parameters for each spot. A treatment plan optimization module was developed to minimize the difference between the prescribed and actual biological effect. We used a gradient-based algorithm to solve the optimization problem. The system was embedded in the Varian Eclipse treatment planning system under a client-server architecture to achieve a user-friendly planning environment. We tested the system with a 1-dimensional homogeneous water case and 3 3-dimensional patient cases.

RESULTS

Our system generated treatment plans with biological spread-out Bragg peaks covering the targeted regions and sparing critical structures. Using 4 NVidia GTX 1080 GPUs, the total computation time, including spot simulation, optimization, and final dose calculation, was 0.6 hour for the prostate case (8282 spots), 0.2 hour for the pancreas case (3795 spots), and 0.3 hour for the brain case (6724 spots). The computation time was dominated by MC spot simulation.

CONCLUSIONS

We built a biological treatment plan optimization system for IMCT that performs simulations using a fast MC engine, goCMC. To the best of our knowledge, this is the first time that full MC-based IMCT inverse planning has been achieved in a clinically viable time frame.

Initial development of goCMC: a GPU-oriented fast cross-platform Monte Carlo engine for carbon ion therapy

Monte Carlo (MC) simulation is considered as the most accurate method for calculation of absorbed dose and fundamental physics quantities related to biological effects in carbon ion therapy. To improve its computational efficiency, we have developed a GPU-oriented fast MC package named goCMC, for carbon therapy. goCMC simulates particle transport in voxelized geometry with kinetic energy up to 450 MeV u-1. Class II condensed history simulation scheme with a continuous slowing down approximation was employed. Energy straggling and multiple scattering were modeled. δ-electrons were terminated with their energy locally deposited. Four types of nuclear interactions were implemented in goCMC, i.e. carbon-hydrogen, carbon-carbon, carbon-oxygen and carbon-calcium inelastic collisions. Total cross section data from Geant4 were used. Secondary particles produced in these interactions were sampled according to particle yield with energy and directional distribution data derived from Geant4 simulation results. Secondary charged particles were transported following the condensed history scheme, whereas secondary neutral particles were ignored. goCMC was developed under OpenCL framework and is executable on different platforms, e.g. GPU and multi-core CPU. We have validated goCMC with Geant4 in cases with different beam energy and phantoms including four homogeneous phantoms, one heterogeneous half-slab phantom, and one patient case. For each case [Formula: see text] carbon ions were simulated, such that in the region with dose greater than 10% of maximum dose, the mean relative statistical uncertainty was less than 1%. Good agreements for dose distributions and range estimations between goCMC and Geant4 were observed. 3D gamma passing rates with 1%/1 mm criterion were over 90% within 10% isodose line except in two extreme cases, and those with 2%/1 mm criterion were all over 96%. Efficiency and code portability were tested with different GPUs and CPUs. Depending on the beam energy and voxel size, the computation time to simulate [Formula: see text] carbons was 9.9-125 s, 2.5-50 s and 60-612 s on an AMD Radeon GPU card, an NVidia GeForce GTX 1080 GPU card and an Intel Xeon E5-2640 CPU, respectively. The combined accuracy, efficiency and portability make goCMC attractive for research and clinical applications in carbon ion therapy.

Multiscale approach predictions for biological outcomes in ion-beam cancer therapy

Ion-beam therapy provides advances in cancer treatment, offering the possibility of excellent dose localization and thus maximising cell-killing within the tumour. The full potential of such therapy can only be realised if the fundamental mechanisms leading to lethal cell damage under ion irradiation are well understood. The key question is whether it is possible to quantitatively predict macroscopic biological effects caused by ion radiation on the basis of physical and chemical effects related to the ion-medium interactions on a nanometre scale. We demonstrate that the phenomenon-based MultiScale Approach to the assessment of radiation damage with ions gives a positive answer to this question. We apply this approach to numerous experiments where survival curves were obtained for different cell lines and conditions. Contrary to other, in essence empirical methods for evaluation of macroscopic effects of ionising radiation, the MultiScale Approach predicts the biodamage based on the physical effects related to ionisation of the medium, transport of secondary particles, chemical interactions, thermo-mechanical pathways of biodamage, and heuristic biological criteria for cell survival. We anticipate this method to give great impetus to the practical improvement of ion-beam cancer therapy and the development of more efficient treatment protocols.

Carbon ion radiotherapy decreases the impact of tumor heterogeneity on radiation response in experimental prostate tumors

OBJECTIVE

To quantitatively study the impact of intrinsic tumor characteristics and microenvironmental factors on local tumor control after irradiation with carbon ((12)C-) ions and photons in an experimental prostate tumor model.

MATERIAL AND METHODS

Three sublines of a syngeneic rat prostate tumor (R3327) differing in grading (highly (-H) moderately (-HI) or anaplastic (-AT1)) were irradiated with increasing single doses of either (12)C-ions or 6 MV photons in Copenhagen rats. Primary endpoint was local tumor control within 300 days. The relative biological effectiveness (RBE) of (12)C-ions was calculated from the dose at 50% tumor control probability (TCD50) of photons and (12)C-ions and was correlated with histological, physiological and genetic tumor parameters.

RESULTS

Experimental findings demonstrated that (i) TCD50-values between the three tumor sublines differed less for (12)C-ions (23.6-32.9 Gy) than for photons (38.2-75.7 Gy), (ii) the slope of the dose-response curve for each tumor line was steeper for (12)C-ions than for photons, and (iii) the RBE increased with tumor grading from 1.62 ± 0.11 (H) to 2.08 ± 0.13 (HI) to 2.30 ± 0.08 (AT1).

CONCLUSION

The response to (12)C-ions is less dependent on resistance factors as well as on heterogeneity between and within tumor sublines as compared to photons. A clear correlation between decreasing differentiation status and increasing RBE was found. (12)C-ions may therefore be a therapeutic option especially in patients with undifferentiated prostate tumors, expressing high resistance against photons.

Differential Superiority of Heavy Charged-Particle Irradiation to X-Rays: Studies on Biological Effectiveness and Side Effect Mechanisms in Multicellular Tumor and Normal Tissue Models

This review is focused on the radiobiology of carbon ions compared to X-rays using multicellular models of tumors and normal mucosa. The first part summarizes basic radiobiological effects, as observed in cancer cells. The second, more clinically oriented part of the review, deals with radiation-induced cell migration and mucositis. Multicellular spheroids from V79 hamster cells were irradiated with X-rays or carbon ions under ambient or restricted oxygen supply conditions. Reliable oxygen enhancement ratios could be derived to be 2.9, 2.8, and 1.4 for irradiation with photons, ¹²C⁺⁶ in the plateau region, and ¹²C⁺⁶ in the Bragg peak, respectively. Similarly, a relative biological effectiveness of 4.3 and 2.1 for ambient pO₂ and hypoxia was obtained, respectively. The high effectiveness of carbon ions was reflected by an enhanced accumulation of cells in G₂/M and a dose-dependent massive induction of apoptosis. These data clearly show that heavy charged particles are more efficient in sterilizing tumor cells than conventional irradiation even under hypoxic conditions. Clinically relevant doses (3 Gy) of X-rays induced an increase in migratory activity of U87 but not of LN229 or HCT116 tumor cells. Such an increase in cell motility following irradiation in situ could be the source of recurrence. In contrast, carbon ion treatment was associated with a dose-dependent decrease in migration with all cell lines and under all conditions investigated. The radiation-induced loss of cell motility was correlated, in most cases, with corresponding changes in β₁ integrin expression. The photon-induced increase in cell migration was paralleled by an elevated phosphorylation status of the epidermal growth factor receptor and AKT-ERK1/2 pathway. Such a hyperphosphorylation did not occur during ¹²C⁺⁶ irradiation under all conditions registered. Comparing the gene toxicity of X-rays with that of particles using the γH2AX technique in organotypic cultures of the oral mucosa, the superior effectiveness of heavy ions was confirmed by a twofold higher number of foci per nucleus. However, proinflammatory signs were similar for both treatment modalities, e.g., the activation of NFκB and the release of IL6 and IL8. The presence of peripheral blood mononuclear cell increased the radiation-induced release of the proinflammatory cytokines by factors of 2–3. Carbon ions are part of the cosmic radiation. Long-term exposure to such particles during extended space flights, as planned by international space agencies, may thus impose a medical and safety risk on the astronauts by a potential induction of mucositis. In summary, particle irradiation is superior to gamma-rays due to a higher radiobiological effectiveness, a reduced hypoxia-induced radioresistance, a multicellular radiosensitization, and the absence of a radiation-induced cell motility. However, the potential of inducing mucositis is similar for both radiation types.

Spatio-Temporal Dynamics of Hypoxia during Radiotherapy

Tumour hypoxia plays a pivotal role in cancer therapy for most therapeutic approaches from radiotherapy to immunotherapy. The detailed and accurate knowledge of the oxygen distribution in a tumour is necessary in order to determine the right treatment strategy. Still, due to the limited spatial and temporal resolution of imaging methods as well as lacking fundamental understanding of internal oxygenation dynamics in tumours, the precise oxygen distribution map is rarely available for treatment planing. We employ an agent-based in silico tumour spheroid model in order to study the complex, localized and fast oxygen dynamics in tumour micro-regions which are induced by radiotherapy. A lattice-free, 3D, agent-based approach for cell representation is coupled with a high-resolution diffusion solver that includes a tissue density-dependent diffusion coefficient. This allows us to assess the space- and time-resolved reoxygenation response of a small subvolume of tumour tissue in response to radiotherapy. In response to irradiation the tumour nodule exhibits characteristic reoxygenation and re-depletion dynamics which we resolve with high spatio-temporal resolution. The reoxygenation follows specific timings, which should be respected in treatment in order to maximise the use of the oxygen enhancement effects. Oxygen dynamics within the tumour create windows of opportunity for the use of adjuvant chemotherapeutica and hypoxia-activated drugs. Overall, we show that by using modelling it is possible to follow the oxygenation dynamics beyond common resolution limits and predict beneficial strategies for therapy and in vitro verification. Models of cell cycle and oxygen dynamics in tumours should in the future be combined with imaging techniques, to allow for a systematic experimental study of possible improved schedules and to ultimately extend the reach of oxygenation monitoring available in clinical treatment.

Determination of the relative biological effectiveness and oxygen enhancement ratio for micronuclei formation using high-LET radiation in solid tumor cells: An in vitro and in vivo study

We determined the relative biological effectiveness (RBE) and oxygen enhancement ratio (OER) of micronuclei (MN) formation in clamped (hypoxic) and non-clamped (normoxic) solid tumors in mice legs following exposure to X-rays and heavy ions. Single-cell suspensions (aerobic) of non-irradiated tumors were prepared in parallel and used directly to determine the radiation response for aerobic cells. Squamous cell carcinoma (SCCVII) cells were transplanted into the right hind legs of syngeneic C3H/He male mice. Irradiation doses with either X-rays or heavy ions at a dose-averaged LET (linear energy transfer) of 14-192keV/μm were delivered to 5-mm diameter tumors and aerobic single-cells in sample-tubes. After irradiation, the tumors were excised and trypsinized to observe MN in single-cells. The single-cell suspensions were used for MN formation assays. The RBE values increased with increasing LET. The maximum RBE values for the three different oxygen conditions; hypoxic tumor, normoxic tumor, and aerobic cells, were 8.18, 5.30, and 3.76 at an LET of 192keV/μm, respectively. After X-irradiation, the OERh/n values (hypoxic tumor/normoxic tumor) were lower than the OERh/a (hypoxic tumor/aerobic cells), and were 1.73 and 2.58, respectively. We found that the OER for the in vivo studies were smaller in comparison to that for the in vitro studies. Both of the OER values at 192keV/μm were small in comparison to those of the X-ray irradiated samples. The OERh/n and OERh/a values at 192keV/μm were 1.12 and 1.19, respectively. Our results suggest that high LET radiation has a large biological effect even if a solid tumor includes substantial numbers of hypoxic cells. To conclude, we found that the RBE values under each oxygen state for non-MN fraction increased with increasing LET and that the OER values for both tumors in vivo and cells in vitro decreased with increasing LET.

Radiobiological description of the LET dependence of the cell survival of oxic and anoxic cells irradiated by carbon ions

Light-ion radiation therapy against hypoxic tumors is highly curative due to reduced dependence on the presence of oxygen in the tumor at elevated linear energy transfer (LET) towards the Bragg peak. Clinical ion beams using spread-out Bragg peak (SOBP) are characterized by a wide spectrum of LET values. Accurate treatment optimization requires a method that can account for influence of the variation in response for a broad range of tumor hypoxia, absorbed doses and LETs. This paper presents a parameterization of the Repairable Conditionally-Repairable (RCR) cell survival model that can describe the survival of oxic and hypoxic cells over a wide range of LET values, and investigates the relationship between hypoxic radiation resistance and LET. The biological response model was tested by fitting cell survival data under oxic and anoxic conditions for V79 cells irradiated with LETs within the range of 30-500 keV/µm. The model provides good agreement with experimental cell survival data for the range of LET investigated, confirming the robustness of the parameterization method. This new version of the RCR model is suitable for describing the biological response of mixed populations of oxic and hypoxic cells and at the same time taking into account the distribution of doses and LETs in the incident beam and its variation with depth in tissue. The model offers a versatile tool for the selection of LET and dose required in the optimization of the therapeutic effect, without severely affecting normal tissue in realistic tumors presenting highly heterogeneous oxic and hypoxic regions.

Modelling of the oxygen enhancement ratio for ion beam radiation therapy

The poor treatment prognosis for tumours with oxygen-deficient areas is usually attributed to the increased radioresistance of hypoxic cells. It can be expressed by the oxygen enhancement ratio (OER), which decreases with increasing linear energy transfer (LET) suggesting a potential clinical advantage of high-LET radiotherapy with heavy ion beams compared to low-LET photon or proton irradiation. The aim of this work is to review the experimental cell survival data from the literature and, based on them, to develop a simple OER model to estimate the clinical impact of OER variations. For this purpose, the standard linear-quadratic model and the Alper-Howard-Flanders model are used. According to our calculations for a carbon ion spread-out Bragg peak at clinically relevant intermediate oxygen levels (0.5-20 mmHg), the advantage of carbon ions might be relatively moderate, with OER values about 1%-15% smaller than for protons. Furthermore, the variations of OER with LET are much smaller in vivo than in vitro due to different oxygen partial pressures used in cell experiments or measured inside tumours. The proposed OER model is a simple tool to quantify the oxygen effect in a practical way and provides the possibility to do hypoxia-based biological optimization in treatment planning.

Relative biological effectiveness of carbon ions for local tumor control of a radioresistant prostate carcinoma in the rat

PURPOSE

To study the relative biological effectiveness (RBE) of carbon ion beams relative to X-rays for local tumor control in a syngeneic rat prostate tumor (Dunning subline R3327-AT1).

METHODS AND MATERIALS

A total of 198 animals with tumors in the distal thigh were treated with increasing single and split doses of either (12)C ions or photons using a 20-mm spread-out Bragg peak. Endpoints of the study were local control (no tumor recurrence within 300 days) and volumetric changes after irradiation. The resulting values for D(50) (dose at 50% tumor control probability) were used to determine RBE values.

RESULTS

The D(50) values for single doses were 32.9 ± 0.9 Gy for (12)C ions and 75.7 ± 1.6 Gy for photons. The respective values for split doses were 38.0 ± 2.3 Gy and 90.6 ± 2.3 Gy. The corresponding RBE values were 2.30 ± 0.08 for single and 2.38 ± 0.16 for split doses. The most prominent side effects were dry and moist desquamation of the skin, which disappeared within weeks.

CONCLUSION

The study confirmed the effectiveness of carbon ion therapy for severely radioresistant tumors. For 1- and 2-fraction photon and (12)C ion radiation, we have established individual D(50) values for local tumor control as well as related RBE values.

Radiobiological characterization of human tumor cell multilayers after conventional and particle irradiation

The goal of this study was to establish planar multilayers from human tumor cells (WiDr and SiHa) as a model for irradiation of solid tumors. In addition to using conventional X rays (250 kV) as a reference standard, multilayers were tested for their suitability in cell survival studies with heavy-ion irradiation ((12)C(6+)) in the plateau and the extended Bragg peak with a scanned ion beam. Multilayers of both cell lines showed decreased survival compared to the corresponding monolayers after both X and heavy-ion irradiation. This multicellular sensitization effect is in contrast to the multicellular resistance or contact effect commonly described in the literature. Flow cytometry measurements showed an arrest of irradiated SiHa cells in G(2)/M phase. In contrast to the transient arrest of the monolayers, the multilayers stayed in a prolonged arrest. After Bragg-peak irradiation of monolayers, the arrest time was increased by 12-24 h, and more cells were arrested than with X rays. For multilayers, there were no differences between G(2) arrest after X rays and heavy ions for the entire observation period.

Resistance of hypoxic cells to ionizing radiation is influenced by homologous recombination status

PURPOSE

To determine the role of DNA repair in hypoxic radioresistance.

METHODS AND MATERIALS

Chinese hamster cell lines with mutations in homologous recombination (XRCC2, XRCC3, BRAC2, RAD51C) or nonhomologous end-joining (DNA-PKcs) genes were irradiated under normoxic (20% oxygen) and hypoxic (<0.1% oxygen) conditions, and the oxygen enhancement ratio (OER) was calculated. In addition, Fanconi anemia fibroblasts (complementation groups C and G) were compared with fibroblasts from nonsyndrome patients. RAD51 foci were studied using immunofluorescence.

RESULTS

All hamster cell lines deficient in homologous recombination showed a decrease in OER (1.5-2.0 vs. 2.6-3.0 for wild-types). In contrast, the OER for the DNA-PKcs-deficient line was comparable to wild-type controls. The two Fanconi anemia cell strains also showed a significant reduction in OER. The OER for RAD51 foci formation at late times after irradiation was considerably lower than that for survival in wild-type cells.

CONCLUSION

Homologous recombination plays an important role in determining hypoxic cell radiosensitivity. Lower OERs have also been reported in cells deficient in XPF and ERCC1, which, similar to homologous recombination genes, are known to play a role in cross-link repair. Because Fanconi anemia cells are also sensitive to cross-linking agents, this strengthens the notion that the capacity to repair cross-links determines hypoxic radiosensitivity.