Efficient coarse simulation of a growing avascular tumor

Data-driven spatio-temporal modelling of glioblastoma

Mathematical oncology provides unique and invaluable insights into tumour growth on both the microscopic and macroscopic levels. This review presents state-of-the-art modelling techniques and focuses on their role in understanding glioblastoma, a malignant form of brain cancer. For each approach, we summarize the scope, drawbacks and assets. We highlight the potential clinical applications of each modelling technique and discuss the connections between the mathematical models and the molecular and imaging data used to inform them. By doing so, we aim to prime cancer researchers with current and emerging computational tools for understanding tumour progression. By providing an in-depth picture of the different modelling techniques, we also aim to assist researchers who seek to build and develop their own models and the associated inference frameworks. Our article thus strikes a unique balance. On the one hand, we provide a comprehensive overview of the available modelling techniques and their applications, including key mathematical expressions. On the other hand, the content is accessible to mathematicians and biomedical scientists alike to accommodate the interdisciplinary nature of cancer research.

Modeling the dynamics of oligodendrocyte precursor cells and the genesis of gliomas

Oligodendrocyte precursor cells (OPCs) have remarkable properties: they represent the most abundant cycling cell population in the adult normal brain and they manage to achieve a uniform and constant density throughout the adult brain. This equilibrium is obtained by the interplay of four processes: division, differentiation or death, migration and active self-repulsion. They are also strongly suspected to be at the origin of gliomas, when their equilibrium is disrupted. In this article, we present a model of the dynamics of OPCs, first in a normal tissue. This model is based on a cellular automaton and its rules are mimicking the ones that regulate the dynamics of real OPCs. The model is able to reproduce the homeostasis of the cell population, with the maintenance of a constant and uniform cell density and the healing of a lesion. We show that there exists a fair quantitative agreement between the simulated and experimental parameters, such as the cell velocity, the time taken to close a lesion, and the duration of the cell cycle. We present three possible scenarios of disruption of the equilibrium: the appearance of an over-proliferating cell, of a deadless/non-differentiating cell, or of a cell that lost any contact-inhibition. We show that the appearance of an over-proliferating cell is sufficient to trigger the growth of a tumor that has low-grade glioma features: an invasive behaviour, a linear radial growth of the tumor with a corresponding growth velocity of less than 2 mm per year, as well a cell density at the center which exceeds the one in normal tissue by a factor of less than two. The loss of contact inhibition leads to a more high-grade-like glioma. The results of our model contribute to the body of evidence that identify OPCs as possible cells of origin of gliomas.

Gliomas are the most common brain tumors and result in more years of life lost than any other tumor. Standard treatments only confer a limited improvement in overall survival, underscoring the need for new therapies. Finding the type of cells at the origin of these tumors could lead to the development of new drugs, specifically targeted towards these cells. The oligodendrocyte precursor cells are suspected to be these cells of origin, because they continue to proliferate through all the adult life. In this article, we present a model of the dynamics of these cells, first in the normal brain, and then we extrapolate our model to the pathological situation. We study several scenarios where, from the normal situation, a cell appears with one property different from those of the normal cells. We show that the alteration of only one of the properties of these cells in the model can lead to the formation of gliomas with different aggressiveness and very similar to real gliomas, reinforcing the suspicion that the precursor cells are at the origin of gliomas.

Simulation of avascular tumor growth by agent-based game model involving phenotype-phenotype interactions

All tumors, both benign and metastatic, undergo an avascular growth stage with nutrients supplied by the surrounding tissue. This avascular growth process is much easier to carry out in more qualitative and quantitative experiments starting from tumor spheroids in vitro with reliable reproducibility. Essentially, this tumor progression would be described as a sequence of phenotypes. Using agent-based simulation in a two-dimensional spatial lattice, we constructed a composite growth model in which the phenotypic behavior of tumor cells depends on not only the local nutrient concentration and cell count but also the game among cells. Our simulation results demonstrated that in silico tumors are qualitatively similar to those observed in tumor spheroid experiments. We also found that the payoffs in the game between two living cell phenotypes can influence the growth velocity and surface roughness of tumors at the same time. Finally, this current model is flexible and can be easily extended to discuss other situations, such as environmental heterogeneity and mutation.