Abstract:Purpose: With the coming era of digital medicine and healthcare technology, mathematical modeling of tumors has become a key step to optimize and realize precision radiation therapy. The purpose of this study was to develop a mathematical model for simulating the change of head and neck (HN) tumor volume during radiation therapy. Methods and materials: A formula was developed to describe the dynamic change of oxygenated compartment within a tumor, which was combined with the lethal lesions model to describe various cell processes during radiation therapy, including potentially lethal lesion repair and misrepair, cell proliferation/loss, and tumor reoxygenation. Parameter sensitivity analysis was performed to evaluate the impacts of lesion- and repair-related biological factors on radiation therapy outcomes. Results: We tested our model on 14 available patients with HN cancer and compared the performance with 3 other models. The mean error of our model for the 12 good fit cases was 12.2%, which is considerably smaller than that of the linear quadratic model (19.7%), the generalized linear quadratic model (19.1%), and a 4-level cell population model (16.6%). Correlation analysis results revealed that for small tumors, there was a positive correlation (correlation coefficient r=0.9416) between hypoxic fraction (hf) and tumor volume, whereas the correlation became negative and not significant (r=-0.4365) for large tumors. It is demonstrated from sensitivity analysis that the production rate of lethal lesions (ηl) has a far greater impact on tumor volume than other parameters. The hf had an insignificant impact on tumor volume but had a notable influence on the volume of surviving cells. The final volume of surviving cells athf=0.5 was almost 8 ×102 times that of hf=0.01. The potentially lethal lesion-related parameters (the production rate of potentially lethal lessions per unit dose ηpl, the rate of correct repair per unit time εpl, and the rate of binary misrepair per unit time ε2pl) had rather small impacts (<1%) on both tumor volume and the volume of surviving cells, which indicates that the repaired and misrepaired sublethal cells only take up a small portion of the total cancer cell population. Conclusions: A population-based tumor-volume model for HN cancer during radiation therapy with a dynamic oxygenated compartment was developed in this study. Comprehensively considering the damage process of tumor cells caused by radiation therapy, the accurate prediction of the volume change of HN tumors during treatment was revealed. Meanwhile, various cell activities and their principles in the process of antitumor treatment were reflected, which has positive clinical reference significance for radiobiology.

A Mathematical Model for Tumor Cell Population Dynamics Based on Target Theory and Tumor Lifespan

Simulating the Dynamic Intra-Tumor Heterogeneity and Therapeutic Responses

Impact of radiation therapy on healthy and cancerous cell dynamics: a Mathematical analysis

Stochastic Model for Tumor Control Probability: Effects of Cell Cycle and (A)symmetric Proliferation

Temporal optimization of radiation therapy to heterogeneous tumour populations and cancer stem cells

A Mathematical Model of Thyroid Disease Response to Radiotherapy

A Mathematical Tumor Model with Immune Resistance and Drug Therapy: An Optimal Control Approach

A multiscale mathematical model of cancer, and its use in analyzing irradiation therapies

Towards Predicting the Response of a Solid Tumour to Chemotherapy and Radiotherapy Treatments: Clinical Insights from a Computational Model

Modeling tumor growth: a simple individual-based model and its analysis

A Population-Based Tumor-Volume Model for Head and Neck Cancer During Radiation Therapy With a Dynamic Oxygenated Compartment

A biophysical model of cell evolution after cytotoxic treatments: damage, repair and cell response

A cellular automata model for avascular solid tumor growth under the effect of therapy

Mathematical modeling and dynamic analysis for cancer resistance incorporating persister cells

Spatio-temporal modelling of phenotypic heterogeneity in tumour tissues and its impact on radiotherapy treatment

Mathematical modeling and dynamic analysis of anti-tumor immune response

Mathematical Modeling of Tumors Growth: Competition based on Gompertz model in Two Dimensions

A kinetic model of continuous radiation damage to populations of cells: Comparison to the LQ model and application to molecular radiotherapy

Predicting Radiotherapy Patient Outcomes with Real-Time Clinical Data Using Mathematical Modelling

Tumor containment: a more general mathematical analysis

A mathematical model with aberrant growth correction in tissue homeostasis and tumor cell growth