Information analysis of Fis binding sites

Spatio-temporal organization of the E. coli chromosome from base to cellular length scales

Escherichia coli has been a vital model organism for studying chromosomal structure, thanks, in part, to its small and circular genome (4.6 million base pairs) and well-characterized biochemical pathways. Over the last several decades, we have made considerable progress in understanding the intricacies of the structure and subsequent function of the E. coli nucleoid. At the smallest scale, DNA, with no physical constraints, takes on a shape reminiscent of a randomly twisted cable, forming mostly random coils but partly affected by its stiffness. This ball-of-spaghetti-like shape forms a structure several times too large to fit into the cell. Once the physiological constraints of the cell are added, the DNA takes on overtwisted (negatively supercoiled) structures, which are shaped by an intricate interplay of many proteins carrying out essential biological processes. At shorter length scales (up to about 1 kb), nucleoid-associated proteins organize and condense the chromosome by inducing loops, bends, and forming bridges. Zooming out further and including cellular processes, topological domains are formed, which are flanked by supercoiling barriers. At the megabase-scale both large, highly self-interacting regions (macrodomains) and strong contacts between distant but co-regulated genes have been observed. At the largest scale, the nucleoid forms a helical ellipsoid. In this review, we will explore the history and recent advances that pave the way for a better understanding of E. coli chromosome organization and structure, discussing the cellular processes that drive changes in DNA shape, and what contributes to compaction and formation of dynamic structures, and in turn how bacterial chromatin affects key processes such as transcription and replication.

Bacterial chromatin proteins, transcription, and DNA topology: Inseparable partners in the control of gene expression

DNA in bacterial chromosomes is organized into higher-order structures by DNA-binding proteins called nucleoid-associated proteins (NAPs) or bacterial chromatin proteins (BCPs). BCPs often bind to or near DNA loci transcribed by RNA polymerase (RNAP) and can either increase or decrease gene expression. To understand the mechanisms by which BCPs alter transcription, one must consider both steric effects and the topological forces that arise when DNA deviates from its fully relaxed double-helical structure. Transcribing RNAP creates DNA negative (-) supercoils upstream and positive (+) supercoils downstream whenever RNAP and DNA are unable to rotate freely. This (-) and (+) supercoiling generates topological forces that resist forward translocation of DNA through RNAP unless the supercoiling is constrained by BCPs or relieved by topoisomerases. BCPs also may enhance topological stress and overall can either inhibit or aid transcription. Here, we review current understanding of how RNAP, BCPs, and DNA topology interplay to control gene expression.

Molecular Information Theory Meets Protein Folding

We propose an application of molecular information theory to analyze the folding of single domain proteins. We analyze results from various areas of protein science, such as sequence-based potentials, reduced amino acid alphabets, backbone configurational entropy, secondary structure content, residue burial layers, and mutational studies of protein stability changes. We found that the average information contained in the sequences of evolved proteins is very close to the average information needed to specify a fold ∼2.2 ± 0.3 bits/(site·operation). The effective alphabet size in evolved proteins equals the effective number of conformations of a residue in the compact unfolded state at around 5. We calculated an energy-to-information conversion efficiency upon folding of around 50%, lower than the theoretical limit of 70%, but much higher than human-built macroscopic machines. We propose a simple mapping between molecular information theory and energy landscape theory and explore the connections between sequence evolution, configurational entropy, and the energetics of protein folding.

The Fis Nucleoid Protein Negatively Regulates the Phase Variation fimS Switch of the Type 1 Pilus Operon in Enteropathogenic Escherichia coli

In Escherichia coli the expression of type 1 pili (T1P) is determined by the site-specific inversion of the fimS ON–OFF switch located immediately upstream of major fimbrial subunit gene fimA. Here we investigated the role of virulence (Ler, GrlR, and GrlA) and global regulators (H-NS, IHF, and Fis) in the regulation of the fimS switch in the human enteropathogenic E. coli (EPEC) O127:H6 strain E2348/69. This strain does not produce detectable T1P and PCR analysis of the fimS switch confirmed that it is locked in the OFF orientation. Among the regulator mutants analyzed, only the ∆fis mutant produced significantly high levels of T1P on its surface and yielded high titers of agglutination of guinea pig erythrocytes. Expression analysis of the fimA, fimB, and fimE promoters using lacZ transcriptional fusions indicated that only PfimA activity is enhanced in the absence of Fis. Collectively, these data demonstrate that Fis is a negative regulator of T1P expression in EPEC and suggest that it is required for the FimE-dependent inversion of the fimS switch from the ON-to-OFF direction. It is possible that a similar mechanism of T1P regulation exists in other intestinal and extra-intestinal pathogenic classes of E. coli.

Blocking, Bending, and Binding: Regulation of Initiation of Chromosome Replication During the Escherichia coli Cell Cycle by Transcriptional Modulators That Interact With Origin DNA

Genome duplication is a critical event in the reproduction cycle of every cell. Because all daughter cells must inherit a complete genome, chromosome replication is tightly regulated, with multiple mechanisms focused on controlling when chromosome replication begins during the cell cycle. In bacteria, chromosome duplication starts when nucleoprotein complexes, termed orisomes, unwind replication origin (oriC) DNA and recruit proteins needed to build new replication forks. Functional orisomes comprise the conserved initiator protein, DnaA, bound to a set of high and low affinity recognition sites in oriC. Orisomes must be assembled each cell cycle. In Escherichia coli, the organism in which orisome assembly has been most thoroughly examined, the process starts with DnaA binding to high affinity sites after chromosome duplication is initiated, and orisome assembly is completed immediately before the next initiation event, when DnaA interacts with oriC’s lower affinity sites, coincident with origin unwinding. A host of regulators, including several transcriptional modulators, targets low affinity DnaA-oriC interactions, exerting their effects by DNA bending, blocking access to recognition sites, and/or facilitating binding of DnaA to both DNA and itself. In this review, we focus on orisome assembly in E. coli. We identify three known transcriptional modulators, SeqA, Fis (factor for inversion stimulation), and IHF (integration host factor), that are not essential for initiation, but which interact directly with E. coli oriC to regulate orisome assembly and replication initiation timing. These regulators function by blocking sites (SeqA) and bending oriC DNA (Fis and IHF) to inhibit or facilitate cooperative low affinity DnaA binding. We also examine how the growth rate regulation of Fis levels might modulate IHF and DnaA binding to oriC under a variety of nutritional conditions. Combined, the regulatory mechanisms mediated by transcriptional modulators help ensure that at all growth rates, bacterial chromosome replication begins once, and only once, per cell cycle.

A horizontally acquired Legionella genomic island encoding a LuxR type regulator and effector proteins displays variation in gene content and regulation

The intracellular pathogen Legionella pneumophila translocates >300 effector proteins into host cells, many of which are regulated at the transcriptional level. Here, we describe a novel L. pneumophila genomic island, which undergoes horizontal gene transfer within the Legionella genus. This island encodes two Icm/Dot effectors: LegK3 and a previously uncharacterized effector which we named CegK3, as well as a LuxR type regulator, which we named RegK3. Analysis of this island in different Legionella species revealed a conserved regulatory element located upstream to the effector-encoding genes in the island. Further analyses, including gene expression analysis, mutagenesis of the RegK3 regulatory element, controlled expression studies, and gel-mobility shift assays, all demonstrate that RegK3 directly activates the expression levels of legK3 and cegK3 effector-encoding genes. Additionally, the expression of all the components of the island is silenced by the Fis repressors. Comparison of expression profiles of these three genes among different Legionella species revealed variability in the activation levels mediated by RegK3, which were positively correlated with the Fis-mediated repression. Furthermore, LegK3 and CegK3 effectors moderately inhibit yeast growth, and importantly, they have a strong synergistic inhibitory effect on yeast growth, suggesting these two effectors are not only co-regulated but also might function together.

Fis protein forms DNA topological barriers to confine transcription-coupled DNA supercoiling in Escherichia coli

Previously, we demonstrated that transcription-coupled DNA supercoiling (TCDS) potently activated or inhibited nearby promoters in Escherichia coli even in the presence of all four DNA topoisomerases, suggesting that DNA topoisomerases are not the only factors regulating TCDS. A different mechanism exists to confine this localized DNA supercoiling. Using an in vivo system containing the TCDS-activated leu-500 promoter (Pleu-500 ), we find that the nucleoid-associated Fis protein potently inhibits the TCDS-mediated activation of Pleu-500 . We also find that deletion of the fis gene significantly enhances TCDS-mediated inhibition of transcription of three genes purH, yieP, and yrdA divergently coupled to different rrn operons in the early log phase. These results suggest that Fis protein forms DNA topological barriers upon binding to its recognition sites, blocks TCDS diffusion, and potently inhibits the TCDS-activated Pleu-500 .

Restriction enzymes use a 24 dimensional coding space to recognize 6 base long DNA sequences

Restriction enzymes recognize and bind to specific sequences on invading bacteriophage DNA. Like a key in a lock, these proteins require many contacts to specify the correct DNA sequence. Using information theory we develop an equation that defines the number of independent contacts, which is the dimensionality of the binding. We show that EcoRI, which binds to the sequence GAATTC, functions in 24 dimensions. Information theory represents messages as spheres in high dimensional spaces. Better sphere packing leads to better communications systems. The densest known packing of hyperspheres occurs on the Leech lattice in 24 dimensions. We suggest that the single protein EcoRI molecule employs a Leech lattice in its operation. Optimizing density of sphere packing explains why 6 base restriction enzymes are so common.

Cooperative DNA binding by proteins through DNA shape complementarity

Localized arrays of proteins cooperatively assemble onto chromosomes to control DNA activity in many contexts. Binding cooperativity is often mediated by specific protein-protein interactions, but cooperativity through DNA structure is becoming increasingly recognized as an additional mechanism. During the site-specific DNA recombination reaction that excises phage λ from the chromosome, the bacterial DNA architectural protein Fis recruits multiple λ-encoded Xis proteins to the attR recombination site. Here, we report X-ray crystal structures of DNA complexes containing Fis + Xis, which show little, if any, contacts between the two proteins. Comparisons with structures of DNA complexes containing only Fis or Xis, together with mutant protein and DNA binding studies, support a mechanism for cooperative protein binding solely by DNA allostery. Fis binding both molds the minor groove to potentiate insertion of the Xis β-hairpin wing motif and bends the DNA to facilitate Xis-DNA contacts within the major groove. The Fis-structured minor groove shape that is optimized for Xis binding requires a precisely positioned pyrimidine-purine base-pair step, whose location has been shown to modulate minor groove widths in Fis-bound complexes to different DNA targets.

Endogenous and Foreign Nucleoid-Associated Proteins of Bacteria: Occurrence, Interactions and Effects on Mobile Genetic Elements and Host's Biology

Mobile Genetic Elements (MGEs) are mosaics of functional gene modules of diverse evolutionary origin and are generally divergent from the hosts´ genetic background. Existing biases in base composition and codon usage of these elements` genes impose transcription and translation limitations that may affect the physical and regulatory integration of MGEs in new hosts. Stable appropriation of the foreign DNA depends on a number of host factors among which are the Nucleoid-Associated Proteins (NAPs). These small, basic, highly abundant proteins bind and bend DNA, altering its topology and folding, thereby affecting all known essential DNA metabolism related processes. Both chromosomally- (endogenous) and MGE- (foreign) encoded NAPs have been shown to exist in bacteria. While the role of host-encoded NAPs in xenogeneic silencing of both episomal (plasmids) and integrative MGEs (pathogenicity islands and prophages) is well acknowledged, less is known about the role of MGE-encoded NAPs in the foreign elements biology or their influence on the host's chromosome expression dynamics. Here we review existing literature on the topic, present examples on the positive and negative effects that endogenous and foreign NAPs exert on global transcriptional gene expression, MGE integrative and excisive recombination dynamics, persistence and transfer to suitable hosts and discuss the nature and relevance of synergistic and antagonizing higher order interactions between diverse types of NAPs.

DNA Sequence Determinants Controlling Affinity, Stability and Shape of DNA Complexes Bound by the Nucleoid Protein Fis

The abundant Fis nucleoid protein selectively binds poorly related DNA sequences with high affinities to regulate diverse DNA reactions. Fis binds DNA primarily through DNA backbone contacts and selects target sites by reading conformational properties of DNA sequences, most prominently intrinsic minor groove widths. High-affinity binding requires Fis-stabilized DNA conformational changes that vary depending on DNA sequence. In order to better understand the molecular basis for high affinity site recognition, we analyzed the effects of DNA sequence within and flanking the core Fis binding site on binding affinity and DNA structure. X-ray crystal structures of Fis-DNA complexes containing variable sequences in the noncontacted center of the binding site or variations within the major groove interfaces show that the DNA can adapt to the Fis dimer surface asymmetrically. We show that the presence and position of pyrimidine-purine base steps within the major groove interfaces affect both local DNA bending and minor groove compression to modulate affinities and lifetimes of Fis-DNA complexes. Sequences flanking the core binding site also modulate complex affinities, lifetimes, and the degree of local and global Fis-induced DNA bending. In particular, a G immediately upstream of the 15 bp core sequence inhibits binding and bending, and A-tracts within the flanking base pairs increase both complex lifetimes and global DNA curvatures. Taken together, our observations support a revised DNA motif specifying high-affinity Fis binding and highlight the range of conformations that Fis-bound DNA can adopt. The affinities and DNA conformations of individual Fis-DNA complexes are likely to be tailored to their context-specific biological functions.

The Nucleoid: an Overview

This review provides a brief review of the current understanding of the structure-function relationship of the Escherichia coli nucleoid developed after the overview by Pettijohn focusing on the physical properties of nucleoids. Isolation of nucleoids requires suppression of DNA expansion by various procedures. The ability to control the expansion of nucleoids in vitro has led to purification of nucleoids for chemical and physical analyses and for high-resolution imaging. Isolated E. coli genomes display a number of individually intertwined supercoiled loops emanating from a central core. Metabolic processes of the DNA double helix lead to three types of topological constraints that all cells must resolve to survive: linking number, catenates, and knots. The major species of nucleoid core protein share functional properties with eukaryotic histones forming chromatin; even the structures are different from histones. Eukaryotic histones play dynamic roles in the remodeling of eukaryotic chromatin, thereby controlling the access of RNA polymerase and transcription factors to promoters. The E. coli genome is tightly packed into the nucleoid, but, at each cell division, the genome must be faithfully replicated, divided, and segregated. Nucleoid activities such as transcription, replication, recombination, and repair are all affected by the structural properties and the special conformations of nucleoid. While it is apparent that much has been learned about the nucleoid, it is also evident that the fundamental interactions organizing the structure of DNA in the nucleoid still need to be clearly defined.

Interpretation of mRNA splicing mutations in genetic disease: review of the literature and guidelines for information-theoretical analysis

The interpretation of genomic variants has become one of the paramount challenges in the post-genome sequencing era. In this review we summarize nearly 20 years of research on the applications of information theory (IT) to interpret coding and non-coding mutations that alter mRNA splicing in rare and common diseases. We compile and summarize the spectrum of published variants analyzed by IT, to provide a broad perspective of the distribution of deleterious natural and cryptic splice site variants detected, as well as those affecting splicing regulatory sequences. Results for natural splice site mutations can be interrogated dynamically with Splicing Mutation Calculator, a companion software program that computes changes in information content for any splice site substitution, linked to corresponding publications containing these mutations. The accuracy of IT-based analysis was assessed in the context of experimentally validated mutations. Because splice site information quantifies binding affinity, IT-based analyses can discern the differences between variants that account for the observed reduced (leaky) versus abolished mRNA splicing. We extend this principle by comparing predicted mutations in natural, cryptic, and regulatory splice sites with observed deleterious phenotypic and benign effects. Our analysis of 1727 variants revealed a number of general principles useful for ensuring portability of these analyses and accurate input and interpretation of mutations. We offer guidelines for optimal use of IT software for interpretation of mRNA splicing mutations.

Interpretation of mRNA splicing mutations in genetic disease: review of the literature and guidelines for information-theoretical analysis

The interpretation of genomic variants has become one of the paramount challenges in the post-genome sequencing era. In this review we summarize nearly 20 years of research on the applications of information theory (IT) to interpret coding and non-coding mutations that alter mRNA splicing in rare and common diseases. We compile and summarize the spectrum of published variants analyzed by IT, to provide a broad perspective of the distribution of deleterious natural and cryptic splice site variants detected, as well as those affecting splicing regulatory sequences. Results for natural splice site mutations can be interrogated dynamically with Splicing Mutation Calculator, a companion software program that computes changes in information content for any splice site substitution, linked to corresponding publications containing these mutations. The accuracy of IT-based analysis was assessed in the context of experimentally validated mutations. Because splice site information quantifies binding affinity, IT-based analyses can discern the differences between variants that account for the observed reduced (leaky) versus abolished mRNA splicing. We extend this principle by comparing predicted mutations in natural, cryptic, and regulatory splice sites with observed deleterious phenotypic and benign effects. Our analysis of 1727 variants revealed a number of general principles useful for ensuring portability of these analyses and accurate input and interpretation of mutations. We offer guidelines for optimal use of IT software for interpretation of mRNA splicing mutations.

DNA recognition by a σ(54) transcriptional activator from Aquifex aeolicus

Transcription initiation by bacterial σ(54)-polymerase requires the action of a transcriptional activator protein. Activators bind sequence-specifically upstream of the transcription initiation site via a DNA-binding domain (DBD). The structurally characterized DBDs from activators all belong to the Fis (factor for inversion stimulation) family of helix-turn-helix DNA-binding proteins. We report here structures of the free and DNA-bound forms of the DBD of NtrC4 (4DBD) from Aquifex aeolicus, a member of the NtrC family of σ(54) activators. Two NtrC4-binding sites were identified upstream (-145 and -85bp) from the start of the lpxC gene, which is responsible for the first committed step in lipid A biosynthesis. This is the first experimental evidence for σ(54) regulation in lpxC expression. 4DBD was crystallized both without DNA and in complex with the -145-binding site. The structures, together with biochemical data, indicate that NtrC4 binds to DNA in a manner that is similar to that of its close homolog, Fis. The greater sequence specificity for the binding of 4DBD relative to Fis seems to arise from a larger number of base-specific contacts contributing to affinity than for Fis.

Two Fis regulators directly repress the expression of numerous effector-encoding genes in Legionella pneumophila

Legionella pneumophila is an intracellular human pathogen that utilizes the Icm/Dot type IVB secretion system to translocate a large repertoire of effectors into host cells. For most of these effectors, there is no information regarding their regulation. Therefore, the aim of this study was to examine the involvement of the three L. pneumophila Fis homologs in the regulation of effector-encoding genes. Deletion mutants constructed in the genes encoding the three Fis regulators revealed that Fis1 (lpg0542 gene) and Fis3 (lpg1743) but not Fis2 (lpg1370) are partially required for intracellular growth of L. pneumophila in Acanthamoeba castellanii. To identify pathogenesis-related genes directly regulated by Fis, we established a novel in vivo system which resulted in the discovery of numerous effector-encoding genes directly regulated by Fis. Further examination of these genes revealed that Fis1 and Fis3 repress the level of expression of effector-encoding genes during exponential phase. Three groups of effector-encoding genes were identified: (i) effectors regulated mainly by Fis1, (ii) effectors regulated mainly by Fis3, and (iii) effectors regulated by both Fis1 and Fis3. Examination of the upstream regulatory region of all of these effector-encoding genes revealed multiple putative Fis regulatory elements, and site-directed mutagenesis confirmed that a few of these sites constitute part of a repressor binding element. Furthermore, gel mobility shift assays demonstrated the direct relation between the Fis1 and Fis3 regulators and these regulatory elements. Collectively, our results demonstrate for the first time that two of the three L. pneumophila Fis regulators directly repress the expression of Icm/Dot effector-encoding genes.

Structural change of DNA induced by nucleoid proteins: growth phase-specific Fis and stationary phase-specific Dps

The effects of nucleoid proteins Fis and Dps of Escherichia coli on the higher order structure of a giant DNA were studied, in which Fis and Dps are known to be expressed mainly in the exponential growth phase and stationary phase, respectively. Fis causes loose shrinking of the higher order structure of a genome-sized DNA, T4 DNA (166 kbp), in a cooperative manner, that is, the DNA conformational transition proceeds through the appearance of a bimodal size distribution or the coexistence of elongated coil and shrunken globular states. The effective volume of the loosely shrunken state induced by Fis is 30-60 times larger than that of the compact state induced by spermidine, suggesting that cellular enzymes can access for DNA with the shrunken state but cannot for the compact state. Interestingly, Dps tends to inhibit the Fis-induced shrinkage of DNA, but promotes DNA compaction in the presence of spermidine. These characteristic effects of nucleotide proteins on a giant DNA are discussed by adopting a simple theoretical model with a mean-field approximation.

Analysis of DevR regulated genes in Mycobacterium tuberculosis

The DevRS two component system of Mycobacterium tuberculosis is responsible for its dormancy in host and becomes operative under hypoxic condition. It is experimentally known that phosphorylated DevR controls the expression of several downstream genes in a complex manner. In the present work we propose a theoretical model to show role of binding sites in DevR mediated gene expression. Individual and collective role of binding sites in regulating DevR mediated gene expression has been shown via modeling. Objective of the present work is twofold. First, to describe qualitatively the temporal dynamics of wild type genes and their known mutants. Based on these results we propose that DevR controlled gene expression follows a specific pattern which is efficient in describing other DevR mediated gene expression. Second, to analyze behavior of the system from information theoretical point of view. Using the tools of information theory we have calculated molecular efficiency of the system and have shown that it is close to the maximum limit of isothermal efficiency.

A fast weak motif-finding algorithm based on community detection in graphs

BACKGROUND

Identification of transcription factor binding sites (also called 'motif discovery') in DNA sequences is a basic step in understanding genetic regulation. Although many successful programs have been developed, the problem is far from being solved on account of diversity in gene expression/regulation and the low specificity of binding sites. State-of-the-art algorithms have their own constraints (e.g., high time or space complexity for finding long motifs, low precision in identification of weak motifs, or the OOPS constraint: one occurrence of the motif instance per sequence) which limit their scope of application.

RESULTS

In this paper, we present a novel and fast algorithm we call TFBSGroup. It is based on community detection from a graph and is used to discover long and weak (l,d) motifs under the ZOMOPS constraint (zero, one or multiple occurrence(s) of the motif instance(s) per sequence), where l is the length of a motif and d is the maximum number of mutations between a motif instance and the motif itself. Firstly, TFBSGroup transforms the (l, d) motif search in sequences to focus on the discovery of dense subgraphs within a graph. It identifies these subgraphs using a fast community detection method for obtaining coarse-grained candidate motifs. Next, it greedily refines these candidate motifs towards the true motif within their own communities. Empirical studies on synthetic (l, d) samples have shown that TFBSGroup is very efficient (e.g., it can find true (18, 6), (24, 8) motifs within 30 seconds). More importantly, the algorithm has succeeded in rapidly identifying motifs in a large data set of prokaryotic promoters generated from the Escherichia coli database RegulonDB. The algorithm has also accurately identified motifs in ChIP-seq data sets for 12 mouse transcription factors involved in ES cell pluripotency and self-renewal.

CONCLUSIONS

Our novel heuristic algorithm, TFBSGroup, is able to quickly identify nearly exact matches for long and weak (l, d) motifs in DNA sequences under the ZOMOPS constraint. It is also capable of finding motifs in real applications. The source code for TFBSGroup can be obtained from http://bioinformatics.bioengr.uic.edu/TFBSGroup/.

Genome-wide analysis of the salmonella Fis regulon and its regulatory mechanism on pathogenicity islands

Fis, one of the most important nucleoid-associated proteins, functions as a global regulator of transcription in bacteria that has been comprehensively studied in Escherichia coli K12. Fis also influences the virulence of Salmonella enterica and pathogenic E. coli by regulating their virulence genes, however, the relevant mechanism is unclear. In this report, using combined RNA-seq and chromatin immunoprecipitation (ChIP)-seq technologies, we first identified 1646 Fis-regulated genes and 885 Fis-binding targets in the S. enterica serovar Typhimurium, and found a Fis regulon different from that in E. coli. Fis has been reported to contribute to the invasion ability of S. enterica. By using cell infection assays, we found it also enhances the intracellular replication ability of S. enterica within macrophage cell, which is of central importance for the pathogenesis of infections. Salmonella pathogenicity islands (SPI)-1 and SPI-2 are crucial for the invasion and survival of S. enterica in host cells. Using mutation and overexpression experiments, real-time PCR analysis, and electrophoretic mobility shift assays, we demonstrated that Fis regulates 63 of the 94 Salmonella pathogenicity island (SPI)-1 and SPI-2 genes, by three regulatory modes: i) binds to SPI regulators in the gene body or in upstream regions; ii) binds to SPI genes directly to mediate transcriptional activation of themselves and downstream genes; iii) binds to gene encoding OmpR which affects SPI gene expression by controlling SPI regulators SsrA and HilD. Our results provide new insights into the impact of Fis on SPI genes and the pathogenicity of S. enterica.

Control of DNA minor groove width and Fis protein binding by the purine 2-amino group

The width of the DNA minor groove varies with sequence and can be a major determinant of DNA shape recognition by proteins. For example, the minor groove within the center of the Fis–DNA complex narrows to about half the mean minor groove width of canonical B-form DNA to fit onto the protein surface. G/C base pairs within this segment, which is not contacted by the Fis protein, reduce binding affinities up to 2000-fold over A/T-rich sequences. We show here through multiple X-ray structures and binding properties of Fis–DNA complexes containing base analogs that the 2-amino group on guanine is the primary molecular determinant controlling minor groove widths. Molecular dynamics simulations of free-DNA targets with canonical and modified bases further demonstrate that sequence-dependent narrowing of minor groove widths is modulated almost entirely by the presence of purine 2-amino groups. We also provide evidence that protein-mediated phosphate neutralization facilitates minor groove compression and is particularly important for binding to non-optimally shaped DNA duplexes.

E. coli Fis protein insulates the cbpA gene from uncontrolled transcription

The Escherichia coli curved DNA binding protein A (CbpA) is a poorly characterised nucleoid associated factor and co-chaperone. It is expressed at high levels as cells enter stationary phase. Using genetics, biochemistry, and genomics, we have examined regulation of, and DNA binding by, CbpA. We show that Fis, the dominant growth-phase nucleoid protein, prevents CbpA expression in growing cells. Regulation by Fis involves an unusual “insulation” mechanism. Thus, Fis protects cbpA from the effects of a distal promoter, located in an adjacent gene. In stationary phase, when Fis levels are low, CbpA binds the E. coli chromosome with a preference for the intrinsically curved Ter macrodomain. Disruption of the cbpA gene prompts dramatic changes in DNA topology. Thus, our work identifies a novel role for Fis and incorporates CbpA into the growing network of factors that mediate bacterial chromosome structure.

Compaction of chromosomal DNA is a fundamental process that impacts on all aspects of cellular biology. However, our understanding of chromosome organisation in bacteria is poorly developed. Since bacteria are amongst the most abundant living organisms on the planet, this represents a startling gap in our knowledge. Despite our lack of understanding, it has long been known that Escherichia coli, and other bacteria, radically re-model their chromosomes in response to environmental stress. This is most notable during periods of starvation, when the E. coli chromosome is super compacted. In dissecting the molecular mechanisms that control this phenomenon, we have found that regulatory cross-talk between DNA–organising proteins plays an essential role. Thus, the major DNA folding protein from growing E. coli inhibits production of the major chromosome organisers in starved cells. Our findings illustrate the highly dynamic nature of bacterial chromosomes. Thus, DNA topology, gene transcription, and chromosome folding proteins entwine to create a web of interactions that define the properties of the chromosome.

Binding of nucleoid-associated protein fis to DNA is regulated by DNA breathing dynamics

Physicochemical properties of DNA, such as shape, affect protein-DNA recognition. However, the properties of DNA that are most relevant for predicting the binding sites of particular transcription factors (TFs) or classes of TFs have yet to be fully understood. Here, using a model that accurately captures the melting behavior and breathing dynamics (spontaneous local openings of the double helix) of double-stranded DNA, we simulated the dynamics of known binding sites of the TF and nucleoid-associated protein Fis in Escherichia coli. Our study involves simulations of breathing dynamics, analysis of large published in vitro and genomic datasets, and targeted experimental tests of our predictions. Our simulation results and available in vitro binding data indicate a strong correlation between DNA breathing dynamics and Fis binding. Indeed, we can define an average DNA breathing profile that is characteristic of Fis binding sites. This profile is significantly enriched among the identified in vivo E. coli Fis binding sites. To test our understanding of how Fis binding is influenced by DNA breathing dynamics, we designed base-pair substitutions, mismatch, and methylation modifications of DNA regions that are known to interact (or not interact) with Fis. The goal in each case was to make the local DNA breathing dynamics either closer to or farther from the breathing profile characteristic of a strong Fis binding site. For the modified DNA segments, we found that Fis-DNA binding, as assessed by gel-shift assay, changed in accordance with our expectations. We conclude that Fis binding is associated with DNA breathing dynamics, which in turn may be regulated by various nucleotide modifications.

Cellular transcription factors (TFs) are proteins that regulate gene expression, and thereby cellular activity and fate, by binding to specific DNA segments. The physicochemical determinants of protein-DNA binding specificity are not completely understood. Here, we report that the propensity of transient opening and re-closing of the double helix, resulting from thermal fluctuations, aka “DNA breathing” or “DNA bubbles,” can be associated with binding affinity in the case of Fis, a well-studied nucleoid-associated protein in Escherichia coli. We found that a particular breathing profile is characteristic of high-affinity Fis binding sites and that DNA fragments known to bind Fis in vivo are statistically enriched for this profile. Furthermore, we used simulations of DNA breathing dynamics to guide design of gel-shift experiments aimed at testing the idea that local breathing influences Fis binding. As a result, we show that via nucleotide modifications but without modifying nucleotides that directly contact Fis, we were able to transform a low-affinity Fis binding site into a high-affinity site and vice versa. The nucleotide modifications were designed only based on DNA breathing simulations. Our study suggests that strong Fis-DNA binding depends on DNA breathing - a novel physicochemical characteristic that could be used for prediction and rational design of TF binding sites.

Chromatin organization, epigenetics and differentiation: an evolutionary perspective

Genome packaging is a universal phenomenon from prokaryotes to higher mammals. Genomic constituents and forces have however, travelled a long evolutionary route. Both DNA and protein elements constitute the genome and also aid in its dynamicity. With the evolution of organisms, these have experienced several structural and functional changes. These evolutionary changes were made to meet the challenging scenario of evolving organisms. This review discusses in detail the evolutionary perspective and functionality gain in the phenomena of genome organization and epigenetics.

Improved predictions of transcription factor binding sites using physicochemical features of DNA

Typical approaches for predicting transcription factor binding sites (TFBSs) involve use of a position-specific weight matrix (PWM) to statistically characterize the sequences of the known sites. Recently, an alternative physicochemical approach, called SiteSleuth, was proposed. In this approach, a linear support vector machine (SVM) classifier is trained to distinguish TFBSs from background sequences based on local chemical and structural features of DNA. SiteSleuth appears to generally perform better than PWM-based methods. Here, we improve the SiteSleuth approach by considering both new physicochemical features and algorithmic modifications. New features are derived from Gibbs energies of amino acid-DNA interactions and hydroxyl radical cleavage profiles of DNA. Algorithmic modifications consist of inclusion of a feature selection step, use of a nonlinear kernel in the SVM classifier, and use of a consensus-based post-processing step for predictions. We also considered SVM classification based on letter features alone to distinguish performance gains from use of SVM-based models versus use of physicochemical features. The accuracy of each of the variant methods considered was assessed by cross validation using data available in the RegulonDB database for 54 Escherichia coli TFs, as well as by experimental validation using published ChIP-chip data available for Fis and Lrp.

Regulation of the vapBC-1 toxin-antitoxin locus in nontypeable Haemophilus influenzae

Nontypeable Haemophilus influenzae (NTHi) are human-adapted commensal bacteria that can cause a number of chronic mucosal infections, including otitis media and bronchitis. One way for these organisms to survive antibiotic therapy and cause recurrent disease is to stop replicating, as most antimicrobials target essential biosynthetic pathways. Toxin-antitoxin (TA) gene pairs have been shown to facilitate entry into a reversible bacteriostatic state. Characteristically, these operons encode a protein toxin and an antitoxin that associate following translation to form a nontoxic complex, which then binds to and regulates the cognate TA promoter. Under stressful conditions, the labile antitoxin is degraded and the complex disintegrates, freeing the stable toxin to facilitate growth arrest. How these events affected the regulation of the TA locus, as well as how the transcription of the operon was subsequently returned to its normal state upon resumption of growth, was not fully understood. Here we show that expression of the NTHi vapBC-1 TA locus is repressed by a complex of VapB-1 and VapC-1 under conditions favorable for growth, and activated by the global transactivator Factor for Inversion Stimulation (Fis) upon nutrient upshift from stationary phase. Further, we demonstrate for the first time that the VapC-1 toxin alone can bind to its cognate TA locus control region and that the presence of VapB-1 directs the binding of the VapBC-1 complex in the transcriptional regulation of vapBC-1.

Silencing of toxic gene expression by Fis

Bacteria and bacteriophages have evolved DNA modification as a strategy to protect their genomes. Mom protein of bacteriophage Mu modifies the phage DNA, rendering it refractile to numerous restriction enzymes and in turn enabling the phage to successfully invade a variety of hosts. A strong fortification, a combined activity of the phage and host factors, prevents untimely expression of mom and associated toxic effects. Here, we identify the bacterial chromatin architectural protein Fis as an additional player in this crowded regulatory cascade. Both in vivo and in vitro studies described here indicate that Fis acts as a transcriptional repressor of mom promoter. Further, our data shows that Fis mediates its repressive effect by denying access to RNA polymerase at mom promoter. We propose that a combined repressive effect of Fis and previously characterized negative regulatory factors could be responsible to keep the gene silenced most of the time. We thus present a new facet of Fis function in Mu biology. In addition to bringing about overall downregulation of Mu genome, it also ensures silencing of the advantageous but potentially lethal mom gene.

"Anticipated" nucleosome positioning pattern in prokaryotes

Linguistic (word count) analysis of prokaryotic genome sequences, by Shannon N-gram extension, reveals that the dominant hidden motifs in A+T rich genomes are T(A)(T)A and G(A)(T)C with uncertain number of repeating A and T. Since prokaryotic sequences are largely protein-coding, the motifs would correspond to amphipathic alpha-helices with alternating lysine and phenylalanine as preferential polar and non-polar residues. The motifs are also known in eukaryotes, as nucleosome positioning patterns. Their existence in prokaryotes as well may serve for binding of histone-like proteins to DNA. In this case the above patterns in prokaryotes may be considered as "anticipated" nucleosome positioning patterns which, quite likely, existed in prokaryotic genomes before the evolutionary separation between eukaryotes and prokaryotes.

Force-driven unbinding of proteins HU and Fis from DNA quantified using a thermodynamic Maxwell relation

Determining numbers of proteins bound to large DNAs is important for understanding their chromosomal functions. Protein numbers may be affected by physical factors such as mechanical forces generated in DNA, e.g. by transcription or replication. We performed single-DNA stretching experiments with bacterial nucleoid proteins HU and Fis, verifying that the force-extension measurements were in thermodynamic equilibrium. We, therefore, could use a thermodynamic Maxwell relation to deduce the change of protein number on a single DNA due to varied force. For the binding of both HU and Fis under conditions studied, numbers of bound proteins decreased as force was increased. Our experiments showed that most of the bound HU proteins were driven off the DNA at 6.3 pN for HU concentrations lower than 150 nM; our HU data were fit well by a statistical-mechanical model of protein-induced bending of DNA. In contrast, a significant amount of Fis proteins could not be forced off the DNA at forces up to 12 pN and Fis concentrations up to 20 nM. This thermodynamic approach may be applied to measure changes in numbers of a wide variety of molecules bound to DNA or other polymers. Force-dependent DNA binding by proteins suggests mechano-chemical mechanisms for gene regulation.

Altered regulation of the OmpF porin by Fis in Escherichia coli during an evolution experiment and between B and K-12 strains

The phenotypic plasticity of global regulatory networks provides bacteria with rapid acclimation to a wide range of environmental conditions, while genetic changes in those networks provide additional flexibility as bacteria evolve across long time scales. We previously identified mutations in the global regulator-encoding gene fis that enhanced organismal fitness during a long-term evolution experiment with Escherichia coli. To gain insight into the effects of these mutations, we produced two-dimensional protein gels with strains carrying different fis alleles, including a beneficial evolved allele and one with an in-frame deletion. We found that Fis controls the expression of the major porin-encoding gene ompF in the E. coli B-derived ancestral strain used in the evolution experiment, a relationship that has not been described before. We further showed that this regulatory connection evolved over two different time scales, perhaps explaining why it was not observed before. On the longer time scale, we showed that this regulation of ompF by Fis is absent from the more widely studied K-12 strain and thus is specific to the B strain. On a shorter time scale, this regulatory linkage was lost during 20,000 generations of experimental evolution of the B strain. Finally, we mapped the Fis binding sites in the ompF regulatory region, and we present a hypothetical model of ompF expression that includes its other known regulators.

Concentration-dependent exchange accelerates turnover of proteins bound to double-stranded DNA

The multistep kinetics through which DNA-binding proteins bind their targets are heavily studied, but relatively little attention has been paid to proteins leaving the double helix. Using single-DNA stretching and fluorescence detection, we find that sequence-neutral DNA-binding proteins Fis, HU and NHP6A readily exchange with themselves and with each other. In experiments focused on the Escherichia coli nucleoid-associated protein Fis, only a small fraction of protein bound to DNA spontaneously dissociates into protein-free solution. However, if Fis is present in solution, we find that a concentration-dependent exchange reaction occurs which turns over the bound protein, with a rate of k_exch = 6 × 10⁴ M⁻¹s⁻¹. The bacterial DNA-binding protein HU and the yeast HMGB protein NHP6A display the same phenomenon of protein in solution accelerating dissociation of previously bound labeled proteins as exchange occurs. Thus, solvated proteins can play a key role in facilitating removal and renewal of proteins bound to the double helix, an effect that likely plays a major role in promoting the turnover of proteins bound to DNA in vivo and, therefore, in controlling the dynamics of gene regulation.

A brief review of molecular information theory

The idea that we could build molecular communications systems can be advanced by investigating how actual molecules from living organisms function. Information theory provides tools for such an investigation. This review describes how we can compute the average information in the DNA binding sites of any genetic control protein and how this can be extended to analyze its individual sites. A formula equivalent to Claude Shannon's channel capacity can be applied to molecular systems and used to compute the efficiency of protein binding. This efficiency is often 70% and a brief explanation for that is given. The results imply that biological systems have evolved to function at channel capacity, which means that we should be able to build molecular communications that are just as robust as our macroscopic ones.

The shape of the DNA minor groove directs binding by the DNA-bending protein Fis

The bacterial nucleoid-associated protein Fis regulates diverse reactions by bending DNA and through DNA-dependent interactions with other control proteins and enzymes. In addition to dynamic nonspecific binding to DNA, Fis forms stable complexes with DNA segments that share little sequence conservation. Here we report the first crystal structures of Fis bound to high- and low-affinity 27-base-pair DNA sites. These 11 structures reveal that Fis selects targets primarily through indirect recognition mechanisms involving the shape of the minor groove and sequence-dependent induced fits over adjacent major groove interfaces. The DNA shows an overall curvature of approximately 65 degrees , and the unprecedented close spacing between helix-turn-helix motifs present in the apodimer is accommodated by severe compression of the central minor groove. In silico DNA structure models show that only the roll, twist, and slide parameters are sufficient to reproduce the changes in minor groove widths and recreate the curved Fis-bound DNA structure. Models based on naked DNA structures suggest that Fis initially selects DNA targets with intrinsically narrow minor grooves using the separation between helix-turn-helix motifs in the Fis dimer as a ruler. Then Fis further compresses the minor groove and bends the DNA to generate the bound structure.

Role of decreased levels of Fis histone-like protein in Crohn's disease-associated adherent invasive Escherichia coli LF82 bacteria interacting with intestinal epithelial cells

The interaction of Crohn's disease (CD)-associated adherent-invasive Escherichia coli (AIEC) strain LF82 with intestinal epithelial cells depends on surface appendages, such as type 1 pili and flagella. Histone-like proteins operate as global regulators to control the expression of these virulence factors. We evaluated the role of histone-like proteins in AIEC reference strain LF82 during infection of intestinal epithelial cells, Intestine-407, and observed that the fis mRNA level was decreased. The role of Fis in AIEC LF82 was determined by studying the phenotype of an LF82 fis::Km mutant. This was the first mutant of strain LF82 that has been described thus far that is unable to express flagellin but still able to produce type 1 pili. The cyclic-di-GMP pathway linking flagella and type 1 pilus expression is not involved in Fis-mediated regulation, and we identified in the present study Fis-binding sites located upstream of the fimE gene and in the intergenic region between fimB and nanC of the fim operon encoding type 1 pili. The major consequence of decreased Fis expression in AIEC bacteria in contact with host cells is a direct downregulation of fimE expression, leading to the preferential ON phase of the fimS element. Thus, by maintaining type 1 pilus expression, AIEC bacteria, which interact with the gut mucosa, have greater ability to colonize and to induce inflammation in CD patients.

The role of FIS in the Rcd checkpoint and stable maintenance of plasmid ColE1

Escherichia coli plasmid ColE1 lacks active partitioning, and copies are distributed randomly to daughter cells at division. The plasmid is maintained stably in the bacterial population as long as its copy number remains high. The accumulation of plasmid dimers and higher multimers depresses copy number, and is an important cause of multicopy plasmid instability. ColE1 dimers are restored to the monomeric state by site-specific recombination, which requires the host-encoded proteins XerCD, ArgR and PepA acting at the plasmid cer site. In addition, a 70 nt RNA expressed from the cer site of plasmid dimers delays the division of dimer-containing cells. Here, we report that the global regulator FIS binds to cer in a sequence-specific manner, close to the Rcd promoter (P cer ). FIS is not required for plasmid dimer resolution, but is essential for repression of P cer in plasmid monomers. Repression also requires the XerCD recombinase, but not ArgR or PepA. We propose a model for monomer–dimer control of P cer in which the promoter is repressed in plasmid monomers by the concerted action of FIS and XerCD. Rcd transcription is triggered in plasmid dimers by the lifting of XerCD-mediated repression in the synaptic complex.

Fis negatively affects binding of Tn4652 transposase by out-competing IHF from the left end of Tn4652

Transposition activity in bacteria is generally maintained at a low level. The activity of mobile DNA elements can be controlled by bacterially encoded global regulators. Regulation of transposition of Tn4652 in Pseudomonas putida is one such example. Activation of transposition of Tn4652 in starving bacteria requires the stationary-phase sigma factor RpoS and integration host factor (IHF). IHF plays a dual role in Tn4652 translocation by activating transcription of the transposase gene tnpA of the transposon and facilitating TnpA binding to the inverted repeats of the transposon. Our previous results have indicated that besides IHF some other P. putida-encoded global regulator(s) might bind to the ends of Tn4652 and regulate transposition activity. In this study, employing a DNase I footprint assay we have identified a binding site of P. putida Fis (factor for inversion stimulation) centred 135 bp inside the left end of Tn4652. Our results of gel mobility shift and DNase I footprint studies revealed that Fis out-competes IHF from the left end of Tn4652, thereby abolishing the binding of TnpA. Thus, the results obtained in this study indicate that the transposition of Tn4652 is regulated by the cellular amount of P. putida global regulators Fis and IHF.

The role of FIS protein in the physiological control of the expression of the Escherichia coli meta-hpa operon

Expression from the Escherichia coli W meta-hpa operon promoter (Pg) is under a strict catabolic repression control mediated by the cAMP-catabolite repression protein (CRP) complex in a glucose-containing medium. The Pg promoter is also activated by the integration host factor (IHF) and repressed by the specific transcriptional regulator HpaR when 4-hydroxyphenylacetate (4HPA) is not present in the medium. Expression from the hpa promoter is also repressed in undefined rich medium such as LB, but the molecular basis of this mechanism is not understood. We present in vitro and in vivo studies to demonstrate the involvement of FIS protein in this catabolic repression. DNase I footprinting experiments show that FIS binds to multiple sites within the Pg promoter. FIS-site I overlaps the CRP-binding site. By using an electromobility shift assay, we demonstrated that FIS efficiently competes with CRP for binding to the Pg promoter, suggesting an antagonist/competitive mechanism. RT-PCR showed that the Pg repression effect is relieved in a FIS deleted strain. The repression role of FIS at Pg was further demonstrated by in vitro transcription assays. These results suggest that FIS contributes to silencing the Pg promoter in the exponential phase of growth in an undefined rich medium when FIS is predominantly expressed.

Downregulation of the Escherichia coli guaB promoter by FIS

The Escherichia coli guaB promoter (P_guaB) regulates transcription of two genes, guaB and guaA, that are required for the synthesis of guanosine 5′-monophosphate (GMP), a precursor for the synthesis of guanine nucleoside triphosphates. Transcription from P_guaB increases as a function of increasing cellular growth rate, and this is referred to as growth rate-dependent control (GRDC). Here we investigated the role of the factor for inversion stimulation (FIS) in the regulation of this promoter. The results showed that there are three binding sites for FIS centred near positions −11, +8 and +29 relative to the guaB transcription start site. Binding of FIS to these sites results in repression of P_guaB in vitro but not in vivo. Deletion of the fis gene results in increased P_guaB activity in vivo, but GRDC of P_guaB is maintained.

Biochemical identification of base and phosphate contacts between Fis and a high-affinity DNA binding site

Fis (factor for inversion stimulation) is a nucleoid-associated protein in Escherichia coli and other bacteria that stimulates certain site-specific DNA recombination events, alters DNA topology, and serves as a global gene regulator. DNA binding is central to the functions of Fis and involves a helix-turn-helix DNA binding motif located in the carboxy-terminal region. Specific DNA binding is observed at a number of sites exhibiting poorly related sequences. Such interactions require four critical base pairs positioned -7, -3, +3, and +7 nucleotides relative to the central nucleotide of a 15-bp core-binding site. To further understand how Fis interacts with DNA, we identified the positions of 14 DNA phosphates (based on ethylation interference assays) that are required for Fis binding. These are the 5' phosphates of the nucleotides at positions -8, -7, -6, +1, +2, +3, and +4 relative to the central nucleotide on both DNA strands. Another five phosphates located in the flanking regions from positions +10 through +14 can serve as additional contact sites. Using a combination of biochemical approaches and various mutant Fis proteins, we probed possible interactions between several key Fis residues and DNA bases or phosphates within a high-affinity binding site. We provide evidence in support of interactions between the R85 Fis residue and a highly conserved guanine at position -7 and between T87 and the critical base pairs at -3 and +3. In addition, we present evidence in support of interactions between N84 and the phosphate 5' to the base at +4, between R89 and the -7 phosphate, between T87 and the +3 and +4 phosphates, and between K90 and the +3 phosphate. This work provides functional evidence for some of the most critical interactions between Fis and DNA required for a high binding affinity and demonstrates the large contribution made by numerous phosphates to the stability of the Fis-DNA complex.

Selective repression by Fis and H-NS at the Escherichia coli dps promoter

Dps is a nucleoid-associated protein that plays a major role in condensation of the Escherichia coli chromosome in stationary phase. Here we show that two other nucleoid-associated proteins, Fis and H-NS, can bind at the dps gene promoter and downregulate its activity. Both Fis and H-NS selectively repress the dps promoter, preventing transcription initiation by RNA polymerase containing sigma(70), the housekeeping sigma factor, but not by RNA polymerase containing sigma(38), the stationary-phase sigma factor. Fis represses by trapping RNA polymerase containing sigma(70) at the promoter. In contrast, H-NS functions by displacing RNA polymerase containing sigma(70), but not RNA polymerase containing sigma(38). Dps levels are known to be very low in exponentially growing cells and rise sharply as cells enter stationary phase. Conversely, Fis levels are high in growing cells but fall to nearly zero in stationary-phase cells. Our data suggest a simple model to explain how the Dps-dependent super-compaction of the folded chromosome is triggered as cell growth ceases.

Genome-wide analysis of Fis binding in Escherichia coli indicates a causative role for A-/AT-tracts

We determined the genome-wide distribution of the nucleoid-associated protein Fis in Escherichia coli using chromatin immunoprecipitation coupled with high-resolution whole genome-tiling microarrays. We identified 894 Fis-associated regions across the E. coli genome. A significant number of these binding sites were found within open reading frames (33%) and between divergently transcribed transcripts (5%). Analysis indicates that A-tracts and AT-tracts are an important signal for preferred Fis-binding sites, and that A(6)-tracts in particular constitute a high-affinity signal that dictates Fis phasing in stretches of DNA containing multiple and variably spaced A-tracts and AT-tracts. Furthermore, we find evidence for an average of two Fis-binding regions per supercoiling domain in the chromosome of exponentially growing cells. Transcriptome analysis shows that approximately 21% of genes are affected by the deletion of fis; however, the changes in magnitude are small. To address the differential Fis bindings under growth environment perturbation, ChIP-chip analysis was performed using cells grown under aerobic and anaerobic growth conditions. Interestingly, the Fis-binding regions are almost identical in aerobic and anaerobic growth conditions-indicating that the E. coli genome topology mediated by Fis is superficially identical in the two conditions. These novel results provide new insight into how Fis modulates DNA topology at a genome scale and thus advance our understanding of the architectural bases of the E. coli nucleoid.

Functional characterization of the Escherichia coli Fis-DNA binding sequence

The Escherichia coli protein Fis is remarkable for its ability to interact specifically with DNA sites of highly variable sequences. The mechanism of this sequence-flexible DNA recognition is not well understood. In a previous study, we examined the contributions of Fis residues to high-affinity binding at different DNA sequences using alanine-scanning mutagenesis and identified several key residues for Fis-DNA recognition. In this work, we investigated the contributions of the 15-bp core Fis binding sequence and its flanking regions to Fis-DNA interactions. Systematic base-pair replacements made in both half sites of a palindromic Fis binding sequence were examined for their effects on the relative Fis binding affinity. Missing contact assays were also used to examine the effects of base removal within the core binding site and its flanking regions on the Fis-DNA binding affinity. The results revealed that: (1) the -7G and +3Y bases in both DNA strands (relative to the central position of the core binding site) are major determinants for high-affinity binding; (2) the C(5) methyl group of thymine, when present at the +4 position, strongly hinders Fis binding; and (3) AT-rich sequences in the central and flanking DNA regions facilitate Fis-DNA interactions by altering the DNA structure and by increasing the local DNA flexibility. We infer that the degeneracy of specific Fis binding sites results from the numerous base-pair combinations that are possible at noncritical DNA positions (from -6 to -4, from -2 to +2, and from +4 to +6), with only moderate penalties on the binding affinity, the roughly similar contributions of -3A or G and +3T or C to the binding affinity, and the minimal requirement of three of the four critical base pairs to achieve considerably high binding affinities.

Effects of Fis on Escherichia coli gene expression during different growth stages

Fis is a nucleoid-associated protein in Escherichia coli that is abundant during early exponential growth in rich medium but is in short supply during stationary phase. Its role as a transcriptional regulator has been demonstrated for an increasing number of genes. In order to gain insight into the global effects of Fis on E. coli gene expression during different stages of growth in rich medium, DNA microarray analyses were conducted in fis and wild-type strains during early, mid-, late-exponential and stationary growth phases. The results uncovered 231 significantly regulated genes that were distributed over 15 functional categories. Regulatory effects were observed at all growth stages examined. Coordinate upregulation was observed for a number of genes involved in translation, flagellar biosynthesis and motility, nutrient transport, carbon compound metabolism, and energy metabolism at different growth stages. Coordinate down-regulation was also observed for genes involved in stress response, amino acid and nucleotide biosynthesis, energy and intermediary metabolism, and nutrient transport. As cells transitioned from the early to the late-exponential growth phase, different functional categories of genes were regulated, and a gradual shift occurred towards mostly down-regulation. The results demonstrate that the growth phase-dependent Fis expression triggers coordinate regulation of 15 categories of functionally related genes during specific stages of growth of an E. coli culture.

Discovery of Fur binding site clusters in Escherichia coli by information theory models

Fur is a DNA binding protein that represses bacterial iron uptake systems. Eleven footprinted Escherichia coli Fur binding sites were used to create an initial information theory model of Fur binding, which was then refined by adding 13 experimentally confirmed sites. When the refined model was scanned across all available footprinted sequences, sequence walkers, which are visual depictions of predicted binding sites, frequently appeared in clusters that fit the footprints ( approximately 83% coverage). This indicated that the model can accurately predict Fur binding. Within the clusters, individual walkers were separated from their neighbors by exactly 3 or 6 bases, consistent with models in which Fur dimers bind on different faces of the DNA helix. When the E. coli genome was scanned, we found 363 unique clusters, which includes all known Fur-repressed genes that are involved in iron metabolism. In contrast, only a few of the known Fur-activated genes have predicted Fur binding sites at their promoters. These observations suggest that Fur is either a direct repressor or an indirect activator. The Pseudomonas aeruginosa and Bacillus subtilis Fur models are highly similar to the E. coli Fur model, suggesting that the Fur-DNA recognition mechanism may be conserved for even distantly related bacteria.

Expression of the Fis protein is sustained in late-exponential- and stationary-phase cultures of Salmonella enterica serovar Typhimurium grown in the absence of aeration

The classic expression pattern of the Fis global regulatory protein during batch culture consists of a high peak in the early logarithmic phase of growth, followed by a sharp decrease through mid-exponential growth phase until Fis is almost undetectable at the end of the exponential phase. We discovered that this pattern is contingent on the growth regime. In Salmonella enterica serovar Typhimurium cultures grown in non-aerated SPI1-inducing conditions, Fis can be detected readily in stationary phase. On the other hand, cultures grown with standard aeration showed the classic Fis expression pattern. Sustained Fis expression in non-aerated cultures was also detected in some Escherichia coli strains, but not in others. This novel pattern of Fis expression was independent of sequence differences in the fis promoter regions of Salmonella and E. coli. Instead, a clear negative correlation between the expression of the Fis protein and of the stress-and-stationary-phase sigma factor RpoS was observed in a variety of strains. An rpoS mutant displayed elevated levels of Fis and had a higher frequency of epithelial cell invasion under these growth conditions. We discuss a model whereby Fis and RpoS levels vary in response to environmental signals allowing the expression and repression of SPI1 invasion genes.

Correlation between binding rate constants and individual information of E. coli Fis binding sites

Individual protein binding sites on DNA can be measured in bits of information. This information is related to the free energy of binding by the second law of thermodynamics, but binding kinetics appear to be inaccessible from sequence information since the relative contributions of the on- and off-rates to the binding constant, and hence the free energy, are unknown. However, the on-rate could be independent of the sequence since a protein is likely to bind once it is near a site. To test this, we used surface plasmon resonance and electromobility shift assays to determine the kinetics for binding of the Fis protein to a range of naturally occurring binding sites. We observed that the logarithm of the off-rate is indeed proportional to the individual information of the binding sites, as predicted. However, the on-rate is also related to the information, but to a lesser degree. We suggest that the on-rate is mostly determined by DNA bending, which in turn is determined by the sequence information. Finally, we observed a break in the binding curve around zero bits of information. The break is expected from information theory because it represents the coding demarcation between specific and nonspecific binding.

Bioinformatic prediction and experimental verification of Fur-regulated genes in the extreme acidophile Acidithiobacillus ferrooxidans

The gamma-proteobacterium Acidithiobacillus ferrooxidans lives in extremely acidic conditions (pH 2) and, unlike most organisms, is confronted with an abundant supply of soluble iron. It is also unusual in that it oxidizes iron as an energy source. Consequently, it faces the challenging dual problems of (i) maintaining intracellular iron homeostasis when confronted with extremely high environmental loads of iron and (ii) of regulating the use of iron both as an energy source and as a metabolic micronutrient. A combined bioinformatic and experimental approach was undertaken to identify Fur regulatory sites in the genome of A. ferrooxidans and to gain insight into the constitution of its Fur regulon. Fur regulatory targets associated with a variety of cellular functions including metal trafficking (e.g. feoPABC, tdr, tonBexbBD, copB, cdf), utilization (e.g. fdx, nif), transcriptional regulation (e.g. phoB, irr, iscR) and redox balance (grx, trx, gst) were identified. Selected predicted Fur regulatory sites were confirmed by FURTA, EMSA and in vitro transcription analyses. This study provides the first model for a Fur-binding site consensus sequence in an acidophilic iron-oxidizing microorganism and lays the foundation for future studies aimed at deepening our understanding of the regulatory networks that control iron uptake, homeostasis and oxidation in extreme acidophiles.