Siphons in Chemical Reaction Networks

Siphons in a chemical reaction system are subsets of the species that have the potential of being absent in a steady state. We present a characterization of minimal siphons in terms of primary decomposition of binomial ideals, we explore the underlying geometry, and we demonstrate the effective computation of siphons using computer algebra software. This leads to a new method for determining whether given initial concentrations allow for various boundary steady states.     

Product-Form Stationary Distributions for Deficiency Zero Chemical Reaction Networks

We consider stochastically modeled chemical reaction systems with mass-action kinetics and prove that a product-form stationary distribution exists for each closed, irreducible subset of the state space if an analogous deterministically modeled system with mass-action kinetics admits a complex balanced equilibrium. Feinberg’s deficiency zero theorem then implies that such a distribution exists so long as the corresponding chemical network is weakly reversible and has a deficiency of zero. The main parameter of the stationary distribution for the stochastically modeled system is a complex balanced equilibrium value for the corresponding deterministically modeled system. We also generalize our main result to some non-mass-action kinetics.     

Evolution of Autocatalytic Sets in Computational Models of Chemical Reaction Networks

Several computational models of chemical reaction networks have been presented in the literature in the past, showing the appearance and (potential) evolution of autocatalytic sets. However, the notion of autocatalytic sets has been defined differently in different modeling contexts, each one having some shortcoming or limitation. Here, we review four such models and definitions, and then formally describe and analyze them in the context of a mathematical framework for studying autocatalytic sets known as RAF theory. The main results are that: (1) RAF theory can capture the various previous definitions of autocatalytic sets and is therefore more complete and general, (2) the formal framework can be used to efficiently detect and analyze autocatalytic sets in all of these different computational models, (3) autocatalytic (RAF) sets are indeed likely to appear and evolve in such models, and (4) this could have important implications for a possible metabolism-first scenario for the origin of life.     

A Stronger Necessary Condition for the Multistationarity of Chemical Reaction Networks

Biochemical reaction networks grow bigger and bigger, fed by the high-throughput data provided by biologists and bred in open repositories of models allowing merging and evolution. Nevertheless, since the available data is still very far from permitting the identification of the increasing number of kinetic parameters of such models, the necessity of structural analyses for describing the dynamics of chemical networks appears stronger every day. Using the structural information, notably from the stoichiometric matrix, of a biochemical reaction system, we state a more strict version of the famous Thomas’ necessary condition for multistationarity. In particular, the obvious cases where Thomas’ condition was trivially satisfied, mutual inhibition due to a multimolecular reaction and mutual activation due to a reversible reaction, can now easily be ruled out. This more strict condition shall not be seen as some version of Thomas’ circuit functionality for the continuous case but rather as related and complementary to the whole domain of the structural analysis of (bio)chemical reaction systems, as pioneered by the chemical reaction network theory.     

Exact Probability Distributions of Selected Species in Stochastic Chemical Reaction Networks

Chemical reactions are discrete, stochastic events. As such, the species’ molecular numbers can be described by an associated master equation. However, handling such an equation may become difficult due to the large size of reaction networks. A commonly used approach to forecast the behaviour of reaction networks is to perform computational simulations of such systems and analyse their outcome statistically. This approach, however, might require high computational costs to provide accurate results. In this paper we opt for an analytical approach to obtain the time-dependent solution of the Chemical Master Equation for selected species in a general reaction network. When the reaction networks are composed exclusively of zeroth and first-order reactions, this analytical approach significantly alleviates the computational burden required by simulation-based methods. By building upon these analytical solutions, we analyse a general monomolecular reaction network with an arbitrary number of species to obtain the exact marginal probability distribution for selected species. Additionally, we study two particular topologies of monomolecular reaction networks, namely (i) an unbranched chain of monomolecular reactions with and without synthesis and degradation reactions and (ii) a circular chain of monomolecular reactions. We illustrate our methodology and alternative ways to use it for non-linear systems by analysing a protein autoactivation mechanism. Later, we compare the computational load required for the implementation of our results and a pure computational approach to analyse an unbranched chain of monomolecular reactions. Finally, we study calcium ions gates in the sarco/endoplasmic reticulum mediated by ryanodine receptors.     

Steady State Detection of Chemical Reaction Networks Using a Simplified Analytical Method

Chemical reaction networks (CRNs) are susceptible to mathematical modelling. The dynamic behavior of CRNs can be investigated by solving the polynomial equations derived from its structure. However, simple CRN give rise to non-linear polynomials that are difficult to resolve. Here we propose a procedure to locate the steady states of CRNs from a formula derived through algebraic geometry methods. We have applied this procedure to define the steady states of a classic CRN that exhibits instability, and to a model of programmed cell death.     

Catalysis in Reaction Networks

We define catalytic networks as chemical reaction networks with an essentially catalytic reaction pathway: one which is “on” in the presence of certain catalysts and “off” in their absence. We show that examples of catalytic networks include synthetic DNA molecular circuits that have been shown to perform signal amplification and molecular logic. Recall that a critical siphon is a subset of the species in a chemical reaction network whose absence is forward invariant and stoichiometrically compatible with a positive point. Our main theorem is that all weakly-reversible networks with critical siphons are catalytic. Consequently, we obtain new proofs for the persistence of atomic event-systems of Adleman et al., and normal networks of Gnacadja. We define autocatalytic networks, and conjecture that a weakly-reversible reaction network has critical siphons if and only if it is autocatalytic.     

Autocatalysis in Reaction Networks

The persistence conjecture is a long-standing open problem in chemical reaction network theory. It concerns the behavior of solutions to coupled ODE systems that arise from applying mass-action kinetics to a network of chemical reactions. The idea is that if all reactions are reversible in a weak sense, then no species can go extinct. A notion that has been found useful in thinking about persistence is that of “critical siphon.” We explore the combinatorics of critical siphons, with a view toward the persistence conjecture. We introduce the notions of “drainable” and “self-replicable” (or autocatalytic) siphons. We show that: Every minimal critical siphon is either drainable or self-replicable; reaction networks without drainable siphons are persistent; and nonautocatalytic weakly reversible networks are persistent. Our results clarify that the difficulties in proving the persistence conjecture are essentially due to competition between drainable and self-replicable siphons.     

Thermodynamics of Random Reaction Networks

Reaction networks are useful for analyzing reaction systems occurring in chemistry, systems biology, or Earth system science. Despite the importance of thermodynamic disequilibrium for many of those systems, the general thermodynamic properties of reaction networks are poorly understood. To circumvent the problem of sparse thermodynamic data, we generate artificial reaction networks and investigate their non-equilibrium steady state for various boundary fluxes. We generate linear and nonlinear networks using four different complex network models (Erdős-Rényi, Barabási-Albert, Watts-Strogatz, Pan-Sinha) and compare their topological properties with real reaction networks. For similar boundary conditions the steady state flow through the linear networks is about one order of magnitude higher than the flow through comparable nonlinear networks. In all networks, the flow decreases with the distance between the inflow and outflow boundary species, with Watts-Strogatz networks showing a significantly smaller slope compared to the three other network types. The distribution of entropy production of the individual reactions inside the network follows a power law in the intermediate region with an exponent of circa −1.5 for linear and −1.66 for nonlinear networks. An elevated entropy production rate is found in reactions associated with weakly connected species. This effect is stronger in nonlinear networks than in the linear ones. Increasing the flow through the nonlinear networks also increases the number of cycles and leads to a narrower distribution of chemical potentials. We conclude that the relation between distribution of dissipation, network topology and strength of disequilibrium is nontrivial and can be studied systematically by artificial reaction networks.     

Dynamical properties of Discrete Reaction Networks

Reaction networks are commonly used to model the dynamics of populations subject to transformations that follow an imposed stoichiometry. This paper focuses on the efficient characterisation of dynamical properties of Discrete Reaction Networks (DRNs). DRNs can be seen as modeling the underlying discrete nondeterministic transitions of stochastic models of reaction networks. In that sense, a proof of non-reachability in a given DRN has immediate implications for any concrete stochastic model based on that DRN, independent of the choice of kinetic laws and constants. Moreover, if we assume that stochastic kinetic rates are given by the mass-action law (or any other kinetic law that gives non-vanishing probability to each reaction if the required number of interacting substrates is present), then reachability properties are equivalent in the two settings. The analysis of two types of global dynamical properties of DRNs is addressed: irreducibility, i.e., the ability to reach any discrete state from any other state; and recurrence, i.e., the ability to return to any initial state. Our results consider both the verification of such properties when species are present in a large copy number, and in the general case. The necessary and sufficient conditions obtained involve algebraic conditions on the network reactions which in most cases can be verified using linear programming. Finally, the relationship of DRN irreducibility and recurrence with dynamical properties of stochastic and continuous models of reaction networks is discussed.     

Evolution of Associative Learning in Chemical Networks

Organisms that can learn about their environment and modify their behaviour appropriately during their lifetime are more likely to survive and reproduce than organisms that do not. While associative learning – the ability to detect correlated features of the environment – has been studied extensively in nervous systems, where the underlying mechanisms are reasonably well understood, mechanisms within single cells that could allow associative learning have received little attention. Here, using in silico evolution of chemical networks, we show that there exists a diversity of remarkably simple and plausible chemical solutions to the associative learning problem, the simplest of which uses only one core chemical reaction. We then asked to what extent a linear combination of chemical concentrations in the network could approximate the ideal Bayesian posterior of an environment given the stimulus history so far? This Bayesian analysis revealed the ‘memory traces’ of the chemical network. The implication of this paper is that there is little reason to believe that a lack of suitable phenotypic variation would prevent associative learning from evolving in cell signalling, metabolic, gene regulatory, or a mixture of these networks in cells.     

Chemical Reaction Systems with Toric Steady States

Mass-action chemical reaction systems are frequently used in computational biology. The corresponding polynomial dynamical systems are often large (consisting of tens or even hundreds of ordinary differential equations) and poorly parameterized (due to noisy measurement data and a small number of data points and repetitions). Therefore, it is often difficult to establish the existence of (positive) steady states or to determine whether more complicated phenomena such as multistationarity exist. If, however, the steady state ideal of the system is a binomial ideal, then we show that these questions can be answered easily. The focus of this work is on systems with this property, and we say that such systems have toric steady states. Our main result gives sufficient conditions for a chemical reaction system to have toric steady states. Furthermore, we analyze the capacity of such a system to exhibit positive steady states and multistationarity. Examples of systems with toric steady states include weakly-reversible zero-deficiency chemical reaction systems. An important application of our work concerns the networks that describe the multisite phosphorylation of a protein by a kinase/phosphatase pair in a sequential and distributive mechanism.     

Uncoupled Analysis of Stochastic Reaction Networks in Fluctuating Environments

The dynamics of stochastic reaction networks within cells are inevitably modulated by factors considered extrinsic to the network such as, for instance, the fluctuations in ribosome copy numbers for a gene regulatory network. While several recent studies demonstrate the importance of accounting for such extrinsic components, the resulting models are typically hard to analyze. In this work we develop a general mathematical framework that allows to uncouple the network from its dynamic environment by incorporating only the environment''s effect onto the network into a new model. More technically, we show how such fluctuating extrinsic components (e.g., chemical species) can be marginalized in order to obtain this decoupled model. We derive its corresponding process- and master equations and show how stochastic simulations can be performed. Using several case studies, we demonstrate the significance of the approach.     

Emergent Robustness in Competition Between Autocatalytic Chemical Networks

The origin of auto-catalytic networks has been proposed as an initial step in pre-biotic evolution. It is possible to derive simple models where auto-catalytic networks naturally arise from simple chemical mixtures. In order for such a system to develop, there needs to be some degree of stability, what is characterised as `robustness'. We demonstrate that competing systems generate this robustness as they create a distributed network of catalytic pathways.     

Hierarchical Feedback Modules and Reaction Hubs in Cell Signaling Networks

Despite much effort, identification of modular structures and study of their organizing and functional roles remain a formidable challenge in molecular systems biology, which, however, is essential in reaching a systematic understanding of large-scale cell regulation networks and hence gaining capacity of exerting effective interference to cell activity. Combining graph theoretic methods with available dynamics information, we successfully retrieved multiple feedback modules of three important signaling networks. These feedbacks are structurally arranged in a hierarchical way and dynamically produce layered temporal profiles of output signals. We found that global and local feedbacks act in very different ways and on distinct features of the information flow conveyed by signal transduction but work highly coordinately to implement specific biological functions. The redundancy embodied with multiple signal-relaying channels and feedback controls bestow great robustness and the reaction hubs seated at junctions of different paths announce their paramount importance through exquisite parameter management. The current investigation reveals intriguing general features of the organization of cell signaling networks and their relevance to biological function, which may find interesting applications in analysis, design and control of bio-networks.     

Oscillations in Biochemical Reaction Networks Arising from Pairs of Subnetworks

Biochemical reaction models show a variety of dynamical behaviors, such as stable steady states, multistability, and oscillations. Biochemical reaction networks with generalized mass action kinetics are represented as directed bipartite graphs with nodes for species and reactions. The bipartite graph of a biochemical reaction network usually contains at least one cycle, i.e., a sequence of nodes and directed edges which starts and ends at the same species node. Cycles can be positive or negative, and it has been shown that oscillations can arise as a result of either a positive cycle or a negative cycle. In earlier work it was shown that oscillations associated with a positive cycle can arise from subnetworks with an odd number of positive cycles. In this article we formulate a similar graph-theoretic condition, which generalizes the negative cycle condition for oscillations. This new graph-theoretic condition for oscillations involves pairs of subnetworks with an even number of positive cycles. An example of a calcium reaction network with generalized mass action kinetics is discussed in detail.