Quantitative profiling of human renal UDP-glucuronosyltransferases and glucuronidation activity: a comparison of normal and tumoral kidney tissues

[Precision Medicine for Patients with Renal Cell Carcinoma Based on Drug-metabolizing Enzyme Expression Levels]

Notable advances have recently been achieved in drug therapies for renal cell carcinoma (RCC). Several tyrosine kinase inhibitors (TKIs) and immune checkpoint inhibitors (ICIs) have been approved for metastatic RCC (mRCC). The current first-line treatment for mRCC involves combination therapies using TKIs and ICIs. However, there is no consensus on which TKI+ICI therapy is best or how to select the appropriate therapy for individual patients with RCC. The kidney expresses various metabolic enzymes, including CYP and uridine diphosphate glucose (UDP)-glucuronosyltransferase (UGT). Although information on CYP and UGT expression in the kidney is limited compared to our understanding of liver expression, the main CYP and UGT subtypes expressed at high levels in the kidney are estimated to be CYP2B6, CYP3A5, CYP4A11, CYP4F2, UGT1A6, UGT1A9, and UGT2B7. In RCC, the expression profiles and levels of these enzymes are somewhat altered compared with normal kidney. The main known subtypes of CYP and UGT in RCC are CYP1B1, CYP3A5, CYP4A11, UGT1A6, UGT1A9, UGT1A10, and UGT2B7. High CYP expression has been reported in several cancers, possibly conferring resistance to anti-cancer drugs including TKIs, due to extensive drug metabolism. Additionally, CYP and UGT expression levels may possibly affect cancer prognosis by metabolizing endogenous substrates, regardless of their role in anti-cancer drug metabolism. In this review, I discuss CYP and UGT expression level profiles in RCC based on previously published papers, including ours, and examine possible relationships between these enzyme expression profiles and treatment outcomes for patients with RCC.

Quantitative Proteomics for Translational Pharmacology and Precision Medicine: State of The Art and Future Outlook

Over the past 20 years, quantitative proteomics has contributed a wealth of protein expression data, which are currently used for a variety of systems pharmacology applications, as a complement or a surrogate for activity of the corresponding proteins. A symposium at the 25th North American International Society for the Study of Xenobiotics meeting, in Boston, in September 2023, was held to explore current and emerging applications of quantitative proteomics in translational pharmacology and strategies for improved integration into model-informed drug development based on practical experience of each of the presenters. A summary of the talks and discussions is presented in this perspective alongside future outlook that was outlined for future meetings. SIGNIFICANCE STATEMENT: This perspective explores current and emerging applications of quantitative proteomics in translational pharmacology and precision medicine and outlines the outlook for improved integration into model-informed drug development.

Clinical Pharmacokinetic and Pharmacodynamic Profile of Vericiguat

Vericiguat is an oral soluble guanylate cyclase stimulator and enhances the cyclic guanosine monophosphate pathway independently of nitric oxide as well as synergistically in normal- and low-nitric oxide conditions. This review describes the pharmacokinetic and pharmacodynamic profile of vericiguat and summarizes the effect of vericiguat on cardiac electrophysiology and population pharmacokinetic/pharmacodynamic relationships. Vericiguat demonstrates virtually complete absorption and increased exposure with food. Vericiguat has high oral bioavailability when taken with food (93.0%) with dose-proportional pharmacokinetics in healthy volunteers. Vericiguat has slightly less than dose-proportional pharmacokinetics with a slight decrease in bioavailability at higher doses in patients with heart failure (HF) with reduced ejection fraction (HFrEF). Vericiguat is a low-clearance drug, with a half-life of approximately 20 h in healthy volunteers and 30 h in patients with HFrEF. Most drug metabolism is achieved by glucuronidation. Vericiguat has pharmacodynamic effects as expected from its pharmacological mechanism of action (i.e., relaxation of the smooth muscles in the vasculature leading to changes in hemodynamics). In the VICTORIA trial (NCT02861534), which enrolled patients with HFrEF, no meaningful exposure-response relationships for the incidence of symptomatic hypotension or syncope were evident. There were no significant imbalances in the incidence of undesirable hemodynamic-related effects (symptomatic hypotension and syncope) in subgroups with HFrEF defined by sex, age, race, and renal impairment. In addition, most patients achieved the 10-mg target dose per the blood pressure-guided titration regimen. No dose adjustments due to body weight, age, sex, race, or hepatic/renal impairment are necessary in adult patients with HFrEF. Observed and predicted changes in vericiguat exposure when co-administered with perpetrator drugs were small and not clinically meaningful. In addition, vericiguat has low potential as a perpetrator to affect exposure and/or pharmacodynamic effects of drugs commonly prescribed in patients with heart failure; therefore, no dose adjustment of these drugs is required in patients taking vericiguat. There is limited experience on the combined use of vericiguat with long-acting nitrates in patients with HFrEF. The ongoing VICTOR trial (NCT05093933), which is investigating vericiguat in patients with HFrEF, permits the co-administration of long-acting nitrates. Combined use of vericiguat and phosphodiesterase type-5 inhibitors has not been studied in patients with HFrEF and is therefore not recommended because of the potential increased risk for symptomatic hypotension. Vericiguat was not associated with electrophysiological abnormalities in preclinical and clinical studies up to the approved dose of 10 mg at steady state. Vericiguat is approved for the treatment of recently decompensated patients with worsening HFrEF. Vericiguat's safety and efficacy profile in patients with HFrEF will be further characterized by the VICTOR trial (NCT05093933) in adults without recent decompensation and in a pediatric population with HF due to left ventricular systolic dysfunction (VALOR trial, NCT05714085).

Identification of UDP-glucuronosyltransferase (UGT) isoforms involved in the metabolism of Chlorophenols (CPs)

Chlorophenols (CPs) are a group of pollutants that pose a great threat to the environment, they are widely used in industrial and agricultural wastes, pesticides, herbicides, textiles, pharmaceuticals and plastics. Among CPs, pentachlorophenol was listed as one of the persistent organic pollutants (POPs) by the Stockholm convention. This study aims to identify the UDP-glucosyltransferase (UGT) isoforms involved in the metabolic elimination of CPs. CPs' mono-glucuronide was detected in the human liver microsomes (HLMs) incubation mixture with co-factor uridine-diphosphate glucuronic acid (UDPGA). HLMs-catalyzed glucuronidation metabolism reaction equations followed Michaelis-Menten or substrate inhibition type. Recombinant enzymes and chemical reagents inhibition experiments were utilized to phenotype the main UGT isoforms involved in the glucuronidation of CPs. UGT1A6 might be the major enzyme in the glucuronidation of mono-chlorophenol isomer. UGT1A1, UGT1A6, UGT1A9, UGT2B4 and UGT2B7 were the most important five UGT isoforms for metabolizing the di-chlorophenol and tri-chlorophenol isomers. UGT1A1 and UGT1A3 were the most important UGT isoforms in the catalysis of tetra-chlorophenol and pentachlorophenol isomers. Species differences were investigated using rat liver microsomes (RLMs), pig liver microsomes (PLMs), dog liver microsomes (DLMs), and monkey liver microsomes (MyLMs). All these results were helpful for elucidating the metabolic elimination and toxicity of CPs.

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.

Quantitative analysis of the UDP-glucuronosyltransferase transcriptome in human tissues

UDP-glucuronosyltransferases (UGTs) are phase II drug metabolizing enzymes that play important roles in the detoxification of endogenous and exogenous substrates. The 22 human UGTs belong to four families (UGT1, UGT2, UGT3, and UGT8) and differ in their expression, substrate specificity, UDP-sugar preference, and physiological functions. Differential expression/activity of the UGTs contributes to interperson variability in drug responses and toxicity, hormone homeostasis, and disease/cancer risks. However, in normal tissues, the tissue-specific expression profiles and transcriptional regulation of the UGTs are still not fully understood. In this study, we comprehensively analyzed the transcriptome of 22 UGTs in 54 human tissues/regions using RNAseq data from GTEx. We then validated the findings in the liver and small intestine samples using real-time PCR. Our results showed large interindividual variability across tissues in the expression of each UGT and the overall composition of UGT pools, consisting of different UGTs and their splice isoforms. Our results also revealed coexpression of the UGTs, Cytochrome P450s, and many transcription factors in the liver, suggesting potential coregulation or functional coordination. Our results provide the groundwork for future studies to detail further the regulation of the expression and activity of the UGTs.

Understanding Statin-Roxadustat Drug-Drug-Disease Interaction Using Physiologically-Based Pharmacokinetic Modeling

A different drug-drug interaction (DDI) scenario may exist in patients with chronic kidney disease (CKD) compared with healthy volunteers (HVs), depending on the interplay between drug-drug and disease (drug-drug-disease interaction (DDDI)). Physiologically-based pharmacokinetic (PBPK) modeling, in lieu of a clinical trial, is a promising tool for evaluating these complex DDDIs in patients. However, the prediction confidence of PBPK modeling in the severe CKD population is still low when nonrenal pathways are involved. More mechanistic virtual disease population and robust validation cases are needed. To this end, we aimed to: (i) understand the implications of severe CKD on statins (atorvastatin, simvastatin, and rosuvastatin) pharmacokinetics (PK) and DDI; and (ii) predict untested clinical scenarios of statin-roxadustat DDI risks in patients to guide suitable dose regimens. A novel virtual severe CKD population was developed incorporating the disease effect on both renal and nonrenal pathways. Drug and disease PBPK models underwent a four-way validation. The verified PBPK models successfully predicted the altered PKs in patients for substrates and inhibitors and recovered the observed statin-rifampicin DDIs in patients and the statin-roxadustat DDIs in HVs within 1.25- and 2-fold error. Further sensitivity analysis revealed that the severe CKD effect on statins PK is mainly mediated by hepatic BCRP for rosuvastatin and OATP1B1/3 for atorvastatin. The magnitude of statin-roxadustat DDI in patients with severe CKD was predicted to be similar to that in HVs. PBPK-guided suitable dose regimens were identified to minimize the risk of side effects or therapeutic failure of statins when co-administered with roxadustat.

Extracellular Vesicles as Surrogates for Drug Metabolism and Clearance: Promise vs. Reality

Drug-metabolizing enzymes (DMEs) and transporters play a major role in drug efficacy and safety. They are regulated at multiple levels and by multiple factors. Estimating their expression and activity could contribute to predicting drug pharmacokinetics and their regulation by drugs or pathophysiological situations. Determining the expression of these proteins in the liver, intestine, and kidney requires the collection of biopsy specimens. Instead, the isolation of extracellular vesicles (EVs), which are nanovesicles released by most cells and present in biological fluids, could deliver this information in a less invasive way. In this article, we review the use of EVs as surrogates for the expression and activity of DMEs, uptake, and efflux transporters. Preliminary evidence has been provided for a correlation between the expression of some enzymes and transporters in EVs and the tissue of origin. In some cases, data obtained in EVs reflect the induction of phase I-DMEs in the tissues. Further studies are required to elucidate to what extent the regulation of other DMEs and transporters in the tissues reflects in the EV cargo. If an association between tissues and their EVs is firmly established, EVs may represent a significant advancement toward precision therapy based on the biotransformation and excretion capacity of each individual.

Quantification of Accurate Composition and Total Abundance of Homologous Proteins by Conserved-Plus-Surrogate Peptide Approach: Quantification of UDP Glucuronosyltransferases in Human Tissues

Characterization of accurate compositions and total abundance of homologous drug-metabolizing enzymes, such as UDP glucuronosyltransferases (UGTs), is important for predicting the fractional contribution of individual isoforms involved in the metabolism of a drug for applications in physiologically based pharmacokinetic (PBPK) modeling. Conventional targeted proteomics utilizes surrogate peptides, which often results in high technical and interlaboratory variability due to peptide-specific digestion leading to data inconsistencies. To address this problem, we developed a novel conserved-plus-surrogate peptide (CPSP) approach for determining the accurate compositions and total or cumulative abundance of homologous UGTs in commercially available pooled human liver microsomes (HLM), human intestinal microsomes (HIM), human kidney microsomes (HKM), and human liver S9 (HLS9) fraction. The relative percent composition of UGT1A and UGT2B isoforms in the human liver was 35:5:36:11:13 for UGT1A1:1A3:1A4:1A6:1A9 and 20:32:22:21:5 for UGT2B4:2B7:2B10:2B15:2B17. The human kidney and intestine also showed unique compositions of UGT1As and UGT2Bs. The reproducibility of the approach was validated by assessing correlations of UGT compositions between HLM and HLS9 (R2> 0.91). The analysis of the conserved peptides also provided the abundance for individual UGT isoforms included in this investigation as well as the total abundance (pmol/mg protein) of UGT1As and UGT2Bs across tissues, i.e., 268 and 342 (HLM), 21 and 92 (HIM), and 138 and 99 (HKM), respectively. The CPSP approach could be used for applications in the in-vitro-to-in-vivo extrapolation of drug metabolism and PBPK modeling. SIGNIFICANCE STATEMENT: We quantified the absolute compositions and total abundance of UDP glucuronosyltransferases (UGTs) in pooled human liver, intestine, and kidney microsomes using a novel conserved-plus-surrogate peptide (CPSP) approach. The CPSP approach addresses the surrogate peptide-specific variability in the determination of the absolute composition of UGTs. The data presented in this manuscript are applicable for the estimation of the fraction metabolized by individual UGTs towards better in vitro-to-in vivo extrapolation of UGT-mediated drug metabolism.

Significance of UGT1A6, UGT1A9, and UGT2B7 genetic variants and their mRNA expression in the clinical outcome of renal cell carcinoma

UDP-glucuronosyltransferase (UGT) metabolizes a number of endogenous and exogenous substrates. Renal cells express high amounts of UGT; however, the significance of UGT in patients with renal cell carcinoma (RCC) remains unknown. In this study, we profile the mRNA expression of UGT subtypes (UGT1A6, UGT1A9, and UGT2B7) and their genetic variants in the kidney tissue of 125 Japanese patients with RCC (Okayama University Hospital, Japan). In addition, we elucidate the association between the UGT variants and UGT mRNA expression levels and clinical outcomes in these patients. The three representative genetic variants, namely, UGT1A6 541A > G, UGT1A9 i399C > T, and UGT2B7-161C > T, were genotyped, and their mRNA expression levels in each tissue were determined. We found that the mRNA expression of the three UGTs (UGT1A6, UGT1A9, and UGT2B7) are significantly downregulated in RCC tissues. Moreover, in patients with RCC, the UGT2B7-161C > T variant and high UGT2B7 mRNA expression are significantly correlated with preferable cancer-specific survival (CSS) and overall survival (OS), respectively. As such, the UGT2B7-161C > T variant and UGT2B7 mRNA expression level were identified as significant independent prognostic factors of CSS and CSS/OS, respectively. Taken together, these findings indicate that UGT2B7 has a role in RCC progression and may, therefore, represent a potential prognostic biomarker for patients with RCC.

The Impact of Suboptimal 25-Hydroxyvitamin D Levels and Cholecalciferol Replacement on the Pharmacokinetics of Oral Midazolam in Control Subjects and Patients With Chronic Kidney Disease

The aim of this study was to investigate the impact of suboptimal 25-hydroxyvitamin D (25-VitD) and cholecalciferol (VitD3 ) supplementation on the pharmacokinetics of oral midazolam (MDZ) in control subjects and subjects with chronic kidney disease (CKD). Subjects with CKD (n = 14) and controls (n = 5) with suboptimal 25-VitD levels (<30 ng/mL) were enrolled in a 2-phase study. In phase 1 (suboptimal), subjects were administered a single oral dose of VitD3 (5000 IU) and MDZ (2 mg). In phase 2 (replete) subjects who achieved 25-VitD repletion after receiving up to 16 weeks of daily cholecalciferol were given the identical single oral doses of VitD3 and MDZ as in phase 1. Concentrations of MDZ and metabolites, 1'-hydroxymidazolam (1'-OHMDZ), and 1'-OHMDZ glucuronide (1'-OHMDZ-G) were measured by liquid chromatography-tandem mass spectrometry and pharmacokinetic analysis was performed. Under suboptimal 25-VitD, reductions in MDZ clearance and renal clearance of 47% and 87%, respectively, and a 72% reduction in renal clearance of 1'-OHMDZ-G were observed in CKD vs controls. In phase 1 versus phase 2, MDZ clearance increased in all control subjects, with a median (interquartile range) increase of 10.5 (0.62-16.7) L/h. No changes in MDZ pharmacokinetics were observed in subjects with CKD between phases 1 and 2. The effects of 25-VitD repletion on MDZ disposition was largely observed in subjects without kidney disease. Impaired MDZ metabolism and/or excretion alterations due to CKD in a suboptimal 25-VitD state may not be reversed by cholecalciferol therapy. Suboptimal 25-VitD may augment the reductions in MDZ and 1'-OHMDZ-G clearance values observed in patients with CKD.

The Somatic Mutation Landscape of UDP-Glycosyltransferase (UGT) Genes in Human Cancers

The human UDP-glycosyltransferase (UGTs) superfamily has a critical role in the metabolism of anticancer drugs and numerous pro/anti-cancer molecules (e.g., steroids, lipids, fatty acids, bile acids and carcinogens). Recent studies have shown wide and abundant expression of UGT genes in human cancers. However, the extent to which UGT genes acquire somatic mutations within tumors remains to be systematically investigated. In the present study, our comprehensive analysis of the somatic mutation profiles of 10,069 tumors from 33 different TCGA cancer types identified 3427 somatic mutations in UGT genes. Overall, nearly 18% (1802/10,069) of the assessed tumors had mutations in UGT genes with huge variations in mutation frequency across different cancer types, ranging from over 25% in five cancers (COAD, LUAD, LUSC, SKCM and UCSC) to less than 5% in eight cancers (LAML, MESO, PCPG, PAAD, PRAD, TGCT, THYM and UVM). All 22 UGT genes showed somatic mutations in tumors, with UGT2B4, UGT3A1 and UGT3A2 showing the largest number of mutations (289, 307 and 255 mutations, respectively). Nearly 65% (2260/3427) of the mutations were missense, frame-shift and nonsense mutations that have been predicted to code for variant UGT proteins. Furthermore, about 10% (362/3427) of the mutations occurred in non-coding regions (5' UTR, 3' UTR and splice sites) that may be able to alter the efficiency of translation initiation, miRNA regulation or the splicing of UGT transcripts. In conclusion, our data show widespread somatic mutations of UGT genes in human cancers that may affect the capacity of cancer cells to metabolize anticancer drugs and endobiotics that control pro/anti-cancer signaling pathways. This highlights their potential utility as biomarkers for predicting therapeutic efficacy and clinical outcomes.

Physiologically based pharmacokinetic combined JAK2 occupancy modelling to simulate PK and PD of baricitinib with kidney transporter inhibitors and in patients with hepatic/renal impairment

PURPOSE

Our aim is to build a physiologically based pharmacokinetic and JAK2 occupancy model (PBPK-JO) to simultaneously predict pharmacokinetic (PK) and pharmacodynamic (PD) changes of baricitinib (BAR) in healthy humans when co-administrated with kidney transporters OAT3 and MATE2-K inhibitors, and in patients with hepatic and renal impairment.

METHODS

Probenecid and vandetanib were selected as OAT3 and MATE2-K competitive inhibitors, respectively. The PBPK-JO model was built using physicochemical and biochemical properties of BAR, and then verified by observed clinical PK. Finally, the model was applied to determine optimal dosing regimens in various clinical situations.

RESULTS

Here, we have successfully simulated PK and JAK2 occupancy profiles in humans by PBPK-JO model. Moreover, this modelling reproduced every observed PK data, and every mean relative deviation (MRD) was below 2. The simulation suggested that PK of BAR had a significant change (2.22-fold increase), however PD only had a slight increase of 1.14-fold. Additionally, the simulation also suggested that vandetanib was almost unlikely to affect the PK and PD of BAR. In simulations of hepatic and renal impairment patients, the predictions suggested that significant changes in the PK and PD of BAR occurred. However, there was a lower fold increase in JAK2 occupancy than in PK in patients relative to healthy individuals.

CONCLUSION

Administration dose adjustment of BAR when co-administrated with OAT3 inhibitors or in patients with hepatic or renal impairment should combine PK and PD changes of BAR, instead of only considering PK alteration.

Physiologically Based Pharmacokinetic (PBPK) Modeling of Clopidogrel and Its Four Relevant Metabolites for CYP2B6, CYP2C8, CYP2C19, and CYP3A4 Drug–Drug–Gene Interaction Predictions

The antiplatelet agent clopidogrel is listed by the FDA as a strong clinical index inhibitor of cytochrome P450 (CYP) 2C8 and weak clinical inhibitor of CYP2B6. Moreover, clopidogrel is a substrate of—among others—CYP2C19 and CYP3A4. This work presents the development of a whole-body physiologically based pharmacokinetic (PBPK) model of clopidogrel including the relevant metabolites, clopidogrel carboxylic acid, clopidogrel acyl glucuronide, 2-oxo-clopidogrel, and the active thiol metabolite, with subsequent application for drug–gene interaction (DGI) and drug–drug interaction (DDI) predictions. Model building was performed in PK-Sim^® using 66 plasma concentration-time profiles of clopidogrel and its metabolites. The comprehensive parent-metabolite model covers biotransformation via carboxylesterase (CES) 1, CES2, CYP2C19, CYP3A4, and uridine 5′-diphospho-glucuronosyltransferase 2B7. Moreover, CYP2C19 was incorporated for normal, intermediate, and poor metabolizer phenotypes. Good predictive performance of the model was demonstrated for the DGI involving CYP2C19, with 17/19 predicted DGI AUC_last and 19/19 predicted DGI C_max ratios within 2-fold of their observed values. Furthermore, DDIs involving bupropion, omeprazole, montelukast, pioglitazone, repaglinide, and rifampicin showed 13/13 predicted DDI AUC_last and 13/13 predicted DDI C_max ratios within 2-fold of their observed ratios. After publication, the model will be made publicly accessible in the Open Systems Pharmacology repository.

Inhibition of UDP-Glucuronosyltransferase Enzymes by Major Cannabinoids and Their Metabolites

The UDP-glucuronosyltransferase (UGT) family of enzymes play a central role in the metabolism and detoxification of a wide range of endogenous and exogenous compounds. UGTs exhibit a high degree of structural similarity and display overlapping substrate specificity, often making estimations of potential drug-drug interactions difficult to fully elucidate. One such interaction yet to be examined may be occurring between UGTs and cannabinoids, as the legalization of recreational and medicinal cannabis and subsequent co-usage of cannabis and therapeutic drugs increases in the United States and internationally. In the present study, the inhibition potential of the major cannabinoids Δ9-tetrahydrocannabinol (THC), cannabidiol (CBD), and cannabinol (CBN), as well as their major metabolites, was determined in microsomes isolated from HEK293 cells overexpressing individual recombinant UGTs and in microsomes from human liver and kidney specimens. The highest inhibition was seen by CBD against the glucuronidation activity of UGTs 1A9, 2B4, 1A6, and 2B7, with binding-corrected IC50 values of 0.12 ± 0.020 µM, 0.22 ± 0.045 µM, 0.40 ± 0.10 µM, and 0.82 ± 0.15 µM, respectively. Strong inhibition of UGT1A9 was also demonstrated by THC and CBN, with binding-corrected IC50 values of 0.45 ± 0.12 μM and 0.51 ± 0.063 μM, respectively. Strong inhibition of UGT2B7 was also observed for THC and CBN; no or weak inhibition was observed with cannabinoid metabolites. This inhibition of UGT activity suggests that in addition to playing an important role in drug-drug interactions, cannabinoid exposure may have important implications in patients with impaired hepatic or kidney function. SIGNIFICANCE STATEMENT: Major cannabinoids found in the plasma of cannabis users inhibit several UDP-glucuronosyltransferase (UGT) enzymes, including UGT1A6, UGT1A9, UGT2B4, and UGT2B7. This study is the first to show the potential of cannabinoids and their metabolites to inhibit all the major kidney UGTs as well as the two most abundant UGTs present in liver. This study suggests that as all three major kidney UGTs are inhibited by cannabinoids, greater drug-drug interaction effects might be observed from co-use of cannabinods and therapeutics that are cleared renally.

The Expression Profiles and Deregulation of UDP-Glycosyltransferase (UGT) Genes in Human Cancers and Their Association with Clinical Outcomes

The human UDP-glycosyltransferase (UGTs) superfamily has 22 functional enzymes that play a critical role in the metabolism of small lipophilic compounds, including carcinogens, drugs, steroids, lipids, fatty acids, and bile acids. The expression profiles of UGT genes in human cancers and their impact on cancer patient survival remains to be systematically investigated. In the present study, a comprehensive analysis of the RNAseq and clinical datasets of 9514 patients from 33 different TCGA (the Genome Cancer Atlas) cancers demonstrated cancer-specific UGT expression profiles with high interindividual variability among and within individual cancers. Notably, cancers derived from drug metabolizing tissues (liver, kidney, gut, pancreas) expressed the largest number of UGT genes (COAD, KIRC, KIRP, LIHC, PAAD); six UGT genes (1A6, 1A9, 1A10, 2A3, 2B7, UGT8) showed high expression in five or more different cancers. Kaplan-Meier plots and logrank tests revealed that six UGT genes were significantly associated with increased overall survival (OS) rates [UGT1A1 (LUSC), UGT1A6 (ACC), UGT1A7 (ACC), UGT2A3 (KIRC), UGT2B15 (BLCA, SKCM)] or decreased OS rates [UGT2B15 (LGG), UGT8 (UVM)] in specific cancers. Finally, differential expression analysis of 611 patients from 12 TCGA cancers identified 16 UGT genes (1A1, 1A3, 1A6, 1A7, 1A8, 1A9, 1A10, 2A1, 2A3, 2B4, 2B7, 2B11, 2B15, 3A1, 3A2, UGT8) that were up/downregulated in at least one cancer relative to normal tissues. In conclusion, our data show widespread expression of UGT genes in cancers, highlighting the capacity for intratumoural drug metabolism through the UGT conjugation pathway. The data also suggests the potentials for specific UGT genes to serve as prognostic biomarkers or therapeutic targets in cancers.

PBPK Modeling as a Tool for Predicting and Understanding Intestinal Metabolism of Uridine 5'-Diphospho-glucuronosyltransferase Substrates

Uridine 5′-diphospho-glucuronosyltransferases (UGTs) are expressed in the small intestines, but prediction of first-pass extraction from the related metabolism is not well studied. This work assesses physiologically based pharmacokinetic (PBPK) modeling as a tool for predicting intestinal metabolism due to UGTs in the human gastrointestinal tract. Available data for intestinal UGT expression levels and in vitro approaches that can be used to predict intestinal metabolism of UGT substrates are reviewed. Human PBPK models for UGT substrates with varying extents of UGT-mediated intestinal metabolism (lorazepam, oxazepam, naloxone, zidovudine, cabotegravir, raltegravir, and dolutegravir) have demonstrated utility for predicting the extent of intestinal metabolism. Drug–drug interactions (DDIs) of UGT1A1 substrates dolutegravir and raltegravir with UGT1A1 inhibitor atazanavir have been simulated, and the role of intestinal metabolism in these clinical DDIs examined. Utility of an in silico tool for predicting substrate specificity for UGTs is discussed. Improved in vitro tools to study metabolism for UGT compounds, such as coculture models for low clearance compounds and better understanding of optimal conditions for in vitro studies, may provide an opportunity for improved in vitro–in vivo extrapolation (IVIVE) and prospective predictions. PBPK modeling shows promise as a useful tool for predicting intestinal metabolism for UGT substrates.

Quantitative Proteomic Map of Enzymes and Transporters in the Human Kidney: Stepping Closer to Mechanistic Kidney Models to Define Local Kinetics

The applications of translational modeling of local drug concentrations in various organs had a sharp increase over the last decade. These are part of the model-informed drug development initiative, adopted by the pharmaceutical industry and promoted by drug regulatory agencies. With respect to the kidney, the models serve as a bridge for understanding animal vs. human observations related to renal drug disposition and any consequential adverse effects. However, quantitative data on key drug-metabolizing enzymes and transporters relevant for predicting renal drug disposition are limited. Using targeted and global quantitative proteomics, we determined the abundance of multiple enzymes and transporters in 20 human kidney cortex samples. Nine enzymes and 22 transporters were quantified (8 for the first time in the kidneys). In addition, > 4,000 proteins were identified and used to form an open database. CYP2B6, CYP3A5, and CYP4F2 showed comparable, but generally low expression, whereas UGT1A9 and UGT2B7 levels were the highest. Significant correlation between abundance and activity (measured by mycophenolic acid clearance) was observed for UGT1A9 (Rs = 0.65, P = 0.004) and UGT2B7 (Rs = 0.70, P = 0.023). Expression of P-gp ≈ MATE-1 and OATP4C1 transporters were high. Strong intercorrelations were observed between several transporters (P-gp/MRP4, MRP2/OAT3, and OAT3/OAT4); no correlation in expression was apparent for functionally related transporters (OCT2/MATEs). This study extends our knowledge of pharmacologically relevant proteins in the kidney cortex, with implications on more prudent use of mechanistic kidney models under the general framework of quantitative systems pharmacology and toxicology.

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.

Population Pharmacokinetics and Pharmacodynamics of Vericiguat in Patients with Heart Failure and Reduced Ejection Fraction

Background

Vericiguat, a stimulator of soluble guanylate cyclase, has been developed as a first-in-class therapy for worsening chronic heart failure in adults with left ventricular ejection fraction < 45%.

Objective

The objective of this article was to characterize the pharmacokinetics and pharmacokinetic variability of vericiguat combined with guideline-directed medical therapy (standard of care), and identify exposure–response relationships for safety (hemodynamics) and pharmacodynamic markers of efficacy (N-terminal pro-B-type natriuretic peptide concentration [NT-proBNP]) in patients with heart failure and left ventricular ejection fraction < 45% in the SOCRATES-REDUCED study (NCT01951625).

Methods

Vericiguat and NT-proBNP plasma concentrations in 454 and 432 patients in SOCRATES-REDUCED, respectively, were analyzed using nonlinear mixed-effects modeling.

Results

Vericiguat pharmacokinetics were well described by a one-compartment model with apparent clearance, apparent volume of distribution, and absorption rate constant. Age, bodyweight, plasma bilirubin, and creatinine clearance were identified as significant covariates on apparent clearance; sex and bodyweight on apparent volume of distribution; and bodyweight and plasma albumin level on absorption rate constant. Pharmacokinetic/pharmacodynamic analysis showed initial minor and transient effects of vericiguat on blood pressure with low clinical impact. There were no changes in heart rate following initial or repeated vericiguat administration. An exposure-dependent and time-dependent turnover pharmacokinetic/pharmacodynamic model for NT-proBNP described production and elimination rates and an demonstrated exposure-dependent reduction in [NT-proBNP] by vericiguat plus standard of care compared with placebo plus standard of care. This effect was dependent on baseline [NT-proBNP].

Conclusions

Vericiguat has predictable pharmacokinetics, with no long-term effects on blood pressure in patients with heart failure and left ventricular ejection fraction < 45%. A pharmacokinetic/pharmacodynamic model described a vericiguat exposure-dependent reduction of NT-proBNP.

Clinical Trial Identifier

NCT01951625.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40262-021-01024-y.

Transit and Metabolic Pathways of Quercetin in Tubular Cells: Involvement of Its Antioxidant Properties in the Kidney

Quercetin is a flavonoid with antioxidant, antiviral, antimicrobial, and anti-inflammatory properties. Therefore, it has been postulated as a molecule with great therapeutic potential. The renoprotective capacity of quercetin against various toxins that produce oxidative stress, in both in vivo and in vitro models, has been shown. However, it is not clear whether quercetin itself or any of its metabolites are responsible for the protective effects on the kidney. Although the pharmacokinetics of quercetin have been widely studied and the complexity of its transit throughout the body is well known, the metabolic processes that occur in the kidney are less known. Because of that, the objective of this review was to delve into the molecular and cellular events triggered by quercetin and/or its metabolites in the tubular cells, which could explain some of the protective properties of this flavonoid against oxidative stress produced by toxin administration. Thus, the following are analyzed: (1) the transit of quercetin to the kidney; (2) the uptake mechanisms of quercetin and its metabolites from plasma to the tubular cells; (3) the metabolic processes triggered in those cells, which affect the accumulation of metabolites in the intracellular space; and (4) the efflux mechanisms of these compounds and their subsequent elimination through urine. Finally, it is discussed whether those processes that are mediated in the tubular cells and that give rise to different metabolites are related to the antioxidant and renoprotective properties observed after the administration of quercetin.

Clinical Investigation of Metabolic and Renal Clearance Pathways Contributing to the Elimination of Fevipiprant Using Probenecid as Perpetrator

Fevipiprant, an oral, nonsteroidal, highly selective, reversible, and competitive prostaglandin D2 receptor 2 antagonist, is eliminated by glucuronidation and by direct renal excretion predominantly via organic anion transporter (OAT) 3. This study aimed to assess the effect of simultaneous UDP-glucuronosyltransferase (UGT) and OAT3 inhibition by probenecid on the pharmacokinetics of fevipiprant and its acyl glucuronide (AG) metabolite to support the dosing recommendation of fevipiprant in the presence of drugs inhibiting these pathways; however, phase III clinical trial results did not support its submission. This was a single-center, open-label, single-sequence, two-period crossover study in healthy subjects. Liquid chromatography with tandem mass spectrometry was used to measure concentrations of fevipiprant and its AG metabolite in plasma and urine. In the presence of probenecid, the mean maximum concentrations of fevipiprant increased approximately 1.7-fold, and the area under the concentration-time curve in plasma increased approximately 2.5-fold, whereas the mean apparent volume of distribution and the AG metabolite:fevipiprant ratio decreased. The apparent systemic clearance decreased by approximately 60% and the renal clearance decreased by approximately 88% in the presence of probenecid. Using these data and those from previous studies, the relative contribution of OAT and UGT inhibition to the overall effect of probenecid was estimated. Furthermore, a general disposition scheme for fevipiprant was developed, showing how a perpetrator drug such as probenecid, which interferes with two key elimination pathways of fevipiprant, causes only a moderate increase in exposure and allows estimation of the drug-drug inhibition when only one of the two pathways is inhibited. SIGNIFICANCE STATEMENT: In this drug-drug interaction (DDI) study, probenecid was used as a tool to inhibit both glucuronidation and active renal secretion of fevipiprant. The combination of plasma and urine pharmacokinetic data from this study with available data allowed the development of a quantitative scheme to describe the fate of fevipiprant in the body, illustrating why the DDI effect on fevipiprant is weak-to-moderate even if a perpetrator drug inhibits several elimination pathways.

A fluorescence-based microplate assay for high-throughput screening and evaluation of human UGT inhibitors

Human UDP-glucuronosyltransferase enzymes (hUGTs), one of the most important classes of conjugative enzymes, are responsible for the glucuronidation and detoxification of a variety of endogenous substances and xenobiotics. Inhibition of hUGTs may cause undesirable effects or adverse drug-drug interactions (DDI) via modulating the glucuronidation rates of endogenous toxins or the drugs that are primarily conjugated by the inhibited hUGTs. Herein, to screen hUGTs inhibitors in a more efficient way, a novel fluorescence-based microplate assay has been developed by utilizing a fluorogenic substrate. Following screening of series of 4-hydroxy-1,8-naphthalimide derivatives, we found that 4-HN-335 is a particularly good substrate for a panel of hUGTs. Under physiological conditions, 4-HN-335 can be readily O-glucuronidated by ten hUGTs, such reactions generate a single O-glucuronide with a high quantum yield (Ф = 0.79) and bring remarkable changes in fluorescence emission. Subsequently, a fluorescence-based microplate assay is developed to simultaneously measure the inhibitory effects of selected compound(s) on ten hUGTs. The newly developed fluorescence-based microplate assay is time- and cost-saving, easy to manage and can be adapted for 96-well microplate format with the Z-factor of 0.92. We further demonstrate the utility of the fluorescence-based assay for high-throughput screening of two compound libraries, resulting in the identification of several potent UGT inhibitors, including natural products and FDA-approved drugs. Collectively, this study reports a novel fluorescence-based microplate assay for simultaneously sensing the residual activities of ten hUGTs, which strongly facilitates the identification and characterization of UGT inhibitors from drugs or herbal constituents and the investigations on UGT-mediated DDI.

Pharmacokinetics of the oral multikinase inhibitor regorafenib and its association with real-world treatment outcomes

Purpose Despite the established activity of regorafenib in metastatic colorectal cancer (CRC), gastrointestinal stromal tumor (GIST), and hepatocellular carcinoma (HCC), its toxicity profile has limited clinical use. We aimed to evaluate the pharmacokinetics of regorafenib and its active metabolites M-2/M-5, and to clarify the relationships between total drug-related exposure and clinical outcomes in real-world practice. Methods Blood samples at steady state were obtained during Cycle 1 from patients treated with regorafenib. Plasma concentrations of regorafenib and its metabolites were measured by liquid chromatography-tandem mass spectrometry. The efficacy and safety endpoints were progression-free survival (PFS) and dose-limiting toxicities (DLTs), respectively. The exposure-response relationships were assessed. Results Thirty-four Japanese patients with advanced cancers were enrolled (CRC, n = 26; GIST and HCC, each n = 4). Nine patients started regorafenib treatment at the recommended dose of 160 mg once daily (3 weeks on / 1 week off), while the other patients received a reduced starting dose to minimize toxicities. The median PFS was significantly longer in patients achieving total trough concentrations (C_trough) of regorafenib and M-2/M-5 ≥2.9 µg/mL than those who did not (112 vs. 57 days; p = 0.044). Furthermore, the cumulative incidence of DLTs during the first 2 cycles was significantly higher in patients with summed C_trough levels ≥4.3 µg/mL than in others (p = 0.0003). Conclusions Dose titration of regorafenib to achieve drug-related C_trough levels between 2.9 and 4.3 µg/mL in Cycle 1 may improve efficacy and safety, warranting further investigation in a larger patient population.Clinical trial registry: Not applicable.

Metabolism of Long-Acting Rilpivirine After Intramuscular Injection: HIV Prevention Trials Network Study 076 (HPTN 076)

A long-acting injectable formulation of rilpivirine (RPV), a non-nucleoside reverse transcriptase inhibitor, is currently under investigation for use in human immunodeficiency virus (HIV) maintenance therapy. We previously characterized RPV metabolism after oral dosing and identified seven metabolites: four metabolites resulting from mono- or dioxygenation of the 2,6-dimethylphenyl ring itself or either of the two methyl groups located on that ring, one N-linked RPV glucuronide conjugate, and two O-linked RPV glucuronides produced via glucuronidation of mono- and dihydroxymethyl metabolites. However, as is true for most drugs, the metabolism of RPV after injection has yet to be reported. The phase II clinical trial HPTN 076 enrolled 136 HIV-uninfected women and investigated the safety and acceptability of long-acting injectable RPV for use in HIV pre-exposure prophylaxis. Through the analysis of plasma samples from 80 of these participants in the active product arm of the study, we were able to detect 2 metabolites after intramuscular injection of long-acting RPV, 2-hydroxymethyl-RPV, and RPV N-glucuronide. Of the total of 80 individuals, 72 participants exhibited detectable levels of 2-hydroxymethyl-RPV in plasma samples whereas RPV N-glucuronide was detectable in plasma samples of 78 participants. In addition, RPV N-glucuronide was detectable in rectal fluid, cervicovaginal fluid, and vaginal tissue. To investigate potential genetic variation in genes encoding enzymes relevant to RPV metabolism, we isolated genomic DNA and performed next-generation sequencing of CYP3A4, CYP3A5, UGT1A1 and UGT1A4. From these analyses, four missense variants were detected for CYP3A4 whereas one missense variant and one frameshift variant were detected for CYP3A5. A total of eight missense variants of UGT1A4 were detected, whereas two variants were detected for UGT1A1; however, these variants did not appear to account for the observed interindividual variability in metabolite levels. These findings provide insight into the metabolism of long-acting RPV and contribute to an overall understanding of metabolism after oral dosing versus injection. ClinicalTrials.gov Identifier: NCT02165202.

Identification of Novel UGT1A1 Variants Including UGT1A1 454C>A through the Genotyping of Healthy Participants of the HPTN 077 Study

Cabotegravir (CAB) is an integrase strand-transfer inhibitor of HIV that has proven effective for HIV treatment and prevention in a long-acting injectable formulation, typically preceded by an oral formulation lead-in phase. Previous in vitro studies have demonstrated that CAB is primarily metabolized via glucuronidation by uridine diphosphate glucuronosyltransferase (UGT) 1A1 and 1A9. In this study, we performed next-generation sequencing of genomic DNA isolated from the HPTN 077 participants to explore the variants within UGT1A1 and UGT1A9. Additionally, to enable correlation of UGT1A1 and UGT1A9 genotypes with plasma CAB-glucuronide levels, we quantified glucuronidated CAB following both oral administration of CAB and intramuscular injection of long-acting CAB. From these studies, 48 previously unreported variants of UGT1A1 and UGT1A9 were detected. Notably, 5/68 individuals carried a UGT1A1 454C>A variant that resulted in amino acid substitution P152T, and the use of in silico tools predicted a deleterious effect of the P152T substitution. Thus, the impact of this mutant on a range of UGT1A1 substrates was tested using a COS-7 cell-based assay. The glucuronide conjugates of CAB, dolutegravir, and raltegravir, were not formed in the COS-7 cells expressing the UGT1A1 P152T mutant. Further, formation of glucuronides of raloxifene and 7-ethyl-10-hydroxycamptothecin were reduced in the cells expressing the UGT1A1 P152T mutant. Using the same approach, we tested the activities of two UGT1A9 mutants, UGT1A9 H217Y and UGT1A9 R464G, and found that these mutations were tolerated and decreased function, respectively. These data provide insight into previously unreported genetic variants of UGT1A1 and UGT1A9.

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.

Parameterization of Microsomal and Cytosolic Scaling Factors: Methodological and Biological Considerations for Scalar Derivation and Validation

Mathematical models that can predict the kinetics of compounds have been increasingly adopted for drug development and risk assessment. Data for these models may be generated from in vitro experimental systems containing enzymes contributing to metabolic clearance, such as subcellular tissue fractions including microsomes and cytosol. Extrapolation from these systems is facilitated by common scaling factors, known as microsomal protein per gram (MPPG) and cytosolic protein per gram (CPPG). Historically, parameterization of MPPG and CPPG has employed the use of recovery factors, commonly benchmarked to cytochromes P450 which work well in some contexts, but could be problematic for other enzymes. Here, we propose absolute quantification of protein content and supplementary assays to evaluate microsomal/cytosolic purity that should be employed. Examples include calculation of microsomal latency by mannose-6-phosphatase activity and immunoblotting of subcellular fractions with fraction-specific markers. Further considerations include tissue source, as disease states can affect enzyme expression and activity, and the methodology used for scalar parameterization. Regional- and organ-specific expression of enzymes, in addition to differences in organ physiology, is another important consideration. Because most efforts have focused on the liver that is, for the most part, homogeneous, derived scalars may not capture the heterogeneity of other major tissues contributing to xenobiotic metabolism including the kidneys and small intestine. Better understanding of these scalars, and how to appropriately derive them from extrahepatic tissues can provide support to the inferences made with physiologically based pharmacokinetic modeling, increase its accuracy in characterizing in vivo drug pharmacokinetics, and improve confidence in go-no-go decisions for clinical trials.

In Vitro Characterization of Ertugliflozin Metabolism by UDP-Glucuronosyltransferase and Cytochrome P450 Enzymes

Ertugliflozin is primarily cleared through UDP-glucurosyltransferase (UGT)-mediated metabolism (86%) with minor oxidative clearance (12%). In vitro phenotyping involved enzyme kinetic characterization of UGTs or cytochrome P450 enzymes catalyzing formation of the major 3-O-β-glucuronide (M5c) and minor 2-O-β-glucuronide (M5a), monohydroxylated ertugliflozin (M1 and M3), and des-ethyl ertugliflozin (M2) metabolites in human liver microsomes (HLMs). Fractional clearance (fCL) from HLM intrinsic clearance (CLint) indicated a major role for glucuronidation (fCL 0.96; CLint 37 µl/min per milligram) versus oxidative metabolism (fCL 0.04; CLint 1.64 µl/min per milligram). Substrate concentration at half-maximal velocity (Km), maximal rate of metabolism (Vmax), and CLint for M5c and M5a formation were 10.8 µM, 375 pmol/min per milligram, and 34.7 µl/min per milligram and 41.7 µM, 94.9 pmol/min per milligram, and 2.28 µl/min per milligram, respectively. Inhibition of HLM CLint with 10 µM digoxin or tranilast (UGT1A9) and 3 µM 16β-phenyllongifolol (UGT2B7/UGT2B4) resulted in fraction metabolism (fm) estimates of 0.81 and 0.19 for UGT1A9 and UGT2B7/UGT2B4, respectively. Relative activity factor scaling of recombinant enzyme kinetics provided comparable fm for UGT1A9 (0.86) and UGT2B7 (0.14). Km and Vmax for M1, M2, and M3 formation ranged 73.0-93.0 µM and 24.3-116 pmol/min per milligram, respectively, and was inhibited by ketoconazole (M1, M2, and M3) and montelukast (M2). In summary, ertugliflozin metabolism in HLMs was primarily mediated by UGT1A9 (78%) with minor contributions from UGT2B7/UGT2B4 (18%), CYP3A4 (3.4%), CYP3A5 (0.4%), and CYP2C8 (0.16%). Considering higher ertugliflozin oxidative metabolism (fCL 0.12) obtained from human mass balance, human systemic clearance is expected to be mediated by UGT1A9 (70%), UGT2B7/UGT2B4 (16%), CYP3A4 (10%), CYP3A5 (1.2%), CYP2C8 (0.5%), and renal elimination (2%). SIGNIFICANCE STATEMENT: This manuscript describes the use of orthogonal approaches (i.e., enzyme kinetics, chemical inhibitors, and recombinant enzymes) to characterize the fraction of ertugliflozin metabolism through various UDP-glucuronosyltransferase (UGT) and cytochrome P450 (CYP) enzyme-mediated pathways. Phenotyping approaches routinely used to characterize CYP hepatic fractional metabolism (fm) to estimate specific enzymes contributing to overall systemic clearance were similarly applied for UGT-mediated metabolism. Defining the in vitro metabolic disposition and fm for ertugliflozin allows risk assessment when considering potential victim-based drug-drug interactions perpetrated by coadministered drugs.

Physiologically Based Pharmacokinetic Models of Probenecid and Furosemide to Predict Transporter Mediated Drug-Drug Interactions

Purpose

To provide whole-body physiologically based pharmacokinetic (PBPK) models of the potent clinical organic anion transporter (OAT) inhibitor probenecid and the clinical OAT victim drug furosemide for their application in transporter-based drug-drug interaction (DDI) modeling.

Methods

PBPK models of probenecid and furosemide were developed in PK-Sim®. Drug-dependent parameters and plasma concentration-time profiles following intravenous and oral probenecid and furosemide administration were gathered from literature and used for model development. For model evaluation, plasma concentration-time profiles, areas under the plasma concentration–time curve (AUC) and peak plasma concentrations (C_max) were predicted and compared to observed data. In addition, the models were applied to predict the outcome of clinical DDI studies.

Results

The developed models accurately describe the reported plasma concentrations of 27 clinical probenecid studies and of 42 studies using furosemide. Furthermore, application of these models to predict the probenecid-furosemide and probenecid-rifampicin DDIs demonstrates their good performance, with 6/7 of the predicted DDI AUC ratios and 4/5 of the predicted DDI C_max ratios within 1.25-fold of the observed values, and all predicted DDI AUC and C_max ratios within 2.0-fold.

Conclusions

Whole-body PBPK models of probenecid and furosemide were built and evaluated, providing useful tools to support the investigation of transporter mediated DDIs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11095-020-02964-z.

Metabolism and Pharmacokinetic Drug-Drug Interaction Profile of Vericiguat, A Soluble Guanylate Cyclase Stimulator: Results From Preclinical and Phase I Healthy Volunteer Studies

BACKGROUND

Vericiguat is a stimulator of soluble guanylate cyclase currently under investigation as a first-in-class therapy for worsening chronic heart failure (NCT02861534). Patients with heart failure often require polypharmacy because of comorbidities. Hence, understanding the clearance mechanisms, elimination, and potential for pharmacokinetic drug-drug interactions of vericiguat is important for dose recommendations in this patient population.

METHODS

Biotransformation and perpetrator properties of vericiguat were characterized in vitro using human hepatocytes, liver microsomes, and recombinant enzymes. This was complemented by a human mass balance study and ten drug-drug interaction studies in healthy volunteers wherein vericiguat was co-administered orally with omeprazole, magnesium/aluminum hydroxide, ketoconazole, rifampicin, mefenamic acid, midazolam, warfarin, digoxin, sacubitril/valsartan, aspirin, or sildenafil.

RESULTS

In the human mass balance study, mean total radioactivity recovered was 98.3% of the dose administered (53.1% and 45.2% excreted via urine and feces, respectively). The main metabolic pathway of vericiguat is glucuronidation via uridine diphosphate-glucuronosyltransferase 1A9 and 1A1. In vitro studies revealed a low risk of vericiguat acting as a perpetrator by inhibiting cytochrome P450s, uridine diphosphate-glucuronosyltransferase isoforms, or major transport proteins, or by inducing cytochrome P450s. These observations were supported by phase I drug-drug interaction studies. Phase I studies that assessed the propensity of vericiguat as a victim drug showed changes in the range that did not warrant recommendations for dose adjustment in phase III.

CONCLUSIONS

A low pharmacokinetic interaction potential of vericiguat was estimated from in vitro data and confirmed in vivo. Thus, vericiguat is suitable for a patient population with multiple comorbidities requiring polypharmacy.

Evidence-based strategies for the characterisation of human drug and chemical glucuronidation in vitro and UDP-glucuronosyltransferase reaction phenotyping

Enzymes of the UDP-glucuronosyltransferase (UGT) superfamily contribute to the elimination of drugs from almost all therapeutic classes. Awareness of the importance of glucuronidation as a drug clearance mechanism along with increased knowledge of the enzymology of drug and chemical metabolism has stimulated interest in the development and application of approaches for the characterisation of human drug glucuronidation in vitro, in particular reaction phenotyping (the fractional contribution of the individual UGT enzymes responsible for the glucuronidation of a given drug), assessment of metabolic stability, and UGT enzyme inhibition by drugs and other xenobiotics. In turn, this has permitted the implementation of in vitro - in vivo extrapolation approaches for the prediction of drug metabolic clearance, intestinal availability, and drug-drug interaction liability, all of which are of considerable importance in pre-clinical drug development. Indeed, regulatory agencies (FDA and EMA) require UGT reaction phenotyping for new chemical entities if glucuronidation accounts for ≥25% of total metabolism. In vitro studies are most commonly performed with recombinant UGT enzymes and human liver microsomes (HLM) as the enzyme sources. Despite the widespread use of in vitro approaches for the characterisation of drug and chemical glucuronidation by HLM and recombinant enzymes, evidence-based guidelines relating to experimental approaches are lacking. Here we present evidence-based strategies for the characterisation of drug and chemical glucuronidation in vitro, and for UGT reaction phenotyping. We anticipate that the strategies will inform practice, encourage development of standardised experimental procedures where feasible, and guide ongoing research in the field.

Epigenetic Regulation of Differentially Expressed Drug-Metabolizing Enzymes in Cancer

Drug metabolism is a biotransformation process of drugs, catalyzed by drug-metabolizing enzymes (DMEs), including phase I DMEs and phase II DMEs. The aberrant expression of DMEs occurs in the different stages of cancer. It can contribute to the development of cancer and lead to individual variations in drug response by affecting the metabolic process of carcinogen and anticancer drugs. Apart from genetic polymorphisms, which we know the most about, current evidence indicates that epigenetic regulation is also central to the expression of DMEs. This review summarizes differentially expressed DMEs in cancer and related epigenetic changes, including DNA methylation, histone modification, and noncoding RNAs. Exploring the epigenetic regulation of differentially expressed DMEs can provide a basis for implementing individualized and rationalized medication. Meanwhile, it can promote the development of new biomarkers and targets for the diagnosis, treatment, and prognosis of cancer. SIGNIFICANCE STATEMENT: This review summarizes the aberrant expression of DMEs in cancer and the related epigenetic regulation of differentially expressed DMEs. Exploring the epigenetic regulatory mechanism of DMEs in cancer can help us to understand the role of DMEs in cancer progression and chemoresistance. Also, it provides a basis for developing new biomarkers and targets for the diagnosis, treatment, and prognosis of cancer.

Analysis of Minor Acidic N-Glycans in Human Serum

Prior investigations by our research group focused on the method development for the simultaneous analysis of sulfated and phosphorylated glycans. Herein, the developed method was applied to analyze minor acidic N-glycans including sulfated and phosphorylated N-glycans in human serum. First, 2-aminobenzoic acid-labeled minor acidic N-glycans were enriched from the serum using a serotonin-immobilized column and were then separated into groups using hydrophilic interaction liquid chromatography, and analyzed by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Phosphorylated hybrid-type and sulfated bi-antennary N-glycans were detected in the serum. In addition, we observed that multiple types of glucuronidated N-glycans were present. These results indicate that the developed method is applicable to the analysis of glucuronidated as well as sulfated and phosphorylated N-glycans. It was also applied to the sera obtained from 17 healthy subjects and 15 pancreatic cancer patients, and the profiles of sulfated, phosphorylated, and glucuronidated N-glycans were compared. The expressed amount of glucuronidated N-glycans was significantly decreased in some pancreatic cancer patients. Numerous examples of the N-glycan analysis in human serum were reported, but phosphorylated and glucuronidated glycans were not investigated. The methods described herein allow the analysis of minor acidic glycans that are typically difficult to detect.

Novel molecular signatures and potential therapeutics in renal cell carcinomas: Insights from a comparative analysis of subtypes

Renal cell carcinomas (RCCs) are among the highest causes of cancer mortality. Although transcriptome profiling studies in the last decade have made significant molecular findings on RCCs, effective diagnosis and treatment strategies have yet to be achieved due to lack of adequate screening and comparative profiling of RCC subtypes. In this study, a comparative analysis was performed on RNA-seq based transcriptome data from each RCC subtype, namely clear cell RCC (KIRC), papillary RCC (KIRP) and kidney chromophobe (KICH), and mutual or subtype-specific reporter biomolecules were identified at RNA, protein, and metabolite levels by the integration of expression profiles with genome-scale biomolecular networks. This approach revealed already-known biomarkers in RCCs as well as novel biomarker candidates and potential therapeutic targets. Our findings also pointed out the incorporation of the molecular mechanisms of KIRC and KIRP, whereas KICH was shown to have distinct molecular signatures. Furthermore, considering the Dipeptidyl Peptidase 4 (DPP4) receptor as a potential therapeutic target specific to KICH, several drug candidates such as ZINC6745464 were identified through virtual screening of ZINC molecules. In this study, we reported valuable data for further experimental and clinical efforts, since the proposed molecules have significant potential for screening and therapeutic purposes in RCCs.

Molecular targeting of renal cell carcinoma by an oral combination

The 5-year survival rate of patients with metastatic renal cell carcinoma (mRCC) is <12% due to treatment failure. Therapeutic strategies that overcome resistance to modestly effective drugs for mRCC, such as sorafenib (SF), could improve outcome in mRCC patients. SF is terminally biotransformed by UDP-glucuronosyltransferase-1A9 (A9) mediated glucuronidation, which inactivates SF. In a clinical-cohort and the TCGA-dataset, A9 transcript and/or protein levels were highly elevated in RCC specimens and predicted metastasis and overall-survival. This suggested that elevated A9 levels even in primary tumors of patients who eventually develop mRCC could be a mechanism for SF failure. 4-methylumbelliferone (MU), a choleretic and antispasmodic drug, downregulated A9 and inhibited SF-glucuronidation in RCC cells. Low-dose SF and MU combinations inhibited growth, motility, invasion and downregulated an invasive signature in RCC cells, patient-derived tumor explants and/or endothelial-RCC cell co-cultures; however, both agents individually were ineffective. A9 overexpression made RCC cells resistant to the combination, while its downregulation sensitized them to SF treatment alone. The combination inhibited kidney tumor growth, angiogenesis and distant metastasis, with no detectable toxicity; A9-overexpressing tumors were resistant to treatment. With effective primary tumor control and abrogation of metastasis in preclinical models, the low-dose SF and MU combinations could be an effective treatment option for mRCC patients. Broadly, our study highlights how targeting specific mechanisms that cause the failure of “old” modestly effective FDA-approved drugs could improve treatment response with minimal alteration in toxicity profile.

Absorption, distribution, metabolism and excretion of molidustat in healthy participants

The absorption, distribution, metabolism and excretion of molidustat were investigated in healthy male participants. In study 1, a mass balance study, radiolabelled molidustat 25 mg (3.57 MBq) was administered as an oral solution (n = 4). Following rapid absorption, molidustat‐related radioactivity was predominantly distributed in plasma rather than in red blood cells. The total recovery of the administered radioactivity was 97.0%, which was mainly excreted renally (90.7%). Metabolite M‐1, produced by N‐glucuronidation, was the dominant component in plasma (80.2% of the area under the concentration‐time curve for total radioactivity) and was primarily excreted via urine (~85% of dose). Only minor amounts of unchanged molidustat were excreted in urine (~4%) and faeces (~6%). Study 2 investigated the absolute bioavailability and pharmacodynamics of molidustat (part 1, n = 12; part 2, n = 16). Orally administered molidustat immediate release tablets had an absolute bioavailability of 59%. Following intravenous administration (1, 5 and 25 mg), total body clearance of molidustat was 28.7‐34.5 L/h and volume of distribution at steady state was 39.3‐50.0 L. All doses of molidustat transiently elevated endogenous erythropoietin levels, irrespective of the route of administration. Molidustat was considered safe and well tolerated at the administered doses.

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.

Multi-Domain Networks Association for Biological Data Using Block Signed Graph Clustering

Multi-domain biological network association and clustering have attracted a lot of attention in biological data integration and understanding, which can provide a more global and accurate understanding of biological phenomenon. In many problems, different domains may have different cluster structures. Due to rapid growth of data collection from different sources, some domains may be strongly or weakly associated with the other domains. A key challenge is how to determine the degree of association among different domains, and to achieve accurate clustering results by data integration. In this paper, we propose an unsupervised learning approach for multi-domain network association by using block signed graph clustering. In particular, with consistency weights calculation, the proposed algorithm automatically identify domains relevant to each other strongly (or weakly) by assigning them larger (or smaller) weights. This approach not only significantly improve clustering accuracy but also understand multi-domain networks association. In each iteration of the proposed algorithm, we update consistency weights based on cluster structure of each domain, and then make use of different sets of eigenvectors to obtain different cluster structures in each domain. Experimental results on both synthetic data sets and real data sets (including neuron activity data and gene expression data) empirically demonstrate the effectiveness of the proposed algorithm in clustering performance and in domain association capability.

Emerging roles for UDP-glucuronosyltransferases in drug resistance and cancer progression

The best-known role of UDP-glucuronosyltransferase enzymes (UGTs) in cancer is the metabolic inactivation of drug therapies. By conjugating glucuronic acid to lipophilic drugs, UGTs impair the biological activity and enhance the water solubility of these agents, driving their elimination. Multiple clinical observations support an expanding role for UGTs as modulators of the drug response and in mediating drug resistance in numerous cancer types. However, accumulating evidence also suggests an influence of the UGT pathway on cancer progression. Dysregulation of the expression and activity of UGTs has been associated with the progression of several cancers, arguing for UGTs as possible mediators of oncogenic pathways and/or disease accelerators in a drug-naive context. The consequences of altered UGT activity on tumour biology are incompletely understood. They might be associated with perturbed levels of bioactive endogenous metabolites such as steroids and bioactive lipids that are inactivated by UGTs or through non-enzymatic mechanisms, thereby eliciting oncogenic signalling cascades. This review highlights the evidence supporting dual roles for the UGT pathway, affecting cancer progression and drug resistance. Pharmacogenomic testing of UGT profiles in patients and the development of therapeutic options that impair UGT actions could provide useful prognostic and predictive biomarkers and enhance the efficacy of anti-cancer drugs.

Sorafenib Activity and Disposition in Liver Cancer Does Not Depend on Organic Cation Transporter 1

Systemic therapy of advanced hepatocellular carcinoma (HCC) with the small-molecule multikinase inhibitor sorafenib is associated with large interindividual pharmacokinetic variability and unpredictable side effects potentially requiring dose reduction or treatment termination. Organic cation transporter (OCT1; gene SLC22A1) has been proposed as a clinical biomarker of HCC response. Because proof is lacking that OCT1 transports sorafenib, we used a combinatorial approach to define how OCT1 contributes to sorafenib transport. Overexpression of functional OCT1 protein in Xenopus laevis oocytes and mammalian cell lines did not facilitate sorafenib transport. Otherwise, sorafenib considerably accumulated in liver cancer cell lines despite negligible OCT1 mRNA and protein levels. Sorafenib pharmacokinetics was independent of OCT1 genotype in mice. Finally, SLC22A1 mRNA expression was significantly reduced by DNA methylation in The Cancer Genome Atlas HCC cohort. These results clearly demonstrate OCT1-independent cellular sorafenib uptake indicating that OCT1 is apparently not a valid biomarker of sorafenib response in HCC.

Factors Affecting Interindividual Variability of Hepatic UGT2B17 Protein Expression Examined Using a Novel Specific Monoclonal Antibody

Accurate quantification of the metabolic enzyme uridine diphospho-glucuronosyltransferase (UGT) UGT2B17 has been hampered by the high sequence identity with other UGT2B enzymes (as high as 94%) and by the lack of a specific antibody. Knowing the significance of the UGT2B17 pathway in drug and hormone metabolism and cancer, we developed a specific monoclonal antibody (EL-2B17mAb), initially validated by the lack of detection in liver microsomes of an individual carrying no UGT2B17 gene copy and in supersomes expressing UGT2B enzymes. Immunohistochemical detection in livers revealed strong labeling of bile ducts and variable labeling of hepatocytes. Expression levels assessed by immunoblotting were highly correlated to mass spectrometry-based quantification (r = 0.93), and three major expression patterns (absent, low, or high) were evidenced. Livers with very low expression were carriers of the functional rs59678213 G variant, located in the binding site for the transcription factor forkhead box A1 (FOXA1) of the UGT2B17 promoter. The highest level of expression was observed for individuals carrying at least one rs59678213 A allele. Multiple regression analysis indicated that the number of gene copies explained only 8% of UGT2B17 protein expression, 49% when adding rs59678213, reaching 54% when including sex. The novel EL-2B17mAb antibody allowed specific UGT2B17 quantification and exposed different patterns of hepatic expression. It further suggests that FOXA1 is a key driver of UGT2B17 expression in the liver. The availability of this molecular tool will help characterize the UGT2B17 level in various disease states and establish more precisely the contribution of the UGT2B17 enzyme to drug and hormone metabolism.

Mechanistic Assessment of Extrahepatic Contributions to Glucuronidation of Integrase Strand Transfer Inhibitors

Integrase strand transfer inhibitor (INSTI)-based regimens dominate initial human immunodeficiency virus treatment. Most INSTIs are metabolized predominantly via UDP-glucuronosyltransferases (UGTs). For drugs predominantly metabolized by UGTs, including INSTIs, in vitro data recovered from human liver microsomes (HLMs) alone often underpredict human oral clearance. While several factors may contribute, extrahepatic glucuronidation may contribute to this underprediction. Thus, we comprehensively characterized the kinetics for the glucuronidation of INSTIs (cabotegravir, dolutegravir, and raltegravir) using pooled human microsomal preparations from liver (HLMs), intestine (HIMs), and kidney (HKMs) tissues; human embryonic kidney 293 cells expressing individual UGTs; and recombinant UGTs. In vitro glucuronidation of cabotegravir (HLMs≈HKMs>>>HIMs), dolutegravir (HLMs>HIMs>>HKMs), and raltegravir (HLMs>HKMs>> HIMs) occurred in hepatic and extrahepatic tissues. The kinetic data from expression systems suggested the major enzymes in each tissue: hepatic UGT1A9 > UGT1A1 (dolutegravir and raltegravir) and UGT1A1 (cabotegravir), intestinal UGT1A3 > UGT1A8 > UGT1A1 (dolutegravir) and UGT1A8 > UGT1A1 (raltegravir), and renal UGT1A9 (dolutegravir and raltegravir). Enzymes catalyzing cabotegravir glucuronidation in the kidney and intestine could not be identified unequivocally. Using data from dolutegravir glucuronidation as a prototype, a "bottom-up" physiologically based pharmacokinetic model was developed in a stepwise approach and predicted dolutegravir oral clearance within 4.5-fold (hepatic data only), 2-fold (hepatic and intestinal data), and 32% (hepatic, intestinal, and renal data). These results suggest clinically meaningful glucuronidation of dolutegravir in tissues other than the liver. Incorporation of additional novel mechanistic and physiologic underpinnings of dolutegravir metabolism along with in silico approaches appears to be a powerful tool to accurately predict the clearance of dolutegravir from in vitro data.

knnAUC: an open-source R package for detecting nonlinear dependence between one continuous variable and one binary variable

BACKGROUND

Testing the dependence of two variables is one of the fundamental tasks in statistics. In this work, we developed an open-source R package (knnAUC) for detecting nonlinear dependence between one continuous variable X and one binary dependent variables Y (0 or 1).

RESULTS

We addressed this problem by using knnAUC (k-nearest neighbors AUC test, the R package is available at https://sourceforge.net/projects/knnauc/ ). In the knnAUC software framework, we first resampled a dataset to get the training and testing dataset according to the sample ratio (from 0 to 1), and then constructed a k-nearest neighbors algorithm classifier to get the yhat estimator (the probability of y = 1) of testy (the true label of testing dataset). Finally, we calculated the AUC (area under the curve of receiver operating characteristic) estimator and tested whether the AUC estimator is greater than 0.5. To evaluate the advantages of knnAUC compared to seven other popular methods, we performed extensive simulations to explore the relationships between eight different methods and compared the false positive rates and statistical power using both simulated and real datasets (Chronic hepatitis B datasets and kidney cancer RNA-seq datasets).

CONCLUSIONS

We concluded that knnAUC is an efficient R package to test non-linear dependence between one continuous variable and one binary dependent variable especially in computational biology area.