The regulation mechanism ofyitJandmetFriboswitches

Regulation of the ThiM riboswitch is facilitated by the trapped structure formed during transcription of the wild-type sequence

The ThiM riboswitch from Escherichia coli is a typical mRNA device that modulates downstream gene expression by sensing TPP. The helix-based RNA folding theory is used to investigate its detailed regulatory behaviors in cells. This RNA molecule is transcriptionally trapped in a state with the unstructured SD sequence in the absence of TPP, which induces downstream gene expression. As a key step to turn on gene expression, formation of this trapped state (the genetic ON state) highly depends on the co-transcriptional folding of its wild-type sequence. Instead of stabilities of the genetic ON and OFF states, the transcription rate, pause, and ligand levels are combined to affect the ThiM riboswitch-mediated gene regulation, which is consistent with a kinetic control model.

RNA-align: quick and accurate alignment of RNA 3D structures based on size-independent TM-scoreRNA

MOTIVATION

Comparison of RNA 3D structures can be used to infer functional relationship of RNA molecules. Most of the current RNA structure alignment programs are built on size-dependent scales, which complicate the interpretation of structure and functional relations. Meanwhile, the low speed prevents the programs from being applied to large-scale RNA structural database search.

RESULTS

We developed an open-source algorithm, RNA-align, for RNA 3D structure alignment which has the structure similarity scaled by a size-independent and statistically interpretable scoring metric. Large-scale benchmark tests show that RNA-align significantly outperforms other state-of-the-art programs in both alignment accuracy and running speed. The major advantage of RNA-align lies at the quick convergence of the heuristic alignment iterations and the coarse-grained secondary structure assignment, both of which are crucial to the speed and accuracy of RNA structure alignments.

AVAILABILITY AND IMPLEMENTATION

https://zhanglab.ccmb.med.umich.edu/RNA-align/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Regulation of the thiamine pyrophosphate (TPP)-sensing riboswitch in NMT1 mRNA from Neurospora crassa

The expression of Neurospora crassa NMT1 involved in thiamine pyrophosphate (TPP) metabolism is regulated at the level of mRNA splicing by a TPP-sensing riboswitch within the precursor NMT1 mRNA. Here, using the systematic helix-based computational method, we investigated the regulation of this riboswitch. We find that the function of the riboswitch does not depend on the transcription process. Whether TPP is present or not, the riboswitch predominately folds into the ON state, while the OFF state aptamer structure does not appear during transcription. Since the transition from the ON state to the aptamer structure is extremely slow, TPP may interact with the RNA before full formation of the aptamer structure, promoting the switch flipping. The potential to fully form helix P0 of the ON state is necessary to restore ligand-dependent gene control by the riboswitch.

Effects of flanking regions on HDV cotranscriptional folding kinetics

Hepatitis delta virus (HDV) ribozyme performs the self-cleavage activity through folding to a double pseudoknot structure. The folding of functional RNA structures is often coupled with the transcription process. In this work, we developed a new approach for predicting the cotranscriptional folding kinetics of RNA secondary structures with pseudoknots. We theoretically studied the cotranscriptional folding behavior of the 99-nucleotide (nt) HDV sequence, two upstream flanking sequences, and one downstream flanking sequence. During transcription, the 99-nt HDV can effectively avoid the trap intermediates and quickly fold to the cleavage-active state. It is different from its refolding kinetics, which folds into an intermediate trap state. For all the sequences, the ribozyme regions (from 1 to 73) all fold to the same structure during transcription. However, the existence of the 30-nt upstream flanking sequence can inhibit the ribozyme region folding into the active native state through forming an alternative helix Alt1 with the segments 70-90. The longer upstream flanking sequence of 54 nt itself forms a stable hairpin structure, which sequesters the formation of the Alt1 helix and leads to rapid formation of the cleavage-active structure. Although the 55-nt downstream flanking sequence could invade the already folded active structure during transcription by forming a more stable helix with the ribozyme region, the slow transition rate could keep the structure in the cleavage-active structure to perform the activity.

Genetic regulation mechanism of the yjdF riboswitch

The yjdF riboswitch resides in potential 5' UTRs of homologues of protein-coding gene yjdF in Firmicutes. Unlike other 30 riboswitch classes previously validated, this riboswitch class, can sense and bind to a broad collection of azaaromatic ligands. Among these compounds, some do activate production of yjdF protein driven by the riboswitch, while others are out of riboswitch-mediated modulation possibly because of the toxicity at high ligand concentrations. By incorporating the structures with pseudoknots and ligand binding kinetics into the co-transcriptional folding theory, we theoretically studied the co-transcriptional folding behaviors of the yjdF riboswitch from Bacillus subtilis at different transcription conditions. Like most riboswitches, the yjdF riboswitch can quickly fold into the aptamer structure without any trapped states during the transcription process. After the aptamer structure is formed, the riboswitch shows two main co-transcriptional folding pathways: aptamer→ON state→OFF state and aptamer → the ligand bound aptamer → the ligand bound ON state. Our results suggested that this translational riboswitch is coupled with the transcription process to exert its biological function and it is kinetically controlled. The threshold concentration for the ligand to activate the riboswitch depends on the transcription rate and the association rate of the ligand binding.

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.

Co-Transcriptional Folding and Regulation Mechanisms of Riboswitches

Riboswitches are genetic control elements within non-coding regions of mRNA. These self-regulatory elements have been found to sense a range of small metabolites, ions, and other physical signals to exert regulatory control of transcription, translation, and splicing. To date, more than a dozen riboswitch classes have been characterized that vary widely in size and secondary structure. Extensive experiments and theoretical studies have made great strides in understanding the general structures, genetic mechanisms, and regulatory activities of individual riboswitches. As the ligand-dependent co-transcriptional folding and unfolding dynamics of riboswitches are the key determinant of gene expression, it is important to investigate the thermodynamics and kinetics of riboswitches both in the presence and absence of metabolites under the transcription. This review will provide a brief summary of the studies about the regulation mechanisms of the pbuE, S_MK, yitJ, and metF riboswitches based on the ligand-dependent co-transcriptional folding of the 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.