Global properties of an age-structured virus model with saturated antibody immune response, multi-target cells and general incidence rate Ángel Cervantes-Pérez, Eric Ávila-Vales arxiv:1712.564v1 [math.ds] 14 Dec 217 Facultad de Matemáticas, Universidad Autónoma de Yucatán, Anillo Periférico Norte, Tablaje 13615, C.P. 97119, Mérida, Yucatán, Mexico Abstract Some viruses, such as human immunodeficiency virus, can infect several types of cell populations. The age of infection can also affect the dynamics of infected cells and production of viral particles. In this work, we study a virus model with infectionage and different types of target cells which takes into account the saturation effect in antibody immune response and a general non-linear infection rate. We construct suitable Lyapunov functionals to show that the global dynamics of the model is completely determined by two critical values: the basic reproduction number of virus and the reproductive number of antibody response. 1 Introduction In recent times, several mathematical models have been proposed in order to try to understand the mechanism of virus infections. These models often describe the changes through time in the concentration of infected and uninfected target cells and viral particles in the blood of an infected individual [1 16]. Modelling the effect of antibody immune response in the neutralization of virus is a very important topic for research since it can provide useful insights into the dynamics of the infection and offer suggestions for clinical treatment. The production of new virions by an infected cell does not occur at a constant rate, and the death rate of infected cells can vary with the time the cell has been infected. Thus, the incorporation of an age structure in models allows us to have a more realistic picture of produced viral particles and the mortality of infected cells; see, for example, the models studied in [1 6]. A recent work by Duan and Yuan [3] introduced a model with age-structure and antibody immune response, described by the system iθ, t θ dt t βt tv t λ dt t 1 + αv t iθ, t + δθiθ, t t dv t P θiθ, t cv t qatv t dat katv t h + At bat, 1
with the initial and boundary conditions i, t βt tv t 1 + αv t, T T, iθ, i θ, V V, A A. In this model, T, i, V, A denote the concentration of uninfected susceptible host cells, infected host cells, free virus particles, and antibody responses released from B cells, respectively. The variable θ is the infection age, i.e., the time that has elapsed since a virion has infected a cell, while iθ, t is the density of infected T cells with infection age θ at time t. The functions P θ and δθ are respectively the time-since-infection structured virion production rate and the death rate of productively infected cells. The authors in [3] showed that the global dynamics of this model is completely determined by the basic reproduction number R and the viral reproductive number R 1. Most mathematical models for human immunodeficiency virus HIV focus on the infection of CD4 + T cells. However, macrophages and dendritic cells are also susceptible to be infected with HIV [4, 5]. Therefore, we need to incorporate in our virus dynamic models a multi-group component to study virus infection in different populations of cells. Viral infection models dealing with the interaction of virus with more than one class of target cells have been studied in [7 13] and references cited therein. When modelling the rate of infection of susceptible target cells by contact with virus, the bilinear incidence term βt V is used frequently in literature [1, 2, 14]. However, there is certain debate that this rate is insufficient to realistically describe the infection process. Some authors have incorporated in their models several types of non-linear incidence functions, including saturated Holling type II incidence βt V/1 + αv [3, 5, 15], Beddington-DeAngelis functional response βt V/1+aT +bv [6, 7], or even more general functions, e.g., of the forms ct fv [16] or ht, V [4, 8, 9]. In this paper, we propose a viral infection model that includes age-structure, multitarget cells, and a general non-linear rate of viral infection. We establish the global properties of its equilibria using semigroup methods, uniform persistence and the construction of Lyapunov functionals. The rest of this paper is organized as follows: in Section 2 we present the equations and assumptions for our model, we calculate the basic reproduction number R and show boundedness and non-negativity of solutions. In Section 3 we rewrite the system as a semilinear Cauchy problem, we prove that the semiflow generated by the system has a global compact attractor and that it is uniformly persistent when R > 1. In Section 4 we determine the steady states of the system and deal with their local stability. In Section 5 we construct several Lyapunov functionals to determine the global stability of equilibria. In Section 6 we present a special case in which the model can be reduced to a system of delay differential equations. In Section 7 we provide some numerical simulations which illustrate the different kinds of global behaviour the model can have. Finally, in Section 8 we present some concluding remarks. 2 Description of the model We consider a within-host viral dynamics model for HIV with multiple target cell populations, which takes into account the saturation effects of antibody immune response and incorporates a general incidence function h j T j, V. The model is given by the following 2
system of differential equations: i j θ, t θ with initial conditions dt j t + i jθ, t t dv t dat λ j T j t h j T j t, V t δ j θi j θ, t P j θi j θ, t cv t qatv t katv t h + At bat 2.1 and boundary conditions T j T j, i j θ, i jθ, V V, A A 2.2 i j, t h j T j t, V t. 2.3 This model considers n classes of target cells, denoted by the subscript j 1,..., n. For each class, T j represents the population of uninfected target cells, while i j θ, t denotes the population of infected cells with infection age θ at time t. We denote by V the concentration of free virus particles, and by A the number of antibody responses released from B cells. We assume that uninfected cells in the j-th class are recruited at a rate λ j and die at a rate. Uninfected cells become infected from contact with the virus according to the incidence function h j T j t, V t. Infected cells in the j-th class die at a rate δ j θ, which depends on the infection age θ. The function P j θ represents the viral production rate of an infected cell from the j-th class with infection age θ, while the parameter c is the viral clearance rate. The virus is removed by the immune system at a rate qa. The antibody responses are generated by contact with the virus at a rate katv t/h + At and vanish at rate b. All parameters are assumed positive. Furthermore, we assume that the functions δ j and P j satisfy δ j, P j L +,, R \ { L } and δ j θ δ min > for all θ, j 1,..., n. Let P max max {ess sup P j j 1,..., n}, where ess sup P j denotes the essential supremum of P j θ for θ,. We make the following hypotheses on the incidence functions h j, j 1,..., n: H1 h j T j, V V h j T j, V with h j, h j C 2 R 2 + R + and hj, V h j T j, for all T j, V. H2 h j T j, V is strictly increasing with respect to both variables when T j, V >. H3 h j T j, V for all T j, V and h j T j, V as V. 2.1 Basic reproduction number and boundedness of solutions For each infected cell in the j-th class, the probability of still being infected after θ time units is given by σ j θ e θ δ jτ dτ, 3
so the total number of virions produced by an infected cell from the j-th class in its entire life span is N j P j θσ j θ. Hence the corresponding basic reproduction number for system 2.1 when the j-th class is the unique class of target cells is given by R j N j c h j λj,, j 1,..., n. Therefore, we can define the basic reproduction number R of virus for the model as R R j 1 c h j λj N j,, 2.4 which represents the expected number of viruses that one virion produces in a fully uninfected cell population. The second equation of 2.1 is a linear transport equation with decay, so it can be solved with the initial conditions 2.2 and boundary conditions 2.3 by the method of integration along the characteristic lines t θ constant. The solution is given by { hj Tj t θ, V t θ σ j θ if θ < t, i j θ, t i j θ t σ jθ σ j θ t if θ t, 2.5 where σ j θ e θ δ jτ dτ. We will now prove that the model is biologically well-posed by showing the nonnegativity and boundedness properties of the system. Theorem 2.1. The solutions of system 2.1 with non-negative initial conditions 2.2 and boundary conditions 2.3 remain non-negative for t and are positively bounded. Proof. Let T 1, i 1,..., T n, i n, V, A be a solution of 2.1 with non-negative initial conditions. Suppose that T j t loses its non-negativity for some j and let t inf{t T j t < }. By continuity of solutions, we have T j t and dt j t. But the first equation of 2.1 implies that dt j t λ j >, a contradiction. Therefore, no such t exists and thus, T j t > for all t. Next, we show that V t > and i j θ, t for all t, θ. Define t 1 min { inf{t V t < }, inf{t i 1, t / L 1 +, },..., inf{t i n, t / L 1 +, } }. Suppose first that t 1 inf{t i j, t / L 1 +, } for some j. Then 2.5 implies that i j θ, t 1 for all θ, which is a contradiction. Hence, t 1 inf{t V t < }. At the time instant t 1, we have dv t 1. However, since i j, t 1 > for all j, the third equation of 2.1 implies that dv t 1 P j θi j θ, t 1 >, which leads to a contradiction. Thus, V t > and i j θ, t for all t, θ. Finally, from the last equation of the model we have At Ae t kv τ h+aτ b dτ. 4
Therefore, the solution remains non-negative on its maximal interval of existence. Next, we consider the boundedness of solutions. From the first equation of 2.1, we have dt j λ j T j for j 1,..., n. This implies that Let Gt lim sup T j t λ j. t T j t + i j θ, t be the total number of cells uninfected and infected at time t. Then dgt dtj t i j θ, t + t [ λ j T j t h j T j t, V t + [ λ j T j t λ j ] δ j θi j θ, t [ T j t + δ min ] i j θ, t, since δ j θ δ min. Let d min{d 1,..., d n, δ min }. Then and thus, dgt λ j d lim sup t n where M λ j /d. so Now, from the equation for dv dv t when t is large enough. Thus, [ T j t + i jθ, t δ j θi j θ, t θ ] i j θ, t λ j dgt, i j θ, t lim sup Gt M, t in 2.1, we have P max dv t lim sup t P j θi j θ, t cv t i j θ, t cv t, P max M cv t V t P maxm. c ] 5
Finally, from the last equation in 2.1, we have for large t, so dat kv t bat kp maxm c lim sup t which completes the proof of the theorem. At kp maxm, bc bat Let Ω { T1, i 1,..., T n, i n, V, A T j λ j, T j + n i j θ [ R + L 1 +,, R ] R 2 + } M, V P maxm, A kp maxm. c bc Then Ω is a positively invariant set for system 2.1. From now on, we will always assume that the initial value T 1, i 1,..., T n, i n, V, A is in Ω. 3 Integrated semigroup formulation We will now reformulate system 2.1 with the initial condition 2.2 as a semilinear Cauchy problem. In order to take into account also the boundary condition 2.3, we need to expand the state space as follows. Let M R R L 1,, R, P R {} L 1,, R, N R {} W 1,1,, R n, X M R R, where W 1,1 denotes a Sobolev space. Consider the linear operator A : DomA X X given by T 1 d 1 T 1 i 1 i 1 i 1 δ 1θi 1.. A T n d n T n, i n i n i n δ n θi n V cv A ba with DomA n i1 P R R. Then n DomA N R R : X i1 i1 6
and we can consider the non-linear map F : X X given by T 1 λ 1 h 1 T 1 t, V t h1 T 1 t, V t i 1 L1.. F T n λ n h n T n t, V t. hn T n t, V t i n L1 V n P j θi j θ, t qatv t A Let katv t h+at T ut T 1 t,,..., T i 1, t n t,, V t, At, i n, t P + R + {} L,, R n, X + P + R + R +. Thus, we can rewrite system 2.1 with the boundary and initial conditions as the following abstract Cauchy problem: i1 dut Aut + F ut for t, with u x X +. 3.1 By applying the results given in [17] and [18], we can conclude the following theorem. Theorem 3.1. System 2.1 generates a unique continuous semiflow {Ut} t on X + that is asymptotically smooth and bounded dissipative. Furthermore, the semiflow {Ut} t has a global compact attractor in X +, which attracts the bounded sets of X +. 3.1 Uniform persistence We will now establish the uniform persistence of system 2.1, using arguments that are highly motivated by those in [6] and [3]. Define M T1,, i 1,..., T n,, i n, V, A X + V + i j θ >, and M X + \ M. For j 1,..., n, let ϑ j sup{θ, P j θ > }. Note that ϑ j can possibly be +. Lemma 3.2. M is positively invariant under the semiflow {Ut} t generated by system 2.1. Moreover, the infection-free equilibrium E λ 1 λ d 1,, 2 λ d 2,,..., n d n,,, is globally asymptotically stable restricted to M. Proof. Let y M and T 1 t,, i 1, t,..., T n t,, i n, t, V t, At Uty. Then we have i j θ, t θ + i jθ, t t δ j θi j θ, t 7
dv t i j, t h j T j t, V t, V V. P j θi j θ, t cv t qatv t i j θ, i jθ, Since T j t λ j for large t, then i j θ, t ĩ j θ, t and V t solution of the system Ṽ t, where ĩ, Ṽ is the ĩ j θ, t θ + ĩ j θ, t δ j θĩ j θ, t t dṽ t P j θĩ j θ, t cṽ t ĩ j, t h j λj, Ṽ t, Ṽ V. ĩ j θ, i jθ 3.2 From the first and third equation of 3.2, we have { λj hj d ĩ j θ, t j, Ṽ t θ σ j θ if θ < t, i j θ t σ jθ σ j θ t if θ t and by substituting this expression in the second equation of 3.2, we obtain where dṽ t t h j λj, Ṽ t θ P j θσ j θ + F j t t P j θi σ j θ jθ t σ j θ t. 3.3 F j t cṽ t, 3.4 For each j, we will prove that F j t for all t. In fact, if t ϑ j, then P j θ for θ > t ϑ j, so F j t t P j θi j θ t. If t < ϑ j, then F j t t P j θi jθ t ϑj P max i jθ t + t ϑ j P j θi jθ t Since y M, then ϑ j t i j θ t V + n i j θ, and also P jθ for θ > ϑ j, so the two terms in the right-hand side of the aforementioned inequality are equal to zero. Thus, Fj t. Accordingly, 3.4 is an autonomous equation and has a unique solution Ṽ t. dat Hence V t Ṽ t and bat for large t, which implies that V t and At approach as t. Now, it follows from 3.3 that ĩ j θ, t for θ t. For t < θ, we have ĩ j θ, t L 1 σ j θ i jθ t σ j θ t e δ mint i j L 1. L 1 Therefore, i j θ, t as t. Since i j θ, t ĩ j θ, t, then Uty approaches E as t and the proof is complete. 8
Theorem 3.3. If R > 1, then the semiflow {Ut} t is uniformly persistent with respect to the pair M, M, i.e., there exists ɛ > such that for each y M, lim inf t d Uty, M ɛ. Furthermore, there exists a compact subset A M that is a global attractor for {Ut} t in M. Proof. Since E is globally asymptotically stable in M, by Theorem 4.2 in [19] we only need to investigate the behaviour of the solutions starting in M in some neighbourhood of E. Then, we will show that W s {E } M, where W s {E } { y X + lim t Uty E }. Assume by contradiction that there exists y W s {E } M. It follows that there exists t > such that V t + n i j θ, t >. Using the same argument as in the proof of Lemma 3.6i in [2], we have that V t > for t and i j θ, t > for any θ, t. Define the function γ j x x P j θe θ x δ jτ dτ. Note that γ j is bounded and satisfies γ j x δ jxγ j x P j x for all x. Consider the function which satisfies dφt Φt cv t γ j θi j θ, t + V t, N j h j T j t, V t cv t qat c 1. Since y W s {E }, we have that T j t λ j, V t, and At as t. Since R > 1, then lim t N j h j T j t, V t cv t qat c 1 N j c h j λj, 1, so Φt is non-decreasing for sufficiently large t. Thus, there exists t 1 > such that Φt Φt 1 for all t > t 1. Since Φt 1 >, this prevents the function V t, i j, t from converging to, L 1 as t, which contradicts that T j t λ j. This completes the proof. 4 Steady states and local stability In this section, we will study the existence and local properties of the steady states for system 2.1. In order to prove the existence of infected equilibria, we first need to prove the following results. 9
Lemma 4.1. For j 1,..., n, T j λ j and V, the equation λ j T j h j T j, V 4.1 has a unique solution given by T j f j V, where f j : R +, λ j ] is a decreasing function with f j λ j. Proof. Let ϕt j, V λ j T j h j T j, V. If V, it is clear that the only solution of 4.1 is T j λ j. If V >, then assumption H2 implies that ϕt j, V is strictly decreasing with respect to T j. Since ϕ, V λ j > and ϕ λ j, V < λ j λ j, then ϕt j, V has a solution T j f j V in the interval, λ j, and it is unique due to monotonicity of ϕ. Thus, the function f j is well defined and positive. Moreover, differentiation of 4.1 gives + h j T j T j, V dt j h j T j, V dv, so f jv dt j dv hj T j, V / + h j T j, V, T j which is negative for V > by assumption H2, Therefore, f j is decreasing. Lemma 4.2. The function ψ j : R + R + defined by ψ j V h j f j V, V is strictly decreasing for V, and ψ j h j λj,. Proof. Let V 1 < V 2. Then f j V 2 < f j V 1 by Lemma 4.1. By H2 this implies that h j f j V 2, V 2 < h j f j V 1, V 2, so h j f j V 2, V 2 < h j f j V 1, V 2. By H3 we also have that h j is non-increasing with respect to the second variable, so h j f j V 1, V 2 h j f j V 1, V 1. Therefore, by transitivity, h j f j V 2, V 2 < h j f j V 1, V 1, i.e., ψ j V 2 < ψ j V 1. For V, we have λj h j f j, V h j d ψ j h j f j, lim lim j, V λj h j, V V V V hence the result. h j λj,, We can now define the viral reproduction number R of model 2.1, which is given by R 1 c N j h j fj bh k, bh k 4.2 and can be interpreted as the average number of antibodies that are activated by the introduction of a single virion within the host, under the condition that there is no previous antibody response. Since bh k >, Lemma 4.2 implies that h j fj bh k, bh k < hj f j, h j λj,, and by using the expression in 2.4 for the basic reproduction number R we can see that R < R. These two numbers, R and R, determine whether system 2.1 has one, two, or three steady states, as asserted in the following theorem. 1
Theorem 4.3. System 2.1 always has an infection-free steady state E T 1,, T 2,,..., T n,,,, where T j λ j/. In addition, when R > 1, system 2.1 has a unique immune-free infected steady state E T 1, i 1 θ,..., T n, i nθ, V,, where T j f j V, i jθ h j fj V, V σ j θ, and V is the unique positive solution of N j h j fj V, V c. 4.3 Also, if R > 1, then system 2.1 has a unique antibody-immune infected steady state Ê ˆT1, î 1 θ,..., ˆTn, î n θ, ˆV, Â, where ˆT j f j ˆV, î j θ h j fj b k h + Â, b k h +  b σ j θ, ˆV h + Â, k and  is the unique positive solution of N j h j fj b k h + Â, b k h +  c qâ. 4.4 Proof. It is clear that for 2.1 the infection-free steady state E always exists. To obtain the immune-free infected steady state E, we have to assume that A and V. We obtain the system λ j T j h j T j, V di j θ δ j θi jθ P j θi jθ cv. 4.5 By Lemma 4.1, the solution to the first equation of this system is T j f jv, while the second equation together with the boundary condition 2.3 yields i j θ h jt j, V σ j θ h j f j V, V σ j θ. Substituting in the third equation of 4.5, we get P j θσ j θh j f j V, V cv and since V, we can divide the above equation by V and write it as F V with F V N j h j f j V, V c. By Lemma 4.2, the function V F V is continuous and strictly decreasing for V, and since lim V h jf j V, V by H3, then lim F V V N j lim h jf j V, V c c <. V 11
This implies that F V has no positive solutions if F, while it has exactly one positive solution if F >. Since h j f j, ψ j h j λj,, then F h j λj N j, c h j λj N j, 1 1, R so the steady state E exists and is positive if and only if R > 1. For the antibody-immune infected steady state Ê, we assume  and ˆV, so we obtain the system λ j ˆTj h j ˆT j, ˆV dî j θ δ j θî j θ P j θî j θ c ˆV qâ ˆV k ˆV h +  b. 4.6 The last equation of this system implies ˆV b k h+â. From the first and second equations, we have ˆT j f j ˆV and î j θ h j ˆT j, ˆV σ j θ h j fj b k h + Â, b k h +  σ j θ. Substituting in the third equation of 4.6, we get P j θσ j θh j fj b k h + Â, b k h +  c ˆV qâ ˆV and since ˆV, this is equivalent to GÂ, where GA N j h j fj b k h + A, b k h + A c qa. By Lemma 4.2 we know that h j f j V, V decreases with respect to V, so A GA is a continuous and decreasing function for A. Since V b k h + A tends to infinity as A and lim h jf j V, V, then V lim GA A N j lim h jf j V, V + lim c qa. V A This implies that GA has no positive solutions if G, while it has exactly one positive solution if G >. We have G N j h j fj bh k, bh k c N j h j fj bh k, bh k so the steady state Ê exists and is positive if and only if R > 1. 1 1 R, The above theorem implies that system 2.1 has only one steady state when R 1, two steady states when R 1 < R and three when 1 < R. We will now analyse their local stability by means of the characteristic equation of the system. 12
Theorem 4.4. The infection-free steady state E of system 2.1 is locally asymptotically stable if R < 1 and it is unstable if R > 1. Proof. Using the expression in 2.5 for i j θ, t, we have P j θi j θ, t t t P j θh j Tj t θ, V t θ σ j θ + k j θh j Tj t θ, V t θ + g j t, where k j θ P j θσ j θ and g j t t P j θi j θ t σ jθ σ j θ t 2.1 as the following system: t P j θi σ j θ jθ t σ j θ t. Thus, we can rewrite dt j t dv t dat λ j T j t h j T j t, V t [ t k j θh j Tj t θ, V t θ ] + g j t cv t qatv t katv t h + At bat. The linearization of this system at E is dt j t T j t h j λj, V t dv t hj λj, t k j θv t θ cv t dat bat, and the characteristic equation is n s + hj λj s + c, k j θe sθ s + b. This equation has the negative roots s b and s, j 1,..., n. The rest of its roots are the solutions of hj λj φs : s + c, k j θe sθ. 4.7 Suppose that 4.7 has a root s with non-negative real part. Then c s + c hj λj, k j θe sθ hj λj, k j θ e s θ, 13
but so k j θ e s θ c k j θ hj λj, N j cr. P j θσ j θ N j, This implies that 1 R. Therefore, when R < 1 all solutions of 4.7 have negative real part and thus, in such case, E is locally asymptotically stable. Otherwise, if R > 1 we have φ c hj λj, k j θ c1 R <, while φs as s. Then φs has at least one positive real root. Therefore, E is unstable when R > 1. We will now study the stability of the immune-free infected steady state E. For that, we define the reproductive number of antibody response R AN as R AN kv bh. We will also need the following result, which shows the relation between R AN and the viral reproduction number R. Lemma 4.5. R < 1 R AN < 1, R 1 R AN 1, and R > 1 R AN > 1. Proof. Suppose that R AN < 1. Then V < bh k. Recall that V satisfies F V N j h j fj V, V c. Since F is strictly decreasing, the condition V < bh k that is, N j h j fj bh k, bh k c <, cr 1 <. implies that F bh k < F V, i.e., Thus, R < 1. The proof that R AN 1 R 1 and that R AN > 1 R > 1 is similar, from which the result follows. Theorem 4.6. The immune-free infected steady state E is unstable when R > 1. 14
Proof. The linearization of system 2.1 at E is dt j t + h j Tj, V T j t h j T j T j, V V t dv t t hj Tj, V k j θt j t θ T j t hj + T j, V k j θv t θ cv t dat kv h b At. Thus, one of the roots of the characteristic equation is s 1 kv h b. If R > 1, then R AN kv bh > 1, which implies that s 1 is a positive root. Hence the theorem. 5 Global stability In order to establish the global stability of the equilibrium E when R 1, we will make use of LaSalle s invariance principle and a Lyapunov functional, similar to those used in [4] for a multi-target cell model with general incidence. Theorem 5.1. If R 1, then the infection-free steady state E is globally asymptotically stable. Proof. Assume R 1. For j 1,..., n, define the function γ j x x P j θe θ x δ jτ dτ. 5.1 We will use the Lyapunov functional W W T 1, i 1,,, T n, i n,, V, A given by W Tj N j T j Tj h j Tj lim, V Tj V h j η j, V dη j + γ j θi j θ, t + V. Then W is non-negative and W E. The time derivative of W along the solutions of 2.1 is dw h j Tj N j 1 lim, V λ j T j h j T j, V V h j T j, V ij θ, t γ j θ + δ j θi j θ, t θ + P j θi j θ, t cv qav h j Tj N j 1 lim, V V h j T j, V T j T j h j Tj N j 1 lim, V h j T j, V V h j T j, V 15
cv qav. Since T j λ j Thus, we have dw γ j θ di j θ, t h j T 1 lim j,v V h j T j,v γj θδ j θ P j θ i j θ, t h j Tj, V for all V, so 1 lim,v V h j T j,v. Tj, then h jt j, V h j Tj Tj T j and also qav, so N j 1 lim V h j T j, V h j T j, V h j T j, V γj θδ j θ P j θ i j θ, t cv N j h j T j, V + γ j θi j θ, t θ θ h j Tj N j h j T j, V lim, V V h j T j, V i j θ, tγ jθ γj θδ j θ P θ i j θ, t cv, γ j θ di j θ, t where we have used integration by parts to expand γ j θ di j θ, t. By 5.1 we know that γ j satisfies γ j x δ jxγ j x P j x, so i j θ, tγ jθ γj θδ j θ P j θ i j θ, t, and since γ j N j, then N j h j T j, V γ j θi j θ, t θ. Thus, dw N j h j T j, V + γ j θi j θ, t h j Tj N j h j T j, V lim, V V h j T j, V cv θ θ h j Tj N j h j T j, V lim, V V h j T j, V γ j θi j θ, t cv θ cv N j h j T j, V h j Tj lim, V cv V h j T j, V 1 γ j θi j θ, t. θ Since V and h j T j, V is non-increasing with respect to V, then We also have h j T j, V h j T j, h j T j, V h j T j, lim h j V V T j,. h j Tj lim, V V h j T j, V lim V h j T j,v h jt j, V h j T j,v h jt j, V 16 h j λj, h j T j,
so N j h j T j, V cv h j Tj lim, V V h j T j, V 1 N j c h jt j, V N j c R 1. h j h j λj, 1 λj, h j T j, 1 Hence we have dw and the equality holds if and only if T j Tj, i j, V. Thus, the largest positively invariant subset of the state space where dw is the set S { T 1, i 1,..., T n, i n, V, A Ω T j T j, i j, V }. By LaSalle s invariance principle, this implies that all solutions with positive initial conditions approach S as t. On the other hand, any solution of 2.1 contained in S satisfies da ba and thus, At as t, that is, the solution approaches E. Therefore, we conclude that E is globally asymptotically stable in the state space Ω. We will now use a Lyapunov functional that is constructed as a combination of those used in [3] and [4] to prove the global stability of E. Theorem 5.2. If R < 1 < R, then the immune-free infected steady state E is globally asymptotically stable. Proof. Consider the Lyapunov functional W 1 + V H Tj N j T j Tj T j h j T j, V h j ξ j, V dξ j + γ j θi jθh i j θ, t i j θ V V + qh k A, 5.2 where Hx x 1 ln x and γ j is defined by 5.1. Then this functional is well-defined based on the uniform persistence of the system. The function H is non-negative and equals zero only at 1, so it is clear that W 1 has a global minimum at the equilibrium E. The derivative of W 1 along the solutions of 2.1 is dw 1 N j 1 h jtj, V λj h j T j, V T j h j T j, V + γ j θ 1 i j θ ij θ, t i j θ, t t + 1 V V + qh kav k h + A ba. P j θi j θ, t cv qav 17
By 4.3 we have 1 V V cv cv V n N jh j Tj, V V V. Also, notice that λ j Tj + h jtj, V and γ j θ 1 i j θ ij θ, t i j θ, t t γ j θi jθ t H i j θ, t i j θ γ j θi i j θ, t θ i j θ, t jθh i j θ + H i θ j θ γ j i i j, t jh i j γ j θi i j θ, t jθh i j θ θ N j h j Tj, V h j T j, V H h j Tj, V γ j θi i j θ, t jθh i j θ P j θi i j θ, t jθh i j θ, γ θi jθ + γ j θ di θ θ P j θi jθh i j θ, t i j θ since γ θi j θ + γ jθ di θ P j θi j θ, γ j N j, and i j, t h j T j, V. Then we obtain dw 1 N j 1 h jtj, V dj h j T j, V Tj T j + h j Tj, V h j T j, V N [ + N j h j Tj, V h j T j, V H h j Tj, V γ j θi i j θ, t jθh i j θ θ ] P j θi i j θ, t jθh i j θ + P j θi j θ, t V P j θi j θ, t V + N j h j Tj, V V N j Tj T j + N j h j Tj, V 1 h jtj, V h j T j, V ln h jt j, V h j T j, V + ln h jt j, V + + qav P j θi jθ h h + A 1 N j h j Tj, V V qav + qav + qhav h + A qhba k γ j θi i j θ, t jθh i j θ θ h jtj, V h j T j, V h jt j, V h j Tj, V + h jt j, V h j T j, V + h jt j, V h j Tj, V h j T j, V + 1 V V V i j θ, t V i j θ + 1 + ln i jθ, t i j θ + qa k kv bh. 18
By H2, we know that h jtj,v h j T j,v < 1 when T j < T j and h jtj,v h j T j,v > 1 when T j > T j, so Tj T j 1 h jtj,v h j T j,v for all T j. It is clear that γ j θi j θh ij θ,t θ i j θ. We also have qav h h+a 1 since h h + A, and qa k kv bh qabh k R AN 1 since the hypothesis R < 1 < R implies that R AN < 1. Thus, dw 1 N j h j Tj, V h jtj, V h j T j, V h jt j, V h j Tj, V + h jt j, V h j T j, V + h jt j, V h j Tj, V ln h jt j, V h j T j, V + ln h jtj, V h j T j, V + 1 V V P j θi V i j θ, t jθ V i j θ 1 ln i jθ, t i j θ. [ ] N j h j Tj, V hj T j, V V H h j T j, V H V P j θi V i j θ, t jθh V i j θ. By H2 and H3, we have that V V V implies V h jt j,v h j T j,v 1, while V V implies 1 h jt j,v h j T j,v V V. Since Hx x 1 ln x is decreasing for x, 1 and increasing for x 1,, then H H V V for all V. Also, hj T j,v h j T j,v non-negativity of H implies that P j θi j θh V i j θ,t V i j θ. Thus, the above argument shows that dw 1 and that the equality holds if and only if T j Tj, i jθ, t i j θ, V V, and A. Therefore, by LaSalle s invariance principle we conclude that E is globally asymptotically stable. Similarly, inspired in [3] and [4], we make a modification to the Lyapunov functional 5.2 in order to prove that Ê is globally asymptotically stable whenever it exists. Theorem 5.3. If R > 1, then the antibody-immune infected steady state Ê is globally asymptotically stable. Proof. Consider the Lyapunov functional W 2 N j T j ˆT j + ˆV H Tj ˆT j h j ˆT j, ˆV h j ξ j, ˆV dξ j + VˆV q A h + ττ  dτ, k  τ + i j θ, t γ j θî j θh î j θ A h+ττ  τ where H and γ j are the same as for 5.2. It can be seen that the term q k dτ is positive for A  and equals zero for A Â, so W 2 has a global minimum at Ê. We calculate the derivative of W 2 along the solutions of 2.1 as follows: dw 2 N j 1 h j ˆT j, ˆV λj h j T j, ˆV T j h j T j, V 19
+ + 1 ˆV V γ j θ 1 îjθ i j θ, t + q kv k h + AA Â h + A b. i j θ, t t P j θi j θ, t cv qav Following the same steps as in the proof of Theorem 5.2, we obtain dw 2 N j ˆT j T j 1 h j ˆT j, ˆV i j θ, t h j T j, ˆV γ j θî j θh î j θ θ + N j h j ˆT j, ˆV h j ˆT j, ˆV h j T j, ˆV h jt j, V h j ˆT j, ˆV + h jt j, V h j T j, ˆV + h jt j, V h j ˆT j, ˆV ln h jt j, V h j T j, ˆV + ln h j ˆT j, ˆV h j T j, ˆV + 1 VˆV + P j θî j θ ˆV i j θ, t V î j θ + 1 + ln i jθ, t î j θ qav + qa ˆV + q kv k h + AA Â h + A b. We have N j ˆT j T j ˆV bh + Â/k, we also have so dw 2 1 h j ˆT j, ˆV h j T j, ˆV qav + qa ˆV + q k h + AA Â kv h + A b and γ j θi j θh ij θ,t i j θ θ. qav + qa ˆV + qv A Â qb qb h + AA + h + AÂ k [ k ] qb A Â h + A q ˆV k qb k A Â2, Since N j h j ˆT j, ˆV h j ˆT j, ˆV h j T j, ˆV h jt j, V h j ˆT j, ˆV + h jt j, V h j T j, ˆV + h jt j, V h j ˆT j, ˆV ln h jt j, V h j T j, ˆV + ln h j ˆT j, ˆV h j T j, ˆV + 1 VˆV ˆV ij θ, t P j θî j θ V î j θ 1 ln i jθ, t î j θ [ ] N j h j ˆT j, ˆV h j T j, V H h j T j, ˆV H VˆV ˆV ij θ, t P j θî j θh V î j θ. Thus, dw 2 and the equality holds if and only if T j ˆT j, i j θ, t î j θ, V ˆV, and A Â. Therefore, by LaSalle s invariance principle we conclude that Ê is globally asymptotically stable. 2
6 Related models In this section, we consider a special case in which the model 2.1 can be reduced to a system of ordinary differential equations with one time delay. For that, we will use a particular form of incidence rate, which is given by the saturated incidence function h j T j, V β jt j tv t 1+α j V t. We will assume that it takes time τ for virus to enter into the target cell and that there is an intracellular delay ω that describes the time required for an infected cell to produce virus, so the death rate of infected cells and viral production rate are given by the functions { if θ < τ, δ j θ 6.1 if θ τ and P j θ δ j { if θ < ω, P j if θ ω. 6.2 We also assume that the initial condition i j θ, i j θ satisfies lim θ i j θ, which means that the initial number of infected cells tends to as the infection age tends to infinity. If I j t i j θ, t is the number of infected cells in the j-th class at time t, then but by 2.5 we have di j t i j θ, t i j θ, t t ij θ, t θ θ θ + δ j θi j θ, t δ j θi j θ, t, so lim i jθ, t lim i σ j θ θ θ jθ t σ j θ t lim θ i jθ te δ j t, di j t i j, t β jt j tv t 1 + α j V t δ j θi j θ, t τ i j θ, t τ δ j i j θ, t. Using the function 6.2 we obtain P j θi j θ, t P j P j ω ω i j θ, t i j θ ω, t ωe δ jθω. Thus, with these assumptions, the model can be reformulated equivalently as the following 21
system: dt j t di j t dv t dat λ j T j t β jt j tv t 1 + α j V t β jt j tv t τ 1 + α j V t i j θ, t δj i j θ, t τ i j θ ω, t ωe δjθω cv t qatv t P j ω katv t h + At bat. 6.3 In this model we have two intracellular time delays, τ and ω. If we consider the case when the death rate of infected cells is constant, i.e., δ j θ δj for all θ, then the above system can be further simplified to the following system of delay differential equations: dt j t di j t dv t dat λ j T j t β jt j tv t 1 + α j V t β jt j tv t 1 + α j V t δ j I j t Pj e δ j ω I j t ω cv t qatv t katv t h + At bat, 6.4 which contains no integral terms and only one time delay. In this case, we have σ j θ e δ j θ and N j Pj e δ j ω /δj, so the basic reproduction numbers of virus and antibodies for 6.4 can be calculated as R 1 c Pj e δ j ω λ j β j δj, R 1 c 7 Numerical examples δ j kpj e δ j ω λ j β j kdj + α j + β j bh. We will now perform some numerical simulations to illustrate the analytical results obtained for the dynamics of model 2.1. 7.1 Model with one class of target cells We first simulate the case n 1 for system 2.1 with the particular forms of δ 1 θ and P 1 θ used in Section 6 to obtain the simplified system 6.4. Hence, we are considering a model with intracellular delay ω in viral production, constant death rate δ for infected cells, and saturated incidence rate. We use the parameters shown in Table 1, which are similar to those used in [7, 21, 22] for models of HIV infection with only one population of target cells. We consider three different values for the infection rate β, given in virions/day, in order to obtain the three possible scenarios for the dynamics of the model. For each case, we show the graphs of solutions with respect to time for several initial conditions, which are given in cells/ml for T t and It, in virions/ml for V t, and in µg for At. 22
Parameter Value Units λ 1 46 cells ml 1 day 1 d 1.46 day 1 α 1.5 δ1.1 day 1 P1 11.459 virions day 1 cell 1 ω.5 days c.25 day 1 q.3 µg 1 day 1 k.15 virion 1 day 1 h.2 b 2.9 day 1 Table 1: Parameters used for numerical simulations of system 6.4 with one class of target cells. Figure 1: Solutions of system 6.4 when R < R < 1 and E is globally asymptotically stable. 23
Figure 2: Solutions of system 6.4 when R < 1 < R and E is globally asymptotically stable. Figure 3: Solutions of system 6.4 when 1 < R < R and stable. Ê is globally asymptotically 24
We first consider the case when β 1 1 8, obtaining R.454 and R.1547, which are less than 1. According to Theorem 5.1, all solutions approach the uninfected equilibrium E as t, as we can see in Figure 1. For β 5 1 8, we have R 2.7238, R.927, so according to Theorem 5.2 all solutions converge to the immune-free infected equilibrium E, as we can see in Figure 2. Lastly, if we take β 5 1 7, then we have R 22.698, R 7.6287, so according to Theorem 5.3 all solutions converge to the immune infected equilibrium Ê, as shown in Figure 3. 7.2 Model with two classes of target cells We will now assume that there are two different populations of target cells: CD4 + T cells denoted by the subscript j 1 and macrophages j 2. For the CD4 + T cells, we consider a viral production kernel P 1 θ similar to that used in [2], which is given by { if θ < θ 1, P 1 θ 1 e rθ θ 1, if θ θ 1, P 1 where P1 6.421 13 virions day 1 cell 1, while r 1 is a saturation parameter and θ 1.25 days is the age at which reverse transcription is completed. We use a death rate for infected cells δ 1 θ of the form 6.1 with d 1.46 day 1, δ1 1.5 day 1 and a saturated infection rate h 1 T 1, V β 1T 1 V 1+α 1 V with α 1.5, β 1 2.4 1 8 virions day 1. For the macrophages, we use a constant death rate δ 2 θ 1/14.1 day 1 and we assume that viral particles at produced at a rate P 2 θ.1 + 1 3.28θ [23]. We assume that the infection rate is given by the function h 2 T 2, V 1.19T 2V [24]. The rest 1 6 +V of the parameters used are given in Table 2 and are based on [2, 7, 23]. Figure 4 shows the solutions for each variable of system 2.1 with the above parameters, using the initial conditions T 1 5, T 2 1 6, i 1 θ, i 2 θ, 1.194 1 4 e 1θ, V.1, A.1. In this case, the basic reproduction numbers are R 1 17.82 for the population of CD4 + T cells and R 2.2963 for macrophages, so the total basic reproduction number of virus is R R 1 + R 2 18.1163. From this, we can see that the relative contribution of T cells to infection is considerably stronger than that of macrophages. The viral reproduction number with antibody response is R 6.3637. Since R > 1 and R > 1, all solutions converge to a positive steady state, which corresponds to the chronic phase of HIV, as can clearly be seen in the simulations. 8 Conclusions In this paper, we have studied an age-structured within-host viral dynamic model that includes multiple populations of target cells. Our model incorporates saturated antibody immune response and a general non-linear incidence rate, so it can be viewed as a generalization of several models that have been studied previously in the literature. This kind of models can be used to evaluate the relative contribution of viral production from different compartments of target cells, such as CD4 + cells, macrophages and dendritic cells, and this can help us improve our understanding of the dynamics of HIV infection. We have extended the results published in [3], where the authors studied a similar model with a saturated incidence that did not include multiple classes of target cells. We proved that the global dynamics of the model is completely determined by the basic 25
Parameter Value Units λ 1 1 4 cells ml 1 day 1 d 1.1 day 1 λ 2 864 cells ml 1 day 1 d 2.24 day 1 c 23 day 1 q.3 µg 1 day 1 k.15 virion 1 day 1 h.2 b 2.9 day 1 Table 2: Parameters used for numerical simulations of system 2.1 with two classes of target cells. Figure 4: Time evolution of the solution of system 2.1 with two classes of target cells. 26
reproduction number and the reproductive number of antibody response. If R is less than one, the only steady state of the model is the infection-free equilibrium and the infection is predicted to die out. When R < 1 < R, an infected equilibrium E appears which is globally asymptotically stable, so infection becomes chronic but with no persistent antibody immune response. Lastly, if R > 1, the antibody-immune infected equilibrium Ê appears and all solutions converge to it; this case corresponds to the scenario when the infection becomes chronic with a persistent antibody immune response. References 1. Nelson, P. W., Gilchrist, M. A., Coombs, D., Hyman, J. M. & Perelson, A. S. An age-structured model of HIV infection that allows for variations in the production rate of viral particles and the death rate of productively infected cells. Math. Biosci. Eng 1, 267 288 24. 2. Rong, L., Feng, Z. & Perelson, A. S. Mathematical analysis of age-structured HIV-1 dynamics with combination antiretroviral therapy. SIAM Journal on Applied Mathematics 67, 731 756 27. 3. Duan, X. & Yuan, S. Global dynamics of an age-structured virus model with saturation effects. Mathematical Methods in the Applied Sciences 4, 1851 1864 217. 4. Wang, X., Lou, Y. & Song, X. Age-Structured Within-Host HIV Dynamics with Multiple Target Cells. Studies in Applied Mathematics 138, 43 76 217. 5. Wang, S., Wu, J. & Rong, L. A note on the global properties of an age-structured viral dynamic model with multiple target cell populations. Mathematical biosciences and engineering: MBE 14, 85 82 217. 6. Yang, Y., Ruan, S. & Xiao, D. Global stability of an age-structured virus dynamics model with Beddington DeAngelis infection function. Math. Biosci. Eng 12, 859 877 215. 7. Tian, X. & Wang, J. Stability analysis for viral infection model with multitarget cells, Beddington-DeAngelis functional response, and humoral immunity. Discrete Dynamics in Nature and Society 215 215. 8. Wang, J., Tian, X. & Wang, X. Stability analysis for delayed viral infection model with multitarget cells and general incidence rate. International Journal of Biomathematics 9, 1657 216. 9. Wang, J., Lang, J. & Li, F. Constructing Lyapunov functionals for a delayed viral infection model with multitarget cells, nonlinear incidence rate, state-dependent removal rate. J Nonlinear Sci Appl 9, 524 536 216. 1. Elaiw, A. M. Global properties of a class of HIV models. Nonlinear Analysis: Real World Applications 11, 2253 2263 21. 11. Elaiw, A. M. Global properties of a class of virus infection models with multitarget cells. Nonlinear Dynamics 69, 423 435 212. 12. Elaiw, A. M. & Azoz, S. A. Global properties of a class of HIV infection models with Beddington DeAngelis functional response. Mathematical Methods in the Applied Sciences 36, 383 394 213. 13. Wang, X., Chen, Y., Liu, S. & Song, X. A class of delayed virus dynamics models with multiple target cells. Computational and Applied Mathematics 32, 211 229 213. 27
14. Culshaw, R. V. & Ruan, S. A delay-differential equation model of HIV infection of CD4+ T-cells. Mathematical biosciences 165, 27 39 2. 15. Li, D. & Ma, W. Asymptotic properties of a HIV-1 infection model with time delay. Journal of Mathematical Analysis and Applications 335, 683 691 27. 16. Georgescu, P. & Hsieh, Y.-H. Global stability for a virus dynamics model with nonlinear incidence of infection and removal. SIAM Journal on Applied Mathematics 67, 337 353 26. 17. Hale, J. K. Asymptotic behavior of dissipative systems 25 American Mathematical Soc., 21. 18. Magal, P. & Thieme, H. R. Eventual compactness for semiflows generated by nonlinear age-structured models. Commun. Pure Appl. Anal 3, 695 727 24. 19. Hale, J. K. & Waltman, P. Persistence in infinite-dimensional systems. SIAM Journal on Mathematical Analysis 2, 388 395 1989. 2. Demasse, R. D. & Ducrot, A. An age-structured within-host model for multistrain malaria infections. SIAM Journal on Applied Mathematics 73, 572 593 213. 21. Pawelek, K. A., Liu, S., Pahlevani, F. & Rong, L. A model of HIV-1 infection with two time delays: mathematical analysis and comparison with patient data. Mathematical Biosciences 235, 98 19 212. 22. Wang, J., Pang, J., Kuniya, T. & Enatsu, Y. Global threshold dynamics in a fivedimensional virus model with cell-mediated, humoral immune responses and distributed delays. Applied Mathematics and Computation 241, 298 316 214. 23. Duffin, R. P. & Tullis, R. H. Mathematical models of the complete course of HIV infection and AIDS. Computational and Mathematical Methods in Medicine 4, 215 221 22. 24. Wasserstein-Robbins, F. A mathematical model of HIV infection: simulating T4, T8, macrophages, antibody, and virus via specific anti-hiv response in the presence of adaptation and tropism. Bulletin of mathematical biology 72, 128 1253 21. 28