Stochastic multi-scale models of competition within heterogeneous cellular populations: Simulation methods and mean-field analysis

Learning differential equation models from stochastic agent-based model simulations

Agent-based models provide a flexible framework that is frequently used for modelling many biological systems, including cell migration, molecular dynamics, ecology and epidemiology. Analysis of the model dynamics can be challenging due to their inherent stochasticity and heavy computational requirements. Common approaches to the analysis of agent-based models include extensive Monte Carlo simulation of the model or the derivation of coarse-grained differential equation models to predict the expected or averaged output from the agent-based model. Both of these approaches have limitations, however, as extensive computation of complex agent-based models may be infeasible, and coarse-grained differential equation models can fail to accurately describe model dynamics in certain parameter regimes. We propose that methods from the equation learning field provide a promising, novel and unifying approach for agent-based model analysis. Equation learning is a recent field of research from data science that aims to infer differential equation models directly from data. We use this tutorial to review how methods from equation learning can be used to learn differential equation models from agent-based model simulations. We demonstrate that this framework is easy to use, requires few model simulations, and accurately predicts model dynamics in parameter regions where coarse-grained differential equation models fail to do so. We highlight these advantages through several case studies involving two agent-based models that are broadly applicable to biological phenomena: a birth-death-migration model commonly used to explore cell biology experiments and a susceptible-infected-recovered model of infectious disease spread.

Approaches to combat hypoxia in cancer therapy and the potential for in silico models in their evaluation

AIM

The negative impact of tumour hypoxia on cancer treatment outcome has been long-known, yet there has been little success combating it. This paper investigates the potential role of in silico modelling to help test emerging hypoxia-targeting treatments in cancer therapy.

METHODS

A Medline search was undertaken on the current landscape of in silico models that simulate cancer therapy and evaluate their ability to test hypoxia-targeting treatments. Techniques and treatments to combat tumour hypoxia and their current challenges are also presented.

RESULTS

Hypoxia-targeting treatments include tumour reoxygenation, hypoxic cell radiosensitization with nitroimidazoles, hypoxia-activated prodrugs and molecular targeting. Their main challenges are toxicity and not achieving adequate delivery to hypoxic regions of the tumour. There is promising research toward combining two or more of these techniques. Different types of in silico therapy models have been developed ranging from temporal to spatial and from stochastic to deterministic models. Numerous models have compared the effectiveness of different radiotherapy fractionation schedules for controlling hypoxic tumours. Similarly, models could help identify and optimize new treatments for overcoming hypoxia that utilize novel hypoxia-targeting technology.

CONCLUSION

Current therapy models should attempt to incorporate more sophisticated modelling of tumour angiogenesis/vasculature and vessel perfusion in order to become more useful for testing hypoxia-targeting treatments, which typically rely upon the tumour vasculature for delivery of additional oxygen, (pro)drugs and nanoparticles.

Age Structure Can Account for Delayed Logistic Proliferation of Scratch Assays

Scratch assays are in vitro methods for studying cell migration. In these experiments, a scratch is made on a cell monolayer and recolonisation of the scratched region is imaged to quantify cell migration rates. Typically, scratch assays are modelled by reaction diffusion equations depicting cell migration by Fickian diffusion and proliferation by a logistic term. In a recent paper (Jin et al. in Bull Math Biol 79(5):1028-1050, 2017), the authors observed experimentally that during the early stage of the recolonisation process, there is a disturbance phase where proliferation is not logistic, and this is followed by a growth phase where proliferation appears to be logistic. The authors did not identify the precise mechanism that causes the disturbance phase but showed that ignoring it can lead to incorrect parameter estimates. The aim of this work is to show that a nonlinear age-structured population model can account for the two phases of proliferation in scratch assays. The model consists of an age-structured cell cycle model of a cell population, coupled with an ordinary differential equation describing the resource concentration dynamics in the substrate. The model assumes a resource-dependent cell cycle threshold age, above which cells are able to proliferate. By studying the dynamics of the full system in terms of the subpopulations of cells that can proliferate and the ones that can not, we are able to find conditions under which the model captures the two-phase behaviour. Through numerical simulations, we are able to show that the interplay between the resource concentration in the substrate and the cell subpopulations dynamics can explain the biphasic dynamics.

A parasitism-mutualism-predation model consisting of crows, cuckoos and cats with stage-structure and maturation delays on crows and cuckoos

In this paper, a parasitism-mutualism-predation model is proposed to investigate the dynamics of multi-interactions among cuckoos, crows and cats with stage-structure and maturation time delays on cuckoos and crows. The crows permit the cuckoos to parasitize their nestlings (eggs) on the crow chicks (eggs). In return, the cuckoo nestlings produce a malodorous cloacal secretion to protect the crow chicks from predation by the cats, which is apparently beneficial to both the crow and cuckoo population. The multi-interactions, i.e., parasitism and mutualism between the cuckoos (nestlings) and crows (chicks), predation between the cats and crow chicks are modeled both by Holling-type II and Beddington-DeAngelis-type functional responses. The existence of positive equilibria of three subsystems of the model are discussed. The criteria for the global stability of the trivial equilibrium are established by the Krein-Rutman theorem and other analysis methods. Moreover, the threshold dynamics for the coexistence and weak persistence of the model are obtained, and we show, both analytically and numerically, that the stabilities of the interior equilibria may change with the increasing maturation time delays. We find there exists an evident difference in the dynamical properties of the parasitism-mutualism-predation model based on whether or not we consider the effects of stage-structure and maturation time delays on cuckoos and crows. Inclusion of stage structure results in many varied dynamical complexities which are difficult to encompass without this inclusion.

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.