Skip to main content

An integrated strategy of biological and physical constraints in biological optimization for cervical carcinoma



For cervical carcinoma cases, this study aimed to evaluate the quality of intensity-modulated radiation therapy (IMRT) plans optimized by biological constraints. Furthermore, a new integrated strategy in biological planning module was proposed and verified.


Twenty patients of advanced stage cervical carcinoma were enrolled in this study. For each patient, dose volume optimization (DVO), biological model optimization (BMO) and integrated strategy optimization (ISO) plans were created using same treatment parameters. Different biological models were also used for organ at risk (OAR) in BMO plans, which include the LKB and Poisson models. Next, BMO plans were compared with their corresponding DVO plans, in order to evaluate BMO plan quality. ISO plans were also compared with DVO and BMO plans, in order to verify the performance of the integrated strategy.


BMO plans produced slightly inhomogeneity and less coverage of planning target volume (PTV) (V95=96.79, HI = 0.10: p < 0.01). However, the tumor control probability (TCP) value, both from DVO and BMO plans, were comparable. For the OARs, BMO plans produced lower normal tissue complication probability (NTCP) of rectum (NTCP = 0.11) and bladder (NTCP = 0.14) than in the corresponding DVO plans (NTCP = 0.19 and 0.18 for rectum and bladder; p < 0.01 for rectum and p = 0.03 for bladder). V95, D98, CI and HI values that were produced by ISO plans (V95 = 98.31, D98 = 54.18Gy, CI = 0.76, HI = 0.09) were greatly better than BMO plans (V95 = 96.79, D98 = 53.42Gy, CI = 0.71, HI = 0.10) with significant differences. Furthermore, ISO plans produced lower NTCP values of rectum (NTCP = 0.14) and bladder (NTCP = 0.16) than DVO plans (NTCP = 0.19 and 0.18 for rectum and bladder, respectively) with significant differences.


BMO plans produced lower NTCP values of OARs compared to DVO plans for cervical carcinoma cases, and resulted in slightly less target coverage and homogeneity. The integrated strategy, proposed in this study, could improve the coverage, conformity and homogeneity of PTV greater than the BMO plans, as well as reduce the NTCP values of OARs greater than the DVO plans.


Intensity-modulated radiation therapy (IMRT) is one of the most common and effective radiotherapy treatment techniques for cervical carcinoma [13]. In conventional IMRT plan designing, physical constraints (dose volume parameters) are usually implemented for optimization and plan evaluation. Recently, biological constraints, such as generalized equivalent uniform dose (gEUD), tumor control probability (TCP) and normal tissue complication probability (NTCP), have been used in IMRT planning process [4]. Biological parameters, which applied in biological models to predict the biological effect of tumor and normal tissues, have more direct correlation with treatment outcome than physical parameters [510]. Research shows that, with similar target coverage, biological optimization could spare more organs at risk (OARs) than physical optimization [1118]. However, there are also some disadvantages in biological optimization; such as highly inhomogeneous target dose distributions [11, 12, 18]. Thus, additional physical constraints involved in the biological optimization have been proposed in several studies in order to improve the target dose homogeneity and conformity [4, 12, 1518]. Kan et. al indicated that the amount of hot spots and the maximum doses within targets could be substantially reduced by adding the physical constraints in biological optimization [18].

The biological planning module in Eclipse System (Varian Medical Systems, Palo Alto, CA) could be used for biological optimization with biological constraints like EUD, TCP and NTCP, as well as physical constraints. Kan et. al found that the target coverage and conformity, produced by biological planning module of Eclipse system, were comparable with physical optimization, but spared more parotid glands [18]. Additionally, the research also outlined the limitations of the biological planning module and stated that it would produce more hot spots in PTV for the current module, for it does not allow the users to specify the priority of physical constraints [18]. In the study, additional physical constraints were implemented in the biological planning module to assistant the biological optimization. There are no further studies on how to effectively integrate the biological constraints and physical constraints in the biological planning module for producing better plan quality.

In this study, the primary purpose is to investigate whether the biological optimization could be an alternative method to the conventional physical optimization for cervical carcinoma, by comparing plan qualities between the dose volume based optimization (DVO, only physical constraints were used in the optimization) and biological model based optimization (BMO, biological constraints were used in the optimization with several assistant physical constraints). Furthermore, the IMRT plan quality produced by BMO were tested and evaluated by different biological models. Based on both dosimetrical and biological endpoints consideration, we propose a novel strategy, integrated strategy based optimization (ISO, physical and biological constraints were motivationally integrated in the optimization) to improve plan quality. And all IMRT plans were evaluated by physical and biological indices.


Patient selection

Twenty patients with cervical carcinoma were selected for this study, which pathology diagnosis ranged from stage IIIB to IV. The computed tomography (CT) images (3 mm slice thickness, 512 × 512 pixels/slice) were used for acquiring volumetric anatomical data and analyzing dosimetry comparison. This study was a retrospective study, and patient treatment outcomes were not involved.

Treatment planning

The Eclipse treatment system (Version 13.5, Varian Medical System Inc, Palo Alto, CA) was used to generate all IMRT plans optimized by biological and physical constraints in this study. Specifically, 6 MV photon beams of a linear accelerator (Trilogy, Varian Medical Systems) with 120 multi-leaf collimator was applied.

The clinical target volume (CTV) of each patient, including the primary tumor area, uterus, and the pelvic and para-aortic lymph nodes, was contoured by an experienced oncologist. The corresponding planning target volume (PTV) was generated by expanding 0.5 cm from CTV symmetrically. The OARs included rectum, bladder, femur-heads and small bowel in this study. In order to improve the target dose conformity, the assistance organ Body-PTV (B-P) was defined as the body volume in the CT data set minus the PTV leaving a 0.3 cm gap. B-P was used in all IMRT optimization in order to standardize the optimized constraints.

The prescribed dose to PTV was set to 54 Gy with 30 fractions. All plans were normalized to a mean dose of PTV (and isodose 95% was set to the prescribed dose) in order to make plan comparisons valid. For all patients, the IMRT plans were created using 7 evenly distributed coplanar fields (every 51 °). The same photon beam settings of DVO, including number and orientations of beams, beam energy and iso-center position, were used for corresponding BMO and ISO. This is to eliminate the differences of plan quality that result from the variation in radiation beam parameters.

The optimization process was applied by two-step algorithm in the Eclipse planning system [19, 20]. All dose calculations were performed with the anisotropic analytical algorithm (AAA) with a calculation grid of 2.5 mm [21]. All plans were generated by a singular, experienced dosimetrist in order to avoid the variation of plan quality.

Dose volume optimization

In DVO plans, the clinical requirements of dose volume for structures were listed in Table 1. According to our institution planning protocol for physical optimization, 95% of PTV should receive at least the prescribed dose, while reducing the dose received by OARs. In this aspect, the higher priorities were to give enough dose coverage to PTV and reduce the maximum doses to B-P; while the dose volume constraints for OARs were given lower priority. The physical constraints were listed in Table 2.

Table 1 Clinical requirements to target and OARs
Table 2 The physical constraints used in DVO

Additionally, in the optimization process, the hot or cold spots within PTV and the high dose region of normal tissue should be re-countered and defined as virtual organs after each optimization cycles. Then, dose volume constrains with appropriate priority were set to these virtual organs, and re-optimization was applied in order to achieve the clinical requirements and produce the desired dose distribution.

Biological model optimization

The BMO plans were generated on biological planning module in Eclipse planning system. In this module, biological models like Poisson-LQ and LKB were involved for biological optimization. The details of these biological models could be found in previous studies [2229]. Additionally, gEUD, a semi-biological constraint, was implemented in this module.

The parameters of the TCP and NTCP functions in the biological planning module, such as the \( {D}_{50} \), n, m, \( \gamma \) and \( \alpha /\beta \) values, could be adjusted manually, and these adjustments would have a significant impact on plan quality. Thus, in this study, we just applied the default parameters, given in the module, which was collected from previous retrospective clinical research [2629].

For qualifying the BMO plans qualities and illuminating the influence caused by different biological models, several biological models (which were pre-set in biological planning module of Eclipse) were implemented. In this research, TCP Poisson-LQ model (cervix-IIIB, IV [26]) was applied for PTV. Different biological models were used for OARs: the rectum used the NTCP Poisson-LQ model (Necrosis/Stenosis [27], as a superscript of “cr” ) and LKB model (Late rectal bleeding, grade > =2 [28], as a superscript of “br”; and Late effects, grade > =3 [29], as a superscript of “ar”), the bladder used the NTCP Poisson-LQ model (Contracture [27], as a superscript of “bb”) and the NTCP LKB model (Late effects, grade > =3 [29], as a superscript of “ab”), and the femur-heads used the NTCP Poisson-LQ model (Necrosis [27]). The pre-set specific values of these biological parameters could be found in Table 3. For comparing different rectum biological models, BMO plans were generated by “ar” “br” and “cr” models respectively with fixed bladder model “ab”. And so did bladder. Fixed the “ar” rectum model, the corresponding BMO plans were generated by “ab” and “bb” bladder models respectively. The small bowel was not considered in the optimization, but the NTCP Poisson-LQ model (Obstruction/Perforation [27]) of small bowel, used for biological evaluation, was also listed in Table 3.

Table 3 The biological constraints and additional physical constrains used in BMO

Due to the similar results of dose distribution from different biological models, only LKB model “ar” for rectum and LKB model “ab” for bladder were implemented in the comparison of DVO, BMO and ISO plans.

Besides the biological constraints, additional physical constraints, such as the Max Dose and Uniformity constraint, were also applied to PTV in BMO plans to reduce hot spots and improve the target dose heterogeneity. A Max Dose constraint was also assigned to the normal tissues represented by B-P in order to improve target conformities. Because the NTCP parameters could not achieve the acceptable doses levels purely, the Max EUD constraints were implemented in optimization process to further reduce the doses to OARs. These additional physical constraints or semi-biological constraints that were used in biological planning module were also listed in Table 3.

Integrated strategy optimization

The biological planning module in the Eclipse not only allowed users to implement physical constraints, but also permitted users to arbitrarily apply either physical or biological constraints for structures. In previous research [4, 12, 1518], physical constraints were only considered as an assistant method to compensate the disadvantage of biological constraints in biological optimization process.

However, based on our comparison of DVO and BMO plans above, a new strategy was proposed: physical constraints were only used for PTV, while biological constraints were only applied for OARs in the biological planning module, to improve the quality of plans generated by the biological planning module. Specifically, in ISO planning, the same physical constraints applied in DVO (constraints for PTV in Table 2) were used for PTV, while the same biological constraints like LKB model “ar” for rectum, LKB model “ab” for bladder, and gEUD constraints applied in BMO (constraints for OARs in Table 3) were used for OARs.

Plan evaluation and statistics analysis

All IMRT plans were evaluated by physical indices and biological indices.

For the dosimetric evaluation, the following indices were reported: V95, the target volume received the prescribe dose as displayed on the cumulative DVH; the maximum and minimum doses, represented by the doses received by 2% (D2) and 98% (D98) of the target volume respectively. Additionally, the homogeneity index (HI) and conformity index (CI) of target were also calculated by the following: \( \mathrm{H}\mathrm{I}=\left({D}_2-{D}_{98}\right)/{D}_{mean} \) [30] and \( \mathrm{C}\mathrm{I}={\left({V}_{PTV* Pre}\right)}^2/\left({V}_{PTV}*{V}_{Pre}\right) \) [31]. For OARs, the mean dose (Dmean) and max dose (Dmax) of rectum and bladder were reported, and the volume of these OARs receiving more than 20, 30, 40 and 50 Gy (V20, V30, V40, V50) were analyzed [32, 33].

The radiobiological evaluation was calculated using the available DVHs. TCP was evaluated for PTV, and the NTCP was used for normal tissues. Specifically, the TCP and NTCP models used in the calculation were listed in Table 3. LKB model “ar” for rectum and LKB model “ab” for bladder were used for plan evaluation in the comparison of DVO, BMO and ISO. While in the comparison of different biological models, the rectum and bladder models was applied respectively.

For statistical analysis, a paired t-test was used to compare the IMRT plans. A p value of p < 0.05 is considered statistically significant. All statistical tests were two-tailed and were performed using R-project software.


Comparison between DVO plans and BMO plans

The physical indices of PTV and OARs extracted from DVO plans and BMO plans were shown in Table 4. For D98 and V95, BMO plans produced lower values (D98 = 53.42 Gy, V95 = 96.79%) with significant differences (p < 0.01) than DVO plans (D98 = 54.57 Gy, V95 = 99.13%). Concerning PTV homogeneity, better results were obtained for DVO plans (HI = 0.08) with significant differences (p < 0.01) in comparison to BMO plans (HI = 0.10). And the same result was observed in the conformity index case. DVO plans produced higher CI value (CI = 0.82) compared to BMO plans (CI = 0.71). Furthermore, there were no regular results between DVO plans and BMO plans in the sparing of rectum and bladder, according to the value of V20, V30, V40 and Dmean. However, Dmax values of rectum and bladder of DVO plans (Dmax = 56.55 Gy and 58.95 Gy for rectum and bladder, respectively) were higher than the corresponding BMO plans (Dmax = 54.29 Gy and 58.22 Gy for rectum and bladder, respectively) with significant differences (p < 0.01 for both rectum and bladder). And DVO plans produced higher V50 values of rectum (V50 = 6.16%) with significant differences (p < 0.01) when compared to BMO plans (V50 = 4.35%).

Table 4 Evaluation of DVO plans, BMO plans and ISO plans

Table 4 also summarized the biological evaluation results for PTV and principal OARs. Considering the TCP value, BMO plans indices were generally higher compared to DVO plans without significant differences. However, the rectum and bladder NTCP values of BMO plans (NTCP = 0.19 and 0.18 for the rectum and bladder, respectively) were lower than DVO plans (NTCP = 0.11 and 0.14 for the rectum and bladder, respectively) with p < 0.01 and p = 0.03 respectively. Additionally, the NTCP values of small bowel were 0 in all plans.

Different biological models

In Table 5, we reported the average NTCP values of OARs optimized by different biological models respectively. The superscripts of NTCP values in the table denote the different biological models of rectum or bladder which listed in Table 3 used in IMRT optimization. In this comparison, when different biological models used to one OAR, we keep the biological models of other OARs fixed. The details of biological models used could be found in the footnote of Table 5. Actually, the different IMRT plans optimized by different NTCP models for one OAR produced very similar dose distribution. But the NTCP values between these biological models shown different behavior. In rectum case, the NTCP value of model ‘ar’ (NTCP = 0.11) was between the model ‘br’ (NTCP = 0.33) and model ‘cr’ (NTCP = 0.01) with significant differences (‘ar’ vs. ‘br’: p < 0.01; ‘ar’ vs. ‘cr’: p < 0.01). In bladder case, the NTCP value of model ‘ab’ (NTCP = 0.14) was greater than that of model ‘bb’ (NTCP = 0.00) with p = 0.03.

Table 5 Evaluation of BMO plans with different biological models

The evaluation of ISO plans

As shown in the Table 4, V95, D98 and CI value produced by ISO plans (V95 = 98.31%, D98 = 54.18 Gy, CI = 0.76) were lower than DVO plans (V95 = 99.13%, D98 = 54.57 Gy, CI = 0.82) with significant differences (p < 0.01 for all indices). ISO plans indices were higher than BMO plans’ (V95 = 96.79%, D98 = 53.42 Gy, CI = 0.71) with significant differences (p < 0.01 for all indices). And these differences were much bigger than the differences between ISO plans and DVO plans. The target homogeneity of ISO plans (HI = 0.09) slightly improved when compared with BMO plans (HI = 0.10), but is slightly worse than DVO plans (HI = 0.08). No significant differences were found in TCP values among these plans.

According to Table 4, ISO plans generate the lowest V20 and Dmean values (V20 = 92.66%, Dmean = 33.49 Gy) of rectum, compared with the other two types plan (V20 = 97.55% and Dmean = 34.62 Gy for DVO plans, respectively; V20 = 95.21% and Dmean = 34.49 Gy for BMO plans, respectively). The same behavior was observed in V30 of bladder.

Additionally, the Dmax value of rectum and bladder generated by ISO plans (Dmax = 56.00 Gy and 58.63 Gy for rectum and bladder, respectively) were generally lower than those generated by DVO plans (Dmax = 56.55 Gy and 58.95 Gy for rectum and bladder, respectively), and higher than those of BMO plans (Dmax = 54.29 Gy and 58.22 Gy for rectum and bladder, respectively) with significant differences (p < 0.01 for all cases). In other words, ISO plans generated an intermediate Dmax value of rectum and bladder between DVO plans and BVO plans. The same situation occurred in the NTCP. Specifically, the NTCP value of rectum produced by ISO plans (NTCP = 0.14) was an intermediate value between DVO plans (NTCP = 0.19) and BMO plans (NTCP = 0.11) with significant differences. Furthermore, the ISO plans produced lower NTCP value (NTCP = 0.16) of bladder compared to the DVO plans (NTCP = 0.18). No significant difference was observed between ISO plans and BMO plans.

The results of small bowel which was not considered in optimization, was also listed in Table 4. There were no significant differences of V20, V45 and NTCP values of small bowel among these three plans. Otherwise, the Dmax value of B-P generated by ISO plans (Dmax = 58.04 Gy) were a median value without significant differences, compared with the one generated by DVO plans (Dmax = 57.59 Gy) and BMO plans (Dmax = 58.39 Gy).

Figure 1 showed the cumulative DVHs, which were averaged over all patients. The PTV curves of these plans were almost overlapped. However, for the rectum, the ISO DVH curve was slightly below the DVO DVH curve and BMO DVH curve in the region from 15Gy to 42Gy. While in the region from 42Gy to 55Gy, the BMO DVH curve was lowest. As for the bladder, the same trend was observed with slightly smaller variance. For B-P, the assistance organ, the ISO DVH curve was lower than the BMO DVH curve, while higher than DVO DVH curve.

Fig. 1
figure 1

The mean DVH averaged over 20 patients of DVO, BMO and ISO plans

Dose distribution at iso-center of one specific patient was shown in Fig. 2, which represents the general condition in all cases. This figure illustrated that the dose distribution produced by BMO plans (Fig. 2.b) had more cold spots (especially in the region the red arrow indicated) in PTV than DVO plans (Fig. 2.a); whereas, in the ISO plans (Fig. 2.c), the cold spots could be reduced as the DVO plans.

Fig. 2
figure 2

Dose distribution of one specific patient. a dose distribution of DVO plan; b dose distribution of BMO plan; c dose distribution of ISO plan


Comparison between DVO plans and BMO plans

Previous studies illustrated that IMRT plans optimized by biological constraints would generally more spare parallel OARs while maintaining the very similar target coverage, compared to those optimized by physical constraints [11, 12, 15, 16, 18]. In this study, compared with DVO, it was reported that using BMO would result in slightly inhomogeneous target volume and lack coverage of PTV, which could be deduced from the HI and \( {V}_{95} \) in Table 4.

The TCP value, both from BMO and DVO plans, were comparable. From the equation of TCP model, it could be deduced that the same prescribed dose (including same numbers of fraction and dose per fraction) was the reason that the comparable TCP value were produced by these two optimizations. But the significant differences were observed in the physical indices, such as V95, D98, CI and HI. This may be due to the TCP model, which was used in the biological planning module, was just effectively maximizing the target-cell killing but not sensitive to the local hot spots, cold spots or the specific dose distribution. Other research studies have produced similar results [11]. On the other hand, since the brachytherapy was not considered in this study, TCP values were slightly smaller (<16) than usual. In other words, the dose received by target was insufficient. The NTCP values of OARs were also very small, which could let the biological constraints slightly effect the optimization.

Some previous studies show that adding certain physical constrains into the biological optimization process might be helpful to improve the target dose homogeneity and conformity [34]. However, Kan et. al reported that the biological planning module produces inferior target dose homogeneity with more hot spots than the DVO plans, even with physical constraints [18]. In our study, two additional physical constraints, namely, the Uniformity constraint and the Max Dose constraint, were applied for PTV in the biological planning module. This was not as effective as hypothesized. It is possible that the priority of constraints could not be modified in the biological planning module of Eclipse system. Similar limitations are discussed in previous research [18].

The NTCP of femur-heads were not referred in the results and discussion section. Because the NCTP values were 0 in all treatment plans.

On the other hand, two major structures were evaluated and compared, including bladder and rectum. Although the advantage of BMO plans was not obvious compared with DVO plans, according to the physical indices (such as V20, V30, V40). BMO plans could produce lower NTCP values than DVO plans (shown in Table 4). Results demonstrated that BMO could more spare rectum and bladder than conventional DVO, similar to previous research of other tumor sites [11, 12, 15, 16].

Based on the comparison of DVO and BMO plans, the min and max dose constraints were sufficient in achieving the clinical requirements in PTV coverage and homogeneity. Compared with BMO plans, these physical constraints had advantage in controlling the dose of local points. Thus, setting the physical constraints for PTV was primarily more efficient and effective. Whereas, for OARs, especially parallel OARs, the setting of physical constraints was more dependent on the planner’s experience and skills. The disadvantage of the constraints was evident when controlling the DVH trend. While biological constraints focusing on the dose received by the whole structure, would be crucial to spare the parallel OARs. Therefore, it is essential for the biological constraints to be the primary consideration for OARs.

Different biological models

Results showed that the different NTCP models for one OAR produce the similar dose distribution. The small NTCP values of OARs (<1), mean the optimization was slightly effected by different biological constraints. Based on this point, only LKB model “ar” for rectum and LKB model “ab” for bladder were implemented in the comparison of DVO, BMO and ISO plans. However, the significant difference of relative NTCP values among different biological models (shown in Table 5) may indicate that different models significantly affect the plan quality from a biological aspect. Actually, it is a challenging work to establish a reliable and consistent biological model. Especially in cervical carcinoma, a combination of fractionated high-dose-rate brachytherapy (HDRB) and external beam radiotherapy (EBRT) should be considered [29]. And for late bladder toxicity after EBRT, there is lack of a clear dose response for whole-bladder radiotherapy [35]. And this is also the main reason that the biological optimization was mostly not implemented in clinical treatment planning. The consistency and reliability of the biological models were not specifically discussed in this work, since our study was just focus on the optimization methods comparison between biological and dose/volume approaches.

The evaluation of ISO plans

In our study, the comparison of DVO plans and BMO plans discussed above showed that the BMO plans produce inferior coverage and homogeneity of target with more cold spots, but also more sparing of OARs. In our previous trail-and-error process, we also evaluated the plan quality that only applied biological constraints to both PTV and OARs. However, the coverage and homogeneity of PTV could not achieve the clinical requirements. Much like the BMO plans, the optimization that involved the additional physical constraints for PTV also did not produced ideal results. As for OARs; by adding physical constraints, biological constraints may be limited to reducing the probability of normal tissue complications.

Based on the above discussion, the novel strategy ISO was proposed and evaluated in this research. This differed from the process in previous studies, which proposed adding certain physical constraints into biological optimization process [4, 12, 1518]. Via ISO strategy, our results of Figs. 1 and 2 and Table 4 show that the ISO generated comparable homogeneity and coverage of PTV as DVO plans, while more sparing OARs NTCP values as BMO plans. However, we must note that the ISO plans produce intermediate values in D98, V95, CI, HI, Dmax and NTCP of rectum, which were directly related with optimized constraints. This demonstrated that the ISO strategy was a tradeoff between physical constraints of PTV and biological constraints of OARs. From another perspective, the least V20 of rectum and V30 of bladder generated by ISO plans indicated that these physical indices were not strongly correlated with the NTCP value. And it could be deduced Dmax, the intermediate value produced by ISO plans, was strongly correlated with NTCP value and should be highly considered in the optimization.

As for the small bowel, no significant differences were observed in all incidences, because the specific constraints for small bowel, or an accurate biology model was not implemented in the optimization process. And this is the main limitation of our work. However, BMO plans produced the highest NTCP value of small bowel, which is contrary to the rectum and bladder. And it is safety to say that the lowest NTCP values of rectum and bladder lead to this highest NTCP value of small bowel in BMO plans. This also could be induced from the similar results of B-P (Fig. 1 and Table 4).

The ISO plans quality was compromised because the priority was not assigned in the biological planning model. The physical optimization, however, was able to modify the priority and dose limits in several repeated optimization cycles [18]. If the priority could be modified for these constraints in ISO, this strategy may produce a better dose distribution and target dose homogeneity. This would be the subject of future studies.

It should be point out that OARs considered in this study were not sufficient. This is contrary to the research by Kan et al. that 7 (5 serial and 2 parallel) OARs were included in both evaluation and optimization process [18]. This is also a limitation of this study. And the ISO approach proposed in this study should be further verified in cases with more OARs and targets.


In this study, BMO plans produce lower NTCP values compared to DVO plans for cervical carcinoma cases, even though this may cause slightly worse target coverage and conformity. DVO plans that apply the min dose and max dose constraints are suitable for PTV, due to its effectively in controlling the local points. On the other hand, BMO plans had an advantage in controlling the entire DVH range for OARs, especially for parallel OARs. We proposed a strategy in this study: set the physical constraints as the PTV primary option and the biological constraints as OAR primary option. By means of this strategy, the indices of PTV in ISO plans improved greater than BMO plans, and the indices of OARs were better than that of DVO plans. The fact that different models significantly affect the plan quality indicates that establishing an accurate biological model of OARs is very necessary and urgent for biological optimization.



Anisotropic analytical algorithm


Computed tomography


Clinical target volume


Generalized equivalent uniform dose


Homogeneity index


Intensity-modulated radiation therapy


Normal tissue complication probability


Organ at risk


Planning target volume


Tumor control probability


  1. Portelance L, Chao KS, Grigsby PW, Bennet H, Low D. Intensity-modulated radiation therapy (IMRT) reduces small bowel, rectum, and bladder doses in patients with cervical cancer receiving pelvic and para-aortic irradiation. Int J Radiat Oncol Biol Phys. 2001;51(1):261–6.

    Article  CAS  PubMed  Google Scholar 

  2. Kavanagh BD, Schefter TE, Wu Q, et al. Clinical application of intensity-modulated radiotherapy for locally advanced cervical cancer. Semin Radiat Oncol. 2002;12(3):260–71.

    Article  PubMed  Google Scholar 

  3. Cihoric N, Tsikkinis A, Tapia C, Aebersold DM, Zlobec I, Lössl K. Dose escalated intensity modulated radiotherapy in the treatment of cervical cancer. Radiat Oncol. 2015;10:240.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Li X, Alber M, Deasy JO, 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.

    Article  Google Scholar 

  5. Lee TF, Fang FM. Quantitative analysis of normal tissue effects in the clinic (QUANTEC) guideline validation using quality of life questionnaire datasets for parotid gland constraints to avoid causing xerostomia during head-and-neck radiotherapy. Radiat Oncol. 2013;106:352–8.

    Article  Google Scholar 

  6. Kierkels RG, Korevaar EW, Steenbakkers RJ, et al. Direct use of multivariable normal tissue complication probability models in treatment plan optimisation for individualised head and neck cancer radiotherapy produces clinically acceptable treatment plans. Radiat Oncol. 2014;112:430–6.

    Article  Google Scholar 

  7. Troeller A, Yan D, Marina O, et al. Comparison and limitations of DVH-based NTCP models derived from 3D-CRT and IMRT data for prediction of gastrointestinal toxicities in prostate cancer patients by using propensity score matched pair analysis. Int J Radiat Oncol Biol Phys. 2015;91:435–43.

    Article  PubMed  Google Scholar 

  8. Lee TF, Yeh SA, Chao PJ, et al. Normal tissue complication probability modeling for cochlea constraints to avoid causing tinnitus after head-and-neck intensity-modulated radiation therapy. Radiat Oncol. 2015;10:194.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Moore KL, Schmidt R, Moiseenko V, et al. Quantifying unnecessary normal tissue complication risks due to suboptimal planning: a secondary study of RTOG 0126. Int J Radiat Oncol Biol Phys. 2015;92:228–35.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Liang X, Penagaricano J, Zheng D, et al. Radiobiological impact of dose calculation algorithms on biologically optimized IMRT lung stereotactic body radiation therapy plans. Radiat Oncol. 2016;11:10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. 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.

    Article  PubMed  Google Scholar 

  12. Wu Q, Djajaputra D, Wu Y, Liu HH, Mohan R. Intensity-modulated radiotherapy optimization with gEUD-guided dose-volume objectives. Phys Med Biol. 2003;48:279–91.

    Article  PubMed  Google Scholar 

  13. Thieke C, Bortfeld T, Niemierko A, Nill S. From physical dose constraints to equivalent uniform dose constraints in inverse radiotherapy planning. Med Phys. 2003;30:2332–9.

    Article  PubMed  Google Scholar 

  14. Stavrev P, Hristov D, Warkentin B, et al. Inverse treatment planning by physically constrained minimization of a biological objective function. Med Phys. 2003;30:2948–58.

    Article  CAS  PubMed  Google Scholar 

  15. Semenenko VA, Reitz B, Day E, et al. Evaluation of a commercial biologically based IMRT treatment planning system. Med Phys. 2008;35:5851–60.

    Article  PubMed  Google Scholar 

  16. Qi XS, Semenenko VA, Li XA. Improved critical structure sparing with biologically-based IMRT optimization. Med Phys. 2009;36:1790–9.

    Article  PubMed  Google Scholar 

  17. Das S. A role for biological optimization within the current treatment planning paradigm. Med Phys. 2009;36:4672–82.

    Article  PubMed  Google Scholar 

  18. Kan MWK, Leung LHT, Yu PKN. The use of biologically related model (Eclipse) for the intensity-modulated radiation therapy planning of nasopharyngeal carcinomas. PLoS One. 2014;9:11.

    Google Scholar 

  19. Ludlum E, Xia P. Comparison of IMRT planning with two-step and one-step optimization: a way to simplify IMRT. Phys Med Biol. 2008;53(3):807–21.

    Article  CAS  PubMed  Google Scholar 

  20. Abate A, Mcbenassi P. Comparison of IMRT planning with two-step and one-step optimization: a strategy for improving therapeutic gain and reducing the integral dose. Phys Med Biol. 2009;54(23):7183–98.

    Article  CAS  PubMed  Google Scholar 

  21. Sievinen J, Ulmer W, Kaissl W. AAA photon dose calculation model in Eclipse. Palo Alto: Varian Medical Systems; 2005. p. 1–18 [RAD #7170B].

    Google Scholar 

  22. Gay HA, Niemierko A. A free program for calculating EUD-based NTCP and TCP in external beam radiotherapy. Phys Med. 2007;23:115–25.

    Article  PubMed  Google Scholar 

  23. Semenenko VA, Li XA. Lyman-Kutcher-Burman NTCP model parameters for radiation pneumonitis and xerostomia based on combined analysis of published clinical data. Phys Med Biol. 2008;53:737–55.

    Article  CAS  PubMed  Google Scholar 

  24. Azzeroni R, Maggio A, Fiorino C, et al. Biological optimization of simultaneous boost on intra-prostatic lesions (DILs): Sensitivity to TCP parameters. Phys Med. 2013;29:592–8.

    Article  CAS  PubMed  Google Scholar 

  25. Zhu J, Zhang ZC, Li BS, et al. Analysis of acute radiation-induced esophagitis in non-small-cell lung cancer patients using the Lyman NTCP model. Radiat Oncol. 2010;97:449–54.

    Article  Google Scholar 

  26. Castro JR, Issa P, Fletcher GH. Carcinoma of the cervix treated by external irradiation alone. Radiology. 1970;95(1):163–6.

    Article  CAS  PubMed  Google Scholar 

  27. Agren-Cronqvist AK, Källman P, Turesson I, Brahme A. Volume and heterogeneity dependence of the dose–response relationship for head and neck tumours. Acta Oncol. 1995;34(6):851–60.

    Article  CAS  PubMed  Google Scholar 

  28. Rancati T, Fiorino C, Gagliardi G, et al. Fitting late rectal bleeding data using different NTCP models: results from an Italian multi-centric study (AIROPROS0101). Radiat Oncol. 2004;73(1):21–32.

    Article  CAS  Google Scholar 

  29. Dale E, Hellebust TP, Skjønsberg A, Høgberg T, Olsen DR. Modeling normal tissue complication probability from repetitive computed tomography scans during fractionated high-dose-rate brachytherapy and external beam radiotherapy of the uterine cervix. Int J Radiat Oncol Biol Phys. 2000;47(4):963–71.

    Article  CAS  PubMed  Google Scholar 

  30. Wu Q, Mohan R, Morris M, et al. Simultaneous integrated boost intensity-modulated radiotherapy for locally advanced head-and-neck squamous cell carcinomas. I: dosimetric results. Int J Radiat Oncol Biol Phys. 2003;56:573–85.

    Article  PubMed  Google Scholar 

  31. van’t Riet A, Mak AC, Moerland MA, et al. A conformation number to quantify the degree of conformality in brachytherapy and external beam irradiation: application to the prostate. Int J Radiat Oncol Biol Phys. 1997;37:731–6.

    Article  Google Scholar 

  32. Stromberger C, Kom Y, Kawgan-Kagan M, et al. Intensity-modulated radiotherapy in patients with cervical cancer. An intra-individual comparison of prone and supine positioning. Radiat Oncol. 2010;5:63.

    Article  PubMed  PubMed Central  Google Scholar 

  33. ZHAI D, YIN Y, GONG G, et al. RapidArc radiotherapy for whole pelvic lymph node in cervical cancer with 6 and 15 MV: a treatment planning comparison with fixed field IMRT. J Radiat Res. 2013;54:166–73.

    Article  PubMed  Google Scholar 

  34. Pow HN, Kwong LW, McMillan A, et al. Xerostomia and quality of life after intensity-modulated radiotherapy vs. conventional radiotherapy for early stage nasopharyngeal carcinoma: Initial report on a randomized controlled clinical trial. Int J Radiat Oncol Biol Phys. 2006;66:981–91.

    Article  PubMed  Google Scholar 

  35. Viswanathan AN, Yorke ED, Marks LB, Eifel PJ, Shipley WU. Radiation dose-volume effects of the urinary bladder. Int J Radiat Oncol Biol Phys. 2010;76(3 Suppl):S116–22.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


We want to thank Dr. Lei Xing and Dr. Yong Yang (Department of Radiation Oncology, Stanford University School of Medicine, Stanford, CA, USA) for the fruitful discussions of biological models and optimization.


This work was jointly supported by National Natural Science Foundation of China (61471226), Natural Science Foundation for Distinguished Young Scholars of Shandong Province (JQ201516), China Postdoctoral Science Foundation (2015 T80739), and research funding from Jinan (201401221). We also thanks to the supporting of Taishan scholar project of Shandong Province.

Availability of data and materials

The data will not be shared.

Authors’ contributions

ZF conceived the idea and DL, ZF, and CT designed the study. ZF and CT collected the data. JZ, JC, GY and SQ performed the analysis. ZF drafted the manuscript with the help of DL, YY, and JZ. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Dengwang Li.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Feng, Z., Tao, C., Zhu, J. et al. An integrated strategy of biological and physical constraints in biological optimization for cervical carcinoma. Radiat Oncol 12, 64 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: