Mathematical Modeling of the Effects of Tumor Heterogeneity on the Efficiency of Radiation Treatment Schedule

Model-based assessment of combination therapies - ranking of radiosensitizing agents in oncology

BACKGROUND

To increase the chances of finding efficacious anticancer drugs, improve development times and reduce costs, it is of interest to rank test compounds based on their potential for human use as early as possible in the drug development process. In this paper, we present a method for ranking radiosensitizers using preclinical data.

METHODS

We used data from three xenograft mice studies to calibrate a model that accounts for radiation treatment combined with radiosensitizers. A nonlinear mixed effects approach was utilized where between-subject variability and inter-study variability were considered. Using the calibrated model, we ranked three different Ataxia telangiectasia-mutated inhibitors in terms of anticancer activity. The ranking was based on the Tumor Static Exposure (TSE) concept and primarily illustrated through TSE-curves.

RESULTS

The model described data well and the predicted number of eradicated tumors was in good agreement with experimental data. The efficacy of the radiosensitizers was evaluated for the median individual and the 95% population percentile. Simulations predicted that a total dose of 220 Gy (5 radiation sessions a week for 6 weeks) was required for 95% of tumors to be eradicated when radiation was given alone. When radiation was combined with doses that achieved at least 8 [Formula: see text] of each radiosensitizer in mouse blood, it was predicted that the radiation dose could be decreased to 50, 65, and 100 Gy, respectively, while maintaining 95% eradication.

CONCLUSIONS

A simulation-based method for calculating TSE-curves was developed, which provides more accurate predictions of tumor eradication than earlier, analytically derived, TSE-curves. The tool we present can potentially be used for radiosensitizer selection before proceeding to subsequent phases of the drug discovery and development process.

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

External beam radiation therapy is a key part of modern cancer treatments which uses high doses of radiation to destroy tumour cells. Despite its widespread usage and extensive study in theoretical, experimental, and clinical works, many questions still remain about how best to administer it. Many mathematical studies have examined optimal scheduling of radiotherapy, and most come to similar conclusions. Importantly though, these studies generally assume intratumoral homogeneity. But in recent years, it has become clear that tumours are not homogeneous masses of cancerous cells, but wildly heterogeneous masses with various subpopulations which grow and respond to treatment differently. One subpopulation of particular importance is cancer stem cells (CSCs) which are known to exhibit higher radioresistence compared with non-CSCs. Knowledge of these differences between cell types could theoretically lead to changes in optimal treatment scheduling. Only a few studies have examined this question, and interestingly, they arrive at apparent conflicting results. However, an understanding of their assumptions reveals a key difference which leads to their differing conclusions. In this paper, we generalize the problem of temporal optimization of dose distribution of radiation therapy to a two cell type model. We do so by creating a mathematical model and a numerical optimization algorithm to find the distribution of dose which leads to optimal cell kill. We then create a data set of optimization solutions and use data analysis tools to learn the relationships between model parameters and the qualitative behaviour of optimization results. Analysis of the model and discussion of biological importance are provided throughout. We find that the key factor in predicting the behaviour of the optimal distribution of radiation is the ratio between the radiosensitivities of the present cell types. These results can provide guidance for treatment in cases where clinicians have knowledge of tumour heterogeneity and of the abundance of CSCs.

A Century of Fractionated Radiotherapy: How Mathematical Oncology Can Break the Rules

Radiotherapy is involved in 50% of all cancer treatments and 40% of cancer cures. Most of these treatments are delivered in fractions of equal doses of radiation (Fractional Equivalent Dosing (FED)) in days to weeks. This treatment paradigm has remained unchanged in the past century and does not account for the development of radioresistance during treatment. Even if under-optimized, deviating from a century of successful therapy delivered in FED can be difficult. One way of exploring the infinite space of fraction size and scheduling to identify optimal fractionation schedules is through mathematical oncology simulations that allow for in silico evaluation. This review article explores the evidence that current fractionation promotes the development of radioresistance, summarizes mathematical solutions to account for radioresistance, both in the curative and non-curative setting, and reviews current clinical data investigating non-FED fractionated radiotherapy.

The impact of random microenvironmental fluctuations on tumor control probability

The tumor control probability (TCP) is a metric used to calculate the probability of controlling or eradicating tumors through radiotherapy. Cancer cells vary in their response to radiation, and although many factors are involved, the tumor microenvironment is a crucial one that determines radiation efficacy. The tumor microenvironment plays a significant role in cancer initiation and propagation, as well as in treatment outcome. We have developed stochastic formulations to study the impact of arbitrary microenvironmental fluctuations on TCP and extinction probability (EP), which is defined as the probability of cancer cells removal in the absence of treatment. Since the derivation of analytical solutions are not possible for complicated cases, we employ a modified Gillespie algorithm to analyze TCP and EP, considering the random variations in cellular proliferation and death rates. Our results show that increasing the standard deviation in kinetic rates initially enhances the probability of tumor eradication. However, if the EP does not reach a probability of 1, the increase in the standard deviation subsequently has a negative impact on probability of cancer cells removal, decreasing the EP over time. The greatest effect on EP has been observed when both birth and death rates are being randomly modified and are anticorrelated. In addition, similar results are observed for TCP, where radiotherapy is included, indicating that increasing the standard deviation in kinetic rates at first enhances the probability of tumor eradication. But, it has a negative impact on treatment effectiveness if the TCP does not reach a probability of 1.

Cellular heterogeneity landscape in laryngeal squamous cell carcinoma

Laryngeal squamous cell carcinoma (LSCC) is a highly malignant tumor originated from respiratory system. Although there have been many improvements in therapy until now, reducing the high mortality remains difficult. Understanding the cellular heterogeneity of LSCC could contribute to improve this problem. Single-cell RNA sequencing was applied to dissect the cell composition and molecular characteristics of LSCC tissues. Immunohistochemistry staining of the LSCC tissues was performed to identify the spatial location of tumor cells. Survival analysis of marker genes was executed in The Cancer Genome Atlas to verify the correlation between each cell clusters and patients' prognosis. The LSCC tissue cells were finely grouped into various clusters, including tumor cells, immune cells, epithelial cells, fibroblasts and endothelial cells. Notably, in tumor cells, keratinocyte-like cells were in the core of tumor while malignant proliferating cells were located at the tumor edge. The malignant proliferating cells were correlated with poor prognosis. In summary, this is the first study to delineate a landscape of the LSCC intratumor heterogeneity. Our work might help researchers have a better understanding for tumor progression.

Three-dimensional tumor growth in time-varying chemical fields: a modeling framework and theoretical study

Background

Contemporary biological observations have revealed a large variety of mechanisms acting during the expansion of a tumor. However, there are still many qualitative and quantitative aspects of the phenomenon that remain largely unknown. In this context, mathematical and computational modeling appears as an invaluable tool providing the means for conducting in silico experiments, which are cheaper and less tedious than real laboratory experiments.

Results

This paper aims at developing an extensible and computationally efficient framework for in silico modeling of tumor growth in a 3-dimensional, inhomogeneous and time-varying chemical environment. The resulting model consists of a set of mathematically derived and algorithmically defined operators, each one addressing the effects of a particular biological mechanism on the state of the system. These operators may be extended or re-adjusted, in case a different set of starting assumptions or a different simulation scenario needs to be considered.

Conclusion

In silico modeling provides an alternative means for testing hypotheses and simulating scenarios for which exact biological knowledge remains elusive. However, finer tuning of pertinent methods presupposes qualitative and quantitative enrichment of available biological evidence. Validation in a strict sense would further require comprehensive, case-specific simulations and detailed comparisons with biomedical observations.

Illuminating the Numbers: Integrating Mathematical Models to Optimize Photomedicine Dosimetry and Combination Therapies

Cancer photomedicine offers unique mechanisms for inducing local tumor damage with the potential to stimulate local and systemic anti-tumor immunity. Optically-active nanomedicine offers these features as well as spatiotemporal control of tumor-focused drug release to realize synergistic combination therapies. Achieving quantitative dosimetry is a major challenge, and dosimetry is fundamental to photomedicine for personalizing and tailoring therapeutic regimens to specific patients and anatomical locations. The challenge of dosimetry is perhaps greater for photomedicine than many standard therapies given the complexity of light delivery and light-tissue interactions as well as the resulting photochemistry responsible for tumor damage and drug-release, in addition to the usual intricacies of therapeutic agent delivery. An emerging multidisciplinary approach in oncology utilizes mathematical and computational models to iteratively and quantitively analyze complex dosimetry, and biological response parameters. These models are parameterized by preclinical and clinical observations and then tested against previously unseen data. Such calibrated and validated models can be deployed to simulate treatment doses, protocols, and combinations that have not yet been experimentally or clinically evaluated and can provide testable optimal treatment outcomes in a practical workflow. Here, we foresee the utility of these computational approaches to guide adaptive therapy, and how mathematical models might be further developed and integrated as a novel methodology to guide precision photomedicine.