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

Mathematical modeling of brain metastases growth and response to therapies: A review

Brain metastases (BMs) are the most common intracranial tumor type and a significant health concern, affecting approximately 10% to 30% of all oncological patients. Although significant progress is being made, many aspects of the metastatic process to the brain and the growth of the resulting lesions are still not well understood. There is a need for an improved understanding of the growth dynamics and the response to treatment of these tumors. Mathematical models have been proven valuable for drawing inferences and making predictions in different fields of cancer research, but few mathematical works have considered BMs. This comprehensive review aims to establish a unified platform and contribute to fostering emerging efforts dedicated to enhancing our mathematical understanding of this intricate and challenging disease. We focus on the progress made in the initial stages of mathematical modeling research regarding BMs and the significant insights gained from such studies. We also explore the vital role of mathematical modeling in predicting treatment outcomes and enhancing the quality of clinical decision-making for patients facing BMs.

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.

A Mathematical Study of the Influence of Hypoxia and Acidity on the Evolutionary Dynamics of Cancer

Hypoxia and acidity act as environmental stressors promoting selection for cancer cells with a more aggressive phenotype. As a result, a deeper theoretical understanding of the spatio-temporal processes that drive the adaptation of tumour cells to hypoxic and acidic microenvironments may open up new avenues of research in oncology and cancer treatment. We present a mathematical model to study the influence of hypoxia and acidity on the evolutionary dynamics of cancer cells in vascularised tumours. The model is formulated as a system of partial integro-differential equations that describe the phenotypic evolution of cancer cells in response to dynamic variations in the spatial distribution of three abiotic factors that are key players in tumour metabolism: oxygen, glucose and lactate. The results of numerical simulations of a calibrated version of the model based on real data recapitulate the eco-evolutionary spatial dynamics of tumour cells and their adaptation to hypoxic and acidic microenvironments. Moreover, such results demonstrate how nonlinear interactions between tumour cells and abiotic factors can lead to the formation of environmental gradients which select for cells with phenotypic characteristics that vary with distance from intra-tumour blood vessels, thus promoting the emergence of intra-tumour phenotypic heterogeneity. Finally, our theoretical findings reconcile the conclusions of earlier studies by showing that the order in which resistance to hypoxia and resistance to acidity arise in tumours depend on the ways in which oxygen and lactate act as environmental stressors in the evolutionary dynamics of cancer cells.