Skip to main content

Calculating the individualized fraction regime in stereotactic body radiotherapy for non-small cell lung cancer based on uncomplicated tumor control probability function

Abstract

Background

To calculate the individualized fraction regime (IFR) in stereotactic body radiotherapy (SBRT) for non-small cell lung cancer (NSCLC) patients using the uncomplicated tumor control probability (UTCP, P+) function.

Methods

Thirty-three patients with peripheral lung cancer or lung metastases who had undergone SBRT were analyzed. Treatment planning was performed using the dose regime of 48 Gy in 4 fractions. Dose volume histogram (DVH) data for the gross tumor volume (GTV), lung, chest wall (CW) and rib were exported and the dose bin was multiplied by a certain percentage of the dose in that bin which ranged from 1 to 200% in steps of 1%. For each dose fraction, P+ values were calculated by considering the tumor control probability (TCP), radiation-induced pneumonitis (RIP), chest wall pain (CWP) and radiation-induced rib fracture (RIRF). UTCP values as a function of physical dose were plotted and the maximum P+ values corresponded to the optimal therapeutic gain. The IFR in 3 fractions was also calculated with the same method by converting the dose using the linear quadratic (LQ) model.

Results

Thirty-three patients attained an IFR using the introduced methods. All the patients achieved a TCP value higher than 92.0%. The IFR ranged from 3 × 10.8 Gy to 3 × 12.5 Gy for 3 fraction regimes and from 4 × 9.2 Gy to 4 × 10.7 Gy for 4 fraction regimes. Four patients with typical tumor characteristics demonstrated that the IFR was patient-specific and could maximize the therapeutic gain. Patients with a large tumor had a lower TCP and UTCP and a smaller fractional dose than patients with a small tumor. Patients with a tumor adjacent to the organ at risk (OAR) or at a high risk of RIP had a lower UTCP and a smaller fractional dose compared with patients with a tumor located distant from the OAR.

Conclusions

The proposed method is capable of predicting the IFR for NSCLC patients undergoing SBRT. Further validation in clinical samples is required.

Background

Stereotactic body radiotherapy (SBRT) has become a standard treatment alternative for patients with medically inoperable early stage non-small cell lung cancer (NSCLC), and for those refusing surgical resection [1,2,3,4]. Recent data have shown that SBRT provides outcomes that are equivalent to those of surgery [5,6,7].

Although SBRT for NSCLC has achieved encouraging outcomes, radiation-induced pneumonitis (RIP), chest wall pain (CWP) and radiation-induced rib fractures (RIRF) are common side effects for NSCLC patients undergoing SBRT. The occurrence of grade ≥ 2 RIP, grade ≥ 2 CWP and symptomatic RIRF ranged between 9.4 and 20.3% [8,9,10,11,12,13,14,15], 10.9 and 39% [16,17,18,19,20] and 12.2 and 17.0% [21,22,23], respectively. Therefore, to develop a method for calculating the individualized fraction regime (IFR) capable of maintaining tumor control probability (TCP) while lowering the risk of normal tissues by considering the tumor size and proximity to the organs at risk (OAR) is a problem to be solved. Recently, two independent studies utilized risk-adapted fraction regimes ranging from 3 to 8 fractions in SBRT treatment for lung cancer and achieved a low incidence of CWP [24, 25]. Unfortunately, the studies considered only the risk of CWP without considering of the tumor size or the toxicities of other OAR.

The current study aimed to develop a method to calculate the patient-specific fraction regime and to maximize the therapeutic gain for peripheral NSCLC patients by incorporating the uncomplicated tumor control probability (UTCP, P+) function.

Methods

Patient eligibility

Computed tomography (CT) simulation data for 33 patients previously diagnosed with primary stage I NSCLC or lung metastases were included in the study. The age of the patients ranged from 51 to 77 years. The basic characteristics of the patients are presented in Table 1.

Table 1 Basic characteristics of the patients

Immobilization and CT scanning

Patients were immobilized in the supine position with a vacuum bag (Medtec Medical, Inc., Buffalo Grove, IL) or a thermoplastic mask (Guangzhou Klarity Medical & Equipment Co., Ltd., Guangzhou, People’s Republic of China). All of the patients were simulated using a Brilliance Big Bore CT (Philips Brilliance CT Big Bore Oncology Configuration, Cleveland, OH, USA) under free breathing conditions. Ten-phase CT images were acquired at a 3-mm slice thickness during scanning using respiratory-correlated four-dimensional computed tomography (4DCT) via a Real-time Position Management System (Varian Medical System, Inc., Palo Alto, CA). Maximum intensity projection (MIP) and average intensity projection (AIP) images were reconstructed. The CT images, including the MIP and AIP images, were transferred to an Eclipse treatment planning system (Version 10.0, Varian Medical System, Inc., Palo Alto, CA) for target delineation, OAR contouring, treatment planning and plan evaluation.

Target defining and OAR contouring

The internal target volume (ITV) was delineated by incorporating the gross tumor volume (GTV) on ten phases of the 4DCT scans under the pulmonary windows. To account for the set-up uncertainties and potential baseline tumor shift, a uniform 5 mm planning target volume (PTV) was created expanding around the ITV. For OAR contouring, the whole lung was limited to the air-inflated lung parenchyma, and the GTV and trachea/ipsilateral bronchus were excluded according to the Radiation Therapy Oncology Group (RTOG) 0915 report [26]. The CW was segmented from the corrected lung edges with a 2 cm expansion in the lateral, anterior, and posterior directions, excluding the lung volume and mediastinal soft tissue [16, 17, 27]. If the 2 cm expansion extended outside the body, then the contour extended only to the external patient surface. To avoid cumbersome delineation of the entire CW, we defined it within a 3 cm limit in the head-to-feet direction from the PTV [27]. To evaluate the incidence of RIRF after SBRT treatment, the rib that was within or closest to the target was delineated under a window level of 750 and a window width of 1400.

Treatment planning

Dose regimes of 4 × 12 Gy were prescribed; 4 × 12 Gy represented 48 Gy in 4 fractions. The treatment was planned on the averaged 4DCT image using Eclipse treatment planning system (Version 10.0). All plans were designed on a TrueBeam LINAC with a 6 MV flattening filter free (FFF) photon beam and a maximum dose rate of 1400 MU/min. Treatment plans were created using dual partial arcs, preventing irradiation of the contralateral lung. The collimator angles for all plans were set to 30° and 330° to minimize the contribution of the tongue-and-groove effect to the dose. Optimization was performed using the progressive resolution optimizer (PRO_10028) algorithm. The objectives were adjusted to ensure a maximum dose higher than 120% of the prescribed dose center in the GTV. The dose was prescribed at 95% of the PTV covered by the prescription dose. Dose calculation was performed using the anisotropic analytical algorithm (AAA_10028) with a grid resolution of 1 mm while accounting for heterogeneity correction. All of the dose constraints and critical organ dose-volume limits should meet the criteria of the RTOG 0915 protocol and other publications [26, 28].

Radiobiological models

The 3-year TCP data was predicted using the Liu et al. model [29] with the isocenter dose as a predictor. The model basically considers the tumor regrowth locally after radiation therapy, and thus can be applied to predict the TCP value as a function of follow-up time. The calculating formula and key modelling parameters were obtained from the original publication, with two respective sets of radiobiological parameters to predict the TCP data for stage T1 and T2 tumors. We employed the Wennberg et al. model [30] to predict the probability of 2-year grade ≥ 2 RIP. The model was Lyman-Kutcher-Burman (LKB) based and showed a dose–response relationship between RIP and the equivalent uniform dose (EUD) of the lung. The Din et al. model [18], a Cox proportional hazard (CPH)-based model, was applied to predict 2-year grade ≥ 2 CWP. The model used a combination of a normalized total dose of 99 Gy (NTD99Gy) and body mass index (BMI) for prediction. The beta coefficients of the model were acquired by taking natural logarithms of the hazard ratio (HR) values shown in Table 2 of the reference. The baseline hazard h0(t) value was obtained from the nomogram by privately contacting the author. The h0(t) value is approximately equal to 0.05 after careful measurement. The risk of 3-year RIRF was predicted using the Stam et al. model [23]. The model was a traditional normal tissue complication probability (NTCP) model in which the time to toxicity was taken into account. To calculate the probability of rib fracture within 3 years, we multiplied the NTCP value by a latency distribution, f(τ), for the time to toxicity using the parameters described in Table 2 of the reference, for which a log-normal distribution is assumed. To estimate the probability of rib fracture within 3 years, the cumulative density function (CDF) was used to calculate the f(τ) value. The calculating process was performed using an in-house developed program on MATLAB 7.0 (MathWorks, USA).

Table 2 The IFR, TCP and UTCP values for the four patients with typical tumor characteristics

Fraction regime individualization

The method described is based on maximizing the P+ function, where P+ = TCP·(1-NTCP) [31, 32]. In this study, P+ = TCP·(1-NTCPlung)·(1-NTCPCW)·(1-NTCPrib). The method investigates the P+ values when the physical dose changes. The dose volume histogram (DVH) data of the GTV (DVHGTV), lung (DVHlung), CW (DVHCW) and rib (DVHrib) were extracted from the treatment planning system at a resolution of 2 cGy. The dose bins of DVHGTV, DVHlung, DVHCW and DVHrib were then multiplied by a certain percentage of the dose in that bin which ranged from 1 to 200% in steps of 1%. Accordingly, 200 groups of DVH data were obtained. The TCP values of GTV (TCPi) and the NTCP values of the lung (NTCPlung,i), CW (NTCPCW,i) and rib (NTCPrib,i) from each dose fraction of the DVH were calculated separately.

The P+i values for each of the 200 sets of DVH data were calculated. The maximum value for P+i correspond to the optimal therapeutic gain. A flow chart of the fraction regime individualization is presented in Fig. 1. As fraction regimes of SBRT are often less than five fractions and single fraction regimes were reported to be unfit by the linear quadratic (LQ) model beyond a fractional dose of up to 20 Gy [33], the DVHGTV, DVHlung, DVHCW and DVHrib data were also converted to 3 fractions regimes using the LQ model. In other words, two groups of P+i values, P+i in 3 fractions (P+i, 3f) and P+i in 4 fractions (P+i, 4f) were obtained.

Fig. 1
figure 1

Process of fraction regime individualization. The P+ values considering the TCP and NTCP values of the lung, CW and rib for each dose bin of the DVH were calculated. The physical dose corresponding to the maximum P+ value is referred to as IFR. Abbreviations: DVH = Dose volume histogram; TCP = Tumor control probability; NTCP=Normal tissue complication probability; P+ = UTCP=Uncomplicated tumor control probability; IFR = Individualized fraction regime

Results

Clinical parameters

The tumor diameter ranges from 1.4–5.0 cm. The GTV and PTV ranges from 1.6–70.6 cc and 13.5–128.9 cc, respectively. Detailed clinical parameters for the patients are presented in Table 1.

Individualized fraction regime

The IFR and corresponding TCP and UTCP values across the 33 patients are shown in Fig. 2. The figure shows that all the patients have patient-specific IFR, TCP and UTCP values. All the patients achieve a TCP value higher than 92.0%. The IFR ranges from 3 × 10.8 Gy to 3 × 12.5 Gy and from 4 × 9.2 Gy to 4 × 10.7 Gy for the 3 and 4 fraction regimes, respectively.

Fig. 2
figure 2

IFR, TCP and UTCP values across 33 patients. a IFR with 3 fractions; b IFR with 4 fractions; c TCP values with 3 fractions; d TCP values with 4 fractions; e UTCP values with 3 fractions; and f UTCP values with 3 fractions. Abbreviations: FD = Fractional dose of the IFR. TCP = Tumor control probability, UTCP = Uncomplicated tumor control probability

Four patients with typical tumor characteristics demonstrate that the method developed is personalized. Patient A has an off-CW small lesion (diameter of 1.9 cm and GTV of 3.4 cc) and patient B has off-CW large lesion (diameter of 3.7 cm and GTV of 26.7 cc). The tumor of patient C is adjacent to the CW and patient D has the largest target to normal lung volume ratio (PTV of 83.8 cc and the normal lung volume of only 1671.7 cc, predicted to be at high risk of RIP). A CT image of the four patients is shown in Fig. 3. Figure 4 shows the UTCP values as a function of physical dose for the four patients. The physical dose corresponding to the maximum UTCP value is referred to as the IFR, regardless of whether 3 or 4 fraction regimes are utilized. Table 2 shows that four patients with typical tumor characteristics possess IFR and individualized UTCP values while maintaining a TCP value higher than 93.0%. The UTCP values and fractional dose of the four patients in descending order are A > B > C > D. Patients with a large tumor has a lower TCP and UTCP and a smaller fractional dose than patients with a small tumor (patient B vs. patient A). Patients with tumor adjacent to the OAR (patient C) or at high risk of RIP (patient D) exhibits lower UTCP values and smaller fractional dose than patients with a tumor located distant from OAR (patient A and B).

Fig. 3
figure 3

CT images of the four patients with typical tumor characteristics. a Patient with an off-CW small tumor; b Patient with an off-CW large tumor; c Patient’s tumor was adjacent to the CW; and d Patient at a high risk of RIP

Fig. 4
figure 4

UTCP values as a function of physical dose for the four patients with typical tumor characteristics. The physical dose corresponding to the maximum UTCP value in the curve is referred to as the IFR. The IFR of 3 fractions and 4 fractions are displayed in red and blue, respectively. a Patient with an off-CW small tumor; b Patient with an off-CW large tumor; c Patient’s tumor was adjacent to the CW; and d Patient at a high risk of RIP. Abbreviations: UTCP = Uncomplicated tumor control probability

Discussion

In the current clinical practice of SBRT treatment for lung cancer, some “rigid” fraction regimes are commonly used. Clinical practice does not take into account tumor diversity, such as the size and proximity to normal tissues, which will lead to an underdose in some patients and to adverse side effects due to overdose in others. In this study, we first developed a method to calculate the IFR to avoid any underdose or overdose for NSCLC patients undergoing SBRT.

Efforts to search for the optimal dose of SBRT for stage I NSCLC are ongoing. Park et al. found that a biologically effective dose (BED) > 100 Gy was required to achieve a > 85% local control rate regardless of tumor size. The optimal dose for small tumors of < 3 cm appeared to be within a range below 150 Gy BED. The escalation of BED to high levels (> 150 Gy) might be required for patients with a tumor size larger than 3 cm [34]. Guckenberger et al. reported that doses of > 100 Gy BED to the CTV based on 4D dose calculation resulted in excellent local control rates [35]; Kestin et al. found a significant dose–response relationship for local control of NSCLC following image-guided SBRT with an optimal PTVmean BED10 > 125 Gy [36]; Lee et al. observed that tumors ≤2 cm had no local recurrence regardless of dose, whereas for tumors > 2 cm, an escalated BED of approximately 150 Gy10 provided significantly higher local tumor control [37]; Koshy et al. concluded that higher doses (> 150 Gy BED) were associated with a significant survival benefit in patients with T2 tumors [38]; and Zhang et al. reported that a medium or medium to high BED (range, 83.2–146 Gy) for SBRT may currently be more beneficial and reasonable in stage I NSCLC [39]. However, the studies above drew inconsistent conclusions for the following reasons. First, the studies were derived from samples at various tumor stages and a wide variety of fraction regimes. Second, the studies mainly concentrated on local control data, and toxicity to the normal tissues was less considered. In other words, we believe that the optimal dose of SBRT for lung cancer is patient specific when the factors of tumor stage and proximity to the OAR are all considered. The so-called“one-size-fits-all” fraction regime does not exist.

The absolute values of LC and NTCP were partially dependent on the radiobiological models used. In the study, the regrowth model was employed to predict the TCP value for the following reasons: (1) It is the one and only TCP-predicting model to separate between stage T1 and T2 tumors. It was also proven to offer a better fit to the clinical data compared with the universal survival curve (USC) and the modified linear quadratic and linear (mLQL) models. (2) It is the only model to estimate local control for a certain follow-up time for primary lung cancer patients. (3) The model is highly in accordance with other published data in which the isocenter dose (also denoted as the maximum dose) was correlated with tumor local control [40, 41]. We used the Wennberg et al. model to calculate the risk of RIP because it uses bilateral lung exclusive of the GTV as the definition of lung volume, which is generally recognized in lung SBRT [8, 13, 15]. The Din et al. model and the Stam et al. model are unique models that calculate the incidence of CWP and RIRF. Moreover, 7 of the 33 patients in the study had pulmonary metastases from different primaries, and we assume a similar dose-response relationship between the primary and secondary lung tumors according to the results of Guckenberger et al. [42].

The applicability of the LQ model for local control modelling has been widely reported. Guckenberger et al. suggested that traditional LQ formalism could model patients with stage I NSCLC undergoing SBRT more accurately than LQ-L formalism based on 395 patients from 13 German and Austrian centers [40]. Shuryak et al. also found that the LQ model provided a significantly better fit to local control data for NSCLC than any of the models requiring extra terms at a high dose range [43]. Santiago analyzed 1975 patients and demonstrated that the LQ model was a robust method for predicting 3-year local control data [41]. Unfortunately, there is limited evidence on the validity of the LQ model for RIP prediction. Only Scheenstra et al. reported that the alpha/beta ratio of 1.3 Gy using the traditional LQ model was applicable for RIP prediction [44]. Whether the USC model characterized by an additional dose modification beyond a certain transitional dose (dT) is more suitable for modelling the RIP prediction has not been well established [45]. The lack of strong evidence for the applicability of the traditional LQ model for RIP prediction prompted us to use the Wennberg model, which is characterized by a USC model for converting the equivalent dose into 2 Gy fractions (EQD2) when the fractional dose of the DVH is greater than 5.8 Gy [30].

Although our study demonstrated that the method based on maximizing the P+ value is able to calculate the IFR in SBRT for lung cancer, there are some limitations. (1) We have to perform the treatment planning and export the dose data before calculating the IFR for each patient, which requires several hours to complete the process of fraction regime individualization. However, we may chart the IFR using the following variables when enough patients are analyzed because tumor size [29, 38], tumor-CW distance [46] and the PTV to total lung volume ratio [47] were widely reported to be correlated with LC, CWP, RIRF and RIP. (2) As the time points of TCP and NTCP data are inconsistent (2 year for CWP and RIP prediction and 3 year for TCP and RIRF prediction, we can’t determine whether the IFR was calculated in 2-year or 3-year follow up. (3) The study used the radiobiological models to calculating the P+ values; however, the correctness of the parameters is a bit questionable. Therefore, clinical validation is required to confirm the results.

Conclusions

The proposed method based on maximizing the P+ value is able to predict the IFR in SBRT for lung cancer patients. However, clinical validation is required to confirm the results.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

4DCT:

Four-dimensional computed tomography

AAA:

Anisotropic analytical algorithm

AIP:

Average intensity projection

BED:

Biologically effective dose

BMI:

Body mass index

CDF:

Cumulative density function

CPH:

Cox proportional hazard

CT:

Computed tomography

CWP:

Chest wall pain

dT :

Transitional dose

DVH:

Dose volume histogram

EQD2 :

Equivalent dose in 2 Gy fractions

EUD:

Equivalent uniform dose

FFF:

Flattening filter free

GTV:

Gross tumor volume

HR:

Hazard ratio

IFR:

Individualized fraction regime

ITV:

Internal target volume

LKB:

Lyman-Kutcher-Burman

LQ:

Linear quadratic

MIP:

Maximum intensity projection

mLQL:

Modified linear quadratic and linear

NSCLC:

Non-small cell lung cancer

NTCP:

Normal tissue complication probability

OAR:

Organs at risk

PRO:

Progressive resolution optimizer

PTV:

Planning target volume

RIP:

Radiation-induced pneumonitis

RIRF:

Radiation-induced rib fractures

RTOG:

Radiation Therapy Oncology Group

SBRT:

Stereotactic body radiotherapy

TCP:

Tumor control probability

USC:

Universal survival curve

UTCP:

P+ uncomplicated tumor control probability

References

  1. Baumann P, Nyman J, Hoyer M, 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]. J Clin Oncol. 2009;27(20):3290–6.

    Article  Google Scholar 

  2. Bonfili P, Di Staso M, Gravina GL, et al. Hypofractionated radical radiotherapy in elderly patients with medically inoperable stage I-II non-small-cell lung cancer [J]. Lung Cancer. 2010;67(1):81–5.

    Article  Google Scholar 

  3. Timmerman R, Paulus R, Galvin J, et al. Stereotactic body radiation therapy for inoperable early stage lung cancer [J]. JAMA. 2010;303(11):1070–6.

    Article  CAS  Google Scholar 

  4. Xia T, Li H, Sun Q, et al. Promising clinical outcome of stereotactic body radiation therapy for patients with inoperable stage I/II non-small-cell lung cancer [J]. Int J Radiat Oncol Biol Phys. 2006;66(1):117–25.

    Article  Google Scholar 

  5. Solda F, Lodge M, Ashley S, et al. Stereotactic radiotherapy (SABR) for the treatment of primary non-small cell lung cancer; systematic review and comparison with a surgical cohort [J]. Radiother Oncol. 2013;109(1):1–7.

    Article  Google Scholar 

  6. Zhang B, Zhu F, Ma X, et al. Matched-pair comparisons of stereotactic body radiotherapy (SBRT) versus surgery for the treatment of early stage non-small cell lung cancer: a systematic review and meta-analysis [J]. Radiother Oncol. 2014;112(2):250–5.

    Article  Google Scholar 

  7. Zheng X, Schipper M, Kidwell K, et al. Survival outcome after stereotactic body radiation therapy and surgery for stage I non-small cell lung cancer: a meta-analysis [J]. Int J Radiat Oncol Biol Phys. 2014;90(3):603–11.

    Article  Google Scholar 

  8. Barriger RB, Forquer JA, Brabham JG, et al. A dose-volume analysis of radiation pneumonitis in non-small cell lung cancer patients treated with stereotactic body radiation therapy [J]. Int J Radiat Oncol Biol Phys. 2012;82(1):457–62.

    Article  Google Scholar 

  9. Baker R, Han G, Sarangkasiri S, et al. Clinical and dosimetric predictors of radiation pneumonitis in a large series of patients treated with stereotactic body radiation therapy to the lung [J]. Int J Radiat Oncol Biol Phys. 2013;85(1):190–5.

    Article  Google Scholar 

  10. Chang JY, Liu H, Balter P, et al. Clinical outcome and predictors of survival and pneumonitis after stereotactic ablative radiotherapy for stage I non-small cell lung cancer [J]. Radiat Oncol. 2012;7:152.

    Article  Google Scholar 

  11. Harder EM, Park HS, Chen ZJ, Decker RH. Pulmonary dose-volume predictors of radiation pneumonitis following stereotactic body radiation therapy [J]. Practical Radiat Oncol. 2016;6(6):e353–9.

    Article  Google Scholar 

  12. Matsuo Y, Shibuya K, Nakamura M, et al. Dose-volume metrics associated with radiation pneumonitis after stereotactic body radiation therapy for lung cancer [J]. Int J Radiat Oncol Biol Phys. 2012;83(4):e545–9.

    Article  Google Scholar 

  13. Nakamura M, Nishimura H, Nakayama M, et al. Dosimetric factors predicting radiation pneumonitis after cyberknife stereotactic body radiotherapy for peripheral lung cancer [J]. Br J Radiol. 2016;89(1068):20160560.

    Article  Google Scholar 

  14. Shi S, Zeng Z, Ye L, Huang Y, He J. Risk factors associated with symptomatic radiation pneumonitis after stereotactic body radiation therapy for stage I non-small cell lung cancer [J]. Technol Cancer Res Treat. 2017;16(3):316–20.

    Article  Google Scholar 

  15. Yamamoto T, Kadoya N, Sato Y, et al. Prognostic value of radiation pneumonitis after stereotactic body radiotherapy: effect of pulmonary emphysema quantitated using CT images [J]. Clin Lung Cancer. 2018;19(1):e85–90.

    Article  Google Scholar 

  16. Dunlap NE, Cai J, Biedermann GB, et al. Chest wall volume receiving >30 Gy predicts risk of severe pain and/or rib fracture after lung stereotactic body radiotherapy [J]. Int J Radiat Oncol Biol Phys. 2010;76(3):796–801.

    Article  Google Scholar 

  17. Mutter RW, Liu F, Abreu A, et al. Dose-volume parameters predict for the development of chest wall pain after stereotactic body radiation for lung cancer [J]. Int J Radiat Oncol Biol Phys. 2012;82(5):1783–90.

    Article  Google Scholar 

  18. Din SU, Williams EL, Jackson A, et al. Impact of fractionation and dose in a multivariate model for radiation-induced chest wall pain [J]. Int J Radiat Oncol Biol Phys. 2015;93(2):418–24.

    Article  Google Scholar 

  19. Murray L, Karakaya E, Hinsley S, et al. Lung stereotactic ablative radiotherapy (SABR): Dosimetric considerations for chest wall toxicity [J]. Br J Radiol. 2016;89(1058):20150628.

    Article  Google Scholar 

  20. Woody NM, Videtic GM, Stephans KL, et al. Predicting chest wall pain from lung stereotactic body radiotherapy for different fractionation schemes [J]. Int J Radiat Oncol Biol Phys. 2012;83(1):427–34.

    Article  Google Scholar 

  21. Thibault I, Chiang A, Erler D, et al. Predictors of chest wall toxicity after lung stereotactic ablative radiotherapy [J]. Clin Oncol (R Coll Radiol). 2016;28(1):28–35.

    Article  CAS  Google Scholar 

  22. Asai K, Shioyama Y, Nakamura K, et al. Radiation-induced rib fractures after hypofractionated stereotactic body radiation therapy: risk factors and dose-volume relationship [J]. Int J Radiat Oncol Biol Phys. 2012;84(3):768–73.

    Article  Google Scholar 

  23. Stam B, van der Bijl E, Peulen H, et al. Dose-effect analysis of radiation induced rib fractures after thoracic SBRT [J]. Radiother Oncol. 2017;123(2):176–81.

    Article  Google Scholar 

  24. Coroller TP, Mak RH, Lewis JH, et al. Low incidence of chest wall pain with a risk-adapted lung stereotactic body radiation therapy approach using three or five fractions based on chest wall dosimetry [J]. PLoS One. 2014;9(4):e94859.

    Article  Google Scholar 

  25. Bongers EM, Haasbeek CJ, Lagerwaard FJ, Slotman BJ, Senan S. Incidence and risk factors for chest wall toxicity after risk-adapted stereotactic radiotherapy for early-stage lung cancer [J]. J Thorac Oncol. 2011;6(12):2052–7.

    Article  Google Scholar 

  26. Videtic GM, Hu C, Singh AK, et al. A randomized phase 2 study comparing 2 stereotactic body radiation therapy schedules for medically inoperable patients with stage I peripheral non-small cell lung cancer: NRG oncology RTOG 0915 (NCCTG N0927) [J]. Int J Radiat Oncol Biol Phys. 2015;93(4):757–64.

    Article  Google Scholar 

  27. Kong FM, Ritter T, Quint DJ, et al. Consideration of dose limits for organs at risk of thoracic radiotherapy: atlas for lung, proximal bronchial tree, esophagus, spinal cord, ribs, and brachial plexus [J]. Int J Radiat Oncol Biol Phys. 2011;81(5):1442–57.

    Article  Google Scholar 

  28. Martin A, Gaya A. Stereotactic body radiotherapy: a review [J]. Clin Oncol (R Coll Radiol). 2010;22(3):157–72.

    Article  CAS  Google Scholar 

  29. Liu F, Tai A, Lee P, et al. Tumor control probability modelling for stereotactic body radiation therapy of early-stage lung cancer using multiple bio-physical models [J]. Radiother Oncol. 2017;122(2):286–94.

    Article  Google Scholar 

  30. Wennberg BM, Baumann P, Gagliardi G, et al. NTCP modelling of lung toxicity after sbrt comparing the universal survival curve and the linear quadratic model for fractionation correction [J]. Acta Oncol. 2011;50(4):518–27.

    Article  Google Scholar 

  31. Agren A, Brahme A, Turesson I. Optimization of uncomplicated control for head and neck tumors [J]. Int J Radiat Oncol Biol Phys. 1990;19(4):1077–85.

    Article  CAS  Google Scholar 

  32. Pizarro F, Hernandez A. Optimization of radiotherapy fractionation schedules based on radiobiological functions [J]. Br J Radiol. 2017;90(1079):20170400.

    Article  Google Scholar 

  33. Ruggieri R, Stavrev P, Naccarato S, et al. Optimal dose and fraction number in SBRT of lung tumours: a radiobiological analysis [J]. Phys Med. 2017;44:188–95.

    Article  Google Scholar 

  34. Park S, Urm S, Cho H. Analysis of biologically equivalent dose of stereotactic body radiotherapy for primary and metastatic lung tumors [J]. Cancer Res Treat. 2014;46(4):403–10.

    Article  Google Scholar 

  35. Guckenberger M, Wulf J, Mueller G, et al. Dose-response relationship for image-guided stereotactic body radiotherapy of pulmonary tumors: relevance of 4D dose calculation [J]. Int J Radiat Oncol Biol Phys. 2009;74(1):47–54.

    Article  Google Scholar 

  36. Kestin L, Grills I, Guckenberger M, et al. Dose-response relationship with clinical outcome for lung stereotactic body radiotherapy (SBRT) delivered via online image guidance [J]. Radiother Oncol. 2014;110(3):499–504.

    Article  Google Scholar 

  37. Lee S, Song SY, Kim SS, et al. Feasible optimization of stereotactic ablative radiotherapy dose by tumor size for stage I non-small-cell lung cancer [J]. Clin Lung Cancer. 2018;19(2):e253–61.

    Article  Google Scholar 

  38. Koshy M, Malik R, Weichselbaum RR, Sher DJ. Increasing radiation therapy dose is associated with improved survival in patients undergoing stereotactic body radiation therapy for stage I non-small-cell lung cancer [J]. Int J Radiat Oncol Biol Phys. 2015;91(2):344–50.

    Article  Google Scholar 

  39. Zhang J, Yang F, Li B, et al. Which is the optimal biologically effective dose of stereotactic body radiotherapy for stage I non-small-cell lung cancer? A meta-analysis [J]. Int J Radiat Oncol Biol Phys. 2011;81(4):e305–16.

    Article  Google Scholar 

  40. Guckenberger M, Klement RJ, Allgauer M, et al. Applicability of the linear-quadratic formalism for modelling local tumor control probability in high dose per fraction stereotactic body radiotherapy for early stage non-small cell lung cancer [J]. Radiother Oncol. 2013;109(1):13–20.

    Article  Google Scholar 

  41. Santiago A, Barczyk S, Jelen U, Engenhart-Cabillic R, Wittig A. Challenges in radiobiological modelling: can we decide between LQ and LQ-L models based on reviewed clinical nsclc treatment outcome data? [J]. Radiat Oncol. 2016;11:67.

    Article  Google Scholar 

  42. Guckenberger M, Klement RJ, Allgauer M, et al. Local tumor control probability modelling of primary and secondary lung tumors in stereotactic body radiotherapy [J]. Radiother Oncol. 2016;118(3):485–91.

    Article  Google Scholar 

  43. Shuryak I, Carlson DJ, Brown JM, Brenner DJ. High-dose and fractionation effects in stereotactic radiation therapy: analysis of tumor control data from 2965 patients [J]. Radiother Oncol. 2015;115(3):327–34.

    Article  Google Scholar 

  44. Scheenstra AE, Rossi MM, Belderbos JS, et al. Alpha/beta ratio for normal lung tissue as estimated from lung cancer patients treated with stereotactic body and conventionally fractionated radiation therapy [J]. Int J Radiat Oncol Biol Phys. 2014;88(1):224–8.

    Article  Google Scholar 

  45. 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 [J]. Int J Radiat Oncol Biol Phys. 2008;70(3):847–52.

    Article  Google Scholar 

  46. Nambu A, Onishi H, Aoki S, et al. Rib fracture after stereotactic radiotherapy for primary lung cancer: prevalence, degree of clinical symptoms, and risk factors [J]. BMC Cancer. 2013;13:68.

    Article  Google Scholar 

  47. Ueyama T, Arimura T, Takumi K, et al. Risk factors for radiation pneumonitis after stereotactic radiation therapy for lung tumours: clinical usefulness of the planning target volume to total lung volume ratio [J]. Br J Radiol. 2018;91(1086):20170453.

    Article  Google Scholar 

Download references

Funding

This work was sponsored by National Natural Science Foundation of China (81602667), Medical Scientific Research Foundation of Guangdong Province (A2015534, B2016048) and the Creative and Facilitating Program of Shantou University.

Author information

Authors and Affiliations

Authors

Contributions

Each author had participated sufficiently in the work to take public responsibility for appropriate portions of the content. JYL wrote the manuscript draft. JYL performed the treatment planning and collected the data. PXL analyzed the data. BTH conceived, designed the experiments and revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Bao-Tian Huang.

Ethics declarations

Ethics approval and consent to participate

The protocol was approved by the Ethics Committees of the Cancer Hospital of Shantou University Medical College. Because this is not a treatment-based study, our institutional review board waived the need for written informed consent from the participants. However, the patient information was kept anonymous to protect their confidentiality.

Consent for publication

Not applicable.

Competing interests

All authors certify that his manuscript has not been published in whole or in part nor is it being considered for publication elsewhere. The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lu, JY., Lin, PX. & Huang, BT. Calculating the individualized fraction regime in stereotactic body radiotherapy for non-small cell lung cancer based on uncomplicated tumor control probability function. Radiat Oncol 14, 111 (2019). https://doi.org/10.1186/s13014-019-1318-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13014-019-1318-9

Keywords