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.

Applying the column generation method to the intensity modulated high dose rate brachytherapy inverse planning problem.

Objective.Intensity modulated high dose rate brachytherapy (IMBT) is a rapidly developing application of brachytherapy where anisotropic dose distributions can be produced at each source dwell position. This technique is made possible by placing rotating metallic shields inside brachytherapy needles or catheters. By dynamically directing the radiation towards the tumours and away from the healthy tissues, a more conformal dose distribution can be obtained. The resulting treatment planning involves optimizing dwell position and shield angle (DPSA). The aim of this study was to investigate the column generation method for IMBT treatment plan optimization.Approach.A column generation optimization algorithm was developed to optimize the dwell times and shield angles. A retrospective study was performed on 10 prostate cases using RapidBrachyMCTPS. At every iteration, the plan was optimized with the chosen DPSA which would best improve the cost function that was added to the plan. The optimization process was stopped when the remaining DPSAs would not add value to the plan to limit the plan complexity.Main results.The average number of DPSAs and voxels were 2270 and 7997, respectively. The column generation approach yielded near-optimal treatment plans by using only 11% of available DPSAs on average in ten prostate cases. The coverage and organs at risk constraints passed in all ten cases.Significance.The column generation method produced high-quality deliverable prostate IMBT plans. The treatment plan quality reached a plateau, where adding more DPSAs had a minimal effect on dose volume histogram parameters. The iterative nature of the column generation method allows early termination of the treatment plan creation process as soon as the dosimetric indices from dose volume histogram satisfy the clinical requirements or if their values stabilize.

Dosimetric impact of a robust optimization approach to mitigate effects from rotational uncertainty in prostate intensity-modulated brachytherapy.

BACKGROUND: Intensity-modulated brachytherapy (IMBT) is an emerging technology for cancer treatment, in which radiation sources are shielded to shape the dose distribution. The rotatable shields provide an additional degree of freedom, but also introduce an additional, directional, type of uncertainty, compared to conventional high-dose-rate brachytherapy (HDR BT). PURPOSE: We propose and evaluate a robust optimization approach to mitigate the effects of rotational uncertainty in the shields with respect to planning criteria. METHODS: A previously suggested prototype for platinum-shielded prostate 169 Yb-based dynamic IMBT is considered. We study a retrospective patient data set (anatomical contours and catheter placement) from two clinics, consisting of six patients that had previously undergone conventional 192 Ir HDR BT treatment. The Monte Carlo-based treatment planning software RapidBrachyMCTPS is used for dose calculations. In our computational experiments, we investigate systematic rotational shield errors of ±10° and ±20°, and the same systematic error is applied to all dwell positions in each scenario. This gives us three scenarios, one nominal and two with errors. The robust optimization approach finds a compromise between the average and worst-case scenario outcomes. RESULTS: We compare dose plans obtained from standard models and their robust counterparts. With dwell times obtained from a linear penalty model (LPM), for 10° errors, the dose to urethra ( D 0.1 c c $D_{0.1cc}$ ) and rectum ( D 0.1 c c $D_{0.1cc}$ and D 1 c c $D_{1cc}$ ) increase with up to 5% and 7%, respectively, in the worst-case scenario, while with the robust counterpart, the corresponding increases were 3% and 3%. For all patients and all evaluated criteria, the worst-case scenario outcome with the robust approach had lower deviation compared to the standard model, without compromising target coverage. We also evaluated shield errors up to 20° and while the deviations increased to a large extent with the standard models, the robust models were capable of handling even such large errors. CONCLUSIONS: We conclude that robust optimization can be used to mitigate the effects from rotational uncertainty and to ensure the treatment plan quality of IMBT.

A novel minimally invasive dynamic-shield, intensity-modulated brachytherapy system for the treatment of cervical cancer.

PURPOSE: To present a novel, MRI-compatible dynamicshield intensity modulated brachytherapy (IMBT) applicator and delivery system using 192 Ir, 75 Se, and 169 Yb radioisotopes for the treatment of locally advanced cervical cancer. Needle-free IMBT is a promising technique for improving target coverage and organs at risk (OAR) sparing. METHODS AND MATERIALS: The IMBT delivery system dynamically controls the rotation of a novel tungsten shield placed inside an MRI-compatible, 6-mm wide intrauterine tandem. Using 36 cervical cancer cases, conventional intracavitary brachytherapy (IC-BT) and intracavitary/interstitial brachytherapy (IC/IS-BT) (10Ci 192 Ir) plans were compared to IMBT (10Ci 192 Ir; 11.5Ci 75 Se; 44Ci 169 Yb). All plans were generated using the Geant4-based Monte Carlo dose calculation engine, RapidBrachyMC. Treatment plans were optimized then normalized to the same high-risk clinical target volume (HR-CTV) D90 and the D2cc for bladder, rectum, and sigmoid in the research brachytherapy planning system, RapidBrachyMCTPS. Plans were renormalized until either of the three OAR reached dose limits to calculate the maximum achievable HR-CTV D90 and D98 . RESULTS: Compared to IC-BT, IMBT with either of the three radionuclides significantly improves the HR-CTV D90 and D98 by up to 5.2% ± 0.3% (P < 0.001) and 6.7% ± 0.5% (P < 0.001), respectively, with the largest dosimetric enhancement when using 169 Yb followed by 75 Se and then 192 Ir. Similarly, D2cc for all OAR improved with IMBT by up to 7.7% ± 0.6% (P < 0.001). For IC/IS-BT cases, needle-free IMBT achieved clinically acceptable plans with 169 Yb-based IMBT further improving HR-CTV D98 by 1.5% ± 0.2% (P = 0.034) and decreasing sigmoid D2cc by 1.9% ± 0.4% (P = 0.048). Delivery times for IMBT are increased by a factor of 1.7, 3.3, and 2.3 for 192 Ir, 75 Se, and 169 Yb, respectively, relative to conventional 192 Ir BT. CONCLUSIONS: Dynamic shield IMBT provides a promising alternative to conventional IC- and IC/IS-BT techniques with significant dosimetric enhancements and even greater improvements with intermediate energy radionuclides. The ability to deliver a highly conformal, OAR-sparing dose without IS needles provides a simplified method for improving the therapeutic ratio less invasively and in a less resource intensive manner.

Fast mixed integer optimization (FMIO) for high dose rate brachytherapy.

The purpose of this work was to develop an efficient quadratic mixed integer programming algorithm for high dose rate (HDR) brachytherapy treatment planning problems and integrate the algorithm into an open-source Monte Carlo based treatment planning software, RapidBrachyMCTPS. The mixed-integer algorithm yields a globally optimum solution to the dose volume histogram (DVH) based problem and, unlike other methods, is not susceptible to local minimum trapping. A hybrid linear-quadratic penalty model coupled to a mixed integer programming model was used to optimize treatment plans for 10 prostate cancer patients. Dose distributions for each dwell position were calculated with RapidBrachyMCTPS with type A uncertainties less than 0.2% in voxels within the planning target volume (PTV). The optimization process was divided into two parts. First, the data was preprocessed, in which the problem size was reduced by eliminating voxels that had negligible impact on the solution (e.g. far from the dwell position). Second, the best combination of dwell times to obtain a plan with the highest score was found. The dwell positions and dose volume constraints were used as input to a commercial mixed integer optimizer (Gurobi Optimization, Inc.). A penalty-based criterion was adopted for the scoring. The voxel-reduction technique successfully reduced the problem size by an average of 91%, without loss of quality. The preprocessing of the optimization process required on average 4 s and solving for the global maximum required on average 33 s. The total optimization time averaged 37 s, which is a substantial improvement over the ∼15 min optimization time reported in published literature. The plan quality was evaluated by evaluating dose volume metrics, including PTV D90, rectum and bladder D1cc and urethra D0.1cc. In conclusion, fast mixed integer optimization is an order of magnitude faster than current mixed-integer approaches for solving HDR brachytherapy treatment planning problems with DVH based metrics.