Research Article | Open Access
Xiaoting Fan, Yi Song, Wencai Zhao, "Modeling Cell-to-Cell Spread of HIV-1 with Nonlocal Infections", Complexity, vol. 2018, Article ID 2139290, 10 pages, 2018. https://doi.org/10.1155/2018/2139290
Modeling Cell-to-Cell Spread of HIV-1 with Nonlocal Infections
This paper is devoted to develop a nonlocal and time-delayed reaction-diffusion model for HIV infection within host cell-to-cell viral transmissions. In a bounded spatial domain, we study threshold dynamics in terms of basic reproduction number for the heterogeneous model. Our results show that if , the infection-free steady state is globally attractive, implying infection becomes extinct, while if , virus will persist in the host environment.
1. Introduction and Model Derivation
Infectious diseases have been threatening human health; infectious diseases such as smallpox, cholera, and acquired immune deficiency syndrome (AIDS) have brought great disaster to the national economy of a country and people’s livelihood. In order to control the spread of infectious diseases, researchers have proposed a great deal of mathematical models to study the evolution behavior of infectious diseases from different views [1–7]. Mathematical models have been confirmed to be an effective and valuable approach to understand the dynamical behavior of infectious diseases [8–11]. Many diseases are caused by a virus, for example, human immunodeficiency virus (HIV) can destroy the immune cells, reduce human immunity, and ultimately lead to AIDS. A lot of mathematical models have been built to explain the virus infection process from different views [12–15]. In recent years, the virus infection dynamical models incorporating spatial dispersion [16, 17] have received widely attentions. For example, under spatial domain equipped with suitable boundary conditions, K. Wang and W. Wang  established the existence of travelling wave solutions via the geometric singular perturbation method and the spatial domain is assumed to be one dimensional, that is, . And then, Brauner et al.  extended the works  to a two-dimensional square domain with periodic boundary conditions by introducing a new parameter which is the largest eigenvalue of some Sturm-Liouville problem and taking the heterogeneous reproductive ratio into account. Recently, Wang et al.  noticed that a realistic domain should be bounded but is typically not a square, under no-flux boundary condition (homogenous Neumann boundary condition) in a general bounded domain with smooth boundary , and threshold dynamics of the model were investigated by appealing to the theory of uniform persistence and the comparison theorem. Some further developments have been performed on the virus infection dynamical models with diffusion term (see, for example, Hattaf and Yousfi , Wang and Ma [22–24], Lai and Zou , Wang and Xu , and Wang et al. ).
It is widely known that virus-to-cell transmission and direct cell-to-cell transmission are two predominant infection modes in within-host environment (see, for example, [28, 29]). The virus-to-cell transmission has iterative process: the binding of virus to a receptor on the surface of CD4 T cells, the fusion of virus with host cells and then the release of genetic materials, the transcription of genetic materials in an infected cell, the assembling of virus inside the infected cells, and budding through the membrane of the infected cells . The cell-to-cell transmission process is due to that a large number of viruses can be transferred from infected cells to uninfected cells through the formation of virally induced structures termed virological synapses . In order to examine the effects of both diffusion and spatial heterogeneity, Wang et al.  proposed the following virus infection dynamical model incorporating cell-to-cell transmission under the homogeneous Neumann boundary condition:
In model (1), is the production rate of newly produced uninfected cells. , , and represent the death rates of uninfected cells, infected cells, and free viruses, respectively. denotes the transmission coefficient for virus-to-cell infection. represents the transmission coefficient for cell-to-cell infection. denotes the production rate of virus due to the lysis of infected cells. is the diffusion coefficient and is Laplacian. For model (1), the virus dynamics are fully determined by an important parameter which is called basic reproduction number , in the following sense: if , the infection-free steady state is globally asymptotically stable, while if , the model is uniformly persistent and the infection steady state is globally asymptotically stable.
The mobility of cells in the latent period will result in a delay term with spatial averaging on the spatial habitat . To formulate this process with the latency properly, we introduce an infection age variable . The infected cells can be divided into two epidemiological categories: latently infected cells () and actively infected cells (). Let represent the density of cells with infection age at time and habitat . From the biological perspective, we assume that all populations remain confined to the region for all time and subject to no-flux boundary condition for : where is the outward normal to the smooth boundary . We adopt the standard model on describing age-structured population with spatial diffusion . Then
In (3), is the gradient of with respect to the spatial variable and represents the divergence of . represents the diffusion rate of cells at age and habitat . is the natural death rate which is independent of the infection age . Let be the average incubation period. Then
In order to make the model mathematically tractable yet without losing the main characteristics, we make some assumptions for function as follows:
Integrating both sides of (3) from to and from to , it then follows that and
For biological reasons, we assume that (see, for example, ), since the recruitment of newly infected cells is divided into two parts: the contact of uninfected cells and virus and the contact of uninfected cells and infected cells. Thus, adopting Beddington-DeAngelis functional response which was firstly proposed by [36, 37] leads to the following condition: where and are the measure of virus interference during infection and and determine how fast the infection rate approaches its saturation value (see, for example, [36, 37]).
In the following, we determine by the method of characteristics. Let with . For , it then follows that
Regarding as a parameter and solving the above equation, we obtain that where is the fundamental solution associated with the partial differential operator , which satisfies , . , and , (see, , Chapter 1).
Since , , we get the formula for for :
It then follows that and
With all these assumptions, the HIV infection dynamics modeling cell-to-cell with nonlocal infections can be described by the following nonlocal and time-delayed reaction-diffusion equations:
Since is decoupled from other equations for model (14), we consider the following dynamical model: for , . We consider a closed environment in the sense that the fluxes for each of these subpopulations are aero, and hence we propose no-flux condition on the boundary: and initial conditions
We now briefly outline the plan of this paper. In Section 2, we establish basic reproduction number . In Section 3, we study threshold dynamics in terms of basic reproduction number for the heterogeneous model. At last, we give discussions and biological implications in Section 4.
2. Basic Reproduction Number
Let be the Banach space with the supremum norm by , , where denotes the Euclidean norm in and denotes the transpose of the vector, and . We easily get that is a closed cone of .
Define with the norm and , then is an ordered Banach space. For any given continuous function for , define by , . Let and . Suppose that , , and are the strongly continuous semigroups associated with , , and subject to homogeneous Neumann boundary conditions, respectively, that is, for any , , and , . and are the Green functions associated with , , and , respectively, subject to homogenous Neumann boundary conditions. According to Section 7.1 and Corollary 7.2.3 in , it follows that is strongly positive and compact for each . We can also obtain that , , is a strongly continuous semigroup. Let be the generator of , respectively . Then, we get that is a semigroup generated by the operator defined on .
In view of Lemma 2.2 in , model (15) exists a unique infection-free steady state Linearizing model (15) at the infection-free steady state , we get the following linearized system: satisfying the following boundary conditions
We consider the following model:
Substituting and into and of model (21) which is called the infectious compartments results in the following eigenvalue problem:
We further consider the nonlocal eigenvalue problem associated with
Lemma 2.1. The eigenvalue problem (23) has a principal eigenvalue with a strictly positive eigenfunction, and for any , has the same sign as .
In order to compute basic reproduction number , we assume that both infected cells and virus are near the infection-free steady state and introduce initial infected cells and virus at time , where the distribution of initial infected cells and virus is described by . From model (15) and as time evolves, those distributions can reach at time . Thus, the total distribution of new infected cells is
Similarly, the total distribution of new virus can be described as
It then follows that where is the next infection operator (see, for example, [41–44]). It maps the initial distribution of infected cells and virus to the total distribution of new infected cells and virus produced during the infection period.
In view of the general results developed in , we have the following observation.
Lemma 2.2. has the same sign as .
3. Threshold Dynamics in a Bounded Spatial Domain
In this section, we study threshold dynamics in terms of for model (15), assuming a bounded spatial domain . We assume that all the parameters , , , , , , , and are spatially dependent. Assume these functions are positive, continuous, and bounded in the domain . Let
For any and , we define by
Obviously, is locally Lipschitz. Then, model (15) can be rewritten as the following abstract functional differential equation: where .
Proof 1. For any and for sufficiently small , we obtain that Thus, . This implies that for all From Corollary 4 in  (see, also , Corollary 8.1.3), we obtain the conclusion stated in Lemma 3.1. The proof is completed.
We are now in the position to study the well-posedness of model (15) in the sense of the following theorem.
Theorem 3.1. For any , model (15) has a unique solution on with . Moreover, the solution semiflow has a compact global attractor for .
Proof 2. It is easy to see that model (15) defines a semiflow by
For any fixed , from the first equation of model (15), we obtain that
From the comparison principle, we easily get that there exists such that for . From the second equation of model (15), we have that
There exists such that for . From the third equation of model (15), we get
There also exists such that for .
Consequently, the existence of solutions of model (15) claimed in Lemma 3.1 is indeed global (i.e., ). The solution semiflow is point dissipative. Moreover, by Theorem 2.2.6 in , we have that is compact for any . Thus, from Theorem 3.4.8 in , we know that has a compact global attractor in for . We complete the proof.
The following results will play an important role in establishing uniform persistence of model (15).
Lemma 3.2. Let be the solution of model (15) with , then we have the following:
(i) If there exists some such that and , then for all , .
(ii) It holds that for any , , and uniformly for .
Proof 3. From the second and third equations of model (15), it follows that
with the following boundary conditions
If and for some , it follows from the strong maximum principle (see, for example,  p. 172, Theorem 4) and the Hopf boundary lemma (see, for example,  p. 170, Theorem 3) that , , and for , , that is, the conclusion holds.
From the first equation of model (15), we get