Abstract
Background/Aim: Stereotactic body radiotherapy for prostate cancer using CyberKnife with circular cone requires a long treatment time. Raster scanning intensity modulated radiotherapy (RS-IMRT) has a potential of improving treatment efficacy, introducing shorter treatment time, better target dose uniformity, and lower organ at risk (OAR) dose. The purpose of the study was to develop a fluence optimization system for RS-IMRT. Patients and Methods: RS-IMRT plans were created for five prostate cancer patients treated with the Novalis system and parameters were compared to the Novalis treatment plans. From 80 nodes available for the CyberKnife, twelve nodes were arbitrarily selected. On the beam’s eye view of each beam, a 100×100 matrix of optimization points was created at the target center plane. The beam fluence map was optimized using the attraction-repulsion model (ARM). The beam fluence maps were converted to the scanning sequence using the ARM and a final dose calculation was performed. Results: For planning target volume (PTV), RS-IMRT plans showed higher dose covering 2% of the volume (D2%) and lower D98% compared to the Novalis plans. However, the homogeneity was within our Institutional clinical protocol. The RS-IMRT plans showed significantly lower OAR dose parameters including bladder volume receiving 100% of prescribed dose (V100%) and dose delivered to 5 cm3 of rectum (D5 cc). Conclusion: We developed a fluence optimization system for RS-IMRT that performs the entire RS-IMRT treatment planning process, including scanning sequence optimization and final dose calculation. The RS-IMRT was capable of generating clinically acceptable plans.
During the last decade, the use of stereotactic body radiotherapy (SBRT) with extreme hypofractionation for prostate cancer has been rapidly increasing owing to excellent clinical outcomes and cost-effectiveness (1-4). Several clinical trials have reported prostate-specific antigen (PSA) control outcomes of approximately 95% and low incidence of grade 3 and higher toxicities for low-risk prostate cancer (5-8). To reduce toxicity of organs at risk (OARs), high conformal dose distribution to the target and accurate beam delivery are essential. CyberKnife robotic radiosurgery system (Accuray Inc., Sunnyvale, CA, USA) has a linear accelerator mounted on a robotic arm to achieve conformal irradiation from multiple noncoplanar beam angles (9). Additionally, the orthogonal kilovoltage imaging system detects fiducial markers implanted in the prostate, and the 6-degree robotic couch allows for precise beam delivery.
The latest version of CyberKnife is available with InCise, a multi-leaf collimator (MLC) device (10, 11). Many studies have reported that InCise has capabilities of treating patients in a shorter treatment time compared to circular cones (12, 13). However, many Institutions still treat patients with the circular cones. Previously, we reported a concept of raster scanning intensity-modulated radiotherapy (RS-IMRT), a new scanning irradiation method for the CyberKnife (14). Although this was a simulation study, the scanning irradiation will be possible because the CyberKnife has great potential of dynamic robotic arm movement and beam delivery to enable motion tracking treatment for lung and liver tumors (15, 16). The scanning sequences were generated by mimicking the beam fluence maps of intensity modulated radiotherapy (IMRT) plans of a Novalis radiosurgery system (BrainLAB, Munich, Germany). The RS-IMRT showed potential in generating dose distribution comparable to the Novalis plans in a shorter treatment time compared to the sequential treatment with the CyberKnife. With small cone size, RS-IMRT showed more homogeneous target dose coverage and lower OAR dose compared to CyberKnife sequential plans (14). This technology has the potential of improving treatment efficacy of the CyberKnife in terms of short treatment time and better dose distribution without additional MLC device.
The purpose of this study was to develop an optimization system of the fluence maps for RS-IMRT. This work aims to deliver a planning system capable of performing entire planning processes of RS-IMRT without interference of commercial treatment planning systems (TPS).
Patients and Methods
Patients and treatment plans. Following Institutional Review Board approval, we analyzed data from five prostate cancer patients treated with the Novalis system, which equips a M3 microMLC with 3 mm leaf width. The patients were treated with a 6 MV photon beams from eight coplanar gantry angles. The treatment plans were developed according to our Institution’s protocol: the prescribed dose was 36.25 Gy/5 fr, normalized by the dose covering 95% of planning target volume (PTV D95%)=prescribed dose. Prescribed dose was targeted at 80% of PTV maximum dose (Dmax), with 70-90% being acceptable. Rectum was limited to dose delivered to 2 cm3 (D2 cc) <35 Gy, D5 cc <30 Gy, and volume receiving 50% of the organ (V50%) <40%. The bladder was limited to D10 cc <35 Gy, V50% <35 cm3, and V100% <5 cm3. Urethral dose was not evaluated in this study because it was not evaluated in the Novalis treatment plans. The prescribed dose of the Novalis plans were 72 Gy/30 fr. To evaluate the ability of the scanning irradiation developed in this study, the clinically used Novalis dose distribution was analyzed by rescaling to PTV D95%=36.25 Gy.
Raster scanning. Figure 1 shows a flow chart for RS-IMRT treatment planning. Targets, OARs, and shells for optimization were generated using iPLAN RT Image (ver. 4.1.2) TPS. Clinical target volume (CTV) was created with a 3 mm margin on the prostate, whereas 1 mm in posterior. PTV was created with a 2 mm margin on CTV. A 5 mm thick shell was created for optimization at 5 mm from the PTV.
Flowchart of treatment planning for the RS-IMRT. RS-IMRT: Raster-scanning intensity modulated radiotherapy; OAR: organ at risk.
The CyberKnife G4 system with a variable aperture Iris collimator (Accuray Inc.) was simulated in this study. For prostate cancer treatment, CyberKnife irradiation is possible from approximately 80 nodes. From these, twelve nodes were arbitrarily selected for treatment planning (Figure 2A). The nodes were selected so that a sufficient angle was opened between each beam. A 100×100-pixel beamlet was created on the beams-eye-view (BEV) of each beam. The resolution was set to 1 mm at the target center plane. An attraction-repulsion model (ARM) (17) was used to optimize the beamlets (Figure 2B). Each of the 100×100 pixels created for the 12 beams works as an optimization point, generating either attraction or repulsion. If the PTV dose is insufficient, the optimized point in the PTV generate an attraction force to increase the beam fluence. If the PTV dose is excessive, the optimization points generate repulsion force to decrease the beam fluence. If the OAR dose is excessive, the optimized points in the OAR generate a repulsive force to decrease the beam fluence.
An example case of treatment planning of the raster scanning IMRT: (A) Selection of nodes used for the treatment. (B) Optimization of fluence map. Red, green, and blue lines represent the differential DVH curves of PTV, rectum, and bladder, respectively. (C) Conversion of the fluence map to the scanning sequence. Spacing of allows represent the scanning speed. IMRT: Intensity-modulated radiotherapy; DVH: dose-volume histogram; PTV: planning target volume.
The obtained fluence map was converted to a scanning beam sequence for circular cone (Figure 2C). In this study, cones with a diameter of 10 mm were used. Control points were equally spaced at 5 mm intervals, half the cone diameter, and the dose rate at each control point was optimized using the ARM. Each optimization point generates an attraction and a repulsion force corresponding to the fluence map, which increases or decreases the dose rate of the scanning irradiation. The detail of the optimization of the scanning sequence is described elsewhere (14).
Three-dimensional dose distributions were calculated with the pencil beam algorithm using a ShioRIS 2.0 software. The calculation method and accuracy are described elsewhere (18). The dose distribution was normalized so that PTV D95%=prescription dose.
Evaluation. Dose-volume histogram (DVH) was calculated for PTV, bladder, and rectum; for PTV, D2%, D98%, and homogeneity index (HI) were evaluated (19). For bladder and rectum, DVH parameters were evaluated as clinical constraints and generalized equivalent uniform dose (gEUD) were calculated. HI and gEUD were calculated by the following equations:
(1)
(2)
where n is a parameter representing the volume dependence of the dose-response relationship for each organ. n values were 0.12 and 0.5 for the rectum and bladder, respectively, as previously reported by Burman et al. (20).
Statistical analysis. Differences between Novalis and RS-IMRT plans were compared using the paired t-test. Statistical significance was defined as a p-Value <0.05.
Results
The RS-IMRT plans were generated for five patients. The PTV were 64.1, 122.7, 74.2, 111.3, and 62.6 cm3 for patient #1, 2, 3, 4, and 5, respectively. Mean±standard deviation of the number of control points were 133,350±34,933. For all cases, the optimization of the fluence map completed within 10 s.
Figure 3A shows an example of dose distribution calculated with RS-IMRT, and Figure 3B shows an DVH of an example case. The solid and dash lines represent the RS-IMRT and Novalis plans, respectively. The Novalis plan showed steeper DVH curve of PTV than that of RS-IMRT, indicating better dose uniformity. On the other hand, in this case, the scanning irradiation plan achieved lower doses in the high-dose areas of the bladder and rectum.
(A) An example dose distribution of RS-IMRT. (B) DVH of an example case (#5). Solid and dash lines represent RS-IMRT and Novalis plans, respectively. RS-IMRT: Raster-scanning intensity-modulated radiotherapy; DVH, dose-volume histogram; PTV: planning target volume.
Figure 4 shows the results for the DVH parameters. The mean (min-max) values for each data are listed in Table I. The RS-IMRT showed significantly lower target homogeneity compared to Novalis. However, the values were within the tolerance of our clinical protocol. RS-IMRT showed significantly lower bladder V100%. and rectum D5 cc. For other parameters, significantly differences were not observed.
DVH parameters of Novalis and RS-IMRT plans calculated for each patient. DVH: Dose-volume histogram; RS-IMRT: raster-scanning intensity-modulated radiotherapy; PTV: planning target volume; Dx% (or x cc): dose delivered to the volume of x% (or x cm3); V100%: volume which receives 100% prescribed dose.
DVH parameters of Novalis and RS-IMRT plans.
Discussion
In this study, we developed an optimization system for RS-IMRT for prostate SBRT. In our previous study, we demonstrated that the RS-IMRT achieved good dose distribution with short treatment time (14). However, in the previous study, the scan sequence was created by mimicking the IMRT fluence map created by iPLAN TPS. The treatment planning could not be performed only with the system alone. In the present study, we simulated beams from nodes that can actually be irradiated with the CyberKnife and created plans with head motion velocities in the range feasible with the CyberKnife. This study enabled a more realistic simulation, and the system developed in this study enabled to optimize the fluence, convert it to a scanning sequence, and perform dose calculation with this system. If this technology is incorporated into the CyberKnife, Institutions without InCise MLC would be able to provide better treatment to patients.
The RS-IMRT plans showed superior OAR dose reduction compared to the Novalis plans. This would be due to the improved dose conformity achieved by using noncoplanar beams. In the previous study, Novalis plans showed lower rectal and bladder doses compared to the RS-IMRT plans, probably because the scanning beam was designed to reproduce the fluence map of Novalis plan and thus could not produce a better plan than the Novalis plan.
On the other hand, the uniformity of PTV is significantly worse than that of Novalis plans. This would be because the beamlets produced by Novalis are made of MLCs and have no joints, whereas the beamlets produced by raster scanning are made of cones overlapping by half their diameter, so the boundaries between beams are not uniform, even though they overlap. Furthermore, because the control point is set at the center plane of the target, the size of the cone is smaller upstream of the beam, which is expected to increase inhomogeneity. However, the dose uniformity of the PTV was within our criteria, and the dose distribution seemed clinically acceptable.
Tomida et al. and Masi et al. previously reported that InCize MLC shortens treatment time by reducing monitor units (MU) and nodes compared to the variable aperture Iris collimator (12, 13). Furthermore, a new optimization system, VOLO, has been reported to reduce the number of nodes and MUs compared to the conventional sequential optimization, enabling high quality treatment in a short time (21-24).
RS-IMRT also achieves good dose distribution with a much smaller number of nodes. In addition, as previously reported, combining a static beam with a large cone diameter significantly shorten the treatment time (14). This technique would be useful for improving treatment efficacy for facilities which do not obtain MLC-equipped machine.
The limitation of this study is that it is a virtual simulation and the current version of the CyberKnife cannot perform such a movement. However, CyberKnife itself has sufficient potential of performing scanning irradiation. The RS-IMRT would become possible if the software is updated to support this.
Conclusion
In this study, we developed an optimization system for RS-IMRT in SBRT for prostate cancer. Our approach allowed the entire planning process of fluence optimization, conversion to scanning sequences, and final dose calculation to be completed within the system. RS-IMRT produced a clinically acceptable treatment plan with fewer nodes than sequential treatment plans. Compared to sequential treatment, the CyberKnife IMRT planning system proposed herein has therapeutic advantages of: (i) shorter treatment times due to fewer nodes, (ii) better uniformity of target dose, and (iii) lower OAR dose. The proposed RS-IMRT technique has the potential to provide a treatment showcasing these benefits without an additional MLC device.
Acknowledgements
This work was supported by JSPS KAKENHI Grant Number JP20K08022.
Footnotes
Authors’ Contributions
Data analysis and investigation: Yuichi Akino and Hiroya Shiomi; Funding: Yuichi Akino; Supervision: Nobuhisa Mabuchi, Norihisa Masai, Ryoong-Jin, and Kazuhiko Ogawa.
Conflicts of Interest
The Author HS is a developer of the ShioRIS 2.0 software.
- Received December 20, 2022.
- Revision received January 6, 2023.
- Accepted January 23, 2023.
- Copyright © 2023 International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved.










