Hidden Markov model analysis of motifs in steroid dehydrogenases and their homologs

Analysis of a large cluster of SLC22 transporter genes, including novel USTs, reveals species-specific amplification of subsets of family members

When the organic anion transporter Oat1 was first identified as NKT (Lopez-Nieto CE, You G, Bush KT, Barros EJ, Beier DR, Nigam SK. J Biol Chem 272: 6471-6478, 1997), it was argued that it, together with Oct1, may be part of a larger subfamily (now known as SLC22) involved in organic ion and xenobiotic transport. The least studied among SLC22 transporters are the so-called unknown substrate transporters (USTs). Here, five novel genes located in a cluster on mouse chromosome 19, immediately between Slc22a8 (Oat3)/Slc22a6 (Oat1) and Slc22a19 (Oat5), were identified as homologs of human USTs. These genes display preferential expression in liver and kidney, and one gene, AB056422, has several splicing variants with differential tissue expression and embryonic expression. Along with Slc22a6, Slc22a8, and Slc22a19, these Usts define the largest known cluster of mammalian Slc22 genes. Given the established functions of Oats, these genes may also be involved in organic anion transport. Usts have characteristic motifs and share a signature residue in the possible active site of transmembrane domain 7, a conserved, positively charged, amino acid, Arg356, possibly a site for interaction with organic anions. In certain species, Oat1 and Oat3 appeared to be highly conserved, whereas the Ust part of this cluster appeared to undergo repeated species-specific amplification, suggesting strong environmental selection pressure, and perhaps providing an explanation for copy number variation in the human locus. One Ust amplification in mouse appears to be recent. This cluster may be coordinately regulated and under selective pressure in a species-specific manner.

Genetic manipulation of 11β-hydroxysteroid dehydrogenases in mice

11β-Hydroxysteroid dehydrogenases (HSDs) interconvert active 11-hydroxy glucocorticoids (cortisol, corticosterone) and their inert 11-keto derivatives (cortisone, 11-dehydrocorticosterone). 11β-HSD type 1 is a predominant reductase that regenerates active glucocorticoids in expressing cells, thus amplifying local glucocorticoid action, whereas 11β-HSD type 2 catalyzes rapid dehydrogenation, potently inactivating intracellular glucocorticoids. Both isozymes thus regulate receptor activation by substrate availability. Spatial and temporal regulation of expression are important determinants of the physiological roles of 11β-HSDs, with each isozyme exhibiting a distinct, tissue-restricted pattern together with dynamic regulation during development and in response to environmental challenges, including diet and stress. Transgenic approaches in the mouse have contributed significantly toward an understanding of the importance of these prereceptor regulatory mechanisms on corticosteroid receptor activity and have highlighted its potential relevance to human health and disease. Here we discuss current ideas of the physiological roles of 11β-HSDs, with emphasis on the key contributions made by studies of 11β-HSD gene manipulation in vivo.

11beta-hydroxysteroid dehydrogenase type 1: a tissue-specific regulator of glucocorticoid response

11beta-hydroxysteroid dehydrogenase type 1 (11beta-HSD1) interconverts inactive cortisone and active cortisol. Although bidirectional, in vivo it is believed to function as a reductase generating active glucocorticoid at a prereceptor level, enhancing glucocorticoid receptor activation. In this review, we discuss both the genetic and enzymatic characterization of 11beta-HSD1, as well as describing its role in physiology and pathology in a tissue-specific manner. The molecular basis of cortisone reductase deficiency, the putative "11beta-HSD1 knockout state" in humans, has been defined and is caused by intronic mutations in HSD11B1 that decrease gene transcription together with mutations in hexose-6-phosphate dehydrogenase, an endoluminal enzyme that provides reduced nicotinamide-adenine dinucleotide phosphate as cofactor to 11beta-HSD1 to permit reductase activity. We speculate that hexose-6-phosphate dehydrogenase activity and therefore reduced nicotinamide-adenine dinucleotide phosphate supply may be crucial in determining the directionality of 11beta-HSD1 activity. Therapeutic inhibition of 11beta-HSD1 reductase activity in patients with obesity and the metabolic syndrome, as well as in glaucoma and osteoporosis, remains an exciting prospect.

The Drosophila Bruno paralogue Bru-3 specifically binds the EDEN translational repression element

We reported in our previous work that the EDEN-dependent translational repression of maternal mRNAs was conserved between Drosophila and Xenopus. In Xenopus, this repression is achieved through the binding of EDEN to the Bruno-like factor, EDEN-BP. We show in the present work that the Drosophila Bruno paralogue, the 45 kDa Bru-3 protein (p45), binds specifically to the EDEN element and acts as a homodimer. We describe for the first time a previously undetected 67 amino acid domain, found in the divergent linker region, the lsm domain (lsm stands for linker-specific motif). We propose that the presence of this domain in a subset of the Bruno-like proteins, including Bru-3, EDEN-BP and CUG-BP but not Bruno nor its other paralogue Bru-2, might be responsible for specific RNA recognition. Interestingly, comparative structural analyses using threaders and molecular modelling suggest that the new domain might be distantly related to the first RNA recognition motif of the Drosophila sex-lethal protein (sxl). The phylogenetic analyses and the experimental data based on its specific binding to the EDEN element support the conclusion that Bru-3 is an EDEN-BP/CUG-BP orthologue.

Mining protein sequences for motifs

We use methods from Data Mining and Knowledge Discovery to design an algorithm for detecting motifs in protein sequences. The algorithm assumes that a motif is constituted by the presence of a "good" combination of residues in appropriate locations of the motif. The algorithm attempts to compile such good combinations into a "pattern dictionary" by processing an aligned training set of protein sequences. The dictionary is subsequently used to detect motifs in new protein sequences. Statistical significance of the detection results are ensured by statistically determining the various parameters of the algorithm. Based on this approach, we have implemented a program called GYM. The Helix-Turn-Helix motif was used as a model system on which to test our program. The program was also extended to detect Homeodomain motifs. The detection results for the two motifs compare favorably with existing programs. In addition, the GYM program provides a lot of useful information about a given protein sequence.

Evolution of adrenal and sex steroid action in vertebrates: a ligand-based mechanism for complexity

Various explanations have been proposed to account for complex differentiation and development in humans, despite the human genome containing only two to three times the number of genes in invertebrates. Ignored are the actions of adrenal and sex steroids-androgens, estrogens, glucocorticoids, mineralocorticoids, and progestins-which act through receptors that arose from an ancestral nuclear receptor in a protochordate. This ligand-based mechanism is unique to vertebrates and was integrated into the already robust network of transcription factors in invertebrates. Adrenal and sex steroids influence almost all aspects of vertebrate differentiation and development. I propose that evolution of this ligand-based mechanism in a primitive vertebrate was an important contribution to vertebrate complexity. Sequencing of genomes from a cephalochordate, such as amphioxus, and from hagfish and lamprey will establish early events in the evolution of steroid hormone signaling, and also allow genetic studies to elucidate how vertebrate complexity depends on steroid hormones.

Hidden Markov models-based system (HMMSPECTR) for detecting structural homologies on the basis of sequential information

HMMSPECTR is a tool for finding putative structural homologs for proteins with known primary sequences. HMMSPECTR contains four major components: a data warehouse with the hidden Markov models (HMM) and alignment libraries; a search program which compares the initial protein sequences with the libraries of HMMs; a secondary structure prediction and comparison program; and a dominant protein selection program that prepares the set of 10-15 "best" proteins from the chosen HMMs. The data warehouse contains four libraries of HMMs. The first two libraries were constructed using different HHM preparation options of the HAMMER program. The third library contains parts ("partial HMM") of initial alignments. The fourth library contains trained HMMs. We tested our program against all of the protein targets proposed in the CASP4 competition. The data warehouse included libraries of structural alignments and HMMs constructed on the basis of proteins publicly available in the Protein Data Bank before the CASP4 meeting. The newest fully automated versions of HMMSPECTR 1.02 and 1.02ss produced better results than the best result reported at CASP4 either by r.m.s.d. or by length (or both) in 64% (HMMSPECTR 1.02) and 79% (HMMSPECTR 1.02ss) of the cases. The improvement is most notable for the targets with complexity 4 (difficult fold recognition cases).

Structure-based phylogenetic analysis of short-chain alcohol dehydrogenases and reclassification of the 17beta-hydroxysteroid dehydrogenase family

Short-chain alcohol dehydrogenases (SCAD) constitute a large and diverse family of ancient origin. Several of its members play an important role in human physiology and disease, especially in the metabolism of steroid substrates (e.g., prostaglandins, estrogens, androgens, and corticosteroids). Their involvement in common human disorders such as endocrine-related cancer, osteoporosis, and Alzheimer disease makes them an important candidate for drug targets. Recent phylogenetic analysis of SCAD is incomplete and does not allow any conclusions on very ancient divergences or on a functional characterization of novel proteins within this complex family. We have developed a 3D structure-based approach to establish the deep-branching pattern within the SCAD family. In this approach, pairwise superpositions of X-ray structures were used to calculate similarity scores as an input for a tree-building algorithm. The resulting phylogeny was validated by comparison with the results of sequence-based algorithms and biochemical data. It was possible to use the 3D data as a template for the reliable determination of the phylogenetic position of novel proteins as a first step toward functional predictions. We were able to discern new patterns in the phylogenetic relationships of the SCAD family, including a basal dichotomy of the 17beta-hydroxysteroid dehydrogenases (17beta-HSDs). These data provide an important contribution toward the development of type-specific inhibitors for 17beta-HSDs for the treatment and prevention of disease. Our structure-based phylogenetic approach can also be applied to increase the reliability of evolutionary reconstructions in other large protein families.

Characterization of multiple Chinese hamster carbonyl reductases

Carbonyl reductase (CR) is an enzyme which can catalyze the oxidoreduction of various carbonyl compounds in the presence of NAD(P)H. With the PCR method, using primers carrying the conserved nucleotide sequence among mammalian CRs, we isolated three different cDNAs (CHCR1, CHCR2 and CHCR3) which encode a unique carbonyl reductase from the Chinese hamster. The PCR products of CHCR1 and CHCR2 were clearly isolated with Bpu1102I, BspEI and XmaI restriction enzymes. The nucleotide-sequence of CHCR3 was completely different from those of CHCR1 and CHCR2. The predicted double-wound betaalphabetaalpha-structures of the CHCRs suggests the presence of a typical NADP(+)-binding motif and is similar to the corresponding region of 3alpha,20beta-hydroxysteroid dehydrogenase and mouse lung tetrameric carbonyl reductase. The deduced amino acid sequence of CHCR1 showed a high homology to CHCR2 (>96%) and the other mammalian CRs (>81%). However, CHCR3 showed a high homology to human CBR3 (>86%) and a relatively lower homology to the other CHCRs (<76%). Bacterial recombinant CHCRs showed typical carbonyl reductase activities towards 4-benzoylpyridine, 4-nitrobenzaldehyde and pyridine 4-carboxyaldehyde. These three CRs showed not only 3-keto reductase of steroids, but also 20-keto reductase. However, these CRs did not show any activity of 17-keto reductase activity. Both CHCR1 and CHCR2 have prostaglandin 9-keto reductase and 15-hydroxyprostaglandin dehydrogenase activities towards PGE(2) and PGF(2alpha) from the analyses of enzymatic reaction products. The results of Western blotting and RT-PCR suggest these CHCRs have a tissue-dependent-distribution in the Chinese hamster.

Evolution of 17beta-hydroxysteroid dehydrogenases and their role in androgen, estrogen and retinoid action

17beta-Hydroxysteroid dehydrogenases (17beta-HSDs) regulate androgen and estrogen concentrations in mammals. By 1995, four distinct enzymes with 17beta-HSD activity had been identified--17beta-HSD-types 1 and 3, which, in vivo, are NADPH-dependent reductases; and 17beta-HSD-types 2 and 4, which, in vivo, are NAD(+)-dependent oxidases. Since then, six additional enzymes with 17beta-HSD activity have been isolated from mammals. With the exception of 17beta-HSD-type 5, which belongs to the aldoketo-reductase (AKR) family, these 17beta-HSDs belong to the short chain dehydrogenase/reductase (SDR) family. Several 17beta-HSDs appear to be examples of convergent evolution. That is, 17beta-HSD activity arose several times from different ancestors. Some 17beta-HSDs share a common ancestor with retinoid oxido-reductases and have retinol dehydrogenase activity. 17beta-HSD-types 2, 6 and 9 appear to have diverged from ancestral retinoid dehydrogenases early in the evolution of deuterostomes during the Cambrian, about 540 million years ago. This coincided with the origin of nuclear receptors for androgens and estrogens suggesting that expression of 17beta-HSDs had an important role in the early evolution of the physiological response to androgens and estrogens.

Prohibitins, stomatins, and plant disease response genes compose a protein superfamily that controls cell proliferation, ion channel regulation, and death

Prohibitins, stomatins, and a group of plant defense response genes are demonstrated to belong to a novel protein superfamily. This superfamily is bound by similar primary and secondary predicted protein structures and hydropathy profiles. A PROSITE-formatted regular expression was generated that is highly predictive for identifying members of this superfamily using PHI-BLAST. The superfamily is named PID (proliferation, ion, and death) because prohibitins are involved in proliferation and cell cycle control, stomatins are involved in ion channel regulation, and the plant defense-related genes are involved in cell death. The plant defense gene family is named HIR (hypersensitive induced reaction) because its members are associated with hypersensitive reactions involving cell death and pathogen resistance. For this study, eight novel maize genes were introduced: four closely related to prohibitins (Zm-phb1, Zm-phb2, Zm-phb3, and Zm-phb4), one to stomatins (Zm-stm1), and three to a gene implicated in plant disease responses (Zm-hir1, Zm-hir2, and Zm-hir3). The maize Zm-hir3 gene transcript is up-regulated in a disease lesion mimic mutation (Les9), supporting a role in maize defense responses. Members of this gene superfamily are involved in diverse functions, but their structural similarity suggests a conserved molecular mechanism, which we postulate to be ion channel regulation.

Sequence analysis of the oxidase/reductase genes upstream of the Rhodococcus erythropolis aldehyde dehydrogenase gene thcA reveals a gene organisation different from Mycobacterium tuberculosis

The sequence of the DNA region upstream of the thiocarbamate-inducible aldehyde dehydrogenase gene thcA of Rhodococcus erythropolis NI86/21 was determined. Most of the predicted ORFs are related to various oxidases/reductases, including short-chain oxidases/reductases, GMC oxidoreductases, alpha-hydroxy acid oxidases (subfamily 1 flavin oxidases/dehydrogenases), and subfamily 2 flavin oxidases/dehydrogenases. One ORF is related to enzymes involved in biosynthesis of PQQ or molybdopterin cofactors. In addition, a putative member of the TetR family of regulatory proteins was identified. The substantial sequence divergence from functionally characterized enzymes precludes a reliable prediction about the probable function of these proteins at this stage. In Mycobacterium tuberculosis H37Rv, most of these ORFs have homologs that are also clustered in the genome, but some striking differences in gene organization were observed between Rhodococcus and Mycobacterium.

Poly-3-hydroxybutyrate degradation in Rhizobium (Sinorhizobium) meliloti: isolation and characterization of a gene encoding 3-hydroxybutyrate dehydrogenase

We have cloned and sequenced the 3-hydroxybutyrate dehydrogenase-encoding gene (bdhA) from Rhizobium (Sinorhizobium) meliloti. The gene has an open reading frame of 777 bp that encodes a polypeptide of 258 amino acid residues (molecular weight 27,177, pI 6.07). The R. meliloti Bdh protein exhibits features common to members of the short-chain alcohol dehydrogenase superfamily. bdhA is the first gene transcribed in an operon that also includes xdhA, encoding xanthine oxidase/dehydrogenase. Transcriptional start site analysis by primer extension identified two transcription starts. S1, a minor start site, was located 46 to 47 nucleotides upstream of the predicted ATG start codon, while S2, the major start site, was mapped 148 nucleotides from the start codon. Analysis of the sequence immediately upstream of either S1 or S2 failed to reveal the presence of any known consensus promoter sequences. Although a sigma54 consensus sequence was identified in the region between S1 and S2, a corresponding transcript was not detected, and a rpoN mutant of R. meliloti was able to utilize 3-hydroxybutyrate as a sole carbon source. The R. meliloti bdhA gene is able to confer upon Escherichia coli the ability to utilize 3-hydroxybutyrate as a sole carbon source. An R. meliloti bdhA mutant accumulates poly-3-hydroxybutyrate to the same extent as the wild type and shows no symbiotic defects. Studies with a strain carrying a lacZ transcriptional fusion to bdhA demonstrated that gene expression is growth phase associated.

Evolution of mammalian 11beta- and 17beta-hydroxysteroid dehydrogenases-type 2 and retinol dehydrogenases from ancestors in Caenorhabditis elegans and evidence for horizontal transfer of a eukaryote dehydrogenase to E. coli

Physiological responses due to steroid hormones and retinoids are regulated by their cognate receptors and dehydrogenases. The origins of either regulatory mechanism are not fully understood. Here we examine the origins of the human 11beta-hydroxysteroid dehydrogenase-type 2, which regulates access of glucocorticoids to cells, and 17beta-hydroxysteroid dehydrogenase-type 2, which regulates access of androgens and estrogens to cells. Sequence comparisons trace their ancestry to homologs in Caenorhabditis elegans. These C. elegans proteins most closely resemble mammalian all-trans and 11-cis-retinol dehydrogenases. The similarity is sufficient -37% to 43% identity to suggest that one or more of the C. elegans homologs metabolizes a retinoid. Receptors for retinoids, but not for androgens, estrogens or glucocorticoids have been identified in C. elegans, suggesting that retinoid-mediated gene transcription is more ancient than that for adrenal and sex steroids. We propose that the hydroxysteroid dehydrogenase-type 2 mechanism for regulating the androgen, estrogen and glucocorticoid concentrations in mammals descended from that for regulating retinoid concentrations. Interestingly, E. coli contains a protein with strong sequence similarity to mammalian retinol dehydrogenases. Sequence comparisons and phylogenetic analysis indicate that the E. coli protein may be an example of horizontal transfer from a eukaryote ancestor.

Spinach CSP41, an mRNA-binding protein and ribonuclease, is homologous to nucleotide-sugar epimerases and hydroxysteroid dehydrogenases

Spinach CSP41 is part of a protein complex that binds to the 3' untranslated region (UTR) of petD precursor-mRNA, a chloroplast gene encoding subunit IV of the cytochrome b6/f complex. CSP41 cleaves the 3'-UTR of petD mRNA within the stem-loop structure, suggesting a key role in the control of chloroplast mRNA stability. We discovered that CSP41 is homologous to nucleotide-sugar epimerases and hydroxysteroid dehydrogenases while seeking distant homologs of these enzymes with a hidden Markov model-based search of Genpept. This analysis identified Synechocystis ORF, Accession 1652543 as a homolog. Subsequent analyses show that spinach CSP41 and Arabidopsis thaliana 2765081 are homologous to the Synechocystis ORF. Information from the solved 3D structures of epimerases and dehydrogenases and our motif analysis of these enzymes is used to predict domains on CSP41 that are important in binding and metabolism of mRNA. Cyanobacteria are among the earliest life forms, indicating that the divergence from a common ancestor of nucleotide-sugar epimerases and an mRNA binding protein with ribonuclease activity was ancient.

Statistical modelling and phylogenetic analysis of a deaminase domain

Deamination reactions are catalyzed by a variety of enzymes including those involved in nucleoside/nucleotide metabolism and cytosine to uracil (C-->U) and adenosine to inosine (A-->I) mRNA editing. The active site of the deaminase (DM) domain in these enzymes contains a conserved histidine (or rarely cysteine), two cysteines and a glutamate proposed to act as a proton shuttle during deamination. Here, a statistical model, a hidden Markov model (HMM), of the DM domain has been created which identifies currently known DM domains and suggests new DM domains in viral, bacterial and eucaryotic proteins. However, no DM domains were identified in the currently predicted proteins from the archaeon Methanococcus jannaschii and possible causes for, and a potential means to ameliorate this situation are discussed. In some of the newly identified DM domains, the glutamate is changed to a residue that could not function as a proton shuttle and in one instance (Mus musculus spermatid protein TENR) the cysteines are also changed to lysine and serine. These may be non-competent DM domains able to bind but not act upon their substrate. Phylogenetic analysis using an HMM-generated alignment of DM domains reveals three branches with clear substructure in each branch. The results suggest DM domains that are candidates for yeast, platyhelminth, plant and mammalian C-->U and A-->I mRNA editing enzymes. Some bacterial and eucaryotic DM domains form distinct branches in the phylogenetic tree suggesting the existence of common, novel substrates.

The Fanconi anemia complementation group A protein contains a peroxidase domain

Computational analysis of the Fanconi anemia (FA) complementation group A protein suggests that it contains a peroxidase domain. FA proteins may be part of a general mechanism that protects cells from oxidative damage.

The proofreading domain of Escherichia coli DNA polymerase I and other DNA and/or RNA exonuclease domains

Prior sequence analysis studies have suggested that bacterial ribonuclease (RNase) Ds comprise a complete domain that is found also in Homo sapiens polymyositis-scleroderma overlap syndrome 100 kDa autoantigen and Werner syndrome protein. This RNase D 3'-->5' exoribonuclease domain was predicted to have a structure and mechanism of action similar to the 3'-->5' exodeoxyibonuclease (proofreading) domain of DNA polymerases. Here, hidden Markov model (HMM) and phylogenetic studies have been used to identify and characterise other sequences that may possess this exonuclease domain. Results indicate that it is also present in the RNase T family; Borrelia burgdorferi P93 protein, an immunodominant antigen in Lyme disease; bacteriophage T4 dexA and Escherichia coli exonuclease I, processive 3'-->5' exodeoxyribonucleases that degrade single-stranded DNA; Bacillus subtilis dinG, a probable helicase involved in DNA repair and possibly replication, and peptide synthase 1; Saccharomyces cerevisiae Pab1p-dependent poly(A) nuclease PAN2 subunit, required for shortening mRNA poly(A) tails; Caenorhabditis elegans and Mus musculus CAF1, transcription factor CCR4-associated factor 1; Xenopus laevis XPMC2, prevention of mitotic catastrophe in fission yeast; Drosophila melanogaster egalitarian, oocyte specification and axis determination, and exuperantia, establishment of oocyte polarity; H.sapiens HEM45, expressed in tumour cell lines and uterus and regulated by oestrogen; and 31 open reading frames including one in Methanococcus jannaschii . Examination of a multiple sequence alignment and two three-dimensional structures of proofreading domains has allowed definition of the core sequence, structural and functional elements of this exonuclease domain.

Statistical modeling and analysis of the LAGLIDADG family of site-specific endonucleases and identification of an intein that encodes a site-specific endonuclease of the HNH family

The LAGLIDADG and HNH families of site-specific DNA endonucleases encoded by viruses, bacteriophages as well as archaeal, eucaryotic nuclear and organellar genomes are characterized by the sequence motifs 'LAGLIDADG' and 'HNH', respectively. These endonucleases have been shown to occur in different environments: LAGLIDADG endonucleases are found in inteins, archaeal and group I introns and as free standing open reading frames (ORFs); HNH endonucleases occur in group I and group II introns and as ORFs. Here, statistical models (hidden Markov models, HMMs) that encompass both the conserved motifs and more variable regions of these families have been created and employed to characterize known and potential new family members. A number of new, putative LAGLIDADG and HNH endonucleases have been identified including an intein-encoded HNH sequence. Analysis of an HMM-generated multiple alignment of 130 LAGLIDADG family members and the three-dimensional structure of the I- Cre I endonuclease has enabled definition of the core elements of the repeated domain (approximately 90 residues) that is present in this family of proteins. A conserved negatively charged residue is proposed to be involved in catalysis. Phylogenetic analysis of the two families indicates a lack of exchange of endonucleases between different mobile elements (environments) and between hosts from different phylogenetic kingdoms. However, there does appear to have been considerable exchange of endonuclease domains amongst elements of the same type. Such events are suggested to be important for the formation of elements of new specficity.