 Research
 Open Access
 Published:
Challenges in radiobiological modeling: can we decide between LQ and LQL models based on reviewed clinical NSCLC treatment outcome data?
Radiation Oncology volume 11, Article number: 67 (2016)
Abstract
Aim
To study the doseresponse of stage I nonsmallcell lung cancer (NSCLC) in terms of longterm local tumor control (LC) after conventional and hypofractionated photon radiotherapy, modeled with the linearquadratic (LQ) and linearquadraticlinear (LQL) approaches and to estimate the clinical α/β ratio within the LQ frame.
Material and methods
We identified studies of curative radiotherapy as single treatment through MedLine search reporting 3year LC as primary outcome of interest. Logistic models coupled with the biologically effective dose (BED) at isocenter and PTV edge according to both the LQ and LQL models with α/β = 10 Gy were fitted. Additionally, α/β was estimated from direct LQ fits.
Results
Thirty one studies were included reporting outcome of 2319 patients. The LQL fit yielded a significant value of 11.0 ± 5.2 Gy for the dose threshold (D_{t}) for BED_{10} at the isocenter. The LQ and LQL fits did not differ substantially. Concerning the estimation of α/β, the value obtained from the direct LQ fit for the complete fractionation range was 3.9 [68 % CI: 2.2–9.0] Gy (p > 0.05).
Conclusion
Both LQ and LQL fits can model local tumor control after conventionally and hypofractionated irradiation and are robust methods for predicting clinical effects. The observed doseeffect for local control in NSCLC is weaker at high doses due to data dispersion. For BED_{10} values of 100–150 Gy in ≥3 fractions, the differences in isoeffects predicted by both models can be neglected.
Introduction
The linearquadratic (LQ) model was developed to describe experimental survival curves of both normal and tumor cells after irradiation. The LQ model fits the cell surviving fraction through a secondorder polynomial on the dose per fraction, with coefficients α and β. The ratio between both coefficients describes the repair capacity of the cells and thus sensitivity to fractionation [1, 2].
The LQ model provides an accurate description of fractionation effects at doses between 1 and 8–10 Gy per fraction [3]. Essentially, this formalism enables isoeffect calculations in current clinical practice, defining the relationships between the biological irradiation effect and key parameters such as dose per fraction, total number of fractions and treatment time. Advancements in this model led to the two most extended, complementary approaches for isoeffect calculation: the biologically effective dose (BED) and the equivalent dose in 2 Gy per fraction (EQD2) [4].
Current treatment of choice for stage I nonsmall cell lung cancer (NSCLC) is surgical tumor extraction. Since it became technically feasible, radiotherapy has been used as an alternative treatment method in inoperable cases. Early approaches used 3Dconformal conventionally fractionated techniques, whereas today, stereotactic body radiotherapy (SBRT) allows for highly precise delivery of radiation, thus enabling hypofractionation to deliver ablative radiation doses in 1–5 fractions [5]. Therefore, SBRT evolved to be the current treatment of choice for earlystage NSCLC in medically inoperable patients and in patients who do not consent to surgery. The high precision of dose delivery facilitates normal tissue sparing, even allowing for dose escalation to potentially improve local control. Despite a growing pool of clinical outcome data, the optimal total dose and fractionation scheme to reach the intended biological effects in terms of both local tumor control and side effects are still under debate [6].
Hypofractionation requires reliable isoeffect calculations. Thus, the relevance of the question has been renewed, whether the description of radiobiological effects based on the LQ formalism is appropriate for hypofractionated treatments [7–9].
If no α/β ratio estimation is available for a specific tumor entity, a generic value of 10 Gy is used for BED calculations, although the precision achievable with such a standard α/β ratio is assumed to be lower. Numerous attempts to calculate clinical α/β values have been made, using available clinical outcome data [10]. Such estimations are specially needed if the α/β of a specific tumor entity is suspected to be lower than 10 Gy. In this case a modified fractionation scheme, which reduces tumor cell recovery between fractions, could increase the therapeutic ratio as is the case e.g. in prostate carcinoma or breast cancer.
Many recent studies aim at outcome review and modeling of the dose response relationship of NSCLC [11–15]. However, studies attempting to estimate the α/β ratio for NSCLC are scarce [11, 15, 16]. It is subject of current debate if the improved outcomes of hypofractionated SBRT are a consequence of an α/β ratio lower than 10 Gy, or even lower than the α/β value of the surrounding normal tissue, which could add a radiobiological rationale to the use of hypofractionation. Alternatively, the improvement could be caused by a reduced repopulation in a shorter overall treatment time.
In addition, in particle radiotherapy, a currently emerging field in radiation oncology, radiobiological considerations are of importance. For proton radiotherapy hypofractionated concepts are aimed for partially as motion management strategy [17–20], so that isoeffect calculations are essential. Apart from isoeffect calculations current treatment planning strategies for light ion therapy also require the attribution of radiobiological properties to both tumor and normal tissues. Specifically, in scannedbeam carbon ion therapy, radiosensitivity is characterized through α/β values obtained from photon irradiation experiments in vitro [21] in one of the mathematical models describing the enhanced biological effect in the Bragg peak, the so called local effect model (LEM), which is implemented in commercial treatment planning systems.
In such situation, clinical long time follow up data is the most valid source of data for modeling approaches, which is however inherently limited by inhomogeneity of treatment parameters and treatment techniques evolving over time. We focused our analysis on the comparison of the LQL versus the LQ model, since most of the mathematical model corrections to the LQ model proposed need additional input parameters [16, 22–24], which are not available for the specific clinical situation.
Therefore, this work aims at:

1)
investigating the doseresponse of NSCLC tumor control data from conventionally fractionated (CF) and stereotactic, hypofractionated radiotherapy treatments (HF), based on a review of published longterm outcome results,

2)
evaluating the validity of the LQ and LQL models for both conventional and SBRT treatments,

3)
and obtaining an estimation of the clinical α/β ratio of NSCLC.
Materials and methods
Study design
We identified inclusion criteria, search strategy, outcome measures of interest and indispensable treatment parameters for the study. The analysis keeps standards of the Preferred Reporting Items for Systematic Reviews and MetaAnalyses (PRISMA) statements [25].
Selection criteria
PubMed was searched (July 2015) without restrictions on the publication date. Abstracts of conference proceedings were excluded and language was restricted to English. Repeated publications based on the same cohort were excluded as were outcome reviews in order to avoid duplicity of cohorts. The selection of studies based on the following criteria was made by two independent researchers.
Study cohorts were eligible only if the following criteria were fulfilled:

1)
Patients with stage I NSCLC (cT1/2, cN0, cM0) with either central or peripheral tumor location.

2)
Treatment with photon radiotherapy with curative intent either 2D or 3Dconformal radiotherapy or SBRT as single modality treatment. Treatment could be delivered with CyberKnife, GammaKnife or linacbased without restrictions on fractionation schemes (e.g. normofractionated, accelerated, hyper or hypofractionated schedules, single fraction irradiation), provided complete information on number of fractions, dose per fraction, total absorbed dose and overall treatment time was available.

3)
Reported outcome of interest (actuarial 3year local control estimations according to KaplanMeier or other methods, in which death was a censored event) with a median followup time in each cohort ≥17 months.

4)
From an original study (i.a. prospective randomized controlled trial or prospective or retrospective observational study or case series).

5)
The reported cohorts must include at least 25 patients.

6)
Patients could be judged to be either operable or inoperable if the clinical decision was made in favor of radiotherapy.

7)
A small amount of heterogeneity in reported dose and patient parameters for a given outcome data point was accepted (e.g. a marginal group of patients treated with a deviating radiation scheme, or a small proportion of tumors included in the report of a stage I population with clinical staging other than I).
Data extraction
The main endpoint of this review was local tumor control (LC) at 3 years. If this value was not explicitly provided in the text, it was extracted from the KaplanMeier diagrams.
The fractionation concept (dose per fraction, number of fractions and total treatment time), and the planning technique were extracted together with further treatment and patientrelated parameters. The range of variation in reported cohort parameters for each local tumor control data point was qualitatively assessed and highly heterogeneous cohorts were excluded. Reported mean or preferably median dose values were used to describe the outcome of each specific patient cohort. All prescription doses were translated to doses at the isocenter and at PTV edge, calculated according to the information provided in each publication (more details in the Additional file 1). Mathematical modeling was performed with both, doses at the isocenter and at the PTV edge. The treatments were classified to be hypofractionated (HF) if 1–10 fractions were delivered with doses per fraction at the isocenter above 6 Gy. Treatments were classified as conventionally fractionated (CF) with a broader definition than in the clinical convention, based on the validity limits of the LQ model according to our current knowledge, namely treatment delivered in more than 10 fractions with fraction doses at isocenter ranging between 1.2 and 6 Gy.
Data analysis and mathematical models
Model parameters were fitted with nonlinear least square optimization methods and confidence intervals were calculated with likelihood profiling. A logistic relationship between tumor control probability (TCP) and the biological effective dose (BED) was assumed, according to the parameterization described in Okunieff et al. [26] and Bentzen et al. [27]. BED was based on the LQ model, calculated from the number of fractions and the dose at the isocenter, taking into account neither repopulation nor hypoxia, according to Eq. 1:
where TCD_{50} is the dose necessary to obtain a local tumor control of 50 % and k is a parameter with dose units that is used to calculate the normalized slope, γ_{50}. This parameter quantifies the change in the expected TCP when a 1 % change in dose occurs, evaluated at the dose level of the TCD_{50}, and represents the maximal slope of the doseresponse relationship. It can be calculated from k and TCD_{50} with the expression [27, 28]:
The same logistic model was implemented with an alternative BED definition, including a transition from linearquadratic dependence for the cell survival to purely linear beyond a certain dose level, the dose threshold D_{t}, as described in [22]:
where n is the number of irradiation fractions, d is the fraction dose, and D_{t} is the threshold for the fraction dose.
The LQ and LQL models were fitted to the joint dataset, and to the CF and HF subsets separately, for BED doses calculated both at the isocenter and the PTV edge. First, the linearquadratic (LQ) model was applied with α/β fixed to 10 Gy, as it is universally accepted for conventional fractionation. The alternative BED definition derived from the LQL model was also applied with α/β equal to 10 Gy, to test if the inclusion of a dose threshold D_{t} would improve the previous fit. Additionally, a study to tentatively estimate the α/β ratio from these clinical data was carried out.
In summary, the following fits were calculated:

1)
LQ model with α/β ratio fixed to 10 Gy on the full dataset, and on the HF and CF datasets separately.

2)
LQL model with α/β ratio fixed to 10 Gy on the full and on the HF datasets.

3)
LQ model with free γ_{50}, TCD_{50} and α/β ratio on the full and the CF datasets.
Finally, in order to compare with the results of the recent analysis of Chi et al. [15] the Spearman’s correlation between 3year LC and BED calculated with different fixed α/β ratios was also investigated on the full and on the HF datasets.
The quality of the fits was assessed with checks of the residuals for normality. Models were ranked according to the Akaike Information Criterion (AIC), and for the LQ and LQL, being nested models, maximum likelihood ratio tests were made. All fitted and calculated values are reported together with their 68 % confidence intervals (CI), whenever possible. Statistical significance was assumed for p values < 0.05. Data handling, statistical analysis, model fitting and graphing were done with the software package R, version 2.15.0 [29].
Results
Selected patient cohorts and description of the studies
In total, 31 studies were identified, which fulfilled the selection criteria, Of those, 8 studies report outcomes after conventionally fractionated treatments of a total of 344 patients [30–37] and 23 studies including 1975 patients reporting on hypofractionated irradiations [38–60]. A total of 34 local control  schedule data points, with doses per fraction ranging from 30 to 1.2 Gy, applied in 1–58 fractions, were collected (see Tables 1 and 2, and Additional file 1 for details of the publication search).
Of all reported tumors, 63.6 % were confirmed to be stage T1, 36.4 % T2 (Table 3). A total of 68.1 % of tumors were histologically confirmed: 45.8 % adenocarcinomas, 34.1 % squamous cell carcinomas, 6.2 % other histologies and 13.9 % carcinoma not otherwise specified (NOS). Of the patients treated with conventional fractionation 86.3% were confirmed medically inoperable, versus 55.2 % of all patients treated with hypofractionated schedules. Median of the reported median ages [age range] was comparable between both groups, namely 72 [range: 35–90] and 75 [range: 29–94] years in the CF and HF groups respectively. Patients, who received conventionally fractionated RT were treated in the time period from 1976 to 2010, whereas patients treated with hypofractionated regimes were irradiated in the time period from 1996 to 2012. In the CF cohort only in one study PETCT was performed for staging in 6 out of 31 patients (Bogart et al. [36]), whereas for many of the HF cohorts PET was a routine procedure; for many of the most recent studies PETstaging was even an inclusion criterion in the retrospective series.
In the 8 series of the CF group, generally a margin of 1–1.5 cm was added around the gross tumor volume (GTV), which was in some cases estimated from port films if no planning computer tomography (CT) scan was available. In the HF series, most frequently no GTVtoCTV (clinical target volume) margins were added, except in 5 out of 23 series. Internal target volume (ITV) concepts were applied in 13/23 studies, based either on addition of GTV from 3DCT scans in expiration/inspiration and free breathing in 4 cases, on slow CT scans in 5 cases, and on 4DCT scans in 4 cases. Nine out of 23 studies did not apply any ITV concept. The most frequently used CTVtoPTV (planning target volume) margins were 0.5 cm in axial and 1 cm in craniocaudal directions. A minimal margin of 0.2 cm was added in one patient cohort treated with the CyberKnife where tumor tracking was used to correct for intrafractional target motion. In 2 out of 23 series the PTV margin definition was patientspecific.
Different dose reporting concepts were found throughout the selected references. Only 5 CF reports out of 8 explicitly mentioned that the dose was prescribed to the isocenter. When not specified, prescription to the isocenter was assumed. In the case of the hypofractionated SBRT data, 9 references reported the prescribed dose to isocenter and 11 to the isodose line encompassing the PTV, which ranged from the 50 to the 100 % isodose, most frequently to the 80 % isodose line. Only one of the SBRT cohorts was treated with IMRT, and in this case the dose was prescribed to the 95 % isodose line enclosing the PTV. More information can be found in the Additional file 1.
Median [min  max] applied BED at the isocenter, calculated with an α/β ratio of 10 Gy (BED_{10}), was 82.0 [67.2–93.5] and 119.8 [76.8–347.5] Gy for the CF and HF groups, respectively. At the PTV edge, values of 77.2 [62.9–87.5] and 89.9 [56.8–211.2] Gy were applied in the CF and HF groups.
Median value of the documented median followup times were 28.5 and 27 months for the CF and HF groups, respectively. Very few of these publications state explicitly the number of patients at risk at each followup time point (four cohorts). In only 5 of the selected studies an estimation of the 95 % CI of the calculated actuarial local control rates was reported and one publication presents the standard error. Therefore, no information could be collected about the precision of the estimated actuarial local control rates, apart from the cohort size and the median followup time.
3 and 5 year clinical outcomes for stage I NSCLC
The median value [range] for the 3year LC was 86 [range: 55–100] % for the HF, and 63 [range: 50–91] %, for the CF series, respectively. Twentytwo out of 26 HF 3year LC data points lie above the 80 % level, while all except two of the CF datasets lie below. All except one HF treatments delivered a BED_{10} of 100 Gy or higher at the isocenter but none of the CF treatments reached a BED_{10} of 100 Gy.
Spearman’s correlation coefficients between dosimetric parameters (total absorbed dose and BED_{10}) and the different clinical outcomes were calculated. For the complete dataset and the total BED_{10} at the isocenter, a significant Spearman’s correlation of 0.716 for the local control with BED_{10} was found. For the BED_{10} evaluated at the PTV edge, a significant correlation of 0.638 was found between LC and BED_{10.} The total absorbed dose at either dose point however did not correlate with LC.
Modeling local control versus BED
Linearquadratic (LQ) model with α/β ratios fixed to 10 Gy (twoparameter fit)
All 3year LC data points for both the HF and CF treatments were fitted to a logistic model coupled with BED_{10} at both isocenter and PTV edge. For the isocenter, the TCD_{50} [68% CI] was 48.3 [23.8–62.4] Gy, k was 44.7 [32.1–64.8] Gy, and the calculated γ_{50} (std error) was 0.27 ± 0.1. A logistic fit under the same assumptions, based on the CF subset alone was made, and resulted in the values: TCD_{50} of 68.9 [50.7–74.4] Gy, k was 20.5 [13.1–50.0] Gy and γ_{50} was 0.84 ± 0.5. The same approach applied to the HF dataset showed a TCD_{50} of −60.2 [−189–3.2] Gy, k of 113.3 [73.4–190.1] Gy, and γ_{50} of −0.13 ± 0.17. These three logistic models are represented in Fig. 1a, and all fit parameter values are summarized in Table 5. The p value was found to be < 0.05 for both TCD_{50} and k simultaneously only for the model based on the full dataset.
The model for the complete dataset based on BED_{10} at PTV edge yielded a TCD_{50} of 28.0 [−0.7–43.1] Gy, and k was equal to 39.7 [28.1–60.5] Gy. These values were 64.2 [48.6–69.3] and 19.5 [12.5–44.6] Gy respectively for the CF, and −19.9 [−93.5–14.5] and 64.4 [42.9–108.1] Gy for the HF group alone. Only the parameter values for the CF dataset yielded significant p values simultaneously. Model summaries and calculated γ_{50} values together with their standard errors obtained with propagation of uncertainties are also presented in Tables 5. These models are represented in Fig. 2a.
Linearquadraticlinear (LQL) model with α/β fixed to 10 Gy (threeparameter fit)
The alternative BED definition including a linear portion in the doseeffect beyond a certain fraction dose was tested, according to Eq. 3, with an α/β ratio fixed to 10 Gy. Fits were performed either based on the full dataset or on the HF subset alone.
For the models based on BED_{10} at the isocenter, the fit based on the full dataset generated a D_{t} value of 11.0 [8.4–16.7] Gy, together with a TCD_{50} value of 44.7 [24.2–58.2] Gy, and k equal to 46.7 [35.1–63.2] Gy. The calculated γ_{50} value was 0.24 ± 0.11. All of the parameter estimates yielded significant p values. 13 out of 34 data points represented a dose per fraction below the estimated D_{t}. This fit is shown in Fig. 1b, together with the LQ fit for comparison, and their respective BED_{LQ} and BED_{LQL} data points.
For the fit based on the HF dataset alone, D_{t} was equal to 9.8 [5.2–15.0] Gy, TCD_{50} was 1.0 [na23.8] Gy, and k equal to 71.4 [54.2–79.6] Gy. Only two data points were found to be below this D_{t}, which was too few to obtain a reliable estimate. Thus, the fit parameters did not produce significant p values. This fit is shown in Fig. 2b, together with the LQ fit for comparison.
When BED_{10} doses at PTV edge were used, similar D_{t} values were found, namely, 12.4 Gy for the complete dataset, and 9.9 Gy for the HF dataset (this fit is shown in Fig. 2b). Likelihood ratio tests showed no difference between LQL and LQ fits, independently of which BED_{10} doses were used, at isocenter or PTV edge (see Table 5).
Correlation of local control with BED
To complete the modeling study, the correlation of the 3year LC with BED under different assumptions for α/β equal to 5, 8.6, 10, 15 and 20 Gy was investigated (see Table 4). For the complete dataset and BED_{10} values at isocenter, a Spearman’s correlation of 0.716 (p < 0.0001) with BED_{10} was found. For all other α/β ratios, correlation values increased marginally from BED_{5} (r = 0.706) to a maximum at BED_{10} and decreased again for BED_{20} (r = 0.706), in all cases being significant. The Spearman’s correlations based on the BED values at the PTV edge decreased with growing α/β values for the complete dataset, from 0.680 to 0.510, all of them being significant and consistently lower than the respective values for the BED_{10} at the isocenter (see Table 4).
In contrast, for the HF subset, the correlation of the LC with the same series of BED_{α/β} at isocenter increased minimally from 0.575 to 0.618, and also the values for BED_{α/β} at PTV edge, from 0.531 to 0.601, all of them being statistically significant.
LQ model with threeparameter logistic fit
Additionally, we attempted to estimate the α/β ratio for NSCLC from the complete data set, fitting all three model parameters simultaneously: TCD_{50}, k and α/β. We obtained for the BED doses at isocenter an α/β value of 3.9 [2.2–9.0] Gy, TCD_{50} of 17.8 [na56.4] Gy, and k of 130.9 [50.1na] Gy, with only TCD_{50} yielding a p value < 0.05. For the CF dataset alone we found a similar α/β of 3.8; however, it was not possible to determine confidence intervals. These fits are represented in Fig. 3b. In order to check plausibility of the fit results and to compare our results with published values we calculated γ_{50} and found a value of 0.0 ± 0.15. For the models based on BED_{10} at PTV edge we found values for alpha/beta of 1.7 [1.3–4.1] Gy for the complete dataset and 4.1 [na] Gy for the CF dataset with no significant p values (Table 5).
Differences in the prediction of isoeffects
In order to translate the impact of the different BED model assumptions on clinical treatment schedules we calculated as an example the doses per fraction which would be necessary to reach selected BED_{10} levels at the isocenter (Fig. 4). Under the assumptions of the LQ model with α/β values of 8, 10, or 15 Gy and the LQL model at the isocenter with α/β equal to 10 Gy and a D_{t} of 11.0 Gy, in order to reach 100 Gy (BED), the maximum differences among models for one and two fractions are 10.5 and 4.5 Gy, respectively. Maximal differences remain below 3.3 Gy for treatments delivered in 3 or more fractions. For 200 Gy (BED), discrepancies between models in fraction size increase to 30.1 and 10.5 Gy for one and two fractions, respectively and remain below 5.6 Gy for 3 fractions and more.
Discussion
Review of clinical outcome data after radiotherapy treatment represents the only possibility to gather longterm information from large numbers of patients, which could serve as basis for statistical analysis for radiobiological modeling. However, this task presents a number of challenges since these datasets are intrinsically heterogeneous. Variability among radiotherapy centers applies to aspects such as target volume definition, dose prescription, planning concepts and delivery techniques with different precision levels. Additional effects hindering precise doseeffect modeling are: the use of multiple and sometimes less appropriate dose calculation algorithms especially in the case of lung tumors (e.g. with limited heterogeneity correction), which may lead to misestimation of the absorbed dose. Additionally, outcome of different tumor stages, dose levels and fractionation schemes are frequently reported together. Furthermore, including historical cohorts implies dealing with changes over time in the standard diagnostic and therapeutic procedures, e.g. for staging, recurrence assessment and radiotherapy image guidance.
To counteract this variability, we applied strict inclusion criteria to the selected publications. The analysis is based on 3year local control in order to depict the doseresponse in mature outcome data, while maintaining a sufficient number of data points. To our knowledge this is one of the largest, most homogeneous patient collectives among similar studies. Through the combination of conventionally fractionated and hypofractionated data, a broad range of doses and fractionation schemes is covered, which is a further requirement in order to achieve conclusive modeling results.
Dosimetric heterogeneity in the PTV can be very pronounced in dose distributions for stereotactic treatments, reaching dose differences between isocenter and PTV edge of up to 50 %. It is not possible to know a priori, which reported dosimetric parameter will describe best the doseeffect relationship: the dose at the PTV edge or the dose at isocenter. Therefore, we calculated models based on both, isocenter and minimum target doses. The TCD_{50} doses estimated for the isocenter doses were in general higher than the ones from the models at PTV edge.
Variability in the estimated doses at PTV edge could arise for instance, from variations in the CTV and PTV margin definitions among institutions, or uncertainties in the dose calculation methods, which in the case of outdated, less accurate dose calculation algorithms for the lung, would produce large dose misestimation and underdosage at the PTV edge.
We calculated the Spearman’s correlations between outcome parameters and the BED_{10} evaluated at both dose specification points and observed that the correlations with BED_{10} at the PTV edge were without exception lower than the corresponding correlations based on BED_{10} at the isocenter. This could be interpreted as an indication of the isocenter doses being more robust than the doses at the PTV edge for retrospective modeling studies, in agreement with previous studies [61].
The BED_{10} fits based on the CF data and the complete dataset differ in both the TCD_{50} and k values. This can be explained by the fact that the information required defining the slope of the doseresponse curve and the TCD_{50}, which together shape the sigmoid region of the logistic function, is provided by the CF data, where generally a lower BED was applied. This explains why TCD_{50} and k are not consistently determined across models, which use different input data: CF + HF, or HF/CF alone. We represented the range of the model input data and the extrapolation regions in all figures, to stress that special care must be taken predicting doses in the extrapolation region.
For comparison, Fig. 1 includes the previously published doseresponse models of Martel et al. [28] and Guckenberger et al. [61]. Martel et al. found a TCD_{50}, which is higher as compared to our findings as they compiled stage I/II data whereas our dataset included stage I tumors only. The values for the slope γ_{50} that we obtained for the complete dataset and the CF subset are lower than the value in Martel et al. [28] and further values in the literature. For instance, Stuschke et al. [11] reported a value of 1.5, and Okunieff et al. [26] found a value of 1.6 on average.
Our fit based on the HF data alone shows a shallower curve than the fits including CF data, since the outcomes after hypofractionation are well above the inflection point of the sigmoidal curve. Our HF fit results for both BED_{10}(isocenter) and BED_{10}(PTV edge) are graphically similar to the findings by Guckenberger et al. [61], which were based on a multicentric compilation of individual patient SBRT data.
The LQL model was proposed to account for the experimental observation that curves of the (log cell survival) versus radiation dose often show a more straightened portion at doses beyond 10 Gy than predicted by the LQ model. This effect is in the majority of cases due to heterogeneity of the radiosensitivity distribution of the cells. This is partly due to differences in the stage of cells in the proliferation cycle, but can also be due to partially hypoxic conditions. Resistant cells tend to survive even at larger doses, which causes the survival curve to become less steep than predicted by the LQ model. This explanation implies that the straightening in the curve is not caused by a fundamental mechanism but by a simple to explain heterogeneity in the distribution of sensitivity. Both models converge at dose of less than about 6 Gy.
Applying stateoftheart fitting methods to compare the performance of the LQ versus the LQL models for different fractionation schemes, we did not find significant differences. Therefore, it was not possible to decide, which model better predicts clinical NSCLC outcome data. For BED corresponding to doses per fraction below 11.0 Gy, the BED points for both concepts overlap, whereas above D_{t} there is a contraction in the BED_{LQL} values, which has no consequence in the fit itself, since it takes place in the region where the TCP approaches 100 %. Clinical consequences of using one model or another are only relevant for highly hypofractionated schedules aiming at delivering BED values well above 100 Gy.
The LQL model was also fitted to the HF datasets alone, although we could not obtain a 68 % CI for the TCD_{50}, and the p values of all three parameters were > 0.05. Fig. 2b clearly demonstrates the similarity of our LQL fit for the HF data subset to the LQ fit previously presented, therefore TCP predictions will be similar with either model. Although our LQL fit based on HF data did not yield significant estimates, the results suggest a D_{t} estimate in the same range of magnitude of 10 Gy. This fit was tested previously also on hypofractionated data alone by Guckenberger et al. [61]. Their dataset had a median dose per fraction at isocenter of 20.8 Gy with range [6–41] Gy. This group found a D_{t} value of 22 Gy with a broad 68 % confidence interval, [14–42] Gy, whose addition to the model did not improve the prediction power.
Estimation of the α/β ratio by fitting three parameters simultaneously on a clinical dataset presenting high dispersion  as is the case of the current work  is challenging. Our LQ model with free α/β did not yield significant values for α/β, nor for TCD_{50} and the slope k, which appears to be too shallow after visual inspection. We concluded that there is no indication for larger α/β values than 10 Gy if the complete range of fractionations is considered. The opposite trend (α/β > 10 Gy) was found for the HF dataset as was also the case in Chi et al. [15], although no significant p values could be obtained in this case, neither.
There are few works aiming to the estimation of a clinical α/β for radiotherapy of NSCLC [10, 11, 15]. Thames et al. [10] published an extremely high α/β value for lung tumors but these authors did not obtain a reliable confidence interval and so, their calculations must be regarded with caution. A similar work was also carried out by Stuschke et al. [11], who found an α/β value of 8.2 Gy. They used a fit similar to ours, but set a constraint to force the model to approach the axes origin by adding a point with low BED and null TCP, with a high fit weight. We also tested this approach (full model information in the Additional file 1), adding a data point at 0 Gy (BED_{α/β}) and 0 % TCP. We observed that this constraint influenced the TCD_{50} value to a small extent only, but could have a strong effect on the slope of the curve, and also on the α/β value, for instance, 3.9 [2.2–9.0] versus 12.6 [10.5–15.0] Gy for the complete dataset and BED_{10}(isocenter). The LQbased fits with α/β of 10 Gy for the CF data alone did not vary much with and without constraint. In contrast, if the fit was based on the HF dataset alone, adding this constraint on the original fit had a pronounced effect on the steepness and TCD_{50} of the curve, which in the constrained fits approached the curves based on CF data (Fig. 4a). For this reason, we think it is preferable not to set a constraint to the model, which largely influences the estimates for k and α/β, and also their standard errors. It seems reasonable to accept that fits done on the HF dataset alone will not reproduce a realistic falloff in LC at low doses nor a plausible, clinical TCD_{50}, since the input doseresponse data are well above that region in this specific case.
The dispersion in the collected data points was high. Due to this fact, in the hypofractionation range the doseeffect relationship appears to be weaker than in other reports [13, 61]. Specifically, for the fraction of HF data above 100 Gy (BED_{10}), the Spearman’s correlation between LC and BED_{10} is low. It can be speculated that in the region of high tumor control probability and highly hypofractionated treatments the relative contribution of nonradiobiological factors to the treatment effect is larger. Such factors are, for instance, subjectivity in target delineation and geographical miss, among others.
The logistic fit was applied also in the study at high dose regions since it is widely used for doseeffect description. However, other functional dependencies of LC with BED might be also appropriate in the high BED range. Concerning the BED concept, it seems likely that more than one approach can fit equivalently well on an inherently noisy dataset like this, for instance the LQ applying a higher α/β ratio, or the LQL model with α/β of 10 Gy, even if these models have contrasting radiobiological implications.
Conclusion
We found a doseeffect relationship in the studied dataset, which in the high BED region was weaker due to considerable dispersion in the data. Both, the LQ and LQL models can be fitted to clinical normo and hypofractionated NSCLC outcome data. The LQL model yielded a significant value for the D_{t} of 11.0 Gy for the model based on BED_{10}(isocenter); however, it produced a comparable TCP fit to the LQ model. For the application of BED_{10} in the range of 100–150 Gy in three fractions or more, the differences in isoeffects predicted by both models can be neglected. Our findings therefore do not allow us to suggest use of the LQL model for an improved fitting compared to the LQ model of local control data in case of hypofractionation. A tentative analysis to establish the optimal α/β ratio in the frame of the LQ model for the full fractionation range did not produce significant estimates, although, it showed a trend for α/β values lower than 10 Gy.
References
 1.
Barendsen GW. Dose fractionation, dose rate and isoeffect relationships for normal tissue responses. Int J Radiat Oncol Biol Phys. 1982;8:1981–97.
 2.
Joiner MC, van der Kogel AJ. Basic clinical radiobiology. 4th ed. London: Hodder Arnold; 2009.
 3.
Herrmann T, Baumann M, Dörr W. Klinische Strahlenbiologie  kurz und bündig. 4th ed. Munich: Elsevier; 2006.
 4.
Bentzen SM, Dörr W, Gahbauer R, Howell RW, Joiner MC, Jones B, et al. Bioeffect modeling and equieffective dose concepts in radiation oncology  terminology, quantities and units. Radiother Oncol. 2012;105:266–8.
 5.
Guckenberger M, Andratschke N, Alheit H, Holy R, Moustakis C, Nestle U, et al. Definition of stereotactic body radiotherapy: principles and practice for the treatment of stage I nonsmall cell lung cancer. Strahlenther Onkol. 2014;190:26–33.
 6.
Guckenberger M, Allgauer M, Appold S, Dieckmann K, Ernst I, Ganswindt U, et al. Safety and efficacy of stereotactic body radiotherapy for stage I nonsmallcell lung cancer in routine clinical practice: a patternsofcare and outcome analysis. J Thorac Oncol. 2013;8:1050–8.
 7.
Fowler JF. Development of radiobiology for oncology  a personal view. Phys Med Biol. 2006;51:R263–86.
 8.
Kirkpatrick JP, Brenner DJ, Orton CG. Point/counterpoint. the linearquadratic model is inappropriate to model high dose per fraction effects in radiosurgery. Med Phys. 2009;36:3381–4.
 9.
Brown JM, Carlson DJ, Brenner DJ. The tumor radiobiology of SRS and SBRT: are more than the 5 Rs involved? Int J Radiat Oncol Biol Phys. 2014;88:254–62.
 10.
Thames HD, Bentzen SM, Turesson I, Overgaard M, Van den Bogaert W. Timedose factors in radiotherapy: a review of the human data. Radiother Oncol. 1990;19:219–35.
 11.
Stuschke M, Pottgen C. Altered fractionation schemes in radiotherapy. Front Radiat Ther Oncol. 2010;42:150–6.
 12.
Zhang J, Yang F, Li B, Li H, Liu J, Huang W, et al. Which is the optimal biologically effective dose of stereotactic body radiotherapy for stage I nonsmallcell lung cancer? a metaanalysis. Int J Radiat Oncol Biol Phys. 2011;81:e305–16.
 13.
Mehta N, King CR, Agazaryan N, Steinberg M, Hua A, Lee P. Stereotactic body radiation therapy and 3dimensional conformal radiotherapy for stage I nonsmall cell lung cancer: a pooled analysis of biological equivalent dose and local control. Pract Radiat Oncol. 2012;2:288–95.
 14.
van Baardwijk A, Bosmans G, Bentzen SM, Boersma L, Dekker A, Wanders R, et al. Radiation dose prescription for nonsmallcell lung cancer according to normal tissue dose constraints: an in silico clinical trial. Int J Radiat Oncol Biol Phys. 2008;71:1103–10.
 15.
Chi A, Wen S, Liao Z, Fowler J, Xu J, Nguyen NP, et al. What would be the most appropriate alpha/beta ratio in the setting of stereotactic body radiation therapy for early stage nonsmall cell lung cancer. Biomed Res Int. 2013;2013:391021.
 16.
Guerrero M, Li XA. Extending the linearquadratic model for large fraction doses pertinent to stereotactic radiotherapy. Phys Med Biol. 2004;49:4825–35.
 17.
Laine AM, Pompos A, Timmerman R, Jiang S, Story MD, Pistenmaa D, Choy H. The role of hypofractionated radiation therapy with photons, protons, and heavy ions for treating extracranial lesions. Front Oncol. 2016;5:302.
 18.
Bush DA, Cheek G, Zaheer S, Wallen J, Mirshahidi H, Katerelos A, et al. Highdose hypofractionated proton beam radiation therapy is safe and effective for central and peripheral earlystage nonsmall cell lung cancer: results of a 12year experience at Loma Linda University Medical Center. Int J Radiat Oncol Biol Phys. 2013;86:964–8.
 19.
Nihei K, Ogino T, Ishikura S, Nishimura H. Highdose proton beam therapy for stage I nonsmallcell lung cancer. Int J Radiat Oncol Biol Phys. 2006;65:107–11.
 20.
Hata M, Tokuuye K, Kagei K, Sugahara S, Nakayama H, Fukumitsu N, et al. Hypofractionated highdose proton beam therapy for stage I nonsmallcell lung cancer: preliminary results of a phase I/II clinical study. Int J Radiat Oncol Biol Phys. 2007;68:786–93.
 21.
Krämer M, Scholz M. Treatment planning for heavyion radiotherapy: calculation and optimization of biologically effective dose. Phys Med Biol. 2000;45:3319–30.
 22.
Astrahan M. Some implications of linearquadraticlinear radiation doseresponse with regard to hypofractionation. Med Phys. 2008;35:4161–72.
 23.
Park C, Papiez L, Zhang S, Story M, Timmerman RD. Universal survival curve and single fraction equivalent dose: useful tools in understanding potency of ablative radiotherapy. Int J Radiat Oncol Biol Phys. 2008;70:847–52.
 24.
Wang JZ, Huang Z, Lo SS, Yuh WTC, Mayr NA. A generalized linearquadratic model for radiosurgery, stereotactic body radiation therapy, and highdose rate brachytherapy. Sci Transl Med. 2010;2:39–48.
 25.
Moher D, Liberati A, Tetzlaff J, Altman DG. Preferred reporting items for systematic reviews and metaanalyses: the PRISMA statement. BMJ. 2009;339:b2535.
 26.
Okunieff P, Morgan D, Niemierko A, Suit HD. Radiation doseresponse of human tumors. Int J Radiat Oncol Biol Phys. 1995;32:1227–37.
 27.
Bentzen SM, Tucker SL. Quantifying the position and steepness of radiation doseresponse curves. Int J Radiat Biol. 1997;71:531–42.
 28.
Martel MK, Ten Haken RK, Hazuka MB, Kessler ML, Strawderman M, Turrisi AT, et al. Estimation of tumor control probability model parameters from 3D dose distributions of nonsmall cell lung cancer patients. Lung Cancer. 1999;24:31–7.
 29.
R version 2.15.0, 20120330, Copyright (C) 2012, The R Foundation for Statistical Computing, ISBN 3900051070.
 30.
Kaskowitz L, Graham MV, Emami B, Halverson KJ, Rush C. Radiation therapy alone for stage I nonsmall cell lung cancer. Int J Radiat Oncol Biol Phys. 1993;27:517–23.
 31.
Jeremic B, Shibamoto Y, Acimovic L, Milisavljevic S. Hyperfractionated radiotherapy alone for clinical stage I nonsmall cell lung cancer. Int J Radiat Oncol Biol Phys. 1997;38:521–5.
 32.
Hayakawa K, Mitsuhashi N, Saito Y, Nakayama Y, Furuta M, Sakurai H, et al. Limited field irradiation for medically inoperable patients with peripheral stage I nonsmall cell lung cancer. Lung Cancer. 1999;26:137–42.
 33.
Cheung PCF, Yeung LTF, Basrur V, Ung YC, Balogh J, Danjoux CE. Accelerated hypofractionation for earlystage nonsmallcell lung cancer. Int J Radiat Oncol Biol Phys. 2002;54:1014–23.
 34.
Langendijk JA, Aaronson NK, de Jong JMA, ten Velde GPM, Muller MJ, Slotman BJ, et al. Quality of life after curative radiotherapy in stage I nonsmallcell lung cancer. Int J Radiat Oncol Biol Phys. 2002;53:847–53.
 35.
Bradley JD, Wahab S, Lockett MA, Perez CA, Purdy JA. Elective nodal failures are uncommon in medically inoperable patients with stage I nonsmallcell lung carcinoma treated with limited radiotherapy fields. Int J Radiat Oncol Biol Phys. 2003;56:342–7.
 36.
Bogart JA, Alpert TE, Kilpatrick MC, Keshler BL, Pohar SS, Shah H, et al. Doseintensive thoracic radiation therapy for patients at high risk with earlystage nonsmallcell lung cancer. Clin Lung Cancer. 2005;6:350–4.
 37.
Zehentmayr F, Wurstbauer K, Deutschmann H, Fussl C, Kopp P, Dagn K, et al. DARTbid: dosedifferentiated accelerated radiation therapy, 1.8 Gy twice daily: high local control in early stage (I/II) nonsmallcell lung cancer. Strahlenther Onkol. 2015;191:256–63.
 38.
Onimaru R, Shirato H, Shimizu S, Kitamura K, Xu B, Fukumoto S, et al. Tolerance of organs at risk in smallvolume, hypofractionated, imageguided radiotherapy for primary and metastatic lung cancers. Int J Radiat Oncol Biol Phys. 2003;56:126–35.
 39.
Xia T, Li H, Sun Q, Wang Y, Fan N, Yu Y, et al. Promising clinical outcome of stereotactic body radiation therapy for patients with inoperable stage I/II nonsmallcell lung cancer. Int J Radiat Oncol Biol Phys. 2006;66:117–25.
 40.
Fritz P, Kraus HJ, Blaschke T, MAhlnickel W, Strauch K, EngelRiedel W, et al. Stereotactic, high singledose irradiation of stage I nonsmall cell lung cancer (NSCLC) using fourdimensional CT scans for treatment planning. Lung Cancer. 2008;60:193–9.
 41.
Onimaru R, Fujino M, Yamazaki K, Onodera Y, Taguchi H, Katoh N, et al. Steep doseresponse relationship for stage I nonsmallcell lung cancer using hypofractionated highdose irradiation by realtime tumortracking radiotherapy. Int J Radiat Oncol Biol Phys. 2008;70:374–81.
 42.
Baumann P, Nyman J, Hoyer M, Wennberg B, Gagliardi G, Lax I, et al. Outcome in a prospective phase II trial of medically inoperable stage I nonsmallcell lung cancer patients treated with stereotactic body radiotherapy. J Clin Oncol. 2009;27:3290–6.
 43.
Brown WT, Wu X, Fayad F, Fowler JF, García S, Monterroso MI, et al. Application of robotic stereotactic radiotherapy to peripheral stage I nonsmall cell lung cancer with curative intent. Clin Oncol (R Coll Radiol). 2009;21:623–31.
 44.
Fakiris AJ, McGarry RC, Yiannoutsos CT, Papiez L, Williams M, Henderson MA, et al. Stereotactic body radiation therapy for earlystage nonsmallcell lung carcinoma: fouryear results of a prospective phase II study. Int J Radiat Oncol Biol Phys. 2009;75:677–82.
 45.
Kopek N, Paludan M, Petersen J, Hansen AT, Grau C, Hoyer M. Comorbidity index predicts for mortality after stereotactic body radiotherapy for medically inoperable earlystage nonsmall cell lung cancer. Radiother Oncol. 2009;93:402–7.
 46.
Stephans KL, Djemil T, Reddy CA, Gajdos SM, Kolar M, Mason D, et al. A comparison of two stereotactic body radiation fractionation schedules for medically inoperable stage I nonsmall cell lung cancer: the Cleveland clinic experience. J Thorac Oncol. 2009;4:976–82.
 47.
Baba F, Shibamoto Y, Ogino H, Murata R, Sugie C, Iwata H, et al. Clinical outcomes of stereotactic body radiotherapy for stage I nonsmall cell lung cancer using different doses depending on tumor size. Radiat Oncol. 2010;5:81.
 48.
Crabtree TD, Denlinger CE, Meyers BF, El Naqa I, Zoole J, Krupnick AS, et al. Stereotactic body radiation therapy versus surgical resection for stage I nonsmall cell lung cancer. J Thorac Cardiovasc Surg. 2010;140:377–86.
 49.
Timmerman R, Paulus R, Galvin J, Michalski J, Straube W, Bradley J, et al. Stereotactic body radiation therapy for inoperable early stage lung cancer. JAMA. 2010;303:1070–6.
 50.
Videtic GMM, Stephans K, Reddy C, Gajdos S, Kolar M, Clouser E, et al. Intensitymodulated radiotherapybased stereotactic body radiotherapy for medically inoperable earlystage lung cancer: excellent local control. Int J Radiat Oncol Biol Phys. 2010;77:344–9.
 51.
Andratschke N, Zimmermann F, Boehm E, Schill S, Schoenknecht C, Thamm R, et al. Stereotactic radiotherapy of histologically proven inoperable stage I nonsmall cell lung cancer: patterns of failure. Radiother Oncol. 2011;101:245–9.
 52.
Hamamoto Y, Kataoka M, Yamashita M, Nogami N, Sugawara Y, Kozuki T, et al. Factors affecting the local control of stereotactic body radiotherapy for lung tumors including primary lung cancer and metastatic lung tumors. Jpn J Radiol. 2012;30:430–4.
 53.
Lagerwaard FJ, Verstegen NE, Haasbeek CJA, Slotman BJ, Paul MA, Smit EF, et al. Outcomes of stereotactic ablative radiotherapy in patients with potentially operable stage I nonsmall cell lung cancer. Int J Radiat Oncol Biol Phys. 2012;83:348–53.
 54.
Shibamoto Y, Hashizume C, Baba F, Ayakawa S, Manabe Y, Nagai A, et al. Stereotactic body radiotherapy using a radiobiologybased regimen for stage I nonsmall cell lung cancer: a multicenter study. Cancer. 2012;118:2078–84.
 55.
Shirata Y, Jingu K, Koto M, Kubozono M, Takeda K, Sugawara T, et al. Prognostic factors for local control of stage I nonsmall cell lung cancer in stereotactic radiotherapy: a retrospective analysis. Radiat Oncol. 2012;7:182.
 56.
Takeda A, Kunieda E, Sanuki N, Aoki Y, Oku Y, Handa H. Stereotactic body radiotherapy (SBRT) for solitary pulmonary nodules clinically diagnosed as lung cancer with no pathological confirmation: comparison with nonsmallcell lung cancer. Lung Cancer. 2012;77:77–82.
 57.
Inoue T, Katoh N, Onimaru R, Shimizu S, Tsuchiya K, Suzuki R, et al. Stereotactic body radiotherapy using gated radiotherapy with realtime tumortracking for stage I nonsmall cell lung cancer. Radiat Oncol. 2013;8:69.
 58.
Takeda A, Sanuki N, Eriguchi T, Kaneko T, Morita S, Handa H, et al. Stereotactic ablative body radiation therapy for octogenarians with nonsmall cell lung cancer. Int J Radiat Oncol Biol Phys. 2013;86:257–63.
 59.
Hamaji M, Chen F, Matsuo Y, Kawaguchi A, Morita S, Ueki N, et al. Videoassisted thoracoscopic lobectomy versus stereotactic radiotherapy for stage I lung cancer. Ann Thorac Surg. 2015;99:1122–9.
 60.
Rwigema JC, Chen AM, Wang PC, Lee JM, Garon E, Lee P. Incidental mediastinal dose does not explain low mediastinal node recurrence rates in patients with earlystage NSCLC treated with stereotactic body radiotherapy. Clin Lung Cancer. 2014;15:287–93.
 61.
Guckenberger M, Klement RJ, Allgauer M, Appold S, Dieckmann K, Ernst I, et al. Applicability of the linearquadratic formalism for modeling local tumor control probability in high dose per fraction stereotactic body radiotherapy for early stage nonsmall cell lung cancer. Radiother Oncol. 2013;109:13–20.
Acknowledgements
The study was financially supported by a grant of the Anneliese Pohl Stiftung. The authors thank Iuliana TomaDasu for her insightful comments about the modeling methods and MarieAnne Chanrion for her useful comments and careful manuscript proofreading.
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors report no conflicts of interest. The authors alone are responsible for the content and writing of the paper.
Authors’ contributions
AW, SB, UJ and RE conceived the study. AS designed the search strategy, collected and assembled the data, analyzed and interpreted the results and drafted the manuscript. AW and SB contributed to define the search strategy and supported the data collection. UJ also contributed to define the search strategy. AW and UJ revised the manuscript and made substantial contributions to the interpretation of the results. All authors read and approved the final manuscript.
Additional file
Additional file 1:
Details of the Pubmed search and of the studies included in this work. (PDF 197 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Santiago, A., Barczyk, S., Jelen, U. et al. Challenges in radiobiological modeling: can we decide between LQ and LQL models based on reviewed clinical NSCLC treatment outcome data?. Radiat Oncol 11, 67 (2016). https://doi.org/10.1186/s1301401606435
Received:
Accepted:
Published:
Keywords
 Nonsmall cell lung cancer
 Doseresponse modeling
 Biologically effective dose
 Linearquadratic model
 Alphabeta ratio