RNA Secondary Structure Elements Research Articles

Predicting functional riboSNitches in the context of alternative splicing

RNAs are the major regulators of gene expression, and their secondary structures play crucial roles at different levels. RiboSNitches are disease-associated SNPs that cause changes in the pre-mRNA secondary structural ensemble. Several riboSNitches have been detected in the 5′ and 3′ untranslated regions and lncRNA. Although cases of secondary structural elements playing a regulatory role in alternative splicing are known, regions specific to splicing events, such as splice junctions have not received much attention. We tested splice-site mutations for their efficiency in disrupting the secondary structure and hypothesized that these could play a crucial role in alternative splicing. Multiple riboSNitch prediction methods were applied to obtain overlapping results that are potentially more reliable. Putative riboSNitches were identified from aberrant 5′ and 3′ splice site mutations, cancer-causing somatic mutations, and genes that harbor the regulatory RNA secondary structural elements. Our workflow for predicting riboSNitches associated with alternative splicing is novel and paves the way for subsequent experimental validation.

RNA Sequence and Structure Determinants of Pol III Transcriptional Termination in Human Cells

The precise mechanism of transcription termination of the eukaryotic RNA polymerase III (Pol III) has been a subject of considerable debate. Although previous studies have clearly shown that multiple uracils at the end of RNA transcripts are required for Pol III termination, the effects of upstream RNA secondary structure in the nascent transcript on transcriptional termination is still unclear. To address this, we developed an in cellulo Pol III transcription termination assay using the recently developed Tornado-Corn RNA aptamer system to create a Pol III-transcribed RNA that produces a detectable fluorescent signal when transcribed in human cells. To study the effects of RNA sequence and structure on Pol III termination, we systematically varied the sequence context upstream of the aptamer and identified sequence characteristics that enhance or diminish termination. For transcription from Pol III type 3 promoters, we found that only poly-U tracts longer than the average length found in the human genome efficiently terminate Pol III transcription without RNA secondary structure elements. We observed that RNA secondary structure elements placed in proximity to shorter poly-U tracts induced termination, and RNA secondary structure by itself was not sufficient to induce termination. For Pol III type 2 promoters, we found that the shorter poly-U tract lengths of 4 uracils were sufficient to induce termination. These findings demonstrate a key role for sequence and structural elements within Pol III-transcribed nascent RNA for efficient transcription termination, and demonstrate a generalizable assay for characterizing Pol III transcription in human cells.

A Mathematical Analysis of RNA Structural Motifs in Viruses

RNA stem-loop structures play an important role in almost every step of the viral replication cycle. In this contribution, a mathematical analysis is performed on a large dataset of RNA secondary structure elements in the coding regions of viruses by using topological indices that capture the Laplacian eigenvalues of the associated RNA graph representations and thereby enable structural classification, supplemented by folding energy and mutational robustness. The application of such an analysis for viral RNA structural motifs is described, being able to extract structural categories such as stem-loop structures of different sizes according to the tree-graph representation of the RNA structure, in our attempt to find novel functional motifs. While the analysis is carried on a large dataset of viral RNA structures, it can be applied more generally to other data that involve RNA secondary structures in biological agents.

Effect of the Cation Size on the Stability of RNA Structures

The folding and subsequent stabilization of functional RNA structures depend strongly on the electrostatic screening of the negative charges on the phosphate group by alkali and alkaline earth metal ions. Such cations can either bind to structurally well-defined RNA sites or form a so-called ion cloud around RNA. However, it is difficult to directly observe these ions experimentally. The common view is that sodium ion is better in stabilizing RNA secondary structural elements than potassium ion that possess the same valance but a larger hydrodynamic radius.

G2-quadruplex in the 3’UTR of IE180 regulates Pseudorabies virus replication by enhancing gene expression

ABSTRACT RNA secondary structure elements in the mRNA 3ʹ-untranslated regions (3ʹUTR) play important roles in post-transcriptional regulation. RNA structure elements in the viral RNA provide valuable model for studying diverse regulation mechanisms. Herpesvirus genomes are double-stranded DNA with GC-rich sequences, which can be transcribed into abundant GC-rich RNAs. It is valuable to explore the structures and function of those GC-rich RNAs. We identified a G2-quadruplex-forming sequence named PQS18-1 in the 3ʹUTR of the unique immediate early gene of Pseudorabies virus (PRV), an important member of Alphaherpesvirinae subfamily. The RNA PQS18-1 was folded into parallel G-quadruplex structure, enhancing gene expression. Both non-G-quadruplex mutant and G3-quadruplex mutant in the 3ʹUTR showed lower gene expression level than the wildtype G2-quadruplex. TMPyP4 destroyed PQS18-1 G2-quadruplex and suppressed gene expression, accordingly reducing PRV replication by one titre in the PK15 cells at 24 h post infection. Our findings indicated that the RNA G2-quadruplex in 3ʹUTR was essential for high expression of IE180 gene, and it could be a specific post-transcription regulation element in response to small molecules or other macromolecules. This study discovers a novel RNA G2-quadruplex in the 3ʹUTR of an immediate early gene of alphaherpesvirus and provides a new nucleic acid target for anti-virus drug design.

3D based on 2D: Calculating helix angles and stacking patterns using forgi 2.0, an RNA Python library centered on secondary structure elements.

We present forgi, a Python library to analyze the tertiary structure of RNA secondary structure elements. Our representation of an RNA molecule is centered on secondary structure elements (stems, bulges and loops). By fitting a cylinder to the helix axis, these elements are carried over into a coarse-grained 3D structure representation. Integration with Biopython allows for handling of all-atom 3D information. forgi can deal with a variety of file formats including dotbracket strings, PDB and MMCIF files. We can handle modified residues, missing residues, cofold and multifold structures as well as nucleotide numbers starting at arbitrary positions. We apply this library to the study of stacking helices in junctions and pseudoknots and investigate how far stacking helices in solved experimental structures can divert from coaxial geometries.

Interplay of primary sequence, position and secondary RNA structure determines alternative splicing of LMNA in a pre-mature aging syndrome.

Aberrant splicing in exon 11 of the LMNA gene causes the premature aging disorder Hutchinson-Gilford Progeria Syndrome. A de novo C1824T mutation activates an internal alternative 5' splice site, resulting in formation of the disease-causing progerin protein. The underlying mechanism for this 5' splice site selection is unknown. Here, we have applied a combination of targeted mutational analysis in a cell-based system and structural mapping by SHAPE-MaP to comprehensively probe the contributions of primary sequence, secondary RNA structure and linear splice site position in determining in vivo mechanisms of splice site choice inLMNA. While splice site choice is in part defined by sequence complementarity to U1 snRNA, we identify RNA secondary structural elements near the alternative 5' splice sites and show that splice site choice is significantly influenced by the structural context of the available splice sites. Furthermore, relative positioning of the competing sites within the primary sequence of the pre-mRNA is a predictor of 5' splice site usage, with the distal position favored over the proximal, regardless of sequence composition. Together, these results demonstrate that 5' splice site selection in LMNA is determined by an intricate interplay among RNA sequence, secondary structure and splice site position.

3D based on 2D: Calculating helix angles and stacking patterns using forgi 2.0, an RNA Python library centered on secondary structure elements.

We present forgi, a Python library to analyze the tertiary structure of RNA secondary structure elements. Our representation of an RNA molecule is centered on secondary structure elements (stems, bulges and loops). By fitting a cylinder to the helix axis, these elements are carried over into a coarse-grained 3D structure representation. Integration with Biopython allows for handling of all-atom 3D information. forgi can deal with a variety of file formats including dotbracket strings, PDB and MMCIF files. We can handle modified residues, missing residues, cofold and multifold structures as well as nucleotide numbers starting at arbitrary positions. We apply this library to the study of stacking helices in junctions and pseudoknots and investigate how far stacking helices in solved experimental structures can divert from coaxial geometries.

3D based on 2D: Calculating helix angles and stacking patterns using forgi 2.0, an RNA Python library centered on secondary structure elements.

We present forgi, a Python library to analyze the tertiary structure of RNA secondary structure elements. Our representation of an RNA molecule is centered on secondary structure elements (stems, bulges and loops). By fitting a cylinder to the helix axis, these elements are carried over into a coarse-grained 3D structure representation. Integration with Biopython allows for handling of all-atom 3D information. forgi can deal with a variety of file formats including dotbracket strings, PDB and MMCIF files. We can handle modified residues, missing residues, cofold and multifold structures as well as nucleotide numbers starting at arbitrary positions. We apply this library to the study of stacking helices in junctions and pseudo knots and investigate how far stacking helices in solved experimental structures can divert from coaxial geometries.

Highly conserved nucleotide motifs present in the 5'UTR of the heme-receptor gene shmR are required for HmuP-dependent expression of shmR in Ensifer meliloti.

Heme may represent a major iron-source for bacteria. In the symbiotic nitrogen-fixing bacterium Ensifer meliloti 1021, iron acquisition from heme depends on the outer-membrane heme-receptor ShmR. Expression of shmR gene is repressed by iron in a RirA dependent manner while under iron-limitation its expression requires the small protein HmuP. In this work, we identified highly conserved nucleotide motifs present upstream the shmR gene. These motifs are widely distributed among Alpha and Beta Proteobacteria, and correlate with the presence of HmuP coding sequences in bacterial genomes. According to data presented in this work, we named these new motifs as HmuP-responsive elements (HPREs). In the analyzed genomes, the HPREs were always present upstream of genes encoding putative heme-receptors. Moreover, in those Alpha and Beta Proteobacteria where transcriptional start sites for shmR homologs are known, HPREs were located in the 5'UTR region. In this work we show that in E. meliloti 1021, HPREs are involved in HmuP-dependent shmR expression. Moreover, we show that changes in sequence composition of the HPREs correlate with changes in a predicted RNA secondary structure element and affect shmR gene expression.

In vivo probing of nascent RNA structures reveals principles of cotranscriptional folding.

Defining the in vivo folding pathway of cellular RNAs is essential to understand how they reach their final native conformation. We here introduce a novel method, named Structural Probing of Elongating Transcripts (SPET-seq), that permits single-base resolution analysis of transcription intermediates’ secondary structures on a transcriptome-wide scale, enabling base-resolution analysis of the RNA folding events. Our results suggest that cotranscriptional RNA folding in vivo is a mixture of cooperative folding events, in which local RNA secondary structure elements are formed as they get transcribed, and non-cooperative events, in which 5′-halves of long-range helices get sequestered into transient non-native interactions until their 3′ counterparts have been transcribed. Together our work provides the first transcriptome-scale overview of RNA cotranscriptional folding in a living organism.

Elucidating the role of essential RNA secondary structural elements in dengue biology and their implication in dengue virulence

Elucidating the role of essential RNA secondary structural elements in dengue biology and their implication in dengue virulence

RNA World research-still evolving.

RNA World research-still evolving.

Evolution of helix formation in the ribosomal Internal Transcribed Spacer 2 (ITS2) and its significance for RNA secondary structures.

Helices are the most common elements of RNA secondary structure. Despite intensive investigations of various types of RNAs, the evolutionary history of the formation of new helices (novel helical structures) remains largely elusive. Here, by studying the nuclear ribosomal Internal Transcribed Spacer 2 (ITS2), a fast-evolving part of the eukaryotic nuclear ribosomal operon, we identify two possible types of helix formation: one type is "dichotomous helix formation"--transition from one large helix to two smaller helices by invagination of the apical part of a helix, which significantly changes the shape of the original secondary structure but does not increase its complexity (i.e., the total length of the RNA). An alternative type is "lateral helix formation"--origin of an extra helical region by the extension of a bulge loop or a spacer in a multi-helix loop of the original helix, which does not disrupt the pre-existing structure but increases RNA size. Moreover, we present examples from the RNA sequence literature indicating that both types of helix formation may have implications for RNA evolution beyond ITS2.

A novel mechanism for the control of translation of specific mRNAs by tumor suppressor protein Pdcd4: inhibition of translation elongation

The tumor suppressor gene Pdcd4 (programmed cell death gene 4) has drawn considerable attention because its downregulation is involved in the development of several types of cancer. Because Pdcd4 interacts with the translation initiation factor eIF4A and inhibits its helicase activity, Pdcd4 has been implicated in the translational suppression of cellular mRNAs containing structured 5'-untranslated regions. However, Pdcd4's role in translation regulation is still poorly understood, because only very few physiological Pdcd4 target mRNAs are known. By using a Pdcd4-deficient clone of the chicken B-cell line DT40, we have discovered that the mRNA of the A-myb proto-oncogene is a novel Pdcd4 target RNA whose translation is suppressed by Pdcd4. Interestingly, the inhibitory effect of Pdcd4 is independent of the Pdcd4-eIF4A interaction, but is dependent on an RNA-binding domain at the N terminus of Pdcd4 and on sequences located within the coding region of A-myb mRNA, indicating that Pdcd4 suppresses A-myb translation by a novel mechanism. Our data show that the Pdcd4 RNA-binding domain preferentially recognizes an RNA secondary structure element formed by the part of the A-myb coding region that mediates Pdcd4-dependent suppression. Previously, we have shown that Pdcd4 also suppresses the translation of the c-myb mRNA by a similar mechanism involving binding of Pdcd4 to RNA secondary structure formed by the c-myb coding region. Surprisingly, our data show that Pdcd4 exerts its inhibitory activity only when the target region of Pdcd4 in A-myb and c-myb mRNA is itself translated, consistent with a mechanism in which Pdcd4 suppresses translation by interfering with translation elongation. Taken together, our work reveals a novel mechanism by which Pdcd4 affects the translational of cellular RNAs. Furthermore, as c-myb and A-myb are members of the Myb proto-oncogene family whose deregulation has been implicated in tumorigenesis, inhibiting their translation might contribute to the tumor-suppressive activity of Pdcd4.

Expansion of the mutually exclusive spliced exome in Drosophila

Mutually exclusive splicing is an important mechanism in a wide range of eukaryotic branches to expand proteome diversity, but the extent of its distribution within a single species and its evolutionary conservation is unknown. Here we present a genome-wide analysis of mutually exclusive spliced exons (MXEs) in Drosophila melanogaster at unprecedented depth. Most of the new MXE candidates are supported by evolutionary conservation, transcriptome data analysis and identification of competing RNA secondary structural elements. The enrichment of the genes with MXEs in transmembrane transporters and ion channel activity is consistent with findings in humans, although the MXEs appeared independently and in non-homologous genes, supporting the idea of a universal benefit of adapting ion channel and receptor properties by tandem exon duplications. The comparison of the mutually exclusive spliced exomes within the Drosophila clade shows high numbers of MXE gain and loss events, suggesting a role of these processes in speciation.

Messenger RNA Sequence Rather than Protein Sequence Determines the Level of Self-synthesis and Antigen Presentation of the EBV-encoded Antigen, EBNA1

Unique purine-rich mRNA sequences embedded in the coding sequences of a distinct group of gammaherpesvirus maintenance proteins underlie the ability of the latently infected cell to minimize immune recognition. The Epstein-Barr virus nuclear antigen, EBNA1, a well characterized lymphocryptovirus maintenance protein has been shown to inhibit in cis antigen presentation, due in part to a large internal repeat domain encoding glycine and alanine residues (GAr) encoded by a purine-rich mRNA sequence. Recent studies have suggested that it is the purine-rich mRNA sequence of this repeat region rather than the encoded GAr polypeptide that directly inhibits EBNA1 self-synthesis and contributes to immune evasion. To test this hypothesis, we generated a series of EBNA1 internal repeat frameshift constructs and assessed their effects on cis-translation and endogenous antigen presentation. Diverse peptide sequences resulting from alternative repeat reading frames did not alleviate the translational inhibition characteristic of EBNA1 self-synthesis or the ensuing reduced surface presentation of EBNA1-specific peptide-MHC class I complexes. Human cells expressing the EBNA1 frameshift variants were also poorly recognized by antigen-specific T-cells. Furthermore, a comparative analysis of the mRNA sequences of the corresponding repeat regions of different viral maintenance homologues highlights the high degree of identity between the nucleotide sequences despite very little homology in the encoded amino acid sequences. Based on these combined observations, we propose that the cis-translational inhibitory effect of the EBNA1 internal repeat sequence operates mechanistically at the nucleotide level, potentially through RNA secondary structural elements, and is unlikely to be mediated through the GAr polypeptide. The demonstration that the EBNA1 repeat mRNA sequence and not the encoded protein sequence underlies immune evasion in this class of virus suggests a novel approach to therapeutic development through the use of anti-sense strategies or small molecules targeting EBNA1 mRNA structure.

Small‐Molecule‐Mediated Cleavage of RNA in Living Cells

Small‐Molecule‐Mediated Cleavage of RNA in Living Cells

Automated 3D structure composition for large RNAs

Understanding the numerous functions that RNAs play in living cells depends critically on knowledge of their three-dimensional structure. Due to the difficulties in experimentally assessing structures of large RNAs, there is currently great demand for new high-resolution structure prediction methods. We present the novel method for the fully automated prediction of RNA 3D structures from a user-defined secondary structure. The concept is founded on the machine translation system. The translation engine operates on the RNA FRABASE database tailored to the dictionary relating the RNA secondary structure and tertiary structure elements. The translation algorithm is very fast. Initial 3D structure is composed in a range of seconds on a single processor. The method assures the prediction of large RNA 3D structures of high quality. Our approach needs neither structural templates nor RNA sequence alignment, required for comparative methods. This enables the building of unresolved yet native and artificial RNA structures. The method is implemented in a publicly available, user-friendly server RNAComposer. It works in an interactive mode and a batch mode. The batch mode is designed for large-scale modelling and accepts atomic distance restraints. Presently, the server is set to build RNA structures of up to 500 residues.

Effect of Neural Network Parameters on RNA Secondary Structure Classification

Helix, Hairpin, Bulge, external loop, internal loop, multi-branch loop are the elements of RNA secondary structure. We have designed a neural network to classify the RNA sequence in to three categories i.e Hairpin, helix, neither of two. This can be extended to classify into all secondary structure elements. If all the elements are predicted then we can determine the entire structure of a RNA family. The parameters of neural network affect the performance of the network. But there are no rules to define the value of these parameters of network. For a given problem the optimal value of parameters can be obtained by performing the experiments on their values. This paper shows the effect on the performance of classification by varying the number of hidden layers, number of neurons and activation functions.