Skip to main content

ORIGINAL RESEARCH article

Front. Artif. Intell., 05 April 2024
Sec. Medicine and Public Health

Artificial neural network-assisted prediction of radiobiological indices in head and neck cancer

Saad Bin Saeed AhmedSaad Bin Saeed Ahmed1Shahzaib NaeemShahzaib Naeem2Agha Muhammad Hammad KhanAgha Muhammad Hammad Khan3Bilal Mazhar QureshiBilal Mazhar Qureshi4Amjad HussainAmjad Hussain5Bulent AydoganBulent Aydogan6Wazir Muhammad
Wazir Muhammad1*
  • 1Department of Physics, Florida Atlantic University, Boca Raton, FL, United States
  • 2Gamma Knife Radiosurgery Center, Dow University of Health Sciences, Karachi, Pakistan
  • 3Radiation Oncology, McGill University, Montreal, QC, Canada
  • 4Radiation Oncology, Aga Khan University Hospital, Karachi, Pakistan
  • 5Tawam Hospital, Al Ain, United Arab Emirates
  • 6Radiation and Cellular Oncology, University of Chicago, Chicago, IL, United States

Background and purpose: We proposed an artificial neural network model to predict radiobiological parameters for the head and neck squamous cell carcinoma patients treated with radiation therapy. The model uses the tumor specification, demographics, and radiation dose distribution to predict the tumor control probability and the normal tissue complications probability. These indices are crucial for the assessment and clinical management of cancer patients during treatment planning.

Methods: Two publicly available datasets of 31 and 215 head and neck squamous cell carcinoma patients treated with conformal radiation therapy were selected. The demographics, tumor specifications, and radiation therapy treatment parameters were extracted from the datasets used as inputs for the training of perceptron. Radiobiological indices are calculated by open-source software using dosevolume histograms from radiation therapy treatment plans. Those indices were used as output in the training of a single-layer neural network. The distribution of data used for training, validation, and testing purposes was 70, 15, and 15%, respectively.

Results: The best performance of the neural network was noted at epoch number 32 with the mean squared error of 0.0465. The accuracy of the prediction of radiobiological indices by the artificial neural network in training, validation, and test phases were determined to be 0.89, 0.87, and 0.82, respectively. We also found that the percentage volume of parotid inside the planning target volume is the significant parameter for the prediction of normal tissue complications probability.

Conclusion: We believe that the model has significant potential to predict radiobiological indices and help clinicians in treatment plan evaluation and treatment management of head and neck squamous cell carcinoma patients.

1 Introduction

Head and neck squamous cell carcinoma (HNSCC) is becoming a significant global health concern, with an estimated 30% increase in incidence rates by 2030 as per Global Cancer Observatory (GLOBOCAN) data, anticipating it to become a major global burden in time (Johnson et al., 2020). Radiation therapy (RT) is one of the modalities used for the management of HNSCC cases, often combined with surgery or chemotherapy making it an important modality in the treatment of HNSCC (Kawashita et al., 2020). The success of RT treatment relies on various factors, including the tumor’s stage, grade, treatment time, nutrition status, and treatment volume, as well as considerations such as tumor hypoxia (Huang and O’Sullivan, 2017). Additionally, the role of biological and radiological markers is evolving with clinical prognostic factors that play a crucial role in influencing treatment success (Konings et al., 2020; Peng et al., 2021).

In oncology, treatment outcomes revolve around two key factors: survival and toxicity. Survival assesses the effectiveness of interventions, measuring overall and progression-free durations (Abola et al., 2014). Simultaneously, managing and mitigating treatment-related toxicity is vital for the patient’s well-being and quality of life (Hunter et al., 2013).

The process of radiotherapy involves several crucial steps, with the planning stage serving as a pivotal checkpoint. During this phase, meticulous consideration is given to optimizing the delivery of radiation, ensuring maximum efficacy in targeting the tumor while minimizing potential harm to surrounding healthy tissues (Gardner et al., 2019). This planning step plays a foundational role in shaping the overall success and precision of radiotherapy treatment. The basic principle of planning is optimizing dose delivery to maximize the dose to the planning target volume (PTV) while minimizing harm to nearby tissues and organs at risk (OAR). Insufficient tumor coverage in radiotherapy can compromise treatment effectiveness, while overly reducing side effects may result in inadequate tumor treatment. Striking the right balance between optimal tumor control and minimizing side effects is crucial in the planning and implementation of radiotherapy. This delicate equilibrium is a key consideration for clinicians to ensure successful outcomes. This emphasizes the need for a nuanced approach to balancing the benefits and potential risks associated with radiotherapy (Pereira et al., 2014).

Like other treatment modalities, RT does come with known specific side effects (Dilalla et al., 2020). Despite leaps of advancements in radiation therapy and its standards, managing OAR toxicity remains a big frontier. Accurate risk estimation is crucial. For this purpose, integrating radiobiological indices like Tumor Control Probability (TCP) and Normal Tissue Complication Probability (NTCP) into treatment planning can quantify the response to irradiation, facilitating informed decision-making. The preference for one over the other creates a seesaw effect, leading to a trade-off where the loss of one aspect results in the gain of another (Yorke, 2023). Although precautionary measures are taken during both the planning and delivery stages, it is challenging to completely avoid toxicity to OAR. Dose-volume histograms (DVHs) conventionally describe 3D-dose distribution but have limitations. Integrating radiobiological indices with DVH, such as TCP and NTCP, provides a more comprehensive evaluation of treatment plans. Presently, the calculation methodologies for Tumor Control Probability (TCP) and Normal Tissue Complication Probability (NTCP) in radiation therapy heavily depend on radiobiological models. Specifically, the poison-based model estimates TCP for a given radiation dose, and the Lyman-Kutcher-Burman (LKB) model utilizes Dose-Volume Histogram (DVH) data to assess the likelihood of normal tissue complications (Warkentin et al., 2004). However, these models exhibit limitations as they oversimplify the intricate biological processes involved, overlooking factors such as tumor heterogeneity, organ- or tissue-specific mechanisms, and patient characteristics. To address these shortcomings, ongoing research direction is toward enhancing these models by integrating uncertainties, thereby promoting more nuanced clinical judgment in the context of treatment planning (Bufacchi et al., 2020).

In recent years, advancements in prediction outcomes have been notable, marked by the emergence of data-driven models (Isaksson et al., 2020). These models integrate dose-volume metrics, advanced bioinformatics tools, and disease- or patient-based prognostic factors. Several publications underscore the efficacy of these prediction models, demonstrating medium to good accuracy in forecasting clinical outcomes for head-and-neck cancer (Parmar et al., 2015; Giraud et al., 2019; Mahmood et al., 2021). These outcomes also encompass HNSCC high-grade RT toxicity including xerostomia, sticky saliva, radiation-induced trismus, and sensorineural hearing loss (Thor et al., 2021; Araújo et al., 2023). Concurrently, ongoing research aims to refine traditional radiobiological models by incorporating uncertainties for improved clinical judgment in treatment planning, addressing the complexity of biological processes and patient-specific factors (Gronberg et al., 2023).

The emergence of the machine learning in healthcare for making informed decisions, leading to effective patient management. The artificial neural network (ANN) showed significant progress in areas of outcome prediction, risk quantification personalized treatments, and diagnosis accuracy for both chronic and infectious diseases (Muhammad et al., 2019; Shafiq et al., 2022; Ataei et al., 2023). This also motivated us to utilize ANN in the current study.

This study aims to develop ANN for the prediction of radiobiological indices in HNSCC. The primary objective is to empower treatment planners with patient-specific information on TCP and NTCP, facilitating a more informed approach to prioritizing PTV and OAR during the treatment planning and evaluation stages. The proposed methodology involves utilizing an artificial neural network, a type of machine learning model, trained on a carefully curated dataset encompassing dose-volume metrics, patient-specific factors, and outcomes. ANN will be designed to predict radiobiological indices, particularly TCP and NTCP. Rigorous validation of the dataset will ensure the robustness of the model. The development of ANN for predicting radiobiological indices holds the potential to revolutionize HNSCC treatment planning. By providing treatment planners with patient-specific information, this approach aims to enhance the overall quality, efficacy, and outcome of radiation therapy, ultimately contributing to improved patient care and management.

2 Materials and methods

2.1 Data sets

Set I: The computed tomography (CT) images of 215 HNSCC patients were acquired from the publicly available repository “The Cancer Imaging Archive” (TCIA) (Clark et al., 2013). Patients were treated with curative intent RT from 2003 to 2013 at the University of Texas MD Anderson Cancer Center (Grossberg et al., 2018). The dataset comprises de-identified DICOM images, RT structures, planning parameters, DVH, and dose deposition data. The collection also consists of patient demographics, tumor grade, and stage, risk factors, recurrence, and survival information.

Set II: The dataset consists of CT scans of 31 HNSCC patients. Three CT scans were acquired for each patient throughout the treatment course. These include treatment planning, mid-treatment, or interim and post-treatment (i.e., follow-up) scans. Patients were treated with volumetric arc therapy (VMAT) at the Department of Radiation Oncology, the University of Miami Miller School of Medicine from 2011 to 2017 (Bejarano et al., 2019). Along with DICOM images, patient demographics, tumor histology, outcomes, RT structures, DVH, and planning parameters were available for public download at TCIA.

2.2 Methodology

The process begins with the acquisition of imaging datasets from TCIA, followed by the upload to the treatment planning system for verification and pre-processing. Both parotid OARs and PTV were selected for DVH data extraction from the treatment planning system and then processed toward a spreadsheet for the calculation of radiobiological indices. Patient particulars along with radiobiological indices were used for the training of the neural network (Figure 1).

Figure 1
www.frontiersin.org

Figure 1. The methodology and basic components utilized for the development of ANN.

2.3 Data acquisition

TCIA is the repository for clinical and imaging data. It hosts various data collections of cancer, and non-cancer patients, and phantom studies. We selected two different sets of HNSCC patients which have RT structures, RT dose, and DVH information along with imaging and patient demographics. An open-source software “NBIA Data Retriever V8.0” was used to download both sets containing DICOM files and spreadsheets of clinical information for all the patients.

2.4 Pre-processing

In Eclipse™ (treatment planning system) V 15.1, the anonymized DICOM files of both sets were uploaded with new IDs, i.e., HN001, and HN002. The CT images of the planning scan, RT structures, and RT dose files were imported to generate a dummy plan for each patient. These plans were evaluated for dose distribution in PTV and OARs. Verification and constancy of downloaded patient data were also ensured with these plans. Patients excluded from further processing were shown in a flow chart (Figure 2). Derived OAR structures were also contoured to extract dose parameters. These parameters along with others were used as input parameters for ANN. The last step in the pre-processing stage is the export of DVH of PTV, OAR, and derived structures in tabular format for the calculation of indices.

Figure 2
www.frontiersin.org

Figure 2. The inclusion/exclusion flow chart of patient datasets.

2.5 Radiobiological indices

RADBIOMOD® (V. beta 0.3) is a spreadsheet-based software tool for the calculation of radiobiological indices and evaluation of RT treatment plans (Chang et al., 2016). It offers multiple radiobiological models for the calculation of TCP and NTCP. We used the Gay & Niemierko (Gay and Niemierko, 2007) model for the calculation of both TCP and NTCP of PTV and parotids, respectively. Gay & Niemierko’s model utilizes the equivalent uniform dose (EUD) mechanistic formulation based on the linear–quadratic cell survival model (Webb and Nahum, 1993). To supplement ANN, we also used the Poisson model (Jackson et al., 1995) and (Zaider and Minerbo, 2000) for the calculation of TCP and the Lyman–Kutcher–Burman model (Lyman, 1985) for NTCPs. The calculation begins with the import of Tabular DVH from each dummy plan in Eclipse to the RADBIOMOD. Dose fractionation, total dose, treatment duration, and α/β ratios, i.e., The proportional contribution of single and double-strand breaks to cell death, were also used as input for the calculation of TCP and NTCP. Then computed radiobiological indices for each patient were routed toward training of an artificial neural network.

2.6 Artificial neural network architecture

A feed-forward neural network with back-propagation training was designed using a machine learning toolbox in Matlab® V.2022a. Under the umbrella of Neural Network Fitting, the network consists of three layers, i.e., input, hidden, and output layer (Figure 3).

Figure 3
www.frontiersin.org

Figure 3. The architecture of feed-forward ANN.

Initially, there were 10 neurons in a hidden layer which were modified as needed during the training stage. Scaled Conjugate Gradient (SCG) was selected as a training function. SCG is a fully automated and fast supervised learning algorithm (Møller, 1993). The SCG algorithm updates the weights W of the neural network (Eq. 1).

W k + 1 = W k α k f W k     (1)
Where W k + 1 is the weight matrix of the neural network at iteration k and α k is step size, which can be adaptively adjusted during the optimization process. The gradient f W k is typically computed using backpropagation, which efficiently calculates the gradient of the error function with respect to each weight in the network. The distribution of training, validation, and testing cohorts are 75, 15, and 15%, respectively. For progress monitoring and validation checks, the maximum number of training epochs used as stopping criteria as per mean squared error (MSE) is defined as follows (Eq. 2).
M S E = 1 N i = 1 i = N y i x i 2     (2)
Where x i is the input value, y i represents the predicted value and N is the total number of datasets. Error histograms and regression plots were used to determine the error distribution and accuracy of ANN at different stages.

2.6.1 Predictors

The input cohort used for the machine learning process was a combination of all eligible patients from Set-I and Set-II, i.e., N = 200. There were 28 patient-specific input variables majorly comprised of demographics, tumor stage, and site, RT-dose, and treatment planning parameters. Subclasses of predictors include the type of concurrent chemotherapy drug administered with RT, radiation dose indices to target and OAR, percentage volumes of parotids lying inside the target, and patients’ weight loss after RT.

2.6.2 Outcomes

The output layer consists of TCP values calculated via Poisson, Zaider-Minerbo, and EUD methods. NTCP calculated for both parotids by Lyman-Kutcher-Burman and EUD methods were also included in the output layer.

3 Results

ANN training begins without the hardware acceleration of the personal computer and it performs efficiently. The training and validation performance evaluated by MSE vs. no. of the epochs plot (Figure 4A). The plot shows a gradual decrease in MSE values versus the number of epochs. The progressive decline in MSE values represents the good accuracy of training and test data. Further, it shows the adjustment of the right weight after every epoch. The MSE reduces around 15 times during the training duration and after reaching the saturation the training stops at epoch number 38. It took 1.06 s to complete the training. The best validation performance was noted at epoch number 32 with an MSE value of 0.0465.

Figure 4
www.frontiersin.org

Figure 4. The performance plot (A) for ANN and error histogram (B) between predicted and target values.

An error histogram was plotted to determine the distribution of errors between the predicted and target values at different instances (Figure 4B). The histogram was also used as a measure of the prediction performance of ANN. Error-values were distributed in 20 bins and placed in columns. The error range for the model lies between −0.9961 and 0.9462. More than 75% of the data points from training, validation, and testing lie in the error range of −0.076 to 0.1258. After the validation checks, the prediction capacity of the model was accessed with 15% of the test cases (Figure 5). Linear regression plots between target and output showed for training, validation, and testing stages, represented by blue, green, and red lines, respectively. Each plot is supported by a line fit equation between the predicted and target values where the target is the independent variable and the output is a dependent variable. The proportionality between target and output is the coefficient of a target for which the value on the stage of training is 0.79 and 0.81, 0.71 for the validation and testing stages, respectively. The constants at the end of equations were errors or residuals and its incorporation scaled the target to make better predictions in output. The error values were 0.098, 0.1, and 0.14 for the training, validation, and testing stages. These equations highlight the performance of ANN. The regression coefficient (R) represents the accuracy of the model. At the final stage, the value of R found for all stages together was 0.877 and 0.89, 0.87, 0.82, respectively, for the training, validation, and testing stages.

Figure 5
www.frontiersin.org

Figure 5. The linear regression plots for the training, validation, and testing stages of the ANN.

4 Discussion

For the prediction of radiobiological indices, datasets were chosen with few special considerations for the development of ANN. (1) The patients were treated with similar conformal RT techniques and narrow integral and per fraction dose variations. (2) Both or ipsilateral parotid must be nearby or involved with the PTV. (3) DVH must be available for both parotids and PTV. The above principles help maintain consistency among input variables and result in smooth progression during the training stage.

The particular significance given to parotid glands as an OAR in this study of HNSCC patients is because even partial sparing of these glands will have a positive effect on patients’ quality of life (QOL) (Deasy et al., 2010). Our ANN model showed the capacity to predict NTCPs for parotid glands and TCPs for target structure. This tool will equip clinicians to make informed decisions with patient-specific information to adapt and modify OAR and PTV structures, which may result in better QOL for patients.

The inclusion of patient data from different institutions strengthens the model and incorporates the inter-variability of treatment planning and contouring practices among the institutes. It also enables the model to adapt further data sets from different regions and origins to enhance the learning and prediction process. Additionally, the architecture of ANN was designed in a way that makes it susceptible to adding more OARs from HNSCC patients.

The data sets used for training and validation were of conventional fractionation, ranging from 2 to 2.12 Gy. This limits the prediction of indices for hypofractionation plans. The non-availability of comprehensive data on the radiobiological effects of higher doses on OARs was one of the limiting factors. However, a substantial portion of HNSCC patients in developing and developed countries are still treated with conventional fractions, which highlights the significance of this prediction model.

Radiobiological indices were calculated from different models to complement the network with the latest available organ-specific radiobiological parameters for HNSCC RT. It enables the prediction model for current and upcoming HNSCC datasets to be tested and comprehended. In contrast to open source available calculators like RADBIOMOD and RBMODELV1 (Patel et al., 2022) which utilize tabular DVH values for the calculation of TCP and NTCP, the proposed ANN model predicts the radiobiological indices by incorporating significant parameters like tumor specifications, RT parameters, DVH and patient demographics. This makes the current model a more comprehensive predictor of TCP and NTCP.

The ANN-driven prediction model with its primary application to assist treatment planners and radiation oncologists with TCP and NTCP values, would potentially minimize the see-saw effect between coverage of PTV and sparing of parotids during RT treatment planning of HNSCC. This could empower radiation oncologists to optimize the margins of PTV near parotids in conjunction with their TCP and NTCP values. Another edge it could provide in adaptive RT by incorporating previous and current treatment plan parameters into the model and producing TCP and NTCP for insight.

For future work, the incorporation of more HNSCC OARs such as the spinal cord, cochlea, and optical structures from current datasets is in process. We are also seeking larger datasets with variable dose fascinations to scale up the performance of the current model.

5 Conclusion

In conclusion, our study has successfully developed ANN with a high prediction accuracy (R = 0.82) for radiobiological indices, offering a valuable tool for HNSCC treatment planning. The observed correlation between parotid volume within the PTV and increased NTCP underscores the importance of minimizing PTV margins to reduce post-irradiation complications. Furthermore, our research expanded predictive factors to include demographics, tumor characteristics, treatment parameters, and additional variables such as concurrent chemotherapy, radiation dose indices, and post-RT weight loss. These enhancements, combined with DVH, TCP, and NTCP models, contribute to the development of personalized and more accurate tools for toxicity prediction, advancing the goal of optimizing cancer care for HNSCC patients.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

SA: Data curation, Formal analysis, Investigation, Methodology, Validation, Writing – original draft, Writing – review & editing. SN: Data curation, Methodology, Writing – review & editing. AK: Data curation, Investigation, Methodology, Writing – review & editing. BQ: Data curation, Investigation, Methodology, Writing – review & editing. AH: Data curation, Formal analysis, Investigation, Methodology, Writing – review & editing. BA: Data curation, Formal analysis, Methodology, Writing – review & editing. WM: Conceptualization, Data curation, Funding acquisition, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. Research reported in this publication was supported by the Proton International under Award Number C-23-262. The content is solely the responsibility of the authors and does not necessarily represent the official views of the Proton International.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Abola, M. V., Prasad, V., and Jena, A. B. (2014). Association between treatment toxicity and outcomes in oncology clinical trials. Ann. Oncol. 25, 2284–2289. doi: 10.1093/annonc/mdu444

PubMed Abstract | Crossref Full Text | Google Scholar

Araújo, A. L. D., Moraes, M. C., Pérez-de-Oliveira, M. E., Silva, V. M. D., Saldivia-Siracusa, C., Pedroso, C. M., et al. (2023). Machine learning for the prediction of toxicities from head and neck cancer treatment: a systematic review with meta-analysis. Oral Oncol. 140:106386. doi: 10.1016/j.oraloncology.2023.106386

PubMed Abstract | Crossref Full Text | Google Scholar

Ataei, A., Deng, J., and Muhammad, W. (2023). Liver cancer risk quantification through an artificial neural network based on personal health data. Acta Oncol. 62, 495–502. doi: 10.1080/0284186X.2023.2213445

Crossref Full Text | Google Scholar

Bejarano, T., De Ornelas-Couto, M., and Mihaylov, I. B. (2019). Longitudinal fan-beam computed tomography dataset for head-and-neck squamous cell carcinoma patients. Med. Phys. 46, 2526–2537. doi: 10.1002/mp.13460

PubMed Abstract | Crossref Full Text | Google Scholar

Bufacchi, A., Caspiani, O., Rambaldi, G., Marmiroli, L., Giovinazzo, G., and Polsoni, M. (2020). Clinical implication in the use of the AAA algorithm versus the AXB in nasopharyngeal carcinomas by comparison of TCP and NTCP values. Radiat. Oncol. 15:150. doi: 10.1186/s13014-020-01591-7

PubMed Abstract | Crossref Full Text | Google Scholar

Chang, J. H., Gehrke, C., Prabhakar, R., Gill, S., Wada, M., Joon, D. L., et al. (2016). RADBIOMOD: a simple program for utilising biological modelling in radiotherapy plan evaluation. Phys. Med. 32, 248–254. doi: 10.1016/j.ejmp.2015.10.091

PubMed Abstract | Crossref Full Text | Google Scholar

Clark, K., Vendt, B., Smith, K., Freymann, J., Kirby, J., Koppel, P., et al. (2013). The Cancer imaging archive (TCIA): maintaining and operating a public information repository. J. Digit. Imaging 26, 1045–1057. doi: 10.1007/s10278-013-9622-7

PubMed Abstract | Crossref Full Text | Google Scholar

Deasy, J. O., Moiseenko, V., Marks, L., Chao, K. C., Nam, J., and Eisbruch, A. (2010). Radiotherapydose–volume effects on salivary gland function. International journal of radiation oncology* biology*. Physics 76, S58–S63. doi: 10.1016/j.ijrobp.2009.06.090

Crossref Full Text | Google Scholar

Dilalla, V., Chaput, G., Williams, T., and Sultanem, K. (2020). Radiotherapy side effects: integrating a survivorship clinical lens to better serve patients. Curr. Oncol. 27, 107–112. doi: 10.3747/co.27.6233

PubMed Abstract | Crossref Full Text | Google Scholar

Gardner, S. J., Kim, J., and Chetty, I. J. (2019). Modern radiation therapy planning and delivery. Hematol. Clin. 33, 947–962. doi: 10.1016/j.hoc.2019.08.005

Crossref Full Text | Google Scholar

Gay, H. A., and Niemierko, A. (2007). A free program for calculating EUD-based NTCP and TCP in external beam radiotherapy. Phys. Med. 23, 115–125. doi: 10.1016/j.ejmp.2007.07.001

PubMed Abstract | Crossref Full Text | Google Scholar

Giraud, P., Giraud, P., Gasnier, A., El Ayachy, R., Kreps, S., Foy, J. P., et al. (2019). Radiomics and machine learning for radiotherapy in head and neck cancers. Front. Oncol. 9:174. doi: 10.3389/fonc.2019.00174

PubMed Abstract | Crossref Full Text | Google Scholar

Gronberg, M. P., Beadle, B. M., Garden, A. S., Skinner, H., Gay, S., Netherton, T., et al. (2023). Deep learning–based dose prediction for automated, individualized quality Assurance of Head and Neck Radiation Therapy Plans. Pract. Radiat. Oncol. 13, e282–e291. doi: 10.1016/j.prro.2022.12.003

PubMed Abstract | Crossref Full Text | Google Scholar

Grossberg, A. J., Mohamed, A. S., Elhalawani, H., Bennett, W. C., Smith, K. E., Nolan, T. S., et al. (2018). Imaging and clinical data archive for head and neck squamous cell carcinoma patients treated with radiotherapy. Sci. Data 5:1. doi: 10.1038/s41597-018-0002-5

PubMed Abstract | Crossref Full Text | Google Scholar

Huang, S. H., and O’Sullivan, B. (2017). Overview of the 8th edition TNM classification for head and neck cancer. Curr. Treat. Options Oncol. 18:40. doi: 10.1007/s11864-017-0484-y

PubMed Abstract | Crossref Full Text | Google Scholar

Hunter, K. U., Schipper, M., Feng, F. Y., Lyden, T., Haxer, M., Murdoch-Kinch, C. A., et al. (2013). Toxicities affecting quality of life after chemo-IMRT of oropharyngeal Cancer: prospective study of patient-reported, observer-rated, and objective outcomes. Int. J. Radiat. Oncol. 85, 935–940. doi: 10.1016/j.ijrobp.2012.08.030

PubMed Abstract | Crossref Full Text | Google Scholar

Isaksson, L. J., Pepa, M., Zaffaroni, M., Marvaso, G., Alterio, D., Volpe, S., et al. (2020). Machine learning-based models for prediction of toxicity outcomes in radiotherapy. Front. Oncol. 10:790. doi: 10.3389/fonc.2020.00790

PubMed Abstract | Crossref Full Text | Google Scholar

Jackson, A., Ten Haken, R. K., Robertson, J. M., Kessler, M. L., Kutcher, G. J., and Lawrence, T. S. (1995). Analysis of clinical complication data for radiation hepatitis using a parallel architecture model. International journal of radiation oncology* biology*. Physics 31, 883–891. doi: 10.1016/0360-3016(94)00471-4

Crossref Full Text | Google Scholar

Johnson, D. E., Burtness, B., Leemans, C. R., Lui, V. W. Y., Bauman, J. E., and Grandis, J. R. (2020). Head and neck squamous cell carcinoma. Nat. Rev. Dis. Primer. 6:92. doi: 10.1038/s41572-020-00224-3

PubMed Abstract | Crossref Full Text | Google Scholar

Kawashita, Y., Soutome, S., Umeda, M., and Saito, T. (2020). Oral management strategies for radiotherapy of head and neck cancer. Jpn. Dent. Sci. Rev. 56, 62–67. doi: 10.1016/j.jdsr.2020.02.001

PubMed Abstract | Crossref Full Text | Google Scholar

Konings, H., Stappers, S., Geens, M., De Winter, B. Y., Lamote, K., Van Meerbeeck, J. P., et al. (2020). A literature review of the potential diagnostic biomarkers of head and neck neoplasms. Front. Oncol. 10:1020. doi: 10.3389/fonc.2020.01020

PubMed Abstract | Crossref Full Text | Google Scholar

Lyman, J. T. (1985). Complication probability as assessed from dose-volume histograms. Radiat. Res. 104, S13–S19. doi: 10.2307/3576626

Crossref Full Text | Google Scholar

Mahmood, H., Shaban, M., Rajpoot, N., and Khurram, S. A. (2021). Artificial intelligence-based methods in head and neck cancer diagnosis: an overview. Br. J. Cancer 124, 1934–1940. doi: 10.1038/s41416-021-01386-x

PubMed Abstract | Crossref Full Text | Google Scholar

Møller, M. F. (1993). A scaled conjugate gradient algorithm for fast supervised learning. Neural Netw. 6, 525–533. doi: 10.1016/S0893-6080(05)80056-5

Crossref Full Text | Google Scholar

Muhammad, W., Hart, G. R., Nartowt, B., Farrell, J. J., Johung, K., Liang, Y., et al. (2019). Pancreatic cancer prediction through an artificial neural network. Front. Artif. Intell. 2:2. doi: 10.3389/frai.2019.00002

PubMed Abstract | Crossref Full Text | Google Scholar

Parmar, C., Grossmann, P., Rietveld, D., Rietbergen, M. M., Lambin, P., and Aerts, H. J. W. L. (2015). Radiomic machine-learning classifiers for prognostic biomarkers of head and neck cancer. Front. Oncol. 5:272. doi: 10.3389/fonc.2015.00272

Crossref Full Text | Google Scholar

Patel, G., Mandal, A., Bharati, A., Choudhary, S., Mishra, R., and Mourya, A. (2022). Development and validation of an indigenous, radiobiological model-based tumor control probability and normal tissue complication probability estimation software for routine plan evaluation in clinics. J. Cancer Res. Ther. 18, 1697–1705. doi: 10.4103/jcrt.JCRT_330_20

PubMed Abstract | Crossref Full Text | Google Scholar

Peng, Z., Wang, Y., Wang, Y., Jiang, S., Fan, R., Zhang, H., et al. (2021). Application of radiomics and machine learning in head and neck cancers. Int. J. Biol. Sci. 17, 475–486. doi: 10.7150/ijbs.55716

PubMed Abstract | Crossref Full Text | Google Scholar

Pereira, G. C., Traughber, M., and Muzic, R. F. (2014). The role of imaging in radiation therapy planning: past, present, and future. Biomed. Res. Int. 2014:e231090. doi: 10.1155/2014/231090

Crossref Full Text | Google Scholar

Shafiq, A., Çolak, A. B., Sindhu, T. N., Lone, S. A., Alsubie, A., and Jarad, F. (2022). Comparative study of artificial neural network versus parametric method in COVID-19 data analysis. Results Phys. 38:105613. doi: 10.1016/j.rinp.2022.105613

PubMed Abstract | Crossref Full Text | Google Scholar

Thor, M., Iyer, A., Jiang, J., Apte, A., Veeraraghavan, H., Allgood, N. B., et al. (2021). Deep learning auto-segmentation and automated treatment planning for trismus risk reduction in head and neck cancer radiotherapy. Phys. Imaging Radiat. Oncol. 19, 96–101. doi: 10.1016/j.phro.2021.07.009

PubMed Abstract | Crossref Full Text | Google Scholar

Warkentin, B., Stavrev, P., Stavreva, N., Field, C., and Fallone, B. G. (2004). A TCP-NTCP estimation module using DVHs and known radiobiological models and parameter sets. J. Appl. Clin. Med. Phys. 5, 50–63. doi: 10.1120/jacmp.v5i1.1970

PubMed Abstract | Crossref Full Text | Google Scholar

Webb, S., and Nahum, A. E. (1993). A model for calculating tumour control probability in radiotherapy including the effects of inhomogeneous distributions of dose and clonogenic cell density. Phys. Med. Biol. 38, 653–666. doi: 10.1088/0031-9155/38/6/001

PubMed Abstract | Crossref Full Text | Google Scholar

Yorke, E. (2023). Modeling clinical outcomes in radiotherapy: NTCP, TCP and the “TECs”. Med. Phys. 50, 122–124. doi: 10.1002/mp.16274

PubMed Abstract | Crossref Full Text | Google Scholar

Zaider, M., and Minerbo, G. N. (2000). Tumour control probability: a formulation applicable to any temporal protocol of dose delivery. Phys. Med. Biol. 45, 279–293. doi: 10.1088/0031-9155/45/2/303

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: tumor control probability, normal tissue complication probability, head and neck cancer, artificial neural network, radiation therapy

Citation: Ahmed SBS, Naeem S, Khan AMH, Qureshi BM, Hussain A, Aydogan B and Muhammad W (2024) Artificial neural network-assisted prediction of radiobiological indices in head and neck cancer. Front. Artif. Intell. 7:1329737. doi: 10.3389/frai.2024.1329737

Received: 01 November 2023; Accepted: 25 March 2024;
Published: 05 April 2024.

Edited by:

Tuan D. Pham, Queen Mary University of London, United Kingdom

Reviewed by:

Andaç Batur Çolak, Istanbul Commerce University, Türkiye
Hasi Hays, University of Arkansas, United States

Copyright © 2024 Ahmed, Naeem, Khan, Qureshi, Hussain, Aydogan and Muhammad. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Wazir Muhammad, wmuhammad@fau.edu; wazirm@hotmail.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.