In silico models and computer simulation are invaluable tools to better understand complex biological processes such as cancer evolution. However, the complexity of the biological environment, with many cell mechanisms in response to changing physical and chemical external stimuli, makes the associated mathematical models highly non-linear and multiparametric. One of the main problems of these models is the determination of the parameters’ values, which are usually fitted for specific conditions, making the conclusions drawn difficult to generalise. We analyse here an important biological problem: the evolution of hypoxia-driven migratory structures in Glioblastoma Multiforme (GBM), the most aggressive and lethal primary brain tumour. We establish a mathematical model considering the interaction of the tumour cells with oxygen concentration in what is called the go or grow paradigm. We reproduce in this work three different experiments, showing the main GBM structures (pseudopalisade and necrotic core formation), only changing the initial and boundary conditions. We prove that it is possible to obtain versatile mathematical tools which, together with a sound parametric analysis, allow to explain complex biological phenomena. We show the utility of this hybrid “biomimetic in vitro-in silico” platform to help to elucidate the mechanisms involved in cancer processes, to better understand the role of the different phenomena, to test new scientific hypotheses and to design new data-driven experiments.

No abstract

Cellular adaptation is the ability of cells to change in response to different stimuli and environmental conditions. It occurs via phenotypic plasticity, that is, changes in gene expression derived from changes in the physiological environment. This phenomenon is important in many biological processes, in particular in cancer evolution and its treatment. Therefore, it is crucial to understand the mechanisms behind it. Specifically, the emergence of the cancer stem cell phenotype, showing enhanced proliferation and invasion rates, is an essential process in tumour progression.We present a mathematical framework to simulate phenotypic heterogeneity in different cell populations as a result of their interaction with chemical species in their microenvironment, through a continuum model using the well-known concept of internal variables to model cell phenotype. The resulting model, derived from conservation laws, incorporates the relationship between the phenotype and the history of the stimuli to which cells have been subjected, together with the inheritance of that phenotype. To illustrate the model capabilities, it is particularised for glioblastoma adaptation to hypoxia. A parametric analysis is carried out to investigate the impact of each model parameter regulating cellular adaptation, showing that it permits reproducing different trends reported in the scientific literature. The framework can be easily adapted to any particular problem of cell plasticity, with the main limitation of having enough cells to allow working with continuum variables. With appropriate calibration and validation, it could be useful for exploring the underlying processes of cellular adaptation, as well as for proposing favorable/unfavourable conditions or treatments.

Modeling and simulation are essential tools for better understanding complex biological processes, such as cancer evolution. However, the resulting mathematical models are often highly non-linear and include many parameters, which, in many cases, are difficult to estimate and present strong correlations. Therefore, a proper parametric analysis is mandatory. Following a previous work in which we modeled the in vitro evolution of Glioblastoma Multiforme (GBM) under hypoxic conditions, we analyze and solve here the problem found of parametric correlation. With this aim, we develop a methodology based on copulas to approximate the multidimensional probability density function of the correlated parameters. Once the model is defined, we analyze the experimental setting to optimize the utility of each configuration in terms of gathered information. We prove that experimental configurations with oxygen gradient and high cell concentration have the highest utility when we want to separate correlated effects in our experimental design. We demonstrate that copulas are an adequate tool to analyze highly-correlated multiparametric mathematical models such as those appearing in Biology, with the added value of providing key information for the optimal design of experiments, reducing time and cost in in vivo and in vitro experimental campaigns, like those required in microfluidic models of GBM evolution.