key: cord-0006350-copt9zw9 authors: Zhang, Tongqian; Wang, Junling; Li, Yuqing; Jiang, Zhichao; Han, Xiaofeng title: Dynamics analysis of a delayed virus model with two different transmission methods and treatments date: 2020-01-06 journal: Adv Differ Equ DOI: 10.1186/s13662-019-2438-0 sha: 45653ee8e35339f3721c3f2ebb922228c5d241ea doc_id: 6350 cord_uid: copt9zw9 In this paper, a delayed virus model with two different transmission methods and treatments is investigated. This model is a time-delayed version of the model in (Zhang et al. in Comput. Math. Methods Med. 2015:758362, 2015). We show that the virus-free equilibrium is locally asymptotically stable if the basic reproduction number is smaller than one, and by regarding the time delay as a bifurcation parameter, the existence of local Hopf bifurcation is investigated. The results show that time delay can change the stability of the endemic equilibrium. Finally, we give some numerical simulations to illustrate the theoretical findings. Using mathematical models to help discover the mechanism of viral transmission to predict the development of infectious diseases has become the mainstream method for controlling and preventing infectious diseases [9] [10] [11] [12] . Therefore, for over a century, lots of mathematical models have been established to explain the evolution of the free virus in a body, and mathematical analysis was implemented to explore the threshold associated with eradication and persistence of the virus; for example, [13] [14] [15] [16] studied the global dynamic behavior of HIV models, [17] [18] [19] [20] [21] [22] [23] [24] analysed the global dynamics of HBV models [25] [26] [27] [28] . A general class of models describing the process of virus invading the target cells and release of the virus due to the infected cell apoptosis has been established and analyzed by Perelson et al. [29, 30] and Nelson et al. [31] as follows: where x, y, and v represent the concentrations of uninfected target cells, infected cells, and virus, respectively. Λ and d are the generation rate and mortality rate of uninfected target cells respectively, α is the infection rate, a is the mortality rate of infected cells, k and u are the generation rate and mortality rate of free virus respectively. However, the above model only considers that free viruses can infect uninfected cells by direct contact with them. Recent studies have shown that virus can be transmitted directly from cell to cell by virological synapses, i.e., cell-to-cell transmission [32] [33] [34] [35] [36] [37] [38] [39] . Spouge et al. [40] built a model to characterize this phenomenon: where C, I, M represent the concentrations of uninfected cells, infected cells, and dead cells, respectively. k I is the rate constant for cell-to-cell spread, r C is the reproductive rate of uninfected cell. μ I , μ C are the rate constants at which uninfected or infected cells die respectively, and the term k I IC represents the cell-to-cell transmission. However, research [41] shows that there is a delay between the time an uninfected cell becomes infected and when it begins to infect other uninfected cells, then Culshaw et al. [42] improved the model by introducing a distributed delay into model (2) and got the following model: where F(u) is the delay kernel. Lai et al. [43] proposed a model containing two different types of infection as follows: where T, T * , and V are the concentrations of uninfected cells, infected cells, and free virus, respectively. β 1 is the infection rate of cell-free virus transmission and β 2 is the infection rate of cell-to-cell transmission. For more details on parameters, please see [43] . The authors proved that a Hopf bifurcation can occur under certain conditions. Recently, Zhang et al. [1] proposed an ordinary differential equations virus model with both two different types of infection and cure rate as follows: where x, y, and v represent the concentrations of uninfected cells, infected cells, and free virus respectively. π is the regeneration rate of uninfected cells, d, a, u are the death rates of three kinds of cells. ρ represents the cure rate, ky is the rate at which infected cells produce free viruses. By constructing suitable Lyapunov function, the authors proved that the equilibria are globally asymptotically stable under some conditions. However, the authors did not consider the time delay in model (3) . In order to understand whether the introduction of time delay or not will change the stability of the equilibria, then motivated by the works [42, 43] and based on [1] , we further consider model (4) by introducing a discrete delay into model (3) as follows: where τ is time delay. βxy is the term of cell-to-cell transmission. αvx represents cellfree virus transmission. For more details on parameters, please see [1] . Considering the biological meanings, we analyze model (4) The paper is organized as follows. Firstly, we summarize some basic results about model (4) in Sect. 2. The local stability of the free equilibrium and Hopf bifurcation of the system are discussed in Sect. 3. We give the properties of Hopf bifurcation in Sect. 4. Finally, we perform some numerical simulations to verify the results in Sect. 5. From [1] , we can conclude some basic results and summarize them in the following theorem in this section. is the basic reproduction number. Then linearization at the original results in a characteristic equation is as follows: Clearly, it has a root λ = -d < 0. Thus we only need to analyze the distribution of roots, which determines the stability of solution of system (5) of the equation where When τ = 0, (7) reduces to Since R = (αk+βu)x 0 (a+ρ)u < 1 implies A + C > 0 and B + D > 0, we know the two roots of (8) always have a negative real part. Next, we assume that equation (7) with τ = 0 has two pure imaginary roots ±i ( > 0), which implies that the equation has at least one positive solution. Obviously, B 2 -D 2 > 0, then we can see that where a + ρ > βx 0 is used. Then we can conclude that when R < 1, equation (9) has no positive real root, which leads to that equation (7) does not have a pure imaginary root. Thus, all the roots of (7) always have negative real parts. Therefore, E 0 is locally stable. When R > 1, it is easy to see B + D < 0, then equation (7) has at least one positive root, thus E 0 is unstable. This proof is completed. Next we discuss the existence of Hopf bifurcation. For this purpose, we let to shift the equilibrium to the original. Then linearization at the original results in a characteristic equation is as follows: When τ = 0, (10) reduces to here (a + ρ)u = (βu + αk)x * is used. And Then all the roots of f (λ) = 0 have negative real parts by using the Routh-Hurwitz criterion [44] , which implies the local asymptotic stability of E * for τ = 0. For τ > 0, let iω (ω > 0) be the root of (10), we get By setting μ = ω 2 , we get where By the method in [45] , we get the following lemmas. If b 0 > c 0 and (H 1 ) hold, suppose z 1 < z 2 , then dG dz 1 < 0 and dG dz 2 > 0. Substituting Lemma 3.2 If (H1) holds, then dα dτ 1 < 0 and dα dτ 2 > 0, where λ(τ ) = α(τ ) + iω(τ ) is the root of (10) and α(τ Since dG dz 1 < 0 and dG dz 2 > 0, then we have dα dτ 1 < 0 and dα dτ 2 > 0. The proof is completed. About the existence of Hopf bifurcation, we have the following theorem. R > 1, b 0 > c 0 , and (H 1 ) hold, then system (4) undergoes a Hopf bifurcation at E * when τ = τ (j) k (k = 1, 2; j = 0, 1, 2, . . .). Furthermore, E * is stable when τ ∈ [0, τ 0 ) and unstable when τ > τ 0 . τ 0 is the Hopf bifurcation value. From Theorem 3.2, we got the sufficient conditions for the Hopf bifurcation to appear. We assume that when τ = τ 0 , system (4) produces a Hopf bifurcation at E * . Next by the normal form theory and the center manifold [46] , we try to establish the explicit formula determining the directions, stability, and period of periodic solutions bifurcating from E * at τ = τ 0 and ω(τ 0 ) = ω 0 . Let ω 0 = ω(t 0 ), τ = τ 0 + ς , ς ∈ R, then ς = 0 is the Hopf bifurcation value of model (4). Let C = C([-τ , 0], R 3 ), we havė where w t (θ ) = w(t + θ ) ∈ C and L ς is given by By using the Riesz representation theorem, we have a function ζ (θ , ς) such that, for ϕ ∈ C, where ζ (θ , ς) is a bounded variation function for [-τ 0 , 0]. And we can choose is the nonlinear part of the right-hand side of system (15) , where For ψ ∈ C 1 ([0, τ 0 ], R 3 * ) and ϕ ∈ C([-τ 0 , 0], R 3 ), define where ζ (θ ) = ζ (θ , 0). We have that H * and H = H(0) are adjoint operators. Then ±iω 0 are eigenvalues of H(0) when τ = τ 0 . Thus they are also eigenvalues of H * . Also, we can get q(θ ) = (1, q 2 , q 3 )e iω 0 θ and q * (s) =D(1, q * 2 , q * 3 )e iω 0 s , which are the eigenvectors of H and H * corresponding to iω 0 and -iω 0 , respectively, and Using the same notations as in [47] , we can obtain -(βh 2 + αh 3 )e 2iω 0 τ 0 (βh 2 + αh 3 )e 2iω 0 τ 0 0 ⎞ ⎟ ⎠ = 2D(βh 2 + αh 3 )e 2iω 0 τ 0 h * 2 -1 , and By substituting E 1 and E 2 into W 20 (θ ) and W 11 (θ ), respectively, g 21 can be expressed by the parameters. Then each g ij can be determined by the parameters. Therefore we get the following expression: Thus, we have from [47] the following theorem. In this paper, we have mainly considered the effect of time delay on the dynamics of a virus model with two different transmission methods and treatments. Our results show that the introduction of time delay has a significant effect on the dynamics of the system. However, it can be seen from [1] that when there is no time delay in system (3), the positive equilibrium E * of system (3) is asymptotically stable if it exists. The appearance of the time delay causes the positive equilibrium of the model (4) to be inverted from stable to unstable, and a periodic solution of small amplitude is generated near the positive equilibrium E * . Biologically, the number of healthy, infected, and free viruses exhibits periodic changes. Next, we take some numerical simulations to validate our main results. We set the basic parameters as follows [48] : d = 0.2, β = 0.000024, α = 0.000024, ρ = 0.2, a = 0.15, k = 150, u = 0.2. Firstly, we set π = 2, direct calculations with Maple 14 show that R = 0.514971 < 1 and E 0 = (10, 0, 0). By Theorem 3.1 the virus-free equilibrium E 0 of the system is stable G z * 1 = -0.008284 < 0, then equation (12) has two positive roots z 1 = 0.0089 and z 2 = 0.3680, and equation (11) has two positive roots 1 = 0.09441603258 and 2 = 0.6066529567. Therefore we have the Hopf bifurcation value τ 0 2 = 3.828340005. If we set τ = 3.8 < 3.828340005, by Theorem 3.2, the virus equilibrium E * is asymptotically stable (see Figs. [5] [6] [7] [8] . If we set τ = 3.9 > 3.828340005, by Theorem 3.2, the virus equilibrium E * is unstable and periodic oscillations occur (see Figs. 9-16 with different initial values). Data sharing not applicable to this article as all data sets are hypothetical during the current study. Global dynamics of a virus dynamical model with cell-to-cell transmission and cure rate WHO: World health organization report: research for universal health coverage Periodic solution and ergodic stationary distribution of stochastic SIRI epidemic systems with nonlinear perturbations Division of STD Prevention: Prevention of Genital HPV Infection and Sequelae: Report of an External Consultants' Meeting. Department of Health and Human Services Condom use and the risk of genital human papillomavirus infection in young women Human papillomavirus, smoking, and sexual practices in the etiology of anal cancer The burden of HPV-related cancers Human papillomavirus and rising oropharyngeal cancer incidence in the United States Influence of asymptomatic infections for the effectiveness of facemasks during pandemic influenza Global analysis of a stochastic TB model with vaccination and treatment Dynamics analysis of stochastic epidemic models with standard incidence Analysis of a stochastic recovery-relapse epidemic model with periodic parameters and media coverage A differential equation model of HIV infection of CD4 + T-cells with cure rate Global stability of a viral infection model with two delays and two types of target cells Global stability of an HIV-1 infection model with saturation infection and intracellular delay Role of cd4 + t-cell proliferation in HIV infection under antiretroviral therapy Viral clearance without destruction of infected cells during acute HBV infection Global properties of an improved hepatitis B virus model Modeling complex decay profiles of hepatitis B virus during antiviral therapy Global stability for a dynamic model of hepatitis B with antivirus treatment Analysis of a stochastic HIV-1 infection model with degenerate diffusion Influence of media intervention on AIDS transmission in MSM groups Modeling inhibitory effect on the growth of uninfected T cells caused by infected T cells: stability and Hopf bifurcation Global behavior of delay differential equations model of HIV infection with apoptosis Extinction and persistence of a stochastic SIRS epidemic model with saturated incidence rate and transfer from infectious to susceptible Dynamics of a stochastic SIS epidemic model with nonlinear incidence rates Dynamics of a stochastic SIR model with both horizontal and vertical transmission Modeling cell-to-cell spread of HIV-1 with nonlocal infection Mathematical analysis of HIV-1 dynamics in vivo HIV-1 dynamics in vivo: virion clearance rate, infected cell life-span, and viral generation time Mathematical analysis of delay differential equation models of HIV-1 infection Recruitment of HIV and its receptors to dendritic cell-T cell junctions Avoiding the void: cell-to-cell spread of human viruses The role of cell-to-cell transmission in HIV infection Cell-to-cell spread of HIV-1 occurs within minutes and may not involve the participation of virus particles Cell-to-cell transmission of retroviruses: innate immunity and interferon-induced restriction factors The immune control and cell-to-cell spread of human T-lymphotropic virus type 1 Cell-to-cell contact via measles virus haemagglutinin-CD46 interaction triggers CD46 downregulation Measles virus: both the haemagglutinin and fusion glycoproteins are required for fusion HIV-1 infection kinetics in tissue cultures Viral dynamics in vivo: limitations on estimates of intracellular delay and virus decay A mathematical model of cell-to-cell spread of HIV-1 that includes a time delay Modeling cell-to-cell spread of HIV-1 with logistic target cell growth On the zeros of a third degree exponential polynomial with applications to a delayed model for the control of testosterone secretion Theory and Applications of Hopf Bifurcation Introduction to Functional Differential Equations Stability and bifurcation analysis for an improved HIV model with time delay and cure rate The authors declare that there is no conflict of interests regarding the publication of this paper. All authors read and approved the final manuscript. Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.