Intensity Modulated Proton Therapy as a Boost Treatment after Prostate Seed Implant: A Treatment Planning Study ()
1. Introduction
Combined radiotherapy of prostate seed implant (PSI) with External Beam Radiation Therapy (EBRT) provides an attractive option for patients with intermediate or high-risk prostate disease as it delivers a higher radiobiological dose to the prostate compared to EBRT alone with similar level of anticipated toxicity [1] [2] . The standard approach of combined therapy includes the deliveries of seed implant and external-beam radiotherapy. The order of these two components to be executed could play an important role in term of treatment outcome. If EBRT is followed by PSI, it would be extremely difficult to compensate sub-optimal dose distribution delivered by PSI, if not impossible. On the other hand, PSI followed by EBRT is an attractive combination as the dose irregularities from PSI may be rectified by the subsequent EBRT [3] .
Traditional 3-dimensional conformal planning has limited capability for dose correction or dose reduction to surrounding normal organs. For instance, it is very difficult to boost the coverage to a low dose region inside the PTV with 3D conformal radiation therapy. However, the wide implementation of Intensity Modulated Radiation Therapy (IMRT) has been proved to overcome the high dose irregularities produced by prostate seed implant (PSI) [4] [5] . The drawback of IMRT is the increased dose inhomogeneity within the planning target volume, which could potentially increase the incidence of long-term side effects when coupled to the dose inhomogeneity inherited from PSI [6] - [8] .
In recent years, Proton Therapy (PT) has seen steadily growing popularity as an emerging external beam radiotherapy technique. Proton centers increased from 3 in 2001 to 13 in 2014 in the U.S. with a few more under construction or in the planning stage. PT is favored because of its sharp dose fall-off beyond the proton Bragg peak and its ability to precisely localize the radiation dosage. Similar to IMRT, intensity modulated proton therapy (IMPT) is a treatment delivery technique based on inverse planning optimization. The treatment planning system (TPS) evaluates the cost function of IMPT plan, adjusts the pencil beam weight, and achieves an optimized delivery plan. Among the IMPT techniques, spot-scanning-based IMPT is the most attractive modality, which modulates the intensity and energy of each small, mono-energetic proton beam spot individually under computer control. Different energy of the protons leads to different pencil beam range and energy spread of the protons determines the sharpness of Bragg Peak. One of the major advantages of spot scanning technique is that beam apertures and field-specific range compensators are not required, resulting in reduced neutron production and higher flexibility for adaptive planning [9] .
In this study, the feasibility of IMPT as a boost for I-125 permanent seed implant patients with localized prostate cancer was investigated. The results from IMPT were compared with those of IMRT plans. Several plan evaluation parameters such as the planning treatment volume (PTV) dose homogeneity, dose conformity and dose to surrounding normal tissues were evaluated to determine if IMPT can be a valuable option as a boost for PSI.
2. Method and Materials
2.1. Prostate Seed Implant
Ten PSI patients were randomly selected and evaluated retrospectively in this study. I-125 prostate seeds were 4.5 mm long and 0.8 mm in diameter, manufactured by Bard Medical, Inc (Covington, Georgia), with activities ranging from 0.3 mCi to 0.4 mCi. I-125 seed contains a titanium tube (0.08 mm thickness), a copper coating, an aluminum layer and a gold core marker. The seed implantation procedure which was previously described [10] used a real-time intra-operative planning approach with a Varian VariSeed treatment planning system and a prescription dose of 110 Gy. Patients were under general anesthesia during the implantation and ultrasound was used to guide needle placements. The dose criteria used in treatment planning included that prostate V100 was larger than 99%, V150 was between 50% and 65%, and D90 was below 125% of the prescribed dose.
2.2. Post-Implant CT Simulation
The post-implant CT scans were obtained with the Siemens SOMATOM CT scanner (Siemens, Erlangen, Germany) approximately one month after the seed implant to allow the decrease of prostate swelling. The slice thickness of the CT scan was 1 mm for the purpose of clear delineation of each radioactive seed. The quality of the PSI was evaluated based on these CT scans.
2.3. Treatment Planning
Both IMPT and IMRT plans were created based on the post-implant CT scans. The prostate and seminal vesicles were contoured as the clinical target volume (CTV). The planning target volume (PTV) was created with a margin of 5 mm around the CTV except that 3-mm margin was used at the prostate rectal interface to reduce the risk of rectal morbidity. Other critical organs such as bladder, rectum, penile bulb, femoral heads and urethra were also contoured. The urethra was defined as a 5-mm diameter volume using a surrogate technique as previously described [11] . The PTV coverage in plan optimization was such that 100% of the prescribed dose was required to cover at least 98% of the PTV. Rectum and bladder were spared as much as possible without compromising PTV coverage. The calculation grids were set at 1 mm and planning priority was given to PTV coverage. To accommodate the high density of I-125 radioactive seeds, the CT number to electron density table and CT number to relative stopping power table were both extended to CT number up to 6000 with linear extrapolation.
Proton treatment plans were created using intensity modulated proton beams from a Varian proton cyclotron accelerator with maximum energy of 250 MeV. Two parallel opposed lateral beams were generated in the Eclipse treatment planning system (RTP) (Varian Medical System, Palo Alto, CA, Varian version 10.0). The proton dose calculation was based on the convolution superposition algorithm with dose Kernels derived from Monte Carlo Calculation. The prescription dose was 45 Gray-equivalents (GyE) in 1.8-GyE fractions with an RBE of 1.1. IMRT boost plans with the same prescription dose and similar dose constraints were generated with seven co-planner 6 MV photon beams of a Trilogy Linac in the Eclipse RTP software Version 10.0 (Varian Medical System, Palo Alto, CA). Several planning-purpose structures were generated to assist the plan optimization to avoid hot spots and improve OARs sparing, including 1) Body minus PTV 0.5 cm expansion; 2) Rectum minus PTV; 3) Bladder minus PTV. Typical dose constraints are listed in Table 1. For comparison purpose, IMRT plans were also normalized to have 98% of the PTV volume covered by 100% of the prescribed dose.
2.4. Plan Evaluation
To compare the dose coverage between IMPT and IMRT plans, dose-volume histograms (DVH) were generated for quantitative evaluations. Specifically, the following parameters were compared between these two sets of plans:
・ PTV coverage: D1% (minimum dose delivered to 1% of PTV) and D99% (minimum dose delivered to 99% of PTV) were used as surrogate markers for maximum and minimum doses within PTV.
・ Conformity Index (CI) was used to evaluate the dose conformity of the two sets of plans compared to prescribed dose. The CI is defined as
(1)
where VPTV is the PTV volume (cc), VTV is the body volume that is covered by the prescribed isodose line, and TVPV is the PTV volume covered by the prescribed dose. In the ideal case, VPTV = VTV = TVPV and CI = 1 [12] [13] .
・ Homogeneity Index (HI) was used to assess the dose uniformity and is defined as the ratio of D5% to D95% in the PTV, where D5% is the maximum dose delivered to 5% of the PTV and D95% is the minimum dose delivered to 95% of the PTV. HI is always equal or greater than 1. Ideally when PTV is uniformly covered with a certain isodose line, HI equals to 1.0.
・ Organs at Risk (OARs), which included rectum, bladder, urethra and femoral heads. Dose distributions to these surrounding normal tissues were compared for both IMPT and IMRT plans.
・ Integral Dose, which is defined as the total dose received by the body volume in the CT scan minus the PTV volume.
Paired T test was performed to analyze the difference between IMPT and IMRT boost plans with p value less than 0.05 indicating statistical significance. SAS (version 9.2, SAS Institute, Cary, NC) statistical software was utilized for all the statistical analysis.
3. Results
Comparison of Dose Distributions between IMPT and IMRT Plans
The typical difference of dose distribution between IMPT and IMRT plans were illustrated in Figure 1 and Figure 2 for one of the patients in the study. The most striking difference was the much larger irradiated volume in the IMRT plan compared to IMPT. It was also obvious from Figure 1 that the low dose distribution to the rectum was significantly smaller in IMPT plans.
Table 1. Typical IMPT optimization criteria.
Figure 3 compared the DVHs of PTV and some Organ at Risk (OARs) between the IMRT and IMPT plan for the same patient.
PTV and Prostate Coverage
For the ten patients selected in this study, both IMRT and IMPT plans demonstrated excellent dosimetric coverage for PTV volumes. On average, 99.91% and 99.84% of the PTV volumes were covered by 95% of the prescribed dose for IMRT and IMPT plans, respectively. The dosimetric comparison of the PTV coverage between IMPT and IMRT plans for all patients were summarized in Table 1. The average value of CI for the ten IMPT plans was estimated to be 1.14, while the corresponding average CI for the IMRT plans was 1.10 (p value = 0.002). On the other hand, IMPT plans demonstrated slightly better homogeneity with an average HI of 1.04 compared to 1.05 for IMRT plans (p value = 0.04). No significant difference was observed in terms of the D99% and D1% values between IMPT and IMRT plans since p values were greater than 0.05.
Rectum
Compared to conventional IMRT boost plans, IMPT boost plans demonstrated improved sparing of rectum, while providing similar PTV coverage. The mean rectal dose in IMRT plans (21.0%) was significantly greater than that in IMPT plans (12.8%) with p < 0.05. The rectum D15%, D25%, D35% and D50% with IMPT plans were 31.2%, 13.7%, 5.9% and 1.6% respectively, compared to 44.1%, 28.2%, 18.1% and 11.4% respectively for the IMRT plans (all p values < 0.05). The D2cc values for the IMPT plans were not significantly different from those of the IMRT plans (p value = 0.643). A detailed comparison for rectum between IMPT and IMRT was listed in Table 2.
Bladder
IMPT also showed better bladder sparing. The mean dose of bladder decreased from 22.7% for IMRT plan to 18.2% (IMPT) (p < 0.05) of the prescribed dose. The bladder D35% and D50% with IMPT plans were 12.8% and 4.3% respectively, which were lower than 21.3% and 12.2% with IMRT plans (both p values < 0.05). D15%, D25% and D2cc were not significantly different between IMPT and IMRT plans (All p value > 0.05). The detailed comparison for bladder between IMPT and IMRT was listed in Table 3.
Urethra
No significant difference was observed for urethra dose between IMPT and IMRT plans. The mean dose of urethra in IMPT plans and IMRT plans were 104% to 103% of the prescribed dose respectively (p value = 0.06). Table 4 illustrated the detailed comparison of urethra and other normal tissues between IMPT and IMRT plans.
Pelvis (Integral Dose)
One of the drawbacks of IMRT is that larger volume of patient body may receive some radiation dose, which may potentially cause secondary cancer [14] . It was noted that the volume of pelvis receiving radiation doses was significantly reduced in IMPT plans (Table 5) compared to IMRT plans. The integral dose comparison was evaluated by the mean dose of body volume outside the PTV, which was also named “Body-PTV”. The mean doses to the volume of “Body-PTV” were estimated as 5.8% and 17.4% of the prescribed dose for IMPT plans and IMRT plans respectively (p value < 0.0001), proving much higher dose drop off of IMPT plans compared to IMRT plans.
Femoral Heads
It was found that the D5% of femoral heads comparison favored IMPT plans since the average D5% dose of IMPT plans was 40.2% of prescribed dose compared to 55.2% for IMRT plans.
Figure 3. Dose volume histogram of plan comparison between IMPT plan and IMRT plan.
Table 2. PTV coverage comparison between IMPT and IMRT.
Note: CI = Conformity Index; HI = Homogeneity Index.
Table 3. Rectum dose comparison between IMPT and IMRT.
Table 4. Bladder dose comparison between IMPT and IMRT.
Table 5. Normal tissue dose comparison between IMPT and IMRT.
4. Discussion
In this study, PSI with an IMPT boost proved to be a valuable option for prostate cancer patients. In combination with PSI, IMPT can deliver a significantly higher dose to gross disease to increase the likelihood of tumor control. Compared to IMRT boost plans, IMPT was able to deliver a comparable or better treatment in terms of normal tissue sparing.
The high density of implanted seeds has been a concern for IMPT as boost treatment. Even though it has been shown that high density markers could cause dose perturbation [15] - [17] , the effect is only limited to the areas around the markers. In the case of prostate seed implants, the under-dose area around seeds in a proton beam is not expected to be dosimetrically significant in overall dose delivery because the extremely high dose in the vicinity of radioactive seeds is more than enough to compensate the under dose in the shadows of the seeds in a proton beam. Analysis of the 10 seed implant cases showed that the mean dose is 199% of the prescribed dose at 2 mm away from each individual seed. The mean dose rate at 5 mm away from the seeds is 33% higher than the prescription dose. The biological equivalent dose (BED) is even higher in the immediate vicinity of the seeds. Therefore, affected dose coverage area was minimum and negligible considering the higher dose from seed implant.
It is also worthy to mention that a linear extrapolation of CT-relative stopping power curve (CT-RSP) was applied in this study. There have been plenty of literatures [18] - [20] exploring the importance of CT-RSP in the accuracy of proton dose calculation. However, for the same reason mentioned previously, the dose contribution to the region around the seed is extremely high from the implanted radioactive seeds. The high Z material of the seeds has very little effect in the combined overall dose distribution even though there may be under-dose region or certain inaccuracy of dose calculation around each seed in a proton beam.
Radiobiologically, a single rounded-off value of RBE of 1.1 for proton relative to megavoltage photon beams has been used in most proton treatment facilities. In this study, the comparison was mainly concerned with the relative dose distributions between IMPT and IMRT plans, which eliminated the issue of radiobiological dose summation and radiobiological difference. The analysis of BED difference between the LDR brachytherapy treatment and proton beam was beyond the scope of this study but could be solved by radiobiological planning in the future. Some commercial software, such as MIM Vista (Cleveland, Ohio), is capable to perform deformable dose evaluation and summation of Brachytherapy plans and EBRT plans using the built-in BED models.
In addition, only two laterally opposed proton beams for IMPT planning was investigated in this study. Further studies could be performed with different gantry angles and multiple beams for better optimization. It could improve the CI, HI and normal tissues sparing at the expense of the setup and delivery time. For patients with hip replacement, IMPT treatment with other gantry angle combination is also necessary due to hip prostheses effect.
5. Conclusion
Prostate Seed Implant (PSI) with an Intensity Modulation Proton Therapy (IMPT) boost can be a valuable option for prostate cancer patients. The IMPT boost delivers comparable or better treatment in terms of normal tissue sparing compared to IMRT plans. Future research will be focused on the implementation of IMPT on other types of cancer diseases.
Acknowledgements
The authors would like to acknowledge Xiaohe Yu for his assistance in this work.
NOTES
*Corresponding author.