Stress-induced DNA duplex destabilization (SIDD) in the E. coli genome: SIDD sites are closely associated with promoters

DNA superhelicity

Closing each strand of a DNA duplex upon itself fixes its linking number L. This topological condition couples together the secondary and tertiary structures of the resulting ccDNA topoisomer, a constraint that is not present in otherwise identical nicked or linear DNAs. Fixing L has a range of structural, energetic and functional consequences. Here we consider how L having different integer values (that is, different superhelicities) affects ccDNA molecules. The approaches used are primarily theoretical, and are developed from a historical perspective. In brief, processes that either relax or increase superhelicity, or repartition what is there, may either release or require free energy. The energies involved can be substantial, sufficient to influence many events, directly or indirectly. Here two examples are developed. The changes of unconstrained superhelicity that occur during nucleosome attachment and release are examined. And a simple theoretical model of superhelically driven DNA structural transitions is described that calculates equilibrium distributions for populations of identical topoisomers. This model is used to examine how these distributions change with superhelicity and other factors, and applied to analyze several situations of biological interest.

DNA supercoiling in bacteria: state of play and challenges from a viewpoint of physics based modeling

DNA supercoiling is central to many fundamental processes of living organisms. Its average level along the chromosome and over time reflects the dynamic equilibrium of opposite activities of topoisomerases, which are required to relax mechanical stresses that are inevitably produced during DNA replication and gene transcription. Supercoiling affects all scales of the spatio-temporal organization of bacterial DNA, from the base pair to the large scale chromosome conformation. Highlighted in vitro and in vivo in the 1960s and 1970s, respectively, the first physical models were proposed concomitantly in order to predict the deformation properties of the double helix. About fifteen years later, polymer physics models demonstrated on larger scales the plectonemic nature and the tree-like organization of supercoiled DNA. Since then, many works have tried to establish a better understanding of the multiple structuring and physiological properties of bacterial DNA in thermodynamic equilibrium and far from equilibrium. The purpose of this essay is to address upcoming challenges by thoroughly exploring the relevance, predictive capacity, and limitations of current physical models, with a specific focus on structural properties beyond the scale of the double helix. We discuss more particularly the problem of DNA conformations, the interplay between DNA supercoiling with gene transcription and DNA replication, its role on nucleoid formation and, finally, the problem of scaling up models. Our primary objective is to foster increased collaboration between physicists and biologists. To achieve this, we have reduced the respective jargon to a minimum and we provide some explanatory background material for the two communities.

Conserved microRNAs and Flipons Shape Gene Expression during Development by Altering Promoter Conformations

The classical view of gene regulation draws from prokaryotic models, where responses to environmental changes involve operons regulated by sequence-specific protein interactions with DNA, although it is now known that operons are also modulated by small RNAs. In eukaryotes, pathways based on microRNAs (miR) regulate the readout of genomic information from transcripts, while alternative nucleic acid structures encoded by flipons influence the readout of genetic programs from DNA. Here, we provide evidence that miR- and flipon-based mechanisms are deeply connected. We analyze the connection between flipon conformation and the 211 highly conserved human miR that are shared with other placental and other bilateral species. The direct interaction between conserved miR (c-miR) and flipons is supported by sequence alignments and the engagement of argonaute proteins by experimentally validated flipons as well as their enrichment in promoters of coding transcripts important in multicellular development, cell surface glycosylation and glutamatergic synapse specification with significant enrichments at false discovery rates as low as 10^-116. We also identify a second subset of c-miR that targets flipons essential for retrotransposon replication, exploiting that vulnerability to limit their spread. We propose that miR can act in a combinatorial manner to regulate the readout of genetic information by specifying when and where flipons form non-B DNA (NoB) conformations, providing the interactions of the conserved hsa-miR-324-3p with RELA and the conserved hsa-miR-744 with ARHGAP5 genes as examples.

The Impact of RNA modifications on the Biology of DNA Virus Infection

Approximately 170 RNA modifications have been identified and these are critical for determining the fate and function of cellular RNAs. Similar to human transcripts, viral RNAs possess an extensive RNA modification landscape. While initial efforts largely focused on investigating the RNA modification landscape in the context of RNA virus infection, a growing body of work has explored the impact of RNA modifications on DNA virus biology. These studies have revealed roles for RNA modifications in DNA virus infection, including gene regulation and viral pathogenesis. In this review, we will discuss the current knowledge on how RNA modifications impact DNA virus biology.

Biochemical and structural basis for YTH domain of human YTHDC1 binding to methylated adenine in DNA

The recently characterized mammalian writer (methyltransferase) and eraser (demethylase) of the DNA N6-methyladenine (N6mA) methyl mark act on single-stranded (ss) and transiently-unpaired DNA. As YTH domain-containing proteins bind N6mA-containing RNA in mammalian cells, we investigated whether mammalian YTH domains are also methyl mark readers of N6mA DNA. Here, we show that the YTH domain of YTHDC1 (known to localize in the nucleus) binds ssDNA containing N6mA, with a 10 nM dissociation constant. This binding is stronger by a factor of 5 than in an RNA context, tested under the same conditions. However, the YTH domains of YTHDF2 and YTHDF1 (predominantly cytoplasmic) exhibited the opposite effect with ∼1.5-2× stronger binding to ssRNA containing N6mA than to the corresponding DNA. We determined two structures of the YTH domain of YTHDC1 in complex with N6mA-containing ssDNA, which illustrated that YTHDC1 binds the methylated adenine in a single-stranded region flanked by duplexed DNA. We discuss the hypothesis that the writer-reader-eraser of N6mA-containining ssDNA is associated with maintaining genome stability. Structural comparison of YTH and SRA domains (the latter a DNA 5-methylcytosine reader) revealed them to be diverse members of a larger family of DNA/RNA modification readers, apparently having originated from bacterial modification-dependent restriction enzymes.

How global DNA unwinding causes non-uniform stress distribution and melting of DNA

DNA unwinding is an important process that controls binding of proteins, gene expression and melting of double-stranded DNA. In a series of all-atom MD simulations on two DNA molecules containing a transcription start TATA-box sequence we demonstrate that application of a global restraint on the DNA twisting dramatically changes the coupling between helical parameters and the distribution of deformation energy along the sequence. Whereas only short range nearest-neighbor coupling is observed in the relaxed case, long-range coupling is induced in the globally restrained case. With increased overall unwinding the elastic deformation energy is strongly non-uniformly distributed resulting ultimately in a local melting transition of only the TATA box segment during the simulations. The deformation energy tends to be stored more in cytidine/guanine rich regions associated with a change in conformational substate distribution. Upon TATA box melting the deformation energy is largely absorbed by the melting bubble with the rest of the sequences relaxing back to near B-form. The simulations allow us to characterize the structural changes and the propagation of the elastic energy but also to calculate the associated free energy change upon DNA unwinding up to DNA melting. Finally, we design an Ising model for predicting the local melting transition based on empirical parameters. The direct comparison with the atomistic MD simulations indicates a remarkably good agreement for the predicted necessary torsional stress to induce a melting transition, for the position and length of the melted region and for the calculated associated free energy change between both approaches.

Whole genome sequence of moderate halophilic marine bacterium Marinobacter litoralis SW-45: Abundance of non-coding RNAs

A report on the de novo Whole Genome Sequence (WGS) of Marinobacter litoralis SW-45, a moderately salt-tolerant bacterium isolated from the seawater in Malaysia is presented. The strain has a genome size of 3.45 Mb and is capable of producing halophilic lipase, protease and esterase enzymes. Computational prediction of non-coding RNA (ncRNA) genes in M. litoralis SW-45 was performed using standalone software known as the non-coding RNA characterization (nocoRNAc). In addition, a phylogenetic tree showing the evolutionary relationship between the strain and other members of the genus Marinobacter was constructed using 16SrRNA sequence information. A total of 385 ncRNA transcripts, 1124 terminator region, and 2350 Stress Induced Duplex Destabilization sites were predicted. The current WGS shotgun project has provided the relevant genetic information that may be useful for the strain's improvement studies. This manuscript gives the first description of M. litoralis with a complete genome.

Hopping into a hot seat: Role of DNA structural features on IS5-mediated gene activation and inactivation under stress

Insertion sequence elements (IS elements) are proposed to play major roles in shaping the genetic and phenotypic landscapes of prokaryotic cells. Recent evidence has raised the possibility that environmental stress conditions increase IS hopping into new sites, and often such hopping has the phenotypic effect of relieving the stress. Although stress-induced targeted mutations have been reported for a number of E. coli genes, the glpFK (glycerol utilization) and the cryptic bglGFB (β-glucoside utilization) systems are among the best characterized where the effects of IS insertion-mediated gene activation are well-characterized at the molecular level. In the glpFK system, starvation of cells incapable of utilizing glycerol leads to an IS5 insertion event that activates the glpFK operon, and enables glycerol utilization. In the case of the cryptic bglGFB operon, insertion of IS5 (and other IS elements) into a specific region in the bglG upstream sequence has the effect of activating the operon in both growing cells, and in starving cells. However, a major unanswered question in the glpFK system, the bgl system, as well as other examples, has been why the insertion events are promoted at specific locations, and how the specific stress condition (glycerol starvation for example) can be mechanistically linked to enhanced insertion at a specific locus. In this paper, we show that a specific DNA structural feature (superhelical stress-induced duplex destabilization, SIDD) is associated with "stress-induced" IS5 insertion in the glpFK, bglGFB, flhDC, fucAO and nfsB systems. We propose a speculative mechanistic model that links specific environmental conditions to the unmasking of an insertional hotspot in the glpFK system. We demonstrate that experimentally altering the predicted stability of a SIDD element in the nfsB gene significantly impacts IS5 insertion at its hotspot.

Environment-directed activation of the Escherichia coliflhDC operon by transposons

The flagellar system in Escherichia coli K12 is expressed under the control of the flhDC-encoded master regulator FlhDC. Transposition of insertion sequence (IS) elements to the upstream flhDC promoter region up-regulates transcription of this operon, resulting in a more rapid motility. Wang and Wood (ISME J 2011;5:1517-1525) provided evidence that insertion of IS5 into upstream activating sites occurs at higher rates in semi-solid agar media in which swarming behaviour is allowed as compared with liquid or solid media where swarming cannot occur. We confirm this conclusion and show that three IS elements, IS1, IS3 and IS5, transpose to multiple upstream sites within a 370 bp region of the flhDC operon control region. Hot spots for IS insertion correlate with positions of stress-induced DNA duplex destabilization (SIDD). We show that IS insertion occurs at maximal rates in 0.24 % agar, with rates decreasing dramatically with increasing or decreasing agar concentrations. In mixed cultures, we show that these mutations preferentially arise from the wild-type parent at frequencies of up to 3×10-3 cell-1 day-1 when the inoculated parental and co-existing IS-activated mutant cells are entering the stationary growth phase. We rigorously show that the apparent increased mutation frequencies cannot be accounted for by increased swimming or by increased growth under the selective conditions used. Thus, our data are consistent with the possibility that appropriate environmental conditions, namely those that permit but hinder flagellar rotation, result in the activation of a mutational pathway that involves IS element insertion upstream of the flhDC operon.

Permanganate/S1 Nuclease Footprinting Reveals Non-B DNA Structures with Regulatory Potential across a Mammalian Genome

DNA in cells is predominantly B-form double helix. Though certain DNA sequences in vitro may fold into other structures, such as triplex, left-handed Z form, or quadruplex DNA, the stability and prevalence of these structures in vivo are not known. Here, using computational analysis of sequence motifs, RNA polymerase II binding data, and genome-wide potassium permanganate-dependent nuclease footprinting data, we map thousands of putative non-B DNA sites at high resolution in mouse B cells. Computational analysis associates these non-B DNAs with particular structures and indicates that they form at locations compatible with an involvement in gene regulation. Further analyses support the notion that non-B DNA structure formation influences the occupancy and positioning of nucleosomes in chromatin. These results suggest that non-B DNAs contribute to the control of a variety of critical cellular and organismal processes.

Significant Role of Segmental Duplications and SIDD Sites in Chromosomal Translocations of Hematological Malignancies: A Multi-parametric Bioinformatic Analysis

Recurrent non-random chromosomal translocations are hallmark characteristics of leukemogenesis, and however, molecular mechanisms underlying these rearrangements are less explored. The fundamental question is, why and how chromosomes break and reunite so precisely in the genome. Meticulous understanding of mechanism leading to chromosomal rearrangement can be achieved by characterizing breakpoints. To address this hypothesis, a novel multi-parametric computational approach for characterization of major leukemic translocations within and around breakpoint region was performed. To best of our knowledge, this bioinformatic analysis is unique in finding the presence of segmental duplications (SDs) flanking breakpoints of all major leukemic translocation. Breakpoint islands (BpIs) were analyzed for stress-induced duplex destabilization (SIDD) sites along with other complex genomic architecture and physicochemical properties. Our study distinctly emphasizes on the probable correlative role of SDs, SIDD sites and various genomic features in the occurrence of breakpoints. Further, it also highlights potential features which may be playing a crucial role in causing double-strand breaks, leading to translocation.

Long-range correlations in the mechanics of small DNA circles under topological stress revealed by multi-scale simulation

It is well established that gene regulation can be achieved through activator and repressor proteins that bind to DNA and switch particular genes on or off, and that complex metabolic networks determine the levels of transcription of a given gene at a given time. Using three complementary computational techniques to study the sequence-dependence of DNA denaturation within DNA minicircles, we have observed that whenever the ends of the DNA are constrained, information can be transferred over long distances directly by the transmission of mechanical stress through the DNA itself, without any requirement for external signalling factors. Our models combine atomistic molecular dynamics (MD) with coarse-grained simulations and statistical mechanical calculations to span three distinct spatial resolutions and timescale regimes. While they give a consensus view of the non-locality of sequence-dependent denaturation in highly bent and supercoiled DNA loops, each also reveals a unique aspect of long-range informational transfer that occurs as a result of restraining the DNA within the closed loop of the minicircles.

Characterization of chromosomal translocation breakpoint sequences in solid tumours: "an in silico analysis"

Chromosomal translocations that results in formation and activation of fusion oncogenes are observed in numerous solid malignancies since years back. Expression of fusion kinases in these cancers drives the initiation & progression that ultimately leads to tumour development and thus comes out to be clinically imperative in terms of diagnosis and treatment of cancer. Nonetheless, molecular mechanisms beneath these translocations remained unexplored consequently limiting our knowledge of carcinogenesis and hence is the current field where further research is required. The issue of prime focus is the precision with which the chromosomes breaks and reunites within genome. Characterization of Genomic sequences located at Breakpoint region may direct us towards the thorough understanding of mechanism leading to chromosomal rearrangement. A unique computational multi-parametric analysis was performed for characterization of genomic sequence within and around breakpoint region. This study turns out to be novel as it reveals the occurrence of Segmental Duplications flanking the breakpoints of all translocation. Breakpoint Islands were also investigated for the presence of other intricate genomic architecture and various physico-chemical parameters. Our study particularly highlights the probable role of SDs and specific genomic features in precise chromosomal breakage. Additionally, it pinpoints the potential features that may be significant for double-strand breaks leading to chromosomal rearrangements.

Distinguishing between productive and abortive promoters using a random forest classifier in Mycoplasma pneumoniae

Distinguishing between promoter-like sequences in bacteria that belong to true or abortive promoters, or to those that do not initiate transcription at all, is one of the important challenges in transcriptomics. To address this problem, we have studied the genome-reduced bacterium Mycoplasma pneumoniae, for which the RNAs associated with transcriptional start sites have been recently experimentally identified. We determined the contribution to transcription events of different genomic features: the –10, extended –10 and –35 boxes, the UP element, the bases surrounding the –10 box and the nearest-neighbor free energy of the promoter region. Using a random forest classifier and the aforementioned features transformed into scores, we could distinguish between true, abortive promoters and non-promoters with good –10 box sequences. The methods used in this characterization of promoters can be extended to other bacteria and have important applications for promoter design in bacterial genome engineering.

Potential non-B DNA regions in the human genome are associated with higher rates of nucleotide mutation and expression variation

While individual non-B DNA structures have been shown to impact gene expression, their broad regulatory role remains elusive. We utilized genomic variants and expression quantitative trait loci (eQTL) data to analyze genome-wide variation propensities of potential non-B DNA regions and their relation to gene expression. Independent of genomic location, these regions were enriched in nucleotide variants. Our results are consistent with previously observed mutagenic properties of these regions and counter a previous study concluding that G-quadruplex regions have a reduced frequency of variants. While such mutagenicity might undermine functionality of these elements, we identified in potential non-B DNA regions a signature of negative selection. Yet, we found a depletion of eQTL-associated variants in potential non-B DNA regions, opposite to what might be expected from their proposed regulatory role. However, we also observed that genes downstream of potential non-B DNA regions showed higher expression variation between individuals. This coupling between mutagenicity and tolerance for expression variability of downstream genes may be a result of evolutionary adaptation, which allows reconciling mutagenicity of non-B DNA structures with their location in functionally important regions and their potential regulatory role.

Torsion-mediated interaction between adjacent genes

DNA torsional stress is generated by virtually all biomolecular processes involving the double helix, in particular transcription where a significant level of stress propagates over several kilobases. If another promoter is located in this range, this stress may strongly modify its opening properties, and hence facilitate or hinder its transcription. This mechanism implies that transcribed genes distant of a few kilobases are not independent, but coupled by torsional stress, an effect for which we propose the first quantitative and systematic model. In contrast to previously proposed mechanisms of transcriptional interference, the suggested coupling is not mediated by the transcription machineries, but results from the universal mechanical features of the double-helix. The model shows that the effect likely affects prokaryotes as well as eukaryotes, but with different consequences owing to their different basal levels of torsion. It also depends crucially on the relative orientation of the genes, enhancing the expression of eukaryotic divergent pairs while reducing that of prokaryotic convergent ones. To test the in vivo influence of the torsional coupling, we analyze the expression of isolated gene pairs in the Drosophila melanogaster genome. Their orientation and distance dependence is fully consistent with the model, suggesting that torsional gene coupling may constitute a widespread mechanism of (co)regulation in eukaryotes.

Role of DNA sequence based structural features of promoters in transcription initiation and gene expression

Regulatory information for transcription initiation is present in a stretch of genomic DNA, called the promoter region that is located upstream of the transcription start site (TSS) of the gene. The promoter region interacts with different transcription factors and RNA polymerase to initiate transcription and contains short stretches of transcription factor binding sites (TFBSs), as well as structurally unique elements. Recent experimental and computational analyses of promoter sequences show that they often have non-B-DNA structural motifs, as well as some conserved structural properties, such as stability, bendability, nucleosome positioning preference and curvature, across a class of organisms. Here, we briefly describe these structural features, the differences observed in various organisms and their possible role in regulation of gene expression.

Transmission of an oxygen availability signal at the Salmonella enterica serovar Typhimurium fis promoter

The nucleoid-associated protein FIS is a global regulator of gene expression and chromosome structure in Escherichia coli and Salmonella enterica. Despite the importance of FIS for infection and intracellular invasion, very little is known about the regulation of S. enterica fis expression. Under standard laboratory growth conditions, fis is highly expressed during rapid growth but is then silenced as growth slows. However, if cells are cultured in non-aerated conditions, fis expression is sustained during stationary phase. This led us to test whether the redox-sensing transcription factors ArcA and FNR regulate S. enterica fis. Deletion of FNR had no detectable effect, whereas deletion of ArcA had the unexpected effect of further elevating fis expression in stationary phase. ArcA required RpoS for induction of fis expression, suggesting that ArcA indirectly affects fis expression. Other putative regulators were found to play diverse roles: FIS acted directly as an auto-repressor (as expected), whereas CRP had little direct effect on fis expression. Deleting regions of the fis promoter led to the discovery of a novel anaerobically-induced transcription start site (Pfis-2) upstream of the primary transcription start site (Pfis-1). Promoter truncation also revealed that the shortest functional fis promoter was incapable of sustained expression. Moreover, fis expression was observed to correlate directly with DNA supercoiling in non-aerated conditions. Thus, the full-length S. enterica fis promoter region may act as a topological switch that is sensitive to stress-induced duplex destabilisation and up-regulates expression in non-aerated conditions.

Twist-DNA: computing base-pair and bubble opening probabilities in genomic superhelical DNA

SUMMARY

Local opening of the DNA double helix is required in many fundamental biological processes and is, in part, controlled by the degree of superhelicity imposed in vivo by the protein machinery. In particular, positions of superhelically destabilized regions correlate with regulatory sites along the genome. Based on a self-consistent linearization of a thermodynamic model of superhelical DNA introduced by Benham, we have developed a program that predicts the locations of these regions by efficiently computing base-pair and bubble opening probabilities in genomic DNA. The program allows visualization of results in standard genome browsers to compare DNA opening properties with other available datasets.

AVAILABILITY AND IMPLEMENTATION

Source codes freely available for download at http://www.cbp.ens-lyon.fr/doku.php?id=developpement:productions:logiciels:twistdna, implemented in Fortran90 and supported on any Unix-based operating system (Linux, Mac OS X).

The genome-wide distribution of non-B DNA motifs is shaped by operon structure and suggests the transcriptional importance of non-B DNA structures in Escherichia coli

Although the right-handed double helical B-form DNA is most common under physiological conditions, DNA is dynamic and can adopt a number of alternative structures, such as the four-stranded G-quadruplex, left-handed Z-DNA, cruciform and others. Active transcription necessitates strand separation and can induce such non-canonical forms at susceptible genomic sequences. Therefore, it has been speculated that these non-B DNA motifs can play regulatory roles in gene transcription. Such conjecture has been supported in higher eukaryotes by direct studies of several individual genes, as well as a number of large-scale analyses. However, the role of non-B DNA structures in many lower organisms, in particular proteobacteria, remains poorly understood and incompletely documented. In this study, we performed the first comprehensive study of the occurrence of B DNA-non-B DNA transition-susceptible sites (non-B DNA motifs) within the context of the operon structure of the Escherichia coli genome. We compared the distributions of non-B DNA motifs in the regulatory regions of operons with those from internal regions. We found an enrichment of some non-B DNA motifs in regulatory regions, and we show that this enrichment cannot be simply explained by base composition bias in these regions. We also showed that the distribution of several non-B DNA motifs within intergenic regions separating divergently oriented operons differs from the distribution found between convergent ones. In particular, we found a strong enrichment of cruciforms in the termination region of operons; this enrichment was observed for operons with Rho-dependent, as well as Rho-independent terminators. Finally, a preference for some non-B DNA motifs was observed near transcription factor-binding sites. Overall, the conspicuous enrichment of transition-susceptible sites in these specific regulatory regions suggests that non-B DNA structures may have roles in the transcriptional regulation of specific operons within the E. coli genome.

Theoretical analysis of competing conformational transitions in superhelical DNA

We develop a statistical mechanical model to analyze the competitive behavior of transitions to multiple alternate conformations in a negatively supercoiled DNA molecule of kilobase length and specified base sequence. Since DNA superhelicity topologically couples together the transition behaviors of all base pairs, a unified model is required to analyze all the transitions to which the DNA sequence is susceptible. Here we present a first model of this type. Our numerical approach generalizes the strategy of previously developed algorithms, which studied superhelical transitions to a single alternate conformation. We apply our multi-state model to study the competition between strand separation and B-Z transitions in superhelical DNA. We show this competition to be highly sensitive to temperature and to the imposed level of supercoiling. Comparison of our results with experimental data shows that, when the energetics appropriate to the experimental conditions are used, the competition between these two transitions is accurately captured by our algorithm. We analyze the superhelical competition between B-Z transitions and denaturation around the c-myc oncogene, where both transitions are known to occur when this gene is transcribing. We apply our model to explore the correlation between stress-induced transitions and transcriptional activity in various organisms. In higher eukaryotes we find a strong enhancement of Z-forming regions immediately 5′ to their transcription start sites (TSS), and a depletion of strand separating sites in a broad region around the TSS. The opposite patterns occur around transcript end locations. We also show that susceptibility to each type of transition is different in eukaryotes and prokaryotes. By analyzing a set of untranscribed pseudogenes we show that the Z-susceptibility just downstream of the TSS is not preserved, suggesting it may be under selection pressure.

The stresses imposed on DNA within organisms can drive the molecule from its standard B-form double-helical structure into other conformations at susceptible sites within the sequence. We present a theoretical method to calculate this transition behavior due to stresses induced by supercoiling. We also develop a numerical algorithm that calculates the transformation probability of each base pair in a user-specified DNA sequence under stress. We apply this method to analyze the competition between transitions to strand separated and left-handed Z-form structures. We find that these two conformations are both competitive under physiological environmental conditions, and that this competition is especially sensitive to temperature. By comparing its results to experimental data we also show that the algorithm properly describes the competition between melting and Z-DNA formation. Analysis of large gene sets from various organisms shows a correlation between sites of stress-induced transitions and locations that are involved in regulating gene expression.

A fundamental regulatory mechanism operating through OmpR and DNA topology controls expression of Salmonella pathogenicity islands SPI-1 and SPI-2

DNA topology has fundamental control over the ability of transcription factors to access their target DNA sites at gene promoters. However, the influence of DNA topology on protein–DNA and protein–protein interactions is poorly understood. For example, relaxation of DNA supercoiling strongly induces the well-studied pathogenicity gene ssrA (also called spiR) in Salmonella enterica, but neither the mechanism nor the proteins involved are known. We have found that relaxation of DNA supercoiling induces expression of the Salmonella pathogenicity island (SPI)-2 regulator ssrA as well as the SPI-1 regulator hilC through a mechanism that requires the two-component regulator OmpR-EnvZ. Additionally, the ompR promoter is autoregulated in the same fashion. Conversely, the SPI-1 regulator hilD is induced by DNA relaxation but is repressed by OmpR. Relaxation of DNA supercoiling caused an increase in OmpR binding to DNA and a concomitant decrease in binding by the nucleoid-associated protein FIS. The reciprocal occupancy of DNA by OmpR and FIS was not due to antagonism between these transcription factors, but was instead a more intrinsic response to altered DNA topology. Surprisingly, DNA relaxation had no detectable effect on the binding of the global repressor H-NS. These results reveal the underlying molecular mechanism that primes SPI genes for rapid induction at the onset of host invasion. Additionally, our results reveal novel features of the archetypal two-component regulator OmpR. OmpR binding to relaxed DNA appears to generate a locally supercoiled state, which may assist promoter activation by relocating supercoiling stress-induced destabilization of DNA strands. Much has been made of the mechanisms that have evolved to regulate horizontally-acquired genes such as SPIs, but parallels among the ssrA, hilC, and ompR promoters illustrate that a fundamental form of regulation based on DNA topology coordinates the expression of these genes regardless of their origins.

DNA is often considered to be a passive carrier of genetic information, but in fact DNA is an active participant in coordinating the expression of the genes it carries. This is because DNA is a dynamic molecule that can assume a wide range of topologies, and this has a direct impact on the formation of the protein–DNA complexes that drive gene expression. In a bacterium, the chromosome is supercoiled to variable levels according to environmental conditions, and supercoiling in turn governs the topology of gene promoters. Thus DNA supercoiling is able to transduce environmental signals to regulate promoter output. A previous study found that the intestinal pathogen Salmonella enterica may use changes in DNA supercoiling to detect when it has entered host immune cells, allowing the bacterium to induce the pathogenicity genes it requires to evade killing by macrophage. In dissecting the underlying molecular mechanisms, we have found that changes in DNA supercoiling also upregulate other key pathogenicity genes, and we have identified the proteins involved in this gene regulatory process. These findings indicate that a fundamental level of gene control arising from the interplay between protein transcription factors and DNA topology regulates Salmonella pathogenicity.

Bubble statistics and positioning in superhelically stressed DNA

We present a general framework to study the thermodynamic denaturation of double-stranded DNA under superhelical stress. We report calculations of position- and size-dependent opening probabilities for bubbles along the sequence. Our results are obtained from transfer-matrix solutions of the Zimm-Bragg model for unconstrained DNA and of a self-consistent linearization of the Benham model for superhelical DNA. The numerical efficiency of our method allows for the analysis of entire genomes and of random sequences of corresponding length (10(6)-10(9) base pairs). We show that, at physiological conditions, opening in superhelical DNA is strongly cooperative with average bubble sizes of 10(2)-10(3) base pairs (bp), and orders of magnitude higher than in unconstrained DNA. In heterogeneous sequences, the mean degree of base-pair opening is self-averaging, while bubble localization and statistics are dominated by sequence disorder. Compared to random sequences with identical GC-content, genomic DNA has a significantly increased probability to open large bubbles under superhelical stress. These bubbles are frequently located directly upstream of transcription start sites.

Bacterial stationary-state mutagenesis and Mammalian tumorigenesis as stress-induced cellular adaptations and the role of epigenetics

Mechanisms of cellular adaptation may have some commonalities across different organisms. Revealing these common mechanisms may provide insight in the organismal level of adaptation and suggest solutions to important problems related to the adaptation. An increased rate of mutations, referred as the mutator phenotype, and beneficial nature of these mutations are common features of the bacterial stationary-state mutagenesis and of the tumorigenic transformations in mammalian cells. We argue that these commonalities of mammalian and bacterial cells result from their stress-induced adaptation that may be described in terms of a common model. Specifically, in both organisms the mutator phenotype is activated in a subpopulation of proliferating stressed cells as a strategy to survival. This strategy is an alternative to other survival strategies, such as senescence and programmed cell death, which are also activated in the stressed cells by different subpopulations. Sustained stress-related proliferative signalling and epigenetic mechanisms play a decisive role in the choice of the mutator phenotype survival strategy in the cells. They reprogram cellular functions by epigenetic silencing of cell-cycle inhibitors, DNA repair, programmed cell death, and by activation of repetitive DNA elements. This reprogramming leads to the mutator phenotype that is implemented by error-prone cell divisions with the involvement of Y family polymerases. Studies supporting the proposed model of stress-induced cellular adaptation are discussed. Cellular mechanisms involved in the bacterial stress-induced adaptation are considered in more detail.

nocoRNAc: characterization of non-coding RNAs in prokaryotes

Background

The interest in non-coding RNAs (ncRNAs) constantly rose during the past few years because of the wide spectrum of biological processes in which they are involved. This led to the discovery of numerous ncRNA genes across many species. However, for most organisms the non-coding transcriptome still remains unexplored to a great extent. Various experimental techniques for the identification of ncRNA transcripts are available, but as these methods are costly and time-consuming, there is a need for computational methods that allow the detection of functional RNAs in complete genomes in order to suggest elements for further experiments. Several programs for the genome-wide prediction of functional RNAs have been developed but most of them predict a genomic locus with no indication whether the element is transcribed or not.

Results

We present NOCORNAc, a program for the genome-wide prediction of ncRNA transcripts in bacteria. NOCORNAc incorporates various procedures for the detection of transcriptional features which are then integrated with functional ncRNA loci to determine the transcript coordinates. We applied RNAz and NOCORNAc to the genome of Streptomyces coelicolor and detected more than 800 putative ncRNA transcripts most of them located antisense to protein-coding regions. Using a custom design microarray we profiled the expression of about 400 of these elements and found more than 300 to be transcribed, 38 of them are predicted novel ncRNA genes in intergenic regions. The expression patterns of many ncRNAs are similarly complex as those of the protein-coding genes, in particular many antisense ncRNAs show a high expression correlation with their protein-coding partner.

Conclusions

We have developed NOCORNAc, a framework that facilitates the automated characterization of functional ncRNAs. NOCORNAc increases the confidence of predicted ncRNA loci, especially if they contain transcribed ncRNAs. NOCORNAc is not restricted to intergenic regions, but it is applicable to the prediction of ncRNA transcripts in whole microbial genomes. The software as well as a user guide and example data is available at http://www.zbit.uni-tuebingen.de/pas/nocornac.htm.

Theoretical analysis of the stress induced B-Z transition in superhelical DNA

We present a method to calculate the propensities of regions within a DNA molecule to transition from B-form to Z-form under negative superhelical stresses. We use statistical mechanics to analyze the competition that occurs among all susceptible Z-forming regions at thermodynamic equilibrium in a superhelically stressed DNA of specified sequence. This method, which we call SIBZ, is similar to the SIDD algorithm that was previously developed to analyze superhelical duplex destabilization. A state of the system is determined by assigning to each base pair either the B- or the Z-conformation, accounting for the dinucleotide repeat unit of Z-DNA. The free energy of a state is comprised of the nucleation energy, the sequence-dependent B-Z transition energy, and the energy associated with the residual superhelicity remaining after the change of twist due to transition. Using this information, SIBZ calculates the equilibrium B-Z transition probability of each base pair in the sequence. This can be done at any physiologically reasonable level of negative superhelicity. We use SIBZ to analyze a variety of representative genomic DNA sequences. We show that the dominant Z-DNA forming regions in a sequence can compete in highly complex ways as the superhelicity level changes. Despite having no tunable parameters, the predictions of SIBZ agree precisely with experimental results, both for the onset of transition in plasmids containing introduced Z-forming sequences and for the locations of Z-forming regions in genomic sequences. We calculate the transition profiles of 5 kb regions taken from each of 12,841 mouse genes and centered on the transcription start site (TSS). We find a substantial increase in the frequency of Z-forming regions immediately upstream from the TSS. The approach developed here has the potential to illuminate the occurrence of Z-form regions in vivo, and the possible roles this transition may play in biological processes.

Promoter prediction in E. coli based on SIDD profiles and Artificial Neural Networks

BACKGROUND

One of the major challenges in biology is the correct identification of promoter regions. Computational methods based on motif searching have been the traditional approach taken. Recent studies have shown that DNA structural properties, such as curvature, stacking energy, and stress-induced duplex destabilization (SIDD) are useful in promoter prediction, as well. In this paper, the currently used SIDD energy threshold method is compared to the proposed artificial neural network (ANN) approach for finding promoters based on SIDD profile data.

RESULTS

When compared to the SIDD threshold prediction method, artificial neural networks showed noticeable improvements for precision, recall, and F-score over a range of values. The maximal F-score for the ANN classifier was 62.3 and 56.8 for the threshold-based classifier.

CONCLUSIONS

Artificial neural networks were used to predict promoters based on SIDD profile data. Results using this technique were an improvement over the previous SIDD threshold approach. Over a wide range of precision-recall values, artificial neural networks were more capable of identifying distinctive characteristics of promoter regions than threshold based methods.

GeneWiz browser: An Interactive Tool for Visualizing Sequenced Chromosomes

We present an interactive web application for visualizing genomic data of prokaryotic chromosomes. The tool (GeneWiz browser) allows users to carry out various analyses such as mapping alignments of homologous genes to other genomes, mapping of short sequencing reads to a reference chromosome, and calculating DNA properties such as curvature or stacking energy along the chromosome. The GeneWiz browser produces an interactive graphic that enables zooming from a global scale down to single nucleotides, without changing the size of the plot. Its ability to disproportionally zoom provides optimal readability and increased functionality compared to other browsers. The tool allows the user to select the display of various genomic features, color setting and data ranges. Custom numerical data can be added to the plot allowing, for example, visualization of gene expression and regulation data. Further, standard atlases are pre-generated for all prokaryotic genomes available in GenBank, providing a fast overview of all available genomes, including recently deposited genome sequences. The tool is available online from http://www.cbs.dtu.dk/services/gwBrowser. Supplemental material including interactive atlases is available online at http://www.cbs.dtu.dk/services/gwBrowser/suppl/.

The functional role of S/MARs in episomal vectors as defined by the stress-induced destabilization profile of the vector sequences

The scaffold/matrix attachment regions (S/MARs) are chromosomal elements that participate in the formation of chromatin domains and have origin of replication support functions. Because of all these functions, in recent years, they have been used as part of episomal vectors for gene transfer. The S/MAR of the human beta-interferon gene has been shown to support efficient episome retention and transgene expression in various mammalian cells. In Jurkat and other cells, DNA plasmid vectors containing Epstein-Barr virus origin of replication (EBV OriP) and the EBV nuclear antigen-1 gene mediate prolonged episome retention in the host cell nucleus, which, however, diminishes over time. In order to enhance retention, we combined this system with an S/MAR element. Unexpectedly, this completely eliminated the capacity of episomes to replicate. Calculation of the stress-induced DNA duplex destabilization profile of the vectors suggested that the S/MAR element had created an increase in molecular stability at the OriP site that may have disturbed replicative potential. In contrast, introduction of an alternative initiation of replication region from the beta-globin gene locus, instead of EBV OriP and the EBV nuclear antigen-1 gene, restored replicative capacity and enhanced episome retention mediated by the S/MAR. These effects were associated with a destabilization profile at the initiation of replication region. These data demonstrate a correlation between S/MAR-mediated vector retention and the presence of an unstable duplex at a replication origin, in this particular setting. We consider that the calculation of stress-induced duplex destabilization may be an informative first step in the design of units that replicate extrachromosomally, particularly as the latter present a safer and, therefore, attractive alternative to integrating viral vectors for gene therapy applications.

A stitch in time: efficient computation of genomic DNA melting bubbles

Background

It is of biological interest to make genome-wide predictions of the locations of DNA melting bubbles using statistical mechanics models. Computationally, this poses the challenge that a generic search through all combinations of bubble starts and ends is quadratic.

Results

An efficient algorithm is described, which shows that the time complexity of the task is O(NlogN) rather than quadratic. The algorithm exploits that bubble lengths may be limited, but without a prior assumption of a maximal bubble length. No approximations, such as windowing, have been introduced to reduce the time complexity. More than just finding the bubbles, the algorithm produces a stitch profile, which is a probabilistic graphical model of bubbles and helical regions. The algorithm applies a probability peak finding method based on a hierarchical analysis of the energy barriers in the Poland-Scheraga model.

Conclusion

Exact and fast computation of genomic stitch profiles is thus feasible. Sequences of several megabases have been computed, only limited by computer memory. Possible applications are the genome-wide comparisons of bubbles with promotors, TSS, viral integration sites, and other melting-related regions.

Superhelical destabilization in regulatory regions of stress response genes

Stress-induced DNA duplex destabilization (SIDD) analysis exploits the known structural and energetic properties of DNA to predict sites that are susceptible to strand separation under negative superhelical stress. When this approach was used to calculate the SIDD profile of the entire Escherichia coli K12 genome, it was found that strongly destabilized sites occur preferentially in intergenic regions that are either known or inferred to contain promoters, but rarely occur in coding regions. Here, we investigate whether the genes grouped in different functional categories have characteristic SIDD properties in their upstream flanks. We report that strong SIDD sites in the E. coli K12 genome are statistically significantly overrepresented in the upstream regions of genes encoding transcriptional regulators. In particular, the upstream regions of genes that directly respond to physiological and environmental stimuli are more destabilized than are those regions of genes that are not involved in these responses. Moreover, if a pathway is controlled by a transcriptional regulator whose gene has a destabilized 5′ flank, then the genes (operons) in that pathway also usually contain strongly destabilized SIDD sites in their 5′ flanks. We observe this statistically significant association of SIDD sites with upstream regions of genes functioning in transcription in 38 of 43 genomes of free-living bacteria, but in only four of 18 genomes of endosymbionts or obligate parasitic bacteria. These results suggest that strong SIDD sites 5′ to participating genes may be involved in transcriptional responses to environmental changes, which are known to transiently alter superhelicity. We propose that these SIDD sites are active and necessary participants in superhelically mediated regulatory mechanisms governing changes in the global pattern of gene expression in prokaryotes in response to physiological or environmental changes.

DNA in vivo experiences regulated amounts of untwisting stress. If sufficiently large, these stresses can destabilize the double helix at specific locations. These sites then become favored locations for strand separations. Gene expression and DNA replication, the two major jobs of DNA, both require the strands of the duplex to be separated. Thus, events that affect the ease of strand separation can regulate the initiation of these processes. Stress-induced DNA duplex destabilization (SIDD) has been implicated in mechanisms regulating several biological processes, including the initiation of gene expression and replication. We have developed computational methods that accurately predict the locations and extents of destabilization within genomic DNA sequences that occur in response to specified stress levels. Here, we report that the easily destabilized sites we find in the Escherichia coli K12 genome are statistically significantly overrepresented in the upstream regions of genes encoding proteins that regulate transcription. In particular, the regions upstream of genes that directly respond to physiological and environmental stimuli are more destabilized than are those regions of genes that are not involved in these responses. These results suggest that strong SIDD sites upstream of participating genes may be involved in transcriptional responses to environmental changes.

Formation of the open complex by bacterial RNA polymerase--a quantitative model

Over the last two decades, a large amount of data on initiation of transcription by bacterial RNA polymerase (RNAP) has been obtained. However, a question of how the open complex is formed still remains open, and several qualitative hypotheses for opening of DNA by RNAP have been proposed. To provide a theoretical framework needed to analyze the assembled experimental data, we here develop the first quantitative model of the open complex formation by bacterial RNAP. We first show that a simple hypothesis (which might follow from recent bioinformatic and experimental results), by which promoter DNA is melted in one step through thermal fluctuations, is inconsistent with experimental data. We next consider a more complex two-step view of the open complex formation. According to this hypothesis, the transcription bubble is formed in the -10 region, and consequently extends to the transcription start site. We derive how the open complex formation rate depends on DNA duplex melting energy and on interaction energies of RNAP with promoter DNA in the closed and open complex. This relationship provides an explicit connection between transcription initiation rate and physical properties of the promoter sequence and promoter-RNAP interactions. We compare our model with both biochemical measurements and genomics data and report a very good agreement with the experiments, with no free parameters used in model testing. This agreement therefore strongly supports both the quantitative model that we propose and the qualitative hypothesis on which the model is based. From a practical point, our results allow efficient estimation of promoter kinetic parameters, as well as engineering of promoter sequences with the desired kinetic properties.

Global features of the Alcanivorax borkumensis SK2 genome

The global feature of the completely sequenced Alcanivorax borkumensis SK2 type strain chromosome is its symmetry and homogeneity. The origin and terminus of replication are located opposite to each other in the chromosome and are discerned with high signal to noise ratios by maximal oligonucleotide usage biases on the leading and lagging strand. Genomic DNA structure is rather uniform throughout the chromosome with respect to intrinsic curvature, position preference or base stacking energy. The orthologs and paralogs of A. borkumensis genes with the highest sequence homology were found in most cases among gamma-Proteobacteria, with Acinetobacter and P. aeruginosa as closest relatives. A. borkumensis shares a similar oligonucleotide usage and promoter structure with the Pseudomonadales. A comparatively low number of only 18 genome islands with atypical oligonucleotide usage was detected in the A. borkumensis chromosome. The gene clusters that confer the assimilation of aliphatic hydrocarbons, are localized in two genome islands which were probably acquired from an ancestor of the Yersinia lineage, whereas the alk genes of Pseudomonas putida still exhibit the typical Alcanivorax oligonucleotide signature indicating a complex evolution of this major hydrocarbonoclastic trait.

Multiple co-regulatory elements and IHF are necessary for the control of fimB expression in response to sialic acid and N-acetylglucosamine in Escherichia coli K-12

Expression of the FimB recombinase, and hence the OFF-to-ON switching of type 1 fimbriation in Escherichia coli, is inhibited by sialic acid (Neu(5)Ac) and by GlcNAc. NanR (Neu(5)Ac-responsive) and NagC (GlcNAc-6P-responsive) activate fimB expression by binding to operators (O(NR) and O(NC1) respectively) located more than 600 bp upstream of the fimB promoter within the large (1.4 kb) nanC-fimB intergenic region. Here it is demonstrated that NagC binding to a second site (O(NC2)), located 212 bp closer to fimB, also controls fimB expression, and that integration host factor (IHF), which binds midway between O(NC1) and O(NC2), facilitates NagC binding to its two operator sites. In contrast, IHF does not enhance the ability of NanR to activate fimB expression in the wild-type background. Neither sequences up to 820 bp upstream of O(NR), nor those 270 bp downstream of O(NC2), are required for activation by NanR and NagC. However, placing the NanR, IHF and NagC binding sites closer to the fimB promoter enhances the ability of the regulators to activate fimB expression. These results support a refined model for how two potentially key indicators of host inflammation, Neu(5)Ac and GlcNAc, regulate type 1 fimbriation.

The recognition and prediction of σ70 promoters in Escherichia coli K-12

Based on the conservation analysis of the 683 latest experimentally verified sigma(70)-promoter sequences of Escherichia coli K-12, it is found that the conservative hexamers segments in different sites play a key role of promoter regions, a novel position-correlation scoring matrix (PCSM) algorithm for predicting sigma(70) promoter is presented. The predictive capacity of the algorithm is tested by 10-cross validation test. The results show that the overall prediction accuracies (sensitivity) and specificity are 91% and 81%, respectively. By selecting the 683 experimentally verified sigma(70) promoters as training set and searching for the complete sequence in E. coli K-12 with 4639221bp. Results show that the 100% of the 683 experimentally verified sigma(70) promoters have been identified and some possible promoters are predicted.

Efflux pump genes of the resistance-nodulation-division family in Burkholderia cenocepacia genome

Background

Burkholderia cenocepacia is recognized as opportunistic pathogen that can cause lung infections in cystic fibrosis patients. A hallmark of B. cenocepacia infections is the inability to eradicate the organism because of multiple intrinsic antibiotic resistance. As Resistance-Nodulation-Division (RND) efflux systems are responsible for much of the intrinsic multidrug resistance in Gram-negative bacteria, this study aims to identify RND genes in the B. cenocepacia genome and start to investigate their involvement into antimicrobial resistance.

Results

Genome analysis and homology searches revealed 14 open reading frames encoding putative drug efflux pumps belonging to RND family in B. cenocepacia J2315 strain.

By reverse transcription (RT)-PCR analysis, it was found that orf3, orf9, orf11, and orf13 were expressed at detectable levels, while orf10 appeared to be weakly expressed in B. cenocepacia. Futhermore, orf3 was strongly induced by chloramphenicol. The orf2 conferred resistance to fluoroquinolones, tetraphenylphosphonium, streptomycin, and ethidium bromide when cloned and expressed in Escherichia coli KAM3, a strain lacking the multidrug efflux pump AcrAB. The orf2-overexpressing E. coli also accumulate low concentrations of ethidium bromide, which was restored to wild type level in the presence of CCCP, an energy uncoupler altering the energy of the drug efflux pump.

Conclusion

The 14 RND pumps gene we have identified in the genome of B. cenocepacia suggest that active efflux could be a major mechanism underlying antimicrobial resistance in this microorganism. We have characterized the ORF2 pump, one of these 14 potential RND efflux systems. Its overexpression in E. coli conferred resistance to several antibiotics and to ethidium bromide but it remains to be determined if this pump play a significant role in the antimicrobial intrinsic resistance of B. cenocepacia. The characterization of antibiotic efflux pumps in B. cenocepacia is an obligatory step prior to the design of specific, potent bacterial inhibitors for the improved control of infectious diseases. Consequently, the topic deserves to be further investigated and future studies will involve systematic investigation on the function and expression of each of the RND efflux pump homologs.

Genome-wide prediction of G4 DNA as regulatory motifs: role in Escherichia coli global regulation

The role of nonlinear DNA in replication, recombination, and transcription has become evident in recent years. Although several studies have predicted and characterized regulatory elements at the sequence level, very few have investigated DNA structure as regulatory motifs. Here, using G-quadruplex or G4 DNA motifs as a model, we have researched the role of DNA structure in transcription on a genome-wide scale. Analyses of >61,000 open reading frames (ORFs) across 18 prokaryotes show enrichment of G4 motifs in regulatory regions and indicate its predominance within promoters of genes pertaining to transcription, secondary metabolite biosynthesis, and signal transduction. Based on this, we predict that G4 DNA may present regulatory signals. This is supported by conserved G4 motifs in promoters of orthologous genes across phylogenetically distant organisms. We hypothesized a regulatory role of G4 DNA during supercoiling stress, when duplex destabilization may result in G4 formation. This is in line with our observations from target site analysis for 55 DNA-binding proteins in Escherichia coli, which reveals significant (P<0.001) association of G4 motifs with target sites of global regulators FIS and Lrp and the sigma factor RpoD (sigma70). These factors together control >1000 genes in the early growth phase and are believed to be induced by supercoiled DNA. We also predict G4 motif-induced supercoiling sensitivity for >30 operons in E. coli, and our findings implicate G4 DNA in DNA-topology-mediated global gene regulation in E. coli.

Promoter prediction and annotation of microbial genomes based on DNA sequence and structural responses to superhelical stress

BACKGROUND

In our previous studies, we found that the sites in prokaryotic genomes which are most susceptible to duplex destabilization under the negative superhelical stresses that occur in vivo are statistically highly significantly associated with intergenic regions that are known or inferred to contain promoters. In this report we investigate how this structural property, either alone or together with other structural and sequence attributes, may be used to search prokaryotic genomes for promoters.

RESULTS

We show that the propensity for stress-induced DNA duplex destabilization (SIDD) is closely associated with specific promoter regions. The extent of destabilization in promoter-containing regions is found to be bimodally distributed. When compared with DNA curvature, deformability, thermostability or sequence motif scores within the -10 region, SIDD is found to be the most informative DNA property regarding promoter locations in the E. coli K12 genome. SIDD properties alone perform better at detecting promoter regions than other programs trained on this genome. Because this approach has a very low false positive rate, it can be used to predict with high confidence the subset of promoters that are strongly destabilized. When SIDD properties are combined with -10 motif scores in a linear classification function, they predict promoter regions with better than 80% accuracy. When these methods were tested with promoter and non-promoter sequences from Bacillus subtilis, they achieved similar or higher accuracies. We also present a strictly SIDD-based predictor for annotating promoter sequences in complete microbial genomes.

CONCLUSION

In this report we show that the propensity to undergo stress-induced duplex destabilization (SIDD) is a distinctive structural attribute of many prokaryotic promoter sequences. We have developed methods to identify promoter sequences in prokaryotic genomes that use SIDD either as a sole predictor or in combination with other DNA structural and sequence properties. Although these methods cannot predict all the promoter-containing regions in a genome, they do find large sets of potential regions that have high probabilities of being true positives. This approach could be especially valuable for annotating those genomes about which there is limited experimental data.