Modeling the within-host dynamics of HIV infection
© Perelson and Ribeiro; licensee BioMed Central Ltd. 2013
Received: 2 April 2013
Accepted: 2 September 2013
Published: 3 September 2013
Skip to main content
© Perelson and Ribeiro; licensee BioMed Central Ltd. 2013
Received: 2 April 2013
Accepted: 2 September 2013
Published: 3 September 2013
The new field of viral dynamics, based on within-host modeling of viral infections, began with models of human immunodeficiency virus (HIV), but now includes many viral infections. Here we review developments in HIV modeling, emphasizing quantitative findings about HIV biology uncovered by studying acute infection, the response to drug therapy and the rate of generation of HIV variants that escape immune responses. We show how modeling has revealed many dynamical features of HIV infection and how it may provide insight into the ultimate cure for this infection.
Since the discovery of HIV as the etiological agent of AIDS, numerous advances have been made in our understanding of the molecular biology, pathogenesis, and epidemiology of the virus, and the host immune response to it. Not least among these has been the knowledge obtained by mathematical analysis and within-host modeling of changes in viral load and T-cell counts after initiation of potent antiretroviral therapy in individual subjects. Indeed, modeling of the kinetics of HIV RNA under drug therapy has led to substantial insights into the dynamics and pathogenesis of HIV-1 [1–6] and the existence of multiple reservoirs that have made eradication of the virus difficult [7, 8]. Through these analyses it has been possible to quantify the rapidity of HIV infection and replication, the rate of virion clearance, the lifespan of productively infected cells [1, 2, 4, 5, 9, 10], and to predict the impact of treatment and the appearance of drug-resistant variants [11–13]. Other modeling efforts have helped clarify controversial issues relating to the mechanism of T-cell depletion in HIV infection  and motivated new experimental and clinical studies. More recent modeling studies have addressed issues such as immune escape and viral evolution, allowing a window into the quantification of the immune mechanisms operating in the setting of HIV infection.
Below we briefly review how quantitative data and modeling have contributed to the understanding of HIV biology.
Only a small fraction of CD4+ T cells in the periphery become infected with HIV  and thus identifying the target cells in this model is not straightforward. However, the model is able to describe the kinetics of T-cell depletion in macaques infected with an X4-tropic virus, such as SHIV89.6P, where most T cells are target cells . Despite its ability to fit data, this model may be too simple in that it does not include any explicit immune response (more on this below). Nonetheless, this simple model and generalizations that include long-lived infected cells and latently infected cells (Figure 1) have proven to be useful and have generated important insights into the biology of HIV [4, 9, 20].
If one follows the response to combination antiretroviral therapy (cART) for more than a week or two, one sees that the rapid first phase of decay of plasma viremia is followed by a slower second phase of decay (Figure 2). This second phase has been attributed to the existence of longer-lived productively infected cells, perhaps resting CD4+ T cells or cells of the macrophage/monocyte lineage . Indeed, there is clear heterogeneity in the cell types that are infected by HIV and in the amount of virus produced by these cell types , consistent with the suggestion that some infected cells may live substantially longer than others. Nonetheless, some modelers have suggested other explanations for this second phase, such as the decline being driven by cytotoxic T lymphocytes (CTLs), which slows as the CTL response declines , or that infected cells have an age-dependent transactivation rate, which slows the generation of virus-producing cells .
Irrespective of the mechanism generating the second phase, with continued cART viral levels decline below the detection limit of clinical assays (50 HIV RNA copies/ml), and with these assays one cannot determine how long the second phase lasts.
Based on the rate of second phase decay, modeling suggested that 3 to 4 years of fully suppressive therapy could eliminate the cells responsible for second phase virus . However, the use of single copy assays (SCAs), which allow one to detect as few as one HIV RNA/ml of plasma , has led to the identification of a third phase of decay with a half-life estimated as 39 weeks in one study  and 69 weeks in another  with overlapping and rather large 95% confidence intervals. Palmer et al., as well as others [29, 32], also suggested a fourth phase or constant level of viremia is attained after very long times on therapy (Figure 2). The existence of these later phases implies that it may not be possible to eliminate HIV with antiviral therapy alone. The sources of third and fourth phase viremia are controversial and may include leakage of bound virus from follicular dendritic cells (FDCs) [33–35], the release of virus from latently infected cells and ongoing viral replication.
HIV incorporates into the host cell genome and can establish a latent form of HIV infection, involving a small fraction of resting memory CD4+ T cells that carry integrated viral genomes [18, 36, 37]. Longitudinal analysis suggested that this latent reservoir could persist in patients for as long as 60 years . While latently infected cells do not produce virus in the resting state, they can do so upon activation . This feature has been explored by modelers to explain the occasional viral ‘blips’ seen in patients who are otherwise well-suppressed [38–42], as well as the low levels of plasma viremia [30, 41, 43–46] detected with research assays that have lower limits of detection in the order of one HIV RNA/ml . Since the latent pool is not depleted by this occasional activation, modeling suggested the possibility that latently infected cells undergo homeostatic proliferation without activation and thus maintain the latent reservoir for decades [41, 45]. This hypothesis was later experimentally confirmed .
Modeling has also examined the role of early therapy in preventing the establishment of latency  and work is underway by both experimental and modeling groups examining the possibility of activating latently infected cells by therapies, such as the use of histone deacetylase inhibitors, to potentially generate a cure for HIV [49, 50]. In assessing such new therapies, models for the activation of latently infected cells [5, 41, 46] may provide critical insights.
So far, we have described the ‘ideal’ scenario for antiretroviral therapy, where treatment is successful in driving the viral load below detection levels. However, viral load sometimes rebounds and new viral strains with mutations that confer resistance to the drugs used in the treatment protocol are observed.
When HIV infects a cell, the viral RNA genome is reverse transcribed into DNA. The reverse transcription process is error prone and results in mutation at an estimated rate of 3 × 10-5 per base per generation, with about two-thirds of these mutations being nucleotide substitutions . Thus, when HIV with its almost 10 kb genome is reverse transcribed there is about a 20% chance that a base substitution occurs. In a chronically infected patient it has been estimated that about 108 cells are infected each day, thus allowing all possible single and many double mutations to be explored each day [3, 11]. Also, drugs may not penetrate all tissues and all cells with equal efficiency, and drug ‘sanctuaries’ may exist . Thus, it is not surprising that resistance to antiretroviral drugs can develop relatively quickly, if combination therapy is not employed. This was explained by modeling within-host viral dynamics [12, 13, 53–58]. Moreover, several quantitative studies indicated that resistance is more likely to pre-exist before therapy is started than to appear de novo during the therapy, especially if therapy is strong enough to curtail replication [12, 54, 59]. Three factors contribute to the development of resistance, and indeed to the observed diversity of HIV: i) its high mutation rate, which is typical of an RNA virus, ii) its fast lifecycle, and iii) the long-term nature of HIV infection. These three factors combine to allow rapid viral evolution and the generation of high diversity.
More recently, studies have shown that some patients derive clinical benefit from continued therapy even when the virus is resistant to the drug protocol, as the resistant virus can be less fit than wild type .
Even though the introduction of antiretroviral treatment was a great success, and improved the quality of life of countless people, it has not become the panacea that was expected. For several reasons, including those alluded to above, treatment has not yet cured HIV. Thus, it has become clear that the development of a vaccine that prevents infection in the first place or ameliorates the course of disease is very important. To this end a better understanding of the clinical progression of HIV and the immune response against it is crucial. Here too, models have been playing a critical role as we discuss below.
HIV is enormously diverse . Thus, it was very surprising to discover that approximately 80% of sexually transmitted infections are the result of a single transmitted/founder virus . This discovery was a perfect example of extraordinary technical developments in assay capabilities pushing the state-of-the-art in modeling to analyze the results generated by those assays. In this case, single genome amplification (SGA) assays allowed an unprecedented look at the phylogenetic structure of HIV early post-infection. To analyze these data, a stochastic model of viral diversity generation was developed [63, 64]. The results conclusively showed that the observed early genetic structure was most compatible with a single virus being transmitted or founding the complete viral population in a large fraction of heterosexual transmissions.
Not every encounter with HIV/SIV results in infection; the transmission probability for HIV has been estimated as 10-3 to 10-2 per coital act . This suggests that when low levels of virus are transmitted the infection may go extinct rather than take off. This possibility, as well as the potential efficacy of pre- and post-exposure prophylaxis, has been explored though the use of stochastic models of early infection [74, 75]. These models suggest that extinction would most likely occur well before virus is detectable even with single copy assays. Thus, the low observed HIV-1 transmission rate may be a consequence of small numbers of virions being transmitted followed by frequent extinction.
Although the modeling studies mentioned so far, dealing with a variety of different processes in HIV infection, have been quite successful in advancing our understanding of HIV biology, a puzzling aspect remains: these models do not necessarily include an explicit immune response.
The basic model given by equation (1), lacks an explicit representation of the immune response and hence has been called a target-cell limited model . Despite this deficit, the model fits viral kinetic data obtained both during natural infection (Figure 3) and while patients are on therapy. However, immune responses against HIV, at both the cellular and humoral level, can be detected and there are considerable data indicating a role for CD8+ T-cell responses during HIV infection, particularly in people called elite controllers [77–79]. Further, depleting CD8+ T cells during acute SIV infection is associated with SIV remaining at high levels rather than reaching a distinct peak and then rapidly falling . Also, HIV-1 tends to start accumulating CTL escape mutations around the time of the peak viremia, supporting the notion that CD8+ T cells play an important role in controlling early viremia . On the other hand, increasing the initial number of HIV-specific CD8+ T cells, by vaccinating macaques prior to infection, did not change the growth rate or decay rate of virus from its peak during primary infection, suggesting a limited role of CD8+ T cells [82, 83]. How to reconcile these various observations and account for the initial immune response in models of acute infection is still a subject of research and debate [84, 85].
While many models have included CD8 responses [68, 86], more than we can review here, they tend to lack comparisons with experimental data leaving the field without good estimates for the parameters that govern CTL effects. An interesting example of these types of models was developed to analyze experiments of SIV in macaques where early drug treatment led to control of the virus in the long term (akin to ‘elite controllers’) [87, 88]. In these models, the interactions between virus, CD4+ T cells and CD8+ T cells were considered , as HIV-induced depletion of CD4+ T cells may affect one’s ability to mount CD8+ T-cell responses .
Modeling antibody responses to HIV is still an undeveloped area, although some work involving data interpretation has been done [91, 92]. Interestingly, in modeling acute SIV infection an improvement in fit to the measured viral loads was attained by allowing the viral infectivity, β in equation (1), to decrease with time, possibly reflecting the effects of antibody that builds up over time in reducing viral infectivity . This effect was suggested by the work of Ma et al. that showed mixing set-point plasma with acute-phase plasma decreased the infectivity of the acute-phase plasma.
While models with only two viral species are easy to analyze, simulation models have looked at much more complex situations in which there are multiple escapes . Further, recent data using single genome amplification and sequencing shows that HIV can escape at multiple T-cell epitopes [81, 102]. Thus, models are being developed that consider the entire HIV genotype and not simply escape at single epitopes that are treated independently .
Some studies quantified the efficiency of cytotoxic T lymphocytes in killing infected cells based on the rate of escape of viral variants at specific epitopes. These results were used to compare the effect of immune pressure in macaques infected with SIV with that of humans . However, precise estimates of escape rates were hampered by a lack of frequent sampling. Recently, a large effort has been expended by the Center for HIV AIDS Vaccine Immunology (CHAVI) to elucidate the early immune response during primary infection. Identifying subjects with a single transmitted/founder sequence and then following the evolution of that sequence in time using either deep sequencing  or single genome amplification techniques has made it possible to study the dynamics of the emergence of both antibody and CTL escape variants in some detail [81, 92, 102]. While early results suggested CTL response against a single epitope only provides a modest amount of pressure , later work using much more frequent patient sampling showed that the CTL response against a single epitope could account for as much as 35% of the killing .
Balamurali et al. measured the post-peak maximal viral decay rates of wild-type and escape mutant virus in macaques and found them to be the same, suggesting that CD8+ T cells may act by non-cytolytic mechanisms. We  and others  tried to address this issue directly by modeling data on the decay rate of virus during cART in the presence and absence of CD8+ T cells . Our results were consistent with CD8+ T cells mainly killing cells before they began producing virus or with CD8s mainly having a non-cytolytic effect, but a later analysis showed a small cytolytic effect on productively infected cells could not be ruled out . Experiments that determine the rate of viral decay during cART only provide information about the fate of productively infected cells. Thus, other modeling studies have examined the possibility that CTLs act by killing cells before they become productively infected [110, 111]. Studies of escape have also been used in other imaginative ways - for example, to estimate the turnover of integrated DNA in resting CD4+ T cells, which is one of the blocks in the elimination of infection .
As AIDS is characterized by a loss of CD4+ T cells, much modeling work has focused on quantifying T-cell turnover during HIV infection using direct labeling as well as cell markers such as T-cell receptor excision circles. Much of this modeling literature has recently been reviewed  and thus will not be discussed here.
Modeling HIV dynamics has been a rich area of study that, we believe, has spearheaded the wider field of modeling in viral and immune system dynamics. Indeed, there is much more research in HIV modeling than we could possibly cover here, but we hope to have given a representative flavor of the most innovative studies.
Looking into the future, there are several areas where we believe modeling can still make important contributions. Clearly, as more quantitative data on the immune response, both cellular and humoral, become available, it will be important to include these in the mechanistic models of HIV infection. This will be especially important to explain the events during primary infection, and in the analysis of vaccine trial data in humans. These models should be able to describe not just the time evolution of the viral load, as do current models, but also the kinetics of the immune response and its effect on virus and infected cells. Another area of interest to modelers is the potential effect of cell-to-cell transmission of HIV-1 [114–119]. If a sizeable proportion of infections in vivo occur this way, the virus may be able to avoid antibody effects, which mostly affect free virus, and thus evade the effect of the humoral immune response. Also, the transfer of multiple viral genomes by cell-to-cell transmission can reduce the efficiency of cART . Finally, as efforts intensify on finding therapies to activate latent cells, models quantifying the effects of such therapies and the balance of re-activation and new infections could help inform best protocols for the clinic.
Overall, we believe that two of the most pressing questions in the HIV field are why primary infection invariably leads to chronicity, and why chronic infection leads to persistence of the virus, even in the face of very potent and durable therapy. It is unquestionable that future modeling studies of novel datasets will help us to understand both of these questions, proposing different mechanisms and hypotheses for these observations. In so doing, we might uncover new intervention strategies to help prevent or eradicate infection. The idea of curing HIV infection either, by viral eradication or functionally curing the infection by controlling it, is gaining traction [121, 122]. Modeling will surely play a role in these endeavors, and in this way modeling will have come full circle from early insights into viral biology to demonstrating that cure is possible.
Mathematical analysis of HIV-1 viral dynamics and immune responses has led to a number of important insights about the dynamics and pathogenesis of HIV infection. Modeling plasma virus decay under therapy demonstrated the fast turnover of virus, explaining the potential for generation of mutants and the development of drug resistance. This early work paved the way for many collaborations between clinicians and modelers to understand the nature of the multiphasic viral decay observed in treated patients, the initial expansion of virus upon infection, the turnover of CD4+ and CD8+ T cells, the probability of single transmitted/founder viruses and many others. It is fair to say that most, if not all, of these insights would not have been possible without close interdisciplinary collaborations allowing quantitative modeling.
Portions of this work were performed under the auspices of the US Department of Energy under contract DE-AC52-06NA25396 and supported by NIH grants P20-GM103452, OD011095, AI028433, and the Center for HIV/AIDS Vaccine Immunology and Immunogen Discovery grant number UM1-AI100645-01. RMR received partial funding from the EU 7th Framework Program under grant no. PCOFUND-GA-2009-246542 and from FCT Portugal and part of his work was done while visiting Instituto de Medicina Molecular, Faculdade de Medicina da Universidade de Lisboa, Lisbon, Portugal.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.