Mathematical modeling is crucial to investigating tthe ongoing coronavirus disease 2019 (COVID-19) pandemic. The primary target area of the SARS-CoV-2 virus is epithelial cells in the human lower respiratory tract. During this viral infection, infected cells can activate innate and adaptive immune responses to viral infection. Immune response in COVID-19 infection can lead to longer recovery time and more severe secondary complications. We formulate a micro-level mathematical model by incorporating a saturation term for SARS-CoV-2-infected epithelial cell loss reliant on infected cell levels. Forward and backward bifurcation between disease-free and endemic equilibrium points have been analyzed. Global stability of both disease-free and endemic equilibrium is provided. We have seen that the disease-free equilibrium is globally stable for R0<1, and endemic equilibrium exists and is globally stable for R0>1. Impulsive application of drug dosing has been applied for the treatment of COVID-19 patients. Additionally, the dynamics of the impulsive system are discussed when a patient takes drug holidays. Numerical simulations support the analytical findings and the dynamical regimes in the systems.
Keywords: antibody response; basic reproduction number; drug holidays; epithelial cell; impulsive control; transcritical bifurcation.