CDEV Subgroup Contributed Talks

Lucy Ham

The University of Melbourne

"Cell fate control in space and time"

Genetically identical cells can adopt distinct, stable states, playing a crucial role in development and tissue organisation. This talk explores the mechanisms driving cell fate decisions, focusing on the interplay between gene regulatory networks and cell-to-cell communication. Using spatial stochastic models that capture fine-scale regulatory dynamics, we demonstrate how feedback loops and paracrine signalling function as switch-like controllers of cell fate, enabling transitions from transient to stable states. We derive mathematical expressions predicting the threshold signalling strength required to trigger phase transitions and establish a fundamental limit on the spatial spread of phenotypic regions. Specifically, we show that the mean region size scales proportionally to the cubic root of the signalling strength, implying that large, stable domains are prohibitively costly to maintain. This trade-off between robustness and signalling precision highlights the constraints organisms must navigate during development to maintain spatial organisation. Our findings provide key insights into the principles governing multicellular patterning and the regulation of tissue structure.

---

Molly Brennan

University College London

"An asymptotic upscaling of transport across bacterial membranes"

Multiscale problems are prevalent in many real world scenarios, especially in biology, where the behaviour of a single microorganism can have considerable impact over lengthscales much larger than its own. In this work we consider the effect of the membrane microstructure of a bacterial cell on the behaviour of concentration profiles of relevant molecules on bacterial and bacterial colony lengthscales. Transport through the outer membrane of gram-negative bacteria is restricted to specific channels and non-specific porins. These provide a size-restricted passageway for small molecules through an otherwise impermeable membrane. The effects of these channels are important, for example, antibiotics must cross the outer membrane in order to effectively target gram-negative bacteria, and quorum sensing molecules must cross the membrane to allow bacterial colonies to coordinate mass phenotypic changes such as the production of virulence factors. In mathematical models this limiting transport mechanism across the membrane is often represented via phenomenological constitutive boundary conditions. In this work, we systematically derive the correct effective boundary conditions to impose across a bacterial membrane in terms of physical channel and porin properties. We use a hybrid mathematical approach, combining multiscale methodology such as asymptotic homogenisation and boundary layer theory with numerical simulations. More broadly, because we consider a generic membrane geometry and do not impose a specific outer problem, the results that we derive have a wide scope of potential applications beyond bacterial membranes, for example, to model water vapour or heat loss through fabrics, or mass transfer through surface coatings in chemical engineering.

---

Augustinas Sukys

The University of Melbourne

"Cell-cycle dependence of bursty gene expression: insights from fitting mechanistic models to single-cell RNA-seq data"

Many genes are expressed in bursts of transcription, associated with alternating active and inactive promoter states. Such transcriptional bursting is characterised by the burst frequency and burst size, which describe how often a burst occurs and how many transcripts are produced per burst. These two burst parameters offer a simple, intuitive and practical quantitative description of bursty gene expression dynamics. However, a transcriptome-wide picture of how the burst frequency and size are modulated due to gene replication and other cell-cycle dependent factors remains missing. To address this, we fit mechanistic models of gene expression to mRNA count data for thousands of mouse genes, obtained by sequencing of single cells whose cell-cycle position has been inferred previously. Although we observe substantial heterogeneity in transcriptional regulation, we find that upon DNA replication, the genome-wide median burst frequency approximately halves, while the median burst size remains mostly unchanged, thus shedding light on the effect of gene dosage compensation. We show that to accurately estimate the bursting kinetics from sequencing data, mechanistic models must explicitly account for gene copy number variation and extrinsic noise due to factors varying across the cell cycle, whereas correcting for technical noise due to imperfect mRNA capture is less critical.

---

Elizabeth Trofimenkoff

University of Lethbridge

"Mathematical modeling of transcription-independent splicing events in human gene expression"

Pre-mRNA often contains introns, which are non-coding sequences that need to be cut out or spliced before translation occurs. The spliceosome, an essential catalyst composed of several proteins with specific sequence affinities, is required for this process. Very long introns must be removed in pieces, a process known as recursive splicing. The experimental literature on the time it takes for the splicing process to occur is inconsistent. Splicing was traditionally believed to be a slow process that could take anywhere from one to tens of minutes per splicing event. However, recent reports suggest that some splicing events occur within a few tens of seconds. We developed the chemical master equation corresponding to the biochemical mechanism of splicing, allowing us to derive the system’s probability distribution, and perform a stability analysis on two conditions based on an unknown association constant parameter associated with the binding step of the scaffolding complex. We also concluded that the distribution of splice times for a single event ranges from a few tens of seconds to a few tens of minutes. Through sensitivity analyses, we have found that the mean splicing time and distribution are almost entirely dependent on the rate at which the spliceosome is activated in the assembly process—i.e. when the U1 and U4 splicing factors dissociate—which confirms that this is the rate limiting step in the catalytic process. Finally, we have examined the distributions of recursive splicing up to six events, and derived analytic solutions for these recursive splicing events in the case where the scaffolding complex strongly binds to the pre-mRNA complex (the condition thought to favor recursive binding), thus providing a model that can be fit to experimental data to in order to evaluate the number of recursive splicing events occurring.

---

Stéphanie Abo

University of Oxford

"Travelling waves in age-structured collective cell migration"

This work examines the interplay between age-structure and migration dynamics in collective cell behaviour. We focus on the integration of cell cycle dynamics with spatial migration, particularly examining the 'go-or-grow' hypothesis in the context of age-dependent processes. Our framework extends classical travelling wave theory to account for the age structure of cell populations, offering new insights into how cell cycle phases influence moving fronts and invasion dynamics. We analyse wave speed characteristics and front dynamics in age-structured systems, addressing a significant gap in current mathematical biology literature. The research provides a novel theoretical foundation for understanding how cell-cycle dependent proliferation and migration behaviours contribute to collective cell dynamics.

---

Gordon R. McNicol

University of Waterloo

"Mechanotransducing structures promote self-driven cell surface patterning"

Cells respond to their local environment through mechanotransduction, converting mechanical signals into a biological response (e.g. cell growth, proliferation or differentiation). The cell cytoskeleton, particularly actomyosin stress fibres (SFs), and focal adhesions (FAs), which bind the cytoskeleton to the extra-cellular matrix (ECM), are central to this process, activating intracellular signalling cascades in response to deformation. We present a novel two-dimensional bio-chemo-mechanical model to describe the development of these structures, coupled through a positive feedback loop, and the associated cell deformation. Building on our previous one-dimensional approach, we similarly employ reaction-diffusion-advection equations to describe the evolution of key scaffolding and signalling proteins, and connect their concentrations to a viscoelastic description of the cell cytoplasm, ECM and adhesions. Further, we now incorporate other key mechanotransducing structures including the stiff cell nucleus, and plasma and cortical membranes. Working in an axisymmetric framework, we employ this model to explain how, dependent upon the mechanical properties of the surrounding ECM, non-uniform patterns of cell striation develop, leading to FA and SF localisation at the cell periphery. Moreover, a linear stability analysis reveals the stability of the axisymmetric configuration to various normal modes of deformation. By identifying non-axisymmetric modes with positive growth rates our model demonstrates a possible mechanism for self-driven surface patterning of cells in vitro.

---

Marc Roussel

University of Lethbridge

"The bacterial dimeric transcription factor NsrR: a case study of a regulatory protein with a large number of states"

In a number of bacteria, nitric oxide (NO) is converted to nitrate by an enzyme called Hmp. In emph{Streptomyces coelicolor}, synthesis of Hmp is in turn controlled by an iron-sulfur protein called NsrR. NsrR represses the transcription of two copies of the emph{hmp} gene in the emph{S. coelicolor} genome, but reaction of NsrR's iron-sulfur cluster with NO causes NsrR to dissociate from the emph{hmp} promoter, thus allowing Hmp to be expressed. While this is a straightforward control mechanism, NsrR is a dimer, and the iron-sulfur cluster in each monomer of NsrR can react with NO several times. Eventually, a repair system restores the NO-damaged iron-sulfur clusters of the dimers. But given that a single reaction with NO is sufficient to cause the NsrR dimer to dissociate from the emph{hmp} promoter, do we need to model the complex chemistry of the dimer, or is a highly simplified model that considers a single NsrR unit and its iron-sulfur cluster sufficient to capture the dynamics of this control system?

---

Paco Castaneda Ruan

The University of Auckland

"Exploring the role of Ca2+ influx in controlling competing oscillatory mechanisms in T cells using ODEs"

Across the spectrum of cell types, the concentration of calcium controls a wide array of cellular functions. These calcium signals, usually in the form of periodic oscillations, play a paramount role in correct cellular activity. T cells are fundamental to the correct behaviour of the immune system. These cells have recently been shown to exhibit two competing oscillatory mechanisms, depending on the influx of extracellular Ca2+. Ca2+ influx is controlled by two molecules, STIM1 and STIM2. When both STIMs are present, T cells showcase sinusoidal Ca2+ oscillations on a raised baseline, but when one of them is absent the nature of the oscillation changes to a mix of Ca2+ spikes and bursting periods. In this talk, we will present an ODE that attempts to explain how these two molecules control the nature of these oscillations in T cells

---