Orientation of the transcription preinitiation complex in Archaea

The basal transcription machinery of Archaea corresponds to the minimal subset of factors required for RNA polymerase II transcription in eukaryotes. Using just two factors, Archaea recruit the RNA polymerase to promoters and define the direction of transcription. Notably, the principal determinant for the orientation of transcription is not the recognition of the TATA box by the TATA-box-binding protein. Instead, transcriptional polarity is governed by the interaction of the archaeal TFIIB homologue with a conserved motif immediately upstream of the TATA box. This interaction yields an archaeal preinitiation complex with the same orientation as the analogous eukaryal complex.

The structural basis for the oriented assembly of a TBP/TFB/promoter complex

Recently the definition of the metazoan RNA polymerase II and archaeal core promoters has been expanded to include a region immediately upstream of the TATA box called the B recognition element (BRE), so named because eukaryal transcription factor TFIIB and its archaeal orthologue TFB interact with the element in a sequence-specific manner. Here we present the 2.4-Å crystal structure of archaeal TBP and the C-terminal core of TFB (TFBc) in a complex with an extended TATA-box-containing promoter that provides a detailed picture of the stereospecific interactions between the BRE and a helix–turn–helix motif in the C-terminal cyclin repeat of TFBc. This interaction is important in determining the level of basal transcription and explicitly defines the direction of transcription.

A functional homolog of a yeast tRNA splicing enzyme is conserved in higher eukaryotes and in Escherichia coli

tRNA splicing in the yeast Saccharomyces cerevisiae requires an endonuclease to excise the intron, tRNA ligase to join the tRNA half-molecules, and 2′-phosphotransferase to transfer the splice junction 2′-phosphate from ligated tRNA to NAD, producing ADP ribose 1′′–2′′ cyclic phosphate (Appr>p). We show here that functional 2′-phosphotransferases are found throughout eukaryotes, occurring in two widely divergent yeasts (Candida albicans and Schizosaccharomyces pombe), a plant (Arabidopsis thaliana), and mammals (Mus musculus); this finding is consistent with a role for the enzyme, acting in concert with ligase, to splice tRNA or other RNA molecules. Surprisingly, functional 2′-phosphotransferase is found also in the bacterium Escherichia coli, which does not have any known introns of this class, and does not appear to have a ligase that generates junctions with a 2′-phosphate. Analysis of the database shows that likely members of the 2′-phosphotransferase family are found also in one other bacterium (Pseudomonas aeruginosa) and two archaeal species (Archaeoglobus fulgidus and Pyrococcus horikoshii). Phylogenetic analysis reveals no evidence for recent horizontal transfer of the 2′-phosphotransferase into Eubacteria, suggesting that the 2′-phosphotransferase has been present there since close to the time that the three kingdoms diverged. Although 2′-phosphotransferase is not present in all Eubacteria, and a gene disruption experiment demonstrates that the protein is not essential in E. coli, the continued presence of 2′-phosphotransferase in Eubacteria over large evolutionary times argues for an important role for the protein.

Chloroplast small heat shock proteins: Evidence for atypical evolution of an organelle-localized protein

Knowledge of the origin and evolution of gene families is critical to our understanding of the evolution of protein function. To gain a detailed understanding of the evolution of the small heat shock proteins (sHSPs) in plants, we have examined the evolutionary history of the chloroplast (CP)-localized sHSPs. Previously, these nuclear-encoded CP proteins had been identified only from angiosperms. This study reveals the presence of the CP sHSPs in a moss, Funaria hygrometrica. Two clones for CP sHSPs were isolated from a F. hygrometrica heat shock cDNA library that represent two distinct CP sHSP genes. Our analysis of the CP sHSPs reveals unexpected evolutionary relationships and patterns of sequence conservation. Phylogenetic analysis of the CP sHSPs with other plant CP sHSPs and eukaryotic, archaeal, and bacterial sHSPs shows that the CP sHSPs are not closely related to the cyanobacterial sHSPs. Thus, they most likely evolved via gene duplication from a nuclear-encoded cytosolic sHSP and not via gene transfer from the CP endosymbiont. Previous sequence analysis had shown that all angiosperm CP sHSPs possess a methionine-rich region in the N-terminal domain. The primary sequence of this region is not highly conserved in the F. hygrometrica CP sHSPs. This lack of sequence conservation indicates that sometime in land plant evolution, after the divergence of mosses from the common ancestor of angiosperms but before the monocot–dicot divergence, there was a change in the selective constraints acting on the CP sHSPs.

Identification and characterization of a single-stranded DNA-binding protein from the archaeon Methanococcus jannaschii

Single-stranded DNA-binding proteins (SSBs) play essential roles in DNA replication, recombination, and repair in bacteria and eukarya. We report here the identification and characterization of the SSB of an archaeon, Methanococcus jannaschii. The M. jannaschii SSB (mjaSSB) has significant amino acid sequence similarity to the eukaryotic SSB, replication protein A (RPA), and contains four tandem repeats of the core single-stranded DNA (ssDNA) binding domain originally defined by structural studies of RPA. Homologous SSBs are encoded by the genomes of other archaeal species, including Methanobacterium thermoautotrophicum and Archaeoglobus fulgidus. The purified mjaSSB binds to ssDNA with high affinity and selectivity. The apparent association constant for binding to ssDNA is similar to that of RPA under comparable experimental conditions, and the affinity for ssDNA exceeds that for double-stranded DNA by at least two orders of magnitude. The binding site size for mjaSSB is ≈20 nucleotides. Given that RPA is related to mjaSSB at the sequence level and to Escherichia coli SSB at the structural level, we conclude that the SSBs of archaea, eukarya, and bacteria share a common core ssDNA-binding domain. This ssDNA-binding domain was presumably present in the common ancestor to all three major branches of life.

A dual-specificity aminoacyl-tRNA synthetase in the deep-rooted eukaryote Giardia lamblia

Cysteinyl-tRNA (Cys-tRNA) is essential for protein synthesis. In most organisms the enzyme responsible for the formation of Cys-tRNA is cysteinyl-tRNA synthetase (CysRS). The only known exceptions are the euryarchaea Methanococcus jannaschii and Methanobacterium thermoautotrophicum, which do not encode a CysRS. Deviating from the accepted concept of one aminoacyl-tRNA synthetase per amino acid, these organisms employ prolyl-tRNA synthetase as the enzyme that carries out Cys-tRNA formation. To date this dual-specificity prolyl-cysteinyl-tRNA synthetase (ProCysRS) is only known to exist in archaea. Analysis of the preliminary genomic sequence of the primitive eukaryote Giardia lamblia indicated the presence of an archaeal prolyl-tRNA synthetase (ProRS). Its proS gene was cloned and the gene product overexpressed in Escherichia coli. By using G. lamblia, M. jannaschii, or E. coli tRNA as substrate, this ProRS was able to form Cys-tRNA and Pro-tRNA in vitro. Cys-AMP formation, but not Pro-AMP synthesis, was tRNA-dependent. The in vitro data were confirmed in vivo, as the cloned G. lamblia proS gene was able to complement a temperature-sensitive E. coli cysS strain. Inhibition studies of CysRS activity with proline analogs (thiaproline and 5′-O-[N-(l-prolyl)-sulfamoyl]adenosine) in a Giardia S-100 extract predicted that the organism also contains a canonical CysRS. This prediction was confirmed by cloning and analysis of the corresponding cysS gene. Like a number of archaea, Giardia contains two enzymes, ProCysRS and CysRS, for Cys-tRNA formation. In contrast, the purified Saccharomyces cerevisiae and E. coli ProRS enzymes were unable to form Cys-tRNA under these conditions. Thus, the dual specificity is restricted to the archaeal genre of ProRS. G. lamblia's archaeal-type prolyl- and alanyl-tRNA synthetases refine our understanding of the evolution and interaction of archaeal and eukaryal translation systems.

Class I release factors in ciliates with variant genetic codes

In eukaryotes with the universal genetic code a single class I release factor (eRF1) most probably recognizes all stop codons (UAA, UAG and UGA) and is essential for termination of nascent peptide synthesis. It is well established that stop codons have been reassigned to amino acid codons at least three times among ciliates. The codon specificities of ciliate eRF1s must have been modified to accommodate the variant codes. In this study we have amplified, cloned and sequenced eRF1 genes of two hypotrichous ciliates, Oxytricha trifallax (UAA and UAG for Gln) and Euplotes aediculatus (UGA for Cys). We also sequenced/identified three protist and two archaeal class I RF genes to enlarge the database of eRF1/aRF1s with the universal code. Extensive comparisons between universal code eRF1s and those of Oxytricha, Euplotes and Tetrahymena, which represent three lineages that acquired variant codes independently, provide important clues to identify stop codon-binding regions in eRF1. Domain 1 in the five ciliate eRF1s, particulary the TASNIKS heptapeptide and its adjacent region, differs significantly from domain 1 in universal code eRF1s. This observation suggests that domain 1 contains the codon recognition site, but that the mechanism of eRF1 codon recognition may be more complex than proposed by Nakamura et al. or Knight and Landweber.

Prediction of operons in microbial genomes

Operon structure is an important organization feature of bacterial genomes. Many sets of genes occur in the same order on multiple genomes; these conserved gene groupings represent candidate operons. This study describes a computational method to estimate the likelihood that such conserved gene sets form operons. The method was used to analyze 34 bacterial and archaeal genomes, and yielded more than 7600 pairs of genes that are highly likely (P ≥ 0.98) to belong to the same operon. The sensitivity of our method is 30–50% for the Escherichia coli genome. The predicted gene pairs are available from our World Wide Web site http://www.tigr.org/tigr-scripts/operons/operons.cgi.

REBASE—restriction enzymes and methylases

REBASE contains comprehensive information about restriction enzymes, DNA methylases and related proteins such as nicking enzymes, specificity subunits and control proteins. It contains published and unpublished references, recognition and cleavage sites, isoschizomers, commercial availability, methy­lation sensitivity, crystal data and sequence data. Homing endonucleases are also included. Most recently, extensive information about the methy­lation sensitivity of restriction enzymes has been added and a new feature contains complete analyses of the putative restriction systems in the sequenced bacterial and archaeal genomes. The data is distributed via email, ftp (ftp.neb.com) and the Web (http://rebase.neb.com).

The Sm domain is an ancient RNA-binding motif with oligo(U) specificity

Sm and Sm-like proteins are members of a family of small proteins that is widespread throughout eukaryotic kingdoms. These proteins form heteromers with one another and bind, as heteromeric complexes, to various RNAs, recognizing primarily short U-rich stretches. Interestingly, completion of several genome projects revealed that archaea also contain genes that may encode Sm-like proteins. Herein, we studied the properties of one Sm-like protein derived from the archaebacterium Archaeoglobus fulgidus and overexpressed in Escherichia coli. This single small protein closely reflects the properties of an Sm or Sm-like protein heteromer. It binds to RNA with a high specificity for oligo(U), and assembles onto the RNA to form a complex that exhibits, as judged by electron microscopy, a ring-like structure similar to the ones observed with the Sm core ribonucleoprotein and the like Sm (LSm) protein heteromer. Importantly, multivariate statistical analysis of negative-stain electron-microscopic images revealed a sevenfold symmetry for the observed ring structure, indicating that the proteins form a homoheptamer. These results support the structural model of the Sm proteins derived from crystallographic studies on Sm heterodimers and demonstrate that the Sm protein family evolved from a single ancestor that was present before the eukaryotic and archaeal kingdoms separated.

A novel site of antibiotic action in the ribosome: Interaction of evernimicin with the large ribosomal subunit

Evernimicin (Evn), an oligosaccharide antibiotic, interacts with the large ribosomal subunit and inhibits bacterial protein synthesis. RNA probing demonstrated that the drug protects a specific set of nucleotides in the loops of hairpins 89 and 91 of 23S rRNA in bacterial and archaeal ribosomes. Spontaneous Evn-resistant mutants of Halobacterium halobium contained mutations in hairpins 89 and 91 of 23S rRNA. In the ribosome tertiary structure, rRNA residues involved in interaction with the drug form a tight cluster that delineates the drug-binding site. Resistance mutations in the bacterial ribosomal protein L16, which is shown to be homologous to archaeal protein L10e, cluster to the same region as the rRNA mutations. The Evn-binding site overlaps with the binding site of initiation factor 2. Evn inhibits activity of initiation factor 2 in vitro, suggesting that the drug interferes with formation of the 70S initiation complex. The site of Evn binding and its mode of action are distinct from other ribosome-targeted antibiotics. This antibiotic target site can potentially be used for the development of new antibacterial drugs.

The complete genome of the crenarchaeon Sulfolobus solfataricus P2

The genome of the crenarchaeon Sulfolobus solfataricus P2 contains 2,992,245 bp on a single chromosome and encodes 2,977 proteins and many RNAs. One-third of the encoded proteins have no detectable homologs in other sequenced genomes. Moreover, 40% appear to be archaeal-specific, and only 12% and 2.3% are shared exclusively with bacteria and eukarya, respectively. The genome shows a high level of plasticity with 200 diverse insertion sequence elements, many putative nonautonomous mobile elements, and evidence of integrase-mediated insertion events. There are also long clusters of regularly spaced tandem repeats. Different transfer systems are used for the uptake of inorganic and organic solutes, and a wealth of intracellular and extracellular proteases, sugar, and sulfur metabolizing enzymes are encoded, as well as enzymes of the central metabolic pathways and motility proteins. The major metabolic electron carrier is not NADH as in bacteria and eukarya but probably ferredoxin. The essential components required for DNA replication, DNA repair and recombination, the cell cycle, transcriptional initiation and translation, but not DNA folding, show a strong eukaryal character with many archaeal-specific features. The results illustrate major differences between crenarchaea and euryarchaea, especially for their DNA replication mechanism and cell cycle processes and their translational apparatus.

Determinants of RNA polymerase alpha subunit for interaction with beta, beta', and sigma subunits: hydroxyl-radical protein footprinting.

Escherichia coli RNA polymerase (RNAP) alpha subunit serves as the initiator for RNAP assembly, which proceeds according to the pathway 2 alpha-->alpha 2-->alpha 2 beta-->alpha 2 beta beta'-->alpha 2 beta beta' sigma. In this work, we have used hydroxyl-radical protein footprinting to define determinants of alpha for interaction with beta, beta', and sigma. Our results indicate that amino acids 30-75 of alpha are protected from hydroxyl-radical-mediated proteolysis upon interaction with beta (i.e., in alpha 2 beta, alpha 2 beta beta', and alpha 2 beta beta' sigma), and amino acids 175-210 of alpha are protected from hydroxyl-radical-mediated proteolysis upon interaction with beta' (i.e., in alpha 2 beta beta' and alpha 2 beta beta' sigma). The protected regions are conserved in the alpha homologs of prokaryotic, eukaryotic, archaeal, and chloroplast RNAPs and contain sites of substitutions that affect RNAP assembly. We conclude that the protected regions define determinants of alpha for direct functional interaction with beta and beta'. The observed maximal magnitude of protection upon interaction with beta and the observed maximal magnitude of protection upon interaction with beta' both correspond to the expected value for complete protection of one of the two alpha protomers of RNAP (i.e., 50% protection). We propose that only one of the two alpha protomers of RNAP interacts with beta and that only one of the two alpha protomers of RNAP interacts with beta'.

A minimal gene set for cellular life derived by comparison of complete bacterial genomes.

The recently sequenced genome of the parasitic bacterium Mycoplasma genitalium contains only 468 identified protein-coding genes that have been dubbed a minimal gene complement [Fraser, C.M., Gocayne, J.D., White, O., Adams, M.D., Clayton, R.A., et al. (1995) Science 270, 397-403]. Although the M. genitalium gene complement is indeed the smallest among known cellular life forms, there is no evidence that it is the minimal self-sufficient gene set. To derive such a set, we compared the 468 predicted M. genitalium protein sequences with the 1703 protein sequences encoded by the other completely sequenced small bacterial genome, that of Haemophilus influenzae. M. genitalium and H. influenzae belong to two ancient bacterial lineages, i.e., Gram-positive and Gram-negative bacteria, respectively. Therefore, the genes that are conserved in these two bacteria are almost certainly essential for cellular function. It is this category of genes that is most likely to approximate the minimal gene set. We found that 240 M. genitalium genes have orthologs among the genes of H. influenzae. This collection of genes falls short of comprising the minimal set as some enzymes responsible for intermediate steps in essential pathways are missing. The apparent reason for this is the phenomenon that we call nonorthologous gene displacement when the same function is fulfilled by nonorthologous proteins in two organisms. We identified 22 nonorthologous displacements and supplemented the set of orthologs with the respective M. genitalium genes. After examining the resulting list of 262 genes for possible functional redundancy and for the presence of apparently parasite-specific genes, 6 genes were removed. We suggest that the remaining 256 genes are close to the minimal gene set that is necessary and sufficient to sustain the existence of a modern-type cell. Most of the proteins encoded by the genes from the minimal set have eukaryotic or archaeal homologs but seven key proteins of DNA replication do not. We speculate that the last common ancestor of the three primary kingdoms had an RNA genome. Possibilities are explored to further reduce the minimal set to model a primitive cell that might have existed at a very early stage of life evolution.

The root of the universal tree and the origin of eukaryotes based on elongation factor phylogeny.

The genes for the protein synthesis elongation factors Tu (EF-Tu) and G (EF-G) are the products of an ancient gene duplication, which appears to predate the divergence of all extant organismal lineages. Thus, it should be possible to root a universal phylogeny based on either protein using the second protein as an outgroup. This approach was originally taken independently with two separate gene duplication pairs, (i) the regulatory and catalytic subunits of the proton ATPases and (ii) the protein synthesis elongation factors EF-Tu and EF-G. Questions about the orthology of the ATPase genes have obscured the former results, and the elongation factor data have been criticized for inadequate taxonomic representation and alignment errors. We have expanded the latter analysis using a broad representation of taxa from all three domains of life. All phylogenetic methods used strongly place the root of the universal tree between two highly distinct groups, the archaeons/eukaryotes and the eubacteria. We also find that a combined data set of EF-Tu and EF-G sequences favors placement of the eukaryotes within the Archaea, as the sister group to the Crenarchaeota. This relationship is supported by bootstrap values of 60-89% with various distance and maximum likelihood methods, while unweighted parsimony gives 58% support for archaeal monophyly.