Physica A: Statistical Mechanics and its Applications
A model for mesoscale patterns in motile populations
Introduction
Cell migration is a multiscale phenomenon where experimental observations are made over a range of scales depending on the system of interest. For example, macroscopic observations of cell densities are made to characterize the invasiveness of cell populations [1], [2], whereas microscopic observations of individual cell trajectories are made to characterize the details of the migration mechanism [3], [4]. Observations of cell migration are also made at intermediate scales, which we shall call the mesoscale. These observations correspond to visual patterns that form within a motile population. Two common mesoscale patterns in cell biology include the formation of motile aggregates of cells [5], [6], as well as the formation of chain-like patterns [7], [8], [9].
Traditional modeling approaches were often interested in describing macroscopic properties of a cell population, such as the spatial distribution of cell densities. For such applications continuum models were appropriate, and various motility mechanisms have been considered [1], [2]. Recent advances in microscopy now provide high quality confocal, time-lapse and magnetic resonance imaging data. This new data provides opportunities to develop models at the individual cell level [10].
There is great interest in the correspondence between the discrete and continuum models of cell motility. The starting point can be a spatially discrete, continuous in time, nearest-neighbor master equation for a probability distribution function [11], [12], [13], [14]. In the appropriate limits, this gives rise to a partial differential equation (PDE) description of cell density. An alternative approach starts with a spatially and temporally discrete exclusion process with biologically realistic discrete cell motility rules. The spatially and temporally discrete model can be averaged to provide a continuum description of the system [15], [16], [17]. This is the approach we take here since the spatially and temporally discrete exclusion process is well suited for the interpretation of time-lapse data.
Here we describe and analyze a lattice-based discrete motility model which takes account of contacts or binding between agents. Motility events are governed, in part, by a binding function , where is the scaled coordination number of the target site and is a function that will bias movements to favor a particular coordination number at the target site. The discrete model, based on a simple exclusion process [18], is averaged to give a continuum description of the expected behavior of the system in terms of: (i) a PDE description of the population density [15], [18], and (ii) differential equations describing the average trajectory, or pathline, of a tagged agent within the population [16]. Using simulation and visual inspection, we demonstrate that the model can produce a range of mesoscale patterns including chains and aggregates. Different patterns correspond to different choices of and different initial conditions.
We find that the macroscopic density is governed by a nonlinear diffusion equation. This is an intriguing outcome. Although standard unbiased exclusion processes involve interacting agents, the population-level density obeys a linear diffusion equation, since the interactions are symmetric and do not appear in the continuum model [15], [18]. In the discrete model presented here, the discrete interactions are asymmetric and now appear in the continuum description through the nonlinear diffusivity , where is the continuum density.
The PDE model developed here admits a range of solutions with complex behaviors depending on the form of . In particular, the PDE with appropriate initial conditions admits both smooth and discontinuous solutions, depending on whether is positive or negative [19], [20]. The discrete model can also generate aggregation patterns when the diffusivity is negative. This behavior is related to the backward heat equation [19], [21]. This contrasts the discrete exclusion process of Deroulers et al. that accounted for maintaining neighbors and gave rise to a positive nonlinear diffusivity [17].
Our analysis suggests that the formation of strong mesoscale patterns (chains, aggregates) and the validity of the continuum models described here depend on the form of . Furthermore our results suggest that strong mesoscale patterns may not have any continuum description. This means that the application of a continuum model to capture a mesoscale chain-like structure may be inappropriate. Although this study is motivated with examples from cell biology, this work has implications for related applications involving collective migration that include ecological applications [22], malignant invasion [23] and pedestrian motility [24].
Section snippets
The discrete model
Simulation data presented in this work will use the two-dimensional square lattice as a specific illustration of the model. However, most of the analysis applies to arbitrary -dimensional periodic lattices with spacing . Agents can be viewed as occupying sites. Alternatively, agents can be regarded as residing in regions, since each site is associated with a spatial region consisting of all points closer to site than to any other. In various contexts these regions are known as Voronoi
Discrete simulation data: Influence of
To demonstrate the influence of different binding functions, we present a suite of simulations on a two-dimensional square lattice of size 100×100. Periodic boundary conditions are imposed along all boundaries. Initially, each site with is occupied with probability 0.5. Snapshots of four simulations, after 200 time steps, are shown in Fig. 2. Results in Fig. 2(a), with , show a symmetric distribution of agents where there is no apparent preference for the local
Deriving a continuous description
To connect the discrete mechanism with a continuum model we define the average occupancy of site , averaged over many statistically identical realizations, as . After averaging, we form a discrete conservation statement describing , which is the change in average occupancy of site during the time interval from to [15], [18]: The positive terms on the right of Eq. (1) represent the
Results
Simulation data will be generated and compared to the solution of the continuum models for two types of problems. First, we consider the evolution of an initially close-packed group of agents. Second, we consider the evolution of a uniformly seeded lattice for which, on average, there are no spatial gradients present initially in the system. These two problems demonstrate a wide range of behaviors supported by the discrete model and demonstrate the relationship between the discrete mechanism
Discussion and conclusion
We have presented and analyzed a discrete motility model capable of producing a range of mesoscale patterns. Our approach enables us to visualize discrete simulations, as well as relating our discrete model to continuum equations describing agent density and agent pathline data. Previous approaches to modeling mesoscale patterns have either taken a discrete or continuum approach separately [13], [37], [38], without any explicit connection between the two. The development and application of
Acknowledgements
This work is supported by the Australian Research Council (ARC). Mat Simpson is an ARC Postdoctoral Fellow and Kerry Landman is an ARC Professorial Fellow. We thank Dr. Heather Young and Dr. Don Newgreen whose experimental work was the motivation for this paper.
References (43)
- et al.
Appl. Math. Lett.
(2004) - et al.
Dev. Biol.
(2007) - et al.
Chem. Eng. Sci.
(1989) - et al.
Dev. Biol.
(2004) - et al.
J. Theoret. Biol.
(2003) - et al.
Physica A
(2009) Appl. Math. Lett.
(1995)Physica A
(2006)- et al.
Phys. Life Rev.
(2005) Physica A
(2002)
J. Theoret. Biol.
J. Theoret. Biol.
J. Neurosci.
Neoplasia
Biotechnol. Bioeng.
Development
Cancer Res.
Phys. Rev. E
SIAM. J. Appl. Math.
Bull. Math. Biol.
J. Math. Biol.
Cited by (48)
Nonlinear stability of shock-fronted travelling waves in reaction-nonlinear diffusion equations
2024, Physica D: Nonlinear PhenomenaReliable and efficient parameter estimation using approximate continuum limit descriptions of stochastic models
2022, Journal of Theoretical BiologyCitation Excerpt :Modelling this kind of experiment is straightforward in a stochastic framework but leads to a more computationally demanding simulation algorithm with more parameters to be estimated (Zhang et al., 2010). Another extension would be to consider modelling how agents interact with each other, for example through modelling cell-to-cell adhesion (Deroulers et al., 2009) or different types of cell-to-cell attraction/repulsion mechanisms (Simpson et al., 2010a). This approach is relevant for modelling various types of experiments involving adhesive cancer cells, and incorporating adhesion into the mathematical model is of interest because describing such mechanisms introduces additional parameters to be estimated.
Shock-fronted travelling waves in a reaction–diffusion model with nonlinear forward–backward–forward diffusion
2021, Physica D: Nonlinear PhenomenaHole-closing model reveals exponents for nonlinear degenerate diffusivity functions in cell biology
2019, Physica D: Nonlinear PhenomenaQuantifying the effect of experimental design choices for in vitro scratch assays
2016, Journal of Theoretical BiologyWhen push comes to shove: Exclusion processes with nonlocal consequences
2015, Physica A: Statistical Mechanics and its ApplicationsCitation Excerpt :Here we study the implications of local movement rules with shoving for the macroscopic description of the population of agents, which may be a single species or multiple species of agents. We shall investigate whether the continuum limit of these probabilistic rules gives rise to a nonlinear diffusion equation and if it does, whether the diffusion coefficient is a polynomial in the average occupancy, as in previous work [9,11,14,15]. Alternatively, since our local movement rules have non-local consequences, might one obtain a different form, for example a higher order diffusion equation [27]?