Coupling volume-excluding compartment-based models of diffusion at different scales: Voronoi and pseudo-compartment approaches

Fock-space methods for diffusion: Capturing volume exclusion via fermionic statistics

Volume exclusion and single-file diffusion play an important role at very small scales, such as those associated with molecular machines, ion channels, and transport in zeolites, while introducing fundamental differences compared to Brownian motion, such as changes to the power-law relation between the mean square displacement and time. In this work we map the chemical master equation for excluded diffusion onto a Schrödinger equation via annihilation and creation ladder operators with fermionic statistics, together with linear and symbolic algebra with the resulting Fock-space representation to describe the effect of volume-exclusion processes in finite one-dimensional chains. We contrast the dynamics with the nonexclusive (bosonic) diffusion for crowded, intermediate, and dilute particle populations. Motivated by shuttling in molecular machines, we proceed to investigate differences in exit time distributions introduced by volume exclusion, incorporating the presence of transport bias. More generally, this study demonstrates how one can analyze volume-excluded transport for small stochastic systems, without the need for stochastic simulation and ensemble averaging, simply by considering anticommutation relations and fermionic statistics in a Fock-space representation of the stochastic dynamics.

Accurate and efficient discretizations for stochastic models providing near agent-based spatial resolution at low computational cost

Understanding how cells proliferate, migrate and die in various environments is essential in determining how organisms develop and repair themselves. Continuum mathematical models, such as the logistic equation and the Fisher-Kolmogorov equation, can describe the global characteristics observed in commonly used cell biology assays, such as proliferation and scratch assays. However, these continuum models do not account for single-cell-level mechanics observed in high-throughput experiments. Mathematical modelling frameworks that represent individual cells, often called agent-based models, can successfully describe key single-cell-level features of these assays but are computationally infeasible when dealing with large populations. In this work, we propose an agent-based model with crowding effects that is computationally efficient and matches the logistic and Fisher-Kolmogorov equations in parameter regimes relevant to proliferation and scratch assays, respectively. This stochastic agent-based model allows multiple agents to be contained within compartments on an underlying lattice, thereby reducing the computational storage compared to existing agent-based models that allow one agent per site only. We propose a systematic method to determine a suitable compartment size. Implementing this compartment-based model with this compartment size provides a balance between computational storage, local resolution of agent behaviour and agreement with classical continuum descriptions.

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.

Simulating biological processes: stochastic physics from whole cells to colonies

The last few decades have revealed the living cell to be a crowded spatially heterogeneous space teeming with biomolecules whose concentrations and activities are governed by intrinsically random forces. It is from this randomness, however, that a vast array of precisely timed and intricately coordinated biological functions emerge that give rise to the complex forms and behaviors we see in the biosphere around us. This seemingly paradoxical nature of life has drawn the interest of an increasing number of physicists, and recent years have seen stochastic modeling grow into a major subdiscipline within biological physics. Here we review some of the major advances that have shaped our understanding of stochasticity in biology. We begin with some historical context, outlining a string of important experimental results that motivated the development of stochastic modeling. We then embark upon a fairly rigorous treatment of the simulation methods that are currently available for the treatment of stochastic biological models, with an eye toward comparing and contrasting their realms of applicability, and the care that must be taken when parameterizing them. Following that, we describe how stochasticity impacts several key biological functions, including transcription, translation, ribosome biogenesis, chromosome replication, and metabolism, before considering how the functions may be coupled into a comprehensive model of a 'minimal cell'. Finally, we close with our expectation for the future of the field, focusing on how mesoscopic stochastic methods may be augmented with atomic-scale molecular modeling approaches in order to understand life across a range of length and time scales.

Coarse-graining and hybrid methods for efficient simulation of stochastic multi-scale models of tumour growth

The development of hybrid methodologies is of current interest in both multi-scale modelling and stochastic reaction-diffusion systems regarding their applications to biology. We formulate a hybrid method for stochastic multi-scale models of cells populations that extends the remit of existing hybrid methods for reaction-diffusion systems. Such method is developed for a stochastic multi-scale model of tumour growth, i.e. population-dynamical models which account for the effects of intrinsic noise affecting both the number of cells and the intracellular dynamics. In order to formulate this method, we develop a coarse-grained approximation for both the full stochastic model and its mean-field limit. Such approximation involves averaging out the age-structure (which accounts for the multi-scale nature of the model) by assuming that the age distribution of the population settles onto equilibrium very fast. We then couple the coarse-grained mean-field model to the full stochastic multi-scale model. By doing so, within the mean-field region, we are neglecting noise in both cell numbers (population) and their birth rates (structure). This implies that, in addition to the issues that arise in stochastic-reaction diffusion systems, we need to account for the age-structure of the population when attempting to couple both descriptions. We exploit our coarse-graining model so that, within the mean-field region, the age-distribution is in equilibrium and we know its explicit form. This allows us to couple both domains consistently, as upon transference of cells from the mean-field to the stochastic region, we sample the equilibrium age distribution. Furthermore, our method allows us to investigate the effects of intracellular noise, i.e. fluctuations of the birth rate, on collective properties such as travelling wave velocity. We show that the combination of population and birth-rate noise gives rise to large fluctuations of the birth rate in the region at the leading edge of front, which cannot be accounted for by the coarse-grained model. Such fluctuations have non-trivial effects on the wave velocity. Beyond the development of a new hybrid method, we thus conclude that birth-rate fluctuations are central to a quantitatively accurate description of invasive phenomena such as tumour growth.