Roosmarijn de Cock

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.

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

Background and ObjectivesDuring 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.MethodsPopulation 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.ResultsPostmenstrual 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.ConclusionsAmikacin 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.

A Neonatal Amikacin Covariate Model Can Be Used to Predict Ontogeny of Other Drugs Eliminated Through Glomerular Filtration in Neonates

ABSTRACTPurposeRecently, 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.MethodsFive 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.ResultsThe 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.ConclusionsThis 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.

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.

Towards rational dosing algorithms for vancomycin in neonates and infants based on population pharmacokinetic modeling

ABSTRACT 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.

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.

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.

PROSPECTIVE VALIDATION OF A MODEL-BASED DOSING REGIMEN FOR AMIKACIN IN NEONATES: FEASIBLE AND RELEVANT

Introduction A neonatal amikacin dosing regimen was previously developed based on a population pharmacokinetic model. The aim of the current study was to prospectively validate this model-derived dosing regimen. Methods First, early (before and after second dose) therapeutic drug monitoring (TDM) observations were evaluated for achieving target trough (<3 mg/L) and peak (>24 mg/L) levels. Secondly, observed concentrations were compared with model-predicted concentrations, whereby the results of an NPDE (normalized prediction distribution error) were considered as well. Subsequently, Monte Carlo simulations were performed. Finally, remaining causes limiting amikacin predictability (prescription errors and disease characteristics of outliers) were explored. Results In 579 neonates [median (range) birth bodyweight 2285 (420–4850) g, postnatal age 2 (1–30) days, gestational age 34 (24–41) weeks], 90.5% of early peak levels reached 24 mg/L and 60.2% of trough levels was <3 mg/L (93.4% ≤5 mg/L). Observations were accurately predicted by the model without bias, which was confirmed by the NPDE. Monte Carlo simulations showed that peak concentrations >24 mg/L were reached in almost all patients. Trough values <3 mg/L were documented in 78–100% and 45–96% of simulated cases, respectively, when ibuprofen was co-administered or not. Suboptimal trough levels were found in patient subgroups with postnatal age <14 days and current weight >2000g. Conclusions Prospective validation of a model-based neonatal amikacin dosing regimen resulted in optimized peak and trough concentrations in almost all patients. Adapted dosing for patients with suboptimal trough levels was proposed. Besides improving dosing individualization, feasibility and relevance of neonatal prospective validation studies was demonstrated.