-
Approximate solutions of a general stochastic velocity-jump process subject to discrete-time noisy observations
Authors:
Arianna Ceccarelli,
Alexander P. Browning,
Ruth E. Baker
Abstract:
Advances in experimental techniques allow the collection of high-space-and-time resolution data that track individual motile entities over time. This poses the question of how to use these data to efficiently and effectively calibrate motion models. However, typical mathematical models often overlook the inherent aspects of data collection, such as the discreteness and the experimental noise of th…
▽ More
Advances in experimental techniques allow the collection of high-space-and-time resolution data that track individual motile entities over time. This poses the question of how to use these data to efficiently and effectively calibrate motion models. However, typical mathematical models often overlook the inherent aspects of data collection, such as the discreteness and the experimental noise of the measured locations. In this paper, we focus on velocity-jump models suitable to describe single-agent motion in one spatial dimension, characterised by successive Markovian transitions between a finite network of $n$ states, each with a specified velocity and a fixed rate of switching to every other state. Since the problem of finding the exact distributions of discrete-time noisy data is generally intractable, we derive a series of approximations for the data distributions and compare them to in-silico data generated by the models using four example network structures. These comparisons suggest that the approximations are accurate given sufficiently infrequent state switching, or equivalently, a sufficiently high data collection frequency. Moreover, for infrequent switching, the PDFs comparisons highlight the importance of accounting for the correlation between subsequent measured locations, due to the likely permanence in the state visited in the previous measurement. The approximate distributions computed can be used for fast parameter inference and model selection between a range of velocity-jump models using single-agent tracking data.
△ Less
Submitted 5 July, 2024; v1 submitted 28 June, 2024;
originally announced June 2024.
-
Quantifying cell cycle regulation by tissue crowding
Authors:
Carles Falcó,
Daniel J. Cohen,
José A. Carrillo,
Ruth E. Baker
Abstract:
The spatiotemporal coordination and regulation of cell proliferation is fundamental in many aspects of development and tissue maintenance. Cells have the ability to adapt their division rates in response to mechanical constraints, yet we do not fully understand how cell proliferation regulation impacts cell migration phenomena. Here, we present a minimal continuum model of cell migration with cell…
▽ More
The spatiotemporal coordination and regulation of cell proliferation is fundamental in many aspects of development and tissue maintenance. Cells have the ability to adapt their division rates in response to mechanical constraints, yet we do not fully understand how cell proliferation regulation impacts cell migration phenomena. Here, we present a minimal continuum model of cell migration with cell cycle dynamics, which includes density-dependent effects and hence can account for cell proliferation regulation. By combining minimal mathematical modelling, Bayesian inference, and recent experimental data, we quantify the impact of tissue crowding across different cell cycle stages in epithelial tissue expansion experiments. Our model suggests that cells sense local density and adapt cell cycle progression in response, during G1 and the combined S/G2/M phases, providing an explicit relationship between each cell cycle stage duration and local tissue density, which is consistent with several experimental observations. Finally, we compare our mathematical model predictions to different experiments studying cell cycle regulation and present a quantitative analysis on the impact of density-dependent regulation on cell migration patterns. Our work presents a systematic approach for investigating and analysing cell cycle data, providing mechanistic insights into how individual cells regulate proliferation, based on population-based experimental measurements.
△ Less
Submitted 24 April, 2024; v1 submitted 16 January, 2024;
originally announced January 2024.
-
Energy translation symmetries and dynamics of separable autonomous two-dimensional ODEs
Authors:
Johannes G. Borgqvist,
Fredrik Ohlsson,
Ruth E. Baker
Abstract:
We study symmetries in the phase plane for separable, autonomous two-state systems of ordinary differential equations (ODEs). We prove two main theoretical results concerning the existence and non-triviality of two orthogonal symmetries for such systems. In particular, we show that these symmetries correspond to translations in the internal energy of the system, and describe their action on soluti…
▽ More
We study symmetries in the phase plane for separable, autonomous two-state systems of ordinary differential equations (ODEs). We prove two main theoretical results concerning the existence and non-triviality of two orthogonal symmetries for such systems. In particular, we show that these symmetries correspond to translations in the internal energy of the system, and describe their action on solution trajectories in the phase plane. In addition, we apply recent results establishing how phase plane symmetries can be extended to incorporate temporal dynamics to these energy translation symmetries. Subsequently, we apply our theoretical results to the analysis of three models from the field of mathematical biology: a canonical biological oscillator model, the Lotka--Volterra (LV) model describing predator-prey dynamics, and the SIR model describing the spread of a disease in a population. We describe the energy translation symmetries in detail, including their action on biological observables of the models, derive analytic expressions for the extensions to the time domain, and discuss their action on solution trajectories.
△ Less
Submitted 15 August, 2023; v1 submitted 20 February, 2023;
originally announced February 2023.
-
Symmetries of systems of first order ODEs: Symbolic symmetry computations, mechanistic model construction and applications in biology
Authors:
Johannes Borgqvist,
Fredrik Ohlsson,
Ruth E. Baker
Abstract:
We discuss the role and merits of symmetry methods for the analysis of biological systems. In particular, we consider systems of first order ordinary differential equations and provide a comprehensive review of the geometrical foundations pertinent to symmetries of such systems. Subsequently, we present an algorithm for finding infinitesimal generators of symmetries for systems with rational react…
▽ More
We discuss the role and merits of symmetry methods for the analysis of biological systems. In particular, we consider systems of first order ordinary differential equations and provide a comprehensive review of the geometrical foundations pertinent to symmetries of such systems. Subsequently, we present an algorithm for finding infinitesimal generators of symmetries for systems with rational reaction terms, and an open-source implementation of the algorithm using symbolic computations. We discuss two complementary perspectives on symmetries in mechanistic modelling; as tools for the analysis of a given model or as a geometrical principle for incorporating biological properties in the construction of new models. Through numerous examples of relevance to modelling in biology we demonstrate the different uses of symmetry methods, and also discuss how to infer symmetries from experimental data.
△ Less
Submitted 10 February, 2022;
originally announced February 2022.
-
Profile likelihood analysis for a stochastic model of diffusion in heterogeneous media
Authors:
Matthew J Simpson,
Alexander P Browning,
Christopher Drovandi,
Elliot J Carr,
Oliver J Maclaren,
Ruth E Baker
Abstract:
We compute profile likelihoods for a stochastic model of diffusive transport motivated by experimental observations of heat conduction in layered skin tissues. This process is modelled as a random walk in a layered one-dimensional material, where each layer has a distinct particle hopping rate. Particles are released at some location, and the duration of time taken for each particle to reach an ab…
▽ More
We compute profile likelihoods for a stochastic model of diffusive transport motivated by experimental observations of heat conduction in layered skin tissues. This process is modelled as a random walk in a layered one-dimensional material, where each layer has a distinct particle hopping rate. Particles are released at some location, and the duration of time taken for each particle to reach an absorbing boundary is recorded. To explore whether this data can be used to identify the hopping rates in each layer, we compute various profile likelihoods using two methods: first, an exact likelihood is evaluated using a relatively expensive Markov chain approach; and, second we form an approximate likelihood by assuming the distribution of exit times is given by a Gamma distribution whose first two moments match the expected moments from the continuum limit description of the stochastic model. Using the exact and approximate likelihoods we construct various profile likelihoods for a range of problems. In cases where parameter values are not identifiable, we make progress by re-interpreting those data with a reduced model with a smaller number of layers.
△ Less
Submitted 9 March, 2021; v1 submitted 6 November, 2020;
originally announced November 2020.
-
Crowded transport within networked representations of complex geometries
Authors:
Daniel B. Wilson,
Francis G. Woodhouse,
Matthew J. Simpson,
Ruth E. Baker
Abstract:
Transport in crowded, complex environments occurs across many spatial scales. Geometric restrictions can hinder the motion of individuals and, combined with crowding between individuals, can have drastic effects on global transport phenomena. However, in general, the interplay between crowding and geometry in complex real-life environments is poorly understood. Existing analytical methodologies ar…
▽ More
Transport in crowded, complex environments occurs across many spatial scales. Geometric restrictions can hinder the motion of individuals and, combined with crowding between individuals, can have drastic effects on global transport phenomena. However, in general, the interplay between crowding and geometry in complex real-life environments is poorly understood. Existing analytical methodologies are not always readily extendable to heterogeneous environments: in these situations predictions of crowded transport behaviour within heterogeneous environments rely on computationally intensive mesh-based approaches. Here, we take a different approach by employing networked representations of complex environments to provide an efficient framework within which the interactions between networked geometry and crowding can be explored. We demonstrate how the framework can be used to: extract detailed information at the level of the whole population or an individual within it; identify the topological features of environments that enable accurate prediction of transport phenomena; and, provide insights into the design of optimal environments.
△ Less
Submitted 10 August, 2021; v1 submitted 24 June, 2020;
originally announced June 2020.
-
Effects of different discretisations of the Laplacian upon stochastic simulations of reaction-diffusion systems on both static and growing domains
Authors:
Bartosz J. Bartmanski,
Ruth E. Baker
Abstract:
By discretising space into compartments and letting system dynamics be governed by the reaction-diffusion master equation, it is possible to derive and simulate a stochastic model of reaction and diffusion on an arbitrary domain. However, there are many implementation choices involved in this process, such as the choice of discretisation and method of derivation of the diffusive jump rates, and it…
▽ More
By discretising space into compartments and letting system dynamics be governed by the reaction-diffusion master equation, it is possible to derive and simulate a stochastic model of reaction and diffusion on an arbitrary domain. However, there are many implementation choices involved in this process, such as the choice of discretisation and method of derivation of the diffusive jump rates, and it is not clear a priori how these affect model predictions. To shed light on this issue, in this work we explore how a variety of discretisations and method for derivation of the diffusive jump rates affect the outputs of stochastic simulations of reaction-diffusion models, in particular using Turing's model of pattern formation as a key example. We consider both static and uniformly growing domains and demonstrate that, while only minor differences are observed for simple reaction-diffusion systems, there can be vast differences in model predictions for systems that include complicated reaction kinetics, such as Turing's model of pattern formation. Our work highlights that care must be taken in using the reaction-diffusion master equation to make predictions as to the dynamics of stochastic reaction-diffusion systems.
△ Less
Submitted 26 November, 2019;
originally announced November 2019.
-
Topology-dependent density optima for efficient simultaneous network exploration
Authors:
Daniel B. Wilson,
Ruth E. Baker,
Francis G. Woodhouse
Abstract:
A random search process in a networked environment is governed by the time it takes to visit every node, termed the cover time. Often, a networked process does not proceed in isolation but competes with many instances of itself within the same environment. A key unanswered question is how to optimise this process: how many concurrent searchers can a topology support before the benefits of parallel…
▽ More
A random search process in a networked environment is governed by the time it takes to visit every node, termed the cover time. Often, a networked process does not proceed in isolation but competes with many instances of itself within the same environment. A key unanswered question is how to optimise this process: how many concurrent searchers can a topology support before the benefits of parallelism are outweighed by competition for space? Here, we introduce the searcher-averaged parallel cover time (APCT) to quantify these economies of scale. We show that the APCT of the networked symmetric exclusion process is optimised at a searcher density that is well predicted by the spectral gap. Furthermore, we find that non-equilibrium processes, realised through the addition of bias, can support significantly increased density optima. Our results suggest novel hybrid strategies of serial and parallel search for efficient information gathering in social interaction and biological transport networks.
△ Less
Submitted 15 March, 2018; v1 submitted 25 September, 2017;
originally announced September 2017.