Head and neck IMPT probabilistic dose accumulation: Feasibility of a 2 mm setup uncertainty setting

Range uncertainty reductions in proton therapy and resulting improvements in quality-adjusted life expectancy (QALE) for head-and-neck cancer patients

Objective. Due to higher dose conformality to the target, proton radiotherapy for cancer has received rapidly-growing interest. However, uncertainties in thein vivoproton range and methods to reduce them remain active areas of research. Based on 20 patients with head-and-neck cancer, this study aims to quantify the benefits of proton range uncertainty reductions in terms of the resulting improvements in quality-adjusted life expectancy (QALE).Approach. For each patient, two different proton therapy treatment plans were created, which assumed a current clinical range uncertainty of approximately 3.5% (IMPT3.5%) and a potentially achievable range uncertainty of 1.0% (IMPT1%). A Markov model considering the probability of tumor control and the development of xerostomia, larynx edema, secondary cancer, and/or metastases as well as death from primary cancer, secondary cancer, metastases, or unrelated causes was constructed, and for every patient and treatment plan, 10 000 simulations of the patient's entire lifetime from the time of treatment until death were performed.Main results.A 3.5%-1% range uncertainty reduction increased QALE by up to 0.4 quality-adjusted life years (QALYs) in the nominal and up to 0.6 QALY in the worst-case scenario, equivalent to 4.8 months and 7.2 months of life in perfect health. This was largely the result of a reduction in healthy tissue toxicity rates, which were reduced by up to 8.5 percentage points (pp) and 10.0 pp in the nominal and worst-case scenario, respectively.Significance. The benefits of a 3.5%-1% range uncertainty reduction in 20 patients with head-and-neck cancer were quantified in terms of the associated improvement in QALE. The highest QALE improvements were observed in patients in the top quartile of youngest patients at the time of treatment, due to the longer potential lifespan over which prevented healthy tissue toxicities would have impacted the patients' quality of life.

Instantaneous in vivo distal edge verification in intensity-modulated proton therapy by means of PET imaging

BACKGROUND

Intensity-modulated proton therapy (IMPT) holds promise for improving outcomes in head-and-neck cancer (HNC) patients by enhancing organ-at-risk (OAR) sparing. A key challenge in IMPT is ensuring an accurate dose delivery at the distal edge of the tumor, where the steep dose gradients make treatment precision highly sensitive to uncertainties in both proton range and patient setup. Thus, IMPT conformality is increased by incorporating robust margins in the treatment optimization. However, an increment in the plan robustness could lead to an OAR overdosing. Therefore, an accurate distal edge verification during dose delivery is crucial to increase IMPT conformality by reducing optimization settings in treatment planning.

PURPOSE

This work aims to evaluate, in a quasi-clinical setting, a novel approach for accurate instantaneous proton beam distal edge verification in IMPT by means of spot-by-spot positron emission tomography (PET) imaging.

METHODS

An anthropomorphic head and neck phantom CIRS-731 HN was irradiated at the head and neck region. The targets were defined as 4 cm diameter spheres. A 60-ms delay was introduced between the proton beam spots in order to enable the spot-by-spot coincidence detection of the 511-keV photons resulting from positron annihilation following the positron emission from very short-lived positron-emitting, mainly 12N (T1/2  = 11.0 ms). Additionally, modified irradiations were carried out using solid water slabs of 2 and 5 mm thickness in the beam path to assess the precision of the approach for detecting range deviations. The positron activity range (PAR) was determined from the 50% distal fall-off position of the 1D longitudinal positron activity profile derived from the 2D image reconstructions. Furthermore, Monte Carlo (MC) simulations were performed using an in-house RayStation/GATE MC framework to predict the positron activity images and verify the PAR measurements.

RESULTS

PAR measurements achieved a precision between 1.5 and 3.6 mm (at 1.5σ clinical level) at the beam spot level within sub-second time scales. Measured PAR shifts of 1.6-2.1  and 4.2--.7 mm were observed with the 2- and 5-mm thickness range shifters, respectively, aligning with the corresponding proton dose range (PDR) shifts of 1.3-1.8 and 3.9-4.3 mm. The simulated PAR agrees with the measured PARs, showing an average range difference of ∼0.4 mm.

CONCLUSION

This study demonstrated the feasibility of instantaneous distal edge verification using PET imaging by introducing beam spot delays during dose delivery. The findings represent a first step toward the clinical implementation of instantaneous in vivo distal edge verification. The approach contributes to the development of real-time range verification aimed at improving IMPT treatments by mitigating range and setup uncertainties, thereby reducing dose to organs-at-risk and ultimately enhancing patient outcomes.

Steenbakkers R, Vogel WV, Both S. Relating proton LETd to biological response of parotid and submandibular glands using PSMA-PET in clinical patients

Background and purpose

A recent study investigated the use of PSMA-PET in monitoring loss of secretory cells in salivary glands of head and neck cancer (HNC) patients. Previously, a dose-effect relation has been formulated to the PSMA-PET uptake in salivary glands. The aim of this study was to derive a proton RBE model from the PSMA-PET uptake in salivary glands after proton therapy of HNC patients.

Materials and methods

Six patients treated with proton therapy were included. These patients received a PET-CT scan using 68Ga (N = 1) or 18F (N = 5) PSMA before treatment (baseline) and one month after the last fraction (follow-up). Physical dose (D), D·LETd and the follow-up PSMA-PET scan were deformed to the baseline PET-CT using deformable image registration. Parotid and submandibular gland delineations were adjusted to include voxels which had an uptake of ≥ 5 g/ml in the baseline PSMA-PET scan.

Results

The average RBE-LET slope was 0.075 [0.009; 0.125] (keV/μm)-1 (mean [95 %CI]) for parotid and submandibular glands combined. When analyzing parotid or submandibular glands separately the RBE-LET curve slope varies with two and five patients showing a positive RBE-LET slope when only analyzing parotid or submandibular glands respectively.

Conclusion

Our study did not find clear evidence of an increased RBE in parotid and submandibular glands with increasing LETd. On average an LETd effect was observed, however our sample size was too small to clearly define an RBE-LET relation. A larger cohort scanned at later time intervals could shed more light on this issue.

Linear approximation of variable relative biological effectiveness models for proton therapy

The McNamara (MCN) and Wedenberg (WED) RBE weighted dose (DRBE), dose and dose-weighted average LET (LETd) were calculated in twenty brain cancer patients. A linear approximation was made for each RBE model to give best agreement to clinically relevant dosimetric parameters. Additional evaluations were done on twenty head and neck and twenty breast cancer patients.The R2 of the fits was ≥0.94 and ≥0.91 for MCN and WED respectively for α/β values ≥1.0 Gy. The graphs derived in this work can be used to convert RBE-LET slopes derived from clinical data to α/β values in the MCN or WED models.

Margin and robustness settings for a library-of-plans IMPT strategy for locally advanced cervical cancer

Objective.This study aims to determine a margin and robustness setting for treating locally advanced cervical cancer (LACC) with a library-of-plans (LoP) based online-adaptive intensity-modulated proton therapy (IMPT).Approach.We analyzed 13 LACC patients with delineated planning and weekly repeat CT scans (reCTs). For each patient, 120 IMPT treatments of 25 fractions were simulated with a LoPs approach. Six different robustness settings (2-7 mm set-up robustness (SR) plus 3% range robustness (RR)) were used to create those 120 IMPT plans. Each fraction was simulated with a weekly reCT, combined with the sampling of inter- and intrafraction treatment uncertainties. The fraction doses were accumulated to obtain a treatment dose to the target volumes, distinguishing between the low-risk clinical target volume (CTV-T-LR) and the elective CTV (CTV-E). If one of the two targets obtained an adequate coverage for more than 90% of the treatments, different anisotropic margins were sampled on top of the robustness setting to the other target to obtain the Pareto-optimal margin in terms of adequate coverage versus increase in target volume.Main results.The percentage of treatments that reach the dose criterionV42.75Gy> 95% for the CTV-T-LR was 22.3%, 28.5%, 51.2%, 73.1%, 85.3%, and 90.0% for 2, 3, 4, 5, 6, and 7 mm SR plus 3% RR and for the CTV-E, this percentage was 60.4%, 73.8%, 86.5%, 92.3%, 96.9%, and 98.5%. The Pareto-optimal margin combined with a 5 mm/3% robustness setting for the CTV-T-LR with an adequate coverage for >90% of the treatments was given by {0, 1, 0, 3, 3, 0} mm in the left, right, anterior, posterior, cranial, caudal direction.Significance.Our study evaluated combinations of robustness and anisotropic margin settings for IMPT for LACC. With 5 mm SR and 3% RR for CTV-E and CTV-T-LR plus a margin to the CTV-T-LR of {0, 1, 0, 3, 3, 0} mm in left, right, anterior, posterior, cranial, and caudal ensured an adequate coverage for >90% of the simulated IMPT treatments.

A probabilistic evaluation of the Dutch robustness and model-based selection protocols for Head-and-Neck IMPT: A multi-institutional study

BACKGROUND AND PURPOSE

In the Netherlands, 2 protocols have been standardized for PT among the 3 proton centers: a robustness evaluation (RE) to ensure adequate CTV dose and a model-based selection (MBS) approach for IMPT patient-selection. This multi-institutional study investigates (i) inter-patient and inter-center variation of target dose from the RE protocol and (ii) the robustness of the MBS protocol against treatment errors for a cohort of head-and-neck cancer (HNC) patients treated in the 3 Dutch proton centers.

MATERIALS AND METHODS

Clinical treatment plans of 100 HNC patients were evaluated. Polynomial Chaos Expansion (PCE) was used to perform a comprehensive robustness evaluation per plan, enabling the probabilistic evaluation of 100,000 complete fractionated treatments. PCE allowed to derive scenario distributions of clinically relevant dosimetric parameters to assess CTV dose (D99.8%/D0.2%, based on a prior photon plan calibration) and tumour control probabilities (TCP) as well as the evaluation of the dose to OARs and normal tissue complication probabilities (NTCP) per center.

RESULTS

For the CTV70.00, doses from the RE protocol were consistent with the clinical plan evaluation metrics used in the 3 centers. For the CTV54.25, D99.8% were consistent with the clinical plan evaluation metrics at center 1 and 2 while, for center 3, a reduction of 1 GyRBE was found on average. This difference did not impact modelled TCP at center 3. Differences between expected and nominal NTCP were below 0.3 percentage point for most patients.

CONCLUSION

The standardization of the RE and MBS protocol lead to comparable results in terms of TCP and the NTCPs. Still, significant inter-patient and inter-center variation in dosimetric parameters remained due to clinical practice differences at each institution. The MBS approach is a robust protocol to qualify patients for PT.

RayStation/GATE Monte Carlo simulation framework for verification of proton therapy based on the12N imaging

Objective.12N, having a half-life of 11 ms, is a highly effective positron emitter that can potentially provide near real-time feedback in proton therapy. There is currently no framework for comparing and validating positron emission imaging of12N. This work describes the development and validation of a Monte Carlo (MC) framework to calculate the images of12N, as well as long-lived isotopes, originating from activation by protons.Approach. The available dual-panel Biograph mCT PET scanner was modeled in GATE and validated by comparing the simulated sensitivity map with the measured one. The distributions of12N and long-lived isotopes were calculated by RayStation and used as the input of GATE simulations. The RayStation/GATE combination was verified using proton beam irradiations of homogeneous phantoms. A 120 MeV pulsed pencil beam with 108protons per pulse was used. Two-dimensional images were created from the GATE output and compared with the images based on the measurements and the 1D longitudinal projection of the full 2D image was used to calculate the12N activity range.Main results. The simulated sensitivity in the center of the FoV (5.44%) agrees well with the measured one (5.41%). The simulated and measured 2D sensitivity maps agree in good detail. The relative difference between the measured and simulated positron activity range for both12N and long-lived isotopes is less than 1%. The broadening of the12N images relative to those of the longer-lived isotopes can be understood in terms of the large positron range of12N.Significance. We developed and validated a MC framework based on RayStation/GATE to support the in-beam PET method for quality assurance of proton therapy. The inclusion of the very short-lived isotope12N makes the framework useful for developing near real-time verification. This represents a significant step towards translating12N real-time in vivo verification to the clinic.

Robustness evaluation of pencil beam scanning proton therapy treatment planning: A systematic review

Compared to conventional radiotherapy using X-rays, proton therapy, in principle, allows better conformity of the dose distribution to target volumes, at the cost of greater sensitivity to physical, anatomical, and positioning uncertainties. Robust planning, both in terms of plan optimization and evaluation, has gained high visibility in publications on the subject and is part of clinical practice in many centers. However, there is currently no consensus on the methods and parameters to be used for robust optimization or robustness evaluation. We propose to overcome this deficiency by following the modified Delphi consensus method. This method first requires a systematic review of the literature. We performed this review using the PubMed and Web Of Science databases, via two different experts. Potential conflicts were resolved by a third expert. We then explored the different methods before focusing on clinical studies that evaluate robustness on a significant number of patients. Many robustness assessment methods are proposed in the literature. Some are more successful than others and their implementation varies between centers. Moreover, they are not all statistically or mathematically equivalent. The most sophisticated and rigorous methods have seen more limited application due to the difficulty of their implementation and their lack of widespread availability.

Target coverage and organs at risk dose in hypofractionated salvage radiotherapy after prostatectomy

Background and purpose

Introducing moderately hypofractionated salvage radiotherapy (SRT) following prostatectomy obligates investigation of its effects on clinical target volume (CTV) coverage and organ-at-risk (OAR) doses. This study assessed interfractional volume and dose changes in OARs and CTV in moderately hypofractionated SRT and evaluated the 8-mm planning target volume (PTV) margin.

Materials and methods

Twenty patients from the PERYTON-trial were included; 10 received conventional SRT (35 × 2 Gy) and 10 hypofractionated SRT (20 × 3 Gy). OARs were delineated on 539 pre-treatment Cone Beam CT (CBCT) scans to compare interfractional OAR volume changes. CTVs for the hypofractionated group were delineated on 199 CBCTs. Dose distributions with 4 and 6 mm PTV margins were generated using voxel-wise minimum robustness evaluation of the original 8-mm PTV plan, and dose changes were assessed.

Results

Median volume changes for bladder and rectum were -26 % and -10 %, respectively. OAR volume changes were not significantly different between the two treatment schedules. The 8-mm PTV margin ensured optimal coverage for prostate bed and vesicle bed CTV (V95 = 100 % in >97 % fractions). However, bladder V60 <25 % was not achieved in 5 % of fractions, and rectum V60 <5 % was unmet in 33 % of fractions. A 6-mm PTV margin resulted in CTV V95 = 100 % in 92 % of fractions for prostate bed, and in 86 % for vesicle bed CTV.

Conclusions

Moderately hypofractionated SRT yielded comparable OAR volume changes to conventionally fractionated SRT. Interfractional changes remained acceptable with a PTV margin of 6 mm for prostate bed and 8 mm for vesicle bed.

Benefit of range uncertainty reduction in robust optimisation for proton therapy of brain, head-and-neck and breast cancer patients

Background and Purpose

The primary cause of range uncertainty in proton therapy is inaccuracy in estimating the stopping-power ratio from computed tomography. This study examined the impact on dose-volume metrics by reducing range uncertainty in robust optimisation for a diverse patient cohort and determined the level of range uncertainty that resulted in a relevant reduction in doses to organs-at-risk (OARs).

Materials and Methods

The effect of reducing range uncertainty on OAR doses was evaluated by robustly optimising six proton plans with varying range uncertainty levels (ranging from 3.5% in the original plan to 1.0%), keeping setup uncertainty fixed. All plans used the initial clinical treatment plan's beam directions and optimisation objectives and were optimised until a clinically acceptable plan was achieved across all setup and range scenarios. The effect of reduced range uncertainty on dose-volume metrics for OARs near the target was evaluated. This study included 30 brain cancer patients, as well as five head-and-neck and five breast cancer patients, investigating the relevance of reducing range uncertainty when different setup uncertainties were used.

Results

Lowering range uncertainty slightly reduced the nominal dose to surrounding tissue. For body volume receiving 80% of the prescribed dose, reducing range uncertainty from 3.5% to 2.0% resulted in a median decrease of 4 cm3 for the brain, 17 cm3 for head-and-neck, and 27 cm3 for breast cancer patients.

Conclusions

Reducing range uncertainty in robust optimisation showed a reduction in dose to OARs. The clinical relevance depends on the affected organs and the clinical dose constraints.

The first investigation of the dosimetric perturbations from the spot position errors in spot-scanning arc therapy (SPArc)

Objective. To quantitatively investigate the impact of spot position error (PE) on the dose distribution in (Spot-scanning arc therapy) SPArc plans compared to Intensity-Modulated Proton Therapy (IMPT).Approach.Twelve representative cases, including brain, lung, liver, and prostate cancers, were retrospectively selected. Spot PEs were simulated during dynamic SPArc treatment delivery. Two types of errors were generated, including random error and systematic error. Two different probability distributions of random errors were used (1) Gaussian distribution (PEran-GS) (2) uniform distribution (PEran-UN). In PEran-UN, four sub-scenarios were considered: 25%, 50%, 75%, and 100% spots were randomly selected in various directions on the scale of 0-1 mm or 0-2 mm of PE. Additionally, systematic error was simulated by shifting all the spot uniformly by 1 or 2 mm in various directions (PEsys). Gamma-index Passing Rate (GPR) is applied to assess the dosimetric perturbation quantitatively.Main results.For PEran-GSin the 1 mm scenario, both SPArc and IMPT are comparable with a GPR exceeding 99%. However, for PEran-GSin 2 mm scenario, SPArc could provide better GPR. As PEsysof 2 mm, SPArc plans have a much better GPR compared to IMPT plans: SPArc's GPR is 99.59 ± 0.47%, 93.82 ± 4.07% and 64.58 ± 15.83% for 3 mm/3%, 2 mm/2% and 1 mm/1% criteria compared to IMPT with 97.49 ± 2.44%, 84.59 ± 4.99% and 42.02 ± 6.31%.Significance.Compared to IMPT, SPArc shows better dosimetric robustness in spot PEs. This study presents the first simulation results and the methodology that serves as a reference to guide future investigations into the accuracy and quality assurance of SPArc treatment delivery.

Optimal timing of re-planning for head and neck adaptive radiotherapy

BACKGROUND AND PURPOSE

Adaptive radiotherapy (ART) relies on re-planning to correct treatment variations, but the optimal timing of re-planning to account for dose changes in head and neck organs at risk (OARs) is still under investigation. We aimed to find out the optimal timing of re-planning in head and neck ART.

MATERIALS AND METHODS

A total of 110 head and neck cancer patients were retrospectively enrolled. A semi auto-segmentation method was applied to obtain the weekly mean dose (Dmean) to OARs. The K-nearest-neighbour method was used for missing data imputation of weekly Dmean. A dose deviation map was built using the planning Dmean and weekly Dmean values and then used to simulate different ART scenarios consisting of 1 to 6 re-plannings. The difference between accumulated Dmean and planning Dmean before re-planning (ΔDmean_acc_noART) and after re-planning (ΔDmean_acc_ART) were evaluated and compared.

RESULTS

Among all the OARs, supraglottic showed the largest ΔDmean_acc_noART (1.23 ± 3.13 Gy) and most cases of ΔDmean_acc_noART > 3 Gy (26 patients). The 3rd week is suggested in the optimal timing of re-planning for 10 OARs. For all the organs except arytenoid, 2 re-plannings were able to guarantee the ΔDmean_acc_ART below 3 Gy while the average |ΔDmean_acc_ART| was below 1 Gy. ART scenarios of 2_4, 3_4, 3_5 (week of re-planning separated with "_") were able to guarantee ΔDmean_acc_ART of 99 % of patients below 3 Gy simultaneously for 19 OARs.

CONCLUSIONS

The optimal timing of re-planning was suggested for different organs at risk in head and neck adaptive radiotherapy. Generic scenarios of timing and frequency for re-planning can be applied to guarantee the increase of accumulated mean dose within 3 Gy simultaneously for multiple organs.

Commissioning and Validation of CT Number to SPR Calibration in Carbon Ion Therapy Facility

Purpose

We performed computed tomography (CT)-stopping power ratio (SPR) calibration in a carbon-ion therapy facility and evaluated SPR estimation accuracy.

Materials and Methods

A polybinary tissue model method was used for the calibration of CT numbers and SPR. As a verification by dose calculation, we created a virtual phantom to which the CT-SPR calibration table was applied. Then, SPR was calculated from the change in the range of the treatment planning beam when changing to 19 different CT numbers, and the accuracy of the treatment planning system (TPS) calculation of SPR values from the CT-SPR calibration table was validated. As a verification by measurement, 5 materials (water, milk, olive oil, ethanol, 40% K2HPO4) were placed in a container, and the SPR was obtained by measurement from the change in the range of the beam that passed through the materials.

Results

The results of the dose calculations of the TPS showed that the results agreed within 1% for the lower CT numbers up to 1000 HU, but there was a difference of 3.0% in the higher CT number volume. The difference between the SPR calculated by TPS and the SPR caused by the difference in the energy of the incident particles agreed within 0.51%. The accuracy of SPR estimation was measured, and the error was within 2% for all materials tested.

Conclusion

These results indicate that the SPR estimation errors are within the range of errors that can be expected in particle therapy. From commissioning and verification results, the CT-SPR calibration table obtained during this commissioning process is clinically applicable.

Impact of and interplay between proton arc therapy and range uncertainties in proton therapy for head-and-neck cancer

Objective. Proton therapy reduces the integral dose to the patient compared to conventional photon treatments. However,in vivoproton range uncertainties remain a considerable hurdle. Range uncertainty reduction benefits depend on clinical practices. During intensity-modulated proton therapy (IMPT), the target is irradiated from only a few directions, but proton arc therapy (PAT), for which the target is irradiated from dozens of angles, may see clinical implementation by the time considerable range uncertainty reductions are achieved. It is therefore crucial to determine the impact of PAT on range uncertainty reduction benefits.Approach. For twenty head-and-neck cancer patients, four different treatment plans were created: an IMPT and a PAT treatment plan assuming current clinical range uncertainties of 3.5% (IMPT3.5%and PAT3.5%), and an IMPT and a PAT treatment plan assuming that range uncertainties can be reduced to 1% (IMPT1%and PAT1%). Plans were evaluated with respect to target coverage and organ-at-risk doses as well as normal tissue complication probabilities (NTCPs) for parotid glands (endpoint: parotid gland flow <25%) and larynx (endpoint: larynx edema).Main results. Implementation of PAT (IMPT3.5%-PAT3.5%) reduced mean NTCPs in the nominal and worst-case scenario by 3.2 percentage points (pp) and 4.2 pp, respectively. Reducing range uncertainties from 3.5% to 1% during use of IMPT (IMPT3.5%-IMPT1%) reduced evaluated NTCPs by 0.9 pp and 2.0 pp. Benefits of range uncertainty reductions subsequently to PAT implementation (PAT3.5%-PAT1%) were 0.2 pp and 1.0 pp, with considerably higher benefits in bilateral compared to unilateral cases.Significance. The mean clinical benefit of implementing PAT was more than twice as high as the benefit of a 3.5%-1% range uncertainty reduction. Range uncertainty reductions are expected to remain beneficial even after PAT implementation, especially in cases with target positions allowing for full leveraging of the higher number of gantry angles during PAT.

Assessment of residual geometrical errors of clinical target volumes and their impact on dose accumulation for head and neck radiotherapy

PURPOSE

To assess the residual geometrical errors (dr) and their impact on the clinical target volumes (CTV) dose coverage for head and neck cancer (HNC) proton therapy patients.

METHODS

We analysed 28 HNC patients treated with 70 Gy (RBE) and 54.25 Gy (RBE) to the therapeutic CTV70 and prophylactic CTV54.25, respectively. Daily cone beam CTs were converted to high quality synthetic CTs (sCTs). The CTVs from the nominal CT were propagated to the corresponding sCTs using a hybrid deformable image registration (propagated CTVs) in RayStation 11B. For 11 patients, all propagated CTVs were reviewed by our HNC radiation oncologist (physician corrected CTVs). The residual geometrical error dr was quantified as a function of the daily CTVs volume overlap with the nominal plan CTV. The errors dr(propagated CTVs) and dr(physician corrected CTVs) and the difference in dice similarity coefficients (ΔDSC) were determined. Using clinical plans, dose coverage and the tumor control probability (TCP) for the nominal, accumulated and voxel-wise minimum scenarios were determined.

RESULTS

The difference in the residual geometrical error dr (propagated CTVs - physician corrected CTVs) and mean DSC (|ΔDSC|mean) were minor: Δdr(CTV70) = 0.16 mm, Δdr(CTV54.25) = 0.26 mm, |ΔDSC|mean < 0.9%. For all 28 patients, dr(CTV70) = 1.91 mm and dr(CTV54.25) = 1.90 mm. However, CTV54.25 above and below the cricoid cartilage differed substantially (1.00 mm c.f. 3.93 mm). The CTV54.25 coverage below the cricoid was then almost always lower, although the TCP of the accumulated dose was higher than the TCP of the voxel-wise minimum dose.

CONCLUSIONS

Setup uncertainty setting of 2 mm is possible. The feasibility of using propagated CTVs for error determination is demonstrated.

Clinical benefit of range uncertainty reduction in proton treatment planning based on dual-energy CT for neuro-oncological patients

OBJECTIVE

Several studies have shown that dual-energy CT (DECT) can lead to improved accuracy for proton range estimation. This study investigated the clinical benefit of reduced range uncertainty, enabled by DECT, in robust optimisation for neuro-oncological patients.

METHODS

DECT scans for 27 neuro-oncological patients were included. Commercial software was applied to create stopping-power ratio (SPR) maps based on the DECT scan. Two plans were robustly optimised on the SPR map, keeping the beam and plan settings identical to the clinical plan. One plan was robustly optimised and evaluated with a range uncertainty of 3% (as used clinically; denoted 3%-plan); the second plan applied a range uncertainty of 2% (2%-plan). Both plans were clinical acceptable and optimal. The dose-volume histogram parameters were compared between the two plans. Two experienced neuro-radiation oncologists determined the relevant dose difference for each organ-at-risk (OAR). Moreover, the OAR toxicity levels were assessed.

RESULTS

For 24 patients, a dose reduction >0.5/1 Gy (relevant dose difference depending on the OAR) was seen in one or more OARs for the 2%-plan; e.g. for brainstem D0.03cc in 10 patients, and hippocampus D40% in 6 patients. Furthermore, 12 patients had a reduction in toxicity level for one or two OARs, showing a clear benefit for the patient.

CONCLUSION

Robust optimisation with reduced range uncertainty allows for reduction of OAR toxicity, providing a rationale for clinical implementation. Based on these results, we have clinically introduced DECT-based proton treatment planning for neuro-oncological patients, accompanied with a reduced range uncertainty of 2%.

ADVANCES IN KNOWLEDGE

This study shows the clinical benefit of range uncertainty reduction from 3% to 2% in robustly optimised proton plans. A dose reduction to one or more OARs was seen for 89% of the patients, and 44% of the patients had an expected toxicity level decrease.

Validation of robust radiobiological optimization algorithms based on the mixed beam model for intensity-modulated carbon-ion therapy

Currently, treatment planning systems (TPSs) that can compute the intensities of intensity-modulated carbon-ion therapy (IMCT) using scanned carbon-ion beams are limited. In the present study, the computational efficacy of the newly designed IMCT algorithms was analyzed for the first time based on the mixed beam model with respect to the physical and biological doses; moreover, the validity and effectiveness of the robust radiobiological optimization were verified. A dose calculation engine was independently generated to validate a clinical dose determined in the TPS. A biological assay was performed using the HSGc-C5 cell line to validate the calculated surviving fraction (SF). Both spot control (SC) and voxel-wise worst-case scenario (WC) algorithms were employed for robust radiobiological optimization followed by their application in a Radiation Therapy Oncology Group benchmark phantom under homogeneous and heterogeneous conditions and a clinical case for range and position errors. Importantly, for the first time, both SC and WC algorithms were implemented in the integrated TPS platform that can compute the intensities of IMCT using scanned carbon-ion beams for robust radiobiological optimization. For assessing the robustness, the difference between the maximum and minimum values of a dose-volume histogram index in the examined error scenarios was considered as a robustness index. The relative biological effectiveness (RBE) determined by the independent dose calculation engine exhibited a -0.6% difference compared with the RBE defined by the TPS at the isocenter, whereas the measured and the calculated SF were similar. Regardless of the objects, compared with the conventional IMCT, the robust radiobiological optimization enhanced the sensitivity of the examined error scenarios by up to 19% for the robustness index. The computational efficacy of the novel IMCT algorithms was verified according to the mixed beam model with respect to the physical and biological doses. The robust radiobiological optimizations lowered the impact of range and position uncertainties considerably in the examined scenarios. The robustness of the WC algorithm was more enhanced compared with that of the SC algorithm. Nevertheless, the SC algorithm can be used as an alternative to the WC IMCT algorithm with respect to the computational cost.

Proton arc therapy increases the benefit of proton therapy for oropharyngeal cancer patients in the model based clinic

BACKGROUND AND PURPOSE

In the model-based approach, patients qualify for proton therapy when the reduction in risk of toxicity (ΔNTCP) obtained with IMPT relative to VMAT is larger than predefined thresholds as defined by the Dutch National Indication Protocol (NIPP). Proton arc therapy (PAT) is an emerging technology which has the potential to further decrease NTCPs compared to IMPT. The aim of this study was to investigate the potential impact of PAT on the number of oropharyngeal cancer (OPC) patients that qualify for proton therapy.

MATERIALS AND METHODS

A prospective cohort of 223 OPC patients subjected to the model-based selection procedure was investigated. 33 (15%) patients were considered unsuitable for proton treatment before plan comparison. When IMPT was compared to VMAT for the remaining 190 patients, 148 (66%) patients qualified for protons and 42 (19%) patients did not. For these 42 patients treated with VMAT, robust PAT plans were generated.

RESULTS

PAT plans provided better or similar target coverage compared to IMPT plans. In the PAT plans, integral dose was significantly reduced by 18% relative to IMPT plans and by 54% relative to VMAT plans. PAT decreased the mean dose to numerous organs-at-risk (OARs), further reducing NTCPs. The ΔNTCP for PAT relative to VMAT passed the NIPP thresholds for 32 out of the 42 patients treated with VMAT, resulting in 180 patients (81%) of the complete cohort qualifying for protons.

CONCLUSION

PAT outperforms IMPT and VMAT, leading to a further reduction of NTCP-values and higher ΔNTCP-values, significantly increasing the percentage of OPC patients selected for proton therapy.

Evaluating the effect of setup uncertainty reduction and adaptation to geometric changes on normal tissue complication probability using online adaptive head and neck intensity modulated proton therapy

Objective. To evaluate the impact of setup uncertainty reduction (SUR) and adaptation to geometrical changes (AGC) on normal tissue complication probability (NTCP) when using online adaptive head and neck intensity modulated proton therapy (IMPT).Approach.A cohort of ten retrospective head and neck cancer patients with daily scatter corrected cone-beam CT (CBCT) was studied. For each patient, two IMPT treatment plans were created: one with a 3 mm setup uncertainty robustness setting and one with no explicit setup robustness. Both plans were recalculated on the daily CBCT considering three scenarios: the robust plan without adaptation, the non-robust plan without adaptation and the non-robust plan with daily online adaptation. Online-adaptation was simulated using an in-house developed workflow based on GPU-accelerated Monte Carlo dose calculation and partial spot-intensity re-optimization. Dose distributions associated with each scenario were accumulated on the planning CT, where NTCP models for six toxicities were applied. NTCP values from each scenario were intercompared to quantify the reduction in toxicity risk induced by SUR alone, AGC alone and SUR and AGC combined. Finally, a decision tree was implemented to assess the clinical significance of the toxicity reduction associated with each mechanism.Main results. For most patients, clinically meaningful NTCP reductions were only achieved when SUR and AGC were performed together. In these conditions, total reductions in NTCP of up to 30.48 pp were obtained, with noticeable NTCP reductions for aspiration, dysphagia and xerostomia (mean reductions of 8.25, 5.42 and 5.12 pp respectively). While SUR had a generally larger impact than AGC on NTCP reductions, SUR alone did not induce clinically meaningful toxicity reductions in any patient, compared to only one for AGC alone.SignificanceOnline adaptive head and neck proton therapy can only yield clinically significant reductions in the risk of long-term side effects when combining the benefits of SUR and AGC.

Automated Robust Planning for IMPT in Oropharyngeal Cancer Patients Using Machine Learning

PURPOSE

The aim of this study was to evaluate an automated treatment planning method for robustly optimized intensity modulated proton therapy (IMPT) plans for oropharyngeal carcinoma patients and to compare the results with manually optimized robust IMPT plans.

METHODS AND MATERIALS

An atlas regression forest-based machine learning (ML) model for dose prediction was trained on CT scans, contours, and dose distributions of robust IMPT plans of 88 oropharyngeal cancer (OPC) patients. The ML model was combined with a robust voxel and dose volume histogram-based dose mimicking optimization algorithm for 21 perturbed scenarios to generate a machine-deliverable plan from the predicted dose distribution. Machine learning optimization (MLO) configuration was performed using a cross-validation approach with 3 × 8 tuning patients and comprised of adjustments to the mimicking optimization, to generate higher-quality MLO plans. Independent testing of the MLO algorithm was performed with another 25 patients. Plan quality of clinical and MLO plans was evaluated by the clinical target volume (D98% voxel-wise minimum dose >94%), organ at risk (OAR) doses, and the normal tissue complication probability (NTCP) (sum (Σ) of grade-2 and grade-3 dysphagia and xerostomia).

RESULTS

Adequate robust target coverage was achieved in 24/25 clinical plans and in 23/25 MLO plans in the primary clinical target volume (CTV). In the elective CTV, 22/25 clinical plans and 24/25 MLO plans passed the robust target coverage evaluation threshold. The MLO average Σgrade 2 and Σgrade 3 NTCPs were comparable to the clinical plans (Σgrade 2 NTCPs: clinical 47.5% vs MLO 48.4%, Σgrade 3 NTCPs: clinical 11.9% vs MLO 12.3%). Significant increases in OAR average doses in MLO plans were found in the pharynx constrictor muscles (163 cGy, P = .002) and cervical esophagus (265 cGy, P = .002). The MLO plans were created within 45 minutes.

CONCLUSION

This study showed that automated MLO planning can generate robustly optimized MLO plans with quality comparable to clinical plans in OPC patients.

Spot scanning proton arc therapy reduces toxicity in oropharyngeal cancer patients

BACKGROUND

Proton arc technology has recently shown dosimetric gains for various treatment indications. The increased number of beams and energy layers (ELs) in proton arc plans, increases the degrees of freedom in plan optimization and thereby flexibility to spare dose in organs at risk (OARs). A relationship exists between dosimetric plan quality, delivery efficiency, the number of ELs -and beams in a proton arc plan.

PURPOSE

This work aims to investigate the effect of the number of beams and ELs in a proton arc plan, on toxicity and delivery time for oropharyngeal cancer patients (OPC) selected for intensity modulated proton therapy (IMPT) based on the Dutch model-based approach.

METHODS

The EL reduction algorithm iteratively selects ELs from beams equidistantly spaced over a 360° arc. The beams in the final plan may contain multiple ELs, making them suited for static delivery on the studied treatment machine. The produced plans can therefore be called "step and shoot" proton arc plans. The number of beams and ELs were varied to determine the relationship with the planning cost function value, normal tissue complication probability (NTCP) and delivery time. Proton arc plans with robust target coverage and optimal energy layer reduction (ELR) settings to reduce NTCP, were generated for 10 OPC patients. Proton arc plans were compared to clinical volumetric modulated arc therapy (VMAT) and IMPT plans in terms of integral dose, OAR dose, NTCP for xerostomia and dysphagia and delivery time. Furthermore, dose-weighted average linear energy transfer (LET_d ) distributions were compared between the IMPT and proton arc plans. A dry run delivery of a plan containing 20 beams and 360 ELs was performed to evaluate delivery time and accuracy.

RESULTS

We found 360 ELs distributed over 30 beams generated proton arc plans with near minimal expected plan toxicity. Relative to corresponding IMPT and VMAT plans, an average reduction of 21 ± 3% and 58 ± 10% in integral dose was observed. D m e a n $_{mean}$ was reduced most in the pharyngeal constrictor muscle (PCM) medius structure, with on average 9.0 ± 4.2 Gy(RBE) (p = 0.0002) compared to the clinical IMPT plans. The average NTCP for grade≥2 and grade≥3 xerostomia at 6 months after treatment significantly decreased with 4.7 ± 1.8% (p = 0.002) and 1.7 ± 0.8% (p = 0.002), respectively, while the average NTCP for grade≥2 and grade≥3 dysphagia decreased with 4.4 ± 2.9% (p = 0.002) and 0.9 ± 0.4% (p = 0.002), respectively, increasing the benefit of protons relative to VMAT. For a "step and shoot" proton arc delivery with auto beam sequencing the estimated delivery time is 11 min, similar to the delivery time of a 6-field IMPT treatment. Gamma analysis between the planned and delivered dose distribution resulted in a 99.99% pass rate using 1mm/1% dose difference/distance to agreement criteria.

CONCLUSIONS

"Step and shoot" proton arc demonstrates potential to further reduce toxicity compared to IMPT and VMAT in OPC treatment. By employing 360 ELs and 30 beams in the proposed ELR method, delivery time can reach clinically acceptable levels without compromising plan toxicity when automatic beam sequencing is available.

An online adaptive plan library approach for intensity modulated proton therapy for head and neck cancer

BACKGROUND AND PURPOSE

In intensity modulated proton therapy (IMPT), the impact of setup errors and anatomical changes is commonly mitigated by robust optimization with population-based setup robustness (SR) settings and offline replanning. In this study we propose and evaluate an alternative approach based on daily plan selection from patient-specific pre-treatment established plan libraries (PLs). Clinical implementation of the PL strategy would be rather straightforward compared to daily online re-planning.

MATERIALS AND METHODS

For 15 head-and-neck cancer patients, the planning CT was used to generate a PL with 5 plans, robustly optimized for increasing SR: 0, 1, 2, 3, 5 mm, and 3% range robustness. Repeat CTs (rCTs) and realistic setup and range uncertainty distributions were used for simulation of treatment courses for the PL approach, treatments with fixed SR (fSR₃) and a trigger-based offline adaptive schedule for 3 mm SR (fSR₃OfA). Daily plan selection in the PL approach was based only on recomputed dose to the CTV on the rCT.

RESULTS

Compared to using fSR₃ and fSR₃OfA, the risk of xerostomia grade ≥ II & III and dysphagia ≥ grade III were significantly reduced with the PL. For 6/15 patients the risk of xerostomia and/or dysphagia ≥ grade II could be reduced by > 2% by using PL. For the other patients, adherence to target coverage constraints was often improved. fSR₃OfA resulted in significantly improved coverage compared to PL for selected patients.

CONCLUSION

The proposed PL approach resulted in overall reduced NTCPs compared to fSR₃ and fSR₃OfA at limited cost in target coverage.

Range uncertainty reductions in proton therapy may lead to the feasibility of novel beam arrangements which improve organ-at-risk sparing

PURPOSE

In proton therapy, dose distributions are currently often conformed to organs at risk (OARs) using the less sharp dose fall-off at the lateral beam edge to reduce the effects of uncertainties in the in vivo proton range. However, range uncertainty reductions may make greater use of the sharper dose fall-off at the distal beam edge feasible, potentially improving OAR sparing. We quantified the benefits of such novel beam arrangements.

METHODS

For each of 10 brain or skull base cases, five treatment plans robust to 2 mm setup and 0%-4% range uncertainty were created for the traditional clinical beam arrangement and a novel beam arrangement making greater use of the distal beam edge to conform the dose distribution to the brainstem. Metrics including the brainstem normal tissue complication probability (NTCP) with the endpoint of necrosis were determined for all plans and all setup and range uncertainty scenarios.

RESULTS

For the traditional beam arrangement, reducing the range uncertainty from the current level of approximately 4% to a potentially achievable level of 1% reduced the brainstem NTCP by up to 0.9 percentage points in the nominal and up to 1.5 percentage points in the worst-case scenario. Switching to the novel beam arrangement at 1% range uncertainty improved these values by a factor of 2, that is, to 1.8 percentage points and 3.2 percentage points, respectively. The novel beam arrangement achieved a lower brainstem NTCP in all cases starting at a range uncertainty of 2%.

CONCLUSION

The benefits of novel beam arrangements may be of the same magnitude or even exceed the direct benefits of range uncertainty reductions. Indirect effects may therefore contribute markedly to the benefits of reducing proton range uncertainties.

Combining inter-observer variability, range and setup uncertainty in a variance-based sensitivity analysis for proton therapy

Margin concepts in proton therapy aim to ensure full dose coverage of the clinical target volume (CTV) in presence of setup and range uncertainty. Due to inter-observer variability (IOV), the CTV itself is uncertain. We present a framework to evaluate the combined impact of IOV, setup and range uncertainty in a variance-based sensitivity analysis (SA). For ten patients with skull base meningioma, the mean calculation time to perform the SA including 1.6 × 10⁴ dose recalculations was 59 min. For two patients in this dataset, IOV had a relevant impact on the estimated CTV D_95% uncertainty.

Technical Note: Does the greater power of pencil beam scanning reduce the need for a proton gantry? A study of head-and-neck and brain tumors

PURPOSE

Proton therapy systems without a gantry can be more compact and less expensive in terms of capital cost, and therefore more available to a larger patient population. Would the advances in pencil beam scanning and robotics make gantry-less treatment possible? In this study, we explore if high-quality treatment plans can be obtained without a gantry.

METHODS AND MATERIALS

We recently showed that proton treatments with the patient in an upright position may be feasible with a new soft robotic immobilization device and imaging which enables multiple possible patient orientations during a treatment. In this study, we evaluate if this new treatment geometry could enable high quality treatment plans without a gantry. We created pencil beam scanning (PBS) treatment plans for seven patients with head-and-neck or brain tumors. Each patient was planned with two scenarios: one with a gantry with the patient in supine position and the other with a gantry-less fixed horizontal beam-line with the patient sitting upright. For the treatment plans, dose-volume-histograms (DVHs), target homogeneity index (HI), mean dose, are reported. A robustness analysis of one plan was performed with 2.5 mm setup errors and 3.5% range uncertainties with nine scenarios.

RESULTS

Most of the PBS-gantry-less plans had similar target HI and OAR mean dose as compared to PBS-gantry plans, and similar robustness with respect to range uncertainties and setup errors.

CONCLUSIONS

Pencil beam scanning provides sufficient power to deliver high quality treatment plans without requiring a gantry for head-and-neck or brain tumors. In combination with the development of the new positioning and immobilization methods required to support this treatment geometry, this work suggests the feasibility of further development of a compact proton therapy system with a fixed horizontal beam-line to treat patients in sitting and reclined positions. This article is protected by copyright. All rights reserved.

Can the mean linear energy transfer of organs be directly related to patient toxicities for current head and neck cancer intensity-modulated proton therapy practice?

BACKGROUND AND PURPOSE

The relative biological effectiveness (RBE) of proton therapy is predicted to vary with the dose-weighted average linear energy transfer (LET_d). However, RBE values may substantially vary for different clinical endpoints. Therefore, the aim of this study was to assess the feasibility of relating mean D⋅LET_d parameters to patient toxicity for HNC patients treated with proton therapy.

MATERIALS AND METHODS

The delivered physical dose (D) and the voxel-wise product of D and LET_d (D⋅LET_d) distributions were calculated for 100 head and neck cancer (HNC) proton therapy patients using our TPS (Raystation v6R). The means and covariance matrix of the accumulated D and D⋅LET_d of all relevant organs-at-risk (OARs) were used to simulate 2.500 data sets of different sizes. For each dataset, an attempt was made to add mean D⋅LET_d parameters to a multivariable NTCP model based on mean D parameters of the same OAR for xerostomia, tube feeding and dysphagia. The likelihood of creating an NTCP model with statistically significant parameters (i.e. power) was calculated as a function of the simulated sample size for various RBE models.

RESULTS

The sample size required to have a power of at least 80% to show an independent effect of mean D⋅LET_d parameters on toxicity is over 15,000 patients for all toxicities.

CONCLUSION

For current clinical practice, it is not feasible to directly model NTCP with both mean D and mean D⋅LET_d of OARs. These findings should not be interpreted as a contradiction of previous evidence for the relationship between RBE and LET_d.

Proton range uncertainty reduction benefits for skull base tumors in terms of normal tissue complication probability (NTCP) and healthy tissue doses

PURPOSE

Proton therapy allows for more conformal dose distributions and lower organ at risk and healthy tissue doses than conventional photon-based radiotherapy, but uncertainties in the proton range currently prevent proton therapy from making full use of these advantages. Numerous developments therefore aim to reduce such range uncertainties. In this work, we quantify the benefits of reductions in range uncertainty for treatments of skull base tumors.

METHODS

The study encompassed 10 skull base patients with clival tumors. For every patient, six treatment plans robust to setup errors of 2 mm and range errors from 0% to 5% were created. The determined metrics included the brainstem and optic chiasm normal tissue complication probability (NTCP) with the endpoints of necrosis and blindness, respectively, as well as the healthy tissue volume receiving at least 70% of the prescription dose.

RESULTS

A range uncertainty reduction from the current level of 4% to a potentially achievable level of 1% reduced the probability of brainstem necrosis by up to 1.3 percentage points in the nominal scenario in which neither setup nor range errors occur and by up to 2.9 percentage points in the worst-case scenario. Such a range uncertainty reduction also reduced the optic chiasm NTCP with the endpoint of blindness by up to 0.9 percentage points in the nominal scenario and by up to 2.2 percentage points in the worst-case scenario. The decrease in the healthy tissue volume receiving at least 70% of the prescription dose ranged from -7.8 to 24.1 cc in the nominal scenario and from -3.4 to 38.4 cc in the worst-case scenario.

CONCLUSION

The benefits quantified as part of this study serve as a guideline of the OAR and healthy tissue dose benefits that range monitoring techniques may be able to achieve. Benefits were observed between all levels of range uncertainty. Even smaller range uncertainty reductions may therefore be beneficial.