Riboswitches and the role of noncoding RNAs in bacterial metabolic control

Enzymes in riboflavin biosynthesis: Potential antibiotic drug targets

The rapid resistance of pathogens to antibiotics has emerged as a major threat to global health. Identification of new antibiotic targets is thus needed for developing alternative drugs. Genes encoding enzymes involved in the biosynthesis of riboflavin and flavin cofactors (FMN/FAD) are attractive targets because these enzymatic reactions are necessary for most bacteria to synthesize flavin cofactors for use in their central metabolic reactions. Moreover, humans lack most of these enzymes because we uptake riboflavin from our diet. This review discusses the current knowledge of enzymes involved in bacterial biosynthesis of riboflavin and other flavin cofactors, as well as the functions of the FMN riboswitch. Here, we highlight recent progress in the structural and mechanistic characterization, and inhibition of GTP cyclohydrolase II (GCH II), lumazine synthase (LS), riboflavin synthase (RFS), FAD synthetase (FADS), and FMN riboswitch, which have been identified as plausible antibiotic targets. As the structures and functions of these enzymes and regulatory systems are not completely understood, they are attractive as subjects for future in-depth biochemical and biophysical analysis.

RLBind: a deep learning method to predict RNA-ligand binding sites

Identification of RNA-small molecule binding sites plays an essential role in RNA-targeted drug discovery and development. These small molecules are expected to be leading compounds to guide the development of new types of RNA-targeted therapeutics compared with regular therapeutics targeting proteins. RNAs can provide many potential drug targets with diverse structures and functions. However, up to now, only a few methods have been proposed. Predicting RNA-small molecule binding sites still remains a big challenge. New computational model is required to better extract the features and predict RNA-small molecule binding sites more accurately. In this paper, a deep learning model, RLBind, was proposed to predict RNA-small molecule binding sites from sequence-dependent and structure-dependent properties by combining global RNA sequence channel and local neighbor nucleotides channel. To our best knowledge, this research was the first to develop a convolutional neural network for RNA-small molecule binding sites prediction. Furthermore, RLBind also can be used as a potential tool when the RNA experimental tertiary structure is not available. The experimental results show that RLBind outperforms other state-of-the-art methods in predicting binding sites. Therefore, our study demonstrates that the combination of global information for full-length sequences and local information for limited local neighbor nucleotides in RNAs can improve the model's predictive performance for binding sites prediction. All datasets and resource codes are available at https://github.com/KailiWang1/RLBind.

Bottom-Up Synthetic Biology Using Cell-Free Protein Synthesis

Technical advances in biotechnology have greatly accelerated the development of bottom-up synthetic biology. Unlike top-down approaches, bottom-up synthetic biology focuses on the construction of a minimal cell from scratch and the application of these principles to solve challenges. Cell-free protein synthesis (CFPS) systems provide minimal machinery for transcription and translation, from either a fractionated cell lysate or individual purified protein elements, thus speeding up the development of synthetic cell projects. In this review, we trace the history of the cell-free technique back to the first in vitro fermentation experiment using yeast cell lysate. Furthermore, we summarized progresses of individual cell mimicry modules, such as compartmentalization, gene expression regulation, energy regeneration and metabolism, growth and division, communication, and motility. Finally, current challenges and future perspectives on the field are outlined.

TBDB: a database of structurally annotated T-box riboswitch:tRNA pairs

T-box riboswitches constitute a large family of tRNA-binding leader sequences that play a central role in gene regulation in many gram-positive bacteria. Accurate inference of the tRNA binding to T-box riboswitches is critical to predict their cis-regulatory activity. However, there is no central repository of information on the tRNA binding specificities of T-box riboswitches, and de novo prediction of binding specificities requires advanced knowledge of computational tools to annotate riboswitch secondary structure features. Here, we present the T-box Riboswitch Annotation Database (TBDB, https://tbdb.io), an open-access database with a collection of 23,535 T-box riboswitch sequences, spanning the major phyla of 3,632 bacterial species. Among structural predictions, the TBDB also identifies specifier sequences, cognate tRNA binding partners, and downstream regulatory targets. To our knowledge, the TBDB presents the largest collection of feature, sequence, and structural annotations carried out on this important family of regulatory RNA.

ABC-F translation factors: from antibiotic resistance to immune response

Energy-dependent translational throttle A (EttA) from Escherichia coli is a paradigmatic ABC-F protein that controls the first step in polypeptide elongation on the ribosome according to the cellular energy status. Biochemical and structural studies have established that ABC-F proteins generally function as translation factors that modulate the conformation of the peptidyl transferase center upon binding to the ribosomal tRNA exit site. These factors, present in both prokaryotes and eukaryotes but not in archaea, use related molecular mechanisms to modulate protein synthesis for heterogenous purposes, ranging from antibiotic resistance and rescue of stalled ribosomes to modulation of the mammalian immune response. Here, we review the canonical studies characterizing the phylogeny, regulation, ribosome interactions, and mechanisms of action of the bacterial ABC-F proteins, and discuss the implications of these studies for the molecular function of eukaryotic ABC-F proteins, including the three human family members.

Cell-Free Protein Synthesis: Chassis toward the Minimal Cell

The quest for a minimal cell not only sheds light on the fundamental principles of life but also brings great advances in related applied fields such as general biotechnology. Minimal cell projects came from the study of a plausible route to the origin of life. Later on, research extended and also referred to the construction of artificial cells, or even more broadly, as in vitro synthetic biology. The cell-free protein synthesis (CFPS) techniques harness the central cellular activity of transcription/translation in an open environment, providing the framework for multiple cellular processes assembling. Therefore, CFPS systems have become the first choice in the construction of the minimal cell. In this review, we focus on the recent advances in the quantitative analysis of CFPS and on its advantage for addressing the bottom-up assembly of a minimal cell and illustrate the importance of systemic chassis behavior, such as stochasticity under a compartmentalized micro-environment.

Bacterial Small RNAs in the Genus Herbaspirillum spp

The genus Herbaspirillum includes several strains isolated from different grasses. The identification of non-coding RNAs (ncRNAs) in the genus Herbaspirillum is an important stage studying the interaction of these molecules and the way they modulate physiological responses of different mechanisms, through RNA⁻RNA interaction or RNA⁻protein interaction. This interaction with their target occurs through the perfect pairing of short sequences (cis-encoded ncRNAs) or by the partial pairing of short sequences (trans-encoded ncRNAs). However, the companion Hfq can stabilize interactions in the trans-acting class. In addition, there are Riboswitches, located at the 5' end of mRNA and less often at the 3' end, which respond to environmental signals, high temperatures, or small binder molecules. Recently, CRISPR (clustered regularly interspaced palindromic repeats), in prokaryotes, have been described that consist of serial repeats of base sequences (spacer DNA) resulting from a previous exposure to exogenous plasmids or bacteriophages. We identified 285 ncRNAs in Herbaspirillum seropedicae (H. seropedicae) SmR1, expressed in different experimental conditions of RNA-seq material, classified as cis-encoded ncRNAs or trans-encoded ncRNAs and detected RNA riboswitch domains and CRISPR sequences. The results provide a better understanding of the participation of this type of RNA in the regulation of the metabolism of bacteria of the genus Herbaspirillum spp.

Computational Methods for Modeling Aptamers and Designing Riboswitches

Riboswitches, which are located within certain noncoding RNA region perform functions as genetic "switches", regulating when and where genes are expressed in response to certain ligands. Understanding the numerous functions of riboswitches requires computation models to predict structures and structural changes of the aptamer domains. Although aptamers often form a complex structure, computational approaches, such as RNAComposer and Rosetta, have already been applied to model the tertiary (three-dimensional (3D)) structure for several aptamers. As structural changes in aptamers must be achieved within the certain time window for effective regulation, kinetics is another key point for understanding aptamer function in riboswitch-mediated gene regulation. The coarse-grained self-organized polymer (SOP) model using Langevin dynamics simulation has been successfully developed to investigate folding kinetics of aptamers, while their co-transcriptional folding kinetics can be modeled by the helix-based computational method and BarMap approach. Based on the known aptamers, the web server Riboswitch Calculator and other theoretical methods provide a new tool to design synthetic riboswitches. This review will represent an overview of these computational methods for modeling structure and kinetics of riboswitch aptamers and for designing riboswitches.

Reversible-Switch Mechanism of the SAM-III Riboswitch

Riboswitches are self-regulatory elements located at the 5' untranslated region of certain mRNAs. The Enterococcus faecalis SAM-III (S_MK) riboswitch regulates downstream gene expression through conformational change by sensing S-adenosylmethionine (SAM) at the translation level. Using the recently developed systematic helix-based computational method, we studied the co-transcriptional folding behavior of the S_MK riboswitch and its shortened construct lacking the first six nucleotides. We find that there are no obvious misfolded structures formed during the transcription and refolding processes for this riboswitch. The full-length riboswitch quickly folds into the ON-state in the absence of SAM, and the coupling between transcription and translation is not required for the riboswitch to function. The potential to form helix P0 is necessary for the riboswitch to function as a switch. For this thermodynamically controlled reversible riboswitch, the fast helix-exchanging transition pathway between the two functional structures guaranteed that this riboswitch can act as a reversible riboswitch.

Riboswitches: From living biosensors to novel targets of antibiotics

Riboswitches are generally located in 5'-UTR region of mRNAs and specifically bind small ligands. Following ligand binding, gene expression is controlled mostly by transcription termination, translation inhibition or mRNA degradation processes. More than 30 classes of known riboswitches have been identified by now. Most riboswitches consist of an aptamer domain and an expression platform. The aptamer domain of each class of riboswitch is a conserved structure and stabilizes specific structures of the expression platforms through binding to specific compounds. In this review, we are highlighting most aspects of riboswitch research including the novel riboswitch discoveries, routine methods for discovering and investigating riboswitches along with newly discovered classes and mechanistic principles of riboswitch-mediated gene expression control. Moreover, we will give an overview about the potential of riboswitches as therapeutic targets for antibiotic design and also their utilization as biosensors for molecular detection.

Genome-Wide Detection of Predicted Non-coding RNAs Related to the Adhesion Process in Vibrio alginolyticus Using High-Throughput Sequencing

The ability of bacteria to adhere to fish mucus can be affected by environmental conditions and is considered to be a key virulence factor of Vibrio alginolyticus. However, the molecular mechanism underlying this ability remains unclear. Our previous study showed that stress conditions such as exposure to Cu, Pb, Hg, and low pH are capable of reducing the adhesion ability of V. alginolyticus. Non-coding RNAs (ncRNAs) play a crucial role in the intricate regulation of bacterial gene expression, thereby affecting bacterial pathogenicity. Thus, we hypothesized that ncRNAs play a key role in the V. alginolyticus adhesion process. To validate this, we combined high-throughput sequencing with computational techniques to detect ncRNA dynamics in samples after stress treatments. The expression of randomly selected novel ncRNAs was confirmed by QPCR. Among the significantly altered ncRNAs, 30 were up-regulated and 2 down-regulated by all stress treatments. The QPCR results reinforced the reliability of the sequencing data. Target prediction and KEGG pathway analysis indicated that these ncRNAs are closely related to pathways associated with in vitro adhesion, and our results indicated that chemical stress-induced reductions in the adhesion ability of V. alginolyticus might be due to the perturbation of ncRNA expression. Our findings provide important information for further functional characterization of ncRNAs during the adhesion process of V. alginolyticus.

Chemical feasibility of the general acid/base mechanism of glmS ribozyme self-cleavage

In numerous Gram-positive bacteria, the glmS ribozyme or catalytic riboswitch regulates the expression of glucosamine-6-phosphate (GlcN6P) synthase via site-specific cleavage of its sugar-phosphate backbone in response to GlcN6P ligand binding. Biochemical data have suggested a crucial catalytic role for an active site guanine (G40 in Thermoanaerobacter tengcongensis, G33 in Bacillus anthracis). We used hybrid quantum chemical/molecular mechanical (QM/MM) calculations to probe the mechanism where G40 is deprotonated and acts as a general base. The calculations suggest that the deprotonated guanine G40(-) is sufficiently reactive to overcome the thermodynamic penalty arising from its rare protonation state, and thus is able to activate the A-1(2'-OH) group toward nucleophilic attack on the adjacent backbone. Furthermore, deprotonation of A-1(2'-OH) and nucleophilic attack are predicted to occur as separate steps, where activation of A-1(2'-OH) precedes nucleophilic attack. Conversely, the transition state associated with the rate-determining step corresponds to concurrent nucleophilic attack and protonation of the G1(O5') leaving group by the ammonium moiety of the GlcN6P cofactor. Overall, our calculations help to explain the crucial roles of G40 (as a general base) and GlcN6P (as a general acid) during glmS ribozyme self-cleavage. In addition, we show that the QM/MM description of the glmS ribozyme self-cleavage reaction is significantly more sensitive to the size of the QM region and the quality of the QM-MM coupling than that of other small ribozymes.

The rise of regulatory RNA

Discoveries over the past decade portend a paradigm shift in molecular biology. Evidence suggests that RNA is not only functional as a messenger between DNA and protein but also involved in the regulation of genome organization and gene expression, which is increasingly elaborate in complex organisms. Regulatory RNA seems to operate at many levels; in particular, it plays an important part in the epigenetic processes that control differentiation and development. These discoveries suggest a central role for RNA in human evolution and ontogeny. Here, we review the emergence of the previously unsuspected world of regulatory RNA from a historical perspective.

Structural insights into the interactions of xpt riboswitch with novel guanine analogues: a molecular dynamics simulation study

Ligand recognition in purine riboswitches is a complex process requiring different levels of conformational changes. Recent efforts in the area of purine riboswitch research have focused on ligand analogue binding studies. In the case of the guanine xanthine phosphoribosyl transferase (xpt) riboswitch, synthetic analogues that resemble guanine have the potential to tightly bind and subsequently influence the genetic expression of xpt mRNA in prokaryotes. We have carried out 25 ns Molecular Dynamics (MD) simulation studies of the aptamer domain of the xpt G-riboswitch in four different states: guanine riboswitch in free form, riboswitch bound with its cognate ligand guanine, and with two guanine analogues SJ1 and SJ2. Our work reveals novel interactions of SJ1 and SJ2 ligands with the binding core residues of the riboswitch. The ligands proposed in this work bind to the riboswitch with greater overall stability and lower root mean square deviations and fluctuations compared to guanine ligand. Reporter gene assay data demonstrate that the ligand analogues, upon binding to the RNA, lower the genetic expression of the guanine riboswitch. Our work has important implications for future ligand design and binding studies in the exciting field of riboswitches.

The promise of riboswitches as potential antibacterial drug targets

Riboswitches represent promising novel RNA structures for developing compounds that artificially regulate gene expression and, thus, bacterial growth. The past years have seen increasing efforts to identify metabolite-analogues which act on riboswitches and which reveal antibacterial activity. Here, we summarize the current inventory of riboswitch-targeting compounds, their characteristics and antibacterial potential.

Computational identification of functional RNA homologs in metagenomic data

A key step toward understanding a metagenomics data set is the identification of functional sequence elements within it, such as protein coding genes and structural RNAs. Relative to protein coding genes, structural RNAs are more difficult to identify because of their reduced alphabet size, lack of open reading frames, and short length. Infernal is a software package that implements “covariance models” (CMs) for RNA homology search, which harness both sequence and structural conservation when searching for RNA homologs. Thanks to the added statistical signal inherent in the secondary structure conservation of many RNA families, Infernal is more powerful than sequence-only based methods such as BLAST and profile HMMs. Together with the Rfam database of CMs, Infernal is a useful tool for identifying RNAs in metagenomics data sets.

Transient RNA structure features are evolutionarily conserved and can be computationally predicted

Functional RNA structures tend to be conserved during evolution. This finding is, for example, exploited by comparative methods for RNA secondary structure prediction that currently provide the state-of-art in terms of prediction accuracy. We here provide strong evidence that homologous RNA genes not only fold into similar final RNA structures, but that their folding pathways also share common transient structural features that have been evolutionarily conserved. For this, we compile and investigate a non-redundant data set of 32 sequences with known transient and final RNA secondary structures and devise a dedicated computational analysis pipeline.

Use of RNA in drug design

This is a review of RNA as a target for small molecules (ribosomes, riboswitches, regulatory RNAs) and RNA-derived oligonucleotides as tools (antisense/small interfering RNA, ribozymes, aptamers/decoy RNA and microRNA). This review highlights the present state of research using RNA as a drug target or as a potential drug candidate and explains at which stage and to what extent rational design could eventually be involved. Special attention has been paid to the recent potential clinical applications of RNA either as drugs or drug targets. The review deals mainly with mechanistic approaches rather than with physicochemical or computational aspects of RNA-based drug design.

Fast alignment and comparison of RNA structures

MOTIVATION

To recognize remote relationships between RNA molecules, one must be able to align structures without regard to sequence similarity. We have implemented a method, which is swift [O(n(2))], sensitive and tolerant of large gaps and insertions. Molecules are broken into overlapping fragments, which are characterized by their memberships in a probabilistic classification based on local geometry and H-bonding descriptors. This leads to a probabilistic similarity measure that is used in a conventional dynamic programming method.

RESULTS

Examples are given of database searching, the detection of structural similarities, which would not be found using sequence based methods, and comparisons with a previously published approach.

AVAILABILITY AND IMPLEMENTATION

Source code (C and perl) and binaries for linux are freely available at www.zbh.uni-hamburg.de/fries.

Multiple conformations of SAM-II riboswitch detected with SAXS and NMR spectroscopy

Riboswitches are a newly discovered large family of structured functional RNA elements that specifically bind small molecule targets out of a myriad of cellular metabolites to modulate gene expression. Structural studies of ligand-bound riboswitches by X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy have provided insights into detailed RNA-ligand recognition and interactions. However, the structures of ligand-free riboswitches remain poorly characterized. In this study, we have used a variety of biochemical, biophysical and computational techniques including small-angle X-ray scattering and NMR spectroscopy to characterize the ligand-free and ligand-bound forms of SAM-II riboswitch. Our data demonstrate that the RNA adopts multiple conformations along its folding pathway and suggest that the RNA undergoes marked conformational changes upon Mg(2+) compaction and S-adenosylmethionine (SAM) metabolite binding. Further studies indicated that Mg(2+) ion is not essential for the ligand binding but can stabilize the complex by facilitating loop/stem interactions. In the presence of millimolar concentration of Mg(2+) ion, the RNA samples a more compact conformation. This conformation is near to, but distinct from, the native fold and competent to bind the metabolite. We conclude that the formation of various secondary and tertiary structural elements, including a pseudoknot, occur to sequester the putative Shine-Dalgarno sequence of the RNA only after metabolite binding.

FINDING NON-CODING RNAs THROUGH GENOME-SCALE CLUSTERING

Non-coding RNAs (ncRNAs) are transcripts that do not code for proteins. Recent findings have shown that RNA-mediated regulatory mechanisms influence a substantial portion of typical microbial genomes. We present an efficient method for finding potential ncRNAs in bacteria by clustering genomic sequences based on homology inferred from both primary sequence and secondary structure. We evaluate our approach using a set of predominantly Firmicutes sequences. Our results showed that, though primary sequence based–homology search was inaccurate for diverged ncRNA sequences, through our clustering method, we were able to infer motifs that recovered nearly all members of most known ncRNA families. Hence, our method shows promise for discovering new families of ncRNA.

Differential analogue binding by two classes of c-di-GMP riboswitches

The ability of bacteria to adapt to a changing environment is essential for their survival. One mechanism bacteria have evolved to sense environmental cues and translate these signals into phenotypic changes uses the second messenger signaling molecule, cyclic diguanosine monophosphate (c-di-GMP). In addition to several classes of protein receptors, two classes of c-di-GMP-binding riboswitches (class I and class II) have been identified as downstream targets of the second messenger in this signaling pathway. The crystal structures of both riboswitch classes bound to c-di-GMP were previously reported. Here, we further investigate the mechanisms that RNA has evolved for recognition and binding of this second messenger. Using a series of c-di-GMP analogues, we probed the interactions made in the RNA-ligand complex for both classes of riboswitches to identify the most critical elements of c-di-GMP for binding. We found that the structural features of c-di-GMP required for binding differ between these two effectors and that the class II riboswitch is much less discriminatory in ligand binding than the class I riboswitch. These data suggest an explanation for the predicted preferential use of the class I motif over the class II motif in the c-di-GMP signaling pathway.

Synthetic circuits, devices and modules

The aim of synthetic biology is to design artificial biological systems for novel applications. From an engineering perspective, construction of biological systems of defined functionality in a hierarchical way is fundamental to this emerging field. Here, we highlight some current advances on design of several basic building blocks in synthetic biology including the artificial gene control elements, synthetic circuits and their assemblies into devices and modules. Such engineered basic building blocks largely expand the synthetic toolbox and contribute to our understanding of the underlying design principles of living cells.

Riboswitches as hormone receptors: hypothetical cytokinin-binding riboswitches in Arabidopsis thaliana

Background

Riboswitches are mRNA elements that change conformation when bound to small molecules. They are known to be key regulators of biosynthetic pathways in both prokaryotes and eukaryotes.

Presentation of the Hypothesis

The hypothesis presented here is that riboswitches function as receptors in hormone perception. We propose that riboswitches initiate or integrate signaling cascades upon binding to classic signaling molecules. The molecular interactions for ligand binding and gene expression control would be the same as for biosynthetic pathways, but the context and the cadre of ligands to consider is dramatically different. The hypothesis arose from the observation that a compound used to identify adenine binding RNA sequences is chemically similar to the classic plant hormone, or growth regulator, cytokinin. A general tenet of the hypothesis is that riboswitch-binding metabolites can be used to make predictions about chemically related signaling molecules. In fact, all cell permeable signaling compounds can be considered as potential riboswitch ligands. The hypothesis is plausible, as demonstrated by a cursory review of the transcriptome and genome of the model plant Arabidopsis thaliana for transcripts that i) contain an adenine aptamer motif, and ii) are also predicted to be cytokinin-regulated. Here, one gene, CRK10 (for Cysteine-rich Receptor-like Kinase 10, At4g23180), contains an adenine aptamer-related sequence and is down-regulated by cytokinin approximately three-fold in public gene expression data. To illustrate the hypothesis, implications of cytokinin-binding to the CRK10 mRNA are discussed.

Testing the hypothesis

At the broadest level, screening various cell permeable signaling molecules against random RNA libraries and comparing hits to sequence and gene expression data bases could determine how broadly the hypothesis applies. Specific cases, such as CRK10 presented here, will require experimental validation of direct ligand binding, altered RNA conformation, and effect on gene expression. Each case will be different depending on the signaling pathway and the physiology involved.

Implications of the hypothesis

This would be a very direct signal perception mechanism for regulating gene expression; rivaling animal steroid hormone receptors, which are frequently ligand dependent transcription initiation factors. Riboswitch-regulated responses could occur by modulating target RNA stability, translatability, and alternative splicing - all known expression platforms used in riboswitches. The specific illustration presented, CRK10, implies a new mechanism for the perception of cytokinin, a classic plant hormone. Experimental support for the hypothesis would add breadth to the growing list of important functions attributed to riboswitches.

Reviewers

This article was reviewed by Anthony Poole, Rob Knight, Mikhail Gelfand.

Protonation states of the key active site residues and structural dynamics of the glmS riboswitch as revealed by molecular dynamics

The glmS catalytic riboswitch is part of the 5'-untranslated region of mRNAs encoding glucosamine-6-phosphate (GlcN6P) synthetase (glmS) in numerous gram-positive bacteria. Binding of the cofactor GlcN6P induces site-specific self-cleavage of the RNA. However, the detailed reaction mechanism as well as the protonation state of the glmS reactive form still remains elusive. To probe the dominant protonation states of key active site residues, we carried out explicit solvent molecular dynamic simulations involving various protonation states of three crucial active site moieties observed in the available crystal structures: (i) guanine G40 (following the Thermoanaerobacter tengcongensis numbering), (ii) the GlcN6P amino/ammonium group, and (iii) the GlcN6P phosphate moiety. We found that a deprotonated G40(-) seems incompatible with the observed glmS active site architecture. Our data suggest that the canonical form of G40 plays a structural role by stabilizing an in-line attack conformation of the cleavage site A-1(2'-OH) nucleophile, rather than a more direct chemical role. In addition, we observe weakened cofactor binding upon protonation of the GlcN6P phosphate moiety, which explains the experimentally observed increase in K(m) with decreasing pH. Finally, we discuss a possible role of cofactor binding and its interaction with the G65 and G1 purines in structural stabilization of the A-1(2'-OH) in-line attack conformation. On the basis of the identified dominant protonation state of the reaction precursor, we propose a hypothesis of the self-cleavage mechanism in which A-1(2'-OH) is activated as a nucleophile by the G1(pro-R(p)) nonbridging oxygen of the scissile phosphate, whereas the ammonium group of GlcN6P acts as the general acid protonating the G1(O5') leaving group.

TRANSAT-- method for detecting the conserved helices of functional RNA structures, including transient, pseudo-knotted and alternative structures

The prediction of functional RNA structures has attracted increased interest, as it allows us to study the potential functional roles of many genes. RNA structure prediction methods, however, assume that there is a unique functional RNA structure and also do not predict functional features required for in vivo folding. In order to understand how functional RNA structures form in vivo, we require sophisticated experiments or reliable prediction methods. So far, there exist only a few, experimentally validated transient RNA structures. On the computational side, there exist several computer programs which aim to predict the co-transcriptional folding pathway in vivo, but these make a range of simplifying assumptions and do not capture all features known to influence RNA folding in vivo. We want to investigate if evolutionarily related RNA genes fold in a similar way in vivo. To this end, we have developed a new computational method, Transat, which detects conserved helices of high statistical significance. We introduce the method, present a comprehensive performance evaluation and show that Transat is able to predict the structural features of known reference structures including pseudo-knotted ones as well as those of known alternative structural configurations. Transat can also identify unstructured sub-sequences bound by other molecules and provides evidence for new helices which may define folding pathways, supporting the notion that homologous RNA sequence not only assume a similar reference RNA structure, but also fold similarly. Finally, we show that the structural features predicted by Transat differ from those assuming thermodynamic equilibrium. Unlike the existing methods for predicting folding pathways, our method works in a comparative way. This has the disadvantage of not being able to predict features as function of time, but has the considerable advantage of highlighting conserved features and of not requiring a detailed knowledge of the cellular environment.

Many non-coding genes exert their function via an RNA structure which starts emerging while the RNA sequence is being transcribed from the genome. The resulting folding pathway is known to depend on a variety of features such as the transcription speed, the concentration of various ions and the binding of proteins and other molecules. Not all of these influences can be adequately captured by the existing computational methods which try to replicate what happens in vivo. So far, it has been challenging to experimentally investigate co-transcriptional folding pathways in vivo and only little data from in vitro experiments exists. In order to investigate if functionally similar RNA sequences from different organisms fold in a similar way, we have developed a new computational method, called Transat, which does not require the detailed computational modeling of the cellular environment. We show in a comprehensive analysis that our method is capable of detecting known structural features and provide evidence that structural features of the in vivo folding pathways have been conserved for several biologically interesting classes of RNA sequences.

Review of microarray studies for host-intracellular pathogen interactions

Obligate intracellular bacteria are privileged soldiers on the battlefield that represent host-pathogen interactions. Microarrays are a powerful technology that can increase our knowledge about how bacteria respond to and interact with their hosts. This review summarizes the limitations inherent to host-pathogen interaction studies and essential strategies to improve microarray investigations of intracellular bacteria. We have compiled the comparative genomic and gene expression analyses of obligate intracellular bacteria currently available from microarrays. In this review we explore ways in which microarrays can be used to identify polymorphisms in different obligate intracellular bacteria such as Coxiella burnetii, Chlamydia trachomatis, Ehrlichia chaffeensis, Rickettsia prowazekii and Tropheryma whipplei. These microarray studies reveal that, while genomic content is highly conserved in obligate intracellular bacteria, genetic polymorphisms can potentially occur to increase bacterial pathogenesis. Additionally, changes in the gene expression of C. trachomatis throughout its life cycle, as well as changes in the gene expression profile of the pathogens R. prowazekii, Rickettsia rickettsii, Rickettsia typhi, T. whipplei and C. trachomatis following environmental changes, are discussed. Finally, an in vivo model of Rickettsia conorii within the skin is discussed. The gene expression analyses highlight the capacity of obligate intracellular bacteria to adapt to environmental changes and potentially to thwart the host response.

Genome-wide detection of predicted non-coding RNAs in Rhizobium etli expressed during free-living and host-associated growth using a high-resolution tiling array

Background

Non-coding RNAs (ncRNAs) play a crucial role in the intricate regulation of bacterial gene expression, allowing bacteria to quickly adapt to changing environments. In the past few years, a growing number of regulatory RNA elements have been predicted by computational methods, mostly in well-studied γ-proteobacteria but lately in several α-proteobacteria as well. Here, we have compared an extensive compilation of these non-coding RNA predictions to intergenic expression data of a whole-genome high-resolution tiling array in the soil-dwelling α-proteobacterium Rhizobium etli.

Results

Expression of 89 candidate ncRNAs was detected, both on the chromosome and on the six megaplasmids encompassing the R. etli genome. Of these, 11 correspond to functionally well characterized ncRNAs, 12 were previously identified in other α-proteobacteria but are as yet uncharacterized and 66 were computationally predicted earlier but had not been experimentally identified and were therefore classified as novel ncRNAs. The latter comprise 17 putative sRNAs and 49 putative cis-regulatory ncRNAs. A selection of these candidate ncRNAs was validated by RT-qPCR, Northern blotting and 5' RACE, confirming the existence of 4 ncRNAs. Interestingly, individual transcript levels of numerous ncRNAs varied during free-living growth and during interaction with the eukaryotic host plant, pointing to possible ncRNA-dependent regulation of these specialized processes.

Conclusions

Our data support the practical value of previous ncRNA prediction algorithms and significantly expand the list of candidate ncRNAs encoded in the intergenic regions of R. etli and, by extension, of α-proteobacteria. Moreover, we show high-resolution tiling arrays to be suitable tools for studying intergenic ncRNA transcription profiles across the genome. The differential expression levels of some of these ncRNAs may indicate a role in adaptation to changing environmental conditions.

Discrimination between closely related cellular metabolites by the SAM-I riboswitch

The SAM-I riboswitch is a cis-acting element of genetic control found in bacterial mRNAs that specifically binds S-adenosylmethionine (SAM). We previously determined the 2.9-A X-ray crystal structure of the effector-binding domain of this RNA element, revealing details of RNA-ligand recognition. To improve this structure, variations were made to the RNA sequence to alter lattice contacts, resulting in a 0.5-A improvement in crystallographic resolution and allowing for a more accurate refinement of the crystallographic model. The basis for SAM specificity was addressed by a structural analysis of the RNA complexed to S-adenosylhomocysteine (SAH) and sinefungin and by measuring the affinity of SAM and SAH for a series of mutants using isothermal titration calorimetry. These data illustrate the importance of two universally conserved base pairs in the RNA that form electrostatic interactions with the positively charged sulfonium group of SAM, thereby providing a basis for discrimination between SAM and SAH.

RNA phosphodiester backbone dynamics of a perdeuterated cUUCGg tetraloop RNA from phosphorus-31 NMR relaxation analysis

We have analyzed the relaxation properties of all (31)P nuclei in an RNA cUUCGg tetraloop model hairpin at proton magnetic field strengths of 300, 600 and 900 MHz in solution. Significant H, P dipolar contributions to R (1) and R (2) relaxation are observed in a protonated RNA sample at 600 MHz. These contributions can be suppressed using a perdeuterated RNA sample. In order to interpret the (31)P relaxation data (R (1), R (2)), we measured the (31)P chemical shift anisotropy (CSA) by solid-state NMR spectroscopy under various salt and hydration conditions. A value of 178.5 ppm for the (31)P CSA in the static state (S (2) = 1) could be determined. In order to obtain information about fast time scale dynamics we performed a modelfree analysis on the basis of our relaxation data. The results show that subnanosecond dynamics detected around the phosphodiester backbone are more pronounced than the dynamics detected for the ribofuranosyl and nucleobase moieties of the individual nucleotides (Duchardt and Schwalbe, J Biomol NMR 32:295-308, 2005; Ferner et al., Nucleic Acids Res 36:1928-1940, 2008). Furthermore, the dynamics of the individual phosphate groups seem to be correlated to the 5' neighbouring nucleobases.

A novel virtual spectrometry: Visualized regulatory motifs on ADM, rPol{beta} and CD83 mRNAs in human-friendly manners

Recently, riboswitches and other structures discovered on mRNAs have been reported as examples of functional RNA structures, motifs. Such motifs were shown to be present as single-form valid structures but they are obscured among other less-valid structures. Here, I present a novel, practical virtual spectrometry (the GenoPoemics Spectrometry) visualizing motifs on mRNA strands as spectra at-a-glance. Every motif along with validity of their existences could be observed on the spectra in human-friendly manners, and whole structures of mRNAs could be overviewed. Therefore, the spectra helped distinguish valid and less valid motifs. The spectrometry was applied to variety of mRNAs such as ADM, rPolbeta and CD83 to identify structures of high validity on them, previously reported functional motifs were successfully revealed. These findings indicate that the structures of mRNAs that may be folded into multiple forms can be further discussed quantitatively based on the visual spectra to discover functional RNA motifs.

Molecular dynamics simulation study of the binding of purine bases to the aptamer domain of the guanine sensing riboswitch

Riboswitches are a novel class of genetic control elements that function through the direct interaction of small metabolite molecules with structured RNA elements. The ligand is bound with high specificity and affinity to its RNA target and induces conformational changes of the RNA's secondary and tertiary structure upon binding. To elucidate the molecular basis of the remarkable ligand selectivity and affinity of one of these riboswitches, extensive all-atom molecular dynamics simulations in explicit solvent (approximately 1 micros total simulation length) of the aptamer domain of the guanine sensing riboswitch are performed. The conformational dynamics is studied when the system is bound to its cognate ligand guanine as well as bound to the non-cognate ligand adenine and in its free form. The simulations indicate that residue U51 in the aptamer domain functions as a general docking platform for purine bases, whereas the interactions between C74 and the ligand are crucial for ligand selectivity. These findings either suggest a two-step ligand recognition process, including a general purine binding step and a subsequent selection of the cognate ligand, or hint at different initial interactions of cognate and noncognate ligands with residues of the ligand binding pocket. To explore possible pathways of complex dissociation, various nonequilibrium simulations are performed which account for the first steps of ligand unbinding. The results delineate the minimal set of conformational changes needed for ligand release, suggest two possible pathways for the dissociation reaction, and underline the importance of long-range tertiary contacts for locking the ligand in the complex.

The Listeria transcriptional landscape from saprophytism to virulence

The bacterium Listeria monocytogenes is ubiquitous in the environment and can lead to severe food-borne infections. It has recently emerged as a multifaceted model in pathogenesis. However, how this bacterium switches from a saprophyte to a pathogen is largely unknown. Here, using tiling arrays and RNAs from wild-type and mutant bacteria grown in vitro, ex vivo and in vivo, we have analysed the transcription of its entire genome. We provide the complete Listeria operon map and have uncovered far more diverse types of RNAs than expected: in addition to 50 small RNAs (<500 nucleotides), at least two of which are involved in virulence in mice, we have identified antisense RNAs covering several open-reading frames and long overlapping 5' and 3' untranslated regions. We discovered that riboswitches can act as terminators for upstream genes. When Listeria reaches the host intestinal lumen, an extensive transcriptional reshaping occurs with a SigB-mediated activation of virulence genes. In contrast, in the blood, PrfA controls transcription of virulence genes. Remarkably, several non-coding RNAs absent in the non-pathogenic species Listeria innocua exhibit the same expression patterns as the virulence genes. Together, our data unravel successive and coordinated global transcriptional changes during infection and point to previously unknown regulatory mechanisms in bacteria.

Design principles for riboswitch function

Scientific and technological advances that enable the tuning of integrated regulatory components to match network and system requirements are critical to reliably control the function of biological systems. RNA provides a promising building block for the construction of tunable regulatory components based on its rich regulatory capacity and our current understanding of the sequence–function relationship. One prominent example of RNA-based regulatory components is riboswitches, genetic elements that mediate ligand control of gene expression through diverse regulatory mechanisms. While characterization of natural and synthetic riboswitches has revealed that riboswitch function can be modulated through sequence alteration, no quantitative frameworks exist to investigate or guide riboswitch tuning. Here, we combined mathematical modeling and experimental approaches to investigate the relationship between riboswitch function and performance. Model results demonstrated that the competition between reversible and irreversible rate constants dictates performance for different regulatory mechanisms. We also found that practical system restrictions, such as an upper limit on ligand concentration, can significantly alter the requirements for riboswitch performance, necessitating alternative tuning strategies. Previous experimental data for natural and synthetic riboswitches as well as experiments conducted in this work support model predictions. From our results, we developed a set of general design principles for synthetic riboswitches. Our results also provide a foundation from which to investigate how natural riboswitches are tuned to meet systems-level regulatory demands.

Riboswitches are RNA-based components that integrate ligand binding and gene regulation to dynamically respond to molecular signals within cells. Natural riboswitches are employed to regulate metabolism and other cellular processes, while synthetic riboswitches have been constructed to expand the sensory and regulatory capabilities exhibited in nature. Characterization studies have revealed that sequence modifications can tune properties of the riboswitch response curve, which links ligand concentration to expression levels. Tunability is critical when matching component properties to the regulatory demands of biological systems; however, the characterization of riboswitch tuning strategies is complicated by the integration of numerous regulatory mechanisms and various processes, such as RNA folding and turnover, that impact riboswitch performance. To develop a generalized framework for examining quantitative aspects of riboswitch tuning, we modeled the kinetics of riboswitch function operating under common regulatory mechanisms. Our results reveal that riboswitch performance is primarily dictated by the competition between reversible and mechanism-specific irreversible rate constants. We demonstrate that practical system restrictions can significantly alter the requirements for riboswitch performance, necessitating a variety of tuning strategies. We developed design principles to guide the construction of synthetic riboswitches and a quantitative framework from which to investigate how natural riboswitches are tuned to meet systems-level regulatory demands.

Riboswitches: ancient and promising genetic regulators

BAIT AND SWITCH: Metabolite-sensing riboswitches make use of RNA structural modulation to regulate gene expression, as illustrated in the scheme, in response to subtle changes in metabolite concentrations. This review describes the current knowledge about naturally occurring riboswitches and their growing potential as antibacterial cellular targets and as molecular biosensors. Newly discovered metabolite-sensing riboswitches have revealed that cellular processes extensively make use of RNA structural modulation to regulate gene expression in response to subtle changes in metabolite concentrations. Riboswitches are involved at various regulation levels of gene expression, such as transcription attenuation, translation initiation, mRNA splicing and mRNA processing. Riboswitches are found in the three kingdoms of life, and in various cases, are involved in the regulation of essential genes, which makes their regulation an essential part of cell survival. Because riboswitches operate without the assistance of accessory proteins, they are believed to be remnants of an ancient time, when gene regulation was strictly based on RNA, from which are left numerous "living molecular fossils", as exemplified by ribozymes, and more spectacularly, by the ribosome. Due to their nature, riboswitches hold high expectations for the manipulation of gene expression and the detection of small metabolites, and also offer an unprecedented potential for the discovery of novel classes of antimicrobial agents.

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.

RNA is remarkably versatile, acting not only as messengers to transfer genetic information from DNA to protein but also as critical structural components and catalytic enzymes in the cell. More intriguingly, RNA elements in messenger RNAs have been widely found in bacteria to control the expression of their downstream genes. The functions of these RNA elements are intrinsically linked to their secondary structures, which are usually conserved across multiple closely related species during evolution and often shared by genes in the same metabolic pathways. We developed a new computational approach to find putative functional RNA elements by looking for conserved RNA secondary structures that are distinguished from random RNA secondary structures in the orthologous RNA sequences from related species. We applied this approach to multiple Shewanella genomes and predicted putative regulatory RNA elements in Shewanella oneidensis, a bacterium that has extraordinary respiratory and metal reducing abilities and great potential for bioremediation and alternative energy generation. Our findings not only recovered many RNA elements that are known or supported by literature evidence but also included exciting novel RNA elements for further exploration.

Large intergenic cruciform-like supermotifs in the Lactobacillus plantarum genome

Twenty-four Lactobacillus plantarum supermotifs (LPSMs) with lengths from approximately 800 to 1,000 nucleotides were identified in the L. plantarum genome. LPSMs were conserved in other L. plantarum strains but not in other species. Secondary structure analysis predicted that LPSMs may fold into cruciform-like structures. Preliminary experiments indicate that the LPSMs are transcribed.

Methods for analysis of ligand-induced RNA conformational changes

Encoded within many RNA sequences is the requisite information for folding of intricate three-dimensional structures. Moreover, many noncoding RNAs can adopt structurally distinct and functionally specialized conformations in response to specific cellular signals. These conformational transitions are often times accompanied by changes in hydrodynamic radii. Therefore, experimental methods that measure changes in hydrodynamic radius can be employed for study of signal-induced RNA conformational changes. Several hydrodynamic methods, including analytical ultracentrifugation, size-exclusion chromatography, and nondenaturing gel electrophoresis, are briefly discussed herein.

Kinetics of riboswitch regulation studied by in vitro transcription

Riboswitches encompass messenger ribonucleic acid transcripts that sense the concentration of small molecule metabolites through binding the target compound and then control the expression of metabolite-related genes in response to the metabolite concentration. While much of the riboswitch-related research has focused on the remarkable ability of different aptamer domains to adopt the intricate structures required to bind a spectrum of biological metabolites with high affinity and specificity, less attention has been paid to the mechanism of riboswitch action. Specifically, the genetic control element of the riboswitch, known as the expression platform, must function cotranscriptionally in the case of transcription termination-controlled riboswitches. By correlating the transcriptional kinetics of the entire switch and the kinetics and thermodynamics of metabolite binding of the aptamer domain, it was found that the FMN-binding riboswitch in the 5' UTR of the Bacillus subtilis ribGBAHT operon functions as a kinetically controlled genetic switch chiefly dependent upon transcriptional pausing and the concentration of the target metabolite. This study has emphasized the importance of studying the switch in its entirety and in the context of an actively transcribing RNA polymerase. Herein I will describe the study of the kinetics of riboswitch transcription and the proposed mechanism for the transcription termination-associated riboswitch control of riboflavin-related genes.

Application of fluorescent measurements for characterization of riboswitch-ligand interactions

Riboswitches are recently discovered messenger RNA motifs involved in gene regulation. They modulate gene expression at various levels, such as transcription, translation, splicing, and mRNA degradation. Because riboswitches exhibit relatively complex structures, they are able to form highly complex ligand-binding sites, which enable the specific recognition of target metabolites in a complex cellular environment. Practically in all studied cases, riboswitches use ligand-induced conformational changes to control gene expression. To monitor the structural reorganization of riboswitches, we use the local fluorescent reporter 2-aminopurine (2AP), which is a structural analog of adenine. The 2AP fluorescence is strongly quenched when the fluorophore is involved in stacking interactions with surrounding bases, and can, therefore, be used to monitor local structural rearrangements. Here, we show specific examples in which 2AP fluorescence can be used to monitor structural changes in the aptamer domain of the S-adenosyl methionine (SAM) riboswitch and where it can be used as a ligand for the guanine riboswitch.

Synthetic biology with RNA motifs

Structural motifs in naturally occurring RNAs and RNPs can be employed as new molecular parts for synthetic biology to facilitate the development of novel devices and systems that modulate cellular functions. In this review, we focus on the following: (i) experimental evolution techniques of RNA molecules in vitro and (ii) their applications for regulating gene expression systems in vivo. For experimental evolution, new artificial RNA aptamers and RNA enzymes (ribozymes) have been selected in vitro. These functional RNA molecules are likely to be applicable in the reprogramming of existing gene regulatory systems. Furthermore, they may be used for designing hypothetical RNA-based living systems in the so-called RNA world. For the regulation of gene expressions in living cells, the development of new riboswitches allows us to modulate the target gene expression in a tailor-made manner. Moreover, recently RNA-based synthetic genetic circuits have been reported by employing functional RNA molecules, expanding the repertory of synthetic biology with RNA motifs.

Riboswitches: emerging themes in RNA structure and function

Riboswitches are RNAs capable of binding cellular metabolites using a diverse array of secondary and tertiary structures to modulate gene expression. The recent determination of the three-dimensional structures of parts of six different riboswitches illuminates common features that allow riboswitches to be grouped into one of two types. Type I riboswitches, as exemplified by the purine riboswitch, are characterized by a single, localized binding pocket supported by a largely pre-established global fold. This arrangement limits ligand-induced conformational changes in the RNA to a small region. In contrast, Type II riboswitches, such as the thiamine pyrophosphate riboswitch, contain binding pockets split into at least two spatially distinct sites. As a result, binding induces both local changes to the binding pocket and global architecture. Similar organizational themes are found in other noncoding RNAs, making it possible to begin to build a hierarchical classification of RNA structure based on the spatial organization of their active sites and associated secondary structural elements.