Multi-objective optimization in radiotherapy: applications to stereotactic radiosurgery and prostate brachytherapy

Artificial intelligence applications in brachytherapy: A literature review

PURPOSE

Artificial intelligence (AI) has the potential to simplify and optimize various steps of the brachytherapy workflow, and this literature review aims to provide an overview of the work done in this field.

METHODS AND MATERIALS

We conducted a literature search in June 2022 on PubMed, Embase, and Cochrane for papers that proposed AI applications in brachytherapy.

RESULTS

A total of 80 papers satisfied inclusion/exclusion criteria. These papers were categorized as follows: segmentation (24), registration and image processing (6), preplanning (13), dose prediction and treatment planning (11), applicator/catheter/needle reconstruction (16), and quality assurance (10). AI techniques ranged from classical models such as support vector machines and decision tree-based learning to newer techniques such as U-Net and deep reinforcement learning, and were applied to facilitate small steps of a process (e.g., optimizing applicator selection) or even automate the entire step of the workflow (e.g., end-to-end preplanning). Many of these algorithms demonstrated human-level performance and offer significant improvements in speed.

CONCLUSIONS

AI has potential to augment, automate, and/or accelerate many steps of the brachytherapy workflow. We recommend that future studies adhere to standard reporting guidelines. We also stress the importance of using larger sample sizes and reporting results using clinically interpretable measures.

Bi-criteria Pareto optimization to balance irradiation time and dosimetric objectives in proton arc therapy

Objective. Proton arc therapy (PAT) is a new delivery technique that exploits the continuous rotation of the gantry to distribute the therapeutic dose over many angular windows instead of using a few static fields, as in conventional (intensity-modulated) proton therapy. Although coming along with many potential clinical and dosimetric benefits, PAT has also raised a new optimization challenge. In addition to the dosimetric goals, the beam delivery time (BDT) needs to be considered in the objective function. Considering this bi-objective formulation, the task of finding a good compromise with appropriate weighting factors can turn out to be cumbersome.Approach. We have computed Pareto-optimal plans for three disease sites: a brain, a lung, and a liver, following a method of iteratively choosing weight vectors to approximate the Pareto front with few points. Mixed-integer programming (MIP) was selected to state the bi-criteria PAT problem and to find Pareto optimal points with a suited solver.Main results. The trade-offs between plan quality and beam irradiation time (staticBDT) are investigated by inspecting three plans from the Pareto front. The latter are carefully picked to demonstrate significant differences in dose distribution and delivery time depending on their location on the frontier. The results were benchmarked against IMPT and SPArc plans showing the strength of degrees of freedom coming along with MIP optimization.Significance. This paper presents for the first time the application of bi-criteria optimization to the PAT problem, which eventually permits the planners to select the best treatment strategy according to the patient conditions and clinical resources available.

Multiobjective optimization challenges in perioperative anesthesia: A review

Physicians use perioperative decision-support tools to mitigate risks and maximize benefits to achieve the most successful outcome for patients. Contemporary risk-assessment practices augment surgeons' judgement and experience with decision-support algorithms driven by big data and machine learning. These algorithms accurately assess risk for a wide range of postoperative complications by parsing large datasets and performing complex calculations that would be cumbersome for busy clinicians. Even with these advancements, large gaps in perioperative risk assessment remain; decision-support algorithms often cannot account for risk-reduction therapies applied during a patient's perioperative course and do not quantify tradeoffs between competing goals of care (eg, balancing postoperative pain control with the risk of respiratory depression or balancing intraoperative volume resuscitation with the risk for complications from pulmonary edema). Multiobjective optimization solutions have been applied to similar problems successfully but have not yet been applied to perioperative decision support. Given the large volume of data available via electronic medical records, including intraoperative data, it is now feasible to successfully apply multiobjective optimization in perioperative care. Clinical application of multiobjective optimization would require semiautomated pipelines for analytics and reporting model outputs and a careful development and validation process. Under these circumstances, multiobjective optimization has the potential to support personalized, patient-centered, shared decision-making with precision and balance.

Multicriteria plan optimization in the hands of physicians: a pilot study in prostate cancer and brain tumors

BACKGROUND

The purpose of this study was to demonstrate the feasibility of physician driven planning in intensity modulated radiotherapy (IMRT) with a multicriteria optimization (MCO) treatment planning system and template based plan optimization. Exploiting the full planning potential of MCO navigation, this alternative planning approach intends to improve planning efficiency and individual plan quality.

METHODS

Planning was retrospectively performed on 12 brain tumor and 10 post-prostatectomy prostate patients previously treated with MCO-IMRT. For each patient, physicians were provided with a template-based generated Pareto surface of optimal plans to navigate, using the beam angles from the original clinical plans. We compared physician generated plans to clinically delivered plans (created by dosimetrists) in terms of dosimetric differences, physician preferences and planning times.

RESULTS

Plan qualities were similar, however physician generated and clinical plans differed in the prioritization of clinical goals. Physician derived prostate plans showed significantly better sparing of the high dose rectum and bladder regions (p(D1) < 0.05; D1: dose received by 1% of the corresponding structure). Physicians' brain tumor plans indicated higher doses for targets and brainstem (p(D1) < 0.05). Within blinded plan comparisons physicians preferred the clinical plans more often (brain: 6:3 out of 12, prostate: 2:6 out of 10) (not statistically significant). While times of physician involvement were comparable for prostate planning, the new workflow reduced the average involved time for brain cases by 30%. Planner times were reduced for all cases. Subjective benefits, such as a better understanding of planning situations, were observed by clinicians through the insight into plan optimization and experiencing dosimetric trade-offs.

CONCLUSIONS

We introduce physician driven planning with MCO for brain and prostate tumors as a feasible planning workflow. The proposed approach standardizes the planning process by utilizing site specific templates and integrates physicians more tightly into treatment planning. Physicians' navigated plan qualities were comparable to the clinical plans. Given the reduction of planning time of the planner and the equal or lower planning time of physicians, this approach has the potential to improve departmental efficiencies.

Application programming in C# environment with recorded user software interactions and its application in autopilot of VMAT/IMRT treatment planning

An autopilot scheme of volumetric‐modulated arc therapy (VMAT)/intensity‐modulated radiation therapy (IMRT) planning with the guidance of prior knowledge is established with recorded interactions between a planner and a commercial treatment planning system (TPS). Microsoft (MS) Visual Studio Coded UI is applied to record some common planner‐TPS interactions as subroutines. The TPS used in this study is a Windows‐based Eclipse system. The interactions of our application program with Eclipse TPS are realized through a series of subroutines obtained by prerecording the mouse clicks or keyboard strokes of a planner in operating the TPS. A strategy to autopilot Eclipse VMAT/IMRT plan selection process is developed as a specific example of the proposed “scripting” method. The autopiloted planning is navigated by a decision function constructed with a reference plan that has the same prescription and similar anatomy with the case at hand. The calculation proceeds by alternating between the Eclipse optimization and the outer‐loop optimization independent of the Eclipse. In the C# program, the dosimetric characteristics of a reference treatment plan are used to assess and modify the Eclipse planning parameters and to guide the search for a clinically sensible treatment plan. The approach is applied to plan a head and neck (HN) VMAT case and a prostate IMRT case. Our study demonstrated the feasibility of application programming method in C# environment with recorded interactions of planner‐TPS. The process mimics a planner's planning process and automatically provides clinically sensible treatment plans that would otherwise require a large amount of manual trial and error of a planner. The proposed technique enables us to harness a commercial TPS by application programming via the use of recorded human computer interactions and provides an effective tool to greatly facilitate the treatment planning process.

PACS number(s): 87.55.D‐, 87.55.kd, 87.55.de

A novel greedy heuristic-based approach to intraoperative planning for permanent prostate brachytherapy

This paper presents a greedy heuristic‐based double iteration and rectification (DIR) approach to intraoperative planning for permanent prostate brachytherapy. The DIR approach adopts a greedy seed selection (GSS) procedure to obtain a preliminary plan. In this process, the potential seeds are evaluated according to their ability to irradiate target while spare organs at risk (OARs), and their impact on dosimetric homogeneity within target volume. A flexible termination condition is developed for the GSS procedure, which guarantees sufficient dose within target volume while avoids overdosing the OARs. The preliminary treatment plan generated by the GSS procedure is further refined by the double iteration (DI) and rectification procedure. The DI procedure removes the needles containing only one seed (single seed) and implements the GSS procedure again to get a temporary plan. The DI procedure terminates until the needles number of the temporary plan does not decrease. This process is guided by constantly removing undesired part rather than imposing extra constrains. Following the DI procedure, the rectification procedure attempts to replace the remaining single seeds with the acceptable ones within the existing needles. The change of dosimetric distribution (DD) after the replacement is evaluated to determine whether to accept or to withdraw the replacement. Experimental results demonstrate that the treatment plans obtained by the DIR approach caters to all clinical considerations. Compared with currently available methods, DIR approach is faster, more reliable, and more suitable for intraoperative treatment planning in the operation room.

PACS number: 87

Prior-knowledge treatment planning for volumetric arc therapy using feature-based database mining

Treatment planning for volumetric arc therapy (VMAT) is a lengthy process that requires many rounds of optimizations to obtain the best treatment settings and optimization constraints for a given patient's geometry. We propose a feature‐selection search engine that explores previously treated cases of similar anatomy, returning the optimal plan configurations and attainable DVH constraints. Using an institutional database of 83 previously treated cases of prostate carcinoma treated with volumetric‐modulated arc therapy, the search procedure first finds the optimal isocenter position with an optimization procedure, then ranks the anatomical similarity as the mean distance between targets. For the best matching plan, the planning information is reformatted to the DICOM format and imported into the treatment planning system to suggest isocenter, arc directions, MLC patterns, and optimization constraints that can be used as starting points in the optimization process. The approach was tested to create prospective treatment plans based on anatomical features that match previously treated cases from the institution database. By starting from a near‐optimal solution and using previous optimization constraints, the best matching test only required simple optimization steps to further decrease target inhomogeneity, ultimately reducing time spend by the therapist in planning arcs' directions and lengths.

PACS number: 87.55.D‐, 87.55.de

Feasibility of identification of gamma knife planning strategies by identification of pareto optimal gamma knife plans

The use of conformity indices to optimize Gamma Knife planning is common, but does not address important tradeoffs between dose to tumor and normal tissue. Pareto analysis has been used for this purpose in other applications, but not for Gamma Knife (GK) planning. The goal of this work is to use computer models to show that Pareto analysis may be feasible for GK planning to identify dosimetric tradeoffs. We define a GK plan A to be Pareto dominant to B if the prescription isodose volume of A covers more tumor but not more normal tissue than B, or if A covers less normal tissue but not less tumor than B. A plan is Pareto optimal if it is not dominated by any other plan. Two different Pareto optimal plans represent different tradeoffs between dose to tumor and normal tissue, because neither plan dominates the other. ‘GK simulator’ software calculated dose distributions for GK plans, and was called repetitively by a genetic algorithm to calculate Pareto dominant plans. Three irregular tumor shapes were tested in 17 trials using various combinations of shots. The mean number of Pareto dominant plans/trial was 59 ± 17 (sd). Different planning strategies were identified by large differences in shot positions, and 70 of the 153 coordinate plots (46%) showed differences of 5mm or more. The Pareto dominant plans dominated other nearby plans. Pareto dominant plans represent different dosimetric tradeoffs and can be systematically calculated using genetic algorithms. Automatic identification of non-intuitive planning strategies may be feasible with these methods.

A new concept for interactive radiotherapy planning with multicriteria optimization: first clinical evaluation

BACKGROUND AND PURPOSE

Currently, inverse planning for intensity-modulated radiotherapy (IMRT) can be a time-consuming trial and error process. This is because many planning objectives are inherently contradictory and cannot reach their individual optimum all at the same time. Therefore in clinical practice the potential of IMRT cannot be fully exploited for all patients. Multicriteria (multiobjective) optimization combined with interactive plan navigation is a promising approach to overcome these problems.

PATIENTS AND METHODS

We developed a new inverse planning system called "Multicriteria Interactive Radiotherapy Assistant (MIRA)". The optimization result is a database of patient specific, Pareto-optimal plan proposals. The database is explored with an intuitive user interface that utilizes both a new interactive element for plan navigation and familiar dose visualizations in form of DVH and isodose projections. Two clinical test cases, one paraspinal meningioma case and one prostate case, were optimized using MIRA and compared with the clinically approved planning program KonRad.

RESULTS

Generating the databases required no user interaction and took approx. 2-3h per case. The interactive exploration required only a few minutes until the best plan was identified, resulting in a significant reduction of human planning time. The achievable plan quality was comparable to KonRad with the additional benefit of having plan alternatives at hand to perform a sensitivity analysis or to decide for a different clinical compromise.

CONCLUSIONS

The MIRA system provides a complete database and interactive exploration of the solution space in real time. Hence, it is ideally suited for the inherently multicriterial problem of inverse IMRT treatment planning.

Optimization of HDR brachytherapy dose distributions using linear programming with penalty costs

Prostate cancer is increasingly treated with high-dose-rate (HDR) brachytherapy, a type of radiotherapy in which a radioactive source is guided through catheters temporarily implanted in the prostate. Clinicians must set dwell times for the source inside the catheters so the resulting dose distribution minimizes deviation from dose prescriptions that conform to patient-specific anatomy. The primary contribution of this paper is to take the well-established dwell times optimization problem defined by Inverse Planning by Simulated Annealing (IPSA) developed at UCSF and exactly formulate it as a linear programming (LP) problem. Because LP problems can be solved exactly and deterministically, this formulation provides strong performance guarantees: one can rapidly find the dwell times solution that globally minimizes IPSA's objective function for any patient case and clinical criteria parameters. For a sample of 20 prostates with volume ranging from 23 to 103 cc, the new LP method optimized dwell times in less than 15 s per case on a standard PC. The dwell times solutions currently being obtained clinically using simulated annealing (SA), a probabilistic method, were quantitatively compared to the mathematically optimal solutions obtained using the LP method. The LP method resulted in significantly improved objective function values compared to SA (P = 1.54 x 10(-7)), but none of the dosimetric indices indicated a statistically significant difference (P < 0.01). The results indicate that solutions generated by the current version of IPSA are clinically equivalent to the mathematically optimal solutions.

Application of influence diagrams to prostate intensity-modulated radiation therapy plan selection

The purpose is to incorporate clinically relevant factors such as patient-specific and dosimetric information as well as data from clinical trials in the decision-making process for the selection of prostate intensity-modulated radiation therapy (IMRT) plans. The approach is to incorporate the decision theoretic concept of an influence diagram into the solution of the multiobjective optimization inverse planning problem. A set of candidate IMRT plans was obtained by varying the importance factors for the planning target volume (PTV) and the organ-at-risk (OAR) in combination with simulated annealing to explore a large part of the solution space. The Pareto set for the PTV and OAR was analysed to demonstrate how the selection of the weighting factors influenced which part of the solution space was explored. An influence diagram based on a Bayesian network with 18 nodes was designed to model the decision process for plan selection. The model possessed nodes for clinical laboratory results, tumour grading, staging information, patient-specific information, dosimetric information, complications and survival statistics from clinical studies. A utility node was utilized for the decision-making process. The influence diagram successfully ranked the plans based on the available information. Sensitivity analyses were used to judge the reasonableness of the diagram and the results. In conclusion, influence diagrams lend themselves well to modelling the decision processes for IMRT plan selection. They provide an excellent means to incorporate the probabilistic nature of data and beliefs into one model. They also provide a means for introducing evidence-based medicine, in the form of results of clinical trials, into the decision-making process.

Biochemotherapy for metastatic melanoma with limited central nervous system involvement

OBJECTIVES

Biochemotherapy outcomes were examined in stage IV melanoma patients with previously treated or active central nervous system (CNS) metastases prior to systemic therapy.

PATIENTS AND METHODS

Patients who received biochemotherapy for metastatic melanoma with active or pretreated CNS metastases were compared to patients without evidence of CNS metastases in terms of response, time to progression (TTP), overall survival (OS), and treatment toxicity.

RESULTS

Twenty-six (16%) of 159 total patients began biochemotherapy with previously treated or active CNS metastases (group I), compared to 133 (84%) who were radiographically free of CNS involvement (group II). A partial or complete response to biochemotherapy was seen in 13 (50%) group I patients, compared to 56 (42%) group II patients (p = 0.243). The median TTP and median survival were 5.5 and 7.0 months, respectively, for group I patients and 6.0 and 9.9 months, respectively, for group II patients (p = 0.222 and 0.434 for TTP and OS, respectively). Five (19%) group I patients survived longer than 24 months. Gamma Knife radiosurgery or surgical resection of CNS disease prior to biochemotherapy improved survival versus delayed treatment (p = 0.017 and 0.005, respectively).

CONCLUSION

Patients with limited CNS metastases and widespread systemic disease can achieve prolonged survival with targeted treatment of CNS lesions and aggressive systemic therapy.