Modeling persistence of motion in a crowded environment: The diffusive limit of excluding velocity-jump processes

BIO-LGCA: A cellular automaton modelling class for analysing collective cell migration

Collective dynamics in multicellular systems such as biological organs and tissues plays a key role in biological development, regeneration, and pathological conditions. Collective tissue dynamics—understood as population behaviour arising from the interplay of the constituting discrete cells—can be studied with on- and off-lattice agent-based models. However, classical on-lattice agent-based models, also known as cellular automata, fail to replicate key aspects of collective migration, which is a central instance of collective behaviour in multicellular systems. To overcome drawbacks of classical on-lattice models, we introduce an on-lattice, agent-based modelling class for collective cell migration, which we call biological lattice-gas cellular automaton (BIO-LGCA). The BIO-LGCA is characterised by synchronous time updates, and the explicit consideration of individual cell velocities. While rules in classical cellular automata are typically chosen ad hoc, rules for cell-cell and cell-environment interactions in the BIO-LGCA can also be derived from experimental cell migration data or biophysical laws for individual cell migration. We introduce elementary BIO-LGCA models of fundamental cell interactions, which may be combined in a modular fashion to model complex multicellular phenomena. We exemplify the mathematical mean-field analysis of specific BIO-LGCA models, which allows to explain collective behaviour. The first example predicts the formation of clusters in adhesively interacting cells. The second example is based on a novel BIO-LGCA combining adhesive interactions and alignment. For this model, our analysis clarifies the nature of the recently discovered invasion plasticity of breast cancer cells in heterogeneous environments.

Pattern formation during embryonic development and pathological tissue dynamics, such as cancer invasion, emerge from individual intercellular interactions. In order to study the impact of single cell dynamics and cell-cell interactions on tissue behaviour, one needs to develop space-time-dependent on- or off-lattice agent-based models (ABMs), which consider the behaviour of individual cells. However, classical on-lattice agent-based models also known as cellular automata fail to replicate key aspects of collective migration, which is a central instance of collective behaviour in multicellular systems. Here, we present the rule- and lattice-based BIO-LGCA modelling class which allows for (i) rigorous derivation of rules from biophysical laws and/or experimental data, (ii) mathematical analysis of collective migration, and (iii) computationally efficient simulations.

Modelling collective cell motion: are on- and off-lattice models equivalent?

Biological processes, such as embryonic development, wound repair and cancer invasion, or bacterial swarming and fruiting body formation, involve collective motion of cells as a coordinated group. Collective cell motion of eukaryotic cells often includes interactions that result in polar alignment of cell velocities, while bacterial patterns typically show features of apolar velocity alignment. For analysing the population-scale effects of these different alignment mechanisms, various on- and off-lattice agent-based models have been introduced. However, discriminating model-specific artefacts from general features of collective cell motion is challenging. In this work, we focus on equivalence criteria at the population level to compare on- and off-lattice models. In particular, we define prototypic off- and on-lattice models of polar and apolar alignment, and show how to obtain an on-lattice from an off-lattice model of velocity alignment. By characterizing the behaviour and dynamical description of collective migration models at the macroscopic level, we suggest the type of phase transitions and possible patterns in the approximative macroscopic partial differential equation descriptions as informative equivalence criteria between on- and off-lattice models. This article is part of the theme issue 'Multi-scale analysis and modelling of collective migration in biological systems'.

Persistent exclusion processes: Inertia, drift, mixing, and correlation

In many biological systems, motile agents exhibit random motion with short-term directional persistence, together with crowding effects arising from spatial exclusion. We formulate and study a class of lattice-based models for multiple walkers with motion persistence and spatial exclusion in one and two dimensions, and use a mean-field approximation to investigate relevant population-level partial differential equations in the continuum limit. We show that this model of a persistent exclusion process is in general well described by a nonlinear diffusion equation. With reference to results presented in the current literature, our results reveal that the nonlinearity arises from the combination of motion persistence and volume exclusion, with linearity in terms of the canonical diffusion or heat equation being recovered in either the case of persistence without spatial exclusion, or spatial exclusion without persistence. We generalize our results to include systems of multiple species of interacting, motion-persistent walkers, as well as to incorporate a global drift in addition to persistence. These models are shown to be governed approximately by systems of nonlinear advection-diffusion equations. By comparing the prediction of the mean-field approximation to stochastic simulation results, we assess the performance of our results. Finally, we also address the problem of inferring the presence of persistence from simulation results, with a view to application to experimental cell-imaging data.

The invasion speed of cell migration models with realistic cell cycle time distributions

Cell proliferation is typically incorporated into stochastic mathematical models of cell migration by assuming that cell divisions occur after an exponentially distributed waiting time. Experimental observations, however, show that this assumption is often far from the real cell cycle time distribution (CCTD). Recent studies have suggested an alternative approach to modelling cell proliferation based on a multi-stage representation of the CCTD. In this paper we investigate the connection between the CCTD and the speed of the collective invasion. We first state a result for a general CCTD, which allows the computation of the invasion speed using the Laplace transform of the CCTD. We use this to deduce the range of speeds for the general case. We then focus on the more realistic case of multi-stage models, using both a stochastic agent-based model and a set of reaction-diffusion equations for the cells' average density. By studying the corresponding travelling wave solutions, we obtain an analytical expression for the speed of invasion for a general N-stage model with identical transition rates, in which case the resulting cell cycle times are Erlang distributed. We show that, for a general N-stage model, the Erlang distribution and the exponential distribution lead to the minimum and maximum invasion speed, respectively. This result allows us to determine the range of possible invasion speeds in terms of the average proliferation time for any multi-stage model.