key: cord-0911410-2s1npse5
authors: Du, Sean Quan; Yuan, Weiming
title: Mathematical modeling of interaction between innate and adaptive immune responses in COVID‐19 and implications for viral pathogenesis
date: 2020-05-13
journal: J Med Virol
DOI: 10.1002/jmv.25866
sha: 6831780ac3a40e59bef5b1fd9b7b59b34702b8a6
doc_id: 911410
cord_uid: 2s1npse5

We have applied mathematical modeling to investigate the infections of the ongoing coronavirus disease‐2019 (COVID‐19) pandemic caused by SARS‐CoV‐2 virus. We first validated our model using the well‐studied influenza viruses and then compared the pathogenesis processes between the two viruses. The interaction between host innate and adaptive immune responses was found to be a potential cause for the higher severity and mortality in COVID‐19 patients. Specifically, the timing mismatch between the two immune responses has a major impact on disease progression. The adaptive immune response of the COVID‐19 patients is more likely to come before the peak of viral load, while the opposite is true for influenza patients. This difference in timing causes delayed depletion of vulnerable epithelial cells in the lungs in COVID‐19 patients while enhancing viral clearance in influenza patients. Stronger adaptive immunity in COVID‐19 patients can potentially lead to longer recovery time and more severe secondary complications. Based on our analysis, delaying the onset of adaptive immune responses during the early phase of infections may be a potential treatment option for high‐risk COVID‐19 patients. Suppressing the adaptive immune response temporarily and avoiding its interference with the innate immune response may allow the innate immunity to more efficiently clear the virus.

vere illness tend to have severe pneumonia, sometimes acute respiratory distress, which could lead to multiple organ failures and death. The hospital stay can be weeks for such patients. [1] [2] [3] [4] There is no proven cure or antiviral drugs available yet.

The influenza virus typically binds to the sialic acid receptors on the surface of epithelial cells in the upper respiratory tracks and the upper divisions of bronchi. 5 In very severe, and often fatal cases of influenza, the infection will spread to the lower lungs. 6, 7 Similarly, the SARS-CoV-2 virus binds to the angiotensin-converting enzyme 2 (ACE2) receptors on the surface epithelial cells. 8 These cells with ACE2 expression are hence vulnerable to virus attack and are called target cells in mathematical models. The density of the target cells has been found to vary significantly in different regions of the respiratory tracks, with the highest in the lungs, followed by the nose, and lastly, the trachea/bronchi tissues. [9] [10] [11] As such pneumonia is relatively common among COVID-19 patients.

The incubation time for influenza is typically 48 hours or less. 12 The eclipse phase is reported to be around 6 hours, with viral load peaks around 2 to 3 days postinfection (dpi). 5 The infection typically resolves in 3 to 5 days. In typical influenza patients, the adaptive immune responses (AIRs) including pathogen-specific antibodies (Abs) and CD8 + cytotoxic T lymphocytes (CTL) are first observed around 5 dpi, peaking around 7 dpi. 13 The mean incubation period for COVID-19 patients is estimated to be 6 days, with a range distribution of 13 days. 14, 15 The median time from symptoms to discharge from the hospital was about 22 days, whereas the median time to death was 18 days. 4 In this paper, we used mathematical modeling to investigate the dynamics of the viral infection/replication inside a human host, in particular, the influenza and the SARS-CoV-2 virus, as well as the interactions of target cells with the innate and AIRs. Our model suggests that most of the differences between the two types of infections can potentially be attributed to the timing mismatch between the two immune responses. More specifically, influenza is a very acute infection; all vulnerable cells are completely depleted and viruses are more or less cleared by the innate immune response, before the adaptive immune response (AIR), which has a transient nature, reaches a significant level. The arrival of AIR may seem futile but it helps to completely clear the remaining viruses. On the other hand, the peak level of the AIR of the COVID-19 patients is more likely to come before the number of infected cells reaches its peak, due to the slower disease progression outlined above. As such the AIR helps to eliminate more viruses in the rapid growth phase and slow down the infection and depletion of vulnerable epithelial cells in the lungs of COVID-19 patients. A strong AIR, in this case, may be bad for the patients, as the disease progression is extended longer while a weaker AIR can be beneficial for the patients due to the rapid depletion of the target cells. One recent observation puzzling for many is the apparent resurgence of the disease for some "cured"

patients. Based on our theory and modeling, this phenomenon can be explained by that of a strong interaction of innate and AIRs in a host, which can sometimes create double peaks in viral load, separated by weeks. In the second part of our study, we applied our modeling to the current bedside treatment regimen. Antiviral drugs and immunosuppressive drugs have been evaluated as potential treatments for high-risk patients. We proposed an immune-suppressing treatment based on the leanings of our modeling study, which is to apply immunosuppressive drugs during the early phase of infection to reduce the AIRs to a level low enough not to interfere with the innate immune response. Once the vulnerable cells are depleted, the drug should be removed to let the AIR run its course to clear the remaining virus. This new approach can potentially lead to fast and complete recovery. This counter-intuitive approach aims to mimic the excellent job human immune systems have done with the influenza infection. Further studies on immune responses of COVID-19 patients and validation of our core assumptions and theory are needed before this can be applied in real patients' situations.

Mathematical models have been developed to understand the dynamics of viral infections. [16] [17] [18] [19] [20] Most of them fall into the category of the target cell-limited model with some variations. The simplest version includes three species: uninfected susceptible target cells (T), which in our case is the surface epithelial cells with ACE2 receptors, located at the respiratory tracks including lungs, nasal and trachea/ bronchial tissues, infected virus-producing cells (I), and the virus particles (V). They can be described by the following set of differential equations: The ability of a virus to grow the infection or to be cleared is determined by the basic reproductive ratio R:

where R represents the number of infected cells newly generated by one infected cell. With R> 1 the infection grows, while with R< 1 infected cells decrease and the virus is cleared out. 21 Here it is important to note that T is a variable. When the target cells are depleted, T goes down, and so does R. This is a very simple yet critically important concept for understanding the virus clearance. Complete depletion of the target cells, namely = T 0, leads to = R 0, hence achieving the best virus clearing ability. R 0 is the native (2).

In their experiment, the authors found that the viral load peaked at 2 to 3 dpi. The level of + CD8 cells and IgM data shows a very narrow peak at about 8 days dpi, while the level of IgG antibody shows a relatively broad peak lasting well over 50 days. We seek to mimic these temporal profiles, using Δ ( ′) t E as an example:

, .

With such a profile, we can choose the peak position ′ t c , and the rising and falling width w r and w f in our simulation to see the effects of these parameters on the viral infection dynamics. Our main assumption is that the AIR is a short term effect, which assumes that the virus is successfully cleared timely and the patient survived. As such, it is not our goal to develop a detailed model for

AIRs, but rather to understand the necessary conditions that are needed for a successful outcome for the patients. Other components like IgA, CD4 + T cells, and a number of other components of the AIR could also play important roles here. 13 To analyze the effect of antiviral drugs that either block infection ϵ k and/or production of viral particles ϵ p , the target cell-limited model is modified as follows: refers to a drug that is 100% effective. With the action of the drug, the reproduction number R is revised as follows:

Finally, we want to add another component to describe the effect of a drug that suppresses the innate immune response,

for reasons that will be discussed in the section of treatments by immunosuppressive drugs. Equations (1), (2), and (3) are changed to the following:

c describes no effect of immunesuppressing drugs while δ ϵ = , 1 c refers to complete suppression.

Mathematical models of influenza infections have been studied extensively, 5,22-25 our interest is to use it as a starting point before we embark on more challenging work on SARS-CoV-2 infections.

We have carefully evaluated the parameter values tabulated by Beauchemin et al 5 and decided to use the values shown in Table 3 .

The choice of the parameters is based on cross-checking multiple sources, and testing many simulation runs with variations of parameters to ensure that the output agrees with what is known about influenza infections. It is not our intent to fit the model to any specific set of data, but rather to be able to reproduce the general features of the influenza infections. The result is shown at Figure 1 . With these parameters, the infected cell count peaks at 2.45 dpi; viral load peaks at 2.98 dpi, and and the R 0 is about 47. γτ , which describes an exponential growth rate (shown in Appendix) is 2.34. As we should expect, the target cells are depleted shortly following the peak of the infected cells, after which the virus count decays exponentially.

The situation with regard to the decay portion is more complicated, and the consideration of the AIR is warranted if one is to understand it beyond 6 dpi. Table 4 shows the parameters for the three components we considered in Equations (2).

These values are chosen so that the temporal profiles resemble what was presented by Miao et al 25 The choice of the peak time is based on Beauchemin et al 5 We should clarify that the temporal profiles are not completely the same as in reference, 25 as some other literature has shown broader profiles. 5 We have done sensitivity analysis by doing multiple simulations, and find that as long as the width is not too narrow, the results are not sensitive to the width or shape of the peak. On the other hand, the location of the peak time is very important for the complete clearing of the virus. The result is shown in Figure 2 .

In Figure 2 , we included the same graph for virus count without the AIR, as shown in Figure As long as the native R 0 is greater than 1, this pattern is the universal behavior of the Equations (1), namely, a rapid exponential rise followed by an exponential decay, then eventually reaching an equilibrium at a low level of infection. Changing parameters in Table 3 will change the time scale, the peak/trough levels and locations of virus count etc., but the general pattern is the same. When the AIR is included, the virus clears out much faster. By day 15, the virus is completely gone. Mathematically virus and infected cells count never actually go to zero. However, we know that the infected cell count has to be an integer. Therefore when this number is below 1, one can safely consider it game-over for the virus. In our modeling, we set this cutoff to be 0.2 infected cell to be conservative.

The temporal profile of the three components of the AIR are also shown in Figure 2 so that the readers can see their peak positions as compared to the actual virus counts over time. After the AIR becomes active, the decay of virus count becomes much faster. We also tried to vary the parameters around what is listed in 

Due to the lack of experimental data, the discussions below are based on our best assumptions and observations. When shifting gears from the influenza virus to the novel corona virus, the most important difference we noticed is the difference in the pace of infection progression. Multiple reports suggested that the viral counts of SARS-CoV-2 do not reach a peak until 2 to 3 weeks postinoculation. [2] [3] [4] Here we will make a bold assumption that, if without the AIR, the peak time is around 8 to 14 days, if not more. As most of the severe cases result in pneumonia and the lungs are significantly larger than the nose, T 0 is estimated to be higher by a factor of 10 4 . Table 5 lists the only two parameters that are changed from Tables 3 and 4, together with the three dependent variables as outputs from the model. Here we aimed to have a peak time for infected cells at around 12 dpi if considering the effect of innate immune response only. Together with a 10 4 times increase of T 0 , these constraints essentially limited our choice for the only other variable to a very narrow range. The result is shown in Figure 3 .

As we can see in Figure 3 , the curve of virus count for an innate immune response only (no AIR) is similar to that of influenza infec- Our initial choice of AIR activity of 0.3 is purely artificial at this point, so we want to explore the effect of different levels of AIR activity. The result is shown in Figure 4 .

As we can see, with the AIR activity level of 0.2-0.4, the virus count reaches a high level, and then the growth is slowed with the activation of the AIR at around 8 dpi. After the effect of the AIR wanes, the virus count continues its ascent, reaching its peak at a later date. At a higher AIR activity level, we start to see an apparent drop in virus count after the onset of AIR, reaching a trough then rising to its peak at a much later day. Overall, an increase in the AIR activity tends to increase the duration of viral activity extensively in this case. This is because early activation of the AIR helps to reduce infection and depletion of target cells before it reaches a peak level, which saves a large number of uninfected target cells for later infection. When the effect of the AIR wanes, the virus takes control again when the effective R rises above 1. Whether the virus count grows or drops during the peak AIR activity is dependent on how strong the AIR effect is, and whether the effective reproductive ratio R at the time is greater or less in Figure 4 with the AIR activity level greater or equal to 0.8 offers a ready explanation for some perplexing clinical observations that some patients can appear to be recovered, but have the virus level resurges at a later time.

What's even more interesting is when we shift the day that the AIR activity peaks, as shown in Figure 5 .

To show the AIR effect on viral replication in a more pronounced fashion, we fix the AIR activity level at 0.5, and vary the day of its peak With additional simulations, we have found that if the AIR peak is after the viral peak, a higher level of AIR activity always helps to clear the virus faster. Contrarily, if the AIR peaks before the viral peak, higher AIR activity can increase the duration of viral activity. We also found that if the AIR level is low enough, for example below 0.1, its effect becomes negligible and the virus replication and clearing process is In essence, our modeling aims to explore potential means to manage the COVID-19 infections so that it can be similar to influenza infections. When the virus can be quickly and completely cleared, the AIR should be transient, similar to the influenza infections. Our modeling then establishes that one of the requirements is that the peak of AIR needs to come after the peak of infected cells, or equivalently, the peak of viral load. This can be viewed as a necessary condition for fast and complete clearing of virus. To achieve rapid clearing of the virus, it is also desirable to have the target cells more or less completely depleted. In reality, when the virus cannot be cleared timely, it is likely that the AIR will also be extended longer, so the conclusions from our modeling related to the later stage development may need to be modified.

When searching for a cure for COVID-19 patients, an effective antiviral drug is highly anticipated, so we decided to model the effect of the antiviral drug on the corona virus infection. We started our base case with a peak of the AIR at 8 dpi, and AIR activity level of 0.3, as shown before. From Equations (3), the expression for R is revised as follows:

Here R 0 is the native reproductive ratio; T T 0 is the percentage of the target cells that are not depleted. The first two brackets show the effect of the antiviral drugs on the reproductive ratio R, which shows that the drug has an equivalent effect regardless of whether it is acting on k or p. Next, we tried to simulate the effect of an antiviral F I G U R E 5 Effect of adaptive immuneresponse (AIR) peak day on viral dynamics: data shows the viral load over time when the day of AIR peak is changed from 6 to 14 days drug with 50% efficacy on k, namely, ϵ = 0.5 k , taken for 15 days consecutively. Taking this drug for a longer time does not necessarily make it more beneficial in this particular scenario. As a comparison, the base scenario is also included with a native R 0 of 7.03, and AIR activity of 0.3 peaking at 8 dpi. The base scenario is shown as a blue line in Figure 6 . When the drug is applied, the effective R is reduced to 3.52, according to the formula above, not considering the effect of target cell depletion.

Similar to the AIR effects shown in Figure 5 , Figure 6 shows that additional drug-induced antiviral activity has a pronounced effect on the duration of viral activity within the host for the three scenarios when the antiviral drug is started on 5 dpi, 10 dpi, or 15 dpi. During the 15 days when the antiviral drug is effective, the virus growth rate is reduced. After the drug is removed, the virus either resumes its growth or remains at a plateau. In all three scenarios, the duration of virus activity is extended substantially, which can lead to detrimental consequences for the patients. On the other hand, if the drug is started 18 or 23 dpi, after the peak of infected cells, the effect of the antiviral drug is to speed up the virus clearing, which is beneficial for the patient. In Figure 6 , the depletion of target cells is also shown for the scenario of the drug started at 18 dpi. As we can see, the depletion is not complete (12% target cells remain), which explains why the decay is slow even with the additional drug activity.

The situation is very different if the antiviral drug efficacy is raised to 90%, namely, ϵ = 0.9 k . When the drug is applied, the effective R becomes 0.7, which leads to virus clearance. We ran the same set of scenarios as above and the results are shown in Figure 7 . The blue line again shows the base scenario with no drug applied. If the drug is applied early, as shown for the scenario with a start date of 5 dpi, the infection can be stopped completely before it makes any real progress. When the drug is started at 10 or 15 dpi, double peaks of viral count are observed. The reason is that when the drug is stopped, R jumps back to 7, so that virus resumes its rapid growth until the target cells are depleted, causing a second peak in virus count. Therefore if the drug is stopped before the virus is cleared out completely, the infection will resurge. This is similar to what we have seen in Figure 4 with the transient antiviral activity of strong AIR, except the peaktrough-peak curve is smoother there. When the drug is applied after the viral peak is reached, it has a small but positive effect of faster recovery, similar to the scenario shown in Figure 6 . The level of target cell depletion is also shown for the scenario of the drug started at 15 dpi. At the first peak of viral load, less than 30% of the target cells are depleted by the infection. This is why a second exponential growth starts after the drug is removed.

To avoid this kind of unwanted effect, the antiviral drug needs to be maintained until the virus is completely cleared. In the Appendix, we showed the relationship between the duration the drug needs to be administered, and the time the drug is started.

The simple rule of thumb is that the later the drug is started, the longer it is required to be applied to ensure complete clearing of the virus.

It is interesting to note that γ-immunoglobulin has been widely used to treat COVID-19 patients in China. In a way, injection of γ-immunoglobulin can be modeled similarly as antiviral drugs F I G U R E 6 Effect of Antiviral drug on viral dynamics: The drug is assumed to be 50% effective on k, and applied to a host with R 0 of 7.03 for 15 days consecutively with various starting dates; viral counts to the left with exponential scale; target cell ∕ T T 0 is shown for the scenario of the drug applied at 18 dpi, to the right with linear scale DU AND YUAN | 9 (assuming that it is able to deliver some antiviral function). It will be interesting to know if some researchers have compiled data to correlate the efficacy of this treatment with the timing of the injections.

Our theory predicts that if the γ-immunoglobulin is used early on, it will most likely delay the depletion of target cells so the duration of the disease will be extended for the patients, hence with poor prognosis. On the other hand, if it is applied after the peak of infected cells, it should be helpful to some extent.

In summary, for an antiviral drug to be useful, the drug needs to be effective enough so that the effective reproductive ratio R falls below 1 (see Equation (5) for the factors affecting R). Under such a condition, it is desirable to apply it early on and must be maintained until the viral load is completely cleared (in the Appendix we provide an estimate for the duration of the regimen). When a drug cannot meet this criterion, applying it before the target cells are fully depleted is typically not beneficial to the patient, and likely will make things worse for the patients. This should have some significant implications for drug companies when they design the clinical trials. At this point, as we do not know much about the potential level of the AIR, R 0 , and the efficacy of the antiviral drug, it is difficult to identify a risk-free way of applying the antiviral drug. A prudent approach is to measure the viral loads frequently post-drug administration. If the viral load does not drop rapidly shortly after, it may indicate that the effective R is not below 1, and the drug should be stopped. If the drug is known to be effective but not effective enough to cause R to go below 1, it might still be useful to apply it after a large percentage of the target cells have been depleted. As a safety measure, it may be desirable to continue the regimen for a few more days after the viral load drops below the detection limit to ensure complete clearance.

Currently, there is already a class of drugs called immunosuppressive drugs (drugs to suppress AIRs). A large body of research has been done on the effects of the drugs on different parts of immune response. 27, 28 The goal of investigating immunosuppressive drugs is to explore their possible applications in COVID-19 disease management.

Immunosuppressants such as corticosteroids have been applied to SARS and COVID-19 patients, with positive results. 29, 30 From the previous discussions, we propose a potential treatment plan by applying immunosuppressive drugs as soon as an infection is detected in the lungs, and remove the drug as soon as the target cells are fully depleted. This approach is simple and only requires the immunosuppressive drugs to be able to reduce the AIR activity to a level low enough not to interfere with the innate immunity. From the discussions around Figure 4 , this approach should clear out the virus quickly if the AIR activity level is below 0.1. If the drug effect is not removed quickly, and the AIR does not rise after the drug is removed, the patients may not be able to completely clear the residue level of the virus. We projected that a short regimen of antiviral drugs 3 to 5 days after the peak of virus counts may help to completely clear the remaining low level of virus load. The purpose of the antiviral drug F I G U R E 7 Effect of antiviral drug on viral dynamics: the drug is assumed to be 90% effective on k, and is applied for 15 days consecutively with various starting dates; viral counts to the left with exponential scale; target cell ∕ T T 0 is for drug started on 15 dpi only, shown with linear scale to the right regimen is to mimic the effect of the AIR in patients with influenza infections.

If the drug effect can be removed immediately after the peak infection, so that the AIR level increases after the drug is removed, we then have a situation that the AIR peaks after the target cells are depleted. The discussions around Figure 5 established that with this scenario, the virus can always be cleared completely, regardless of the level of the AIR. With this approach, the novel coronavirus infections are essentially reduced to influenza-like infections in adults, except with slower disease progression.

One possible concern with our proposed approach is that the immunosuppressive drugs may inevitably have some limited side effects on the innate immune response, in addition to the intended effect on the AIR. This is considered by Equations (4), which gives us the following (see Appendix): To address this possible concern, we did some simulations with the following: the immunosuppressive drug is assumed to have either a 50% or 90% effect on ϵ δ , the drug with 50% effect on ϵ δ is applied from 1 dpi up to 12 dpi, while the second drug is applied from 1 dpi up to 10 dpi. The choice to apply the drug from 1 dpi is merely to simplify our simulation, and may not be necessary, as our main goal is to suppress the AIR. The AIR peak day is assumed to be the same as the base scenario at 8 dpi, and the AIR activity level is reduced to 0.1.

The same base scenario without drug is also shown as the blue line in Figure 8 . With both scenarios, we observed a big reduction of the duration of viral activity. Not only the date of viral peak is moved from 16 to 11 dpi, but the decay from the peak is also much faster.

This rapid decay is beneficial for the patients as well, and it's due to the complete depletion of the target cells. The overall virus replication curve with the immunosuppressive drug resembles closely the curve for adult influenza infections (shown in Figure 2 ) except with a longer time scale. When compared to the base scenario with AIR (the blue line in Figure 8 ), we find that the maximum viral load is increased by a factor of 3.6-4.7. However, if we compare it to the base scenario of innate response alone without AIR, the maximum viral load varies by less than 15%. In summary, we think that the benefit of shortening the disease duration outweighs the cost of increased viral load for a short period.

We see multiple potential advantages of this novel approach.

There are a number of well-studied, FDA-approved drug F I G U R E 8 Effect of immunosuppressive drug on viral dynamics: The drug is assumed to be 50 or 90% effective so ϵ δ = 0.5 or 0.9; drug is applied from 1 dpi until 10 or 12 dpi (with the end day shown on legend); viral counts shown to the left with exponential scale; target cell ∕ T T 0 is shown for 10 dpi for both 50% and 90% drugs, with linear scale to the right DU AND YUAN | 11 candidates to choose from. Animal models can help to define the range for the main parameters quickly. Finally, as long as the drug is applied for enough days so that the suppression of the AIR is beyond the peak of virus load, this approach always resulted in a quick clearance of the virus. Removing the drug in time should enable the AIR to be fully stimulated, so the patient most likely can have antibody-based immunity for the virus in the near future, as in the case of patients recovered from influenza infections. In our discussion so far, we have assumed that the drug effect is immediate, namely, the drug effect shows up immediately after its application and the effect goes away immediately after the drug is removed. This is clearly not realistic, so one of the key considerations when selecting the right immunosuppressive drugs is its pharmacokinetics; a drug with short elimination and absorption half time is more desirable, considering the time-sensitive nature of the effect. In a recent study by Wang et al, 29 the authors reviewed the effect of a low-dose and short-term application of corticosteroid to COVID-19 patients at an early stage, and found some clear benefits of for the patients as compared to the patients without the treatment, including the shortening the disease course. This is consistent with our findings. In their study, the drug corticosteroid was applied for 11 days or more, which according to our theory, might be a little too long for optimal results, considering its half time of 12 to 36 days.

Previous studies have investigated the RNA expression of ACE2 in 72 human tissues. 9 Two more recent studies have shown that the percentage of epithelial cells with ACE2 receptors in bronchial is much lower than in turbinate and lung tissues (0.2% vs 2-3%). 10, 11 Considering the difference in surface area between the lungs and the nasal cavity, [31] [32] [33] we arrived at a rough estimate of 10 4 as the ratio of target cells in these two organs. A quick simulation using the much lower T 0 suggests that the peak of infected cells is changed from 12 to 8 dpi in the nose with the innate immune response alone. If infection is localized at the nasal cavity only, it is then much more likely for the peak of the AIR to come after the peak of infected cells. In this case, most likely the infection will resolve itself just like influenza infections, except the disease may take a few more days. This may explain why a lot of people do not develop severe symptoms or remain asymptomatic if the infection is localized in the nasal area without reaching and propagating in the lungs. If our analysis can be verified clinically, methods and techniques (for example antiviral nasal sprays) to enhance trachea/ bronchial tubes as a barrier to stop the spread of the virus from upper respiratory tracks to the lungs can be effective in reducing the percentage of patients with pneumonia.

Another important parameter we have not discussed is the cell regeneration rate, τ d , of the epithelial cells. In Table 3 , we have listed it as 2.9 × 10 −4 , which corresponds to 10 −3 day −1 for d. This is based on a study by Rawlins and Hogan 34 on the half-life of epithelial cells in the mouse lung. This parameter, if it is higher, will make the depletion of the target cells more difficult. This is because as the target cells are depleted, new target cells are generated. The higher the regeneration rate, the more new target cells are available for viruses to infect and stay alive. Our simulation results remain essentially the same if this parameter is lower or increased by no more than a factor of 20. Beyond that, depletion of the target cells becomes substantially more difficult to achieve. Thereby, clearing the virus also becomes harder to accomplish. Could this be the reason for the observed higher severity rate among older patients? We reason that this is not likely, as the cell regeneration rate should be higher for younger patients, while older people are known to have more difficulty repairing lung damage. 35 In this article, we have applied the target cell-limited model to study the novel corona virus infections mostly in human lungs.

Through comparison between SARS-CoV-2 and influenza viruses, and by analyzing the interactions among various players (target cells, innate and AIRs, and different types of drugs), we propose that the main reason for the higher probability of severe symptoms, prolonged hospital stays, and even fatal outcomes for COVID-19 patients 2,3 may be that the peak of infected cells and virus counts often comes after the peak of the AIR in these patients. This mismatch of timing and the resulted interference with innate immunity by adaptive immunity, leads to incomplete depletion of the target cells, thus providing uninfected target cells for continued infection.

Prolonged infection can induce overactive immune responses, secondary complications, and sometimes fatal outcomes.

We also discussed the prospects of antiviral drugs and immunosuppressive drugs in combating COVID-19 infections. Our proposal to apply immunosuppressive drugs at an early stage to reduce the interference from adaptive immunity so that the innate immunity and the depletion of target cells can together achieve fast elimination of the virus seems to be unconventional. However, it is particularly interesting to note that, as we were working on our mathematical modeling, the recent study by Wang et al 29 demonstrated that an early stage low-dose and short-term application of corticosteroid treatment in patients with severe COVID-19 pneumonia was beneficial and essentially validated our findings. We envision that new insights from our analysis and modeling will encourage more work in this direction.

Clinical features of patients infected with 2019 novel coronavirus in Wuhan

Clinical Characteristics of 138 Hospitalized Patients With 2019 Novel Coronavirus-Infected Pneumonia in Wuhan, China

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study

Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study

A review of mathematical models of influenza A infections within a host or cell culture: lessons learned and challenges ahead

H5N1 virus attachment to lower respiratory tract

Tropism of avian influenza A (H5N1) in the upper and lower respiratory tract

Evolution of the novel coronavirus from the ongoing Wuhan outbreak and modeling of its spike protein for risk of human transmission

Tissue distribution of ACE2 protein, the functional receptor for SARS coronavirus. A first step in understanding SARS pathogenesis

Single-cell RNA expression profiling of ACE2, the putative receptor of Wuhan 2019-nCov

Single-cell RNA expression profiling of ACE2, the putative receptor of Wuhan 2019-nCoV

Principles and Practice of Clinical Virology

Defense mechanisms against influenza virus infection in the respiratory tract mucosa

Incubation period of 2019 novel coronavirus (2019-nCoV) infections among travellers from Wuhan, China

Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia

Mathematical analysis of viral replication dynamics and antiviral treatment strategies: from basic models to age-based multi-scale modeling

Population dynamics of immune responses to persistent viruses

Viral dynamics in hepatitis B virus infection

Virus dynamics and drug therapy

Mathematical models of HIV pathogenesis and treatment

Virus Dynamics: Mathematical Principles of Immunology and Virology

Influenza and the challenge for immunology

Cell-mediated protection in influenza infection

Pathogenesis of emerging Avian influenza viruses in mammals and the host innate immune response

Quantifying the early immune response and adaptive immune response kinetics in mice infected with influenza a virus

Hepatitis C viral dynamics in vivo and the antiviral efficacy of interferon-alpha therapy

Effects of antirejection drugs on innate immune cells after Kidney transplantation

Immunobiology: The Immune System in Health and Disease

Early, low-dose and short-term application of corticosteroid treatment in patients with severe COVID-10 pneumonia: single-center experience from Wuhan, China. medRxiv

Description and clinical treatment of an early outbreak of severe acute respiratory syndrome (SARS) in Guangzhou, PR China

Histology a Text & Atlas

Anatomical and histological factors affecting intranasal drug and vaccine delivery

Ciliated epithelial cell lifespan in the mouse trachea and lung

Regeneration of the aging lung: a mini-review

The authors would like to express their appreciation for Professor Yong-Zhong Qian for some productive discussions. Weiming Yuan is supported by the NIH grant (R21 AI 149365) and partially supported by NIH grant P30CA014089 to the University of Southern California Norris Comprehensive Cancer Center from the National Cancer Center.

The authors declare that there are no conflict of interests.

http://orcid.org/0000-0002-4780-7157

Solving the above equation, we find the following:The time it takes for the infected cells to reach the peak, t r , is as follows, assuming complete depletion of the target cells: Notice that here we used I max and I min , because to clear virus completely we need to get to a level lower than the initial level of infected cells. Immediately we can see that when = R 0, we have the lowest t c , which means the fastest clearing of virus. If we assume that, the above result can be further simplified: