Computing weakly reversible realizations of chemical reaction networks

Gabor Szederkenyi,Katalin M. Hangos,Zsolt Tuza

DOI: https://doi.org/10.48550/arXiv.1103.2892

2011-03-15

Dynamical Systems

Abstract:An algorithm is given in this paper for the computation of dynamically equivalent weakly reversible realizations with the maximal number of reactions, for chemical reaction networks (CRNs) with mass action kinetics.

What problem does this paper attempt to address?

Finding weakly reversible realizations of chemical reaction networks using optimization

Gabor Szederkenyi,Katalin M. Hangos,Zsolt Tuza

DOI: https://doi.org/10.48550/arXiv.1103.4741

2011-03-24

Abstract:An algorithm is given in this paper for the computation of dynamically equivalent weakly reversible realizations with the maximal number of reactions, for chemical reaction networks (CRNs) with mass action kinetics. The original problem statement can be traced back at least 30 years ago. The algorithm uses standard linear and mixed integer linear programming, and it is based on elementary graph theory and important former results on the dense realizations of CRNs. The proposed method is also capable of determining if no dynamically equivalent weakly reversible structure exists for a given reaction network with a previously fixed complex set.

Dynamical Systems,Molecular Networks
Computing weakly reversible linearly conjugate chemical reaction networks with minimal deficiency

Matthew D. Johnston,David Siegel,Gábor Szederkényi

DOI: https://doi.org/10.1016/j.mbs.2012.09.008

IF: 3.935

2013-01-01

Mathematical Biosciences

Abstract:Mass-action kinetics is frequently used in systems biology to model the behavior of interacting chemical species. Many important dynamical properties are known to hold for such systems if their underlying networks are weakly reversible and have a low deficiency. In particular, the Deficiency Zero and Deficiency One Theorems guarantee strong regularity with regards to the number and stability of positive equilibrium states. It is also known that chemical reaction networks with distinct reaction structure can admit mass-action systems with the same qualitative dynamics. The theory of linear conjugacy encapsulates the cases where this relationship is captured by a linear transformation. In this paper, we propose a mixed-integer linear programming algorithm capable of determining the minimal deficiency weakly reversible reaction network which admits a mass-action system which is linearly conjugate to a given reaction network.

biology,mathematical & computational biology
On the steady states of weakly reversible chemical reaction networks

Jian Deng,Christopher Jones,Martin Feinberg,Adrian Nachman

DOI: https://doi.org/10.48550/arXiv.1111.2386

2011-11-10

Quantitative Methods

Abstract:A natural condition on the structure of the underlying chemical reaction network, namely weak reversibility, is shown to guarantee the existence of an equilibrium (steady state) in each positive stoichiometric compatibility class for the associated mass-action system. Furthermore, an index formula is given for the set of equilibria in a given stoichiometric compatibility class.
Differentiable Programming of Chemical Reaction Networks

Alexander Mordvintsev,Ettore Randazzo,Eyvind Niklasson

DOI: https://doi.org/10.48550/arXiv.2302.02714

2023-02-06

Abstract:We present a differentiable formulation of abstract chemical reaction networks (CRNs) that can be trained to solve a variety of computational tasks. Chemical reaction networks are one of the most fundamental computational substrates used by nature. We study well-mixed single-chamber systems, as well as systems with multiple chambers separated by membranes, under mass-action kinetics. We demonstrate that differentiable optimisation, combined with proper regularisation, can discover non-trivial sparse reaction networks that can implement various sorts of oscillators and other chemical computing devices.

Molecular Networks,Machine Learning
Optimality-preserving Reduction of Chemical Reaction Networks

Kim G. Larsen,Daniele Toller,Mirco Tribastone,Max Tschaikowski,Andrea Vandin

DOI: https://doi.org/10.48550/arXiv.2301.08553

2023-01-20

Abstract:Across many disciplines, chemical reaction networks (CRNs) are an established population model defined as a system of coupled nonlinear ordinary differential equations. In many applications, for example, in systems biology and epidemiology, CRN parameters such as the kinetic reaction rates can be used as control inputs to steer the system toward a given target. Unfortunately, the resulting optimal control problem is nonlinear, therefore, computationally very challenging. We address this issue by introducing an optimality-preserving reduction algorithm for CRNs. The algorithm partitions the original state variables into a reduced set of macro-variables for which one can define a reduced optimal control problem from which one can exactly recover the solution of the original control problem. Notably, the reduction algorithm runs with polynomial time complexity in the size of the CRN. We use this result to reduce reachability and control problems of large-scale protein-interaction networks and vaccination models with hundreds of thousands of state variables.

Systems and Control
Programming and training rate-independent chemical reaction networks

Marko Vasić,Cameron Chalk,Austin Luchsinger,Sarfraz Khurshid,David Soloveichik

DOI: https://doi.org/10.1073/pnas.2111552119

2022-06-14

Abstract:Embedding computation in biochemical environments incompatible with traditional electronics is expected to have a wide-ranging impact in synthetic biology, medicine, nanofabrication, and other fields. Natural biochemical systems are typically modeled by chemical reaction networks (CRNs) which can also be used as a specification language for synthetic chemical computation. In this paper, we identify a syntactically checkable class of CRNs called noncompetitive (NC) whose equilibria are absolutely robust to reaction rates and kinetic rate law, because their behavior is captured solely by their stoichiometric structure. In spite of the inherently parallel nature of chemistry, the robustness property allows for programming as if each reaction applies sequentially. We also present a technique to program NC-CRNs using well-founded deep learning methods, showing a translation procedure from rectified linear unit (ReLU) neural networks to NC-CRNs. In the case of binary weight ReLU networks, our translation procedure is surprisingly tight in the sense that a single bimolecular reaction corresponds to a single ReLU node and vice versa. This compactness argues that neural networks may be a fitting paradigm for programming rate-independent chemical computation. As proof of principle, we demonstrate our scheme with numerical simulations of CRNs translated from neural networks trained on traditional machine learning datasets, as well as tasks better aligned with potential biological applications including virus detection and spatial pattern formation.
Uniqueness of weakly reversible and deficiency zero realizations of dynamical systems

Gheorghe Craciun,Jiaxin Jin,Polly Y. Yu

DOI: https://doi.org/10.1016/j.mbs.2021.108720

IF: 3.935

2021-12-01

Mathematical Biosciences

Abstract:A reaction network together with a choice of rate constants uniquely gives rise to a system of differential equations, according to the law of mass-action kinetics. On the other hand, different networks can generate the same dynamical system under mass-action kinetics. Therefore, the problem of identifying "the" underlying network of a dynamical system is not well-posed, in general. Here we show that the problem of identifying an underlying weakly reversible deficiency zero network is well-posed, in the sense that the solution is unique whenever it exists. This can be very useful in applications because from the perspective of both dynamics and network structure, a weakly reversible deficiency zero (WR<sub>0</sub>) realization is the simplest possible one. Moreover, while mass-action systems can exhibit practically any dynamical behavior, including multistability, oscillations, and chaos, WR<sub>0</sub> systems are remarkably stable for any choice of rate constants: they have a unique positive steady state within each invariant polyhedron, and cannot give rise to oscillations or chaotic dynamics. We also prove that both of our hypotheses (i.e., weak reversibility and deficiency zero) are necessary for uniqueness.

biology,mathematical & computational biology
CRNs Exposed: Systematic Exploration of Chemical Reaction Networks

Marko Vasic,David Soloveichik,Sarfraz Khurshid

DOI: https://doi.org/10.48550/arXiv.1912.06197

2019-12-12

Emerging Technologies

Abstract:Formal methods have enabled breakthroughs in many fields, such as in hardware verification, machine learning and biological systems. The key object of interest in systems biology, synthetic biology, and molecular programming is chemical reaction networks (CRNs) which formalizes coupled chemical reactions in a well-mixed solution. CRNs are pivotal for our understanding of biological regulatory and metabolic networks, as well as for programming engineered molecular behavior. Although it is clear that small CRNs are capable of complex dynamics and computational behavior, it remains difficult to explore the space of CRNs in search for desired functionality. We use Alloy, a tool for expressing structural constraints and behavior in software systems, to enumerate CRNs with declaratively specified properties. We show how this framework can enumerate CRNs with a variety of structural constraints including biologically motivated catalytic networks and metabolic networks, and seesaw networks motivated by DNA nanotechnology. We also use the framework to explore analog function computation in rate-independent CRNs. By computing the desired output value with stoichiometry rather than with reaction rates (in the sense that $X \to Y+Y$ computes multiplication by $2$), such CRNs are completely robust to the choice of reaction rates or rate law. We find the smallest CRNs computing the max, minmax, abs and ReLU (rectified linear unit) functions in a natural subclass of rate-independent CRNs where rate-independence follows from structural network properties.
Enhancing reliability of composed non-output-oblivious chemical reaction networks

Yu, Sihai,Mao, Qianfei

DOI: https://doi.org/10.1007/s11047-024-09999-w

2024-10-22

Natural Computing

Abstract:A chemical reaction network (CRN) is a molecular computing paradigm which aims to exploit the dynamic reactions among various chemical substances as computational resources. The well-mixed nature of chemical molecules, however, tends to obstruct the connection of multiple networks into a larger and more complex network. Such poor modularity is mainly attributed to the non-output-oblivious characteristic of most CRNs because their outputs can only be stabilized after all reactions have completed, whereby feeding the occasional outputs of one CRN before the CRN stabilizes to other connected CRN will cause unexpected effects and errors on both networks. In this paper, we propose an extended clock module to intermediate the connection between two non-output-oblivious CRNs. Especially, the clock yields a reasonable delay in-between the two CRNs in a natural way, which attempts to stabilize the former network before its outputs can be processed by the connected network. Numerical experiments show the clock can significantly reduce the occurrence of errors when composing two CRNs to conduct more complex tasks, thereby improving the modularity of non-output-oblivious CRNs. To this end, we illustrate how to use this method to combine the chemical reaction decider and the chemical reaction computer for predicate and function computations. Finally, a synchronization-waiting mechanism is explored for composing multiple non-output-oblivious CRNs.

computer science, artificial intelligence, theory & methods
Recurrent neural chemical reaction networks that approximate arbitrary dynamics

Alexander Dack,Benjamin Qureshi,Thomas E. Ouldridge,Tomislav Plesa

2024-06-06

Abstract:Many important phenomena in chemistry and biology are realized via dynamical features such as multi-stability, oscillations, and chaos. Construction of novel chemical systems with such finely-tuned dynamics is a challenging problem central to the growing field of synthetic biology. In this paper, we address this problem by putting forward a molecular version of a recurrent artificial neural network, which we call a recurrent neural chemical reaction network (RNCRN). We prove that the RNCRN, with sufficiently many auxiliary chemical species and suitable fast reactions, can be systematically trained to achieve any dynamics. This approximation ability is shown to hold independent of the initial conditions for the auxiliary species, making the RNCRN more experimentally feasible. To demonstrate the results, we present a number of relatively simple RNCRNs trained to display a variety of biologically-important dynamical features.

Molecular Networks,Dynamical Systems
Source-only realizations, weakly reversible deficiency one networks and dynamical equivalence

Abhishek Deshpande

DOI: https://doi.org/10.48550/arXiv.2205.00801

2022-05-02

Abstract:Reaction networks can display a wide array of dynamics. However, it is possible for different reaction networks to display the same dynamics. This phenomenon is called dynamical equivalence and makes network identification a hard problem to solve. We show that to find a strongly endotactic, endotactic, consistent, or conservative realization that is dynamically equivalent to the mass-action system generated by a given network, it suffices to consider only the source vertices of the given network. In addition, we show that weakly reversible deficiency one realizations are not unique. We also present a characterization of the dynamical relationships that exist between several types of weakly reversible deficiency one networks.

Dynamical Systems
Stochastic Chemical Reaction Networks with Discontinuous Limits and AIMD processes

Lucie Laurence,Philippe Robert

2024-06-18

Abstract:In this paper we study a class of stochastic chemical reaction networks (CRNs) for which chemical species are created by a sequence of chain reactions. We prove that under some convenient conditions on the initial state, some of these networks exhibit a discrete-induced transitions (DIT) property: isolated, random, events have a direct impact on the macroscopic state of the process. If this phenomenon has already been noticed in several CRNs, in auto-catalytic networks in the literature of physics in particular, there are up to now few rigorous studies in this domain. A scaling analysis of several cases of such CRNs with several classes of initial states is achieved. The DIT property is investigated for the case of a CRN with four nodes. We show that on the normal timescale and for a subset of (large) initial states and for convenient Skorohod topologies, the scaled process converges in distribution to a Markov process with jumps, an Additive Increase/Multiplicative Decrease (AIMD) process. This asymptotically discontinuous limiting behavior is a consequence of a DIT property due to random, local, blowups of jumps occurring during small time intervals. With an explicit representation of invariant measures of AIMD processes and time-change arguments, we show that, with a speed-up of the timescale, the scaled process is converging in distribution to a continuous deterministic function. The DIT analyzed in this paper is connected to a simple chain reaction between three chemical species and is therefore likely to be a quite generic phenomenon for a large class of CRNs.

Probability
Asymptotology of Chemical Reaction Networks

A. N. Gorban,O. Radulescu,A. Y. Zinovyev

DOI: https://doi.org/10.1016/j.ces.2009.09.005

2010-06-15

Abstract:The concept of the limiting step is extended to the asymptotology of multiscale reaction networks. Complete theory for linear networks with well separated reaction rate constants is developed. We present algorithms for explicit approximations of eigenvalues and eigenvectors of kinetic matrix. Accuracy of estimates is proven. Performance of the algorithms is demonstrated on simple examples. Application of algorithms to nonlinear systems is discussed.

Chemical Physics,Computational Physics
Reachability Problems for Continuous Chemical Reaction Networks

Adam Case,Jack H. Lutz,D. M. Stull

DOI: https://doi.org/10.48550/arXiv.1508.04125

2015-09-03

Abstract:Chemical reaction networks (CRNs) model the behavior of molecules in a well-mixed system. The emerging field of molecular programming uses CRNs not only as a descriptive tool, but as a programming language for chemical computation. Recently, Chen, Doty and Soloveichik introduced a new model of chemical kinetics, rate-independent continuous CRNs (CCRNs), to study the chemical computation of continuous functions. A fundamental question of a CRN is whether a state of the system is reachable through a sequence of reactions in the network. This is known as the reachability problem. In this paper, we investigate CCRN-REACH, the reachability problem for this model of chemical reaction networks. We show that, for continuous CRNs, constructing a path to a state of the network is computable in polynomial time. We also prove that a related problem, Sub-CCRN-REACH, is NP-complete.

Computational Complexity
Universal Slow Dynamics of Chemical Reaction Networks

Masanari Shimada,Pegah Behrad,Eric De Giuli

DOI: https://doi.org/10.1103/PhysRevE.109.044105

2024-03-14

Abstract:Understanding the emergent behavior of chemical reaction networks (CRNs) is a fundamental aspect of biology and its origin from inanimate matter. A closed CRN monotonically tends to thermal equilibrium, but when it is opened to external reservoirs, a range of behaviors is possible, including transition to a new equilibrium state, a non-equilibrium state, or indefinite growth. This study shows that slowly driven CRNs are governed by the conserved quantities of the closed system, which are generally far fewer in number than the species. Considering both deterministic and stochastic dynamics, a universal slow dynamics equation is derived with singular perturbation methods, and is shown to be thermodynamically consistent. The slow dynamics is highly robust against microscopic details of the network, which may be unknown in practical situations. In particular, non-equilibrium states of realistic large CRNs can be sought without knowledge of bulk reaction rates. The framework is successfully tested against a suite of networks of increasing complexity and argued to be relevant in the treatment of open CRNs as chemical machines.

Molecular Networks,Statistical Mechanics,Adaptation and Self-Organizing Systems
Weakly reversible deficiency one realizations of polynomial dynamical systems: an algorithmic perspective

Gheorghe Craciun,Abhishek Deshpande,Jiaxin Jin

2023-03-17

Abstract:Given a dynamical system with polynomial right-hand side, can it be generated by a reaction network that possesses certain properties? This question is important because some network properties may guarantee specific dynamical properties, such as existence or uniqueness of equilibria, persistence, permanence, or global stability. Here we focus on this problem in the context of weakly reversible deficiency one networks. In particular, we describe an algorithm for deciding if a polynomial dynamical system admits a weakly reversible deficiency one realization, and identifying one if it does exist. In addition, we show that weakly reversible deficiency one realizations can be partitioned into mutually exclusive Type I and Type II realizations, where Type I realizations guarantee existence and uniqueness of positive steady states, while Type II realizations are related to stoichiometric generators, and therefore to multistability.

Dynamical Systems
Computing Threshold Circuits with Void Reactions in Step Chemical Reaction Networks

Rachel Anderson,Alberto Avila,Bin Fu,Timothy Gomez,Elise Grizzell,Aiden Massie,Gourab Mukhopadhyay,Adrian Salinas,Robert Schweller,Evan Tomai,Tim Wylie

2024-02-13

Abstract:We introduce a new model of \emph{step} Chemical Reaction Networks (step CRNs), motivated by the step-wise addition of materials in standard lab procedures. Step CRNs have ordered reactants that transform into products via reaction rules over a series of steps. We study an important subset of weak reaction rules, \emph{void} rules, in which chemical species may only be deleted but never changed. We demonstrate the capabilities of these simple limited systems to simulate threshold circuits and compute functions using various configurations of rule sizes and step constructions, and prove that without steps, void rules are incapable of these computations, which further motivates the step model. Additionally, we prove the coNP-completeness of verifying if a given step CRN computes a function, holding even for $O(1)$ step systems.

Molecular Networks,Emerging Technologies
Translated Chemical Reaction Networks

Matthew D. Johnston

DOI: https://doi.org/10.1007/s11538-014-9947-5

2013-06-12

Abstract: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 network's 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.

Dynamical Systems,Molecular Networks
Deficiency, kinetic invertibility, and catalysis in stochastic chemical reaction networks

Shesha Gopal Marehalli Srinivas,Matteo Polettini,Massimiliano Esposito,Francesco Avanzini

DOI: https://doi.org/10.1063/5.0147283

IF: 4.304

2023-05-22

The Journal of Chemical Physics

Abstract:Stochastic chemical processes are described by the chemical master equation satisfying the law of mass-action. We first ask whether the dual master equation, which has the same steady state as the chemical master equation, but with inverted reaction currents, satisfies the law of mass-action and, hence, still describes a chemical process. We prove that the answer depends on the topological property of the underlying chemical reaction network known as deficiency. The answer is yes only for deficiency-zero networks. It is no for all other networks, implying that their steady-state currents cannot be inverted by controlling the kinetic constants of the reactions. Hence, the network deficiency imposes a form of non-invertibility to the chemical dynamics. We then ask whether catalytic chemical networks are deficiency-zero. We prove that the answer is no when they are driven out of equilibrium due to the exchange of some species with the environment.

chemistry, physical,physics, atomic, molecular & chemical

Computing weakly reversible realizations of chemical reaction networks

Finding weakly reversible realizations of chemical reaction networks using optimization

Computing weakly reversible linearly conjugate chemical reaction networks with minimal deficiency

On the steady states of weakly reversible chemical reaction networks

Differentiable Programming of Chemical Reaction Networks

Optimality-preserving Reduction of Chemical Reaction Networks

Programming and training rate-independent chemical reaction networks

Uniqueness of weakly reversible and deficiency zero realizations of dynamical systems

CRNs Exposed: Systematic Exploration of Chemical Reaction Networks

Enhancing reliability of composed non-output-oblivious chemical reaction networks

Recurrent neural chemical reaction networks that approximate arbitrary dynamics

Source-only realizations, weakly reversible deficiency one networks and dynamical equivalence

Stochastic Chemical Reaction Networks with Discontinuous Limits and AIMD processes

Asymptotology of Chemical Reaction Networks

Reachability Problems for Continuous Chemical Reaction Networks

Universal Slow Dynamics of Chemical Reaction Networks

Weakly reversible deficiency one realizations of polynomial dynamical systems: an algorithmic perspective

Computing Threshold Circuits with Void Reactions in Step Chemical Reaction Networks

Translated Chemical Reaction Networks

Deficiency, kinetic invertibility, and catalysis in stochastic chemical reaction networks