Beam configuration selection for robust intensity-modulated proton therapy in cervical cancer using Pareto front comparison

Dosimetric impact of bone marrow sparing for robustly optimized IMPT for locally advanced cervical cancer

BACKGROUND AND PURPOSE

To investigate the trade-off between bone marrow sparing (BMS) and dose to organs at risk (OARs) for intensity modulated proton therapy (IMPT) for women with locally advanced cervical cancer (LACC).

MATERIALS AND METHODS

Twenty LACC patients were retrospectively included. IMPT plans were created for each patient using automated treatment planning. These plans progressively reduced bone marrow mean doses by steps of 1 GyRBE, while constraining target coverage and conformality. The relation between bone marrow dose and bladder, small bowel, rectum, and sigmoid doses was evaluated.

RESULTS

A total of 140 IMPT plans were created. Plans without BMS had an average [range] bone marrow mean dose of 17.3 [14.7-21.6] GyRBE , which reduced to 12.0 [10.0-14.0] GyRBE with maximum BMS. The mean OAR dose [range] increased modestly for 1 GyRBE BMS: 0.2 [0.0 - 0.6] GyRBE for bladder, 0.3 [-0.2 - 0.7] GyRBE for rectum, 0.4 [0.1 - 0.8] GyRBE for small bowel, and 0.2 [-0.2 - 0.4] GyRBE for sigmoid. Moreover, for maximum BMS, mean OAR doses [range] escalated by 3.3 [0.1 - 6.7] GyRBE for bladder, 5.8 [1.8 - 12.4] GyRBE for rectum, 3.9 [1.6 - 5.9] GyRBE for small bowel, and 2.7 [0.6 - 5.9] GyRBE for sigmoid.

CONCLUSION

Achieving 1 GyRBE BMS for IMPT is feasible for LACC patients with limited dosimetric impact on other OARs. While further bone marrow dose reduction is possible for some patients, it may increase OAR doses substantially for others. Hence, we recommend a personalized approach when introducing BMS into clinical IMPT treatment planning to carefully assess individual patient benefits and risks.

Reflections on beam configuration optimization for intensity-modulated proton therapy

Presumably, intensity-modulated proton radiotherapy (IMPT) is the most powerful form of proton radiotherapy. In the current state of the art, IMPT beam configurations (i.e. the number of beams and their directions) are, in general, chosen subjectively based on prior experience and practicality. Beam configuration optimization (BCO) for IMPT could, in theory, significantly enhance IMPT’s therapeutic potential. However, BCO is complex and highly computer resource-intensive. Some algorithms for BCO have been developed for intensity-modulated photon therapy (IMRT). They are rarely used clinically mainly because the large number of beams typically employed in IMRT renders BCO essentially unnecessary. Moreover, in the newer form of IMRT, volumetric modulated arc therapy, there are no individual static beams. BCO is of greater importance for IMPT because it typically employs a very small number of beams (2-4) and, when the number of beams is small, BCO is critical for improving plan quality. However, the unique properties and requirements of protons, particularly in IMPT, make BCO challenging. Protons are more sensitive than photons to anatomic changes, exhibit variable relative biological effectiveness along their paths, and, as recently discovered, may spare the immune system. Such factors must be considered in IMPT BCO, though doing so would make BCO more resource intensive and make it more challenging to extend BCO algorithms developed for IMRT to IMPT. A limited amount of research in IMPT BCO has been conducted; however, considerable additional work is needed for its further development to make it truly effective and computationally practical. This article aims to provide a review of existing BCO algorithms, most of which were developed for IMRT, and addresses important requirements specific to BCO for IMPT optimization that necessitate the modification of existing approaches or the development of new effective and efficient ones.

Dosimetric Impact of Intrafraction Motion in Online-Adaptive Intensity Modulated Proton Therapy for Cervical Cancer

PURPOSE

A method was recently developed for online-adaptive intensity modulated proton therapy (IMPT) in patients with cervical cancer. The advantage of this approach, relying on the use of tight margins, is challenged by the intrafraction target motion. The purpose of this study was to evaluate the dosimetric effect of intrafraction motion on the target owing to changes in bladder filling in patients with cervical cancer treated with online-adaptive IMPT.

METHODS AND MATERIALS

In 10 patients selected to have large uterus motion induced by bladder filling, the intrafraction anatomic changes were simulated for several prefraction durations for online (automated) contouring and planning. For each scenario, the coverage of the primary target was evaluated with margins of 2.5 and 5 mm.

RESULTS

Using a 5- mm planning target volume margin, median accumulated D98% was greater than 42.75 Gy_RBE1.1 (95% of the prescribed dose) in the case of a prefraction duration of 5 and 10 minutes. For a prefraction duration of 15 minutes, this parameter deteriorated to 42.6 Gy_RBE1.1. When margins were reduced to 2.5 mm, only a 5-minute duration resulted in median target D98% above 42.75 Gy_RBE1.1. In addition, smaller bladders were found to be associated with larger dose degradations compared with larger bladders.

CONCLUSIONS

This study indicates that intrafraction anatomic changes can have a substantial dosimetric effect on target coverage in an online-adaptive IMPT scenario for patients subject to large uterus motion. A margin of 5 mm was sufficient to compensate for the intrafraction motion due to bladder filling for up to 10 minutes of prefraction time. However, compensation for the uncertainties that were disregarded in this study, by using margins or robust optimization, is also required. Furthermore, a large bladder volume restrains intrafraction target motion and is recommended for treating patients in this scenario. Assuming that online-adaptive IMPT remains beneficial as long as narrow margins are used (5 mm or below), this study demonstrates its feasibility with regard to intrafraction motion.

Inter-fraction motion robustness and organ sparing potential of proton therapy for cervical cancer

PURPOSE

Large-field photon radiotherapy is current standard in the treatment of cervical cancer patients. However, with the increasing availability of Pencil Beam Scanning Proton Therapy (PBS-PT) and robust treatment planning techniques, protons may have significant advantages for cervical cancer patients in the reduction of toxicity. In this study, PBS-PT and photon Volumetric Modulated Arc Therapy (VMAT) were compared, examining target coverage and organ at risk (OAR) dose, taking inter- and intra-fraction motion into account.

MATERIALS AND METHODS

Twelve cervical cancer patients were included in this in-silico planning study. In all cases, a planning CT scan, five weekly repeat CT scans (reCTs) and an additional reCT 10 min after the first reCT were available. Two-arc VMAT and robustly optimised two- and four-field (2F and 4F) PBS-PT plans were robustly evaluated on planCTs and reCTs using set-up and range uncertainty. Nominal OAR doses and voxel-wise minimum target coverage robustness were compared.

RESULTS

Average voxel-wise minimum accumulated doses for pelvic target structures over all patients were adequate for both photon and proton treatment techniques (D98 > 95%, [91.7-99.3%]). Average accumulated dose of the para-aortic region was lower than the required 95%, D98 > 94.4% [91.1-98.2%]. With PBS-PT 4F, dose to all OARs was significantly lower than with VMAT. Major differences were observed for mean bowel bag V_15Gy: 60% [39-70%] for VMAT vs 30% [10-52%] and 32% [9-54%] for PBS-PT 2F and 4F and for mean bone marrow V_10Gy: 88% [82-97%] for VMAT vs 66% [60-73%] and 67% [60-75%] for PBS-PT 2F and 4F.

CONCLUSION

Robustly optimised PBS-PT for cervical cancer patients shows equivalent target robustness against inter- and intra-fraction variability compared to VMAT, and offers significantly better OAR sparing.

A Pareto-based beam orientation optimization method for spot scanning intensity-modulated proton therapy

PURPOSE

To provide a proof of principle of a Pareto-based method to automatically generate optimal intensity-modulated proton therapy (IMPT) plans for various noncoplanar beam orientations.

METHODS

A novel multicriteria beam orientation optimization (MCBOO) method was developed to generate Pareto database of optimal plans. The MCBOO method automatically explores the beam orientations and the scalarization parameters of the IMPT plans simultaneously. The MCBOO method is based on multicriteria bilevel optimization (i.e., hierarchical optimization with two nested levels, named the upper and lower level optimization). In MCBOO, the upper level optimization explores the noncoplanar beam orientation space, while the lower level explores the scalarization parameters for a given beam orientation. Differential evolution method was used in both levels, and the Pareto optimal plans were aggregated from the bilevel optimizations to construct the Pareto database. The MCBOO method was implemented on a multinode multi-GPU cluster, and it was tested on three brain tumor patient cases. The Pareto database of the three patients was generated for a set of DVH-based objectives. A statistical analysis was performed between a selected set of MCBOO plans and the manual plan (plan with manually selected beam orientation based on the clinical experience and optimized with the same single plan iterative optimizer used in the MCBOO). The selected set of MCBOO plans consisted of plans that matched the performance of the manual plan [i.e., MCBOO plans that have the same target coverage (within 2%) as the manual plan or better and achieved the same dose (within 2%) or lower to all of the organs at risks (OARs) but one OAR]. Additionally, a dosimetric comparison between of one of the selected MCBOO plans vs the manual plan was conducted.

RESULTS

The multicriteria beam orientation optimization algorithm automatically generated Pareto plans for the three noncoplanar brain tumor cases. The MCBOO plans provided an alternative objective trade-offs to the manual plan. The selected MCBOO plans showed a reduction in dose to multiple organs at risk vs the manual plan with a maximum value which ranged between 10.8 and 12.9 Gy for the three patients. The trade-off of the OAR dose reduction resulted in higher dose to no more than one OAR for each of the selected MCBOO plans vs the manual plan. The maximum dose increase in the MCBOO plans over the manual plan ranged from 7.8 to 11.8 Gy.

CONCLUSIONS

A novel multicriteria beam orientation optimization method was developed and tested on three IMPT patient cases. The method automatically generates Pareto plans database by exploring the noncoplanar beam orientations. The method was able to identify beam orientations with Pareto optimal plans that are comparable to the manually created plans with varying objective trade-offs.

Adaptive method for multicriteria optimization of intensity-modulated proton therapy

PURPOSE

Provide an adaptive multicriteria optimization (MCO) method for intensity-modulated proton therapy (IMPT) utilizing GPU technology. Previously described limitations of MCO such as Pareto approximation and limitation on the number of objectives were addressed.

METHODS

The treatment planning process for IMPT must account for multiple objectives, which requires extensive treatment planning resources. Often a large number of objectives (>10) are required. Hence the need for an MCO algorithm that can handle large number of objectives. The novelty of the MCO method presented here lies on the introduction of the adaptive weighting scheme that can generate a well-distributed and dense representation of the Pareto surface for a large number of objectives in an efficient manner. In our approach the generated Pareto surface is constructed for a set of DVH objectives. The MCO algorithm is based on the augmented weighted Chebychev metric (AWCM) method with an adaptive weighting scheme. This scheme uses the differential evolution (DE) method to generate a set of well-distributed Pareto points. The quality of the Pareto points' distribution in the objective space was assessed quantitatively using the Pareto sampling metric. The MCO algorithm was developed to perform multiple parallel searches to achieve a rapid mapping of the Pareto surface, produce clinically deliverable plans, and was implemented on a GPU cluster. The MCO algorithm was tested on two clinical cases with 10 and 18 objectives. For each case one of the MCO-generated plans was selected for comparison with the clinically generated plan. The MCO plan was randomly selected out of the set of MCO plans that had target coverage similar to the clinically generated plan and the same or better sparing of the organs at risk (OAR). Additionally, a validation study of the AWCM method vs the weighted sum method was performed.

RESULTS

The adaptive MCO algorithm generated Pareto points on the Pareto hypersurface in a fast (2-3 hr) and efficient manner for 2 cases with 10 and 18 objectives. The MCO algorithm generated a dense and well-distributed set of Pareto points on the objective space, and was able to achieve minimization of the Pareto sampling metric. The selected MCO plan showed an improvement of the DVH objectives in comparison to the clinically optimized plan in both cases. For case one, the MCO plan showed a 48% reduction of the 50% dose to OARs and a 16% reduction of the 1% dose to OARs. For case 2, the MCO plan showed a 72% reduction of the 50% dose to OARs and a 42% reduction of the 1% dose to OARs. The comparison of AWCM to WS showed that the AWCM method has a dosimetric advantage over WS for both patient cases.

CONCLUSION

We introduced an adaptive MCO algorithm for IMPT accelerated using GPUs. The algorithm is based on an adaptive method for generating Pareto plans in the objective space. We have shown that the algorithm can provide rapid and efficient mapping of the multicriteria Pareto surface with clinically deliverable plans.

Automation in intensity modulated radiotherapy treatment planning-a review of recent innovations

Radiotherapy treatment planning of complex radiotherapy techniques, such as intensity modulated radiotherapy and volumetric modulated arc therapy, is a resource-intensive process requiring a high level of treatment planner intervention to ensure high plan quality. This can lead to variability in the quality of treatment plans and the efficiency in which plans are produced, depending on the skills and experience of the operator and available planning time. Within the last few years, there has been significant progress in the research and development of intensity modulated radiotherapy treatment planning approaches with automation support, with most commercial manufacturers now offering some form of solution. There is a rapidly growing number of research articles published in the scientific literature on the topic. This paper critically reviews the body of publications up to April 2018. The review describes the different types of automation algorithms, including the advantages and current limitations. Also included is a discussion on the potential issues with routine clinical implementation of such software, and highlights areas for future research.

Target tailoring and proton beam therapy to reduce small bowel dose in cervical cancer radiotherapy : A comparison of benefits

PURPOSE

The aim of the study was to investigate the potential clinical benefit from both target tailoring by excluding the tumour-free proximal part of the uterus during image-guided adaptive radiotherapy (IGART) and improved dose conformity based on intensity-modulated proton therapy (IMPT).

METHODS

The study included planning CTs from 11 previously treated patients with cervical cancer with a >4-cm tumour-free part of the proximal uterus on diagnostic magnetic resonance imaging (MRI). IGART and robustly optimised IMPT plans were generated for both conventional target volumes and for MRI-based target tailoring (where the non-invaded proximal part of the uterus was excluded), yielding four treatment plans per patient. For each plan, the V15Gy, V30Gy, V45Gy and Dmean for bladder, sigmoid, rectum and bowel bag were compared, and the normal tissue complication probability (NTCP) for ≥grade 2 acute small bowel toxicity was calculated.

RESULTS

Both IMPT and MRI-based target tailoring resulted in significant reductions in V15Gy, V30Gy, V45Gy and Dmean for bladder and small bowel. IMPT reduced the NTCP for small bowel toxicity from 25% to 18%; this was further reduced to 9% when combined with MRI-based target tailoring. In four of the 11 patients (36%), NTCP reductions of >10% were estimated by IMPT, and in six of the 11 patients (55%) when combined with MRI-based target tailoring. This >10% NTCP reduction was expected if the V45Gy for bowel bag was >275 cm3 and >200 cm3, respectively, during standard IGART alone.

CONCLUSIONS

In patients with cervical cancer, both proton therapy and MRI-based target tailoring lead to a significant reduction in the dose to surrounding organs at risk and small bowel toxicity.

Impact of bowel gas and body outline variations on total accumulated dose with intensity-modulated proton therapy in locally advanced cervical cancer patients

INTRODUCTION

Density changes occurring during fractionated radiotherapy in the pelvic region may degrade proton dose distributions. The aim of the study was to quantify the dosimetric impact of gas cavities and body outline variations.

MATERIAL AND METHODS

Seven patients with locally advanced cervical cancer (LACC) were analyzed through a total of 175 daily cone beam computed tomography (CBCT) scans. Four-beams intensity-modulated proton therapy (IMPT) dose plans were generated targeting the internal target volume (ITV) composed of: primary tumor, elective and pathological nodes. The planned dose was 45 Gy [Relative-Biological-Effectiveness-weighted (RBE)] in 25 fractions and simultaneously integrated boosts of pathologic lymph nodes were 55-57.5 Gy (RBE). In total, 475 modified CTs were generated to evaluate the effect of: 1/gas cavities, 2/outline variations and 3/the two combined. The anatomy of each fraction was simulated by propagating gas cavities contours and body outlines from each daily CBCT to the pCT. Hounsfield units corresponding to gas and fat were assigned to the propagated contours. D98 (least dose received by the hottest 98% of the volume) and D99.9 for targets and V43Gy(RBE) (volume receiving ≥43 Gy(RBE)) for organs at risk (OARs) were recalculated on each modified CT, and total dose was evaluated through dose volume histogram (DVH) addition across all fractions.

RESULTS

Weight changes during radiotherapy were between -3.1% and 1.2%. Gas cavities and outline variations induced a median [range] dose degradation for ITV45 of 1.0% [0.5-3.5%] for D98 and 2.1% [0.8-6.4%] for D99.9. Outline variations had larger dosimetric impact than gas cavities. Worst nodal dose degradation was 2.0% for D98 and 2.3% for D99.9. The impact on bladder, bowel and rectum was limited with V43Gy(RBE) variations ≤3.5 cm³.

CONCLUSION

Bowel gas cavities and outline variations had minor impact on accumulated dose in targets and OAR of four-field IMPT in a LACC population of moderate weight changes.

Potential proton and photon dose degradation in advanced head and neck cancer patients by intratherapy changes

Purpose

Evaluation of dose degradation by anatomic changes for head‐and‐neck cancer (HNC) intensity‐modulated proton therapy (IMPT) relative to intensity‐modulated photon therapy (IMRT) and identification of potential indicators for IMPT treatment plan adaptation.

Methods

For 31 advanced HNC datasets, IMPT and IMRT plans were recalculated on a computed tomography scan (CT) taken after about 4 weeks of therapy. Dose parameter changes were determined for the organs at risk (OARs) spinal cord, brain stem, parotid glands, brachial plexus, and mandible, for the clinical target volume (CTV) and the healthy tissue outside planning target volume (PTV). Correlation of dose degradation with target volume changes and quality of rigid CT matching was investigated.

Results

Recalculated IMPT dose distributions showed stronger degradation than the IMRT doses. OAR analysis revealed significant changes in parotid median dose (IMPT) and near maximum dose (D_1ml) of spinal cord (IMPT, IMRT) and mandible (IMPT). OAR dose parameters remained lower in IMPT cases. CTV coverage (V_95%) and overdose (V_107%) deteriorated for IMPT plans to (93.4 ± 5.4)% and (10.6 ± 12.5)%, while those for IMRT plans remained acceptable. Recalculated plans showed similarly decreased PTV conformity, but considerable hotspots, also outside the PTV, emerged in IMPT cases. Lower CT matching quality was significantly correlated with loss of PTV conformity (IMPT, IMRT), CTV homogeneity and coverage (IMPT). Target shrinkage correlated with increased dose in brachial plexus (IMRT, IMPT), hotspot generation outside the PTV (IMPT) and lower PTV conformity (IMRT).

Conclusions

The study underlines the necessity of precise positioning and monitoring of anatomy changes, especially in IMPT which might require adaptation more often. Since OAR doses remained typically below constraints, IMPT plan adaptation will be indicated by target dose degradations.

Dosimetric advantages of proton therapy compared with photon therapy using an adaptive strategy in cervical cancer

Background Image-guided adaptive proton therapy (IGAPT) can potentially be applied to take into account interfraction motion while limiting organ at risk (OAR) dose in cervical cancer radiation therapy (RT). In this study, the potential dosimetric advantages of IGAPT compared with photon-based image-guided adaptive RT (IGART) were investigated. Material and methods For 13 cervical cancer patients, full and empty bladder planning computed tomography (CT) images and weekly CTs were acquired. Based on both primary clinical target volumes (pCTVs) [i.e. gross tumor volume (GTV), cervix, corpus-uterus and upper part of the vagina] on planning CTs, the pretreatment observed full range primary internal target volume (pITV) was interpolated to derive pITV subranges. Given corresponding ITVs (i.e. pITVs including lymph nodes), patient-specific photon and proton plan libraries were generated. Using all weekly CTs, IGART and IGAPT treatments were simulated by selecting library plans and recalculating the dose. For each recalculated IGART and IGAPT fraction, CTV (i.e. pCTV including lymph nodes) coverage was assessed and differences in fractionated substitutes of dose-volume histogram (DVH) parameters (V15Gy, V30Gy, V45Gy, Dmean, D2cc) for bladder, bowel and rectum were tested for significance (Wilcoxon signed-rank test). Also, differences in toxicity-related DVH parameters (rectum V30Gy, bowel V45Gy) were approximated based on accumulated dose distributions. Results In 92% (96%) of all recalculated IGAPT (IGART) fractions adequate CTV coverage (V95% >98%) was obtained. All dose parameters for bladder, bowel and rectum, except the fractionated substitute for rectum V45Gy, were improved using IGAPT. Also, IGAPT reduced the mean dose to bowel, bladder and rectum significantly (p < 0.01). In addition, an average decrease of rectum V30Gy and bowel V45Gy indicated reductions in toxicity probabilities when using IGAPT. Conclusion This study demonstrates the feasibility of IGAPT in cervical cancer using a plan-library based plan-of-the-day approach. Compared to photon-based IGART, IGAPT maintains target coverage while significant dose reductions for the bladder, bowel and rectum can be achieved.