sRNAscanner: A Computational Tool for Intergenic Small RNA Detection in Bacterial Genomes

Background:Bacterial non-coding small RNAs (sRNAs) have attracted considerable attention due to their ubiquitous nature and contribution to numerous cellular processes including survival, adaptation and pathogenesis. Existing computational approaches for identifying bacterial sRNAs demonstrate varying levels of success and there remains considerable room for improvement. Methodology/Principal Findings: Here we have proposed a transcriptional signal-based computational method to identify intergenic sRNA transcriptional units (TUs) in completely sequenced bacterial genomes. Our sRNAscanner tool uses position weight matrices derived from experimentally defined E. coli K-12 MG1655 sRNA promoter and rho-independent terminator signals to identify intergenic sRNA TUs through sliding window based genome scans. Analysis of genomes representative of twelve species suggested that sRNAscanner demonstrated equivalent sensitivity to sRNAPredict2, the best performing bioinformatics tool available presently. However, each algorithm yielded substantial numbers of known and uncharacterized hits that were unique to one or the other tool only. sRNAscanner identified 118 novel putative intergenic sRNA genes in Salmonella enterica Typhimurium LT2, none of which were flagged by sRNAPredict2. Candidate sRNA locations were compared with available deep sequencing libraries derived from Hfq-co-immunoprecipitated RNA purified from a second Typhimurium strain (Sittka et al. (2008) PLoS Genetics 4: e1000163). Sixteen potential novel sRNAs computationally predicted and detected in deep sequencing libraries were selected for experimental validation by Northern analysis using total RNA isolated from bacteria grown under eleven different growth conditions. RNA bands of expected sizes were detected in Northern blots for six of the examined candidates. Furthermore, the 5'-ends of these six Northern-supported sRNA candidates were successfully mapped using 5'-RACE analysis. Conclusions/Significance: We have developed, computationally examined and experimentally validated the sRNAscanner algorithm. Data derived from this study has successfully identified six novel S. Typhimurium sRNA genes. In addition, the computational specificity analysis we have undertaken suggests that similar to 40% of sRNAscanner hits with high cumulative sum of scores represent genuine, undiscovered sRNA genes. Collectively, these data strongly support the utility of sRNAscanner and offer a glimpse of its potential to reveal large numbers of sRNA genes that have to date defied identification. sRNAscanner is available from: http://bicmku.in:8081/sRNAscanner or http://cluster.physics.iisc.ernet.in/sRNAscanner/.

High-quality annotation of promoter regions for 913 bacterial genomes

Motivation: The number of bacterial genomes being sequenced is increasing very rapidly and hence, it is crucial to have procedures for rapid and reliable annotation of their functional elements such as promoter regions, which control the expression of each gene or each transcription unit of the genome. The present work addresses this requirement and presents a generic method applicable across organisms. Results: Relative stability of the DNA double helical sequences has been used to discriminate promoter regions from non-promoter regions. Based on the difference in stability between neighboring regions, an algorithm has been implemented to predict promoter regions on a large scale over 913 microbial genome sequences. The average free energy values for the promoter regions as well as their downstream regions are found to differ, depending on their GC content. Threshold values to identify promoter regions have been derived using sequences flanking a subset of translation start sites from all microbial genomes and then used to predict promoters over the complete genome sequences. An average recall value of 72% (which indicates the percentage of protein and RNA coding genes with predicted promoter regions assigned to them) and precision of 56% is achieved over the 913 microbial genome dataset.

Promiscuous restriction is a cellular defense strategy that confers fitness advantage to bacteria

Most bacterial genomes harbor restriction-modification systems, encoding a REase and its cognate MTase. On attack by a foreign DNA, the REase recognizes it as nonself and subjects it to restriction. Should REases be highly specific for targeting the invading foreign DNA? It is often considered to be the case. However, when bacteria harboring a promiscuous or high-fidelity variant of the REase were challenged with bacteriophages, fitness was maximal under conditions of catalytic promiscuity. We also delineate possible mechanisms by which the REase recognizes the chromosome as self at the noncanonical sites, thereby preventing lethal dsDNA breaks. This study provides a fundamental understanding of how bacteria exploit an existing defense system to gain fitness advantage during a host-parasite coevolutionary ``arms race.''

Polyketide Quinones Are Alternate Intermediate Electron Carriers during Mycobacterial Respiration in Oxygen-Deficient Niches

Mycobacterium tuberculosis (Mtb) adaptation to hypoxia is considered crucial to its prolonged latent persistence in humans. Mtb lesions are known to contain physiologically heterogeneous microenvironments that bring about differential responses from bacteria. Here we exploit metabolic variability within biofilm cells to identify alternate respiratory polyketide quinones (PkQs) from both Mycobacterium smegmatis (Msmeg) and Mtb. PkQs are specifically expressed in biofilms and other oxygen-deficient niches to maintain cellular bioenergetics. Under such conditions, these metabolites function as mobile electron carriers in the respiratory electron transport chain. In the absence of PkQs, mycobacteria escape from the hypoxic core of biofilms and prefer oxygenrich conditions. Unlike the ubiquitous isoprenoid pathway for the biosynthesis of respiratory quinones, PkQs are produced by type III polyketide synthases using fatty acyl-CoA precursors. The biosynthetic pathway is conserved in several other bacterial genomes, and our study reveals a redox-balancing chemicocellular process in microbial physiology.

CsrA interacting small RNAs in Haemophilus spp genomes: a theoretical analysis

The csrA is a carbon storage regulator gene that encodes a protein with multiple RNA interaction sites. Bacterial non-coding small RNAs like csrB, csrC and their counterparts in diverse bacterial genus are identified to control the regulatory activities of CsrA and its orthologs. An attempt has been made in this study to identify 'novel' non-coding small RNAs that are involved in the regulatory activities of csrA gene. All CsrA-interacting small RNAs are computationally fingerprinted to have multiple occurrence of 7-nucleotide CsrA interacting repeats [CAGGA(U/A/C)G] along with a 18-nucleotide upstream binding site. However, in several of the genomes like Haemophilus spp, the upstream binding site is not identified. The current methodology overcomes this difficulty by identifying small RNA-specific orphan transcriptional units within the intergenic regions of the genome. The results could identify all known CsrA-interacting small RNAs in E. coli, Vibrio cholerae and Pseudomonas aeruginosa genomes and additionally has picked six new possible CsrA-interacting small RNA regions in E. coli. Our computational analysis indicates that known rygD and rprA sRNAs in E. coli could possibly interact with CsrA proteins. The study is extended to three of the Haemophilus genomes that could identify seven new possible CsrA interacting small RNAs.

Distribution of simple repetitive (TG/CA)n and (CT/AG)n sequences in human and rodent genomes

Sixteen million nucleotide sequence of genome of various organisms have been analysed to detect and study the extent of occurrence of simple repetitive sequences. Two sequence motifs (TG/CA)n and (CT/AG)n capable of adopting unusual DNA structures, left handed Z-conformation and triple-helical conformation respectively, are found to be abundant in rodent and human genomes, but almost completely absent in bacterial genome. (TG/CA)n and (CT/AG)n sequences are present mostly in the intron or 5'/3' flanking regions of the genes. The presence of such repeat motifs in genomic sequence of higher eukaryotes has been correlated with their possible functional significance in nucleosome organization, recombination and gene expression.

Functional correlation of bacterial LuxS with their quaternary associations: interface analysis of the structure networks

Background The genome of a wide variety of prokaryotes contains the luxS gene homologue, which encodes for the protein S-ribosylhomocysteinelyase (LuxS). This protein is responsible for the production of the quorum sensing molecule, AI-2 and has been implicated in a variety of functions such as flagellar motility, metabolic regulation, toxin production and even in pathogenicity. A high structural similarity is present in the LuxS structures determined from a few species. In this study, we have modelled the structures from several other species and have investigated their dimer interfaces. We have attempted to correlate the interface features of LuxS with the phenotypic nature of the organisms. Results The protein structure networks (PSN) are constructed and graph theoretical analysis is performed on the structures obtained from X-ray crystallography and on the modelled ones. The interfaces, which are known to contain the active site, are characterized from the PSNs of these homodimeric proteins. The key features presented by the protein interfaces are investigated for the classification of the proteins in relation to their function. From our analysis, structural interface motifs are identified for each class in our dataset, which showed distinctly different pattern at the interface of LuxS for the probiotics and some extremophiles. Our analysis also reveals potential sites of mutation and geometric patterns at the interface that was not evident from conventional sequence alignment studies. Conclusion The structure network approach employed in this study for the analysis of dimeric interfaces in LuxS has brought out certain structural details at the side-chain interaction level, which were elusive from the conventional structure comparison methods. The results from this study provide a better understanding of the relation between the luxS gene and its functional role in the prokaryotes. This study also makes it possible to explore the potential direction towards the design of inhibitors of LuxS and thus towards a wide range of antimicrobials.

Purification of a bacterial organophosphate-hydrolysing phosphatase by cibacron 3GA-Sepharose affinity chromatography

A phosphatase catalysing the hydrolysis of organophosphorus pesticides was purified to homogeneity using Cibacron 3GA-Sepharose CL 6B affinity chromatography. The enzyme which is localized in the periplasm of the bacterium Image NC5 was extracted by treating with 0.2M MgCl2, pH 8.4. The enzyme was adsorbed to the Cibacron-Sepharose at pH 7.0 and eluted with Tris-HCl buffer at pH 8.0, with 47 per cent recovery. The enzyme thus obtained was electrophoretically homogeneous. This simple affinity purification procedure enhances the potential for its use in large scale detoxification systems.

Future switching technology: bacterial protein-based programmable all-optical switch for integrated optics applications?

We report on the bacterial protein-based all-optical switches which operate at low laser power, high speed and fulfil most of the requirements to be an ideal all-optical switch without any moving parts involved. This consists of conventional optical waveguides coated with bacteriorhodopsin films at switching locations. The principle of operation of the switch is based on the light-induced refractive index change of bacteriorhodopsin. This approach opens the possibility of realizing proteinbased all-optical switches for communication network, integrated optics and optical computers.

The BglG group of antiterminators: a growing family of bacterial regulators

The product of the bglG gene of Escherichia coli was among the first bacterial antiterminators to be identified and characterized. Since the elucidation ten years ago of its role in the regulation of the bgl operon of E. coli,a large number of homologies have been discovered in both Gram-positive and Gram-negative bacteria. Often the homologues of BglG in other organisms are also involved in regulating β-glucoside utilization. Surprisingly, in many cases, they mediate antitermination to regulate a variety of other catabolic functions. Because of the high degree of conservation of the cis-acting regulatory elements, antiterminators from one organism can function in another. Generally the antiterminator protein itself is negatively regulated by phosphorylation by a component of the phosphotransferase system. This family of proteins thus represents a highly evolved regulatory system that is conserved across evolutionarily distant genuses.

Relative stability of DNA as a generic criterion for promoter prediction: whole genome annotation of microbial genomes with varying nucleotide base composition

The rapid increase in genome sequence information has necessitated the annotation of their functional elements, particularly those occurring in the non-coding regions, in the genomic context. Promoter region is the key regulatory region, which enables the gene to be transcribed or repressed, but it is difficult to determine experimentally. Hence an in silico identification of promoters is crucial in order to guide experimental work and to pin point the key region that controls the transcription initiation of a gene. In this analysis, we demonstrate that while the promoter regions are in general less stable than the flanking regions, their average free energy varies depending on the GC composition of the flanking genomic sequence. We have therefore obtained a set of free energy threshold values, for genomic DNA with varying GC content and used them as generic criteria for predicting promoter regions in several microbial genomes, using an in-house developed tool `PromPredict'. On applying it to predict promoter regions corresponding to the 1144 and 612 experimentally validated TSSs in E. coli (50.8% GC) and B. subtilis (43.5% GC) sensitivity of 99% and 95% and precision values of 58% and 60%, respectively, were achieved. For the limited data set of 81 TSSs available for M. tuberculosis (65.6% GC) a sensitivity of 100% and precision of 49% was obtained.

A Mycobacterial Cyclic AMP Phosphodiesterase That Moonlights as a Modifier of Cell Wall Permeability

Mycobacterium tuberculosis utilizes many mechanisms to establish itself within the macrophage, and bacterially derived cAMP is important in modulating the host cellular response. Although the genome of M. tuberculosis is endowed with a number of mammalian-like adenylyl cyclases, only a single cAMP phosphodiesterase has been identified that can decrease levels of cAMP produced by the bacterium. We present the crystal structure of the full-length and sole cAMP phosphodiesterase, Rv0805, found in M. tuberculosis, whose orthologs are present only in /the genomes of slow growing and pathogenic mycobacteria. The dimeric core catalytic domain of Rv0805 adopts a metallophosphoesterase fold, and the C-terminal region builds the active site and contributes to multiple substrate utilization.Localization of Rv0805 to the cell wall is dependent on its C terminus, and expression of either wild type or mutationally inactivated Rv0805 in M. smegmatis alters cell permeability to hydrophobic cytotoxic compounds. Rv0805 may therefore play a key role in the pathogenicity of mycobacteria, not only by hydrolyzing bacterial cAMP, but also by moonlighting as a protein that can alter cell wall functioning.

Inducible nitric oxide synthase and control of intracellular bacterial pathogens

Inducible nitric oxide synthase (iNOS) has important functions in innate immunity and regulation of immune functions. Here, the role of iNOS in the pathogenesis of various intracellular bacterial infections is discussed. These pathogens have also evolved a broad array of strategies to repair damage by reactive nitrogen intermediates, and to suppress or inhibit functions of iNOS.

A Single Mammalian Mitochondrial Translation Initiation Factor Functionally Replaces Two Bacterial Factors

The mechanism of translation in eubacteria and organelles is thought to be similar. In eubacteria, the three initiation factors IF1, IF2, and IF3 are vital. Although the homologs of IF2 and IF3 are found in mammalian mitochondria, an IF1 homolog has never been detected. Here, we show that bovine mitochondrial IF2 (IF2mt) complements E. coli containing a deletion of the IF2 gene (E. coli ΔinfB). We find that IF1 is no longer essential in an IF2mt-supported E. coli ΔinfB strain. Furthermore, biochemical and molecular modeling data show that a conserved insertion of 37 amino acids in the IF2mt substitutes for the function of IF1. Deletion of this insertion from IF2mt supports E. coli for the essential function of IF2. However, in this background, IF1 remains essential. These observations provide strong evidence that a single factor (IF2mt) in mammalian mitochondria performs the functions of two eubacterial factors, IF1 and IF2.