A linear optimization model for high dose rate brachytherapy using a novel distance metric

Purpose.We propose a linear network-based optimization model (LNBM) for high dose rate brachytherapy (HDR-BT) that uses a novel distance metric to measure the discrepancy between the dose delivered and the prescription. Unlike models in the literature, LNBM takes advantage of the adjacency structure of the patients' voxels by formalizing them into a network.Methods.We apply LNBM to a set of 7 cervical cancer cases treated with HDR-BT. State-of-the-art commercial optimization software solves LNBM to global optimality. The results of LNBM are compared with those of inverse planning by simulated annealing (IPSA) based on tumor coverage, dosimetric indices for the critical organs at risk (OARs), isodose contour plots, and two metrics of homogeneity new to this work (hot-spots volumes and diameters).Results.LNBM produces plans with improved tumor coverage and with improved isodose contour plots and dosimetric indices for OARs that receive highest dose (bladder and rectum in this study) when compared with IPSA. Using new metrics of homogeneity, we also demonstrate that LNBM produces more homogeneous plans on these cases. An analysis of the solutions of LNBM shows that they use a significant part of the voxel network structure, providing evidence that the plans produced are different from those created using traditional penalty approaches and are more directly guided by the geometry of the patients' anatomy.Conclusions.The proposed linear network-based optimization model efficiently generates more homogeneous high quality treatment plans for HDR-BT.

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.

A new optimization algorithm for HDR brachytherapy that improves DVH-based planning: Truncated Conditional Value-at-Risk (TCVaR)

Purpose:To introduce a new optimization algorithm that improves DVH results and is designed for the type of heterogeneous dose distributions that occur in brachytherapy.Methods:The new optimization algorithm is based on a prior mathematical approach that uses mean doses of the DVH metric tails. The prior mean dose approach is referred to as conditional value-at-risk (CVaR), and unfortunately produces noticeably worse DVH metric results than gradient-based approaches. We have improved upon the CVaR approach, using the so-called Truncated CVaR (TCVaR), by excluding the hottest or coldest voxels in the structure from the calculations of the mean dose of the tail. Our approach applies an iterative sequence of convex approximations to improve the selection of the excluded voxels. Data Envelopment Analysis was used to quantify the sensitivity of TCVaR results to parameter choice and to compare the quality of a library of 256 TCVaR plans created for each of prostate, breast, and cervix treatment sites with commercially-generated plans.Results:In terms of traditional DVH metrics, TCVaR outperformed CVaR and the improvements increased monotonically as more iterations were used to identify and exclude the hottest/coldest voxels from the optimization problem. TCVaR also outperformed the Eclipse-Brachyvision TPS, with an improvement in PTVD95% (for equivalent organ-at-risk doses) of up to 5% (prostate), 3% (breast), and 1% (cervix).Conclusions:A novel optimization algorithm for HDR treatment planning produced plans with superior DVH metrics compared with a prior convex optimization algorithm as well as Eclipse-Brachyvision. The algorithm is computationally efficient and has potential applications as a primary optimization algorithm or quality assurance for existing optimization approaches.

An extended dose-volume model in high dose-rate brachytherapy - Using mean-tail-dose to reduce tumor underdosage

Purpose

High dose–rate brachytherapy is a method of radiotherapy for cancer treatment in which the radiation source is placed within the body. In addition to give a high enough dose to a tumor, it is also important to spare nearby healthy organs [organs at risk (OAR)]. Dose plans are commonly evaluated using the so‐called dosimetric indices; for the tumor, the portion of the structure that receives a sufficiently high dose is calculated, while for OAR it is instead the portion of the structure that receives a sufficiently low dose that is of interest. Models that include dosimetric indices are referred to as dose–volume models (DVMs) and have received much interest recently. Such models do not take the dose to the coldest (least irradiated) volume of the tumor into account, which is a distinct weakness since research indicates that the treatment effect can be largely impaired by tumor underdosage even to small volumes. Therefore, our aim is to extend a DVM to also consider the dose to the coldest volume.

Methods

An improved DVM for dose planning is proposed. In addition to optimizing with respect to dosimetric indices, this model also takes mean dose to the coldest volume of the tumor into account.

Results

Our extended model has been evaluated against a standard DVM in ten prostate geometries. Our results show that the dose to the coldest volume could be increased, while also computing times for the dose planning were improved.

Conclusion

While the proposed model yields dose plans similar to other models in most aspects, it fulfils its purpose of increasing the dose to cold tumor volumes. An additional benefit is shorter solution times, and especially for clinically relevant times (of minutes) we show major improvements in tumour dosimetric indices.