Hybrid proton-photon inverse optimization with uniformity-regularized proton and photon target dose

TEAM: Triangular-mEsh Adaptive and Multiscale proton spot generation method

BACKGROUND

Proton therapy is preferred for its dose conformality to spare normal tissues and organs-at-risk (OAR) via Bragg peaks with negligible exit dose. However, proton dose conformality can be further optimized: (1) the spot placement is based on the structured (e.g., Cartesian) grid, which may not offer conformal shaping to complex tumor targets; (2) the spot sampling pattern is uniform, which may be insufficient at the tumor boundary to provide the sharp dose falloff, and at the same time may be redundant at the tumor interior to provide the uniform dose coverage, for example, due to multiple Coulomb scattering (MCS); and (3) the lateral spot penumbra increases with respect to the depth due to MCS, which blurs the lateral dose falloff. On the other hand, while (1) the deliverable spots are subject to the minimum-monitor-unit (MMU) constraint, and (2) the dose rate is proportional to the MMU threshold, the current spot sampling method is sensitive to the MMU threshold and can fail to provide satisfactory plan quality for a large MMU threshold (i.e., high-dose-rate delivery).

PURPOSE

This work will develop a novel Triangular-mEsh-based Adaptive and Multiscale (TEAM) proton spot generation method to address these issues for optimizing proton dose conformality and plan delivery efficiency.

METHODS

Compared to the standard clinically-used spot placement method, three key elements of TEAM are as follows: (1) a triangular mesh instead of a structured grid: the triangular mesh is geometrically more conformal to complex target shapes and therefore more efficient and accurate for dose shaping inside and around the target; (2) adaptive sampling instead of uniform sampling: the adaptive sampling consists of relatively dense sampling at the tumor boundary to create the sharp dose falloff, which is more accurate, and coarse sampling at the tumor interior to uniformly cover the target, which is more efficient; and (3) depth-dependent sampling instead of depth-independent sampling: the depth-dependent sampling is used to compensate for MCS, that is, with increasingly dense sampling at the tumor boundary to improve dose shaping accuracy, and increasingly coarse sampling at the tumor interior to improve dose shaping efficiency, as the depth increases. In the TEAM method the spot locations are generated for each energy layer and layer-by-layer in the multiscale fashion; and then the spot weights are derived by solving the IMPT problem of dose-volume planning objectives, MMU constraints, and robustness optimization with respect to range and setup uncertainties.

RESULTS

Compared to the standard clinically-used spot placement method UNIFORM, TEAM achieved (1) better plan quality using <60% number of spots of UNIFORM; (2) better robustness to the number of spots; (3) better robustness to a large MMU threshold. Furthermore, TEAM provided better plan quality with fewer spots than other adaptive methods (Cartesian-grid or triangular-mesh).

CONCLUSIONS

A novel triangular-mesh-based proton spot placement method called TEAM is proposed, and it is demonstrated to improve plan quality, robustness to the number of spots, and robustness to the MMU threshold, compared to the clinically-used spot placement method and other adaptive methods.

Simultaneous dose and dose rate optimization via dose modifying factor modeling for FLASH effective dose

BACKGROUND

Although the FLASH radiotherapy (FLASH) can improve the sparing of organs-at-risk (OAR) via the FLASH effect, it is generally a tradeoff between the physical dose coverage and the biological FLASH coverage, for which the concept of FLASH effective dose (FED) is needed to quantify the net improvement of FLASH, compared to the conventional radiotherapy (CONV).

PURPOSE

This work will develop the first-of-its-kind treatment planning method called simultaneous dose and dose rate optimization via dose modifying factor modeling (SDDRO-DMF) for proton FLASH that directly optimizes FED.

METHODS

SDDRO-DMF models and optimizes FED using FLASH dose modifying factor (DMF) models, which can be classified into two categories: (1) the phenomenological model of the FLASH effect, such as the FLASH effectiveness model (FEM); (2) the mechanistic model of the FLASH radiobiology, such as the radiolytic oxygen depletion (ROD) model. The general framework of SDDRO-DMF will be developed, with specific DMF models using FEM and ROD, as a demonstration of general applicability of SDDRO-DMF for proton FLASH via transmission beams (TB) or Bragg peaks (BP) with single-field or multi-field irradiation. The FLASH dose rate is modeled as pencil beam scanning dose rate. The solution algorithm for solving the inverse optimization problem of SDDRO-DMF is based on iterative convex relaxation method.

RESULTS

SDDRO-DMF is validated in comparison with IMPT and a state-of-the-art method called SDDRO, with demonstrated efficacy and improvement for reducing the high dose and the high-dose volume for OAR in terms of FED. For example, in a SBRT lung case of the dose-limiting factor that the max dose of brachial plexus should be no more than 26 Gy, only SDDRO-DMF met this max dose constraint; moreover, SDDRO-DMF completely eliminated the high-dose (V70%) volume to zero for CTV10mm (a high-dose region as a 10 mm ring expansion of CTV).

CONCLUSION

We have proposed a new proton FLASH optimization method called SDDRO-DMF that directly optimizes FED using phenomenological or mechanistic models of DMF, and have demonstrated the efficacy of SDDO-DMF in reducing the high-dose volume or/and the high-dose value for OAR, compared to IMPT and a state-of-the-art method SDDRO.

A novel treatment planning method via scissor beams for uniform-target-dose proton GRID with peak-valley-dose-ratio optimization

BACKGROUND

Proton spatially fractionated RT (SFRT) can potentially synergize the unique advantages of using proton Bragg peak and SFRT peak-valley dose ratio (PVDR) to reduce the radiation-induced damage for normal tissues. Uniform-target-dose (UTD) proton GRID is a proton SFRT modality that can be clinically desirable and conveniently adopted since its UTD resembles target dose distribution in conventional proton RT (CONV). However, UTD proton GRID is not used clinically, which is likely due to the lack of an effective treatment planning method.

PURPOSE

This work will develop a novel treatment planning method using scissor beams (SB) for UTD proton GRID, with the joint optimization of PVDR and dose objectives.

METHODS

The SB method for spatial dose modulation in normal tissues with UTD has two steps: (1) a primary beam (PB) is halved with interleaved beamlets, to generate spatial dose modulation in normal tissues; (2) a complementary beam (CB) is added to fill in previously valley-dose positions in the target to generate UTD, while the CB is angled slightly from the PB, to maintain spatial dose modulation in normal tissues. A treatment planning method with PVDR optimization via the joint total variation and L1 (TVL1) regularization is developed to jointly optimize PVDR and dose objectives. The plan optimization solution is obtained using an iterative convex relaxation algorithm.

RESULTS

The new methods SB and SB-TVL1 were validated in comparison with CONV. Compared to CONV of relatively homogeneous dose distribution, SB had modulated spatial dose pattern in normal tissues with UTD and comparable plan quality. Compared to SB, SB-TVL1 further maximized PVDR, with comparable dose-volume parameters.

CONCLUSIONS

A novel SB method is proposed that can generate modulated spatial dose pattern in normal tissues to achieve UTD proton GRID. A treatment planning method with PVDR optimization capability via TVL1 regularization is developed that can jointly optimize PVDR and dose objectives for proton GRID.

A new approach to combined proton-photon therapy for metastatic cancer patients

Objective.Proton therapy is a limited resource and is typically not available to metastatic cancer patients. Combined proton-photon therapy (CPPT), where most fractions are delivered with photons and only few with protons, represents an approach to distribute proton resources over a larger patient population. In this study, we consider stereotactic radiotherapy of multiple brain or liver metastases, and develop an approach to optimally take advantage of a single proton fraction by optimizing the proton and photon dose contributions to each individual metastasis.Approach.CPPT treatments must balance two competing goals: (1) deliver a larger dose in the proton fractions to reduce integral dose, and (2) fractionate the dose in the normal tissue between metastases, which requires using the photon fractions. Such CPPT treatments are generated by simultaneously optimizing intensity modulated proton therapy (IMPT) and intensity modulated radiotherapy (IMRT) plans based on their cumulative biologically effective dose (BEDα/β). The dose contributions of the proton and photon fractions to each individual metastasis are handled as additional optimization variables in the optimization problem. The method is demonstrated for two patients with 29 and 30 brain metastases, and two patients with 4 and 3 liver metastases.Main results.Optimized CPPT plans increase the proton dose contribution to most of the metastases, while using photons to fractionate the dose around metastases which are large or located close to critical structures. On average, the optimized CPPT plans reduce the mean brain BED2by 29% and the mean liver BED4by 42% compared to IMRT-only plans. Thereby, the CPPT plans approach the dosimetric quality of IMPT-only plans, for which the mean brain BED2and mean liver BED4are reduced by 28% and 58%, respectively, compared to IMRT-only plans.Significance.CPPT with optimized proton and photon dose contributions to individual metastases may benefit selected metastatic cancer patients without tying up major proton resources.

A novel fast robust optimization algorithm for intensity-modulated proton therapy with minimum monitor unit constraint

BACKGROUND

Intensity-modulated proton therapy (IMPT) optimizes spot intensities and position, providing better conformability. However, the successful application of IMPT is dependent upon addressing the challenges posed by range and setup uncertainties. In order to address the uncertainties in IMPT, robust optimization is essential.

PURPOSE

This study aims to develop a novel fast algorithm for robust optimization of IMPT with minimum monitor unit (MU) constraint.

METHODS AND MATERIALS

The study formulates a robust optimization problem and proposes a novel, fast algorithm based on the alternating direction method of multipliers (ADMM) framework. This algorithm enables distributed computation and parallel processing. Ten clinical cases were used as test scenarios to evaluate the performance of the proposed approach. The robust optimization method (RBO-NEW) was compared with plans that only consider nominal optimization using CTV (NMO-CTV) without handling uncertainties and PTV (NMO-PTV) to handle the uncertainties, as well as with conventional robust-optimized plans (RBO-CONV). Dosimetric metrics, including D95, homogeneity index, and Dmean, were used to evaluate the dose distribution quality. The area under the root-mean-square dose (RMSD)-volume histogram curves (AUC) and dose-volume histogram (DVH) bands were used to evaluate the robustness of the treatment plan. Optimization time cost was also assessed to measure computational efficiency.

RESULTS

The results demonstrated that the RBO plans exhibited better plan quality and robustness than the NMO plans, with RBO-NEW showing superior computational efficiency and plan quality compared to RBO-CONV. Specifically, statistical analysis results indicated that RBO-NEW was able to reduce the computational time from389.70 ± 207.40 $389.70\pm 207.40$ to228.60 ± 123.67 $228.60\pm 123.67$ s (p < 0.01 $p<0.01$ ) and reduce the mean organ-at-risk (OAR) dose from9.38 ± 12.80 $9.38\pm 12.80$ % of the prescription dose to9.07 ± 12.39 $9.07\pm 12.39$ % of the prescription dose (p < 0.05 $p<0.05$ ) compared to RBO-CONV.

CONCLUSION

This study introduces a novel fast robust optimization algorithm for IMPT treatment planning with minimum MU constraint. Such an algorithm is not only able to enhance the plan's robustness and computational efficiency without compromising OAR sparing but also able to improve treatment plan quality and reliability.

Cardinality-constrained plan-quality and delivery-time optimization method for proton therapy

BACKGROUND

While minimizing plan delivery time is beneficial for proton therapy in terms of motion management, patient comfort, and treatment throughput, it often poses a tradeoff with optimizing plan quality. A key component of plan delivery time is the energy switching time, which is approximately proportional to the number of energy layers, that is, the cardinality.

PURPOSE

This work aims to develop a novel optimization method that can efficiently compute the pareto surface between plan quality and energy layer cardinality, for the planner to navigate through this quality-and-efficiency tradeoff and select the appropriate plan of a balanced tradeoff.

METHODS

A new IMPT method CARD is proposed that (1) explicitly incorporates the minimization of energy layer cardinality as an optimization objective, and (2) automatically generates a set of plans sequentially with a descending order in number of energy layers. The energy layer cardinality is penalized through the l1,0-norm regularization with an upper bound, and the upper bound is monotonically decreased to compute a series of treatment plans with gradually decreased energy layer cardinality on the quality-and-efficiency pareto surface. For any given treatment plan, the plan optimality is enforced using dose-volume planning objectives and the plan deliverability is imposed through minimum-monitor-unit (MMU) constraints, with optimization solution algorithm based on iterative convex relaxation.

RESULTS

The new method CARD was validated in comparison with the benchmark plan of all energy layers (P0), and a state-of-the-art method called MMSEL, using prostate, head-and-neck (HN), lung, pancreas, liver and brain cases. While labor-intensive and time-consuming manual parameter tuning was needed for MMSEL to generate plans of predefined energy layer cardinality, CARD automatically and efficiently computed all plans with sequentially decreasing predefined energy layer cardinality all at once. With the acceptable plan quality (i.e., no more than 110% of total optimization objective value from P0), CARD achieved the reduction of number of energy layers to 52% (from 77 to 40), 48% (from 135 to 65), 59% (from 85 to 50), 67% (from 52 to 35), 80% (from 50 to 40), and 30% (from 66 to 20), for prostate, HN, lung, pancreas, liver, and brain cases, respectively, compared to P0, with overall better plan quality than MMSEL. Moreover, due to the nonconvexity of the MMU constraint, CARD provided the similar or even smaller optimization objective than P0, at the same time with fewer number of energy layers, that is, 55 versus 77, 85 versus 135, 45 versus 52, and 25 versus 66 for prostate, HN, pancreas, and brain cases, respectively.

CONCLUSIONS

We have developed a novel optimization algorithm CARD that can efficiently and automatically compute a series of treatment plans of any given energy layer sequentially, which allows the planner to navigate through the plan-quality and energy-layer-cardinality tradeoff and select the appropriate plan of a balanced tradeoff.

Biological optimization for hybrid proton-photon radiotherapy

Objective.Hybrid proton-photon radiotherapy (RT) is a cancer treatment option to broaden access to proton RT. Additionally, with a refined treatment planning method, hybrid RT has the potential to offer superior plan quality compared to proton-only or photon-only RT, particularly in terms of target coverage and sparing organs-at-risk (OARs), when considering robustness to setup and range uncertainties. However, there is a concern regarding the underestimation of the biological effect of protons on OARs, especially those in close proximity to targets. This study seeks to develop a hybrid treatment planning method with biological dose optimization, suitable for clinical implementation on existing proton and photon machines, with each photon or proton treatment fraction delivering a uniform target dose.Approach.The proposed hybrid biological dose optimization method optimized proton and photon plan variables, along with the number of fractions for each modality, minimizing biological dose to the OARs and surrounding normal tissues. To mitigate underestimation of hot biological dose spots, proton biological dose was minimized within a ring structure surrounding the target. Hybrid plans were designed to be deliverable separately and robustly on existing proton and photon machines, with enforced uniform target dose constraints for the proton and photon fraction doses. A probabilistic formulation was utilized for robust optimization of setup and range uncertainties for protons and photons. The nonconvex optimization problem, arising from minimum monitor unit constraint and dose-volume histogram constraints, was solved using an iterative convex relaxation method.Main results.Hybrid planning with biological dose optimization effectively eliminated hot spots of biological dose, particularly in normal tissues surrounding the target, outperforming proton-only planning. It also provided superior overall plan quality and OAR sparing compared to proton-only or photon-only planning strategies.Significance.This study presents a novel hybrid biological treatment planning method capable of generating plans with reduced biological hot spots, superior plan quality to proton-only or photon-only plans, and clinical deliverability on existing proton and photon machines, separately and robustly.

Treatment planning comparison for head and neck cancer between photon, proton, and combined proton-photon therapy - From a fixed beam line to an arc

BACKGROUND AND PURPOSE

This study investigates whether combined proton-photon therapy (CPPT) improves treatment plan quality compared to single-modality intensity-modulated radiation therapy (IMRT) or intensity-modulated proton therapy (IMPT) for head and neck cancer (HNC) patients. Different proton beam arrangements for CPPT and IMPT are compared, which could be of specific interest concerning potential future upright-positioned treatments. Furthermore, it is evaluated if CPPT benefits remain under inter-fractional anatomical changes for HNC treatments.

MATERIAL AND METHODS

Five HNC patients with a planning CT and multiple (4-7) repeated CTs were studied. CPPT with simultaneously optimized photon and proton fluence, single-modality IMPT, and IMRT treatment plans were optimized on the planning CT and then recalculated and reoptimized on each repeated CT. For CPPT and IMPT, plans with different degrees of freedom for the proton beams were optimized. Fixed horizontal proton beam line (FHB), gantry-like, and arc-like plans were compared.

RESULTS

The target coverage for CPPT without adaptation is insufficient (average V95%=88.4 %), while adapted plans can recover the initial treatment plan quality for target (average V95%=95.5 %) and organs-at-risk. CPPT with increased proton beam flexibility increases plan quality and reduces normal tissue complication probability of Xerostomia and Dysphagia. On average, Xerostomia NTCP reductions compared to IMRT are -2.7 %/-3.4 %/-5.0 % for CPPT FHB/CPPT Gantry/CPPT Arc. The differences for IMPT FHB/IMPT Gantry/IMPT Arc are + 0.8 %/-0.9 %/-4.3 %.

CONCLUSION

CPPT for HNC needs adaptive treatments. Increasing proton beam flexibility in CPPT, either by using a gantry or an upright-positioned patient, improves treatment plan quality. However, the photon component is substantially reduced, therefore, the balance between improved plan quality and costs must be further determined.

Combined BNCT-CIRT treatment planning for glioblastoma using the effect-based optimization

Objective. Boron neutron capture therapy (BNCT) and carbon ion radiotherapy (CIRT) are emerging treatment modalities for glioblastoma. In this study, we investigated the methodology and feasibility to combine BNCT and CIRT treatments. The combined treatment plan illustrated how the synergistic utilization of BNCT's biological targeting and CIRT's intensity modulation capabilities could lead to optimized treatment outcomes.Approach. The Monte Carlo toolkit, TOPAS, was employed to calculate the dose distribution for BNCT, while matRad was utilized for the optimization of CIRT. The biological effect-based approach, instead of the dose-based approach, was adopted to develop the combined BNCT-CIRT treatment plans for six patients diagnosed with glioblastoma, considering the different radiosensitivity and fraction. Five optional combined treatment plans with specific BNCT effect proportions for each patient were evaluated to identify the optimal treatment that minimizes damage on normal tissue.Main results. Individual BNCT exhibits a significant effect gradient along with the beam direction in the large tumor, while combined BNCT-CIRT treatments can achieve uniform effect delivery within the clinical target volume (CTV) through the effect filling with reversed gradient by the CIRT part. In addition, the increasing BNCT effect proportion in combined treatments can reduce damage in the normal brain tissue near the CTV. Besides, the combined treatments effectively minimize damage to the skin compared to individual BNCT treatments.Significance. The initial endeavor to combine BNCT and CIRT treatment plans is achieved by the effect-based optimization. The observed advantages of the combined treatment suggest its potential applicability for tumors characterized by pleomorphic, infiltrative, radioresistant and voluminous features.

Lattice position optimization for LATTICE therapy

BACKGROUND

LATTICE radiation therapy delivers 3D heterogenous dose of high peak-to-valley dose ratio (PVDR) to the tumor target, with peak dose at lattice vertices inside the target and valley dose for the rest of the target. Although the lattice vertex positions can impact PVDR inside the target and sparing of organs-at-risk (OAR), they are fixed as constants and not optimized during treatment planning in current clinical practice.

PURPOSE

This work proposes a new LATTICE plan optimization method that can optimize lattice vertex positions during LATTICE treatment planning, which is the first lattice position optimization study to the best of our knowledge.

METHODS

The new LATTICE treatment planning method optimizes lattice vertex positions as well as other plan variables (e.g., photon fluences or proton spot weights), with optimization objectives for target PVDR and OAR sparing. To satisfy mathematical differentiability, the lattice vertices are approximated in sigmoid functions. For geometric feasibility, proper geometry constraints are enforced onto lattice vertex positions. The lattice position optimization problem is solved by iterative convex relaxation (ICR) method and alternating direction method of multipliers (ADMM), and lattice vertex positions and photon/proton plan variables are jointly updated via the Quasi-Newton method.

RESULTS

Both photon and proton LATTICE RT were considered, and the optimal lattice vertex positions in terms of plan objectives were found by solving all possible combinations on given discrete positions via exhaustive searching based on standard IMRT/IMPT, which served as the ground truth for validating the new LATTICE method. The results show that the new method indeed provided the optimal lattice vertex positions with the smallest optimization objective, the largest target PVDR, and the best OAR sparing.

CONCLUSIONS

A new LATTICE treatment planning method is proposed and validated that can optimize lattice vertex positions as well as other photon or proton plan variables for improving target PVDR and OAR sparing.

An orthogonal matching pursuit optimization method for solving minimum-monitor-unit problems: Applications to proton IMPT, ARC and FLASH

BACKGROUND

The intensities (i.e., number of protons in monitor unit [MU]) of deliverable proton spots need to be either zero or meet a minimum-MU (MMU) threshold, which is a nonconvex problem. Since the dose rate is proportionally associated with the MMU threshold, higher-dose-rate proton radiation therapy (RT) (e.g., efficient intensity modulated proton therapy (IMPT) and ARC proton therapy, and high-dose-rate-induced FLASH effect needs to solve the MMU problem with larger MMU threshold, which however makes the nonconvex problem more difficult to solve.

PURPOSE

This work will develop a more effective optimization method based on orthogonal matching pursuit (OMP) for solving the MMU problem with large MMU thresholds, compared to state-of-the-art methods, such as alternating direction method of multipliers (ADMM), proximal gradient descent method (PGD), or stochastic coordinate descent method (SCD).

METHODS

The new method consists of two essential components. First, the iterative convex relaxation (ICR) method is used to determine the active sets for dose-volume planning constraints and decouple the MMU constraint from the rest. Second, a modified OMP optimization algorithm is used to handle the MMU constraint: the non-zero spots are greedily selected via OMP to form the solution set to be optimized, and then a convex constrained subproblem is formed and can be conveniently solved to optimize the spot weights restricted to this solution set via OMP. During this iterative process, the new non-zero spots localized via OMP will be adaptively added to or removed from the optimization objective.

RESULTS

The new method via OMP is validated in comparison with ADMM, PGD and SCD for high-dose-rate IMPT, ARC, and FLASH problems of large MMU thresholds, and the results suggest that OMP substantially improved the plan quality from PGD, ADMM and SCD in terms of both target dose conformality (e.g., quantified by max target dose and conformity index) and normal tissue sparing (e.g., mean and max dose). For example, in the brain case, the max target dose for IMPT/ARC/FLASH was 368.0%/358.3%/283.4% respectively for PGD, 154.4%/179.8%/150.0% for ADMM, 134.5%/130.4%/123.0% for SCD, while it was <120% in all scenarios for OMP; compared to PGD/ADMM/SCD, OMP improved the conformity index from 0.42/0.52/0.33 to 0.65 for IMPT and 0.46/0.60/0.61 to 0.83 for ARC.

CONCLUSIONS

A new OMP-based optimization algorithm is developed to solve the MMU problems with large MMU thresholds, and validated using examples of IMPT, ARC, and FLASH with substantially improved plan quality from ADMM, PGD, and SCD.

Beam angle optimization for proton therapy via group-sparsity based angle generation method

BACKGROUND

In treatment planning, beam angle optimization (BAO) refers to the selection of a subset with a given number of beam angles from all available angles that provides the best plan quality. BAO is a NP-hard combinatorial problem. Although exhaustive search (ES) can exactly solve BAO by exploring all possible combinations, ES is very time-consuming and practically infeasible.

PURPOSE

To the best of our knowledge, (1) no optimization method has been demonstrated that can provide the exact solution to BAO, and (2) no study has validated an optimization method for solving BAO by benchmarking with the optimal BAO solution (e.g., via ES), both of which will be addressed by this work.

METHODS

This work considers BAO for proton therapy, for example, the selection of 2-4 beam angles for IMPT. The optimal BAO solution is obtained via ES and serves as the ground truth. A new BAO algorithm, namely angle generation (AG) method, is proposed, and demonstrated to provide nearly-exact solutions for BAO in reference to the ES solution. AG iteratively optimizes the angular set via group-sparsity (GS) regularization, until the planning objective does not decrease further.

RESULTS

Since GS alone can also solve BAO, AG was validated and compared with GS for 2-angle brain, 3-angle lung, and 4-angle brain cases, in reference to the optimal BAO solutions obtained by ES: the AG solution had the rank (1/276, 1/2024, 4/10 626), while the GS solution had the rank (42/276, 279/2024, 4328/10 626).

CONCLUSIONS

A new BAO algorithm called AG is proposed and shown to provide substantially improved accuracy for BAO from current methods with nearly-exact solutions to BAO, in reference to the ground truth of optimal BAO solution via ES.

Global democratisation of proton radiotherapy

Proton radiotherapy is an advanced treatment option compared with conventional x-ray treatment, delivering much lower doses of radiation to healthy tissues surrounding the tumour. However, proton therapy is currently not widely available. In this Review, we summarise the evolution of proton therapy to date, together with the benefits to patients and society. These developments have led to an exponential growth in the number of hospitals using proton radiotherapy worldwide. However, the gap between the number of patients who should be treated with proton radiotherapy and those who have access to it remains large. We summarise the ongoing research and development that is contributing to closing this gap, including the improvement of treatment efficiency and efficacy, and advances in fixed-beam treatments that do not require an enormously large, heavy, and costly gantry. The ultimate goal of decreasing the size of proton therapy machines to fit into standard treatment rooms appears to be within reach, and we discuss future research and development opportunities to achieve this goal.

A treatment plan optimization method with direct minimization of number of energy jumps for proton arc therapy

Objective. The optimization of energy layer distributions is crucial to proton arc therapy: on one hand, a sufficient number of energy layers is needed to ensure the plan quality; on the other hand, an excess number of energy jumps (NEJ) can substantially slow down the treatment delivery. This work will develop a new treatment plan optimization method with direct minimization of (NEJ), which will be shown to outperform state-of-the-art methods in both plan quality and delivery efficiency.Approach. The proposed method jointly optimizes the plan quality and minimizes the NEJ. To minimize NEJ, (1) the proton spotsxis summed per energy layer to form the energy vectory; (2)yis binarized via sigmoid transform intoy1; (3)y1is multiplied with a predefined energy order vector via dot product intoy2; (4)y2is filtered through the finite-differencing kernel intoy3in order to identify NEJ; (5) only the NEJ ofy3is penalized, whilexis optimized for plan quality. The solution algorithm to this new method is based on iterative convex relaxation.Main results. The new method is validated in comparison with state-of-the-art methods called energy sequencing (ES) method and energy matrix (EM) method. In terms of delivery efficiency, the new method had fewer NEJ, less energy switching time, and generally less total delivery time. In terms of plan quality, the new method had smaller optimization objective values, lower normal tissue dose, and generally better target coverage.Significance. We have developed a new treatment plan optimization method with direct minimization of NEJ, and demonstrated that this new method outperformed state-of-the-art methods (ES and EM) in both plan quality and delivery efficiency.

TVL1-IMPT: Optimization of Peak-to-Valley Dose Ratio Via Joint Total-Variation and L1 Dose Regularization for Spatially Fractionated Pencil-Beam-Scanning Proton Therapy

PURPOSE

Proton minibeam radiation therapy (pMBRT) is a novel proton modality of spatially fractionated RT. pMBRT can reduce the radiation damage to normal tissues via biological dose sparing of high peak-to-valley dose ratio (PVDR). This work will develop a new spatially fractionated IMPT treatment planning method for pMBRT that jointly optimizes the plan quality and maximizes the PVDR.

METHODS

The new optimization method simultaneously maximizes the normal-tissue PVDR and optimizes the dose distribution at tumor targets and organs at risk. The PVDR maximization is through the joint total variation (TV) and L1 regularization with respect to the normal-tissue dose. That is, the beam-eye view projects dose slices of several depths for each beam angle; the TV of dose is maximized, corresponding to the PVDR maximization; and the L1 of dose is minimized, corresponding to the minimization of the organs-at-risk dose and maximization of survival fraction (SF).

RESULTS

The new IMPT method with TV and L1 regularization was validated in comparison with the conventional IMPT method for pMBRT in several clinical cases. The results show that TVL1 provided larger PVDR and SF than the conventional IMPT method for biological sparing of normal tissues, with preserved plan quality in terms of physical dose distribution.

CONCLUSIONS

A new spatially fractionated IMPT treatment planning method was developed for pMBRT that can optimize and improve normal-tissue PVDR and SF by incorporating TV and L1 dose regularization with properly chosen regularization parameters into IMPT.

An iterative convex relaxation method for proton LET optimization

Objective:A constant relative biological effectiveness of 1.1 in current clinical practice of proton radiotherapy (RT) is a crude approximation and may severely underestimate the biological dose from proton RT to normal tissues, especially near the treatment target at the end of Bragg peaks that exhibits high linear energy transfer (LET). LET optimization can account for biological effectiveness of protons during treatment planning, for minimizing biological proton dose and hot spots to normal tissues. However, the LET optimization is usually nonlinear and nonconvex to solve, for which this work will develop an effective optimization method based on iterative convex relaxation (ICR).Approach: In contrast to the generic nonlinear optimization method, such as Quasi-Newton (QN) method, that does not account for specific characteristics of LET optimization, ICR is tailored to LET modeling and optimization in order to effectively and efficiently solve the LET problem. Specifically, nonlinear dose-averaged LET term is iteratively linearized and becomes convex during ICR, while nonconvex dose-volume constraint and minimum-monitor-unit constraint are also handled by ICR, so that the solution for LET optimization is obtained by solving a sequence of convex and linearized convex subproblems. Since the high LET mostly occurs near the target, a 1 cm normal-tissue expansion of clinical target volume (CTV) (excluding CTV), i.e. CTV1cm, is defined to as an auxiliary structure during treatment planning, where LET is minimized.Main results: ICR was validated in comparison with QN for abdomen, lung, and head-and-neck cases. ICR was effective for LET optimization, as ICR substantially reduced the LET and biological dose in CTV1cm the ring, with preserved dose conformality to CTV. Compared to QN, ICR had smaller LET, physical and biological dose in CTV1cm, and higher conformity index values; ICR was also computationally more efficient, which was about 3 times faster than QN.Significance: A LET-specific optimization method based on ICR has been developed for solving proton LET optimization, which has been shown to be more computationally efficient than generic nonlinear optimizer via QN, with better plan quality in terms of LET, biological and physical dose conformality.

Fraction optimization for hybrid proton-photon treatment planning

BACKGROUND

Hybrid proton-photon radiotherapy (RT) can provide better plan quality than proton or photon only RT, in terms of robustness of target coverage and sparing of organs-at-risk (OAR).

PURPOSE

This work develops a hybrid treatment planning method that can optimize the number of proton and photon fractions as well as proton and photon plan variables, so that the hybrid plans can be clinically delivered day-to-day using either proton or photon machine.

METHODS

In the new hybrid treatment planning method, the total dose distribution (sum of proton dose and photon dose) is optimized for robust target coverage and optimal OAR sparing, by jointly optimizing proton spots and photon fluences, while the target dose uniformity is also enforced individually on both proton dose and photon dose, so that the hybrid plans can be separately and robustly delivered on proton or photon machine. To ensure the target dose uniformity for proton and photon plans, the number of proton and photon fractions is explicitly modeled and optimized, so that the target dose for proton and photon dose components is constrained to be a constant fraction of the total prescription dose while the plan quality based on total dose is optimized. The feasibility of hybrid planning using the proposed method is validated with representative clinical cases including abdomen, lung, head-and-neck (HN), and brain.

RESULTS

For all cases, hybrid plans provided better overall plan quality and OAR sparing than proton-only or photon-only plans, better target dose uniformity and robustness than proton-only plans, quantified by treatment planning objectives and dosimetric parameters. Moreover, for HN and brain cases, hybrid plans also had better target coverage than photon-only plans.

CONCLUSIONS

We have developed a new hybrid treatment planning method that optimizes number of proton and photon fractions as well as proton spots and photon fluences, for generating hybrid plans that can be separately and robustly delivered on proton or photon machines. Preliminary results have demonstrated that hybrid plans generated by the new method have better plan quality than proton-only or photon-only plans.

Efficiency enhancements of a Monte Carlo beamlet based treatment planning process: implementation and parameter study

Objective.The computational effort to perform beamlet calculation, plan optimization and final dose calculation of a treatment planning process (TPP) generating intensity modulated treatment plans is enormous, especially if Monte Carlo (MC) simulations are used for dose calculation. The goal of this work is to improve the computational efficiency of a fully MC based TPP for static and dynamic photon, electron and mixed photon-electron treatment techniques by implementing multiple methods and studying the influence of their parameters.Approach.A framework is implemented calculating MC beamlets efficiently in parallel on each available CPU core. The user can specify the desired statistical uncertainty of the beamlets, a fractional sparse dose threshold to save beamlets in a sparse format and minimal distances to the PTV surface from which 2 × 2 × 2 = 8 (medium) or even 4 × 4 × 4 = 64 (large) voxels are merged. The compromise between final plan quality and computational efficiency of beamlet calculation and optimization is studied for several parameter values to find a reasonable trade-off. For this purpose, four clinical and one academic case are considered with different treatment techniques.Main results.Setting the statistical uncertainty to 5% (photon beamlets) and 15% (electron beamlets), the fractional sparse dose threshold relative to the maximal beamlet dose to 0.1% and minimal distances for medium and large voxels to the PTV to 1 cm and 2 cm, respectively, does not lead to substantial degradation in final plan quality compared to using 2.5% (photon beamlets) and 5% (electron beamlets) statistical uncertainty and no sparse format nor voxel merging. Only OAR sparing is slightly degraded. Furthermore, computation times are reduced by about 58% (photon beamlets), 88% (electron beamlets) and 96% (optimization).Significance.Several methods are implemented improving computational efficiency of beamlet calculation and plan optimization of a fully MC based TPP without substantial degradation in final plan quality.

A novel stochastic optimization method for handling misalignments of proton and photon doses in combined treatments

Objective. Combined proton–photon treatments, where most fractions are delivered with photons and only a few are delivered with protons, may represent a practical approach to optimally use limited proton resources. It has been shown that, when organs at risk (OARs) are located within or near the tumor, the optimal multi-modality treatment uses protons to hypofractionate parts of the target volume and photons to achieve near-uniform fractionation in dose-limiting healthy tissues, thus exploiting the fractionation effect. These plans may be sensitive to range and setup errors, especially misalignments between proton and photon doses. Thus, we developed a novel stochastic optimization method to directly incorporate these uncertainties into the biologically effective dose (BED)-based simultaneous optimization of proton and photon plans. Approach. The method considers the expected value E b and standard deviation σ b of the cumulative BED b in every voxel of a structure. For the target, a piecewise quadratic penalty function of the form b min − E b − 2 σ b + 2 is minimized, aiming for plans in which the expected BED minus two times the standard deviation exceeds the prescribed BED b min . Analogously, E b + 2 σ b − b max + 2 is considered for OARs. Main results. Using a spinal metastasis case and a liver cancer patient, it is demonstrated that the novel stochastic optimization method yields robust combined treatment plans. Tumor coverage and a good sparing of the main OARs are maintained despite range and setup errors, and especially misalignments between proton and photon doses. This is achieved without explicitly considering all combinations of proton and photon error scenarios. Significance. Concerns about range and setup errors for safe clinical implementation of optimized proton–photon radiotherapy can be addressed through an appropriate stochastic planning method.

Energy layer optimization via energy matrix regularization for proton spot-scanning arc therapy

PURPOSE

Spot-scanning arc therapy (SPArc) is an emerging proton modality that can potentially offer a combination of advantages in plan quality and delivery efficiency, compared to traditional IMPT of a few beam angles. Unlike IMPT, frequent low-to-high energy layer switching (so called switch-up (SU)) can degrade delivery efficiency for SPArc. However, it is a tradeoff between the minimization of SU times and the optimization of plan quality. This work will consider the energy layer optimization (ELO) problem for SPArc and develop a new ELO method via energy matrix (EM) regularization to improve plan quality and delivery efficiency.

METHODS

The major innovation of EM method for ELO is to design an energy matrix that encourages desirable energy-layer map with minimal SU during SPArc, and then incorporate this energy matrix into the SPArc treatment planning to simultaneously minimize the number of SU and optimize plan quality. The EM method is solved by the fast iterative shrinkage-thresholding algorithm and validated in comparison with a state-of-the-art method, so-called energy sequencing (ES).

RESULTS

EM is validated and compared with ES using representative clinical cases. In terms of delivery efficiency, EM had fewer SU than ES with an average of 35% reduction of SU. In terms of plan quality, compared to ES, EM had smaller optimization objective values and better target dose conformality, and generally lower dose to organs-at-risk and lower integral dose to body. In terms of computational efficiency, EM was substantially more efficient than ES by at least ten-fold.

CONCLUSION

We have developed a new ELO method for SPArc using energy matrix regularization, and shown that this new method EM can improve both delivery efficiency and plan quality, with substantially reduced computational time, compared to ES. This article is protected by copyright. All rights reserved.

Simultaneous dose and dose rate optimization (SDDRO) of the FLASH effect for pencil-beam-scanning proton therapy

PURPOSE

Compared to CONV-RT (with conventional dose rate), FLASH-RT (with ultra-high dose rate) can provide biological dose sparing for organs-at-risk (OARs) via the so-called FLASH effect, in addition to physical dose sparing. However, the FLASH effect only occurs, when both dose and dose rate meet certain minimum thresholds. This work will develop a simultaneous dose and dose rate optimization (SDDRO) method accounting for both FLASH dose and dose rate constraints during treatment planning for pencil-beam-scanning proton therapy.

METHODS

SDDRO optimizes the FLASH effect (specific to FLASH-RT) as well as the dose distribution (similar to CONV-RT). The nonlinear dose rate constraint is linearized, and the reformulated optimization problem is efficiently solved via iterative convex relaxation powered by alternating direction method of multipliers. To resolve and quantify the generic tradeoff of FLASH-RT between FLASH and dose optimization, we propose the use of FLASH effective dose based on dose modifying factor (DMF) owing to the FLASH effect.

RESULTS

FLASH-RT via transmission beams (TB) (IMPT-TB or SDDRO) and CONV-RT via Bragg peaks (BP) (IMPT-BP) were evaluated for clinical prostate, lung, head-and-neck (HN), and brain cases. Despite the use of TB, which is generally suboptimal to BP for normal tissue sparing, FLASH-RT via SDDRO considerably reduced FLASH effective dose for high-dose OAR adjacent to the target. For example, in the lung SBRT case, the max esophageal dose constraint 27 Gy was only met by SDDRO (24.8 Gy), compared to IMPT-BP (35.3 Gy) or IMPT-TB (36.6 Gy); in the brain SRS case, the brain constraint V12Gy≤15cc was also only met by SDDRO (13.7cc), compared to IMPT-BP (43.9cc) or IMPT-TB (18.4cc). In addition, SDDRO substantially improved the FLASH coverage from IMPT-TB, e.g., an increase from 37.2% to 67.1% for lung, from 39.1% to 58.3% for prostate, from 65.4% to 82.1% for HN, from 50.8% to 73.3% for the brain.

CONCLUSIONS

Both FLASH dose and dose rate constraints are incorporated into SDDRO for FLASH-RT that jointly optimizes the FLASH effect and physical dose distribution. FLASH effective dose via FLASH DMF is introduced to reconcile the tradeoff between physical dose sparing and FLASH sparing, and quantify the net effective gain from CONV-RT to FLASH-RT.

Minimum-monitor-unit optimization via a stochastic coordinate descent method

Objective. Deliverable proton spots are subject to the minimum monitor-unit (MMU) constraint. The MMU optimization problem with relatively large MMU threshold remains mathematically challenging due to its strong nonconvexity. However, the MMU optimization is fundamental to proton radiotherapy (RT), including efficient IMPT and proton arc delivery (ARC). This work aims to develop a new optimization algorithm that is effective in solving the MMU problem.Approach.Our new algorithm is primarily based on stochastic coordinate decent (SCD) method. It involves three major steps: first to decouple the determination of active sets for dose-volume-histogram (DVH) planning constraints from the MMU problem via iterative convex relaxation method; second to handle the nonconvexity of the MMU constraint via SCD to localize the index set of nonzero spots; third to solve convex subproblems projected to this convex set of nonzero spots via projected gradient descent method.Main results.Our new method SCD is validated and compared with alternating direction method of multipliers (ADMM) for IMPT and ARC. The results suggest SCD had better plan quality than ADMM, e.g. the improvement of conformal index (CI) from 0.56 to 0.69 during IMPT, and from 0.28 to 0.80 during ARC for the lung case. Moreover, SCD successfully handled the nonconvexity from large MMU threshold that ADMM failed to handle, in the sense that (1) the plan quality from ARC was worse than IMPT (e.g. CI was 0.28 with IMPT and 0.56 with ARC for the lung case), when ADMM was used; (2) in contrast, with SCD, ARC achieved better plan quality than IMPT (e.g. CI was 0.69 with IMPT and 0.80 with ARC for the lung case), which is compatible with more optimization degrees of freedom from ARC compared to IMPT.Significance. To the best of our knowledge, our new MMU optimization method via SCD can effectively handle the nonconvexity from large MMU threshold that none of the current methods can solve. Therefore, we have developed a unique MMU optimization algorithm via SCD that can be used for efficient IMPT, proton ARC, and other particle RT applications where large MMU threshold is desirable (e.g. for the delivery of high dose rates or/and a large number of spots).

An adaptive spot placement method on Cartesian grid for pencil beam scanning proton therapy

Pencil beam scanning proton radiotherapy (RT) offers flexible proton spot placement near treatment targets for delivering tumoricidal radiation dose to tumor targets while sparing organs-at-risk. Currently the spot placement is mostly based on a non-adaptive sampling (NS) strategy on a Cartesian grid. However, the spot density or spacing during NS is a constant for the Cartesian grid that is independent of the geometry of tumor targets, and thus can be suboptimal in terms of plan quality (e.g. target dose conformality) and delivery efficiency (e.g. number of spots). This work develops an adaptive sampling (AS) spot placement method on the Cartesian grid that fully accounts for the geometry of tumor targets. Compared with NS, AS places (1) a relatively fine grid of spots at the boundary of tumor targets to account for the geometry of tumor targets and treatment uncertainties (setup and range uncertainty) for improving dose conformality, and (2) a relatively coarse grid of spots in the interior of tumor targets to reduce the number of spots for improving delivery efficiency and robustness to the minimum-minitor-unit (MMU) constraint. The results demonstrate that (1) AS achieved comparable plan quality with NS for regular MMU and substantially improved plan quality from NS for large MMU, using merely about 10% of spots from NS, where AS was derived from the same Cartesian grid as NS; (2) on the other hand, with similar number of spots, AS had better plan quality than NS consistently for regular and large MMU.

Optimal Allocation of Proton Therapy Slots in Combined Proton-Photon Radiation Therapy

PURPOSE

Proton therapy is a limited resource that is not available to all patients who may benefit from it. We investigated combined proton-photon treatments, in which some fractions are delivered with protons and the remaining fractions with photons, as an approach to maximize the benefit of limited proton therapy resources at a population level.

METHODS AND MATERIALS

To quantify differences in normal-tissue complication probability (NTCP) between protons and photons, we considered a cohort of 45 patients with head and neck cancer for whom intensity modulated radiation therapy and intensity modulated proton therapy plans were previously created, in combination with NTCP models for xerostomia and dysphagia considered in the Netherlands for proton patient selection. Assuming limited availability of proton slots, we developed methods to optimally assign proton fractions in combined proton-photon treatments to minimize the average NTCP on a population level. The combined treatments were compared with patient selection strategies in which patients are assigned to single-modality proton or photon treatments.

RESULTS

There is a benefit of combined proton-photon treatments compared with patient selection, owing to the nonlinearity of NTCP functions; that is, the initial proton fractions are the most beneficial, whereas additional proton fractions have a decreasing benefit when a flatter part of the NTCP curve is reached. This effect was small for the patient cohort and NTCP models considered, but it may be larger if dose-response relationships are better known. In addition, when proton slots are limited, patient selection methods face a trade-off between leaving slots unused and blocking slots for future patients who may have a larger benefit. Combined proton-photon treatments with flexible proton slot assignment provide a method to make optimal use of all available resources.

CONCLUSIONS

Combined proton-photon treatments allow for better use of limited proton therapy resources. The benefit over patient selection schemes depends on the NTCP models and the dose differences between protons and photons.

SDDRO-joint: simultaneous dose and dose rate optimization with the joint use of transmission beams and Bragg peaks for FLASH proton therapy

Cancer radiotherapy (RT) with the irradiation at ultra-high dose rates, namely FLASH-RT, can substantially reduce radiation-induced normal tissue toxicities while maintaining tumor response. Currently, clinical FLASH-RT on deep-seated tumors can only be performed with proton beams. One way to achieve ultra-high dose rates at depth is through the use of high-energy transmission beams (TB), where the Bragg peaks (BP) fall outside the body. However, planning with TB alone does not fully leverage the degrees of freedom for dose shaping as traditional intensity modulated proton therapy (IMPT) which uses the BP of multi-energy proton beams at the tumor target. This work will develop a simultaneous dose and dose rate optimization (SDDRO) method with the joint use of TB and BP, namely SDDRO-Joint. Specifically, BP are placed inside tumor targets to improve the target dose conformality and sparse the normal-tissue dose, while TB primarily cover the tumor boundary to achieve ultra-high dose rate coverage of organs-at-risk (OAR) close to tumor targets. The sparing of OAR and other normal tissues via SDDRO-Joint is jointly by TB and BP, i.e. the FLASH sparing by TB and the dose sparing by BP. The results suggest that the addition of BP substantially increased the target dose conformality for SDDRO. Noticeably SDDRO-Joint also provided slightly higher conformal index values than the conventional IMPT method with BP alone.

Joint Optimization of Photon-Carbon Ion Treatments for Glioblastoma

PURPOSE

Carbon ions are radiobiologically more effective than photons and are beneficial for treating radioresistant gross tumor volumes (GTV). However, owing to a reduced fractionation effect, they may be disadvantageous for treating infiltrative tumors, in which healthy tissue inside the clinical target volume (CTV) must be protected through fractionation. This work addresses the question: What is the ideal combined photon-carbon ion fluence distribution for treating infiltrative tumors given a specific fraction allocation between photons and carbon ions?

METHODS AND MATERIALS

We present a method to simultaneously optimize sequentially delivered intensity modulated photon (IMRT) and carbon ion (CIRT) treatments based on cumulative biological effect, incorporating both the variable relative biological effect of carbon ions and the fractionation effect within the linear quadratic model. The method is demonstrated for 6 glioblastoma patients in comparison with the current clinical standard of independently optimized CIRT-IMRT plans.

RESULTS

Compared with the reference plan, joint optimization strategies yield inhomogeneous photon and carbon ion dose distributions that cumulatively deliver a homogeneous biological effect distribution. In the optimal distributions, the dose to CTV is mostly delivered by photons and carbon ions are restricted to the GTV with variations depending on tumor size and location. Improvements in conformity of high-dose regions are reflected by a mean EQD2 reduction of 3.29 ± 1.22 Gy in a dose fall-off margin around the CTV. Carbon ions may deliver higher doses to the center of the GTV, and photon contributions are increased at interfaces with CTV and critical structures. This results in a mean EQD2 reduction of 8.3 ± 2.28 Gy, in which the brain stem abuts the target volumes.

CONCLUSIONS

We have developed a biophysical model to optimize combined photon-carbon ion treatments. For 6 glioblastoma patient cases, we show that our approach results in a more targeted application of carbon ions that (1) reduces dose in normal tissues within the target volume, which can only be protected through fractionation; and (2) boosts central target volume regions to reduce integral dose. Joint optimization of IMRT-CIRT treatments enable the exploration of a new spectrum of plans that can better address physical and radiobiological treatment planning challenges.

Simultaneous dose and dose rate optimization (SDDRO) for FLASH proton therapy

PURPOSE

FLASH radiotherapy (RT) can potentially reduce normal tissue toxicity while preserving tumoricidal effectiveness to improve the therapeutic ratio. The key of FLASH for sparing normal tissues is to irradiate tissues with an ultra-high dose rate (i.e., ≥40 Gy/s), for which proton RT can be used. However, currently available treatment plan optimization method only optimizes the dose distribution and does not directly optimize the dose rate. The contribution of this work to FLASH proton RT is the development of a novel treatment optimization method, that is, simultaneous dose and dose rate optimization (SDDRO), to optimize tissue-receiving dose rate distribution as well as dose distribution.

METHODS

Distinguished from existing methods, SDDRO accounts for dose rate constraint and optimizes dose rate distribution. In terms of mathematical formulation, SDDRO is a constrained optimization problem with dose-volume constraint on dose distribution, minimum dose rate constraint on dose-averaged tissue-receiving dose rates, minimum monitor unit constraint on spot weight, and maximum intensity constraint on beam intensity. In terms of optimization algorithm, SDDRO is solved by iterative convex relaxation and alternating direction method of multipliers. SDDRO algorithms are presented for both scenarios with either constant or variable beam intensity.

RESULTS

SDDRO was compared with intensity modulated proton therapy (IMPT) (dose optimization alone, and no dose rate optimization) using three lung cases. SDDRO substantially improved the dose rate distribution compared to IMPT, for example, increasing of the region-of-interest (ROI) volume (ROI = CTV_10mm: the ring sandwiched by 10 mm outer and inner expansion of CTV boundary) receiving at least 40 Gy/s from ~30-50% to at least 98%, and the lung volume receiving at least 40 Gy/s from ~30-40% to ~70-90%. Moreover, both dose and dose rate distributions from SDDRO were further considerably improved via the combined use of hypofractionation and multiple beams.

CONCLUSIONS

We have developed a joint dose and dose rate optimization method for FLASH proton RT, namely SDDRO, which is first-of-its-kind to the best of our knowledge. The results suggest that (a) SDDRO can substantially improve the FLASH-dose rate coverage (e.g., in terms of dose rate volume histogram) compared to IMPT for the purpose of normal tissue sparing while preserving the dose distribution and (b) the combination of hypofractionation and multiple beams can further considerably improve the SDDRO plan quality in terms of both dose and dose rate distribution.

A treatment planning study of combined carbon ion-beam plus photon intensity-modulated radiotherapy

BACKGROUND AND PURPOSE

Combined photon intensity-modulated radiotherapy (IMRT) and sequential dose-escalated carbon ion beam therapy (IBT) is a technically advanced treatment option for head and neck malignancies. We proposed and evaluated an integrated planning strategy as opposed to an established and largely separated planning workflow.

MATERIALS AND METHODS

Ten patients with representative malignancies of the head and neck region underwent combined carbon-photon radiotherapy (RT) in our facilities. Clinical plans were created according to the separated workflow with independent optimization stages for both modalities. Experimental plans incorporated the existing carbon IBT dose distribution into the optimization stage of a step-and-shoot photon IMRT (bias dose planning).

RESULTS

Cumulative dose distributions showed statistically significant differences between the two planning strategies and were predominantly in favor of the integrated approach. As such, target irradiation was generally maintained or even improved in a subset of metrics, while normal tissue sparing was widely enhanced; for instance, in the ipsilateral temporal lobe with median Dmean of -16% (p < 0.001). Maximum doses D1% (with adjustment for different fractionation) fell below thresholds for toxicity risk in a minority of instances, where they were previously exceeded. Integral dose did not differ significantly.

CONCLUSIONS

Our findings indicate that combination planning of carbon-photon RT for head and neck malignancies may benefit from a proposed bias dose method, yielding favorable dose distribution characteristics and a streamlined planning workflow with fewer plan revisions. Further research is necessary to validate these observations in terms of robustness and their potential for higher tumor control.

Accounting for Range Uncertainties in the Optimization of Combined Proton-Photon Treatments Via Stochastic Optimization

PURPOSE

Proton treatment slots are a limited resource. Combined proton-photon treatments, in which most fractions are delivered with photons and only a few with protons, may represent a practical solution to optimize the allocation of proton resources over the patient population. We demonstrate how a limited number of proton fractions can be optimally used in multimodality treatments and address the issue of the robustness of combined treatments against proton range uncertainties.

METHODS AND MATERIALS

Combined proton-photon treatments are planned by simultaneously optimizing intensity modulated radiation therapy and proton therapy plans while accounting for the fractionation effect through the biologically effective dose model. The method was investigated for different tumor sites (a spinal metastasis, a sacral chordoma, and an atypical meningioma) in which organs at risk (OARs) were located within or near the tumor. Stochastic optimization was applied to mitigate range uncertainties.

RESULTS

In optimal combinations, proton and photon fractions deliver similar doses to OARs overlaying the target volume to protect these dose-limiting normal tissues through fractionation. Meanwhile, parts of the tumor are hypofractionated with protons. Thus, the total dose delivered with photons is reduced compared with simple combinations in which each modality delivers the prescribed dose per fraction to the target volume. The benefit of optimal combinations persists when range errors are accounted for via stochastic optimization.

CONCLUSIONS

Limited proton resources are optimally used in combined treatments if parts of the tumor are hypofractionated with protons and near-uniform fractionation is maintained in serial OARs. Proton range uncertainties can be efficiently accounted for through stochastic optimization and are not an obstacle for clinical application.