Toward a Transportable Cell Culture Platform for Evaluating Radiotherapy Dose Modifying Factors

The current tools for validating dose delivery and optimizing new radiotherapy technologies in radiation therapy do not account for important dose modifying factors (DMFs), such as variations in cellular repair capability, tumor oxygenation, ultra-high dose rates and the type of ionizing radiation used. These factors play a crucial role in tumor control and normal tissue complications. To address this need, we explored the feasibility of developing a transportable cell culture platform (TCCP) to assess the relative biological effectiveness (RBE) of ionizing radiation. We measured cell recovery, clonogenic viability and metabolic viability of MDA-MB-231 cells over several days at room temperature in a range of concentrations of fetal bovine serum (FBS) in medium-supplemented gelatin, under both normoxic and hypoxic oxygen environments. Additionally, we measured the clonogenic viability of the cells to characterize how the duration of the TCCP at room temperature affected their radiosensitivity at doses up to 16 Gy. We found that (78±2)% of MDA-MB-231 cells were successfully recovered after being kept at room temperature for three days in 50% FBS in medium-supplemented gelatin at hypoxia (0.4±0.1)% pO2, while metabolic and clonogenic viabilities as measured by ATP luminescence and colony formation were found to be (58±5)% and (57±4)%, respectively. Additionally, irradiating a TCCP under normoxic and hypoxic conditions yielded a clonogenic oxygen enhancement ratio (OER) of 1.4±0.6 and a metabolic OER of 1.9±0.4. Our results demonstrate that the TCCP can be used to assess the RBE of a DMF and provides a feasible platform for assessing DMFs in radiation therapy applications.

Combined RBE and OER optimization in proton therapy with FLUKA based on EF5-PET

INTRODUCTION

Tumor hypoxia is associated with poor treatment outcome. Hypoxic regions are more radioresistant than well-oxygenated regions, as quantified by the oxygen enhancement ratio (OER). In optimization of proton therapy, including OER in addition to the relative biological effectiveness (RBE) could therefore be used to adapt to patient-specific radioresistance governed by intrinsic radiosensitivity and hypoxia.

METHODS

A combined RBE and OER weighted dose (ROWD) calculation method was implemented in a FLUKA Monte Carlo (MC) based treatment planning tool. The method is based on the linear quadratic model, with α and β parameters as a function of the OER, and therefore a function of the linear energy transfer (LET) and partial oxygen pressure (pO2 ). Proton therapy plans for two head and neck cancer (HNC) patients were optimized with pO2 estimated from [18 F]-EF5 positron emission tomography (PET) images. For the ROWD calculations, an RBE of 1.1 (RBE1.1,OER ) and two variable RBE models, Rørvik (ROR) and McNamara (MCN), were used, alongside a reference plan without incorporation of OER (RBE1.1 ).

RESULTS

For the HNC patients, treatment plans in line with the prescription dose and with acceptable target ROWD could be generated with the established tool. The physical dose was the main factor modulated in the ROWD. The impact of incorporating OER during optimization of HNC patients was demonstrated by the substantial difference found between ROWD and physical dose in the hypoxic tumor region. The largest physical dose differences between the ROWD optimized plans and the reference plan was 12.2 Gy.

CONCLUSION

The FLUKA MC based tool was able to optimize proton treatment plans taking the tumor pO2 distribution from hypoxia PET images into account. Independent of RBE-model, both elevated LET and physical dose were found in the hypoxic regions, which shows the potential to increase the tumor control compared to a conventional optimization approach.

Event-by-event approach to the oxygen-effect-incorporated stochastic microdosimetric kinetic model for hypofractionated multi-ion therapy

An oxygen-effect-incorporated stochastic microdosimetric kinetic (OSMK) model was previously developed to estimate the survival fraction of cells exposed to charged-particle beams with wide dose and linear energy transfer (LET) ranges under various oxygen conditions. In the model, hypoxia-induced radioresistance was formulated based on the dose-averaged radiation quality. This approximation may cause inaccuracy in the estimation of the biological effectiveness of the radiation with wide variation in energy deposited to a sensitive volume per event, such as spread-out Bragg peak (SOBP) beams. The purpose of this study was to apply an alternative approach so as to consider the energy depositions on an event-by-event basis. The production probability of radiation-induced lesions per energy was formulated with oxygen partial pressure to account for the hypoxia-induced radioresistance. The reduction in the oxygen enhancement ratio for high-LET radiations was modeled by reducing the sensitive-volume size and increasing the saturation energy in microdosimetry. The modified OSMK model was tested against the reported survival data of three cell lines exposed to six species of ions with wide dose and LET ranges under aerobic and hypoxic conditions. The model reasonably reproduced the reported cell survival data. To evaluate the event-by-event approach, survival distributions of Chinese hamster ovary cells exposed to SOBP beams were estimated using the original and modified OSMK models. The differences in the estimated survival distributions between the models were marginal even under extreme hypoxia. The event-by-event approach improved the theoretical validity of the OSMK model. However, the original OSMK model can still provide an accurate estimation of the biological effectiveness of therapeutic radiations.

Hypoxia adapted relative biological effectiveness models for proton therapy: a simulation study

In proton therapy, a constant relative biological effectiveness (RBE) factor of 1.1 is applied although the RBE has been shown to depend on factors including the Linear Energy Transfer (LET). The biological effectiveness of radiotherapy has also been shown to depend on the level of oxygenation, quantified by the oxygen enhancement ratio (OER). To estimate the biological effectiveness across different levels of oxygenation the RBE-OER-weighted dose (ROWD) can be used. To investigate the consistency between different approaches to estimate ROWD, we implemented and compared OER models in a Monte Carlo (MC) simulation tool. Five OER models were explored: Wenzl and Wilkens 2011 (WEN), Tinganelliet al2015 (TIN), Strigariet al2018 (STR), Dahleet al2020 (DAH) and Meinet al2021 (MEI). OER calculations were combined with a proton RBE model and the microdosimetric kinetic model for ROWD calculations. ROWD and OER were studied for a water phantom scenario and a head and neck cancer case using hypoxia PET data for the OER calculation. The OER and ROWD estimates from the WEN, MEI and DAH showed good agreement while STR and TIN gave higher OER values and lower ROWD. The WEN, STR and DAH showed some degree of OER-LET dependency while this was negligible for the MEI and TIN models. The ROWD for all implemented models is reduced in hypoxic regions with an OER of 1.0-2.1 in the target volume. While some variations between the models were observed, all models display a large difference in the estimated dose from hypoxic and normoxic regions. This shows the potential to increase the dose or LET in hypoxic regions or reduce the dose to normoxic regions which again could lead to normal tissue sparing. With reliable hypoxia imaging, RBE-OER weighting could become a useful tool for proton therapy plan optimization.

Mechanistic modelling of oxygen enhancement ratio of radiation via Monte Carlo simulation-based DNA damage calculation

Objective.Oxygen plays an important role in affecting the cellular radio-sensitivity to ionizing radiation. The objective of this study is to build a mechanistic model to compute oxygen enhancement ratio (OER) using a GPU-based Monte Carlo (MC) simulation package gMicroMC for microscopic radiation transport simulation and DNA damage calculation.Approach.We first simulated the water radiolysis process in the presence of DNA and oxygen for 1 ns and recorded the produced DNA damages. In this process, chemical reactions among oxygen, water radiolysis free radicals and DNA molecules were considered. We then applied a probabilistic approach to model the reactions between oxygen and indirect DNA damages for a maximal reaction time oft0. Finally, we defined two parametersP0andP1, representing probabilities for DNA damages without and with oxygen fixation effect not being restored in the repair process, to compute the final DNA double strand breaks (DSBs). As cell survival fraction is mainly determined by the number of DSBs, we assumed that the same numbers of DSBs resulted in the same cell survival rates, which enabled us to compute the OER as the ratio of doses producing the same number of DSBs without and with oxygen. We determined the three parameters (t0,P0andP1) by fitting the OERs obtained in our computation to a set of published experimental data under x-ray irradiation. We then validated the model by performing OER studies under proton irradiation and studied model sensitivity to parameter values.Main results.We obtained the model parameters ast0= 3.8 ms,P0= 0.08, andP1= 0.28 with a mean difference of 3.8% between the OERs computed by our model and that obtained from experimental measurements under x-ray irradiation. Applying the established model to proton irradiation, we obtained OERs as functions of oxygen concentration, LET, and dose values, which generally agreed with published experimental data. The parameter sensitivity analysis revealed that the absolute magnitude of the OER curve relied on the values ofP0andP1, while the curve was subject to a horizontal shift when adjustingt0.Significance.This study developed a mechanistic model that fully relies on microscopic MC simulations to compute OER.

Oxygen enhancement ratios of cancer cells after exposure to intensity modulated x-ray fields: DNA damage and cell survival

Hypoxic cancer cells within solid tumours show radio-resistance, leading to malignant progression in fractionated radiotherapy. When prescribing dose to tumours under heterogeneous oxygen pressure with intensity-modulated radiation fields, intercellular signalling could have an impact on radiosensitivity between in-field and out-of-field (OF) cells. However, the impact of hypoxia on radio-sensitivity under modulated radiation intensity remains to be fully clarified. Here, we investigate the impact of hypoxia on in-field and OF radio-sensitivities using two types of cancer cells, DU145 and H1299. Using a nBIONIX hypoxic culture kit and a shielding technique to irradiate 50% of a cell culture flask, oxygen enhancement ratios for double-strand breaks (DSB) and cell death endpoints were determined. Thesein vitromeasurements indicate that hypoxia impacts OF cells, although the hypoxic impacts on OF cells for cell survival were dose-dependent and smaller compared to those for in-field and uniformly irradiated cells. These decreased radio-sensitivities of OF cells were shown as a consistent tendency for both DSB and cell death endpoints, suggesting that radiation-induced intercellular communication is of importance in advanced radiotherapy dose-distributions such as with intensity-modulated radiotherapy.

Patient-Derived Nasopharyngeal Cancer Organoids for Disease Modeling and Radiation Dose Optimization

Effective radiation treatment (RT) for recurrent nasopharyngeal cancers (NPC), featuring an intrinsic hypoxic sub-volume, remains a clinical challenge. Lack of disease‐specific in-vitro models of NPC, together with difficulties in establishing patient derived xenograft (PDX) models, have further hindered development of personalized therapeutic options. Herein, we established two NPC organoid lines from recurrent NPC PDX models and further characterized and compared these models with original patient tumors using RNA sequencing analysis. Organoids were cultured in hypoxic conditions to examine the effects of hypoxia and radioresistance. These models were then utilized to determine the radiobiological parameters, such as α/β ratio and oxygen enhancement ratio (OER), characteristic to radiosensitive normoxic and radioresistant hypoxic NPC, using simple dose-survival data analytic tools. The results were further validated in-vitro and in-vivo, to determine the optimal boost dose and fractionation regimen required to achieve effective NPC tumor regression. Despite the differences in tumor microenvironment due to the lack of human stroma, RNA sequencing analysis revealed good correlation of NPC PDX and organoid models with patient tumors. Additionally, the established models also mimicked inter-tumoral heterogeneity. Hypoxic NPC organoids were highly radioresistant and had high α/β ratio compared to its normoxic counterparts. In-vitro and in-vivo fractionation studies showed that hypoxic NPC was less sensitive to RT fractionation scheme and required a large bolus dose or 1.4 times of the fractionated dose that was effective against normoxic cells in order to compensate for oxygen deficiency. This study is the first direct experimental evidence to predict optimal RT boost dose required to cause sufficient damage to recurrent hypoxic NPC tumor cells, which can be further used to develop dose-painting algorithms in clinical practice.

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.

RBE and OER within the spread-out Bragg peak for proton beam therapy: in vitro study at the Proton Medical Research Center at the University of Tsukuba

There are few reports on the biological homogeneity within the spread-out Bragg peak (SOBP) of proton beams. Therefore, to evaluate the relative biological effectiveness (RBE) and the oxygen enhancement ratio (OER), human salivary gland tumor (HSG) cells were irradiated at the plateau position (position A) and three different positions within a 6-cm-wide SOBP (position B, 26 mm proximal to the middle; position C, middle; position D, 26 mm distal to the middle) using 155-MeV/n proton beams under both normoxic and hypoxic conditions at the Proton Medical Research Center, University of Tsukuba, Japan. The RBE to the plateau region (RBE(plateau)) and the OER value were calculated from the doses corresponding to 10% survival data. Under the normoxic condition, the RBE(plateau) was 1.00, 0.99 and 1.09 for positions B, C and D, respectively. Under the hypoxic condition, the RBE(plateau) was 1.10, 1.06 and 1.12 for positions B, C and D, respectively. The OER was 2.84, 2.60, 2.63 and 2.76 for positions A, B, C and D, respectively. There were no significant differences in either the RBE(plateau) or the OER between these three positions within the SOBP. In conclusion, biological homogeneity need not necessarily be taken into account for treatment planning for proton beam therapy at the University of Tsukuba.

Theoretical analysis of the dose dependence of the oxygen enhancement ratio and its relevance for clinical applications

BACKGROUND

The increased resistance of hypoxic cells to ionizing radiation is usually believed to be the primary reason for treatment failure in tumors with oxygen-deficient areas. This oxygen effect can be expressed quantitatively by the oxygen enhancement ratio (OER). Here we investigate theoretically the dependence of the OER on the applied local dose for different types of ionizing irradiation and discuss its importance for clinical applications in radiotherapy for two scenarios: small dose variations during hypoxia-based dose painting and larger dose changes introduced by altered fractionation schemes.

METHODS

Using the widespread Alper-Howard-Flanders and standard linear-quadratic (LQ) models, OER calculations are performed for T1 human kidney and V79 Chinese hamster cells for various dose levels and various hypoxic oxygen partial pressures (pO2) between 0.01 and 20 mmHg as present in clinical situations in vivo. Our work comprises the analysis for both low linear energy transfer (LET) treatment with photons or protons and high-LET treatment with heavy ions. A detailed analysis of experimental data from the literature with respect to the dose dependence of the oxygen effect is performed, revealing controversial opinions whether the OER increases, decreases or stays constant with dose.

RESULTS

The behavior of the OER with dose per fraction depends primarily on the ratios of the LQ parameters alpha and beta under hypoxic and aerobic conditions, which themselves depend on LET, pO2 and the cell or tissue type. According to our calculations, the OER variations with dose in vivo for low-LET treatments are moderate, with changes in the OER up to 11% for dose painting (1 or 3 Gy per fraction compared to 2 Gy) and up to 22% in hyper-/hypofractionation (0.5 or 20 Gy per fraction compared to 2 Gy) for oxygen tensions between 0.2 and 20 mmHg typically measured clinically in hypoxic tumors. For extremely hypoxic cells (0.01 mmHg), the dose dependence of the OER becomes more pronounced (up to 36%). For high LET, OER variations up to 4% for the whole range of oxygen tensions between 0.01 and 20 mmHg were found, which were much smaller than for low LET.

CONCLUSIONS

The formalism presented in this paper can be used for various tissue and radiation types to estimate OER variations with dose and help to decide in clinical practice whether some dose changes in dose painting or in fractionation can bring more benefit in terms of the OER in the treatment of a specific hypoxic tumor.