Advertisement

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Modeling the Slow CD4+ T Cell Decline in HIV-Infected Individuals

  • Sunpeng Wang,

    Affiliation Department of Biology, New York University, New York, New York, United States of America

    ⨯
  • Patricia Hottz,

    Affiliation Departamento de Medicina Preventiva, Faculdade de Medicina, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil

    ⨯
  • Mauro Schechter,

    Affiliations Departamento de Medicina Preventiva, Faculdade de Medicina, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil, Projeto Praça Onze, Hospital Escola São Francisco de Assis, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil

    ⨯
  • Libin Rong

    rong2@oakland.edu

    Affiliation Department of Mathematics and Statistics, and Center for Biomedical Research, Oakland University, Rochester, Michigan, United States of America

    ⨯

Abstract

The progressive loss of CD4+ T cell population is the hallmark of HIV-1 infection but the mechanism underlying the slow T cell decline remains unclear. Some recent studies suggested that pyroptosis, a form of programmed cell death triggered during abortive HIV infection, is associated with the release of inflammatory cytokines, which can attract more CD4+ T cells to be infected. In this paper, we developed mathematical models to study whether this mechanism can explain the time scale of CD4+ T cell decline during HIV infection. Simulations of the models showed that cytokine induced T cell movement can explain the very slow decline of CD4+ T cells within untreated patients. The long-term CD4+ T cell dynamics predicted by the models were shown to be consistent with available data from patients in Rio de Janeiro, Brazil. Highly active antiretroviral therapy has the potential to restore the CD4+ T cell population but CD4+ response depends on the effectiveness of the therapy, when the therapy is initiated, and whether there are drug sanctuary sites. The model also showed that chronic inflammation induced by pyroptosis may facilitate persistence of the HIV latent reservoir by promoting homeostatic proliferation of memory CD4+ cells. These results improve our understanding of the long-term T cell dynamics in HIV-1 infection, and support that new treatment strategies, such as the use of caspase-1 inhibitors that inhibit pyroptosis, may maintain the CD4+ T cell population and reduce the latent reservoir size.

Author Summary

The CD4+ T cell population within HIV-infected individuals declines slowly as disease progresses. When CD4+ cells drop to below 200 cells/ul, the infection is usually considered to enter the late stage, i.e., acquired immune deficiency syndrome (AIDS). CD4+ T cell depletion can take many years but the biological events underlying such slow decline are not well understood. Some studies showed that the majority of infected T cells in lymph nodes die by pyroptosis, a form of programmed cell death, which can release inflammatory signals attracting more CD4+ T cells to be infected. We developed mathematical models to describe this process and explored whether they can generate the long-term CD4+ T cell decline. We showed that pyroptosis induced cell movement can explain the slow time scale of CD4+ T cell depletion and that pyroptosis may also contribute to the persistence of latently infected cells, which represent a major obstacle to HIV eradication. The modeling prediction agrees with patient data in Rio de Janeiro, Brazil. These results suggest that a combination of current treatment regimens and caspase-1 inhibitor that can inhibit pyroptosis might provide a new way to maintain the CD4+ T cell population and eradicate the HIV latent reservoir.

Citation: Wang S, Hottz P, Schechter M, Rong L (2015) Modeling the Slow CD4+ T Cell Decline in HIV-Infected Individuals. PLoS Comput Biol 11(12): e1004665. https://doi.org/10.1371/journal.pcbi.1004665

Editor: Becca Asquith, Imperial College London, UNITED KINGDOM

Received: March 24, 2015; Accepted: November 17, 2015; Published: December 28, 2015

Copyright: © 2015 Wang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: LR is supported by NSF grants DMS-1122290 and DMS-1349939. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

HIV-1 progression to the AIDS stage within untreated patients usually takes many years. As HIV-1 infection progresses, the CD4+ T cell population declines slowly and the infected individual becomes progressively more susceptible to certain opportunistic infections and neoplasms. These are particularly common when CD4+ T cells reach a level below 200 cells/ul, which defines AIDS [1–7]. How HIV-1 infection induces progressive CD4+ T cell depletion is unclear [8]. One explanation is that the turnover rate of CD4+ T cells is significantly increased in HIV or simian immunodeficiency virus (SIV) infected subjects [9,10]. Therefore, massive activation of CD4+ T cells, which leads to more viral infection and cell death, might outrun the regeneration of T cells and cause progressive depletion. Another explanation is the failure of CD4+ memory T cell homeostasis during progressive HIV infection. This is possibly due to the destruction of the microenvironment of organs and tissues supporting T cell regeneration [3,11–14]. It remains unclear whether the impaired conformation of T cell regenerative tissues leads to the regeneration failure or it is merely a pathogenic reformation caused by HIV to promote viral replication.

Mathematical models may shed light on how the complex interplay between the immune response and viral infection leads to overt immunodeficiency. Matrajt et al. used a model to analyze the simian-human immunodeficiency virus (SHIV) infection data in macaques [15]. They found that uninfected or bystander cell death accounts for the majority of CD4+ T cell death [15]. Mohri et al. studied the turnover of CD4+ T cells and found that T cell depletion is primarily induced by increased cellular destruction rather than decreased cellular production [16]. Kovacs et al. also showed that HIV does not impair CD4+ T cell production but increases T cell proliferation [17]. Using a model including the activation of resting CD4+ T cells, Ribeiro et al. found that HIV infection increases both the activation rate of resting CD4+ T cells and the rates of death and proliferation of activated CD4+ T cells [18]. Chan et al. showed that the rapid proliferation of CD4+ T cells provides more targets for infection and that preservation of CD4+ T cells in natural host monkeys is due to the limited CD4+ T cell proliferation [19]. Thus, CD4+ T cell depletion may be caused by the massive immune activation during chronic infection. However, a model by Yates et al. suggested that if immune activation drives T cell decline, then the predicted decline would be very fast, which is not consistent with the time scale of T cell depletion during chronic infection [20]. The above observations and analyses may explain T cell depletion but the long-term dynamics of CD4+ T cells have been neither simulated by models nor compared with patient data. In a recent study, Hernandez-Vargas and Middleton [21] developed a model including the infection of macrophages to explain the three stages of HIV infection. Fast infection of CD4+ T cells can explain the CD4+ T cell and viral load dynamics in the early stages, while slow infection of macrophages may explain the dynamics in the advanced stages of infection. Whether macrophages form a long-term reservoir causing T cell depletion and viral explosion in the later stages of infection needs further experimental investigation.

Different from apoptosis, a programmed process that results in non-inflammatory cell death, pyroptosis is a form of programmed cell death associated with antimicrobial responses during inflammation [22]. During HIV infection, Doitsh et al. [23,24] found that when virus enters a CD4+ T cell that is non-permissive to viral infection, the caspase-1 pathway is triggered to induce pyroptosis, which can secrete inflammatory cytokines such as IL-1β. These cytokines establish a chronic inflammation state and attract more CD4+ T cells to the inflamed sites, resulting in more infection and cell death. Thus, pyroptosis generates a vicious cycle in which dying CD4+ T cells secrete inflammatory signals that attract more CD4+ T cells to be infected and die [23]. These findings suggest that HIV-1 may use the intrinsic feature of the immune system to seek targets of infection, establish productive viral replication, and meanwhile destroy the CD4+ T cell population.

Here we developed mathematical models incorporating the effect of pyroptosis to study whether it can explain the very slow T cell depletion during HIV-1 infection. Using the models we explored if highly active antiretroviral therapy (HAART) can preserve the CD4+ T cell population. We studied the effect of CD4+ T cell proliferation and CD8+ T cell response on CD4+ decline. We also compared our modeling prediction with clinical data obtained from patients in Rio de Janeiro, Brazil [25–28]. At last, we probed the possible contribution of chronic inflammation associated with pyroptosis to the HIV latent reservoir persistence.

Methods

Patient data

The patient data were obtained from seroconverters in 3 cohorts [25–28]. One cohort consists of high-risk, HIV-seronegative homosexual and bisexual men who did not report injection drug use, were enrolled between July 1995 and June 1998 and seroconverted during follow-up [26]. The other cohorts consist of seroconverters from high-risk HIV-seronegative homosexual and bisexual men patients who were enrolled from December 1998 to May 2001 in a study designed to evaluate the behavior impact of post-exposure prophylaxis [27], and participants from the control arm of SPARTAC, a randomized trial designed to evaluate the impact of short term antiretroviral therapy on the course of primary HIV infection [28]. The median of the CD4+ T cell data was derived from these cohort studies. The median of the Current Study Multicenter AIDS Cohort Study (MACS) was obtained from the study [29]. These patient data and medians were compared with modeling prediction.

One-compartment model

Inflammatory cytokines released by abortively HIV-infected cells can attract more CD4+ T cells to be infected. In the following one-compartment model, to minimize the number of variables and parameters we described the effect of pyroptosis by use of an enhanced viral infection rate because of increased availability of CD4+ T cells attracted by cytokines to the inflamed sites.

The variable T represents the population of uninfected CD4+ T cells. They are generated at the rate λ. Proliferation of target cells will be considered later. The infection rate is modeled by a mass action term kVT, which is enhanced by the inflammatory cytokine (C) with a factor γi. Uninfected T cells die at a per capita rate d1. T* is the population of productively infected T cells and their death rate is d2. A fraction (f) of new infection is assumed to be abortively infected. The death rate of abortively infected T cells (M*) is d3. Virus (V) is generated by productively infected T cells with a viral production rate pv and is cleared at a rate d4. Inflammatory cytokines are released with a burst size (Nc) when an abortively infected cell dies. Thus, Ncd3 represents the generation rate of cytokines per abortively infected cell. The decay rate of cytokines is assumed to be d5. The schematic diagram of this model is shown in Fig 1. Parameters and values are listed in Table 1.

thumbnail
Download:
Fig 1. Schematic diagram of the one-compartment model.

Inflammatory cytokines released during cell death by pyroptosis attract more CD4+ T cells (T) to be infected. The term γiC·kVT represents cytokine enhanced viral infection due to increased CD4+ T cell availability.

https://doi.org/10.1371/journal.pcbi.1004665.g001

thumbnail
Download:
Table 1. A summary of all parameters and values used in the models.

https://doi.org/10.1371/journal.pcbi.1004665.t001

In the above one-compartment model, we described the consequence of pyroptosis but did not explicitly model the cytokine-induced attraction of CD4+ T cells from elsewhere to the place where abortive infection occurs. Below we developed another model with two compartments to include cytokine-induced T cell movement explicitly. The model is more complicated and contains more parameters.

Two-compartment model

In the model there are two compartments: one represents the blood (T1) and the other represents human lymphoid tissues (T2) such as lymph nodes in which abortive infection takes place on a large scale [23]. CD4+ T cells in compartment I (or II) can transport to compartment II (or I) at a rate σ1 (or σ2). In blood, cytokines released during abortive infection cannot accumulate as in lymphoid tissues. They cannot attract other immune cells to fight the infection and contribute to inflammation. Thus, pyroptosis is assumed to take place only in lymphoid tissues (compartment II), as observed in ref. [23]. The transportation rate σ1 from the blood to tissues is assumed to be enhanced by a factor (1+γrC) due to inflammatory cytokines (C) released during pyroptosis in compartment II. Viruses (V1 and V2) can also transport between two compartments with the rates D2(V1-V2) and D1(V2-V1), which depend on the difference of viral load in the two compartments. Because the dynamics of the virus are much faster than those of infected cells, it is reasonable to assume that they are proportional to each other. Thus, we only included the transportation of virus between compartments. In the Supporting Information (S1 Text and S7 Fig), we added the transportation of infected cells to the model and found that the model prediction is similar to the case without infected cell transportation. All the other variables and parameters (summarized in Table 1) can be defined similarly as those in the one-compartment model (Fig 1). The schematic diagram of the two-compartment model is shown in Fig 2.

thumbnail
Download:
Fig 2. Schematic diagram of the two-compartment model.

Parameters σ1 and σ2 represent the transportation rates of uninfected T cells between blood and lymph nodes. Cytokines (C) attract uninfected T cells (T1) from blood into lymph nodes at an enhanced rate γrC. Parameters D1 and D2 represent the transportation rates of virus between two compartments.

https://doi.org/10.1371/journal.pcbi.1004665.g002

Model parameters

For model simulation, we fixed most of parameters based on existing experimental data and our previous modeling studies [30–33]. Because the CD4+ T cell level within an uninfected individual ranges normally from 500 cells/μl to 1500 cells/μl, we changed the unit to cells/ml and assumed CD4+ T cells to be 106 cells/ml before infection [34]. The death rate (d1) of uninfected CD4+ T cells is assumed to be 0.01 day-1 [35]. Thus, from the steady state of target cells before infection, we obtained that the generation rate (λ) of target cells is 106(0.01) = 104 cells ml-1 day-1. The viral infection rate k is assumed to be 2.4×10−8 ml virion-1 day-1 [30]. The death rate of infected T cells is d2 = 1 day-1 [36]. We chose the parameter γi to be 2×10−4 ml molecule-1. The viral production rate of productively infected T cells in the one-compartment model is chosen to be 2.5×104 virions cell-1 day-1 [37]. As described by Doitsh et al. [23,24], abortive infection accounts for 95% of the total infection. Thus, we chose f to be 0.95. Because abortive infection mainly takes place in non-permissive quiescent T cells, we chose their death rate (d3) to be 0.001 day-1 [31,32]. The burst size of cytokines is fixed to Nc = 15 molecules. The half-life of IL-1β is about 2.5 hours [38]. Thus, we chose the decay rate of cytokines to be d5 = 6.6 day-1. We also performed sensitivity tests of the modeling prediction on a number of parameters.

Data fitting

We fit both the one-compartment and two-compartment models to subjects with more than 10 data points [25–29]. The root mean square (RMS) between model prediction and patient data is minimized for each patient. RMS is calculated using the following formula where T(ti)+T*(ti) represents the CD4+ T cell population level in blood at time ti predicted by the model, is the corresponding patient data at ti. We used T1(ti)+T1*(ti) in the fitting for the two-compartment model. Parameter estimates are based on the best fit that achieves the minimum RMS. Data fitting is performed using the R programming language.

Model comparison by AIC

In order to statistically compare the best fits of using the two models, we calculated the Akaike information criterion (AIC). The model with a lower AIC value fits the data better from a statistical viewpoint. The AIC is calculated using the following formula where n is the number of observations (i.e. number of data points) and m is the number of fitted parameters. RSS is the residual sum of squares. T(ti), T*(ti) and are the same as those defined in the calculation of RMS.

Confidence interval

We obtained the 95% confidence intervals for fitted parameters using a bootstrap method [39], where the residuals to the best fit were re-sampled 200 times.

Results

Slow depletion of CD4+ T cells

Using the parameter values listed in Methods and initial values V(0) = 1×10−3 RNA copies/ml, T(0) = 103 cells/μl, T*(0) = 0, M*(0) = 0, and C(0) = 0 in the one-compartment model, we showed that the population of CD4+ count declines from 103 cells/μl to about 200 cells/μl around the 6th year after infection (Fig 3A). This is consistent with the slow time scale of T cell decline during HIV infection. The entire T cell depletion course consists of two major phases. The first massive depletion phase is rapid, followed by a slower chronic depletion phase (Fig 3A). The first-phase T cell decline is due to the substantial viral infection during the early stage. If there is no infection (k = 0), then the T cell level would stabilize at the initial level (Fig 3A). The slow second-phase T cell decline is due to pyroptosis enhanced viral infection. Without the effect of inflammatory cytokines released during pyroptosis (i.e. γi = 0 or no inflammation in Fig 3A), a balance between T cell generation and viral infection is reached and the T cell population is maintained at a steady state level. This agrees with the prediction of most viral dynamics models without treatment. Because of pyroptosis, cytokine-enhanced viral infection breaks the balance between cellular production and viral infection, which makes the T cell level decline at a very low rate and approach the immune-deficient level after several years (Fig 3A).

thumbnail
Download:
Fig 3. CD4+ T cell decline and viral load dynamics predicted by the one-compartment model.

(A) CD4+ T cell decline consists of two major phases. A rapid and massive decline is caused by enormous viral infection during the early stage. It is followed by a progressive depletion phase, which is driven by pyroptosis enhanced viral infection. (B) Viral load dynamics generated by the one-compartment model with and without chronic inflammation.

https://doi.org/10.1371/journal.pcbi.1004665.g003

The viral load change was plotted in Fig 3B. Without the effect of inflammatory cytokines, the viral load reaches a steady state level. When there is cytokine enhanced viral infection, viral load increases very slowly during the phase of chronic infection (Fig 3B).

Using a constant λ is a simple way to approximate the generation of target cells. We included the proliferation of target cells in the model (S1 Text). Simulation with different proliferation rates is shown in S1 Fig. As the proliferation rate increases, the decline of CD4+ T cells becomes faster. This is because more target cells lead to more abortive infection, which releases more cytokines attracting more CD4+ T cells to be infected and die. This prediction is consistent with the observation that the level of T cell proliferation in non-pathogenic infection (e.g. SIV infection in natural host monkeys such as sooty mangabeys or mandrills that do not develop AIDS-like diseases) was much lower than in pathogenic infection, e.g., SIV in rhesus macaques [40,41]. This provides an additional support to the view that an attenuated rather than effective adaptive immune response preserves immune function in natural host monkeys [42].

We performed sensitivity analysis of the CD4+ T cell decline for a number of parameters. Fig 4 shows that the sensitivity tests on parameters k, λ, pv and γi. S2–S5 Figs show the tests on parameters Nc, d3, d5, and f, respectively. We found that the model is robust in generating the slow decline of CD4+ T cells, although the model prediction is more sensitive to three parameters k, pv and f (see Figs 4A, 4C and S5).

thumbnail
Download:
Fig 4. CD4+ T cell dynamics predicted by the one-compartment model with different parameter values.

(A) Sensitivity test on the infection rate k. The value of k increases from 1.44×10−8 (orange line) to 3.36 ×10−8 ml virion-1 day-1 (black line) with four equal increments. (B) Sensitivity test on the target cell generate rate λ. The value of λ increases from 6000 (orange line) to 14000 cell ml-1day-1 (black line) with four equal increments. (C) Sensitivity test on the viral production rate pv. The value of pv increases from 15000 (orange line) to 35000 virions cell-1 day-1 (black line) with four equal increments. (D) Sensitivity test on γi, the effect of cytokines on infection. The value of γi increases from 1.2×10−4 (orange line) to 2.8 ×10−4 ml molecule-1 (black line) with four equal increments.

https://doi.org/10.1371/journal.pcbi.1004665.g004

In the above simulation, we assumed that the viral infection enhancement parameter γi is a constant. When the concentration of inflammatory cytokines is low, they may not be able to trigger the attraction of CD4+ T cells from elsewhere. Thus, we simulated a scenario in which enhanced viral infection is triggered only when the level of cytokines is above a threshold value. We chose γi to be the following step function. It is zero if the level of cytokines is below a certain level.

The threshold value was chosen to be 2000 or 4000 molecules/ml in Fig 5A. CD4+ T cells do not decline until the level of cytokines reaches the corresponding threshold (Fig 5B).

thumbnail
Download:
Fig 5. Simulation with a concentration-dependent cytokine enhanced function γi.

(A) The dynamics of cytokines. The black and brown dashed lines represent two threshold values we chose in simulation. (B) CD4+ T cell dynamics generated by the one-compartment model with a step function for the parameter γi. The brown line corresponds to the threshold of 2000 molecules/ml and the black corresponds to the threshold of 4000 molecules/ml. (C) CD4+ T cell dynamics generated by the one-compartment model when the parameter γi is an exponential function of the concentration of cytokines.

https://doi.org/10.1371/journal.pcbi.1004665.g005

A more realistic scenario is that γi increases gradually when the concentration of cytokines is above the threshold. We chose γi(C) to be the following exponential function.

The hill coefficient ρ determines how fast γi(C) increases from 0 to its maximum value γi. Both ρ and γi were fixed to 2×10−4 ml molecule-1. With a non-constant parameter γi(C), we found that CD4+ T cells also undergo a slow decline to below 200 cells/μl (Fig 5B and 5C). Using an exponential function for γi(C), the decline of CD4+ T cells is smoother than the case using a step function.

Influence of HAART

Using the one-compartment model, we studied if HAART can rescue the CD4+ T cell population. During HAART we assumed that the viral infection rate k is reduced by a factor (1-ε), where ε is the overall drug efficacy of the treatment [32]. The simulation shows that if the treatment effectiveness is very high, then CD4+ count can rebound to its pre-infection level (Fig 6A) no matter when HAART is initiated. For lower treatment effectiveness (e.g. ε = 0.6 in Fig 6B), the patient needs a relatively long time to restore the CD4+ T cell population. The later HAART starts, the longer it takes for CD4+ T cell restoration (Fig 6B). When the treatment effectiveness is further lower, CD4+ T cell depletion could not be prevented. These results suggest that HAART has the potential to rescue CD4+ T cell population, but CD4+ response depends on the effectiveness of the therapy and when the therapy is initiated.

thumbnail
Download:
Fig 6. Influence of HAART on CD4+ T cell dynamics.

(A) If treatment is very effective, then the CD4+ T cell population is predicted to be restored. (B) If drug efficacy is lower, then the later HAART is initiated the longer it takes for the T cell population to increase.

https://doi.org/10.1371/journal.pcbi.1004665.g006

Model with CD8+ T cell response

We included CD8+ T cells in the one-compartment model to study the interaction between CD4+ T cell decline and CD8+ T cell response. CD8+ T cells (E) are assumed to kill infected T cells at a rate αET*. The activation rate of CD8+ T cells depends on the level of infected cells with a half-maximal saturation constant θ. pE is the maximum activation rate. CD4+ T cells play an important role in activating the adaptive immune response. We used another saturation function T/(T+η) to account for this influence. The T* and E equations are given below.

The simulation of the model with CD8+ T cell response is shown in Fig 7. Parameter values are listed in Table 1. For comparison, we plotted the predicted T cell dynamics with and without the influence of CD4+ T cells. In column A of Fig 7, we performed the simulation without T/(T+η). CD4+ T cells decline slowly and CD8+ T cells reach a steady state level. Column B shows the simulation with the term T/(T+η). CD8+ T cell response becomes weaker than in column A because of the slow depletion of CD4+ T cells. CD4+ T cells decline faster because of the incapability of CD8+ T cells to control viral infection.

thumbnail
Download:
Fig 7. Simulation of the one-compartment model with CD8+ T cell response.

Column A: predicted T cell dynamics assuming that CD8 activation is not regulated by CD4+ T cells (i.e. without the term T/(T+η). Column B: predicted T cell dynamics assuming that CD8 activation is regulated by CD4+ T cells. Parameter values are listed in Table 1.

https://doi.org/10.1371/journal.pcbi.1004665.g007

Two-compartment model

Inflammatory signals released during pyroptosis induce the movement of CD4+ T cells from circulation in blood to inflamed lymph nodes [43–46]. We developed a more comprehensive model by including two cell compartments (Fig 2). One is the blood compartment and the other is the compartment of lymphoid tissues where pyroptosis takes place. Simulation of the two-compartment model shows that CD4+ count in blood declines from 103 cells/ul to 200 cells/ul over a long time period (Fig 8A). The viral load change in blood is also similar to that shown in Fig 3B except that T cell and viral load dynamics generated by the two-compartment model have less oscillation than by the one-compartment model.

thumbnail
Download:
Fig 8. Simulation of the two-compartment model.

(A) CD4+ T cell dynamics predicted by the two-compartment model. (B) Simulation of the model with ε1 = ε2 = 0.9. The red and black lines are the simulation when HAART is initiated 1 year and 3 years, respectively, after infection. (C) Long-term CD4+ T cell dynamics predicted with ε1 = 0.9 and ε2 = 0.4 when HAART is initiated at the beginning of infection. Because approximately 98% of viral replication takes place within lymph nodes [2,46], in the simulation we fixed λ1 to be 104 cell ml-1 day-1 and λ2 to be 50 times of λ1. Using the equilibrium before infection, we chose the rate σ1 to be 50 times of σ2 (σ1 = 0.01 day-1 and σ2 = 0.0002 day-1). The other parameters were chosen to be γr = 5×10−6 ml molecule-1, D1 = 0.1 day-1, D2 = 0.2 day-1, pv1 = 1000 day-1, and pv2 = 2000 day-1 [32].

https://doi.org/10.1371/journal.pcbi.1004665.g008

Using the two-compartment model we also tested if HAART can rescue CD4+ T cell population. We assumed that the drug efficacies of HAART within blood and lymph node are different (i.e., the viral infection rate k in compartment I is reduced by 1-ε1 and k in compartment II is reduced by 1-ε2). We found that if the drug efficacies in both compartments are high, then CD4+ T cell depletion can be prevented (Fig 8B). The time for CD4+ restoration also depends on when HAART is initiated. However, if the drug efficacy in compartment II is relatively low (e.g. ε2 = 0.4) compared with the high efficacy in compartment I (e.g. ε1 = 0.9), then CD4+ T cells decline even when HAART is initiated at the beginning of viral infection (Fig 8C). In the simulation, CD4+ T cells stabilize at 230 cells/ul after more than 30 years (Fig 8C). This result suggests that even if some lymphoid tissues might be difficult for drug's penetration (i.e. drug sanctuary sites), CD4+ T cells can be maintained at a higher level in treated patients than in untreated patients. This may explain the increased life expectancies in HIV patients treated with combination therapy [47–51]. However, because of the CD4+ cell decline (Fig 8C), life expectancy should be lower in patients with lower baseline CD4+ cell counts than in those with higher baseline counts. This is consistent with the reported life expectancy of individuals on combination therapy in a collaborative analysis of 14 cohort studies [47].

Comparison with long-term CD4+ T cell data

We compared modeling predictions with the CD4+ T cell data shown in [25–29]. Using the one-compartment model, we fit parameters k, γi, λ, pv and fix the other parameters for each patient. We also fit the model to the median data calculated from all the patients in the two cohort study [25] and the median data of the Current Study Multicenter AIDS Cohort Study (MACS) [29]. Using the two-compartment model, we fit parameters k, γr, λ1, pv1 to the same patient and median data. Figs 9 and 10 show that both models provide a good fit to the long-term CD4+ T cell data in untreated HIV-1 patients. The fit to the median data is better than the fit to individual patients based on the calculated error between modeling prediction and data. These data fits suggest that pyroptosis induced CD4+ T cell movement during abortive infection can explain the progressive CD4+ T cell depletion observed in untreated HIV-1 patients.

thumbnail
Download:
Fig 9. Fitting of the one-compartment model to patient data.

Patient median was derived from two cohorts of studies in ref. [25] and MACS median was derived from the Multicenter AIDS Cohort Study in ref. [29]. Parameters values based on the best fits, 95% confidence intervals, and the AIC values of the fitting are listed in Table 2.

https://doi.org/10.1371/journal.pcbi.1004665.g009

thumbnail
Download:
Fig 10. Fitting of the two-compartment model to the same patient data shown in Fig 9.

Parameters values based on the best fits, 95% confidence intervals, and the AIC values of the fitting are listed in Table 3.

https://doi.org/10.1371/journal.pcbi.1004665.g010

Parameter estimates and their 95% confidence intervals based on the fits to the one-compartment and two-compartment models are listed in Tables 2 and 3, respectively. The estimate of the viral production rate pv in the one-compartment model is higher than the viral production rate pv1 in blood of the two-compartment model (pv2 = 2000 virions per cell per day is fixed during fitting). This is because in the one-compartment model 95% of infection is assumed to be abortive and only 5% of infection produces virus. Thus, a higher value of viral production rate is needed to generate viral load with reasonable magnitude. In the second-compartment model, although only 5% of infection produces viruses in lymphoid tissues, the target cell level is much higher in lymphoid tissues than in blood (i.e. λ2 >> λ1). Thus, the viral production rates in the two compartments are on the same order of magnitude.

thumbnail
Download:
Table 2. Parameter estimates, 95% confidence intervals, and AIC values when fitting the one-compartment model to patient data.

https://doi.org/10.1371/journal.pcbi.1004665.t002

thumbnail
Download:
Table 3. Parameter estimates, 95% confidence intervals, and AIC values when fitting the two-compartment model to patient data.

https://doi.org/10.1371/journal.pcbi.1004665.t003

The Akaike information criterion (AIC) value is calculated to compare data fitting using the two models (Tables 2 and 3). We found that for patients 11, 38, 44, and median of patient data, the AIC value of using the second model is less than that of using the first model. This suggests that the two-compartment model provides a better fit to the data for these patients from a statistical viewpoint.

Latent reservoir persistence

IL-7 plays an important role in latently infected CD4+ T cell proliferation [52]. It has been observed to be over expressed in inflamed tissues [53,54]. Inflammatory cytokines released during cell death by pyroptosis may promote the establishment and persistence of the latent reservoir in HIV patients. Here we included the population of latently infected CD4+ T cell (L) into the one-compartment model. Latently infected CD4+ cells are produced with a fraction μ during HIV-1 infection. They can also be maintained by proliferation which is assumed to rely on the cytokine level (see the term 1+φC in the following equation where φ is fixed to 10−2 ml molecule-1). We chose the base proliferation rate pL to be 0.001 day-1 [32], which represents a limited proliferation capacity in the absence of inflammatory cytokines. The carrying capacity of latently infected cells (Lmax) is fixed at 100 cells/ml [32]. The other parameter values are listed in Table 1. The equations of L and T* are given below and the other equations are the same as those in the one-compartment model.

If there is no chronic inflammation (i.e. φ = 0 in the L equation), then latently infected cells undergo a slow decline (Fig 11A). However, if the proliferation is enhanced by cytokines released during cell death by pyroptosis, then the latent reservoir can be maintained at a higher level (Fig 11B). This result suggests that inflammatory cytokines generated during abortive infection might contribute to the establishment of the latent reservoir and the maintenance of its size. We also performed sensitivity test of latently infected cells on the parameter φ, the effectiveness of cytokines promoting latently infected cell proliferation. The modeling prediction is robust to this parameter (S6 Fig).

thumbnail
Download:
Fig 11. Simulation of the one-compartment model with latently infected cells.

(A) Population of latently infected T cells without cytokine enhanced homeostatic proliferation. (B) Population of latently infected T cells with cytokine enhanced homeostatic proliferation.

https://doi.org/10.1371/journal.pcbi.1004665.g011

Latently infected cells can be activated by relevant antigens and become productively infected cells. In S1 Text, we included the activation of latently infected cells in the one-compartment model. Simulation with different values of the activation rate is shown in S8 Fig. As the activation rate increases, the size of the latent reservoir decreases.

Discussion

The mechanisms underlying the slow time scale of CD4+ T cell decline in untreated HIV-1 patients remain unclear. HIV-mediated cell death can contribute to the loss of CD4+ T cells, but quantitative image analysis suggested that infection-induced cell death could be compensated by upregulated T cell division [55,56]. Some studies suggested that the destruction of bystander non-infected cells may account for the CD4+ T cell decline during disease progression [57–60]. Immune activation might be the reason of bystander cellular demise [60]. It drives uninfected CD4+ T cells into several rounds of division and cells are susceptible to activation-induced death [61,62]. However, a mathematical model showed that the decline of CD4+ T cells would be very rapid if immune activation drives T cell depletion [20]. Another possible reason of T cell decline might be the regeneration failure of CD4+ T cells during disease progression [3,11–14]. A recent study found that about 95% of CD4+ T cells within lymph nodes die from pyroptosis and release inflammatory signals that attract more CD4+ T cells from elsewhere to be infected [23]. HIV-1 may use this vicious infection cycle to promote disease progression and chronic T cell depletion.

In this paper, we developed mathematical models to explore whether cell death induced by pyroptosis can explain the slow time scale of CD4+ T cell decline in untreated HIV patients. In the first model, we assumed that increased availability of target cells due to attraction by inflammatory cytokines facilitates viral infection, which drains the CD4+ T cell population slowly during chronic infection. In the second model, we explicitly included the movement of CD4+ T cells from blood to lymphoid tissues where pyroptosis occurs. Both models generate a very slow decline of CD4+ T cells in plasma (Figs 3 and 8), and agree with the long-term CD4+ T cell data from untreated HIV patients in several cohorts in Brazil (Figs 9 and 10).

We found that the entire CD4+ T cell decline consists of two major phases (Fig 3). The first-phase decline is very rapid. This decline is due to the enormous virus infection and virus-induced cell death during primary infection. Following the first phase, CD4+ T cells partially recover because of cell regeneration and viral control by immune responses. However, a balance cannot be established between cell generation and viral infection. Chronic inflammatory cytokines released during pyroptosis can attract CD4+ T cells from other places to inflamed lymphoid tissues. These cells are infected and die, resulting in a slow decline of CD4+ T cells in plasma. These results suggest that HIV-mediated cell death causes the dramatic decline of CD4+ T cells during primary infection and that persistent chronic inflammation acts like an erosive force which gradually drains the CD4+ T cell population in plasma during chronic infection. HAART was shown to have the potential to restore the CD4+ T cell population (Figs 6 and 8), which agrees with the robust and sustained CD4 recovery among patients remaining on therapy [63] and a normal life expectancy in patients with a good CD4 response and undetectable viral load [50]. However, CD4 response depends on the effectiveness of the therapy, when the therapy is initiated, and whether there exist drug sanctuary sites (Figs 6 and 8). This may explain the considerable variability in the increase of life expectancy in patients treated with combination therapy between 1996 and 2005 [47].

Our model has limitations. First, it does not account for the spatial effect of CD4+ T cells. Although we used a two-compartment model to describe the transportation of cells and virus between blood and lymphoid tissues, release of cytokines during cell death by pyroptosis and attraction of CD4+ T cells are mainly constrained to occur locally. Ordinary differential equation models could not capture these features. It would be valuable to develop spatial models that can describe the vicious cycle within lymphoid tissues. Spatial models require precise description and parameterization of diffusion of cytokines and attraction of CD4+ T cells, and are also computationally demanding in studying T cell dynamics within blood and different lymphoid tissues. The second limitation of our model is that we did not consider a detailed inflammatory signal transduction cascade between T cells and relevant tissues. Recruitment of T cells to the inflamed tissue goes through several steps of immunological reaction. Upon secretion of IL-1β, expression of adhesion molecules such as E/P-selectin and ICAM-1 on the vascular endothelium is upregulated [64]. Binding to these molecules facilitates T cell's attachment to vascular endothelium. After attachment T cells undergo conformational changes and penetrate into the inflamed tissue [65,66]. In our models, we used a very simple factor multiplied by the concentration of cytokines to describe the effect of inflammatory cytokines. A more comprehensive model requires a detailed description of intracellular processes underlying the inflammatory signal cascade and related data for model verification. The third limitation is that our model cannot generate viral load explosion in the later stages of HIV infection. Assuming that all parameters are constant and that only one cell population produces virus, our model cannot describe viral explosion. However, as CD4+ T cells drop to very low levels, the immune system cannot kill infected cells or neutralize virus effectively. This leads to a reduction in the death rate of infected cells or viral clearance rate, and may explain the viral explosion. Infection of other cell populations such as macrophages (as suggested by Hernandez-Vargas and Middleton in ref. [21]) or other viral reservoirs may also explain the dramatic viral load increase during the AIDS stage.

Our simulation shows that the latent reservoir may be maintained by chronic inflammation. How inflammation promotes the latent reservoir persistence is not fully understood. Some results suggested that caspase-1 can promote cellular survival. For example, epithelial cells activate caspase-1 to enhance membrane repair in response to the pore-forming toxins to prevent proteolysis [67]. Whether latently infected T cells can use this caspase-1 pathway to promote their survival remains unknown. Another possibility is through the dysregulated action of IL-7 or IL-15 that can stimulate homeostatic proliferation of latently infected cells. Stromal cells are located in secondary lymph organs such as lymph node trabeculae, lymph vessels, and conduits [68]. IL-7 is observed to be significantly expressed by stromal cells within inflamed lymph nodes [69]. It would be valuable to explore whether HIV-1 can use the caspase-1 pathway to persist in latent cells and whether IL-7 production can be inhibited in the inflamed microenvironment.

The results suggest that cell death by pyroptosis plays an important role in driving slow CD4+ T cell depletion. If pyroptosis can be inhibited, then CD4+ T cells might be maintained. VX-765 is a caspase-1 inhibitor [70–73] used to treat chronic epilepsy and psoriasis. It was found to be safe and well tolerated in humans in a phase IIa trial of epilepsy [74]. Doitsh et al. showed that VX-765 can inhibit secretion of IL-1β and also block cleavage of caspase-1 in HIV-infected tonsillar and splenic lymphoid tissues [23]. However, the active form of VX-765 cannot effectively inhibit cell death, which may be due to reduced cellular permeability [70]. It remains unclear whether the pro-drug VX-765 can efficiently block cell death in vivo. We showed that if antiretroviral drugs cannot effectively block viral replication in lymphoid tissues, then HIV-1 can still establish chronic inflammation in these sites. This is consistent with the observation of persistent inflammation in patients under long-term antiretroviral treatment [75,76]. When drug sanctuary sites exist, CD4+ T cells undergo a very slow depletion or stabilize at a low level (Fig 8). In this case, the immune system would be vulnerable to various opportunistic infections and neoplasms. This may explain the shorter extension of life expectancy in treated patients who had a low CD4+ cell nadir [47–49,51]. If antiretroviral drugs and caspase-1 inhibitors can be effectively delivered to human lymphoid tissues via some transporters [77,78], then CD4+ T cell depletion might be prevented and life expectancy of treated patients might be further extended.

Supporting Information

S1 Text. Additional models and sensitivity test.

https://doi.org/10.1371/journal.pcbi.1004665.s001

(DOCX)

S1 Fig. CD4+ T cells predicted by the one-compartment model with target cell proliferation.

The proliferation rate p is chosen to be 0.001 day-1 (red), 0.01 day-1 (blue), 0.1 day-1 (black), and 1 day-1 (brown) in the simulation.

https://doi.org/10.1371/journal.pcbi.1004665.s002

(EPS)

S2 Fig. Sensitivity test of CD4+ T cells on the cytokine burst size Nc.

The value of Nc increases from 9 molecules (orange line) to 21 molecules (black line) with four equal increments.

https://doi.org/10.1371/journal.pcbi.1004665.s003

(EPS)

S3 Fig. Sensitivity test of CD4+ T cells on the death rate (d3) of abortively infected cells.

The value of d3 increases from 0.0006 day-1 (orange line) to 0.0014 day-1 (black line) with four equal increments.

https://doi.org/10.1371/journal.pcbi.1004665.s004

(EPS)

S4 Fig. Sensitivity test of CD4+ T cells on the cytokine decay rate d5.

The value of d5 increases from 3.96 day-1 (orange line) to 9.24 day-1 (black line) with four equal increments.

https://doi.org/10.1371/journal.pcbi.1004665.s005

(EPS)

S5 Fig. Sensitivity test of CD4+ T cells on the fraction (f) of abortive infection.

The value of f increases from 93.1% (orange line) to 96.9% (black line) with four equal increments.

https://doi.org/10.1371/journal.pcbi.1004665.s006

(EPS)

S6 Fig. Sensitivity test of latently infected cells on the effect (φ) of cytokines promoting latently infected cell proliferation.

The value of φ increases from 0.006 ml molecule-1 (orange line) to 0.014 ml molecule-1 (black line) with four equal increments.

https://doi.org/10.1371/journal.pcbi.1004665.s007

(EPS)

S7 Fig. CD4+ T cells predicted by the two-compartment model with transportation of productively infected cells between compartments.

In the simulation, the value of D1* is fixed to 0.2 day-1 and D2* is fixed to 0.1 day-1.

https://doi.org/10.1371/journal.pcbi.1004665.s008

(EPS)

S8 Fig. Simulation of latently infected cells with different rates of activation aL.

The value of aL increases from 0.01 day-1 (red line) to 0.05 day-1 (orange line) with four equal increments.

https://doi.org/10.1371/journal.pcbi.1004665.s009

(EPS)

Author Contributions

Conceived and designed the experiments: SW LR. Performed the experiments: SW LR. Analyzed the data: SW LR. Contributed reagents/materials/analysis tools: PH MS. Wrote the paper: SW LR.

Reference

  1. 1. Epstein FH, Pantaleo G, Graziosi C, Fauci AS (1993) The immunopathogenesis of human immunodeficiency virus infection. New Engl J Med 328: 327–335. pmid:8093551
  2. 2. Fauci AS, Schnittman SM, Poli G, Koenig S, Pantaleo G (1991) Immunopathogenic mechanisms in human immunodeficiency virus (HIV) infection. Ann Intern Med 114: 678–693. pmid:1672248
  3. 3. Okoye AA, Picker LJ (2013) CD4+ T-cell depletion in HIV infection: mechanisms of immunological failure. Immunol Rev 254: 54–64. pmid:23772614
  4. 4. Quagliarello V (1982) The Acquired Immunodeficiency Syndrome: current status. Yale J Biol Med 55: 443–452. pmid:6134399
  5. 5. Gottlieb MS, Schroff R, Schanker HM, Weisman JD, Fan PT, et al. (1981) Pneumocystis carinii pneumonia and mucosal candidiasis in previously healthy homosexual men: evidence of a new acquired cellular immunodeficiency. New Engl J Med 305: 1425–1431. pmid:6272109
  6. 6. Masur H, Michelis MA, Greene JB, Onorato I, Vande Stouwe RA, et al. (1981) An outbreak of community-acquired Pneumocystis carinii pneumonia: initial manifestation of cellular immune dysfunction. New Engl J Med 305: 1431–1438. pmid:6975437
  7. 7. Muñoz A, Wang M-C, Bass S, Talyor JMG, Kingsley LA, et al. (1989) Acquired immunodeficiency syndrome (AIDS)-free time after human immunodeficiency virus type 1 (HIV-1) seroconversion in homosexual men. Am J Epidemiol 130: 530–539. pmid:2669471
  8. 8. De Boer RJ (2007) Time scales of CD4+ T cell depletion in HIV infection. PLoS Med 4: e193. pmid:17518517
  9. 9. Ho DD, Neumann AU, Perelson AS, Chen W, Leonard JM, et al. (1995) Rapid turnover of plasma virions and CD4 lymphocytes in HIV-1 infection. Nature 373: 123–126. pmid:7816094
  10. 10. Mohri H, Bonhoeffer S, Monard S, Perelson AS, Ho DD (1998) Rapid turnover of T lymphocytes in SIV-infected rhesus macaques. Science 279: 1223–1227. pmid:9469816
  11. 11. McCune JM (2001) The dynamics of CD4+ T-cell depletion in HIV disease. Nature 410: 974–979. pmid:11309627
  12. 12. Okoye A, Meier-Schellersheim M, Brenchley JM, Hagen SI, Walker JM, et al. (2007) Progressive CD4+ central–memory T cell decline results in CD4+ effector–memory insufficiency and overt disease in chronic SIV infection. J Exp Med 204: 2171–2185. pmid:17724130
  13. 13. Zeng M, Haase AT, Schacker TW (2012) Lymphoid tissue structure and HIV-1 infection: life or death for T cells. Trends Immunol 33: 306–314. pmid:22613276
  14. 14. Zeng M, Paiardini M, Engram JC, Beilman GJ, Chipman JG, et al. (2012) Critical role of CD4 T cells in maintaining lymphoid tissue structure for immune cell homeostasis and reconstitution. Blood 120: 1856–1867. pmid:22613799
  15. 15. Matrajt L, Younan PM, Kiem H-P, Schiffer JT (2014) The Majority of CD4+ T-Cell Depletion during Acute Simian-Human Immunodeficiency Virus SHIV89. 6P Infection Occurs in Uninfected Cells. J Virol 88: 3202–3212. pmid:24390339
  16. 16. Mohri H, Perelson AS, Tung K, Ribeiro RM, Ramratnam B, et al. (2001) Increased turnover of T lymphocytes in HIV-1 infection and its reduction by antiretroviral therapy. J Exp Med 194: 1277–1287. pmid:11696593
  17. 17. Kovacs JA, Lempicki RA, Sidorov IA, Adelsberger JW, Herpin B, et al. (2001) Identification of dynamically distinct subpopulations of T lymphocytes that are differentially affected by HIV. J Exp Med 194: 1731–1741. pmid:11748275
  18. 18. Ribeiro RM, Mohri H, Ho DD, Perelson AS (2002) In vivo dynamics of T cell activation, proliferation, and death in HIV-1 infection: why are CD4+ but not CD8+ T cells depleted? Proc Natl Acad Sci U S A 99: 15572–15577. pmid:12434018
  19. 19. Chan ML, Petravic J, Ortiz AM, Engram J, Paiardini M, et al. (2010) Limited CD4+ T cell proliferation leads to preservation of CD4+ T cell counts in SIV-infected sooty mangabeys. Proc Biol Sci 277: 3773–3781. pmid:20591864
  20. 20. Yates A, Stark J, Klein N, Antia R, Callard R (2007) Understanding the slow depletion of memory CD4+ T cells in HIV infection. PLoS Med 4: e177. pmid:17518516
  21. 21. Hernandez-Vargas EA, Middleton RH (2013) Modeling the three stages in HIV infection. J Theor Biol 320: 33–40. pmid:23238280
  22. 22. Miao EA, Rajan JV, Aderem A (2011) Caspase-1-induced pyroptotic cell death. Immunol Rev 243: 206–214. pmid:21884178
  23. 23. Doitsh G, Galloway NL, Geng X, Yang Z, Monroe KM, et al. (2014) Cell death by pyroptosis drives CD4 T-cell depletion in HIV-1 infection. Nature 505: 509–514. pmid:24356306
  24. 24. Doitsh G, Cavrois M, Lassen KG, Zepeda O, Yang Z, et al. (2010) Abortive HIV infection mediates CD4 T cell depletion and inflammation in human lymphoid tissue. Cell 143: 789–801. pmid:21111238
  25. 25. Djomand G, Duerr A, Faulhaber JC, Struchiner CJ, Pacheco AG, et al. (2006) Viral load and CD4 count dynamics after HIV-1 seroconversion in homosexual and bisexual men in Rio de Janeiro, Brazil. J Acquir Immune Defic Syndr 43: 401–404.
  26. 26. Harrison LH, do Lago RF, Friedman RK, Rodrigues J, Santos EM, et al. (1999) Incident HIV infection in a high-risk, homosexual, male cohort in Rio de Janeiro, Brazil. J Acquir Immune Defic Syndr 21: 408–412.
  27. 27. Schechter M, do Lago RF, Mendelsohn AB, Moreira RI, Moulton LH, et al. (2004) Behavioral impact, acceptability, and HIV incidence among homosexual men with access to postexposure chemoprophylaxis for HIV. J Acquir Immune Defic Syndr 35: 519–525.
  28. 28. Fidler S, Porter K, Ewings F, Frater J, Ramjee G, et al. (2013) Short-course antiretroviral therapy in primary HIV infection. New Engl J Med 368: 207–217. pmid:23323897
  29. 29. Lyles RH, Muñoz A, Yamashita TE, Bazmi H, Detels R, et al. (2000) Natural history of human immunodeficiency virus type 1 viremia after seroconversion and proximal to AIDS in a large cohort of homosexual men. J Infect Dis 181: 872–880. pmid:10720507
  30. 30. Wang S, Rong L (2014) Stochastic population switch may explain the latent reservoir stability and intermittent viral blips in HIV patients on suppressive therapy. J Theor Biol 360: 137–148. pmid:25016044
  31. 31. Rong L, Perelson AS (2009) Modeling latently infected cell activation: viral and latent reservoir persistence, and viral blips in HIV-infected patients on potent therapy. PLoS Comput Biol 5: e1000533.
  32. 32. Rong L, Perelson AS (2009) Asymmetric division of activated latently infected cells may explain the decay kinetics of the HIV-1 latent reservoir and intermittent viral blips. Math Biosci 217: 77–87. pmid:18977369
  33. 33. Rong L, Perelson AS (2009) Modeling HIV persistence, the latent reservoir, and viral blips. J Theor Biol 260: 308–331. pmid:19539630
  34. 34. Bofill M, Janossy G, Lee C, MacDonald-Burns D, Phillips A, et al. (1992) Laboratory control values for CD4 and CD8 T lymphocytes. Implications for HIV-1 diagnosis. Clin Exp Immunol 88: 243–252. pmid:1349272
  35. 35. Perelson AS, Kirschner DE, De Boer R (1993) Dynamics of HIV infection of CD4+ T cells. Math Biosci 114: 81–125. pmid:8096155
  36. 36. Markowitz M, Louie M, Hurley A, Sun E, Di Mascio M, et al. (2003) A novel antiviral intervention results in more accurate assessment of human immunodeficiency virus type 1 replication dynamics and T-cell decay in vivo. J Virol 77: 5037–5038. pmid:12663814
  37. 37. Chen HY, Di Mascio M, Perelson AS, Ho DD, Zhang L (2007) Determination of virus burst size in vivo using a single-cycle SIV in rhesus macaques. Proc Natl Acad Sci U S A 104: 19079–19084. pmid:18025463
  38. 38. Hazuda D, Lee J, Young P (1988) The kinetics of interleukin 1 secretion from activated monocytes. Differences between interleukin 1 alpha and interleukin 1 beta. J Biol Chem 263: 8473–8479. pmid:3259579
  39. 39. Efron B, Tibshirani R (1986) Bootstrap methods for standard errors, confidence intervals, and other measures of statistical accuracy. Statist Sci: 54–75.
  40. 40. Silvestri G, Fedanov A, Germon S, Kozyr N, Kaiser WJ, et al. (2005) Divergent host responses during primary simian immunodeficiency virus SIVsm infection of natural sooty mangabey and nonnatural rhesus macaque hosts. J Virol 79: 4043–4054. pmid:15767406
  41. 41. Silvestri G, Sodora DL, Koup RA, Paiardini M, O'Neil SP, et al. (2003) Nonpathogenic SIV infection of sooty mangabeys is characterized by limited bystander immunopathology despite chronic high-level viremia. Immunity 18: 441–452. pmid:12648460
  42. 42. Shedlock DJ, Silvestri G, Weiner DB (2009) Monkeying around with HIV vaccines: using rhesus macaques to define 'gatekeepers' for clinical trials. Nat Rev Immunol 9: 717–728. pmid:19859066
  43. 43. Rosenberg YJ, Janossy G (1999) The importance of lymphocyte trafficking in regulating blood lymphocyte levels during HIV and SIV infections. Semin Immunol 11: 139–154. pmid:10329500
  44. 44. Hengel RL, Jones BM, Kennedy MS, Hubbard MR, McDougal JS (1999) Markers of lymphocyte homing distinguish CD4 T cell subsets that turn over in response to HIV-1 infection in humans. J Immunol 163: 3539–3548. pmid:10477629
  45. 45. Kirschner D, Webb G, Cloyd M (2000) Model of HIV-1 disease progression based on virus-induced lymph node homing and homing-induced apoptosis of CD4+ lymphocytes. J Acquir Immune Defic Syndr 24: 352–362.
  46. 46. Embretson J, Zupancic M, Ribas JL, Burke A, Racz P, et al. (1993) Massive covert infection of helper T lymphocytes and macrophages by HIV during the incubation period of AIDS. Nature 362: 359–362. pmid:8096068
  47. 47. Collaboration ATC (2008) Life expectancy of individuals on combination antiretroviral therapy in high-income countries: a collaborative analysis of 14 cohort studies. Lancet 372: 293–299. pmid:18657708
  48. 48. Fang C, Chang Y, Hsu H, Twu S, Chen K-T, et al. (2007) Life expectancy of patients with newly-diagnosed HIV infection in the era of highly active antiretroviral therapy. QJM 100: 97–105. pmid:17277317
  49. 49. Lloyd-Smith E, Brodkin E, Wood E, Kerr T, Tyndall MW, et al. (2006) Impact of HAART and injection drug use on life expectancy of two HIV-positive cohorts in British Columbia. AIDS 20: 445–450. pmid:16439879
  50. 50. May MT, Gompels M, Delpech V, Porter K, Orkin C, et al. (2014) Impact on life expectancy of HIV-1 positive individuals of CD4+ cell count and viral load response to antiretroviral therapy. AIDS 28: 1193–1202.
  51. 51. Mills EJ, Bakanda C, Birungi J, Chan K, Ford N, et al. (2011) Life expectancy of persons receiving combination antiretroviral therapy in low-income countries: a cohort analysis from Uganda. Ann Intern Med 155: 209–216. pmid:21768555
  52. 52. Geginat J, Sallusto F, Lanzavecchia A (2001) Cytokine-driven proliferation and differentiation of human naive, central memory, and effector memory CD4+ T cells. J Exp Med 194: 1711–1719. pmid:11748273
  53. 53. Bikker A, Erik Hack C, PJG Lafeber F, AG van Roon J (2012) Interleukin-7: a key mediator in T cell-driven autoimmunity, inflammation, and tissue destruction. Curr Pharm Des 18: 2347–2356. pmid:22390698
  54. 54. Hartgring S, Willis CR, Bijlsma J, Lafeber F, van Roon J (2012) Interleukin-7-aggravated joint inflammation and tissue destruction in collagen-induced arthritis is associated with T-cell and B-cell activation. Arthritis Res Ther 14: R137. pmid:22676399
  55. 55. Haase A (1999) Population biology of HIV-1 infection: viral and CD4+ T cell demographics and dynamics in lymphatic tissues. Annu Rev Immunol 17: 625–656. pmid:10358770
  56. 56. Nowak MA, Lloyd AL, Vasquez GM, Wiltrout TA, Wahl LM, et al. (1997) Viral dynamics of primary viremia and antiretroviral therapy in simian immunodeficiency virus infection. J Virol 71: 7518–7525. pmid:9311831
  57. 57. Chun T-W, Carruth L, Finzi D, Shen X, DiGiuseppe JA, et al. (1997) Quantification of latent tissue reservoirs and total body viral load in HIV-1 infection. Nature 387: 183–188.
  58. 58. Harper ME, Marselle LM, Gallo RC, Wong-Staal F (1986) Detection of lymphocytes expressing human T-lymphotropic virus type III in lymph nodes and peripheral blood from infected individuals by in situ hybridization. Proc Natl Acad Sci U S A 83: 772–776. pmid:3003749
  59. 59. Nardelli B, Gonzalez CJ, Schechter M, Valentine FT (1995) CD4+ blood lymphocytes are rapidly killed in vitro by contact with autologous human immunodeficiency virus-infected cells. Proc Natl Acad Sci U S A 92: 7312–7316. pmid:7638187
  60. 60. Sousa AE, Carneiro J, Meier-Schellersheim M, Grossman Z, Victorino RM (2002) CD4 T cell depletion is linked directly to immune activation in the pathogenesis of HIV-1 and HIV-2 but only indirectly to the viral load. J Immunol 169: 3400–3406. pmid:12218162
  61. 61. Fleury S, Rizzardi G, Chapuis A, Tambussi G, Knabenhans C, et al. (2000) Long-term kinetics of T cell production in HIV-infected subjects treated with highly active antiretroviral therapy. Proc Natl Acad Sci U S A 97: 5393–5398. pmid:10805798
  62. 62. Hellerstein M, Hanley M, Cesar D, Siler S, Papageorgopoulos C, et al. (1999) Directly measured kinetics of circulating T lymphocytes in normal and HIV-1-infected humans. Nat Med 5: 83–89. pmid:9883844
  63. 63. Nash D, Katyal M, Brinkhof MW, Keiser O, May M, et al. (2008) Long-term immunologic response to antiretroviral therapy in low-income countries: a collaborative analysis of prospective studies. AIDS 22: 2291–2302. pmid:18981768
  64. 64. Osborn L (1997) Endothelium and inflammation: Molecules and therapeutic prospects. In: Hamann A, editor. Adhesion Molecules and Chemokines in Lymphocyte Trafficking Harwood. Academic Publishers. pp. 217–236.
    • 65. Epstein FH, Gabay C, Kushner I (1999) Acute-phase proteins and other systemic responses to inflammation. New Engl J Med 340: 448–454. pmid:9971870
    • 66. Kunkel EJ, Butcher EC (2002) Chemokines and the tissue-specific migration of lymphocytes. Immunity 16: 1–4. pmid:11825560
    • 67. Gurcel L, Abrami L, Girardin S, Tschopp J, van der Goot FG (2006) Caspase-1 activation of lipid metabolic pathways in response to bacterial pore-forming toxins promotes cell survival. Cell 126: 1135–1145. pmid:16990137
    • 68. Onder L, Narang P, Scandella E, Chai Q, Iolyeva M, et al. (2012) IL-7–producing stromal cells are critical for lymph node remodeling. Blood 120: 4675–4683. pmid:22955921
    • 69. Malhotra D, Fletcher AL, Astarita J, Lukacs-Kornek V, Tayalia P, et al. (2012) Transcriptional profiling of stroma from inflamed and resting lymph nodes defines immunological hallmarks. Nat Immunol 13: 499–510. pmid:22466668
    • 70. Boxer MB, Quinn AM, Shen M, Jadhav A, Leister W, et al. (2010) A highly potent and selective caspase 1 inhibitor that utilizes a key 3-cyanopropanoic acid moiety. ChemMedChem 5: 730–738. pmid:20229566
    • 71. Randle JC, Harding MW, Ku G, Schönharting M, Kurrle R (2001) ICE/Caspase-1 inhibitors as novel anti-inflammatory drugs. Expert Opin Investig Drugs 10: 1207–1209. pmid:11772244
    • 72. Stack JH, Beaumont K, Larsen PD, Straley KS, Henkel GW, et al. (2005) IL-converting enzyme/caspase-1 inhibitor VX-765 blocks the hypersensitive response to an inflammatory stimulus in monocytes from familial cold autoinflammatory syndrome patients. J Immunol 175: 2630–2634. pmid:16081838
    • 73. Maroso M, Balosso S, Ravizza T, Iori V, Wright CI, et al. (2011) Interleukin-1β biosynthesis inhibition reduces acute seizures and drug resistant chronic epileptic activity in mice. Neurotherapeutics 8: 304–315. pmid:21431948
    • 74. Vezzani A, Balosso S, Maroso M, Zardoni D, Noé F, et al. (2010) ICE/caspase 1 inhibitors and IL-1beta receptor antagonists as potential therapeutics in epilepsy. Curr Opin Investig Drugs 11: 43–50.
    • 75. Fletcher CV, Staskus K, Wietgrefe SW, Rothenberger M, Reilly C, et al. (2014) Persistent HIV-1 replication is associated with lower antiretroviral drug concentrations in lymphatic tissues. Proc Natl Acad Sci U S A 111: 2307–2312. pmid:24469825
    • 76. Rönsholt FF, Ullum H, Katzenstein TL, Gerstoft J, Ostrowski SR (2013) Persistent inflammation and endothelial activation in HIV-1 infected patients after 12 years of antiretroviral therapy. PLoS ONE 8: e65182. pmid:23755191
    • 77. Sakaeda T, Nakamura T, Okumura K (2004) Pharmacogenetics of drug transporters and its impact on the pharmacotherapy. Curr Top Med Chem 4: 1383–1398.
    • 78. Löscher W, Potschka H (2005) Role of drug efflux transporters in the brain for drug disposition and treatment of brain diseases. Prog Neurobiol 76: 22–76. pmid:16011870