Abstract
In this investigation, we propose and analyze a virus dynamics model with multi-stages of infected cells. The model incorporates the effect of both humoral and cell-mediated immune responses. We consider two modes of transmissions, virus-to-cell and cell-to-cell. Multiple intracellular discrete-time delays have been integrated into the model. The incidence rate of infection as well as the generation and removal rates of all compartments are described by general nonlinear functions. We derive five threshold parameters which determine the existence of the equilibria of the model under consideration. A set of conditions on the general functions has been established which is sufficient to investigate the global stability of the five equilibria of the model. The global asymptotic stability of all equilibria is proven by utilizing Lyapunov function and LaSalle's invariance principle. The theoretical results are illustrated by numerical simulations of the model with specific forms of the general functions.