Relative biological effectiveness for cell-killing effect on various human cell lines irradiated with heavy-ion medical accelerator in Chiba (HIMAC) carbon-ion beams

LET-based approximation of the microdosimetric kinetic model for proton radiotherapy

BACKGROUND

Phenomenological relative biological effectiveness (RBE) models for proton therapy, based on the dose-averaged linear energy transfer (LET), have been developed to address the apparent RBE increase towards the end of the proton range. The results of these phenomenological models substantially differ due to varying empirical assumptions and fitting functions. In contrast, more theory-based approaches are used in carbon ion radiotherapy, such as the microdosimetric kinetic model (MKM). However, implementing microdosimetry-based models in LET-based proton therapy treatment planning systems poses challenges.

PURPOSE

This work presents a LET-based version of the MKM that is practical for clinical use in proton radiotherapy.

METHODS

At first, we derived an approximation of the Mayo Clinic Florida (MCF) MKM for relatively-sparsely ionizing radiation such as protons. The mathematical formalism of the proposed model is equivalent to the original MKM, but it maintains some key features of the MCF MKM, such as the determination of model parameters from measurable cell characteristics. Subsequently, we carried out Monte Carlo calculations with PHITS in different simulated scenarios to establish a heuristic correlation between microdosimetric quantities and the dose averaged LET of protons.

RESULTS

A simple allometric function was found able to describe the relationship between the dose-averaged LET of protons and the dose-mean lineal energy, which includes the contributions of secondary particles. The LET-based MKM was used to model the in vitro clonogenic survival RBE of five human and rodent cell lines (A549, AG01522, CHO, T98G, and U87) exposed to pristine and spread-out Bragg peak (SOBP) proton beams. The results of the LET-based MKM agree well with the biological data in a comparable or better way with respect to the other models included in the study. A sensitivity analysis on the model results was also performed.

CONCLUSIONS

The LET-based MKM integrates the predictive theoretical framework of the MCF MKM with a straightforward mathematical description of the RBE based on the dose-averaged LET, a physical quantity readily available in modern treatment planning systems for proton therapy.

Understanding Relative Biological Effectiveness and Clinical Outcome of Prostate Cancer Therapy Using Particle Irradiation: Analysis of Tumor Control Probability With the Modified Microdosimetric Kinetic Model

PURPOSE

Recent experimental studies and clinical trial results might indicate that-at least for some indications-continued use of the mechanistic model for relative biological effectiveness (RBE) applied at carbon ion therapy facilities in Europe for several decades (LEM-I) may be unwarranted. We present a novel clinical framework for prostate cancer treatment planning and tumor control probability (TCP) prediction based on the modified microdosimetric kinetic model (mMKM) for particle therapy.

METHODS AND MATERIALS

Treatment plans of 91 patients with prostate tumors (proton: 46, carbon ions: 45) applying 66 GyRBE [RBE = 1.1 for protons and LEM-I, (α/β)x = 2.0 Gy, for carbon ions] in 20 fractions were recalculated using mMKM [(α/β)x = 3.1 Gy]). Based solely on the response data of photon-irradiated patient groups stratified according to risk and usage of androgen deprivation therapy, we derived parameters for an mMKM-based Poisson-TCP model. Subsequently, new carbon and helium ion plans, adhering to prescribed biological dose criteria, were generated. These were systematically compared with the clinical experience of Japanese centers employing an analogous fractionation scheme and existing proton plans.

RESULTS

mMKM predictions suggested significant biological dose deviation between the proton and carbon ion arms. Patients irradiated with protons received (3.25 ± 0.08)  GyRBEmMKM/Fx, whereas patients treated with carbon ions received(2.51 ± 0.05)  GyRBEmMKM/Fx. TCP predictions were (86 ± 3)% for protons and (52 ± 4)% for carbon ions, matching the clinical outcome of 85% and 50%. Newly optimized carbon ion plans, guided by the mMKM/TCP model, effectively replicated clinical data from Japanese centers. Using mMKM, helium ions exhibited similar target coverage as proton and carbon ions and improved rectum and bladder sparing compared with proton.

CONCLUSIONS

Our mMKM-based model for prostate cancer treatment planning and TCP prediction was validated against clinical data for proton and carbon ion therapy, and its application was extended to helium ion therapy. Based on the data presented in this work, mMKM seems to be a good candidate for clinical biological calculations in carbon ion therapy for prostate cancer.

Research Trends in the Study of the Relative Biological Effectiveness: A Bibliometric Study

The relative biological effectiveness is a mathematical quantity first defined in the 1950s. This has resulted in more than 4,000 scientific papers published to date. Yet defining the correct value of the RBE to use in clinical practice remains a challenge. A scientific analysis in the radiation research literature can provide an understanding of how this mathematical quantity has evolved. The purpose of this study is to investigate documents published since 1950 using bibliometric indicators and network visualization. This analysis seeks to provide an assessment of global research activities, key themes, and RBE research within the radiation-related field. It strives to highlight top-performing authors, organizations, and nations that have made major contributions to this research domain, as well as their interactions. The Scopus Collection was searched for articles and reviews pertaining to RBE in radiation research from 1950 through 2023. Scopus and Bibiometrix analytic tools were used to investigate the most productive countries, researchers, collaboration networks, journals, along with the citation analysis of references and keywords. A total of 4,632 documents were retrieved produced by authors originating from 71 countries. Publication trends could be separated in 20-year groupings beginning with slow accrual from 1950 to 1970, an early rise from 1970-1990, followed by a sharp increase in the years 1990s-2010s that matches the development of charged particle therapy in clinics worldwide and opened discussion on the true value of the RBE in proton beam therapy. Since the 2010s, a steady 200 papers, on average, have been published per year. The United States produced the most publications overall (N = 1,378) and Radiation Research was the most likely journal to have published articles related to the RBE (606 publications during this period). J. Debus was the most prolific author (112 contributions, with 2,900 citations). The RBE has captured the research interest of over 7,000 authors in the past decade alone. This study supports that notion that the growth of the body of evidence surrounding the RBE, which started 75 years ago, is far from reaching its end. Applications to medicine have continuously dominated the field, with physics competing with Biochemistry, Genetics and Molecular Biology for second place over the decades. Future research can be predicted to continue.

Carbon-ion radiotherapy for hepatocellular carcinoma with major vascular invasion: a retrospective cohort study

BACKGROUND

Macroscopic vascular invasion (MVI) significantly impacts survival in patients with hepatocellular carcinoma (HCC), warranting systemic therapy over locoregional therapy. Despite novel approaches, HCC with MVI has a poor prognosis compared to early-to intermediate-stage HCC. This study aimed to evaluate the safety and efficacy of carbon-ion radiotherapy (C-ion RT) for HCC characterized by MVI.

METHODS

This retrospective cohort study evaluated HCC patients with MVI treated using C-ion RT with a dose of 45.0-48.0 Gy/2 fractions or 52.8-60.0 Gy/4 fractions between 1995 and 2020 at our institution in Japan. We analyzed the prognostic factors and rates of local recurrence, survival, and adverse events. The local recurrence rate was determined using the cumulative incidence function, with death as a competing event. Survival rates were determined using the Kaplan-Meier method. The log-rank test for univariate analysis and the Cox proportional hazards model for multivariate analysis were used to compare subgroups.

RESULTS

In total, 76 patients with a median age of 71 years (range, 45-86 years) were evaluated. Among them, 68 had Child-Pugh grade A while eight had grade B disease. In 17 patients, the vascular tumor thrombus reached the inferior vena cava or main trunk of the portal vein. Over a median follow-up period of 27.9 months (range, 1.5-180.4 months), the 2-year overall survival, progression-free survival, and local recurrence rates were 70.0% (95% confidence interval [CI]: 57.7-79.4%), 32.7% (95% CI: 22.0-43.8%), and 8.9% (95% CI: 1.7-23.5%), respectively. A naïve tumor and a single lesion were significant prognostic factors for overall survival in the univariate analysis. Albumin-bilirubin grade 1 and a single lesion were independent prognostic factors in the multivariate analysis. Overall, four patients (5%) experienced grade 3 late adverse events, with no observed grade 4 or 5 acute or late adverse events.

CONCLUSIONS

C-ion RT for HCC with MVI showed favorable local control and survival benefits with minimal toxicity.

Most DNA repair defects do not modify the relationship between relative biological effectiveness and linear energy transfer in CRISPR-edited cells

BACKGROUND

Cancer is a highly heterogeneous disease, driven by frequent genetic alterations which have significant effects on radiosensitivity. However, radiotherapy for a given cancer type is typically given with a standard dose determined from population-level trials. As a result, a proportion of patients are under- or over-dosed, reducing the clinical benefit of radiotherapy. Biological optimization would not only allow individual dose prescription but also a more efficient allocation of limited resources, such as proton and carbon ion therapy. Proton and ion radiotherapy offer an advantage over photons due to their elevated Relative Biological Effectiveness (RBE) resulting from their elevated Linear Energy Transfer (LET). Despite significant interest in optimizing LET by tailoring radiotherapy plans, RBE's genetic dependence remains unclear.

PURPOSE

The aim of this study is to better define the RBE/LET relationship in a panel of cell lines with different defects in DSB repair pathways, but otherwise identical biological features and genetic background to isolate these effects.

METHODS

Normal human cells (RPE1), genetically modified to introduce defects in DNA double-strand break (DSB) repair genes, ATM, BRCA1, DCLRE1C, LIG4, PRKDC and TP53, were used to map the RBE-LET relationship. Cell survival was measured with clonogenic assays after exposure to photons, protons (LET 1 and 12 keV/µm) and alpha particles (129 keV/µm). Gene knockout sensitizer enhancement ratio (SER) values were calculated as the ratio of the mean inactivation dose (MID) of wild-type cells to repair-deficient cells, and RBE values were calculated as the ratio of the MID of X-ray and particle irradiated cells. 53BP1 foci were used to quantify radiation-induced DSBs and their repair following irradiation.

RESULTS

Deletion of NHEJ genes had the greatest impact on photon sensitivity (ATM-/- SER = 2.0 and Lig4-/- SER = 1.8), with genes associated with HR having smaller effects (BRCA1-/- SER = 1.2). Wild-type cells showed RBEs of 1.1, 1.3, 5.0 for low- and high-LET protons and alpha particles respectively. SERs for different genes were independent of LET, apart from NHEJ knockouts which proved to be markedly hypersensitive across all tested LETs. Due to this hypersensitivity, the impact of high LET was reduced in cell models lacking the NHEJ repair pathway. HR-defective cells had moderately increased sensitivity across all tested LETs, but, notably, the contribution of HR pathway to survival appeared independent of LET. Analysis of 53BP1 foci shows that NHEJ-defective cells had the least DSB repair capacity after low LET exposure, and no visible repair after high LET exposure. HR-defective cells also had slower repair kinetics, but the impact of HR defects is not as severe as NHEJ defects.

CONCLUSIONS

DSB repair defects, particularly in NHEJ, conferred significant radiosensitivity across all LETs. This sensitization appeared independent of LET, suggesting that the contribution of different DNA repair pathways to survival does not depend on radiation quality.

Enhancing Tissue Equivalence in 7Li Heavy Ion Therapy with MC Algorithm Optimized Polymer-Based Bioinks

The unique physical properties of heavy ion beams, particularly their distinctive depth-dose distribution and sharp lateral dose reduction profiles, have led to their widespread adoption in tumor therapy worldwide. However, the physical properties of heavy ion beams must be investigated to deliver a sufficient dose to tumors without damaging organs at risk. These studies should be performed on phantoms made of biomaterials that closely mimic human tissue. Polymers can serve as soft tissue substitutes and are suitable materials for building radiological phantoms due to their physical, mechanical, biological, and chemical properties. Extensive research, development, and applications of polymeric biomaterials have been encouraged due to these properties. In this study, we investigated the ionization, recoils, phonon release, collision events, and lateral straggle properties of polymeric biomaterials that closely resemble soft tissue using lithium-ion beams and Monte Carlo Transport of Ions in Matter simulation. The results indicated that the Bragg peak position closest to soft tissue was achieved with a 7.3% difference in polymethylmethacrylate, with an average recoils value of 10.5%. Additionally, average values of 33% were observed in collision events and 22.6% in lateral straggle. A significant contribution of this study to the existing literature lies in the exploration of secondary interactions alongside the assessment of linear energy transfer induced by the 7Li beam used for treatment. Furthermore, we analyzed the tissue-equivalent properties of polymer biomaterials using heavy ion beams, taking into account phonon release resulting from ionization, recoils, lateral straggle, and all other interactions. This approach allows for the evaluation of the most suitable polymeric biomaterials for heavy ion therapy while considering the full range of interactions involved.

Current evidence and the potential role of proton beam therapy for hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a leading cause of cancer-related death, and external beam radiation therapy has emerged as a promising approach for managing HCC. Proton beam therapy (PBT) offers dosimetric advantages over X-ray therapy, with superior physical properties known as the Bragg peak. PBT holds promise for reducing hepatotoxicity and allowing safe dose-escalation to the tumor. It has been tried in various clinical conditions and has shown promising local tumor control and survival outcomes. A recent phase III trial demonstrated the non-inferiority of PBT in local tumor control compared to current standard radiofrequency ablation in early-stage HCC. PBT also tended to show more favorable outcomes compared to transarterial chemoembolization in the intermediate stage, and has proven effective in-field disease control and safe toxicity profiles in advanced HCC. In this review, we discuss the rationale, clinical studies, optimal indication, and future directions of PBT in HCC treatment.

Chondrosarcoma Resistance to Radiation Therapy: Origins and Potential Therapeutic Solutions

Chondrosarcoma is a malignant cartilaginous tumor that is particularly chemoresistant and radioresistant to X-rays. The first line of treatment is surgery, though this is almost impossible in some specific locations. Such resistances can be explained by the particular composition of the tumor, which develops within a dense cartilaginous matrix, producing a resistant area where the oxygen tension is very low. This microenvironment forces the cells to adapt and dedifferentiate into cancer stem cells, which are described to be more resistant to conventional treatments. One of the main avenues considered to treat this type of tumor is hadrontherapy, in particular for its ballistic properties but also its greater biological effectiveness against tumor cells. In this review, we describe the different forms of chondrosarcoma resistance and how hadrontherapy, combined with other treatments involving targeted inhibitors, could help to better treat high-grade chondrosarcoma.

Treatment planning of carbon ion radiotherapy for prostate cancer based on cellular experiments with PC3 human prostate cancer cells

[Purpose] Treatment plans for carbon ion radiotherapy (CIRT) in Japan are designed to uniformly deliver the prescribed clinical dose based on the radiosensitivity of human salivary gland (HSG) cells to the planning target volume (PTV). However, sensitivity to carbon beams varies between cell lines, that is, it should be checked that the clinical dose distribution based on the cell radiosensitivity of the treatment site is uniform within the PTV. [Methods] We modeled the linear energy transfer (LET) dependence of the linear-quadratic (LQ) coefficients specific to prostate cancer, which accounts for the majority of CIRT. This was achieved by irradiating prostate cancer cells (PC3) with X-rays from a 4 MV-Linac and carbon beams with different LETs of 11.1-214.3 keV/μm. By using the radiosensitivity of PC3 cells derived from cellular experiments, we reconstructed prostate-cancer-specific clinical dose distributions on patient computed tomography (CT). [Results] The LQ coefficient, α, of PC3 cells was larger than that of HSG cells at low (<50 keV/μm) LET and smaller at high (>50 keV/μm) LET, which was validated by cellular experiments performed on rectangular SOBPs. The reconstructed dose distribution on patient CT was sloped when 1 fraction incident from the one side of the patient was considered, but remained uniform from the sum of 12 fractions of the left-right opposing beams (as is used in clinical practice). [Conclusion] Our study reveals the inhomogeneity of clinical doses in single-field plans calculated using the PC3 radiosensitivity data. However, this inhomogeneity is compensated by using the combination of left-right opposing beams.

On the effectiveness of proton boron fusion therapy (PBFT) at cellular level

The present work introduced a framework to investigate the effectiveness of proton boron fusion therapy (PBFT) at the cellular level. The framework consisted of a cell array generator program coupled with PHITS Monte Carlo package with a dedicated terminal-based code editor that was developed in this work. The framework enabled users to model large cell arrays with normal, all boron, and random boron filled cytoplasm, to investigate the underlying mechanism of PBFT. It was found that alpha particles and neutrons could be produced in absence of boron mainly because of nuclear reaction induced by proton interaction with ¹⁶O, ¹²C and ¹⁴N nuclei. The effectiveness of PBFT is highly dependent on the incident proton energy, source size, cell array size, buffer medium thickness layer, concentration and distribution of boron in the cell array. To quantitatively assess the effectiveness of PBFT, of the total energy deposition by alpha particle for different cases were determined. The number of alpha particle hits in cell cytoplasm and nucleus for normal and 100 ppm boron were determined. The obtained results and the developed tools would be useful for future development of PBFT to objectively determine the effectiveness of this treatment modality.

The Mayo Clinic Florida Microdosimetric Kinetic Model of Clonogenic Survival: Application to Various Repair-Competent Rodent and Human Cell Lines

The relative biological effectiveness (RBE) calculations used during the planning of ion therapy treatments are generally based on the microdosimetric kinetic model (MKM) and the local effect model (LEM). The Mayo Clinic Florida MKM (MCF MKM) was recently developed to overcome the limitations of previous MKMs in reproducing the biological data and to eliminate the need for ion-exposed in vitro data as input for the model calculations. Since we are considering to implement the MCF MKM in clinic, this article presents (a) an extensive benchmark of the MCF MKM predictions against corresponding in vitro clonogenic survival data for 4 rodent and 10 cell lines exposed to ions from 1H to 238U, and (b) a systematic comparison with published results of the latest version of the LEM (LEM IV). Additionally, we introduce a novel approach to derive an approximate value of the MCF MKM model parameters by knowing only the animal species and the mean number of chromosomes. The overall good agreement between MCF MKM predictions and in vitro data suggests the MCF MKM can be reliably used for the RBE calculations. In most cases, a reasonable agreement was found between the MCF MKM and the LEM IV.

An empirical model of proton RBE based on the linear correlation between x-ray and proton radiosensitivity

BACKGROUND

Proton relative biological effectiveness (RBE) is known to depend on physical factors of the proton beam, such as its linear energy transfer (LET), as well as on cell-line specific biological factors, such as their ability to repair DNA damage. However, in a clinical setting, proton RBE is still considered to have a fixed value of 1.1 despite the existence of several empirical models that can predict proton RBE based on how a cell's survival curve (linear-quadratic model [LQM]) parameters α and β vary with the LET of the proton beam. Part of the hesitation to incorporate variable RBE models in the clinic is due to the great noise in the biological datasets on which these models are trained, often making it unclear which model, if any, provides sufficiently accurate RBE predictions to warrant a departure from RBE = 1.1.

PURPOSE

Here, we introduce a novel model of proton RBE based on how a cell's intrinsic radiosensitivity varies with LET, rather than its LQM parameters.

METHODS AND MATERIALS

We performed clonogenic cell survival assays for eight cell lines exposed to 6 MV x-rays and 1.2, 2.6, or 9.9 keV/µm protons, and combined our measurements with published survival data (n = 397 total cell line/LET combinations). We characterized how radiosensitivity metrics of the form DSF% , (the dose required to achieve survival fraction [SF], e.g., D10% ) varied with proton LET, and calculated the Bayesian information criteria associated with different LET-dependent functions to determine which functions best described the underlying trends. This allowed us to construct a six-parameter model that predicts cells' proton survival curves based on the LET dependence of their radiosensitivity, rather than the LET dependence of the LQM parameters themselves. We compared the accuracy of our model to previously established empirical proton RBE models, and implemented our model within a clinical treatment plan evaluation workflow to demonstrate its feasibility in a clinical setting.

RESULTS

Our analyses of the trends in the data show that DSF% is linearly correlated between x-rays and protons, regardless of the choice of the survival level (e.g., D10% , D37% , or D50% are similarly correlated), and that the slope and intercept of these correlations vary with proton LET. The model we constructed based on these trends predicts proton RBE within 15%-30% at the 68.3% confidence level and offers a more accurate general description of the experimental data than previously published empirical models. In the context of a clinical treatment plan, our model generally predicted higher RBE-weighted doses than the other empirical models, with RBE-weighted doses in the distal portion of the field being up to 50.7% higher than the planned RBE-weighted doses (RBE = 1.1) to the tumor.

CONCLUSIONS

We established a new empirical proton RBE model that is more accurate than previous empirical models, and that predicts much higher RBE values in the distal edge of clinical proton beams.

Geometrical Properties of the Nucleus and Chromosome Intermingling Are Possible Major Parameters of Chromosome Aberration Formation

Ionizing radiation causes chromosome aberrations, which are possible biomarkers to assess space radiation cancer risks. Using the Monte Carlo codes Relativistic Ion Tracks (RITRACKS) and Radiation-Induced Tracks, Chromosome Aberrations, Repair and Damage (RITCARD), we investigated how geometrical properties of the cell nucleus, irradiated with ion beams of linear energy transfer (LET) ranging from 0.22 keV/μm to 195 keV/μm, influence the yield of simple and complex exchanges. We focused on the effect of (1) nuclear volume by considering spherical nuclei of varying radii; (2) nuclear shape by considering ellipsoidal nuclei of varying thicknesses; (3) beam orientation; and (4) chromosome intermingling by constraining or not constraining chromosomes in non-overlapping domains. In general, small nuclear volumes yield a higher number of complex exchanges, as compared to larger nuclear volumes, and a higher number of simple exchanges for LET < 40 keV/μm. Nuclear flattening reduces complex exchanges for high-LET beams when irradiated along the flattened axis. The beam orientation also affects yields for ellipsoidal nuclei. Reducing chromosome intermingling decreases both simple and complex exchanges. Our results suggest that the beam orientation, the geometry of the cell nucleus, and the organization of the chromosomes within are important parameters for the formation of aberrations that must be considered to model and translate in vitro results to in vivo risks.

Biological Effectiveness of Ion Beam for In Vitro Cell Irradiations

Despite numerous ion beam irradiation of cell experiments performed over the past five decades, the relationship between the biological effectiveness of ion beams and the physical characteristics of the ion beam remains unclear. Using 1,118 sets of in vitro cell survival experiments with ion beam irradiation, compiled by the Particle Irradiation Data Ensemble (PIDE) project, the relationship between cell survival and the fluence and linear energy transfer (LET) of the ion beam was established. Unlike previous studies, the closed-form analytical function is independent of photon irradiation and takes a universal form across all ion and cell species. A new understanding of the biological effectiveness of ion beams is crucial for predicting tumor response and toxicities in ion beam radiation therapy, along with radiation protection for high-LET ion beams with low fluence.

Biological effectiveness and relative biological effectiveness of ion beams for in-vitro cell irradiation

Biological effectiveness and relative biological effectiveness are critical for proton and ion beam radiotherapy. However, the relationship between the two quantities and physical character of ion beams is not well established. By analyzing 1188 sets of in‐vitro cell irradiation experiments using ion beams ranging from protons to ²³⁸U, compiled by the Particle Irradiation Data Ensemble (PIDE) project, the biological effectiveness of the ion beams, with cell survival fractionation (SF) as the endpoint, was found to be dependent on the fluence and linear energy transfer (LET) of the ion beam. Consequently, the relative biological effectiveness of the ion beam to photon beam was also established as a function of LET. A common form of relationship among SF, fluence, and LET was found to be valid for all ion beam experiments. The close form relationship could be used for proton and ion beam radiotherapy applications.

A Consistent Protocol Reveals a Large Heterogeneity in the Biological Effectiveness of Proton and Carbon-Ion Beams for Various Sarcoma and Normal-Tissue-Derived Cell Lines

Simple Summary

Using a consistent experimental protocol, we found a large heterogeneity in the relative biological effectiveness (RBE) values of both proton and carbon-ion beams in various sarcomas and normal-tissue-derived cell lines. Our data suggest that proton beam therapy may be more beneficial for some types of tumors. In carbon-ion therapy, for some types of tumors, large heterogeneity in RBE should prompt consideration of dose reduction or an increased dose per fraction. In particular, a higher RBE value in normal tissues requires caution. Specific dose evaluations for tumor and normal tissues are needed for both proton and carbon-ion therapies.

Abstract

This study investigated variations in the relative biological effectiveness (RBE) values among various sarcoma and normal-tissue-derived cell lines (normal cell line) in proton beam and carbon-ion irradiations. We used a consistent protocol that specified the timing of irradiation after plating cells and detailed the colony formation assay. We examined the cell type dependence of RBE for proton beam and carbon-ion irradiations using four human sarcoma cell lines (MG63 osteosarcoma, HT1080 fibrosarcoma, SW872 liposarcoma, and SW1353 chondrosarcoma) and three normal cell lines (HDF human dermal fibroblast, hTERT-HME1 mammary gland, and NuLi-1 bronchus epithelium). The cells were irradiated with gamma rays, proton beams at the center of the spread-out Bragg peak, or carbon-ion beams at 54.4 keV/μm linear energy transfer. In all sarcoma and normal cell lines, the average RBE values in proton beam and carbon-ion irradiations were 1.08 ± 0.11 and 2.08 ± 0.36, which were consistent with the values of 1.1 and 2.13 used in current treatment planning systems, respectively. Up to 34% difference in the RBE of the proton beam was observed between MG63 and HT1080. Similarly, a 32% difference in the RBE of the carbon-ion beam was observed between SW872 and the other sarcoma cell lines. In proton beam irradiation, normal cell lines had less variation in RBE values (within 10%), whereas in carbon-ion irradiation, RBE values differed by up to 48% between hTERT-HME1 and NuLi-1. Our results suggest that specific dose evaluations for tumor and normal tissues are necessary for treatment planning in both proton and carbon-ion therapies.

Normal Tissue Complication Probability Model for Acute Radiation Dermatitis in Patients with Head and Neck Cancer Treated with Carbon Ion Radiotherapy

PURPOSE

This study aimed to explore the prognostic factors associated with acute radiation dermatitis (ARD). A normal tissue complication probability (NTCP) model for ARD in patients with head and neck cancer (HNC) treated with carbon ion radiotherapy (CIRT) was developed.

MATERIALS AND METHODS

A total of 187 patients were included in the analysis, and the endpoint was ≥grade 2 ARD. The biological and physical dose-surface parameters associated with ARD were used in the logistic regression model. The mean areas under the receiver operating characteristic curve (AUC) in the internal cross-validation and Akaike's corrected Information Criterion (AICc) were examined for model evaluation and selection. The multivariate logistic regression NTCP models were established based on factors with weak correlation.

RESULTS

Tumour volume, planning target volume to the skin, radiation technique and all dose-surface parameters were significantly associated with ARD (P < 0.05). Models with high performance for grade 2-3 ARD were constructed. The most significant prognostic predictors were S_{40 Gy (relative biological effectiveness, RBE)} and S_{20 Gy} [absolute surface area receiving RBE-weighted dose of 40 Gy (RBE) or physical dose of 20 Gy]. The internal cross-validation-based AUCs for models with S_{40 Gy (RBE)} and S_{20 Gy} were 0.78 and 0.77, respectively. The biological and physical dose-surface parameters had similar performance at various dose levels. However, the performance of the multivariate NTCP models based on two factors was not better than that of the univariate models.

CONCLUSIONS

NTCP models for ARD may provide a basis for the development of individualised treatment strategies and reduce the incidence of severe ARD in patients with HNC receiving CIRT. Furthermore, biological and physical dose-surface parameter-based models are comparable. However, further validation with more evaluation parameters is warranted.

Lethal DNA damage caused by ion-induced shock waves in cells

The elucidation of fundamental mechanisms underlying ion-induced radiation damage of biological systems is crucial for advancing radiotherapy with ion beams and for radiation protection in space. The study of ion-induced biodamage using the phenomenon-based multiscale approach (MSA) to the physics of radiation damage with ions has led to the prediction of nanoscale shock waves created by ions in a biological medium at the high linear energy transfer (LET). The high-LET regime corresponds to the keV and higher-energy losses by ions per nanometer, which is typical for ions heavier than carbon at the Bragg peak region in biological media. This paper reveals that the thermomechanical stress of the DNA molecule caused by the ion-induced shock wave becomes the dominant mechanism of complex DNA damage at the high-LET ion irradiation. Damage of the DNA molecule in water caused by a projectile-ion-induced shock wave is studied by means of reactive molecular dynamics simulations. Five projectile ions (carbon, oxygen, silicon, argon, and iron) at the Bragg peak energies are considered. For the chosen segment of the DNA molecule and the collision geometry, the number of DNA strand breaks is evaluated for each projectile ion as a function of the bond dissociation energy and the distance from the ion's path to the DNA strands. Simulations reveal that argon and especially iron ions induce the breakage of multiple bonds in a DNA double convolution containing 20 DNA base pairs. The DNA damage produced in segments of such size leads to complex irreparable lesions in a cell. This makes the shock-wave-induced thermomechanical stress the dominant mechanism of complex DNA damage at the high-LET ion irradiation. A detailed theory for evaluating the DNA damage caused by ions at high-LET is formulated and integrated into the MSA formalism. The theoretical analysis reveals that a single ion hitting a cell nucleus at high-LET is sufficient to produce highly complex, lethal damages to a cell by the shock-wave-induced thermomechanical stress. Accounting for the shock-wave-induced thermomechanical mechanism of DNA damage provides an explanation for the "overkill" effect observed experimentally in the dependence of cell survival probabilities on the radiation dose delivered with iron ions. This important observation provides strong experimental evidence of the ion-induced shock-wave effect and the related mechanism of radiation damage in cells.

Evaluation of Threshold Dose of Damaged Hepatic Tissue After Carbon-Ion Radiation Therapy Using Gd-EOB-DTPA-Enhanced Magnetic Resonance Imaging

Purpose

To evaluate the threshold dose and associated factors using signal-intensity changes in the irradiated area after carbon-ion radiation therapy (C-ion RT) for patients with liver cancer.

Methods and Materials

Patients treated for the first time with C-ion RT for malignant liver tumors and followed up with 3-Tesla gadoxetic acid (Gd-EOB-DTPA)–enhanced magnetic resonance imaging (MRI) 3 months after treatment completion were retrospectively enrolled. The volume of focal liver reaction (FLR), a low-intensity area in the hepatobiliary phase of Gd-EOB-DTPA after treatment, was measured. Corrected FLR (cFLR) volume, defined as FLR corrected for changes in tumor volume from before to after treatment, was calculated, and the threshold dose was determined by applying the cFLR volume in the dose-volume histogram. To evaluate potential mismatch in fusion images of planning computed tomography and follow-up MRI, the concordance coefficient (CC) was measured, and patients with a CC < 0.7 were excluded. Sixty patients were included. Multiple regression analysis was performed with the threshold dose as the objective variable and the age, dose, number of fractionations, Child-Pugh score, pretreatment liver volume, and pretreatment tumor volume as explanatory variables. The Student t test or Mann-Whitney U test was used as required.

Results

The median threshold doses for each number of dose fractionations (4 fractions, 12 fractions, and overall) were 51.6, 51.9, and 51.8 Gy (relative biological effectiveness [RBE]), respectively, in patients categorized as Child-Pugh class A and 27.0, 28.8, and 27.0 Gy (RBE), respectively, in patients categorized as Child-Pugh class B. In the multiple-regression analysis, only the Child-Pugh score was significant (P < .001). The number of dose fractionations was not statistically significant.

Conclusions

Although few patients in the study had decreased liver function, baseline liver function was the only factor significantly associated with the median threshold dose. These findings facilitate appropriate patient selection to receive C-ion RT for malignant hepatic tumors.

An Analysis of 20 Cases of Radiation-Associated Sarcoma, Including 4 Cases Treated by Carbon Ion Radiotherapy

Introduction Radiation-associated sarcoma (RAS) is one of the most life-threatening complications associated with the treatment of malignant neoplasms. Because all RAS patients have a history of radiotherapy, there have been no effective treatment options when RAS is not completely resected. Methods We retrospectively reviewed 20 RAS patients, including 4 unresectable cases treated by carbon ion radiotherapy (CIRT). Results The primary diseases targeted by radiotherapy included malignant lymphoma (n=4), cervical cancer (n=3), pharyngeal cancer (n=3), breast cancer (n=2), lung cancer (n=1), rectal cancer (n=1), maxillary cancer (n=1), synovial sarcoma (n=1), and benign neoplasms (n=4). The histological diagnoses of RAS included osteosarcoma (n=8), leiomyosarcoma (n=3), undifferentiated pleomorphic sarcoma (n=3), rhabdomyosarcoma (n=1), angiosarcoma (n=1), malignant peripheral nerve sheath tumor (n=1), spindle cell sarcoma NOS (n=1), and sarcoma not further specified (n=2). The median survival time from the diagnosis of RAS was 26 months. Eleven patients underwent surgery. Five of these patients achieved a continuous disease free status or showed no evidence disease. Four patients underwent CIRT. One of these patients with leiomyosarcoma achieved a continuous disease free status, and the other patient with osteosarcoma achieved a partial response. On the other hand, 2 patients experienced Grade 3 toxicities that required surgical treatment. Conclusion RAS originates from various types of diseases that are treated by radiotherapy and shows diverse pathological features. Complete resection achieves a good prognosis. CIRT can be an effective and feasible option for unresectable RAS.

Performance Evaluation for Repair of HSGc-C5 Carcinoma Cell Using Geant4-DNA

Track-structure Monte Carlo simulations are useful tools to evaluate initial DNA damage induced by irradiation. In the previous study, we have developed a Gean4-DNA-based application to estimate the cell surviving fraction of V79 cells after irradiation, bridging the gap between the initial DNA damage and the DNA rejoining kinetics by means of the two-lesion kinetics (TLK) model. However, since the DNA repair performance depends on cell line, the same model parameters cannot be used for different cell lines. Thus, we extended the Geant4-DNA application with a TLK model for the evaluation of DNA damage repair performance in HSGc-C5 carcinoma cells which are typically used for evaluating proton/carbon radiation treatment effects. For this evaluation, we also performed experimental measurements for cell surviving fractions and DNA rejoining kinetics of the HSGc-C5 cells irradiated by 70 MeV protons at the cyclotron facility at the National Institutes for Quantum and Radiological Science and Technology (QST). Concerning fast- and slow-DNA rejoining, the TLK model parameters were adequately optimized with the simulated initial DNA damage. The optimized DNA rejoining speeds were reasonably agreed with the experimental DNA rejoining speeds. Using the optimized TLK model, the Geant4-DNA simulation is now able to predict cell survival and DNA-rejoining kinetics for HSGc-C5 cells.

Comparison between the results of a recently-developed biological weighting function (V79-RBE10BWF) and thein vitroclonogenic survival RBE10of other repair-competent asynchronized normoxic mammalian cell lines and ions not used for the development of the model

728 simulated microdosimetric lineal energy spectra (26 different ions between 1H and 238U, 28 energy points from 1 to 1000 MeV/n) were used in combination with a recently-developed biological weighting function (Parisi et al., 2020) and 571 published in vitro clonogenic survival curves in order to: 1) assess prediction intervals for the in silico results by deriving an empirical indication of the experimental uncertainty from the dispersion in the in vitro hamster lung fibroblast (V79) data used for the development of the biophysical model; 2) explore the possibility of modeling the relative biological effectiveness (RBE) of the 10% clonogenic survival of asynchronized normoxic repair-competent mammalian cell lines other than the one used for the development of the model (V79); 3) investigate the predictive power of the model through a comparison between in silico results and in vitro data for 10 ions not used for the development of the model. At first, different strategies for the assessment of the in silico prediction intervals were compared. The possible sources of uncertainty responsible for the dispersion in the in vitro data were also shortly reviewed. Secondly, also because of the relevant scatter in the in vitro data, no statistically-relevant differences were found between the RBE10 of the investigated different asynchronized normoxic repair-competent mammalian cell lines. The only exception (Chinese Hamster peritoneal fibroblasts, B14FAF28), is likely due to the limited dataset (all in vitro ion data were extracted from a single publication), systematic differences in the linear energy transfer (LET) calculations for the employed very-heavy ions, and the use of reference photon survival curves extracted from a different publication. Finally, the in silico predictions for the 10 ions not used for the model development were in good agreement with the corresponding in vitro data.

DNA double-strand breaks in cancer cells as a function of proton linear energy transfer and its variation in time

PURPOSE

The complex relationship between linear energy transfer (LET) and cellular response to radiation is not yet fully elucidated. To better characterize DNA damage after irradiations with therapeutic protons, we monitored formation and disappearance of DNA double-strand breaks (DNA DSB) as a function of LET and time. Comparisons with conventional γ-rays and high LET carbon ions were also performed.

MATERIALS AND METHODS

In the present work, we performed immunofluorescence-based assay to determine the amount of DNA DSB induced by different LET values along the 62 MeV therapeutic proton Spread out Bragg peak (SOBP) in three cancer cell lines, i.e. HTB140 melanoma, MCF-7 breast adenocarcinoma and HTB177 non-small lung cancer cells. Time dependence of foci formation was followed as well. To determine irradiation positions, corresponding to the desired LET values, numerical simulations were carried out using Geant4 toolkit. We compared γ-H2AX foci persistence after irradiations with protons to that of γ-rays and carbon ions.

RESULTS

With the rise of LET values along the therapeutic proton SOBP, the increase of γ-H2AX foci number is detected in the three cell lines up to the distal end of the SOBP, while there is a decrease on its distal fall-off part. With the prolonged incubation time, the number of foci gradually drops tending to attain the residual level. For the maximum number of DNA DSB, irradiation with protons attain higher level than that of γ-rays. Carbon ions produce more DNA DSB than protons but not substantially. The number of residual foci produced by γ-rays is significantly lower than that of protons and particularly carbon ions. Carbon ions do not produce considerably higher number of foci than protons, as it could be expected due to their physical properties.

CONCLUSIONS

In situ visualization of γ-H2AX foci reveal creation of more lesions in the three cell lines by clinically relevant proton SOBP than γ-rays. The lack of significant differences in the number of γ-H2AX foci between the proton and carbon ion-irradiated samples suggests an increased complexity of DNA lesions and slower repair kinetics after carbon ions compared to protons. For all three irradiation types, there is no major difference between the three cell lines shortly after irradiations, while later on, the formation of residual foci starts to express the inherent nature of tested cells, therefore increasing discrepancy between them.

Definitive carbon ion radiotherapy for tracheobronchial adenoid cystic carcinoma: a preliminary report

Background

Tracheobronchial adenoid cystic carcinoma (TACC) is a rare tumour. About one-third of patients miss their chance of surgery or complete resection as it is mostly detected in the advanced stage; hence, photon radiotherapy (RT) is used. However, the outcomes of photon RT remain unsatisfactory. Carbon ion radiotherapy (CIRT) is thought to improve the therapeutic gain ratio; however, the outcomes of CIRT in TACC are unclear. Therefore, we aimed to assess the effects and toxicities of CIRT in patients with TACC.

Methods

The inclusion criteria were as follows: 1) age 18–80 years; 2) Eastern Cooperative Oncology Group Performance Status 0–2; 3) histologically confirmed TACC; 4) stage III–IV disease; 5) visible primary tumour; and 6) no previous RT history. The planned prescription doses of CIRT were 66–72.6 GyE/22–23 fractions. The rates of overall survival (OS), local control (LC), and progression-free survival (PFS) were calculated using the Kaplan-Meier method. Treatment-induced toxicities and tumour response were scored according to the Common Terminology Criteria for Adverse Events and Response Evaluation Criteria in Solid Tumors, respectively.

Results

Eighteen patients with a median age of 48 (range 30–73) years were enrolled. The median follow-up time was 20.7 (range 5.8–44.1) months. The overall response rate was 88.2%. Five patients developed lung metastasis after 12.2–41.0 months and one of them experienced local recurrence at 31.9 months after CIRT. The rates of 2-year OS, LC, and PFS were 100, 100, and 61.4%, respectively. Except for one patient who experienced grade 4 tracheal stenosis, which was relieved after stent implantation, no other ≥3 grade toxicities were observed.

Conclusions

CIRT might be safe and effective in the management of TACC based on a short observation period. Further studies with more cases and longer observation are warranted.

Improved Asymptotic Expansions in High- and Low-Dose Ranges for Generalized Multi-Hit Model of Radiation-Induced Cell Survival

By considering an upper bound on the number of radiation-induced potential lethal damages that can be repaired in a cell, we have proposed the generalized multi-hit (GMH) model with a closed-form solution, which can better fit various radiation-induced cell survival curves. Recent analysis shows that the asymptotic expansions that we gave before can be used to approximate the generalized single-hit single-target (GSHST) model rather than the GMH model. To illustrate the asymptotic trends of radiation-induced cell survival curves, in this study, we improve the asymptotic expansions of the GMH model in low- and high-dose ranges based on the limit formula of the incomplete gamma function in the corresponding dose ranges. When the upper limit of the number of radiation-induced potential lethal damages is one, the improved expansions of the GMH model can be reduced to the previous expansions of the GSHST model, and the improved asymptotic expansions of the GMH model also indicate that the GMH model has the generalized linear-quadratic-linear (LQL) feature. The numerical simulations indicate that the improved asymptotic expansions in high- and low-dose ranges agree well with the non-linear fitting of the GMH model in six kinds of cell lines under the corresponding dose ranges. In addition, we analyze the relative errors of the improved expansions of the GMH model in high- and low-dose ranges to demonstrate the accuracy and effectiveness of the improved expansions. Based on the error analysis, we further give the reasonable ranges of radiation dose applicable to the improved asymptotic expansions of the GMH model.

Cell lines of the same anatomic site and histologic type show large variability in intrinsic radiosensitivity and relative biological effectiveness to protons and carbon ions

PURPOSE

To show that intrinsic radiosensitivity varies greatly for protons and carbon (C) ions in addition to photons, and that DNA repair capacity remains important in governing this variability.

METHODS

We measured or obtained from the literature clonogenic survival data for a number of human cancer cell lines exposed to photons, protons (9.9 keV/μm), and C-ions (13.3-77.1 keV/μm). We characterized their intrinsic radiosensitivity by the dose for 10% or 50% survival (D_10% or D_50% ), and quantified the variability at each radiation quality by the coefficient of variation (COV) in D_10% and D_50% . We also treated cells with DNA repair inhibitors prior to irradiation to assess how DNA repair capacity affects their variability.

RESULTS

We found no statistically significant differences in the COVs of D_10% or D_50% between any of the radiation qualities investigated. The same was true regardless of whether the cells were treated with DNA repair inhibitors, or whether they were stratified into histologic subsets. Even within histologic subsets, we found remarkable differences in radiosensitivity for high LET C-ions that were often greater than the variations in RBE, with brain cancer cells varying in D_10% (D_50% ) up to 100% (131%) for 77.1 keV/μm C-ions, and non-small cell lung cancer and pancreatic cancer cell lines varying up to 55% (76%) and 51% (78%), respectively, for 60.5 keV/μm C-ions. The cell lines with modulated DNA repair capacity had greater variability in intrinsic radiosensitivity across all radiation qualities.

CONCLUSIONS

Even for cell lines of the same histologic type, there are remarkable variations in intrinsic radiosensitivity, and these variations do not differ significantly between photon, proton or C-ion radiation. The importance of DNA repair capacity in governing the variability in intrinsic radiosensitivity is not significantly diminished for higher LET radiation.

A systematic review on the usage of averaged LET in radiation biology for particle therapy

Linear Energy Transfer (LET) is widely used to express the radiation quality of ion beams, when characterizing the biological effectiveness. However, averaged LET may be defined in multiple ways, and the chosen definition may impact the resulting reported value. We review averaged LET definitions found in the literature, and quantify which impact using these various definitions have for different reference setups. We recorded the averaged LET definitions used in 354 publications quantifying the relative biological effectiveness (RBE) of hadronic beams, and investigated how these various definitions impact the reported averaged LET using a Monte Carlo particle transport code. We find that the kind of averaged LET being applied is, generally, poorly defined. Some definitions of averaged LET may influence the reported averaged LET values up to an order of magnitude. For publications involving protons, most applied dose averaged LET when reporting RBE. The absence of what target medium is used and what secondary particles are included further contributes to an ill-defined averaged LET. We also found evidence of inconsistent usage of averaged LET definitions when deriving LET-based RBE models. To conclude, due to commonly ill-defined averaged LET and to the inherent problems of LET-based RBE models, averaged LET may only be used as a coarse indicator of radiation quality. We propose a more rigorous way of reporting LET values, and suggest that ideally the entire particle fluence spectra should be recorded and provided for future RBE studies, from which any type of averaged LET (or other quantities) may be inferred.

A Mathematical Radiobiological Model (MRM) to Predict Complex DNA Damage and Cell Survival for Ionizing Particle Radiations of Varying Quality

Predicting radiobiological effects is important in different areas of basic or clinical applications using ionizing radiation (IR); for example, towards optimizing radiation protection or radiation therapy protocols. In this case, we utilized as a basis the ‘MultiScale Approach (MSA)’ model and developed an integrated mathematical radiobiological model (MRM) with several modifications and improvements. Based on this new adaptation of the MSA model, we have predicted cell-specific levels of initial complex DNA damage and cell survival for irradiation with ¹¹Β, ¹²C, ¹⁴Ν, ¹⁶Ο, ²⁰Νe, ⁴⁰Αr, ²⁸Si and ⁵⁶Fe ions by using only three input parameters (particle’s LET and two cell-specific parameters: the cross sectional area of each cell nucleus and its genome size). The model-predicted survival curves are in good agreement with the experimental ones. The particle Relative Biological Effectiveness (RBE) and Oxygen Enhancement Ratio (OER) are also calculated in a very satisfactory way. The proposed integrated MRM model (within current limitations) can be a useful tool for the assessment of radiation biological damage for ions used in hadron-beam radiation therapy or radiation protection purposes.

High Doses of Photons and Carbon Ions Comparably Increase Vascular Permeability in R3327-HI Prostate Tumors: A Dynamic Contrast-Enhanced MRI Study

Carbon- (12C-) ion radiotherapy exhibits enhanced biological effectiveness compared to photon radiotherapy, however, the contribution of its interaction with the vasculature remains debatable. The effect of high-dose 12C-ion and photon irradiation on vascular permeability in moderately differentiated rat prostate tumors was compared using dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI). Syngeneic R3327-HI rat prostate tumors were irradiated with a single dose of either 18 or 37 Gy 12C ions, or 37 or 75 Gy 6-MV photons (sub-curative and curative dose levels, respectively). DCE-MRI was performed one day prior to and 3, 7, 14 and 21 days postirradiation. Voxel-based tumor concentration-time curves were clustered based on their curve shape and treatment response was assessed as the longitudinal changes in the relative abundance per cluster. Radiation-induced vascular damage and increased permeability occurred at day 7 postirradiation for all treatment groups except for the 75 Gy photon-irradiated group, where the onset of vascular damage was delayed until day 14. No differences between irradiation modalities were found. Therefore, early vascular damage cannot explain the higher effectiveness of 12C ions relative to photons in terms of local tumor control for this moderately differentiated prostate tumor and the applied single high doses.

The relative biological effectiveness of carbon ion radiation therapy for early stage lung cancer

BACKGROUND AND PURPOSE

Carbon ion radiation therapy (CIRT) is recognized as an effective alternative treatment modality for early stage lung cancer, but a quantitative understanding of relative biological effectiveness (RBE) compared to photon therapy is lacking. In this work, a mechanistic tumor response model previously validated for lung photon radiotherapy was used to estimate the RBE of CIRT compared to photon radiotherapy, as a function of dose and the fractionation schedule.

MATERIALS AND METHODS

Clinical outcome data of 9 patient cohorts (394 patients) treated with CIRT for early stage lung cancer, representing all published data, were included. Fractional dose, number of fractions, treatment schedule, and local control rates were used for model simulations relative to standard photon outcomes. Four parameters were fitted: α, α/β, and the oxygen enhancement ratios of cells either accessing only glucose, not oxygen (OER_I), or cells dying from starvation (OER_H). The resulting dose-response relationship of CIRT was compared with the previously determined dose-response relationship of photon radiotherapy for lung cancer, and an RBE of CIRT was derived.

RESULTS

Best-fit CIRT parameters were: α = 1.12 Gy^-1 [95%-CI: 0.97-1.26], α/β = 23.9 Gy [95%-CI: 8.9-38.9], and the oxygen induced radioresistance of hypoxic cell populations were characterized by OER_I = 1.08 [95%-CI: 1.00-1.41] (cells lacking oxygen but not glucose), and OER_H = 1.01 [95%-CI: 1.00-1.44] (cells lacking oxygen and glucose). Depending on dose and fractionation, the derived RBE ranges from 2.1 to 1.5, with decreasing values for larger fractional dose and fewer number of fractions.

CONCLUSION

Fitted radiobiological parameters were consistent with known carbon in vitro radiobiology, and the resulting dose-response curve well-fitted the reported data over a wide range of dose-fractionation schemes. The same model, with only a few fitted parameters of clear mechanistic meaning, thus synthesizes both photon radiotherapy and CIRT clinical experience with early stage lung tumors.

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.

High LET Radiation Overcomes In Vitro Resistance to X-Rays of Chondrosarcoma Cell Lines

Chondrosarcomas are malignant tumors of the cartilage that are chemoresistant and radioresistant to X-rays. This restricts the treatment options essential to surgery. In this study, we investigated the sensitivity of chondrosarcoma to X-rays and C-ions in vitro. The sensitivity of 4 chondrosarcoma cell lines (SW1353, CH2879, OUMS27, and L835) was determined by clonogenic survival assays and cell cycle progression. In addition, biomarkers of DNA damage responses were analyzed in the SW1353 cell line. Chondrosarcoma cells showed a heterogeneous sensitivity toward irradiation. Chondrosarcoma cell lines were more sensitive to C-ions exposure compared to X-rays. Using D10 values, the relative biological effectiveness of C-ions was higher (relative biological effectiveness = 5.5) with cells resistant to X-rays (CH2879) and lower (relative biological effectiveness = 3.7) with sensitive cells (L835). C-ions induced more G2 phase blockage and micronuclei in SW1353 cells as compared to X-rays with the same doses. Persistent unrepaired DNA damage was also higher following C-ions irradiation. These results indicate that chondrosarcoma cell lines displayed a heterogeneous response to conventional radiation treatment; however, treatment with C-ions irradiation was more efficient in killing chondrosarcoma cells, compared to X-rays.

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.

Significance of concurrent use of weekly cisplatin in carbon-ion radiotherapy for locally advanced adenocarcinoma of the uterine cervix: A propensity score-matched analysis

BACKGROUND

Although carbon-ion radiotherapy (C-ion RT) with concurrent chemotherapy (chemo-C-ion RT) is a promising treatment for adenocarcinoma (AC) of the uterine cervix, its long-term efficacy remains unclear. We evaluated the long-term significance of concurrent weekly cisplatin and C-ion RT for locally advanced AC of the uterine cervix.

METHODS

We performed a pooled analysis of patients with stage IIB-IVA AC of the uterine cervix who underwent C-ion RT alone or chemo-C-ion RT between September 2007 and December 2018 at our institution. Patients received 74.4 Gy (relative biological effectiveness) with or without cisplatin (40 mg/m² per week for up to 5 weeks), underwent no prior pelvic RT or systemic therapy, and had a performance status of 0-2. Propensity score matching was based on the year of diagnosis, regional lymph node metastasis, and stage.

RESULTS

The matched cohort contained 26 patients who underwent C-ion RT and 26 who underwent chemo-C-ion RT. The median age and follow-up period were 57 (range, 28-79) years and 34 (range, 2-126) months, respectively. The 5-year overall survival rate was significantly better in the chemo-C-ion RT group (72%) than in the C-ion RT group (46%; P = .041). The 5-year distant metastatic-free rate was also significantly better in the chemo-C-ion RT group (66%) than in the C-ion RT group (41%; P = .048). The incidence of grade ≥ 3 late toxicities was comparable between the two groups.

CONCLUSIONS

Chemo-C-ion RT for locally advanced AC of the uterine cervix is associated with a long-term survival benefit.

Rare entities in head-and-neck cancer: salvage re-irradiation with carbon ions

Background

The objective of this investigation is to evaluate the outcomes and toxicity of carbon-ion re-irradiation (CIR) in patients with rare head and neck cancers (HNC). There is a paucity of data regarding treatment approaches in this patient cohort, which we aim to address in this work.

Methods

Thirty-two (n = 32) consecutive patients with uncommon HNC treated between 2010 and 2017 were retrospectively analyzed in terms of clinical outcomes, patterns of failure, and toxicity.

Results

Mucoepidermoid carcinoma (MEC) was the most common histology (22%). Patients received a median cumulative dose equivalent in 2 Gy fractions (EQD₂) after CIR of 128.6 Gy (range, 105.8–146.5 Gy). The local and distant control rates 1 year after CIR were 66 and 72%. No serious acute or late toxicity (≥ grade 3) after CIR was observed.

Conclusions

CIR may represent an effective and safe treatment alternative to palliative systemic therapies in these rare indications.

RBE-weighted dose conversions for carbon ionradiotherapy between microdosimetric kinetic model and local effect model for the targets and organs at risk in prostate carcinoma

BACKGROUND AND PURPOSE

The aim of this study was to establish curves for the conversion of RBE-weighted doses for targets and organs at risk (OARs) from the microdosimetric kinetic model (MKM) calculation to that of the local effect model I (LEM) for carbon ion radiotherapy (CIRT) for prostate carcinoma (PCA).

MATERIALS AND METHODS

This study was performed in the experimental treatment planning system (eTPS, V8A, Raystation, Sweden), which incorporates both MKM and LEM. CIRT plans from 10 PCA patients were collected. There were 5 steps to establish the curves: (1) design MKM plans in eTPS; (2) recalculate the physical doses from MKM to LEM and create a LEM plan in eTPS; (3) plot the RBE-weighted MKM to LEM conversion curves; (4) convert the MKM rectum constraint dose volume histogram (DVH) from NIRS to a LEM DVH; and (5) compare patients' rectum DVHs and follow-up with the converted constraint DVH.

RESULTS

The conversion factors for MKM doses of 0.18 Gy (RBE) to 4.55 Gy (RBE) per fraction to LEM doses were 2.72-1.06. For fraction sizes of >1 Gy (RBE), the conversion factors matched Fossati's curve and for fraction sizes of <1.00 Gy (RBE) the values were on the extrapolated Fossati's curve. A LEM rectum constraint DVH was established. Ten patients' rectum DVHs were all lower than LEM constraint DVHs. No complications were reported clinically.

CONCLUSION

For PCA receiving CIRT, the RBE-weighted doses using MKM for targets and OARs could be converted to LEM doses using conversion curves.

Distinct modes of death in human neural stem and glioblastoma cells irradiated with carbon-ion radiation and gamma-rays

Purpose: Accumulated damage in neural stem cells (NSCs) during brain tumor radiotherapy causes cognitive dysfunction to the patients. Carbon-ion radiotherapy can reduce undesired irradiation of normal tissues more efficiently than conventional photon radiotherapy. This study elucidates the responses of NSCs to carbon-ion radiation.Methods: Human NSCs and glioblastoma A-172 cells were irradiated with carbon-ion radiation and γ-rays, which have different linear-energy-transfer (LET) values of 108 and 0.2 keV/μm, respectively. After irradiation, growth rates were measured, apoptotic cells were detected by flow cytometry, and DNA synthesizing cells were immunocytochemically visualized.Results: Growth rates of NSCs and A-172 cells were decreased after irradiation. The percentages of apoptotic cells were remarkably increased in NSCs but not in A-172 cells. In contrast, the fractions of DNA synthesizing A-172 cells were decreased in a dose-dependent manner. These results indicate that apoptosis induction and DNA synthesis inhibition contribute to the growth inhibition of NSCs and glioblastoma cells, respectively. In addition, high-LET carbon ions induced more profound effects than low-LET γ-rays.Conclusions: Apoptosis is an important clinical target to protect NSCs during brain tumor radiotherapy using carbon-ion radiation as well as conventional X-rays.

Targeting the Hedgehog pathway in combination with X‑ray or carbon ion radiation decreases migration of MCF‑7 breast cancer cells

The use of carbon ion therapy for cancer treatment is becoming more widespread due to the advantages of carbon ions compared with X‑rays. Breast cancer patients may benefit from these advantages, as the surrounding healthy tissues receive a lower dose, and the increased biological effectiveness of carbon ions can better control radioresistant cancer cells. Accumulating evidence indicates that the Hedgehog (Hh) pathway is linked to the development and progression of breast cancer, as well as to resistance to X‑irradiation and the migratory capacity of cancer cells. Hence, there is an increasing interest in targeting the Hh pathway in combination with radiotherapy. Several studies have already investigated this treatment strategy with conventional radiotherapy. However, to the best of our knowledge, the combination of Hh inhibitors with particle therapy has not yet been explored. The aim of the present study was to investigate the potential of the Hh inhibitor GANT61 as an effective modulator of radiosensitivity and migration potential in MCF‑7 breast cancer cells, and compare potential differences between carbon ion irradiation and X‑ray exposure. Although Hh targeting was not able to radiosensitise cells to any radiation type used, the combination of GANT61 with X‑rays or carbon ions (energy: 95 MeV/n; linear energy transfer: 73 keV/µm) was more effective in decreasing MCF‑7 cell migration compared with either radiation type alone. Gene expression of the Hh pathway was affected to different degrees in response to X‑ray and carbon ion irradiation, as well as in response to the combination of GANT61 with irradiation. In conclusion, combining Hh inhibition with radiation (X‑rays or carbon ions) more effectively decreased breast cancer cell migration compared with radiation treatment alone.

Sensitization of chondrosarcoma cells with PARP inhibitor and high-LET radiation

Chondrosarcoma is a malignant tumor that arises from cartilaginous tissue and is radioresistant and chemoresistant to conventional treatments. The preferred treatment consists of surgical resection, which might cause severe disabilities for the patient; in addition, this procedure might be impossible for inoperable locations, such as the skull base. Carbon ion irradiation (hadron therapy) has been proposed as an alternative treatment, primarily due to its greater biological effectiveness and improved ballistic properties compared with conventional radiotherapy with X-rays. The goal of this study was to characterize the genetic mutations of a grade III chondrosarcoma cell line (CH2879) and examine the cellular responses to conventional radiotherapy (X-rays) and hadron therapy (proton and carbon ions) in the presence of the PARP inhibitor Olaparib. To better understand PARP inhibition, we first analyzed the formation of poly-ADP ribose chains by western blot; we observed an increase in its signal after irradiation, which disappeared on addition of the PARP inhibitor. PARPi enhanced ratio of approximately 1.3, 1.8, and 1.5 following irradiation of cells with X-rays, protons, and C-ions, respectively, as detected by clonogenic assay. The decrease in cell survival was confirmed by proliferation assay. The radiosensitivity of CH2879 cells was associated with mutations in homologous recombination repair genes, such as RAD50, SMARCA2 and NBN. This study demonstrates the capacity of the PARP inhibitor Olaparib to radiosensitize mutated chondrosarcoma cells to conventional photon irradiation, proton and carbon ion irradiation.

Predicting In Vitro Cancer Cell Survival Based on Measurable Cell Characteristics

Variation in cellular characteristics may determine tumor response and, consequently, patient survival in radiation therapy. However, patient-specific prediction of cellular radiation response is currently unavailable for treatment planning. Thus, the importance of developing a novel approach based on clinically accessible parameters prior to treatment (e.g., by biopsy) is high. The goal of this study was to predict in vitro cancer cell survival through the p53mutation status and the number of chromosomes (NoC). To predict cell survival, we modified a mechanistic radiation response model incorporating DNA repair and cell death, originally designed for normal human cells. Cell-specific parameters of 24 cell lines originating from two laboratories (OncoRay, Dresden, Germany and HIMAC, Chiba, Japan) were considered for modeling. In a first step, we obtained estimates of the only unknown model input parameter genome size (GS) by fitting cell survival simulations onto experimental data. We then analyzed measured and published input model parameters (NoC, p53-mutation status and cell-cycle distribution) to assess their impact on measured and simulated parameters (modeled GS, and measured α, β, SF₂ and γ-H2AX). The resulting data suggested a linear correlation between NoC and modeled GS (R² > 0.93) allowing for estimating GS based on NoC. Applying the estimated GS resulted in predicted cell survival that matched measured data mostly within the experimental uncertainty. The measured radiobiological value β increased quadratically with the cell's modeled GS irrespective of other cell-specific parameters. The measured α and SF₂ split into two groups, depending on the cells' p53-mutation status, both linearly increasing and decreasing, respectively, with modeled GS. Model predictions of foci numbers were, on average, in agreement with published γ-H2AX measurement data. In conclusion, knowledge of clinically accessible parameters (p53-mutation status and NoC) may support patient stratification in radiotherapy based on cell-specific survival prediction testable in prospective clinical trials.

Bystander effectors of chondrosarcoma cells irradiated at different LET impair proliferation of chondrocytes

While the dose-response relationship of radiation-induced bystander effect (RIBE) is controversial at low and high linear energy transfer (LET), mechanisms and effectors of cell-to-cell communication stay unclear and highly dependent of cell type. In the present study, we investigated the capacity of chondrocytes in responding to bystander factors released by chondrosarcoma cells irradiated at different doses (0.05 to 8 Gy) with X-rays and C-ions. Following a medium transfer protocol, cell survival, proliferation and DNA damages were quantified in bystander chondrocytes. The bystander factors secreted by chondrosarcoma cells were characterized. A significant and major RIBE response was observed in chondrocyte cells (T/C-28a2) receiving conditioned medium from chondrosarcoma cells (SW1353) irradiated with 0.1 Gy of X-rays and 0.05 Gy of C-ions, resulting in cell survivals of 36% and 62%, respectively. Micronuclei induction in bystander cells was observed from the same low doses. The cell survival results obtained by clonogenic assays were confirmed using impedancemetry. The bystander activity was vanished after a heat treatment or a dilution of the conditioned media. The cytokines which are well known as bystander factors, TNF-α and IL-6, were increased as a function of doses and LET according to an ELISA multiplex analysis. Together, the results demonstrate that irradiated chondrosarcoma cells can communicate stress factors to non-irradiated chondrocytes, inducing a wide and specific bystander response related to both doses and LET.

The radiobiological effects of He, C and Ne ions as a function of LET on various glioblastoma cell lines

The effects of the charged ion species 4He, 12C and 20Ne on glioblastoma multiforme (GBM) T98G, U87 and LN18 cell lines were compared with the effects of 200 kVp X-rays (1.7 keV/μm). These cell lines have different genetic profiles. Individual GBM relative biological effectiveness (RBE) was estimated in two ways: the RBE10 at 10% survival fraction and the RBE2Gy after 2 Gy doses. The linear quadratic model radiosensitivity parameters α and β and the α/β ratio of each ion type were determined as a function of LET. Mono-energetic 4He, 12C and 20Ne ions were generated by the Heavy Ion Medical Accelerator at the National Institute of Radiological Sciences in Chiba, Japan. Colony-formation assays were used to evaluate the survival fractions. The LET of the various ions used ranged from 2.3 to 100 keV/μm (covering the depth-dose plateau region to clinically relevant LET at the Bragg peak). For U87 and LN18, the RBE10 increased with LET and peaked at 85 keV/μm, whereas T98G peaked at 100 keV/μm. All three GBM α parameters peaked at 100 keV/μm. There is a statistically significant difference between the three GBM RBE10 values, except at 100 keV/μm (P < 0.01), and a statistically significant difference between the α values of the GBM cell lines, except at 85 and 100 keV/μm. The biological response varied depending on the GBM cell lines and on the ions used.

Radiobiological Characterization of Canine Malignant Melanoma Cell Lines with Different Types of Ionizing Radiation and Efficacy Evaluation with Cytotoxic Agents

Canine malignant melanoma (CMM) is a locally and systemically aggressive cancer that shares many biological and clinical characteristics with human mucosal melanoma. Hypofractionated radiation protocols have been used to treat CMM but little is known about its radiation biology. This pilot study is designed to investigate response of CMM cell lines to various ionizing radiations and cytotoxic agents to better understand this canine cancer. Four CMM cell lines were evaluated by clonogenic survival assay under aerobic and hypoxic conditions and parameters such as alpha beta (α/β) ratio, oxygen enhancement ratio (OER), and relative biological effectiveness (RBE) were calculated after ¹³⁷Cs, 6 megavoltage (MV) photon, or carbon ion irradiation. Six cytotoxic agents (cisplatin, camptothecin, mitomycin C, bleomycin, methtyl methanesulfonate and etoposide) were also assessed for their efficacy. Under aerobic condition with 6 MV photon, the α/β ratio of the four cell lines ranged from 0.3 to >100, indicating a wide variation of cellular sensitivity. The ratio increased under hypoxic condition compared to aerobic condition and this was more dramatic in ¹³⁷Cs and 6 MV photon treatments. OER of carbon was lower than ¹³⁷Cs at D₁₀ in 3 of the 4 cell lines. The RBE values generally increased with the increase of LET. Different cell lines showed sensitivity/resistance to different cytotoxic agents. This study revealed that CMM has a wide range of radiosensitivity and that hypoxia can reduce it, indicating that widely used hypofractionated protocols may not be optimal for all CMM patients. Several cytotoxic agents that have never been clinically assessed can improve treatment outcome.

Potential lethal damage repair in glioblastoma cells irradiated with ion beams of various types and levels of linear energy transfer

Glioblastoma (GBM), a Grade IV brain tumour, is a well-known radioresistant cancer. To investigate one of the causes of radioresistance, we studied the capacity for potential lethal damage repair (PLDR) of three altered strains of GBM: T98G, U87 and LN18, irradiated with various ions and various levels of linear energy transfer (LET). The GBM cells were exposed to 12C and 28Si ion beams with LETs of 55, 100 and 200 keV/μm, and with X-ray beams of 1.7 keV/μm. Mono-energetic 12C ions and 28Si ions were generated by the Heavy Ion Medical Accelerator at the National Institute of Radiological Science, Chiba, Japan. Clonogenic assays were used to determine cell inactivation. The ability of the cells to repair potential lethal damage was demonstrated by allowing one identical set of irradiated cells to repair for 24 h before subplating. The results show there is definite PLDR with X-rays, some evidence of PLDR at 55 keV/μm, and minimal PLDR at 100 keV/μm. There is no observable PLDR at 200 keV/μm. This is the first study, to the authors' knowledge, demonstrating the capability of GBM cells to repair potential lethal damage following charged ion irradiations. It is concluded that a GBM's PLDR is dependent on LET, dose and GBM strain; and the more radioresistant the cell strain, the greater the PLDR.

Single fraction carbon ion radiotherapy for colorectal cancer liver metastasis: A dose escalation study

Prognosis is usually grim for those with liver metastasis from colorectal cancer (CRC) who cannot receive resection. Radiation therapy can be an option for those unsuitable for resection, with carbon ion radiotherapy (CIRT) being more effective and less toxic than X-ray due to its physio-biological characteristics. The objective of this study is to identify the optimal dose of single fraction CIRT for colorectal cancer liver metastasis. Thirty-one patients with liver metastasis from CRC were enrolled in the present study. Twenty-nine patients received a single-fraction CIRT, escalating the dose from 36 Gy (RBE) in 5% to 10% increments until unacceptable incidence of dose-limiting toxicity was observed. Dose-limiting toxicity was defined as grade ≥3 acute toxicity attributed to radiotherapy. The prescribed doses were as follows: 36 Gy (RBE) (3 cases), 40 Gy (2 cases), 44 Gy (4 cases), 46 Gy (6 cases), 48 Gy (3 cases), 53 Gy (8 cases) and 58 Gy (3 cases). Dose-limiting toxicity was not observed, but late grade 3 liver toxicity due to biliary obstruction was observed in 2 patients at 53 Gy (RBE). Both cases had lesions close to the hepatic portal region, and, therefore, the dose was escalated to 58 Gy (RBE), limited to peripheral lesions. The 3-year actuarial overall survival rate of all 29 patients was 78%, and the median survival time was 65 months. Local control improved significantly at ≥53 Gy (RBE), with a 3-year actuarial local control rate of 82%, compared to 28% in lower doses. Treatment for CRC liver metastasis with single-fraction CIRT appeared to be safe up to 58 Gy (RBE) as long as the central hepatic portal region was avoided.