Modelling the movement of interacting cell populations: A moment dynamics approach

Using a probabilistic approach to derive a two-phase model of flow-induced cell migration

Interstitial fluid flow is a feature of many solid tumors. In vitro experiments have shown that such fluid flow can direct tumor cell movement upstream or downstream depending on the balance between the competing mechanisms of tensotaxis (cell migration up stress gradients) and autologous chemotaxis (downstream cell movement in response to flow-induced gradients of self-secreted chemoattractants). In this work we develop a probabilistic-continuum, two-phase model for cell migration in response to interstitial flow. We use a kinetic description for the cell velocity probability density function, and model the flow-dependent mechanical and chemical stimuli as forcing terms that bias cell migration upstream and downstream. Using velocity-space averaging, we reformulate the model as a system of continuum equations for the spatiotemporal evolution of the cell volume fraction and flux in response to forcing terms that depend on the local direction and magnitude of the mechanochemical cues. We specialize our model to describe a one-dimensional cell layer subject to fluid flow. Using a combination of numerical simulations and asymptotic analysis, we delineate the parameter regime where transitions from downstream to upstream cell migration occur. As has been observed experimentally, the model predicts downstream-oriented chemotactic migration at low cell volume fractions, and upstream-oriented tensotactic migration at larger volume fractions. We show that the locus of the critical volume fraction, at which the system transitions from downstream to upstream migration, is dominated by the ratio of the rate of chemokine secretion and advection. Our model also predicts that, because the tensotactic stimulus depends strongly on the cell volume fraction, upstream, tensotaxis-dominated migration occurs only transiently when the cells are initially seeded, and transitions to downstream, chemotaxis-dominated migration occur at later times due to the dispersive effect of cell diffusion.

Agent-Based and Continuum Models for Spatial Dynamics of Infection by Oncolytic Viruses

The use of oncolytic viruses as cancer treatment has received considerable attention in recent years, however the spatial dynamics of this viral infection is still poorly understood. We present here a stochastic agent-based model describing infected and uninfected cells for solid tumours, which interact with viruses in the absence of an immune response. Two kinds of movement, namely undirected random and pressure-driven movements, are considered: the continuum limit of the models is derived and a systematic comparison between the systems of partial differential equations and the individual-based model, in one and two dimensions, is carried out. In the case of undirected movement, a good agreement between agent-based simulations and the numerical and well-known analytical results for the continuum model is possible. For pressure-driven motion, instead, we observe a wide parameter range in which the infection of the agents remains confined to the center of the tumour, even though the continuum model shows traveling waves of infection; outcomes appear to be more sensitive to stochasticity and uninfected regions appear harder to invade, giving rise to irregular, unpredictable growth patterns. Our results show that the presence of spatial constraints in tumours' microenvironments limiting free expansion has a very significant impact on virotherapy. Outcomes for these tumours suggest a notable increase in variability. All these aspects can have important effects when designing individually tailored therapies where virotherapy is included.

Regulated bacterial interaction networks: A mathematical framework to describe competitive growth under inclusion of metabolite cross-feeding

When bacterial species with the same resource preferences share the same growth environment, it is commonly believed that direct competition will arise. A large variety of competition and more general 'interaction' models have been formulated, but what is currently lacking are models that link monoculture growth kinetics and community growth under inclusion of emerging biological interactions, such as metabolite cross-feeding. In order to understand and mathematically describe the nature of potential cross-feeding interactions, we design experiments where two bacterial species Pseudomonas putida and Pseudomonas veronii grow in liquid medium either in mono- or as co-culture in a resource-limited environment. We measure population growth under single substrate competition or with double species-specific substrates (substrate 'indifference'), and starting from varying cell ratios of either species. Using experimental data as input, we first consider a mean-field model of resource-based competition, which captures well the empirically observed growth rates for monocultures, but fails to correctly predict growth rates in co-culture mixtures, in particular for skewed starting species ratios. Based on this, we extend the model by cross-feeding interactions where the consumption of substrate by one consumer produces metabolites that in turn are resources for the other consumer, thus leading to positive feedback in the species system. Two different cross-feeding options were considered, which either lead to constant metabolite cross-feeding, or to a regulated form, where metabolite utilization is activated with rates according to either a threshold or a Hill function, dependent on metabolite concentration. Both mathematical proof and experimental data indicate regulated cross-feeding to be the preferred model to constant metabolite utilization, with best co-culture growth predictions in case of high Hill coefficients, close to binary (on/off) activation states. This suggests that species use the appearing metabolite concentrations only when they are becoming high enough; possibly as a consequence of their lower energetic content than the primary substrate. Metabolite sharing was particularly relevant at unbalanced starting cell ratios, causing the minority partner to proliferate more than expected from the competitive substrate because of metabolite release from the majority partner. This effect thus likely quells immediate substrate competition and may be important in natural communities with typical very skewed relative taxa abundances and slower-growing taxa. In conclusion, the regulated bacterial interaction network correctly describes species substrate growth reactions in mixtures with few kinetic parameters that can be obtained from monoculture growth experiments.

On predicting heterogeneity in nanoparticle dosage

Nanoparticles are increasingly employed as a vehicle for the targeted delivery of therapeutics to specific cell types. However, much remains to be discovered about the fundamental biology that dictates the interactions between nanoparticles and cells. Accordingly, few nanoparticle-based targeted therapeutics have succeeded in clinical trials. One element that hinders our understanding of nanoparticle-cell interactions is the presence of heterogeneity in nanoparticle dosage data obtained from standard experiments. It is difficult to distinguish between heterogeneity that arises from stochasticity in nanoparticle-cell interactions, and that which arises from heterogeneity in the cell population. Mathematical investigations have revealed that both sources of heterogeneity contribute meaningfully to the heterogeneity in nanoparticle dosage. However, these investigations have relied on simplified models of nanoparticle internalisation. Here we present a stochastic mathematical model of nanoparticle internalisation that incorporates a suite of relevant biological phenomena such as multistage internalisation, cell division, asymmetric nanoparticle inheritance and nanoparticle saturation. Critically, our model provides information about nanoparticle dosage at an individual cell level. We perform model simulations to examine the influence of specific biological phenomena on the heterogeneity in nanoparticle dosage in the absence of heterogeneity in the cell population. Under certain modelling assumptions, we derive analytic approximations of the nanoparticle dosage distribution. We demonstrate that the analytic approximations are accurate, and show that nanoparticle dosage can be described by a Poisson mixture distribution with rate parameters that are a function of Beta-distributed random variables. We discuss the implications of the analytic results with respect to parameter estimation and model identifiability from standard experimental data. Finally, we highlight extensions and directions for future research.

Predicting population extinction in lattice-based birth-death-movement models

The question of whether a population will persist or go extinct is of key interest throughout ecology and biology. Various mathematical techniques allow us to generate knowledge regarding individual behaviour, which can be analysed to obtain predictions about the ultimate survival or extinction of the population. A common model employed to describe population dynamics is the lattice-based random walk model with crowding (exclusion). This model can incorporate behaviour such as birth, death and movement, while including natural phenomena such as finite size effects. Performing sufficiently many realizations of the random walk model to extract representative population behaviour is computationally intensive. Therefore, continuum approximations of random walk models are routinely employed. However, standard continuum approximations are notoriously incapable of making accurate predictions about population extinction. Here, we develop a new continuum approximation, the state-space diffusion approximation, which explicitly accounts for population extinction. Predictions from our approximation faithfully capture the behaviour in the random walk model, and provides additional information compared to standard approximations. We examine the influence of the number of lattice sites and initial number of individuals on the long-term population behaviour, and demonstrate the reduction in computation time between the random walk model and our approximation.

A Brownian dynamics tumor progression simulator with application to glioblastoma

Tumor progression modeling offers the potential to predict tumor-spreading behavior to improve prognostic accuracy and guide therapy development. Common simulation methods include continuous reaction-diffusion (RD) approaches that capture mean spatio-temporal tumor spreading behavior and discrete agent-based (AB) approaches which capture individual cell events such as proliferation or migration. The brain cancer glioblastoma (GBM) is especially appropriate for such proliferation-migration modeling approaches because tumor cells seldom metastasize outside of the central nervous system and cells are both highly proliferative and migratory. In glioblastoma research, current RD estimates of proliferation and migration parameters are derived from computed tomography or magnetic resonance images. However, these estimates of glioblastoma cell migration rates, modeled as a diffusion coefficient, are approximately 1-2 orders of magnitude larger than single-cell measurements in animal models of this disease. To identify possible sources for this discrepancy, we evaluated the fundamental RD simulation assumptions that cells are point-like structures that can overlap. To give cells physical size (~10 μm), we used a Brownian dynamics approach that simulates individual single-cell diffusive migration, growth, and proliferation activity via a gridless, off-lattice, AB method where cells can be prohibited from overlapping each other. We found that for realistic single-cell parameter growth and migration rates, a non-overlapping model gives rise to a jammed configuration in the center of the tumor and a biased outward diffusion of cells in the tumor periphery, creating a quasi-ballistic advancing tumor front. The simulations demonstrate that a fast-progressing tumor can result from minimally diffusive cells, but at a rate that is still dependent on single-cell diffusive migration rates. Thus, modeling with the assumption of physically-grounded volume conservation can account for the apparent discrepancy between estimated and measured diffusion of GBM cells and provide a new theoretical framework that naturally links single-cell growth and migration dynamics to tumor-level progression.

Continuum approximations for lattice-free multi-species models of collective cell migration

Cell migration within tissues involves the interaction of many cells from distinct subpopulations. In this work, we present a discrete model of collective cell migration where the motion of individual cells is driven by random forces, short range repulsion forces to mimic crowding, and longer range attraction forces to mimic adhesion. This discrete model can be used to simulate a population of cells that is composed of K ≥ 1 distinct subpopulations. To analyse the discrete model we formulate a hierarchy of moment equations that describe the spatial evolution of the density of agents, pairs of agents, triplets of agents, and so forth. To solve the hierarchy of moment equations we introduce two forms of closure: (i) the mean field approximation, which effectively assumes that the distributions of individual agents are independent; and (ii) a moment dynamics description that is based on the Kirkwood superposition approximation. The moment dynamics description provides an approximate way of incorporating spatial patterns, such as agent clustering, into the continuum description. Comparing the performance of the two continuum descriptions confirms that both perform well when adhesive forces are sufficiently weak. In contrast, the moment dynamics description outperforms the mean field model when adhesive forces are sufficiently large. This is a first attempt to provide an accurate continuum description of a lattice-free, multi-species model of collective cell migration.

Co-operation, Competition and Crowding: A Discrete Framework Linking Allee Kinetics, Nonlinear Diffusion, Shocks and Sharp-Fronted Travelling Waves

Invasion processes are ubiquitous throughout cell biology and ecology. During invasion, individuals can become isolated from the bulk population and behave differently. We present a discrete, exclusion-based description of the birth, death and movement of individuals. The model distinguishes between individuals that are part of, or are isolated from, the bulk population by imposing different rates of birth, death and movement. This enables the simulation of various co-operative or competitive mechanisms, where there is either a positive or negative benefit associated with being part of the bulk population, respectively. The mean-field approximation of the discrete process gives rise to 22 different classes of partial differential equation, which can include Allee kinetics and nonlinear diffusion. Here we examine the ability of each class of partial differential equation to support travelling wave solutions and interpret the long time behaviour in terms of the individual-level parameters. For the first time we show that the strong Allee effect and nonlinear diffusion can result in shock-fronted travelling waves. We also demonstrate how differences in group and individual motility rates can influence the persistence of a population and provide conditions for the successful invasion of a population.

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.