article: The kinetics of HIV-1 decay under treatment depends on the class of antiretrovirals used. Mathematical models are useful to interpret the different profiles, providing quantitative information about the kinetics of virus replication and the cell populations contributing to viral decay. We modeled proviral integration in short- and long-lived infected cells to compare viral kinetics under treatment with and without the integrase inhibitor raltegravir (RAL). We fitted the model to data obtained from participants treated with RAL-containing regimes or with a four-drug regimen of protease and reverse transcriptase inhibitors. Our model explains the existence and quantifies the three phases of HIV-1 RNA decay in RAL-based regimens vs. the two phases observed in therapies without RAL. Our findings indicate that HIV-1 infection is mostly sustained by short-lived infected cells with fast integration and a short viral production period, and by long-lived infected cells with slow integration but an equally short viral production period. We propose that these cells represent activated and resting infected CD4+ T-cells, respectively, and estimate that infection of resting cells represent ~4% of productive reverse transcription events in chronic infection. RAL reveals the kinetics of proviral integration, showing that in short-lived cells the pre-integration population has a half-life of ~7 hours, whereas in long-lived cells this half-life is ~6 weeks. We also show that the efficacy of RAL can be estimated by the difference in viral load at the start of the second phase in protocols with and without RAL. Overall, we provide a mechanistic model of viral infection that parsimoniously explains the kinetics of viral load decline under multiple classes of antiretrovirals.
  Viral dynamics analysis is a powerful tool to probe the lifecycle of viral infections [1–3]. In the case of HIV-1, treatment with reverse transcriptase inhibitors (RTI) and protease inhibitors (PI) resulted in a two-phase decline in plasma viral load [4–6]. Mathematical modeling attributed the first phase to the loss of productively infected cells with a short half-life (t1/2~0. 7 days), presumably activated CD4+ T-cells [7]. The second phase was attributed to the loss of “long-lived” infected cells with a slower loss rate (t1/2 ~14 days) [4]. The nature of the cells responsible for this second phase is still debated, and possibilities include, among others, macrophages or resting CD4+ T cells [4,8, 9]. In these early studies, the first phase lasted for about 6–11 days, and the second phase started after a viral load drop of 93%-99% from baseline, implying that the long-lived infected cells contributed ~1–7% of the total virus before treatment [4]. Different classes of drugs may lead to different patterns of viral decline, depending on where in the viral lifecycle the drug acts [10,11]. Modeling of one of the first clinical trials of HIV-1 treatment with an integrase strand transfer inhibitor (InSTI) showed that the viral load decline had first and second phase rates similar to those seen with InSTI-free regimens, with half-lives of ~1. 2 days and ~15. 5 days, respectively [12]. Although the data was sparse for the second phase, a unique characteristic of the viral load decline was that this phase started at much lower viral load levels than previously seen with InSTI-free regimens [12]. A recent study, with more frequent data, comparing InSTI-free and InSTI-containing regimens found that the rates of these two phases of decay were somewhat slower in the presence of an InSTI [13]. This study also confirmed that for InSTI-based regimens the viral load at the start of the second phase is ~1 log below the viral load at the start of the second phase for InSTI-free regimens [13] (Fig 1A). Another characteristic of viral decay under an InSTI is the separation of the first phase of decline into two sub-phases (phase 1a and phase 1b, Fig 1B) [10,14,15]. Analyzing the viral decay during the first 10 days of therapy with InSTI-containing regimens, we estimated [14] the half-life of the cells contributing to phase 1a as ~0. 8 days, similar to the previously estimated half-life of short-lived infected cells (~0. 7 days) [16]. In the same study, the half-life of the population of cells contributing to phase 1b was estimated as ~4. 6 days, which is much shorter than the previously estimated half-life of the cells contributing to the second phase (~14 days) [4]. Several mechanistic models were proposed to explain the differences in the profile of viral decay between InSTI-free and InSTI-containing regimens [9,10,12,14,15,17], but a thorough comparison with appropriate data is still lacking. Here, we analyze in detail several data sets under InSTI-based and InSTI-free therapies using a unified approach. Our experimental data includes more frequent sampling throughout the first month of the viral load decline under InSTI-based therapy than any previous modeling study. We develop mathematical models that separate the pre-integration and post-integration state in both short-lived and long-lived infected cells. We then use a mixed-effect modeling approach to simultaneously fit the viral load data in InSTI-free combination drug therapy and therapies including raltegravir (RAL), an InSTI [13,14,16,18]. We discovered that the best explanation for all the available kinetic data is that short-lived infected cells have a fast integration rate and a short viral production period, while long-lived infected cells have slow integration, but an equally short viral production period. Thus, we attribute the long life span of some infected cells to slow proviral integration. We propose that both types of cells (short- and long-lived) represent mostly infected CD4+ T-cells but in activated and resting states, respectively.
  We analyzed a total of 47 participants in three HIV-1 treatment data sets (see Methods) [13,16,18]. Eight participants were treated with a 4-drug combination without InSTI (“quad-therapy) [16]; eleven participants were treated with a combination of RAL and two RTIs (“RAL-combination”); and twenty-eight people were treated with RAL monotherapy for 9 days (“RAL-monotherapy”). Fig 1A presents the plasma viral load of all the participants from the quad-therapy and RAL-combination protocols, as well as the respective medians. This allows a clear visualization of the different patterns of decay under the two regimens (Fig 1B). In supplementary Figs B-D in S1 Text, we present the viral load data for all 47 individuals analyzed. Examination of these data shows that the second phase of viral load decline in individuals treated without RAL starts before than in participants under RAL-therapy (~ 5 days post-treatment initiation vs. ~10 days). At the start of the second phase, the viral load level in the RAL-combination participants is ~1-log lower than the viral load at the start of the second phase in the quad-therapy group (denoted by Δ in Fig 1A). We observe also that participants treated with RAL seem to have an intermediate phase of decline (phase 1b), which is slower than the first phase in the quad-therapy group, but faster than their second phase (Fig 1B). Finally, the slope of the last phase (phase 2) for RAL-combination therapy seems similar to the second phase slope for quad-therapy. Thus, these data sets with more frequent measurements confirm the differences in viral load decay profiles between InSTI-containing and InSTI-free regimens presented before [12–14]. These observations of the differences between the two types of drug regimens, lead to four questions: i) Are these empirical observations borne out by rigorous statistical analyses? ii) Could the early second phase in the quad-therapy be due to the same mechanisms (i. e. , decay of the same population of infected cells) as phase 1b in the InSTI-based treatment, which starts also at about ~5 days? iii) If not, what populations of infected cells contribute to each phase of decay in each type of treatment? And iv) Why do the second phases of decay seem have the same slope, but occur at very different viral load levels in the two types of regimens? We now answer these questions in turn. To put the observations of the previous section in a rigorous footing, we used a mixed-effects approach [19] using two empirical models, allowing for two or three exponential decays, respectively (see Methods). Two exponentials correspond to the two phases of decay (as it seems to be the case in the quad-therapy) and three exponentials correspond to three phases of decay (phase 1a, 1b and 2). We fitted these two models to the viral load of the three data sets simultaneously to assess which model is more consistent with the data. We found that the three-exponentials model did not fit the quad-therapy statistically better than the two-exponentials model (p = 0. 89), while for the RAL-containing regimens the reverse is true (p<0. 0001). Moreover, the initial decay rate was not significantly different between the treatment protocols, but the second decay rate was significantly different between the RAL-containing regimens (phase 1b) and the quad-therapy (phase 2) (p = 0. 0007), with estimates of ~0. 15 day-1 and ~0. 05 day-1, respectively (Table 1). The estimate for the third decay rate in the RAL-combination arm (phase 2) was ~0. 045 day-1 (Table 1), which is very similar to the phase 2 decay rate in the quad-therapy (~0. 05 day-1). These statistical results, thus, indicate that it is likely phase 1b for the RAL-containing treatments has a different viral source from phase two in the treatment without RAL. To explore the potential different populations of infected cells contributing to the viral load, we next analyzed the data in the context of mechanistic models of viral dynamics. We generalized a model developed to analyze the effects of RAL [14] by including the effect of reverse transcriptase and protease inhibitors and a compartment of long-lived infected cells, as illustrated in Fig 2A (details in Methods). To understand whether phase 2 (without RAL) and phase 1b (with RAL) have the same origin, we used the mechanistic model described by Eq (1) in Methods, without long-lived cells [10,14] (see section 2 in S1 Text). The model shows that there are two conditions that must be satisfied to observe two early phases of decay (viz. phase 1a and 1b): i) RAL must be present and ii) the efficacy of RAL must be above a critical threshold (ω>ξ (δ1+k) / (δ2+k), where ω is the efficacy of RAL, ξ is the combined efficacy of RTIs and PIs, δ1 and δ2 are the loss rates of short-lived infected cells before and after integration respectively, and k is the integration rate in these cells), but less than 100%. In the case of RAL monotherapy, i. e. , ξ = 0, the value of the critical threshold is 0, and thus with monotherapy one should always see two early phases of decay (section 2 in S1 Text) [15]. In the case of therapy without an InSTI, ω = 0 and the model predicts (see section 2. 2 in S1 Text) that after a short delay (so-called shoulder phase) only one exponential-like decay is observed early on, i. e. , during phase one [9]. The interpretation of these results is that in the presence of RAL, the loss of cells in the pre-integration stage is uncovered (phase 1b), whereas without RAL, integration proceeds so fast that it is not possible to observe the loss of those cells (i. e. , there is no phase 1b). Thus, to explain the second phase one needs another source of infected cells (e. g. long-lived HIV-infected cells). In the RAL-containing regimens, phase 1 (phases 1a and 1b together) lasts longer (~10 days) than phase 1 during quad-therapy (~5 days). Consequently, for RAL-containing regimens phase 2 starts at a later time and at a lower viral load (Fig 1 and [12,13]). But why does this happen? To investigate this issue, we used the model presented in Eq (1) (Fig 2A) and fitted it to the three datasets simultaneously. We found for the best fits that the rate of integration in the “long-lived cells” was very slow (k1<0. 05 day-1) when compared to that of “short-lived cells” (k = 2. 6 day-1). In addition, the estimate of the loss rate of long-lived cells producing virus (δM2 ~0. 90 day-1) was very close to the corresponding rate in short-lived cells (δ2 = 0. 86 day-1). This suggests the hypothesis that “long-lived cells” have slow integration leading to their long lifespan, but after becoming productively infected they are lost at a rate equal to that of short-lived cells, i. e. , δM2 = δ2. Indeed making this assumption gave the best fitting model (with lowest AIC, see Table C in S1 Text). Thus, in the best model to describe the data, we have one population of productively infected cells with a short lifespan, which is generated from two populations, one with fast integration events and another with slow integration events (Fig 2B). In what follows we present analyses of the data using this model, which we call the SRI model (Slow and Rapid Integration model) to differentiate from the standard model in Fig 2A. The SRI model has three compartments (Fig 2B), corresponding to cells with fast integration (I1) that quickly become productively infected (I2) and cells with slow integration (M1), which eventually also become productively infected (I2). In this new interpretation, the previous M2 compartment in the standard model is simply the subpopulation of productively infected cells (I2) that were generated from slow integration events (i. e. , from M1), but it does not correspond to a physiologically separate population of cells. Fig 3 shows the predicted population viral load profiles for each treatment protocol resulting from fitting the SRI model to all the data simultaneously, and Table 1 summarizes the corresponding parameters. In Fig 4, we show best fits of the SRI model to the viral load data for representative individuals from the three treatments protocols. Individual best fits for all participants are shown in Figs B-D in S1 Text and the estimated parameters in Table E in S1 Text. Overall, we found that the first phase of decline, which is phase 1 or 1a without or with RAL, respectively, corresponds to the loss of short-lived productively infected cells (~δ2, t½~ 20 h). Phase 1b only exists in RAL-based therapy and corresponds to the decay of short-lived cells pre-integration (at rate ~δ1+ (1-ω) k, t½~ 1. 8 days). Phase 2 corresponds to the decay of long-lived cells pre-integration (at rate ~δM1+ (1-ω) k1), which in the absence of RAL (ω = 0) is slightly faster (t½ ~ 19 days) than in the presence of RAL (t½ ~ 33 days), as RAL slows the loss of these cells by integration. Note that this estimated half-life in the absence of RAL is consistent with previous estimates (~14 days, range 6 to 24 days) [4,13]. Also, the difference in the second phase slope between the treatment regimens in the absence and presence of RAL (0. 037 day-1 vs. 0. 021 day-1) is sufficiently small that over 30 days of treatment it can’t be discerned in the figures (Fig 3). From the best-fit parameters, we can calculate the relative proportion of infection events leading to infected cells with fast and slow integration (see section 3. 1 in S1 Text). We estimate that ~96% of infections result in cells with fast integration (I1) and only 4% generate infected cells with slow integration (M1). Yet, due to the long-lived nature of the latter population, before treatment the pool of cells with slow integration (M1) represents about 40% of total infected cells in the SRI model (i. e. , M1 is 40% of I1+M1+I2). These results offer for the first time a consistent picture for the three treatment scenarios based on the same mechanistic model. However, we still need an explanation for the observed lower viral load during the second phase decline in InSTI-based therapies. Based on the SRI model, we can directly calculate the reduction in viral load during the second phase with RAL treatment with respect to the viral load in the second phase with quad therapy. We find that this reduction is approximately 1− (1−ω) eωk1t (section 3. 2 in S1 Text). For example, at 12. 5 days post-treatment, i. e. the beginning of phase two in the RAL-containing regimen, the viral load in the RAL-combination therapy is reduced by ~93% with respect to the viral load in the quad-therapy group (Fig 3). Interestingly, if we linearly extrapolate the second phase of decay for therapy without and with RAL back to time 0, then this reduction in viral load is given by ~ω (Fig 3). Therefore, our model provides a way to directly estimate the efficacy of RAL in blocking integration, ω, when compared to a highly-potent therapy without RAL, such as the quad-treatment. Thus, the SRI model intrinsically explains this difference in viral loads, based on the dynamics of different infected cell populations, and without the need for any extra assumptions. Still, it would be important to have some intuitive understanding why there is a lower level of virus at the beginning of phase 2 in RAL-based regimens. In the SRI model, productively infected cells (I2) generate the observed virus. The different phases of virus decay correspond to loss of different cell populations contributing to generating those productively infected cells. In more detail and as shown in Movie S1, viral decay is observed as follows. Initially productively infected cells (I2) are lost quickly (rate δ2), corresponding to phase 1 or phase 1a, depending on the type of treatment. At the same time, cells progressing through integration (I1 and M1) replenish the productively infected cells (I2). In the presence of RTIs and PIs without InSTIs, the pool of cells with fast integration (I1) is quickly exhausted. In this case, the second phase can be observed when the main contribution to productively infected cells (I2) comes from the conversion of slowly integrating cells (M1) at rate k1 into productively infected cells. The productively infected cells then have an effective decay rate given by (~δM1+k1), which is the rate limiting step. Paradoxically, in the presence of an InSTI, the pool of cells that can integrate fast (I1) is exhausted more slowly, because integration is blocked/slowed but not prevented completely. This results in phase 1b with slope (~δ1+k (1-ω) ) due to the slow decay of cells in I1. Together phases 1a and 1b last longer than phase 1, because it takes longer to lose most cells in I1. Again, as in the quad-therapy, the second phase is observed when the main contribution to productively infected cells comes from the slow integrating cells (M1), which in the presence of an InSTI occurs later, because integration is slowed in this compartment too. This explains the lower viral load level at the start of the second phase with an InSTI regimen. All the infected cell populations are decaying from early on at the start of treatment, but the observed viral load reflects the decline of different populations in turn, as the populations with faster turnover are lost. However, it raises another question, why does the second phase in the quad-treatment start so early (~5 days), even earlier than what had been observed before in RTI+PI treatments? In the quad-therapy data set the first phase lasted only ~5 days, whereas early studies with less potent RTIs and PIs indicated that this phase lasted for about 6 to 11 days [4]. Using the SRI model, we found that the higher the effectiveness of the RTI+PI regimen (i. e. , the larger the effectiveness of the RTI, η and of the PI, ε), the earlier the switch from the first to the second phase (Fig E (a) in S1 Text). In fact, we can calculate (Eq. S15) that the difference in the times of transition from ~8 days in the earlier studies to 5 days with the quad-regime implies that the latter had an efficacy ~1. 6-fold higher, i. e. , early therapies had an efficacy ξ ~ 0. 62. This result compares favorably with a previous study that indicated that early therapies had an efficacy ~75% of the quad-treatment protocol [20]. In the quad regimen with higher efficacy, the transition from the first to second phase occurs earlier and at a slightly lower viral load, because of faster viral decay (Fig E (b) in S1 Text), as discussed in [16]. Thus, the SRI model clearly explains the early transition to the second phase in the quad-therapy regimen is due to the high efficacy of the combination treatment.
  Analysis of viral dynamics under different treatment regimens provides insight into the possible virus sources for the observed phases of plasma viral load decay. Under RAL-containing regimens the original phase 1 decay described for PI and RTI containing regimes [4,5] is replaced by two phases, called 1a and 1b [10,14,15], which are then followed by a second phase that starts at lower viral load levels than the start of the second phase with InSTI-free regimens due to the presence of phase 1b. Here we presented for the first time a detailed quantification of the dynamics of cell subpopulations contributing to the observed viral kinetics over the first month of treatment. There are other slower phases of decay at later times, which are not studied here [3,21,22]. The mechanistic reasons for the peculiar viral load decay profile under InSTI regimens have been the subject of previous studies. Murray et al. [12] first developed models to analyze the viral load decay kinetics under RAL. The models presented in [12] are different from our SRI model, they were not fit to individual patient data, and they cannot recapitulate the phase 1b observed in the data. However, the models in [12] that could qualitatively explain the data had the interesting characteristic that the second phase virus was produced from the same cell population as the first phase virus. The difference was in the infection kinetics (e. g. , proviral integration), which could be fast (in short-lived cells) or slow (in long-lived cells). This idea echoes our final model where we found that the death rate of short-lived cells was the same as the death rate of long-lived cells after they both start producing virus (i. e. , δM2 = δ2). Sedaghat et al. [9,17,23], Gilmore et al. [10] and Wang et al. [15] presented detailed mathematical analyses of models similar to ours. However, they did not analyze clinical data and/or had only very sparse viral load data, and, thus, had to make various assumptions about parameters and their relationships. For example, Sedaghat et al. [9,17,23] studied without fitting data the features of the second phase of viral decline. They concluded that in the most likely scenario the loss rate of long-lived infected cells pre-integration (~δM1+k1) should be greater than the loss rate of long-lived productively infected cells (~δM2); and that the latter should be the slope of the second phase (equal) in both treatment cases (with and without integrase inhibitors). By quantitative analyses of viral load data and fitting, we rather found that δM2>δM1+k1 even when using the conclusions in Sedaghat et al. [9,17,23] as initial guesses in the fitting procedure. Thus, we now show that the second phase slope is the loss rate of long-lived cells pre-integration, i. e. , ~δM1+ (1-ω) k1, and that this slope is slightly different with and without an integrase inhibitor, as also found previously by simple linear regression [13]. Interestingly, Sedaghat et al. briefly considered an “intriguing possibility” consistent with our results that the loss rate of long-lived productively infected cells is the same as that for short-lived productively infected cells, but dismissed it [17]. Our detailed data-driven comparison between treatment regimens with and without RAL leads to a novel interpretation of the profile of HIV-1 decay under treatment. The existence of an early (fast) and a late (slow) phase of viral decay is explained, as before, by decay of two “different” cell populations. However, the best fit of the SRI model supports the idea that the dynamics of proviral integration distinguishes these two cell populations, but they both could be CD4+ T-cells. The early fast decay in viral load (including phase 1a and 1b) is most likely due to infection of activated CD4+ T-cells, which quickly progress to viral production. The slow late decay (phase 2) is likely due to infection of resting CD4+ T-cells, which then progress slowly to provirus integration. This interpretation is consistent with data indicating that in resting CD4+ T-cells integration proceeds slowly [24–26]. In the SRI model a fraction of these cells eventually complete integration and become productively infected cells. The resting cells have a slow loss rate before integration, thus they are long-lived, but we found that after becoming productively infected they have a high loss rate, similar to short-lived productively infected cells. One explanation for this result could be that some stimulus activates these cells, which then quickly integrate and start producing virus. Although this interpretation of the data may be surprising, with hindsight, one recognizes that the differences observed in viral load decay profiles between treatments with and without RAL must be due to some process involving proviral integration, which is the step of the lifecycle blocked by the drug. We note that this interpretation is unlikely to depend on the assumptions of the SRI model, which we tested by examining other models and parameter ranges without finding any important effect (section 4. 4 in S1 Text). In particular, alternative models, where the second phase is due to cells with virus already integrated that are lost slowly, does not generate viral decay profiles consistent with those observed. Indeed, we started our analysis with this assumption (Fig 2A) and found that the resulting differences in viral loads in the second phase between treatments with or without RAL were too small compared with the data (see section 4. 3 and Fig A in S1 Text). We estimate that before treatment, the vast majority of virus is produced by short-lived cells [4], which is indicated by the large decay in viral load during the first phase, and only a small fraction of virus is produced from cells with slow proviral integration. Several recent reports have indicated that non-integrated HIV DNA in resting CD4+ T-cell is capable of integration and generating productive infection upon activation [24,27–32]. Indeed, this phenomenon has been dubbed “pre-integration latency” [7]. Here we also find that the loss of cells in this pre-integration stage in the absence of InSTI therapy is very slow (t1/2 ~ 19 days), consistent with a study measuring the decay of unintegrated HIV DNA in chronically infected patients under successful HAART that reported a half-life of 26 days [33]. Some studies have described a very rapid decay of the pre-integration complex (PIC) with half-life of ~1 day in resting cells in vitro [24,28]. However, there are other studies that report stable accumulation of HIV DNA transcripts over longer timespans [34] with inducible pre-integration latency over times up to 28 days in vitro [27,30]. It has also been suggested that in vivo the half-life of pre-integration latent cells can be longer than in vitro [35]. Moreover, recent evidence has shown that 2-LTR circles, which are very stable, can be cleaved by the viral integrase forming a viable substrate for integration and rescuing viral production [29]. Altogether, it is possible that in vivo resting infected cells in the pre-integration stage have a broad distribution of half-lives. Thus, although the half-life of “pre-integration latency” is somewhat controversial, we think that our results add a new piece of evidence to indicate that in vivo on average this half-life is long and physiologically relevant, since these cells can be activated into producing virus. An alternative explanation invoked for the second phase of viral decay is that the long-lived cells are infected macrophages [4,9, 23]. Some results indicate that HIV DNA integration is slow in primary blood monocytes [36], however other studies have measured fast integration (~3. 4 hours) [23]. In any case, it is thought that infected macrophages are long-lived even after starting to produce virus [7], which is inconsistent with our findings. Moreover, it is unlikely that infected monocytes with unintegrated HIV DNA form as large a pool of infected cells as estimated here [37]. However, it is possible that in addition to the cell compartments in the SRI model, a compartment of long-lived productively infected cells, such as macrophages (Fig 2A), exists. The contribution of these cells to the second phase in viral load should be minimal and, thus, it would not affect the results presented by the SRI model. Altogether, the most parsimonious interpretation of the analyses of these detailed data sets is that the slow second phase viral decay corresponds to the loss of infected cells with unintegrated provirus (possibly resting CD4+ T-cells) by death, by loss of the pre-integration intermediates, or by integration of the HIV DNA to generate productively infected cells. This description is valid both in regimens with and without an InSTI, and thus alters the interpretation of the second phase of viral decay presented in previous studies [4,9, 17]. Furthermore, our results explain why the second phase decay is slightly slower in the RAL containing regimen than in the RAL-free regime, as found previously [13]. This difference is small because it depends on the integration rate k1, which is slow. In conclusion, we have proposed a new model that explains the main differences of viral load decline under InSTI therapy (as prototyped here by RAL) when compared with InSTI-free regimens: i) the two early phases 1a and 1b of viral decay; ii) the ~1 log lower viral load at the start of the second phase; and iii) the nearly identical second phase slope with and without an InSTI. Fitting this model to frequently sampled viral load data indicates that the second phase of viral decay is due to a subset of cells completing integration very slowly, which most likely are infected resting CD4+ T-cells that eventually complete integration to produce virus and then die quickly.
  The first set of data is from a trial with a highly active antiretroviral quad-regimen [16], where nine chronically HIV-1-infected individuals were treated with the boosted PI lopinavir-ritonavir (1,066 and 266 mg/day, respectively), and the RTIs efavirenz (600 mg/day), lamivudine (300 mg/day), and tenofovir DF (300 mg/day). Plasma viral loads were measured by RT-PCR (Roche Amplicor Ul- trasensitive Cobas 1. 5) with a limit of quantification (LoQ) of 50 copies/ml at 6 h intervals during the first 72 h, then daily until day 10, and then weekly until day 28 [16]. One participant could not be analyzed, because the precise time of each viral load quantification was not available. We refer to this data set as the “quad treatment” data. The second data set is from 28 HIV-1-infected participants treated with RAL monotherapy for 9 days [18]. Participants received 100,200,400, or 600 mg RAL twice daily, and had plasma HIV RNA measured before the first dose, 6 and 12 hours post-dose, and on days 1,2, 3,4, 7,9. Plasma was assayed for HIV-1 RNA using the Amplicor HIV-1 Monitor Assay, version 1. 5 with a LoQ = 400 copies/ml. If the result was below that level, the UltraSensitive HIV-1 Monitor Assay (LoQ = 50 copies/ml) was used (Roche Molecular Diagnostics, Alameda, CA). We treated plasma HIV RNA data below the limits of quantification as censored at the corresponding values. Because no difference in the pattern of decay or efficiency between the different doses was found [14,18], we analyze all doses together. We refer to this data set as the “RAL-monotherapy” data. The third set of data was obtained from 11 HIV-1-infected participants treated with a combination of FTC/TDF 200 mg/300 mg daily plus RAL 400 mg twice daily enrolled in the intensive viral dynamics A5249s substudy of ACTG protocol A5248 [13]. Plasma HIV-1 RNA was measured at baseline, and at 2,4, 6,12,18,24,30,36,42, and 48 hours and days 3,4, 7,10,14,21 and 28 after treatment initiation. Plasma was assayed for HIV-1 RNA using the Amplicor HIV-1 Monitor, version 1. 5, UltraSensitive protocol (LoQ = 50 copies/mL; Roche Molecular Systems, Branchburg, NJ). We refer to this data set as the “RAL-combination therapy” data. The study protocols were approved by the respective institutional review board at each of the participating clinical research sites, and all subjects provided signed informed consent (see details in [13,16,18]). All data analyzed were anonymized. We first analyzed the data with double and triple exponential decay curves, corresponding to sums of (j = 2 or j = 3) terms of the form Cje−αjt, to see if the data was consistent with two or three phases of decay and what the respective slopes (αj) were. We then generalized our model developed to analyze the effects of RAL [14] by including the effect of protease inhibitors and a compartment of long-lived infected cells, to account for the long-term (up to 30 days) follow-up. The model is a modification of the standard model of virus dynamics [1,3] separating the compartment of infected cells in pre-integration and post-integration states in both short-lived and long-lived infected cells. A schematic of the model is shown in Fig 2A, and the full system of differential equations describing the dynamics is given in the S1 Text (section 1). We assumed that each drug has the same efficacy in short and long-lived infected cells and that the dynamics of virus is much faster than that of infected cells. Both of these assumptions have been widely used before [4,5, 9,14,15,17], and result in the simplified model dI^1dt= (1−ξ) T^V−δ1I^1−k (1−ω) I^1dM^1dt= (1−ξ) M^V−δM1M^1−k1 (1−ω) M^1dVIdt=kc (1−ω) I^1−δ2VIdVMdt=k1c (1−ω) M^1−δM2VM (1) where VI and VM denote virus produced by short- and long-lived infected cells, I and M, respectively, total virus is V = VI + VM, ω is the effectiveness of the InSTI and ξ is the combined effectiveness of the PIs (ε) and RTIs (ε) given by (1-ξ) = (1-η) (1-η). See section 1 in S1 Text for further details about the definitions of the variables and parameters. In particular, note that without loss of generality, we rescaled the variables so that the rate of infection, which in our model includes reverse transcription, and the viral production rate no longer appear in the equations, although they can be different between short-lived, I, and long-lived, M, cells. We fitted the model to the plasma HIV-1 RNA data using non-linear mixed effects (NLME) models. In this approach, we represent the parameters for each individual (i) as μi = θeφi, where θ is the median value of the parameter in the population, and φi the random terms, which are normally distributed with zero mean and a variance to be estimated. We fitted the model in Eq (1) simultaneously to the three sets of data using the software MONOLIX (www. lixoft. eu), to estimate the population parameters (by maximum likelihood) and the variances of the random effects. In total, we fit 613 data points of 47 participants simultaneously to estimate between 4 and 6 parameters, and their corresponding variances in various versions of the model. For each model fit, we estimate the log-likelihood (log L) and compute the Akaike Information Criteria (AIC = -2log L+2m, where m is the number of parameters estimated) [38]. We compared the AIC for different model assumptions and computed ΔAIC, the difference in AIC between the best model and other model scenarios analyzed. We assumed models had similar support if their ΔAIC<2 [38]. P-values were based on the log-likelihood ratio test and significance was assessed at the α = 0. 05 level. To implement the fitting, we first explored the possible range of values for the new parameters in our model, δM1, δM2 and k1 describing long-lived infected cells (see Fig 2A), keeping fixed the remaining parameters based on estimates from previous studies [14,16,39], and obtaining the values of T^ and M^ from steady state assumptions (see Eq. S. 4). We performed Latin-hypercube sampling resulting in 64,000 (40 values for each parameter) simulations of the model in Eq (1) for different values of the parameters [40]. We looked for combinations of parameters that showed viral decay properties consistent with the observations in the data (section 4. 1 in S1 Text). Namely, a reduction in the viral load greater than 70% at the start of the second phase of decline under RAL-combination therapy relative to start of the second phase in the quad-based treatment; and similar second phase slopes. We found that these biological conditions are only satisfied for small values of δM1 (<0. 07 day-1) and k1 (<0. 08 day-1), along with values of δM2 greater than 0. 15 day-1 (Fig A in S1 Text). These results provided starting guesses for these parameters in the viral load fits. To simplify the fitting procedure and obtain convergence, we fixed five parameters (ξ, k, ω, c, and the fraction of virus produced from long-lived infected cells, see section 4. 2 in S1 Text) at values previously estimated by us and others [14,39]. In addition, the values of T^ and M^ in Eq (1) were obtained from the other parameters based on pre-therapy steady state assumptions (see Eq. S. 4 in S1 Text). Still, it was difficult to estimate independently both δM1 and k1, since when one increased the other tended to decrease (see Tables A and B in S1 Text). Therefore, we performed the fitting procedure for different fixed values of δM1, and found the best fits when δM1 = 0. 01–0. 02 day-1, based on AIC (see Table B in S1 Text). We then estimated the remaining parameters: V0, δ1, δ2, δM2 and k1.
summary: Antiretroviral therapy in HIV-1 infection leads to characteristic viral load decay profiles. Interpretation of this profile with mathematical models has provided important insights into the dynamics of the viral lifecycle (e. g. , the turnover of infected cells). Here we develop a unified model and analyze viral load data from HIV-infected participants under treatment with and without an integrase inhibitor. Our model offers a new explanation for the observed differences in the decay profiles, and strongly supports the hypothesis that the second phase of viral decay is due to the loss of long-lived resting CD4+ T-cells with slow proviral integration. We estimate parameters of this decay and estimate the efficacy of an integrase inhibitor.
keywords: blood cells
  medicine and health sciences
  radiochemistry
  immune cells
  pathology and laboratory medicine
  viral transmission and infection
  pathogens
  immunology
  microbiology
  enzymology
  mathematical models
  nuclear decay
  retroviruses
  viruses
  immunodeficiency viruses
  rna viruses
  enzyme inhibitors
  viral load
  research and analysis methods
  white blood cells
  animal cells
  medical microbiology
  hiv
  mathematical and statistical techniques
  microbial pathogens
  t cells
  chemistry
  hiv-1
  physics
  biochemistry
  macrophages
  cell biology
  nuclear physics
  virology
  viral pathogens
  biology and life sciences
  cellular types
  physical sciences
  protease inhibitors
  lentivirus
  organisms
title: Treatment with integrase inhibitor suggests a new interpretation of HIV RNA decay curves that reveals a subset of cells with slow integration