Proton therapy of prostate cancer by anterior-oblique beams: implications of setup and anatomy variations.

Proton therapy of prostate by anterior beams could offer an attractive option for treating patients with hip prosthesis and limiting the high-dose exposure to the rectum. We investigated the impact of setup and anatomy variations on the anterior-oblique (AO) proton plan dose, and strategies to manage these effects via range verification and adaptive delivery. Ten patients treated by bilateral (BL) passive-scattering proton therapy (79.2 Gy in 44 fractions) who underwent weekly verification CT scans were selected. Plans with AO beams were additionally created. To isolate the effect of daily variations, initial AO plans did not include range uncertainty margins. The use of fixed planning margins and adaptive range adjustments to manage these effects was investigated. For each case, the planned dose was recalculated on weekly CTs, and accumulated on the simulation CT using deformable registration to approximate the delivered dose. Planned and accumulated doses were compared for each scenario to quantify dose deviations induced by variations. The possibility of estimating the necessary range adjustments before each treatment was explored by simulating the procedure of a diode-based in vivo range verification technique, which would potentially be used clinically. The average planned rectum, penile bulb and femoral heads mean doses were smaller for initial AO compared to BL plans (by 8.3, 16.1 and 25.9 Gy, respectively). After considering interfractional variations in AO plans, the target coverage was substantially reduced. The maximum reduction of V 79.2/D 95/D mean/EUD for AO (without distal margins) (25.3%/10.7/1.6/4.9 Gy, respectively) was considerably larger than BL plans. The loss of coverage was mainly related to changes in water equivalent path length of the prostate after fiducial-based setup, caused by discrepancies in patient anterior surface and bony-anatomy alignment. Target coverage was recovered partially when using fixed planning margins, and fully when applying adaptive range adjustments. The accumulated organs-at-risk dose for AO beams after range adjustment demonstrated full sparing of femoral heads and superior sparing of penile bulb and rectum compared to the conventional BL cases. Our study indicates that using AO beams makes prostate treatment more susceptible to target underdose induced by interfractional variations. Adaptive range verification/adjustment may facilitate the use of anterior beam approaches, and ensure adequate target coverage in every fraction of the treatment.

Clinical commissioning of an in vivo range verification system for prostate cancer treatment with anterior and anterior oblique proton beams.

The purpose of this work is the clinical commissioning of a recently developed in vivo range verification system (IRVS) for treatment of prostate cancer by anterior and anterior oblique proton beams. The IRVS is designed to perform a complete workflow for pre-treatment range verification and adjustment. It contains specifically designed dosimetry and electronic hardware and a specific software for workflow control with database connection to the treatment and imaging systems. An essential part of the IRVS system is an array of Si-diode detectors, designed to be mounted to the endorectal water balloon routinely used for prostate immobilization. The diodes can measure dose rate as function of time from which the water equivalent path length (WEPL) and the dose received are extracted. The former is used for pre-treatment beam range verification and correction, if necessary, while the latter is to monitor the dose delivered to patient rectum during the treatment and serves as an additional verification. The entire IRVS workflow was tested for anterior and 30 degree inclined proton beam in both solid water and anthropomorphic pelvic phantoms, with the measured WEPL and rectal doses compared to the treatment plan. Gafchromic films were also used for measurement of the rectal dose and compared to IRVS results. The WEPL measurement accuracy was in the order of 1 mm and after beam range correction, the dose received by the rectal wall were 1.6% and 0.4% from treatment planning, respectively, for the anterior and anterior oblique field. We believe the implementation of IRVS would make the treatment of prostate with anterior proton beams more accurate and reliable.

Assessment of uncertainties in radiation-induced cancer risk predictions at clinically relevant doses.

PURPOSE: Theoretical dose-response models offer the possibility to assess second cancer induction risks after external beam therapy. The parameters used in these models are determined with limited data from epidemiological studies. Risk estimations are thus associated with considerable uncertainties. This study aims at illustrating uncertainties when predicting the risk for organ-specific second cancers in the primary radiation field illustrated by choosing selected treatment plans for brain cancer patients. METHODS: A widely used risk model was considered in this study. The uncertainties of the model parameters were estimated with reported data of second cancer incidences for various organs. Standard error propagation was then subsequently applied to assess the uncertainty in the risk model. Next, second cancer risks of five pediatric patients treated for cancer in the head and neck regions were calculated. For each case, treatment plans for proton and photon therapy were designed to estimate the uncertainties (a) in the lifetime attributable risk (LAR) for a given treatment modality and (b) when comparing risks of two different treatment modalities. RESULTS: Uncertainties in excess of 100% of the risk were found for almost all organs considered. When applied to treatment plans, the calculated LAR values have uncertainties of the same magnitude. A comparison between cancer risks of different treatment modalities, however, does allow statistically significant conclusions. In the studied cases, the patient averaged LAR ratio of proton and photon treatments was 0.35, 0.56, and 0.59 for brain carcinoma, brain sarcoma, and bone sarcoma, respectively. Their corresponding uncertainties were estimated to be potentially below 5%, depending on uncertainties in dosimetry. CONCLUSIONS: The uncertainty in the dose-response curve in cancer risk models makes it currently impractical to predict the risk for an individual external beam treatment. On the other hand, the ratio of absolute risks between two modalities is less sensitive to the uncertainties in the risk model and can provide statistically significant estimates.

Validation of a deformable image registration technique for cone beam CT-based dose verification.

PURPOSE: As radiation therapy evolves toward more adaptive techniques, image guidance plays an increasingly important role, not only in patient setup but also in monitoring the delivered dose and adapting the treatment to patient changes. This study aimed to validate a method for evaluation of delivered intensity modulated radiotherapy (IMRT) dose based on multimodal deformable image registration (dir) for prostate treatments. METHODS: A pelvic phantom was scanned with CT and cone-beam computed tomography (CBCT). Both images were digitally deformed using two realistic patient-based deformation fields. The original CT was then registered to the deformed CBCT resulting in a secondary deformed CT. The registration quality was assessed as the ability of the dir method to recover the artificially induced deformations. The primary and secondary deformed CT images as well as vector fields were compared to evaluate the efficacy of the registration method and it's suitability to be used for dose calculation. plastimatch, a free and open source software was used for deformable image registration. A B-spline algorithm with optimized parameters was used to achieve the best registration quality. Geometric image evaluation was performed through voxel-based Hounsfield unit (HU) and vector field comparison. For dosimetric evaluation, IMRT treatment plans were created and optimized on the original CT image and recomputed on the two warped images to be compared. The dose volume histograms were compared for the warped structures that were identical in both warped images. This procedure was repeated for the phantom with full, half full, and empty bladder. RESULTS: The results indicated mean HU differences of up to 120 between registered and ground-truth deformed CT images. However, when the CBCT intensities were calibrated using a region of interest (ROI)-based calibration curve, these differences were reduced by up to 60%. Similarly, the mean differences in average vector field lengths decreased from 10.1 to 2.5 mm when CBCT was calibrated prior to registration. The results showed no dependence on the level of bladder filling. In comparison with the dose calculated on the primary deformed CT, differences in mean dose averaged over all organs were 0.2% and 3.9% for dose calculated on the secondary deformed CT with and without CBCT calibration, respectively, and 0.5% for dose calculated directly on the calibrated CBCT, for the full-bladder scenario. Gamma analysis for the distance to agreement of 2 mm and 2% of prescribed dose indicated a pass rate of 100% for both cases involving calibrated CBCT and on average 86% without CBCT calibration. CONCLUSIONS: Using deformable registration on the planning CT images to evaluate the IMRT dose based on daily CBCTs was found feasible. The proposed method will provide an accurate dose distribution using planning CT and pretreatment CBCT data, avoiding the additional uncertainties introduced by CBCT inhomogeneity and artifacts. This is a necessary initial step toward future image-guided adaptive radiotherapy of the prostate.

Dosimetric feasibility of real-time MRI-guided proton therapy.

PURPOSE: Magnetic resonance imaging (MRI) is a prime candidate for image-guided radiotherapy. This study was designed to assess the feasibility of real-time MRI-guided proton therapy by quantifying the dosimetric effects induced by the magnetic field in patients' plans and identifying the associated clinical consequences. METHODS: Monte Carlo dose calculation was performed for nine patients of various treatment sites (lung, liver, prostate, brain, skull-base, and spine) and tissue homogeneities, in the presence of 0.5 and 1.5 T magnetic fields. Dose volume histogram (DVH) parameters such as D95, D5, and V20 as well as equivalent uniform dose were compared for the target and organs at risk, before and after applying the magnetic field. The authors further assessed whether the plans affected by clinically relevant dose distortions could be corrected independent of the planning system. RESULTS: By comparing the resulting dose distributions and analyzing the respective DVHs, it was determined that despite the observed lateral beam deflection, for magnetic fields of up to 0.5 T, neither was the target coverage jeopardized nor was the dose to the nearby organs increased in all cases except for prostate. However, for a 1.5 T magnetic field, the dose distortions were more pronounced and of clinical concern in all cases except for spine. In such circumstances, the target was severely underdosed, as indicated by a decrease in D95 of up to 41% of the prescribed dose compared to the nominal situation (no magnetic field). Sites such as liver and spine were less affected due to higher tissue homogeneity, typically smaller beam range, and the choice of beam directions. Simulations revealed that small modifications to certain plan parameters such as beam isocenter (up to 19 mm) and gantry angle (up to 10°) are sufficient to compensate for the magnetic field-induced dose disturbances. The authors' observations indicate that the degree of required corrections strongly depends on the beam range and direction relative to the magnetic field. This method was also applicable to more heterogeneous scenarios such as skull-base tumors. CONCLUSIONS: This study confirmed the dosimetric feasibility of real-time MRI-guided proton therapy and delivering a clinically acceptable dose to patients with various tumor locations within magnetic fields of up to 1.5 T. This work could serve as a guide and encouragement for further efforts toward clinical implementation of hybrid MRI-proton gantry systems.

SU-E-T-258: Assessment of Radiation Induced Second Cancer Risks in Proton Therapy and IMRT for Organs inside the Main Radiation Field.

PURPOSE: Radiation therapy can potentially cause a second malignancy. There is clinical evidence that those occur typically within the beam path in the medium/high dose region. The purpose of this study was to assess the risk for developing a radiation induced tumor within the treated volume and to compare this risk for proton therapy and IMRT. METHODS: Fully contoured age and gender specific whole body phantoms (4-year and 14-year old) were uploaded into a treatment planning system and typical tumor volumes were contoured based on patients treated for optic glioma and vertebral body Ewing's sarcoma. Lifetime attributable risks (LARs) for developing a second malignancy were calculated using a risk model incorporating factors for cell kill, mutations, repopulation, and inhomogeneous organ doses. RESULTS: For standard fractionation schemes, the LAR for developing a second malignancy from radiation therapy alone were found to be up to 2.7% for a 4-year old optic glioma patient treated with IMRT considering a soft tissue carcinoma risk model only. Sarcoma risks were found to be below 1% in all cases. For the 14-year old, risks were found to be about a factor of 2 lower. For the Ewing's sarcoma cases the risks based on the sarcoma model were typically higher than the carcinoma risks, i.e. up to 1.3% LAR for soft tissue sarcoma. Generally, the risk from proton therapy turned out to be lower by a factor of 2 to 10. However, comparison of a 3-field and 4-field proton plan shows that the distribution of the dose, i.e., the particular treatment plan, plays a role as well. CONCLUSIONS: In general, proton therapy can significantly reduce the risk for developing an in-field second malignancy. Risk analysis based on our formalism could be applied within treatment planning programs to guide treatment plans for pediatric patients. Federal Share of program income earned by Massachusetts General Hospital on NIH/NCI C06 CA059267.

SU-E-T-470: Comparison of Proton Treatment Planning and Monte Carlo Calculation Using TOPAS for Liver Cancer.

PURPOSE: To compare Monte Carlo (MC) calculated and planned dose distributions (pencil beam algorithm) for patients with liver cancer treated with proton radiation therapy. METHODS: Six patients with unresectable Hepatocellular carcinoma were chosen from the institutional protocol list. We applied the newly developed TOPAS (Tool for Particle Simulation) Monte Carlo (MC) tool and an in-house (mcauto) program, which connects the planning system with the MC. Two beams, typically right lateral (RL) and anterior-posterior (AP), were simulated for each patient with a total prescribed dose of 58 Gy. The calculated absolute dose was determined by separately simulating an SOBP dose in a water phantom for normalization to the prescription dose. The difference between MC and planned dose were calculated and Dose Volume Histograms (DVHs) for the critical organs with non-negligible dose (whole liver, heart, small and large bowel and chest wall) were analyzed. RESULTS: The resulting dose distributions were in quite good agreement. The main discrepancy in all cases was observed in the lateral penumbrae. These discrepancies can mainly result from the range compensator gradient and tissue composition. The Dose Volume Histograms (DVHs) also presented good agreement between doses for the CTV as well as all the OARs. The difference in D95 ranged from 0.7-1.5 Gy that is translated to 1.3-2.5% of the prescribed dose. CONCLUSIONS: TOPAS Monte Carlo tool presented an efficient and accurate method for dose calculation in liver and to validate clinical treatment planning. Discrepancies with doses calculated using the pencil beam algorithm were seen but were generally quite small.

Monte Carlo patient study on the comparison of prompt gamma and PET imaging for range verification in proton therapy.

The purpose of this work was to compare the clinical adaptation of prompt gamma (PG) imaging and positron emission tomography (PET) as independent tools for non-invasive proton beam range verification and treatment validation. The PG range correlation and its differences with PET have been modeled for the first time in a highly heterogeneous tissue environment, using different field sizes and configurations. Four patients with different tumor locations (head and neck, prostate, spine and abdomen) were chosen to compare the site-specific behaviors of the PG and PET images, using both passive scattered and pencil beam fields. Accurate reconstruction of dose, PG and PET distributions was achieved by using the planning computed tomography (CT) image in a validated GEANT4-based Monte Carlo code capable of modeling the treatment nozzle and patient anatomy in detail. The physical and biological washout phenomenon and decay half-lives for PET activity for the most abundant isotopes such as (11)C, (15)O, (13)N, (30)P and (38)K were taken into account in the data analysis. The attenuation of the gamma signal after traversing the patient geometry and respective detection efficiencies were estimated for both methods to ensure proper comparison. The projected dose, PG and PET profiles along many lines in the beam direction were analyzed to investigate the correlation consistency across the beam width. For all subjects, the PG method showed on average approximately 10 times higher gamma production rates than the PET method before, and 60 to 80 times higher production after including the washout correction and acquisition time delay. This rate strongly depended on tissue density and elemental composition. For broad passive scattered fields, it was demonstrated that large differences exist between PG and PET signal falloff positions and the correlation with the dose distribution for different lines in the beam direction. These variations also depended on the treatment site and the particular subject. Thus, similar to PET, direct range verification with PG in passive scattering is not easily viable. However, upon development of an optimized 3D PG detector, indirect range verification by comparing measured and simulated PG distributions (currently being explored for the PET method) would be more beneficial because it can avoid the inherent biological challenges of the PET imaging. The improved correlation of PG and PET with dose when using pencil beams was evident. PG imaging was found to be potentially advantageous especially for small tumors in the presence of high tissue heterogeneities. Including the effects of detector acceptance and efficiency may hold PET superior in terms of the amplitude of the detected signal (depending on the future development of PG detection technology), but the ability to perform online measurements and avoid signal disintegration (due to washout) with PG are important factors that can outweigh the benefits of higher detection sensitivity.

Probing the high momentum component of the deuteron at high Q2.

The (2)H(e,e'p)n cross section at a momentum transfer of 3.5 (GeV/c)(2) was measured over a kinematical range that made it possible to study this reaction for a set of fixed missing momenta as a function of the neutron recoil angle θ(nq) and to extract missing momentum distributions for fixed values of θ(nq) up to 0.55 GeV/c. In the region of 35°≤θ(nq)≤45° recent calculations, which predict that final-state interactions are small, agree reasonably well with the experimental data. Therefore, these experimental reduced cross sections provide direct access to the high momentum component of the deuteron momentum distribution in exclusive deuteron electrodisintegration.