Technical Note: Improving the workflow in a carbon ion therapy center with custom software for enhanced patient care

Carbon-ion radiation therapy (CIRT) is an up-and-coming modality for cancer treatment. Implementation of CIRT requires collaboration among specialists like radiation oncologists, medical physicists, and other healthcare professionals. Effective communication among team members is necessary for the success of CIRT. However, the current workflows involving data management, treatment planning, scheduling, and quality assurance (QA) can be susceptible to errors, leading to delays and decreased efficiency. With the aim of addressing these challenges, a team of medical physicists developed an in-house workflow management software using FileMaker Pro. This tool has streamlined the workflow and improved the efficiency and quality of patient care.

Technical note: A comparison of in-house 3D-printed and commercially available patient-specific skin collimators for use with electron beam therapy

PURPOSE

Skin collimation is a useful tool in electron beam therapy (EBT) to decrease the penumbra at the field edge and minimize dose to nearby superficial organs at risk (OARs), but manually fabricating these collimation devices in the clinic to conform to the patient's anatomy can be a difficult and time intensive process. This work compares two types of patient-specific skin collimation (in-house 3D printed and vendor-provided machined brass) using clinically relevant metrics.

METHODS

Attenuation measurements were performed to determine the thickness of each material needed to adequately shield both 6 and 9 MeV electron beams. Relative and absolute dose planes at various depths were measured using radiochromic film to compare the surface dose, flatness, and penumbra of the different skin collimation materials.

RESULTS

Clinically acceptable thicknesses of each material were determined for both 6 and 9 MeV electron beams. Field width, flatness, and penumbra results between the two systems were very similar and significantly improved compared to measurements performed with no surface collimation.

CONCLUSION

Both skin collimation methods investigated in this work generate sharp penumbras at the field edge and can minimize dose to superficial OARs compared to treatment fields with no surface collimation. The benefits of skin collimation are greatest for lower energy electron beams, and the benefits decrease as the measurement depth increases. Using bolus with skin collimation is recommended to avoid surface dose enhancement seen with collimators placed on the skin surface. Ultimately, the appropriate choice of material will depend on the desire to create these devices in-house or outsource the fabrication to a vendor.

Impact of proton PBS machine operating parameters on the effectiveness of layer rescanning for interplay effect mitigation in lung SBRT treatment

BACKGROUND

Rescanning is a common technique used in proton pencil beam scanning to mitigate the interplay effect. Advances in machine operating parameters across different generations of particle therapy systems have led to improvements in beam delivery time (BDT). However, the potential impact of these improvements on the effectiveness of rescanning remains an underexplored area in the existing research.

METHODS

We systematically investigated the impact of proton machine operating parameters on the effectiveness of layer rescanning in mitigating interplay effect during lung SBRT treatment, using the CIRS phantom. Focused on the Hitachi synchrotron particle therapy system, we explored machine operating parameters from our institution's current (2015) and upcoming systems (2025A and 2025B). Accumulated dynamic 4D dose were reconstructed to assess the interplay effect and layer rescanning effectiveness.

RESULTS

Achieving target coverage and dose homogeneity within 2% deviation required 6, 6, and 20 times layer rescanning for the 2015, 2025A, and 2025B machine parameters, respectively. Beyond this point, further increasing the number of layer rescanning did not further improve the dose distribution. BDTs without rescanning were 50.4, 24.4, and 11.4 s for 2015, 2025A, and 2025B, respectively. However, after incorporating proper number of layer rescanning (six for 2015 and 2025A, 20 for 2025B), BDTs increased to 67.0, 39.6, and 42.3 s for 2015, 2025A, and 2025B machine parameters. Our data also demonstrated the potential problem of false negative and false positive if the randomness of the respiratory phase at which the beam is initiated is not considered in the evaluation of interplay effect.

CONCLUSION

The effectiveness of layer rescanning for mitigating interplay effect is affected by machine operating parameters. Therefore, past clinical experiences may not be applicable to modern machines.

Beam delivery characteristics of the Hitachi carbon ion scanning system at Osaka Heavy Ion Medical Accelerator in Kansai (HIMAK)

BACKGROUND

Using the pencil beam raster scanning method employed at most carbon beam treatment facilities, spots can be moved without interrupting the beam, allowing for the delivery of a dose between spots (move dose). This technique is also known as Dose-Driven-Continuous-Scanning (DDCS). To minimize its impact on HIMAK patient dosimetry, there's an upper limit to the move dose. Spots within a layer are grouped into sets, or "break points," allowing continuous irradiation. The beam is turned off when transitioning between sets or at the end of a treatment layer or spill. The control system beam-off is accomplished by turning off the RF Knockout (RFKO) extraction and after a brief delay the High Speed Steering Magnet (HSST) redirects the beam transport away from isocenter to a beam dump.

PURPOSE

The influence of the move dose and beam on/off control on the dose distribution and irradiation time was evaluated by measurements never before reported and modelled for Hitachi Carbon DDCS.

METHOD

We conducted fixed-point and scanning irradiation experiments at three different energies, both with and without breakpoints. For fixed-point irradiation, we utilized a 2D array detector and an oscilloscope to measure beam intensity over time. The oscilloscope data enabled us to confirm beam-off and beam-on timing due to breakpoints, as well as the relative timing of the RFKO signal, HSST signal, and dose monitor (DM) signals. From these measurements, we analyzed and modelled the temporal characteristics of the beam intensity. We also developed a model for the spot shape and amplitude at isocenter occurring after the beam-off signal which we called flap dose and its dependence on beam intensity. In the case of scanning irradiation, we measured move doses using the 2D array detector and compared these measurements with our model.

RESULT

We observed that the most dominant time variation of the beam intensity was at 1 kHz and its harmonic frequencies. Our findings revealed that the derived beam intensity cannot reach the preset beam intensity when each spot belongs to different breakpoints. The beam-off time due to breakpoints was approximately 100 ms, while the beam rise time and fall time (tdecay ) were remarkably fast, about 10 ms and 0.2 ms, respectively. Moreover, we measured the time lag (tdelay ) of approximately 0.2 ms between the RFKO and HSST signals. Since tdelay ≈ tdecay at HIMAK then the HSST is activated after the residual beam intensity, resulting in essentially zero flap dose at isocenter from the HSST. Our measurements of the move dose demonstrated excellent agreement with the modelled move dose.

CONCLUSION

We conducted the first move dose measurement for a Hitachi Carbon synchrotron, and our findings, considering beam on/off control details, indicate that Hitachi's carbon synchrotron provides a stable beam at HIMAK. Our work suggests that measuring both move dose and flap dose should be part of the commissioning process and possibly using our model in the Treatment Planning System (TPS) for new facilities with treatment delivery control systems with higher beam intensities and faster beam-off control.

Characterizing Proton-Induced Biological Effects in a Mouse Spinal Cord Model: A Comparison of Bragg Peak and Entrance Beam Response in Single and Fractionated Exposures

PURPOSE

Proton relative biological effectiveness (RBE) is a dynamic variable influenced by factors like linear energy transfer (LET), dose, tissue type, and biological endpoint. The standard fixed proton RBE of 1.1, currently used in clinical planning, may not accurately represent the true biological effects of proton therapy (PT) in all cases. This uncertainty can contribute to radiation-induced normal tissue toxicity in patients. In late-responding tissues such as the spinal cord, toxicity can cause devastating complications. This study investigated spinal cord tolerance in mice subjected to proton irradiation and characterized the influence of fractionation on proton- induced myelopathy at entrance (ENT) and Bragg peak (BP) positions.

METHODS AND MATERIALS

Cervical spinal cords of 8-week-old C57BL/6J female mice were irradiated with single- or multi-fractions (18x) using lateral opposed radiation fields at 1 of 2 positions along the Bragg curve: ENT (dose-mean LET = 1.2 keV/μm) and BP (LET = 6.9 keV/μm). Mice were monitored over 1 year for changes in weight, mobility, and general health, with radiation-induced myelopathy as the primary biological endpoint. Calculations of the RBE of the ENT and BP curve (RBEENT/BP) were performed.

RESULTS

Single-fraction RBEENT/BP for 50% effect probability (tolerance dose (TD50), grade II paresis, determined using log-logistic model fitting) was 1.10 ± 0.06 (95% CI) and for multifraction treatments it was 1.19 ± 0.05 (95% CI). Higher incidence and faster onset of paralysis were seen in mice treated at the BP compared with ENT.

CONCLUSIONS

The findings challenge the universally fixed RBE value in PT, indicating up to a 25% mouse spinal cord RBEENT/BP variation for multifraction treatments. These results highlight the importance of considering fractionation in determining RBE for PT. Robust characterization of proton-induced toxicity, aided by in vivo models, is paramount for refining clinical decision-making and mitigating potential patient side effects.

Variable RBE in proton radiotherapy: a comparative study with the predictive Mayo Clinic Florida microdosimetric kinetic model and phenomenological models of cell survival

Objectives. (1) To examine to what extent the cell- and exposure- specific information neglected in the phenomenological proton relative biological effectiveness (RBE) models could influence the computed RBE in proton therapy. (2) To explore similarities and differences in the formalism and the results between the linear energy transfer (LET)-based phenomenological proton RBE models and the microdosimetry-based Mayo Clinic Florida microdosimetric kinetic model (MCF MKM). (3) To investigate how the relationship between the RBE and the dose-mean proton LET is affected by the proton energy spectrum and the secondary fragments.Approach. We systematically compared six selected phenomenological proton RBE models with the MCF MKM in track-segment simulations, monoenergetic proton beams in a water phantom, and two spread-out Bragg peaks. A representative comparison within vitrodata for human glioblastoma cells (U87 cell line) is also included.Main results. Marked differences were observed between the results of the phenomenological proton RBE models, as reported in previous studies. The dispersion of these models' results was found to be comparable to the spread in the MCF MKM results obtained by varying the cell-specific parameters neglected in the phenomenological models. Furthermore, while single cell-specific correlation between RBE and the dose-mean proton LET seems reasonable above 2 keVμm-1, caution is necessary at lower LET values due to the relevant contribution of secondary fragments. The comparison within vitrodata demonstrates comparable agreement between the MCF MKM predictions and the results of the phenomenological models.Significance. The study highlights the importance of considering cell-specific characteristics and detailed radiation quality information for accurate RBE calculations in proton therapy. Furthermore, these results provide confidence in the use of the MCF MKM for clonogenic survival RBE calculations in proton therapy, offering a more mechanistic approach compared to phenomenological models.

Implementing dispersion measurement as part of scanning proton therapy commissioning and quality assurance

Introduction. Dispersion in an accelerator quantifies the deviation of the proton trajectory when there is a momentum deviation. We present for the first time a safe method of measuring dispersion in the clinic, using a scintillator detector and the momentum deviations within a spill. This is an important accelerator quantity as we found that this is the reason behind the large dose fluctuation in our absolute dosimetry measurement.Methods. Dispersions are measured for nine energies in a Hitachi ProBeat system at three locations (isocenter and at two profile monitors) and at two gantry angles (0 and 90 degrees) by first measuring the spot position and momentum drift within a spill. The spot position drift is measured by the XRV-4000 at the isocenter, and by the two profile monitors located at 0.57 and 2.27 m from the isocenter. The momentum drift is calculated from the intra-spill range drift which is measured using the Ranger accessory. The dispersion at isocenter and its gradient are calculated using the weighted least square regression on the measured dispersions at the three locations. A constraint is formulated on the dispersion and its gradient to ensure minimal intra-spill spot position deviation around the isocenter.Results. The measured intra-spill range and spot positional drift at isocenter are less than0.25mmand0.7mmrespectively. The momentum spread calculated from the range drift are less than 0.08%. The dispersion at the isocenter ranged from0.50to4.30mand the zero-crossing happens upstream of isocenter for all energies. 2 of the 9 energies (168.0 and 187.5 MeV) violated the constraint and has an intra-spill spot positional deviation greater than1.0within5cmfrom the isocenter.Conclusion. This measurement is recommended as part of commissioning and annual quality assurance for accelerator monitoring and to ensure intra-spill spot deviations remain low.

The effect of fitting the reference photon dose-response on the clonogenic survival predicted with the Mayo Clinic Florida microdosimetric kinetic model in case of accelerated ions

The Mayo Clinic Florida microdosimetric kinetic model (MCF MKM) is a recently developed clonogenic survival model. Since the MCF MKM relies on novel strategies to a priori determine the cell-specific model parameters, the only experiment-specific input values are the α and ß terms of the linear-quadratic model (LQM) of clonogenic survival for the reference photon exposure. Because the two LQM terms are anti-correlated, the fitting process of the reference photon survival curve was found to significantly influence the MCF MKM calculations. This article reports this effect for two clinically relevant cell lines (human brain glioblastoma A-172, human healthy foreskin fibroblasts AG01522) and ions (1H and 12C ions).

A high-resolution cone beam computed tomography (HRCBCT) reconstruction framework for CBCT-guided online adaptive therapy

BACKGROUND

Kilo-voltage cone-beam computed tomography (CBCT) is a prevalent modality used for adaptive radiotherapy (ART) due to its compatibility with linear accelerators and ability to provide online imaging. However, the widely-used Feldkamp-Davis-Kress (FDK) reconstruction algorithm has several limitations, including potential streak aliasing artifacts and elevated noise levels. Iterative reconstruction (IR) techniques, such as total variation (TV) minimization, dictionary-based methods, and prior information-based methods, have emerged as viable solutions to address these limitations and improve the quality and applicability of CBCT in ART.

PURPOSE

One of the primary challenges in IR-based techniques is finding the right balance between minimizing image noise and preserving image resolution. To overcome this challenge, we have developed a new reconstruction technique called high-resolution CBCT (HRCBCT) that specifically focuses on improving image resolution while reducing noise levels.

METHODS

The HRCBCT reconstruction technique builds upon the conventional IR approach, incorporating three components: the data fidelity term, the resolution preservation term, and the regularization term. The data fidelity term ensures alignment between reconstructed values and measured projection data, while the resolution preservation term exploits the high resolution of the initial Feldkamp-Davis-Kress (FDK) algorithm. The regularization term mitigates noise during the IR process. To enhance convergence and resolution at each iterative stage, we applied Iterative Filtered Backprojection (IFBP) to the data fidelity minimization process.

RESULTS

We evaluated the performance of the proposed HRCBCT algorithm using data from two physical phantoms and one head and neck patient. The HRCBCT algorithm outperformed all four different algorithms; FDK, Iterative Filtered Back Projection (IFBP), Compressed Sensing based Iterative Reconstruction (CSIR), and Prior Image Constrained Compressed Sensing (PICCS) methods in terms of resolution and noise reduction for all data sets. Line profiles across three line pairs of resolution revealed that the HRCBCT algorithm delivered the highest distinguishable line pairs compared to the other algorithms. Similarly, the Modulation Transfer Function (MTF) measurements, obtained from the tungsten wire insert on the CatPhan 600 physical phantom, showed a significant improvement with HRCBCT over traditional algorithms.

CONCLUSION

The proposed HRCBCT algorithm offers a promising solution for enhancing CBCT image quality in adaptive radiotherapy settings. By addressing the challenges inherent in traditional IR methods, the algorithm delivers high-definition CBCT images with improved resolution and reduced noise throughout each iterative step. Implementing the HR CBCT algorithm could significantly impact the accuracy of treatment planning during online adaptive therapy.

Investigation of scan path optimization in improving proton pencil beam scanning continuous delivery

Objective. To investigate the impact of scan path optimization on the dose accuracy and beam delivery time (BDT) of proton pencil beam scanning in the dose-driven continuous scanning (DDCS).Approach. A diverse set of six clinical plans, representing various spot patterns and treatment sites, was used to evaluate the effectiveness of scan time optimization and scan length optimization. The DDCS dose discrepancy and BDT with optimized scan paths was compared to the default serpentine scan path.Main results. Both scan time optimization and scan path optimization were able to reduce the DDCS dose discrepancy compared to the default serpentine scan path. All plans, except for the layer repainting lung plan, achieved a 2%/2 mm gamma pass rate of over 99% and less than 1% PTV DVH root mean square error (RMSE) through scan path optimization. In the case of the layer repainting lung plan, when compared to the default serpentine scan path, the 2%/2 mm gamma pass rate showed improvements from 91.3% to 93.1% and 95.8%, while the PTV DVH RMSE decreased from 2.1% to 1.7% and 1.1% for scan time optimization and scan length optimization, respectively. Although scan time optimization resulted in shorter total scan times for all plans compared to the default scan path and scan length optimization tended to have longer total scan times. However, due to the short total scan times and their minimal contribution to the total BDT, the impact of scan path optimization on the total BDT was practically negligible.Significance. Both scan time optimization and scan length optimization proved to be effective in minimizing DDCS dose discrepancy. No definitive winner can be determined between these two optimization approaches. Both scan time and scan length optimization had minimal effect on the total BDT.

Requirements for dose calculation on an active scanned proton beamline for small, shallow fields

INTRODUCTION

A growing interest in using proton pencil beam scanning in combination with collimators for the treatment of small, shallow targets, such as ocular melanoma or pre-clinical research emerged recently. This study aims at demonstrating that the dose of a synchrotron-based PBS system with a dedicated small, shallow field nozzle can be accurately predicted by a commercial treatment planning system (TPS) following appropriate tuning of both, nozzle and TPS.

MATERIALS

A removable extension to the clinical nozzle was developed to modify the beam shape passively. Five circular apertures with diameters between 5 to 34mm, mounted 72cm downstream of a range shifter were used. For each collimator treatment plans with spread-out Bragg peaks (SOBP) with a modulation of 3 to 30mm were measured and calculated with GATE/Geant4 and the research TPS RayStation (RS11B-R). The dose grid, multiple coulomb scattering and block discretization resolution were varied to find the optimal balance between accuracy and performance.

RESULTS

For SOBPs deeper than 10mm, the dose in the target agreed within 1% between RS11B-R, GATE/Geant4 and measurements for aperture diameters between 8 to 34mm, but deviated up to 5% for smaller apertures. A plastic taper was introduced reducing scatter contributions to the patient (from the pipe) and improving the dose calculation accuracy of the TPS to a 5% level in the entrance region for large apertures.

CONCLUSION

The commercial TPS and GATE/Geant4 can accurately calculate the dose for shallow, small proton fields using a collimator and pencil beam scanning.

Selecting Optimal Proton Pencil Beam Scanning Plan Parameters to Reduce Dose Discrepancy between Discrete Spot Plan and Continuous Scanning: A Proof-of-Concept Study

Pencil beam scanning delivered with continuous scanning has several advantages over conventional discrete spot scanning. Such advantages include improved beam delivery efficiency and reduced beam delivery time. However, a move dose is delivered between consecutive spots with continuous scanning, and current treatment planning systems do not take this into account. Therefore, continuous scanning and discrete spot plans have an inherent dose discrepancy. Using the operating parameters of the state-of-the-art particle therapy system, we conducted a proof-of-concept study in which we systematically generated 28 plans for cubic targets with different combinations of plan parameters and simulated the dose discrepancies between continuous scanning and a planned one. A nomograph to guide the selection of plan parameters was developed to reduce the dose discrepancy. The effectiveness of the nomograph was evaluated with two clinical cases (one prostate and one liver). Plans with parameters guided by the nomograph decreased dose discrepancy than those used standard plan parameters. Specifically, the 2%/2 mm gamma passing rate increased from 96.3% to 100% for the prostate case and from 97.8% to 99.7% for the liver case. The CTV DVH root mean square error decreased from 2.2% to 0.2% for the prostate case and from 1.8% to 0.9% for the liver case. The decreased dose discrepancy may allow the relaxing of the delivery constraint for some cases, leading to greater benefits in continuous scanning. Further investigation is warranted.

Technical note: Progressive deep learning: An accelerated training strategy for medical image segmentation

BACKGROUND

Recent advancements in Deep Learning (DL) methodologies have led to state-of-the-art performance in a wide range of applications especially in object recognition, classification, and segmentation of medical images. However, training modern DL models requires a large amount of computation and long training times due to the complex nature of network structures and the large number of training datasets involved. Moreover, it is an intensive, repetitive manual process to select the optimized configuration of hyperparameters for a given DL network.

PURPOSE

In this study, we present a novel approach to accelerate the training time of DL models via the progressive feeding of training datasets based on similarity measures for medical image segmentation. We term this approach Progressive Deep Learning (PDL).

METHODS

The two-stage PDL approach was tested on the auto-segmentation task for two imaging modalities: CT and MRI. The training datasets were ranked according to similarity measures between each sample based on Mean Square Error (MSE), Peak Signal-to-Noise Ratio (PSNR), Structural Similarity Index (SSIM), and the Universal Quality Image Index (UQI) values. At the start of the training process, a relatively coarse sampling of training datasets with higher ranks was used to optimize the hyperparameters of the DL network. Following this, the samples with higher ranks were used in step 1 to yield accelerated loss minimization in early training epochs and the total dataset was added in step 2 for the remainder of training.

RESULTS

Our results demonstrate that the PDL approach can reduce the training time by nearly half (∼49%) and can predict segmentations (CT U-net/DenseNet dice coefficient: 0.9506/0.9508, MR U-net/DenseNet dice coefficient: 0.9508/0.9510) without major statistical difference (Wilcoxon signed-rank test) compared to the conventional DL approach. The total training times with a fixed cutoff at 0.95 DSC for the CT dataset using DenseNet and U-Net architectures, respectively, were 17 h, 20 min and 4 h, 45 min in the conventional case compared to 8 h, 45 min and 2 h, 20 min with PDL. For the MRI dataset, the total training times using the same architectures were 2 h, 54 min and 52 min in the conventional case and 1 h, 14 min and 25 min with PDL.

CONCLUSION

The proposed PDL training approach offers the ability to substantially reduce the training time for medical image segmentation while maintaining the performance achieved in the conventional case.

Technical note: Delivery benefit and dosimetric implication of synchrotron-based proton pencil beam scanning using continuous scanning mode

BACKGROUND

Discrete spot scanning (DSS) is the commonly used method for proton pencil beam scanning (PBS). There is lack of data on the dose-driven continuous scanning (DDCS).

PURPOSE

To investigate delivery benefits and dosimetric implications of DDCS versus DSS for PBS systems.

METHODS

The irradiation duty factor, beam delivery time (BDT), and dose deviation were simulated for eight treatment plans in prostate, head and neck, liver, and lung, with both conventional fractionation and hypofractionation schemes. DDCS results were compared with those of DSS.

RESULTS

The DDCS irradiation duty factor (range, 11%-41%) was appreciably improved compared to DSS delivery (range, 4%-14%), within which, hypofractionation schemes had greater improvement than conventional fractionation. With decreasing stop ratio constraints, the DDCS BDT reduction was greater, but dose deviation also increased. With stop ratio constraints of 2, 1, 0.5, and 0, DDCS BDT reduction reached to 6%, 10%, 12%, and 15%, respectively, and dose deviation reached to 0.6%, 1.7%, 3.0%, and 5.2% root mean square error in PTV DVH, respectively. The 3%/2-mm gamma passing rate was greater than 99% with stop ratio constraints of 2 and 1, and greater than 95% with a stop ratio of 0.5. When the stop ratio constraint was removed, five of the eight treatment plans had a 3%/2-mm gamma passing rate greater than 95%, and the other three plans had a 3%/2-mm gamma passing rate between 90% and 95%.

CONCLUSIONS

The irradiation duty factor was considerably improved with DDCS. Smaller stop ratio constraints led to shorter BDTs, but with the cost of larger dose deviations. Our finding suggested that a stop ratio of 1 constraint seems to yield acceptable DDCS dose deviation.

The effect of spill change on reliable absolute dosimetry in a synchrotron proton spot scanning system

BACKGROUND

Absolute dosimetry measurement is an integral part of Treatment Planning System (TPS) commissioning and it involves measuring the integrated absorbed dose to water for all energies in a pencil beam scanning delivery system. During the commissioning of Singapore's first proton therapy center, a uniform scanned field with an Advanced Markus chamber method was employed for this measurement, and a large dose fluctuation of at least 5% was observed for 10% of the energy layers during repeated measurements.

PURPOSE

This study aims to understand the root cause of this fluctuation by relating the actual delivered spot information in the log file with the charge measurement by the ion chambers.

METHODS

A dedicated pencil beam dose algorithm was developed, taking into account the log file parameters, to calculate the dose for a single energy layer in a homogeneous water phantom. Three energies, 70.2, 182.7, and 228.7 MeV were used in this study, with the 182.7 MeV energy exhibiting large dose fluctuation. The dose fluctuation was investigated as a function of detector's sizes (pinpoint 3D, Advanced Markus, PTW 34070, and PTW 34089) and water depth (2 , 6, and 20 cm). Twelve ion chambers measurements were performed for each chamber and depth. The comparison of the theoretically predicted integrated dose and the charge measurement served as a validation of the algorithm.

RESULTS

About 5.9% and 9.6% dose fluctuation were observed in Advanced Markus and pinpoint 3D measurements at 2 cm depth for 182.7 MeV, while fluctuation of 1.6% and 1.1% were observed in Advanced Markus with 228.7 and 70.2 MeV at similar depth. Fluctuation of less than 0.1% was observed for PTW34070 and PTW 34089 for all energies. The fluctuation was found to diminish with larger spot size at 20 cm depth to 1.3% for 182.7 MeV. The theoretical and measured charge comparison showed a high linear correlation of R 2 > 0.80 ${R^2} > 0.80$ for all datasets, indicating the fluctuation originated from the delivered spot characteristics. The cause of fluctuation was identified to be due to the spill change occurring close to the detector, and since the spot positional deviation profiles were different between two spills, this resulted in local hot spots between columns of spots. The actual position of spill change varies randomly during measurement, which led to a random occurrence of hot spot within the detector's sensitive volume and a fluctuating dose measurement.

CONCLUSION

This is the first report of a dose fluctuation greater than 5% in absolute dosimetry measurement with a uniform scanned field and the cause of the fluctuation has been conclusively determined. It is important to choose the MU and scanning pattern carefully to avoid spill change happening when the spot delivery is near the detector.

A methodology to abridge microdosimetric distributions without a significant loss of the spectral information needed for the RBE computation in carbon ion therapy

BACKGROUND

In order to compute the relative biological effectiveness (RBE) of ion radiation therapy with the Mayo Clinic Florida microdosimetric kinetic model (MCF MKM), it is necessary to process entire microdosimetric distributions. Therefore, a posteriori RBE recalculations (i.e., for a different cell line or another biological endpoint) would require whole spectral information. It is currently not practical to compute and store all this data for each clinical voxel.

PURPOSE

To develop a methodology that allows to store a limited amount of physical information without losing accuracy in the RBE calculations nor the possibility of a posteriori RBE recalculations.

METHODS

Computer simulations for four monoenergetic ¹² C ion beams and a ¹² C ion spread-out Bragg peak (SOBP) were performed to assess lineal energy distributions as a function of the depth within a water phantom. These distributions were used in combination with the MCF MKM to compute the in vitro clonogenic survival RBE for human salivary gland tumor cells (HSG cell line) and human skin fibroblasts (NB1RGB cell line). The RBE values were also calculated with a new abridged microdosimetric distribution methodology (AMDM) and compared with the reference RBE calculations using the entire distributions.

RESULTS

The maximum relative deviation between the RBE values computed using the entire distributions and the AMDM was 0.61% (monoenergetic beams) and 0.49% (SOBP) for the HSG cell line, while 0.45% (monoenergetic beams) and 0.26% (SOBP) for the NB1RGB cell line.

CONCLUSION

The excellent agreement between the RBE values computed using the entire lineal energy distributions and the AMDM represents a milestone for the clinical implementation of the MCF MKM.

Validation of clinical acceptability of deep-learning-based automated segmentation of organs-at-risk for head-and-neck radiotherapy treatment planning

Introduction

Organ-at-risk segmentation for head and neck cancer radiation therapy is a complex and time-consuming process (requiring up to 42 individual structure, and may delay start of treatment or even limit access to function-preserving care. Feasibility of using a deep learning (DL) based autosegmentation model to reduce contouring time without compromising contour accuracy is assessed through a blinded randomized trial of radiation oncologists (ROs) using retrospective, de-identified patient data.

Methods

Two head and neck expert ROs used dedicated time to create gold standard (GS) contours on computed tomography (CT) images. 445 CTs were used to train a custom 3D U-Net DL model covering 42 organs-at-risk, with an additional 20 CTs were held out for the randomized trial. For each held-out patient dataset, one of the eight participant ROs was randomly allocated to review and revise the contours produced by the DL model, while another reviewed contours produced by a medical dosimetry assistant (MDA), both blinded to their origin. Time required for MDAs and ROs to contour was recorded, and the unrevised DL contours, as well as the RO-revised contours by the MDAs and DL model were compared to the GS for that patient.

Results

Mean time for initial MDA contouring was 2.3 hours (range 1.6-3.8 hours) and RO-revision took 1.1 hours (range, 0.4-4.4 hours), compared to 0.7 hours (range 0.1-2.0 hours) for the RO-revisions to DL contours. Total time reduced by 76% (95%-Confidence Interval: 65%-88%) and RO-revision time reduced by 35% (95%-CI,-39%-91%). All geometric and dosimetric metrics computed, agreement with GS was equivalent or significantly greater (p<0.05) for RO-revised DL contours compared to the RO-revised MDA contours, including volumetric Dice similarity coefficient (VDSC), surface DSC, added path length, and the 95%-Hausdorff distance. 32 OARs (76%) had mean VDSC greater than 0.8 for the RO-revised DL contours, compared to 20 (48%) for RO-revised MDA contours, and 34 (81%) for the unrevised DL OARs.

Conclusion

DL autosegmentation demonstrated significant time-savings for organ-at-risk contouring while improving agreement with the institutional GS, indicating comparable accuracy of DL model. Integration into the clinical practice with a prospective evaluation is currently underway.

Pragmatic, Prospective Comparative Effectiveness Trial of Carbon Ion Therapy, Surgery, and Proton Therapy for the Management of Pelvic Sarcomas (Soft Tissue/Bone) Involving the Bone: The PROSPER Study Rationale and Design

Surgical treatment of pelvic sarcoma involving the bone is the standard of care but is associated with several sequelae and reduced functional quality of life (QOL). Treatment with photon and proton radiotherapy is associated with relapse. Carbon ion radiotherapy (CIRT) may reduce both relapse rates and treatment sequelae. The PROSPER study is a tricontinental, nonrandomized, prospective, three-arm, pragmatic trial evaluating treatments of pelvic sarcoma involving the bone. Patients aged at least 15 years are eligible for inclusion. Participants must have an Eastern Cooperative Oncology Group Performance Status score of two or less, newly diagnosed disease, and histopathologic confirmation of pelvic chordoma, chondrosarcoma, osteosarcoma, Ewing sarcoma with bone involvement, rhabdomyosarcoma (RMS) with bone involvement, or non-RMS soft tissue sarcoma with bone involvement. Treatment arms include (1) CIRT (n = 30) delivered in Europe and Asia, (2) surgical treatment with or without adjuvant radiotherapy (n = 30), and (3) proton therapy (n = 30). Arms two and three will be conducted at Mayo Clinic campuses in Arizona, Florida, and Minnesota. The primary end point is to compare the 1-year change in functional QOL between CIRT and surgical treatment. Additional comparisons among the three arms will be made between treatment sequelae, local control, and other QOL measures.

A Geant4 shielding design for the first US carbon multi-ion hybrid synchrotron facility

The Mayo Clinic Florida Integrated Oncology Building will be the home of the first spot-scanning only carbon/proton hybrid therapy system by Hitachi, Ltd. It will provide proton beams up to kinetic energies of 230 MeV and carbon beams up to 430 MeV n^-1for clinical deployment. To provide adequate radiation protection, the Geant4 (v10.6) Monte Carlo toolkit was utilized to quantify the ambient dose equivalent at a 10 mm depth (H*(10)) for photons and neutrons. To perform accurate calculations of the ambient dose equivalent, three-dimensional computer-aided design files of the entire planned facility were imported into Geant4, as well as certain particle system components such as the bending magnets, fast Faraday cup, and gantry. Particle fluence was scored using 60 cm diameter spheres, which were strategically placed throughout areas of interests. Analytical calculations were performed as first-pass design checks. Major shielding slabs were optimized using Geant4 simulations iteratively, with more than 20 alternative designs evaluated within Geant4. The 430 MeV n^-1carbon beams played the most significant role in concrete thickness Requirements. The primary wall thickness for the carbon fixed beam room is 4 meters. The presence of the proton gantry structure in the simulation caused the ambient dose equivalent to increase by around 67% at the maze entrance, but a decrease in the high energy beam transport corridor. All shielding primary and secondary goals for clinical operations were met per state regulation and national guidelines.

The Mayo Clinic Florida microdosimetric kinetic model of clonogenic survival: formalism and first benchmark against in vitro and in silico data

Objective. To develop a new model (Mayo Clinic Florida microdosimetric kinetic model, MCF MKM) capable of accurately describing the in vitro clonogenic survival at low and high linear energy transfer (LET) using single-event microdosimetric spectra in a single target. Methodology. The MCF MKM is based on the ‘post-processing average’ implementation of the non-Poisson microdosimetric kinetic model and includes a novel expression to compute the particle-specific quadratic-dependence of the cell survival with respect to dose (β of the linear-quadratic model). A new methodology to a priori calculate the mean radius of the MCF MKM subnuclear domains is also introduced. Lineal energy spectra were simulated with the Particle and Heavy Ion Transport code System (PHITS) for 1H, 4He, 12C, 20Ne, 40Ar, 56Fe, and 132Xe ions and used in combination with the MCF MKM to calculate the ion-specific LET-dependence of the relative biological effectiveness (RBE) for Chinese hamster lung fibroblasts (V79 cell line) and human salivary gland tumor cells (HSG cell line). The results were compared with in vitro data from the Particle Irradiation Data Ensemble (PIDE) and in silico results of different models. The possibility of performing experiment-specific predictions to explain the scatter in the in vitro RBE data was also investigated. Finally, a sensitivity analysis on the model parameters is also included. Main results. The RBE values predicted with the MCF MKM were found to be in good agreement with the in vitro data for all tested conditions. Though all MCF MKM model parameters were determined a priori, the accuracy of the MCF MKM was found to be comparable or superior to that of other models. The model parameters determined a priori were in good agreement with the ones obtained by fitting all available in vitro data. Significance. The MCF MKM will be considered for implementation in cancer radiotherapy treatment planning with accelerated ions.

Exploring the Biological and Physical Basis of Boron Neutron Capture Therapy (BNCT) as a Promising Treatment Frontier in Breast Cancer

Simple Summary

BNCT is a biologically targeted, densely ionizing form of radiation therapy that allows for increased tumor cell kill, while reducing toxicity to surrounding normal tissues. Although BNCT has been investigated in the treatment of head and neck cancers and recurrent brain tumors, its applicability to breast cancer has not been previoulsy investigated. In this review we discuss the physical and biological properties of various boronated compounds, and advantages and challenges associated with the potential use of BNCT in the treatment of breast cancer.

Abstract

BNCT is a high LET radiation therapy modality that allows for biologically targeted radiation delivery to tumors while reducing normal tissue impacts. Although the clinical use of BNCT has largely been limited to phase I/II trials and has primarily focused on difficult-to-treat malignancies such as recurrent head and neck cancer and recurrent gliomas, recently there has been a renewed interest in expanding the use of BNCT to other disease sites, including breast cancer. Given its high LET characteristics, its biologically targeted and tumor specific nature, as well as its potential for use in complex treatment settings including reirradiation and widespread metastatic disease, BNCT offers several unique advantages over traditional external beam radiation therapy. The two main boron compounds investigated to date in BNCT clinical trials are BSH and BPA. Of these, BPA in particular shows promise in breast cancer given that is taken up by the LAT-1 amino acid transporter that is highly overexpressed in breast cancer cells. As the efficacy of BNCT is directly dependent on the extent of boron accumulation in tumors, extensive preclinical efforts to develop novel boron delivery agents have been undertaken in recent years. Preclinical studies have shown promise in antibody linked boron compounds targeting ER/HER2 receptors, boron encapsulating liposomes, and nanoparticle-based boron delivery systems. This review aims to summarize the physical and biological basis of BNCT, the preclinical and limited clinical data available to date, and discuss its potential to be utilized for the successful treatment of various breast cancer disease states.

On the calculation of the relative biological effectiveness of ion radiation therapy using a biological weighting function, the microdosimetric kinetic model (MKM) and subsequent corrections (non-Poisson MKM and modified MKM)

Objective. To investigate similarities and differences in the formalism, processing, and the results of relative biological effectiveness (RBE) calculations with a biological weighting function (BWF), the microdosimetric kinetic model (MKM) and subsequent modifications (non-Poisson MKM, modified MKM). This includes: (a) the extension of the V79-RBE10% BWF to model the RBE for other clonogenic survival levels; (b) a novel implementation of MKMs as weighting functions; (c) a benchmark against Chinese Hamster lung fibroblast (V79) in vitro data; (d) a study on the effect of pre- or post- processing the average biophysical quantities used for the RBE calculations; (e) a possible modification of the modified MKM parameters to improve the model accuracy at high linear energy transfer (LET). Methodology. Lineal energy spectra were simulated for two spherical targets (diameter = 0.464 or 1.0 μm) using PHITS for 1H, 4He, 12C, 20Ne, 40Ar, 56Fe and 132Xe ions. The results of the in silico calculations were compared with published in vitro data. Main results. All models appear to underestimate the RBE α of hydrogen ions. All MKMs generally overestimate the RBE50%, RBE10% and RBE1% for ions with an LET greater than ∼200 keV μm−1. This overestimation is greater for small surviving fractions and is likely due to the assumption of a radiation-independent quadratic term of clonogenic survival (ß). The overall RBE trends seem to be best described by the novel ‘post-processing average’ implementation of the non-Poisson MKM. In case of calculations with the non-Poisson MKM, pre- or post- processing the average biophysical quantities affects the computed RBE values significantly. Significance. This study presents a systematic analysis of the formalism and results of widely used microdosimetric models of clonogenic survival for ions relevant for cancer particle therapy and space radiation protection. Points for improvements were highlighted and will contribute to the development of upgraded biophysical models.

RBE-weighted dose and its impact on the risk of acute coronary event for breast cancer patients treated with intensity modulated proton therapy

PURPOSE

To evaluate the relative biological effectiveness (RBE)-weighted dose to the heart and to estimate RBE uncertainties when assuming a constant RBE of 1.1, for breast cancer patients receiving intensity-modulated proton therapy (IMPT). Further, to study the impact of RBE uncertainties on the risk of an acute coronary event (ACE).

MATERIAL AND METHODS

We analyzed 20 patients who received IMPT to either the left breast (n = 10) or left chest wall (n = 10) and regional lymph nodes. The Monte Carlo simulation engine, MCsquare, was used to simulate the dose-averaged linear energy transfer (LETd) map. The RBE-weighted dose to the heart and its substructures was calculated using three different RBE models. The risk of ACE was estimated per its linear relationship with mean heart dose (MHD) as established by Darby et al.

RESULTS

The median MHD increased from 1.33 GyRBE assuming an RBE of 1.1 to 1.64, 1.87, and 1.99 GyRBE when using the RBE-weighted dose models. The median values (and ranges) of the excess absolute risk of ACE were 0.4% (0.1%-0.8%) when assuming an RBE of 1.1, and 0.6% (0.2%-1.0%), 0.6% (0.2%-1.1%), and 0.7% (0.2%-1.1%) with the RBE-weighted models. For our patient cohort, the maximum excess absolute risk of ACE increased by 0.3% with the RBE-weighted doses compared to the constant RBE of 1.1, reaching an excess absolute ACE risk of 1.1%. The interpatient LETd variation was small for the relevant high-dose regions of the heart.

CONCLUSION

All three RBE models predicted a higher biological dose compared to the clinical standard dose assuming a constant RBE of 1.1. An underestimation of the biological dose results in underestimation of the ACE risk. Analyzing the voxel-by-voxel biological dose and the LET map alongside clinical outcomes is warranted in the development of a more accurate normal-tissue complication probability model.

Design of a 3D patient-specific collision avoidance virtual framework for half-gantry proton therapy system

Introduction

This study presents a comprehensive collision avoidance framework based on three‐dimension (3D) computer‐aided design (CAD) modeling, a graphical user interface (GUI) as peripheral to the radiation treatment planning (RTP) environment, and patient‐specific plan parameters for intensity‐modulated proton therapy (IMPT).

Methods

A stand‐alone software application was developed leveraging the Varian scripting application programming interface (API) for RTP database object accessibility. The Collision Avoider software models the Hitachi ProBeat‐V half gantry design and the Kuka robotic couch with triangle mesh structures. Patient‐specific plan parameters are displayed in the collision avoidance software for potential proximity evaluation. The external surfaces of the patients and the immobilization devices are contoured based on computed tomography (CT) images. A “table junction‐to‐CT‐origin” (JCT) measurement is made for every patient at the time of CT simulation to accurately provide reference location of the patient contours to the treatment couch. Collision evaluations were performed virtually with the program during treatment planning to prevent four major types of collisional events: collisions between the gantry head and the treatment couch, gantry head and the patient's body, gantry head and the robotic arm, and collisions between the gantry head and the immobilization devices.

Results

The Collision Avoider software was able to accurately model the proton treatment delivery system and the robotic couch position. Commonly employed clinical beam configuration and JCT values were investigated. Brain and head and neck patients require more complex gantry and patient positioning system configurations. Physical measurements were performed to validate 3D CAD model geometry. Twelve clinical proton treatment plans were used to validate the accuracy of the software. The software can predict all four types of collisional events in our clinic since its full implementation in 2020.

Conclusion

A highly efficient patient‐specific collision prevention program for scanning proton therapy has been successfully implemented. The graphical program has provided accurate collision detection since its inception at our institution.

Nuclear Fragmentation Imaging for Carbon-Ion Radiation Therapy Monitoring: an In Silico Study

Purpose

This article presents an in vivo imaging technique based on nuclear fragmentation of carbon ions in irradiated tissues for potential real-time monitoring of carbon-ion radiation therapy (CIRT) treatment delivery and quality assurance purposes in clinical settings.

Materials and Methods

A proof-of-concept imaging and monitoring system (IMS) was devised to implement the technique. Monte Carlo simulations were performed for a prospective pencil-beam scanning CIRT nozzle. The development IMS benchmark considered a 5×5-cm² pixelated charged-particle detector stack positioned downstream from a target phantom and list-mode data acquisition. The abundance and production origins, that is, vertices, of the detected fragments were studied. Fragment trajectories were approximated by straight lines and a beam back-projection algorithm was built to reconstruct the vertices. The spatial distribution of the vertices was then used to determine plan relevant markers.

Results

The IMS technique was applied for a simulated CIRT case, a primary brain tumor. Four treatment plan monitoring markers were conclusively recovered: a depth dose distribution correlated profile, ion beam range, treatment target boundaries, and the beam spot position. Promising millimeter-scale (3-mm, ≤10% uncertainty) beam range and submillimeter (≤0.6-mm precision for shifts <3 cm) beam spot position verification accuracies were obtained for typical therapeutic energies between 150 and 290 MeV/u.

Conclusions

This work demonstrated a viable online monitoring technique for CIRT treatment delivery. The method's strong advantage is that it requires few signal inputs (position and timing), which can be simultaneously acquired with readily available technology. Future investigations will probe the technique's applicability to motion-sensitive organ sites and patient tissue heterogeneities. In-beam measurements with candidate detector-acquisition systems are ultimately essential to validate the IMS benchmark performance and subsequent deployment in the clinic.

Dosimetric Assessment of a High Precision System for Mouse Proton Irradiation to Assess Spinal Cord Toxicity

The uncertainty associated with the relative biological effectiveness (RBE) in proton therapy, particularly near the Bragg peak (BP), has led to the shift towards biological-based treatment planning. Proton RBE uncertainty has recently been reported as a possible cause for brainstem necrosis in pediatric patients treated with proton therapy. Despite this, in vivo studies have been limited due to the complexity of accurate delivery and absolute dosimetry. The purpose of this investigation was to create a precise and efficient method of treating the mouse spinal cord with various portions of the proton Bragg curve and to quantify associated uncertainties for the characterization of proton RBE. Mice were restrained in 3D printed acrylic boxes, shaped to their external contour, with a silicone insert extending down to mold around the mouse. Brass collimators were designed for parallel opposed beams to treat the spinal cord while shielding the brain and upper extremities of the animal. Up to six animals may be accommodated for simultaneous treatment within the restraint system. Two plans were generated targeting the cervical spinal cord, with either the entrance (ENT) or the BP portion of the beam. Dosimetric uncertainty was measured using EBT3 radiochromic film with a dose-averaged linear energy transfer (LETd) correction. Positional uncertainty was assessed by collecting a library of live mouse scans (n = 6 mice, two independent scans per mouse) and comparing the following dosimetric statistics from the mouse cervical spinal cord: Volume receiving 90% of the prescription dose (V90); mean dose to the spinal cord; and LETd. Film analysis results showed the dosimetric uncertainty to be ±1.2% and ±5.4% for the ENT and BP plans, respectively. Preliminary results from the mouse library showed the V90 to be 96.3 ± 4.8% for the BP plan. Positional uncertainty of the ENT plan was not measured due to the inherent robustness of that treatment plan. The proposed high-throughput mouse proton irradiation setup resulted in accurate dose delivery to mouse spinal cords positioned along the ENT and BP. Future directions include adapting the setup to account for weight fluctuations in mice undergoing fractionated irradiation.

Functional avoidance-based intensity modulated proton therapy with 4DCT derived ventilation imaging for lung cancer

The primary objective is to evaluate the potential dosimetric gains of performing functional avoidance‐based proton treatment planning using 4DCT derived ventilation imaging. 4DCT data of 31 patients from a prospective functional avoidance clinical trial were evaluated with intensity modulated proton therapy (IMPT) plans and compared with clinical volumetric modulated arc therapy (VMAT) plans. Dosimetric parameters were compared between standard and functional plans with IMPT and VMAT with one‐way analysis of variance and post hoc paired student t‐test. Normal Tissue Complication Probability (NTCP) models were employed to estimate the risk of two toxicity endpoints for healthy lung tissues. Dose degradation due to proton motion interplay effect was evaluated. Functional IMPT plans led to significant dose reduction to functional lung structures when compared with functional VMAT without significant dose increase to Organ at Risk (OAR) structures. When interplay effect is considered, no significant dose degradation was observed for the OARs or the clinical target volume (CTV) volumes for functional IMPT. Using fV20 as the dose metric and Grade 2+ pneumonitis as toxicity endpoint, there is a mean 5.7% reduction in Grade 2+ RP with the functional IMPT and as high as 26% in reduction for individual patient when compared to the standard IMPT planning. Functional IMPT was able to spare healthy lung tissue to avoid excess dose to normal structures while maintaining satisfying target coverage. NTCP calculation also shows that the risk of pulmonary complications can be further reduced with functional based IMPT.

Inhibition of ATM Induces Hypersensitivity to Proton Irradiation by Upregulating Toxic End Joining

Proton Bragg peak irradiation has a higher ionizing density than conventional photon irradiation or the entrance of the proton beam profile. Whether targeting the DNA damage response (DDR) could enhance vulnerability to the distinct pattern of damage induced by proton Bragg peak irradiation is currently unknown. Here, we performed genetic or pharmacologic manipulation of key DDR elements and evaluated DNA damage signaling, DNA repair, and tumor control in cell lines and xenografts treated with the same physical dose across a radiotherapy linear energy transfer spectrum. Radiotherapy consisted of 6 MV photons and the entrance beam or Bragg peak of a 76.8 MeV spot scanning proton beam. More complex DNA double-strand breaks (DSB) induced by Bragg peak proton irradiation preferentially underwent resection and engaged homologous recombination (HR) machinery. Unexpectedly, the ataxia-telangiectasia mutated (ATM) inhibitor, AZD0156, but not an inhibitor of ATM and Rad3-related, rendered cells hypersensitive to more densely ionizing proton Bragg peak irradiation. ATM inhibition blocked resection and shunted more DSBs to processing by toxic ligation through nonhomologous end-joining, whereas loss of DNA ligation via XRCC4 or Lig4 knockdown rescued resection and abolished the enhanced Bragg peak cell killing. Proton Bragg peak monotherapy selectively sensitized cell lines and tumor xenografts with inherent HR defects, and the repair defect induced by ATM inhibitor coadministration showed enhanced efficacy in HR-proficient models. In summary, inherent defects in HR or administration of an ATM inhibitor in HR-proficient tumors selectively enhances the relative biological effectiveness of proton Bragg peak irradiation. SIGNIFICANCE: Coadministration of an ATM inhibitor rewires DNA repair machinery to render cancer cells uniquely hypersensitive to DNA damage induced by the proton Bragg peak, which is characterized by higher density ionization.See related commentary by Nickoloff, p. 3156.

A High-Precision Method for In Vitro Proton Irradiation

Purpose

Although proton therapy has become a well-established radiation modality, continued efforts are needed to improve our understanding of the molecular and cellular mechanisms occurring during treatment. Such studies are challenging, requiring many resources. The purpose of this study was to create a phantom that would allow multiple in vitro experiments to be irradiated simultaneously with a spot-scanning proton beam.

Materials and Methods

The setup included a modified patient-couch top coupled with a high-precision robotic arm for positioning. An acrylic phantom was created to hold 4 6-well cell-culture plates at 2 different positions along the Bragg curve in a reproducible manner. The proton treatment plan consisted of 1 large field encompassing all 4 plates with a monoenergetic 76.8-MeV posterior beam. For robust delivery, a mini pyramid filter was used to broaden the Bragg peak (BP) in the depth direction. Both a Markus ionization chamber and EBT3 radiochromic film measurements were used to verify absolute dose.

Results

A treatment plan for the simultaneous irradiation of 2 plates irradiated with high linear energy transfer protons (BP, 7 keV/μm) and 2 plates irradiated with low linear energy transfer protons (entrance, 2.2 keV/μm) was created. Dose uncertainty was larger across the setup for cell plates positioned at the BP because of beam divergence and, subsequently, variable proton-path lengths. Markus chamber measurements resulted in uncertainty values of ±1.8% from the mean dose. Negligible differences were seen in the entrance region (<0.3%).

Conclusion

The proposed proton irradiation setup allows 4 plates to be simultaneously irradiated with 2 different portions (entrance and BP) of a 76.8-MeV beam. Dosimetric uncertainties across the setup are within ±1.8% of the mean dose.

The Importance of Verification CT-QA Scans in Patients Treated with IMPT for Head and Neck Cancers

Purpose

To understand how verification computed tomography-quality assurance (CT-QA) scans influenced clinical decision-making to replan patients with head and neck cancer and identify predictors for replanning to guide intensity-modulated proton therapy (IMPT) clinical practice.

Patients and Methods

We performed a quality-improvement study by prospectively collecting data on 160 consecutive patients with head and neck cancer treated using spot-scanning IMPT who underwent weekly verification CT-QA scans. Kaplan-Meier estimates were used to determine the cumulative probability of a replan by week. Predictors for replanning were determined with univariate (UVA) and multivariate (MVA) Cox model hazard ratios (HRs). Logistic regression was used to determine odds ratios (ORs). P < .05 was considered statistically significant.

Results

Of the 160 patients, 79 (49.4%) had verification CT-QA scans, which prompted a replan. The cumulative probability of a replan by week 1 was 13.7% (95% confidence interval [CI], 8.82-18.9), week 2, 25.0% (95% CI, 18.0-31.4), week 3, 33.1% (95% CI, 25.4-40.0), week 4, 45.6% (95% CI, 37.3-52.8), and week 5 and 6, 49.4% (95% CI, 41.0-56.6). Predictors for replanning were sinonasal disease site (UVA: HR, 1.82, P = .04; MVA: HR, 3.64, P = .03), advanced stage disease (UVA: HR, 4.68, P < .01; MVA: HR, 3.10, P < .05), dose > 60 Gy equivalent (GyE; relative biologic effectiveness, 1.1) (UVA: HR, 1.99, P < .01; MVA: HR, 2.20, P < .01), primary disease (UVA: HR, 2.00 versus recurrent, P = .01; MVA: HR, 2.46, P = .01), concurrent chemotherapy (UVA: HR, 2.05, P < .01; MVA: not statistically significant [NS]), definitive intent treatment (UVA: HR, 1.70 versus adjuvant, P < .02; MVA: NS), bilateral neck treatment (UVA: HR, 2.07, P = .03; MVA: NS), and greater number of beams (5 beam UVA: HR, 5.55 versus 1 or 2 beams, P < .02; MVA: NS). Maximal weight change from baseline was associated with higher odds of a replan (≥3 kg: OR, 1.97, P = .04; ≥ 5 kg: OR, 2.13, P = .02).

Conclusions

Weekly verification CT-QA scans frequently influenced clinical decision-making to replan. Additional studies that evaluate the practice of monitoring IMPT-treated patients with weekly CT-QA scans and whether that improves clinical outcomes are warranted.

Technical Delivery Parameters of 2000 Proton Treatment Courses

Purpose

To summarize the technical delivery parameters of proton plans delivered at the Mayo Clinic in Rochester, Minnesota.

Materials and Methods

The database of treated patient proton plans was queried to extract field parameters such as gantry angle, patient support angle, minimum and maximum water-equivalent depth (WED) treated, number of layers, field size, patient orientation, and monitor units. The plans were analyzed in aggregate, by disease site, and by fractionation.

Results

There were 2963 proton plans for 2023 distinct treatment sites delivered between June 2015 and September 2018. The mean number of fields per plan was 2.8. The mean number of energy layers per field was 51.9. The mean monitor unit per field was 117.4. The median maximum field dimension was 12.4 cm; 95% of the fields had a maximum dimension < 28.7 cm, and the maximum field dimension was 39.8 cm. The median maximum field WED was 16.4 cm; 95% of the fields reached a maximum WED of ≤ 26.4 cm, and the maximum field WED was 32.4 cm.

Conclusion

A large variety of disease sites were treated using the maximum field size (40 cm) and WED (32.4 cm) capabilities of our half-gantry system.

A Pareto-based beam orientation optimization method for spot scanning intensity-modulated proton therapy

PURPOSE

To provide a proof of principle of a Pareto-based method to automatically generate optimal intensity-modulated proton therapy (IMPT) plans for various noncoplanar beam orientations.

METHODS

A novel multicriteria beam orientation optimization (MCBOO) method was developed to generate Pareto database of optimal plans. The MCBOO method automatically explores the beam orientations and the scalarization parameters of the IMPT plans simultaneously. The MCBOO method is based on multicriteria bilevel optimization (i.e., hierarchical optimization with two nested levels, named the upper and lower level optimization). In MCBOO, the upper level optimization explores the noncoplanar beam orientation space, while the lower level explores the scalarization parameters for a given beam orientation. Differential evolution method was used in both levels, and the Pareto optimal plans were aggregated from the bilevel optimizations to construct the Pareto database. The MCBOO method was implemented on a multinode multi-GPU cluster, and it was tested on three brain tumor patient cases. The Pareto database of the three patients was generated for a set of DVH-based objectives. A statistical analysis was performed between a selected set of MCBOO plans and the manual plan (plan with manually selected beam orientation based on the clinical experience and optimized with the same single plan iterative optimizer used in the MCBOO). The selected set of MCBOO plans consisted of plans that matched the performance of the manual plan [i.e., MCBOO plans that have the same target coverage (within 2%) as the manual plan or better and achieved the same dose (within 2%) or lower to all of the organs at risks (OARs) but one OAR]. Additionally, a dosimetric comparison between of one of the selected MCBOO plans vs the manual plan was conducted.

RESULTS

The multicriteria beam orientation optimization algorithm automatically generated Pareto plans for the three noncoplanar brain tumor cases. The MCBOO plans provided an alternative objective trade-offs to the manual plan. The selected MCBOO plans showed a reduction in dose to multiple organs at risk vs the manual plan with a maximum value which ranged between 10.8 and 12.9 Gy for the three patients. The trade-off of the OAR dose reduction resulted in higher dose to no more than one OAR for each of the selected MCBOO plans vs the manual plan. The maximum dose increase in the MCBOO plans over the manual plan ranged from 7.8 to 11.8 Gy.

CONCLUSIONS

A novel multicriteria beam orientation optimization method was developed and tested on three IMPT patient cases. The method automatically generates Pareto plans database by exploring the noncoplanar beam orientations. The method was able to identify beam orientations with Pareto optimal plans that are comparable to the manually created plans with varying objective trade-offs.

Incorporation of Biologic Response Variance Modeling Into the Clinic: Limiting Risk of Brachial Plexopathy and Other Late Effects of Breast Cancer Proton Beam Therapy

Purpose:

The relative biologic effectiveness (RBE) rises with increasing linear energy transfer toward the end of proton tracks. Presently, there is no consensus on how RBE heterogeneity should be accounted for in breast cancer proton therapy treatment planning. Our purpose was to determine the dosimetric consequences of incorporating a brachial plexus (BP) biologic dose constraint and to describe other clinical implications of biologic planning.

Methods and Materials:

We instituted a biologic dose constraint for the BP in the context of MC1631, a randomized trial of conventional versus hypofractionated postmastectomy intensity modulated proton therapy (IMPT). IMPT plans of 13 patients treated before the implementation of the biologic dose constraint (cohort A) were compared with IMPT plans of 38 patients treated on MC1631 after its implementation (cohort B) using (1) a commercially available Eclipse treatment planning system (RBE = 1.1); (2) an in-house graphic processor unit-based Monte Carlo physical dose simulation (RBE = 1.1); and (3) an in-house Monte Carlo biologic dose (MCBD) simulation that assumes a linear relationship between RBE and dose-averaged linear energy transfer (product of RBE and physical dose = biologic dose).

Results:

Before implementation of a BP biologic dose constraint, the Eclipse mean BP D0.01 cm³ was 107%, and the MCBD estimate was 128% (ie, 64 Gy [RBE = biologic dose] in 25 fractions for a 50-Gy [RBE = 1.1] prescription), compared with 100.0% and 116.0%, respectively, after the implementation of the constraint. Implementation of the BP biologic dose constraint did not significantly affect clinical target volume coverage. MCBD plans predicted greater internal mammary node coverage and higher heart dose than Eclipse plans.

Conclusions:

Institution of a BP biologic dose constraint may reduce brachial plexopathy risk without compromising target coverage. MCBD plan evaluation provides valuable information to physicians that may assist in making clinical judgments regarding relative priority of target coverage versus normal tissue sparing.

Adaptive method for multicriteria optimization of intensity-modulated proton therapy

PURPOSE

Provide an adaptive multicriteria optimization (MCO) method for intensity-modulated proton therapy (IMPT) utilizing GPU technology. Previously described limitations of MCO such as Pareto approximation and limitation on the number of objectives were addressed.

METHODS

The treatment planning process for IMPT must account for multiple objectives, which requires extensive treatment planning resources. Often a large number of objectives (>10) are required. Hence the need for an MCO algorithm that can handle large number of objectives. The novelty of the MCO method presented here lies on the introduction of the adaptive weighting scheme that can generate a well-distributed and dense representation of the Pareto surface for a large number of objectives in an efficient manner. In our approach the generated Pareto surface is constructed for a set of DVH objectives. The MCO algorithm is based on the augmented weighted Chebychev metric (AWCM) method with an adaptive weighting scheme. This scheme uses the differential evolution (DE) method to generate a set of well-distributed Pareto points. The quality of the Pareto points' distribution in the objective space was assessed quantitatively using the Pareto sampling metric. The MCO algorithm was developed to perform multiple parallel searches to achieve a rapid mapping of the Pareto surface, produce clinically deliverable plans, and was implemented on a GPU cluster. The MCO algorithm was tested on two clinical cases with 10 and 18 objectives. For each case one of the MCO-generated plans was selected for comparison with the clinically generated plan. The MCO plan was randomly selected out of the set of MCO plans that had target coverage similar to the clinically generated plan and the same or better sparing of the organs at risk (OAR). Additionally, a validation study of the AWCM method vs the weighted sum method was performed.

RESULTS

The adaptive MCO algorithm generated Pareto points on the Pareto hypersurface in a fast (2-3 hr) and efficient manner for 2 cases with 10 and 18 objectives. The MCO algorithm generated a dense and well-distributed set of Pareto points on the objective space, and was able to achieve minimization of the Pareto sampling metric. The selected MCO plan showed an improvement of the DVH objectives in comparison to the clinically optimized plan in both cases. For case one, the MCO plan showed a 48% reduction of the 50% dose to OARs and a 16% reduction of the 1% dose to OARs. For case 2, the MCO plan showed a 72% reduction of the 50% dose to OARs and a 42% reduction of the 1% dose to OARs. The comparison of AWCM to WS showed that the AWCM method has a dosimetric advantage over WS for both patient cases.

CONCLUSION

We introduced an adaptive MCO algorithm for IMPT accelerated using GPUs. The algorithm is based on an adaptive method for generating Pareto plans in the objective space. We have shown that the algorithm can provide rapid and efficient mapping of the multicriteria Pareto surface with clinically deliverable plans.

Knowledge of endoscopic ultrasound-delivered fiducial composition and dimension necessary when planning proton beam radiotherapy

BACKGROUND AND STUDY AIMS

Little consideration has been given to selection of endoscopic ultrasound-guided fiducials for proton radiotherapy and the resulting perturbations in the therapy dose and pattern. Our aim was to assess the impact of perturbations caused by six fiducials of different composition and dimensions in a phantom gel model.

MATERIALS AND METHODS

 The phantom was submerged in a water bath and irradiated with a uniform 10 cm × 10 cm field of 119.7 MeV monoenergetic spot scanning protons delivered through a 45 mm range shifter. The proton "Bragg Peak" was evaluated.

RESULTS

 Dose perturbations manifesting as dose reductions up to 30 % were observed. A carbon composite (1 × 5 mm) and gold (0.4 × 10 mm) fiducial with backload potential rather than dedicated EUS pre-loaded gold fiducial needles had the best performance in terms of minimizing the dose perturbation.

CONCLUSIONS

 Our data demonstrate that a carbon composite fiducial has a less untoward effect on proton therapy dose distribution than dedicated EUS pre-loaded gold fiducial needles. Such information is important to consider when selecting fiducials specifically for proton therapy.

Clinical efficacy and safety of a highly conformal, supine, hybrid forward and inverse planned intensity modulated radiation therapy technique for craniospinal irradiation

PURPOSE

To demonstrate the clinical efficacy and safety of a highly conformal, supine, hybrid forward and inverse planned intensity modulated radiation therapy (IMRT) technique for photon craniospinal irradiation (CSI).

METHODS

Patients who received supine, hybrid IMRT CSI from 2009 to 2014 were included in this retrospective review. Clinical target volume (CTV) was defined as intracranial contents and thecal sac, including nerve roots. Dose was prescribed such that >99% of CTV received >99% of prescription and >95% of the planning target volume received >95% of prescription, with no attempt to include vertebral bodies. Lateral fields were utilized at the cranium and upper cervical spine. Spine fields were either single posterior or 2-3 obliques. Plans were generated with a hybrid of forward and inverse planned IMRT. Inferior borders of the cranium fields and superior border of the lower spine field were designed with 6-15 cm long, gradual dose gradients by sequential closing of multileaf collimator leaves using forward planned multiple static segment IMRT delivery. The sliding window upper spine IMRT field was created by the inverse planning system to match gradients of the brain and lower spine fields. The lower spine field gradient was similarly completed.

RESULTS

The cohort consisted of 34 patients. Median CSI dose was 36 Gy (range: 18-39.6 Gy). With a median follow up of 59.4 months, there were no isolated recurrences or spinal myelopathies at CTV margins or field gradients. Eleven patients had recurrence, all of which were intracranial.

CONCLUSIONS

Our hybrid forward and inverse planned IMRT supine CSI technique did not result in any isolated recurrences or myelopathies at CTV margins or field gradients. This suggests our target volumes and blended gradients are appropriate for highly conformal three-dimensional planning.

Biological Model for Predicting Toxicity in Head and Neck Cancer Patients Receiving Proton Therapy

PURPOSE

To use a linear energy transfer (LET) dependent formula for relative biological effectiveness (RBE) to generate a biological model that can be used to predict toxicity in patients treated with proton therapy for cancer of the head and neck.

PATIENTS AND METHODS

Patients treated with protons to a dose of 60 to 70 Gy (RBE = 1.1) for head and neck cancer were eligible to participate in this study. Treatment plans were developed using graphics processing unit Monte Carlo calculations. The equation, RBE = (1.1)[0.08(LETd)+0.88], was the biological model. The physical model assumes RBE = 1.1. Tumor volumes and organs at risk (OARs) were contoured, and isodose lines were created for 105%-120% of the prescribed dose. Dose to volume of OARs was calculated for both models for comparative purposes. Physician-reported toxicity was graded from 0 to 5 using the Common Terminology Criteria for Adverse Events, version 4.03. Patient-reported outcomes were obtained using the Promis10 and European Organisation for Research and Treatment of Cancer's QLQ-H&N35 instruments.

RESULTS

Eleven patients were included in this study. In each case the biological model revealed an increased dose to several OARs compared with the physical model. For selected OARs, the volume receiving >105% of the target dose was 2-fold to 15-fold greater in the biological model than the volume predicted by the physical model. Patients experienced toxicity that was consistent with the dose to OARs predicted by the biological model. Furthermore, 1 patient developed mucosal ulceration and another developed osteoradionecrosis at the location of a biological hot spot. In each case, the biological hot spot was located 2 mm inside the clinical target volume.

CONCLUSION

The results suggest that increases in dose predicted by the biological model are clinically relevant and that LET and RBE corrections and optimization should be a component of the treatment-planning process in proton therapy.

Implication of spot position error on plan quality and patient safety in pencil-beam-scanning proton therapy

PURPOSE

To quantitatively and systematically assess dosimetric effects induced by spot positioning error as a function of spot spacing (SS) on intensity-modulated proton therapy (IMPT) plan quality and to facilitate evaluation of safety tolerance limits on spot position.

METHODS

Spot position errors (PE) ranging from 1 to 2 mm were simulated. Simple plans were created on a water phantom, and IMPT plans were calculated on two pediatric patients with a brain tumor of 28 and 3 cc, respectively, using a commercial planning system. For the phantom, a uniform dose was delivered to targets located at different depths from 10 to 20 cm with various field sizes from 2(2) to 15(2) cm(2). Two nominal spot sizes, 4.0 and 6.6 mm of 1 σ in water at isocenter, were used for treatment planning. The SS ranged from 0.5 σ to 1.5 σ, which is 2-6 mm for the small spot size and 3.3-9.9 mm for the large spot size. Various perturbation scenarios of a single spot error and systematic and random multiple spot errors were studied. To quantify the dosimetric effects, percent dose error (PDE) depth profiles and the value of percent dose error at the maximum dose difference (PDE [ΔDmax]) were used for evaluation.

RESULTS

A pair of hot and cold spots was created per spot shift. PDE[ΔDmax] is found to be a complex function of PE, SS, spot size, depth, and global spot distribution that can be well defined in simple models. For volumetric targets, the PDE [ΔDmax] is not noticeably affected by the change of field size or target volume within the studied ranges. In general, reducing SS decreased the dose error. For the facility studied, given a single spot error with a PE of 1.2 mm and for both spot sizes, a SS of 1σ resulted in a 2% maximum dose error; a SS larger than 1.25 σ substantially increased the dose error and its sensitivity to PE. A similar trend was observed in multiple spot errors (both systematic and random errors). Systematic PE can lead to noticeable hot spots along the field edges, which may be near critical structures. However, random PE showed minimal dose error.

CONCLUSIONS

Dose error dependence for PE was quantitatively and systematically characterized and an analytic tool was built to simulate systematic and random errors for patient-specific IMPT. This information facilitates the determination of facility specific spot position error thresholds.

Design and characterization of an economical (192)Ir hemi-brain small animal irradiator

PURPOSE

To describe the design and dosimetric characterization of a simple and economical small animal irradiator.

MATERIALS AND METHODS

A high dose rate (HDR) (192)Ir brachytherapy source from a commercially available afterloader was used with a 1.3 cm thick tungsten collimator to provide sharp beam penumbra suitable for hemi-brain irradiation of mice. The unit was equipped with continuous gas anesthesia to allow robust animal immobilization. Dosimetric characterization of the device was performed with Gafchromic film measurements.

RESULTS

The tungsten collimator provided a sharp penumbra suitable for hemi-brain irradiation, and dose rates on the order of 200 cGy/minute were achieved. The sharpness of the penumbra attainable with this device compares favorably to those measured experimentally for 6 MV photons, and 6 and 20 MeV electron beams from a linear accelerator, and was comparable to those measured for a 300 kVp orthovoltage beam and a Monte Carlo simulated 90 MeV proton beam.

CONCLUSIONS

Due to its simplicity and low cost, the apparatus described is an attractive alternative for small animal irradiation experiments requiring steep dose gradients.

Feasibility of proton transmission-beam stereotactic ablative radiotherapy versus photon stereotactic ablative radiotherapy for lung tumors: a dosimetric and feasibility study

Stereotactic ablative radiotherapy is being increasingly adopted in the treatment of lung tumors. The use of proton beam therapy can further reduce dose to normal structures. However, uncertainty exists in proton-based treatment plans, including range uncertainties, large sensitivity to position uncertainty, and calculation of dose deposition in heterogeneous areas. This study investigated the feasibility of proton transmission beams, i.e. without the Bragg peak, to treat lung tumors with stereotactic ablative radiotherapy. We compared three representative treatment plans using proton transmission beams versus conformal static-gantry photon beams. It was found that proton treatment plans using transmission beams passing through the patient were feasible and demonstrated lower dose to normal structures and markedly reduced treatment times than photon plans. This is the first study to demonstrate the feasibility of proton-based stereotactic ablative radiotherapy planning for lung tumors using proton transmission beams alone. Further research using this novel approach for proton-based planning is warranted.