Abstract
Background/Aim: This study aimed to evaluate whether the predictive performance of tumor control probability (TCP) by carbon ion radiotherapy (CIRT) could be improved by incorporating cancer type-specific carbon ion sensitivity data obtained from non-small cell lung cancer (NSCLC).
Materials and Methods: A TCP model based on the linear-quadratic (LQ) formalism and fractionated irradiation was employed, with σ defined as an index of inter-tumor heterogeneity on radiosensitivity. The LQ model parameters α and β for NSCLC were obtained from 13 cell lines subjected to carbon ion irradiation under clinically relevant conditions, followed by clonogenic assays. The human salivary gland (HSG) cell line was used as a control. The α and β values were incorporated into the TCP model, and agreement with the clinical data (i.e., the local control rates for 48 NSCLCs treated with CIRT, as previously reported) was evaluated using the coefficient of determination (R2) derived from least-squares fitting.
Results: When σ was set to 0.20, the NSCLC-derived TCP curve showed better agreement with the clinical data than the HSG cell-derived TCP curve (R2: 0.72 vs. 0.53, respectively). Similarly, when σ was set to 0.15, the NSCLC-derived curve showed better agreement with the clinical data (R2: 0.56 vs. 0.38, respectively). Furthermore, when weighting factors of 1, 3, and 10 were applied to the two data points with TCP=1.0, the NSCLC-derived TCP curve showed consistently better agreement with the clinical data than the HSG cell-derived TCP curve (R2: 0.74 vs. 0.55, 0.78 vs. 0.63, and 0.79 vs. 0.65, respectively).
Conclusion: Incorporation of cancer type-specific carbon ion sensitivity data can improve the predictive performance of TCP modeling compared with the conventional HSG cell-based approach.
- Carbon ion radiotherapy
- tumor control probability
- relative biological effectiveness
- linear-quadratic model
- non-small cell lung cancer
Introduction
Carbon ion radiotherapy (CIRT) is a promising radiotherapy modality that offers superior dose conformality and a greater cell-killing effect than photons (1). The cell-killing effect of carbon ions at an equivalent physical dose is influenced by variations in linear energy transfer (LET) along the beam depth (2). In Japan, the effect of LET differences on the cell-killing effect of carbon ions is addressed during dose prescription by using the clinical dose, which is defined as the physical dose weighted by the relative biological effectiveness (RBE) (3). The RBE values are derived from measurement of clonogenic survival of human salivary gland tumor (HSG) cells irradiated with carbon ions (3). In practice, the clonogenic survival data are fitted to the linear-quadratic (LQ) model, and the resulting α and β parameters, which are determined for each beam energy, are incorporated into the treatment planning systems. In that sense, clinical dose prescriptions for CIRT are performed using the HSG-derived α and β values, regardless of cancer type.
Results from basic research have demonstrated that cancer cells exhibit a wide range of sensitivities to carbon ions (4-6). This finding suggests that the use of cancer type-specific carbon ion sensitivity data may allow for more accurate prediction of antitumor effects; however, this remains to be elucidated. In fact, previous studies have evaluated the tumor control probability (TCP) of CIRT in non-small cell lung cancer (NSCLC) (7) and prostate cancer (8) using the α and β values derived from HSG cells, or even by treating these values as free parameters for fitting. To address this gap, we aimed to investigate whether the predictive performance of TCP in CIRT could be improved through incorporation of cancer type-specific carbon ion sensitivity data.
Materials and Methods
Study overview. A TCP model based on the concepts of the LQ model and fractionated irradiation was used. The LQ model parameters α and β were obtained from cell experiments and incorporated into the TCP model. Then, the agreement between the TCP model and the clinical data was evaluated. NSCLC was chosen as a model for analysis. The carbon ion sensitivity data for 13 NSCLC cell lines were obtained, and data from HSG cells were used as a reference. For the clinical data, the local control rates of 48 NSCLCs treated with CIRT, as reported by Miyamoto et al., were utilized (9).
The TCP model. It is assumed that the relative frequency of the number of surviving clonogenic cells in a given tumor after treatment follows a Poisson distribution (10). TCP is the probability that no surviving clonogenic cells remain in a tumor, which corresponds to the zero-order term of the Poisson distribution, expressed as follows: (10, 11):
(1)
where Ns denotes the average number of surviving clonogenic cells per tumor after treatment.
For fractionated radiotherapy, in which the LQ model and a time factor accounting for cancer cell proliferation are taken into consideration, Ns is expressed as follows (12):
(2)
where N0 is the number of clonogenic cells in the tumor prior to treatment; n is the number of fractions; d is the fractionated dose (Gy); Td is the average doubling time of the tumor cells (days); and T is the overall treatment time (days).
Individual tumors within a given cohort should exhibit different radiosensitivities. To account for this, and to achieve a better fit with clinically observed data, the distribution of α among individual tumors within the cohort was incorporated (13). When gi is the probability density of tumors having α=αi, the overall mean TCP for the cohort is expressed as follows (11):
(3)
Assuming that αi follows a Gaussian distribution, gi is expressed as follows:
(4)
where αm and σ are the mean and standard deviation of α, respectively. Thus, the TCP for the cohort is obtained as follows:
(5)
Using Equation 5, the TCP was calculated by performing numerical integration of αi over the range −5σ to 5σ, with a step size of 0.01. The following values were used: N0=109, Td=7, (7), and T=42 (9). Fitting the TCP model to the clinical data was performed using least-squares fitting (7). For several analyses, fitting was performed by assigning a weighting factor to the clinical data points, with a local control rate of 1.0, as reported by Kanai et al. (7). The extent of agreement between the TCP model and the clinical data was evaluated using the coefficient of determination (R2), derived from least-squares fitting. These calculations were performed using Python (version 3.7.3, Python Software Foundation, Beaverton, OR, USA).
Cell experiments. The 13 human NSCLC lines were cultured in RPMI-1640 (Sigma-Aldrich, St. Louis, MO, USA) supplemented with 10% fetal bovine serum (Life Technologies, Carlsbad, CA, USA; Table I). Clonogenic assays were performed as previously described (14). Briefly, cells were seeded in 6-well plates (n=4) and cultured at 37°C until they attached to the plates (6-12 h, depending on cell line). The cells were exposed to carbon ions, at 1, 2, 3, or 4 Gy, at the Gunma University Heavy Ion Medical Center. The beam specification was as follows: 290 MeV/nucleon; and an average LET at the center of a 6 cm spread-out Bragg peak (SOBP) of approximately 50 keV/μm (4). After incubation for an additional 10-14 days, the cells were fixed with methanol and stained with crystal violet. Colonies comprising at least 50 cells were counted under an inverted microscope. The surviving fractions were calculated after normalization to the corresponding un-irradiated controls. The survival data were fitted to the LQ model, and the α and β values were obtained. In cases where the β values were extremely small positive, zero, or negative, the α/β ratio becomes difficult to interpret (i.e., extremely large, undefined, or negative, respectively). As no consensus exists regarding the appropriate handling of such situations, these cases were treated as missing values in the present study.
Linear quadratic (LQ) model parameters for the cell lines used in this study.
In CIRT in Japan, the HSG cell line is used for standardization of clinical doses (3); however, this cell line has been reported to be cross-contaminated with HeLa cells in multiple laboratories. In addition, authentication of uncontaminated HSG cells is theoretically impossible due to the absence of genotypic data for the original line at the time of their establishment (15). From this standpoint, the present study adopted the α and β values for HSG cells (α=0.76 and β=0.076) obtained at the National Institutes for Quantum Science and Technology (formerly the National Institute of Radiological Sciences, Chiba, Japan), where the standardization of clinical carbon ion doses was first established, and under irradiation settings consistent with those used in the aforementioned experiments with NSCLC cells (7).
Clinical data. As a reference for the CIRT outcomes of NSCLC, clinical data reported by Miyamoto et al. were employed (9). Their study presents the results of a dose-escalation trial for stage I NSCLC, in which the prescribed dose ranged from 59.4 Gy (RBE) to 95.4 Gy (RBE), delivered in 18 fractions (Table II). Among the 48 tumors analyzed, 26 were adenocarcinomas and 22 were squamous cell carcinomas. In the present study, the tumor control was defined according to the local control rates reported by Miyamoto et al. The local control rate from the clinical results was plotted as a function of the physical dose (Gy) instead of the clinical dose; this was done by applying an RBE value of 2.38 at the center of a 60-mm SOBP, as previously described (7).
Clinical data for carbon ion radiotherapy of non-small cell lung cancer used for tumor control probability modeling (9).
Results
The α and β values of the 13 NSCLC cell lines were obtained from clonogenic assays following carbon ion irradiation (Table I). TCP curves for each individual NSCLC cell line were generated based on Equation 1 using the cell line-specific α and β values, since the distribution of α (considered in Equation 3-5) can be neglected in this context (Figure 1). As expected, the TCP curves showed substantial shifts along the dose axis, reflecting the diversity of radiosensitivity among cell lines. Nevertheless, no histology-specific trends in the curve shape were identified. Based on these findings, in the subsequent analyses the mean α and β values across all examined NSCLC cell lines were used as the representative α and β values for NSCLC cells.
Tumor control probability (TCP) curves generated using the α and β values derived from each of the 13 non-small cell lung cancer cell lines. The lower right-hand panel compares the TCP curves by histological subtype. Adeno: Adenocarcinoma. Cell lines with the β values <0.001 are not shown.
Next, we compared the predictive performance of the TCP model between the setting using the α and β values derived from HSG cells with that using the values derived from NSCLC cells on the clinical outcome of CIRT for 48 stage I NSCLCs (9). In Equation 5, which describes the TCP model, σ represents the extent of inter-tumor heterogeneity in radiosensitivity (7). Given that the true σ value for the population is unknown, we employed 0.20, the σ value estimated from least-squares fitting of the TCP model to clinical outcomes and in which α, β, and σ were treated as free parameters. The α and β values for HSG cells, or those for NSCLC cells, were substituted for αm and β, respectively, in Equation 5. As a result, the R2 value for the NSCLC cell-derived TCP curve outperformed that for the HSG cell-derived TCP curve (R2: 0.72 vs. 0.53, respectively; Figure 2).
Tumor control probability (TCP) curves generated using the α and β values for human salivary gland (HSG) cells or for non-small cell lung cancer (NSCLC) cells, with σ fixed at 0.20. The R2 values obtained from least-squares fitting of the TCP model to the clinical data are shown.
The study by Kanai et al., which analyzed the fitting of the HSG cell-derived TCP curve to CIRT outcomes using the same clinical dataset used in the present study, concluded that the σ value that best described the clinical outcomes was 0.15±0.05 (7). Therefore, we next compared the predictive performance of the TCP model using the cell type-specific α and β values but employing a σ value of 0.15. As a result, the R2 value for the NSCLC cell-derived TCP curve outperformed that for the HSG cell-derived TCP curve (R2: 0.56 vs. 0.38, respectively; Figure 3).
Tumor control probability (TCP) curves generated using the α and β values for human salivary gland (HSG) cells or non-small cell lung cancer cells, with σ fixed at 0.15. The R2 values obtained from least-squares fitting of the TCP model to the clinical data are shown.
The study by Kanai et al. also performed least-squares fitting of the TCP model to clinical outcomes using the HSG cell-derived α and β values, treating σ as a free parameter (7). Therefore, we next compared the predictive performance of the TCP model using cell type-specific α and β values. In that analysis, Kanai et al. performed curve fitting by assigning “large weighting factors” to the two data points showing TCP=1.0, based on the assumption that these points were reliable due to saturation of the local control rate at 100%; however, the actual values for the weighting factors used were not provided in the manuscript. Therefore, we performed curve fitting by employing three different weighting factors: 1, 3, and 10. As a result, the R2 value for the NSCLC cell-derived TCP curve outperformed that for the HSG cell-derived TCP curve for all weighting factors analyzed (R2: 0.74 vs. 0.55, 0.78 vs. 0.63, and 0.79 vs. 0.65, for the weighting factors 1, 3, and 10, respectively; Figure 4). Taken together, these data indicate that the predictive performance of the TCP during CIRT could be improved by incorporating cancer type-specific carbon ion sensitivity data.
Tumor control probability (TCP) curves generated using the α and β values for human salivary gland (HSG) cells or non-small cell lung cancer (NSCLC) cells by applying a weighting factor of 1, 3, or 10 to data points showing TCP=1.0. The σ value was treated as a free parameter for least-squares fitting. The R2 values obtained from least-squares fitting of the TCP model to the clinical data are shown.
Discussion
We found that incorporating in vitro carbon ion sensitivity data derived from 13 NSCLC cell lines into the TCP model improved concordance with the clinical outcomes of CIRT for NSCLC. As a reference control, we used the study by Kanai et al. (7), which examined the fit of TCP, obtained using the sensitivity data derived from HSG cells, to the same clinical dataset. In that study, the authors demonstrated that the TCP based on HSG cell data was generally consistent with clinical outcomes, and so they concluded that the parameters that best accounted for the results were as follows: α=0.75±0.5, σ=0.15±0.05 (assuming β=0.076 and Td=7). In the present study, we evaluated the fit of the TCP model to clinical outcomes under multiple assumptions regarding the parameters, including those concluded by Kanai et al. Across all settings, the TCP based on NSCLC cell line-derived data consistently demonstrated a superior fit over the TCP based on HSG cell line-derived data. These findings indicate that the predictive performance of the TCP for CIRT can be improved by incorporating cancer type-specific carbon ion sensitivity data. The results suggest that the use of cell line-derived, cancer type-specific α and β values may facilitate prediction of optimal clinical dose-fractionation schemes for cancer types with limited clinical data. Moreover, cancer type-specific α and β values could also be implemented to LET-modulated CIRT that has been attracting interest in recent years (16, 17). To further validate these findings, additional studies involving other cancer types are warranted.
For decades, attempts have been made to predict the clinical outcomes of photon therapy using a TCP model based on the LQ parameters derived from cultured cells. Yaes et al. investigated TCP by incorporating α and β values from cell lines derived from six cancer types (18). In that study, the authors demonstrated that, for head and neck squamous cell carcinoma, the TCP curve based on cell line data was considerably steeper than that derived from clinical data. This observation is consistent with our findings, as shown in Figure 1. As a potential explanation for this discrepancy, Peters et al. suggested inter-tumor heterogeneity of radiosensitivity among clinical tumors (19). In the present study, we also assumed inter-tumor heterogeneity of radiosensitivity and so introduced σ as the parameter representing this variability. Since the true value of σ in the population is unknown, we performed multiple validations, each with σ set to a different value. In all cases, incorporation of σ brought the steepness of the cell line-derived TCP curve closer to that of the clinical data-derived TCP. In this regard, it can be said that the notion proposed by Peters et al. was confirmed by the present study. In addition, when σ was treated as a free parameter when fitting the TCP model to the clinical data, the resulting σ value was comparable with the median absolute deviation of the α values obtained from the NSCLC cell lines (i.e., approximately 0.25). Although further validation is warranted, this finding indicates the possibility that σ for a given cancer type cohort may be estimable from in vitro cell experiments.
As another potential explanation for the discrepancy between cell line- and clinical data-derived TCP, Weichselbaum et al. suggested intra-tumor heterogeneity of radiosensitivity within a clinical tumor, in which the presence of a small subset of radioresistant clones is assumed (20). It has been considered that this discrepancy is resolved by assuming that, within a tumor consisting of 109 cancer cells (defined as N0 in the present study), approximately 200 cells are intrinsically radioresistant. Weichselbaum et al. further speculated that the cause of this intrinsic cancer cell radioresistance is an enhanced DNA damage repair capacity. If this is the case, incorporating variations of β into TCP modeling may improve the predictive performance of clinical outcomes, since repair capacity is considered to be represented by β in the clonogenic survival curves approximated by the LQ model. In the present study, variation in the β term was not taken into account in the TCP model, which used the settings employed in previous studies (7, 13). Supporting this approach, Kanai et al. reported that variation in the β value of HSG cells with respect to varying carbon ion LET was modest (7); however, other studies have shown that cancer cells exhibit a wide range of β values in response to carbon ions under clinically relevant experimental settings (6, 21), suggesting the need for considering incorporating the β variation into the TCP model in the future. Nevertheless, doing so would render TCP modeling more complex, warranting further investigation.
When attempting to predict clinical TCP using cancer cell-derived radiosensitivity data, attention should be paid to the potential discrepancy in radiosensitivity between cultured cell lines and clinical tumors within a given cancer type. For example, it is well established that the α/β value for prostate cancer in the clinical setting is as low as 1-2 (22), whereas that for prostate cancer cell lines is typically much higher (23). This discrepancy may be attributable to the fact that prostate cancer cell lines commonly used in radiation research are derived from metastatic lesions and lack androgen sensitivity in most cases; thus, they exhibit biological characteristics distinct from those of primary prostate cancers, which are the targets of definitive radiotherapy in clinical practice. For this reason, direct application of cell line-derived radiosensitivity data to clinical TCP prediction is unlikely to be feasible in the context of this cancer type. Indeed, Kang et al. performed TCP modeling for CIRT in prostate cancer using equations based on the LQ model; however, their analysis treated α and β as free parameters (8). From this standpoint, clinical TCP prediction based on cell line-derived data requires careful evaluation of the biological properties of the cell lines employed.
The present study has the following limitations. First, the clinical dataset did not include physical dose data; therefore, direct comparison in units of Gy were not feasible. Second, the TCP model incorporated multiple parameters for which accurate estimation of true values is inherently difficult (i.e., σ, N0, and Td). These parameters constitute sources of uncertainty regarding the predictive performance of the TCP model. Third, as described in the Methods section, negative β values were treated as missing values due to the lack of consensus regarding their interpretation, which may have influenced the TCP estimation results. Lastly, in this study, differences in the weighting of TCP across clinical data plots, arising from variations in the number of cases per plot, were not taken into account. This decision was made because such differences have not been addressed in previous studies and incorporating them would substantially increase the complexity of the model.
In summary, the present study demonstrates that incorporation of cancer type-specific carbon ion sensitivity data can improve the predictive performance of TCP modeling compared with the conventional HSG cell-based approach. These findings suggest that cell type-specific biological parameters should be considered carefully when estimating clinical TCP in CIRT. Further validation with larger datasets and across additional cancer types is warranted to establish the generalizability of this approach, and to support its translation into clinical treatment planning.
Acknowledgements
This work was supported by Takeda Science Foundation, Shimadzu Science Foundation, Nishikawa Medical Foundation, the Japanese Society for Radiation Oncology, the Headquarters of Gunma University, and GHMC.
Footnotes
Authors’ Contributions
E.Y. performed experiments, analyzed data, and wrote the manuscript. T. Oike conceptualized the project, performed experiments, and finalized the manuscript. M.S. and M.T. supervised TCP analysis. T. Ohno supervised entire project and acquired funding. All Authors read and agreed with the final version.
Conflicts of Interest
The Authors declare no conflicts of interest in relation to this study.
Artificial Intelligence (AI) Disclosure
No artificial intelligence (AI) tools, including large language models or machine learning software, were used in the preparation, analysis, or presentation of this manuscript.
- Received September 26, 2025.
- Revision received October 7, 2025.
- Accepted October 20, 2025.
- Copyright © 2026 The Author(s). Published by the International Institute of Anticancer Research.
This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.










