INTRODUCTION

TOP ABSTRACT INTRODUCTION RESULTS AND DISCUSSION METHODS WEB SITE REFERENCES REFERENCES

Biotin (vitamin H) is an essential cofactor for a class of important metabolic enzymes, biotin carboxylases and decarboxylases (Perkins and Pero 2001). The biotin biosynthetic pathway is widespread among microorganisms. The well-studied systems of biotin biosynthesis from Escherichia coli, Bacillus subtilis, and Bacillus sphaericus differ in the first step of biosynthesis. B. subtilis and B. sphaericus use pimeloyl-CoA synthase encoded by the bioW gene to synthesize pimeloyl-CoA from pimelic acid. In addition, pimelic acid formation in B. subtilis has been proposed to use cytochrome P450 encoded by bioI (Stok and De Voss 2000). In E. coli, pimeloyl-CoA is synthesized from L-alanine and/or acetate via acetyl-CoA, instead of pimelic acid (Ifuku et al. 1994), and products of the bioC and bioH genes are required for pimeloyl-CoA synthesis in E. coli. The pathway from pimeloyl-CoA to biotin is similar in E. coli and bacilli and uses products of the bioF, bioD, bioA, and bioB genes (Fig. 1). Genes encoding biotin transporters have not been identified in bacteria until now, but E. coli can uptake biotin by active transport (Piffeteau and Gaudry 1985), and a gene for biotin transport, bioP, has been mapped on the E. coli chromosome (Eisenberg 1985).

View larger version (37K): [in this window] [in a new window]   Figure 1   The biotin biosynthesis pathway in bacteria.

The operon organization of the biotin biosynthetic genes differs between E. coli and bacilli. E. coli has bioBFCD operon located divergently with the bioA gene and single bioH gene (DeMoll 1994). In contrast, B. subtilis has the single bioWAFDBI operon (Perkins et al. 1996). Two unlinked biotin biosynthetic operons, bioDAYB and bioXWF, were described in B. sphaericus (Gloeckler et al. 1990). The functions of two new biotin-related genes, bioX and bioY, are presently unknown; however, it has been proposed that BioX of B. sphaericus and BioC of E. coli may function as acyl carrier proteins involved in the pimeloyl-CoA synthesis (Lemoine et al. 1996). Recently, four biotin biosynthetic gene clusters, orf1-bioDA, orf2-bioFB, bioH-orf3, and bioFIIHIIC, were characterized in Gram-positive bacterium Kurthia sp. (Kiyasu et al. 2001). The authors of this study suggested that, in contrast to B. subtilis and B. sphaericus, Kurthia sp. produces pimeloyl-CoA by a pathway similar to that of E. coli.

The biotin operon of E. coli is negatively regulated by biotin and the bifunctional protein BirA (DeMoll 1994). The biotin-protein ligase BirA mediates biotinylation of acetyl-CoA carboxylase via a two-step reaction. Firstly, the adenylate of biotin is synthesized from substrates biotin and ATP and, at the second step, transferred to a unique lysine residue on carboxylase. When biotin is unclaimed, two generated BirA-biotinyl-5'-AMP monomers bind cooperatively to the bioO operator between the divergent bioA and bioBCDF operons and repress transcription in both directions. The BirA protein is composed of the N-terminal DNA-binding (D-b) domain containing a helix-turn-helix (HTH) structure, the central domain, and the C-terminal domain. The central catalytic domain contains the binding site for biotinyl-5'-AMP and also is required for transcriptional regulation (Kwon et al. 2000). The BirA protein of B. subtilis has a similar structure and also can act as a repressor of the bioWAFDBI operon (Bower et al. 1996). Recently, two new BirA-regulated operons of unknown function, yhfUST and yuiG, were detected in B. subtilis by expression microarray analysis (Lee et al. 2001). Imperfect palindromic sequences, which are partially similar to the bioO operator from E. coli, were found upstream of the BirA-regulated operons from B. subtilis, B. sphaericus, and Kurthia sp. (Gloeckler et al. 1990; Kiyasu et al. 2001; Lee et al. 2001).

The large number of complete genomes now available provides an opportunity to perform global comparison of whole metabolic pathways and regulons in a variety of bacteria. The comparative analysis of binding sites for transcriptional regulators in bacterial genomes is a powerful approach to functional annotation of genomes (for review, see Gelfand 1999). The general assumption in such studies is that true sites mostly occur upstream of orthologous genes, whereas false positives are scattered at random in the genome. In addition, analysis of gene clustering on the chromosome allows one to detect functionally coupled genes (Overbeek et al. 1999).

Here, we report the comparative study of the biotin regulon and metabolic pathway in all available prokaryotic genomes. It is shown that birA is the most widely distributed biotin-related gene in bacteria. However, only a fraction of BirA orthologs possess the N-terminal D-b domain with the HTH motif (D-b-BirA). Presence of D-b-BirA in a genome coincides with occurrence of potential BirA sites upstream of biotin-related genes. The BirA-mediated regulation was found in such diverse bacterial lineages as proteobacteria, low-GC Gram-positive bacteria, and archaea. At that, BirA is the only transcriptional regulator with the binding signal conserved in eubacteria and archaea. On the practical side, this analysis allowed us to predict new members of biotin regulons, to assign biotin-transport function to BioY, and to detect nonorthologous displacement of bioH in several lineages and individual genomes.

	RESULTS AND DISCUSSION

TOP ABSTRACT INTRODUCTION RESULTS AND DISCUSSION METHODS WEB SITE REFERENCES REFERENCES

Orthologs of birA and biotin biosynthetic genes (BBS) from E. coli and B. subtilis were identified in all available bacterial genomes by similarity search (Table 1). The biotin-protein ligase BirA is widely distributed in eubacteria and archaea. Only Buchnera sp., Borrelia burgdorferi, Aeropyrum pernix, thermoplasmas, and mycoplasmas have neither the BBS genes nor birA, which is consistent with the lack of biotin-dependent carboxylases in the genomes of these microorganisms. The BBS genes are less widespread than birA: among all complete genomes, Sinorhizobium meliloti, Rickettsia prowazekii, Deinococcus radiodurans, Thermotoga maritima, Treponema pallidum, most archaea, and Gram-positive pathogens from the Bacillus/Clostridium group lack the BBS genes, but have birA. Among archaeal genomes, only Methanococcus jannaschii has a cluster of the BBS genes. Phylogenetic analysis of the BBS proteins shows that this archaeal BBS gene cluster may be the result of possible horizontal gene transfer from bacilli. The detailed phylogenetic and positional analysis of the BBS genes is given below.

BirA Regulon

To analyze possible transcriptional regulation of the BBS genes, we started with identification of the N-terminal regulatory domains in the detected BirA proteins. Using multiple alignment, we compiled the list of 46 sequences of the BirA N-terminal domains that have the same length as the known regulatory domain of E. coli BirA. To determine the significance of the possible helix-turn-helix (HTH) regulatory motif in each of the collected sequences, the HTH motif prediction program (Dodd and Egan 1990) was used (Fig. 2). After that, eight sequences without HTH motifs were removed, and 38 BirA proteins with the predicted DNA-binding regulatory domains (D-b-BirA) were retained (Table 1). We also retained the BirA protein from Bacillus cereus, although it was predicted to contain no HTH motif. This looks like a false-negative prediction. Indeed, not only is BirA highly conserved among bacilli, but the B. cereus genome has several strong BirA sites upstream of biotin-related operons. To support the selection of D-b-BirA, the phylogenetic tree of 50 BirA N-terminal domains was constructed (Fig. 3). It shows that each sequence without a potential HTH motif is highly diverged from the D-b-BirA sequences and looks like an outgroup in this tree.

View larger version (137K): [in this window] [in a new window]   Figure 2   Multiple alignment of the BirA N-terminal domains and identification of the HTH motif. The known secondary structure of the Escherichia coli BirA is shown in the first row. The 2 and 3 helices form the helix-turn-helix (HTH) structure. The score and the probability of the candidate HTH motif are given. A score of <2.5 is not significant. Non-HTH proteins are boxed, except BirA from Bacillus cereus, which is a false-negative prediction (see text). The genome abbreviations are listed in Table 1.

View larger version (13K): [in this window] [in a new window]   Figure 3   Maximum likelihood tree of the N-terminal domains of BirA. Domains containing the regulatory HTH motif are shown in solid lines. Other N-terminal domains of BirA (without HTH) are shown as outgroups by broken lines. The proteobacterial and nonproteobacterial subtrees are separated by arrowtail signs. The genome abbreviations are listed in Table 1.

D-b-BirA is widely distributed in the Bacillus/Clostridium group, gamma-proteobacteria, and archaea. In addition, it was found in Nitrosomonas europaea, Methylobacillus flagellatus, Magnetococcus sp., and Thermus thermophilus. The N-terminal domains of BirA from the Pasteurellaceae family of gamma-proteobacteria possibly have lost their regulatory function. The genomes of Clostridium acetobutylicum, Lactococcus lactis, Halobacterium sp., Pyrococcus abyssii, and Pyrococcus furiosus have two BirA paralogs, with and without the N-terminal regulatory domain. The phylogenetic analysis of the catalytic BirA domains shows that paralogous BirA in the first three genomes could result from a recent duplication. In P. abyssii and P. furiosus, BirA without the N-terminal regulatory domain is close to the other archaeal BirA, whereas the second BirA (D-b-BirA) has a weakly conserved catalytic domain and a well-conserved N-terminal regulatory domain.

Based on the phylogenetic analysis of the D-b domains, all D-b-BirAs were divided into two major groups, proteobacterial and nonproteobacterial (Fig. 3). Consistent with this, two different recognition rules (profiles) for the BirA sites were constructed using the sets of upstream regions of the BBS genes from various genomes. The BirA profile for proteobacteria (with consensus 5'-tTGTaAACC-N14 ... 16-GGTTtACAa-3', where strongly conserved positions are shown in capitals) is more strict than that for other bacteria (5'wwTGTtAAC-N14 ... 16-GTTaACAww-3', where `w' stands for A or T). The constructed profiles were used to detect new candidate members of the BirA regulons in the genomes containing D-b-BirA. Proteobacteria possess only one strong BirA site per genome occurring upstream of the BBS operon. However, most Gram-positive bacteria and some archaea have multiple BirA sites located upstream of BBS genes and new genes of the BirA regulon (Table 1). For a control, we checked the genomes without D-b-BirA for the existence of BirA sites upstream of the BBS operons, and found none.

After comparison of the BirA regulons from numerous bacteria, we predicted several new biotin-regulated genes. A gene of unknown function, bioY (so named by Gloeckler et al. 1990), is widely distributed in bacteria and often clusters with genes of biotin metabolism. The homologs of BioY form a unique protein family (InterPro entry IPR003784), and have no significant similarity to any gene of known function. Analysis of the BirA sites showed that bioY is always under regulation of the biotin repressor in genomes containing regulatory D-b-BirA. The existence of the BirA-regulated bioY in several complete genomes that have no BBS genes indicates that bioY is probably not involved in biotin biosynthesis. On the other hand, proteins of the BioY family have six candidate transmembrane segments, an arrangement typical for prokaryotic transporters. The phylogenetic tree of the BioY protein family consists of several branches, and within each branch most members are positionally linked to BBS genes, or have upstream candidate BirA-binding sites, or both (Fig. 4A). Taken together, these observations strongly imply that all BioY paralogs are transporters of biotin or some biotin precursor.

View larger version (37K): [in this window] [in a new window]   Figure 4   Maximum likelihood trees of the predicted biotin-related transporter BioY (A), the hypothetical long-chain-fatty acid-CoA ligase YhfT (B), and the hypothetical acetyl-CoA-acetyltransferase YhfS (C). Genes predicted to be regulated by BirA are boxed and shown in bold. The co-occurrence of the bioY, yhfS, and yhfT genes in one genome is shown by thick lines. Background colors signify: (black) single bioY gene; (blue) bioY from the biotin biosynthetic operon; (red) bioY in one operon with cbiO-cbiQ; (yellow) bioY in one operon with lspA; (magenta) bioY in the fadH-fabZ-fabK-bioY-fabD operon; (green) bioY positionally linked to the yhfS-yhfT gene pair. The bioY genes positionally linked to birA are shown by broken lines. The genome abbreviations are listed in Table 1.

Another gene pair of unknown function, yhfS-yhfT, has been detected in several bacteria from the Bacillus/Clostridium group and in S. meliloti. Except for the latter genome, the yhfS-yhfT genes are always under predicted regulation by BirA. YhfT and YhfS are homologous to numerous long-chain fatty acid-CoA ligases and acetyl-CoA-acetyltransferases, respectively. Each of them forms a separate branch on the phylogenetic tree for the corresponding protein family (Fig. 4B,C). One of the bioY paralogs from B. subtilis, yhfU, belongs to the yhfUST operon, and transcription of this operon is repressed by BirA (Lee et al. 2001). In addition, yhfU and yhfS-yhfT are clustered in the genomes of B. cereus, Lactococcus lactis, Clostridium difficile, and S. meliloti; whereas Streptococcus pyogenes, Streptococcus equi, and Staphylococcus aureus have separate BirA-regulated yhfST and yhfU operons. Surprisingly, all YhfU paralogs except one from C. difficile form a separate branch in the phylogenetic tree of the BioY family (Fig. 4A). Again, occurrence of the positionally linked yhfU-yhfS-yhfT genes in complete genomes without BBS genes rules out their involvement in the first steps of biotin biosynthesis. A plausible hypothesis is that the YhfS-YhfT proteins are involved in fatty acid metabolism, the pathway that requires biotin at one of the early steps (cf. clustering of bioY with fatty acid biosynthetic genes in T. maritima; see below).

Positional Analysis of Biotin Genes

To reveal new biotin-related genes, we analyzed putative operon structures and chromosomal clustering of the BBS, birA, and bioY genes. In some eubacterial and archaeal genomes, bioY is clustered with a hypothetical two-component ABC cassette that encodes ATPase and permease components from the CbiO and CbiQ families, respectively (Table 1; Fig. 4A). The cbiN-cbiO-cbiQ operon of Salmonella typhimurium encodes the permease, ATPase, and the second permease components, respectively, of a putative cobalt transporter (Roth et al. 1993). Analysis of the phylogenetic trees for the CbiO and CbiQ protein families shows the existence of separate tree branches for the bioY-linked CbiO and CbiQ components of putative ABC transporters from S. meliloti, R. capsulatus, Agrobacterium tumefaciens, Bordetella pertussis, Thermomonospora fusca, two corynebacteria, and D. radiodurans (data not shown). The bioY genes from T. pallidum, Halobacterium sp., and Archaeoglobus fulgidus form possible operons with cbiO homologs and hypothetical transmembrane proteins (with six predicted TMS) that are not similar to any known protein. Both Methanosarcina genomes have BirA-regulated bioY-cbiO1-cbiO2-cbiQ operons encoding two paralogous ATPase components from the CbiO family. Computational approaches alone cannot explain the possible functional link between the predicted biotin transporter BioY and the putative ABC transporter CbiO-CbiQ, but the obtained data seem to be sufficiently strong to warrant experimental analysis.

Another interesting finding is that bioY from T. maritima was found in one operon with genes involved in fatty acid biosynthesis (Table 1). One logical explanation of this linkage is that fatty acid biosynthesis requires biotin as a coenzyme for a hypothetical biotin carboxylase. In addition, positional linkage of the bioY gene with a hypothetical signal peptidase lspA was observed in all cyanobacteria; the functional meaning of this observation is unclear.

Some differences in the gene organization and BirA-mediated regulation of the bioY genes were observed in three Pyrococcus genomes. Strong BirA sites in the common regulatory regions of divergently transcribed bioY and birA genes were predicted in the genomes of P. abyssii and P. furiosus. Besides the regulatory birA gene, these two genomes also contain the second birA gene, encoding BirA without the regulatory domain. In contrast, Pyrococcus horikoshii has no regulatory birA gene, and BirA sites were not found in this genome.

We predicted possible coregulation of various biotin-dependent carboxylases and BirA in some genomes (Table 1). The pycA and pycB genes encoding the biotin-dependent pyruvate carboxylase were found in one candidate operon with birA in two Methanosarcina genomes. These Methanosarcina operons and the single pycA gene from A. fulgidus are preceded by weak BirA sites. The genes encoding subunits of putative propionyl-CoA carboxylase (pccA and pccB) are clustered on the chromosome with the birA gene in all actinobacteria and Halobacterium sp. Finally, in R. capsulatus, birA is located within a long gene cluster encoding components of the malonate decarboxylase Na⁺ pump. The BirA-regulated gene clusters from C. acetobutylicum, L. lactis, and some archaea contain the birA gene itself; therefore, the biotin repressors from these bacteria can be autoregulated.

The bioC-bioH gene pair is required for the synthesis of pimeloyl-CoA in E. coli. The bioC gene is widely distributed in bacteria, whereas bioH was not found in many bioC-containing bacterial genomes. Instead, we predict several nonorthologous gene displacements of bioH in some of these genomes. It was recently shown that the bioZ gene from the bioABFDZ operon of Mesorhizobium loti can complement bioH of E. coli (Sullivan et al. 2001). The orthologs of bioZ with the same gene organization were found in A. tumefaciens and Brucella melitensis.

Using comparative analysis, we have detected displacement of bioH by another gene, named here bioG, in some proteobacteria (including all Pasteurellaceae), the CFB group of bacteria, and Fusobacterium nucleatum (Table 1). The bioG gene always forms an operon with bioC and other BBS genes in these genomes; furthermore, in Bacteroides fragilis there is a single gene encoding a fused protein BioC-BioG. Interestingly, all gamma-proteobacteria except Pasteurellaceae possess the bioC-bioH gene pair, whereas all Pasteurellaceae have bioC-bioG. Neisseria meningitidis has both bioC-bioH and bioC-bioG gene pairs, and the latter likely has been acquired from Haemophilus influenzae or a closely related bacterium, as the respective genes are highly similar. The phylogenetic tree of the BioC family has a separate branch for the proteins associated with BioG (Fig. 5).

View larger version (39K): [in this window] [in a new window]   Figure 5   Maximum likelihood tree of BioC. The proteins predicted to be associated with (blue) BioH, (red) BioG, (yellow) BioZ, and (green) BioK. The genome abbreviations are listed in Table 1.

Another bioC-linked gene, named bioK, was found in two cyanobacteria, Synechococcus sp. and Prochlorococcus marinus. The genomes of these bacteria contain the bioFKCDA operon and the bioB gene. Two other cyanobacteria, Synechocystis sp. and Nostoc sp., have all biotin biosynthetic genes except bioC and bioK. Therefore, they possibly use a different pathway for pimeloyl-CoA synthesis.

Using similarity search, we detected that BioC possesses an S-adenosylmethionine binding motif (InterPro entry IPR000379) and belongs to the methyltransferase superfamily. BioK and BioG are not similar to any known protein. The BioZ protein is similar to the 3-oxoacyl-[acyl-carrier-protein] synthase FabH involved in fatty acid biosynthesis in bacteria. Another BioC-linked protein, BioH, possesses the active-site serine of a wide variety of enzymes including esterases, lipases, and peptidases (InterPro entry IPR000379) and is similar to arylesterase EstE from Pseudomonas fluorescens (26% identity). All bioK and bioG genes, as well as most bioH genes, are located immediately upstream of the bioC gene in the biotin operon.

The observed diversity of enzymes for the first step of biotin biosynthesis can reflect either frequent nonorthologous gene displacements, or possible use of different substrates for biotin biosynthesis. In contrast, B. subtilis, S. aureus, Corynebacterium diphtheriae, Aquifex aeolicus, and M. jannaschii possess pimeloyl-CoA synthase encoded by the bioW gene and can use pimelate as a biotin precursor (Table 1).

It remains unclear why the comparative analysis of regulation and operon structures failed to identify missing BBS genes in the complete genomes of Clostridium perfringens and C. acetobutylicum. The former has no the bioF and bioA counterparts, whereas the latter lacks only bioF. However, these bacteria possess the predicted biotin transporter BioY. It would be interesting to check if these bacteria can synthesize biotin de novo, and if they can, to search for genes missing in their incomplete BBS pathways.

Conclusions

The biotin-protein ligase BirA is a ubiquitous enzyme in bacteria. In addition, BirA can act as a repressor of transcription when it has the N-terminal DNA-binding domain. Using a global analysis of BirA proteins and DNA-binding sites in available bacterial genomes, we have found that the BirA regulon is widely distributed in eubacteria and archaea. A correlation exists between the presence of D-b-BirA and finding of the BirA sites in bacterial genomes. Conservation of the BirA binding sites across large phylogenetic distances allows us to suggest that D-b-BirA is the first example of an ancient DNA-binding transcriptional factor common to eubacteria and archaea. It is unlikely that numerous BirA regulons in various archaea result from mass gene transfer from bacteria, as this scenario would involve many similar, but independent events (although some cases of horizontal transfer are very clear). In contrast, analysis of regulatory systems for biosynthesis of riboflavin and thiamin showed that they are operated by conserved RNA elements, the RFN element (Vitreschak et al. 2002) and the Thi-box (Miranda-Rios et al. 2001), respectively. These unique regulatory elements are widely distributed in eubacteria and, in addition, several Thi-boxes have been found in archaeal genomes (Vitreschak et al. 2002). Thus, it seems very likely that, in general, the regulatory systems for vitamin biosynthesis are ancient.

Comparative analysis of the biotin regulon in complete genomes resulted in new functional assignments for the bioY, yhfS, and yhfT genes. The first of them, bioY, widely distributed in eubacteria and archaea, is a member of the BirA regulon in all genomes containing D-b-BirA, and it has been predicted to encode a transporter for biotin or biotin-related compounds. Proteins YhfS and YhfT, associated with BioY, can be involved in the metabolic pathway that requires biotin as a coenzyme. The systematic comparison of putative operon structures revealed the conserved gene string bioY-cbiO-cbiQ in some bacterial genomes. Such functional linkage between the putative ABC transporter CbiO-CbiQ and the biotin transporter BioY is enigmatic.

Positional analysis resulted in dissection of novel interesting examples of coregulation of biotin-related genes. Positional linkage between birA and genes encoding biotin-dependent carboxylases was found in Actinobacteria and some archaea, and a fraction of these genes were predicted to be regulated by the biotin repressor. Several genomes have divergently transcribed birA and bioY genes with predicted BirA sites in their common regulatory region. Another example of coregulation of bioY with genes of fatty acid biosynthesis in T. maritima can be easily explained, as biotin is a required cofactor of carboxylase, the latter being involved in the first step of fatty acid biosynthesis.

The enzymes mediating the first step of the biotin biosynthetic pathway are diverse. BioW and BioC represent two major types of enzymes involved in the synthesis of pimeloyl-CoA, a biotin precursor. Moreover, another type of pimeloyl-CoA synthetase, namely, PauA, was found recently in Pseudomonas mendocina (Binieda et al. 1999). In contrast to BioW, PauA belongs to the newly recognized superfamily of acyl-CoA synthetases (Sanchez et al. 2000) and is involved in catabolism rather than biosynthesis. The most interesting observation is that various bacteria have different BioC-associated proteins (BioH, BioG, BioK, or BioZ). It can be explained either by utilization of different sources for biotin biosynthesis or by nonorthologous displacements of the BioC-linked proteins.

This report once again shows the power of comparative genomics for prediction of regulatory sites and functional annotation of genomes, especially when experimental data are limited. In particular, this approach is a powerful tool for prediction of missing transport genes, shown by this study and in the analysis of riboflavin (Vitreschak et al. 2002) and thiamin (A. Vitreschak, D. Rodionov, A. Mironov, and M. Gelfand, in prep.) regulons.

	METHODS

TOP ABSTRACT INTRODUCTION RESULTS AND DISCUSSION METHODS WEB SITE REFERENCES REFERENCES

Complete and partial bacterial genomes were downloaded from GenBank (Benson et al. 2000). Preliminary sequence data were also obtained from the Web sites of the Institute for Genomic Research (http://www.tigr.org), the University of Oklahoma's Advanced Center for Genome Technology (http://www.genome.ou.edu/), the Wellcome Trust Sanger Institute (http://www.sanger.ac.uk/), the DOE Joint Genome Institute (http://jgi.doe.gov), and the ERGO database (Overbeek et al. 2000; http://ergo.integratedgenomics.com/ERGO/). The gene identifiers from the ERGO database and GenBank are used throughout.

The existence of BirA with an N-terminal DNA-binding domain (D-b-BirA) is a prerequisite to the comparative analysis of the BirA regulons in bacteria. Therefore, the bacterial genomes containing D-b-BirA were selected and divided into two major groups, proteobacterial and nonproteobactertial including archaeal, according to the phylogenetic tree of the DNA-binding domains of D-b-BirA (Fig. 3). Two training sets were composed; each of them included the upstream regions of the biotin biosynthetic genes (operons) from one of the above genomic groups.

For construction of the BirA profiles, we used the "inverted repeat" option in the SignalX program (Mironov et al. 2000) with a 14-16-bp spacer between two 9-bp units of the inverted repeat. The positional nucleotide weights in the profile were defined as
where N(b,k) is the count of nucleotide b in position k (Mironov et al. 1999). The score of a candidate site was calculated as the sum of the respective positional nucleotide weights:
where L is the length of the site. All genomes containing D-b-BirA were scanned using the constructed profiles, and the genes with candidate regulatory sites in the upstream regions were selected.

Protein alignment was performed using the Smith-Waterman algorithm implemented in the GenomeExplorer program (Mironov et al. 2000). Orthologous proteins were defined by the best-bidirectional-hits criterion (Tatusov et al. 2000). Distant homologs were identified using PSI-BLAST (Altschul et al. 1997). Multiple sequence alignments were constructed using CLUSTALX (Thompson et al. 1997). Phylogenetic trees were created by the maximum likelihood method implemented in PHYLIP (Felsenstein 1981) and drawn using the GeneMaster program (A.A. Mironov, unpubl.). Prediction of potential transmembrane segments in protein sequences was done using TMpred (http://www.ch.embnet.org/software/TMPRED_form.html). Helix-turn-helix (HTH) DNA-binding motifs were analyzed using the weight matrix method (Dodd and Egan 1990; http://npsa-pbil.ibcp.fr/). The significance of a candidate HTH motif in a given sequence was estimated using the HTH score and probability reported by the above program. In addition, the InterPro database (Apweiler et al. 2000; http://www.ebi.ac.uk/interpro/) was used to verify the protein functional and structural annotation.

	WEB SITE REFERENCES

TOP ABSTRACT INTRODUCTION RESULTS AND DISCUSSION METHODS WEB SITE REFERENCES REFERENCES

http://ergo.integratedgenomics.com/ERGO/; ERGO database.

http://jgi.doe.gov; DOE Joint Genome Institute.

http://npsa-pbil.ibcp.fr; Network Protein Sequence Analysis server.

http://www.ch.embnet.org/software/TMPRED_form.html; TMpred Server.

http://www.ebi.ac.uk/interpro/; InterPro database.

http://www.genome.ou.edu; University of Oklahoma's Advanced Center for Genome Technology.

http://www.sanger.ac.uk; Wellcome Trust Sanger Institute.

http://www.tigr.org/; Institute for Genomic Research.