|Year : 2016 | Volume
| Issue : 2 | Page : 123-128
Estimation of the effects of normal tissue sparing using equivalent uniform dose-based optimization
K Senthilkumar1, KJ Maria Das2, K Balasubramanian3, AC Deka3, BR Patil3
1 Department of Medical Physics, Karnataka Cancer Therapy and Research Institute, Hubli, Karnataka; Research and Development Centre, Bharathiar University, Coimbatore, Tamil Nadu, India
2 Department of Radiotherapy, Sanjay Gandhi Post Graduate Institute of Medical Sciences, Lucknow, Uttar Pradesh, India
3 Department of Medical Physics, Karnataka Cancer Therapy and Research Institute, Hubli, Karnataka, India
|Date of Submission||30-Jan-2016|
|Date of Decision||04-Mar-2016|
|Date of Acceptance||28-Mar-2016|
|Date of Web Publication||3-May-2016|
Department of Medical Physics, Karnataka Cancer Therapy and Research Institute, Navanagar, Hubli - 580 025, Karnataka
Source of Support: None, Conflict of Interest: None
| Abstract|| |
In this study, we intend to estimate the effects of normal tissue sparing between intensity modulated radiotherapy (IMRT) treatment plans generated with and without a dose volume (DV)-based physical cost function using equivalent uniform dose (EUD). Twenty prostate cancer patients were retrospectively selected for this study. For each patient, two IMRT plans were generated (i) EUD-based optimization with a DV-based physical cost function to control inhomogeneity (EUDWith DV) and (ii) EUD-based optimization without a DV-based physical cost function to allow inhomogeneity (EUDWithout DV). The generated plans were prescribed a dose of 72 Gy in 36 fractions to planning target volume (PTV). Mean dose, D30%, and D5%were evaluated for all organ at risk (OAR). Normal tissue complication probability was also calculated for all OARs using BioSuite software. The average volume of PTV for all patients was 103.02 ± 27 cm3. The PTV mean dose for EUDWith DVplans was 73.67 ± 1.7 Gy, whereas for EUDWithout DVplans was 80.42 ± 2.7 Gy. It was found that PTV volume receiving dose more than 115% of prescription dose was negligible in EUDWith DV plans, whereas it was 28% in EUDWithout DV plans. In almost all dosimetric parameters evaluated, dose to OARs in EUDWith DVplans was higher than in EUDWithout DVplans. Allowing inhomogeneous dose (EUDWithout DV) inside the target would achieve better normal tissue sparing compared to homogenous dose distribution (EUDWith DV). Hence, this inhomogeneous dose could be intentionally dumped on the high-risk volume to achieve high local control. Therefore, it was concluded that EUD optimized plans offer added advantage of less OAR dose as well as selectively boosting dose to gross tumor volume.
Keywords: Biological optimization; equivalent uniform dose; in-homogeneity; intensity modulated radiotherapy
|How to cite this article:|
Senthilkumar K, Maria Das K J, Balasubramanian K, Deka A C, Patil B R. Estimation of the effects of normal tissue sparing using equivalent uniform dose-based optimization. J Med Phys 2016;41:123-8
|How to cite this URL:|
Senthilkumar K, Maria Das K J, Balasubramanian K, Deka A C, Patil B R. Estimation of the effects of normal tissue sparing using equivalent uniform dose-based optimization. J Med Phys [serial online] 2016 [cited 2019 Jun 24];41:123-8. Available from: http://www.jmp.org.in/text.asp?2016/41/2/123/181631
| Introduction|| |
With the advancement of radiotherapy, treatment plan should be optimized to produce desired dose distribution inside the tumor with reduced normal tissue dose. This could be achieved using intensity modulated radiotherapy (IMRT). Nowadays, most of the IMRT treatment planning systems (TPS) incorporate the dose volume (DV)-based physical cost functions for IMRT optimization. The major drawback associated with the DV-based physical cost function used in IMRT optimization is that it does not represent the nonlinear response of tumor or normal tissues. Further IMRT plan score does not get affected by small cold spot inside the tumor when using DV-based physical cost functions in optimization. On the other hand, plan score for IMRT plan based on equivalent uniform dose (EUD) cost functions would be significantly diminished if there is cold spot inside the tumor. Furthermore, a single DV-based physical cost function for tumor does not represent the real nature of dose response of tumor. However, it can be argued that adding multiple DV-based physical cost functions for tumor would represent the dose response nature of the tumor but to a lesser extent. Although EUD-based cost functions are highly degenerative, it might also pose some clinically unaccepted problems. For instance, it may lead to very inhomogeneous target dose distribution. This target dose inhomogeneity was controlled by adding DV-based cost functions by limiting the real potential of EUD dose distributions in clinical practice. In this study, we intend to compare the effects of dose distributions in organ at risk (OAR) estimated from IMRT treatment plans generated with and without inhomogeneous dose distributions using EUD-based cost functions.
| Materials and Methods|| |
Twenty patients treated for prostate cancer were retrospectively selected for this study and all these patients underwent radiotherapy computed tomography (CT) scans of 3 mm slice thickness, extending from the second lumbar vertebrae to proximal third of femoral diaphysis. T2 weighted magnetic resonance imaging scans were fused to the CT images for delineation of gross target volume (GTV) which included visible prostate and clinical target volume (CTV). CTV includes prostate plus seminal vesicle. Planning target volume (PTV) was expanded nonuniformly from the CTV. The contours for OARs such as rectum, bladder, and femoral heads were delineated from the CT data.
For each patient, two IMRT plans were generated (i) EUD-based optimization with a DV-based physical cost function to control inhomogeneity (EUDWith DV), a combination of cost functions for target which is widely used in day to day clinical practice and (ii) EUD-based optimization without a DV-based physical cost function to allow inhomogeneity (EUDWithout DV), a combination of cost functions for target which we are proposing to estimate the effects of dose distributions in OAR. The generated plans were prescribed a dose of 72 Gy in 36 fractions to PTV. The treatment plan acceptance criterion was to deliver 95% of the prescribed dose to 95% volume of PTV. Treatment planning was generated using Monaco TPS version 5.0 (CMS Inc., St. Louis, MO, USA). Serial, parallel, and physical cost functions were used for OARs in all plans. To make both treatment plans identical, same kind of cost functions were used for both the plans for particular OAR. Both treatment plans used same beam energy, number of beams, beam angles, and isocenter. Monaco used 3 mm grid spacing for dose calculation. Monaco TPS is based on constrained optimization method. For all Monaco plans, X-ray Volume Monte Carlo (XVMC) algorithm with 3% variance was used in the segment shape optimization phase. All plans were corrected for tissue heterogeneities. The bladder and rectum dose reporting were done for entire organ. Iterative adjustment of isoconstraint values followed until the mandated dosimetric criteria was achieved.
Biological optimization in Monaco TPS
The EUD concept was introduced by Niemierko for tumor and normal tissues as the biological EUD, if given uniformly, would result in same biological effect as the actual nonuniform dose distribution. The phenomenological form of EUD,
The above equation applies to both tumors and normal tissues. In this equation, N is the number of voxels in the anatomical structure of interest, Di is the dose in the ith voxel, and a is the tumor or normal tissue specific parameter that describes the DV effect. This EUD formalism is based on the power law dependence which stimulus the response of complex biological system. The above expression is generalized mean of nonuniform dose distribution. For a = ∞, the EUD is equal to the maximum dose, and for a = −∞, the EUD is equal to minimum dose. For a = 1, the EUD is equal to arithmetic mean, and for a = 0, it is equal to geometric mean.
As stated by Wu et al., EUD-based cost function produces better normal structure sparing over DV-based cost function for the same minimal target dose in IMRT plans. The same has been demonstrated by several authors.,, Further, EUD-based cost functions are insensitive to hot spots inside the tumor which leads to highly inhomogeneous target dose, if used alone. In clinical practice, a physical cost function (DV-based) is added with EUD cost function to achieve homogenous dose distribution inside the tumor by accepting higher OAR dose.
Monaco TPS is the first commercial IMRT TPS that incorporated biological-based optimization features. It offers three biological-based cost functions namely; Poisson statistics cell kill model for target; serial complication model and parallel complication model for OAR. Monaco also offers several physical DV-based cost functions. The biological cost functions incorporated into Monaco TPS were developed by Alber and Reemtsen. Detailing the full mathematics of their work is beyond the scope of this paper. For each cost function, a three-dimensional dose distribution is reduced to a single index called isoeffect. On the other hand, clinical goals specified by the user are referred to as isoconstraint.
Isoeffect for target calculated using Poisson cell kill model is as follows:
Where α′ is the average cell sensitivity, ρ′ is the average clonogen density, V is the total volume of the organ, and is a biological response function given by,
Where is the local density of clonogenic tumor cells, is the cell sensitivity for particular voxel, and is the absorbed dose in the particular voxel. At present, user can only specify the parameter cell sensitivity ranges from 0.1 to 1.0 Gy −1. Equation (3) conceptually represents the EUD formalism which was discussed earlier.
Physical dose evaluation indices
The cumulative DV histograms (DVH) parameters were reported for the following:
- Rectum - Mean dose, D30%, and D5%
- Bladder - Mean dose and D5%
- Left femur head - D5%
- Right femur head - D5%.
The treatment time and total monitor units (MU) were also compared.
Biological dose evaluation indices
Normal tissue complication probability (NTCP) values were calculated using BioSuite software proposed by Nahum and Uzan. Lyman-Kutcher-Burman NTCP model was used for calculation of NTCP., Although BioSuite software offers a list of endpoints with default parameters extracted from literature, it is possible for users to use their own experimental data. We used the default BioSuite end points for all our cases to calculate NTCP. Differential DVH for OAR for both the plans was exported from Monaco TPS and was converted to a BioSuite compatible format. BioSuite plan was generated with same planning parameters used in Monaco TPS including total dose, fraction size, and number of fractions. Differential DVHs were imported into the BioSuite software and corresponding endpoints were associated with respective DVHs. NTCP values were calculated for OARs for all patients.
To determine the statistical significance, two-tailed paired t-tests were performed with P < 0.05 considered to be statistically significant. All calculations were performed using the online statistical packages software called VassarStats (Vassar College, Poughkeepsie, NY, USA).
| Results|| |
The average volume (mean ± standard deviation) of PTV for all patients was 103.02 ± 27.03 cm 3. The PTV mean dose for EUDWith DV plans was 73.67 ± 1.7 Gy, whereas for EUDWithout DV plans was 80.42 ± 2.7 Gy. The PTV volume receiving dose more than 107% of the prescription dose in EUDWith DV plans was 2.1 ± 0.6 cm 3 while for EUDWithout DV plans was 40.6 ± 3.8 cm 3 (40.8%). Similarly, PTV volume receiving dose more that 115% of the prescription dose in EUDWith DV plans was negligible, whereas for EUDWithout DV plans was 27.9 ± 1.7 cm 3 (27%). [Table 1] shows the target volume (X) receiving the dose (Y). [Figure 1] (a) to (e) shows the cumulative DVH curves of a typical prostate cancer patient.
|Figure 1: Comparison of the dose volume histograms between EUDWith DVplan and EUDWithout DVplan for (a) PTV_72Gy (b) Rectum (c) Bladder (d) Left femoral head and (e) Right femoral head of a typical prostate cancer patient|
Click here to view
The mean dose for rectum in EUDWith DV plans was 5.6 Gy higher than EUDWithout DV plans. The same correlation continues in D30% also. For bladder, the mean dose in EUDWith DV plan was higher than EUDWithout DV plans. For rectum and bladder, D5% dose difference was not statistically significant between both the plans. The D5% dose for right femur in EUDWith DV was 5.01 Gy higher than EUDWithout DV plans (p<0.002). Similarly, the D5% dose for left femur in EUDWith DV was 4.8 Gy higher than EUDWithout DV plans (p<0.002). [Table 2] summarizes the dosimetric values of all OARs evaluated.
|Table 2: Comparison of evaluated dosimetric values for all organs at risk|
Click here to view
The calculated NTCP for rectum was4.4 ± 0.25% for EUDWith DV, whereas for EUDWithout DV, it was 3.3 ± 0.24% with P = 0.0028. For bladder, the calculated NTCP was 1.25 ± 0.19% for EUDWith DV plans, whereas for EUDWithout DV was 0.84 ± 0.16% with P = 0.0017. The calculated NTCP for left femur was 0.28 ± 0.11% for EUDWithout DV plans, whereas for EUDWith DV was0.70 ± 0.15% with P = 0.001. For right femur, the calculated NTCP was 0.27 ± 0.10% for EUDWithout DV plans, whereas for EUDWith DV was 0.64 ± 0.15% with P = 0.001. The calculated NTCP values are shown in [Table 3].
|Table 3: Calculated normal tissue complication probability values for all organs at risk|
Click here to view
The mean treatment delivery time was 7.2 ± 2.3 min for homogenous plans (EUDWith DV), whereas for inhomogeneous (EUDWithout DV) plans, it was 9.12 ± 1.9 min. The calculated mean MU was 456 ± 39 for homogenous plans (EUDWith DV) and 610 ± 42 for inhomogeneous plans (EUDWithout DV).
| Discussion|| |
Plans generated without DV cost function in the EUD-based optimization resulted in inhomogeneous dose distribution inside the PTV. It was found that the PTV volume receiving dose more than 107% of prescription dose was around 40% and the PTV volume receiving dose more than 115% of the prescription dose was around 28% in EUDWithout DV plans. This inhomogeneous dose could be deliberately dumped to high-risk volume to achieve high local control. This opens up the possibility of selectively boosting the substantial volume of the tumor using EUD-based cost function for IMRT optimization.
Mean dose for both rectum and bladder was high in EUDWith DV plans compared to EUDWithout DV plans. For both femoral heads, D5% was high in EUDWith DV plans compared to EUDWithout DV plans. It was found that by allowing inhomogeneous dose (EUDWithout DV plans) inside the target, it was possible to achieve reduced OAR dose compared to homogenous dose distribution (EUDWith DV plans). Even calculated NTCP values for all OARs substantiate the reduced OAR dose with inhomogeneous dose inside the target. Therefore, it was noted that EUD-based IMRT optimized plans offer added advantage of less OAR dose as well as selectively boosting dose to gross tumor volume. It was already stated by Wu et al. that introducing a DV-based physical cost function with EUD-based optimization to get homogenous dose would degrade the dose distributions in OARs. In this study, we estimated the effects of dose distributions in OARs using EUD-based IMRT optimization.
At the same time, for both rectum and bladder, there was no statistically significant difference between EUDWith DV and EUDWithout DV plans in D5 doses. It is obvious that D5% dose was evaluated for full organ, i.e., some part of the rectum and bladder was inside the PTV. Although we tried to prevent the overlapping OAR volume from receiving more than the prescription dose, still for EUDWithout DV plans, D5% dose was slightly high compared to EUDWith DV because of increased dose inhomogeneity inside the target.
Total treatment time was significantly high in inhomogeneous dose (EUDWithout DV) plans because of more total energy required to deliver such high dose. Consequently, treatment time was also high to deliver inhomogeneous dose to the target. One should note that the absolute treatment time may vary depending on how efficient the sequencing of the segments and therefore be TPS dependent.
In clinical practice, there is a concern for accepting inhomogeneous dose inside the tumor. Goitein and Niemierko  stated that inhomogeneous dose can be accepted if it is not due to treatment delivery methods. However, further study is needed on how much tumor volume can be selectively boosted by EUD-based cost functions optimized plans. The hypothesis of this study should also be tested where we practice simultaneous integrated boosts such as in the treatment of head and neck cancer.
| Conclusion|| |
This study demonstrated that by allowing inhomogeneous dose (EUDWithout DV) inside the target one can achieve better normal tissue sparing as compared to homogenous dose distribution (EUDWith DV). Hence, this inhomogeneous dose could be intentionally dumped on the high-risk volume to achieve high local control. Therefore, it is concluded that EUD optimized plans offer added advantage of less OAR dose as well as selectively boosting dose to gross tumor volume.
Financial support and sponsorship
Conflicts of interest
There are no conflicts of interest.
| References|| |
Niemierko A. Reporting and analyzing dose distributions: A concept of equivalent uniform dose. Med Phys 1997;24:103-10.
Wu Q, Mohan R, Niemierko A, Schmidt-Ullrich R. Optimization of intensity-modulated radiotherapy plans based on the equivalent uniform dose. Int J Radiat Oncol Biol Phys 2002;52:224-35.
Bertsekas DP. Contrained Optimization and Lagrange Multiplier Methods. Belmont, Massachusetts: Athena Scientific; 1996.
Fippel M. Fast Monte Carlo dose calculation for photon beams based on the VMC electron algorithm. Med Phys 1999;26:1466-75.
Abramowitz M, Stegun IA. Handbook of Mathematical Functions. Washington, D.C: National Bureau of Standards; 1968.
Semenenko VA, Reitz B, Day E, Qi XS, Miften M, Li XA. Evaluation of a commercial biologically based IMRT treatment planning system. Med Phys 2008;35:5851-60.
Qi XS, Semenenko VA, Li XA. Improved critical structure sparing with biologically based IMRT optimization. Med Phys 2009;36:1790-9.
Anderson N, Lawford C, Khoo V, Rolfo M, Joon DL, Wada M. Improved normal tissue sparing in head and neck radiotherapy using biological cost function based-IMRT. Technol Cancer Res Treat 2011;10:575-83.
Alber M, Reemtsen R. Intensity modulated radiation therapy planning by use of a barrier-penalty multiplier method. Optim Methods Softw 2007;22:391-411.
Allen Li X, Alber M, Deasy JO, Jackson A, Ken Jee KW, Marks LB, et al.
The use and QA of biologically related models for treatment planning: Short report of the TG-166 of the Therapy Physics Committee of the AAPM. Med Phys 2012;39:1386-409.
Nahum AE, Uzan J. (Radio) biological optimization of external-beam radiotherapy. Comput Math Methods Med 2012;2012:329214.
Lyman JT, Wolbarst AB. Optimization of radiation therapy, III: A method of assessing complication probabilities from dose-volume histograms. Int J Radiat Oncol Biol Phys 1987;13:103-9.
Lyman JT. Complication probability as assessed from dose-volume histograms. Radiat Res Suppl 1985;8:S13-9.
Goitein M, Niemierko A. Intensity modulated therapy and inhomogeneous dose to the tumor: A note of caution. Int J Radiat Oncol Biol Phys 1996;36:519-22.
[Table 1], [Table 2], [Table 3]