SARS-CoV-2 Nsp5 Protein Causes Acute Lung Inflammation: a Dynamical Mathematical Model

In the present work we propose a dynamical mathematical model of the lung cells inflammation process in response to SARS-CoV-2 infection. In this scenario the main protease Nsp5 enhances the inflammatory process, increasing the levels of NF  B, IL-6, Cox2, and PGE2 with respect to a reference state without the virus. In presence of the virus the translation rates of NF  B and I  B arise to a high constant value, and when the translation rate of IL-6 also increases above the threshold value of 7 pg mL -1 s -1 the model predicts a persistent over stimulated immune state with high levels of the cytokine IL-6. Our model shows how such over stimulated immune state becomes autonomous of the signals from other immune cells such as macrophages and lymphocytes, and does not shut down by itself. We also show that in the context of the dynamical model presented here, Dexamethasone or Nimesulide have little effect on such inflammation state of the infected lung cell, and the only form to suppress it is with the inhibition of the activity of the viral protein Nsp5.To that end, our model suggest that drugs like Saquinavir may be useful. In this form, our model suggests that Nsp5 is effectively a central node underlying the severe acute lung inflammation during SARS-CoV-2 infection. The persistent production of IL-6 by lung cells can be one of the causes of the cytokine storm observed in critical patients with COVID19. Nsp5 seems to be the switch to start inflammation, the consequent overproduction of the ACE2 receptor, and an important underlying cause of the most severe


Introduction
Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) virus is an intracellular infectious agent whose replication cycle depends on the host´s cell structures and functions.In particular, it uses the translational apparatus of different types of infected cells to express its proteins (Nakagawa et al., 2016).SARS-CoV-2 causes the Coronavirus Disease 2019 (COVID19) that has infected over 117,000,000 persons and killed over 2,400,000 worldwide since the end of 2019.No specific and effective therapeutic drugs to defeat SARS-CoV-2 infection have been developed yet (Díaz, 2020a), although several effective vaccines have been developed and are being applied.
When SARS-CoV-2 virion infects the organism, S protein binds with high affinity to the surface receptor Angiotensin-Converting Enzyme 2 (ACE2), highly abundant in the lung alveolar type II cells and other cell types (Gordon et al., 2020;Hamming et al, 2004), and forms a molecular complex that begins the process of fusion of the virion envelope with the host cell membrane.
Finally, viral (+)ssRNA is released into the host cytoplasm (Letko et al., 2020) (Figure 1).Infection begins when the protein S of the virion binds with high affinity to the cells surface receptor ACE2 (Angiotensin-Converting Enzyme 2) that is highly abundant in the lung alveolar type II cells.The formation of the complex S-ACE2 initiates the process of fusion between the virion envelope and the cells membrane leading to the liberation of the nucleocapside with the viral genome into the cytoplasm .SARS-CoV-2 genome consists of a positive-sense nonsegmented single stranded mRNA ((+)ssRNA) of ~ 30 kb.The open reading frames 1a (orf1a) and 1b (orf1b) are located near the 5'UTR of the (+)ssRNA and they code for the polyproteins pp1a and pp1ab.Maturation of these polyproteins results in 11 nonstructural proteins (Nsp) from the orf1a segment (Nsp1 to Nsp11) and 5 nonstructural proteins from the orf1b segment (Nsp12 to Nsp16).Nsp proteins form the replicationtranscription complex (RTC) a set of nested subgenomic minus-strands of RNA ((-)sgRNA) are synthesized in a process of discontinuous transcription.These (-) sgRNAs serve as the templates for the production of subgenomic mRNAs from which the structural proteins E, M, N and S, together with the accessory proteins orf 3a, orf6, orf7a, orf7b, orf8, orf9b, orf9c and orf10 are synthesized.The presence of viral (+)ssRNA initiates the TLR3 and TLR4 mediated immune response, which increases the concentration of free NF B in cytoplasm.Cytoplasmic NF B enters the nucleus where promotes the transcription of a series of target genes.Viral protein Nsp5 inhibits the inhibition exerted by HDAC2 on nuclear NF B transcription factor, enhancing the synthesis of inflammatory cytokines like IL-6.(Figure drawn with Biorender template https://app.biorender.com/).
The process of translation occurs in the cytoplasm and produces a set of viral polyproteins whose maturation results in 11 nonstructural proteins (Nsp) from the orf1a segment (Nsp1 to Nsp11) and 5 nonstructural proteins from the orf1b segment (Nsp12 to Nsp16).Nsp proteins form the replicationtranscription complex (RTC) in a double-membrane vesicle where a set of nested subgenomic minus-strands of RNA ((-)sgRNA) are synthesized in a process of discontinuous transcription.These (-) sgRNAs serve as the templates for the production of subgenomic mRNAs from which the structural proteins E, M, N and S, together with the accessory proteins orf 3a, orf6, orf7a, orf7b, orf8, orf9b, orf9c and orf10 are synthesized (Díaz, 2020a;Sevajol et al., 2014;Dongwan et al., 2020).SARS-CoV-2 uses the host translational machinery to redirect it to viral protein synthesis and replication, while host mRNA translation is inhibited (Nakagawa et al., 2016).SARS-CoV and SARS-CoV-2 genomes have ~ 79% of homology (Forster et al., 2020), and most of the set of structural and nonstructural proteins are practically the same.However, the virus species differ in the accessory proteins orf8, orf8a, orf 8b, orf9c and orf10 (Bartlam et al., 2005;Gordon et al., 2020).
Viral proteins are inserted into the host molecular machinery to modify and redirect a great number of host cell functions towards the production of additional viral particles (Masters, 2006;Wu et al., 2020).Experimental analysis of the interaction of viral and host proteins, or interactome, by Gordon and collaborators (Gordon et al., 2020) has been a fundamental contribution to understand the form in which SARS-CoV-2 virus takes control of the host molecular network to produce new virions and propagate the infection (Hekman et al., 2020).Gordon and collaborators (2020) cloned, tagged and expressed 26 viral proteins in human cells using affinity-purification mass spectrometry to identify the human proteins physically associated with each other.They found around 332 SARS-CoV-2-human protein-protein interactions that constitute the virus interactome.From these results, the construction of the network representation of the interactome is possible.In the particular case of the Gordon interactome, the analysis of its undirected network model indicates a modular freescale hierarchical type of structure in which proteins orf8, N and Nsp7 are the main hubs (Díaz, 2020a).
Despite all the molecular knowledge rapidly accumulated on SARS-CoV-2 itself and the infection or inflammation processes in which it is involved in human cells, we still lack dynamical models to understand which are the key nodes underlying severe inflammation processes in some cases.Mathematical and computational dynamic network models are useful tools to integrate data, find experimental holes and also propose or understand the underlying dynamics mechanisms involved in the complex viral/human networks involved during SARS-CoV-2 infection (Breitling, 2010;Weinstein et al., 2020;Enciso et al., 2020;Díaz, 2020a).The number of nodes, the number of connections of each node to its neighbours, and the distribution of these connections in the network determine its complexity, structure and dynamical properties (Kumar et al., 2015).
In this network, the viral main protease Nsp5 has a special role in lung inflammation process through its link with Histone Deacetylase 2 (HDAC2) (Díaz, 2020a;Gordon et al., 2020;Hekman et al., 2020)   The change from a hub (N protein in SARS-CoV) to a poor connected protein (Nsp5 in SARS-CoV-2) as the underlying cause of severe acute respiratory syndrome could be an adaptation that made SARS-CoV-2 more pathogenic.Nsp5 is also the main protease that cleavages the SARS-CoV-2 polyprotein at 11 conserved sites, and its malfunction or deletion could stop the viral replication cycle and lung inflammation (Cohen et al., 2000;Xu et al., 2020).Thus, Nsp5 is a molecular switch that shows a low rate of mutation (Vilar e Isom, 2020), and its malfunctioning can decrease the strength of inflammation and the persistence of ACE2 in the lung epithelial cells giving rise to a less intense inflammatory response (Long et al., 2020).
Activation of the NF B family of proteins (p65/RelA, p105/p50, p100/p52, RelB, and c-Rel) is mediated by the phosphorylation of the inhibitor of NF B proteins (IBs) by IB kinases (IKKs) and the posterior ubiquitinylation and degradation of the phosphorylated IBs (Figure 1).The released NF κB proteins enter the nucleus to activate specific genes (Wagner et al., 2015).In the nucleus, HDAC2 blocks the transcriptional activity of NF B, inhibiting the production and release of the cytokine Interleukin 6 (IL-6) and IB, and interfering with the functioning of the positive feedback circuit between IL-6 and NF B (Rahman, and MacNee, 1998) (Figure 3a).As we mentioned before, the sustained inactivation of HDAC2 by Nsp5 produces a sustained increase in the transcriptional activity of NF κB (Figures1 and 2) that leads to a high production of IL-6 in the lung epithelium cells.Release of this excess of IL-6 to the vascular system contributes to the "cytokines storm" observed in critical patients (Magro, 2020;Leisman et al., 2020), although the average serum concentration of IL-6 in critical COVID19 patients are lower with respect to other respiratory syndromes (Leisman et al., 2020) .However, the uninterrupted production of IL-6 is closely associated to a persistent presence of the ACE2 receptor in the alveolar type 2 cells due to enhanced ACE2 gene expression mediated by the JAK-STAT pathway (Hennighausen and Lee, 2020) (Figure 2).
NF B also promotes transcription of the Cox2 gene leading to the expression and increased enzymatic activity of Cox2 (Figures 2 and 3B).However, NF B, Cox2 and IL-6 also form a positive feedback circuit in which NF B promotes Cox2 enzymatic activity that catalyses the conversion of Arachidonic acid (AA) into Prostaglandin G2 (PGG2), which is in turn modified by the peroxidase moiety of the Cox2 enzyme to produce prostaglandin H2 (PGH2) that is converted to prostaglandin E2 (PGE2) (Alexanian and Sorokin, 2017).PGE2 is then released to the vascular system, binds to its membrane EP2 and EP4 prostanoid receptors and promotes the cleavage of the IB-NF B complex increasing the amount of free NF B, which increases the production of IL-6 (Figure 3B) (Cho et al., 2014;Bouffi et al., 2010).IL-6, in turn, increases NF B activity (Wang et al., 2003) (Figure 2).These regulatory circuits has three therapeutic targets: IL-6 production (Tocilizumab) (Magro 2020), Cox2 (Nimesulide) (Alexanian and Sorokin, 2017) and NF B (Dexamethasone) (Newton et al., 1998;Aghai et al, 2006).However, none of these have shown to be really effective to prevent or revert acute inflammatory cases after SARS-CoV-2 infection.In this work, we propose a dynamical mathematical model that suggests that the three feedback loops mentioned above constitute the minimal circuit for the inflammatory process in the epithelial lung cells during SARS-CoV-2 infection (Figure 2), and we use an ordinary differential equations (ODEs) continuous model to explore the effect of Nsp5 in the qualitative dynamics of this circuit in which this viral protein acts as an enhancing perturbation in the phase space of this dynamical system (Diaz, 2020b) (Figure 2).For the model of the circuit, we chose as main molecular components the proteins NF B, IB, IL-6, Cox2, and PGE2.In this model, Nsp5 prevents the movement of HDAC2 into the nucleus allowing an increase in NF B transcriptional activity.The input that turns on the circuit in lung cells is a signal from monocytes, lymphocytes, and other immune cells that activates IL-6 production.We propose that this model represents the basic set of interactions that settles on the inflammation reaction to the invasion of the lung cells by SARS-CoV-2, and that it is a first approximation to understand the main dynamical features of this response.
Our hypothesis is that Nsp5 leads to the sustained overproduction of IL-6, Cox2, PGE2, and NF B, in the epithelial lung cells, which is a necessary but not a sufficient condition for a cytokine storm.
In the particular case of epithelial lung cells, we found that the computational and mathematical analysis of the model supports our hypothesis, and that Nsp5 effectively enhances NF B, IL-6, Cox2 and PGE2 production during the time that remains bounded to HDAC2.Furthermore, the model predicts the existence of an over stimulated immune state (OSIS) in which the amount of IL-6 produced by the infected lung cells is high, and such state becomes autonomous and robust, independent of any stimulus from external sources like monocytes, and lymphocytes.Hence this OSIS cannot be completely shut down by classical anti-inflammatory drugs like Nimesulide and Dexamethasone, and can be possibly beaten only with specific inhibitors of Nsp5 like Saquinavir (Xu et al., 2020).

Model
Immune response to infections involves complex spatio-temporal processes involving multiple molecular components, non-linear interactions and cell-cell interactions.In consequence, the immune response to infections is a multidimensional, multiscale and multilevel dynamical process in which innate and acquired immunity mechanisms are involved to protect the organism against pathogens (Eftimie et al, 2016).In this scenario, the number of variables that become both sufficient and necessary for a detailed enough description of the spatio-temporal dynamics of the immune system is tremendously high.Some computational integrative efforts have been made in this direction (Danos et al., 2007), but imply very large systems that cannot be modeled with continuous quantitative models that are, in turn, necessary if we want to understand which are key nodes driving the system towards the OSIS.The lack of quantitative experimental data and knowledge of the values of most of the model parameters are additional challenges to achieve such dynamic models of the immune system in response to specific infectious diseases.Thus, most models of the immune system have been qualitative (Eftimie et al., 2016;Martínez-Sánchez et al., 2018).
A possible approach implies the use of ordinary differential equations (ODEs) qualitative modeling that is oriented to determine all the possible trajectories of a dynamical system along its ndimensional phase space (with n > 0), when the system is subject to a set of initial conditions and parameters values.High dimensional nonlinear dynamical systems can exhibit a variety of coexisting dynamical behaviors that determine the structure of their phase space (Strogatz, 2015).However, the structure of the phase space of biological nonlinear dynamical systems is highly dependent on the set of parameters values.Variations in the value of one or more parameters can drive a drastic change in phase space generating a different qualitative dynamics of the system.This process is known as bifurcation, and it is important to find the possible bifurcations in the system and to identify which parameter or parameters are the responsible of these changes (Strogatz, 2015)(See Supplementary Material 2 (SM 2)).
In order to avoid the problem due to the high dimensionality and complexity of the immune response to SARS-CoV-2 infection, and to use the minimum number of parameters in the model that we are proposing in this work, we assume that the interaction between nodes IL-6, NF B, IB, Cox2 and PGE2 settles on the main dynamical features of the inflammation process in lung epithelial cells due to the presence of the Nsp5 viral protein.In the proposed circuit shown in Figure 2, IL-6, NF B, Cox2 and PGE2 define three feedback loops (Figures 3A-3B): a positive feedback loop between NF B and IL-6 (Figure 3A); a positive feedback loop between NF B and IB (Figure 3A); and a negative feedback loop between PGE2 and IB (Figure 3B).
These feedback loops are distributed between three compartments: a) IL-6 is secreted into the extracellular medium where it acts in paracrine and exocrine form.IL-6 binds to its receptor IL6R at the cell surface.Once the complex IL6-IL6R is formed, it activates IKK in the cytoplasm (Wang et al, 2003).PGE2 is also secreted into this compartment where it binds to EP2 and EP4 receptors at the cell surface, activating Protein Kinase B (Akt) in the cytoplasm (Cho et al., 2014).
b) IKK, free NF B, and the complex IB-NF B are located in the cytoplasm, where IB is degraded by ubiquitination.Protein Cox2 is located in the endoplasmic reticulum (ER).
c) Free NF B enters the nucleus where induce the transcription of IL6, Cox2, IKB and NF B genes.Transcriptional activity of nuclear NF B is inhibited by HDAC2, and enhanced by Nsp5.
In this form, the model cell consists of three compartments: extracellular medium, cytoplasm and nucleus.We assume that the ratio of the cytoplasmatic volume to external volume (Vcyt / Vout ) is one, and that the radio of the cytoplasmatic volume to the nuclear volume (Vcyt / Vnuc) is 2. In each compartment the concentration of the molecules are measured in pg mL -1 .
In the circuit of Figure 2, the rate of production of free cytoplasmic NF B depends on the rate of its release from the IB-NF B complex, which is proportional to the product of the concentrations of IB and IKK (Figure 3A), its own rate of production due to nuclear NF B transcriptional activity inhibited by HDAC2 (Figures 1 and 3A), the rate at which is recaptured by IB (Figure 3A), and the rate of transport into the nucleus: where nf is the concentration of free cytoplasmic NF B, max nf V is the maximum rate of translation of the NF B gene,

 
on nf p t is the probability that the NF B gene is in its active state (on), ikk is the concentration of IKK, ibnf is the concentration of IB-NF B complex, ka and k2 are rate constants and D is the transport constant (Table 1). .
We assume that the concentration of IKK activated by IL-6 is proportional to the external concentration of the cytokine: where 1 is a constant of proportionality.In a similar form, PGE2 increases the rate of cleavage of the complex IB-NF B by phosphorylation of the IKK subunit by Akt (Bai et al, 2009)(Figure 3B).Thus, we assume that the concentration of Akt is proportional to the external concentration of PGE2: where k1 and k 11 are the new rate constants (Table 1).
The rate at which the complex IB-NF B is broken in presence of IKK depends on the rate at which IKK and Akt remove IB from the complex under the action of IL-6 and PGE2, and the rate at which the complex is formed again.: where ibnf is the concentration of the complex IB-NF B, and k4 and k5 are rate constants (Table 1).
The rate at which NF B enter the nucleus depends on the rate of transport of the molecule into de nucleus, adjusted for the change in volume between the cytoplasmic and nuclear compartments, and on the rate of inhibition of nuclear NK B by its inhibitors: where nf* is the nuclear concentration of NF B and kdeg is a rate constant.
The probabilities of activation of genes IL6, Cox2, IB and NF B (Figures 3A and 3B) are given by: see Supplementary Material I (SM I) for the mathematical deduction of these equations.For Eq. 5-7, the variable , is the probability that the gene i is expressed at time t, hd* is the amount of HDCA2 in nucleus, and k2, k-2, k3, k-3, k7, k-7, k15, k-15 are rate constants (Table 1).
We assume that free cytoplasmic HDAC2 is produced according to the rate equation: where hd is the amount of free cytoplasmic HDAC2 that is produced at a constant rate r, D2 is the rate of transport of free HDAC2 into the nucleus.In this equation, nsp5 is the amount of Nsp5 in the infected cell that sequester HDAC2 in cytoplasm (El Baba & Herbein, 2020), and  is a constant (Table 1).Free cytoplasmic HDCA2 enters the nucleus according to the equation: where hd* is the amount of HDAC2 in the nucleus, and kdeg2 is a rate constant (Table 1).
Nsp5 is a protein with only one link in the virus interactome (Gordon et al., 2020;Díaz. 2020a;Hekman et al, 2020), which is an input to the circuit of Figure 2.However, the real kinetic mechanism of production of this viral protein is unknown, thus we assume a plausible simple kinetic mechanism that leads Nsp5 to have a maximum constant concentration of 20 nM.We set the parameters  and k18 to the values shown in Table 1: where  is the constant rate of production of Nsp5, and k18 is a rate constant (Table 1).The effect of the variation of parameter  on the dynamics of the circuit of Figure 2 is analyzed in SM 2.
During the early stage of infection by SARS-CoV-2, the rate of variation of the concentration of IL-6 outside the lung cell depends on the rate of production of IL-6 by monocytes, leucocytes and other cells under stimulation by Tall-like Receptors (TLRs) (denoted by (t)) (Magro, 2020;Jafarzadeh et al., 2020), on the rate of production by lung cells due to nuclear NF B transcriptional activity enhanced by Nsp5 (Figure 3A), and on the rate of degradation of IL-6 by different mechanisms.
Finally, we assume that the concentration of IKK activated by IL-6 is proportional to the external concentration of the cytokine (see Eq. 2): where max 6 IL V is the maximum rate of translation of IL6,   6 on IL p t is the probability that the IL6 gene is activated at time t, and k5 and k6 are rate constants (Table 1).
The rate of variation of the amount of free IB in the system is the balance between the rate of translation of gene IB and the rate of formation of new IB-NF B complexes: where max I B V  is the maximum rate of translation of IB, and is the probability that the gene IB is activated at time t (Table 1).
The rate of variation of the amount of the enzyme Cox2 in the system depends on the rate of translation of the Cox2 gene, and its rate of inhibition by different mechanisms: where cox2 is the amount of the enzyme in the system, max 2 Cox V is the maximum rate of variation,   2 on Cox p t is the probability that the gene Cox2 is activated at time t, k8, k9 and k10 are rate constants (Table 1).
The rate of variation of the amount of PGE2 in the system depends on the amount of Cox2, on the amount of IL-6, and on its rate of degradation: where k11, k12 and k17 are rate constants (Table 1).In the reference state the translation velocities are: max In this section we present in detail and focus on the results that are relevant to understand the role of Nsp5 in the process of acute lung inflammation.The rest of the analyses and results related to the mathematical aspects of the model are reported in SM 2.
We solved the model using the Euler predictor-corrector method with time step of 0.05 s and 720,000 integration steps to simulate a period of 10 hrs post infection, which is the estimated time for the first bursting of new virions into the extracellular medium (Bar-On et al., 2020;Kumar et al., 2020).We used the parameter values shown in Table 1 to generate an arbitrary reference state (ARS) of the system for the initial conditions ibnf = 5 pg mL -1 , and the rest of the variables equal to 0 at time t = 0.In this case, the input (t) = 10 pg mL -1 s -1 and nsp5 = 0 for   0, t   .In this ARS, the translation velocities of all genes are set to 2 pg mL -1 s -1 (Hauser et al., 2019).Figure 4A shows that the output of the model is the transitory activation of NF B, Cox2 and PGE2, in contrast with a sustained high activation of IL-6.In order to clarify the source of the steady activation of IL-6, we solved the model with the same initial conditions and (t) = 0.In this case, the system does not turn on during all the time of simulation (data not shown), indicating that the signal from monocytes and other cells is a necessary condition to start the IL-6 sustained production in lung cells (Figure 4A). Figure 4B shows that when nsp5  0, (t) = 0 and equal gene translation rates, the response of the system is a high sustained level of activation of NF B, Cox2, and PGE2 with respect to the ARS.In this case, the concentration of IL-6 in the cell was only slightly affected by the virus.Figures from 4A to 4D show the first five minutes of simulation during which the circuit of Figure 2 reaches its steady state.This result indicates that the presence of the virus induce a sustained and persistent production of NF B, Cox2, and PGE2 in the circuit of Figure 2. In absence of the viral protein Nsp5, the response to an external input a(t) = 10 pg mL -1 s -1 is a transitory activation of NF B, Cox2, and PGE2.In this case, IL-6 shows a steady high concentration.In this case, all the gene translation rates are set to 2 pg mL -1 s -1 .B) In presence of Nsp5 the circuit shows steady higher concentrations of NF B, Cox2, and PGE2 (red line) with respect to ARS (black line).IL-6 steady concentration (black line) is unaffected by the presence of the viral protein (red line).In this panel, a(t) = 10 pg mL -1 s -1 , and all the gene translation rates are set to 2 pg mL -1 s -1 .C) When NF B and IB translation rates are max 10 nf V  pg mL -1 s -1 and max 15

Preprints
pg mL -1 s -1 , respectively, and Cox2 and IL-6 translation rates are 2 nM s -1 the circuit shows a weak sustained activation of NF B, Cox2, and PGE2 (black line) in response to the external input (t) = 10 pg mL -1 s -1 and in absence of Nsp5.In this case, IL-6 steady concentration does not change (black line).In presence of Nsp5 the circuit shows steady higher concentrations of NF B, Cox2, and PGE2 (red line), although the concentration of IL-6 is lower than is absence of the virus.D) In the over stimulated immune state (OSIS), the presence of Nsp5 produces an enhanced production of NF B, Cox 2, PGE2 and IL-5 (red lines).This particular state is obtained when (t) = 10 pg mL -1 s -1 , max 10 nf V  pg mL -1 s -1 , max 15 Translation rates of the set of genes of the model rarely have the same value, even when their activation is promoted by the same transcription factor (NF B in this case).As shown in Table 1, the results of Figures 4A and 4B are obtained when the rates of translation are the same for all genes.However, in real cells this is not the case, and Cox2, Il-6, NF B and IB are translated at different rates.In Homo sapiens the estimated average rate of translation ( max V ) of a gene is about 10 4 proteins per hour (Hausser et al., 2019).However, this average value can be increased according to the state of activity of the cell.In the present work, we arbitrary chose the value of max 1 1 2 pg mL s V    for all genes in the ARF (Table 1).However, we made parameter variation to know how different values of this set of parameters affect the qualitative behavior of the model (SM 2). Figure 4C shows that when in the system max 1 1 10 pg mL s nf V  


, max 1 1 15 pg mL s , and max 1 1 6 2 pg mL s , (t) = 10 pg mL -1 s -1 and nsp5 = 0, the circuit tends to a steady state in which NF B, Cox2, and PGE2 exhibit a constant low concentration during the time of simulation, while IL-6 steady concentration is unaffected with respect to the value shown in Figure 4A.When nsp5  0, the steady concentration of NF B, and Cox2, are increased ~ 10 fold with respect to their maximum concentration in the ARS, while PGE2 concentration is increased ~ 4 fold; Nsp5 produces a ~ 4 fold decrease in IL-6 steady concentration with respect to the reference concentration shown in Figure 4A.This result indicates that a the difference in translation rates between NF B and IB (the intensity of the inhibition of NF B by IB, Figure 3A) in absence of the virus with respect to their values in ARS induces a low sustained production of NF B, Cox2, and PGE2 without affecting IL-6 production.However the presence of Nsp5 induces a high production of NF B and Cox2, and a moderate increase in PGE2 concentration that decreases IL-6 concentration, indicating that a higher rate of IB translation with respect to NF B translation can decrease the concentration of the cytokine IL-6 in presence of the virus probably due to a lower concentration of PGE2in the circuit (Figures 2 and 3A).


the circuit tends to a steady state in which NF B, Cox2 and PGE2 exhibit a constant low concentration during the time of simulation, while IL-6 is unaffected with respect to its steady concentration in Figure 4A.However, when nsp5  0 the dynamical behavior of the circuit changes because NF B and Cox2 show the same steady concentration that in Figure 4C, but PGE2 has a higher steady concentration.In this case, IL-6 shows an increase of ~ 2 fold with respect to the level shown in Figure 4C.Thus, a high rate of IL-6 translation can overcome the effect of an increased rate of IB translation in presence of the virus, allowing an overproduction of the cytokine.
All these results suggest an important dynamical property of the circuit shown in Figure 2: a high value of the parameter max 6 IL V produces an over production of NF B, Cox2, PGE2, and IL-6 in the circuit in presence of Nsp5 with respect to the ARS (Figures 4A and 4D).We name over stimulated immune state (OSIS) to the state of the system shown in Figure 4D, for which: = 2 pg mL -1 s -1 , max 6 IL V = 20 pg mL -1 s -1 , nsp5  0 and (t) = 10 pg mL -1 s -1 .
Figure 5 shows the response of the circuit to a step function defined as:   10 0 10 0 10 The latter simulates an external signal of IL-6 initially generated at a rate of 10 pg mL -1 s -1 during 10 minutes, and switched off after the 10 minutes.Figure 5A is the response of the system in ARS to the step signal, in absence and presence of Nsp5.Interestingly, in the absence of Nsp5, the circuit exhibits a low concentration of NF B, Cox2 and PGE2, and a high concentration of IL-6, all which are zero when the external signal is switched off.In presence of Nsp5 the ARS is perturbed by (t), and reaches a new steady state with high concentrations of NF B, Cox2, and PGE2 that remain even when the external signal is switched off.Furthermore, the high IL-6 concentration reached when the circuit is turned on falls to a lower concentration different from zero when (t) = 0.
Figure 5B is the response of the circuit in OSIS to the step function.In absence of Nsp5, the circuit exhibits a low concentration of NF B, Cox2 and PGE2, and a high concentration of IL-6, all which are zero when the external signal is switched off.However, in presence of Nsp5 the OSIS is perturbed by (t), and reaches a new steady state with high concentrations of NF B, Cox2, and PGE2 that remain even when the external signal is switched off.Furthermore, the higher IL-6 concentration reached when the circuit is turned on only falls ~ 50% of its maximum concentration when (t) = 0.These results suggest that the OSIS becomes autonomous of the external signal generated by monocytes and other immune cells in presence of Nsp5, i.e., the OSIS becomes a persistent state in presence of SARS-CoV-2, and never shut down by itself or by removing external signaling during the 10 hrs of simulation.16) of main text and is applied to the ARS in absence of Nsp5 (black line) the level of cytokine IL-6 initially increases to a maximum value of 20 pg mL -1 , and becomes zero when the external signal is switched off.There is not any effect on NF B, Cox2 and PGE2 concentrations in the circuit.However, in presence of Nsp5 (red line) the levels of NF B, Cox2, PGE2 and IL-6 are increased in response to (t) and never become zero when the external signal is switched off.B) When the step function of Eq. ( 16) is applied to the circuit in absence of Nsp5 (black line) with max 10 nf V  pg mL -1 s -1 , max 15 and max 6 20 IL V  pg mL -1 s -1 , there is a low increase in NF B, Cox2 and PGE2 concentrations that becomes zero when the external signal is switched off.IL-6 initially increases to a maximum value of 20 pg mL -1 , and becomes zero when the external signal is also zero.In contrast, when the circuit is in OSIS in presence of Nsp5 (red line) the concentrations of NF B, Cox2, PGE2 and IL-6 are increased in response to (t) and remain at a high value even the external signal is switched off.
In order to reduce the effect of Nsp5 on the circuit of Figure 2, we tested the effect of Nimesulide on NF B production.Nimesulide, a Cox2 inhibitor, acts inhibiting the synthesis of Cox2 and lipooxigenase enzyme and their products (Suleyman et al., 2008).We modified Eq. ( 14) to introduce an inhibition term that decreases the rate of production of Cox2 in presence of Nimesulide: where nim is the amount of Nimesulide at time t, whose rate of variation of its concentration in the cell is: where r2 = 40 pg mL -1 s -1 and k19 = 3 s -1 .Solution of the modified model (Eq.17) is presented in Figure 6.In Figure 6A, Nimesulide is applied to the circuit in presence of Nsp5, and with all gene translation rates set to 2 pg mL -1 s -1 .In this case, Nimesulide has no effect on NF B, a slightly effect on IL-6, and a significant decrease in Cox2 and PGE2 concentrations to values near to zero in presence of 14 pg mL -1 of the drug.The same effect is observed when the circuit is in OSIS (Figure 6B).Thus, Nimesulide is not an effective anti-inflammatory agent in the context of the dynamic circuit modeled here, and does not impact the concentration of IL-6 in the presence of Nsp5, even at high concentrations.
Dexamethasone (DX), is a potent synthetic glucocorticoid that inhibits NF B transcription, and is another anti-inflammatory drug that is being used to ameliorate acute inflammatory states due to SARS-CoV-2 infection.In order to test the effect of DX on the inflammatory process in the circuit of Figure 2, we modified Eq. ( 8) to introduce an inhibition term that decreases the probability of activation of NF B: where DX is the amount of DX at time t, whose rate of variation of its concentration in the cell is: where r3 = 30 pg mL -1 s -1 and k20 = 3 s -1 .
Figure 6.-Effect of Nimesulide on lung cell inflammation.A) Effect of 14 pg mL -1 of Nimesulide on the levels of NF B, Cox2, PGE2 and IL-6 when (t) = 10 pg mL -1 s -1 , and all gene translation rates are 2 pg mL -1 s -1 in presence of Nsp5 (red line).Nimesulide decreases the concentration of Cox2 and PGE2, but has not effect on NF B and IL-6 (blue line).B) Effect of 14 pg mL -1 of Nimesulide on the OSIS (red line).Nimesulide significantly decreases the concentration of Cox2 and PGE2, but has not effect on NF B, The drug has a weak effect on IL-6 concentration (blue line).
Figure 7 shows the effect of DX on the circuit.Figure 7A shows that when DX is applied to the circuit in presence of Nsp5, and when all genes have the same translation rate of 2 pg mL -1 s -1 , the drug eliminates NF B, Cox2 and PGE2 from the circuit, but has a null effect on the concentration of IL-6.However, when the circuit is in OSIS (Figure 7B), DX has little effect on the immune response in presence of Nsp5. Figure 7C shows that when all translation rates are the same (2 pg mL -1 s -1 ), DX can shut down the immune response in presence of the virus, when the circuit input (t) is the function of Eq. ( 16).However, DX cannot shut down the effects of Nsp5 when the circuit is in OSIS, having a little anti-inflammatory effect (Figure 7D).Hence, our results suggest that DX is an effective anti-inflammatory drug during mild inflammatory states, but is less effective during critical acute lung inflammation due to SARS-CoV-2 (Horby et al., 2021).This dynamical behavior of the circuit of Figure 2 is unaffected by higher doses of DX.When DX and Nimesulide are applied together (with 35 pg mL -1 of DX and 40 pg mL -1 of Nimesulide) to the circuit in OSIS, in which (t) is the function of Eq. ( 16), there is also a limited anti-inflammatory effect (Figure 8A).16) of main text and is applied to the circuit in presence of Nsp5 (red line), with all gene translation rates equal to 2 pg mL -1 s -1 , 14 pg mL -1 of DX decreases the concentration of NF B, Cox2, PGE2 and IL-6 to zero (black line).D) When the external signal a(t) is the step function of Eq. ( 16) of main text and 14 pg mL -1 of DX are applied to the circuit in the OSIS, the drug has a practically null effect on the inflammation process.16), there is also a little anti-inflammatory effect.B) The molecular docking of Nsp5 completely shuts down the inflammation process when the circuit is in OSIS, with (t) given by Eq. ( 16), suggesting that Nsp5 is effectively the cause of the acute lung inflammation during SARS-CoV-2 infection.
In contrast, the predicted limited effect of Nimesulide and DX on reverting the OSIS, the molecular docking of Nsp5 with a hypothetical drug like Saquinavir (Xu et al., 2020), is able to completely eliminate inflammation.We tested this hypothesis modifying Nsp5 production in Eq. ( 11) with an inhibitory term: where drug is the amount of the hypothetical drug that can remove Nsp5 from the circuit, whose concentration at time t is given by: with r4 = 30 pg mL -1 s -1 and k21 = 3 s -1 .
Figure 8B shows that when the hypothetical drug is applied to the circuit in OSIS, with (t) given by Eq. ( 16), the inflammatory process is completely shut down, suggesting that Nsp5 is effectively the key cause of the acute lung inflammation during SARS-CoV-2 infection, and a fundamental potential target for its treatment.

Discussion
The complexity of the immune response to infections implies a complex set of different kinds of cells that involve complex regulatory networks and molecules.Together, they underlie specific responses that may lead to resolution or to inflammatory responses whose dynamical properties depend on the networks structure (i.e., the number of nodes, links, and the distribution of links between the nodes) occurring at different scales and levels (Eftimie et al., 2016;Morel et al., 2006).
The structure of such complex and dynamical network is not random due to the specific nature of its nodes and links, so that if the network has a hierarchical structure its hubs determine the overall network dynamics and the structure of its phase space (Eftimie et al., 2016;Díaz. 2020b;Cessac, 2009).In high dimensional phase spaces, steady points of diverse kinds can coexist with isolated closed curves (limit cycles), strange attractors of fractal structure, and other limit sets giving rise to a variety of biological dynamical behaviors like molecular switches, periodic and quasi-periodic oscillations, and bursting (among others).Two central problems in the analysis of this kind of dynamical systems are: a) the structure of the phase space is strongly dependent of the value of the set of parameters of the dynamical system, and b) the nonlinearity characteristic of biological systems such as the one treated here (see SM 2) (Nayak et al, 2018).
In most of the ODEs-based biological models so far studied, there is limited knowledge of the real values of the parameters, and this leads to uncertainty concerning the extent to which the model explains the real process under study.In this case, the validation of the model relied on experimental data.Nonetheless, with the exception of some particular cases (Hodgkin and Huxley, 1952), qualitative models have to be proposed and quantitative validations remain ahead (Eftimie et al., 2016).This challenge increases with the dimensionality of the model.Thus, is necessary to choose some particular properties of the network under study (inflammation network in this case) that incorporate key dynamical features of the whole system under study, and that can also aid to reduce the number of parameters of the model to the minimum possible.
In this sense, we limited the scope of this work to the study of the inflammation process of epithelial lung cells during SARS-CoV-2 infection.The epithelial lung cells form a sub-network that have two highly connected nodes: the internal signal integrator NF B, and the external signal integrator IL-6 Receptor (IL6R) (Magro, 2020), which is activated by external IL-6 (Figures 2 and  3).Both integrators trigger the activation of two important intermediate processes: Cox2 and PGE2 synthesis.Thus, we propose, as a first approximation, that the dynamical behavior of this these four nodes of the circuit of Figure 2 reflect with certain accuracy the dynamical properties of the complex inflammation sub-network of lung cells.Thus, the dynamics of IL-6, NF B, Cox2 and PGE2 nodes is described by a model that consist of 12 nonlinear coupled ordinary differential equations (Table 1), that we assume represents the central dynamical core involved in the lung cells inflammation process in response to SarsCoV-2 infection.Such core includes three coupled feedback loops shown in Figure 3.We used this low-dimensionality model to explore the effect of Nsp5 in the qualitative dynamics of the inflammation circuit of Figure 2 (Diaz, 2020b;Kumar et al, 2020).
SARS-CoV-2 main protease Nsp5 has a role as an epigenetic regulator of host DNA expression in lung cells; and thus its physical association with HADC2 (Gordon et al., 2020;Kumar et al., 2020) avoids the movement of the deacetylase from the cytoplasm to the nucleus, increasing the probability of the binding of NF B to the promoter site of its target genes (Eqs.5-8) (El Baba & Herbein, 2020).NF B is a central integrator of the signals that initiate the inflammation process (Figure 1), which include IL-6 secreted by monocytes, leucocytes, macrophages, among other cells, together with signals from TLR4 and TLR7 that suppress the IB-NF B complex in a MyD88 independent manner (Vidya et al., 2017) .NF B enhanced transcriptional activity includes the transcription of IL-6, Cox2, IB and NF B itself that are key genes in the process of inflammation (Figure 3).
Results obtained from our model suggest that Nsp5 can effectively transform a weak inflammation response into a persistent sustained one that has a relative high concentration of the cytokine IL-6 (Figure 4B).In the example of Figure 4B, all genes have the same translation rate, a condition that have low probably of occurrence in the real human immune system in which every gene is subject to different processes of post transcriptional regulation (Hausser et al., 2019).When this restriction is removed, and IL-6 has a high translation rate, Nsp5 boosts the concentration of the four main proteins of the circuit with respect to their respective values in ARS, and the inflammation response becomes strongly persistent with a high IL-6 concentration.Furthermore, in this state that we named OSIS, the dynamics of the circuit of Figure 2 becomes independent of any kind of external signal and never turns off, leading to the deregulation of the production of the cytokine IL-6 (Figure 5B).In this form, Nsp5 changes the qualitative dynamical behavior of the system, in which the trajectories that initially span around the ARS stable node now span around the OSIS stable node (see SM 2) (Díaz, 2020b).
An important biological consequence of this change in the qualitative dynamics of the circuit is that ACE2 receptor synthesis also becomes persistent and autonomous of any external signal to the lung cells (Figure 2).Furthermore, neither Nimesulide nor Dexamethasone or both can eliminate the OSIS suggesting that this dynamical behavior could be the form in which the virus assures the persistence of its reproductive cycle without any perturbation from the natural defenses of the body.This can be also a possible clue about the possible adaptive substitution of N protein used in SARS-CoV for the use of the main protease Nsp5 in SARS-CoV-2 as the switch for the inflammation process, which is necessary for ACE2 sustained production, taking into consideration that SARS-CoV N protein targets Cox2 while SARS-CoV-2 Nsp5 protein targets the master integrator NF B.
If this hypothesis is true, this is the probable cause of the increased pathogenicity of SARS-CoV-2 with respect to SARS-CoV.
It is of interest that the persistence of the OSIS is supported by a high translational rate of IL-6, which in the model must have a value of max 6 7 IL V  pg mL -1 s -1 .This result pictures a probable hypothetical scenario where the external signal (t) triggers the initial inflammation response with a low translation rate of IL-6 that can be shut down with DX.As viral infection continues, IL-6, NF B, and IB increase their rate of translation ( max 2 Cox V value has a little weight in this process) and the dynamical behavior of the circuit of Figure 2 becomes independent of the value of (t), and selfsustained, i.e., the control of the production of the cytokine IL-6 is translated from the external immune cells to the lung cells, and cannot be stopped even with high doses of DX.In this form, in the model, max In this extreme situation (OSIS), the unique treatment possible suggested by our model is an inhibitor of Nsp5 like Saquinavir (Xun et al., 2020).In this case, the OSIS is completely shut down, which indirectly implies the down production of ACE2.Saquinavir is an anti-retroviral drug used against the main protease of HIV-1 (Bensussen et al., 2018;Xun et al., 2020), with some undesirable effects like diarrhea, abdominal pain, and nausea.We suggest that this drug could be a plausible treatment against the effects of SARS-CoV-2 in acute lung inflammation, and as a complement to the new vaccine against this coronavirus.
In a previous work (Díaz, 2020a) the role of the viral protein Orf8 as a hub of the SARS-CoV-2 infection was suggested.This highly connected protein has as main targets those processes related to the vesicular trafficking required for the ensemble of new virions (Figure 1) (Gordon et al., 2020).Thus, the inhibition of the effects of Orf8 in the host cells could block the reproduction cycle of the virus but not the inflammation process because Nsp5 is not directly linked to Orf8 (Gordon et al., 2020).Rapamycin has been suggested as a drug against the effects of Orf8 but produces severe immunosuppressant effects, which, according to the results of our model, will not be of care in cases of severe lung inflammation because Nsp5 uncouple the circuit of Figure 2 from the process of viral replication (Figure 1 and Figure 5B).In this case, a treatment with Rapamycin and Saquinavir can be an alternative for patients with severe lung inflammation.This suggestion needs experimental verification and clinical evaluation.

Conclusions
In this work we propose a dynamical model that puts forward a novel hypothetical mechanisms underlying lung cells inflammation process in response to SARS-CoV-2.In this scenario the main protease Nsp5 effectively enhances the inflammatory process, increasing the levels of NF B, IL-6, Cox2, and PGE2 with respect to the reference state.When the translation rates of NF B and IB are increased to a high constant value, and the translation rate of IL-6 is increased above the threshold value of 7 pg mL -1 s -1 the circuit enters in a persistent over stimulated immune state (OSIS) with high levels of the cytokine IL-6.The OSIS never shuts down by itself, and becomes autonomous of the signals from other immune cells like macrophages and lymphocytes.Our model explains why DX or Nimesulide have little effect on the OSIS, and the only means to suppress such acute and sustained inflammation state is by inhibiting the protein Nsp5 with drugs such as Saquinavir.
In this form, the model suggests, in accordance to our hypothesis, that Nsp5 is effectively a key node underlying severe acute lung inflammation during SARS-CoV-2 infection.The persistent production of IL-6 by lung cells during infection can be one of the causes of the cytokine storm observed in critical patients with COVID19.From an evolutionary point of view, the use of Nsp5 as the switch to start inflammation, and the consequent overproduction of the ACE2 receptor, could have driven the increased pathogenicity of SARS-CoV-2 with respect to SARS-CoV.
The present theoretical work is part of the project "Dynamics of the SARS-CoV-2 network" and is based on the experimental data reported in the literature at the moment.

Figure 1 .
Figure 1.-Process of infection of lung cells by SARS-CoV-2.Infection begins when the protein S of the virion binds with high affinity to the cells surface receptor ACE2 (Angiotensin-Converting Enzyme 2) that is highly abundant in the lung alveolar type II cells.The formation of the complex S-ACE2 initiates the process of fusion between the virion envelope and the cells membrane leading to the liberation of the nucleocapside with the viral genome into the cytoplasm .SARS-CoV-2 genome consists of a positive-sense nonsegmented single stranded mRNA ((+)ssRNA) of ~ 30 kb.The open reading frames 1a (orf1a) and 1b (orf1b) are located near the 5'UTR of the (+)ssRNA and they code for the polyproteins pp1a and pp1ab.Maturation of (Figures1 and 2); suggesting that Nsp5 represses HDAC2 inhibitory action on NF B (Figure1), and thus enhances the transcription of the pro-inflammatory genes targeted by the Nuclear Factor B (NF B) (Figures1 and 2)(Wagner et al., 2015).Such role distinguishes SARS-CoV and SARS-CoV-2.The SARS-CoV virus uses N protein to promote the sustained transcription of the Cyclooxygenase 2 (Cox2) enzyme by binding directly to the NF-B transcriptional regulatory elements, and to the CCAAT/enhancer binding proteins of the gene Cox2(Yan et al., 2006).Hence, SARS-CoV and SARS-CoV-2 produce severe acute respiratory syndrome through different molecular mechanisms.

Figure 2 .
Figure 2.-Circuit representation of inflammation.IL-6, NF B, Cox2, and PGE2 are the basic nodes of the circuit that represents the basic features of inflammation of lung cells during SARS-CoV-2 infection.In this circuit the arrows represent activation and the bars inhibition.In the figure are shown two positive feedback loops between IL-6 and NF B, and between NF B and Cox2.IL-6 also induces the expression of the ACE2 receptor in the lung cell membrane.NF B is inhibited by HDAC2, which in turn is inhibited by the viral protein Nsp5.The circuit has two targets for the anti-inflammatory drugs Nimesulide and Dexamethasone.

Figure 3 .
Figure 3.-Feedback loops of the model of inflammation.The mathematical model proposed in this work is based on the dynamical features of three feedback loops: A) two positive feedback loops between IL-6 and NF B, and NF B and IB; B) a negative feedback loop between PGE2 and IB.In Figure A, (t) is a external signal that represents the amount of IL-6 produced by immune cells like monocytes and lymphocytes, among others.

Figure 4 .
Figure 4.-Effect of Nsp5 on IL-6 production.A) Arbitrary reference state (ARS) of the circuit of Figure 2.In absence of the viral protein Nsp5, the response to an external input a(t) = 10 pg mL -1 s -1 is a transitory activation of NF B, Cox2, and PGE2.In this case, IL-6 shows a steady high concentration.In this case, all the gene translation rates are set to 2 pg mL -1 s -1 .B) In presence of Nsp5 the circuit shows steady higher concentrations of NF B, Cox2, and PGE2 (red line) with respect to ARS (black line).IL-6 steady concentration (black line) is unaffected by the presence of the viral protein (red line).In this panel, a(t) = 10 pg mL -1 s -1 , and all the gene translation rates are set to 2 pg mL -1 s -1 .C) When NF B and IB translation

Figure 5 .
Figure 5.-Response of the circuit to a step function.A) When the external signal (t) is the step function of Eq. (16) of main text and is applied to the ARS in absence of Nsp5 (black line) the level of cytokine IL-6 initially increases to a maximum value of 20 pg mL -1 , and becomes zero when the external signal is switched off.There is not any effect on NF B, Cox2 and PGE2 concentrations in the circuit.However, in presence of

Figure 7 .
Figure7.-Effect of Dexamethasone on lung cell inflammation.A) Effect of 14 pg mL -1 of Dexamethasone (DX) on the levels of NF B, Cox2, PGE2 and IL-6 when a(t) = 10 pg mL -1 s -1 , and all gene translation rates are 2 pg mL -1 s -1 in presence of Nsp5 (red line).DX effectively decreases the concentration of NF B, Cox2 and PGE2 to cero, but has not effect on IL-6 (black line).B) Effect of 14 pg mL -1 of DX on the OSIS (red line).DX has a weak effect on the concentration of NF B, Cox2, PGE2, and IL-6 (red line).C) When external signal a(t) is the step function of Eq. (16) of main text and is applied to the circuit in presence of Nsp5 (red line), with all gene translation rates equal to 2 pg mL -1 s -1 , 14 pg mL -1 of DX decreases the concentration of NF B, Cox2, PGE2 and IL-6 to zero (black line).D) When the external signal a(t) is the step function of Eq. (16) of main text and 14 pg mL -1 of DX are applied to the circuit in the OSIS, the drug has a practically null effect on the inflammation process.

Figure 8 .
Figure 8.-Effect of the molecular docking of Nsp5 on lung inflammation.A) When DX and Nimesulide are applied together (with 35 pg mL -1 of DX and 40 pg mL -1 of Nimesulide) to the circuit in OSIS (red line), in which (t) is the function of Eq. (16), there is also a little anti-inflammatory

Table 1 .
-Equations and parameters of the basic model of inflammation