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

Molecular cloning and functional expression of a protein-serine/threonine phosphatase from the hyperthermophilic archaeon Pyrodictium abyssi TAG11

An open reading frame coding for a putative protein-serine/threonine phosphatase was identified in the hyperthermophilic archaeon Pyrodictium abyssi TAG11 and named Py-PP1. Py-PP1 was expressed in Escherichia coli, purified from inclusion bodies, and biochemically characterized. The phosphatase gene is part of an operon which may provide, for the first time, insight into a physiological role for archaeal protein phosphatases in vivo.

Functional analysis of the secretory precursor processing machinery of Bacillus subtilis: identification of a eubacterial homolog of archaeal and eukaryotic signal peptidases

Approximately 47% of the genes of the Gram-positive bacterium Bacillus subtilis belong to paralogous gene families. The present studies were aimed at the functional analysis of the sip gene family of B. subtilis, consisting of five chromosomal genes, denoted sipS, sipT, sipU, sipV, and sipW. All five sip genes specify type I signal peptidases (SPases), which are actively involved in the processing of secretory preproteins. Interestingly, strains lacking as many as four of these SPases could be obtained. As shown with a temperature-sensitive SipS variant, only cells lacking both SipS and SipT were not viable, which may be caused by jamming of the secretion machinery with secretory preproteins. Thus, SipS and SipT are of major importance for protein secretion. This conclusion is underscored by the observation that only the transcription of the sipS and sipT genes is temporally controlled via the DegS-DegU regulatory system, in concert with the transcription of most genes for secretory preproteins. Notably, the newly identified SPase SipW is highly similar to SPases from archaea and the ER membrane of eukaryotes, suggesting that these enzymes form a subfamily of the type I SPases, which is conserved in the three domains of life.

A novel aminopeptidase associated with the 60 kDa chaperonin in the thermophilic archaeon Sulfolobus solfataricus

The chaperonins are high-molecular-weight protein complexes having a characteristic double-ring toroidal shape; they are thought to aid the folding of denatured or newly synthesized polypeptides. These proteins exist as two functionally similar but distantly related families, one including the bacterial and organellar chaperonins and the other (termed the CCT-TRiC family) including the chaperonins of the Archaea and the eukaryotes. The CCT-TRiC chaperonins, particularly their archeal members, are less well known than their bacterial counterparts, and their main cellular function is still doubtful. In this work, we report that the chaperonin of the thermophilic archaeon Sulfolobus solfataricus interacts with several polypeptides other than the two subunits that constitute the 18-mer double-ring structure. We have cloned and sequenced the gene encoding one 90 kDa chaperonin-associated protein and have shown, using biochemical assays, that the product is an enzyme belonging to the family of zinc-dependent aminopeptidases. The Sulfolobus protein shows maximal homology to eukaryotic (yeast and mouse) aminopeptidases. It contains a leucine zipper motif and can be phosphorylated by an unidentified kinase present in the cell extracts. The possible significance of an association between an aminopeptidase and a chaperonin is discussed.

Molecular biology of hyperthermophilic Archaea

The sequences of a number of archaeal genomes have recently been completed, and many more are expected shortly. Consequently, the research of Archaea in general and hyperthermophiles in particular has entered a new phase, with many exciting discoveries to be expected. The wealth of sequence information has already led, and will continue to lead to the identification of many enzymes with unique properties, some of which have potential for industrial applications. Subsequent functional genomics will help reveal fundamental matters such as details concerning the genetic, biochemical and physiological adaptation of extremophiles, and hence give insight into their genomic evolution, polypeptide structure-function relations, and metabolic regulation. In order to optimally exploit many unique features that are now emerging, the development of genetic systems for hyperthermophilic Archaea is an absolute requirement. Such systems would allow the application of this class of Archaea as so-called "cell factories": (i) expression of certain archaeal enzymes for which no suitable conventional (mesophilic bacterial or eukaryal) systems are available, (ii) selection for thermostable variants of potentially interesting enzymes from mesophilic origin, and (iii) the development of in vivo production systems by metabolic engineering. An overview is given of recent insight in the molecular biology of hyperthermophilic Archaea, as well as of a number of promising developments that should result in the generation of suitable genetic systems in the near future.

C1 transfer enzymes and coenzymes linking methylotrophic bacteria and methanogenic Archaea

Methanogenic and sulfate-reducing Archaea are considered to have an energy metabolism involving C1 transfer coenzymes and enzymes unique for this group of strictly anaerobic microorganisms. An aerobic methylotrophic bacterium, Methylobacterium extorquens AM1, was found to contain a cluster of genes that are predicted to encode some of these enzymes and was shown to contain two of the enzyme activities and one of the methanogenic coenzymes. Insertion mutants were all unable to grow on C1 compounds, suggesting that the archaeal enzymes function in aerobic C1 metabolism. Thus, methylotrophy and methanogenesis involve common genes that cross the bacterial/archaeal boundaries.

Molecular characterization and postsplicing fate of three introns within the single rRNA operon of the hyperthermophilic archaeon Aeropyrum pernix K1

The single rRNA operon (arnS-arnL) of the hyperthermophilic archaeon Aeropyrum pernix K1 was sequenced. The DNA sequence data and detailed RNA analyses disclosed an unusual feature: the presence of three introns at hitherto undescribed insertion positions within the rRNA genes. The 699-nucleotide (nt) intron Ialpha was located at position 908 (Escherichia coli numbering [H. F. Noller, Annu. Rev. Biochem. 53:119-162, 1984]) of the 16S rRNA, while the 202-nt intron Ibeta and 575-nt intron Igamma were located at positions 1085 and 1927 (E. coli numbering), respectively, of the 23S rRNA. They were located within highly conserved sites which have been implicated as crucial for rRNA function in E. coli. All three introns were remarkably AT rich (41.5 to 43.1 mol% G+C) compared with the mature rRNAs (67.7 and 69.2 mol% G+C for 16S and 23S rRNAs, respectively). No obvious primary sequence similarities were detected among them. After splicing from rRNA transcripts in vivo, a large quantity of intronic RNAs were stably retained in the linear monomeric form, whereas a trace of topoisomeric RNA molecules also appeared, as characterized by their behavior in two-dimensional gel electrophoresis. Secondary structural models of the Ialpha-, Ibeta-, and Igamma-containing rRNA precursors agree with the bulge-helix-bulge motif. Two of the introns, Ialpha and Igamma, contained open reading frames whose protein translation exhibited no overall similarity with proteins reported so far. However, both share a LAGLI-DADG motif characteristic of homing endonucleases.

Clustered genes encoding the methyltransferases of methanogenesis from monomethylamine

Coenzyme M (CoM) is methylated during methanogenesis from monomethyamine in a reaction catalyzed by three proteins. Using monomethylamine, a 52-kDa polypeptide termed monomethylamine methyltransferase (MMAMT) methylates the corrinoid cofactor bound to a second polypeptide, monomethylamine corrinoid protein (MMCP). Methylated MMCP then serves as a substrate for MT2-A, which methylates CoM. The genes for these proteins are clustered on 6.8 kb of DNA in Methanosarcina barkeri MS. The gene encoding MMCP (mtmC) is located directly upstream of the gene encoding MMAMT (mtmB). The gene encoding MT2-A (mtbA) was found 1.1 kb upstream of mtmC, but no obvious open reading frame was found in the intergenic region between mtbA and mtmC. A single monocistronic transcript was found for mtbA that initiated 76 bp from the translational start. Separate transcripts of 2.4 and 4.7 kb were detected, both of which carried mtmCB. The larger transcript also encoded mtmP, which is homologous to the APC family of cationic amine permeases and may therefore encode a methylamine permease. A single transcriptional start site was found 447 bp upstream of the translational start of mtmC. MtmC possesses the corrinoid binding motif found in corrinoid proteins involved in dimethylsulfide- and methanol-dependent methanogenesis, as well as in methionine synthase. The open reading frame of mtmB was interrupted by a single in-frame, midframe, UAG codon which was also found in mtmB from M. barkeri NIH. A mechanism that circumvents UAG-directed termination of translation must operate during expression of mtmB in this methanogen.

N-terminal processing: the methionine aminopeptidase and N alpha-acetyl transferase families

Removal of the initiator methionine and/or acetylation of the alpha-amino group are among the earliest possible chemical modifications that occur during protein synthesis in eukaryotes. These events are catalyzed by methionine aminopeptidase and N alpha-acetyltransferase, respectively. Recent advances in the isolation and characterization of these enzymes indicate that they exist as isoforms that vary in cellular location, function, and evolutionary origins.

Aeromonas hydrophila adenylyl cyclase 2: a new class of adenylyl cyclases with thermophilic properties and sequence similarities to proteins from hyperthermophilic archaebacteria

Complementation of an Escherichia coli cya mutant with a genomic library from Aeromonas hydrophila allowed isolation of clones containing two different cya genes. Whereas one of these genes (cyaA) coded for an adenylyl cyclase (AC1) belonging to the previously described class I adenylyl cyclases (ACs), the second one (cyaB) coded for a protein (AC2) that did not match any previously characterized protein when compared to protein sequence databases. In particular, it did not align with any of members of the three known classes of ACs. The purified AC2 enzyme exhibited remarkable biochemical characteristics, namely, an optimum activity at a high temperature (65 degrees C) and at an alkalinic pH (9.5). In order to investigate the functions of both cyclases in A. hydrophila, each gene was inactivated in the chromosome and the resulting mutant strains were examined for physiological alterations. It was shown that, in contrast to cyaA, the cyaB gene was not expressed under usual laboratory growth conditions. However, introduction of a plasmid harboring the cyaB gene in a cyaA mutant, as well as in a cyaA cyaB mutant, allowed cyclic AMP production. AC2 is the first member of a new class of previously unrecognized ACs, and to date, no functional counterpart has been demonstrated in other organisms. However, scanning databases revealed a significant similarity between AC2 and the gene product of three hyperthermophilic archaebacteria: Methanobacterium thermoautotrophicum, Archaeglobus fulgidus, and Methanococcus jannaschii. The possibility of a gene transfer between such phylogenetically divergent bacteria is discussed.

The A1A0 ATPase from Methanosarcina mazei: cloning of the 5' end of the aha operon encoding the membrane domain and expression of the proteolipid in a membrane-bound form in Escherichia coli

Three additional ATPase genes, clustered in the order ahaH, ahaI, and ahaK, were found upstream of the previously characterized genes ahaECFABDG coding for the archaeal A1A0 ATPase from Methanosarcina mazei. ahaH, the first gene in the cluster, is preceded by a conserved promoter sequence. Northern blot analysis revealed that the clusters ahaHIK and ahaECFABDG are transcribed as one message. AhaH is a hydrophilic polypeptide and is similar to peptides of previously unassigned function encoded by genes preceding postulated ATPase genes in Methanobacterium thermoautotrophicum and Methanococcus jannaschii. AhaI has a two-domain structure with a hydrophilic domain of 39 kDa and a hydrophobic domain with seven predicted transmembrane alpha helices. It is similar to the 100-kDa polypeptide of V1V0 ATPases and is therefore suggested to participate in proton transport. AhaK is a hydrophobic polypeptide with two predicted transmembrane alpha helices and, on the basis of sequence comparisons and immunological studies, is identified as the proteolipid, a polypeptide which is essential for proton translocation. However, it is only one-half and one-third the size of the proteolipids from M. thermoautotrophicum and M. jannaschii, respectively. ahaK is expressed in Escherichia coli, and it is incorporated into the cytoplasmic membrane despite the different chemical natures of lipids from archaea and bacteria. This is the first report on the expression and incorporation into E. coli lipids of a membrane integral enzyme from a methanogens, which will facilitate analysis of the structure and function of the membrane domain of the methanoarchaeal ATPase.

Conserved sequence motifs among bacterial, eukaryotic, and archaeal phosphatases that define a new phosphohydrolase superfamily

Members of a new molecular family of bacterial nonspecific acid phosphatases (NSAPs), indicated as class C, were found to share significant sequence similarities to bacterial class B NSAPs and to some plant acid phosphatases, representing the first example of a family of bacterial NSAPs that has a relatively close eukaryotic counterpart. Despite the lack of an overall similarity, conserved sequence motifs were also identified among the above enzyme families (class B and class C bacterial NSAPs, and related plant phosphatases) and several other families of phosphohydrolases, including bacterial phosphoglycolate phosphatases, histidinol-phosphatase domains of the bacterial bifunctional enzymes imidazole-glycerolphosphate dehydratases, and bacterial, eukaryotic, and archaeal phosphoserine phosphatases and threalose-6-phosphatases. These conserved motifs are clustered within two domains, separated by a variable spacer region, according to the pattern [FILMAVT]-D-[ILFRMVY]-D-[GSNDE]-[TV]-[ILVAM]-[AT S VILMC]-X-¿YFWHKR)-X-¿YFWHNQ¿-X( 102,191)-¿KRHNQ¿-G-D-¿FYWHILVMC¿-¿QNH¿-¿FWYGP¿-D -¿PSNQYW¿. The dephosphorylating activity common to all these proteins supports the definition of this phosphatase motif and the inclusion of these enzymes into a superfamily of phosphohydrolases that we propose to indicate as "DDDD" after the presence of the four invariant aspartate residues. Database searches retrieved various hypothetical proteins of unknown function containing this or similar motifs, for which a phosphohydrolase activity could be hypothesized.

Biosynthesis of 2-aceto-2-hydroxy acids: acetolactate synthases and acetohydroxyacid synthases

Two groups of enzymes are classified as acetolactate synthase (EC 4. 1.3.18). This review deals chiefly with the FAD-dependent, biosynthetic enzymes which readily catalyze the formation of acetohydroxybutyrate from pyruvate and 2-oxobutyrate, as well as of acetolactate from two molecules of pyruvate (the ALS/AHAS group). These enzymes are generally susceptible to inhibition by one or more of the branched-chain amino acids which are ultimate products of the acetohydroxyacids, as well as by several classes of herbicides (sulfonylureas, imidazolinones and others). Some ALS/AHASs also catalyze the (non-physiological) oxidative decarboxylation of pyruvate, leading to peracetic acid; the possible relationship of this process to oxygen toxicity is considered. The bacterial ALS/AHAS which have been well characterized consist of catalytic subunits (around 60 kDa) and smaller regulatory subunits in an alpha2beta2 structure. In the case of Escherichia coli isozyme III, assembly and dissociation of the holoenzyme has been studied. The quaternary structure of the eukaryotic enzymes is less clear and in plants and yeast only catalytic polypeptides (homologous to those of bacteria) have been clearly identified. The presence of regulatory polypeptides in these organisms cannot be ruled out, however, and genes which encode putative ALS/AHAS regulatory subunits have been identified in some cases. A consensus sequence can be constructed from the 21 sequences which have been shown experimentally to represent ALS/AHAS catalytic polypeptides. Many other sequences fit this consensus, but some genes identified as putative 'acetolactate synthase genes' are almost certainly not ALS/AHAS. The solution of the crystal structures of several thiamin diphosphate (ThDP)-dependent enzymes which are homologous to ALS/AHAS, together with the availability of many amino acid sequences for the latter enzymes, has made it possible for two laboratories to propose similar, reasonable models for a dimer of catalytic subunits of an ALS/AHAS. A number of characteristics of these enzymes can now be better understood on the basis of such models: the nature of the herbicide binding site, the structural role of FAD and the binding of ThDP-Mg2+. The models are also guides for experimental testing of ideas concerning structure-function relationships in these enzymes, e.g. the nature of the substrate recognition site. Among the important remaining questions is how the enzyme suppresses alternative reactions of the intrinsically reactive hydroxyethylThDP enamine formed by the decarboxylation of the first substrate molecule and specifically promotes its condensation with 2-oxobutyrate or pyruvate.

Derivation of the three-dimensional architecture of bacterial ribonuclease P RNAs from comparative sequence analysis

The secondary structure of bacterial RNase P RNA, a ribozyme responsible for the maturation of the 5' end of tRNAs, is well established on the basis of sequence comparison analysis. RNase P RNA secondary structures fall into two types, A and B, which share a common core formed by the assembly of two main folding domains, but differ in their peripheral elements.A revised alignment of 137 available sequences reveals new covariations allowing for the refinement of both types of secondary structures. Phylogenetic evidence is thus provided for the extension of stems P11, P14, P19, P10.1 and P15.1 through further canonical base-pairs or GAellipsisGA mismatches. These refinements led in turn to a new organization of the catalytic core, with coaxial stackings of helices P2 and P19 as well as P1 and P4. New inter-domain tertiary interactions involve loop L9 and helix P1 and loop L8 with helix P4. These features were incorporated into atomic-scale 3D models of RNase P RNA for representatives of each structural type, namely Escherichia coli and Bacillus subtilis. In each model, the juxtaposition of the core helices creates a cradle onto which the pre-tRNA substrate binds with most evolutionarily conserved residues converging towards the cleavage site. The inner cores of both types are stabilized similarly, albeit by different peripheral elements, emphasizing the modular and hierarchical organisation of the architecture of RNase P RNAs. Similarities are thus apparent between the type A modules, P16/P17/P6 and P13/P14, and their type B analogs, P5.1/P15.1 and P10. 1/P10.1a, respectively. Other noteworthy features of these models include compactness and good agreement with published crosslinking data.

Beyond complete genomes: from sequence to structure and function

Computer analysis of complete prokaryotic genomes shows that microbial proteins are in general highly conserved--approximately 70% of them contain ancient conserved regions. This allows us to delineate families of orthologs across a wide phylogenetic range and, in many cases, predict protein functions with considerable precision. Sequence database searches using newly developed, sensitive algorithms result in the unification of such orthologous families into larger superfamilies sharing common sequence motifs. For many of these superfamilies, prediction of the structural fold and specific amino acid residues involved in enzymatic catalysis is possible. Taken together, sequence and structure comparisons provide a powerful methodology that can successfully complement traditional experimental approaches.

Novel evolutionary histories and adaptive features of proteins from hyperthermophiles

The hyperthermophiles include both bacteria and archaea, although the majority of isolates growing above 100 degreesC are archaea. Newly described adaptive features of hyperthermophiles include proteins stable to 200 degreesC, nucleosomes, chaperonins and high-capacity DNA modifying enzymes. The ongoing release of genomic sequence data from hyperthermophiles will continue to accelerate the discovery of novel proteins.

Genomics and computational molecular biology

There has been a dramatic increase in the number of completely sequenced bacterial genomes during the past two years as a result of the efforts both of public genome agencies and the pharmaceutical industry. The availability of completely sequenced genomes permits more systematic analyses of genes, evolution and genome function than was otherwise possible. Using computational methods - which are used to identify genes and their functions including statistics, sequence similarity, motifs, profiles, protein folds and probabilistic models - it is possible to develop characteristic genome signatures, assign functions to genes, identify pathogenic genes, identify metabolic pathways, develop diagnostic probes and discover potential drug-binding sites. All of these directions are critical to understanding bacterial growth, pathogenicity and host-pathogen interactions.

Measuring genome evolution

The determination of complete genome sequences provides us with an opportunity to describe and analyze evolution at the comprehensive level of genomes. Here we compare nine genomes with respect to their protein coding genes at two levels: (i) we compare genomes as "bags of genes" and measure the fraction of orthologs shared between genomes and (ii) we quantify correlations between genes with respect to their relative positions in genomes. Distances between the genomes are related to their divergence times, measured as the number of amino acid substitutions per site in a set of 34 orthologous genes that are shared among all the genomes compared. We establish a hierarchy of rates at which genomes have changed during evolution. Protein sequence identity is the most conserved, followed by the complement of genes within the genome. Next is the degree of conservation of the order of genes, whereas gene regulation appears to evolve at the highest rate. Finally, we show that some genomes are more highly organized than others: they show a higher degree of the clustering of genes that have orthologs in other genomes.

Analyzing genomes with cumulative skew diagrams

A novel method of cumulative diagrams shows that the nucleotide composition of a microbial chromosome changes at two points separated by about a half of its length. These points coincide with sites of replication origin and terminus for all bacteria where such sites are known. The leading strand is found to contain more guanine than cytosine residues. This fact is used to predict origin and terminus locations in other bacterial and archaeal genomes. Local changes, visible as diagram distortions, may represent recent genome rearrangements, as demonstrated for two strains of Escherichia coli . Analysis of the diagrams of viral and mitochondrial genomes suggests a link between the base composition bias and the time spent by DNA in a single stranded state during replication.

Cytoskeletal proteins: the evolution of cell division

The prokaryotic cell division protein FtsZ and eukaryotic tubulin have been shown to have very similar structures and are most likely homologs. The evolutionary transition from FtsZ to tubulin could provide a window into the transition from prokaryotic cells to eukaryotic cells.

RadA protein is an archaeal RecA protein homolog that catalyzes DNA strand exchange

With the discovery that the Saccharomyces cerevisiae Rad51 protein is both structurally and functionally similar to the Escherichia coli RecA protein, the RecA paradigm for homologous recombination was extended to the Eucarya. The ubiquitous presence of RecA and Rad51 protein homologs raises the question of whether this archetypal protein exists within the third domain of life, the Archaea. Here we present the isolation of a Rad51/RecA protein homolog from the archaeon Sulfolobus solfataricus, and show that this protein, RadA, possesses the characteristics of a DNA strand exchange protein: The RadA protein is a DNA-dependent ATPase, forms a nucleoprotein filament on DNA, and catalyzes DNA pairing and strand exchange.

Domain structure and RNA annealing activity of the Escherichia coli regulatory protein StpA

The Escherichia coli regulatory protein StpA bears striking similarity to the chromatin-associated protein H-NS. These two proteins have many structural, functional and mechanistic parallels. Although H-NS is more abundant in the cell, both proteins act as transcriptional regulators, both bind to curved DNA and both restrain DNA supercoils. However, StpA is better able to promote RNA annealing and trans-splicing in vitro. In this study, phylogenetic analyses and experiments to examine the protease sensitivity of StpA and H-NS suggest a similar structure for the two proteins. Both proteins consist of two structured domains separated by an exposed protease-sensitive linker. The N-terminal (StpA-NterL) and C-terminal (StpA-CterL) domains of StpA, as well as the full-length StpA and H-NS proteins, were cloned, overproduced in E. coli and purified to homogeneity. StpA-CterL, but not StpA-NterL, promotes strand annealing of complementary RNA oligonucleotides and in vitro trans-splicing of a model group I intron. Both StpA and StpA-CterL exhibited stronger RNA-modulating activity than H-NS. Phylogenetic analyses showed that the N-terminal and C-terminal domains can exist autonomously. The phylogenetic and experimental data are compatible with a two-domain model for StpA and H-NS, with independently functioning modules joined by a non-conserved linker, and with the observed RNA-related activities residing entirely within the C-terminal domain.

A transcription frame-based analysis of the genomic DNA sequence of a hyper-thermophilic archaeon for the identification of genes, pseudo-genes and operon structures

An algorithm for identifying transcription units, independently regulated genes and operons, and pseudo-genes that are not expected to be expressed, has been developed by combining a system for predicting transcription and translation signals, and a system for scoring the triplet periodicity in ORF candidates. By using the algorithm, the 1.09 Mb sequence that covers approximately 60% of the genome of Pyrococcus sp. OT3 has been analyzed. The identified ORFs show the expected biological and physical characteristics, while the rejected ORF candidates do not. Frequent use of operon structures for transcription, and gene duplication followed by mutation or termination of the duplicated genes, are discussed.

Genomics: re-evaluation of translation machinery evolution

Experiments based on genome sequence analysis have revealed unexpected complexity in the evolution of the translation apparatus, including concerted evolution of Gln-tRNA synthetase and Glu-tRNAGln amidotransferase, and a novel, class I Lys-tRNA synthetase shared by archaea and spirochaetes.

PPi-dependent phosphofructokinase from Thermoproteus tenax, an archaeal descendant of an ancient line in phosphofructokinase evolution

Flux into the glycolytic pathway of most cells is controlled via allosteric regulation of the irreversible, committing step catalyzed by ATP-dependent phosphofructokinase (PFK) (ATP-PFK; EC 2.7.1.11), the key enzyme of glycolysis. In some organisms, the step is catalyzed by PPi-dependent PFK (PPi-PFK; EC 2.7.1.90), which uses PPi instead of ATP as the phosphoryl donor, conserving ATP and rendering the reaction reversible under physiological conditions. We have determined the enzymic properties of PPi-PFK from the anaerobic, hyperthermophilic archaeon Thermoproteus tenax, purified the enzyme to homogeneity, and sequenced the gene. The approximately 100-kDa PPi-PFK from T. tenax consists of 37-kDa subunits; is not regulated by classical effectors of ATP-PFKs such as ATP, ADP, fructose 2,6-bisphosphate, or metabolic intermediates; and shares 20 to 50% sequence identity with known PFK enzymes. Phylogenetic analyses of biochemically characterized PFKs grouped the enzymes into three monophyletic clusters: PFK group I represents only classical ATP-PFKs from Bacteria and Eucarya; PFK group II contains only PPi-PFKs from the genus Propionibacterium, plants, and amitochondriate protists; whereas group III consists of PFKs with either cosubstrate specificity, i.e., the PPi-dependent enzymes from T. tenax and Amycolatopsis methanolica and the ATP-PFK from Streptomyces coelicolor. Comparative analyses of the pattern of conserved active-site residues strongly suggest that the group III PFKs originally bound PPi as a cosubstrate.

A novel DNA polymerase family found in Archaea

One of the most puzzling results from the complete genome sequence of the methanogenic archaeon Methanococcus jannaschii was that the organism may have only one DNA polymerase gene. This is because no other DNA polymerase-like open reading frames (ORFs) were found besides one ORF having the typical alpha-like DNA polymerase (family B). Recently, we identified the genes of DNA polymerase II (the second DNA polymerase) from the hyperthermophilic archaeon Pyrococcus furiosus, which has also at least one alpha-like DNA polymerase (T. Uemori, Y. Sato, I. Kato, H. Doi, and Y. Ishino, Genes Cells 2:499-512, 1997). The genes in M. jannaschii encoding the proteins that are homologous to the DNA polymerase II of P. furiosus have been located and cloned. The gene products of M. jannaschii expressed in Escherichia coli had both DNA polymerizing and 3'-->5' exonuclease activities. We propose here a novel DNA polymerase family which is entirely different from other hitherto-described DNA polymerases.

Strand compositional asymmetry in bacterial and large viral genomes

Several bacterial genomes exhibit preference for G over C on the DNA leading strand extending from the origin of replication to the ter-region in the genomes of Escherichia coli, Mycoplasma genitalium, Bacillus subtilis, and marginally in Haemophilus influenzae, Mycoplasma pneumoniae, and Helicobacter pylori. Strand compositional asymmetry is not observed in the cyanobacterium Synechocystis sp. genome nor in the archaeal genomes of Methanococcus jannaschii, Methanobacterium thermoautotrophicum, and Archaeoglobus fulgidus. A strong strand compositional asymmetry is observed in beta-type but not alpha- or gamma-type human herpesviruses featuring G > C downstream of oriL and C > G upstream of oriL. Dinucleotide relative abundances (i.e., dinucleotide representations normalized by the component nucleotide frequencies) are consonant with respect to the leading and lagging strands. Strand compositional asymmetry may reflect on differences in replication synthesis of the leading versus lagging strand, on differences between template and coding strand associated with transcription-coupled repair mechanisms, on differences in gene density between the two strands, on differences in residue and codon biases in relation to gene function, expression level, or operon organization, or on differences in single or context-dependent base mutational rates. The absence of strand asymmetry in the archaeal genomes may reflect the presence of multiple origins of replication.

Archaeal translation initiation revisited: the initiation factor 2 and eukaryotic initiation factor 2B alpha-beta-delta subunit families

As the amount of available sequence data increases, it becomes apparent that our understanding of translation initiation is far from comprehensive and that prior conclusions concerning the origin of the process are wrong. Contrary to earlier conclusions, key elements of translation initiation originated at the Universal Ancestor stage, for homologous counterparts exist in all three primary taxa. Herein, we explore the evolutionary relationships among the components of bacterial initiation factor 2 (IF-2) and eukaryotic IF-2 (eIF-2)/eIF-2B, i.e., the initiation factors involved in introducing the initiator tRNA into the translation mechanism and performing the first step in the peptide chain elongation cycle. All Archaea appear to posses a fully functional eIF-2 molecule, but they lack the associated GTP recycling function, eIF-2B (a five-subunit molecule). Yet, the Archaea do posses members of the gene family defined by the (related) eIF-2B subunits alpha, beta, and delta, although these are not specifically related to any of the three eukaryotic subunits. Additional members of this family also occur in some (but by no means all) Bacteria and even in some eukaryotes. The functional significance of the other members of this family is unclear and requires experimental resolution. Similarly, the occurrence of bacterial IF-2-like molecules in all Archaea and in some eukaryotes further complicates the picture of translation initiation. Overall, these data lend further support to the suggestion that the rudiments of translation initiation were present at the Universal Ancestor stage.

A tRNA(Glu) gene from the hyperthermophilic archaeon Pyrococcus furiosus contains the 3'-terminal CCA sequence of the mature tRNA

We cloned a gene encoding tRNA(Glu) of the hyperthermophilic archaeon Pyrococcus furiosus. This gene contains the CCA sequence corresponding to the 3'-terminus of the mature tRNA. It is known that, like in eukaryal tRNAs, the CCA-termini of archaeal tRNAs are generally not encoded. Therefore, we analyzed all tRNA genes in the genome of Methanococcus jannaschii estimated by its whole genome sequence. Twenty-one of 37 listed tRNA genes contained the 3'-terminal CCA sequence. The corresponding M. jannaschii tRNA(Glu) gene does not contain the CCA sequence, although the tRNA sequences of the M. jannaschii and P. furiosus tRNA(Glu) genes are 86% identical.

Archaeal genomics: do archaea have a mixed heritage?

A third complete archaeal genome sequence, replete with eukaryote-like genes for replication, transcription and translation, has appeared. The sequence also shows bacteria-like features. It is time to come to grips with this evidence for a mixed heritage.

GeneMark.hmm: new solutions for gene finding

The number of completely sequenced bacterial genomes has been growing fast. There are computer methods available for finding genes but yet there is a need for more accurate algorithms. The GeneMark. hmm algorithm presented here was designed to improve the gene prediction quality in terms of finding exact gene boundaries. The idea was to embed the GeneMark models into naturally derived hidden Markov model framework with gene boundaries modeled as transitions between hidden states. We also used the specially derived ribosome binding site pattern to refine predictions of translation initiation codons. The algorithm was evaluated on several test sets including 10 complete bacterial genomes. It was shown that the new algorithm is significantly more accurate than GeneMark in exact gene prediction. Interestingly, the high gene finding accuracy was observed even in the case when Markov models of order zero, one and two were used. We present the analysis of false positive and false negative predictions with the caution that these categories are not precisely defined if the public database annotation is used as a control.

Biosynthesis of the ansamycin antibiotic rifamycin: deductions from the molecular analysis of the rif biosynthetic gene cluster of Amycolatopsis mediterranei S699

BACKGROUND

The ansamycin class of antibiotics are produced by various Actinomycetes. Their carbon framework arises from the polyketide pathway via a polyketide synthase (PKS) that uses an unusual starter unit. Rifamycin (rif), produced by Amycolatopsis mediterranei, is the archetype ansamycin and it is medically important. Although its basic precursors (3-amino-5-hydroxy benzoic acid AHBA, and acetic and propionic acids) had been established, and several biosynthetic intermediates had been identified, very little was known about the origin of AHBA nor had the PKS and the various genes and enzymes that modify the initial intermediate been characterized.

RESULTS

A set of 34 genes clustered around the rifK gene encoding AHBA synthase were defined by sequencing all but 5 kilobases (kb) of a 95 kb contiguous region of DNA from A. mediterranei. The involvement of some of the genes in the biosynthesis of rifamycin B was examined. At least five genes were shown to be essential for the synthesis of AHBA, five genes were determined to encode the modular type I PKS that uses AHBA as the starter unit, and 20 or more genes appear to govern modification of the polyketide-derived framework, and rifamycin resistance and export. Putative regulatory genes were also identified. Disruption of the PKS genes at the end of rifA abolished rifamycin B production and resulted in the formation of P8/1-OG, a known shunt product of rifamycin biosynthesis, whereas disruption of the orf6 and orf9 genes, which may encode deoxysugar biosynthesis enzymes, had no apparent effect.

CONCLUSIONS

Rifamycin production in A. mediterranei is governed by a single gene cluster consisting of structural, resistance and export, and regulatory genes. The genes characterized here could be modified to produce novel forms of the rifamycins that may be effective against rifamycin-resistant microorganisms.

Universally conserved translation initiation factors

The process by which translation is initiated has long been considered similar in Bacteria and Eukarya but accomplished by a different unrelated set of factors in the two cases. This not only implies separate evolutionary histories for the two but also implies that at the universal ancestor stage, a translation initiation mechanism either did not exist or was of a different nature than the extant processes. We demonstrate herein that (i) the "analogous" translation initiation factors IF-1 and eIF-1A are actually related in sequence, (ii) the "eukaryotic" translation factor SUI1 is universal in distribution, and (iii) the eukaryotic/archaeal translation factor eIF-5A is homologous to the bacterial translation factor EF-P. Thus, the rudiments of translation initiation would seem to have been present in the universal ancestor stage. However, significant development and refinement subsequently occurred independently on both the bacterial lineage and on the archaeal/eukaryotic line.

Modular organization of inteins and C-terminal autocatalytic domains

Analysis of the conserved sequence features of inteins (protein "introns") reveals that they are composed of three distinct modular domains. The N-terminal (N) and C-terminal (C) domains are predicted to perform different parts of the autocatalytic protein splicing reaction. An optional endonuclease domain (EN) is shown to correspond to different types of homing endonucleases in different inteins. The N domain contains motifs predicted to catalyze the first steps of protein splicing, leading to the cleavage of the intein N terminus from its protein host. Intein N domain motifs are also found in C-terminal autocatalytic domains (CADs) present in hedgehog and other protein families. Specific residues in the N domain of intein and CADs are proposed to form a charge relay system involved in cleaving their N-termini. The intein C domain is apparently unique to inteins and contains motifs that catalyze the final protein splicing steps: ligation of the intein flanks and cleavage of its C terminus to release the free intein and spliced host protein. All intein EN domains known thus far have dodecapeptide (DOD, LAGLI-DADG) type homing endonuclease motifs. This work identifies an EN domain with an HNH homing-endonuclease motif and two new small inteins with no EN domains. One of these small inteins might be inactive or a "pseudo intein." The results suggest a modular architecture for inteins, clarify their origin and relationship to other protein families, and extend recent experimental findings on the functional roles of intein N, C, and EN motifs.

Archaeal-type lysyl-tRNA synthetase in the Lyme disease spirochete Borrelia burgdorferi

Lysyl-tRNAs are essential for protein biosynthesis by ribosomal mRNA translation in all organisms. They are synthesized by lysyl-tRNA synthetases (EC 6.1.1.6), a group of enzymes composed of two unrelated families. In bacteria and eukarya, all known lysyl-tRNA synthetases are subclass IIc-type aminoacyl-tRNA synthetases, whereas some archaea have been shown to contain an unrelated class I-type lysyl-tRNA synthetase. Examination of the preliminary genomic sequence of the bacterial pathogen Borrelia burgdorferi, the causative agent of Lyme disease, indicated the presence of an open reading frame with over 55% similarity at the amino acid level to archaeal class I-type lysyl-tRNA synthetases. In contrast, no coding region with significant similarity to any class II-type lysyl-tRNA synthetase could be detected. Heterologous expression of this open reading frame in Escherichia coli led to the production of a protein with canonical lysyl-tRNA synthetase activity in vitro. Analysis of B. burgdorferi mRNA showed that the lysyl-tRNA synthetase-encoding gene is highly expressed, confirming that B. burgdorferi contains a functional class I-type lysyl-tRNA synthetase. The detection of an archaeal-type lysyl-tRNA synthetase in B. burgdorferi and other pathogenic spirochetes, but not to date elsewhere in bacteria or eukarya, indicates that the gene that encodes this enzyme has a common origin with its orthologue from the archaeal kingdom. This difference between the lysyl-tRNA synthetases of spirochetes and their hosts may be readily exploitable for the development of anti-spirochete therapeutics.

Archaea: what can we learn from their sequences?

Gene-by-gene and traditional biochemical approaches continue to reveal surprising molecular features in the archaeal domain. In addition, the complete sequencing of several archaeal genomes has further confirmed the phenotypic coherence of these micro-organisms at the molecular level. Nevertheless, the phylogeny of Archaea and the nature of the last universal common ancestor are still matters for debate.