The P450 superfamily: update on new sequences, gene mapping, and recommended nomenclature

We provide here a list of 154 P450 genes and seven putative pseudogenes that have been characterized as of October 20, 1990. These genes have been described in a total of 23 eukaryotes (including nine mammalian and one plant species) and six prokaryotes. Of 27 gene families so far described, 10 exist in all mammals. These 10 families comprise 18 subfamilies, of which 16 and 14 have been mapped in the human and mouse genomes, respectively; to date, each subfamily appears to represent a cluster of tightly linked genes. We propose here a modest revision of the initially proposed (Nebert et al., DNA 6, 1-11, 1987) and updated (Nebert et al., DNA 8, 1-13, 1989) nomenclature system based on evolution of the superfamily. For the gene we recommend that the italicized root symbol CYP for human (Cyp for mouse), representing cytochrome P450, be followed by an Arabic number denoting the family, a letter designating the subfamily (when two or more exist), and an Arabic numeral representing the individual gene within the subfamily. A hyphen should precede the final number in mouse genes. We suggest that the human nomenclature system be used for other species. This system is consistent with our earlier proposed nomenclature for P450 of all eukaryotes and prokaryotes, except that we are discouraging the future use of cumbersome Roman numerals.     

Disruption of the ugt1 locus in mice resembles human Crigler-Najjar type I disease

The 9 UDP-glucuronosyltranferases (UGTs) encoded by the UGT1 locus in humans are key enzymes in the metabolism of most drugs as well as endogenous substances such as bile acids, fatty acids, steroids, hormones, neurotransmitters, and bilirubin. Severe unconjugated hyperbilirubinemia in humans that suffer from Crigler-Najjar type I disease results from lesions in the UGT1A1 gene and is often fatal. To examine the physiological importance of the Ugt1 locus in mice, this locus was rendered non-functional by interrupting exon 4 to create Ugt1(-/-) mice. Because UGT1A1 in humans is responsible for 100% of the conjugated bilirubin, it followed that newborn Ugt1(-/-) mice developed serum levels of unconjugated bilirubin that were 40-60 times higher than Ugt1(+/-) or wild-type mice. The result of extreme unconjugated bilirubin in Ugt1(-/-) mice, comparable to the induced levels noted in patients with Crigler-Najjar type 1 disease, is fatal in neonatal Ugt1(-/-) mice within 2 weeks following birth. The extreme jaundice is present as a phenotype in skin color after 8 h. Neonatal Ugt1(-/-) mice exhibit no detectable UGT1A-specific RNA, which corresponds to a complete absence of UGT1A proteins in liver microsomes. Conserved glucuronidation activity attributed to the Ugt1 locus can be defined in Ugt1(-/-) mice, because UGT2-dependent glucuronidation activity is unaffected. Remarkably, the loss of UGT1A functionality in liver results in significant alterations in cellular metabolism as investigated through changes in gene expression. Thus, the loss of UGT1A function in Ugt1(-/-) mice leads to a metabolic syndrome that can serve as a model to further investigate the toxicities associated with unconjugated bilirubin and the impact of this disease in humans.     

Steroid UDP glucuronosyltransferases

The glucuronidation of steroids is a major process necessary for their elimination in the bile and urine. In general, steroid glucuronides are biologically less reactive than their parent steroids. However, in some cases often associated with disease and steroid therapy, more reactive or toxic glucuronides may be formed. The concentrations of specific steroid glucuronides in the blood may also indicate hormonal imbalances and may funnction as diagnostic markers of genetic defects in steroid synthesis and metabolism. In this review, the forms of UDP glucuronosyltransferase involved in steroid glucuronidation are described in terms of their specificities, functional domains and regulation. The available evidence suggests that steroid glucuronidation is mainly carried out by members of the UGT2B subfamily which are encoded by genes containing 6 exons. Members of this subfamily exhibit a regioselectively in their glucuronidation of steroids that is mediated by domains in the amino-terminal half on the protein encoded by exons 1 and 2. Although much of this review will describe studies in the rat, preliminary evidence indicates that a similar situation may exist in humans.     

Murine aldo-keto reductase family 1 subfamily B: identification of AKR1B8 as an ortholog of human AKR1B10

Aldo-keto reductase family 1 member B10 (AKR1B10), over-expressed in multiple human cancers, might be implicated in cancer development and progression via detoxifying cytotoxic carbonyls and regulating fatty acid synthesis. In the present study, we investigated the ortholog of AKR1B10 in mice, an ideal modeling organism greatly contributing to human disease investigations. In the mouse, there are three aldo-keto reductase family 1 subfamily B (AKR1B) members, i.e., AKR1B3, AKR1B7, and AKR1B8. Among them, AKR1B8 has the highest similarity to human AKR1B10 in terms of amino acid sequence, computer-modeled structures, substrate spectra and specificity, and tissue distribution. More importantly, similar to human AKR1B10, mouse AKR1B8 associates with murine acetyl-CoA carboxylase-α and mediates fatty acid synthesis in colon cancer cells. Taken together, our data suggest that murine AKR1B8 is the ortholog of human AKR1B10.     

Adaptive evolution of multiple-variable exons and structural diversity of drug-metabolizing enzymes

Background  

The human genome contains a large number of gene clusters with multiple-variable-first exons, including the drug-metabolizing UDP glucuronosyltransferase (UGT1) and I-branching β-1,6-N-acetylglucosaminyltransferase (GCNT2, also known as IGNT) clusters, organized in a tandem array, similar to that of the protocadherin (PCDH), immunoglobulin (IG), and T-cell receptor (TCR) clusters. To gain insight into the evolutionary processes that may have shaped their diversity, we performed comprehensive comparative analyses for vertebrate multiple-variable-first-exon clusters.     

A novel mouse kinesin of the UNC-104/KIF1 subfamily encoded by the Kif1b gene

Kinesin and kinesin-related proteins are microtubule-dependent motor proteins that transport organelles. We have cloned and sequenced a full-length 9924 bp mouse cDNA for a new kinesin of the UNC-104/KIF1 subfamily. Northern blot analysis of mouse RNAs detected high levels of a 10 kb mRNA in brain and eye, but lower levels in other tissues. Human RNA dot-blot analysis detected this mRNA in all tissues examined, although at different levels. The overall structure of the new kinesin (predicted size 204 kDa) was most similar to mouse KIF1A; however, 2.1 kb of the 5' portion of the cDNA were identical to the published sequence for KIF1B (Nangaku, M., Sato-Yoshitake, R., Okada, Y., Noda, Y., Takemura, R., Yamazaki, H., Hirokawa, N., 1994. KIF1B, a novel microtubule plus end-directed monomeric motor protein for transport of mitochondria. Cell 79, 1209-1220). We localized the Kif1b gene to the distal end of mouse Chromosome 4 by haplotype analysis of an interspecific backcross from The Jackson Laboratory. We had previously mapped the gene for the novel kinesin to the same location (Gong, T.-W.L., Burmeister, M., Lomax, M.I., 1996b. The novel gene D4Mille maps to mouse Chromosome 4 and human Chromosome 1p36. Mamm. Genome 7, 790-791). We conclude, therefore, that the Kif1b gene generates two major kinesin isoforms by alternative splicing. The shorter 7.8 kb mRNA encodes a 130 kDa kinesin, KIF1Bp130, whereas the 10 kb mRNA encodes a 204 kDa kinesin, KIF1Bp204. In addition, alternative splicing of two exons in the conserved region adjacent to the motor domain generates four different isoforms of each kinesin, leading to eight kinesin isoforms derived from the Kif1b gene.     

An updated nomenclature for keratin-associated proteins (KAPs)

Most protein in hair and wool is of two broad types: keratin intermediate filament-forming proteins (commonly known as keratins) and keratin-associated proteins (KAPs). Keratin nomenclature was reviewed in 2006, but the KAP nomenclature has not been revised since 1993. Recently there has been an increase in the number of KAP genes (KRTAPs) identified in humans and other species, and increasingly reports of variation in these genes. We therefore propose that an updated naming system is needed to accommodate the complexity of the KAPs. It is proposed that the system is founded in the previous nomenclature, but with the abbreviation sp-KAPm-nL*x for KAP proteins and sp-KRTAPm-n(p/L)*x for KAP genes. In this system "sp" is a unique letter-based code for different species as described by the protein knowledge-based UniProt. "m" is a number identifying the gene or protein family, "n" is a constituent member of that family, "p" signifies a pseudogene if present, "L" if present signifies "like" and refers to a temporary "place-holder" until the family is confirmed and "x" signifies a genetic variant or allele. We support the use of non-italicised text for the proteins and italicised text for the genes. This nomenclature is not that different to the existing system, but it includes species information and also describes genetic variation if identified, and hence is more informative. For example, GenBank sequence JN091630 would historically have been named KRTAP7-1 for the gene and KAP7-1 for the protein, but with the proposed nomenclature would be SHEEP-KRTAP7-1*A and SHEEP-KAP7-1*A for the gene and protein respectively. This nomenclature will facilitate more efficient storage and retrieval of data and define a common language for the KAP proteins and genes from all mammalian species.     

Trace amine-associated receptors form structurally and functionally distinct subfamilies of novel G protein-coupled receptors

Trace amines are endogenous compounds structurally related to classical biogenic amines that have been studied for decades, triggered by their link to psychiatric conditions of high epidemiological and economical relevance. The understanding of their pharmacology on the molecular level was hampered until the recent discovery of trace-amine-specific receptors. We completed the identification of all members of this novel GPCR family in human, chimpanzee, rat, and mouse and observed remarkable interspecies differences, even between human and chimpanzee. The analysis of the chromosomal localizations, phylogenetic relationships, and ligand pocket vectors reveals three distinct receptor subfamilies. As most of these receptors do not respond to trace amines, each subfamily will presumably have a distinct pharmacological profile, which remains to be identified. We propose a uniform nomenclature describing this novel GPCR family in all mammalian species as trace-amine-associated receptors (TAARs), which resolves the ambiguities and contradictions of the previous naming.     

Cloning and characterization of a canine UDP-glucuronosyltransferase.

UDP-glucuronosyltransferases (UGTs) are a major family of enzymes catalyzing the transfer of glucuronic acid to a range of endogenous compounds and xenobiotics facilitating their elimination in either urine or bile. Although the dog is commonly used in drug metabolism studies, relatively little is known about the expression and activity of UGTs in this species. This report describes the molecular cloning and functional characterization of the first dog UGT, UGT1A6. The cloned protein is composed of 528 amino acids with the variable region demonstrating a 67-72% identity with the variable regions of mouse, rat, and human UGT1A6. The enzyme expressed stably in V79 cells predominantly catalyzed the glucuronidation of simple, planar phenols (e.g., for 1-naphthol, K(m) = 41 microM, V(max) = 0.07 nmol/min/mg protein), a class of compounds extensively glucuronidated by human UGT1A6. Based on sequence homology and common catalytic activity, this dog UGT1A protein appears to be the canine orthologue of human UGT1A6.     

Cloning and functional expression of UGT genes encoding sterol glucosyltransferases from Saccharomyces cerevisiae, Candida albicans, Pichia pastoris, and Dictyostelium discoideum

Sterol glucosides, typical membrane-bound lipids of many eukaryotes, are biosynthesized by a UDP-glucose:sterol glucosyltransferase (EC 2. 4.1.173). We cloned genes from three different yeasts and from Dictyostelium discoideum, the deduced amino acid sequences of which all showed similarities with plant sterol glucosyltransferases (Ugt80A1, Ugt80A2). These genes from Saccharomyces cerevisiae (UGT51 = YLR189C), Pichia pastoris (UGT51B1), Candida albicans (UGT51C1), and Dictyostelium discoideum (ugt52) were expressed in Escherichia coli. In vitro enzyme assays with cell-free extracts of the transgenic E. coli strains showed that the genes encode UDP-glucose:sterol glucosyltransferases which can use different sterols such as cholesterol, sitosterol, and ergosterol as sugar acceptors. An S. cerevisiae null mutant of UGT51 had lost its ability to synthesize sterol glucoside but exhibited normal growth under various culture conditions. Expression of either UGT51 or UGT51B1 in this null mutant under the control of a galactose-induced promoter restored sterol glucoside synthesis in vitro. Lipid extracts of these cells contained a novel glycolipid. This lipid was purified and identified as ergosterol-beta-D-glucopyranoside by nuclear magnetic resonance spectroscopy. These data prove that the cloned genes encode sterol-beta-D-glucosyltransferases and that sterol glucoside synthesis is an inherent feature of eukaryotic microorganisms.     

Functional specialization of UDP‐glycosyltransferase 73P12 in licorice to produce a sweet triterpenoid saponin,glycyrrhizin

Glycyrrhizin, a sweet triterpenoid saponin found in the roots and stolons of Glycyrrhiza species (licorice), is an important active ingredient in traditional herbal medicine. We previously identified two cytochrome P450 monooxygenases, CYP88D6 and CYP72A154, that produce an aglycone of glycyrrhizin, glycyrrhetinic acid, in Glycyrrhiza uralensis. The sugar moiety of glycyrrhizin, which is composed of two glucuronic acids, makes it sweet and reduces its side‐effects. Here, we report that UDP‐glycosyltransferase (UGT) 73P12 catalyzes the second glucuronosylation as the final step of glycyrrhizin biosynthesis in G. uralensis; the UGT73P12 produced glycyrrhizin by transferring a glucuronosyl moiety of UDP‐glucuronic acid to glycyrrhetinic acid 3‐O‐monoglucuronide. We also obtained a natural variant of UGT73P12 from a glycyrrhizin‐deficient (83‐555) strain of G. uralensis. The natural variant showed loss of specificity for UDP‐glucuronic acid and resulted in the production of an alternative saponin, glucoglycyrrhizin. These results are consistent with the chemical phenotype of the 83‐555 strain, and suggest the contribution of UGT73P12 to glycyrrhizin biosynthesis in planta. Furthermore, we identified Arg32 as the essential residue of UGT73P12 that provides high specificity for UDP‐glucuronic acid. These results strongly suggest the existence of an electrostatic interaction between the positively charged Arg32 and the negatively charged carboxy group of UDP‐glucuronic acid. The functional arginine residue and resultant specificity for UDP‐glucuronic acid are unique to UGT73P12 in the UGT73P subfamily. Our findings demonstrate the functional specialization of UGT73P12 for glycyrrhizin biosynthesis during divergent evolution, and provide mechanistic insights into UDP‐sugar selectivity for the rational engineering of sweet triterpenoid saponins.     

Enantioselective Inhibition of Carprofen Towards UDP‐glucuronosyltransferase (UGT) 2B7

UDP‐glucuronosyltransferases (UGTs)‐catalyzed glucuronidation conjugation reaction plays an important role in the elimination of many important clinical drugs and endogenous substances. The present study aims to investigate the enantioselective inhibition of carprofen towards UGT isoforms. In vitro a recombinant UGT isoforms‐catalyzed 4‐methylumbelliferone (4‐MU) glucuronidation incubation mixture was used to screen the inhibition potential of (R)‐carprofen and (S)‐carprofen towards multiple UGT isoforms. The results showed that (S)‐carprofen exhibited stronger inhibition potential than (R)‐carprofen towards UGT2B7. However, no significant difference was observed for the inhibition of (R)‐carprofen and (S)‐carprofen towards other UGT isoforms. Furthermore, the inhibition kinetic behavior was compared for the inhibition of (S)‐carprofen and (R)‐carprofen towards UGT2B7. A Lineweaver–Burk plot showed that both (S)‐carprofen and (R)‐carprofen exhibited competitive inhibition towards UGT2B7‐catalyzed 4‐MU glucuronidation. The inhibition kinetic parameter (K_i) was calculated to be 7.0 μM and 31.1 μM for (S)‐carprofen and (R)‐carprofen, respectively. Based on the standard for drug–drug interaction, the threshold for (S)‐carprofen and (R)‐carprofen to induce a drug–drug interaction is 0.7 μM and 3.1 μM, respectively. In conclusion, enantioselective inhibition of carprofen towards UDP‐glucuronosyltransferase (UGT) 2B7 was demonstrated in the present study. Using the in vitro inhibition kinetic parameter, the concentration threshold of (S)‐carprofen and (R)‐carprofen to possibly induce the drug–drug interaction was obtained. Therefore, clinical monitoring of the plasma concentration of (S)‐carprofen is more important than (R)‐carprofen to avoid a possible drug–drug interaction between carprofen and the drugs mainly undergoing UGT2B7‐catalyzed metabolism. Chirality 27:189–193, 2015. © 2014 Wiley Periodicals, Inc.     

Inter and intra-ethnic differences in the distribution of the molecular variants of TPMT, UGT1A1 and MDR1 genes in the South Indian population

Molecular variants of polymorphic drug metabolizing enzymes and drug transporters are attributed to differences in individual's therapeutic response and drug toxicity in different populations. We sought to determine the genotype and allele frequencies of polymorphisms for major phase II drug-metabolizing enzymes (TPMT, UGT1A1) and drug transporter (MDR1) in South Indians. Allelic variants of TPMT (*2,*3A,*3B,*3C & *8), UGT1A1 (TA)6>7 and MDR1 (2677G>T/A & 3435C>T) were evaluated in 450-608 healthy South Indian subjects. Genomic DNA was extracted by phenol-chloroform method and genotype was determined by PCR-RFLP, qRT-PCR, allele specific PCR, direct sequencing and SNaPshot techniques. The frequency distributions of TPMT, UGT1A1 and MDR1 gene polymorphisms were compared between the individual 4 South Indian populations viz., Tamilian, Kannadiga, Andhrite and Keralite. The combined frequency distribution of the South Indian populations together, was also compared with that of other major populations. The allele frequencies of TPMT*3C, UGT1A1 (TA)7, MDR1 2677T, 2677A and 3435T were 1.2, 39.8, 60.3, 3.7, and 61.6% respectively. The other variant alleles such as TPMT*2, *3A, *3B and *8 were not identified in the South Indian population. Sub-population analysis showed that the distribution of UGT1A1 (TA)6>7 and MDR1 allelic variants differed between the four ethnic groups. However, the frequencies of TPMT*3C allele were similar in the four South Indian populations. The distribution of TPMT, UGT1A1 and MDR1 gene polymorphisms of the South Indian population was significantly different from other populations.     

Genomic structure of murine macrophage inflammatory protein-1 alpha and conservation of potential regulatory sequences with a human homolog, LD78

The gene for a murine macrophage inflammatory cytokine, MIP-1 alpha, belongs to a newly recognized superfamily encoding small, inducible peptides shown to be up-regulated in association with cellular activation or transformation (tentatively designated the scy, or small cytokine, gene family). Secreted scy family peptides as a group, and MIP-1 alpha in particular, have inflammatory and mitogenic activities, and the family has been divided into CXC and CC subfamilies according to the spacing of conserved cysteine residues in the primary amino acid sequences. We have isolated and characterized a genomic clone encoding the CC subfamily member MIP-1 alpha. The organization of the murine MIP-1 alpha gene into three exons interrupted by two introns is identical to that found for other members of the CC subfamily (e.g., huLD78, muJE, huJE/MCP-1, muTCA3, and hul-309), which has been taken as evidence of evolution from a common ancestral gene. With the exception of the ratPF4 gene, which shares the two-intron/three-exon pattern typical of the CC subfamily, sequenced genes encoding CXC subfamily peptides (e.g., hulL-8 and hulP-10) include an additional intervening sequence that creates a fourth exon. Genomic nucleotide sequences 5' of the MIP-1 alpha cap site are highly homologous to corresponding regions of the human gene encoding a CC peptide variously designated as LD78/GOS19/pAT464, including consensus regulatory motifs in common, reinforcing the contention that MIP-1 alpha and LD78 may be interspecies homologs.     

A new and unified nomenclature for male fertility restorer (RF) proteins in higher plants

The male fertility restorer (RF) proteins belong to extended protein families associated with the cytoplasmic male sterility in higher plants. Up till now, there is no devised nomenclature for naming the RF proteins. The systematic sequencing of new plant species in recent years has uncovered the existence of several novel RF genes and their encoded proteins. Their naming has been simply arbitrary and could not be adequately handled in the context of comparative functional genomics. We propose in this study a unified nomenclature for the RF extended protein families across all plant species. This new and unified nomenclature relies upon previously developed nomenclature for the first ever characterized RF gene, RF2A/ALDH2B2, a member of ALDH gene superfamily, and adheres to the guidelines issued by the ALDH Genome Nomenclature Committees. The proposed nomenclature reveals that RF gene superfamily encodes currently members of 51 families. This unified nomenclature accommodates functional RF genes and pseudogenes, and offers the flexibility needed to incorporate additional RFs as they become available in future. In addition, we provide a phylogenetic relationship between the RF extended families and use computational protein modeling to demonstrate the high divergence of RF functional specializations through specific structural features of selected members of RF superfamily.     

Glucuronidation of the steroid enantiomers ent-17β-estradiol, ent-androsterone and ent-etiocholanolone by the human UDP-glucuronosyltransferases

Steroids enantiomers are interesting compounds for detailed exploration of drug metabolizing enzymes, such as the UDP-glucuronosyltransferases (UGTs). We have now studied the glucuronidation of the enantiomers of estradiol, androsterone and etiocholanolone by the 19 human UGTs of subfamilies 1A, 2A and 2B. The results reveal that the pattern of human UGTs of subfamily 2B that glucuronidate ent-17β-estradiol, particularly 2B15 and 2B17, resembles the glucuronidation of epiestradiol (17α-estradiol) rather than 17β-estradiol, the main physiological estrogen. The UGTs of subfamilies 1A and 2A exhibit higher degree of regioselectivity than enantioselectivity in the conjugation of these estradiols, regardless of whether the activity is primarily toward the non-chiral site, 3-OH (UGT1A1, UGT1A3, UGT1A7, UGT1A8 and, above all, UGT1A10), or the 17-OH (UGT1A4). In the cases of etiocholanolone and androsterone, glucuronidation of the ent-androgens, like the conjugation of the natural androgens, is mainly catalyzed by UGTs of subfamilies 2A and 2B. Nevertheless, the glucuronidation of ent-etiocholanolone and ent-androsterone by both UGT2B7 and UGT2B17 differs considerably from their respective activity toward the corresponding endogenous androgens, whereas UGT2A1-catalyzed conjugation is much less affected by the stereochemistry differences. Kinetic analyses reveal that the K(m) value of UGT2A1 for ent-estradiol is much higher than the corresponding value in the other two high activity enzymes, UGT1A10 and UGT2B7. Taken together, the results highlight large enantioselectivity differences between individual UGTs, particularly those of subfamily 2B.