The UDP-Glycosyltransferase (UGT) Superfamily: New Members, New Functions, and Novel Paradigms.

UDP-glycosyltransferases (UGTs) catalyze the covalent addition of sugars to a broad range of lipophilic molecules. This biotransformation plays a critical role in elimination of a broad range of exogenous chemicals and by-products of endogenous metabolism, and also controls the levels and distribution of many endogenous signaling molecules. In mammals, the superfamily comprises four families: UGT1, UGT2, UGT3, and UGT8. UGT1 and UGT2 enzymes have important roles in pharmacology and toxicology including contributing to interindividual differences in drug disposition as well as to cancer risk. These UGTs are highly expressed in organs of detoxification (e.g., liver, kidney, intestine) and can be induced by pathways that sense demand for detoxification and for modulation of endobiotic signaling molecules. The functions of the UGT3 and UGT8 family enzymes have only been characterized relatively recently; these enzymes show different UDP-sugar preferences to that of UGT1 and UGT2 enzymes, and to date, their contributions to drug metabolism appear to be relatively minor. This review summarizes and provides critical analysis of the current state of research into all four families of UGT enzymes. Key areas discussed include the roles of UGTs in drug metabolism, cancer risk, and regulation of signaling, as well as the transcriptional and posttranscriptional control of UGT expression and function. The latter part of this review provides an in-depth analysis of the known and predicted functions of UGT3 and UGT8 enzymes, focused on their likely roles in modulation of levels of endogenous signaling pathways.

Trace amine associated receptor 1: predicted effects of single nucleotide variants on structure-function in geographically diverse populations.

Trace Amine Associated Receptor 1 (TAAR1) is a novel pharmaceutical target under investigation for the treatment of several neuropsychiatric conditions. TAAR1 single nucleotide variants (SNV) have been found in patients with schizophrenia and metabolic disorders. However, the frequency of variants in geographically diverse populations and the functional effects of such variants are unknown. In this study, we aimed to characterise the distribution of TAAR1 SNVs in five different WHO regions using the Database of Genotypes and Phenotypes (dbGaP) and conducted a critical computational analysis using available TAAR1 structural data to identify SNVs affecting ligand binding and/or functional regions. Our analysis shows 19 orthosteric, 9 signalling and 16 micro-switch SNVs hypothesised to critically influence the agonist induced TAAR1 activation. These SNVs may non-proportionally influence populations from discrete regions and differentially influence the activity of TAAR1-targeting therapeutics in genetically and geographically diverse populations. Notably, our dataset presented with orthosteric SNVs D1033.32N (found only in the South-East Asian Region and Western Pacific Region) and T1945.42A (found only in South-East Asian Region), and 2 signalling SNVs (V1253.54A/T2526.36A, found in African Region and commonly, respectively), all of which have previously demonstrated to influence ligand induced functions of TAAR1. Furthermore, bioinformatics analysis using SIFT4G, MutationTaster 2, PROVEAN and MutationAssessor predicted all 16 micro-switch SNVs are damaging and may further influence the agonist activation of TAAR1, thereby possibly impacting upon clinical outcomes. Understanding the genetic basis of TAAR1 function and the impact of common mutations within clinical populations is important for the safe and effective utilisation of novel and existing pharmacotherapies.

The molecular basis of dapsone activation of CYP2C9-catalyzed nonsteroidal anti-inflammatory drug oxidation.

Positive heterotropic cooperativity, or "activation," results in an instantaneous increase in enzyme activity in the absence of an increase in protein expression. Thus, cytochrome P450 (CYP) enzyme activation presents as a potential drug-drug interaction mechanism. It has been demonstrated previously that dapsone activates the CYP2C9-catalyzed oxidation of a number of nonsteroidal anti-inflammatory drugs in vitro. Here, we conducted molecular dynamics simulations (MDS) together with enzyme kinetic investigations and site-directed mutagenesis to elucidate the molecular basis of the activation of CYP2C9-catalyzed S-flurbiprofen 4'-hydroxylation and S-naproxen O-demethylation by dapsone. Supplementation of incubations of recombinant CYP2C9 with dapsone increased the catalytic efficiency of flurbiprofen and naproxen oxidation by 2.3- and 16.5-fold, respectively. MDS demonstrated that activation arises predominantly from aromatic interactions between the substrate, dapsone, and the phenyl rings of Phe114 and Phe476 within a common binding domain of the CYP2C9 active site, rather than involvement of a distinct effector site. Mutagenesis of Phe114 and Phe476 abrogated flurbiprofen and naproxen oxidation, and MDS and kinetic studies with the CYP2C9 mutants further identified a pivotal role of Phe476 in dapsone activation. MDS additionally showed that aromatic stacking interactions between two molecules of naproxen are necessary for binding in a catalytically favorable orientation. In contrast to flurbiprofen and naproxen, dapsone did not activate the 4'-hydroxylation of diclofenac, suggesting that the CYP2C9 active site favors cooperative binding of nonsteroidal anti-inflammatory drugs with a planar or near-planar geometry. More generally, the work confirms the utility of MDS for investigating ligand binding in CYP enzymes.

Binding of SEP-363856 within TAAR1 and the 5HT1A receptor: implications for the design of novel antipsychotic drugs.

Current medications for schizophrenia typically modulate dopaminergic neurotransmission. While affecting positive symptoms, antipsychotic drugs have little clinical effect on negative symptoms and cognitive impairment. Moreover, newer 'atypical' antipsychotic drugs also have significant metabolic adverse-effects. The recent positive clinical trial of the novel drug candidate SEP-363856, which targets non-dopamine receptors (trace amine-associated receptor and the 5HT1A receptor), is a potentially promising development for the management of schizophrenia. In this perspective, we briefly overview the role of TAAR1 and the 5HT1A receptor in schizophrenia and explore the specific binding characteristics of SEP-363856 at these receptors. Molecular dynamics simulations (MDS) indicate that SEP-363856 interacts with a small, common set of conserved residues within the TAAR1 and 5HT1A ligand-binding domain. The primary interaction of SEP-363856 involves binding to the negatively charged aspartate residue (Asp1033.32, TAAR1; Asp1163.32, 5HT1A). In general, the binding of SEP-363856 within TAAR1 involves a greater number of aromatic contacts compared to 5HT1A. MDS provides important insights into the molecular basis of binding site interactions of SEP-363856 with TAAR1 and the 5HT1A receptor, which will be beneficial for understanding the pharmacological uniqueness of SEP-363856 and for the design of novel drug candidates for these newly targeted receptors in the treatment of schizophrenia and related disorders.

Chemical similarities and differences among inhibitors of nitric oxide synthase, arginase and dimethylarginine dimethylaminohydrolase-1: Implications for the design of novel enzyme inhibitors modulating the nitric oxide pathway.

Nitric oxide (NO) is a signalling molecule that controls a multitude of regulatory functions including neurotransmission, vascular tone, immune response, and angiogenesis. Regulating NO concentrations in cells using small molecules is an active area of research in the treatment of several pathologies such as cardiovascular disease, cancer, and inflammatory conditions. Small molecule-inhibition of critical NO regulatory enzymes, NO synthase (NOS), arginase, and dimethylarginine dimethyaminohydrolase-1 (DDAH1), has shown therapeutic benefits as well as limitations and is a focus of current research.In recent years, DDAH1 has been explored as a potential target to indirectly regulate NO in diseases characterized by excessive NO production. This review discusses the biological and pathophysiological role of the NO pathway, the existing inhibitors of NOS, arginase and DDAH1, and the conventional and structure-guided structure-activity relationship studies involved in their discovery. The key structural elements of amino acid-derived inhibitors responsible for selective inhibition of each enzyme, and the chemical features responsible for dual enzyme inhibition are also discussed. Finally, a synthetic scheme for developing both selective and dual inhibitors using common starting materials is provided, offering unique insights in the quest for the rational design of novel NO pathway inhibitors.

Binding of clozapine to the GABAB receptor: clinical and structural insights.

Clozapine is the gold-standard agent for treatment resistant schizophrenia but its mechanism of action remains unclear. There is emerging evidence of the potential role of the GABAB receptor in the pathogenesis of schizophrenia. It has been hypothesised that clozapine can mediate its actions via the GABAB receptor. Baclofen is currently recognised as the prototype GABAB receptor agonist. There are some potential clinical similarities between clozapine and baclofen. Indeed, baclofen has been previously proposed for use as an antipsychotic agent. Our analysis of the X-ray crystal structure of GABAB receptor along with molecular docking calculations, suggests that clozapine could directly bind to the GABAB receptor similar to that of baclofen. This finding could lead to a better understanding of the pharmacological uniqueness of clozapine, potential development of a biomarker for treatment resistant schizophrenia and the development of more targeted treatments leading to personalisation of treatment.

Structural Insights to Human Immunodeficiency Virus (HIV-1) Targets and Their Inhibition.

Human immunodeficiency virus (HIV) is a deadly virus that attacks the body's immune system, subsequently leading to AIDS (acquired immunodeficiency syndrome) and ultimately death. Currently, there is no vaccine or effective cure for this infection; however, antiretrovirals that act at various phases of the virus life cycle have been useful to control the viral load in patients. One of the major problems with antiretroviral therapies involves drug resistance. The three-dimensional structure from crystallography studies are instrumental in understanding the structural basis of drug binding to various targets. This chapter provides key insights into different targets and drugs used in the treatment from a structural perspective. Specifically, an insight into the binding characteristics of drugs at the active and allosteric sites of different targets and the importance of targeting allosteric sites for design of new-generation antiretrovirals to overcome complex and resistant forms of the virus has been reviewed.

Arginine-259 of UGT2B7 Confers UDP-Sugar Selectivity.

Enzymes of the human UDP-glycosyltransferase (UGT) superfamily typically catalyze the covalent addition of the sugar moiety from a UDP-sugar cofactor to relatively low-molecular weight lipophilic compounds. Although UDP-glucuronic acid (UDP-GlcUA) is most commonly employed as the cofactor by UGT1 and UGT2 family enzymes, UGT2B7 and several other enzymes can use both UDP-GlcUA and UDP-glucose (UDP-Glc), leading to the formation of glucuronide and glucoside conjugates. An investigation of UGT2B7-catalyzed morphine glycosidation indicated that glucuronidation is the principal route of metabolism because the binding affinity of UDP-GlcUA is higher than that of UDP-Glc. Currently, it is unclear which residues in the UGT2B7 cofactor binding domain are responsible for the preferential binding of UDP-GlcUA. Here, molecular dynamics (MD) simulations were performed together with site-directed mutagenesis and enzyme kinetic studies to identify residues within the UGT2B7 binding site responsible for the selective cofactor binding. MD simulations demonstrated that Arg259, which is located within the N-terminal domain, specifically interacts with UDP-GlcUA, whereby the side chain of Arg259 H-bonds and forms a salt bridge with the carboxylate group of glucuronic acid. Consistent with the MD simulations, substitution of Arg259 with Leu resulted in the loss of morphine, 4-methylumbelliferone, and zidovudine glucuronidation activity, but morphine glucosidation was preserved. SIGNIFICANCE STATEMENT: Despite the importance of uridine diphosphate glycosyltransferase (UGT) enzymes in drug and chemical metabolism, cofactor binding interactions are incompletely understood, as is the molecular basis for preferential glucuronidation by UGT1 and UGT2 family enzymes. The study demonstrated that long timescale molecular dynamics (MD) simulations with a UGT2B7 homology model can be used to identify critical binding interactions of a UGT protein with UDP-sugar cofactors. Further, the data provide a basis for the application of MD simulations to the elucidation of UGT-aglycone interactions.

Slow-, Tight-Binding Inhibition of CYP17A1 by Abiraterone Redefines Its Kinetic Selectivity and Dosing Regimen.

Substantial evidence underscores the clinical efficacy of inhibiting CYP17A1-mediated androgen biosynthesis by abiraterone for treatment of prostate oncology. Previous structural analysis and in vitro assays revealed inconsistencies surrounding the nature and potency of CYP17A1 inhibition by abiraterone. Here, we establish that abiraterone is a slow-, tight-binding inhibitor of CYP17A1, with initial weak binding preceding the subsequent slow isomerization to a high-affinity CYP17A1-abiraterone complex. The in vitro inhibition constant of the final high-affinity CYP17A1-abiraterone complex ( ( K i * = 0.39 nM )yielded a binding free energy of -12.8 kcal/mol that was quantitatively consistent with the in silico prediction of -14.5 kcal/mol. Prolonged suppression of dehydroepiandrosterone (DHEA) concentrations observed in VCaP cells after abiraterone washout corroborated its protracted CYP17A1 engagement. Molecular dynamics simulations illuminated potential structural determinants underlying the rapid reversible binding characterizing the two-step induced-fit model. Given the extended residence time (42 hours) of abiraterone within the CYP17A1 active site, in silico simulations demonstrated sustained target engagement even when most abiraterone has been eliminated systemically. Subsequent pharmacokinetic-pharmacodynamic (PK-PD) modeling linking time-dependent CYP17A1 occupancy to in vitro steroidogenic dynamics predicted comparable suppression of downstream DHEA-sulfate at both 1000- and 500-mg doses of abiraterone acetate. This enabled mechanistic rationalization of a clinically reported PK-PD disconnect, in which equipotent reduction of downstream plasma DHEA-sulfate levels was achieved despite a lower systemic exposure of abiraterone. Our novel findings provide the impetus for re-evaluating the current dosing paradigm of abiraterone with the aim of preserving PD efficacy while mitigating its dose-dependent adverse effects and financial burden. SIGNIFICANCE STATEMENT: With the advent of novel molecularly targeted anticancer modalities, it is becoming increasingly evident that optimal dose selection must necessarily be predicated on mechanistic characterization of the relationships between target exposure, drug-target interactions, and pharmacodynamic endpoints. Nevertheless, efficacy has always been perceived as being exclusively synonymous with affinity-based measurements of drug-target binding. This work demonstrates how elucidating the slow-, tight-binding inhibition of CYP17A1 by abiraterone via in vitro and in silico analyses was pivotal in establishing the role of kinetic selectivity in mediating time-dependent CYP17A1 engagement and eventually downstream efficacy outcomes.

Computational Prediction of the Site(s) of Metabolism and Binding Modes of Protein Kinase Inhibitors Metabolized by CYP3A4.

Protein kinase inhibitors (KIs), which are mainly biotransformed by CYP3A4-catalyzed oxidation, represent a rapidly expanding class of drugs used primarily for the treatment of cancer. Ligand- and structure-based methods were applied here to investigate whether computational approaches may be used to predict the site(s) of metabolism (SOM) of KIs and to identify amino acids within the CYP3A4 active site involved in KI binding. A data set of the experimentally determined SOMs of 31 KIs known to undergo biotransformation by CYP3A4 was collated. The structure-based (molecular docking) approach employed three CYP3A4 X-ray crystal structures to account for structural plasticity of this enzyme. Docking pose and SOM predictivity were influenced by the X-ray crystal template used for docking and the scoring function used for ranking binding poses. The best prediction of SOM (77%) was achieved using the substrate (bromoergocryptine)-bound X-ray crystal template together with the potential of mean force score. Binding interactions of KIs with CYP3A4 active site residues were generally similar to those observed for other substrates of this enzyme. The ligand-based molecular superposition approach, using bromoergocryptine from the X-ray cocrystal structure as a template, poorly predicted (42%) the SOM of KIs, although predictivity improved to 71% when the docked conformation of sorafenib was used as the template. Among the web-based approaches examined, all web servers provided excellent predictivity, with one web server predicting the SOM of 87% of the data set molecules. Computational approaches may be used to predict the SOM of KIs, and presumably other classes of CYP3A4 substrates, but predictivity varies between methods.

Warfarin resistance associated with genetic polymorphism of VKORC1: linking clinical response to molecular mechanism using computational modeling.

The variable response to warfarin treatment often has a genetic basis. A protein homology model of human vitamin K epoxide reductase, subunit 1 (VKORC1), was generated to elucidate the mechanism of warfarin resistance observed in a patient with the Val66Met mutation. The VKORC1 homology model comprises four transmembrane (TM) helical domains and a half helical lid domain. Cys132 and Cys135, located in the N-terminal end of TM-4, are linked through a disulfide bond. Two distinct binding sites for warfarin were identified. Site-1, which binds vitamin K epoxide (KO) in a catalytically favorable orientation, shows higher affinity for S-warfarin compared with R-warfarin. Site-2, positioned in the domain occupied by the hydrophobic tail of KO, binds both warfarin enantiomers with similar affinity. Displacement of Arg37 occurs in the Val66Met mutant, blocking access of warfarin (but not KO) to Site-1, consistent with clinical observation of warfarin resistance.

Cytochrome P450 structure-function: insights from molecular dynamics simulations.

Cytochrome P450 (CYP) family 1, 2, and 3 enzymes play an essential role in the metabolic clearance and detoxification of a myriad of structurally and chemically diverse drugs and non-drug xenobiotics. The individual CYP enzymes exhibit distinct substrate and inhibitor selectivities, and hence differing patterns of inhibitory drug-drug interactions. In addition, CYP enzymes differ in terms of regulation of expression, genetic polymorphism, and environmental factors that alter activity. The availability of three-dimensional structures from X-ray crystallography have been invaluable for understanding the structural basis of the ligand selectivity of CYP enzymes. Moreover, the X-ray crystal structures demonstrate that CYP proteins exhibit marked flexibility, particularly around the active site, and the principle of ligand-induced conformational changes is now well accepted. Recent studies have demonstrated that molecular dynamics simulations (MDS) provide an additional approach for modeling the structural flexibility of CYP enzymes, both in the presence and absence of bound ligand, and understanding the functional consequences of plasticity. However, most of the MDS studies reported to date have utilized short simulation time scales, and few have validated the computationally-generated data experimentally (e.g. by site-directed mutagenesis and enzyme kinetic approaches). Although modeling approaches require further development and validation, MDS has the potential to provide a deeper understanding of CYP structure-function than is available from experimental techniques such as X-ray crystallography alone.

Human UDP-Glucuronosyltransferase (UGT) 2B10: Validation of Cotinine as a Selective Probe Substrate, Inhibition by UGT Enzyme-Selective Inhibitors and Antidepressant and Antipsychotic Drugs, and Structural Determinants of Enzyme Inhibition.

Although there is evidence for an important role of UGT2B10 in the N-glucuronidation of drugs and other xenobiotics, the inhibitor selectivity of this enzyme is poorly understood. This study sought primarily to characterize the inhibition selectivity of UGT2B10 by UDP-glucuronosyltransferase (UGT) enzyme-selective inhibitors used for reaction phenotyping, and 34 antidepressant and antipsychotic drugs that contain an amine functional group. Initial studies demonstrated that cotinine is a highly selective substrate of human liver microsomal UGT2B10. The kinetics of cotinine N-glucuronidation by recombinant UGT and human liver microsomes (± bovine serum albumin) were consistent with the involvement of a single enzyme. Of the UGT enzyme-selective inhibitors employed for reaction phenotyping, only the UGT2B4/7 inhibitor fluconazole reduced recombinant UGT2B10 activity to an appreciable extent. The majority of antidepressant and antipsychotic drugs screened for effects on UGT2B10 inhibited enzyme activity with IC50 values <100 µM. The most potent inhibition was observed with the tricyclic antidepressants amitriptyline and doxepin and the tetracyclic antidepressant mianserin, and the structurally related compounds desloratadine and loratadine. Molecular modeling using a ligand-based approach indicated that hydrophobic and charge interactions are involved in inhibitor binding, whereas spatial features influence the potency of UGT2B10 inhibition. Respective mean Ki,u (± S.D.) values for amitriptyline, doxepin, and mianserin inhibition of human liver microsomal UGT2B10 were 0.61 ± 0.05, 0.95 ± 0.18, and 0.43 ± 0.01 µM. In vitro-in vivo extrapolation indicates that these drugs may perpetrate inhibitory drug-drug interactions when coadministered with compounds that are cleared predominantly by UGT2B10.

A Fragment-Based Approach for the Computational Prediction of the Nonspecific Binding of Drugs to Hepatic Microsomes.

Correction for the nonspecific binding (NSB) of drugs to liver microsomes is essential for the accurate measurement of the kinetic parameters Km and Ki, and hence in vitro-in vivo extrapolation to predict hepatic clearance and drug-drug interaction potential. Although a number of computational approaches for the estimation of drug microsomal NSB have been published, they generally rely on compound lipophilicity and charge state at the expense of other physicochemical and chemical properties. In this work, we report the development of a fragment-based hologram quantitative structure activity relationship (HQSAR) approach for the prediction of NSB using a database of 132 compounds. The model has excellent predictivity, with a noncross-validated r2 of 0.966 and cross-validated r2 of 0.680, with a predictive r2 of 0.748 for an external test set comprising 34 drugs. The HQSAR method reliably predicted the fraction unbound in incubations of 95% of the training and test set drugs, excluding compounds with a steroid or morphinan 4,5-epoxide nucleus. Using the same data set of compounds, performance of the HQSAR method was superior to a model based on logP/D as the sole descriptor (predictive r2 for the test set compounds, 0.534). Thus, the HQSAR method provides an alternative approach to laboratory-based procedures for the prediction of the NSB of drugs to liver microsomes, irrespective of the drug charge state (acid, base, or neutral).

A novel function for UDP glycosyltransferase 8: galactosidation of bile acids.

The human UDP glycosyltransferase (UGT) superfamily comprises four families of enzymes that catalyze the addition of sugar residues to small lipophilic chemicals. The UGT1 and UGT2 enzymes use UDP-glucuronic acid, and UGT3 enzymes use UDP-N-acetylglucosamine, UDP-glucose, and UDP-xylose to conjugate xenobiotics, including drugs and endobiotics such as metabolic byproducts, hormones, and signaling molecules. This metabolism renders the substrate more polar and more readily excreted from the body and/or functionally inactive. The fourth UGT family, called UGT8, contains only one member that, unlike other UGTs, is considered biosynthetic. UGT8 uses UDP galactose to galactosidate ceramide, a key step in the synthesis of brain sphingolipids. To date other substrates for this UGT have not been identified and there has been no suggestion that UGT8 is involved in metabolism of endo- or xenobiotics. We re-examined the functions of UGT8 and discovered that it efficiently galactosidates bile acids and drug-like bile acid analogs. UGT8 conjugates bile acids ∼60-fold more efficiently than ceramide based on in vitro assays with substrate preference deoxycholic acid > chenodeoxycholic acid > cholic acid > hyodeoxycholic acid > ursodeoxycholic acid. Activities of human and mouse UGT8 are qualitatively similar. UGT8 is expressed at significant levels in kidney and gastrointestinal tract (intestine, colon) where conjugation of bile acids is likely to be metabolically significant. We also investigate the structural determinants of UDP-galactose selectivity. Our novel findings suggest a new role for UGT8 as a modulator of bile acid homeostasis and signaling.

Insights into the UDP-sugar selectivities of human UDP-glycosyltransferases (UGT): a molecular modeling perspective.

Enzymes of the human uridine diphosphate (UDP)-glycosyltransferase (UGT) superfamily typically catalyze the covalent addition of a sugar from UDP-sugar cofactors to relatively small lipophilic compounds. The sugar conjugates are often biologically less active with improved water-solubility, facilitating more effective elimination from the body. Experimental data indicate that UGT proteins exhibit differing selectivities toward various UDP-sugars. Although, three-dimensional (3D) structures of UGT proteins are required to provide insights into the UDP-sugar selectivities observed for the various UGT proteins, there are currently, no experimental structures available for human UGTs bound to UDP-sugar(s). Thus, the absence of 3D structures poses a major challenge for analyzing UDP-sugar selectivity at an atomic level. In this commentary, we highlight the application of comparative homology modeling for understanding the UDP-sugar selectivities of UGT proteins. Homology models of the C-terminal (CT) domain indicate a highly conserved structural fold across the UGT family with backbone root mean-squared deviations (rmsds) between 0.066 and 0.079 Å with respect to the UGT2B7-CT X-ray crystal structure. The models show that four residues in the terminal portion of the CT signature sequence play an important role in UDP-sugar selectivity. The N-terminal domain is less likely to be associated with UDP-sugar selectivity, although, a conserved residue, Arg-259 (UGT2B7 numbering) in the UGT 1 and 2 families may influence UDP-sugar selectivity. Overall, the models demonstrate excellent agreement with experimental observations in predicting the key residues that influence the selectivity of UDP-sugar binding.