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

Pharmacokinetic and Pharmacodynamic Approaches to Optimize Antibiotic Use in Neonates

Newborn infants (particularly those born preterm) are frequently exposed to empiric antibiotics at birth, and antibiotics are among the most commonly prescribed medications in neonatal intensive care units. Challenges in optimizing neonatal antibiotic dosing include: technical and ethical barriers to neonatal pharmacoanalytic study design and sampling, difficulty in extrapolating adult and pediatric data due to unique neonatal physiology, and a lack of validated pharmacodynamic targets specific to neonatal populations. In this review, we summarize basic concepts in pharmacokinetics (PK) and pharmacodynamics (PD), describe pharmacometric strategies utilized in contemporary PK/PD analyses, and review the evolution of PK/PD data guiding neonatal dosing among 3 commonly used antibiotics.

PK/PD-Guided Strategies for Appropriate Antibiotic Use in the Era of Antimicrobial Resistance

Antimicrobial resistance (AMR) poses a critical global health threat, necessitating the optimal use of existing antibiotics. Pharmacokinetic/pharmacodynamic (PK/PD) principles provide a scientific framework for optimizing antimicrobial therapy, particularly to respond to evolving resistance patterns. This review examines PK/PD strategies for antimicrobial dosing optimization, focusing on three key aspects. First, we discuss the importance of drug concentration management for enhancing efficacy while preventing toxicity, considering various patient populations, including pediatric and elderly patients with their unique physiological characteristics. Second, we analyze different PK modeling approaches: the classic top-down approach exemplified by population PK analysis, the bottom-up approach represented by physiologically based PK modeling, and hybrid models combining both approaches for enhanced predictive performance. Third, we explore clinical applications, including nomogram-based dosing strategies, Bayesian estimation, and emerging artificial intelligence applications, for real-time dose optimization. Critical challenges in implementing PK/PD simulation are addressed, particularly the selection of appropriate PK models, the optimization of PK/PD indices, and considerations concerning antimicrobial concentrations at infection sites. Understanding these principles and challenges is crucial for optimizing antimicrobial therapy and combating AMR through improved dosing strategies.

A Quantitative Examination and Comparison of the Ability of Australian Gentamicin Dosing Guidelines to Achieve Target Therapeutic Concentrations in Neonates

Background: Effective gentamicin dosing is crucial to the survival of neonates with suspected sepsis but requires a careful balance between attaining both effective peak and safe trough concentrations. We aimed to systematically compare existing gentamicin dosing guidelines for neonates in Australia to determine the extent to which they reach therapeutic targets. Methods: Simulations of a single gentamicin dose to a virtual representative neonatal population according to each Australian guideline were performed using population pharmacokinetic modelling. We determined the proportion of neonates who would achieve peak gentamicin concentrations of ≥5 or ≥10 mg/L and trough concentrations of ≤1 or ≤2 mg/L. We calculated the probability of target attainment (PTA) according to gestation at birth (22 to 40 weeks) and postnatal age (1-7, 8-14, 15-21, 22-28 days). Results: Five unique dosing guidelines were identified. Guidelines varied considerably with respect to dose (4.5 to 7 mg/kg), dosing interval (24 to 48 h), and characteristics used to individualise dosing regimens (e.g., gestation at birth and postnatal age). Guidelines were satisfactory in routinely achieving effective peak concentrations ≥ 5 mg/L, but PTAs for effective peak concentrations ≥ 10 mg/L varied considerably from 5% to 100% based on dose, gestation, and postnatal age. PTAs for trough concentrations ≤ 1 mg/L ranged from 0% to 100%, being lowest among extremely preterm infants. Conclusions: Current Australian gentamicin guidelines demonstrate significant variability in their ability to achieve defined therapeutic targets, necessitating efforts to improve standardisation of dosing recommendations. Further research to define optimal pharmacodynamic targets in neonates with respect to clinical outcomes are also urgently warranted.

Postdiscontinuation Antibiotic Exposure in Hospitalized Infants at Risk for Late-onset Sepsis in the Neonatal Intensive Care Unit

BACKGROUND

In the neonatal intensive care unit, infants are at risk for late-onset sepsis. When blood cultures are negative, antibiotic stewardship efforts encourage stopping antibiotics, yet the duration of therapeutic exposure after the last dose is unknown.

METHODS

This retrospective cohort study of simulated antibiotic exposures used published population pharmacokinetic models within drug-specific neonatal intensive care unit cohorts of preterm and term infants, postnatal age 7-60 days and exposed to cefepime, piperacillin-tazobactam or tobramycin. Monte Carlo simulations (NONMEM 7.3) were used to predict steady-state exposures after a 72-hour antibiotic course per Neofax dosing. Exposure was assessed relative to drug-specific minimum inhibitory concentration (MIC) targets between 1 and 16 mcg/mL for Pseudomonas and Enterobacteriaceae species. Postdiscontinuation antibiotic exposure (PDAE) was defined as the time from the last dose to when antibiotic concentration decreased below a specific MIC.

RESULTS

Piperacillin-tazobactam, cefepime and tobramycin cohorts included infants with median gestation age 29, 32 and 32 weeks and postnatal age 17, 19 and 15 days, respectively. The mean PDAE was 19-68 hours, depending on the specific antibiotic/MIC combination. PDAE was longer for infants <28 days old and preterm (vs. term) infants. Cefepime exhibited the longest mean PDAE of 68 hours for Enterobacteriaceae MIC 1. Piperacillin mean PDAE was 25 hours for Enterobacteriaceae MIC 8. Tobramycin had a short mean PDAE of 19 hours.

CONCLUSIONS

Piperacillin and cefepime exposures remained therapeutic long after the expected 8- to 12-hour dosing interval. PDAE is an important consideration for antibiotic stewardship among hospitalized infants, particularly premature infants and those within 1 month postbirth.

Personalized Dosing of Medicines for Children: A Primer on Pediatric Pharmacometrics for Clinicians

The widespread use of drugs for unapproved purposes remains common in children, primarily attributable to practical, ethical, and financial constraints associated with pediatric drug research. Pharmacometrics, the scientific discipline that involves the application of mathematical models to understand and quantify drug effects, holds promise in advancing pediatric pharmacotherapy by expediting drug development, extending applications, and personalizing dosing. In this review, we delineate the principles of pharmacometrics, and explore its clinical applications and prospects. The fundamental aspect of any pharmacometric analysis lies in the selection of appropriate methods for quantifying pharmacokinetics and pharmacodynamics. Population pharmacokinetic modeling is a data-driven method ('top-down' approach) to approximate population-level pharmacokinetic parameters, while identifying factors contributing to inter-individual variability. Model-informed precision dosing is increasingly used to leverage population pharmacokinetic models and patient data, to formulate individualized dosing recommendations. Physiologically based pharmacokinetic models integrate physicochemical drug properties with biological parameters ('bottom-up approach'), and is particularly valuable in situations with limited clinical data, such as early drug development, assessing drug-drug interactions, or adapting dosing for patients with specific comorbidities. The effective implementation of these complex models hinges on strong collaboration between clinicians and pharmacometricians, given the pivotal role of data availability. Promising advancements aimed at improving data availability encompass innovative techniques such as opportunistic sampling, minimally invasive sampling approaches, microdialysis, and in vitro investigations. Additionally, ongoing research efforts to enhance measurement instruments for evaluating pharmacodynamics responses, including biomarkers and clinical scoring systems, are expected to significantly bolster our capacity to understand drug effects in children.

Population pharmacokinetic analysis for dose regimen optimization of vancomycin in Southern Chinese children

Changes in physiological factors may result in large pharmacokinetic variability of vancomycin in pediatric patients, thereby leading to either supratherapeutic or subtherapeutic exposure and potentially affecting clinical outcomes. This study set out to characterize the disposition of vancomycin, quantify the exposure target and establish an optimal dosage regimen among the Southern Chinese pediatric population. Routine therapeutic drug monitoring data of 453 patients were available. We performed a retrospective population pharmacokinetic analysis of hospitalized children prescribed intravenous vancomycin using NONMEM® software. A one-compartment PPK model of vancomycin with body weight and renal functions as covariates based on a cutoff of 2 years old children was proposed in this study. Both internal and external validation showing acceptable and robust predictive performance of the model to estimate PK parameters. The value of area under the curve over 24 h to minimum inhibitory concentration ratio (AUC0-24/MIC) ≥ 260 was a significant predictor for therapeutic efficacy. Monte Carlo simulations served as a model-informed precision dosing approach and suggested that different optimal dose regimens in various scenarios should be considered rather than flat dosing. The evaluation of vancomycin exposure-efficacy relationship indicated that lower target level of AUC0-24/MIC may be needed to achieve clinical effectiveness in children, which was used to derive the recommended dosing regimen. Further prospective studies will be needed to corroborate and elucidate these results.

Virtual Patient Simulation Using Copula Modeling

Virtual patient simulation is increasingly performed to support model-based optimization of clinical trial designs or individualized dosing strategies. Quantitative pharmacological models typically incorporate individual-level patient characteristics, or covariates, which enable the generation of virtual patient cohorts. The individual-level patient characteristics, or covariates, used as input for such simulations should accurately reflect the values seen in real patient populations. Current methods often make unrealistic assumptions about the correlation between patient's covariates or require direct access to actual data sets with individual-level patient data, which may often be limited by data sharing limitations. We propose and evaluate the use of copulas to address current shortcomings in simulation of patient-associated covariates for virtual patient simulations for model-based dose and trial optimization in clinical pharmacology. Copulas are multivariate distribution functions that can capture joint distributions, including the correlation, of covariate sets. We compare the performance of copulas to alternative simulation strategies, and we demonstrate their utility in several case studies. Our work demonstrates that copulas can reproduce realistic patient characteristics, both in terms of individual covariates and the dependence structure between different covariates, outperforming alternative methods, in particular when aiming to reproduce high-dimensional covariate sets. In conclusion, copulas represent a versatile and generalizable approach for virtual patient simulation which preserve relationships between covariates, and offer an open science strategy to facilitate re-use of patient data sets.

When will the Glomerular Filtration Rate in Former Preterm Neonates Catch up with Their Term Peers?

AIMS

Whether and when glomerular filtration rate (GFR) in preterms catches up with term peers is unknown. This study aims to develop a GFR maturation model for (pre)term-born individuals from birth to 18 years of age. Secondarily, the function is applied to data of different renally excreted drugs.

METHODS

We combined published inulin clearance values and serum creatinine (Scr) concentrations in (pre)term born individuals throughout childhood. Inulin clearance was assumed to be equal to GFR, and Scr to reflect creatinine synthesis rate/GFR. We developed a GFR function consisting of GFRbirth (GFR at birth), and an Emax model dependent on PNA (with GFRmax, PNA50 (PNA at which half ofGFR max is reached) and Hill coefficient). The final GFR model was applied to predict gentamicin, tobramycin and vancomycin concentrations.

RESULT

In the GFR model, GFRbirth varied with birthweight linearly while in the PNA-based Emax equation, GA was the best covariate for PNA50, and current weight for GFRmax. The final model showed that for a child born at 26 weeks GA, absolute GFR is 18%, 63%, 80%, 92% and 96% of the GFR of a child born at 40 weeks GA at 1 month, 6 months, 1 year, 3 years and 12 years, respectively. PopPK models with the GFR maturation equations predicted concentrations of renally cleared antibiotics across (pre)term-born neonates until 18 years well.

CONCLUSIONS

GFR of preterm individuals catches up with term peers at around three years of age, implying reduced dosages of renally cleared drugs should be considered below this age.

Physiologically Based Pharmacokinetics Modeling in the Neonatal Population-Current Advances, Challenges, and Opportunities

Physiologically based pharmacokinetic (PBPK) modeling is an approach to predicting drug pharmacokinetics, using knowledge of the human physiology involved and drug physiochemical properties. This approach is useful when predicting drug pharmacokinetics in under-studied populations, such as pediatrics. PBPK modeling is a particularly important tool for dose optimization for the neonatal population, given that clinical trials rarely include this patient population. However, important knowledge gaps exist for neonates, resulting in uncertainty with the model predictions. This review aims to outline the sources of variability that should be considered with developing a neonatal PBPK model, the data that are currently available for the neonatal ontogeny, and lastly to highlight the data gaps where further research would be needed.

External validation of population pharmacokinetic models of gentamicin in paediatric population from preterm newborns to adolescents

The aim of this study was to externally validate the predictive performance of published population pharmacokinetic models of gentamicin in all paediatric age groups, from preterm newborns to adolescents. We first selected published population pharmacokinetic models of gentamicin developed in the paediatric population with a wide age range. The parameters of the literature models were then re-estimated using the PRIOR subroutine in NONMEM^®. The predictive ability of the literature and the tweaked models was evaluated. Retrospectively collected data from a routine clinical practice (512 concentrations from 308 patients) were used for validation. The models with covariates characterising developmental changes in clearance and volume of distribution had better predictive performance, which improved further after re-estimation. The tweaked model by Wang 2019 performed best, with suitable accuracy and precision across the complete paediatric population. For patients treated in the intensive care unit, a lower proportion of patients would be expected to reach the target trough concentration at standard dosing. The selected model could be used for model-informed precision dosing in clinical settings where the entire paediatric population is treated. However, for use in clinical practice, the next step should include additional analysis of the impact of intensive care treatment on gentamicin pharmacokinetics, followed by prospective validation.

Current knowledge, challenges and innovations in developmental pharmacology: A combined conect4children Expert Group and European Society for Developmental, Perinatal and Paediatric Pharmacology White Paper

Developmental pharmacology describes the impact of maturation on drug disposition (pharmacokinetics, PK) and drug effects (pharmacodynamics, PD) throughout the paediatric age range. This paper, written by a multidisciplinary group of experts, summarizes current knowledge, and provides suggestions to pharmaceutical companies, regulatory agencies and academicians on how to incorporate the latest knowledge regarding developmental pharmacology and innovative techniques into neonatal and paediatric drug development. Biological aspects of drug absorption, distribution, metabolism and excretion throughout development are summarized. Although this area made enormous progress during the last two decades, remaining knowledge gaps were identified. Minimal risk and burden designs allow for optimally informative but minimally invasive PK sampling, while concomitant profiling of drug metabolites may provide additional insight in the unique PK behaviour in children. Furthermore, developmental PD needs to be considered during drug development, which is illustrated by disease- and/or target organ-specific examples. Identifying and testing PD targets and effects in special populations, and application of age- and/or population-specific assessment tools are discussed. Drug development plans also need to incorporate innovative techniques such as preclinical models to study therapeutic strategies, and shift from sequential enrolment of subgroups, to more rational designs. To stimulate appropriate research plans, illustrations of specific PK/PD-related as well as drug safety-related challenges during drug development are provided. The suggestions made in this joint paper of the Innovative Medicines Initiative conect4children Expert group on Developmental Pharmacology and the European Society for Developmental, Perinatal and Paediatric Pharmacology, should facilitate all those involved in drug development.

Evaluating and Improving Neonatal Gentamicin Pharmacokinetic Models Using Aggregated Routine Clinical Care Data

Model-informed precision dosing (MIPD) can aid dose decision-making for drugs such as gentamicin that have high inter-individual variability, a narrow therapeutic window, and a high risk of exposure-related adverse events. However, MIPD in neonates is challenging due to their dynamic development and maturation and by the need to minimize blood sampling due to low blood volume. Here, we investigate the ability of six published neonatal gentamicin population pharmacokinetic models to predict gentamicin concentrations in routine therapeutic drug monitoring from nine sites in the United State (n = 475 patients). We find that four out of six models predicted with acceptable levels of error and bias for clinical use. These models included known important covariates for gentamicin PK, showed little bias in prediction residuals over covariate ranges, and were developed on patient populations with similar covariate distributions as the one assessed here. These four models were refit using the published parameters as informative Bayesian priors or without priors in a continuous learning process. We find that refit models generally reduce error and bias on a held-out validation data set, but that informative prior use is not uniformly advantageous. Our work informs clinicians implementing MIPD of gentamicin in neonates, as well as pharmacometricians developing or improving PK models for use in MIPD.

Use of Antibiotics in Preterm Newborns

Due to complex maturational and physiological changes that characterize neonates and affect their response to pharmacological treatments, neonatal pharmacology is different from children and adults and deserves particular attention. Although preterms are usually considered part of the neonatal population, they have physiological and pharmacological hallmarks different from full-terms and, therefore, need specific considerations. Antibiotics are widely used among preterms. In fact, during their stay in neonatal intensive care units (NICUs), invasive procedures, including central catheters for parental nutrition and ventilators for respiratory support, are often sources of microbes and require antimicrobial treatments. Unfortunately, the majority of drugs administered to neonates are off-label due to the lack of clinical studies conducted on this special population. In fact, physiological and ethical concerns represent a huge limit in performing pharmacokinetic (PK) studies on these subjects, since they limit the number and volume of blood sampling. Therapeutic drug monitoring (TDM) is a useful tool that allows dose adjustments aiming to fit plasma concentrations within the therapeutic range and to reach specific drug target attainment. In this review of the last ten years’ literature, we performed Pubmed research aiming to summarize the PK aspects for the most used antibiotics in preterms.

Development and application of a pediatric mechanistic kidney model

Pediatric physiologically‐based pharmacokinetic (P‐PBPK) models have been used to predict age related changes in the pharmacokinetics (PKs) of renally cleared drugs mainly in relation to changes in glomerular filtration rate. With emerging data on ontogeny of renal transporters, mechanistic models of renal clearance accounting for the role of active and passive secretion should be developed and evaluated. Data on age‐related physiological changes and ontogeny of renal transporters were applied into a mechanistic kidney within a P‐PBPK model. Plasma concentration–time profile and PK parameters of cimetidine, ciprofloxacin, metformin, tenofovir, and zidovudine were predicted in subjects aged 1 day to 18 years. The predicted and observed plasma concentration–time profiles and PK parameters were compared. The predicted concentration–time profile means and 5th and 95th percent intervals generally captured the observed data and variability in various studies. Overall, based on drugs and age bands, predicted to observed clearance were all within two‐fold and in 11 of 16 cases within 1.5‐fold. Predicted to observed area under the curve (AUC) and maximum plasma concentration (C_max) were within two‐fold in 12 of 14 and 12 of 15 cases, respectively. Predictions in neonates and early infants (up to 14 weeks postnatal age) were reasonable with 15–20 predicted PK parameters within two‐fold of the observed. ciprofloxacin but not zidovudine PK predictions were sensitive to basal kidney uptake transporter ontogeny. The results indicate that a mechanistic kidney model accounting for physiology and ontogeny of renal processes and transporters can predict the PK of renally excreted drugs in children. Further data especially in neonates are required to verify the model and ontogeny profiles.

Implementation of Vancomycin Therapeutic Monitoring Guidelines: Focus on Bayesian Estimation Tools in Neonatal and Pediatric Patients

BACKGROUND

The 2020 consensus guidelines for vancomycin therapeutic monitoring recommend using Bayesian estimation targeting the ratio of the area under the curve over 24 hours to minimum inhibitory concentration as an optimal approach to individualize therapy in pediatric patients. To support institutional guideline implementation in children, the objective of this study was to comprehensively assess and compare published population-based pharmacokinetic (PK) vancomycin models and available Bayesian estimation tools, specific to neonatal and pediatric patients.

METHODS

PubMed and Embase databases were searched from January 1994 to December 2020 for studies in which a vancomycin population PK model was developed to determine clearance and volume of distribution in neonatal and pediatric populations. Available Bayesian software programs were identified and assessed from published articles, software program websites, and direct communication with the software company. In the present review, 14 neonatal and 20 pediatric models were included. Six programs (Adult and Pediatric Kinetics, BestDose, DoseMeRx, InsightRx, MwPharm++, and PrecisePK) were evaluated.

RESULTS

Among neonatal models, Frymoyer et al and Capparelli et al used the largest PK samples to generate their models, which were externally validated. Among the pediatric models, Le et al used the largest sample size, with multiple external validations. Of the Bayesian programs, DoseMeRx, InsightRx, and PrecisePK used clinically validated neonatal and pediatric models.

CONCLUSIONS

To optimize vancomycin use in neonatal and pediatric patients, clinicians should focus on selecting a model that best fits their patient population and use Bayesian estimation tools for therapeutic area under the -curve-targeted dosing and monitoring.

Prediction of glomerular filtration rate maturation across preterm and term neonates and young infants using inulin as marker

Describing glomerular filtration rate (GFR) maturation across the heterogeneous population of preterm and term neonates and infants is important to predict the clearance of renally cleared drugs. This study aims to describe the GFR maturation in (pre)term neonates and young infants (PNA < 90 days) using individual inulin clearance data (CL_inulin). To this end, published GFR maturation models were evaluated by comparing their predicted GFR with CL_inulin retrieved from literature. The best model was subsequently optimized in NONMEM V7.4.3 to better fit the CL_inulin values. Our study evaluated seven models and collected 381 individual CL_inulin values from 333 subjects with median (range) birthweight (BWb) 1880 g (580-4950), gestational age (GA) 34 weeks (25-43), current weight (CW) 1890 g (480-6200), postnatal age (PNA) 3 days (0-75), and CL_inulin 2.20 ml/min (0.43-17.90). The De Cock 2014 model (covariates: BWb and PNA) performed the best in predicting CL_inulin, followed by the Rhodin 2009 model (covariates: CW and postmenstrual age). The final optimized model shows that GFR at birth is determined by BWb, thereafter the maturation rate of GFR is dependent on PNA and GA, with a higher GA showing an overall faster maturation. To conclude, using individual CL_inulin data, we found that a model for neonatal GFR requires a distinction between prenatal maturation quantified by BWb and postnatal maturation. To capture postnatal GFR maturation in (pre)term neonates and young infants, we developed an optimized model in which PNA-related maturation was dependent on GA.

Pre- and Postnatal Maturation are Important for Fentanyl Exposure in Preterm and Term Newborns: A Pooled Population Pharmacokinetic Study

BACKGROUND AND OBJECTIVE

Fentanyl is an opioid commonly used to prevent and treat severe pain in neonates; however, its use is off label and mostly based on bodyweight. Given the limited pharmacokinetic information across the entire neonatal age range, we characterized the pharmacokinetics of fentanyl across preterm and term neonates to individualize dosing.

METHODS

We pooled data from two previous studies on 164 newborns with a median gestational age of 29.0 weeks (range 23.9-42.3), birthweight of 1055 g (range 390-4245), and postnatal age (PNA) of 1 day (range 0-68). In total, 673 plasma samples upon bolus dosing (69 patients; median dose 2.1 μg/kg, median 2 boluses per patient) or continuous infusions (95 patients; median dose 1.1 μg/kg/h for 30 h) with and without boluses were used for population pharmacokinetic modeling in NONMEM^® 7.4.

RESULTS

Clearance in neonates with birthweight of 2000 and 3000 g was 2.8- and 5.0-fold the clearance in a neonate with birthweight of 1000 g, respectively. Fentanyl clearance at PNA of 7, 14, and 21 days was 2.7-fold, 3.8-fold, and 4.6-fold the clearance at 1 day, respectively. Bodyweight-based dosing resulted in large differences in fentanyl concentrations. Depending on PNA and birthweight, fentanyl concentrations increased slowly after the start of therapy for both intermittent boluses and continuous infusion and reached a maximum concentration at 12-48 h.

CONCLUSIONS

As both prenatal and postnatal maturation are important for fentanyl exposure, we propose a birthweight- and PNA-based dosage regimen. To provide rapid analgesia in the first 24 h of treatment, additional loading doses need to be considered.

Pharmacokinetics and Safety of Ceftobiprole in Pediatric Patients

BACKGROUND

Ceftobiprole, the active moiety of the prodrug ceftobiprole medocaril, is an advanced-generation, broad-spectrum, intravenous cephalosporin, which is currently approved for the treatment of adults with hospital-acquired or community-acquired pneumonia.

METHODS

Noncompartmental pharmacokinetics and safety were analyzed from 2 recently completed pediatric studies, a single-dose, phase 1 study in neonates and infants up to 3 months of age (7.5 mg/kg) and a phase 3 study in patients 3 months to 17 years of age with pneumonia (10-20 mg/kg with a maximum of 500 mg per dose every 8 hours for up to 14 days).

RESULTS

Total ceftobiprole plasma concentrations peaked at the end of infusion. Half life (median ranging from 1.9 to 2.9 hours) and overall exposure (median AUC ranging from 66.6 to 173 μg•h/mL) were similar to those in adults (mean ± SD, 3.3 ± 0.3 hours and 102 ± 11.9 μg•h/mL, respectively). Calculated free-ceftobiprole concentrations in the single-dose study remained above a minimum inhibitory concentration (MIC) of 4 mg/L (fT > MIC of 4 mg/L) for a mean of 5.29 hours after dosing. In the pneumonia study, mean fT > MIC of 4 mg/L was ≥5.28 hours in all dose groups. Ceftobiprole was well tolerated in both studies.

CONCLUSIONS

Pharmacokinetic parameters of ceftobiprole characterized in the pediatric population were within the range of those observed in adults. In the pneumonia study, the lowest percentage of the dosing interval with fT > MIC of 4 mg/L was 50.8%, which suggests that pharmacokinetic-pharmacodynamic target attainment can be sufficient in pediatric patients. Ceftobiprole was well tolerated.

An Update on Population Pharmacokinetic Analyses of Vancomycin, Part II: In Pediatric Patients

Vancomycin is widely used in pediatric patients, however, large inter- and intraindividual variability are observed in vancomycin pharmacokinetics, affecting proper therapeutic monitoring. This review aimed at providing a comprehensive synthesis of the population pharmacokinetic models of vancomycin in pediatric patients and identifying potential factors responsible for the variability observed in various subpopulations. We conducted a literature search of the PubMed and EMBASE databases to obtain population pharmacokinetic studies for vancomycin published between January 2011 and January 2020, which resulted in a total of 33 studies. Vancomycin pharmacokinetics were generally characterized using a one-compartment model (n = 27), while a two-compartment model was used in six studies. The median (interquartile range) of the typical vancomycin clearance (CL) and the total volume of distribution adjusted to the median or mean body weight of the respective study was 0.103 L/h/kg (0.071-0.125) and 0.64 L/kg (0.59-1.03), respectively. Median weight-adjusted CL between different child age groups, such as infants and adolescents, did not appear to vary significantly, although the sample size for many age groups was very small. Examples of the conditions with relatively abnormal vancomycin pharmacokinetic values include renal insufficiency, sepsis, hematological and solid malignancy, and hypothermia treatment. Factors influencing pediatric vancomycin pharmacokinetics after adjusting for size and maturation include various renal function descriptors and some case-specific variables such as dialysate flow rate, ultrafiltrate output, and hypothermia. This review was able to document possible variables explaining the high variability observed in certain subpopulations and contrast vancomycin pharmacokinetics in different pediatric subpopulations.

Pharmacokinetics of Antibiotics in Pediatric Intensive Care: Fostering Variability to Attain Precision Medicine

Children show important developmental and maturational changes, which may contribute greatly to pharmacokinetic (PK) variability observed in pediatric patients. These PK alterations are further enhanced by disease-related, non-maturational factors. Specific to the intensive care setting, such factors include critical illness, inflammatory status, augmented renal clearance (ARC), as well as therapeutic interventions (e.g., extracorporeal organ support systems or whole-body hypothermia [WBH]). This narrative review illustrates the relevance of both maturational and non-maturational changes in absorption, distribution, metabolism, and excretion (ADME) applied to antibiotics. It hereby provides a focused assessment of the available literature on the impact of critical illness—in general, and in specific subpopulations (ARC, extracorporeal organ support systems, WBH)—on PK and potential underexposure in children and neonates. Overall, literature discussing antibiotic PK alterations in pediatric intensive care is scarce. Most studies describe antibiotics commonly monitored in clinical practice such as vancomycin and aminoglycosides. Because of the large PK variability, therapeutic drug monitoring, further extended to other antibiotics, and integration of model-informed precision dosing in clinical practice are suggested to optimise antibiotic dose and exposure in each newborn, infant, or child during intensive care.

Optimizing Ceftobiprole Dosage in Pediatric Patients: A Model-Based Approach

Ceftobiprole is an advanced-generation cephalosporin for intravenous administration with activity against Gram-positive and Gram-negative organisms. A population pharmacokinetic (PK) model characterizing the disposition of ceftobiprole in plasma using data from patients in three pediatric studies was developed. Model-based simulations were subsequently performed to assist in dose optimization for the treatment of pediatric patients with hospital-acquired or community-acquired pneumonia. The population PK data set comprised 518 ceftobiprole plasma concentrations from 107 patients from 0 (birth) to 17 years of age. Ceftobiprole PK was well described by a three-compartment model with linear elimination. Ceftobiprole clearance was modeled as a function of glomerular filtration rate; other PK parameters were scaled to body weight. The final population PK model provided a robust and reliable description of the PK of ceftobiprole in the pediatric study population. Model-based simulations using the final model suggested that a ceftobiprole dose of 15 mg/kg of body weight infused over 2 h and administered every 12 h in neonates and infants <3 months of age or every 8 h in older pediatric patients would result in a ceftobiprole exposure consistent with that in adults and good pharmacokinetic-pharmacodynamic target attainment. The dose should be reduced to 10 mg/kg every 12 h in neonates and infants <3 months of age who weigh <4 kg to avoid high exposures. Extended intervals and reduced doses may be required for pediatric patients older than 3 months of age with renal impairment.

Optimizing gentamicin dosing in different pediatric age groups using population pharmacokinetics and Monte Carlo simulation

Introduction

The use of once daily dosing of aminoglycosides in pediatrics is increasing but studies on dose optimization targeting the pediatric population are limited. This study aimed to derive a population pharmacokinetic model of gentamicin and apply it to design optimal dosing regimens in pediatrics.

Methods

Population pharmacokinetics of gentamicin in pediatrics was described from a retrospective chart review of plasma gentamicin concentration data (peak/ trough levels) of pediatric patients (1 month − 12 years), admitted to non-critically ill pediatrics. Monte Carlo simulations were performed on the resulting pharmacokinetic model to assess the probability of achieving a C_max/MIC target of 10 mg/L over a range of gentamicin MICs of 0.5–2 mg/L and once daily gentamicin dosing regimens. Results: A two-compartment model with additive residual error best described the model with weight incorporated as a significant covariate for both clearance and volume of distribution. Monte Carlo simulations demonstrated a good probability of target attainment even at a MIC of 2 mg/L, where neonates required doses of 6-7 mg/kg/day and older pediatrics required lower daily doses of 4–5 mg/kg/day while maintaining trough gentamicin concentration below the toxicity limit of 1 mg/L. Conclusion: Once daily dosing is a reasonable option in pediatrics that allows target attainment while maintaining trough gentamicin level below the limits of toxicity.

Population Pharmacokinetic Models of Vancomycin in Paediatric Patients: A Systematic Review

BACKGROUND

Vancomycin is commonly used to treat gram-positive bacterial infections in the paediatric population, but dosing can be challenging. Population pharmacokinetic (popPK) modelling can improve individualization of dosing regimens. The primary objective of this study was to describe popPK models of vancomycin and factors that influence pharmacokinetic (PK) variability in paediatric patients.

METHODS

Systematic searches were conducted in the Cochrane Central Register of Controlled Trials, MEDLINE, EMBASE, International Pharmaceutical Abstracts and the grey literature without language or publication status restrictions from inception to 17 August 2020. Observational studies that described the development of popPK models of vancomycin in paediatric patients (< 18 years of age) were included. Risk of bias was assessed using the National Heart, Lung and Blood Institute Study Quality Assessment Tool for Case Series Studies.

RESULTS

Sixty-four observational studies (1 randomized controlled trial, 13 prospective studies and 50 retrospective studies of 9019 patients with at least 25,769 serum vancomycin concentrations) were included. The mean age was 2.5 years (range 1 day-18 years), serum creatinine was 47.1 ± 33.6 µmol/L, and estimated creatinine clearance was 97.4 ± 76 mL/min/1.73m². Most studies found that vancomycin PK was best described by a one-compartment model (71.9%). There was a wide range of clearance and volume of distribution (V_d) values (range 0.014-0.27 L/kg/h and 0.43-1.46 L/kg, respectively) with interindividual variability as high as 49.7% for clearance and 136% for V_d, proportional residual variability up to 37.5% and additive residual variability up to 17.5 mg/L. The most significant covariates for clearance were weight, age, and serum creatinine or creatinine clearance, and weight for V_d. Variable dosing recommendations were suggested.

CONCLUSION

Numerous popPK models of vancomycin were derived, however external validation of suggested dosing regimens and analyses in subgroup paediatric populations such as dialysis patients are still needed before a popPK model with best predictive performance can be applied for dosing recommendations. Significant intraindividual and interindividual PK variability was present, which demonstrated the need for ongoing therapeutic drug monitoring and derivation of PK models for vancomycin for certain subgroup populations, such as dialysis patients.

Estimation of Ontogeny Functions for Renal Transporters Using a Combined Population Pharmacokinetic and Physiology-Based Pharmacokinetic Approach: Application to OAT1,3

To date, information on the ontogeny of renal transporters is limited. Here, we propose to estimate the in vivo functional ontogeny of transporters using a combined population pharmacokinetic (popPK) and physiology-based pharmacokinetic (PBPK) modeling approach called popPBPK. Clavulanic acid and amoxicillin were used as probes for glomerular filtration, combined glomerular filtration, and active secretion through OAT1,3, respectively. The predictive value of the estimated OAT1,3 ontogeny function was assessed by PBPK predictions of renal clearance (CL_R) of other OAT1,3 substrates: cefazolin and piperacillin. Individual CL_Rpost-hoc values, obtained from a published popPK model on the concomitant use of clavulanic acid and amoxicillin in critically ill children between 1 month and 15 years, were used as dependent variables in the popPBPK analysis. CL_R was re-parameterized according to PBPK principles, resulting in the estimation of OAT1,3-mediated intrinsic clearance (CL_{int,OAT1,3,invivo}) and its ontogeny. CL_{int,OAT1,3,invivo} ontogeny was described by a sigmoidal function, reaching half of adult level around 7 months of age, comparable to findings based on renal transporter-specific protein expression data. PBPK-based CL_R predictions including this ontogeny function were reasonably accurate for piperacillin in a similar age range (2.5 months–15 years) as well as for cefazolin in neonates as compared to published data (%RMSPE of 21.2 and 22.8%, respectively and %PE within ±50%). Using this novel approach, we estimated an in vivo functional ontogeny profile for CL_{int,OAT1,3,invivo} that yields accurate CL_R predictions for different OAT1,3 substrates across different ages. This approach deserves further study on functional ontogeny of other transporters.

Dosing Recommendations for Vancomycin in Children and Adolescents with Varying Levels of Obesity and Renal Dysfunction: a Population Pharmacokinetic Study in 1892 Children Aged 1-18 Years

Vancomycin is an effective but potentially nephrotoxic antibiotic commonly used for severe infections. Dosing guidelines for vancomycin in obese children and adolescents with or without renal impairment are currently lacking. This study describes the pharmacokinetics of vancomycin in a large pediatric cohort with varying degrees of obesity and renal function to design practical dosing guidelines for this population. A multi-center retrospective population pharmacokinetic study was conducted using data from patients aged 1−18 years who received >1 dose of vancomycin and had ≥1 vancomycin concentration measured between January 2006 and December 2012. Besides pharmacokinetic data, age, gender, body weight, creatinine clearance (CL_cr, bedside Schwartz equation), ward, race, and neutropenic status were collected. Population pharmacokinetic analysis and simulations were performed using NONMEM7.4. A total of 1892 patients (5524 samples) were included, with total body weight (TBW) ranging 6−188 kg (1344 normal weight, 247 overweight, and 301 obese patients) and CL_cr down to 8.6 mL/min/1.73 m². The two-compartment model, with clearance (CL) significantly increasing with TBW and CL_cr, central and peripheral volume of distribution and inter-compartmental clearance increasing with TBW, performed well for all age, weight, and renal function ranges. A dosing guideline is proposed that integrates body weight and CL_cr resulting in effective and safe exposures across all ages, body weight, and renal functions in the pediatric population. We have characterized the full pharmacokinetic profile of vancomycin in obese children and adolescents aged 1−18 years and propose a practical dosing guideline that integrates both body weight and renal function.

Recent advances in the ontogeny of drug disposition

Developmental changes that occur throughout childhood have long been known to impact drug disposition. However, pharmacokinetic studies in the paediatric population have historically been limited due to ethical concerns arising from incorporating children into clinical trials. As such, much of the early work in the field of developmental pharmacology was reliant on difficult-to-interpret in vitro and in vivo animal studies. Over the last 2 decades, our understanding of the mechanistic processes underlying age-related changes in drug disposition has advanced considerably. Progress has largely been driven by technological advances in mass spectrometry-based methods for quantifying proteins implicated in drug disposition, and in silico tools that leverage these data to predict age-related changes in pharmacokinetics. This review summarizes our current understanding of the impact of childhood development on drug disposition, particularly focusing on research of the past 20 years, but also highlighting select examples of earlier foundational research. Equally important to the studies reviewed herein are the areas that we cannot currently describe due to the lack of research evidence; these gaps provide a map of drug disposition pathways for which developmental trends still need to be characterized.

Simplified Dosing Regimens for Gentamicin in Neonatal Sepsis

Background: The effectiveness of antibiotics for the treatment of severe bacterial infections in newborns in resource-limited settings has been determined by empirical evidence. However, such an approach does not warrant optimal exposure to antibiotic agents, which are known to show different disposition characteristics in this population. Here we evaluate the rationale for a simplified regimen of gentamicin taking into account the effect of body size and organ maturation on pharmacokinetics. The analysis is supported by efficacy data from a series of clinical trials in this population.

Methods: A previously published pharmacokinetic model was used to simulate gentamicin concentration vs. time profiles in a virtual cohort of neonates. Model predictive performance was assessed by supplementary external validation procedures using therapeutic drug monitoring data collected in neonates and young infants with or without sepsis. Subsequently, clinical trial simulations were performed to characterize the exposure to intra-muscular gentamicin after a q.d. regimen. The selection of a simplified regimen was based on peak and trough drug levels during the course of treatment.

Results: In contrast to current World Health Organization guidelines, which recommend gentamicin doses between 5 and 7.5 mg/kg, our analysis shows that gentamicin can be used as a fixed dose regimen according to three weight-bands: 10 mg for patients with body weight <2.5 kg, 16 mg for patients with body weight between 2.5 and 4 kg, and 30 mg for those with body weight >4 kg.

Conclusion: The choice of the dose of an antibiotic must be supported by a strong scientific rationale, taking into account the differences in drug disposition in the target patient population. Our analysis reveals that a simplified regimen is feasible and could be used in resource-limited settings for the treatment of sepsis in neonates and young infants with sepsis aged 0–59 days.

The Neonatal and Juvenile Pig in Pediatric Drug Discovery and Development

Pharmacotherapy in pediatric patients is challenging in view of the maturation of organ systems and processes that affect pharmacokinetics and pharmacodynamics. Especially for the youngest age groups and for pediatric-only indications, neonatal and juvenile animal models can be useful to assess drug safety and to better understand the mechanisms of diseases or conditions. In this respect, the use of neonatal and juvenile pigs in the field of pediatric drug discovery and development is promising, although still limited at this point. This review summarizes the comparative postnatal development of pigs and humans and discusses the advantages of the juvenile pig in view of developmental pharmacology, pediatric diseases, drug discovery and drug safety testing. Furthermore, limitations and unexplored aspects of this large animal model are covered. At this point in time, the potential of the neonatal and juvenile pig as nonclinical safety models for pediatric drug development is underexplored.

Population Pharmacokinetics of Unbound and Total Teicoplanin in Critically Ill Pediatric Patients

BACKGROUND AND OBJECTIVES

Teicoplanin is a highly protein-bound antibiotic, increasingly used to treat serious Gram-positive infections in critically ill children. Maturational and pathophysiological intensive care unit-related changes often lead to altered pharmacokinetics. In this study, the objectives were to develop a pediatric population-pharmacokinetic model of unbound and total teicoplanin concentrations, to investigate the impact of plasma albumin levels and renal function on teicoplanin pharmacokinetics, and to evaluate the efficacy of the current weight-based dosing regimen.

METHODS

An observational pharmacokinetic study was performed and blood samples were collected for quantification of unbound and total concentrations of teicoplanin after the first dose and in assumed steady-state conditions. A population-pharmacokinetic analysis was conducted using a standard sequential approach and Monte Carlo simulations were performed for a probability of target attainment analysis using previously published pharmacokinetic-pharmacodynamic targets.

RESULTS

A two-compartment model with allometric scaling of pharmacokinetic parameters and non-linear plasma protein binding best described the data. Neither the inclusion of albumin nor the renal function significantly improved the model and no other covariates were supported for inclusion in the final model. The probability of target attainment analysis showed that the standard dosing regimen does not satisfactory attain the majority of the proposed targets.

CONCLUSIONS

We successfully characterized the pharmacokinetics of unbound and total teicoplanin in critically ill pediatric patients. The highly variable unbound fraction of teicoplanin could not be predicted using albumin levels, which may support the use of therapeutic drug monitoring of unbound concentrations. Poor target attainment was shown for the most commonly used dosing regimen, regardless of the pharmacokinetic-pharmacodynamic target evaluated.

The Predictive Value of Glomerular Filtration Rate-Based Scaling of Pediatric Clearance and Doses for Drugs Eliminated by Glomerular Filtration with Varying Protein-Binding Properties

INTRODUCTION

For drugs eliminated by glomerular filtration (GF), clearance (CL) is determined by GF rate (GFR) and the unbound fraction of the drug. When predicting CL of GF-eliminated drugs in children, instead of physiologically based pharmacokinetic (PBPK) methods that consider changes in both GFR and protein binding, empiric bodyweight-based methods are often used. In this article, we explore the predictive value of scaling using a GFR function, and compare the results with linear and allometric scaling methods for drugs with different protein-binding properties.

METHODS

First, different GFR maturation functions were compared to identify the GFR function that would yield the most accurate GFR predictions across the pediatric age range compared with published pediatric inulin/mannitol CL values. Subsequently, the accuracy of pediatric CL scaling using this GFR maturation function was assessed and compared with PBPK CL predictions for hypothetical drugs binding, to varying extents, to serum albumin or α-acid glycoprotein across the pediatric age range. Additionally, empiric bodyweight-based methods were assessed.

RESULTS

The published GFR maturation functions yielded comparable maturation profiles, with the function reported by Salem et al. leading to the most accurate predictions. On the basis of this function, GFR-based scaling yields reasonably accurate (percentage prediction error ≤ 50%) pediatric CL values for all drugs, except for some drugs highly bound to AGP in neonates. Overall, this method was more accurate than linear or 0.75 allometric bodyweight-based scaling.

CONCLUSION

When scaling CL and dose by GFR function, maturational changes in plasma protein concentrations impact GF minimally, making this method a superior alternative to empiric bodyweight-based scaling.

The role of antibiotic pharmacokinetic studies performed post-licensing

Post-licensing pharmacometric studies can provide a better understanding of the pharmacokinetic (PK) alterations in special patient populations and may lead to better clinical outcomes. Some patient populations exhibit markedly different pathophysiology to general ward patients or healthy individuals. This may be developmental (paediatric patients), a manifestation of an underlying disease pathology (patients with obesity or haematological malignancies) or due to medical interventions (critically ill patients receiving extracorporeal therapies). This paper outlines the factors that affect the PK of special patient populations and describes some novel methods of antimicrobial administration that may increase antimicrobial concentrations at the site of infection and improve treatment of severe infection.

Population pharmacokinetic meta-analysis of individual data to design the first randomized efficacy trial of vancomycin in neonates and young infants

OBJECTIVES

In the absence of consensus, the present meta-analysis was performed to determine an optimal dosing regimen of vancomycin for neonates.

METHODS

A 'meta-model' with 4894 concentrations from 1631 neonates was built using NONMEM, and Monte Carlo simulations were performed to design an optimal intermittent infusion, aiming to reach a target AUC0-24 of 400 mg·h/L at steady-state in at least 80% of neonates.

RESULTS

A two-compartment model best fitted the data. Current weight, postmenstrual age (PMA) and serum creatinine were the significant covariates for CL. After model validation, simulations showed that a loading dose (25 mg/kg) and a maintenance dose (15 mg/kg q12h if <35 weeks PMA and 15 mg/kg q8h if ≥35 weeks PMA) achieved the AUC0-24 target earlier than a standard 'Blue Book' dosage regimen in >89% of the treated patients.

CONCLUSIONS

The results of a population meta-analysis of vancomycin data have been used to develop a new dosing regimen for neonatal use and to assist in the design of the model-based, multinational European trial, NeoVanc.

Does "Birth" as an Event Impact Maturation Trajectory of Renal Clearance via Glomerular Filtration? Reexamining Data in Preterm and Full-Term Neonates by Avoiding the Creatinine Bias

Glomerular filtration rate (GFR) is an important measure of renal function. Various models for its maturation have recently been compared; however, these have used markers, which are subject to different renal elimination processes. Inulin clearance data (a purer probe of GFR) collected from the literature were used to determine age‐related changes in GFR aspects of renal drug excretion in pediatrics. An ontogeny model was derived using a best‐fit model with various combinations of covariates such as postnatal age, gestational age at birth, and body weight. The model was applied to the prediction of systemic clearance of amikacin, gentamicin, vancomycin, and gadobutrol. During neonatal life, GFR increased as a function of both gestational age at birth and postnatal age, hence implying an impact of birth and a discrepancy in GFR for neonates with the same postmenstrual age depending on gestational age at birth (ie, neonates who were outside the womb longer had higher GFR, on average). The difference in GFR between pre‐term and full‐term neonates with the same postmenstrual age was negligible from beyond 1.25 years. Considering both postnatal age and gestational age at birth in GFR ontogeny models is important because postmenstrual age alone ignores the impact of birth. Most GFR models use covariates of body size in addition to age. Therefore, prediction from these models will also depend on the change in anthropometric characteristics with age. The latter may not be similar in various ethnic groups, and this makes the head‐to‐head comparison of models very challenging.

The Utility of Pharmacometric Models in Clinical Pharmacology Research in Infants

Purpose of commentary

Acquiring knowledge on drug disposition and action in infant is challenging because of the problem of sparse and unbalanced data obtained for each individual infant due to the limited blood volume as well as the issue of extensive inter-subject and intra-subject variability in drug exposure and response due to the fast growth and dynamic maturation changes in infants. This commentary highlights the importance of using population-based pharmacometric models to improve knowledge on drug disposition and action in infants.

Recent findings

Pharmacometric modeling remains to be critical in clinical pharmacology research in infants. Many pediatric covariate models developed for scaling of drug clearance use a combination of allometric weight scaling to account for size change and a sigmoid function of antenatal development and postnatal maturation to characterize the age-related maturation. To expedite the development of safe and effective dosing regimens in infants, a number of strategies have been proposed recently, including the use of pediatric covariate model obtained from one drug for extrapolation to other drugs undergoing similar elimination pathways, as well as the combination of opportunistic clinical studies and population-based pharmacometrics models.

Summary

Population-based pharmacometric modeling plays a pivotal role in clinical pharmacology research in infants. Most of the covariate models reported so far focus on antibiotics undergoing renal elimination. Novel modeling strategies have been proposed recently to facilitate clinical pharmacology research and expedite the dose optimization process in infants.

Toward precision medicine in pediatric population using cytochrome P450 phenotyping approaches and physiologically based pharmacokinetic modeling

The activity of drug-metabolizing enzymes (DME) shows high inter- and intra-individual variability. Genetic polymorphisms, exposure to drugs, and environmental toxins are known to significantly alter DME expression. In addition, the activity of these enzymes is highly age-dependent due to maturation processes that occur during development. Currently, there is a vast choice of phenotyping methods in adults using exogenous probes to characterize the activity of these enzymes. However, this can hardly be applied to children since it requires the intake of non-therapeutic xenobiotics. In addition, sampling may be challenging in the pediatric population for a variety of reasons: limited volume (e.g., blood), inappropriate sampling methods for age (e.g., urine), and metric requiring invasive or multiple blood samples. This review covers the main existing methods that can be used in the pediatric population to determine DME activity, with a particular focus on cytochrome P450 enzymes. Less invasive tools are described, including phenotyping using endogenous probes. Finally, the potential of pediatric model-informed precision dosing using physiologically based pharmacokinetic modeling is discussed.

Augmented renal clearance in pediatric intensive care: are we undertreating our sickest patients?

Many critically ill patients display a supraphysiological renal function with enhanced renal perfusion and glomerular hyperfiltration. This phenomenon described as augmented renal clearance (ARC) may result in enhanced drug elimination through renal excretion mechanisms. Augmented renal clearance seems to be triggered by systemic inflammation and therapeutic interventions in intensive care. There is growing evidence that ARC is not restricted to the adult intensive care population, but is also prevalent in critically ill children. Augmented renal clearance is often overlooked due to the lack of reliable methods to assess renal function in critically ill children. Standard equations to calculate glomerular filtration rate (GFR) are developed for patients who have a steady-state creatinine production and a stable renal function. Those formulas are not reliable in critically ill patients with acutely changing GFR and tend to underestimate true GFR in patients with ARC. Tools for real-time, continuous, and non-invasive measurement of fluctuating GFR are most needed to identify changes in kidney function during critical illness and therapeutic interventions. Such devices are currently being validated and hold a strong potential to become the standard of practice. In the meantime, urinary creatinine clearance is considered the most reliable method to detect ARC in critically ill patients. Augmented renal clearance is clearly associated with subtherapeutic antimicrobial concentrations and subsequent therapeutic failure. This warrants the need for adjusted dosing regimens to optimize pharmacokinetic and pharmacodynamic target attainment. This review aims to summarize current knowledge on ARC in critically ill children, to give insight into its possible pathophysiological mechanism, to evaluate screening methods for ARC in the pediatric intensive care population, and to illustrate the effect of ARC on drug exposure, therapeutic efficacy, and clinical outcome.

External Validation of Model-Based Dosing Guidelines for Vancomycin, Gentamicin, and Tobramycin in Critically Ill Neonates and Children: A Pragmatic Two-Center Study

BACKGROUND

The Dutch Pediatric Formulary (DPF) increasingly bases its guidelines on model-based dosing simulations from pharmacokinetic studies. This resulted in nationwide dose changes for vancomycin, gentamicin, and tobramycin in 2015.

OBJECTIVE

We aimed to evaluate target attainment of these altered, model-based doses in critically ill neonates and children.

METHODS

This was a retrospective cohort study in neonatal intensive care unit (NICU) and pediatric ICU (PICU) patients receiving vancomycin, gentamicin, or tobramycin between January 2015 and March 2017 in two university hospitals. The first therapeutic drug monitoring concentration for each patient was collected, as was clinical and dosing information. Vancomycin and tobramycin target trough concentrations were 10-15 and ≤ 1 mg/L, respectively. Target gentamicin trough and peak concentrations were < 1 and 8-12 mg/L, respectively.

RESULTS

In total, 482 patients were included (vancomycin [PICU] n = 62, [NICU] n = 102; gentamicin [NICU] n = 97; tobramycin [NICU] n = 221). Overall, median trough concentrations were within the target range for all cohorts but showed large interindividual variability, causing nontarget attainment. Trough concentrations were outside the target range in 66.1%, 60.8%, 14.7%, and 23.1% of patients in these four cohorts, respectively. Gentamicin peak concentrations were outside the range in 69% of NICU patients (term neonates 87.1%, preterm infants 57.1%). Higher creatinine concentrations were associated with higher vancomycin and tobramycin trough concentrations.

CONCLUSION

This study illustrates the need to validate model-based dosing advice in the real-world setting as both sub- and supratherapeutic concentrations of vancomycin, gentamicin, and tobramycin were very prevalent. Our data underline the necessity for further individualization by addressing the high interindividual variability to improve target attainment.

Optimizing Antibiotic Drug Therapy in Pediatrics: Current State and Future Needs

The selection of the right antibiotic and right dose necessitates clinicians understand the contribution of pharmacokinetic variability stemming from age-related physiologic maturation and the pharmacodynamics to optimize drug exposure for clinical response. The complexity of selecting the right dose arises from the multiplicity of pediatric age groups, from premature neonates to adolescents. Body size and age (which relate to organ function) must be incorporated to optimize antibiotic dosing in this vulnerable population. In the effort to optimize and individualize drug dosing regimens, clinical pharmacometrics that incorporate population-based pharmacokinetic modeling, Bayesian estimation, and Monte Carlo simulations are utilized as a quantitative approach to understanding and predicting the pharmacology and clinical and microbiologic efficacy of antibiotics. In addition, opportunistic study designs and alternative blood sampling strategies can serve as practical approaches to ensure successful conduct of pediatric studies. This review article examines relevant literature on optimization of antibiotic pharmacotherapy in pediatric populations published within the last decade. Specific pediatric antibiotic data, including beta-lactam antibiotics, aminoglycosides, and vancomycin, are critically evaluated.

Pharmacokinetic modelling and Bayesian estimation-assisted decision tools to optimize vancomycin dosage in neonates: only one piece of the puzzle

Introduction: Vancomycin is commonly administered to neonates, while observational data on therapeutic drug monitoring (TDM, trough levels) suggest that vancomycin exposure and dosage remain substandard. Area covered: Data on vancomycin pharmacokinetics (PK) and its covariates are abundant. Consequently, modeling is an obvious tool to improve targeted exposure, with a shift from TDM trough levels to area under the curve (AUC_24h) targets, as in adults. Continuous administration appeared as a practice to facilitate AUC_24h target attainment, while Bayesian model-supported targeting emerged as a novel tool. However, the AUC_24h/MIC (minimal inhibitory concentration) target itself should consider neonate-specific aspects (bloodstream infections, coagulase-negative staphylococci, protein binding, underexplored causes of variability, like assays, preparation and administration inaccuracies, or missing covariates). Expert opinion: To improve targeted exposure in neonates, initial vancomycin prescription should be based on 'a priori model-based individual dosing' using validated dosing regimens, followed by further tailoring by dosing optimization applying Bayesian estimation-assisted TDM. Future research should focus on the feasibility to integrate these tools (individualized dosing, Bayesian models) in clinical practice, and to perform PK/PD studies in the relevant animal models and human neonatal setting (coagulase-negative staphylococci, bloodstream infections).

Population Pharmacokinetic Modeling of Gentamicin in Pediatrics

The primary objective of this work was to characterize the pharmacokinetics (PK) of gentamicin across the whole pediatric age spectrum from premature neonates to young adults with a single model by identifying significant clinical predictors. A nonlinear mixed-effect population PK model was developed with retrospective therapeutic drug-monitoring data. A total of 6459 drug concentration measurements from 3370 hospitalized patients were collected for model building (n = 2357) and evaluation (n = 1013). In agreement with previously reported models, a 2-compartment model with first-order elimination best described the drug PK. Patient-specific factors significantly impacting gentamicin clearance included fat-free mass, postmenstrual age, and serum creatinine (SCr). Based on our model, the deviation of the individual SCr from the age-dependent expected mean SCr value (SCrM) can result in a 40% lower clearance in a patient with renal impairment than that in a patient with normal kidney function, with SCrM:SCr ratios between 0.16 and 3.2 in this study. Consistent with the known age-dependent changes of the proportion of extracellular water in body weight, the inclusion of the impact of extracellular water maturation on the central volume of distribution was found to improve the model fitting significantly. In comparison with other published models, model evaluation suggested the developed model was the least biased and physiologically most representative. These results will be used to inform individualized initial dosing strategies and serve as a prior PK model for Bayesian updating and forecasting as individual clinical observations become available.

Tobramycin Clearance Is Best Described by Renal Function Estimates in Obese and Non-obese Individuals: Results of a Prospective Rich Sampling Pharmacokinetic Study

Purpose

Tobramycin is an aminoglycoside antibiotic of which the 24 h exposure correlates with efficacy. Recently, we found that clearance of the aminoglycoside gentamicin correlates with total body weight (TBW). In this study, we investigate the full pharmacokinetic profile of tobramycin in obese and non-obese individuals with normal renal function.

Methods

Morbidly obese individuals (n = 20) undergoing bariatric surgery and non-obese healthy volunteers (n = 8), with TBW ranging 57–194 kg, received an IV dose of tobramycin with plasma concentrations measured over 24 h (n = 10 per individual). Statistical analysis, modelling and simulations were performed using NONMEM.

Results

In a two-compartment model, TBW was the best predictor for central volume of distribution (p < 0.001). For clearance, MDRD (de-indexed for body surface area) was identified as best covariate (p < 0.001), and was superior over TBW ((p < 0.05). Other renal function estimates (24 h urine GFR and de-indexed CKD-EPI) led to similar results as MDRD (all p < 0.001)).

Conclusions

In obese and non-obese individuals with normal renal function, renal function estimates such as MDRD were identified as best predictors for tobramycin clearance, which may imply that other processes are involved in clearance of tobramycin versus gentamicin. To ensure similar exposure across body weights, we propose a MDRD-based dosing nomogram for obese patients.

Electronic supplementary material

The online version of this article (10.1007/s11095-019-2651-2) contains supplementary material, which is available to authorized users.

A review of population pharmacokinetic models of gentamicin in paediatric patients

WHAT IS KNOWN AND OBJECTIVES

Gentamicin is often used for the treatment of Gram-negative infections. Due to pharmacokinetic variability in paediatric patients, appropriate dosing of gentamicin in the paediatric population is challenging. This article reviews published population pharmacokinetic models of gentamicin in paediatric patients, identifies covariates that significantly influence gentamicin pharmacokinetics, and determines whether there is a consensus on proposed dosing for intravenous gentamicin in this population.

METHODS

The PubMed database was searched for articles published until the end of 2017. If the articles described population pharmacokinetic models of gentamicin in the paediatric population (after intravenous administration of gentamicin), the following data were extracted: type of study, year of publication, population characteristics and number of patients, gentamicin dosing, total number of gentamicin (serum and/or plasma) concentrations, type of population modelling approach, developed model with pharmacokinetic parameters and covariates included.

RESULTS AND DISCUSSION

In most of the studies, one- or two-compartment modelling was applied. The mean estimated gentamicin clearance for newborns, infants and the complete paediatric population was 0.048, 0.13 and 0.067 L/h/kg, respectively, and the mean predicted volume of distribution was 0.475, 0.35 and 0.33 L/kg, respectively. The values reflect differences in body composition and kidney maturation within the different paediatric populations. Gentamicin pharmacokinetics were most influenced by age, body size and renal function.

WHAT IS NEW AND CONCLUSION

Based on our review, the authors agree on a prolonged dosing interval for preterm and term newborns (up to 48 hours). However, there was no agreement on proposed dosing with respect to gestational age. In general, the proposed daily doses were lower compared to those initially applied for preterm newborns and comparable to those for term newborns. For infants and children, the dosing interval remained unchanged (24 hours), but the proposed daily doses were higher than actually applied. When differences in the paediatric population are considered and an appropriate population PK model with applicable covariates is applied, dosing can be individualized. In the future, studies of gentamicin pharmacokinetics in paediatric patients should focus on currently underestimated covariates, such as fat-free mass, concomitantly administered drugs, body temperature and critical illness because these can change gentamicin PK considerably. Consequently, different dosing is required and TDM becomes even more important.

Optimizing Gentamicin Dosing in Pediatrics Using Monte Carlo Simulations

Gentamicin is known to have concentration-dependent bactericidal activity, and its nephrotoxic effect is well described. We developed a population pharmacokinetic/pharmacodynamic model to optimize gentamicin dosing in pediatrics. Data were retrospectively collected for pediatric patients 1 month to 12 years of age, admitted to general pediatric wards or intensive care units and received gentamicin for suspected or proven Gram-negative infections at King Saud University Medical City, Riyadh, Saudi Arabia. A total of 306 gentamicin peak and trough concentrations sets from 107 patients were analyzed with mean (±standard deviation) patient age and weight of 4.5 ± 3.5 years and 16.7 ± 10.8 kg, respectively. Gentamicin pharmacokinetics were adequately described with a one compartment system (R = 0.82, bias = 1.75% and precision = 88% for population predictions and R = 0.94, bias = 5% and precision = 29% for individual predictions). The gentamicin pharmacokinetic parameters were as follows: volume of distribution = 8.9 L, total body clearance = 2.8 L/h for a 20-kg patient. Monte Carlo simulations showed that doses of 5-6 mg/kg/dose once daily are adequate only to treat infections with Gram-negative organisms having minimal inhibitory concentration less than 1 µg/mL. While, at minimal inhibitory concentration of 1 µg/mL, higher doses (7-8 mg/kg/dose once daily) are needed to maximize the efficacy of gentamicin. However, at minimal inhibitory concentration of 2 µg/mL, even a 10 mg/kg dose showed poor target attainment (52%). The finding of this study highlights the need to reevaluate the current breakpoints of gentamicin and also to assess the safety of higher doses of gentamicin in pediatrics.

Vancomycin pharmacokinetics in critically ill neonates receiving extracorporeal membrane oxygenation

Objective

This study aimed to characterise vancomycin pharmacokinetics in critically ill neonates undergoing extracorporeal membrane oxygenation.

Methods

In a retrospective analysis, the pharmacokinetics of vancomycin were determined in 25 full-term neonates receiving extracorporeal membrane oxygenation and compared with those of matched controls (n = 25) not receiving extracorporeal membrane oxygenation.

Results

The half-life of vancomycin in the neonates undergoing extracorporeal membrane oxygenation was significantly prolonged compared with that in the controls (17.45 ± 11.01 hour vs 5.92 ± 2.70 hour, P<0.001). Clearance decreased significantly in the extracorporeal membrane oxygenation group relative to the control group (0.03 ± 0.02 L/kg/hr vs 0.08 ± 0.05 L/kg/hr, P<0.001). No significant difference was found in the volume of distribution between the two groups (0.63 ± 0.30 L/kg in the extracorporeal membrane oxygenation group vs 0.57 ± 0.14 L/kg/hr in control, P=0.596). Clearance values were significantly correlated with serum creatinine (r = - 0.528, P<0.001). In the subgroup analysis using patients with serum creatinine < 0.5 mg/dL, similar results were obtained including significantly prolonged half-life (11.52 ± 6.31 hour vs 5.44 ± 2.36 hour, P<0.001) and decreased clearance (0.05 ± 0.02 L/kg/hr vs 0.09 ± 0.05 L/kg/hr, P<0.001) in the extracorporeal membrane oxygenation group relative to the control group.

Conclusions

Vancomycin clearance decreased significantly in the neonates undergoing extracorporeal membrane oxygenation compared with the controls. Dosing adjustments of vancomycin and close therapeutic drug monitoring are required for the safe and effective management of neonates during extracorporeal membrane oxygenation.

Physiologically based pharmacokinetic (PBPK) modeling and simulation in neonatal drug development: how clinicians can contribute

Introduction: Legal initiatives to stimulate neonatal drug development should be accompanied by development of valid research tools. Physiologically based (PB)-pharmacokinetic (PK) modeling and simulation are established tools, accepted by regulatory authorities. Consequently, PBPK holds promise to be a strong research tool to support neonatal drug development. Area covered: The currently available PBPK models still have poor predictive performance in neonates. Using an illustrative approach on distinct PK processes of absorption, distribution, metabolism, excretion, and real-world data in neonates, we provide evidence on the need to further refine available PBPK system parameters through generation and integration of new knowledge. This necessitates cross talk between clinicians and modelers to integrate knowledge (PK datasets, system knowledge, maturational physiology) or test and refine PBPK models. Expert opinion: Besides refining these models for 'small molecules', PBPK model development should also be more widely applied for therapeutic proteins and to determine exposure through breastfeeding. Researchers should also be aware that PBPK modeling in combination with clinical observations can also be used to elucidate age-related changes that are almost impossible to study based on in vivo or in vitro data. This approach has been explored for hepatic biliary excretion, renal tubular activity, and central nervous system exposure.