Comparative modelling of the human UDP-glucuronosyltransferases: insights into structure and mechanism

Computational and artificial intelligence-based approaches for drug metabolism and transport prediction

Drug metabolism and transport, orchestrated by drug-metabolizing enzymes (DMEs) and drug transporters (DTs), are implicated in drug-drug interactions (DDIs) and adverse drug reactions (ADRs). Reliable and precise predictions of DDIs and ADRs are critical in the early stages of drug development to reduce the rate of drug candidate failure. A variety of experimental and computational technologies have been developed to predict DDIs and ADRs. Recent artificial intelligence (AI) approaches offer new opportunities for better predicting and understanding the complex processes related to drug metabolism and transport. We summarize the role of major DMEs and DTs, and provide an overview of current progress in computational approaches for the prediction of drug metabolism, transport, and DDIs, with an emphasis on AI including machine learning (ML) and deep learning (DL) modeling.

1H, 13C, 15N Backbone and sidechain chemical shift assignments of the C-terminal domain of human UDP-glucuronosyltransferase 2B17 (UGT2B17-C)

UDP-glucuronosyltransferases are the principal enzymes involved in the glucuronidation of metabolites and xenobiotics for physiological clearance in humans. Though glucuronidation is an indispensable process in the phase II metabolic pathway, UGT-mediated glucuronidation of most prescribed drugs (> 55%) and clinical evidence of UGT-associated drug resistance are major concerns for therapeutic development. While UGTs are highly conserved enzymes, they manifest unique substrate and inhibitor specificity which is poorly understood given the dearth of experimentally determined full-length structures. Such information is important not only to conceptualize their specificity but is central to the design of inhibitors specific to a given UGT in order to avoid toxicity associated with pan-UGT inhibitors. Here, we provide the 1H, 13C and 15N backbone (~ 90%) and sidechain (~ 62%) assignments for the C-terminal domain of UGT2B17, which can be used to determine the molecular binding sites of inhibitor and substrate, and to understand the atomic basis for inhibitor selectivity between UGT2B17 and other members of the UGT2B subfamily. Given the physiological relevance of UGT2B17 in the elimination of hormone-based cancer drugs, these assignments will contribute towards dissecting the structural basis for substrate specificity, selective inhibitor recognition and other aspects of enzyme activity with the goal of selectively overcoming glucuronidation-based drug resistance.

Machine learning and structure-based modeling for the prediction of UDP-glucuronosyltransferase inhibition

UDP-glucuronosyltransferases (UGTs) are responsible for 35% of the phase II drug metabolism. In this study, we focused on UGT1A1, which is a key UGT isoform. Strong inhibition of UGT1A1 may trigger adverse drug/herb-drug interactions, or result in disorders of endobiotic metabolism. Most of the current machine learning methods predicting the inhibition of drug metabolizing enzymes neglect protein structure and dynamics, both being essential for the recognition of various substrates and inhibitors. We performed molecular dynamics simulations on a homology model of the human UGT1A1 structure containing both the cofactor- (UDP-glucuronic acid) and substrate-binding domains to explore UGT conformational changes. Then, we created models for the prediction of UGT1A1 inhibitors by integrating information on UGT1A1 structure and dynamics, interactions with diverse ligands, and machine learning. These models can be helpful for further prediction of drug-drug interactions of drug candidates and safety treatments.

•

UGTs are responsible for 35% of the phase II drug metabolism reactions

•

We created machine learning models for prediction of UGT1A1 inhibitors

•

Our simulations suggested key residues of UGT1A1 involved in the substrate binding

Predicting Regioselectivity of AO, CYP, FMO, and UGT Metabolism Using Quantum Mechanical Simulations and Machine Learning

Unexpected metabolism in modification and conjugation phases can lead to the failure of many late-stage drug candidates or even withdrawal of approved drugs. Thus, it is critical to predict the sites of metabolism (SoM) for enzymes, which interact with drug-like molecules, in the early stages of the research. This study presents methods for predicting the isoform-specific metabolism for human AOs, FMOs, and UGTs and general CYP metabolism for preclinical species. The models use semi-empirical quantum mechanical simulations, validated using experimentally obtained data and DFT calculations, to estimate the reactivity of each SoM in the context of the whole molecule. Ligand-based models, trained and tested using high-quality regioselectivity data, combine the reactivity of the potential SoM with the orientation and steric effects of the binding pockets of the different enzyme isoforms. The resulting models achieve κ values of up to 0.94 and AUC of up to 0.92.

Discovery and Characterization of the Key Constituents in Ginkgo biloba Leaf Extract With Potent Inhibitory Effects on Human UDP-Glucuronosyltransferase 1A1

Human UDP-glucuronosyltransferase 1A1 (hUGT1A1) is one of the most essential phase II enzymes in humans. Dysfunction or strong inhibition of hUGT1A1 may result in hyperbilirubinaemia and clinically relevant drug/herb-drug interactions (DDIs/HDIs). Recently, a high-throughput fluorescence-based assay was constructed by us to find the compounds/herbal extracts with strong inhibition against intracellular hUGT1A1. Following screening of over one hundred of herbal products, the extract of Ginkgo biloba leaves (GBL) displayed the most potent hUGT1A1 inhibition in HeLa-UGT1A1 cells (Hela cells overexpressed hUGT1A1). Further investigations demonstrated that four biflavones including bilobetin, isoginkgetin, sciadopitysin and ginkgetin, are key constituents responsible for hUGT1A1 inhibition in living cells. These biflavones potently inhibit hUGT1A1 in both human liver microsomes (HLM) and living cells, with the IC₅₀ values ranging from 0.075 to 0.41 μM in living cells. Inhibition kinetic analyses and docking simulations suggested that four tested biflavones potently inhibit hUGT1A1-catalyzed NHPN-O-glucuronidation in HLM via a mixed inhibition manner, showing the K_i values ranging from 0.07 to 0.74 μM. Collectively, our findings uncover the key constituents in GBL responsible for hUGT1A1 inhibition and decipher their inhibitory mechanisms against hUGT1A1, which will be very helpful for guiding the rational use of GBL-related herbal products in clinical settings.

Inhibition of Human UGT1A1-Mediated Bilirubin Glucuronidation by the Popular Flavonoids Baicalein, Baicalin and Hyperoside is responsible for Herbs (Shuang-huang-lian) -Induced Jaundice

Bilirubin-related adverse drug reactions (ADRs) or malady (e.g., jaundice) induced by some herbs rich in certain flavonoids have been widely reported. However, the causes and mechanisms of the ADRs are not well understood. The aim of this paper was to explore the mechanism of Shuang-huang-lian injections (SHL) and its major constituents-induced jaundice via inhibiting human UDP-glucuronosyltransferases1A1 (hUGT1A1)-mediated bilirubin glucuronidation. The inhibitory effects of SHL and its major constituents in the herbal medicine including baicalein (BAI), baicalin (BA) and hyperoside (HYP) on bilirubin glucuronidation were investigated. This study indicated that the average formation rates of bilirubin glucuronides (i.e., BMG1, BMG2, BDG) displayed significant differences (P <0.05), specially, the formation of mono-glucuronides (BMGs) was favored regardless whether an inhibitor was absent or presence. SHL, BAI, BA and HYP dose-dependently inhibit bilirubin glucuronidation, showing the IC₅₀ values against total bilirubin glucuronidation (TBG) were in the range of (7.69 {plus minus} 0.94) μg/mL - (37.09 {plus minus} 2.03) μg/mL, (4.51 {plus minus} 0.27) μM - (20.84 {plus minus} 1.99) μM, (22.36 {plus minus} 5.74) μM - (41.35 {plus minus} 2.40) μM, and (15.16 {plus minus} 1.12) μM - (42.80 {plus minus} 2.63) μM for SHL, BAI, BA, and HYP, respectively. Both inhibition kinetics assays and molecular docking simulations suggested that SHL, BAI, BA, and HYP significantly inhibited hUGT1A1-mediated bilirubin glucuronidation via a mixed-type inhibition. Collectively, some naturally occurring flavonoids (BAI, BA and HYP) in SHL have been identified as the inhibitors against hUGT1A1-mediated bilirubin glucuronidation, which well-explains the bilirubin-related ADRs or malady triggered by SHL in clinical settings. Significance Statement Herbal products and their components (e.g., flavonoids), which been widely used in the whole world, may cause liver injury. As a commonly used herbal products rich in flavonoids, Shuang-huang-lian injections (SHL), easily lead to symptoms of liver injury (e.g., jaundice) owing to significant inhibition of hUGT1A1-mediated bilirubin glucuronidation by its flavonoid components (i.e., baicalein, baicalin, hyperoside). Herbs-induced bilirubin-related ADRs and its associated clinical significance should be seriously considered.

Methylophiopogonanone A is a naturally occurring broad-spectrum inhibitor against human UDP-glucuronosyltransferases: Inhibition behaviours and implication in herb-drug interactions

Methylophiopogonanone A (MOA) is an abundant homoisoflavonoid in the Chinese herb Ophiopogonis Radix. Recent investigations revealed that MOA inhibited several human cytochrome P450 enzymes (CYPs) and stimulated OATP1B1. However, the inhibitory effects of MOA on phase II drug-metabolizing enzymes, such as human UDP-glucuronosyltransferases (hUGTs), have not been well investigated. Herein, the inhibition potentials of MOA on hUGTs were assessed. The results clearly demonstrated that MOA dose-dependently inhibited all tested hUGTs including UGT1A1 (IC₅₀  = 1.23 μM), one of the most important detoxification enzymes in humans. Further investigations showed that MOA strongly inhibited UGT1A1-catalysed NHPH-O-glucuronidation in a range of biological settings including hUGT1A1, human liver microsomes (HLM) and HeLa cells overexpressing UGT1A1. Inhibition kinetic analyses demonstrated that MOA competitively inhibited UGT1A1-catalysed NHPH-O-glucuronidation in both hUGT1A1 and HLM, with K_i values of 0.52 and 1.22 μM, respectively. Collectively, our findings expanded knowledge of the interactions between MOA and human drug-metabolizing enzymes, which would be very helpful for guiding the use of MOA-related herbal products in clinical settings.

Predicting reactivity to drug metabolism: beyond P450s-modelling FMOs and UGTs

We present a study based on density functional theory calculations to explore the rate limiting steps of product formation for oxidation by Flavin-containing Monooxygenase (FMO) and glucuronidation by the UDP-glucuronosyltransferase (UGT) family of enzymes. FMOs are responsible for the modification phase of metabolism of a wide diversity of drugs, working in conjunction with Cytochrome P450 (CYP) family of enzymes, and UGTs are the most important class of drug conjugation enzymes. Reactivity calculations are important for prediction of metabolism by CYPs and reactivity alone explains around 70-85% of the experimentally observed sites of metabolism within CYP substrates. In the current work we extend this approach to propose model systems which can be used to calculate the activation energies, i.e. reactivity, for the rate-limiting steps for both FMO oxidation and glucuronidation of potential sites of metabolism. These results are validated by comparison with the experimentally observed reaction rates and sites of metabolism, indicating that the presented models are suitable to provide the basis of a reactivity component within generalizable models to predict either FMO or UGT metabolism.

Establishment of the experimental procedure for prediction of conjugation capacity in mutant UGT1A1

UDP-glucuronosyltransferase 1A1 (UGT1A1) is an enzyme that is found in the endoplasmic reticulum membrane and can reportedly have a large number of amino acid substitutions that result in the reduction of glucuronidation capacity. For example, adverse drug reactions when patients receive CPT-11 (irinotecan) such as in cancer chemotherapy are caused by amino acid substitutions in UGT1A1. We previously found that the extent of the docking when the hydroxyl residue of bilirubin was oriented toward UDP-glucuronic acid correlated with in vitro conjugation capacity. In this study, we analyzed the conformation of mutant UGT1A1s by means of structural optimization with water and lipid bilayers instead of the optimization in vacuo that we used in our previous study. We then derived a mathematical model that can predict the conjugation capacities of mutant UGT1A1s by using results of substrate docking in silico and results of in vitro analysis of glucuronidation of acetaminophen and 17β-estradiol by UGT1A1s. This experimental procedure showed that the in silico conjugation capacities of other mutant UGT1A1s with bilirubin or SN-38 were similar to reported in vitro conjugation capacities. Our results suggest that this experimental procedure described herein can correctly predict the conjugation capacities of mutant UGT1A1s and any substrate.

The evolution of UDP-glycosyl/glucuronosyltransferase 1E (UGT1E) genes in bird lineages is linked to feeding habits but UGT2 genes is not

UDP-glycosyltransferase (UGT) catalyzes the transfer of glycosyl groups (e.g., glucuronic acid) to exogenous or endogenous chemicals and plays an important role in conjugation reactions. In vertebrates, UGT genes are divided into 5 families: UGT1, UGT2, UGT3, UGT5, and UGT8. Among these UGT enzymes, UGT1 and UGT2 enzymes are known to be important xenobiotic metabolizing enzymes in mammals. However, little is known about UGT1 and UGT2 genes in avian species. In this study, we therefore aimed to classify avian UGT1 and UGT2 genes based on their evolutionary relationships. We also investigated the association between UGT molecular evolution and ecological factors, specifically feeding habits, habitat, and migration. By examining the genomes of 43 avian species with differing ecology, we showed that avian UGT1E genes are divided into 6 groups and UGT2 genes into 3 groups. Correlations between UGT gene count and ecological factors suggested that the number of UGT1E genes is decreasing in carnivorous species. Estimates of selection pressure also support the hypothesis that diet influenced avian UGT1E gene evolution, similar to mammalian UGT1A and UGT2B genes.

Chemical Probes for Human UDP-Glucuronosyltransferases: A Comprehensive Review

UGTs play crucial roles in the metabolism and detoxification of both endogenous and xenobiotic compounds. The key roles of UGTs in human health have garnered great interest in the design and development of specific probes for human UGTs. However, in contrast to other human enzymes, the probe substrates for human UGTs are rarely reported, owing to the highly overlapping substrate specificities of UGTs and the lack of the integrated crystal structures of UGTs. Over the past decades, many efforts are made to develop specific probe substrates for UGTs and use them in both basic research and drug discovery. This review focuses on recent progress in the development of probe substrates for UGTs and their biomedical applications. A long list of chemical probes for UGTs, including non-fluorescent and fluorescent probes along with their structural information and kinetic parameters, are prepared and analyzed. Additionally, challenges and future directions in this field are highlighted in the final section. All information and knowledge presented in this review provide practical tools/methods for measuring UGT activities in complex biological samples, which will be very helpful for rapid screening and characterization of UGT modulators, and for exploring the relevance of UGT enzymes to human diseases.

Revisiting the Latency of Uridine Diphosphate-Glucuronosyltransferases (UGTs)-How Does the Endoplasmic Reticulum Membrane Influence Their Function?

Uridine diphosphate-glucuronosyltransferases (UGTs) are phase 2 conjugation enzymes mainly located in the endoplasmic reticulum (ER) of the liver and many other tissues, and can be recovered in artificial ER membrane preparations (microsomes). They catalyze glucuronidation reactions in various aglycone substrates, contributing significantly to the body's chemical defense mechanism. There has been controversy over the last 50 years in the UGT field with respect to the explanation for the phenomenon of latency: full UGT activity revealed by chemical or physical disruption of the microsomal membrane. Because latency can lead to inaccurate measurements of UGT activity in vitro, and subsequent underprediction of drug clearance in vivo, it is important to understand the mechanisms behind this phenomenon. Three major hypotheses have been advanced to explain UGT latency: compartmentation, conformation, and adenine nucleotide inhibition. In this review, we discuss the evidence behind each hypothesis in depth, and suggest some additional studies that may reveal more information on this intriguing phenomenon.

Inter-isoform Hetero-dimerization of Human UDP-Glucuronosyltransferases (UGTs) 1A1, 1A9, and 2B7 and Impacts on Glucuronidation Activity

Human UDP-glucuronosyltransferases (UGTs) play a pivotal role in phase II metabolism by catalyzing the glucuronidation of endobiotics and xenobiotics. The catalytic activities of UGTs are highly impacted by both genetic polymorphisms and oligomerization. The present study aimed to assess the inter-isoform hetero-dimerization of UGT1A1, 1A9, and 2B7, including the wild type (1A1*1, 1A9*1, and 2B7*1) and the naturally occurring (1A1*1b, 1A9*2/*3/*5, and 2B7*71S/*2/*5) variants. The related enzymes were double expressed in Bac-to-Bac systems. The fluorescence resonance energy transfer (FRET) technique and co-immunoprecipitation (Co-IP) revealed stable hetero-dimerization of UGT1A1, 1A9, and 2B7 allozymes. Variable FRET efficiencies and donor-acceptor distances suggested that genetic polymorphisms resulted in altered affinities to the target protein. In addition, the metabolic activities of UGTs were differentially altered upon hetero-dimerization via double expression systems. Moreover, protein interactions also changed the regioselectivity of UGT1A9 for querectin glucuronidation. These findings provide in-depth understanding of human UGT dimerization as well as clues for complicated UGT dependent metabolism in humans.

Structure and Protein-Protein Interactions of Human UDP-Glucuronosyltransferases

Mammalian UDP-glucuronosyltransferases (UGTs) catalyze the transfer of glucuronic acid from UDP-glucuronic acid to various xenobiotics and endobiotics. Since UGTs comprise rate-limiting enzymes for metabolism of various compounds, co-administration of UGT-inhibiting drugs and genetic deficiency of UGT genes can cause an increased blood concentration of these compounds. During the last few decades, extensive efforts have been made to advance the understanding of gene structure, function, substrate specificity, and inhibition/induction properties of UGTs. However, molecular mechanisms and physiological importance of the oligomerization and protein–protein interactions of UGTs are still largely unknown. While three-dimensional structures of human UGTs can be useful to reveal the details of oligomerization and protein–protein interactions of UGTs, little is known about the protein structures of human UGTs due to the difficulty in solving crystal structures of membrane-bound proteins. Meanwhile, soluble forms of plant and bacterial UGTs as well as a partial domain of human UGT2B7 have been crystallized and enabled us to predict three-dimensional structures of human UGTs using a homology-modeling technique. The homology-modeled structures of human UGTs do not only provide the detailed information about substrate binding or substrate specificity in human UGTs, but also contribute with unique knowledge on oligomerization and protein–protein interactions of UGTs. Furthermore, various in vitro approaches indicate that UGT-mediated glucuronidation is involved in cell death, apoptosis, and oxidative stress as well. In the present review article, recent understandings of UGT protein structures as well as physiological importance of the oligomerization and protein–protein interactions of human UGTs are discussed.

In silico pharmacology: Drug membrane partitioning and crossing

Over the past decade, molecular dynamics (MD) simulations have become particularly powerful to rationalize drug insertion and partitioning in lipid bilayers. MD simulations efficiently support experimental evidences, with a comprehensive understanding of molecular interactions driving insertion and crossing. Prediction of drug partitioning is discussed with respect to drug families (anesthetics; β-blockers; non-steroidal anti-inflammatory drugs; antioxidants; antiviral drugs; antimicrobial peptides). To accurately evaluate passive permeation coefficients turned out to be a complex theoretical challenge; however the recent methodological developments based on biased MD simulations are particularly promising. Particular attention is paid to membrane composition (e.g., presence of cholesterol), which influences drug partitioning and permeation. Recent studies concerning in silico models of membrane proteins involved in drug transport (influx and efflux) are also reported here. These studies have allowed gaining insight in drug efflux by, e.g., ABC transporters at an atomic resolution, explicitly accounting for the mandatory forces induced by the surrounded lipid bilayer. Large-scale conformational changes were thoroughly analyzed.

Substrate selectivity of human intestinal UDP-glucuronosyltransferases (UGTs): in silico and in vitro insights

The current drug development process aims to produce safe, effective drugs within a reasonable time and at a reasonable cost. Phase II metabolism (glucuronidation) can affect drug action and pharmacokinetics to a considerable extent and so its studies and prediction at initial stages of drug development are very imperative. Extensive glucuronidation is an obstacle to oral bioavailability because the first-pass glucuronidation [or premature clearance by UDP-glucuronosyltransferases (UGTs)] of orally administered agents frequently results in poor oral bioavailability and lack of efficacy. Modeling of new chemical entities/drugs for UGTs and their kinetic data can be useful in understanding the binding patterns to be used in the design of better molecules. This review concentrates on first-pass glucuronidation by intestinal UGTs, including their topology, expression profile, and pharmacogenomics. In addition, recent advances are discussed with respect to substrate selectivity at the binding pocket, structural requirements, and mechanism of enzyme actions.

Genetic variations and haplotype diversity of the UGT1 gene cluster in the Chinese population

Vertebrates require tremendous molecular diversity to defend against numerous small hydrophobic chemicals. UDP-glucuronosyltransferases (UGTs) are a large family of detoxification enzymes that glucuronidate xenobiotics and endobiotics, facilitating their excretion from the body. The UGT1 gene cluster contains a tandem array of variable first exons, each preceded by a specific promoter, and a common set of downstream constant exons, similar to the genomic organization of the protocadherin (Pcdh), immunoglobulin, and T-cell receptor gene clusters. To assist pharmacogenomics studies in Chinese, we sequenced nine first exons, promoter and intronic regions, and five common exons of the UGT1 gene cluster in a population sample of 253 unrelated Chinese individuals. We identified 101 polymorphisms and found 15 novel SNPs. We then computed allele frequencies for each polymorphism and reconstructed their linkage disequilibrium (LD) map. The UGT1 cluster can be divided into five linkage blocks: Block 9 (UGT1A9), Block 9/7/6 (UGT1A9, UGT1A7, and UGT1A6), Block 5 (UGT1A5), Block 4/3 (UGT1A4 and UGT1A3), and Block 3′ UTR. Furthermore, we inferred haplotypes and selected their tagSNPs. Finally, comparing our data with those of three other populations of the HapMap project revealed ethnic specificity of the UGT1 genetic diversity in Chinese. These findings have important implications for future molecular genetic studies of the UGT1 gene cluster as well as for personalized medical therapies in Chinese.

First-pass metabolism via UDP-glucuronosyltransferase: a barrier to oral bioavailability of phenolics

Glucuronidation mediated by UDP-glucuronosyltransferases (UGTs) is a significant metabolic pathway that facilitates efficient elimination of numerous endobiotics and xenobiotics, including phenolics. UGT genetic deficiency and polymorphisms or inhibition of glucuronidation by concomitant use of drugs are associated with inherited physiological disorders or drug-induced toxicities. Moreover, extensive glucuronidation can be a barrier to oral bioavailability as the first-pass glucuronidation (or premature clearance by UGTs) of orally administered agents usually results in the poor oral bioavailability and lack of efficacies. This review focused on the first-pass glucuronidation of phenolics including natural polyphenols and pharmaceuticals. The complexity of UGT-mediated metabolism of phenolics is highlighted with species-, gender-, organ- and isoform-dependent specificity, as well as functional compensation between UGT1A and 2B subfamily. In addition, recent advances are discussed with respect to the mechanisms of enzymatic actions, including the important properties such as binding pocket size and phosphorylation requirements.

The crystal structure of human UDP-glucuronosyltransferase 2B7 C-terminal end is the first mammalian UGT target to be revealed: the significance for human UGTs from both the 1A and 2B families

Human UDP-glucuronosyltransferases (EC 2.4.1.17) (UGTs) are major phase II metabolism enzymes that detoxify a multitude of endo- and xenobiotics through the covalent addition of a glucuronic acid moiety. UGTs are promiscuous enzymes that regulate the levels of numerous important endobiotics in a range of tissues, and inactivate most therapeutic compounds in concert with phase I enzymes. In spite of the importance of these enzymes, we have only a limited understanding of the molecular mechanisms governing their substrate specificity and catalytic activity. Until recently, no three-dimensional structural information was available for any mammalian UGT. The 1.8-å resolution apo crystal structure of the UDP-glucuronic acid binding domain of human UGT2B7 (2B7CT) is the only structure of a mammalian UGT target determined to date. In this review, we summarize what has been learned about human UGT function from the analysis of this and other related glycosyltransferase (GT) crystal structures.

A molecular model of the human UDP-glucuronosyltransferase 1A1, its membrane orientation, and the interactions between different parts of the enzyme

The vertebrate UDP-glucuronosyltransferases (UGTs) are membrane-bound enzymes of the endoplasmic reticulum that process both endogenous and exogenous substrates. The human UGTs are well known biologically, but biophysical understanding is scarce, largely because of problems in purification. The one resolved crystal structure covers the C-terminal domain of the human UGT2B7. Here, we present a homology model of the complete monomeric human UGT1A1, the enzyme that catalyzes bilirubin glucuronidation. The enzyme can be seen as composed of four different domains: two large ones, the N- and C-terminal domains, and two small ones, the "envelope" helices and the transmembrane segment that includes the cytoplasmic tail. The hydrophobic core of the N-terminal domain and the two envelope helices that connect the large domains are shown to be structurally well conserved even among distant homologs and can thus be modeled with good certainty according to plant and bacterial structures. We consider alternative solutions for the highly variable N-terminal regions that probably contribute to substrate binding. The bilirubin binding site, known pathological mutations in UGT1A1, and other specific residues have been examined in the context of the model with regard to available experimental data. A putative orientation of the protein relative to the membrane has been derived from the location of predicted N-glycosylation sites. The model presents extensive interactions between the N- and C-terminal domains, the two envelope helices, and the membrane. Together, these interactions could allow for a concerted large-scale conformational change during catalysis.

In silico and in vitro approaches to elucidate the thermal stability of human UDP-glucuronosyltransferase (UGT) 1A9

UDP-Glucuronosyltransferases (UGTs) are predominant drug metabolizing enzymes in the liver and extrahepatic tissues. Human UGT1A9 is uniquely stable against heat treatment. To understand the unique properties of UGT1A9, the three-dimensional structure was constructed by homology modeling using a crystal structure of TDP-epi-vancosaminyltransferase as template. Sequence alignment analysis revealed that 13 amino acid residues (Arg42, Lys91, Ala92, Tyr106, Gly111, Tyr113, Asp115, Asn152, Leu173, Leu219, His221, Arg222, and Glu241) are unique to UGT1A9 as compared with UGT1A7, UGT1A8 and UGT1A10. To examine the roles of these residues in the conformational stability of UGT1A9, molecular dynamics simulation of the structures was carried out at 310 K and 360 K in aqueous solution for 3.0 nanoseconds. Root mean square deviation analyses revealed that Arg42, Leu173, Leu219, His221 and Arg222 were responsible for the thermal stability. Root mean square fluctuation analyses and a dynamical cross correlation map revealed that Lys91, Ala92, Tyr106, Gly111, Tyr113, Asp115, Leu219, His221, Arg222 and Glu241 were responsible for the thermal stability. In vitro study using mutants of these residues demonstrated that all these amino acids may be collectively involved in the thermal stability of UGT1A9. The results presented here provide a molecular basis for the thermal stability of human UGT1A9.

Influence of N-terminal domain histidine and proline residues on the substrate selectivities of human UDP-glucuronosyltransferase 1A1, 1A6, 1A9, 2B7, and 2B10

An N-terminal domain histidine [corresponding to position 39 of UDP-glucuronosyltransferase (UGT) 1A1] is conserved in all UGT1A and UGT2B subfamily proteins except UGT1A4 (Pro-40) and UGT2B10 (Leu-34). Unlike most UGT1A and UGT2B xenobiotic-metabolizing enzymes, UGT1A4 and UGT2B10 lack the ability to glucuronidate 4-methylumbelliferone (4MU) and 1-naphthol (1NP), both planar phenols, and naproxen (a carboxylic acid). However, only UGT1A4 glucuronidates the tertiary amines lamotrigine (LTG) and trifluoperazine (TFP). In this study, we sought to elucidate the influence of specific N-terminal histidine and proline residues on UGT enzyme substrate selectivity. The conserved N-terminal domain histidine of UGT1A1, UGT1A6, UGT1A9, and UGT2B7 was mutated to proline and leucine 34 of UGT2B10 was substituted with histidine, and the capacity of the wild-type and mutant proteins to glucuronidate 4MU, 1NP, LTG, TFP, and naproxen was characterized. Whereas UGT1A1(H39P), UGT1A6(H38P), and UGT1A9(H37P) lacked the ability to metabolize 4MU, 1NP, and naproxen, all glucuronidated LTG. K(m) values for UGT1A1(H39P) and UGT1A9(H37P) were 774 and 3812 microM, respectively, compared with 1579 microM for UGT1A4. UGT1A1(H39P) also glucuronidated TFP with a V(max)/K(m) value comparable to that of UGT1A4. In contrast to the wild-type enzyme, UGT2B10(L34H) glucuronidated 4MU and 1NP with respective K(m) values of 260 and 118 microM. UGT2B7(H35P) lacked activity toward all substrates. The data confirm a pivotal role for an N-terminal domain proline in the glucuronidation of the tertiary amines LTG and TFP by UGT1A subfamily proteins, whereas glucuronidation reactions involving proton abstraction generally, although not invariably, require a histidine at the equivalent position in both UGT1A and UGT2B enzymes.

Glucuronidation of polychlorinated biphenylols and UDP-glucuronic acid concentrations in channel catfish liver and intestine

Polychlorinated biphenylols (OH-PCBs) are potentially toxic polychlorinated biphenyl metabolites that can be eliminated by glucuronidation, catalyzed by UDP-glucuronosyltransferases (UGTs). OH-PCBs with a 3,5-dichloro-4-hydroxy substitution pattern have been detected in blood from humans and wildlife, suggesting slow elimination. In this study we assessed the glucuronidation of 4-OH-PCBs with zero, one, or two chlorine atoms flanking the 4-hydroxyl group and zero to four chlorine atoms in the aphenolic ring in microsomes from channel catfish liver and proximal intestine. Product formation was quantitated with [(14)C]UDP-glucuronic acid (UDPGA). Physiological concentrations of UDPGA were measured in preparations of liver and intestine. When the OH-PCB concentrations were varied in the presence of saturating UDPGA concentrations, glucuronidation V(max) values were higher in hepatic than in intestinal microsomes (0.40-3.4 and 0.12-0.78 nmol/min/mg of protein, respectively), whereas the K(m) values were generally lower for intestine (0.042-0.47 mM) than for liver (0.11-1.64 mM). In both tissues V(max) values with 3,5-dichloro-4-OH-PCBs were lower than with the corresponding 3-chloro-4-OH-PCBs. Varying the UDPGA concentrations in the presence of saturating concentrations of OH-PCB showed that the K(m) for UDPGA was lower in intestine (27 microM) than in liver (690 microM). The measured concentration of UDPGA in catfish liver (246-377 nmol/g) was lower than the K(m) for UDPGA, suggesting that in vivo rates of glucuronidation may be suboptimal, whereas in intestine the measured UDPGA concentration (71-258 nmol/g) was higher than the K(m) for UDPGA. Although liver has a greater glucuronidation capacity than proximal intestine, the properties of intestinal UGTs in channel catfish enable them to efficiently glucuronidate low concentrations of OH-PCBs.

Critical Roles of Residues 36 and 40 in the Phenol and Tertiary Amine Aglycone Substrate Selectivities of UDP-Glucuronosyltransferases 1A3 and 1A4

Despite high sequence identity, UGT1A3 and UGT1A4 differ in terms of substrate selectivity. UGT1A3 glucuronidates the planar phenols 1-naphthol (1-NP) and 4-methylumbelliferone (4-MU), whereas UGT1A4 converts the tertiary amines lamotrigine (LTG) and trifluoperazine (TFP) to quaternary ammonium glucuronides. Residues 45 to 154 (which incorporate 21 of the 35 amino acid differences) and 45 to 535 were exchanged between UGT1A3 and UGT1A4 to generate UGT1A3-4((45-535)), UGT1A3-4((45-154))-3, UGT1A4-3((45-535)), and UGT1A4-3((45-154))-4 hybrid proteins. Although differences in kinetic parameters were observed between the parent enzymes and chimeras, UGT1A4-3((45-535)) and UGT1A4-3((45-154))-4 [but not UGT1A3-4((45-535)) and UGT1A3-4((45-154))-3] retained the capacity to glucuronidate LTG and TFP. Likewise, UGT1A3-4((45-535)) and UGT1A3-4((45-154))-3 retained the capacity to glucuronidate 1-NP and 4-MU, but UGT1A4-3((45-535)) and UGT1A4-3((45-154))-4 exhibited low or absent activity. Within the first 44 residues, UGT1A3 and UGT1A4 differ in sequence at positions 36 and 40. "Reciprocal" mutagenesis was performed to generate the UGT1A3(I36T), UGT1A3(H40P), UGT1A4(T36I), and UGT1A4 (P40H) mutants. The T36I and P40H mutations in UGT1A4 reduced in vitro clearances for LTG and TFP glucuronidation by >90%. Conversely, the I36T and H40P mutations in UGT1A3 reduced the in vitro clearances for 1-NP and 4-MU glucuronidation by >90%. Introduction of the single H40P mutation in UGT1A3 conferred LTG and TFP glucuronidation, whereas the single T36I mutation in UGT1A4 conferred 1-NP and 4-MU glucuronidation. Thus, residues 36 and 40 of UGT1A3 and UGT1A4 are pivotal for the respective selectivities of these enzymes toward planar phenols and tertiary amines, although other regions of the proteins influence binding affinity and/or turnover.