Generic Workflow to Predict Medicine Concentrations in Human Milk Using Physiologically-Based Pharmacokinetic (PBPK) Modelling-A Contribution from the ConcePTION Project.

Women commonly take medication during lactation. Currently, there is little information about the exposure-related safety of maternal medicines for breastfed infants. The aim was to explore the performance of a generic physiologically-based pharmacokinetic (PBPK) model to predict concentrations in human milk for ten physiochemically diverse medicines. First, PBPK models were developed for "non-lactating" adult individuals in PK-Sim/MoBi v9.1 (Open Systems Pharmacology). The PBPK models predicted the area-under-the-curve (AUC) and maximum concentrations (Cmax) in plasma within a two-fold error. Next, the PBPK models were extended to include lactation physiology. Plasma and human milk concentrations were simulated for a three-months postpartum population, and the corresponding AUC-based milk-to-plasma (M/P) ratios and relative infant doses were calculated. The lactation PBPK models resulted in reasonable predictions for eight medicines, while an overprediction of human milk concentrations and M/P ratios (>2-fold) was observed for two medicines. From a safety perspective, none of the models resulted in underpredictions of observed human milk concentrations. The present effort resulted in a generic workflow to predict medicine concentrations in human milk. This generic PBPK model represents an important step towards an evidence-based safety assessment of maternal medication during lactation, applicable in an early drug development stage.

Current Practices, Gap Analysis, and Proposed Workflows for PBPK Modeling of Cytochrome P450 Induction: An Industry Perspective.

The International Consortium for Innovation and Quality (IQ) Physiologically Based Pharmacokinetic (PBPK) Modeling Induction Working Group (IWG) conducted a survey across participating companies around general strategies for PBPK modeling of induction, including experience with its utility to address various questions, regulatory interactions, and regulatory acceptance. The results highlight areas where PBPK modeling is used with high confidence and identifies opportunities where confidence is lower and further evaluation is needed. To enhance the survey results, the PBPK-IWG also collected case studies and analyzed recent literature examples where PBPK models were applied to predict CYP3A induction-mediated drug-drug interactions. PBPK modeling of induction has evolved and progressed significantly, proving to have great potential to accelerate drug discovery and development. With the aim of enabling optimal use for new molecular entities that are either substrates and/or inducers of CYP3A, the PBPK-IWG proposes initial workflows for PBPK application, discusses future trends, and identifies gaps that need to be addressed.

Characterization of Hepatic UDP-Glucuronosyltransferase Enzyme Abundance-Activity Correlations and Population Variability Using a Proteomics Approach and Comparison with Cytochrome P450 Enzymes.

The expression of ten major drug-metabolizing UDP-glucuronosyltransferase (UGT) enzymes in a panel of 130 human hepatic microsomal samples was measured using a liquid chromatography-tandem mass spectrometry-based approach. Simultaneously, ten cytochromes P450 and P450 reductase were also measured, and activity-expression relationships were assessed for comparison. The resulting data sets demonstrated that, with the exception of UGT2B17, 10th to 90th percentiles of UGT expression spanned 3- to 8-fold ranges. These ranges were small relative to ranges of reported mean UGT enzyme expression across different laboratories. We tested correlation of UGT expression with enzymatic activities using selective probe substrates. A high degree of abundance-activity correlation (Spearman's rank correlation coefficient > 0.6) was observed for UGT1As (1A1, 3, 4, 6) and cytochromes P450. In contrast, protein abundance and activity did not correlate strongly for UGT1A9 and UGT2B enzymes (2B4, 7, 10, 15, and 17). Protein abundance was strongly correlated for UGTs 2B7, 2B10, and 2B15. We suggest a number of factors may contribute to these differences including incomplete selectivity of probe substrates, correlated expression of these UGT2B isoforms, and the impact of splice and polymorphic variants on the peptides used in proteomics analysis, and exemplify this in the case of UGT2B10. Extensive correlation analyses identified important criteria for validating the fidelity of proteomics and enzymatic activity approaches for assessing UGT variability, population differences, and ontogenetic changes. SIGNIFICANCE STATEMENT: Protein expression data allow detailed assessment of interindividual variability and enzyme ontogeny. This study has observed that expression and enzyme activity are well correlated for hepatic UGT1A enzymes and cytochromes P450. However, for the UGT2B family, caution is advised when assuming correlation of expression and activity as is often done in physiologically based pharmacokinetic modeling. This can be due to incomplete probe substrate specificities, but may also be related to presence of inactive UGT protein materials and the effect of splicing variations.

Pediatric Pharmacokinetics and Dose Predictions: A Report of a Satellite Meeting to the 10th Juvenile Toxicity Symposium.

On April 24, 2019, a symposium on Pediatric Pharmacokinetics and Dose Predictions was held as a satellite meeting to the 10th Juvenile Toxicity Symposium. This symposium brought together scientists from academia, industry, and clinical research organizations with the aim to update each other on the current knowledge on pediatric drug development. Through more knowledge on specific ontogeny profiles of drug metabolism and transporter proteins, integrated into physiologically-based pharmacokinetic (PBPK) models, we have gained a more integrated understanding of age-related differences in pharmacokinetics (PKs), Relevant examples were presented during the meeting. PBPK may be considered the gold standard for pediatric PK prediction, but still it is important to know that simpler methods, such as allometry, allometry combined with maturation function, functions based on the elimination pathway, or linear models, also perform well, depending on the age range or the mechanisms involved. Knowledge from different methods and information sources should be combined (e.g., microdosing can reveal early read-out of age-related differences in exposure), and such results can be a value to verify models. To further establish best practices for dose setting in pediatrics, more in vitro and in vivo research is needed on aspects such as age-related changes in the exposure-response relationship and the impact of disease on PK. New information coupled with the refining of model-based drug development approaches will allow faster targeting of intended age groups and allow more efficient design of pediatric clinical trials.

Coexpression of Human Hepatic Uridine Diphosphate Glucuronosyltransferase Proteins: Implications for Ontogenetic Mechanisms and Isoform Coregulation.

Uridine diphosphate glucuronosyltransferases (UGTs) catalyze glucuronidation to facilitate systemic and local clearance of numerous chemicals and drugs. To investigate whether UGT expression is coregulated in human liver, we analyzed the protein expression of UGTs 1A1, 1A3, 1A4, 1A6, 1A9, 2B7, 3A1, and 3A2 using western blots from 164 healthy human liver samples, comparing expression with age and sex. UGT1A6 levels were significantly higher in children than adults, and UGT3A1 and 3A2 expression significantly increased with age from childhood to age >65 yearas. In children aged <18 years, UGT1A4/1A9 protein expression was significantly correlated, but not for adults aged >18 years. UGT1A3 expression was always significantly correlated with other UGT1A isoforms in all adults aged >18 years. In individuals aged ≥12 years, expression of UGT1A1/1A4, UGT1A1/1A6, UGT1A1/1A9, and UGT1A4/1A6 significantly correlated, which was not observed in children aged <12 years. In contrast, UGT1A4/2B7 showed significant correlation in children aged <12 years, but not in individuals aged ≥12 years, and this was observed in female but not male individuals. Expression of UGT1A6/1A9 and UGT3A1/3A2 correlated in the entire sample population, but UGT3As did not correlate with other UGTs. These correlations were sex dependent, as UGT1A3/1A1, UGT1A4/2B7 and UGT3A1/3A2 correlated more highly in male than female individuals, while UGT1A4/1A6 protein correlated more significantly in female than male individuals. This is the first report on the ontogeny of UGT3A isoforms, showing maximal expression in the elderly, and is the first demonstration that UGT isoforms commonly coexpress in vivo, in both age-dependent and sex-dependent manners.

Use of Phenotypically Poor Metabolizer Individual Donor Human Liver Microsomes To Identify Selective Substrates of UGT2B10.

UDP-glucuronosyltransferase (UGT)1A4 and UGT2B10 are the human UGT isoforms most frequently involved in N-glucuronidation of drugs. UGT2B10 exhibits higher affinity than UGT1A4 for numerous substrates, making it potentially the more important enzyme for metabolism of these compounds in vivo. Clinically relevant UGT2B10 polymorphisms, including a null activity splice site mutation common in African populations, can lead to large exposure differences for UGT2B10 substrates that may limit their developability as marketed drugs. UGT phenotyping approaches using recombinantly expressed UGTs are limited by low enzyme activity and lack of validation of scaling to in vivo. In this study, we describe the use of an efficient experimental protocol for identification of UGT2B10-selective substrates (i.e., those with high fraction metabolized by UGT2B10), which exploits the activity difference between pooled human liver microsomes (HLM) and HLM from a phenotypically UGT2B10 poor metabolizer donor. Following characterization of the approach with eight known UGT2B10 substrates, we used ligand-based virtual screening and literature precedents to select 24 potential UGT2B10 substrates of 140 UGT-metabolized drugs for testing. Of these, dothiepin, cidoxepin, cyclobenzaprine, azatadine, cyproheptadine, bifonazole, and asenapine were indicated to be selective UGT2B10 substrates that have not previously been described. UGT phenotyping experiments and tests comparing conjugative and oxidative clearance were then used to confirm these findings. These approaches provide rapid and sensitive ways to evaluate whether a potential drug candidate cleared via glucuronidation will be sensitive to UGT2B10 polymorphisms in vivo. SIGNIFICANCE STATEMENT: The role of highly polymorphic UDP-glucuronosyltransferase (UGT)2B10 is likely to be underestimated currently for many compounds cleared via N-glucuronidation due to high test concentrations often used in vitro and low activity of UGT2B10 preparations. The methodology described in this study can be combined with the assessment of UGT versus oxidative in vitro metabolism to rapidly identify compounds likely to be sensitive to UGT2B10 polymorphism (high fraction metabolized by UGT2B10), enabling either chemical modification or polymorphism risk assessment before candidate selection.

Characterization of the Ontogeny of Hepatic UDP-Glucuronosyltransferase Enzymes Based on Glucuronidation Activity Measured in Human Liver Microsomes.

An understanding of the postnatal development of hepatic UDP-glucuronosyltransferase (UGT) enzymes is required for accurate prediction of the age-dependent changes in pharmacokinetics of many drugs used in children. However, the maturation rate of hepatic UGT isoforms remains a major knowledge gap. This study aimed to establish the age-associated changes in glucuronidation activity of 10 major hepatic UGT isoforms in humans, namely, UGT1A1, UGT1A3, UGT1A4, UGT1A6, UGT1A9, UGT2B4, UGT2B7, UGT2B10, UGT2B15, and UGT2B17. Human liver microsomes from pediatric and adult donors were incubated under optimized incubation conditions to assess the activity rates of hepatic UGT isoforms using a panel of 19 in vitro UGT probe substrates and clinically used drugs. Statistically strong correlations of glucuronidation activities allowed the ontogeny of UGT1A1, UGT1A4, UGT2B7, UGT2B10, and UGT2B15 to be established using multiple selective UGT substrates and matched human liver microsome samples. The postnatal development of hepatic UGTs is isoform-dependent using either individual or cross-correlated selective isoform substrates. Maximal adult activity was reached at different times ranging from within a month (UGT1A1, UGT2B4, UGT2B7, UGT2B10, and UGT2B15), during infancy (UGT1A3, UGT1A4, and UGT1A9), to adolescence (UGT1A6 and UGT2B17). This study provides an extensive characterization of the postnatal ontogeny profiles of hepatic UGT enzymes that are instrumental for predicting drug disposition via in vitro-in vivo extrapolation algorithms and verifying pharmacokinetic predictions against in vivo observations via pediatric physiologically based pharmacokinetic modeling in pediatric patients.

The Ontogeny of UDP-glucuronosyltransferase Enzymes, Recommendations for Future Profiling Studies and Application Through Physiologically Based Pharmacokinetic Modelling.

Limited understanding of drug pharmacokinetics in children is one of the major challenges in paediatric drug development. This is most critical in neonates and infants owing to rapid changes in physiological functions, especially in the activity of drug-metabolising enzymes. Paediatric physiologically based pharmacokinetic models that integrate ontogeny functions for cytochrome P450 enzymes have aided our understanding of drug exposure in children, including those under the age of 2 years. Paediatric physiologically based pharmacokinetic models have consequently been recognised by the European Medicines Agency and the US Food and Drug Administration as innovative tools in paediatric drug development and regulatory decision making. However, little is currently known about age-related changes in UDP-glucuronosyltransferase-mediated metabolism, which represents the most important conjugation reaction for xenobiotics. Therefore, the objective of the review was to conduct a thorough literature survey to summarise our current understanding of age-related changes in UDP-glucuronosyltransferases as well as associated clinical and experimental sources of variance. Our findings indicate that there are distinct differences in UDP-glucuronosyltransferase expression and activity between isoforms for different age groups. In addition, there is substantial variability between individuals and laboratories reported for human liver microsomes, which results in part from a lack of standardised experimental conditions. Therefore, we provide a number of best practice recommendations for experimental conditions, which ultimately may help improve the quality of data used for quantitative clinical pharmacology approaches, and thus for safe and effective pharmacotherapy in children.

Optimization of Experimental Conditions of Automated Glucuronidation Assays in Human Liver Microsomes Using a Cocktail Approach and Ultra-High Performance Liquid Chromatography-Tandem Mass Spectrometry.

UDP-glucuronosyltransferase (UGT)-mediated metabolism is possibly the most important conjugation reaction for marketed drugs. However, there are currently no generally accepted standard incubation conditions for UGT microsomal assays, and substantial differences in experimental design and methodology between laboratories hinder cross-study comparison of in vitro activities. This study aimed to define optimal experimental conditions to determine glucuronidation activity of multiple UGT isoforms simultaneously using human liver microsomes. Hepatic glucuronidation activities of UGT1A1, UGT1A3, UGT1A4, UGT1A6, UGT1A9, UGT2B4, UGT2B7, UGT2B10, UGT2B15, and UGT2B17 were determined using cocktail incubations of 10 UGT probe substrates. Buffer components and cosubstrates were assessed over a range of concentrations including magnesium chloride (MgCl2; 0-10 mM) and uridine 5'-diphosphoglucuronic acid (UDPGA; 1-25 mM) with either Tris-HCl or potassium phosphate buffer (100 mM, pH 7.4). Greater microsomal glucuronidation activity by different hepatic UGT isoforms was obtained using 10 mM MgCl2 and 5 mM UDPGA with 100 mM Tris-HCl buffer. The influence of bovine serum albumin (BSA; 0.1%-2% w/v) on glucuronidation activity was also assessed. Enzyme- and substrate-dependent effects of BSA were observed, resulting in decreased total activity of UGT1A1, UGT1A3, and UGT2B17 and increased total UGT1A9 and UGT2B7 activity. The inclusion of BSA did not significantly reduce the between-subject variability of UGT activity. Future in vitro UGT profiling studies under the proposed optimized experimental conditions would allow high-quality positive control data to be generated across laboratories, with effective control of a high degree of between-donor variability for UGT activity and for chemical optimization toward lower-clearance drug molecules in a pharmaceutical drug discovery setting.

Meta-analysis of expression of hepatic organic anion-transporting polypeptide (OATP) transporters in cellular systems relative to human liver tissue.

Organic anion-transporting polypeptide (OATP)1B1, OATP1B3, and OATP2B1 transporters play an important role in hepatic drug disposition. Recently, an increasing number of studies have reported proteomic expression data for OATP transporters. However, systematic analysis and understanding of the actual differences in OATP expression between liver tissue and commonly used cellular systems is lacking. In the current study, meta-analysis was performed to assess the protein expression of OATP transporters reported in hepatocytes relative to liver tissue and to identify any potential correlations in transporter expression levels in the same individual. OATP1B1 was identified as the most abundant uptake transporter at 5.9 ± 8.3, 5.8 ± 3.3, and 4.2 ± 1.7 fmol/µg protein in liver tissue, sandwich-cultured human hepatocytes (SCHH), and cryopreserved suspended hepatocytes, respectively. The rank order in average expression in liver tissue and cellular systems was OATP1B1 > OATP1B3 ≈ OATP2B1. Abundance levels of the OATP transporters investigated were not significantly different between liver and cellular systems, with the exception of OATP2B1 expression in SCHH relative to liver tissue. Analysis of OATP1B1, OATP1B3, and OATP2B1 liver expression data in the same individuals (n = 86) identified weak (OATP1B1-OATP2B1) to moderately (OATP1B3-OATP2B1) significant correlations. A significant weak correlation was noted between OATP1B1 abundance and age of human donors, whereas expression of the OATPs investigated was independent of sex. Implications of the current analysis on the in vitro-in vivo extrapolation of transporter-mediated drug disposition using physiologically based pharmacokinetic models are discussed.