Abstract:Complex systems can be advantageously modeled by formal reaction systems (RS), a.k.a. chemical reaction networks in chemistry. Reaction-based models can indeed be interpreted in a hierarchy of semantics, depending on the question at hand, most notably by Ordinary Differential Equations (ODEs), Continuous Time Markov Chains (CTMCs), discrete Petri nets and asynchronous Boolean transition systems. The last three semantics can be easily related in the framework of abstract interpretation. The first two are classically related by Kurtz's limit theorem which states that if reactions are density-dependent families, then, as the volume goes to infinity, the mean reactant concentrations of the CTMC tends towards the solution of the ODE. In the more realistic context of bounded volumes, it is easy to show, by moment closure, that the restriction to reactions with at most one reactant ensures similarly that the mean of the CTMC trajectories is equal to the solution of the ODE at all time points. In this paper, we generalize that result in presence of polyreactant reactions, by introducing the Stoichiometric Influence and Modification Graph (SIMG) of an RS, and by showing that the equality between the two interpretations holds for the variables that belong to distinct SIMG ancestors of polyreactant reactions. We illustrate this approach with several examples. Evaluation on BioModels reveals that the condition for all variables is satisfied on models with no polymolecular reaction only. However, our theorem can be applied selectively to certain variables of the model to provide insights into their behaviour within more complex systems. Interestingly, we also show that the equality holds for a basic oscillatory RS implementing the sine and cosine functions of time.

Mathematical Analysis of Chemical Reaction Systems

Mathematical Formalism of Nonequilibrium Thermodynamics for Nonlinear Chemical Reaction Systems with General Rate Law

Translated Chemical Reaction Networks

Discrete models for chemically reacting systems

Probabilistic and Maximum Entropy Modeling of Chemical Reaction Systems: Characteristics and Comparisons to Mass Action Kinetic Models

Graphical Conditions ensuring Equality between Differential and Mean Stochastic Dynamics

Stochastic analysis of chemical reactions in multi-component interacting systems at criticality

Explicit exponential convergence to equilibrium for nonlinear reaction-diffusion systems with detailed balance condition

Mathematical Analysis of Some Reaction Networks Inducing Biological Growth/Decay Functions.

Multiscale Analysis of Reaction Networks

Quasichemical Models of Multicomponent Nonlinear Diffusion

A Graphic Formulation of Nonisothermal Chemical Reaction Systems and the Analysis of Detailed Balanced Networks.

Stability analysis of irreversible chemical reaction-diffusion systems with boundary equilibria

Positive Steady-State Varieties of Small Chemical Reaction Networks

Analysis of Steady State Multiplicity of Chemical Reaction Systems

Fixed Points for Stochastic Open Chemical Systems

Model reduction in chemical dynamics: slow invariant manifolds, singular perturbations, thermodynamic estimates, and analysis of reaction graph

Convergence to equilibrium of renormalised solutions to nonlinear chemical reaction-diffusion systems

Comprehensive Numerical Analysis of Time-Fractional Reaction–Diffusion Models with Applications to Chemical and Biological Phenomena

Sensitivity of chemical reaction networks: a structural approach 3. Regular multimolecular systems

Modeling of chemical reaction systems with detailed balance using gradient structures