Chemical Reaction Network Theory elucidates sources of multistability in interferon signaling.
Ontology highlight
ABSTRACT:
SUBMITTER: Otero-Muras I
PROVIDER: S-EPMC5400276 | biostudies-literature | 2017 Apr
REPOSITORIES: biostudies-literature
Publications
Chemical Reaction Network Theory elucidates sources of multistability in interferon signaling.
Otero-Muras Irene I Yordanov Pencho P Stelling Joerg J
PLoS computational biology 20170403 4
Bistability has important implications in signaling pathways, since it indicates a potential cell decision between alternative outcomes. We present two approaches developed in the framework of the Chemical Reaction Network Theory for easy and efficient search of multiple steady state behavior in signaling networks (both with and without mass conservation), and apply them to search for sources of bistability at different levels of the interferon signaling pathway. Different type I interferon subt ...[more]
Similar Datasets
Project description:BackgroundReaction networks (RNs) comprise a set X of species and a set [Formula: see text] of reactions [Formula: see text], each converting a multiset of educts [Formula: see text] into a multiset [Formula: see text] of products. RNs are equivalent to directed hypergraphs. However, not all RNs necessarily admit a chemical interpretation. Instead, they might contradict fundamental principles of physics such as the conservation of energy and mass or the reversibility of chemical reactions. The consequences of these necessary conditions for the stoichiometric matrix [Formula: see text] have been discussed extensively in the chemical literature. Here, we provide sufficient conditions for [Formula: see text] that guarantee the interpretation of RNs in terms of balanced sum formulas and structural formulas, respectively.ResultsChemically plausible RNs allow neither a perpetuum mobile, i.e., a "futile cycle" of reactions with non-vanishing energy production, nor the creation or annihilation of mass. Such RNs are said to be thermodynamically sound and conservative. For finite RNs, both conditions can be expressed equivalently as properties of the stoichiometric matrix [Formula: see text]. The first condition is vacuous for reversible networks, but it excludes irreversible futile cycles and-in a stricter sense-futile cycles that even contain an irreversible reaction. The second condition is equivalent to the existence of a strictly positive reaction invariant. It is also sufficient for the existence of a realization in terms of sum formulas, obeying conservation of "atoms". In particular, these realizations can be chosen such that any two species have distinct sum formulas, unless [Formula: see text] implies that they are "obligatory isomers". In terms of structural formulas, every compound is a labeled multigraph, in essence a Lewis formula, and reactions comprise only a rearrangement of bonds such that the total bond order is preserved. In particular, for every conservative RN, there exists a Lewis realization, in which any two compounds are realized by pairwisely distinct multigraphs. Finally, we show that, in general, there are infinitely many realizations for a given conservative RN.Conclusions"Chemical" RNs are directed hypergraphs with a stoichiometric matrix [Formula: see text] whose left kernel contains a strictly positive vector and whose right kernel does not contain a futile cycle involving an irreversible reaction. This simple characterization also provides a concise specification of random models for chemical RNs that additionally constrain [Formula: see text] by rank, sparsity, or distribution of the non-zero entries. Furthermore, it suggests several interesting avenues for future research, in particular, concerning alternative representations of reaction networks and infinite chemical universes.
| S-EPMC9484159 | biostudies-literature
Project description:BACKGROUND:The transitions between epithelial (E) and mesenchymal (M) cell phenotypes are essential in many biological processes like tissue development and cancer metastasis. Previous studies, both modeling and experimental, suggested that in addition to E and M states, the network responsible for these phenotypes exhibits intermediate phenotypes between E and M states. The number and importance of such states is subject to intense discussion in the epithelial-mesenchymal transition (EMT) community. RESULTS:Previous modeling efforts used traditional bifurcation analysis to explore the number of the steady states that correspond to E, M and intermediate states by varying one or two parameters at a time. Since the system has dozens of parameters that are largely unknown, it remains a challenging problem to fully describe the potential set of states and their relationship across all parameters. We use the computational tool DSGRN (Dynamic Signatures Generated by Regulatory Networks) to explore the intermediate states of an EMT model network by computing summaries of the dynamics across all of parameter space. We find that the only attractors in the system are equilibria, that E and M states dominate across parameter space, but that bistability and multistability are common. Even at extreme levels of some of the known inducers of the transition, there is a certain proportion of the parameter space at which an E or an M state co-exists with other stable steady states. CONCLUSIONS:Our results suggest that the multistability is broadly present in the EMT network across parameters and thus response of cells to signals may strongly depend on the particular cell line and genetic background.
| S-EPMC7041120 | biostudies-literature
Project description:An elastic network model (ENM), usually C? coarse-grained one, has been widely used to study protein dynamics as an alternative to classical molecular dynamics simulation. This simple approach dramatically saves the computational cost, but sometimes fails to describe a feasible conformational change due to unrealistically excessive spring connections. To overcome this limitation, we propose a mass-weighted chemical elastic network model (MWCENM) in which the total mass of each residue is assumed to be concentrated on the representative alpha carbon atom and various stiffness values are precisely assigned according to the types of chemical interactions. We test MWCENM on several well-known proteins of which both closed and open conformations are available as well as three ?-helix rich proteins. Their normal mode analysis reveals that MWCENM not only generates more plausible conformational changes, especially for closed forms of proteins, but also preserves protein secondary structures thus distinguishing MWCENM from traditional ENMs. In addition, MWCENM also reduces computational burden by using a more sparse stiffness matrix.
| S-EPMC3649262 | biostudies-literature
Project description:The organization and management of large amounts of data has become a major point in almost all areas of human knowledge. In this context, semantic approaches propose a structure for the target data, defining ontologies that state the types of entities on a certain field and how these entities are interrelated. In this work, we introduce OntoRXN, a novel ontology describing the reaction networks constructed from computational chemistry calculations. Under our paradigm, these networks are handled as undirected graphs, without assuming any traversal direction. From there, we propose a core class structure including reaction steps, network stages, chemical species, and the lower-level entities for the individual computational calculations. These individual calculations are founded on the OntoCompChem ontology and on the ioChem-BD database, where information is parsed and stored in CML format. OntoRXN is introduced through several examples in which knowledge graphs based on the ontology are generated for different chemical systems available on ioChem-BD. Finally, the resulting knowledge graphs are explored through SPARQL queries, illustrating the power of the semantic approach to standardize the analysis of intricate datasets and to simplify the development of complex workflows.
| S-EPMC9153116 | biostudies-literature
Project description:Interferon gamma (IFN-γ), is a cytokine, which is an important regulator of host defense system by mediating both innate and adaptive immune responses. IFN-γ signaling is primarily associated with inflammation and cell-mediated immune responses. IFN-γ is also represented as antitumor cytokine which facilitates immunosurveillance in tumor cells. In addition, IFN-γ mediated signaling also elicits pro-tumorigenic transformations and promotes tumor progression. Impact of IFN-γ signaling in mammalian cells has been widely studied which indicate that IFN-γ orchestrates distinct cellular functions including immunomodulation, leukocyte trafficking, apoptosis, anti-microbial, and both anti- and pro-tumorigenic role. However, a detailed network of IFN-γ signaling pathway is currently lacking. Therefore, we systematically curated the literature information pertaining to IFN-γ signaling and develop a comprehensive signaling network to facilitate better understanding of IFN-γ mediated signaling. A total of 124 proteins were catalogued that were experimentally proven to be involved in IFN-γ signaling cascade. These 124 proteins were found to participate in 81 protein-protein interactions, 94 post-translational modifications, 20 translocation events, 54 activation/inhibiton reactions. Further, 236 differential expressed genes were also documented in IFN-γ mediated signaling. IFN-γ signaling pathway is made freely available to scientific audience through NetPath at ( http://www.netpath.org/pathways?path_id=NetPath_32 ). We believe that documentation of reactions pertaining to IFN-γ signaling and development of pathway map will facilitate further research in IFN-γ associated human diseases including cancer.
| S-EPMC6235777 | biostudies-literature
Project description:Early fMRI studies suggested that brain areas processing self-related and other-related information were highly overlapping. Hypothesising functional localisation of the cortex, researchers have tried to locate "self-specific" and "other-specific" regions within these overlapping areas by subtracting suspected confounding signals in task-based fMRI experiments. Inspired by recent advances in whole-brain dynamic modelling, we instead explored an alternative hypothesis that similar spatial activation patterns could be associated with different processing modes in the form of different synchronisation patterns. Combining an automated synthesis of fMRI data with a presumption-free diffusion spectrum image (DSI) fibre-tracking algorithm, we isolated a network putatively composed of brain areas and white matter tracts involved in self-other processing. We sampled synchronisation patterns from the dynamical systems of this network using various combinations of physiological parameters. Our results showed that the self-other processing network, with simulated gamma-band activity, tended to stabilise at a number of distinct synchronisation patterns. This phenomenon, termed "multistability," could serve as an alternative model in theorising the mechanism of processing self-other information.
| S-EPMC5335699 | biostudies-literature
Project description:The stability of synchronised networked systems is a multi-faceted challenge for many natural and technological fields, from cardiac and neuronal tissue pacemakers to power grids. For these, the ongoing transition to distributed renewable energy sources leads to a proliferation of dynamical actors. The desynchronisation of a few or even one of those would likely result in a substantial blackout. Thus the dynamical stability of the synchronous state has become a leading topic in power grid research. Here we uncover that, when taking into account physical losses in the network, the back-reaction of the network induces new exotic solitary states in the individual actors and the stability characteristics of the synchronous state are dramatically altered. These effects will have to be explicitly taken into account in the design of future power grids. We expect the results presented here to transfer to other systems of coupled heterogeneous Newtonian oscillators.
| S-EPMC6992754 | biostudies-literature
Project description:BackgroundSignaling networks are designed to sense an environmental stimulus and adapt to it. We propose and study a minimal model of signaling network that can sense and respond to external stimuli of varying strength in an adaptive manner. The structure of this minimal network is derived based on some simple assumptions on its differential response to external stimuli.MethodologyWe employ stochastic differential equations and probability distributions obtained from stochastic simulations to characterize differential signaling response in our minimal network model. Gillespie's stochastic simulation algorithm (SSA) is used in this study.Conclusions/significanceWe show that the proposed minimal signaling network displays two distinct types of response as the strength of the stimulus is decreased. The signaling network has a deterministic part that undergoes rapid activation by a strong stimulus in which case cell-to-cell fluctuations can be ignored. As the strength of the stimulus decreases, the stochastic part of the network begins dominating the signaling response where slow activation is observed with characteristic large cell-to-cell stochastic variability. Interestingly, this proposed stochastic signaling network can capture some of the essential signaling behaviors of a complex apoptotic cell death signaling network that has been studied through experiments and large-scale computer simulations. Thus we claim that the proposed signaling network is an appropriate minimal model of apoptosis signaling. Elucidating the fundamental design principles of complex cellular signaling pathways such as apoptosis signaling remains a challenging task. We demonstrate how our proposed minimal model can help elucidate the effect of a specific apoptotic inhibitor Bcl-2 on apoptotic signaling in a cell-type independent manner. We also discuss the implications of our study in elucidating the adaptive strategy of cell death signaling pathways.
| S-EPMC2920308 | biostudies-literature
Project description:BackgroundWhile most game theoretical models assume that individuals randomly interact with all other group members, strong evidence indicates that individuals tend to preferentially interact with some of them. The position of an individual in a network affects, among other factors related to survival, its predation risk and competitive success. Here I then modified the Hawk-Dove game to explore the effect of social network structure on competitive strategy of individuals that differ in their fighting ability and may adjust their use of the Hawk, Dove and Assessor tactics to maximize their foraging success when they meet opponents they are connected with.ResultsFrom randomly generated networks, I demonstrate that phenotypic assortment by fighting ability reduces individuals' aggressiveness and, as such, favours cooperative interactions. Furthermore, the success of individuals with the weakest fighting ability is usually highest within networks where they most frequently meet opponents with the same fighting ability as their own, suggesting they might benefit from breaking connections with strong contestants. This might be the case when strong contestants systematically rely on the aggressive Hawk tactic or the risk of being predated is low and independent of the number of neighbours. Thus, I extended the model and built a dynamic model to allow individuals not only to adjust their behaviour to local conditions but also to modify the structure of the social network. The number of connections and degree of phenotypic assortment are then affected by ecological factors (e.g. resources value and predation risk), but above all by whether individuals can reliably assess the competitive ability of their opponents and adjust their behaviour accordingly.ConclusionsThese findings provide strong evidence that behaviour can play a key role in shaping network structure and highlight the importance of considering the coevolution of network and behaviour to apprehend its consequences on population dynamics.
| S-EPMC11365163 | biostudies-literature
Project description:The goal of many single-cell studies on eukaryotic cells is to gain insight into the biochemical reactions that control cell fate and state. In this paper we introduce the concept of Effective Stoichiometric Spaces (ESS) to guide the reconstruction of biochemical networks from multiplexed, fixed time-point, single-cell data. In contrast to methods based solely on statistical models of data, the ESS method leverages the power of the geometric theory of toric varieties to begin unraveling the structure of chemical reaction networks (CRN). This application of toric theory enables a data-driven mapping of covariance relationships in single-cell measurements into stoichiometric information, one in which each cell subpopulation has its associated ESS interpreted in terms of CRN theory. In the development of ESS we reframe certain aspects of the theory of CRN to better match data analysis. As an application of our approach we process cytomery- and image-based single-cell datasets and identify differences in cells treated with kinase inhibitors. Our approach is directly applicable to data acquired using readily accessible experimental methods such as Fluorescence Activated Cell Sorting (FACS) and multiplex immunofluorescence.
| S-EPMC6919632 | biostudies-literature