Modelling cell adaptation using internal variables: Accounting for cell plasticity in continuum mathematical biology.

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 favourable/unfavourable conditions or treatments.