1. Introduction
Non-small cell lung cancer (NSCLC) continues to be the leading cause of cancer-related mortality [1] . Dose escalation has been proposed to improve local control [2] -[5] and capitalize on the therapeutic ratio of thoracic radiotherapy, which is mainly limited by the potential toxicity to organs at risk (OARs). In brief, NSCLC is the most commonly diagnosed type of lung cancer that accounts for 85% of all cases characterized by slow growth and spread compared to other small cell lung cancer [6] . NSCLC encompasses a number of different types of lung cancer, which are categorized as “squamous” or “non-squamous”. Smoking, ionizing radiation, occupational risks, environment, indoor air pollution and diseases as risk factors for lung cancer are the ultimate causes and risk factors. NSCLC is associated with persistent cough, shortness of breath, coughing up phlegm (sputum) with signs of blood, aches or pains when breathing or coughing, loss of appetite, fatigue, and loss of weight. NSCLC is detected using several options that include cytology, bronchoscopy, and needle biopsy. The Prognosis based on three year survival rates for lung cancer patient who received stereotactic body radiation therapy reveals a survival rate of 55.8% [7] .
Stereotactic body radiation therapy (SBRT) is a technique that can attain a high cumulative biologically effective dose (BED) inside the tumor with a steep gradient outside the tumor, preventing dose to OARs. Achieving a high BED has been shown to improve patient survival and local control rates for many tumors [8] . However, SBRT confronts several issues including patient positioning, breathing patterns, and high dose escalation. All these issues can be detrimental to the treatment outcome.
Advancements in imaging and radiation treatment delivery systems have provided radiation oncologists with better target definition. For this reason, SBRT has recently emerged as a favored treatment modality for early-stage NSCLC. Modern SBRT has yielded durable local control results, comparable to surgery [9] , and excellent overall survival outcomes even in inoperable patients.
Several authors have carried out studies on lung SBRT, examining the impacts of beam arrangement, image guidance, radiation dose constraints, dose prescription, and treatment planning algorithm on patient outcomes. For example, Papiez and Timmerman [10] discussed dose shaping and dose delivery techniques using standard beam arrangements with coplanar and noncoplanar beams. They established that SBRT treatment can be a new standard therapy for inoperable lung cancer. Chang and Timmerman [11] used SBRT along with recent advancements in diagnostic and real-time imaging to make precise measurements of tumor motion during treatment, and showed that the beams can be adjusted in real time to maintain the dose on the tumor. Through clinical trials of lung, liver, spine, pancreas, kidney, and prostate tumors, they proved that the effective doses delivered by SBRT are a comparable and noninvasive alternative to surgery.
However, much work remains to determine the best treatment parameters. Ding et al. [12] published a rationale with the Radiation Therapy Oncology Group (RTOG 0236) demonstrating the need for tissue heterogeneity corrections in the lung dose calculation. According to them, heterogeneity corrections have a major impact on the SBRT technique for lung tumors and normal lung tissues. This warning has been echoed by Herman et al. [13] . Another group directed by Timmerman et al. [14] uses tumor motion tracking methods such as breath-hold to achieve the lowest possible radiation dose to OARs. This group has suggested adopting a dose-response relationship for tumor control using SBRT, with increasing doses. Similarly, Xiao et al. [15] conducted a retrospective study to determine dose constraints and prescriptions for SBRT lung protocols and establish a heterogeneity correction for dose calculations. Likewise, Takeda et al. [16] performed a retrospective study of clinical outcomes in SBRT lung patients with Stages 1A and 1B NSCLC, and concluded that delivering a total of 50 Gy in 5 fractions at the periphery of the PTV is a viable treatment option. Finally, Takahashi et al. [17] compared the dose distributions calculated by three algorithms and revealed the importance of choosing an adequate algorithm for SBRT.
This research contributes to the literature by determining useful radiobiological parameters of SBRT treatments, such as the NTCP of organs at risk and the TCP of the tumor, in a sample of 39 patients. The results are calculated using BioSuite software, and will be compared to published data as quality assurance [18] - [20] .
2. Methods and Materials
Ethical statement: Hospital Sisters Health System institutional review board approval was obtained for the treatment.
2.1. Patients
A cohort of thirty-nine patients with stage IA/IB NSCLC underwent definitive thoracic radiation as a single- mode treatment. All of them were treated by SBRT, but they were prescribed a variety of treatment plans. The demographics and stages of the patients are listed in Table 1. Volumetric modulated arc therapy (VMAT) and intensity modulated radiation therapy (IMRT) are the methods of choice for dose delivery; Figure 1 illustrates the VMAT method applied to one patient.
SBRT provided an adequate solution for challenging central and advanced lung tumors, because of the close proximity of organs at risk. VMAT dispenses the treatment more rapidly than IMRT, making the procedure
(a)(b)(c)
Figure 1. A 72-year-old woman with lung cancer was treated with SBRT for a total dose of 50 Gy in 5 fractions. (a) (top) VMAT with 3 arcs; (b) (middle) isodose lines; (c) (bottom) dose color “wash” normalized to 50 Gy.
more bearable for the patient. The specific advantages of its shorter duration are reduced intrafraction motion by the patient and less loss of biological effect.
2.2. Dosimetry
The patients were placed supine with their arms above their heads, while the tumor motion was assessed with 4DCT. The modalities of this procedure are described elsewhere [21] . Due to the treatment length for SBRT, close attention was paid to comfort as well as the reproducibility and stability of their position. Briefly, the computed tomography (CT) scans of all cases were acquired with 512 × 512 pixels and 0.25 cm slice spacing on the flat tabletop of a GE LightSpeed CT Scanner (Milwaukee, USA). The planning target volume (PTV) was created using a 5 mm wide isotropic expansion of the clinical target volume (CTV). The OARs, such as lungs, heart, chest, skin, esophagus and spinal cord, were delineated using axial CT images. All cases were planned using the VMAT (RapidArc) technique or IMRT technique, in 6X mode at 600 MU/min dose rate. The treatment plans were inversely optimized using the progressive resolution optimizer (PRO) in Eclipse TPS version 11 (Varian Medical Systems, Palo Alto, CA), such that at least 95% of the PTV received the prescription dose of 50 Gy in 4 - 5 fractions. The optimized plans were calculated using an analytical anisotropic algorithm (AAA) with tissue heterogeneity correction. Each calculated plan was then normalized such that 100% of the prescribed dose covered 95% of the PTV. The calculation grid was set to 2.5 mm for all cases. The DVHs (Figure 2) from the approved plans of the thirty-nine patients, for both tumors and OARs, were imported into BioSuite [22] to calculate important radiobiological parameters.
2.3. Radiobiology Parameters
Several radiobiological parameters are used to describe causes and effects when tissue is subjected to radiation ionization. Among these are the equivalent uniform dose (EUD), normal tissue complication probability (NTCP), and tumor control probability (TCP). All three can be analyzed with variation of model analysis. These parameters and their modes of calculation have been described elsewhere [23] [24] . Briefly, Niemierko characterized EUD as the theoretical dose which would cause the same biological effect as the actual inhomogeneous dose, if the former were distributed uniformly throughout the volume of the structure of interest [25] [26] . According to Niemierko, EUD is given by
(1)
In this equation, vi is unitless and denotes the ith partial volume receiving dose Di in Gy. The value a is a unitless model parameter specific to the OAR or tumor of interest. This formula has seen extensive use because of its simplicity. NTCP and TCP can be calculated using dose volume histograms (DVHs), combined
Figure 2. Dose volume histogram of lung SBRT treated at 50 Gy in 5 fractions.
with additional parameters portraying the reactions of various tissues to ionizing radiation, including the α/β ratio, a parameter n that controls dose volume effects, TD50 or dose sensitivity, a parameter m that controls the slope of dose response curve, and the reference volume v [27] - [30] . This study uses Biosuite software to model NTCP and TCP, based in part on inputs from the Lyman-Kutcher-Burman (LKB) model parameters. NTCP was calculated for OARs and linear quadratic (LQ) Poisson “Marsden” model [31] , to estimate TCP for lung tumors.
3. Results
Data from thirty-nine patients who underwent SBRT lung treatment were used for this study (Figure 3(a), Figure 3(b)).
Their prescribed doses were either 5 daily fractions of 10 Gy or 4 daily fractions of 12.5 Gy, for a total dose of 50 Gy using SBRT. The patient characteristics are summarized in Table 1.
The primary organs involved were the lung, heart, esophagus, spinal cords, chest wall, and skin. The median tumor size was 3.7 cm (range 1.6 to 6.0 cm). Eleven patients (28.20%) had a tumor size of <3 cm, and 28 patients (71.80%) had a tumor size of ≥3 cm. The tumors were located in different areas of the lung: upper (13), middle (17), or lower (9). Among the 39 patients, four received high biologically effective doses of radiation (BED > 151.2 Gy), 11 received medium-high doses (BED, 112.5 Gy), five received medium-low doses (105.6 Gy), and 19 received low doses (BED = 100 Gy). Table 2 shows the characteristics of the four BED groups.
According to the LQ model, the median BED is 100 and the range is 100 - 151.2 Gy. The mean PTV volume is 31.60 ± 25.55 cc. The mean EUDs in Gy to the different organs at risk are 5.19 ± 2.84 (lungs), 5.66 ± 4.95 (spinal cord), 61.45 ± 29.18 (chest wall), 3.31 ± 5.92 (heart), 6.45 ± 5.18 (esophagus), and 12.22 ± 5.94 (skin). On average, the heart received the lowest EUD, while the chest and ribs received the highest (61.45 ± 29.18 Gy). Figure 4 summarizes all observed EUDs for the six OARs and thirty-nine cases.
(a)(b)
Figure 3. A 70-year-old woman with lung cancer was treated with SBRT for a total dose of 50 Gy given in 5 fractions.
Figure 4. EUDs for organs at risk in lung SBRT. (The y-axis is in cGy and the x-axis is the patient number).
The mean NTCPs were 3.75% ± 2.61%, 36.25% ± 36.42%, and 0.59% ± 1.48% for the lungs, chest, and esophagus, respectively. The NTCP values for cords, heart, and skin were 0.00%. All NTCP values are shown in Figure 5.
The mean TCP was 99.72% ± 0.44%. Using the PTV encompassing dose or maximum (isocentric) dose, our data indicate that the mean isocenter dose is 53.08 ± 2.41 Gy. Figure 6 displays the BED value at the isocenter, the isocenter dose, and the TCP for all patients.
4. Discussion
This research employs radiobiological modeling to calculate the biological effects of radiation therapy based on the linear-quadratic (LQ) model. In this study, we analyzed radiobiological parameters for thirty-nine patients who underwent lung SBRT. Several predictive factors such as TCP, BED, EUD, and NTCP were characterized using the LKB model. All input parameters of this study are derived from the reports of Burman et al. [32] and Luxton et al. [33] .
For patients who undergo lung SBRT, the most obvious toxicities are pneumonitis, esophagitis, and chest- wall pain; pulmonary and cardiovascular complications are also very common. Accordingly, within our sample the NTCP values for the cords, heart, and skin are 0.00%. The highest observed NTCP value is for the chest, while smaller but still significant NTCP values are observed for the total lung, treated lung, and esophagus. The large value of NTCP for the chest can be attributed to several factors, including tumor location and some issues that are patient-dependent. For example, it is believed that female gender and a small tumor-chest wall separation are both risk factors for rib fracture [34] .
Uematsu et al. [35] described a 5-year experience with 50 patients who underwent SBRT, using 50 - 60 Gy in 5 - 10 fractions with a median follow-up of 36 months. In their sample only one patient (2%) had a rib fracture detected by follow-up CT, and this one did not require medical treatment. Uematsu’s study is supported by the work of Nambu et al. [36] , who states that the frequency and degree of chest wall pain is greatly dependent on the total dose or fraction. They suggested that a BED10 ranging from 96 to 119 Gy at the isocenter is safe enough to prevent any chest wall adverse reaction. However, different frequencies are expected at the horizons as we follow-up our patients.
We used the LQM model with α/β = 10 Gy for our BED calculations and compared the results to published reports with a variety of SBRT schemes. In addition, we examine the dose-response relationship for local tumor
Figure 5. NTCPs for organs at risk in lung SBRT (The y-axis is a % pro- bability, and the x-axis is the patient number).
Figure 6. TCP, BED at the isocenter, and isocenter dose for lung SBRT The y-axis is in Gy and the x-axis represents the patient number).
control. To achieve better local tumor control, the consensus in the literature is that the BED should be over 100 Gy. The mean BED value for our study, estimated at the isocenter, was 119.57 Gy. This is consistent with reviews by Onishi et al. [37] and Baumann et al. [38] which state that higher BED values result in better outcomes. They showed, considering a variety of treatment schemes, that BED ≥ 100 Gy results in a recurrence rate of only 8.4%, while patients receiving BED < 100 Gy have a recurrence rate of 26.4%.
According to Zhang et al. [39] , BED > 146 Gy resulted in a local control rate of 100%, two years after treatment. The BEDs in our study range from 100 - 151.2 Gy. It is believed that in order for a lower BED to maximize tumor control, it is paramount to avoid complications and effectively treat tumors in central locations.
In this study, the SBRT treatments were performed with either IMRT or VMAT. Both techniques are susceptible to interplay between the multileaf collimator (MLC) and tumor respiratory motion. Li et al. [40] evaluated the interplay effect for intra-fraction cumulative dose and concluded that less than 1% discrepancy in the PTV and ITV minimum doses was observed using an energy mapping algorithm. However, in a SBRT scenario, where 3 - 5 dose fractions are delivered, the effect of this interplay on dose distributions is not well defined.
Examining the TCP and NTCP results, dose calculation with an analytical anisotropic algorithm (AAA) may become more and more clinically significant for lung cancer. It is paramount to use the most accurate dose calculation algorithm possible, in order to avoid underestimation of the required dose and increasing the probability of local disease recurrence. On the other hand, overestimation of the dose will increase toxicity to normal tissues. For example, Latifi et al. [41] show a higher recurrence rate for patients whose treatments were planned with Pencil-Beam (PB) compared to collapsed cone convolution (CCC) algorithms, even though the prescribed nominal dose and constraints were identical.
This study has several limitations, including the lack of patient follow-up in order to establish a survival curve. Only randomized assessments among patients undergoing the SBRT regimen will avoid the disparities and selection biases present in these comparisons.
Another limitation is that LQM was used to calculate radiobiological parameters. Timmerman et al. [42] [43] showed that for SBRT with a dose greater than 50 Gy, cell survival curves seem linear rather than downward bending, as LQM would suggest. For this reason, we believe that LQM may overestimate the cytotoxic effects of SBRT. A new method called the universal survival curve (USC) model has been introduced which associates the LQM for the low dose region with a multitarget model asymptote for high dose region [44] . In addition, inconsistencies on endpoints due to variation in breathing patterns of the lung, making more difficult to rely solely on DVH-based parameters. EUD, NTCP, and TCP can be calculated directly from the DVH. These parameters are therefore dependent on of the treatment planning system (TPS) employed, which includes interpolation, binning, and normalization methods. The quantitative evaluation of EUD, NTCP, and TCP reveals a strong dependence on the actual clinical treatment situation. Finally, SBRT radiation dose and fractionation depends on the oncologist’s experience and also varies due to a lack of consensus in the literature regarding recommended BED. In addition, tailored treatment planning has to account for organ motion and real-time tomographic acquisition.
5. Conclusion
NTCP, EUD, and TCP may serve as a useful tool for prospective and retrospective analysis of radiotherapy treatment plans in lung SBRT.