Multiple stable periodic oscillations in a mathematical model of CTL response to HTLV-I infection

Delay induced stability switch in a mathematical model of CD8 T-cell response to SARS-CoV-2 mediated by receptor ACE2

The pathogen SARS-CoV-2 binds to the receptor angiotensin-converting enzyme 2 (ACE2) of the target cells and then replicates itself through the host, eventually releasing free virus particles. After infection, the CD8 T-cell response is triggered and appears to play a critical role in the defense against virus infections. Infected cells and their activated CD8 T-cells can cause tissue damage. Here, we established a mathematical model of within-host SARS-CoV-2 infection that incorporates the receptor ACE2, the CD8 T-cell response, and the damaged tissues. According to this model, we can get the basic reproduction number R0 and the immune reproduction number R1. We provide the theoretical proof for the stability of the disease-free equilibrium, immune-inactivated equilibrium, and immune-activated equilibrium. Finally, our numerical simulations show that the time delay in CD8 T-cell production can induce complex dynamics such as stability switching. These results provide insights into the mechanisms of SARS-CoV-2 infection and may help in the development of effective drugs against COVID-19.

Stability and Hopf bifurcation of an HIV infection model with two time delays

This work focuses on an HIV infection model with intracellular delay and immune response delay, in which the former delay refers to the time it takes for healthy cells to become infectious after infection, and the latter delay refers to the time when immune cells are activated and induced by infected cells. By investigating the properties of the associated characteristic equation, we derive sufficient criteria for the asymptotic stability of the equilibria and the existence of Hopf bifurcation to the delayed model. Based on normal form theory and center manifold theorem, the stability and the direction of the Hopf bifurcating periodic solutions are studied. The results reveal that the intracellular delay cannot affect the stability of the immunity-present equilibrium, but the immune response delay can destabilize the stable immunity-present equilibrium through the Hopf bifurcation. Numerical simulations are provided to support the theoretical results.

Dynamics and control strategy for a delayed viral infection model

In this paper, we derive a delayed epidemic model to describe the characterization of cytotoxic T lymphocyte (CTL)-mediated immune response against virus infection. The stability of equilibria and the existence of Hopf bifurcation are analysed. Theoretical results reveal that if the basic reproductive number is greater than 1, the positive equilibrium may lose its stability and the bifurcated periodic solution occurs when time delay is taken as the bifurcation parameter. Furthermore, we investigate an optimal control problem according to the delayed model based on the available therapy for hepatitis B infection. With the aim of minimizing the infected cells and viral load with consideration for the treatment costs, the optimal solution is discussed analytically. For the case when periodic solution occurs, numerical simulations are performed to suggest the optimal therapeutic strategy and compare the model-predicted consequences.

Stability analysis of general delayed HTLV-I dynamics model with mitosis and CTL immunity

This paper formulates and analyzes a general delayed mathematical model which describe the within-host dynamics of Human T-cell lymphotropic virus class I (HTLV-I) under the effect Cytotoxic T Lymphocyte (CTL) immunity. The models consist of four components: uninfected CD$ 4^{+} $T cells, latently infected cells, actively infected cells and CTLs. The mitotic division of actively infected cells are modeled. We consider general nonlinear functions for the generation, proliferation and clearance rates for all types of cells. The incidence rate of infection is also modeled by a general nonlinear function. These general functions are assumed to be satisfy some suitable conditions. To account for series of events in the infection process and activation of latently infected cells, we introduce two intracellular distributed-time delays into the models: (ⅰ) delay in the formation of latently infected cells, (ⅱ) delay in the activation of latently infected cells. We determine a bounded domain for the system's solutions. We calculate two threshold numbers, the basic reproductive number $ R_{0} $ and the CTL immunity stimulation number $ R_{1} $. We determine the conditions for the existence and global stability of the equilibrium points. We study the global stability of all equilibrium points using Lyapunov method. We prove the following: (a) if $ R_{0}\leq 1 $, then the infection-free equilibrium point is globally asymptotically stable (GAS), (b) if $ R_{1}\leq 1 < R_{0} $, then the infected equilibrium point without CTL immunity is GAS, (c) if $ R_{1} > 1 $, then the infected equilibrium point with CTL immunity is GAS. We present numerical simulations for the system by choosing special shapes of the general functions. The effects of proliferation of CTLs and time delay on the HTLV-I progression is investigated. We noted that the CTL immunity does not play the role in clearing the HTLV-I from the body, but it has an important role in controlling and suppressing the viral infection. On the other hand, we observed that, increasing the time delay intervals can have similar influences as drug therapies in removing viruses from the body. This gives some impression to develop two types of treatments, the first type aims to extend the intracellular delay periods, while the second type aims to activate and stimulate the CTL immune response.

Threshold dynamics of a viral infection model with defectively infected cells

In this paper, we investigate the global dynamics of a viral infection model with defectively infected cells. The explicit expression of the basic reproduction number of virus is obtained by using the next generation matrix approach, where each term has a clear biological interpretation. We show that the basic reproduction number serves as a threshold parameter. The virus dies out if the basic reproduction number is not greater than unity, otherwise the virus persists and the viral load eventually approaches a positive number. The result is established by Lyapunov's direct method. Our novel arguments for the stability of the infection equilibrium not only simplify the analysis (compared with some traditional ones in the literature) but also demonstrate some correlation between the two Lyapunov functions for the infection-free and infection equilibria.

Analysis of an HTLV/HIV dual infection model with diffusion

In the literature, several HTLV-I and HIV single infections models with spatial dependence have been developed and analyzed. However, modeling HTLV/HIV dual infection with diffusion has not been studied. In this work we derive and investigate a PDE model that describes the dynamics of HTLV/HIV dual infection taking into account the mobility of viruses and cells. The model includes the effect of Cytotoxic T lymphocytes (CTLs) immunity. Although HTLV-I and HIV primarily target the same host, CD4⁺T cells, via infected-to-cell (ITC) contact, however the HIV can also be transmitted through free-to-cell (FTC) contact. Moreover, HTLV-I has a vertical transmission through mitosis of active HTLV-infected cells. The well-posedness of solutions, including the existence of global solutions and the boundedness, is justified. We derive eight threshold parameters which govern the existence and stability of the eight steady states of the model. We study the global stability of all steady states based on the construction of suitable Lyapunov functions and usage of Lyapunov-LaSalle asymptotic stability theorem. Lastly, numerical simulations are carried out in order to verify the validity of our theoretical results.

Stability of an HTLV-HIV coinfection model with multiple delays and CTL-mediated immunity

In the literature, several mathematical models have been formulated and developed to describe the within-host dynamics of either human immunodeficiency virus (HIV) or human T-lymphotropic virus type I (HTLV-I) monoinfections. In this paper, we formulate and analyze a novel within-host dynamics model of HTLV-HIV coinfection taking into consideration the response of cytotoxic T lymphocytes (CTLs). The uninfected CD 4 + T cells can be infected via HIV by two mechanisms, free-to-cell and infected-to-cell. On the other hand, the HTLV-I has two modes for transmission, (i) horizontal, via direct infected-to-cell touch, and (ii) vertical, by mitotic division of active HTLV-infected cells. It is well known that the intracellular time delays play an important role in within-host virus dynamics. In this work, we consider six types of distributed-time delays. We investigate the fundamental properties of solutions. Then, we calculate the steady states of the model in terms of threshold parameters. Moreover, we study the global stability of the steady states by using the Lyapunov method. We conduct numerical simulations to illustrate and support our theoretical results. In addition, we discuss the effect of multiple time delays on stability of the steady states of the system.

Modeling and analysis of a within-host HIV/HTLV-I co-infection

Human immunodeficiency virus (HIV) and human T-lymphotropic virus type I (HTLV-I) are two retroviruses that attack the CD4 + T cells and impair their functions. Both HIV and HTLV-I can be transmitted between individuals through direct contact with certain body fluids from infected individuals. Therefore, a person can be co-infected with both viruses. HIV causes acquired immunodeficiency syndrome (AIDS), while HTLV-I is the causative agent for adult T-cell leukemia (ATL) and HTLV-I-associated myelopathy/tropical spastic paraparesis (HAM/TSP). Several mathematical models have been developed in the literature to describe the within-host dynamics of HIV and HTLV-I mono-infections. However, modeling a within-host dynamics of HIV/HTLV-I co-infection has not been involved. The present paper is concerned with the formulation and investigation of a new HIV/HTLV-I co-infection model under the effect of Cytotoxic T lymphocytes (CTLs) immune response. The model describes the interaction between susceptible CD4 + T cells, silent HIV-infected cells, active HIV-infected cells, silent HTLV-infected cells, Tax-expressing HTLV-infected cells, free HIV particles, HIV-specific CTLs and HTLV-specific CTLs. The HIV can spread by virus-to-cell transmission. On the other side, HTLV-I has two modes of transmission, (i) horizontal transmission via direct cell-to-cell contact through the virological synapse, and (ii) vertical transmission through the mitotic division of Tax-expressing HTLV-infected cells. The well-posedness of the model is established by showing that the solutions of the model are nonnegative and bounded. We define a set of threshold parameters which govern the existence and stability of all equilibria of the model. We explore the global asymptotic stability of all equilibria by utilizing Lyapunov function and Lyapunov-LaSalle asymptotic stability theorem. We have presented numerical simulations to justify the applicability and effectiveness of the theoretical results. In addition, we evaluate the effect of HTLV-I infection on the HIV dynamics and vice versa.

Modeling and stability analysis of HIV/HTLV-I co-infection

Human immunodeficiency virus (HIV) and human T-lymphotropic virus type I (HTLV-I) are two retroviruses that infect the susceptible CD[Formula: see text]T cells. It is known that HIV and HTLV-I have in common a way of transmission through direct contact with certain body fluids related to infected individuals. Therefore, it is not surprising that a mono-infected person with one of these viruses can be co-infected with the other virus. In the literature, a great number of mathematical models has been presented to describe the within-host dynamics of HIV or HTLV-I mono-infection. However, the within-host dynamics of HIV/HTLV-I co-infection has not been modeled. In this paper, we develop a new within-host HIV/HTLV-I co-infection model. The model includes the impact of Cytotoxic T lymphocytes (CTLs) immune response, which is important to control the progression of viral co-infection. The model describes the interaction between susceptible CD[Formula: see text]T cells, silent HIV-infected cells, active HIV-infected cells, silent HTLV-infected cells, Tax-expressing HTLV-infected cells, free HIV particles, HIV-specific CTLs and HTLV-specific CTLs. We first show the nonnegativity and boundedness of the model’s solutions and then we calculate all possible equilibria. We derive the threshold parameters which govern the existence and stability of all equilibria of the model. We prove the global asymptotic stability of all equilibria by utilizing Lyapunov function and LaSalle’s invariance principle. We have presented numerical simulations to illustrate the effectiveness of our main results. In addition, we discuss the effect of HTLV-I infection on the HIV-infected patients and vice versa.

Analysis of a within-host HIV/HTLV-I co-infection model with immunity

Human immunodeficiency virus (HIV) and human T-lymphotropic virus type I (HTLV-I) are two retroviruses that attack the immune cells and impair their functions. Both HIV and HTLV-I can be transmitted between individuals through direct contact with certain body fluids from infected individuals. Therefore, a person can be co-infected with both viruses. HIV causes acquired immunodeficiency syndrome, while HTLV-I is the causative agent for adult T-cell leukemia and HTLV-I-associated myelopathy/tropical spastic paraparesis. Several mathematical models have been developed in the literature to describe the within-host dynamics of HIV and HTLV-I mono-infections. However, modeling a within-host dynamics of HIV/HTLV-I co-infection has not been involved. In the present paper, we are concerned to formulate and analyze a new HIV/HTLV co-infection model under the effect of Cytotoxic T lymphocytes (CTLs) immune response. The model describes the interaction between susceptible CD4⁺T cells, silent HIV-infected cells, active HIV-infected cells, silent HTLV-infected cells, Tax-expressing HTLV-infected cells, free HIV particles, HIV-specific CTLs and HTLV-specific CTLs. The HIV can spread by two routes of transmission, virus-to-cell and cell-to-cell. On the other side, HTLV-I has two modes of transmission, (i) horizontal transmission via direct cell-to-cell contact, and (ii) vertical transmission through mitotic division of Tax-expressing HTLV-infected cells. The well-posedness of the model is established by showing that the solutions of the model are nonnegative and bounded. We define a set of threshold parameters which govern the existence and stability of all equilibria of the model. We explore the global asymptotic stability of all equilibria by utilizing Lyapunov function and LaSalle's invariance principle. We have presented numerical simulations to justify the applicability and effectiveness of the theoretical results. In addition, we evaluate the effect of HTLV-I infection on the HIV dynamics and vice versa.

Stability and bifurcation analysis for a delayed viral infection model with full logistic proliferation

In this paper, we study a delayed viral infection model with cellular infection and full logistic proliferations for both healthy and infected cells. The global asymptotic stabilities of the equilibria are studied by constructing Lyapunov functionals. Moreover, we investigated the existence of Hopf bifurcation at the infected equilibrium by regarding the possible combination of the two delays as bifurcation parameters. The results show that time delays may destabilize the infected equilibrium and lead to Hopf bifurcation. Finally, numerical simulations are carried out to illustrate the main results and explore the dynamics including Hopf bifurcation and stability switches.

Stability and bifurcation of a delayed diffusive predator-prey system with food-limited and nonlinear harvesting

Based on ecological significance, a delayed diffusive predator-prey system with food-limited and nonlinear harvesting subject to the Neumann boundary conditions is investigated in this paper. Firstly, the sufficient conditions of the stability of nonnegative constant steady state solutions of system are derived. The existence of Hopf bifurcation is obtained by analyzing the associated characteristic equation and the conditions of Turing instability are derived when the system has no delay. Furthermore, the occurrence conditions the Hopf bifurcation are discussed by regarding delay expressing the gestation time of the predator as the bifurcation parameter. Secondly, by using upper-lower solution method, the global asymptotical stability of a unique positive constant steady state solution of system is investigated. Moreover, we also give the detailed formulas to determine the direction, stability of Hopf bifurcation by applying the normal form theory and center manifold reduction. Finally, numerical simulations are carried out to demonstrate our theoretical results.

The stationary distribution and extinction of a double thresholds HTLV-I infection model with nonlinear CTL immune response disturbed by white noise

This paper investigates the stochastic HTLV-I infection model with CTL immune response, and the corresponding deterministic model has two basic reproduction numbers. We consider the nonlinear CTL immune response for the interaction between the virus and the CTL immune cells. Firstly, for the theoretical needs of system dynamical behavior, we prove that the stochastic model solution is positive and global. In addition, we obtain the existence of ergodic stationary distribution by stochastic Lyapunov functions. Meanwhile, sufficient condition for the extinction of the stochastic system is acquired. Reasonably, the dynamical behavior of deterministic model is included in our result of stochastic model when the white noise disappears.

Inhibitory killer cell immunoglobulin-like receptors strengthen CD8+ T cell-mediated control of HIV-1, HCV, and HTLV-1

Killer cell immunoglobulin-like receptors (KIRs) are expressed predominantly on natural killer cells, where they play a key role in the regulation of innate immune responses. Recent studies show that inhibitory KIRs can also affect adaptive T cell-mediated immunity. In mice and in human T cells in vitro, inhibitory KIR ligation enhanced CD8+ T cell survival. To investigate the clinical relevance of these observations, we conducted an extensive immunogenetic analysis of multiple independent cohorts of HIV-1-, hepatitis C virus (HCV)-, and human T cell leukemia virus type 1 (HTLV-1)-infected individuals in conjunction with in vitro assays of T cell survival, analysis of ex vivo KIR expression, and mathematical modeling of host-virus dynamics. Our data suggest that functional engagement of inhibitory KIRs enhances the CD8+ T cell response against HIV-1, HCV, and HTLV-1 and is a significant determinant of clinical outcome in all three viral infections.

Dynamics of an infection model with two delays

In this paper, a HTLV-I infection model with two delays is considered. It is found that the dynamics of this model are determined by two threshold parameters R0 and R1, basic reproduction numbers for viral infection and for CTL response, respectively. If R0 < 1, the infection-free equilibrium P0 is globally asymptotically stable. If R1 < 1 < R0, the asymptomatic-carrier equilibrium P1 is globally asymptotically stable. If R1 > 1, there exists a unique HAM/TSP equilibrium P2. The stability of P2 is changed when the second delay τ2 varies, that is there exist stability switches for P2.

Stability and multi-parametric Hopf bifurcation analyses of viral infection model with time delay

Ever since HIV was first diagnosed in human, a great number of scientific works have been undertaken to explore the biological mechanisms involved in the infection and progression of the disease. This paper deals with stability and bifurcation analyses of mathematical model that represents the dynamics of HIV infection of thymus. The existence and stability of the equilibria are investigated. The model is described by a system of delay differential equations with logistic growth term, cure rate and discrete type of time delay. Choosing the time delay as a bifurcation parameter, the analysis is mainly focused on the Hopf bifurcation problem to predict the existence of a limit cycle bifurcating from the infected steady state. Further, using center manifold theory and normal form method we derive explicit formulae to determine the stability and direction of the limit cycles. Moreover the mitosis rate r also plays a vital role in the model, so we fix it as second bifurcation parameter in the incidence of viral infection. Our analysis shows that, while both the bifurcation parameters can destabilize the equilibrium E* and cause limit cycles. Numerical simulations are performed to investigate the qualitative behaviors of the inherent model.

HIV evolution and progression of the infection to AIDS

In this paper, we propose and discuss a possible mechanism, which, via continuous mutations and evolution, eventually enables HIV to break from immune control. In order to investigate this mechanism, we employ a simple mathematical model, which describes the relationship between evolving HIV and the specific CTL response and explicitly takes into consideration the role of CD4(+)T cells (helper T cells) in the activation of the CTL response. Based on the assumption that HIV evolves towards higher replication rates, we quantitatively analyze the dynamical properties of this model. The model exhibits the existence of two thresholds, defined as the immune activation threshold and the immunodeficiency threshold, which are critical for the activation and persistence of the specific cell-mediated immune response: the specific CTL response can be established and is able to effectively control an infection when the virus replication rate is between these two thresholds. If the replication rate is below the immune activation threshold, then the specific immune response cannot be reliably established due to the shortage of antigen-presenting cells. Besides, the specific immune response cannot be established when the virus replication rate is above the immunodeficiency threshold due to low levels of CD4(+)T cells. The latter case implies the collapse of the immune system and beginning of AIDS. The interval between these two thresholds roughly corresponds to the asymptomatic stage of HIV infection. The model shows that the duration of the asymptomatic stage and progression of the disease are very sensitive to variations in the model parameters. In particularly, the rate of production of the naive lymphocytes appears to be crucial.

Joint effects of mitosis and intracellular delay on viral dynamics: two-parameter bifurcation analysis

To understand joint effects of logistic growth in target cells and intracellular delay on viral dynamics in vivo, we carry out two-parameter bifurcation analysis of an in-host model that describes infections of many viruses including HIV-I, HBV and HTLV-I. The bifurcation parameters are the mitosis rate r of the target cells and an intracellular delay τ in the incidence of viral infection. We describe the stability region of the chronic-infection equilibrium E* in the two-dimensional (r, τ) parameter space, as well as the global Hopf bifurcation curves as each of τ and r varies. Our analysis shows that, while both τ and r can destabilize E* and cause Hopf bifurcations, they do behave differently. The intracellular delay τ can cause Hopf bifurcations only when r is positive and sufficiently large, while r can cause Hopf bifurcations even when τ = 0. Intracellular delay τ can cause stability switches in E* while r does not.