No Arabic abstract
We study how the properties of allowing multiple positive nondegenerate equilibria (MPNE) and multiple positive linearly stable equilibria (MPSE) are inherited in chemical reaction networks (CRNs). Specifically, when is it that we can deduce that a CRN admits MPNE or MPSE based on analysis of its subnetworks? Using basic techniques from analysis we are able to identify a number of situations where MPNE and MPSE are inherited as we build up a network. Some of these modifications are known while others are new, but all results are proved using the same basic framework, which we believe will yield further results. The results are presented primarily for mass action kinetics, although with natural, and in some cases immediate, generalisation to other classes of kinetics.
The goal of this paper is to gather and develop some necessary and sufficient criteria for injectivity and multistationarity in vector fields associated with a chemical reaction network under a variety of more or less general assumptions on the nature of the network and the reaction rates. The results are primarily linear algebraic or matrix-theoretic, with some graph-theoretic results also mentioned. Several results appear in, or are close to, results in the literature. Here, we emphasise the connections between the results, and where possible, present elementary proofs which rely solely on basic linear algebra and calculus. A number of examples are provided to illustrate the variety of subtly different conclusions which can be reached via different computations. In addition, many of the computations are implemented in a web-based open source platform, allowing the reader to test examples including and beyond those analysed in the paper.
Many biochemical and industrial applications involve complicated networks of simultaneously occurring chemical reactions. Under the assumption of mass action kinetics, the dynamics of these chemical reaction networks are governed by systems of polynomial ordinary differential equations. The steady states of these mass action systems have been analysed via a variety of techniques, including elementary flux mode analysis, algebraic techniques (e.g. Groebner bases), and deficiency theory. In this paper, we present a novel method for characterizing the steady states of mass action systems. Our method explicitly links a networks capacity to permit a particular class of steady states, called toric steady states, to topological properties of a related network called a translated chemical reaction network. These networks share their reaction stoichiometries with their source network but are permitted to have different complex stoichiometries and different network topologies. We apply the results to examples drawn from the biochemical literature.
A chemical reaction network (CRN) is composed of reactions that can be seen as interactions among entities called species, which exist within the system. Endowed with kinetics, CRN has a corresponding set of ordinary differential equations (ODEs). In Chemical Reaction Network Theory, we are interested with connections between the structure of the CRN and qualitative properties of the corresponding ODEs. One of the results in Decomposition Theory of CRNs is that the intersection of the sets of positive steady states of the subsystems is equal to the set of positive steady states of the whole system, if the decomposition is independent. Hence, computational approach using independent decompositions can be used as an efficient tool in studying large systems. In this work, we provide a necessary and sufficient condition for the existence of a nontrivial independent decomposition of a CRN, which leads to a novel step-by-step method to obtain such decomposition, if it exists. We also illustrate these results using real-life examples. In particular, we show that a CRN of a popular model of anaerobic yeast fermentation pathway has a nontrivial independent decomposition, while a particular biological system, which is a metabolic network with one positive feedforward and a negative feedback has none. Finally, we analyze properties of steady states of reaction networks of specific influenza virus models.
We are concerned with polynomial ordinary differential systems that arise from modelling chemical reaction networks. For such systems, which may be of high dimension and may depend on many parameters, it is frequently of interest to obtain a reduction of dimension in certain parameter ranges. Singular perturbation theory, as initiated by Tikhonov and Fenichel, provides a path toward such reductions. In the present paper we discuss parameter values that lead to singular perturbation reductions (so-called Tikhonov-Fenichel parameter values, or TFPVs). An algorithmic approach is known, but it is feasible for small dimensions only. Here we characterize conditions for classes of reaction networks for which TFPVs arise by turning off reactions (by setting rate parameters to zero), or by removing certain species (which relates to the classical quasi-steady state approach to model reduction). In particular, we obtain definitive results for the class of complex balanced reaction networks (of deficiency zero) and first order reaction networks.
We study chemical reaction networks with discrete state spaces, such as the standard continuous time Markov chain model, and present sufficient conditions on the structure of the network that guarantee the system exhibits an extinction event. The conditions we derive involve creating a modified chemical reaction network called a domination-expanded reaction network and then checking properties of this network. We apply the results to several networks including an EnvZ-OmpR signaling pathway in Escherichia coli. This analysis produces a system of equalities and inequalities which, in contrast to previous results on extinction events, allows algorithmic implementation. Such an implementation will be investigated in a companion paper where the results are applied to 458 models from the European Bioinformatics Institutes BioModels database.