Network reconstruction and validation of the Snf1/AMPK pathway in baker's yeast based on a comprehensive literature review

AMPK Activates Cellulase Secretion in Penicillium funiculosum by Downregulating P-HOG1 MAPK Levels

Cellulase production for hydrolyzing plant cell walls is energy-intensive in filamentous fungi during nutrient scarcity. AMP-activated protein kinase (AMPK), encoded by snf1, is known to be the nutrient and energy sensor in eukaryotes. Previous studies on AMPK identified its role in alternate carbon utilization in pathogenic fungi. However, the precise role of AMPK in cellulase production remains elusive. In the present study, we employed gene-deletion analysis, quantitative proteomics and chemical-genetic approaches to investigate the role of AMPK in cellulase synthesis in Penicillium funiculosum. Gene-deletion analysis revealed that AMPK does not promote transcription and translation but is essential for cellulase secretion in a calcium-dependent manner. Proteomic analysis of the snf1-deleted (Δsnf1) strain confirmed trapped cellulase inside the mycelia and identified HOG1 MAPK activation as the most significant Ca2+-induced signaling event during carbon stress in Δsnf1. Western blot analysis analysis revealed that the phosphorylated HOG1 (P-HOG1)/HOG1 MAPK ratio maintained by Ca2+-signaling/Ca2+-activated AMPK, respectively, forms a secretion checkpoint for cellulases, and disturbing this equilibrium blocks cellulase secretion. The proteomic analysis also indicated a massive increase in mTORC1-activated anabolic pathways during carbon stress in Δsnf1. Our study suggests that AMPK maintains homeostasis by acting as a global repressor during carbon stress.

Data integration strategies for whole-cell modeling

Data makes the world go round-and high quality data is a prerequisite for precise models, especially for whole-cell models (WCM). Data for WCM must be reusable, contain information about the exact experimental background, and should-in its entirety-cover all relevant processes in the cell. Here, we review basic requirements to data for WCM and strategies how to combine them. As a species-specific resource, we introduce the Yeast Cell Model Data Base (YCMDB) to illustrate requirements and solutions. We discuss recent standards for data as well as for computational models including the modeling process as data to be reported. We outline strategies for constructions of WCM despite their inherent complexity.

A detailed molecular network map and model of the NLRP3 inflammasome

The NLRP3 inflammasome is a key regulator of inflammation that responds to a broad range of stimuli. The exact mechanism of activation has not been determined, but there is a consensus on cellular potassium efflux as a major common denominator. Once NLRP3 is activated, it forms high-order complexes together with NEK7 that trigger aggregation of ASC into specks. Typically, there is only one speck per cell, consistent with the proposal that specks form - or end up at - the centrosome. ASC polymerisation in turn triggers caspase-1 activation, leading to maturation and release of IL-1β and pyroptosis, i.e., highly inflammatory cell death. Several gain-of-function mutations in the NLRP3 inflammasome have been suggested to induce spontaneous activation of NLRP3 and hence contribute to development and disease severity in numerous autoinflammatory and autoimmune diseases. Consequently, the NLRP3 inflammasome is of significant clinical interest, and recent attention has drastically improved our insight in the range of involved triggers and mechanisms of signal transduction. However, despite recent progress in knowledge, a clear and comprehensive overview of how these mechanisms interplay to shape the system level function is missing from the literature. Here, we provide such an overview as a resource to researchers working in or entering the field, as well as a computational model that allows for evaluating and explaining the function of the NLRP3 inflammasome system from the current molecular knowledge. We present a detailed reconstruction of the molecular network surrounding the NLRP3 inflammasome, which account for each specific reaction and the known regulatory constraints on each event as well as the mechanisms of drug action and impact of genetics when known. Furthermore, an executable model from this network reconstruction is generated with the aim to be used to explain NLRP3 activation from priming and activation to the maturation and release of IL-1β and IL-18. Finally, we test this detailed mechanistic model against data on the effect of different modes of inhibition of NLRP3 assembly. While the exact mechanisms of NLRP3 activation remains elusive, the literature indicates that the different stimuli converge on a single activation mechanism that is additionally controlled by distinct (positive or negative) priming and licensing events through covalent modifications of the NLRP3 molecule. Taken together, we present a compilation of the literature knowledge on the molecular mechanisms on NLRP3 activation, a detailed mechanistic model of NLRP3 activation, and explore the convergence of diverse NLRP3 activation stimuli into a single input mechanism.

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.

Modelling of glucose repression signalling in yeast Saccharomyces cerevisiae

Saccharomyces cerevisiae has a sophisticated signalling system that plays a crucial role in cellular adaptation to changing environments. The SNF1 pathway regulates energy homeostasis upon glucose derepression; hence, it plays an important role in various processes, such as metabolism, cell cycle and autophagy. To unravel its behaviour, SNF1 signalling has been extensively studied. However, the pathway components are strongly interconnected and inconstant; therefore, elucidating its dynamic behaviour based on experimental data only is challenging. To tackle this complexity, systems biology approaches have been successfully employed. This review summarizes the progress, advantages and disadvantages of the available mathematical modelling frameworks covering Boolean, dynamic kinetic, single-cell models, which have been used to study processes and phenomena ranging from crosstalks to sources of cell-to-cell variability in the context of SNF1 signalling. Based on the lessons from existing models, we further discuss how to develop a consensus dynamic mechanistic model of the entire SNF1 pathway that can provide novel insights into the dynamics of nutrient signalling.

D-Xylose Sensing in Saccharomyces cerevisiae: Insights from D-Glucose Signaling and Native D-Xylose Utilizers

Extension of the substrate range is among one of the metabolic engineering goals for microorganisms used in biotechnological processes because it enables the use of a wide range of raw materials as substrates. One of the most prominent examples is the engineering of baker's yeast Saccharomyces cerevisiae for the utilization of d-xylose, a five-carbon sugar found in high abundance in lignocellulosic biomass and a key substrate to achieve good process economy in chemical production from renewable and non-edible plant feedstocks. Despite many excellent engineering strategies that have allowed recombinant S. cerevisiae to ferment d-xylose to ethanol at high yields, the consumption rate of d-xylose is still significantly lower than that of its preferred sugar d-glucose. In mixed d-glucose/d-xylose cultivations, d-xylose is only utilized after d-glucose depletion, which leads to prolonged process times and added costs. Due to this limitation, the response on d-xylose in the native sugar signaling pathways has emerged as a promising next-level engineering target. Here we review the current status of the knowledge of the response of S. cerevisiae signaling pathways to d-xylose. To do this, we first summarize the response of the native sensing and signaling pathways in S. cerevisiae to d-glucose (the preferred sugar of the yeast). Using the d-glucose case as a point of reference, we then proceed to discuss the known signaling response to d-xylose in S. cerevisiae and current attempts of improving the response by signaling engineering using native targets and synthetic (non-native) regulatory circuits.

Multiscale models quantifying yeast physiology: towards a whole-cell model

The yeast Saccharomyces cerevisiae is widely used as a cell factory and as an important eukaryal model organism for studying cellular physiology related to human health and disease. Yeast was also the first eukaryal organism for which a genome-scale metabolic model (GEM) was developed. In recent years there has been interest in expanding the modeling framework for yeast by incorporating enzymatic parameters and other heterogeneous cellular networks to obtain a more comprehensive description of cellular physiology. We review the latest developments in multiscale models of yeast, and illustrate how a new generation of multiscale models could significantly enhance the predictive performance and expand the applications of classical GEMs in cell factory design and basic studies of yeast physiology.

Transcriptional regulatory proteins in central carbon metabolism of Pichia pastoris and Saccharomyces cerevisiae

System-wide interactions in living cells and discovery of the diverse roles of transcriptional regulatory proteins that are mediator proteins with catalytic domains and regulatory subunits and transcription factors in the cellular pathways have become crucial for understanding the cellular response to environmental conditions. This review provides information for future metabolic engineering strategies through analyses on the highly interconnected regulatory networks in Saccharomyces cerevisiae and Pichia pastoris and identifying their components. We discuss the current knowledge on the carbon catabolite repression (CCR) mechanism, interconnecting regulatory system of the central metabolic pathways that regulate cell metabolism based on nutrient availability in the industrial yeasts. The regulatory proteins and their functions in the CCR signalling pathways in both yeasts are presented and discussed. We highlight the importance of metabolic signalling networks by signifying ways on how effective engineering strategies can be designed for generating novel regulatory circuits, furthermore to activate pathways that reconfigure the network architecture. We summarize the evidence that engineering of multilayer regulation is needed for directed evolution of the cellular network by putting the transcriptional control into a new perspective for the regulation of central carbon metabolism of the industrial yeasts; furthermore, we suggest research directions that may help to enhance production of recombinant products in the widely used, creatively engineered, but relatively less studied P. pastoris through de novo metabolic engineering strategies based on the discovery of components of signalling pathways in CCR metabolism. KEY POINTS: • Transcriptional regulation and control is the key phenomenon in the cellular processes. • Designing de novo metabolic engineering strategies depends on the discovery of signalling pathways in CCR metabolism. • Crosstalk between pathways occurs through essential parts of transcriptional machinery connected to specific catalytic domains. • In S. cerevisiae, a major part of CCR metabolism is controlled through Snf1 kinase, Glc7 phosphatase, and Srb10 kinase. • In P. pastoris, signalling pathways in CCR metabolism have not yet been clearly known yet. • Cellular regulations on the transcription of promoters are controlled with carbon sources.

An abscisic acid-responsive protein interaction network for sucrose non-fermenting related kinase1 in abiotic stress response

Yeast Snf1 (Sucrose non-fermenting1), mammalian AMPK (5′ AMP-activated protein kinase) and plant SnRK1 (Snf1-Related Kinase1) are conserved heterotrimeric kinase complexes that re-establish energy homeostasis following stress. The hormone abscisic acid (ABA) plays a crucial role in plant stress response. Activation of SnRK1 or ABA signaling results in overlapping transcriptional changes, suggesting these stress pathways share common targets. To investigate how SnRK1 and ABA interact during stress response in Arabidopsis thaliana, we screened the SnRK1 complex by yeast two-hybrid against a library of proteins encoded by 258 ABA-regulated genes. Here, we identify 125 SnRK1- interacting proteins (SnIPs). Network analysis indicates that a subset of SnIPs form signaling modules in response to abiotic stress. Functional studies show the involvement of SnRK1 and select SnIPs in abiotic stress responses. This targeted study uncovers the largest set of SnRK1 interactors, which can be used to further characterize SnRK1 role in plant survival under stress.

Carianopol et al. construct a detailed protein interaction network for the SnRK1 kinase complex to investigate the interaction of SnRK1 and ABA during stress response. They identify 125 proteins that interact with SnRK1, which can be used further to characterise the role of SnRK1 in plant survival under stress.

A scalable method for parameter-free simulation and validation of mechanistic cellular signal transduction network models

The metabolic modelling community has established the gold standard for bottom-up systems biology with reconstruction, validation and simulation of mechanistic genome-scale models. Similar methods have not been established for signal transduction networks, where the representation of complexes and internal states leads to scalability issues in both model formulation and execution. While rule- and agent-based methods allow efficient model definition and execution, respectively, model parametrisation introduces an additional layer of uncertainty due to the sparsity of reliably measured parameters. Here, we present a scalable method for parameter-free simulation of mechanistic signal transduction networks. It is based on rxncon and uses a bipartite Boolean logic with separate update rules for reactions and states. Using two generic update rules, we enable translation of any rxncon model into a unique Boolean model, which can be used for network validation and simulation-allowing the prediction of system-level function directly from molecular mechanistic data. Through scalable model definition and simulation, and the independence of quantitative parameters, it opens up for simulation and validation of mechanistic genome-scale models of signal transduction networks.

Robustness of Nutrient Signaling Is Maintained by Interconnectivity Between Signal Transduction Pathways

Systems biology approaches provide means to study the interplay between biological processes leading to the mechanistic understanding of the properties of complex biological systems. Here, we developed a vector format rule-based Boolean logic model of the yeast S. cerevisiae cAMP-PKA, Snf1, and the Snf3-Rgt2 pathway to better understand the role of crosstalk on network robustness and function. We identified that phosphatases are the common unknown components of the network and that crosstalk from the cAMP-PKA pathway to other pathways plays a critical role in nutrient sensing events. The model was simulated with known crosstalk combinations and subsequent analysis led to the identification of characteristics and impact of pathway interconnections. Our results revealed that the interconnections between the Snf1 and Snf3-Rgt2 pathway led to increased robustness in these signaling pathways. Overall, our approach contributes to the understanding of the function and importance of crosstalk in nutrient signaling.

Using rxncon to Develop Rule-Based Models

We present a protocol for building, validating, and simulating models of signal transduction networks. These networks are challenging modeling targets due to the combinatorial complexity and sparse data, which have made it a major challenge even to formalize the current knowledge. To address this, the community has developed methods to model biomolecular reaction networks based on site dynamics. The strength of this approach is that reactions and states can be defined at variable resolution, which makes it possible to adapt the model resolution to the empirical data. This improves both scalability and accuracy, making it possible to formalize large models of signal transduction networks. Here, we present a method to build and validate large models of signal transduction networks. The workflow is based on rxncon, the reaction-contingency language. In a five-step process, we create a mechanistic network model, convert it into an executable Boolean model, use the Boolean model to evaluate and improve the network, and finally export the rxncon model into a rule-based format. We provide an introduction to the rxncon language and an annotated, step-by-step protocol for the workflow. Finally, we create a small model of the insulin signaling pathway to illustrate the protocol, together with some of the challenges-and some of their solutions-in modeling signal transduction.

Snf1 cooperates with the CWI MAPK pathway to mediate the degradation of Med13 following oxidative stress

Eukaryotic cells, when faced with unfavorable environmental conditions, mount either pro-survival or pro-death programs. The conserved cyclin C-Cdk8 kinase plays a key role in this decision. Both are members of the Cdk8 kinase module that, along with Med12 and Med13, associate with the core Mediator complex of RNA polymerase II. In Saccharomyces cerevisiae, oxidative stress triggers Med13 destruction, which releases cyclin C into the cytoplasm to promote mitochondrial fission and programmed cell death. The SCF^Grr1 ubiquitin ligase mediates Med13 degradation dependent on the cell wall integrity pathway, MAPK Slt2. Here we show that the AMP kinase Snf1 activates a second SCF^Grr1 responsive degron in Med13. Deletion of Snf1 resulted in nuclear retention of cyclin C and failure to induce mitochondrial fragmentation. This degron was able to confer oxidative-stress-induced destruction when fused to a heterologous protein in a Snf1 dependent manner. Although snf1∆ mutants failed to destroy Med13, deleting the degron did not prevent destruction. These results indicate that the control of Med13 degradation following H₂O₂ stress is complex, being controlled simultaneously by CWI and MAPK pathways.