Toward a pre-clinical irradiator using clinical infrastructure

Evaluation of plan quality and treatment efficiency in cranial stereotactic radiosurgery treatment plans with a variable source-to-axis distance

INTRODUCTION

Radiotherapy deliveries with dynamic couch motions that shorten the source-to-axis distance (SAD) on a C-arm linac have the potential to increase treatment efficiency through the increase of the effective dose rate. In this investigation, we convert clinically deliverable volumetric modulated arc therapy (VMAT) and dynamic conformal arc (DCA) plans for cranial radiosurgery into virtual isocenter plans through implementation of couch trajectories that maintain the target at a shortened but variable SAD throughout treatment.

MATERIALS AND METHODS

A randomly sampled population of patients treated with cranial radiosurgery from within the last three years were separated into groups with one, two, and three lesions. All plans had a single isocenter (regardless of number of targets), and a single prescription dose. Patient treatment plans were converted from their original delivery at a standard isocenter to a dynamic virtual isocenter in MATLAB. The virtual isocenter plan featured a variable isocenter position based upon the closest achievable source-to-target distance (referred to herein as a virtual source-to-axis distance [vSAD]) which avoided collision zones on a TrueBeam STx platform. Apertures were magnified according to the vSAD and monitor units at a given control point were scaled based on the inverse square law. Doses were calculated for the plans with a virtual isocenter in the Eclipse (v13.6.23) treatment planning system (TPS) and were compared with the clinical plans. Plan metrics (MU, Paddick conformity index, gradient index, and the volume receiving 12 Gy or more), normal brain dose-volume differences, as well as maximum doses received by organs at risk (OARs) were assessed. The values were compared between standard and virtual isocenter plans with Wilcoxon Sign Ranked tests to determine significance. A subset of the plans were mapped to the MAX-HD anthropomorphic phantom which contained an insert housing EBT3 GafChromic film and a PTW 31010 microion chamber for dose verification on a linac.

RESULTS

Delivering plans at a virtual isocenter resulted in an average reduction of 20.9% (p = 3×10^-6 ) and 20.6% (p = 3.0×10^-6 ) of MUs across all VMAT and all DCA plans, respectively. There was no significant change in OAR max doses received by plans delivered at a virtual isocenter. The low dose wash (1.0-2.0 Gy or 5-11% of the prescription dose) was increased (by approximately 20 cc) for plans with three lesions. This was equivalent to a 2.7%-3.8% volumetric increase in normal tissue receiving the respective dose level when comparing the plan with a virtual isocenter to a plan with a standard isocenter. Gamma pass rates with a 5%/1mm analysis criteria were 96.40% ± 2.90% and 95.07% ± 3.10% for deliveries at standard and virtual isocenter, respectively. Absolute point dose agreements were within -0.36% ± 3.45% and -0.55% ± 3.39% for deliveries at a standard and virtual isocenter, respectively. Potential time savings per arc were found to have linear relationship with the monitor units delivered per arc (savings of 0.009 s/MU with an r² = 0.866 when fit to plans with a single lesion).

CONCLUSIONS

Converting clinical plans at standard isocenter to a virtual isocenter design did not show any losses to plan quality while simultaneously improving treatment efficiency through MU reductions.

Region-of-interest intra-arc MV imaging to facilitate sub-mm positional accuracy with minimal imaging dose during treatment deliveries of small cranial lesions

Purpose

To automate the generation of region‐of‐interest (ROI) apertures for use with megavoltage imaging for online positional corrections during cranial stereotactic radiosurgery.

Materials and methods

Digitally reconstructed radiographs (DRRs) were created for a 3D‐printed skull phantom at 5 degree gantry angle increments for a three‐arc beam arrangement. At each angle, 3000 random rectangular apertures were generated, and 100 shifts on a grid were applied to the anatomy within the frame. For all shifts, the mutual information (MI) between the shifted and unshifted DRR was calculated to derive an average MI gradient. The top 10% of apertures that minimized registration errors were overlaid and discretely thresholded to generate imaging plans. Imaging was acquired with the skull while implementing simulated patient motion on a linac. Control point‐specific couch motions were derived to align the skull to its planned positioning.

Results

Apertures with a range of repositioning errors less than 0.1 mm possessed a 42% larger average MI gradient when compared with apertures with a range greater than 1 mm. Dose calculations with Monte Carlo exhibited an 84% reduction in the dose received by 50% of the skull with the 50% thresholded plan when compared to a constant 22 × 22 cm² imaging plan. For all different imaging plans (with and without motion), the calculated median 3D‐errors with respect to the tracking of a metal‐BB fiducial positioned at isocenter in the skull were sub‐mm except for the 80% thresholded plan.

Conclusions

Sub‐mm positional errors are achievable with couch motions derived from control point–specific ROI imaging. Smaller apertures that conform to an anatomical ROI can be utilized to minimize the imaging dose incurred at the expense of larger errors.

Investigating the impacts of intrafraction motion on dosimetric outcomes when treating small targets with virtual cones

Purpose

Intrafraction patient motion is a well‐documented phenomenon in radiation therapy. In stereotactic radiosurgery applications in which target sizes can be very small and dose gradients very steep, patient motion can significantly impact the magnitude and positional accuracy of the delivered dose. This work investigates the impact of intrafraction motion on dose metrics for small targets when treated with a virtual cone.

Materials and Methods

Monte Carlo simulations were performed to calculate dose kernels for treatment apertures ranging from 1 × 2.5 mm² to 10 × 10 mm². The phantom was an 8.2‐cm diameter sphere and isotropic voxels had lengths of 0.25 mm. Simulated treatments consisted of 3 arcs: 1 axial arc (360° gantry rotation, couch angle 0°) and 2 oblique arcs (180° gantry rotation, couch angle ±45°). Dose distributions were calculated via superposition of the rotated kernels. Two different collimator orientations were considered to create a virtual cone: (a) each treatment arc was delivered twice, once each with a static collimator angle of ±45°, and (b) each treatment arc was delivered once, with dynamic collimator rotation throughout the arc. Two different intrafraction motion patterns were considered: (a) constant linear motion and (b) sudden, persistent motion. The impact of motion on dose distributions for target sizes ranging from 1 to 10 mm diameter spheres was quantified as a function of the aperture size used to treat the lesions.

Results

The impact of motion on both the target and the surrounding tissue was a function of both aperture shape and target size. When a 0.5‐mm linear drift along each dimension occurred during treatment, targets ≥5 mm saw less than a 10% decrease in coverage by the prescription dose. Smaller apertures accrued larger penalties with respect to dosimetric hotspots seen in the tissues surrounding the target volume during intrafraction motion. For example, treating a 4‐mm‐sized target that undergoes 2.60 mm (3D vector) of continuous linear motion, the D₅ in the concentric shells that extend 1, 2, and 3 mm from the surface of the target was 39%, 24%, and 14% smaller, respectively when comparing the delivery of a larger aperture (6 × 10 mm²) to a smaller aperture (2 × 5 mm²). Using a static collimator for shaping a virtual cone during treatment minimized the dosimetric impact of motion in the majority of cases. For example, the volume that is covered by 70% or more of the prescription dose is smaller in 60.4% of cases when using the static collimator. The volume covered by 50, and 30% or more of the prescription dose is also smaller when treating with a static collimator, but the clinical significance of this finding is unknown.

Conclusions

In this work, the dosimetric trade‐offs between aperture size and target size when irradiating with virtual cones has been demonstrated. These findings provide information about the tradeoffs between target coverage and normal tissue sparing that may help inform clinical decision making when treating smaller targets with virtual cones.

Small Animal IMRT Using 3D-Printed Compensators

PURPOSE

Preclinical radiation replicating clinical intensity modulated radiation therapy (IMRT) techniques can provide data translatable to clinical practice. For this work, treatment plans were created for oxygen-guided dose-painting in small animals using inverse-planned IMRT. Spatially varying beam intensities were achieved using 3-dimensional (3D)-printed compensators.

METHODS AND MATERIALS

Optimized beam fluence from arbitrary gantry angles was determined using a verified model of the XRAD225Cx treatment beam. Compensators were 3D-printed with varied thickness to provide desired attenuation using copper/polylactic-acid. Spatial resolution capabilities were investigated using printed test-patterns. Following American Association of Physicists in Medicine TG119, a 5-beam IMRT plan was created for a miniaturized (∼1/8th scale) C-shape target. Electron paramagnetic resonance imaging of murine tumor oxygenation guided simultaneous integrated boost (SIB) plans conformally treating tumor to a base dose (Rx₁) with boost (Rx₂) based on tumor oxygenation. The 3D-printed compensator intensity modulation accuracy and precision was evaluated by individually delivering each field to a phantom containing radiochromic film and subsequent per-field gamma analysis. The methodology was validated end-to-end with composite delivery (incorporating 3D-printed tungsten/polylactic-acid beam trimmers to reduce out-of-field leakage) of the oxygen-guided SIB plan to a phantom containing film and subsequent gamma analysis.

RESULTS

Resolution test-patterns demonstrate practical printer resolution of ∼0.7 mm, corresponding to 1.0 mm bixels at the isocenter. The miniaturized C-shape plan provides planning target volume coverage (V_95% = 95%) with organ sparing (organs at risk D_max < 50%). The SIB plan to hypoxic tumor demonstrates the utility of this approach (hypoxic tumor V_95%,Rx2 = 91.6%, normoxic tumor V_95%,Rx1 = 95.7%, normal tissue V_100%,Rx1 = 7.1%). The more challenging SIB plan to boost the normoxic tumor rim achieved normoxic tumor V_95%,Rx2 = 90.9%, hypoxic tumor V_95%,Rx1 = 62.7%, and normal tissue V_100%,Rx2 = 5.3%. Average per-field gamma passing rates using 3%/1.0 mm, 3%/0.7 mm, and 3%/0.5 mm criteria were 98.8% ± 2.8%, 96.6% ± 4.1%, and 90.6% ± 5.9%, respectively. Composite delivery of the hypoxia boost plan and gamma analysis (3%/1 mm) gave passing results of 95.3% and 98.1% for the 2 measured orthogonal dose planes.

CONCLUSIONS

This simple and cost-effective approach using 3D-printed compensators for small-animal IMRT provides a methodology enabling preclinical studies that can be readily translated into the clinic. The presented oxygen-guided dose-painting demonstrates that this methodology will facilitate studies driving much needed biologic personalization of radiation therapy for improvements in patient outcomes.