Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, April 2004, p. 1159-1167, Vol. 48, No. 4
0066-4804/04/$08.00+0 DOI: 10.1128/AAC.48.4.1159-1167.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Department of Pharmacy,1 Division of Neonatology, Kitasato University Hospital,3 Division of Infectious Disease, Kitasato University School of Medicine, 1-15-1 Kitasato, Sagamihara-shi, Kanagawa 228-8555,2 Laboratory of Analytical Chemistry, School of Parmaceutical Sciences, Kitasato University, 108-8641 Shirogane, Minatoku, Tokyo, Japan4
Received 27 January 2003/ Returned for modification 29 August 2003/ Accepted 20 December 2003
|
|
|---|
33 weeks, Varbekacin = 0.54 liters/kg, CLvancomycin = 0.0250 x BW/serum creatinine level for PCAs of <34 weeks and CLvancomycin = 0.0323 x BW/serum creatinine level for PCAs of
34 weeks, Vvancomycin = 0.66 liters/kg, CLpanipenem = 0.0832 for PCAs of <33 weeks and CLpanipenem = 0.179 x BW for PCAs of
33 weeks, and Vpanipenem = 0.53 liters/kg. When the CL of each drug was evaluated by the nonlinear mixed-effect model, we found that the mean CL for subjects with PCAs of <33 to 34 weeks was significantly smaller than those with PCAs of
33 to 34 weeks, and CL showed an exponential increase with PCA. Many antibiotics are excreted by glomerular filtration, and maturation of glomerular filtration is the most important factor for estimation of antibiotic clearance. Clinicians should consider PCA, serum creatinine level, BW, and chemical features when determining the initial antibiotic dosing regimen for neonates. |
|
|---|
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Patient demographic dataa
|
Blood samples for retrieval of serum were obtained by capillary heel stick just prior to infusion. Peaks concentrations were measured 30 min or 1 h (2 to 3 h for vancomycin) after the end of infusion, and concentrations for evaluation of the elimination phase were measured 6 h after the end of infusion or just before administration of the next dose by the method of Sawchuk-Zaske (41). In our preliminary studies, trough levels were not detectable in some patients, and therefore, sampling every 6 h was also selected for measurement of the values for the elimination phase. Samples were obtained on days 1, 3, and 6 of the course of antibiotic treatment.
Serum was prepared by centrifugation and stored at -20°C prior to data analysis. Serum arbekacin concentrations were determined by high-performance liquid chromatography with a combination of reverse-phase, ion-pair chromatography, postcolumn derivatization with o-phthalaldehyde, and fluorescence detection, as described by Kubo et al. (28), with the following equipment and settings: detector, S-FL-330 (Soma Optics, Tokyo, Japan); wavelength, 440 nm; and column, Radial-PAK C18 (Nihon Waters, Tokyo, Japan). The lower limit of quantification was 0.3 mg/liter. The coefficient of variation was less than 3.0%. The recovery of arbekacin was 98.5%. Arbekacin was a gift from Meiji Seika. Serum vancomycin concentrations were analyzed by a fluorescence polarization immunoassay with a TDx Vancomycin Dinapack kit (Abbott, Osaka, Japan). The coefficient of variation was less than 3.0%. Serum panipenem concentrations were determined by high-performance liquid chromatography with a Multisolvent Delivery System 600 (Nihon Waters) with the following equipment and settings: detector, UV-970 (Jusco, Tokyo, Japan); wavelength, 300 nm; column, ODS-2 (GL Science, Tokyo, Japan); mobile phase, 35% (vol/vol) methanol in 5 mM sodium phosphate buffer containing 5 mM sodium dodecyl sulfate (pH 5.8); and flow rate, 0.8 ml/min. The lower limit of quantification was 0.5 mg/liter. The coefficient of variation was less than 2.7%. The recovery of panipenem was 98.4%. Panipenem was a gift from Sankyo.
Pharmacokinetic analysis.
The simultaneous analysis of all concentration-time and patient physiologic data was performed by using the NONMEM program with first-order methods (double precision, version V, level 1.0), a computer program developed for PPK analysis on a PCAT computer running under a Windows system (5). The concentration-time courses of each antibiotic were described by use of a one-compartment model (in accordance with the available post-distribution phase data) with a short infusion and first-order elimination. Fixed- and random-effect parameters were estimated by use of the NONMEM program. The basic pharmacokinetic parameters of CL and the volume of distribution (V) corresponding to the proposed model were determined for each patient by using the PREDPP subroutines supplied with ADVAN1/TRANS3 from the NONMEM library. In the first phase of the analysis, a basic model with no covariates on CL or V was used. Additive, proportional, and log-linear error models were compared for determination of the interindividual variability in CL and V and for the residual variability in the drug concentration. Proportional interindividual variability models were invoked for each
as follows:
![]() |
![]() |
are the population parameters.
CL and
V are random variables with a mean of zero and a variance of
2. CLi and Vi are the individualized estimated parameters.
The additive error model was used for residual variability:
![]() |
ij is the independent identically distributed statistical error with a mean of zero and a variance of
2.
Fixed-effects modeling.
In the fitting process, the following patient demographic and biochemical data were used as covariables in the population model: PCA, GA, PNA, Apgar score, BW, and Cr. Covariates were investigated on the basis of the CL and V values for each antibiotic. The candidate covariate was screened in turn by adding it to the basic model: CL =
1 +
2 x covariate and V =
1 +
2 x covariate, where
1 and
2 are the intercept and slope parameters, respectively. An objective function value (the negative value of twice the log-likelihood difference [-2 l.l.d.]) is produced by the NONMEM program for each model in this regression. Comparisons among the different models were based on the differences in the minimum value of the objective function. Changes in the objective function greater than 6.635 indicate a statistically significant (P < 0.01) improvement in the fit of the data on the basis of a
2 distribution with 1 degree of freedom. On the basis of these results, all significant factors were used to construct the full regression equation. A stepwise procedure was used to determine the final model. The final model was obtained by removing covariates from the full model. After deletion of each factor in the full model, the objective function value of this reduced model was compared with that of the full model. At that time, a more restrictive criterion was adopted, and -2 l.l.d. of more than 7.8 was required to maintain covariates for the final model (P < 0.005). The final estimates of the PPK parameters were defined. Individual estimates of the CL and V values for each antibiotic were generated by Bayesian feedback with the NONMEM program after the population analysis.
Model validation. In the present study the bootstrap resampling technique was applied as an internal validation. The distribution of the pharmacokinetic parameters was confirmed by selecting 500 data sets by sampling and replacing data sets as needed. Finally, the mean values for the parameters from the bootstrap replicates were compared with the mean values made with the original data set.
To evaluate the predictive performance of the population parameters, we also obtained individual Bayesian estimates of CL and V of each antibiotic for these patients with the POSTHOC option in the NONMEM software package. The pharmacokinetic parameters derived in this study were evaluated by determination of their distribution against the linear regression of the predicted versus the observed concentrations in serum.
|
|
|---|
The mean CL and V values were as follows: CLarbekacin = 0.107 liters/h, Varbekacin = 0.91 liters, CLvancomycin = 0.105 liters/h, Vvancomycin = 0.91 liters, CLpanipenem = 0.175 liters/h, and Vpanipenem = 0.55 liters. Table 2 summarizes the main covariate models that were tested by use of the NONMEM program. In the fitting process, the NONMEM analysis with single covariates identified PCA, GA, PNA (except with vancomycin), BW, inverted Cr, exponential PCA, the complex factor of BW, and inverted Cr (BW/Cr) as having an influence on CL. The Apgar score was not significantly correlated. The values of -2 l.l.d. for the effects of different PCAs on CL were minimal for PCAs of 33 to 34 weeks. When evaluating each CL, we found that the mean CL for the subjects with PCAs of <33 to 34 weeks was significantly smaller than the mean CL for those with PCAs of
33 to 34 weeks, and for all three antibiotics CL showed an exponential increase with PCA. The influence of PCAs of <33 to 34 weeks, 33 to 34 weeks, or >33 to 34 weeks on the CL was chosen for the next analysis step.
|
View this table: [in a new window] |
TABLE 2. Hypothesis testing for influence of fixed effects on pharmacokinetic parameters for three types of antibiotics
|
1 +
2 x BW x (1 +
3/Cr) (PCAs, <33 or 34 weeks), CL =
4 +
5 x BW x (1 +
6/Cr) (PCAs,
33 or 34 weeks), and V =
7+
8 x BW, where
3 to
8 are the intercept or slope parameters.
Each parameter was eliminated from the full model (Table 3). The final model with PCA, Cr, and BW was examined; and the final regression equations and parameter estimates are shown in Table 4. The final formulas for the PPK parameters were as follows: CLarbekacin = 0.0238 x BW/Cr for PCAs of <33 weeks and CLarbekacin = 0.0367 x BW/Cr for PCAs of
33 weeks, Varbekacin = 0.54 liters/kg, CLvancomycin = 0.0250 x BW/Cr for PCAs of <34 weeks and CLvancomycin = 0.0323 x BW/Cr for PCAs of
34 weeks, Vvancomycin = 0.66 liters/kg, CLpanipenem = 0.0832 for PCAs of <33 weeks and CLpanipenem = 0.179 x BW for PCAs of
33 weeks, and Vpanipenem = 0.53 liters/kg.
|
View this table: [in a new window] |
TABLE 3. Hypothesis testing according to intersubject variability by using reduced models of the full model
|
|
View this table: [in a new window] |
TABLE 4. Final estimates for the population pharmacokinetic parameters for antibiotics in neonates
|
The standard deviations in intraindividual variability were determined to be 2.43, 3.22, and 4.2 mg/liter for arebkacin, vancomycin, and panipenem, respectively. The highly significant relationships between pharmacokinetic parameters and physiologic factors are presented in Fig. 1 to 4. These parameters were estimated by Bayesian feedback with the NONMEM program.
![]() View larger version (18K): [in a new window] |
FIG. 1. Relationship between CLarbekacin normalized by weight and PCA. CLarbekacin was estimated by the Bayesian method. The horizontal lines show the average CLarbekacin for neonates with PCAs of <33 and 33 weeks. CLarbekacin was normalized by weight but depends on PCA and Cr (Scr).
|
![]() View larger version (23K): [in a new window] |
FIG. 4. Scatterplot of Varbekacin, Vvancomycin, and Vpanipenem versus PCA. Each value of V was estimated by the Bayesian method. There was no significant relationship between V and time.
|
Figure 5 shows a scatterplot of the observed concentrations and the estimated concentrations. The mean values for the population analyzed showed a good predictive performance. The distribution of the plot against the line for the formula for CL displays a bilateral symmetry around the regression, so the observed population mean was assumed to be an unprejudiced parameter.
![]() View larger version (25K): [in a new window] |
FIG. 5. Scatterplot of observed concentrations versus estimated concentrations of three kinds of antibiotics. Estimated concentrations were calculated by using our population model, according to the covariate results for each individual from NONMEM analysis.
|
|
|
|---|
In the present study, PCA, GA, PNA, Cr, and BW were important factors that correlated with individual estimates of CL (Table 2); but the correlation for PNA was not clearly significant. Many investigators have reported a high correlation between GA and antibiotic CL (7, 8, 23, 45, 46). Clinically, renal function is the most important factor for elimination because many antibiotics are completely excreted in urine. Many studies have evaluated the changes in renal function during development in neonates according to their GAs. Although PNA must be considered along with GA to evaluate renal function correctly, it is difficult to identify the significant covariate in the process of determining population parameters (17). The use of a large number of patients of the same GA with various PNAs is necessary for analysis of a small population. For this reason, the use of PCA is an obvious choice for the evaluation of renal function in neonates. Our results showed that the mean CL of antibiotics in patients with PCAs of <33 to 34 weeks was significantly lower than that in patients with PCAs of
33 to 34 weeks, and CL showed a logarithmic rise with PCA. Many investigators (2, 3, 18, 29) have reported that changes in the development of the renal function in neonates are correlated with PCA and that renal function increases after 34 to 35 weeks. These associations suggest that antibiotic CL steeply increases in neonates with PCAs of 33 to 36 weeks.
If patients with a wider range of renal function had been available, it is likely that creatinine CL would have been a significant predictor of CL for antibiotics excreted mainly by glomerular filtration. Schwartz et al. (42) estimate GFR (in millimeters per minute per 1.73 m2) as follows: GFR = k · HT/Cr, where HT is the height (in centimeters) and k is a constant proportional to HT/Cr. Although HT or body surface area might be more accurate for estimation of GFR, many dosages have previously been normalized per kilogram of BW. After the inclusion of Cr and BW in the regression model, BW/Cr was a powerful predictor of the CL of each antibiotic in the early stages of the model-building process, but for panipenem it was not significant in the final model.
In this study, significant covariates related to CLpanipenem were found, and it was confirmed that CLpanipenem does not depend on the maturation process. Panipenem and other carbapenems are metabolized by the enzyme dehydropeptidase-1 and are very unstable in solution (34). The renal function of low-birth-weight babies is extremely low; therefore, CL by extrarenal pathways is relatively greater than that by the renal pathway according to the level of maturation of the neonate. Our results indicate that CLpanipenem (0.179 liter/h) was almost identical to GFR in full-term babies (0.144 liter/h; GA, 38 to 41 weeks) (24), and similar CLs were observed for other carbapenems, such as meropenem (0.157 liter/h/kg) (46) and impenem (0.150 liter/h/kg) (39).
The renal CL of antibiotics is affected by glomerular filtration, tubular secretion, and tubular reabsorption (35). Antibiotics with high degrees of binding to serum proteins show a markedly lower renal CL (6). We reviewed 14 reports in the literature about the relationship between plasma protein binding and renal CL for seven antibiotics in neonates, and the results are shown in Fig. 6. With some exceptions, protein binding ratios provide good theoretical values of renal CL (renal CL = fu x GFR) in neonates (fu is the unbound fraction of antibiotic). A similar result was shown in the present study. Because of the high degree of interpatient variability, very small CL, and low body protein levels, CLvancomycin and CLarbekacin were similar in neonates with PCAs of less than 33 weeks. However, in the group with PCAs of >33 weeks, CLvancomycin was lower than CLarbekacin. The levels of protein binding of vancomycin and arbekacin are less than 10 to 82% (1) and less than 15% (31), respectively.
![]() View larger version (15K): [in a new window] |
FIG. 6. Correlation between total and renal CL of antibiotics with protein binding ratios reported in the literature. Fourteen antibiotics were reviewed in 16 studies. The values of antibiotic CL from the literature were divided into two groups: circles indicate PCA of 33 weeks or BW of 2,000 g for amikacin (36), amoxicillin (11, 23), aztreonam (13), cefoperazone (38), cefotaxime (21), ceftazidime (40, 44), ceftizoxime (26), ceftriaxone (32), meropenem (46), and piperacillin (25); squares indicate PCA of <33 weeks or BW of <2,000 for cefoperazone (8), flucloxacillin (10), piperacillin (25), ticarcillin (11), and vancomycin (12). The figure shows that CL is associated with protein binding ratios and that the relationship between renal CL and protein binding is stronger than that between total CL and protein binding.
|
![]() View larger version (21K): [in a new window] |
FIG. 7. Correlation between PCA and CL of arbekacin, vancomycin, and panipenem. Antibiotic CL was estimated by the Bayesian method.
|
In our study the mean Vvancomycin for the neonates was smaller than that for adults, and conversely, the mean Varbekacin and Vpanipenem for the neonates were larger than those for adults. Those findings confirm the findings of other investigators with the same group of antibiotics (4, 17, 27, 49). We suggest that the discrepancy may be due to estimation of the distribution property in adults. Although some antibiotics are extensively bound to tissue and plasma proteins and have very large V in adults, many antibiotics show small apparent V (0.1 to 0.3 liters/kg) due to high levels of aqueous solubility. Antibiotic diffusion from blood to extracellular fluid depends on many factors, such as tissue binding, free drug concentration (plasma protein binding), pKa, molecular weight, and lipid solubility (37, 47). Plasma protein binding may interfere with the extravascular distribution (9). The distribution profiles of 12 antibiotics with a wide range of protein binding values in neonates and the relationship between V in neonates and V in adults are reviewed and presented in Fig. 8. V in neonates is difficult to estimate from values for adults (Fig. 8), but it depends on the level of protein binding and the amount of extracellular fluid. While protein binding is an important factor for estimation of pharmacokinetics parameters, we need to consider the significant differences between neonates and adults. In our study Varbekacin and Vpanipenem were double those for adults because of low plasma protein binding levels. Because Vvancomycin is large (more than 0.6 liters/kg) in adults, the drug is likely extensively distributed in tissue, in addition to extracellular fluid. The apparent Vvancomycin in neonates tends to be small compared to that in adults, even though the antibiotic is hydrophilic; the reason may be that increased BW is the result of increased extracellular fluid levels in neonates and, thus, Vvancomycin is relatively small because vancomycin does not distribute to extracellular fluid. These data on hydrophilic antibiotics assume that the value of V in neonates, which is less than 0.4 liters/kg in adults, depends on the protein binding ratio. For antibiotics for which the values of V are >0.6 liters/kg in adults, the values of V are smaller in neonates.
![]() View larger version (15K): [in a new window] |
FIG. 8. (A) Correlation between apparent V for antibiotics in adults and neonates reported in the literature. Apparent V for amikacin (36), amoxicillin (11, 23), cefazolin (14), cefoperazone (8, 38), cefotaxime (21), ceftazidime (40, 44), ceftizoxime (26), ceftriaxone (30, 32), flucloxacillin (22), meropenem (46), piperacillin (25), and ticarcillin (10) are low (less than 0.4 liters/kg) in adults due to aqueous solubility. (B) Correlation between protein binding ratio and V. The apparent V was strongly associated with the protein binding ratio. As the protein binding ratio increases in serum (fu decreases), the apparent V is decreased.
|
![]() View larger version (17K): [in a new window] |
FIG. 2. Relationship between CLvancomycin normalized by weight and PCA. CLvancomycin was estimated by the Bayesian method. The horizontal lines shows the average CLvancomycin for neonates with PCAs of <34 and 34 weeks. CLvancomycin normalized by weight still depends on PCA and Cr (Scr).
|
![]() View larger version (17K): [in a new window] |
FIG. 3. Relationship between CLpanipenem and PCA. CLpanipenem for neonates with PCAs of <33 weeks was not dependent on maturation but was constant for extrarenal CL; CLpanipenem for neonates with PCAs of 34 weeks increased by PCA.
|
|
|
|---|
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»