Human variability in isoform-specific UDP-glucuronosyltransferases: markers of acute and chronic exposure, polymorphisms and uncertainty factors

Association of UGT1A6 gene polymorphisms with sodium valproate-induced tremor in patients with epilepsy

BACKGROUND

Individual susceptibility to sodium valproate (VPA)-induced tremors may be due to genetic polymorphisms in the gene encoding the uridine diphosphate glucuronosyltransferase (UGT) enzyme, which affec the drug's clinical efficacy and cause toxic side effects. This study aimed to investigate the association between UGT1A6 polymorphisms and VPA-induced tremors in patients with epilepsy.

METHODS

In total, 128 patients with epilepsy were enrolled. Patients with epilepsy who received VPA were divided into tremor and non-tremor groups. Polymerase chain reaction-restriction fragment length polymorphism was used to investigate the genotype of UGT1A6 polymorphisms.

RESULTS

Carriers of the UGT1A6 A541G mutant genotype conferred a higher risk of tremor than wild-type carriers (odds ratio 2.128, P = 0.045). Logistic regression analysis showed that the A541G mutant genotype was a significant genetic risk factor for VPA-induced tremors. This suggests that individual susceptibility to VPA-induced tremors may result, at least partially, from genetic variation in UGT1A6 A541G.

CONCLUSIONS

Patients with epilepsy carrying the UGT1A6 A541G mutant genotype may have VPA-induced tremors, and early detection of this genotype will help guide the clinical individualizsation of VPA treatment.

Population-scale variability of the human UDP-glycosyltransferase gene family

Human UDP-glycosyltransferases (UGTs) are responsible for the glucuronidation of a wide variety of endogenous substrates and commonly prescribed drugs. Different genetic polymorphisms in UGT genes are implicated in interindividual differences in drug response and cancer risk. However, the genetic complexity beyond these variants has not been comprehensively assessed. We here leveraged whole-exome and whole-genome sequencing data from 141,456 unrelated individuals across 7 major human populations to provide a comprehensive profile of genetic variability across the human UGT gene family. Overall, 9666 exonic variants were observed of which 98.9% were rare. To interpret the functional impact of UGT missense variants, we developed a gene family-specific variant effect predictor. This algorithm identified a total of 1208 deleterious variants, most of which were found in African and South Asian populations. Structural analysis corroborated the predicted effects for multiple variations in substrate binding sites. Combined, our analyses provide a systematic overview of UGT variability, which can yield insights into interindividual differences in phase 2 metabolism and facilitate the translation of sequencing data into personalized predictions of UGT substrate disposition.

Applying Physiologically Based Pharmacokinetic Modeling to Interpret Carbamazepine's Nonlinear Pharmacokinetics and Its Induction Potential on Cytochrome P450 3A4 and Cytochrome P450 2C9 Enzymes

Carbamazepine (CBZ) is commonly prescribed for epilepsy and frequently used in polypharmacy. However, concerns arise regarding its ability to induce the metabolism of other drugs, including itself, potentially leading to the undertreatment of co-administered drugs. Additionally, CBZ exhibits nonlinear pharmacokinetics (PK), but the root causes have not been fully studied. This study aims to investigate the mechanisms behind CBZ's nonlinear PK and its induction potential on CYP3A4 and CYP2C9 enzymes. To achieve this, we developed and validated a physiologically based pharmacokinetic (PBPK) parent-metabolite model of CBZ and its active metabolite Carbamazepine-10,11-epoxide in GastroPlus®. The model was utilized for Drug-Drug Interaction (DDI) prediction with CYP3A4 and CYP2C9 victim drugs and to further explore the underlying mechanisms behind CBZ's nonlinear PK. The model accurately recapitulated CBZ plasma PK. Good DDI performance was demonstrated by the prediction of CBZ DDIs with quinidine, dolutegravir, phenytoin, and tolbutamide; however, with midazolam, the predicted/observed DDI AUClast ratio was 0.49 (slightly outside of the two-fold range). CBZ's nonlinear PK can be attributed to its nonlinear metabolism caused by autoinduction, as well as nonlinear absorption due to poor solubility. In further applications, the model can help understand DDI potential when CBZ serves as a CYP3A4 and CYP2C9 inducer.

Application of the Key Characteristics of Carcinogens to Bisphenol A

The ten key characteristics (KCs) of carcinogens are based on characteristics of known human carcinogens and encompass many types of endpoints. We propose that an objective review of the large amount of cancer mechanistic evidence for the chemical bisphenol A (BPA) can be achieved through use of these KCs. A search on metabolic and mechanistic data relevant to the carcinogenicity of BPA was conducted and web-based software tools were used to screen and organize the results. We applied the KCs to systematically identify, organize, and summarize mechanistic information for BPA, and to bring relevant carcinogenic mechanisms into focus. For some KCs with very large data sets, we utilized reviews focused on specific endpoints. Over 3000 studies for BPA from various data streams (exposed humans, animals, in vitro and cell-free systems) were identified. Mechanistic data relevant to each of the ten KCs were identified, with receptor-mediated effects, epigenetic alterations, oxidative stress, and cell proliferation being especially data rich. Reactive and bioactive metabolites are also associated with a number of KCs. This review demonstrates how the KCs can be applied to evaluate mechanistic data, especially for data-rich chemicals. While individual entities may have different approaches for the incorporation of mechanistic data in cancer hazard identification, the KCs provide a practical framework for conducting an objective examination of the available mechanistic data without a priori assumptions on mode of action. This analysis of the mechanistic data available for BPA suggests multiple and inter-connected mechanisms through which this chemical can act.

Uridine 5'-Diphospho-glucuronosyltransferase 1A3 (UGT1A3) Prediction of Hepatic Clearance of Organic Anion Transporting Polypeptide 1B3 (OATP1B3) Substrate Telmisartan by Glucuronidation Using In Vitro-In Vivo Extrapolation (IVIVE)

BACKGROUND AND OBJECTIVE

The prediction of pharmacokinetic parameters for drugs metabolised by cytochrome P450 enzymes has been the subject of active research for many years, while the application of in vitro-in vivo extrapolation (IVIVE) techniques for non-cytochrome P450 enzymes has not been thoroughly evaluated. There is still no established quantitative method for predicting hepatic clearance of drugs metabolised by uridine 5'-diphospho-glucuronosyltransferases (UGTs), not to mention those which undergo hepatic uptake. The objective of the study was to predict the human hepatic clearance for telmisartan based on in vitro metabolic stability and hepatic uptake results.

METHODS

Telmisartan was examined in liver systems, allowing to estimate intrinsic clearance (CLint, in vitro) based on the substrate disappearance rate with the use of liquid chromatography tandem mass spectrometry (LC-MS/MS) technique. Obtained CLint, in vitro values were corrected for corresponding unbound fractions. Prediction of human hepatic clearance was made from scaled unbound CLint, in vitro data with the use of the well-stirred model, and finally referenced to the literature value of observed clearance in humans, allowing determination of the essential scaling factors.

RESULTS

The in vitro scaled CLint, in vitro by UGT1A3 was assessed using three systems, human hepatocytes, liver microsomes, and recombinant enzymes. Obtained values were scaled and hepatic metabolism clearance was predicted, resulting in significant clearance underprediction. Utilization of the extended clearance concept (ECC) and hepatic uptake improved prediction of hepatic metabolism clearance. The scaling factors for hepatocytes, assessing the in vitro-in vivo difference, changed from sixfold difference to only twofold difference with the application of the ECC.

CONCLUSIONS

The study showed that taking into consideration hepatic uptake of a drug allows us to obtain satisfactory scaling factors, hence enabling the prediction of in vivo hepatic glucuronidation from in vitro data.

Pharmacokinetics, Pharmacodynamics, and Side Effects of Midazolam: A Review and Case Example

Midazolam, a short-acting benzodiazepine, is widely used to alleviate patient anxiety, enhance compliance, and aid in anesthesia. While its side effects are typically dose-dependent and manageable with vigilant perioperative monitoring, serious cardiorespiratory complications, including fatalities and permanent neurological impairment, have been documented. Prolonged exposure to benzodiazepines, such as midazolam, has been associated with neurological changes in infants. Despite attempts to employ therapeutic drug monitoring for optimal sedation dosing, its efficacy has been limited. Consequently, efforts are underway to identify alternative predictive markers to guide individualized dosing and mitigate adverse effects. Understanding these factors is crucial for determining midazolam's suitability for future administration, particularly after a severe adverse reaction. This article aims to elucidate the factors influencing midazolam's pharmacokinetics and pharmacodynamics, potentially leading to adverse events. Finally, a case study is presented to exemplify the complex investigation into the causative factors of midazolam-related adverse events.

Assessing city-wide pharmaceutical emissions to wastewater via modelling and passive sampling

With increasing numbers of chemicals used in modern society, assessing human and environmental exposure to them is becoming increasingly difficult. Recent advances in wastewater-based epidemiology enable valuable insights into public exposure to data-poor compounds. However, measuring all >26,000 chemicals registered under REACH is not just technically unfeasible but would also be incredibly expensive. In this paper, we argue that estimating emissions of chemicals based on usage data could offer a more comprehensive, systematic and efficient approach than repeated monitoring. Emissions of 29 active pharmaceutical ingredients (APIs) to wastewater were estimated for a medium-sized city in the Netherlands. Usage data was collected both on national and local scale and included prescription data, usage in health-care institutions and over-the-counter sales. Different routes of administration were considered as well as the excretion and subsequent in-sewer back-transformation of conjugates into respective parent compounds. Results suggest model-based emission estimation on a city-level is feasible and in good agreement with wastewater measurements obtained via passive sampling. Results highlight the need to include excretion fractions in the conceptual framework of emission estimation but suggest that the choice of an appropriate excretion fraction has a substantial impact on the resulting model performance.

Cannabinoid-Induced Stereoselective Inhibition of R-S-Oxazepam Glucuronidation: Cannabinoid-Oxazepam Drug Interactions

Benzodiazepines (BZDs) such as oxazepam are commonly prescribed depressant drugs known for their anxiolytic, hypnotic, muscle relaxant, and anticonvulsant effects and are frequently used in conjunction with other illicit drugs including cannabis. Oxazepam is metabolized in an enantiomeric-specific manner by glucuronidation, with S-oxazepam metabolized primarily by UGT2B15 and R-oxazepam glucuronidation mediated by both UGT 1A9 and 2B7. The goal of the present study was to evaluate the potential inhibitory effects of major cannabinoids, Δ9-tetrahydrocannabinol (THC) and cannabidiol (CBD), and major THC metabolites, 11-hydroxy-Δ9-tetrahydrocannabinol (11-OH-THC) and 11-nor-9-carboxy-Δ9-tetrahydrocannabinol (11-COOH-THC), on the UGT-mediated metabolism of R- and S-oxazepam. The cannabinoids and metabolites were screened as inhibitors of R- and S-oxazepam glucuronidation in microsomes isolated from HEK293 cells overexpressing individual UGT enzymes (rUGTs). The IC50 values were determined in human liver microsomes (HLM), human kidney microsomes (HKM), and rUGTs and utilized to estimate the nonspecific, binding-corrected Ki (Ki,u) values and predict the area under the concentration-time curve ratio (AUCR). The estimated Ki,u values observed in HLM for S- and R-oxazepam glucuronidation by CBD, 11-OH-THC, and THC were in the micromolar range (0.82 to 3.7 µM), with the Ki,u values observed for R-oxazepam glucuronidation approximately 2- to 5-fold lower as compared to those observed for S-oxazepam glucuronidation. The mechanistic static modeling predicted a potential clinically significant interaction between oral THC and CBD with oxazepam, with the AUCR values ranging from 1.25 to 3.45. These data suggest a pharmacokinetic drug-drug interaction when major cannabinoids like CBD or THC and oxazepam are concurrently administered.

Establishment of UGT1A1-knockout human iPS-derived hepatic organoids for UGT1A1-specific kinetics and toxicity evaluation

Uridine diphosphate glucuronosyltransferases (UGTs) are highly expressed in the liver and are involved in the metabolism of many drugs. In particular, UGT1A1 has a genetic polymorphism that causes decreased activity, leading to drug-induced hepatotoxicity. Therefore, an in vitro evaluation system that accurately predicts the kinetics of drugs involving UGT1A1 is required. However, there is no such evaluation system because of the absence of the UGT1A1-selective inhibitor. Here, using human induced pluripotent stem (iPS) cells, genome editing technology, and organoid technology, we generated UGT1A1-knockout human iPS hepatocyte-derived liver organoids (UGT1A1-KO i-HOs) as a model for UGT1A1-specific kinetics and toxicity evaluation. i-HOs showed higher gene expression of many drug-metabolizing enzymes including UGT1A1 than human iPS cell-derived hepatocyte-like cells (iPS-HLCs), suggesting that hepatic organoid technology improves liver functions. Wild-type (WT) i-HOs showed similar levels of UGT1A1 activity to primary human (cryopreserved) hepatocytes, while UGT1A1-KO i-HOs completely lost the activity. Additionally, to evaluate whether this model can be used to predict drug-induced hepatotoxicity, UGT1A1-KO i-HOs were exposed to SN-38, the active metabolite of irinotecan, an anticancer drug, and acetaminophen and confirmed that these cells could predict UGT1A1-mediated toxicity. Thus, we succeeded in generating model cells that enable evaluation of UGT1A1-specific kinetics and toxicity.

Bisphenol-A and phthalate metabolism in children with neurodevelopmental disorders

BACKGROUND

The etiology of autism spectrum (ASD) and Attention Deficit/Hyperactivity (ADHD) disorders are multifactorial. Epidemiological studies have shown associations with environmental pollutants, such as plasticizers. This study focused on two of these compounds, the Bisphenol-A (BPA) and Diethylhexyl Phthalate (DEHP). The major pathway for BPA and DEHP excretion is via glucuronidation. Glucuronidation makes insoluble substances more water-soluble allowing for their subsequent elimination in urine.

HYPOTHESIS

Detoxification of these two plasticizers is compromised in children with ASD and ADHD. Consequently, their tissues are more exposed to these two plasticizers.

METHODS

We measured the efficiency of glucuronidation in three groups of children, ASD (n = 66), ADHD (n = 46) and healthy controls (CTR, n = 37). The children were recruited from the clinics of Rutgers-NJ Medical School. A urine specimen was collected from each child. Multiple mass spectrometric analyses including the complete metabolome were determined and used to derive values for the efficiency of glucuronidation for 12 varied glucuronidation pathways including those for BPA and MEHP.

RESULTS

(1) Both fold differences and metabolome analyses showed that the three groups of children were metabolically different from each other. (2) Of the 12 pathways examined, only the BPA and DEHP pathways discriminated between the three groups. (3) Glucuronidation efficiencies for BPA were reduced by 11% for ASD (p = 0.020) and 17% for ADHD (p<0.001) compared to controls. DEHP showed similar, but not significant trends.

CONCLUSION

ASD and ADHD are clinically and metabolically different but share a reduction in the efficiency of detoxification for both BPA and DEHP with the reductions for BPA being statistically significant.

Drug-drug interactions that alter the exposure of glucuronidated drugs: Scope, UDP-glucuronosyltransferase (UGT) enzyme selectivity, mechanisms (inhibition and induction), and clinical significance

Drug-drug interactions (DDIs) arising from the perturbation of drug metabolising enzyme activities represent both a clinical problem and a potential economic loss for the pharmaceutical industry. DDIs involving glucuronidated drugs have historically attracted little attention and there is a perception that interactions are of minor clinical relevance. This review critically examines the scope and aetiology of DDIs that result in altered exposure of glucuronidated drugs. Interaction mechanisms, namely inhibition and induction of UDP-glucuronosyltransferase (UGT) enzymes and the potential interplay with drug transporters, are reviewed in detail, as is the clinical significance of known DDIs. Altered victim drug exposure arising from modulation of UGT enzyme activities is relatively common and, notably, the incidence and importance of UGT induction as a DDI mechanism is greater than generally believed. Numerous DDIs are clinically relevant, resulting in either loss of efficacy or an increased risk of adverse effects, necessitating dose individualisation. Several generalisations relating to the likelihood of DDIs can be drawn from the known substrate and inhibitor selectivities of UGT enzymes, highlighting the importance of comprehensive reaction phenotyping studies at an early stage of drug development. Further, rigorous assessment of the DDI liability of new chemical entities that undergo glucuronidation to a significant extent has been recommended recently by regulatory guidance. Although evidence-based approaches exist for the in vitro characterisation of UGT enzyme inhibition and induction, the availability of drugs considered appropriate for use as 'probe' substrates in clinical DDI studies is limited and this should be research priority.

Meat tenderization using acetaminophen (paracetamol/APAP): A review on deductive biochemical mechanisms, toxicological implications and strategies for mitigation

Meats consist of edible portions originating from domestic and wild animals. Meat's palatability and sensory accessibility largely depend on its tenderness to consumers. Although many factors influence meat tenderness, the cooking method cannot be neglected. Different chemical, mechanical, and natural means of meat tenderization have been considered healthy and safe for consumers. However, many households, food vendors, and bars in developing countries engage in the unhealthy use of acetaminophen (paracetamol/APAP) in meat tenderization due to the cost reduction it offers in the overall cooking process. Acetaminophen (paracetamol/APAP) is one of the most popular, relatively cheap, and ubiquitous over-the-counter drugs that induce serious toxicity challenges when misused. It is important to note that acetaminophen during cooking is hydrolyses into a toxic compound known as 4-aminophenol, which damages the liver and kidney and results in organ failure. Despite the reports on the increase in the use of acetaminophen for meat tenderizing in many web reports, there have not been any serious scientific publications on this subject. This study adopted classical/traditional methodology to review relevant literature retrieved from Scopus, PubMed, and ScienceDirect using relevant key terms (Acetaminophen, Toxicity, Meat tenderization, APAP, paracetamol, mechanisms) and Boolean operators (AND and OR). This paper provides in-depth information on the hazard and health implications of consuming acetaminophen tenderized meat via genetic and metabolic pathways deductions. Understanding these unsafe practices will promote awareness and mitigation strategies.

Drug Metabolism of Hepatocyte-like Organoids and Their Applicability in In Vitro Toxicity Testing

Emerging advances in the field of in vitro toxicity testing attempt to meet the need for reliable human-based safety assessment in drug development. Intrahepatic cholangiocyte organoids (ICOs) are described as a donor-derived in vitro model for disease modelling and regenerative medicine. Here, we explored the potential of hepatocyte-like ICOs (HL-ICOs) in in vitro toxicity testing by exploring the expression and activity of genes involved in drug metabolism, a key determinant in drug-induced toxicity, and the exposure of HL-ICOs to well-known hepatotoxicants. The current state of drug metabolism in HL-ICOs showed levels comparable to those of PHHs and HepaRGs for CYP3A4; however, other enzymes, such as CYP2B6 and CYP2D6, were expressed at lower levels. Additionally, EC50 values were determined in HL-ICOs for acetaminophen (24.0−26.8 mM), diclofenac (475.5−>500 µM), perhexiline (9.7−>31.5 µM), troglitazone (23.1−90.8 µM), and valproic acid (>10 mM). Exposure to the hepatotoxicants showed EC50s in HL-ICOs comparable to those in PHHs and HepaRGs; however, for acetaminophen exposure, HL-ICOs were less sensitive. Further elucidation of enzyme and transporter activity in drug metabolism in HL-ICOs and exposure to a more extensive compound set are needed to accurately define the potential of HL-ICOs in in vitro toxicity testing.

UGT2B7 c.-161C>T polymorphism frequency in Croatian population

Uridine diphosphate glucuronosyltransferase-2B7 (UGT2B7), enzyme responsible for the elimination of a number of xenobiotics through glucuronidation, is expressed in the gut, kidneys, intestines, and brain. However, data on the frequency of UGT2B7 polymorphisms in the Croatian population are limited. The aim of this study was to assess the frequency of the UGT2B7 c.-161C>T (rs7668258) polymorphism in the Croatian population and to compare it with reported frequencies in other populations. This polymorphism is in complete linkage disequilibrium with the UGT2B7 c.802C>T (UGT2B7*2, rs7439366) variant, which is important in clinical medicine. The study reports data of 501 participants from University Hospital Centre Zagreb. All data were collected and analysed retrospectively. Genotyping was performed by real-time polymerase chain reaction (PCR) using the TaqMan® Drug Metabolism Genotyping Assay for UGT2B7 c.-161C>T (rs7668258). We found that 120 (23.95 %) participants were carriers of the UGT2B7 c.-161CC genotype and 255 (50.9 %) were heterozygous carriers (UGT2B7 c.-161CT), while 126 (25.15 %) were homozygous carriers of the variant allele (UGT2B7 c.-161TT). The frequency of the variant UGT2B7 c.-161C>T allele in this study was T=0.506. The frequency of the UGT2B7 c.-161C>T allelic variants and genotypes in the Croatian population is similar to other European populations.

Dysregulated hepatic UDP-glucuronosyltransferases and flavonoids glucuronidation in experimental colitis

Glucuronidation catalyzed by UDP-glucuronosyltransferases (UGTs) is one of the most important phase II mechanisms, facilitating drug clearance via conjugation of glucuronic acid with polar groups of xenobiotics. Accumulating evidence suggests that IBDs impact drug disposition, but whether and how IBDs regulate UGTs and drug glucuronidation remains undefined. In this study, we aim to investigate the expression of UGTs and drug glucuronidation in experimental colitis. Given that glucuronidation occurs primarily in the liver, we analyzed the mRNA changes in hepatic UGTs with a DSS-induced mouse colitis model. Twelve UGTs were downregulated in the liver of colitis mice including UGT1A1 and UGT1A9 (two representative UGTs). Colitis in mice downregulated UGT1A1 and UGT1A9 in the liver but not in small intestine, colon, and kidney. We also established that the downregulation of UGTs was attributed to the disease itself rather than the DSS compound. Moreover, colitis-reduced UGT1A1 and UGT1A9 lead to dampened baicalein and puerarin glucuronidation. PXR was the only UGT regulator significantly downregulated in colitis mice, suggesting dysregulation of PXR is associated with the downregulation of UGT1A1 and UGT1A9, thereby potentially resulting in dysfunction of baicalein and puerarin glucuronidation. Collectively, we establish that UGTs and glucuronidation are dysregulated in colitis, and this effect may cause variation in drug responsiveness in IBDs.

Identification of naturally occurring inhibitors in Xian-Ling-Gu-Bao capsule against the glucuronidation of estrogens

Xian-Ling-Gu-Bao (XLGB) capsule, a well-known traditional Chinese medicine prescription, is widely used for the treatment of osteoporosis. It could significantly increase the levels of estrogen in ovariectomized rats and mice. However, this working mechanism has not been well elucidated. Considering that UDP-glucuronosyltransferase (UGT) enzymes are the important enzymes that inactivate and regulate estrogen activity in vivo, this study aimed to identify the bioactive compounds from XLGB against the glucuronidation of estrogens. First, thirty compounds were considered as candidate bioactive compounds based on our previous studies including pharmacological evaluation, chemical profiles, and metabolic profiles. Second, the characteristics of estrogen glucuronidation by uridine diphosphate glucuronic acid (UDPGA)-supplemented human liver microsomes (HLM), human intestine microsomes (HIM), and expressed UGT enzymes were determined, and the incubation systems of their key UGT enzymes were optimized. Then, inhibitory effects and mechanisms of XLGB and its main compounds toward the key UGT isozymes were further investigated. As a result, estrogen underwent efficient glucuronidation by HLM and HIM. UGT1A10, 1A1, and 2B7 were mainly responsible for the glucuronidation of estrone, β-estradiol, and estriol, respectively. For E1 and E2, UGT1A10 and 1A1 tended to mediate estrogen-3-O-glucuronidation, while UGT2B7 preferred catalyzing estrogen-16-O-glucuronidation. Furthermore, the incubation system for active UGT isoforms was optimized including Tris-HCl buffer, detergents, MgCl2 concentration, β-glucuronidase inhibitors, UDPGA concentration, protein concentration, and incubation time. Based on optimal incubation conditions, eleven, nine, and nine compounds were identified as the potent inhibitors for UGT1A10, 1A1, and 2B7, respectively (IC50 &lt; 4.97 μM and Ki &lt; 3.35 μM). Among them, six compounds (bavachin, isobavachin, isobavachalcone, neobavaisoflavone, corylifol A, and icariside II) simultaneously demonstrated potent inhibitory effects against these three active enzymes. Prenylated flavanols from Epimedium brevicornu Maxim., prenylated flavonoids from Psoralea corylifolia L., and salvianolic acids from Salvia miltiorrhiza Bge. were characterized as the most important and effective compounds. The identification of potent natural inhibitors of XLGB against the glucuronidation of estrogen laid an important foundation for the pharmacodynamic material basis.

A Computational Understanding of Inter-Individual Variability in CYP2D6 Activity to Investigate the Impact of Missense Mutations on Ochratoxin A Metabolism

Cytochrome P-450 (CYP) enzymes have a key role in the metabolism of xenobiotics of food origin, and their highly polymorphic nature concurs with the diverse inter-individual variability in the toxicokinetics (TK) and toxicodynamics (TD) of food chemicals. Ochratoxin A is a well-known mycotoxin which contaminates a large variety of food and is associated with food safety concerns. It is a minor substrate of CYP2D6, although the effects of CYP2D6 polymorphisms on its metabolism may be overlooked. Insights on this aspect would provide a useful mechanistic basis for a more science-based hazard assessment, particularly to integrate inter-individual differences in CYP2D6 metabolism. This work presents a molecular modelling approach for the analysis of mechanistic features with regard to the metabolic capacity of CYP2D6 variants to oxidise a number of substrates. The outcomes highlighted that a low-frequency CYP2D6 variant (CYP2D6*110) is likely to enhance ochratoxin A oxidation with possible consequences on TK and TD. It is therefore recommended to further analyse such TK and TD consequences. Generally speaking, we propose the identification of mechanistic features and parameters that could provide a semi-quantitative means to discriminate ligands based on the likelihood to undergo transformation by CYP2D6 variants. This would support the development of a fit-for-purpose pipeline which can be extended to a tool allowing for the bulk analysis of a large number of compounds. Such a tool would ultimately include inter-phenotypic differences of polymorphic xenobiotic-metabolising enzymes in the hazard assessment and risk characterisation of food chemicals.

PF-PLC micelles ameliorate cholestatic liver injury via regulating TLR4/MyD88/NF-κB and PXR/CAR/UGT1A1 signaling pathways in EE-induced rats

Paeoniflorin (PF) has a certain therapeutic effect on cholestasis liver injury. To further improve the bioavailability of PF and play its pharmacological role in liver protection, PF-phospholipid complex micelles (PF-PLC micelles) were prepared based on our previous research on PF-PLC. The protective effects of PF and PF-PLC micelles on cholestasis liver injury induced by 17α-ethynylestradiol (EE) were compared, and the possible mechanisms were further explored. Herein, we showed that PF-PLC micelles effectively improved liver function, alleviated liver pathological damage, and localized infiltration of inflammatory cells. Mechanism studies indicated that PF-PLC micelles treatment could suppress the TLR4/MyD88/NF-κB pathway, and further reduce the levels of pro-inflammatory factors. Meanwhile, our experimental results demonstrated that the beneficial effect of PF-PLC micelles on EE-induced cholestasis may be achieved by the upregulation of nuclear receptors and metabolic enzymes (PXR/CAR/UGT1A1). All these results indicate that PF-PLC micelles have great potential in the treatment of cholestatic liver disease.

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.

Investigating the interaction between organic anion transporter 1 and ochratoxin A: An in silico structural study to depict early molecular events of substrate recruitment and the impact of single point mutations

Organic anion transporters (OATs) belong to a subgroup of the solute carrier 22 transporter family. OATs have a central role in xenobiotic disposition affecting the toxicokinetics of its substrates and inter-individual differences in their expression, activity and function impact both toxicokinetics and toxicodynamics. Amongst OATs, OAT1 (solute carrier family 22 member 6) is involved in the urinary excretion of many xenobiotics bringing substrates into renal proximal tubular cells which can then be secreted across the apical membrane into the tubule lumen. The mycotoxin ochratoxin A has been shown to have a high affinity for OAT1, which is an important renal transporter involved in its urinary excretion. Nowadays, molecular modeling techniques are widely applied to assess protein-ligand interactions and may provide a tool to depict the mechanic of xenobiotic action be it toxicokinetics or toxicodynamics. This work provides a structured pipeline consisting of docking and molecular dynamic simulations to study OAT1-ligand interactions and the impact of OAT1 polymorphisms on such interactions. Such a computational structure-based analytical framework allowed to: i) model OAT1-substrate complex formation and depict the features correlating its sequence, structure and its capability to recruit substrates; and ii) investigate the impact of OAT1 missense mutations on substrate recruitment. Perspectives on applying such a structured pipeline to xenobiotic-metabolising enzymes are discussed.

Human variability in polymorphic CYP2D6 metabolism: Implications for the risk assessment of chemicals in food and emerging designer drugs

The major human cytochrome P450 CYP2D6 isoform enzyme plays important roles in the liver and in the brain with regards to xenobiotic metabolism. Xenobiotics as CYP2D6 substrates include a whole range of pharmaceuticals, pesticides and plant alkaloids to cite but a few. In addition, a number of endogenous compounds have been shown to be substrates of CYP2D6 including trace amines in the brain such as tyramine and 5-methoxytryptamine as well as anandamide and progesterone. Because of the polymorphic nature of CYP2D6, considerable inter-phenotypic and inter-ethnic differences in the pharmaco/toxicokinetics (PK/TK) and metabolism of CYP2D6 substrates exist with potential consequences on the pharmacology and toxicity of chemicals. Here, large extensive literature searches have been performed to collect PK data from published human studies for a wide range of pharmaceutical probe substrates and investigate human variability in CYP2D6 metabolism. The computed kinetic parameters resulted in the largest open source database, quantifying inter-phenotypic differences for the kinetics of CYP2D6 probe substrates in Caucasian and Asian populations, to date. The database is available in supplementary material (CYPD6 DB) and EFSA knowledge junction (DOI to added). Subsequently, meta-analyses using a hierarchical Bayesian model for markers of chronic oral exposure (oral clearance, area under the plasma concentration time curve) and acute oral exposure (maximum plasma concentration (Cmax) provided estimates of inter-phenotypic differences and CYP2D6-related uncertainty factors (UFs) for chemical risk assessment in Caucasian and Asian populations classified as ultra-rapid (UM), extensive (EMs), intermediate (IMs) and poor metabolisers (PMs). The model allowed the integration of inter-individual (i.e. inter-phenotypic and inter-ethnic), inter-compound and inter-study variability together with uncertainty in each PK parameter. Key findings include 1. Higher frequencies of PMs in Caucasian populations compared to Asian populations (>8% vs 1-2%) for which EM and IM were the most frequent phenotype. 2. Large inter-phenotypic differences in PK parameters for Caucasian EMs (coefficients of variation (CV) > 50%) compared with Caucasian PMs and Asian EMs and IMs (i.e CV < 40%). 3. Inter-phenotypic PK differences between EMs and PMs in Caucasian populations increase with the quantitative contribution of CYP2D6 for the metabolism (fm) for a range of substrates (fm_CYP2D6 range: 20-95% of dose) (range: 1-54) to a much larger extent than those for Asian populations (range: 1-4). 4. Exponential meta-regressions between Fm_CYP2D6 in EMs and inter-phenotypic differences were also shown to differ between Caucasian and Asian populations as well as CYP2D6-related UFs. Finally, implications of these results for the risk assessment of food chemicals and emerging designer drugs of public health concern, as CYP2D6 substrates, are highlighted and include the integration of in vitro metabolism data and CYP2D6-variability distributions for the development of quantitative in vitro in vivo extrapolation models.

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.

OpenCYP: An open source database exploring human variability in activities and frequencies of polymophisms for major cytochrome P-450 isoforms across world populations

The open source database "OpenCYP database" has been developed based on the results of extensive literature searches from the peer-reviewed literature. OpenCYP provides data on human variability on baseline of activities and polymophism frequencies for selected cytochrome P-450 isoforms (CYP1A2, CYP2A6, CYP2D6, CYP3A4/3A5 and CYP3A7) in healthy adult populations from world populations. CYP enzymatic activities were generally expressed as the metabolic ratio (MR) between an unchanged probe drug and its metabolite(s) in urine or plasma measured in healthy adults. Data on other age groups were very limited and fragmented, constituting an important data gap. Quantitative comparisons were often hampered by the different experimental conditions used. However, variability was quite limited for CYP1A2, using caffeine as a probe substrate, with a symmetrical distribution of metabolic activity values. For CYP3A4, human variability was dependent on the probe substrate itself with low variability when data considering the dextromethorphan/demethilathed metabolite MR were used and large variability when the urinary 6β-hydroxycortisol/cortisol ratio was used. The largest variability in CYP activity was shown for CYP2D6 activity, after oral dosing of dextromethorphan, for which genetic polymorphisms are well characterised and constitute a significant source of variability. It is foreseen that the OpenCYP database can contribute to promising tools to support the further development of QIVIVE and PBK models for human risk assessment of chemicals particularly when combined with information on isoform-specific content in cells using proteomic approaches.

Human Variability in Carboxylesterases and carboxylesterase-related Uncertainty Factors for Chemical Risk Assessment

Carboxylesterases (CES) are an important class of enzymes involved in the hydrolysis of a range of chemicals and show large inter-individual variability in vitro. An extensive literature search was performed to identify in vivo probe substrates for CES1 and CES2 together with their protein content and enzymatic activity. Human pharmacokinetic (PK) data on Cmax, clearance, and AUC were extracted from 89 publications and Bayesian meta-analysis was performed using a hierarchical model to derive CES-related variability distributions and related uncertainty factors (UF). The CES-related variability indicated that 97.5% of healthy adults are covered by the kinetic default UF (3.16), except for clopidogrel and dabigatran etexilate. Clopidogrel is metabolised for a small amount by the polymorphic CYP2C19, which can have an impact on the overall pharmacokinetics, while the variability seen for dabigatran etexilate might be due to differences in the absorption, since this can be influenced by food intake. The overall CES-related variability was moderate to high in vivo (<CV 50%), which might be due to possible polymorphism in the enzyme but also to the small sample size available per chemical. The presented CES-related variability can be used in combination with in vitro data to derive pathway-specific distributions.

Stratification of Volunteers According to Flavanone Metabolite Excretion and Phase II Metabolism Profile after Single Doses of 'Pera' Orange and 'Moro' Blood Orange Juices

Large interindividual variations in the biological response to citrus flavanones have been observed, and this could be associated with high variations in their bioavailability. The aim of this study was to identify the main determinants underlying interindividual differences in citrus flavanone metabolism and excretion. In a randomized cross-over study, non-obese and obese volunteers, aged 19-40 years, ingested single doses of Pera and Moro orange juices, and urine was collected for 24 h. A large difference in the recovery of the urinary flavanone phase II metabolites was observed, with hesperetin-sulfate and hesperetin-sulfo-O-glucuronide being the major metabolites. Subjects were stratified according to their total excretion of flavanone metabolites as high, medium, and low excretors, but the expected correlation with the microbiome was not observed at the genus level. A second stratification was proposed according to phase II flavanone metabolism, whereby participants were divided into two excretion groups: Profiles A and B. Profile B individuals showed greater biotransformation of hesperetin-sulfate to hesperetin-sulfo-O-glucuronide, as well as transformation of flavanone-monoglucuronide to the respective diglucuronides, suggestive of an influence of polymorphisms on UDP-glucuronosyltransferase. In conclusion, this study proposes a new stratification of volunteers based on their metabolic profiles. Gut microbiota composition and polymorphisms of phase II enzymes may be related to the interindividual variability of metabolism.

Inter-phenotypic differences in CYP2C9 and CYP2C19 metabolism: Bayesian meta-regression of human population variability in kinetics and application in chemical risk assessment

Quantifying variability in pharmacokinetics (PK) and toxicokinetics (TK) provides a science-based approach to refine uncertainty factors (UFs) for chemical risk assessment. In this context, genetic polymorphisms in cytochromes P450 (CYPs) drive inter-phenotypic differences and may result in reduction or increase in metabolism of drugs or other xenobiotics. Here, an extensive literature search was performed to identify PK data for probe substrates of the human polymorphic isoforms CYP2C9 and CYP2C19. Relevant data from 158 publications were extracted for markers of chronic exposure (clearance and area under the plasma concentration-time curve) and analysed using a Bayesian meta-regression model. Enzyme function (EF), driven by inter-phenotypic differences across a range of allozymes present in extensive and poor metabolisers (EMs and PMs), and fraction metabolised (Fm), were identified as exhibiting the highest impact on the metabolism. The Bayesian meta-regression model provided good predictions for such inter-phenotypic differences. Integration of population distributions for inter-phenotypic differences and estimates for EF and Fm allowed the derivation of CYP2C9- and CYP2C19-related UFs which ranged from 2.7 to 12.7, and were above the default factor for human variability in TK (3.16) for PMs and major substrates (Fm >60%). These results provide population distributions and pathway-related UFs as conservative in silico options to integrate variability in CYP2C9 and CYP2C19 metabolism using in vitro kinetic evidence and in the absence of human data. The future development of quantitative extrapolation models is discussed with particular attention to integrating human in vitro and in vivo PK or TK data with pathway-related variability for chemical risk assessment.

Human variability in glutathione-S-transferase activities, tissue distribution and major polymorphic variants: Meta-analysis and implication for chemical risk assessment

The input into the QIVIVE and Physiologically-Based kinetic and dynamic models of drug metabolising enzymes performance and their inter-individual differences significantly improve the modelling performance, supporting the development and integration of alternative approaches to animal testing. Bayesian meta-analyses allow generating and integrating statistical distributions with human in vitro metabolism data for quantitative in vitro-in vivo extrapolation. Such data are lacking on glutathione-S-transferases (GSTs). This paper reports for the first time results on the human variability of GST activities in healthy individuals, their tissue localisation and the frequencies of their major polymorphic variants by means of extensive literature search, data collection, data base creation and meta-analysis. A limited number of papers focussed on in vivo GST inter-individual differences in humans. Ex-vivo total GST activity without discriminating amongst isozymes is generally reported, resulting in a high inter-individual variability. The highest levels of cytosolic GSTs in humans are measured in the kidney, liver, adrenal glands and blood. The frequencies of GST polymorphisms for cytosolic isozymes in populations of different geographical ancestry were also presented. Bayesian meta-analyses to derive GST-related uncertainty factors provided uncertain estimates, due to the limited database. Considering the relevance of GST activities and their pivotal role in cellular adaptive response mechanisms to chemical stressors, further studies are needed to identify GST probe substrates for specific isozymes and quantify inter-individual differences.