Stochastic modeling of aging cells reveals how damage accumulation, repair, and cell-division asymmetry affect clonal senescence and population fitness

Balance between asymmetric allocation and repair of somatic damage in unicellular life forms as an ancient form of r/K selection

Over the course of multiple divisions, cells accumulate diverse nongenetic, somatic damage including misfolded and aggregated proteins and cell wall defects. If the rate of damage accumulation exceeds the rate of dilution through cell growth, a dedicated mitigation strategy is required to prevent eventual population collapse. Strategies for somatic damage control can be divided into two categories, asymmetric allocation and repair, which are not, in principle, mutually exclusive. We explore a mathematical model to identify the optimal strategy, maximizing the total cell number, over a wide range of environmental and physiological conditions. The optimal strategy is primarily determined by extrinsic, damage-independent mortality and the physiological model for damage accumulation that can be either independent (linear) or increasing (exponential) with respect to the prior accumulated damage. Under the linear regime, the optimal strategy is either exclusively repair or asymmetric allocation, whereas under the exponential regime, the optimal strategy is a combination of asymmetry and repair. Repair is preferred when extrinsic mortality is low, whereas at high extrinsic mortality, asymmetric damage allocation becomes the strategy of choice. We hypothesize that at an early stage of life evolution, optimization over repair and asymmetric allocation of somatic damage gave rise to r and K selection strategists.

Sis2 regulates yeast replicative lifespan in a dose-dependent manner

Application of microfluidic platforms facilitated high-precision measurements of yeast replicative lifespan (RLS); however, comparative quantification of lifespan across strain libraries has been missing. Here we microfluidically measure the RLS of 307 yeast strains, each deleted for a single gene. Despite previous reports of extended lifespan in these strains, we found that 56% of them did not actually live longer than the wild-type; while the remaining 44% showed extended lifespans, the degree of extension was often different from what was previously reported. Deletion of SIS2 gene led to the largest RLS increase observed. Sis2 regulated yeast lifespan in a dose-dependent manner, implying a role for the coenzyme A biosynthesis pathway in lifespan regulation. Introduction of the human PPCDC gene in the sis2Δ background neutralized the lifespan extension. RNA-seq experiments revealed transcriptional increases in cell-cycle machinery components in sis2Δ background. High-precision lifespan measurement will be essential to elucidate the gene network governing lifespan.

Unicellular life balances asymmetric allocation and repair of somatic damage representing the origin of r/K selection

Over the course of multiple divisions, cells accumulate diverse non-genetic, somatic damage including misfolded and aggregated proteins and cell wall defects. If the rate of damage accumulation exceeds the rate of dilution through cell growth, a dedicated mitigation strategy is required to prevent eventual population collapse. Strategies for somatic damage control can be divided into two categories, asymmetric allocation and repair, which are not, in principle, mutually exclusive. Through mathematical modelling, we identify the optimal strategy, maximizing the total cell number, over a wide range of environmental and physiological conditions. The optimal strategy is primarily determined by extrinsic (damage-independent) mortality and the physiological model for damage accumulation that can be either independent (linear) or increasing (exponential) with respect to the prior accumulated damage. Under the linear regime, the optimal strategy is either exclusively repair or asymmetric allocation whereas under the exponential regime, the optimal strategy is mixed. Repair is preferred when extrinsic mortality is low, whereas at high extrinsic mortality, asymmetric damage allocation becomes the strategy of choice. We hypothesize that optimization over somatic damage repair and asymmetric allocation in early cellular life forms gave rise to the r and K selection strategies.

The choice of the objective function in flux balance analysis is crucial for predicting replicative lifespans in yeast

Flux balance analysis (FBA) is a powerful tool to study genome-scale models of the cellular metabolism, based on finding the optimal flux distributions over the network. While the objective function is crucial for the outcome, its choice, even though motivated by evolutionary arguments, has not been directly connected to related measures. Here, we used an available multi-scale mathematical model of yeast replicative ageing, integrating cellular metabolism, nutrient sensing and damage accumulation, to systematically test the effect of commonly used objective functions on features of replicative ageing in budding yeast, such as the number of cell divisions and the corresponding time between divisions. The simulations confirmed that assuming maximal growth is essential for reaching realistic lifespans. The usage of the parsimonious solution or the additional maximisation of a growth-independent energy cost can improve lifespan predictions, explained by either increased respiratory activity using resources otherwise allocated to cellular growth or by enhancing antioxidative activity, specifically in early life. Our work provides a new perspective on choosing the objective function in FBA by connecting it to replicative ageing.

Changes in DNA double-strand break repair during aging correlate with an increase in genomic mutations

A double -strand break (DSB) is one of the most deleterious forms of DNA damage. In eukaryotic cells, two main repair pathways have evolved to repair DSBs, homologous recombination (HR) and non-homologous end-joining (NHEJ). HR is the predominant pathway of repair in the unicellular eukaryotic organism, S. cerevisiae. However, during replicative aging the relative use of HR and NHEJ shifts in favor of end-joining repair. By monitoring repair events in the HO-DSB system, we find that early in replicative aging there is a decrease in the association of long-range resection factors, Dna2-Sgs1 and Exo1 at the break site and a decrease in DNA resection. Subsequently, as aging progressed, the recovery of Ku70 at DSBs decreased and the break site associated with the nuclear pore complex at the nuclear periphery, which is the location where DSB repair occurs through alternative pathways that are more mutagenic. End-bridging remained intact as HR and NHEJ declined, but eventually it too became disrupted in cells at advanced replicative age. In all, our work provides insight into the molecular changes in DSB repair pathway during replicative aging. HR first declined, resulting in a transient increase in the NHEJ. However, with increased cellular divisions, Ku70 recovery at DSBs and NHEJ subsequently declined. In wild type cells of advanced replicative age, there was a high frequency of repair products with genomic deletions and microhomologies at the break junction, events not observed in young cells which repaired primarily by HR.

Multi-scale model suggests the trade-off between protein and ATP demand as a driver of metabolic changes during yeast replicative ageing

The accumulation of protein damage is one of the major drivers of replicative ageing, describing a cell’s reduced ability to reproduce over time even under optimal conditions. Reactive oxygen and nitrogen species are precursors of protein damage and therefore tightly linked to ageing. At the same time, they are an inevitable by-product of the cell’s metabolism. Cells are able to sense high levels of reactive oxygen and nitrogen species and can subsequently adapt their metabolism through gene regulation to slow down damage accumulation. However, the older or damaged a cell is the less flexibility it has to allocate enzymes across the metabolic network, forcing further adaptions in the metabolism. To investigate changes in the metabolism during replicative ageing, we developed an multi-scale mathematical model using budding yeast as a model organism. The model consists of three interconnected modules: a Boolean model of the signalling network, an enzyme-constrained flux balance model of the central carbon metabolism and a dynamic model of growth and protein damage accumulation with discrete cell divisions. The model can explain known features of replicative ageing, like average lifespan and increase in generation time during successive division, in yeast wildtype cells by a decreasing pool of functional enzymes and an increasing energy demand for maintenance. We further used the model to identify three consecutive metabolic phases, that a cell can undergo during its life, and their influence on the replicative potential, and proposed an intervention span for lifespan control.

The effect of stress on biophysical characteristics of misfolded protein aggregates in living Saccharomyces cerevisiae cells

Aggregation of misfolded or damaged proteins is often attributed to numerous metabolic and neurodegenerative disorders. To reveal underlying mechanisms and cellular responses, it is crucial to investigate protein aggregate dynamics in cells. Here, we used super-resolution single-molecule microscopy to obtain biophysical characteristics of individual aggregates of a model misfolded protein ∆ssCPY* labelled with GFP. We demonstrated that oxidative and hyperosmotic stress lead to increased aggregate stoichiometries but not necessarily the total number of aggregates. Moreover, our data suggest the importance of the thioredoxin peroxidase Tsa1 for the controlled sequestering and clearance of aggregates upon both conditions. Our work provides novel insights into the understanding of the cellular response to stress via revealing the dynamical properties of stress-induced protein aggregates.

Modeling aging and its impact on cellular function and organismal behavior

Aging is a complex phenomenon of functional decay in a biological organism. Although the effects of aging are readily recognizable in a wide range of organisms, the cause(s) of aging are ill defined and poorly understood. Experimental methods on model organisms have driven significant insight into aging as a process, but have not provided a complete model of aging. Computational biology offers a unique opportunity to resolve this gap in our knowledge by generating extensive and testable models that can help us understand the fundamental nature of aging, identify the presence and characteristics of unaccounted aging factor(s), demonstrate the mechanics of particular factor(s) in driving aging, and understand the secondary effects of aging on biological function. In this review, we will address each of the above roles for computational biology in aging research. Concurrently, we will explore the different applications of computational biology to aging in single-celled versus multicellular organisms. Given the long history of computational biogerontological research on lower eukaryotes, we emphasize the key future goals of gradually integrating prior models into a holistic map of aging and translating successful models to higher-complexity organisms.

Senolytics and the compression of late-life mortality

Senescent cells play an important role in mammalian ageing and in the etiology of age-related diseases. Treatment of mice with senolytics - drugs that selectively remove senescent cells - causes an extension of median lifespan but has little effect on maximum lifespan. Postponement of some mortality to later ages, without a corresponding increase in maximum mortality, can be termed 'compression of mortality'. When we fit the standard Gompertz mortality model to the survival data following senolytic treatment, we find an increase in the slope parameter, commonly described as the 'actuarial ageing rate'. These observations raise important questions about the actions of senolytic treatments and their effects on health and survival, which are not yet sufficiently understood. To explore how the survival data from senolytics experiments might be explained, we combine a recent exploration of the evolutionary basis of cellular senescence with theoretical consideration of the molecular processes that might be involved. We perform numerical simulations of senescent cell accumulation and senolytic treatment in an ageing population. The simulations suggest that while senolytics diminish the burden of senescent cells, they may also impair the general repair capacity of the organism, leading to a faster accumulation post-treatment of new senescent cells. Our results suggest a framework to address the benefits and possible side effects of senolytic therapies, with the potential to aid in the design of optimal treatment regimens.

Damage Repair versus Aging in an Individual-Based Model of Biofilms

Damage is an inevitable consequence of life. For unicellular organisms, this leads to a trade-off between allocating resources into damage repair or into growth coupled with segregation of damage upon cell division, i.e., aging and senescence. Few studies considered repair as an alternative to senescence. None considered biofilms, where the majority of unicellular organisms live, although fitness advantages in well-mixed systems often turn into disadvantages in spatially structured systems such as biofilms. We compared the fitness consequences of aging versus an adaptive repair mechanism based on sensing damage, using an individual-based model of a generic unicellular organism growing in biofilms. We found that senescence is not beneficial provided that growth is limited by substrate availability. Instead, it is useful as a stress response to deal with damage that failed to be repaired when (i) extrinsic mortality was high; (ii) a degree of multicellularity was present; and (iii) damage segregation was effective.

The extent of senescence due to damage accumulation—or aging—is evidently evolvable as it differs hugely between species and is not universal, suggesting that its fitness advantages depend on life history and environment. In contrast, repair of damage is present in all organisms studied. Despite the fundamental trade-off between investing resources into repair or into growth, repair and segregation of damage have not always been considered alternatives. For unicellular organisms, unrepaired damage could be divided asymmetrically between daughter cells, leading to senescence of one and rejuvenation of the other. Repair of “unicells” has been predicted to be advantageous in well-mixed environments such as chemostats. Most microorganisms, however, live in spatially structured systems, such as biofilms, with gradients of environmental conditions and cellular physiology as well as a clonal population structure. To investigate whether this clonal structure might favor senescence by damage segregation (a division-of-labor strategy akin to the germline-soma division in multicellular organisms), we used an individual-based computational model and developed an adaptive repair strategy where cells respond to their current intracellular damage levels by investing into repair machinery accordingly. Our simulations showed that the new adaptive repair strategy was advantageous provided that growth was limited by substrate availability, which is typical for biofilms. Thus, biofilms do not favor a germline-soma-like division of labor between daughter cells in terms of damage segregation. We suggest that damage segregation is beneficial only when extrinsic mortality is high, a degree of multicellularity is present, and an active mechanism makes segregation effective.

IMPORTANCE Damage is an inevitable consequence of life. For unicellular organisms, this leads to a trade-off between allocating resources into damage repair or into growth coupled with segregation of damage upon cell division, i.e., aging and senescence. Few studies considered repair as an alternative to senescence. None considered biofilms, where the majority of unicellular organisms live, although fitness advantages in well-mixed systems often turn into disadvantages in spatially structured systems such as biofilms. We compared the fitness consequences of aging versus an adaptive repair mechanism based on sensing damage, using an individual-based model of a generic unicellular organism growing in biofilms. We found that senescence is not beneficial provided that growth is limited by substrate availability. Instead, it is useful as a stress response to deal with damage that failed to be repaired when (i) extrinsic mortality was high; (ii) a degree of multicellularity was present; and (iii) damage segregation was effective.

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.

Astaxanthin as a Putative Geroprotector: Molecular Basis and Focus on Brain Aging

In recent years, the scientific interest in natural compounds with geroprotective activities has grown exponentially. Among the various naturally derived molecules, astaxanthin (ASX) represents a highly promising candidate geroprotector. By virtue of the central polyene chain, ASX acts as a scavenger of free radicals in the internal membrane layer and simultaneously controls oxidation on the membrane surface. Moreover, several studies have highlighted ASX’s ability to modulate numerous biological mechanisms at the cellular level, including the modulation of transcription factors and genes directly linked to longevity-related pathways. One of the main relevant evolutionarily-conserved transcription factors modulated by astaxanthin is the forkhead box O3 gene (FOXO3), which has been recognized as a critical controller of cell fate and function. Moreover, FOXO3 is one of only two genes shown to robustly affect human longevity. Due to its tropism in the brain, ASX has recently been studied as a putative neuroprotective molecule capable of delaying or preventing brain aging in different experimental models of brain damage or neurodegenerative diseases. Astaxanthin has been observed to slow down brain aging by increasing brain-derived neurotrophic factor (BDNF) levels in the brain, attenuating oxidative damage to lipids, protein, and DNA and protecting mitochondrial functions. Emerging data now suggest that ASX can modulate Nrf2, FOXO3, Sirt1, and Klotho proteins that are linked to longevity. Together, these mechanisms provide support for a role of ASX as a potential geroneuroprotector.

Optimum Threshold Minimizes Noise in Timing of Intracellular Events

How the noisy expression of regulatory proteins affects timing of intracellular events is an intriguing fundamental problem that influences diverse cellular processes. Here we use the bacteriophage λ to study event timing in individual cells where cell lysis is the result of expression and accumulation of a single protein (holin) in the Escherichia coli cell membrane up to a critical threshold level. Site-directed mutagenesis of the holin gene generated phage variants that vary in their lysis times from 30 to 190 min. Observation of the lysis times of single cells reveals an intriguing finding-the noise in lysis timing first decreases with increasing lysis time to reach a minimum and then sharply increases at longer lysis times. A mathematical model with stochastic expression of holin together with dilution from cell growth was sufficient to explain the non-monotonic noise profile and identify holin accumulation thresholds that generate precision in lysis timing.

Characterization of the impact of GMP/GDP synthesis inhibition on replicative lifespan extension in yeast

Slowing down aging-associated accumulation of molecular damage or its prevention represents a promising therapeutic paradigm to combat aging-related disease and death. While several chemical compounds extend lifespan in model organisms, their mechanism of action is often unknown, reducing their therapeutic potential. Using a systematic approach, here we characterize the impact of the GMP pathway on yeast lifespan and elucidate GMP synthesis inhibition as a lifespan extension mechanism. We further discover that proteasome activation extends lifespan in part through the GMP pathway. GMP synthesis inhibition exerts its lifespan extension effect independently of the canonical nutrient-sensing pathway regulating lifespan. Exposing longitudinally aging yeast cells to GMP pathway inhibition in an age-dependent manner, we demonstrate that the lifespan extension is facilitated by slowing, rather than reversing, the aging process in cells. Using a GUK1 mutant with lower GMP-to-GDP conversion activity, we observe lifespan extension, suggesting that reduced GDP level by itself can also extend yeast lifespan. These findings elucidate the involvement of nucleotide metabolism in the aging process. The existence of clinically-approved GMP pathway inhibitors elicits the potential of a new class of therapeutics for aging-related disorders.