Next Article in Journal
Biochemical and Biophysical Characterization of the dsDNA Packaging Motor from the Lactococcus lactis Bacteriophage Asccphi28
Next Article in Special Issue
viromeBrowser: A Shiny App for Browsing Virome Sequencing Analysis Results
Previous Article in Journal
Beyond the Surface: Endocytosis of Mosquito-Borne Flaviviruses
Previous Article in Special Issue
Characterization of a Novel Mitovirus of the Sand Fly Lutzomyia longipalpis Using Genomic and Virus–Host Interaction Signatures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Structure and Hierarchy of SARS-CoV-2 Infection Dynamics Models Revealed by Reaction Network Analysis

1
Department of Fundamental Sciences, Ernst-Abbe University of Applied Sciences Jena, Carl-Zeiss-Promenade 2, 07745 Jena, Germany
2
Bio Systems Analysis Group, Department of Mathematics and Computer Science, University of Jena, Ernst-Abbe-Platz 2, 07743 Jena, Germany
3
Department of Mathematics and Natural Sciences, Centre for Applied Mathematics and Bioinformatics, Gulf University for Science and Technology, 32093 Hawally, Kuwait
4
European Virus Bioinformatics Center, Leutragraben 1, 07743 Jena, Germany
*
Authors to whom correspondence should be addressed.
Submission received: 2 November 2020 / Revised: 8 December 2020 / Accepted: 16 December 2020 / Published: 23 December 2020
(This article belongs to the Special Issue Virus Bioinformatics 2020)

Abstract

:
This work provides a mathematical technique for analyzing and comparing infection dynamics models with respect to their potential long-term behavior, resulting in a hierarchy integrating all models. We apply our technique to coupled ordinary and partial differential equation models of SARS-CoV-2 infection dynamics operating on different scales, that is, within a single organism and between several hosts. The structure of a model is assessed by the theory of chemical organizations, not requiring quantitative kinetic information. We present the Hasse diagrams of organizations for the twelve virus models analyzed within this study. For comparing models, each organization is characterized by the types of species it contains. For this, each species is mapped to one out of four types, representing uninfected, infected, immune system, and bacterial species, respectively. Subsequently, we can integrate these results with those of our former work on Influenza-A virus resulting in a single joint hierarchy of 24 models. It appears that the SARS-CoV-2 models are simpler with respect to their long term behavior and thus display a simpler hierarchy with little dependencies compared to the Influenza-A models. Our results can support further development towards more complex SARS-CoV-2 models targeting the higher levels of the hierarchy.

1. Introduction

The current SARS-CoV-2 pandemic has required huge efforts from global society and the scientific community to track, understand, and combat its proliferation. Models of the infection dynamics can help understand SARS-CoV-2 pathogenesis, develop optimal treatments, and introduce appropriate measures to prevent the spread of the virus. There are a multitude of modeling approaches with different properties, applications and aims that can be classed into categories of in-host models (e.g., [1,2,3,4,5,6,7,8]) versus host-to-host models (such as [9,10,11,12]), discrete versus continuous models and ODE versus PDE models (for an overview we refer to [13,14,15,16]). There is an accumulating body of literature on SARS-CoV-2 infection dynamics that make use of these various tools and provide datasets that can be analyzed retrospectively once consensus modeling strategies have been derived [17,18]. The aforementioned models have in common that they rely on an identifiable reaction network, for instance, a set of species and a set of reactions that describe the possible interactions of these species. We have shown that for Influenza A virus infection dynamics [19] reaction network analysis (especially COT [20]) provides metrics to understand, analyze, and categorize different in-host ODE models. In our current work, we applied this to SARS-CoV-2 infection dynamics by extending our previous approach [19] in several directions: We incorporate in-host, host-to-host, and linked models consisting either of ODEs or PDEs. Finally, we combine the models of SARS-CoV-2 with Influenza A in order to compare the dynamics for both viruses. In the next chapter (Materials and Methods) we introduce the herein used method of reaction network analysis by applying it to an example model from Almocera [1]. In the Results chapter, we describe the structural analysis of a set of representative models. For each of these models, we then derive what we call signature. The signature of a model gives a brief overview of its potential dynamical behavior, which allows for relating several models to each other and combining them in a hierarchy. Finally, we link the respective hierarchy of SARS-CoV-2 with that from Influenza A from [19]. In the Conclusions, we summarize our findings and the benefits of the technique applied herein. This novel method can be used as an instrument for a deeper understanding of infection dynamics models and further for an appropriate construction of future virus infection models.
All models, as well as the software tools used to do the analysis, can be found on Github (https://github.com/stephanpeter/orgs-covid).

2. Materials and Methods: Procedure for the Organizational Analysis

Our method is an extension of that used in [19] to analyze Influenza A virus modeling. We will briefly describe our method for the example of the in-host ODE model from Almocera [1] (see Figure 1).
The method we apply in this work aims at deriving the signature of each model we analyze. Therefore, for each model, the set of organizations must be computed. The computation of the organizations can be divided into two steps:
  • Step 1—Deriving the set of reactions: Each summand of each ODE (or PDE) is translated into a reaction as illustrated by the transition from Subfigures (a) to (b) in Figure 1. On the left-hand side of each reaction formula, there is a set of species, the so-called support of a reaction. The support of a reaction is the unique set of species that are needed to run the reaction. If only one of the species of the support of the reaction is missing then that reaction is not active. The term (of the ODE (or PDE)) that belongs to that reaction must be zero if and only if the concentration of at least one of the species in the support of that reaction is zero. The number of the appearance of each species of a reaction on the right-hand side of a reaction is bigger or less than the number on the left-hand side depending on whether the regarding term has a positive or negative sign in the ODE (or PDE)) of the regarding species. As an example we consider reaction R 3 . The corresponding summand is c V E V . It is zero if and only if the concentration of at least one of E or V is zero. Thus the support of R 3 contains exactly the species E and V. On the left-hand side of the reaction equation of R 3 the species E resp. V appear only to the power of one because of the power of E resp. V is one in c V E V . Since the summand c V E V appears only in the ODE of V, namely with a negative sign, the right-hand side of the reaction equation of R 3 contains one less of V than the left-hand side. The number of E is equal on both sides of the reaction equation since the amount of E is not affected by the reaction R 3 .
  • Step 2—Calculating the organizations from the set of reactions: The second step is to compute the organizations (as defined in [21]) from the derived reactions. Each organization consists of a subset of species that is
    • closed and
    • self-maintaining.
Before explaining closedness and self-maintenance we have to introduce some other terms regarding reaction networks. Firstly, we discuss reaction equations shortly. As an example we take the reaction
R 6 : E + V r 2 V .
As for every reaction, it contains a reaction arrow, that separates its left-hand side from its right-hand side. The left-hand side, which we call the support of R 6 (or shortly s u p p ( R 6 ) ), includes one entity of each species E resp. V linked by a plus. The right-hand side of the reaction equation of R 6 contains two entities of species E.
Now we explain the stoichiometric matrix N of a given reaction network, i.e., of a set of species together with a set of reactions. The stoichiometric matrix N Z n × m consists of n lines and m columns. Thereby, n is the number of species of the reaction network and m the number of reactions. e.g., the sixth column of the reaction network of the example Almocera In-host Model can be determined from the reaction equation of R 6 . For each of the two species, the difference of its number of occurrence on the right-hand side minus the number of occurrences on the left-hand side must be calculated. For viruses V we have no appearance on the right-hand side of the reaction equation of R 6 and one appearance on the left-hand side. Thus the element in the first line and sixth column of N equals 0 1 = 1 . Similarly for E we get for the second line and sixth column of N the value 2 1 = 1 .
For a given subset S S of species, a vector v R + m of m non-negative real numbers v r from
R + { x R : x 0 } ,
one for each reaction r R , is called feasible flux if and only if for all reactions r R
v r > 0 s u p p ( r ) S .
We are now able to explain what it means for a subset of species to be closed resp. self-maintaining. These concepts of closedness and self-maintenance stem from chemical organization theory (COT) [23,24].
  • A subset S S of species is closed if and only if for each reaction with its support contained in that subset, also all species appearing on the right-hand side of the reaction equation are contained in that subset. In other words, no reaction that is active on the subset S produces a species that is not contained in that subset.
  • A subset S S of species is self-maintaining if and only if there is a feasible flux for S v R + m for S such that
    N · v 0 R + n .
All organizations of a given reaction network can be arranged in a so-called Hasse diagram. For the In-host Almocera example model the Hasse diagram is shown in Figure 1c. From the bottom to the top the organizations have increasing size, indicating an increasing number of species. A line is drawn between two organizations if and only if one is a subset of the other and there is no organization between them. Thus, there is a line between the organizations { E } and { V , E } .
In [21] it was proven that for appropriate kinetics, for example, mass-action kinetics, every fixed-point of an ODE is represented by exactly one organization. More precisely, the set of species, that have a strictly positive concentration in the fixed-point, is an organization. Later on, this result was generalized from fixed-points to all kinds of attractors of such dynamical systems [22]. This means that the Hasse diagram of organizations gives an overview of all possible attractors of a dynamical system on the abstract level of the underlying reaction network, regardless of the parameters like reaction constants that were used for the ODEs (or PDEs). Furthermore, COT allows for statements about transitions between those attractors. See Figure 2 for a summary of the fore-mentioned explanations.
Finally, different models with different reaction networks can be analyzed, compared and related to each other in a hierarchy as it was done for Influenza-A virus (IAV) infections [19] and also incorporating the PDEs analysis [25].
In this work, we apply the technique described above to models of SARS-CoV-2 infection dynamics and extend it by including not only dynamical systems
  • consisting of ODEs but also of PDEs,
  • describing in-host dynamics but also host-to-host and mixed (in-host and host-to-host) models.
  • Analyzed infection dynamics of SARS-CoV-2 but also compared to Influenza models
We also show that it is possible to analyze the interrelations between models of different application cases, in this case by ordering Influenza A, SARS-CoV-2, and general virus infection models together in one single hierarchy.

3. Analysis of the Models

In this chapter we provide the analysis of all the models listed in Table 1. We provide the respective ODEs and PDEs and the Hasse diagrams of organizations. The reactions that we derive from the models are listed in the Appendix A. We will start with the models basing on ODEs. These models, be they in-host or host-to-host models, were constructed especially for modeling the SARS-CoV-2 infection dynamics. After having analyzed the ODE models, we did the same for the PDE models. Note, that the latter were not solely built to model SARS-CoV-2 infection but rather viral infection dynamics in general. Finally, we analyze one ODE model (the Almocera model [1,2]) linking an in-host scenario (the example model from the Introduction) with a host-to-host scenario of virus infection dynamics.

3.1. In-Host Models

Here we firstly present four in-host ODE models with increasing numbers of species (see Figure 3, Figure 4, Figure 5 and Figure 6). These models describe the spread of the infection within a host, in this case humans. All models contain a virus species but the models differ in terms of the identity of the species.
Now we analyze two in-host PDE models (see Figure 7 and Figure 8). Contrary to ODE models they are able to deal with spatial inhomogeneities of viral infection processes in the host. They were designed for general viral infections. Thus it is recommended for future SARS-CoV-2 infection dynamics modeling to adapt these models to capture the specifics of this new virus.
Within the set of in-host models we observe a principal difference between the Abuin Model and the other models: The Abuin model does not have an organization with regard to the viral species causing the infection. Thus the Abuin model implicitly assumes the vanishing of the infection over time. The other models do not share this property and thus contain no assumptions regarding viral persistence, which may confer an advantage to these models since it is unclear what the extent of SARS-CoV-2 persistence is. The Su Model exhibits the most complex Hasse diagram. This is the model with the biggest number of species and the only model that explicitly focuses on the genetic aspects of SARS-CoV-2 infection dynamics. Interestingly the Su Model has a distributed organization that is not an organization. This emphasizes the role of the distribution of the species in space or time which is the subject of current research which in turn was recently initiated [25]. Since the Su Model only has ODEs that do not allow for modeling spatial inhomogeneities this is an indication that adapting this model to PDEs may improve the model quality.

3.2. Host-To-Host Models

In this section we first analyze three different host-to-host ODE models describing SARS-CoV-2 infection as it spreads in a human population from one individual to the next (see Figure 9, Figure 10 and Figure 11). Thus these models have three host species in common: susceptible, uninfected individuals (S) and infected individuals (I).
The host-to-host PDE models we subsequently analyzed (see Figure 12 and Figure 13) have thus far only been applied to general viral infections. Because of the importance of the spatial dimension in SARS-CoV-2 transmissions, through interventions such as social distancing, it is pertinent to apply this approach to the current outbreak.
The Wu Model and the Fitzgibbon-I Model are the only ones that cover the infection dynamics on the level of organizations. However the Nesteruk Model, the Bai Model, and the Fitgibbon-II Model have solely organizations without species representing infection. Thus these latter (three) models implicitly assume a vanishing of the infection in the long-term. It is currently unclear whether this assumption is justified for SARS-CoV-2.

3.3. A Linked In-Host/Host-To-Host Model

Here we analyze a model that we called “linked model” (see Figure 14) as it includes in-host as well as host-to-host dynamics, both described by ODEs. This model is designed for viral infections in general and its application to SARS-CoV-2 was deemed of interest because of its bigger focus compared to the previous models. For the analysis of the in-host part of the model we refer the reader to Figure 1 in the Materials and Methods Section 2. In the following we analyzed the linked model, where the in-host model is incorporated into a host-to-host model.

3.4. Hierarchy of Models

In this section, we present the hierarchies of all models. First, we show the hierarchy of SARS-CoV-2 models (see Figure 15) and then the merged hierarchy of SARS-CoV-2 in addition to IAV models (see Figure 16).

4. Conclusions

In this work, we revealed the mathematical structure of different SARS-CoV-2 infection dynamics models operating at different scales of in-host and host-to-host. We have shown how spatial PDE models can be coherently considered as well, resulting in an integrated hierarchical overview of models (Figure 15). This overview can be merged with those of other virus species (here, Influenza-A, Figure 16), which can serve as a beneficial instrument supporting further development of SARS-CoV-2 models.
Additional models similar to those studied here can be easily analyzed and added to the overview, because chemical organization theory can directly be applied to the model’s reaction network structure without requiring further kinetic details. However, complex multi-agent simulation models are problematic, because a corresponding reaction network cannot be derived in a straight forward way.
For each analyzed model we obtained a qualitative description of all possible long-term dynamics. That is, we distinguish sets of species that can persist (organizations) and sets of species that can definitely not persist, no matter what kind of quantitative kinetics are chosen.
We found a number of similarities among the models, for instance, the Abuin, Nesteruk, and Bai model are all in the same group regarding their organizations. Nevertheless, we also found a surprisingly high diversity of models with respect to their long-term qualitative behavior (Figure 15). Interestingly, there is only a small overlap between SARS-CoV-2 and Influenza A models. Furthermore, compared to the Influenza-A models the SARS-CoV-2 models appear to be simpler (mostly level 2) and thus display a simpler hierarchy, i.e., with only one inclusion relation between the Wu and the Su model (Figure 16).
Finally, there is only one model (the Wu Model) implying the unconditional persistence of infected cells. In stark contrast with this, three host-to-host (Nesteruk, Bai, Fitzgibbon-II) and one in-host model (Abuin) imply the definite extinction of any infection in the long-term.

Author Contributions

S.P. performed the implementation and computations. All authors analyzed the results and did the final conclusions. S.P. and B.I. wrote the paper with critical input from P.D. All authors reviewed and approved the final version of the paper. P.D. and B.I. supervised the project and B.I. conceived the study. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Zeiss Foundation, project “Eine virtuelle Werkstatt für die Digitalisierung in den Wissenschaften” (Durchbrüche 2018).

Acknowledgments

The authors are grateful to Hendrik Huthoff, Head of Education at Jena School for Microbial Communication, for critically reading the manuscript.

Conflicts of Interest

The authors declare no conflict of interest

Appendix A. List of the Reactions of All Models with Reactions Constants in Brackets

Figure A1. Lists of reactions of the Vargas Models.
Figure A1. Lists of reactions of the Vargas Models.
Viruses 13 00014 g0a1
Figure A2. Lists of reactions of the Abuin Model and the Su Model.
Figure A2. Lists of reactions of the Abuin Model and the Su Model.
Viruses 13 00014 g0a2
Figure A3. Lists of reactions of the Nesteruk Model, the Wu Model, and the Bai Model.
Figure A3. Lists of reactions of the Nesteruk Model, the Wu Model, and the Bai Model.
Viruses 13 00014 g0a3
Figure A4. Lists of reactions of the Bocharov Models.
Figure A4. Lists of reactions of the Bocharov Models.
Viruses 13 00014 g0a4
Figure A5. Lists of reactions of the Fitzgibbon Models.
Figure A5. Lists of reactions of the Fitzgibbon Models.
Viruses 13 00014 g0a5
Figure A6. Lists of reactions of the Almocera Model.
Figure A6. Lists of reactions of the Almocera Model.
Viruses 13 00014 g0a6

References

  1. Almocera, A.E.S.; Hernandez-Vargas, E.A. Multiscale model within-host and between-host for viral infectious diseases. J. Math. Biol. 2018, 77, 1035–1057. [Google Scholar] [CrossRef] [PubMed]
  2. Boianelli, A.; Nguyen, V.K.; Ebensen, T.; Schulze, K.; Wilk, E.; Sharma, N.; Stegemann-Koniszewski, S.; Bruder, D.; Toapanta, F.R.; Guzman, C.A.; et al. Modeling Influenza Virus Infection: A Roadmap for Influenza Research. Viruses 2015, 7, 5274–5304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Hernandez-Vargas, E.A.; Velasco-Hernandez, J.X. In-host Mathematical Modelling of COVID-19 in Humans. Annu. Rev. Control 2020, 50, 448–456. [Google Scholar] [CrossRef] [PubMed]
  4. Du, S.Q.; Yuan, W. Mathematical modeling of interaction between innate and adaptive immune responses in COVID-19 and implications for viral pathogenesis. J. Med. Virol. 2020, 92, 1615–1628. [Google Scholar] [CrossRef]
  5. Vargas, E.A.H.; Velasco-Hernandez, J.X. In-host modelling of covid-19 kinetics in humans. medRxiv 2020, 44487. [Google Scholar] [CrossRef] [Green Version]
  6. Tasevich, A.L.; Bocharov, G.A.; Vol’pert, V.A. Reaction-diffusion equations in immunology. Zhurnal Vychislitel’noi Matematiki i Matematicheskoi Fiziki 2018, 58, 2048–2059. [Google Scholar] [CrossRef]
  7. Abuin, P.; Anderson, A.; Ferramosca, A.; Hernandez-Vargas, E.A.; Gonzalez, A.H. Characterization of SARS-CoV-2 Dynamics in the Host. arXiv 2020, arXiv2006.08447. [Google Scholar]
  8. Su, Z.; Wu, Y. A multiscale and comparative model for receptor binding of 2019 novel coronavirus and the implication of its life cycle in host cells. BioRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  9. Nesteruk, I. Statistics-Based Predictions of Coronavirus Epidemic Spreading in Mainland China. Innov. Biosyst. Bioeng. 2020. [Google Scholar] [CrossRef] [Green Version]
  10. Wu, J.T.; Leung, K.; Leung, G.M. Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: A modelling study. Lancet 2020, 395, 689–697. [Google Scholar] [CrossRef] [Green Version]
  11. Fitzgibbon, W.; Morgan, J.; Webb, G.; Wu, Y. Analysis of a reaction–diffusion epidemic model with asymptomatic transmission. J. Biol. Syst. 2020, 28, 561–587. [Google Scholar] [CrossRef]
  12. Krishna, M.V.; Prakash, J. Mathematical modelling on phase based transmissibility of Coronavirus. Infect. Dis. Model. 2020, 5, 375–385. [Google Scholar] [CrossRef] [PubMed]
  13. Ewald, J.; Sieber, P.; Garde, R.; Lang, S.N.; Schuster, S.; Ibrahim, B. Trends in mathematical modeling of host–pathogen interactions. Cell. Mol. Life Sci. 2020, 77, 467–480. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Eker, S. Validity and usefulness of COVID-19 models. Humanit. Soc. Sci. Commun. 2020, 7, 1–5. [Google Scholar] [CrossRef]
  15. Ibrahim, B.; Henze, R.; Gruenert, G.; Egbert, M.; Huwald, J.; Dittrich, P. Spatial rule-based modeling: A method and its application to the human mitotic kinetochore. Cells 2013, 2, 506–544. [Google Scholar] [CrossRef]
  16. Wang, J. Mathematical models for COVID-19: Applications, limitations, and potentials. J. Public Health Emerg. 2020. [Google Scholar] [CrossRef]
  17. Hufsky, F.; Lamkiewicz, K.; Almeida, A.; Aouacheria, A.; Arighi, C.; Bateman, A.; Baumbach, J.; Beerenwinkel, N.; Brandt, C.; Cacciabue, M.; et al. Computational Strategies to Combat COVID-19: Useful Tools to Accelerate SARS-CoV-2 and Coronavirus Research. Briefings Bioinformat. 2020. [Google Scholar] [CrossRef]
  18. Ibrahim, B.; McMahon, D.P.; Hufsky, F.; Beer, M.; Deng, L.; Le Mercier, P.; Palmarini, M.; Thiel, V.; Marz, M. A new era of virus bioinformatics. Virus Res. 2018, 251, 86–90. [Google Scholar] [CrossRef] [Green Version]
  19. Peter, S.; Hölzer, M.; Lamkiewicz, K.; Di Fenizio, P.S.; Al Hwaeer, H.; Marz, M.; Schuster, S.; Dittrich, P.; Ibrahim, B. Structure and hierarchy of influenza virus models revealed by reaction network analysis. Viruses 2019, 11, 449. [Google Scholar] [CrossRef] [Green Version]
  20. Speroni di Fenizio, P.; Dittrich, P. Chemical Organizations at Different Spatial Scales. In Advances in Artificial Life; Springer: Berlin/Heidelberg, Germany, 2007; pp. 1–11. [Google Scholar]
  21. Di Fenizio, P.S.; Dittrich, P.; Banzhaf, W.; Ziegler, J. Towards a theory of organizations. In German Workshop on Artificial Life (GWAL 2000); DUV: Bayreuth, Germany, 2000. [Google Scholar]
  22. Peter, S.; Dittrich, P. On the Relation between Organizations and Limit Sets in Chemical Reaction Systems. Adv. Complex Syst. 2011, 14, 77–96. [Google Scholar] [CrossRef]
  23. Dittrich, P.; Speroni di Fenizio, P. Chemical Organization Theory. Bull. Math. Biol. 2007, 69, 1199–1231. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Kreyssig, P.; Wozar, C.; Peter, S.; Veloz, T.; Ibrahim, B.; Dittrich, P. Effects of small particle numbers on long-term behaviour in discrete biochemical systems. Bioinformatics 2014, 30, 475–481. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Peter, S.; Ghanim, F.; Dittrich, P.; Ibrahim, B. Organizations in reaction-diffusion systems: Effects of diffusion and boundary conditions. Ecol. Complex. 2020, 43, 100855. [Google Scholar] [CrossRef]
  26. Bai, S. Simulations of COVID-19 spread by spatial agent-based model and ordinary differential equations. Int. J. Simul. Process. Model. 2020, 15, 268–277. [Google Scholar] [CrossRef]
  27. Wölfel, R.; Corman, V.M.; Guggemos, W.; Seilmaier, M.; Zange, S.; Müller, M.A.; Niemeyer, D.; Jones, T.C.; Vollmar, P.; Rothe, C.; et al. Virological assessment of hospitalized patients with COVID-2019. Nature 2020, 581, 465–469. [Google Scholar] [CrossRef] [Green Version]
  28. Hernandez-Vargas, E.A. Modeling and Control of Infectious Diseases in the Host: With MATLAB and R; Academic Press: New York, NY, USA, 2019. [Google Scholar]
  29. Perelson, A.S. Modelling viral and immune system dynamics. Nat. Rev. Immunol. 2002, 2, 28–36. [Google Scholar] [CrossRef]
  30. Ciupe, S.M.; Heffernan, J.M. In-host modeling. Infect. Dis. Model. 2017, 2, 188–202. [Google Scholar] [CrossRef]
  31. Kreyssig, P.; Escuela, G.; Reynaert, B.; Veloz, T.; Ibrahim, B.; Dittrich, P. Cycles and the qualitative evolution of chemical systems. PLoS ONE 2012, 7, e45772. [Google Scholar] [CrossRef]
  32. Murray, J.D. Mathematical biology: I. An introduction. Interdisciplinary applied mathematics. In Mathematical Biology; Springer: New York, NY, USA, 2002. [Google Scholar]
  33. Bailey, N.T. The Mathematical Theory of Epidemics; Technical Report; Griffin: London, UK, 1957. [Google Scholar]
  34. Weitz, J.S.; Dushoff, J. Modeling post-death transmission of Ebola: Challenges for inference and opportunities for control. Sci. Rep. 2015, 5, 8751. [Google Scholar] [CrossRef]
  35. Shao, P.; Shan, Y. Beware of asymptomatic transmission: Study on 2019-nCoV prevention and control measures based on extended SEIR model. BioRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  36. Anderson, R.M.; Anderson, B.; May, R.M. Infectious Diseases of Humans: Dynamics and Control; Oxford University Press: New York, NY, USA, 1992. [Google Scholar]
  37. Handel, A.; Longini, I.M., Jr.; Antia, R. Towards a quantitative understanding of the within-host dynamics of influenza A infections. J. R. Soc. Interface 2009, 7, 35–47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Smith, A.M.; Smith, A.P. A Critical, Nonlinear Threshold Dictates Bacterial Invasion and Initial Kinetics During Influenza. Sci. Rep. 2016, 6, 38703. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Hancioglu, B.; Swigon, D.; Clermont, G. A dynamical model of human immune response to influenza A virus infection. J. Theor. Biol. 2007, 246, 70–86. [Google Scholar] [CrossRef] [PubMed]
  40. Lee, H.Y.; Topham, D.J.; Park, S.Y.; Hollenbaugh, J.; Treanor, J.; Mosmann, T.R.; Jin, X.; Ward, B.M.; Miao, H.; Holden-Wiltse, J.; et al. Simulation and prediction of the adaptive immune response to influenza A virus infection. J. Virol. 2009, 83, 7151–7165. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Cao, P.; Yan, A.W.; Heffernan, J.M.; Petrie, S.; Moss, R.G.; Carolan, L.A.; Guarnaccia, T.A.; Kelso, A.; Barr, I.G.; McVernon, J.; et al. Innate Immunity and the Inter-exposure Interval Determine the Dynamics of Secondary Influenza Virus Infection and Explain Observed Viral Hierarchies. PLoS Comput. Biol. 2015, 11, e1004334. [Google Scholar] [CrossRef] [PubMed]
  42. Miao, H.; Hollenbaugh, J.A.; Zand, M.S.; Holden-Wiltse, J.; Mosmann, T.R.; Perelson, A.S.; Wu, H.; Topham, D.J. Quantifying the early immune response and adaptive immune response kinetics in mice infected with influenza A virus. J. Virol. 2010, 84, 6687–6698. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Bocharov, G.A.; Romanyukha, A.A. Mathematical model of antiviral immune response. III. Influenza A virus infection. J. Theor. Biol. 1994, 167, 323–360. [Google Scholar] [CrossRef]
  44. Hernandez-Vargas, A.E.; Meyer-Hermann, M. Innate immune system dynamics to influenza virus. Ifac Proc. Vol. 2012, 45, 260–265. [Google Scholar] [CrossRef] [Green Version]
  45. Baccam, P.; Beauchemin, C.; Macken, C.A.; Hayden, F.G.; Perelson, A.S. Kinetics of influenza A virus infection in humans. J. Virol. 2006, 80, 7590–7599. [Google Scholar] [CrossRef] [Green Version]
  46. Pawelek, K.A.; Huynh, G.T.; Quinlivan, M.; Cullinane, A.; Rong, L.; Perelson, A.S. Modeling within-host dynamics of influenza virus infection including immune responses. PLoS Comput. Biol. 2012, 8, e1002588. [Google Scholar] [CrossRef] [Green Version]
  47. Saenz, R.A.; Quinlivan, M.; Elton, D.; Macrae, S.; Blunden, A.S.; Mumford, J.A.; Daly, J.M.; Digard, P.; Cullinane, A.; Grenfell, B.T.; et al. Dynamics of influenza virus infection and pathology. J. Virol. 2010, 84, 3974–3983. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. The in-host part of the Almocera Model [1,2] has two variables resp. ODEs (see Subfigure (a)): viruses ( V ) and T-cells ( E ). This is the starting point of our method consisting of three steps briefly described below. Step 1: We derive from the ODE system a set of six reactions (see Subfigure (b)): R 1 , , R 6 . Step 2: We compute from the set of reactions the set of organizations (an organization is a subset of species with specific properties as described below in this Chapter): { V , E } and { E } . We arrange the organizations in a Hasse diagram (see Subfigure (c)), where organizations get bigger from bottom to top and are linked by a line, if the lower organization is a subset of the upper one. Step 3: We derive from the set of organizations the signature of the model (see Subfigure (d)). For our example, the signature is X , X X , where X represents the organization { E } and X X represents { V , E } . The signature tells us via colored Xs, which of the types ( uninfected or susceptible cells or individuals , viruses or infected cells or individuals, or immune system, e.g., T-cells) of species are contained in the organizations of the model. We maintain this coloring throughout this work. Note that we use underlining X ̲ to tag host-to-host species in contrast to in-host species. One should understand the following aspects concerning the method presented above: The long-term behavior of simulations of the dynamics of the model can be easily estimated from the signature: We know from [22] that there is always an organization representing the species persisting in the long-run. Thus species that are not contained in any organization will go extinct for sure after a sufficiently long time period. On the other hand, species that are contained in all organizations of a model, will persist in the long-run for sure. If a species is contained in some organizations of a model but not in all, it has the potential to persist but also to go extinct. It depends on the applied kinetic laws, the initial conditions, and the reaction constants, which case occurs. A hierarchy of a set of models can be constructed relying on their signatures. Like the set of organizations of one model, it can again be visualized as a Hasse diagram. The more combinations of colors a signature contains, the higher is its position in the Hasse diagram and the bigger is the variety of its potential dynamical behavior. If all combinations of colors of one signature are present in a second signature, then the models can be linked by a line. Kinetic laws: Note that except for reaction R 6 , where Michaelis–Menten kinetics are applied, all the other reactions of this example follow mass-action kinetics. The technique of computing and analyzing chemical organizations used in this work applies to both these kinetic laws. Distributed organizations: When the species V and E are separated (we say “distributed”), the reaction R 3 is inactive. Then, due to the remaining for reactions, the set { V , E } is still self-maintaining and closed and thus some kind of an organization. We write { V | E } (instead of { V , E } ) to denote this and call { V | E } a “distributed organization.”
Figure 1. The in-host part of the Almocera Model [1,2] has two variables resp. ODEs (see Subfigure (a)): viruses ( V ) and T-cells ( E ). This is the starting point of our method consisting of three steps briefly described below. Step 1: We derive from the ODE system a set of six reactions (see Subfigure (b)): R 1 , , R 6 . Step 2: We compute from the set of reactions the set of organizations (an organization is a subset of species with specific properties as described below in this Chapter): { V , E } and { E } . We arrange the organizations in a Hasse diagram (see Subfigure (c)), where organizations get bigger from bottom to top and are linked by a line, if the lower organization is a subset of the upper one. Step 3: We derive from the set of organizations the signature of the model (see Subfigure (d)). For our example, the signature is X , X X , where X represents the organization { E } and X X represents { V , E } . The signature tells us via colored Xs, which of the types ( uninfected or susceptible cells or individuals , viruses or infected cells or individuals, or immune system, e.g., T-cells) of species are contained in the organizations of the model. We maintain this coloring throughout this work. Note that we use underlining X ̲ to tag host-to-host species in contrast to in-host species. One should understand the following aspects concerning the method presented above: The long-term behavior of simulations of the dynamics of the model can be easily estimated from the signature: We know from [22] that there is always an organization representing the species persisting in the long-run. Thus species that are not contained in any organization will go extinct for sure after a sufficiently long time period. On the other hand, species that are contained in all organizations of a model, will persist in the long-run for sure. If a species is contained in some organizations of a model but not in all, it has the potential to persist but also to go extinct. It depends on the applied kinetic laws, the initial conditions, and the reaction constants, which case occurs. A hierarchy of a set of models can be constructed relying on their signatures. Like the set of organizations of one model, it can again be visualized as a Hasse diagram. The more combinations of colors a signature contains, the higher is its position in the Hasse diagram and the bigger is the variety of its potential dynamical behavior. If all combinations of colors of one signature are present in a second signature, then the models can be linked by a line. Kinetic laws: Note that except for reaction R 6 , where Michaelis–Menten kinetics are applied, all the other reactions of this example follow mass-action kinetics. The technique of computing and analyzing chemical organizations used in this work applies to both these kinetic laws. Distributed organizations: When the species V and E are separated (we say “distributed”), the reaction R 3 is inactive. Then, due to the remaining for reactions, the set { V , E } is still self-maintaining and closed and thus some kind of an organization. We write { V | E } (instead of { V , E } ) to denote this and call { V | E } a “distributed organization.”
Viruses 13 00014 g001
Figure 2. Overview of the relations between organizations and virus dynamics research. The in vivo and in vitro measurements and observations (left) lead to quantitative data as well as models (center) concerning the virus dynamics. Those models in turn exist at different scales of abstraction, such as non spatial (like ODEs) or spatial modeling (PDEs) but also in-host (upper left) and host-to-host modeling. If a system of partial differential equations can be solved numerically (right), which can be challenging, COT provides a tool to further study the dynamics (upper right). When solving the system of partial differential equations becomes too difficult, COT allows for describing the principal components of the system and provides conclusions about its dynamical behavior in simple ways (right). Additionally, COT works without any kind of kinetic data, such as reaction or diffusion rates.
Figure 2. Overview of the relations between organizations and virus dynamics research. The in vivo and in vitro measurements and observations (left) lead to quantitative data as well as models (center) concerning the virus dynamics. Those models in turn exist at different scales of abstraction, such as non spatial (like ODEs) or spatial modeling (PDEs) but also in-host (upper left) and host-to-host modeling. If a system of partial differential equations can be solved numerically (right), which can be challenging, COT provides a tool to further study the dynamics (upper right). When solving the system of partial differential equations becomes too difficult, COT allows for describing the principal components of the system and provides conclusions about its dynamical behavior in simple ways (right). Additionally, COT works without any kind of kinetic data, such as reaction or diffusion rates.
Viruses 13 00014 g002
Figure 3. The Vargas-I Model [5,27] with two variables: Exponential growth viruses (Vg) and decay of viruses (Vd). There are two organizations: the empty set and the single species set {Vg}. The set of all species is not an organization since Vd decays but is not produced by any reaction. The signature of this model is: , X .
Figure 3. The Vargas-I Model [5,27] with two variables: Exponential growth viruses (Vg) and decay of viruses (Vd). There are two organizations: the empty set and the single species set {Vg}. The set of all species is not an organization since Vd decays but is not produced by any reaction. The signature of this model is: , X .
Viruses 13 00014 g003
Figure 4. The Vargas-II Model [1,2,5] with two variables: viruses V and T-cells T, and two organizations: { V , T } and { T } . The empty set is not an organization for this model since T has an inflow reaction with reaction constant s T and thus does not go extinct. The organization { V , T } exists also as distributed organization { V | T } . So, if V and T are separated the two reactions with reaction constants c T and r are inactive, but this does neither destroy the self-maintenance nor the closedness. The signature of this model is: X , X X . Note that, by replacing T by E in this model, we get almost the same reactions and the same stoichiometric matrix as for the in-host Almocera Model we introduced in the Materials and Methods Section 2.
Figure 4. The Vargas-II Model [1,2,5] with two variables: viruses V and T-cells T, and two organizations: { V , T } and { T } . The empty set is not an organization for this model since T has an inflow reaction with reaction constant s T and thus does not go extinct. The organization { V , T } exists also as distributed organization { V | T } . So, if V and T are separated the two reactions with reaction constants c T and r are inactive, but this does neither destroy the self-maintenance nor the closedness. The signature of this model is: X , X X . Note that, by replacing T by E in this model, we get almost the same reactions and the same stoichiometric matrix as for the in-host Almocera Model we introduced in the Materials and Methods Section 2.
Viruses 13 00014 g004
Figure 5. The Abuin Model [5,7,28,29,30] with three variables: susceptible host cells U, infected host cells I, and viral particles V. There are two organizations: the empty set { } and { U } . None of them includes an active reaction; thus we say, that they are “non-reactive”. Note that for this model, a principal part of the infection dynamics, concerning I and V, does not take place within an organization. Thus, from the role organizations play in the long-run of dynamical systems [22,31] we know that this model induces a vanishing of I and V in the long run. The signature of this model is: , X .
Figure 5. The Abuin Model [5,7,28,29,30] with three variables: susceptible host cells U, infected host cells I, and viral particles V. There are two organizations: the empty set { } and { U } . None of them includes an active reaction; thus we say, that they are “non-reactive”. Note that for this model, a principal part of the infection dynamics, concerning I and V, does not take place within an organization. Thus, from the role organizations play in the long-run of dynamical systems [22,31] we know that this model induces a vanishing of I and V in the long run. The signature of this model is: , X .
Viruses 13 00014 g005
Figure 6. The Su Model [8] with 7 variables: healthy cells H, bound virus Vb , RNA genome m, proteins and replicated RNA packaged together in cytoplasm P, cytokines stimulating inflammatory responses S, infected cells In , coronavirus V. Here we have three organizations and one distributed organization { H , m , S | In } that is not an organization. So, the species H, m, S, and In can only survive together, if S and In are separated. This deactivates the reaction with the reaction constant d I and thus In is able to persist. The signature of this model is: X , X X .
Figure 6. The Su Model [8] with 7 variables: healthy cells H, bound virus Vb , RNA genome m, proteins and replicated RNA packaged together in cytoplasm P, cytokines stimulating inflammatory responses S, infected cells In , coronavirus V. Here we have three organizations and one distributed organization { H , m , S | In } that is not an organization. So, the species H, m, S, and In can only survive together, if S and In are separated. This deactivates the reaction with the reaction constant d I and thus In is able to persist. The signature of this model is: X , X X .
Viruses 13 00014 g006
Figure 7. The Bocharov-I Model [6] with one variable: the virus concentration v. There are no boundary conditions specified. Thus, we assume Neumann boundary conditions for simplicity and in the style of the other PDE models analyzed in this work. There is the maximum number of two organizations for a model with one species here. This represents all possible long-term dynamics of the infection, i.e., its persistence as well as its extinction. The signature of this model is: , X .
Figure 7. The Bocharov-I Model [6] with one variable: the virus concentration v. There are no boundary conditions specified. Thus, we assume Neumann boundary conditions for simplicity and in the style of the other PDE models analyzed in this work. There is the maximum number of two organizations for a model with one species here. This represents all possible long-term dynamics of the infection, i.e., its persistence as well as its extinction. The signature of this model is: , X .
Viruses 13 00014 g007
Figure 8. The Bocharov-II Model [5] from the year 2018 with two variables: the virus concentration v and immune cell concentration c. The functions ϕ ( v ) and ψ ( v ) are assumed to be strictly positive if and only if v > 0 . As for the Bocharov-I Model (see Figure 7) we assume Neumann boundary conditions. There are three organizations. Only one subset of species is not an organization, i.e., the set { c } . Thus the model provides a relatively big variety of possible long-term behaviors. The signature of this model is: , X , X X .
Figure 8. The Bocharov-II Model [5] from the year 2018 with two variables: the virus concentration v and immune cell concentration c. The functions ϕ ( v ) and ψ ( v ) are assumed to be strictly positive if and only if v > 0 . As for the Bocharov-I Model (see Figure 7) we assume Neumann boundary conditions. There are three organizations. Only one subset of species is not an organization, i.e., the set { c } . Thus the model provides a relatively big variety of possible long-term behaviors. The signature of this model is: , X , X X .
Viruses 13 00014 g008
Figure 9. The Nesteruk Model [9,32,33] with 3 variables (SIR): the number of susceptible persons S, infected (sick and infection-spreading) persons I, and removed (sum of isolated, recovered, and dead) persons R. There are four organizations that only contain healthy individuals. None of them contains infected individuals and all are non-reactive, i.e., no reaction is active for these organizations. Thus, as for the in-host Abuin Model from the previous section, the whole infection dynamics takes place outside the organizations. This model inherently assumes that all infected individuals I will vanish finally. The signature of this model is: , X ̲ .
Figure 9. The Nesteruk Model [9,32,33] with 3 variables (SIR): the number of susceptible persons S, infected (sick and infection-spreading) persons I, and removed (sum of isolated, recovered, and dead) persons R. There are four organizations that only contain healthy individuals. None of them contains infected individuals and all are non-reactive, i.e., no reaction is active for these organizations. Thus, as for the in-host Abuin Model from the previous section, the whole infection dynamics takes place outside the organizations. This model inherently assumes that all infected individuals I will vanish finally. The signature of this model is: , X ̲ .
Viruses 13 00014 g009
Figure 10. The Wu Model [10] (a SEIR model) with 4 variables: the number of susceptible S, latent E, infectious I, and removed R individuals, and only one organization: { S , E , I . R } , that contains all species. Thus this model implicitly assumes the infection to persist forever once it occurs which is a totally contrary assumption compared to the models assuming the vanishing of infection in the like the previously analyzed Nesteruk model. The signature of this model is: X ̲ X ̲ .
Figure 10. The Wu Model [10] (a SEIR model) with 4 variables: the number of susceptible S, latent E, infectious I, and removed R individuals, and only one organization: { S , E , I . R } , that contains all species. Thus this model implicitly assumes the infection to persist forever once it occurs which is a totally contrary assumption compared to the models assuming the vanishing of infection in the like the previously analyzed Nesteruk model. The signature of this model is: X ̲ X ̲ .
Viruses 13 00014 g010
Figure 11. The Bai Model [26,34,35] with 5 variables (SEIRD): the number of susceptible S, exposed E, infected I, recovered R, and dead D individuals. This model has a similar structure in terms of organizations to the Nesteruk model: there is no organization containing species representing the infection. Thus infection is implicitly assumed to vanish in the long-term. The remaining multitude of organizations exists simply due to the fact that recovered individuals R and dead individuals D can be combined arbitrarily with each other and with susceptible individuals S to form organizations. The signature of this model is: , X ̲ .
Figure 11. The Bai Model [26,34,35] with 5 variables (SEIRD): the number of susceptible S, exposed E, infected I, recovered R, and dead D individuals. This model has a similar structure in terms of organizations to the Nesteruk model: there is no organization containing species representing the infection. Thus infection is implicitly assumed to vanish in the long-term. The remaining multitude of organizations exists simply due to the fact that recovered individuals R and dead individuals D can be combined arbitrarily with each other and with susceptible individuals S to form organizations. The signature of this model is: , X ̲ .
Viruses 13 00014 g011
Figure 12. The Fitzgibbon-I Model [11] with one variable ρ representing the current strength of the infection. ρ obeys Neumann boundary conditions. The reaction network structure of this model is almost equal to that of the in-host PDE Bocharov-I model (see Figure 7) from the previous section. Thus there is the maximum number of two organizations. All long-term dynamics of the infection, i.e., its persistence as well as its extinction, are possible. The signature of this model is: , X ̲ .
Figure 12. The Fitzgibbon-I Model [11] with one variable ρ representing the current strength of the infection. ρ obeys Neumann boundary conditions. The reaction network structure of this model is almost equal to that of the in-host PDE Bocharov-I model (see Figure 7) from the previous section. Thus there is the maximum number of two organizations. All long-term dynamics of the infection, i.e., its persistence as well as its extinction, are possible. The signature of this model is: , X ̲ .
Viruses 13 00014 g012
Figure 13. The Fitzgibbon-II Model [11] with two variables: susceptible individuals S ( t ) and infected individuals I ( x , t ) . Note that only the infected individuals I are modeled dependent not only of time but also of space. It follows Neumann boundary condition. As for the Nesteruk Model and the Bai Model we here have no organization with any species representing the infection. Thus this model implies the infection to go extinct. The signature of this model is: , X ̲ .
Figure 13. The Fitzgibbon-II Model [11] with two variables: susceptible individuals S ( t ) and infected individuals I ( x , t ) . Note that only the infected individuals I are modeled dependent not only of time but also of space. It follows Neumann boundary condition. As for the Nesteruk Model and the Bai Model we here have no organization with any species representing the infection. Thus this model implies the infection to go extinct. The signature of this model is: , X ̲ .
Viruses 13 00014 g013
Figure 14. The linked Almocera Model [1,36] has the two variables of the in-host model above (see Figure 1), i.e., viruses (V) and T-cells (E), plus two further variables: susceptible (S) and infected (I) individuals. Note that viruses V as a part of the in-host model influence the host-to-host dynamics via the reaction with reaction constant r W . There are three organizations: All of them contain E and S, since they are produced by the two inflow reactions with the reaction constants N E and N S , respectively. If V is added to { E , S } , we get the organization { V , E , S } that can exist as a distributed organization by separation of V and E in the same way as it was the case for the in-host model (see Figure 1). If, furthermore, I is added, then we get the full organization where still V and E can be separated, but I must not be separated from V and S since then it could not regenerate vie the reaction with the reaction constant r W . Interestingly, in all the organizations of this model in-host species are mixed with host-to-host species. We find that in the long-term the species representing healthy, uninfected hosts, i.e., E and S always persist. Contrarily, V and I might persist too, but might also go extinct. If this is the case for the virus V than also the infected individuals I go extinct. This must not be the case the other way around. Lastly, the COT analysis shows that the model assumes that the T-cells E can exist independently of the virus, but infected individuals can only exist persistently if in contact with healthy individuals S and viruses V. The signature of this model is: , X ̲ . The signature of this model is: X X ̲ , X X X ̲ , X X X X ̲ .
Figure 14. The linked Almocera Model [1,36] has the two variables of the in-host model above (see Figure 1), i.e., viruses (V) and T-cells (E), plus two further variables: susceptible (S) and infected (I) individuals. Note that viruses V as a part of the in-host model influence the host-to-host dynamics via the reaction with reaction constant r W . There are three organizations: All of them contain E and S, since they are produced by the two inflow reactions with the reaction constants N E and N S , respectively. If V is added to { E , S } , we get the organization { V , E , S } that can exist as a distributed organization by separation of V and E in the same way as it was the case for the in-host model (see Figure 1). If, furthermore, I is added, then we get the full organization where still V and E can be separated, but I must not be separated from V and S since then it could not regenerate vie the reaction with the reaction constant r W . Interestingly, in all the organizations of this model in-host species are mixed with host-to-host species. We find that in the long-term the species representing healthy, uninfected hosts, i.e., E and S always persist. Contrarily, V and I might persist too, but might also go extinct. If this is the case for the virus V than also the infected individuals I go extinct. This must not be the case the other way around. Lastly, the COT analysis shows that the model assumes that the T-cells E can exist independently of the virus, but infected individuals can only exist persistently if in contact with healthy individuals S and viruses V. The signature of this model is: , X ̲ . The signature of this model is: X X ̲ , X X X ̲ , X X X X ̲ .
Viruses 13 00014 g014
Figure 15. Hierarchy of SARS-CoV-2 models (tagged with CoV()) and general virus models with respect to their long-term behavior identified by their signature. There are three different levels of increasing signatures. The higher the position of a model in the hierarchy the more diverse is its potential dynamical behavior. By clicking on the model names you are directed to the respective part of this work where the model is analyzed. PDE model names are written in italic.
Figure 15. Hierarchy of SARS-CoV-2 models (tagged with CoV()) and general virus models with respect to their long-term behavior identified by their signature. There are three different levels of increasing signatures. The higher the position of a model in the hierarchy the more diverse is its potential dynamical behavior. By clicking on the model names you are directed to the respective part of this work where the model is analyzed. PDE model names are written in italic.
Viruses 13 00014 g015
Figure 16. Merged Hasse-diagram of hierarchy of SARS-CoV-2 (tagged with CoV()), IAV models (tagged with I A V ( ) ), and general virus models. They are positioned at five different levels according to the size of their signatures. The higher the level of a model, the bigger is the number of components of its signature and thus the diversity of its potential dynamical behavior. By clicking on the model names you are directed to the respective part of this work where the models are analyzed [37,38,39,40,41,42,43,44,45,46,47]. PDE model names are written in italic. Note, that contrary to many Influenza-A infection models (third level and above) all SARS-CoV-2 infection models are on the second level (counted from bottom to top) and thus have a maximum number of two different species in their signature. This means that they are less complex in terms of their organizations than most of the Influenza-A infection models. Note that the Smith Influenza-A model is the only one that considers co-infection by bacteria X. For more information about the COT analysis of the Influenza-A infection models see [19].
Figure 16. Merged Hasse-diagram of hierarchy of SARS-CoV-2 (tagged with CoV()), IAV models (tagged with I A V ( ) ), and general virus models. They are positioned at five different levels according to the size of their signatures. The higher the level of a model, the bigger is the number of components of its signature and thus the diversity of its potential dynamical behavior. By clicking on the model names you are directed to the respective part of this work where the models are analyzed [37,38,39,40,41,42,43,44,45,46,47]. PDE model names are written in italic. Note, that contrary to many Influenza-A infection models (third level and above) all SARS-CoV-2 infection models are on the second level (counted from bottom to top) and thus have a maximum number of two different species in their signature. This means that they are less complex in terms of their organizations than most of the Influenza-A infection models. Note that the Smith Influenza-A model is the only one that considers co-infection by bacteria X. For more information about the COT analysis of the Influenza-A infection models see [19].
Viruses 13 00014 g016
Table 1. Overview of all models analyzed in this work each named by its first author and followed by the names of the variables used in the models. By clicking on the model names or the model types (left) you are directed to the part of this work where the respective model is analyzed. The model names tagged with CoV() in the beginning are explicitly for SARS-CoV-2 infection whereas the others are designed for viral infections in general. All models except for the two models from Bocharov and Almocera (both published in 2018) were published in 2020.
Table 1. Overview of all models analyzed in this work each named by its first author and followed by the names of the variables used in the models. By clicking on the model names or the model types (left) you are directed to the part of this work where the respective model is analyzed. The model names tagged with CoV() in the beginning are explicitly for SARS-CoV-2 infection whereas the others are designed for viral infections in general. All models except for the two models from Bocharov and Almocera (both published in 2018) were published in 2020.
Model TypeODEPDE
In-host[5] CoVVargas-I: V g , V d [6] Bocharov-I: v
[5] CoVVargas-II: V , T [6] Bocharov-II: v , c
[7]CoVAbuin: U , I , V
[8] CoVSu: H , Vb , m , P , S , In , V
Host-to-host[9] CoVNesteruk: S , I , R [11] CoVFitzgibbon-I: ρ
[10] CoVWu: S , E , I , R [11] CoVFitzgibbon-II: S , I
[26] CoVBai: S , E , I , R , D
Linked[1] Almocera:
In-Host: V , E ; Linked: V , E , S , I
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Peter, S.; Dittrich, P.; Ibrahim, B. Structure and Hierarchy of SARS-CoV-2 Infection Dynamics Models Revealed by Reaction Network Analysis. Viruses 2021, 13, 14. https://0-doi-org.brum.beds.ac.uk/10.3390/v13010014

AMA Style

Peter S, Dittrich P, Ibrahim B. Structure and Hierarchy of SARS-CoV-2 Infection Dynamics Models Revealed by Reaction Network Analysis. Viruses. 2021; 13(1):14. https://0-doi-org.brum.beds.ac.uk/10.3390/v13010014

Chicago/Turabian Style

Peter, Stephan, Peter Dittrich, and Bashar Ibrahim. 2021. "Structure and Hierarchy of SARS-CoV-2 Infection Dynamics Models Revealed by Reaction Network Analysis" Viruses 13, no. 1: 14. https://0-doi-org.brum.beds.ac.uk/10.3390/v13010014

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop