The adaptive responses of a living cell to internal and external signals are controlled by networks of proteins whose interactions are so complex that the functional integration of the network cannot be comprehended by intuitive reasoning alone. Mathematical modeling, based on biochemical rate equations, provides a rigorous and reliable tool for unraveling the complexities of molecular regulatory networks. The budding yeast cell cycle is a challenging test case for this approach, because the control system is known in exquisite detail and its function is constrained by the phenotypic properties of >100 genetically engineered strains. We show that a mathematical model built on a consensus picture of this control system is largely successful in explaining the phenotypes of mutants described so far. A few inconsistencies between the model and experiments indicate aspects of the mechanism that require revision. In addition, the model allows one to frame and critique hypotheses about how the division cycle is regulated in wild-type and mutant cells, to predict the phenotypes of new mutant combinations, and to estimate the effective values of biochemical rate constants that are difficult to measure directly in vivo.

Cytokines such as TNF and FASL can trigger death or survival depending on cell lines and cellular conditions. The mechanistic details of how a cell chooses among these cell fates are still unclear. The understanding of these processes is important since they are altered in many diseases, including cancer and AIDS. Using a discrete modelling formalism, we present a mathematical model of cell fate decision recapitulating and integrating the most consistent facts extracted from the literature. This model provides a generic high-level view of the interplays between NFκB pro-survival pathway, RIP1-dependent necrosis, and the apoptosis pathway in response to death receptor-mediated signals. Wild type simulations demonstrate robust segregation of cellular responses to receptor engagement. Model simulations recapitulate documented phenotypes of protein knockdowns and enable the prediction of the effects of novel knockdowns. In silico experiments simulate the outcomes following ligand removal at different stages, and suggest experimental approaches to further validate and specialise the model for particular cell types. We also propose a reduced conceptual model implementing the logic of the decision process. This analysis gives specific predictions regarding cross-talks between the three pathways, as well as the transient role of RIP1 protein in necrosis, and confirms the phenotypes of novel perturbations. Our wild type and mutant simulations provide novel insights to restore apoptosis in defective cells. The model analysis expands our understanding of how cell fate decision is made. Moreover, our current model can be used to assess contradictory or controversial data from the literature. Ultimately, it constitutes a valuable reasoning tool to delineate novel experiments. Activation of death receptors (TNFR and Fas) can trigger either survival or cell death according to the cell type and the cellular conditions. In other words, the same signal can have antagonist responses. On one hand, the cell can survive by activating the NFκB signalling pathway. On the other hand, it can die by apoptosis or necrosis. Apoptosis is a suicide mechanism, i.e., an orchestrated way to disrupt cellular components and pack them into specialized vesicles that can be easily removed from the environment, whereas necrosis is a type of death that involves release of intracellular components in the surrounding tissues, possibly causing inflammatory response and severe injury. We, biologists and theoreticians, have recapitulated and integrated known biological data from the literature into an influence diagram describing the molecular events leading to each possible outcome. The diagram has been translated into a dynamical Boolean model. Simulations of wild type, mutant cells and drug treatments qualitatively match current data, and predict several novel mutant phenotypes, along with general characteristics of the cell fate decision mechanism: transient activation of some key proteins in necrosis, mutual inhibitory cross-talks between the three pathways. Our model can further be used to assess contradictory data and address specific biological questions through in silico experiments.

Summary: BIOCHAM (the BIOCHemical Abstract Machine) is a software environment for modeling biochemical systems. It is based on two aspects: (1) the analysis and simulation of boolean, kinetic and stochastic models and (2) the formalization of biological properties in temporal logic. BIOCHAM provides tools and languages for describing protein networks with a simple and straightforward syntax, and for integrating biological properties into the model. It then becomes possible to analyze, query, verify and maintain the model with respect to those properties. For kinetic models, BIOCHAM can search for appropriate parameter values in order to reproduce a specific behavior observed in experiments and formalized in temporal logic. Coupled with other methods such as bifurcation diagrams, this search assists the modeler/biologist in the modeling process. Availability: BIOCHAM (v. 2.5) is a free software available for download, with example models, at Contact:Sylvain.Soliman@inria.fr

The Mitogen-Activated Protein Kinase (MAPK) network consists of tightly interconnected signalling pathways involved in diverse cellular processes, such as cell cycle, survival, apoptosis and differentiation. Although several studies reported the involvement of these signalling cascades in cancer deregulations, the precise mechanisms underlying their influence on the balance between cell proliferation and cell death (cell fate decision) in pathological circumstances remain elusive. Based on an extensive analysis of published data, we have built a comprehensive and generic reaction map for the MAPK signalling network, using CellDesigner software. In order to explore the MAPK responses to different stimuli and better understand their contributions to cell fate decision, we have considered the most crucial components and interactions and encoded them into a logical model, using the software GINsim. Our logical model analysis particularly focuses on urinary bladder cancer, where MAPK network deregulations have often been associated with specific phenotypes. To cope with the combinatorial explosion of the number of states, we have applied novel algorithms for model reduction and for the compression of state transition graphs, both implemented into the software GINsim. The results of systematic simulations for different signal combinations and network perturbations were found globally coherent with published data. In silico experiments further enabled us to delineate the roles of specific components, cross-talks and regulatory feedbacks in cell fate decision. Finally, tentative proliferative or anti-proliferative mechanisms can be connected with established bladder cancer deregulations, namely Epidermal Growth Factor Receptor (EGFR) over-expression and Fibroblast Growth Factor Receptor 3 (FGFR3) activating mutations. Depending on environmental conditions, strongly intertwined cellular signalling pathways are activated, involving activation/inactivation of proteins and genes in response to external and/or internal stimuli. Alterations of some components of these pathways can lead to wrong cell behaviours. For instance, cancer-related deregulations lead to high proliferation of malignant cells enabling sustained tumour growth. Understanding the precise mechanisms underlying these pathways is necessary to delineate efficient therapeutical approaches for each specific tumour type. We particularly focused on the Mitogen-Activated Protein Kinase (MAPK) signalling network, whose involvement in cancer is well established, although the precise conditions leading to its positive or negative influence on cell proliferation are still poorly understood. We tackled this problem by first collecting sparse published biological information into a comprehensive map describing the MAPK network in terms of stylised chemical reactions. This information source was then used to build a dynamical Boolean model recapitulating network responses to characteristic stimuli observed in selected bladder cancers. Systematic model simulations further allowed us to link specific network components and interactions with proliferative/anti-proliferative cell responses.

Modeling of signaling pathways is an important step towards the understanding and the treatment of diseases such as cancers, HIV or auto-immune diseases. MaBoSS is a software that allows to simulate populations of cells and to model stochastically the intracellular mechanisms that are deregulated in diseases. MaBoSS provides an output of a Boolean model in the form of time-dependent probabilities, for all biological entities (genes, proteins, phenotypes, etc.) of the model. We present a new version of MaBoSS (2.0), including an updated version of the core software and an environment. With this environment, the needs for modeling signaling pathways are facilitated, including model construction, visualization, simulations of mutations, drug treatments and sensitivity analyses. It offers a framework for automated production of theoretical predictions. MaBoSS software can be found at https://maboss.curie.fr, including tutorials on existing models and examples of models. Supplementary data are available at Bioinformatics online.