Skip to main content
Download PDF

Clinical implications in the use of the PBC algorithm versus the AAA by comparison of different NTCP models/parameters

Radiation Oncology volume 8, Article number: 164 (2013) Cite this article

Abstract

Purpose

Retrospective analysis of 3D clinical treatment plans to investigate qualitative, possible, clinical consequences of the use of PBC versus AAA.

Methods

The 3D dose distributions of 80 treatment plans at four different tumour sites, produced using PBC algorithm, were recalculated using AAA and the same number of monitor units provided by PBC and clinically delivered to each patient; the consequences of the difference on the dose-effect relations for normal tissue injury were studied by comparing different NTCP model/parameters extracted from a review of published studies. In this study the AAA dose calculation is considered as benchmark data. The paired Student t-test was used for statistical comparison of all results obtained from the use of the two algorithms.

Results

In the prostate plans, the AAA predicted lower NTCP value (NTCPAAA) for the risk of late rectal bleeding for each of the seven combinations of NTCP parameters, the maximum mean decrease was 2.2%. In the head-and-neck treatments, each combination of parameters used for the risk of xerostemia from irradiation of the parotid glands involved lower NTCPAAA, that varied from 12.8% (sd=3.0%) to 57.5% (sd=4.0%), while when the PBC algorithm was used the NTCPPBC’s ranging was from 15.2% (sd=2.7%) to 63.8% (sd=3.8%), according the combination of parameters used; the differences were statistically significant. Also NTCPAAA regarding the risk of radiation pneumonitis in the lung treatments was found to be lower than NTCPPBC for each of the eight sets of NTCP parameters; the maximum mean decrease was 4.5%. A mean increase of 4.3% was found when the NTCPAAA was calculated by the parameters evaluated from dose distribution calculated by a convolution-superposition (CS) algorithm. A markedly different pattern was observed for the risk relating to the development of pneumonitis following breast treatments: the AAA predicted higher NTCP value. The mean NTCPAAA varied from 0.2% (sd = 0.1%) to 2.1% (sd = 0.3%), while the mean NTCPPBC varied from 0.1% (sd = 0.0%) to 1.8% (sd = 0.2%) depending on the chosen parameters set.

Conclusions

When the original PBC treatment plans were recalculated using AAA with the same number of monitor units provided by PBC, the NTCPAAA was lower than the NTCPPBC, except for the breast treatments. The NTCP is strongly affected by the wide-ranging values of radiobiological parameters.

Background

As a result of the increased sophistication of treatment techniques and delivery methods, the accuracy of highly conformal radiotherapy has improved rapidly with technological advances in recent years. However more demanding modern treatment techniques require better modeling of treatment beams and more sophisticated modeling in the presence of inhomogeneities in order to guarantee accuracy in the calculation of dose distribution. In the clinical routine, calculations of dose to the tumor are performed by commercial treatment planning systems (TPS) and the majority of these systems include dose calculation algorithms with a limited ability to account for the effects of electron transport [1]. The Pencil Beam Convolution (PBC) algorithm is commonly used in clinical practice. However it is well known that it has shortcomings regarding the presence of inhomogeneities, particularly in those regions where charged particle equilibrium does not exist [2–5]. The introduction of convolution-superposition (CS) algorithms that better account for electron transport, have enabled improved calculation of dose distribution, principally in the absence of electronic equilibrium [6–11]. In the Eclipse TPS (Varian Medical Systems) the Anisotropic Analytical Algorithm (AAA) is implemented; it is a 3D pencil-beam kernel-based superposition algorithm [12]. The AAA includes separately modeled contributions from three sources: primary photons, extra-focal photons and contaminating electrons; each of these has an associated fluence, an energy deposition density function and a scatter kernel. A better consideration of inhomogeneities is obtained when the AAA is used. The higher accuracy of the AAA, compared to the Eclipse’s Pencil Beam Convolution (PBC) algorithm is well-established [13–16]. Monte Carlo (MC) simulation is considered to be a gold standard in dose calculation, and it is therefore used to evaluate other dose calculation algorithms. Sterpin et al. [17] investigated the accuracy of the AAA in two studies. First the AAA was compared both with MC and measurements in an inhomogeneous phantom. Second, the AAA and MC were compared with four Intensity Modulated Radiation Therapy (IMRT) treatment plans in the presence of inhomogeneous tissue. They showed good agreement between the AAA and MC and evaluated the improved accuracy of the AAA compared to the PBC algorithm.

Studies Fogliata et al. [18] carried out, show how Collapsed Cone (CC) and AAA manifested a high degree of consistency compared to the MC method, when the impact of photon dose calculation algorithms on expected dose distribution in lungs under different respiratory phases was investigated. PBC proved to be severely defective in calculations, particularly for cases where specific respiratory phases (e.g. deep inspiration breath hold) were assumed for treatment.

In the study of Bragg et al. [19] compared to the PBC algorithm, the AAA was not found to significantly alter the quality of IMRT treatments plans for prostate, parotid or nasopharynx. While its more accurate modeling of lateral electron transport demonstrates significant increases in the volume of PTV being underdosed in IMRT non-small cell lung cancer (NSCLC) treatments plans.

Nielsen et al. [20] investigated the differences in calculated dose distributions and NTCP values between six different dose calculation algorithms for NSCLC treatments. The study showed how the calculated NTCP values for pneumonitis were more sensitive to the choice of algorithm than mean lung dose and V20 which are commonly used for plan evaluation. Furthermore, NTCP for the lungs was calculated using two different model parameter sets; within each dose calculation algorithm, large differences were found between the calculated NTCP values.

The aim of our study was the retrospective analysis of 3D clinical treatment plans to investigate qualitative, possible clinical consequences of the use of PBC versus AAA, which was considered as benchmark data. The 3D dose distributions of 80 treatment plans at four different tumor sites, produced using PBC algorithm, were recalculated using AAA and the same number of monitor units provided by PBC and clinically delivered to each patient. Similarly to the study of Nielsen et al., in addition to the information relating to the normal tissue/target dosimetry and the Tumor Control Probability (TCP), the comparison was performed investigating the consequences on the dose-effect relations for normal tissue injury comparing different Normal Tissue Complication Probability (NTCP) model parameters extracted from a review of published studies.

Several authors have proposed studies to compute NTCP from a survey on clinical tolerance data [21, 22]. After the first parameterization of dose-volume effects reported by Burman et al. [23] and based on the experience from the 2D radiotherapy era, the availability of 3D dose-volume information strongly increased the amount of quantitative data. Consequently much effort was dedicated both to the proposal of reliable dose-volume constraints that demonstrated capability to reduce toxicity, and the development of normal tissue complication probability (NTCP) models properly fitting clinical data. However NTCP models have been re-evaluated using 3D dose data calculated by rather simple computation algorithms; as it is difficult to quantify the clinical consequence of approximate dose calculations, the value of the reported NTCP parameters remains questionable. De Jaeger et al. [24] compared the results of dose calculations in lung tissue using an approximate dose calculation algorithm (the equivalent-pathlength model, EPL) with calculations for a convolution-superposition (CS) algorithm, and the consequences with respect to the estimation of normal lung tissue injury in a group of patients with NSCLC was evaluated [25]. The study demonstrated that when more accurate dose data is available, a re-evaluation of NTCP model parameters is necessary to avoid NTCP being grossly over or underestimated.

The degree to which the NTCP is strongly affected by the wide-ranging values of radiobiological parameters is shown in this study. At present, the evaluation of the optimal values of the radiobiological parameters is difficult; absolute NTCP values are not reliable enough to be considered for evaluating a treatment plan. However the NTCP values have the attractive feature of synthesizing in only one value the whole dose distribution throughout the organ of interest and along with the dose-volume parameters are useful tools for comparing rival plans or for defining dose escalation strategies.

Methods

Patient data, treatment planning and delivery technique

3D clinical treatment plans of 80 patients were reviewed for this study. 20/80 were irradiated for left breast cancer with two tangential fields, 20/80 were treated for non-small cell lung cancer (NSCLC) by three-dimensional conformal radition therapy (3DCRT); as the PTV locations varied widely, the beam angles were adjusted for each individual patient to meet both the set dose-volume constraints for Organs At Risk (OAR) and acceptable dose distributions. The tumor was situated in the middle or lower lobe, eight plans required four fields; all others used five fields. 20/80 patients were irradiated with intensity-modulated radiotherapy for head-and-neck cancer using seven coplanar fields arrangement [26]; finally 20/80 prostate cancer 3DCRT treatments plans (five coplanar field technique 0°,45°,90°,270°,315°) were re-evaluated.

The target volumes were defined in accordance with the 1993 International Commission on Radiation Units and Measurements Report 50 (ICRU Report 50). The gross tumor volumes (GTV) included all known gross disease as determined by imaging and clinical findings. GTVs were expanded to yield corresponding clinical target volumes (CTVs) according to clinical assessment in each case.

For breast cancer, the CTV was glandular breast tissue and the PTV was generated by expanding the CTV by 0.7 cm isotropically, except in the direction of the skin surface. For lung cancer, the stage was IIIA/IIIB (T3) with a broad scope of disease; the GTV was similar to the CTV. The radiation oncologist identified the GTV, and using a margin of 0.3 cm, the CTV was delineated. A margin of 0.7 cm for middle lobe tumors and 1.0 cm for tumors in the lower lobes was added to create the PTV. For head-and-neck cancer, the margins were adjusted to 1.0 cm beyond the GTV to obtain the CTV; the CTV was expanded symmetrically by 0.3 cm in all directions to account for patient setup and motion within the thermoplastic mask. For prostate cancer the CTV was considered to be the prostate plus seminal vesicles; the planning treatment volume (PTV) was obtained by expanding in 3D the CTV by 1.0 cm and 0.7 cm on the prostate–rectum interface to avoid excessive rectal wall involvement.

All patients, except for five head-and- neck cases, were treated with one fraction per day, 5 days a week, with the fraction dose equal to 2Gy at the ICRU reference point [27]. Five head-and-neck treatments receiving 69.96 Gy to PTV1 and 59.40 Gy to PTV2 with simultaneous integrated boost in 33 fractions, the remaining ones received 70.0 Gy, according to clinical risk. Respectively for lung, prostate and breast cancer the prescription dose was 60.0, 76.0 and 50.0 Gy.

The treatment plans were developed using Eclipse 8.1 TPS (Treatment Planning System); the dose distributions of the clinical treatment plans initially performed using the PBC algorithm were recalculated with AAA using the same number of monitor units provided by PBC.

The paired Student t-test was used for statistical comparison of all results obtained from the use of the two algorithms. All tests were two-tailed with a p value of < 0.05 considered the threshold for statistical significance.

For the validation of both the algorithms implemented in the TPS, the tests, the analyses, and the acceptability criteria were in large part based on the report of the AAPM Report 55 [28], and other documents such as the technical report by IAEA [29] were consulted. For the AAA, the outcomes of some tests were comparable to those provided by Van Esch et al. [14].

Dose analysis

For the PTV we evaluated D95%, D2% dose levels on the dose volume histogram (DVH) above which lay 95% and 2% of the volume of the PTV; they were used as a surrogate for dose minimum and dose maximum, respectively. The mean dose to the PTV was also considered.

To describe the degree of the PTV dose inhomogeneity, the Inhomogeneity Index (II) was used and it has been calculated as (D2% − D95%)/D median , II is equal to 0 if no intra-target inhomogeneity is observed.

For breast treatments, D15%, D2% and the mean dose, Dmean, to the heart and homolateral lung were assessed; D2%, D20%, D60%, V20 (relative volume of the lung receiving at least 20.0 Gy) and Dmean to the lung-CTV for NSCLC treatments were considered; mean parotid glands dose for head-and-neck treatment; finally D2%, D50%, D95% and Dmean to the rectum and D2%, D50%, D80% to femoral heads for prostate treatments were recorded.

NTCP and TCP analysis

The NTCP were evaluated by applying two different radiobiological models: the Lyman-Kutcher-Burman (LKB) model [21–23, 30] and the relative seriality (RS) model [31]. For the five head-and-neck simultaneous integrated boost treatments, the linear quadratic equivalent doses at 2 Gy per fraction (EQD2) DVHs were calculated; α∕β = 3 Gy was assumed.

The NTCP and TCP calculation was performed by a home-made software; the dose calculation in Eclipse was calculated using the minimum available grid size, 0.25 cm, and the step size for differential DVH export was chosen to be 25 cGy, it was a compromise between time necessary for TCP/NTCP calculation and the accuracy of the computed values.

Six different combinations of NTCP model/parameters were available for the heart, three quantifying the risk for pericarditis, and three the excess risk for cardiac mortality. For the lung we applied eight different model parameter sets relating to the development of pneumonitis. The lung-GTV is the volume considered from some of the used NTCP parameter sets (e.g. De Jaeger et al. [24]) and using these parameters on the lung-CTV volume would generally be incorrect, as any differences in dose distribution from the PBC or AAA is most noticeable at the lung/tumour interface. However because the tumor was at an advanced stage and the GTV was similar to the CTV, lung-GTV and lung-CTV were only slightly different. Therefore using the NTCP model parameters for lung-GTV on the lung-CTV may be suitable.

For quantifying the risk of xerostomia from irradiation of the parotid glands we used three different combinations of NTCP model/parameters. Seven sets of parameters were used for the prediction of late rectal bleeding in the patients treated for prostate cancer and only one set of parameters relating to the necrosis of femoral heads. The applied parameters are listed in Tables 1, 2, 3, 4, and 5.

Table 1 Summary of NTCP modeling studies of cardiac toxicity
Full size table
Table 2 Summary of NTCP modeling studies of lung toxicity
Full size table
Table 3 Summary of NTCP modeling studies of parotid glands toxicity
Full size table
Table 4 Summary of NTCP modeling studies of rectal toxicity
Full size table
Table 5 NTCP modeling study of femoral heads toxicity
Full size table

Using the LQ model the TCP was also calculated from the DVH of the CTV. For head-and-neck, breast and lung tumor α/β ratio of 10 Gy was used. At present much debate is going on about the value of α/β for prostate cancer. Modeling studies suggest that α/β could be as low as 1.49 Gy, while other studies show higher α/β [32–37]. Because valid data were not available, we decided to do our investigation with two extreme values; a α/β ratio of 1.49 Gy [32] and 10 Gy [37] were chosen, subsequently, these sets of parameters are referred to as α/β_1.49 and α/β_10 respectively.

Results

The results of the comparison of the treatment plans as calculated by two algorithms, PBC and AAA, are summarized in Tables 6, 7, 8, and 9 for breast, lung, head-and-neck and prostate treatment respectively.

Table 6 Summary of differences between treatment plan and radiobiological parameters from the two algorithms for breast treatment
Full size table
Table 7 Summary of differences between treatment plan and radiobiological parameters from the two algorithms for lung treatment
Full size table
Table 8 Summary of differences between treatment plan and radiobiological parameters from the two algorithms for head-and-neck treatment
Full size table
Table 9 Summary of differences between treatment plan and radiobiological parameters from the two algorithms for prostate treatment
Full size table

Subsequently, NTCP calculated with the AAA and PBC algorithm are referred to as NTCPAAA and NTCPPBC, respectively; the NTCP values less than 0.1% are assumed to be zero.

The breast treatment

When AAA was used, the maximum percentage difference was −3.3% for D95% and a increase of 2.0% for II was found (Table 6). The poorer coverage of the PTV was reflected in the TCP, which was significantly lower when the AAA was used, the mean value was 77.3% (sd = 7.7%) and 85.1% (sd = 4.3%) for PBC (p < 0.001). For the ipsilateral lung while mean D2% decreased when the AAA was applied, the mean D15% and Dmean increased by 3.0 Gy and 1.8 Gy respectively; the mean NTCPAAA values were higher than NTCPPBC (Figures 1 and 2). The mean NTCPAAA varied from 0.2% (sd = 0.1%) to 2.1% (sd = 0.3%), while the mean NTCPPBC varied from 0.1% (sd = 0.0%) to 1.8% (sd = 0.2%) depending on the chosen parameters set.

Figure 1
figure 1

Comparison of NTCP for risk of developing pneumonitis following breast treatment computed with the AAA (ordinate) and the PBC algorithm (abscissa) for NTCP models/parameters sets from Table2. Each symbol represents data of an individual patient. The dotted line indicates the line of identity.

Full size image
Figure 2
figure 2

Example of a comparative DVH for a breast plan. The curves calculated by the PBC algorithm are depicted by solid lines and those calculated by the AAA by dotted lines.

Full size image

A mean increase of 0.1% was observed on the NTCPAAA value when it was estimated by De Jeager et al. [24] parameters and re-evaluated using convolution-superposition (CS) algorithm for the dose calculation. When the Seriality Model was applied, the NTCP value estimated by Seppenwoolde et al. [38] was lower than had been predicted by the LKB model.

The incidence of radiation pneumonitis predicted by both Emami et al. [39] / Burman et al. [23] and Emami et al. [39] / Agren-Cronqvist [40] parameters was zero.

Following left breast cancer treatments, the risk of excess cardiac mortality was found to be low; there was no statistically significant difference between the AAA and PBC. The risk for developing pericarditis was zero for all the considered NTCP parameters. These results might be due to little cardiac tissue having been exposed to radiation beam, but on the other hand it is also necessary to consider that presently there are insufficient clinical data on the dose–response characteristics of cardiac tissue on which to base reliable estimates of radiobiological parameters.

The lung treatment

Changes were observed specifically in the PTV coverage; there was a mean percentage difference of about −3% for the D95%, D2% and Dmean when the AAA was used and as a result, a mean decrease of 8.3% for TCP was observed (Table 7).

For the normal lung dose parameters, the AAA predicted a mean reduction of 3.6 Gy for D2%, while the differences for Dmean and D20% were not statistically significant. D60% and V20% were found to be slightly higher when calculated with the AAA than had been predicted by the PBC algorithm. The resultant NTCPAAA values were lower than the NTCPPBC (Figure 3). For both Seppenwoolde et al. [38] and Emami et al. [39] / Burman et al. [23] parameters, the comparison of NTCP values as calculated by two models, LKB and RS, involved higher values if the LKB model was applied.

Figure 3
figure 3

Comparison of NTCP for risk of developing pneumonitis following NSCLC treatment computed with the AAA (ordinate) and the PBC algorithm (abscissa) for NTCP models/parameters sets from Table2. Each symbol represents data of an individual patient. The dotted line indicates the line of identity.

Full size image

When De Jaeger et al. [24] parameters were used we found a mean increase of 4.3%: 14.9% (sd=7.2%) versus 10.6% (sd=4.8%), if the NTCPAAA was calculated with the parameters evaluated from dose distribution based on a CS algorithm.

The head-and-neck treatment

−3.7%, -2.9% and −3.4% were the percentage difference between AAA and PBC for D95%, D2% and Dmean respectively, and the differences were statistically significant (Table 8). A reduction of 4% was observed in the mean TCP with the use of the AAA.

The mean NTCPAAA of parotid gland toxicity values were lower than NTCPPBC. Furthermore the NTCP value obtained using Emami et al. [39]/Burman et al. [23] parameters was lower than the value predicted by the other sets of NTCP parameters. Using the Eisbruck et al. [41] parameters, the risk of a decrease in the salivary flow to 25% of the pre-treatment flow at 1 year post treatment was much higher (see Table 8) than the risk calculated by Roesink et al. [42] parameters which considered the same endpoint.

The prostate treatment

This type of treatment showed a markedly different pattern. The present study showed no clinically significant differences for any of the evaluated PTV dose parameters (Table 9). The type of algorithm did not affect the II (Figure 4). Across the two algorithms, the TCP was within 2.0%. Regarding the rectum, the AAA predicted slightly higher values for the dose parameters (though statistically significant) except D2%. But on the other hand when the PBC was used, low percentage of rectum volume was exposed to a higher dose than had been obtained with the AAA and the NTCPPBC was higher than NTCPAAA (Figure 5).

Figure 4
figure 4

Comparison of inhomogeneity index for treatment plans computed with the AAA (ordinate) and the PBC algorithm (abscissa) for the breast, NSCLC, head-and-neck and prostate treatments. Each symbol represents data of an individual patient. The dotted line indicates the line of identity.

Full size image
Figure 5
figure 5

Example of a comparative DVH for a prostate plan. The curves calculated by the PBC algorithm are depicted by solid lines and those calculated by the AAA by dotted lines.

Full size image

The tissue-bone interfaces were encountered by the lateral fields (90°, 270°); the slight shift towards higher doses of the femoral heads DVH curve was found when the AAA was used (Figure 5). Consequently, the only available combination of parameters (Emami et al. [39] / Burman et al. [23]) predicted a slightly higher risk of necrosis, even if of limited clinical significance.

Discussion

Since the AAA is considered to be a more accurate dose computation algorithm, the comparison between the AAA and PBC dose distributions by the analyzed dose indices, provides an indication of the difference between the dose predicted by the PBC and that actually delivered. In our study the consequences of the difference on the dose-effect relations for normal tissue injury were analyzed, comparing different NTCP model/parameters.

The influence of the low density lung in proximity to the target volume was the factor largely responsible for the differences between the algorithms’ dose distributions for the breast plans. When AAA was used, the increase of II was the result of the broadening of the penumbra which was nicely modeled by the AAA; the TCP was significantly lower.

Using a more accurate algorithm the radiation in the lung was clearly transported further away; when the AAA was used the mean D2% to the lung decreased by 8.1%. This is in accordance with the findings by Knöös et al. [16], which in their comparison of dose calculation algorithms divided the algorithms in two groups based on how changes in electron transport were accounted for; for the breast plans, they found that the D5% to the adjacent lung decreased by 9.5% or more when accurate algorithms were used. While mean D2% decreased when the AAA was applied, the mean D15% and Dmean increased; as result, the mean NTCPAAA values were higher because the available NTCP parameter sets described the lung by a prevalent parallel architecture for the analyzed endpoint. The incidence of radiation pneumonitis predicted by both Emami et al. [39] / Burman et al. [23] and Emami et al. [39] / Agren-Cronqvist [40] parameters was zero. It is worth remembering that the Emami et al. [39] study was derived from the clinical experience of the 2D planning era, without any individually assessed dose-volume information; that most external radiation therapy was delivered with opposing fields and the normal tissue was irradiated with a fairly uniform dose fraction size.

It would be interesting to know the NTCP values for the parameters estimated from Emami’s study. However more reliable parameterizations, fitting radiobiological models to clinical and dosimetric data, have been published and have to be considered to predict the risk of normal tissue toxicity.

In the lung treatment, the extremely inhomogeneous lung region resulted in the major difference of the two algorithms’ abilities to account for inhomogeneities on the final dose distribution. The factor responsible was the widening of the penumbra, in fact looking closer at the isodoses in proximity to the target revealed that they had a larger separation when the more accurate AAA was applied and was worse for those treatments where a large amount of lung tissue is involved in the PTV. These results are in agreement with the report of Bragg et al. [13], who found that lung plans generated the most significant issues in PTV coverage. When the AAA was used, a reduction for D2% was observed, while the differences for Dmean and D20% were not statistically significant, contrary to what was observed in the breast plans. The NTCPAAA values were lower than the NTCPPBC. The highest NTCP values and the maximum differences between NTCPAAA and NTCPPBC were found for Emami et al. [39] parameters. One can see from Table 2 how much these parameters are markedly different than the other parameterizations, whether for the LKB model or for the Seriality model.

The results obtained show how the use of NTCP parameters based on more accurate dose calculation should be recommended to avoid underestimating the calculated values of NTCP. When De Jaeger et al. [24] parameters were used the underestimate was not significant when the lung tissue volume involved was small, such as for breast treatments, but it was relevant when the volume was major such as for lung treatments; we found a mean increase of 4.3%, when the NTCPAAA was calculated with the parameters evaluated from dose distribution based on a CS algorithm.

In the head-and-neck treatments, as every beam passed through a region of low density material the resulting dose distributions of the two algorithms were dissimilar. A reduction of the mean TCP and of the mean NTCP values, quantifying parotid gland toxicity, was observed with the use of the AAA.

Moreover using the Eisbruck et al. [41] parameters, the NTCP value was much higher than the risk calculated by Roesink et al. [42] parameters which considered the same endpoint. This can be taken as a warning to radiation oncologists: before introducing a predictive model into clinical practice, it is necessary to assess if its predictions “make sense” in regard to that clinic’s treatment plans and experience [43, 44].

The dose distribution across the two algorithms was found to be very similar in the prostate plans; no clinically significant differences for any of the evaluated PTV dose parameters were observed. Such small differences between the dose distributions were found because only two of the five beams used for these treatments were such that there was involvement of heterogeneities.

When the PBC was used, a low percentage of rectum volume was exposed to a higher dose than had been obtained with the AAA and the NTCPPBC was higher than NTCPAAA. That is consistent with a prevalent serial architecture of the rectum for the analyzed endpoint. The effect of the variability of the NTCP parameters on the NTCP value is shown in Table 9. The major differences between the two algorithms were for those set of parameters with the n-parameter closer to 0; the highest NTCP values were with the parameters proposed by Tucker et al. [45] and Rancati et al. [46]. It is interesting to note how a slight difference for the m-value and D50-value between Tucker et al. [45] and Söhn et al. [47] parameters resulted in an important difference in the corresponding NTCP values.

The comparison of the two algorithms in the present study is in accordance with the literature; the differences are of minor clinical significance in many situations such as for prostate treatments and probably for other lesions in the pelvic area. The adoption of the AAA into clinical treatment planning practice requires one to fully understand its effect and its potential consequences so as to re-evaluate an assessment of dose-effect relationships and of parameters used in treatment planning decisions [48]. Similarly the introduction of a predictive model into clinical practice has to be prudent as it is necessary to assess if it is based on calculations and treatments similar to those for which the NTCP has to be calculated. The results found in this study show how the NTCP is strongly affected by the wide-ranging values of radiobiological parameters and the differences between the dose distributions of the two tested algorithms yield statistically significant differences in the NTCP values.

Conclusions

In this study, we have tried to investigate qualitative, possible clinical consequences of the use of PBC versus AAA (keeping the same number of monitor units provided by PBC and clinically delivered to each patient) by comparing different NTCP model/parameters. As general result, the NTCPAAA was lower than the NTCPPBC, except for the breast treatments.

The difference in NTCP between PBC and AAA treatment plans could be clinically significant and it may be the subject of a future prospective study.

Moreover we have observed how much the NTCP value depends strongly on the choice of radiobiological parameters. Radiobiological modeling can play an important role in high quality radiotherapy, however uncritical reliance on model results may compromise treatment outcomes and patient safety. It is important to use NTCP parameter sets based on calculations and treatments similar to those for which the NTCP has to be calculated; additionally, it is necessary to improve models and obtain more robust radiobiological parameters [49–56].

Abbreviations

AAA:

Anisotropic analytical algorithm

DVH:

Dose volume histogram

CC:

Collapsed cone

CS:

Convolution-superposition

CTVs:

Clinical target volumes

Dx%:

Dose level on the

DVH:

above which lay x% of the observed volume

3DCRT:

Three-dimensional conformal radition therapy

GTV:

Gross tumor volume

IMRT:

Intensity modulated radiation therapy

LKB:

Lyman-kutcher-burman model

MC:

Monte Carlo

NSCLC:

Non-small cell lung cancer

NTCP:

Normal tissue complication probability

OARs:

Organs at risk

PBC:

Pencil beam convolution algorithm

PTV:

Planning treatment volume

RS:

Relative seriality model

TPS:

Treatment planning system.

References

  1. Arnfield MR, Siantar CH, Siembers J, Garmon P, Cox L, Mohan R: The impact of electron transport on the accuracy of computed dose. Med Phys 2000, 27: 1266-74. 10.1118/1.599004

    Article  CAS  PubMed  Google Scholar 

  2. Carrasco P, Jornet N, Duch MA, Weber L, Ginjaume M, Eudaldo T, Jurado D, Ruiz A, Ribas M: Comparison of dose calculation algorithms in phantoms with lung equivalent heterogeneities under conditions of lateral electronic disequilibrium. Med Phys 2004, 31: 2899-911. 10.1118/1.1788932

    Article  CAS  PubMed  Google Scholar 

  3. Shaine BH, Al-Ghazi MSAL, El-Khatib E: Experimental evaluation of interface doses in the presence of air cavities compared with treatment planning algorithms. Med Phys 1999, 26: 350-5. 10.1118/1.598526

    Article  Google Scholar 

  4. Engelsman M, Damen EMF, Koken PW, Van’t Veld AA, Van Ingen KM, Mijnheer BJ: Impact of simple tissue inhomogneity correction algorithms on conformal radiotherapy of lung tumours. Radiother Oncol 2001, 60: 299-09. 10.1016/S0167-8140(01)00387-5

    Article  CAS  PubMed  Google Scholar 

  5. AAPM: Tissue inhomogeneity corrections for megavoltage photon beams. New York: Report 85; 2004.

    Google Scholar 

  6. Mackie TR, Scrimger JW, Battista JJ: A convolution method of calculating dose for 15 MV x rays. Med Phys 1985, 12: 188-96. 10.1118/1.595774

    Article  CAS  PubMed  Google Scholar 

  7. Boyer AL, Mok EC: A photon dose distribution model employing convolution calculations. Med Phys 1985, 12: 169-77. 10.1118/1.595772

    Article  CAS  PubMed  Google Scholar 

  8. Boyer AL, Mok EC: Calculation of photon dose distribution in an inhomogeneous medium using convolution. Med Phys 1986, 13: 503-09. 10.1118/1.595964

    Article  CAS  PubMed  Google Scholar 

  9. Mohan R, Chui C, Lidofsky L: Differential pencil beam dose computation model for photons. Med Phys 1986, 13: 64-73. 10.1118/1.595924

    Article  CAS  PubMed  Google Scholar 

  10. Woo MK, Cunningham JR: The validity of the density scaling method in primary electron transport for photon and electron beams. Med Phys 1990, 17: 187-94. 10.1118/1.596497

    Article  CAS  PubMed  Google Scholar 

  11. Aspradakis MM, Morrison RH, Richmond ND, Steele A: Experimental verification of convolution/superposition photon dose calculation for radiotherapy treatment planning. Phys Med Biol 2003, 48: 2873-93. 10.1088/0031-9155/48/17/309

    Article  PubMed  Google Scholar 

  12. Ulmer W, Pyyry J, Kaissl W: 3D photon superposition/convolution algorithm and its foundation on results of Monte Carlo calculations. Phys Med Biol 2005, 50: 1767-90. 10.1088/0031-9155/50/8/010

    Article  CAS  PubMed  Google Scholar 

  13. Bragg CM, Wingate K, Conway J: Clinical implications of the anisotropic analytical algorithm for IMRT treatment planning and verification. Radiother Oncol 2008, 86: 276-84. 10.1016/j.radonc.2008.01.011

    Article  PubMed  Google Scholar 

  14. Van Esch A, Tillikainen L, Pyykkonen J, Tenhunen M, Helminen H, Siljamäki S, Alakuijala J, Paiusco M, Lori M, Huyskens DP: Testing of the Analytical Anisotropic Algorithm for photon dose calculation. Med Phys 2006, 33: 4130-48. 10.1118/1.2358333

    Article  PubMed  Google Scholar 

  15. Fogliata A, Nicolini G, Vanetti E, Clivio A, Cozzi L: Dosimetric validation of the anisotropic analytical algorithm for photon dose calculation: fundamental characterization in water. Phys Med Biol 2006, 51: 1421-38. 10.1088/0031-9155/51/6/004

    Article  PubMed  Google Scholar 

  16. Knöös T, Wieslander E, Cozzi L, Brink C, Fogliata A, Albers D, Nyström H, Lassen S: Comparison of dose calculation algorithms for treatment planning in external photon beam therapy for clinical situations. Phys Med Biol 2006, 51: 5785-07. 10.1088/0031-9155/51/22/005

    Article  PubMed  Google Scholar 

  17. Sterpin E, Tomsej M, De Smedt B, Reynaert N, Vynckier S: Monte Carlo evaluation of the AAA treatment planning algorithm in a heterogeneous multilayer phantom and IMRT clinical treatments for an Electa Sl25 Linear Accelerator. Med Phys 2007, 34: 1665-77. 10.1118/1.2727314

    Article  CAS  PubMed  Google Scholar 

  18. Fogliata A, Nicolini G, Vanetti E, Clivio A, Winkler P, Cozzi L: The impact of photon dose calculation algorithms on expected dose dostributions in lungs under different respiratory phases. Phys Med Biol 2008, 53: 2375-90. 10.1088/0031-9155/53/9/011

    Article  PubMed  Google Scholar 

  19. Bragg CM, Conway J: Dosimetric verification of the Anisotropic Analytical Algorithm for radiotherapy treatment planning. Radiother Oncol 2006, 81: 315-23. 10.1016/j.radonc.2006.10.020

    Article  PubMed  Google Scholar 

  20. Nielsen TB, Wieslander E, Fogliata A, Nielsen M, Hansen O, Brink C: Influence of dose calculation algorithms on the predicted dose distributions and NTCP values for NSCLC patients. Med Phys 2011, 38: 2412-18. 10.1118/1.3575418

    Article  PubMed  Google Scholar 

  21. Webb S, Nahum AE: A model for calculating tumor control probability in radiotherapy including the effects of inhomogeneous disributions of dose and clonogenic cell density. Phys Med Biol 1993, 39: 653-66.

    Article  Google Scholar 

  22. Kutcher GJ, Burman C, Brewster L, Goitein M, Mohan R: Histogram reduction method for calculating complication probabilities for three-dimensional treatment planning evaluations. J Radiat Oncol Biol Phys 1991, 21: 137-46.

    Article  CAS  Google Scholar 

  23. Burman C, Kutcher GJ, Emami B, Goiten M: Fitting of normal tissue tolerance data to an analytic function. Int J Radiat Oncol Biol Phys 1991, 21: 123-35.

    Article  CAS  PubMed  Google Scholar 

  24. De Jaeger K, Hoogeman MS, Engelsman M, Seppenwoolde Y, Damen EM, Mijnheer BJ, Boersma LJ, Lebesque JV: Incorporating an improved dose-calculation algorithm in conformal radiotherapy of lung cancer: re-evaluation of dose in normal lung tissue. Radiother Oncol 2003, 69: 1-10. 10.1016/S0167-8140(03)00195-6

    Article  PubMed  Google Scholar 

  25. Belderbos JS, De Jaeger K, Heemsbergen WD, Seppenwoolde Y, Baas P, Boersma LJ, Lebesque JV: First results of a phase I/II dose escalation trial in non-small cell lung cancer using three-dimensional conformal radiotherapy. Radiother Oncol 2003, 66: 119-26. 10.1016/S0167-8140(02)00377-8

    Article  PubMed  Google Scholar 

  26. Martin E, Deville C, Bonnetain F, Bosset M, Créhange G, Truc G, Naudy S, Maingon P: Intensity-modulated radiation therapy in head and neck cancer : prescribed dose, clinical challenges and results. Radiother Oncol 2007, 85: 392-08. 10.1016/j.radonc.2007.10.045

    Article  PubMed  Google Scholar 

  27. ICRU report 62: Prescribing, recording and reporting photon beam therapy. 1999. Supplement to ICRU report 50

    Google Scholar 

  28. AAPM Report 55: Radiation treatment planning dosimetry verification. New York; 1995.

    Google Scholar 

  29. TRS-430: Commissioning and quality assurance of computerized planning systems for radiation treatment of cancer. Vienna: IAEA; 2004.

    Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  31. Kallman P, Agren A, Brahme A: Tumor and normal tissue responses to fractionated non-uniform dose delivery. Int J Radiat Biol 1992, 62: 249-62. 10.1080/09553009214552071

    Article  CAS  PubMed  Google Scholar 

  32. Fowler JF, Chappel R, Ritter M: Is α/β for prostate tumors really low? Int J Radiat Oncol Biol Phys 2001, 50: 1021-31. 10.1016/S0360-3016(01)01607-8

    Article  CAS  PubMed  Google Scholar 

  33. Kal HB, van Gellekom MPR: How low is the α/β ratio for prostate cancer? Int Radiat Oncol Biol Phys 2003, 57: 1116-21. 10.1016/S0360-3016(03)01455-X

    Article  Google Scholar 

  34. Nahum AE, Movsas B, Horwitz EM, Stobbe CC, Chapman JD: Incorporating clinical measurements of hypoxia into tumor local control modelling for prostate cancer: implications for the a/b ratio. Int Radiat Oncol Biol Phys 2003, 57: 391-401. 10.1016/S0360-3016(03)00534-0

    Article  Google Scholar 

  35. Valdagni R, Italia C, Montanaro P, Lanceni A, Lattuada P, Magnani T, Fiorino C, Nahum A: Is the alpha-beta ratio for prostate cancer really low? A prospective, non randomized trial comparing standard and hypofractionated conformal radiation therapy. Radiather Oncol 2005, 75: 74-82. 10.1016/j.radonc.2004.12.019

    Article  Google Scholar 

  36. Fowler JF, Nahum AE, Orton CG: The best radiotherapy for the treatment of prostate cancer involves hypofractionation. Med Phys 2006, 33: 3081-84. 10.1118/1.2179008

    Article  PubMed  Google Scholar 

  37. Nahum AE, Sanchez N: Tumor control probability modelling: basic principles and application in treatment planning. Phys Med 2001, 17: 13-22.

    Google Scholar 

  38. Seppenwoolde Y, Lebesque JV, de Jaeger K, Belderbos JS, Boersma LJ, Schilstra C, Henning GT, Hayman JA, Martel MK, Ten Haken RK: Comparing different NTCP models that predict the incidence of radiation pneumonitis. Int J Radiat Oncol Biol Phys 2003, 55: 724-35. 10.1016/S0360-3016(02)03986-X

    Article  PubMed  Google Scholar 

  39. Emami B, Lyman J, Brown A, Coia L, Goitein M, Munzenrider JE, Shank B, Solin LJ, Wesson M: Tolerance of normal tissue to therapeutic irradiation. Int J Radiat Oncol Biol Phys 1991, 21: 109-12.

    Article  CAS  PubMed  Google Scholar 

  40. Ågren Cronqvist A-K: Quantification of the response of heterogeneous tumors and organized normal tissue to fractionated radiotherapy. Thesis. Stockholm University: Departement of Medical radiation Physics; 1995.

    Google Scholar 

  41. Eisbruch A, Ten Haken RK, Kim HM, Marsh LH, Ship JA: Dose, volume and function relationship in parotid salivary glands following conformal and intensity-modulated irradiation of head and neck cancer. Int J Radiat Oncol Biol Phys 1999, 45: 577-87. 10.1016/S0360-3016(99)00247-3

    Article  CAS  PubMed  Google Scholar 

  42. Roesink JM, Moerland MA, Hoekstra A, Van Rijk PP, Terhaard CH: Scintigraphic assessment of early and late parotid gland function after radiotherapy for head and neck cancer: A prospective study of dose volume response relationship. Int J Radiat Oncol Biol Phys 2004, 58: 1451-60. 10.1016/j.ijrobp.2003.09.021

    Article  PubMed  Google Scholar 

  43. Lawrence BM, Ellen DY, Andrew J, Randall KTH, Louis SC, Avraham E, Søren MB, Jiho N, Joseph OD: Use of normal tissue complication probability models in the clinic. Int Radiat Oncol Biol Phys 2010, 76: S10-S19. 10.1016/j.ijrobp.2009.07.1754

    Article  Google Scholar 

  44. Kong FM, Pan C, Eisbruch A, Haken RKT: Physical Models and simpler dosimetric descriptors of radiation late toxicity. Semin Radiat Oncol 2007, 17: 108-120. 10.1016/j.semradonc.2006.11.007

    Article  PubMed  Google Scholar 

  45. Tucker SL, Dong L, Bosch WR, Michalski J, Winter K, Lee AK, Cheung MR, Kuban DA, Cox JD, Mohan R: Fit of a generalized Lyman normal tissue complication probability (NTCP) model to grade ≥ 2 late rectal toxicity data from patients treated on protocol RTOG 94–06. Int J Radiat Oncol Biol Phys 2007, 69: S8-9.

    Article  Google Scholar 

  46. Rancati T, Fiorino C, Vavassori V, Baccolini M, Bianchi C, Foppiano F, Menegotti L, Monti A, Pasquino M, Stasi M, Fellin G, Valdagni R: Late rectal bleeding after conformal radiotherapy for prostate cancer: NTCP modeling. Radiother Oncol 2008, 88: S332-3.

    Google Scholar 

  47. Söhn M, Yan D, Liang J, Meldolesi E, Vargas C, Alber M: Incidence of late rectal bleeding in high-dose conformal radiotherapy of prostate cancer using equivalent uniform dose-based and dose-volume-based normal tissue complication probability models. Int J Radiat Oncol Biol Phys 2007, 67: 1066-73. 10.1016/j.ijrobp.2006.10.014

    Article  PubMed  Google Scholar 

  48. Morgan AM, Knöös T, McNee SG, Evans CJ, Thwaites DI: Clinical implications of the implementation of advanced treatment planning algorithms for thoracic treatments. Radiother Oncol 2008, 86: 48-54. 10.1016/j.radonc.2007.11.033

    Article  PubMed  Google Scholar 

  49. Schilstra C, Meertens H: Calculation of the uncertainty in complication probability for various dose–response models, applied to the parotid gland. Int J Radiat Oncol Biol Phys 2001, 50: 147-158. 10.1016/S0360-3016(00)01553-4

    Article  CAS  PubMed  Google Scholar 

  50. Brahme A: Optimized radiation therapy based on radiobiological objectives. Semin Radiat Oncol 1999, 9: 35-47. 10.1016/S1053-4296(99)80053-8

    Article  CAS  PubMed  Google Scholar 

  51. Niemierko A, Urie M, Goiten M: Optimization of 3D radiation therapy with both physical and biological end points and constraints. Int J Radiat Oncol Biol Phys 1992, 23: 99-107. 10.1016/0360-3016(92)90548-V

    Article  CAS  PubMed  Google Scholar 

  52. Stewart RD, Li XA: BGRT: Biologically guided radiation therapy – The future is fast approaching! Med Phys 2007, 34: 3739-51. 10.1118/1.2779861

    Article  PubMed  Google Scholar 

  53. Bentzen SM: Radiation therapy: Intensity modulated, image guided, biologically optimized and evidence based. Radiather Oncol 2005, 77: 227-30. 10.1016/j.radonc.2005.11.001

    Article  Google Scholar 

  54. Begg A, van der Kogel A: Clinical radiobiology in 2008. Radiother Oncol 2008, 86: 295-299. 10.1016/j.radonc.2008.02.004

    Article  PubMed  Google Scholar 

  55. Chavaudra N, Bourhis J, Foray N: Quantified relationship between cellular radiosensitivity, DNA repair defects and chromatic relaxation: a study of 19 human tumor cell lines from different origin. Radiother Oncol 2004, 73: 373-82. 10.1016/j.radonc.2004.07.016

    Article  CAS  PubMed  Google Scholar 

  56. Hill RP, Kasper P, Griffin AM, O’Sullivan B, Catton C, Alasti H, Abbas A, Heydarian M, Ferguson P, Wunder JS, Bell RS: Studies of the in vivo radiosensitivity of human skin fibroblasts. Radiother Oncol 2007, 84: 75-83. 10.1016/j.radonc.2007.05.025

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Martel MK, Sahijdak WM, Ten Haken RK, Kessler ML, Tumsi AT: Fraction size and dose parameters related to the incidence of pericardial effusions. Int J Radiat Oncol Biol Phys 1998, 40: 155-61. 10.1016/S0360-3016(97)00584-1

    Article  CAS  PubMed  Google Scholar 

  58. Gagliardi G, Lax I, Ottolenghi A, Rutquist LE: Long-term cardiac mortality after radiotherapy of breast cancer–application of the relative seriality model. Br J Radiol 1996, 69: 839-46. 10.1259/0007-1285-69-825-839

    Article  CAS  PubMed  Google Scholar 

  59. Eriksson F, Gagliardi G, Liedberg A, Lax I, Lee C, Levitt S, Lind B, Rutqvist LE: Long-term cardiac mortality following radiation therapy for Hodgkin’s disease: analysis with the relative seriality model. Radiother Oncol 2000, 55: 153-62. 10.1016/S0167-8140(00)00166-3

    Article  CAS  PubMed  Google Scholar 

  60. Kwa SL, Lebesgue JV, Theuws JCM, Marks LB, Munley MT, Bentel G, Oetzel D, Spahn U, Graham MV, Drzymala RE, Purdy JA, Lichter AS, Martel MK, Ten Haken RK: Radiation pneumonitis as a function of mean dose: an analysis of pooled data of 540 patients. Int J Radiat Oncol Biol Phys 1998, 42: 1-9.

    Article  CAS  PubMed  Google Scholar 

  61. Gagliardi G, Bjöhle J, Lax I, Ottolenghi A, Eriksson F, Liedberg A, Lind P, Rutqvist LE: Radiation pneumonitis after breast cancer irradiation: Analysis of the complication probability using the relative seriality model. Int J Radiat Oncol Biol Phys 2000, 46: 373-81. 10.1016/S0360-3016(99)00420-4

    Article  CAS  PubMed  Google Scholar 

  62. Rancati T, Fiorino C, Gagliardi G, Cattaneo GM, Sanguineti G: Casanova Borca V, Cozzarini C, Fellin G, Foppiano F, Girelli G, Menegotti L, Piazzolla A, Vavassori V, Valdagni R: Fitting late rectal bleeding data using different NTCP models: results from an Italian multi-centric study (AIROPROS0101). Radiother Oncol 2004, 73: 21-32. 10.1016/j.radonc.2004.08.013

    Article  CAS  PubMed  Google Scholar 

  63. Peeters ST, Hoogeman MS, Heemsbergen WD, Hart AA, Koper PC, Lebesque JV: Rectal bleeding, fecal incontinence and high stool frequency after conformal radiotherapy for prostate cancer: normal tissue complication probability modeling. Int J Radiat Oncol Biol Phys 2006, 66: 11-19. 10.1016/j.ijrobp.2006.03.034

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We wish to thank Mr. Thomas Dalton for the English revision of the manuscript.

Author information

Authors and Affiliations

  1. Medical Physics, PioXI Clinic and U.O.C. Medical Physics, S. Giovanni Calibita Fatebenefratelli Hospital, Rome, Italy

    Antonella Bufacchi

  2. UPMC San Pietro FBF Advanced Radiotherapy Center, Rome, Italy

    Barbara Nardiello

  3. U.O.C. Medical Physics, S. Giovanni Calibita Fatebenefratelli Hospital, Rome, Italy

    Roberto Capparella & Luisa Begnozzi

Authors
  1. Antonella Bufacchi
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Barbara Nardiello
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Roberto Capparella
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Luisa Begnozzi
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Antonella Bufacchi.

Additional information

Competing interests

The authors declare they have no competing interests.

Authors’ contributions

AB designed the study, performed the analyses and wrote the manuscript; BN carried out the data from DVH of the treatment plans and helped to perform the analyses; RC helped to make graphs and tables; LB read the different manuscript’s versions. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Bufacchi, A., Nardiello, B., Capparella, R. et al. Clinical implications in the use of the PBC algorithm versus the AAA by comparison of different NTCP models/parameters. Radiat Oncol 8, 164 (2013). https://doi.org/10.1186/1748-717X-8-164

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1748-717X-8-164

Keywords

Radiation Oncology

ISSN: 1748-717X

Contact us