RANGE ADAPTIVE PROTON THERAPY FOR PROSTATE CANCER

Purpose: The rapid distal falloff of a proton beam allows for sparing of normal tissues distal to the target. However proton beams that aim directly towards critical structures are avoided due to concerns of range uncertainties, such as CT number conversion and anatomy variations. We propose to eliminate range uncertainty and enable prostate treatment with a single anterior beam by detecting the proton’s range at the prostate-rectal interface and adaptively adjusting the range in vivo and in real-time. Materials and Methods: A prototype device, consisting of an endorectal liquid scintillation detector and dual-inverted Lucite wedges for range compensation, was designed to test the feasibility and accuracy of the technique. Liquid scintillation filled volume was fitted with optical fiber and placed inside the rectum of an anthropomorphic pelvic phantom. Photodiode-generated current signal was generated as a function of proton beam distal depth, and the spatial resolution of this technique was calculated by relating the variance in detecting proton spills to its maximum penetration depth. The relative water-equivalent thickness of the wedges was measured in a water phantom and prospectively tested to determine the accuracy of range corrections. Treatment simulation studies were performed to test the potential dosimetric benefit in sparing the rectum. Results: The spatial resolution of the detector in phantom measurement was 0.5 mm. The precision of the range correction was 0.04 mm. The residual margin to ensure CTV coverage was 1.1 mm. The composite distal margin for 95% treatment confidence was 2.4 mm. Planning studies based on a previously estimated 2mm margin (90% treatment confidence) for 27 patients showed a rectal sparing up to 51% at 70 Gy and 57% at 40 Gy relative to IMRT and bilateral proton treatment. Conclusion: We demonstrated the feasibility of our design. Use of this technique allows for proton treatment using a single anterior beam, significantly reducing the rectal dose.

Implementation of an Anthropomorphic Phantom for the Evaluation of Proton Therapy Treatment Procedures

With an increasing number of institutions offering proton therapy, the number of multi-institutional clinical trials involving proton therapy will also increase in the coming years. The Radiological Physics Center monitors sites involved in clinical trials through the use of site visits and remote auditing with thermoluminescent dosimeters (TLD) and mailable anthropomorphic phantoms. Currently, there are no heterogeneous phantoms that have been commissioned to evaluate proton therapy. It was hypothesized that an anthropomorphic pelvis phantom can be designed to audit treatment procedures (patient simulation, treatment planning and treatment delivery) at proton facilities to confirm agreement between the measured dose and calculated dose within 5%/3mm with a reproducibility of 3%. A pelvis phantom originally designed for use with photon treatments was retrofitted for use in proton therapy. The relative stopping power (SP) of each phantom material was measured. Hounsfield Units (HU) for each phantom material were measured with a CT scanner and compared to the relative stopping power calibration curve. The tissue equivalency for each material was calculated. Two proton treatment plans were created; one which did not correct for material SP differences (Plan 1) and one plan which did correct for SP differences (Plan 2). Film and TLD were loaded into the phantom and the phantom was irradiated 3 times per plan. The measured values were compared to the HU-SP calibration curve and it was found that the stopping powers for the materials could be underestimated by 5-10%. Plan 1 passed the criteria for the TLD and film margins with reproducibility under 3% between the 3 trials. Plan 2 failed because the right-left film dose profile average displacement was -9.0 mm on the left side and 6.0 mm on the right side. Plan 2 was intended to improve the agreements and instead introduced large displacements along the path of the beam. Plan 2 more closely represented the actual phantom composition with corrected stopping powers and should have shown an agreement between the measured and calculated dose within 5%/3mm. The hypothesis was rejected and the pelvis phantom was found to be not suitable to evaluate proton therapy treatment procedures.

MODIFICATION AND IMPLEMENTATION OF THE RPC HETEROGENEOUS THORAX PHANTOM FOR VERIFICATION OF PROTON THERAPY TREATMENT PROCEDURES

The Radiological Physics Center (RPC) provides heterogeneous phantoms that are used to evaluate radiation treatment procedures as part of a comprehensive quality assurance program for institutions participating in clinical trials. It was hypothesized that the existing RPC heterogeneous thorax phantom can be modified to assess lung tumor proton beam therapy procedures involving patient simulation, treatment planning, and treatment delivery, and could confirm agreement between the measured dose and calculated dose within 5%/3mm with a reproducibility of 5%. The Hounsfield Units (HU) for lung equivalent materials (balsa wood and cork) was measured using a CT scanner. The relative linear stopping power (RLSP) of these materials was measured. The linear energy transfer (LET) of Gafchromic EBT2 film was analyzed utilizing parallel and perpendicular orientations in a water tank and compared to ion chamber readings. Both parallel and perpendicular orientations displayed a quenching effect underperforming the ion chamber, with the parallel orientation showing an average 31 % difference and the perpendicular showing an average of 15% difference. Two treatment plans were created that delivered the prescribed dose to the target volume, while achieving low entrance doses. Both treatment plans were designed using smeared compensators and expanded apertures, as would be utilized for a patient in the clinic. Plan 1a contained two beams that were set to orthogonal angles and a zero degree couch kick. Plan 1b utilized two beams set to 10 and 80 degrees with a 15 degree couch kick. EBT2 film and TLD were inserted and the phantom was irradiated 3 times for each plan. Both plans passed the criteria for the TLD measurements where the TLD values were within 7% of the dose calculated by Eclipse. Utilizing the 5%/3mm criteria, the 3 trial average of overall pass rate was 71% for Plan 1a. The 3 trial average for the overall pass rate was 76% for Plan 1b. The trials were then analyzed using RPC conventional lung treatment guidelines set forth by the RTOG: 5%/5mm, and an overall pass rate of 85%. Utilizing these criteria, only Plan 1b passed for all 3 trials, with an average overall pass rate of 89%.

DUAL ENERGY COMPUTED TOMOGRAPHY FOR PROTON THERAPY TREATMENT PLANNING

Proton radiation therapy is gaining popularity because of the unique characteristics of its dose distribution, e.g., high dose-gradient at the distal end of the percentage-depth-dose curve (known as the Bragg peak). The high dose-gradient offers the possibility of delivering high dose to the target while still sparing critical organs distal to the target. However, the high dose-gradient is a double-edged sword: a small shift of the highly conformal high-dose area can cause the target to be substantially under-dosed or the critical organs to be substantially over-dosed. Because of that, large margins are required in treatment planning to ensure adequate dose coverage of the target, which prevents us from realizing the full potential of proton beams. Therefore, it is critical to reduce uncertainties in the proton radiation therapy. One major uncertainty in a proton treatment is the range uncertainty related to the estimation of proton stopping power ratio (SPR) distribution inside a patient. The SPR distribution inside a patient is required to account for tissue heterogeneities when calculating dose distribution inside the patient. In current clinical practice, the SPR distribution inside a patient is estimated from the patient’s treatment planning computed tomography (CT) images based on the CT number-to-SPR calibration curve. The SPR derived from a single CT number carries large uncertainties in the presence of human tissue composition variations, which is the major drawback of the current SPR estimation method. We propose to solve this problem by using dual energy CT (DECT) and hypothesize that the range uncertainty can be reduced by a factor of two from currently used value of 3.5%. A MATLAB program was developed to calculate the electron density ratio (EDR) and effective atomic number (EAN) from two CT measurements of the same object. An empirical relationship was discovered between mean excitation energies and EANs existing in human body tissues. With the MATLAB program and the empirical relationship, a DECT-based method was successfully developed to derive SPRs for human body tissues (the DECT method). The DECT method is more robust against the uncertainties in human tissues compositions than the current single-CT-based method, because the DECT method incorporated both density and elemental composition information in the SPR estimation. Furthermore, we studied practical limitations of the DECT method. We found that the accuracy of the DECT method using conventional kV-kV x-ray pair is susceptible to CT number variations, which compromises the theoretical advantage of the DECT method. Our solution to this problem is to use a different x-ray pair for the DECT. The accuracy of the DECT method using different combinations of x-ray energies, i.e., the kV-kV, kV-MV and MV-MV pair, was compared using the measured imaging uncertainties for each case. The kV-MV DECT was found to be the most robust against CT number variations. In addition, we studied how uncertainties propagate through the DECT calculation, and found general principles of selecting x-ray pairs for the DECT method to minimize its sensitivity to CT number variations. The uncertainties in SPRs estimated using the kV-MV DECT were analyzed further and compared to those using the stoichiometric method. The uncertainties in SPR estimation can be divided into five categories according to their origins: the inherent uncertainty, the DECT modeling uncertainty, the CT imaging uncertainty, the uncertainty in the mean excitation energy, and SPR variation with proton energy. Additionally, human body tissues were divided into three tissue groups – low density (lung) tissues, soft tissues and bone tissues. The uncertainties were estimated separately because their uncertainties were different under each condition. An estimate of the composite range uncertainty (2s) was determined for three tumor sites – prostate, lung, and head-and-neck, by combining the uncertainty estimates of all three tissue groups, weighted by their proportions along typical beam path for each treatment site. In conclusion, the DECT method holds theoretical advantages in estimating SPRs for human tissues over the current single-CT-based method. Using existing imaging techniques, the kV-MV DECT approach was capable of reducing the range uncertainty from the currently used value of 3.5% to 1.9%-2.3%, but it is short to reach our original goal of reducing the range uncertainty by a factor of two. The dominant source of uncertainties in the kV-MV DECT was the uncertainties in CT imaging, especially in MV CT imaging. Further reduction in beam hardening effect, the impact of scatter, out-of-field object etc. would reduce the Hounsfeld Unit variations in CT imaging. The kV-MV DECT still has the potential to reduce the range uncertainty further.

The Development and Implementation of an Anthropomorphic Head Phantom for the Assessment of Proton Therapy Treatment Procedures

Proton therapy has become an increasingly more common method of radiation therapy, with the dose sparing to distal tissue making it an appealing option, particularly for treatment of brain tumors. This study sought to develop a head phantom for the Radiological Physics Center (RPC), the first to be used for credentialing of institutions wishing to participate in clinical trials involving brain tumor treatment of proton therapy. It was hypothesized that a head phantom could be created for the evaluation of proton therapy treatment procedures (treatment simulation, planning, and delivery) to assure agreement between the measured dose and calculated dose within ±5%/3mm with a reproducibility of ±3%. The relative stopping power (RSP) and Hounsfield Units (HU) were measured for potential phantom materials and a human skull was cast in tissue-equivalent Alderson material (RLSP 1.00, HU 16) with anatomical airways and a cylindrical hole for imaging and dosimetry inserts drilled into the phantom material. Two treatment plans, proton passive scattering and proton spot scanning, were created. Thermoluminescent dosimeters (TLDs) and film were loaded into the phantom dosimetry insert. Each treatment plan was delivered three separate times. Each treatment plan passed our 5%/3mm criteria, with a reproducibility of ±3%. The hypothesis was accepted and the phantom was found to be suitable for remote audits of proton therapy treatment facilities.

The effect of surgical titanium rods on proton therapy delivered for cervical bone tumors: experimental validation using an anthropomorphic phantom.

To investigate the effect of metal implants in proton radiotherapy, dose distributions of different, clinically relevant treatment plans have been measured in an anthropomorphic phantom and compared to treatment planning predictions. The anthropomorphic phantom, which is sliced into four segments in the cranio-caudal direction, is composed of tissue equivalent materials and contains a titanium implant in a vertebral body in the cervical region. GafChromic® films were laid between the different segments to measure the 2D delivered dose. Three different four-field plans have then been applied: a Single-Field-Uniform-Dose (SFUD) plan, both with and without artifact correction implemented, and an Intensity-Modulated-Proton-Therapy (IMPT) plan with the artifacts corrected. For corrections, the artifacts were manually outlined and the Hounsfield Units manually set to an average value for soft tissue. Results show a surprisingly good agreement between prescribed and delivered dose distributions when artifacts have been corrected, with > 97% and 98% of points fulfilling the gamma criterion of 3%/3 mm for both SFUD and the IMPT plans, respectively. In contrast, without artifact corrections, up to 18% of measured points fail the gamma criterion of 3%/3 mm for the SFUD plan. These measurements indicate that correcting manually for the reconstruction artifacts resulting from metal implants substantially improves the accuracy of the calculated dose distribution.

Prophylactic use of bevacizumab to avoid anterior segment neovascularization following proton therapy for uveal melanoma.

PURPOSE To investigate whether the prophylactic use of bevacizumab reduces the rate of rubeosis after proton therapy for uveal melanoma and improves the possibility to treat ischemic, reapplicated retina with laser photocoagulation. DESIGN Comparative retrospective case series. METHODS Uveal melanoma patients with ischemic retinal detachment and treated with proton therapy were included in this institutional study. Twenty-four eyes received prophylactic intravitreal bevacizumab injections and were compared with a control group of 44 eyes without bevacizumab treatment. Bevacizumab injections were performed at the time of tantalum clip insertion and were repeated every 2 months during 6 months, and every 3 months thereafter. Ultra-widefield angiography allowed determination of the extent of retinal ischemia, which was treated with laser photocoagulation after retinal reapplication. Main outcome measures were the time to rubeosis, the time to retinal reattachment, and the time to laser photocoagulation of ischemic retina. RESULTS Baseline characteristics were balanced between the groups, except for thicker tumors and larger retinal detachments in the bevacizumab group, potentially to the disadvantage of the study group. Nevertheless, bevacizumab prophylaxis significantly reduced the rate of iris rubeosis from 36% to 4% (log-rank test P = .02) and tended to shorten the time to retinal reapplication until laser photocoagulation of the nonperfusion areas could be performed. CONCLUSIONS Prophylactic intravitreal bevacizumab in patients treated with proton therapy for uveal melanoma with ischemic retinal detachment prevented anterior segment neovascularization, until laser photocoagulation to the reapplied retina could be performed.

Long-term outcomes and prognostic factors of skull-base chondrosarcoma patients treated with pencil-beam scanning proton therapy at the Paul Scherrer Institute

BACKGROUND Skull-base chondrosarcoma (ChSa) is a rare disease, and the prognostication of this disease entity is ill defined. METHODS We assessed the long-term local control (LC) results, overall survival (OS), and prognostic factors of skull-base ChSa patients treated with pencil beam scanning proton therapy (PBS PT). Seventy-seven (male, 35; 46%) patients with histologically confirmed ChSa were treated at the Paul Scherrer Institute. Median age was 38.9 years (range, 10.2-70.0y). Median delivered dose was 70.0 GyRBE (range, 64.0-76.0 GyRBE). LC, OS, and toxicity-free survival (TFS) rates were calculated using the Kaplan Meier method. RESULTS After a mean follow-up of 69.2 months (range, 4.6-190.8 mo), 6 local (7.8%) failures were observed, 2 of which were late failures. Five (6.5%) patients died. The actuarial 8-year LC and OS were 89.7% and 93.5%, respectively. Tumor volume > 25 cm(3) (P = .02), brainstem/optic apparatus compression at the time of PT (P = .04) and age >30 years (P = .08) were associated with lower rates of LC. High-grade (≥3) radiation-induced toxicity was observed in 6 (7.8%) patients. The 8-year high-grade TFS was 90.8%. A higher rate of high-grade toxicity was observed for older patients (P = .073), those with larger tumor volume (P = .069), and those treated with 5 weekly fractions (P = .069). CONCLUSIONS This is the largest PT series reporting the outcome of patients with low-grade ChSa of the skull base treated with PBS only. Our data indicate that protons are both safe and effective. Tumor volume, brainstem/optic apparatus compression, and age were prognosticators of local failures.

Pencil Beam Scanning Proton Therapy for Pediatric Parameningeal Rhabdomyosarcomas: Clinical Outcome of Patients Treated at the Paul Scherrer Institute

BACKGROUND Parameningeal rhabdomyosarcomas (PM-RMSs) represent approximately 25% of all rhabdomyosarcoma (RMS) cases. These tumors are associated with early recurrence and poor prognosis. This study assessed the clinical outcome and late toxicity of pencil beam scanning (PBS) proton therapy (PT) in the treatment of children with PM-RMS. PROCEDURES Thirty-nine children with PM-RMS received neoadjuvant chemotherapy followed by PBS-PT at the Paul Scherrer Institute, with concomitant chemotherapy. The median age was 5.8 years (range, 1.2-16.1). Due to young age, 25 patients (64%) required general anesthesia during PT. The median time from the start of chemotherapy to PT was 13 weeks (range, 3-23 weeks). Median prescription dose was 54 Gy (relative biologic effectiveness, RBE). RESULTS With a mean follow-up of 41 months (range, 9-106 months), 10 patients failed. The actuarial 5-year progression-free survival (PFS) was 72% (95% CI, 67-94%) and the 5-year overall survival was 73% (95% CI, 69-96%). On univariate analysis, a delay in the initiation of PT (>13 weeks) was a significant detrimental factor for PFS. Three (8%) patients presented with grade 3 radiation-induced toxicity. The estimated actuarial 5-year toxicity ≥grade 3 free survival was 95% (95% CI, 94-96%). CONCLUSIONS Our data contribute to the growing body of evidence demonstrating the safety and effectiveness of PT for pediatric patients with PM-RMS. These preliminary results are encouraging and in line with other combined proton-photon and photons series; observed toxicity was acceptable.

Treatment planning aspects and Monte Carlo methods in proton therapy

Over the last years, the interest in proton radiotherapy is rapidly increasing. Protons provide superior physical properties compared with conventional radiotherapy using photons. These properties result in depth dose curves with a large dose peak at the end of the proton track and the finite proton range allows sparing the distally located healthy tissue. These properties offer an increased flexibility in proton radiotherapy, but also increase the demand in accurate dose estimations. To carry out accurate dose calculations, first an accurate and detailed characterization of the physical proton beam exiting the treatment head is necessary for both currently available delivery techniques: scattered and scanned proton beams. Since Monte Carlo (MC) methods follow the particle track simulating the interactions from first principles, this technique is perfectly suited to accurately model the treatment head. Nevertheless, careful validation of these MC models is necessary. While for the dose estimation pencil beam algorithms provide the advantage of fast computations, they are limited in accuracy. In contrast, MC dose calculation algorithms overcome these limitations and due to recent improvements in efficiency, these algorithms are expected to improve the accuracy of the calculated dose distributions and to be introduced in clinical routine in the near future.

Robustness of the Voluntary Breath-Hold Approach for the Treatment of Peripheral Lung Tumors Using Hypofractionated Pencil Beam Scanning Proton Therapy.

PURPOSE The safe clinical implementation of pencil beam scanning (PBS) proton therapy for lung tumors is complicated by the delivery uncertainties caused by breathing motion. The purpose of this feasibility study was to investigate whether a voluntary breath-hold technique could limit the delivery uncertainties resulting from interfractional motion. METHODS AND MATERIALS Data from 15 patients with peripheral lung tumors previously treated with stereotactic radiation therapy were included in this study. The patients had 1 computed tomographic (CT) scan in voluntary breath-hold acquired before treatment and 3 scans during the treatment course. PBS proton treatment plans with 2 fields (2F) and 3 fields (3F), respectively, were calculated based on the planning CT scan and subsequently recalculated on the 3 repeated CT scans. Recalculated plans were considered robust if the V95% (volume receiving ≥95% of the prescribed dose) of the gross target volume (GTV) was within 5% of what was expected from the planning CT data throughout the simulated treatment. RESULTS A total of 14/15 simulated treatments for both 2F and 3F met the robustness criteria. Reduced V95% was associated with baseline shifts (2F, P=.056; 3F, P=.008) and tumor size (2F, P=.025; 3F, P=.025). Smaller tumors with large baseline shifts were also at risk for reduced V95% (interaction term baseline/size: 2F, P=.005; 3F, P=.002). CONCLUSIONS The breath-hold approach is a realistic clinical option for treating lung tumors with PBS proton therapy. Potential risk factors for reduced V95% are small targets in combination with large baseline shifts. On the basis of these results, the baseline shift of the tumor should be monitored (eg, through image guided therapy), and appropriate measures should be taken accordingly. The intrafractional motion needs to be investigated to confirm that the breath-hold approach is robust.

Dose-painting intensity-modulated proton therapy for intermediate- and high-risk meningioma.

BACKGROUND Newly diagnosed WHO grade II-III or any WHO grade recurrent meningioma exhibit an aggressive behavior and thus are considered as high- or intermediate risk tumors. Given the unsatisfactory rates of disease control and survival after primary or adjuvant radiation therapy, optimization of treatment strategies is needed. We investigated the potential of dose-painting intensity-modulated proton beam-therapy (IMPT) for intermediate- and high-risk meningioma. MATERIAL AND METHODS Imaging data from five patients undergoing proton beam-therapy were used. The dose-painting target was defined using [68]Ga-[1,4,7,10-tetraazacyclododecane tetraacetic acid]- d-Phe(1),Tyr(3)-octreotate ([68]Ga-DOTATATE)-positron emission tomography (PET) in target delineation. IMPT and photon intensity-modulated radiation therapy (IMRT) treatment plans were generated for each patient using an in-house developed treatment planning system (TPS) supporting spot-scanning technology and a commercial TPS, respectively. Doses of 66 Gy (2.2 Gy/fraction) and 54 Gy (1.8 Gy/fraction) were prescribed to the PET-based planning target volume (PTVPET) and the union of PET- and anatomical imaging-based PTV, respectively, in 30 fractions, using simultaneous integrated boost. RESULTS Dose coverage of the PTVsPET was equally good or slightly better in IMPT plans: dose inhomogeneity was 10 ± 3% in the IMPT plans vs. 13 ± 1% in the IMRT plans (p = 0.33). The brain Dmean and brainstem D50 were small in the IMPT plans: 26.5 ± 1.5 Gy(RBE) and 0.002 ± 0.0 Gy(RBE), respectively, vs. 29.5 ± 1.5 Gy (p = 0.001) and 7.5 ± 11.1 Gy (p = 0.02) for the IMRT plans, respectively. The doses delivered to the optic structures were also decreased with IMPT. CONCLUSIONS Dose-painting IMPT is technically feasible using currently available planning tools and resulted in dose conformity of the dose-painted target comparable to IMRT with a significant reduction of radiation dose delivered to the brain, brainstem and optic apparatus. Dose escalation with IMPT may improve tumor control and decrease radiation-induced toxicity.

Monte Carlo and analytical dose calculations for ocular proton therapy

Uveal melanoma is a rare but life-threatening form of ocular cancer. Contemporary treatment techniques include proton therapy, which enables conservation of the eye and its useful vision. Dose to the proximal structures is widely believed to play a role in treatment side effects, therefore, reliable dose estimates are required for properly evaluating the therapeutic value and complication risk of treatment plans. Unfortunately, current simplistic dose calculation algorithms can result in errors of up to 30% in the proximal region. In addition, they lack predictive methods for absolute dose per monitor unit (D/MU) values. ^ To facilitate more accurate dose predictions, a Monte Carlo model of an ocular proton nozzle was created and benchmarked against measured dose profiles to within ±3% or ±0.5 mm and D/MU values to within ±3%. The benchmarked Monte Carlo model was used to develop and validate a new broad beam dose algorithm that included the influence of edgescattered protons on the cross-field intensity profile, the effect of energy straggling in the distal portion of poly-energetic beams, and the proton fluence loss as a function of residual range. Generally, the analytical algorithm predicted relative dose distributions that were within ±3% or ±0.5 mm and absolute D/MU values that were within ±3% of Monte Carlo calculations. Slightly larger dose differences were observed at depths less than 7 mm, an effect attributed to the dose contributions of edge-scattered protons. Additional comparisons of Monte Carlo and broad beam dose predictions were made in a detailed eye model developed in this work, with generally similar findings. ^ Monte Carlo was shown to be an excellent predictor of the measured dose profiles and D/MU values and a valuable tool for developing and validating a broad beam dose algorithm for ocular proton therapy. The more detailed physics modeling by the Monte Carlo and broad beam dose algorithms represent an improvement in the accuracy of relative dose predictions over current techniques, and they provide absolute dose predictions. It is anticipated these improvements can be used to develop treatment strategies that reduce the incidence or severity of treatment complications by sparing normal tissue. ^