The complete genome sequence of the gram-positive bacterium Bacillus subtilis

Analogous enzymes: independent inventions in enzyme evolution

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

The Escherichia coli citrate carrier CitT: a member of a novel eubacterial transporter family related to the 2-oxoglutarate/malate translocator from spinach chloroplasts

Under anoxic conditions in the presence of an oxidizable cosubstrate such as glucose or glycerol, Escherichia coli converts citrate to acetate and succinate. Two enzymes are specifically required for the fermentation of the tricarboxylic acid, i.e., a citrate uptake system and citrate lyase. Here we report that the open reading frame (designated citT) located at 13.90 min on the E. coli chromosome between rna and the citrate lyase genes encodes a citrate carrier. E. coli transformed with a plasmid expressing citT was capable of aerobic growth on citrate, which provides convincing evidence for a function of CitT as a citrate carrier. Transport studies with cell suspensions of the transformed strain indicated that CitT catalyzes a homologous exchange of citrate or a heterologous exchange against succinate, fumarate, or tartrate. Since succinate is the end product of citrate fermentation in E. coli, it is likely that CitT functions in vivo as a citrate/succinate antiporter. Analysis of the primary sequence showed that CitT (487 amino acids, 53.1 kDa) is a highly hydrophobic protein with 12 putative transmembrane helices. Sequence comparisons revealed that CitT is related to the 2-oxoglutarate/malate translocator (SODiT1 gene product) from spinach chloroplasts and five bacterial gene products, none of which has yet been functionally characterized. It is suggested that the E. coli CitT protein is a member of a novel family of eubacterial transporters involved in the transport of di- and tricarboxylic acids.

Substitution of an alanine residue for glycine 146 in TMP kinase from Escherichia coli is responsible for bacterial hypersensitivity to bromodeoxyuridine

The wild-type TMP kinases from Escherichia coli and from a strain hypersensitive to 5-bromo-2'-deoxyuridine were characterized comparatively. The mutation at codon 146 causes the substitution of an alanine residue for glycine in the enzyme, which is accompanied by changes in the relative affinities for 5-Br-UMP and TMP compared to those of the wild-type TMP kinase. Plasmids carrying the wild-type tmk gene from Escherichia coli or Bacillus subtilis, but not the defective tmk gene, restored the resistance to bromodeoxyuridine of an E. coli mutant strain.

Identification and expression of the Bacillus subtilis fructose-1, 6-bisphosphatase gene (fbp)

The Bacillus subtilis fbp gene encoding fructose-1,6-bisphosphatase (FBPase) was originally identified as yydE. The fbp gene was expressed at a fairly constant level in cells undergoing glycolysis or gluconeogenesis. fbp transcription was initiated 94 bp upstream of the translation initiation codon, resulting in a 2.4-kb monocistronic transcript. Interestingly, B. subtilis FBPase exhibited no significant similarity to other FBPases in protein sequence databases.

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.

Cell surface localization and processing of the ComG proteins, required for DNA binding during transformation of Bacillus subtilis

The comG operon of Bacillus subtilis encodes seven proteins essential for the binding of transforming DNA to the competent cell surface. We have explored the processing of the ComG proteins and the cellular localization of six of them. All of the proteins were found to be membrane associated. The four proteins with N-terminal sequence motifs typical of type 4 pre-pilins (ComGC, GD, GE and GG) are processed by a pathway that requires the product of comC, also an essential competence gene. The unprocessed forms of ComGC and GD behave like integral membrane proteins. Pre-ComGG differs from pre-ComGC and pre-ComGD, in that it is accessible to proteolysis only from the cytoplasmic face of the membrane and at least a portion of it behaves like a peripheral membrane protein. The mature forms of these proteins are translocated to the outer face of the membrane and are liberated when peptidoglycan is hydrolysed by lysozyme or mutanolysin. ComGG exists in part as a disulphide-cross-linked homodimer in vivo. ComGC was found to possess an intramolecular disulphide bond. The previously identified homodimer form of this protein is not stabilized by disulphide bond formation. ComGF behaves as an integral membrane protein, while ComGA, a putative ATPase, is located on the inner face of the membrane as a peripheral membrane protein. Possible roles of the ComG proteins in DNA binding to the competent cell surface are discussed in the light of these and other results.

Novel approaches to the discovery of antimicrobial agents

With the completion of numerous bacterial genome sequences, the discovery of antibacterial drugs has fully entered the genomic era. The strategies for effectively using genomic information for target identification, target characterization, screen development and compound evaluation are emerging, and have greatly increased the number of antibacterial targets available for screening. Fortunately, simultaneous efforts in improving miniaturization, robotics and database tools are underway so that the potential of genomics can be realized.

Bacillus subtilis RNase III cleaves both 5'- and 3'-sites of the small cytoplasmic RNA precursor

Bacillus subtilis small cytoplasmic RNA (scRNA) is a member of the signal recognition particle RNA family. It is transcribed as a 354-nucleotide primary transcript and processed to a 271-nucleotide mature scRNA. In the precursor, the 5'- and 3'-flanking regions form a stable double-stranded structure based on their complementary sequence. This structure is similar to those of substrates for the double-stranded RNA processing enzyme, RNase III. The B. subtilis enzyme that has similar activity to Escherichia coli RNase III has been purified and is designated Bs-RNase III. Recently, B. subtilis rncS has been shown to encode Bs-RNase III (Wang, W., and Bechhofer, D. H. (1997) J. Bacteriol. 179, 7379-7385). We show here that Bs-RNase III and the purified His-tagged product of rncS cleave pre-scRNA at both 5'- and 3'-sites to produce an intermediate scRNA (scRNA-275), although processing at the 3'-site is less efficient. The 5'-end of scRNA-275 was identical to that of the mature scRNA, whereas it contains four excess nucleotides at the 3'-end. Bs-RNase III cleavage yields a two-base 3'-overhang, which is consistent with the manner in which E. coli RNase III cleaves. We also show that truncation of the rncS gene affected processing, and significant amounts of an intermediate scRNA (scRNA-275) were found to accumulate in the rncS-truncated mutant. It is concluded that Bs-RNase III is an enzyme that processes pre-scRNA.

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.

Cloning and characterization of a novel binding factor (GMEB-2) of the glucocorticoid modulatory element

The 21-base pair glucocorticoid modulatory element (GME) of the rat tyrosine aminotransferase gene is the only cis-acting element known to modulate the transcriptional activity of receptors bound to glucocorticoid response elements. Specifically, the GME increases the activity of complexes bound both by physiological concentrations of glucocorticoids, due to a left shift in the dose-response curve, and by saturating concentrations of anti-glucocorticoids. For this reason, the nuclear protein(s) that has been demonstrated to bind to the GME is of major interest as a possible transcription factor with hitherto undescribed properties. Subsequent studies indicated that not one but two proteins of 88 and 67 kDa (= GMEB-1 and -2, respectively) formed a heteromeric complex with double-stranded GME oligonucleotides in gel shift assays and participated in the expression of GME activity (Oshima, H., Szapary, D., and Simons, S. S., Jr. (1995) J. Biol. Chem. 270, 21893-21910). Here, we report the use of polymerase chain reaction of degenerate oligonucleotides and 5'- and 3'-rapid amplification of cDNA ends to clone two cDNAs of 2. 0 and 1.9 kilobase pairs that probably result from alternative splicing. Both cDNAs encoded open reading frames containing all four previously sequenced peptides. The longer 2.0-kilobase pair cDNA encoded an open reading frame for an acidic, 529-amino acid protein and afforded a major 67-kDa and a minor 58-kDa protein after in vitro transcription/translation. Both proteins were recognized by a mono-epitopic antibody raised against a peptide of GMEB-2. The in vitro translated protein bound to GME DNA in gel shift assays. However, the binding to GME DNA increased markedly after mixing with authentic GMEB-1 to give a gel-shifted complex that was similar to that derived from HTC cell cytosol. GMEB-2 shares a unique domain (KDWKR) with proteins derived from diverse organisms as follows: Drosophila (DEAF-I), rat (Suppressin), and Caenorhabditis elegans (three unknown open reading frames). Collectively, these data suggest that the 67-kDa GMEB-2 not only is an important factor for the modulation of glucocorticoid receptor bound to glucocorticoid response elements but also may belong to a novel family of transcription factors.

Replication terminator protein-based replication fork-arrest systems in various Bacillus species

The replication terminator protein (RTP) of Bacillus subtilis interacts with its cognate DNA terminators to cause replication fork arrest, thereby ensuring that the forks approaching one another at the conclusion of a round of replication meet within a restricted terminus region. A similar situation exists in Escherichia coli, but it appears that the fork-arrest systems in these two organisms have evolved independently of one another. In the present work, RTP homologs in four species closely related to B. subtilis (B. atrophaeus, B. amyloliquefaciens, B. mojavensis, and B. vallismortis) have been identified and characterized. An RTP homolog could not be identified in another closely related species, B. licheniformis. The nucleotide and amino acid changes from B. subtilis among the four homologs are consistent with the recently established phylogenetic tree for these species. The GC contents of the rtp genes raise the possibility that these organisms arose within this branch of the tree by horizontal transfer into a common ancestor after their divergence from B. licheniformis. Only 5 amino acid residue positions were changed among the four homologs, despite an up to 17.2% change in the nucleotide sequence, a finding that highlights the importance of the precise folded structure to the functioning of RTP. The absence of any significant change in the proposed DNA-binding region of RTP emphasizes the importance of its high affinity for the DNA terminator in its functioning. By coincidence, the single change (E30K) found in the B. mojavensis RTP corresponds exactly to that purposefully introduced by others into B. subtilis RTP to implicate a crucial role for E30 in the fork-arrest mechanism. The natural occurrence of this variant is difficult to reconcile with such an implication, and it was shown directly that RTP.E30K functions normally in fork arrest in B. subtilis in vivo. Additional DNA terminators were identified in the new RTP homolog-containing strains, allowing the definition of a Bacillus terminator consensus and identification of two more terminators in the B. subtilis 168 genome sequence to bring the total to nine.

Genetic recombination in Bacillus subtilis 168: effects of recU and recS mutations on DNA repair and homologous recombination

Bacillus subtilis recombination-deficient mutants were constructed by inserting a selectable marker (cat gene) into the yppB and ypbC coding regions. The yppB:cat and ypbC:cat null alleles rendered cells sensitive to DNA-damaging agents, impaired plasmid transformation (25- and 100-fold), and moderately affected chromosomal transformation when present in an otherwise Rec+ B. subtilis strain. The yppB gene complemented the defect of the recG40 strain. yppB and ypbC and their respective null alleles were termed "recU" and "recU1" (recU:cat) and "recS" and "recS1" (recS:cat), respectively. The recU and recS mutations were introduced into rec-deficient strains representative of the alpha (recF), beta (addA5 addB72), gamma (recH342), and epsilon (recG40) epistatic groups. The recU mutation did not modify the sensitivity of recH cells to DNA-damaging agents, but it did affect inter- and intramolecular recombination in recH cells. The recS mutation did not modify the sensitivity of addAB cells to DNA-damaging agents, and it marginally affected recF, recH, and recU cells. The recS mutation markedly reduced (about 250-fold) intermolecular recombination in recH cells, and there were reductions of 10- to 20-fold in recF, addAB, and recU cells. Intramolecular recombination was blocked in recS recF, recS addAB, and recS recU cells. RecU and RecS have no functional counterparts in Escherichia coli. Altogether, these data indicate that the recU and recS proteins are required for DNA repair and intramolecular recombination and that the recF (alpha epistatic group), addAB (beta), recH (gamma), recU (epsilon), and recS genes provide overlapping activities that compensate for the effects of single mutation. We tentatively placed recS within a new group, termed "zeta".

The ins and outs of peptide signaling

Many bacteria communicate by secreting and responding to extracellular peptides (pheromones). Some peptide pheromones act via receptors on the cell surface, which are often membrane-bound histidine protein kinases. Other peptide pheromones are transported into the cell by an oligopeptide permease and interact with intracellular receptors to modulate gene expression.

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.

The complete genome of Bacillus subtilis: from sequence annotation to data management and analysis

The completion of the entire 4.2-Mb genome sequence of the gram-positive bacterium Bacillus subtilis has been a milestone for biological studies on this model organism. This paper describes bioinformatics work related to this joint European and Japanese project: methods and strategies for gene annotation and detection of sequencing errors, using an integrated cooperative computer environment (Imagene); construction of a specialized database for data management and a WWW server for data retrieval (SubtiList); DNA sequence analysis, yielding striking results on oligonucleotide bias, repeated sequences, and codon usage, all landmarks of evolutionary events shaping the B. subtilis genome.

Oligonucleotide bias in Bacillus subtilis: general trends and taxonomic comparisons

We present a general analysis of oligonucleotide usage in the complete genome of Bacillus subtilis . Several datasets were built in order to assign various biological contexts to the biased use of words and to reveal local asymmetries in word usage that may be coupled with replication, the control of gene expression and the restriction/modification system. This analysis was complemented by cross-comparisons with the complete genomes of Escherichia coli , Haemophilus influenzae and Methanococcus jannaschii . We have observed a large number of biased oligonucleotides for words of size up to 8, throughout the datasets and species, indicating that such long strict words play an important role as biological signals. We speculate that some of them are involved in interactions with DNA and/or RNA polymerases. An extensive analysis of palindrome abundances and distributions provides the surprising result that prophage-like elements embedded in the genome exhibit a smaller avoidance of restriction sites. This may reinforce a recently proposed hypothesis of a selfish gene phenomena in the transfer of restriction/modification systems in bacteria.

Combining diverse evidence for gene recognition in completely sequenced bacterial genomes

Analysis of a newly sequenced bacterial genome starts with identification of protein-coding genes. Functional assignment of proteins requires the exact knowledge of protein N-termini. We present a new program ORPHEUS that identifies candidate genes and accurately predicts gene starts. The analysis starts with a database similarity search and identification of reliable gene fragments. The latter are used to derive statistical characteristics of protein-coding regions and ribosome-binding sites and to predict the complete set of genes in the analyzed genome. In a test on Bacillus subtilis and Escherichia coli genomes, the program correctly identified 93.3% (resp. 96.3%) of experimentally annotated genes longer than 100 codons described in the PIR-International database, and for these genes 96.3% (83.9%) of starts were predicted exactly. Furthermore, 98.9% (99.1%) of genes longer than 100 codons annotated in GenBank were found, and 92.9% (75.7%) of predicted starts coincided with the feature table description. Finally, for the complete gene complements of B.subtilis and E.coli , including genes shorter than 100 codons, gene prediction accuracy was 88.9 and 87.1%, respectively, with 94.2 and 76.7% starts coinciding with the existing annotation.

Mammalian mitochondrial methionyl-tRNA transformylase from bovine liver. Purification, characterization, and gene structure

The mammalian mitochondrial methionyl-tRNA transformylase (MTFmt) was partially purified 2,200-fold from bovine liver mitochondria using column chromatography. The polypeptide responsible for MTFmt activity was excised from a sodium dodecyl sulfate-polyacrylamide gel and the amino acid sequences of several peptides were determined. The cDNA encoding bovine MTFmt was obtained and its nucleotide sequence was determined. The deduced amino acid sequence of the mature form of MTFmt consists of 357 amino acid residues. This sequence is about 30% identical to the corresponding Escherichia coli and yeast mitochondrial MTFs. Kinetic parameters governing the formylation of various tRNAs were obtained. Bovine MTFmt formylates its homologous mitochondrial methionyl-tRNA and the E. coli initiator methionyl-tRNA (Met-tRNAfMet) with essentially equal efficiency. The E. coli elongator methionyl-tRNA (Met-tRNAmMet) was also formylated although with somewhat less favorable kinetics. These results suggest that the substrate specificity of MTFmt is not as rigid as that of the E. coli MTF which clearly discriminates between the bacterial initiator and elongator Met-tRNAs. These observations are discussed in terms of the presence of a single tRNAMet gene in mammalian mitochondria.

Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence

Countless millions of people have died from tuberculosis, a chronic infectious disease caused by the tubercle bacillus. The complete genome sequence of the best-characterized strain of Mycobacterium tuberculosis, H37Rv, has been determined and analysed in order to improve our understanding of the biology of this slow-growing pathogen and to help the conception of new prophylactic and therapeutic interventions. The genome comprises 4,411,529 base pairs, contains around 4,000 genes, and has a very high guanine + cytosine content that is reflected in the biased amino-acid content of the proteins. M. tuberculosis differs radically from other bacteria in that a very large portion of its coding capacity is devoted to the production of enzymes involved in lipogenesis and lipolysis, and to two new families of glycine-rich proteins with a repetitive structure that may represent a source of antigenic variation.

The ABC family of multidrug transporters in microorganisms

Multidrug transporters are membrane proteins that are able to expel a broad range of toxic molecules from the cell. In humans, the overexpression of the multidrug resistance P-glycoprotein (Pgp) and the multidrug resistance-associated protein MRP1 (MRP) is a principal cause of resistance of cancers to chemotherapy. These multidrug transporters belong to the ATP-binding cassette (ABC) family of transport proteins that utilize the energy of ATP hydrolysis for activity. In microorganisms, multidrug transporters play an important role in conferring antibiotic resistance on pathogens. In the last decade, homologs of human Pgp and MRP have been found in microorganisms such as Plasmodium falciparum, Candida albicans, Saccharomyces cerevisiae and, more recently, in Lactococcus lactis. In this review, we will summarize the current state of knowledge on three major aspects of microbial ABC-type multidrug transporters: (i) the functional and structural similarities among these proteins in prokaryotic and eukaryotic cells, (ii) the molecular mechanism of these transporters, and (iii) their potential physiological role.

Characterization of yhcN, a new forespore-specific gene of Bacillus subtilis

A new Bacillus subtilis sporulation-specific gene, yhcN, has been identified, the expression of which is dependent on the forespore-specific sigma factor sigmaG and to a much lesser extent on sigmaF. A translational yhcN-lacZ fusion is expressed at a very high level in the forespore, and the protein encoded by yhcN was detected in the inner spore membrane. A yhcN mutant sporulates normally and yhcN spores have identical resistance properties to wild-type spores. However, the outgrowth of yhcN spores is slower than that of wild-type spores.

The Bacillus stearothermophilus argCJBD operon harbours a strong promoter as evaluated in Escherichia coli cells

We have shown that the B. stearothermophilus argCJBD genes form a single operon. In B. stearothermophilus, a specific repressor governs operon expression by binding to the argCo operator site overlapping the Parg promoter sequence (Dion et al., 1997). Therefore, the enzymatic and transcriptional analyses performed in this work did not reflect the potential strength of Parg in the native host. For evaluation of the Parg promoter strength, E. coli was used as a host since its own ArgR repressor does not interact with the B. stearothermophilus heterologous operator. Parg-promoted argC gene expression dramatically increased, reaching up to 38% of the total protein in E. coli cells. An AT-rich sequence upstream of a -35 site of Parg was found to be indispensable for the promoter strength. Plasmids carrying the B. stearothermophilus argCJBD operon linked with its Parg/argCo region were unstable in E. coli. Stabilization of plasmids was achieved by repression of B. stearothermophilus arg genes through the action of the B. subtilis AhrC repressor.

SMC protein complexes and higher-order chromosome dynamics

The structural maintenance of chromosome (SMC) family of proteins represents an expanding group of chromosomal ATPases that are highly conserved among Bacteria, Archaea and Eukarya. During the past year, significant progress has been made towards understanding the cellular functions and molecular activities of this new class of proteins. Emerging evidence suggests that eukaryotic SMC proteins form large protein complexes with non-SMC subunits and act as key components for a wide variety of higher-order chromosome dynamics.

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.

The yvsA-yvqA (293 degrees-289 degrees) region of the Bacillus subtilis chromosome containing genes involved in metal ion uptake and a putative sigma factor

The region between yvsA (293 degrees) and yvqA (289 degrees) of the Bacillus subtilis chromosome has been sequenced within the framework of the B. subtilis 168 international sequencing programme. A primary analysis of the 42 ORFs identified in this 43 kb region is presented. The region included a high proportion of genes that did not show homology with genes in other bacteria. The identified ORFs showed homology to proteins involved in the transport of metal ions, two-component signal transducers, ATP-binding-cassette-type transporters and a sigma factor.

The glucose kinase of Bacillus subtilis

The open reading frame yqgR (now termed glcK), which had been sequenced as part of the genome project, encodes a glucose kinase of Bacillus subtilis. A 1.1-kb DNA fragment containing glcK complemented an Escherichia coli strain deficient in glucose kinase activity. Insertional mutagenesis of glcK resulted in a complete inactivation of glucose kinase activity in crude protein extracts, indicating that B. subtilis contains one major glucose kinase. The glcK gene encodes a 321-residue protein with a molecular mass of 33.5 kDa. The glucose kinase was overexpressed as a fusion protein to a six-His affinity tag and purified to homogeneity. The enzyme had K(m) values for ATP and glucose of 0.77 and 0.24 mM, respectively, and a Vmax of 93 mumol min-1 mg-1. A B. subtilis strain deficient for glucose kinase grew at the same rate on different carbon sources tested, including disaccharides such as maltose, trehalose, and sucrose.

Systems for categorizing functions of gene products

As the sequencing of the total DNA of many organisms continues, attention is turning next to the interpretation of the function of all of the genes and gene products, with the aim of learning the full meaning of the entire genetic blueprint of a sequenced organism in concrete terms. To set the stage for accomplishing this, systems need to be constructed for the expression of the functions of gene products in systematic yet rich ways. As much as possible, the same systems should be applicable to all organisms, so that when comparability exists among organisms the connections will become clear.

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.

PRD--a protein domain involved in PTS-dependent induction and carbon catabolite repression of catabolic operons in bacteria

Several operon-specific transcriptional regulators, including antiterminators and activators, contain a duplicated conserved domain, the PTS regulation domain (PRD). These duplicated domains modify the activity of the transcriptional regulators both positively and negatively. PRD-containing regulators are very common in Gram-positive bacteria. In contrast, antiterminators controlling beta-glucoside utilization are the only functionally characterized members of this family from gram-negative bacteria. PRD-containing regulators are controlled by PTS-dependent phosphorylation with different consequences: (i) In the absence of inducer, the phosphorylated EIIB component of the sugar permease donates its phosphate to a PRD, thereby inactivating the regulator. In the presence of the substrate, the regulator is dephosphorylated, and the phosphate is transferred to the sugar, resulting in induction of the operon. (ii) In gram-positive bacteria, a novel mechanism of carbon catabolite repression mediated by PRD-containing regulators has been demonstrated. In the absence of PTS substrates, the HPr protein is phosphorylated by enzyme I at His-15. This form of HPr can, in turn, phosphorylate PRD-containing regulators and stimulate their activity. In the presence of rapidly metabolizable carbon sources, ATP-dependent phosphorylation of HPr at Ser-46 by HPr kinase inhibits phosphorylation by enzyme I, and PRD-containing regulators cannot, therefore, be stimulated and are inactive. All regulators of this family contain two copies of PRD, which are functionally specialized in either induction or catabolite repression.

A spore coat protein, CotS, of Bacillus subtilis is synthesized under the regulation of sigmaK and GerE during development and is located in the inner coat layer of spores

The spore coat of Bacillus subtilis has a unique morphology and consists of polypeptides of different sizes, whose synthesis and assembly are precisely regulated by a cascade of transcription factors and regulatory proteins. We examined the factors that regulate cotS gene expression and CotS assembly into the coat layer of B. subtilis by Northern blot and Western blot analysis. Transcription of cotS mRNA was not detected in sporulating cells of sigmaK and gerE mutants by Northern blot analysis. By Western blot analysis using anti-CotS antibody, CotS was first detected in protein samples solubilized from wild-type cells at 5 h after the start of sporulation. CotS was not detected in the vegetative cells and spores of a gerE mutant or in the spores of mutants deficient in sigmaE, sigmaF, sigmaG, or sigmaK. CotS was detected in the sporangium but not in the spores of a cotE mutant. The sequence of the promoter region of cotS was similar to the consensus sequences for binding of sigmaK and GerE. These results demonstrate that sigmaK and GerE are required for cotS expression and that CotE is essential for the assembly of CotS in the coat. Immunoelectron microscopic observation using anti-CotS antibody revealed that CotS is located within the spore coat, in particular in the inner coats of dormant spores.

HrpW of Erwinia amylovora, a new Hrp-secreted protein

Erwinia amylovora strain CFBP1430 secretes a protein called HrpW in a Hrp-dependent manner. HrpW was detected in culture supernatant of the wild-type strain grown on solid inducing hrp medium. This protein shares structural similarities with elicitors of the hypersensitive response such as HrpN of Erwinia amylovora and PopA of Ralstonia solanacearum. Furthermore, the C-terminal region of HrpW is homologous to class III pectate lyases. An hrpW mutant is as aggressive as the wild-type strain on pear and apple seedlings. It elicits the hypersensitive response on tobacco at a lower concentration than the wild-type strain.

Identification of target promoters for the Bacillus subtilis sigma X factor using a consensus-directed search

The promoter selectivity of RNA polymerase (RNAP) can be altered by the association with alternative sigma subunits. Bacillus subtilis hosts a multitude of sigma factors, several of which coordinate the complex developmental program culminating in endospore formation. Genome sequencing has revealed an unanticipated seven new sigma factors of the highly divergent extracytoplasmic function (ECF) sub-family. Virtually nothing is known regarding either the promoter selectivity or the target genes for these newly identified sigma factors. We have used saturation mutagenesis to define a promoter consensus for recognition by one such ECF sigma factor, sigma X. The resulting consensus sequence was used to identify candidate sigma X target sites. Three newly identified sigma X-dependent promoters precede genes encoding regulatory proteins: an AbrB homolog (Abh), a putative response regulator aspartate phosphatase (RapD), and a regulator of autolysin expression (LytR). sigma X also contributes to the expression of CsbB, a putative membrane-bound glucosyl transferase that is partially controlled by the sigma B stress response sigma factor. Since LytR modulates the expression of the major autolytic amidase and CsbB may function in peptidoglycan synthesis or modification, we suggest that sigma X participates in the regulation of peptidoglycan synthesis and turnover.

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.

Sequence specificity of illegitimate plasmid recombination in Bacillus subtilis: possible recognition sites for DNA topoisomerase I

Previous work in our group indicated that structural plasmid instability in Bacillus subtilis is often caused by illegitimate recombination between non-repeated sequences, characterized by a relatively high AT content. Recently we developed a positive selection vector for analysis of plasmid recombination events in B. subtilis which enables measurement of recombination frequencies without interference of selective growth differences of cells carrying wild-type or deleted plasmids. Here we have used this system to further analyse the sequence specificity of illegitimate plasmid recombination events and to assess the role of the host-encoded DNA topoisomerase I enzyme in this process. Several lines of evidence suggest that single-strand DNA nicks introduced by DNA topoisomerase I are a major source of plasmid deletions in pGP100. First, strains overproducing DNA topoisomerase I showed increased levels of plasmid deletion. Second, these deletions occurred predominantly (>90% of the recombinants) between non-repeated DNA sequences, the majority of which resemble potential DNA topoisomerase I target sites. Sequence alignment of 66 deletion end-points confirmed the previously reported high AT content and, most importantly, revealed a highly conserved C residue at position -4 relative to the site of cleavage at both deletion termini. Based on these genetic data we propose the following putative consensus cleavage site for DNA topoisomerase I of B.subtilis: 5'-A/TCATA/TTAA/TA/TA-3'.

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.

The cbaAB genes for bo3-type cytochrome c oxidase in Bacillus stearothermophilus

Structural genes were cloned for cytochrome bo3-type cytochrome c oxidase recently isolated from a Gram-positive thermophile Bacillus stearothermophilus. Sequencing and Northern blotting analyses indicated that the two genes cbaA and cbaB composed an operon encoding for subunits I and II, respectively, and that the oxidase was SoxB-type. They are the first genes for a SoxB-type cytochrome c oxidase whose natural substrate is known.

Identification and enzymatic characterization of the maltose-inducible alpha-glucosidase MalL (sucrase-isomaltase-maltase) of Bacillus subtilis

A gene coding for a putative alpha-glucosidase has been identified in the open reading frame yvdL (now termed malL), which was sequenced as part of the Bacillus subtilis genome project. The enzyme was overproduced in Escherichia coli and purified. Further analyses indicate that MalL is a specific oligo-1,4-1,6-alpha-glucosidase (sucrase-maltase-isomaltase). MalL expression in B. subtilis requires maltose induction and is subject to carbon catabolite repression by glucose and fructose. Insertional mutagenesis of malL resulted in a complete inactivation of the maltose-inducible alpha-glucosidase activity in crude protein extracts and a Mal- phenotype.

Integrated physical and genetic mapping of Bacillus cereus and other gram-positive bacteria based on IS231A transposition vectors

The genome structure of Bacillus cereus is relatively complex, its DNA being modulated between a size-varying chromosome and large plasmids. To study the genetic organization of the B. cereus type strain ATCC 14579, thermosensitive transposition vectors were designed on the basis of IS231A-derived cassettes containing uncommon restriction sites. A highly preferred insertion site for IS231A was detected in the chromosome by Southern blotting and pulsed-field gel electrophoresis (PFGE) analyses of independent insertion mutants. However, once this insertional hot spot was occupied, secondary IS231A insertions occurred randomly, as demonstrated by isolation of independent B. cereus auxotrophs at a frequency of approximately 0.6%. The hot-spot site, as well as several auxotrophic mutations, were mapped by using NotI, SfiI, and AscI PFGE restriction profiles. It was confirmed by sequencing that one of the insertions, generating an Ade- phenotype, had disrupted a gene of the purine synthesis pathway. These results showed that combined PFGE and sequencing analyses of mini-IS231A insertions enable the construction of integrated physical and genetic maps of B. cereus type strain. Moreover, the presence of the ultrarare I-SceI restriction site in the mini-IS231A allowed the isolation, in double-insertion mutants, of contiguous and nonoverlapping large chromosomal fragments, convenient for direct sequencing. The system detailed in this report is therefore a powerful tool for comparative genetic studies among members of the B. cereus group (i.e., B. cereus, B. thuringiensis, B. mycoides, and B. anthracis) and could also be applied to more distantly related gram-positive bacteria.

Gas vesicle genes identified in Bacillus megaterium and functional expression in Escherichia coli

Gas vesicles are intracellular, protein-coated, and hollow organelles found in cyanobacteria and halophilic archaea. They are permeable to ambient gases by diffusion and provide buoyancy, enabling cells to move upwards in liquid to access oxygen and/or light. In halobacteria, gas vesicle production is encoded in a 9-kb cluster of 14 genes (4 of known function). In cyanobacteria, the number of genes involved has not been determined. We now report the cloning and sequence analysis of an 8,142-bp cluster of 15 putative gas vesicle genes (gvp) from Bacillus megaterium VT1660 and their functional expression in Escherichia coli. Evidence includes homologies by sequence analysis to known gas vesicle genes, the buoyancy phenotype of E. coli strains that carry this gvp gene cluster, the presence of pressure-sensitive, refractile bodies in phase-contrast microscopy, structural details in phase-contrast microscopy, structural details in direct interference-contrast microscopy, and shape and size revealed by transmission electron microscopy. In B. megaterium, the gvp region carries a cluster of 15 putative genes arranged in one orientation; they are open reading frame 1 and gvpA, -P, -Q, -B, -R, -N, -F, -G, -L, -S, -K, -J, -T, and -U, of which the last 11 genes, in a 5.7-kb gene cluster, are the maximum required for gas vesicle synthesis and function in E. coli. To our knowledge, this is the first example of a functional gas vesicle gene cluster in nonaquatic bacteria and the first example of the interspecies transfer of genes resulting in the synthesis of a functional organelle.

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.

Antibacterial agents that inhibit two-component signal transduction systems

A class of antibacterials has been discovered that inhibits the growth of Gram-positive pathogenic bacteria. RWJ-49815, a representative of a family of hydrophobic tyramines, in addition to being a potent bactericidal Gram-positive antibacterial, inhibits the autophosphorylation of kinase A of the KinA::Spo0F two-component signal transduction system in vitro. Analogs of RWJ-49815 vary greatly in their ability to inhibit growth of bacteria and this ability correlates directly with their activity as kinase A inhibitors. Compared with the potent quinolone, ciprofloxacin, RWJ-49815 exhibits reduced resistance emergence in a laboratory passage experiment. Inhibition of the histidine protein kinase::response regulator two-component signal transduction pathways may present an opportunity to depress chromosomal resistance emergence by targeting multiple proteins with a single inhibitor in a single bacterium. Such inhibitors may represent a class of antibacterials that potentially may represent a breakthrough in antibacterial therapy.

Bipolar localization of Bacillus subtilis topoisomerase IV, an enzyme required for chromosome segregation

In Bacillus subtilis, parE and parC were shown to be essential genes for the segregation of replicated chromosomes. Disruption of either one of these genes resulted in failure of the nucleoid to segregate. Purified ParE and ParC proteins reconstituted to form topoisomerase IV (topo IV), which was highly proficient for ATP-dependent superhelical DNA relaxation and decatenation of interlocked DNA networks. By immunofluorescence microscopy and by directly visualizing fluorescence by using green fluorescence protein fusions, we determined that ParC is localized at the poles of the bacteria in rapidly growing cultures. The bipolar localization of ParC required functional ParE, suggesting that topo IV activity is required for the localization. ParE was found to be distributed uniformly throughout the cell. On the other hand, fluorescence microscopy showed that the GyrA and GyrB subunits of gyrase were associated with the nucleoid. Our results provide a physiologic distinction between DNA gyrase and topo IV. The subcellular localization of topo IV provides physical evidence that it may be part of the bacterial segregation machinery.

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.

Microbial genome analyses: global comparisons of transport capabilities based on phylogenies, bioenergetics and substrate specificities

We have conducted genome sequence analyses of seven prokaryotic microorganisms for which completely sequenced genomes are available (Escherichia coli, Haemophilus influenzae, Helicobacter pylori, Bacillus subtilis, Mycoplasma genitalium, Synechocystis PCC6803 and Methanococcus jannaschii). We report the distribution of encoded known and putative polytopic cytoplasmic membrane transport proteins within these genomes. Transport systems for each organism were classified according to (1) putative membrane topology, (2) protein family, (3) bioenergetics, and (4) substrate specificities. The overall transport capabilities of each organism were thereby estimated. Probable function was assigned to greater than 90% of the putative transport proteins identified. The results show the following: (1) Numbers of transport systems in eubacteria are approximately proportional to genome size and correspond to 9.7 to 10.8% of the total encoded genes except for H. pylori (5.4%), Synechocystis (4.7%) and M. jannaschii (3.5%) which exhibit substantially lower proportions. (2) The distribution of topological types is similar in all seven organisms. (3) Transport systems belonging to 67 families were identified within the genomes of these organisms, and about half of these families are also found in eukaryotes. (4) 12% of these families are found exclusively in Gram-negative bacteria, but none is found exclusively in Gram-positive bacteria, cyanobacteria or archaea. (5) Two superfamilies, the ATP-binding cassette (ABC) and major facilitator (MF) superfamilies account for nearly 50% of all transporters in each organism, but the relative representation of these two transporter types varies over a tenfold range, depending on the organism. (6) Secondary, pmf-dependent carriers are 1.5 to threefold more prevalent than primary ATP-dependent carriers in E. coli, H. influenzae, H. pylori and B. subtilis while primary carriers are about twofold more prevalent in M. genitalium and Synechocystis. M. jannaschii exhibits a slight preference for secondary carriers. (7) Bioenergetics of transport generally correlate with the primary forms of energy generated via available metabolic pathways but ecological niche and substrate availability may also be determining factors. (8) All organisms display a similar range of transport specificities with quantitative differences presumably reflective of disparate ecological niches. (9) M. jannaschii and Synechocystis have a two to threefold increased proportion of transporters for inorganic ions with a concomitant decrease in transporters for organic compounds. (10) 6 to 18% of all transporters in these bacteria probably function as drug export systems showing that these systems are prevalent in non-pathogenic as well as pathogenic organisms. (11) All seven prokaryotes examined encode proteins homologous to known channel proteins, but none of the channel types identified occurs in all of these organisms. (12) The phosphoenolpyruvate:sugar phosphotransferase system is prevalent in the large genome organisms, E. coli and B. subtilis, and is present in the small genome organisms, H. influenzae and M. genitalium, but is totally lacking in H. pylori, Synechocystis and M. jannaschii. Details of the information summarized in this article are available on our web sites, and this information will be periodically updated and corrected as new sequence and biochemical data become available.

A Bacillus subtilis gene induced by cold shock encodes a membrane phospholipid desaturase

Bacillus subtilis grown at 37 degrees C synthesizes saturated fatty acids with only traces of unsaturated fatty acids (UFAs). However, when cultures growing at 37 degrees C are transferred to 20 degrees C, UFA synthesis is induced. We report the identification and characterization of the gene encoding the fatty acid desaturase of B. subtilis. This gene, called des, was isolated by complementation of Escherichia coli strains with mutations in either of two different genes of UFA synthesis. The des gene encodes a polypeptide of 352 amino acid residues containing the three conserved histidine cluster motifs and two putative membrane-spanning domains characteristic of the membrane-bound desaturases of plants and cyanobacteria. Expression of the des gene in E. coli resulted in desaturation of palmitic acid moieties of the membrane phospholipids to give the novel mono-UFA cis-5-hexadecenoic acid, indicating that the B. subtilis des gene product is a delta5 acyl-lipid desaturase. The des gene was disrupted, and the resulting null mutant strains were unable to synthesize UFAs upon a shift to low growth temperatures. The des null mutant strain grew as well as its congenic parent at 20 or 37 degrees C but showed severely reduced survival during stationary phase. Analysis of operon fusions in which the des promoter directed the synthesis of a lacZ reporter gene showed that des expression is repressed at 37 degrees C, but a shift of cultures from 37 to 20 degrees C resulted in a 10- to 15-fold increase in transcription. This is the first report of a membrane phospholipid desaturase in a nonphotosynthetic organism and the first direct evidence for cold induction of a desaturase.

The Bacillus subtilis galE gene is essential in the presence of glucose and galactose

Bacillus subtilis is unable to grow by consuming galactose because it is unable to transport it into the cell. The transcription of galE is not influenced by galactose but is repressed by glucose. Galactose is toxic for galE-negative bacteria because it results in elevated levels of metabolic intermediates. These negative effects are reduced in galK and galT mutants. Glucose is also toxic for galE-negative strains.

Genome-wide analysis of integral membrane proteins from eubacterial, archaean, and eukaryotic organisms

We have carried out detailed statistical analyses of integral membrane proteins of the helix-bundle class from eubacterial, archaean, and eukaryotic organisms for which genome-wide sequence data are available. Twenty to 30% of all ORFs are predicted to encode membrane proteins, with the larger genomes containing a higher fraction than the smaller ones. Although there is a general tendency that proteins with a smaller number of transmembrane segments are more prevalent than those with many, uni-cellular organisms appear to prefer proteins with 6 and 12 transmembrane segments, whereas Caenorhabditis elegans and Homo sapiens have a slight preference for proteins with seven transmembrane segments. In all organisms, there is a tendency that membrane proteins either have many transmembrane segments with short connecting loops or few transmembrane segments with large extra-membraneous domains. Membrane proteins from all organisms studied, except possibly the archaeon Methanococcus jannaschii, follow the so-called "positive-inside" rule; i.e., they tend to have a higher frequency of positively charged residues in cytoplasmic than in extra-cytoplasmic segments.

The 172 kb prkA-addAB region from 83 degrees to 97 degrees of the Bacillus subtilis chromosome contains several dysfunctional genes, the glyB marker, many genes encoding transporter proteins, and the ubiquitous hit gene

A 171812 bp nucleotide sequence between prkA and addAB (83 degrees to 97 degrees) on the genetic map of the Bacillus subtilis 168 chromosome was determined and analysed. An accurate physical/genetic map of this previously poorly described chromosomal region was constructed. One hundred and seventy open reading frames (ORFs) were identified on the DNA fragment. These include the previously described genes cspB, glpPFKD, spoVR, phoAIV, papQ, citRA, sspB, prsA, hpr, pbpF, hemEHY, aprE, comK and addAB. ORF yhaF in this region corresponds to the glyB marker. Among the striking features of this region are: an abundance of genes encoding (putative) transporter proteins, several dysfunctional genes, the ubiquitous hit gene, and five multidrug-resistance-like genes. These analyses have also revealed the existence of numerous paralogues of ORFs in this region: about two-thirds of the putative genes seem to have at least one paralogue in the B. subtilis genome.