Meta-analysis of hepatic cytochrome P450 ontogeny to underwrite the prediction of pediatric pharmacokinetics using physiologically based pharmacokinetic modeling

Dosing of Venetoclax in Pediatric Patients with Relapsed Acute Myeloid Leukemia: Analysis of Developmental Pharmacokinetics and Exposure-Response Relationships

PURPOSE

This work aimed to characterize the pharmacokinetics and exposure-response relationships of venetoclax in pediatric patients with relapsed or refractory (R/R) acute myeloid leukemia (AML) to identify venetoclax doses to be administered to pediatric patients in the phase 3 study.

METHODS

Data from 121 patients across three phase 1 studies enrolling pediatric patients with R/R malignancies were utilized to develop a population pharmacokinetic model to describe venetoclax pharmacokinetics in pediatric patients. Individual patient average venetoclax plasma concentration up to the event of interest, derived based on the population pharmacokinetics analysis, was used to evaluate the exposure-response relationships to efficacy (complete response) and safety (neutropenia and thrombocytopenia) endpoints for patients with AML who received venetoclax in combination with azacitidine, decitabine, or cytarabine (n = 36). The population pharmacokinetic model was then used to simulate exposures in pediatric age- and weight-based subgroups to identify the venetoclax doses for pediatric patients.

FINDINGS

The pharmacokinetic data were adequately described by the two-compartment population pharmacokinetic model with first-order absorption and elimination. The model accounted for cytochrome P450 3A developmental changes using a maturation function and incorporated allometric scaling to account for growth and body size effect. Weight was identified as a statistically significant covariate on clearance and volume of distribution and retained in the final model. Population pharmacokinetic estimates were comparable to previously reported estimates in adults. Exposure-response analyses suggested that the clinical efficacy of venetoclax in combination with high-dose cytarabine (HDAC) is maximized at 600 mg adult-equivalent, and higher doses are unlikely to enhance clinical efficacy. Venetoclax 600 mg adult-equivalent was selected for further development in combination with HDAC. Additionally, venetoclax 400 mg adult-equivalent was selected for bridging/maintenance therapy in combination with azacitidine. Flat exposure-response relationships were observed with Grade ≥3 neutropenia and thrombocytopenia. Doses were selected based on weight (allometric scaling) for children aged ≥2 years old and based on weight and CYP3A ontogeny for children aged <2 years. The selected age- and weight-based dosing scheme of venetoclax is projected to achieve venetoclax exposures in pediatric subgroups comparable to those observed in adults receiving venetoclax 400 mg or 600 mg.

IMPLICATIONS

This work characterized the pharmacokinetics and exposure-response relationships of venetoclax in pediatric patients and guided the selection of pediatric dosing regimens in support of the venetoclax phase 3 trial in pediatric AML (NCT05183035).

CLINICAL STUDIES

NCT03236857, NCT03181126, and NCT03194932.

Age-Dependent Abundance of CYP450 Enzymes Involved in Metronidazole Metabolism: Application to Pediatric PBPK Modeling

The expression of cytochrome P450 (CYP) enzymes is highly variable and associated with factors, such as age, genotype, sex, and disease states. In this study, quantification of metronidazole metabolizing CYP isoforms (CYP2A6, CYP2E1, CYP3A4, CYP3A5, and CYP3A7) in human liver microsomes from 115 children and 35 adults was performed using a quantitative proteomics method. The data confirmed age-dependent increase in CYP2A6, CYP2E1, and CYP3A4 abundance, whereas, as expected, CYP3A7 abundance showed postnatal decrease with age. In particular, the fold difference (neonatal to adulthood levels) in the protein abundance of CYP2A6, CYP2E1, and CYP3A4 was 14, 11, and 20, respectively. In contrast, protein abundance of CYP3A7 was > 125-fold higher in the liver microsomes of neonates than of adults. The abundance of CYP2A6 and CYP3A5 was associated with genotypes, rs4803381 and rs776746, respectively. A proteomics-informed physiologically based pharmacokinetic (PBPK) model was developed to describe the pharmacokinetics of metronidazole and its primary metabolite, 2-hydroxymethylmetronidazole. The model revealed an increase in the metabolite-to-parent ratio with age and showed a strong correlation between CYP2A6 abundance and metabolite formation (r 2 = 0.75). Notably, the estimated contribution of CYP3A7 was ~ 75% in metronidazole clearance in neonates. These data suggest that variability in CYP2A6 and CYP3A7 in younger children poses the risk of variable pharmacokinetics of metronidazole and its active metabolite with a potential impact on drug efficacy and safety. No sex-dependent difference was observed in the protein abundance of the studied CYPs. The successful integration of hepatic CYP ontogeny data derived from a large liver bank into the pediatric PBPK model of metronidazole can be extended to other drugs metabolized by the studied CYPs.

Development and Evaluation of Ontogeny Functions of the Major UDP-Glucuronosyltransferase Enzymes to Underwrite Physiologically Based Pharmacokinetic Modeling in Pediatric Populations

Uridine 5'-diphospho-glucuronosyltransferases (UGTs) demonstrate variable expression in the pediatric population. Thus, understanding of age-dependent maturation of UGTs is critical for accurate pediatric pharmacokinetics (PK) prediction of drugs that are susceptible for glucuronidation. Ontogeny functions of major UGTs have been previously developed and reported. However, those ontogeny functions are based on in vitro data (i.e., enzyme abundance, in vitro substrate activity, and so on) and therefore, may not translate to in vivo maturation of UGTs in the clinical setting. This report describes meta-analysis of the literature to develop and compare ontogeny functions for 8 primary UGTs (UGT1A1, UGT1A4, UGT1A6, UGT1A9, UGT2B7, UGT2B10, UGT2B15, and UGT2B17) based on published in vitro and in vivo studies. Once integrated with physiologically based pharmacokinetics modeling models, in vivo activity-based ontogeny functions demonstrated somewhat greater prediction accuracy (mean squared error, MSE: 0.05) compared to in vitro activity (MSE: 0.104) and in vitro abundance-based ontogeny functions (MSE: 0.129).

The potential of enteroids derived from children and adults to study age-dependent differences in intestinal CYP3A4/5 metabolism

Drug metabolism in the intestinal wall affects bioavailability of orally administered drugs and is influenced by age. Hence, it is important to fully understand the drug metabolizing capacity of the gut to predict systemic exposure. The aim of this study was to investigate the potential of enteroids as a tool to study CYP3A4/5 -mediated metabolism in both children and adults. Bioconversion of midazolam, a CYP3A4/5 model substrate, was studied using enteroid monolayers as well as tissue explants in the Ussing chamber, both derived from pediatric [median (range age): 54 weeks (2 days - 13 years), n = 21] and adult (n = 5) tissue. Caco-2 cellular monolayers were employed as controls. In addition, mRNA expression of CYP3A4 was determined in enteroid monolayers (n = 11), tissue (n = 23) and Caco-2 using RT-qPCR. Midazolam metabolism was successfully detected in all enteroid monolayers, as well as in all tissue explants studied in the Ussing chamber, whereas Caco-2 showed no significant metabolite formation. The extracted fraction of midazolam was similar between enteroid monolayers and tissue. The fraction of midazolam extracted increased with age in enteroid monolayers derived from 0 to 70 week old donors. No statistically significant correlation was observed in tissue likely due to high variability observed and the smaller donor numbers included in the study. At the level of gene expression, CYP3A4 increased with age in tissues (n = 32), while this was not reflected in enteroid monolayers (n = 16). Notably, asymmetric metabolite formation was observed in enteroids and tissue, with higher metabolite formation on the luminal side of the barrier. In summary, we demonstrated that enteroids can be used to measure CYP3A4/5 midazolam metabolism, which we show is similar as observed in fresh isolated tissue. This was the case both in children and adults, indicating the potential of enteroids to predict intestinal metabolism. This study provides promising data to further develop enteroids to study drug metabolism in vitro and potentially predict oral absorption for special populations as an alternative to using fresh tissue.

Prediction of Pediatric Pharmacokinetics for CYP3A4 Metabolized Drugs: Comparison of the Performance of Two Hepatic Ontogeny Within a Physiologically Based Pharmacokinetic Model

The rapid growth in the use of pediatric physiologically based pharmacokinetic (PBPK) models, particularly for regulatory applications, has focused emphasis on model verification and ensuring system parameters are robust, including how these change with age. Uncertainty remains regarding the ontogeny of some enzymes and transporters, in this study 2 published ontogeny profiles for hepatic CYP3A4 were compared. Clinical pharmacokinetic data on 4 intravenously administered CYP3A4 substrates (alfentanil, fentanyl, midazolam, and sildenafil) used across the pediatric age range was collected from the literature. The PBPK models were verified in the adult population and then used to compare the Salem and a modified Upreti ontogeny profiles for CYP3A4 in terms of parent drug clearance and area under the curve from birth onward. Overall, the modified Upreti ontogeny profile resulted in 15 out of 17 age-related predictions within 2-fold and 12 out of 17 predictions within 1.5-fold ranges of observed values, for the Salem ontogeny these values were 12 out of 17 and 8 out of 17, respectively. The Upreti ontogeny profile performed better than Salem, average fold error and absolute average fold error were 1.14 and 1.35 compared to 1.56 and 1.90, respectively. Identifying the optimal CYP3A4 ontogeny is important for regulatory use of PBPK especially given the number of drugs cleared by this enzyme. This study broadens the evidence from previous studies that Upreti is more favorable than Salem, but further work is needed especially in the neonatal and early infant age range.

Evaluation of Empiric Voriconazole Dosing and Therapeutic Drug Monitoring in Hospitalized Pediatric Patients

SUMMARY

Invasive fungal infections are a significant cause of morbidity and mortality in children with immunodeficiencies. Current dosing recommendations for voriconazole often result in subtherapeutic exposure in pediatric patients. In this single-center retrospective study, we reviewed hospitalized pediatric patients receiving voriconazole with at least one inpatient serum trough concentration measured. Patient characteristics and voriconazole dosing courses with associated trough concentrations were summarized for all patients as well as grouped by age (0 to 1 y, 2 to 11 y, and 12 to 18 y). Of 106 included patients, the median age was 9 years (range, 29 d to 18 y). Five hundred ninety courses of voriconazole were administered with 365 associated troughs. Most troughs were subtherapeutic (49%) and 30% of patients never attained a therapeutic trough. The median oral daily dose associated with a therapeutic trough was higher in younger age groups: 21.6 mg/kg 0 to 1 year, 17.9 mg/kg 2 to 11, and 9.5 mg/kg 12 to 18 years ( P <0.001). Patients younger than 2 years had the largest proportion of subtherapeutic troughs and variability in dosing. Attainment of therapeutic voriconazole concentrations was challenging across all pediatric age groups. Higher starting doses for patients younger than 2 years are likely needed.

Getting the dose right using physiologically-based pharmacokinetic modeling: dexamethasone to prevent post-extubation stridor in children as proof of concept

Introduction

Critically ill patients show large variability in drug disposition due to e.g., age, size, disease and treatment modalities. Physiologically-based pharmacokinetic (PBPK) models can be used to design individualized dosing regimens taking this into account. Dexamethasone, prescribed for the prevention post-extubation stridor (PES), is metabolized by the drug metabolizing enzyme CYP3A. As CYP3A4 undergoes major changes during childhood, we aimed to develop age-appropriate dosing recommendations for children of dexamethasone for PES, as proof of concept for PBPK modeling to individualize dosing for critically ill patients.

Methods

All simulations were conducted in Simcyp™ v21 (a population-based PBPK modeling platform), using an available dexamethasone compound model and pediatric population model in which CYP3A4 ontogeny is incorporated. Published pharmacokinetic (PK) data was used for model verification. Evidence for the dose to prevent post-extubation stridor was strongest for 2-6 year old children, hence simulated drug concentrations resulting from this dose from this age group were targeted when simulating age-appropriate doses for the whole pediatric age range.

Results

Dexamethasone plasma concentrations upon single and multiple intravenous administration were predicted adequately across the pediatric age range. Exposure-matched predictions of dexamethasone PK indicated that doses (in mg/kg) for the 2-6 years olds can be applied in 3 month-2 year old children, whereas lower doses are needed in children of other age groups (60% lower for 0-2 weeks, 40% lower for 2-4 weeks, 20% lower for 1-3 months, 20% lower for 6-12 year olds, 40% lower for 12-18 years olds).

Discussion

We show that PBPK modeling is a valuable tool that can be used to develop model-informed recommendations using dexamethasone to prevent PES in children. Based on exposure matching, the dose of dexamethasone should be reduced compared to commonly used doses, in infants <3 months and children ≥6 years, reflecting age-related variation in drug disposition. PBPK modeling is an promising tool to optimize dosing of critically ill patients.

Combining data on the bioavailability of midazolam and physiologically-based pharmacokinetic modeling to investigate intestinal CYP3A4 ontogeny

Pediatric physiologically-based modeling in drug development has grown in the past decade and optimizing the underlying systems parameters is important in relation to overall performance. In this study, variation of clinical oral bioavailability of midazolam as a function of age is used to assess the underlying ontogeny models for intestinal CYP3A4. Data on midazolam bioavailability in adults and children and different ontogeny patterns for intestinal CYP3A4 were first collected from the literature. A pediatric PBPK model was then used to assess six different ontogeny models in predicting bioavailability from preterm neonates to adults. The average fold error ranged from 0.7 to 1.38, with the rank order of least to most biased model being No Ontogeny < Upreti = Johnson < Goelen < Chen < Kiss. The absolute average fold error ranged from 1.17 to 1.64 with the rank order of most to least precise being Johnson > Upreti > No Ontogeny > Goelen > Kiss > Chen. The optimal ontogeny model is difficult to discern when considering the possible influence of CYP3A5 and other population variability; however, this study suggests that from term neonates and older a faster onset Johnson model with a lower fraction at birth may be close to this. For inclusion in other PBPK models, independent verification will be needed to confirm these results. Further research is needed in this area both in terms of age-related changes in midazolam and similar drug bioavailability and intestinal CYP3A4 ontogeny.

Clinical pharmacokinetics of leriglitazone and a translational approach using PBPK modeling to guide the selection of the starting dose in children

Leriglitazone is a unique peroxisome proliferator-activated receptor-gamma (PPARγ) agonist that crosses the blood-brain barrier in humans and clinical trials have shown evidence of efficacy in neurodegenerative diseases. At clinical doses which are well-tolerated, leriglitazone reaches the target central nervous system (CNS) concentrations that are needed for PPARγ engagement and efficacy; PPARγ engagement is also supported by clinical and anti-inflammatory biomarker changes in the Cerebrospinal fluid in the CNS. Plasma pharmacokinetics (PK) of leriglitazone were determined in a phase 1 study in male healthy volunteers comprising a single ascending dose (SAD) and a multiple ascending dose (MAD) at oral doses of 30, 90, and 270 mg and 135 and 270 mg, respectively. Leriglitazone was rapidly absorbed with no food effect on overall exposure and showed a linear PK profile with dose-exposure correlation. A physiologically based pharmacokinetic (PBPK) model was developed for leriglitazone based on phase 1 data (SAD part) and incorporated CYP3A4 (fmCYP3A4 = 24%) and CYP2C8-mediated (fmCYP2C8 = 45%) metabolism, as well as biliary clearance (feBIL = 19.5%) derived from in vitro data, and was verified by comparing the observed versus predicted concentration-time profiles from the MAD part. The PBPK model was prospectively applied to predict the starting pediatric doses and was preliminarily verified with data from five pediatric patients.

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.

Activity-Based Protein Profiling to Probe Relationships between Cytochrome P450 Enzymes and Early-Age Metabolism of Two Polycyclic Aromatic Hydrocarbons (PAHs): Phenanthrene and Retene

A growing body of literature has linked early-life exposures to polycyclic aromatic hydrocarbons (PAH) with adverse neurodevelopmental effects. Once in the body, metabolism serves as a powerful mediator of PAH toxicity by bioactivating and detoxifying PAH metabolites. Since enzyme expression and activity vary considerably throughout human development, we evaluated infant metabolism of PAHs as a potential contributing factor to PAH susceptibility. We measured and compared rates of phenanthrene and retene (two primary PAH constituents of woodsmoke) metabolism in human hepatic microsomes from individuals ≤21 months of age to a pooled sample (n = 200) consisting primarily of adults. We used activity-based protein profiling (ABPP) to characterize cytochrome P450 enzymes (CYPs) in the same hepatic microsome samples. Once incubated in microsomes, phenanthrene demonstrated rapid depletion. Best-fit models for phenanthrene metabolism demonstrated either 1 or 2 phases, depending on the sample, indicating that multiple enzymes could metabolize phenanthrene. We observed no statistically significant differences in phenanthrene metabolism as a function of age, although samples from the youngest individuals had the slowest phenanthrene metabolism rates. We observed slower rates of retene metabolism compared with phenanthrene also in multiple phases. Rates of retene metabolism increased in an age-dependent manner until adult (pooled) metabolism rates were achieved at ∼12 months. ABPP identified 28 unique CYPs among all samples, and we observed lower amounts of active CYPs in individuals ≤21 months of age compared to the pooled sample. Phenanthrene metabolism correlated to CYPs 1A1, 1A2, 2C8, 4A22, 3A4, and 3A43 and retene metabolism correlated to CYPs 1A1, 1A2, and 2C8 measured by ABPP and vendor-supplied substrate marker activities. These results will aid efforts to determine human health risk and susceptibility to PAHs exposure during early life.

A Mechanistic Physiologically-Based Pharmacokinetic Platform Model to Guide Adult and Pediatric Intravenous and Subcutaneous Dosing for Bispecific T Cell Engagers

Bispecific T cell engagers (Bi-TCEs) have revolutionized the treatment of oncology indications across both liquid and solid tumors. Bi-TCEs are rapidly evolving from conventional intravenous (i.v.) to more convenient subcutaneous (s.c.) administrations and extending beyond adults to also benefit pediatric patients. Leveraging clinical development experience across three generations of Bi-TCE molecules across both liquid and solid tumor indications from i.v./s.c. dosing in adults and pediatric subjects, we developed a mechanistic-physiologically-based pharmacokinetic (PBPK) platform model for Bi-TCEs. The model utilizes a full PBPK model framework and was successfully validated for PK predictions following i.v. and s.c. dosing across both liquid and solid tumor space in adults for eight Bi-TCEs. After refinement to incorporate physiological ontogeny, the model was successfully validated to predict pediatric PKs in 1 month - < 2 years, 2-11 years, and 12-17 years old subjects following i.v. dosing. Following s.c. dosing in pediatric subjects, the model predicted similar bioavailability, however, a shorter time to maximum concentration (Tmax ) for the three age groups compared with adults. The model was also applied to guide the dosing strategy for first generation of Bi-TCEs for organ impairment, specifically renal impairment, and was able to accurately predict the impact of renal impairment on PK for these relatively small-size Bi-TCEs. This work highlights a novel mechanistic platform model for accurately predicting the PK in adult and pediatric patients across liquid and solid tumor indications from i.v./s.c. dosing and can be used to guide optimal dose and dosing regimen selection and accelerating the clinical development for Bi-TCEs.

Challenges of pediatric pharmacotherapy: A narrative review of pharmacokinetics, pharmacodynamics, and pharmacogenetics

PURPOSE

Personalized pharmacotherapy, including for the pediatric population, provides optimal treatment and has emerged as a major trend owing to advanced drug therapeutics and diversified drug selection. However, it is essential to understand the growth and developmental characteristics of this population to provide appropriate drug therapy. In recent years, clinical pharmacogenetics has accumulated knowledge in pediatric pharmacotherapy, and guidelines from professional organizations, such as the Clinical Pharmacogenetics Implementation Consortium, can be consulted to determine the efficacy of specific drugs and the risk of adverse effects. However, the existence of a large knowledge gap hinders the use of these findings in clinical practice.

METHODS

We provide a narrative review of the knowledge gaps in pharmacokinetics (PK) and pharmacodynamics (PD) in the pediatric population, focusing on the differences from the perspective of growth and developmental characteristics. In addition, we explored PK/PD in relation to pediatric clinical pharmacogenetics.

RESULTS

The lack of direct and indirect biomarkers for more accurate assessment of the effects of drug administration limits the current knowledge of PD. In addition, incorporating pharmacogenetic insights as pivotal covariates is indispensable in this comprehensive synthesis for precision therapy; therefore, we have provided recommendations regarding the current status and challenges of personalized pediatric pharmacotherapy. The integration of clinical pharmacogenetics with the health care system and institution of educational programs for health care providers is necessary for its safe and effective implementation. A comprehensive understanding of the physiological and genetic complexities of the pediatric population will facilitate the development of effective and personalized pharmacotherapeutic strategies.

PBPK modeling to predict the pharmacokinetics of pantoprazole in different CYP2C19 genotypes

Pantoprazole is used to treat gastroesophageal reflux disease (GERD), maintain healing of erosive esophagitis (EE), and control symptoms related to Zollinger-Ellison syndrome (ZES). Pantoprazole is mainly metabolized by cytochrome P450 (CYP) 2C19, converting to 4'-demethyl pantoprazole. CYP2C19 is a genetically polymorphic enzyme, and the genetic polymorphism affects the pharmacokinetics and/or pharmacodynamics of pantoprazole. In this study, we aimed to establish the physiologically based pharmacokinetic (PBPK) model to predict the pharmacokinetics of pantoprazole in populations with various CYP2C19 metabolic activities. A comprehensive investigation of previous reports and drug databases was conducted to collect the clinical pharmacogenomic data, physicochemical data, and disposition properties of pantoprazole, and the collected data were used for model establishment. The model was evaluated by comparing the predicted plasma concentration-time profiles and/or pharmacokinetic parameters (AUC and Cmax) with the clinical observation results. The predicted plasma concentration-time profiles in different CYP2C19 phenotypes properly captured the observed profiles. All fold error values for AUC and Cmax were included in the two-fold range. Consequently, the minimal PBPK model for pantoprazole related to CYP2C19 genetic polymorphism was properly established and it can predict the pharmacokinetics of pantoprazole in different CYP2C19 phenotypes. The present model can broaden the insight into the individualized pharmacotherapy for pantoprazole.

Physiologically based mechanistic insight into differential risk of valproate hepatotoxicity between children and adults: A focus on ontogeny impact

The anticonvulsant valproic acid (VPA) despite complex pharmacokinetics has been in clinical use for nearly 6 decades. Previous reports indicated neonates, infants, and toddlers/preschoolers had higher risk of valproate hepatotoxicity than adults. However, dosing recommendations for those less than 10 years of age are lacking. To decipher clinical puzzles, physiologically-based pharmacokinetic (PBPK) models of VPA and its hepatotoxic metabolite 4-ene-VPA were constructed and simulated with particularly integrated information of drug-metabolizing enzyme ontogeny. Adult and pediatric PK data of VPA (n = 143 subjects) and 4-ene-VPA (n = 8 subjects) collected from previous reports were used for model development and validation. Sensitivity analyses were performed to characterize ontogeny impacts of CYP2C9 and UGT2B7 on dispositions of VPA and 4-ene-VPA across age groups. Optimal VPA dosing for each pediatric age group was also predicted and objectively judged by ensuring VPA efficacy and avoiding 4-ene-VPA hepatotoxicity. The study revealed UGT2B7 ontogeny was quite influential on VPA clearance even in neonates and small children. Intrinsic clearance of CYP2C9 was the most prominent determinant for areas under the concentration-time curve of VPA and 4-ene-VPA in infants, and toddlers/preschoolers, reflecting higher hepatotoxicity risk due to noxious 4-ene-VPA accumulation in these groups. The ontogeny-based PBPK approach complements conventional allometric methods in dosing estimation for the young by providing more mechanistic insight of the processes changing with age. The established ontogeny-based PBPK approach for VPA therapy deserves further corroboration by real-world therapeutic data to affirm its clinical applicability.

Physiologically Based Pharmacokinetics Modeling in the Neonatal Population-Current Advances, Challenges, and Opportunities

Physiologically based pharmacokinetic (PBPK) modeling is an approach to predicting drug pharmacokinetics, using knowledge of the human physiology involved and drug physiochemical properties. This approach is useful when predicting drug pharmacokinetics in under-studied populations, such as pediatrics. PBPK modeling is a particularly important tool for dose optimization for the neonatal population, given that clinical trials rarely include this patient population. However, important knowledge gaps exist for neonates, resulting in uncertainty with the model predictions. This review aims to outline the sources of variability that should be considered with developing a neonatal PBPK model, the data that are currently available for the neonatal ontogeny, and lastly to highlight the data gaps where further research would be needed.

Physiologically-Based Pharmacokinetic Modeling for Drug Dosing in Pediatric Patients: A Tutorial for a Pragmatic Approach in Clinical Care

It is well-accepted that off-label drug dosing recommendations for pediatric patients should be based on the best available evidence. However, the available traditional evidence is often low. To bridge this gap, physiologically-based pharmacokinetic (PBPK) modeling is a scientifically well-founded tool that can be used to enable model-informed dosing (MID) recommendations in children in clinical practice. In this tutorial, we provide a pragmatic, PBPK-based pediatric modeling workflow. For this approach to be successfully implemented in pediatric clinical practice, a thorough understanding of the model assumptions and limitations is required. More importantly, careful evaluation of an MID approach within the context of overall benefits and the potential risks is crucial. The tutorial is aimed to help modelers, researchers, and clinicians, to effectively use PBPK simulations to support pediatric drug dosing.

Escitalopram and Sertraline Population Pharmacokinetic Analysis in Pediatric Patients

BACKGROUND AND OBJECTIVE

Escitalopram and sertraline are commonly prescribed for anxiety and depressive disorders in children and adolescents. The pharmacokinetics (PK) of these medications have been evaluated in adults and demonstrate extensive variability, but studies in pediatric patients are limited. Therefore, we performed a population PK analysis for escitalopram and sertraline in children and adolescents to characterize the effects of demographic, clinical, and pharmacogenetic factors on drug exposure.

METHODS

A PK dataset was generated by extracting data from the electronic health record and opportunistic sampling of escitalopram- and sertraline-treated psychiatrically hospitalized pediatric patients aged 5-18 years. A population PK analysis of escitalopram and sertraline was performed using NONMEM. Concentration-time profiles were simulated using MwPharm++ to evaluate how covariates included in the final models influence medication exposure and compared to adult therapeutic reference ranges.

RESULTS

The final escitalopram cohort consisted of 315 samples from 288 patients, and the sertraline cohort consisted of 265 samples from 255 patients. A one-compartment model with a proportional residual error model best described the data for both medications. For escitalopram, CYP2C19 phenotype and concomitant CYP2C19 inhibitors affected apparent clearance (CL/F), and normalizing CL/F and apparent volume of distribution (V/F) to body surface area (BSA) improved estimations. The final escitalopram model estimated CL/F and V/F at 14.2 L/h/1.73 m2 and 428 L/1.73 m2, respectively. For sertraline, CYP2C19 phenotype and concomitant CYP2C19 inhibitors influenced CL/F, and empirical allometric scaling of patient body weight on CL/F and V/F was significant. The final sertraline model estimated CL/F and V/F at 124 L/h/70 kg and 4320 L/70 kg, respectively. Normalized trough concentrations (Ctrough) for CYP2C19 poor metabolizers taking escitalopram were 3.98-fold higher compared to normal metabolizers (151.1 ng/mL vs 38.0 ng/mL, p < 0.0001), and normalized Ctrough for CYP2C19 poor metabolizers taking sertraline were 3.23-fold higher compared to normal, rapid, and ultrarapid metabolizers combined (121.7 ng/mL vs 37.68 ng/mL, p < 0.0001). Escitalopram- and sertraline-treated poor metabolizers may benefit from a dose reduction of 50-75% and 25-50%, respectively, to normalize exposure to other phenotypes.

CONCLUSION

To our knowledge, this is the largest population PK analysis of escitalopram and sertraline in pediatric patients. Significant PK variability for both medications was observed and was largely explained by CYP2C19 phenotype. Slower CYP2C19 metabolizers taking escitalopram or sertraline may benefit from dose reductions given increased exposure.

Tacrolimus pharmacokinetics are influenced by CYP3A5, age, and concomitant fluconazole in pediatric kidney transplant patients

Tacrolimus, the most common immunosuppressant for organ transplant, has a narrow therapeutic range and is metabolized by CYP3A4/5. Trough concentration monitoring and dosing adjustments are used to reach a therapeutic range. CYP3A5 intermediate and normal metabolizers (*1 allele carriers; IM/NM) demonstrate faster tacrolimus metabolism than poor metabolizers (PM). We analyzed the electronic health records of 93 patients aged <21 years for the first 8 weeks after a kidney transplant between January 2010 and December 2021. The target tacrolimus trough was 10-15 ng/mL in the first 4 weeks and 7-10 ng/mL in the next 4 weeks. Banked DNA was collected and genotyped for CYP3A5*3, *6, *7, and *8 alleles. We found that CYP3A5 IM/NM (n = 21) took longer than PM (n = 72) to reach the therapeutic range (7 vs. 4 days, p = 0.048). IM/NM had more dose adjustments (8 vs. 6, p = 0.025) and needed >150% of the required daily dose compared with PM. The concentration/dose ratio was influenced by age and concomitant fluconazole (p = 0.0003, p = 0.034, respectively) and the average daily dose decreases with age in CYP3A5 PM (p = 0.001). Tremors were more common in patients who ever had a trough concentration >15 ng/mL compared with those who never had a trough concentration >15 ng/mL (OR 3.31, 95% CI 1.03-8.98, p = 0.038). Using standard dosing, CYP3A5 IM/NM took longer to reach the goal range and require more dose adjustments and higher doses than PM. Preemptive genotyping could decrease the number of dose changes necessary to reach a therapeutic dose. We have implemented pre-transplant CYP3A5 testing at our institution.

Development and Verification of a Japanese Pediatric Physiologically Based Pharmacokinetic Model with Emphasis on Drugs Eliminated by Cytochrome P450 or Renal Excretion

Physiologically based pharmacokinetic (PBPK) models are useful in bridging drug exposure in different ethnic groups, and there is increasing regulatory application of this approach in adults. Reported pediatric PBPK models tend to focus on the North European population, with few examples in other ethnic groups. This study describes the development and verification of a Japanese pediatric PBPK population. The development of the model was based on the existing North European pediatric population. Japanese systems and clinical data were collated from public databases and the literature, and the underlying demographics and equations were optimized so that physiological outputs represented the Japanese pediatric population. The model was tested using 14 different small molecule drugs, eliminated by a variety of pathways, including cytochrome P450 3A4 (CYP3A4) metabolism and renal excretion. Given the limitations of the clinical data, the overall performance of the model was good, with 44/62 predictions for PK parameters (area under the plasma drug concentration-time curve, AUC; maximum serum concentration, Cmax ; clearance, CL) being within 0.8- to 1.25-fold, 56/62 within 0.67- to 1.5-fold, and 61/62 within 0.5- to 2.0-fold of the observed values. Specific results for the 5 CYP3A4 substrates showed 20/31 cases were predicted within 0.8- to 1.25-fold, 27/31 within 0.67- to 1.5-fold, and all were within 0.5- to 2.0-fold of the observed values. Given the increased regulatory use of pediatric PBPK in drug development, expanding these models to other ethnic groups are important. Considering qualifying these models based on the context of use, there is a need to expand on the current research to include a larger range of drugs with different elimination pathways. Collaboration among academic, industry, model providers, and regulators will facilitate further development.

Investigating Tacrolimus Disposition in Paediatric Patients with a Physiologically Based Pharmacokinetic Model Incorporating CYP3A4 Ontogeny, Mechanistic Absorption and Red Blood Cell Binding

Tacrolimus is a crucial immunosuppressant for organ transplant patients, requiring therapeutic drug monitoring due to its variable exposure after oral intake. Physiologically based pharmacokinetic (PBPK) modelling has provided insights into tacrolimus disposition in adults but has limited application in paediatrics. This study investigated age dependency in tacrolimus exposure at the levels of absorption, metabolism, and distribution. Based on the literature data, a PBPK model was developed to predict tacrolimus exposure in adults after intravenous and oral administration. This model was then extrapolated to the paediatric population, using a unique reference dataset of kidney transplant patients. Selecting adequate ontogeny profiles for hepatic and intestinal CYP3A4 appeared critical to using the model in children. The best model performance was achieved by using the Upreti ontogeny in both the liver and intestines. To mechanistically evaluate the impact of absorption on tacrolimus exposure, biorelevant in vitro solubility and dissolution data were obtained. A relatively fast and complete release of tacrolimus from its amorphous formulation was observed when mimicking adult or paediatric dissolution conditions (dose, fluid volume). In both the adult and paediatric PBPK models, the in vitro dissolution profiles could be adequately substituted by diffusion-layer-based dissolution modelling. At the level of distribution, sensitivity analysis suggested that differences in blood plasma partitioning of tacrolimus may contribute to the variability in exposure in paediatric patients.

Integration of Biorelevant Pediatric Dissolution Methodology into PBPK Modeling to Predict In Vivo Performance and Bioequivalence of Generic Drugs in Pediatric Populations: a Carbamazepine Case Study

This study investigated the impact of gastro-intestinal fluid volume and bile salt (BS) concentration on the dissolution of carbamazepine (CBZ) immediate release (IR) 100 mg tablets and to integrate these in vitro biorelevant dissolution profiles into physiologically based pharmacokinetic modelling (PBPK) in pediatric and adult populations to determine the biopredictive dissolution profile. Dissolution profiles of CBZ IR tablets (100 mg) were generated in 50-900 mL biorelevant adult fasted state simulated gastric and intestinal fluid (Ad-FaSSGF and Ad-FaSSIF), also in three alternative compositions of biorelevant pediatric FaSSGF and FaSSIF medias at 200 mL. This study found that CBZ dissolution was poorly sensitive to changes in the composition of the biorelevant media, where dissimilar dissolution (F2 = 46.2) was only observed when the BS concentration was changed from 3000 to 89 μM (Ad-FaSSIF vs Ped-FaSSIF 50% 14 BS). PBPK modeling demonstrated the most predictive dissolution volume and media composition to forecast the PK was 500 mL of Ad-FaSSGF/Ad-FaSSIF media for adults and 200 mL Ped-FaSSGF/FaSSIF media for pediatrics. A virtual bioequivalence simulation was conducted by using Ad-FaSSGF and/or Ad-FaSSIF 500 mL or Ped-FaSSGF and/or Ped-FaSSIF 200 mL dissolution data for CBZ 100 mg (reference and generic test) IR product. The CBZ PBPK models showed bioequivalence of the product. This study demonstrates that the integration of biorelevant dissolution data can predict the PK profile of a poorly soluble drug in both populations. Further work using more pediatric drug products is needed to verify biorelevant dissolution data to predict the in vivo performance in pediatrics.

Estimation of FMO3 Ontogeny by Mechanistic Population Pharmacokinetic Modelling of Risdiplam and Its Impact on Drug-Drug Interactions in Children

BACKGROUND AND OBJECTIVE

Spinal muscular atrophy (SMA) is a progressive neuromuscular disease caused by insufficient levels of survival motor neuron (SMN) protein. Risdiplam (EvrysdiTM) increases SMN protein and is approved for the treatment of SMA. Risdiplam has high oral bioavailability and is primarily eliminated through hepatic metabolism by flavin-containing monooxygenase3 (FMO3) and cytochrome P450 (CYP) 3A, by 75% and 20%, respectively. While the FMO3 ontogeny is critical input data for the prediction of risdiplam pharmacokinetics (PK) in children, it was mostly studied in vitro, and robust in vivo FMO3 ontogeny is currently lacking. We derived in vivo FMO3 ontogeny by mechanistic population PK modelling of risdiplam and investigated its impact on drug-drug interactions in children.

METHODS

Population and physiologically based PK (PPK and PBPK) modelling conducted during the development of risdiplam were integrated into a mechanistic PPK (Mech-PPK) model to estimate in vivo FMO3 ontogeny. A total of 10,205 risdiplam plasma concentration-time data from 525 subjects aged 2 months-61 years were included. Six different structural models were examined to describe the in vivo FMO3 ontogeny. Impact of the newly estimated FMO3 ontogeny on predictions of drug-drug interaction (DDI) in children was investigated by simulations for dual CYP3A-FMO3 substrates including risdiplam and theoretical substrates covering a range of metabolic fractions (fm) of CYP3A and FMO3 (fmCYP3A:fmFMO3 = 10%:90%, 50%:50%, 90%:10%).

RESULTS

All six models consistently predicted higher FMO3 expression/activity in children, reaching a maximum at the age of 2 years with an approximately threefold difference compared with adults. Different trajectories of FMO3 ontogeny in infants < 4 months of age were predicted by the six models, likely due to limited observations for this age range. Use of this  in vivo FMO3 ontogeny function improved prediction of risdiplam PK in children compared to in vitro FMO3 ontogeny functions. The simulations of theoretical dual CYP3A-FMO3 substrates predicted comparable or decreased CYP3A-victim DDI propensity in children compared to adults across the range of fm values. Refinement of FMO3 ontogeny in the risdiplam model had no impact on the previously predicted low CYP3A-victim or -perpetrator DDI risk of risdiplam in children.

CONCLUSION

Mech-PPK modelling successfully estimated in vivo FMO3 ontogeny from risdiplam data collected from 525 subjects aged 2 months-61 years. To our knowledge, this is the first investigation of in vivo FMO3 ontogeny by population approach using comprehensive data covering a wide age range. Derivation of a robust in vivo FMO3 ontogeny function has significant implications on the prospective prediction of PK and DDI in children for other FMO3 substrates in the future, as illustrated in the current study for FMO3 and/or dual CYP3A-FMO3 substrates.

CLINICAL TRIAL REGISTRY NUMBERS

NCT02633709, NCT03032172, NCT02908685, NCT02913482, NCT03988907.

Minimal physiologically-based hybrid model of pharmacokinetics in pregnant women: Application to antenatal corticosteroids

Minimal physiologically-based pharmacokinetic (mPBPK) models are an alternative to full physiologically-based pharmacokinetic (PBPK) models as they offer reduced complexity while maintaining the physiological interpretation of key model components. Full PBPK models have been developed for pregnancy, but a mPBPK model eases the ability to perform a "top-down" meta-analysis melding all available pharmacokinetic (PK) data in the mother and fetus. Our hybrid mPBPK model consists of mPBPK models for the mother and fetus with connection by the placenta. This model was applied to describe the rich PK data of antenatal corticosteroid betamethasone (BET) jointly with the limited data for dexamethasone (DEX) in the mother and fetus. Physiologic model parameters were obtained from the literature while drug-dependent parameters were estimated by the simultaneous fitting of all available data for DEX and BET. Maternal clearances of DEX and BET confirmed the literature values, and the expected fetal-to-maternal plasma ratios ranged from 0.3 to 0.4 for both drugs. Simulations of maternal plasma concentrations for the dosing regimens of BET and DEX recommended by the World Health Organization based on our findings revealed up to 60% lower exposures than found in nonpregnant women and offers a means of devising alternative dosing regimens. Our hybrid mPBPK model and meta-analysis approach could facilitate assessment of other classes of drugs indicated for the treatment of pregnant women.

Precision sirolimus dosing in children: The potential for model-informed dosing and novel drug monitoring

The mTOR inhibitor sirolimus is prescribed to treat children with varying diseases, ranging from vascular anomalies to sporadic lymphangioleiomyomatosis to transplantation (solid organ or hematopoietic cell). Precision dosing of sirolimus using therapeutic drug monitoring (TDM) of sirolimus concentrations in whole blood drawn at the trough (before the next dose) time-point is the current standard of care. For sirolimus, trough concentrations are only modestly correlated with the area under the curve, with R 2 values ranging from 0.52 to 0.84. Thus, it should not be surprising, even with the use of sirolimus TDM, that patients treated with sirolimus have variable pharmacokinetics, toxicity, and effectiveness. Model-informed precision dosing (MIPD) will be beneficial and should be implemented. The data do not suggest dried blood spots point-of-care sampling of sirolimus concentrations for precision dosing of sirolimus. Future research on precision dosing of sirolimus should focus on pharmacogenomic and pharmacometabolomic tools to predict sirolimus pharmacokinetics and wearables for point-of-care quantitation and MIPD.

Another Step Toward Qualification of Pediatric Physiologically Based Pharmacokinetic Models to Facilitate Inclusivity and Diversity in Pediatric Clinical Studies

Robust prediction of pharmacokinetics (PKs) in pediatric subjects of diverse ages, ethnicities, and morbidities is critical. Qualification of pediatric physiologically-based pharmacokinetic (P-PBPK) models is an essential step toward enabling precision dosing of these vulnerable groups. Twenty-two manuscripts involving P-PBPK predictions and corresponding observed PK data (e.g., area under the curve and clearance) for 22 small-molecule compounds metabolized by CYP (3A4, 1A2, and 2C9), UGT (1A9 and 2B7), FMO3, renal, non-renal, and complex routes were identified; ratios of mean predicted/observed (P/O) PK parameters were calculated. Seventy-eight of 115 mean predicted PK parameters were within 0.8 to 1.25-fold of observed data, 98 within 0.67 to 1.5-fold, 109 within 2-fold, and only 6 P/O ratios were outside of these bounds. A set of 12 CYP3A4-metabolized compounds and a set of 6 metabolized by other enzymes, CYP1A2 (1 compound), CYP2C9 (2 compounds), UGT1A9 (1 compound) and UGT2B7 (2 compounds) had 56 of 59 and 22 of 25 mean P/O ratios, respectively, that fell within the > 0.5 and < 2.0-fold boundaries. For compounds covering renal, non-renal, complex, and FM03 routes of elimination, 29 of 31 mean P/O ratios fell within the 0.67 to 1.5-fold bounds, including 4 of 5 P/O ratios from newborns. P-PBPK modeling and simulation is a strategic component of the complement of precision dosing methods and has a vital role to play in dose adjustment in vulnerable pediatric populations, such as those with disease or in different ethnic groups. Qualification of such models is an essential step toward acceptance of this methodology by regulators.

FDA approval summary for lonafarnib (Zokinvy) for the treatment of Hutchinson-Gilford progeria syndrome and processing-deficient progeroid laminopathies

The U.S. Food and Drug Administration recently approved lonafarnib as the first treatment for Hutchinson-Gilford progeria syndrome (HGPS) and processing-deficient progeroid laminopathies. This approval was primarily based on a comparison of patients with HGPS treated with lonafarnib in 2 open-label trials with an untreated patient cohort. With up to 11 years of follow-up, it was found that the lonafarnib treated patients with HGPS had a survival benefit of 2.5 years compared with the untreated patients with HGPS. This large treatment effect on the objective endpoint of mortality using a well-matched comparator group mitigated potential sources of bias and together with other evidence, established compelling evidence of a drug effect with benefits that outweighed the risks. This approval is an example of U.S. Food and Drug Administration's regulatory flexibility for a rare disease while ensuring that standards for drug approval are met.

Sleep Quality and Duration in Children That Consume Caffeine: Impact of Dose and Genetic Variation in ADORA2A and CYP1A

Caffeine is the most consumed drug in the world, and it is commonly used by children. Despite being considered relatively safe, caffeine can have marked effects on sleep. Studies in adults suggest that genetic variants in the adenosine A2A receptor (ADORA2A, rs5751876) and cytochrome P450 1A (CYP1A, rs2472297, rs762551) loci are correlated with caffeine-associated sleep disturbances and caffeine intake (dose), but these associations have not been assessed in children. We examined the independent and interaction effects of daily caffeine dose and candidate variants in ADORA2A and CYP1A on the sleep quality and duration in 6112 children aged 9-10 years who used caffeine and were enrolled in the Adolescent Brain Cognitive Development (ABCD) study. We found that children with higher daily caffeine doses had lower odds of reporting > 9 h of sleep per night (OR = 0.81, 95% CI = 0.74-0.88, and p = 1.2 × 10^-6). For every mg/kg/day of caffeine consumed, there was a 19% (95% CI = 12-26%) decrease in the odds of children reporting > 9 h of sleep. However, neither ADORA2A nor CYP1A genetic variants were associated with sleep quality, duration, or caffeine dose. Likewise, genotype by caffeine dose interactions were not detected. Our findings suggest that a daily caffeine dose has a clear negative correlation with sleep duration in children, but this association is not moderated by the ADORA2A or CYP1A genetic variation.

Feasibility of a Pragmatic PBPK Modeling Approach: Towards Model-Informed Dosing in Pediatric Clinical Care

BACKGROUND AND OBJECTIVE

More than half of all drugs are still prescribed off-label to children. Pharmacokinetic (PK) data are needed to support off-label dosing, however for many drugs such data are either sparse or not representative. Physiologically-based pharmacokinetic (PBPK) models are increasingly used to study PK and guide dosing decisions. Building compound models to study PK requires expertise and is time-consuming. Therefore, in this paper, we studied the feasibility of predicting pediatric exposure by pragmatically combining existing compound models, developed e.g. for studies in adults, with a pediatric and preterm physiology model.

METHODS

Seven drugs, with various PK characteristics, were selected (meropenem, ceftazidime, azithromycin, propofol, midazolam, lorazepam, and caffeine) as a proof of concept. Simcyp® v20 was used to predict exposure in adults, children, and (pre)term neonates, by combining an existing compound model with relevant virtual physiology models. Predictive performance was evaluated by calculating the ratios of predicted-to-observed PK parameter values (0.5- to 2-fold acceptance range) and by visual predictive checks with prediction error values.

RESULTS

Overall, model predicted PK in infants, children and adolescents capture clinical data. Confidence in PBPK model performance was therefore considered high. Predictive performance tends to decrease when predicting PK in the (pre)term neonatal population.

CONCLUSION

Pragmatic PBPK modeling in pediatrics, based on compound models verified with adult data, is feasible. A thorough understanding of the model assumptions and limitations is required, before model-informed doses can be recommended for clinical use.

Effects of inflammation on voriconazole levels: a systematic review

AIMS

This study aimed to review the studies evaluating the effect of the inflammatory state on voriconazole (VRZ) levels.

METHODS

The study included randomized clinical trials, cohort studies, and case-control studies that focused on the influence of the inflammatory state on VRZ levels. Following the preferred reporting items for systematic reviews and meta-analyses (PRISMA) guidelines, relevant articles published until 2021 were searched in several databases, including PubMed, Embase, Web of Science, and the Cochrane Library.

RESULTS

Twenty studies were included in this review, of which 15 described adult populations, 3 described pediatric populations, and 2 included both adult and pediatric populations. Seventeen studies used C-reactive protein (CRP) as an indicator of inflammation, 6 described a dose-response relationship for the effect of inflammation represented by CRP on VRZ concentrations, and 4 examined the effect of CRP on the metabolic rate of VRZ.

CONCLUSIONS

Our findings showed that the level of inflammation can significantly affect VRZ levels. However, the effect of inflammation on VRZ concentrations in children is controversial and must be analyzed along with age. Clinicians dosing VRZ should take into account the patient's inflammatory state. The impact of inflammation on genotype-based dosing decisions requires further study to explain the high pharmacokinetic variability of VRZ.

A database of pediatric drug effects to evaluate ontogenic mechanisms from child growth and development

BACKGROUND

Adverse drug effects (ADEs) in children are common and may result in disability and death, necessitating post-marketing monitoring of their use. Evaluating drug safety is especially challenging in children due to the processes of growth and maturation, which can alter how children respond to treatment. Current drug safety-signal-detection methods do not account for these dynamics.

METHODS

We recently developed a method called disproportionality generalized additive models (dGAMs) to better identify safety signals for drugs across child-development stages.

FINDINGS

We used dGAMs on a database of 264,453 pediatric adverse-event reports and found 19,438 ADEs signals associated with development and validated these signals against a small reference set of pediatric ADEs. Using our approach, we can hypothesize on the ontogenic dynamics of ADE signals, such as that montelukast-induced psychiatric disorders appear most significant in the second year of life. Additionally, we integrated pediatric enzyme expression data and found that pharmacogenes with dynamic childhood expression, such as CYP2C18 and CYP27B1, are associated with pediatric ADEs.

CONCLUSIONS

We curated KidSIDES, a database of pediatric drug safety signals, for the research community and developed the Pediatric Drug Safety portal (PDSportal) to facilitate evaluation of drug safety signals across childhood growth and development.

FUNDING

This study was supported by grants from the National Institutes of Health (NIH).

Cytochrome P450 enzymes in the pediatric population: Connecting knowledge on P450 expression with pediatric pharmacokinetics

Cytochrome P450 enzymes play an important role in the pharmacokinetics, efficacy, and toxicity of drugs. Age-dependent changes in P450 enzyme expression have been studied based on several detection systems, as well as by deconvolution of in vivo pharmacokinetic data observed in pediatric populations. The age-dependent changes in P450 enzyme expression can be important determinants of drug disposition in childhood, in addition to the changes in body size and the other physiological parameters, and effects of pharmacogenetics and disease on organ functions. As a tool incorporating drug-specific and body-specific factors, physiologically-based pharmacokinetic (PBPK) models have become increasingly used to characterize and explore mechanistic insights into drug disposition. Thus, PBPK models can be a bridge between findings from basic science and utilization in predictive science. Pediatric PBPK models incorporate additional system specific information on developmental physiology and ontogeny and have been used to predict pharmacokinetic parameters from preterm neonates onwards. These models have been advocated by regulatory authorities in order to support pediatric clinical trials. The purpose of this chapter is to highlight accumulated knowledge and findings from basic research focusing on P450 enzymes, as well as the current status and future challenges of expanding the utilization of pediatric PBPK modeling.

Challenges in the Treatment of Invasive Aspergillosis in Immunocompromised Children

Invasive aspergillosis (IA) is associated with significant morbidity and mortality. Voriconazole remains the drug of choice for the treatment of IA in children; however, the complex kinetics of voriconazole in children make dosing challenging and therapeutic drug monitoring (TDM) essential for treatment success. The overarching goal of this review is to discuss the role of voriconazole, posaconazole, isavuconazole, liposomal amphotericin B, echinocandins, and combination antifungal therapy for the treatment of IA in children. We also provide a detailed discussion of antifungal TDM in children.

Development and application of a pediatric mechanistic kidney model

Pediatric physiologically‐based pharmacokinetic (P‐PBPK) models have been used to predict age related changes in the pharmacokinetics (PKs) of renally cleared drugs mainly in relation to changes in glomerular filtration rate. With emerging data on ontogeny of renal transporters, mechanistic models of renal clearance accounting for the role of active and passive secretion should be developed and evaluated. Data on age‐related physiological changes and ontogeny of renal transporters were applied into a mechanistic kidney within a P‐PBPK model. Plasma concentration–time profile and PK parameters of cimetidine, ciprofloxacin, metformin, tenofovir, and zidovudine were predicted in subjects aged 1 day to 18 years. The predicted and observed plasma concentration–time profiles and PK parameters were compared. The predicted concentration–time profile means and 5th and 95th percent intervals generally captured the observed data and variability in various studies. Overall, based on drugs and age bands, predicted to observed clearance were all within two‐fold and in 11 of 16 cases within 1.5‐fold. Predicted to observed area under the curve (AUC) and maximum plasma concentration (C_max) were within two‐fold in 12 of 14 and 12 of 15 cases, respectively. Predictions in neonates and early infants (up to 14 weeks postnatal age) were reasonable with 15–20 predicted PK parameters within two‐fold of the observed. ciprofloxacin but not zidovudine PK predictions were sensitive to basal kidney uptake transporter ontogeny. The results indicate that a mechanistic kidney model accounting for physiology and ontogeny of renal processes and transporters can predict the PK of renally excreted drugs in children. Further data especially in neonates are required to verify the model and ontogeny profiles.

The Development of a PBPK Model for Atomoxetine Using Levels in Plasma, Saliva and Brain Extracellular Fluid in Patients with Normal and Deteriorated Kidney Function

BACKGROUND

Atomoxetine is a treatment for attention-deficit hyperactivity disorder. It inhibits Norepinephrine Transporters (NET) in the brain. Renal impairment can reduce hepatic CYP2D6 activity and atomoxetine elimination which may increase its body exposure. Atomoxetine can be secreted in saliva.

OBJECTIVE

The objective of this work was to test the hypothesis that atomoxetine saliva levels (sATX) can be used to predict ATX brain Extracellular Fluid (bECF) levels and their pharmacological effects in healthy subjects and those with End-Stage Renal Disease (ESRD).

METHODS

The pharmacokinetics of atomoxetine after intravenous administration to rats with chemically induced acute and chronic renal impairments were investigated. A physiologically-based pharmacokinetic (PBPK) model was built and verified in rats using previously published measured atomoxetine levels in plasma and brain tissue. The rat PBPK model was then scaled to humans and verified using published measured atomoxetine levels in plasma, saliva, and bECF.

RESULTS

The rat PBPK model predicted the observed reduced atomoxetine clearance due to renal impairment in rats. The PBPK model predicted atomoxetine exposure in human plasma, sATX and bECF. Additionally, it predicted that ATX bECF levels needed to inhibit NET are achieved at 80 mg dose. In ESRD patients, the developed PBPK model predicted that the previously reported 65% increase in plasma exposure in these patients can be associated with a 63% increase in bECF. The PBPK simulations showed that there is a significant correlation between sATX and bECF in human.

CONCLUSION

Saliva levels can be used to predict atomoxetine pharmacological response.

Applications, Challenges, and Outlook for PBPK Modeling and Simulation: A Regulatory, Industrial and Academic Perspective

Several regulatory guidances on the use of physiologically based pharmacokinetic (PBPK) analyses and physiologically based biopharmaceutics model(s) (PBBM(s)) have been issued. Workshops are routinely held, demonstrating substantial interest in applying these modeling approaches to address scientific questions in drug development. PBPK models and PBBMs have remarkably contributed to model-informed drug development (MIDD) such as anticipating clinical PK outcomes affected by extrinsic and intrinsic factors in general and specific populations. In this review, we proposed practical considerations for a "base" PBPK model construction and development, summarized current status, challenges including model validation and gaps in system models, and future perspectives in PBPK evaluation to assess a) drug metabolizing enzyme(s)- or drug transporter(s)- mediated drug-drug interactions b) dosing regimen prediction, sampling timepoint selection and dose validation in pediatric patients from newborns to adolescents, c) drug exposure in patients with renal and/or and hepatic organ impairment, d) maternal-fetal drug disposition during pregnancy, and e) pH-mediated drug-drug interactions in patients treated with proton pump inhibitors/acid-reducing agents (PPIs/ARAs) intended for gastric protection. Since PBPK can simulate outcomes in clinical studies with enrollment challenges or ethical issues, the impact of PBPK models on waivers and how to strengthen study waiver is discussed.

Increasing application of pediatric physiologically based pharmacokinetic models across academic and industry organizations

There has been a significant increase in the use of physiologically based pharmacokinetic (PBPK) models during the past 20 years, especially for pediatrics. The aim of this study was to give a detailed overview of the growth and areas of application of pediatric PBPK (P‐PBPK) models. A total of 181 publications and publicly available regulatory reviews were identified and categorized according to year, author affiliation, platform, and primary application of the P‐PBPK model (in clinical settings, drug development or to advance pediatric model development in general). Secondary application areas, including dose selection, biologics, and drug interactions, were also assessed. The growth rate for P‐PBPK modeling increased 33‐fold between 2005 and 2020; this was mainly attributed to growth in clinical and drug development applications. For primary applications, 50% of articles were classified under clinical, 18% under drug development, and 33% under model development. The most common secondary applications were dose selection (75% drug development), pharmacokinetic prediction and covariate identification (47% clinical), and model parameter identification (68% model development), respectively. Although population PK modeling remains the mainstay of approaches supporting pediatric drug development, the data presented here demonstrate the widespread application of P‐PBPK models in both drug development and clinical settings. Although applications for pharmacokinetic and drug–drug interaction predictions in pediatrics is advocated, this approach remains underused in areas such as assessment of pediatric formulations, toxicology, and trial design. The increasing number of publications supporting the development and refinement of the pediatric model parameters can only serve to enhance optimal use of P‐PBPK models.

Physiologically based pharmacokinetic modeling of tadalafil to inform pediatric dose selection in children with pulmonary arterial hypertension

Tadalafil, a phosphodiesterase 5 inhibitor, is being investigated as a treatment for pulmonary arterial hypertension (PAH) in children aged 6 months to less than 18 years. Tadalafil pharmacokinetic (PK) data in children less than 2 years old are unavailable, therefore a physiologically based pharmacokinetic (PBPK) model was developed to enable estimation of tadalafil doses in children less than 2 years old. The model was verified in adults and extended for use in children by modifying CYP3A‐mediated intrinsic clearance to include CYP3A7. To account for co‐dosing of the commonly prescribed moderate CYP3A4 inducer bosentan, predicted exposures were increased by a factor of 1.54 based on changes in exposure in adults with PAH. This factor was predictable using a bosentan PBPK model. The tadalafil model was verified in children aged greater than or equal to 2 years by comparing predicted and observed exposures. Tadalafil doses for children less than 2 years old were calculated as target area under the concentration curve from zero to 24 h (AUC_0–24)/predicted AUC_0–24, with target AUC_0–24 of 10,000 ng*h/ml based on adult 40 mg single dose exposures determined in patients without bosentan background treatment. These doses were 2 mg, 3 mg, 4 mg, and 6 mg, respectively, for children aged birth to less than 1 month, 1 month to less than 6 months, 6 months to less than 1 year, and 1 to less than 2 years. Due to uncertainties in CYP maturation, a nonmechanistic steady‐state volume scalar, and lack of PK data in children less than 2 years old, accumulation of tadalafil to steady‐state in children less than 2 years was not verifiable. Safety of proposed doses is supported by postmarketing research and investigator‐led trials.

Physiologically-Based Pharmacokinetic Modeling to Investigate the Effect of Maturation on Buprenorphine Pharmacokinetics in Newborns with Neonatal Opioid Withdrawal Syndrome

Neonatal opioid withdrawal syndrome (NOWS) is a major public health concern whose incidence has paralleled the opioid epidemic in the United States. Sublingual buprenorphine is an emerging treatment for NOWS, but given concerns about long-term adverse effects of perinatal opioid exposure, precision dosing of buprenorphine is needed. Buprenorphine pharmacokinetics (PK) in newborns, however, is highly variable. To evaluate underlying sources of PK variability, a neonatal physiologically-based pharmacokinetic (PBPK) model of sublingual buprenorphine was developed using Simcyp (version 19.1). The PBPK model included metabolism by cytochrome P450 (CYP) 3A4, CYP2C8, UDP-glucuronosyltransferase (UGT) 1A1, UGT1A3, UGT2B7, and UGT2B17, with additional biliary excretion. Maturation of metabolizing enzymes was incorporated, and default CYP2C8 and UGT2B7 ontogeny profiles were updated according to recent literature. A biliary clearance developmental profile was outlined using clinical data from neonates receiving sublingual buprenorphine as NOWS treatment. Extensive PBPK model validation in adults demonstrated good predictability, with geometric mean (95% confidence interval (CI)) predicted/observed ratios (P/O ratios) of area under the curve from zero to infinity (AUC0-∞ ), peak concentration (Cmax ), and time to reach peak concentration (Tmax ) equaling 1.00 (0.74-1.33), 1.04 (0.84-1.29), and 0.95 (0.72-1.26), respectively. In neonates, the geometric mean (95% CI) P/O ratio of whole blood concentrations was 0.75 (95% CI 0.64-0.87). PBPK modeling and simulation demonstrated that variability in biliary clearance, sublingual absorption, and CYP3A4 abundance are likely important drivers of buprenorphine PK variability in neonates. The PBPK model could be used to guide development of improved buprenorphine starting dose regimens for the treatment of NOWS.

Performance Evaluation of Montelukast Pediatric Formulations: Part II - a PBPK Modelling Approach

This study aimed to build a physiologically based pharmacokinetic (PBPK) model coupled with age-appropriate in vitro dissolution data to describe drug performance in adults and pediatric patients. Montelukast sodium was chosen as a model drug. Two case studies were investigated: case study 1 focused on the description of formulation performance from adults to children; case study 2 focused on the description of the impact of medicine co-administration with vehicles on drug exposure in infants. The PBPK model for adults and pediatric patients was developed in Simcyp^® v18.2 informed by age-appropriate in vitro dissolution results obtained in a previous study. Oral administration of montelukast was simulated with the ADAM™ model. For case study 1, the developed PBPK model accurately described montelukast exposure in adults and children populations after the administration of montelukast chewable tablets. Two-stage dissolution testing in simulated fasted gastric to intestinal conditions resulted in the best description of in vivo drug performance in adults and children. For case study 2, a good description of in vivo drug performance in infants after medicine co-administration with vehicles was achieved by incorporating in vitro drug dissolution (under simulated fasted gastric to fed intestinal conditions) into a fed state PBPK model with consideration of the in vivo dosing conditions (mixing of formulation with applesauce or formula). The case studies presented demonstrate how a PBPK absorption modelling strategy can facilitate the description of drug performance in the pediatric population to support decision-making and biopharmaceutics understanding during pediatric drug development.

Population Pharmacokinetics of Ponatinib in Healthy Adult Volunteers and Patients With Hematologic Malignancies and Model-Informed Dose Selection for Pediatric Development

The BCR-ABL1 inhibitor ponatinib is approved for the treatment of adults with chronic myeloid leukemia or Philadelphia chromosome-positive acute lymphoblastic leukemia, including those with the T315I mutation. We report a population pharmacokinetic model-based analysis for ponatinib and its application to inform dose selection for pediatric development. Plasma concentration-time data were collected from 260 participants (86 healthy volunteers; 174 patients with hematologic malignancies) enrolled across 7 clinical trials. Data were analyzed using nonlinear mixed-effects modeling. Ponatinib pharmacokinetics were described by a 2-compartment model with first-order elimination from the central compartment. The final model included body weight and age as covariates on the apparent central volume of distribution; however, exposure variability explained by these covariates was small compared with overall variability in the population. None of the covariates evaluated, including sex, age (19-85 years), race, body weight (40.7-152.0 kg), total bilirubin (0.1-3.16 mg/dL), alanine aminotransferase (6-188 U/L), albumin (23.0-52.5 g/L), and creatinine clearance (≥28 mL/min) had clinically meaningful effects on apparent oral clearance. Simulations based on the final model predicted that daily doses of 15 to 45 mg result in steady-state average concentrations that are in the pharmacological range for BCR-ABL1 inhibition and approximate or exceed concentrations associated with suppression of T315I mutant clones. The final model was adapted using allometric scaling to inform dose selection for pediatric development. Clinicaltrials.gov identifier: NCT00660920; NCT01667133; NCT01650805.

Application of PBPK Modeling and Simulation for Regulatory Decision Making and Its Impact on US Prescribing Information: An Update on the 2018-2019 Submissions to the US FDA's Office of Clinical Pharmacology

Since 2016, results from physiologically based pharmacokinetic (PBPK) analyses have been routinely found in the clinical pharmacology section of regulatory applications submitted to the US Food and Drug Administration (FDA). In 2018, the Food and Drug Administration's Office of Clinical Pharmacology published a commentary summarizing the application of PBPK modeling in the submissions it received between 2008 and 2017 and its impact on prescribing information. In this commentary, we provide an update on the application of PBPK modeling in submissions received between 2018 and 2019 and highlight a few notable examples.

Leveraging Physiologically Based Pharmacokinetic Modeling and Experimental Data to Guide Dosing Modification of CYP3A-Mediated Drug-Drug Interactions in the Pediatric Population

Solithromycin is a novel fluoroketolide antibiotic that is both a substrate and time-dependent inhibitor of CYP3A. Solithromycin has demonstrated efficacy in adults with community-acquired bacterial pneumonia and has also been investigated in pediatric patients. The objective of this study was to develop a framework for leveraging physiologically based pharmacokinetic (PBPK) modeling to predict CYP3A-mediated drug-drug interaction (DDI) potential in the pediatric population using solithromycin as a case study. To account for age, we performed in vitro metabolism and time-dependent inhibition studies for solithromycin for CYP3A4, CYP3A5, and CYP3A7. The PBPK model included CYP3A4 and CYP3A5 metabolism and time-dependent inhibition, glomerular filtration, P-glycoprotein transport, and enterohepatic recirculation. The average fold error of simulated and observed plasma concentrations of solithromycin in both adults (1966 plasma samples) and pediatric patients from 4 days to 17.9 years (684 plasma samples) were within 0.5- to 2.0-fold. The geometric mean ratios for the simulated area under the concentration versus time curve (AUC) extrapolated to infinity were within 0.75- to 1.25-fold of observed values in healthy adults receiving solithromycin with midazolam or ketoconazole. DDI potential was simulated in pediatric patients (1 month to 17 years of age) and adults. Solithromycin increased the simulated midazolam AUC 4- to 6-fold, and ketoconazole increased the simulated solithromycin AUC 1- to 2-fold in virtual subjects ranging from 1 month to 65 years of age. This study presents a systematic approach for incorporating CYP3A in vitro data into adult and pediatric PBPK models to predict pediatric CYP3A-mediated DDI potential. SIGNIFICANCE STATEMENT: Using solithromycin, this study presents a framework for investigating and incorporating CYP3A4, CYP3A5, and CYP3A7 in vitro data into adult and pediatric physiologically based pharmacokinetic models to predict CYP3A-mediated DDI potential in adult and pediatric subjects during drug development. In this study, minor age-related differences in inhibitor concentration resulted in differences in the magnitude of the DDI. Therefore, age-related differences in DDI potential for substrates metabolized primarily by CYP3A4 can be minimized by closely matching adult and pediatric inhibitor concentrations.

Impact of Inflammation on Cytochromes P450 Activity in Pediatrics: A Systematic Review

Background and Objective

Cytochromes P450 (CYP) are the major enzymes involved in hepatic metabolism of drugs. Personalization of treatment in pediatrics is a major challenge, as it must not only take into account genetic, environmental, and physiological factors but also ontogeny. Published data in adults show that inflammation had an isoform-specific impact on CYP activities and we aimed to evaluate this impact in the pediatric population.

Methods

Articles listed in PubMed through 7 January, 2021 that studied the impact of inflammation on CYP activities in pediatrics were included in this systematic review. Sources of inflammation, victim drugs (CYP involved), effect of drug–disease interactions, number and age of subjects, and study design were extracted.

Results

Twenty-seven studies and case reports were included. The impact of inflammation on CYP activities appeared to be age dependent and isoform-specific, with some drug–disease interactions having significant pharmacokinetic and clinical impact. For example, midazolam clearance decreases by 70%, while immunosuppressant and theophylline concentrations increase three-fold and two-fold with intensive care unit admission and infection. Cytochrome P450 activity appears to return to baseline level when the disease is resolved.

Conclusions

Studies that have assessed the impact of inflammation on CYP activity are lacking in pediatrics, yet it is a major factor to consider to improve drug efficacy or safety. The scarce current data show that the impact of inflammation is isoform and age dependent. An effort must be made to improve the understanding of the impact of inflammation on CYP activities in children to better individualize treatment.

Ontogeny of Drug-Metabolizing Enzymes

Almost 50% of prescription drugs lack age-appropriate dosing guidelines and therefore are used "off-label." Only ~10% drugs prescribed to neonates and infants have been studied for safety or efficacy. Immaturity of drug metabolism in children is often associated with drug toxicity. This chapter summarizes data on the ontogeny of major human metabolizing enzymes involved in oxidation, reduction, hydrolysis, and conjugation of drugs. The ontogeny data of individual drug-metabolizing enzymes are important for accurate prediction of drug pharmacokinetics and toxicity in children. This information is critical for designing clinical studies to appropriately test pharmacological hypotheses and develop safer pediatric drugs, and to replace the long-standing practice of body weight- or surface area-normalized drug dosing. The application of ontogeny data in physiologically based pharmacokinetic model and regulatory submission are discussed.

CYP2C19 Phenotype and Body Weight-Guided Voriconazole Initial Dose in Infants and Children after Hematopoietic Cell Transplantation

Prophylactic voriconazole use is recommended for children undergoing hematopoietic cell transplantation (HCT). Dosing considerations are essential, due to the narrow therapeutic window of voriconazole. Known covariates do not sufficiently explain the large interindividual pharmacokinetic (PK) variability of voriconazole. Moreover, knowledge of voriconazole PK for age <2 years is limited. We investigated genetic and clinical covariate associations with voriconazole interindividual PK variability and subsequently simulated dosing regimens in children. This study was conducted as part of a single-institution, phase I study of intravenous voriconazole therapy for children undergoing HCT. We conducted a population PK analysis and tested covariate effects on voriconazole PK, including 67 genetic variants and clinical variables. We analyzed plasma voriconazole and N-oxide metabolite concentrations from 58 children <21 years of age (including 12 children <2 years of age). A two-compartment parent mixed linear/nonlinear model best described our data. The CYP2C19 phenotype and body weight were significant covariates (P < 0.05 for both). Our model performance for age <2 years was comparable to that for other age groups. Simulation of the final model suggested the following doses to attain target steady-state trough concentrations of 1.5 to 5.0 mg/liter for the CYP2C19 normal phenotype: 16 mg/kg (weight of <15 kg), 12 mg/kg (weight of 15 to 30 kg), or 10 mg/kg (weight of >30 kg); doses were 33 to 50% lower for CYP2C19 poor/intermediate phenotypes and 25 to 50% higher for CYP2C19 rapid/ultrarapid phenotypes. We propose a new starting-dose regimen, combined with therapeutic drug monitoring, for intravenous voriconazole therapy in children of all ages. Future studies should validate this dosing regimen.

Chlorpyrifos and Δ9 Tetrahydrocannabinol exposure and effects on parameters associated with the endocannabinoid system and risk factors for obesity

Marilyn Silva. Retired from a career in toxicology and risk assessment. Increased childhood and adult obesity are associated with chlorpyrifos (CPF), an organophosphate pesticide. Cannabis (Δ9Tetrahydrocannabinol: Δ9THC) use has increased globally with legalization. CPF applications on cannabis crops lacks federally regulated tolerances and may pose health risks through exposure during development and in adulthood. Both CPF and Δ9THC affect the endocannabinoid system (eCBS), a regulator of appetite, energy balance, and gut microbiota, which, if disrupted, increases risk for obesity and related diseases. CPF inhibits eCB metabolism and Δ9THC is a partial agonist/antagonist at the cannabinoid receptor (CB1R). Effects of each on obesogenic parameters were examined via literature search. Male rodents with CPF exposure showed increased body weights, dysbiosis, inflammation and oxidative stress, potentially associated with increased eCBs acting through the gut-microbiota-adipose-brain regulatory loop. Δ9THC generally decreased body weights via partial agonism at the CB1R, lowering levels of eCBs. Dysbiosis and/or oxidative stress associated inflammation occurred with CPF, but these parameters were not tested with Δ9THC. Database deficiencies included limited endpoints to compare between chemicals/age-groups, inter-study variables (dose ranges, dosing vehicle, rodent strain, treatment duration, etc.). CPF and Δ9THC were not tested together, but human co-chemical effects would depend on exposure ratio, subject age, exposure duration, and health status, among others. An overriding concern is that both chemicals are well-documented developmental neurotoxins in addition to their low dose effects on energy balance. A co-exposure risk assessment is warranted with increased use and lack of federal CPF regulation on cannabis.

Model-Based Drug-Drug Interaction Extrapolation Strategy From Adults to Children: Risdiplam in Pediatric Patients With Spinal Muscular Atrophy

Risdiplam (Evrysdi) improves motor neuron function in patients with spinal muscular atrophy (SMA) and has been approved for the treatment of patients ≥2 months old. Risdiplam exhibits time‐dependent inhibition of cytochrome P450 (CYP) 3A in vitro. While many pediatric patients receive risdiplam, a drug–drug interaction (DDI) study in pediatric patients with SMA was not feasible. Therefore, a novel physiologically‐based pharmacokinetic (PBPK) model‐based strategy was proposed to extrapolate DDI risk from healthy adults to children with SMA in an iterative manner. A clinical DDI study was performed in healthy adults at relevant risdiplam exposures observed in children. Risdiplam caused an 1.11‐fold increase in the ratio of midazolam area under the curve with and without risdiplam (AUCR)), suggesting an 18‐fold lower in vivo CYP3A inactivation constant compared with the in vitro value. A pediatric PBPK model for risdiplam was validated with independent data and combined with a validated midazolam pediatric PBPK model to extrapolate DDI from adults to pediatric patients with SMA. The impact of selected intestinal and hepatic CYP3A ontogenies on the DDI susceptibility in children relative to adults was investigated. The PBPK analysis suggests that primary CYP3A inhibition by risdiplam occurs in the intestine rather than the liver. The PBPK‐predicted risdiplam CYP3A inhibition risk in pediatric patients with SMA aged 2 months–18 years was negligible (midazolam AUCR of 1.09–1.18) and included in the US prescribing information of risdiplam. Comprehensive evaluation of the sensitivity of predicted CYP3A DDI on selected intestinal and hepatic CYP3A ontogeny functions, together with PBPK model‐based strategy proposed here, aim to guide and facilitate DDI extrapolations in pediatric populations.

Physiologically Based Pharmacokinetic Models Are Effective Support for Pediatric Drug Development

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

Physiologically Based Pharmacokinetic Modeling and Allometric Scaling in Pediatric Drug Development: Where Do We Draw the Line?

Developing medicines for children is now established in legislation in both the United States and Europe; new drugs require pediatric study or investigation plans as part of their development. Particularly in early age groups, many developmental processes are not reflected by simple scalars such as body weight or body surface area, and even projecting doses based on simple allometric scaling can lead to significant overdoses in certain age groups. Modeling and simulation methodology, including physiologically based modeling, has evolved as part of the drug development toolkit and is being increasingly applied to various aspects of pediatric drug development. Pediatric physiologically based pharmacokinetic (PBPK) models account for the development of organs and the ontogeny of specific enzymes and transporters that determine the age-related pharmacokinetic profiles. However, when should this approach be used, and when will simpler methods such as allometric scaling suffice in answering specific problems? The aim of this review article is to illustrate the application of allometric scaling and PBPK in pediatric drug development and explore the optimal application of the latter approach with reference to case examples. In reality, allometric scaling included as part of population pharmacokinetic and PBPK approaches are all part of a model-informed drug development toolkit helping with decision making during the process of drug discovery and development; to that end, they should be viewed as complementary.