Complete genome sequence of Methanobacterium thermoautotrophicum deltaH: functional analysis and comparative genomics

Adenylosuccinate synthetase genes: molecular cloning and phylogenetic analysis of a highly conserved archaeal gene

Adenylosuccinate synthetase (PurA) catalyzes the first step in the de novo AMP synthesis and has been extensively studied in both Bacteria and Eukarya. We cloned the purA gene from the hyperthermophilic archaeon, Pyrococcus furiosus. The gene appears to be individually transcribed and encodes a protein of 339 amino acids. The amino acid sequence comparison with other archael PurAs found from recent genome analyses indicated that two deletions, one central and the other C-terminal, are a common feature of archaeal PurAs. None of the 21 PurA homologues analyzed from Eukarya and Bacteria exhibited this feature. Amino acid sequences of PurAs in Archaea showed 64% average identities which were significantly higher than the 50% and 55% calculated for Bacteria and Eukarya, respectively. Several residues conserved in PurAs of both Eukarya and Bacteria and shown to be of catalytic importance are missing in the archaeal PurAs. Phylogenetic analysis using PurA as the marker grouped life into 3 domains, hence it was consistent with results derived from 16-18S ribosomal RNA sequences. The topology within the three domains, in general, portrayed the hitherto accepted evolutionary relationship among the organisms utilized. PurA can, thus, serve as an additional marker to evaluate phylogenetic inferences drawn from sequence data from rRNA and other conserved genes. The presence of two unique deletions in both euryarchaeal and crenarchaeal PurAs, but not in those of Bacteria and Eukarya, is a strong evidence confirming the common lineage of these two subdomains of Archaea.

The euryarchaeotes, a subdomain of Archaea, survive on a single DNA polymerase: fact or farce?

Archaea is now recognized as the third domain of life. Since their discovery, much effort has been directed towards understanding the molecular biology and biochemistry of Archaea. The objective is to comprehend the complete structure and the depth of the phylogenetic tree of life. DNA replication is one of the most important events in living organisms and DNA polymerase is the key enzyme in the molecular machinery which drives the process. All archaeal DNA polymerases were thought to belong to family B. This was because all of the products of pol genes that had been cloned showed amino acid sequence similarities to those of this family, which includes three eukaryal DNA replicases and Escherichia coli DNA polymerase II. Recently, we found a new heterodimeric DNA polymerase from the hyperthermophilic archaeon, Pyrococcus furiosus. The genes coding for the subunits of this DNA polymerase are conserved in the euryarchaeotes whose genomes have been completely sequenced. The biochemical characteristics of the novel DNA polymerase family suggest that its members play an important role in DNA replication within euryarchaeal cells. We review here our current knowledge on DNA polymerases in Archaea with emphasis on the novel DNA polymerase discovered in Euryarchaeota.

A heterodimeric DNA polymerase: evidence that members of Euryarchaeota possess a distinct DNA polymerase

We describe here a DNA polymerase family highly conserved in Euryarchaeota, a subdomain of Archaea. The DNA polymerase is composed of two proteins, DP1 and DP2. Sequence analysis showed that considerable similarity exists between DP1 and the second subunit of eukaryotic DNA polymerase delta, a protein essential for the propagation of Eukarya, and that DP2 has conserved motifs found in proteins with nucleotide-polymerizing activity. These results, together with our previous biochemical analyses of one of the members, DNA polymerase II (DP1 + DP2) from Pyrococcus furiosus, implicate the DNA polymerases of this family in the DNA replication process of Euryarchaeota. The discovery of this DNA-polymerase family, aside from providing an opportunity to enhance our knowledge of the evolution of DNA polymerases, is a significant step toward the complete understanding of DNA replication across the three domains of life.

FKBP-type peptidyl-prolyl cis-trans isomerase from a sulfur-dependent hyperthermophilic archaeon, Thermococcus sp. KS-1

A gene encoding an FK506 binding protein (FKBP)-type peptidyl-prolyl cis-trans isomerase (PPIase) was cloned from a hyperthermophilic archaeon, Thermococcus sp. KS-1, and sequenced. This gene encoded an FKBP with 159 amino-acid residues with a molecular mass of 17.6kDa. Two insertion sequences with 13 and 44 amino acids were found in the regions corresponding to the bulge and flap regions of human FKBP-12, respectively. Comparison with other archaeal FKBP sequences obtained from reported genome sequences revealed that the insertion sequences in the bulge and flap regions were common to archaeal FKBPs. It was also revealed that archaeal FKBPs are classified into two groups: one is approx. 17kDa and the other 27kDa. This Thermococcus FKBP (TcFK) belonged to the smaller archaeal FKBP. In this TcFK, 9 out of 15 amino acid residues forming the FK506 binding pocket of human FKBP12 were found. This gene was expressed in Escherichia coli and the recombinant protein was purified. The purified protein showed PPIase activity and its activity was inhibited by FK506 with an IC50 of 7 microM. This enzyme showed high kinetic stability with a half-life of 40 min at 100 degrees C. Catalytic efficiency of this recombinant PPIase was 1.2-times higher with the substrate N-succinyl-A-L-P-F-p-nitroanilide than with N-succinyl-A-A-P-F-p-nitroanilide.

Human asparaginyl-tRNA synthetase: molecular cloning and the inference of the evolutionary history of Asx-tRNA synthetase family

We have cloned and sequenced a cDNA encoding human cytoplasmic asparaginyl-tRNA synthetase (AsnRS). The N-terminal appended domain of 112 amino acid represents the signature sequence for the eukaryotic AsnRS and is absent from archaebacterial or eubacterial enzymes. The canonical ortholog for AsnRS is absent from most archaebacterial and some eubacterial genomes, indicating that in those organisms, formation of asparaginyl-tRNA is independent of the enzyme. The high degree of sequence conservation among asparaginyl- and aspartyl-tRNA synthetases (AsxRS) made it possible to infer the evolutionary paths of the two enzymes. The data show the neighbor relationship between AsnRS and eubacterial aspartyl-tRNA synthetase, and support the occurrence of AsnRS early in the course of evolution, which is in contrast to the proposed late occurrence of glutaminyl-tRNA synthetase.

Aminoacylation of tRNA in the evolution of an aminoacyl-tRNA synthetase

Aminoacyl-tRNA synthetases catalyze aminoacylation of tRNAs by joining an amino acid to its cognate tRNA. The selection of the cognate tRNA is jointly determined by separate structural domains that examine different regions of the tRNA. The cysteine-tRNA synthetase of Escherichia coli has domains that select for tRNAs containing U73, the GCA anticodon, and a specific tertiary structure at the corner of the tRNA L shape. The E. coli enzyme does not efficiently recognize the yeast or human tRNACys, indicating the evolution of determinants for tRNA aminoacylation from E. coli to yeast to human and the coevolution of synthetase domains that interact with these determinants. By successively modifying the yeast and human tRNACys to ones that are efficiently aminoacylated by the E. coli enzyme, we have identified determinants of the tRNA that are important for aminoacylation but that have diverged in the course of evolution. These determinants provide clues to the divergence of synthetase domains. We propose that the domain for selecting U73 is conserved in evolution. In contrast, we propose that the domain for selecting the corner of the tRNA L shape diverged early, after the separation between E. coli and yeast, while that for selecting the GCA-containing anticodon loop diverged late, after the separation between yeast and human.

The CCA-adding enzyme has a single active site

The CCA-adding enzyme (tRNA nucleotidyltransferase) synthesizes and repairs the 3'-terminal CCA sequence of tRNA. The eubacterial, eukaryotic, and archaeal CCA-adding enzymes all share a single active-site signature motif, which identifies these enzymes as belonging to the nucleotidyltransferase superfamily. Here we show that mutations at Asp-53 or Asp-55 of the Sulfolobus shibatae signature sequence abolish addition of both C and A, demonstrating that a single active site is responsible for addition of both nucleotides. Mutations at Asp-106 (and to a lesser extent, at Glu-173 and Asp-215) selectively impaired addition of A, but not C. We have previously demonstrated that the tRNA acceptor stem remains fixed on the surface of the CCA-adding enzyme during C and A addition (Shi, P.-Y., Maizels, N., and Weiner, A. M. (1998) EMBO J. 17, 3197-3206). Taken together with this new evidence that there is a single active site for catalysis, our data suggest that specificity of nucleotide addition is determined by a process of collaborative templating: as the single active site catalyzes addition of each nucleotide, the growing 3'-end of the tRNA would progressively refold to create a binding pocket for addition of the next nucleotide.

Predicting function: from genes to genomes and back

Predicting function from sequence using computational tools is a highly complicated procedure that is generally done for each gene individually. This review focuses on the added value that is provided by completely sequenced genomes in function prediction. Various levels of sequence annotation and function prediction are discussed, ranging from genomic sequence to that of complex cellular processes. Protein function is currently best described in the context of molecular interactions. In the near future it will be possible to predict protein function in the context of higher order processes such as the regulation of gene expression, metabolic pathways and signalling cascades. The analysis of such higher levels of function description uses, besides the information from completely sequenced genomes, also the additional information from proteomics and expression data. The final goal will be to elucidate the mapping between genotype and phenotype.

Isolation of a minD-like gene in the hyperthermophilic archaeon Pyrococcus AL585, and phylogenetic characterization of related proteins in the three domains of life

The region upstream of the dinF-like gene of the hyperthermophilic archaeon Pyrococcus strain AL585 has been cloned and sequenced. This region contains an open reading frame (ORF) that encodes a polypeptide with a high similarity to MinD proteins and their Mrp paralogues. Transcripts of the dinF-like and the minD-like genes were detected by RT-PCR, indicating that they are both expressed in vivo. The MinD and MinD-like proteins belong to a broad family of ATPases involved in chromosome and plasmid partitioning. MinD-like proteins can be defined by specific amino-acid sequence signatures. A systematic search for proteins sharing these signatures in current databases and newly sequenced genomes show that MinD-like proteins are present in all archaeal genomes sequenced so far, often in several copies. Phylogenetic analysis identifies two groups of MinD-like proteins which are also characterized by more conserved amino-acid motifs. A first group, which includes the Escherichia coli MinD and the Pyrococcus AL585 MinDL protein, contains only procaryotic proteins. This group can be further divided into a subgroup of archaeal proteins and two subgroups of bacterial proteins. A second group includes proteins more related to the E. coli Mrp protein and contains representants of the three domains of life. The conservation of MinD-like proteins in the three domains of life suggests that these proteins play a central role in cellular metabolism.

Isolation of carbon- and nitrogen-deprivation-induced loci of Sinorhizobium meliloti 1021 by Tn5-luxAB mutagenesis

Soil bacteria, such as Sinorhizobium meliloti, are subject to variation in environmental conditions, including carbon- and nitrogen-deprivation. The ability of bacteria to sense changes in their environment and respond accordingly is of vital importance to their survival and persistence in the soil and rhizosphere. A derivative of Tn5 which creates transcriptional fusions to the promoterless luxAB genes was used to mutagenize S. meliloti 1021 and 5000 insertion mutants were subsequently screened for gene fusions induced by selected environmental stresses. The isolation of 21 gene fusions induced by nitrogen-deprivation and 12 induced by carbon-deprivation is described. Cloning and partial DNA sequence analysis of the transposon-tagged loci revealed a variety of novel genes, as well as S. meliloti genes with significant similarity to known bacterial loci. In addition, nodule occupancy studies were carried out with selected Tn5-luxAB insertion mutants to examine the role of the tagged genes in competition.

Snapshot of a large dynamic replicon in a halophilic archaeon: megaplasmid or minichromosome?

Extremely halophilic archaea, which flourish in hypersaline environments, are known to contain a variety of large dynamic replicons. Previously, the analysis of one such replicon, pNRC100, in Halobacterium sp. strain NRC-1, showed that it undergoes high-frequency insertion sequence (IS) element-mediated insertions and deletions, as well as inversions via recombination between 39-kb-long inverted repeats (IRs). Now, the complete sequencing of pNRC100, a 191,346-bp circle, has shown the presence of 27 IS elements representing eight families. A total of 176 ORFs or likely genes of 850-bp average size were found, 39 of which were repeated within the large IRs. More than one-half of the ORFs are likely to represent novel genes that have no known homologs in the databases. Among ORFs with previously characterized homologs, three different copies of putative plasmid replication and four copies of partitioning genes were found, suggesting that pNRC100 evolved from IS element-mediated fusions of several smaller plasmids. Consistent with this idea, putative genes typically found on plasmids, including those encoding a restriction-modification system and arsenic resistance, as well as buoyant gas-filled vesicles and a two-component regulatory system, were found on pNRC100. However, additional putative genes not expected on an extrachromosomal element, such as those encoding an electron transport chain cytochrome d oxidase, DNA nucleotide synthesis enzymes thioredoxin and thioredoxin reductase, and eukaryotic-like TATA-binding protein transcription factors and a chromosomal replication initiator protein were also found. A multi-step IS element-mediated process is proposed to account for the acquisition of these chromosomal genes. The finding of essential genes on pNRC100 and its property of resistance to curing suggest that this replicon may be evolving into a new chromosome.

How to interpret an anonymous bacterial genome: machine learning approach to gene identification

In this report we address the problem of accurate statistical modeling of DNA sequences, either coding or noncoding, for a bacterial species whose genome (or a large portion) was sequenced but not yet characterized experimentally. Availability of these models is critical for successful solution of the genome annotation task by statistical methods of gene finding. We present the method, GeneMark-Genesis, which learns the parameters of Markov models of protein-coding and noncoding regions from anonymous bacterial genomic sequence. These models are subsequently used in the GeneMark and GeneMark.hmm gene-finding programs. Although there is basically one model of a noncoding region for a given genome, several models of protein-coding region are automatically obtained by GeneMark-Genesis. The diversity of protein-coding models reflects the diversity of oligonucleotide compositions, particularly the diversity of codon usage strategies observed in genes from one and the same genome. In the simplest and the most important case, there are just two gene models-typical and atypical ones. We show that the atypical model allows one to predict genes that escape identification by the typical model. Many genes predicted by the atypical model appear to be horizontally transferred genes. The early versions of GeneMark-Genesis were used for annotating the genomes of Methanoccocus jannaschii and Helicobacter pylori. We report the results of accuracy testing of the full-scale version of GeneMark-Genesis on 10 completely sequenced bacterial genomes. Interestingly, the GeneMark.hmm program that employed the typical and atypical models defined by GeneMark-Genesis was able to predict 683 new atypical genes with 176 of them confirmed by similarity search.

Identification of a gene for a rubrerythrin/nigerythrin-like protein in Spirillum volutans by using amino acid sequence data from mass spectrometry and NH2-terminal sequencing

A hydrogen peroxide-resistant mutant of the catalase-negative microaerophile, Spirillum volutans, constitutively expresses a 21.5 kDa protein that is undetectable and non-inducible in the wild-type cells. Part of the gene that encodes the protein was cloned using amino acid sequence data obtained by both mass spectrometry and NH2-terminal sequencing. The deduced 158 amino acid polypeptide shows high relatedness to rubrerythrin and nigerythrin previously described in the anaerobes Clostridium perfringens and Desulfovibrio vulgaris. The protein also shows high similarity to putative rubrerythrin proteins found in the anaerobic archeons Archaeoglobus fulgidus, Methanococcus jannaschii and Methanobacterium thermoautotrophicum. This is the first report of this type of protein in an organism that must respire with oxygen. It seems likely that the novel combination of methodologies used in this study could be applied to the rapid cloning of other genes in bacteria for which no genomic library yet exists.

Nucleotide excision repair in the third kingdom

Nucleotide excision repair, a general repair mechanism for removing DNA damage, is initiated by dual incisions bracketing the lesion. In procaryotes, the dual incisions result in excision of the damage in 12- to 13-nucleotide-long oligomers, and in eucaryotes they result in excision of the damage in the form of 24- to 32-nucleotide-long oligomers. We wished to find out if Archaea perform excision repair. Using cell extracts from Methanobacterium thermoautotrophicum, we found that this organism removes UV-induced (6-4) photoproducts in the form of 10- to 11-mers by incising the sixth to seventh phosphodiester bond 5' to the damage and the fourth phosphodiester bond 3' to the damage.

Cloning, sequencing, and expression of the gene encoding cyclic 2, 3-diphosphoglycerate synthetase, the key enzyme of cyclic 2, 3-diphosphoglycerate metabolism in Methanothermus fervidus

Cyclic 2,3-diphosphoglycerate synthetase (cDPGS) catalyzes the synthesis of cyclic 2,3-diphosphoglycerate (cDPG) by formation of an intramolecular phosphoanhydride bond in 2,3-diphosphoglycerate. cDPG is known to be accumulated to high intracellular concentrations (>300 mM) as a putative thermoadapter in some hyperthermophilic methanogens. For the first time, we have purified active cDPGS from a methanogen, the hyperthermophilic archaeon Methanothermus fervidus, sequenced the coding gene, and expressed it in Escherichia coli. cDPGS purification resulted in enzyme preparations containing two isoforms differing in their electrophoretic mobility under denaturing conditions. Since both polypeptides showed the same N-terminal amino acid sequence and Southern analyses indicate the presence of only one gene coding for cDPGS in M. fervidus, the two polypeptides originate from the same gene but differ by a not yet identified modification. The native cDPGS represents a dimer with an apparent molecular mass of 112 kDa and catalyzes the reversible formation of the intramolecular phosphoanhydride bond at the expense of ATP. The enzyme shows a clear preference for the synthetic reaction: the substrate affinity and the Vmax of the synthetic reaction are a factor of 8 to 10 higher than the corresponding values for the reverse reaction. Comparison with the kinetic properties of the electrophoretically homogeneous, apparently unmodified recombinant enzyme from E. coli revealed a twofold-higher Vmax of the enzyme from M. fervidus in the synthesizing direction.

Glutamyl-tRNA(Gln) amidotransferase in Deinococcus radiodurans may be confined to asparagine biosynthesis

Asparaginyl-tRNA (Asn-tRNA) and glutaminyl-tRNA (Gln-tRNA) are essential components of protein synthesis. They can be formed by direct acylation by asparaginyl-tRNA synthetase (AsnRS) or glutaminyl-tRNA synthetase (GlnRS). The alternative route involves transamidation of incorrectly charged tRNA. Examination of the preliminary genomic sequence of the radiation-resistant bacterium Deinococcus radiodurans suggests the presence of both direct and indirect routes of Asn-tRNA and Gln-tRNA formation. Biochemical experiments demonstrate the presence of AsnRS and GlnRS, as well as glutamyl-tRNA synthetase (GluRS), a discriminating and a nondiscriminating aspartyl-tRNA synthetase (AspRS). Moreover, both Gln-tRNA and Asn-tRNA transamidation activities are present. Surprisingly, they are catalyzed by a single enzyme encoded by three ORFs orthologous to Bacillus subtilis gatCAB. However, the transamidation route to Gln-tRNA formation is idled by the inability of the discriminating D. radiodurans GluRS to produce the required mischarged Glu-tRNAGln substrate. The presence of apparently redundant complete routes to Asn-tRNA formation, combined with the absence from the D. radiodurans genome of genes encoding tRNA-independent asparagine synthetase and the lack of this enzyme in D. radiodurans extracts, suggests that the gatCAB genes may be responsible for biosynthesis of asparagine in this asparagine prototroph.

PCR-based subtractive hybridization and differences in gene content among strains of Helicobacter pylori

Genes that are characteristic of only certain strains of a bacterial species can be of great biologic interest. Here we describe a PCR-based subtractive hybridization method for efficiently detecting such DNAs and apply it to the gastric pathogen Helicobacter pylori. Eighteen DNAs specific to a monkey-colonizing strain (J166) were obtained by subtractive hybridization against an unrelated strain whose genome has been fully sequenced (26695). Seven J166-specific clones had no DNA sequence match to the 26695 genome, and 11 other clones were mixed, with adjacent patches that did and did not match any sequences in 26695. At the protein level, seven clones had homology to putative DNA restriction-modification enzymes, and two had homology to putative metabolic enzymes. Nine others had no database match with proteins of assigned function. PCR tests of 13 unrelated H. pylori strains by using primers specific for 12 subtracted clones and complementary Southern blot hybridizations indicated that these DNAs are highly polymorphic in the H. pylori population, with each strain yielding a different pattern of gene-specific PCR amplification. The search for polymorphic DNAs, as described here, should help identify previously unknown virulence genes in pathogens and provide new insights into microbial genetic diversity and evolution.

Chlorella virus PBCV-1 encodes functional glutamine: fructose-6-phosphate amidotransferase and UDP-glucose dehydrogenase enzymes

DNA sequence analysis of the 330-kb Chlorella virus PBCV-1 genome unexpectedly revealed several open reading frames which encode proteins that are homologous to sugar-manipulating enzymes including glutamine:fructose-6-phosphate amidotransferase (GFAT), UDP-glucose dehydrogenase (UDP-GlcDH), and hyaluronan synthase (HAS). PBCV-1 genes encoding the putative GFAT and UDP-GlcDH enzymes were expressed in Escherichia coli, and both recombinant proteins have the predicted enzyme activity in cell free extracts. These same two genes are transcribed early in PBCV-1 infection, and both genes are widespread among the Chlorella viruses. The products of the reactions catalyzed by these two enzymes are precursors in the biosynthesis of hyaluronan polysaccharide. Previous experiments established that, like the GFAT and UDP-GlcDH genes, the HAS gene is transcribed early and encodes a functional enzyme (DeAngelis, P. L., Jing. W., Graves, M. V., Burbank, D. E., and Van Etten, J. L. (1997) Science 278, 1800-1803). Interestingly, the predicted amino-acid sequences of the PBCV-1 GFAT and UDP-GlcDH enzymes are more similar to bacterial GFAT and UDP-GlcDH enzymes than to their eukaryotic counterparts. In contrast, the amino-acid sequence of the PBCV-1 HAS enzyme more closely resembles eukaryotic enzymes.

Group II chaperonin in a thermophilic methanogen, Methanococcus thermolithotrophicus. Chaperone activity and filament-forming ability

A gene encoding 544 amino acids for a subunit of group II chaperonin (thermosome) was cloned from a thermophilic methanogen, Methanococcus thermolithotrophicus. The deduced amino acid sequence showed 66.5, 56.1, and 20.1% similarities to those of Methanopyrus kandleri and Thermoplasma acidophilum and group I chaperonin of Escherichia coli, respectively. We call this chaperonin MTTS (M. thermolithotrophicus thermosome). The MTTS gene was expressed in E. coli. The purified recombinant MTTS seemed to be monomeric on gel filtration in the absence of Mg2+ and ATP. The monomer assembled to an oligomer (complex) in the presence of 50 mM MgCl2, 0.25 mM ATP, and 0.3 M (NH4)2SO4. It was eluted immediately before the elution volume of E. coli GroEL tetradecamer on gel filtration with a TSKgel G3000SWXL column. This reconstructed MTTS complex showed the cylindrical structure with two stacked rings in electron microscopy. The MTTS complex formed filamentous structures in the presence of Mg2+ and ATP at the protein concentration above 3.0 mg/ml. This filament formation was reversible. The MTTS filament was dissociated to the complex by dilution to the protein concentration of 0.2 mg/ml, even in the presence of Mg2+ and ATP. The MTTS complex exhibited weak ATPase activity with the hydrolysis rate of 74 mol of ATP hydrolysis/mol of MTTS complex/min at 70 degreesC. The MTTS complex promoted the refolding of chemically denatured thermophilic archaeal citrate synthase and glucose dehydrogenase at 50 degreesC in an ATP-dependent fashion. The analysis of nucleotide specificity of chaperone activity of MTTS suggested that it was coupled with hydrolysis of ATP, CTP, or UTP.

Constructing multigenome views of whole microbial genomes

We have designed and implemented a system to carry out cross-genome comparisons of open reading frames (ORFs) from multiple genomes. This implementation includes a genome profiling system that allows us to explore pairwise comparisons at different levels of match similarity and ask biologically motivated queries involving number and identity of ORFs, their function, functional category, distribution in genomes or in biological domains, and statistics on their matches and match families. This analysis required precise definition of new classification terms and concepts. We define the terms genomic signature, summary signature, biologic domain signature, domain class, match level, match family, and extended match family, then use these terms to define concepts, including genomically universal proteins and proteins characteristics of sets of genomes. We initiate an analysis based on automated FASTA (Pearson, 1996) comparison of 22,419 conceptually translated protein sequences from nine microbial genomes.

The 1.8 A crystal structure of the ycaC gene product from Escherichia coli reveals an octameric hydrolase of unknown specificity

BACKGROUND

The ycaC gene comprises a 621 base pair open reading frame in Escherichia coli. The ycaC gene product (ycaCgp) is uncharacterized and has no assigned function. The closest sequence homologs with an assigned function belong to a family of bacterial hydrolases that catalyze isochorismatase-like reactions, but these have only low sequence similarity to ycaCgp (approximately 20% amino acid identity). The ycaCgp was obtained and identified during crystallization trials of an unrelated E. coli protein with which it co-purified.

RESULTS

The 1.8 A crystal structure of ycaCgp reveals an octameric complex comprised of two tetrameric rings. A large three-layer (alphabetaalpha) sandwich domain and a small helical domain form the folded structure of the monomeric unit. Comparisons with sequence and structure databases suggest that ycaCgp belongs to a diverse family of bacterial hydrolases. The most closely related three-dimensional structure is that of the D2 tetrameric N-carbamoylsarcosine amidohydrolase (CSHase) from an Arthrobacter species. A conspicuous cleft between two ycaCgp subunits contains several conserved residues including Cys118, which we propose to be catalytic. In the active site, a nonprolyl cis peptide bond precedes Val114 and coincides with a cis peptide bond in CSHase in a region of dissimilar sequence. The crystal structure reveals a probable error or mutation relative to the reported genomic sequence.

CONCLUSIONS

Although the specific function of ycaCgp is not yet known, structural studies solidify the relationship of this protein to other hydrolases and illuminate its active site and key elements of the catalytic mechanism.

Electrochemical and spectroscopic properties of the iron-sulfur flavoprotein from Methanosarcina thermophila

An iron-sulfur flavoprotein (Isf) from the methanoarchaeaon Methanosarcina thermophila, which participates in electron transfer reactions required for the fermentation of acetate to methane, was characterized by electrochemistry and EPR and Mössbauer spectroscopy. The midpoint potential (Em) of the FMN/FMNH2 couple was -0.277 V. No flavin semiquinone was observed during potentiometric titrations; however, low amounts of the radical were observed when Isf was quickly frozen after reaction with CO and the CO dehydrogenase/acetyl-CoA synthase complex from M. thermophila. Isf contained a [4Fe-4S]2+/1+ cluster with g values of 2.06 and 1.93 and an unusual split signal with g values at 1.86 and 1.82. The unusual morphology was attributed to microheterogeneity among Isf molecules. The Em value for the 2+/1+ redox couple of the cluster was -0.394 V. Extracts from H2-CO2-grown Methanobacterium thermoautotrophicum cells catalyzed either the H2- or CO-dependent reduction of M. thermophila Isf. In addition, Isf homologs were found in the genomic sequences of the CO2-reducing methanoarchaea M. thermoautotrophicum and Methanococcus jannaschii. These results support a general role for Isf in electron transfer reactions of both acetate-fermenting and CO2-reducing methanoarchaea. It is suggested that Isf functions to couple electron transfer from ferredoxin to membrane-bound electron carriers, such as methanophenazine and/or b-type cytochromes.

A comparative analysis of ABC transporters in complete microbial genomes

The ABC transporter is a major class of cellular translocation machinery in all bacterial species encoded in the largest set of paralogous genes. The operon structure is frequently found for the genes of three molecular components: the ATP-binding protein, the membrane protein, and the substrate-binding protein. Here, we developed an "ortholog group table" by comparison and classification of known and putative ABC transporters in the complete genomes of seven microorganisms. Our procedure was to first search and classify the most conserved ATP-binding protein components by the sequence similarity and then to classify the entire transporter units by examining the similarity of the other components and the conservation of the operon structure. The resulting 25 ortholog groups of ABC transporters were well correlated with known functions. Through the analysis, we could assign substrate specificity to hypothetical transporters, predict additional transporter operons, and identify novel types of putative transporters. The ortholog group table was also used as a reference data set for functional assignment in four additional genomes. In general, the ABC transporter operons were strongly conserved despite the extensive shuffling of gene locations in bacterial evolution. In Synechocystis, however, the tendency of forming operons was clearly diminished. Our result suggests that the ancestral ABC transporter operons may have arisen early in evolution before the speciation of bacteria and archaea.

The serine, threonine, and/or tyrosine-specific protein kinases and protein phosphatases of prokaryotic organisms: a family portrait

Inspection of the genomes for the bacteria Bacillus subtilis 168, Borrelia burgdorferi B31, Escherichia coli K-12, Haemophilus influenzae KW20, Helicobacter pylori 26695, Mycoplasma genitalium G-37, and Synechocystis sp PCC 6803 and for the archaeons Archaeoglobus fulgidus VC-16 DSM4304, Methanobacterium thermoautotrophicum delta H, and Methanococcus jannaschii DSM2661 revealed that each contains at least one ORF whose predicted product displays sequence features characteristic of eukaryote-like protein-serine/threonine/tyrosine kinases and protein-serine/threonine/tyrosine phosphatases. Orthologs for all four major protein phosphatase families (PPP, PPM, conventional PTP, and low molecular weight PTP) were present in the bacteria surveyed, but not all strains contained all types. The three archaeons surveyed lacked recognizable homologs of the PPM family of eukaryotic protein-serine/threonine phosphatases; and only two prokaryotes were found to contain ORFs for potential phosphatases from all four major families. Intriguingly, our searches revealed a potential ancestral link between the catalytic subunits of microbial arsenate reductases and the protein-tyrosine phosphatases; they share similar ligands (arsenate versus phosphate) and features of their catalytic mechanism (formation of arseno-versus phospho-cysteinyl intermediates). It appears that all prokaryotic organisms, at one time, contained the genetic information necessary to construct protein phosphorylation-dephosphorylation networks that target serine, threonine, and/or tyrosine residues on proteins. However, the potential for functional redundancy among the four protein phosphatase families has led many prokaryotic organisms to discard one, two, or three of the four.

Genomic analysis reveals chromosomal variation in natural populations of the uncultured psychrophilic archaeon Cenarchaeum symbiosum

Molecular phylogenetic surveys have recently revealed an ecologically widespread crenarchaeal group that inhabits cold and temperate terrestrial and marine environments. To date these organisms have resisted isolation in pure culture, and so their phenotypic and genotypic characteristics remain largely unknown. To characterize these archaea, and to extend methodological approaches for characterizing uncultivated microorganisms, we initiated genomic analyses of the nonthermophilic crenarchaeote Cenarchaeum symbiosum found living in association with a marine sponge, Axinella mexicana. Complex DNA libraries derived from the host-symbiont population yielded several large clones containing the ribosomal operon from C. symbiosum. Unexpectedly, cloning and sequence analysis revealed the presence of two closely related variants that were consistently found in the majority of host individuals analyzed. Homologous regions from the two variants were sequenced and compared in detail. The variants exhibit >99.2% sequence identity in both small- and large-subunit rRNA genes and they contain homologous protein-encoding genes in identical order and orientation over a 28-kbp overlapping region. Our study not only indicates the potential for characterizing uncultivated prokaryotes by genome sequencing but also identifies the primary complication inherent in the approach: the widespread genomic microheterogeneity in naturally occurring prokaryotic populations.

Life's third domain (Archaea): an established fact or an endangered paradigm?

The three-domain proposal of Woese et al. (Proc. Natl. Acad. Sci. USA 87, 4576 (1990)) divides all living organisms into three primary groups or domains named Archaea (or archaebacteria), Bacteria (or eubacteria), and Eucarya (or eukaryotes), with Eucarya being relatives (or descendants) of Archaea. Although this proposal is currently widely accepted, sequence features and phylogenies derived from many highly conserved proteins are inconsistent with it and point to a close and specific relationship between archaebacteria and gram-positive bacteria, whereas gram-negative bacteria are indicated to be phylogenetically distinct. A closer relationship of archaebacteria to gram-positive bacteria in comparison to gram-negative bacteria is generally seen for the majority of the available gene/protein sequences. To account for these results, and the fact that both archaebacteria and gram-positive bacteria are prokaryotes surrounded by a single cell membrane, I propose that the primary division within prokaryotes is between Monoderm prokaryotes (surrounded by a single membrane) and Diderm prokaryotes (i.e., all true gram-negative bacteria containing both an inner cytoplasmic membrane and an outer membrane). This proposal is consistent with both cell morphology and signature sequences in different proteins. Protein phylogenies and signature sequences also show that all eukaryotic cells have received significant gene contributions from both an archaebacterium and a gram- negative eubacterium. Thus, the hypothesis that archaebacteria and eukaryotes shared a common ancestor exclusive of eubacteria, or that the ancestral eukaryotic cell directly descended from an archaea, is erroneous. These results call into question the validity of the currently popular three-domain proposal and the assignment of a domain status to archaebacteria. A new classifica- tion of organisms consistent with phenotype and macromolecular sequence data is proposed.

Molecular analysis of Methanobacterium phage psiM2

The methanogenic archaeon Methanobacterium thermoautotrophicum Marburg is infected by the double-stranded DNA phage psiM2. The complete phage genome sequence of 26 111 bp was established. Thirty-one open reading frames (orfs), all of them organized in the same direction of transcription, were identified. On the basis of comparison of the deduced amino acid sequences to known proteins and by searching for conserved motifs, putative functions were assigned to the products of six orfs. These included three proteins involved in packaging DNA into the capsid, two putative phage structural proteins and a protein related to the Int family of site-specific recombinases. Analysis of the N-terminal amino acid sequences of three phage-encoded proteins led to the identification of two genes encoding structural proteins and of peiP, the structural gene of pseudomurein endoisopeptidase. This enzyme is involved in the lysis of host cells, and it appears to belong to a novel enzyme family. peiP was overexpressed in Escherichia coli, and its product was shown to catalyse the in vitro lysis of M. thermoautotrophicum cells. Comparison of the phage psiM2 DNA sequence with parts of the sequence of the wild-type phage psiM1 suggests that psiM2 is a deletion derivative, which formed by homologous recombination between two copies of a direct repeat.

Findings emerging from complete microbial genome sequences

Sixteen microorganisms, including one eukaryote, four archaeons, and 11 eubacteria, have been completely sequenced and published. More than 50 genomes are scheduled to be completed by the year 2000. This explosive growth of information is forcing change in many scientific disciplines (e.g. bioinformatics and molecular genetics), spawning new fields, and even changing the way scientific information is used and shared. Novel, global genome sequence comparisons seem slow to appear but the infrastructure for these projects is being built, and we expect exciting developments in the near future.

Sulfolobus genome: from genomics to biology

Major progress in sequencing the genome of Sulfolobus solfataricus has been closely concerted with the characterization and sequencing of many extrachromosomal genetic elements, including viruses, cryptic plasmids and conjugative plasmids, as well as mobile archaeal introns and transposons. The latter have provided a basis for developing the first generation of vectors that are now being used to study the genetics of Sulfolobus and other Archaea.

Comparing genomes in terms of protein structure: surveys of a finite parts list

We give an overview of the emerging field of structural genomics, describing how genomes can be compared in terms of protein structure. As the number of genes in a genome and the total number of protein folds are both quite limited, these comparisons take the form of surveys of a finite parts list, similar in respects to demographic censuses. Fold surveys have many similarities with other whole-genome characterizations, e.g., analyses of motifs or pathways. However, structure has a number of aspects that make it particularly suitable for comparing genomes, namely the way it allows for the precise definition of a basic protein module and the fact that it has a better defined relationship to sequence similarity than does protein function. An essential requirement for a structure survey is a library of folds, which groups the known structures into 'fold families.' This library can be built up automatically using a structure comparison program, and we described how important objective statistical measures are for assessing similarities within the library and between the library and genome sequences. After building the library, one can use it to count the number of folds in genomes, expressing the results in the form of Venn diagrams and 'top-10' statistics for shared and common folds. Depending on the counting methodology employed, these statistics can reflect different aspects of the genome, such as the amount of internal duplication or gene expression. Previous analyses have shown that the common folds shared between very different microorganisms, i.e., in different kingdoms, have a remarkably similar structure, being comprised of repeated strand-helix-strand super-secondary structure units. A major difficulty with this sort of 'fold-counting' is that only a small subset of the structures in a complete genome are currently known and this subset is prone to sampling bias. One way of overcoming biases is through structure prediction, which can be applied uniformly and comprehensively to a whole genome. Various investigators have, in fact, already applied many of the existing techniques for predicting secondary structure and transmembrane (TM) helices to the recently sequenced genomes. The results have been consistent: microbial genomes have similar fractions of strands and helices even though they have significantly different amino acid composition. The fraction of membrane proteins with a given number of TM helices falls off rapidly with more TM elements, approximately according to a Zipf law. This latter finding indicates that there is no preference for the highly studied 7-TM proteins in microbial genomes. Continuously updated tables and further information pertinent to this review are available over the web at http://bioinfo.mbb.yale.edu/genome.

A phylogenomic study of the MutS family of proteins

The MutS protein of Escherichia coli plays a key role in the recognition and repair of errors made during the replication of DNA. Homologs of MutS have been found in many species including eukaryotes, Archaea and other bacteria, and together these proteins have been grouped into the MutS family. Although many of these proteins have similar activities to the E.coli MutS, there is significant diversity of function among the MutS family members. This diversity is even seen within species; many species encode multiple MutS homologs with distinct functions. To better characterize the MutS protein family, I have used a combination of phylogenetic reconstructions and analysis of complete genome sequences. This phylogenomic analysis is used to infer the evolutionary relationships among the MutS family members and to divide the family into subfamilies of orthologs. Analysis of the distribution of these orthologs in particular species and examination of the relationships within and between subfamilies is used to identify likely evolutionary events (e.g. gene duplications, lateral transfer and gene loss) in the history of the MutS family. In particular, evidence is presented that a gene duplication early in the evolution of life resulted in two main MutS lineages, one including proteins known to function in mismatch repair and the other including proteins known to function in chromosome segregation and crossing-over. The inferred evolutionary history of the MutS family is used to make predictions about some of the uncharacterized genes and species included in the analysis. For example, since function is generally conserved within subfamilies and lineages, it is proposed that the function of uncharacterized proteins can be predicted by their position in the MutS family tree. The uses of phylogenomic approaches to the study of genes and genomes are discussed.

Toprim--a conserved catalytic domain in type IA and II topoisomerases, DnaG-type primases, OLD family nucleases and RecR proteins

Iterative profile searches and structural modeling show that bacterial DnaG-type primases, small primase-like proteins from bacteria and archaea, type IA and type II topoisomerases, bacterial and archaeal nucleases of the OLD family and bacterial DNA repair proteins of the RecR/M family contain a common domain, designated Toprim (topoisomerase-primase) domain. The domain consists of approximately 100 amino acids and has two conserved motifs, one of which centers at a conserved glutamate and the other one at two conserved aspartates (DxD). Examination of the structure of Topo IA and Topo II and modeling of the Toprim domains of the primases reveal a compact beta/alpha fold, with the conserved negatively charged residues juxtaposed, and inserts seen in Topo IA and Topo II. The conserved glutamate may act as a general base in nucleotide polymerization by primases and in strand rejoining by topoisomerases and as a general acid in strand cleavage by topoisomerases and nucleases. The role of this glutamate in catalysis is supported by site-directed mutagenesis data on primases and Topo IA. The DxD motif may coordinate Mg2+that is required for the activity of all Toprim-containing enzymes. The common ancestor of all life forms could encode a prototype Toprim enzyme that might have had both nucleotidyl transferase and polynucleotide cleaving activity.

Skewed oligomers and origins of replication

The putative origin of replication in prokaryotic genomes can be located by a new method that finds short oligomers whose orientation is preferentially skewed around the origin. The skewed oligomer method is shown to work for all bacterial genomes and one of three archaeal genomes sequences to date, confirming known or predicted origins in most cases and in three cases (H. pylori, M. thermoautotrophicum, and Synechocystis sp.), suggesting origins that were previously unknown. In many cases, the presence of conserved genes and nucleotide motifs confirms the predictions. An algorithm for finding these skewed seven-base and eight-base sequences is described, along with a method for combining evidence from multiple skewed oligomers to accurately locate the replication origin. Possible explanations for the phenomenon of skewed oligomers are discussed. Explanations are presented for why some bacterial genomes contain hundreds of highly skewed oligomers, whereas others contain only a handful.

Biochemistry of methanogenesis: a tribute to Marjory Stephenson. 1998 Marjory Stephenson Prize Lecture

Max-Planck-Institut für terrestrische Mikrobiologie, Karl-von-Frisch-Straße, D-35043 Marburg, and Laboratorium für Mikrobiologie, Fachbereich Biologie, Philipps-Universität, Karl-von-Frisch-Straße, D-35032 Marburg, Germany

In 1933, Stephenson & Stickland (1933a) published that they had isolated from river mud, by the single cell technique, a methanogenic organism capable of growth in an inorganic medium with formate as the sole carbon source.

Replicational and transcriptional selection on codon usage in Borrelia burgdorferi

With more than 10 fully sequenced, publicly available prokaryotic genomes, it is now becoming possible to gain useful insights into genome evolution. Before the genome era, many evolutionary processes were evaluated from limited data sets and evolutionary models were constructed on the basis of small amounts of evidence. In this paper, I show that genes on the Borrelia burgdorferi genome have two separate, distinct, and significantly different codon usages, depending on whether the gene is transcribed on the leading or lagging strand of replication. Asymmetrical replication is the major source of codon usage variation. Replicational selection is responsible for the higher number of genes on the leading strands, and transcriptional selection appears to be responsible for the enrichment of highly expressed genes on these strands. Replicational-transcriptional selection, therefore, has an influence on the codon usage of a gene. This is a new paradigm of codon selection in prokaryotes.

Class-directed structure determination: foundation for a protein structure initiative

The recent sequencing of many complete genomes, combined with the development of methods that allow rapid structure determination for many proteins, has changed the way in which protein structure determinations can be approached. One-by-one determinations of individual protein structures will soon be augmented by class-directed structure analyses in which a group of proteins is targeted and structures of representative members are determined and used to represent the entire group. Such a shift in approach would be the foundation for a broad protein structure initiative targeting classes of proteins important for biotechnology and for a fundamental understanding of protein function.

DNA repair in Mycobacterium tuberculosis. What have we learnt from the genome sequence?

The genome sequence of Mycobacterium tuberculosis was analysed by searching for homologues of genes known to be involved in the reversal or repair of DNA damage in Escherichia coli and related organisms. Genes necessary to perform nucleotide excision repair (NER), base excision repair (BER), recombination, and SOS repair and mutagenesis were identified. In particular, all of the genes known to be directly involved in the repair of oxidative and alkylative damage are present in M. tuberculosis. In contrast, we failed to identify homologues of genes involved in mismatch repair. This finding has potentially significant implications with respect to genome stability, strain variability at repeat loci and the emergence of chromosomally encoded drug resistance mutations.

The deoxyxylulose phosphate pathway of terpenoid biosynthesis in plants and microorganisms

Recent studies have uncovered the existence of an alternative, non-mevalonate pathway for the formation of isopentenyl pyrophosphate and dimethylallyl pyrophosphate, the two building blocks of terpene biosynthesis.

Self-identification of protein-coding regions in microbial genomes

A new method for predicting protein-coding regions in microbial genomic DNA sequences is presented. It uses an ab initio iterative Markov modeling procedure to automatically perform the partition of genomic sequences into three subsets shown to correspond to coding, coding on the opposite strand, and noncoding segments. In contrast to current methods, such as GENEMARK [Borodovsky, M. & McIninch, J. D. (1993) Comput. Chem. 17, 123-133], no training set or prior knowledge of the statistical properties of the studied genome are required. This new method tolerates error rates of 1-2% and can process unassembled sequences. It is thus ideal for the analysis of genome survey and/or fragmented sequence data from uncharacterized microorganisms. The method was validated on 10 complete bacterial genomes (from four major phylogenetic lineages). The results show that protein-coding regions can be identified with an accuracy of up to 90% with a totally automated and objective procedure.

The tRNA(guanine-26,N2-N2) methyltransferase (Trm1) from the hyperthermophilic archaeon Pyrococcus furiosus: cloning, sequencing of the gene and its expression in Escherichia coli

The structural gene pfTRM1 (GenBank accession no. AF051912), encoding tRNA(guanine-26, N 2- N 2) methyltransferase (EC 2.1.1.32) of the strictly anaerobic hyperthermophilic archaeon Pyrococcus furiosus, has been identified by sequence similarity to the TRM1 gene of Saccharomyces cerevisiae (YDR120c). The pfTRM1 gene in a 3.0 kb restriction DNA fragment of P.furiosus genomic DNA has been cloned by library screening using a PCR probe to the 5'-part of the corresponding ORF. Sequence analysis revealed an entire ORF of 1143 bp encoding a polypeptide of 381 residues (calculated molecular mass 43.3 kDa). The deduced amino acid sequence of this newly identified gene shares significant similarity with the TRM1- like genes of three other archaea (Methanococcus jannaschii, Methanobacterium thermoautotrophicum and Archaeoglobus fulgidus), one eukaryon (Caenorhabditis elegans) and one hyperthermophilic eubacterium (Aquifex aeolicus). Two short consensus motifs for S-adenosyl-l-methionine binding are detected in the sequence of pfTrm1p. Cloning of the P.furiosus TRM1 gene in an Escherichia coli expression vector allowed expression of the recombinant protein (pfTrm1p) with an apparent molecular mass of 42 kDa. A protein extract from the transformed E.coli cells shows enzymatic activity for the quantitative formation of N 2, N 2-dimethylguanosine at position 26 in a transcript of yeast tRNAPhe used as substrate. The recombinant enzyme was also shown to modify bulk E.coli tRNAs in vivo.

Ambineela, an unusual blue protein isolated from the archaeon Acidianus ambivalens

A novel blue protein, named ambineela, was isolated from the soluble extract of the thermoacidophilic archaeon Acidianus ambivalens. In solution, the purified protein is a monomer with 50 kDa and has a basic character (pI approximately 8.7). The electronic spectrum shows two bands, centred at 395 and 625 nm (A625/A395 = 0.7). The protein does not contain any transition metal; its blue colour is due to an unidentified non-fluorescent cofactor, covalently bound to it. Ambineela N-terminal sequence exhibits a consensus ADP-binding region, suggesting that its unknown cofactor may comprise this molecule or an analogue.

A superfamily of metalloenzymes unifies phosphopentomutase and cofactor-independent phosphoglycerate mutase with alkaline phosphatases and sulfatases

Sequence analysis of the probable archaeal phosphoglycerate mutase resulted in the identification of a superfamily of metalloenzymes with similar metal-binding sites and predicted conserved structural fold. This superfamily unites alkaline phosphatase, N-acetylgalactosamine-4-sulfatase, and cerebroside sulfatase, enzymes with known three-dimensional structures, with phosphopentomutase, 2,3-bisphosphoglycerate-independent phosphoglycerate mutase, phosphoglycerol transferase, phosphonate monoesterase, streptomycin-6-phosphate phosphatase, alkaline phosphodiesterase/nucleotide pyrophosphatase PC-1, and several closely related sulfatases. In addition to the metal-binding motifs, all these enzymes contain a set of conserved amino acid residues that are likely to be required for the enzymatic activity. Mutational changes in the vicinity of these residues in several sulfatases cause mucopolysaccharidosis (Hunter, Maroteaux-Lamy, Morquio, and Sanfilippo syndromes) and metachromatic leucodystrophy.

A glycyl radical solution: oxygen-dependent interconversion of pyruvate formate-lyase

Pyruvate formate-lyase (PFL) catalyses the non-oxidative dissimilation of pyruvate to formate and acetyl-CoA using a radical-chemical mechanism. The enzyme is enzymically interconverted between inactive and active forms, the active form contains an organic free radical located on a glycyl residue in the C-terminal portion of the polypeptide chain. Introduction of the radical into PFL only occurs anaerobically, and the activating enzyme responsible is an iron-sulphur protein that uses S-adenosyl methionine as cofactor and reduced flavodoxin as reductant. As the radical form of PFL is inactivated by molecular oxygen it is safeguarded during the transition to aerobiosis by conversion back to the radical-free, oxygen-stable form. This reaction is catalysed by the anaerobically induced multimeric enzyme alcohol dehydrogenase. The genes encoding PFL and its activating enzyme are adjacent on the chromosome but form discrete transcriptional units. This genetic organization is highly conserved in many, but not all, organisms that have PFL. Recent studies have shown that proteins exhibiting significant similarity to PFL and its activating enzyme are relatively widespread in facultative and obligate anaerobic eubacteria, as well as archaea. The physiological function of many of these PFL-like enzymes remains to be established. It is becoming increasingly apparent that glycyl radical enzymes are more prevalent than previously surmised. They represent a class of enzymes with unusual biochemistry and probably predate the appearance of molecular oxygen.

Deletion of the rbo gene increases the oxygen sensitivity of the sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough

The rbo gene of Desulfovibrio vulgaris Hildenborough encodes rubredoxin oxidoreductase (Rbo), a 14-kDa iron sulfur protein; forms an operon with the gene for rubredoxin; and is preceded by the gene for the oxygen-sensing protein DcrA. We have deleted the rbo gene from D. vulgaris with the sacB mutagenesis procedure developed previously (R. Fu and G. Voordouw, Microbiology 143:1815-1826, 1997). The absence of the rbo-gene in the resulting mutant, D. vulgaris L2, was confirmed by PCR and protein blotting with Rbo-specific polyclonal antibodies. D. vulgaris L2 grows like the wild type under anaerobic conditions. Exposure to air for 24 h caused a 100-fold drop in CFU of L2 relative to the wild type. The lag times of liquid cultures of inocula exposed to air were on average also greater for L2 than for the wild type. These results demonstrate that Rbo, which is not homologous with superoxide dismutase or catalase, acts as an oxygen defense protein in the anaerobic, sulfate-reducing bacterium D. vulgaris Hildenborough and likely also in other sulfate-reducing bacteria and anaerobic archaea in which it has been found.

Archaea and the cell cycle

Sequence similarity data suggest that archaeal chromosome replication is eukaryotic in character. Putative nucleoid-processing proteins display similarities to both eukaryotic and bacterial counterparts, whereas cell division may occur through a predominantly bacterial mechanism. Insights into the organization of the archaeal cell cycle are therefore of interest, not only for understanding archaeal biology, but also for investigating how components from the other two domains interact and work in concert within the same cell; in addition, archaea may have the potential to provide insights into eukaryotic initiation of chromosome replication.

Analogous enzymes: independent inventions in enzyme evolution

It is known that the same reaction may be catalyzed by structurally unrelated enzymes. We performed a systematic search for such analogous (as opposed to homologous) enzymes by evaluating sequence conservation among enzymes with the same enzyme classification (EC) number using sensitive, iterative sequence database search methods. Enzymes without detectable sequence similarity to each other were found for 105 EC numbers (a total of 243 distinct proteins). In 34 cases, independent evolutionary origin of the suspected analogous enzymes was corroborated by showing that they possess different structural folds. Analogous enzymes were found in each class of enzymes, but their overall distribution on the map of biochemical pathways is patchy, suggesting multiple events of gene transfer and selective loss in evolution, rather than acquisition of entire pathways catalyzed by a set of unrelated enzymes. Recruitment of enzymes that catalyze a similar but distinct reaction seems to be a major scenario for the evolution of analogous enzymes, which should be taken into account for functional annotation of genomes. For many analogous enzymes, the bacterial form of the enzyme is different from the eukaryotic one; such enzymes may be promising targets for the development of new antibacterial drugs.

Cloning of the trkAH gene cluster and characterization of the Trk K(+)-uptake system of Vibrio alginolyticus

K(+)-uptake genes of Vibrio alginolyticus were identified by cloning chromosomal DNA fragments of this organism into plasmids, followed by electroporation and selection for growth at low K+ concentrations of cells of an Escherichia coli strain defective in K+ uptake. A 4.1 kb DNA fragment contained a cluster of three ORFs on the same DNA strand: the previously identified trkA gene, a gene similar to E. coli trkH (V. alginolyticus trkH) and a new gene, orf1, whose function is not clear. Products of V. alginolyticus trkA and orf1 were detected in E. coli minicells. trkA and trkH from V. alginolyticus restored growth at low K+ concentrations of an E. coli delta trkA and an E. coli delta trkG delta trkH strain, respectively, suggesting that these V. alginolyticus genes can functionally replace their E. coli counterparts. In addition, a plasmid containing V. alginolyticus trkAH permitted growth of an E. coli delta sapABCDF (delta trkE) strain at low K+ concentrations. This effect was mainly due to V. alginolyticus trkH and was enhanced by trkA from this organism. Measurements of net K(+)-uptake rates indicated that the presence of these genes in E. coli renders the Trk systems independent of products from the E. coli sapABCDF (trkE) operon.

Novel homologs of replication protein A in archaea: implications for the evolution of ssDNA-binding proteins

In Bacteria and Eukarya, ssDNA-binding proteins are central to most aspects of DNA metabolism. Until recently, however, no counterpart of an ssDNA-binding protein had been identified in the third domain of life, Archaea. Here, we report the discovery of a novel type of ssDNA-binding protein in the genomes of several archaeons. These proteins, in contrast to all known members of this protein family, possess four conserved DNA-binding sites within a single polypeptide or, in one case, two polypeptides. This peculiar structural organization allows us to propose a model for the evolution of this class of proteins.