Abstract
Objective Investigate the impact of using antiviral drugs to control local outbreaks of COVID-19.
Methods Using a simulation-based model of viral transmission we tested the impact of different intervention measures for the control of COVID-19.
Results The use of an antiviral drug, in combination with contact tracing, quarantine and isolation, results in a significant decrease of the mean final size and the peak incidence of local outbreaks of COVID-19, provided delays in contact tracing are small.
Conclusions Integrating antiviral drugs together with contact tracing and quarantine is predicted through this model to be an effective tool for the control of local outbreaks of COVID-19.
1 Introduction
To control local outbreaks of COVID-19 we investigate the use of contact tracing and isolation in combination with an antiviral compound. Even when perfect isolation is in place, it may not be sufficient to contain a local COVID-19 outbreak [5]. Therefore, in the absence of a vaccine, an antiviral drug in addition to isolation could be used to contain the current COVID-19 epidemic. Today there are no Corona-specific drugs and development of potent and safe drugs typically takes years. However, there are a number of drugs, originally targeted towards other viral infections, in clinical trials for their ability to control SARS-CoV-2 infection. Most of these drugs inhibit key components of the coronavirus infection life cycle including viral entry, replication, RNA synthesis and protein synthesis. We therefore assume that an antiviral compound will reduce the viral load of an infected individual with COVID-19. For our modeling experiments we searched for an experimental drug for which viral load data were available to inform our model. Remdesivir, was recently shown to inhibit SARS-CoV-2 in vitro [13] and is currently under evaluation in clinical trials. Remdesivir is an investigational broad-spectrum antiviral agent that has been developed against the Ebola-virus and which functions as a nucleotide-analog inhibitor of the viral RNA-dependent RNA polymerase and has activity against a wide range of RNA-viruses [14]. In addition, it is also active against SARS-CoV and MERS-CoV, which is expected given the high level of similarity between the genes that are involved in the replication cycle of these corona viruses. For the aforementioned reasons we inform our model with data on the control of MERS-CoV viral load by Remdesivir in a translational murine model [11]. The animal model that was utilized was specifically developed to better approximate the pharmacokinetics and drug exposure profile in humans. Therefore, the measure of viral titers in lung tissue at different time-points in this model serves as a reasonable proxy for viral dynamics upon compound exposure in the controlled setting of a viral challenge. To this end, we calibrate the model to represent the viral load decrease thereof.
In this manuscript we first present the effect of isolation, considering both home quarantine (for individuals that are part of a contact trace network and for infected individuals with mild symptoms) and hospital isolation (for severe cases). We argue that when an individual is quarantined at home, this will only result in a partial reduction of contacts, accounting for contacts with household members and other isolation imperfectness. To compensate for this imperfect isolation, we consider the use of an antiviral compound. We test these different control measures in a simulation study that aims at representing, given the available information, the current COVID-19 epidemic.
2 Methods
2.1 Epidemiological model
The disease dynamics are depicted in the left panel of Figure 1. The possible transitions between epidemic classes are described by the arrows.
Possible transitions among the different epidemic compartments (left panel). A measure describing how infectiousness is distributed over time since infection (right panel).
Individuals are initially susceptible (S) and once infected, they enter the exposed class (E). The infection, that is at first asymptomatic (Ia), can lead to the onset of mild (Im) or severe symptoms (Is). Symptomatic individuals are hospitalized (H), where they are isolated, or are confined in home quarantine (Q), based on the severity of symptoms. Ultimately, all infectives are assumed to either recover from infection or die (R). Individuals that are hospitalized are immediately isolated; therefore they can no longer transmit the disease. The quarantined individuals, can still make contacts, although at a decreased rate.
The transition from the susceptible to the exposed class is governed by a stochastic process based on the notion of infectious contact processes [10]. First, contacts between individuals are generated. When such contacts generated between susceptible and infectious people, these can result in an infection event according to a Bernoulli probability value based on the time since infection. This probability is computed, at a precise time point, as the product of two components: the infectiousness measure, ν(t), which quantifies the level of infectiousness over time, and the total amount of infectivity q, i.e. the number of infections over the contact rate [12]. The function ν(t) is defined over the exposed and infectious period, or analogously over the incubation and symptomatic period, along which it integrates to one. This function is scaled to have a similar shape among different infectives, based on their lengths of exposed and infectious period. According to this framework, an infectious individual makes effective contacts at a rate, r(t) given by:
where λ is the contact rate. The mean number of effective contacts is an approximation of the reproduction number. The two quantities are identical in an infinite and homogeneous population, where the probability of making two effective contacts with the same person is zero. For the considered population size, the probability of this event is extremely low. Therefore, throughout the manuscript, we approximate the reproduction number with the mean number of secondary cases.
In the considered framework, isolation/quarantine is implemented by reducing the contact rate λ at the time of diagnosis. To date, little is known about the difference in viral load among severe and mild cases. Zhou et al. [19] indicate that in nasal and throat swabs the viral load is higher in mild cases. However, virus was detected also in blood and stool sample [9]. Due to this uncertainty, we assume that the same curve is defined for all the infected individuals. In this work, we assume the population to be homogeneous, closed and finite population. The former assumption was chosen because of the limited knowledge on the COVID-19 determinants as mentioned above. The two latter assumptions relate to the control measures currently in place, e.g. in Italy, aiming at containing immigration and emigration in a country with an ongoing outbreak.
2.2 Simulation parameters and distributions
In Table 1 we report the parameters and distributions that were utilized in the simulation study. Where distributions are not reported, the parameters are assumed to be constant. In the last column, when available, we report the references to the literature that justifies the choice of the parameter value, or distribution, we use.
Model Parameters
We assume infectious individuals to make, on average, between two and three effective contacts. This value, is set accordingly to the current estimate of the basic reproduction number [2; 7; 16; 17]. The infectiousness measure is set to represent the viral load observations reported in [9; 18; 19], the peak of which is reached within a few days after symptoms onset. The total length is modeled as the convolution of incubation and symptomatic period. Incubation and symptomatic period, are set, respectively, to have a mean length of 5.2 days and 18 days. There is no precise estimate of the length of the symptomatic period to date, thus we choose it according to the aggregated nasal and throat swabs data [19]. The time to hospitalization is estimated from the data presented in [1]. We assume that the time to hospitalization coincides with the time of diagnosis. At this time-point, depending on the severity of the symptoms, individuals are isolated (severe cases) or quarantined (mild cases). The population size is set to 500 to represent a localised outbreak of COVID-19. Furthermore, we assume that the contact tracing starts when individuals are diagnosed.
2.3 Contact tracing and isolation
In order to implement contact tracing we keep track of a contact history Hi for each individual i for all contacts made since the time of infection. When an individual i is found to be infected with SARS-CoV-2, a contact tracing procedure is started. We assume that each contact in Hi will be traced back successfully with probability η. Depending on the considered scenario, traced-back individuals will be monitored, or put in quarantine/isolation if positive for SARS-Covid-2, or injected with the antiviral drug and put in quarantine/isolation if positive for SARS-Covid-2. We assume that traced-back individuals that test positive are isolated in the 16% of cases. The quarantine will result in a decreased contact rate (i.e. imperfect isolation), λq, while in case of perfect isolation the contact rate is set to zero. Similarly, diagnosed individuals will also be quarantined: at home (mild symptoms), with a decreased contact rate λq, or in the hospital (severe symptoms), where we assume that perfect isolation is possible. We assume that 16% of infections is severe and require hospitalisation [4].
2.4 Antiviral compounds
To compensate for imperfect isolation we investigate the use of antiviral compounds to reduce the infectiousness of an infected individual. We assume that, once the antiviral compound has been administered, the infectiousness measure will exponentially decay according to an inverse Malthusian growth model (shown in Figure 2) [8]. The rate of this decay is set to represent the reduction in viral load, due to Remdesevir, as reported in [11] for the MERS coronavirus.
The blue and the orange lines describe the infectiousness measure, respectively, before (dashed blue) and after antiviral injection (solid yellow). The red arrows indicate the injection times.
3 Scenarios
We assume the following parameters for the reduction of contacts because of home quarantine: λq = 0.1λ, 0.25λ, 0.5λ and for the probability of tracing back a contact in history Hi: η = 0.25, 0.5, 0.75 [6]. Scenarios: In all the considered scenarios we assume that individuals are isolated, or quarantined, when diagnosed. Moreover, we assume that contact tracing starts at the time of diagnosis.
IAS: Traced individuals are monitored for two weeks, and isolated/quarantined if they show symptoms during this period. This scenario is similar to the baseline scenario described by Hellewell et al. [5] with the exception that in our description only severe cases are isolated while the mild are home quarantined. This scenario reflects more realistically the current practice of containment.
IBS: Traced individuals are isolated/quarantined, as soon as they test positive for SARS-COV-2. We assume that an individual that is infected tests positive 2 days after infection. Therefore, a traced individual is tested immediately when traced, and, if this test was negative, we test the individual again two days later.
IBTBS: Diagnosed patient is immediately treated with the antiviral drug. Furthermore, traced individuals are isolated/quarantined and injected with the antiviral drug, as soon as they test positive for SARS-COV-2. We assume that an individual that is infected tests positive 2 days after infection. Therefore, a traced individual is tested immediately when traced, and, if this test was negative, we test the individual again two days later.
For each scenario we run 5000 simulations. Among these, we compute the mean final size and the cases at peak for the one in which at least the 10% of individuals have been infected. Doing this, we only account for outbreaks that are most challenging to contain.
4 Results
Quarantine, isolation and antiviral treatments lead, in different levels, to the mitigation of the outbreak by reducing the final size as well as by reducing the number of cases at the peak of the epidemic. The containment performance depends, among all the scenarios, on the probability to successfully trace contacts and on the reduction in contact rate due to quarantine (Figure 3). Isolation and quarantine lead to a substantial decrease in mean final size and peak incidence. When performed prior to symptom onset their efficacy increases, which is important, as there is a positive probability of presymptomatic infection [3]. The antiviral treatment is shown to have a substantial impact and, together with quarantine and isolation, significantly reduces the mean final size, the peak incidence and the number of outbreaks that are most challenging to contain.
Left panel: distributions of the mean final size value for Scenario IAS (yellow), Scenario IBS (green) and Scenario IBTBS (blue) when the quarantine contact rate is λq = 0.25λ together with the probability that a simulation leads to a number of cases smaller than the 10% of the population (purple asterisks). Right panel: Peak incidence together with the 2.5 % and the 97.5 % percentiles.
5 Discussion
We assume that we have sufficient antiviral drugs doses to treat all individuals that are encountered via the contact tracing procedure. This is motivated by the fact that we consider an emerging outbreak and the required number of doses will thus be limited.
Furthermore, we assume that all the individuals will show symptoms, sooner or later, during their infectious period, and therefore all infected individuals will be diagnosed. Due to the awareness of COVID-19 given by media and government officials, individuals are more likely to act upon even mild symptoms. This assumption is in line with the work of Hellewell et al. [5]. Moreover, we assume that all infected individuals will eventually be diagnosed. While this is a limitation of our study, we argue that in the scenarios where we aggressively trace and treat the contacts of individuals, we are more likely to find (and constrain) cases that would otherwise go undetected.
In the hospital, we assume perfect isolation, meaning that infected individuals cannot spread the infection. However, health care workers are at risk and they could be infected by infected individuals in isolation [5].
In the IBS and IBTBS scenarios we assume that the traced individuals that test positive are isolated in the 16% of cases, even before showing actual severe symptoms. This is a simplification; a more realistic implementation would be to first quarantine cases and isolated them when the severe symptoms onset. We expect that our implementation may only slightly affect the outcomes, without influencing the conclusion that antivirals, in addition to isolation and quarantine, are effective to control an epidemic. Furthermore, in the sensitivity we tested analysis a scenario in which infection can be detected only after 4 days since infection, close to the incubation period, limiting the impact of a early isolation. In this case, we also highlight the importance of antiviral drugs.
Although this model is informed with data on the control of MERS-CoV viral load using prophylaxis with Remdesivir, it stands to reason that different classes of viral inhibitors control the viral load in different ways. Additionally, despite the sequence similarity of MERS-CoV and SARS-COV-2 it remains to be established whether the impact of Remdesivir (or other antivirals) on the viral load is similar. To this end, longitudinal data of the viral load on COVID-19 infected patients treated with different viral inhibitors will be informative. Furthermore, Sheahan et al. [11] demonstrated that the degree of the clinical benefit of Remdesivir for MERS depends on the viral dose and also on the timing of the treatment of the viral inhibitor. A lethal viral dose and delaying the initiation of antiviral treatment failed to fully prevent viral pathogenesis. Although Remdesivir proved effective at reducing the viral load also in these conditions, the argument for an early start of antiviral treatment is evident. Presumably, reducing the viral load with an antiviral compound loses its efficacy in advanced disease as the tissue damage is sustained by inflammatory processes in absence of the viral initiator. In our implementation, antiviral injections are immediately administered to successfully traced back individuals, mostly in their asymptomatic phase, and to the diagnosed patient. Therefore, we believe that the assumptions on the use of this drug, in the considered scenario, are reasonable.
6 Conclusion
The ongoing epidemic of COVID-19 threatens the health system of many countries. Control strategies are fundamental measures that public health officials must implement to contain this epidemic. Isolation and quarantine can be of great importance, however, relying on their exclusive use could fail to contain an ongoing outbreak. In addition, when several infected individuals need hospital care, the number of cases at peak should be minimized as much as possible to avoid that regions run out of hospital capacity. With this study we highlight the impact of a potential antiviral compound that reduces the viral load and, consequently, the infectiousness of infectives. We describe the action of Remdesivir, however, the proposed investigation can be extended to other antivirals, with the caveat that the efficacy of administering an antiviral compound, in addition to isolation and quarantine, depends on the effectiveness of the respective drug. We show that the use of this compound leads to a substantial reduction of the mean final size and the peak incidence. In addition, the number of outbreaks that are most challenging to contain decreases when the antiviral is administered to diagnosed and traced individuals. Therefore, the administration of an antiviral drug, together with isolation and quarantine, is expected to have a major impact in the control of local COVID-19 outbreaks.
We remain hopeful that the ongoing clinical trials will reveal an antiviral compound that can be used as a treatment and prophylaxis. Yet, our work shows that such compounds have a great potential to mitigate pandemics and research towards new drugs targeting different virus families with pandemic potential is thus warranted.
Data Availability
Our study concerns the modelling of antiviral compounds to contain COVID-19 epidemics. Our modelling code has been posted to github.
Author contributions
A.T. contributed to the conceptualization of the study, the construction and implementation of the mathematical model, the experimental setting and the writing of the manuscript. P.L. contributed to the conceptualization of the study, the construction of the mathematical model, the experimental setting and the writing of the manuscript. J.V. contributed to the construction of the mathematical model, the experimental setting and the writing of the manuscript. A.V. contributed to the conceptualization of the study and the writing of the manuscript. J.N. contributed to the conceptualization of the study and the writing of the manuscript. N.H. contributed to the conceptualization of the study, the construction of the mathematical model, the experimental setting, the writing of the manuscript and the project supervision.
Competing interests
Besides his employment at the Hasselt University, JV is employed at Bioqube Ventures. Bioqube Ventures was not involved in this work, nor does it prosper financially as a result of the current study. The other authors declare that they have no competing interests.
Materials and Correspondence
Correspondence and material requests should be addressed to Prof. Dr. Niel Hens. Source code of our model and experiments can be found on https://github.com/AndreaTorneri/ViralTransm.
7 Sensitivity Analysis
We report in Figure 4 and 5 the sensitivity analysis for the quarantine contact rate: λq = 0.1λ, 0.5λ. The introduction of an antiviral compound substantially contributes in reducing the mean final size, the peak incidence and the probability of a challenging outbreak in all the considered settings. This decrease, compared to the scenario in which only isolation/quarantine is implemented, increases when quarantine is less effective (left panels). In Figure 6 and 7, we vary the reproduction number that is set, respectively, to R0 = 2 and R0 = 3. The effect of the antiviral drug, in addition to isolation and quarantine, increases when the reproduction number increases. In case of R0 = 3, the peak incidence decreases when the antiviral compound is used, compared to control measures based only on isolation and quarantine after symptoms onset.
Distributions of the mean final size value for Scenario IAS (yellow), Scenario IBS (green) and Scenario IBTBS (blue) when the quarantine contact rate is λq = 0.5 (left panel) and λq = 0.1 (right panel) together with the probability that a simulation leads to a number of cases smaller than the 10% of the population (purple asterisks)
Mean peak incidence for Scenario IAS (yellow), Scenario IBS (green), Scenario IATBS (red) and Scenario IBTBS (blue) when the quarantine contact rate is λq = 0.1λ (left panel) and λq = 0.5λ (right panel) together with 2.5% and 97.5% percentiles.
Mean final size distribution (left panel) and mean peak incidence for Scenario IAS (yellow), Scenario IBS (green) and Scenario IBTBS (blue) when the quarantine contact rate is λq = 0.25λ, and R0 = 2. In the left panel, for each scenario we report the probability that a simulation leads to a number of cases smaller than the 10% of the population (purple asterisks). In the right panel, together with the point estimates we report the 2.5% and 97.5 % percentiles.
Mean final size distribution (left panel) and mean peak incidence for Scenario IAS (yellow), Scenario IBS (green) and Scenario IBTBS (blue) when the quarantine contact rate is λq = 0.25λ, and R0 = 3. In the left panel, for each scenario we report the probability that a simulation leads to a number of cases smaller than the 10% of the population (purple asterisks). In the right panel, together with the point estimates we report the 2.5% and 97.5 % percentiles.
In Figure 8 we investigate the effect of a longer time needed for the test to detect an infectious individual.We assume the test is positive when performed on an infectious individual who has been infected since at least 4 days. Simulations show a substantial increase, both in the mean final size and the peak incidence for the IBS scenario. Instead, the use of an antiviral drug results also in this case of remarkable impact in both the final size and the peak incidence.
Mean final size distribution (left panel) and mean peak incidence for Scenario IAS (yellow), Scenario IBS (green) and Scenario IBTBS (blue) when the quarantine contact rate is λq = 0.25λ, R0 = 2.5 and the test detect positively an infectious individual after 4 days since infection. In the left panel, for each scenario we report the probability that a simulation leads to a number of cases smaller than the 10% of the population (purple asterisks). In the right panel, together with the point estimates we report the 2.5% and 97.5 % percentiles.
Acknowledgements
This work is funded by the Epipose project from the European Union’s SC1-PHE-CORONAVIRUS-2020 programme, project number 101003688 and by the European Union’s Horizon 2020 research and innovation programme (grant agreement 682540 — TransMID). A.T. acknowledges support from the special research fund of the University of Antwerp. P.L. was supported by funding from the Flemish Government under the “Onderzoeksprogramma Artificiële Intelligentie (AI) Vlaanderen” programme.
Footnotes
↵† joint first author
