Kinetic Modeling for BT200 to Predict the Level of Plasma-Derived Coagulation Factor VIII in Humans.

Lack of Factor VIII (FVIII) concentrates is one of limiting factors for Hemophilia A prophylaxis in resource-limited countries. Rondaptivon pegol (BT200) is a pegylated aptamer and has been shown to elevate the level of von Willebrand Factor (VWF) and FVIII in previous studies. A population pharmacokinetic model for BT200 was built and linked to the kinetic models of VWF and FVIII based on reasonable assumptions. The developed PK/PD model for BT200 described the observed kinetic of BT200, VWF, and FVIII in healthy volunteers and patients with mild-to-moderate hemophilia A from two clinical trials. The developed model was evaluated using an external dataset in patients with severe hemophilia A taking recombinant FVIII products. The developed and evaluated PK/PD model was able to describe and predict concentration-time profiles of BT200, VWF, and FVIII in healthy volunteers and patients with hemophilia A. Concentration-time profiles of FVIII were then predicted following coadministration of plasma-derived FVIII concentrate and BT200 under various dosing scenarios in virtual patients with severe hemophilia A. Plasma-derived products, that contain VWF, are more accessible in low-resource countries as compared to their recombinant counterparts. The predicted time above 1 and 3 IU/dL FVIII in one week was compared between scenarios in the absence and presence of BT200. A combination dose of 6 mg BT200 once weekly plus 10 IU/kg plasma-derived FVIII twice weekly maintained similar coverage to a 30 IU/kg FVIII thrice weekly dose in absence of BT200, representing only 22% of the FVIII dose per week.

Estimated prophylactic dose required to achieve 3% trough as a function of age and concentrate class in multi-country severe WAPPS-Hemo haemophilia patients.

INTRODUCTION: Web-Accessible Population-Pharmacokinetic Service-Haemophilia (WAPPS-Hemo) data are available to study factor-concentrate usage, defined as the required weekly dose to achieve a 3% trough (WD3T), across standard and extended half-life (SHL/EHL) products. AIM: To provide baseline usage data including (i) differences across plasma-derived (pdSHL) versus recombinant (rSHL) products, (ii) SHL versus EHL, and (iii) effect of age and positive inhibitor history. METHODS: PK profiles (n = 14,416 patients, 0.3-85.2 years) and linear mixed effects models were used to estimate usage versus age, controlling for significant factors, using 95% confidence intervals to perform comparisons across all ages and posthoc tests to assess the differences. RESULTS: Average usage was significantly higher for pdSHL versus rSHL in patients with a positive inhibitor history (PIH; 1.9-2.5 times higher), for SHL versus EHL (4-10 times), and was significantly associated with age. CONCLUSION: Baseline usage patterns from 2017 to early 2023 provide a benchmark for assessing the impact of emerging technologies in haemophilia.

Health utilities in adults with hemophilia A: A retrospective cohort study.

INTRODUCTION: Haemophilia A negatively affects a patient's quality of life. There is a limited amount of health utility data (a measure of health-related quality of life) available for patients with haemophilia A. This information is crucial for cost-effectiveness analysis for haemophilia A treatment. OBJECTIVES: The goal of this project is to elicit the health utilities and factors impacting utility values for haemophilia A patients in Canada. METHODS: This is a population-based, cross-sectional, retrospective study of health utilities in patients with haemophilia A using Patient Report Outcomes Burdens and Experiences (PROBE) components from the Canadian Bleeding Disorders Registry (CBDR). A review of the mean utilities for three severity states, defined by clotting factor VIII level, was completed. A multiple linear regression analysis was completed to examine the determinants of health utilities including age, treatment type, chronic pain status, number of limited joints, and bleed rate. RESULTS: The average utility values (and standard deviations) for patients with haemophilia A in Canada are .79(.17), .76(.20), and .77(.19) for patients with severe, moderate, and mild haemophilia. The regression showed chronic pain status and the number of additional comorbidities as major significant factors (p-value < .001) in haemophilia A utility. Haemophilia severity was shown to be a major factor with smaller p-value (p-value < .05). CONCLUSIONS: Haemophilia A patients have lower utility than the general population. Chronic pain was shown to be a significant, major factor in health-related quality of life. Our study is essential for valuing health outcomes in haemophilia A-related cost-effectiveness analysis.

Canadian clinical experience on switching from standard half-life recombinant factor VIII (rFVIII), octocog alfa, to extended half-life rFVIII, damoctocog alfa pegol, in persons with haemophilia A ≥ 12 years followed in a Comprehensive Hemophilia Care Program in Canada.

INTRODUCTION: Damoctocog alfa pegol (BAY 94-9027, Jivi®) is an extended half-life recombinant factor (F)VIII replacement, indicated for the treatment of haemophilia A in patients aged ≥12 years. Following introduction of damoctocog alfa pegol in Canada in 2020, there have been no reports on routine clinical effectiveness and satisfaction, when switching from a previous FVIII product in Canada. AIM: To report changes in pharmacokinetics, effectiveness, utilization and patient satisfaction when switching to damoctocog alfa pegol prophylaxis from previous standard half-life octocog alfa (BAY 81-8973, Kovaltry®) treatment. METHODS: A single-centre, intra-patient comparison of pharmacokinetics and clinical outcomes was performed. Blood samples drawn once pre-dose and ≥2 times post-dose were measured by a one-stage assay to assess pharmacokinetic parameters including area under the curve (AUC, primary endpoint). Patient-reported outcomes data were collected using the Patient-Reported Outcomes, Burdens and Experiences questionnaire (PROBE). Clinical outcomes included annualized bleeding rate (ABR) and factor utilization. RESULTS: Dose-normalized AUC was significantly increased after switch to damoctocog alfa pegol from octocog alfa. Median (quartile [Q]1; Q3) annualized bleeding rates were 0.67 (0.00; 1.33) with damoctocog alfa pegol and 1.33 (0.00; 2.67) with octocog alfa. Half of the patients receiving damoctocog alfa pegol prophylaxis experienced zero bleeds (n = 9, 50.0%) versus 38.9% (n = 7) of patients treated with octocog alfa. Patients' good quality of life was maintained. CONCLUSION: This study provides routine clinical evidence supporting the benefits of switching from octocog alfa to damoctocog alfa pegol for patients with severe haemophilia A.

Maternal Ezetimibe Concentrations Measured in Breast Milk and Its Use in Breastfeeding Infant Exposure Predictions.

BACKGROUND: Lactating mothers taking ezetimibe, an antihyperlipidemic agent, may be hesitant to breastfeed despite the known benefit of breastfeeding to both mother and infant. Currently, no data exist on the presence or concentration of ezetimibe and its main active metabolite, ezetimibe-glucuronide (EZE-glucuronide), in human breast milk. METHODS: Voluntary breast milk samples containing ezetimibe and EZE-glucuronide were attained from lactating mothers taking ezetimibe as part of their treatment. An assay was developed and validated to measure ezetimibe and EZE-glucuronide concentrations in breast milk. A workflow that utilized a developed and evaluated pediatric physiologically based pharmacokinetic (PBPK) model, the measured concentrations in milk, and weight-normalized breast milk intake volumes was applied to predict infant exposures and determine the upper area under the curve ratio (UAR). RESULTS: Fifteen breast milk samples from two maternal-infant pairs were collected. The developed liquid chromatography-tandem mass spectrometry (LC-MS/MS) assay showed an analytical range of 0.039-5.0 ng/mL and 0.39-50.0 ng/mL for ezetimibe and EZE-glucuronide, respectively. The measured concentrations in the breast milk samples were 0.17-1.02 ng/mL and 0.42-2.65 ng/mL of ezetimibe and EZE-glucuronide, respectively. The evaluated pediatric PBPK model demonstrated minimal exposure overlap in adult therapeutic dose and breastfed infant simulated area under the concentration-time curve from time zero to 24 h (AUC24). Calculated UAR across infant age groups ranged from 0.0015 to 0.0026. CONCLUSIONS: PBPK model-predicted ezetimibe and EZE-glucuronide exposures and UAR suggest that breastfeeding infants would receive non-therapeutic exposures. Future work should involve a 'mother-infant pair study' to ascertain breastfed infant plasma ezetimibe and EZE-glucuronide concentrations to confirm the findings of this work.

Maximum likelihood estimation of renal transporter ontogeny profiles for pediatric PBPK modeling.

Optimal treatment of infants with many renally cleared drugs must account for maturational differences in renal transporter (RT) activity. Pediatric physiologically-based pharmacokinetic (PBPK) models may incorporate RT activity, but this requires ontogeny profiles for RT activity in children, especially neonates, to predict drug disposition. Therefore, RT expression measurements from human kidney postmortem cortical tissue samples were normalized to represent a fraction of mature RT activity. Using these data, maximum likelihood estimated the distributions of RT activity across the pediatric age spectrum, including preterm and term neonates. PBPK models of four RT substrates (acyclovir, ciprofloxacin, furosemide, and meropenem) were evaluated with and without ontogeny profiles using average fold error (AFE), absolute average fold error (AAFE), and proportion of observations within the 5-95% prediction interval. Novel maximum likelihood profiles estimated ontogeny distributions for the following RT: OAT1, OAT3, OCT2, P-gp, URAT1, BCRP, MATE1, MRP2, MRP4, and MATE-2 K. Profiles for OAT3, P-gp, and MATE1 improved infant furosemide and neonate meropenem PBPK model AFE from 0.08 to 0.70 and 0.53 to 1.34 and model AAFE from 12.08 to 1.44 and 2.09 to 1.36, respectively, and improved the percent of data within the 5-95% prediction interval from 48% to 98% for neonatal ciprofloxacin simulations, respectively. Even after accounting for other critical population-specific maturational differences, novel RT ontogeny profiles substantially improved neonatal PBPK model performance, providing validated estimates of maturational differences in RT activity for optimal dosing in children.

Anakinra Removal by Continuous Renal Replacement Therapy: An Ex Vivo Analysis.

OBJECTIVES: Patients with sepsis are at significant risk for multiple organ dysfunction, including the lungs and kidneys. To manage the morbidity associated with kidney impairment, continuous renal replacement therapy (CRRT) may be required. The extent of anakinra pharmacokinetics in CRRT remains unknown. The objectives of this study were to investigate the anakinra-circuit interaction and quantify the rate of removal from plasma. DESIGN: The anakinra-circuit interaction was evaluated using a closed-loop ex vivo CRRT circuit. CRRT was performed in three phases based on the method of solute removal: 1) hemofiltration, 2) hemodialysis, and 3) hemodiafiltration. Standard control samples of anakinra were included to assess drug degradation. SETTING: University research laboratory. PATIENTS: None. INTERVENTIONS: Anakinra was administered to the CRRT circuit and serial prefilter blood samples were collected along with time-matched control and hemofiltrate samples. Each circuit was run in triplicate to assess inter-run variability. Concentrations of anakinra in each reference fluid were measured by enzyme-linked immunosorbent assay. Transmembrane filter clearance was estimated by the product of the sieving coefficient/dialysate saturation constant and circuit flow rates. MEASUREMENTS AND MAIN RESULTS: Removal of anakinra from plasma occurred within minutes for each CRRT modality. Average drug remaining (%) in plasma following anakinra administration was lowest with hemodiafiltration (34.9%). The average sieving coefficient was 0.34, 0.37, and 0.41 for hemodiafiltration, hemofiltration, and hemodialysis, respectively. Transmembrane clearance was fairly consistent across each modality with the highest during hemodialysis (5.53 mL/min), followed by hemodiafiltration (4.99 mL/min), and hemofiltration (3.94 mL/min). Percent drug remaining within the control samples (93.1%) remained consistent across each experiment, indicating negligible degradation within the blood. CONCLUSIONS: The results of this analysis are the first to demonstrate that large molecule therapeutic proteins such as anakinra, are removed from plasma with modern CRRT technology. Current dosing recommendations for patients with severe renal impairment may result in subtherapeutic anakinra concentrations in those receiving CRRT.

Cannabidiol Exposure Through Maternal Marijuana Use: Predictions in Breastfed Infants.

BACKGROUND AND OBJECTIVE: Knowledge about exposure to cannabidiol (CBD) in breastfed infants can provide an improved understanding of potential risk. The aim was to predict CBD exposure in breastfed infants from mothers taking CBD and CBD-containing products. METHODS: Cannabidiol concentrations in milk previously attained from data collected through an existing human milk research biorepository were used to simulate infant doses and identify subgroups. A developed pediatric physiologically based pharmacokinetic model produced virtual breastfed infants administered the simulated CBD doses. Predicted breastfed infant exposures and upper area under the curve ratios were compared to the lowest therapeutic dose for approved indications in children. RESULTS: The existing human milk research biorepository contained 200 samples from 181 unique breastfeeding mothers for whom self-reported administration data and CBD concentrations had previously been measured. Samples that were above the lower limit of quantification with only one maternal administration type revealed that administration type, i.e., joint/blunt or edible versus oil or pipe, resulted in significantly different subgroups in terms of milk concentrations. Resulting simulated infant doses (ng/kg) were described by lognormal distributions with geometric means and geometric standard deviations: 0.61 ± 2.41 all concentrations, 0.10 ± 0.37 joint/blunt or edible, and 2.23 ± 8.15 oil or pipe. Doses administered to breastfed infants had exposures magnitudes lower than exposures in children aged 4-11 years administered the lowest therapeutic dose for approved indications, and low upper area under the curve ratios. CONCLUSIONS: Based on real-world use, breastfeeding infants are predicted to receive very small exposures of CBD through milk. Studies examining adverse reactions will provide further insight into potential risk.

Addressing maternal medication use during breastfeeding using clinical resources and a novel physiologically based pharmacokinetic model-derived metric: A qualitative study.

Introduction: Breastfeeding has major benefits to the maternal-infant dyad and yet healthcare providers have expressed uncertainty about advocating breastfeeding when mothers are taking medications. The tendency for some providers to be more cautious in their advising approach is likely a consequence of limited, unfamiliar, and unreliable existing information on medication use during lactation. A novel risk metric termed the Upper Area Under the Curve Ratio (UAR) was developed to overcome existing resource shortcomings. However, the perception and use of the UAR in practice by providers is not known. The aim of this study was to understand existing resource use and potential UAR use in practice, their advantages and disadvantages, and areas of improvement for the UAR. Methods: Healthcare providers mainly practicing in California with experience advising on medication use during lactation were recruited. One-on-one semi-structured interviews that included questions on current practices when advising medication use during breastfeeding, and approaches to a given a scenario with and without information about the UAR were conducted. The Framework Method was applied for data analysis to construct themes and codes. Results: Twenty-eight providers representing multiple professions and disciplines were interviewed. Six main themes emerged: (1) Current Practice Approaches, (2) Advantages of Existing Resources, (3) Disadvantages of Existing Resources, (4) Advantages of the UAR, (5) Disadvantages of the UAR, and (6) Strategies to Improve the UAR. Overall, 108 codes were identified that illustrated theme topics ranging from a general lack of metric use to the realities of advising. A workflow describing current practice approaches connected all other themes. Almost all disadvantages of existing resources could be overcome by advantages of other resources and the UAR. Several improvements to the UAR were identified to address its shortcomings. Conclusion: Through interviews with providers who use resources to advise on medication use during breastfeeding, an improved understanding of current practice approaches and accessed resources was ascertained. Ultimately, it was found that the UAR would confer multiple benefits over existing resources, and improvements of the UAR were identified. Future work should focus on implementing the suggested recommendations to ensure optimal uptake of the UAR to improve advising practices.

Verifying in vitro-determined enzyme contributions to cannabidiol clearance for exposure predictions in human through physiologically-based pharmacokinetic modeling.

Cannabidiol (CBD) is approved for treatment of seizures associated with two forms of epilepsy that become apparent in infancy or early childhood. To consider an adult physiologically-based pharmacokinetic (PBPK) model for pediatric scaling, we assessed in vitro-derived cytochrome P450 (CYP) and uridine 5'-diphospho-glucuronosyltransferase (UGT) enzyme contributions to CBD clearance in human. An i.v. PBPK model was constructed using CBD physicochemical properties and knowledge of disposition. The i.v. datasets were used for model building and evaluation. Oral PBPK models for CBD administered in fasted and fed states were developed using single dose oral datasets and parameters optimized from the i.v. model and evaluated with multiple dose datasets. Relative contributions of CBD metabolizing enzymes were partitioned according to in vitro studies. Clinical drug-drug interaction (DDI) studies were simulated using CBD fed state, itraconazole, fluconazole, and rifampicin PBPK models. Linear mixed effect modeling was used to estimate area under the concentration-time curve from zero to infinity (AUC0-∞ ) perpetrator + CBD versus CBD alone. The i.v. and oral datasets used in model evaluation produced acceptable average fold error (AFE) of 1.28 and absolute AFE of 1.65. Relative contributions of drug-metabolizing enzymes to CBD clearance were proposed from in vitro data: UGT1A7 4%, UGT1A9 16%, UGT2B7 10%, CYP3A4 38%, CYP2C19 21%, and CYP2C9 11%. The simulated DDI studies using the in vitro-derived values produced AUC0-∞ treatment ratios comparable to observed: itraconazole 1.24 versus 1.07, fluconazole 1.45 versus 1.22, and rifampicin 0.49 versus 0.69. The constructed CBD PBPK models can predict adult exposures and have potential for use in pediatrics where exposure estimates are limited.

Cost-utility analysis of emicizumab for the treatment of severe hemophilia A patients in Canada.

INTRODUCTION: EHL FVIII products and emicizumab provide clinicians with other prophylactic options for treating hemophilia A, however, it is unclear if emicizumab is a cost-saving option. The objective of this study is to estimate the health and economic effects of using prophylactic EHL FVIII, SHL FVIII, and emicizumab in severe haemophilia A patients. MATERIALS AND METHODS: A state-transition Markov model evaluated the cost-effectiveness of prophylactic SHL FVIII, EHL FVIII, and emicizumab in a cohort of 2-year-old male patients over a lifetime horizon in the form of a cost-utility analysis using a Canadian provincial ministry of health payer perspective. The transition probabilities, costs, and utilities were obtained from literature and the Canadian Bleeding Disorders Registry. Probabilistic sensitivity and scenario analyses were performed to test the robustness of the model. RESULTS: The base-case analysis, over a lifetime horizon, resulted in a total cost and utilities per person for SHL FVIII, EHL FVIII, and emicizumab of $27.2 million (M), $36.7 M, and $26.2 M, respectively, and 31.30, 31.16, and 31.61 quality-adjusted life years, respectively. Emicizumab treatment resulted in 29 and 16 less bleeds in a lifetime compared to SHL FVIII and EHL FVIII, respectively. Probabilistic sensitivity analysis showed that emicizumab was cost-saving 100% of the time compared to SHL FVIII and EHL FVIII. CONCLUSION: The cost-utility analysis showed that emicizumab is more effective and may be less costly than FVIII for Canadian haemophilia A patients, conditional on drug cost assumptions. Our model indicates that emicizumab may be a potentially favourable treatment option for minimising healthcare costs and providing higher effectiveness.

A personalized limited sampling approach to better estimate terminal half-life of FVIII concentrates.

INTRODUCTION: Hemophilia A is a bleeding disorder characterized by a deficiency of a coagulation factor VIII and optimally treated using pharmacokinetics (PK)-guided prophylactic replacement therapy. To decrease patient burden, PK can be estimated from sparse sampling leveraging population PK modeling. However, recommendations for sampling times meant for patients with hemophilia A as a group may not be optimal at the individual level. OBJECTIVE: To evaluate a personalized limited sampling approach (Personalized LSA) that suggests a next sampling time point that would provide a more accurate estimation of terminal half-life of FVIII concentrates when using a population PK approach. METHODS: 331 PK studies with rich sampling were extracted from the WAPPS-Hemo database. Two sampling approaches were evaluated and compared: 974 PK studies consisting of two samples were built from the rich sampling data including one sample selected using the personalized LSA prediction; 974 PK studies consisting of two samples were built from the rich sampling data including one sample selected randomly. Half-life values were estimated on the sparse data and compared within patients to the estimates obtained on the rich data for assessing the error on half-life values. RESULTS: Relative errors between estimates from sparse sampling data using personalized LSA and from rich sampling data were always lower than 20% and significantly lower than the comparative approach that used random sampling (median-95th percentile were 3.8%-13.1% vs. 7.0%-23.5%, respectively, p-value < 10-10 ). Moreover, less than 4% of the samples suggested by the personalized LSA were below the limit of quantification. CONCLUSIONS: Identifying the most informative sampling points for PK assessment using a Personalized LSA approach that accounts for individual differences in PK improves the precision of FVIII terminal half-life estimates in sparse sampling.

Use of physiologically-based pharmacokinetic modeling to inform dosing of the opioid analgesics fentanyl and methadone in children with obesity.

Obesity is an increasingly alarming public health threat, with nearly 20% of children classified as obese in the United States today. Children with obesity are commonly prescribed the opioids fentanyl and methadone, and accurate dosing is critical to reducing the risk of serious adverse events associated with overexposure. However, pharmacokinetic studies in children with obesity are challenging to conduct, so there is limited information to guide fentanyl and methadone dosing in these children. To address this clinical knowledge gap, physiologically-based pharmacokinetic models of fentanyl and methadone were developed in adults and scaled to children with and without obesity to explore the interplay of obesity, age, and pharmacogenomics. These models included key obesity-induced changes in physiology and pharmacogenomic effects. Model predictions captured observed concentrations in children with obesity well, with an overall average fold error of 0.72 and 1.08 for fentanyl and methadone, respectively. Model simulations support a reduced fentanyl dose (1 vs. 2 µg/kg/h) starting at an earlier age (6 years) in virtual children with obesity, highlighting the importance of considering both age and obesity status when selecting an infusion rate most likely to achieve steady-state concentrations within the target range. Methadone dosing simulations highlight the importance of considering genotype in addition to obesity status when possible, as cytochrome P450 (CYP)2B6*6/*6 virtual children with obesity required half the dose to match the exposure of wildtype children without obesity. This physiologically-based pharmacokinetic modeling approach can be applied to explore dosing of other critical drugs in children with obesity.

Use of Real-World Data and Physiologically-Based Pharmacokinetic Modeling to Characterize Enoxaparin Disposition in Children With Obesity.

Dosing guidance for children with obesity is often unknown despite the fact that nearly 20% of US children are classified as obese. Enoxaparin, a commonly prescribed low-molecular-weight heparin, is dosed based on body weight irrespective of obesity status to achieve maximum concentration within a narrow therapeutic or prophylactic target range. However, whether children with and without obesity experience equivalent enoxaparin exposure remains unclear. To address this clinical question, 2,825 anti-activated factor X (anti-Xa) surrogate concentrations were collected from the electronic health records of 596 children, including those with obesity. Using linear mixed-effects regression models, we observed that 4-hour anti-Xa concentrations were statistically significantly different in children with and without obesity, even for children with the same absolute dose (P = 0.004). To further mechanistically explore obesity-associated differences in anti-Xa concentration, a pediatric physiologically-based pharmacokinetic (PBPK) model was developed in adults, and then scaled to children with and without obesity. This PBPK model incorporated binding of enoxaparin to antithrombin to form anti-Xa and elimination via heparinase-mediated metabolism and glomerular filtration. Following scaling, the PBPK model predicted real-world pediatric concentrations well, with an average fold error (standard deviation of the fold error) of 0.82 (0.23) and 0.87 (0.26) in children with and without obesity, respectively. PBPK model simulations revealed that children with obesity have at most 20% higher 4-hour anti-Xa concentrations under recommended, total body weight-based dosing compared to children without obesity owing to reduced weight-normalized clearance. Enoxaparin exposure was better matched across age groups and obesity status using fat-free mass weight-based dosing.

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

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

Predicting Individual Changes in Terminal Half-Life After Switching to Extended Half-Life Concentrates in Patients With Severe Hemophilia.

Predicting individual effects of switching from standard half-life (SHL) to extended half-life (EHL) FVIII/FIX concentrates is pivotal in clinical care, but large-scale individual data are scarce. The aim of this study was to assess individual changes in terminal half-life (THL) after switching to EHL concentrates and identifying determinants of a clinically relevant THL extension in people with severe hemophilia. Data from participants with pharmacokinetic studies on both SHL and EHL were extracted from the Web-Accessible Population Pharmacokinetics Service (WAPPS) database and stratified according to hemophilia type and age groups (children/adults). A 30% increase in THL was considered clinically relevant. Predictors of a relevant increase were identified using logistic regression. Data from 688 persons with severe hemophilia (2174 infusions) were included: 89% hemophilia A; median age: 21.7 (interquartile range [IQR]: 11.5-37.7); positive inhibitor history: 11.7%. THL increased by 38% (IQR: 17%-67%) and 212% (139%-367%) for hemophilia A and B, respectively. All EHL-FIX concentrate users showed clinically relevant THL extension. However, 40% (242/612) of people with hemophilia A showed limited extension or decrease in THL after switching. Relevant FVIII-THL extension was predicted by short baseline THL and blood group non-O in both children and adults. In conclusion, clinically relevant THL extension was observed in all 75/76 participants switching to EHL-FIX, and in 60% of 612 switching to EHL-FVIII. Short THL on SHL-FVIII and blood group non-O were identified as predictors for a relevant THL increase after switching to EHL-FVIII. Individualized pharmacokinetic assessment may guide clinical decision-making when switching from SHL to EHL-FVIII.

Physiologically Based Pharmacokinetic Modeling of Metformin in Children and Adolescents With Obesity.

Childhood obesity continues to rise in the United States and, with it, the off-label use of metformin for weight loss. The influence of age and obesity on the drug's disposition and exposure has not previously been studied using a mechanistic framework. Here, an adult physiologically based pharmacokinetic (PBPK) model of metformin was scaled to pediatric populations without obesity, with overweight/obesity, and with severe obesity; a published virtual population of children and adolescents with obesity was leveraged during model evaluation. When the pediatric model was simulated in groups aged 10 to 18 years, oral clearance following 1000 mg of metformin was higher (≈1200 mL/min) in those with obesity and severe obesity compared to the groups without and with overweight (≈1000 mL/min). In addition, simulated area under the concentration-time curve in older children and adolescents with obesity and severe obesity was comparable to that in adults with a similar dose-exposure relationship. Overall, simulations using the pediatric PBPK model support the use of adult doses of metformin in older children and adolescents with obesity. Moreover, the virtual population of children and adolescents with obesity offers a valuable tool to facilitate development of pediatric PBPK models for studying populations with obesity and, in turn, contribute information to inform drug labeling in this special population.

Model-Based Assessment of the Contribution of Monocytes and Macrophages to the Pharmacokinetics of Monoclonal Antibodies.

PURPOSE: We have hypothesized that a high concentration of circulating monocytes and macrophages may contribute to the fast weight-based clearance of monoclonal antibodies (mAbs) in young children. Exploring this hypothesis, this work uses modeling to clarify the role of monocytes and macrophages in the elimination of mAbs. METHODS: Leveraging pre-clinical data from mice, a minimal physiologically-based pharmacokinetic model was developed to characterize mAb uptake and FcRn-mediated recycling in circulating monocytes, macrophages, and endothelial cells. The model characterized IgG disposition in complex scenarios of site-specific FcRn deletion and variable endogenous IgG levels. Evaluation was performed for predicting IgG disposition with co-administration of high dose IVIG. A one-at-a-time sensitivity analysis quantified the role of relevant cellular parameters on IgG elimination in various scenarios. RESULTS: The plasma AUC of mAbs was highly sensitive to endothelial cell parameters, but had near-nil sensitivity to monocyte and macrophage parameters, even in scenarios with 90% loss of FcRn expression/activity. In mice with normal FcRn expression, simulations suggest that less than 2% of an IV dose is eliminated in macrophages, while endothelial cells are predicted to dominate mAb elimination. CONCLUSIONS: The model suggests that the role of monocytes and macrophages in IgG homeostasis includes extensive uptake and highly efficient FcRn-mediated protection, but not appreciable degradation when FcRn is present. Therefore, it is very unlikely that a high concentration of circulating monocytes can contribute to explaining the fast weight-based clearance of mAbs in very young children, even if FcRn expression/activity was 90% lower in children than in adults.

Development and Evaluation of an In Silico Dermal Absorption Model Relevant for Children.

The higher skin surface area to body weight ratio in children and the prematurity of skin in neonates may lead to higher chemical exposure as compared to adults. The objectives of this study were: (i) to provide a comprehensive review of the age-dependent anatomical and physiological changes in pediatric skin, and (ii) to construct and evaluate an age-dependent pediatric dermal absorption model. A comprehensive review was conducted to gather data quantifying the differences in the anatomy and physiology of child and adult skin. Maturation functions were developed for model parameters that were found to be age-dependent. A pediatric dermal absorption model was constructed by updating a MoBi implementation of the Dancik et al. 2013 skin permeation model with these maturation functions. Using a workflow for adult-to-child model extrapolation, the predictive performance of the model was evaluated by comparing its predicted rates of flux of diamorphine, phenobarbital and buprenorphine against experimental observations using neonatal skin. For diamorphine and phenobarbital, the model provided reasonable predictions. The ratios of predicted:observed flux in neonates for diamorphine ranged from 0.55 to 1.40. For phenobarbital, the ratios ranged from 0.93 to 1.26. For buprenorphine, the model showed acceptable predictive performance. Overall, the physiologically based pediatric dermal absorption model demonstrated satisfactory prediction accuracy. The prediction of dermal absorption in neonates using a model-based approach will be useful for both drug development and human health risk assessment.

Development and Evaluation of a Virtual Population of Children with Obesity for Physiologically Based Pharmacokinetic Modeling.

BACKGROUND AND OBJECTIVE: While one in five children in the USA are now obese, and more than three-quarters receive at least one drug during childhood, there is limited dosing guidance for this vulnerable patient population. Physiologically based pharmacokinetic modeling can bridge the gap in the understanding of how pharmacokinetics, including drug distribution and clearance, changes with obesity by incorporating known obesity-related physiological changes in children. The objective of this study was to develop a virtual population of children with obesity to enable physiologically based pharmacokinetic modeling, then use the novel virtual population in conjunction with previously developed models of clindamycin and trimethoprim/sulfamethoxazole to better understand dosing of these drugs in children with obesity. METHODS: To enable physiologically based pharmacokinetic modeling, a virtual population of children with obesity was developed using national survey, electronic health record, and clinical trial data, as well as data extracted from the literature. The virtual population accounts for key obesity-related changes in physiology relevant to pharmacokinetics, including increased body size, body composition, organ size and blood flow, plasma protein concentrations, and glomerular filtration rate. The virtual population was then used to predict the pharmacokinetics of clindamycin and trimethoprim/sulfamethoxazole in children with obesity using previously developed physiologically based pharmacokinetic models. RESULTS: Model simulations predicted observed concentrations well, with an overall average fold error of 1.09, 1.24, and 1.53 for clindamycin, trimethoprim, and sulfamethoxazole, respectively. Relative to children without obesity, children with obesity experienced decreased clindamycin and trimethoprim/sulfamethoxazole weight-normalized clearance and volume of distribution, and higher absolute doses under recommended pediatric weight-based dosing regimens. CONCLUSIONS: Model simulations support current recommended weight-based dosing in children with obesity for clindamycin and trimethoprim/sulfamethoxazole, as they met target exposure despite these changes in clearance and volume of distribution.