Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, November 2006, p. 3548-3555, Vol. 50, No. 11
0066-4804/06/$08.00+0 doi:10.1128/AAC.00943-05
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Université Paris-Descartes, Faculté de Médecine,1 Assistance Publique-Hôpitaux de Paris,2 Pharmacologie Clinique,3 Gynécologie-Obstétrique 1,4 Groupe Hospitalier Cochin-Saint-Vincent de Paul,5 INSERM,6 Virologie,7 Unité d'Immunologie-Hématologie Pédiatrique,8 Groupe Hospitalier Necker-Enfants Malades,9 Centre René-Huguenin Saint-Cloud,10 EA3620, Paris, France,11
Received 25 July 2005/ Returned for modification 23 November 2005/ Accepted 10 August 2006
|
|
|---|
|
|
|---|
|
|
|---|
Analytical methods. The assay for lopinavir was performed by using high-pressure liquid chromatography with UV detection. Briefly, 500 µl of sodium carbonate (0.2 M) and 6 ml of ethyl acetate-hexane (50:50 vol/vol) was added to 200 µl of plasma. After 20 min of mixing, the supernatant was evaporated at 30°C under a stream of nitrogen. Dry residues were then reconstituted with 200 µl of the mobile phase (tetramethylammonium perchlorate [0.01 M] in trifluoroacetic acid [0.01%]-acetonitrile [39:61; vol/vol]), and 100 µl of this mixture was injected into the chromatographic system. The separation was performed on a Nucleosil C8 column (125 by 4.6 mm, 3 µm) at a flow rate of 1 ml/min. Detection was performed at 205 nm. The quantification limit of the method was 0.10 mg/liter, with the interassay precision and bias being less than 7% in the calibration range of 0.1 to 20 mg/liter.
Population pharmacokinetic modeling.
Concentration-time data were analyzed by use of the first-order conditional estimation method with interaction of the nonlinear mixed effects modeling program NONMEM (version V, level 1.1, double precision) (1). Several structural pharmacokinetic models were investigated. Classical one- and two-compartment models with first- and zero-order absorption were first evaluated. A one-compartment pharmacokinetic model in which a single rate constant (k) was used for both the absorption and the elimination processes was also tried (31). The equation for this model was as follows:
![]() |
is the time that elapsed between the administration of two doses (i.e., 12 h), C is concentration, t is time, F is bioavailability, D is lopinavir dose, and V is volume. The explicit solution for this pharmacokinetic model was coded in the $PRED section of the control stream, and its parameters were the apparent total clearance (CL/F) and the volume of distribution (V/F) of lopinavir. Several error models (i.e., proportional, exponential, and additive random effects model) were also investigated as means of describing interpatient and residual variabilities.
Systematic testing for the influence of continuous covariates on the pharmacokinetic parameters was done by use of a generalized model, according to the following equation, by using, for example, CL/F and BW: CL/F = TV(CL/F) · (BW/median BW)
BW, where TV(CL/F) is the typical value of the apparent clearance for a patient with the median covariate value, and
BW is the influential factor for body weight. Binary covariates (gender, combined treatment, dosage form) were investigated as follows: CL/F = TV(CL/F) ·
SEX, where SEX is equal to 1 for males and 0 for females.
The possible effect of a covariate on lopinavir bioavailability was investigated as described above by estimation of the same influential factor on both CL/F and V/F.
The significance of a relationship between a pharmacokinetic parameter and a covariate was assessed by use of the chi-square test of the difference between the objective functions of the basic model (without the covariate) and the model with the covariate. A covariate was retained in the model if it produced a minimum decrease in the objective function of 4 units (P < 0.05, 1 degree of freedom) and if one of the following criteria was satisfied: (i) it led to a reduction of the interindividual variability (
) of the associated pharmacokinetic parameter or (ii) if its effect was biologically plausible. An intermediate multivariate model that included all selected covariates was then obtained. A covariate was retained in the final multivariate model if its deletion from the intermediate model led to a 7-point increase in the objective function (P < 0.01, 1 degree of freedom). At each step, the goodness of fit was evaluated by use of a graph of the weighted residuals versus time after administration of the dose (time) or by use of a graph of the weighted residuals versus the predicted concentration.
The accuracy and robustness of the final population model were assessed by a bootstrap method, which consisted of repeated random sampling with replacement from the original data set. This resampling was repeated 1,000 times, and the values of the parameters estimated from the bootstrap set were compared to the estimates obtained from the original data set. The entire procedure was performed in an automated fashion by using Wings for NONMEM (24).
Individual Bayesian estimates of the pharmacokinetic parameters were used to calculate the individual area under the concentration-time curve from time zero to 12 h (AUC0-12) and the trough concentration (Ctrough). The AUC0-12 achieved for the recommended pediatric dosage regimen was also calculated.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Characteristics of the 157 children (67 girls, 90 boys) at sampling time
|
![]() View larger version (14K): [in a new window] |
FIG. 1. Age distribution at the time of retrieval of the first blood sample in the population of the study.
|
(where
is the interindividual variability of V/F), suggesting that these parameters could not be estimated (data not shown). Because of these results, a limited form of this model that corresponded to a one-compartment model with first-order absorption when ka and kel are not distinguishable was used (31). This model provided the best fit, with a further 20-point decrease in the objective function compared to the value obtained by use of the classical one-compartment model. The one-compartment model with zero-order absorption led to pharmacokinetic parameter estimates that were in disagreement with previously published values (data not shown), and the two-compartment model seemed overparameterized for our data, as it systematically led to convergence failure. The graph of the observed concentrations as a function of time after dosing along with the typical pharmacokinetic curve is shown in Fig. 2. Interpatient variability was described by an exponential error model, whereas residual variability was described by a combined exponential and additive error model. A significant covariance term was found between the
values for CL/F and V/F.
![]() View larger version (20K): [in a new window] |
FIG. 2. Observed lopinavir concentrations and typical pharmacokinetic curve.
|
|
View this table: [in a new window] |
TABLE 2. Effect of the tested covariates on the objective functiona
|
Table 3 summarizes the population parameter estimates. The goodness of fit was also evaluated graphically from the distribution of the points on graphs of the weighted residuals versus time and weighted residuals predicted concentration (Fig. 3A and B).
|
View this table: [in a new window] |
TABLE 3. Population pharmacokinetic parameters of lopinavir in 157 children and bootstrap validationa
|
![]() View larger version (17K): [in a new window] |
FIG. 3. Goodness of fit visualized on the weighted residuals (WRES) versus time after dose (A) and weighted residuals versus predicted plasma lopinavir concentrations (PRED) (B).
|
Evolution of lopinavir exposure with respect to BW, age, and sex. The different values of the influential factor of body weight on CL/F and V/F indicated the nonlinear pattern of the lopinavir half-life (t1/2) with respect to body weight. These nonlinear relationships between CL/F, weight-normalized CL/F, t1/2, and body weight are displayed in Fig. 4. As seen in Fig. 4, the weight-normalized CL/F was even higher in the neonatal subgroup (i.e., age <8 days), with a mean value of 0.44 liter/h/kg (six neonates). Eleven children, including the six neonates, had BWs less than 7 kg at the time of sampling and had a mean age and a mean BW of 2.3 months and 3.7 kg, respectively. The mean CL/F, V/F, and t1/2 values for these children were 1.15 liters/h (i.e., 0.35 liter/h/kg), 6.70 liters, and 3.9 h, respectively. In this subpopulation, the 12-mg/kg lopinavir dose that is recommended for children with BWs between 7 and 14 kg would provide mean AUC0-12, maximum concentration, and Ctrough values of 38.6 mg · h/liter, 3.68 mg/liter, and 2.16 mg/liter, respectively. If the pharmacokinetic target is the mean AUC0-12 calculated for our population (i.e., 108 mg · h/liter), a 34-mg/kg dose would be necessary. The corresponding calculated mean Ctrough would be 6.11 mg/liter.
![]() View larger version (15K): [in a new window] |
FIG. 4. Evolution of lopinavir plasma clearance (A), weight-normalized clearance (B), and terminal half-life (C) according to BW.
|
![]() View larger version (29K): [in a new window] |
FIG. 5. Calculated (A) mean lopinavir AUC0-12 and (B) steady-state Ctrough values achieved for the recommended pediatric dosage regimen in boys and girls with respect to BW (calculated values and mean curves).
|
|
View this table: [in a new window] |
TABLE 4. Lopinavir efficacy after 12 months of treatment in pretreated children by age and sexa
|
|
|
|---|
For neonates and infants, our model also showed the decrease in the lopinavir half-life in children with BWs of less than about 10 kg that was recently reported for 12 infants whose ages were between 14 days and 6 months (E. G. Chadwick, J. H. Rodman, P. Palumbo, D. Persaud, J. Chen, J. Gardella, K. Luzuriaga, R. Yogev, P. Emmanuel, M. Rathore, and the PACTG 1030 Study Team, Abstr. 12th Conf. Retrovir. Opportunistic Infect., abstr. 766, 2005). Indeed, for a median BW of 5.1 kg, that study reported median CL/F and V/F values of 1.48 liters/h and 7.74 liters, respectively. If a one-compartment model is assumed, these parameters allow calculation of an elimination half-life of 3.6 h, which is consistent with our own results (Fig. 4). Such a result is surprising, as lopinavir is extensively metabolized by cytochrome P450 3A4 (CYP 3A4) (15), the activity of which is known to be markedly reduced in the first weeks of life (18). A first possible explanation would be a decrease in lopinavir bioavailability secondary to absorption disorders. However, a change in bioavailability would have the same consequences on V/F and CL/F, and the elimination half-life would remain constant throughout the BW range of the study. Another possibility could be the presence, at least during the first weeks of life, of CYP 3A7 (4). This enzyme, which has a 90% amino acid sequence similarity with CYP 3A4 (16), could be involved in lopinavir metabolism; and the pharmacokinetic enhancement by ritonavir could be weaker with this cytochrome than with CYP 3A4. The pharmacokinetic enhancement by ritonavir could also be reduced by a decrease in ritonavir bioavailability in children younger than 6 months of age. Indeed, a previously published study investigating the pharmacokinetics of ritonavir in 48 children between the ages of 6 months and 14 years found that the two children who were 6 months of age had a ritonavir CL/F almost 10 times higher than the CL/F found in older children and imputed this result to a decrease in ritonavir bioavailability (23). Although this result was obtained with high ritonavir doses, this phenomenon could also be observed with the low ritonavir dose of the lopinavir/ritonavir combination and explain the increase in lopinavir CL/F that we observed for children younger than 6 months of age. The metabolic immaturity could also be balanced by an increase in the lopinavir free fraction (fu) in the youngest children. Indeed, lopinavir is known to be highly bound to
1-acid glycoprotein (AAG) (5), the concentration of which is markedly reduced at birth (about 0.2 to 0.3 g/liter) and gradually increases to reach adult values (0.7 to 1.0 g/liter) after 9 to 12 months of life (22). Ritonavir is also highly bound to plasma proteins, including AAG (14); but as the plasma ritonavir concentration is much lower than the lopinavir plasma concentration (most of the plasma ritonavir concentrations were less than 0.5 mg/liter in our population), a decrease in the AAG concentration could have more consequences on lopinavir binding than on ritonavir binding. Therefore, more free lopinavir than free ritonavir could be available for CYP 3A4, which could partially counteract the competitive inhibition by ritonavir.
So, it seems likely that this increased ability to eliminate lopinavir in these children is, in fact, a consequence of a decrease in the pharmacokinetic enhancement provided by ritonavir. Such a phenomenon could therefore be found for the pharmacokinetics of other boosted PIs in neonates and infants. Unfortunately, to our knowledge, no information about the pharmacokinetics of boosted PIs in this particular population is available. Such studies are therefore warranted to confirm this hypothesis. These changes in lopinavir pharmacokinetics in neonates and infants could have consequences on treatment efficacy in this subpopulation. Indeed, because of the decrease in the elimination half-life, it will be difficult to achieve the Ctroughs or genotypic inhibitory quotients that are usually accepted as target values (6, 17). The 12-mg/kg BID lopinavir dose that is recommended for children with BWs greater than 7 kg could be inappropriate for the youngest children, as this dose provides a mean Ctrough of 2.5 mg/liter, which was found to be the limit below which treatment failure is likely, whatever the lopinavir mutation score at the baseline is (6). Our results indicate that a 34-mg/kg BID dose would provide a lopinavir exposure close to the mean reported exposure in children. However, this dose is at least three times higher than the recommended pediatric dose and could therefore raise a safety issue. Furthermore, this theoretical 34-mg/kg dose was calculated by assuming that the pharmacokinetics of lopinavir are linear with the dose, which is probably not the case. Finally, because a change in fu is thought to have no clinically relevant consequence for metabolized and orally administered drugs, as the concentration of the free drug (i.e., the active form) remains constant (2), a dose increase could be irrelevant if the change in the pharmacokinetics of lopinavir are a result of modified protein binding, as suggested above. It therefore seems important to investigate the active form of lopinavir and to specifically determine the optimal lopinavir dosing regimen for neonates and infants in further studies. For example, an alternative, given the decrease in the lopinavir half-life, could be the administration of lopinavir three times a day instead of BID.
The sex-related change in the lopinavir CL/F that was found in adolescents has an important consequence on lopinavir exposure, as, for the current recommended lopinavir regimen, the AUC0-12 calculated for the oldest boys (i.e., 50 mg · h/liter) is approximately half the AUC0-12 calculated for the oldest girls (i.e., 100 mg · h/liter). Poorer compliance by boys than by girls could have been a possible bias that would explain our results. The omission of the previous doses is indeed impossible to investigate when a single blood sample is drawn after the drug was taken during the consultation. However, such a link between sex and compliance would have provided a relationship between sex and V/F similar to the one that was found between sex and CL/F, which was not the case (Table 2). The previously published pediatric pharmacokinetic studies did not find such a sex-related difference, probably because of the lack of children older than 12 years of age (27) or the small number of patients older than 12 years of age included (3). A previously published study performed with adults (18 men and 7 women) also investigated the sex-related differences in the pharmacokinetics of lopinavir (25). Even though the median lopinavir AUC0-12 seemed greater in women (90.6 mg · h/liter) than in men (74.1 mg · h/liter), the difference was not significant. As stated above, lopinavir is extensively metabolized at the hepatic level by CYP 3A4, a metabolic pathway that is not known to display clinically relevant sex-related differences in adults (8, 28). Our results could therefore indicate that adolescence could be characterized by pharmacokinetic specificities compared to the pharmacokinetics in both younger children and adults, so that the current pharmacokinetic assumption that places adults and adolescents in the same category should be reconsidered, at least for metabolized drugs. This was previously suggested for CYP 1A2, as the rate of demethylation of caffeine decreased in early puberty for girls and late puberty for boys (20). However, it must be noticed that similar sex-related differences in the pharmacokinetics of indinavir (10) and saquinavir (25), two other PIs metabolized by CYP 3A4, have been reported in adults. Unfortunately, the pharmacokinetic studies performed with children did not include enough adolescents to investigate the possible sex-based differences in indinavir and saquinavir pharmacokinetics in this population (7, 12, 13, 19). Furthermore, a result that could support a sex-related difference in the lopinavir CL/F in adults is that our model seemed to predict the values of the lopinavir pharmacokinetic parameters in adults quite accurately. Indeed, the study cited above reported median lopinavir AUC0-12s of 90.6 and 74.1 mg · h/liter for women and men, respectively (25). As the lopinavir dose was 400 mg BID, these AUC0-12 values corresponded to CL/F values of 4.41 and 5.40 liters/h for women and men, respectively. The inclusion of the mean BW of 70 kg that was reported in this study in our CL/F equations gave values of 4 liters/h for women and 5.6 liters/h for men. Another previously published pharmacokinetic study performed with adults reported mean CL/F and V/F values of 5.73 liters/h and 61.6 liters, respectively (9). As the population in that study was characterized by a mean BW of 72 kg and comprised 85% men, it was possible to calculate respective CL/F and V/F values of 5.6 liters/h (by using the CL/F equation for boys older than age 12 years) and 49.8 liters, respectively, by the use of our model. So, the consistency between the values calculated by the use of our model and the values reported in previous studies could suggest a sex-related difference in lopinavir pharmacokinetics in adults. Besides, it also indicates that our model accurately predicts the lopinavir CL/F and V/F over a very large BW range (i.e., from the BWs of neonates to the BWs adults), which supports the structural pharmacokinetic model used.
Our data did not allow us to investigate rigorously the possible consequence of this increase in the lopinavir CL/F on treatment efficacy, so the need to give higher lopinavir doses to boys than to girls should be evaluated in further studies. Such a study seems more especially important, as many countries do not practice therapeutic drug monitoring with antiretroviral drugs, which could balance the pharmacokinetic discrepancies between patients. However, the interest in therapeutic drug monitoring for children may be hampered by the high intraindividual variability that we observed. This variability is probably overestimated, as our study design (i.e., a single sample was obtained per occasion) did not allow us to estimate the interoccasion variability. This interoccasion variability could also be important, as the median delay between the first and the last occasion was 0.9 year (range, 0.03 to 2.4 years), a delay during which the pharmacokinetic parameters values naturally evolve, especially in young children.
In conclusion, this study has shown that the pharmacokinetics of lopinavir vary greatly during childhood and identified children younger than age 6 months and boys older than age 12 years as subpopulations in which the levels of exposure to lopinavir may be decreased. The virological and clinical relevance of this pharmacokinetic modification and the usefulness of a higher dosage in these children should be further evaluated.
Published ahead of print on 28 August 2006. ![]()
|
|
|---|
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»