Multiphase modelling of glioma pseudopalisading under acidosis

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.

[1]  H M Byrne,et al.  A mathematical model of the stress induced during avascular tumour growth , 2000, Journal of mathematical biology.

[2]  F. Feuerhake,et al.  Why one-size-fits-all vaso-modulatory interventions fail to control glioma invasion: in silico insights , 2016, Scientific Reports.

[3]  C. Engwer,et al.  Effective equations for anisotropic glioma spread with proliferation: a multiscale approach and comparisons with previous settings. , 2016, Mathematical medicine and biology : a journal of the IMA.

[4]  K. Painter,et al.  Mathematical modelling of glioma growth: the use of Diffusion Tensor Imaging (DTI) data to predict the anisotropic pathways of cancer invasion. , 2013, Journal of theoretical biology.

[5]  Daniel J Brat,et al.  Pseudopalisades in Glioblastoma Are Hypoxic, Express Extracellular Matrix Proteases, and Are Formed by an Actively Migrating Cell Population , 2004, Cancer Research.

[7]  Christina Surulescu,et al.  Mathematical modeling of glioma invasion: acid- and vasculature mediated go-or-grow dichotomy and the influence of tissue anisotropy , 2020, Appl. Math. Comput..

[8]  S. Evje,et al.  A multiphase model for exploring tumor cell migration driven by autologous chemotaxis , 2018, Chemical Engineering Science.

[9]  Ansgar Jüngel,et al.  ENTROPY STRUCTURE OF A CROSS-DIFFUSION TUMOR-GROWTH MODEL , 2012 .

[10]  Pawan Kumar,et al.  Multiscale modeling of glioma pseudopalisades: contributions from the tumor microenvironment , 2020, Journal of Mathematical Biology.

[11]  H. Frieboes,et al.  Three-dimensional multispecies nonlinear tumor growth--I Model and numerical method. , 2008, Journal of theoretical biology.

[12]  R. E. Craine,et al.  CONTINUUM THEORIES OF MIXTURES: BASIC THEORY AND HISTORICAL DEVELOPMENT , 1976 .

[13]  Helen M. Byrne,et al.  A two-phase model of solid tumour growth , 2003, Appl. Math. Lett..

[14]  B. Scheithauer,et al.  Histopathology, classification, and grading of gliomas , 1995, Glia.

[15]  R K Jain,et al.  Noninvasive measurement of interstitial pH profiles in normal and neoplastic tissue using fluorescence ratio imaging microscopy. , 1994, Cancer research.

[16]  C. Engwer,et al.  A multiscale model for glioma spread including cell-tissue interactions and proliferation. , 2015, Mathematical biosciences and engineering : MBE.

[17]  Michael Berens,et al.  A mathematical model of glioblastoma tumor spheroid invasion in a three-dimensional in vitro experiment. , 2007, Biophysical journal.

[18]  Matthew P. Jacobson,et al.  Dysregulated pH: a perfect storm for cancer progression , 2011, Nature Reviews Cancer.

[19]  Luigi Preziosi,et al.  Multiphase modeling of tumor growth with matrix remodeling and fibrosis , 2009, Math. Comput. Model..

[20]  Cristian Barbarosie,et al.  REPRESENTATION OF DIVERGENCE-FREE VECTOR FIELDS , 2011 .

[21]  Helen M Byrne,et al.  A mechanical model of tumor encapsulation and transcapsular spread. , 2002, Mathematical biosciences.

[22]  E. Kostelich,et al.  Virtual glioblastoma: growth, migration and treatment in a three‐dimensional mathematical model , 2009, Cell proliferation.

[23]  A. Klar,et al.  Modeling glioma invasion with anisotropy- and hypoxia-triggered motility enhancement: from subcellular dynamics to macroscopic PDEs with multiple taxis , 2020, Mathematical Models and Methods in Applied Sciences.

[24]  Ansgar Jungel,et al.  The boundedness-by-entropy principle for cross-diffusion systems , 2014, 1403.5419.

[25]  Víctor M. Pérez-García,et al.  Hypoxic Cell Waves Around Necrotic Cores in Glioblastoma: A Biomathematical Model and Its Therapeutic Implications , 2012, Bulletin of Mathematical Biology.

[26]  Harald Garcke,et al.  A multiphase Cahn--Hilliard--Darcy model for tumour growth with necrosis , 2017, 1701.06656.

[27]  Luigi Preziosi,et al.  Multiphase modelling of tumour growth and extracellular matrix interaction: mathematical tools and applications , 2009, Journal of mathematical biology.

[28]  J. Lennerz,et al.  Neuropathology for the neuroradiologist: palisades and pseudopalisades. , 2006, AJNR. American journal of neuroradiology.

[29]  D. Drew,et al.  Theory of Multicomponent Fluids , 1998 .

[30]  T. Hillen,et al.  Glioma follow white matter tracts: a multiscale DTI-based model , 2015, Journal of mathematical biology.

[31]  C. Surulescu,et al.  Modeling multiple taxis: Tumor invasion with phenotypic heterogeneity, haptotaxis, and unilateral interspecies repellence , 2020, Discrete & Continuous Dynamical Systems - B.

[32]  Leonardo F. Jurado,et al.  Histopathology , 2019, Fungal Infections of the Central Nervous System.

[33]  Christina Surulescu,et al.  A Flux-Limited Model for Glioma Patterning with Hypoxia-Induced Angiogenesis , 2020, Symmetry.

[34]  Christina Surulescu,et al.  A Multiscale Modeling Approach to Glioma Invasion with Therapy , 2017 .

[35]  J R King,et al.  Multiphase modelling of cell behaviour on artificial scaffolds: effects of nutrient depletion and spatially nonuniform porosity. , 2007, Mathematical medicine and biology : a journal of the IMA.

[36]  H. Byrne,et al.  The role of cell-cell interactions in a two-phase model for avascular tumour growth , 2002, Journal of mathematical biology.

[37]  D. Brat,et al.  Malignant Glioma Physiology: Cellular Response to Hypoxia and Its Role in Tumor Progression , 2003, Annals of Internal Medicine.

[38]  Cass T. Miller,et al.  A multiphase model for three-dimensional tumor growth , 2013, New journal of physics.