Why one-size-fits-all vaso-modulatory interventions fail to control glioma invasion: in silico insights

There is an ongoing debate on the therapeutic potential of vaso-modulatory interventions against glioma invasion. Prominent vasculature-targeting therapies involve functional tumour-associated blood vessel deterioration and normalisation. The former aims at tumour infarction and nutrient deprivation medi- ated by vascular targeting agents that induce occlusion/collapse of tumour blood vessels. In contrast, the therapeutic intention of normalising the abnormal structure and function of tumour vascular net- works, e.g. via alleviating stress-induced vaso-occlusion, is to improve chemo-, immuno- and radiation therapy efficacy. Although both strategies have shown therapeutic potential, it remains unclear why they often fail to control glioma invasion into the surrounding healthy brain tissue. To shed light on this issue, we propose a mathematical model of glioma invasion focusing on the interplay between the mi- gration/proliferation dichotomy (Go-or-Grow) of glioma cells and modulations of the functional tumour vasculature. Vaso-modulatory interventions are modelled by varying the degree of vaso-occlusion. We discovered the existence of a critical cell proliferation/diffusion ratio that separates glioma invasion re- sponses to vaso-modulatory interventions into two distinct regimes. While for tumours, belonging to one regime, vascular modulations reduce the tumour front speed and increase the infiltration width, for those in the other regime the invasion speed increases and infiltration width decreases. We show how these in silico findings can be used to guide individualised approaches of vaso-modulatory treatment strategies and thereby improve success rates

Hypoxic Cell Waves around Necrotic Cores in Glioblastoma: A Biomathematical Model and its Therapeutic Implications

Glioblastoma is a rapidly evolving high-grade astrocytoma that is distinguished pathologically from lower grade gliomas by the presence of necrosis and microvascular hiperplasia. Necrotic areas are typically surrounded by hypercellular regions known as "pseudopalisades" originated by local tumor vessel occlusions that induce collective cellular migration events. This leads to the formation of waves of tumor cells actively migrating away from central hypoxia. We present a mathematical model that incorporates the interplay among two tumor cell phenotypes, a necrotic core and the oxygen distribution. Our simulations reveal the formation of a traveling wave of tumor cells that reproduces the observed histologic patterns of pseudopalisades. Additional simulations of the model equations show that preventing the collapse of tumor microvessels leads to slower glioma invasion, a fact that might be exploited for therapeutic purposes.Comment: 29 pages, 9 figure

Modeling Three-dimensional Invasive Solid Tumor Growth in Heterogeneous Microenvironment under Chemotherapy

A systematic understanding of the evolution and growth dynamics of invasive solid tumors in response to different chemotherapy strategies is crucial for the development of individually optimized oncotherapy. Here, we develop a hybrid three-dimensional (3D) computational model that integrates pharmacokinetic model, continuum diffusion-reaction model and discrete cell automaton model to investigate 3D invasive solid tumor growth in heterogeneous microenvironment under chemotherapy. Specifically, we consider the effects of heterogeneous environment on drug diffusion, tumor growth, invasion and the drug-tumor interaction on individual cell level. We employ the hybrid model to investigate the evolution and growth dynamics of avascular invasive solid tumors under different chemotherapy strategies. Our simulations reproduce the well-established observation that constant dosing is generally more effective in suppressing primary tumor growth than periodic dosing, due to the resulting continuous high drug concentration. In highly heterogeneous microenvironment, the malignancy of the tumor is significantly enhanced, leading to inefficiency of chemotherapies. The effects of geometrically-confined microenvironment and non-uniform drug dosing are also investigated. Our computational model, when supplemented with sufficient clinical data, could eventually lead to the development of efficient in silico tools for prognosis and treatment strategy optimization.Comment: 41 pages, 8 figure

Modelling the response of vascular tumours to chemotherapy: A multiscale approach

An existing multiscale model is extended to study the response of a vascularised tumour to treatment with chemotherapeutic drugs which target proliferating cells. The underlying hybrid cellular automaton model couples tissue-level processes (e.g. blood flow, vascular adaptation, oxygen and drug transport) with cellular and subcellular phenomena (e.g. competition for space, progress through the cell cycle, natural cell death and drug-induced cell kill and the expression of angiogenic factors). New simulations suggest that, in the absence of therapy, vascular adaptation induced by angiogenic factors can stimulate spatio-temporal oscillations in the tumour's composition.\ud \ud Numerical simulations are presented and show that, depending on the choice of model parameters, when a drug which kills proliferating cells is continuously infused through the vasculature, three cases may arise: the tumour is eliminated by the drug; the tumour continues to expand into the normal tissue; or, the tumour undergoes spatio-temporal oscillations, with regions of high vascular and tumour cell density alternating with regions of low vascular and tumour cell density. The implications of these results and possible directions for future research are also discussed

Computer simulation of glioma growth and morphology

Despite major advances in the study of glioma, the quantitative links between intra-tumor molecular/cellular properties, clinically observable properties such as morphology, and critical tumor behaviors such as growth and invasiveness remain unclear, hampering more effective coupling of tumor physical characteristics with implications for prognosis and therapy. Although molecular biology, histopathology, and radiological imaging are employed in this endeavor, studies are severely challenged by the multitude of different physical scales involved in tumor growth, i.e., from molecular nanoscale to cell microscale and finally to tissue centimeter scale. Consequently, it is often difficult to determine the underlying dynamics across dimensions. New techniques are needed to tackle these issues. Here, we address this multi-scalar problem by employing a novel predictive three-dimensional mathematical and computational model based on first-principle equations (conservation laws of physics) that describe mathematically the diffusion of cell substrates and other processes determining tumor mass growth and invasion. The model uses conserved variables to represent known determinants of glioma behavior, e.g., cell density and oxygen concentration, as well as biological functional relationships and parameters linking phenomena at different scales whose specific forms and values are hypothesized and calculated based on in vitro and in vivo experiments and from histopathology of tissue specimens from human gliomas. This model enables correlation of glioma morphology to tumor growth by quantifying interdependence of tumor mass on the microenvironment (e.g., hypoxia, tissue disruption) and on the cellular phenotypes (e.g., mitosis and apoptosis rates, cell adhesion strength). Once functional relationships between variables and associated parameter values have been informed, e.g., from histopathology or intra-operative analysis, this model can be used for disease diagnosis/prognosis, hypothesis testing, and to guide surgery and therapy. In particular, this tool identifies and quantifies the effects of vascularization and other cell-scale glioma morphological characteristics as predictors of tumor-scale growth and invasion