Renal Clearance in Newborns and Infants: Predictive Performance of Population-Based Modeling for Drug Development

Optimization of Vancomycin Initial Dosing Regimen in Neonates Using an Externally Evaluated Population Pharmacokinetic Model

BACKGROUND

Vancomycin therapeutic monitoring guidelines were revised in March 2020, and a population pharmacokinetics-guided Bayesian approach to estimate the 24-hour area under the concentration-time curve to the minimum inhibitory concentration ratio has since been recommended instead of trough concentrations. To comply with these latest guidelines, we evaluated published population pharmacokinetic models of vancomycin using an external dataset of neonatal patients and selected the most predictive model to develop a new initial dosing regimen.

METHODS

The models were identified from the literature and tested using a retrospective dataset of Canadian neonates. Their predictive performance was assessed using prediction- and simulation-based diagnostics. Monte Carlo simulations were performed to develop the initial dosing regimen with the highest probability of therapeutic target attainment.

RESULTS

A total of 144 vancomycin concentrations were derived from 63 neonates in the external population. Five of the 28 models retained for evaluation were found predictive with a bias of 15% and an imprecision of 30%. Overall, the Grimsley and Thomson model performed best, with a bias of -0.8% and an imprecision of 20.9%; therefore, it was applied in the simulations. A novel initial dosing regimen of 15 mg/kg, followed by 11 mg/kg every 8 hours should favor therapeutic target attainment.

CONCLUSIONS

A predictive population pharmacokinetic model of vancomycin was identified after an external evaluation and used to recommend a novel initial dosing regimen. The implementation of these model-based tools may guide physicians in selecting the most appropriate initial vancomycin dose, leading to improved clinical outcomes.

Pediatric Pharmacology for the Primary Care Provider: Advances and Limitations

Despite >1 in 5 children taking prescription drugs in the United States, off-label drug use is common. To increase the study of drugs in children, regulatory bodies have enacted legislation to incentivize and require pediatric drug studies. As a result of this legislation, novel trial approaches, and an increase in personnel with pediatric expertise, there have been numerous advancements in pediatric drug development. With this review, we aim to highlight developments in pediatric pharmacology over the past 6 years for the most common disease processes that may be treated pharmacologically by the pediatric primary care provider. Using information extracted from label changes between 2018 and 2023, the published literature, and Clinicaltrials.gov, we discuss advances across multiple therapeutic areas relevant to the pediatric primary care provider, including asthma, obesity and related disorders, mental health disorders, infections, and dermatologic conditions. We highlight instances in which new drugs have been developed on the basis of a deeper mechanistic understanding of illness and instances in which labels have been expanded in older drugs on the basis of newly available data. We then consider additional factors that affect pediatric drug use, including cost and nonpharmacologic therapies. Although there is work to be done, efforts focused on pediatric-specific drug development will increase the availability of evidence-based, labeled guidance for commonly prescribed drugs and improve outcomes through the safe and effective use of drugs in children.

Bayesian Vancomycin Model Selection for Therapeutic Drug Monitoring in Neonates

BACKGROUND AND OBJECTIVE

Pharmacokinetic models can inform drug dosing of vancomycin in neonates to optimize therapy. However, the model selected needs to describe the intended population to provide appropriate dose recommendations. Our study aims to identify the population pharmacokinetic (PopPK) model(s) with the best performance to predict vancomycin exposure in neonates in our hospital.

METHODS

Relevant published PopPK models for vancomycin in neonates were selected based on demographics and vancomycin dosing strategy. The predictive performance of the models was evaluated in Tucuxi using a local cohort of 69 neonates. Mean absolute error (MAE), relative bias (rBias) and relative root mean square error (rRMSE) were used to quantify the accuracy and precision of the predictive performance of each model for three different approaches: a priori, a posteriori, and Bayesian forecasting for the next course of therapy based on the previous course predictions. A PopPK model was considered clinically acceptable if rBias was between ± 20 and 95% confidence intervals included zero.

RESULTS

A total of 25 PopPK models were identified and nine were considered suitable for further evaluation. The model of De Cock et al. 2014 was the only clinically acceptable model based on a priori [MAE 0.35 mg/L, rBias 0.8 % (95% confidence interval (CI) - 7.5, 9.1%), and rRMSE 8.9%], a posteriori [MAE 0.037 mg/L, rBias - 0.23% (95% CI - 1.3, 0.88%), and rRMSE 6.02%] and Bayesian forecasting for the next courses [MAE 0.89 mg/L, rBias 5.45% (95% CI - 8.2, 19.1%), and rRMSE 38.3%) approaches.

CONCLUSIONS

The De Cock model was selected based on a comprehensive approach of model selection to individualize vancomycin dosing in our neonates.

Developmental Pharmacokinetics of Antibiotics Used in Neonatal ICU: Focus on Preterm Infants

Neonatal Infections are among the most common reasons for admission to the intensive care unit. Neonatal sepsis (NS) significantly contributes to mortality rates. Empiric antibiotic therapy of NS recommended by current international guidelines includes benzylpenicillin, ampicillin/amoxicillin, and aminoglycosides (gentamicin). The rise of antibacterial resistance precipitates the growth of the use of antibiotics of the Watch (second, third, and fourth generations of cephalosporines, carbapenems, macrolides, glycopeptides, rifamycins, fluoroquinolones) and Reserve groups (fifth generation of cephalosporines, oxazolidinones, lipoglycopeptides, fosfomycin), which are associated with a less clinical experience and higher risks of toxic reactions. A proper dosing regimen is essential for effective and safe antibiotic therapy, but its choice in neonates is complicated with high variability in the maturation of organ systems affecting drug absorption, distribution, metabolism, and excretion. Changes in antibiotic pharmacokinetic parameters result in altered efficacy and safety. Population pharmacokinetics can help to prognosis outcomes of antibiotic therapy, but it should be considered that the neonatal population is heterogeneous, and this heterogeneity is mainly determined by gestational and postnatal age. Preterm neonates are common in clinical practice, and due to the different physiology compared to the full terms, constitute a specific neonatal subpopulation. The objective of this review is to summarize the evidence about the developmental changes (specific for preterm and full-term infants, separately) of pharmacokinetic parameters of antibiotics used in neonatal intensive care units.

Physiologically-based pharmacokinetic model to investigate the effect of pregnancy on risperidone and paliperidone pharmacokinetics: Application to a pregnant woman and her neonate

This study aimed to determine the effects of pregnancy and ontogeny on risperidone and paliperidone pharmacokinetics by assessing their serum concentrations in two subjects and constructing a customized physiologically-based pharmacokinetic (PBPK) model. Risperidone and paliperidone serum concentrations were determined in a pregnant woman and her newborn. PBPK models for risperidone and paliperidone in adults, pediatric, and pregnant populations were developed and verified using the Simcyp simulator. These models were then applied to our two subjects, generating their "virtual twins." Effects of pregnancy on both drugs were examined using models with fixed pharmacokinetic parameters. In the neonatal PBPK simulation, 10 different models for estimating the renal function of neonates were evaluated. Risperidone was not detected in the serum of both pregnant woman and her newborn. Maternal and neonatal serum paliperidone concentrations were between 2.05-3.80 and 0.82-1.03 ng/ml, respectively. Developed PBPK models accurately predicted paliperidone's pharmacokinetics, as shown by minimal bias and acceptable precision across populations. The individualized maternal model predicted all observed paliperidone concentrations within the 90% prediction interval. Fixed-parameter simulations showed that CYP2D6 activity largely affects risperidone and paliperidone pharmacokinetics during pregnancy. The Flanders metadata equation showed the lowest absolute bias (mean error: 22.3% ± 6.0%) and the greatest precision (root mean square error: 23.8%) in predicting paliperidone plasma concentration in the neonatal population. Our constructed PBPK model can predict risperidone and paliperidone pharmacokinetics in pregnant and neonatal populations, which could help with precision dosing using the PBPK model-informed approach in special populations.

Population Pharmacokinetic Modeling for Ceftazidime-Avibactam Renal Dose Adjustments in Pediatric Patients 3 months and Older

Ceftazidime-avibactam is a novel β-lactam/β-lactamase inhibitor combination developed to treat serious Gram-negative bacterial infections; approved indications include complicated urinary tract infection, complicated intra-abdominal infection, and hospital-acquired pneumonia including ventilator-associated pneumonia in patients ≥ 3 months old. Because of the predominantly renal clearance of ceftazidime and avibactam, dose adjustments (reductions) are required for patients with estimated creatinine clearance (CrCL) ≤ 50 mL/min. We describe the application of combined adult and pediatric population pharmacokinetic models in developing ceftazidime-avibactam dose recommendations for pediatric patients ≥ 2 to < 18 years old with body surface area-normalized CrCL ≤ 50 mL/min/1.73 m² , including moderate, severe, or very severe renal impairment, or end-stage renal disease requiring hemodialysis, and for patients ≥ 3 months to < 2 years old with mild, moderate, or severe renal impairment. Models included allometric scaling for all subjects and simulations (1,000 subjects per age group, renal function group, and indication) were performed nonparametrically using post hoc random effects. Doses were selected based on simulated pediatric patients achieving steady-state exposures similar to adults and high probability of target attainment (using a simultaneous joint target for both ceftazidime and avibactam). Because there were few children with renal impairment in the ceftazidime-avibactam clinical trials, selected pediatric doses were guided by extrapolation and matching of adult exposures associated with efficacy and within established safety margins. The recommended doses for pediatric patients with estimated CrCL ≤ 50 mL/min/1.73 m² use equivalent adjustments in dose quantity and/or administration interval (vs. the corresponding age group with normal renal function) as those for adults.

Optimal Teicoplanin Dosing Regimen in Neonates and Children Developed by Leveraging Real-World Clinical Information

BACKGROUND

Teicoplanin is a glycopeptide antibiotic used for the treatment of methicillin-resistant Staphylococcus aureus infections. To ensure successful target attainment, therapeutic drug monitoring-informed dosage adjustment is recommended. However, it relies on the experience of the clinician and the frequency of drug measurements. This study aimed to design a new optimal dosing regimen of teicoplanin with a maintenance dosing strategy for neonates and children based on their physiological characteristics.

METHODS

Data from teicoplanin-treated patients (n = 214) were collected from electronic medical records. Covariate analyses were performed using population pharmacokinetic (PK) modeling with 399 serum teicoplanin concentrations from 48 neonates and 166 children. Multiple PK simulations were conducted to explore optimal dosing regimens that would allow control of the trough concentration to the target of 15-30 mg/L quicker than the current standard regimen.

RESULTS

Allometrically scaled body weight, postmenstrual age (PMA), renal function, and serum albumin were implemented as substantial covariates for teicoplanin clearance in a two-compartment PK model. Covariate analyses and comprehensive simulation assessments recommended the following modifications to the current regimen: (1) decreased dose for premature babies (PMA ≤28 weeks), (2) decreased dose for children with renal dysfunction, and (3) increased dose for children (0.5-11 years) with an estimated glomerular filtration rate of ≥90 mL/min/1.73 m2.

CONCLUSIONS

This study leverages real-world clinical information and proposes new optimal dosing regimens for teicoplanin in neonates and children through PK modeling and simulation analyses, taking into account the age, including PMA, and renal function of patients.

Fluconazole Population Pharmacokinetics after Fosfluconazole Administration and Dosing Optimization in Extremely Low-Birth-Weight Infants

A prospective single-center study was conducted to characterize the pharmacokinetics (PK) of fluconazole (FLCZ) in extremely low-birth-weight infants (ELBWIs) who received fosfluconazole (F-FLCZ). Intravenous F-FLCZ was administered at a dose of 3 mg/kg of body weight every 72 h during the first 2 weeks of life, every 48 h during the third and fourth weeks of life, and every 24 h after 5 weeks of life. Blood samples from ELBWIs treated with F-FLCZ were collected using scavenged samples. The concentration of FLCZ was determined using liquid chromatography-tandem mass spectrometry. The population pharmacokinetic model was established using Phenix NLME 8.2 software. In total, 18 ELBWIs were included in this analysis. Individual PK parameters were determined by a one-compartment analysis with first-order conversion. Postmenstrual age (PMA), serum creatinine (SCr), and alkaline phosphatase were considered covariates for clearance (CL). The mean population CL and the volume of distribution were 0.011 L/h/kg^0.75 and 0.95 L/kg, respectively. Simulation assessments with the final model revealed that the current regimen (3 mg/kg every 72 h) could achieve the proposed target FLCZ trough concentration (>2 μg/mL) in 43.3% and 72.2% of infants with a PMA of ≥37 and 30 to 36 weeks, respectively, and an SCr level of <0.5 mg/dL. Shortened dosing intervals (every 48 or 24 h) might improve the probability of target attainment. This study was the first to assess the PK of F-FLCZ in ELBWI, as well as the first to provide fundamental information about FLCZ exposure after F-FLCZ administration, with the goal of facilitating dose optimization in the ELBWI population.

IMPORTANCE Invasive fungal infection is an important cause of mortality and morbidity in very preterm or very-low-birth-weight infants. In order to limit the risk of invasive fungal infections in this population, the administration of fluconazole is generally recommended for extremely low-birth-weight infants admitted to a neonatal intensive care unit with a Candida species colonization prevalence rate of >10%, under the guidelines of the Infectious Diseases Society of America. Fosfluconazole can reduce the volume of solution required for intravenous therapy compared to fluconazole because it has increased solubility, which is a major advantage for infants undergoing strict fluid management. To date, no study has demonstrated the fluconazole pharmacokinetics after fosfluconazole administration in neonates and infants, and this needs to be clarified. Here, we characterized the pharmacokinetics of fluconazole in extremely low-birth-weight infants who received F-FLCZ and explored the appropriate dosage in this patient population.

Simultaneous pharmacokinetic/pharmacodynamic (PKPD) assessment of ampicillin and gentamicin in the treatment of neonatal sepsis

Objectives

This study aimed to simultaneously investigate the pharmacokinetics of ampicillin and gentamicin, currently the WHO standard of care for treating neonatal sepsis.

Methods

Pharmacokinetic data were collected in 59 neonates receiving ampicillin and gentamicin for suspected or proven sepsis in the NeoFosfo trial (NCT03453177). A panel of 23 clinical Escherichia coli isolates from neonates with sepsis, resistant to either ampicillin, gentamicin or both, were tested for susceptibility using chequerboards. Pharmacokinetic/pharmacodynamic (PKPD) modelling and simulations were used to compare single-agent (EUCAST MIC) and combination (chequerboard MIC) target attainment with standard dosing regimens.

Results

A model was established that simultaneously estimated parameters of a one-compartment ampicillin model and a two-compartment gentamicin model. A common clearance for both drugs was used (6.89 L/h/70 kg) relating to glomerular filtration (CL_GFR), with an additional clearance term added for ampicillin (5.3 L/h/70 kg). Covariate modelling included a priori allometric weight and post-menstrual age scaling of clearance. Further covariate relationships on renal clearance were postnatal age and serum creatinine.

Simulation-based PKPD assessments suggest good Gram-positive (MIC ≤ 0.25 mg/L) cover. However, less than one-quarter of neonates were predicted to receive efficacious coverage against Enterobacterales (MIC ≤ 2 mg/L). The benefit of the ampicillin/gentamicin combination was limited, with only 2/23 E. coli clinical strains showing FIC index < 0.5 (synergy) and most in the range 0.5–1 (suggesting additivity). Simulations showed that feasible dosing strategies would be insufficient to cover resistant strains.

Conclusions

PKPD simulations showed ampicillin and gentamicin combination therapy was insufficient to cover Enterobacterales, suggesting the need for alternative empirical treatment options for neonatal sepsis.

An update on the use of allometric and other scaling methods to scale drug clearance in children: towards decision tables

INTRODUCTION

When pediatric data are not available for a drug, allometric and other methods are applied to scale drug clearance across the pediatric age-range from adult values. This is applied when designing first-in-child studies, but also for off-label drug prescription.

AREAS COVERED

This review provides an overview of the systematic accuracy of allometric and other pediatric clearance scaling methods compared to gold-standard PBPK predictions. The findings are summarized in decision tables to provide a priori guidance on the selection of appropriate pediatric clearance scaling methods for both novel drugs for which no pediatric data are available and existing drugs in clinical practice.

EXPERT OPINION

While allometric scaling principles are commonly used to scale pediatric clearance, there is no universal allometric exponent (i.e., 1, 0.75 or 0.67) that can accurately scale clearance for all drugs from adults to children of all ages. Therefore, pediatric scaling decision tables based on age, drug elimination route, binding plasma protein, fraction unbound, extraction ratio, and/or isoenzyme maturation are proposed to a priori select the appropriate (allometric) clearance scaling method, thereby reducing the need for full PBPK-based clearance predictions. Guidance on allometric scaling when estimating pediatric clearance values is provided as well.

Prediction of Pregnancy-Induced Changes in Secretory and Total Renal Clearance of Drugs Transported by Organic Anion Transporters

Pregnancy can significantly change the pharmacokinetics of drugs, including those renally secreted by organic anion transporters (OATs). Quantifying these changes in pregnant women is logistically and ethically challenging. Hence, predicting the in vivo plasma renal secretory clearance (CLsec) and renal CL (CLrenal) of OAT drugs in pregnancy is important to design correct dosing regimens of OAT drugs. Here, we first quantified the fold-change in renal OAT activity in pregnant versus nonpregnant individual using available selective OAT probe drug CLrenal data (training dataset; OAT1: tenofovir, OAT2: acyclovir, OAT3: oseltamivir carboxylate). The fold-change in OAT1 activity during the 2nd and 3rd trimester was 2.9 and 1.0 compared with nonpregnant individual, respectively. OAT2 activity increased 3.1-fold during the 3rd trimester. OAT3 activity increased 2.2, 1.7 and 1.3-fold during the 1st, 2nd, and 3rd trimester, respectively. Based on these data, we predicted the CLsec, CLrenal and total clearance ((CLtotal) of drugs in pregnancy, which are secreted by multiple OATs (verification dataset; amoxicillin, pravastatin, cefazolin and ketorolac, R-ketorolac, S-ketorolac). Then, the predicted clearances (CLs) were compared with the observed values. The predicted/observed CLsec, CLrenal, and CLtotal of drugs in pregnancy of all verification drugs were within 0.80-1.25 fold except for CLsec of amoxicillin in the 3rd trimester (0.76-fold) and cefazolin in the 2nd trimester (1.27-fold). Overall, we successfully predicted the CLsec, CLrenal, and CLtotal of drugs in pregnancy that are renally secreted by multiple OATs. This approach could be used in the future to adjust dosing regimens of renally secreted OAT drugs which are administered to pregnant women. SIGNIFICANCE STATEMENT: To the authors' knowledge, this is the first report to successfully predict renal secretory clearance and renal clearance of multiple OAT substrate drugs during pregnancy. The data presented here could be used in the future to adjust dosing regimens of renally secreted OAT drugs in pregnancy. In addition, the mechanistic approach used here could be extended to drugs transported by other renal transporters.

Physiologically Based Pharmacokinetic Models Are Effective Support for Pediatric Drug Development

Pediatric drug development faces many difficulties. Traditionally, pediatric drug doses are simply calculated linearly based on the body weight, age, and body surface area of adults. Due to the ontogeny of children, this simple linear scaling may lead to drug overdose in pediatric patients. The physiologically based pharmacokinetic (PBPK) model, as a mathematical model, contributes to the research and development of pediatric drugs. An example of a PBPK model guiding drug dose selection in pediatrics has emerged and has been approved by the relevant regulatory agencies. In this review, we discuss the principle of the PBPK model, emphasize the necessity of establishing a pediatric PBPK model, introduce the absorption, distribution, metabolism, and excretion of the pediatric PBPK model, and understand the various applications and related prospects of the pediatric PBPK model.

IV and oral fosfomycin pharmacokinetics in neonates with suspected clinical sepsis

BACKGROUND

Fosfomycin has the potential to be re-purposed as part of a combination therapy to treat neonatal sepsis where resistance to current standard of care (SOC) is common. Limited data exist on neonatal fosfomycin pharmacokinetics and estimates of bioavailability and CSF/plasma ratio in this vulnerable population are lacking.

OBJECTIVES

To generate data informing the appropriate dosing of IV and oral fosfomycin in neonates using a population pharmacokinetic analysis of plasma and CSF data.

METHODS

The NeoFosfo study (NCT03453177) was a randomized trial that examined the safety and pharmacokinetics of fosfomycin comparing SOC versus SOC plus fosfomycin. Sixty-one neonates received fosfomycin (100 mg/kg IV q12h for 48 h) and then they converted to oral therapy at the same dose. Two plasma pharmacokinetic samples were taken following the first IV and oral doses, sample times were randomized to cover the whole pharmacokinetic profile and opportunistic CSF pharmacokinetic samples were collected. A population pharmacokinetic model was developed in NONMEM and simulations were performed.

RESULTS

In total, 238 plasma and 15 CSF concentrations were collected. A two-compartment disposition model, with an additional CSF compartment and first-order absorption, best described the data. Bioavailability was estimated as 0.48 (95% CI = 0.347-0.775) and the CSF/plasma ratio as 0.32 (95% CI = 0.272-0.409). Allometric weight and postmenstrual age (PMA) scaling was applied; additional covariates included postnatal age (PNA) on clearance and CSF protein on CSF/plasma ratio.

CONCLUSIONS

Through this analysis a population pharmacokinetic model has been developed that can be used alongside currently available pharmacodynamic targets to select a neonatal fosfomycin dose based on an infant's PMA, PNA and weight.

Application of Physiologically Based Pharmacokinetic-Pharmacodynamic Modeling in Preterm Neonates to Guide Gentamicin Dosing Decisions and Predict Antibacterial Effect

Clinical studies in preterm neonates are rarely performed due to ethical concerns and difficulties associated with trials and recruitment. Consequently, dose selection in this population is primarily empirical. Scaling neonatal doses from adult doses does not account for developmental changes and may not accurately predict drug kinetics. This is especially important for gentamicin, a narrow therapeutic index aminoglycoside antibiotic. While gentamicin's bactericidal effect is associated with its peak plasma concentration, keeping trough concentrations below 1 µg/mL prevents toxicity and also helps to counteract adaptive resistance in bacteria such as Escherichia coli. In this study, physiologically based pharmacokinetic-pharmacodynamic (PBPK-PD) modeling was used to support and/or guide dosing decisions and to predict the antibacterial effect in preterm neonates. A gentamicin PBPK model was successfully verified in healthy adults and preterm neonates across all gestational ages. Clinical data from a neonatal intensive care unit at NYU Langone Hospital-Long Island was used to identify dosing regimens associated with increased incidence of elevated gentamicin trough concentrations in different preterm patient cohorts. Model predictions demonstrated that a higher dose with an extended-dosing interval (every 36 hours) in neonates with a postmenstrual age of 30 to 34 weeks and ≥35 weeks, with postnatal age 8 to 28 days and 0 to 7 days, respectively, were more likely to have a trough <1 µg/mL when compared with once-daily (every 24 hours) dosing. PBPK-PD modeling suggested that a higher dose administered every 36 hours may provide effective antibacterial therapy.

Drug Clearance in Neonates: A Combination of Population Pharmacokinetic Modelling and Machine Learning Approaches to Improve Individual Prediction

BACKGROUND

Population pharmacokinetic evaluations have been widely used in neonatal pharmacokinetic studies, while machine learning has become a popular approach to solving complex problems in the current era of big data.

OBJECTIVE

The aim of this proof-of-concept study was to evaluate whether combining population pharmacokinetic and machine learning approaches could provide a more accurate prediction of the clearance of renally eliminated drugs in individual neonates.

METHODS

Six drugs that are primarily eliminated by the kidneys were selected (vancomycin, latamoxef, cefepime, azlocillin, ceftazidime, and amoxicillin) as 'proof of concept' compounds. Individual estimates of clearance obtained from population pharmacokinetic models were used as reference clearances, and diverse machine learning methods and nested cross-validation were adopted and evaluated against these reference clearances. The predictive performance of these combined methods was compared with the performance of two other predictive methods: a covariate-based maturation model and a postmenstrual age and body weight scaling model. Relative error was used to evaluate the different methods.

RESULTS

The extra tree regressor was selected as the best-fit machine learning method. Using the combined method, more than 95% of predictions for all six drugs had a relative error of < 50% and the mean relative error was reduced by an average of 44.3% and 71.3% compared with the other two predictive methods.

CONCLUSION

A combined population pharmacokinetic and machine learning approach provided improved predictions of individual clearances of renally cleared drugs in neonates. For a new patient treated in clinical practice, individual clearance can be predicted a priori using our model code combined with demographic data.

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.

Serum Creatinine and Serum Cystatin C are Both Relevant Renal Markers to Estimate Vancomycin Clearance in Critically Ill Neonates

Purpose: Serum creatinine (SCr) is used as a marker of kidney function to guide dosing of renally eliminated drugs. Serum Cystatin C (S-CysC) has been suggested as a more reliable kidney marker than SCr in adults and children. Purpose of this study was to investigate S-CysC as alternative renal marker to SCr for estimating vancomycin clearance in neonates undergoing intensive care.

Methods: Vancomycin pharmacokinetics (PK), SCr and S-CysC data were collected in patients undergoing vancomycin treatment in the neonatal intensive care unit of Robert Debré Hospital - Paris. A population PK analysis was performed utilizing routine therapeutic drug monitoring samples. S-CysC and SCr were compared as covariates on vancomycin clearance using stepwise covariate modeling (forward inclusion [p < 0.05] and backward elimination [p < 0.01]). Model performance was evaluated by graphical and statistical criteria.

Results: A total of 108 vancomycin concentrations from 66 patients (postmenstrual age [PMA] of 26–46 weeks) were modeled with an allometric one-compartment model. The median (range) values for SCr and S-CysC were 41 (12–153) µmol/l and 1.43 (0.95–2.83) mg/l, respectively. Following stepwise covariate model building, SCr was retained as single marker of kidney function (after accounting for weight and PMA) in the final model. Compared to the final model based on SCr, the alternative model based on S-CysC showed very similar performance (e.g. BIC of 578.3 vs. 576.4) but included one additional covariate: impact of mechanical ventilation on vancomycin clearance, in addition to the effects of size and maturation.

Conclusion: ill neonates. However, if using S-CysC for this purpose mechanical ventilation needs to be taken into account.

Glomerular Filtration Rate in Asphyxiated Neonates Under Therapeutic Whole-Body Hypothermia, Quantified by Mannitol Clearance

Background

Therapeutic hypothermia (TH) is an established intervention to improve the outcome of neonates with moderate-to-severe hypoxic-ischemic encephalopathy resulting from perinatal asphyxia. Despite this beneficial effect, TH may further affect drug elimination pathways such as the glomerular filtration rate.

Objectives

The objective of this study was to quantify the effect of TH in addition to asphyxia on mannitol clearance as a surrogate for the glomerular filtration rate.

Methods

The effect of asphyxia and TH (mild vs moderate/severe) on mannitol clearance was assessed using a population approach, based on mannitol observations collected in the ALBINO (ALlopurinol in addition to TH for hypoxic-ischemic Brain Injury on Neurocognitive Outcome) trial, as some were exposed to a second dose of 10 mg/kg intravenous mannitol as placebo to ensure blinding. Pharmacokinetic analysis and model development were conducted using NONMEM version 7.4.

Results

Based on 77 observations from 17 neonates (TH = 13), a one-compartment model with first-order linear elimination best described the observed data. To account for prenatal glomerular filtration rate maturation, both birthweight and gestational age were implemented as clearance covariates using an earlier published three-quarters power function and a sigmoid hyperbolic function. Our final model predicted a mannitol clearance of 0.15 L/h for a typical asphyxia neonate (39.5 weeks, birthweight 3.25 kg, no TH), lower than the reported value of 0.33 L/h for a healthy neonate of similar age and weight. By introducing TH as a binary covariate on clearance, the additional impact of TH on mannitol clearance was quantified (60% decrease).

Conclusions

Mannitol clearance was decreased by approximately 60% in neonates undergoing TH, although this is likely confounded with asphyxia severity.

Trial Registration

ClinicalTrials.gov identifier NCT03162653.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40262-021-00991-6.

Dosing Recommendations for Pediatric Patients With Renal Impairment

A treatment gap exists for pediatric patients with renal impairment. Alterations in renal clearance and metabolism of drugs render standard dosage regimens inappropriate and may lead to drug toxicity, but these studies are not routinely conducted during drug development. The objective of this study was to examine the clinical evidence behind current renal impairment dosage recommendations for pediatric patients in a standard pediatric dosing handbook. The sources of recommendations and comparisons included the pediatric dosing handbook (Lexicomp), the U.S. Food and Drug Administration-approved manufacturer's labels, and published studies in the literature. One hundred twenty-six drugs in Lexicomp had pediatric renal dosing recommendations. Only 14% (18 of 126) of Lexicomp pediatric renal dosing recommendations referenced a pediatric clinical study, and 15% of manufacturer's labels (19 of 126) described specific dosing regimens for renally impaired pediatric patients. Forty-two products had published information on pediatric renal dosing, but 19 (45%) were case studies. When pediatric clinical studies were not referenced in Lexicomp, the renal dosing recommendations followed the adult and pediatric dosing recommendations on the manufacturer's label. Clinical evidence in pediatric patients does not exist for most renal dosing recommendations in a widely used pediatric dosing handbook, and the adult renal dosing recommendations from the manufacturer's label are currently the primary source of pediatric renal dosing information.

Application of physiologically based pharmacokinetic modeling for sertraline dosing recommendations in pregnancy

Pregnancy is a period of significant change that impacts physiological and metabolic status leading to alterations in the disposition of drugs. Uncertainty in drug dosing in pregnancy can lead to suboptimal therapy, which can contribute to disease exacerbation. A few studies show there are increased dosing requirements for antidepressants in late pregnancy; however, the quantitative data to guide dose adjustments are sparse. We aimed to develop a physiologically based pharmacokinetic (PBPK) model that allows gestational-age dependent prediction of sertraline dosing in pregnancy. A minimal physiological model with defined gut, liver, plasma, and lumped placental-fetal compartments was constructed using the ordinary differential equation solver package, ‘mrgsolve’, in R. We extracted data from the literature to parameterize the model, including sertraline physicochemical properties, in vitro metabolism studies, disposition in nonpregnant women, and physiological changes during pregnancy. The model predicted the pharmacokinetic parameters from a clinical study with eight subjects for the second trimester and six subjects for the third trimester. Based on the model, gestational-dependent changes in physiology and metabolism account for increased clearance of sertraline (up to 143% at 40 weeks gestational age), potentially leading to under-dosing of pregnant women when nonpregnancy doses are used. The PBPK model was converted to a prototype web-based interactive dosing tool to demonstrate how the output of a PBPK model may translate into optimal sertraline dosing in pregnancy. Quantitative prediction of drug exposure using PBPK modeling in pregnancy will support clinically appropriate dosing and increase the therapeutic benefit for pregnant women.

The Influence of Haemostatic System Maturation on the Dose-Response Relationship of Unfractionated Heparin

BACKGROUND

Unfractionated heparin (UFH) dosing and monitoring guidelines for children are often extrapolated from adult data. This practice is suboptimal given the inherent differences in haemostatic maturation and drug handling in children compared with adults.

OBJECTIVE

The aim of this work was to investigate the impact of haemostatic system maturation on the dose-response relationship of UFH in children.

METHODS

A quantitative model for haemostasis in adults was adapted to account for maturation in UFH pharmacokinetic (PK) parameters with and without age-related changes in coagulation factor concentrations. The adult and adapted models were used to predict the time courses of anti-factor Xa activity (aXa) and activated partial thromboplastin time (aPTT) in patients receiving UFH infusion. Predictions from both models were compared with observed aXa and aPTT measurements from 31 paediatric patients receiving UFH during extracorporeal membrane oxygenation (ECMO).

RESULTS

The model with maturation for both UFH PK and the haemostatic system had an improved aXa and aPTT predictive performance compared with maturation in UFH PK only and the original adult model. Despite the minor effect of haemostatic system maturation on baseline aPTT, it led to substantial changes in the time course of aPTT sensitivity to UFH. This finding suggests that between-subject variability in clotting factors concentrations is potentially a major contributor to the overall variability of aPTT response to UFH. In addition, time-varying clotting factors concentrations may explain within-subject changes in aPTT sensitivity to UFH.

CONCLUSION

We developed the first quantitative systems pharmacology (QSP) model that provides a mechanistic and quantitative basis for linking physiological and pharmacological maturation to UFH effect and response biomarkers. After appropriate clinical validation, the model could be useful for the development of paediatric-specific individualised UFH dosing recommendations.

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.

β-Lactam antimicrobial pharmacokinetics and target attainment in critically ill patients aged 1 day to 90 years: the ABDose study

Background

The pharmacokinetics of β-lactam antibiotics in critical illness remain poorly characterized, particularly in neonates, children and the elderly. We undertook a pharmacokinetic study of commonly used β-lactam antibiotics in critically ill patients of all ages. The aims were to produce a whole-life β-lactam pharmacokinetic model and describe the extent to which standard doses achieve pharmacokinetic/pharmacodynamic targets associated with clinical cure.

Patients and methods

A total of 212 critically ill participants with an age range from 1 day (gestational age 24 weeks) to 90 years were recruited from a UK hospital, providing 1339 pharmacokinetic samples. Population pharmacokinetic analysis was undertaken using non-linear mixed-effects modelling (NONMEM) for each drug. Pooled data were used to estimate maturation and decline of β-lactam pharmacokinetics throughout life.

Results

Pharmacokinetic models for eight drugs were described, including what is thought to be the first benzylpenicillin model in critically ill adults. We estimate that 50% of adult β-lactam clearance is achieved by 43 weeks post-menstrual age (chronological plus gestational age). Fifty percent of decline from peak adult clearance occurs by 71 years. Paediatric participants were significantly less likely than adults to achieve pharmacokinetic/pharmacodynamic targets with standard antibiotic doses (P &lt; 0.01).

Conclusions

We believe this to be the first prospective whole-life antibiotic pharmacokinetic study in the critically ill. The study provides further evidence that standard antibiotic doses fail to achieve pharmacokinetic/pharmacodynamic targets associated with clinical success in adults, children and neonates. Maturation and decline parameters estimated from this study could be adopted as a standard for future prospective studies.

Predictive Pediatric Modeling and Simulation Using Ontogeny Information

Food and Drug Administration submissions of physiologically based pharmacokinetic (PBPK) modeling and simulation of small-molecule drugs document the relevance of pediatric drug development and, in particular, information on dosing strategies in children. The most relevant prerequisite for reliable PBPK-based translation of adult pharmacokinetics of a small molecule to children is knowledge of the drug-specific absorption, distribution, metabolism, and elimination (ADME) processes in adults together with existing information about ontogeny of ADME processes relevant for the drug. All mechanisms driving a drug's clearance are of specific importance. For other drug modalities, our knowledge of ADME processes and ontogeny is still limited. More research is required, for example, to understand why some therapeutic proteins show complex differences in pharmacokinetics between adults and children, whereas other proteins seem to follow simple allometric scaling rules. Ontogeny information originates from various sources, such as (semi)quantitative mRNA expression, in vitro activity data, and deconvolution of in vivo pharmacokinetic data. The workflow for pediatric predictions is well described in several articles documenting successful translation from adults to children. The technical hurdles for PBPK modeling are low. State-of-the-art PBPK modeling software tools provide integrated pediatric translation workflows. For example, PK-Sim and MoBi are freely available as fully transparent open-source software via Open Systems Pharmacology (OSP). With the latest 2019 software release, version 8.0, OSP even provides a fully integrated technical framework for the qualification (and requalification) of any specific intended PBPK use in line with Food and Drug Administration and European Medicines Agency PBPK guidance. Qualification packages for pediatric translation are available on the OSP platform.

Pediatric Renal Ontogeny and Applications in Drug Development

The clinical applications of renal ontogeny mainly include renal function evaluation and optimal dosing of renally eliminated drugs in pediatric patients, which rely on pharmacometric models and/or bedside estimated glomerular filtration rate equations. However, these applications in drug development are based on an understanding of renal function development, especially when considering premature infants, and the renal biomarkers that can be used for renal function assessment. This review provides a general overview on (1) renal function development, (2) the biomarkers that are used to assess renal function, and (3) the practical application of this knowledge to drug dosing for renally eliminated drugs during pediatric development. While pharmacometric approaches for estimating renal function during development have improved considerably, the number of drug development programs that have studied premature infants is small and suggests that caution should be taken in estimating doses for renally eliminated drugs during periods of rapid change in renal function.

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.

Creatinine-Based Renal Function Assessment in Pediatric Drug Development: An Analysis Using Clinical Data for Renally Eliminated Drugs

The estimated glomerular filtration rate (eGFR) equations based on serum creatinine (SCR) have been used for pediatric dose adjustment in drug labeling. This study evaluated the performance of those equations in estimating individual clearance of drugs that are predominantly eliminated by glomerular filtration, using clinical data from the renally eliminated drugs gadobutrol, gadoterate, amikacin, and vancomycin. The eGFR was compared with the observed drug clearance (CL) in 352 pediatric patients from birth to 12 years of age. Multiple eGFR equations overestimated the drug CL on average, including the original and bedside Schwartz equations, which showed an average eGFR/CL ratio between 1 and 3. Further analysis with bedside Schwartz equation showed a higher eGFR/CL ratio in the subjects with a lower SCR or CL. Supraphysiological eGFR as high as 380 mL/min/1.73 m² was obtained using the bedside Schwartz equation for some of the subjects, most of whom are children < 2 years of age with SCR < 0.2 mg/dL. Excluding the subjects with supraphysiological eGFR from the analysis did not change the overall trend of overestimation. In conclusion, Schwartz equations led to an overestimation of drug clearance for the drugs evaluated. When greater precision is required in predicting eGFR for pediatric patients, such as in drug dosing, revised k constants for the Schwartz equation or new methods of glomerular filtration rate estimation may be necessary.