Ancient Systems of Sodium/Potassium Homeostasis as Predecessors of Membrane Bioenergetics.

Cell cytoplasm of archaea, bacteria, and eukaryotes contains substantially more potassium than sodium, and potassium cations are specifically required for many key cellular processes, including protein synthesis. This distinct ionic composition and requirements have been attributed to the emergence of the first cells in potassium-rich habitats. Different, albeit complementary, scenarios have been proposed for the primordial potassium-rich environments based on experimental data and theoretical considerations. Specifically, building on the observation that potassium prevails over sodium in the vapor of inland geothermal systems, we have argued that the first cells could emerge in the pools and puddles at the periphery of primordial anoxic geothermal fields, where the elementary composition of the condensed vapor would resemble the internal milieu of modern cells. Marine and freshwater environments generally contain more sodium than potassium. Therefore, to invade such environments, while maintaining excess of potassium over sodium in the cytoplasm, primordial cells needed means to extrude sodium ions. The foray into new, sodium-rich habitats was the likely driving force behind the evolution of diverse redox-, light-, chemically-, or osmotically-dependent sodium export pumps and the increase of membrane tightness. Here we present a scenario that details how the interplay between several, initially independent sodium pumps might have triggered the evolution of sodium-dependent membrane bioenergetics, followed by the separate emergence of the proton-dependent bioenergetics in archaea and bacteria. We also discuss the development of systems that utilize the sodium/potassium gradient across the cell membranes.

MHYT, a new integral membrane sensor domain.

MHYT, a new conserved protein domain with a likely signaling function, is described. This domain consists of six transmembrane segments, three of which contain conserved methionine, histidine, and tyrosine residues that are projected to lie near the outer face of the cytoplasmic membrane. In Synechocystis sp. PCC6803, this domain forms the N-terminus of the sensor histidine kinase Slr2098. In Pseudomonas aeruginosa and several other organisms, the MHYT domain forms the N-terminal part of a three-domain protein together with previously described GGDEF and EAL domains, both of which have been associated with signal transduction due to their presence in likely signaling proteins. In Bacillus subtilis YkoW protein, an additional PAS domain is found between the MHYT and GGDEF domains. A ykoW null mutant of B. subtilis did not exhibit any growth alterations, consistent with a non-essential, signaling role of this protein. A model of the membrane topology of the MHYT domain indicates that its conserved residues could coordinate one or two copper ions, suggesting a role in sensing oxygen, CO, or NO.

Conserved core structure and active site residues in alkaline phosphatase superfamily enzymes.

Cofactor-independent phosphoglycerate mutase (iPGM) has been previously identified as a member of the alkaline phosphatase (AlkP) superfamily of enzymes, based on the conservation of the predicted metal-binding residues. Structural alignment of iPGM with AlkP and cerebroside sulfatase confirmed that all these enzymes have a common core structure and revealed similarly located conserved Ser (in iPGM and AlkP) or Cys (in sulfatases) residues in their active sites. In AlkP, this Ser residue is phosphorylated during catalysis, whereas in sulfatases the active site Cys residues are modified to formylglycine and sulfatated. Similarly located Thr residue forms a phosphoenzyme intermediate in one more enzyme of the AlkP superfamily, alkaline phosphodiesterase/nucleotide pyrophosphatase PC-1 (autotaxin). Using structure-based sequence alignment, we identified homologous Ser, Thr, or Cys residues in other enzymes of the AlkP superfamily, such as phosphopentomutase, phosphoglycerol transferase, phosphonoacetate hydrolase, and GPI-anchoring enzymes (glycosylphosphatidylinositol phosphoethanolamine transferases) MCD4, GPI7, and GPI13. We predict that catalytical cycles of all the enzymes of AlkP superfamily include phosphoenzyme (or sulfoenzyme) intermediates.

Independent evolution of heavy metal-associated domains in copper chaperones and copper-transporting atpases.

Copper chaperones are small cytoplasmic proteins that bind intracellular copper (Cu) and deliver it to Cu-dependent enzymes such as cytochrome oxidase, superoxide dismutase, and amine oxidase. Copper chaperones are similar in sequence and structure to the Cu-binding heavy metal-associated (HMA) domains of Cu-transporting ATPases (Cu-ATPases), and the genes for copper chaperones and Cu-ATPases are often located in the same operon. Phylogenetic analysis shows that Cu chaperones and HMA domains of Cu-ATPases represent ancient and distinct lineages that have evolved largely independently since their initial separation. Copper chaperone-Cu-ATPase operons appear to have evolved independently in different prokaryotic lineages, probably due to a strong selective pressure for coexpression of these genes.

Novel domains of the prokaryotic two-component signal transduction systems.

The archetypal two-component signal transduction systems include a sensor histidine kinase and a response regulator, which consists of a receiver CheY-like domain and a DNA-binding domain. Sequence analysis of the sensor kinases and response regulators encoded in complete bacterial and archaeal genomes revealed complex domain architectures for many of them and allowed the identification of several novel conserved domains, such as PAS, GAF, HAMP, GGDEF, EAL, and HD-GYP. All of these domains are widely represented in bacteria, including 19 copies of the GGDEF domain and 17 copies of the EAL domain encoded in the Escherichia coli genome. In contrast, these novel signaling domains are much less abundant in bacterial parasites and in archaea, with none at all found in some archaeal species. This skewed phyletic distribution suggests that the newly discovered complexity of signal transduction systems emerged early in the evolution of bacteria, with subsequent massive loss in parasites and some horizontal dissemination among archaea. Only a few proteins containing these domains have been studied experimentally, and their exact biochemical functions remain obscure; they may include transformations of novel signal molecules, such as the recently identified cyclic diguanylate. Recent experimental data provide the first direct evidence of the participation of these domains in signal transduction pathways, including regulation of virulence genes and extracellular enzyme production in the human pathogens Bordetella pertussis and Borrelia burgdorferi and the plant pathogen Xanthomonas campestris. Gene-neighborhood analysis of these new domains suggests their participation in a variety of processes, from mercury and phage resistance to maintenance of virulence plasmids. It appears that the real picture of the complexity of phosphorelay signal transduction in prokaryotes is only beginning to unfold.

Sodium ion cycle in bacterial pathogens: evidence from cross-genome comparisons.

Analysis of the bacterial genome sequences shows that many human and animal pathogens encode primary membrane Na+ pumps, Na+-transporting dicarboxylate decarboxylases or Na+ translocating NADH:ubiquinone oxidoreductase, and a number of Na+ -dependent permeases. This indicates that these bacteria can utilize Na+ as a coupling ion instead of or in addition to the H+ cycle. This capability to use a Na+ cycle might be an important virulence factor for such pathogens as Vibrio cholerae, Neisseria meningitidis, Salmonella enterica serovar Typhi, and Yersinia pestis. In Treponema pallidum, Chlamydia trachomatis, and Chlamydia pneumoniae, the Na+ gradient may well be the only energy source for secondary transport. A survey of preliminary genome sequences of Porphyromonas gingivalis, Actinobacillus actinomycetemcomitans, and Treponema denticola indicates that these oral pathogens also rely on the Na+ cycle for at least part of their energy metabolism. The possible roles of the Na+ cycling in the energy metabolism and pathogenicity of these organisms are reviewed. The recent discovery of an effective natural antibiotic, korormicin, targeted against the Na+ -translocating NADH:ubiquinone oxidoreductase, suggests a potential use of Na+ pumps as drug targets and/or vaccine candidates. The antimicrobial potential of other inhibitors of the Na+ cycle, such as monensin, Li+ and Ag+ ions, and amiloride derivatives, is discussed.

The COG database: new developments in phylogenetic classification of proteins from complete genomes.

The database of Clusters of Orthologous Groups of proteins (COGs), which represents an attempt on a phylogenetic classification of the proteins encoded in complete genomes, currently consists of 2791 COGs including 45 350 proteins from 30 genomes of bacteria, archaea and the yeast Saccharomyces cerevisiae (http://www.ncbi.nlm.nih. gov/COG). In addition, a supplement to the COGs is available, in which proteins encoded in the genomes of two multicellular eukaryotes, the nematode Caenorhabditis elegans and the fruit fly Drosophila melanogaster, and shared with bacteria and/or archaea were included. The new features added to the COG database include information pages with structural and functional details on each COG and literature references, improvements of the COGNITOR program that is used to fit new proteins into the COGs, and classification of genomes and COGs constructed by using principal component analysis.

Conserved 'hypothetical' proteins: new hints and new puzzles.

Conserved hypothetical proteins, i.e. conserved proteins whose functions are still unknown, pose a challenge not just to functional genomics but also to general biology. For many conserved proteins, computational analysis provides only a general prediction of biochemical function; their exact biological functions have to be established through direct experimentation. In the few cases when this has been accomplished, the results were remarkable, revealing the deoxyxylulose pathway and a new essential enzyme, the ITP pyrophosphatase. Comparative genome analysis is also instrumental in illuminating unsolved problems in biology, e.g. the mechanism of FtsZ-independent cell division in Chlamydia, Ureaplasma and Aeropyrum or the role of uncharacterized conserved domains in signal transduction.

The synthetase domains of cobalamin biosynthesis amidotransferases cobB and cobQ belong to a new family of ATP-dependent amidoligases, related to dethiobiotin synthetase.

Phosphotransacetylases of Escherichia coli and several other bacteria contain an additional 350-aa N-terminal fragment that is not required for phosphotransacetylase activity. Sequence analysis of this fragment revealed that it is closely related to a family of ATP-dependent enzymes that also includes dethiobiotin synthetase and the synthetase domains of two amidotransferases involved in cobalamin biosynthesis, cobyrinic acid a,c-diamide synthase (CobB) and cobyric acid synthase (CobQ). Further database searches showed that this enzyme family is also related to the MinD family of ATPases involved in regulation of cell division in bacteria and archaea. Analysis of sequence conservation in the members of this enzyme family using the structure of dethiobiotin synthetase active site as a guide allowed us to suggest a model for the interaction of CobB and CobQ with their respective substrates. CobB and CobQ were also found to contain unusual Triad family (class I) glutamine amidotransferase domains with conserved Cys and His residues, but lacking the Glu residue of the catalytic triad. These results should help in understanding the enzymology of cobalamin biosynthesis and in resolving the role of phosphotransacetylase in regulation of the carbon flow to and from acetate.

Who's your neighbor? New computational approaches for functional genomics.

Several recently developed computational approaches in comparative genomics go beyond sequence comparison. By analyzing phylogenetic profiles of protein families, domain fusions, gene adjacency in genomes, and expression patterns, these methods predict many functional interactions between proteins and help deduce specific functions for numerous proteins. Although some of the resultant predictions may not be highly specific, these developments herald a new era in genomics in which the benefits of comparative analysis of the rapidly growing collection of complete genomes will become increasingly obvious.

Aldolases of the DhnA family: a possible solution to the problem of pentose and hexose biosynthesis in archaea.

Sequence analysis of the recently identified class I aldolase of Escherichia coli (dhnA gene product) helped to identify its homologs in Chlamydia trachomatis, Chlamydiophyla pneumoniae and in each of the completely sequenced archaeal genomes. Iterative database searches revealed sequence similarities between the DhnA-family enzymes, deoxyribose phosphate aldolases and bacterial (class II) fructose bisphosphate aldolases and allowed prediction of similar three-dimensional structures (TIM-barrel fold) in all these enzymes. The Schiff base-forming lysyl residues of DhnA and deoxyribose phosphate aldolase are conserved in all members of the DhnA and deoxyribose phosphate aldolase families, indicating that these enzymes share common features with both class I and class II aldolases. The DhnA-family enzymes are predicted to possess an aldolase activity and to play a critical role in sugar biosynthesis in archaea.

Acetyl-CoA synthetase from the amitochondriate eukaryote Giardia lamblia belongs to the newly recognized superfamily of acyl-CoA synthetases (Nucleoside diphosphate-forming).

The gene coding for the acetyl-CoA synthetase (ADP-forming) from the amitochondriate eukaryote Giardia lamblia has been expressed in Escherichia coli. The recombinant enzyme exhibited the same substrate specificity as the native enzyme, utilizing acetyl-CoA and adenine nucleotides as preferred substrates and less efficiently, propionyl- and succinyl-CoA. N- and C-terminal parts of the G. lamblia acetyl-CoA synthetase sequence were found to be homologous to the alpha- and beta-subunits, respectively, of succinyl-CoA synthetase. Sequence analysis of homologous enzymes from various bacteria, archaea, and the eukaryote, Plasmodium falciparum, identified conserved features in their organization, which allowed us to delineate a new superfamily of acyl-CoA synthetases (nucleoside diphosphate-forming) and its signature motifs. The representatives of this new superfamily of thiokinases vary in their domain arrangement, some consisting of separate alpha- and beta-subunits and others comprising fusion proteins in alpha-beta or beta-alpha orientation. The presence of homologs of acetyl-CoA synthetase (ADP-forming) in such human pathogens as G. lamblia, Yersinia pestis, Bordetella pertussis, Pseudomonas aeruginosa, Vibrio cholerae, Salmonella typhi, Porphyromonas gingivalis, and the malaria agent P. falciparum suggests that they might be used as potential drug targets.

Towards understanding the first genome sequence of a crenarchaeon by genome annotation using clusters of orthologous groups of proteins (COGs).

BACKGROUND: Standard archival sequence databases have not been designed as tools for genome annotation and are far from being optimal for this purpose. We used the database of Clusters of Orthologous Groups of proteins (COGs) to reannotate the genomes of two archaea, Aeropyrum pernix, the first member of the Crenarchaea to be sequenced, and Pyrococcus abyssi. RESULTS: A. pernix and P. abyssi proteins were assigned to COGs using the COGNITOR program; the results were verified on a case-by-case basis and augmented by additional database searches using the PSI-BLAST and TBLASTN programs. Functions were predicted for over 300 proteins from A. pernix, which could not be assigned a function using conventional methods with a conservative sequence similarity threshold, an approximately 50% increase compared to the original annotation. A. pernix shares most of the conserved core of proteins that were previously identified in the Euryarchaeota. Cluster analysis or distance matrix tree construction based on the co-occurrence of genomes in COGs showed that A. pernix forms a distinct group within the archaea, although grouping with the two species of Pyrococci, indicative of similar repertoires of conserved genes, was observed. No indication of a specific relationship between Crenarchaeota and eukaryotes was obtained in these analyses. Several proteins that are conserved in Euryarchaeota and most bacteria are unexpectedly missing in A. pernix, including the entire set of de novo purine biosynthesis enzymes, the GTPase FtsZ (a key component of the bacterial and euryarchaeal cell-division machinery), and the tRNA-specific pseudouridine synthase, previously considered universal. A. pernix is represented in 48 COGs that do not contain any euryarchaeal members. Many of these proteins are TCA cycle and electron transport chain enzymes, reflecting the aerobic lifestyle of A. pernix. CONCLUSIONS: Special-purpose databases organized on the basis of phylogenetic analysis and carefully curated with respect to known and predicted protein functions provide for a significant improvement in genome annotation. A differential genome display approach helps in a systematic investigation of common and distinct features of gene repertoires and in some cases reveals unexpected connections that may be indicative of functional similarities between phylogenetically distant organisms and of lateral gene exchange.

Using the COG database to improve gene recognition in complete genomes.

A complete understanding of the biology of an organism necessarily starts with knowledge of its genetic makeup. Proteins encoded in a genome must be identified and characterized, and the presence or absence of specific sets of proteins must be noted in order to determine the possible biochemical pathways or functional systems utilized by that organism. The COG database presents a set of tools suited to these purposes, including the ability to select protein families (COGs) that contain proteins from a specified set of species. The selection is based upon a phylogenetic pattern, which is a shorthand representation of the presence or absence of a particular species in a COG. Here we present the use of phylogenetic patterns as a means to perform targeted searches for undetected protein-coding genes in complete genomes.

The COG database: a tool for genome-scale analysis of protein functions and evolution.

Rational classification of proteins encoded in sequenced genomes is critical for making the genome sequences maximally useful for functional and evolutionary studies. The database of Clusters of Orthologous Groups of proteins (COGs) is an attempt on a phylogenetic classification of the proteins encoded in 21 complete genomes of bacteria, archaea and eukaryotes (http://www. ncbi.nlm. nih.gov/COG). The COGs were constructed by applying the criterion of consistency of genome-specific best hits to the results of an exhaustive comparison of all protein sequences from these genomes. The database comprises 2091 COGs that include 56-83% of the gene products from each of the complete bacterial and archaeal genomes and approximately 35% of those from the yeast Saccharomyces cerevisiae genome. The COG database is accompanied by the COGNITOR program that is used to fit new proteins into the COGs and can be applied to functional and phylogenetic annotation of newly sequenced genomes.

Searching for drug targets in microbial genomes.

Comparative analysis of the complete genome sequences of 10 bacterial pathogens available in the public databases offers the first insights into the drug discovery approaches of the near future. Genes that are conserved in different genomes often turn out to be essential, which makes them attractive targets for new broad-spectrum antibiotics. Subtractive genome analysis reveals the genes that are conserved in all or most of the pathogenic bacteria but not in eukaryotes; these are the most obvious candidates for drug targets. Species-specific genes, on the other hand, may offer the possibility to design drugs against a particular, narrow group of pathogens.

Comparative genomics of the Archaea (Euryarchaeota): evolution of conserved protein families, the stable core, and the variable shell.

Comparative analysis of the protein sequences encoded in the four euryarchaeal species whose genomes have been sequenced completely (Methanococcus jannaschii, Methanobacterium thermoautotrophicum, Archaeoglobus fulgidus, and Pyrococcus horikoshii) revealed 1326 orthologous sets, of which 543 are represented in all four species. The proteins that belong to these conserved euryarchaeal families comprise 31%-35% of the gene complement and may be considered the evolutionarily stable core of the archaeal genomes. The core gene set includes the great majority of genes coding for proteins involved in genome replication and expression, but only a relatively small subset of metabolic functions. For many gene families that are conserved in all euryarchaea, previously undetected orthologs in bacteria and eukaryotes were identified. A number of euryarchaeal synapomorphies (unique shared characters) were identified; these are protein families that possess sequence signatures or domain architectures that are conserved in all euryarchaea but are not found in bacteria or eukaryotes. In addition, euryarchaea-specific expansions of several protein and domain families were detected. In terms of their apparent phylogenetic affinities, the archaeal protein families split into bacterial and eukaryotic families. The majority of the proteins that have only eukaryotic orthologs or show the greatest similarity to their eukaryotic counterparts belong to the core set. The families of euryarchaeal genes that are conserved in only two or three species constitute a relatively mobile component of the genomes whose evolution should have involved multiple events of lineage-specific gene loss and horizontal gene transfer. Frequently these proteins have detectable orthologs only in bacteria or show the greatest similarity to the bacterial homologs, which might suggest a significant role of horizontal gene transfer from bacteria in the evolution of the euryarchaeota.