A multi-criteria optimization approach for HDR prostate brachytherapy: I. Pareto surface approximation

Penalty weight tuning in high dose rate brachytherapy using multi-objective Bayesian optimization

Objective.Treatment plan optimization in high dose rate brachytherapy often requires manual fine-tuning of penalty weights for each objective, which can be time-consuming and dependent on the planner's experience. To automate this process, this study used a multi-criteria approach called multi-objective Bayesian optimization with q-noisy expected hypervolume improvement as its acquisition function (MOBO-qNEHVI).Approach.The treatment plans of 13 prostate cancer patients were retrospectively imported to a research treatment planning system, RapidBrachyMTPS, where fast mixed integer optimization (FMIO) performs dwell time optimization given a set of penalty weights to deliver 15 Gy to the target volume. MOBO-qNEHVI was used to find patient-specific Pareto optimal penalty weight vectors that yield clinically acceptable dose volume histogram metrics. The relationship between the number of MOBO-qNEHVI iterations and the number of clinically acceptable plans per patient (acceptance rate) was investigated. The performance time was obtained for various parameter configurations.Main results.MOBO-qNEHVI found clinically acceptable treatment plans for all patients. With increasing the number of MOBO-qNEHVI iterations, the acceptance rate grew logarithmically while the performance time grew exponentially. Fixing the penalty weight of the tumour volume to maximum value, adding the target dose as a parameter, initiating MOBO-qNEHVI with 25 parallel sampling of FMIO, and running 6 MOBO-qNEHVI iterations found solutions that delivered 15 Gy to the hottest 95% of the clinical target volume while respecting the dose constraints to the organs at risk. The average acceptance rate for each patient was 89.74% ± 8.11%, and performance time was 66.6 ± 12.6 s. The initiation took 22.47 ± 7.57 s, and each iteration took 7.35 ± 2.45 s to find one Pareto solution.Significance.MOBO-qNEHVI combined with FMIO can automatically explore the trade-offs between treatment plan objectives in a patient specific manner within a minute. This approach can reduce the dependency of plan quality on planner's experience and reduce dose to the organs at risk.

Catheters and dose optimization using a modified CVT algorithm and multi-criteria optimization in prostate HDR brachytherapy

Currently, in HDR brachytherapy planning, the catheter's positions are often selected by the planner which involves the planner's experience. The catheters are then inserted using a template which helps to guide the catheters. For certain applications, it is of interest to choose the optimal location and number of catheters needed for dose coverage and potential decrease of the treatment's toxicity. Hence, it is of great importance to develop patient-specific algorithms for catheters and dose optimization. A modified Centroidal Voronoi tessellation (CVT) algorithm is implemented and merged with a GPU-based multi-criteria optimization algorithm (gMCO). The CVT algorithm optimizes the catheters' positions, and the gMCO algorithm optimizes the dwell times and dwell positions. The CVT algorithm can be used simultaneously for insertion with or without a template. Some improvements to the CVT algorithm are presented such as a new way of considering the area that needs to be covered. One hundred and eight previously treated prostates HDR cases using real-time ultrasound (US) are used to evaluate the different optimization procedures. The plan robustness is evaluated using two types of errors; deviations (random) in the insertion and deviation (systematic) in the reconstruction of the catheters. Using gMCO on clinically inserted catheter increases the acceptance rate by 37% for RTOG criteria. Our results show that all the patients respect RTOG criteria with 11 catheters using CVT+gMCO with a template of 5 mm. The number of catheters needed for all patients to respect RTOG criteria with the freehand technique is 10 catheters using CVT+gMCO. When deviations are introduced, using a template, the acceptance rate goes to 85% with 3 mm deviations using 11 catheters. This decrease is less significant when the number of catheters is higher, decreasing by less than 5% with a 3 mm deviation using 13 catheters or more. In conclusion, it is feasible to decrease the number of catheters needed to treat most patients. Some cases still need a high number of catheters to reach the plan's criteria. Using gMCO allows an increase in the plan quality while using CVT reduces the number of catheters. A higher number of catheters equates to plans that are more robust to deviations. This article is protected by copyright. All rights reserved.

Emerging technologies in brachytherapy

Brachytherapy is a mature treatment modality. The literature is abundant in terms of review articles and comprehensive books on the latest established as well as evolving clinical practices. The intent of this article is to part ways and look beyond the current state-of-the-art and review emerging technologies that are noteworthy and perhaps may drive the future innovations in the field. There are plenty of candidate topics that deserve a deeper look, of course, but with practical limits in this communicative platform, we explore four topics that perhaps is worthwhile to review in detail at this time. First, intensity modulated brachytherapy (IMBT) is reviewed. The IMBT takes advantage ofanisotropicradiation profile generated through intelligent high-density shielding designs incorporated onto sources and applicators such to achieve high quality plans. Second, emerging applications of 3D printing (i.e. additive manufacturing) in brachytherapy are reviewed. With the advent of 3D printing, interest in this technology in brachytherapy has been immense and translation swift due to their potential to tailor applicators and treatments customizable to each individual patient. This is followed by, in third, innovations in treatment planning concerning catheter placement and dwell times where new modelling approaches, solution algorithms, and technological advances are reviewed. And, fourth and lastly, applications of a new machine learning technique, called deep learning, which has the potential to improve and automate all aspects of brachytherapy workflow, are reviewed. We do not expect that all ideas and innovations reviewed in this article will ultimately reach clinic but, nonetheless, this review provides a decent glimpse of what is to come. It would be exciting to monitor as IMBT, 3D printing, novel optimization algorithms, and deep learning technologies evolve over time and translate into pilot testing and sensibly phased clinical trials, and ultimately make a difference for cancer patients. Today's fancy is tomorrow's reality. The future is bright for brachytherapy.

Optimization in treatment planning of high dose-rate brachytherapy - Review and analysis of mathematical models

Treatment planning in high dose-rate brachytherapy has traditionally been conducted with manual forward planning, but inverse planning is today increasingly used in clinical practice. There is a large variety of proposed optimization models and algorithms to model and solve the treatment planning problem. Two major parts of inverse treatment planning for which mathematical optimization can be used are the decisions about catheter placement and dwell time distributions. Both these problems as well as integrated approaches are included in this review. The proposed models include linear penalty models, dose-volume models, mean-tail dose models, quadratic penalty models, radiobiological models, and multiobjective models. The aim of this survey is twofold: (i) to give a broad overview over mathematical optimization models used for treatment planning of brachytherapy and (ii) to provide mathematical analyses and comparisons between models. New technologies for brachytherapy treatments and methods for treatment planning are also discussed. Of particular interest for future research is a thorough comparison between optimization models and algorithms on the same dataset, and clinical validation of proposed optimization approaches with respect to patient outcome.

An Inverse Dose Optimization Algorithm for Three-Dimensional Brachytherapy

PURPOSE

To investigate an implementation method and the results of an inverse dose optimization algorithm, Gradient Based Planning Optimization (GBPO), for three-dimensional brachytherapy.

METHODS

The GBPO used a quadratic objective function, and a dwell time modulation item was added to the objective function to restrict the dwell time variance. We retrospectively studied 4 cervical cancer patients using different applicators and 15 cervical cancer patients using the Fletcher applicator. We assessed the plan quality of GBPO by isodose lines for the patients using different applicators. For the 15 patients using the Fletcher applicator, we utilized dose-volume histogram (DVH) parameters of HR-CTV (D100%, V150%) and organs at risk (OARs) (D0.1cc, D1cc, D2cc) to evaluate the difference between the GBPO plans and the IPSA (Inverse Planning Simulated Annealing) plans, as well as the GBPO plans and the Graphic plans.

RESULTS

For the 4 patients using different applicators, the dose distributions are conformable. For the 15 patients using the Fletcher applicator, when the dwell time modulation factor (DTMF) is less than 20, the dwell time deviation reduces quickly; however, after the DTMF increased to 100, the dwell time deviation has no remarkable change. The difference in dosimetric parameters between the GBPO plans and the IPSA plans is not statistically significant (P>0.05). The GBPO plans have a higher D100% (3.57 ± 0.36, 3.38 ± 0.34; P<0.01) and a lower V150% (55.73 ± 4.06, 57.75 ± 3.79; P<0.01) than those of the Graphic plans. The differences in other DVH parameters are negligible between the GBPO plans and the Graphic plans.

CONCLUSIONS

The GBPO plans have a comparable quality as the IPSA plans and the Graphic plans for the studied cervical cancer cases. The GBPO algorithm could be integrated into a three-dimensional brachytherapy treatment planning system after studying more sites.

Bi-objective optimization of catheter positions for high-dose-rate prostate brachytherapy

Purpose

Bi‐objective simultaneous optimization of catheter positions and dwell times for high‐dose‐rate (HDR) prostate brachytherapy, based directly on dose‐volume indices, has shown promising results. However, optimization with the state‐of‐the‐art evolutionary algorithm MO‐RV‐GOMEA so far required several hours of runtime, and resulting catheter positions were not always clinically feasible. The aim of this study is to extend the optimization model and apply GPU parallelization to achieve clinically acceptable computation times. The resulting optimization procedure is compared with a previously introduced method based solely on geometric criteria, the adapted Centroidal Voronoi Tessellations (CVT) algorithm.

Methods

Bi‐objective simultaneous optimization was performed with a GPU‐parallelized version of MO‐RV‐GOMEA. This optimization of catheter positions and dwell times was retrospectively applied to the data of 26 patients previously treated with HDR prostate brachytherapy for 8–16 catheters (steps of 2). Optimization of catheter positions using CVT was performed in seconds, after which optimization of only the dwell times using MO‐RV‐GOMEA was performed in 1 min.

Results

Simultaneous optimization of catheter positions and dwell times using MO‐RV‐GOMEA was performed in 5 min. For 16 down to 8 catheters (steps of 2), MO‐RV‐GOMEA found plans satisfying the planning‐aims for 20, 20, 18, 14, and 11 out of the 26 patients, respectively. CVT achieved this for 19, 17, 13, 9, and 2 patients, respectively. The P‐value for the difference between MO‐RV‐GOMEA and CVT was 0.023 for 16 catheters, 0.005 for 14 catheters, and <0.001 for 12, 10, and 8 catheters.

Conclusions

With bi‐objective simultaneous optimization on a GPU, high‐quality catheter positions can now be obtained within 5 min, which is clinically acceptable, but slower than CVT. For 16 catheters, the difference between MO‐RV‐GOMEA and CVT is clinically irrelevant. For 14 catheters and less, MO‐RV‐GOMEA outperforms CVT in finding plans satisfying all planning‐aims.

Evaluating the impact of real-time multicriteria optimizers integrated with interactive plan navigation tools for HDR brachytherapy

PURPOSE

Currently in high-dose-rate (HDR) brachytherapy planning, manual fine-tuning of an objective function is a common practice. Furthermore, automated planning approaches such as multicriteria optimization (MCO) are still limited to the automatic generation of a single treatment plan. This study aims to quantify planning efficiency gains when using a graphics processing unit-based MCO (gMCO) algorithm combined with a novel graphical user interface (gMCO-GUI) that integrates efficient automated and interactive plan navigation tools.

METHODS AND MATERIALS

The gMCO algorithm was used to generate 1000 Pareto optimal plans per case for 379 prostate cases. gMCO-GUI was developed to allow plan navigation through all plans. gMCO-GUI integrates interactive parameter selection tools directly with the optimization algorithm to allow plan navigation. The quality of each plan was evaluated based on the Radiation Treatment Oncology Group 0924 protocol and a more stringent institutional protocol (INSTp). gMCO-GUI allows real-time time display of the dose-volume histogram indices, the dose-volume histogram curves, and the isodose lines during the plan navigation.

RESULTS

Over the 379 cases, the fraction of Radiation Treatment Oncology Group 0924 protocol valid plans with target coverage greater than 95% was 90.8%, compared with 66.0% for clinical plans. The fraction of INSTp valid plans with target coverage greater than 95% was 81.8%, compared with 62.3% for clinical plans. The average time to compute 1000 deliverable plans with gMCO was 12.5 s, including the full computation of the 3D dose distributions.

CONCLUSIONS

Combining the gMCO algorithm with automated and interactive plan navigation tools resulted in simultaneous gains in both plan quality and planning efficiency.