Turing pattern formation on the sphere is robust to the removal of a hole.

The formation of buds on the cell membrane of budding yeast cells is thought to be driven by reactions and diffusion involving the protein Cdc42. These processes can be described by a coupled system of partial differential equations known as the Schnakenberg system. The Schnakenberg system is known to exhibit diffusion-driven pattern formation, thus providing a mechanism for bud formation. However, it is not known how the accumulation of bud scars on the cell membrane affect the ability of the Schnakenberg system to form patterns. We have approached this problem by modelling a bud scar on the cell membrane with a hole on the sphere. We have studied how the spectrum of the Laplace-Beltrami operator, which determines the resulting pattern, is affected by the size of the hole, and by numerically solving the Schnakenberg system on a sphere with a hole using the finite element method. Both theoretical predictions and numerical solutions show that pattern formation is robust to the introduction of a bud scar of considerable size, which lends credence to the hypothesis that bud formation is driven by diffusion-driven instability.

Occam's razor gets a new edge: the use of symmetries in model selection.

We demonstrate the power of using symmetries for model selection in the context of mechanistic modelling. We analyse two different models called the power law model (PLM) and the immunological model (IM) describing the increase in cancer risk with age, due to mutation accumulation or immunosenescence, respectively. The IM fits several cancer types better than the PLM implying that it would be selected based on minimizing residuals. However, recently a symmetry-based method for model selection has been developed, which has been successfully used in an in silico setting to find the correct model when traditional model fitting has failed. Here, we apply this method in a real-world setting to investigate the mechanisms of carcinogenesis. First, we derive distinct symmetry transformations of the two models and then we select the model which not only fits the original data but is also invariant under transformations by its symmetry. Contrary to the initial conclusion, we conclude that the PLM realistically describes the mechanism underlying the colon cancer dataset. These conclusions agree with experimental knowledge, and this work demonstrates how a model selection criterion based on biological properties can be implemented using symmetries.

Cell polarisation in a bulk-surface model can be driven by both classic and non-classic Turing instability.

The GTPase Cdc42 is the master regulator of eukaryotic cell polarisation. During this process, the active form of Cdc42 is accumulated at a particular site on the cell membrane called the pole. It is believed that the accumulation of the active Cdc42 resulting in a pole is driven by a combination of activation-inactivation reactions and diffusion. It has been proposed using mathematical modelling that this is the result of diffusion-driven instability, originally proposed by Alan Turing. In this study, we developed, analysed and validated a 3D bulk-surface model of the dynamics of Cdc42. We show that the model can undergo both classic and non-classic Turing instability by deriving necessary conditions for which this occurs and conclude that the non-classic case can be viewed as a limit case of the classic case of diffusion-driven instability. Using three-dimensional Spatio-temporal simulation we predicted pole size and time to polarisation, suggesting that cell polarisation is mainly driven by the reaction strength parameter and that the size of the pole is determined by the relative diffusion.

The synergy of damage repair and retention promotes rejuvenation and prolongs healthy lifespans in cell lineages.

Damaged proteins are inherited asymmetrically during cell division in the yeast Saccharomyces cerevisiae, such that most damage is retained within the mother cell. The consequence is an ageing mother and a rejuvenated daughter cell with full replicative potential. Daughters of old and damaged mothers are however born with increasing levels of damage resulting in lowered replicative lifespans. Remarkably, these prematurely old daughters can give rise to rejuvenated cells with low damage levels and recovered lifespans, called second-degree rejuvenation. We aimed to investigate how damage repair and retention together can promote rejuvenation and at the same time ensure low damage levels in mother cells, reflected in longer health spans. We developed a dynamic model for damage accumulation over successive divisions in individual cells as part of a dynamically growing cell lineage. With detailed knowledge about single-cell dynamics and relationships between all cells in the lineage, we can infer how individual damage repair and retention strategies affect the propagation of damage in the population. We show that damage retention lowers damage levels in the population by reducing the variability across the lineage, and results in larger population sizes. Repairing damage efficiently in early life, as opposed to investing in repair when damage has already accumulated, counteracts accelerated ageing caused by damage retention. It prolongs the health span of individual cells which are moreover less prone to stress. In combination, damage retention and early investment in repair are beneficial for healthy ageing in yeast cell populations.

Symmetry structures in dynamic models of biochemical systems.

Understanding the complex interactions of biochemical processes underlying human disease represents the holy grail of systems biology. When processes are modelled in ordinary differential equation (ODE) fashion, the most common tool for their analysis is linear stability analysis where the long-term behaviour of the model is determined by linearizing the system around its steady states. However, this asymptotic behaviour is often insufficient for completely determining the structure of the underlying system. A complementary technique for analysing a system of ODEs is to consider the set of symmetries of its solutions. Symmetries provide a powerful concept for the development of mechanistic models by describing structures corresponding to the underlying dynamics of biological systems. To demonstrate their capability, we consider symmetries of the nonlinear Hill model describing enzymatic reaction kinetics and derive a class of symmetry transformations for each order of the model. We consider a minimal example consisting of the application of symmetry-based methods to a model selection problem, where we are able to demonstrate superior performance compared to ordinary residual-based model selection. Moreover, we demonstrate that symmetries reveal the intrinsic properties of a system of interest based on a single time series. Finally, we show and propose that symmetry-based methodology should be considered as the first step in a systematic model building and in the case when multiple time series are available it should complement the commonly used statistical methodologies.

Synergistic effects of repair, resilience and retention of damage determine the conditions for replicative ageing.

Accumulation of damaged proteins is a hallmark of ageing, occurring in organisms ranging from bacteria and yeast to mammalian cells. During cell division in Saccharomyces cerevisiae, damaged proteins are retained within the mother cell, resulting in an ageing mother while a new daughter cell exhibits full replicative potential. The cell-specific features determining the ageing remain elusive. It has been suggested that the replicative ageing is dependent on the ability of the cell to repair and retain pre-existing damage. To deepen the understanding of how these factors influence the life of individual cells, we developed and experimentally validated a dynamic model of damage accumulation accounting for replicative ageing on the single cell level. The model includes five essential properties: cell growth, damage formation, damage repair, cell division and cell death, represented in a theoretical framework describing the conditions allowing for replicative ageing, starvation, immortality or clonal senescence. We introduce the resilience to damage, which can be interpreted as the difference in volume between an old and a young cell. We show that the capacity to retain damage deteriorates with high age, that asymmetric division allows for retention of damage, and that there is a trade-off between retention and the resilience property. Finally, we derive the maximal degree of asymmetry as a function of resilience, proposing that asymmetric cell division is beneficial with respect to replicative ageing as it increases the lifespan of a given organism. The proposed model contributes to a deeper understanding of the ageing process in eukaryotic organisms.

Single-cell study links metabolism with nutrient signaling and reveals sources of variability.

BACKGROUND: The yeast AMPK/SNF1 pathway is best known for its role in glucose de/repression. When glucose becomes limited, the Snf1 kinase is activated and phosphorylates the transcriptional repressor Mig1, which is then exported from the nucleus. The exact mechanism how the Snf1-Mig1 pathway is regulated is not entirely elucidated. RESULTS: Glucose uptake through the low affinity transporter Hxt1 results in nuclear accumulation of Mig1 in response to all glucose concentrations upshift, however with increasing glucose concentration the nuclear localization of Mig1 is more intense. Strains expressing Hxt7 display a constant response to all glucose concentration upshifts. We show that differences in amount of hexose transporter molecules in the cell could cause cell-to-cell variability in the Mig1-Snf1 system. We further apply mathematical modelling to our data, both general deterministic and a nonlinear mixed effect model. Our model suggests a presently unrecognized regulatory step of the Snf1-Mig1 pathway at the level of Mig1 dephosphorylation. Model predictions point to parameters involved in the transport of Mig1 in and out of the nucleus as a majorsource of cell to cell variability. CONCLUSIONS: With this modelling approach we have been able to suggest steps that contribute to the cell-to-cell variability. Our data indicate a close link between the glucose uptake rate, which determines the glycolytic rate, and the activity of the Snf1/Mig1 system. This study hence establishes a close relation between metabolism and signalling.