key: cord-0994591-jyhtp18g authors: Oduro, Bismark; Magagula, Vusi Mpendulo title: COVID-19 intervention models: An initial aggressive treatment strategy for controlling the infection date: 2021-01-23 journal: Infect Dis Model DOI: 10.1016/j.idm.2021.01.007 sha: 87c7d24f75140ba988da49b044f21dd467dd090f doc_id: 994591 cord_uid: jyhtp18g The novel coronavirus (COVID-19) outbreak emerged in December 2019. The disease has caused loss of many lives and has become an unprecedented threat to public health worldwide. We develop simple COVID-19 epidemic models to study treatment strategies to control the pandemic. The results show that eradication of the disease is possible if the efficacy of treatment is perfect. We also investigate the existence of a dual-rate effect. Conditions under which the effect occurs are derived. When the effect is present, a tactic to control the infection might be to initially treat infected individuals aggressively at a relatively high rate to drive the prevalence to a lower region that can be maintained in the long run at relatively moderate rate and cost. Coronavirus is a pathogen that targets the human respiratory system. It is mostly refereed to as COVID-19. The World Health Organization (WHO) officially named it the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The WHO declared the outbreak of SARS-CoV-2 a Public Health Emergency of International Concern on 30 January 2020 [1] . Compared with the SARS-CoV that caused an outbreak of SARS in 2003, SARS-CoV-2 has a stronger transmission capacity. It belongs to the same family of previous outbreaks of coronaviruses (CoVs), namely, severe acute respiratory syndrome (SARS)-CoV and the Middle East respiratory syndrome (MERS)-CoV. These are characterized as agents that are a great public health threat. Patients were admitted to hospitals with an initial diagnosis of pneumonia of an unknown etiology in late December 2019 Wuhan, Hubei Province, China. These patients were epidemiologically linked to a seafood and wet animal wholesale market in China, in the Hubei Province [2, 3] . This seafood and wet animal wholesale market sold live animals such as poultry, bats, marmots, and snakes. The coronavirus has caused a lot of damage economically, health wise and otherwise. Service-oriented economies will be particularly negatively affected, and have more jobs at risk. Countries that are more reliant on tourism (more than 15% of GDP) will be more affected by this crisis. This current crisis is generating spillover effects throughout supply chains. Therefore, countries highly dependent on foreign trade are more negatively affected. Research study suggests that on average, each additional month of crisis costs 2.5-3% [4] of global GDP. Since the beginning of the virus, we have about seven million cases and approximately four hundred thousands deaths [5] with the United States having the highest cases as of June 2, 2020. Almost 50% of the cases have recovered worldwide. To contain the spread of the coronavirus, the WHO suggested that people should wash their hands with soap for more than two minutes, maintain social distancing, wear masks in public spaces, effective mass testing, lockdowns towns where the infection rates are higher. These lockdowns were meant to prepare governments to buy adequate health equipments in order to minimize casualties in terms of death. Many countries embarked on effective testing, had total lockdowns for almost two months in order to minimize the spread of the coronavirus. Some managed to fumigate public places, effectively tested large percentages of their populations. Countries like South Korea, Taiwan, New Zealand managed to bring down a lot of their cases in a short space of time due to their effective mass testing strategies. However, despite all these efforts, the coronaviruses cases spike up once the economy is re-opened. The best solution in containing the virus is getting some treatment drugs. Several studies have been done on medications for treating COVID-19 patients. Cai et al. [6] examined the effects of favipiravir (FPV) versus lopinavir (LPV)/ritonavir (RTV) for the treatment of COVID-19. They discovered that favipiravir showed better therapeutic responses on COVID-19 in terms of disease progression and viral clearance and hence one of the potential drugs for treatment. Many researchers [7, 8] also worked on the remdesivir drug. Remdesivir has been recently a strong candidate for the treatment of COVID-19 and exhibits broad-spectrum antiviral activity against RNA viruses. Zhai et al. [7] worked analyzing the different treatment drugs that have been tried in different clinical trials. Felsenstein et al. [8] also analyzed different treatment drugs and their efficacy on reducing fatalities. They observed that the following anti-viral treatment drugs are potential treatment drugs for COVID-19. One of them is chloroquine which has antiviral activity against Influenza, Chikungunya virus, seasonal CoVs, and SARS. Another anti-viral treatment drug which is known as azithromycin is sometimes used for its immunomodulatory properties, especially in patients with chronic pulmonary disorders [9, 10] . Felsenstein et al. [8] reported that convalescent plasma, i.e. plasma from individuals following COVID-19 resolution and rich in immunoglobulins directed against SARS-CoV2, is being entertained as possible treatment option. However, many clinical trials are being undertaken to determine the best treatment drug. These medications once developed will be useful in mitigating the spread of the coronavirus. However, the question is, what should be the optimal way of administering the medications so that the virus is contained? In this paper, we try to answer this question by using the idea of dual-rate effect analysis. This effect has extensively studied in [11, 12, 13] and has shown to occur several models, including the re-infestation models in [11, 12, 13] and black pod disease model in [14] . In [13] , Oduro et al. stated that "the dual-rate effect has important implications for control efforts, as in the presence of this effect a strategy of adopting a highly aggressive short-term intervention at the beginning of a sustained control program may be more cost-effective in the long run than treatment at a fixed rate without a heightened initial effort". We study compartmental COVID-19 models that incorporate treatment as a control measure with the ultimate goal of assessing optimal treatment strategies and deriving recommendations. These models were inspired by efforts to eradicate or contain the COVID-19 pandemic. The results reported here are purely mathematical and more widely applicable. They may be of particular relevance to the virus control efforts in lower middle-income countries (LMIC) where funds for the COVID-19 control programs are severely limited and often in competition for funds with one another as well as with other pressing health care needs and sectors such as education and infrastructure. The rest of the paper is organized as follows, in Section 2, we formulate and analyze simple SEIR models. Dual-rate effect analysis are done in Sections 3. We then summarize the results in Section 4. In this section we present two mathematical models, the first one with 100% efficacy of treatment and the second one with less than 100% efficacy of treatment. 2.1.1 First Scenario: 100% efficacy of treatment. In this model, we assume that the population can be subdivided into four compartments. These are the susceptible (S) class, exposed (E) class, infected (I) class and lastly the recovered (R) class. The susceptible class recruits new population at the rate b. In all four classes, the natural death rate is given by µ and the mortality due to COVID-19 is given by δ. The transmission rate is given by β and the progression rate of exposed individuals to infected individuals is given by p. The infected individuals are treated with a recommended COVID-19 medication at rate v. These infected individuals can also be given other medical therapy at the rate q. The combination of the treatment and other medical therapy ensures that the infected individuals recover. The dynamics explained in this subsection are displayed in Figure 1 . The dynamics in Figure 1 can be represented as a system of nonlinear ordinary differential equations given by We summarize the existence of equilibrium in the following theorem. . (ii) When R c ≤ 1, then the disease-free equilibrium DF E = b µ , 0, 0, 0 is the only biologically feasible equilibrium. J o u r n a l P r e -p r o o f (iii) When R c > 1, then in addition to the DF E, there exists a unique endemic equilibrium EE = (S * , E * , I * , R * ) with 2.1.2 Second Scenario: Less 100% efficacy of treatment. In this model, we assume that the population can be subdivided into five compartments. These are the susceptible (S) class, exposed (E) class, infected (I 1 ) class, infected (I 2 ) class after receiving ineffective treatment, and lastly the recovered (R) class. The recruitment rate is denoted by b, the natural mortality rate by µ and the mortality rate due to COVID-19 by δ. The transmission rate is given by β and the progression rate of exposed individuals to infected individuals is given by p. The infected individuals are treated at rate v. These infected individuals can also be given other medical therapy at the rate q. However, if the treatment is ineffective, then the individuals from I 1 move to I 2 at the rate (1 − k)v where k is the fraction of individuals who have been successfully treated with recommended COVID-19 medication. Thus those who have been successfully treated progress to the recovered class from the infected class I 1 . The combination of the COVID-19 medication and other medical therapy ensures that the infected individuals recover. Figure 2 shows the dynamics explained in this subsection. The dynamics in Figure 2 can be represented as a system of nonlinear ordinary differ-5 J o u r n a l P r e -p r o o f ential equations of the form: We summarize the existence of equilibrium in the following theorem. Theorem 2.2 Consider the model 2 (i) The control reproduction number (R c ) is given by . (ii) When R c ≤ 1, then the disease-free equilibrium DF E = b µ , 0, 0, 0, 0 is the only biologically feasible equilibrium. (iii) When R c > 1, then in addition to the DF E, there exists a unique endemic equilibrium EE = (S * , E * , I * 1 , I * 2 , R * ) with We use Pennsylvania (USA) state-level time series for cumulative active COVID-19 cases compiled by Worldometer [5] , from March 15, 2020, through May 31, 2020 to estimate the model parameters. The sources of the data provided by Worldometer were checked for confirmation. The model fitting was performed using the standard nonlinear least squares approach using the fmincon Optimization Toolbox embedded in MATLAB. We provide the estimated parameters in Table 1 . The recruitment rate was estimated from the population data for Pennsylvania, by US Census Bureau [15] and the natural death rate was adopted from the National Center for Health Statistics[22]. Since the treatment rate, v, is our man control parameter, let us consider the effect of this parameter on the control reproduction number for both models. Thus, the R c is non-increasing with respect to v. Also This suggests that eradication will be achieved with sufficiently high treatment rate, the control reproduction number decreases significantly as the treatment rate increases. We have The R c decreases with v but asymptotically approaches the right-hand side of 4. This suggests that with some parameter settings, high treatment rate may not be sufficient to eradicate the disease. Figures 3 -5 show the effect of varying the treatment rate, v, on the infectious class in scenario 1 and 2 respectively. The values of parameters used are as given in Table 1 unless otherwise stated. First, increasing the treatment rate decreases the population of infected individuals as shown in Figures 3 and 4 . We also note that decreasing the treatment rate increases the population of infected individuals. The treatment intervention delays the peak of the epidemic curve and thus reduces the number of infested individuals. Second, Figure 5 shows the effect of treatment on the infected population that still remain infected after being treated (Scenario 2). Of course, the potency of the treatment will play a major factor, fixing k = 0.8 and increasing the treatment rate leads to increase in the I 2 (t) infected population as shown in Figure 5 . Other medical therapies may be needed to contain individuals in the I 2 (t) class if the COVID-19 treatment is not 100% effective. In Figure 6 , we show the effect of k (fraction of successfully COVID-19 treatment) on I 2 (t). Increasing k, decreases the I 2 (t) population. For example, increasing the baseline k = 0.5 by 50% reduces the population in I 2 (t) from about 7 × 10 5 to below 4 × 10 5 individuals. Thus, the efficacy of the treatment has direct impact on this infected population. Baseline k 20% reduction in k 40% reduction in k 50% increase in k 80% increase in k 100% increase in k Figure 6 : Impact of k on I 2 (t) in the scenario 2 model, baseline k = 0.5 In this section, we investigate long-term cost effective strategies for controlling the COVID-19 infection. Assuming there is a recommended medication treatment, we show the existence of a dual-rate effect. When this effect is present, an initial aggressive treatment rate may plunge the prevalence to a level that can be controlled at a lower treatment rate and relatively moderate cost [11, 12, 13]. In what follows, we provide some useful definitions with notations that are needed in subsequent section. is the overall cases of COVID-19 between T − and T + . Following the approach in [11, 12, 13] , we conceptualize the cost (C) as the total amount of treatment needed over the interval and it is defined as where v is the COVID-19 treatment rate. We established in Section 2.1 that when v is kept constant, the susceptible, exposed, infested, and recovered individuals approach an equilibrium, with I * (v) denoting the the number of infected individuals at equilibrium. If an endemic equilibrium can be reached, then the long-term cost of COVID-19 treatment can be conceptualized as follows: It is important to note that, the equilibrium value I * (v) is a monotone decreasing function of v. However, under certain conditions the dependence of C on v is not monotone. ii. I * (v 1 ) < I * (v 2 ). We will use the following observation. Proposition 3.1 If the cost function C(v) has a strict local maximum at some v ∈ (0, ∞), then the dual-rate effect occurs. Proof. Let v c be the critical point of C(v). Then there exist two distinct treatment Since I * depends on the treatment rate and decreases monotonically with v, a low treatment rate 0 < v l < v crit will produce a relatively high I * and a high treatment rate v c < v h , will yield a relative low I * . That is, I * (v l ) > I * (v h ). This implies that there is a cost C(v l ) = C(v h ) that corresponds to two different values of I * (v). We present results for dual-rate rate effect. Theorem 3.1 Consider the model with perfect treatment efficacy and let β > µ(p+µ)(q+δ+µ) pb . Then the dual-rate effect occurs whenever the endemic equilibrium exists. Proof. Using the I * (v) of the endemic equilibrium in Theorem 2.1 and the definition in equation (7), we have It follows that This yields a critical treatment rate v * = βpb(q+δ+µ) µ(p+µ) By the second derivative test, the cost function C(v) has a local maximum value at v * . The rest of the proof follows directly from Proposition 3.1. In Figures 7 and 8 , we present an illustration of the presence of the dual-rate effect. For example, (with b = 9213, β = 1.386e − 06, µ = 0.008638, p = 0.2710, q = 0.4655, δ = 0.014067) a cost of 500 corresponds to two different treatment rates, approximately v 1 = 0.05 and v 2 = 0.8. Since the I * is a monotone decreasing function of v, I * (0.05) = 612 individuals, a relatively low infection level than I * (0.8) = 10419 individuals. Here, an optimal strategy is to initially treat at rate v = 0.8 to push the system into a region with low endemic level of infection that can be maintained in the long run at relatively low rate v = 0.05. Then the dual-rate effect occurs whenever the endemic equilibrium exists. Proof. Similarly to the proof of Theorem 3.1, let The critical treatment rate occurs at For example, choosing β = µ(q+δ+µ)(p+µ) for all 0 < k < 1. The cost function has a strict local maximum at v * , the existence of the dual-rate effect follows from Proposition 3.1. As illustrated in Figures 9 and 10 , the dual rate occurs in this scenario where the efficacy of treatment is not perfect. In this paper, we designed and analyzed two mathematical models with treatment, one with 100% efficacy of the treatment and the other one with less than 100% efficacy of the treatment. We derived conditions for both disease-free and endemic equilibria to exist, and argued that the infection can be eradicated in the scenario with perfect efficacy of treatment if the treatment rate is sufficiently high. On the other hand, treatment may not be sufficient for eradication if the efficacy is not prefect. The prevalence can only be driven to arbitrarily low levels by sufficiently large treatment rate. We also investigated the long-time cost effective strategies for controlling the COVID-19 infection. In particular, we showed the existence of the dual-rate effect in both models. When this effect is present, the optimal strategy might be an initially very aggressive treatment to plunge the prevalence to a level that can be controlled at a lower treatment rate and in the long run at relatively moderate cost. The high initial treatment rate may require a substantial initial investment (cost) which may payoffs in the long run. Without the initially very high treatment, the same average cost would only allow a much less significant reduction in long-term prevalence level. COVID-19 and the cardiovascular system Pneumonia of Unknown Etiology in Wuhan, China: Potential for International Spread Via Commercial Air Travel Outbreak of Pneumonia of Unknown Etiology in Wuhan China: the Mystery and the Miracle Economic Effects of Coronavirus Outbreak (COVID-19) on the World Economy The epidemiology, diagnosis and treatment of COVID-19 COVID-19: immunology and treatment options Azithromycin polarizes macrophages to an M2 phenotype via inhibition of the STAT1 and NF-KB signaling pathways Predictive screening of M1 and M2 macrophages reveals the immunomodulatory effectiveness of post spinal cord injury azithromycin treatment Models of disease vector control: When can aggressive initial intervention lower long-term cost A model of insect control with imperfect treatment Assessing the effect of fungicide treatment on Cocoa black pod disease in Ghana: Insight from mathematical modeling Mathematical Tools for Understanding Infectious Disease Dynamics. Kindle Edition Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission A brief history of R0 and a recipe for its calculation Australia's been asked to make a coronavirus vaccine at "unprecedented speed Modeling Infectious diseases in Humans and Animals The first coronavirus drug candidate is set for testing in China Ethical approval This article does not contain any studies with human participants performed by any of the authors.