Global dynamics of a mathematical model for HTLV-I infection of CD4 T cells with delayed CTL response

Human T-cell leukaemia virus type I (HTLV-I) preferentially infects the CD4⁺ T cells. The HTLV-I infection causes a strong HTLV-I specific immune response from CD8⁺ cytotoxic T cells (CTLs). The persistent cytotoxicity of the CTL is believed to contribute to the development of a progressive neurologic disease, HTLV-I associated myelopathy/tropical spastic paraparesis (HAM/TSP). We investigate the global dynamics of a mathematical model for the CTL response to HTLV-I infection in vivo. To account for a series of immunological events leading to the CTL response, we incorporate a time delay in the response term. Our mathematical analysis establishes that the global dynamics are determined by two threshold parameters R₀ and R₁, basic reproduction numbers for viral infection and for CTL response, respectively. If R₀≤1, the infection-free equilibrium P₀ is globally asymptotically stable, and the HTLV-I viruses are cleared. If R₁≤1<R₀, the asymptomatic-carrier equilibrium P₁ is globally asymptotically stable, and the HTLV-I infection becomes chronic but with no persistent CTL response. If R₁>1, a unique HAM/TSP equilibrium P₂ exists, at which the HTLV-I infection is chronic with a persistent CTL response. We show that the time delay can destabilize the HAM/TSP equilibrium, leading to Hopf bifurcations and stable periodic oscillations. Implications of our results to the pathogenesis of HTLV-I infection and HAM/TSP development are discussed.     

Dynamics of an HIV‐1 virus model with both virus‐to‐cell and cell‐to‐cell transmissions,general incidence rate,intracellular delay,and CTL immune responses

Direct cell‐to‐cell transmission of HIV‐1 is a more efficient means of virus infection than virus‐to‐cell transmission. In this paper, we incorporate both these transmissions into an HIV‐1 virus model with nonlinear general incidence rate, intracellular delay, and cytotoxic T lymphocyte (CTL) immune responses. This model admits three types of equilibria: infection‐free equilibrium, CTL‐inactivated equilibrium, and CTL‐activated equilibrium. By using Lyapunov functionals and LaSalle invariance principle, it is verified that global threshold dynamics of the model can be explicitly described by the basic reproduction numbers.     

Global properties of a class of HIV infection models with Beddington–DeAngelis functional response

In this paper, the global properties of a class of human immunodeficiency virus (HIV) models with Beddington–DeAngelis functional response are investigated. Lyapunov functions are constructed to establish the global asymptotic stability of the uninfected and infected steady states of three HIV infection models. The first model considers the interaction process of the HIV and the CD4 ⁺ T cells and takes into account the latently and actively infected cells. The second model describes two co‐circulation populations of target cells, representing CD4 ⁺ T cells and macrophages. The third model is a two‐target‐cell model taking into account the latently and actively infected cells. We have proven that if the basic reproduction number R₀ is less than unity, then the uninfected steady state is globally asymptotically stable, and if R₀ > 1, then the infected steady state is globally asymptotically stable. Copyright © 2012 John Wiley & Sons, Ltd.     

Stability of a general adaptive immunity virus dynamics model with multistages of infected cells and two routes of infection

This paper studies an (n+4)-dimensional nonlinear virus dynamics model that characterizes the interactions of the viruses, susceptible host cells, n-stages of infected cells, B cells and cytotoxic T lymphocyte (CTL) cells. Both viral and cellular infections have been incorporated into the model. The infected-susceptible and virus-susceptible infection rates as well as the generation and removal rates of all compartments are described by general nonlinear functions. Five threshold parameters are computed, which insure 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 dynamics 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.     

Bifurcation analysis of the HIV‐1 within host model

In this paper, a bifurcation solution's analysis is proposed for an HIV‐1 within the host model around its chronic equilibrium point, this is carried out based on Lyapunov–Schmidt approach. It is shown that the coefficient b, which represents the healthy CD4⁺ T‐cells growth rate, is a bifurcation parameter; this means that the rate of multiplication of healthy cells can have serious effects on the qualitative dynamical properties and structural stability of the infection evolution dynamics. Copyright © 2015 John Wiley & Sons, Ltd.     

On global stability of an HIV pathogenesis model with cure rate

In this paper, applying both Lyapunov function techniques and monotone iterative techniques, we establish new sufficient conditions under which the infected equilibrium of an HIV pathogenesis model with cure rate is globally asymptotically stable. By giving an explicit expression for eventual lower bound of the concentration of susceptible CD4⁺ T cells, we establish an affirmative partial answer to the numerical simulations investigated in the recent paper [Liu, Wang, Hu and Ma, Global stability of an HIV pathogenesis model with cure rate, Nonlinear Analysis RWA (2011) 12 : 2947–2961]. Our monotone iterative techniques are applicable for the small and large growth rate in logistic functions for the proliferation rate of healthy and infected CD4⁺ T cells. Copyright © 2014 John Wiley & Sons, Ltd.     

Global stability of delay-distributed viral infection model with two modes of viral transmission and B-cell impairment

This paper formulates a virus dynamics model with impairment of B-cell functions. The model incorporates two modes of viral transmission: cell-free and cell-to-cell. The cell-free and cell-cell incidence rates are modeled by general functions. The model incorporates both, latently and actively, infected cells as well as three distributed time delays. Nonnegativity and boundedness properties of the solutions are proven to show the well-posedness of the model. The model admits two equilibria that are determined by the basic reproduction number R₀. The global stability of each equilibrium is proven by utilizing Lyapunov function and LaSalle's invariance principle. The theoretical results are illustrated by numerical simulations. The effect of impairment of B-cell functions and time delays on the virus dynamics are studied. We have shown that if the functions of B-cell is impaired, then the concentration of viruses is increased in the plasma. Moreover, we have observed that increasing the time delay will suppress the viral replication.     

A fractional-order model of HIV infection with drug therapy effect

In this paper, the fractional-order model that describes HIV infection of CD4⁺ T cells with therapy effect is given. Generalized Euler Method (GEM) is employed to get numerical solution of such problem. The fractional derivatives are described in the Caputo sense.     

Long term dynamics in a mathematical model of HIV-1 infection with delay in different variants of the basic drug therapy model

Infection with HIV-1, degrading the human immune system and recent advances of drug therapy to arrest HIV-1 infection, has generated considerable research interest in the area. Bonhoeffer et al. (1997) [1], introduced a population model representing long term dynamics of HIV infection in response to available drug therapies. We consider a similar type of approximate model incorporating time delay in the process of infection on the healthy T cells which, in turn, implies inclusion of a similar delay in the process of viral replication. The model is studied both analytically and numerically. We also include a similar delay in the killing rate of infected CD4⁺ T cells by Cytotoxic T-Lymphocyte (CTL) and in the stimulation of CTL and analyse two resulting models numerically.The models with no time delay present have two equilibria: one where there is no infection and a non-trivial equilibrium where the infection can persist. If there is no time delay then the non-trivial equilibrium is locally asymptotically stable. Both our analytical results (for the first model) and our numerical results (for all three models) indicate that introduction of a time delay can destabilize the non-trivial equilibrium. The numerical results indicate that such destabilization occurs at realistic time delays and that there is a threshold time delay beneath which the equilibrium with infection present is locally asymptotically stable and above which this equilibrium is unstable and exhibits oscillatory solutions of increasing amplitude.     

A delay differential equation model of HIV infection of CD4+ T-cells with cure rate

A delay differential equation as a mathematical model that described HIV infection of CD4⁺ T-cells is analyzed. When the constant death rate of infected but not yet virus-producing cells is equal to zero, the stability of the non-negative equilibria and the existence of Hopf bifurcation are investigated. A stability switch in the system due to variation of delay parameter has been observed, so is the phenomena of Hopf bifurcation and stable limit cycle. The estimation of the length of delay to preserve stability has been calculated. Further, when the constant death rate of infected but not yet virus-producing cells is not equal to zero, by using the geometric stability switch criterion in the delay differential system with delay dependent parameters, we present that stable equilibria become unstable as the time delay increases. Numerical simulations are carried out to explain the mathematical conclusions.     

Global and Bifurcation Analysis of an HIV Pathogenesis Model with Saturated Reverse Function

In this paper, an HIV dynamics model with the proliferation of CD4 T cells is proposed. The authors consider nonnegativity, boundedness, global asymptotic stability of the solutions and bifurcation properties of the steady states. It is proved that the virus is cleared from the host under some conditions if the basic reproduction number R_0 is less than unity. Meanwhile, the model exhibits the phenomenon of backward bifurcation. We also obtain one equilibrium is semi-stable by using center manifold theory. It is proved that the endemic equilibrium is globally asymptotically stable under some conditions if R_0 is greater than unity. It also is proved that the model undergoes Hopf bifurcation from the endemic equilibrium under some conditions. It is novelty that the model exhibits two famous bifurcations,backward bifurcation and Hopf bifurcation. The model is extended to incorporate the specific Cytotoxic T Lymphocytes(CTLs) immune response. Stabilities of equilibria and Hopf bifurcation are considered accordingly. In addition, some numerical simulations for justifying the theoretical analysis results are also given in paper.     

Stabilization of HIV/AIDS model by receding horizon control

This work concerns the stabilization of uninfected steady state of an ordinary differential equation system modeling the interaction of the HIV virus and the immune system of the human body. The control variable is the drug dose, which, in turn, affects the rate of infection of CD4⁺ T cells by HIV virus. The feedback controller is constructed by a variant of the receding horizon control (RHC) method. Simulation results are discussed.     

Analysis of the dynamics of a delayed HIV pathogenesis model

In this paper, considering full Logistic proliferation of CD4⁺ T cells, we study an HIV pathogenesis model with antiretroviral therapy and HIV replication time. We first analyze the existence and stability of the equilibrium, and then investigate the effect of the time delay on the stability of the infected steady state. Sufficient conditions are given to ensure that the infected steady state is asymptotically stable for all delay. Furthermore, we apply the Nyquist criterion to estimate the length of delay for which stability continues to hold, and investigate the existence of Hopf bifurcation by using a delay τ as a bifurcation parameter. Finally, numerical simulations are presented to illustrate the main results.     

Bifurcation analysis of a multidelayed HIV model in presence of immune response and understanding of in‐host viral dynamics

In this paper, we proposed a multidelayed in‐host HIV model to represent the interaction between human immunodeficiency virus and immune response. One delay was considered to incorporate the time required by the virus for various intracellular events to make a host cell productively infective, and the second delay was introduced to take into account the time required for adaptive immune system to respond against infection. We extensively analyzed this multidelayed model analytically and numerically. We show that delay may have both destabilizing and stabilizing effects even when the system contains a single immune response delay. It happens when there exists two sequences of critical values of this delay. If the system starts with stable state in absence of delay, then the smallest value of these critical delays causes instability and the next higher value causes stability. The system may also show stability switching for different values of the virus replication factor. These results demonstrate the possible reasons of intrapatients and interpatients variability of CD4⁺ T cells and virus counts in HIV‐infected patients.     

Study on a HIV/AIDS model with application to Yunnan province,China

This paper presents an epidemic model aiming at the prevalence of HIV/AIDS in Yunnan, China. The total population in the model is restricted within high risk population. By the epidemic characteristics of HIV/AIDS in Yunnan province, the population is divided into two groups: injecting drug users (IDUs) and people engaged in commercial sex (PECS) which includes female sex workers (FSWs), and clients of female sex workers (C). For a better understanding of HIV/AIDS transmission dynamics, we do some necessary mathematical analysis. The conditions and thresholds for the existence of four equilibria are established. We compute the reproduction number for each group independently, and show that when both the reproduction numbers are less than unity, the disease-free equilibrium is globally stable. The local stabilities for other equilibria including two boundary equilibria and one positive equilibrium are figured out. When we omit the infectivity of AIDS patients, global stability of these equilibria are obtained. For the simulation, parameters are chosen to fit as much as possible prevalence data publicly available for Yunnan. Increasing strength of the control measure on high risk population is necessary to reduce the HIV/AIDS in Yunnan.     

一类HIV病毒动力学模型的全局稳定性与周期解

建立了一类较广泛的HIV感染CD4+T细胞病毒动力学模型,给出了一个感染细胞在其整个感染期内产生的病毒的平均数(基本再生数)R0的表达式,运用Lyapunov原理和Routh-Hurwitz判据得到了该模型的未感染平衡点与感染平衡点的存在性与稳定性条件.同时也得到了模型存在轨道渐近稳定周期解和系统持续生存的条件,并通过数值模拟验证了所得到的结果.