Increased regulatory demands for pediatric drug development research have fostered interest in the use of modeling and simulation among industry and academia. Physiologically based pharmacokinetic (PBPK) modeling offers a unique modality to incorporate multiple levels of information to estimate age-specific pharmacokinetics. This tutorial will serve to provide the reader with a basic understanding of the procedural steps to developing a pediatric PBPK model and facilitate a discussion of the advantages and limitations of this modeling technique.

Abstract. The use of physiologically based pharmacokinetic (PBPK) models in the field of pediatric drug development has garnered much interest of late due to a recent Food and Drug Administration recommendation. The purpose of this study is to illustrate the developmental processes involved in creation of a pediatric PBPK model incorporating existing adult drug data. Lorazepam, a benzodiazepine utilized in both adults and children, was used as an example. A population-PBPK model was developed in PK-Sim v4.2® and scaled to account for age-related changes in size and composition of tissue compartments, protein binding, and growth/maturation of elimination processes. Dose (milligrams per kilogram) requirements for children aged 0-18 years were calculated based on simulations that achieved targeted exposures based on adult references. Predictive accuracy of the PBPK model for producing comparable plasma concentrations among 63 pediatric subjects was assessed using average-fold error (AFE). Estimates of clearance (CL) and volume of distribution (V ss ) were compared with observed values for a subset of 15 children using fold error (FE). Pediatric dose requirements in young children (1-3 years) exceeded adult levels on a linear weight-adjusted (milligrams per kilogram) basis. AFE values for model-derived concentration estimates were within 1.5-and 2-fold deviation from observed values for 73% and 92% of patients, respectively. For CL, 60% and 80% of predictions were within 1.5 and 2 FE, respectively. Comparatively, predictions of V ss were more accurate with 80% and 100% of estimates within 1.5 and 2 FE, respectively. Using the presented workflow, the developed pediatric model estimated lorazepam pharmacokinetics in children as a function of age.

Developmental changes in GI fluid composition can result in relevant discrepancies in luminal compound solubility between children and adults.

; for the Best Pharmaceuticals for Children Act-Pediatric Trials Network Steering Committee IMPORTANCE Children of all ages appear susceptible to severe acute respiratory syndrome coronavirus 2 infection. To support pediatric clinical studies for investigational treatments of coronavirus disease 2019 (COVID-19), pediatric-specific dosing is required. OBJECTIVE To define pediatric-specific dosing regimens for hydroxychloroquine and remdesivir for COVID-19 treatment. DESIGN, SETTING, AND PARTICIPANTS Pharmacokinetic modeling and simulation were used to extrapolate investigated adult dosages toward children (March 2020-April 2020). Physiologically based pharmacokinetic modeling was used to inform pediatric dosing for hydroxychloroquine. For remdesivir, pediatric dosages were derived using allometric-scaling with age-dependent exponents. Dosing simulations were conducted using simulated pediatric and adult participants based on the demographics of a white US population. INTERVENTIONS Simulated drug exposures following a 5-day course of hydroxychloroquine (400 mg every 12 hours × 2 doses followed by 200 mg every 12 hours × 8 doses) and a single 200-mg intravenous dose of remdesivir were computed for simulated adult participants. A simulation-based dose-ranging study was conducted in simulated children exploring different absolute and weight-normalized dosing strategies. MAIN OUTCOMES AND MEASURES The primary outcome for hydroxychloroquine was average unbound plasma concentrations for 5 treatment days. Additionally, unbound interstitial lung concentrations were simulated. For remdesivir, the primary outcome was plasma exposure (area under the curve, 0 to infinity) following single-dose administration. RESULTS For hydroxychloroquine, the physiologically based pharmacokinetic model analysis included 500 and 600 simulated white adult and pediatric participants, respectively, and supported weight-normalized dosing for children weighing less than 50 kg. Geometric mean-simulated average unbound plasma concentration values among children within different developmental age groups (32-35 ng/mL) were congruent to adults (32 ng/mL). Simulated unbound hydroxychloroquine concentrations in lung interstitial fluid mirrored those in unbound plasma and were notably lower than in vitro concentrations needed to mediate antiviral activity. For remdesivir, the analysis included 1000 and 6000 simulated adult and pediatric participants, respectively. The proposed pediatric dosing strategy supported weight-normalized dosing for participants weighing less than 60 kg. Geometric mean-simulated plasma area under the time curve 0 to infinity values among children within different developmental age-groups (4315-5027 ng × h/mL) were similar to adults (4398 ng × h/mL). CONCLUSIONS AND RELEVANCE This analysis provides pediatric-specific dosing suggestions for hydroxychloroquine and remdesivir and raises concerns regarding hydroxychloroquine use for COVID-19 treatment because concentrations were less than those needed to mediate an antiviral effect.

The small intestine represents the region where the majority of drug and nutrient absorption transpires. Among adults, small intestinal transit kinetics is well delineated; however, the applicability of these values toward children remains unclear. This article serves to examine the relationship between age and mean small intestinal transit time (SITT) based on the available literature. In addition, the influence of alterations in intestinal transit time was explored among children using a model-based approach. Primary literature sources depicting SITT from children to adults were ascertained via the PubMed database. Data were limited to subjects without pathologies that could influence intestinal motility. Random-effect meta-regression models with between-study variability were employed to assess the influence of age on SITT. Three separate models with age as a linear or higher-order (i.e., second-and third-order polynomial) regressor were implemented to assess for the potential of both linear and curvilinear relationships. Examination of the influence of altered intestinal transit kinetics on the absorption of a sustained release theophylline preparation was explored among children between 8 and 14 years using physiologically based pharmacokinetic (PBPK) modeling. Age was not found to be a significant modulator of small intestinal transit within either the linear or higher-order polynomial meta-regression models. PBPK simulations indicated a lack of influence of variations in SITT on the absorption of theophylline from the examined sustained release formulation in older children. Based on the current literature, there is no evidence to suggest that mean SITT differs between children and adults.

Differences in fentanyl pharmacokinetics (PK) between obese and nonobese adults have previously been reported; however, the impact of childhood obesity on fentanyl PK is relatively unknown. We developed a population pharmacokinetic (PopPK) model using opportunistically collected samples from a cohort of predominately obese children receiving fentanyl per the standard of care. Using a probability‐based approach, we evaluated the ability of different continuous infusion strategies to provide steady‐state concentrations (Css) within an analgesic concentration range (1‐3 ng/mL). Fifty‐three samples from 32 children were used for PopPK model development. Median (range) age and body weight of study participants were 13 years (2‐19 years) and 52 kg (16‐164 kg), respectively. The majority of children (94%) were obese. A 2‐compartment model allometrically scaled by total body weight provided an appropriate fit to the data. Estimated typical clearance was 32.5 L/h (scaled to 70 kg). A fixed dose rate infusion of 1 µg/kg/h was associated with probabilities between 49% and 58% for achieving Css within target; however, the risk of achieving Css > 3 ng/mL increased with increasing body weight (15% at 16 kg vs 43% at 164 kg). A proposed model‐based infusion strategy maintained consistent probabilities across the examined weight range for achieving Css within (58%) and above (20%) target. Use of an allometric relationship between weight and clearance was appropriate for describing the PK of intravenous fentanyl in our cohort of predominately obese children. Our proposed model‐derived continuous infusion strategy maximized the probability of achieving target Css in children of varying weights.

Objective To characterize hydroxychloroquine exposure in patients with rheumatic disease receiving long-term hydroxychloroquine compared to target concentrations with reported antiviral activity against the 2019 coronavirus SARS-CoV-2. Methods We evaluated total hydroxychloroquine concentrations in serum and plasma from published literature values, frozen serum samples from a pediatric lupus trial, and simulated concentrations using a published pharmacokinetic model during pregnancy. For each source, we compared observed or predicted hydroxychloroquine concentrations to target concentrations with reported antiviral activity against SARS-CoV-2. Results The average total serum/plasma hydroxychloroquine concentrations were below the lowest SARS-CoV-2 target of 0.48 mg/L in all studies. Assuming the highest antiviral target exposure (total plasma concentration of 4.1 mg/L), all studies had approximately one-tenth the necessary concentration for in-vitro viral inhibition. Pharmacokinetic model simulations confirmed that pregnant adults receiving common dosing for rheumatic diseases did not achieve target exposures; however, the models predict that a dosage of 600 mg once a day during pregnancy would obtain the lowest median target exposure for most patients after the first dose. Conclusion We found that the average patient receiving treatment with hydroxychloroquine for rheumatic diseases, including children and non-pregnant/pregnant adults, are unlikely to achieve total serum or plasma concentrations shown to inhibit SARS-CoV-2 in-vitro. Nevertheless, patients receiving hydroxychloroquine long-term may have tissue concentrations far exceeding that of serum/plasma. Because the therapeutic window for hydroxychloroquine in the setting of SARS-CoV-2 is unknown, well-designed clinical trials that include patients with rheumatic disease are urgently needed to characterize the efficacy, safety, and target exposures for hydroxychloroquine.

The current investigation depicts a proficient modality for estimation of protein binding in pediatrics and will, therefore, aid in reducing uncertainty associated with pediatric pharmacokinetic predictions.