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.

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.

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.

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.

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.

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.

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.

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.

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 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.

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 10 - 18 y of age, oral clearance (CL/F) following 1,000 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 AUC 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. This article is protected by copyright. All rights reserved.

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.

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.

Physiologically Based Pharmacokinetic Modeling of Meropenem in Preterm and Term Infants

BACKGROUND

Meropenem is a broad-spectrum carbapenem antibiotic approved by the US Food and Drug Administration for use in pediatric patients, including treating complicated intra-abdominal infections in infants < 3 months of age. The impact of maturation in glomerular filtration rate and tubular secretion by renal transporters on meropenem pharmacokinetics, and the effect on meropenem dosing, remains unknown. We applied physiologically based pharmacokinetic (PBPK) modeling to characterize the disposition of meropenem in preterm and term infants.

METHODS

An adult meropenem PBPK model was developed in PK-Sim® (Version 8) and scaled to infants accounting for renal transporter ontogeny and glomerular filtration rate maturation. The PBPK model was evaluated using 645 plasma concentrations from 181 infants (gestational age 23-40 weeks; postnatal age 1-95 days). The PBPK model-based simulations were performed to evaluate meropenem dosing in the product label for infants < 3 months of age treated for complicated intra-abdominal infections.

RESULTS

Our model predicted plasma concentrations in infants in agreement with the observed data (average fold error of 0.90). The PBPK model-predicted clearance in a virtual infant population was successfully able to capture the post hoc estimated clearance of meropenem in this population, estimated by a previously published model. For 90% of virtual infants, a 4-mg/L target plasma concentration was achieved for > 50% of the dosing interval following product label-recommended dosing.

CONCLUSIONS

Our PBPK model supports the meropenem dosing regimens recommended in the product label for infants <3 months of age.

A Physiological Approach to Pharmacokinetics in Chronic Kidney Disease

The conventional approach to approximating the pharmacokinetics of drugs in patients with chronic kidney disease (CKD) only accounts for changes in the estimated glomerular filtration rate. However, CKD is a systemic and multifaceted disease that alters many body systems. Therefore, the objective of this exercise was to develop and evaluate a whole-body mechanistic approach to predicting pharmacokinetics in patients with CKD. Physiologically based pharmacokinetic models were developed in PK-Sim v8.0 (www.open-systems-pharmacology.org) to mechanistically represent the disposition of 7 compounds in healthy human adults. The 7 compounds selected were eliminated by glomerular filtration and active tubular secretion by the organic cation transport system to varying degrees. After a literature search, the healthy adult models were adapted to patients with CKD by numerically accounting for changes in glomerular filtration rate, kidney volume, renal perfusion, hematocrit, plasma protein concentrations, and gastrointestinal transit. Literature-informed interindividual variability was applied to the physiological parameters to facilitate a population approach. Model performance in CKD was evaluated against pharmacokinetic data from 8 clinical trials in the literature. Overall, integration of the CKD parameterization enabled exposure predictions that were within 1.5-fold error across all compounds and patients with varying stages of renal impairment. Notable improvement was observed over the conventional approach to scaling exposure, which failed in all but 1 scenario in patients with advanced CKD. Further research is required to qualify its use for first-in-CKD dose selection and clinical trial planning for a wider selection of renally eliminated compounds, including those subject to anion transport.

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.

Terminal half-life of FVIII and FIX according to age, blood group and concentrate type: Data from the WAPPS database

BACKGROUND

Real-life data on pharmacokinetics of factor (F) VIII/IX concentrates, especially extended half-life (EHL), concentrates in large cohorts of persons with hemophilia are currently lacking.

OBJECTIVES

This cross-sectional study aimed to establish reference values for terminal half-life (THL) for FVIII/IX concentrates according to concentrate type, age, blood group and inhibitor history.

PATIENTS/METHODS

Data were extracted from the Web-Accessible Population Pharmacokinetics Service database. Groups were compared by nonparametric tests. THL was modelled according to patient characteristics and concentrate type.

RESULTS

Infusion data (n = 8022) were collected from 4832 subjects (including 2222 children) with severe hemophilia (age: 1 month-85 years; 89% hemophilia A; 34% using EHL concentrates, 9.8% with history of inhibitors). THL of FVIII-EHL was longer than of FVIII standard half-life (SHL; median 15.1 vs. 11.1 h). FVIII-THL was dependent on age, concentrate type, blood group, and inhibitor history. THL of FIX-EHL was longer than of FIX-SHL (median 106.9 vs. 36.5 h). FIX-THL increased with age until 30 years and remained stable thereafter. FVIII-THL was shorter in subjects with blood group O. THL was decreased by 1.3 h for FVIII and 22 h for FIX in subjects with a positive inhibitor history.

CONCLUSIONS

We established reference values for FVIII/IX concentrates according to patient characteristics and concentrate type in a large database of hemophilia patients. These reference values may inform clinical practice (e.g., assessment of immune tolerance success), economic implications of procurement processes and value attribution of novel treatments (e.g., mimetics, gene therapy).

Correction to: Incorporating Breastfeeding-Related Variability with Physiologically Based Pharmacokinetic Modeling to Predict Infant Exposure to Maternal Medication Through Breast Milk: a Workflow Applied to Lamotrigine

A Correction to this paper has been published: 10.1208/s12248-021-00615-8

Correction to: Incorporating Breastfeeding-Related Variability with Physiologically Based Pharmacokinetic Modeling to Predict Infant Exposure to Maternal Medication Through Breast Milk: a Workflow Applied to Lamotrigine

A Correction to this paper has been published: https://doi.org/10.1208/s12248-021-00614-9.

Pediatric Dose Selection for Therapeutic Proteins

In selecting optimal dosing regimens in support of the clinical use of monoclonal antibodies and other therapeutic proteins in pediatric indications, the unique pharmacokinetic properties of this class of biologics, as well as the underlying physiologic and pathophysiologic processes and their modulation by childhood growth and development, needs to be appreciated. During drug development, first-in-pediatric dose selection is a capstone event in the pediatric investigation plan that relies heavily on extrapolation of pharmacokinetic and pharmacodynamic data from adult to pediatric populations. It is facilitated by combinations of pharmacometric approaches, including allometry, physiologically based pharmacokinetic modeling, and population pharmacokinetic analyses, although data on reliability and qualification of some of these tools in the context of therapeutic proteins are still limited but emerging. Presented data suggest nonlinear relationships between body weight and both clearance and volume of distribution for therapeutic proteins in pediatric populations, with allometric exponents of 0.75 and 0.8, respectively. For newborns and infants (<1 year), even higher nonlinearity seems to occur. Translation of the quantitative characterization of the pediatric pharmacokinetics of therapeutic proteins into dosing regimens for the drug label requires compromising between precision dosing and clinical practicability, with tiered dosing algorithms based on size or age strata being the currently most frequently applied methodology.

Incorporating Breastfeeding-Related Variability with Physiologically Based Pharmacokinetic Modeling to Predict Infant Exposure to Maternal Medication Through Breast Milk: a Workflow Applied to Lamotrigine

Current methods to assess risk in infants exposed to maternal medication through breast milk do not specifically account for infants most vulnerable to high drug exposure. A workflow applied to lamotrigine incorporated variability in infant anatomy and physiology, milk intake volume, and milk concentration to predict infant exposure. An adult physiologically based pharmacokinetic model of lamotrigine was developed and evaluated. The model was scaled to account for growth and maturation of a virtual infant population (n=100). Daily infant doses were simulated using milk intake volume and concentration models described by a nonlinear equation of weight-normalized intake across infant age and a linear function on the relationship of observed milk concentrations and maternal doses, respectively. Average infant plasma concentration at steady state was obtained through simulation. Models were evaluated by comparing observed to simulated infant plasma concentrations from breastfeeding infants based on a 90% prediction interval (PI). Upper AUC ratio (UAR) was defined as a novel risk metric. Twenty-five paired (milk concentrations measured) and 18 unpaired (milk concentrations unknown) infant plasma samples were retrieved from the literature. Forty-four percent and 11% of the paired and unpaired infant plasma concentrations were outside of the 90% PI, respectively. Over all ages (0-7 months), unpaired predictions captured more observed infant plasma concentrations within 90% PI than paired. UAR was 0.18-0.44 when mothers received 200 mg lamotrigine, suggesting that infants can receive 18-44% of the exposure per dose as compared to adults. UARs determined for further medications could reveal trends to better classify at-risk mother-infant pairs.

Pharmacokinetic implications of dosing emicizumab based on vial size: A simulation study

INTRODUCTION

Emicizumab is dosed as mg/kg and, according to the label, any unused drug left in the vial(s) must be discarded, thereby wasting expensive resources. The aim of this study was to use population pharmacokinetics to illustrate the implications of changing the dosing interval to avoid wastage.

METHODS

We used a previously published emicizumab PopPK model after extending its validation to children. We simulated PK parameters for labelled dosing regimens and for regimens using full vials with infusion frequency varied to keep the steady-state drug concentration unchanged. Cost and drug savings were calculated.

RESULTS

The model evaluation was successful. When rounding up, the average individual below 53, 47 and 39 has a time-to-trough increase of up to 5.7, 7.9 and 5.8 days for the QW, Q2 W and Q4 W regimen, respectively. This resulted in an annual cost reduction of up to $173,136, $75,747 and $61,319 USD per patient. At higher body weights, rounding down the dose to the nearest vial resulted in negligible changes in the steady state concentration and cost savings of up to $93,781, $46,891 and $23,446 USD per patient, respectively.

CONCLUSION

Individuals with a lower body weight may benefit from increasing dose intervals and rounding up dose up to the nearest vial, and individuals with a higher body weight from maintaining the injection frequency and rounding dose down to the nearest vial without significant change in emicizumab levels. Administering the entire vial may result in a reduction of vials used annually and potential cost savings.

Evaluation of models for predicting pediatric fraction unbound in plasma for human health risk assessment

Pediatric physiologically based pharmacokinetic (PBPK) models facilitate the prediction of PK parameters in children under specific exposure conditions. Pharmacokinetic outcomes are highly sensitive to fraction unbound in plasma (fup) as incorporated into PBPK models. Rarely is fup in children (fup_child) experimentally derived and prediction is based upon fup in adults (fup_adult) as well as a ratio of plasma protein concentrations between children and adults. The objectives were to (i) evaluate protein concentration vs. age profile derived from ontogeny models, (ii) assess predictive performances of fup ontogeny models, and (iii) determine overall uncertainty in fup_child prediction resulting from a combination of quantitative structure-property relationship (QSPR) model and ontogeny models. The plasma albumin and alpha-acid glycoprotein (AAG) concentration data for pediatrics and fup_child and fup_adult data were obtained from literature. The protein concentration vs. age profile derived from ontogeny models were compared to observed levels. Fup_child values were calculated according to ontogeny models using both observed and QSPR-predicted fup_adult as inputs and predictive performances of ontogeny models assessed by comparing predicted fup_child to observed values. Protein concentrations vs. age profiles derived from non-linear equations were more congruent with observed albumin levels than linear or step-wise models. When observed fup_adult values were used as input, the fup_child data were under-predicted with average fold error (AFE) amounts ranging 0.79-0.81 and 0.77-0.97 for albumin and AAG ontogeny models, respectively. When QSPR-predicted fup_adult values were used as input, AFE of fup_child ranged 1.2-1.35 and 0.98-1.2 for albumin and AAG models, respectively. The choice of ontogeny model with respect to prediction accuracy is more important for AAG, highly bound compounds and infants. For these compounds and scenarios, experimental determination of fup_child for inclusion into a pediatric PBPK model is necessary to have confidence in PBPK model outputs.

A comparison of methods for prediction of pharmacokinetics when switching to extended half-life products in hemophilia A patients

INTRODUCTION

Hemophilia A is a genetic bleeding disorder resulting from a lack of clotting factor VIII. Where extended half-life products are available, people with hemophilia may stop their current drug regimen and switch to a EHL product providing a more convenient dosing regimen. While most factor VIII concentrate regimens are started prophylactically based on international units per weight, this "one-size-fits-all" approach does not account for the large pharmacokinetic variability between individuals.

AIMS

We explored methods to predict individual PK of an EHL product by using population pharmacokinetic models and eta-values (η), a value that quantifies how individuals deviate from a population for any PK parameter, derived from a prior product. In addition, we wanted to investigate which individuals would benefit from this method compared to using a PopPK model alone.

METHODS

PK data from subjects (n = 39) who have taken both Adynovate and Eloctate was collected from clinical trial data and from the Web-Accessible Population Pharmacokinetic Service - Hemophilia (WAPPS-Hemo) database. In addition, PK data from subjects (n = 200) who switched from a standard half-life product to Eloctate was also extracted from the WAPPS-Hemo database. Two methods to estimate individual PK outcomes of the second product were compared. The PopPK method used the Eloctate PopPK model published from WAPPS-Hemo, while the η-method incorporated individually scaled η from the prior product's PopPK model. Both methods were assessed for its performance in predicting PK outcomes. Absolute percent differences were calculated between the predicted and observed PK outcomes. Infusions were parsed into subgroups based on number of samples and individual η-percentiles for analysis.

RESULTS

For the three switching protocols (Adynovate to Eloctate, Eloctate to Adynovate, and SHL FVIII to Eloctate), the η-method resulted in a relative difference reduction in mean absolute percent difference of 27.8% (range 1-59%), 4.9% (range 0-129%), and 18.0% (0-79%) in half-life compared to the PopPK method respectively. With some exceptions (in particular central volume), the η-method produced relative difference reduction in mean absolute percent differences up to 33% lower compared to the PopPK method. When individuals were parsed based on their η-values (either CL or V1), the two methods differentiate up to 64% in terms of half-life and time to 0.02 IU/mL predictions for individuals with a low (0th to 20th percentile) η_CL or η_V1 on the first product. Individuals with higher number of observations per infusion on the first product resulted in better predictions in PK parameter estimates when using the η-method.

CONCLUSION

The use of prior knowledge by implementing η-values into PopPK models may provide clinicians with a safer and more effective method to choose a dosing regimen for patients with hemophilia A switching from one factor concentrate to another. However, the η-method was unable to better predict an increase or decrease in half-life of a future product compared to the PopPK method, and thus supports the conclusion that most individuals would still benefit from a trial on the EHL and subsequent estimation of their individual PK profile from sparse measurements on the EHL.

Physiologically Based Pharmacokinetic Modeling for Trimethoprim and Sulfamethoxazole in Children

OBJECTIVE

The aims of this study were to (1) determine whether opportunistically collected data can be used to develop physiologically based pharmacokinetic (PBPK) models in pediatric patients; and (2) characterize age-related maturational changes in drug disposition for the renally eliminated and hepatically metabolized antibiotic trimethoprim (TMP)-sulfamethoxazole (SMX).

METHODS

We developed separate population PBPK models for TMP and SMX in children after oral administration of the combined TMP-SMX product and used sparse and opportunistically collected plasma concentration samples to validate our pediatric model. We evaluated predictability of the pediatric PBPK model based on the number of observed pediatric data out of the 90% prediction interval. We performed dosing simulations to target organ and tissue (skin) concentrations greater than the methicillin-resistant Staphylococcus aureus (MRSA) minimum inhibitory concentration (TMP 2 mg/L; SMX 9.5 mg/L) for at least 50% of the dosing interval.

RESULTS

We found 67-87% and 71-91% of the observed data for TMP and SMX, respectively, were captured within the 90% prediction interval across five age groups, suggesting adequate fit of our model. Our model-rederived optimal dosing of TMP at the target tissue was in the range of recommended dosing for TMP-SMX in children in all age groups by current guidelines for the treatment of MRSA.

CONCLUSION

We successfully developed a pediatric PBPK model of the combination antibiotic TMP-SMX using sparse and opportunistic pediatric pharmacokinetic samples. This novel and efficient approach has the potential to expand the use of PBPK modeling in pediatric drug development.

Quantifying breast milk intake by term and preterm infants for input into paediatric physiologically based pharmacokinetic models

Despite the many benefits of breast milk, mothers taking medication are often uncertain about the risks of drug exposure to their infants and decide not to breastfeed. Physiologically based pharmacokinetic models can contribute to drug‐in‐milk safety assessments by predicting the infant exposure and subsequently, risk for toxic effects that would result from continuous breastfeeding. This review aimed to quantify breast milk intake feeding parameters in term and preterm infants using literature data for input into paediatric physiologically based pharmacokinetic models designed for drug‐in‐milk risk assessment. Ovid MEDLINE and Embase were searched up to July 2, 2019. Key study reference lists and grey literature were reviewed. Title, abstract and full text were screened in nonduplicate. Daily weight‐normalized human milk intake (WHMI) and feeding frequency by age were extracted. The review process retrieved 52 studies. A nonlinear regression equation was constructed to describe the WHMI of exclusively breastfed term infants from birth to 1 year of age. In all cases, preterm infants fed with similar feeding parameters to term infants on a weight‐normalized basis. Maximum WHMI was 152.6 ml/kg/day at 19.7 days, and weighted mean feeding frequency was 7.7 feeds/day. Existing methods for approximating breast milk intake were refined by using a comprehensive set of literature data to describe WHMI and feeding frequency. Milk feeding parameters were quantified for preterm infants, a vulnerable population at risk for high drug exposure and toxic effects. A high‐risk period of exposure at 2–4 weeks of age was identified and can inform future drug‐in‐milk risk assessments.

Clinical application of Web Accessible Population Pharmacokinetic Service-Hemophilia (WAPPS-Hemo): Patterns of blood sampling and patient characteristics among clinician users

BACKGROUND

Use of population pharmacokinetics (PopPK) to facilitate PK-informed prophylaxis in clinical practice has gained momentum among haemophilia providers due to the accessibility of tools such as the Web Accessible Population Pharmacokinetic Service-Hemophilia (WAPPS-Hemo) and availability of extended half-life (EHL) factor concentrates. It is unknown how clinicians implement PopPK.

AIM

To investigate the evolution of PopPK use in clinical practice by comparing blood sampling strategies, patient features, and factor group between initial and recent periods of WAPPS-Hemo availability.

METHODS

PK data for haemophilia A and haemophilia B patients from two time periods were extracted from the WAPPS-Hemo database: early availability (10/2015-09/2016) and recent use (10/2017-09/2018). We compared patient characteristics (age, body weight, haemophilia type), product type and dose, and blood sampling times between the time frames.

RESULTS

Over 1900 eligible infusions were submitted to WAPPS-Hemo during the periods studied, with 85% representing FVIII concentrates. In the recent cohort, PK profiles were requested for younger patients (median age 18 vs 26 years), with increased proportional EHL FVIII use (29% vs 14% of infusions). High-use centres generally submitted fewer blood samples per infusion than non-high-use centres, although the number of samples collected by non-high-use centres decreased significantly over time. During both periods, blood sample timing was generally consistent with ISTH recommended windows.

CONCLUSION

The use of WAPPS-Hemo by haemophilia providers grew by over threefold between the time periods investigated. While sampling times have included key time points proposed first by Björkman since early WAPPS-Hemo usage, a trend towards minimizing sampling was observed.

A comparison of methods for prediction of pharmacokinetics across factor concentrate switching in hemophilia patients

INTRODUCTION

This study proposes a method to predict individual pharmacokinetics of a future product by using the individual pharmacokinetic profile on the current product and the PopPK models of the current and future product.

METHODS

Individual dense data was collected from two PK crossover studies, one enrolling 29 patients switching from Advate to Eloctate and one enrolling 15 patients switching from Advate to Novoeight. Three methods were designed to predict the second product's individual PK parameters (CL, V1, Q, and V2). Method 1 used the second product's typical population value of PK parameters from its PopPK model. Method 2 used the second product's calculated PK parameters based on individual covariates and its PopPK model. Method 3 used method 2, along with the predicted η-values of CL and V1 from the first product and its PopPK model. Each method was used to assess PK prediction during switching from Advate to Novoeight, Novoeight to Advate, and Advate to Eloctate.

RESULTS

The three methods produced different outcomes. The mean absolute relative errors for half-life were lowest for method 3 for each study (11.6%, 13.1%, 13.6%). The regression line between predicted and observed half-life for method 3 was closest to the line of identity for each study (0.84, 0.67, 0.66).

CONCLUSION

Taking into account individual PK from a previous clotting factor product was shown to provide better means of estimating individual PK for a new product. This may improve regimen design across switches and reduce the time to tailor optimal dose of FVIII products.

Predicting Escitalopram Exposure to Breastfeeding Infants: Integrating Analytical and In Silico Techniques

BACKGROUND

Escitalopram is used for post-partum depression; however, there are limited pharmacokinetic data of escitalopram in milk and plasma of infants breastfed by women taking the drug.

OBJECTIVE

The objective of this study was to apply physiologically-based pharmacokinetic (PBPK) modelling to predict infant drug exposure (plasma area under the curve from time zero to infinity [AUC_∞]) based on drug monitoring data of escitalopram in breast milk.

METHODS

Using a newly developed liquid chromatography-tandem mass spectrometry (LC-MS/MS) method, we quantified escitalopram concentrations in milk samples of 18 breastfeeding women with escitalopram therapy at steady state, collected at three to five time points. The escitalopram concentrations in breast milk were used with infant feeding parameters from the literature to simulate infant daily dose. We used PK-Sim^® to develop an adult PBPK model for escitalopram and extrapolated it to a population of 1600 infants up to 12 months of age. An integration of the simulated infant daily dose and the virtual infants with variable physiological-pharmacological parameters was used to predict drug exposure (plasma AUC_∞) distribution in the population of infants breastfed by women receiving escitalopram 20 mg/day.

RESULTS

Escitalopram concentrations in milk were 50 ± 17 ng/mL (mean ± standard deviation). The simulated infant plasma AUC_∞ following escitalopram exposure through breast milk was low, with a median of 1.7% (range 0.5-5.9%) of the corresponding maternal plasma AUC_∞, indicating no substantial exposure.

CONCLUSIONS

Infant exposure levels to escitalopram in breast milk are low. A PBPK modeling approach can be used to translate data on drug monitoring in milk into a population distribution of infant plasma levels for drug safety assessment.

Physiologically-Based Pharmacokinetic Modeling vs. Allometric Scaling for the Prediction of Infliximab Pharmacokinetics in Pediatric Patients

The comparative performances of physiologically‐based pharmacokinetic (PBPK) modeling and allometric scaling for predicting the pharmacokinetics (PKs) of large molecules in pediatrics are unknown. Therefore, both methods were evaluated for accuracy in translating knowledge of infliximab PKs from adults to children. PBPK modeling was performed using the base model for large molecules in PK‐Sim version 7.4 with modifications in Mobi. Eight population PK models from literature were reconstructed and scaled by allometry to pediatrics. Evaluation data included seven pediatric studies (~4–18 years). Both methods performed comparably with 66.7% and 68.6% of model‐predicted concentrations falling within twofold of the observed concentrations for PBPK modeling and allometry, respectively. Considerable variability was noted among the allometric models. Therefore, pediatric clinical trial planning would benefit from using approaches that require predictions depending on the specific question i.e., PBPK modeling and allometry.

Model qualification of the PK-Sim® pediatric module for pediatric exposure assessment of CYP450 metabolized compounds

Pediatric physiologically based pharmacokinetic (PBPK) models facilitate the estimation of pharmacokinetic (PK) parameters in children under specific exposure conditions. In human health risk assessment, PBPK modeling has been used to determine a chemical-specific human kinetic adjustment factor (HKAF). Due to increased demands in regulatory assessment, model evaluation and qualification have gained growing attention. The aim of this study was to undertake model qualification of pediatric PBPK models for compounds that are primarily metabolized by cytochrome P450 (CYP) enzymes. The objectives were to determine the appropriateness of the virtual individual creating algorithm in PK-Sim® in predicting PK parameters and their variability in children and identify critical system-specific inputs. PBPK models in adults were constructed for several pharmaceuticals (grouped by major clearance process such as CYP3A4). Several age groups of virtual individuals were created to represent children in pediatric clinical studies. The mean and variance of clearance (CL) from virtual populations were compared to observed values. Sensitivity analysis on area under the curve (AUC) was performed. System-specific parameters of virtual children that contribute to inter-individual PK properties were assessed. Eighty-one percent of the comparisons between simulated and observed clearance values were within twofold error. The mean fold errors were 1.1, 1, 0.7 and 1.8 in adolescents, children, infants and neonates, respectively. CL variability was reasonably predicted for 70% of the comparisons with comparable coefficients of variation between observed and predicted. The sensitivity analysis revealed that fraction unbound in plasma, parameters related to CYP enzyme-mediated metabolism and liver volumewere most important in the estimation of pediatric exposure. A comparison of variabilities in weight, height and liver volume in virtual children showed reliable agreement with observed data. The presented results of predictive performance and properties of virtual populations provide confidence in the use of PK-Sim for pediatric PBPK modeling in toxicological applications including PBPK-based-HKAF derivation.

Comparative pharmacokinetics of two extended half-life FVIII concentrates (Eloctate and Adynovate) in adolescents with hemophilia A: Is there a difference?

Essentials The PK parameters of Eloctate vs Adynovate were compared using one-stage and chromogenic assays in 25 boys (12-18 years). The FVIII levels were taken at 3, 24, 48, and 72 hours following a dose of either FVIII; levels analyzed by WAPPS PK program. The PK profiles (half-life, clearance, and time to 5%, 3%, and 1%) were not statistically different for the two EHL FVIIIs. The significant interpatient variability in PK is mainly related to VWF levels (and blood group).

BACKGROUND

A head-to-head comparison of the pharmokinetcs (PK) of extended half-life (EHL) factor VIII (FVIII) concentrates in the same subjects has not been reported. Recently, boys (ages 12-18 years) with hemophilia A in Canada were required to switch from Eloctate to Adynovate.

OBJECTIVES

Compare the PK profiles of Eloctate vs Adynovate in the same boys.

METHODS

Boys switching from Eloctate to Adynovate prophylaxis had FVIII levels sampled at 3, 24, 48, and 72 hours following a regular prophylactic infusion of Eloctate and then 1-3 months later, of Adynovate. Testing was done by one-stage assay (OSA) and chromogenic assay (CA). The PK parameters were determined with the Web Accessible Population Pharmacokinetic Service (WAPPS)-Hemo PK tool.

RESULTS

Twenty-five boys (mean age 15.3 years; range: 12.1-18.4; 9 O blood group) underwent switching. Mean (range) terminal half-lives with the OSA were 16.1 hours (10.4 to 23.4; Eloctate) and 16.7 hours (11.0 to 23.6; Adynovate) (NS). With the CA, these were 18.0 hours (12.0 to 25.5; Eloctate) and 16.0 hours (10.3 to 22.9; Adynovate) (P = 0.001). There were no significant differences between the two EHL-FVIIIs in clearance, area under the concentration vs time curve (AUC), Vss, or time for FVIII levels to drop to 5%, 3%, and 1%. At the 72-h time point, mean observed FVIII levels following a mean dose of 39.3 IU/kg of Eloctate were 4.4% (OSA) and 4.4% (CA). For Adynovate, these were 5.1% (OSA) and 5.3% (CA) following similar doses. There was considerable interpatient variation in PK, mainly explained by differences in blood group/von Willebrand factor (VWF) levels.

CONCLUSIONS

Eloctate and Adynovate have almost identical PK parameters. When switching from one to another no prophylaxis regimen change is needed.

Using pharmacokinetics for tailoring prophylaxis in people with hemophilia switching between clotting factor products: A scoping review

The objective of this scoping review is to summarize the current use of pharmacokinetics for tailoring prophylaxis in hemophilia patients switching between clotting factor products. Patients with hemophilia may require switching of clotting factor concentrates due to a variety of factors, but there have been perceived risks associated with switching, such as inhibitor development or suboptimal protection due to inadequate dosing while titrating treatment. Studies that look at patients switching from one clotting factor concentrate to another are categorized in terms of their primary and/or secondary objectives, notably biosimilarity and comparative pharmacokinetic studies and inhibitor development studies. Research on how best to switch concentrates with respect to dosing regimen are lacking, and currently a trial‐and‐error approach is used for dosing the new factor concentrate. In the future, studies looking at the predictability of pharmacokinetics (PK) of a new factor concentrate based on individual PK knowledge of the original factor concentrate may offer clinical benefit by providing a safer switching approach and protocol.

Biodistribution and Physiologically-Based Pharmacokinetic Modeling of Gold Nanoparticles in Mice with Interspecies Extrapolation

Gold nanoparticles (AuNPs) are a focus of growing medical research applications due to their unique chemical, electrical and optical properties. Because of uncertain toxicity, "green" synthesis methods are emerging, using plant extracts to improve biological and environmental compatibility. Here we explore the biodistribution of green AuNPs in mice and prepare a physiologically-based pharmacokinetic (PBPK) model to guide interspecies extrapolation. Monodisperse AuNPs were synthesized and capped with epigallocatechin gallate (EGCG) and curcumin. 64 CD-1 mice received the AuNPs by intraperitoneal injection. To assess biodistribution, groups of six mice were sacrificed at 1, 7, 14, 28 and 56 days, and their organs were analyzed for gold content using inductively coupled plasma mass spectrometry (ICP-MS). A physiologically-based pharmacokinetic (PBPK) model was developed to describe the biodistribution data in mice. To assess the potential for interspecies extrapolation, organism-specific parameters in the model were adapted to represent rats, and the rat PBPK model was subsequently evaluated with PK data for citrate-capped AuNPs from literature. The liver and spleen displayed strong uptake, and the PBPK model suggested that extravasation and phagocytosis were key drivers. Organ predictions following interspecies extrapolation were successful for rats receiving citrate-capped AuNPs. This work lays the foundation for the pre-clinical extrapolation of the pharmacokinetics of AuNPs from mice to larger species.

Impact of Adopting Population Pharmacokinetics for Tailoring Prophylaxis in Haemophilia A Patients: A Historically Controlled Observational Study

BACKGROUND

 Performing individual pharmacokinetics (PK) studies in clinical practice can be simplified by adopting population PK-based profiling on limited post-infusion samples. The objective of this study was to assess the impact of population PK in tailoring prophylaxis in patients with haemophilia A.

PATIENTS AND METHODS

 Individual weekly treatment plans were developed considering predicted plasma factor activity levels and patients' lifestyle. Patients were trained using a visual traffic-light scheme to help modulate their level of physical activity with respect to factor infusions timing. Annualized joint bleeding rate (ABJR), haemophilia-specific quality of life questionnaire for adults (Haemo-QoL-A) and factor utilization were measured for 12 months before and after tailoring, compared within patients and analysed separately for those previously on prophylaxis (P), situational prophylaxis (SP) or on-demand (OD).

RESULTS

 Sixteen patients previously on P, 10 on SP and 10 on OD were enrolled in the study. The median (lower, upper quartile) ABJR changed from 2.0 (0, 4.0) to 0 (0, 1.6) for P (p = 0.003), from 2.0 (2.0, 13.6) to 3.0 (1.4, 7.2) for SP (p = 0.183) and from 16.0 (13.0, 25.0) to 2.3 (0, 5.0) for OD (p = 0.003). The Haemo-QoL-A total score improved for 58% of P, 50% of SP and 29% of OD patients. Factor utilization (IU/kg/patient/year) increased by 2,400 (121; 2,586) for P, 1,052 (308; 1,578) for SP and 2,086 (1,498; 2,576) for OD. One of 138 measurements demonstrated a factor activity level below the critical threshold of 0.03 IU/mL while the predicted level was above the threshold.

CONCLUSION

 Implementing tailored prophylaxis using a Bayesian forecasting approach in a routine clinical practice setting may improve haemophilia clinical outcomes.

Physiologically-Based Pharmacokinetic model for Ciprofloxacin in children with complicated Urinary Tract Infection

In a recent multicenter population pharmacokinetic study of ciprofloxacin administered to children suffering from complicated urinary tract infection (cUTI), the apparent volume of distribution (V) and total plasma clearance (CL) were decreased by 83.6% and 41.5% respectively, compared to healthy children. To understand these differences, a physiologically-based pharmacokinetic model (PBPK) for ciprofloxacin was developed for cUTI children. First, a PBPK model in adults was developed, modified incorporating age-dependent functions and evaluated with paediatric data generated from a published model in healthy children. Then, the model was then adapted to a cUTI paediatric population according to the degree of renal impairment (KF) affecting renal clearance (CL_Renal,) and CYP1A2 clearance (CL_CYP1A2). Serum and urine samples obtained from 22 cUTI children were used for model evaluation. Lastly, a parameter sensitivity analysis identified the most influential parameters on V and CL. The PBPK model predicted the ciprofloxacin exposure in adults and children, capturing age-related pharmacokinetic changes. Plasma concentrations and fraction excreted unchanged in urine (f_e) predictions improved in paediatric cUTI patients once CL_renal and CL_CYP1A2 were corrected by KF. The presented PBPK model for ciprofloxacin demonstrates its adequacy to simulate different dosing scenarios to obtain PK predictions in a healthy population from 3 months old onwards. Model adaptation of CL_Renal and CL_CYP1A2 according to KF explained partially the differences seen in the plasma drug concentrations and f_e vs time profiles between healthy and cUTI children. Nevertheless, it is necessary to further investigate the disease-related changes in cUTI to improve model predictions.

Physiologically Based Pharmacokinetic Approach to Determine Dosing on Extracorporeal Life Support: Fluconazole in Children on ECMO

Extracorporeal life support (e.g., dialysis, extracorporeal membrane oxygenation (ECMO)) can affect drug disposition, placing patients at risk for therapeutic failure. In this population, dose selection to achieve safe and effective drug exposure is difficult. We developed a novel and flexible approach that uses physiologically based pharmacokinetic (PBPK) modeling to translate results from ECMO ex vivo experiments into bedside dosing recommendations. To determine fluconazole dosing in children on ECMO, we developed a PBPK model, which was validated using fluconazole pharmacokinetic (PK) data in adults and critically ill infants. Next, an ECMO compartment was added to the PBPK model and parameterized using data from a previously published ex vivo study. Simulations using the final ECMO PBPK model reasonably characterized observed PK data in infants on ECMO, and the model was used to derive dosing in children on ECMO across the pediatric age spectrum. This approach can be generalized to other forms of extracorporeal life support (ECLS), such as dialysis.

Performing and interpreting individual pharmacokinetic profiles in patients with Hemophilia A or B: Rationale and general considerations

OBJECTIVES

In a separate document, we have provided specific guidance on performing individual pharmacokinetic (PK) studies using limited samples in persons with hemophilia with the goal to optimize prophylaxis with clotting factor concentrates. This paper, intended for clinicians, aims to describe how to interpret and apply PK properties obtained in persons with hemophilia.

METHODS

The members of the Working Party on population PK (PopPK) of the ISTH SSC Subcommittee on Factor VIII and IX and rare bleeding disorders, together with additional hemophilia and PK experts, completed a survey and ranking exercise whereby key areas of interest in the field were identified. The group had regular web conferences to refine the manuscript's scope and structure, taking into account comments from the external feedback to the earlier document.

RESULTS

Many clinical decisions in hemophilia are based on some form of explicit or implicit PK assessment. Individual patient PK profiles can be analyzed through traditional or PopPK methods, with the latter providing the advantage of fewer samples needing to be collected on any prophylaxis regimen, and without the need the for a washout period. The most useful presentation of PK results for clinical decision making are a curve of the factor activity level over time, the time to achieve a certain activity level, or related parameters like half-life or exposure (AUC). Software platforms have been developed to deliver this information to clinicians at the point of care. Key characteristics of studies measuring average PK parameters were reviewed, outlining what makes a credible head-to-head comparison among different concentrates. Large data collections of PK and treatment outcomes currently ongoing will advance care in the future.

CONCLUSIONS

Traditionally used to compare different concentrates, PK can support tailoring of hemophilia treatment by individual profiling, which is greatly simplified by adopting a PopPK/Bayesian method and limited sampling protocol.

Improving Pediatric Protein Binding Estimates: An Evaluation of α1-Acid Glycoprotein Maturation in Healthy and Infected Subjects

BACKGROUND

Differences in plasma protein levels observed between children and adults can alter the extent of xenobiotic binding in plasma, resulting in divergent patterns of exposure.

OBJECTIVE

This study aims to quantify the ontogeny of α1-acid glycoprotein in both healthy and infected subjects.

METHODS

Data pertaining to α1-acid glycoprotein from healthy subjects were compiled over 26 different publications. For subjects diagnosed or suspected of infection, α1-acid glycoprotein levels were obtained from 214 individuals acquired over three clinical investigations. The analysis evaluated the use of linear, power, exponential, log-linear, and sigmoid E max models to describe the ontogeny of α1-acid glycoprotein. Utility of the derived ontogeny equation for estimation of pediatric fraction unbound was evaluated using average-fold error and absolute average-fold error as measures of bias and precision, respectively. A comparison to fraction unbound estimates derived using a previously proposed linear equation was also instituted.

RESULTS

The sigmoid E max model provided the comparatively best depiction of α1-acid glycoprotein ontogeny in both healthy and infected subjects. Despite median α1-acid glycoprotein levels in infected subjects being more than two-fold greater than those observed in healthy subjects, a similar ontogeny pattern was observed when levels were normalized toward adult levels. For estimation of pediatric fraction unbound, the α1-acid glycoprotein ontogeny equation derived from this work (average fold error 0.99; absolute average fold error 1.24) provided a superior predictive performance in comparison to the previous equation (average fold error 0.74; absolute average fold error 1.45).

CONCLUSION

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

Development of a Pediatric Physiologically-Based Pharmacokinetic Model of Clindamycin Using Opportunistic Pharmacokinetic Data

Physiologically-based pharmacokinetic (PBPK) modeling is a powerful tool used to characterize maturational changes in drug disposition to inform dosing across childhood; however, its use is limited in pediatric drug development. Access to pediatric pharmacokinetic data is a barrier to widespread application of this model, which impedes its development and optimization. To support the development of a pediatric PBPK model, we sought to leverage opportunistically-collected plasma concentrations of the commonly used antibiotic clindamycin. The pediatric PBPK model was optimized following development of an adult PBPK model that adequately described literature data. We evaluated the predictability of the pediatric population PBPK model across four age groups and found that 63-93% of the observed data were captured within the 90% prediction interval of the model. We then used the pediatric PBPK model to optimize intravenous clindamycin dosing for a future prospective validation trial. The optimal dosing proposed by this model was 9 mg/kg/dose in children ≤5 months of age, 12 mg/kg/dose in children >5 months-6 years of age, and 10 mg/kg/dose in children 6-18 years of age, all administered every 8 h. The simulated exposures achieved with the dosing regimen proposed were comparable with adult plasma and tissue exposures for the treatment of community-acquired methicillin-resistant Staphylococcus aureus infections. Our model demonstrated the feasibility of using opportunistic pediatric data to develop pediatric PBPK models, extending the reach of this powerful modeling tool and potentially transforming the pediatric drug development field.

Modeling of Body Weight Metrics for Effective and Cost-Efficient Conventional Factor VIII Dosing in Hemophilia A Prophylaxis

The total body weight-based dosing strategy currently used in the prophylactic treatment of hemophilia A may not be appropriate for all populations. The assumptions that guide weight-based dosing are not valid in overweight and obese populations, resulting in overdosing and ineffective resource utilization. We explored different weight metrics including lean body weight, ideal body weight, and adjusted body weight to determine an alternative dosing strategy that is both safe and resource-efficient in normal and overweight/obese adult patients. Using a validated population pharmacokinetic model, we simulated a variety of dosing regimens using different doses, weight metrics, and frequencies; we also investigated the implications of assuming various levels of endogenous factor production. Ideal body weight performed the best across all of the regimens explored, maintaining safety while moderating resource consumption for overweight and obese patients.

Point/Counterpoint: Are Outstanding Leaders Born or Made?

The question of whether outstanding leaders are born or made has been debated for years. There are numerous examples of historical figures that came naturally to leadership, while others developed their leadership skills through tenacity and experience. To understand leadership, both nature (the genetic component) and nurture (the environmental influences) must be considered. This article represents the work of two Academic Leadership Fellows Program groups who debated each position at the 2016 American Association of Colleges of Pharmacy (AACP) Interim Meeting in Tampa, Fla., in February 2016.

Population PBPK modelling of trastuzumab: a framework for quantifying and predicting inter-individual variability

In this work we proposed a population physiologically-based pharmacokinetic (popPBPK) framework for quantifying and predicting inter-individual pharmacokinetic variability using the anti-HER2 monoclonal antibody (mAb) trastuzumab as an example. First, a PBPK model was developed to account for the possible mechanistic sources of variability. Within the model, five key factors that contribute to variability were identified and the nature of their contribution was quantified with local and global sensitivity analyses. The five key factors were the concentration of membrane-bound HER2 ([Formula: see text]), the convective flow rate of mAb through vascular pores ([Formula: see text]), the endocytic transport rate of mAb through vascular endothelium ([Formula: see text]), the degradation rate of mAb-HER2 complexes ([Formula: see text]) and the concentration of shed HER2 extracellular domain in circulation ([Formula: see text]). [Formula: see text] was the most important parameter governing trastuzumab distribution into tissues and primarily affected variability in the first 500 h post-administration. [Formula: see text] was the most significant contributor to variability in clearance. These findings were used together with population generation methods to accurately predict the observed variability in four experimental trials with trastuzumab. To explore anthropometric sources of variability, virtual populations were created to represent participants in the four experimental trials. Using populations with only their expected anthropometric diversity resulted in under-prediction of the observed inter-individual variability. Adapting the populations to include literature-based variability around the five key parameters enabled accurate predictions of the variability in the four trials. The successful application of this framework demonstrates the utility of popPBPK methods to understand the mechanistic underpinnings of pharmacokinetic variability.