Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, March 2002, p. 746-754, Vol. 46, No. 3
0066-4804/02/$04.00+0 DOI: 10.1128/AAC.46.3.746-754.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Worldwide Clinical Pharmacology, GlaxoSmithKline Inc., Research Triangle Park, North Carolina
Received 12 January 2001/ Returned for modification 1 August 2001/ Accepted 26 November 2001
|
|
|---|
|
|
|---|
A strong correlation between antiviral activity and the concentrations of PIs in plasma has been demonstrated in several studies (1, 7, 8, 18, 26, 28). Subtherapeutic concentrations of PIs in plasma can be associated with viral rebound and the development of drug resistance, indicating the importance of maintaining consistent drug concentrations in plasma (1, 7, 8, 18, 26, 28). A regimen that results in increased trough exposures could potentially provide therapy for patients with drug-resistant viruses and make once-daily (q.d.) dosing a feasible therapeutic option.
Amprenavir (Agenerase; formerly 141W94) is a recently introduced HIV PI. Like the other approved PIs, amprenavir is metabolized primarily by the hepatic cytochrome P450 enzyme system, specifically, CYP3A4 (J. Woolley, S. Studenberg, C. Boehlert, G. Bowers, A. Sinhababu, and P. Adams, Abstr. 37th Intersci. Conf. Antimicrob. Agents Chemother., abstr. A-60, 1997). All HIV PIs are substrates and inhibitors of CYP3A4, but the HIV PI ritonavir is the most potent inhibitor of CYP3A4 (Woolley et al., 37th ICAAC). Previous studies have indicated significant increases in drug exposure from coadministration of ritonavir with indinavir, saquinavir, and nelfinavir (12, 13, 18, 19, 22, 23, 30; D. Burger, P. Hugen, H. Hofstede, P. Koopmans, M. Stek, Jr., P. Reiss, and J. Lange, Abstr. 37th Intersci. Conf. Antimicrob. Agents Chemother., abstr. 321, 1999; Woolley et al., 37th ICAAC). The magnitude of pharmacokinetic interaction between ritonavir and amprenavir on coadministration was recently evaluated in a series of three randomized, two-sequence, multiple-dose studies with healthy subjects (25). Those studies demonstrated that ritonavir is a potent inhibitor of amprenavir metabolism and that addition of ritonavir (100 or 300 mg every 12 h) to amprenavir (450 or 900 mg every 12 h) resulted in clinically and statistically significant increases in the amprenavir area under the plasma concentration-versus-time curve (AUC) from time zero to 24 h (AUC0-24) and minimum concentration in plasma (Cmin).
In this report, we define a model using data from the first two studies (450 mg of amprenavir with 100 or 300 mg of ritonavir) and NONMEM for inhibition of amprenavir elimination by ritonavir. We then validated this model by prospectively estimating the pharmacokinetics of the combination of 900 mg of amprenavir plus 100 mg of ritonavir. Amprenavir data from all three interaction studies were used to develop the final model. This model was then used to assess different dosing regimens based on simulated Cmins by using Monte Carlo simulation for twice-daily (b.i.d.) and q.d. regimens and to determine regimens for clinical investigation with subject populations whose viruses have different susceptibilities to HIV PIs.
(This work was presented in part at the 7th Conference on Retroviruses and Opportunistic Infections, San Francisco, Calif., 30 January to 2 February 2000 [B. M. Sadler, P. J. Piliero, S. L. Preston, and D. S. Stein, Abstr. 7th Conf. Retrovir. Opportunistic Infect., abstr. G302e, 2000] and the 1st International Workshop on Clinical Pharmacology of HIV Therapy, Noordwijk, The Netherlands, March 2000 [M. Sole, D. S. Stein, and B. Sadler, Abstr. 1st Int. Workshop Clin. Pharmacol. HIV Ther., abstr. 5.2, 2000].)
|
|
|---|
In each study, 18 subjects were randomly assigned to receive one of two treatment sequences (nine subjects per sequence) in a two-period, one-way crossover design. In the first study (PRO10017), 450 mg of amprenavir every 12 h (b.i.d.) alone (sequence 1) or 300 mg of ritonavir b.i.d. alone (sequence 2) was administered for 7 days. For the following 7 days, all subjects received both amprenavir and ritonavir at these same doses. In the subsequent two studies, the study design remained the same but a regimen of 450 mg of amprenavir b.i.d. and 100 mg of ritonavir b.i.d. (PRO10022) or a regimen of 900 mg of amprenavir b.i.d. and 100 mg of ritonavir b.i.d. (PRO10023) was investigated. All data from both sequences were used in the analysis.
All subjects provided written informed consent to participate in the trials. All study protocols were reviewed and approved by the institutional review board at each of the study sites.
Blood sampling and amprenavir and ritonavir assays. Blood samples for the determination of plasma amprenavir and ritonavir concentrations were collected on days 1, 7, and 14. Seventeen samples were collected between predosing and 24 h postdosing on day 1, and 15 samples were collected between predosing and 12 h postdosing on days 7 and 14. Details of the sampling and analysis methods have been described previously (25).
Data analysis. Traditional noncompartmental analyses of the data and clinical safety evaluations are reported in a separate publication (25).
Pharmacokinetic modeling. Plasma concentration-time data from each treatment sequence in each study were used to fit a standard linear compartmental model by using the nonlinear mixed-effects model software package NONMEM (version 5) (27). NONMEM was used to describe the interaction between amprenavir and ritonavir and simulate selected amprenavir pharmacokinetic parameters (maximum concentration of drug in plasma [Cmax], AUC0-24, and Cmin) for b.i.d. and q.d. dosing regimens of amprenavir in combination with ritonavir. NONMEM was developed independently for amprenavir and ritonavir. Models were built by forward addition at the P < 0.05 significance level and backward elimination at the P < 0.01 significance level by a log likelihood test.
Amprenavir model.
A standard two-compartment, first-order absorption model was found to describe the data well. The final pharmacokinetic model for amprenavir included the effects of the observed ritonavir concentration on the elimination rate. The central volume of distribution (V/F) was found to be a function of weight and
1-acid glycoprotein (AAG) concentration. All parameters included log-normal interindividual error terms. The absorption rate constant (Ka), V/F, and the transfer rate constants between the central and peripheral compartments (K23 and K32) were modeled separately for days 1 and 14. It cannot be determined if these effects are dependent on time or on the administration of ritonavir. The effect on Ka, V/F, K23, and K32 was small. The equations for the parameters for the amprenavir model are given below:
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
1,
2, and
3 are log-normal interindividual error terms. Weight (WT) was expressed in kilograms. The maximum effect value (Emax;
6) was constrained to be less than or equal to 0.99.
Ritonavir model.
Because ritonavir concentrations were found to be a significant factor in the final amprenavir model, a population-based model of ritonavir concentrations was also developed and provided input to the amprenavir model during Monte Carlo simulations. Pharmacokinetic data obtained after 14 days of dosing only were used to develop the ritonavir model. This should reflect steady-state dosing. In the initial model developed with data from the first two studies, ritonavir pharmacokinetics were also well described by a two-compartment linear model with first-order absorption and an absorption lag time. In the final model (data from all three studies), a second compartment was not statistically supported by the data. The final pharmacokinetic model for ritonavir included the effects of dose on Ka and bioavailability (F1; relative bioavailability of the 300-mg dose compared to that of the 100-mg dose), the effects of age on K, and the effects of sex and weight on the volume of distribution (V). Ka, K, and V included log-normal interindividual error terms (
1,
2, and
3, respectively). It should be noted that the effect of dose on F could as well have been applied to V/F. However, it was believed to be biologically more reasonable to have a dose effect on F1 than on V/F. The results of the subsequent simulation would have been identical in either case. The equations for the ritonavir model are given below:
![]() |
![]() |
![]() |
![]() |
![]() |
) was best expressed in terms of a constant coefficient of variation (
1) of the predicted concentration () plus a simple additive error term (
2). Monte Carlo simulation. Simulations were performed by using a composition of the amprenavir and ritonavir models. The ritonavir model was used to simulate ritonavir concentrations, which where then used to calculate the amprenavir K. Weight and AAG concentration were simulated from distributions derived from the study data. Sex was equally distributed between men and women. The value of Dose300 was 0 for 100 mg, 0.5 for 200 mg, and 1 for 300 mg. Pharmacokinetic parameters from day 14 were used for the amprenavir model.
Sensitivity analyses. The concentration of ritonavir that results in 50% of the maximum inhibition (Km) predicted by the model was 0.0122 µg/ml. This is below the smallest observed concentration of ritonavir (0.028 µg/ml) in any of the three studies, suggesting that the inhibition of elimination of amprenavir was essentially saturated in all subjects for the duration of the b.i.d. dosing interval. Because of the lack of data near the Km, this value was in question. Log likelihood plots were constructed to test the range of values for Km that were consistent with the observed data.
Model validation. The initial model was developed with data from studies PRO10017 and PRO10022. These studies included doses of 450 mg of amprenavir plus 100 mg of ritonavir b.i.d. and 450 mg of amprenavir plus 300 mg of ritonavir b.i.d.. The results from this analysis were then used to predict the outcome of study PRO10023 (900 mg of amprenavir plus 100 mg of ritonavir b.i.d.). When the simulations and prediction for PRO10023 were done, the study had not yet been completed and no data were available. The two endpoints of interest were AUC0-24 and Cmin. The AUC0-24 predicted by the simulations was within 3.9% of the observed value (48.95 ng · h/ml predicted versus 47.71 ng · h/ml observed), and the Cmin was within 19.4% (2.95 ng/ml predicted versus 2.47 ng/ml observed). The final simulations were based on the analysis of data from all three studies.
Cmin-to-protein binding-adjusted IC50 ratio central tendency estimates. Amprenavir Cmins obtained from simulated pharmacokinetic parameters for the b.i.d. and q.d. dosing regimens of amprenavir in combination with ritonavir were compared with protein binding-adjusted 50% inhibitory concentrations (IC50s) for clinical viral isolates from PI-naïve and -experienced patients to determine if the amprenavir Cmin would exceed the IC50s for the viral isolates.
The distribution of the ratios of Cmin to IC50 in the population was calculated by resampling from the observed IC50s for the viral isolates and the observed Cmins (for 1,200 mg b.i.d.) or the simulated Cmins (for 600 mg of amprenavir plus 100 mg of ritonavir b.i.d. and 1,200 mg of amprenavir plus 200 mg of ritonavir q.d.). The ratios were adjusted for the
90% protein binding of amprenavir. Phenotypic assays for IC50 were performed by Virco NV (Mechelen, Belgium) by the method of Hertogs et al. (11) to determine the amprenavir IC50 for the viral isolates. The three populations used were wild-type isolates (n = 334; IC50 = 14.6 ± 12.5 ng/ml), isolates from patients clinically described as first-time amprenavir failures (n = 36; IC50 = 61.4 ± 83.9 ng/ml), and isolates from patients described clinically as multiple PI failures (n = 328; IC50 = 90.3 ± 84.6 ng/ml).
|
|
|---|
Pharmacokinetic modeling.
The initial model, based on data from the first two studies, was revised with the data from the third study. In the full data set, a two-compartment model for ritonavir was not statistically supported by a log likelihood test, and so a one-compartment model was used. A two-compartment model was once again supported for amprenavir. The pharmacokinetic parameter estimates of the models for amprenavir and ritonavir are presented in Table 1.
In the amprenavir model, weight and AAG were found to have an effect on V. K was found to be dependent on the (observed) ritonavir concentration by the (Emax) relationship:
![]() |
|
View this table: [in a new window] |
TABLE 1. Pharmacokinetic parameter estimates
|
|
View this table: [in a new window] |
TABLE 2. Interindividual error terms (OMEGAs)
|
![]() View larger version (9K): [in a new window] |
FIG. 1. Likelihood plot for values of Km in the amprenavir model from all three studies. The dotted horizontal line indicates the line where P is equal to 0.05. For values of Km in which the log likelihood is above this line, the observed data are less than 5% as likely to occur as they would be with the minimum value of Km.
|
![]() View larger version (16K): [in a new window] |
FIG. 2. Sensitivity plots for different values of Km versus the predicted Cmin for 900 mg of amprenavir q.d. (A) and 1,200 mg of amprenavir q.d. (B) in combination with 100, 200, or 300 mg of ritonavir q.d., based on the model from all three studies.
|
![]() View larger version (16K): [in a new window] |
FIG. 3. Geometric mean simulated amprenavir (APV) concentration-versus-time curves at steady state for 450, 600, and 900 mg of amprenavir b.i.d. in combination with 100 mg of ritonavir (RTV) b.i.d.
|
![]() View larger version (12K): [in a new window] |
FIG. 4. Geometric mean simulated amprenavir concentration-versus-time curves at steady state for 900 and 1,200 mg of amprenavir (AMP) q.d. in combination with 200 mg of ritonavir (RTV) q.d.
|
|
View this table: [in a new window] |
TABLE 3. Simulated and observed exposures from a two-compartment NONMEM of amprenavir and ritonavir or amprenavir alonea
|
![]() View larger version (32K): [in a new window] |
FIG. 5. Central tendency ratios of amprenavir (APV) Cmin to protein binding-adjusted IC50 for amprenavir b.i.d. alone (APV) or in combination with ritonavir (APV/R) for wild-type and resistant clinical isolates. Viral isolates of the wild type (WT; n = 334), from patients with first-time amprenavir failures (APV failures; n = 36), and from patients with multiple PI failures (multi PI failures; n = 284) were tested. The central line is the median of the range, the box goes from the 25th to the 75th percentiles of the range, and the lines with bars represent the 10th and 90th percentiles of the range.
|
|
|
|---|
The Km estimated in our model for ritonavir was 0.0122 µg/ml (total concentration). This is in reasonable agreement with the range of Km or Ki values reported for ritonavir in in vitro investigations under various experimental conditions (14). Koudriakova et al. (14) calculated a ritonavir Km of 0.068 µM (total concentration, 0.034 ng/ml) for expressed CYP3A4 and a value of 0.063 ± 0.045 µM in microsomes. They noted that ritonavir's inhibitory effect lasted well after its removal, indicating that ritonavir inactivated CYP3A4. Kumar et al. (16) calculated a Km of 0.08 µM (0.040 ng/ml) for the formation of ritonavir's primary CYP3A4 metabolite, M1. Ritonavir has a Ki of 0.013 µM (0.0066 ng/ml) for inhibition of metabolism of ABT-378, another HIV PI (15), and has Ki values of 0.017 µM (0.008 ng/ml) for terfenadine (16) and 0.019 µM (0.009 ng/ml) for testosterone (6), which are classic CYP3A4 substrates. All concentrations are given as total concentrations. The level of protein binding of ritonavir in plasma is high (
98%). The protein binding in the medium used for these in vitro estimates of Km and Ki is not known. However, differences in free drug concentration could bias these comparisons.
As the calculated Km affects the model predictions of different regimens, we conducted a sensitivity analysis to further explore this source of variability. This sensitivity analysis indicated that for a wide range of Km values the greatest effect (sensitivity to Km) occurred in the simulated q.d. regimens. It is for this reason that we would not currently recommend the use of q.d. regimens with less than 200 mg of ritonavir.
The pharmacokinetic profiles of PIs show variability among different patients, leading to low Cmins in some patients (4, 7, 20). Multiple studies have found a relationship between Cmin or average exposure (i.e., AUC) and antiviral activity or time to the development of resistance for PIs. Prior studies of indinavir, saquinavir, and amprenavir have all demonstrated that the relationships of AUC or Cmin to antiviral activity are well described by a sigmoid Emax curve. For both indinavir and amprenavir the plateau of the relationship occurs at average exposures at or above the protein binding-adjusted average IC50 for the viral population. This suggests that maintenance of a concentration above the IC50 could lead to a better overall antiviral effect and possibly better long-term suppression of virus. The lower limit of the 95% confidence interval for the amprenavir concentration in the regimen of 600 mg of amprenavir plus 100 mg of ritonavir b.i.d. was 0.264 µg/ml. This value is above the mean (protein binding-adjusted) IC50 for wild-type virus examined in this study (0.146 µg/ml).
In the traditional analysis of the three pharmacokinetic drug interaction studies (25), the addition of 100 or 300 mg of ritonavir b.i.d. to 450 mg of amprenavir b.i.d. resulted in statistically significant increases in the amprenavir Cmin and AUC compared with those achieved when amprenavir was administered alone. The modeling clearly showed that this is due to a decrease in elimination of amprenavir. In the present study, sensitivity analyses with simulations of the effect of the predicted amprenavir Cmin on variable ritonavir Km values calculated that the maximum reduction in amprenavir clearance was 66%. A 66% magnitude of inhibition was predicted to be maintained throughout a 24-h dosing interval if the Km value is <0.001, but no inhibition would occur throughout this interval if Km was >0.1. Between these values of Km (which includes the point estimate for Km), the 200-mg dose of ritonavir would be expected to sustain a greater and longer duration of inhibition of amprenavir clearance than the 100-mg dose of ritonavir, but it would be less than that seen with the 300-mg dose. This expectation is highly dependent on the estimate of the concentration at the site of effect that causes 50% of the maximal effect.
Several sources of potential bias exist in this analysis. First, the assumption of the model was an Emax inhibition of the elimination of amprenavir by ritonavir. This model fit the data best. However, in vitro work has suggested a persistent or even irreversible inhibition of cytochrome P450 CYP3A4 activity. Attempts to construct a model based on this assumption were unsuccessful. These models were unsuccessful either because they were not numerically stable or the model did not fit the data as well as a traditional Emax model did. It should be noted that if the effect of ritonavir were persistent rather than instantaneous, as assumed in this model, the results presented here would underestimate the magnitude of the effect for the q.d. dosing. The actual Cmin would be expected to be even higher than the concentration suggested by this model, presumably with more antiviral activity. Preliminary clinical data obtained with a regimen of 1,200 mg of amprenavir and 200 mg of ritonavir q.d. in HIV-infected subjects indicate that Cmins are higher than those predicted by the conservative assumptions of our model, suggesting some degree of persistence of the effect of ritonavir on amprenavir clearance in vivo (R. Wood, C. Trepo, J. M. Livrozet, K. Arasteh, J. Eron, P. Kaur, O. Naderer, and M. B. Wire, 8th Conf. Retrovir. Opportunistic Infect., poster 332, 2001).
Another potential source of bias is the possibility that amprenavir affects the metabolism of ritonavir. Given the study design, this effect could not be distinguished from a time effect on ritonavir clearance, which is known to occur. However, if the Km of the interaction is as small as is suggested by these analyses, this effect should be small since the effect of ritonavir on amprenavir is already nearly saturated. Finally, the clinical studies suggested that the 100- and 300-mg b.i.d doses of ritonavir had very similar effects on amprenavir clearance. In fact, study PRO10017 suggested that the concentrations of amprenavir when given with 100 mg of ritonavir b.i.d. would be higher than those seen when amprenavir was given with 300 mg of ritonavir. The model used here (Emax inhibition of clearance) would not permit a lower dose (100 mg) to have a larger effect than the higher dose (300 mg). Instead, the model suggests that the concentrations achieved with either dose nearly saturate the inhibition. Hence, in the final simulations, the concentrations simulated with the 300-mg dose of ritonavir are slightly higher than those seen with the 100-mg dose.
It should be noted that healthy volunteers have a fairly constant AAG concentration, and that concentration is significantly lower than those in HIV-positive patients. The model suggested that AAG influences V. This is likely related to the binding of amprenavir to AAG, reducing the free fraction available for diffusion into tissues outside the plasma. No effect of AAG is seen on elimination, which is where the effect of ritonavir is seen. Therefore, we believe that differences in AAG concentrations should not preclude extrapolation of these results to HIV-positive patients. Other studies have shown that the AAG concentration does not affect the free drug concentration (but does affect the total drug concentration and free fraction [24]).
In summary, NONMEM analysis and Monte Carlo simulation were used to predict the pharmacokinetic parameters of amprenavir in the presence of ritonavir. A combination of 600 mg of amprenavir plus 100 mg of ritonavir b.i.d. is predicted to result in clinically and statistically significant increases in amprenavir Cmins. This combination has the potential to reduce the pill count from 16 to 10 capsules per day, which may reduce the adherence problems observed with some patients. In addition, the potential exists to enhance efficacy in salvage regimens, particularly against PI-resistant virus, as a result of the increase in Cmin above that obtained with the 1,200-mg dose of amprenavir b.i.d. alone. In the case of patients who have failed a traditional amprenavir regimen, these simulations suggest that a regimen of 600 mg of amprenavir with 100 mg of ritonavir b.i.d. would restore the Cmin-to-IC50 ratio to a value similar to that achieved with 1,200 mg of amprenavir b.i.d. alone for wild-type viruses. Regimens of 900 or 1,200 mg of amprenavir q.d. in combination with 200 mg of ritonavir q.d. were also predicted to give trough concentrations comparable to those achieved with the approved 1,200-mg dose of amprenavir b.i.d. alone. These regimens could offer the potential of a further reduction in pill burden and the advantage of less frequent dosing. Clinical studies are in progress to evaluate the pharmacokinetics, efficacies, and safety profiles of different doses of amprenavir in combination with ritonavir for b.i.d. and q.d. regimens. Clinicians should exercise appropriate caution if they are basing their dosing regimens for HIV-positive subjects on a model derived from studies conducted with healthy volunteers.
|
|
|---|
1-acid glycoprotein on pharmacokinetics of amprenavir, a human immunodeficiency virus protease inhibitor. Antimicrob. Agents Chemother. 45:852-856.
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»