Dose averaged linear energy transfer optimization for large sacral chordomas in carbon ion therapy

BACKGROUND

Carbon ion beams are well accepted as densely ionizing radiation with a high linear energy transfer (LET). However, the current clinical practice does not fully exploit the highest possible dose-averaged LET (LETd) and, consequently, the biological potential in the target. This aspect becomes worse in larger tumors for which inferior clinical outcomes and corresponding lower LETd was reported.

PURPOSE

The vicinity to critical organs in general and the inferior overall survival reported for larger sacral chordomas treated with carbon ion radiotherapy (CIRT), makes the treatment of such tumors challenging. In this work it was aimed to increase the LETd in large volume tumors while maintaining the relative biological effectiveness (RBE)-weighted dose, utilizing the LETd optimization functions of a commercial treatment planning system (TPS).

METHODS

Ten reference sequential boost carbon ion treatment plans, designed to mimic clinical plans for large sacral chordoma tumors, were generated. High dose clinical target volumes (CTV-HD) larger than250 cm 3 $250 \,{\rm cm}^{3}$ were considered as large targets. The total RBE-weighted median dose prescription with the local effect model (LEM) wasD RBE , 50 % = 73.6 Gy $\textrm {D}_{\rm RBE, 50\%}=73.6 \,{\rm Gy}$ in 16 fractions (nine to low dose and seven to high dose planning target volume). No LETd optimization was performed in the reference plans, while LETd optimized plans used the minimum LETd (Lmin) optimization function in RayStation 2023B. Three different Lmin values were investigated and specified for the seven boost fractions:L min = 60 keV / μ m $\textrm {L}_{\rm min}=60 \,{\rm keV}/{\umu }{\rm m}$ ,L min = 80 keV / μ m $\textrm {L}_{\rm min}=80 \,{\rm keV}/{\umu }{\rm m}$ andL min = 100 keV / μ m $\textrm {L}_{\rm min}=100 \,{\rm keV}/{\umu }{\rm m}$ . To compare the LETd optimized against reference plans, LETd and RBE-weighted dose based goals similar to and less strict than clinical ones were specified for the target. The goals for the organs at risk (OAR) remained unchanged. Robustness evaluation was studied for eight scenarios (± 3.5 % $\pm 3.5\%$ range uncertainty and± 3 mm $\pm 3 \,{\rm mm}$ setup uncertainty along the main three axes).

RESULTS

The optimization method withL min = 60 keV / μ m $\textrm {L}_{\rm min}=60 \,{\rm keV}/{\umu }{\rm m}$ resulted in an optimal LETd distribution with an average increase ofLET d , 98 % ${\rm {LET}}_{{\rm {d,}}98\%}$ (andLET d , 50 % ${\rm {LET}}_{{\rm {d,}}50\%}$ ) in the CTV-HD by8.9 ± 1.5 keV / μ m $8.9\pm 1.5 \,{\rm keV}/{\umu }{\rm m}$ (27 % $27\%$ ) (and6.9 ± 1.3 keV / μ m $6.9\pm 1.3 \,{\rm keV}/{\umu }{\rm m}$ (17 % $17\%$ )), without significant difference in the RBE-weighted dose. By allowing± 5 % $\pm 5\%$ over- and under-dosage in the target, theLET d , 98 % ${\rm {LET}}_{{\rm {d,}}98\%}$ (andLET d , 50 % ${\rm {LET}}_{{\rm {d,}}50\%}$ ) can be increased by11.3 ± 1.2 keV / μ m $11.3\pm 1.2 \,{\rm keV}/{\umu }{\rm m}$ (34 % $34\%$ ) (and11.7 ± 3.4 keV / μ m $11.7\pm 3.4 \,{\rm keV}/{\umu }{\rm m}$ (29 % $29\%$ )), using the optimization parametersL min = 80 keV / μ m $\textrm {L}_{\rm min}=80 \,{\rm keV}/{\umu }{\rm m}$ . The pass rate for the OAR goals in the LETd optimized plans was in the same level as the reference plans. LETd optimization lead to less robust plans compared to reference plans.

CONCLUSIONS

Compared to conventionally optimized treatment plans, the LETd in the target was increased while maintaining the RBE-weighted dose using TPS LETd optimization functionalities. Regularly assessing RBE-weighted dose robustness and acquiring more in-room images remain crucial and inevitable aspects during treatment.

Robustness of carbon-ion radiotherapy against DNA damage repair associated radiosensitivity variation based on a biophysical model

BACKGROUND

Interpatient variation of tumor radiosensitivity is rarely considered during the treatment planning process despite its known significance for the therapeutic outcome.

PURPOSE

To apply our mechanistic biophysical model to investigate the biological robustness of carbon ion radiotherapy (CIRT) against DNA damage repair interference (DDRi) associated patient-to-patient variability in radiosensitivity and its potential clinical advantages against conventional radiotherapy approaches.

METHODS AND MATERIALS

The "UNIfied and VERSatile bio response Engine" (UNIVERSE) was extended by carbon ions and its predictions were compared to a panel of in vitro and in vivo data including various endpoints and DDRi settings within clinically relevant dose and linear energy transfer (LET) ranges. The implications of UNIVERSE predictions were then assessed in a clinical patient scenario considering DDRi variance.

RESULTS

UNIVERSE tests well against the applied benchmarks. While in vitro survival curves were predicted with an R2 > 0.92, deviations from in vivo RBE data were less than 5.6% The conducted paradigmatic patient plan study implies a markedly reduced significance of DDRi based radiosensitivity variability in CIRT (13% change ofD 50 ${{D}_{50}}$ in target) compared to conventional radiotherapy (62%) and that boosting the LET within the target further amplifies this robustness of CIRT (8%). In the case of heightened tumor radiosensitivity, a dose de-escalation strategy for photons allows a reduction of the maximum effective dose within the normal tissue (NT) from aD 2 ${{D}_2}$ of 2.65 to 1.64 Gy, which lies below the level found for CIRT (D 2 ${{D}_2}$  = 2.41 Gy) for the analyzed plan and parameters. However, even after de-escalation, the integral effective dose in the NT is found to be substantially higher for conventional radiotherapy in comparison to CIRT (D m e a n ${{D}_{mean}}$ of 0.75, 0.46, and 0.24 Gy for the conventional plan, its de-escalation and CIRT, respectively).

CONCLUSIONS

The framework offers adequate predictions of in vitro and in vivo radiation effects of CIRT while allowing the consideration of DRRi based solely on parameters derived from photon data. The results of the patient planning study underline the potential of CIRT to minimize important sources of interpatient divergence in therapy outcome, especially when combined with techniques that allow to maximize the LET within the tumor. Despite the potential of de-escalation strategies for conventional radiotherapy to reduce the maximum effective dose in the NT, CIRT appears to remain a more favorable option due to its ability to reduce the integral effective dose within the NT.

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.

Investigation on the physical dose filtered by linear energy transfer for treatment plan evaluation in carbon ion therapy

BACKGROUND

Large tumor size has been reported as a predicting factor for inferior clinical outcome in carbon ion radiotherapy (CIRT). Besides the clinical factors accompanied with such tumors, larger tumors receive typically more low linear energy transfer (LET) contributions than small ones which may be the underlying physical cause. Although dose averaged LET is often used as a single parameter descriptor to quantify the beam quality, there is no evidence that this parameter is the optimal clinical predictor for the complex mixed radiation fields in CIRT.

PURPOSE

Purpose of this study was to investigate on a novel dosimetric quantity, namely high-LET-dose (D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ , the physical dose filtered based on an LET threshold) as a single parameter estimator to differentiate between carbon ion treatment plans (cTP) with a small and large tumor volume.

METHODS

Ten cTPs with a planning target volume,PTV ≥ 500 cm 3 $\mathrm{PTV}\ge {500}\,{{\rm cm}^{3}}$ (large) and nine with aPTV < 500 cm 3 $\mathrm{PTV}<{500}\,{{\rm cm}^{3}}$ (small) were selected for this study. To find a reasonable LET threshold (L thr $\textrm {L}_{\textrm {thr}}$ ) that results in a significant difference in terms ofD > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ , the voxel based normalized high-LET-dose (D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ ) distribution in the clinical target volume (CTV) was studied on a subset (12 out of 19 cTPs) for 18 LET thresholds, using standard distribution descriptors (mean, variance and skewness). The classical dose volume histogram concept was used to evaluate theD > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ andD ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ distributions within the target of all 19 cTPs at the before determinedL thr $\textrm {L}_{\textrm {thr}}$ . Statistical significance of the difference between the two groups in terms of meanD > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ andD ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ volume histogram parameters was evaluated by means of (two-sided) t-test or Mann-Whitney-U-test. In addition, the minimum target coverage at the above determinedL thr $\textrm {L}_{\textrm {thr}}$ was compared and validated against three other thresholds to verify its potential in differentiation between small and large volume tumors.

RESULTS

AnL thr $\textrm {L}_{\textrm {thr}}$ of approximately30 keV / μ m ${30}\,{\rm keV/}\umu {\rm m}$ was found to be a reasonable threshold to classify the two groups. At this threshold, theD > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ andD ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ were significantly larger (p < 0.05 $p<0.05$ ) in small CTVs. For the small tumor group, the near-minimum and medianD > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ (andD ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ ) in the CTV were in average9.3 ± 1.5 Gy $9.3\pm {1.5}\,{\rm Gy}$ (0.31 ± 0.08) and13.6 ± 1.6 Gy $13.6\pm {1.6}\,{\rm Gy}$ (0.46 ± 0.06), respectively. For the large tumors, these parameters were6.6 ± 0.2 Gy $6.6\pm {0.2}\,{\rm Gy}$ (0.20 ± 0.01) and8.6 ± 0.4 Gy $8.6\pm {0.4}\,{\rm Gy}$ (0.28 ± 0.02). The difference between the two groups in terms of mean near-minimum and medianD > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ (D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ ) was 2.7 Gy (11%) and 5.0 Gy (18%), respectively.

CONCLUSIONS

The feasibility of high-LET-dose based evaluation was shown in this study where a lowerD > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ was found in cTPs with a large tumor size. Further investigation is needed to draw clinical conclusions. The proposed methodology in this work can be utilized for future high-LET-dose based studies.

The LET trilemma: Conflicts between robust target coverage, uniform dose, and dose-averaged LET in carbon therapy

BACKGROUND

Linear energy transfer (LET) is closely related to the biological effect of ionizing radiation. Increasing the dose-averaged LET (LETd ) within the target volume has been proposed as a means to improve clinical outcome for hypoxic tumors. However, doing so can lead to reduced robustness to range uncertainty.

PURPOSE

To quantify the relationship between robust target coverage, target dose uniformity, and LETd , we employ robust optimization using dose-based and LETd -based functions and allow varying amounts of target non-uniformity.

METHODS AND MATERIALS

Robust carbon therapy optimization is used to create plans for phantom cases with increasing target sizes (radii 1, 3, and 5 cm). First, the influence of respectively range and setup uncertainty on the LETd in the target is studied. Second, we employ strategies allowing overdosage in the clinical target volume (CTV) or gross tumor volume (GTV), which enable increased LETd in the target. The relationship between robust target coverage and LETd in the target is illustrated by tradeoff curves generated by optimization using varying weights for the LETd -based functions.

RESULTS

As the range uncertainty used in the robust optimization increased from 0% to 5%, the near-minimum nominal LETd decreased by 17%-29% (9-21 keV/µm) for the different target sizes. The effect of increasing setup uncertainty was marginal. Allowing 10% overdosage in the CTV enabled 9%-29% (6-12 keV/µm) increased near-minimum worst case LETd for the different target sizes, compared to uniform dose plans. When 10% overdosage was allowed in the GTV only, the increase was 1%-20% (1-8 keV/µm).

CONCLUSIONS

There is an inherent conflict between range uncertainty robustness and high LETd in the target, which is aggravated with increasing target size. For large tumors, it is possible to simultaneously achieve two of the three qualities range robustness, uniform dose, and high LETd in the target.

Treatment planning of scanned proton beams in RayStation

The use of scanned proton beams in external beam radiation therapy has seen a rapid development over the past decade. This technique places new demands on treatment planning, as compared to conventional photon-based radiation therapy. In this article, several proton specific functions as implemented in the treatment planning system RayStation are presented. We will cover algorithms for energy layer and spot selection, basic optimization including the handling of spot weight limits, optimization of the linear energy transfer (LET) distribution, robust optimization including the special case of 4D optimization, proton arc planning, and automatic planning using deep learning. We will further present the Monte Carlo (MC) proton dose engine in RayStation to some detail, from the material interpretation of the CT data, through the beam model parameterization, to the actual MC transport mechanism. Useful tools for plan evaluation, including robustness evaluation, and the versatile scripting interface are also described. The overall aim of the paper is to give an overview of some of the key proton planning functions in RayStation, with example usages, and at the same time provide the details about the underlying algorithms that previously have not been fully publicly available.

Partitioning of discrete proton arcs into interlaced subplans can bring proton arc advances to existing proton facilities

BACKGROUND

Proton arcs have shown potential to reduce the dose to organs at risks (OARs) by delivering the protons from many different directions. While most previous studies have been focused on dynamic arcs (delivery during rotation), an alternative approach is discrete arcs, where step-and-shoot delivery is used over a large number of beam directions. The major advantage of discrete arcs is that they can be delivered at existing proton facilities. However, this advantage comes at the expense of longer treatment times.

PURPOSE

To exploit the dosimetric advantages of proton arcs, while achieving reasonable delivery times, we propose a partitioning approach where discrete arc plans are split into subplans to be delivered over different fractions in the treatment course.

METHODS

For three oropharyngeal cancer patients, four different arc plans have been created and compared to the corresponding clinical IMPT plan. The treatment plans are all planned to be delivered in 35 fractions, but with different delivery approaches over the fractions. The first arc plan (1×30) has 30 directions to be delivered every fraction, while the others are partitioned into subplans with 10 and 6 beam directions, each to be delivered every third (3×10), fifth fraction (5×6), or seventh fraction (7×10). All plans are assessed with respect to delivery time, target robustness over the treatment course, doses to OARs and NTCP for dysphagia and xerostomia.

RESULTS

The delivery time (including an additional delay of 30 s between the discrete directions to simulate manual interaction with the treatment control system) is reduced from on average 25.2 min for the 1×30 plan to 9.2 min for the 3×10 and 7×10 plans and 5.7 min for the 5×6 plans. The delivery time for the IMPT plan is 7.9 min. When accounting for the combination of delivery time, target robustness, OAR sparing, and NTCP reduction, the plans with 10 directions in each fraction are the preferred choice. Both the 3×10 and 7×10 plans show improved target robustness compared to the 1×30 plans, while keeping OAR doses and NTCP values at almost as low levels as for the 1×30 plans. For all patients the NTCP values for dysphagia are lower for the partitioned plans with 10 directions compared to the IMPT plans. NTCP reduction for xerostomia compared to IMPT is seen in two of the three patients. The best results are seen for the first patient, where the NTCP reductions for the 7×10 plan are 1.6 p.p. (grade 2 xerostomia) and 1.5 p.p. (grade 2 dysphagia). The corresponding NTCP reductions for the 1×30 plan are 2.7 p.p. (xerostomia, grade 2) and 2.0 p.p. (dysphagia, grade 2).

CONCLUSIONS

Discrete proton arcs can be implemented at any proton facility with reasonable treatment times using a partitioning approach. The technique also makes the proton arc treatments more robust to changes in the patient anatomy.

Technical note: In silico benchmarking of the linear energy transfer-based functionalities for carbon ion beams in a commercial treatment planning system

BACKGROUND

The increasing number of studies dealing with linear energy transfer (LET)-based evaluation and optimization in the field of carbon ion radiotherapy (CIRT) indicates the rising demand for LET implementation in commercial treatment planning systems (TPS). Benchmarking studies could play a key role in detecting (and thus preventing) computation errors prior implementing such functionalities in a TPS.

PURPOSE

This in silico study was conducted to benchmark the following two LET-related functionalities in a commercial TPS against Monte Carlo simulations: (1) dose averaged LET (LET_d ) scoring and (2) physical dose filtration based on LET for future LET-based treatment plan evaluation and optimization studies.

METHODS

The LET_d scoring and LET-based dose filtering (in which the deposited dose can be separated into the dose below and above the user specified LET threshold) functionalities for carbon ions in the research version RayStation (RS) 9A-IonPG TPS (RaySearch Laboratories, Sweden) were benchmarked against GATE/Geant4 simulations. Pristine Bragg peaks (BPs) and cuboid targets, positioned at different depths in a homogeneous water phantom and a setup with heterogeneity were used for this study.

RESULTS

For all setups (homogeneous and heterogeneous), the mean absolute (and relative) LET_d difference was less than 1 keV/ μ $\umu$ m (3.5%) in the plateau and target and less than 2 keV/ μ $\umu$ m (8.3%) in the fragmentation tail. The maximum local differences were 4 and 6 keV/ μ $\umu$ m, respectively. The mean absolute (and relative) physical dose differences for both low- and high-LET doses were less than 1 cGy (1.5%) in the plateau, target and tail with a maximum absolute difference of 2 cGy.

CONCLUSIONS

No computation error was found in the tested functionalities except for LET_d in lateral direction outside the target, showing the limitation of the implemented monochrome model in the tested TPS version.

Technical note: Impact of beamline-specific particle energy spectra on clinical plans in carbon ion beam therapy

PURPOSE

The Local Effect Model version one (LEM I) is applied clinically across Europe to quantify the relative biological effectiveness (RBE) of carbon ion beams. It requires the full particle fluence spectrum differential in energy in each voxel as input parameter. Treatment planning systems (TPSs) use beamline-specific look-up tables generated with Monte Carlo (MC) codes. In this study, the changes in RBE weighted dose were quantified using different levels of details in the simulation or different MC codes.

METHODS

The particle fluence differential in energy was simulated with FLUKA and Geant4 at 500 depths in water in 1-mm steps for 58 initial carbon ion energies (between 120.0 and 402.8 MeV/u). A dedicated beam model was applied, including the full description of the Nozzle using GATE-RTionV1.0 (Geant4.10.03p03). In addition, two tables generated with FLUKA were compared. The starting points of the FLUKA simulations were phase space (PhS) files from, firstly, the Geant4 nozzle simulations, and secondly, a clinical beam model where an analytic approach was used to mimic the beamline. Treatment plans (TPs) were generated with RayStation 8B (RaySearch Laboratories AB, Sweden) for cubic targets in water and 10 clinical patient cases using the clinical beam model. Subsequently, the RBE weighted dose was re-computed using the two other fluence tables (FLUKA PhS or Geant4).

RESULTS

The fluence spectra of the primary and secondary particles simulated with Geant4 and FLUKA generally agreed well for the primary particles. Differences were mainly observed for the secondary particles. Interchanging the two energy spectra (FLUKA vs. GEANT4) to calculate the RBE weighted dose distributions resulted in average deviations of less than 1% in the entrance up to the end of the target region, with a maximum local deviation at the distal edge of the target. In the fragment tail, larger discrepancies of up to 5% on average were found for deep-seated targets. The patient and water phantom cases demonstrated similar results.

CONCLUSION

RBE weighted doses agreed well within all tested setups, confirming the clinical beam model provided by the TPS vendor. Furthermore, the results showed that the open source and generally available MC code Geant4 (in particular using GATE or GATE-RTion) can also be used to generate basic beam data required for RBE calculation in carbon ion therapy.

Fast robust optimization of proton PBS arc therapy plans using early energy layer selection and spot assignment

Objective. Proton pencil-beam scanning arcs (PBS arcs) have gained much attention during the past years, due to its potential for increased clinical benefit compared to conventional proton therapy. Previous studies on PBS arcs have primarily been focused on plan quality, and lately efforts have been made to reduce the delivery time. However, the methods presented so far suffer from slow optimization processes. Approach. We present a new method for fast robust optimization of PBS arc plans. The new method assigns a single energy layer per discretized direction prior to spot weight optimization and reduces the number of initial spots considerably compared to conventional methods. We used the new method for three prostate cancer patients with a prescribed dose to the CTV of 77 GyRBE in 35 fractions. For each of the patients, four plans were created: 2-beam IMPT (2IMPT), 1-beam PBS arc (1Arc), 1-beam PBS arc without focus on reducing upward energy jumps (1Arc_unseq) and two-beam PBS arc (2Arc). Main results. All PBS arc plans show a reduced integral dose compared to their respective 2IMPT plans. In the nominal case, the average CTV D98 and D2 metrics over the three patients were best for the 2Arc, followed by 2IMPT ( D 98 ¯ / D 2 ¯ : 7523/7986 cGyRBE (2IMPT), 7478/7984 cGy (1Arc), 7486/7951 cGy (1Arc_unseq), 7531/7951 cGyRBE (2Arc)). The average robust target coverage in terms of V95 of the voxelwise minimum dose distribution (evaluated over 42 scenarios) was: 98.0% (2IMPT), 88.6% (1Arc), 92.5% (1Arc_unseq), 97.3% (2Arc). The optimization time, including spot selection and spot dose computation, is longest for the 2Arc plan, but is below 6 min for all patients. The maximum estimated delivery time for all types of arc plans is just above 5 min Significance. The ability for efficient treatment planning constitutes an important step towards clinical introduction of proton PBS arcs.

4D robust optimization including uncertainties in time structures can reduce the interplay effect in proton pencil beam scanning radiation therapy

PURPOSE

Interplay effects in proton radiotherapy can create large distortions in the dose distribution and severely degrade the plan quality. Standard methods to mitigate these effects include abdominal compression, gating, and rescanning. We propose a new method to include the time structures of the delivery and organ motion in the framework of four-dimensional (4D) robust optimization to generate plans that are robust against interplay effects.

METHODS

The method considers multiple scenarios reflecting the uncertainties in the delivery and in the organ motion. In each scenario, the pencil beam scanning spots are distributed to different phases of the breathing cycle according to each individual spot time stamp, and a partial beam dose is calculated for each phase. The partial beam doses are accumulated on a reference phase through deformable image registrations. Minimax optimization is performed to take all scenarios into account simultaneously. For simplicity, the uncertainties in this proof of concept study are limited to variations in the breathing pattern. The method is evaluated for three different nonsmall cell lung cancer patients and compared to plans using conventional 4D robust optimization both with and without rescanning. We assess the ability of the method to mitigate distortions from the interplay effect over multiple evaluation scenarios using 4D dose calculations. This interplay evaluation is performed in an experimentally validated framework, which is independent of the optimization in the plan generation step.

RESULTS

For the three studied patients, 4D optimization including time structures is efficient, especially for large tumor motions, where rescanning of conventional 4D robustly optimized plans is not sufficient to mitigate the interplay effect. The most efficient approach of the new method is achieved when it is combined with rescanning. For the patient with the largest motion, the mean V95% is 99.2% and mean V107% is 3.65% for the best rescanned 4D plan optimized with time structure. This can be compared to conventional 4D optimized plans with mean V95% of 92.7% and mean V107% of 13.1%.

CONCLUSIONS

The current study shows the potential of reducing interplay effects in proton pencil beam scanning radiotherapy by incorporating organ motion and delivery characteristics in a 4D robust optimization.

Effectiveness of different rescanning techniques for scanned proton radiotherapy in lung cancer patients

Non-small cell lung cancer (NSCLC) is a tumour type thought to be well-suited for proton radiotherapy. However, the lung region poses many problems related to organ motion and can for actively scanned beams induce severe interplay effects. In this study we investigate four mitigating rescanning techniques: (1) volumetric rescanning, (2) layered rescanning, (3) breath-sampled (BS) layered rescanning, and (4) continuous breath-sampled (CBS) layered rescanning. The breath-sampled methods will spread the layer rescans over a full breathing cycle, resulting in an improved averaging effect at the expense of longer treatment times. In CBS, we aim at further improving the averaging by delivering as many rescans as possible within one breathing cycle. The interplay effect was evaluated for 4D robustly optimized treatment plans (with and without rescanning) for seven NSCLC patients in the treatment planning system RayStation. The optimization and final dose calculation used a Monte Carlo dose engine to account for the density heterogeneities in the lung region. A realistic treatment delivery time structure given from the IBA ScanAlgo simulation tool served as basis for the interplay evaluation. Both slow (2.0 s) and fast (0.1 s) energy switching times were simulated. For all seven studied patients, rescanning improves the dose conformity to the target. The general trend is that the breath-sampled techniques are superior to layered and volumetric rescanning with respect to both target coverage and variability in dose to OARs. The spacing between rescans in our breath-sampled techniques is set at planning, based on the average breathing cycle length obtained in conjunction with CT acquisition. For moderately varied breathing cycle lengths between planning and delivery (up to 15%), the breath-sampled techniques still mitigate the interplay effect well. This shows the potential for smooth implementation at the clinic without additional motion monitoring equipment.