Programmed cell death (PCD) is an essential cellular mechanism that is evolutionary conserved, mediated through various pathways and acts by integrating different stimuli. Many diseases such as neurodegenerative diseases and cancers are found to be caused by, or associated with, regulations in the cell death pathways. Yeast Saccharomyces cerevisiae, is a unicellular eukaryotic organism that shares with human cells components and pathways of the PCD and is therefore used as a model organism. Boolean modeling is becoming promising approach to capture qualitative behavior and describe essential properties of such complex networks. Here we present large literature-based and to our knowledge first Boolean model that combines pathways leading to apoptosis (a type of PCD) in yeast. Analysis of the yeast model confirmed experimental findings of anti-apoptotic role of Bir1p and pro-apoptotic role of Stm1p and revealed activation of the stress protein kinase Hog proposing the maximal level of activation upon heat stress. In addition we extended the yeast model and created an in silico humanized yeast in which human pro- and anti-apoptotic regulators Bcl-2 family and Valosin-contain protein (VCP) are included in the model. We showed that accumulation of Bax in silico humanized yeast shows apoptotic markers and that VCP is essential target of Akt Signaling. The presented Boolean model provides comprehensive description of yeast apoptosis network behavior. Extended model of humanized yeast gives new insights of how complex human disease like neurodegeneration can initially be tested.

The interplay between nutrient-induced signaling and metabolism plays an important role in maintaining homeostasis and its malfunction has been implicated in many different human diseases such as obesity, type 2 diabetes, cancer, and neurological disorders. Therefore, unraveling the role of nutrients as signaling molecules and metabolites together with their interconnectivity may provide a deeper understanding of how these conditions occur. Both signaling and metabolism have been extensively studied using various systems biology approaches. However, they are mainly studied individually and in addition, current models lack both the complexity of the dynamics and the effects of the crosstalk in the signaling system. To gain a better understanding of the interconnectivity between nutrient signaling and metabolism in yeast cells, we developed a hybrid model, combining a Boolean module, describing the main pathways of glucose and nitrogen signaling, and an enzyme-constrained model accounting for the central carbon metabolism ofSaccharomyces cerevisiae, using a regulatory network as a link. The resulting hybrid model was able to capture a diverse utalization of isoenzymes and to our knowledge outperforms constraint-based models in the prediction of individual enzymes for both respiratory and mixed metabolism. The model showed that during fermentation, enzyme utilization has a major contribution in governing protein allocation, while in low glucose conditions robustness and control are prioritized. In addition, the model was capable of reproducing the regulatory effects that are associated with the Crabtree effect and glucose repression, as well as regulatory effects associated with lifespan increase during caloric restriction. Overall, we show that our hybrid model provides a comprehensive framework for the study of the non-trivial effects of the interplay between signaling and metabolism, suggesting connections between the Snf1 signaling pathways and processes that have been related to chronological lifespan of yeast cells.

Systems Biology is an approach to biology and medicine that has the potential to lead to a better understanding of how biological properties emerge from the interaction of genes, proteins, molecules, cells and organisms. The approach aims at elucidating how these interactions govern biological function by employing experimental data, mathematical models and computational simulations. As Systems Biology is inherently multidisciplinary, education within this field meets numerous hurdles including departmental barriers, availability of all required expertise locally, appropriate teaching material and example curricula. As university education at the Bachelor’s level is traditionally built upon disciplinary degrees, we believe that the most effective way to implement education in Systems Biology would be at the Master’s level, as it offers a more flexible framework. Our team of experts and active performers of Systems Biology education suggest here (i) a definition of the skills that students should acquire within a Master’s programme in Systems Biology, (ii) a possible basic educational curriculum with flexibility to adjust to different application areas and local research strengths, (iii) a description of possible career paths for students who undergo such an education, (iv) conditions that should improve the recruitment of students to such programmes and (v) mechanisms for collaboration and excellence spreading among education professionals. With the growing interest of industry in applying Systems Biology approaches in their fields, a concerted action between academia and industry is needed to build this expertise. Here we present a reflection of the European situation and expertise, where most of the challenges we discuss are universal, anticipating that our suggestions will be useful internationally. We believe that one of the overriding goals of any Systems Biology education should be a student’s ability to phrase and communicate research questions in such a manner that they can be solved by the integration of experiments and modelling, as well as to communicate and collaborate productively across different experimental and theoretical disciplines in research and development.

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. 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. 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.

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.