

ORIGINAL ARTICLE 



Year : 2016  Volume
: 41
 Issue : 2  Page : 123128 

Estimation of the effects of normal tissue sparing using equivalent uniform dosebased optimization
K Senthilkumar^{1}, KJ Maria Das^{2}, K Balasubramanian^{3}, AC Deka^{3}, BR Patil^{3}
^{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  30Jan2016 
Date of Decision  04Mar2016 
Date of Acceptance  28Mar2016 
Date of Web Publication  3May2016 
Correspondence Address: K Senthilkumar Department of Medical Physics, Karnataka Cancer Therapy and Research Institute, Navanagar, Hubli  580 025, Karnataka India
Source of Support: None, Conflict of Interest: None  Check 
DOI: 10.4103/09716203.181631
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) EUDbased optimization with a DVbased physical cost function to control inhomogeneity (EUD_{With DV}) and (ii) EUDbased optimization without a DVbased physical cost function to allow inhomogeneity (EUD_{Without DV}). The generated plans were prescribed a dose of 72 Gy in 36 fractions to planning target volume (PTV). Mean dose, D_{30%}, and D_{5%}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 cm^{3}. The PTV mean dose for EUD_{With DV}plans was 73.67 ± 1.7 Gy, whereas for EUD_{Without DV}plans was 80.42 ± 2.7 Gy. It was found that PTV volume receiving dose more than 115% of prescription dose was negligible in EUD_{With DV} plans, whereas it was 28% in EUD_{Without DV} plans. In almost all dosimetric parameters evaluated, dose to OARs in EUD_{With DV}plans was higher than in EUD_{Without DV}plans. Allowing inhomogeneous dose (EUD_{Without DV}) inside the target would achieve better normal tissue sparing compared to homogenous dose distribution (EUD_{With DV}). Hence, this inhomogeneous dose could be intentionally dumped on the highrisk 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; inhomogeneity; 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 dosebased optimization. J Med Phys 2016;41:1238 
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 dosebased optimization. J Med Phys [serial online] 2016 [cited 2021 Mar 5];41:1238. Available from: https://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 DVbased 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 DVbased 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.^{[1]} Furthermore, a single DVbased physical cost function for tumor does not represent the real nature of dose response of tumor. However, it can be argued that adding multiple DVbased physical cost functions for tumor would represent the dose response nature of the tumor but to a lesser extent. Although EUDbased cost functions are highly degenerative, it might also pose some clinically unaccepted problems.^{[2]} For instance, it may lead to very inhomogeneous target dose distribution. This target dose inhomogeneity was controlled by adding DVbased 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 EUDbased cost functions.
Materials and Methods   
Planning techniques
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) EUDbased optimization with a DVbased physical cost function to control inhomogeneity (EUD_{With DV}), a combination of cost functions for target which is widely used in day to day clinical practice and (ii) EUDbased optimization without a DVbased physical cost function to allow inhomogeneity (EUD_{Without 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.^{[3]} For all Monaco plans, Xray Volume Monte Carlo (XVMC) algorithm with 3% variance was used in the segment shape optimization phase.^{[4]} 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.^{[1]} 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, D_{i} is the dose in the i^{th} 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.^{[5]} 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.,^{[2]} EUDbased cost function produces better normal structure sparing over DVbased cost function for the same minimal target dose in IMRT plans. The same has been demonstrated by several authors.^{[6],[7],[8]} Further, EUDbased cost functions are insensitive to hot spots inside the tumor which leads to highly inhomogeneous target dose, if used alone.^{[2]} In clinical practice, a physical cost function (DVbased) 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 biologicalbased optimization features. It offers three biologicalbased cost functions namely; Poisson statistics cell kill model for target; serial complication model and parallel complication model for OAR. Monaco also offers several physical DVbased cost functions. The biological cost functions incorporated into Monaco TPS were developed by Alber and Reemtsen.^{[9]} Detailing the full mathematics of their work is beyond the scope of this paper. For each cost function, a threedimensional dose distribution is reduced to a single index called isoeffect.^{[10]} 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, D_{30%}, and D_{5%}
 Bladder  Mean dose and D_{5%}
 Left femur head  D_{5%}
 Right femur head  D_{5%}.
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.^{[11]} LymanKutcherBurman NTCP model was used for calculation of NTCP.^{[12],[13]} 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.
Statistical analysis
To determine the statistical significance, twotailed paired ttests 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 EUD_{With DV} plans was 73.67 ± 1.7 Gy, whereas for EUD_{Without DV} plans was 80.42 ± 2.7 Gy. The PTV volume receiving dose more than 107% of the prescription dose in EUD_{With DV} plans was 2.1 ± 0.6 cm ^{3} while for EUD_{Without DV} plans was 40.6 ± 3.8 cm ^{3} (40.8%). Similarly, PTV volume receiving dose more that 115% of the prescription dose in EUD_{With DV} plans was negligible, whereas for EUD_{Without 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 EUD_{With DV}plan and EUD_{Without DV}plan 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 EUD_{With DV} plans was 5.6 Gy higher than EUD_{Without DV} plans. The same correlation continues in D_{30%} also. For bladder, the mean dose in EUD_{With DV} plan was higher than EUD_{Without DV} plans. For rectum and bladder, D_{5%} dose difference was not statistically significant between both the plans. The D_{5%} dose for right femur in EUD_{With DV} was 5.01 Gy higher than EUD_{Without DV} plans (p<0.002). Similarly, the D_{5%} dose for left femur in EUD_{With DV} was 4.8 Gy higher than EUD_{Without 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 EUD_{With DV,} whereas for EUD_{Without DV}, it was 3.3 ± 0.24% with P = 0.0028. For bladder, the calculated NTCP was 1.25 ± 0.19% for EUD_{With DV} plans, whereas for EUD_{Without DV} was 0.84 ± 0.16% with P = 0.0017. The calculated NTCP for left femur was 0.28 ± 0.11% for EUD_{Without DV} plans, whereas for EUD_{With DV} was0.70 ± 0.15% with P = 0.001. For right femur, the calculated NTCP was 0.27 ± 0.10% for EUD_{Without DV} plans, whereas for EUD_{With 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 (EUD_{With DV}), whereas for inhomogeneous (EUD_{Without DV}) plans, it was 9.12 ± 1.9 min. The calculated mean MU was 456 ± 39 for homogenous plans (EUD_{With DV}) and 610 ± 42 for inhomogeneous plans (EUD_{Without DV}).
Discussion   
Plans generated without DV cost function in the EUDbased 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 EUD_{Without DV} plans. This inhomogeneous dose could be deliberately dumped to highrisk volume to achieve high local control. This opens up the possibility of selectively boosting the substantial volume of the tumor using EUDbased cost function for IMRT optimization.
Mean dose for both rectum and bladder was high in EUD_{With DV} plans compared to EUD_{Without DV} plans. For both femoral heads, D_{5%} was high in EUD_{With DV} plans compared to EUD_{Without DV} plans. It was found that by allowing inhomogeneous dose (EUD_{Without DV} plans) inside the target, it was possible to achieve reduced OAR dose compared to homogenous dose distribution (EUD_{With 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 EUDbased 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.^{[2]} that introducing a DVbased physical cost function with EUDbased 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 EUDbased IMRT optimization.
At the same time, for both rectum and bladder, there was no statistically significant difference between EUD_{With DV} and EUD_{Without DV} plans in D_{5} doses. It is obvious that D_{5%} 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 EUD_{Without DV} plans, D_{5%} dose was slightly high compared to EUD_{With DV} because of increased dose inhomogeneity inside the target.
Total treatment time was significantly high in inhomogeneous dose (EUD_{Without 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 ^{[14]} 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 EUDbased 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 (EUD_{Without DV}) inside the target one can achieve better normal tissue sparing as compared to homogenous dose distribution (EUD_{With DV}). Hence, this inhomogeneous dose could be intentionally dumped on the highrisk 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
Nil.
Conflicts of interest
There are no conflicts of interest.
References   
1.  Niemierko A. Reporting and analyzing dose distributions: A concept of equivalent uniform dose. Med Phys 1997;24:10310. 
2.  Wu Q, Mohan R, Niemierko A, SchmidtUllrich R. Optimization of intensitymodulated radiotherapy plans based on the equivalent uniform dose. Int J Radiat Oncol Biol Phys 2002;52:22435. 
3.  Bertsekas DP. Contrained Optimization and Lagrange Multiplier Methods. Belmont, Massachusetts: Athena Scientific; 1996. 
4.  Fippel M. Fast Monte Carlo dose calculation for photon beams based on the VMC electron algorithm. Med Phys 1999;26:146675. 
5.  Abramowitz M, Stegun IA. Handbook of Mathematical Functions. Washington, D.C: National Bureau of Standards; 1968. 
6.  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:585160. 
7.  Qi XS, Semenenko VA, Li XA. Improved critical structure sparing with biologically based IMRT optimization. Med Phys 2009;36:17909. 
8.  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 basedIMRT. Technol Cancer Res Treat 2011;10:57583. 
9.  Alber M, Reemtsen R. Intensity modulated radiation therapy planning by use of a barrierpenalty multiplier method. Optim Methods Softw 2007;22:391411. 
10.  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 TG166 of the Therapy Physics Committee of the AAPM. Med Phys 2012;39:1386409. 
11.  Nahum AE, Uzan J. (Radio) biological optimization of externalbeam radiotherapy. Comput Math Methods Med 2012;2012:329214. 
12.  Lyman JT, Wolbarst AB. Optimization of radiation therapy, III: A method of assessing complication probabilities from dosevolume histograms. Int J Radiat Oncol Biol Phys 1987;13:1039. [ PUBMED] 
13.  Lyman JT. Complication probability as assessed from dosevolume histograms. Radiat Res Suppl 1985;8:S139. [ PUBMED] 
14.  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:51922. [ PUBMED] 
[Figure 1]
[Table 1], [Table 2], [Table 3]
