Skip to main content
Advertisement

Main menu

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Methods & Resources
    • Author Interviews
    • Archive
    • Subjects
  • Collections
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Journal of Human Immunity
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research

User menu

  • My alerts

Search

  • Advanced search
Life Science Alliance
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Journal of Human Immunity
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research
  • My alerts
Life Science Alliance

Advanced Search

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Methods & Resources
    • Author Interviews
    • Archive
    • Subjects
  • Collections
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Follow LSA on Bluesky
  • Follow lsa Template on Twitter
Research Article
Source Data
Transparent Process
Open Access

Incomplete antiviral treatment may induce longer durations of viral shedding during SARS-CoV-2 infection

View ORCID ProfileKwang Su Kim, View ORCID ProfileShoya Iwanami, Takafumi Oda, Yasuhisa Fujita, View ORCID ProfileKeiji Kuba, Taiga Miyazaki, View ORCID ProfileKeisuke Ejima  Correspondence email, View ORCID ProfileShingo Iwami  Correspondence email
Kwang Su Kim
1Interdisciplinary Biology Laboratory (iBLab), Division of Biological Science, Graduate School of Science, Nagoya University, Nagoya, Japan
Roles: Data curation, Formal analysis, Investigation, Visualization, Methodology, Writing—original draft, review, and editing
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Kwang Su Kim
Shoya Iwanami
1Interdisciplinary Biology Laboratory (iBLab), Division of Biological Science, Graduate School of Science, Nagoya University, Nagoya, Japan
Roles: Data curation, Formal analysis, Investigation, Visualization, Methodology
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Shoya Iwanami
Takafumi Oda
2Department of Biology, Faculty of Sciences, Kyushu University, Fukuoka, Japan
Roles: Data curation, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing—original draft, review, and editing
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yasuhisa Fujita
1Interdisciplinary Biology Laboratory (iBLab), Division of Biological Science, Graduate School of Science, Nagoya University, Nagoya, Japan
Roles: Data curation, Formal analysis, Methodology
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Keiji Kuba
3Department of Biochemistry and Metabolic Science, Akita University Graduate School of Medicine, Akita, Japan
Roles: Conceptualization, Investigation, Writing—original draft, review, and editing
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Keiji Kuba
Taiga Miyazaki
4Department of Infectious Diseases, Nagasaki University Graduate School of Biomedical Sciences, Nagasaki, Japan
Roles: Conceptualization, Investigation, Writing—original draft, review, and editing
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Keisuke Ejima
5Department of Epidemiology and Biostatistics, Indiana University School of Public Health-Bloomington, Bloomington, IN, USA
Roles: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Writing—original draft, review, and editing
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Keisuke Ejima
  • For correspondence: kejima@iu.edu
Shingo Iwami
1Interdisciplinary Biology Laboratory (iBLab), Division of Biological Science, Graduate School of Science, Nagoya University, Nagoya, Japan
6Institute of Mathematics for Industry, Kyushu University, Fukuoka, Japan
7Institute for the Advanced Study of Human Biology (ASHBi), Kyoto University, Kyoto, Japan
8NEXT-Ganken Program, Japanese Foundation for Cancer Research, Tokyo, Japan
9Interdisciplinary Theoretical and Mathematical Sciences Program (iTHEMS), RIKEN, Saitama, Japan
10Science Groove Inc., Fukuoka, Japan
Roles: Conceptualization, Supervision, Funding acquisition, Investigation, Methodology, Project administration, Writing—original draft, review, and editing
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Shingo Iwami
  • For correspondence: iwami.iblab@bio.nagoya-u.ac.jp
Published 3 August 2021. DOI: 10.26508/lsa.202101049
  • Article
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF
Loading

Abstract

The duration of viral shedding is determined by a balance between de novo infection and removal of infected cells. That is, if infection is completely blocked with antiviral drugs (100% inhibition), the duration of viral shedding is minimal and is determined by the length of virus production. However, some mathematical models predict that if infected individuals are treated with antiviral drugs with efficacy below 100%, viral shedding may last longer than without treatment because further de novo infections are driven by entry of the virus into partially protected, uninfected cells at a slower rate. Using a simple mathematical model, we quantified SARS-CoV-2 infection dynamics in non-human primates and characterized the kinetics of viral shedding. We counterintuitively found that treatments initiated early, such as 0.5 d after virus inoculation, with intermediate to relatively high efficacy (30–70% inhibition of virus replication) yield a prolonged duration of viral shedding (by about 6.0 d) compared with no treatment.

Introduction

Two main processes are involved in the pathogenesis of coronavirus disease 2019 (COVID-19): the disease primarily develops through replication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and is then driven by exaggerated host immune/inflammatory responses to the virus, leading to various tissue damage. Thus, antiviral agents are considered to likely be more beneficial in the earlier stages of COVID-19, whereas immunosuppressive/anti-inflammatory therapies may have the greatest effects in later stages of the disease (1).

Remdesivir (RDV) is currently the only antiviral drug approved by the Food and Drug Administration for the treatment of COVID-19 and is recommended for hospitalized patients who require supplemental oxygen. However, it is not routinely recommended for patients at advanced stages of the disease, such as for patients maintained on artificial respiration by means of mechanical ventilators, owing to the lack of data showing benefit (2, 3, 4, 5). The corticosteroid dexamethasone has been reported to improve survival in hospitalized patients who require supplemental oxygen, especially in patients who require mechanical ventilation (6, 7, 8, 9).

Although mortality is an important and ultimate clinical outcome at both the individual and the population level (10), the efficacy of antiviral drugs that reduce infection with the virus or viral replication may also be evaluated primarily by using viral load kinetics (5, 11, 12, 13). Although viral load was similar among symptomatic and asymptomatic outpatients in a previous study (14), high SARS-CoV-2 viral load on admission is associated with in-hospital mortality in COVID-19 patients who require hospitalization (15, 16). In addition, reporting viral load has an important role in infection prevention practices, given that viral load correlates with levels of infectiousness (17, 18, 19), even though detection of viral RNA by PCR does not necessarily indicate the infectivity of SARS-CoV-2 in the late phase of disease.

The currently available data suggest that reducing the amount of viral load and/or the duration of viral shedding might be beneficial for suppressing early inflammation in patients or for preventing further human-to-human transmission (20). An effective antiviral drug for COVID-19 is urgently needed, and the duration of viral shedding would be an important clinical outcome for assessing drug efficacy. Indeed, to date, several clinical studies have used the duration of viral shedding as an outcome to evaluate the efficacy of candidate drugs (11, 13).

In general, antiviral drugs that have limited efficacy for inhibiting viral replication are not expected to substantially affect virus infection dynamics in the clinical setting regardless of merit. Exhaustive simulations of mathematical models for antiviral treatment; however, sometimes uncover unnatural behaviors even under realistic parameter ranges. For example, it is observed that viral shedding may paradoxically last longer with treatment than without treatment if the antiviral efficacy is below 100% (21, 22). Although even simple limited models of target cells with antiviral treatment blocking virus production (or virus entry) can explain this paradoxical phenomenon, it has not yet been reported in both a clinical and an experimental setting. In a recent article (23), these paradoxical phenomena were observed for the first time in a SARS-CoV-2 infection non-human primate model treated with RDV, which is known as a nucleoside analogue. In the present study, we used a simple mathematical model and its extended models to quantify SARS-CoV-2 infection dynamics in the non-human primate and to characterize viral shedding to better understand the viral infection dynamics behind the paradoxically longer duration of viral shedding with antiviral treatment.

Results

Characterizing SARS-CoV-2 infection dynamics in different specimens from rhesus macaque

We quantitatively analyzed time-course data on viral load from nose and throat swabs of rhesus macaques infected with SARS-CoV-2 (reported in references 23 and 24) using a viral dynamics model (see the Materials and Methods section). To consider interindividual variations in viral loads, we used a nonlinear mixed-effect modelling approach to estimate parameters (see the Materials and Methods section). The estimated fixed effect parameters for nose and throat swabs are listed in Table 1, and the individual estimated parameters for each macaque are listed in Table S1. When we compared the population parameters between specimens, the maximum rate constant for viral replication was statistically significantly larger for nose swabs than for throat swabs (P = 4.12 × 10−10 by the Wald test). The viral load curves based on the individual parameters from the rhesus macaques are depicted for nose (curves in purple) and throat (curves in orange) swabs in Fig 1A.

View this table:
  • View inline
  • View popup
Table 1.

Estimated parameters (fixed effect) for SARS-CoV-2 infection in nose and throat.

Table S1 Individual estimated parameters and initial values in nose and throat swabs.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1. Comparison of viral load trajectory between nose and throat swabs in SARS-CoV-2– infected rhesus macaques without and with treatment.

Viral loads were measured using nasal (purple) and throat (orange) swabs for SARS-CoV-2–infected macaques. (A) The closed dots and curves correspond to the observed and estimated viral loads in the macaques without treatment (A). (B, C) The expected viral load and uninfected target cell proportion in nose, throat, and lung are calculated by the population parameter and shown in (B, C) (without treatment). (D) and (E, F) are same as (A) and (B, C), respectively, but for the group with treatment. (B, C, E, F) The black dashed curves in (E, F) represent the expected viral load and uninfected target cell proportion shown in (B, C). The black horizontal dotted and gray vertical lines show the detection limit of viral load and the timing of initiation of the treatments (0.5 d), respectively.

Source data are available online for this figure.

Source Data for Figure 1[LSA-2021-01049_SdataF1.xlsx]

Expected viral loads in nose and throat swabs (with estimated fixed effect parameters) are shown in Fig 1B and C, respectively. The viral RNA load in nose swabs peaked earlier than in throat swabs: at 0.680 d after virus inoculation compared with 3.09 d after inoculation in throat swabs. The peak timing in the nose was numerically earlier than in the throat, although the difference among specimens was not significant. We further calculated and compared the following quantities, which are composed of the estimated parameters (Table 1): the within-host basic reproduction number (R0 = γ/δ), which is the average number of newly infected cells produced by any single infected cell when target cells are fully uninfected (22), and the Malthusian parameter (M = γ − δ), which is an indicator of the initial growth rate of the viral load (25, 26). Both R0 and M of SARS-CoV-2 were significantly larger in nose swabs than in throat swabs (P = 1.98 × 10−10 and 1.54 × 10−11 by the Wald test, respectively) (Table 1). Thus, both the larger R0 and M and the earlier peak viral load in nose swabs suggest that the virus replicates and spreads more effectively in the nose than in the throat, implying that the infection dynamics of SARS-CoV-2 differs in different organs. These organ-specific infection dynamics may thus affect treatment outcomes in different organs. In the next section, we provide a detailed analysis of anti-SARS-CoV-2 treatment with RDV in infected rhesus macaques.

Evaluation of antiviral treatment with remdesivir in SARS-CoV-2–infected rhesus macaque

In addition to our finding of the organ specificity of virus infection dynamics, other differences between organs may affect drug efficacy, such as target cell constitutions, drug distributions in plasma and tissue compartments, etc. (27, 28). To quantitatively evaluate how these organ-specific heterogeneities affect treatment outcomes, we further analyzed time-course data on viral load from the same two specimens of infected rhesus macaques treated with RDV. Note that the study of RDV-treated rhesus macaque reported that there was no significant reduction in viral load in the nose or throat at each time point although RDV suppressed virus infection in the lower respiratory tract (23). Here we used the above parameter values obtained from the SARS-CoV-2–infected rhesus macaques without treatment, and then estimated the antiviral efficacy of RDV against SARS-CoV-2 in the nose and throat (Table 1). The viral load curves are depicted with observed data in Fig 1D for nose and throat swabs.

The fixed effect of the antiviral efficacy in throat is estimated to be ε =0.618 (95% CI: 0.373–0.863), although we failed to estimate individual antiviral efficacy in nose (see below). This is because the treatment with RDV was initiated at day 0.5, which was close to the peak viral load in nose (i.e., day 0.680; Fig 1B). Because only a very small fraction of target cells remain uninfected after the viral load peak (29, 30), the number of ongoing de novo infections that the RDV can potentially interrupt is limited. Therefore, even if RDV blocked virus production with relatively high efficacy up to 100%, the viral load decay was not significantly influenced (Figs 1E and S1A). We also quantitatively analyzed time-course data on viral load from bronchoalveolar lavage (BAL) samples (Fig S2A–D). The viral load decay was not influenced by treatment (Fig S2D) owing to the late initiation of RDV therapy (i.e., day 0.5), which was after the viral load peak in BAL fluid (i.e., 0.360 d; Fig S1B). In fact, there was little difference in the duration of viral shedding between rhesus macaques without and with RDV treatment (Figs 1E and S2D). On the other hand, because the initiation of RDV (i.e., day 0.5) was before the viral load peak in throat (i.e., day 3.09; Fig 1C), blocking virus production affected the viral load and the duration of viral shedding (Fig 1F). We further investigated the treatment effect of varying the treatment efficacy and the timing of treatment initiation.

Figure S1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S1. Different treatment outcomes in nose when drug efficacy is varied.

(A, B) The expected viral load and uninfected target cell proportion in (A) nose and (B) lung are calculated under the condition of antiviral treatment blocking virus production with different efficacies (ε = 0.00, 0.30, 0.70, 0.90, and 1.00 corresponding to the blue, light blue, green, yellow, and red solid lines, respectively) at 0.5 d after inoculation. The black horizontal and gray vertical lines show the detection limit of the viral load and the timing of initiation of the treatments, respectively.

Source data are available online for this figure.

Source Data for Figure S1[LSA-2021-01049_SdataFS1.xlsx]

Figure S2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S2. Viral load trajectory in lung of SARS-CoV-2–infected rhesus macaques without and with treatment.

Viral loads were measured using BAL for SARS-CoV-2–infected macaques. (A) The closed dots and curves correspond to the observed and estimated viral loads in the macaques without treatment (A). (B) The expected viral load and uninfected target cell proportion in lung are calculated by the population parameters and shown in (B) (without treatment). (C) and (D) are same as (A) and (B), respectively, but for the group with treatment. (B, D) The black dashed curves in (D) represent the expected viral load and uninfected target cell proportion shown in (B). The black horizontal dotted and gray vertical lines show the detection limit of the viral load and the timing of initiation of treatments (0.5 d), respectively. Comparing population parameters between specimens, the maximum rate constant for viral replication in lung was statistically significantly larger than that in nose (P = 1.97 × 10−3 by the Wald test). On the other hand, the rate constant for virus infection was statistically significantly smaller than that in nose (P = 3.57 × 10−2 by the Wald test). Also, both R0 and M of SARS-CoV-2 in lung were statistically significantly larger than those in nose (P = 7.35 × 10−4 and 2.55 × 10−4 by the Wald test, respectively). The estimated parameters (fixed effect and individual parameter) and initial values are summarized in Tables S2 and S3.

Source data are available online for this figure.

Source Data for Figure S2[LSA-2021-01049_SdataFS2.xlsx]

Table S2 Estimated parameters (fixed effect) for SARS-CoV-2 infection in the lung.

Table S3 Individual estimated parameters and initial values in the lung.

Paradoxically longer durations of viral shedding with antiviral treatment

In our calculations using the nose data (see Fig 1E), the duration of viral shedding was not influenced much by treatment because ongoing de novo infection was limited when treatment was initiated. By contrast, in the throat, the duration of viral shedding was determined by a balance between de novo infection and removal of infected cells. Surprisingly, as shown in Fig 1F, we found a trend for viral shedding to last longer in the macaques treated with RDV than in those without treatment despite more than 60% of de novo infection being blocked.

To elucidate a viral infection dynamics behind this counterintuitive longer duration of viral shedding in the throat, we further analyzed the target cell limited model for throat by varying drug efficacy and the timing of treatment initiation (Fig 2). If the efficacy of treatment for blocking virus production was close to 100% (Fig 2A and B), the duration of viral shedding was shortened because all ongoing de novo infections were blocked. However, with an intermediate efficacy of treatment (10–70%), the duration became longer because uninfected cells were temporally protected by the treatment but some ongoing de novo infections were not blocked. Thus, the protected cells were gradually infected and drove further but slower infection than in the case without treatment (Fig 2A and B). For low-efficacy treatments (<10%), the treatments no longer protected uninfected cells and de novo infections occurred at almost same level without treatment (i.e., the duration became shorter again) (Fig 2A and B). This trend was highlighted when the treatment was initiated early. As we discussed above, limited uninfected target cells remained after the viral load hit its peak. Thus, the treatment impact on the duration of viral shedding was limited when the treatment was initiated after the viral load peak (Fig 2C, see also Fig 1E). We further compared the model we used in this study and other extended models, all of which have been used to describe virus dynamics of SARS-CoV-2 and other viruses, to show our conclusion is consistent without loss of generality (see Figs S3A–H, S4A–H, and S5A–H).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2. Different treatment outcomes in throat as the result of varying drug efficacy and the timing of initiation.

(A, B, C) The expected viral load and uninfected target cell proportion are calculated under the condition of antiviral treatment blocking virus production with different efficacies (ε = 0.00, 0.10, 0.30, 0.70, 0.90, and 1.00, corresponding to the blue, light blue, green, yellow, orange, and red solid lines, respectively) and with initiation at 0.5, 2.0, or 5.0 d after inoculation using the estimated model parameter values for throat shown in (A, B, C), respectively. The black horizontal dotted and gray vertical lines show the detection limit of viral load and the timing of initiation of the treatments, respectively.

Source data are available online for this figure.

Source Data for Figure 2[LSA-2021-01049_SdataF2.xlsx]

Figure S3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S3. Comparison of viral load trajectory derived from estimated parameters of the model with eclipse phase between nose and throat swabs in rhesus macaques infected with SARS-CoV-2 without and with treatment.

We used a virus dynamics model considering an eclipse phase of infection used for SARS-CoV-2 and other viruses (51). The following is derived from the eclipse phase model and applied to data fitting here: dfT(t)dt=−βfT(t)V(t),dfE(t)dt=βfT(t)V(t)−kfE(t),dV(t)/dt=(1−εH(t))γfE(t)−δV(t) where fT(t) is the relative fraction of the uninfected target cell population at day t to that at day 0 (i.e., fT(0) = 1), fE(t) is that of eclipse phase cells, and V(t) is the amount of virus at day t, respectively. Note that time 0 corresponds to the day of inoculation. The parameters γ, β, and δ are the maximum rate constant for viral replication, the rate constant for virus infection, and the death rate of infected cells, respectively. 1/k is the average duration of the eclipse phase. Here, we fixed k as 3 d−1 (52) and assumed fE(0) as 0. Random effect and categorical covariate are assumed for γ, β, δ, ε, V(0) and γ, which provides the lowest Bayesian information criterion, respectively. (A) The closed dots and curves correspond to the observed and estimated viral loads in the macaques without treatment (A). (B, C, D) The expected viral load and uninfected target cell proportion in nose, throat, and lung are calculated by the population parameter and shown in (B, C, D) (without treatment). (E) and (F, G, H) are same as (A) and (B, C, D), respectively, but for the group with treatment. (B, C, D, F, G, H) The black dashed curves in (F, G, H) represent the expected viral load and uninfected target cell proportion shown in (B, C, D). The black horizontal dotted and gray vertical lines show the detection limit of viral load and the timing of initiation of the treatments (0.5 d), respectively. The estimated parameter values are listed in Table S4.

Source data are available online for this figure.

Source Data for Figure S3[LSA-2021-01049_SdataFS3.xlsx]

Table S4 Estimated parameters (fixed effect) for SARS-CoV-2 infection in nose and throat by fitting the viral dynamics with eclipse phase.

Figure S4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S4. Comparison of viral load trajectory derived from estimated parameters of the model with multiple target cells between nose and throat swabs in rhesus macaques infected with SARS-CoV-2 without and with treatment.

We used a virus dynamics model considering multiple target cells and applied to data fitting here (53): df1(t)dt=−β1f1(t)V(t),df2(t)dt=−β2f2(t)V(t), dV(t)/dt=(1−εH(t))γ1f1(t)V(t)+(1−εH(t))γ2f2(t)V(t)−δV(t), where f1(t) and f2(t) are the relative fraction of the first and second uninfected target cell population at day t (i.e., f1(0) + f2(0) = 1), and V(t) is the amount of virus at day t, respectively. Note that time 0 corresponds to the day of inoculation. The parameters γ1 and γ2 are the maximum rate constants for viral replication for first and second target cells, respectively. β1 and β2 represent the rate constants for virus infection for first and second target cells, respectively. δ is the death rate of infected cells. Random effect and categorical covariate are used to γ1, γ2, β1, β2, δ, ε, V(0), f1(0), f2(0) and γ1, β1, which provide the lowest Bayesian information criterion, respectively. (A) The closed dots and curves correspond to the observed and estimated viral loads in the macaques without treatment (A). (B, C, D) The expected viral load and uninfected target cell proportion (f1(t) + f2(t)) in nose, throat, and lung are calculated by the population parameter and shown in (B, C, D) (without treatment). (E) and (F, G, H) are same as (A) and (B, C, D), respectively, but for the group with treatment. (B, C, D, F, G, H) The black dashed curves in (F, G, H) represent the expected viral load and uninfected target cell proportion shown in (B, C, D). The black horizontal dotted and gray vertical lines show the detection limit of viral load and the timing of initiation of the treatments (0.5 d), respectively. The estimated parameter values are listed in Table S5.

Source data are available online for this figure.

Source Data for Figure S4[LSA-2021-01049_SdataFS4.xlsx]

Table S5 Estimated parameters (fixed effect) for SARS-CoV-2 infection in nose and throat by fitting the viral dynamics with multiple target cell.

Figure S5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S5. Comparison of viral load trajectory derived from estimated parameters of the model with innate immune response between nose and throat swabs in rhesus macaques infected with SARS-CoV-2 without and with treatment.

We used a virus dynamics model considering innate immune response by IFNs and applied to data fitting here (51): df(t)dt=−βf(t)V(t)1+ηF(t),dV(t)dt=(1−εH(t))γf(t)V(t)1+ηF(t)−δV(t),dF(t)dt=sV(t)−αF(t), where f(t) = T(t)/T(0) is the relative fraction of the uninfected target cell population at day t, V(t) is the amount of virus at day t, and F(t) is the concentration of IFNs produced from infected cells, respectively. Note that time 0 corresponds to the day of inoculation. The parameters γ and β are the maximum rate constants for viral replication and the rate constant for virus infection, respectively. δ is the death rate of infected cells. IFN is assumed to be secreted proportional to the amount of virus at a rate of s = 1.0 × 10−5 and cleared at a rate of α = 0.4. We also fixed F(0) = 1, η = 0.1, and V(0) = 1 without loss of generality. Random effect and categorical covariate are used to γ, β, δ and γ, β, respectively. (A) The closed dots and curves correspond to the observed and estimated viral loads in the macaques without treatment (A). (B, C, D) The expected viral load and uninfected target cell proportion in nose, throat, and lung are calculated by the population parameter and shown in (B, C, D) (without treatment). (E) and (F, G, H) are same as (A) and (B, C, D), respectively, but for the group with treatment. (B, C, D, F, G, H) The black dashed curves in (F, G, H) represent the expected viral load and uninfected target cell proportion shown in (B, C, D). The black horizontal dotted and gray vertical lines show the detection limit of viral load and the timing of initiation of the treatments (0.5 d), respectively. The estimated parameter values are listed in Table S6.

Source data are available online for this figure.

Source Data for Figure S5[LSA-2021-01049_SdataFS5.xlsx]

Table S6 Estimated parameters (fixed effect) for SARS-CoV-2 infection in nose and throat by fitting the viral dynamics with innate immune response.

Taken together, under incomplete antiviral treatments, our findings show that uninfected target cells remain longer and infection slowly continues depending on target cell availability. The availability of target cells might explain the paradoxical phenomenon in different organs first observed in the study of SARS-CoV-2 infection in the non-human primate model.

Discussion

It is well known that entry of SARS-CoV-2 into cells requires interaction of the virus spike protein with host angiotensin-converting enzyme 2 (ACE2) receptor, and cleavage and activation of spike by a serine protease, TMPRSS2 (31). Analyses of single-cell sequencing datasets have demonstrated that ACE2 receptor is abundantly expressed in multiple organs, such as nose, lung, eye, and intestine (27), which suggests that the target organ of SARS-CoV-2 is not limited to respiratory tracts. Given that the expression level of ACE2 varies in different organs, virus infection dynamics is expected to differ as well. In fact, although the differences might be due to the sensitivity of the quantification methods used, different viral loads of SARS-CoV-2 in different specimens have been reported (32, 33, 34, 35). Whereas interferon signalling upon SARS-CoV-2 infection is suggested to induce ACE2 expression in single-cell levels (27), cellular entry of SARS coronavirus leads to down-regulation of bulk amounts of ACE2 protein in the lungs and heart (36, 37), which suggests a complicated mechanism for virus spread in the tissue microenvironment during the course of de novo infection of SARS-CoV-2. Thus, different dynamics of the virus and regulation of host factors in the organs might contribute to a paradoxically prolonged duration of virus shedding with incomplete RDV treatment.

In this study, we quantitatively analyzed time-course data on viral load from different specimens, such as nasal swabs, throat swabs, and BAL fluid, of SARS-CoV-2–infected rhesus macaques without and with RDV treatment (23, 24). Because direct and multi-route inoculations with a high viral titer were used in the infection experiments in the rhesus macaques, we observed extremely rapid infections in multiple organs compared with natural infections in humans via respiratory transmission. However, in these infection experiments, because the multi-route inoculations led to effective and rapid infections across organs, we found that the peak viral load appeared faster and that R0 was larger in the nose than in the throat (Fig 1), implying different treatment outcomes in different organs. Because the initiation of the RDV treatment in the nose and lung (i.e., in BAL fluid) was close to the peak of viral load, the RDV had little effect on altering the viral load decay even though it completely blocked virus production; therefore, we could not observe an antiviral effect on the duration of virus shedding (Figs 1E and S2D). In contrast, interestingly, we demonstrated that in the throat an efficacy of RDV of more than 60% induced longer durations of viral shedding (Fig 1F). This paradoxically longer duration could be explained by our simple target cell limited mathematical model: because an incomplete antiviral treatment reduces ongoing infections and protects uninfected target cells, progeny virus that escape from treatment further induce de novo infections into protected, uninfected cells but at a slower rate (Fig 2). We quantitatively revealed the effect of RDV on SARS-CoV-2 replication in throat in treated rhesus macaque and explained why no clear effect of treatment in nose and lung is observed by using a mathematical model.

In many clinical trials, the duration of viral shedding has been used as a primary outcome to evaluate the antiviral effect of treatments with lopinavir and ritonavir (11), hydroxychloroquine (38, 39), and meplazumab (40). Because the viral load of SARS-CoV-2 peaks on or before symptom onset in many patients (18), and the mean interval between symptom onset and hospitalization is estimated to be 4.64 d (41), it is expected that antiviral treatments initiated after symptom onset will not shorten the duration of viral shedding. However, prophylactic use of antiviral drugs or contact-tracing-based (i.e., identifying patients before symptom onset) treatments may shorten the duration of viral shedding.

There are various studies using mathematical models that can explain SARS-CoV-2 infection dynamics and give insight into treatment strategy (29, 42, 43, 44). Another mathematical study of nasal viral load also showed a longer duration of virus shedding under RDV treatment; however, that study did not directly estimate treatment efficacy (45). Moreover, the previous discussion of the difference in SARS-CoV-2 viral dynamics in lung and nose used a more detailed mathematical model (46 Preprint). Because our mathematical model (two-dimensional ordinary differential equations, or ODEs) is derived from a basic virus dynamics model (three-dimensional ODEs [21]) by assuming quasi-steady state, we can reduce the number of parameters needing to be estimated. Note that there might be a discrepancy in predicted virus infection dynamics between the two-dimensional and the three-dimensional ODEs if the antiviral efficacy is close to 100%. However, because the estimated antiviral efficacy is around 60%, the two-dimensional ODE well approximates the virus infection dynamics for the macaques treated with RDV. Ours is the first study to reveal the treatment conditions that induce longer durations of viral shedding based on data collected during treatment.

The limited number of rhesus macaques in the SARS-CoV-2 infection experiments (28 and 17 animals with and without treatment, respectively, in references 23 and 24) used in our analysis might be a limitation. For example, we could not obtain statistically significant differences for the time from virus inoculation to the viral load peak between nose and throat swabs because of large variance due to the small sample size. Once more experimental datasets become available, our mathematical model may be able to more precisely differentiate the virus infection dynamics of different organs. In addition, although relevant data and evidence for prophylactic use of antiviral drugs for SARS-CoV-2 infection in humans are not available so far, future studies should verify our findings by using human data. Thus, collecting longitudinal data and quantifying the viral dynamics with mathematical models are important to evaluate outcomes of antiviral treatments.

Materials and Methods

Study data

The longitudinal viral load was extracted from the published studies of SARS-CoV-2 infection in a non-human primate model (23). The rhesus macaques were inoculated with a total dose of 2.6 × 106 TCID50 through intranasal, oral, ocular, and intratracheal routes. After the inoculations, the viral load was measured every day. Also, BAL was performed using 10 ml of sterile saline to extract viral RNA in BAL fluid. The rhesus macaques were assigned to receive either an intravenous dose of RDV (10 mg/kg loading dose on 0.5 d after inoculation with SARS-CoV-2, followed by 5 mg/kg daily) or an equal volume of vehicle solution (2 ml/kg loading dose, followed by 1 ml/kg daily). We also added viral loads from nose and throat of rhesus macaques infected with SARS-CoV-2 with the same experimental protocol in (24). Data from macaques with less than three data points above the detection limit were excluded. Finally, we used data from 28 and 17 animals without and with treatment, respectively.

Mathematical model

The following mathematical model describes SARS-CoV-2 viral dynamics (29, 30, 47):df(t)dt=−βf(t)V(t),(1)dV(t)dt=(1−εH(t))γf(t)V(t)−δV(t).(2)

The variables f(t) and V(t) are the relative fraction of uninfected target cells at time t to those at time 0, and the amount of virus at time, respectively. The parameters β, γ, and δ represent the rate constant for virus infection, the maximum rate constant for viral replication, and the death rate of virus-producing cells, respectively. To describe the off- and on-treatments, we employed a Heaviside function: H(t) = 0 if t <t* (i.e., before treatment initiation); otherwise H(t) = 1. ε is the efficacy of antiviral treatment blocking virus production (0 < ε ≤ 1). If ε = 1, the virus replication from the infected cells is perfectly inhibited (i.e., the antiviral efficacy is 100%). We evaluated the expected antiviral effect of the treatment on the duration of virus shedding under different inhibition rates (ε) and initiation timings (t*) [Equations (1) and (2)].

The nonlinear mixed-effect model

All viral load data were fitted using a nonlinear mixed-effect modelling approach, which estimates population parameters while accounting for interindividual variation in virus dynamics. The model included both a fixed effect (constant across rhesus macaques) and a random effect (different between rhesus macaques) for each parameter. The parameter values for rhesus macaque i can be expressed as ϑi(=ϑ×eπi), where ϑ is a fixed effect and eπi is a random effect; πi is normally distributed as N(0, Ω). We estimated fixed effects and random effects using the stochastic approximation expectation-maximization algorithm and the empirical Bayes’s method, respectively. In addition, we used location of specimens for measuring viral load as a categorical covariate in estimating γ and β which provide the lowest Bayesian information criterion. Fitting was carried out using MONOLIX 2019R2. To handle data under the detection limit, the likelihood was constructed assuming that the data are within the interval from 0 to α, where α is the lowest observed viral load (above the detection limit) around the day when the viral load is censored (48). The viral load trajectories were depicted using the best fit parameter estimates for individual data in Figs 1 and S1. The estimated parameters (fixed effect and individual parameter) and initial values are summarized in Tables 1 and S1.

Statistical test

To evaluate statistical differences for each of the estimated parameters (γ, β, R0, and M), we applied the Wald test. The time from virus inoculation to viral load peak was calculated at the population level by running the model using estimated fixed effects and the initial values. To compare the timing of viral load peak in the nose, throat, and lung, the differences of these values were tested by the Jackknife test (49, 50).

Acknowledgements

This study was supported in part by the Basic Science Research Program through the National Research Foundation of Korea funded by the Ministry of Education 2019R1A6A3A12031316 (to KS Kim); Grants-in-Aid for JSPS Research Fellow 19J12319 (to S Iwanami); Scientific Research (KAKENHI) B 18KT0018 (to S Iwami), 18H01139 (to S Iwami), 16H04845 (to S Iwami), Scientific Research in Innovative Areas 20H05042 (to S Iwami), 19H04839 (to S Iwami), 18H05103 (to S Iwami); AMED CREST 19gm1310002 (to S Iwami); AMED Japan Program for Infectious Diseases Research and Infrastructure, 20wm0325007h0001, 20wm0325004s0201, 20wm0325012s0301, 20wm0325015s0301 (to S Iwami); AMED Research Program on HIV/AIDS 19fk0410023s0101 (to S Iwami); AMED Research Program on Emerging and Re-emerging Infectious Diseases 19fk0108050h0003 (to S Iwami), 19fk0108156h0001 (to S Iwami), 20fk0108140s0801 (to S Iwami) and 20fk0108413s0301 (to S Iwami); AMED Program for Basic and Clinical Research on Hepatitis 19fk0210036h0502 (to S Iwami); AMED Program on the Innovative Development and the Application of New Drugs for Hepatitis B 19fk0310114h0103 (to S Iwami); JST MIRAI (to S Iwami); Moonshot R&D grant Number JPMJMS2021 (to S Iwami) and JPMJMS2025 (to S Iwami); Mitsui Life Social Welfare Foundation (to S Iwami); Shin-Nihon of Advanced Medical Research (to S Iwami); Suzuken Memorial Foundation (to S Iwami); Life Science Foundation of Japan (to S Iwami); SECOM Science and Technology Foundation (to S Iwami); The Japan Prize Foundation (to S Iwami); Daiwa Securities Health Foundation (to S Iwami).

Author Contributions

  • KS Kim: data curation, formal analysis, investigation, visualization, methodology, and writing—original draft, review, and editing.

  • S Iwanami: data curation, formal analysis, investigation, visualization, and methodology.

  • T Oda: data curation, formal analysis, validation, investigation, visualization, methodology, and writing—original draft, review, and editing.

  • Y Fujita: data curation, formal analysis, and methodology.

  • K Kuba: conceptualization, investigation, and writing—original draft, review, and editing.

  • T Miyazaki: conceptualization, investigation, and writing—original draft, review, and editing.

  • K Ejima: conceptualization, data curation, formal analysis, investigation, methodology, project administration, and writing—original draft, review, and editing.

  • S Iwami: conceptualization, supervision, funding acquisition, investigation, methodology, project administration, and writing—original draft, review, and editing.

Conflict of Interest Statement

The authors declare that they have no conflict of interest.

  • Received February 5, 2021.
  • Revision received July 23, 2021.
  • Accepted July 23, 2021.
  • © 2021 Kim et al.
Creative Commons logoCreative Commons logohttps://creativecommons.org/licenses/by/4.0/

This article is available under a Creative Commons License (Attribution 4.0 International, as described at https://creativecommons.org/licenses/by/4.0/).

References

  1. 1.↵
    1. Wiersinga WJ,
    2. Rhodes A,
    3. Cheng AC,
    4. Peacock SJ,
    5. Prescott HC
    (2020) Pathophysiology, transmission, diagnosis, and treatment of coronavirus disease 2019 (COVID-19): A review. JAMA 324: 782–793. doi:10.1001/jama.2020.12839
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Wang Y,
    2. Zhang D,
    3. Du G,
    4. Du R,
    5. Zhao J,
    6. Jin Y,
    7. Fu S,
    8. Gao L,
    9. Cheng Z,
    10. Lu Q, et al.
    (2020) Remdesivir in adults with severe COVID-19: A randomised, double-blind, placebo-controlled, multicentre trial. Lancet 395: 1569–1578. doi:10.1016/S0140-6736(20)31022-9
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Spinner CD,
    2. Gottlieb RL,
    3. Criner GJ,
    4. Arribas López JR,
    5. Cattelan AM,
    6. Soriano Viladomiu A,
    7. Ogbuagu O,
    8. Malhotra P,
    9. Mullane KM,
    10. Castagna A, et al.
    (2020) Effect of remdesivir vs standard care on clinical status at 11 days in patients with moderate COVID-19: A randomized clinical trial. JAMA 324: 1048–1057. doi:10.1001/jama.2020.16349
    OpenUrlCrossRefPubMed
  4. 4.↵
    1. Goldman JD,
    2. Lye DCB,
    3. Hui DS,
    4. Marks KM,
    5. Bruno R,
    6. Montejano R,
    7. Spinner CD,
    8. Galli M,
    9. Ahn MY,
    10. Nahass RG, et al.
    (2020) Remdesivir for 5 or 10 days in patients with severe covid-19. N Engl J Med 383: 1827–1837. doi:10.1056/NEJMoa2015301
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. Beigel JH,
    2. Tomashek KM,
    3. Dodd LE,
    4. Mehta AK,
    5. Zingman BS,
    6. Kalil AC,
    7. Hohmann E,
    8. Chu HY,
    9. Luetkemeyer A,
    10. Kline S, et al.
    (2020) Remdesivir for the treatment of covid-19: Final report. N Engl J Med 383: 1813–1826. doi:10.1056/NEJMoa2007764
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Horby P,
    2. Lim WS,
    3. Emberson JR,
    4. Mafham M,
    5. Bell JL,
    6. Linsell L,
    7. Staplin N,
    8. Brightling C,
    9. Ustianowski A,
    10. Elmahi E, et al., RECOVERY Collaborative Group
    , (2020) Dexamethasone in hospitalized patients with covid-19. N Engl J Med 384: 693–704. doi:10.1056/NEJMoa2021436
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Jeronimo CMP,
    2. Brito-Sousa JD,
    3. Borba MGS,
    4. Baía-da-Silva DC,
    5. Val FFA,
    6. Monteiro WM,
    7. Melo GC,
    8. Sampaio VS,
    9. Alexandre MAA,
    10. Xavier M, et al.
    (2020) Methylprednisolone as adjunctive therapy for patients hospitalized with Coronavirus Disease 2019 (COVID-19; metcovid): A randomized, double-blind, phase IIb, placebo-controlled trial. Clin Infect Dis 72: e373–e381. doi:10.1093/cid/ciaa1177
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Tomazini BM,
    2. Maia IS,
    3. Cavalcanti AB,
    4. Berwanger O,
    5. Rosa RG,
    6. Veiga VC,
    7. Avezum A,
    8. Lopes RD,
    9. Bueno FR,
    10. Silva MVAO, et al.
    (2020) Effect of dexamethasone on days alive and ventilator-free in patients with moderate or severe acute respiratory distress syndrome and COVID-19: The CoDEX randomized clinical trial. JAMA 324: 1307–1316. doi:10.1001/jama.2020.17021
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Angus DC,
    2. Derde L,
    3. Al-Beidh F,
    4. Annane D,
    5. Arabi Y,
    6. Beane A,
    7. van Bentum-Puijk W,
    8. Berry L,
    9. Bhimani Z,
    10. Bonten M, et al.
    (2020) Effect of hydrocortisone on mortality and organ support in patients with severe COVID-19: The REMAP-CAP COVID-19 corticosteroid domain randomized clinical trial. JAMA 324: 1317–1329. doi:10.1001/jama.2020.17022
    OpenUrlCrossRefPubMed
  10. 10.↵
    1. Pan H,
    2. Peto R,
    3. Henao-Restrepo AM,
    4. Preziosi MP,
    5. Sathiyamoorthy V,
    6. Abdool Karim Q,
    7. Alejandria MM,
    8. Hernández García C,
    9. Kieny MP,
    10. Malekzadeh R, et al., WHO Solidarity Trial Consortium
    , (2020) Repurposed antiviral drugs for covid-19: Interim WHO solidarity trial results. N Engl J Med 384: 497–511. doi:10.1056/NEJMoa2023184
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Li Y,
    2. Xie Z,
    3. Lin W,
    4. Cai W,
    5. Wen C,
    6. Guan Y,
    7. Mo X,
    8. Wang J,
    9. Wang Y,
    10. Peng P, et al.
    (2020) Efficacy and safety of lopinavir/ritonavir or arbidol in adult patients with mild/moderate COVID-19: An exploratory randomized controlled trial. Med (N Y) 1: 105–113.e4. doi:10.1016/j.medj.2020.04.001
    OpenUrlCrossRef
  12. 12.↵
    1. RECOVERY Collaborative Group
    (2020) Lopinavir-ritonavir in patients admitted to hospital with COVID-19 (RECOVERY): A randomised, controlled, open-label, platform trial. Lancet 396: 1345–1352. doi:10.1016/S0140-6736(20)32013-4
    OpenUrlCrossRefPubMed
  13. 13.↵
    1. Doi Y,
    2. Hibino M,
    3. Hase R,
    4. Yamamoto M,
    5. Kasamatsu Y,
    6. Hirose M,
    7. Mutoh Y,
    8. Homma Y,
    9. Terada M,
    10. Ogawa T, et al.
    (2020) A prospective, randomized, open-label trial of early versus late favipiravir therapy in hospitalized patients with COVID-19. Antimicrob Agents Chemother 64: e01897-20. doi:10.1128/AAC.01897-20
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    1. Lee S,
    2. Kim T,
    3. Lee E,
    4. Lee C,
    5. Kim H,
    6. Rhee H,
    7. Park SY,
    8. Son HJ,
    9. Yu S,
    10. Park JW, et al.
    (2020) Clinical course and molecular viral shedding among asymptomatic and symptomatic patients with SARS-CoV-2 infection in a community treatment center in the Republic of Korea. JAMA Intern Med 180: 1447–1452. doi:10.1001/jamainternmed.2020.3862
    OpenUrlCrossRef
  15. 15.↵
    1. Westblade LF,
    2. Brar G,
    3. Pinheiro LC,
    4. Paidoussis D,
    5. Rajan M,
    6. Martin P,
    7. Goyal P,
    8. Sepulveda JL,
    9. Zhang L,
    10. George G, et al.
    (2020) SARS-CoV-2 viral load predicts mortality in patients with and without cancer who are hospitalized with COVID-19. Cancer Cell 38: 661–671.e2. doi:10.1016/j.ccell.2020.09.007
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Magleby R,
    2. Westblade LF,
    3. Trzebucki A,
    4. Simon MS,
    5. Rajan M,
    6. Park J,
    7. Goyal P,
    8. Safford MM,
    9. Satlin MJ, et al.
    (2020) Impact of SARS-CoV-2 viral load on risk of intubation and mortality among hospitalized patients with coronavirus disease 2019. Clin Infect Dis 10: ciaa851. doi:10.1093/cid/ciaa851
    OpenUrlCrossRef
  17. 17.↵
    1. Bullard J,
    2. Dust K,
    3. Funk D,
    4. Strong JE,
    5. Alexander D,
    6. Garnett L,
    7. Boodman C,
    8. Bello A,
    9. Hedley A,
    10. Schiffman Z, et al.
    (2020) Predicting infectious severe acute respiratory syndrome coronavirus 2 from diagnostic samples. Clin Infect Dis 71: 2663–2666. doi:10.1093/cid/ciaa638
    OpenUrlCrossRefPubMed
  18. 18.↵
    1. He X,
    2. Lau EHY,
    3. Wu P,
    4. Deng X,
    5. Wang J,
    6. Hao X,
    7. Lau YC,
    8. Wong JY,
    9. Guan Y,
    10. Tan X, et al.
    (2020) Temporal dynamics in viral shedding and transmissibility of COVID-19. Nat Med 26: 672–675. doi:10.1038/s41591-020-0869-5
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Marks M,
    2. Millat-Martinez P,
    3. Ouchi D,
    4. Roberts CH,
    5. Alemany A,
    6. Corbacho-Monné M,
    7. Ubals M,
    8. Tobias A,
    9. Tebé C,
    10. Ballana E, et al.
    (2021) Transmission of COVID-19 in 282 clusters in catalonia, Spain: A cohort study. Lancet Infect Dis 21: 629–636. doi:10.1016/S1473-3099(20)30985-3
    OpenUrlCrossRefPubMed
  20. 20.↵
    1. Mitjà O,
    2. Clotet B
    (2020) Use of antiviral drugs to reduce COVID-19 transmission. Lancet Glob Health 8: e639–e640. doi:10.1016/S2214-109X(20)30114-5
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. Nowak M,
    2. May RM
    (2000) Virus Dynamics: Mathematical Principles of Immunology and Virology: Mathematical Principles of Immunology and Virology. UK: Oxford University Press.
  22. 22.↵
    1. Perelson AS
    (2002) Modelling viral and immune system dynamics. Nat Rev Immunol 2: 28–36. doi:10.1038/nri700
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Williamson BN,
    2. Feldmann F,
    3. Schwarz B,
    4. Meade-White K,
    5. Porter DP,
    6. Schulz J,
    7. van Doremalen N,
    8. Leighton I,
    9. Yinda CK,
    10. Pérez-Pérez L, et al.
    (2020) Clinical benefit of remdesivir in rhesus macaques infected with SARS-CoV-2. Nature 585: 273–276. doi:10.1101/2020.04.15.043166
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. Munster VJ,
    2. Feldmann F,
    3. Williamson BN,
    4. van Doremalen N,
    5. Pérez-Pérez L,
    6. Schulz J,
    7. Meade-White K,
    8. Okumura A,
    9. Callison J,
    10. Brumbaugh B, et al.
    (2020) Respiratory disease in rhesus macaques inoculated with SARS-CoV-2. Nature 585: 268–272. doi:10.1038/s41586-020-2324-7
    OpenUrlCrossRefPubMed
  25. 25.↵
    1. Iwanami S,
    2. Kakizoe Y,
    3. Morita S,
    4. Miura T,
    5. Nakaoka S,
    6. Iwami S
    (2017) A highly pathogenic simian/human immunodeficiency virus effectively produces infectious virions compared with a less pathogenic virus in cell culture. Theor Biol Med Model 14: 9. doi:10.1186/s12976-017-0055-8
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Iwami S,
    2. Takeuchi JS,
    3. Nakaoka S,
    4. Mammano F,
    5. Clavel F,
    6. Inaba H,
    7. Kobayashi T,
    8. Misawa N,
    9. Aihara K,
    10. Koyanagi Y, et al.
    (2015) Cell-to-cell infection by HIV contributes over half of virus infection. ELife 4: e08150. doi:10.7554/eLife.08150
    OpenUrlCrossRefPubMed
  27. 27.↵
    1. Ziegler CGK,
    2. Allon SJ,
    3. Nyquist SK,
    4. Mbano IM,
    5. Miao VN,
    6. Tzouanas CN,
    7. Cao Y,
    8. Yousif AS,
    9. Bals J,
    10. Hauser BM, et al.
    (2020) SARS-CoV-2 receptor ACE2 is an interferon-stimulated gene in human airway epithelial cells and is detected in specific cell subsets across tissues. Cell 181: 1016–1035.e19. doi:10.1016/j.cell.2020.04.035
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. Eichler H-G,
    2. Müller M
    (1998) Drug distribution. Clin Pharmacokinet 34: 95–99. doi:10.2165/00003088-199834020-00001
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Iwanami S,
    2. Ejima K,
    3. Kim KS,
    4. Noshita K,
    5. Fujita Y,
    6. Miyazaki T,
    7. Kohno S,
    8. Miyazaki Y,
    9. Morimoto S,
    10. Nakaoka S, et al.
    (2021) Detection of significant antiviral drug effects on COVID-19 with reasonable sample sizes in randomized controlled trials: A modeling study. PLoS Med 18: e1003660. doi:10.1371/journal.pmed.1003660
    OpenUrlCrossRef
  30. 30.↵
    1. Kim KS,
    2. Ejima K,
    3. Iwanami S,
    4. Fujita Y,
    5. Ohashi H,
    6. Koizumi Y,
    7. Asai Y,
    8. Nakaoka S,
    9. Watashi K,
    10. Aihara K, et al.
    (2021) A quantitative model used to compare within-host SARS-CoV-2, MERS-CoV, and SARS-CoV dynamics provides insights into the pathogenesis and treatment of SARS-CoV-2. PLoS Biol 19: e3001128. doi:10.1371/journal.pbio.3001128
    OpenUrlCrossRef
  31. 31.↵
    1. Hoffmann M,
    2. Kleine-Weber H,
    3. Schroeder S,
    4. Krüger N,
    5. Herrler T,
    6. Erichsen S,
    7. Schiergens TS,
    8. Herrler G,
    9. Wu NH,
    10. Nitsche A, et al.
    (2020) SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor. Cell 181: 271–280.e8. doi:10.1016/j.cell.2020.02.052
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. Zou L,
    2. Ruan F,
    3. Huang M,
    4. Liang L,
    5. Huang H,
    6. Hong Z,
    7. Yu J,
    8. Kang M,
    9. Song Y,
    10. Xia J, et al.
    (2020) SARS-CoV-2 viral load in upper respiratory specimens of infected patients. New Engl J Med 382: 1177–1179. doi:10.1056/NEJMc2001737
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Wölfel R,
    2. Corman VM,
    3. Guggemos W,
    4. Seilmaier M,
    5. Zange S,
    6. Müller MA,
    7. Niemeyer D,
    8. Jones TC,
    9. Vollmar P,
    10. Rothe C, et al.
    (2020) Virological assessment of hospitalized patients with COVID-2019. Nature 581: 465–469. doi:10.1038/s41586-020-2196-x
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Pan Y,
    2. Zhang D,
    3. Yang P,
    4. Poon LLM,
    5. Wang Q
    (2020) Viral load of SARS-CoV-2 in clinical samples. Lancet Infect Dis 20: 411–412. doi:10.1016/S1473-3099(20)30113-4
    OpenUrlCrossRefPubMed
  35. 35.↵
    1. Zheng S,
    2. Fan J,
    3. Yu F,
    4. Feng B,
    5. Lou B,
    6. Zou Q,
    7. Xie G,
    8. Lin S,
    9. Wang R,
    10. Yang X, et al.
    (2020) Viral load dynamics and disease severity in patients infected with SARS-CoV-2 in Zhejiang province, China, January-march 2020: Retrospective cohort study. BMJ 10: m1443. doi:10.1136/bmj.m1443
    OpenUrlCrossRef
  36. 36.↵
    1. Kuba K,
    2. Imai Y,
    3. Rao S,
    4. Gao H,
    5. Guo F,
    6. Guan B,
    7. Huan Y,
    8. Yang P,
    9. Zhang Y,
    10. Deng W, et al.
    (2005) A crucial role of angiotensin converting enzyme 2 (ACE2) in SARS coronavirus-induced lung injury. Nat Med 11: 875–879. doi:10.1038/nm1267
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Oudit GY,
    2. Kassiri Z,
    3. Jiang C,
    4. Liu PP,
    5. Poutanen SM,
    6. Penninger JM,
    7. Butany J
    (2009) SARS-coronavirus modulation of myocardial ACE2 expression and inflammation in patients with SARS. Eur J Clin Invest 39: 618–625. doi:10.1111/j.1365-2362.2009.02153.x
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Gautret P,
    2. Lagier JC,
    3. Parola P,
    4. Hoang VT,
    5. Meddeb L,
    6. Mailhe M,
    7. Doudier B,
    8. Courjon J,
    9. Giordanengo V,
    10. Vieira VE, et al.
    (2020) Hydroxychloroquine and azithromycin as a treatment of COVID-19: Results of an open-label non-randomized clinical trial. Int J Antimicrob Agents 56: 105949. doi:10.1016/j.ijantimicag.2020.105949
    OpenUrlCrossRefPubMed
  39. 39.↵
    1. Tang W,
    2. Cao Z,
    3. Han M,
    4. Wang Z,
    5. Chen J,
    6. Sun W,
    7. Wu Y,
    8. Xiao W,
    9. Liu S,
    10. Chen E, et al.
    (2020) Hydroxychloroquine in patients with mainly mild to moderate coronavirus disease 2019: Open label, randomised controlled trial. BMJ 369: m1849. doi:10.1136/bmj.m1849
    OpenUrlAbstract/FREE Full Text
  40. 40.↵
    1. Bian H,
    2. Zheng Z-H,
    3. Wei D,
    4. Wen A,
    5. Zhang Z,
    6. Lian J-Q,
    7. Kang W-Z,
    8. Hao C-Q,
    9. Wang J,
    10. Xie R-H, et al.
    (2021) Safety and efficacy of meplazumab in healthy volunteers and COVID-19 patients: A randomized phase 1 and an exploratory phase 2 trial. Signal Transduct Target Ther 6: 194. doi:10.1038/s41392-021-00603-6
    OpenUrlCrossRef
  41. 41.↵
    1. Bi Q,
    2. Wu Y,
    3. Mei S,
    4. Ye C,
    5. Zou X,
    6. Zhang Z,
    7. Liu X,
    8. Wei L,
    9. Truelove SA,
    10. Zhang T, et al.
    (2020) Epidemiology and transmission of COVID-19 in 391 cases and 1286 of their close contacts in shenzhen, China: A retrospective cohort study. Lancet Infect Dis 10: 1016–1025. doi:10.1016/S1473-3099(20)30287-5
    OpenUrlCrossRefPubMed
  42. 42.↵
    1. Goyal A,
    2. Cardozo-Ojeda EF,
    3. Schiffer JT
    (2020) Potency and timing of antiviral therapy as determinants of duration of SARS-CoV-2 shedding and intensity of inflammatory response. Sci Adv 6: eabc7112. doi:10.1126/sciadv.abc7112
    OpenUrlFREE Full Text
  43. 43.↵
    1. Neant N,
    2. Lingas G,
    3. Le Hingrat Q,
    4. Ghosn J,
    5. Engelmann I,
    6. Lepiller Q,
    7. Gaymard A,
    8. Ferré V,
    9. Hartard C,
    10. Plantier JC, et al.
    (2021) Modeling SARS-CoV-2 viral kinetics and association with mortality in hospitalized patients from the French COVID cohort. Proc Natl Acad Sci U S A 118: e2017962118. doi:10.1073/pnas.2017962118
    OpenUrlAbstract/FREE Full Text
  44. 44.↵
    1. Padmanabhan P,
    2. Desikan R,
    3. Dixit NM
    (2020) Targeting TMPRSS2 and Cathepsin B/L together may be synergistic against SARS-CoV-2 infection. PLoS Comput Biol 16: e1008461. doi:10.1371/journal.pcbi.1008461
    OpenUrlCrossRef
  45. 45.↵
    1. Dobrovolny HM
    (2020) Quantifying the effect of remdesivir in rhesus macaques infected with SARS-CoV-2. Virology 550: 61–69. doi:10.1016/j.virol.2020.07.015
    OpenUrlCrossRef
  46. 46.↵
    1. Goyal A,
    2. Duke ER,
    3. Cardozo-Ojeda EF,
    4. Schiffer JT
    (2020) Mathematical modeling explains differential SARS CoV-2 kinetics in lung and nasal passages in remdesivir treated rhesus macaques. BioRxiv doi:10.1101/2020.06.21.163550(Preprint posted June 22, 2020).
    OpenUrlAbstract/FREE Full Text
  47. 47.↵
    1. Ikeda H,
    2. Nakaoka S,
    3. de Boer RJ,
    4. Morita S,
    5. Misawa N,
    6. Koyanagi Y,
    7. Aihara K,
    8. Sato K,
    9. Iwami S
    (2016) Quantifying the effect of Vpu on the promotion of HIV-1 replication in the humanized mouse model. Retrovirology 13: 23. doi:10.1186/s12977-016-0252-2
    OpenUrlCrossRef
  48. 48.↵
    1. Samson A,
    2. Lavielle M,
    3. Mentré F
    (2006) Extension of the SAEM algorithm to left-censored data in nonlinear mixed-effects model: Application to HIV dynamics model. Comput Stat Data Anal 51: 1562–1574. doi:10.1016/j.csda.2006.05.007
    OpenUrlCrossRef
  49. 49.↵
    1. Efron B,
    2. Stein C
    (1981) The Jackknife estimate of variance. Ann Statist 9: 586–596. doi:10.1214/aos/1176345462
    OpenUrlCrossRef
  50. 50.↵
    1. Shao J,
    2. Wu CFJ
    (1989) A general theory for Jackknife variance estimation. Ann Statist 17: 1176–1197. doi:10.1214/aos/1176347263
    OpenUrlCrossRef
  51. 51.↵
    1. Baccam P,
    2. Beauchemin C,
    3. Macken CA,
    4. Hayden FG,
    5. Perelson AS
    (2006) Kinetics of influenza A virus infection in humans. J Virol 80: 7590–7599. doi:10.1128/JVI.01623-05
    OpenUrlAbstract/FREE Full Text
  52. 52.↵
    1. Goncalves A,
    2. Bertrand J,
    3. Ke R,
    4. Comets E,
    5. de Lamballerie X,
    6. Malvy D,
    7. Pizzorno A,
    8. Terrier O,
    9. Rosa Calatrava M,
    10. Mentré F, et al.
    (2020) Timing of antiviral treatment initiation is critical to reduce SARS-CoV-2 viral load. CPT Pharmacometrics Syst Pharmacol 9: 509–514. doi:10.1002/psp4.12543
    OpenUrlCrossRef
  53. 53.↵
    1. Callaway DS,
    2. Perelson AS
    (2002) HIV-1 infection and low steady state viral loads. Bull Math Biol 64: 29–64. doi:10.1006/bulm.2001.0266
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Download PDF
Email Article

Thank you for your interest in spreading the word on Life Science Alliance.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Incomplete antiviral treatment may induce longer durations of viral shedding during SARS-CoV-2 infection
(Your Name) has sent you a message from Life Science Alliance
(Your Name) thought you would like to see the Life Science Alliance web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Incomplete antiviral treatment on SARS-CoV-2 infection
Kwang Su Kim, Shoya Iwanami, Takafumi Oda, Yasuhisa Fujita, Keiji Kuba, Taiga Miyazaki, Keisuke Ejima, Shingo Iwami
Life Science Alliance Aug 2021, 4 (10) e202101049; DOI: 10.26508/lsa.202101049

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Incomplete antiviral treatment on SARS-CoV-2 infection
Kwang Su Kim, Shoya Iwanami, Takafumi Oda, Yasuhisa Fujita, Keiji Kuba, Taiga Miyazaki, Keisuke Ejima, Shingo Iwami
Life Science Alliance Aug 2021, 4 (10) e202101049; DOI: 10.26508/lsa.202101049
Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
Issue Cover

In this Issue

Volume 4, No. 10
October 2021
  • Table of Contents
  • Cover (PDF)
  • About the Cover
  • Masthead (PDF)
Advertisement

Jump to section

  • Article
    • Abstract
    • Introduction
    • Results
    • Discussion
    • Materials and Methods
    • Acknowledgements
    • References
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF

Subjects

  • Systems & Computational Biology

Related Articles

  • No related articles found.

Cited By...

  • No citing articles found.
  • Google Scholar

More in this TOC Section

  • Human GIGYF2 and GRB10 do not interact
  • hnRNP A1 binds a novel ESS element inducing CFTR missplicing
  • Mitochondria form inhibition granules
Show more Research Article

Similar Articles

EMBO Press LogoRockefeller University Press LogoCold Spring Harbor Logo

Content

  • Home
  • Newest Articles
  • Current Issue
  • Archive
  • Subject Collections

For Authors

  • Submit a Manuscript
  • Author Guidelines
  • License, copyright, Fee

Other Services

  • Alerts
  • Bluesky
  • X/Twitter
  • RSS Feeds

More Information

  • Editors & Staff
  • Reviewer Guidelines
  • Feedback
  • Licensing and Reuse
  • Privacy Policy

ISSN: 2575-1077
© 2025 Life Science Alliance LLC

Life Science Alliance is registered as a trademark in the U.S. Patent and Trade Mark Office and in the European Union Intellectual Property Office.