# A Model for SARS-CoV-2 Infection with Treatment

Authors: Amar Nath Chatterjee; Fahad Al Basir
Journal: Computational and Mathematical Methods in Medicine
 (2020)

---

## Abstract

The current emergence of coronavirus (SARS-CoV-2) puts the world in threat. The structural research on the receptor recognition by SARS-CoV-2 has identified the key interactions between SARS-CoV-2 spike protein and its host (epithelial cell) receptor, also known as angiotensin-converting enzyme 2 (ACE2). It controls both the cross-species and human-to-human transmissions of SARS-CoV-2. In view of this, we propose and analyze a mathematical model for investigating the effect of CTL responses over the viral mutation to control the viral infection when a postinfection immunostimulant drug (pidotimod) is administered at regular intervals. Dynamics of the system with and without impulses have been analyzed using the basic reproduction number. This study shows that the proper dosing interval and drug dose both are important to eradicate the viral infection.


---

## Body

## 1. Introduction

A novel coronavirus named SARS-CoV-2 (an interim name proposed by WHO (World Health Organization)) became a pandemic since December 2019. The first infectious respiratory syndrome was recognized in Wuhan, Hubei province of China. Dedicated virologists identified and recognized the virus within a short time [1]. The SARS-CoV-2 is a single-stranded RNA virus genome which is closely related to severe acute respiratory syndrome- (SARS-) CoV [2]. The infection of SARS-CoV-2 is associated with a SARS-CoV-like a disease with a fatality rate of 3.4% [3]. The World Health Organization (WHO) have named the disease as COVID-19 and declared it as a public health emergency worldwide [4].The common symptoms of COVID-19 are fever, fatigue, dry cough, and myalgia. Also, some patients suffer from headaches, abdominal pain, diarrhea, nausea, and vomiting. In the acute phase of infection, the disease may lead to respiratory failure which leads to death also. From clinical observation, within 1-2 days after patient symptoms, the patient becomes morbid after 4-6 days and the infection may clear within 18 days [5] depending on the immune system. Thus, appropriate quarantine measure for a minimum of two weeks is taken by the public health authorities for inhibiting community spread [6].In [1], Zhou et al. identified the respiratory tract as the principal infection site for COVID-19 infection. SARS-CoV-2 infects primary human airway epithelial cells. The angiotensin-converting enzyme 2 (ACE2) receptor of epithelial cells plays an important role in cellular entry [1, 7]. It has been observed that ACE2 could be expressed in the oral cavity. ACE2 receptors are higher in the tongue than buccal and gingival tissues. These findings imply that the mucosa of the oral cavity may be a potentially high-risk route of COVID-19 infection. Thus, epithelial cells of the tongue are the major routes of entry for COVID-19. Zhou et al. [1] also reported that SARS-CoV-2 spikes (S) bind with the ACE2 receptor of epithelial cells with high affinity. The bonding between S (spikes) of SARS-CoV-2 and ACE2 [7] results from the fusion between the viral envelope and the target cell membrane, and the epithelial cells become infected. The S protein plays a major role in the induction of protective immunity during the infection of SARS-CoV-2 by eliciting neutralization antibody and T cell responses [8]. The S protein is not only capable of neutralizing antibody, but it also contains several immunogenic T cell epitopes. Some of the epitopes are found in either the S1 or S2 domain. These proteins are useful for SARS-CoV-2 drug development [9].We know that virus clearance after acute infection is associated with strong antibody responses. Antibody responses have the potential to control the infection [10]. Also, CTL responses help to resolve infection and virus persistence caused by weak CTL responses [11]. Antibody responses against SARS-CoV-2 play an important role in preventing the viral entry process [8]. Hsueh et al. [2] found that antibodies block viral entry by binding to the S glycoprotein of SARS-CoV-2. To fight against the pathogen SARS-CoV-2, the body requires SARS-CoV-2-specific CD4+ T helper cells for developing this specific antibody [8]. Antibody-mediated immunity protection helps the anti-SARS-CoV serum to neutralize COVID-19 infection. Besides that, the role of T cell responses in COVID-19 infection is very much important. Cytotoxic T lymphocyte (CTL) responses are important for recognizing and killing infected cells, particularly in the lungs [8]. But the kinetics of the CTL responses and antibody responses during SARS-CoV-2 infection is yet to be explored. Our study will focus on the role of CTL and its possible implication on treatment and drug development. The drug that stimulates the CTL responses represents the best hope for control of COVID-19. Here, we have modeled the situation where CTLs can effectively control the viral infection when the postinfection drug is administered at regular intervals.Mathematical modeling with real data can help in predicting the dynamics and control of an infectious disease [12, 13]. A four-dimensional dynamical model for a viral infection is proposed by Tang et al. [14] for MERS-CoV mediated by DPP4 receptors. In the case of SARS-CoV-2, the infection process is almost similar with MERS-CoV and SARS-CoV. For SARS-CoV-2 infection, the ACE2 receptors of epithelial cells are the major target area.Since the dynamics of the disease transmission of SARS-CoV-2 in the cellular level is yet to be explored, we investigate the system in the light of the previous literature of [14–18] to formulate the dynamic model which plays a significant role in describing the interaction between uninfected cells, free virus, and CTL responses. We propose a novel deterministic model which describes the cell biological infection of SARS-CoV-2 with epithelial cells and the role of the ACE2 receptor.We explained the dynamics in the acute infection stage. It has been observed that CTLs proliferate and differentiate antibody production after they encounter antigens. Here, we investigate the effect of CTL responses over the viral mutation to control viral infection when a postinfection drug is administered at regular intervals by a mathematical perspective.It is clinically evident that immunostimulants play a crucial role in the case of respiratory disease. Among the currently available immunostimulants, pidotimod is the most effective for the respiratory disease [19]. Pidotimod increases the level of immunoglobulins (IgA, IgM, and IgG) and activates the CTL responses to fight against the disease.In this article, we have considered the infection dynamics of SARS-CoV-2 infection in the acute stage. We have used impulsive differential equations to study the immunostimulant drug dynamics and the effects of perfect drug adherence. In recent years, the effects of perfect adherence have been studied by using impulsive differential equations in [20–26]. With the help of impulsive differential equations, the effect of maximal acceptable drug holidays and optimal dosage can be found more precisely [20, 26].The article is organized as follows. The very next section contains the formulation of the impulsive mathematical model. Dynamics of the system without impulses has been provided in Section3. The system with impulses has been analyzed in Section 4. Numerical simulations, on the basis of the outcomes of Sections 3 and 4, have been included in Section 5. Discussion in Section 6 concludes the paper.


## 2. Model Formulation

As discussed in the previous section, we propose a model considering the interaction between epithelial cells and SARS-CoV-2 virus along with lytic CTL responses over the infected cells. We consider five populations, namely, the uninfected epithelial cellsTt, infected cells It, ACE2 receptor of the epithelial cells Et, SARS-CoV-2 virus V(t) and CTLs against the pathogen Ct.In this model, we consider which represents the concentration of ACE2 on the surface of uninfected cells, which can be recognized by the surface spike (S) protein of SARS-CoV-2 [27].It is assumed that the susceptible cells are produced at a rateλ1 from the precursor cells and die at a rate dT. The susceptible cells become infected at a rate βEtVtTt. The constant dI is the death rate of the infected cells. Infected cells are also cleared by the body’s defensive CTLs at a rate p.The infected cells produce new viruses at the ratemdI during their life, and dV is the death rate of new virions, where m is any positive integer. It is also assumed that ACE2 is produced from the surface of uninfected cells at the constant rate λ2 and the ACE2 is destroyed, when free viruses try to infect uninfected cells, at the rate θβEtVtTt and is hydrolyzed at the rate dEE.CTL proliferation in the presence of infected cells is described by the term(1)αIC1−CCmax,which shows the antigen-dependent proliferation. Here, we consider the logistic growth of CTL with Cmax as the maximum concentration of CTL, and dc is its rate of decay.With the above assumptions, we have the following mathematical model characterizing the SARS-CoV-2 dynamics:(2)dTdt=λ1−βEVT−dTT,dIdt=βEVT−dII−pIC,dVdt=mdII−dVV,dEdt=λ2−θβEVT−dEE,dCdt=αIC1−CCmax−dcC.A short description of the model parameters and their values is shown in Table1. We now modify the above model by incorporating pulse periodic drug dosing using impulsive differential equations [28, 29].Table 1
Set of parameter values used of numerical simulations.
ParameterExplanationAssigned valueλ1Production rate of uninfected cell5λ2Production rate of ACE21βDisease transmission rate0.0001θBonding rate of ACE20.3dTDeath rate of uninfected cells0.1dIDeath rate of infected cells0.1dVRemoval rate of virus0.1dEHydrolyzing rate of epithelial cells0.1dcDecay rate of CTL0.1pKilling rate of infected cells by CTL0.01mNumber of new virions produced10-100αProliferation rate of CTL0.22CmaxMaximum proliferation of CTL100We consider the perfect adherence behavior of the immunostimulant drug for SARS-CoV-2-infected patients at fixed drug dosing timestk, k∈ℕ.We assume that CTL cells increase by a fixed amountω, which is proportional to the total number of CTLs that the drug can stimulate. Thus, the above model takes the following form:
(3)dTdt=λ1−βEVT−dTT,dIdt=βEVT−dII−pIC,dVdt=mdII−dVV,dEdt=λ2−θβEVT−dEE,dCdt=αIC1−CCmax−dcC,t≠tk,Ctk+=ω+Ctk−,t=tk.Here,Ctk− denotes the CTL cell concentration immediately before the impulse, Ctk+ denotes the concentration after the impulse, and ω is the fixed amount which is proportional to the total number of CTLs the drug stimulates at each impulse time tk, k∈ℕ.Remark 1.
It can be noted that when there is no drug application in the system, model (3) becomes model (2).


## 3. Analysis of the System without the Drug

In this section, we analyze the dynamics of the system without impulses, i.e., system (1). We have derived the basic reproduction number for the system. Stability of equilibria is discussed using the number.


### 3.1. Existence of Equilibria

Model (2) has three steady states, namely, (i) the disease-free equilibrium E1λ1/dT,0,0,λ2/dE,0; (ii) with E¯>dTdV/βλ1m, there is a CTL response-free equilibrium, E2T¯,I¯,V¯,E¯,0, where
(4)T¯=dVβmE¯,I¯=βλ1mE¯−dTdVβdImE¯,V¯=βλ1mE¯−dTdVβdVE¯,E¯=−θβλ1m−βλ2m+θβλ1m−βλ2m2+4βmdTdVdEθ2βdEm;and (iii) the endemic equilibrium E∗ which is given by
(5)T∗=λ1α−dIαI∗−pdcdTα,V∗=dImI∗dV,E∗=λ2α−θαdII∗−θpdcdEα,(6)C∗=αI∗−dcCmaxαI∗,where I∗ is the positive root of the cubic equation
(7)L0I3+L1I2+L2I+L3=0,with
(8)L0=−α2θβdI3m,L1=−2αθβdI2dcmp+α2θβdI2λ1m+α2βdI2λ2m,L2=α2dTdIdVdE+αθβdIdcλ1mp+αβdIdcλ2mp−α2βdIλ1λ2m−θβdIdc2mp2,L3=αdTdVdEdcp.Remark 2.
Note thatL0<0 and L3>0. Thus, equation (7) has at least one positive real root. If L1>0 and L2<0, then (3) can have two positive roots. For a feasible endemic equilibrium, we also need
(9)minλ1α−pdcdIα,λ2α−θpdcθαdI>I∗>dcα.


### 3.2. Stability of Equilibria

In this section, the characteristic equation at any equilibrium is determined for the local stability of system (2). Linearizing system (2) at any equilibrium ET,I,V,E,C yields the characteristic equation
(10)Δξ=∣ξIn−A∣=0,where In is the identity matrix and A=aij is the following 5×5 matrix given by
(11)A=−βEV−dT0−βET−βVT0βEV−dI−pCβETβVT−pI0dIm−dv00−θβEV0−θβET−θβVT−dE00αC1−CCmax00a55,with a55=αI1−2C/Cmax−dc. We finally get the characteristic equation as
(12)ψξ=ξ5+A1ξ4+A2ξ3+A3ξ2+A4ξ+A5=0.The coefficientsAi, i=1,2,⋯,5, are given in the appendix.Looking at stability of any equilibriumE, the Routh-Hurwitz criterion gives that all roots of this characteristic equation (12) have negative real parts, provided the following conditions hold
(13)A5>0,A1A2−A3>0,A3A1A2−A3−A1A1A4−A5>0,A1A2−A3A3A4−A2A5−A1A4−A52>0.Let us define the basic reproduction number as(14)R0=mβλ1λ2dTdEdV.Then, using (5), we can derived the following result.Theorem 3.
Disease-free equilibriumE1λ1/dT,0,0,λ2/dE,0 of model (2) is stable for R0<1 and unstable for R0>1.AtE2, one eigenvalue is −dc and the rest of the eigenvalues satisfy the following equation:
(15)ξ4+B1ξ3+B2ξ2+B3ξ+B4=0.The coefficientsBi, i=1,2,⋯,5, are given in the appendix.Using the Routh-Hurwitz criterion, we have the following theorem:Theorem 4.
The CTL-free equilibrium,E2T¯,I¯,V¯,E¯,0, is asymptotically stable if and only if the following conditions are satisfied:
(16)B1>0,B2>0,B3>0,B4>0,B1B2−B3>0,B1B2−B3B3−B12B4>0.DenotingAi∗=AiE∗ and using (5), we have the following theorem establishing the stability of coexisting equilibrium E∗.Theorem 5.
The coexisting equilibriumE∗ is asymptotically stable if and only if the following conditions are satisfied:
(17)A5∗>0,A1∗A2∗−A3∗>0,A3∗A1∗A2∗−A3∗−A1∗A1∗A4∗−A5∗>0,A1∗A2∗−A3∗A3∗A4∗−A2∗A5∗−A1∗A4∗−A5∗2>0.


## 3.1. Existence of Equilibria

Model (2) has three steady states, namely, (i) the disease-free equilibrium E1λ1/dT,0,0,λ2/dE,0; (ii) with E¯>dTdV/βλ1m, there is a CTL response-free equilibrium, E2T¯,I¯,V¯,E¯,0, where
(4)T¯=dVβmE¯,I¯=βλ1mE¯−dTdVβdImE¯,V¯=βλ1mE¯−dTdVβdVE¯,E¯=−θβλ1m−βλ2m+θβλ1m−βλ2m2+4βmdTdVdEθ2βdEm;and (iii) the endemic equilibrium E∗ which is given by
(5)T∗=λ1α−dIαI∗−pdcdTα,V∗=dImI∗dV,E∗=λ2α−θαdII∗−θpdcdEα,(6)C∗=αI∗−dcCmaxαI∗,where I∗ is the positive root of the cubic equation
(7)L0I3+L1I2+L2I+L3=0,with
(8)L0=−α2θβdI3m,L1=−2αθβdI2dcmp+α2θβdI2λ1m+α2βdI2λ2m,L2=α2dTdIdVdE+αθβdIdcλ1mp+αβdIdcλ2mp−α2βdIλ1λ2m−θβdIdc2mp2,L3=αdTdVdEdcp.Remark 2.
Note thatL0<0 and L3>0. Thus, equation (7) has at least one positive real root. If L1>0 and L2<0, then (3) can have two positive roots. For a feasible endemic equilibrium, we also need
(9)minλ1α−pdcdIα,λ2α−θpdcθαdI>I∗>dcα.


## 3.2. Stability of Equilibria

In this section, the characteristic equation at any equilibrium is determined for the local stability of system (2). Linearizing system (2) at any equilibrium ET,I,V,E,C yields the characteristic equation
(10)Δξ=∣ξIn−A∣=0,where In is the identity matrix and A=aij is the following 5×5 matrix given by
(11)A=−βEV−dT0−βET−βVT0βEV−dI−pCβETβVT−pI0dIm−dv00−θβEV0−θβET−θβVT−dE00αC1−CCmax00a55,with a55=αI1−2C/Cmax−dc. We finally get the characteristic equation as
(12)ψξ=ξ5+A1ξ4+A2ξ3+A3ξ2+A4ξ+A5=0.The coefficientsAi, i=1,2,⋯,5, are given in the appendix.Looking at stability of any equilibriumE, the Routh-Hurwitz criterion gives that all roots of this characteristic equation (12) have negative real parts, provided the following conditions hold
(13)A5>0,A1A2−A3>0,A3A1A2−A3−A1A1A4−A5>0,A1A2−A3A3A4−A2A5−A1A4−A52>0.Let us define the basic reproduction number as(14)R0=mβλ1λ2dTdEdV.Then, using (5), we can derived the following result.Theorem 3.
Disease-free equilibriumE1λ1/dT,0,0,λ2/dE,0 of model (2) is stable for R0<1 and unstable for R0>1.AtE2, one eigenvalue is −dc and the rest of the eigenvalues satisfy the following equation:
(15)ξ4+B1ξ3+B2ξ2+B3ξ+B4=0.The coefficientsBi, i=1,2,⋯,5, are given in the appendix.Using the Routh-Hurwitz criterion, we have the following theorem:Theorem 4.
The CTL-free equilibrium,E2T¯,I¯,V¯,E¯,0, is asymptotically stable if and only if the following conditions are satisfied:
(16)B1>0,B2>0,B3>0,B4>0,B1B2−B3>0,B1B2−B3B3−B12B4>0.DenotingAi∗=AiE∗ and using (5), we have the following theorem establishing the stability of coexisting equilibrium E∗.Theorem 5.
The coexisting equilibriumE∗ is asymptotically stable if and only if the following conditions are satisfied:
(17)A5∗>0,A1∗A2∗−A3∗>0,A3∗A1∗A2∗−A3∗−A1∗A1∗A4∗−A5∗>0,A1∗A2∗−A3∗A3∗A4∗−A2∗A5∗−A1∗A4∗−A5∗2>0.


## 4. Dynamics of the System with Impulsive Drug Dosing

In this section, we consider the model system (3). Before analyzing the system, we first discuss the one-dimensional impulse system as follows:
(18)dCdt=αIC1−CCmax−dcC,t≠tk,Ctk+=ω+Ctk−,t=tk.Ctk− denotes the CTL responses immediately before the impulse drug dosing, Ctk+ denotes the concentration after the impulse, and ω is the dose that is taken at each impulse time tk, k∈ℕ.We now consider the following linear system:(19)dCdt=−dcC,t≠tk,ΔC=ω,t=tk,where Δ=Ctk+−Ctk−. Let τ=tk+1−tk be the period of the campaign. The solution of system (10) is
(20)Ct=Ctk+e−dct−tk,fortk<t≤tk+1.In presence of impulsive dosing, we can get the recursion relation at the moments of impulse as(21)Ctk+=Ctk−+ω.Thus, the amount of CTL before and after the impulse is obtained as(22)Ctk+=ω1−e−kτdc1−e−τdc,Ctk+1−=ω1−e−kτdce−τdc1−e−τdc.Thus, the limiting case of the CTL amount before and after one cycle is as follows:(23)limk→∞Ctk+=ω1−e−τdc,limk→∞Ctk+1−=ωe−τdc1−e−τdc,Ctk+1+=ωe−τdc1−e−τdc+ω=ω1−e−τdc.Definition 6.
LetΛ≡Su,Sa,I,C and B0=B:R+4→R+; then, we say that B belong to class B0 if the following conditions hold:
(i)
B is continuous on tk,tk+1×R+3, n∈N, and for all Λ∈R4, limt,μ→tk+,ΛBt,μ=Btk+,Λ exists(ii)
B is locally Lipschitzian in ΛWe now recall some results for our analysis from [28, 29].Lemma 7.
LetZt be a solution of system (9) with Z0+≥0. Then, Zit≥0, i=1,⋯,4, for all t≥0. Moreover, Zit>0, i=1,⋯,4, for all t>0 if Zi0+>0, i=1,⋯,4.Lemma 8.
There exists a constantγ such that Tt≤γ, It≤γ, Vt≤γEt≤γ, and Ct≤γ for each and every solution Zt of system (9) for all sufficiently large t.Lemma 9.
LetB∈B0 and also consider that
(24)D+Bt,Z≤jt,Bt,Zt,t≠tk,Bt,Zt+≤ΦnBt,Zt,t=tk,where j:R+×R+→R is continuous in tk,tk+1 for e∈R+2, n∈N, the limit limt,V→tk+jt,g=jtk+,x exists, and Φnii=1,2:R+→R+ is nondecreasing. Let yt be a maximal solution of the following impulsive differential equation:
(25)dxtdt=jt,xt,t≠tk,xt+=Φnxt,t=tk,x0+=x0,existing on 0+,∞. Then, B0+,Z0≤x0 implies that Bt,Zt≤yt,t≥0, for any solution Zt of system (9). If j satisfies additional smoothness conditions to ensure the existence and uniqueness of solutions for (12), then yt is the unique solution of (12).We now consider the following subsystem:(26)dCtdt=−dcC,t≠tk,Ctk+=Ctk+ω,C0+=C0.The lemma provided above gives the following result.Lemma 10.
System (13) has a unique positive periodic solution C~t with period τ and given by
(27)C~t=ωexp−dct−tk1−exp−τdc,tk<t≤tk+1,C~0+=dc1−exp−τdc.We use this result to derive the following theorem.Theorem 11.
The disease-free periodic orbitT~,0,0,E~,C~ of system (2) is locally asymptotically stable if
(28)R~0<1,where
(29)R~0=mdIβdTdEdVτ∫0τT~E~dI+pC~dt.Proof.
Let the solution of system (9) without infected people be denoted by T~,0,0,E~,C~, where
(30)C~t=ωexp−dct−tk1−exp−τdc,tk<t≤tk+1,with initial condition C0+ as in Lemma 10. We now test the stability of the equilibria. The variational matrix at T~,0,0,E~,C~ is given by
(31)Mt=mij=−dT0m13000−dI+pC~βE~T~000mdI−dv0000m43−dE00m5200−dc.The monodromy matrixℙ of the variational matrix Mt is
(32)ℙτ=Inexp∫0τMtdt,where In is the identity matrix. Note that m13, m43, and m52 are not required for this analysis; therefore, we have not mentioned their expressions.We can writeℙτ=diagσ1,σ2,σ3,σ4,σ5, where σi, i=1,2,3,4,5, are the Floquet multipliers and they are determined as
(33)σ1=exp−dTτ,σ2,3=exp∫0τ12−A±A2−4Bdt,σ4=exp−dEτ,σ5=exp−dcτ.Here,A=dI+dV+pC~ and B=dVdI+pC~−mdIβE~T~. Clearly, λ1,4,5<1. It is easy to check that A2−4B>0, and if B≥0 and hold, then we have λ2,3<1. Thus, according to Floquet theory, the periodic solution T~,0,0,E~,C~ of system (9) is locally asymptotically stable if the conditions given in (14) hold.


## 5. Numerical Results and Discussion

In this section, we have observed the dynamical behaviors of the system without the drug (Figures1 and 2) and with impulsive effect of the drug dose (Figures 3 and 4) through numerical simulations taking the parameters mainly from [14, 19, 30].Figure 1
Existence and stability of equilibria is shown with respect toR0. Parameter values used in this figure are taken from Table 1 and m=10. We have varied the value of β in 0.00001,0.0001.
(a)(b)Figure 2
In the absence of the drug, the effect of the growth rate of CTL, i.e.,α on the steady-state values of model population, is shown. Parameter values used in this figure are the same as Figure 1 except α.
(a)(b)(c)(d)(e)Figure 3
Numerical solution of the model system with and without the drug dose is shown taking parameters as in Figure1. In this figure, τ=2 and ω=50.
(a)(b)(c)(d)(e)Figure 4
Numerical solution of the model system for different rates of drug dosing and different intervals of impulses.
(a)(b)(c)(d)(e)We have mainly focused on the role of CTL and its possible implication on the treatment and drug development. The drug that stimulates the CTL responses represents the best hope for control of COVID-19. Here, we have determined the situation where CTLs can effectively control the viral infection when the postinfection drug is administered at regular intervals.Existence of equilibria of the system without the drug dose is shown for different values of basic reproduction numberR0. In plotting Figure 1, we have varied the value of infection rate β. It is observed that for the lower infection rate (that corresponds to R0<1), disease-free equilibrium E1 is stable (corroborated with Theorem 3). It becomes unstable and ensures the existence of the CTL-free equilibrium E2 which is stable if R0<2.957 (which corresponds to β=0.00005963) and unstable otherwise. (This satisfies Theorem 4.) Again, we see that when E2 is unstable, E∗ is feasible. Also, whenever E∗ exists, it is stable which verified the Theorem 5.The effect of the immune response rateα is plotted in Figure 2. We observe that in the absence of the drug, the CTL count and ACE2 increase with increasing value of α. The steady-state value of infected cell I∗ and virus V∗ decreases significantly as α increases.Due to the impulsive nature of the drugs, there are no equilibria of the system; i.e., population does not reach towards the equilibrium point, rather approach a periodic orbit. Hence, we evaluate equilibrium-like periodic orbits. There are two periodic orbits of system (3), namely, the disease-free periodic orbit and endemic periodic orbit. Here, our aim is to find the stability of the disease-free periodic orbit.Figure3 compares the system without and with impulse drug effect. In the absence of the drug, we observe that the CTL count approaches a stable equilibrium. Under regular drug dosing, the CTL count oscillates in an impulsive periodic orbit. Assuming perfect adherence, if the drug is sufficiently strong, both infected cell and virus population approach towards extinction. In this case, the total number of uninfected cells reaches its maximum level which implies that the system approaches towards its infection-free state (Theorem 11).If we take sufficiently large impulsive intervalτ=5 days (keeping rate ω=50 fixed, as in Figure 3) or lower dosage effect ω=20 (keeping interval τ=2 fixed, as in Figure 3), in both the cases, infection remains present in the system. Thus, the proper dosage of drug and optimal dosing interval are important for infection management.


## 6. Conclusion

In this article, the role of the immunostimulant drug (mainly pidotimod) during interactions between SARS-CoV-2 spike protein and epithelial cell receptor ACE2 in COVID-19 infection has been studied as a possible drug dosing policy. To reactivate the CTL responses during the acute infection period, immune activator drugs are delivered to the host system in an impulsive mode.When the immunostimulant drug is administered, the best possible CTL responses can act against the infected or virus-producing cells to neutralize infection. This particular situation can keep the infected cell population at a very low level. In the proposed mathematical model, we have analyzed the optimal dosing regimen for which infection can be controlled.From this study, it has been observed that when the basic reproduction ratio lies below one, we expect the system to attain its disease-free state. However, the system switches from the disease-free state to the CTL-free equilibrium state when1<R0<2.957. If R0>2.957, the CTL-free equilibrium moves to an endemic state (Figure 1).Here, we have explored the immunostimulant drug dynamics by the help of impulsive differential equations. With the help of impulsive differential equations, we have studied how the effect of the maximal acceptable optimal dosage can be found more precisely. The impulsive system shows that the proper dosage and dosing intervals are important for the eradication of the infected cells and virus population which results in the control of the pandemic (Figure3).It has also been observed that the length of the dosing interval and the drug dose play a very decisive role to control and eradicate the infection. The most interesting prediction of this model is that effective therapy can often be achieved, even for low adherence, if the dosing regimen is adjusted appropriately (Figure4). Also, if the treatment regimen is not adjusted properly, the therapy is not effective at all. This approach might also be applicable to a combination of antiviral therapy.Future extension work of the combination of drug therapy should also include more realistic patterns of nonadherence (random drug holidays, imperfect timing of successive doses) and more accurate intracellular pharmacokinetics which leads towards better estimates of drug dosage and drug dosing intervals.We end the paper with the quotation: “This outbreak is a test of political, financial and scientific solidarity for the world to fight a common enemy that does not respect borders..., what matters now is stopping the outbreak and saving lives,” by Dr. Tedros, Director General, WHO [31].



---

*Source: 1352982-2020-09-01.xml*