Abstract:In this essay, we investigate some relations between Chemical Reaction Networks (CRN) and Mathematical Epidemiology (ME) and report on several pleasant surprises which we had simply by putting these two topics together. Firstly, we propose a definition of ME models as a subset of CRN models. Secondly, we review a fundamental stability result for boundary points, known in ME as the NGM method since it replaces the investigation of the Jacobian by that of a matrix whose origins lie in probability (the theory of branching processes). This important result seems to be little known outside of ME; even in ME, it has not been made clear before that the method gets sometimes the right answer, even though the conditions of the NGM theorem are not all satisfied. Thus, beyond the theorem, there is a heuristic approach, the validity conditions for which are not sufficiently understood. Thirdly, we show that some simple CRN models with absolute concentration robustness (ACR), are close qualitatively to simple ME models, in the sense that they have an unique disease free equilibrium, and a unique interior fixed point, and the latter enters the positive domain and becomes stable precisely when $R_0:=s_{dfe} \mathcal{R}=\frac{s_{dfe}}{s_e}>1.$ (where $s$ denotes the "ACR species"). Thus, for these "ME type models", a "relay phenomena" takes place: precisely when the DFE loses stability, a new fixed point enters the domain, and takes over. Last but not least, we offer in the associated GitHub repository <a class="link-external link-https" href="https://github.com/adhalanay/epidemiology_crns" rel="external noopener nofollow">this https URL</a> a Mathematica package, Epid-CRN, which is addressed to researchers of both disciplines, and provide illustrative notebooks, which in particular solve a few minor open ME and CRN problems. Our package may also be used to study easy cases of analogue continuous time Markov chain (CTMC) ME and CRN models.

Advancing Mathematical Epidemiology and Chemical Reaction Network Theory via Synergies Between Them

Advancing Mathematical Epidemiology and Chemical Reaction Network Theory via Synergies Between Them

Investigating synergies between chemical reaction networks (CRN) and mathematical epidemiology (ME), using the Mathematica package Epid-CRN

Advancing Mathematical Epidemic Modeling via synergies with Chemical Reaction Network Theory and Lagrange-Hamilton Geometry

Finding bifurcations in mathematical epidemiology via reaction network methods

From Genome-Scale Data to Models of Infectious Disease: A Bayesian Network-Based Strategy to Drive Model Development

CERENA: ChEmical REaction Network Analyzer—A Toolbox for the Simulation and Analysis of Stochastic Chemical Kinetics

A Scalable Moment-Closure Approximation for Large-Scale Biochemical Reaction Networks

Modeling of Stochastic Biological Processes with Non-Polynomial Propensities Using Non-Central Conditional Moment Equation

EpiModel: An R Package for Mathematical Modeling of Infectious Disease over Networks

Bayesian Parameter Estimation for Biochemical Reaction Networks Using Region-Based Adaptive Parallel Tempering

Epidemionics: from the host-host interactions to the systematic analysis of the emergent macroscopic dynamics of epidemic networks

Scalable Parameter Estimation for Genome-Scale Biochemical Reaction Networks

The two-step exponential decay reaction network: analysis of the solutions and relation to epidemiological SIR models with logistic and Gompertz type infection contact patterns

Poincaré and SimBio: a versatile and extensible Python ecosystem for modeling systems

Scaling methods for accelerating kinetic Monte Carlo simulations of chemical reaction networks

Differential and graphical approaches to multistability theory for chemical reaction networks

Polyhedral geometry and combinatorics of an autocatalytic ecosystem

A Rule-Based Epidemiological Modelling Framework

Performing global sensitivity analysis on simulations of a continuous-time Markov chain model motivated by epidemiology