Spatial-Stochastic modelling of synthetic gene regulatory networks

Beyond microtubules: The cellular environment at the endoplasmic reticulum attracts proteins to the nucleus, enabling nuclear transport

All proteins are translated in the cytoplasm, yet many, including transcription factors, play vital roles in the nucleus. While previous research has concentrated on molecular motors for the transport of these proteins to the nucleus, recent observations reveal perinuclear accumulation even in the absence of an energy source, hinting at alternative mechanisms. Here, we propose that structural properties of the cellular environment, specifically the endoplasmic reticulum (ER), can promote molecular transport to the perinucleus without requiring additional energy expenditure. Specifically, physical interaction between proteins and the ER impedes their diffusion and leads to their accumulation near the nucleus. This result explains why larger proteins, more frequently interacting with the ER membrane, tend to accumulate at the perinucleus. Interestingly, such diffusion in a heterogeneous environment follows Chapman's law rather than the popular Fick's law. Our findings suggest a novel protein transport mechanism arising solely from characteristics of the intracellular environment.

Reliable and efficient parameter estimation using approximate continuum limit descriptions of stochastic models

Stochastic individual-based mathematical models are attractive for modelling biological phenomena because they naturally capture the stochasticity and variability that is often evident in biological data. Such models also allow us to track the motion of individuals within the population of interest. Unfortunately, capturing this microscopic detail means that simulation and parameter inference can become computationally expensive. One approach for overcoming this computational limitation is to coarse-grain the stochastic model to provide an approximate continuum model that can be solved using far less computational effort. However, coarse-grained continuum models can be biased or inaccurate, particularly for certain parameter regimes. In this work, we combine stochastic and continuum mathematical models in the context of lattice-based models of two-dimensional cell biology experiments by demonstrating how to simulate two commonly used experiments: cell proliferation assays and barrier assays. Our approach involves building a simple statistical model of the discrepancy between the expensive stochastic model and the associated computationally inexpensive coarse-grained continuum model. We form this statistical model based on a limited number of expensive stochastic model evaluations at design points sampled from a user-chosen distribution, corresponding to a computer experiment design problem. With straightforward design point selection schemes, we show that using the statistical model of the discrepancy in tandem with the computationally inexpensive continuum model allows us to carry out prediction and inference while correcting for biases and inaccuracies due to the continuum approximation. We demonstrate this approach by simulating a proliferation assay, where the continuum limit model is the well-known logistic ordinary differential equation, as well as a barrier assay where the continuum limit model is closely related to the well-known Fisher-KPP partial differential equation. We construct an approximate likelihood function for parameter inference, both with and without discrepancy correction terms. Using maximum likelihood estimation, we provide point estimates of the unknown parameters, and use the profile likelihood to characterise the uncertainty in these estimates and form approximate confidence intervals. For the range of inference problems considered, working with the continuum limit model alone leads to biased parameter estimation and confidence intervals with poor coverage. In contrast, incorporating correction terms arising from the statistical model of the model discrepancy allows us to recover the parameters accurately with minimal computational overhead. The main tradeoff is that the associated confidence intervals are typically broader, reflecting the additional uncertainty introduced by the approximation process. All algorithms required to replicate the results in this work are written in the open source Julia language and are available at GitHub.

Mathematical Modelling of p53 Signalling during DNA Damage Response: A Survey

No gene has garnered more interest than p53 since its discovery over 40 years ago. In the last two decades, thanks to seminal work from Uri Alon and Ghalit Lahav, p53 has defined a truly synergistic topic in the field of mathematical biology, with a rich body of research connecting mathematic endeavour with experimental design and data. In this review we survey and distill the extensive literature of mathematical models of p53. Specifically, we focus on models which seek to reproduce the oscillatory dynamics of p53 in response to DNA damage. We review the standard modelling approaches used in the field categorising them into three types: time delay models, spatial models and coupled negative-positive feedback models, providing sample model equations and simulation results which show clear oscillatory dynamics. We discuss the interplay between mathematics and biology and show how one informs the other; the deep connections between the two disciplines has helped to develop our understanding of this complex gene and paint a picture of its dynamical response. Although yet more is to be elucidated, we offer the current state-of-the-art understanding of p53 response to DNA damage.

MeSCoT: The tool for quantitative trait simulation through the mechanistic modelling of genes' regulatory interactions

This work represents a novel mechanistic approach to simulate and study genomic networks with accompanying regulatory interactions and complex mechanisms of quantitative trait formation. The approach implemented in MeSCoT software is conceptually based on the omnigenic genetic model of quantitative (complex) trait, and closely imitates the basic in vivo mechanisms of quantitative trait realization. The software provides a framework to study molecular mechanisms of gene-by-gene and gene-by-environment interactions underlying quantitative trait’s realization and allows detailed mechanistic studies of impact of genetic and phenotypic variance on gene regulation. MeSCoT performs a detailed simulation of genes’ regulatory interactions for variable genomic architectures and generates complete set of transcriptional and translational data together with simulated quantitative trait values. Such data provide opportunities to study, for example, verification of novel statistical methods aiming to integrate intermediate phenotypes together with final phenotype in quantitative genetic analyses or to investigate novel approaches for exploiting gene-by-gene and gene-by-environment interactions.