Examination of the structural and functional versatility of glmS ribozymes by using in vitro selection

Self-cleaving ribozymes associated with the glmS genes of many Gram-positive bacteria are activated by binding to glucosamine-6-phosphate (GlcN6P). Representatives of the glmS ribozyme class function as metabolite-sensing riboswitches whose self-cleavage activities down-regulate the expression of GlmS enzymes that synthesizes GlcN6P. As with other riboswitches, natural glmS ribozyme isolates are highly specific for their target metabolite. Other small molecules closely related to GlcN6P, such as glucose-6-phosphate, cannot activate self-cleavage. We applied in vitro selection methods in an attempt to identify variants of a Bacillus cereus glmS ribozyme that expand the range of compounds that induce self-cleavage. In addition, we sought to increase the number of variant ribozymes of this class to further examine the proposed secondary structure model. Although numerous variant ribozymes were obtained that efficiently self-cleave, none exhibited changes in target specificity. These findings are consistent with the hypothesis that GlcN6P is used by the ribozyme as a coenzyme for RNA cleavage, rather than an allosteric effector

Folding of the cocaine aptamer studied by EPR and fluorescence spectroscopies using the bifunctional spectroscopic probe Ç

The cocaine aptamer is a DNA molecule that binds cocaine at the junction of three helices. The bifunctional spectroscopic probe Ç was incorporated independently into three different positions of the aptamer and changes in structure and dynamics upon addition of the cocaine ligand were studied. Nucleoside Ç contains a rigid nitroxide spin label and can be studied directly by electron paramagnetic resonance (EPR) spectroscopy and fluorescence spectroscopy after reduction of the nitroxide to yield the fluoroside Çf. Both the EPR and the fluorescence data for aptamer 2 indicate that helix III is formed before cocaine binding. Upon addition of cocaine, increased fluorescence of a fully base-paired Çf, placed at the three-way junction in helix III, was observed and is consistent with a helical tilt from a coaxial stack of helices II and III. EPR and fluorescence data clearly show that helix I is formed upon addition of cocaine, concomitant with the formation of the Y-shaped three-way helical junction. The EPR data indicate that nucleotides in helix I are more mobile than nucleotides in regular duplex regions and may reflect increased dynamics due to the short length of helix I

Comparative genomics of metabolic capacities of regulons controlled by cis-regulatory RNA motifs in bacteria

BACKGROUND: In silico comparative genomics approaches have been efficiently used for functional prediction and reconstruction of metabolic and regulatory networks. Riboswitches are metabolite-sensing structures often found in bacterial mRNA leaders controlling gene expression on transcriptional or translational levels. An increasing number of riboswitches and other cis-regulatory RNAs have been recently classified into numerous RNA families in the Rfam database. High conservation of these RNA motifs provides a unique advantage for their genomic identification and comparative analysis. RESULTS: A comparative genomics approach implemented in the RegPredict tool was used for reconstruction and functional annotation of regulons controlled by RNAs from 43 Rfam families in diverse taxonomic groups of Bacteria. The inferred regulons include ~5200 cis-regulatory RNAs and more than 12000 target genes in 255 microbial genomes. All predicted RNA-regulated genes were classified into specific and overall functional categories. Analysis of taxonomic distribution of these categories allowed us to establish major functional preferences for each analyzed cis-regulatory RNA motif family. Overall, most RNA motif regulons showed predictable functional content in accordance with their experimentally established effector ligands. Our results suggest that some RNA motifs (including thiamin pyrophosphate and cobalamin riboswitches that control the cofactor metabolism) are widespread and likely originated from the last common ancestor of all bacteria. However, many more analyzed RNA motifs are restricted to a narrow taxonomic group of bacteria and likely represent more recent evolutionary innovations. CONCLUSIONS: The reconstructed regulatory networks for major known RNA motifs substantially expand the existing knowledge of transcriptional regulation in bacteria. The inferred regulons can be used for genetic experiments, functional annotations of genes, metabolic reconstruction and evolutionary analysis. The obtained genome-wide collection of reference RNA motif regulons is available in the RegPrecise database (http://regprecise.lbl.gov/)

Discovering cis-Regulatory RNAs in Shewanella Genomes by Support Vector Machines

An increasing number of cis-regulatory RNA elements have been found to regulate gene expression post-transcriptionally in various biological processes in bacterial systems. Effective computational tools for large-scale identification of novel regulatory RNAs are strongly desired to facilitate our exploration of gene regulation mechanisms and regulatory networks. We present a new computational program named RSSVM (RNA Sampler+Support Vector Machine), which employs Support Vector Machines (SVMs) for efficient identification of functional RNA motifs from random RNA secondary structures. RSSVM uses a set of distinctive features to represent the common RNA secondary structure and structural alignment predicted by RNA Sampler, a tool for accurate common RNA secondary structure prediction, and is trained with functional RNAs from a variety of bacterial RNA motif/gene families covering a wide range of sequence identities. When tested on a large number of known and random RNA motifs, RSSVM shows a significantly higher sensitivity than other leading RNA identification programs while maintaining the same false positive rate. RSSVM performs particularly well on sets with low sequence identities. The combination of RNA Sampler and RSSVM provides a new, fast, and efficient pipeline for large-scale discovery of regulatory RNA motifs. We applied RSSVM to multiple Shewanella genomes and identified putative regulatory RNA motifs in the 5′ untranslated regions (UTRs) in S. oneidensis, an important bacterial organism with extraordinary respiratory and metal reducing abilities and great potential for bioremediation and alternative energy generation. From 1002 sets of 5′-UTRs of orthologous operons, we identified 166 putative regulatory RNA motifs, including 17 of the 19 known RNA motifs from Rfam, an additional 21 RNA motifs that are supported by literature evidence, 72 RNA motifs overlapping predicted transcription terminators or attenuators, and other candidate regulatory RNA motifs. Our study provides a list of promising novel regulatory RNA motifs potentially involved in post-transcriptional gene regulation. Combined with the previous cis-regulatory DNA motif study in S. oneidensis, this genome-wide discovery of cis-regulatory RNA motifs may offer more comprehensive views of gene regulation at a different level in this organism. The RSSVM software, predictions, and analysis results on Shewanella genomes are available at http://ural.wustl.edu/resources.html#RSSVM

Folding of the lysine riboswitch: importance of peripheral elements for transcriptional regulation

The Bacillus subtilis lysC lysine riboswitch modulates its own gene expression upon lysine binding through a transcription attenuation mechanism. The riboswitch aptamer is organized around a single five-way junction that provides the scaffold for two long-range tertiary interactions (loop L2–loop L3 and helix P2–loop L4)—all of this for the creation of a specific lysine binding site. We have determined that the interaction P2–L4 is particularly important for the organization of the ligand-binding site and for the riboswitch transcription attenuation control. Moreover, we have observed that a folding synergy between L2–L3 and P2–L4 allows both interactions to fold at lower magnesium ion concentrations. The P2–L4 interaction is also critical for the close juxtaposition involving stems P1 and P5. This is facilitated by the presence of lysine, suggesting an active role of the ligand in the folding transition. We also show that a previously uncharacterized stem–loop located in the expression platform is highly important for the riboswitch activity. Thus, folding elements located in the aptamer and the expression platform both influence the lysine riboswitch gene regulation

Comparative genomics and evolution of bacterial regulatory systems

Abstract: Recent advances in genome sequencing and development of comparative genomics techniques allow one to study evolution of regulation in prokaryotes at several different levels: microevolution of orthologous regulatory sites, changes in regulon content, evolution of interacting regulatory systems, and co-evolution of transcription factors and their binding signals. Regulatory interactions appear to be very dynamic in some cases and surprisingly stable in others. The review presents several examples where comparative analysis uncovered plausible scenarios of evolution of regulatory systems

Two groups of phenylalanine biosynthetic operon leader peptides genes: a high level of apparently incidental frameshifting in decoding Escherichia coli pheL

The bacterial pheL gene encodes the leader peptide for the phenylalanine biosynthetic operon. Translation of pheL mRNA controls transcription attenuation and, consequently, expression of the downstream pheA gene. Fifty-three unique pheL genes have been identified in sequenced genomes of the gamma subdivision. There are two groups of pheL genes, both of which are short and contain a run(s) of phenylalanine codons at an internal position. One group is somewhat diverse and features different termination and 5′-flanking codons. The other group, mostly restricted to Enterobacteria and including Escherichia coli pheL, has a conserved nucleotide sequence that ends with UUC_CCC_UGA. When these three codons in E. coli pheL mRNA are in the ribosomal E-, P- and A-sites, there is an unusually high level, 15%, of +1 ribosomal frameshifting due to features of the nascent peptide sequence that include the penultimate phenylalanine. This level increases to 60% with a natural, heterologous, nascent peptide stimulator. Nevertheless, studies with different tRNAPro mutants in Salmonella enterica suggest that frameshifting at the end of pheL does not influence expression of the downstream pheA. This finding of incidental, rather than utilized, frameshifting is cautionary for other studies of programmed frameshifting