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.

What problem does this paper attempt to address?

This paper attempts to address the problem of the synergetic relationship between Chemical Reaction Networks (CRN) and Mathematical Epidemiology (ME). Specifically, the author proposes the following main research objectives: 1. **Define ME models as a subset of CRN models**: The author proposes to define mathematical epidemiology models as a subset of chemical reaction network models (see Open Problem 1). This definition helps to unify the theoretical frameworks of these two fields and reveals their intrinsic connections. 2. **Review and extend the Next - Generation Matrix method (NGM)**: The paper reviews the Next - Generation Matrix method (NGM) used in mathematical epidemiology to study the stability of boundary points. This method replaces the study of the Jacobian matrix through the theory of branching processes in probability theory. The author points out that although some conditions of the NGM theorem are not satisfied, in some cases, this method can still give the correct answer. In addition, the author emphasizes the lack of understanding of the validity conditions of the NGM method and proposes some heuristic methods (see Section 4.5). 3. **Explore the similarities between CRN models with Absolute Concentration Robustness (ACR) and simple ME models**: The author shows that simple CRN models with absolute concentration robustness are qualitatively very close to simple ME models. These models have a unique Disease - Free Equilibrium (DFE) and a unique internal equilibrium, and when \[ R_0:=\frac{s_d\beta}{\gamma_i + \mu}>1 \] the internal equilibrium enters the positive domain and becomes stable. This phenomenon is called the "relay effect": that is, when the DFE loses its stability, a new equilibrium enters the domain and takes over the dominant position. 4. **Develop and provide the Mathematica package Epid - CRN**: To promote the cooperation between researchers in these two fields, the author has developed a Mathematica package named Epid - CRN and provided several example notebooks. These tools can help researchers solve some small - scale open problems and provide computational support for more complex models. In conclusion, this paper aims to reveal the deep - seated connections between CRN and ME by integrating their methods and to provide new tools and perspectives for future research.

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

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

Finding bifurcations in mathematical epidemiology via reaction network methods

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

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

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

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

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

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

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

CRNs Exposed: Systematic Exploration of Chemical Reaction Networks

Predicting epidemics on weighted networks

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

The Critical Role of Networks to Describe Disease Spreading Dynamics in Social Systems: A Perspective

A HYBRID INTEGRO-DIFFERENTIAL EQUATION AND NETWORK BASED MODEL OF EPIDEMICS

Reactmine: a statistical search algorithm for inferring chemical reactions from time series data

bayesnec: An R Package for Concentration-Response Modeling and Estimation of Toxicity Metrics

Synthetic likelihood method for reaction network inference

A Coupled Spatial-Network Model: A Mathematical Framework for Applications in Epidemiology