Towards Rational Dosing Algorithms for Vancomycin in Neonates and Infants Based on Population Pharmacokinetic Modeling.

Because of the recent awareness that vancomycin doses should aim to meet a target area under the concentration-time curve (AUC) instead of trough concentrations, more aggressive dosing regimens are warranted also in the pediatric population. In this study, both neonatal and pediatric pharmacokinetic models for vancomycin were externally evaluated and subsequently used to derive model-based dosing algorithms for neonates, infants, and children. For the external validation, predictions from previously published pharmacokinetic models were compared to new data. Simulations were performed in order to evaluate current dosing regimens and to propose a model-based dosing algorithm. The AUC/MIC over 24 h (AUC24/MIC) was evaluated for all investigated dosing schedules (target of >400), without any concentration exceeding 40 mg/liter. Both the neonatal and pediatric models of vancomycin performed well in the external data sets, resulting in concentrations that were predicted correctly and without bias. For neonates, a dosing algorithm based on body weight at birth and postnatal age is proposed, with daily doses divided over three to four doses. For infants aged <1 year, doses between 32 and 60 mg/kg/day over four doses are proposed, while above 1 year of age, 60 mg/kg/day seems appropriate. As the time to reach steady-state concentrations varies from 155 h in preterm infants to 36 h in children aged >1 year, an initial loading dose is proposed. Based on the externally validated neonatal and pediatric vancomycin models, novel dosing algorithms are proposed for neonates and children aged <1 year. For children aged 1 year and older, the currently advised maintenance dose of 60 mg/kg/day seems appropriate.

Novel model-based dosing guidelines for gentamicin and tobramycin in preterm and term neonates.

OBJECTIVES: In the heterogeneous group of preterm and term neonates, gentamicin and tobramycin are mainly dosed according to empirical guidelines, after which therapeutic drug monitoring and subsequent dose adaptation are applied. In view of the variety of neonatal guidelines available, the purpose of this study was to evaluate target concentration attainment of these guidelines, and to propose a new model-based dosing guideline for these drugs in neonates. METHODS: Demographic characteristics of 1854 neonates (birth weight 390-5200 g, post-natal age 0-27 days) were extracted from earlier studies and sampled to obtain a test dataset of 5000 virtual patients. Monte Carlo simulations on the basis of validated models were undertaken to evaluate the attainment of target peak (5-12 mg/L) and trough (<0.5 mg/L) concentrations, and cumulative AUC, with the existing and proposed guidelines. RESULTS: Across the entire neonatal age and weight range, the Dutch National Formulary for Children, the British National Formulary for Children, Neofax and the Red Book resulted in adequate peak but elevated trough concentrations (63%-90% above target). The proposed dosing guideline (4.5 mg/kg gentamicin or 5.5 mg/kg tobramycin) with a dosing interval based on birth weight and post-natal age leads to adequate peak concentrations with only 33%-38% of the trough concentrations above target, and a constant AUC across weight and post-natal age. CONCLUSIONS: The proposed neonatal dosing guideline for gentamicin and tobramycin results in improved attainment of target concentrations and should be prospectively evaluated in clinical studies to evaluate the efficacy and safety of this treatment.

Simultaneous pharmacokinetic modeling of gentamicin, tobramycin and vancomycin clearance from neonates to adults: towards a semi-physiological function for maturation in glomerular filtration.

PURPOSE: Since glomerular filtration rate (GFR) is responsible for the elimination of a large number of water-soluble drugs, the aim of this study was to develop a semi-physiological function for GFR maturation from neonates to adults. METHODS: In the pharmacokinetic analysis (NONMEM VI) based on data of gentamicin, tobramycin and vancomycin collected in 1,760 patients (age 1 day-18 years, bodyweight 415 g-85 kg), a distinction was made between drug-specific and system-specific information. Since the maturational model for clearance is considered to contain system-specific information on the developmental changes in GFR, one GFR maturational function was derived for all three drugs. RESULTS: Simultaneous analysis of these three drugs showed that maturation of GFR mediated clearance from preterm neonates to adults was best described by a bodyweight-dependent exponent (BDE) function with an exponent varying from 1.4 in neonates to 1.0 in adults (ClGFR = Cldrug*(BW/4 kg)(BDE) with BDE = 2.23*BW(-0.065)). Population clearance values (Cldrug) for gentamicin, tobramycin and vancomycin were 0.21, 0.28 and 0.39 L/h for a full term neonate of 4 kg, respectively. DISCUSSION: Based on an integrated analysis of gentamicin, tobramycin and vancomycin, a semi-physiological function for GFR mediated clearance was derived that can potentially be used to establish evidence based dosing regimens of renally excreted drugs in children.

Low but inducible contribution of renal elimination to clearance of propylene glycol in preterm and term neonates.

BACKGROUND: Despite limited information being available on the pharmacokinetics of excipients, propylene glycol (PG) is often used as an excipient in both adults and children. The aim of this study is to characterize the renal and hepatic elimination of PG in preterm and term neonates. METHODS: The pharmacokinetic analysis of PG was performed in NONMEM 6.2. on the basis of PG concentrations in plasma and/or urine samples for a total of 69 (pre)term neonates (birth weight 630-3980 g, gestational age 24-41 weeks, postnatal age 1-29 days) who received PG coadministered with intravenous paracetamol (5-10 mg/kg per 6 hours), phenobarbital (5 mg·kg(-1)·d(-1)), or both. To capture the time-dependent trend in the renal excretion of PG, different models based on time after the first dose, urine volume, and creatinine amount in urine were tested. RESULTS: A one-compartment model parameterized in terms of renal clearance, hepatic clearance, and volume of distribution was found to adequately describe the observations in both plasma and urine. After the first dose was administered, the renal elimination of PG was 15% of total clearance, which increased over time to 25% at 24 hours after the first dose of PG. This increase was best described using a hyperbolic function based on time after the first dose. CONCLUSIONS: Renal elimination of PG in (pre)term neonates is low, particularly compared with the reported percentage of 45% in adults, but it may increase with time after the first dose of PG. To study whether this increase is caused by an autoinduced increase in the renal secretion or a reduction of tubular reabsorption of PG, further research is needed.

A neonatal amikacin covariate model can be used to predict ontogeny of other drugs eliminated through glomerular filtration in neonates.

PURPOSE: Recently, a covariate model characterizing developmental changes in clearance of amikacin in neonates has been developed using birth bodyweight and postnatal age. The aim of this study was to evaluate whether this covariate model can be used to predict maturation in clearance of other renally excreted drugs. METHODS: Five different neonatal datasets were available on netilmicin, vancomycin, tobramycin and gentamicin. The extensively validated covariate model for amikacin clearance was used to predict clearance of these drugs. In addition, independent reference models were developed based on a systematic covariate analysis. RESULTS: The descriptive and predictive properties of the models developed using the amikacin covariate model were good, and fairly similar to the independent reference models (goodness-of-fit plots, NPDE). Moreover, similar clearance values were obtained for both approaches. Finally, the same covariates as in the covariate model of amikacin, i.e. birth bodyweight and postnatal age, were identified on clearance in the independent reference models. CONCLUSIONS: This study shows that pediatric covariate models may contain physiological information since information derived from one drug can be used to describe other drugs. This semi-physiological approach may be used to optimize sparse data analysis and to derive individualized dosing algorithms for drugs in children.

Developmental pharmacokinetics of propylene glycol in preterm and term neonates.

AIM: Propylene glycol (PG) is often applied as an excipient in drug formulations. As these formulations may also be used in neonates, the aim of this study was to characterize the pharmacokinetics of propylene glycol, co-administered intravenously with paracetamol (800 mg PG/1000 mg paracetamol) or phenobarbital (700 mg PG/200 mg phenobarbital) in preterm and term neonates. METHODS: A population pharmacokinetic analysis was performed based on 372 PG plasma concentrations from 62 (pre)term neonates (birth weight (bBW) 630-3980 g, postnatal age (PNA) 1-30 days) using NONMEM 6.2. The model was subsequently used to simulate PG exposure upon administration of paracetamol or phenobarbital in neonates (gestational age 24-40 weeks). RESULTS: In a one compartment model, birth weight and PNA were both identified as covariates for PG clearance using an allometric function (CL(i) = 0.0849 × {(bBW/2720)(1.69) × (PNA/3)(0.201)}). Volume of distribution scaled allometrically with current bodyweight (V(i) = 0.967 × {(BW/2720)(1.45)}) and was estimated 1.77 times higher when co-administered with phenobarbital compared with paracetamol. By introducing these covariates a large part of the interindividual variability on clearance (65%) as well as on volume of distribution (53%) was explained. The final model shows that for commonly used dosing regimens, the population mean PG peak and trough concentrations range between 33-144 and 28-218 mg l(-1) (peak) and 19-109 and 6-112 mg l(-1) (trough) for paracetamol and phenobarbital formulations, respectively, depending on birth weight and age of the neonates. CONCLUSION: A pharmacokinetic model was developed for PG co-administered with paracetamol or phenobarbital in neonates. As such, large variability in PG exposure may be expected in neonates which is dependent on birth weight and PNA.

Is indirect hyperbilirubinemia a useful biomarker of reduced propofol clearance in neonates?

AIM: Large interindividual variability in neonatal propofol clearance is documented which, in part, can be explained by postmenstrual age (PMA) and postnatal age (PNA). We aimed to document whether indirect bilirubin, instead of or in addition to PNA, could improve predictability of propofol clearance and serve as a useful biomarker of reduced propofol clearance in neonates. METHODS: Indirect serum bilirubin was introduced as a dichotomous or continuous variable (both age-normalized) in a previously developed three-compartment pharmacokinetic model, based on 235 concentration-time points obtained in 25 neonates after single bolus administration of propofol. For pharmacokinetic analysis, nonlinear mixed effect modeling 6.2 was used. RESULTS: The covariates PMA and PNA explained 67% of the interindividual variability compared with 45% in the model with PMA and bilirubin. CONCLUSION: Age, reflected by PMA and PNA, is a more relevant clinical predictor of neonatal propofol clearance compared with PMA and raised indirect hyperbilirubinemia.

Maturation of the glomerular filtration rate in neonates, as reflected by amikacin clearance.

BACKGROUND AND OBJECTIVES: During the newborn period and early infancy, renal function matures, resulting in changes in the glomerular filtration rate (GFR). This study was performed to quantify developmental changes in the GFR in (pre)term neonates by use of amikacin clearance as proof of concept. The model was used to derive a rational dosing regimen in comparison with currently used dosing regimens for amikacin. METHODS: Population pharmacokinetic modelling was performed in nonlinear mixed-effect modelling software (NONMEM version 6.2) using data from 874 neonates obtained from two previously published datasets (gestational age 24-43 weeks; postnatal age 1-30 days; birthweight 385-4650 g). The influence of different age-related, weight-related and other covariates was investigated. The model was validated both internally and externally. RESULTS: Postmenstrual age was identified as the most significant covariate on clearance. However, the combination of birthweight and postnatal age proved to be superior to postmenstrual age alone. Birthweight was best described using an allometric function with an exponent of 1.34. Postnatal age was identified using a linear function with a slope of 0.2, while co-administration of ibuprofen proved to be a third covariate. Current bodyweight was the most important covariate for the volume of distribution, using an allometric function. The external evaluation supported the prediction of the final pharmacokinetic model. This analysis illustrated clearly that the currently used dosing regimens for amikacin in reference handbooks may possibly increase the risk of toxicities and should be revised. Consequently, a new model-based dosing regimen based on current bodyweight and postnatal age was derived. CONCLUSIONS: Amikacin clearance, reflecting the GFR in neonates, can be predicted by birthweight representing the antenatal state of maturation of the kidney, postnatal age representing postnatal maturation, and co-administration of ibuprofen. Finally, the model reflects maturation of the GFR, allowing for adjustments of dosing regimens for other renally excreted drugs in preterm and term neonates.

The role of population PK-PD modelling in paediatric clinical research.

Children differ from adults in their response to drugs. While this may be the result of changes in dose exposure (pharmacokinetics [PK]) and/or exposure response (pharmacodynamics [PD]) relationships, the magnitude of these changes may not be solely reflected by differences in body weight. As a consequence, dosing recommendations empirically derived from adults dosing regimens using linear extrapolations based on body weight, can result in therapeutic failure, occurrence of adverse effect or even fatalities. In order to define rational, patient-tailored dosing schemes, population PK-PD studies in children are needed. For the analysis of the data, population modelling using non-linear mixed effect modelling is the preferred tool since this approach allows for the analysis of sparse and unbalanced datasets. Additionally, it permits the exploration of the influence of different covariates such as body weight and age to explain the variability in drug response. Finally, using this approach, these PK-PD studies can be designed in the most efficient manner in order to obtain the maximum information on the PK-PD parameters with the highest precision. Once a population PK-PD model is developed, internal and external validations should be performed. If the model performs well in these validation procedures, model simulations can be used to define a dosing regimen, which in turn needs to be tested and challenged in a prospective clinical trial. This methodology will improve the efficacy/safety balance of dosing guidelines, which will be of benefit to the individual child.