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

Age-Specific ADME Gene Expression in Infant Intestinal Enteroids

In childhood, developmental changes and environmental interactions highly affect orally dosed drug disposition across the age range. To optimize dosing regimens and ensure safe use of drugs in pediatric patients, understanding this age-dependent biology is necessary. In this proof-of-concept study, we aimed to culture age-specific enteroids from infant tissue which represent its original donor material, specifically for drug transport and metabolism. Enteroid lines from fresh infant tissues (n = 8, age range: 0.3-45 postnatal weeks) and adult tissues (n = 3) were established and expanded to 3D self-organizing enteroids. The gene expression of drug transporters P-gp (ABCB1), BCRP (ABCG2), MRP2 (ABCC2), and PEPT1 (SLC15A1) and drug metabolizing enzymes CYP3A4, CYP2C18, and UGT1A1 was determined with RT-qPCR in fresh tissue and its derivative differentiated enteroids. Expression levels of P-gp, BCRP, MRP2, and CYP3A4 were similar between tissues and enteroids. PEPT1 and CYP2C18 expression was lower in enteroids compared to that in the tissue. The expression of UGT1A1 in the tissue was lower than that in enteroids. The gene expression did not change with the enteroid passage number for all genes studied. Similar maturational patterns in tissues and enteroids were visually observed for P-gp, PEPT1, MRP2, CYP3A4, CYP2C18, and VIL1. In this explorative study, interpatient variability was high, likely due to the diverse patient characteristics of the sampled population (e.g., disease, age, and treatment). To summarize, maturational patterns of clinically relevant ADME genes in tissue were maintained in enteroids. These findings are an important step toward the potential use of pediatric enteroids in pediatric drug development, which in the future may lead to improved pediatric safety predictions during drug development. We reason that such an approach can contribute to a potential age-specific platform to study and predict drug exposure and intestinal safety in pediatrics.

Retrospective Data Analysis Reveals Unusual Metabolism Pattern of Ethanol in Pediatrics as Compared to Adult and Geriatric Populations

BACKGROUND

About 95% of consumed ethanol is metabolized by oxidative pathways. Less than 1% is metabolized via nonoxidative pathways: glucuronidation, sulfation, and the formation of fatty acid esters of ethanol. In neonates, the glucuronidation pathway has been reported to be underdeveloped but matures with age. This work compared the test results of patients' random urine samples submitted to our facility for ethyl glucuronide (EtG) and ethyl sulfate (EtS) measurements across pediatric and adult populations.

METHODS

Test results (n = 63 498) from urine samples tested for EtG and EtS by quantitative liquid chromatography-tandem mass spectrometry at our facility were utilized for this study. EtG and EtS concentrations were compared across the age partitions 0 to 17 years (pediatric), 18 to 80 years (adult), and 81 to 100 years (geriatric). Eight pediatric patients from a tertiary academic hospital contributed clinical context via abstracted clinical information.

RESULTS

Across the individual age partitions, 60% to 65% of patients had both EtG and EtS present in urine. Approximately 5% to 10% of patients had only EtG, and 25% to 35% had neither metabolite present. The lowest percentages (<1.5%) had EtS present in the absence of EtG. Markedly, no pediatric patients had only EtS present; compared to the adult population, this was statistically significant (Fisher exact test, P = 0.025).

CONCLUSIONS

From the data presented in this work, EtG is more prevalent relative to EtS in urine samples of patients assessed for ethanol exposure.

Update of the scientific opinion on tetrabromobisphenol A (TBBPA) and its derivatives in food

The European Commission asked EFSA to update its 2011 risk assessment on tetrabromobisphenol A (TBBPA) and five derivatives in food. Neurotoxicity and carcinogenicity were considered as the critical effects of TBBPA in rodent studies. The available evidence indicates that the carcinogenicity of TBBPA occurs via non-genotoxic mechanisms. Taking into account the new data, the CONTAM Panel considered it appropriate to set a tolerable daily intake (TDI). Based on decreased interest in social interaction in male mice, a lowest observed adverse effect level (LOAEL) of 0.2 mg/kg body weight (bw) per day was identified and selected as the reference point for the risk characterisation. Applying the default uncertainty factor of 100 for inter- and intraspecies variability, and a factor of 3 to extrapolate from the LOAEL to NOAEL, a TDI for TBBPA of 0.7 μg/kg bw per day was established. Around 2100 analytical results for TBBPA in food were used to estimate dietary exposure for the European population. The most important contributors to the chronic dietary LB exposure to TBBPA were fish and seafood, meat and meat products and milk and dairy products. The exposure estimates to TBBPA were all below the TDI, including those estimated for breastfed and formula-fed infants. Accounting for the uncertainties affecting the assessment, the CONTAM Panel concluded with 90%-95% certainty that the current dietary exposure to TBBPA does not raise a health concern for any of the population groups considered. There were insufficient data on the toxicity of any of the TBBPA derivatives to derive reference points, or to allow a comparison with TBBPA that would support assignment to an assessment group for the purposes of combined risk assessment.

Population pharmacokinetics of oxcarbazepine active metabolite in Chinese paediatric patients with epilepsy: Model-based dose optimization

The pharmacological activity of oxcarbazepine (OXC) is primarily exerted through its active 10-monohydroxy metabolite (MHD). Nonetheless, there is limited pharmacokinetic information available regarding paediatric patients with epilepsy treated with OXC, especially in infants and toddlers. Concurrently, this drug exhibits substantial variability in pharmacokinetics and therapeutic response across different individuals. We aimed to develop a model to quantitatively investigate factors that affect MHD pharmacokinetics to formulate a dosage guideline for OXC in Chinese paediatric patients. A total of 297 MHD trough concentrations were obtained from 287 epileptic children. Six body weight (BW)-based allometric models were used for population pharmacokinetic modelling, while investigating the impact of other covariates on the apparent clearance. The one-compartment model and age cut-off model for the apparent clearance (CL/F) were established to describe the pharmacokinetics of MHD. The probability to obtain target trough concentration ranges (TTCRs) of MHD between 3 and 35 mg/L was determined by Monte Carlo simulations for doses ranging from 8 to 90 mg/kg/day. A new dose optimization strategy combining the dosage guidelines and Bayesian method provides a tailored approach for Chinese paediatric epileptic patients based on their individual BW and desired TTCRs of MHD, and also supports current dose recommendations, with the exception of children weighing ≤5 kg.

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

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

Significant Effects of Renal Function on Mycophenolic Acid Total Clearance in Pediatric Kidney Transplant Recipients with Population Pharmacokinetic Modeling

BACKGROUND AND OBJECTIVES

Mycophenolic acid (MPA) is an immunosuppressant commonly prescribed in pediatric kidney transplantation to prevent graft rejection. Large variabilities in MPA plasma exposures have been observed in this population, which could result in severe adverse effects. The majority of the MPA pharmacokinetic data have been reported in adult populations, whereas information in pediatric patients is still very limited. The objective of this study was to establish a novel, nonlinear mixed-effects model for MPA and investigate the clinical variables affecting MPA population pharmacokinetics in pediatric kidney transplant recipients.

METHODS

Data were collected retrospectively from pediatric kidney transplant patients (≤ 18 years when MPA concentrations were initially collected; on oral administration of mycophenolate mofetil) in Calgary, Alberta, Canada. Nonlinear mixed-effect modeling was conducted using stochastic approximation expectation-maximization in Monolix 2021R2 (Lixoft SAS, France) to determine population pharmacokinetic estimates, interindividual variabilities, and interoccasional variabilities. Covariate models were constructed using the Model Proposal function in Monolix in conjunction with a systematic stepwise inclusion/elimination protocol. The best model was selected based on objective function values, relative standard errors, goodness-of-fit plots, prediction-corrected visual predictive checks, and numerical predictive checks.

RESULTS

A total of 50 pediatric kidney transplant patients (25 female) with 219 MPA plasma concentration-time profiles were included. The average age (± standard deviation) and posttransplant time for the sample population were 12.8 ± 4.8 years and 762 ± 1160 days, respectively. The majority of study subjects (i.e., > 85% based on all occasions) were co-administered tacrolimus. A two-compartment, first-order absorption with lag time and linear elimination structural model with lognormal distributed proportional residual errors best described the MPA concentration-time data. The absorption rate constant (2.52 h-1 or 0.042 min-1), lag time (0.166 h or 9.96 min), volumes of distributions of the central (22.8 L) and peripheral (216 L) compartments, and intercompartment clearance (17.6 L h-1 or 0.293 L min-1) were consistent with literature values; whereas total MPA clearance (0.72 L h-1 or 0.012 L min-1) was relatively reduced, likely due to the general lack of cyclosporine interactions and the stabilized graft functions from significantly longer posttransplant time in our sample population. Of the clinical variables tested, only estimated glomerular filtration rate (eGFR) was identified a significant covariate affecting total MPA clearance with a positive, exponential relationship. The final population pharmacokinetic model was successfully evaluated/validated using a variety of complementary methods.

CONCLUSION

We have successfully constructed and validated a novel population pharmacokinetic model of MPA in pediatric kidney transplant patients. A positive, nonlinear relationship between eGFR and total MPA clearance identified in our model is likely attributed to multiple concurrent mechanisms, which warrant further systematic investigations.

Nervonic acid improves liver inflammation in a mouse model of Parkinson's disease by inhibiting proinflammatory signaling pathways and regulating metabolic pathways

BACKGROUND

Nervonic acid (NA) - a type of bioactive fatty acid that is found in natural sources - can inhibit inflammatory reactions and regulate immune system balance. Therefore, the use of NA for the treatment of neurodegenerative diseases has received considerable attention. Our previous study found that NA inhibited inflammatory responses in the brain of Parkinson's disease (PD) mouse models. In addition to the brain, PD is also associated with visceral organ dysfunction, especially impaired liver function. Thus, studying the role of NA in PD-mediated inflammation of the liver is particularly important.

METHODS

A combined transcriptome and metabolomic approach was utilized to investigate the anti-inflammatory effects of NA on the liver of PD mice. Inflammatory signaling molecules and metabolic pathway-related genes were examined in the liver using real-time PCR and western blotting.

RESULTS

Liver transcriptome analysis revealed that NA exerted anti-inflammatory effects by controlling several pro-inflammatory signaling pathways, such as the down-regulation of the tumor necrosis factor and nuclear factor kappa B signaling pathways, both of which were essential in the development of inflammatory disease. In addition, liver metabolomic results revealed that metabolites related to steroid hormone biosynthesis, arachidonic acid metabolism, and linoleic acid metabolism were up-regulated and those related to valine, leucine, and isoleucine degradation pathways were down-regulated in NA treatment groups compared with the PD model. The integration of metabolomic and transcriptomic results showed NA significantly exerted its anti-inflammatory function by regulating the transcription and metabolic pathways of multiple genes. Particularly, linoleic acid metabolism, arachidonic acid metabolism, and steroid hormone biosynthesis were the crucial pathways of the anti-inflammatory action of NA. Key genes in these metabolic pathways and key molecules in inflammatory signaling pathways were also verified, which were consistent with transcriptomic results.

CONCLUSION

These findings provide novel insights into the liver protective effects of NA against PD mice. This study also showed that NA could be a useful dietary element for improving and treating PD-induced liver inflammation.

Contribution of UGT Enzymes to Human Drug Metabolism Stereoselectivity: A Case Study of Medetomidine, RO5263397, Propranolol, and Testosterone

The enantiomeric forms of chiral compounds have identical physical properties but may vary greatly in their metabolism by individual enzymes. Enantioselectivity in UDP-glucuronosyl transferase (UGT) metabolism has been reported for a number of compounds and with different UGT isoforms involved. However, the impact of such individual enzyme results on overall clearance stereoselectivity is often not clear. The enantiomers of medetomidine, RO5263397, and propranolol and the epimers testosterone and epitestosterone exhibit more than a 10-fold difference in glucuronidation rates by individual UGT enzymes. In this study, we examined the translation of human UGT stereoselectivity to hepatic drug clearance considering the combination of multiple UGTs to overall glucuronidation, the contribution of other metabolic enzymes such as cytochrome P450s (P450s), and the potential for differences in protein binding and blood/plasma partitioning. For medetomidine and RO5263397, the high individual enzyme (UGT2B10) enantioselectivity translated into ∼3- to >10-fold differences in predicted human hepatic in vivo clearance. For propranolol, the UGT enantioselectivity was irrelevant in the context of high P450 metabolism. For testosterone, a complex picture emerged due to differential epimeric selectivity of various contributing enzymes and potential for extrahepatic metabolism. Quite different patterns of P450- and UGT-mediated metabolism were observed across species, as well as differences in stereoselectivity, indicating that extrapolation from human enzyme and tissue data are essential when predicting human clearance enantioselectivity. SIGNIFICANCE STATEMENT: Individual enzyme stereoselectivity illustrates the importance of three-dimensional drug-metabolizing enzyme-substrate interactions and is essential when considering the clearance of racemic drugs. However, translation from in vitro to in vivo can be challenging as contributions from multiple enzymes and enzyme classes must be combined with protein binding and blood/plasma partitioning data to estimate the net intrinsic clearance for each enantiomer. Preclinical species may be misleading as enzyme involvement and metabolism stereoselectivity can differ substantially.

Ontogeny of CYP3A and UGT activity in preterm piglets: a translational model for drug metabolism in preterm newborns

Despite considerable progress in understanding drug metabolism in the human pediatric population, data remains scarce in preterm neonates. Improving our knowledge of the ADME properties in this vulnerable age group is of utmost importance to avoid suboptimal dosing, which may lead to adverse drug reactions. The juvenile (mini)pig is a representative model for hepatic drug metabolism in human neonates and infants, especially phase I reactions. However, the effect of prematurity on the onset of hepatic phase I and phase II enzyme activity has yet to be investigated in this animal model. Therefore, the aim of this study was to assess the ontogeny of CYP3A and UGT enzyme activity in the liver of preterm (gestational day 105-107) and term-born (gestational day 115-117) domestic piglets. In addition, the ontogeny pattern between the preterm and term group was compared to examine whether postconceptional or postnatal age affects the onset of enzyme activity. The following age groups were included: preterm postnatal day (PND) 0 (n = 10), PND 5 (n = 10), PND 11 (n = 8), PND 26 (n = 10) and term PND 0 (n = 10), PND 5 (n = 10), PND 11 (n = 8), PND 19 (n = 18) and PND 26 (n = 10). Liver microsomes were extracted, and the metabolism of CYP3A and UGT-specific substrates assessed enzyme activity. Preterm CYP3A activity was only detectable at PND 26, whereas term CYP3A activity showed a gradual postnatal increase from PND 11 onwards. UGT activity gradually increased between PND 0 and PND 26 in preterm and term-born piglets, albeit, being systematically lower in the preterm group. Thus, postconceptional age is suggested as the main driver affecting porcine CYP3A and UGT enzyme ontogeny. These data are a valuable step forward in the characterization of the preterm piglet as a translational model for hepatic drug metabolism in the preterm human neonate.

Developmental pharmacokinetics of indomethacin in preterm neonates: Severely decreased drug clearance in the first week of life

Indomethacin is used commonly in preterm neonates for the prevention of intracranial hemorrhage and closure of an abnormally open cardiac vessel. Due to biomedical advances, the infants who receive this drug in the neonatal intensive care unit setting have become younger, smaller, and less mature (more preterm) at the time of treatment. To develop a pharmacokinetics (PK) model to aid future dosing, we designed a prospective cohort study to characterize indomethacin PK in a dynamically changing patient population. A population PK base model was created using NONMEM, and a covariate model was developed in a primary development cohort and subsequently was tested for accuracy in a validation cohort. Postnatal age was a significant covariate for hepatic clearance (CL_H ) and renal clearance (CL_R ). The typical value of the total clearance (CL, the sum of CL_R and CL_H ) was 3.09 ml/h and expressed as CL/WT_median  = 3.96 ml/h/kg, where WT_median is the median body weight. The intersubject variability of CL_R and CL_H were 61% and 207%, respectively. The typical value of the volume of distribution V_p  = 366 ml (V_p /WT_median  = 470 ml/kg), and its intersubject variability was 38.8%. Half-life was 82.1 h. Compared with more mature and older preterm populations studied previously, indomethacin CL is considerably lower in this contemporary population. Model-informed precision dosing incorporating important covariates other than weight alone offers an opportunity to individualize dosing in a susceptible patient undergoing rapid change.

Patent Ductus Arteriosus: A Contemporary Perspective for the Pediatric and Adult Cardiac Care Provider

The burden of patent ductus arteriosus (PDA) continues to be significant. In view of marked differences in preterm infants versus more mature, term counterparts (viewed on a continuum with adolescent and adult patients), mechanisms regulating ductal patency, genetic contributions, clinical consequences, and diagnostic and treatment thresholds are discussed separately, when appropriate. Among both preterm infants and older children and adults, a range of hemodynamic profiles highlighting the markedly variable consequences of the PDA are provided. In most contemporary settings, transcatheter closure is preferable over surgical ligation, but data on longer-term outcomes, particularly among preterm infants, are lacking. The present review provides recommendations to identify gaps in PDA diagnosis, management, and treatment on which subsequent research can be developed. Ultimately, the combination of refined diagnostic thresholds and expanded treatment options provides the best opportunities to address the burden of PDA. Although fundamental gaps remain unanswered, the present review provides pediatric and adult cardiac care providers with a contemporary framework in PDA care to support the practice of evidence-based medicine.

Analysis of influencing factors on monohydroxylated derivative of oxcarbazepine plasma concentration in children with epilepsy

PURPOSE

This study aimed to investigate the factors affecting the plasma concentration of monohydroxylated derivative (MHD) of oxcarbazepine (OXC) in children with epilepsy.

METHODS

We recruited 125 children with epilepsy who received OXC monotherapy. Among them, 16 single nucleotide polymorphisms were detected by MassARRAY genotyping technology to evaluate the influence of related factors on the plasma concentration of OXC monotherapy. MHD is the main active metabolite of OXC, and its plasma concentration was measured by high-performance liquid chromatography (HPLC).

RESULTS

Bivariate correlation analysis revealed that concentration-dose ratio (CDR) increased with weight, and the corresponding maintenance dose decreased with weight (r=0.317, P=0.001 for CDR; r=-0.285, P=0.000 for OXC maintenance dose). The duration of seizure was found to be associated with CDR (0.90 ± 0.36 vs 0.74 ± 0.26 μg·kg/mg/mL for ≥6 years vs <1 year, P=0.028; 0.90 ± 0.36 vs 0.64 ± 0.21 μg·kg/mg/mL for ≥6 years vs 1-3 years, P=0.004; 0.90 ± 0.36 vs 0.69 ± 0.18 μg·kg/mg/mL for ≥6 years vs 3-6 years, P=0.031). The CDR of patients with ABCB1 rs1045642 mutation homozygous GG type is higher than heterozygous AG type (0.79 ± 0.30 vs 0.68 ± 0.20 μg·kg/mg/mL for AG vs GG, P=0.032).

CONCLUSION

This study clarified the association of weight, duration of seizure, and gene polymorphisms of ABCB1 rs1045642 with MHD plasma concentration in children with epilepsy.

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

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

A potential implication of UDP-glucuronosyltransferase 2B10 in the detoxification of drugs used in pediatric hematopoietic stem cell transplantation setting: an in silico investigation

Background

Sinusoidal occlusion syndrome (SOS) is a potentially severe complication following hematopoietic stem cell transplantation (HSCT) in pediatric patients. Treatment related risk factors such as intensity of conditioning, hepatotoxic co-medication and patient related factors such as genetic variants predispose individuals to develop SOS. The variant allele for SNP rs17146905 in UDP-glucuronosyl transferase 2B10 (UGT2B10) gene was correlated with the occurrence of SOS in an exome-wide association study. UGT2B10 is a phase II drug metabolizing enzyme involved in the N-glucuronidation of tertiary amine containing drugs.

Methods

To shed light on the functionality of UGT2B10 enzyme in the metabolism of drugs used in pediatric HSCT setting, we performed in silico screening against custom based library of putative ligands. First, a list of potential substrates for in silico analysis was prepared using a systematic consensus-based strategy. The list comprised of drugs and their metabolites used in pediatric HSCT setting. The three-dimensional structure of UGT2B10 was not available from the Research Collaboratory Structural Bioinformatics - Protein Data Bank (RCSB - PDB) repository and thus we predicted the first human UGT2B10 3D model by using multiple template homology modeling with MODELLER Version 9.2 and molecular docking calculations with AutoDock Vina Version 1.2 were implemented to quantify the estimated binding affinity between selected putative substrates or ligands and UGT2B10. Finally, we performed molecular dynamics simulations using GROMACS Version 5.1.4 to confirm the potential UGT2B10 ligands prioritized after molecular docking (exhibiting negative free binding energy).

Results

Four potential ligands for UGT2B10 namely acetaminophen, lorazepam, mycophenolic acid and voriconazole n-oxide intermediate were identified. Other metabolites of voriconazole satisfied the criteria of being possible ligands of UGT2B10. Except for bilirubin and 4-Hydroxy Voriconazole, all the ligands (particularly voriconazole and hydroxy voriconazole) are oriented in substrate binding site close to the co-factor UDP (mean ± SD; 0.72 ± 0.33 nm). Further in vitro screening of the putative ligands prioritized by in silico pipeline is warranted to understand the nature of the ligands either as inhibitors or substrates of UGT2B10.

Conclusions

These results may indicate the clinical and pharmacological relevance UGT2B10 in pediatric HSCT setting. With this systematic computational methodology, we provide a rational-, time-, and cost-effective way to identify and prioritize the interesting putative substrates or inhibitors of UGT2B10 for further testing in in vitro experiments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-021-00402-5.

Characterizations of Human UDP-Glucuronosyltransferase Enzymes in the Conjugation of p-Cresol

p-Cresol is a uremic toxin that is formed by intestinal microbiota and extensively conjugated by first-pass metabolism. p-Cresol glucuronide exerts various forms of cellular toxicity in vitro and is accumulated in the plasma of subjects with kidney disease, where associations with adverse cardiovascular and renal outcomes are evident. The objective of this study was to determine the contributions of human UDP-glucuronosyltransferase (UGT) enzymes in the formation of p-cresol glucuronide. Utilizing commonly expressed hepatic or renal human recombinant UGTs (ie, hrUGT1A1, 1A3, 1A4, 1A6, 1A7, 1A8, 1A9, 1A10, 2B4, 2B7, 2B10, 2B15, and 2B17), hrUGT1A6 and hrUGT1A9 exhibited the highest catalytic activities in the generation of p-cresol glucuronide. The kinetics of p-cresol glucuronide formation in hrUGT1A6 and pooled human liver microsomes were best described by the Hill equation and in hrUGT1A9 and pooled human kidney microsomes by substrate inhibition. Using inhibitory and selective UGT inhibitors (ie, acetaminophen or amentoflavone for UGT1A6 and niflumic acid for UGT1A9), UGT1A6 was identified the predominant enzyme responsible for p-cresol glucuronide production in pooled human liver (78.4%-81.3% contribution) and kidney (54.3%-62.9%) microsomes, whereas UGT1A9 provided minor contributions (2.8% and 35.5%, respectively). The relative contributions of UGT1A6 (72.6 ± 11.3%, mean ± SD) and UGT1A9 (5.7 ± 4.1%) in individual human liver microsomes from 12 adult donors were highly variable, where an inverse association (R = -.784, p = .003) between UGT1A6 contribution and UGT1A9 probe substrate activity (ie, mycophenolic acid) was evident. Our novel findings provide valuable tools for conducting further mechanistic studies and for designing clinical interventions to mitigate the toxicities associated with p-cresol glucuronide.

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.

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.

Physiologically Based Pharmacokinetic Modeling Framework to Predict Neonatal Pharmacokinetics of Transplacentally Acquired Emtricitabine, Dolutegravir, and Raltegravir

BACKGROUND AND OBJECTIVE

Little is understood about neonatal pharmacokinetics immediately after delivery and during the first days of life following intrauterine exposure to maternal medications. Our objective was to develop and evaluate a novel, physiologically based pharmacokinetic modeling workflow for predicting perinatal and postnatal disposition of commonly used antiretroviral drugs administered prenatally to pregnant women living with human immunodeficiency virus.

METHODS

Using previously published, maternal-fetal, physiologically based pharmacokinetic models for emtricitabine, dolutegravir, and raltegravir built with PK-Sim/MoBi^®, placental drug transfer was predicted in late pregnancy. The total drug amount in fetal compartments at term delivery was estimated and subsequently integrated as initial conditions in different tissues of a whole-body, neonatal, physiologically based pharmacokinetic model to predict drug concentrations in the neonatal elimination phase after birth. Neonatal elimination processes were parameterized according to published data. Model performance was assessed by clinical data.

RESULTS

Neonatal physiologically based pharmacokinetic models generally captured the initial plasma concentrations after delivery but underestimated concentrations in the terminal phase. The mean percentage error for predicted plasma concentrations was - 71.5%, - 33.8%, and 76.7% for emtricitabine, dolutegravir, and raltegravir, respectively. A sensitivity analysis suggested that the activity of organic cation transporter 2 and uridine diphosphate glucuronosyltransferase 1A1 during the first postnatal days in term newborns is ~11% and ~30% of that in adults, respectively.

CONCLUSIONS

These findings demonstrate the general feasibility of applying physiologically based pharmacokinetic models to predict washout concentrations of transplacentally acquired drugs in newborns. These models can increase the understanding of pharmacokinetics during the first postnatal days and allow the prediction of drug exposure in this vulnerable population.

Enzyme Kinetics of Uridine Diphosphate Glucuronosyltransferases (UGTs)

Glucuronidation, catalyzed by uridine diphosphate glucuronosyltransferases (UGTs), is an important process for the metabolism and clearance of many lipophilic chemicals, including drugs, environmental chemicals, and endogenous compounds. Glucuronidation is a bisubstrate reaction that requires the aglycone and the cofactor, UDP-GlcUA. Accumulating evidence suggests that the bisubstrate reaction follows a compulsory-order ternary mechanism. To simplify the kinetic modeling of glucuronidation reactions in vitro, UDP-GlcUA is usually added to incubations in large excess. Many factors have been shown to influence UGT activity and kinetics in vitro, and these must be accounted for during experimental design and data interpretation. While the assessment of drug-drug interactions resulting from UGT inhibition has been challenging in the past, the increasing availability of UGT enzyme-selective substrate and inhibitor "probes" provides the prospect for more reliable reaction phenotyping and assessment of drug-drug interaction potential. Although extrapolation of the in vitro intrinsic clearance of a glucuronidated drug often underpredicts in vivo clearance, careful selection of in vitro experimental conditions and inclusion of extrahepatic glucuronidation may improve the predictivity of in vitro-in vivo extrapolation. Physiologically based pharmacokinetic (PBPK) modeling has also shown to be of value for predicting PK of drugs eliminated by glucuronidation.

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.

Pharmacokinetics and Pharmacodynamics of Posaconazole

Posaconazole is typically used for preventing invasive yeast and mold infections such as invasive aspergillosis in high-risk immunocompromised patients. The oral suspension was the first released formulation and many pharmacokinetic and pharmacodynamic studies of this formulation have been published. Erratic absorption profiles associated with this formulation were widely reported. Posaconazole exposure was found to be significantly influenced by food and many gastrointestinal conditions, including pH and motility. As a result, low posaconazole plasma concentrations were obtained in large groups of patients. These issues of erratic absorption urged the development of the subsequently marketed delayed-release tablet, which proved to be associated with higher and more stable exposure profiles. Shortly thereafter, an intravenous formulation was released for patients who are not able to take oral formulations. Both new formulations require a loading dose on day 1 to achieve high posaconazole concentrations more quickly, which was not possible with the oral suspension. So far, there appears to be no evidence of increased toxicity correlated to the higher posaconazole exposure achieved with the regimen for these formulations. The higher systemic availability of posaconazole for the delayed-release tablet and intravenous formulation have resulted in these two formulations being preferable for both prophylaxis and treatment of invasive fungal disease. This review aimed to integrate the current knowledge on posaconazole pharmacokinetics, pharmacodynamics, major toxicity, existing resistance, clinical experience in special populations, and new therapeutic strategies in order to get a clear understanding of the clinical use of this drug.

Infant exposure to dolutegravir through placental and breastmilk transfer: a population pharmacokinetic analysis of DolPHIN-1

BACKGROUND

Rapid reduction of HIV viral load is paramount to prevent peripartum transmission in women diagnosed late in pregnancy. We investigated dolutegravir population pharmacokinetics in maternal plasma, cord, breastmilk and infant plasma of DolPHIN-1 participants (NCT02245022) presenting with untreated HIV late in pregnancy (28-36 weeks gestation).

METHODS

Pregnant women from Uganda and South Africa were randomised (1:1) to daily dolutegravir (50 mg) or efavirenz-based therapy. Dolutegravir pharmacokinetic sampling (0-24 hours) was undertaken 14 days after treatment initiation and within 1-3 weeks of delivery, with matched maternal and cord samples at delivery. Mothers switched to efavirenz and maternal and infant plasma and breastmilk samples taken 24, 48 or 72 hours post-switch. Nonlinear mixed effects (NONMEM v. 7.4) was used to describe dolutegravir in all matrices and to evaluate covariates.

RESULTS

Twenty-eight women and 22 infants were included. Maternal dolutegravir was described by a two-compartment model linked to a fetal and breastmilk compartment. Cord and breastmilk to maternal plasma ratios were 1.279 (1.209-1.281) and 0.033 (0.021-0.050), respectively. Infant dolutegravir was described by breastmilk-to-infant and infant elimination rate constants. No covariate effects were observed. Predicted infant dolutegravir half-life and time to protein adjusted-IC90 (0.064 mg/L) for those above this threshold were 37.9 hours (22.1-63.5) and 108.9 hours [(18.6-129.6); 4.5 days (0.8-5.4); n=13].

CONCLUSIONS

Breastfeeding contributed relatively little to infant plasma exposures but a median of 4.5 days additional prophylaxis to some of the breastfed infants was observed following maternal dolutegravir cessation (3-15 days postpartum), which waned with time postpartum as transplacental dolutegravir cleared.

Use of computational toxicology (CompTox) tools to predict in vivo toxicity for risk assessment

Computational Toxicology tools were used to predict toxicity for three pesticides: propyzamide (PZ), carbaryl (CB) and chlorpyrifos (CPF). The tools used included: a) ToxCast/Tox21 assays (AC₅₀ s μM: concentration 50% maximum activity); b) in vitro-to-in vivo extrapolation (IVIVE) using ToxCast/Tox21 AC₅₀s to predict administered equivalent doses (AED: mg/kg/d) to compare to known in vivo Lowest-Observed-Effect-Level (LOEL)/Benchmark Dose (BMD); c) high throughput toxicokinetics population based (HTTK-Pop) using AC₅₀s for endpoints associated with the mode of action (MOA) to predict age-adjusted AED for comparison with in vivo LOEL/BMDs. ToxCast/Tox21 active-hit-calls for each chemical were predictive of targets associated with each MOA, however, assays directly relevant to the MOAs for each chemical were limited. IVIVE AEDs were predictive of in vivo LOEL/BMD₁₀s for all three pesticides. HTTK-Pop was predictive of in vivo LOEL/BMD₁₀s for PZ and CPF but not for CB after human age adjustments 11-15 (PZ) and 6-10 (CB) or 6-10 and 11-20 (CPF) corresponding to treated rat ages (in vivo endpoints). The predictions of computational tools are useful for risk assessment to identify targets in chemical MOAs and to support in vivo endpoints. Data can also aid is decisions about the need for further studies.

Genetic and Non-genetic Factors Contributing to the Significant Variation in the Plasma Trough Concentration-to-Dose Ratio of Valproic Acid in Children With Epilepsy

Objective: This study was conducted to evaluate the potential genetic and non-genetic factors contributing to plasma trough concentration-to-dose (C ₀/D) ratio of valproic acid (VPA) in pediatric patients with epilepsy. Study Design: A single-center, retrospective cohort study was performed by collecting data from 194 children aged 1-14 years between May 2018 and November 2018. The oral solution (n = 135) group and the sustained-release (SR) tablet group (n = 59) were defined, and the plasma VPA C ₀ was measured. Twenty-six single-nucleotide polymorphisms (SNPs) were chosen for genotyping with the MassARRAY system. A multiple logistic regression model was used for data analysis. Results: Body weight (BW) and age were positively correlated with the C ₀/D ratio in 194 patients, but the positive correlation disappeared after the patients were divided into oral solution and SR tablet subgroups. The average C ₀/D ratio was significantly increased by 2.11-fold (P = 0.000) in children who took VPA SR tablets compared with children who were administered VPA oral solutions. No significant association between genetic variants and the C ₀/D ratio was found, even for the five well-studied SNPs, namely UGT2B7 G211T, C802T, C161T, T125C, and CYP2C9 ^* 3 A1075C. However, a significant association between the C ₀/D ratio and UGT1A6/9 Del>A (rs144486213) was observed in the VPA oral solution group, but not in the VPA SR tablet group. Conclusions: The dosage forms of sodium valproate, rather than BW, age, or genetic polymorphisms, significantly affected the VPA C ₀/D ratios in pediatric patients with epilepsy. Based on our findings, switching the dosage form between solution and SR tablet should be performed cautiously. Total daily dose adjustment should be considered, and the plasma concentration, seizure-control effect, and adverse drug reaction should also be monitored very closely.