Gliomas are primary brain tumors with a high invasive potential and infiltrative spread. Among them, glioblastoma multiforme (GBM) exhibits microvascular hyperplasia and pronounced necrosis triggered by hypoxia. Histological samples showing garland-like hypercellular structures (so-called pseudopalisades) centered around the occlusion site of a capillary are typical for GBM and hint on poor prognosis of patient survival. We propose a multiscale modeling approach in the kinetic theory of active particles framework and deduce by an upscaling process a reaction-diffusion model with repellent pH-taxis. We prove existence of a unique global bounded classical solution for a version of the obtained macroscopic system and investigate the asymptotic behavior of the solution. Moreover, we study two different types of scaling and compare the behavior of the obtained macroscopic PDEs by way of simulations. These show that patterns (not necessarily of Turing type), including pseudopalisades, can be formed for some parameter ranges, in accordance with the tumor grade. This is true when the PDEs are obtained via parabolic scaling (undirected tissue), while no such patterns are observed for the PDEs arising by a hyperbolic limit (directed tissue). This suggests that brain tissue might be undirected - at least as far as glioma migration is concerned. We also investigate two different ways of including cell level descriptions of response to hypoxia and the way they are related .

We propose a model for glioma patterns in a microlocal tumor environment under the influence of acidity, angiogenesis, and tissue anisotropy. The bottom-up model deduction eventually leads to a system of reaction–diffusion–taxis equations for glioma and endothelial cell population densities, of which the former infers flux limitation both in the self-diffusion and taxis terms. The model extends a recently introduced (Kumar, Li and Surulescu, 2020) description of glioma pseudopalisade formation with the aim of studying the effect of hypoxia-induced tumor vascularization on the establishment and maintenance of these histological patterns which are typical for high-grade brain cancer. Numerical simulations of the population level dynamics are performed to investigate several model scenarios containing this and further effects.

Objective. The distribution of hypoxia within tissues plays a critical role in tumor diagnosis and prognosis. Recognizing the significance of tumor oxygenation and hypoxia gradients, we introduce mathematical frameworks grounded in mechanistic modeling approaches for their quantitative assessment within a tumor microenvironment. By utilizing known blood vasculature, we aim to predict hypoxia levels across different tumor types. Approach. Our approach offers a computational method to measure and predict hypoxia using known blood vasculature. By formulating a reaction-diffusion model for oxygen distribution, we derive the corresponding hypoxia profile. Main results. The framework successfully replicates observed inter- and intra-tumor heterogeneity in experimentally obtained hypoxia profiles across various tumor types (breast, ovarian, pancreatic). Additionally, we propose a data-driven method to deduce partial differential equation models with spatially dependent parameters, which allows us to comprehend the variability of hypoxia profiles within tissues. The versatility of our framework lies in capturing diverse and dynamic behaviors of tumor oxygenation, as well as categorizing states of vascularization based on the dynamics of oxygen molecules, as identified by the model parameters. Significance. The proposed data-informed mechanistic method quantitatively assesses hypoxia in the tumor microenvironment by integrating diverse histopathological data and making predictions across different types of data. The framework provides valuable insights from both modeling and biological perspectives, advancing our comprehension of spatio-temporal dynamics of tumor oxygenation.

We propose a multiphase modeling approach to describe glioma pseudopalisade patterning under the influence of acidosis. The phases considered at the model onset are glioma, normal tissue, necrotic matter, and interstitial fluid in a void-free volume with acidity represented by proton concentration. We start from mass and momentum balance to characterize the respective volume fractions and deduce reaction-cross diffusion equations for the space-time evolution of glioma, normal tissue, and necrosis. These are supplemented with a reaction-diffusion equation for the acidity dynamics and lead to formation of patterns which are typical for high grade gliomas. Unlike previous works, our deduction also works in higher dimensions and involves less restrictions. We also investigate the existence of weak solutions to the obtained system of equations and perform numerical simulations to illustrate the solution behavior and the pattern occurrence.

The distribution of hypoxia within tissues plays a critical role in tumor diagnosis and prognosis. Recognizing the significance of tumor oxygenation and hypoxia gradients, we introduce mathematical frameworks grounded in mechanistic modeling approaches for their quantitative assessment within a tumor microenvironment. Our approach provides a non-invasive method to measure and predict hypoxia using known blood vasculature. Formulating a reaction-diffusion model for oxygen distribution, we apply it to derive the corresponding hypoxia profile. The modeling and simulations successfully replicate the observed inter- and intra-tumor heterogeneity in experimentally obtained hypoxia profiles across various tumor tissues (breast, ovarian, and pancreatic) in our dataset. Employing a data-driven approach, we propose a method to deduce partial differential equation (PDE) models with spatially dependent parameters, enabling us to comprehend the variability of hypoxia profiles within a tissue. The versatility of our framework lies not only in capturing diverse and dynamic behaviors of tumor oxygenation but also in categorizing states of vascularization. These categories are distinguished based on the dynamics of oxygen molecules, as identified by the model parameters.