key: cord-1026193-cp6dp5ud authors: Hattaf, Khalid; Dutta, Hemen title: Modeling the dynamics of viral infections in presence of latently infected cells date: 2020-06-04 journal: Chaos Solitons Fractals DOI: 10.1016/j.chaos.2020.109916 sha: 15cb8cdb5ff19ca2dc34036339aad3b8bdfadbee doc_id: 1026193 cord_uid: cp6dp5ud The study aims to develop a new mathematical model in order to explain the dynamics of viral infections in vivo such as HIV infection. The model includes three classes of cells, takes into account the cure of infected cells in latent period and also incorporates three modes of transmission. The mention modes are modeled by three general incidence functions covering several special cases available in the literature. The basic properties of the model as well as its stability analysis have been carried out rigorously. Further, an application is given and also numerical simulation results have been incorporated supporting the analytical results. Dangerous viruses have been appearing from time to time attacking the human body and causing several health problems, deaths and socio-economic damage worldwide. For instance, a novel coronavirus (nCoV) was identified on January 7, 2020 responsible for the initial cases of pneumonia occurred in Wuhan, China. This virus has recently been reported as extreme acute respiratory syndrome coronavirus 2 (SARS-CoV-2) by the Coronavirus Study Group of the International Committee on Taxonomy of Viruses [7] . Because its genetic sequence is very similar to that of severe acute respiratory syndrome (SARS), another coronavirus that appeared for first time in 2002 and caused a pandemic with more than 80 0 0 infected people and 800 deaths. Another coronavirus that causes severe disease in humans is middle east respiratory syndrome (MERS), identified for the first time in 2012 in the Middle East having connection with camels. As of 24 March 2020, more than 372,757 cases of SARS-CoV-2 have been confirmed in 196 countries and territories, including 195,511 cases in European region, and the number of deaths from the virus has reached 16,231 globally, according to the World Health Organization (WHO) [27] . Also, already existing viruses like hepatitis B virus (HBV) and hu- * Corresponding author. E-mail addresses: k.hattaf@yahoo.fr (K. Hattaf), hemen_dutta08@rediffmail.com (H. Dutta). man immunodeficiency virus (HIV) are still causing large numbers of deaths worldwide. As per an estimation in 2015, 257 million people were living with chronic hepatitis B infection (defined as hepatitis B surface antigen positive), and death of 887,0 0 0 people, mostly due to cirrhosis and hepatocellular carcinoma (i.e. primary liver cancer) [28] . Also, 770,0 0 0 people died from HIV-related causes and 1.7 million people were newly infected in 2018 [29] . The rapid spread of viruses has been attracting the attention of researchers to carry out new researches to understand various associated aspects. It has been shown biologically that viruses can spread by both modes, viz., virus-to-cell mode via the extracellular environment and through cell-to-cell transmission mode by direct cell-cell contact [13, 14, 23] . During cell-to-cell transmission, it is possible for several viral particles to be transferred from infected to uninfected cells via structures called infectious or virological synapses [13] . In the case of HIV infection, cell-to-cell transmission may be 10-10 0 0 times more effective than the passive dissemination of virons via the extracellular environment [2, 4, 5] . In [30] , HIV is used as a model virus in understanding the processes that consequence in viral spread via virus-to-cell or cell contact dependent mode. In the same work, it is demonstrated that HIV can spread by virus-to-cell mode if all steps of the viral replication cycle are efficient. If any of these steps is impaired, then co-cultures of HIV donor cells with specific target cells allow effective viral spreading. Latently infected cells are those cells which include the viruses but do not producing them due to time lag between both the moment of infection and the moment when the infected cell becomes active to create new infectious viruses. For some viruses like herpesvirus, latency is an essential mechanism for immune evasion and viral persistence [16, 17] . One of the major hurdles in eliminating HIV infection completely is the presence of a small reservoir of latently infected cells that can survive through long-term suppression of viral replication by extremely active anti-retroviral therapy (HAART) [6, 21, 22] . A recent study found that cell-to-cell transmission contributes to the production of latently infected cells [1] . This mode of viral transmission is found to contribute to both the spread and pathogenesis of the virus in vivo and to the production and maintenance of latent reservoirs that represent a serious barrier to the curing of the infection [19] . The above biological studies inspire to propose the following model: where S ( t ), L ( t ), A ( t ) and V ( t ) respectively indicate the concentration of susceptible cells, latently infected cells, actively infected cells and free virus particles at time t . It is also considered that susceptible cells are generated at a constant rate of λ, die at a rate of dS and become infected either through free virus at a rate of p ( S, V ) V , or through direct contact with actively infected cells at a rate of q ( S, A ) A , or by direct contact with latently infected cells at a rate of h ( S, L ) L . The term p(S, V ) V + q (S, A ) A + h (S, L ) L then describes the total infection rate of susceptible cells. It is assumed a fraction η of infected cells to be latently infected cells and the remaining (1 − η) becomes actively infected cells, where 0 < η ≤ 1. We take e and a as the parameters for the death rates of latently and actively infected cells. Also, r, k and μ are the parameters respectively for the production rate of virus from actively infected cells, the rate at which the latently infected cells are transmitted to actively infected cells, and the clearance rate of virus. Fig. 1 is an illustration of the schematic diagram of model (1) . As in Hattaf and Yousfi [10] , 11 ], the general incidences p ( S, V ), q ( S, A ) and h ( S, L ) are assumed to be continuously differentiable in the interior of IR 2 + and also satisfy the following hypotheses: It should be noted that model (1) extends as well as generalizes many special cases available in the literature. As example, one can get the latent infection model of Pankavich [18] with p(S, V ) = β 1 S, q (S, A ) = 0 , h (S, L ) = 0 and ρ = 0 , where β 1 is the virus-to-cell infection rate. Also, the model of Callaway and Perelson [3] with an overall drug efficacy. Furthermore, the model (1) includes the class of HIV infection models for different rates of infection and latently infected cells as investigated in Wang et al. [26] . The model (1) reduces to when all of infected cells become latently infected cells (i.e. η = 1 ). This system includes the following viral models: • Model of Rong et al. [20] with p(S, The paper has been organized as follows: In Section 2 , the existence of equilibria and the well-posedness of the model including nonnegativity and boundedness of the solutions have been established. The stability analysis of the equilibria has been presented in Section 3 . An application and illustrative numerical simulations have been given in Section 4 . The Section 5 is for conclusion in terms of both mathematical and biological perspectives of the findings. This section first establishes the well-posedness of our model (1) by proving the nonnegativity and boundedness of solutions. Then we derive the threshold parameter for the existence of equilibria. Then IR 5 + is positively invariant with regard to (1) . Next, we prove the boundedness of solutions. Denote It is obvious that model (1) has always one infection-free equilibrium Then we define the basic reproduction number of (1) as follows: . This number can be rewritten as follows The other equilibrium of model (1) satisfies the following system kA − μV = 0 . . This section deals with the stability analysis of two equilibria E p and E * . For more details, the local stability is obtained by analyzing the characteristic equation. However, the global stability is established by using two techniques, the first is based on the definition which consists to prove that the solution of the model converges to the equilibrium point independently of the initial conditions. The second technique is based on the direct Lyapunov method which aims to construct an appropriate Lyapunov functional. where It is obvious that ξ = −d is a negative root of Eq. (8) . The other roots are determined by the following equation: If R 0 < 1 , then a 1 > 0, a 2 > 0, a 3 > 0 and By the Routh-Hurwitz Theorem [8] , we know that all roots of (9) have negative real parts. Thus, the infection-free equilibrium E p is locally asymptotically stable if R 0 < 1 . If R 0 > 1 , then a 3 < 0. So, Eq. (9) has at least one positive root. This implies that the infection-free equilibrium E p becomes unstable when R 0 > 1 . The following theorem characterizes the global stability of the infection-free equilibrium E p when R 0 < 1 . Proof. We have the following inequality from the first equation of (1) Then for arbitrary > 0, there exist a t 1 > 0 such that From the last three equations of (1) , for t ≥ t 1 , we have Consider the following comparison system This system can be rewritten as follows The eigenvalues of the matrix M satisfy the following equation: It follows from the Routh-Hurwitz Theorem that all eigenvalues of M have negative real parts. Therefore, From the first equation of (1) , for t ≥ t 1 , we can get This completes the proof of Theorem 3.2 . Next, we investigate the global asymptotic stability of the chronic infection equilibrium E * ( S * , L * , A * , V * ) by assuming that R 0 > 1 and the incidence functions p, q and h satisfy, ∀ S, L, A, V > 0, the following further hypothesis: Since the incidence function p ( S, V ) is strictly monotonically increasing with respect to S , we deduce that Since ( x ) ≥ 0, we have dW dt ≤ 0 with equality if and only if S = S * , L = L * , A = A * and V = V * . From LaSalle's invariance principle [12] , we conclude that the chronic infection equilibrium E * is globally asymptotically stable when R 0 > 1 . In this section, we first apply our main results to the following HIV infection model: Then the assumption ( H 4 ) is also satisfied. In this special case, the basic reproduction number has the following form: According to Theorems 3.2 and 3.3 , we deduce the following result. Corollary 4.1. Let R 0 defined by (11) . (i) If R 0 < 1 , then the infection-free equilibrium E p of model (10) is globally asymptotically stable. (ii) If R 0 > 1 , then the infection-free equilibrium E p becomes unstable and the chronic infection equilibrium E * of model (10) is globally asymptotically stable for ρ = 0 and η = 1 . Now, we carry out numerical simulations to confirm the above analytical results concerning the dynamics of model with the parameter values given in Table 1 . Firstly, we choose β 1 = 5 . 2 × 10 −5 , ρ = 0 . 01 and η = 0 . 7 . By a simple computation, we have R 0 = 0 . 9733 < 1 . Hence, model (10) has an infection-free equilibrium E p (719.4245, 0, 0, 0). Fig. 2 shows that the trajectories of (10) starting from different initial values tend to E p . This supports the global stability result given by Corollary 4.1 (i). Secondly, we choose β 1 = 2 . 4 × 10 −4 , ρ = 0 and η = 1 . In this case, we have R 0 = 6 . 7313 > 1 . Fig. 3 shows that the trajectories of model (10) Thirdly, we choose β 1 = 2 . 4 × 10 −4 , ρ = 0 . 01 and η = 0 . 7 . Then R 0 = 7 . 6442 > 1 . Fig. 4 shows that E * (618.5765, 5.4708, 3.5012, 58.3445) is globally asymptotically stable despite ρ = 0 and η = 1. Thus, the conditions ρ = 0 and η = 1 are not necessary for the global asymptotic stability of the chronic infection equilibrium E * . We will prove analytically this open problem in future work. The work developed a model by using tools and techniques of mathematics for infection by viruses with latently infected cells and cell-to-cell transmission. The well-posedness of the model, including nonnegativity and boundedness of solutions have been established. It also identified two equilibria, viz. an infection-free equilibrium E p and chronic infection equilibrium E * . It is proved that E p is globally asymptotically stable if the basic reproduction number is less than one, which implies that the virus is eradicable and the infection is curable. And, when the basic reproduction number is greater than one, E p becomes unstable and E * is globally asymptotically stable for ρ = 0 and η = 1 . In this case, the virus persists in the host. On the other hand, the models and results presented in Callaway and Perelson [3] , Hu et al. [9] , Maziane et al. [15] , Pankavich [18] , Rong et al. [20] , Sun et al. [24] , Wang et al. [25] , 26 ] are improved and generalized. It is observed from the explicit formula of the basic reproduction number R 0 given in (3) that R 0 is as sum of the basic reproduction numbers for the classical virus-to-cell infection mode R 01 , cell-to-cell transmission mode via actively infected cells R 02 and cell-to-cell transmission mode via latently infected cells R 03 . Hence, ignoring the cell-to-cell transmission in viral infections dynamics will underestimate R 0 . Further, the study of this paper can be of help for the biologists to develop effective treatments against viral infections in order to reduce the patients' basic reproduction number to less than one and eliminate the latent reservoirs which represent the main obstacle to the eradication of virus from the human body. Author 1 discussed the research problem with Author 2 first. Author 2 gave his opinion on the formulation of the model and possible solutions. All the authors contributed in developing the main results and their proofs. All authors contributed in each section of the paper and approved the paper for submission. The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. HIV-1-infected CD4 + T cells facilitate latent infection of resting CD4 + T cells through cell-cell contact Rapid and efficient cell-to-cell transmission of human immunodeficiency virus infection from monocyte-derived macrophages to peripheral blood lymphocytes HIV-1 infection and low steady state viral loads Predominant mode of human immunodeficiency virus transfer between T cells is mediated by sustained Env-dependent neutralization-resistant virological synapses Quantitation of human immunodeficiency virus type 1 infection kinetics Identification of a reservoir for HIV-1 in patients on highly active antiretroviral therapy The species severe acute respiratory syndromerelated coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2 Routh-Hurwitz theorem, in tables of integrals, series, and products Analysis of a CD4 + T cell viral infection model with a class of saturated infection rate A generalized virus dynamics model with cell-to-cell transmission and cure rate Qualitative analysis of a generalized virus dynamics model with both modes of transmission and distributed delays The stability of dynamical systems Virus entry: open sesame Virus cell-to-cell transmission Dynamics of a class of HIV infection models with cure of infected cells in eclipse stage Viral latency and its regulation: lessons from the g-herpesviruses Towards an understanding of the herpes simplex virus type 1 latency-reactivation cycle The effects of latent infection on the dynamics of HIV Mechanisms of HIV-1 cell-to-cell transmission and the establishment of the latent reservoir Modeling within-host HIV-1 dynamics and the evolution of drug resistance: trade-offs between viral enzyme function and drug susceptibility An integrated overview of HIV-1 latency Heterogeneous clearance rates of long-lived lymphocytes infected with HIV: intrinsic stability predicts lifelong persistence Avoiding the void: cell-to-cell spread of human viruses Global stability of infection-free state and endemic infection state of a modified human immunodeficiency virus infection model Global dynamics for viral infection model with Beddington-DeAngelis functional response and an eclipse stage of infected cells Global dynamics of a class of HIV-1 infection models with latently infected cells COVID-19) Situation Report-64 Cell-to-cell transmission can overcome multiple donor and target cell barriers imposed on cell-free HIV The authors would like to thank Karmina K Ali, University of Zakho for reading the manuscript and help in formatting the paper. The authors are also grateful to the editor and anonymous reviewers for their careful reading and helpful suggestions.