Modelling Directional Guidance and Motility Regulation in Cell Migration

In vivo imaging in transgenic songbirds reveals superdiffusive neuron migration in the adult brain

Neuron migration is a key phase of neurogenesis, critical for the assembly and function of neuronal circuits. In songbirds, this process continues throughout life, but how these newborn neurons disperse through the adult brain is unclear. We address this question using in vivo two-photon imaging in transgenic zebra finches that express GFP in young neurons and other cell types. In juvenile and adult birds, migratory cells are present at a high density, travel in all directions, and make frequent course changes. Notably, these dynamic migration patterns are well fit by a superdiffusive model. Simulations reveal that these superdiffusive dynamics are sufficient to disperse new neurons throughout the song nucleus HVC. These results suggest that superdiffusive migration may underlie the formation and maintenance of nuclear brain structures in the postnatal brain and indicate that transgenic songbirds are a useful resource for future studies into the mechanisms of adult neurogenesis.

Learning anisotropic interaction rules from individual trajectories in a heterogeneous cellular population

Interacting particle system (IPS) models have proven to be highly successful for describing the spatial movement of organisms. However, it is challenging to infer the interaction rules directly from data. In the field of equation discovery, the weak-form sparse identification of nonlinear dynamics (WSINDy) methodology has been shown to be computationally efficient for identifying the governing equations of complex systems from noisy data. Motivated by the success of IPS models to describe the spatial movement of organisms, we develop WSINDy for the second-order IPS to learn equations for communities of cells. Our approach learns the directional interaction rules for each individual cell that in aggregate govern the dynamics of a heterogeneous population of migrating cells. To sort a cell according to the active classes present in its model, we also develop a novel ad hoc classification scheme (which accounts for the fact that some cells do not have enough evidence to accurately infer a model). Aggregated models are then constructed hierarchically to simultaneously identify different species of cells present in the population and determine best-fit models for each species. We demonstrate the efficiency and proficiency of the method on several test scenarios, motivated by common cell migration experiments.

Population dynamics with spatial structure and an Allee effect

Population dynamics including a strong Allee effect describe the situation where long-term population survival or extinction depends on the initial population density. A simple mathematical model of an Allee effect is one where initial densities below the threshold lead to extinction, whereas initial densities above the threshold lead to survival. Mean-field models of population dynamics neglect spatial structure that can arise through short-range interactions, such as competition and dispersal. The influence of non-mean-field effects has not been studied in the presence of an Allee effect. To address this, we develop an individual-based model that incorporates both short-range interactions and an Allee effect. To explore the role of spatial structure we derive a mathematically tractable continuum approximation of the IBM in terms of the dynamics of spatial moments. In the limit of long-range interactions where the mean-field approximation holds, our modelling framework recovers the mean-field Allee threshold. We show that the Allee threshold is sensitive to spatial structure neglected by mean-field models. For example, there are cases where the mean-field model predicts extinction but the population actually survives. Through simulations we show that our new spatial moment dynamics model accurately captures the modified Allee threshold in the presence of spatial structure.

Identifying density-dependent interactions in collective cell behaviour

Scratch assays are routinely used to study collective cell behaviour in vitro. Typical experimental protocols do not vary the initial density of cells, and typical mathematical modelling approaches describe cell motility and proliferation based on assumptions of linear diffusion and logistic growth. Jin et al. (Jin et al. 2016 J. Theor. Biol. 390, 136-145 (doi:10.1016/j.jtbi.2015.10.040)) find that the behaviour of cells in scratch assays is density-dependent, and show that standard modelling approaches cannot simultaneously describe data initiated across a range of initial densities. To address this limitation, we calibrate an individual-based model to scratch assay data across a large range of initial densities. Our model allows proliferation, motility, and a direction bias to depend on interactions between neighbouring cells. By considering a hierarchy of models where we systematically and sequentially remove interactions, we perform model selection analysis to identify the minimum interactions required for the model to simultaneously describe data across all initial densities. The calibrated model is able to match the experimental data across all densities using a single parameter distribution, and captures details about the spatial structure of cells. Our results provide strong evidence to suggest that motility is density-dependent in these experiments. On the other hand, we do not see the effect of crowding on proliferation in these experiments. These results are significant as they are precisely the opposite of the assumptions in standard continuum models, such as the Fisher-Kolmogorov equation and its generalizations.

A computational model of bidirectional axonal growth in micro-tissue engineered neuronal networks (micro-TENNs)

Micro-Tissue Engineered Neural Networks (Micro-TENNs) are living three-dimensional constructs designed to replicate the neuroanatomy of white matter pathways in the brain and are being developed as implantable micro-tissue for axon tract reconstruction, or as anatomically-relevant in vitro experimental platforms. Micro-TENNs are composed of discrete neuronal aggregates connected by bundles of long-projecting axonal tracts within miniature tubular hydrogels. In order to help design and optimize micro-TENN performance, we have created a new computational model including geometric and functional properties. The model is built upon the three-dimensional diffusion equation and incorporates large-scale uni- and bi-directional growth that simulates realistic neuron morphologies. The model captures unique features of 3D axonal tract development that are not apparent in planar outgrowth and may be insightful for how white matter pathways form during brain development. The processes of axonal outgrowth, branching, turning and aggregation/bundling from each neuron are described through functions built on concentration equations and growth time distributed across the growth segments. Once developed we conducted multiple parametric studies to explore the applicability of the method and conducted preliminary validation via comparisons to experimentally grown micro-TENNs for a range of growth conditions. Using this framework, the model can be applied to study micro-TENN growth processes and functional characteristics using spiking network or compartmental network modeling. This model may be applied to improve our understanding of axonal tract development and functionality, as well as to optimize the fabrication of implantable tissue engineered brain pathways for nervous system reconstruction and/or modulation.

Rhodoptilometrin, a Crinoid-Derived Anthraquinone, Induces Cell Regeneration by Promoting Wound Healing and Oxidative Phosphorylation in Human Gingival Fibroblast Cells

Gingival recession (GR) potentially leads to the exposure of tooth root to the oral cavity microenvironment and increases susceptibility to dental caries, dentin hypersensitivity, and other dental diseases. Even though many etiological factors were reported, the specific mechanism of GR is yet to be elucidated. Given the species richness concerning marine biodiversity, it could be a treasure trove for drug discovery. In this study, we demonstrate the effects of a marine compound, (+)-rhodoptilometrin from crinoid, on gingival cell migration, wound healing, and oxidative phosphorylation (OXPHOS). Experimental results showed that (+)-rhodoptilometrin can significantly increase wound healing, migration, and proliferation of human gingival fibroblast cells, and it does not have effects on oral mucosa fibroblast cells. In addition, (+)-rhodoptilometrin increases the gene and protein expression levels of focal adhesion kinase (FAK), fibronectin, and type I collagen, changes the intracellular distribution of FAK and F-actin, and increases OXPHOS and the expression levels of complexes I~V in the mitochondria. Based on our results, we believe that (+)-rhodoptilometrin might increase FAK expression and promote mitochondrial function to affect cell migration and promote gingival regeneration. Therefore, (+)-rhodoptilometrin may be a promising therapeutic agent for GR.

Pattern formation in multiphase models of chemotactic cell aggregation

We develop a continuum model for the aggregation of cells cultured in a nutrient-rich medium in a culture well. We consider a 2D geometry, representing a vertical slice through the culture well, and assume that the cell layer depth is small compared with the typical lengthscale of the culture well. We adopt a continuum mechanics approach, treating the cells and culture medium as a two-phase mixture. Specifically, the cells and culture medium are treated as fluids. Additionally, the cell phase can generate forces in response to environmental cues, which include the concentration of a chemoattractant that is produced by the cells within the culture medium. The model leads to a system of coupled nonlinear partial differential equations for the volume fraction and velocity of the cell phase, the culture medium pressure and the chemoattractant concentration, which must be solved subject to appropriate boundary and initial conditions. To gain insight into the system, we consider two model reductions, appropriate when the cell layer depth is thin compared to the typical length scale of the culture well: a (simple) 1D and a (more involved) thin-film extensional flow reduction. By investigating the resulting systems of equations analytically and numerically, we identify conditions under which small amplitude perturbations to a homogeneous steady state (corresponding to a spatially uniform cell distribution) can lead to a spatially varying steady state (pattern formation). Our analysis reveals that the simpler 1D reduction has the same qualitative features as the thin-film extensional flow reduction in the linear and weakly nonlinear regimes, motivating the use of the simpler 1D modelling approach when a qualitative understanding of the system is required. However, the thin-film extensional flow reduction may be more appropriate when detailed quantitative agreement between modelling predictions and experimental data is desired. Furthermore, full numerical simulations of the two model reductions in regions of parameter space when the system is not close to marginal stability reveal significant differences in the evolution of the volume fraction and velocity of the cell phase, and chemoattractant concentration.

Scalable population-level modelling of biological cells incorporating mechanics and kinetics in continuous time

The processes taking place inside the living cell are now understood to the point where predictive computational models can be used to gain detailed understanding of important biological phenomena. A key challenge is to extrapolate this detailed knowledge of the individual cell to be able to explain at the population level how cells interact and respond with each other and their environment. In particular, the goal is to understand how organisms develop, maintain and repair functional tissues and organs. In this paper, we propose a novel computational framework for modelling populations of interacting cells. Our framework incorporates mechanistic, constitutive descriptions of biomechanical properties of the cell population, and uses a coarse-graining approach to derive individual rate laws that enable propagation of the population through time. Thanks to its multiscale nature, the resulting simulation algorithm is extremely scalable and highly efficient. As highlighted in our computational examples, the framework is also very flexible and may straightforwardly be coupled with continuous-time descriptions of biochemical signalling within, and between, individual cells.

Distinguishing cell shoving mechanisms

Motivated by in vitro time–lapse images of ovarian cancer spheroids inducing mesothelial cell clearance, the traditional agent–based model of cell migration, based on simple volume exclusion, was extended to include the possibility that a cell seeking to move into an occupied location may push the resident cell, and any cells neighbouring it, out of the way to occupy that location. In traditional discrete models of motile cells with volume exclusion such a move would be aborted. We introduce a new shoving mechanism which allows cells to choose the direction to shove cells that expends the least amount of shoving effort (to account for the likely resistance of cells to being pushed). We call this motility rule ‘smart shoving’. We examine whether agent–based simulations of different shoving mechanisms can be distinguished on the basis of single realisations and averages over many realisations. We emphasise the difficulty in distinguishing cell mechanisms from cellular automata simulations based on snap–shots of cell distributions, site–occupancy averages and the evolution of the number of cells of each species averaged over many realisations. This difficulty suggests the need for higher resolution cell tracking.

Quantifying the roles of random motility and directed motility using advection-diffusion theory for a 3T3 fibroblast cell migration assay stimulated with an electric field

Background

Directed cell migration can be driven by a range of external stimuli, such as spatial gradients of: chemical signals (chemotaxis); adhesion sites (haptotaxis); or temperature (thermotaxis). Continuum models of cell migration typically include a diffusion term to capture the undirected component of cell motility and an advection term to capture the directed component of cell motility. However, there is no consensus in the literature about the form that the advection term takes. Some theoretical studies suggest that the advection term ought to include receptor saturation effects. However, others adopt a much simpler constant coefficient. One of the limitations of including receptor saturation effects is that it introduces several additional unknown parameters into the model. Therefore, a relevant research question is to investigate whether directed cell migration is best described by a simple constant tactic coefficient or a more complicated model incorporating saturation effects.

Results

We study directed cell migration using an experimental device in which the directed component of the cell motility is driven by a spatial gradient of electric potential, which is known as electrotaxis. The electric field (EF) is proportional to the spatial gradient of the electric potential. The spatial variation of electric potential across the experimental device varies in such a way that there are several subregions on the device in which the EF takes on different values that are approximately constant within those subregions. We use cell trajectory data to quantify the motion of 3T3 fibroblast cells at different locations on the device to examine how different values of the EF influences cell motility. The undirected (random) motility of the cells is quantified in terms of the cell diffusivity, D, and the directed motility is quantified in terms of a cell drift velocity, v. Estimates D and v are obtained under a range of four different EF conditions, which correspond to normal physiological conditions. Our results suggest that there is no anisotropy in D, and that D appears to be approximately independent of the EF and the electric potential. The drift velocity increases approximately linearly with the EF, suggesting that the simplest linear advection term, with no additional saturation parameters, provides a good explanation of these physiologically relevant data.

Conclusions

We find that the simplest linear advection term in a continuum model of directed cell motility is sufficient to describe a range of different electrotaxis experiments for 3T3 fibroblast cells subject to normal physiological values of the electric field. This is useful information because alternative models that include saturation effects involve additional parameters that need to be estimated before a partial differential equation model can be applied to interpret or predict a cell migration experiment.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0413-5) contains supplementary material, which is available to authorized users.

Phase separation in driven granular gases: exploring the elusive character of nonequilibrium steady states

The emergence of patterns and phase separation in many-body systems far from thermal equilibrium is discussed using the example of driven granular gases. It is shown that phase separation follows a similar mechanism as in the systems of active Brownian particles. Depending on the quantities chosen for observation, it may or may not be easy to find functionals analogous to the free energy in equilibrium statistical physics. We argue that although such functionals can always be derived from the dynamics, it is of only limited value for predicting relevant aspects of the nonequilibrium steady state of the system. Consequently, although there is indeed a 'principle' governing the selection of collective nonequilibrium steady states (and the corresponding large deviation functional can be identified), it is not generally useful for predicting the behaviour of the system.

Spatial moment dynamics for collective cell movement incorporating a neighbour-dependent directional bias

The ability of cells to undergo collective movement plays a fundamental role in tissue repair, development and cancer. Interactions occurring at the level of individual cells may lead to the development of spatial structure which will affect the dynamics of migrating cells at a population level. Models that try to predict population-level behaviour often take a mean-field approach, which assumes that individuals interact with one another in proportion to their average density and ignores the presence of any small-scale spatial structure. In this work, we develop a lattice-free individual-based model (IBM) that uses random walk theory to model the stochastic interactions occurring at the scale of individual migrating cells. We incorporate a mechanism for local directional bias such that an individual's direction of movement is dependent on the degree of cell crowding in its neighbourhood. As an alternative to the mean-field approach, we also employ spatial moment theory to develop a population-level model which accounts for spatial structure and predicts how these individual-level interactions propagate to the scale of the whole population. The IBM is used to derive an equation for dynamics of the second spatial moment (the average density of pairs of cells) which incorporates the neighbour-dependent directional bias, and we solve this numerically for a spatially homogeneous case.

Spatial structure arising from neighbour-dependent bias in collective cell movement

Mathematical models of collective cell movement often neglect the effects of spatial structure, such as clustering, on the population dynamics. Typically, they assume that individuals interact with one another in proportion to their average density (the mean-field assumption) which means that cell-cell interactions occurring over short spatial ranges are not accounted for. However, in vitro cell culture studies have shown that spatial correlations can play an important role in determining collective behaviour. Here, we take a combined experimental and modelling approach to explore how individual-level interactions give rise to spatial structure in a moving cell population. Using imaging data from in vitro experiments, we quantify the extent of spatial structure in a population of 3T3 fibroblast cells. To understand how this spatial structure arises, we develop a lattice-free individual-based model (IBM) and simulate cell movement in two spatial dimensions. Our model allows an individual's direction of movement to be affected by interactions with other cells in its neighbourhood, providing insights into how directional bias generates spatial structure. We consider how this behaviour scales up to the population level by using the IBM to derive a continuum description in terms of the dynamics of spatial moments. In particular, we account for spatial correlations between cells by considering dynamics of the second spatial moment (the average density of pairs of cells). Our numerical results suggest that the moment dynamics description can provide a good approximation to averaged simulation results from the underlying IBM. Using our in vitro data, we estimate parameters for the model and show that it can generate similar spatial structure to that observed in a 3T3 fibroblast cell population.

Random walks in nonuniform environments with local dynamic interactions

We consider a class of lattice random walk models in which the random walker is initially confined to a finite connected set of allowed sites but has the opportunity to enlarge this set by colliding with its boundaries, each such collision having a given probability of breaking through. The model is motivated by an analogy to cell motility in tissue, where motile cells have the ability to remodel extracellular matrix, but is presented here as a generic model for stochastic erosion. For the one-dimensional case, we report some exact analytic results, some mean-field type analytic approximate results and simulations. We compute exactly the mean and variance of the time taken to enlarge the interval from a single site to a given size. The problem of determining the statistics of the interval length and the walker's position at a given time is more difficult and we report several interesting observations from simulations. Our simulations include the case in which the initial interval length is random and the case in which the initial state of the lattice is a random mixture of allowed and forbidden sites, with the walker placed at random on an allowed site. To illustrate the extension of these ideas to higher-dimensional systems, we consider the erosion of the simple cubic lattice commencing from a single site and report simulations of measures of cluster size and shape and the mean-square displacement of the walker.

Multi-dimensional, mesoscopic Monte Carlo simulations of inhomogeneous reaction-drift-diffusion systems on graphics-processing units

For many biological applications, a macroscopic (deterministic) treatment of reaction-drift-diffusion systems is insufficient. Instead, one has to properly handle the stochastic nature of the problem and generate true sample paths of the underlying probability distribution. Unfortunately, stochastic algorithms are computationally expensive and, in most cases, the large number of participating particles renders the relevant parameter regimes inaccessible. In an attempt to address this problem we present a genuine stochastic, multi-dimensional algorithm that solves the inhomogeneous, non-linear, drift-diffusion problem on a mesoscopic level. Our method improves on existing implementations in being multi-dimensional and handling inhomogeneous drift and diffusion. The algorithm is well suited for an implementation on data-parallel hardware architectures such as general-purpose graphics processing units (GPUs). We integrate the method into an operator-splitting approach that decouples chemical reactions from the spatial evolution. We demonstrate the validity and applicability of our algorithm with a comprehensive suite of standard test problems that also serve to quantify the numerical accuracy of the method. We provide a freely available, fully functional GPU implementation. Integration into Inchman, a user-friendly web service, that allows researchers to perform parallel simulations of reaction-drift-diffusion systems on GPU clusters is underway.

Building macroscale models from microscale probabilistic models: a general probabilistic approach for nonlinear diffusion and multispecies phenomena

A discrete agent-based model on a periodic lattice of arbitrary dimension is considered. Agents move to nearest-neighbor sites by a motility mechanism accounting for general interactions, which may include volume exclusion. The partial differential equation describing the average occupancy of the agent population is derived systematically. A diffusion equation arises for all types of interactions and is nonlinear except for the simplest interactions. In addition, multiple species of interacting subpopulations give rise to an advection-diffusion equation for each subpopulation. This work extends and generalizes previous specific results, providing a construction method for determining the transport coefficients in terms of a single conditional transition probability, which depends on the occupancy of sites in an influence region. These coefficients characterize the diffusion of agents in a crowded environment in biological and physical processes.

A random cell motility gradient downstream of FGF controls elongation of an amniote embryo

Vertebrate embryos are characterized by an elongated antero-posterior (AP) body axis, which forms by progressive cell deposition from a posterior growth zone in the embryo. Here, we used tissue ablation in the chicken embryo to demonstrate that the caudal presomitic mesoderm (PSM) has a key role in axis elongation. Using time-lapse microscopy, we analysed the movements of fluorescently labelled cells in the PSM during embryo elongation, which revealed a clear posterior-to-anterior gradient of cell motility and directionality in the PSM. We tracked the movement of the PSM extracellular matrix in parallel with the labelled cells and subtracted the extracellular matrix movement from the global motion of cells. After subtraction, cell motility remained graded but lacked directionality, indicating that the posterior cell movements associated with axis elongation in the PSM are not intrinsic but reflect tissue deformation. The gradient of cell motion along the PSM parallels the fibroblast growth factor (FGF)/mitogen-activated protein kinase (MAPK) gradient, which has been implicated in the control of cell motility in this tissue. Both FGF signalling gain- and loss-of-function experiments lead to disruption of the motility gradient and a slowing down of axis elongation. Furthermore, embryos treated with cell movement inhibitors (blebbistatin or RhoK inhibitor), but not cell cycle inhibitors, show a slower axis elongation rate. We propose that the gradient of random cell motility downstream of FGF signalling in the PSM controls posterior elongation in the amniote embryo. Our data indicate that tissue elongation is an emergent property that arises from the collective regulation of graded, random cell motion rather than by the regulation of directionality of individual cellular movements.

Characterisation of human bone marrow stromal cell heterogeneity for skeletal regeneration strategies using a two-stage colony assay and computational modelling

Skeletal regeneration and tissue engineering strategies rely critically on the efficient expansion of progenitor cell populations whilst simultaneously preserving multipotentiality and the ability to induce differentiation towards bone and cartilage. Cell population heterogeneity has a significant impact on this process, but is currently poorly quantified, hampering the interpretation of experimental results and the design of optimised expansion protocols. The objective of this study was to characterise individual human bone marrow stromal cell heterogeneity in terms of colony expansion potential. For this purpose, a novel two-stage CFU-F assay was developed in which cells from primary single cell-derived colonies were detached and reseeded again at clonal density as single cells to form new secondary colonies. This clearly demonstrated how secondary colony growth potential varies markedly both between and within primary colonies. Depending on the primary colony, cells either generated small secondary colonies only, or else a wide range of colony sizes. Using computational modelling it was shown how such colony heterogeneity could arise from hierarchical progenitor cell populations and what the limits of such a population structure were in explaining the experimental data. In addition the model demonstrated the significant potential impact of cell mobility on expansion potential and its implications for inducing population heterogeneity. This combined experimental-computational approach will ascertain the impact of cell culture protocols on the expansion potential and functional composition of heterogeneous progenitor populations. Such insights are likely to be of crucial importance for the success of skeletal regeneration strategies.

A framework for modeling the growth and development of neurons and networks

The development of neural tissue is a complex organizing process, in which it is difficult to grasp how the various localized interactions between dividing cells leads relentlessly to global network organization. Simulation is a useful tool for exploring such complex processes because it permits rigorous analysis of observed global behavior in terms of the mechanistic axioms declared in the simulated model. We describe a novel simulation tool, CX3D, for modeling the development of large realistic neural networks such as the neocortex, in a physical 3D space. In CX3D, as in biology, neurons arise by the replication and migration of precursors, which mature into cells able to extend axons and dendrites. Individual neurons are discretized into spherical (for the soma) and cylindrical (for neurites) elements that have appropriate mechanical properties. The growth functions of each neuron are encapsulated in set of pre-defined modules that are automatically distributed across its segments during growth. The extracellular space is also discretized, and allows for the diffusion of extracellular signaling molecules, as well as the physical interactions of the many developing neurons. We demonstrate the utility of CX3D by simulating three interesting developmental processes: neocortical lamination based on mechanical properties of tissues; a growth model of a neocortical pyramidal cell based on layer-specific guidance cues; and the formation of a neural network in vitro by employing neurite fasciculation. We also provide some examples in which previous models from the literature are re-implemented in CX3D. Our results suggest that CX3D is a powerful tool for understanding neural development.

Modeling tumor cell migration: From microscopic to macroscopic models

It has been shown experimentally that contact interactions may influence the migration of cancer cells. Previous works have modelized this thanks to stochastic, discrete models (cellular automata) at the cell level. However, for the study of the growth of real-size tumors with several million cells, it is best to use a macroscopic model having the form of a partial differential equation (PDE) for the density of cells. The difficulty is to predict the effect, at the macroscopic scale, of contact interactions that take place at the microscopic scale. To address this, we use a multiscale approach: starting from a very simple, yet experimentally validated, microscopic model of migration with contact interactions, we derive a macroscopic model. We show that a diffusion equation arises, as is often postulated in the field of glioma modeling, but it is nonlinear because of the interactions. We give the explicit dependence of diffusivity on the cell density and on a parameter governing cell-cell interactions. We discuss in detail the conditions of validity of the approximations used in the derivation, and we compare analytic results from our PDE to numerical simulations and to some in vitro experiments. We notice that the family of microscopic models we started from includes as special cases some kinetically constrained models that were introduced for the study of the physics of glasses, supercooled liquids, and jamming systems.

Random walk models in biology

Mathematical modelling of the movement of animals, micro-organisms and cells is of great relevance in the fields of biology, ecology and medicine. Movement models can take many different forms, but the most widely used are based on the extensions of simple random walk processes. In this review paper, our aim is twofold: to introduce the mathematics behind random walks in a straightforward manner and to explain how such models can be used to aid our understanding of biological processes. We introduce the mathematical theory behind the simple random walk and explain how this relates to Brownian motion and diffusive processes in general. We demonstrate how these simple models can be extended to include drift and waiting times or be used to calculate first passage times. We discuss biased random walks and show how hyperbolic models can be used to generate correlated random walks. We cover two main applications of the random walk model. Firstly, we review models and results relating to the movement, dispersal and population redistribution of animals and micro-organisms. This includes direct calculation of mean squared displacement, mean dispersal distance, tortuosity measures, as well as possible limitations of these model approaches. Secondly, oriented movement and chemotaxis models are reviewed. General hyperbolic models based on the linear transport equation are introduced and we show how a reinforced random walk can be used to model movement where the individual changes its environment. We discuss the applications of these models in the context of cell migration leading to blood vessel growth (angiogenesis). Finally, we discuss how the various random walk models and approaches are related and the connections that underpin many of the key processes involved.

Experimental characterization and computational modelling of two-dimensional cell spreading for skeletal regeneration

Limited cell ingrowth is a major problem for tissue engineering and the clinical application of porous biomaterials as bone substitutes. As a first step, migration and proliferation of an interacting cell population can be studied in two-dimensional culture. Mathematical modelling is essential to generalize the results of these experiments and to derive the intrinsic parameters that can be used for predictions. However, a more thorough evaluation of theoretical models is hampered by limited experimental observations. In this study, experiments and image analysis methods were developed to provide a detailed spatial and temporal picture of how cell distributions evolve. These methods were used to quantify the migration and proliferation of skeletal cell types including MG63 and human bone marrow stromal cells (HBMSCs). The high level of detail with which the cell distributions were mapped enabled a precise assessment of the correspondence between experimental results and theoretical model predictions. This analysis revealed that the standard Fisher equation is appropriate for describing the migration behaviour of the HBMSC population, while for the MG63 cells a sharp front model is more appropriate. In combination with experiments, this type of mathematical model will prove useful in predicting cell ingrowth and improving strategies and control of skeletal tissue regeneration.

Simulating invasion with cellular automata: connecting cell-scale and population-scale properties

Interpretive and predictive tools are needed to assist in the understanding of cell invasion processes. Cell invasion involves cell motility and proliferation, and is central to many biological processes including developmental morphogenesis and tumor invasion. Experimental data can be collected across a wide range of scales, from the population scale to the individual cell scale. Standard continuum or discrete models used in isolation are insufficient to capture this wide range of data. We develop a discrete cellular automata model of invasion with experimentally motivated rules. The cellular automata algorithm is applied to a narrow two-dimensional lattice and simulations reveal the formation of invasion waves moving with constant speed. The simulation results are averaged in one dimension-these data are used to identify the time history of the leading edge to characterize the population-scale wave speed. This allows the relationship between the population-scale wave speed and the cell-scale parameters to be determined. This relationship is analogous to well-known continuum results for Fisher's equation. The cellular automata algorithm also produces individual cell trajectories within the invasion wave that are analogous to cell trajectories obtained with new experimental techniques. Our approach allows both the cell-scale and population-scale properties of invasion to be predicted in a way that is consistent with multiscale experimental data. Furthermore we suggest that the cellular automata algorithm can be used in conjunction with individual data to overcome limitations associated with identifying cell motility mechanisms using continuum models alone.

Coalescence of interacting cell populations

We analyse the coalescence of invasive cell populations by studying both the temporal and steady behaviour of a system of coupled reaction-diffusion equations. This problem is relevant to recent experimental observations of the dynamics of opposingly directed invasion waves of cells. Two cell types, u and v, are considered with the cell motility governed by linear or nonlinear diffusion. The cells proliferate logistically so that the long-term total cell density, u+v approaches a carrying capacity. The steady-state solutions for u and v are denoted u(s) and v(s). The steady solutions are spatially invariant and satisfy u(s)+v(s)=1. However, this expression is underdetermined so the relative proportion of each cell type u(s) and v(s) cannot be determined a priori. Various properties of this model are studied, such as how the relative proportion of u(s) and v(s) depends on the relative motility and relative proliferation rates. The model is analysed using a combination of numerical simulations and a comparison principle. This investigation unearths some novel outcomes regarding the role of overcrowding and cell death in this type of cell migration assay. These observations have relevance to experimental design and interpretation regarding the identification and parameterisation of mechanisms involved in cell invasion.

Multi-scale modeling of a wound-healing cell migration assay

A continuum model and a discrete model are developed to capture the population-scale and cell-scale behavior in a wound-healing cell migration assay created from a scrape wound in a confluent cell monolayer. During wound closure, the cell population forms a sustained traveling wave, with close contact between cells behind the wavefront. Cells exhibit contact inhibition of migration and contact-limited proliferation. The continuum model includes the two dominant mechanisms and characteristics of cell migration and proliferation, using a cell diffusivity function that decreases with cell density and a logistic proliferative growth term. The discrete model arises naturally from the continuum model. Individual cells are simulated as continuous-time random walkers with nearest-neighbor transitions, together with a birth/death process. The migration and proliferation parameters are determined by analysing individual mice 3T3 fibroblast cell trajectories obtained during the development of a confluent cell monolayer and in a wound healing assay. The population-scale model successfully predicts the shape and speed of the traveling wave, while the discrete model is also successful in capturing the contact inhibition of migration effects.

Looking inside an invasion wave of cells using continuum models: proliferation is the key

Recently, a suite of cell migration assays were conducted to investigate the migration of neural crest (NC) cells along the gut during the development of the enteric nervous system (ENS). The NC cells colonise the gastro-intestinal tract as a rostro-caudal wave. Local behaviour was shown to be controlled by position relative to the leading edge of the wavefront. The assays involved chick-quail grafting techniques allowing the total invading population to be considered as a two-species system. A two-species continuum model with logistic proliferation and a migration mechanism is developed here to simulate the chick-quail graft experiments and provide a means of looking at the processes occurring within the invasion wave. Five migration mechanisms are considered--linear diffusion, two cases of nonlinear diffusion, chemokinesis and chemotaxis. The model results agree with the experimental observations, regardless of the specific type of migration mechanism. The results show that NC cell invasion is driven by proliferation and cell motility at the leading edge of the wave. Furthermore, logistic proliferation exerts the dominant control on the system. This observation is confirmed by analysing some simplified invasion models. Once the basic experiments were mathematically replicated, the mathematical models were used in turn to make some predictions that were yet to be experimentally tested. This involved conducting a sensitivity analysis of the system by interrupting the proliferation and/or migration ability of the leading edge. Numerical results show that the system is stable against these changes. Of the three experiments suggested, one was carried out and the experimental results were concordant with the theoretical predictions. The outcome of two other suggested experiments are predicted and left for future experimental validation.