Skip to main page content

• HOME
• Current Issue
• Archive
• Alerts
• About ASM
• Contact us
• Tech Support
• Journals.ASM.org

This Article Abstract Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Citation Map Services E-mail this article to a friend Similar articles in this journal Similar articles in ASM journals Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Reprints and Permissions Copyright Information Books from ASM Press MicrobeWorld Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Vinks, A. A. Articles by Mouton, J. W. Search for Related Content PubMed PubMed Citation Articles by Vinks, A. A. Articles by Mouton, J. W.

 Previous Article  |  Next Article 

Antimicrobial Agents and Chemotherapy, February 2003, p. 541-547, Vol. 47, No. 2
0066-4804/03/$08.00+0     DOI: 10.1128/AAC.47.2.541-547.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.

Population Pharmacokinetic Analysis of Nonlinear Behavior of Piperacillin during Intermittent or Continuous Infusion in Patients with Cystic Fibrosis

Alexander A. Vinks,^1* Jan G. den Hollander,² Shelley E. Overbeek,³ Roger W. Jelliffe,⁴ and Johan W. Mouton^2,5

Division of Pharmacology Research and Clinical Pharmacology, Cincinnati Children's Hospital Medical Center, Cincinnati, Ohio,¹ Department of Medical Microbiology and Infectious Diseases, Canisius-Wilhelmina Hospital, Nijmegen,⁵ Department of Clinical Microbiology and Infectious Diseases,² Department of Pulmonary Medicine, Erasmus University Medical Center, Rotterdam, The Netherlands,³ Laboratory of Applied Pharmacokinetics, Keck School of Medicine, University of Southern California, Los Angeles, California⁴

Received 25 October 2001/ Returned for modification 12 August 2002/ Accepted 31 October 2002

Top Abstract Introduction Materials and Methods Results Discussion References

 
The purpose of this study was to describe the nonlinear pharmacokinetics of piperacillin observed during intermittent infusion and continuous infusion by using a nonparametric population modeling approach. Data were 120 serum piperacillin concentration measurements from eight adult cystic fibrosis (CF) patients. Individual pharmacokinetic parameter estimates during intermittent infusion or continuous infusion were calculated by noncompartmental analysis and with a maximum iterative two-stage Bayesian estimator. To simultaneously describe concentration-time data during intermittent infusion and continuous infusion, nonlinear models were parameterized as two-compartment Michaelis-Menten models. Models were fit to the data with the nonparametric expectation maximization algorithm. The calculations were executed on a remote supercomputer. Nonlinear models were evaluated by log-likelihood estimates, residual plots, and R² values, and predictive performance was based on bias (mean weighted error [MWE]) and precision (mean weighted square error [MWSE]). A linear pharmacokinetic model could not describe combined intermittent infusion and continuous infusion data well. A good population model fit to the intermittent infusion and continuous infusion data was obtained with the constructed nonlinear models. Maximum a posteriori probability (MAP) Bayesian R² values for the nonlinear models were 0.96 to 0.97. Median parameter estimates for the best nonlinear model were as follows: K_m, 58 ± 75 mg/liter (mean and standard deviation); V_max, 1,904 ± 1,009 mg/h; volume of distribution of the central compartment, 14.1 ± 3.0 liters; k₁₂, 0.63 ± 0.41 h^-1; and k₂₁, 0.37 ± 0.19 h^-1. The median bias (MWE) and precision (MWSE) values for MAP Bayesian estimation with the Michaelis-Menten model were 0.05 and 4.6 mg/liters, respectively. The developed nonlinear pharmacokinetic models can be used to optimize piperacillin therapy administered via continuous infusion in patients with CF and have distinct advantages over conventional linear models.

Top Abstract Introduction Materials and Methods Results Discussion References

 
Piperacillin-tazobactam is a parenterally administered combination of a ß-lactam antibiotic with a ß-lactamase inhibitor (4). This drug combination is frequently used in the treatment of pulmonary exacerbations in patients with cystic fibrosis (CF) because of its good bactericidal activity against Pseudomonas aeruginosa and gram-positive microorganisms (20). Piperacillin, like other ß-lactam antibiotics, exhibits time-dependent killing, and the parameter for the bactericidal effect is the time during which the concentration of the drug exceeds the MIC (5). As shown for ceftazidime (23), continuous infusion of piperacillin-tazobactam may result in superior efficacy with equal or lower total daily doses. Administration by continuous infusion also has distinct advantages in CF home treatment programs (24). A recent report on continuous intravenous administration of vancomycin also emphasized this point (27).

Data on the pharmacokinetics of piperacillin in patients with CF are scarce (12; J. S. Bertino, M. D. Reed, C. Meyers, and J. L. Blumer, ediatr. Res. 16:1210A, abstr. 254, 1982), while data for the combination are lacking. Piperacillin-tazobactam is commonly administered by intermittent infusion. In non-CF patient populations and with this mode of administration, the pharmacokinetics of the ß-lactamase inhibitor- antibiotic combination appeared to be linear (1), although nonlinear analysis was not addressed. Some authors, however, have argued for nonlinear behavior, and studies have documented evidence for decreasing clearance with increasing concentrations (2, 3, 19). Piperacillin is, for the most part, excreted through the kidneys by active tubular secretion and by glomerular filtration. A typical characteristic of the active secretion process is a limitation in the capacity of this process. This behavior can be described in terms of Michaelis-Menten (MM) kinetics, where the rate of active transport increases toward a maximum value (V_max) with increasing plasma drug concentrations (11). The concentration at which the rate is 50% of its maximum is the K_m. At concentrations well below the K_m, the rate of excretion and the concentration vary in direct proportion. Drug clearance decreases in a nonlinear fashion with increasing concentrations, most notably at concentrations exceeding the K_m. In patients with CF, drugs that are extensively excreted by renal tubules (such as penicillins) show enhanced renal clearance (20). Using multiple constant-rate infusions of ticarcillin, Wang et al. (25) showed that renal clearance was enhanced due to increased affinity of the renal secretory system. Typically, secretory capacity was unaltered, but K_m values were significantly lower than those of controls.

Recently, while studying the pharmacokinetics and efficacy of high-dose piperacillin-tazobactam delivered via continuous infusion in patients with CF, den Hollander et al. observed significantly higher piperacillin clearance during continuous infusion than during intermittent infusion (J. G. den Hollander, S. E. Overbeek, A. A. Vinks, H. A. Verbrugh, and J. W. Mouton, Abstr. Proc. 20th Int. Congr. Chemother., abstr. 3071, 1997). This observation led us to hypothesize that the observed differences may be explained by a mechanism similar to that described for ticarcillin. As nonlinearity is a concentration-dependent phenomenon, it will be influenced by the mode of administration, especially when dosing regimens lead to very different concentrations. The high concentrations in serum during intermittent infusion may saturate the tubular secretion process, while the relatively low concentrations during continuous infusion (which remain below the K_m) may explain the observed differences in clearance. This scenario may have important clinical implications for the dosing of piperacillin administered via continuous infusion.

Until recently, it was not possible to analyze such nonlinear data, as appropriate software was lacking. With the development of a resource for the modeling of nonlinear pharmacokinetic-pharmacodynamic systems via the Internet, such an analysis has become feasible. The purpose of the present investigation was to analyze the nonlinear pharmacokinetics of piperacillin during intermittent infusion and continuous infusion in patients with CF by using a nonparametric population modeling approach. As tazobactam concentrations did not have nonlinear characteristics, these data were not part of the study.

Top Abstract Introduction Materials and Methods Results Discussion References

 
Patients. Data were 120 serum concentration measurements obtained from eight adult CF patients (six women and two men) during intermittent infusion and continuous infusion of a piperacillin-tazobactam combination (Tazocin; AHP, Hoofddorp, The Netherlands). Patients participated in a clinical trial that was conducted at University Hospital Dijkzigt, Rotterdam, The Netherlands, and that was designed to evaluate the clinical efficacy, safety, and pharmacokinetics of piperacillin-tazobactam in CF patients during intermittent infusion and continuous infusion. Patients were treated for an acute exacerbation of their chronic P. aeruginosa infection with high doses of piperacillin-tazobactam. The total daily doses were 300 mg of piperacillin/kg of body weight/day (16 g) and 37.5 mg of tazobactam/kg/day (2 g). Patients received the piperacillin-tazobactam combination (piperacillin, 4,000 mg; tazobactam, 500 mg) in a crossover fashion either via intermittent infusion (piperacillin, 4 g; tazobactam, 0.5 g) every 6 h or via continuous infusion (piperacillin, 16 g; tazobactam, 2 g) over 24 h with an infusion pump. Prior to treatment, patients were randomized to either one of the regimens and were switched to the alternative regimen after 48 h.

Venous blood samples for drug assays were obtained from the arm opposite that used for the infusion of the antibiotic by using an indwelling needle. During intermittent infusion, samples were taken on day 2 prior to the sixth dose (time zero) and at 5, 30, 45, and 60 min and 1.5, 2, 3, 4, 5, and 6 h after the start of infusion. The duration of intermittent infusion was 30 min. During continuous infusion, samples were taken on the second day of continuous infusion treatment, just prior to the changing of the infusion reservoir (time zero), and at 2, 4, and 6 h into the infusion regimen. After samples were obtained, the blood was allowed to clot on ice for 20 min. It was then centrifuged, and the serum was stored at -70°C until analysis.

The study was approved by the Institutional Review Board of the University Hospital Dijkzigt, and written informed consent was obtained from each patient.

Drug assay. Piperacillin concentrations in serum were analyzed by high-performance liquid chromatography by a modification of the method of Reed et al. (14). For piperacillin, briefly, 0.2 ml of a methanol solution containing 50 mg of methicillin/liter as an internal standard was added to an equal volume of serum sample, and the mixture was centrifuged. The supernatant was injected into a C₁₈ reversed-phase column (Chrompack, Middelburg, The Netherlands) with 0.05 M KH₂PO₄-0.05 M K₂HPO₄ (pH 6.9) containing 16% (vol/vol) acetonitrile as the mobile phase. Detection was done at 215 nm. The method was validated by using published guidelines (22). The lower limit of quantitation was 0.5 mg/liter and was determined by the method of Kucharczyk (10). All high-performance liquid chromatography measurements were obtained at 4°C.

In order to weight each serum drug concentration correctly by the reciprocal of its variance (Fisher information), the assay standard deviation (SD) was determined over the working range of the assay as described by Jelliffe et al. (9). The assay variance was estimated by regression modeling on the basis of four different concentrations over the working range. The assay error pattern was fitted best by the polynomial equation SD = 0.71470 + (0.030110 x C) + (0.000245 x C²), where the SD for the assay is given in milligrams per liter, C is the measured piperacillin concentration, and C² is the square of C. This polynomial equation for the assay error was also incorporated into the appropriate sections of the iterative two-stage Bayesian (IT2B) program and the nonparametric expectation maximization algorithm (NPEM3) program for population analysis (User Manual for Version 10.7 of the USC_*PACK Collection of PC Programs, Laboratory of Applied Pharmacokinetics, School of Medicine, University of Southern California, Los Angeles).

Pharmacokinetic parameter estimates. Serum piperacillin concentration-time data (intermittent infusion and continuous infusion data) were analyzed by noncompartmental analysis (KINFIT, MW/PHARM, version 3.50; MediWare, Groningen, The Netherlands) (13) and by IT2B estimation (front part of NPEM3; intermittent infusion data) (User Manual for Version 10.7 of the USC_*PACK Collection of PC Programs). IT2B modeling included one- and two-compartment open-model analyses with zero-order input and first-order elimination from the central compartment. Parameters in the model were volume of distribution of the central compartment (V_c) (in liters), the elimination rate constant k_el (in hour^-1), and the intercompartmental rate constants k₁₂ and k₂₁ (in hour^-1), if applicable. The inverse of the estimated assay variance was used as the weighting in the pharmacokinetic modeling throughout. The goodness of fit was assessed by visual inspection of the predicted concentration-time profiles, log-likelihood estimates, residual analyses, and coefficients of determination. Model discrimination was performed by using Akaike's information criterion (AIC) (28). The Wilcoxon signed rank test was used to compare individual pharmacokinetic parameter estimates during intermittent infusion and continuous infusion. A P value of <0.05 was considered significant.

The IT2B individual pharmacokinetic parameter estimates were used in the simulation module of MW/PHARM to obtain piperacillin concentration estimates during continuous infusion.

Population pharmacokinetic analysis. Based on the IT2B analysis and to mimic the renal elimination pathways for piperacillin, tubular secretion and glomerular filtration (25), models were parameterized as two-compartment models with elimination from the central compartment modeled as an MM process, as an MM process in combination with first-order elimination, or as a first-order elimination process only (Fig. 1).

View larger version (8K): [in this window] [in a new window]

FIG. 1. Two-compartment open model used to fit piperacillin data. IV, zero-order infusion rate; Vc, central compartment; Vp, peripheral compartment; Ke, first-order rate constant; Km, MM constant; Vmax, maximum excretion rate; k12 and k21, intercompartmental rate constants. Structural models evaluated in this study had either first-order elimination (Ke), MM elimination, or a combination of MM and first-order elimination.

The general differential equation for the models was as follows:

where K_m is the MM constant (in milligrams per liter), V_max is the maximum elimination rate (in milligrams per hour), C_p is the plasma drug concentration (in milligrams per liter), C₂ is the concentration in the peripheral compartment, k_el is the first-order elimination rate constant (in hour^-1), k₁₂ and k₂₁ are distribution rate constants, and r(t) is the zero-order drug input rate in the central compartment (in milligrams per hour). In one model, the elimination rate constant k_e was defined as a function of the creatinine clearance (CL_CR) by the equation k_el = K_i + (K_s x CL_CR), where K_i is the intercept (representing residual nonrenal elimination) and K_s is the slope factor.

Differential equations for the models were defined in BOXES (USC_*PACK PC programs, version 10.7) (User Manual for Version 10.7 of the USC_*PACK Collection of PC Programs). Four different models were used in the final analysis (Fig. 1). (i) The MM model includes no first-order elimination from the central compartment. (ii) The MM-k_el combination model includes first-order elimination (k_el) from the central compartment, defined as k_el = 0.01 + (K_s x CL_CR), with K_s fixed at 0.0052 h^-1, according to published data (18). The value chosen for K_i corresponds to the residual k_el in dialysis patients (CL_CR, 0; half-life, 70 h). K_s is the mean slope value, as reported for people with normal renal function (CL_CR, 120 ml/min) (18). This model was evaluated to include CL_CR as a potential useful clinical descriptor. (ii) The MM-Ke combination model includes first-order elimination (k_el) from the central compartment, with no descriptors for k_el. (iv) The linear k_el model includes no MM elimination.

All models include V_c, k₁₂, and k₂₁ as additional parameters. In the modeling, the ranges for K_m and V_max were set at 0 to 750 mg/liter and 0 to 4,000 mg/h, respectively. Models were fit to the combined intermittent infusion and continuous infusion data with the Big-NPEM program (NPEM3) developed by R. W. Jelliffe and A. Schumitzky (16). To enable comparison with the IT2B analysis of intermittent infusion data, the best nonlinear model (MM model) was also fit to the intermittent infusion data. NPEM3 resides on the Cray T3E supercomputer at the San Diego Supercomputer Center and can be accessed over the Internet through the file transmission protocol. This resource can also be accessed more easily from its own web site (https://gridport.npaci.edu/LAPK/). Similar but smaller parallel computing facilities can now also be accessed at the University of Southern California Laboratory of Applied Pharmacokinetics modeling resource at www.lapk.org. Calculations were executed on the remote mainframe computer by up- and downloading of data and results via the Big-NPEM protocol and through the file transmission protocol. The inverse of the estimated assay variance was used as the weighting. With NPEM3, the total observation variance is approximated by the assay variance. Models were evaluated by residual analyses, coefficients of determination, AIC, and log-likelihood values. Finally, maximum a posteriori probability (MAP) Bayesian estimation was used to generate parameter estimates for each patient. Predictive performance evaluation was based on bias (mean weighted error [MWE]) and precision (mean weighted square error [MWSE]) (17; User Manual for Version 10.7 of the USC_*PACK Collection of PC Programs).

Top Abstract Introduction Materials and Methods Results Discussion References

 
Pharmacokinetic parameter estimates. The demographic data and genotypes for the patients are given in Tables 1 and 2, respectively.

View this table: [in this window] [in a new window]

TABLE 1. Demographic and clinical data for eight adult patients with CF and participating in this study

View this table: [in this window] [in a new window]

TABLE 2. Genotypes for patients in this study

In the IT2B modeling, a two-compartment model best described the data for most of the patients (six of eight), based on the AIC, and was therefore chosen as the most appropriate for the population analysis. Mean piperacillin pharmacokinetic values based on separate intermittent infusion and continuous infusion noncompartmental data analyses are summarized in Table 3. Piperacillin clearance during continuous infusion (24.4 ± 11.7 liters/h) (mean and SD) was significantly higher than that during intermittent infusion (13.1 ± 2.3 liters/h) (P < 0.02). The mean and median parameter estimates for the IT2B modeling with intermittent infusion data only are shown in Table 4. These estimates were used to predict piperacillin concentrations during continuous infusion. However, as shown in Fig. 2, the observed piperacillin concentrations during continuous infusion (32.8 ± 14.4 mg/liter) (mean and SD) were significantly lower than those predicted (54.3 ± 12.1 mg/liter) (P < 0.001) by the individual IT2B pharmacokinetic values.

View this table: [in this window] [in a new window]

TABLE 3. Pharmacokinetic parameter estimates obtained by noncompartmental analysisa

View this table: [in this window] [in a new window]

TABLE 4. Piperacillin IT2B pharmacokinetic parameter estimatesa obtained with a linear pharmacokinetic model

View larger version (22K): [in this window] [in a new window]

FIG. 2. Mean (and SD) piperacillin concentration-time profile during intermittent infusion (piperacillin, 4 g; tazobactam, 0.5 g; every 6 h) and continuous infusion (piperacillin, 16 g; tazobactam, 2 g; over 24 h) with the linear model. The solid line represents the fit with the linear model based on the intermittent pharmacokinetic data. For clarity, drug regimen crossover data are presented similarly (intermittent infusion followed by continuous infusion).

Population modeling. The median and mean population parameter estimates and associated dispersion factors generated by NPEM3 are shown in Table 5. Parameter values were not significantly different among nonlinear models. In contrast, mean and median estimates for volume of distribution in the linear model (e.g., median V_c, 23.1 ± 13.5 liters) were significantly higher than the nonlinear model estimates (e.g., median V_c, 13.7 ± 3.3 liters for the MM model) (P < 0.04).

View this table: [in this window] [in a new window]

TABLE 5. Piperacillin population pharmacokinetic parameter estimates obtained by NPEM3 analysisa

Evaluation of predictive performance. Table 6 summarizes log-likelihood values, coefficients of determination (R²), and predictive performance of the three nonlinear models and the linear model, determined with mean error (measure of bias) and mean square error (measure of precision). Based on the log-likelihood values, the MM model describes the data best. The bias and precision estimates of the nonlinear models were in the same range, but all were significantly better (P < 0.01) than those of the linear model, as is evident from the smaller MWE and MWSE values.

View this table: [in this window] [in a new window]

TABLE 6. Evaluation of the predictive performance of piperacillin population modelsa

The MM model based on intermittent infusion data fit the data better than the linear model (Fig. 3). This fit was almost identical to that obtained with the MM model based on combined intermittent infusion and continuous infusion data (Tables 5 and 6). Figure 4 shows the relationship between the observed piperacillin concentrations (y axis) and the MAP Bayesian-predicted concentrations obtained with the linear and MM models (x axis) for all 120 observations in the eight patients. In the Bayesian estimation, median values were used for each parameter as the measure of central tendency. This strategy has been found to improve predictive performance (R. W. Jelliffe, personal communication). As determined by weighted linear regression analysis, the line of best fit (solid line) for the linear k_el model differed significantly from the line of identity (broken line) (Fig. 4A). The line of best fit (solid line) for the MM model did not differ from the line of identity (broken line), and the intercept did not differ significantly from 0.0 (Fig. 4A). The coefficient of determination (R²) was 0.973. Predictive performance relationships for the other nonlinear combination models were comparable to those for the MM model (R², 0.956 and 0.960).

View larger version (22K): [in this window] [in a new window]

FIG. 3. Mean (and SD) piperacillin concentration-time profile during intermittent infusion (piperacillin, 4 g; tazobactam, 0.5 g; every 6 h) and continuous infusion (piperacillin, 16 g; tazobactam, 2 g; over 24 h) with the MM model. The solid line represents the fit with the MM model based on the intermittent pharmacokinetic data.

View larger version (21K): [in this window] [in a new window]

FIG. 4. Scatter plot of observed versus MAP Bayesian-predicted concentrations obtained with the population model median parameter estimates for the linear model and the best nonlinear model (MM model). (A) For the kel model, the slope and the intercept of the line of best fit (solid line) are significantly different from 1.0 and 0.0 (broken line), respectively. The R2 value is 0.88. (B) For the MM model, the slope and the intercept of the line of best fit (solid line) are not different from identity and 0.0 (broken line), respectively. The R2 value is 0.97. pred. conc., predicted concentration; post., a posteriori.

Top Abstract Introduction Materials and Methods Results Discussion References

 
The purpose of this study was to describe the nonlinear behavior of piperacillin observed when the drug was administered via two distinct dosing methods, intermittent infusion and continuous infusion (Den Hollander et al., Abstr. Proc. 20th Int. Congr. Chemother.). During continuous infusion of the same total daily dose, piperacillin steady-state concentrations were, on average, 40% (range, 6 to 69%) lower than expected based on intermittent infusion pharmacokinetic parameter estimates. The results indicate that a nonlinear model better describes piperacillin pharmacokinetics in patients with CF than does a linear model (Fig. 4 and Table 6). With the increasing use of ß-lactam antibiotics administered by continuous infusion, this finding has important clinical implications. Most dosing guidelines for continuous infusion of antibiotics are not based on clinical efficacy studies but are typically derived from the total daily intermittent infusion dose (24). As the renal clearance of ticarcillin, a related acylureido-penicillin derivative, was increased in patients with CF as a result of increased tubular secretion compared to that in healthy volunteers (25), we hypothesized the existence of a similarly affected renal anion excretion system for piperacillin. Recently, more theoretical evidence for this phenomenon was provided. Investigators studying CF were able to simulate an increased rate of renal organic anion secretion in the renal tubules and hence, due to enhanced tubular secretion, lower steady-state serum drug concentrations (26). An enhanced affinity for the secretory site(s) similar to that shown for ticarcillin in patients with CF and lower K_m values may very well explain our observations for piperacillin. In fact, with our nonlinear modeling approach, the K_m and V_max estimates obtained in the present study correspond well to those reported for ticarcillin. In the study of Wang et al. (25), unbound ticarcillin K_m and V_max estimates were 33.7 mg/liter (coefficient of variation, 36%) and 950 mg/h (coefficient of variation, 48%), respectively. However, when corrected for the concentration-dependent protein binding (37 to 57%) of ticarcillin, the corresponding K_m and V_max estimates for total ticarcillin are in the ranges of 91 mg/liter and 1,667 to 2,568 mg/h, respectively, and are in good agreement with our data for piperacillin (Table 5).

Over the past decades, there has been controversy over the proposed nonlinear elimination of piperacillin. Interestingly, for healthy volunteers, a recent population pharmacokinetic analysis claimed linear pharmacokinetics for piperacillin after multiple doses (1). An explanation for the inconsistent findings may be related to study design. In the present study, we used the same dose administered via two different modes of administration, intermittent infusion and continuous infusion, in a crossover fashion. Pharmacokinetic modeling of piperacillin after intermittent administration did not show obvious nonlinear or "hockey stick" behavior indicative of nonlinear elimination. As the elimination of piperacillin is, for the most part, governed by two processes, nonlinear secretion and linear filtration, changes in tubular secretion at the point at which concentrations approach saturation may very well be compensated for by increases in the filtration clearance of unbound drug. Only at concentrations well below the K_m may these effects on clearance become most obvious.

Although the objective of this study was to model the nonlinear behavior of piperacillin, we also estimated clearance and other pharmacokinetic parameters after multiple doses by linear modeling. The calculated pharmacokinetic parameter values after intermittent infusion are different from earlier observations in young patients with CF. Reed et al. (15) and Bertino et al. (Pediatr. Res. 16:121A, abstr. 254), using high-dose therapy (600 to 900 mg/kg/day), found mean total body clearance values of 305.6 and 341.6 ml/min/1.73 m², respectively, values that are higher than the 226.2 ml/min/1.73 m² found in this study. The main reason for this discrepancy is most likely the older age of our patients, as with age, renal elimination gradually diminishes. Besides age, the CF genotype may play an important role in potential pharmacokinetic differences, as was recently shown for tobramycin pharmacokinetics (P. M. Beringer, A. A. Vinks, R. W. Jelliffe, H. G. M. Heijerman, and B. J. Shapiro, Pediatr. Pulmonol., abstr. 342, 28(S19):261, 1999). In fact, the oldest patient in our study was 31 years old but did not show a large difference between predicted and observed steady-state piperacillin concentrations. Although the patient had a positive sweat test, the genotype was not known, and it is tempting to speculate that this patient may have had a mild mutation. This mild mutation may have resulted in less of an effect or a different effect on tubular piperacillin processing.

We chose to model our data in accordance with the major renal elimination processes, tubular secretion and filtration. Reabsorption is assumed to be negligible on the basis of the pK_a values of piperacillin, which indicate that most of the drug is in its ionized form at normal urinary pHs, thus preventing reabsorption. The models in which the rate of excretion was determined by tubular secretion and glomerular filtration also performed well.

A comparison by linear regression of observed versus predicted piperacillin concentrations resulted in nearly identical plots, such as the one shown in Fig. 4, for all three nonlinear models. CL_CR is an attractive descriptor in the model that can be useful clinically. Unfortunately, only a small range of CL_CR values was available, making further validation with a larger population over a broad range of clearance values necessary.

Despite the good ability of the MAP Bayesian estimator to predict observed concentrations, as indicated by low bias and relatively good precision, the parameter estimates showed large interindividual variability. This result may have been partly due to genotypic variability, as the highest K_m and V_max values were associated with non-F508 mutations (Table 2, patients 2 and 6). Further studies are necessary to determine whether a direct relationship exists between genotype and piperacillin clearance.

In summary, for this group of CF patients, the nonlinear population modeling approach could describe piperacillin data well and was superior to the linear model. The population pharmacokinetic models can be used to design effective continuous infusion dosage regimens for the treatment of patients with CF. Models can also be used as a priori information in Bayesian feedback programs for further piperacillin dosage optimization. The modeling method using a remote supercomputer provides a new tool and a research resource for describing the behavior of large nonlinear systems. Recent new methods for access to the system and for improving the computations themselves have greatly improved the ease of access and the speed of the computations (A. Botnen, R. W. Jelliffe, M. Thomas, and N. Hoem, Abstr. 14th IEEE Symp. Comput. Based Med. Syst., p. 163, 2001; R. Leary, R. W. Jelliffe, A. Schumitzky, and M. Van Guilder, Abstr. 14th IEEE Symp. Comput. Based Med. Syst., p. 389, 2001).

 
This study was funded in part by a grant from Wyeth-Lederle, Hoofddorp, The Netherlands, and by travel grants from the Dutch Scientific Organization (NWO) and the Scientific Fund of the Royal Dutch Pharmaceutical Association (KNMP) to Alexander A. Vinks. The modeling resource used was supported by NIH grant RR11526 to the University of Southern California Laboratory of Applied Pharmacokinetics.

 
* Corresponding author. Mailing address: Pharmacology Research Center, MLC 7025, Cincinnati Children's Hospital Medical Center, 3333 Burnet Ave., Cincinnati, OH 45229. Phone: (513) 636-0159. Fax: (513) 636-0168. E-mail: sander.vinks{at}chmcc.org.

Top Abstract Introduction Materials and Methods Results Discussion References

 

1
1. Auclair, B., and M. P. Ducharme. 1999. Piperacillin and tazobactam exhibit linear pharmacokinetics after multiple standard clinical doses. Antimicrob Agents Chemother. 43:1465-1468.[Abstract/Free Full Text]
2
2. Batra, V. K., J. A. Morrison, K. C. Lasseter, and V. A. Joy. 1979. Piperacillin kinetics. Clin. Pharmacol. Ther. 26:41-53.[Medline]
3
3. Bergan, T., and J. D. Williams. 1982. Dose dependence of piperacillin pharmacokinetics. Chemotherapy 28:153-159.[Medline]
4
4. Bryson, H. M., and R. N. Brogden. 1994. Piperacillin/tazobactam. A review of its antibacterial activity, pharmacokinetic properties and therapeutic potential. Drugs 47:506-535.[Medline]
5
5. Craig, W. A., and S. C. Ebert. 1992. Continuous infusion of beta-lactam antibiotics. Antimicrob. Agents Chemother. 36:2577-2583.[Free Full Text]
6
6. DuBois, D., and E. F. DuBois. 1916. A formula to estimate the appropriate surface area if height and weight are known. Arch. Intern. Med. 17:1065-1072.
7
7. Hallynck, T. H., H. H. Soep, J. A. Thomis, J. Boelaert, R. Daneels, and L. Dettli. 1981. Should clearance be normalised to body surface or to lean body mass? Br. J. Clin. Pharmacol. 11:523-526.[Medline]
8
8. Jelliffe, R. W., and S. M. Jelliffe. 1971. Estimation of creatinine clearance from changing serum-creatinine levels. Lancet 2:710.
9
9. Jelliffe, R. W., A. Schumitzky, M. Van Guilder, M. Liu, L. Hu, P. Maire, P. Gomis, X. Barbaut, and B. Tahani. 1993. Individualizing drug dosage regimens: roles of population pharmacokinetic and dynamic models, Bayesian fitting, and adaptive control. Ther. Drug Monit. 15:380-393.[Medline]
10
10. Kucharczyk, N. 1993. Estimation of the lower limit of quantitation, a method detection performance parameter for biomedical assays, from calibration curves. J. Chromatogr. 612:71-76.[Medline]
11
11. Ludden, T. M. 1991. Nonlinear pharmacokinetics: clinical implications. Clin. Pharmacokinet. 20:429-446.[Medline]
12
12. McCarty, J. M., S. J. Tilden, P. Black, J. C. Craft, J. Blumer, W. Waring, and N. A. Halsey. 1988. Comparison of piperacillin alone versus piperacillin plus tobramycin for treatment of respiratory infections in children with cystic fibrosis. Pediatr. Pulmonol. 4:201-204.[Medline]
13
13. Proost, J. H., and D. K. Meijer. 1992. MW/Pharm, an integrated software package for drug dosage regimen calculation and therapeutic drug monitoring. Comput. Biol. Med. 22:155-163.[CrossRef][Medline]
14
14. Reed, M. D., J. Goldfarb, T. S. Yamashita, E. Lemon, and J. L. Blumer. 1994. Single-dose pharmacokinetics of piperacillin and tazobactam in infants and children. Antimicrob. Agents Chemother. 38:2817-2826.[Abstract/Free Full Text]
15
15. Reed, M. D., R. C. Stern, C. M. Myers, J. D. Klinger, T. S. Yamashita, and J. L. Blumer. 1987. Therapeutic evaluation of piperacillin for acute pulmonary exacerbations in cystic fibrosis. Pediatr. Pulmonol. 3:101-109.[Medline]
16
16. Schumitzky, A. 1991. Nonparametric EM algorithms for estimating prior distributions. Appl. Math. Comput. 45:141-157.
17
17. Sheiner, L. B., and S. L. Beal. 1981. Some suggestions for measuring predictive performance. J. Pharmacokinet. Biopharm. 9:503-512.[CrossRef][Medline]
18
18. Sorgel, F., and M. Kinzig. 1993. The chemistry, pharmacokinetics and tissue distribution of piperacillin/tazobactam. J. Antimicrob. Chemother. 31(Suppl. A):39-60.[Abstract/Free Full Text]
19
19. Tjandramaga, T. B., A. Mullie, R. Verbesselt, P. J. De Schepper, and L. Verbist. 1978. Piperacillin: human pharmacokinetics after intravenous and intramuscular administration. Antimicrob. Agents Chemother. 14:829-837.[Abstract/Free Full Text]
20
20. Touw, D. J. 1998. Clinical pharmacokinetics of antimicrobial drugs in cystic fibrosis. Pharm. World Sci. 20:149-160.[CrossRef][Medline]
21
21. Touw, D. J., A. A. Vinks, F. Jacobs, H. G. Heijerman, and W. Bakker. 1996. Creatinine clearance as predictor of tobramycin elimination in adult patients with cystic fibrosis. Ther. Drug Monit. 18:562-569.[CrossRef][Medline]
22
22. Van der Vaart, F. J. 1992. Guidelines for validation of analytical procedures in pharmaceutical and biopharmaceutical analyis. Pharm. Weekbl. 46:1229-1235.
23
23. Vinks, A. A., R. W. Brimicombe, H. G. Heijerman, and W. Bakker. 1997. Continuous infusion of ceftazidime in cystic fibrosis patients during home treatment: clinical outcome, microbiology and pharmacokinetics. J. Antimicrob. Chemother. 40:125-133.[Abstract/Free Full Text]
24
24. Vinks, A. A., J. W. Mouton, H. G. M. Heijerman, and W. Bakker. 1999. Pharmacodynamic and pharmacokinetic optimisation of treatment with ß-lactam antibiotics in cystic fibrosis. Rev. Med. Microbiol. 10:89-107.
25
25. Wang, J. P., J. D. Unadkat, S. M. al-Habet, T. A. O'Sullivan, J. Williams-Warren, A. L. Smith, and B. Ramsey. 1993. Disposition of drugs in cystic fibrosis. IV. Mechanisms for enhanced renal clearance of ticarcillin. Clin. Pharmacol. Ther. 54:293-302.[Medline]
26
26. Woodland, C., D. Blowey, S. Ito, M. Spino, and G. Koren. 1998. Hypothetical framework for enhanced renal tubular secretion of drugs in cystic fibrosis. Med. Hypotheses 51:489-491.[CrossRef][Medline]
27
27. Wysocki, M., F. Delatour, F. Faurisson, A. Rauss, Y. Pean, B. Misset, F. Thomas, J. F. Timsit, T. Similowski, H. Mentec, L. Mier, and D. Dreyfuss. 2001. Continuous versus intermittent infusion of vancomycin in severe staphylococcal infections: prospective multicenter randomized study. Antimicrob. Agents Chemother. 45:2460-2467.[Abstract/Free Full Text]
28
28. Yamaoka, K., T. Nakagawa, and T. Uno. 1978. Application of Akaike's information criterion (AIC) in the evaluation of linear pharmacokinetic equations. J. Pharmacokinet. Biopharm. 6:165-175.

---

Antimicrobial Agents and Chemotherapy, February 2003, p. 541-547, Vol. 47, No. 2
0066-4804/03/$08.00+0     DOI: 10.1128/AAC.47.2.541-547.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.

This article has been cited by other articles:

• Owens, R. C. Jr., Shorr, A. F. (2009). Rational dosing of antimicrobial agents: Pharmacokinetic and pharmacodynamic strategies. Am J Health Syst Pharm 66: S23-S30 [Abstract] [Full Text]  
• Bulitta, J. B., Duffull, S. B., Kinzig-Schippers, M., Holzgrabe, U., Stephan, U., Drusano, G. L., Sorgel, F. (2007). Systematic Comparison of the Population Pharmacokinetics and Pharmacodynamics of Piperacillin in Cystic Fibrosis Patients and Healthy Volunteers. Antimicrob. Agents Chemother. 51: 2497-2507 [Abstract] [Full Text]  
• Gavin, P. J., Suseno, M. T., Thomson, R. B. Jr., Gaydos, J. M., Pierson, C. L., Halstead, D. C., Aslanzadeh, J., Brecher, S., Rotstein, C., Brossette, S. E., Peterson, L. R. (2006). Clinical Correlation of the CLSI Susceptibility Breakpoint for Piperacillin- Tazobactam against Extended-Spectrum-{beta}-Lactamase-Producing Escherichia coli and Klebsiella Species. Antimicrob. Agents Chemother. 50: 2244-2247 [Abstract] [Full Text]  
• Li, C., Kuti, J. L., Nightingale, C. H., Mansfield, D. L., Dana, A., Nicolau, D. P. (2005). Population pharmacokinetics and pharmacodynamics of piperacillin/tazobactam in patients with complicated intra-abdominal infection. J Antimicrob Chemother 56: 388-395 [Abstract] [Full Text]  
• Lodise, T. P. Jr., Lomaestro, B., Rodvold, K. A., Danziger, L. H., Drusano, G. L. (2004). Pharmacodynamic Profiling of Piperacillin in the Presence of Tazobactam in Patients through the Use of Population Pharmacokinetic Models and Monte Carlo Simulation. Antimicrob. Agents Chemother. 48: 4718-4724 [Abstract] [Full Text]  

This Article Abstract Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Citation Map Services E-mail this article to a friend Similar articles in this journal Similar articles in ASM journals Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Reprints and Permissions Copyright Information Books from ASM Press MicrobeWorld Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Vinks, A. A. Articles by Mouton, J. W. Search for Related Content PubMed PubMed Citation Articles by Vinks, A. A. Articles by Mouton, J. W.