Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, October 2007, p. 3720-3725, Vol. 51, No. 10
0066-4804/07/$08.00+0 doi:10.1128/AAC.00318-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Wil H. F. Goessens,5
Johan W. Mouton,6
Meindert Danhof,2,3 and
John N. van den Anker4,7
Medical Centre Haaglanden, Department of Obstetrics and Gynecology, Lijnbaan 32, 2512 VA The Hague, The Netherlands,1 LAP&P Consultants BV, Archimedesweg 31, 2333 CM Leiden, The Netherlands,2 Leiden-Amsterdam Center for Drug Research, Leiden University, Division of Pharmacology, P.O. Box 9502, 2300 RA Leiden, The Netherlands,3 Erasmus MC-Sophia, Sophia Children's Hospital, Dr Molewaterplein 60, 3015 GJ Rotterdam, The Netherlands,4 Erasmus MC, University Medical Centre Rotterdam, Department of Medical Microbiology, Dr Molewaterplein 60, 3015 GJ Rotterdam, The Netherlands,5 Canisius Wilhelmina Hospital, Department of Clinical Microbiology and Infectious Diseases, Nijmegen, The Netherlands,6 Children's National Medical Center, Division of Pediatric Clinical Pharmacology, 111 Michigan Avenue, N.W., Washington, DC 200107
Received 7 March 2007/ Returned for modification 10 June 2007/ Accepted 15 July 2007
|
|
|---|
4 mg/liter, simulating the current dosing regimen of 50,000 U/kg every 12 h. This regimen is therefore adequate for the treatment of common infections in neonates on the third day of life. |
|
|---|
Infectious complications during the immediate postnatal period are not uncommon and require prompt antibiotic treatment (1, 5, 26, 28). Differences in body composition and organ function can significantly affect the pharmacokinetics of drugs in neonates. In very premature neonates (i.e., neonates with gestational ages of less than 32 weeks), the disposition of antibiotics may differ from that in full-term neonates as a result of differences in absorption, distribution, biotransformation, and excretion (6, 17, 27). As a result, dose estimation on the basis of body size or allometric scaling may be inadequate. Specifically, in very premature neonates the processes of maturation of the organs may influence the relevant pharmacokinetic parameters.
The efficacies of the penicillins are primarily correlated to the percentage of time that the concentrations of unbound drug remained above the MIC (fT > MIC) (9, 18, 30). In general, the therapeutic goal for the cure of infections caused by gram-positive organisms is an fT > MIC of the antimicrobial of at least 40%, which corresponds to an in vivo static effect in animal studies (8). Studies showing a clear relationship between exposure and efficacy in premature neonates are not available. Because premature neonates must be regarded as immunocompromised, in our opinion it is reasonable that in this patient group the exposure should correspond to exposures that correlate to a 1- to 2-log drop in the numbers of CFU in various models and thus be bactericidal rather than bacteriostatic. Thus, for penicillins, this percentage should be at least 50% (11). Since the goal of treatment is to attain this target for every individual in the population, the dosing regimen in this age group should be defined by taking the interindividual variability in pharmacokinetics into account. Monte Carlo simulation (MCS) is a technique that is commonly used to determine the probability of achieving therapeutic concentrations on the basis of population pharmacokinetic parameter estimates and their measures of dispersion (3, 12, 13, 21-23). We investigated the pharmacokinetics of penicillin G and the adequacy of the dosing regimen in very premature neonates on the third day of life to allow us to construct a population pharmacokinetic model and to use the parameter estimates to perform MCSs with this specific age group. We then used that information to define the optimal dosing regimens.
|
|
|---|
Penicillin G was administered as an intravenous bolus injection of 50,000 U/kg every 12 h. For each subject, the leukocyte count, platelet count, and blood and superficial specimen cultures were performed as part of the routine workup.
Pharmacokinetics study. The pharmacokinetics of penicillin G were studied on day 3 of life. Blood samples (200 µl) were taken from an indwelling arterial line just before the administration of an intravenous bolus dose and at 0.03, 0.5, 1, 2.5, 4, 8, and 12 h after administration. A sample was taken at 24 h from those patients who did not receive a subsequent dose. Samples were immediately centrifuged at 3,000 x g in a microcentrifuge (Merck-type Eppendorf 5414) for 1 min, and the serum was stored at –70°C.
Penicillin G high-pressure liquid chromatography assay. Chromatographic analysis was performed with a glass-prepacked column (100 by 3 mm) containing ODS-2 Chromospher Spherisorb beads (5-µm-diameter particle size; Chrompack, Middelburg, The Netherlands) combined with a guard column. A Bio LC pump (model 410; Perkin-Elmer, Norwalk, Conn.) was used to deliver the eluent, which consisted of 16% (vol/vol) acetonitrile and 50 mM sodium phosphate buffer (pH 6.9), at a flow rate of 0.8 ml/min. The separations were carried out at room temperature. The eluate was monitored with a Perkin-Elmer LC-95 UV/visible spectrophotometer detector at a wavelength of 215 nm. As an internal standard 25 µg/ml methicillin in 100% (vol/vol) methanol was used. Briefly, 100 µl of the internal standard was added to a 100-µl aliquot of the serum sample. This mixture was immediately vortexed for 30 s. Subsequently, the sample was kept at –20°C for 10 min, again vortexed for 30 s, and finally, centrifuged at 1,500 x g for 10 min at room temperature. The supernatant was filtered (Millipore), and 10 µl was injected onto the column.
High-pressure liquid chromatography-grade acetonitrile was purchased from Rathburn (Walkerurb, Scotland). The other chemicals were purchased from Aldrich-Chemie (Steinheim, Germany). All chemicals applied were of the highest grade commercially available.
The lower limit of detection of penicillin was 0.5 µg/ml. The coefficients of interassay variation determined at concentrations of 100 and 20 µg/ml were 2.6% and 2.3%, respectively. The intra-assay values were 0.75% and 1.05%, respectively.
Pharmacokinetic analysis. Pharmacokinetic parameters were estimated by means of nonlinear mixed effect (population) modeling (NONMEM). This approach estimates the structural pharmacokinetic parameters by considering both the interindividual variability within the population and the intraindividual (i.e., residual) variability. The model was implemented in the NONMEM ADVAN5 subroutine, and the analysis was performed by using the first-order conditional estimation method with the interaction option. All fitting procedures were performed with the use of the Compaq Visual FORTRAN, standard edition 6.6 (Compaq Computer Cooperation, Euston, TX), and NONMEM, version V (NONMEM project group; University of California, San Francisco).
To determine the basic structural pharmacokinetic parameters, several models were evaluated. One-, two-, and three-compartment models were tested and evaluated for their goodness of fit. Model selection and the identification of variability were based on the likelihood ratio test, pharmacokinetic parameter point estimates, and their respective confidence intervals and goodness-of-fit plots. For the likelihood ratio test for differences between two models, the objective function value with a prespecified level of significance of a P value of <0.001 was used. NONMEM minimizes an objective function in performing nonlinear regression analysis. To detect systematic deviations in the model fits, the goodness-of-fit plots were visually inspected. The data of individual observations versus individual or population predictions should be randomly distributed around the line of identity. The weighted residuals versus time or population predictions should be randomly distributed around zero.
The stochastic part of the model was selected to describe the interindividual variability in the pharmacokinetic parameters, and a lognormal distribution of all model parameters over the population was assumed. Therefore, an exponential distribution model was used to account for interindividual variability, as follows: Pi =
· exp(
i), where Pi is the value of the model parameter P for individual i,
is the population estimate for parameter P, and
i is the normally distributed interindividual random variable with mean zero and variance
2. Selection of an appropriate residual error model was based on the likelihood ratio test and inspection of the goodness-of-fit plots. The model was modified to objectively account for unexplained inconsistencies in the data. To reduce the influence of neonates with large unexplained inconsistencies in the concentration-time profile on the population estimates, these neonates were objectively determined by means of the residual error and were weighted less in the estimates for the population.
To refine the stochastic model, covariate analysis was also performed. The estimated pharmacokinetic parameters were plotted independently against the covariates gestational age, birth weight, gender, and the presence of a dosing history to determine whether these influenced the pharmacokinetics. The effects of the covariates were tested for statistical significance by the likelihood ratio test, and the residual intra- and interindividual variabilities were visually evaluated. A covariate was retained in the model if it produced a decrease in objective function of >10.8 (P < 0.001). In addition, we investigated whether there were significant differences in the pharmacokinetics between neonates with birth weights of less than 1,000 g and neonates with birth weights of more than 1,000 g. The volume of distribution at steady state (Vss) and the half-life (t1/2) were calculated by standard procedures (15).
Estimation of fT > MIC and MCSs. The estimates of the pharmacokinetic parameters and measures of dispersion were used to simulate various dosing regimens and obtain the percent fT > MIC as a function of the MIC (23). Protein binding was estimated to be 40% ± 2.5% (14). The protein binding of penicillin G in premature neonates is unknown. However, protein binding in neonates is generally less than that in adults. It is therefore likely that the estimated protein binding of 40% overestimates the nonactive protein-bound fraction of penicillin G in these neonates. The use of 40% is therefore a conservative estimate. MCS was performed by using the MICLAB (version 2.36) program (Medimatics, Maastricht, The Netherlands) by simulating 10,000 subjects for each regimen. The program allows inclusion of the covariance matrix (or correlation matrix) of the parameter estimates used in the simulations. The output consisted of a probability distribution, a cumulative probability distribution, and specific confidence intervals over user-defined MIC and percent fT > MIC ranges.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Demographic, laboratory, and clinical parameters for 20 patients studied on day 3 after birth
|
![]() View larger version (18K): [in a new window] |
FIG. 1. Individual plots for the 20 very preterm neonates. The black dots correspond to the individual datum points. The line represents the individual estimate, and the dotted line represents the population estimate. Neonates 1, 3, and 6 were weighted less in the population estimation. ID, neonate number.
|
![]() View larger version (17K): [in a new window] |
FIG. 2. Plot of individual predicted versus observed concentrations of penicillin G for 20 patients. The correlation coefficient was 0.801. The individual datum points for the entire population and the x-y line are also shown.
|
|
View this table: [in a new window] |
TABLE 2. Pharmacokinetic parameters for penicillin G in 20 preterm neonates
|
![]() View larger version (11K): [in a new window] |
FIG. 3. Observed relation between body weight and clearance.
|
4 mg/liter. Figure 4 shows the percent fT > MICs for the dose of 50,000 U/kg and three different dosing intervals on the basis of the mean population parameter estimates and the correlation matrix of these estimates.
![]() View larger version (12K): [in a new window] |
FIG. 4. fT > MIC for penicillin G, based on the pharmacokinetic estimates and the correlation matrix of the parameter estimates, as a function of the MICs for three regimens.
|
|
|
|---|
Most studies on the pharmacokinetics of penicillin G in neonates have been performed after intramuscular administration (19, 20, 25). However, the intramuscular route should be avoided, because this may result in erratic absorption in a sick, infected newborn with a restricted blood supply to the extremities (31). Mulhall (25) performed a study on the pharmacokinetics of penicillin G after intravenous administration in four neonates with gestational ages ranging from 27 to 40 weeks and found a CL of 0.12 ± 0.07 liter/h/kg (mean ± standard deviation), a V of 0.61 ± 0.28 liter/kg (mean ± standard deviation), and a t1/2 of 3.8 h. These results are similar to our data.
Other penicillins have larger V and longer t1/2ß values in premature neonates, especially compared with those in adults (11). V varied from 0.3 liter/kg for ampicillin to 0.41 to 0.68 liter/kg for amoxicillin in neonates (2, 16, 32), whereas V was 0.45 liter/kg for penicillin G in our study. The t1/2ß of penicillin G in neonates with gestational ages of less than 32 weeks was longer than the value of 0.5 h reported for healthy adults (14), but it is in the same range as the values of between 2 and 9.5 h that have been reported for other penicillins in neonates (2, 7, 10, 11, 16, 32).
Among the limited number of samples available from patients in this age group, there was some unexplained inconsistency in the data. This might have been caused by the subcutaneous administration of penicillin G, erratic sampling times, or the accidental exchange of samples. Large unexplained inconsistencies in the data were found for three neonates, based objectively on their residual errors, and therefore, the data were weighted less in the population estimates. By using this method, none of the neonates were excluded from the study and all data were used in the analysis. Consequently, there is a deviation between the line of identity and the regression line of the observed versus the predicted concentrations, indicating that the description of the pharmacokinetics in this age group can be improved. To this end, more data are needed.
The presence of a third elimination phase (t1/2
) for penicillin G has been described previously both in animal models and in adult humans (14, 29). Our data could not be described by a three-compartment model. The t1/2ß found in our study, however, was comparable to the t1/2
of 3.1 h in human adults (14). Both the limited number of samples taken in the initial distribution phase and the unique body composition of the neonate, which comprises approximately 75% water, complicate the distinction between the initial distribution and the second elimination phase. It is possible that the initial phase in our study represents both the initial and second phases found in the study of Ebert et al. (14). Thus, the slow elimination that we found may represent the third elimination phase for penicillin G determined in that study.
Interindividual variability was partly explained by variations in CL, V2, and the infusion rate. The variability in the infusion rate represents not only the variation in the rate of manual administration of the intravenous bolus injection but also the variation in the sampling times between the first samples. Especially for the first samples, the exact sampling times are crucial.
Growth and development are major aspects in infants; therefore, both size and gestational age may have an impact on the prediction of CL (4). The maturation of CL begins before birth, suggesting that gestational age would be a physiologically appropriate covariate to explain the time course of changes in CL (4). In our data, changes in CL were best explained by differences in body weight, probably because both gestational age and other factors influencing the development of renal function are represented by an adequate increase in body weight in time. Furthermore, the range of gestational ages of the nenonates included in our study was relatively small. While CL slightly increased as a function of birth weight for the entire group, differences between the subgroup with birth weights of more than 1,000 g and the subgroup with birth weights of less than 1,000 g could not be demonstrated. The analysis of correlations between covariates and pharmacokinetic parameter estimates is more sensitive when the data for all neonates are included in the study.
The enhanced interindividual pharmacokinetic variability in premature neonates complicates the calculation of the therapeutic dosing regimen. The dosing regimen should be adequate for the entire population. We concluded that the 100% PTA obtained with simulation of the recommended regimen was adequate for the treatment of infections in neonates on the third day of life. In case of meningeal involvement, effective concentrations in the cerebrospinal fluid (CSF) are required. Little is known about the pharmacokinetics of penicillin G in the CSF of premature neonates. We do not know what percentage of penicillin G penetrates the CSF in premature neonates with meningitis. Given the low penicillin MIC for group B streptococci (up to 0.12 mg/liter), with the use of the currently recommended dosing regimen, it is sufficient when the level of penetration of penicillin into the CSF is at least 3%.
Shortening of the dosing interval to 8 h will not have additional value for the treatment of infections in neonates with gestational ages of less than 32 weeks. When these infections are caused by microorganisms with low MICs, like Streptococcus agalactiae, a regimen with a prolonged dosing interval of 24 h is also likely to have clinical success. For empirical therapy, however, this regimen is suboptimal.
Published ahead of print on 23 July 2007. ![]()
Present address: Wilhelmina Hospital Assen, Department of Pediatrics, Europaweg-Zuid 1, 9401 RK Assen, The Netherlands. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»