Inference of the High-Level Interaction Topology between the Metabolic and Cell-Cycle Oscillators from Single-Cell Dynamics

Newton's cradle: Cell cycle regulation by two mutually inhibitory oscillators

The cell division cycle is a fundamental physiological process displaying a great degree of plasticity during the course of multicellular development. This plasticity is evident in the transition from rapid and stringently-timed divisions of the early embryo to subsequent size-controlled mitotic cycles. Later in development, cells may pause and restart proliferation in response to myriads of internal or external signals, or permanently exit the cell cycle following terminal differentiation or senescence. Beyond this, cells can undergo modified cell division variants, such as endoreplication, which increases their ploidy, or meiosis, which reduces their ploidy. This wealth of behaviours has led to numerous conceptual analogies intended as frameworks for understanding the proliferative program. Here, we aim to unify these mechanisms under one dynamical paradigm. To this end, we take a control theoretical approach to frame the cell cycle as a pair of arrestable and mutually-inhibiting, doubly amplified, negative feedback oscillators controlling chromosome replication and segregation events, respectively. Under appropriate conditions, this framework can reproduce fixed-period oscillations, checkpoint arrests of variable duration, and endocycles. Subsequently, we use phase plane and bifurcation analysis to explain the dynamical basis of these properties. Then, using a physiologically realistic, biochemical model, we show that the very same regulatory structure underpins the diverse functions of the cell cycle control network. We conclude that Newton's cradle may be a suitable mechanical analogy of how the cell cycle is regulated.

Improving pathway prediction accuracy of constraints-based metabolic network models by treating enzymes as microcompartments

Metabolic network models have become increasingly precise and accurate as the most widespread and practical digital representations of living cells. The prediction functions were significantly expanded by integrating cellular resources and abiotic constraints in recent years. However, if unreasonable modeling methods were adopted due to a lack of consideration of biological knowledge, the conflicts between stoichiometric and other constraints, such as thermodynamic feasibility and enzyme resource availability, would lead to distorted predictions. In this work, we investigated a prediction anomaly of EcoETM, a constraints-based metabolic network model, and introduced the idea of enzyme compartmentalization into the analysis process. Through rational combination of reactions, we avoid the false prediction of pathway feasibility caused by the unrealistic assumption of free intermediate metabolites. This allowed us to correct the pathway structures of l-serine and l-tryptophan. A specific analysis explains the application method of the EcoETM-like model and demonstrates its potential and value in correcting the prediction results in pathway structure by resolving the conflict between different constraints and incorporating the evolved roles of enzymes as reaction compartments. Notably, this work also reveals the trade-off between product yield and thermodynamic feasibility. Our work is of great value for the structural improvement of constraints-based models.

Global quantitative understanding of non-equilibrium cell fate decision-making in response to pheromone

Cell-cycle arrest and polarized growth are commonly used to characterize the response of yeast to pheromone. However, the quantitative decision-making processes underlying time-dependent changes in cell fate remain unclear. In this study, we conducted single-cell level experiments to observe multidimensional responses, uncovering diverse fates of yeast cells. Multiple states are revealed, along with the kinetic switching rates and pathways among them, giving rise to a quantitative landscape of mating response. To quantify the experimentally observed cell fates, we developed a theoretical framework based on non-equilibrium landscape and flux theory. Additionally, we performed stochastic simulations of biochemical reactions to elucidate signal transduction and cell growth. Notably, our experimental findings have provided the first global quantitative evidence of the real-time synchronization between intracellular signaling, physiological growth, and morphological functions. These results validate the proposed underlying mechanism governing the emergence of multiple cell fate states. This study introduces an emerging mechanistic approach to understand non-equilibrium cell fate decision-making in response to pheromone.

Cytoplasmic division cycles without the nucleus and mitotic CDK/cyclin complexes

Cytoplasmic divisions are thought to rely on nuclear divisions and mitotic signals. We demonstrate in Drosophila embryos that cytoplasm can divide repeatedly without nuclei and mitotic CDK/cyclin complexes. Cdk1 normally slows an otherwise faster cytoplasmic division cycle, coupling it with nuclear divisions, and when uncoupled, cytoplasm starts dividing before mitosis. In developing embryos where CDK/cyclin activity can license mitotic microtubule (MT) organizers like the spindle, cytoplasmic divisions can occur without the centrosome, a principal organizer of interphase MTs. However, centrosomes become essential in the absence of CDK/cyclin activity, implying that the cytoplasm can employ either the centrosome-based interphase or CDK/cyclin-dependent mitotic MTs to facilitate its divisions. Finally, we present evidence that autonomous cytoplasmic divisions occur during unperturbed fly embryogenesis and that they may help extrude mitotically stalled nuclei during blastoderm formation. We postulate that cytoplasmic divisions occur in cycles governed by a yet-to-be-uncovered clock mechanism autonomous from CDK/cyclin complexes.

Network regulation meets substrate modification chemistry

Biochemical networks are at the heart of cellular information processing. These networks contain distinct facets: (i) processing of information from the environment via cascades/pathways along with network regulation and (ii) modification of substrates in different ways, to confer protein functionality, stability and processing. While many studies focus on these factors individually, how they interact and the consequences for cellular systems behaviour are poorly understood. We develop a systems framework for this purpose by examining the interplay of network regulation (canonical feedback and feed-forward circuits) and multisite modification, as an exemplar of substrate modification. Using computational, analytical and semi-analytical approaches, we reveal distinct and unexpected ways in which the substrate modification and network levels combine and the emergent behaviour arising therefrom. This has important consequences for dissecting the behaviour of specific signalling networks, tracing the origins of systems behaviour, inference of networks from data, robustness/evolvability and multi-level engineering of biomolecular networks. Overall, we repeatedly demonstrate how focusing on only one level (say network regulation) can lead to profoundly misleading conclusions about all these aspects, and reveal a number of important consequences for experimental/theoretical/data-driven interrogations of cellular signalling systems.

Temporal segregation of biosynthetic processes is responsible for metabolic oscillations during the budding yeast cell cycle

Many cell biological and biochemical mechanisms controlling the fundamental process of eukaryotic cell division have been identified; however, the temporal dynamics of biosynthetic processes during the cell division cycle are still elusive. Here, we show that key biosynthetic processes are temporally segregated along the cell cycle. Using budding yeast as a model and single-cell methods to dynamically measure metabolic activity, we observe two peaks in protein synthesis, in the G1 and S/G2/M phase, whereas lipid and polysaccharide synthesis peaks only once, during the S/G2/M phase. Integrating the inferred biosynthetic rates into a thermodynamic-stoichiometric metabolic model, we find that this temporal segregation in biosynthetic processes causes flux changes in primary metabolism, with an acceleration of glucose-uptake flux in G1 and phase-shifted oscillations of oxygen and carbon dioxide exchanges. Through experimental validation of the model predictions, we demonstrate that primary metabolism oscillates with cell-cycle periodicity to satisfy the changing demands of biosynthetic processes exhibiting unexpected dynamics during the cell cycle.

Spatial-temporal regulation of fatty alcohol biosynthesis in yeast

BACKGROUND

Construction of efficient microbial cell factories is one of the core steps for establishing green bio-manufacturing processes. However, the complex metabolic regulation makes it challenging in driving the metabolic flux toward the product biosynthesis. Dynamically coupling the biosynthetic pathways with the cellular metabolism at spatial-temporal manner should be helpful for improving the production with alleviating the cellular stresses.

RESULTS

In this study, we observed the mismatch between fatty alcohol biosynthesis and cellular metabolism, which compromised the fatty alcohol production in Saccharomyces cerevisiae. To enhance the fatty alcohol production, we spatial-temporally regulated fatty alcohol biosynthetic pathway by peroxisomal compartmentalization (spatial) and dynamic regulation of gene expression (temporal). In particular, fatty acid/acyl-CoA responsive promoters were identified by comparative transcriptional analysis, which helped to dynamically regulate the expression of acyl-CoA reductase gene MaFAR1 and improved fatty alcohol biosynthesis by 1.62-fold. Furthermore, enhancing the peroxisomal supply of acyl-CoA and NADPH further improved fatty alcohol production to 282 mg/L, 2.52 times higher than the starting strain.

CONCLUSIONS

This spatial-temporal regulation strategy partially coordinated fatty alcohol biosynthesis with cellular metabolism including peroxisome biogenesis and precursor supply, which should be applied for production of other products in microbes.

Genome-Wide Analysis of Yeast Metabolic Cycle through Metabolic Network Models Reveals Superiority of Integrated ATAC-seq Data over RNA-seq Data

Saccharomyces cerevisiae undergoes robust oscillations to regulate its physiology for adaptation and survival under nutrient-limited conditions. Environmental cues can induce rhythmic metabolic alterations in order to facilitate the coordination of dynamic metabolic behaviors. Of such metabolic processes, the yeast metabolic cycle enables adaptation of the cells to varying nutritional status through oscillations in gene expression and metabolite production levels. In this process, yeast metabolism is altered between diverse cellular states based on changing oxygen consumption levels: quiescent (reductive charging [RC]), growth (oxidative [OX]), and proliferation (reductive building [RB]) phases. We characterized metabolic alterations during the yeast metabolic cycle using a variety of approaches. Gene expression levels are widely used for condition-specific metabolic simulations, whereas the use of epigenetic information in metabolic modeling is still limited despite the clear relationship between epigenetics and metabolism. This prompted us to investigate the contribution of epigenomic information to metabolic predictions for progression of the yeast metabolic cycle. In this regard, we determined altered pathways through the prediction of regulated reactions and corresponding model genes relying on differential chromatin accessibility levels. The predicted metabolic alterations were confirmed via data analysis and literature. We subsequently utilized RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) data sets in the contextualization of the yeast model. The use of ATAC-seq data considerably enhanced the predictive capability of the model. To the best of our knowledge, this is the first attempt to use genome-wide chromatin accessibility data in metabolic modeling. The preliminary results showed that epigenomic data sets can pave the way for more accurate metabolic simulations.

IMPORTANCE Dynamic chromatin organization mediates the emergence of condition-specific phenotypes in eukaryotic organisms. Saccharomyces cerevisiae can alter its metabolic profile via regulation of genome accessibility and robust transcriptional oscillations under nutrient-limited conditions. Thus, both epigenetic information and transcriptomic information are crucial in the understanding of condition-specific metabolic behavior in this organism. Based on genome-wide alterations in chromatin accessibility and transcription, we investigated the yeast metabolic cycle, which is a remarkable example of coordinated and dynamic yeast behavior. In this regard, we assessed the use of ATAC-seq and RNA-seq data sets in condition-specific metabolic modeling. To our knowledge, this is the first attempt to use chromatin accessibility data in the reconstruction of context-specific metabolic models, despite the extensive use of transcriptomic data. As a result of comparative analyses, we propose that the incorporation of epigenetic information is a promising approach in the accurate prediction of metabolic dynamics.

A frequency-amplitude coordinator and its optimal energy consumption for biological oscillators

Biorhythm including neuron firing and protein-mRNA interaction are fundamental activities with diffusive effect. Their well-balanced spatiotemporal dynamics are beneficial for healthy sustainability. Therefore, calibrating both anomalous frequency and amplitude of biorhythm prevents physiological dysfunctions or diseases. However, many works were devoted to modulate frequency exclusively whereas amplitude is usually ignored, although both quantities are equally significant for coordinating biological functions and outputs. Especially, a feasible method coordinating the two quantities concurrently and precisely is still lacking. Here, for the first time, we propose a universal approach to design a frequency-amplitude coordinator rigorously via dynamical systems tools. We consider both spatial and temporal information. With a single well-designed coordinator, they can be calibrated to desired levels simultaneously and precisely. The practical usefulness and efficacy of our method are demonstrated in representative neuronal and gene regulatory models. We further reveal its fundamental mechanism and optimal energy consumption providing inspiration for biorhythm regulation in future.

Stress-induced growth rate reduction restricts metabolic resource utilization to modulate osmo-adaptation time

A near-constant feature of stress responses is a downregulation or arrest of the cell cycle, resulting in transient growth slowdown. To investigate the role of growth slowdown in the hyperosmotic shock response of S. cerevisiae, we perturbed the G1/S checkpoint protein Sic1 to enable osmo-stress response activation with diminished growth slowdown. We document that in this mutant, adaptation to stress is accelerated rather than delayed. This accelerated recovery of the mutant proceeds by liquidation of internal glycogen stores, which are then shunted into the osmo-shock response. Therefore, osmo-adaptation in wild-type cells is delayed because growth slowdown prevents full accessibility to cellular glycogen stores. However, faster adaptation comes at the cost of acute sensitivity to subsequent osmo-stresses. We suggest that stress-induced growth slowdown acts as an arbiter to regulate the resources devoted to osmo-shock, balancing short-term adaptation with long-term robustness.

Multiple Layers of Phospho-Regulation Coordinate Metabolism and the Cell Cycle in Budding Yeast

The coordination of metabolism and growth with cell division is crucial for proliferation. While it has long been known that cell metabolism regulates the cell division cycle, it is becoming increasingly clear that the cell division cycle also regulates metabolism. In budding yeast, we previously showed that over half of all measured metabolites change concentration through the cell cycle indicating that metabolic fluxes are extensively regulated during cell cycle progression. However, how this regulation is achieved still remains poorly understood. Since both the cell cycle and metabolism are regulated to a large extent by protein phosphorylation, we here decided to measure the phosphoproteome through the budding yeast cell cycle. Specifically, we chose a cell cycle synchronization strategy that avoids stress and nutrient-related perturbations of metabolism, and we grew the yeast on ethanol minimal medium to force cells to utilize their full biosynthetic repertoire. Using a tandem-mass-tagging approach, we found over 200 sites on metabolic enzymes and transporters to be phospho-regulated. These sites were distributed among many pathways including carbohydrate catabolism, lipid metabolism, and amino acid synthesis and therefore likely contribute to changing metabolic fluxes through the cell cycle. Among all one thousand sites whose phosphorylation increases through the cell cycle, the CDK consensus motif and an arginine-directed motif were highly enriched. This arginine-directed R-R-x-S motif is associated with protein-kinase A, which regulates metabolism and promotes growth. Finally, we also found over one thousand sites that are dephosphorylated through the G1/S transition. We speculate that the phosphatase Glc7/PP1, known to regulate both the cell cycle and carbon metabolism, may play an important role because its regulatory subunits are phospho-regulated in our data. In summary, our results identify extensive cell cycle dependent phosphorylation and dephosphorylation of metabolic enzymes and suggest multiple mechanisms through which the cell division cycle regulates metabolic signaling pathways to temporally coordinate biosynthesis with distinct phases of the cell division cycle.

Differential scaling between G1 protein production and cell size dynamics promotes commitment to the cell division cycle in budding yeast

In the unicellular eukaryote Saccharomyces cerevisiae, Cln3-cyclin-dependent kinase activity enables Start, the irreversible commitment to the cell division cycle. However, the concentration of Cln3 has been paradoxically considered to remain constant during G1, due to the presumed scaling of its production rate with cell size dynamics. Measuring metabolic and biosynthetic activity during cell cycle progression in single cells, we found that cells exhibit pulses in their protein production rate. Rather than scaling with cell size dynamics, these pulses follow the intrinsic metabolic dynamics, peaking around Start. Using a viral-based bicistronic construct and targeted proteomics to measure Cln3 at the single-cell and population levels, we show that the differential scaling between protein production and cell size leads to a temporal increase in Cln3 concentration, and passage through Start. This differential scaling causes Start in both daughter and mother cells across growth conditions. Thus, uncoupling between two fundamental physiological parameters drives cell cycle commitment.