Choosing the Allometric Exponent in Covariate Model Building

Pharmacokinetics and safety of rifapentine in children: dosing for latent tuberculosis infection

OBJECTIVES

To assess the safety of 4-week daily rifapentine-isoniazid regimen in latent tuberculosis for Chinese children, and to provide paediatric-specific evidence for extrapolating adult dosing strategies to children.

METHODS

An open-label, prospective, single-arm clinical trial was conducted among eligible patients (aged <10 years old). Rifapentine concentrations and laboratory safety biomarker (total bilirubin) were analysed and used for population pharmacokinetic-toxicity model development. Simulations were performed to compare efficacy and safety of weight-based and flat-dosing strategy.

RESULTS

Once-daily rifapentine treatment was well tolerated: 310 samples (rifapentine n = 139; total bilirubin n = 171) from 36 children (age range 0.89-10 years) were captured well by a joint one-compartment pharmacokinetic with time-varying clearance and an indirect response model. The model adequately described rifapentine autoinduction, reaching a plateau after 21 days and increasing clearance by 70.4%. Simulation suggested that weight-based dosing may cause underexposure in children under 14.5 kg. A flat-dosing strategy could ensure plasma levels within the therapeutic windows. Rifapentine's impact on total bilirubin was within a 2-fold range, and the effect subsided within 5 days after discontinuation.

CONCLUSIONS

Our study suggested that a flat-dosing strategy of rifapentine was potentially safe and effective for latent tuberculosis infection treatment in Chinese children aged 1 to 10 years old.

All You Need to Know About Allometric Scaling: An Integrative Review on the Theoretical Basis, Empirical Evidence, and Application in Human Pharmacology

Scaling approaches are used to describe or predict clearance for paediatric or obese populations from normal-weight adult values. Theoretical allometry assumes the existence of a universal bodyweight-based scaling relationship. Although theoretical allometry is highly disputed, it is commonly applied in pharmacological data analyses and clinical practice. The aim of the current review is to (1) increase pharmacologists' understanding of theoretical allometry to better understand the (implicit) assumptions and (dis)advantages and (2) highlight important methodological considerations with the application of this methodology. Theoretical allometry originated in an empirical, and later debated, observation by Kleiber of a scaling exponent of 0.75 between basal metabolic rate and body mass of mammals. The mathematical framework of West, Brown, and Enquist provides one possible explanation for this value. To date, multiple key assumptions of this framework have been disputed or disproven, and an increasing body of evidence is emerging against the existence of one universal allometric exponent. The promise of ease and universality of use that comes with theoretical approaches may be the reason they are so strongly sought after and defended. However, ecologists have suggested that the theory should move from a 'Newtonian approach', in which physical explanations are sought for a universal law and variability is of minor importance, to a 'Darwinian approach', in which variability is considered of primary importance for which evolutionary explanations can be found. No scientific support was found for the application of allometry for within-species scaling, so the application of basal metabolic rate-based scaling principles to clearance scaling remains unsubstantiated. Recent insights from physiologically based modelling approaches emphasise the interplay between drugs with different properties and physiological variables that underlie drug clearance, which drives the variability in the allometric scaling exponent in the field of pharmacology. To deal with this variability, drug-specific or patient-specific adaptations to theoretical allometric scaling are proposed, that introduce empiric elements and reduce the universality of the theory. The use of allometric scaling with an exponent of 0.75 may hold empirical merit for paediatric populations, except for the youngest individuals (aged ≤ 5 years). Nevertheless, biological interpretations and extrapolation potential attributed to models based on 0.75 allometric scaling are theoretically unfounded, and merits of the empirical application of this function should, as for all models, always be supported by appropriate model validation procedures. In this respect, it is not the value of the allometric exponent but the description and prediction of individual clearance values and drug concentrations that are of primary interest.

Population Pharmacokinetic Model of Intravenous Immunoglobulin in Patients Treated for Various Immune System Disorders

PURPOSE

Intravenous immunoglobulin (IVIG) is used to treat various immune system disorders, but the factors influencing its disposition are not well understood. This study aimed to estimate the population pharmacokinetic parameters of IVIG and to investigate the effect of genetic polymorphism of the FCGRT gene encoding the neonatal Fc receptor (FcRn) and clinical variability on the pharmacokinetic properties of IVIG in patients with immune system disorders.

METHODS

Patients were recruited from 4 hospitals in Malaysia. Clinical data were recorded, and blood samples were taken for pharmacokinetic and genetic studies. Population pharmacokinetic parameters were estimated by nonlinear mixed-effects modeling in Monolix. Age, weight, baseline immunoglobulin G concentration, ethnicity, sex, genotype, disease type, and comorbidity were investigated as potential covariates. Models were evaluated using the difference in objective function value, goodness-of-fit plots, visual predictive checks, and bootstrap analysis.

FINDINGS

A total of 292 blood samples were analyzed from 79 patients. The IVIG concentrations were best described by a 2-compartment model with linear elimination. Weight was found to be an important covariate for volume of distribution in the central compartment (Vc), volume of distribution in the peripheral compartment (Vp), and clearance in the central compartment, whereas disease type was found to be an important covariate for Vp. Goodness-of-fit plots indicated that the model fit the data adequately. Genetic polymorphism of the FCGRT gene encoding the neonatal Fc receptor did not affect the pharmacokinetic properties of IVIG.

IMPLICATIONS

This study supports the use of dosage based on weight as per current practice. The study findings highlight that Vp is significantly influenced by the type of disease being treated with IVIG. This relationship suggests that different disease types, particularly inflammatory and autoimmune conditions, may alter tissue permeability and fluid distribution due to varying degrees of inflammation. Increased inflammation can lead to enhanced permeability and retention of IVIG in peripheral tissues, reflecting higher Vp values.

Development and Quality Control of a Population Pharmacokinetic Model Library for Caspofungin

BACKGROUND

Caspofungin is an echinocandin antifungal agent commonly used as the first-line therapy for invasive candidiasis, salvage therapy for invasive aspergillosis, and empirical therapy for presumed fungal infections. Pharmacokinetic (PK) variabilities and suboptimal exposure have been reported for caspofungin, increasing the risk of insufficient efficacy.

OBJECTIVE

This work aimed to develop a caspofungin population pharmacokinetic (popPK) library and demonstrate its utility by assessing the probability of target attainment across diverse settings.

METHODS

We established a caspofungin popPK model library following a rigorous literature review, re-implementing selected models in R with rxode2. Quality control procedures included a comparison of different studies and assessing covariate impacts. Model libraries were primarily used to perform Monte Carlo simulations to estimate target attainment and guide personalized dosing in Candida infections.

RESULTS

A total of 13 models, one- or two-compartment models, were included. The most significant covariates were body size (weight and body surface area), liver function, and albumin level. The results show that children and adults showed considerable differences in pharmacokinetics. For C. albicans and C. parapsilosis, none of the populations achieved a PTA of ≥90% at their respective susceptible MIC values. In contrast, for C. glabrata, 70% of the adult studies reached a PTA of ≥90%, while all pediatric studies achieved the same PTA level.

CONCLUSION

At the recommended dosage, adult patients showed notably lower exposure to caspofungin compared to pediatric patients. Considering body size, liver function, and serum albumin is crucial when determining caspofungin dosage regimens. Furthermore, further research is required to comprehensively understand the pharmacokinetics of caspofungin in pediatric patients.

Tofacitinib pharmacokinetics in children and adolescents with juvenile idiopathic arthritis

These analyses characterized tofacitinib pharmacokinetics (PKs) in children and adolescents with juvenile idiopathic arthritis (JIA). Data were pooled from phase I (NCT01513902), phase III (NCT02592434), and open-label, long-term extension (NCT01500551) studies of tofacitinib tablet/solution (weight-based doses administered twice daily [b.i.d.]) in patients with JIA aged 2 to less than 18 years. Population PK modeling used a nonlinear mixed-effects approach, with covariates identified using stepwise forward-inclusion backward-deletion procedures. Simulations were performed to derive dosing recommendations for children and adolescents with JIA. Two hundred forty-six pediatric patients were included in the population PK model. A one-compartment model with first-order elimination and absorption with body weight as a covariate for oral clearance and apparent volume of distribution sufficiently described the data. Oral solution was associated with comparable average concentration (Cavg) and slightly higher (113.9%) maximum concentration (Cmax) versus tablet, which was confirmed by a subsequent randomized, open-label, bioavailability study conducted in healthy adult participants (n = 12) by demonstrating adjusted geometric mean ratios (90% confidence interval) between oral solution and tablet of 1.04 (1.00-1.09) and 1.10 (1.00-1.21) for area under the curve extrapolated to infinity and Cmax, respectively (NCT04111614). A dosing regimen of 3.2 mg b.i.d. solution in patients 10 to less than 20 kg, 4 mg b.i.d. solution in patients 20 to less than 40 kg, and 5 mg b.i.d. tablet/solution in patients greater than or equal to 40 kg, irrespective of age, was proposed to achieve constant Cavg across weight groups. In summary, population PK characterization informed a simplified tofacitinib dosing regimen that has been implemented in pediatric patients with JIA.

Integrative population pharmacokinetic/pharmacodynamic analysis of nemonoxacin capsule in Chinese patients with community-acquired pneumonia

Introduction: Nemonoxacin is an innovative quinolone antibiotic for treatment of community-acquired pneumonia (CAP). As more data are available from clinical studies, it is necessary to perform an integrative pharmacokinetic/pharmacodynamic (PK/PD) analysis to support and justify the optimal dosing regimen of nemonoxacin in clinical practice. Methods and Results: We developed a population PK model using non-linear mixed effect model based on the data of 195 Chinese subjects receiving nemonoxacin in phase I to III clinical trials. The base model was a standard two-compartment PK model defined by clearance (12 L/h) and central volume of distribution (86 L). Covariates included creatinine clearance (CL_cr), body weight (BW), sex, disease status and food. Compared to the subject with BW 60 kg, C_max and A U C 0 - 24 ,   ss reduced by 24% and 19% in the subject with BW 80 kg, respectively. Compared to the subject with CL_cr 150 ml/min, A U C 0 - 24 ,   ss and T_1/2 increased by 28% and 24%, respectively in the subject with CL_cr 30 ml/min. Compared to the fasted status, T_max of nemonoxacin increased by 1.2 h in the subject with fed status. Effects of sex and disease status on PK parameters were small (change of PK parameters ≤19%). AUC_0-24/MIC and %T > MIC were identified as the optimal PK/PD indices for predicting clinical efficacy. The AUC_0-24/MIC target was 63.3, 97.8, and 115.7 against Streptococcus pneumoniae, Staphylococcus aureus, and Haemophilus influenzae, respectively. The %T > MIC target was 7.96% against Klebsiella pneumoniae. Monte Carlo simulation showed that treatment with nemonoxacin 500 mg q24 h could attain a PK/PD cutoff value higher than the MIC₉₀ against S. pneumoniae and S. aureus. The corresponding cumulative fraction of response (CFR) was greater than 93%, while nemonoxacin 750 mg q24 h would provide higher PK/PD cutoff value against Haemophilus parainfluenzae, and higher CFR (83%) than 500 mg q24 h. Conclusion: Integrative PK/PD analysis justifies the reliable clinical and microbiological efficacy of nemonoxacin 500 mg q24 h in treating CAP caused by S. pneumoniae, S. aureus, and K. pneumoniae, irrespective of patient sex, mild renal impairment, empty stomach or not. However, nemonoxacin 750 mg q24 h would provide better efficacy than 500 mg q24 h for the CAP caused by H. parainfluenzae in terms of CFR.

Population Pharmacokinetic Modelling of Intravenous Immunoglobulin Treatment in Patients with Guillain-Barré Syndrome

Background and Objective

Intravenous immunoglobulin (IVIg) at a standard dosage is the treatment of choice for Guillain–Barré syndrome. The pharmacokinetics, however, is highly variable between patients, and a rapid clearance of IVIg is associated with poor recovery. We aimed to develop a model to predict the pharmacokinetics of a standard 5-day IVIg course (0.4 g/kg/day) in patients with Guillain–Barré syndrome.

Methods

Non-linear mixed-effects modelling software (NONMEM^®) was used to construct a pharmacokinetic model based on a model-building cohort of 177 patients with Guillain–Barré syndrome, with a total of 589 sequential serum samples tested for total immunoglobulin G (IgG) levels, and evaluated on an independent validation cohort that consisted of 177 patients with Guillain–Barré syndrome with 689 sequential serum samples.

Results

The final two-compartment model accurately described the daily increment in serum IgG levels during a standard IVIg course; the initial rapid fall and then a gradual decline to steady-state levels thereafter. The covariates that increased IgG clearance were a more severe disease (as indicated by the Guillain–Barré syndrome disability score) and concomitant methylprednisolone treatment. When the current dosing regimen was simulated, the percentage of patients who reached a target ∆IgG > 7.3 g/L at 2 weeks decreased from 74% in mildly affected patients to only 33% in the most severely affected and mechanically ventilated patients (Guillain–Barré syndrome disability score of 5).

Conclusions

This is the first population-pharmacokinetic model for standard IVIg treatment in Guillain–Barré syndrome. The model provides a new tool to predict the pharmacokinetics of alternative regimens of IVIg in Guillain–Barré syndrome to design future trials and personalise treatment.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40262-022-01136-z.

Characterizing Pharmacokinetics in Children With Obesity-Physiological, Drug, Patient, and Methodological Considerations

Childhood obesity is an alarming public health problem. The pediatric obesity rate has quadrupled in the past 30 years, and currently nearly 20% of United States children and 9% of children worldwide are classified as obese. Drug distribution and elimination processes, which determine drug exposure (and thus dosing), can vary significantly between patients with and without obesity. Obesity-related physiological changes, such as increased tissue volume and perfusion, altered blood protein concentrations, and tissue composition can greatly affect a drug's volume of distribution, which might necessitate adjustment in loading doses. Obesity-related changes in the drug eliminating organs, such as altered enzyme activity in the liver and glomerular filtration rate, can affect the rate of drug elimination, which may warrant an adjustment in the maintenance dosing rate. Although weight-based dosing (i.e., in mg/kg) is commonly practiced in pediatrics, choice of the right body size metric (e.g., total body weight, lean body weight, body surface area, etc.) for dosing children with obesity still remains a question. To address this gap, the interplay between obesity-related physiological changes (e.g., altered organ size, composition, and function), and drug-specific properties (e.g., lipophilicity and elimination pathway) needs to be characterized in a quantitative framework. Additionally, methodological considerations, such as adequate sample size and optimal sampling scheme, should also be considered to ensure accurate and precise top-down covariate selection, particularly when designing opportunistic studies in pediatric drug development. Further factors affecting dosing, including existing dosing recommendations, target therapeutic ranges, dose capping, and formulations constraints, are also important to consider when undergoing dose selection for children with obesity. Opportunities to bridge the dosing knowledge gap in children with obesity include modeling and simulating techniques (i.e., population pharmacokinetic and physiologically-based pharmacokinetic [PBPK] modeling), opportunistic clinical data, and real world data. In this review, key considerations related to physiology, drug parameters, patient factors, and methodology that need to be accounted for while studying the influence of obesity on pharmacokinetics in children are highlighted and discussed. Future studies will need to leverage these modeling opportunities to better describe drug exposure in children with obesity as the childhood obesity epidemic continues.

Allometric Scaling in Pharmacokinetic Studies in Anesthesiology

A clinical review is presented of basic allometric scaling theory and its application to pharmacokinetic models in anesthesia and other fields in the biologic sciences.

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

INTRODUCTION

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

AREAS COVERED

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

EXPERT OPINION

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

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.

Pharmacokinetic modeling of R and S-Methadone and their metabolites to study the effects of various covariates in post-operative children

Methadone is a synthetic opioid used as an analgesic and for the treatment of opioid abuse disorder. The analgesic dose in the pediatric population is not well-defined. The pharmacokinetics (PKs) of methadone is highly variable due to the variability in alpha-1 acid glycoprotein (AAG) and genotypic differences in drug-metabolizing enzymes. Additionally, the R and S enantiomers of methadone have unique PK and pharmacodynamic properties. This study aims to describe the PKs of R and S methadone and its metabolite 2-ethylidene-1,5-dimethyl-3,3-diphenylpyrrolidine (EDDP) in pediatric surgical patients and to identify sources of inter- and intra-individual variability. Children aged 8-17.9 years undergoing orthopedic surgeries received intravenous methadone 0.1 mg/kg intra-operatively followed by oral methadone 0.1 mg/kg postoperatively every 12 h. Pharmacokinetics of R and S methadone and EDDP were determined using liquid chromatography tandem mass spectrometry assays and the data were modeled using nonlinear mixed-effects modeling in NONMEM. R and S methadone PKs were well-described by two-compartment disposition models with first-order absorption and elimination. EDDP metabolites were described by one compartment disposition models with first order elimination. Clearance of both R and S methadone were allometrically scaled by bodyweight. CYP2B6 phenotype was a determinant of the clearance of both the enantiomers in an additive gene model. The intronic CYP3A4 single-nucleotide polymorphism (SNP) rs2246709 was associated with decreased clearance of R and S methadone. Concentrations of AAG and the SNP of AAG rs17650 independently increased the volume of distribution of both the enantiomers. The knowledge of these important covariates will aid in the optimal dosing of methadone in children.

Population pharmacokinetics of unbound valproic acid in pediatric epilepsy patients in China: a protein binding model

PURPOSE

The purpose of this study was to establish a protein binding model of unbound valproic acid (VPA) based on Chinese pediatric patients with epilepsy and provide a reference for clinical medication.

METHODS

A total of 313 patients were included and both their total and unbound VPA concentrations (375 pairs of concentrations) were measured. NONMEM software was used for population pharmacokinetic modeling. The stepwise method was used to screen the potential covariates. Goodness-of-fit plot, bootstrap, and visual predictive check were used for model evaluation. In addition, dose recommendations for typical patients aged 0 to 16 years were proposed by Monte Carlo simulations.

RESULTS

A one-compartment model of first-order absorption and first-order elimination was used to describe the pharmacokinetic characteristics of unbound VPA, and the linear non-saturable binding equation was introduced to describe the protein binding. Body weight, age-based maturation, and co-medicated with lamotrigine could affect the CL/F of unbound and bound VPA. Model evaluation showed satisfactory robustness of the final model. The dosing regimens for children aged 0 to 16 years were proposed based on the final established model.

CONCLUSION

We developed a population pharmacokinetic model of unbound and bound VPA that took account of protein binding. The VPA dosing regimen in pediatric patients with epilepsy needs to be optimized by the body weight, age, and co-medications.

Population Pharmacokinetics of Clotting Factor Concentrates and Desmopressin in Hemophilia

Hemophilia A and B are bleeding disorders caused by a deficiency of clotting factor VIII and IX, respectively. Patients with severe hemophilia (< 0.01 IU mL⁻¹) and some patients with moderate hemophilia (0.01–0.05 IU mL⁻¹) administer clotting factor concentrates prophylactically. Desmopressin (d-amino d-arginine vasopressin) can be applied in patients with non-severe hemophilia A. The aim of administration of factor concentrates or desmopressin is the prevention or cessation of bleeding. Despite weight-based dosing, it has been demonstrated that factor concentrates still exhibit considerable pharmacokinetic variability. Population pharmacokinetic analyses, in which this variability is quantified and explained, are increasingly performed in hemophilia research. These analyses can assist in the identification of important patient characteristics and can be applied to perform patient-tailored dosing. This review aims to present and discuss the population pharmacokinetic analyses that have been conducted to develop population pharmacokinetic models describing factor levels after administration of factor VIII or factor IX concentrates or d-amino d-arginine vasopressin. In total, 33 publications were retrieved from the literature. Two approaches were applied to perform population pharmacokinetic analyses, the standard two-stage approach and non-linear mixed-effect modeling. Using the standard two-stage approach, four population pharmacokinetic models were established describing factor VIII levels. In the remaining 29 analyses, the non-linear mixed-effect modeling approach was applied. NONMEM was the preferred software to establish population pharmacokinetic models. In total, 18 population pharmacokinetic analyses were conducted on the basis of data from a single product. From all available population pharmacokinetic analyses, 27 studies also included data from pediatric patients. In the majority of the population pharmacokinetic models, the population pharmacokinetic parameters were allometrically scaled using actual body weight. In this review, the available methods used for constructing the models, key features of these models, patient population characteristics, and established covariate relationships are described in detail.

Population pharmacokinetics of sildenafil in extremely premature infants

AIMS

To characterize the population pharmacokinetics (PK) of sildenafil and its active metabolite, N-desmethyl sildenafil (DMS), in premature infants.

METHODS

We performed a multicentre, open-label trial to characterize the PK of sildenafil in infants ≤28 weeks gestation and < 365 postnatal days (cohort 1) or < 32 weeks gestation and 3-42 postnatal days (cohort 2). In cohort 1, we obtained PK samples from infants receiving sildenafil as ordered per the local standard of care (intravenous [IV] or enteral). In cohort 2, we administered a single IV dose of sildenafil and performed PK sampling. We performed a population PK analysis and dose-exposure simulations using the software NONMEM®.

RESULTS

We enrolled 34 infants (cohort 1 n = 25; cohort 2 n = 9) and collected 109 plasma PK samples. Sildenafil was given enterally (0.42-2.09 mg/kg) in 24 infants in cohort 1 and via IV (0.125 or 0.25 mg/kg) in all infants in cohort 2. A 2-compartment PK model for sildenafil and 1-compartment model for DMS, with presystemic conversion of sildenafil to DMS, characterized the data well. Coadministration of fluconazole (n = 4), a CYP3A inhibitor, resulted in an estimated 59% decrease in sildenafil clearance. IV doses of 0.125, 0.5 and 1 mg/kg every 8 hours (in the absence of fluconazole) resulted in steady-state maximum sildenafil concentrations that were generally within the range of those reported to inhibit phosphodiesterase type 5 activity in vitro.

CONCLUSIONS

We successfully characterized the PK of sildenafil and DMS in premature infants and applied the model to inform dosing for a follow-up, phase II study.

Population Pharmacokinetic/Pharmacodynamic Modeling of Methylprednisolone in Neonates Undergoing Cardiopulmonary Bypass

Methylprednisolone is used in neonates to modulate cardiopulmonary bypass (CPB)–induced inflammation, but optimal dosing and exposure are unknown. We used plasma methylprednisolone and interleukin (IL)‐6 and IL‐10 concentrations from neonates enrolled in a randomized trial comparing one vs. two doses of methylprednisolone to develop indirect response population pharmacokinetic/pharmacodynamic models characterizing the exposure–response relationships. We applied the models to simulate methylprednisolone dosages resulting in the desired IL‐6 and ‐10 exposures, known mediators of CPB‐induced inflammation. A total of 64 neonates (median weight 3.2 kg, range 2.2–4.3) contributed 290 plasma methylprednisolone concentrations (range 1.07–12,700 ng/mL) and IL‐6 (0–681 pg/mL) and IL‐10 (0.1–1125 pg/mL). Methylprednisolone plasma exposure following a single 10 mg/kg intravenous dose inhibited IL‐6 and stimulated IL‐10 production when compared with placebo. Higher (30 mg/kg) or more frequent (twice) dosing did not confer additional benefit. Clinical efficacy studies are needed to evaluate the effect of optimized dosing on outcomes.

Population Pharmacokinetics of Doxycycline in Children

Doxycycline is a tetracycline-class antimicrobial labeled by the United States (U.S.) Food and Drug Administration for children >8 years of age for many common childhood infections. Doxycycline is not labeled for children ≤8 years of age, due to the association between tetracycline class antibiotics and tooth staining, although doxycycline may be used off-label in severe conditions. Accordingly, there is a paucity of pharmacokinetic (PK) data to guide dosing in children 8 years and younger. We leveraged opportunistically-collected plasma samples after intravenous (IV) and oral doxycycline doses received per standard of care to characterize the PK of doxycycline in children of different ages, and evaluated the effect of obesity and fasting status on PK parameters.We developed a population PK model of doxycycline using data collected from 47 patients 0-18 years of age, including 14 participants ≤8 years. We developed a 1 compartment PK model and found doxycycline clearance to be 3.32 L/h/70 kg and volume to be 96.8 L/70kg for all patients; comparable to values reported in adults. We estimated a bioavailability of 89.6%, also consistent with adult data. Allometrically scaled clearance and volume of distribution did not differ between children 2 to ≤8 years of age and children >8 to ≤18 years of age, suggesting that younger children may be given the same per kg dosing. Obese and fasting status were not selected for inclusion in the final model. Additional doxycycline PK samples collected in future studies may be used to improve model performance and maximize its clinical value.