Treatment planning for heavy-ion radiotherapy: calculation and optimization of biologically effective dose

NRG Oncology White Paper on the Relative Biological Effectiveness in Proton Therapy

This position paper, led by the NRG Oncology Particle Therapy Work Group, focuses on the concept of relative biologic effect (RBE) in clinical proton therapy, with the goal of providing recommendation for the next generation clinical trials with proton therapy on the best practice of investigating and using RBE, which could deviate from the current standard proton RBE value of 1.1 relative to photons. In part one, current clinical utilization and practice are reviewed, giving the context and history of RBE. Evidence for variation in RBE is presented along with the concept of linear energy transfer (LET). The intertwined nature of tumor radiobiology, normal tissue constraints, and treatment planning with LET and RBE considerations is then reviewed. Part two summarizes current and past clinical data and then suggests the next steps to explore and employ tools for improved dynamic models for RBE. In part three, approaches and methods for the next generation of prospective clinical trials are explored, with the goal of optimizing RBE to be both more reflective of clinical reality and also deployable in trials to allow clinical validation and inter-patient comparisons. These concepts provide the foundation for personalized biologic treatments reviewed in part four. Finally, we conclude with a summary including short- and long-term scientific focus points for clinical proton therapy. The practicalities and capacity to use RBE in treatment planning are reviewed and considered with more biological data in hand. The intermediate step of LET optimization is summarized and proposed as a potential bridge to the ultimate goal of case-specific RBE planning that can be achieved as a hypothesis generating tool in near-term proton trials.

Up modulation of dose-averaged linear energy transfer by simultaneous integrated boost in carbon-ion radiotherapy for pancreatic carcinoma

BACKGROUND

Local recurrence in locally advanced pancreatic cancer (LAPC) after carbon-ion radiotherapy (CIRT) may partly attribute to low dose-averaged linear energy transfer (LETd), despite high CIRT dose.

PURPOSE

This study aimed to investigate the approaches to up-modulate the CIRT LETd and to evaluate the corresponding oxygen enhancement ratio (OER) reduction.

METHODS

10 LAPCs that had been irradiated by CIRT with 67.5 Gy (RBE) in 15 fractions were selected. Their original plans were taken as the control plan for the LETd and OER investigations. Our considerations for up-modulating LETd were: (1) to deliver high doses to gross tumor volume core (GTVcore), while keeping dose constraints of the gastrointestinal (GI) tract in tolerance; (2) to put more Bragg-peak (BP) within the modulated targets; (3) to increase the BP density, high doses were necessary; (4) CIRT LETd could be effectively increased to small volumes; and (5) simultaneous integrated boost technique (SIB) could achieve the aforementioned tasks. The LETd and the corresponding OER distributions of each type of SIB plan were evaluated.

RESULTS

We delivered up to 100 Gy (RBE) to GTVcore using SIB. The mean LETd of GTV increased significantly by 21.3% from 47.8 to 58.0 keV/μm (p < 0.05). Meanwhile, the mean OER of GTVcore decreased by 6.6%, from 1.51 to 1.41 (p < 0.05). The GI LETdS in all modulated plans were not more than those in the original plans.

CONCLUSIONS

SIB could effectively increase CIRT LETd to LAPC, thus producing reduced OER, which may effectively overcome the radioresistance of LAPCs.

Biological dose optimization incorporating intra-tumoural cellular radiosensitivity heterogeneity in ion-beam therapy treatment planning

Objective.Treatment plans of ion-beam therapy have been made under an assumption that all cancer cells within a tumour equally respond to a given radiation dose. However, an intra-tumoural cellular radiosensitivity heterogeneity clearly exists, and it may lead to an overestimation of therapeutic effects of the radiation. The purpose of this study is to develop a biological model that can incorporate the radiosensitivity heterogeneity into biological optimization for ion-beam therapy treatment planning.Approach.The radiosensitivity heterogeneity was modeled as the variability of a cell-line specific parameter in the microdosimetric kinetic model following the gamma distribution. To validate the developed intra-tumoural-radiosensitivity-heterogeneity-incorporated microdosimetric kinetic (HMK) model, a treatment plan with H-ion beams was made for a chordoma case, assuming a radiosensitivity heterogeneous region within the tumour. To investigate the effects of the radiosensitivity heterogeneity on the biological effectiveness of H-, He-, C-, O-, and Ne-ion beams, the relative biological effectiveness (RBE)-weighted dose distributions were planned for a cuboid target with the stated ion beams without considering the heterogeneity. The planned dose distributions were then recalculated by taking the heterogeneity into account.Main results. The cell survival fraction and corresponding RBE-weighted dose were formulated based on the HMK model. The first derivative of the RBE-weighted dose distribution was also derived, which is needed for fast biological optimization. For the patient plan, the biological optimization increased the dose to the radiosensitivity heterogeneous region to compensate for the heterogeneity-induced reduction in biological effectiveness of the H-ion beams. The reduction in biological effectiveness due to the heterogeneity was pronounced for low linear energy transfer (LET) beams but moderate for high-LET beams. The RBE-weighted dose in the cuboid target decreased by 7.6% for the H-ion beam, while it decreased by just 1.4% for the Ne-ion beam.Significance.Optimal treatment plans that consider intra-tumoural cellular radiosensitivity heterogeneity can be devised using the HMK model.

Modeling of the resensitization effect on carbon-ion radiotherapy for stage I non-small cell lung cancer

Objective. To investigate the effect of redistribution and reoxygenation on the 3-year tumor control probability (TCP) of patients with stage I non-small cell lung cancer (NSCLC) treated with carbon-ion radiotherapy.Approach. A meta-analysis of published clinical data of 233 NSCLC patients treated by carbon-ion radiotherapy under 18-, 9-, 4-, and single-fraction schedules was conducted. The linear-quadratic (LQ)-based cell-survival model incorporating the radiobiological 5Rs, radiosensitivity, repopulation, repair, redistribution, and reoxygenation, was developed to reproduce the clinical TCP data. Redistribution and reoxygenation were regarded together as a single phenomenon and termed 'resensitization' in the model. The optimum interval time between fractions was investigated for each fraction schedule using the determined model parameters.Main results.The clinical TCP data for 18-, 9-, and 4-fraction schedules were reasonably reproduced by the model without the resensitization effect, whereas its incorporation was essential to reproduce the TCP data for all fraction schedules including the single fraction. The curative dose for the single-fraction schedule was estimated to be 49.0 Gy (RBE), which corresponds to the clinically adopted dose prescription of 50.0 Gy (RBE). For 18-, 9-, and 4-fraction schedules, a 2-to-3-day interval is required to maximize the resensitization effect during the time interval. In contrast, the single-fraction schedule cannot benefit from the resensitization effect, and the shorter treatment time is preferable to reduce the effect of sub-lethal damage repair during the treatment.Significance.The LQ-based cell-survival model incorporating the radiobiological 5Rs was developed and used to evaluate the effect of the resensitization on clinical results of NSCLC patients treated with hypo-fractionated carbon-ion radiotherapy. The incorporation of the resensitization into the cell-survival model improves the reproducibility to the clinical TCP data. A shorter treatment time is preferable in the single-fraction schedule, while a 2-to-3-day interval between fractions is preferable in the multi-fraction schedules for effective treatments.

A digital male pelvis phantom series showing anatomical variations over the course of fractionated radiotherapy treatment

BACKGROUND

Daily IGRT images show day-to-day anatomical variations in patients undergoing fractionated prostate radiotherapy. This is of particular importance in particle beam treatments.

PURPOSE

To develop a digital phantom series showing variation in pelvic anatomy for evaluating treatment planning and IGRT procedures in particle radiotherapy.

METHODS

A pelvic phantom series was developed from the planning MRI and kVCT (planning CT) images along with six of the daily serial MVCT images taken of a single patient treated with a full bladder on a Tomotherapy unit. The selected patient had clearly visible yet unexceptional internal anatomy variation. Prostate, urethra, bladder, rectum, bowel, bowel gas, bone and soft tissue were contoured and a single Hounsfield Unit was assigned to each region. Treatment plans developed on the kVCT for photon, proton and carbon beams were recalculated on each phantom to demonstrate a clinical application of the series. Proton plans were developed with and without robust optimization.

RESULTS

Limited to axial slices with prostate, the bladder volume varied from 6 to 46 cm3, the rectal volume (excluding gas) from 22 to 52 cm3, and rectal gas volume from zero to 18 cm3. The water equivalent path length to the prostate varied by up to 1.5 cm . The variations resulted in larger changes in the RBE-weighted Dose Volume Histograms of the non-robust proton plan and the carbon plan compared to the robust proton plan, the latter similar to the photon plan. The prostate coverage (V100%) decreased by an average of 18% in the carbon plan, 16% in the non-robust proton plan, 1.8% in the robust proton plan, and 4.4% in the photon plan. The volume of rectum receiving 75% of the prescription dose (V75%) increased by an average of 3.7 cm3, 4.7 cm3, 1.9 cm3, and 0.6 cm3 in those four plans, respectively.

CONCLUSIONS

The digital pelvic phantom series provides for quantitative investigation of IGRT procedures and new methods for improving accuracy in particle therapy and may be used in cross-institutional comparisons for clinical trial quality assurance.

The rationale for a carbon ion radiation therapy facility in Australia

Australia has taken a collaborative nationally networked approach to achieve particle therapy capability. This supports the under-construction proton therapy facility in Adelaide, other potential proton centres and an under-evaluation proposal for a hybrid carbon ion and proton centre in western Sydney. A wide-ranging overview is presented of the rationale for carbon ion radiation therapy, applying observations to the case for an Australian facility and to the clinical and research potential from such a national centre.

Comment on "Reproducibility study of Monte Carlo simulations for nanoparticle dose enhancement and biological modeling of cell survival curves" by Veltenet al[Biomed Phys Eng Express 2023;9:045004]

This comment highlights two methodological issues with the recent article by Velten et al [Biomed Phys Eng Express 2023;9:045004]: First, the approach taken in this work with a local effect model (LEM) in 2D leads to a significant overstimation of the number of radiation-induced lesions. This results in order of magnitude smaller predicted survival rates compared to the conventional LEM. Second, the dose without nanoparticles is used as the 'macroscopic dose' against which cell survival is plotted. However, for the considered gold concentrations, the average absorbed dose under secondary particle equilibrium is between 2 and 20 times higher with nanoparticles than without.

Intensity-modulated proton and carbon-ion radiotherapy using a fixed-beam system for locally advanced lung cancer: dosimetric comparison with x-ray radiotherapy and normal tissue complication probability (NTCP) evaluation

Objective. To assess the dosimetric consequences and the normal tissue complication probability (NTCP) for the organs at risk (OARs) in intensity-modulated particle radiotherapy of proton (IMPT) and carbon-ion (IMCT) using a fixed-beam delivery system when compared with intensity-modulated photon radiotherapy (IMRT) for locally advanced small-cell lung cancer.Approach. The plans were all designed under the same total relative biological effectiveness (RBE)-weighted prescription dose, in which the planning target volume (PTV) of the internal gross target volume(IGTV) and the PTV of the clinical target volume was irradiated with 69.3 Gy (RBE) and 63 Gy (RBE), respectively, using a simultaneously integrated boosting (SIB) technique. NTCPs were estimated for heart, lung, esophagus and spinal cord by Lyman-Kutcher-Burman (LKB) and logistic models. Dose escalation was simulated under the desired NTCP values (0.05, 0.10 and 0.50) of the three radiation techniques.Main results. Under the similar target coverage, almost all OARs were significantly better spared (p< 0.05) when using the particle radiotherapy except for D1cc (the dose to 1 cm3of the volume) of the proximal bronchial tree (p> 0.05). At least 57.6% of mean heart dose, 28.8% of mean lung dose and 19.1% of mean esophageal dose were reduced compared with IMRT. The mean NTCP of radiation-induced pneumonitis (RP) in the ipsilateral lung was 0.39 ± 0.33 (0.39 ± 0.31) in IMPT plans and 0.36 ± 0.32 (0.35 ± 0.30) in IMCT plans compared with 0.66 ± 0.30 (0.64 ± 0.28) in IMRT plans by LKB (logistic) models. The target dose could be escalated to 78.3/76.9 Gy (RBE) in IMPT/IMCT plans compared with 61.7 Gy (RBE) in IMRT plans when 0.50 of NTCP in terms of RP in the ipsilateral lung was applied.Significance. This study presents the potential of better control of the side effects and improvement of local control originating from the dosimetric advantage with the application of IMPT and IMCT with the SIB technique for locally advanced lung cancer, even with limited beam directions.

Effects of nuclear interaction corrections and trichrome fragment spectra modelling on dose and linear energy transfer distributions in carbon ion radiotherapy

Background and Purpose

Nuclear interaction correction (NIC) and trichrome fragment spectra modelling improve relative biological effectiveness-weighted dose (DRBE) and dose-averaged linear energy transfer (LETd) calculation for carbon ions. The effect of those novel approaches on the clinical dose and LET distributions was investigated.

Materials and Methods

The effect of the NIC and trichrome algorithm was assessed, creating single beam plans for a virtual water phantom with standard settings and NIC + trichrome corrections. Reference DRBE and LETd distributions were simulated using FLUKA version 2021.2.9. Thirty clinically applied scanned carbon ion treatment plans were recalculated applying NIC, trichrome and NIC + trichrome corrections, using the LEM low dose approximation and compared to clinical plans (base RS). Four treatment sites were analysed: six prostate adenocarcinoma, ten head and neck, nine locally advanced pancreatic adenocarcinoma and five sacral chordoma. The FLUKA and clinical plans were compared in terms of DRBE deviations for D98%, D50%, D2% for the clinical target volume (CTV) and D50% in ring-like dose regions retrieved from isodose curves in base RS plans. Additionally, region-based median LETd deviations and global gamma parameters were evaluated.

Results

Dose deviations comparing base RS and evaluation plans were within ± 1% supported by γ-pass rates over 97% for all cases. No significant LETd deviations were reported in the CTV, but significant median LETd deviations were up to 80% for very low dose regions.

Conclusion

Our results showed improved accuracy of the predicted DRBE and LETd. Considering clinically relevant constraints, no significant modifications of clinical protocols are expected with the introduction of NIC + trichrome.

Computational approaches in the estimation of radiobiological damage for human-malignant cells irradiated with clinical proton and carbon beams

PURPOSE

The use of Monte Carlo (MC) simulations capable of reproducing radiobiological effects of ionising radiation on human cell lines is of great importance, especially for cases involving protons and heavier ion beams. In the latter, huge uncertainties can arise mainly related to the effects of the secondary particles produced in the beam-tissue interaction. This paper reports on a detailed MC study performed using Geant4-based approach on three cancer cell lines, the HTB-177, CRL-5876 and MCF-7, that were previously irradiated with therapeutic proton and carbon ion beams.

METHODS

A Geant4-based approach used jointly with analytical calculations has been developed to provide a more realistic estimation of the radiobiological damage produced by proton and carbon beams in tissues, reproducing available data obtained from in vitro cell irradiations. The MC "Hadrontherapy" Geant4 application and the Local Effect Model: LEM I, LEM II and LEM III coupled with the different numerical approaches: RapidRusso (RR) and RapidScholz (RS) were used in the study.

RESULTS

Experimental survival curves are compared with those evaluated using the highlighted Geant4 MC-based approach via chi-square statistical analysis, for the combinations of radiobiological models and numerical approaches, as outlined above.

CONCLUSION

This study has presented a comparison of the survival data from MC simulations to experimental survival data for three cancer cell lines. An overall best level of agreement was obtained for the HTB-177 cells.

Dosimetric rationale and preliminary experience in proton plus carbon-ion radiotherapy for esophageal carcinoma: a retrospective analysis

BACKGROUND

Concurrent chemoradiotherapy has been standard of care for unresectable esophageal carcinoma. There were no reports on proton radiotherapy (PRT) plus carbon-ion radiotherapy (CIRT) with pencil beam scanning (PBS) for esophageal carcinoma. This study evaluated the tolerability and efficiency of proton and sequential carbon-ion boost radiotherapy for esophageal carcinoma.

METHODS

From April 2017 to July 2020, 20 patients with primary esophageal carcinoma at stages II-IV were treated with PRT plus sequential CIRT with PBS. A median relative biological effectiveness-weighted PRT dose of 50 Gy in 25 fractions, and a sequential CIRT dose of 21 Gy in 7 fractions were delivered. Respiratory motion management was used if the tumor moved > 5 mm during the breathing cycle. A dosimetric comparison of photon intensity-modulated radiotherapy (IMRT), PRT, and CIRT was performed. The median times and rates of survivals were estimated using the Kaplan-Meier method. Comparison of the dose-volume parameters of the organs at risk employed the Wilcoxon matched-pairs test.

RESULTS

Twenty patients (15 men and 5 women, median age 70 years) were included in the analysis. With a median follow-up period of 25.0 months, the 2-year overall survival and progression-free survival rates were 69.2% and 57.4%, respectively. The patients tolerated radiotherapy and chemotherapy well. Grades 1, 2, 3, and 4 acute hematological toxicities were detected in 25%, 30%, 10%, and 30% of patients, respectively. Grades 3-5 acute non-hematological toxicities were not observed. Late toxicity events included grades 1, 2, and 3 in 50%, 20%, and 10% (pulmonary and esophageal toxicity in each) of patients. Grades 4-5 late toxicities were not noted. PRT or CIRT produced lower doses to organs at risk than did photon IMRT, especially the maximum dose delivered to the spinal cord and the mean doses delivered to the lungs and heart.

CONCLUSIONS

PRT plus CIRT with PBS appears to be a safe and effective treatment for esophageal carcinoma. PRT and CIRT delivered lower doses to organs at risk than did photon IMRT. Further investigation is warranted.

Monte Carlo simulations of cell survival in proton SOBP

Objective. The objective of this study is to develop a multi-scale modeling approach that accurately predicts radiation-induced DNA damage and survival fraction in specific cell lines.Approach. A Monte Carlo based simulation framework was employed to make the predictions. The FLUKA Monte Carlo code was utilized to estimate absorbed doses and fluence energy spectra, which were then used in the Monte Carlo Damage Simulation code to compute DNA damage yields in Chinese hamster V79 cell lines. The outputs were converted into cell survival fractions using a previously published theoretical model. To reduce the uncertainties of the predictions, new values for the parameters of the theoretical model were computed, expanding the database of experimental points considered in the previous estimation. Simulated results were validated against experimental data, confirming the applicability of the framework for proton beams up to 230 MeV. Additionally, the impact of secondary particles on cell survival was estimated.Main results. The simulated survival fraction versus depth in a glycerol phantom is reported for eighteen different configurations. Two proton spread out Bragg peaks at several doses were simulated and compared with experimental data. In all cases, the simulations follow the experimental trends, demonstrating the accuracy of the predictions up to 230 MeV.Significance. This study holds significant importance as it contributes to the advancement of models for predicting biological responses to radiation, ultimately contributing to more effective cancer treatment in proton therapy.

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.

Evaluation of proton and carbon ion beam models in TReatment Planning for Particles 4D (TRiP4D) referring to a commercial treatment planning system

PURPOSE

To investigate the accuracy of the treatment planning system (TPS) TRiP4D in reproducing doses computed by the clinically used TPS SyngoRT.

METHODS

Proton and carbon ion beam models in TRiP4D were converted from SyngoRT. Cubic plans with different depths in a water-tank phantom (WP) and previously treated and experimentally verified patient plans from SyngoRT were recalculated in TRiP4D. The target mean dose deviation (ΔD_mean,T) and global gamma index (2%-2 mm for the absorbed dose and 3%-3mm for the RBE-weighted dose with 10% threshold) were evaluated.

RESULTS

The carbon and proton absorbed dose gamma passing rates (γ-PRs) were ≥99.93% and ΔD_mean,T smaller than -0.22%. On average, the RBE-weighted dose D_mean,T was -1.26% lower for TRiP4D than SyngoRT for cubic plans. In TRiP4D, the faster analytical 'low dose approximation' (Krämer, 2006) was used, while SyngoRT used a stochastic implementation (Krämer, 2000). The average ΔD_{mean, T} could be reduced to -0.59% when applying the same biological effect calculation algorithm. However, the dose recalculation time increased by a factor of 79-477. ΔD_mean,T variation up to -2.27% and -2.79% was observed for carbon absorbed and RBE-weighted doses in patient plans. The γ-PRs were ≥93.92% and ≥91.83% for patient plans, except for one proton beam with a range shifter (γ-PR of 64.19%).

CONCLUSION

The absorbed dose between TRiP4D and SyngoRT were identical for both proton and carbon ion plans in the WP. Compared to SyngoRT, TRiP4D underestimated the target RBE-weighted dose; however more efficient in RBE-weighted dose calculation. Large variation for proton beam with range shifter was observed. TRiP4D will be used to evaluate doses delivered to moving targets. Uncertainties inherent to the 4D-dose reconstruction calculation are expected to be significantly larger than the dose errors reported here. For this reason, the residual differences between TRiP4D and SyngoRT observed in this study are considered acceptable. The study was approved by the Institutional Research Board of Shanghai Proton and Heavy Ion Center (approval number SPHIC-MP-2020-04, RS).

Creation, evolution, and future challenges of ion beam therapy from a medical physicist's viewpoint (Part 2). Chapter 2. Biophysical model, treatment planning system and image guided radiotherapy

When an ion beam penetrates deeply into the body, its kinetic energy decreases, and its biological effect increases due to the change of the beam quality. To give a uniform biological effect to the target, it is necessary to reduce the absorbed dose with the depth. A bio-physical model estimating the relationship between ion beam quality and biological effect is necessary to determine the relative biological effectiveness (RBE) of the ion beam that changes with depth. For this reason, Lawrence Berkeley Laboratory, National Institute of Radiological Sciences (NIRS) and GSI have each developed their own model at the starting of the ion beam therapy. Also, NIRS developed a new model at the starting of the scanning irradiation. Although the Local Effect Model (LEM) at the GSI and the modified Microdosimetric Kinetic Model (MKM) at the NIRS, the both are currently used, can similarly predict radiation quality-induced changes in surviving fraction of cultured cell, the clinical RBE-weighted doses for the same absorbed dose are different. This is because the LEM uses X-rays as a reference for clinical RBE, whereas the modified MKM uses carbon ion beam as a reference and multiplies it by a clinical factor of 2.41. Therefore, both are converted through the absorbed dose. In PART 2, I will describe the development of such a bio-physical model, as well as the birth and evolution of a treatment planning system and image guided radiotherapy.

State-of-the-art and potential of experimental microdosimetry in ion-beam therapy

In radiotherapy, radiation-quality should be an expression of the biological and physical characteristics of ionizing radiation such as spatial distribution of ionization or energy deposition. Linear energy transfer (LET) and lineal energy (y) are two descriptors used to quantify the radiation quality. These two quantities are connected and exhibit similar features. In ion-beam therapy (IBT), lineal energy can be measured with microdosimeters, which are specifically designed to cope with the high fluence of particles in clinical beams, while the quantification of LET is generally based on calculations. In pre-clinical studies, microdosimetric spectra are used for the indirect determination of relative biological effectiveness (RBE), e.g., using the microdosimetric kinetic model (MKM) or biophysical response functions. In this context it is important to consider saturation effects, which occur when the highest values of y become less biologically relevant compared to the relative contribution they make to the physical dose. Recent clinical data suggests that local tumor control and normal tissue effects can be linked to macroscopic and microscopic dosimetry parameters. In particular, positive clinical outcomes have been correlated to the highest LET values in the density distribution, and there is no evident link to the saturation discussed above. A systematic collection of microdosimetric information in combination with clinical data in retrospective studies may clarify the role of radiation quality at the highest LET. In the clinical setting, microdosimetry is not widely used yet, despite its potential to be linked with LET by experimentally-determined y values. Through this connection, both play an important role in complex therapy techniques such as intensity modulated particle therapy (IMPT), LET-painting and multi-ion optimization. This review summarizes the current state of microdosimetry for IBT and its potential, as well as research and development needed to make experimental microdosimetry a mature procedure in a clinical context.

Formation of spread-out Bragg peak for helium-ion beam using microdosimetric kinetic model

PURPOSE

To evaluate the applicability of microdosimetric kinetic model (MKM) to helium-ion therapy by forming a spread-out Bragg peak (SOBP) of a helium-ion beam using the MKM developed for carbon-ion radiotherapy and confirming the predictions in biological experiments.

METHODS

Using a ridge filter, a 90-mm wide SOBP for a 210 MeV/u helium-ion beam was created in a broad beam delivery system. The ridge filter was designed such that a uniform biological response was achieved with a cell survival rate of 7% over the SOBP region. Biological experiments were then performed using the SOBP beam in a human salivary gland (HSG) cell line to measure the cell survival rates.

RESULTS

The biological responses were uniform in the SOBP region, as expected by the MKM; however, the mean of the measured cell survival rates was (11.2 ± 0.6) % in the SOBP region, which was 60% higher than the designed rate. When investigating the biological parameters of the HSG cell line used in the experiments, we found that they were altered slightly from the MKM parameters used for carbon-ion radiotherapy. The new β parameter reproduced the measured survival rates within 6.5% in the SOBP region.

CONCLUSION

We produced biologically uniform SOBP using MKM for carbon-ion radiotherapy. The measured survival rates in the SOBP region were higher than expected, and the survival rates were reproduced by modifying the MKM parameter. This study was limited to one SOBP, and further investigations are required to prove that MKM is generally applicable to helium-ion radiotherapy.

Accumulation of sublethal radiation damage and its effect on cell survival

Objective.Determine the extent of sublethal radiation damage (SRD) in a cell population that received a given dose of radiation and the impact of this damage on cell survival.Approach.We developed a novel formalism to account for accumulation of SRD with increasing dose. It is based on a very general formula for cell survival that correctly predicts the basic properties of cell survival curves, such as the transition from the linear-quadratic to a linear dependence at high doses. Using this formalism we analyzed extensive experimental data for photons, protons and heavy ions to evaluate model parameters, quantify the extent of SRD and its impact on cell survival.Main results.Significant accumulation of SRD begins at doses below 1 Gy. As dose increases, so does the number of damaged cells and the amount of SRD in individual cells. SRD buildup in a cell increases the likelihood of complex irrepairable damage. For this reason, during a dose fraction delivery, each dose increment makes cells more radiosensitive. This gradual radosensitization is evidenced by the increasing slope of survival curves observed experimentally. It continues until the fraction is delivered, unless radiosensitivity reaches its maximum first. The maximum radiosensitivity is achieved when SRD accumulated in most cells is the maximum damage they can repair. After this maximum is reached, the slope of a survival curve, logarithm of survival versus dose, becomes constant, dose independent. The survival curve becomes a straight line, as experimental data at high doses show. These processes are random. They cause large cell-to-cell variability in the extent of damage and radiosensitivity of individual cells.Significance.SRD is in effect a radiosensitizer and its accumulation is a significant factor affecting cell survival, especially at high doses. We developed a novel formalism to study this phenomena and reported pertinent data for several particle types.

Calibration and evaluation of the relative biological effectiveness for carbon-ion radiotherapy in a new relative to a clinically applied treatment planning system

BACKGROUND

The study objective was to validate the relative biological effectiveness (RBE) in RayStation for carbon-ion radiotherapy (CIRT) using the Syngo treatment planning system as reference.

METHODS

Local effect model I was established in RayStation (Ray-LEM) with the same parameters as in LEM I in Syngo (Syngo-LEM). Three cube plans covering most of the tumors treated at our center were generated with Syngo-LEM. Ray-LEM re-calculated the Syngo plans and compared the RBEs to the Syngo counterparts. The results showed that RayStation RBE was smaller than Syngo RBE. To ensure that Ray-LEM reproduced Syngo RBE, the observed deviations were used to scale the maximum RBE (RBEmax) in Ray-LEM. After this calibration, we further compared the RayStation RBE to Syngo RBE using additional plans in both homogeneous phantoms and patients, to ensure that the calibrated Ray-LEM reproduced Syngo RBE even with more complex planning features.

RESULTS

The calibration increased the RBEmax by 2.3% to raise the Ray-LEM RBE. The target mean RBE deviations in the phantom evaluation plans were median: 0.0 (minimum: - 1.1 to maximum: 0.7) %, and the target mean RBE deviations of the clinical target volumes of 16 patient cases were - 0.4 (- 1.5 to 0.2) %.

CONCLUSIONS

The residual RBE difference between RayStation and Syngo was found to be ≤ 1.0%. Thus, we can propose to use RayStation for clinical CIRT treatment planning. However, the potential differences due to the absorbed beam model warrants further exploration.

A Dosiomics Analysis Based on Linear Energy Transfer and Biological Dose Maps to Predict Local Recurrence in Sacral Chordomas after Carbon-Ion Radiotherapy

Carbon Ion Radiotherapy (CIRT) is one of the most promising therapeutic options to reduce Local Recurrence (LR) in Sacral Chordomas (SC). The aim of this work is to compare the performances of survival models fed with dosiomics features and conventional DVH metrics extracted from relative biological effectiveness (RBE)-weighted dose (D_RBE) and dose-averaged Linear Energy Transfer (LET_d) maps, towards the identification of possible prognostic factors for LR in SC patients treated with CIRT. This retrospective study included 50 patients affected by SC with a focus on patients that presented a relapse in a high-dose region. Survival models were built to predict both LR and High-Dose Local Recurrencies (HD-LR). The models were evaluated through Harrell Concordance Index (C-index) and patients were stratified into high/low-risk groups. Local Recurrence-free Kaplan-Meier curves were estimated and evaluated through log-rank tests. The model with highest performance (median(interquartile-range) C-index of 0.86 (0.22)) was built on features extracted from LET_d maps, with D_RBE models showing promising but weaker results (C-index of 0.83 (0.21), 0.80 (0.21)). Although the study should be extended to a wider patient population, LET_d maps show potential as a prognostic factor for SC HD-LR in CIRT, and dosiomics appears to be the most promising approach against more conventional methods (e.g., DVH-based).

Calculating dose-averaged linear energy transfer in an analytical treatment planning system for carbon-ion radiotherapy

BACKGROUND

Compelling evidence shows the association between the relative biological effectiveness (RBE) of carbon-ion radiotherapy (CIRT) and the dose averaged linear energy transfer (LETd). However, the ability to calculate the LETd in commercially available treatment planning systems (TPS) is lacking.

PURPOSE

This study aims to develop a method of calculating the LETd of CIRT plans that could be robustly carried out in RayStation (V10B, Raysearch, Sweden).

METHODS

The calculation used the fragment spectra in RayStation for the CIRT treatment planning. The dose-weighted averaging procedure was supported by the microdosimetric kinetic model (MKM). The MKM-based pencil beam dose engine (PBA, v4.2) for calculating RBE-weighted doses was reformulated to become a LET-weighted calculating engine. A separate module was then configured to inversely calculate the LETd from the absorbed dose of a plan and the associated fragment spectra. In this study, the ion and energy-specific LET table in the LETd module was further matched with the values decoded from the baseline data of the Syngo TPS (V13C, Siemens, Germany). The LETd distributions of several monoenergetic and modulated beams were calculated and validated against the values derived from the Syngo TPS and the published data.

RESULTS

The differences in LETds of the monoenergetic beams between the new method and the traditional method were within 3% in the entrance and Bragg-peak regions. However, a larger difference was observed in the distal region. The results of the modulated beams were in good agreement with the works from the published literature.

CONCLUSIONS

The method presented herein reformulates the MKM dose engine in the RayStation TPS to inversely calculate LETds. The robustness and accuracy were demonstrated.

CONVERSION OF DOSE DISTRIBUTION TO CELL SURVIVAL FRACTION THROUGH DNA DAMAGE: A MONTE CARLO STUDY

Ionizing radiation plays an important role in cancer treatment. Radiation is able to damage the genetic material of cells, blocking their ability to divide and proliferate further. Since radiation affects both healthy and malignant tissues, for all radiation treatments, the design of an accurate treatment plan is fundamental. Usually, weight factors, such as the relative biological effectiveness, are applied to estimate the impact of the kind of radiation and the irradiated medium on the dose deposition. However, these factors can only provide a partial estimation of the real effect on tissues. In this work, a flexible system that is able to predict cell survival fractions according to the planned dose distribution is presented. Dose deposition and subsequent DNA damage were simulated with a multi-scale modeling approach by first applying the FLUKA Monte Carlo (MC) code to estimate the absorbed doses and fluence energy spectra and then using the MC Damage Simulation code to compute the DNA damage yields. Lastly, the results are converted into cell survival fraction using a theoretical model. The comparisons between the simulated survival fractions with experimental data are reported for a proton spread out Bragg peak at several doses. The presented approach helps to elucidate radiobiological responses along the Bragg curve and has the flexibility to be extended to a wide range of situations of clinical interest.

3D high resolution clonogenic survival measurement of xrs-5 cells in low-dose region of carbon ion plans

PURPOSE

In this study, we performed biological verification measurements of cell survival of a ¹²C ion irradiation plan employing a high-resolution 3D culture setup. This allowed, in particular, to access the cell inactivation in the low-dose regions close to the target area.

MATERIALS AND METHODS

We established the protocol for a 3D culture setup where xrs-5 cells were grown inside a layered matrigel structure in 384-well plates. Their radiosensitivity to conventional and ¹²C ion radiation was evaluated by irradiating them either with 250 kV X-rays at GSI or with monoenergetic ¹²C beams of 110 MeV/u at MIT, and compared with those of monolayers. A treatment plan for a rectangular target was prepared using the GSI research treatment planning system TRiP98. xrs-5 cells were seeded in the matrigel-based setup and irradiated in dose fall-off regions using active scanning ¹²C ion beams. In addition, film dosimetry utilizing radiochromic EBT3 film has been performed to assess the field homogeneity downstream of 384-well V-bottom plates with or without additional agarose coating of the well plate bottom.

RESULTS

Dose response curves following X-ray and ¹²C ion irradiation had linear shape and showed a significant decrease in survival fraction at even moderate doses. Survival measurements in the low-dose regions of the plan for the extended target showed good agreement to the predicted survival fraction. The irradiated film profiles yielded a flat dose distribution without apparent artifacts or inhomogeneities for well plates both with and without agarose coating, confirming the suitability of the experimental setup.

CONCLUSIONS

We conclude that the V-bottom 384-well plates in combination with the radiation-sensitive xrs-5 cell line constitute a suitable radiobiological verification tool which can be used especially for low doses. Furthermore, the measured survival of xrs-5 cells show a good agreement with the expected survival in the low-dose out-of-field regions, both laterally and downstream of the target.

Spot-scanning hadron arc (SHArc) therapy: A proof of concept using single and multi-ion strategies with helium, carbon, oxygen and neon ions

PURPOSE

To present particle arc therapy treatments using single and multi-ion therapy optimization strategies with helium (⁴ He), carbon (¹² C), oxygen (¹⁶ O) and neon (²⁰ Ne) ion beams.

METHODS AND MATERIALS

An optimization procedure and workflow were devised for spot-scanning hadron arc therapy (SHArc) treatment planning in the PRECISE (PaRticle thErapy using single and Combined Ion optimization StratEgies) treatment planning system (TPS). Physical and biological beam models were developed for helium, carbon, oxygen and neon ions via FLUKA MC simulation. SHArc treatments were optimized using both single ion (¹² C, ¹⁶ O, or ²⁰ Ne) and multi-ion therapy (¹⁶ O+⁴ He or ²⁰ Ne+⁴ He) applying variable relative biological effectiveness (RBE) modeling using a modified microdosimetric kinetic model (mMKM) with (α/β)_x values of 2Gy, 5Gy and 3.1Gy respectively, for glioblastoma, pancreatic adenocarcinoma, and prostate adenocarcinoma patient cases. Dose, effective dose, linear energy transfer (LET) and RBE were computed with the GPU-accelerated dose engine FRoG and dosimetric/biophysical attributes were evaluated in the context of conventional particle and photon-based therapies (e.g., volumetric modulated arc therapy [VMAT]).

RESULTS

All SHArc plans met the target optimization goals (3GyRBE) and demonstrated increased target conformity and substantially lower low-dose bath to surrounding normal tissues than VMAT. SHArc plans using a single ion species (¹² C, ¹⁶ O, or ²⁰ Ne) exhibited favorable LET distributions with the highest-LET components centralized in the target volume, with values ranging from ∼80-170keV/μm, ∼130-220keV/μm and ∼180-350keV/μm, for ¹² C, ¹⁶ O, or ²⁰ Ne, respectively, exceeding mean target LET of conventional particle therapy (¹² C:∼60, ¹⁶ O:∼78 ²⁰ Ne:∼100 keV/μm). Multi-ion therapy with SHArc delivery (SHArc_MIT ) provided a similar level of target LET enhancement as SHArc compared to conventional planning, however, with additional benefits of homogenous physical dose and RBE distributions.

CONCLUSION

Here, we demonstrate that arc delivery of light and heavy ion beams, using either a single ion species (¹² C, ¹⁶ O, or ²⁰ Ne) or combining two ions in a single fraction (¹⁶ O+⁴ He or ²⁰ Ne+⁴ He), affords enhanced physical and biological distributions (e.g., LET) compared with conventional delivery with photons or particle beams. SHArc marks the first single and multi-ion arc therapy treatment optimization approach using light and heavy ions. This article is protected by copyright. All rights reserved.

Including Volume Effects in Biological Treatment Plan Optimization for Carbon Ion Therapy: Generalized Equivalent Uniform Dose-Based Objective in TRiP98

We describe a way to include biologically based objectives in plan optimization specific for carbon ion therapy, beyond the standard voxel-dose-based criteria already implemented in TRiP98, research planning software for ion beams. The aim is to account for volume effects—tissue architecture-dependent response to damage—in the optimization procedure, using the concept of generalized equivalent uniform dose (gEUD), which is an expression to convert a heterogeneous dose distribution (e.g., in an organ at risk (OAR)) into a uniform dose associated with the same biological effect. Moreover, gEUD is closely related to normal tissue complication probability (NTCP). The multi-field optimization problem here takes also into account the relative biological effectiveness (RBE), which in the case of ion beams is not factorizable and introduces strong non-linearity. We implemented the gEUD-based optimization in TRiP98, allowing us to control the whole dose–volume histogram (DVH) shape of OAR with a single objective by adjusting the prescribed gEUD₀ and the volume effect parameter a, reducing the volume receiving dose levels close to mean dose when a = 1 (large volume effect) while close to maximum dose for a >> 1 (small volume effect), depending on the organ type considered. We studied the role of gEUD₀ and a in the optimization, and we compared voxel-dose-based and gEUD-based optimization in chordoma cases with different anatomies. In particular, for a plan containing multiple OARs, we obtained the same target coverage and similar DVHs for OARs with a small volume effect while decreasing the mean dose received by the proximal parotid, thus reducing its NTCP by a factor of 2.5. Further investigations are done for this plan, considering also the distal parotid gland, obtaining a NTCP reduction by a factor of 1.9 for the proximal and 2.9 for the distal one. In conclusion, this novel optimization method can be applied to different OARs, but it achieves the largest improvement for organs whose volume effect is larger. This allows TRiP98 to perform a double level of biologically driven optimization for ion beams, including at the same time RBE-weighted dose and volume effects in inverse planning. An outlook is presented on the possible extension of this method to the target.

Management of Motion and Anatomical Variations in Charged Particle Therapy: Past, Present, and Into the Future

The major aim of radiation therapy is to provide curative or palliative treatment to cancerous malignancies while minimizing damage to healthy tissues. Charged particle radiotherapy utilizing carbon ions or protons is uniquely suited for this task due to its ability to achieve highly conformal dose distributions around the tumor volume. For these treatment modalities, uncertainties in the localization of patient anatomy due to inter- and intra-fractional motion present a heightened risk of undesired dose delivery. A diverse range of mitigation strategies have been developed and clinically implemented in various disease sites to monitor and correct for patient motion, but much work remains. This review provides an overview of current clinical practices for inter and intra-fractional motion management in charged particle therapy, including motion control, current imaging and motion tracking modalities, as well as treatment planning and delivery techniques. We also cover progress to date on emerging technologies including particle-based radiography imaging, novel treatment delivery methods such as tumor tracking and FLASH, and artificial intelligence and discuss their potential impact towards improving or increasing the challenge of motion mitigation in charged particle therapy.

Estimate of the Biological Dose in Hadrontherapy Using GATE

Simple Summary

This study presents the implementation of a biological dose module using the Monte Carlo software, GATE. Both mMKM and NanOx biophysics models of cell survival predictions were used as input. The code was validated in terms of biological dose, relative biological effectiveness and cell survival against experimental data from the HIMBC (Hyogo, Japan) ion beam line.

Abstract

For the evaluation of the biological effects, Monte Carlo toolkits were used to provide an RBE-weighted dose using databases of survival fraction coefficients predicted through biophysical models. Biophysics models, such as the mMKM and NanOx models, have previously been developed to estimate a biological dose. Using the mMKM model, we calculated the saturation corrected dose mean specific energy z1D* (Gy) and the dose at 10% D₁₀ for human salivary gland (HSG) cells using Monte Carlo Track Structure codes LPCHEM and Geant4-DNA, and compared these with data from the literature for monoenergetic ions. These two models were used to create databases of survival fraction coefficients for several ion types (hydrogen, carbon, helium and oxygen) and for energies ranging from 0.1 to 400 MeV/n. We calculated α values as a function of LET with the mMKM and the NanOx models, and compared these with the literature. In order to estimate the biological dose for SOBPs, these databases were used with a Monte Carlo toolkit. We considered GATE, an open-source software based on the GEANT4 Monte Carlo toolkit. We implemented a tool, the BioDoseActor, in GATE, using the mMKM and NanOx databases of cell survival predictions as input, to estimate, at a voxel scale, biological outcomes when treating a patient. We modeled the HIBMC 320 MeV/u carbon-ion beam line. We then tested the BioDoseActor for the estimation of biological dose, the relative biological effectiveness (RBE) and the cell survival fraction for the irradiation of the HSG cell line. We then tested the implementation for the prediction of cell survival fraction, RBE and biological dose for the HIBMC 320 MeV/u carbon-ion beamline. For the cell survival fraction, we obtained satisfying results. Concerning the prediction of the biological dose, a 10% relative difference between mMKM and NanOx was reported.

Dose- and LET-dependent changes in mouse skin contracture up to a year after either single dose or fractionated doses of carbon ion or gamma rays

Time dependence of relative biological effectiveness (RBE) of carbon ions for skin damage was investigated to answer the question of whether the flat distribution of biological doses within a Spread-Out Bragg peak (SOBP) which is designed based on in vitro cell kill could also be flat for in vivo late responding tissue. Two spots of Indian ink intracutaneously injected into the legs of C3H mice were measured by calipers. An equieffective dose to produce 30% skin contraction was calculated from a dose-response curve and used to calculate the RBE of carbon ion beams. We discovered skin contraction progressed after irradiation and then reached a stable/slow progression phase. Equieffective doses decreased with time and the decrease was most prominent for gamma rays and least prominent for 100 keV/μm carbon ions. Survival parameter of alpha but not beta in the linear-quadratic model is closely related to the RBE of carbon ions. Biological doses within the SOBP increased with time but their distribution was still flat up to 1 year after irradiation. The outcomes of skin contraction studies suggest that (i) despite the higher RBE for skin contracture after carbon ions compared to gamma rays, gamma rays can result in a more severe late effect of skin contracture. This is due to the carbon effect saturating at a lower dose than gamma rays, and (ii) the biological dose distribution throughout the SOBP remains approximately the same even one year after exposure.

The history of ion beam therapy in Germany

The advantageous depth dose profile of ion beams together with state of the art beam delivery and treatment planning systems allow for highly conformal tumor treatments in patients. First treatments date back to 1954 at the Lawrence Berkeley Laboratory (LBL) and in Europe, ion beam therapy started in the mid-1990s at the Paul-Scherrer Institute (PSI) with protons and at the Helmholtz Center for Heavy Ion Research (GSI) with carbon ions, followed by the Heidelberg Ion Therapy Center (HIT) in Heidelberg. This review describes the historical development of ion beam therapy in Germany based on the pioneering work at LBL and in the context of simultaneous developments in other countries as well as recent developments.

In Silico Feasibility Study of Carbon Ion Radiotherapy With Simultaneous Integrated Boost for Head and Neck Adenoid Cystic Carcinoma

Purpose

In carbon ion radiotherapy (CIRT), a simultaneous integrated boost (SIB) approach has not been fully exploited so far. The feasibility of a CIRT-SIB strategy for head and neck adenoid cystic carcinoma (ACC) patients was investigated in order to improve treatment planning dose distributions.

Methods and Materials

CIRT plans of 10 ACC patients treated at the National Center for Oncological Hadrontherapy (CNAO, Pavia, Italy) with sequential boost (SEQ) irradiation and prescription doses of 41.0 Gy [relative biological effectiveness (RBE)]/10 fractions to low-risk (LR) clinical target volume (CTV) plus 24.6 Gy(RBE)/6 fractions to the high-risk (HR) CTV were re-planned with two SIB dose levels to the LR-CTV, namely, 48.0 Gy(RBE) and 54.4 Gy(RBE). While planning with SIB, the HR-CTV coverage had higher priority, with fixed organ-at-risk dose constraints among the SIB and SEQ plans. The homogeneity and conformity indexes were selected for CTV coverage comparison. The biologically effective dose (BED) was calculated to compare the different fractionation schemes.

Results

Comparable HR-CTV coverage was achieved with the treatment approaches, while superior conformality and homogeneity were obtained with the SIB technique in both CTVs. With the SEQ, SIB_48.0, and SIB_54.4, the LR-CTV median doses were respectively 50.3%, 11.9%, and 6.0% higher than the prescriptions. Significant reductions of the median and near-maximum BEDs were achieved with both SIB dose levels in the LR-CTV.

Conclusions

The SIB approach resulted in highly conformal dose distributions with the reduction of the unintended dose to the LR-CTV. A prescription dose range for the LR-CTV will be clinically defined to offer tailored personalized treatments, according to the clinical and imaging characteristics of the patients.

Therapeutic Efficacy of Variable Biological Effectiveness of Proton Therapy in U-CH2 and MUG-Chor1 Human Chordoma Cell Death

Simple Summary

Chordoma is a rare, slow-growing cancer of the spinal cord. Photon radiation therapy and surgery are the standard of care for chordoma. Proton radiation therapy has become an increasingly common treatment in comparison to photon radiation therapy due to the ability to reduce off-target radiation dose. However, there is still an increased risk of toxicity to the surrounding critical structures that lead to poor treatment outcomes. Moreover, the biologic effectiveness of protons to sterilize chordoma cells remains uncertain and likely varies according to the proton energy spectrum throughout the proton field. We aim to investigate the tumoricidal properties of proton radiation therapy at the middle and end of the proton radiation field and elucidate variations in the relative biological effectiveness for chordoma cells. Our study helps quantify the therapeutic value of treating chordoma near the end of the proton field, where linear energy transfer is relatively high.

Abstract

Background: Chordoma is a cancer of spinal cord, skull base, and sacral area. Currently, the standard of care to treat chordoma is resection followed by radiation therapy. Since, chordoma is present in the spinal cord and these are very sensitive structures and often complete removal by surgery is not possible. As a result, chordoma has a high chance of recurrence and developing resistance to radiation therapy. In addition, treatment of chordoma by conventional radiation therapy can also damage normal tissues surrounding chordoma. Thus, current therapeutic options to treat chordoma are insufficient and novel therapies are desperately needed to treat locally advanced and metastatic chordoma. (2) Methods: In the present investigation, human chordoma cell lines of sacral origin MUG-Chor1 and U-CH2 were cultured and irradiated with Proton Beam Radiation using the clinical superconducting cyclotron and pencil-beam (active) scanning at Middle and End of the Spread-Out Bragg Peak (SOBP). Proton radiation was given at the following doses: Mug-Chor1 at 0, 1, 2, 4, and 8 Gy and U-CH2 at 0, 4, 8, 12, and 16 Gy. These doses were selected based on a pilot study in our lab and attempted to produce approximate survival fractions in the range of 1, 0.9, 0.5, 0.1, and 0.01, respectively, chosen for linear quadratic model fitting of the dose response. (3) Results: In this study, we investigated relative biological effectiveness (RBE) of proton radiation at the end of Spread Out Bragg Peak assuming that the reference radiation is a proton radiation in the middle of the SOBP. We observed differences in the survival of both Human chordoma cell lines, U-CH2 and MUG-Chor1. The data showed that there was a significantly higher cell death at the end of the Bragg peak as compared to middle of the Bragg peak. Based on the linear quadratic (LQ) fit for cell survival we calculated the RBE between M-SOBP and E-SOBP at 95% CI level and it was observed that RBE was higher than 1 at E-SOBP and caused significantly higher cell killing. Proton field at E-SOBP caused complex DNA damage in comparison to M-EOBP and the genes such as DNA topoisomerase 1, GTSE1, RAD51B were downregulated in E-SOBP treated cells. Thus, we conclude that there seems to be substantial variation in RBE (1.3–1.7) at the E-SOBP compared with the M-SOBP.

Compensating for beam modulation due to microscopic lung heterogeneities in carbon ion therapy treatment planning

PURPOSE

To predict and mitigate for the degradation in physical and biologically effective dose distributions of particle beams caused by microscopic heterogeneities in lung tissue.

MATERIALS AND METHODS

The TRiP98 treatment planning system was adapted to account for the beam-modulating effect of heterogeneous lung tissue in physical and biological inverse treatment planning. The implementation employs an analytical model that derives the degradation from the established "modulation power" parameter P mod and the total water-equivalent thickness of lung parenchyma traversed by the beam. Beam modulation was reproduced through an on-the-fly convolution of the reference Bragg curve with Gaussian kernels depending on the modulation power of lung tissue (upstream). For biological doses, the degradation was determined by modulating dose-averaged α , β , and LET distributions. Carbon SOBP measurements behind lung substitute material were performed to validate the code. The implementation was then applied to a phantom and patient case.

RESULTS

Experimental results show the passage through a 20-cm Gammex LN300 slab led to a decrease in target coverage and broadening of the SOBP distal fall-off. However, dose coverage was regained through optimization. A good agreement between calculated and measured SOBPs was also found. In addition, a patient case study revealed a 3.2% decrease in D 95 from degradation ( P mod = 450 μ m), which was reduced to a 0.4% difference after optimization. Furthermore, widening of the RBE distribution beyond the target distal edge was observed. This implies an increased degradation in the biological dose, which could be harmful to healthy tissues distal to the target.

CONCLUSIONS

This is the first implementation capable of compensating for lung dose perturbations, which is more effective than margin extensions. A larger patient study is needed to examine the observed modulation in the RBE distribution and judge the clinical relevance also in IMPT, where margins might prove insufficient to recover target coverage.

Dosimetric Impact of Inter-Fraction Anatomical Changes in Carbon Ion Boost Treatment for High-Risk Prostate Cancer (AIRC IG 14300)

Rectum and bladder volumes play an important role in the dose distribution reproducibility in prostate cancer adenocarcinoma (PCa) radiotherapy, especially for particle therapy, where density variation can strongly affect the dose distribution. We investigated the reliability and reproducibility of our image-guided radiotherapy (IGRT) and treatment planning protocol for carbon ion radiotherapy (CIRT) within the phase II mixed beam study (AIRC IG 14300) for the treatment of high-risk PCa. In order to calculate the daily dose distribution, a set of synthetic computed tomography (sCT) images was generated from the cone beam computed tomography (CBCT) images acquired in each treatment session. Planning target volume (PTV) together with rectum and bladder volume variation was evaluated with sCT dose-volume histogram (DVH) metric deviations from the planning values. The correlations between the bladder and rectum volumes, and the corresponding DVH metrics, were also assessed. No significant difference in the bladder, rectum, and PTV median volumes between the planning computed tomography (pCT) and the sCT was found. In addition, no significant difference was assessed when comparing the average DVHs and median DVH metrics between pCT and sCT. Dose deviations determined by bladder and rectum filling variations demonstrated that dose distributions were reproducible in terms of both target coverage and organs at risk (OARs) sparing.

Technical Note: Field size analysis of patient-specific quality assurance in scanned carbon ion radiotherapy

PURPOSE

To evaluate the dose difference between measurement and double Gaussian beam model prediction according to the field size and correct the measurements in patient-specific quality assurance (QA).

METHODS

The field size dependence of the dose was evaluated with volumes of 20 × 20 × 80 mm³ , 40 × 40 × 80 mm³ , 60 × 60 × 80 mm³ , and 80 × 80 × 80 mm³ of 1 Gy uniform dose at three depths. Additional two 80 × 80 × 80 mm³ volumes of nonuniform fields were created: one high-dose field was given 1 Gy at the central 40 × 40 mm² and 0.5 Gy in its surrounding, and the other low-dose field was given 0.5 Gy in the middle and 1 Gy at the periphery. The dose in the center of the spread-out Bragg peak (SOBP) was measured in a water phantom and compared with the treatment planning system (TPS) predication. A field factor based on the two-dimensional (2D) dose distribution was proposed to estimate the field size. The field factor was first evaluated against the dose difference in the square fields, and then used to analyze and correct the patient-specific QA results.

RESULTS

TPS overestimated dose for fields smaller than 80 × 80 mm² . A practically positive correlation was observed between the measured dose and the field factor. In the patient-specific QA, measured doses were lower than the TPS predication as they were calculated a relatively small field factor. The corrected dose differences were no longer field factor dependent.

CONCLUSIONS

Using the proposed field factor, we have shown that all the measurements with a large dose deviation were due to the small-sized field. It is clinically relevant to take into consideration the field size in the QA analysis as long as the double Gaussian beam model being used for the dose calculation. Correction to the measurement can be made based on the field factor.

Adaptation of stochastic microdosimetric kinetic model to hypoxia for hypo-fractionated multi-ion therapy treatment planning

For hypo-fractionated multi-ion therapy (HFMIT), the stochastic microdosimetric kinetic (SMK) model had been developed to estimate the biological effectiveness of radiation beams with wide linear energy transfer (LET) and dose ranges. The HFMIT will be applied to radioresistant tumors with oxygen-deficient regions. The response of cells to radiation is strongly dependent on the oxygen condition in addition to radiation type, LET and absorbed dose. This study presents an adaptation of the SMK model to account for oxygen-pressure dependent cell responses, and develops the oxygen-effect-incorporated stochastic microdosimetric kinetic (OSMK) model. In the model, following assumptions were made: the numbers of radiation-induced sublethal lesions (double-strand breaks) are reduced due to lack of oxygen, and the numbers of oxygen-mediated lesions are reduced for radiation with high LET. The model parameters were determined by fitting survival data under aerobic and anoxic conditions for human salivary gland tumor cells and V79 cells exposed to helium-, carbon-, and neon-ion beams over the LET range of 18.5-654.0 keVμm^-1. The OSMK model provided good agreement with the experimental survival data of the cells with determination coefficients >0.9. In terms of oxygen enhancement ratio, the OSMK model reproduced the experimental data behavior, including slight dependence on particle type at the same LET. The OSMK model was then implemented into the in-house treatment planning software for the HFMIT to validate its applicability in clinical practice. A treatment plan with helium- and neon-ion beams was made for a pancreatic cancer case assuming an oxygen-deficient region within the tumor. The biological optimization based on the OSMK model preferentially placed the neon-ion beam to the hypoxic region, while it placed both helium- and neon-ion beams to the surrounding normoxic region. The OSMK model offered the accuracy and usability required for hypoxia-based biological optimization in HFMIT treatment planning.

Biological effectiveness of He-3 and He-4 ion beams for cancer hadrontherapy: a study based on the BIANCA biophysical model

While cancer therapy with protons and C-ions is continuously spreading, in the near future patients will be also treated with He-ions which, in comparison to photons, combine the higher precision of protons with the higher relative biological effectiveness (RBE) of C-ions. Similarly to C-ions, also for He-ions the RBE variation along the beam must be known as precisely as possible, especially for active beam delivery systems. In this framework the BIANCA biophysical model, which has already been applied to calculate the RBE along proton and C-ion beams, was extended to⁴He-ions and, following interface with the FLUKA code, was benchmarked against cell survival data on CHO normal cells and Renca tumour cells irradiated at different positions along therapeutic-like⁴He-ion beams at the Heidelberg Ion-beam Therapy centre, where the first He-ion patient will be treated soon. Very good agreement between simulations and data was obtained, showing that BIANCA can now be used to predict RBE following irradiation with all ion types that are currently used, or will be used soon, for hadrontherapy. Thanks to the development of a reference simulation database describing V79 cell survival for ion and photon irradiation, these predictions can be cell-type specific because analogous databases can be produced, in principle, for any cell line. Furthermore, survival data on CHO cells irradiated by a He-3 beam were reproduced to compare the biophysical properties of He-4 and He-3 beams, which is currently an open question. This comparison showed that, at the same depth, He-4 beams tend to have a higher RBE with respect to He-3 beams, and that this difference is also modulated by the considered physical dose, as well as the cell radiosensitivity. However, at least for the considered cases, no significant difference was found for the ratio between the RBE-weighted dose in the SOBP and that in the entrance plateau.

Dosimetric Validation of a System to Treat Moving Tumors Using Scanned Ion Beams That Are Synchronized With Anatomical Motion

Purpose

The purpose of this study was to validate the dosimetric performance of scanned ion beam deliveries with motion-synchronization to heterogenous targets.

Methods

A 4D library of treatment plans, comprised of up to 10 3D sub-plans, was created with robust and conventional 4D optimization methods. Each sub-plan corresponded to one phase of periodic target motion. The plan libraries were delivered to a test phantom, comprising plastic slabs, dosimeters, and heterogenous phantoms. This phantom emulated range changes that occur when treating moving tumors. Similar treatment plans, but without motion synchronization, were also delivered to a test phantom with a stationary target and to a moving target; these were used to assess how the target motion degrades the quality of dose distributions and the extent to which motion synchronization can improve dosimetric quality. The accuracy of calculated dose distributions was verified by comparison with corresponding measurements. Comparisons utilized the gamma index analysis method. Plan quality was assessed based on conformity, dose coverage, overdose, and homogeneity values, each extracted from calculated dose distributions.

Results

High pass rates for the gamma index analysis confirmed that the methods used to calculate and reconstruct dose distributions were sufficiently accurate for the purposes of this study. Calculated and reconstructed dose distributions revealed that the motion-synchronized and static deliveries exhibited similar quality in terms of dose coverage, overdose, and homogeneity for all deliveries considered. Motion-synchronization substantially improved conformity in deliveries with moving targets. Importantly, measurements at multiple locations within the target also confirmed that the motion-synchronized delivery system satisfactorily compensated for changes in beam range caused by the phantom motion. Specifically, the overall planning and delivery approach achieved the desired dose distribution by avoiding range undershoots and overshoots caused by tumor motion.

Conclusions

We validated a dose delivery system that synchronizes the movement of the ion beam to that of a moving target in a test phantom. Measured and calculated dose distributions revealed that this system satisfactorily compensated for target motion in the presence of beam range changes due to target motion. The implication of this finding is that the prototype system is suitable for additional preclinical research studies, such as irregular anatomic motion.

Biological Impact of Target Fragments on Proton Treatment Plans: An Analysis Based on the Current Cross-Section Data and a Full Mixed Field Approach

Simple Summary

Proton therapy is now an established external radiotherapy modality for cancer treatment. Clinical routine currently neglects the radiobiological impact of nuclear target fragments even if experimental evidences show a significant enhancement in cell-killing effect due to secondary particles. This paper quantifies the contribution of proton target fragments of different charge in different irradiation scenarios and compares the computationally predicted corrections to the overall biological dose with experimental data.

Abstract

Clinical routine in proton therapy currently neglects the radiobiological impact of nuclear target fragments generated by proton beams. This is partially due to the difficult characterization of the irradiation field. The detection of low energetic fragments, secondary protons and fragments, is in fact challenging due to their very short range. However, considering their low residual energy and therefore high LET, the possible contribution of such heavy particles to the overall biological effect could be not negligible. In this context, we performed a systematic analysis aimed at an explicit assessment of the RBE (relative biological effectiveness, i.e., the ratio of photon to proton physical dose needed to achieve the same biological effect) contribution of target fragments in the biological dose calculations of proton fields. The TOPAS Monte Carlo code has been used to characterize the radiation field, i.e., for the scoring of primary protons and fragments in an exemplary water target. TRiP98, in combination with LEM IV RBE tables, was then employed to evaluate the RBE with a mixed field approach accounting for fragments’ contributions. The results were compared with that obtained by considering only primary protons for the pristine beam and spread out Bragg peak (SOBP) irradiations, in order to estimate the relative weight of target fragments to the overall RBE. A sensitivity analysis of the secondary particles production cross-sections to the biological dose has been also carried out in this study. Finally, our modeling approach was applied to the analysis of a selection of cell survival and RBE data extracted from published in vitro studies. Our results indicate that, for high energy proton beams, the main contribution to the biological effect due to the secondary particles can be attributed to secondary protons, while the contribution of heavier fragments is mainly due to helium. The impact of target fragments on the biological dose is maximized in the entrance channels and for small α/β values. When applied to the description of survival data, model predictions including all fragments allowed better agreement to experimental data at high energies, while a minor effect was observed in the peak region. An improved description was also obtained when including the fragments’ contribution to describe RBE data. Overall, this analysis indicates that a minor contribution can be expected to the overall RBE resulting from target fragments. However, considering the fragmentation effects can improve the agreement with experimental data for high energy proton beams.

How LEM-based RBE and dose-averaged LET affected clinical outcomes of sacral chordoma patients treated with carbon ion radiotherapy

PURPOSE/OBJECTIVE

To understand the role of relative biological effectiveness (RBE) and dose-averaged linear energy transfer (LET_d) distributions in the treatment of sacral chordoma (SC) patients with carbon ion radiotherapy (CIRT).

MATERIAL/METHODS

Clinical plans of 50 SC patients consecutively treated before August 2018 with a local effect model-based optimization were recalculated with the modified microdosimetric kinetic RBE model (mMKM). Twenty-six patients were classified as progressive disease and the relapse volume was contoured on the corresponding follow-up diagnostic sequence. The remaining 24 patients populated the control group. Target prescription dose (D_RBE|50%), near-to-minimum- (D_RBE|95%) and near-to-maximum- (D_RBE|2%) doses were compared between the two cohorts in both RBE systems. LET_d distribution was evaluated for in-field relapsed cases with respect to the control group.

RESULTS

Target D_MKM|50% and D_MKM|95% were respectively 10% and 18% lower than what we aimed at. Dosimetric evaluators showed no significant difference, in neither of the RBE frameworks, between relapsed and control sets. Half of the relapse volumes were located in a well-covered high dose region. On average, over these cases, median target LET_d was significantly lower than the control cohort mean value (27 vs 30 keV/μm). Most notably, the volume receiving dose from high-LET particles (>50 keV/μm) lay substantially below recently reported data in the literature.

CONCLUSION

A combined multi model RBE- and LET-based optimization could play a key role in the enhancement of the therapeutic ratio of CIRT for large radioresistant tumors such as sacral chordomas.

Interaction of therapeutic12C ions with bone-like targets: physical characterization and dosimetric effect at material interfaces

Carbon therapy is a promising treatment option for cancer. The physical and biological properties of carbon ions can theoretically allow for the delivery of curative doses to the tumor, while simultaneously limiting risks of toxicity to adjacent healthy structures. The treatment effectiveness can be further improved by decreasing the uncertainties stemming from several sources, including the modeling of tissue heterogeneity. Current treatment plans employ density-based conversion methods to translate patient-specific anatomy into a water system, where dose distribution is calculated. This approach neglects differences in nuclear interactions stemming from the elemental composition of each tissue. In this work, we investigated the interaction of therapeutic carbon ions with bone-like materials. The study concentrated on nuclear interactions and included attenuation curves of 200 and 400 AMeV beams in different types of bones, as well as kinetic energy spectra of all charged fragments produced up to 29 degrees from the beam direction. The comparison between measurements and calculations of the treatment planning system TRiP98 indicated that bone tissue causes less fragmentation of carbon ions than water. Overall, hydrogen and helium particles were found to be the most abundant species, while heavier fragments were mostly detected within 5 degrees from the beam direction. We also investigated how the presence of a soft tissue-bone interface could affect the depth-dose profile. The results revealed a dose spike in the transition region, that extended from the entry channel to the target volume. The findings of this work indicated that the tissue-to-water conversion method based only on density considerations can result in dose inaccuracies. Tissue heterogeneity regions containing bones can potentially produce dose spikes, whose magnitude will depend on the patient anatomy. Dose uncertainties can be decreased by modeling nuclear interactions directly in bones, without applying the tissue-to-water conversion.

Carbon ion radiotherapy as definitive treatment in locally recurrent pancreatic cancer

Purpose

Data on management of locally recurrent pancreatic cancer (LRPC) after primary resection are limited. Recently, surprisingly high overall survival rates were reported after irradiation with carbon ions. Here, we report on our clinical experience using carbon ion radiotherapy as definitive treatment in LRPC at the Heidelberg Ion-Beam Therapy Center (HIT).

Methods

Between 2015 and 2019, we treated 13 patients with LRPC with carbon ions with a median total dose of 48 Gy (RBE) in 12 fractions using an active raster-scanning technique at a rotating gantry. No concomitant chemotherapy was administered. Overall survival, local control, and toxicity rates were evaluated 18 months after the last patient finished radiotherapy.

Results

With a median follow-up time of 9.5 months, one patient is still alive (8%). Median OS was 12.7 months. Ten patients (77%) developed distant metastases. Additionally, one local recurrence (8%) and two regional tumor recurrences (15%) were observed. The estimated 1‑year local control and locoregional control rates were 87.5% and 75%, respectively. During radiotherapy, we registered one gastrointestinal bleeding CTCAE grade III (8%) due to gastritis. The bleeding was sufficiently managed with conservative therapy. No further higher-grade acute or late toxicities were observed.

Conclusion

We demonstrate high local control rates in a rare cohort of LRPC patients treated with carbon ion radiotherapy. The observed median overall survival rate was not improved compared to historical in-house data using photon radiotherapy. This is likely due to a high rate of distant tumor progression, highlighting the necessity of additional chemotherapy.

FLASH radiotherapy with carbon ion beams

FLASH radiotherapy is considered a new potential breakthrough in cancer treatment. Ultra-high dose rates (>40 Gy/s) have been shown to reduce toxicity in the normal tissue without compromising tumor control, resulting in a widened therapeutic window. These high dose rates are more easily achievable in the clinic with charged particles, and clinical trials are, indeed, ongoing using electrons or protons. FLASH could be an attractive solution also for heavier ions such as carbon and could even enhance the therapeutic window. However, it is not yet known whether the FLASH effect will be the same as for sparsely ionizing radiation when densely ionizing carbons ions are used. Here we discuss the technical challenges in beam delivery and present a promising solution using 3D range-modulators in order to apply ultra-high dose rates (UHDR) compatible with FLASH with carbon ions. Furthermore, we will discuss the possible outcome of C-ion therapy at UHDR on the level of the radiobiological and radiation chemical effects.

Carbon Ion Dose Constraints in the Head and Neck and Skull Base: Review of MedAustron Institutional Protocols

Background

Dose constraints are of paramount importance for the outcome of any radiotherapy treatment. In this article, we report dose-volume constraints as well as currently used fractionation schedules for carbon ion radiotherapy as applied in MedAustron (Wiener Neustadt, Austria).

Materials and Methods

For fractionation schedules, both German and Japanese regimes were used. From the clinical experience of National Institute of Radiological Sciences (Chiba, Japan) and Heidelberg Ion Therapy (Heidelberg, Germany; formerly GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany) and the work by colleagues in Centro Nazionale Adroterapia Oncologica (Pavia, Italy) recalculating the dose from the microdosimetric kinetic model to the local effect model, we have set the dose constraints for critical organs of the head and neck area. Where no clinical data was available, an educated guess was made, based on data available from photon and proton series.

Results

We report the constraints for the optic nerve and chiasm, brainstem, spinal cord, cochlea, brain parenchyma, salivary gland, eye and adnexa, and mandibular/maxillary bone; constraints are grouped based on a fractionation scheme (German versus Japanese) and the risk of toxicity (safe, low to middle, and middle to high).

Conclusion

We think validation of dose constraints should present a relevant part of the activity of any carbon ion radiotherapy facility, and we anticipate future multicentric, joint evaluations.

Joint Optimization of Photon-Carbon Ion Treatments for Glioblastoma

PURPOSE

Carbon ions are radiobiologically more effective than photons and are beneficial for treating radioresistant gross tumor volumes (GTV). However, owing to a reduced fractionation effect, they may be disadvantageous for treating infiltrative tumors, in which healthy tissue inside the clinical target volume (CTV) must be protected through fractionation. This work addresses the question: What is the ideal combined photon-carbon ion fluence distribution for treating infiltrative tumors given a specific fraction allocation between photons and carbon ions?

METHODS AND MATERIALS

We present a method to simultaneously optimize sequentially delivered intensity modulated photon (IMRT) and carbon ion (CIRT) treatments based on cumulative biological effect, incorporating both the variable relative biological effect of carbon ions and the fractionation effect within the linear quadratic model. The method is demonstrated for 6 glioblastoma patients in comparison with the current clinical standard of independently optimized CIRT-IMRT plans.

RESULTS

Compared with the reference plan, joint optimization strategies yield inhomogeneous photon and carbon ion dose distributions that cumulatively deliver a homogeneous biological effect distribution. In the optimal distributions, the dose to CTV is mostly delivered by photons and carbon ions are restricted to the GTV with variations depending on tumor size and location. Improvements in conformity of high-dose regions are reflected by a mean EQD2 reduction of 3.29 ± 1.22 Gy in a dose fall-off margin around the CTV. Carbon ions may deliver higher doses to the center of the GTV, and photon contributions are increased at interfaces with CTV and critical structures. This results in a mean EQD2 reduction of 8.3 ± 2.28 Gy, in which the brain stem abuts the target volumes.

CONCLUSIONS

We have developed a biophysical model to optimize combined photon-carbon ion treatments. For 6 glioblastoma patient cases, we show that our approach results in a more targeted application of carbon ions that (1) reduces dose in normal tissues within the target volume, which can only be protected through fractionation; and (2) boosts central target volume regions to reduce integral dose. Joint optimization of IMRT-CIRT treatments enable the exploration of a new spectrum of plans that can better address physical and radiobiological treatment planning challenges.

High-dose hypofractionated pencil beam scanning carbon ion radiotherapy for lung tumors: Dosimetric impact of different spot sizes and robustness to interfractional uncertainties

PURPOSE

The robustness against setup and motion uncertainties of gated four-dimensional restricted robust optimization (4DRRO) was investigated for hypofractionated carbon ion radiotherapy (CIRT) of lung tumors.

METHODS

CIRT plans of 9 patients were optimized using 4DRRO strategy with 3 mm setup errors, 3% density errors and 3 breathing phases related to the gate window. The prescription was 60 Gy(RBE) in 4 fractions. Standard spots (SS) were compared to big spots (BS). Plans were recalculated on multiple 4DCTs acquired within 3 weeks from treatment simulation and rigidly registered with planning images using bone matching. Warped dose distributions were generated using deformable image registration and accumulated on the planning 4DCTs. Target coverage (D98%, D95% and V95%) and dose to lung were evaluated in the recalculated and accumulated dose distributions.

RESULTS

Comparable target coverage was obtained with both spot sizes (p = 0.53 for D95%). The mean lung dose increased of 0.6 Gy(RBE) with BS (p = 0.0078), still respecting the dose constraint of a 4-fraction stereotactic treatment for the risk of radiation pneumonitis. Statistically significant differences were found in the recalculated and accumulated D95% (p = 0.048 and p = 0.024), with BS showing to be more robust. Using BS, the average degradations of the D98%, D95% and V95% in the accumulated doses were -2.7%, -1.6% and -1.5%.

CONCLUSIONS

Gated 4DRRO was highly robust against setup and motion uncertainties. BS increased the dose to healthy tissues but were more robust than SS. The selected optimization settings guaranteed adequate target coverage during the simulated treatment course with acceptable risk of toxicity.

Mechanistic Modelling of Slow and Fast NHEJ DNA Repair Pathways Following Radiation for G0/G1 Normal Tissue Cells

Mechanistic in silico models can provide insight into biological mechanisms and highlight uncertainties for experimental investigation. Radiation-induced double-strand breaks (DSBs) are known to be toxic lesions if not repaired correctly. Non-homologous end joining (NHEJ) is the major DSB-repair pathway available throughout the cell cycle and, recently, has been hypothesised to consist of a fast and slow component in G0/G1. The slow component has been shown to be resection-dependent, requiring the nuclease Artemis to function. However, the pathway is not yet fully understood. This study compares two hypothesised models, simulating the action of individual repair proteins on DSB ends in a step-by-step manner, enabling the modelling of both wild-type and protein-deficient cell systems. Performance is benchmarked against experimental data from 21 cell lines and 18 radiation qualities. A model where resection-dependent and independent pathways are entirely separated can only reproduce experimental repair kinetics with additional restraints on end motion and protein recruitment. However, a model where the pathways are entwined was found to effectively fit without needing additional mechanisms. It has been shown that DaMaRiS is a useful tool when analysing the connections between resection-dependent and independent NHEJ repair pathways and robustly matches with experimental results from several sources.

Proton and carbon ion radiotherapy in skull base chordomas: a prospective study based on a dual particle and a patient-customized treatment strategy

BACKGROUND

The aim of this study is to evaluate results in terms of local control (LC), overall survival (OS), and toxicity profile and to better identify factors influencing clinical outcome of skull base chordoma treated with proton therapy (PT) and carbon ion radiotherapy (CIRT).

METHODS

We prospectively collected and analyzed data of 135 patients treated between November 2011 and December 2018. Total prescription dose in the PT group (70 patients) and CIRT group (65 patients) was 74 Gy relative biological effectiveness (RBE) delivered in 37 fractions and 70.4 Gy(RBE) delivered in 16 fractions, respectively (CIRT in unfavorable patients). LC and OS were evaluated using the Kaplan-Meier method. Univariate and multivariate analyses were performed, to identify prognostic factors on clinical outcomes.

RESULTS

After a median follow-up of 44 (range, 6-87) months, 14 (21%) and 8 (11%) local failures were observed in CIRT and PT group, respectively. Five-year LC rate was 71% in CIRT cohort and 84% in PT cohort. The estimated 5-year OS rate in the CIRT and PT group was 82% and 83%, respectively. On multivariate analysis, gross tumor volume (GTV), optic pathways, and/or brainstem compression and dose coverage are independent prognostic factors of local failure risk. High rate toxicity grade ≥3 was reported in 11% of patients.

CONCLUSIONS

Particle radiotherapy is an effective treatment for skull base chordoma with acceptable late toxicity. GTV, optic pathways, and/or brainstem compression and target coverage were independent prognostic factors for LC.

KEY POINTS

• Proton and carbon ion therapy are effective and safe in skull base chordoma.• Prognostic factors are GTV, organs at risk compression, and dose coverage.• Dual particle therapy and customized strategy was adopted.

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.