Abstract
Despite the development of safe and effective vaccines, effective treatments for COVID-19 disease are still desperately needed. Several antiviral drugs have shown to be effective in reducing hospitalizations, and it is expected Despite the development of safe and effective vaccines, effective treatments for COVID-19 disease are still desperately needed. Several antiviral drugs have shown to be effective in reducing hospitalizations, and are to be widely available in the coming months. In the present work, we use an agent-based mathematical model to assess the potential population impact of the use of antiviral treatments in locations in four countries with different demographic structure and current levels of vaccination coverage: Kenya, Mexico, United States (US) and Belgium. For each location, we varied daily treatment initiation rate (DTIR) and antiviral effect in reducing transmission (AVT). Irrespective of location and AVT, widespread antiviral treatment of symptomatic adult infections (≥ 20% DTIR) prevented the majority of COVID-19 deaths, and recruiting 6% of all adult symptomatic infections daily reduced mortality by a third in all countries. Furthermore, we show that targeting antiviral treatment to the oldest age group (65 years old and older, DTIR of 20%) can prevent over 47% of deaths. Our results suggest that early antiviral treatment (within two days of symptoms) is needed to mitigate transmission, preventing 50% more infections compared to late treatment (started 3 to 5 days after symptoms onset). Our results highlight the synergistic effect of vaccination and antiviral treatment: as vaccination rate increased, antiviral treatment had a bigger impact on overall transmission. These results suggest that antiviral treatments can become a strategic tool that, in combination with vaccination, can significantly reduce COVID-19 hospitalizations and deaths and can help control SASRS-CoV-2 transmission.
Introduction
With over 5 million deaths worldwide [1], the COVID-19 pandemic has proven difficult to contain. Despite the development of many safe and effective vaccines, controlling SARS-CoV-2 transmission has shown to be elusive for several reasons, including vaccine supply shortages in low- and middle-income countries [2], vaccine hesitancy [3], and the emergence of new variants [4]. Indeed, the delta and omicron variants, that emerged in summer and fall of 2021, quickly became the predominant strains and have caused large epidemic outbreaks, even in highly vaccinated regions [1]. Rapidly producing such COVID-19 vaccines has been an amazing scientific endeavor, but effective tools to treat COVID-19 disease are still desperately needed. Several antiviral treatments have been found to reduce hospitalizations by 30% to 89% [5–8] when taken within the first five days after developing symptoms, and others are being developed [9]. Seven of these are currently approved for COVID-19 treatment in the United States (US) [10]. Further, earlier studies have shown that these antiviral treatments significantly reduced the amount of infectious virus in treated individuals [11]. Experience with other infectious diseases, such as influenza and HIV, has shown that potent antiviral treatments can be an effective means of prevention and epidemic control [12, 13]. Hence, the advent of effective antiviral drugs raises the possibility that in treating infected individuals we may be prevented onward transmission (indirect population benefit) while also protecting the treated person from severe disease (direct benefit). In this work, we use an agent-based mathematical model to evaluate the potential population impact of widespread use of antiviral treatments reducing hospitalization risk and viral load. We explored the use of antiviral treatments in four locations in different countries (Kenya, Mexico, US and Belgium) with very different demographic composition and vastly different proportions of vaccinated individuals under epidemic waves with characteristics similar to the delta and omicron waves. Additionally, we analyzed the impact of antiviral treatment on overall transmission by assuming that antiviral treatment would reduce viral transmission of treated individuals by 25, 50, 75 or 100%. We showed that the synergistic use of vaccine and antiviral treatments can significantly reduce the burden of COVID-19 and help restore normalcy. Further, our model suggested that targeted use of antiviral treatments can be used to prevent the majority of deaths.
Results
Briefly, we used COVASIM [14], a previously developed agent-based model of SARS-CoV-2 spread. Our model simulates a population of 500,000 people interacting through a network over the course of a 6 month time horizon, where each individual in the population is an agent. Every day, individuals get in contact with others in four possible locations: “home”, “school”, “work” or “community”. At a given point in time, individuals can be susceptible, infected asymptomatic, pre-symptomatic or symptomatic, and recovered. Symptomatic infected individuals can develop a mild, severe or critical infection after symptoms onset. Age-specific proportions of infections are assumed to progress without developing symptoms being 30% less infectious than symptomatic ones (Table S1). We assumed that 40% of the population has been infected in previous epidemic waves and protected from re-infection over the study duration, and compared these results to scenarios assuming 20 and 60% pre-existing immunity. We simulated synthetic populations of equal size at four different locations in four different countries: Kenya, Mexico, US and Belgium. For each location, the model uses country-specific demographics to inform population structure and household sizes. Because we are interested in investigating wide-ranging population-level effects of antiviral treatment, we explored two main scenarios: first, we assumed deployment would take place under the omicron wave, representing a high transmissible variant for which vaccines might not very effective at preventing infections. Then, we repeated the analysis assuming antiviral deployment under the delta wave, a less transmissible variant for which vaccines remained moderately effective in preventing infections. For each of these waves, we used variant-specific parameters for viral transmission, vaccination coverage during that wave and vaccine effectiveness (Methods for full details).
We assumed that the antiviral treatment would have two primary effects. First, in line with results from clinical trials, we assumed an antiviral effect on hospitalization, denoted by AVH, by which a course of antiviral treatment would reduce the rate of hospitalization of symptomatic infected individuals by 88% (main results) or 30% (sensitivity analysis) [5–8]. Second, because clinical trials have reported an antiviral effect on viral load, [5, 11], we assumed that the antiviral treatment would have an effect in overall transmission, denoted by AVT. We explored reductions of 25, 50 75 or 100% on transmission for treated individuals.
Because all but one antiviral treatments are currently approved for only treatment [10], we limited our analysis to treatment interventions mimicking the conditions under which the antiviral products were mostly tested in clinical trials: we assumed that antiviral treatment would be given to adult symptomatically infected individuals (denoted below as eligible symptomatic infections) within five days from symptoms onset. We explored scenarios with treatment initiated within five days, within two days or between days 3 and 5 from symptoms onset, separately to study the impact of treatment timing on the results. We also varied the rate at which symptomatic adults could be identified and recruited for treatment, considering scenarios with a daily treatment initiation rate (DTIR) of 2, 4, 6, 8, 10, 20, … 100% of eligible symptomatic adults. Throughout the text, we present the median percentage deaths and infections averted compared to base-case scenario assuming that existing non-pharmaceutical measures remain in place for the next 6 months without additional vaccination across 100 simulations. Cumulative deaths and infections (and respective confidence bounds) are presented in the Supplemental Material.
Antiviral treatment deployment under the omicron wave
Targeted use of antiviral treatment can avert large numbers of deaths
We first examine the use of antiviral treatment in all adults (18 years or older). For all the countries considered, irrespective of the antiviral effect on viral transmission, large numbers of deaths were averted for all scenarios considered, even if the DTIR was low. A 6% DTIR among all adult symptomatic infections averted over 20% of deaths when compared to no antiviral use in all countries (24, 23, 23 and 24% deaths averted for Kenya, Mexico, US and Belgium respectively assuming AVT=25%). Over 58% of deaths would be averted if 20% of all adult symptomatic infections were identified and recruited for treatment daily (58, 58, 59 and 59% of deaths averted for Kenya, Mexico, US and Belgium respectively assuming AVT=25%). If the DTIR increased to 50% or more, over 85% of deaths would be averted in all countries. We observed minimal differences between countries and between assumed antiviral reduction in viral transmission, pointing to the fact that these reductions are a result of direct rather than indirect antiviral treatment protection (Figs. 1 and S1).
Percentage of deaths averted (compared to a baseline of no antiviral treatment) for A) Kenya, B) Mexico, C) United States and D) Belgium. Here, we assumed an epidemic wave with parameters similar to those of the omicron epidemic wave (transmissibility, vaccine effectiveness, and vaccination coverage). For each location, the colors represent four possible values of AVT (25, 50, 75 or 100% reduction in viral transmission in treated symptomatic individuals) and a daily treatment initiation rate (DTIR) of 2-100% of adult symptomatic individuals within the first 5 days of symptoms.
Because identifying large numbers of all the eligible symptomatic adults (individuals aged 18 or older) might be difficult to implement in practice, especially in countries where testing is not widely available, we then evaluated additional strategies where we restricted antiviral treatment to adults over 30, 50 or 65 years of age and considered a 20% DTIR among these populations. Concentrating antiviral treatment in the older age groups can be extremely valuable: despite being the group less populated and with higher vaccination rates (table S4), identifying and treating 20% of the daily symptomatic infections among older adults (aged 65 years old and older) averted 32, 43, 47 and 47% of all deaths in Kenya, Mexico, US and Belgium respectively (assuming 25% reduction in viral load). Adults aged 50 and older represent 26, 40 and 44% of the adult populations in Mexico, US and Belgium respectively, but recruiting 20% of the daily symptomatic infections in this age group averted the majority of deaths in these countries (53, 54 and 56% deaths averted in Mexico, US and Belgium respectively).
More impressively, while this age group represents only 13% of the adult population in Kenya, initiating antiviral treatment in 20% of the symptomatic infections in this age group daily resulted in 47% of deaths averted. Extending the use of antiviral treatment to younger adults resulted in a marginal gain in Mexico, US and Belgium (a maximum 7, 4 and 6% more deaths averted respectively) when compared with targeted treatment in adults older than 50 (Fig. 2).
Percentage of deaths averted (compared to a baseline of no antiviral treatment) for A) Kenya, B) Mexico, C) United States and D) Belgium. Here, we assumed an epidemic wave with parameters similar to those of the omicron epidemic wave (transmissibility, vaccine effectiveness, and vaccination coverage). For each location, the colors represent four possible values of AVT (25, 50, 75 or 100% reduction in viral transmission in treated symptomatic individuals) and targeting the antiviral treatment to symptomatic adults older than 18, 30, 50 or 65 years of age with a daily treatment initiation rate of 20%.
Synergistic effect of antiviral treatment and vaccine can mitigate SARS-CoV-2 transmission
Our model projected that antiviral treatment will have a very limited impact in mitigating transmission in the presence of a variant as transmissible as omicron in countries like Kenya or Mexico, and will have some impact in highly vaccinated countries like the US and Belgium, provided that a large proportion of symptomatic infections are identified and recruited for treatment and the antiviral effect on transmission is high (Figs. 3 and S2). Further, we found a clear synergy of combining antiviral treatment and vaccination. While Belgium and the US have similar population structure, vaccination coverage in Belgium has been much higher than in the US (tables S5 and S6). When combined with antiviral treatment, this difference had a big effect in mitigating transmission: for example, in the most optimistic scenario, assuming that the antiviral treatment completely blocked transmission, a 30% DTIR of adult symptomatic infections averted 9% of overall infections in the US, and twice as many infections were averted in Belgium (18%). In contrast, for this scenario, only 7% of overall infections were averted in Mexico, which had a high vaccination coverage but had not boosted the population at that time, and 5% of the overall infections were averted in Kenya. Indeed with very limited vaccine supply, the effective reproductive number in Kenya is much larger than in the other countries, resulting in epidemics that are more difficult to control (Figs. 3 and 4). As expected, the antiviral effect on viral transmission played a bigger role in controlling SARS-CoV-2 transmission. With a low AVT (25%), a very limited number of infections were averted, irrespective of treatment initiation rate or location, with a maximum of 6% of infections averted, but a maximum of 35% of infections were averted with high AVT (Belgium, Fig. 3D). Note that, for AVT= 25, 50 or 75%, the number of infections averted plateaued for medium and high antiviral treatment coverage. For example, in the US, a DTIR of 50% of adult symptomatic infections with an antiviral treatment with AVT=50% averted 5% of infections, and daily recruiting an additional 20% of eligible symptomatic infections (70% DTIR) resulted in only 2% more infections averted (Fig. 3).
Percentage of infections averted (compared to a baseline of no antiviral treatment) for A) Kenya, B) Mexico, C) United States and D) Belgium. Here, we assumed an epidemic wave with parameters similar to those of the omicron epidemic wave (transmissibility, vaccine effectiveness, and vaccination coverage). For each location, the colors represent four possible values of AVT (25, 50, 75 or 100% reduction in viral transmission in treated symptomatic individuals) and a daily treatment initiation rate (DTIR) of 2-100% of adult symptomatic individuals within the first 5 days of symptoms.
Daily new infections assuming no antiviral treatment (Baseline) or assuming coverage of 10-100% of eligible symptomatic individuals in A) Kenya, B) Mexico, C) United States and D) Belgium. Here, we assumed an epidemic wave with parameters similar to those of the omicron epidemic wave (transmissibility, vaccine effectiveness, and vaccination coverage). For each location, each column represents a different value of AVT (25, 50, 75 or 100% reduction in viral load).
Early treatment is needed to mitigate transmission
We analyzed the difference between early and late treatment by considering treating eligible symptomatic infections within the first two days of symptoms or alternatively between days 3 and 5 after symptom onset. For locations in Kenya and Mexico, late treatment had no impact on transmission (Fig. 5A and B). For locations in the US or Belgium, starting antiviral treatment during the first two days of symptom onset prevented roughly twice as many subsequent infections compared to late treatment (Fig. 5C and D). For example, with a DTIR of 50% of symptomatic infections, and AVT = 75%, 14% and 8% of overall infections were averted in Belgium with early and late treatment respectively (Fig. 5D).
Percentage of deaths averted (compared to a baseline of no antiviral treatment) for A) Kenya, B) Mexico, C) United States and D) Belgium for early (within the first two days of symptoms, left column) or late (between days 2 and 5 of symptoms, right column) treatment. Here, we assumed an epidemic wave with parameters similar to those of the omicron epidemic wave (transmissibility, vaccine effectiveness, and vaccination coverage). For each location, the colors represent four possible values of AVT (25, 50, 75 or 100% reduction in viral transmission in treated symptomatic individuals) and covering 10-100% of eligible symptomatic individuals.
Delta epidemic wave
In this section we present the results assuming antiviral treatment was given under similar conditions as those experienced during the delta wave: we considered vaccination coverages as of October 12th, 2021, more effective vaccines (table S2), and a virus as transmissible as the delta variant (full details in Methods). Under this scenario, the number of deaths averted was slightly higher than for the Omicron wave with an average 2% more deaths averted in all countries (Figs. 1 and S3). Additionally, there were some differences in the number of deaths averted depending on the assumed antiviral effect on transmission, with a maximum of 3, 7, 9 and 12% more deaths averted in Kenya, Mexico, US and Belgium respectively for an antiviral treatment fully blocking transmission compared with one reducing transmission by 25%, suggesting a small indirect effect on the number of deaths averted as a result of reduced transmission. More importantly, under this scenario, we observed that antiviral treatment markedly cut overall viral transmission, and the synergy with vaccination was more apparent. For example, assuming a DTIR of 30% of symptomatic adults and AVT=50%, 5, 12, 16 and 32% of infections were averted in Kenya, Mexico, US and Belgium respectively (compared to 1, 1, 4 and 8% under an omicron wave respectively), Figs. 3 and S4.
Sensitivity analysis
In this section we analyze the sensitivity of our results to key model parameters. All the results presented here assumed an omicron-type epidemic.
Increase or decreased infection-acquired pre-existing immunity
Different places, even within a country, have experienced the COVID-19 pandemic in different ways, depending on local non-pharmaceutical interventions, culture (e.g. vaccine hesitancy, adherence to mask use) and population demographics [15, 16]. Hence, we considered the effect of pre-existing infection-acquired immunity on our results, assuming a lower (20%) or higher (60%) fraction of the population previously infected and currently immune. We observed minimal differences in the number of deaths averted for lower or higher assumed infection-acquired pre-existing immunity. This is expected, as the majority of deaths averted in our model are the result of direct effect of the antiviral treatment. However, the model projected differences in transmission assuming lower or higher infection-acquired pre-existing immunity. Lower infection-acquired immunity resulted in bigger epidemic waves and lower proportions of infections averted with antiviral treatment compared to the main results for all countries, all coverages and all AVEs (maximum 4, 6, 12 and 21% infections averted in Kenya, Mexico, US and Belgium respectively), Figs. S5 left column, S6. In contrast, higher infection-acquired immunity resulted in lower epidemic curves and more infections averted (maximum 13, 18, 41 and 82% infections averted in Kenya, Mexico, US and Belgium respectively), Figs. S5 right column and S7.
Reduced antiviral effectiveness against hospitalization, AVH = 30%
One antiviral treatment (molnupiravir) was shown to have a 30% reduction in hospitalizations and deaths [5], so we repeated our analysis assuming this lower value for AVH. As expected, an antiviral treatment which reduces the hospitalization risk of symptomatically infected by 30% averted considerably less deaths than our main scenario (with AVH = 88%). For example, a DTIR of 20% among adult symptomatic infections resulted in over 20% of the deaths averted for all four countries, a ∼40% reduction compared to the main scenario. Further, under this scenario, there was a clear vaccination effect that was compounded with the antiviral effect on transmission: locations with higher vaccination rates had more deaths averted and as the antiviral effect on transmission increased, the number of deaths averted also increased, Fig. S8.
Decrease in asymptomatic infectiousness
The relative infectivity of asymptomatic individuals remains uncertain [17]. We explored scenarios in which asymptomatic individuals are 50% less infectious than symptomatic ones. The overall antiviral effect on mortality was not sensitive to this parameter: our model projected nearly identical deaths averted for each scenario considered, Fig. S9. However, the number of infections averted assuming a lower infectivity of asymptomatic infections was higher, with higher differences observed in locations with higher vaccination coverages. For example, under this scenario, 20% of the infections were averted in the US assuming a DTIR of 50% and an antiviral treatment fully reducing transmission, compared with 12% infections averted in the main scenario, Fig. S10.
Discussion
Despite the existence of effective vaccines thousands of people continue to suffer from severe COVID infections around the world due to issues with supply, hesitancy and logistics. This highlights the need of effective and affordable treatment options. Several antiviral treatments have shown to be highly efficacious against COVID-19 hospitalization and death, provided that treatment is started early, within 5 days of symptoms onset. In the US, several antiviral treatments have been given emergency use authorization, and 20 million courses of Paxlovid will be delivered to US government in 2022. The Medicines Patent Pool and the manufacturers of molnupiravir and paxlovid (Merck and Pfizer respectively) have announced license agreements to facilitate global access for these drugs [18, 19]. Hence, it is likely that in the next few months antiviral treatments will become widely available globally. In the present work we analyzed the potential population impact of antiviral treatments for reducing SARS-CoV-2 transmission and COVID-19 related hospitalizations and deaths. We considered four equal-sized populations from four different countries (Kenya, Mexico, US and Belgium), representing different population structures and vaccination coverage. We further explored the impact of this intervention under epidemic waves with parameters similar to the delta and omicron waves. Our results suggest that irrespective of location, AVT and variant, widespread use of antiviral treatment (daily treatment initiation rate of 20% of the daily adult symptomatic infections) could prevent over 60% of deaths. Even a small fraction of all adult symptomatic infections can still have a large impact on the expected number of COVID deaths. Indeed, a daily treatment initiation rate of 6% of symptomatic infections averted over 20% of deaths in all countries. Because it might be challenging to identify and test symptomatic infections among all adults, we further explored scenarios restricting antiviral treatment to older age groups. We demonstrated that targeting 20% of the daily symptomatic infections among those 65 years old and older could avert over 45% of deaths in all countries. This highlights the need to better use and promote the use of currently available and potential new antiviral treatments. We projected that under an omicron-like epidemic wave, antiviral treatment will have very limited impact on transmission. However, if newer variants behaved more like the delta variant, then we projected that antiviral treatment can have a larger impact on curbing overall infections. Furthermore, our model suggests that there is a synergistic effect of combining antiviral treatment and vaccination: countries with larger proportion of their populations vaccinated are expected to benefit more from adding antiviral treatment to their pandemic control toolbox. This emphasizes the need of continuous effort to improve vaccination coverage especially in settings where it is extremely low. Finally, our simulations showed that early treatment, when viral transmission is highest [20, 21], is important for mitigating transmission and could be used as a prevention tool, especially in the presence of variants as transmissible as delta, but our model suggests that the use of antiviral treatment to curb transmission will be extremely limited in the presence of more transmissible variants, such as omicron. This is in agreement with previous work [22–26] that has shown the potential use of antiviral therapies to reduce COVID-19 related mortality.
Our model, like all mathematical models, is subject to several limitations. We did not consider the development of antiviral resistance, yet it is possible that if antiviral treatments are widely used, resistance could rapidly develop. Further, as new variants are continuously emerging, antivirals that are highly effective against current variants might become ineffective against future ones [27, 28]. We analyzed the use of antiviral treatment in symptomatic individuals, and did not explore its use in asymptomatic infections or as a prophylaxis. Further, we assumed that an antiviral effect reducing transmission was independent of when treatment was started (as long as it started within the first five days of treatment). In reality, it is possible that antiviral treatments might have different effects in reducing overall transmission depending on when during the course of an infection they are taken, e.g. reducing overall viral load and hence transmission if taken early on but having only a modest effect as time goes by. We assumed that the vaccine had no effect in reducing infectiousness, but if vaccination does reduce infectiousness then the synergy between antiviral and vaccine might be less than reported here. We assumed no further vaccination during the period of our analysis, which corresponds to the current situation in most middle- and high-vaccinated countries, where the adult population who are willing to be vaccinated and boosted has mostly been immunized. Nevertheless, vaccines for children younger than 5 years old are currently in clinical trials, and are expected to be available during summer 2022. Vaccinating this age group will increase the overall proportion of the population vaccinated in high-income countries, hence reducing further on-going transmission. In addition, WHO recently released an ambitious plan to vaccinate 70% of the global population by mid 2022 [29]. For simplicity, we did not include waning immunity, but explored scenarios where a highler or lower proportions of the population are currently immune. While it is extremely hard to predict when and how different populations will become vaccinated and the percentage of pre-existing immunity in each population, we believe that the scenarios that we have considered here, with different combinations of vaccination coverage and proportions of the population with pre-existing infection-induced immunity, are sufficient to capture the potential population-level effect of the potential use of antiviral treatment on transmission and COVID-19-related severe outcomes.
Taken together, our results suggest that antiviral treatment can be an extremely useful tool to reduce COVID-19 related deaths and to alleviate the COVID-19 burden on overwhelmed and exhausted healthcare systems. Further, our results show that the population level impact of antiviral treatments are enhanced by their synergistic use in combination with vaccination.
Methods
Main model
We adapted a previously developed agent-based transmission model, COVASIM [14]. Namely, we extended the model to include the use of antiviral treatments in the population. We briefly describe here the main features of this model (as were used in the present work) and we describe in full detail the adaptations we made. We refer to the original article by Kerr et al. for the full details of the model implementation. This is an agent-based model that simulates SARS-CoV-2 transmission and interventions that was originally calibrated to data for Seattle, WA, US. Each individual in the population is modeled as an agent in a network, with 500,000 agents. The model has four possible contact layers: “home”, “school”, “work” or “community”. For each location, the model uses country-specific demographics and household sizes. The model simulates infections, interventions antiviral treatment and vaccination. We did 100 runs for each scenario. We report the median of these runs and the 10th and 90th percentiles for the lower and upper bounds.
Disease dynamics
Individuals in the network can be susceptible, exposed (infected but not infectious), asymptomatic, presymptomatic or symptomatic. Symptomatic individuals have one of three fates: they develop a mild, severe, or critical disease. All mild infected individuals recover, and infected individuals reaching a critical state can recover or die. The latent period is sampled from a log-normal distributions with a mean of 4.6 days. The length of time between becoming infectious and developing symptoms is sampled from a log-normal distribution with mean 1.1 days. The times to develop severe symptoms, to progress to become critically ill and to death are sampled from lognormal distributions with means 6.6, 1.5 and 10.7 days respectively. Asymptomatic and mild cases recover on average on 8 days (sampled from log-normal distributions), while severe and critically ill cases recover on average on 18.1 days. Infectiousness was assumed to be linearly correlated to viral load [30]. As with the original model [14], we modeled viral load having two modes: first, a high mode where viral load is highest, around or before symptom onset, and a low mode having a longer duration but a lower viral load (50% lower than during the high mode). We assumed that children are less likely to develop symptoms than adults but equally likely to become infected [31, 32], Table S3. Further, asymptomatic individuals are assumed to be 30% less infectious than symptomatic individuals [33].
Viral transmission
We assumed that antiviral treatment would be deployed under epidemic waves with characteristics similar to the delta or the omicron waves. For the delta wave, we assumed an increased transmissibility of 97% with respect to the ancestral variant [4]. For the omicron wave, we assumed an increased transmissibility of 66% with respect to the delta variant [34]. These assumptions resulted in basic reproduction number R0 for a Delta-like variant of approximately 4.95-5.68 and of 7.75-8.63 for an omicron-like variant. We assumed that 0.5% of the population is infected at the beginning of the simulation and that 40% of the population (selected at random) in each location was previously infected and is immune for the duration of the study (20 and 60% were also explored in Sensitivity analysis). Table S1 summarizes all the parameters used in the model.
Antiviral treatment
As described above, we assumed that antiviral treatment had two main effects. First, antiviral treatment reduced the probability of symptomatic infected individuals becoming hospitalized. Second, we reduced viral load in treated individuals by 25, 50, 75 or 100%. Because in our model, viral load linearly correlates with infectiousness, this assumption resulted in equal reductions in the overall transmissibility of treated infected individuals. Each day of the simulation, the model identified eligible individuals (symptomatic adults over 18, 30, 50 or 65 years of age) whose symptoms onset was within the time frame studied (5 days, 2 days or within 3-5 days of symptoms onset). It then randomly selected a percentage of these individuals (excluding those who are already in treatment) for antiviral treatment initiation. Once an infected symptomatic individual initiated treatment, we assumed that antiviral effects would last for the subsequent days of his or her infection.
Vaccination
We considered a leaky vaccine (that is, a vaccine that confers partial protection to all vaccinated individuals) [35] having three effects on vaccinated individuals: to reduce the probability of acquiring a SARS-CoV-2 infection (denoted byVESUS), to reduce the probability of developing COVID-19 symptoms after infection (denoted by VESYMP), and to reduce the probability of hospitalization conditioned on symptomatic infection (denoted by VEH). Then it follows that
where VEDIS and VEHOSP are the unconditional vaccine effectiveness on symptomatic infection and hospitalization respectively. Importantly, under this model vaccination does not influence onward transmission of infection except through a reduction in SARS-CoV-2 infection. Vaccine effectiveness estimates for each modeled wave can be found in table S2. For Mexico, the US and Belgium, we used country- and age-specific vaccination rates (vaccination rates as of October 12th 2021 and January 3rd, 2022 for the delta and omicron waves respectively), tables S5-S7 [36–39]. For the US and Belgium, vaccinated individuals were further split between those having high (boosted individuals) or low protection for the omicron wave. Since Kenya had fully vaccinated only a very small proportion of its population –1.5% and 7.7% of its population for each of these dates– and targeted front-line health workers, teachers, police and military, we distributed its vaccine randomly among all adults [2, 39, 40].
Data Availability
Code will be available upon publication.
Funding
Scientific Computing Infrastructure at Fred Hutch was funded by ORIP grant S10OD028685. HJ was supported by R56AI143418 and R01CA152089 from the National Institutes of Health. LM, DD and HJ were supported by UM1 AI068635 from the National Institutes of Health. ERB was supported by the Infectious Diseases Clinical Research Consortium through the National Institute of Allergy and Infectious Diseases, part of the National Institutes of Health, under award number UM1AI148684. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health
Author contributions
LM, ERB and HJ conceived the study. LM conducted the analysis and wrote the first draft of the manuscript. HJ, DD, and ERB contributed to writing the manuscript.
Competing interests
The authors declare no competing interests.
Code availability
Code will be available at: https://github.com/lulelita
Supplemental Information
Parameters used in the model.
Vaccine effectiveness values used in the model during the delta and the omicron waves. For the delta wave, vaccine effectiveness assumes full coverage. For the omicron wave, we assumed that boosted vaccinated individuals will get the same protection as that given by full coverage during the delta wave.
Age-specific parameters for disease progression.
Demographic distribution of the adult population
Cumulative deaths over next 6 months in A) Kenya, B) Mexico, C) United States and D) Belgium. Here, we assumed an epidemic wave with parameters similar to those of the omicron epidemic wave (transmissibility, vaccine effectiveness, and vaccination coverage). For each location, the colors represent four possible values of AVT (25, 50, 75 or 100% reduction in viral transmission in treated symptomatic individuals) and a daily treatment initiation rate (DTIR) of 2-100% of adult symptomatic individuals within the first 5 days of symptoms. Gray bars represent baseline cumulative deaths in absence of antiviral treatment.
Cumulative infections over next 6 months for A) Kenya, B) Mexico, C) United States and D) Belgium. Here, we assumed an epidemic wave with parameters similar to those of the omicron epidemic wave (transmissibility, vaccine effectiveness, and vaccination coverage). For each location, the colors represent four possible values of AVT (25, 50, 75 or 100% reduction in viral transmission in treated symptomatic individuals) and a daily treatment initiation rate (DTIR) of 2-100% of adult symptomatic individuals within the first 5 days of symptoms. Gray bars represent baseline cumulative infections in absence of antiviral treatment.
Percentage of deaths averted (compared to a baseline of no antiviral treatment) for A) Kenya, B) Mexico, C) United States and D) Belgium. Here, we assumed an epidemic wave with parameters similar to those of the delta epidemic wave (transmissibility, vaccine effectiveness, and vaccination coverage). For each location, the colors represent four possible values of AVT (25, 50, 75 or 100% reduction in viral transmission in treated symptomatic individuals) and a daily treatment initiation rate (DTIR) of 2-100% of adult symptomatic individuals within the first 5 days of symptoms.
Percentage of infections averted (compared to a baseline of no antiviral treatment) for A) Kenya, B) Mexico, C) United States and D) Belgium. Here, we assumed an epidemic wave with parameters similar to those of the delta epidemic wave (transmissibility, vaccine effectiveness, and vaccination coverage). For each location, the colors represent four possible values of AVT (25, 50, 75 or 100% reduction in viral transmission in treated symptomatic individuals) and a daily treatment initiation rate (DTIR) of 2-100% of adult symptomatic individuals within the first 5 days of symptoms.
Percentage of deaths averted for A) Kenya, B) Mexico, C) United States and D) Belgium assuming 20 (left), 40 (middle) or 60% (right) of the population has been previously infected and is currently immune. For each location, the colors represent four possible values of AVT (25, 50, 75 or 100% reduction in viral transmission in treated symptomatic individuals) and a daily treatment initiation rate (DTIR) of 2-100% of adult symptomatic individuals within the first 5 days of symptoms.
Daily new infections assuming no antiviral treatment (Baseline) or assuming a daily treatment initiation rate of of 10-100% of eligible symptomatic individuals in A) Kenya, B) Mexico, C) United States and D) Belgium assuming 20% of the population has been previously infected and is now recovered. For each location, each column represents a different value of AVT (25, 50, 75 or 100% reduction in viral transmission in treated symptomatic individuals).
Daily new infections assuming no antiviral treatment (Baseline) or assuming a daily treatment initiation rate of of 10-100% of eligible symptomatic individuals in A) Kenya, B) Mexico, C) United States and D) Belgium assuming 60% of the population has been previously infected and is now recovered. For each location, each column represents a different value of AVT 25, 50, 75 or 100% reduction in viral transmission in treated symptomatic individuals).
Percentage of deaths averted (compared to a baseline of no antiviral treatment) for A) Kenya, B) Mexico, C) United States and D) Belgium assuming antiviral treatment would reduce hospitalizations by 30%. Here, we assumed an epidemic wave with parameters similar to those of the delta epidemic wave (transmissibility, vaccine effectiveness, and vaccination coverage). For each location, the colors represent four possible values of AVT (25, 50, 75 or 100% reduction in viral transmission in treated symptomatic individuals) and a daily treatment initiation rate (DTIR) of 2-100% of adult symptomatic individuals within the first 5 days of symptoms.
Percentage of deaths averted (compared to a baseline of no antiviral treatment) for A) Kenya, B) Mexico, C) United States and D) Belgium assuming asymptomatic infections are 50% less infectious. For each location, the colors represent four possible values of AVT (25, 50, 75 or 100% reduction in viral transmission in treated symptomatic individuals) and a daily treatment initiation rate (DTIR) of 2-100% of adult symptomatic individuals within the first 5 days of symptoms.
Percentage of infections averted (compared to a baseline of no antiviral treatment) for A) Kenya, B) Mexico, C) United States and D) Belgium assuming asymptomatic infections are 50% less infectious. For each location, the colors represent four possible values of AVT (25, 50, 75 or 100% reduction in viral transmission in treated symptomatic individuals) and a daily treatment initiation rate (DTIR) of 2-100% of adult symptomatic individuals within the first 5 days of symptoms.
Footnotes
New runs with lower coverage and under a viral transmission similar to the Omicron variant were added.
References
- [1].↵
- [2].↵
- [3].↵
- [4].↵
- [5].↵
- [6].
- [7].
- [8].↵
- [9].↵
- [10].↵
- [11].↵
- [12].↵
- [13].↵
- [14].↵
- [15].↵
- [16].↵
- [17].↵
- [18].↵
- [19].↵
- [20].↵
- [21].↵
- [22].↵
- [23].
- [24].
- [25].
- [26].↵
- [27].↵
- [28].↵
- [29].↵
- [30].↵
- [31].↵
- [32].↵
- [33].↵
- [34].↵
- [35].↵
- [36].↵
- [37].↵
- [38].↵
- [39].↵
- [40].↵
- [41].
- [42].
- [43].
- [44].
- [45].
- [46].
- [47].
- [48].
- [49].
- [50].
- [51].
- [52].
- [53].
- [54].
- [55].
- [56].
- [57].