The importance of dead material within a tumour on the dynamics in response to radiotherapy

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

Longitudinal tumour volume data from head-and-neck cancer patients show that tumours of comparable pre-treatment size and stage may respond very differently to the same radiotherapy fractionation protocol. Mathematical models are often proposed to predict treatment outcome in this context, and have the potential to guide clinical decision-making and inform personalised fractionation protocols. Hindering effective use of models in this context is the sparsity of clinical measurements juxtaposed with the model complexity required to produce the full range of possible patient responses. In this work, we present a compartment model of tumour volume and tumour composition, which, despite relative simplicity, is capable of producing a wide range of patient responses. We then develop novel statistical methodology and leverage a cohort of existing clinical data to produce a predictive model of both tumour volume progression and the associated level of uncertainty that evolves throughout a patient's course of treatment. To capture inter-patient variability, all model parameters are patient specific, with a bootstrap particle filter-like Bayesian approach developed to model a set of training data as prior knowledge. We validate our approach against a subset of unseen data, and demonstrate both the predictive ability of our trained model and its limitations.

Simulating tumor volume dynamics in response to radiotherapy: Implications of model selection

From the beginning of the usage of radiotherapy (RT) for cancer treatment, mathematical modeling has been integral to understanding radiobiology and for designing treatment approaches and schedules. There has been extensive modeling of response to RT with the inclusion of various degrees of biological complexity. In this study, we compare three models of tumor volume dynamics: (1) exponential growth with RT directly reducing tumor volume, (2) logistic growth with direct tumor volume reduction, and (3) logistic growth with RT reducing the tumor carrying capacity with the objective of understanding the implications of model selection and informing the process of model calibration and parameterization. For all three models, we: examined the rates of change in tumor volume during and RT treatment course; performed parameter sensitivity and identifiability analyses; and investigated the impact of the parameter sensitivity on the tumor volume trajectories. In examining the tumor volume dynamics trends, we coined a new metric - the point of maximum reduction of tumor volume (MRV) - to quantify the magnitude and timing of the expected largest impact of RT during a treatment course. We found distinct timing differences in MRV, dependent on model selection. The parameter identifiability and sensitivity analyses revealed the interdependence of the different model parameters and that it is only possible to independently identify tumor growth and radiation response parameters if the underlying tumor growth rate is sufficiently large. Ultimately, the results of these analyses help us to better understand the implications of model selection while simultaneously generating falsifiable hypotheses about MRV timing that can be tested on longitudinal measurements of tumor volume from pre-clinical or clinical data with high acquisition frequency. Although, our study only compares three particular models, the results demonstrate that caution is necessary in selecting models of response to RT, given the artifacts imposed by each model.

Modelling the Tumour Microenvironment, but What Exactly Do We Mean by "Model"?

The Oxford English Dictionary includes 17 definitions for the word "model" as a noun and another 11 as a verb. Therefore, context is necessary to understand the meaning of the word model. For instance, "model railways" refer to replicas of railways and trains at a smaller scale and a "model student" refers to an exemplary individual. In some cases, a specific context, like cancer research, may not be sufficient to provide one specific meaning for model. Even if the context is narrowed, specifically, to research related to the tumour microenvironment, "model" can be understood in a wide variety of ways, from an animal model to a mathematical expression. This paper presents a review of different "models" of the tumour microenvironment, as grouped by different definitions of the word into four categories: model organisms, in vitro models, mathematical models and computational models. Then, the frequencies of different meanings of the word "model" related to the tumour microenvironment are measured from numbers of entries in the MEDLINE database of the United States National Library of Medicine at the National Institutes of Health. The frequencies of the main components of the microenvironment and the organ-related cancers modelled are also assessed quantitatively with specific keywords. Whilst animal models, particularly xenografts and mouse models, are the most commonly used "models", the number of these entries has been slowly decreasing. Mathematical models, as well as prognostic and risk models, follow in frequency, and these have been growing in use.

An integrated and fast imaging quality assurance phantom for a 0.35 T magnetic resonance imaging linear accelerator

Purpose

Periodic imaging quality assurance (QA) of magnetic resonance imaging linear accelerator (MRL) is critical. The feasibility of a new MRL imaging phantom used for QA in the low field was evaluated with automated image analysis of various parameters for accuracy and reproducibility.

Methods and materials

The new MRL imaging phantom was scanned across every 30 degrees of the gantry, having the on/off state of the linac in a low-field MRL system using three magnetic resonance imaging sequences: true fast imaging with steady-state precession (TrueFISP), T1 weighted (T1W), and T2 weighted (T2W). The DICOM files were used to calculate the imaging parameters: geometric distortion, uniformity, resolution, signal-to-noise ratio (SNR), and laser alignment. The point spread function (PSF) and edge spread function (ESF) were also calculated for resolution analysis.

Results

The phantom data showed a small standard deviation - and high consistency for each imaging parameter. The highest variability in data was observed with the true fast imaging sequence at the calibration angle, which was expected because of low resolution and short scan time (25 sec). The mean magnitude of the largest distortion measured within 200 mm diameter with TrueFISP was 0.31 ± 0.05 mm. The PSF, ESF, signal uniformity, and SNR measurements remained consistent. Laser alignment traditional offsets and angular deviation remained consistent.

Conclusions

The new MRL imaging phantom is reliable, reproducible, time effective, and easy to use for a 0.35 T MRL system. The results promise a more streamlined, time-saving, and error-free QA process for low-field MRL adapted in our clinical setting.

Mathematical modelling of the dynamics of image-informed tumor habitats in a murine model of glioma

Tumors exhibit high molecular, phenotypic, and physiological heterogeneity. In this effort, we employ quantitative magnetic resonance imaging (MRI) data to capture this heterogeneity through imaging-based subregions or "habitats" in a murine model of glioma. We then demonstrate the ability to model and predict the growth of the habitats using coupled ordinary differential equations (ODEs) in the presence and absence of radiotherapy. Female Wistar rats (N = 21) were inoculated intracranially with 10⁶ C6 glioma cells, a subset of which received 20 Gy (N = 5) or 40 Gy (N = 8) of radiation. All rats underwent diffusion-weighted and dynamic contrast-enhanced MRI at up to seven time points. All MRI data at each visit were subsequently clustered using k-means to identify physiological tumor habitats. A family of four models consisting of three coupled ODEs were developed and calibrated to the habitat time series of control and treated rats and evaluated for predictive capability. The Akaike Information Criterion was used for model selection, and the normalized sum-of-square-error (SSE) was used to evaluate goodness-of-fit in model calibration and prediction. Three tumor habitats with significantly different imaging data characteristics (p < 0.05) were identified: high-vascularity high-cellularity, low-vascularity high-cellularity, and low-vascularity low-cellularity. Model selection resulted in a five-parameter model whose predictions of habitat dynamics yielded SSEs that were similar to the SSEs from the calibrated model. It is thus feasible to mathematically describe habitat dynamics in a preclinical model of glioma using biology-based ODEs, showing promise for forecasting heterogeneous tumor behavior.

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

We present a mathematical model that describes how tumour heterogeneity evolves in a tissue slice that is oxygenated by a single blood vessel. Phenotype is identified with the stemness level of a cell and determines its proliferative capacity, apoptosis propensity and response to treatment. Our study is based on numerical bifurcation analysis and dynamical simulations of a system of coupled, non-local (in phenotypic "space") partial differential equations that link the phenotypic evolution of the tumour cells to local tissue oxygen levels. In our formulation, we consider a 1D geometry where oxygen is supplied by a blood vessel located on the domain boundary and consumed by the tumour cells as it diffuses through the tissue. For biologically relevant parameter values, the system exhibits multiple steady states; in particular, depending on the initial conditions, the tumour is either eliminated ("tumour-extinction") or it persists ("tumour-invasion"). We conclude by using the model to investigate tumour responses to radiotherapy, and focus on identifying radiotherapy strategies which can eliminate the tumour. Numerical simulations reveal how phenotypic heterogeneity evolves during treatment and highlight the critical role of tissue oxygen levels on the efficacy of radiation protocols that are commonly used in the clinic.

Mathematical modeling of radiotherapy and its impact on tumor interactions with the immune system

Radiotherapy is a primary therapeutic modality widely utilized with curative intent. Traditionally tumor response was hypothesized to be due to high levels of cell death induced by irreparable DNA damage. However, the immunomodulatory aspect of radiation is now widely accepted. As such, interest into the combination of radiotherapy and immunotherapy is increasing, the synergy of which has the potential to improve tumor regression beyond that observed after either treatment alone. However, questions regarding the timing (sequential vs concurrent) and dose fractionation (hyper-, standard-, or hypo-fractionation) that result in improved anti-tumor immune responses, and thus potentially enhanced tumor inhibition, remain. Here we discuss the biological response to radiotherapy and its immunomodulatory properties before giving an overview of pre-clinical data and clinical trials concerned with answering these questions. Finally, we review published mathematical models of the impact of radiotherapy on tumor-immune interactions. Ranging from considering the impact of properties of the tumor microenvironment on the induction of anti-tumor responses, to the impact of choice of radiation site in the setting of metastatic disease, these models all have an underlying feature in common: the push towards personalized therapy.

Tumor Volume Regression during and after Radiochemotherapy: A Macroscopic Description

Tumor volume regression during and after chemo and radio therapy is a useful information for clinical decisions. Indeed, a quantitative, patient oriented, description of the response to treatment can guide towards the modification of the scheduled doses or the evaluation of the best time for surgery. We propose a macroscopic algorithm which permits to follow quantitatively the time evolution of the tumor volume during and after radiochemotherapy. The method, initially validated with different cell-lines implanted in mice, is then successfully applied to the available data for partially responding and complete recovery patients.

Simulation of Gamma-Ray Transmission Buildup Factors for Stratified Spherical Layers

Deterministic particle transport codes usually take into account scattered photons with correct attenuation laws and application of buildup factor to incident beam. Transmission buildup factors for adipose, bone, muscle, and skin human tissues, as well as for various combinations of these media for point isotropic photon source with energies of .15, 1.5 and 15 MeV, for different thickness of layers, were carried out using Geant4 (version 10.5) simulation toolkit. Also, we performed the analysis of existing multilayered shield fitting models (Lin and Jiang, Kalos, Burke and Beck) of buildup factor and the proposition of a new model. We found that the model combining those of Burke and Beck, for low atomic number (Z) followed by high Z materials and Kalos 1 for high Z followed by low Z materials, accurately reproduces simulation results with approximated deviation of 3 ± 3%, 2 ± 2%, and 3 ± 2% for 2, 3, and 4 layers, respectively. Since buildup factors are the key parameter for point kernel calculations, a correct study can be of great interest to the large community of radiation physicists, in general, and to medical imaging and radiotreatment physicists, especially.

Forecasting Individual Patient Response to Radiation Therapy in Head and Neck Cancer With a Dynamic Carrying Capacity Model

Purpose:

To model and predict individual patient responses to radiation therapy.

Methods and Materials:

We modeled tumor dynamics as logistic growth and the effect of radiation as a reduction in the tumor carrying capacity, motivated by the effect of radiation on the tumor microenvironment. The model was assessed on weekly tumor volume data collected for 2 independent cohorts of patients with head and neck cancer from the H. Lee Moffitt Cancer Center (MCC) and the MD Anderson Cancer Center (MDACC) who received 66 to 70 Gy in standard daily fractions or with accelerated fractionation. To predict response to radiation therapy for individual patients, we developed a new forecasting framework that combined the learned tumor growth rate and carrying capacity reduction fraction (δ) distribution with weekly measurements of tumor volume reduction for a given test patient to estimate δ, which was used to predict patient-specific outcomes.

Results:

The model fit data from MCC with high accuracy with patient-specific δ and a fixed tumor growth rate across all patients. The model fit data from an independent cohort from MDACC with comparable accuracy using the tumor growth rate learned from the MCC cohort, showing transferability of the growth rate. The forecasting framework predicted patient-specific outcomes with 76% sensitivity and 83% specificity for locoregional control and 68% sensitivity and 85% specificity for disease-free survival with the inclusion of 4 on-treatment tumor volume measurements.

Conclusions:

These results demonstrate that our simple mathematical model can describe a variety of tumor volume dynamics. Furthermore, combining historically observed patient responses with a few patient-specific tumor volume measurements allowed for the accurate prediction of patient outcomes, which may inform treatment adaptation and personalization.

A Mathematical Model of Breast Tumor Progression Based on Immune Infiltration

Breast cancer is the most prominent type of cancer among women. Understanding the microenvironment of breast cancer and the interactions between cells and cytokines will lead to better treatment approaches for patients. In this study, we developed a data-driven mathematical model to investigate the dynamics of key cells and cytokines involved in breast cancer development. We used gene expression profiles of tumors to estimate the relative abundance of each immune cell and group patients based on their immune patterns. Dynamical results show the complex interplay between cells and molecules, and sensitivity analysis emphasizes the direct effects of macrophages and adipocytes on cancer cell growth. In addition, we observed the dual effect of IFN-γ on cancer proliferation, either through direct inhibition of cancer cells or by increasing the cytotoxicity of CD8+ T-cells.

Phenotypic variation modulates the growth dynamics and response to radiotherapy of solid tumours under normoxia and hypoxia

In cancer, treatment failure and disease recurrence have been associated with small subpopulations of cancer cells with a stem-like phenotype. In this paper, we develop and investigate a phenotype-structured model of solid tumour growth in which cells are structured by a stemness level, which varies continuously between stem-like and terminally differentiated behaviours. Cell evolution is driven by proliferation and death, as well as advection and diffusion with respect to the stemness structure variable. Here, the magnitude and sign of the advective flux are allowed to vary with the oxygen level. We use the model to investigate how the environment, in particular oxygen levels, affects the tumour's population dynamics and composition, and its response to radiotherapy. We use a combination of numerical and analytical techniques to quantify how under physiological oxygen levels the cells evolve to a differentiated phenotype and under low oxygen level (i.e., hypoxia) they de-differentiate. Under normoxia, the proportion of cancer stem cells is typically negligible and the tumour may ultimately become extinct whereas under hypoxia cancer stem cells comprise a dominant proportion of the tumour volume, enhancing radio-resistance and favouring the tumour's long-term survival. We then investigate how such phenotypic heterogeneity impacts the tumour's response to treatment with radiotherapy under normoxia and hypoxia. Of particular interest is establishing how the presence of radio-resistant cancer stem cells can facilitate a tumour's regrowth following radiotherapy. We also use the model to show how radiation-induced changes in tumour oxygen levels can give rise to complex re-growth dynamics. For example, transient periods of hypoxia induced by damage to tumour blood vessels may rescue the cancer cell population from extinction and drive secondary regrowth.

Data Driven Mathematical Model of Colon Cancer Progression

Every colon cancer has its own unique characteristics, and therefore may respond differently to identical treatments. Here, we develop a data driven mathematical model for the interaction network of key components of immune microenvironment in colon cancer. We estimate the relative abundance of each immune cell from gene expression profiles of tumors, and group patients based on their immune patterns. Then we compare the tumor sensitivity and progression in each of these groups of patients, and observe differences in the patterns of tumor growth between the groups. For instance, in tumors with a smaller density of naive macrophages than activated macrophages, a higher activation rate of macrophages leads to an increase in cancer cell density, demonstrating a negative effect of macrophages. Other tumors however, exhibit an opposite trend, showing a positive effect of macrophages in controlling tumor size. Although the results indicate that for all patients the size of the tumor is sensitive to the parameters related to macrophages, such as their activation and death rate, this research demonstrates that no single biomarker could predict the dynamics of tumors.

Bayesian Information-Theoretic Calibration of Radiotherapy Sensitivity Parameters for Informing Effective Scanning Protocols in Cancer

With new advancements in technology, it is now possible to collect data for a variety of different metrics describing tumor growth, including tumor volume, composition, and vascularity, among others. For any proposed model of tumor growth and treatment, we observe large variability among individual patients’ parameter values, particularly those relating to treatment response; thus, exploiting the use of these various metrics for model calibration can be helpful to infer such patient-specific parameters both accurately and early, so that treatment protocols can be adjusted mid-course for maximum efficacy. However, taking measurements can be costly and invasive, limiting clinicians to a sparse collection schedule. As such, the determination of optimal times and metrics for which to collect data in order to best inform proper treatment protocols could be of great assistance to clinicians. In this investigation, we employ a Bayesian information-theoretic calibration protocol for experimental design in order to identify the optimal times at which to collect data for informing treatment parameters. Within this procedure, data collection times are chosen sequentially to maximize the reduction in parameter uncertainty with each added measurement, ensuring that a budget of n high-fidelity experimental measurements results in maximum information gain about the low-fidelity model parameter values. In addition to investigating the optimal temporal pattern for data collection, we also develop a framework for deciding which metrics should be utilized at each data collection point. We illustrate this framework with a variety of toy examples, each utilizing a radiotherapy treatment regimen. For each scenario, we analyze the dependence of the predictive power of the low-fidelity model upon the measurement budget.