Capturing the Dynamics of a Hybrid Multiscale Cancer Model with a Continuum Model

Cancer is a complex disease involving processes at spatial scales from subcellular, like cell signalling, to tissue scale, such as vascular network formation. A number of multiscale models have been developed to study the dynamics that emerge from the coupling between the intracellular, cellular and tissue scales. Here, we develop a continuum partial differential equation model to capture the dynamics of a particular multiscale model (a hybrid cellular automaton with discrete cells, diffusible factors and an explicit vascular network). The purpose is to test under which circumstances such a continuum model gives equivalent predictions to the original multiscale model, in the knowledge that the system details are known, and differences in model results can be explained in terms of model features (rather than unknown experimental confounding factors). The continuum model qualitatively replicates the dynamics from the multiscale model, with certain discrepancies observed owing to the differences in the modelling of certain processes. The continuum model admits travelling wave solutions for normal tissue growth and tumour invasion, with similar behaviour observed in the multiscale model. However, the continuum model enables us to analyse the spatially homogeneous steady states of the system, and hence to analyse these waves in more detail. We show that the tumour microenvironmental effects from the multiscale model mean that tumour invasion exhibits a so-called pushed wave when the carrying capacity for tumour cell proliferation is less than the total cell density at the tumour wave front. These pushed waves of tumour invasion propagate by triggering apoptosis of normal cells at the wave front. Otherwise, numerical evidence suggests that the wave speed can be predicted from linear analysis about the normal tissue steady state.