Tips and traps analyzing pediatric PK data

Evaluation and optimization of sample size of neonates and infants for pediatric clinical studies on cefiderocol using a model-based approach

When planning pediatric clinical trials, optimizing the sample size of neonates/infants is essential because it is difficult to enroll these subjects. In this simulation study, we evaluated the sample size of neonates/infants using a model-based optimal approach for identifying their pharmacokinetics for cefiderocol. We assessed the usefulness of data for estimation performance (accuracy and variance of parameter estimation) from adults and the impact of data from very young subjects, including preterm neonates. Stochastic simulation and estimation were utilized to assess the impact of sample size allocation for age categories in estimation performance for population pharmacokinetic parameters in pediatrics. The inclusion of adult pharmacokinetic information improved the estimation performance of population pharmacokinetic parameters as the coefficient of variation (CV) range of parameter estimation decreased from 4.9%-593.7% to 2.3%-17.3%. When sample size allocation was based on the age groups of gestational age and postnatal age, the data showed 15 neonates/infants would be necessary to appropriately estimate pediatric pharmacokinetic parameters (<20%CV). By using the postmenstrual age (PMA), which is theoretically considered to be associated with the maturation of organs, the number of neonates/infants required for appropriate parameter estimation could be reduced to seven (one and six with <32 and >32 weeks PMA, respectively) to nine (three and six with <37 and >37 weeks PMA, respectively) subjects. The model-based optimal design approach allowed efficient evaluation of the sample size of neonates/infants for estimation of pediatric pharmacokinetic parameters. This approach to assessment should be useful when designing pediatric clinical trials, especially those including young children.

Model-informed dose optimization of mycophenolic acid in pediatric kidney transplant patients

PURPOSE

We aimed to develop and evaluate a population PK model of mycophenolic acid (MPA) in pediatric kidney transplant patients to aid MPA dose optimization.

METHODS

Data were collected from pediatric kidney transplant recipients from a Dutch academic hospital (Radboudumc, the Netherlands). Pharmacokinetic model-building and model-validation analyses were performed using NONMEM. Subsequently, we externally evaluated the final model using data from another academic hospital. The final model was used to develop an optimized dosing regimen.

RESULTS

Thirty pediatric patients were included of whom 266 measured MPA plasma concentrations, including 20 full pharmacokinetic (PK) curves and 24 limited sampling curves, were available. A two-compartment model with a transition compartment for Erlang-type absorption best described the data. The final population PK parameter estimates were Ktr (1.48 h-1; 95% CI, 1.15-1.84), CL/F (16.0 L h-1; 95% CI, 10.3-20.4), Vc/F (24.9 L; 95% CI, 93.0-6.71E25), Vp/F (1590 L; 95% CI, 651-2994), and Q/F (36.2 L h-1; 95% CI, 9.63-74.7). The performance of the PK model in the external population was adequate. An optimized initial dose scheme based on bodyweight was developed. With the licensed initial dose, 35% of patients were predicted to achieve the target AUC, compared to 42% using the optimized scheme.

CONCLUSION

We have successfully developed a pharmacokinetic model for MPA in pediatric renal transplant patients. The optimized dosing regimen is expected to result in better target attainment early in treatment. It can be used in combination with model-informed follow-up dosing to further individualize the dose when PK samples become available.

A population pharmacokinetics model of balovaptan to support dose selection in adult and pediatric populations

Balovaptan is a brain-penetrating vasopressin receptor 1a antagonist previously investigated for the core symptoms of autism spectrum disorder (ASD). A population pharmacokinetic (PK) model of balovaptan was developed, initially to assist clinical dosing for adult and pediatric ASD studies and subsequently for new clinical indications including malignant cerebral edema (MCE) and post-traumatic stress disorder. The final model incorporates one-compartment disposition and describes time- and dose-dependent non-linear PK through empirical drug binding and a gut extraction component with turnover. An age effect on clearance observed in children was modeled by an asymptotic function that predicts adult-equivalent exposures at 40% of the adult dose for children aged 2-4 years, 70% for 5-9 years, and at the full adult dose for ≥ 10 years. The model was adapted for intravenous (IV) balovaptan dosing and combined with in vitro and ex vivo pharmacodynamic data to simulate brain receptor occupancy as a guide for dosing in a phase II trial of MCE prophylaxis after acute ischemic stroke. A sequence of three stepped-dose daily infusions of 50, 25 and 15 mg over 30 or 60 min was predicted to achieve a target occupancy of ≥ 80% in ≥ 95% of patients over a 3-day period. This model predicts both oral and IV balovaptan exposure across a wide age range and will be a valuable tool to analyze and predict its PK in new indications and target populations, including pediatric patients.

The MPRINT Hub Data, Model, Knowledge and Research Coordination Center: Bridging the gap in maternal-pediatric therapeutics research through data integration and pharmacometrics

Maternal and pediatric populations have historically been considered "therapeutic orphans" due to their limited inclusion in clinical trials. Physiologic changes during pregnancy and lactation and growth and maturation of children alter pharmacokinetics (PK) and pharmacodynamics (PD) of drugs. Precision therapy in these populations requires knowledge of these effects. Efforts to enhance maternal and pediatric participation in clinical studies have increased over the past few decades. However, studies supporting precision therapeutics in these populations are often small and, in isolation, may have limited impact. Integration of data from various studies, for example through physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) modeling or bioinformatics approaches, can augment the value of data from these studies, and help identify gaps in understanding. To catalyze research in maternal and pediatric precision therapeutics, the Obstetric and Pediatric Pharmacology and Therapeutics Branch of the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD) established the Maternal and Pediatric Precision in Therapeutics (MPRINT) Hub. Herein, we provide an overview of the status of maternal-pediatric therapeutics research and introduce the Indiana University-Ohio State University MPRINT Hub Data, Model, Knowledge and Research Coordination Center (DMKRCC), which aims to facilitate research in maternal and pediatric precision therapeutics through the integration and assessment of existing knowledge, supporting pharmacometrics and clinical trials design, development of new real-world evidence resources, educational initiatives, and building collaborations among public and private partners, including other NICHD-funded networks. By fostering use of existing data and resources, the DMKRCC will identify critical gaps in knowledge and support efforts to overcome these gaps to enhance maternal-pediatric precision therapeutics.

Update on Therapeutic Drug Monitoring of Beta-Lactam Antibiotics in Critically Ill Patients—A Narrative Review

Beta-lactam antibiotics remain one of the most preferred groups of antibiotics in critical care due to their excellent safety profiles and their activity against a wide spectrum of pathogens. The cornerstone of appropriate therapy with beta-lactams is to achieve an adequate plasmatic concentration of a given antibiotic, which is derived primarily from the minimum inhibitory concentration (MIC) of the specific pathogen. In a critically ill patient, the plasmatic levels of drugs could be affected by many significant changes in the patient’s physiology, such as hypoalbuminemia, endothelial dysfunction with the leakage of intravascular fluid into interstitial space and acute kidney injury. Predicting antibiotic concentration from models based on non-critically ill populations may be misleading. Therapeutic drug monitoring (TDM) has been shown to be effective in achieving adequate concentrations of many drugs, including beta-lactam antibiotics. Reliable methods, such as high-performance liquid chromatography, provide the accurate testing of a wide range of beta-lactam antibiotics. Long turnaround times remain the main drawback limiting their widespread use, although progress has been made recently in the implementation of different novel methods of antibiotic testing. However, whether the TDM approach can effectively improve clinically relevant patient outcomes must be proved in future clinical trials.

Pharmacokinetic parameters over time during sepsis and the association of target attainment and outcomes in critically ill children and young adults receiving ceftriaxone

INTRODUCTION

Early sepsis results in pharmacokinetic (PK) changes due to physiologic alterations. PK changes can lead to suboptimal drug target attainment, risking inadequate coverage from antibiotics like ceftriaxone. Little is known about how ceftriaxone PK and target attainment quantitatively change over time in patients with sepsis or the association between target attainment and outcomes in critically ill children and young adults.

METHODS

A retrospective analysis of a prospective study was conducted in a single-center pediatric intensive care unit. Septic patients given at least one ceftriaxone dose (commonly as 50 mg/kg every 12 h) and who had blood obtained in both the first 48 h of therapy (early) and afterwards (late) were included. Normalized clearance and central volume were estimated and compared in both sepsis phases. We evaluated target attainment, defined as concentrations above 1× or 4× the minimum inhibitory concentration (MIC) for 100% of dosing intervals, and investigated the association between target attainment and clinical outcomes.

RESULTS

Fifty-five septic patients (median age: 7.5 years) were included. Normalized clearance and central volume were similar in both phases (6.18 ± 1.48 L/h/70 kg early vs. 6.10 ± 1.61 L/h/70 kg late, p = 0.60; 26.6 [IQR 22.3, 31.3] L/70 kg early vs. 24.5 [IQR 22.0, 29.4] L/70 kg late, p = 0.18). Individual percent differences in normalized clearance and central volume between sepsis phases ranged from -39% to 276% and -51% to 212% (reference, late sepsis), respectively. Fewer patients attained the 1× MIC target in late sepsis (82% late vs. 96% early, p = 0.013), which was associated with transition to once daily dosing, typically done due to transfer from the pediatric intensive care unit (PICU) to a lower acuity unit. Failure to attain either target in late sepsis was associated with antibiotic broadening.

CONCLUSION

Ceftriaxone PK parameters were similar between early and late sepsis, but there were large individual differences. Fewer patients attained MIC targets in late sepsis and all who did not attain the less stringent target received once daily dosing during this period. The failure to attain targets in late sepsis was associated with antibiotic broadening and could be an area for antibiotic stewardship intervention.

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

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

Gentamicin Pharmacokinetics and Optimal Dosage in Infant Patients: A Case Report and Literature Review

Gentamicin is an aminoglycoside antibiotic that is mostly used for the pediatric population. While the pediatric population is classified into neonates, infants, children, and adolescents based on developmental or maturational changes, infants are often overlooked in research. Three infant cases receiving gentamicin are presented to illustrate the pharmacokinetics and optimum dosage of gentamicin. Three infant patients received gentamicin (5.6-7.5 mg/kg/day) for urinary tract infections (UTIs) or bacteremia caused by Enterobacter aerogenes. The trough (Cmin) and peak (Cpeak) concentrations of gentamicin were 0.2-1.8 and 8.9 mg/L, respectively. The Cmin of a patient receiving gentamicin at 9.0 mg/kg/day was 3.3 mg/L, and the patient showed a decrease in urinary volume. The other two patients fully recovered from the infection and did not experience any adverse events. Additionally, we reviewed three studies regarding infant patients receiving gentamicin. The studies used gentamicin therapy for Gram-negative pathogen infections and UTIs caused by Escherichia coli and Enterococcus faecalis. The Cmin and Cpeak of patients receiving gentamicin at 2.2-7.5 mg/kg/day were 0.58-2.15 mg/kg and 4.67-8.88 mg/L, respectively. All patients were cured without any adverse events. Gentamicin dosages below 7.5 mg/kg/day may be effective and safe for use in infant patients. However, the optimal dosing regimen of gentamicin in infant patients is controversial, and limited data are available.

Interactions of the protease inhibitor, ritonavir, with common anesthesia drugs

The protease inhibitor, ritonavir, is a strong inhibitor of CYP 3A. The drug is used for management of the human immunovirus and is currently part of an oral antiviral drug combination (nirmatrelvir-ritonavir) for the early treatment of SARS-2 COVID-19-positive patients aged 12 years and over who have recognized comorbidities. The CYP 3A enzyme system is responsible for clearance of numerous drugs used in anesthesia (e.g., alfentanil, fentanyl, methadone, rocuronium, bupivacaine, midazolam, ketamine). Ritonavir will have an impact on drug clearances that are dependent on ritonavir concentration, anesthesia drug intrinsic hepatic clearance, metabolic pathways, concentration-response relationship, and route of administration. Drugs with a steep concentration-response relationship (ketamine, midazolam, rocuronium) are mostly affected because small changes in concentration have major changes in effect response. An increase in midazolam concentration is observed after oral administration because CYP 3A in the gastrointestinal wall is inhibited, causing a large increase in relative bioavailability. Fentanyl infusion may be associated with a modest increase in plasma concentration and effect, but the large between subject variability of pharmacokinetic and pharmacodynamic concentration changes suggests it will have little impact on an individual patient, especially when used with adverse effect monitoring. It has been proposed that drugs that have no or only a small metabolic pathway involving the CYP 3A enzyme be used during anesthesia, for example, propofol, atracurium, remifentanil, and the volatile agents. That anesthesia approach denies children of drugs with considerable value. It is better that the inhibitory changes in clearance of these drugs are understood so that rational drug choices can be made to tailor drug use to the individual patient. Altered drug dose, anticipation of duration of effect, timing of administration, use of reversal agents and perioperative monitoring would better behoove children undergoing anesthesia.

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.

Population pharmacokinetic analysis of dexmedetomidine in children using real-world data from electronic health records and remnant specimens

AIMS

Our objectives were to perform a population pharmacokinetic analysis of dexmedetomidine in children using remnant specimens and electronic health records (EHRs) and explore the impact of patient's characteristics and pharmacogenetics on dexmedetomidine clearance.

METHODS

Dexmedetomidine dosing and patient data were gathered from EHRs and combined with opportunistically sampled remnant specimens. Population pharmacokinetic models were developed using nonlinear mixed-effects modelling. Stage 1 developed a model without genotype variables; Stage 2 added pharmacogenetic effects.

RESULTS

Our final study population included 354 post-cardiac surgery patients aged 0-22 years (median 16 mo). The data were best described with a 2-compartment model with allometric scaling for weight and Hill maturation function for age. Population parameter estimates and 95% confidence intervals were 27.3 L/h (24.0-31.1 L/h) for total clearance, 161 L (139-187 L) for central compartment volume of distribution, 26.0 L/h (22.5-30.0 L/h) for intercompartmental clearance and 7903 L (5617-11 119 L) for peripheral compartment volume of distribution. The estimate for postmenstrual age when 50% of adult clearance is achieved was 42.0 weeks (41.5-42.5 weeks) and the Hill coefficient estimate was 7.04 (6.99-7.08). Genotype was not statistically or clinically significant.

CONCLUSION

Our study demonstrates the use of real-world EHR data and remnant specimens to perform a population pharmacokinetic analysis and investigate covariate effects in a large paediatric population. Weight and age were important predictors of clearance. We did not find evidence for pharmacogenetic effects of UGT1A4 or UGT2B10 genotype or CYP2A6 risk score.

Effects of Posaconazole on Tacrolimus Population Pharmacokinetics and Initial Dose in Children With Crohn’s Disease Undergoing Hematopoietic Stem Cell Transplantation

The present study explored the effects of posaconazole on tacrolimus population pharmacokinetics (PPK) in children with Crohn’s disease (CD) undergoing hematopoietic stem cell transplantation (HSCT). Tacrolimus concentrations, physiological and biochemical factors, and concomitant medications from 51 CD children undergoing HSCT were used to establish a PPK model based on a nonlinear mixed-effect model. Steady-state concentrations of tacrolimus for children weighing less than 20 kg treated with different dose regimens were simulated by the Monte Carlo method. Weight and concomitant medications were included as covariates. At the same weight, the relative tacrolimus clearance was 1:0.43 in children without or with posaconazole. Compared to children not receiving posaconazole, the simulated tacrolimus steady-state concentrations at different doses for different body weights were all higher in children receiving posaconazole (p &lt; 0.01). Furthermore, in children not receiving posaconazole, the dosage regimen with the best probability of achieving the target concentration was 0.6 mg/kg/day for children weighing 5–8.2 kg and 0.5 mg/kg/day for children weighing 8.2–20 kg, while for children receiving posaconazole, the best probability of reaching the target concentration of tacrolimus was a dosage regimen of 0.5 mg/kg/day for children weighing 5–20 kg. In conclusion, the PPK for tacrolimus was determined in children with CD undergoing HSCT for the first time. Co-treatment with posaconazole significantly increased tacrolimus concentrations, and we recommend a specific initial dose regimen for tacrolimus.

Amikacin Therapy in Japanese Pediatric Patients: Narrative Review

Children show a very wide range of physical development processes. These changes impact pharmacokinetic (PK) variability in pediatric patients. Most PK studies have been conducted on the Caucasian population. Therefore, whether current evidence of how developmental change affects PK and exposure-response relationships applies to Japanese pediatric patients remains unclear. This narrative review focuses on amikacin therapy in Japanese pediatric patients and shows the relationship between amikacin concentrations and efficacy/toxicity. Ten relevant articles were identified. Of these, nine articles were published in the 1980s. All studies reported a maximum concentration (Cmax) and minimum concentration (Cmin) of amikacin. Overall, articles reporting PK/pharmacodynamic (PD) indices and minimum inhibitory concentration (MIC) of isolated bacteria in Japanese pediatric patients is lacking, whereas all patients recovered from an infection state and showed negative cultures. Five of the included studies reported the association between Cmin and toxicity. The Cmin in three of four patients who developed toxicity was above 10 mg/L. This narrative review shows that further PK study of amikacin in Japanese pediatric patients is necessary. In particular, the pursuit of knowledge of Cmax/MIC ratio is vital. On the other hand, this review demonstrates that the optimal Cmin for Japanese pediatric patients is below 10 mg/L as a candidate concentration. However, it is noted that the number of patients who developed toxicity is very small.

R-praziquantel integrated population pharmacokinetics in preschool- and school-aged African children infected with Schistosoma mansoni and S. haematobium and Lao adults infected with Opisthorchis viverrini

Racemic praziquantel (PZQ) is the standard treatment for schistosomiasis and liver fluke infections (opisthorchiasis and clonorchiasis). The development of an optimal pediatric formulation and dose selection would benefit from a population pharmacokinetic (popPK) model. A popPK model was developed for R-PZQ, the active enantiomer of PZQ, in 664 subjects, 493 African children (2-15 years) infected with Schistosoma mansoni and S. haematobium, and 171 Lao adults (15-78 years) infected with Opisthorchis viverrini. Racemate tablets were administered as single doses of 20, 40 and 60 mg/kg in children and 30, 40 and 50 mg/kg in 129 adults, and as 3 × 25 mg/kg apart in 42 adults. Samples collected by the dried-blood-spot technique were assayed by LC-MS/MS. A two-compartment disposition model, with allometric scaling and dual first-order and transit absorption, was developed using Phoenix™ software. Inversely parallel functions of age described the apparent oral bioavailability (BA) and clearance maturation in children and ageing in adults. BA decreased slightly in children with dose increase, and by 35% in adults with multiple dosing. Crushing tablets for preschool-aged children increased the first-order absorption rate by 64%. The mean transit absorption time was 70% higher in children. A popPK model for R-PZQ integrated African children over 2 years of age with schistosomiasis and Lao adults with opisthorchiasis, and should be useful to support dose optimization in children. In vitro hepatic and intestinal metabolism data would help refining and validating the model in younger children as well as in target ethnic pediatric and adult groups.

Pediatric pharmacokinetic considerations and implications for drug dosing

Optimizing the dosing of medicines for pediatric patients in routine clinical practice and determining the dose for clinical trials is still a challenging task. Children differ from adults in their response to drugs due to inherent differences in pharmacokinetics and/or pharmacodynamics, and responses may also vary among pediatric patients of different ages. However, the greatest disparities compared to adult pharmacokinetic profiles are observed in children below 2 years of age. The maturation of the liver and the kidneys, as well as the variation in body composition, are considered to be the main sources of pharmacokinetic variability. Hence, besides specific pharmacodynamic features, understanding age-related changes in drug absorption, distribution, and elimination is fundamental for optimizing drug efficacy and avoiding toxicity. This paper summarizes the pharmacokinetic changes throughout the childhood, along with the effect of developmental changes on drug dosage calculation. In clinical practice, age and body weight-based dosing regimens are usually used. In spite of dosing recommendations based on age and/or body weight, variabilities in pharmacokinetics and pharmacodynamic response remain, implying a need to monitor patients and optimize the dosing regimen according to physiological characteristics, disease characteristics and therapy.

Dosing strategies of imipenem in neonates based on pharmacometric modelling and simulation

OBJECTIVES

Imipenem is a broad-spectrum antibacterial agent used in critically ill neonates after failure of first-line treatments. Few studies have described imipenem disposition in this population. The objectives of our study were: (i) to characterize imipenem population pharmacokinetics (PK) in a cohort of neonates; and (ii) to conduct model-based simulations to evaluate the performance of six different dosing regimens aiming at optimizing PK target attainment.

METHODS

A total of 173 plasma samples from 82 neonates were collected over 15 years at the Lausanne University Hospital, Switzerland. The majority of study subjects were preterm neonates with a median gestational age (GA) of 27 weeks (range: 24-41), a postnatal age (PNA) of 21 days (2-153) and a body weight (BW) of 1.16 kg (0.5-4.1). PK data were analysed using non-linear mixed-effect modelling (NONMEM).

RESULTS

A one-compartment model best characterized imipenem disposition. Population PK parameters estimates of CL and volume of distribution were 0.21 L/h and 0.73 L, with an interpatient variability (CV%) of 20.1% on CL in a representative neonate (GA 27 weeks, PNA 21 days, BW 1.16 kg, serum creatinine, SCr 46.6 μmol/L). GA and PNA exhibited the greatest impact on PK parameters, followed by SCr. These covariates explained 36% and 15% of interindividual variability in CL, respectively.Simulated regimens using a dose of 20-25 mg/kg every 6-12 h according to postnatal age led to the highest PTA (T>MIC over 100% of time).

CONCLUSIONS

Dosing adjustment according to BW, GA and PNA optimizes imipenem exposure in neonates.

Population Pharmacokinetic Modeling of Total and Free Ceftriaxone in Critically Ill Children and Young Adults and Monte Carlo Simulations Support Twice Daily Dosing for Target Attainment

Critical illness, including sepsis, causes significant pathophysiologic changes that alter the pharmacokinetics (PK) of antibiotics. Ceftriaxone is one of the most prescribed antibiotics in patients admitted to the pediatric intensive care unit (PICU). We sought to develop population PK models of both total ceftriaxone and free ceftriaxone in children admitted to a single-center PICU using a scavenged opportunistic sampling approach. We tested if the presence of sepsis and phase of illness (before or after 48 hours of antibiotic treatment) altered ceftriaxone PK parameters. We performed Monte Carlo simulations to evaluate whether dosing regimens commonly used in PICUs in the United States (50 mg/kg every 12 hours vs. 24 hours) resulted in adequate antimicrobial coverage. We found that a two-compartment model best described both total and free ceftriaxone concentrations. For free concentrations, the population clearance value is 6.54 L/h/70 kg, central volume is 25.4 L/70 kg and the peripheral volume is 19.6 L/70kg. For both models, we found that allometric weight scaling, post-menstrual age, creatinine clearance and daily highest temperature had significant effects on clearance. Presence of sepsis or phase of illness did not have a significant effect on clearance or volume of distribution. Monte Carlo simulations demonstrated that to achieve free concentrations above 1 μg/mL for 100% of the dosing intervals, a dosing regimen of 50 mg/kg every 12 hours is recommended for most patients. A continuous infusion could be considered if the target is to maintain free concentrations four times above the minimum inhibitory concentrations (4 μg/mL).

Wide size dispersion and use of body composition and maturation improves the reliability of allometric exponent estimates

To evaluate study designs and the influence of dispersion of body size, body composition and maturation of clearance or reliable estimation of allometric exponents. Non-linear mixed effects modeling and parametric bootstrap were employed to assess how the study sample size, number of observations per subject, between subject variability (BSV) and dispersion of size distribution affected estimation bias and uncertainty of allometric exponents. The role of covariate model misspecification was investigated using a large data set ranging from neonates to adults. A decrease in study sample size, number of observations per subject, an increase in BSV and a decrease in dispersion of size distribution, increased the uncertainty of allometric exponent estimates. Studies conducted only in adults with drugs exhibiting normal (30%) BSV in clearance may need to include at least 1000 subjects to be able to distinguish between allometric exponents of 2/3 and 1. Nevertheless, studies including both children and adults can distinguish these exponents with only 100 subjects. A marked bias of 45% (95%CI 41-49%) in the estimate of the allometric exponent of clearance was obtained when maturation and body composition were ignored in infants. A wide dispersion of body size (e.g. infants, children and adults) is required to reliably estimate allometric exponents. Ignoring differences in body composition and maturation of clearance may bias the exponent for clearance. Therefore, pharmacometricians should avoid estimating allometric exponent parameters without suitable designs and covariate models. Instead, they are encouraged to rely on the well-developed theory and evidence that clearance and volume parameters in humans scale with theory-based exponents.

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

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

Model-Informed Pediatric Drug Development: Application of Pharmacometrics to Define the Right Dose for Children

One of the biggest challenges in pediatric drug development is defining a safe and effective dose in pediatric populations, which span across a wide age and development range from neonates to adolescents. Model-informed drug development approaches are particularly suited to address knowledge gaps including data leveraging to increase the success of pediatric studies. Considering the often limited number of patients available for study and logistic difficulties to collect the necessary data in pediatric populations, the application of pharmacometrics and modeling and simulation techniques can improve clinical trial efficiency, increase the probability of regulatory success, and optimize therapeutic individualization in support of dedicated trials. This review describes the state of pediatric model-informed drug development to define the right dose for children and provides suggestions for future development.

Test-dose pharmacokinetics guided melphalan dose adjustment in reduced intensity conditioning allogeneic transplant for non-malignant disorders

AIMS

We studied melphalan pharmacokinetics (PK) and feasibility of melphalan full-dose adjustment based on test-dose PK in children and young adults with non-malignant disorders (NMD) undergoing allogeneic hematopoietic cell transplantation (HCT) using reduced intensity conditioning (RIC) containing alemtuzumab, fludarabine and melphalan.

METHODS

Patients received test-dose melphalan (10% of planned full-dose) prior to conditioning. Blood samples for PK were obtained around test and full-dose melphalan (140 mg/m² or 4.7 mg/kg in patients <10 kg). Melphalan concentration was measured by liquid chromatography electrospray ionization tandem mass-spectrometry assay and data were analysed using a population-PK model and Bayesian estimation. Test and full-dose melphalan clearance estimates were evaluated by pairwise Wilcoxon test and Bland-Altman plot.

RESULTS

Twenty-four patients undergoing 25 transplants were included in the final analysis. Patients received standard full-dose melphalan in 17 transplants, with median area under the concentration-time curve (AUC) of 5.5 mg*h/L (range, 3.0-9.5 mg*h/L). Patients received test-dose melphalan in 23 transplants with a test-dose PK predicted full-dose AUC range of 2.9-16.8 mg*h/L. In seven transplants where patients had baseline organ impairment, test-dose PK predicted higher exposure for standard full-dose (median AUC 13.8 mg*h/L). Melphalan full-dose was adjusted in these patients, with achievement of desired target AUC (3.6-5.4 mg*h/L) and no excess toxicity. Mean ratio of test-dose clearance to full-dose clearance was 1.03. Twenty of 22 patients (91%) were within the 95% confidence intervals of the clearance ratio.

CONCLUSION

Melphalan test-dose PK reliably predicts full-dose PK and allows for accurate adjustment of full-dose melphalan in RIC-HCT for NMD. This approach can avoid excess toxicity from increased systemic exposure, especially in patients with organ impairment.

Why do We Use the Concepts of Adult Anesthesia Pharmacology in Developing Brains? Will It Have an Impact on Outcomes? Challenges in Neuromonitoring and Pharmacology in Pediatric Anesthesia

BACKGROUND

Pediatric sedation and anesthesia techniques have plenty of difficulties and challenges. Data on the pharmacologic, electroencephalographic, and neurologic response to anesthesia at different brain development times are only partially known. New data in neuroscience, pharmacology, and intraoperative neuromonitoring will impact changing concepts and clinical practice. In this article, we develop a conversation to guide the debate and search for a view more attuned to the updated knowledge in neurodevelopment, electroencephalography, and clinical pharmacology for the anesthesiologic practice in the pediatric population.

Initial dose recommendation for sirolimus in paediatric kaposiform haemangioendothelioma patients based on population pharmacokinetics and pharmacogenomics

Objective

Sirolimus has been used to treat paediatric kaposiform haemangioendothelioma patients. However, there is considerable pharmacokinetic variability among individuals, and it is difficult to develop an initial dosing regimen. The goal of the present study is to recommend an initial sirolimus dose in paediatric kaposiform haemangioendothelioma patients based on population pharmacokinetics and pharmacogenomics.

Methods

This was a retrospective clinical study. A population pharmacokinetics model was established and population characteristics, laboratory test results, drug combinations, and pharmacogenomics were considered as potential covariates. The Monte Carlo method was used to simulate the optimal initial dosage.

Results

The final covariates that affect sirolimus clearance include weight and the CYP3A5 genotype. The initial dosage of sirolimus for individuals with CYP3A5*3/*3 was 0.20 mg/kg split into two doses for 5 to 60 kg body weight. For individuals with CYP3A5*1, the initial dose was 0.23 mg/kg split into two doses for 5 to 30 kg body weight and 0.20 mg/kg split into two doses for 30 to 60 kg body weight.

Conclusion

The recommendation for the initial sirolimus dose in paediatric kaposiform haemangioendothelioma patients was based on population pharmacokinetics and pharmacogenomics. This study may provide practical value for sirolimus clinical use in paediatric kaposiform haemangioendothelioma patients.

Contribution of Population Pharmacokinetics of Glycopeptides and Antifungals to Dosage Adaptation in Paediatric Onco-hematological Malignancies: A Review

The response to medications in children differs not only in comparison to adults but also between children of the different age groups and according to the disease. This is true for anti-infectives that are widely prescribed in children with malignancy. In the absence of pharmacokinetic/pharmacodynamic paediatric studies, dosage is frequently based on protocols adapted to adults. After a short presentation of the drugs, we reviewed the population pharmacokinetic studies available for glycopeptides (vancomycin and teicoplanin, n = 5) and antifungals (voriconazole, posaconazole, and amphotericin B, n = 9) currently administered in children with onco-hematological malignancies. For each of them, we reported the main study characteristics including identified covariates affecting pharmacokinetics and proposed paediatric dosage recommendations. This review highlighted the very limited amount of data available, the lack of consensus regarding PK/PD targets used for dosing optimization and regarding dosage recommendations when available. Additional PK studies are urgently needed in this specific patient population. In addition to pharmacokinetics, efficacy may be altered in immunocompromised patients and prospective clinical evaluation of new dosage regimen should be provided as they are missing in most cases.

Clinical pharmacology and dosing regimen optimization of neonatal opioid withdrawal syndrome treatments

In this paper, we review the management of neonatal opioid withdrawal syndrome (NOWS) and clinical pharmacology of primary treatment agents in NOWS, including morphine, methadone, buprenorphine, clonidine, and phenobarbital. Pharmacologic treatment strategies in NOWS have been mostly empirical, and heterogeneity in dosing regimens adds to the difficulty of extrapolating study results to broader patient populations. As population pharmacokinetics (PKs) of pharmacologic agents in NOWS become more well‐defined and knowledge of patient‐specific factors affecting treatment outcomes continue to accumulate, PK/pharmacodynamic modeling and simulation will be powerful tools to aid the design of optimal dosing regimens at the patient level. Although there is an increasing number of clinical trials on the comparative efficacy of treatment agents in NOWS, here, we also draw attention to the importance of optimizing the dosing regimen, which can be arguably equally important at identifying the optimal treatment agent.

External evaluation of the predictive performance of seven population pharmacokinetic models for phenobarbital in neonates

AIM

Several studies have reported population pharmacokinetic models for phenobarbital (PB), but the predictive performance of these models has not been well documented. This study aims to do external evaluation of the predictive performance in published pharmacokinetic models.

METHODS

Therapeutic drug monitoring data collected in neonates and young infants treated with PB for seizure control was used for external evaluation. A literature review was conducted through PubMed to identify population pharmacokinetic models. Prediction- and simulation-based diagnostics, and Bayesian forecasting were performed for external evaluation. The incorporation of allometric scaling for body size and maturation factors into the published models was also tested for prediction improvement.

RESULTS

A total of 79 serum concentrations from 28 subjects were included in the external dataset. Seven population pharmacokinetic studies of PB were identified as relevant in the literature search and included for our evaluation. The model by Voller et al showed the best performance concerning prediction-based evaluation. In simulation-based analyses, the normalized prediction distribution error of two models (those of Shellhaas et al and Marsot et al) obeyed a normal distribution. Bayesian forecasting with more than one observation improved predictive capability. Incorporation of both allometric size scaling and maturation function generally enhanced the predictive performance, with improvement as observed in the model of Vucicevic et al. CONCLUSIONS: The predictive performance of published pharmacokinetic models of PB was diverse. Bayesian forecasting and incorporation of both size and maturation factors could improve the predictability of the models for neonates.

Effect of Body Weight and Age on the Pharmacokinetics of Dihydroartemisinin: FDA Basis for Dose Determination of Artesunate for Injection in Pediatric Patients with Severe Malaria

For treatment of severe malaria, the WHO recommends 3 mg/kg intravenous artesunate in pediatric patients weighing less than 20 kg. Here we describe FDA's rationale for selecting 2.4 mg/kg in pediatric patients weighing less than 20 kg based on literature review and independent analyses.

Population Pharmacokinetic Analysis and Dose Regimen Optimization in Japanese Infants with an Extremely Low Birth Weight

Vancomycin is a synthetic antibiotic effective against Gram-positive pathogens. Although the clinical applicability of vancomycin for infants has been increasing, the pharmacokinetic data for vancomycin in extremely low-birth-weight infants are limited. The aim of this study was to construct a population pharmacokinetics model for vancomycin in extremely-low-birth-weight infants and establish an optimal dosage regimen. We enrolled children aged less than 1 year with a birth weight of less than 1,000 g and body weight at vancomycin prescription of less than 1,500 g. Pharmacokinetic data from 19 patients were analyzed, and a population pharmacokinetics model was developed using nonlinear mixed-effects modeling software. Goodness-of-fit plots, a nonparametric bootstrap analysis, and a prediction-corrected visual predictive check were employed to evaluate the final model. The dosage regimen was optimized based on the final model. The pharmacokinetic data fit a one-compartment model with first-order elimination, and body weight and estimated serum creatinine level were used as significant covariates. In a simulation using the final model, the optimal dosage regimen, especially when the serum creatinine level (>0.6 mg/dl) was high, was 5.0 to 7.5 mg/kg of body weight twice a day every 12 h; this was required to reduce the dosage compared with that in previous studies. The recommended doses based on the current target time course concentration curves may not be appropriate for extremely-low-birth-weight infants.

Opioid Treatment for Neonatal Opioid Withdrawal Syndrome: Current Challenges and Future Approaches

Chronic intrauterine exposure to psychoactive drugs often results in neonatal opioid withdrawal syndrome (NOWS). When nonpharmacologic measures are insufficient in controlling NOWS, morphine, methadone, and buprenorphine are first-line medications commonly used to treat infants with NOWS because of in utero exposure to opioids. Research suggests that buprenorphine may be the leading drug therapy used to treat NOWS when compared with morphine and methadone. Currently, there are no consensus or standardized treatment guidelines for medications prescribed for NOWS. Opioids used to treat NOWS exhibit large interpatient variability in pharmacokinetics (PK) and pharmacodynamic (PD) response in neonates. Organ systems undergo rapid maturation after birth that may alter drug disposition and exposure for any given dose during development. Data regarding the PK and PD of opioids in neonates are sparse. Pharmacometric methods such as physiologically based pharmacokinetic and population pharmacokinetic modeling can be used to explore factors predictive of some of the variability associated with the PK/PD of opioids in newborns. This review discusses the utility of pharmacometric techniques for enhancing precision dosing in infants requiring opioid treatment for NOWS. Applying these approaches may contribute to optimizing the outcome by reducing cumulative drug exposure, mitigating adverse drug effects, and reducing the burden of NOWS in neonates.

A Pilot Randomized, Controlled, Double-Blind Trial of Bumetanide to Treat Neonatal Seizures

OBJECTIVE

In the absence of controlled trials, treatment of neonatal seizures has changed minimally despite poor drug efficacy. We tested bumetanide added to phenobarbital to treat neonatal seizures in the first trial to include a standard-therapy control group.

METHODS

A randomized, double-blind, dose-escalation design was employed. Neonates with postmenstrual age 33 to 44 weeks at risk of or with seizures were eligible. Subjects with electroencephalography (EEG)-confirmed seizures after ≥20 and <40mg/kg phenobarbital were randomized to receive additional phenobarbital with either placebo (control) or 0.1, 0.2, or 0.3mg/kg bumetanide (treatment). Continuous EEG monitoring data from ≥2 hours before to ≥48 hours after study drug administration (SDA) were analyzed for seizures.

RESULTS

Subjects were randomized to treatment (n = 27) and control (n = 16) groups. Pharmacokinetics were highly variable among subjects and altered by hypothermia. The only statistically significant adverse event was diuresis in treated subjects (48% vs 13%, p = 0.02). One treated (4%) and 3 control subjects died (19%, p = 0.14). Among survivors, 2 of 26 treated subjects (8%) and 0 of 13 control subjects had hearing impairment, as did 1 nonrandomized subject. Total seizure burden varied widely, with much higher seizure burden in treatment versus control groups (median = 3.1 vs 1.2 min/h, p = 0.006). There was significantly greater reduction in seizure burden 0 to 4 hours and 2 to 4 hours post-SDA (both p < 0.01) compared with 2-hour baseline in treatment versus control groups with adjustment for seizure burden.

INTERPRETATION

Although definitive proof of efficacy awaits an appropriately powered phase 3 trial, this randomized, controlled, multicenter trial demonstrated an additional reduction in seizure burden attributable to bumetanide over phenobarbital without increased serious adverse effects. Future trials of bumetanide and other drugs should include a control group and balance seizure severity. ANN NEUROL 2021;89:327-340.

A POETIC Phase II study of continuous oral everolimus in recurrent, radiographically progressive pediatric low-grade glioma

BACKGROUND

To evaluate efficacy, pharmacokinetics (PK) and pharmacodynamics of single-agent everolimus in pediatric patients with radiographically progressive low-grade glioma (LGG).

METHODS

Everolimus was administered at 5 mg/m² once daily as a tablet or liquid for a planned 48-week duration or until unacceptable toxicity or disease progression. Patients with neurofibromatosis type 1 were excluded. PK and pharmacodynamic endpoints were assessed in consenting patients.

RESULTS

Twenty-three eligible patients (median age 9.2 years) were enrolled. All patients received prior chemotherapy (median number of prior regimens two) and/or radiotherapy (two patients). By week 48, two patients had a partial response, 10 stable disease, and 11 clinical or radiographic progression; two discontinued study prior to 1 year (toxicity: 1, physician determination: 1). With a median follow up of 1.8 years (range 0.2-6.7 years), the 2-, 3-, and 5-year progression-free survivals (PFS) were 39 ± 11%, 26 ± 11%, and 26 ± 11%, respectively; two patients died of disease. The 2-, 3-, and 5-year overall survival (OS) were all 93 ± 6%. Grade 1 and 2 toxicities predominated; two definitively related grade 3 toxicities (mucositis and neutropenia) occurred. Grade 4 elevation of liver enzymes was possibly related in one patient. Predose blood levels showed substantial variability between patients with 45.5% below and 18.2% above the target range of 5-15 ng/mL. Pharmacodynamic analysis demonstrated significant inhibition in phospho-S6, 4E-BP1, and modulation of c-Myc expression.

CONCLUSION

Daily oral everolimus provides a well-tolerated, alternative treatment for multiple recurrent, radiographically progressive pediatric LGG. Based on these results, everolimus is being investigated further for this patient population.

Model-Informed Precision Dosing of Antibiotics in Pediatric Patients: A Narrative Review

Optimal pharmacotherapy in pediatric patients with suspected infections requires understanding and integration of relevant data on the antibiotic, bacterial pathogen, and patient characteristics. Because of age-related physiological maturation and non-maturational covariates (e.g., disease state, inflammation, organ failure, co-morbidity, co-medication and extracorporeal systems), antibiotic pharmacokinetics is highly variable in pediatric patients and difficult to predict without using population pharmacokinetics models. The intra- and inter-individual variability can result in under- or overexposure in a significant proportion of patients. Therapeutic drug monitoring typically covers assessment of pharmacokinetics and pharmacodynamics, and concurrent dose adaptation after initial standard dosing and drug concentration analysis. Model-informed precision dosing (MIPD) captures drug, disease, and patient characteristics in modeling approaches and can be used to perform Bayesian forecasting and dose optimization. Incorporating MIPD in the electronic patient record system brings pharmacometrics to the bedside of the patient, with the aim of a consisted and optimal drug exposure. In this narrative review, we evaluated studies assessing optimization of antibiotic pharmacotherapy using MIPD in pediatric populations. Four eligible studies involving amikacin and vancomycin were identified from 418 records. Key articles, independent of year of publication, were also selected to highlight important attributes of MIPD. Although very little research has been conducted until this moment, the available data on vancomycin indicate that MIPD is superior compared to conventional dosing strategies with respect to target attainment. The utility of MIPD in pediatrics needs to be further confirmed in frequently used antibiotic classes, particularly aminoglycosides and beta-lactams.

Hydromorphone population pharmacokinetics in pediatric surgical patients

BACKGROUND

Hydromorphone is an opioid agonist used for pediatric analgesia. Due to lack of data, pediatric dosing (based on adult pharmacokinetic models) is not optimal.

AIM

This study characterizes hydromorphone population pharmacokinetics in pediatric surgical patients.

METHODS

In this prospective observational study, 34 children (4-18 years, bodyweight 23-89.6 kg) received multiple intravenous hydromorphone boluses followed by postoperative hydromorphone patient-controlled analgesia. Arterial blood samples were collected before and at 3, 10, 30, and 90 (and few samples at 1350) minutes after the first dose. Hydromorphone concentrations were measured by validated LC-MS/MS assay. Nonlinear mixed-effects modeling was used for pharmacokinetic model development. The final population pharmacokinetic model was evaluated by visual predictive check and bootstrap analysis. Monte Carlo simulations based on the final pharmacokinetic model determined optimal patient-controlled analgesia parameters to achieve a target of 20 ng/mL (as the median effective analgesic concentration), using minimum effective analgesic concentration of 4 ng/mL as a proxy for patient-controlled analgesia dose demand, and not exceeding the defined safe upper threshold of 40 ng/mL.

RESULTS

Hydromorphone pharmacokinetic profiles were adequately described by a two-compartmental model with first-order elimination. Bodyweight was found to be a significant covariate for hydromorphone clearance. Allometrically scaledpharmacokinetic parameter estimates (per 70 kg), systemic clearance (0.748 L/min), volume of distribution (33 L), peripheral clearance (1.57 L/min), and peripheral volume of distribution (146 L) were similar to reported adult parameter estimates. Sex, race, age, and type of surgery were not identified as significant covariates. To identify optimal patient-controlled analgesia dosing parameters, we simulated several initial loading doses, demand doses, and lockout intervals. Our simulations support an initial patient-controlled analgesia loading dose of 15 µg/kg followed by a demand dose of 6 µg/kg with lockout intervals of 20 minutes.

CONCLUSIONS

After intravenous hydromorphone, plasma pharmacokinetic profiles in children undergoing different surgeries were well described by a two-compartment population allometric pharmacokinetic model using bodyweight as the size descriptor. Model informed simulations identified patient-controlled analgesia parameters to inform initial settings, with adjustments as needed based on observed individual effects.

Initial dosage optimization of ciclosporin in pediatric Chinese patients who underwent bone marrow transplants based on population pharmacokinetics

Bone marrow transplants (BMT) are an established therapeutic strategy for patients with severe aplastic anemia, acute lymphoblastic leukemia, acute myeloid leukemia or chronic myeloid leukemia. However, the successful application of BMT is limited by graft-vs.-host disease (GVHD). Ciclosporin has been widely used for treating GVHD in pediatric patients who underwent BMT. The present study aimed to optimize the dosage of ciclosporin for safety and effectiveness based on population pharmacokinetics. A non-linear mixed-effects model was used to analyze the clinical data of pediatric patients who underwent BMT between September 2016 and September 2019 at the Children's Hospital of Fudan University. Monte Carlo simulations were used to identify the optimal dose of ciclosporin. The final population pharmacokinetic model indicated that body weight and days post-transplant influenced the clearance of ciclosporin in pediatric patients who underwent BMT. The present study indicated that the optimal initial dose of ciclosporin for pediatric patients weighing 5-30 kg who underwent BMT was 6 mg/kg/day split into 2 doses.

Pharmacokinetics and population pharmacokinetics in pediatric oncology

Pharmacokinetic research has become increasingly important in pediatric oncology as it can have direct clinical implications and is a crucial component in individualized medicine. Population pharmacokinetics has become a popular method especially in children, due to the potential for sparse sampling, flexible sampling times, computing of heterogeneous data, and identification of variability sources. However, population pharmacokinetic reports can be complex and difficult to interpret. The aim of this article is to provide a basic explanation of population pharmacokinetics, using clinical examples from the field of pediatric oncology, to facilitate the translation of pharmacokinetic research into the daily clinic.

External Validation of a Pharmacokinetic Model of Propofol for Target-Controlled Infusion in Children under Two Years Old

BACKGROUND

Previously, a linked pharmacokinetic-pharmacodynamic model (the Kim model) of propofol with concurrent infusion of remifentanil was developed for children aged 2-12 years. There are few options for pharmacokinetic-pharmacodynamic model of propofol for children under two years old. We performed an external validation of the Kim model for children under two years old to evaluate whether the model is applicable to this age group.

METHODS

Twenty-four children were enrolled. After routine anesthetic induction, a continuous infusion of 2% propofol and remifentanil was commenced using the Kim model. The target effect-site concentration of propofol was set as 2, 3, 4, and 5 μg/mL, followed by arterial blood sampling after 10 min of each equilibrium. Population estimates of four parameters-pooled bias, inaccuracy, divergence, and wobble-were used to evaluate the performance of the Kim model.

RESULTS

A total of 95 plasma concentrations were used for evaluation of the Kim model. The population estimate (95% confidence interval) of bias was -0.96% (-8.45%, 6.54%) and that of inaccuracy was 21.0% (15.0%-27.0%) for the plasma concentration of propofol.

CONCLUSION

The pooled bias and inaccuracy of the pharmacokinetic predictions are clinically acceptable. Therefore, our external validation of the Kim model indicated that the model can be applicable to target-controlled infusion of propofol in children younger than 2 years, with the recommended use of actual bispectral index monitoring in clinical settings that remifentanil is present.

TRIAL REGISTRATION

Clinical Research Information Service Identifier: KCT0001752.

Pharmacokinetic and pharmacodynamic considerations of general anesthesia in pediatric subjects

Introduction: The target concentration strategy uses PKPD information for dose determination. Models have also quantified exposure-response relationships, improved understanding of developmental pharmacokinetics, rationalized dose prescription, provided insight into the importance of covariate information, explained drug interactions and driven decision-making and learning during drug development.Areas covered: The prime PKPD consideration is parameter estimation and quantification of variability. The main sources of variability in children are age (maturation) and weight (size). Model use is mostly confined to pharmacokinetics, partly because anesthesia effect models in the young are imprecise. Exploration of PK and PD covariates and their variability hold potential to better individualize treatment.Expert opinion: The ability to model drugs using computer-based technology is hindered because covariate data required to individualize treatment using these programs remain lacking. Target concentration intervention strategies remain incomplete because covariate information that might better predict individualization of dose is absent. Pharmacogenomics appear a valuable area for investigation for pharmacodynamics and pharmacodynamics. Effect measures in the very young are imprecise. Assessment of the analgesic component of anesthesia is crude. While neuromuscular monitoring is satisfactory, depth of anaesthesia EEG interpretation is inadequate. Closed loop anesthesia is possible with better understanding of EEG changes.

Implementation of a Physiologically Based Pharmacokinetic Modeling Approach to Guide Optimal Dosing Regimens for Imatinib and Potential Drug Interactions in Paediatrics

Long-term use of imatinib is effective and well-tolerated in children with chronic myeloid leukaemia (CML) yet defining an optimal dosing regimen for imatinib in younger patients is a challenge. The potential interactions between imatinib and coadministered drugs in this "special" population also remains largely unexplored. This study implements a physiologically based pharmacokinetic (PBPK) modeling approach to investigate optimal dosing regimens and potential drug interactions with imatinib in the paediatric population. A PBPK model for imatinib was developed in the Simcyp Simulator (version 17) utilizing in silico, in vitro drug metabolism, and in vivo pharmacokinetic data and verified using an independent set of published clinical pharmacokinetic data. The model was then extrapolated to children and adolescents (aged 2-18 years) by incorporating developmental changes in organ size and maturation of drug-metabolising enzymes and plasma protein responsible for imatinib disposition. The PBPK model described imatinib pharmacokinetics in adult and paediatric populations and predicted drug interaction with carbamazepine, a cytochrome P450 (CYP)3A4 and 2C8 inducer, with a good accuracy (evaluated by visual inspections of the simulation results and predicted pharmacokinetic parameters that were within 1.25-fold of the clinically observed values). The PBPK simulation suggests that the optimal dosing regimen range for imatinib is 230-340 mg/m2/d in paediatrics, which is supported by the recommended initial dose for treatment of childhood CML. The simulations also highlighted that children and adults being treated with imatinib have similar vulnerability to CYP modulations. A PBPK model for imatinib was successfully developed with an excellent performance in predicting imatinib pharmacokinetics across age groups. This PBPK model is beneficial to guide optimal dosing regimens for imatinib and predict drug interactions with CYP modulators in the paediatric population.

Multicenter-Based Population Pharmacokinetic Analysis of Ciclosporin in Hematopoietic Stem Cell Transplantation Patients

PURPOSE

To explore the contribution of physiological characteristics to variability in ciclosporin pharmacokinetics in hematopoietic stem cell transplantation patients.

METHODS

Clinical data from 563 patients were collected from centers in three regions. Ciclosporin concentrations were measured using immunoassays. The patients' demographics, hematological and biological indicators, coadministered drugs, region, and disease diagnosis were recorded from medical records. Data analysis was performed using NONMEM based on a one-compartment model to describe the pharmacokinetics of ciclosporin. The reliability and stability of the final model were evaluated using bootstrap resampling, goodness-of-fit plots, and prediction-corrected visual predictive checks.

RESULTS

The population estimate of the clearance (CL) was 30.4 L/h, the volume of distribution (V) was 874.0 L and the bioavailability (F) was 81.1%. The between-subject variability in these parameters was 26.3, 68.0, and 110.8%, respectively. Coadministration of fluconazole, itraconazole, or voriconazole decreased CL by 17.6%, 28.4%, and 29.2%, respectively. Females' CL increased by approximately 12.0%. In addition, CL and V decreased with hematocrit, total protein, and uric acid increase, and CL also decreased with age and aspartate aminotransferase increase. However, CL increased with creatinine clearance increase.

CONCLUSIONS

A multicenter-based population pharmacokinetic model of ciclosporin was established. The pharmacokinetics of ciclosporin exhibited discrepancies among different regions.

Antithrombin Population Pharmacokinetics in Pediatric Ventricular Assist Device Patients

OBJECTIVES

Describe the pharmacokinetics of antithrombin in pediatric patients undergoing ventricular assist device therapy and provide dosing recommendations for antithrombin in this population.

DESIGN

A retrospective population pharmacokinetic study was designed.

SETTING

Large tertiary care children's hospital Subject inclusion criteria consisted of less than 19 years old.

PATIENTS

Subjects less than 19 years old undergoing therapy with a HeartWare ventricular assist device (HeartWare, Framingham, MA) or Berlin EXCOR ventricular assist device (Berlin GmbH, Berlin, Germany), who received a dose of antithrombin with a postdose antithrombin activity level from January 1, 2011, to June 30, 2017.

INTERVENTIONS

Population pharmacokinetic analysis and simulation using NONMEM v.7.4 (Icon, PLC, Dublin, Ireland).

MEASUREMENTS AND MAIN RESULTS

A total of 41 patients met study criteria (median age, 5.8 years [interquartile range, 1.6-9.9 yr]), and 53.7% underwent therapy with the pulsatile Berlin EXCOR pediatric ventricular assist device (Berlin Heart GmbH, Berlin, Germany). All patients received unfractionated heparin continuous infusion at a mean ± SD dose of 29 ± 14 U/kg/hr. A total of 181 antithrombin doses (44.1 ± 24.6 U/kg/dose) were included, and baseline antithrombin activity levels were 77 ± 12 U/dL. Antithrombin activity levels were drawn a median 19.9 hours (interquartile range, 8.8-41.6 hr) after antithrombin dose. A one-compartment proportional error model best fit the data, with allometric scaling of fat-free mass providing a better model fit than actual body weight. Unfractionated heparin and baseline antithrombin were identified as significant covariates. A 50 U/kg dose of antithrombin had a simulated half-life 13.2 ± 6.6 hours.

CONCLUSIONS

Antithrombin should be dosed on fat-free mass in pediatric ventricular assist device patients. Unfractionated heparin dose and baseline antithrombin activity level should be considered when dosing antithrombin in pediatric ventricular assist device patients.

Population pharmacokinetics of tacrolimus in pediatric patients with systemic-onset juvenile idiopathic arthritis: Initial dosage recommendations

Pediatric patients with systemic-onset juvenile idiopathic arthritis (SOJIA) may be treated with tacrolimus. However, the therapeutic range for tacrolimus is narrow with considerable inter- and intra-individual variability, making it difficult to formulate an ideal dosage regimen for personalized treatment. The purpose of the present study was to set up a population pharmacokinetics (PPK) model of tacrolimus treatment for SOJIA to determine the optimal initial dosage. Patients with SOJIA were analyzed using non-linear mixed-effects modeling. Different regimens were analyzed using Monte Carlo simulation with concentration profiles. A first-order absorption and elimination one-compartment model was selected as the most appropriate model for SOJIA. Based on initial dosage recommendations, the regimen of 0.5 mg every 24 h (q24h) appeared to be most suitable for subjects with a body weight of 5 kg, while the 0.5 mg q12h regimen was most suitable for subjects with a body weight of 15–25 kg, the 1/0.5 mg q24h regimen was appropriate for the 26–35 kg group and the 1 mg q12h regimen was suitable for the subjects with a body weight of 36–50 kg. To the best of our knowledge, the present study established the first PPK model of tacrolimus treatment that may be used for the selection of the initial dose based on body weight of pediatric patients with SOJIA.

Neurological excitation in a 6-week-old male infant after morphine overdose

While respiratory depression is a known complication of morphine overdose, the neuro-excitatory side effect of the morphine metabolite morphine-3-glucuronide is less widely known. Here, we report the case of an infant with neurological excitation after morphine overdose. The neuro-excitation in this infant was probably induced by an elevated morphine-3-glucuronide concentration.

A Population Pharmacokinetic Model of Intravenous Dexmedetomidine for Mechanically Ventilated Children after Neurosurgery

Dexmedetomidine is a selective alpha-2 adrenergic agonist with concurrent sedative and analgesic effects, and it is being increasingly used in pediatric anesthesia and intensive care. This study aimed to investigate the pharmacokinetics of intravenous dexmedetomidine in mechanically ventilated children in the intensive care unit (ICU) after neurosurgery. Pediatric patients aged 2–12 years, who were mechanically ventilated in ICU after neurosurgery, were allocated into a low-dose (n = 15) or high-dose (n = 14) group. The low-dose group received dexmedetomidine at a loading dose of 0.25 µg/kg for 10 min, followed by a maintenance dose of 0.25 µg/kg/h for 50 min, whereas the high-dose group received dexmedetomidine at a loading dose of 0.5 µg/kg for 10 min, followed by a maintenance dose of 0.5 µg/kg/h for 50 min. Serial blood samples were collected for a pharmacokinetic analysis up to 480 min after the end of the infusion. The sedative effect of dexmedetomidine was assessed using the Bispectral Index and University of Michigan Sedation Scale. Adverse reactions, electrocardiography findings, and vital signs were monitored for a safety assessment. A population pharmacokinetic analysis was performed using non-linear mixed effects modeling. Dexmedetomidine induced a moderate-to-deep degree of sedation during infusion in both groups. The pharmacokinetics of dexmedetomidine were best described by a two-compartment disposition model with first-order elimination kinetics. The parameters were standardized for a body weight of 70 kg using an allometric power model. The population estimates (95% confidence interval) per 70 kg body weight were as follows: clearance of 81.0 (72.9–90.9) L/h, central volume of distribution of 64.2 (50.6–81.0) L, intercompartment clearance of 116.4 (90.6–156.0) L/h, and peripheral volume of distribution of 167 (132–217) L. No serious adverse reactions or hemodynamic changes requiring the discontinuation of dexmedetomidine were observed. Dexmedetomidine had increased clearance and volume of distribution in mechanically ventilated children in ICU after neurosurgery, thereby indicating the need to adjust the dosage to obtain a target plasma concentration.

Pharmacometrics of clobazam in pediatrics: Prediction of effective clobazam doses for Dravet syndrome

OBJECTIVE

To describe the use of a population pharmacokinetic (PopPK) model incorporating weight and ontogeny to identify effective clobazam (CLB) dosing for use in a clinical trial in pediatric patients with Dravet syndrome.

METHODS

Pharmacokinetic data were combined from 3 CLB trials (OV-1012, OV-1017, and study 301) and a simulated study (study 401) for a total of 1306 CLB and 1305 N-desmethyl clobazam (N-CLB) samples from 193 Lennox-Gastaut syndrome patients and healthy subjects aged 6 months to 45 years. A structured approach based on US Food and Drug Administration guidance and pharmacometric knowledge discovery was developed using a nonlinear mixed-effects approach. Graphing and fitting using logistical weight regression were used to identify covariates for inclusion in the final model, which was evaluated using goodness-of-fit criteria and validated using prediction-corrected visual predictive check (pcVPC). Using the final PopPK model, a simulation study determined CLB and N-CLB distributions after 4 weeks of 1.5 and 2.0 mg/kg CLB.

RESULTS

The parameters of the final PopPK model were similar to previous reports. Fixed-effect parameters were precisely estimated, with no significant increase in NONMEM objective function value. Intersubject variability estimates were similar to previous reports, with <35% shrinkage associated with parameter variability, except for intercompartmental clearance and apparent volumes of distribution of peripheral compartments. Goodness-of-fit plots and pcVPC show that the model adequately described CLB and N-CLB data. The CLB/N-CLB ratio in virtual study subjects aged <3 years was 0.23 for 1.5 and 2.0 mg/kg and was 0.14 for subjects aged ≥3 years, which is 2 to 3 times those reported in a previous stiripentol/CLB/valproate study in which seizure improvement was reported.

SIGNIFICANCE

The PopPK model dosing parameters of 1.5 and 2.0 mg/kg are likely to result in efficacious concentrations of CLB and N-CLB in pediatric patients as young as 16 months. Dosages exceeding 1.5 mg/kg should be monitored for tolerability, particularly in patients aged <2 years, as there may be a higher incidence of sedation.

Model extrapolation to a real-world dataset: evaluation of tacrolimus population pharmacokinetics and drug interaction in pediatric liver transplantation patients

1. Numerous tacrolimus population pharmacokinetic (PPK) models in pediatric liver transplantation patients have been established to define an optimal dose schedule. However, the applicability of extrapolating these PPK models to our clinical center remains unknown. The goals of the present study was to evaluate model external predictiveness and establish a new model applicable to traditional therapeutic drug monitoring data.2. Published PPK models were collected from the literature and assessed using our real-world dataset including 41 pediatric liver transplantation patients via the individual prediction error method. The establishment of a new model was characterized using non-linear mixed-effects modeling.3. Nine published pediatric liver transplantation PPK models were identified, three of which could be applied to our real-world dataset. However, these models were dissatisfactory in terms of individual prediction error and hence, inadequate for extrapolation. Finally, a new model applicable to our real-world dataset was established as follows: CL/F = 22.9 × (WT/70)^0.75 × (1 - WZ × 0.264) × (1 - FCZ × 0.338) × (1 + ASPI × 0.281) × (POD/41)^0.0486 L/h; V/F = 906 × (WT/70) L. Where WT, WZ, FCZ, ASPI and POD were weight, Wuzhi capsule, fluconazole, aspirin and post-transplantation day, respectively. In conclusion, published models were inadequate for application to our real-world dataset. The present study produced a new model applicable to our real-world study data.

Population Pharmacokinetics and Dosing of Milrinone After Patent Ductus Arteriosus Ligation in Preterm Infants

OBJECTIVES

The postoperative course of patent ductus arteriosus ligation is often complicated by postligation cardiac syndrome, occurring in 10-45% of operated infants. Milrinone might prevent profound hemodynamic instability and improve the recovery of cardiac function in this setting. The present study aimed to describe the population pharmacokinetics of milrinone in premature neonates at risk of postligation cardiac syndrome and give dosing recommendations.

DESIGN

A prospective single group open-label pharmacokinetics study.

SETTINGS

Two tertiary care neonatal ICUs: Tallinn Children's Hospital and Tartu University Hospital, Estonia.

PATIENTS

Ten neonates with postmenstrual age of 24.6-30.1 weeks and postnatal age of 5-27 days undergoing patent ductus arteriosus ligation and at risk of postligation cardiac syndrome, based on echocardiographic assessment of left ventricular output of less than 200 mL/kg/min 1 hour after the surgery.

INTERVENTIONS

Milrinone at a dose of 0.73 μg/kg/min for 3 hours followed by 0.16 μg/kg/min for 21 hours. Four blood samples from each patient for milrinone plasma concentration measurements were collected.

MEASUREMENTS AND MAIN RESULTS

Concentration-time data of milrinone were analyzed with nonlinear mixed-effects modeling software (NONMEM Version 7.3 [ICON Development Solutions, Ellicott City, MD]). Probability of target attainment simulations gave a dosing schedule that maximally attains concentration targets of 150-250 μg/L. Milrinone pharmacokinetics was described by a one-compartmental linear model with allometric scaling to bodyweight and an age maturation function of glomerular filtration rate. Parameter estimates for a patient with the median weight were 0.350 (L/hr) for clearance and 0.329 (L) for volume of distribution. The best probability of target attainment was achieved with a loading dose of 0.50 μg/kg/min for 3 hours followed by 0.15 μg/kg/min (postmenstrual age < 27 wk) or 0.20 μg/kg/min (postmenstrual age ≥ 27 wk).

CONCLUSIONS

Population pharmacokinetic modeling and simulations suggest a slow loading dose followed by maintenance infusion to reach therapeutic milrinone plasma concentrations within the timeframe of the postligation cardiac syndrome.

Application of Size and Maturation Functions to Population Pharmacokinetic Modeling of Pediatric Patients

Traditionally, dosage for pediatric patients has been optimized using simple weight-scaled methods, but these methods do not always meet the requirements of children. To overcome this discrepancy, population pharmacokinetic (PK) modeling of size and maturation functions has been proposed. The main objective of the present study was to evaluate a new modeling method for pediatric patients using clinical data from three different clinical studies. To develop the PK models, a nonlinear mixed effect modeling method was employed, and to explore PK differences in pediatric patients, size with allometric and maturation with Michaelis–Menten type functions were evaluated. Goodness of fit plots, visual predictive check and bootstrap were used for model evaluation. Single application of size scaling to PK parameters was statistically significant for the over one year old group. On the other hand, simultaneous use of size and maturation functions was statistically significant for infants younger than one year old. In conclusion, population PK modeling for pediatric patients was successfully performed using clinical data. Size and maturation functions were applied according to established criteria, and single use of size function was applicable for over one year ages, while size and maturation functions were more effective for PK analysis of neonates and infants.

Population Pharmacokinetic and Exposure-Response Analyses of Prasugrel in Pediatric Patients with Sickle Cell Anemia

BACKGROUND AND OBJECTIVE

Prasugrel, a P2Y₁₂ adenosine diphosphate (ADP) receptor antagonist, inhibits ADP-mediated platelet activation and aggregation in patients with sickle cell anemia (SCA). We developed a population pharmacokinetic (popPK) model in pediatric patients from 2 to <18 years of age with SCA, and performed exposure-response evaluations to characterize the effects of prasugrel in a subset of these patients who weighed 19 kg or more and experienced at least two episodes of vaso-occlusive crises (VOC) in the past year.

METHODS

A three-compartment popPK model adapted from that used in adults with acute coronary syndrome was used to describe the relationship between plasma concentrations of prasugrel's active metabolite (Pras-AM) and time using data from phase II and III clinical studies in children. A VOC event rate model was developed from the phase III study to explore the exposure-response relationship between Pras-AM exposure and VOC, and included evaluation of covariates.

RESULTS

The final popPK model for children with SCA provided a reasonable fit to Pras-AM plasma concentrations over time, with estimates of apparent clearance (CL/F) (172 L/h) and apparent volume of distribution (Vd/F) (51.7 L) that were comparable to previous studies in adults. The final model included weight as a covariate on both CL/F and Vd/F, and age as a covariate on CL/F. Analyses of safety (bleeding events requiring medical intervention) and efficacy (VOC event rate) variables showed no apparent relationship to model-predicted Pras-AM exposure quartiles, and no statistically significant effects of intrinsic or extrinsic factors on the VOC event rate were identified in the VOC event rate model. The effect of post hoc exposures on the VOC event rate did not reach statistical significance.

CONCLUSIONS

A popPK model was developed that provided reasonable parameter estimates, goodness-of-fit diagnostics, and visual predictive checks when applied to Pras-AM plasma concentrations in pediatric patients with SCA. Post hoc exposures obtained from this model did not correlate with measures of VOC or bleeding events in this population.

Pharmacokinetics of Fentanyl and Its Derivatives in Children: A Comprehensive Review

Fentanyl and its derivatives sufentanil, alfentanil, and remifentanil are potent opioids. A comprehensive review of the use of fentanyl and its derivatives in the pediatric population was performed using the National Library of Medicine PubMed. Studies were included if they contained original pharmacokinetic parameters or models using established routes of administration in patients younger than 18 years of age. Of 372 retrieved articles, 44 eligible pharmacokinetic studies contained data of 821 patients younger than 18 years of age, including more than 46 preterm infants, 64 full-term neonates, 115 infants/toddlers, 188 children, and 28 adolescents. Underlying diagnoses included congenital heart and pulmonary disease and abdominal disorders. Routes of drug administration were intravenous, epidural, oral-transmucosal, intranasal, and transdermal. Despite extensive use in daily clinical practice, few studies have been performed. Preterm and term infants have lower clearance and protein binding. Pharmacokinetics was not altered by chronic renal or hepatic disease. Analyses of the pooled individual patients' data revealed that clearance maturation relating to body weight could be best described by the Hill function for sufentanil (R ² = 0.71, B _max 876 mL/min, K ₅₀ 16.3 kg) and alfentanil (R ² = 0.70, B _{max (fixed)} 420 mL/min, K ₅₀ 28 kg). The allometric exponent for estimation of clearance of sufentanil was 0.99 and 0.75 for alfentanil clearance. Maturation of remifentanil clearance was described by linear regression to bodyweight (R ² = 0.69). The allometric exponent for estimation of remifentanil clearance was 0.76. For fentanyl, linear regression showed only a weak correlation between clearance and bodyweight in preterm and term neonates (R ² = 0.22) owing to a lack of data in older age groups. A large heterogeneity regarding study design, clinical setting, drug administration, laboratory assays, and pharmacokinetic estimation was observed between studies introducing bias into the analyses performed in this review. A limitation of this review is that pharmacokinetic data, based on different modes of administration, dosing schemes, and parameter estimation methods, were combined.