A Preterm Physiologically Based Pharmacokinetic Model. Part I: Physiological Parameters and Model Building

Salivary Therapeutic Monitoring of Buprenorphine in Neonates After Maternal Sublingual Dosing Guided by Physiologically Based Pharmacokinetic Modeling

BACKGROUND

Opioid use disorder (OUD) during pregnancy is associated with high mortality rates and neonatal opioid withdrawal syndrome (NOWS). Buprenorphine, an opioid, is used to treat OUD and NOWS. Buprenorphine active metabolite (norbuprenorphine) can cross the placenta and cause neonatal respiratory depression (EC 50 = 35 ng/mL) at high brain extracellular fluid (bECF) levels. Neonatal therapeutic drug monitoring using saliva decreases the likelihood of distress and infections associated with frequent blood sampling.

METHODS

An adult physiologically based pharmacokinetic model for buprenorphine and norbuprenorphine after intravenous and sublingual administration was constructed, vetted, and scaled to newborn and pregnant populations. The pregnancy model predicted that buprenorphine and norbuprenorphine doses would be transplacentally transferred to the newborns. The newborn physiologically based pharmacokinetic model was used to estimate the buprenorphine and norbuprenorphine levels in newborn plasma, bECF, and saliva after these doses.

RESULTS

After maternal sublingual administration of buprenorphine (4 mg/d), the estimated plasma concentrations of buprenorphine and norbuprenorphine in newborns exceeded the toxicity thresholds for 8 and 24 hours, respectively. However, the norbuprenorphine bECF levels were lower than the respiratory depression threshold. Furthermore, the salivary buprenorphine threshold levels in newborns for buprenorphine analgesia, norbuprenorphine analgesia, and norbuprenorphine hypoventilation were observed to be 22, 2, and 162 ng/mL.

CONCLUSIONS

Using neonatal saliva for buprenorphine therapeutic drug monitoring can facilitate newborn safety during the maternal treatment of OUD using sublingual buprenorphine. Nevertheless, the suitability of using adult values of respiratory depression EC 50 for newborns must be confirmed.

Model-informed drug development in pediatric, pregnancy and geriatric drug development: States of the art and future

The challenges of drug development in pediatric, pregnant and geriatric populations are a worldwide concern shared by regulatory authorities, pharmaceutical companies, and healthcare professionals. Model-informed drug development (MIDD) can integrate and quantify real-world data of physiology, pharmacology, and disease processes by using modeling and simulation techniques to facilitate decision-making in drug development. In this article, we reviewed current MIDD policy updates, reflected on the integrity of physiological data used for MIDD and the effects of physiological changes on the drug PK, as well as summarized current MIDD strategies and applications, so as to present the state of the art of MIDD in pediatric, pregnant and geriatric populations. Some considerations are put forth for the future improvements of MIDD including refining regulatory considerations, improving the integrity of physiological data, applying the emerging technologies, and exploring the application of MIDD in new therapies like gene therapies for special populations.

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.

Interdisciplinary Collaboration on Real World Data to Close the Knowledge Gap: A Reflection on "De Sutter et al. Predicting Volume of Distribution in Neonates: Performance of Physiologically Based Pharmacokinetic Modelling"

This commentary further reflects on the paper of De Sutter et al. on predicting volume of distribution in neonates, and the performance of physiologically based pharmacokinetic models We hereby stressed the add on value to collaborate on real world data to further close this knowledge gap. We illustrated this by weight distribution characteristics in breastfed (physiology) and in asphyxiated (pathophysiology), with additional reflection on their kidney and liver function.

Pragmatic physiologically-based pharmacokinetic modeling to support clinical implementation of optimized gentamicin dosing in term neonates and infants: proof-of-concept

Introduction

Modeling and simulation can support dosing recommendations for clinical practice, but a simple framework is missing. In this proof-of-concept study, we aimed to develop neonatal and infant gentamicin dosing guidelines, supported by a pragmatic physiologically-based pharmacokinetic (PBPK) modeling approach and a decision framework for implementation.

Methods

An already existing PBPK model was verified with data of 87 adults, 485 children and 912 neonates, based on visual predictive checks and predicted-to-observed pharmacokinetic (PK) parameter ratios. After acceptance of the model, dosages now recommended by the Dutch Pediatric Formulary (DPF) were simulated, along with several alternative dosing scenarios, aiming for recommended peak (i.e., 8-12 mg/L for neonates and 15-20 mg/L for infants) and trough (i.e., <1 mg/L) levels. We then used a decision framework to weigh benefits and risks for implementation.

Results

The PBPK model adequately described gentamicin PK. Simulations of current DPF dosages showed that the dosing interval for term neonates up to 6 weeks of age should be extended to 36-48 h to reach trough levels <1 mg/L. For infants, a 7.5 mg/kg/24 h dose will reach adequate peak levels. The benefits of these dose adaptations outweigh remaining uncertainties which can be minimized by routine drug monitoring.

Conclusion

We used a PBPK model to show that current DPF dosages for gentamicin in term neonates and infants needed to be optimized. In the context of potential uncertainties, the risk-benefit analysis proved positive; the model-informed dose is ready for clinical implementation.

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.

Predicting Volume of Distribution in Neonates: Performance of Physiologically Based Pharmacokinetic Modelling

The volume of distribution at steady state (Vss) in neonates is still often estimated through isometric scaling from adult values, disregarding developmental changes beyond body weight. This study aimed to compare the accuracy of two physiologically based pharmacokinetic (PBPK) Vss prediction methods in neonates (Poulin & Theil with Berezhkovskiy correction (P&T+) and Rodgers & Rowland (R&R)) with isometrical scaling. PBPK models were developed for 24 drugs using in-vitro and in-silico data. Simulations were done in Simcyp (V22) using predefined populations. Clinical data from 86 studies in neonates (including preterms) were used for comparison, and accuracy was assessed using (absolute) average fold errors ((A)AFEs). Isometric scaling resulted in underestimated Vss values in neonates (AFE: 0.61), and both PBPK methods reduced the magnitude of underprediction (AFE: 0.82-0.83). The P&T+ method demonstrated superior overall accuracy compared to isometric scaling (AAFE of 1.68 and 1.77, respectively), while the R&R method exhibited lower overall accuracy (AAFE: 2.03). Drug characteristics (LogP and ionization type) and inclusion of preterm neonates did not significantly impact the magnitude of error associated with isometric scaling or PBPK modeling. These results highlight both the limitations and the applicability of PBPK methods for the prediction of Vss in the absence of clinical data.

Evaluating the Pharmacokinetics of Fentanyl in the Brain Extracellular Fluid, Saliva, Urine, and Plasma of Newborns from Transplacental Exposure from Parturient Mothers Dosed with Epidural Fentanyl Utilizing PBPK Modeling

BACKGROUND AND OBJECTIVE

Fentanyl can mitigate the mother and newborn complications resulting from labor pain. However, fentanyl shows a narrow therapeutic index between its respiratory depressive and analgesic effects. Thus, prenatally acquired high fentanyl levels in the newborn brain extracellular fluid (bECF) may induce respiratory depression which requires therapeutic drug monitoring (TDM). TDM using saliva and urine in newborns can reduce the possibility of infections and distress associated with TDM using blood. The objective of this study was to develop  a physiologically based pharmacokinetic (PBPK) model to predict fentanyl concentrations in different newborn tissues due to intrauterine exposure.

METHODS

A fentanyl PBPK model in adults after intravenous and epidural administration was built, validated, and scaled to pregnancy and newborn populations. The dose that the newborn received transplacentally at birth was calculated using the pregnancy model. Then, the newborn bECF, saliva, plasma, and urine concentrations after such a dose were predicted using the newborn PBPK model.

RESULTS

After a maternal epidural dose of fentanyl 245 µg, the predicted newborn plasma and bECF levels were below the toxicity thresholds. Furthermore, the salivary threshold levels in newborns for fentanyl analgesic and respiratory depression effects were estimated to be 0.39 and 14.7-18.2 ng/ml, respectively.

CONCLUSION

The salivary TDM of fentanyl in newborns can be useful in newborns exposed to intrauterine exposure from parturient females dosed with epidural fentanyl. However, newborn-specific values of µ-opioid receptors IC50 for respiratory depression are needed.

Infant Exposure to Antituberculosis Drugs via Breast Milk and Assessment of Potential Adverse Effects in Breastfed Infants: Critical Review of Data

Infants of mothers treated for tuberculosis might be exposed to drugs via breast milk. The existing information on the exposure of breastfed infants lacks a critical review of the published data. We aimed to evaluate the quality of the existing data on antituberculosis (anti-TB) drug concentrations in the plasma and milk as a methodologically sound basis for the potential risk of breastfeeding under therapy. We performed a systematic search in PubMed for bedaquiline, clofazimine, cycloserine/terizidone, levofloxacin, linezolid, pretomanid/pa824, pyrazinamide, streptomycin, ethambutol, rifampicin and isoniazid, supplemented with update references found in LactMed^®. We calculated the external infant exposure (EID) for each drug and compared it with the recommended WHO dose for infants (relative external infant dose) and assessed their potential to elicit adverse effects in the breastfed infant. Breast milk concentration data were mainly not satisfactory to properly estimate the EID. Most of the studies suffer from limitations in the sample collection, quantity, timing and study design. Infant plasma concentrations are extremely scarce and very little data exist documenting the clinical outcome in exposed infants. Concerns for potential adverse effects in breastfed infants could be ruled out for bedaquiline, cycloserine/terizidone, linezolid and pyrazinamide. Adequate studies should be performed covering the scenario in treated mothers, breast milk and infants.

The Analysis of Pethidine Pharmacokinetics in Newborn Saliva, Plasma, and Brain Extracellular Fluid After Prenatal Intrauterine Exposure from Pregnant Mothers Receiving Intramuscular Dose Using PBPK Modeling

BACKGROUND AND OBJECTIVE

Pethidine (meperidine) can decrease labor pain-associated mother's hyperventilation and high cortisol-induced newborn complications. However, prenatal transplacentally acquired pethidine can cause side effects in newborns. High pethidine concentrations in the newborn brain extracellular fluid (bECF) can cause a serotonin crisis. Therapeutic drug monitoring (TDM) in newborns' blood distresses them and increases infection incidence, which can be overcome by using salivary TDM. Physiologically based pharmacokinetic (PBPK) modeling can predict drug concentrations in newborn plasma, saliva, and bECF after intrauterine pethidine exposure.

METHODS

A healthy adult PBPK model was constructed, verified, and scaled to newborn and pregnant populations after intravenous and intramuscular pethidine administration. The pregnancy PBPK model was used to predict the newborn dose received transplacentally at birth, which was used as input to the newborn PBPK model to predict newborn plasma, saliva, and bECF pethidine concentrations and set correlation equations between them.

RESULTS

Pethidine can be classified as a Salivary Excretion Classification System class II drug. The developed PBPK model predicted that, after maternal pethidine intramuscular doses of 100 mg and 150 mg, the newborn plasma and bECF concentrations were below the toxicity thresholds. Moreover, it was estimated that newborn saliva concentrations of 4.7 µM, 11.4 µM, and 57.7 µM can be used as salivary threshold concentrations for pethidine analgesic effects, side effects, and the risk for serotonin crisis, respectively, in newborns.

CONCLUSION

It was shown that saliva can be used for pethidine TDM in newborns during the first few days after delivery to mothers receiving pethidine.

Physiologically Based Pharmacokinetic (PBPK) Model-Informed Dosing Guidelines for Pediatric Clinical Care: A Pragmatic Approach for a Special Population

Physiologically based pharmacokinetic (PBPK) modeling can be an attractive tool to increase the evidence base of pediatric drug dosing recommendations by making optimal use of existing pharmacokinetic (PK) data. A pragmatic approach of combining available compound models with a virtual pediatric physiology model can be a rational solution to predict PK and hence support dosing guidelines for children in real-life clinical care, when it can also be employed by individuals with little experience in PBPK modeling. This comes within reach as user-friendly PBPK modeling platforms exist and, for many drugs and populations, models are ready for use. We have identified a list of drugs that can serve as a starting point for pragmatic PBPK modeling to address current clinical dosing needs.

General clinical and methodological considerations on the extrapolation of pharmacokinetics and optimization of study protocols for small molecules and monoclonal antibodies in children

Pharmacometric modelling plays a key role in both the design and analysis of regulatory trials in paediatric drug development. Studies in adults provide a rich source of data to inform the paediatric investigation plans, including knowledge on drug pharmacokinetics (PK), safety and efficacy. In children, drug disposition differs widely from birth to adolescence but extrapolating adult to paediatric PK, safety and efficacy either with pharmacometric or physiologically based approaches can help design or in some cases reduce the need for clinical studies. Aspects to consider when extrapolating PK include the maturation of drug metabolizing enzyme expression, glomerular filtration, drug excretory systems, and the expression and activity of specific transporters in conjunction with other drug properties such as fraction unbound. Knowledge of these can be used to develop extrapolation tools such as allometric scaling plus maturation functions or physiologically based PK. PK/pharmacodynamic approaches and well-designed clinical trials in children are of key importance in paediatric drug development. In this white paper, state-of-the-art of current methods used for paediatric extrapolation will be discussed. This paper is part of a conect4children implementation of innovative methodologies including pharmacometric and physiologically based PK modelling in clinical trial design/paediatric drug development through dissemination of expertise and expert advice. The suggestions arising from this white paper should define a minimum set of standards in paediatric modelling and contribute to the regulatory science.

Population pharmacokinetic modeling of caffeine in preterm infants with apnea of prematurity: New findings from concomitant erythromycin and AHR genetic polymorphisms

Current standard-dose caffeine therapy results in significant intersubject variability. The aims of this study were to develop and evaluate population pharmacokinetic (PPK) models of caffeine in preterm infants through comprehensive screening of covariates and then to propose model-informed precision dosing of caffeine for this population. A total of 129 caffeine concentrations from 96 premature neonates were incorporated into this study. Comprehensive medical record and genotype data of these neonates were collected for analysis. PPK modeling was performed by a nonlinear mixed effects modeling program (NONMEM). Final models based on the current weight (CW) or body surface area (BSA) were evaluated via multiple graphic and statistical methods. The model-informed dosing regimen was performed through Monte Carlo simulations. In addition to CW or BSA, postnatal age, coadministration with erythromycin (ERY), and aryl hydrocarbon receptor coding gene (AHR) variant (rs2158041) were incorporated into the final PPK models. Multiple evaluation results showed satisfactory prediction performance and stability of the CW- and BSA-based models. Monte Carlo simulations demonstrated that trough concentrations of caffeine in preterm infants would be affected by concomitant ERY therapy and rs2158041 under varying dose regimens. For the first time, ERY and rs2158041 were found to be associated with the clearance of caffeine in premature infants. Similar predictive performance and stability were obtained for both CW- and BSA-based PPK models. These findings provide novel insights into caffeine precision therapy for preterm infants.

Current and future physiologically based pharmacokinetic (PBPK) modeling approaches to optimize pharmacotherapy in preterm neonates

INTRODUCTION

There is a need for structured approaches to inform on pharmacotherapy in preterm neonates. With their proven track record up to regulatory acceptance, physiologically based pharmacokinetic (PBPK) modeling and simulation provide such a structured approach, and hold the promise to support drug development in preterm neonates.

AREAS COVERED

Compared to the general and pediatric use of PBPK modeling, its use to inform pharmacotherapy in preterms is limited. Using a systematic search (PBPK + preterm), we retained 25 records (20 research papers, 2 letters, 3 abstracts). We subsequently collated the published information on PBPK software packages (PK-Sim®, Simcyp®), and their applications and optimization efforts in preterm neonates. It is encouraging that these applications cover a broad range of scenarios (pharmacokinetic-dynamic analyses, drug-drug interactions, developmental pharmacogenetics, lactation related exposure) and compounds (small molecules, proteins). Furthermore, specific compartments (cerebrospinal fluid, tissue) or (patho)physiologic processes (cardiac output, biliary excretion, first pass metabolism) are considered.

EXPERT OPINION

Knowledge gaps exist, giving rise to various levels of model uncertainty in PBPK applications in preterm neonates. To improve this setting, we need cross talk between clinicians and modelers to generate and integrate knowledge (PK datasets, system knowledge, maturational physiology and pathophysiology) to further refine PBPK models.

Optimal Teicoplanin Dosing Regimen in Neonates and Children Developed by Leveraging Real-World Clinical Information

BACKGROUND

Teicoplanin is a glycopeptide antibiotic used for the treatment of methicillin-resistant Staphylococcus aureus infections. To ensure successful target attainment, therapeutic drug monitoring-informed dosage adjustment is recommended. However, it relies on the experience of the clinician and the frequency of drug measurements. This study aimed to design a new optimal dosing regimen of teicoplanin with a maintenance dosing strategy for neonates and children based on their physiological characteristics.

METHODS

Data from teicoplanin-treated patients (n = 214) were collected from electronic medical records. Covariate analyses were performed using population pharmacokinetic (PK) modeling with 399 serum teicoplanin concentrations from 48 neonates and 166 children. Multiple PK simulations were conducted to explore optimal dosing regimens that would allow control of the trough concentration to the target of 15-30 mg/L quicker than the current standard regimen.

RESULTS

Allometrically scaled body weight, postmenstrual age (PMA), renal function, and serum albumin were implemented as substantial covariates for teicoplanin clearance in a two-compartment PK model. Covariate analyses and comprehensive simulation assessments recommended the following modifications to the current regimen: (1) decreased dose for premature babies (PMA ≤28 weeks), (2) decreased dose for children with renal dysfunction, and (3) increased dose for children (0.5-11 years) with an estimated glomerular filtration rate of ≥90 mL/min/1.73 m2.

CONCLUSIONS

This study leverages real-world clinical information and proposes new optimal dosing regimens for teicoplanin in neonates and children through PK modeling and simulation analyses, taking into account the age, including PMA, and renal function of patients.

From immature pharmacotherapy towards pharmacotherapy of the immature

To truly attain effective and safe pharmacotherapy, the similarities and dissimilarities in physiology between micro-preemies and extreme preterm infants should be explored. The higher incidence of pulmonary hypertension and presence of adrenal insufficiency of prematurity in micro-preemies hereby serve as illustrations. The current limited data on pharmacokinetics, -dynamics and safety reflect the obvious need to collect such data, and to tailor modelling tools to their physiology and needs. Drug utilization hereby mirrors different needs and practices and may serve to guide prioritization decisions. Physiological data, combined with even limited observations on pharmacokinetics and -dynamics can be translated to effective modelling tools to attain effective and safe pharmacotherapy. We therefore discuss how valid research tools in pharmacology like physiology-based pharmacokinetic models can be developed, and how clinicians can contribute to such efforts, with the overarching aim to enable this shift from immature pharmacotherapy to pharmacotherapy for the immature.

Prediction of glomerular filtration rate maturation across preterm and term neonates and young infants using inulin as marker

Describing glomerular filtration rate (GFR) maturation across the heterogeneous population of preterm and term neonates and infants is important to predict the clearance of renally cleared drugs. This study aims to describe the GFR maturation in (pre)term neonates and young infants (PNA < 90 days) using individual inulin clearance data (CL_inulin). To this end, published GFR maturation models were evaluated by comparing their predicted GFR with CL_inulin retrieved from literature. The best model was subsequently optimized in NONMEM V7.4.3 to better fit the CL_inulin values. Our study evaluated seven models and collected 381 individual CL_inulin values from 333 subjects with median (range) birthweight (BWb) 1880 g (580-4950), gestational age (GA) 34 weeks (25-43), current weight (CW) 1890 g (480-6200), postnatal age (PNA) 3 days (0-75), and CL_inulin 2.20 ml/min (0.43-17.90). The De Cock 2014 model (covariates: BWb and PNA) performed the best in predicting CL_inulin, followed by the Rhodin 2009 model (covariates: CW and postmenstrual age). The final optimized model shows that GFR at birth is determined by BWb, thereafter the maturation rate of GFR is dependent on PNA and GA, with a higher GA showing an overall faster maturation. To conclude, using individual CL_inulin data, we found that a model for neonatal GFR requires a distinction between prenatal maturation quantified by BWb and postnatal maturation. To capture postnatal GFR maturation in (pre)term neonates and young infants, we developed an optimized model in which PNA-related maturation was dependent on GA.

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

INTRODUCTION

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

AREAS COVERED

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

EXPERT OPINION

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

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.

Application of a Physiologically Based Pharmacokinetic Approach to Predict Theophylline Pharmacokinetics Using Virtual Non-Pregnant, Pregnant, Fetal, Breast-Feeding, and Neonatal Populations

Perinatal pharmacology is influenced by a myriad of physiological variables that are changing dynamically. The influence of these covariates has not been assessed systemically. The objective of this work was to use theophylline as a model drug and to predict its pharmacokinetics before, during (including prediction of the umbilical cord level), and after pregnancy as well as in milk (after single and multiple doses) and in neonates using a physiological-based pharmacokinetic (PBPK) model. Neonatal theophylline exposure from milk consumption was projected in both normal term and preterm subjects. Predicted infant daily doses were calculated using theophylline average and maximum concentration in the milk as well as an estimate of milk consumption. Predicted concentrations and parameters from the PBPK model were compared to the observed data. PBPK predicted theophylline concentrations in non-pregnant and pregnant populations at different gestational weeks were within 2-fold of the observations and the observed concentrations fell within the 5th-95th prediction interval from the PBPK simulations. The PBPK model predicted an average cord-to-maternal plasma ratio of 1.0, which also agrees well with experimental observations. Predicted postpartum theophylline concentration profiles in milk were also in good agreement with observations with a predicted milk-to-plasma ratio of 0.68. For an infant of 2 kg consuming 150 ml of milk per day, the lactation model predicted a relative infant dose (RID) of 12 and 17% using predicted average (C_avg,ss) and maximum (C_max,ss) concentration in milk at steady state. The maximum RID of 17% corresponds to an absolute infant daily dose of 1.4 ± 0.5 mg/kg/day. This dose, when administered as 0.233 mg/kg every 4 h, to resemble breastfeeding frequency, resulted in plasma concentrations as high as 3.9 (1.9-6.8) mg/L and 2.8 (1.3-5.3) (5th-95th percentiles) on day 7 in preterm (32 GW) and full-term neonatal populations.

Prediction of cytochromes P450 3A and 2C19 modulation by both inflammation and drug interactions using physiologically based pharmacokinetics

Xenobiotics can interact with cytochromes P450 (CYPs), resulting in drug-drug interactions, but CYPs can also contribute to drug-disease interactions, especially in the case of inflammation, which downregulates CYP activities through pretranscriptional and posttranscriptional mechanisms. Interleukin-6 (IL-6), a key proinflammatory cytokine, is mainly responsible for this effect. The aim of our study was to develop a physiologically based pharmacokinetic (PBPK) model to foresee the impact of elevated IL-6 levels in combination with drug interactions with esomeprazole on CYP3A and CYP2C19. Data from a cohort of elective hip surgery patients whose CYP3A and CYP2C19 activities were measured before and after surgery were used to validate the accurate prediction of the developed models. Successive steps were to fit models for IL-6, esomeprazole, and omeprazole and its metabolite from the literature and to validate them. The models for midazolam and its metabolite were obtained from the literature. When appropriate, a correction factor was applied to convert drug concentrations from whole blood to plasma. Mean ratios between simulated and observed areas under the curve for omeprazole/5-hydroxy omeprazole, esomeprazole, and IL-6 were 1.53, 1.06, and 0.69, respectively, indicating an accurate prediction of the developed models. The impact of IL-6 and esomeprazole on the exposure to CYP3A and CYP2C19 probe substrates and respective metabolites were correctly predicted. Indeed, the ratio between predicted and observed mean concentrations were <2 for all observations (ranging from 0.51 to 1.7). The impact of IL-6 and esomeprazole on CYP3A and CYP2C19 activities after a hip surgery were correctly predicted with the developed PBPK models.

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.

Physiologically based pharmacokinetic modelling of acetaminophen in preterm neonates-The impact of metabolising enzyme ontogeny and reduced cardiac output

In preterm neonates, physiologically based pharmacokinetic (PBPK) models are suited for studying the effects of maturational and non-maturational factors on the pharmacokinetics of drugs with complex age-dependent metabolic pathways like acetaminophen (APAP). The aim of this study was to determine the impact of drug metabolising enzymes ontogeny on the pharmacokinetics of APAP in preterm neonates and to study the effect of reduced cardiac output (CO) on its PK using PBPK modelling. A PBPK model for APAP was first developed and validated in adults and then scaled to paediatric age groups to account for the effect of enzyme ontogeny. In preterm neonates, CO was reduced by 10%, 20%, and 30% to determine how this might affect APAP PK in preterm neonates. In all age groups, the predicted concentration-time profiles of APAP were within 5th and 95th percentile of the clinically observed concentration-time profiles and the predicted Cmax and AUC were within 2-folds of the reported parameters in clinical studies. Sulfation accounted for most of APAP metabolism in children, with the highest contribution of 68% in preterm neonates. A reduction in CO by up to 30% did not significantly alter the clearance of APAP in preterm neonates. The model successfully incorporated the ontogeny of drug metabolising enzymes involved in APAP metabolism and adequately predicted the PK of APAP in preterm neonates. A reduction in hepatic perfusion as a result of up to 30% reduction in CO has no effect on the PK of APAP in preterm neonates.

Progress in Drug Development-Pediatric Dose Selection: Workshop Summary

The "Pediatric Dose Selection" workshop was held in October 2020 and sponsored by the U.S. Food and Drug Administration and the University of Maryland Center for Excellence in Regulatory Science and Innovation. A summary of the presentations in the context of pediatric drug development is summarized in this article.

Considerations for Drug Dosing in Premature Infants

In premature infants, effective and safe drug therapy depends on optimal dose selection and requires a thorough understanding of the underlying disease(s) of these fragile infants as well as the pharmacokinetics and pharmacodynamics of the drugs selected to treat their diseases. Differences in gestational and postnatal age or weight are the major determinants of the observed variability in drug disposition and effect in these infants. This article presents an outline on how to translate the results of a population pharmacokinetic/pharmacodynamic study into rational dosing regimens, and how physiologically based pharmacokinetic modeling, electronic health records, and the abundantly available data of vital functions of premature infants during their stay in the neonatal intensive care unit for evaluation of their pharmacotherapy can be used to tailor the most safe and effective dose in these infants.

Preterm Physiologically Based Pharmacokinetic Model. Part II: Applications of the Model to Predict Drug Pharmacokinetics in the Preterm Population

BACKGROUND

Preterm neonates are usually not part of a traditional drug development programme, however they are frequently administered medicines. Developing modelling and simulation tools, such as physiologically based pharmacokinetic (PBPK) models that incorporate developmental physiology and maturation of drug metabolism, can be used to predict drug exposure in this group of patients, and may help to optimize drug dose adjustment.

OBJECTIVE

The aim of this study was to assess and verify the predictability of a preterm PBPK model using compounds that undergo diverse renal and/or hepatic clearance based on the knowledge of their disposition in adults.

METHODS

A PBPK model was developed in the Simcyp Simulator V17 to predict the pharmacokinetics (PK) of drugs in preterm neonates. Drug parameters for alfentanil, midazolam, caffeine, ibuprofen, gentamicin and vancomycin were collated from the literature. Predicted PK parameters and profiles were compared against the observed data.

RESULTS

The preterm PBPK model predicted the PK changes of the six compounds using ontogeny functions for cytochrome P450 (CYP) 1A2, CYP2C9 and CYP3A4 after oral and intravenous administrations. For gentamicin and vancomycin, the maturation of renal function was able to predict the exposure of these two compounds after intravenous administration. All PK parameter predictions were within a twofold error criteria.

CONCLUSION

While the developed preterm model for the prediction of PK behaviour in preterm patients is not intended to replace clinical studies, it can potentially help with deciding on first-time dosing in this population and study design in the absence of clinical data.

Development of a Virtual Chinese Pediatric Population Physiological Model Targeting Specific Metabolism and Kidney Elimination Pathways

Background: Physiologically based pharmacokinetic (PBPK) modeling and simulating may be a powerful tool in predicting drug behaviors in specific populations. It is a mathematical model that relates the pharmacokinetic (PK) profile of a compound with human anatomical characteristics, physiological characteristics, and biochemical parameters. Predictions using PBPK models offer a promising way to guide drug development and can be used to optimize clinical dosing regimens. However, PK data of new drugs in the pediatric population are too limited to guide clinical therapy, which may lead to frequent adverse events or insufficient efficacy for pediatric patients, particularly in neonates and infants.

Objective: The objective of this study was to establish a virtual Chinese pediatric population based on the physiological parameters of Chinese children that could be utilized in PBPK models.

Methods: A Chinese pediatric PBPK model was developed in Simcyp Simulator by collecting published Chinese pediatric physiological and anthropometric data to use as system parameters. This pediatric population model was then evaluated in the Chinese pediatric population by predicting the pharmacokinetic characteristics of four probe drugs: theophylline (major CYP1A2 substrate), fentanyl (major CYP3A4 substrate), vancomycin, and ceftazidime (renal-eliminated).

Results: The predicted maximum concentration (C_max), area under the curve of concentration-time (AUC), and clearance (CL) for theophylline (CYP1A2 metabolism pathway) and fentanyl (CYP3A4 metabolism pathway) were within two folds of the observed data. For drugs mainly eliminated by renal clearance (vancomycin and ceftazidime) in the Chinese pediatric population, the ratio of prediction to observation for major PK parameters was within a 2-fold error range.

Conclusion: The model is a supplement to the previous Chinese population PBPK model. We anticipate the model to be a better representative of the pediatric Chinese population for drugs PK, offering greater clinical precision for medication given to the pediatric population, ultimately advancing clinical development of pediatric drugs. We can refine this model further by collecting more physiological parameters of Chinese children.

Physiologically based pharmacokinetic/pharmacodynamic model for the prediction of morphine brain disposition and analgesia in adults and children

Morphine is a widely used opioid analgesic, which shows large differences in clinical response in children, even when aiming for equivalent plasma drug concentrations. Age-dependent brain disposition of morphine could contribute to this variability, as developmental increase in blood-brain barrier (BBB) P-glycoprotein (Pgp) expression has been reported. In addition, age-related pharmacodynamics might also explain the variability in effect. To assess the influence of these processes on morphine effectiveness, a multi-compartment brain physiologically based pharmacokinetic/pharmacodynamic (PB-PK/PD) model was developed in R (Version 3.6.2). Active Pgp-mediated morphine transport was measured in MDCKII-Pgp cells grown on transwell filters and translated by an in vitro-in vivo extrapolation approach, which included developmental Pgp expression. Passive BBB permeability of morphine and its active metabolite morphine-6-glucuronide (M6G) and their pharmacodynamic parameters were derived from experiments reported in literature. Model simulations after single dose morphine were compared with measured and published concentrations of morphine and M6G in plasma, brain extracellular fluid (ECF) and cerebrospinal fluid (CSF), as well as published drug responses in children (1 day– 16 years) and adults. Visual predictive checks indicated acceptable overlays between simulated and measured morphine and M6G concentration-time profiles and prediction errors were between 1 and -1. Incorporation of active Pgp-mediated BBB transport into the PB-PK/PD model resulted in a 1.3-fold reduced brain exposure in adults, indicating only a modest contribution on brain disposition. Analgesic effect-time profiles could be described reasonably well for older children and adults, but were largely underpredicted for neonates. In summary, an age-appropriate morphine PB-PK/PD model was developed for the prediction of brain pharmacokinetics and analgesic effects. In the neonatal population, pharmacodynamic characteristics, but not brain drug disposition, appear to be altered compared to adults and older children, which may explain the reported differences in analgesic effect.

Developmental processes in children can affect pharmacokinetics: “what the body does to the drug” as well as pharmacodynamics: “what the drug does to the body”. A typical example is morphine, of which the analgesic response is variable and particularly neonates suffer more often from respiratory depression, even when receiving doses corrected for differences in elimination. One way to mathematically incorporate developmental processes is by employing physiologically based pharmacokinetic/pharmacodynamic (PB-PK/PD) models, where physiological differences between individuals are incorporated. In this study, we developed a morphine PB-PK/PD model to predict brain drug disposition as well as analgesic response in adults and children, as both processes could potentially contribute to developmental variability in the effect of morphine. We found that age-related variation in BBB expression of the main morphine efflux transporter P-glycoprotein was not responsible for differences in brain exposure. In contrast, pharmacodynamic modelling suggested an increased sensitivity to morphine in neonates.

Does "Birth" as an Event Impact Maturation Trajectory of Renal Clearance via Glomerular Filtration? Reexamining Data in Preterm and Full-Term Neonates by Avoiding the Creatinine Bias

Glomerular filtration rate (GFR) is an important measure of renal function. Various models for its maturation have recently been compared; however, these have used markers, which are subject to different renal elimination processes. Inulin clearance data (a purer probe of GFR) collected from the literature were used to determine age‐related changes in GFR aspects of renal drug excretion in pediatrics. An ontogeny model was derived using a best‐fit model with various combinations of covariates such as postnatal age, gestational age at birth, and body weight. The model was applied to the prediction of systemic clearance of amikacin, gentamicin, vancomycin, and gadobutrol. During neonatal life, GFR increased as a function of both gestational age at birth and postnatal age, hence implying an impact of birth and a discrepancy in GFR for neonates with the same postmenstrual age depending on gestational age at birth (ie, neonates who were outside the womb longer had higher GFR, on average). The difference in GFR between pre‐term and full‐term neonates with the same postmenstrual age was negligible from beyond 1.25 years. Considering both postnatal age and gestational age at birth in GFR ontogeny models is important because postmenstrual age alone ignores the impact of birth. Most GFR models use covariates of body size in addition to age. Therefore, prediction from these models will also depend on the change in anthropometric characteristics with age. The latter may not be similar in various ethnic groups, and this makes the head‐to‐head comparison of models very challenging.

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

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