Information analysis of Fis binding sites

Differential expression of two bc1 complexes in the strict acidophilic chemolithoautotrophic bacterium Acidithiobacillus ferrooxidans suggests a model for their respective roles in iron or sulfur oxidation

Three strains of the strict acidophilic chemolithoautotrophic Acidithiobacillus ferrooxidans, including the type strain ATCC 23270, contain a petIIABC gene cluster that encodes the three proteins, cytochrome c1, cytochrome b and a Rieske protein, that constitute a bc1 electron-transfer complex. RT-PCR and Northern blotting show that the petIIABC cluster is co-transcribed with cycA, encoding a cytochrome c belonging to the c4 family, sdrA, encoding a putative short-chain dehydrogenase, and hip, encoding a high potential iron-sulfur protein, suggesting that the six genes constitute an operon, termed the petII operon. Previous results indicated that A. ferrooxidans contains a second pet operon, termed the petI operon, which contains a gene cluster that is similarly organized except that it lacks hip. Real-time PCR and Northern blot experiments demonstrate that petI is transcribed mainly in cells grown in medium containing iron, whereas petII is transcribed in cells grown in media containing sulfur or iron. Primer extension experiments revealed possible transcription initiation sites for the petI and petII operons. A model is presented in which petI is proposed to encode the bc1 complex, functioning in the uphill flow of electrons from iron to NAD(P), whereas petII is suggested to be involved in electron transfer from sulfur (or formate) to oxygen (or ferric iron). A. ferrooxidans is the only organism, to date, to exhibit two functional bc1 complexes.

Anatomy of Escherichia coli sigma70 promoters

Information theory was used to build a promoter model that accounts for the -10, the -35 and the uncertainty of the gap between them on a common scale. Helical face assignment indicated that base -7, rather than -11, of the -10 may be flipping to initiate transcription. We found that the sequence conservation of sigma70 binding sites is 6.5 +/- 0.1 bits. Some promoters lack a -35 region, but have a 6.7 +/- 0.2 bit extended -10, almost the same information as the bipartite promoter. These results and similarities between the contacts in the extended -10 binding and the -35 suggest that the flexible bipartite sigma factor evolved from a simpler polymerase. Binding predicted by the bipartite model is enriched around 35 bases upstream of the translational start. This distance is the smallest 5' mRNA leader necessary for ribosome binding, suggesting that selective pressure minimizes transcript length. The promoter model was combined with models of the transcription factors Fur and Lrp to locate new promoters, to quantify promoter strengths, and to predict activation and repression. Finally, the DNA-bending proteins Fis, H-NS and IHF frequently have sites within one DNA persistence length from the -35, so bending allows distal activators to reach the polymerase.

Mechanism of chromosome compaction and looping by the Escherichia coli nucleoid protein Fis

Fis, the most abundant DNA-binding protein in Escherichia coli during rapid growth, has been suspected to play an important role in defining nucleoid structure. Using bulk-phase and single-DNA molecule experiments, we analyze the structural consequences of non-specific binding by Fis to DNA. Fis binds DNA in a largely sequence-neutral fashion at nanomolar concentrations, resulting in mild compaction under applied force due to DNA bending. With increasing concentration, Fis first coats DNA to form an ordered array with one Fis dimer bound per 21 bp and then abruptly shifts to forming a higher-order Fis-DNA filament, referred to as a low-mobility complex (LMC). The LMC initially contains two Fis dimers per 21 bp of DNA, but additional Fis dimers assemble into the LMC as the concentration is increased further. These complexes, formed at or above 1 microM Fis, are able to collapse large DNA molecules via stabilization of DNA loops. The opening and closing of loops on single DNA molecules can be followed in real time as abrupt jumps in DNA extension. Formation of loop-stabilizing complexes is sensitive to high ionic strength, even under conditions where DNA bending-compaction is unaltered. Analyses of mutants indicate that Fis-mediated DNA looping does not involve tertiary or quaternary changes in the Fis dimer structure but that a number of surface-exposed residues located both within and outside the helix-turn-helix DNA-binding region are critical. These results suggest that Fis may play a role in vivo as a domain barrier element by organizing DNA loops within the E. coli chromosome.

Dynamic state of DNA topology is essential for genome condensation in bacteria

In bacteria, Dps is one of the critical proteins to build up a condensed nucleoid in response to the environmental stresses. In this study, we found that the expression of Dps and the nucleoid condensation was not simply correlated in Escherichia coli, and that Fis, which is an E. coli (gamma-Proteobacteria)-specific nucleoid protein, interfered with the Dps-dependent nucleoid condensation. Atomic force microscopy and Northern blot analyses indicated that the inhibitory effect of Fis was due to the repression of the expression of Topoismerase I (Topo I) and DNA gyrase. In the Deltafis strain, both topA and gyrA/B genes were found to be upregulated. Overexpression of Topo I and DNA gyrase enhanced the nucleoid condensation in the presence of Dps. DNA-topology assays using the cell extract showed that the extracts from the Deltafis and Topo I-/DNA gyrase-overexpressing strains, but not the wild-type extract, shifted the population toward relaxed forms. These results indicate that the topology of DNA is dynamically transmutable and that the topology control is important for Dps-induced nucleoid condensation.

Roles for DNA supercoiling and the Fis protein in modulating expression of virulence genes during intracellular growth of Salmonella enterica serovar Typhimurium

Adaptation of bacterial pathogens to an intracellular environment requires resetting of the expression levels of a wide range of both virulence and housekeeping genes. We investigated the possibility that changes in DNA supercoiling could modulate the expression of genes known to be important in the intracellular growth of the pathogen Salmonella enterica serovar Typhimurium. Our data show that DNA becomes relaxed when Salmonella grows in murine macrophage but not in epithelial cells, indicating that DNA supercoiling plays a role in discrimination between two types of intracellular environment. The ssrA regulatory gene within the SPI-2 pathogenicity island that is required for survival in macrophage was found to be upregulated by DNA relaxation. This enhancement of expression also required the Fis nucleoid-associated protein. Manipulating the level of the Fis protein modulated both the level of DNA supercoiling and ssrA transcription. We discuss a model of bacterial intracellular adaptation in which Fis and DNA supercoiling collaborate to fine-tune virulence gene expression.

Recovering motifs from biased genomes: application of signal correction

A significant problem in biological motif analysis arises when the background symbol distribution is biased (e.g. high/low GC content in the case of DNA sequences). This can lead to overestimation of the amount of information encoded in a motif. A motif can be depicted as a signal using information theory (IT). We apply two concepts from IT, distortion and patterned interference (a type of noise), to model genomic and codon bias respectively. This modeling approach allows us to correct a raw signal to recover signals that are weakened by compositional bias. The corrected signal is more likely to be discriminated from a biased background by a macromolecule. We apply this correction technique to recover ribosome-binding site (RBS) signals from available sequenced and annotated prokaryotic genomes having diverse compositional biases. We observed that linear correction was sufficient for recovering signals even at the extremes of these biases. Further comparative genomics studies were made possible upon correction of these signals. We find that the average Euclidian distance between RBS signal frequency matrices of different genomes can be significantly reduced by using the correction technique. Within this reduced average distance, we can find examples of class-specific RBS signals. Our results have implications for motif-based prediction, particularly with regards to the estimation of reliable inter-genomic model parameters.

CorreLogo: an online server for 3D sequence logos of RNA and DNA alignments

We present an online server that generates a 3D representation of properties of user-submitted RNA or DNA alignments. The visualized properties are information of single alignment columns, mutual information of two alignment positions as well as the position-specific fraction of gaps. The nucleotide composition of both single columns and column pairs is visualized with the help of color-coded 3D bars labeled with letters. The server generates both VRML and JVX output that can be viewed with a VRML viewer or the JavaView applet, respectively. We show that combining these different features of an alignment into one 3D representation is helpful in identifying correlations between bases and potential RNA and DNA base pairs. Significant known correlations between the tRNA 3' anticodon cardinal nucleotide and the extended anticodon were observed, as were correlations within the amino acid acceptor stem and between the cardinal nucleotide and the acceptor stem. The online server can be accessed using the URL http://correlogo.abcc.ncifcrf.gov.

Fis is required for proper regulation of ssaG expression in Salmonella enterica serovar Typhimurium

Salmonella pathogenicity island 2 (SPI2) encodes a type III secretion system (TTSS) necessary for bacterial survival and replication in intracellular environment of host cells. SPI2 genes are transcribed preferentially after Salmonella enters the host cells. Transcriptional regulation of ssaG encoding the component of SPI2-TTSS apparatus was studied in vivo and in vitro. Fis, one of the major components of bacterial nucleoid, activated the stationary phase-specific expression of ssaG when Salmonella was grown in LB media. Gel-shift and footprinting analysis showed Fis bound to four distinct sites of the ssaG promoter region with different affinities. All four Fis-binding sites were required for timely transcription activation of ssaG after Salmonella entered macrophage cells. Gentamicin protection experiments using bacteria grown to stationary phase prior to infection showed that the ability of the fis mutant strain to replicate within the RAW264.7 macrophage cells was lower than the wild type. These observations confirm that Fis plays an important role in regulations of SPI2 as well as SPI1 for an efficient regulation of the virulence genes.

Common and variable contributions of Fis residues to high-affinity binding at different DNA sequences

Fis is a nucleoid-associated protein that interacts with poorly related DNA sequences with a high degree of specificity. A difference of more than 3 orders of magnitude in apparent Kd values was observed between specific (Kd, approximately 1 to 4 nM) and nonspecific (Kd, approximately 4 microM) DNA binding. To examine the contributions of Fis residues to the high-affinity binding at different DNA sequences, 13 alanine substitutions were generated in or near the Fis helix-turn-helix DNA binding motif, and the resulting proteins were purified. In vitro binding assays at three different Fis sites (fis P II, hin distal, and lambda attR) revealed that R85, T87, R89, K90, and K91 played major roles in high-affinity DNA binding and that R85, T87, and K90 were consistently vital for binding to all three sites. Other residues made variable contributions to binding, depending on the binding site. N84 was required only for binding to the lambda attR Fis site, and the role of R89 was dramatically altered by the lambda attR DNA flanking sequence. The effects of Fis mutations on fis P II or hin distal site binding in vitro generally correlated with their abilities to mediate fis P repression or DNA inversion in vivo, demonstrating that the in vitro DNA-binding effects are relevant in vivo. The results suggest that while Fis is able to recognize a minimal common set of DNA sequence determinants at different binding sites, it is also equipped with a number of residues that contribute to the binding strength, some of which play variable roles.

Repressor of Phage 16 - 3 with Altered Binding Specificity Indicates Spatial Differences in Repressor-Operator Complexes

The C repressor protein of phage 16-3 , which is required for establishing and maintaining lysogeny, recognizes structurally different operators which differ by 2 bp in the length of the spacer between the conserved palindromic sequences. A “rotationally flexible protein homodimers” model has been proposed in order to explain the conformational adaptivity of the 16-3 repressor. In this paper, we report on the isolation of a repressor mutant with altered binding specificity which was used to identify a residue-base pair contact and to monitor the spatial relationship of the recognition helix of C repressor to the contacting major groove of DNA within the two kinds of repressor-operator complexes. Our results indicate spatial differences at the interface which may reflect different docking arrangements in recognition of the structurally different operators by the 16-3 repressor.

A mutational analysis of the ColE1-encoded cell cycle regulator Rcd confirms its role in plasmid stability

Multimers of multicopy plasmids cause instability. They arise by homologous recombination and accumulate by over-replication in a process known as the dimer catastrophe. Dimers are resolved to monomers by site-specific recombination systems such as Xer-cer of plasmid ColE1. In addition, the Rcd checkpoint hypothesis proposes that a short transcript (Rcd) coded within ColE1 cer delays the division of multimer-containing cells. The crucial observation underpinning the checkpoint hypothesis is that when the Rcd promoter (P(cer)) is inactivated by mutation of its invariant T, the plasmid becomes unstable. Recently, we discovered that this mutation also alters a potential Fis binding site in cer. ColE1-like plasmids are less stable in fis mutant hosts and it is conceivable that instability caused by the mutation is due to altered Fis binding, rather than the loss of Rcd expression per se. We have therefore undertaken an independent test of the role of P(cer)-Rcd in multicopy plasmid stability. We have generated a series of loss-of-function mutants of Rcd and detailed analysis of two of these shows that they cause a level of instability indistinguishable from P(cer) inactivation. This result is consistent with the predictions of the checkpoint hypothesis and confirms the role of Rcd in plasmid stability.

Twenty Years of Delila and Molecular Information Theory: The Altenberg-Austin Workshop in Theoretical Biology Biological Information, Beyond Metaphor: Causality, Explanation, and Unification Altenberg, Austria, 11-14 July 2002

A brief personal history is given about how information theory can be applied to binding sites of genetic control molecules on nucleic acids. The primary example used is ribosome binding sites in Escherichia coli. Once the sites are aligned, the information needed to describe the sites can be computed using Claude Shannon's method. This is displayed by a computer graphic called a sequence logo. The logo represents an average binding site, and the mathematics easily allows one to determine the components of this average. That is, given a set of binding sites, the information for individual binding sites can also be computed. One can go further and predict the information of sites that are not in the original data set. Information theory also allows one to model the flexibility of ribosome binding sites, and this led us to a simple model for ribosome translational initiation in which the molecular components fit together only when the ribosome is at a good ribosome binding site. Since information theory is general, the same mathematics applies to human splice junctions, where we can predict the effect of sequence changes that cause human genetic diseases and cancer. The second example given is the Pribnow 'box' which, when viewed by the information theory method, reveals a mechanism for initiation of both transcription and DNA replication. Replication, transcription, splicing, and translation into protein represent the central dogma, so these examples show how molecular information theory is contributing to our knowledge of basic biology.

The seven E. coli ribosomal RNA operon upstream regulatory regions differ in structure and transcription factor binding efficiencies

Ribosomal RNAs in E. coli are transcribed from seven operons, which are highly conserved in their organization and sequence. However, the upstream regulatory DNA regions differ considerably, suggesting differences in regulation. We have therefore analyzed the conformation of all seven DNA elements located upstream of the major E. coli rRNA P1 promoters. As judged by temperature-dependent gel electrophoresis with isolated DNA fragments comprising the individual P1 promoters and the complete upstream regulatory regions, all seven rRNA upstream sequences are intrinsically curved. The degree of intrinsic curvature was highest for the rrnB and rrnD fragments and less pronounced for the rrnA and rrnE operons. Comparison of the experimentally determined differences in curvature with programs for the prediction of DNA conformation revealed a generally high degree of conformity. Moreover, the analysis showed that the center of curvature is located at about the same position in all fragments. The different upstream regions were analyzed for their capacity to bind the transcription factors FIS and H-NS, which are known as antagonists in the regulation of rRNA synthesis. Gel retardation experiments revealed that both proteins interact with the upstream promoter regions of all seven rDNA fragments, with the affinities of the different DNA fragments for FIS and H-NS and the structure of the resulting complexes deviating considerably. FIS binding was non-cooperative, and at comparable protein concentrations the occupancy of the different DNA fragments varied between two and four binding sites. In contrast, H-NS was shown to bind cooperatively and intermediate states of occupancy could not be resolved for each fragment. The different gel electrophoretic mobilities of the individual DNA/protein complexes indicate variable structures and topologies of the upstream activating sequence regulatory complexes. Our results are highly suggestive of differential regulation of the individual rRNA operons.

The role of nucleoid-associated proteins in the organization and compaction of bacterial chromatin

The bacterial chromosomal DNA is folded into a compact structure called nucleoid. The shape and size of this 'body' is determined by a number of factors. Major players are DNA supercoiling, macromolecular crowding and architectural proteins, associated with the nucleoid, which are the topic of this MicroReview. Although many of these proteins were identified more than 25 years ago, the molecular mechanisms involved in the organization and compaction of DNA have only started to become clear in recent years. Many of these new insights can be attributed to the use of recently developed biophysical techniques.

Fis regulates transcriptional induction of RpoS in Salmonella enterica

The sigma factor RpoS is known to regulate at least 60 genes in response to environmental sources of stress or during growth to stationary phase (SP). Accumulation of RpoS relies on integration of multiple genetic controls, including regulation at the levels of transcription, translation, protein stability, and protein activity. Growth to SP in rich medium results in a 30-fold induction of RpoS, although the mechanism of this regulation is not understood. We characterized the activity of promoters serving rpoS in Salmonella enterica serovar Typhimurium and report that regulation of transcription during growth into SP depends on Fis, a DNA-binding protein whose abundance is high during exponential growth and very low in SP. A fis mutant of S. enterica serovar Typhimurium showed a ninefold increase in expression from the major rpoS promoter (PrpoS) during exponential growth, whereas expression during SP was unaffected. Increased transcription from PrpoS in the absence of Fis eliminated the transcriptional induction as cells enter SP. The mutant phenotype can be complemented by wild-type fis carried on a single-copy plasmid. Fis regulation of rpoS requires the presence of a Fis site positioned at -50 with respect to PrpoS, and this site is bound by Fis in vitro. A model is presented in which Fis binding to this site allows repression of rpoS specifically during exponential growth, thus mediating transcriptional regulation of rpoS.

Hierarchical gene regulators adapt to its host milieus

The facultative intracellular pathogen Salmonella enterica serovar Typhimurium possesses an elaborate set of virulence genes that enables the bacterium successfully to move between and adapt to the environment, different host organisms and various micro-niches within a given host. Expression of virulence attributes is by no means constitutive. Rather, the regulation of virulence determinants is highly coordinated and integrated into normal bacterial physiological responses. By integrating discriminating virulence gene regulators with conserved housekeeping regulatory processes, the bacteria can sense alterations in the repertoire of environmental cues, and translate the sensing events into a pragmatic and coordinated expression of virulence genes. While the description of transmissible genetic elements that import global gene regulatory factors into a cell brings conceptual problems into the established regulatory network, the existence of mobile gene regulators may actually enable the bacteria to further modulate virulence expression.

A global role for Fis in the transcriptional control of metabolism and type III secretion in Salmonella enterica serovar Typhimurium

Fis is a key DNA-binding protein involved in nucleoid organization and modulation of many DNA transactions, including transcription in enteric bacteria. The regulon of genes whose expression is influenced by Fis inSalmonella entericaserovar Typhimurium (S. typhimurium) has been defined by DNA microarray analysis. These data suggest that Fis plays a central role in coordinating the expression of both metabolic and type III secretion factors. The genes that were most strongly up-regulated by Fis were those involved in virulence and located in the pathogenicity islands SPI-1, SPI-2, SPI-3 and SPI-5. Similarly, motility and flagellar genes required Fis for full expression. This was shown to be a direct effect as purified Fis protein bound to the promoter regions of representative flagella and SPI-2 genes. Genes contributing to aspects of metabolism known to assist the bacterium during survival in the mammalian gut were also Fis-regulated, usually negatively. This category included components of metabolic pathways for propanediol utilization, biotin synthesis, vitamin B12transport, fatty acids and acetate metabolism, as well as genes for the glyoxylate bypass of the tricarboxylic acid cycle. Genes found to be positively regulated by Fis included those for ethanolamine utilization. The data reported reveal the central role played by Fis in coordinating the expression of both housekeeping and virulence factors required byS. typhimuriumduring life in the gut lumen or during systemic infection of host cells.

Molecular flip-flops formed by overlapping Fis sites

The DNA-binding protein Fis frequently uses pairs of sites 7 or 11 base pairs (bp) apart. Two overlapping Fis sites separated by 11 bp are found in the Escherichia coli origin of chromosomal replication. Only one of these sites is bound by Fis at a time, so the structure is a molecular flip-flop that could direct alternative firing of replication complexes in opposite directions. Alternatively, the flip-flop could represent part of an on-off switch for replication. Because they can be used to create precise switched states, molecular flip-flops could be used as the basis of a novel molecular computer.

The genome and proteome of coliphage T1

The genome of enterobacterial phage T1 has been sequenced, revealing that its 50.7-kb terminally redundant, circularly permuted sequence contains 48,836 bp of nonredundant nucleotides. Seventy-seven open reading frames (ORFs) were identified, with a high percentage of small genes located at the termini of the genomes displaying no homology to existing phage or prophage proteins. Of the genes showing homologs (47%), we identified those involved in host DNA degradation (three endonucleases) and T1 replication (DNA helicase, primase, and single-stranded DNA-binding proteins) and recombination (RecE and Erf homologs). While the tail genes showed homology to those from temperate coliphage N15, the capsid biosynthetic genes were unique. Phage proteins were resolved by 2D gel electrophoresis, and mass spectrometry was used to identify several of the spots including the major head, portal, and tail proteins, thus verifying the annotation.

Size and sequence heterogeneity in the ITS1 of Xylaria hypoxylon isolates

During a survey of 375 strains of the Xylariales, one isolate (F127076) was observed to have an ITS1 size of 833 bp. This size exceeds the average ITS1 size in the Xylariales (mean = 209 +/- 57 bp). Comparison of the DNA sequence with GenBank and with a proprietary DNA database revealed low homology with Xylaria hypoxylon ATCC 42768, and with one undescribed Xylaria species. When the ITS2 sequence was compared, these isolates were 96-98% homologous. Sequences of other variable genes confirmed the relatedness among these strains. A closer observation of the ITS sequence of this isolate revealed the presence of three repeated domains of 250 bp plus one truncated domain, showing 52-75% homology. Sequence similarity suggests that the repeated domain is derived from the fusion of the ITS1 with a DNA fragment derived from the gamma-glutamyltranspeptidase gene. This example suggests that the rate of evolution of ITS1 can be independent of the rate of evolution of other genes, even when this variability is not a result of slipped strand misspairing events like in other Xylariales. This observation also indicates that recombination with other nuclear genes could participate in the evolution of the internal transcribed spacer.

Modulation of CRP-dependent transcription at the Escherichia coli acsP2 promoter by nucleoprotein complexes: anti-activation by the nucleoid proteins FIS and IHF

acs encodes acetyl-coenzyme A synthetase, a high-affinity enzyme that allows cells to scavenge for acetate during carbon starvation. CRP activates acs transcription by binding tandem DNA sites located upstream of the major promoter, acsP2. Here, we used electrophoretic mobility shift assays and DNase I footprint analyses to demonstrate that the nucleoid proteins FIS and IHF each bind multiple sites within the acs regulatory region, that FIS competes successfully with CRP for binding to their overlapping and neighbouring sites and that IHF binds independently of either FIS or CRP. Using in vitro transcription assays, we demonstrated that FIS and IHF independently reduce CRP-dependent acs transcription. Using in vivo reporter assays, we showed that disruption of DNA sites for FIS or deletion of DNA sites for IHF increases acs transcription. We propose that FIS and IHF each function directly as anti-activators of CRP, each working independently at different times during growth to set the levels of CRP-dependent acs transcription.

Genome-wide prediction, display and refinement of binding sites with information theory-based models

BACKGROUND

We present Delila-genome, a software system for identification, visualization and analysis of protein binding sites in complete genome sequences. Binding sites are predicted by scanning genomic sequences with information theory-based (or user-defined) weight matrices. Matrices are refined by adding experimentally-defined binding sites to published binding sites. Delila-Genome was used to examine the accuracy of individual information contents of binding sites detected with refined matrices as a measure of the strengths of the corresponding protein-nucleic acid interactions. The software can then be used to predict novel sites by rescanning the genome with the refined matrices.

RESULTS

Parameters for genome scans are entered using a Java-based GUI interface and backend scripts in Perl. Multi-processor CPU load-sharing minimized the average response time for scans of different chromosomes. Scans of human genome assemblies required 4-6 hours for transcription factor binding sites and 10-19 hours for splice sites, respectively, on 24- and 3-node Mosix and Beowulf clusters. Individual binding sites are displayed either as high-resolution sequence walkers or in low-resolution custom tracks in the UCSC genome browser. For large datasets, we applied a data reduction strategy that limited displays of binding sites exceeding a threshold information content to specific chromosomal regions within or adjacent to genes. An HTML document is produced listing binding sites ranked by binding site strength or chromosomal location hyperlinked to the UCSC custom track, other annotation databases and binding site sequences. Post-genome scan tools parse binding site annotations of selected chromosome intervals and compare the results of genome scans using different weight matrices. Comparisons of multiple genome scans can display binding sites that are unique to each scan and identify sites with significantly altered binding strengths.

CONCLUSIONS

Delila-Genome was used to scan the human genome sequence with information weight matrices of transcription factor binding sites, including PXR/RXRalpha, AHR and NF-kappaB p50/p65, and matrices for RNA binding sites including splice donor, acceptor, and SC35 recognition sites. Comparisons of genome scans with the original and refined PXR/RXRalpha information weight matrices indicate that the refined model more accurately predicts the strengths of known binding sites and is more sensitive for detection of novel binding sites.

Selective regulation of ptsG expression by Fis. Formation of either activating or repressing nucleoprotein complex in response to glucose

Transcription of ptsG encoding glucose-specific permease, enzyme IICB(Glc), in Escherichia coli is initiated from two promoters, P1 and P2. ptsG transcription is repressed by Mlc, a glucose-inducible regulator of carbohydrate metabolism. The regulation of ptsG P1 transcription is also under positive control by cyclic AMP receptor protein and cyclic AMP complex (CRP.cAMP) as observed in other Mlc regulon. We report here that Fis, one of the nucleoid-associated proteins, plays a key role in glucose induction of Mlc regulon. ptsG transcription was induced when wild-type cells were grown in the presence of glucose. However, in a fis mutant, the basal level of ptsG transcription was higher but decreased when cells were grown in the presence of glucose, which implies the possibility of regulatory interactions among Fis, Mlc, and CRP.cAMP. Footprinting experiments with various probes and transcription assays revealed that Fis assists both Mlc repression and CRP.cAMP activation of ptsG P1 through the formation of Fis.CRP.Mlc or Fis.CRP nucleoprotein complexes at ptsG P1 promoter depending on the availability of glucose in the growth medium. ptsG P2 transcription was inhibited by Fis and Mlc. Tighter Mlc repression and enhanced CRP.cAMP activation of ptsG P1 by Fis enable cells to regulate Mlc regulon efficiently by selectively controlling the concentration of enzyme IICB(Glc) that modulates Mlc activity.

The phi29 transcriptional regulator contacts the nucleoid protein p6 to organize a repression complex

The nucleoid protein p6 of Bacillus subtilis phage phi29 binds to DNA, recognizing a structural feature rather than a specific sequence. Upon binding to the viral DNA ends, p6 generates an extended nucleoprotein complex that activates the initiation of phi29 DNA replication. Protein p6 also participates in transcription regulation, repressing the early C2 promoter and assisting the viral regulatory protein p4 in controlling the switch from early to late transcription. Proteins p6 and p4 bind cooperatively to an approximately 200 bp DNA region located between the late A3 and the early A2c promoters, generating an extended nucleoprotein complex that helps to repress the early A2c promoter and to activate the late A3 promoter. We show that stable assembly of this complex requires interaction between protein p6 and the C-terminus of protein p4. Therefore, at this DNA region, stable polymerization of protein p6 relies on p4-specified signals in addition to the structural features of the DNA. Protein p4 would define the phase and boundaries of the p6-DNA complex.

Architecture of Fis-activated transcription complexes at the Escherichia coli rrnB P1 and rrnE P1 promoters

The transcription factor Fis activates the Escherichia coli rRNA promoters rrnB P1 and rrnE P1 by binding to sites centered at -71 and -72, respectively, and interacting with the C-terminal domain of the alpha subunit of RNA polymerase (RNAP alphaCTD). To understand the mechanism of activation by Fis at these promoters, we used oriented alpha-heterodimeric RNAPs and heterodimers of Fis to determine whether one or both subunits of alpha and Fis participate in the alphaCTD-Fis interaction. Our results imply that only one alphaCTD in the alpha dimer and only one activation-proficient subunit in the Fis dimer are required for activation by Fis. A library of alanine substitutions in alpha was used to identify the alphaCTD determinants required for Fis-dependent transcription at rrnB P1 and rrnE P1. We propose that the transcriptional activation region of the promoter-proximal subunit of the Fis dimer interacts with a determinant that includes E273 of one alphaCTD to activate transcription. We further suggest that the Fis contact to alphaCTD results in alphaCTD interactions with DNA that differ somewhat from those that occur at UP elements in the absence of Fis. The accompanying paper shows that the 273 determinant on alphaCTD is also targeted by Fis at the proP P2 promoter where the activator binds overlapping the -35 hexamer. Thus, similar Fis-alphaCTD interactions are used for activation of transcription when the activator is bound at very different positions on the DNA.

Consensus sequence Zen

Consensus sequences are widely used in molecular biology but they have many flaws. As a result, binding sites of proteins and other molecules are missed during studies of genetic sequences and important biological effects cannot be seen. Information theory provides a mathematically robust way to avoid consensus sequences. Instead of using consensus sequences, sequence conservation can be quantitatively presented in bits of information by using sequence logo graphics to represent the average of a set of sites, and sequence walker graphics to represent individual sites.

Strong minor groove base conservation in sequence logos implies DNA distortion or base flipping during replication and transcription initiation

The sequence logo for DNA binding sites of the bacteriophage P1 replication protein RepA shows unusually high sequence conservation ( approximately 2 bits) at a minor groove that faces RepA. However, B-form DNA can support only 1 bit of sequence conservation via contacts into the minor groove. The high conservation in RepA sites therefore implies a distorted DNA helix with direct or indirect contacts to the protein. Here I show that a high minor groove conservation signature also appears in sequence logos of sites for other replication origin binding proteins (Rts1, DnaA, P4 alpha, EBNA1, ORC) and promoter binding proteins (sigma(70), sigma(D) factors). This finding implies that DNA binding proteins generally use non-B-form DNA distortion such as base flipping to initiate replication and transcription.

The P1 phage replication protein RepA contacts an otherwise inaccessible thymine N3 proton by DNA distortion or base flipping

The RepA protein from bacteriophage P1 binds DNA to initiate replication. RepA covers one face of the DNA and the binding site has a completely conserved T that directly faces RepA from the minor groove at position +7. Although all four bases can be distinguished through contacts in the major groove of B-form DNA, contacts in the minor groove cannot easily distinguish between A and T bases. Therefore the 100% conservation at this position cannot be accounted for by direct contacts approaching into the minor groove of B-form DNA. RepA binding sites with modified base pairs at position +7 were used to investigate contacts with RepA. The data show that RepA contacts the N3 proton of T at position +7 and that the T=A hydrogen bonds are already broken in the DNA before RepA binds. To accommodate the N3 proton contact the T(+7 )/A(+7)((')) base pair must be distorted. One possibility is that T(+7) is flipped out of the helix. The energetics of the contact allows RepA to distinguish between all four bases, accounting for the observed high sequence conservation. After protein binding, base pair distortion or base flipping could initiate DNA melting as the second step in DNA replication.

Contributions of UP elements and the transcription factor FIS to expression from the seven rrn P1 promoters in Escherichia coli

The high activity of the rrnB P1 promoter in Escherichia coli results from a cis-acting DNA sequence, the UP element, and a trans-acting transcription factor, FIS. In this study, we examine the effects of FIS and the UP element at the other six rrn P1 promoters. We find that UP elements are present at all of the rrn P1 promoters, but they make different relative contributions to promoter activity. Similarly, FIS binds upstream of, and activates, all seven rrn P1 promoters but to different extents. The total number of FIS binding sites, as well as their positions relative to the transcription start site, differ at each rrn P1 promoter. Surprisingly, the FIS sites upstream of site I play a much larger role in transcription from most rrn P1 promoters compared to rrnB P1. Our studies indicate that the overall activities of the seven rrn P1 promoters are similar, and the same contributors are responsible for these high activities, but these inputs make different relative contributions and may act through slightly different mechanisms at each promoter. These studies have implications for the control of gene expression of unlinked multigene families.

Anatomy of Escherichia coli ribosome binding sites

During translational initiation in prokaryotes, the 3' end of the 16S rRNA binds to a region just upstream of the initiation codon. The relationship between this Shine-Dalgarno (SD) region and the binding of ribosomes to translation start-points has been well studied, but a unified mathematical connection between the SD, the initiation codon and the spacing between them has been lacking. Using information theory, we constructed a model that treats these three components uniformly by assigning to the SD and the initiation region (IR) conservations in bits of information, and by assigning to the spacing an uncertainty, also in bits. To build the model, we first aligned the SD region by maximizing the information content there. The ease of this process confirmed the existence of the SD pattern within a set of 4122 reviewed and revised Escherichia coli gene starts. This large data set allowed us to show graphically, by sequence logos, that the spacing between the SD and the initiation region affects both the SD site conservation and its pattern. We used the aligned SD, the spacing, and the initiation region to model ribosome binding and to identify gene starts that do not conform to the ribosome binding site model. A total of 569 experimentally proven starts are more conserved (have higher information content) than the full set of revised starts, which probably reflects an experimental bias against the detection of gene products that have inefficient ribosome binding sites. Models were refined cyclically by removing non-conforming weak sites. After this procedure, models derived from either the original or the revised gene start annotation were similar. Therefore, this information theory-based technique provides a method for easily constructing biologically sensible ribosome binding site models. Such models should be useful for refining gene-start predictions of any sequenced bacterial genome.

Growth phase and growth rate regulation of the rapA gene, encoding the RNA polymerase-associated protein RapA in Escherichia coli

The Escherichia coli rapA gene encodes the RNA polymerase (RNAP)-associated protein RapA, which is a bacterial member of the SWI/SNF helicase-like protein family. We have studied the rapA promoter and its regulation in vivo and determined the interaction between RNAP and the promoter in vitro. We have found that the expression of rapA is growth phase dependent, peaking at the early log phase. The growth phase control of rapA is determined at least by one particular feature of the promoter: it uses CTP as the transcription-initiating nucleotide instead of a purine, which is used for most E. coli promoters. We also found that the rapA promoter is subject to growth rate regulation in vivo and that it forms intrinsic unstable initiation complexes with RNAP in vitro. Furthermore, we have shown that a GC-rich or discriminator sequence between the -10 and +1 positions of the rapA promoter is responsible for its growth rate control and the instability of its initiation complexes with RNAP.

Involvement of FIS in the H-NS-mediated regulation of virF gene of Shigella and enteroinvasive Escherichia coli

The mechanism of pathogenicity in Shigella and enteroinvasive Escherichia coli (EIEC) requires the co-ordinated expression of several genes located on both the virulence plasmid and the chromosome. We found that cells lacking a functional FIS protein (factor for inversion stimulation) are partially impaired in expressing the virulence genes and that full expression is totally restored when Shigella wild-type fis gene is offered in trans. We also identified virF, among the virulence genes, as a target of FIS-mediated activation and showed that FIS binds to four specific sites in the promoter region of virF. Previous studies have demonstrated that the expression of VirF, the first positive activator of a multistep regulatory cascade, is subject to temperature-dependent regulation by H-NS, one of the main nucleoid-associated proteins. We now demonstrate that two of the four FIS sites overlap one of the two H-NS sites responsible for thermoregulation (H-NS site I). FIS was found to exercise a direct positive transcriptional control at permissive temperature (37 degrees C), when H-NS fails to repress virF, as well as an indirect effect by partially counteracting H-NS inhibition at the transition temperature (32 degrees C). Our data indicate that FIS may be relevant for the rapid increase in virF expression after penetration of bacteria into the host.

Control of transcription by nucleoid proteins

Nucleoid proteins are a group of abundant DNA binding proteins that modulate the structure of the bacterial chromosome. They have been recruited as specific negative and positive regulators of gene transcription and their fluctuating patterns of expression are often exploited to impart an additional level of control with respect to environmental conditions.

Genome organisation and chromatin structure in Escherichia coli

We have analysed the complete sequence of the Escherichia coli K12 isolate MG1655 genome for chromatin-associated protein binding sites, and compared the predicted location of predicted sites with experimental expression data from 'DNA chip' experiments. Of the dozen proteins associated with chromatin in E. coli, only three have been shown to have significant binding preferences: integration host factor (IHF) has the strongest binding site preference, and FIS sites show a weak consensus, and there is no clear consensus site for binding of the H-NS protein. Using hidden Markov models (HMMs), we predict the location of 608 IHF sites, scattered throughout the genome. A subset of the IHF sites associated with repeats tends to be clustered around the origin of replication. We estimate there could be roughly 6000 FIS sites in E. coli, and the sites tend to be localised in two regions flanking the replication termini. We also show that the regions upstream of genes regulated by H-NS are more curved and have a higher AT content than regions upstream of other genes. These regions in general would also be localised near the replication terminus.

Transposon Tn7 gene insertion into an evolutionarily conserved human homolog of Escherichia coli attTn7

Escherichia coli transposon Tn7 can integrate into its target DNA sequence, attTn7 at the 3' end of glmS, with high specificity and efficiency. Remarkably, the insertional recognition sequence in the E. coli genome displays a high degree of identity with the corresponding region at the 3' end of the corresponding human gene for glutamine-fructose-6-phosphate transaminase (GFPT), located at 2p13. It was therefore of interest to determine whether Tn7 could recognize the corresponding human sequence, and transpose at that site. Strains of E. coli DH5alpha were prepared carrying the tnsA-E genes on one plasmid, and attTn7 or the human equivalent on a second recipient plasmid within the alpha-complementation fragment of the lacZ gene. Each strain was transformed with a donor plasmid carrying a gentamycin resistance gene within the Tn7L and Tn7R cassettes. Restriction mapping and sequence analysis of recipient plasmids isolated from white colonies demonstrated that Tn7 inserted the gentamycin resistance gene both into the E. coli attTn7 sequence, and into its human counterpart. No nonspecific insertion was observed in a control plasmid containing only the lacZ fragment. These results provide a basis to investigate whether TnsA-D proteins can mediate gene insertion into comparably conserved sites in eukaryotic chromosomes.

Information analysis of human splice site mutations

Splice site nucleotide substitutions can be analyzed by comparing the individual information contents (Ri, bits) of the normal and variant splice junction sequences [Rogan and Schneider, 1995]. In the present study, we related splicing abnormalities to changes in Ri values of 111 previously reported splice site substitutions in 41 different genes. Mutant donor and acceptor sites have significantly less information than their normal counterparts. With one possible exception, primary mutant sites with <2.4 bits were not spliced. Sites with Ri values > or = 2.4 bits but less than the corresponding natural site usually decreased, but did not abolish splicing. Substitutions that produced small changes in Ri probably do not impair splicing and are often polymorphisms. The Ri values of activated cryptic sites were generally comparable to or greater than those of the corresponding natural splice sites. Information analysis revealed preexisting cryptic splice junctions that are used instead of the mutated natural site. Other cryptic sites were created or strengthened by sequence changes that simultaneously altered the natural site. Comparison between normal and mutant splice site Ri values distinguishes substitutions that impair splicing from those which do not, distinguishes null alleles from those that are partially functional, and detects activated cryptic splice sites.

Antagonistic control of the Escherichia coli bgl promoter by FIS and CAP in vitro

The wild-type Escherichia coli bgl promoter is silent in vivo but active in vitro. Silencing in vivo is directed by silencer sequences that flank the promoter, and requires nucleoid-associated protein H-NS and other unidentified cellular factors. Here we show that the DNA bending protein FIS is a repressor of the bgl promoter. Two FIS binding sites, centred at positions -52 and -27, overlap the CAP binding site and the -35 box respectively. FIS efficiently competes with CAP for binding to the wild-type promoter. However, FIS does not prevent binding of RNA polymerase. It interferes with the formation of a heparin-resistant complex and represses transcription initiation up to 40-fold. The presence of CAP has very little effect on the FIS-mediated repression of the wild-type bgl promoter in vitro. However, when a bgl promoter allele was tested that carries an improved CAP binding site (which leads to activation in vivo) CAP effectively counteracted repression by FIS in vitro. These results suggest that FIS contributes to silencing of the wild-type bgl promoter in vivo, presumably in the early exponential phase when FIS is predominantly expressed.

Activation of Escherichia coli leuV transcription by FIS

The transcription factor FIS has been implicated in the regulation of several stable RNA promoters, including that for the major tRNALeu species in Escherichia coli, tRNA1Leu. However, no evidence for direct involvement of FIS in tRNA1Leu expression has been reported. We show here that FIS binds to a site upstream of the leuV promoter (centered at -71) and that it directly stimulates leuV transcription in vitro. A mutation in the FIS binding site reduces transcription from a leuV promoter in strains containing FIS but has no effect on transcription in strains lacking FIS, indicating that FIS contributes to leuV expression in vivo. We also find that RNA polymerase forms an unusual heparin-sensitive complex with the leuV promoter, having a downstream protection boundary of approximately -7, and that the first two nucleotides of the transcript, GTP and UTP, are required for formation of a heparin-stable complex that extends downstream of the transcription start site. These studies have implications for the regulation of leuV transcription.

In vitro selection of integration host factor binding sites

Integration host factor (IHF) is a bacterial protein that binds and severely bends a specific DNA target. IHF binding sites are approximately 30 to 35 bp long and are apparently divided into two domains. While the 3' domain is conserved, the 5' domain is degenerate but is typically AT rich. As a result of physical constraints that IHF must impose on DNA in order to bind, it is believed that this 5' domain must possess structural characteristics conducive for both binding and bending with little regard for specific contacts between the protein and the DNA. We have examined the sequence requirements of the 5' binding domain of the IHF binding target. Using a SELEX procedure, we randomized and selected variants of a natural IHF site. We then analyzed these variants to determine how the 5' binding domain affects the structure, affinity, and function of an IHF-DNA complex in a native system. Despite finding individual sequences that varied over 100-fold in affinity for IHF, we found no apparent correlation between affinity and function.

Using sequence logos and information analysis of Lrp DNA binding sites to investigate discrepancies between natural selection and SELEX

In vitro experiments that characterize DNA-protein interactions by artificial selection, such as SELEX,are often performed with the assumption that the experimental conditions are equivalent to natural ones. To test whether SELEX gives natural results, we compared sequence logos composed from naturally occurring leucine-responsive regulatory protein (Lrp) binding sites with those composed from SELEX-generated binding sites. The sequence logos were significantly different, indicating that the binding conditions are disparate. A likely explanation is that the SELEX experiment selected for a dimeric or trimeric Lrp complex bound to DNA. In contrast, natural sites appear to be bound by a monomer. This discrepancy suggests that in vitro selections do not necessarily give binding site sets comparable with the natural binding sites.

Evolution of new protein function: recombinational enhancer Fis originated by horizontal gene transfer from the transcriptional regulator NtrC

New protein function is thought to evolve mostly by gene duplication and divergence. Here we present phylogenetic evidence that the multifunctional protein Fis of the gamma proteobacterial species derived from the COOH-terminal domain of an ancestral alpha proteobacterial NtrC transcriptional regulatory protein. All of the known enterobacterial fis genes are preceded by an open reading frame, named yhdG, that is highly similar to nifR3, a gene that forms an operon with ntrC in several alpha proteobacterial species. Thus, we propose that yhdG and fis were acquired by a lineage ancestral to the gamma proteobacteria in a single horizontal gene transfer event, and later diverged to their present functions.