Skip Navigation

Virologic Response to Potent Antiretroviral Therapy and Modeling of HIV Dynamics in Early Pediatric Infection

  1. Paul Palumbo1,a,
  2. Hulin Wu2,a,
  3. Ellen Chadwick4,
  4. Ping Ruan2,
  5. Katherine Luzuriaga3,
  6. John Rodman5,b and
  7. Ram Yogev4

    for the Pediatric AIDS Clinical Trials Group 345 Investigatorsc c

  1. 1 Department of Pediatrics, University of Medicine and Dentistry-New Jersey Medical School, Newark, New Jersey
  2. 2 Statistical and Data Analysis Center, Harvard School of Public Health, Boston
  3. 3 Department of Pediatrics, Program in Molecular Medicine, University of Massachusetts
  4. 4 Children's Memorial Hospital, Northwestern University, Chicago, Illinois
  5. 5 St. Jude Children's Research Hospital, Memphis, Tennessee
  1. Reprints or correspondence: Dr. Paul Palumbo Section of Infectious Diseases and International Health Dartmouth-Hitchcock Medical Center and Dartmouth Medical School One Medical Center Drive Lebanon NH 03756 (paul.palumbo{at}dartmouth.edu).
  • Present affiliations: Section of Infectious Diseases and International Health, Dartmouth-Hitchcock Medical Center and Dartmouth Medical School, Lebanon, New Hampshire (P.P.); Department of Biostatistics and Computational Biology, University of Rochester School of Medicine and Dentistry, Rochester, New York (H.W.).

  • Deceased. We dedicate this article to the memory of our long-standing colleague and team member, John Rodman.

Abstract

Background. Human immunodeficiency virus (HIV) infection in infancy features a persistently high viral load and elevated antiretroviral drug clearance rates, which pose significant therapeutic challenges to the clinician. Viral and cellular kinetic analyses performed in HIV-infected adults have yielded significant insights into the dynamic setting of this viral infection. Similar studies are needed in pediatric populations, in whom differing dynamics might translate into age-specific treatment approaches.

Methods. Viral and cellular kinetic analyses were performed using a nonlinear mixed-effects model in a cohort of 48 infants 1–24 months of age enrolled in a trial of ritonavir-based highly active antiretroviral therapy (HAART).

Results. Infected cell compartment kinetics were comparable with reported adult values, with no age-specific differences demonstrated—suggesting the ability to suppress viral replication in infants receiving HAART. Comparisons between 2 ritonavir dosing schedules revealed significant improvement in phase 1/2 decay constants in favor of the higher dose. A negative correlation was established between plasma RNA levels and phase 1 decay rates, which has worrisome implications for infant therapeutics given high infant pretreatment plasma virus levels.

Conclusions. Ritonavir-based HAART regimens in infancy result in HIV decay constants comparable to those reported in adults, without age-specific variability. Despite higher plasma HIV levels and CD4 lymphocyte counts in infancy, HAART can result in timely, effective control of viral replication.

HIV infection in the first few years of life possesses some unique features and poses significant therapeutic challenges to the clinician. It has been well documented in natural history studies of HIV infection that, although peak plasma virus levels during acute infection among adults and infants are similar, subsequent declines in plasma virus levels require considerably longer time in infancy, resulting in higher systemic viral exposure during the first 2 years of life [1, 2]. These sustained high virus levels are temporally related to observed increased rates of disease progression in early infancy and are likely to be a causative factor. This phenomenon occurs in the context of multiple variables that may impact virus-cell homeostasis during infancy: proportionately higher numbers of CD4 lymphocytes relative to body size, a developing and relatively immature immune system, significant somatic cell proliferation, and pervasive immune system activation through natural and iatrogenic antigen exposure.

Viral and cellular kinetic analyses performed in HIV-infected adults have yielded significant insights into the dynamic setting of this viral infection [3, 4]. Implications of these findings include remarkably high levels of viral production and cell turnover (supporting the need for early, aggressive therapeutic intervention) and the risk of rapid emergence of viral resistance to antiretroviral agents [57]. Similar studies are needed in pediatric populations, especially among young infants whose unique features suggest the possibility of different cellular and viral dynamics in comparison with adults. In fact, preliminary data from a small number of infants receiving a potent, non- protease inhibitor-containing antiretroviral regimen suggest age-specific viral turnover rates [8].

The analyses reported here were performed on a cohort of 48 infants <2 years of age who were enrolled in Pediatric AIDS Clinical Trials Group (PACTG) 345, a phase 1/2 trial of the safety and antiviral effects of a potent 3-drug regimen (zidovudine [AZT], lamivudine [3TC], and ritonavir [RTV]) [9]. Key questions addressed here include the following: (1) Will HIV-infected infants experience “adultlike” responses to highly active therapy? (2) Are viral/cellular production and turnover rates similar when comparing infants with adults? (3) What are the relative viral and cellular kinetics of 2 dosing regimens?

Methods

Study population. Forty-eight HIV-infected infants between the ages of 1 and 24 months were enrolled and stratified into 3 age groups: 1–3 months, 3–6 months, and 6–24 months (table 1) [9]. All infants were naive to 3TC and protease inhibitor antiretroviral therapy. On enrollment, all infants were treated with the following antiretroviral regimen: AZT at 160 mg/m23 times a day; 3TC at 4 mg/m2 every 12 h; and RTV at 350 mg/m2 every 12 h in cohort 1 and at 450 mg/m2 every 12 h in cohort 2. At the inception of this trial, pediatric dosing for children >2 years of age was 350–400 mg/m2every 12 h. For cohort 1, RTV monotherapy was administered for 7 days, after which triple therapy with the addition of AZT and 3TC was instituted. Infants in cohort 2 received triple therapy from the outset. Infants were enrolled into cohort 1 between July 1997 and March 1998, and cohort 2 infants were enrolled between February 1998 and March 2001. This clinical trial was conducted with the approval of the institutional review boards at each participating clinical trial site, and informed consent was obtained from the parent or guardian of each infant.

Table 1.

Baseline characteristics.

Monitoring. Intensive blood sampling for determining plasma RNA levels was performed before entry; at baseline; at days 1, 3, 7, 14, and 28; monthly from 1 to 18 months; and every 2 months thereafter. The NucliSens assay (Organon Teknika) was used to determine quantitative plasma RNA levels, with a lower limit of quantitation of 400 RNA copies/mL [10]. All assays were performed in a single laboratory certified by the Division of AIDS (DAIDS), National Institute of Allergy and Infectious Diseases, and incorporated quantitative standards supplied by the Viral Quality Assurance Laboratory of DAIDS [11]. Plasma samples were stored at −70°C until the assay was performed. Samples from the first 28 days of the study were assayed in batch format to avoid interassay variability, whereas subsequent samples were assayed in real time. Virologic failure was defined as (1) failure to achieve at least a 1log10 decrease in plasma RNA level by 28 days of therapy; (2) failure to achieve a nondetectable plasma RNA level (<400 copies/mL) by 3 months of therapy; or (3) rebound to detectable plasma RNA levels after achieving nondetectable status. All virologic failures were based on repeat plasma RNA testing for confirmation of initial failure values.

Model of viral/cellular dynamics and statistical analyses. Biphasic plasma RNA data from baseline through 100 days of treatment were fitted into 2 viral dynamic models developed by Wu and Ding [12, 13]. These comprised a nonlinear mixedeffects (NLME) model with and without a cohort covariate. The model accommodates missing values and ignores data on the shoulder during the first day of treatment [5, 14] in an attempt to reduce possible estimation bias. Specifically, the NLME model estimates viral decay rates with the viral RNA data at baseline and day 1 removed. This model also does not require the assumption of pretreatment steady state, which may not be valid during early pediatric HIV infection and does not assume complete control of viral replication at the initiation of treatment (see Wu and Ding [13] for more details). The biphasic viral dynamic model can be written as Formula where d1 and d2 are the first-and second-phase decay rates and P1 and P2 are macroparameters (see Ding and Wu [15] and Wu and Ding [13] for more details).

For plasma RNA values below the limit of quantitation (400 copies/mL), a value of 200 copies/mL was imputed. Only the first in any series of values below quantitation was included in the model. Phase 1 decay rates (d1) for productively infected cells, phase 2 decay rates (d2) for long-lived cells, and the respective compartment half-lives (H1 and H2) were calculated. Assuming complete suppression of viral replication by the treatment regimen, the decay rates d1 and d2 are equivalent to the death rates of productively infected cells and long-lived/latently infected cells, respectively [6]. If this assumption is invalid, interpretation of the decay rates becomes more complicated [12]. The half-life of the 2 decay phases can be calculated as H1 = log(2)/d1 and H2 = log(2)/d2). The proportion of plasma virus originating from the 2 compartments can be approximated by R1 = P1/(P1+ P2) and R2 = P2/(P1+ P2). The transition time between the 2 phases is defined as tT = (log P1 log P)/(d1-d2).

The nonlinear least squares method was used to fit individual plasma RNA data. The Wilcoxon 2-sided test was used for 2sample comparisons (P values were not adjusted for multiple comparisons in this exploratory analysis). Comparisons of this method with other methodologies have been published [15, 16].

Results

Cohort characteristics and virologic response. Seventeen infants were enrolled into cohort 1, who received RTV at a dose of 350 mg/m2 every 12 h, and 31 were enrolled into cohort 2, who received RTV at 450 mg/m2 every 12 h (table 1). The mean baseline log10 plasma RNA level ranged from 4.7 to 5.6, whereas the mean baseline CD4 lymphocyte count was in the low normal range for age, from 1473 to 2261 cells/mm3.

Figure 1 A-1 C shows plasma RNA levels over the initial 6-month treatment period, comparing the 2 RTV dosing cohorts for each of the 3 enrollment age groups. All infants experienced a dramatic early (or phase 1) decline in plasma RNA level over the first 7–14 days of therapy. This was observed regardless of ultimate virologic success or failure. Declines in plasma RNA levels of 2–3 orders of magnitude were common. Subsequent to the acute phase 1 decline in plasma RNA level, a much slower phase 2 decline over several months was observed in those with sustained virologic suppression. The majority of such infants reached nondetectability by 2–4 months of age.

Table 2.

Phase 1 kinetics of productively infected cells.

Table 3.

Phase 2 kinetics of long-lived/latently infected cells.

Figure 1.

Plasma HIV RNA levels over the initial 6-month treatment period, comparing the 2 ritonavir dosing cohorts for each of the 3 enrollment age groups-6–24 months (A), 3–6 months (B), and 1–3 months (C). The X-axis represents the infant age at which samples were obtained, in months. Each line represents an individual infant.

Viral and cellular dynamics. Viral decay curves were fit to the data (NLME model), resulting in biphasic patterns. Analyses using a viral dynamic model that does not assume steady state at baseline allowed the calculation of first-and second-phase decay constants representing the death rate of productively infected (phase 1) and long lived/latently infected (phase 2) cells. Tables 2 and 3 list a series of calculated variables by age group for phases 1 and 2, respectively. Data were sufficient to allow phase 1 analyses for all 48 infants. Although the NLME model allowed phase 2 calculations for all infants, 13 of 48 had limited data due to early events such as plasma RNA rebound, reaching nondetectability before 100 days, or study dropout. Some analyses of phase 2 kinetics were performed both with and without data on these 13 infants.

Phase 1 decay constants ranged from 0.353 to 0.405 among those receiving the lower RTV dose, with no statistically significant difference between the age groups (table 2). The halflife of productively infected cells was short, calculated as being between 1.7 and 2.0. Almost all plasma virus at baseline, estimated at 98%-99%, originated from this compartment (i.e., productively infected cells). A consistent trend of a higher phase 1 decline (range, 0.566–0.593) and a shorter half-life of productively infected cells (1.2 days) was observed for the higherdose RTV cohort; the trend was statistically significant when all age groups were combined (P = .0397 for the entire cohort [n = 48]; P = .0202 for the subcohort [n = 35] after the exclusion of the 13 subjects with limited longitudinal data). The time to achieve the transition between phases 1 and 2 was ∼7 days, decreased from 11–14 days for cohort 1.

Phase 2 decay constants for the lower-dose RTV cohort ranged from 0.012 to 0.018, with relatively prolonged half-lives (range, 38–47 days) for this long-lived/latently infected cell compartment (table 3). As with phase 1 decay constants, no age-specific differences were observed. When compared with the lower-dose RTV cohort, there was a consistent trend toward higher phase 2 decay rates (range, 0.045–0.053) and shorter half-lives (range, 13–15.6 days) for the higher-dose RTV cohort; this trend was statistically significant when all age groups were combined (P = .0202 for n = 35; P = .077 for n = 48).

Correlation between viral decay rates and baseline covariates. The Spearman correlation rank test was used to test for significant correlations between individual cohort 1 and 2 decay rates (NLME model 1) and the following baseline variables: plasma RNA level, CD4 and CD8 lymphocyte counts, weight, and age. Only baseline plasma RNA level demonstrated a significant (negative) correlation with phase 1 decay rates (for cohort 1, r = −0.534 and P = .032 ; for cohort 2, r = −0.605 and P = .0009). A significant correlation between the cohort 1 phase 2 decay constant and plasma RNA was also demonstrated (r = 0.882; P = .0004), but not for the higher-dose RTV cohort (r = 0.237 ; P = .195). Once again, no age-specific effect was demonstrated for the decay rate for either phase. When both cohorts were combined in NLME model 2, which incorporated an adjustment for treatment, baseline plasma RNA level remained the only covariate with a significant correlation with the decay rates (for phase 1, r = −0.411 and P = .0048; for phase 2, r = 0.455 and P = .0018).

Discussion

The triple antiretroviral drug combination of AZT, 3TC, and RTV was well tolerated, safe, and highly active when administered to infants enrolled in PACTG 345 between 1 and 24 months of age [9]. This trial, with its intensive early monitoring, allowed for comprehensive analyses of infected cell compartment kinetics through the marker of plasma virus decay in response to this highly active antiretroviral regimen.

The kinetic parameters of viral decay (or the kinetics of infected cell compartment decay) in this study appear to be similar to values reported in adults and recently in infants [36, 8, 17]. This observation suggests that the ability to suppress viral replication in infants with potent antiretroviral regimens is comparable to that in adults. In fact, the 43% rate of durable viral suppression (defined as <400 copies/mL at 96 weeks) demonstrated in this study [9] is comparable to that in many early adult trials with protease inhibitor-containing regimens, notwithstanding the higher plateau plasma virus levels observed in infancy. The success rate of durable viral suppression and the calculation of half-lives of the long-lived cell population of 2 weeks and 3–4 weeks for the higher-and lower-dose RTV cohorts, respectively, have to be tempered by the concept of a truly latent viral reservoir. Several groups have demonstrated the persistence of replication-competent viral genomes within an additional cell compartment-quiescent memory cell populations-which makes complete viral eradication all but unattainable and will require unique therapeutic approaches [1822].

No significant age-specific differences in viral or infected cell dynamics/clearance were detected in this cohort. This is in contrast to a previous report of an age-specific decline in viral (or infected cell compartment) clearance [8]. There are several possible explanations for these divergent results: (1) both studies enrolled a relatively small number of infants for kinetic analyses; (2) potential pharmacokinetic and pharmacodynamic interpatient variability may directly affect these initial analyses; (3) somewhat different models were used in the 2 studies; and (4) the present study used a protease inhibitor-containing regimen, compared with a non-protease inhibitor regimen in the previous study. In support of the last hypothesis, pharmacokinetic data from this study (Rodman et al., unpublished data) also suggested that the RTV dose of 350 mg/m2 twice daily may be suboptimal for infants with an observed increase in drug clearance with decreasing age, leading to the enrollment of cohort 2 with a RTV dose of 450 mg/m2. Although phase 1 plasma virus decay of 3 orders of magnitude in the first 7–14 days was universally observed, regardless of ultimate outcome, up to 4 months was required to achieve undetectable plasma virus, suggesting the need for more intensive therapy (i.e., higher doses and/or additional antiretroviral agents) to accelerate the process. Several studies in adults have demonstrated that optimal long-term outcome is inversely related to the rate of achieving complete suppression of viral replication. In addition, the slower phase 2 decay in children receiving the lower dose of RTV (350 mg/m2twice daily) is probably due to the lower drug levels found in those children with good adherence, suggesting that periodic evaluation of drug levels may improve the virologic outcome.

Analyses comparing viral and infected cell dynamics between the 2 RTV dose cohorts demonstrated a significant improvement in rates of decline for both phases 1 and 2 in favor of the higher dose. Specifically, the transition time between the 2 phases was reduced by 150%-from 14 to 7 days-and the half-life of phase 2 was reduced from a range of 38 to 47 days to 13 to 16 days. Recent studies in adults have documented increasingly potent highly active antiretroviral therapy regimens that result in “better” decay rates and predict that even the most aggressive antiretroviral regimens may yet be less than “ideal” [23, 24]. Long-term follow-up and evaluation of potential low-level viral replication is indicated to determine whether the predicted benefit of more-rapid clearance of infected cell compartments translates into clinical benefit.

An intriguing study result was the negative correlation between baseline plasma RNA levels and phase 1 decay rates. This has also been demonstrated in adult cohorts and has particular significance for HIV-infected infants with high plasma RNA levels persisting through the first 1–2 years of life. The potential effect of even a small decrease in phase 1 decay for infants with high plasma virus levels is depicted in figure 2. The solid lines represent optimal (or “best case”) decay parameters for 2 individuals with high and low plasma RNA starting points. The phase 1 half-life is 1 day, with a transition to phase 2 at 7 days. The dashed line represents a 20% slower phase 1 decay constant, consistent with the negative correlation between baseline plasma RNA level and phase 1 decay. Phase 2 decay constants were assumed to be identical in all cases, with a half-life of 14 days. A likely mechanism underlying the negative correlation between baseline RNA and the phase 1 decay constant is suboptimal therapeutics-either less than maximum dosing or difficulty with adherence, or both-which may be particularly problematic in the presence of high baseline plasma RNA levels.

Figure 2.

Graphic depiction of the effect of 2 distinct phase 1 decay rates on the time to achieving nondetectable plasma virus with highly active antiretroviral treatment regimens. The solid lines represent optimal (or “best case”) decay parameters for 2 individuals with high and low plasma HIV RNA starting points. The phase 1 half-life is 1 day, with a transition to phase 2 at 7 days. The dashed line represents a 20% slower phase 1 decay constant, consistent with the negative correlation between baseline plasma RNA level and phase 1 decay. Phase 2 decay constants were assumed to be identical in all cases, with a half-life of 14 days. The horizontal dashed line represents 50 RNA copies/mL, the limit of detection.

The higher value for R2 (and related lower R1) for the higher-dose versus lower-dose cohort is intriguing but remains unexplained (tables 2 and 3). Note that R1 = P1/(P1 + P2) and R2 P2/(P1 + P2); thus, R2 = 1−R2. One possible reason is that the more-potent higher-dose treatment may result not only in faster decay rates (i.e., higher d1 and d2) but also in the conversion of some productively infected cells into long-lived infected cells that are slowly producing virus (i.e., higher R1 and lower R2).

Further intensive efforts are needed to elucidate the viruscell-drug interactions in young infants. Despite some daunting challenges (e.g., high plasma virus levels and a developing immune system), these early treatment results and kinetic analyses suggest a hopeful environment for the control of viral replication in HIV-infected infants.

Pediatric Aids Clinical Trials Group (Pactg) 345 Investigators

Protocol team members are Paula Britto and Jane Lindsey (Harvard School of Public Health); Carol Elgie (Frontier Science and Technology Research Foundation); Debbie Cloutier (Children's Memorial Hospital, Chicago); Wende Levy (PACTG Operations Center); Lynette Purdue and Elizabeth Smith (Division of AIDS, National Institute of Allergy and Infectious Diseases, National Institutes of Health [NIH]); Lynne Mofenson (National Institute of Child Health and Human Development, NIH); Nancy Pugh (Abbott Laboratories); and Neil Graham and Kathleen Zerbe (Glaxo Wellcome). Clinical trials site investigators are S. K. Burchett (Boston); P. Palumbo (Newark); R. Yogev (Chicago); K. Boyer (Chicago); J. A. Englund (Chicago); A. A. Gershon (New York City); W. Borkowsky (New York City); S. A. Spector (San Diego); V. R. Bonagura (Long Island); E. J. Abrams (New York City); S. Pahwa (Long Island);

S. R. Rana (Washington, DC); A. Kovacs (Los Angeles); E. J. McFarland (Denver); S. R. Lavoie (Richmond); J. H. Rodman (Memphis); R. Rutstein (Philadelphia); and S. S. Bakshi (New York City).

Footnotes

  • Potential conflicts of interest: E.C. has received grant support from Glaxo-SmithKline and Abbott Laboratories, participates on the speakers' bureau for Abbott, and is a shareholder in Abbott. All other authors report no potential conflicts of interest.

  • Financial support: Public Health Service (grants UO1 AI25883, AI055290, AI38855, and AI41110).

  • Study group members are listed after the text.

  • Received January 13, 2006.
  • Accepted January 29, 2006.

References

  1. 1.
  2. 2.
  3. 3.
  4. 4.
  5. 5.
  6. 6.
  7. 7.
  8. 8.
  9. 9.
  10. 10.
  11. 11.
  12. 12.
  13. 13.
  14. 14.
  15. 15.
  16. 16.
  17. 17.
  18. 18.
  19. 19.
  20. 20.
  21. 21.
  22. 22.
  23. 23.
  24. 24.
| Table of Contents