Population pharmacokinetic studies in pediatrics: issues in design and analysis

Population pharmacokinetic models of anti-PD-1 mAbs in patients with multiple tumor types: A systematic review

Aims and background

A number of population pharmacokinetic (PPK) models of anti-programmed cell death-1 (PD-1) monoclonal antibodies (mAbs) in multiple tumor types have been published to characterize the influencing factors of their pharmacokinetics. This review described PPK models of anti-PD-1 mAbs that investigate the magnitude and types of covariate effects in PK parameters, provide a reference for building PPK models of other anti-PD-1 mAbs, and identify areas requiring additional research to facilitate the application of PPK models.

Methods

A systematic search for analyses of PPK models of eleven anti-PD-1 mAbs on the market that were carried out in humans was conducted using PubMed, Embase, and the Cochrane Library. The search covered the period from the inception of the databases to April 2022.

Results

Currently, there are fourteen analyses on PPK models of anti-PD-1 mAbs summarized in this review, including seven models that refer to nivolumab, four referring to pembrolizumab, one referring to cemiplimab, one referring to camrelizumab, and one referred to dostarlimab. Most analyses described the pharmacokinetics of anti-PD-1 mAbs with a two-compartment model with time-varying clearance (CL) and a sigmoidal maximum effect. The estimated CL and volume of distribution in the central (V_C) ranged from 0.179 to 0.290 L/day and 2.98 to 4.46 L, respectively. The median (range) of interindividual variability (IIV) for CL and V_C was 30.9% (8.7%–50.8%) and 29.0% (4.32%–40.7%), respectively. The commonly identified significant covariates were body weight (BW) on CL and V_C, and albumin (ALB), tumor type, sex, and performance status (PS) on CL. Other less assessed significant covariates included lactate dehydrogenase (LDH), immunoglobulin G (IgG), ipilimumab coadministration (IPICO) on CL, and body mass index (BMI), malignant pleural mesothelioma (MESO) on V_C.

Conclusion

This review provides detailed information about the characteristics of PPK models of anti-PD-1 mAbs, the effects of covariates on PK parameters, and the current status of the application of the models. ALB, BW, specific tumor type, sex, and PS should be considered for the future development of the PPK model of anti-PD-1 mAbs. Other potential covariates that were assessed less frequently but still have significance (e.g., LDH, IgG, and IPICO) should not be ignored. Thus, further research and thorough investigation are needed to assess new or potential covariates, which will pave the way for personalized anti-PD-1 mAbs therapy.

A General Biphasic Bodyweight Model for Scaling Basal Metabolic Rate, Glomerular Filtration Rate, and Drug Clearance from Birth to Adulthood

The objective of this study is to propose a unified, continuous, and bodyweight-only equation to quantify the changes of human basal metabolic rate (BMR), glomerular filtration rate (GFR), and drug clearance (CL) from infancy to adulthood. The BMR datasets were retrieved from a comprehensive historical database of male and female subjects (0.02 to 64 years). The CL datasets for 17 drugs and the GFR dataset were generated from published maturation and growth models with reported parameter values. A statistical approach was used to simulate the model-generated CL and GFR data for a hypothetical population with 26 age groups (from 0 to 20 years). A biphasic equation with two power-law functions of bodyweight was proposed and evaluated as a general model using nonlinear regression and dimensionless analysis. All datasets universally reveal biphasic curves with two distinct linear segments on log-log plots. The biphasic equation consists of two reciprocal allometric terms that asymptotically determine the overall curvature. The fitting results show a superlinear scaling phase (asymptotic exponent >1; ca. 1.5-3.5) and a sublinear scaling phase (asymptotic exponent <1; ca. 0.5-0.7), which are separated at the phase transition bodyweight ranging from 5 to 20 kg with a mean value of 10 kg (corresponding to 1 year of age). The dimensionless analysis generalizes and offers quantitative realization of the maturation and growth process. In conclusion, the proposed mixed-allometry equation is a generic model that quantitatively describes the phase transition in the human maturation process of diverse human functions.

Knowledge Gaps in the Pharmacokinetics of Therapeutic Proteins in Pediatric Patients

Therapeutic proteins such as monoclonal antibodies and their derivatives, fusions proteins, hormone analogs and enzymes for replacement therapy are an ever-growing mainstay in our pharmacopoeia. While a growing number of these medications are developed for and used in younger and younger pediatric patients, knowledge gaps in the basic understanding of the molecular and physiologic processes governing the disposition of these compounds in the human body and their modulation by age and childhood development are a hindrance to the effective and timely development and clinical use of these compounds, especially in very young pediatric patient populations. This is particularly the case for the widespread lack of information on the ontogeny and age-associated expression and function of receptor systems that are involved in the molecular processes driving the pharmacokinetics of these compounds. This article briefly highlights three receptor systems as examples, the neonatal Fc receptor, the asialoglycoprotein receptor, and the mannose receptor. It furthermore provides suggestions on how these gaps should be addressed and prioritized to provide the field of pediatric clinical pharmacology the urgently needed tools for a more effective development and clinical utilization of this important class of drugs with rapidly evolving importance as cornerstone in pediatric pharmacotherapy.

The incidence rate and influence factors of hemolysis, lipemia, icterus in fasting serum biochemistry specimens

Objective

Hemolysis, icterus, and lipemia (HIL) of blood samples have been a concern in hospitals because they reflect pre-analytical processes’ quality control. However, very few studies investigate the influence of patients’ gender, age, and department, as well as sample-related turnaround time, on the incidence rate of HIL in fasting serum biochemistry specimens.

Methods

A retrospective, descriptive study was conducted to investigate the incidence rate of HIL based on the HIL index in 501,612 fasting serum biochemistry specimens from January 2017 to May 2018 in a tertiary university hospital with 4,200 beds in Sichuan, southwest China. A subgroup analysis was conducted to evaluate the differences in the HIL incidence rate by gender, age and department of patients, and turnaround time of specimens.

Results

The incidence rate of hemolysis, lipemia and icterus was 384, 53, and 612 per 10,000 specimens. The male patients had a significantly elevated incidence of hemolysis (4.13% vs. 3.54%), lipemia (0.67% vs. 0.38%), and icterus (6.95% vs. 5.43%) than female patients. Hemolysis, lipemia, and icterus incidence rate were significantly associated with the male sex with an odds ratio (OR) of 1.174 [95% confidence interval (CI), 1.140–1.208], 1.757 (95%CI: 1.623–1.903), and 1.303 (95%CI: 1.273–1.333), respectively, (P<0.05). The hospitalized patients had a higher incidence of hemolysis (4.03% vs. 3.54%), lipemia (0.63% vs. 0.36%), and icterus (7.10% vs. 4.75%) than outpatients (P<0.001). Specimens with relatively longer transfer time and/or detection time had a higher HIL incidence (P<0.001). The Pediatrics had the highest incidence of hemolysis (16.2%) with an adjusted OR (AOR) of 4.93 (95%CI, 4.59–5.29, P<0.001). The Neonatology department had the highest icterus incidence (30.1%) with an AOR of 4.93 (95%CI: 4.59–5.29, P<0.001). The Neonatology department (2.32%) and Gastrointestinal Surgery (2.05%) had the highest lipemia incidence, with an AOR of 1.17 (95%CI: 0.91–1.51) and 4.76 (95%CI: 4.70–5.53), both P-value <0.001. There was an increasing tendency of hemolysis and icterus incidence for children under one year or adults aged more than 40.

Conclusion

Evaluation of HIL incidence rate and HIL-related influence factors in fasting serum biochemistry specimens are impartment to interpret the results more accurately and provide better clinical services to patients.

Bioinformatics and Computer Simulation approaches to the discovery and analysis of Bioactive Peptides

The traditional process of separating and purifying bioactive peptides is laborious and time-consuming. Using a traditional process to identify is difficult, and there is a lack of fast and accurate activity evaluation methods. How to extract bioactive peptides quickly and efficiently is still the focus of bioactive peptides research. In order to improve the present situation of the research, bioinformatics techniques and peptidome methods are widely used in this field. At the same time, bioactive peptides have their own specific pharmacokinetic characteristics, so computer simulation methods have incomparable advantages in studying the pharmacokinetics and pharmacokinetic-pharmacodynamic correlation models of bioactive peptides. The purpose of this review is to summarize the combined applications of bioinformatics and computer simulation methods in the study of bioactive peptides, with focuses on the role of bioinformatics in simulating the selection of enzymatic hydrolysis and precursor proteins, activity prediction, molecular docking, physicochemical properties, and molecular dynamics. Our review shows that new bioactive peptide molecular sequences with high activity can be obtained by computer-aided design. The significance of the pharmacokinetic-pharmacodynamic correlation model in the study of bioactive peptides is emphasized. Finally, some problems and future development potential of bioactive peptides binding new technologies are prospected.

Model-Based Exposure-Response Assessment for Spectinamide 1810 in a Mouse Model of Tuberculosis

Despite decades of research, tuberculosis remains a leading cause of death from a single infectious agent. Spectinamides are a promising novel class of antituberculosis agents, and the lead spectinamide 1810 has demonstrated excellent efficacy, safety, and drug-like properties in numerous in vitro and in vivo assessments in mouse models of tuberculosis. In the current dose ranging and dose fractionation study, we used 29 different combinations of dose level and dosing frequency to characterize the exposure-response relationship for spectinamide 1810 in a mouse model of Mycobacterium tuberculosis infection and in healthy animals. The obtained data on 1810 plasma concentrations and counts of CFU in lungs were analyzed using a population pharmacokinetic/pharmacodynamic (PK/PD) approach as well as classical anti-infective PK/PD indices. The analysis results indicate that there was no difference in the PK of 1810 in infected compared to healthy, uninfected animals. The PK/PD index analysis showed that bacterial killing of 1810 in mice was best predicted by the ratio of maximum free drug concentration to MIC (fC_max/MIC) and the ratio of the area under the free concentration-time curve to the MIC (fAUC/MIC) rather than the cumulative percentage of time that the free drug concentration is above the MIC (f%T_MIC). A novel PK/PD model with consideration of postantibiotic effect could adequately describe the exposure-response relationship for 1810 and supports the notion that the in vitro observed postantibiotic effect of this spectinamide also translates to the in vivo situation in mice. The obtained results and pharmacometric model for the exposure-response relationship of 1810 provide a rational basis for dose selection in future efficacy studies of this compound against M. tuberculosis.

Rupatadine Oral Solution Titration by Body Weight in Paediatric Patients Suffering from Allergic Rhinitis: A Population Pharmacokinetic Study

Background

Allergic rhinitis (AR) and chronic urticaria, both are treated in children with doses of second generation of antihistamines that have been mostly based on extrapolation of data obtained in adults. The objectives of this work were to develop a model to explain the pharmacokinetics (PK) of rupatadine, a second generation antihistamine, administered to children 2−11 years old and to calculate the non-compartmental PK parameters for two groups of age (2–5 and 6–11 years old) based on the individual Bayesian estimates from the selected model.

Methods

Data from two PK studies with rupatadine oral solution (1 mg/mL) were pooled: Study A, an extensive blood sampling study performed in 11 children (6–11 years old) who received a single oral dose of rupatadine; and Study B, a sparse blood sampling study in 40 children (2–5 years old) receiving multiple oral doses. A simultaneous population PK model was developed using data available for all children. Using individual Bayesian estimates from the selected model, steady-state plasma concentrations for both studies were simulated and the non-parametric PK parameters were calculated for two age groups: 2–5 years (subgroup I) and 6–11 years (subgroup II).

Results

A two-compartment model with first-order absorption and elimination with clearance depending on body weight, better described the PK of rupatadine for 2–11 year old children. The plasma clearance dependence on weight was linear. The mean (SD) non-compartment PK parameters calculated using simulated plasma profiles at steady state were: C_max, 2.54 (1.26) vs 1.96 (0.52) ng/mL; AUC_0-24h, 10.74 (3.09) vs 10.38 (4.31) ng/mL/h; and t_1/2, 12.28 (3.09) vs 15.94 (4.09) h, for children 6–11 and 2–5 years old, respectively.

Conclusions

The PK of rupatadine depends on the weight of paediatric patients but not on their age. The dosage strategy adjusted by body weight in children 2–11 years old (2.5 mL if weight 10–25 kg, and 5 mL if ≥ 25 kg) provides similar exposure between the two groups of age, and to that obtained in adults with the 10 mg dose tablet formulation.

Pediatric Dose Selection for Therapeutic Proteins

In selecting optimal dosing regimens in support of the clinical use of monoclonal antibodies and other therapeutic proteins in pediatric indications, the unique pharmacokinetic properties of this class of biologics, as well as the underlying physiologic and pathophysiologic processes and their modulation by childhood growth and development, needs to be appreciated. During drug development, first-in-pediatric dose selection is a capstone event in the pediatric investigation plan that relies heavily on extrapolation of pharmacokinetic and pharmacodynamic data from adult to pediatric populations. It is facilitated by combinations of pharmacometric approaches, including allometry, physiologically based pharmacokinetic modeling, and population pharmacokinetic analyses, although data on reliability and qualification of some of these tools in the context of therapeutic proteins are still limited but emerging. Presented data suggest nonlinear relationships between body weight and both clearance and volume of distribution for therapeutic proteins in pediatric populations, with allometric exponents of 0.75 and 0.8, respectively. For newborns and infants (<1 year), even higher nonlinearity seems to occur. Translation of the quantitative characterization of the pediatric pharmacokinetics of therapeutic proteins into dosing regimens for the drug label requires compromising between precision dosing and clinical practicability, with tiered dosing algorithms based on size or age strata being the currently most frequently applied methodology.

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.

Antibiotic Resistance in Pediatric Infections: Global Emerging Threats, Predicting the Near Future

Antibiotic resistance is a public health threat of the utmost importance, especially when it comes to children: according to WHO data, infections caused by multidrug resistant bacteria produce 700,000 deaths across all ages, of which around 200,000 are newborns. This surging issue has multipronged roots that are specific to the pediatric age. For instance, the problematic overuse and misuse of antibiotics (for wrong diagnoses and indications, or at wrong dosage) is also fueled by the lack of pediatric-specific data and trials. The ever-evolving nature of this age group also poses another issue: the partly age-dependent changes of a developing system of cytochromes determine a rather diverse population in terms of biochemical characteristics and pharmacokinetics profiles, hard to easily codify in an age- or weight-dependent dosage. The pediatric population is also penalized by the contraindications of tetracyclines and fluoroquinolones, and by congenital malformations which often require repeated hospitalizations and pharmacological and surgical treatments from a very young age. Emerging threats for the pediatric age are MRSA, VRSA, ESBL-producing Enterobacteriaceae, carbapenem-resistant Enterobacteriaceae and the alarming colistin resistance. Urgent actions need to be taken in order to step back from a now likely post-antibiotic era, where simple infections might cause infant death once again.

Use of Real-World Data and Pharmacometric Modeling in Support of Lacosamide Dosing in Pediatric Patients Under 4 Years of Age

The antiepileptic drug lacosamide (LCM) is approved in the United States and the European Union as monotherapy as well as adjunctive therapy for the treatment of focal seizures in children ≥4 years of age and adults. Using real-world therapeutic drug monitoring data, we performed a pharmacometric analysis for 315 pediatric patients (>1 month to <18 years of age) who received lacosamide as both monotherapy and adjunctive therapy. Population pharmacokinetic modeling was performed using nonlinear mixed-effects modeling with a 1-compartment structural model with linear elimination, where clearance and volume of distribution were allometrically scaled for body weight, with no further need for age-associated maturation functions. A covariate analysis for age, sex, race, and coadministration of other antiepileptic drugs identified phenobarbital and felbamate to significantly increase lacosamide clearance (1.71- and 1.46-fold, respectively). Based on the developed population pharmacokinetic model, simulations were performed in virtual pediatric patients to explore age-associated dose requirements to match lacosamide exposure in patient groups of different age with the exposure achieved in children ≥4 year of age with the weight-based dosing recommendations provided by the US Food and Drug Administration. Based on this approach, our analysis suggested that children ≥3 years of age needed the same dose as recommended by the US Food and Drug Administration for children ≥4 years of age (12 mg/kg/d), while children 1 to 3 years of age may need 13 to 14 mg/kg/d and infants between 1 month and 1 year of age may need 15 to 18 mg/kg/d (based on their actual age) to match the exposure seen in children ≥4 years of age.

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.

Population Pharmacokinetics of Cefoxitin Administered for Pediatric Cardiac Surgery Prophylaxis

BACKGROUND

Available data about pharmacokinetics (PK) of antimicrobials administered as surgical prophylaxis to children undergoing cardiac surgery with cardiopulmonary bypass (CPB) showed that drug concentrations during CPB may be supra or subtherapeutic. The aim of this study was to determine the population PK and pharmacodynamic target attainment (PTA) of cefoxitin during pediatric CPB surgery.

METHODS

A prospective interventional study was conducted. Cefoxitin (40 mg/kg, up to max 1000 mg) was administered before skin incision. Blood samples were obtained in the operatory room throughout surgery. Population PK, PTA, and safety of cefoxitin were evaluated in neonates, infants, children <10 and >10 years old.

RESULTS

Forty patients were enrolled. Cefoxitin levels correlated with time from bolus administration (r = -0.6, P = 0.0001) and, after 240 minutes from bolus, drug values below the target (8 mg/L) were shown. Cefoxitin concentrations were best described by a one-compartment model with first order elimination. A significant relationship was identified between body weight, age, body mass index, and serum creatinine on drug clearance and age, body weight, and body mass index on cefoxitin volume of distribution. The PTA for free drug concentration being above the minimum inhibitory concentration of 8 mg/L for at least 240 minutes was >90% in all age groups except in patients >10 years of age (PTA = 62%).

CONCLUSIONS

Cefoxitin PK appears to be significantly influenced by CPB with generally reduced drug clearance. The PTA was adequately achieved in the majority of patients except in patients >10 years old or longer surgeries.

The pharmacokinetics, safety, and tolerability of mirabegron in children and adolescents with neurogenic detrusor overactivity or idiopathic overactive bladder and development of a population pharmacokinetic model-based pediatric dose estimation

INTRODUCTION

Mirabegron, a selective β3-adrenoreceptor agonist, is a well-established alternative to antimuscarinics in adults with overactive bladder (OAB) symptoms and is under development for use in pediatric patients. Understanding drug pharmacokinetics (PK) in pediatric patients is needed to determine appropriate dosing. Conducting these studies is ethically complex, particularly as regulatory guidance requires that PK is assessed in pediatric patients with a therapeutic need for the drug. It is also vital to evaluate the safety/tolerability and palatability/acceptability of pediatric formulations.

PURPOSE

The purpose of the study was to characterize the PK of mirabegron in pediatric patients with neurogenic detrusor overactivity or idiopathic OAB, to provide a basis for a weight-based dosing algorithm, and to evaluate the safety, tolerability, and palatability/acceptability of the formulations.

MATERIALS AND METHODS

A preliminary population PK model constructed from adult data with allometric scaling was used to predict single weight-adjusted mirabegron doses. This was developed to achieve exposures in pediatric patients in two phase 1 studies that were consistent with steady state in adults following once-daily 25 mg ('low dose') and 50 mg ('high dose') dosing. In study 1, adolescents (12-<18 years) and children (5-<12 years) received a single tablet under fed or fasted conditions. In study 2, children (3-<12 years) received a single oral suspension dose under fed conditions. The PK data were used to assess the predictive value of the preliminary PK model and to update it to analyze mirabegron PK in pediatric patients. The safety/tolerability and palatability/acceptability of the formulations were evaluated.

RESULTS

Forty-three patients comprised six study cohorts: adolescents, low-dose tablets, fed (n = 7); children, low-dose tablets, fed (n = 7); adolescents, high-dose tablets, fed (n = 8); children, high-dose tablets, fed (n = 6); children, high-dose tablets, fasted (n = 6); and children, high-dose oral suspension, fed (n = 9). The population PK model-based doses for tablets and oral suspension achieved exposures that were typically consistent with steady state in adults. The final population PK model was used to describe the PK for mirabegron in pediatric patients (Table). Both formulations were well tolerated, and there were no reports of bad taste or swallowing difficulties for the tablets, although some found the oral suspension unpleasant.

CONCLUSIONS

The single, weight-adjusted pediatric mirabegron doses were successfully predicted by population PK modeling to achieve drug exposures comparable with steady state in adults. The finalized PK model used to characterize the pediatric PK of mirabegron will be utilized to develop a weight-based dosing algorithm. The single mirabegron doses were well tolerated.

Use of Model-Based Compartmental Analysis and a Super-Child Design to Study Whole-Body Retinol Kinetics and Vitamin A Total Body Stores in Children from 3 Lower-Income Countries

BACKGROUND

Model-based compartmental analysis has been used to describe and quantify whole-body vitamin A metabolism and estimate total body stores (TBS) in animals and humans.

OBJECTIVES

We applied compartmental modeling and a super-child design to estimate retinol kinetic parameters and TBS for young children in Bangladesh, Guatemala, and the Philippines.

METHODS

Children ingested [13C10]retinyl acetate and 1 or 2 blood samples were collected from each child from 6 h to 28 d after dosing. Temporal data for fraction of dose in plasma [13C10]retinol were modeled using WinSAAM software and a 6-component model with vitamin A intake included as weighted data.

RESULTS

Model-predicted TBS was 198, 533, and 1062 μmol for the Bangladeshi (age, 9-17 mo), Filipino (12-18 mo), and Guatemalan children (35-65 mo). Retinol kinetics were similar for Filipino and Guatemalan groups and generally faster for Bangladeshi children, although fractional transfer of plasma retinol to a larger exchangeable storage pool was the same for the 3 groups. Recycling to plasma from that pool was ∼2.5 times faster in the Bangladeshi children compared with the other groups and the recycling number was 2-3 times greater. Differences in kinetics between groups are likely related to differences in vitamin A stores and intakes (geometric means: 352, 727, and 764 μg retinol activity equivalents/d for the Bangladeshi, Filipino, and Guatemalan children, respectively).

CONCLUSIONS

By collecting 1 or 2 blood samples from each child to generate a composite plasma tracer data set with a minimum of 5 children/time, group TBS and retinol kinetics can be estimated in children by compartmental analysis; inclusion of vitamin A intake data increases confidence in model predictions. The super-child modeling approach is an effective technique for comparing vitamin A status among children from different populations. These trials were registered at www.clinicaltrials.gov as NCT03000543 (Bangladesh), NCT03345147 (Guatemala), and NCT03030339 (Philippines).

Pharmacokinetics and Clinical Pharmacology of Monoclonal Antibodies in Pediatric Patients

Monoclonal antibodies (mAbs) and their derivatives are increasingly used in pediatric pharmacotherapy, and the number of antibody-based drug products with approved pediatric indications is continuously growing. In most instances, pediatric use is being pursued after the efficacy and safety of novel antibody medications have been established in adult indications. The pediatric extrapolation exercise that is frequently used in this context to bridge efficacy and safety from adults to children is oftentimes challenged through uncertainties and knowledge gaps in how to reliably extrapolate pharmacokinetics and clinical pharmacology of mAbs to different pediatric age groups, and how to derive age-appropriate dosing regimens that strike a balance between precision dosing and practicability. The article highlights some of the pharmacokinetic and clinical pharmacology challenges with regard to therapeutic use of mAbs and antibody derivatives in children, including immunogenicity events. Although considering body size-based differences in drug disposition can account for many of the perceived and actual differences in the distribution and elimination of antibody-based therapeutics between children and adults, increasing evidence suggests potential or actual age-associated differences beyond size differences, especially for young pediatric patients such as newborns and infants. To overcome age-associated differences in antibody disposition, various different dosing approaches have been applied to ensure safe and efficacious antibody exposure for pediatric populations of different ages. The development of such dosing regimens and the associated pathway to pediatric indication approval is illustrated in more detail for two antibody-based biologics, the fusion protein abatacept and the mAb tocilizumab.

Antiplatelet Effect of a Pulaimab [Anti-GPIIb/IIIa F(ab)2 Injection] Evaluated by a Population Pharmacokinetic-pharmacodynamic Model

BACKGROUND

Cardiovascular disease has one of the highest mortality rates among all the diseases. Platelets play an important role in the pathogenesis of cardiovascular diseases. Platelet membrane glycoprotein GPIIb/IIIa antagonists are the most effective antiplatelet drugs, and pulaimab is one of these. The study aims to promote individual medication of pulaimab [anti-GPIIb/IIIa F(ab)2 injection] by discovering the pharmacological relationship among the dose, concentration, and effects. The goal of this study is to establish a population pharmacokineticpharmacodynamic model to evaluate the antiplatelet effect of intravenous pulaimab injection.

METHODS

Data were collected from 59 healthy subjects who participated in a Phase-I clinical trial. Plasma concentration was used as the pharmacokinetic index, and platelet aggregation inhibition rate was used as the pharmacodynamic index. The basic pharmacokinetics model was a two-compartment model, whereas the basic pharmacodynamics model was a sigmoid-EMAX model with a direct effect. The covariable model was established by a stepwise method. The final model was verified by a goodness-of-fit method, and predictive performance was assessed by a Bootstrap (BS) method.

RESULTS

In the final model, typical population values of the parameters were as follows: central distribution Volume (V1), 183 L; peripheral distribution Volume (V2), 349 L; Central Clearance (CL), 31 L/h; peripheral clearance(Q), 204 L/h; effect compartment concentration reaching half of the maximum effect (EC50), 0.252 mg/L; maximum effect value (EMAX), 54.0%; and shape factor (γ), 0.42. In the covariable model, thrombin time had significant effects on CL and EMAX. Verification by the goodness-of-fit and BS methods showed that the final model was stable and reliable.

CONCLUSION

A model was successfully established to evaluate the antiplatelet effect of intravenous pulaimab injection that could provide support for the clinical therapeutic regimen.

Addition of Vitamin A Intake Data during Compartmental Modeling of Retinol Kinetics in Theoretical Humans Leads to Accurate Prediction of Vitamin A Total Body Stores and Kinetic Parameters in Studies of Reasonable Duration

BACKGROUND

Mathematical modeling of theoretical data has been used to validate experimental protocols and methods in several fields.

OBJECTIVES

We hypothesized that adding dietary vitamin A intake data as an input during compartmental modeling of retinol kinetics would lead to accurate prediction of vitamin A total body stores (TBS) at 2 specified study lengths and would reduce study duration required to accurately define the system.

METHODS

We generated reference values for state variables (including TBS and intake) and kinetic parameters for 12 theoretical individuals (4 each of children, younger adults, and older adults) based on modeling plasma retinol tracer data for 365 d. We compared TBS predictions using data to 28 d (children) or 56 d (adults) without and with intake included in the model to reference values for each subject. Then, by truncating data sets from 365 d, we determined the shortest study duration required to accurately define the system without and with inclusion of vitamin A intake.

RESULTS

Reference values for TBS ranged from 30 to 3023 µmol. Study durations of 28 and 56 d were sufficient to accurately predict TBS for 6 of the 12 subjects without intake; adding intake resulted in accurate predictions of TBS for all individuals. When intake was not included as a modeling input, durations of 35-310 d were required to define the system; inclusion of intake data substantially reduced the time required to 10-42 d.

CONCLUSIONS

Inclusion of vitamin A intake as additional data input when modeling vitamin A kinetics allows investigators to accurately predict TBS and define the vitamin A system in studies of reasonable length (4 wk in children and 8 wk in adults). Because it is generally possible to obtain estimates/measures of intake, including such data increases confidence in model predictions while also making studies more feasible.

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.

Population Pharmacokinetics and Pharmacodynamics of Dexmedetomidine in Children Undergoing Ambulatory Surgery

BACKGROUND

Dexmedetomidine (DEX) is an α-2 adrenergic agonist with sedative and analgesic properties. Although not approved for pediatric use by the Food and Drug Administration, DEX is increasingly used in pediatric anesthesia and critical care. However, very limited information is available regarding the pharmacokinetics of DEX in children. The aim of this study was to investigate DEX pharmacokinetics and pharmacodynamics (PK-PD) in Mexican children 2-18 years of age who were undergoing outpatient surgical procedures.

METHODS

Thirty children 2-18 years of age with American Society of Anesthesiologists physical status score of I/II were enrolled in this study. DEX (0.7 μg/kg) was administered as a single-dose intravenous infusion. Venous blood samples were collected, and plasma DEX concentrations were analyzed with a combination of high-performance liquid chromatography and electrospray ionization-tandem mass spectrometry. Population PK-PD models were constructed using the Monolix program.

RESULTS

A 2-compartment model adequately described the concentration-time relationship. The parameters were standardized for a body weight of 70 kg by using an allometric model. Population parameters estimates were as follows: mean (between-subject variability): clearance (Cl) (L/h × 70 kg) = 20.8 (27%); central volume of distribution (V1) (L × 70 kg) = 21.9 (20%); peripheral volume of distribution (V2) (L × 70 kg) = 81.2 (21%); and intercompartmental clearance (Q) (L/h × 70 kg) = 75.8 (25%). The PK-PD model predicted a maximum mean arterial blood pressure reduction of 45% with an IC50 of 0.501 ng/ml, and a maximum heart rate reduction of 28.9% with an IC50 of 0.552 ng/ml.

CONCLUSIONS

Our results suggest that in Mexican children 2-18 years of age with American Society of Anesthesiologists score of I/II, the DEX dose should be adjusted in accordance with lower DEX clearance.

Population pharmacokinetics of trabectedin in adolescent patients with cancer

PURPOSE

To characterize the trabectedin population pharmacokinetics in children and adolescent patients with cancer and compare it with the trabectedin pharmacokinetics in adults.

METHODS

Plasma concentrations from ten adolescent and three children with cancer (age range 4.0-17.0 years) treated with trabectedin at doses ranging from 1.1 to 1.7 mg/m², administered as a 24-h continuous intravenous infusion every 3 weeks, were available for the analysis. An external model evaluation was performed to verify whether a previously developed adult population pharmacokinetic model was predictive of the pediatric plasma concentrations of trabectedin. The maximum a posteriori estimation of the individual pharmacokinetic parameters for pediatric patients was conducted, after successful completion of the external evaluation step. The relationships between pharmacokinetic parameters and body size were evaluated.

RESULTS

External evaluation methods showed no major differences between the adult population and children and adolescent patients of this study. The mean ± standard deviation (SD) of the individual estimated clearance and central volume of distribution in these children/adolescent patients was 36.4 ± 16.1 L/h and 13.2 ± 6.54 L, respectively. These values were similar to the typical values reported for adult patients-37.6 L/h and 13.9 L (for females) and 16.1 L (for males). The median area under the plasma concentration versus time curve (AUC) in children/adolescent patients was 55.1 µg h/L, while in the adult population the median AUC was 61.3 µg h/L, both administered a 1.5 mg/m² dose regimen with mean (range) BSA for adults = 1.86 (0.90-2.80) vs children/adolescent patients = 1.49 (0.66-2.54).

CONCLUSIONS

The adult population pharmacokinetic model adequately described the trabectedin plasma concentrations and its variability in the pediatric population of patients involved in this assessment that mostly comprised adolescents. The trabectedin systemic exposure achieved in this population was comparable (within 12%) to the exposure obtained in adult population when the same dose, expressed in mg/m², was administered.

[Pharmacokinetic and pharmacodynamic modelling in paediatric drug development with a focus on physiology-based pharmacokinetic simulations]

Pharmacokinetic and pharmacodynamic models are mandatory for dosing and the safe use of drugs in the paediatric population. Different modelling methods allow for the development of dosing regimens for children requiring only a small number of blood samples or none at all. The medicines regulatory authorities recommend using these methods for paediatric drug development programs. Taking sildenafil as an example, the least invasive method of physiology-based pharmacokinetic simulation for the development of dosing regimens in the paediatric population is presented.

Population Pharmacokinetics and Pharmacodynamics of Bezlotoxumab in Adults with Primary and Recurrent Clostridium difficile Infection

The fully human monoclonal antibody bezlotoxumab is indicated for preventing the recurrence of Clostridioides difficile (formerly Clostridium difficile) infection (CDI) in adults who receive antibacterial treatment for CDI and who are at high risk for a CDI recurrence. The efficacy and safety of 10-mg/kg of body weight bezlotoxumab were demonstrated in two phase 3 trials: the MODIFY I (ClinicalTrials.gov registration number NCT01241552) and MODIFY II (ClinicalTrials.gov registration number NCT01513239) trials. Here, a population pharmacokinetic (popPK) analysis, performed using data from the MODIFY I and II trials (n = 1,515) and from three phase 1 trials (n = 72) to characterize bezlotoxumab pharmacokinetics (PK) in phase 3 clinical trial participants and in healthy subjects, is reported. A stepwise covariate search was conducted to identify factors influencing PK. Post hoc-estimated bezlotoxumab exposures from the popPK model were used to conduct an exposure-response analysis for CDI recurrence. Bezlotoxumab PK were described by a two-compartment model with linear elimination and allometric scaling for clearance and the volume of distribution by body weight. Although the final popPK model included gender, ethnicity (Japanese descent), race (black versus nonblack), and albumin level as significant covariates, the impact of these factors was not clinically meaningful, based on the totality of the PK and clinical experience. Exposure-response analysis of CDI recurrence demonstrated a similar low rate of CDI recurrence over the entire range of exposures achieved in the phase 3 trials, indicating that exposures were on the maximal response plateau of the exposure-response curve. Overall, the analyses confirmed the appropriateness of the 10-mg/kg dose across the phase 3 trial population with no dose adjustments necessary over a broad demographic background.

Model-informed pediatric development applied to bilastine: Analysis of the clinical PK data and confirmation of the dose selected for the target population

Bilastine is a non-sedating second-generation H₁ antihistamine approved for treatment of allergic rhinoconjunctivitis (AR) and urticaria (U) in adults at the oral (p.o.) dose of 20 mg once daily (OD). Optimal attributes can be anticipated for its clinical use in pediatrics due to its favorable safety and tolerability and age-independent PD profile. The aim of this work was to characterize bilastine PK in children through population modeling of data from a limited sampling confirmatory clinical trial in children with AR or U. The objective was also to ascertain whether the proposed dose (10 mg/day) in the target pediatric subset aged 2-<12 years matches the systemic exposure seen in adults at the 20 mg/day dose. A popPK model characterizing bilastine PK behavior in children aged from 4 to <12 years treated with 10 mg oral bilastine daily was successfully developed and qualified. No relationship was found between bilastine PK and age or weight; stopping rules pre-stablished to finalize the trial, i.e., model completeness and no dependence of exposure on decreasing age, were thus fulfilled. On a second step, the popPK model in children was linked to the PD model in adults assuming the same PD as described in adults and used to compare the PD outcome between both populations. Finally, an allometric scaling method and a physiological approximation were used to evaluate the suitability of the selected dose in the youngest children, showing that children from 2 years were deemed to belong to the same population as well. The achievement of comparable PK (i.e., within the range) to that observed in adults after the therapeutic dose of 20 mg, together with the achievement of similar PD and additional integrative analysis, served to confirm the validity of the 10 mg daily dose for the target pediatric subset (2 to <12 years).

A Population-Based (Super-Child) Approach for Predicting Vitamin A Total Body Stores and Retinol Kinetics in Children Is Validated by the Application of Model-Based Compartmental Analysis to Theoretical Data

BACKGROUND

Public health nutritionists need accurate and feasible methods to assess vitamin A status and to evaluate efficacy of interventions, especially in children. The application of population-based designs to tracer kinetic data is an effective approach that reduces sample burden for each child.

OBJECTIVES

Objectives of the study were to use theoretical data to validate a population-based (super-child) approach for estimating group mean vitamin A total body stores (TBS) and retinol kinetics in children and to use population-based data to improve individual TBS predictions using retinol isotope dilution (RID).

METHODS

We generated plasma retinol kinetic data from 6 h to 56 d for 50 theoretical children with high vitamin A intakes, assigning values within physiologically reasonable ranges for state variables and kinetic parameters ("known values"). Mean data sets for all subjects at extensive (n = 36) and reduced (n = 11) sampling times, plus 5 data sets for reduced numbers (5/time, except all at 4 d) and times, were analyzed using Simulation, Analysis and Modeling software. Results were compared with known values; population RID coefficients were used to calculate TBS for individuals.

RESULTS

For extensive and reduced data sets including all subjects, population TBS predictions were within 1% of the known value. For 5 data sets reflecting numbers and times being used in ongoing super-child studies, predictions were within 1-17% of the known group value. Using RID equation coefficients from population modeling, TBS predictions at 4 d were within 25% of the known value for 66-80% of subjects and reflected the range of assigned values; when ranked, predicted and assigned values were significantly correlated (Rs  = 0.93, P < 0.0001). Results indicate that 7 d may be better than 4 d for applying RID in children. For all data sets, predictions for kinetic parameters reflected the range of known values.

CONCLUSION

The population-based (super-child) approach provides a feasible experimental design for quantifying retinol kinetics, accurately estimating group mean TBS, and predicting TBS for individuals reasonably well.

Population Pharmacokinetics of Gentamicin in Neonates with Hypoxemic-Ischemic Encephalopathy Receiving Controlled Hypothermia

OBJECTIVE

Identify population pharmacokinetics and pharmacodynamic target attainment of gentamicin in neonates with hypoxic-ischemic encephalopathy (HIE) undergoing controlled hypothermia (CH).

DESIGN

Prospective open-label pharmacokinetic study. Gentamicin concentrations were modeled and dosing regimens simulated for a 5000-patient neonatal population with HIE receiving CH using PMetrics, a nonparametric, pharmacometric modeling, and simulation package for R.

SETTING

A 189-bed children's tertiary care teaching hospital.

RESULTS

Twelve patients, 5 (42%) females and 7 (58%) males, met inclusion criteria with a median gestation age of 39.9 weeks (interquartile range [IQR] 38.5-40.2 wks) and a median birthweight (BW) of 3.3 kg (IQR 3.1-3.7 kg). Gentamicin concentrations were best described by a two-compartment model with first-order elimination with BW as a covariate on volume of distribution (Vd). The mean total body population clearance (CL) was 2.2 ± 0.7 ml/minute/kg, and the volume of the central compartment was 0.44 ± 0.06 L/kg. The R² , bias, and precision for the observed versus population predicted model were 0.917, 1.15, and 10.9 μg/ml; the R² , bias, and precision for the observed versus individual predicted model were 0.982, -0.132, and 0.932 μg/ml, respectively. The calculated mean population estimate for the total Vd was 0.96 ± 0.4 L/kg. The dosing regimen that most consistently produced a maximum concentration (C_max ) in the range of 10-12 mg/L with a minimum concentration (C_min ) level less than 2 mg/L was 5 mg/kg/dose given every 36 hours.

CONCLUSION

These data suggest the population pharmacokinetics of gentamicin in neonates with HIE receiving CH have an increase in gentamicin CL and are different from previous reports in neonates with HIE not receiving CH and/or neonates without HIE. This analysis suggests a dosing regimen of 5 mg/kg/dose every 36 hours results in a gentamicin C_max within the range of 10-12 mg/L with a C_min lower than 2 mg/L, which is appropriate for treating susceptible gram-negative organisms with minimum inhibitory concentrations of 1 mg/L or lower.

Evaluation of the Adequacy of WHO Revised Dosages of the First-Line Antituberculosis Drugs in Children with Tuberculosis Using Population Pharmacokinetic Modeling and Simulations

Optimal doses for antituberculosis (anti-TB) drugs in children have yet to be established. In 2010, the World Health Organization (WHO) recommended revised dosages of the first-line anti-TB drugs for children. Pharmacokinetic (PK) studies that investigated the adequacy of the WHO revised dosages to date have yielded conflicting results. We performed population PK modeling using data from one of these studies to identify optimal dosage ranges. Ghanaian children with tuberculosis on recommended therapy with rifampin (RIF), isoniazid (INH), pyrazinamide (PZA), and ethambutol (EMB) for at least 4 weeks had blood samples collected predose and at 1, 2, 4, and 8 hours postdose. Drug concentrations were determined by validated liquid chromatography-mass spectrometry methods. Nonlinear mixed-effects models were applied to describe the population PK of those drugs using MonolixSuite2016R1 (Lixoft, France). Bayesian estimation was performed, the correlation coefficient, bias, and precision between the observed and predicted areas under the concentration-time curve (AUCs) were calculated, and Bland-Altman plots were analyzed. The population PK of RIF and PZA was described by a one-compartment model and that for INH and EMB by a two-compartment model. Plasma maximum concentration (C_max) and AUC targets were based on published results for children from India. The lowest target values for pediatric TB patients were attainable at the WHO-recommended dosage schedule for RIF and INH, except for N-acetyltransferase 2 non-slow acetylators (rapid and intermediate acetylators) in the lower-weight bands. However, higher published adult targets were not attainable for RIF and INH. The targets were not achieved for PZA and EMB. (This study has been registered at ClinicalTrials.gov under identifier NCT01687504.).

Model Informed Pediatric Development Applied to Bilastine: Ontogenic PK Model Development, Dose Selection for First Time in Children and PK Study Design

PURPOSE

Bilastine is an H₁ antagonist whose pharmacokinetics (PK) and pharmacodynamics (PD) have been resolved in adults with a therapeutic oral dose of 20 mg/day. Bilastine has favorable characteristics for use in pediatrics but the PK/PD and the optimal dose in children had yet to be clinically explored. The purpose is to: (1) Develop an ontogenic predictive model of bilastine PK linked to the PD in adults by integrating current knowledge; (2) Use the model to design a PK study in children; (3) Confirm the selected dose and the study design through the evaluation of model predictability in the first recruited children; (4) Consider for inclusion the group of younger children (< 6 years).

METHODS

A semi-mechanistic approach was applied to predict bilastine PK in children assuming the same PD as described in adults. The model was used to simulate the time evolution of plasma levels and wheal and flare effects after several doses and design an adaptive PK trial in children that was then confirmed using data from the first recruits by comparing observations with model predictions.

RESULTS

PK/PD simulations supported the selection of 10 mg/day in 2 to <12 year olds. Results from the first interim analysis confirmed the model predictions and design hence trial continuation.

CONCLUSION

The model successfully predicted bilastine PK in pediatrics and optimally assisted the selection of the dose and sampling scheme for the trial in children. The selected dose was considered suitable for younger children and the forthcoming safety study in children aged 2 to <12 years.

Pharmacokinetics and Model-Based Dosing to Optimize Fludarabine Therapy in Pediatric Hematopoietic Cell Transplant Recipients

A prospective multicenter study was conducted to characterize the pharmacokinetics (PK) and pharmacodynamics (PD) of fludarabine plasma (f-ara-a) and intracellular triphosphate (f-ara-ATP) in children undergoing hematopoietic cell transplantation (HCT) and receiving fludarabine with conditioning. Plasma and peripheral blood mononuclear cells (PBMCs) were collected over the course of therapy for quantitation of f-ara-a and f-ara-ATP. Nonlinear mixed-effects modeling was used to develop the PK model, including identification of covariates impacting drug disposition. Data from a total of 133 children (median age, 5 years; range, .2 to 17.9) undergoing HCT for a variety of malignant and nonmalignant disorders were available for PK-PD modeling. The implementation of allometric scaling of PK parameters alone was insufficient to describe drug clearance, particularly in very young children. Renal impairment was predicted to increase drug exposure across all ages. The rate of f-ara-a entry into PBMCs (expressed in pmoles per million cells) decreased over the course of therapy, resulting in 78% lower f-ara-ATP after the fourth dose (1.7 pmoles/million cells [range, .2 to 7.2]) compared with first dose (7.9 pmoles/million cells [range, .7 to 18.2]). The overall incidence of treatment-related mortality (TRM) was low at 3% and 8% at days 60 and 360, respectively, and no association with f-ara-a exposure and TRM was found. In the setting of malignancy, disease-free survival was highest at 1 year after HCT in subjects achieving a systemic f-ara-a cumulative area under the curve (cAUC) greater than 15 mg*hour/L compared to patients with a cAUC less than 15 mg*hour/L (82.6% versus 52.8% P = .04). These results suggest that individualized model-based dosing of fludarabine in infants and young children may reduce morbidity and mortality through improved rates of disease-free survival and limiting drug-related toxicity. ClinicalTrials.gov Identifier: NCT01316549.

Population Pharmacokinetics and Pharmacodynamic Target Attainment of Vancomycin in Neonates on Extracorporeal Life Support

OBJECTIVES

To evaluate the population pharmacokinetics and pharmacodynamic target attainment of vancomycin in neonates with a contemporary ¼-inch extracorporeal life support circuit with a Quadrox-iD Pediatric oxygenator (Maquet Cardiovascular, LLC, Wayne, NJ).

DESIGN

Retrospective medical record review.

SETTING

Two free-standing tertiary/quaternary pediatric children's hospitals.

PATIENTS

Neonates receiving either veno-arterial or veno-venous extracorporeal life support and vancomycin for empiric or definitive therapy with resulting serum concentrations.

INTERVENTIONS

None.

MEASUREMENTS AND MAIN RESULTS

Twelve patients with a median gestations age of 39 weeks (range 36-41 wk) and a median postnatal age of 9.5 days (range 0-28 d) accounted for 14 courses of vancomycin therapy while on extracorporeal life support and were included in the analysis. The median weight was 3.1 kg (range 2.2-4.41 kg) with five of 12 patients (41.7%) being female. Vancomycin concentrations were best described by an one-compartment model incorporating allometric scaling of estimated glomerular filtration rate on clearance. The mean total body clearance (mL/min/kg) for the population was 3.48 ± 1.31 mL/min/kg, and the mean total volume of distribution (L/kg) for the population was 1.2 ± 0.4 L/kg. The intermittent and continuous infusion dosing regimens that provided for the highest percentage of trough concentrations in the range of 10-20 mg/L were the 10 mg/kg/dose IV q8h, 12.5 mg/kg/dose IV q8-12h, 15 mg/kg/dose IV q12h, and 20 mg/kg/dose IV q12h, and the 20, 25, and 30 mg/kg/d continuous infusion regimens, respectively. All regimens allowed for an area under the concentration:minimum inhibitory concentration ratio of 400:1 for minimum inhibitory concentrations of less than or equal to 0.5 mg/L for a 90% PTA. None of the simulated regimens had a greater than 90% probability of achieving an area under the concentration:minimum inhibitory concentration ratio of 400:1 for vancomycin minimum inhibitory concentrations greater than or equal to 1 mg/L while maintaining trough concentrations in the range of 10-20 mg/L.

CONCLUSIONS

To our knowledge, this is the first pharmacokinetic and pharmacodynamic study of neonates receiving vancomycin with a contemporary ¼-inch extracorporeal life support circuit including the Quadrox-iD Pediatric oxygenator (Maquet Cardiovascular, LLC). The data suggest differences in vancomycin pharmacokinetics compared with previous extracorporeal life support data, notably a more rapid clearance, which could result in lower vancomycin concentrations. Considering this, a more aggressive initial dosing regimen may need to be employed in infants on extracorporeal life support.

Potential Sources of Inter-Subject Variability in Monoclonal Antibody Pharmacokinetics

Understanding inter-subject variability in drug pharmacokinetics and pharmacodynamics is important to ensure that all patients attain suitable drug exposure to achieve efficacy and avoid toxicity. Inter-subject variability in the pharmacokinetics of therapeutic monoclonal antibodies (mAbs) is generally moderate to high; however, the factors responsible for the high inter-subject variability have not been comprehensively reviewed. In this review, the extent of inter-subject variability for mAb pharmacokinetics is presented and potential factors contributing to this variability are explored and summarised. Disease status, age, sex, ethnicity, body size, genetic polymorphisms, concomitant medication, co-morbidities, immune status and multiple other patient-specific details have been considered. The inter-subject variability for mAb pharmacokinetics most likely depends on the complex interplay of multiple factors. However, studies aimed at investigating the reasons for the inter-subject variability are sparse. Population pharmacokinetic models and physiologically based pharmacokinetic models are useful tools to identify important covariates, aiding in the understanding of factors contributing to inter-subject variability. Further understanding of inter-subject variability in pharmacokinetics should aid in development of dosing regimens that are more appropriate.

Population Pharmacokinetics and Pharmacodynamic Target Attainment of Meropenem in Critically Ill Young Children

OBJECTIVE

This study aims to describe the population pharmacokinetics and pharmacodynamic target attainment of meropenem in critically ill children.

METHODS

The study involved a retrospective medical record review from a 189-bed, freestanding children's tertiary care teaching hospital of patients ages 1 to 9 years who received meropenem with concurrent therapeutic drug monitoring.

RESULTS

There were 9 patients ages 1 to 9 years (mean age, 3.1 ± 2.9 years) with a mean weight of 17.1 ± 11.9 kg who met the inclusion/exclusion criteria and were included in the pharmacokinetic analysis. Meropenem concentrations were best described by a 2-compartment model with first-order elimination, with an R² and bias of 0.91 and 13.2 mg/L, respectively, for the observed versus population predicted concentrations, and an R², bias, and imprecision of 1, 0.0675, and 1 mg/L, respectively, for the observed versus individual predicted concentrations. The mean total body drug clearance for the population was 6.99 ± 2.5 mL/min/kg, and V_c was 0.57 ± 0.47 L/kg. The calculated population estimate for the total volume of distribution was 0.78 ± 0.73 L/kg. Standard 0.5-hour meropenem infusions did not provide for appropriate pharmacodynamic exposures of 40% free time > minimum inhibitory concentration (40% fT > MIC) for Gram-negative organisms with susceptible MICs. Dosage regimens employing prolonged and continuous infusion regimens did provide appropriate pharmacodynamic exposures of 40% fT > MIC for Gram-negative organisms up to the break point for Pseudomonas aeruginosa of 4 mg/L.

CONCLUSION

These data suggest the reference dosage regimens for meropenem (20-40 mg/kg per dose every 8 hours) do not meet an appropriate pharmacodynamic target attainment in critically ill children ages 1 to 9 years. Based on these data, only the 3- to 4-hour prolonged infusion and 24-hour continuous infusion regimens were able to achieve an optimal probability of target attainment against all susceptible Gram-negative bacteria in critically ill children for 40% fT > MIC. Dosage regimens of 120 and 160 mg/kg/day as continuous infusion regimens may be necessary to achieve an optimal probability of target attainment against all susceptible Gram-negative bacteria in critically ill children for 80% fT > MIC. Based on these findings, confirmation with a larger, prospective investigation in critically ill children is warranted.

Data Analysis Protocol for the Development and Evaluation of Population Pharmacokinetic Models for Incorporation Into the Web-Accessible Population Pharmacokinetic Service - Hemophilia (WAPPS-Hemo)

BACKGROUND

Hemophilia is an inherited bleeding disorder caused by a deficiency in a specific clotting factor. This results in spontaneous bleeding episodes and eventual arthropathy. The mainstay of hemophilia treatment is prophylactic replacement of the missing factor, but an optimal regimen remains to be determined. Rather, individualized prophylaxis has been suggested to improve both patient safety and resource utilization. However, uptake of this approach has been hampered by the demanding sampling schedules and complex calculations required to obtain individual estimates of pharmacokinetic (PK) parameters. The use of population pharmacokinetics (PopPK) can alleviate this burden by reducing the number of plasma samples required for accurate estimation, but few tools incorporating this approach are readily available to clinicians.

OBJECTIVE

The Web-accessible Population Pharmacokinetic Service - Hemophilia (WAPPS-Hemo) project aims to bridge this gap by providing a Web-accessible service for the reliable estimation of individual PK parameters from only a few patient samples. This service is predicated on the development of validated brand-specific PopPK models.

METHODS

We describe the data analysis plan for the development and evaluation of each PopPK model to be incorporated into the WAPPS-Hemo platform. The data sources and structure of the dataset are discussed first, followed by the procedures for handling both data below limit of quantification (BLQ) and absence of such BLQ data. Next, we outline the strategies for building the appropriate structural and covariate models, including the possible need for a process algorithm when PK behavior varies between subjects or significant covariates are not provided. Prior to use in a prospective manner, the models will undergo extensive evaluation using a variety of techniques such as diagnostic plots, bootstrap analysis and cross-validation. Finally, we describe the incorporation of a validated PopPK model into the Bayesian post hoc model to produce individualized estimates of PK parameters.

RESULTS

Dense PK data has been collected for more than 20 brands of factor concentrate from both industry-sponsored and investigator-driven studies. The model development process is underway for the majority of molecules, with refinement and validation to be completed in 2017. Further, the WAPPS-Hemo co-investigator network has contributed more than 300 PK assessments for use in model development and evaluation. This constitutes the largest repository of this type of PK data globally.

CONCLUSIONS

The WAPPS-Hemo service aims to eliminate barriers to the uptake of individualized PK-tailored hemophilia treatment. By incorporating this tool into routine practice, clinicians can implement a personalized dosing strategy without performing rigorous sampling or complex calculations. This service is centred on validated models developed according to the robust approach to PopPK modeling described herein.

CLINICALTRIAL

ClinicalTrials.gov NCT02061072; https://clinicaltrials.gov/ct2/show/NCT02061072 (Archived by WebCite at http://www.webcitation.org/6mRIXJh55).

Bottom-up Meets Top-down: Complementary Physiologically Based Pharmacokinetic and Population Pharmacokinetic Modeling for Regulatory Approval of a Dosing Algorithm of Valganciclovir in Very Young Children

Population pharmacokinetic (PopPK) and physiologically based pharmacokinetic (PBPK) models are frequently used to support pediatric drug development. Both methods have strengths and limitations and we used them complementarily to support the regulatory approval of a dosing algorithm for valganciclovir (VGCV) in children <4 months old. An existing pediatric PBPK model was extended to neonates and showed that potential physiological differences compared with older children are minor. The PopPK model was used to simulate ganciclovir (GCV) exposures in children with population typical combinations of body size and renal function and to assess the effectiveness of an alternative dosing algorithm suggested by the US Food and Drug Administration. PBPK and PopPK confirmed that the proposed VGCV dosing algorithm achieves similar GCV exposures in children of all ages and that the alternative dosing algorithm leads to underexposure in a substantial fraction of patients. Our approach raised the confidence in the VGCV dosing algorithm for children <4 months old and supported the regulatory approval.

Development of a population pharmacokinetic model of prucalopride in children with functional constipation

A recent phase 3 trial of prucalopride in children with functional constipation (SPD555-303 ClinicalTrials.gov Identifier: NCT01330381) reported negative efficacy results. Here, we developed a population pharmacokinetic (PK) model of prucalopride in children to assess prucalopride exposure in SPD555-303. An initial population PK model in children was developed based on sampled single-dose data from a phase 1 study (PRU-USA-12). This model was subsequently updated with sampled data from SPD555-303 and used to simulate plasma concentration-time profiles for children aged 6 months to 18 years who were treated once daily with prucalopride 0.02, 0.04, or 0.06 mg kg-1 (maximum dose, 2 mg). Simulated PK profiles were compared with those of adults at the recommended dose of 2 mg once daily. Data were available from 38 patients (median age, 8.5 years) in PRU-USA-12 and 137 patients (median age, 7.9 years) in SPD555-303. Mean (range) area under the plasma concentration-time curve (AUC) at steady state was 62.3 (40.5-82.7) ng mL-1 h (dose, 0.03 mg kg-1) in PRU-USA-12 and 100.3 (22.7-286.0) ng mL-1 h (dose, 0.04 mg kg-1; maximum, 2 mg) in SPD555-303. Prucalopride 0.04 mg kg-1 once daily in children produced similar maximum plasma concentrations and approximately 10% lower AUC compared with adults receiving 2 mg once daily. This population PK analysis indicates that the PK profile of prucalopride in children in SPD555-303 was similar to that observed in adults. The negative efficacy results of SPD555-303 cannot be explained by differences in prucalopride exposure between children and adults.

Children in clinical trials: towards evidence-based pediatric pharmacotherapy using pharmacokinetic-pharmacodynamic modeling

INTRODUCTION

In pediatric pharmacotherapy, many drugs are still used off-label, and their efficacy and safety is not well characterized. Different efficacy and safety profiles in children of varying ages may be anticipated, due to developmental changes occurring across pediatric life.

AREAS COVERED

Beside pharmacokinetic (PK) studies, pharmacodynamic (PD) studies are urgently needed. Validated PKPD models can be used to derive optimal dosing regimens for children of different ages, which can be evaluated in a prospective study before implementation in clinical practice. Strategies should be developed to ensure that formularies update their drug dosing guidelines regularly according to the most recent advances in research, allowing for clinicians to integrate these guidelines in daily practice. Expert commentary: We anticipate a trend towards a systems-level approach in pediatric modeling to optimally use the information gained in pediatric trials. For this approach, properly designed clinical PKPD studies will remain the backbone of pediatric research.

Accurately Achieving Target Busulfan Exposure in Children and Adolescents With Very Limited Sampling and the BestDose Software

BACKGROUND

Busulfan dose adjustment is routinely guided by plasma concentration monitoring using 4-9 blood samples per dose adjustment, but a pharmacometric Bayesian approach could reduce this sample burden.

METHODS

The authors developed a nonparametric population model with Pmetrics. They used it to simulate optimal initial busulfan dosages, and in a blinded manner, they compared dosage adjustments using the model in the BestDose software to dosage adjustments calculated by noncompartmental estimation of area under the time-concentration curve at a national reference laboratory in a cohort of patients not included in model building.

RESULTS

Mean (range) age of the 53 model-building subjects was 7.8 years (0.2-19.0 years) and weight was 26.5 kg (5.6-78.0 kg), similar to nearly 120 validation subjects. There were 16.7 samples (6-26 samples) per subject to build the model. The BestDose cohort was also diverse: 10.2 years (0.25-18 years) and 46.4 kg (5.2-110.9 kg). Mean bias and imprecision of the 1-compartment model-predicted busulfan concentrations were 0.42% and 9.2%, and were similar in the validation cohorts. Initial dosages to achieve average concentrations of 600-900 ng/mL were 1.1 mg/kg (≤12 kg, 67% in the target range) and 1.0 mg/kg (>12 kg, 76% in the target range). Using all 9 concentrations after dose 1 in the Bayesian estimation of dose requirements, the mean (95% confidence interval) bias of BestDose calculations for the third dose was 0.2% (-2.4% to 2.9%, P = 0.85), compared with the standard noncompartmental method based on 9 concentrations. With 1 optimally timed concentration 15 minutes after the infusion (calculated with the authors' novel MMopt algorithm) bias was -9.2% (-16.7% to -1.5%, P = 0.02). With 2 concentrations at 15 minutes and 4 hours bias was only 1.9% (-0.3% to 4.2%, P = 0.08).

CONCLUSIONS

BestDose accurately calculates busulfan intravenous dosage requirements to achieve target plasma exposures in children up to 18 years of age and 110 kg using only 2 blood samples per adjustment compared with 6-9 samples for standard noncompartmental dose calculations.

Pharmacometric Approaches to Personalize Use of Primarily Renally Eliminated Antibiotics in Preterm and Term Neonates

Sepsis remains a major cause of mortality and morbidity in neonates, and, as a consequence, antibiotics are the most frequently prescribed drugs in this vulnerable patient population. Growth and dynamic maturation processes during the first weeks of life result in large inter- and intrasubject variability in the pharmacokinetics (PK) and pharmacodynamics (PD) of antibiotics. In this review we (1) summarize the available population PK data and models for primarily renally eliminated antibiotics, (2) discuss quantitative approaches to account for effects of growth and maturation processes on drug exposure and response, (3) evaluate current dose recommendations, and (4) identify opportunities to further optimize and personalize dosing strategies of these antibiotics in preterm and term neonates. Although population PK models have been developed for several of these drugs, exposure-response relationships of primarily renally eliminated antibiotics in these fragile infants are not well understood, monitoring strategies remain inconsistent, and consensus on optimal, personalized dosing of these drugs in these patients is absent. Tailored PK/PD studies and models are useful to better understand relationships between drug exposures and microbiological or clinical outcomes. Pharmacometric modeling and simulation approaches facilitate quantitative evaluation and optimization of treatment strategies. National and international collaborations and platforms are essential to standardize and harmonize not only studies and models but also monitoring and dosing strategies. Simple bedside decision tools assist clinical pharmacologists and neonatologists in their efforts to fine-tune and personalize the use of primarily renally eliminated antibiotics in term and preterm neonates.

Improving evidence on anticoagulant therapies for venous thromboembolism in children: key challenges and opportunities

Venous thromboembolism (VTE) is increasingly diagnosed in pediatric patients, and anticoagulant use in this population has become common, despite the absence of US Food and Drug Administration (FDA) approval for this indication. Guidelines for the use of anticoagulants in pediatrics are largely extrapolated from large randomized controlled trials (RCTs) in adults, smaller dose-finding and observational studies in children, and expert opinion. The recently FDA-approved direct oral anticoagulants (DOACs), such as dabigatran, rivaroxaban, apixaban, and edoxaban, provide potential advantages over oral vitamin K antagonists and subcutaneous low-molecular-weight heparins (LMWHs). However, key questions arise regarding their potential off-label clinical application in pediatric thromboembolic disease. In this Perspective, we provide background on the use of LMWHs such as enoxaparin as the mainstay of treatment of pediatric provoked VTE; identify key questions and challenges with regard to DOAC trials and future DOAC therapy in pediatric VTE; and discuss applicable lessons learned from the recent pilot/feasibility phase of a large multicenter RCT of anticoagulant duration in pediatric VTE. The challenges and lessons learned present opportunities to improve evidence for anticoagulant therapies in pediatric VTE through future clinical trials.

Population pharmacokinetics of busulfan in pediatric and young adult patients undergoing hematopoietic cell transplant: a model-based dosing algorithm for personalized therapy and implementation into routine clinical use

BACKGROUND

Population pharmacokinetic (PK) studies of busulfan in children have shown that individualized model-based algorithms provide improved targeted busulfan therapy when compared with conventional dose guidelines. The adoption of population PK models into routine clinical practice has been hampered by the tendency of pharmacologists to develop complex models too impractical for clinicians to use. The authors aimed to develop a population PK model for busulfan in children that can reliably achieve therapeutic exposure (concentration at steady state) and implement a simple model-based tool for the initial dosing of busulfan in children undergoing hematopoietic cell transplantation.

PATIENTS AND METHODS

Model development was conducted using retrospective data available in 90 pediatric and young adult patients who had undergone hematopoietic cell transplantation with busulfan conditioning. Busulfan drug levels and potential covariates influencing drug exposure were analyzed using the nonlinear mixed effects modeling software, NONMEM. The final population PK model was implemented into a clinician-friendly Microsoft Excel-based tool and used to recommend initial doses of busulfan in a group of 21 pediatric patients prospectively dosed based on the population PK model.

RESULTS

Modeling of busulfan time-concentration data indicates that busulfan clearance displays nonlinearity in children, decreasing up to approximately 20% between the concentrations of 250-2000 ng/mL. Important patient-specific covariates found to significantly impact busulfan clearance were actual body weight and age. The percentage of individuals achieving a therapeutic concentration at steady state was significantly higher in subjects receiving initial doses based on the population PK model (81%) than in historical controls dosed on conventional guidelines (52%) (P = 0.02).

CONCLUSIONS

When compared with the conventional dosing guidelines, the model-based algorithm demonstrates significant improvement for providing targeted busulfan therapy in children and young adults.

Effect of Kidney Function on Drug Kinetics and Dosing in Neonates, Infants, and Children

Neonates, infants, and children differ from adults in many aspects, not just in age, weight, and body composition. Growth, maturation and environmental factors affect drug kinetics, response and dosing in pediatric patients. Almost 80% of drugs have not been studied in children, and dosing of these drugs is derived from adult doses by adjusting for body weight/size. As developmental and maturational changes are complex processes, such simplified methods may result in subtherapeutic effects or adverse events. Kidney function is impaired during the first 2 years of life as a result of normal growth and development. Reduced kidney function during childhood has an impact not only on renal clearance but also on absorption, distribution, metabolism and nonrenal clearance of drugs. 'Omics'-based technologies, such as proteomics and metabolomics, can be leveraged to uncover novel markers for kidney function during normal development, acute kidney injury, and chronic diseases. Pharmacometric modeling and simulation can be applied to simplify the design of pediatric investigations, characterize the effects of kidney function on drug exposure and response, and fine-tune dosing in pediatric patients, especially in those with impaired kidney function. One case study of amikacin dosing in neonates with reduced kidney function is presented. Collaborative efforts between clinicians and scientists in academia, industry, and regulatory agencies are required to evaluate new renal biomarkers, collect and share prospective pharmacokinetic, genetic and clinical data, build integrated pharmacometric models for key drugs, optimize and standardize dosing strategies, develop bedside decision tools, and enhance labels of drugs utilized in neonates, infants, and children.

Challenges and considerations for development of therapeutic proteins in pediatric patients

Target specificity and generally good tolerability of therapeutic proteins (TPs) present desirable treatment opportunities for pediatric patients. However, little is known on the ontogeny of processes related to the pharmacokinetics (PK) and disposition of TPs. The science, regulatory requirements and strategy of developing TPs for children are evolving. Our current review of TPs, (with focus on monoclonal antibodies and fusion proteins) that were approved for pediatric use indicates that dose-selection for pediatric pivotal studies is often based on adult PK information alone. This approach might not be sufficient if more complex PK properties than simple linear PK are present. Body weight-based dosing for pediatric patients directly scaled down from adult dosing can lead to under-exposure in young pediatric patients who are usually in the lowest body-weight range. Tiered-fixed dosing can be reasonably effective for TPs in achieving comparable exposure in children over a wide age range. The uniqueness of the pediatric population, the practical challenges in conducting clinical studies in this population, as well as regulations from health authorities warrant including pharmacometrics as an integral component of pediatric drug development. We propose a framework distinct from previous proposals, to guide clinical pharmacology strategy for pediatric drug development specifically for TPs.

Influence of covariate distribution on the predictive performance of pharmacokinetic models in paediatric research

AIMS

The accuracy of model-based predictions often reported in paediatric research has not been thoroughly characterized. The aim of this exercise is therefore to evaluate the role of covariate distributions when a pharmacokinetic model is used for simulation purposes.

METHODS

Plasma concentrations of a hypothetical drug were simulated in a paediatric population using a pharmacokinetic model in which body weight was correlated with clearance and volume of distribution. Two subgroups of children were then selected from the overall population according to a typical study design, in which pre-specified body weight ranges (10-15 kg and 30-40 kg) were used as inclusion criteria. The simulated data sets were then analyzed using non-linear mixed effects modelling. Model performance was assessed by comparing the accuracy of AUC predictions obtained for each subgroup, based on the model derived from the overall population and by extrapolation of the model parameters across subgroups.

RESULTS

Our findings show that systemic exposure as well as pharmacokinetic parameters cannot be accurately predicted from the pharmacokinetic model obtained from a population with a different covariate range from the one explored during model building. Predictions were accurate only when a model was used for prediction in a subgroup of the initial population.

CONCLUSIONS

In contrast to current practice, the use of pharmacokinetic modelling in children should be limited to interpolations within the range of values observed during model building. Furthermore, the covariate point estimate must be kept in the model even when predictions refer to a subset different from the original population.

Population pharmacokinetics of meropenem administered as a prolonged infusion in children with cystic fibrosis

OBJECTIVES

Meropenem is frequently used to treat pulmonary exacerbations in children with cystic fibrosis (CF) in the USA. Prolonged-infusion meropenem improves the time that free drug concentrations remain above the MIC (fT> MIC) in adults, but data in CF children are sparse. We describe the population pharmacokinetics, tolerability and treatment burden of prolonged-infusion meropenem in CF children.

METHODS

Thirty children aged 6-17 years with a pulmonary exacerbation received 40 mg/kg meropenem every 8 h; each dose was administered as a 3 h infusion. Pharmacokinetics were determined using population methods in Pmetrics. Monte Carlo simulation was employed to compare 0.5 with 3 h infusions to estimate the probability of pharmacodynamic target attainment (PTA) at 40% fT> MIC. NCT#01429259.

RESULTS

A two-compartment model fitted the data best with clearance and volume predicted by body weight. Clearance and volume of the central compartment were 0.41 ± 0.23 L/h/kg and 0.30 ± 0.17 L/kg, respectively. Half-life was 1.11 ± 0.38 h. At MICs of 1, 2 and 4 mg/L, PTAs for the 0.5 h infusion were 87.6%, 70.1% and 35.4%, respectively. The prolonged infusion increased PTAs to >99% for these MICs and achieved 82.8% at 8 mg/L. Of the 30 children, 18 (60%) completed treatment with prolonged infusion; 5 did so at home without any reported burden. Nine patients were changed to a 0.5 h infusion when discharged home.

CONCLUSIONS

In these CF children, meropenem clearance was greater compared with published values from non-CF children. Prolonged infusion provided an exposure benefit against pathogens with MICs ≥1 mg/L, was well tolerated and was feasible to administer in the hospital and home settings, the latter depending on perception and family schedule.

Population Pharmacokinetic Model of Sublingual Buprenorphine in Neonatal Abstinence Syndrome

OBJECTIVE

Neonatal abstinence syndrome (NAS)--a clinical entity of infants from in utero exposure to psychoactive xenobiotic and buprenorphine--has been successfully used to treat NAS. However, nothing is known about the pharmacokinetics (PK) of buprenorphine in neonates with NAS. To our knowledge, this is the first study to investigate the population pharmacokinetic of sublingual buprenorphine in neonates with NAS.

DESIGN

A retrospective population PK analysis of: (1) neonates with NAS treated with sublingual buprenorphine in randomized, double blinded clinical study and (2) data from healthy adults from a previously published pharmacokinetic study.

SETTING

Neonatal intensive care unit and general clinical research unit.

PATIENTS

Twenty-four neonates with NAS and five healthy adults.

INTERVENTIONS

All participants received sublingual buprenorphine per study protocol.

MEASUREMENTS AND MAIN RESULTS

A total of 303 PK data from 29 neonates and adults were used for model development. A population pharmacokinetic analysis was conducted using a first order conditional estimation with interaction in the NONMEM software program. A two-compartment linear PK model with first-order absorption process best described the pharmacokinetics of sublingual buprenorphine in neonates. The apparent clearance (CL) of buprenorphine was linearly related to body weight and matured with increasing age via two distinct saturated pathways. A typical neonate with NAS (body weight, 2.9 kg; postnatal age; 5.4 days) had a CL of 3.5 L/kg/hour and elimination half-life of 11 hours. Phenobarbital did not affect the clearance of buprenorphine compared to neonates of similar age and weight.

CONCLUSIONS

This is the first study to investigate the population PK of sublingual buprenorphine in neonatal NAS. To our knowledge, this is also the first report to describe the age-dependent changes of buprenorphine PK in this patient population. No buprenorphine dose adjustment is needed for neonates with NAS treated with buprenorphine and concurrent phenobarbital.