A model of HIV drug resistance driven by heterogeneities in host immunity and adherence patterns
- Anna Bershteyn^{1}Email author and
- Philip A Eckhoff^{1}
https://doi.org/10.1186/1752-0509-7-11
© Bershteyn and Eckhoff; licensee BioMed Central Ltd. 2013
Received: 5 June 2012
Accepted: 16 January 2013
Published: 4 February 2013
Abstract
Background
Population transmission models of antiretroviral therapy (ART) and pre-exposure prophylaxis (PrEP) use simplistic assumptions – typically constant, homogeneous rates – to represent the short-term risk and long-term effects of drug resistance. In contrast, within-host models of drug resistance allow for more detailed dynamics of host immunity, latent reservoirs of virus, and drug PK/PD. Bridging these two levels of modeling detail requires an understanding of the “levers” – model parameters or combinations thereof – that change only one independent observable at a time. Using the example of accidental tenofovir-based pre-exposure prophyaxis (PrEP) use during HIV infection, we will explore methods of implementing host heterogeneities and their long-term effects on drug resistance.
Results
We combined and extended existing models of virus dynamics by incorporating pharmacokinetics, pharmacodynamics, and adherence behavior. We identified two “levers” associated with the host immune pressure against the virus, which can be used to independently modify the setpoint viral load and the shape of the acute phase viral load peak. We propose parameter relationships that can explain differences in acute and setpoint viral load among hosts, and demonstrate their influence on the rates of emergence and reversion of drug resistance. The importance of these dynamics is illustrated by modeling long-lived latent reservoirs of virus, through which past intervals of drug resistance can lead to failure of suppressive drug regimens. Finally, we analyze assumptions about temporal patterns of drug adherence and their impact on resistance dynamics, finding that with the same overall level of adherence, the dwell times in drug-adherent versus not-adherent states can alter the levels of drug-resistant virus incorporated into latent reservoirs.
Conclusions
We have shown how a diverse range of observable viral load trajectories can be produced from a basic model of virus dynamics using immunity-related “levers”. Immune pressure, in turn, influences the dynamics of drug resistance, with increased immune activity delaying drug resistance and driving more rapid return to dominance of drug-susceptible virus after drug cessation. Both immune pressure and patterns of drug adherence influence the long-term risk of drug resistance. In the case of accidental PrEP use during infection, rapid transitions between adherence states and/or weak immunity fortifies the “memory” of previous PrEP exposure, increasing the risk of future drug resistance. This model framework provides a means for analyzing individual-level risks of drug resistance and implementing heterogeneities among hosts, thereby achieving a crucial prerequisite for improving population-level models of drug resistance.
Background
Quantitative analysis of the risk of HIV drug resistance is important on both an individual and a population level, especially when the “real-world” use of drugs may differ significantly from the more ideal setting of randomized controlled trials. Examples of non-ideal “real-world” situations include poor adherence, late entry into therapy, drug stockouts, unauthorized re-distribution of antiretrovirals, and interruption of multi-drug regimens in which one drug has a longer half-life, which creates an interval of “effective monotherapy” when only this drug is present at high levels [1].
Interest in models of drug resistance has increased recently due to the addition of pre-exposure prophylaxis (PrEP) to the world’s HIV prevention toolkit [2]. Clinical trials have demonstrated that PrEP with tenofovir disoproxil fumarate (TDF), combined TDF and emtricitabine (FTC), and tenofovir vaginal gel can reduce the risk of HIV acquisition [3–6]. No PrEP regimen prevented transmission entirely, though increased adherence correlated with increased protection against HIV. This creates a risk of accidental PrEP use during “breakthough” infections, until the individual is diagnosed with HIV and PrEP is discontinued. Additionally, there is risk of accidental initiation of PrEP by infected individuals due to faulty testing or early “window period” testing. In addition to their use in PrEP, TDF and FTC are also found in two popular single-pill, once-daily fixed-dose combination therapies (Atripla and Complera) as well as the fixed-dose “quad” pill recently approved by the FDA. Thus, there is concern as to whether resistance caused by PrEP could threaten the ability to use two of the safest available drugs available in a convenient, single-dose, once-daily regimen.
Mathematical models of ART and PrEP have been used to assess the risk of drug resistance on the individual and population level [7]. However, state-of-the-art population models have failed to capture heterogeneities in the risk of drug resistance among individuals, due to the disparity in model detail between population-level and within-host models.
Population-level models of PrEP and resistance tend to focus on specific conditions of HIV transmission, with a majority of recent oral PrEP models focusing on heterosexual generalized epidemics [8–10], and others on concentrated epidemics among men who have sex with men [11]. Of the small number of population PrEP models that include drug resistance, only one has considered reversion of resistance by assuming that resistance and reversion occur at fixed rates for the treated and post-treated subpopulations, respectively [10, 11]. Recent models of combination ART that account for resistance and reversion treat these similarly [12, 13]; though one model used stochastic numerical methods as a proxy for variability in acquisition and transmission of resistance, these are still assumed to occur at a fixed rate [13] with no consideration for heterogeneous host biology and behavior.
Within-host models have provided some mechanistic insights into HIV progression, drug efficacy, and the risk of resistance. Nowak and Bangham formulated a model of virus dynamics that includes mutation [14], and has inspired dozens of variations, such as mechanisms of viral escape arising from competition among quasispecies-specific cytotoxic T lymphocyte (CTL) [15] and antibody [16] responses. A common application of these models is predicting the relationship between dosage and acquisition of resistance: a quantity not well-characterized in clinical studies [17]. Very recently, a detailed within-host model similar to the one presented here was used to propose a mechanism for drug resistance during HAART that includes a protease inhibitor, proposing a relationship between adherence level and risk of treatment failure [18]. This and other models [19] have assumed that adherence is primarily characterized by the percentage, not the pattern, of doses taken. Other models have incorporated time-correlated field data [20] and to explore the impact of “clustering” dropped doses on the drug resistance outcome [21]. Measuring the true patterns of drug adherence can be challenging [22–25], but models can represent a range of possible behaviors in silico to understand the effects of different possible patterns.
The host immune response is another source of heterogeneity that, as we will show, contributes substantially to the dynamics of resistance and reversion. Though a wealth of biological studies have revealed the complex relationship between immunity and virologic control, few models have explored the effect of viral replication and immunity on drug efficacy [26, 27].
Ultimately, population models must capture the important outcomes of within-host heterogeneities and react appropriately to a range of assumptions about viral fitness, host immune response, drug pharmacology, and adherence. Here, we use the example of tenofovir-based PrEP to explore the insights that can be gained from considering such heterogeneities in a combined model of within-host virus dynamics, drug pharmacokinetics/pharmacodynamics (PK/PD), and patterns of drug adherence.
Methods
We constructed a model that combines the pharmacokinetic/pharmacodynamic outcome of time-varying drug adherence with the dynamics of competition between drug-resistant mutants. Hypotheses or data about drug adherence patterns are provided as input for the simulation, and the model predicts the time-varying populations of CD4+ T cells and virus (WT and drug-resistant) in the form of plasma RNA and integrated DNA. The model is deterministic, but can be driven stochastically through the choice of adherence pattern.
Next, a detailed pharmacokinetics model is used to calculate the time-varying concentration of the active form of the drug in the relevant compartment, given the sequence of doses that were taken over time. For oral TDF, the model includes bioavailability, partitioning across the cell membrane, and the two phosphorylation steps required to create the active diphosphate form of the drug.
A pharmacodynamics model is used to translate the drug concentration into an effect on the replicative capacity of the WT and mutant viral quasispecies.
Finally, a dynamic model of viral replication responds to the changes in replicative capacity caused by fluctuating drug concentrations. Here, uninfected CD4+ target cells (denoted x) are infected by a given viral quasispecies (v_{ i } for quasispecies i, e.g., v_{ WT } and v_{ mutant }) to form a productively infected cell (y_{ i }). Infected cells can be eliminated by virus-specific CTLs (z). In an alternative version of this model (Figure 1b), infection can yield long-lived latently infected cells (w_{ i }) with probability g, or productively infected cells (y_{ i }) with probability 1-g. Latently infected cells produce no virus and are not susceptible to CTLs until they become activated into productively infected cells at a slow rate, h.
Model parameters
Parameter values
Parameter | Value | Units | Ref | Parameter | Value | Units | Ref |
v _{ SS } | 51,000 | copies mL ^{ -1 } | [32] | m _{ WT } | 1.139 | Unitless | [46] |
λ | 10^{5} | cells mL^{ -1 } day^{ -1 } | [32] | m _{ K65R } | 1.126 | Unitless | [46] |
b | 0.027 | day ^{ -1 } | [32] | IC _{ 50 ,WT } | 0.050* | μg mL ^{ -1 } | [46] |
d | 0.1 | day ^{ -1 } | [32] | IC _{ 50 ,K65R } | 2.2* | μg mL ^{ -1 } | [46] |
u | 5 | day ^{ -1 } | [32] | K _{ a } | 14.64 | day ^{ -1 } | [41] |
a | 0.5 | day ^{ -1 } | [32] | K _{ e } | 9.6 | day ^{ -1 } | [41] |
c _{ 1 } | 1000 | cells mL^{ -1 } day^{ -1 } | [32] | K _{ 1f } | 9.6 | day ^{ -1 } | [41] |
β | 2f × 10^{5} | mL count^{ -1 } day^{ -1 } | [32] | K _{ 1b } | 30.3 | day ^{ -1 } | [41] |
f | Variable (Equation 9) | Unitless | K _{ 2f } | 270.7 | day ^{ -1 } | [41] | |
h | 0.002 | day ^{ -1 } | [41] | K _{ 2b } | 95.5 | day ^{ -1 } | [41] |
w _{ ss } /y _{ ss } (for g ) | 0.39 | Unitless | [72] | K _{ m } | 24000 | day ^{ -1 } | [41] |
k _{ WT } | 100/f | day ^{ -1 } | [32] | K _{ c } | 1.1 | day ^{ -1 } | [41] |
k _{ K65R } | 61.054/f | day ^{ -1 } | [46] | f _{ plasma-bound } | 0.07 | Unitless | [41] |
H | 1800 | Unitless | [42] | * Concentration in plasma. |
Adherence patterns
A hypothesis or data set related to adherence is converted to a bit string representing whether the dose is taken or missed at each time in the dosing schedule. The sequence could potentially be informed by field data, e.g., if a time-correlated data set such as electronic monitoring of bottle opening were used [24, 25].
where J_{ i } represents the set of possible states that could be entered from i, r_{i → j} is the rate of transitioning from state i to state j, and rand is a random number uniformly distributed between 0 and 1. After time ∆t, the next state was chosen among possible states J_{ i } proportionally to the relative magnitude of the transition rates r_{i → j}. Finally, the states were discretized into a binary sequence of taken or missed doses by determining whether the individual was in an ART-taking or ART-missing state at the time of each dose.
Virus dynamics and immune response
We assumed that the population of free viruses and cell-incorporated viruses interact with a common pool of uninfected target cells, and that immune pressure p was exerted equally on all infected cells. (This assumption would break down if the drug resistance mutation occurred in an epitope targeted by two CTLs that may be present at significantly different levels or subject to different levels of immune regulation.) The steady-state frequency of HIV-specific CTLs z_{ss} = c/b was set to 37 cells/uL, consistent with phenotypic analysis of HIV-infected human peripheral blood [35, 36]. We assumed z grows at rate c independently of z, y, or the product zy, the dynamics of which are described elsewhere [32–34]. We explored this system under varying p because this is a known heterogeneity among hosts: viremia is elevated [37] or suppressed [38] in the presence of certain HLA-restricted CTL subtypes.
whereas the growth of target cells, free virions, and CTLs were still governed by Equations 2, 4, and 5, respectively.
We investigated the effect of changing p with and without a compensating change in f that preserves v_{ ss } at the mean value measured in discordant heterosexual couples in Zambia, 51,000 copies/mL [39].
In the latency simulations shown, we set the latently infected cell reactivation rate h to 0.002 days^{-1}, representing the long-lived component of the latently infected cell reservoir [41]. The simplifying assumptions made by this implementation of latency include (1) a constant low rate of gradual reactivation, in lieu of the wide range of re-activation rates possible from different cell subtypes, and (2) the production of a single activated cell by each latently infected cell, neglecting proliferation in the process of re-activation.
Pharmacokinetics
The rate constants for this model, also taken from Dixit and Perelson [41], are listed in Table 1.
Pharmacodynamics
The pharmacodynamics component of the model takes as input the time-varying drug concentration in its active form (C_{ cpp }) and determines its time-varying effect on viral replication. We modeled only the single mutants known to confer drug resistance using recently reported dose–response data from single-cycle replication studies [46]. Mutation is assumed to incur a fitness penalty that reduces the virus production rate k. In the case of the K65R mutant, this rate is reduced by 39% at zero drug concentration. Multiple-position mutants were not modeled.
where IC_{ 50 } is the median-effect dose and m is the Hill coefficient. We used values recently reported from an in vitro study [46] and multiplied the IC_{ 50 } by the partition coefficient of TDF across the cell membrane in the pharmacokinetic model [42] to reflect the intracellular activity of TDF-DP.
Results
Heterogeneities in infection and immunity
The early transient (acute phase) and steady-state (chronic phase) trajectories of viral load are measurable and vary among individuals, potentially reflecting heterogeneity in underlying viral replication parameter values. To understand the insights that may be gained by including these heterogeneities in our model of adherence/resistance (Figure 1a), we began by identifying the “levers” within the virus dynamics model (Figure 1a and b) that can be used to reproduce these clinically observable heterogeneities. We hypothesized that the immune pressure parameter p, alone and in combination with other parameters, could be leveraged to effect these heterogeneities.
Observations of viremia during early HIV infection have revealed an early peak in viral load that is 2 logs (±1 log) higher than the setpoint viral load [48–52]. The time of the acute phase peak has been reported in the range of 12–31 days [53] or 5–19 days [54]. Some models lacking an immune compartment attribute the peak [55] and setpoint [56] viremia exclusively to the depletion of uninfected target cells rather than immune pressure. While such models could still achieve a 2-log reduction in viremia from peak to setpoint, the depletion of uninfected target cells is extreme when this drop is driven by target cell depletion alone, with nearly all CD4+ T cells ending up in the infected state. An immune response component was required to preserve a substantial fraction of uninfected target cells while producing a characteristic acute phase, both of which are features of realistic infection dynamics. Furthermore, correlations between CTLs and viremia in humans [57–59], as well as direct studies of CD8+ T cell depletion in primates [60–62], provide evidence for an important role for CTLs in bringing viral load down to setpoint. Therefore, we chose a model wherein both target cell depletion and CTL-mediated killing contribute to setpoint and acute phase virus dynamics. We then hypothesized that variability in immune response among hosts could, even in a simple model, reproduce the range of possible viral load trajectories among individuals.
A second, independent “lever” was required to modulate the acute phase viral load trajectory. This is because viral load during the acute phase is not predictive of setpoint viral load [53, 63] nor rate of AIDS progression [50, 63], although AIDS progression is correlated with setpoint viral load [64–66] and symptoms such as fever, vomiting, diarrhea, and headache during acute infection [50].
Frequent longitudinal sampling of patients with an estimated date of infection has revealed three equally common patterns of viral load trajectories: rapid decline to setpoint, a slow “shoulder” with gradual decline to setpoint, and an initial interval of viral suppression followed by a rebound of viremia, which takes over 90 days to return to setpoint [53]. With this range of biological variability in mind, we searched for a “lever” that can vary the shape (height, timing, and “shoulder”) of the acute phase viral load trajectory without changing setpoint viral load.
Sensitivity analysis of the virus dynamics model (Additional file 1: Table S1) found the acute phase peak time to be most sensitive to β and k, and the acute phase peak height to be most sensitive to k (closely followed by λ and u). Dividing k and multiplying β by a factor f preserves the number of virions produced with each replication cycle, but modulates the setpoint viral load. To keep the setpoint viral load constant, we used f to offset the change in setpoint viral load caused by p, changing the two parameters together as per Equation 9. We found that these counterbalanced parameters modulate the peak time, magnitude, and settling time of the peak of acute phase viremia while maintaining a constant setpoint viral load (Figure 2c-d). A biological interpretation of f might be ascribed to variability in non-CTL components of the host immune response, such as the relative magnitude of the neutralizing antibody response, which reduces β, and innate immune responses such as antiviral factors, which reduce k.
Both “levers” described here rely on the immune pressure parameter p, an expected source of heterogeneity among hosts. We next asked how p influences the response to antiretroviral therapy and the risk of resistance, using the context of misused tenofovir-based PrEP as a relevant scenario.
Resistance and reversion of mutations: effect of immune response
Human studies of TDF monotherapy have not agreed on a characteristic timescale for the emergence of resistance. Two human studies of daily TDF, lasting 21 days [67] and 28 days [68] respectively, did not find evidence of resistance. In the rhesus macaque SHIV model, TDF monotherapy elicited K65R mutations with associated viral rebound in all animals within 2 to 12 weeks (median 4 weeks) of therapy initiation [69]. This implies that the human trials were likely too short to detect the emergence of resistance. We therefore used our model to simulate longer timescales of drug exposure.
To investigate the effect of immune pressure on resistance and reversion rates, we next analyzed the effect of varying p, balanced with f (per Equation 9 and similarly to Figure 2c), on resistance and reversion times for a fixed interval of TDF monotherapy. We defined resistance as the first time K65R viral load reaches 50 copies/mL, and reversion as the crossover point between K65R and WT viremia. We analyzed a range of p that permits growth of the K65R mutant according to its R_{ 0 } at peak and trough TDF-DP concentrations, shown in Figure 3c. Note that the pharmacokinetic model with intracellular TDF-DP ensures that the peak-to-trough change in active drug concentration only brings R_{ 0 } of the WT above 1 for a very narrow range of p, illustrating the pharmacological “forgiveness” of daily TDF for most levels of immune pressure.
Increasing p lengthened the delay time until resistance from 29 days to ever-increasing values as R_{ 0 } of K65R declined toward 1, as shown in Figure 3d. The time from monotherapy cessation until reversion to a WT-dominated infection declined with increasing p over the same range. Reversion required a maximum of 37 days in the absence of CTL activity with p set to zero. As increasing p drove the R_{ 0 } of K65R toward 1, this duration dropped toward a minimum of 15 days (Figure 3e). Thus, increasing immune pressure led to slower appearance and more rapid reversion of the drug-resistant mutant. Intuitively, this result can be construed as a balance between immune pressure, which favors the WT due to the mutant’s fitness cost, and pharmacological pressure, which favors the mutant. In the drug-free condition shown in Figure 2b, a sufficiently high value of p drives the mutant R_{ 0 } below 1, only allowing the WT to survive in the absence of drug. With daily TDF dosing, it is the mutant that can survive for an intermediate range of p, and this range occurs at overall lower values of p compared to the drug-free state as a result of pharmacological pressure exerted on both quasispecies.
The dependence of resistance/reversion rates on immune strength may have further implications on the effect of late entry into therapy, or, in the case of PrEP, whether the accidental misuse of PrEP occurred during initial or early infection (due to breakthrough infection or false negative HIV tests during the “window period” of early infection) or well into infection (due to unauthorized drug use or failure of HIV testing).
Long-term effect of monotherapy: latent reservoirs
The net effect of increasing p in the short term is a shorter interval of resistance-dominated infection, i.e., fewer viral replication cycles involving the drug-resistant mutant. However, the model as used thus far allows the system to rapidly return to equilibrium after cessation of monotherapy, yielding no long-term implications of transient drug resistance. To investigate these long-term implications, we included another key feature of HIV infection: the incorporation of virus into long-lived latently infected cells capable of re-activating after months or years of suppressed viral replication.
For each mutant, we included a latently infected cell compartment w_{ i }, as shown in Equations 6–7 and illustrated in Figure 1c. We hypothesized that slow re-activation rates would endow this compartment with a “memory” of the virus subtypes that have undergone prior replication cycles, such that more intense and prolonged infection with drug-resistant virus would increase the fraction of drug-resistant mutants available for future re-activation. During non-suppressed viremia, the contribution of this compartment would be negligible compared to rates of forward- and back-mutation. However, during fully suppressive therapy, the latent compartment would become the sole source of new virions, and therefore would govern the production rate of drug-resistant mutant virions.
The level of drug-resistant mutants present in this compartment (w_{ mutant }) grows during monotherapy when replication of the mutant is favored over that of the wild-type, and declines in the absence of monotherapy while being gradually replaced with wild-type latently infected cells (w_{ WT }) until it reaches its equilibrium concentration. At equilibrium, both the latent and the active infected cell populations are primarily maintained through forward mutation of the wild-type virus due to its error rate q. When only 30 days of monotherapy are used to initially grow w_{ mutant }, this equilibrium level is reached within a decade, as shown in Figure 4b. When monotherapy is applied for longer intervals, such as 90 days (Figure 4c) or 1 year (Figure 4d), w_{ mutant } requires over a decade to reach equilibrium. Thus, as expected and shown in Figure 4e (graphed at the values achieved after 2 years of suppressive therapy), the fraction of the latent reservoir harboring drug resistance during the final period of suppressive therapy increases with increasing duration of initial monotherapy, except for short monotherapy exposures (<50 days) for which equilibrium is reached during the intervening decade. This threshold time would be shorter for shorter intervening spans between monotherapy and suppressive therapy. The fraction of drug-resistant free virus produced from the latent reservoir exhibits the same trend as the fraction contained in the latent reservoir, but is lower due to the reduced replicative capacity of the mutant, as shown in Figure 4f.
The previously discussed examples were conducted with immune pressure p, set to 1.5 × 10^{-5}. As a final step in analyzing the behavior of our extended model with latency, we varied p, which reduced the fraction of drug-resistant mutant in the latent compartment (Figure 4g) and in the free virus produced by the latent compartment (Figure 4h). As with the monotherapy duration, the trends were the same for latently infected cells and free virus produced, but the fraction of mutant virus was lower in the free virus pool due to the reduced replicative capacity of the mutant. At a sufficiently large p, no further reduction in the fraction of mutant was observed. Similar to the minimum duration of monotherapy required for variation in Figure 4e-f, this threshold occurs because back-mutation of the wild-type virus became the primary source of drug-resistant mutants.
We hypothesize that incorporation of drug-resistant virus into the long-lived latent reservoir could explain why past exposure to subclinical drug concentrations (e.g., single-dose nevirapine administration for prevention of mother-to-child-transmission) can lead to increased risk of treatment failure long after the initial drug exposure has waned and any rapidly-replicating viral populations should have returned to equilibrium. Thus, minimizing the total exposure to drug-resistant viremia would reduce the incorporation of these mutants into the latent pool.
Patterns of adherence: effect on drug resistance dynamics
We have already learned that increased immune pressure can reduce exposure to resistance by delaying its onset and hastening its reversion. However, the pattern of adherence to a drug regimen, including accidental PrEP during infection, is another important source of heterogeneities that may influence the total exposure the drug resistance and thus influence the long-term outcome. We therefore finish with an exploration of drug adherence patterns and their effect on total exposure to drug-resistant viremia.
The patterns of adherence – particularly, whether reversion can occur within the intervals of non-adherence – can therefore influence the overall exposure to drug resistance even for the same total fraction of doses missed. Knowing that immune pressure influences the required time for development and reversion of resistance, the critical dwell times in the dose-taking and dose-missing states would vary among individuals, such that populations with different immune phenotypes may exhibit different responses to particular patterns of drug adherence. Through our latency model, we further hypothesize that these changes in short-term resistance dynamics may influence treatment outcomes long after the exposure to drug resistance.
Discussion
Modern population models of HIV transmission are limited in their ability to predict the impact of drug resistance on individuals and the epidemic, because their representation of the development of drug resistance, if any, is simplistic. This simplification is chosen in part for lack of a way to model within-host dynamics in a manner that is consistent with observables such as viral load, which is heterogeneous across individuals, populations, and time periods [70]. We have shown how observable quantities such as setpoint and acute peak viral load can be independently tuned in such a model over the range of values observed in patients [53, 63]. Modulation of the parameters p and f in our model can recapitulate common patterns of early viremia: rapid decline to setpoint, a slow “shoulder” with gradual decline to setpoint, or an initial interval of viral suppression followed by a rebound of viremia, which approaches setpoint after 90 days or more [53].
The kinetics of tenofovir resistance due to K65R mutation are not fully characterized in humans [67–69]. By accounting for adherence patterns, pharmacokinetics, and the competition between quasispecies at time-varying drug concentrations, we were able to model the dependence of resistance and reversion rates on heterogeneities in host immunity and patterns of drug adherence. At constant overall levels of adherence, we showed the relationship between the duration of dose-taking and dose-skipping intervals and the resulting proportion of time that K65R dominated the infection. We also examined the extremes of periodic versus exponentially distributed transitions between drug-taking and drug-missing states. Reality is likely somewhere in-between: for example, lack of a prescription refill may cause monthly transitions into a drug-missing state, yet the transition back to a drug-taking state may be more variable, depending on the accessibility of the pharmacy, availability of drugs, transportation resources, and other factors.
Without a long-lived latent compartment, our model predicted rapid reversion of resistance with no longer-term impact on the host. With fully suppressive therapy, the populations of both WT and mutant virus drop toward extinction, as back-mutation of one quasispecies can no longer contribute to sustaining the others. In this state, re-activation of long-since-incorporated latent virus becomes a significant contributor to the viral pool. By including a latent compartment as has been described [32, 33], we were able to observe increasing proportions of K65R mutant in both the latent reservoir and the free virus population as the duration of TDF exposure increased, and as immune pressure p decreased.
This model feature provides a potential link between past monotherapy exposure (e.g., misuse of PrEP or single-dose nevirapine administration for prevention of mother-to-child-transmission) and future susceptibility to treatment failure. It explains why past exposure to monotherapy may increase the risk of treatment failure many years later. Further, it allows for host heterogeneities in immune response to be considered in evaluating the risk of this outcome.
Our analysis has thus far been limited to single-position mutations. In part, this is driven by the availability of detailed in vitro dose–response data for single-position mutants [46]. Additionally, the mutation rate of HIV is high enough that any given single-point mutant is likely to exist within a patient by the end of the acute phase. In contrast, fewer than one-third of double mutants are present during a WT-dominated infection, and a given triple-mutant is unlikely to already exist somewhere in the body [71]. In such situations, a stochastic simulation with discrete viral counts would be required to account for small populations such as rare mutants.
In the present analysis, we deliberately chose scenarios in which our continuous deterministic model would agree with the outcome of a stochastic model. We confirmed this through stochastic simulation of a subset of the experiment points. Although we found that simulations of a single milliliter did exhibit species fade-out, simulations of approximately 1% of the total system volume fully approached deterministic model results. Application of stochastic modeling to new areas in which the results could deviate from deterministic models, such as multi-position mutants, is an area of ongoing work.
Conclusions
To date, no model has linked the population-level dynamics of HIV transmission with a representation of within-host interactions of HIV virus dynamics, host immune response, and drug PK/PD [7]. We have shown that host immunity and patterns of drug adherence are important drivers of individual-level risks and dynamics of drug adherence. Understanding how within-host models can capture among-host heterogeneities is an important step toward the milestone of bridging within-host and population-level models of HIV drug resistance.
Declarations
Acknowledgements
We gratefully acknowledge Dr. C. Scott Wylie, Dr. Joshua L. Proctor, Dr. Daniel J. Klein, Dr. Erez Leiberman Aiden, and Dr. Aviva Presser Aiden for helpful discussions during the development of this model. The authors thank Bill and Melinda Gates for their active support of this work and their sponsorship through the Global Good Fund.
Authors’ Affiliations
References
- Taylor S, Boffito M, Khoo S, Smit E, Back D: Stopping antiretroviral therapy. AIDS. 2007, 21: 1673-1682. 10.1097/QAD.0b013e3281c61394.PubMedView ArticleGoogle Scholar
- Hurt CB, Eron JJ, Cohen MS: Pre-Exposure Prophylaxis and Antiretroviral Resistance: HIV Prevention at a Cost?. Clin Infect Dis. 2011, 53: 1265-1270. 10.1093/cid/cir684.PubMedPubMed CentralView ArticleGoogle Scholar
- Karim QA, Karim SSA, Frohlich JA, Grobler AC, Baxter C, Mansoor LE, Kharsany ABM, Sibeko S, Mlisana KP, Omar Z, Gengiah TN, Maarschalk S, Arulappan N, Mlotshwa M, Morris L, Taylor D, On behalf of the CAPRISA 004 Trial Group: Effectiveness and Safety of Tenofovir Gel, an Antiretroviral Microbicide, for the Prevention of HIV Infection in Women. Science. 2010, 329: 1168-1174. 10.1126/science.1193748.View ArticleGoogle Scholar
- Van Damme L, Corneli A, Ahmed K, Agot K, Lombaard J, Kapiga S, Malahleha M, Owino F, Manongi R, Onyango J, Temu L, Monedi MC, Mak’Oketch P, Makanda M, Reblin I, Makatu SE, Saylor L, Kiernan H, Kirkendale S, Wong C, Grant R, Kashuba A, Nanda K, Mandala J, Fransen K, Deese J, Crucitti T, Mastro TD, Taylor D: Preexposure Prophylaxis for HIV Infection among African Women. N Eng J Med. 2012, 367: 411-422. 10.1056/NEJMoa1202614.View ArticleGoogle Scholar
- Thigpen MC, Kebaabetswe PM, Paxton LA, Smith DK, Rose CE, Segolodi TM, Henderson FL, Pathak SR, Soud FA, Chillag KL, Mutanhaurwa R, Chirwa LI, Kasonde M, Abebe D, Buliva E, Gvetadze RJ, Johnson S, Sukalac T, Thomas VT, Hart C, Johnson JA, Malotte CK, Hendrix CW, Brooks JT: Antiretroviral Preexposure Prophylaxis for Heterosexual HIV Transmission in Botswana. N Eng J Med. 2012, 367: 423-434. 10.1056/NEJMoa1110711.View ArticleGoogle Scholar
- Baeten JM, Donnell D, Ndase P, Mugo NR, Campbell JD, Wangisi J, Tappero JW, Bukusi EA, Cohen CR, Katabira E, Ronald A, Tumwesigye E, Were E, Fife KH, Kiarie J, Farquhar C, John-Stewart G, Kakia A, Odoyo J, Mucunguzi A, Nakku-Joloba E, Twesigye R, Ngure K, Apaka C, Tamooh H, Gabona F, Mujugira A, Panteleeff D, Thomas KK, Kidoguchi L: Antiretroviral Prophylaxis for HIV Prevention in Heterosexual Men and Women. N Eng J Med. 2012, 367: 399-410. 10.1056/NEJMoa1108524.View ArticleGoogle Scholar
- Baggaley RF, Powers KA, Boily MC: What do mathematical models tell us about the emergence and spread of drug-resistant HIV?. Curr Opin HIV AIDS. 2011, 6: 131-10.1097/COH.0b013e328343ad03.PubMedPubMed CentralView ArticleGoogle Scholar
- Abbas UL, Anderson RM, Mellors JW: Potential Impact of Antiretroviral Chemoprophylaxis on HIV-1 Transmission in Resource-Limited Settings. PLoS One. 2007, 2: e875-10.1371/journal.pone.0000875.PubMedPubMed CentralView ArticleGoogle Scholar
- Abbas UL, Hood G, Wetzel AW, Mellors JW: Factors Influencing the Emergence and Spread of HIV Drug Resistance Arising from Rollout of Antiretroviral Pre-Exposure Prophylaxis (PrEP). PLoS One. 2011, 6: e18165-10.1371/journal.pone.0018165.PubMedPubMed CentralView ArticleGoogle Scholar
- Supervie V, Barrett M, Kahn JS, Musuka G, Moeti TL, Busang L, Blower S: Modeling dynamic interactions between pre-exposure prophylaxis interventions & treatment programs: predicting HIV transmission & resistance. Sci Rep. 2011, 1: 11-View ArticleGoogle Scholar
- Supervie V, García-Lerma JG, Heneine W, Blower S: HIV, transmitted drug resistance, and the paradox of preexposure prophylaxis. Proc Natl Acad Sci. 2010, 107: 12381-12386. 10.1073/pnas.1006061107.PubMedPubMed CentralView ArticleGoogle Scholar
- Hoare A, Kerr SJ, Ruxrungtham K, Ananworanich J, Law MG, Cooper DA, Phanuphak P, Wilson DP: Hidden Drug Resistant HIV to Emerge in the Era of Universal Treatment Access in Southeast Asia. PLoS One. 2010, 5: e10981-10.1371/journal.pone.0010981.PubMedPubMed CentralView ArticleGoogle Scholar
- Marks AJ, Pillay D, McLean AR: The effect of intrinsic stochasticity on transmitted HIV drug resistance patterns. J Theor Biol. 2010, 262: 1-13. 10.1016/j.jtbi.2009.09.017.PubMedView ArticleGoogle Scholar
- Nowak MA, Bangham CRM: Population Dynamics of Immune Responses to Persistent Viruses. Science. 1996, 272: 74-79. 10.1126/science.272.5258.74.PubMedView ArticleGoogle Scholar
- Wang G, Deem MW: Physical Theory of the Competition that Allows HIV to Escape from the Immune System. Phys Rev Lett. 2006, 97: 188106-PubMedView ArticleGoogle Scholar
- Ciupe SM, De Leenheer P, Kepler TB: Paradoxical suppression of poly-specific broadly neutralizing antibodies in the presence of strain-specific neutralizing antibodies following HIV infection. J Theor Biol. 2011, 277: 55-66. 10.1016/j.jtbi.2011.01.050.PubMedView ArticleGoogle Scholar
- Bangsberg DR, Moss AR, Deeks SG: Paradoxes of adherence and drug resistance to HIV antiretroviral therapy. J Antimicrob Chemother. 2004, 53: 696-699. 10.1093/jac/dkh162.PubMedView ArticleGoogle Scholar
- Rosenbloom DIS, Hill AL, Rabi SA, Siliciano RF, Nowak MA: Antiretroviral dynamics determines HIV evolution and predicts therapy outcome. Nat Med. 2012, 18: 1378-1385. 10.1038/nm.2892.PubMedPubMed CentralView ArticleGoogle Scholar
- Wu H, Huang Y, Acosta EP, Rosenkranz SL, Kuritzkes DR, Eron JJ, Perelson AS, Gerber JG: Modeling Long-Term HIV Dynamics and Antiretroviral Response: Effects of Drug Potency, Pharmacokinetics, Adherence, and Drug Resistance. J Acquir Immune Defic Syndr. 2005, 39: 272-283. 10.1097/01.qai.0000165907.04710.da.PubMedView ArticleGoogle Scholar
- Huang Y, Wu H, Holden-Wiltse J, Acosta EP: A dynamic Bayesian nonlinear mixed-effects model of HIV response incorporating medication adherence, drug resistance and covariates. Ann Appl Stat. 2011, 5: 551-577. 10.1214/10-AOAS376.PubMedPubMed CentralView ArticleGoogle Scholar
- Wahl LM, Nowak MA: Adherence and drug resistance: predictions for therapy outcome. Proc R Soc Lond B Biol Sci. 2000, 267: 835-843. 10.1098/rspb.2000.1079.View ArticleGoogle Scholar
- Cramer JA: How often is medication taken as prescribed? A novel assessment technique. JAMA. 1989, 261: 3273-3277. 10.1001/jama.1989.03420220087032.PubMedView ArticleGoogle Scholar
- Kenna LA, Labbé L, Barrett JS, Pfister M: Modeling and simulation of adherence: Approaches and applications in therapeutics. AAPS J. 2005, 7: E390-E407. 10.1208/aapsj070240.PubMedPubMed CentralView ArticleGoogle Scholar
- Arnsten JH, Demas PA, Farzadegan H, Grant RW, Gourevitch MN, Chang C-J, Buono D, Eckholdt H, Howard AA, Schoenbaum EE: Antiretroviral Therapy Adherence and Viral Suppression in HIV-Infected Drug Users: Comparison of Self-Report and Electronic Monitoring. Clin Infect Dis. 2001, 33: 1417-1423. 10.1086/323201.PubMedPubMed CentralView ArticleGoogle Scholar
- Liu H, Golin CE, Miller LG, Hays RD, Beck CK, Sanandaji S, Christian J, Maldonado T, Duran D, Kaplan AH, Wenger NS: A Comparison Study of Multiple Measures of Adherence to HIV Protease Inhibitors. Ann Intern Med. 2001, 134: 968-977. 10.7326/0003-4819-134-10-200105150-00011.PubMedView ArticleGoogle Scholar
- Sedaghat AR, Wilke CO: Kinetics of the viral cycle influence pharmacodynamics of antiretroviral therapy. Biol Direct. 2011, 6: 42-10.1186/1745-6150-6-42.PubMedPubMed CentralView ArticleGoogle Scholar
- Yi S, Nelson PW, Ulsoy AG: ASME Dynamic Systems and Control Conference (DSCC2008) October 20–22, 2008. 2008, Ann Arbor, Michigan, USA, 1549-1557.View ArticleGoogle Scholar
- Elliott E: Estimates of error rates for codes on burst-noise channels. Bell Syst Tech J. 1963, 42: 1977-1997.View ArticleGoogle Scholar
- Gilbert EN: Capacity of a burst-noise channel. Bell Syst Tech J. 1960, 39: 1253-1265.View ArticleGoogle Scholar
- Gillespie DT: Exact stochastic simulation of coupled chemical reactions. J Phys Chem. 1977, 81: 2340-2361. 10.1021/j100540a008.View ArticleGoogle Scholar
- Gillespie DT: Stochastic Simulation of Chemical Kinetics. Annu Rev Phys Chem. 2007, 58: 35-55. 10.1146/annurev.physchem.58.032806.104637.PubMedView ArticleGoogle Scholar
- Nowak MA, May RMC: Virus dynamics: mathematical principles of immunology and virology. 2000, USA: Oxford University PressGoogle Scholar
- Perelson AS: Modelling viral and immune system dynamics. Nat Rev Immunol. 2002, 2: 28-36. 10.1038/nri700.PubMedView ArticleGoogle Scholar
- Wodarz D, Nowak MA: Correlates of cytotoxic T-lymphocyte-mediated virus control: implications for immuno-suppressive infections and their treatment. Philos Trans R Soc B Biol Sci. 2000, 355: 1059-1070. 10.1098/rstb.2000.0643.View ArticleGoogle Scholar
- Grant MD, Smaill FM, Singal DP, Rosenthal KL: The influence of lymphocyte counts and disease progression on circulating and inducible anti-HIV-1 cytotoxic T-cell activity in HIV-1-infected subjects. AIDS. 1992, 6: 1085-10.1097/00002030-199210000-00004.PubMedView ArticleGoogle Scholar
- Gea-Banacloche JC, Migueles SA, Martino L, Shupert WL, McNeil AC, Sabbaghian MS, Ehler L, Prussin C, Stevens R, Lambert L, Altman J, Hallahan CW, De Quiros JCLB, Connors M: Maintenance of Large Numbers of Virus-Specific CD8+ T Cells in HIV-Infected Progressors and Long-Term Nonprogressors. J Immunol. 2000, 165: 1082-1092.PubMedView ArticleGoogle Scholar
- Ogg GS, Jin X, Bonhoeffer S, Dunbar PR, Nowak MA, Monard S, Segal JP, Cao Y, Rowland-Jones SL, Cerundolo V, Hurley A, Markowitz M, Ho DD, Nixon DF, McMichael AJ: Quantitation of HIV-1-Specific Cytotoxic T Lymphocytes and Plasma Load of Viral RNA. Science. 1998, 279: 2103-2106. 10.1126/science.279.5359.2103.PubMedView ArticleGoogle Scholar
- The International HIV Controllers Study: The Major Genetic Determinants of HIV-1 Control Affect HLA Class I Peptide Presentation. Science. 2010, 330: 1551-1557.PubMed CentralView ArticleGoogle Scholar
- Fideli ÜS, Allen SA, Musonda R, Trask S, Hahn BH, Weiss H, Mulenga J, Kasolo F, Vermund SH, Aldrovandi GM: Virologic and immunologic determinants of heterosexual transmission of human immunodeficiency virus type 1 in Africa. AIDS Res Hum Retroviruses. 2001, 17: 901-910. 10.1089/088922201750290023.PubMedPubMed CentralView ArticleGoogle Scholar
- Chun TW, Carruth L, Finzi D, Shen X, DiGiuseppe JA, Taylor H, Hermankova M, Chadwick K, Margolick J, Quinn TC: Quantification of latent tissue reservoirs and total body viral load in HIV-1 infection. Nature. 1997, 387: 183-188. 10.1038/387183a0.PubMedView ArticleGoogle Scholar
- Rong L, Perelson AS: Modeling Latently Infected Cell Activation: Viral and Latent Reservoir Persistence, and Viral Blips in HIV-infected Patients on Potent Therapy. PLoS Comput Biol. 2009, 5: e1000533-10.1371/journal.pcbi.1000533.PubMedPubMed CentralView ArticleGoogle Scholar
- Dixit NM, Perelson AS: Complex patterns of viral load decay under antiretroviral therapy: influence of pharmacokinetics and intracellular delay. J Theor Biol. 2004, 226: 95-109. 10.1016/j.jtbi.2003.09.002.PubMedView ArticleGoogle Scholar
- Anderson PL, Kiser JJ, Gardner EM, Rower JE, Meditz A, Grant RM: Pharmacological considerations for tenofovir and emtricitabine to prevent HIV infection. J Antimicrob Chemother. 2011, 66: 240-250. 10.1093/jac/dkq447.PubMedPubMed CentralView ArticleGoogle Scholar
- Hawkins T, Veikley W, St Claire RL, Guyer B, Clark N, Kearney BP: Intracellular pharmacokinetics of tenofovir diphosphate, carbovir triphosphate, and lamivudine triphosphate in patients receiving triple-nucleoside regimens. J Acquir Immune Defic Syndr. 2005, 39: 406-411. 10.1097/01.qai.0000167155.44980.e8.PubMedView ArticleGoogle Scholar
- Jullien V, Tréluyer J-M, Rey E, Jaffray P, Krivine A, Moachon L, Lillo-Le Louet A, Lescoat A, Dupin N, Salmon D, Pons G, Urien S: Population Pharmacokinetics of Tenofovir in Human Immunodeficiency Virus-Infected Patients Taking Highly Active Antiretroviral Therapy. Antimicrob Agents Chemother. 2005, 49: 3361-3366. 10.1128/AAC.49.8.3361-3366.2005.PubMedPubMed CentralView ArticleGoogle Scholar
- Sampah MES, Shen L, Jilek BL, Siliciano RF: Dose–response curve slope is a missing dimension in the analysis of HIV-1 drug resistance. Proc Natl Acad Sci. 2011Google Scholar
- Chou T-C: Relationships between inhibition Constants and Fractional Inhibition in Enzyme-Catalyzed Reactions with Different Numbers of Reactants, Different Reaction Mechanisms, and Different Types and Mechanisms of Inhibition. Curr Mol Pharmacol. 1974, 10: 235-247.Google Scholar
- Abu-Raddad LJ, Longini IA: No HIV stage is dominant in driving the HIV epidemic in sub-Saharan Africa. AIDS. 2008, 22: 1055-1061. 10.1097/QAD.0b013e3282f8af84.PubMedView ArticleGoogle Scholar
- Pilcher C, Tien H, Eron J, Vernazza P, Leu S, Stewart P, Goh L, Cohen M: Brief but efficient: acute HIV infection and the sexual transmission of HIV. J Infect Dis. 2004, 189: 1785-1792. 10.1086/386333.PubMedView ArticleGoogle Scholar
- Lavreys L, Baeten JM, Chohan V, McClelland RS, Hassan WM, Richardson BA, Mandaliya K, Ndinya-Achola JO, Overbaugh J: Higher Set Point Plasma Viral Load and More-Severe Acute HIV Type 1 (HIV-1) Illness Predict Mortality among High-Risk HIV-1–Infected African Women. Clin Infect Dis. 2006, 42: 1333-1339. 10.1086/503258.PubMedView ArticleGoogle Scholar
- Lavreys L, Baeten JM, Overbaugh J, Panteleeff DD, Chohan BH, Richardson BA, Mandaliya K, Ndinya-Achola JO, Kreiss JK: Virus Load during Primary Human Immunodeficiency Virus (HIV) Type 1 Infection Is Related to the Severity of Acute HIV Illness in Kenyan Women. Clin Infect Dis. 2002, 35: 77-81. 10.1086/340862.PubMedView ArticleGoogle Scholar
- Kumwenda JJ, Makanani B, Taulo F, Nkhoma C, Kafulafula G, Li Q, Kumwenda N, Taha TE: Natural History and Risk Factors Associated with Early and Established HIV Type 1 Infection among Reproductive-Age Women in Malawi. Clin Infect Dis. 2008, 46: 1913-1920. 10.1086/588478.PubMedView ArticleGoogle Scholar
- Kaufmann GR, Cunningham P, Kelleher AD, Zaunders J, Carr A, Vizzard J, Law M, Cooper DA: Patterns of viral dynamics during primary human immunodeficiency virus type 1 infection. J Infect Dis. 1998, 178: 1812-10.1086/314480.PubMedView ArticleGoogle Scholar
- Little SJ, McLean AR, Spina CA, Richman DD, Havlir DV: Viral Dynamics of Acute HIV-1 Infection. J Exp Med. 1999, 190: 841-850. 10.1084/jem.190.6.841.PubMedPubMed CentralView ArticleGoogle Scholar
- Stafford MA, Corey L, Cao Y, Daar ES, Ho DD, Perelson AS: Modeling Plasma Virus Concentration during Primary HIV Infection. J Theor Biol. 2000, 203: 285-301. 10.1006/jtbi.2000.1076.PubMedView ArticleGoogle Scholar
- De Boer RJ, Perelson AS: Target cell limited and immune control models of HIV infection: a comparison. J Theor Biol. 1998, 190: 201-214. 10.1006/jtbi.1997.0548.PubMedView ArticleGoogle Scholar
- Streeck H, Jolin JS, Qi Y, Yassine-Diab B, Johnson RC, Kwon DS, Addo MM, Brumme C, Routy J-P, Little S, Jessen HK, Kelleher AD, Hecht FM, Sekaly R-P, Rosenberg ES, Walker BD, Carrington M, Altfeld M: Human Immunodeficiency Virus Type 1-Specific CD8+ T-Cell Responses during Primary Infection Are Major Determinants of the Viral Set Point and Loss of CD4+ T Cells. J Virol. 2009, 83: 7641-7648. 10.1128/JVI.00182-09.PubMedPubMed CentralView ArticleGoogle Scholar
- Fellay J, Shianna KV, Ge D, Colombo S, Ledergerber B, Weale M, Zhang K, Gumbs C, Castagna A, Cossarizza A, Cozzi-Lepri A, De Luca A, Easterbrook P, Francioli P, Mallal S, Martinez-Picado J, Miro JM, Obel N, Smith JP, Wyniger J, Descombes P, Antonarakis SE, Letvin NL, McMichael AJ, Haynes BF, Telenti A, Goldstein DB: A Whole-Genome Association Study of Major Determinants for Host Control of HIV-1. Science. 2007, 317: 944-947. 10.1126/science.1143767.PubMedPubMed CentralView ArticleGoogle Scholar
- Altfeld M, Addo MM, Rosenberg ES, Hecht FM, Lee PK, Vogel M, Yu XG, Draenert R, Johnston MN, Strick D: Influence of HLA-B57 on clinical presentation and viral control during acute HIV-1 infection. AIDS. 2003, 17: 2581-10.1097/00002030-200312050-00005.PubMedView ArticleGoogle Scholar
- Sodora DL, Allan JS, Apetrei C, Brenchley JM, Douek DC, Else JG, Estes JD, Hahn BH, Hirsch VM, Kaur A, Kirchhoff F, Muller-Trutwin M, Pandrea I, Schmitz JE, Silvestri G: Toward an AIDS vaccine: lessons from natural simian immunodeficiency virus infections of African nonhuman primate hosts. Nat Med. 2009, 15: 861-865. 10.1038/nm.2013.PubMedPubMed CentralView ArticleGoogle Scholar
- Jin X, Bauer DE, Tuttleton SE, Lewin S, Gettie A, Blanchard J, Irwin CE, Safrit JT, Mittler J, Weinberger L, Kostrikis LG, Zhang L, Perelson AS, Ho DD: Dramatic Rise in Plasma Viremia after CD8+ T Cell Depletion in Simian Immunodeficiency Virus–infected Macaques. J Exp Med. 1999, 189: 991-998. 10.1084/jem.189.6.991.PubMedPubMed CentralView ArticleGoogle Scholar
- Madden LJ, Zandonatti MA, Flynn CT, Taffe MA, Cecilia M, Marcondes G, Schmitz JE, Reimann KA, Henriksen SJ, Fox HS: CD8+ cell depletion amplifies the acute retroviral syndrome. J Neurovirol. 2004, 10: 58-66.PubMedView ArticleGoogle Scholar
- Lindbäck S, Karlsson AC, Mittler J, Blaxhult A, Carlsson M, Briheim G, Sönnerborg A, Gaines H: Viral dynamics in primary HIV-1 infection. AIDS. 2000, 14: 2283-2291. 10.1097/00002030-200010200-00009.PubMedView ArticleGoogle Scholar
- De Wolf F, Spijkerman I, Schellekens PT, Langendam M, Kuiken C, Bakker M, Roos M, Coutinho R, Miedema F, Goudsmit J: AIDS prognosis based on HIV-1 RNA, CD4+ T-cell count and function: markers with reciprocal predictive value over time after seroconversion. AIDS. 1997, 11: 1799-10.1097/00002030-199715000-00003.PubMedView ArticleGoogle Scholar
- Mellors JW, Rinaldo CR, Gupta P, White RM, Todd JA, Kingsley LA: Prognosis in HIV-1 Infection Predicted by the Quantity of Virus in Plasma. Science. 1996, 272: 1167-1170. 10.1126/science.272.5265.1167.PubMedView ArticleGoogle Scholar
- O’Brien TR, Blattner WA, Waters D, Eyster ME, Hilgartner MW, Cohen AR, Luban N, Hatzakis A, Aledort LM, Rosenberg PS, Miley WJ, Kroner BL, Goedert JJ: Serum HIV-1 RNA Levels and Time to Development of AIDS in the Multicenter Hemophilia Cohort Study. JAMA. 1996, 276: 105-110. 10.1001/jama.1996.03540020027025.PubMedView ArticleGoogle Scholar
- Louie M, Hogan C, Hurley A, Simon V, Chung C, Padte N, Lamy P, Flaherty J, Coakley D, Di Mascio M: Determining the antiviral activity of tenofovir disoproxil fumarate in treatment-naive chronically HIV-1-infected individuals. AIDS. 2003, 17: 1151-10.1097/00002030-200305230-00006.PubMedView ArticleGoogle Scholar
- Barditch-Crovo P, Deeks SG, Collier A, Safrin S, Coakley DF, Miller M, Kearney BP, Coleman RL, Lamy PD, Kahn JO: Phase I/II trial of the pharmacokinetics, safety, and antiretroviral activity of tenofovir disoproxil fumarate in human immunodeficiency virus-infected adults. Antimicrob Agents Chemother. 2001, 45: 2733-10.1128/AAC.45.10.2733-2739.2001.PubMedPubMed CentralView ArticleGoogle Scholar
- Van Rompay KK, Johnson JA, Blackwood EJ, Singh RP, Lipscomb J, Matthews TB, Marthas ML, Pedersen NC, Bischofberger N, Heneine W, North TW: Sequential emergence and clinical implications of viral mutants with K70E and K65R mutation in reverse transcriptase during prolonged tenofovir monotherapy in rhesus macaques with chronic RT-SHIV infection. Retrovirology. 2007, 4: 25-10.1186/1742-4690-4-25.PubMedPubMed CentralView ArticleGoogle Scholar
- Fraser C, Hollingsworth TD, Chapman R, De Wolf F, Hanage WP: Variation in HIV-1 set-point viral load: Epidemiological analysis and an evolutionary hypothesis. Proc Natl Acad Sci. 2007, 104: 17441-17446. 10.1073/pnas.0708559104.PubMedPubMed CentralView ArticleGoogle Scholar
- Dorman KS, Zhou H: A branching process model of drug resistant HIV. Deterministic And Stochastic Models Of AIDS Epidemics And HIV Infections With Intervention. 2005, WORLD SCIENTIFIC, 457-496.Google Scholar
- Chun T-W, Carruth L, Finzi D, Shen X, DiGiuseppe JA, Taylor H, Hermankova M, Chadwick K, Margolick J, Quinn TC, Kuo Y-H, Brookmeyer R, Zeiger MA, Barditch-Crovo P, Siliciano RF: Quantification of latent tissue reservoirs and total body viral load in HIV-1 infection. Nature. 1997, 387: 183-188. 10.1038/387183a0.PubMedView ArticleGoogle Scholar
Copyright
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.