Variability of α/β ratios for prostate cancer with the fractionation schedule: caution against using the linear-quadratic model for hypofractionated radiotherapy

Background Prostate cancer (PCa) is known to be suitable for hypofractionated radiotherapy due to the very low α/β ratio (about 1.5–3 Gy). However, several randomized controlled trials have not shown the superiority of hypofractionated radiotherapy over conventionally fractionated radiotherapy. Besides, in vivo and in vitro experimental results show that the linear-quadratic (LQ) model may not be appropriate for hypofractionated radiotherapy, and we guess it may be due to the influence of fractionation schedules on the α/β ratio. Therefore, this study attempted to estimate the α/β ratio in different fractionation schedules and evaluate the applicability of the LQ model in hypofractionated radiotherapy. Methods The maximum likelihood principle in mathematical statistics was used to fit the parameters: α and β values in the tumor control probability (TCP) formula derived from the LQ model. In addition, the fitting results were substituted into the original TCP formula to calculate 5-year biochemical relapse-free survival for further verification. Results Information necessary for fitting could be extracted from a total of 23,281 PCa patients. A total of 16,442 PCa patients were grouped according to fractionation schedules. We found that, for patients who received conventionally fractionated radiotherapy, moderately hypofractionated radiotherapy, and stereotactic body radiotherapy, the average α/β ratios were 1.78 Gy (95% CI 1.59–1.98), 3.46 Gy (95% CI 3.27–3.65), and 4.24 Gy (95% CI 4.10–4.39), respectively. Hence, the calculated α/β ratios for PCa tended to become higher when the dose per fraction increased. Among all PCa patients, 14,641 could be grouped according to the risks of PCa in patients receiving radiotherapy with different fractionation schedules. The results showed that as the risk increased, the k (natural logarithm of an effective target cell number) and α values decreased, indicating that the number of effective target cells decreased and the radioresistance increased. Conclusions The LQ model appeared to be inappropriate for high doses per fraction owing to α/β ratios tending to become higher when the dose per fraction increased. Therefore, to convert the conventionally fractionated radiation doses to equivalent high doses per fraction using the standard LQ model, a higher α/β ratio should be used for calculation. Supplementary Information The online version contains supplementary material available at 10.1186/s13014-022-02010-9.

Background With the development of high-precision radiotherapy, fractionation schedules to treat various tumors are changing [1]. Hypofractionated radiotherapy is being increasingly employed in clinics as stereotactic body radiotherapy (SBRT) and moderately hypofractionated intensity-modulated radiotherapy (IMRT), which have become valuable therapeutic approaches for a variety of tumors owing to the improved dose distribution. In addition, for prostate cancer (PCa), hypofractionated IMRT and SBRT seem to have radiobiological advantages based on the linear-quadratic (LQ) model estimation.
Since definitive hypofractionated radiotherapy is a relatively novel treatment, optimal dose fractionation schedules often need to be inferred from mathematical calculation, and an LQ model-based formula is frequently used to convert the conventionally fractionated radiation doses to high doses per fraction by clinicians due to its convenience and simplicity [2]. Recently, however, several investigators demonstrated that the standard LQ model may not be applicable to hypofractionated radiotherapy especially in SBRT [2][3][4][5][6], although other researchers insist that the LQ model can be used to estimate the antitumor effects of hypofractionated radiotherapy [7,8].
The α/β ratio is a key factor in the LQ model [2]. Basically, the α/β ratio of a tumor is obtained from an in vitro dose-survival curve of tumor cells [2,6], but this method cannot be applied to human tumors in patients. The α/β ratio can also be obtained from in vivo tumor or normal tissue responses to different fractionation schedules, and following this in vivo method, a mathematical method was elaborated to estimate the α/β ratio from clinical data employing various fractionation schedules [9,10]. Using the method, the α/β ratios for various tumors have been reported, and PCa was found to have a low α/β ratio [9,10], which was lower than the α/β ratio for normal tissue late reactions. Accordingly, moderately hypofractionated IMRT and SBRT are being increasingly used in the treatment of PCa. However, since the reliability of the LQ model in SBRT was questioned in recent studies [4][5][6], it may be necessary to re-evaluate the validity of converting conventionally fractionated doses to hypofractionated doses with the LQ model. Previous studies only tried to demonstrate that PCa has a low α/β ratio [9][10][11][12], and variability of the α/β ratios with the dose per fraction has not been investigated. Therefore, we carried out an analysis using an established mathematical calculation method to estimate the variation in the α/β ratio of PCa according to the daily fractional dose and to verify the applicability of the LQ model in hypofractionated radiotherapy.

Clinical data collection
We searched for relevant articles in PubMed with key words of "radiotherapy" or "radiation therapy" and "prostate cancer" or "prostatic carcinoma". The inclusion criteria were as follows: (1) patients with PCa undergoing conventionally fractionated radiotherapy, moderately hypofractionated radiotherapy, or SBRT, with or without androgen deprivation therapy (ADT) and (2) 5-year biochemical relapse-free survival (5y-bRFS), number of patients, total dose, and fraction number or dose per fraction available from the articles. Articles that lacked the necessary fitting data or that used other fractionation schedules, such as hyperfractionated radiotherapy, were excluded. Conventionally fractionated radiotherapy was defined as that using 1.8-2.1 Gy per fraction. Moderately hypofractionated radiotherapy was defined as that using 2.19-3.5 Gy in our study. SBRT was defined as that using 6.5 Gy per fraction or greater according to National Comprehensive Cancer Network (NCCN) guidelines.
Five-year bRFS according to the ASTRO or Phoenix definition was evaluated. The ASTRO definition of biochemical relapse is three consecutive rises in prostate-specific antigen (PSA) from the nadir [13]. The Phoenix definition of biochemical relapse is a rise of PSA over 2 ng/mL from the nadir [14]. Risk stratification of PCa was mostly made according to the NCCN guidelines risk group classification and in part D' Amico's classification.

Estimation of the α/β ratios
Statistical analyses were carried out exactly following the method of Miralbell and coworkers [15,16]. Briefly, standard LQ models for tumor control at 5 years of the form: or were fitted to the obtained data. In the formula, P is interpreted as the tumor control probability (5y-bRFS); doses to equivalent high doses per fraction using the standard LQ model, a higher α/β ratio should be used for calculation. Keywords: LQ model, α/β ratio, Hypofractionated radiotherapy, High dose per fraction, Prostate cancer D is the total dose; N is the number of fractions during the whole radiotherapy; d is the dose per fraction; and k represents the natural logarithm of an effective target cell number. α represents unrepairable lethal damage caused by a one-track action and β represents repairable sublethal damage caused by a two-track action in the DNA damage repair kinetics [17]. So, the α/β ratio can be considered as the balance between the two forms of damage.
In order to reduce the errors of fitting parameters, we performed jackknife method. We removed one group every time from the 22 groups of data (taking conventionally fractionated radiotherapy group as an example); we could get 22 groups of sample data, which is C21 22in mathematics. Using these 22 groups of data, the average and standard error of the sample can be obtained, and then the population mean and the confidence interval of fitting parameters can be obtained. Moreover, we go a step further and calculate the Bias-Corrected and accelerated(BCa) intervals according to the Jung et al. 's methods [19]. For the parameters of the LQ model, significant figures were rounded to the 2nd decimal place.
In order to verify the accuracy of the results, we substituted the fitting results: k, α and α/β values and known parameters, i.e., total dose (D), single dose (d) and total fraction number (N), into the original TCP formula (1) or (2) and got a calculated P (5y-bRFS). Using goodness of fit test by chi-square test, the formula was , where T i means theoretical frequency (original P) and O i means observed frequency (calculated P). Then, we put calculated X 2 into the critical value table and obtained corresponding P value. This method was checked whether there was statistical difference between the calculated P and the P from the original study. The statistical software SPSS 22.0 was used.

Results
There were 45 articles (23,281 PCa patients) incorporated in this study published during 2003 to February 2021. Detailed information on the 45 articles is shown in Table 1. Among them, 5 articles reported on more than 1000 patients [20][21][22][23][24], and 3 reported on 500-1000 patients [25][26][27]. Among all 45 studies, a total of 38 articles (16,442 PCa patients) could be grouped according to fractionation schedules of radiotherapy and complete information could be extracted based on the inclusion and exclusion criteria; 15 were enrolled in the conventionally fractionated radiotherapy group, 24 were enrolled in the moderately hypofractionated radiotherapy group, and 8 were enrolled in the SBRT group. Nine articles were duplicated since the studies investigated both conventional fractionation and moderate hypofractionation. Detailed data from each article are shown in Tables 2, 3 and 4 according to the three different regimens of radiotherapy. To explore the relationship between α/β ratios and risks of PCa, we also divided each group into three subgroups by the risks in patients receiving radiotherapy with different fractionation schedules. Of all the 45 studies, 21 studies (14,641 PCa patients) could be grouped and the characteristics are shown in Additional file 1: Tables S1 to S3.
The preliminary results of verification of fitting results were shown in Table 7. The X 2 were all < 1 in all three risk groups and the P values were all > 0.995 that meant there was no statistical difference between the calculated TCP and the TCP from the original study. In other words, our fitting was accurate.
In summary, for PCa patients receiving conventionally fractionated radiotherapy, moderately hypofractionated radiotherapy, and SBRT, the mean α/β ratios were 1.78, 3.46, and 4.24 Gy, respectively. Meanwhile, as the risks of PCa increased, the k and α values decreased.

Discussion
The α/β ratio proposed in the early 1970's derives from the LQ models [28,29]. Factors that can influence α and/ or β independently increase or decrease the α/β ratio. The major influencing factors are internal factors from cells themselves and external factors from physical or chemical effects [17,30]. The internal factors include cell cycle regulation, cell repopulation, and DNA damage repair after irradiation. The external physical factors include temperature (hyperthermia), oxygenation (hypoxia), characteristics of radioactive rays-like linear energy transfer, and the dose rate. The external chemical factors are some anticancer drugs such as cisplatin, EGFR inhibitors, and PARP1 inhibitors. Thus, there are multiple factors that affect the α/β ratio and modify the radiosensitivity of tumors. Our study showed that the α/β ratio tended to become higher when the dose per fraction increased. The α/β ratios may increase also dynamically during treatment, from approximately 4 Gy for 'short' fractionation schedules to about 1.5 Gy for long schedules, which probably reflects the process of accelerated repopulation in normal acute skin reactions [31,32]. For late-responding tissues and slow-growing tumors like PCa, however, there may be no repopulation during radiotherapy [31], and the α/β ratio increase is not due to tumor cells repopulation. Also the time factors should not be considered in late-responding tissues [33,34]. Thus, we did not take the time factor into consideration when converting doses using the standard LQ model.
Recent randomized trials demonstrated that hypofractionated radiotherapy was not superior to conventional radiotherapy in PCa. In the Radiation Therapy Oncology Group (RTOG) 0415 [35], Hypofractionated Irradiation for Prostate Cancer trial (HYPRO) [25], and the Fox Chase trial (ClinicalTrials.gov identifier: NCT00062309) [36], biological effective doses (BEDs) in hypofractionated vs. conventionally fractionated radiotherapy groups were calculated as 186.7 vs. 162.4 Gy, 211.0 vs 182.0 Gy, and 196.6 vs 177.3 Gy, respectively, using an α/β ratio of 1.5 Gy. All BEDs in the hypofractionated groups were 19-29 Gy higher than BEDs of the conventionally fractionated groups. Nevertheless, the higher BEDs did not lead to satisfactory improvements in the outcome. This may be attributable to the inaccurate conversion of radiation doses using the LQ model. In our study, the α/β ratio tended to become higher when the dose per fraction increased. When the doses in the three trials were    [3,37,38]. Thus, several models were proposed modifying the standard LQ model to reasonably convert conventionally fractionated doses to equivalent single or hypofractionated doses. The lethalpotentially-lethal (LPL) model considered DNA lesion repair and could explain very effectively the shoulder on survival curves [39]; The modified LQ (MLQ) model made a better fit to the iso-effect data than the LQ model in a single high dose [40]. The "universal survival curve"(USC) model proposed by Park et al. combined two classical radiobiological models: the multitarget model and the standard LQ model that provide superior approximation of survival curves in the high-dose range. [41], and generalized LQ (gLQ) model encompasses the full dose range of possible dose delivery patterns and special radiotherapy schemes. [3]. Characteristics of these models have already been described [5]. Wang et al. [3] demonstrated that the problems in the LQ model derived from the amount of sublethal damage were reduced owing to conversion to lethal damage at a single high dose; if sublethal damage is converted to lethal damage, then the α/β ratio is elevated with a single high dose according to the definition of the α/β ratio. Other studies also revealed that cell death at high doses exceeded the probability of intracellular cell repair, and higher α/β ratios were shown with a linear survival curve [39,42]. An in vivo study involving a murine tumor model demonstrated that an equivalent single high dose converted from fractionated radiotherapy was lower than the actual dose. However, when a higher α/β ratio was used, the discrepancy became smaller [38]. Our data agreed with their results. At different fractional doses, the α/β ratio tended to be higher when the dose per fraction increased (1.78 Gy for conventional fractionation, 3.46 Gy for moderate hypofractionation, and 4.24 Gy for SBRT). Especially in the SBRT groups, the high α/β ratio was marked.
We also found that α and k values decreased with risk elevation in the conventional fractionation and moderate hypofractionation groups. These results were similar to those in the previous study [15]. A decrease in the α values with escalation of the risk group can be attributed to higher radio-resistance of tumor cells in higher risk patients. k represents the natural logarithm of an effective target cell number, and a decrease in k values means that the effective target cell number is reduced with escalation of the risk group.
A recent publication was based on the dose distribution delivered to patients and provided another method for fitting parameter [43]. This method used the linearquadratic Poisson TCP model with dose distribution like GTV (prostate gland), mpMRI-GTV, D 50 and dose volume histograms (DVH) to estimate α/β ratio and obtain solution space, initial parameter values, and optimal solution by optimizer. Our research approach was different from this method and used maximum likelihood principle in mathematical statistics to fit the α/β ratio according to the Miralbell model with tumor control probability (5y-bRFS), total dose, and number of fractions or dose per fraction; this method was much easier than the methods mentioned above.
The α/β ratio in the low-risk patients of the SBRT group was in the negative range. A study using external beam radiation therapy alone also had negative α/β  ratios [44]. Repeated measures of PSA at 6 institutions were analyzed and data from 3 institutions including RTOG showed negative α/β ratios. In the Peter Mac-Callum Cancer Center, the α/β ratio was − 2.05 Gy (95% CI −∞ − +∞ ). Another study found that the α/β ratio of arteriovenous malformation obliteration after radiosurgery was markedly negative (α/β = − 49.3 ± 5.3) [45]. However, neither study explained why the α/β ratio was negative. These results as well as ours suggest a limitation of this calculation method in that it could possibly yield unrealistic α/β ratios, especially when the patient number is small. Although, we imposed non-negativity restrictions on the α/β ratio of the low-risk group (− 10.7 Gy to 0.032 Gy), it is merely a mathematical statistics method and may affect real result. There are several limitations in our study. Since we divided the whole group into three subgroups according to the fractionation schedule, the dose ranges per fraction were relatively narrow in each fractionation group. This may increase the variability of the estimated α/β ratios, but we tried to solve this problem by including as many patients as possible. Subtle variations in patient evaluation including the definition of PSA failure and treatment including the dose prescription method among respective studies would also contribute to variability in the estimated α/β ratios; this problem is common to all studies of this kind, and is considered to be ameliorated by including a large number of patients. An analysis of over 14,000 patients showed that the derived α/β ratios were not different between studies using the ASTRO definition and those using the Phoenix definition [46]. The patient number in our study was larger than in any other studies investigating the α/β ratio for PCa. Also, the influence of ADT was not considered, as was the case with other previous studies, since a previous study indicated the minimal influence of ADT [15]. Another limitation is that our analysis was based on the prescription dose and the treatment outcome of the incorporated studies. Thus, 3D dose description or the DVHs were not incorporated into the analysis, because, for fitting parameters, we only needed prescription dose and treatment outcome according to TCP formula. Furthermore, the incorporated studies all used prescription dose which was widely used in clinical practice rather than dose distribution.
In conclusion, our study using mathematical statistics with 5y-bRFS data in PCa patients demonstrated that the α/β ratio was dependent on the fractionation schedule. In SBRT, the estimated α/β ratio was > 4 Gy. Therefore, to convert conventionally fractionated radiation doses to an equivalent single high dose, it may be necessary to use either a modified formula or a higher α/β ratio with the standard LQ model.