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.

kboolnet: a toolkit for the verification, validation, and visualization of reaction-contingency (rxncon) models.

BACKGROUND: Computational models of cell signaling networks are extremely useful tools for the exploration of underlying system behavior and prediction of response to various perturbations. By representing signaling cascades as executable Boolean networks, the previously developed rxncon ("reaction-contingency") formalism and associated Python package enable accurate and scalable modeling of signal transduction even in large (thousands of components) biological systems. The models are split into reactions, which generate states, and contingencies, that impinge on reactions; this avoids the so-called "combinatorial explosion" of system size. Boolean description of the biological system compensates for the poor availability of kinetic parameters which are necessary for quantitative models. Unfortunately, few tools are available to support rxncon model development, especially for large, intricate systems. RESULTS: We present the kboolnet toolkit ( https://github.com/Kufalab-UCSD/kboolnet , complete documentation at https://github.com/Kufalab-UCSD/kboolnet/wiki ), an R package and a set of scripts that seamlessly integrate with the python-based rxncon software and collectively provide a complete workflow for the verification, validation, and visualization of rxncon models. The verification script VerifyModel.R checks for responsiveness to repeated stimulations as well as consistency of steady state behavior. The validation scripts TruthTable.R, SensitivityAnalysis.R, and ScoreNet.R provide various readouts for the comparison of model predictions to experimental data. In particular, ScoreNet.R compares model predictions to a cloud-stored MIDAS-format experimental database to provide a numerical score for tracking model accuracy. Finally, the visualization scripts allow for graphical representations of model topology and behavior. The entire kboolnet toolkit is cloud-enabled, allowing for easy collaborative development; most scripts also allow for the extraction and analysis of individual user-defined "modules". CONCLUSION: The kboolnet toolkit provides a modular, cloud-enabled workflow for the development of rxncon models, as well as their verification, validation, and visualization. This will enable the creation of larger, more comprehensive, and more rigorous models of cell signaling using the rxncon formalism in the future.

A yeast cell cycle model integrating stress, signaling, and physiology.

The cell division cycle in eukaryotic cells is a series of highly coordinated molecular interactions that ensure that cell growth, duplication of genetic material, and actual cell division are precisely orchestrated to give rise to two viable progeny cells. Moreover, the cell cycle machinery is responsible for incorporating information about external cues or internal processes that the cell must keep track of to ensure a coordinated, timely progression of all related processes. This is most pronounced in multicellular organisms, but also a cardinal feature in model organisms such as baker's yeast. The complex and integrative behavior is difficult to grasp and requires mathematical modeling to fully understand the quantitative interplay of the single components within the entire system. Here, we present a self-oscillating mathematical model of the yeast cell cycle that comprises all major cyclins and their main regulators. Furthermore, it accounts for the regulation of the cell cycle machinery by a series of external stimuli such as mating pheromones and changes in osmotic pressure or nutrient quality. We demonstrate how the external perturbations modify the dynamics of cell cycle components and how the cell cycle resumes after adaptation to or relief from stress.

Identification of periodic attractors in Boolean networks using a priori information.

Boolean networks (BNs) have been developed to describe various biological processes, which requires analysis of attractors, the long-term stable states. While many methods have been proposed to detection and enumeration of attractors, there are no methods which have been demonstrated to be theoretically better than the naive method and be practically used for large biological BNs. Here, we present a novel method to calculate attractors based on a priori information, which works much and verifiably faster than the naive method. We apply the method to two BNs which differ in size, modeling formalism, and biological scope. Despite these differences, the method presented here provides a powerful tool for the analysis of both networks. First, our analysis of a BN studying the effect of the microenvironment during angiogenesis shows that the previously defined microenvironments inducing the specialized phalanx behavior in endothelial cells (ECs) additionally induce stalk behavior. We obtain this result from an extended network version which was previously not analyzed. Second, we were able to heuristically detect attractors in a cell cycle control network formalized as a bipartite Boolean model (bBM) with 3158 nodes. These attractors are directly interpretable in terms of genotype-to-phenotype relationships, allowing network validation equivalent to an in silico mutagenesis screen. Our approach contributes to the development of scalable analysis methods required for whole-cell modeling efforts.

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.

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.

A comprehensive, mechanistically detailed, and executable model of the cell division cycle in Saccharomyces cerevisiae.

Understanding how cellular functions emerge from the underlying molecular mechanisms is a key challenge in biology. This will require computational models, whose predictive power is expected to increase with coverage and precision of formulation. Genome-scale models revolutionised the metabolic field and made the first whole-cell model possible. However, the lack of genome-scale models of signalling networks blocks the development of eukaryotic whole-cell models. Here, we present a comprehensive mechanistic model of the molecular network that controls the cell division cycle in Saccharomyces cerevisiae. We use rxncon, the reaction-contingency language, to neutralise the scalability issues preventing formulation, visualisation and simulation of signalling networks at the genome-scale. We use parameter-free modelling to validate the network and to predict genotype-to-phenotype relationships down to residue resolution. This mechanistic genome-scale model offers a new perspective on eukaryotic cell cycle control, and opens up for similar models-and eventually whole-cell models-of human cells.

Toward Community Standards and Software for Whole-Cell Modeling.

OBJECTIVE: Whole-cell (WC) modeling is a promising tool for biological research, bioengineering, and medicine. However, substantial work remains to create accurate comprehensive models of complex cells. METHODS: We organized the 2015 Whole-Cell Modeling Summer School to teach WC modeling and evaluate the need for new WC modeling standards and software by recoding a recently published WC model in the Systems Biology Markup Language. RESULTS: Our analysis revealed several challenges to representing WC models using the current standards. CONCLUSION: We, therefore, propose several new WC modeling standards, software, and databases. SIGNIFICANCE: We anticipate that these new standards and software will enable more comprehensive models.

The MYpop toolbox: Putting yeast stress responses in cellular context on single cell and population scales.

Systems biology holds the promise to integrate multiple sources of information in order to build ever more complete models of cellular function. To do this, the field must overcome two significant challenges. First, the current strategy to model average cells must be replaced with population based models accounting for cell-to-cell variability. Second, models must be integrated with each other and with basic cellular function. This requires a core model of cellular physiology as well as a multiscale simulation platform to support large-scale simulation of culture or tissues from single cells. Here, we present such a simulation platform with a core model of yeast physiology as scaffold to integrate and simulate SBML models. The software automates this integration helping users simulate their model of choice in context of the cell division cycle. We benchmark model merging, simulation and analysis by integrating a minimal model of osmotic stress into the core model and analyzing it. We characterize the effect of single cell differences on the dynamics of osmoadaptation, estimating when normal cell growth is resumed and obtaining an explanation for experimentally observed glycerol dynamics based on population dynamics. Hence, the platform can be used to reconcile single cell and population level data.

Strategies for structuring interdisciplinary education in Systems Biology: an European perspective.

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.

Stochastic simulation of Boolean rxncon models: towards quantitative analysis of large signaling networks.

BACKGROUND: Cellular decision-making is governed by molecular networks that are highly complex. An integrative understanding of these networks on a genome wide level is essential to understand cellular health and disease. In most cases however, such an understanding is beyond human comprehension and requires computational modeling. Mathematical modeling of biological networks at the level of biochemical details has hitherto relied on state transition models. These are typically based on enumeration of all relevant model states, and hence become very complex unless severely--and often arbitrarily--reduced. Furthermore, the parameters required for genome wide networks will remain underdetermined for the conceivable future. Alternatively, networks can be simulated by Boolean models, although these typically sacrifice molecular detail as well as distinction between different levels or modes of activity. However, the modeling community still lacks methods that can simulate genome scale networks on the level of biochemical reaction detail in a quantitative or semi quantitative manner. RESULTS: Here, we present a probabilistic bipartite Boolean modeling method that addresses these issues. The method is based on the reaction-contingency formalism, and enables fast simulation of large networks. We demonstrate its scalability by applying it to the yeast mitogen-activated protein kinase (MAPK) network consisting of 140 proteins and 608 nodes. CONCLUSION: The probabilistic Boolean model can be generated and parameterized automatically from a rxncon network description, using only two global parameters, and its qualitative behavior is robust against order of magnitude variation in these parameters. Our method can hence be used to simulate the outcome of large signal transduction network reconstruction, with little or no overhead in model creation or parameterization.

Bud-Localization of CLB2 mRNA Can Constitute a Growth Rate Dependent Daughter Sizer.

Maintenance of cellular size is a fundamental systems level process that requires balancing of cell growth with proliferation. This is achieved via the cell division cycle, which is driven by the sequential accumulation and destruction of cyclins. The regulatory network around these cyclins, particularly in G1, has been interpreted as a size control network in budding yeast, and cell size as being decisive for the START transition. However, it is not clear why disruptions in the G1 network may lead to altered size rather than loss of size control, or why the S-G2-M duration also depends on nutrients. With a mathematical population model comprised of individually growing cells, we show that cyclin translation would suffice to explain the observed growth rate dependence of cell volume at START. Moreover, we assess the impact of the observed bud-localisation of the G2 cyclin CLB2 mRNA, and find that localised cyclin translation could provide an efficient mechanism for measuring the biosynthetic capacity in specific compartments: The mother in G1, and the growing bud in G2. Hence, iteration of the same principle can ensure that the mother cell is strong enough to grow a bud, and that the bud is strong enough for independent life. Cell sizes emerge in the model, which predicts that a single CDK-cyclin pair per growth phase suffices for size control in budding yeast, despite the necessity of the cell cycle network around the cyclins to integrate other cues. Size control seems to be exerted twice, where the G2/M control affects bud size through bud-localized translation of CLB2 mRNA, explaining the dependence of the S-G2-M duration on nutrients. Taken together, our findings suggest that cell size is an emergent rather than a regulatory property of the network linking growth and proliferation.

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

BACKGROUND/OBJECTIVES: The SNF1/AMPK protein kinase has a central role in energy homeostasis in eukaryotic cells. It is activated by energy depletion and stimulates processes leading to the production of ATP while it downregulates ATP-consuming processes. The yeast SNF1 complex is best known for its role in glucose derepression. METHODS: We performed a network reconstruction of the Snf1 pathway based on a comprehensive literature review. The network was formalised in the rxncon language, and we used the rxncon toolbox for model validation and gap filling. RESULTS: We present a machine-readable network definition that summarises the mechanistic knowledge of the Snf1 pathway. Furthermore, we used the known input/output relationships in the network to identify and fill gaps in the information transfer through the pathway, to produce a functional network model. Finally, we convert the functional network model into a rule-based model as a proof-of-principle. CONCLUSIONS: The workflow presented here enables large scale reconstruction, validation and gap filling of signal transduction networks. It is analogous to but distinct from that established for metabolic networks. We demonstrate the workflow capabilities, and the direct link between the reconstruction and dynamic modelling, with the Snf1 network. This network is a distillation of the knowledge from all previous publications on the Snf1/AMPK pathway. The network is a knowledge resource for modellers and experimentalists alike, and a template for similar efforts in higher eukaryotes. Finally, we envisage the workflow as an instrumental tool for reconstruction of large signalling networks across Eukaryota.

A Monte Carlo evaluation of beam characteristics for total body irradiation at extended treatment distances.

The aim is to study beam characteristics at large distances when focusing on the electron component. In particular, to investigate the utility of spoilers with various thicknesses as an electron source, as well as the effect of different spoiler-to-surface distances (STSD) on the beam characteristics and, consequently, on the dose in the superficial region. A MC model of a 15 MV Varian accelerator, validated earlier by experimental data at isocenter and extended distances used in large-field total body irradiation, is applied to evaluate beam characteristics at distances larger than 400 cm. Calculations are carried out using BEAMnrc/DOSXYZnrc code packages and phase space data are analyzed by the beam data processor BEAMdp. The electron component of the beam is analyzed at isocenter and extended distances, with and without spoilers as beam modifiers, assuming vacuum or air surrounding the accelerator head. Spoiler thickness of 1.6 cm is found to be optimal compared to thicknesses of 0.8 cm and 2.4 cm. The STSD variations should be taken into account when treating patients, in particular when the treatment protocols are based on a fixed distance to the patient central sagittal plane, and also, in order to maintain high dose in the superficial region.

Bridging the gaps in systems biology.

Systems biology aims at creating mathematical models, i.e., computational reconstructions of biological systems and processes that will result in a new level of understanding-the elucidation of the basic and presumably conserved "design" and "engineering" principles of biomolecular systems. Thus, systems biology will move biology from a phenomenological to a predictive science. Mathematical modeling of biological networks and processes has already greatly improved our understanding of many cellular processes. However, given the massive amount of qualitative and quantitative data currently produced and number of burning questions in health care and biotechnology needed to be solved is still in its early phases. The field requires novel approaches for abstraction, for modeling bioprocesses that follow different biochemical and biophysical rules, and for combining different modules into larger models that still allow realistic simulation with the computational power available today. We have identified and discussed currently most prominent problems in systems biology: (1) how to bridge different scales of modeling abstraction, (2) how to bridge the gap between topological and mechanistic modeling, and (3) how to bridge the wet and dry laboratory gap. The future success of systems biology largely depends on bridging the recognized gaps.

Glucose de-repression by yeast AMP-activated protein kinase SNF1 is controlled via at least two independent steps.

The AMP-activated protein kinase, AMPK, controls energy homeostasis in eukaryotic cells but little is known about the mechanisms governing the dynamics of its activation/deactivation. The yeast AMPK, SNF1, is activated in response to glucose depletion and mediates glucose de-repression by inactivating the transcriptional repressor Mig1. Here we show that overexpression of the Snf1-activating kinase Sak1 results, in the presence of glucose, in constitutive Snf1 activation without alleviating glucose repression. Co-overexpression of the regulatory subunit Reg1 of the Glc-Reg1 phosphatase complex partly restores glucose regulation of Snf1. We generated a set of 24 kinetic mathematical models based on dynamic data of Snf1 pathway activation and deactivation. The models that reproduced our experimental observations best featured (a) glucose regulation of both Snf1 phosphorylation and dephosphorylation, (b) determination of the Mig1 phosphorylation status in the absence of glucose by Snf1 activity only and (c) a regulatory step directing active Snf1 to Mig1 under glucose limitation. Hence it appears that glucose de-repression via Snf1-Mig1 is regulated by glucose via at least two independent steps: the control of activation of the Snf1 kinase and directing active Snf1 to inactivating its target Mig1.

Reaction-contingency based bipartite Boolean modelling.

BACKGROUND: Intracellular signalling systems are highly complex, rendering mathematical modelling of large signalling networks infeasible or impractical. Boolean modelling provides one feasible approach to whole-network modelling, but at the cost of dequantification and decontextualisation of activation. That is, these models cannot distinguish between different downstream roles played by the same component activated in different contexts. RESULTS: Here, we address this with a bipartite Boolean modelling approach. Briefly, we use a state oriented approach with separate update rules based on reactions and contingencies. This approach retains contextual activation information and distinguishes distinct signals passing through a single component. Furthermore, we integrate this approach in the rxncon framework to support automatic model generation and iterative model definition and validation. We benchmark this method with the previously mapped MAP kinase network in yeast, showing that minor adjustments suffice to produce a functional network description. CONCLUSIONS: Taken together, we (i) present a bipartite Boolean modelling approach that retains contextual activation information, (ii) provide software support for automatic model generation, visualisation and simulation, and (iii) demonstrate its use for iterative model generation and validation.

Initiation of the transcriptional response to hyperosmotic shock correlates with the potential for volume recovery.

The control of activity and localization of transcription factors is critical for appropriate transcriptional responses. In eukaryotes, signal transduction components such as mitogen-activated protein kinase (MAPK) shuttle into the nucleus to activate transcription. It is not known in detail how different amounts of nuclear MAPK over time affect the transcriptional response. In the present study, we aimed to address this issue by studying the high osmolarity glycerol (HOG) system in Saccharomyces cerevisiae. We employed a conditional osmotic system, which changes the period of the MAPK Hog1 signal independent of the initial stress level. We determined the dynamics of the Hog1 nuclear localization and cell volume by single-cell analysis in well-controlled microfluidics systems and compared the responses with the global transcriptional output of cell populations. We discovered that the onset of the initial transcriptional response correlates with the potential of cells for rapid adaptation; cells that are capable of recovering quickly initiate the transcriptional responses immediately, whereas cells that require longer time to adapt also respond later. This is reflected by Hog1 nuclear localization, Hog1 promoter association and the transcriptional response, but not Hog1 phosphorylation, suggesting that a presently uncharacterized rapid adaptive mechanism precedes the Hog1 nuclear response. Furthermore, we found that the period of Hog1 nuclear residence affects the amplitude of the transcriptional response rather than the spectrum of responsive genes.

Information content and scalability in signal transduction network reconstruction formats.

One of the first steps towards holistic understanding of cellular networks is the integration of the available information in a human and machine readable format. This network reconstruction process is well established for metabolic networks, and numerous genome wide metabolic reconstructions are already available. Extending these strategies to signalling networks has proven difficult, primarily due to the combinatorial nature of regulatory modifications. The combinatorial nature of possible protein-protein interactions and post translational modifications affects both network size and the correspondence between the reconstructed network and the underlying empirical data. Here, we discuss different approaches to reconstruction of signal transduction networks. We divide the current approaches into topological, specific state based and reaction-contingency based, and discuss their different information content and scalability. The discussion focusses on graphical formats but the points are in general applicable also to mathematical models and databases. While the formats have complementary strengths especially for small networks, reaction-contingency based formats have a number of advantages in the light of global network reconstruction. In particular, they minimise the need for assumptions, maximise the congruence with empirical data, and scale efficiently with network size.

Biographer: web-based editing and rendering of SBGN compliant biochemical networks.

MOTIVATION: The rapid accumulation of knowledge in the field of Systems Biology during the past years requires advanced, but simple-to-use, methods for the visualization of information in a structured and easily comprehensible manner. RESULTS: We have developed biographer, a web-based renderer and editor for reaction networks, which can be integrated as a library into tools dealing with network-related information. Our software enables visualizations based on the emerging standard Systems Biology Graphical Notation. It is able to import networks encoded in various formats such as SBML, SBGN-ML and jSBGN, a custom lightweight exchange format. The core package is implemented in HTML5, CSS and JavaScript and can be used within any kind of web-based project. It features interactive graph-editing tools and automatic graph layout algorithms. In addition, we provide a standalone graph editor and a web server, which contains enhanced features like web services for the import and export of models and visualizations in different formats. AVAILABILITY: The biographer tool can be used at and downloaded from the web page http://biographer.biologie.hu-berlin.de/. The different software packages, including a server-independent version as well as a web server for Windows and Linux based systems, are available at http://code.google.com/p/biographer/ under the open-source license LGPL