Challenges in radiobiological modeling: can we decide between LQ and LQ-L models based on reviewed clinical NSCLC treatment outcome data?
© Santiago et al. 2016
Received: 12 February 2016
Accepted: 27 April 2016
Published: 6 May 2016
To study the dose-response of stage I non-small-cell lung cancer (NSCLC) in terms of long-term local tumor control (LC) after conventional and hypofractionated photon radiotherapy, modeled with the linear-quadratic (LQ) and linear-quadratic-linear (LQ-L) 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 3-year 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 LQ-L models with α/β = 10 Gy were fitted. Additionally, α/β was estimated from direct LQ fits.
Thirty one studies were included reporting outcome of 2319 patients. The LQ-L fit yielded a significant value of 11.0 ± 5.2 Gy for the dose threshold (Dt) for BED10 at the isocenter. The LQ and LQ-L 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).
Both LQ and LQ-L fits can model local tumor control after conventionally and hypofractionated irradiation and are robust methods for predicting clinical effects. The observed dose-effect for local control in NSCLC is weaker at high doses due to data dispersion. For BED10 values of 100–150 Gy in ≥3 fractions, the differences in isoeffects predicted by both models can be neglected.
KeywordsNon-small cell lung cancer Dose-response modeling Biologically effective dose Linear-quadratic model Alpha-beta ratio
The linear-quadratic (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 second-order 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 . 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) .
Current treatment of choice for stage I non-small 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 3D-conformal 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 . Therefore, SBRT evolved to be the current treatment of choice for early-stage 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 .
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 . 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 scanned-beam carbon ion therapy, radiosensitivity is characterized through α/β values obtained from photon irradiation experiments in vitro  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 LQ-L 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.
investigating the dose-response of NSCLC tumor control data from conventionally fractionated (CF) and stereotactic, hypofractionated radiotherapy treatments (HF), based on a review of published long-term outcome results,
evaluating the validity of the LQ and LQ-L models for both conventional and SBRT treatments,
and obtaining an estimation of the clinical α/β ratio of NSCLC.
Materials and methods
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 Meta-Analyses (PRISMA) statements .
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.
Patients with stage I NSCLC (cT1/2, cN0, cM0) with either central or peripheral tumor location.
Treatment with photon radiotherapy with curative intent either 2D or 3D-conformal radiotherapy or SBRT as single modality treatment. Treatment could be delivered with CyberKnife, GammaKnife or linac-based 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.
Reported outcome of interest (actuarial 3-year local control estimations according to Kaplan-Meier or other methods, in which death was a censored event) with a median follow-up time in each cohort ≥17 months.
From an original study (i.a. prospective randomized controlled trial or prospective or retrospective observational study or case series).
The reported cohorts must include at least 25 patients.
Patients could be judged to be either operable or inoperable if the clinical decision was made in favor of radiotherapy.
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).
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 Kaplan-Meier 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 patient-related 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
The LQ and LQ-L 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 linear-quadratic (LQ) model was applied with α/β fixed to 10 Gy, as it is universally accepted for conventional fractionation. The alternative BED definition derived from the LQ-L model was also applied with α/β equal to 10 Gy, to test if the inclusion of a dose threshold Dt would improve the previous fit. Additionally, a study to tentatively estimate the α/β ratio from these clinical data was carried out.
LQ model with α/β ratio fixed to 10 Gy on the full dataset, and on the HF and CF datasets separately.
LQ-L model with α/β ratio fixed to 10 Gy on the full and on the HF datasets.
LQ model with free γ50, TCD50 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.  the Spearman’s correlation between 3-year 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 LQ-L, 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 .
Selected patient cohorts and description of the studies
Characteristics of included studies with conventionally fractionated treatment regimes. Studies published between 1993 and 2015
No. of pats. with stage T1 - T2
BED10@ isoc [Gy]
BED10@ PTV edge [Gy]
Dose calculation algorithm
Follow-up Median (range) [m]
Kaskowitz 1993 
Jeremic 1997 
1.2 (2× day)
Hayakawa 1999 
no dens corr
Cheung 2002 
Langendijk 2002 
Bradley 2003 
no dens corr
Zehentmayr 2015 
19 (Ia)–21 (Ib)
1.8 (2× day)
Characteristics of included studies with hypofractionated treatment regimes. Studies published between 2003 and 2015
No. of pats. with stage T1 - T2
BED10@PTV edge [Gy]
Dose calculation algorithm
Follow-up median (range) [m]
Onimaru 2003 
dens inhom corr
Xia 2006 
Fritz 2008 
Onimaru 2008 
Baumann 2009 
PB, dens inhom corr
Brown 2009 
Fakiris 2009 
no dens inhom corr
Kopek 2009 
Helax-TMS/ Eclipse, ns
Stephans 2009 
dens inhom corr
Baba 2010 
PB convol with Batho
26 (7–66) (living pats)
Crabtree 2010 
Timmerman 2010 
dens inhom corr
Videtic 2010 
dens inhom corr
Andratschke 2011 
dens inhom corr
Hamamoto 2012 
PB, no dens inhom corr
Lagerwaard 2012 
12 20 7.5
Shibamoto 2012 
PB convo, Batho
Shirata 2012 
PB convolution Batho
Takeda 2012 
Inoue 2013 
dens inhom corr
Takeda 2013 
Hamaji 2015 
PB convol, Batho
Summary of cohort characteristics and clinical follow up for conventionally fractionated and hypofractionated datasets
Total No. pats.
% histology unknown
Median age (range) [y]
Median follow-up (range) [m]
Median BED10@isoc (range) [Gy]
Median No. Pats (range)
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 GTV-to-CTV (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 3D-CT scans in expiration/inspiration and free breathing in 4 cases, on slow CT scans in 5 cases, and on 4D-CT scans in 4 cases. Nine out of 23 studies did not apply any ITV concept. The most frequently used CTV-to-PTV (planning target volume) margins were 0.5 cm in axial and 1 cm in cranio-caudal 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 patient-specific.
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 (BED10), 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 follow-up 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 follow-up 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 follow-up time.
3- and 5- year clinical outcomes for stage I NSCLC
The median value [range] for the 3-year LC was 86 [range: 55–100] % for the HF, and 63 [range: 50–91] %, for the CF series, respectively. Twenty-two out of 26 HF 3-year 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 BED10 of 100 Gy or higher at the isocenter but none of the CF treatments reached a BED10 of 100 Gy.
Spearman’s correlation coefficients between dosimetric parameters (total absorbed dose and BED10) and the different clinical outcomes were calculated. For the complete dataset and the total BED10 at the isocenter, a significant Spearman’s correlation of 0.716 for the local control with BED10 was found. For the BED10 evaluated at the PTV edge, a significant correlation of 0.638 was found between LC and BED10. The total absorbed dose at either dose point however did not correlate with LC.
Modeling local control versus BED
Linear-quadratic (LQ) model with α/β ratios fixed to 10 Gy (two-parameter fit)
Linear-quadratic-linear (LQ-L) model with α/β fixed to 10 Gy (three-parameter fit)
The alternative BED definition including a linear portion in the dose-effect 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 BED10 at the isocenter, the fit based on the full dataset generated a Dt value of 11.0 [8.4–16.7] Gy, together with a TCD50 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 Dt. This fit is shown in Fig. 1b, together with the LQ fit for comparison, and their respective BEDLQ and BEDLQ-L data points.
For the fit based on the HF dataset alone, Dt was equal to 9.8 [5.2–15.0] Gy, TCD50 was 1.0 [na-23.8] Gy, and k equal to 71.4 [54.2–79.6] Gy. Only two data points were found to be below this Dt, 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 BED10 doses at PTV edge were used, similar Dt 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 LQ-L and LQ fits, independently of which BED10 doses were used, at isocenter or PTV edge (see Table 5).
Correlation of local control with BED
Spearman’s correlations of the local control with the biologically effective dose at the isocenter and the planning target volume edge, calculated with different α/β values, for the complete dataset and exclusively for the hypofractionated dataset (all of them significant, p > 0.05)
vs BED @ Isocenter
vs BED @ PTV edge
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 three-parameter logistic fit
Summary of the models for both BED doses calculated at isocenter and PTV edge; all fit parameter values are provided with standard errors (and 68 % CI). This table includes the maximum likelihood ratio tests for comparison between the corresponding LQ and LQ-L models with α/β = 10 Gy
Model concept and dataset
α/β (std error) [Gy]
Dt (std error) [Gy]
TCD50CI 68 % [Gy]
k CI 68 % [Gy]
γ50(std error) [%/%]
Likelihood ratio test LQ vs LQL: Dataset, Df, LogLik, Df, Chisq, Pr (>Chisq)
LQ fixed α/β
CF + HF
α/β = 10
Isocenter, all data 4, 2223.8, 1, 3.42, 0.064
α/β = 10
α/β = 10
Isocenter, HF 4, 1316.4, 1, 1.68, 0.195
LQ-L fixed α/β
CF + HF
α/β = 10
α/β = 10
LQ: free α/β, TCD50 and k
LQ fixed α/β
CF + HF
α/β = 10
PTV Edge, all data 4, 2218.3, 1, 0.614, 0.433
α/β = 10
α/β = 10
PTV Edge, HF 4, 1315.9, 1, 0.009, 0.925
LQ-L fixed α/β
CF + HF
α/β = 10
α/β = 10
LQ: free α/β, TCD50 and k
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 BED10 levels at the isocenter (Fig. 4). Under the assumptions of the LQ model with α/β values of 8, 10, or 15 Gy and the LQ-L model at the isocenter with α/β equal to 10 Gy and a Dt 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.
Review of clinical outcome data after radiotherapy treatment represents the only possibility to gather long-term 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 dose-effect 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 mis-estimation 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 3-year local control in order to depict the dose-response 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 dose-effect 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 TCD50 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 mis-estimation and underdosage at the PTV edge.
We calculated the Spearman’s correlations between outcome parameters and the BED10 evaluated at both dose specification points and observed that the correlations with BED10 at the PTV edge were without exception lower than the corresponding correlations based on BED10 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 .
The BED10 fits based on the CF data and the complete dataset differ in both the TCD50 and k values. This can be explained by the fact that the information required defining the slope of the dose-response curve and the TCD50, 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 TCD50 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 dose-response models of Martel et al.  and Guckenberger et al. . Martel et al. found a TCD50, 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.  and further values in the literature. For instance, Stuschke et al.  reported a value of 1.5, and Okunieff et al.  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 BED10(isocenter) and BED10(PTV edge) are graphically similar to the findings by Guckenberger et al. , which were based on a multicentric compilation of individual patient SBRT data.
The LQ-L 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 state-of-the-art fitting methods to compare the performance of the LQ versus the LQ-L 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 Dt there is a contraction in the BEDLQ-L 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 LQ-L model was also fitted to the HF datasets alone, although we could not obtain a 68 % CI for the TCD50, and the p values of all three parameters were > 0.05. Fig. 2b clearly demonstrates the similarity of our LQ-L fit for the HF data subset to the LQ fit previously presented, therefore TCP predictions will be similar with either model. Although our LQ-L fit based on HF data did not yield significant estimates, the results suggest a Dt estimate in the same range of magnitude of 10 Gy. This fit was tested previously also on hypofractionated data alone by Guckenberger et al. . Their dataset had a median dose per fraction at isocenter of 20.8 Gy with range [6–41] Gy. This group found a Dt 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 TCD50 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. , although no significant p values could be obtained in this case, neither.
The dispersion in the collected data points was high. Due to this fact, in the hypofractionation range the dose-effect relationship appears to be weaker than in other reports [13, 61]. Specifically, for the fraction of HF data above 100 Gy (BED10), the Spearman’s correlation between LC and BED10 is low. It can be speculated that in the region of high tumor control probability and highly hypofractionated treatments the relative contribution of non-radiobiological 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 dose-effect 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 LQ-L model with α/β of 10 Gy, even if these models have contrasting radiobiological implications.
We found a dose-effect relationship in the studied dataset, which in the high BED region was weaker due to considerable dispersion in the data. Both, the LQ and LQ-L models can be fitted to clinical normo- and hypofractionated NSCLC outcome data. The LQ-L model yielded a significant value for the Dt of 11.0 Gy for the model based on BED10(isocenter); however, it produced a comparable TCP fit to the LQ model. For the application of BED10 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 LQ-L 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.
The study was financially supported by a grant of the Anneliese Pohl Stiftung. The authors thank Iuliana Toma-Dasu for her insightful comments about the modeling methods and Marie-Anne Chanrion for her useful comments and careful manuscript proof-reading.
Open AccessThis 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.
- Barendsen GW. Dose fractionation, dose rate and iso-effect relationships for normal tissue responses. Int J Radiat Oncol Biol Phys. 1982;8:1981–97.View ArticlePubMedGoogle Scholar
- Joiner MC, van der Kogel AJ. Basic clinical radiobiology. 4th ed. London: Hodder Arnold; 2009.Google Scholar
- Herrmann T, Baumann M, Dörr W. Klinische Strahlenbiologie -- kurz und bündig. 4th ed. Munich: Elsevier; 2006.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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 non-small cell lung cancer. Strahlenther Onkol. 2014;190:26–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Guckenberger M, Allgauer M, Appold S, Dieckmann K, Ernst I, Ganswindt U, et al. Safety and efficacy of stereotactic body radiotherapy for stage I non-small-cell lung cancer in routine clinical practice: a patterns-of-care and outcome analysis. J Thorac Oncol. 2013;8:1050–8.View ArticlePubMedGoogle Scholar
- Fowler JF. Development of radiobiology for oncology -- a personal view. Phys Med Biol. 2006;51:R263–86.View ArticlePubMedGoogle Scholar
- Kirkpatrick JP, Brenner DJ, Orton CG. Point/counterpoint. the linear-quadratic model is inappropriate to model high dose per fraction effects in radiosurgery. Med Phys. 2009;36:3381–4.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Thames HD, Bentzen SM, Turesson I, Overgaard M, Van den Bogaert W. Time-dose factors in radiotherapy: a review of the human data. Radiother Oncol. 1990;19:219–35.View ArticlePubMedGoogle Scholar
- Stuschke M, Pottgen C. Altered fractionation schemes in radiotherapy. Front Radiat Ther Oncol. 2010;42:150–6.View ArticlePubMedGoogle Scholar
- 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 non-small-cell lung cancer? a meta-analysis. Int J Radiat Oncol Biol Phys. 2011;81:e305–16.View ArticlePubMedGoogle Scholar
- Mehta N, King CR, Agazaryan N, Steinberg M, Hua A, Lee P. Stereotactic body radiation therapy and 3-dimensional conformal radiotherapy for stage I non-small cell lung cancer: a pooled analysis of biological equivalent dose and local control. Pract Radiat Oncol. 2012;2:288–95.View ArticlePubMedGoogle Scholar
- van Baardwijk A, Bosmans G, Bentzen SM, Boersma L, Dekker A, Wanders R, et al. Radiation dose prescription for non-small-cell lung cancer according to normal tissue dose constraints: an in silico clinical trial. Int J Radiat Oncol Biol Phys. 2008;71:1103–10.View ArticlePubMedGoogle Scholar
- 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 non-small cell lung cancer. Biomed Res Int. 2013;2013:391021.PubMedPubMed CentralGoogle Scholar
- Guerrero M, Li XA. Extending the linear-quadratic model for large fraction doses pertinent to stereotactic radiotherapy. Phys Med Biol. 2004;49:4825–35.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Bush DA, Cheek G, Zaheer S, Wallen J, Mirshahidi H, Katerelos A, et al. High-dose hypofractionated proton beam radiation therapy is safe and effective for central and peripheral early-stage non-small cell lung cancer: results of a 12-year experience at Loma Linda University Medical Center. Int J Radiat Oncol Biol Phys. 2013;86:964–8.View ArticlePubMedGoogle Scholar
- Nihei K, Ogino T, Ishikura S, Nishimura H. High-dose proton beam therapy for stage I non-small-cell lung cancer. Int J Radiat Oncol Biol Phys. 2006;65:107–11.View ArticlePubMedGoogle Scholar
- Hata M, Tokuuye K, Kagei K, Sugahara S, Nakayama H, Fukumitsu N, et al. Hypofractionated high-dose proton beam therapy for stage I non-small-cell lung cancer: preliminary results of a phase I/II clinical study. Int J Radiat Oncol Biol Phys. 2007;68:786–93.View ArticlePubMedGoogle Scholar
- Krämer M, Scholz M. Treatment planning for heavy-ion radiotherapy: calculation and optimization of biologically effective dose. Phys Med Biol. 2000;45:3319–30.View ArticlePubMedGoogle Scholar
- Astrahan M. Some implications of linear-quadratic-linear radiation dose-response with regard to hypofractionation. Med Phys. 2008;35:4161–72.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Wang JZ, Huang Z, Lo SS, Yuh WTC, Mayr NA. A generalized linear-quadratic model for radiosurgery, stereotactic body radiation therapy, and high-dose rate brachytherapy. Sci Transl Med. 2010;2:39–48.View ArticleGoogle Scholar
- Moher D, Liberati A, Tetzlaff J, Altman DG. Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. BMJ. 2009;339:b2535.View ArticlePubMedPubMed CentralGoogle Scholar
- Okunieff P, Morgan D, Niemierko A, Suit HD. Radiation dose-response of human tumors. Int J Radiat Oncol Biol Phys. 1995;32:1227–37.View ArticlePubMedGoogle Scholar
- Bentzen SM, Tucker SL. Quantifying the position and steepness of radiation dose-response curves. Int J Radiat Biol. 1997;71:531–42.View ArticlePubMedGoogle Scholar
- Martel MK, Ten Haken RK, Hazuka MB, Kessler ML, Strawderman M, Turrisi AT, et al. Estimation of tumor control probability model parameters from 3-D dose distributions of non-small cell lung cancer patients. Lung Cancer. 1999;24:31–7.View ArticlePubMedGoogle Scholar
- R version 2.15.0, 2012-03-30, Copyright (C) 2012, The R Foundation for Statistical Computing, ISBN 3-900051-07-0.Google Scholar
- Kaskowitz L, Graham MV, Emami B, Halverson KJ, Rush C. Radiation therapy alone for stage I non-small cell lung cancer. Int J Radiat Oncol Biol Phys. 1993;27:517–23.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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 non-small cell lung cancer. Lung Cancer. 1999;26:137–42.View ArticlePubMedGoogle Scholar
- Cheung PCF, Yeung LTF, Basrur V, Ung YC, Balogh J, Danjoux CE. Accelerated hypofractionation for early-stage non-small-cell lung cancer. Int J Radiat Oncol Biol Phys. 2002;54:1014–23.View ArticlePubMedGoogle Scholar
- 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 non-small-cell lung cancer. Int J Radiat Oncol Biol Phys. 2002;53:847–53.View ArticlePubMedGoogle Scholar
- Bradley JD, Wahab S, Lockett MA, Perez CA, Purdy JA. Elective nodal failures are uncommon in medically inoperable patients with stage I non-small-cell lung carcinoma treated with limited radiotherapy fields. Int J Radiat Oncol Biol Phys. 2003;56:342–7.View ArticlePubMedGoogle Scholar
- Bogart JA, Alpert TE, Kilpatrick MC, Keshler BL, Pohar SS, Shah H, et al. Dose-intensive thoracic radiation therapy for patients at high risk with early-stage non-small-cell lung cancer. Clin Lung Cancer. 2005;6:350–4.View ArticlePubMedGoogle Scholar
- Zehentmayr F, Wurstbauer K, Deutschmann H, Fussl C, Kopp P, Dagn K, et al. DART-bid: dose-differentiated accelerated radiation therapy, 1.8 Gy twice daily: high local control in early stage (I/II) non-small-cell lung cancer. Strahlenther Onkol. 2015;191:256–63.View ArticlePubMedGoogle Scholar
- Onimaru R, Shirato H, Shimizu S, Kitamura K, Xu B, Fukumoto S, et al. Tolerance of organs at risk in small-volume, hypofractionated, image-guided radiotherapy for primary and metastatic lung cancers. Int J Radiat Oncol Biol Phys. 2003;56:126–35.View ArticlePubMedGoogle Scholar
- 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 non-small-cell lung cancer. Int J Radiat Oncol Biol Phys. 2006;66:117–25.View ArticlePubMedGoogle Scholar
- Fritz P, Kraus HJ, Blaschke T, MAhlnickel W, Strauch K, Engel-Riedel W, et al. Stereotactic, high single-dose irradiation of stage I non-small cell lung cancer (NSCLC) using four-dimensional CT scans for treatment planning. Lung Cancer. 2008;60:193–9.View ArticlePubMedGoogle Scholar
- Onimaru R, Fujino M, Yamazaki K, Onodera Y, Taguchi H, Katoh N, et al. Steep dose-response relationship for stage I non-small-cell lung cancer using hypofractionated high-dose irradiation by real-time tumor-tracking radiotherapy. Int J Radiat Oncol Biol Phys. 2008;70:374–81.View ArticlePubMedGoogle Scholar
- 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 non-small-cell lung cancer patients treated with stereotactic body radiotherapy. J Clin Oncol. 2009;27:3290–6.View ArticlePubMedGoogle Scholar
- Brown WT, Wu X, Fayad F, Fowler JF, García S, Monterroso MI, et al. Application of robotic stereotactic radiotherapy to peripheral stage I non-small cell lung cancer with curative intent. Clin Oncol (R Coll Radiol). 2009;21:623–31.View ArticleGoogle Scholar
- Fakiris AJ, McGarry RC, Yiannoutsos CT, Papiez L, Williams M, Henderson MA, et al. Stereotactic body radiation therapy for early-stage non-small-cell lung carcinoma: four-year results of a prospective phase II study. Int J Radiat Oncol Biol Phys. 2009;75:677–82.View ArticlePubMedGoogle Scholar
- Kopek N, Paludan M, Petersen J, Hansen AT, Grau C, Hoyer M. Co-morbidity index predicts for mortality after stereotactic body radiotherapy for medically inoperable early-stage non-small cell lung cancer. Radiother Oncol. 2009;93:402–7.View ArticlePubMedGoogle Scholar
- 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 non-small cell lung cancer: the Cleveland clinic experience. J Thorac Oncol. 2009;4:976–82.View ArticlePubMedGoogle Scholar
- Baba F, Shibamoto Y, Ogino H, Murata R, Sugie C, Iwata H, et al. Clinical outcomes of stereotactic body radiotherapy for stage I non-small cell lung cancer using different doses depending on tumor size. Radiat Oncol. 2010;5:81.View ArticlePubMedPubMed CentralGoogle Scholar
- 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 non-small cell lung cancer. J Thorac Cardiovasc Surg. 2010;140:377–86.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Videtic GMM, Stephans K, Reddy C, Gajdos S, Kolar M, Clouser E, et al. Intensity-modulated radiotherapy-based stereotactic body radiotherapy for medically inoperable early-stage lung cancer: excellent local control. Int J Radiat Oncol Biol Phys. 2010;77:344–9.View ArticlePubMedGoogle Scholar
- Andratschke N, Zimmermann F, Boehm E, Schill S, Schoenknecht C, Thamm R, et al. Stereotactic radiotherapy of histologically proven inoperable stage I non-small cell lung cancer: patterns of failure. Radiother Oncol. 2011;101:245–9.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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 non-small cell lung cancer. Int J Radiat Oncol Biol Phys. 2012;83:348–53.View ArticlePubMedGoogle Scholar
- Shibamoto Y, Hashizume C, Baba F, Ayakawa S, Manabe Y, Nagai A, et al. Stereotactic body radiotherapy using a radiobiology-based regimen for stage I non-small cell lung cancer: a multicenter study. Cancer. 2012;118:2078–84.View ArticlePubMedGoogle Scholar
- Shirata Y, Jingu K, Koto M, Kubozono M, Takeda K, Sugawara T, et al. Prognostic factors for local control of stage I non-small cell lung cancer in stereotactic radiotherapy: a retrospective analysis. Radiat Oncol. 2012;7:182.View ArticlePubMedPubMed CentralGoogle Scholar
- 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 non-small-cell lung cancer. Lung Cancer. 2012;77:77–82.View ArticlePubMedGoogle Scholar
- Inoue T, Katoh N, Onimaru R, Shimizu S, Tsuchiya K, Suzuki R, et al. Stereotactic body radiotherapy using gated radiotherapy with real-time tumor-tracking for stage I non-small cell lung cancer. Radiat Oncol. 2013;8:69.View ArticlePubMedPubMed CentralGoogle Scholar
- Takeda A, Sanuki N, Eriguchi T, Kaneko T, Morita S, Handa H, et al. Stereotactic ablative body radiation therapy for octogenarians with non-small cell lung cancer. Int J Radiat Oncol Biol Phys. 2013;86:257–63.View ArticlePubMedGoogle Scholar
- Hamaji M, Chen F, Matsuo Y, Kawaguchi A, Morita S, Ueki N, et al. Video-assisted thoracoscopic lobectomy versus stereotactic radiotherapy for stage I lung cancer. Ann Thorac Surg. 2015;99:1122–9.View ArticlePubMedGoogle Scholar
- 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 early-stage NSCLC treated with stereotactic body radiotherapy. Clin Lung Cancer. 2014;15:287–93.View ArticlePubMedGoogle Scholar
- Guckenberger M, Klement RJ, Allgauer M, Appold S, Dieckmann K, Ernst I, et al. Applicability of the linear-quadratic formalism for modeling local tumor control probability in high dose per fraction stereotactic body radiotherapy for early stage non-small cell lung cancer. Radiother Oncol. 2013;109:13–20.View ArticlePubMedGoogle Scholar