THUMP--a predicted RNA-binding domain shared by 4-thiouridine, pseudouridine synthases and RNA methylases

Sequence profile searches were used to identify an ancient domain in ThiI-like thiouridine synthases, conserved RNA methylases, archaeal pseudouridine synthases and several uncharacterized proteins. We predict that this domain is an RNA-binding domain that adopts an alpha/beta fold similar to that found in the C-terminal domain of translation initiation factor 3 and ribosomal protein S8.     

Structural insights into the interaction of the Nip7 PUA domain with polyuridine RNA

The conserved protein Nip7 is involved in ribosome biogenesis, being required for proper 27S pre-rRNA processing and 60S ribosome subunit assembly in Saccharomyces cerevisiae. Yeast Nip7p interacts with nucleolar proteins and with the exosome subunit Rrp43p, but its molecular function remains to be determined. Solution of the Pyrococcus abyssi Nip7 (PaNip7) crystal structure revealed a monomeric protein composed by two alpha-beta domains. The N-terminal domain is formed by a five-stranded antiparallel beta-sheet surrounded by three alpha-helices and a 310 helix while the C-terminal, a mixed beta-sheet domain composed by strands beta8 to beta12, one alpha-helix, and a 310 helix, corresponds to the conserved PUA domain (after Pseudo-Uridine synthases and Archaeosine-specific transglycosylases). By combining structural analyses and RNA interaction assays, we assessed the ability of both yeast and archaeal Nip7 orthologues to interact with RNA. Structural alignment of the PaNip7 PUA domain with the RNA-interacting surface of the ArcTGT (archaeosine tRNA-guanine transglycosylase) PUA domain indicated that in the archaeal PUA domain positively charged residues (R151, R152, K155, and K158) are involved in RNA interaction. However, equivalent positions are occupied by mostly hydrophobic residues (A/G160, I161, F164, and A167) in eukaryotic Nip7 orthologues. Both proteins can bind specifically to polyuridine, and RNA interaction requires specific residues of the PUA domain as determined by site-directed mutagenesis. This work provides experimental verification that the PUA domain mediates Nip7 interaction with RNA and reveals that the preference for interaction with polyuridine sequences is conserved in Archaea and eukaryotic Nip7 proteins.     

Evolutionary Comparisons of RecA-Like Proteins Across All Major Kingdoms of Living Organisms

Protein sequences with similarities to Escherichia coli RecA were compared across the major kingdoms of eubacteria, archaebacteria, and eukaryotes. The archaeal sequences branch monophyletically and are most closely related to the eukaryotic paralogous Rad51 and Dmc1 groups. A multiple alignment of the sequences suggests a modular structure of RecA-like proteins consisting of distinct segments, some of which are conserved only within subgroups of sequences. The eukaryotic and archaeal sequences share an N-terminal domain which may play a role in interactions with other factors and nucleic acids. Several positions in the alignment blocks are highly conserved within the eubacteria as one group and within the eukaryotes and archaebacteria as a second group, but compared between the groups these positions display nonconservative amino acid substitutions. Conservation within the RecA-like core domain identifies possible key residues involved in ATP-induced conformational changes. We propose that RecA-like proteins derive evolutionarily from an assortment of independent domains and that the functional homologs of RecA in noneubacteria comprise an array of RecA-like proteins acting in series or cooperatively. Received: 25 October 1996 / Accepted: 31 December 1996     

Close Evolutionary Relatedness of α-Amylases from Archaea and Plants

The amino acid sequences of 22 α-amylases from family 13 of glycosyl hydrolases were analyzed with the aim of revealing the evolutionary relationships between the archaeal α-amylases and their eubacterial and eukaryotic counterparts. Two evolutionary distance trees were constructed: (i) the first one based on the alignment of extracted best-conserved sequence regions (58 residues) comprising β2, β3, β4, β5, β7, and β8 strand segments of the catalytic (α/β)₈-barrel and a short conserved stretch in domain B protruding out of the barrel in the β3 →α3 loop, and (ii) the second one based on the alignment of the substantial continuous part of the (α/β)₈-barrel involving the entire domain B (consensus length: 386 residues). With regard to archaeal α-amylases, both trees compared brought, in fact, the same results; i.e., all family 13 α-amylases from domain Archaea were clustered with barley pI isozymes, which represent all plant α-amylases. The enzymes from Bacillus licheniformis and Escherichia coli, representing liquefying and cytoplasmic α-amylases, respectively, seem to be the further closest relatives to archaeal α-amylases. This evolutionary relatedness clearly reflects the discussed similarities in the amino acid sequences of these α-amylases, especially in the best-conserved sequence regions. Since the results for α-amylases belonging to all three domains (Eucarya, Eubacteria, Archaea) offered by both evolutionary trees are very similar, it is proposed that the investigated conserved sequence regions may indeed constitute the ``sequence fingerprints' of a given α-amylase. Received: 3 June 1998 / Accepted: 20 August 1998     

The PUA domain - a structural and functional overview

The pseudouridine synthase and archaeosine transglycosylase (PUA) domain is a compact and highly conserved RNA-binding motif that is widespread among diverse types of proteins from the three kingdoms of life. Its three-dimensional architecture is well established, and the structures of several PUA-RNA complexes reveal a common RNA recognition surface, but also some versatility in the way in which the motif binds to RNA. The PUA domain is often part of RNA modification enzymes and ribonucleoproteins, but it has also been unexpectedly found fused to enzymes involved in proline biosynthesis, where it plays an unknown role. The functional impact of the domain varies with the protein studied, ranging from minor to essential effects. PUA motifs are involved in dyskeratosis congenita and cancer, pointing to links between RNA metabolism and human diseases.     

An HflX-type GTPase from Sulfolobus solfataricus binds to the 50S ribosomal subunit in all nucleotide-bound states

HflX GTPases are found in all three domains of life, the Bacteria, Archaea, and Eukarya. HflX from Escherichia coli has been shown to bind to the 50S ribosomal subunit in a nucleotide-dependent manner, and this interaction strongly stimulates its GTPase activity. We recently determined the structure of an HflX ortholog from the archaeon Sulfolobus solfataricus (SsoHflX). It revealed the presence of a novel HflX domain that might function in RNA binding and is linked to a canonical G domain. This domain arrangement is common to all archaeal, bacterial, and eukaryotic HflX GTPases. This paper shows that the archaeal SsoHflX, like its bacterial orthologs, binds to the 50S ribosomal subunit. This interaction does not depend on the presence of guanine nucleotides. The HflX domain is sufficient for ribosome interaction. Binding appears to be restricted to free 50S ribosomal subunits and does not occur with 70S ribosomes engaged in translation. The fingerprint (1)H-(15)N heteronuclear correlation nuclear magnetic resonance (NMR) spectrum of SsoHflX reveals a large number of well-resolved resonances that are broadened upon binding to the 50S ribosomal subunit. The GTPase activity of SsoHflX is stimulated by crude fractions of 50S ribosomal subunits, but this effect is lost with further high-salt purification of the 50S ribosomal subunits, suggesting that the stimulation depends on an extrinsic factor bound to the 50S ribosomal subunit. Our results reveal common properties but also marked differences between archaeal and bacterial HflX proteins.     

Eukaryotic ribonucleases P/MRP: the crystal structure of the P3 domain

Ribonuclease (RNase) P is a site‐specific endoribonuclease found in all kingdoms of life. Typical RNase P consists of a catalytic RNA component and a protein moiety. In the eukaryotes, the RNase P lineage has split into two, giving rise to a closely related enzyme, RNase MRP, which has similar components but has evolved to have different specificities. The eukaryotic RNases P/MRP have acquired an essential helix‐loop‐helix protein‐binding RNA domain P3 that has an important function in eukaryotic enzymes and distinguishes them from bacterial and archaeal RNases P. Here, we present a crystal structure of the P3 RNA domain from Saccharomyces cerevisiae RNase MRP in a complex with RNase P/MRP proteins Pop6 and Pop7 solved to 2.7 Å. The structure suggests similar structural organization of the P3 RNA domains in RNases P/MRP and possible functions of the P3 domains and proteins bound to them in the stabilization of the holoenzymes' structures as well as in interactions with substrates. It provides the first insight into the structural organization of the eukaryotic enzymes of the RNase P/MRP family.     

Identifying bacterial and archaeal homologs of pentameric ligand-gated ion channel (pLGIC) family using domain-based and alignment-based approaches

Identification of bacterial and archaeal counterparts to eukaryotic ion channels has greatly facilitated studies of structural biophysics of the channels. Often, searches based only on sequence alignment tools are inadequate for discovering such distant bacterial and archaeal counterparts. We address the discovery of bacterial and archaeal members of the Pentameric Ligand-Gated Ion Channel (pLGIC) family by a combination of four computational methods. One domain-based method involves retrieval of proteins with pLGIC-relevant domains by matching those domains to previously established domain templates in the InterPro family of databases. The second domain-based method involves searches using ungapped de-novo motifs discovered by MEME which were trained with well characterized members of the pLGIC family. The third and fourth methods involve the use of two sequence alignment search algorithms BLASTp and psiBLAST respectively. The sequences returned from all methods were screened by having the correct topology for pLGIC's, and by returning an annotated member of this family as one of the first ten hits using BLASTp against a comprehensive database of eukaryotic proteins. We found the domain based searches to have high specificity but low sensitivity, while the sequence alignment methods have higher sensitivity but lower specificity. The four methods together discovered 69 putative bacterial and archaeal members of the pLGIC family. We ranked and divide the 69 proteins into groups according to the similarity of their domain compositions with known eukaryotic pLGIC's. One especially notable group is more closely related to eukaryotic pLGIC's than to any other known protein family, and has the overall topology of pLGIC's, but the functional domains they contain are sufficiently different from those found in known pLGIC's that they do not score very well against the pLGIC domain templates. We conclude that multiple methods used in a coordinated fashion outperform any single method for identifying likely distant bacterial and archaeal proteins that may provide useful models for important eukaryotic channel function. We note also that the methods used here are largely standard and readily accessible. The novelty is in the effectiveness of a strategy that combines these methods for identifying bacterial and archea relatives of this family. Therefore the paper may serve as a template for a broad group of workers to reliably identify bacterial and archaeal counterparts to eukaryotic proteins.     

Symbiosis Between Methanogenic Archaea and δ-Proteobacteria as the Origin of Eukaryotes: The Syntrophic Hypothesis

We present a novel hypothesis for the origin of the eukaryotic cell, or eukaryogenesis, based on a metabolic symbiosis (syntrophy) between a methanogenic archaeon (methanobacterial-like) and a δ-proteobacterium (an ancestral sulfate-reducing myxobacterium). This syntrophic symbiosis was originally mediated by interspecies H₂ transfer in anaerobic, possibly moderately thermophilic, environments. During eukaryogenesis, progressive cellular and genomic cointegration of both types of prokaryotic partners occurred. Initially, the establishment of permanent consortia, accompanied by extensive membrane development and close cell–cell interactions, led to a highly evolved symbiotic structure already endowed with some primitive eukaryotic features, such as a complex membrane system defining a protonuclear space (corresponding to the archaeal cytoplasm), and a protoplasmic region (derived from fusion of the surrounding bacterial cells). Simultaneously, bacterial-to-archaeal preferential gene transfer and eventual replacement took place. Bacterial genome extinction was thus accomplished by gradual transfer to the archaeal host, where genes adapted to a new genetic environment. Emerging eukaryotes would have inherited archaeal genome organization and dynamics and, consequently, most DNA-processing information systems. Conversely, primordial genes for social and developmental behavior would have been provided by the ancient myxobacterial symbiont. Metabolism would have been issued mainly from the versatile bacterial organotrophy, and progressively, methanogenesis was lost. Received: 5 January 1998 / Accepted: 18 March 1998     

A Protein Related to Eucaryal and Bacterial DNA-Motor Proteins in the Hyperthermophilic Archaeon Sulfolobus acidocaldarius

We have isolated a new gene encoding a putative 103-kDa protein from the hyperthermophilic archaeon Sulfolobus acidocaldarius. Analysis of the deduced amino-acid sequence shows an extended central domain, predicted to form coiled-coil structures, and two terminal domains that display purine NTPase motifs. These features are reminiscent of mechanochemical motor proteins which use the energy of ATP hydrolysis to move specific cellular components. Comparative analysis of the amino-acid sequence of the terminal domains and predicted structural organization of this putative purine NTPase show that it is related both to eucaryal proteins from the ``SMC family' involved in the condensation of chromosomes and to several bacterial and eucaryal proteins involved in DNA recombination/repair. Further analyses revealed that these proteins are all members of the so called ``UvrA-related NTP-binding proteins superfamily' and form a large subgroup of motor-like NTPases involved in different DNA processing mechanisms. The presence of such protein in Archaea, Bacteria, and Eucarya suggests an early origin of DNA-motor proteins that could have emerged and diversified by domain shuffling. Received: 29 June 1996 / Accepted: 28 February 1997     

DNA Repair Systems in Archaea: Mementos from the Last Universal Common Ancestor?

DNA repair in the Archaea is relevant to the consideration of genome maintenance and replication fidelity in the last universal common ancestor (LUCA) from two perspectives. First, these prokaryotes embody a mix of bacterial and eukaryal molecular features. Second, DNA repair proteins would have been essential in LUCA to maintain genome integrity, regardless of the environmental temperature. Yet we know very little of the basic molecular mechanisms of DNA damage and repair in the Archaea in general. Many studies on DNA repair in archaea have been conducted with hyperthermophiles because of the additional stress imposed on their macromolecules by high temperatures. In addition, of the six complete archaeal genome sequences published so far, five are thermophilic archaea. We have recently shown that the hyperthermophile Pyrococcus furiosus has an extraordinarily high capacity for repair of radiation-induced double-strand breaks and we have identified and sequenced several genes involved in DNA repair in P. furiosus. At the sequence level, only a few genes share homology with known bacterial repair genes. For instance, our phylogenetic analysis indicates that archaeal recombinases occur in two paralogous gene families, one of which is very deeply branched, and both recombinases are more closely related to the eukaryotic RAD51 and Dmc1 gene families than to the Escherichia coli recA gene. We have also identified a gene encoding a repair endo/exonuclease in the genomes of several Archaea. The archaeal sequences are highly homologous to those of the eukaryotic Rad2 family and they cluster with genes of the FEN-1 subfamily, which are known to be involved in DNA replication and repair in eukaryotes. We argue that there is a commonality of mechanisms and protein sequences, shared between prokaryotes and eukaryotes for several modes of DNA repair, reflecting diversification from a minimal set of genes thought to represent the genome of the LUCA.     

Evolution of DNA Polymerase Families: Evidences for Multiple Gene Exchange Between Cellular and Viral Proteins

A phylogenetic analysis of the five major families of DNA polymerase is presented. Viral and plasmid sequences are included in this compilation along with cellular enzymes. The classification by Ito and Braithwaite (Ito and Braithwaite 1991) of the A, B, C, D, and X families has been extended to accommodate the ``Y family' of DNA polymerases that are related to the eukaryotic RAD30 and the bacterial UmuC gene products. After analysis, our data suggest that no DNA polymerase family was universally conserved among the three biological domains and no simple evolutionary scenario could explain that observation. Furthermore, viruses and plasmids carry a remarkably diverse set of DNA polymerase genes, suggesting that lateral gene transfer is frequent and includes non-orthologous gene displacements between cells and viruses. The relationships between viral and host genes appear very complex. We propose that the gamma DNA polymerase of the mitochondrion replication apparatus is of phage origin and that this gene replaced the one in the bacterial ancestor. Often there was no obvious relation between the viral and the host DNA polymerase, but an interesting exception concerned the family B enzymes: in which ancient gene exchange can be detected between the viruses and their hosts. Additional evidence for horizontal gene transfers between cells and viruses comes from an analysis of the small damage-inducible DNA polymerases. Taken together, these findings suggest a complex evolutionary history of the DNA replication apparatus that involved significant exchanges between viruses, plasmids, and their hosts.     

Crystal structure of the catalytic domain of RluD, the only rRNA pseudouridine synthase required for normal growth of Escherichia coli

Escherichia coli pseudouridine synthase RluD makes pseudouridines 1911, 1915, and 1917 in the loop of helix 69 in 23S RNA. These are the most highly conserved ribosomal pseudouridines known. Of 11 pseudouridine synthases in E. coli, only cells lacking RluD have severe growth defects and abnormal ribosomes. We have determined the 2.0 A structure of the catalytic domain of RluD (residues 77-326), the first structure of an RluA family member. The catalytic domain folds into a mainly antiparallel beta-sheet flanked by several loops and helices. A positively charged cleft that presumably binds RNA leads to the conserved Asp 139. The RluD N-terminal S4 domain, connected by a flexible linker, is disordered in our structure. RluD is very similar in both catalytic domain structure and active site arrangement to the pseudouridine synthases RsuA, TruB, and TruA. We identify five sequence motifs, two of which are novel, in the RluA, RsuA, TruB, and TruA families, uniting them as one superfamily. These results strongly suggest that four of the five families of pseudouridine synthases arose by divergent evolution. The RluD structure also provides insight into its multisite specificity.     

Identifying bacterial and archaeal homologs of pentameric ligand-gated ion channel (pLGIC) family using domain-based and alignment-based approaches

Identification of bacterial and archaeal counterparts to eukaryotic ion channels has greatly facilitated studies of structural biophysics of the channels. Often, searches based only on sequence alignment tools are inadequate for discovering such distant bacterial and archaeal counterparts. We address the discovery of bacterial and archaeal members of the Pentameric Ligand-Gated Ion Channel (pLGIC) family by a combination of four computational methods. One domain-based method involves retrieval of proteins with pLGIC-relevant domains by matching those domains to previously established domain templates in the InterPro family of databases. The second domain-based method involves searches using ungapped de-novo motifs discovered by MEME which were trained with well characterized members of the pLGIC family. The third and fourth methods involve the use of two sequence alignment search algorithms BLASTp and psiBLAST respectively. The sequences returned from all methods were screened by having the correct topology for pLGIC's, and by returning an annotated member of this family as one of the first ten hits using BLASTp against a comprehensive database of eukaryotic proteins. We found the domain based searches to have high specificity but low sensitivity, while the sequence alignment methods have higher sensitivity but lower specificity. The four methods together discovered 69 putative bacterial and archaeal members of the pLGIC family. We ranked and divide the 69 proteins into groups according to the similarity of their domain compositions with known eukaryotic pLGIC's. One especially notable group is more closely related to eukaryotic pLGIC's than to any other known protein family, and has the overall topology of pLGIC's, but the functional domains they contain are sufficiently different from those found in known pLGIC's that they do not score very well against the pLGIC domain templates. We conclude that multiple methods used in a coordinated fashion outperform any single method for identifying likely distant bacterial and archaeal proteins that may provide useful models for important eukaryotic channel function. We note also that the methods used here are largely standard and readily accessible. The novelty is in the effectiveness of a strategy that combines these methods for identifying bacterial and archea relatives of this family. Therefore the paper may serve as a template for a broad group of workers to reliably identify bacterial and archaeal counterparts to eukaryotic proteins.     

Families of Proteins Forming Transmembrane Channels

Channel-forming proteins/peptides fall into over 100 currently recognized families, most of which are restricted to prokaryotes or eukaryotes, but a few of which are ubiquitous. These proteins fall into three major currently recognized classes: (i) α-helix-type channels present in bacterial, archaeal and eukaryotic cytoplasmic and organellar membranes, (ii) β-barrel-type porins present in the outer membranes of Gram-negative bacterial cells, mitochondria and chloroplasts, and (iii) protein/peptide toxins targeted to the cytoplasmic membranes of cells other than those that synthesize the toxins. High-resolution 3-dimensional structural data are available for representative proteins/peptides of all three of these channel-forming types. Each type exhibits distinctive features that distinguish them from the other channel protein types and from carriers. Structural, functional, and evolutionary aspects of transmembrane channel-formers are discussed. Received: 10 September 1999/Revised: 11 February 2000     

Taxonomic distribution,repeats, and functions of the S1 domain‐containing proteins as members of the OB‐fold family

Proteins of the nucleic acid‐binding proteins superfamily perform such functions as processing, transport, storage, stretching, translation, and degradation of RNA. It is one of the 16 superfamilies containing the OB‐fold in protein structures. Here, we have analyzed the superfamily of nucleic acid‐binding proteins (the number of sequences exceeds 200,000) and obtained that this superfamily prevalently consists of proteins containing the cold shock DNA‐binding domain (ca. 131,000 protein sequences). Proteins containing the S1 domain compose 57% from the cold shock DNA‐binding domain family. Furthermore, we have found that the S1 domain was identified mainly in the bacterial proteins (ca. 83%) compared to the eukaryotic and archaeal proteins, which are available in the UniProt database. We have found that the number of multiple repeats of S1 domain in the S1 domain‐containing proteins depends on the taxonomic affiliation. All archaeal proteins contain one copy of the S1 domain, while the number of repeats in the eukaryotic proteins varies between 1 and 15 and correlates with the protein size. In the bacterial proteins, the number of repeats is no more than 6, regardless of the protein size. The large variation of the repeat number of S1 domain as one of the structural variants of the OB‐fold is a distinctive feature of S1 domain‐containing proteins. Proteins from the other families and superfamilies have either one OB‐fold or change slightly the repeat numbers. On the whole, it can be supposed that the repeat number is a vital for multifunctional activity of the S1 domain‐containing proteins. Proteins 2017; 85:602–613. © 2016 Wiley Periodicals, Inc.