Reformulated McNamara RBE-weighted beam orientation optimization for intensity modulated proton therapy

Hypoxia-informed RBE-weighted beam orientation optimization for intensity modulated proton therapy

BACKGROUND

Variable relative biological effectiveness (RBE) models in treatment planning have been proposed to optimize the therapeutic ratio of proton therapy. It has been reported that proton RBE decreases with increasing tumor oxygen level, offering an opportunity to address hypoxia-related radioresistance with RBE-weighted optimization.

PURPOSE

Here, we obtain a voxel-level estimation of partial oxygen pressure to weigh RBE values in a single biologically informed beam orientation optimization (BOO) algorithm.

METHODS

Three glioblastoma patients with [18 F]-fluoromisonidazole (FMISO)-PET/CT images were selected from the institutional database. Oxygen values were derived from tracer uptake using a nonlinear least squares curve fitting. McNamara RBE, calculated from proton dose, was then weighed using oxygen enhancement ratios (OER) for each voxel and incorporated into the dose fidelity term of the BOO algorithm. The nonlinear optimization problem was solved using a split-Bregman approach, with FISTA as the solver. The proposed hypoxia informed RBE-weighted method (HypRBE) was compared to dose fidelity terms using the constant RBE of 1.1 (cRBE) and the normoxic McNamara RBE model (RegRBE). Tumor homogeneity index (HI), maximum biological dose (Dmax), and D95%, as well as OAR therapeutic index (TI = gEUDCTV /gEUDOAR ) were evaluated along with worst-case statistics after normalization to normal tissue isotoxicity.

RESULTS

Compared to [cRBE, RegRBE], HypRBE increased tumor HI, Dmax, and D95% across all plans by on average [31.3%, 31.8%], [48.6%, 27.1%], and [50.4%, 23.8%], respectively. In the worst-case scenario, the parameters increase on average by [12.5%, 14.7%], [7.3%,-8.9%], and [22.3%, 2.1%]. Despite increased OAR Dmean and Dmax by [8.0%, 3.0%] and [13.1%, -0.1%], HypRBE increased average TI by [22.0%, 21.1%]. Worst-case OAR Dmean, Dmax, and TI worsened by [17.9%, 4.3%], [24.5%, -1.2%], and [9.6%, 10.5%], but in the best cases, HypRBE escalates tumor coverage significantly without compromising OAR dose, increasing the therapeutic ratio.

CONCLUSIONS

We have developed an optimization algorithm whose dose fidelity term accounts for hypoxia-informed RBE values. We have shown that HypRBE selects bE:\Alok\aaeams better suited to deliver high physical dose to low RBE, hypoxic tumor regions while sparing the radiosensitive normal tissue.

Impact of proton therapy on the DNA damage induction and repair in hematopoietic stem and progenitor cells

Proton therapy is of great interest to pediatric cancer patients because of its optimal depth dose distribution. In view of healthy tissue damage and the increased risk of secondary cancers, we investigated DNA damage induction and repair of radiosensitive hematopoietic stem and progenitor cells (HSPCs) exposed to therapeutic proton and photon irradiation due to their role in radiation-induced leukemia. Human CD34+ HSPCs were exposed to 6 MV X-rays, mid- and distal spread-out Bragg peak (SOBP) protons at doses ranging from 0.5 to 2 Gy. Persistent chromosomal damage was assessed with the micronucleus assay, while DNA damage induction and repair were analyzed with the γ-H2AX foci assay. No differences were found in induction and disappearance of γ-H2AX foci between 6 MV X-rays, mid- and distal SOBP protons at 1 Gy. A significantly higher number of micronuclei was found for distal SOBP protons compared to 6 MV X-rays and mid- SOBP protons at 0.5 and 1 Gy, while no significant differences in micronuclei were found at 2 Gy. In HSPCs, mid-SOBP protons are as damaging as conventional X-rays. Distal SOBP protons showed a higher number of micronuclei in HSPCs depending on the radiation dose, indicating possible changes of the in vivo biological response.

Analytical HDR prostate brachytherapy planning with automatic catheter and isotope selection

BACKGROUND

High dose rate (HDR) brachytherapy is commonly used to treat prostate cancer. Existing HDR planning systems solve the dwell time problem for predetermined catheters and a single energy source.

PURPOSE

Additional degrees of freedom can be obtained by relaxing the catheters' pre-designation and introducing more source types, and may have a dosimetric benefit, particularly in improving conformality to spare the urethra. This study presents a novel analytical approach to solving the corresponding HDR planning problem.

METHODS

The catheter and dual-energy source selection problem was formulated as a constrained optimization problem with a non-convex group sparsity regularization. The optimization problem was solved using the fast-iterative shrinkage-thresholding algorithm (FISTA). Two isotopes were considered. The dose rates for the HDR 4140 Ytterbium (Yb-169) source and the Elekta Iridium (Ir-192) HDR Flexisource were modeled according to the TG-43U1 formalism and benchmarked accordingly. Twenty-two retrospective HDR prostate brachytherapy patients treated with Ir-192 were considered. An Ir-192 only (IRO), Yb-169 only (YBO), and dual-source (DS) plan with optimized catheter location was created for each patient with N catheters, where N is the number of catheters used in the clinically delivered plans. The DS plans jointly optimized Yb-169 and Ir-192 dwell times. All plans and the clinical plans were normalized to deliver a 15 Gy prescription (Rx) dose to 95% of the clinical treatment volume (CTV) and evaluated for the CTV D90%, V150%, and V200%, urethra D0.1cc and D1cc, bladder V75%, and rectum V75%. Dose-volume histograms (DVHs) were generated for each structure.

RESULTS

The DS plans ubiquitously selected Ir-192 as the only treatment source. IRO outperformed YBO in organ at risk (OARs) OAR sparing, reducing the urethra D0.1cc and D1cc by 0.98% (p = 2.22 ∗ 10 - 9 $p\ = \ 2.22*{10^{ - 9}}$ ) and 1.09% (p = 1.22 ∗ 10 - 10 $p\ = \ 1.22*{10^{ - 10}}$ ) of the Rx dose, respectively, and reducing the bladder and rectum V75% by 0.09 (p = 0.0023 $p\ = \ 0.0023$ ) and 0.13 cubic centimeters (cc) (p = 0.033 $p\ = \ 0.033$ ), respectively. The YBO plans delivered a more homogenous dose to the CTV, with a smaller V150% and V200% by 3.20 (p = 4.67 ∗ 10 - 10 $p\ = \ 4.67*{10^{ - 10}}$ ) and 1.91 cc (p = 5.79 ∗ 10 - 10 $p\ = \ 5.79*{10^{ - 10}}$ ), respectively, and a lower CTV D90% by 0.49% (p = 0.0056 $p\ = \ 0.0056$ ) of the prescription dose. The IRO plans reduce the urethral D1cc by 2.82% (p = 1.38 ∗ 10 - 4 $p\ = \ 1.38*{10^{ - 4}}$ ) of the Rx dose compared to the clinical plans, at the cost of increased bladder and rectal V75% by 0.57 (p = 0.0022 $p\ = \ 0.0022$ ) and 0.21 cc (p = 0.019 $p\ = \ 0.019$ ), respectively, and increased CTV V150% by a mean of 1.46 cc (p = 0.010 $p\ = \ 0.010$ ) and CTV D90% by an average of 1.40% of the Rx dose (p = 8.80 ∗ 10 - 8 $p\ = \ 8.80*{10^{ - 8}}$ ). While these differences are statistically significant, the clinical differences between the plans are minimal.

CONCLUSIONS

The proposed analytical HDR planning algorithm integrates catheter and isotope selection with dwell time optimization for varying clinical goals, including urethra sparing. The planning method can guide HDR implants and identify promising isotopes for specific HDR clinical goals, such as target conformality or OAR sparing.

A unified path seeking algorithm for IMRT and IMPT beam orientation optimization

Objective. Fully automated beam orientation optimization (BOO) for intensity-modulated radiotherapy and intensity modulated proton therapy (IMPT) is gaining interest, since achieving optimal plan quality for an unknown number of fixed beam arrangements is tedious. Fast group sparsity-based optimization methods have been proposed to find the optimal orientation, but manual tuning is required to eliminate the exact number of beams from a large candidate set. Here, we introduce a fast, automated gradient descent-based path-seeking algorithm (PathGD), which performs fluence map optimization for sequentially added beams, to visualize the dosimetric benefit of one added field at a time.Approach. Several configurations of 2-4 proton and 5-15 photon beams were selected for three head-and-neck patients using PathGD, which was compared to group sparsity-regularized BOO solved with the fast iterative shrinkage-thresholding algorithm (GS-FISTA), and manually selected IMPT beams or one coplanar photon VMAT arc (MAN). Once beams were chosen, all plans were compared on computational efficiency, dosimetry, and for proton plans, robustness.Main results. With each added proton beam, Clinical Target Volume (CTV) and organs at risk (OAR) dosimetric cost improved on average across plans by [1.1%, 13.6%], and for photons, [0.6%, 2.0%]. Comparing algorithms, beam selection for PathGD was faster than GS-FISTA on average by 35%, and PathGD matched the CTV coverage of GS-FISTA plans while reducing OAR mean and maximum dose in all structures by an average of 13.6%. PathGD was able to improve CTV [Dmax, D95%] by [2.6%, 5.2%] and reduced worst-case [max, mean] dose in OARs by [11.1%, 13.1%].Significance. The benefit of a path-seeking algorithm is the beam-by-beam analysis of dosimetric cost. PathGD was shown to be most efficient and dosimetrically desirable amongst group sparsity and manual BOO methods, and highlights the sensitivity of beam addition for IMPT in particular.

Dose and dose rate objectives in Bragg peak and shoot-through beam orientation optimization for FLASH proton therapy

PURPOSE

The combined use of Bragg peak (BP) and shoot-through (ST) beams has previously been shown to increase the normal tissue volume receiving FLASH dose rates while maintaining dose conformality compared to conventional intensity-modulated proton therapy (IMPT) methods. However, the fixed beam optimization method has not considered the effects of beam orientation on the dose and dose rates. To maximize the proton FLASH effect, here, we incorporate dose rate objectives into our beam orientation optimization framework.

METHODS

From our previously developed group-sparsity dose objectives, we add upper and lower dose rate terms using a surrogate dose-averaged dose rate definition and solve using the fast-iterative shrinking threshold algorithm. We compare the dosimetry for three head-and-neck cases between four techniques: (1) spread-out BP IMPT (BP), (2) dose rate optimization using BP beams only (BP-DR), (3) dose rate optimization using ST beams only (ST-DR), and (4) dose rate optimization using combined BP and ST (BPST-DR), with the goal of sparing organs at risk without loss of tumor coverage and maintaining high dose rate within a 10 mm region of interest (ROI) surrounding the clinical target volume (CTV).

RESULTS

For BP, BP-DR, ST-DR, and BPST-DR, CTV homogeneity index and Dmax were found to be on average 0.886, 0.867, 0.687, and 0.936 and 107%, 109%, 135%, and 101% of prescription, respectively. Although ST-DR plans were not able to meet dosimetric standards, BPST-DR was able to match or improve either maximum or mean dose in the right submandibular gland, left and right parotids, constrictors, larynx, and spinal cord compared to BP plans. Volume of ROIs receiving greater than 40 Gy/s ( V γ 0 ) ${V_{\gamma 0}})$ was 51.0%, 91.4%, 95.5%, and 92.1% on average.

CONCLUSIONS

The dose rate techniques, particularly BPST-DR, were able to significantly increase dose rate without compromising physical dose compared with BP. Our algorithm efficiently selects beams that are optimal for both dose and dose rate.