On some models for cancer cell migration through tissue networks

Modelling non-local cell-cell adhesion: a multiscale approach

Cell-cell adhesion plays a vital role in the development and maintenance of multicellular organisms. One of its functions is regulation of cell migration, such as occurs, e.g. during embryogenesis or in cancer. In this work, we develop a versatile multiscale approach to modelling a moving self-adhesive cell population that combines a careful microscopic description of a deterministic adhesion-driven motion component with an efficient mesoscopic representation of a stochastic velocity-jump process. This approach gives rise to mesoscopic models in the form of kinetic transport equations featuring multiple non-localities. Subsequent parabolic and hyperbolic scalings produce general classes of equations with non-local adhesion and myopic diffusion, a special case being the classical macroscopic model proposed in Armstrong et al. (J Theoret Biol 243(1): 98-113, 2006). Our simulations show how the combination of the two motion effects can unfold. Cell-cell adhesion relies on the subcellular cell adhesion molecule binding. Our approach lends itself conveniently to capturing this microscopic effect. On the macroscale, this results in an additional non-linear integral equation of a novel type that is coupled to the cell density equation.

Risk of lung cancer due to external environmental factor and epidemiological data analysis

Lung cancer is a cancer with the fastest growth in the incidence and mortality all over the world, which is an extremely serious threat to human's life and health. Evidences reveal that external environmental factors are the key drivers of lung cancer, such as smoking, radiation exposure and so on. Therefore, it is urgent to explain the mechanism of lung cancer risk due to external environmental factors experimentally and theoretically. However, it is still an open issue regarding how external environment factors affect lung cancer risk. In this paper, we summarize the main mathematical models involved the gene mutations for cancers, and review the application of the models to analyze the mechanism of lung cancer and the risk of lung cancer due to external environmental exposure. In addition, we apply the model described and the epidemiological data to analyze the influence of external environmental factors on lung cancer risk. The result indicates that radiation can cause significantly an increase in the mutation rate of cells, in particular the mutation in stability gene that leads to genomic instability. These studies not only can offer insights into the relationship between external environmental factors and human lung cancer risk, but also can provide theoretical guidance for the prevention and control of lung cancer.

Multiscale modeling of glioma pseudopalisades: contributions from the tumor microenvironment

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 .

Estimating the extent of glioblastoma invasion : Approximate stationalization of anisotropic advection-diffusion-reaction equations in the context of glioblastoma invasion

Glioblastoma Multiforme is a malignant brain tumor with poor prognosis. There have been numerous attempts to model the invasion of tumorous glioma cells via partial differential equations in the form of advection–diffusion–reaction equations. The patient-wise parametrization of these models, and their validation via experimental data has been found to be difficult, as time sequence measurements are mostly missing. Also the clinical interest lies in the actual (invisible) tumor extent for a particular MRI/DTI scan and not in a predictive estimate. Therefore we propose a stationalized approach to estimate the extent of glioblastoma (GBM) invasion at the time of a given MRI/DTI scan. The underlying dynamics can be derived from an instationary GBM model, falling into the wide class of advection-diffusion-reaction equations. The stationalization is introduced via an analytic solution of the Fisher-KPP equation, the simplest model in the considered model class. We investigate the applicability in 1D and 2D, in the presence of inhomogeneous diffusion coefficients and on a real 3D DTI-dataset.

A multiscale model for glioma spread including cell-tissue interactions and proliferation

Glioma is a broad class of brain and spinal cord tumors arising from glia cells, which are the main brain cells that can develop into neoplasms. They are highly invasive and lead to irregular tumor margins which are not precisely identifiable by medical imaging, thus rendering a precise enough resection very difficult. The understanding of glioma spread patterns is hence essential for both radiological therapy as well as surgical treatment. In this paper we propose a multiscale model for glioma growth including interactions of the cells with the underlying tissue network, along with proliferative effects. Our current accounting for two subpopulations of cells to accomodate proliferation according to the go-or-grow dichtomoty is an extension of the setting in [16]. As in that paper, we assume that cancer cells use neuronal fiber tracts as invasive pathways. Hence, the individual structure of brain tissue seems to be decisive for the tumor spread. Diffusion tensor imaging (DTI) is able to provide such information, thus opening the way for patient specific modeling of glioma invasion. Starting from a multiscale model involving subcellular (microscopic) and individual (mesoscale) cell dynamics, we perform a parabolic scaling to obtain an approximating reaction-diffusion-transport equation on the macroscale of the tumor cell population. Numerical simulations based on DTI data are carried out in order to assess the performance of our modeling approach.

Glioma follow white matter tracts: a multiscale DTI-based model

Gliomas are a class of rarely curable tumors arising from abnormal glia cells in the human brain. The understanding of glioma spreading patterns is essential for both radiological therapy as well as surgical treatment. Diffusion tensor imaging (DTI) allows to infer the white matter fibre structure of the brain in a noninvasive way. Painter and Hillen (J Theor Biol 323:25-39, 2013) used a kinetic partial differential equation to include DTI data into a class of anisotropic diffusion models for glioma spread. Here we extend this model to explicitly include adhesion mechanisms between glioma cells and the extracellular matrix components which are associated to white matter tracts. The mathematical modelling follows the multiscale approach proposed by Kelkel and Surulescu (Math Models Methods Appl Sci 23(3), 2012). We use scaling arguments to deduce a macroscopic advection-diffusion model for this process. The tumor diffusion tensor and the tumor drift velocity depend on both, the directions of the white matter tracts as well as the binding dynamics of the adhesion molecules. The advanced computational platform DUNE enables us to accurately solve our macroscopic model. It turns out that the inclusion of cell binding dynamics on the microlevel is an important factor to explain finger-like spread of glioma.

Mathematical Analysis and Numerical Simulations for a System Modeling Acid-Mediated Tumor Cell Invasion

This work is concerned with the mathematical analysis of a model proposed by Gatenby and Gawlinski (1996) in order to support the hypothesis that tumor-induced alteration of microenvironmental pH may provide a simple but comprehensive mechanism to explain cancer invasion. We give an intuitive proof for the existence of a solution under general initial conditions upon using an iterative approach. Numerical simulations are also performed, which endorse the predictions of the model when compared with experimentally observed qualitative facts.