2D meets 4G: G-quadruplexes in RNA secondary structure prediction

Thermodynamics and Kinetics-Directed Regulation of Nucleic Acid-Based Molecular Recognition

Nucleic acid-based molecular recognition plays crucial roles in various fields like biosensing and disease diagnostics. To achieve optimal detection and analysis, it is essential to regulate the response performance of nucleic acid probes or switches to match specific application requirements by regulating thermodynamics and kinetics properties. However, the impacts of thermodynamics and kinetics theories on recognition performance are sometimes obscure and the relative conclusions are not intuitive. To promote the thorough understanding and rational utilization of thermodynamics and kinetics theories, this review focuses on the landmarks and recent advances of nucleic acid thermodynamics and kinetics and summarizes the nucleic acid thermodynamics and kinetics-based strategies for regulation of nucleic acid-based molecular recognition. This work hopes such a review can provide reference and guidance for the development and optimization of nucleic acid probes and switches in the future, as well as for advancements in other nucleic acid-related fields.

In Silico Identification of Potential Quadruplex Forming Sequences in LncRNAs of Cervical Cancer

Long non-coding RNAs (lncRNAs) have emerged as auxiliary regulators of gene expression influencing tumor microenvironment, metastasis and radio-resistance in cancer. The presence of lncRNA in extracellular fluids makes them promising diagnostic markers. LncRNAs deploy higher-order structures to facilitate a complex range of functions. Among such structures, G-quadruplexes (G4s) can be detected or targeted by small molecular probes to drive theranostic applications. The in vitro identification of G4 formation in lncRNAs can be a tedious and expensive proposition. Bioinformatics-driven strategies can provide comprehensive and economic alternatives in conjunction with suitable experimental validation. We propose a pipeline to identify G4-forming sequences, protein partners and biological functions associated with dysregulated lncRNAs in cervical cancer. We identified 17 lncRNA clusters which possess transcripts that can fold into a G4 structure. We confirmed in vitro G4 formation in the four biologically active isoforms of SNHG20, MEG3, CRNDE and LINP1 by Circular Dichroism spectroscopy and Thioflavin-T-assisted fluorescence spectroscopy and reverse-transcriptase stop assay. Gene expression data demonstrated that these four lncRNAs can be potential prognostic biomarkers of cervical cancer. Two approaches were employed for identifying G4 specific protein partners for these lncRNAs and FMR2 was a potential interacting partner for all four clusters. We report a detailed investigation of G4 formation in lncRNAs that are dysregulated in cervical cancer. LncRNAs MEG3, CRNDE, LINP1 and SNHG20 are shown to influence cervical cancer progression and we report G4 specific protein partners for these lncRNAs. The protein partners and G4s predicted in lncRNAs can be exploited for theranostic objectives.

Mono-valent salt corrections for RNA secondary structures in the ViennaRNA package

BACKGROUND

RNA features a highly negatively charged phosphate backbone that attracts a cloud of counter-ions that reduce the electrostatic repulsion in a concentration dependent manner. Ion concentrations thus have a large influence on folding and stability of RNA structures. Despite their well-documented effects, salt effects are not handled consistently by currently available secondary structure prediction algorithms. Combining Debye-Hückel potentials for line charges and Manning's counter-ion condensation theory, Einert et al. (Biophys J 100: 2745-2753, 2011) modeled the energetic contributions of monovalent cations on loops and helices.

RESULTS

The model of Einert et al. is adapted to match the structure of the dynamic programming recursion of RNA secondary structure prediction algorithms. An empirical term describing the salt dependence of the duplex initiation energy is added to improve co-folding predictions for two or more RNA strands. The slightly modified model is implemented in the ViennaRNA package in such way that only the energy parameters but not the algorithmic structure is affected. A comparison with data from the literature show that predicted free energies and melting temperatures are in reasonable agreement with experiments.

CONCLUSION

The new feature in the ViennaRNA package makes it possible to study effects of salt concentrations on RNA folding in a systematic manner. Strictly speaking, the model pertains only to mono-valent cations, and thus covers the most important parameter, i.e., the NaCl concentration. It remains a question for future research to what extent unspecific effects of bi- and tri-valent cations can be approximated in a similar manner.

AVAILABILITY

Corrections for the concentration of monovalent cations are available in the ViennaRNA package starting from version 2.6.0.

G-QINDER Tool: Bioinformatically Predicted Formation of Different Four-Stranded DNA Motifs from (GT)n and (GA)n Repeats

The recently introduced semi-orthogonal system of nucleic acid imaging offers a greatly improved method of identifying DNA sequences that are capable of adopting noncanonical structures. This paper uses our newly developed G-QINDER tool to identify specific repeat sequences that adopt unique structural motifs in DNA: TG and AG repeats. The structures were found to adopt a left-handed G-quadruplex form under extreme crowding conditions and a unique tetrahelical motif under certain other conditions. The tetrahelical structure likely consists of stacked AGAG-tetrads but, unlike G-quadruplexes, their stability does not appear to be dependent on the type of monovalent cation present. The occurrence of TG and AG repeats in genomes is not rare, and they are also found frequently in the regulatory regions of nucleic acids, so it is reasonable to assume that putative structural motifs, like other noncanonical forms, could play an important regulatory role in cells. This hypothesis is supported by the structural stability of the AGAG motif; its unfolding can occur even at physiological temperatures since the melting temperature is primarily dependent on the number of AG repeats in the sequence.

Small molecule-based detection of non-canonical RNA G-quadruplex structures that modulate protein translation

Tandem repeats of guanine-rich sequences in RNA often form thermodynamically stable four-stranded RNA structures. Such RNA G-quadruplexes have long been considered to be linked to essential biological processes, yet their physiological significance in cells remains unclear. Here, we report a approach that permits the detection of RNA G-quadruplex structures that modulate protein translation in mammalian cells. The approach combines antibody arrays and RGB-1, a small molecule that selectively stabilizes RNA G-quadruplex structures. Analysis of the protein and mRNA products of 84 cancer-related human genes identified Nectin-4 and CapG as G-quadruplex-controlled genes whose mRNAs harbor non-canonical G-quadruplex structures on their 5′UTR region. Further investigations revealed that the RNA G-quadruplex of CapG exhibits a structural polymorphism, suggesting a possible mechanism that ensures the translation repression in a KCl concentration range of 25–100 mM. The approach described in the present study sets the stage for further discoveries of RNA G-quadruplexes.

A Key Molecular Regulator, RNA G-Quadruplex and Its Function in Plants

RNA structure plays key roles in plant growth, development, and adaptation. One of the complex RNA structures is the RNA G-quadruplex (RG4) where guanine-rich sequences are folded into two or more layers of G-quartets. Previous computational predictions of RG4 revealed that it is widespread across the whole transcriptomes in many plant species, raising the hypothesis that RG4 is likely to be an important regulatory motif in plants. Recently, with the advances in both high-throughput sequencing and cell imaging technologies, RG4 can be detected in living cells as well as at the genome-wide scale. Here, we provide a comprehensive review of recent developments in new methods for detecting RG4 in plants. We also summarize the new functions of RG4 in regulating plant growth and development. We then discuss the possible role of RG4 in adapting to environmental conditions along with evolutionary perspectives.

RNA G-quadruplexes (rG4s): genomics and biological functions

G-quadruplexes (G4s) are non-classical DNA or RNA secondary structures that have been first observed decades ago. Over the years, these four-stranded structural motifs have been demonstrated to have significant regulatory roles in diverse biological processes, but challenges remain in detecting them globally and reliably. Compared to DNA G4s (dG4s), the study of RNA G4s (rG4s) has received less attention until recently. In this review, we will summarize the innovative high-throughput methods recently developed to detect rG4s on a transcriptome-wide scale, highlight the many novel and important functions of rG4 being discovered in vivo across the tree of life, and discuss the key biological questions to be addressed in the near future.

Topology-based classification of tetrads and quadruplex structures

Motivation

Quadruplexes attract the attention of researchers from many fields of bio-science. Due to a specific structure, these tertiary motifs are involved in various biological processes. They are also promising therapeutic targets in many strategies of drug development, including anticancer and neurological disease treatment. The uniqueness and diversity of their forms cause that quadruplexes show great potential in novel biological applications. The existing approaches for quadruplex analysis are based on sequence or 3D structure features and address canonical motifs only.

Results

In our study, we analyzed tetrads and quadruplexes contained in nucleic acid molecules deposited in Protein Data Bank. Focusing on their secondary structure topology, we adjusted its graphical diagram and proposed new dot-bracket and arc representations. We defined the novel classification of these motifs. It can handle both canonical and non-canonical cases. Based on this new taxonomy, we implemented a method that automatically recognizes the types of tetrads and quadruplexes occurring as unimolecular structures. Finally, we conducted a statistical analysis of these motifs found in experimentally determined nucleic acid structures in relation to the new classification.

Availability and implementation

https://github.com/tzok/eltetrado/

Supplementary information

Supplementary data are available at Bioinformatics online.

Fundamentals of G-quadruplex biology

Several decades elapsed between the first descriptions of G-quadruplex nucleic acid structures (G4s) assembled in vitro and the emergence of experimental findings indicating that such structures can form and function in living systems. A large body of evidence now supports roles for G4s in many aspects of nucleic acid biology, spanning processes from transcription and chromatin structure, mRNA processing, protein translation, DNA replication and genome stability, and telomere and mitochondrial function. Nonetheless, it must be acknowledged that some of this evidence is tentative, which is not surprising given the technical challenges associated with demonstrating G4s in biology. Here I provide an overview of evidence for G4 biology, focusing particularly on the many potential pitfalls that can be encountered in its investigation, and briefly discuss some of broader biological processes that may be impacted by G4s including cancer.

G-quadruplexes in mRNA: A key structure for biological function

The last several years have seen exciting advances in the understanding of the structure and function of higher-order structures of RNA. Expression levels of some specific genes were shown to be directly regulated by environmentally-responsive formation of certain secondary structures such as stem-loops and pseudoknots. Even among these noncanonical structures, RNA G-quadruplexes, which form on the regions of guanine-rich sequences in mRNA, are highly stable structures that are involved in a variety of biological processes. However, many questions regarding the biological significance of RNA G-quadruplexes remain unsettled, mainly because it is difficult to locate the structures in mRNA. This review focuses on emerging methods that locate RNA G-quadruplexes in mRNA by computational and biochemical techniques. In addition, recent reports on the biological functions of RNA G-quadruplexes are also covered to highlight their various roles in cells, such as in regulating mRNA processing and translation.

A guide to computational methods for G-quadruplex prediction

Guanine-rich nucleic acids can fold into the non-B DNA or RNA structures called G-quadruplexes (G4). Recent methodological developments have allowed the characterization of specific G-quadruplex structures in vitro as well as in vivo, and at a much higher throughput, in silico, which has greatly expanded our understanding of G4-associated functions. Typically, the consensus motif G3+N1-7G3+N1-7G3+N1-7G3+ has been used to identify potential G-quadruplexes from primary sequence. Since, various algorithms have been developed to predict the potential formation of quadruplexes directly from DNA or RNA sequences and the number of studies reporting genome-wide G4 exploration across species has rapidly increased. More recently, new methodologies have also appeared, proposing other estimates which consider non-canonical sequences and/or structure propensity and stability. The present review aims at providing an updated overview of the current open-source G-quadruplex prediction algorithms and straightforward examples of their implementation.

Stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells

This review provides the biophysicochemical background and recent advances in stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells.

Stable G-quadruplex enabling sequences are selected against by the context-dependent codon bias

The distributions of secondary structural elements appear to differ between coding regions (CDS) of mRNAs compared to the untranslated regions (UTRs), presumably as a mechanism to fine-tune gene expression, including efficiency of translation. However, a systematic and comprehensive analysis of secondary structure avoidance because of potential bias in codon usage is difficult as some of the common secondary structures, such as, hairpins can be formed by numerous sequence combinations. Using G-quadruplex (GQ) as the model secondary structure we studied the impact of codon bias on GQs within the CDS. Because GQs can be predicted using specific consensus sequence motifs, they provide an excellent platform for investigation of the selectivity of such putative structures at the codon level. Using a bioinformatics approach, we calculated the frequencies of putative GQs within the CDS of a variety of species. Our results suggest that the most stable GQs appear to be significantly underrepresented within the CDS, through the use of specific synonymous codon combinations. Furthermore, we identified many peptide sequence motifs in which silent mutations can potentially alter translation via stable GQ formation. This work not only provides a comprehensive analysis on how stable secondary structures appear to be avoided within the CDS of mRNA, but also broadens the current understanding of synonymous codon usage as they relate to the structure-function relationship of RNA.

G-quadruplexes formation in the 5'UTRs of mRNAs associated with colorectal cancer pathways

RNA G-quadruplexes (rG4) are stable non-canonical secondary structures composed of G-rich sequences. Many rG4 structures located in the 5'UTRs of mRNAs act as translation repressors due to their high stability which is thought to impede ribosomal scanning. That said, it is not known if these are mRNA-specific examples, or if they are indicative of a global expression regulation mechanism of the mRNAs involved in a common pathway based on structure folding recognition. Gene-ontology analysis of mRNAs bearing a predicted rG4 motif in their 5'UTRs revealed an enrichment for mRNAs associated with the colorectal cancer pathway. Bioinformatic tools for rG4 prediction, and experimental in vitro validations were used to confirm and compare the folding of the predicted rG4s of the mRNAs associated with dysregulated pathways in colorectal cancer. The rG4 folding was confirmed for the first time for 9 mRNAs. A repressive effect of 3 rG4 candidates on the expression of a reporter gene was also measured in colorectal cancer cell lines. This work highlights the fact that rG4 prediction is not yet accurate, and that experimental characterization is still essential in order to identify the precise rG4 folding sequences and the possible common features shared between the rG4 overrepresented in important biological pathways.

Motif independent identification of potential RNA G-quadruplexes by G4RNA screener

Motivation

G-quadruplex structures in RNA molecules are known to have regulatory impacts in cells but are difficult to locate in the genome. The minimal requirements for G-quadruplex folding in RNA (G_≥3N_1-7G_≥3N_1-7G_≥3N_1-7G_≥3) is being challenged by observations made on specific examples in recent years. The definition of potential G-quadruplex sequences has major repercussions on the observation of the structure since it introduces a bias. The canonical motif only describes a sub-population of the reported G-quadruplexes. To address these issues, we propose an RNA G-quadruplex prediction strategy that does not rely on a motif definition.

Results

We trained an artificial neural network with sequences of experimentally validated G-quadruplexes from the G4RNA database encoded using an abstract definition of their sequence. This artificial neural network, G4NN, evaluates the similarity of a given sequence to known G-quadruplexes and reports it as a score. G4NN has a predictive power comparable to the reported G richness and G/C skewness evaluations that are the current state-of-the-art for the identification of potential RNA G-quadruplexes. We combined these approaches in the G4RNA screener, a program designed to manage and evaluate the sequences to identify potential G-quadruplexes.

Availability and implementation

G4RNA screener is available for download at http://gitlabscottgroup.med.usherbrooke.ca/J-Michel/g4rna_screener.

Supplementary information

Supplementary data are available at Bioinformatics online.

Detecting RNA G-Quadruplexes (rG4s) in the Transcriptome

RNA G-quadruplex (rG4) secondary structures are proposed to play key roles in fundamental biological processes that include the modulation of transcriptional, co-transcriptional, and posttranscriptional events. Recent methodological developments that include predictive algorithms and structure-based sequencing have enabled the detection and mapping of rG4 structures on a transcriptome-wide scale at high sensitivity and resolution. The data generated by these studies provide valuable insights into the potentially diverse roles of rG4s in biology and open up a number of mechanistic hypotheses. Herein we highlight these methodologies and discuss the associated findings in relation to rG4-related biological mechanisms.

G4RNA screener web server: User focused interface for RNA G-quadruplex prediction

Though RNA G-quadruplexes became a focus of study over a decade ago, the main challenge associated with the identification of new potential G-quadruplexes remains a bottleneck step. It slows the study of these non-canonical structures in nucleic acids, and thus the understanding of their significance. The G4RNA screener is an accurate tool for the prediction of RNA G-quadruplexes but its deployment has brought to light an issue with its accessibility to G-quadruplex experts and biologists. G4RNA screener web server is a platform that provides a much needed interface to manage the input, parameters and result display of the main command-line ready tool. It is accessible at http://scottgroup.med.usherbrooke.ca/G4RNA_screener/.

Single-Stranded Condensation Stochastically Blocks G-Quadruplex Assembly in Human Telomeric RNA

TERRA is an RNA molecule transcribed from human subtelomeric regions toward chromosome ends potentially involved in regulation of heterochromatin stability, semiconservative replication, and telomerase inhibition, among others. TERRA contains tandem repeats of the sequence GGGUUA, with a strong tendency to fold into a four-stranded arrangement known as a parallel G-quadruplex. Here, we demonstrate by using single-molecule force spectroscopy that this potential is limited by the inherent capacity of RNA to self-associate randomly and further condense into entropically more favorable structures. We stretched RNA constructions with more than four and less than eight hexanucleotide repeats, thus unable to form several G-quadruplexes in tandem, flanked by non-G-rich overhangs of random sequence by optical tweezers on a one by one basis. We found that condensed RNA stochastically blocks G-quadruplex folding pathways with a near 20% probability, a behavior that is not found in DNA analogous molecules.

Machine learning model for sequence-driven DNA G-quadruplex formation

We describe a sequence-based computational model to predict DNA G-quadruplex (G4) formation. The model was developed using large-scale machine learning from an extensive experimental G4-formation dataset, recently obtained for the human genome via G4-seq methodology. Our model differentiates many widely accepted putative quadruplex sequences that do not actually form stable genomic G4 structures, correctly assessing the G4 folding potential of over 700,000 such sequences in the human genome. Moreover, our approach reveals the relative importance of sequence-based features coming from both within the G4 motifs and their flanking regions. The developed model can be applied to any DNA sequence or genome to characterise sequence-driven intramolecular G4 formation propensities.

RNA structure prediction: from 2D to 3D

We summarize different levels of RNA structure prediction, from classical 2D structure to extended secondary structure and motif-based research toward 3D structure prediction of RNA. We outline the importance of classical secondary structure during all those levels of structure prediction.

A new cis-acting motif is required for the axonal SMN-dependent Anxa2 mRNA localization

Spinal muscular atrophy (SMA) is caused by mutations and/or deletions of the survival motor neuron gene (SMN1). Besides its function in the biogenesis of spliceosomal snRNPs, SMN might possess a motor neuron specific role and could function in the transport of axonal mRNAs and in the modulation of local protein translation. Accordingly, SMN colocalizes with axonal mRNAs of differentiated NSC-34 motor neuron-like cells. We recently showed that SMN depletion gives rise to a decrease in the axonal transport of the mRNAs encoding Annexin A2 (Anxa2). In this work, we have characterized the structural features of the Anxa2 mRNA required for its axonal targeting by SMN. We found that a G-rich motif located near the 3'UTR is essential for axonal localization of the Anxa2 transcript. We also show that mutations in the motif sequence abolish targeting of Anxa2 reporter mRNAs in axon-like structures of differentiated NSC-34 cells. Finally, localization of both wild-type and mutated Anxa2 reporters is restricted to the cell body in SMN-depleted cells. Altogether, our studies show that this G-motif represents a novel and essential determinant for axonal localization of the Anxa2 mRNA mediated by the SMN complex.

Cancer-associated noncoding mutations affect RNA G-quadruplex-mediated regulation of gene expression

Cancer is a multifactorial disease driven by a combination of genetic and environmental factors. Many cancer driver mutations have been characterised in protein-coding regions of the genome. However, mutations in noncoding regions associated with cancer have been less investigated. G-quadruplex (G4) nucleic acids are four-stranded secondary structures formed in guanine-rich sequences and prevalent in the regulatory regions. In this study, we used published whole cancer genome sequence data to find mutations in cancer patients that overlap potential RNA G4-forming sequences in 5′ UTRs. Using RNAfold, we assessed the effect of these mutations on the thermodynamic stability of predicted RNA G4s in the context of full-length 5′ UTRs. Of the 217 identified mutations, we found that 33 are predicted to destabilise and 21 predicted to stabilise potential RNA G4s. We experimentally validated the effect of destabilising mutations in the 5′ UTRs of BCL2 and CXCL14 and one stabilising mutation in the 5′ UTR of TAOK2. These mutations resulted in an increase or a decrease in translation of these mRNAs, respectively. These findings suggest that mutations that modulate the G4 stability in the noncoding regions could act as cancer driver mutations, which present an opportunity for early cancer diagnosis using individual sequencing information.

RNA G-Quadruplexes as Key Motifs of the Transcriptome

G-Quadruplexes are non-canonical secondary structures that can be adopted under physiological conditions by guanine-rich DNA and RNA molecules. They have been reported to occur, and to perform multiple biological functions, in the genomes and transcriptomes of many species, including humans. This chapter focuses specifically on RNA G-quadruplexes and reviews the most recent discoveries in the field, as well as addresses the upcoming challenges researchers studying these structures face.

Multiple G-quartet structures in pre-edited mRNAs suggest evolutionary driving force for RNA editing in trypanosomes

Mitochondrial transcript maturation in African trypanosomes requires a U-nucleotide specific RNA editing reaction. In its most extreme form hundreds of U's are inserted into and deleted from primary transcripts to generate functional mRNAs. Unfortunately, both origin and biological role of the process have remained enigmatic. Here we report a so far unrecognized structural feature of pre-edited mRNAs. We demonstrate that the cryptic pre-mRNAs contain numerous clustered G-nt, which fold into G-quadruplex (GQ) structures. We identified 27 GQ's in the different pre-mRNAs and demonstrate a positive correlation between the steady state abundance of guide (g)RNAs and the sequence position of GQ-elements. We postulate that the driving force for selecting G-rich sequences lies in the formation of DNA/RNA hybrid G-quadruplex (HQ) structures between the pre-edited transcripts and the non-template strands of mitochondrial DNA. HQ's are transcription termination/replication initiation sites and thus guarantee an unperturbed replication of the mt-genome. This is of special importance in the insect-stage of the parasite. In the transcription-on state, the identified GQ's require editing as a GQ-resolving activity indicating a link between replication, transcription and RNA editing. We propose that the different processes have coevolved and suggest the parasite life-cycle and the single mitochondrion as evolutionary driving forces.

RNA folding with hard and soft constraints

Background

A large class of RNA secondary structure prediction programs uses an elaborate energy model grounded in extensive thermodynamic measurements and exact dynamic programming algorithms. External experimental evidence can be in principle be incorporated by means of hard constraints that restrict the search space or by means of soft constraints that distort the energy model. In particular recent advances in coupling chemical and enzymatic probing with sequencing techniques but also comparative approaches provide an increasing amount of experimental data to be combined with secondary structure prediction.

Results

Responding to the increasing needs for a versatile and user-friendly inclusion of external evidence into diverse flavors of RNA secondary structure prediction tools we implemented a generic layer of constraint handling into the ViennaRNA Package. It makes explicit use of the conceptual separation of the “folding grammar” defining the search space and the actual energy evaluation, which allows constraints to be interleaved in a natural way between recursion steps and evaluation of the standard energy function.

Conclusions

The extension of the ViennaRNA Package provides a generic way to include diverse types of constraints into RNA folding algorithms. The computational overhead incurred is negligible in practice. A wide variety of application scenarios can be accommodated by the new framework, including the incorporation of structure probing data, non-standard base pairs and chemical modifications, as well as structure-dependent ligand binding.

Electronic supplementary material

The online version of this article (doi:10.1186/s13015-016-0070-z) contains supplementary material, which is available to authorized users.

Predicting RNA secondary structures from sequence and probing data

RNA secondary structures have proven essential for understanding the regulatory functions performed by RNA such as microRNAs, bacterial small RNAs, or riboswitches. This success is in part due to the availability of efficient computational methods for predicting RNA secondary structures. Recent advances focus on dealing with the inherent uncertainty of prediction by considering the ensemble of possible structures rather than the single most stable one. Moreover, the advent of high-throughput structural probing has spurred the development of computational methods that incorporate such experimental data as auxiliary information.

The RNA chaperone activity of the Trypanosoma brucei editosome raises the dynamic of bound pre-mRNAs

Mitochondrial transcript maturation in African trypanosomes requires an RNA editing reaction that is characterized by the insertion and deletion of U-nucleotides into otherwise non-functional mRNAs. The reaction is catalyzed by editosomes and requires guide (g)RNAs as templates. Recent data demonstrate that the binding of pre-edited mRNAs to editosomes is followed by a chaperone-type RNA remodeling reaction. Here we map the changes in RNA folding using selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE). We demonstrate that pre-mRNAs in their free state adopt intricately folded, highly stable 2D-structures. Editosome binding renders the pre-mRNAs to adopt 2D-conformations of reduced stabilities. On average about 30% of the nucleotides in every pre-mRNA are affected with a prevalence for U-nucleotides. The data demonstrate that the chaperone activity acts by increasing the flexibility of U-residues to lower their base-pairing probability. This results in a simplified RNA folding landscape with a reduced energy barrier to facilitate the binding of gRNAs. The data provide a first rational for the enigmatic U-specificity of the editing reaction.

Aven recognition of RNA G-quadruplexes regulates translation of the mixed lineage leukemia protooncogenes

G-quadruplexes (G4) are extremely stable secondary structures forming stacks of guanine tetrads. DNA G4 structures have been extensively studied, however, less is known about G4 motifs in mRNAs, especially in their coding sequences. Herein, we show that Aven stimulates the mRNA translation of the mixed lineage leukemia (MLL) proto-oncogene in an arginine methylation-dependent manner. The Aven RGG/RG motif bound G4 structures within the coding regions of the MLL1 and MLL4 mRNAs increasing their polysomal association and translation, resulting in the induction of transcription of leukemic genes. The DHX36 RNA helicase associated with the Aven complex and was required for optimal translation of G4 mRNAs. Depletion of Aven led to a decrease in synthesis of MLL1 and MLL4 proteins resulting in reduced proliferation of leukemic cells. These findings identify an Aven-centered complex that stimulates the translation of G4 harboring mRNAs, thereby promoting survival of leukemic cells.

DOI:http://dx.doi.org/10.7554/eLife.06234.001

To make a protein, the DNA sequence that encodes it is first copied to make a molecule of messenger RNA (or mRNA for short). The mRNA is then used as a set of instructions to assemble a protein in a process called translation. Both DNA and RNA molecules can fold into particular shapes. One such structure is known as a G-quartet and involves the DNA or RNA folding back on itself to form a highly stable planar structure. Stacks of G-quartets can form structures known as G-quadruplexes, but little is known about the G-quadruplexes that form in mRNA molecules.

Leukemia affects cells in the bone marrow and causes blood cells to develop abnormally. A protein called Aven is often found in increased amounts in certain types of leukemic cells, but it was not clear how Aven affects how leukemia develops. Thandapani, Song et al. have now found that in leukemic cells, Aven binds to G-quadruplexes found in two mRNA molecules that encode proteins that are linked to leukemia. This binding increases the translation of these mRNAs, with translation occurring most efficiently when a particular enzyme called a helicase—which remodels RNA—also bound to Aven. Reducing the amount of Aven in cells caused fewer of the leukemic proteins to be produced, which also reduced the growth and multiplcation of leukemic cells.

These findings raise the possibility that drugs that disrupt how Aven works could form part of treatments for leukemia. The next challenge will be to identify the signaling pathways that communicate with Aven and to define all the G-quadruplex mRNAs that are regulated by Aven.

DOI:http://dx.doi.org/10.7554/eLife.06234.002

G4RNA: an RNA G-quadruplex database

G-quadruplexes (G4) are tetrahelical structures formed from planar arrangement of guanines in nucleic acids. A simple, regular motif was originally proposed to describe G4-forming sequences. More recently, however, formation of G4 was discovered to depend, at least in part, on the contextual backdrop of neighboring sequences. Prediction of G4 folding is thus becoming more challenging as G4 outlier structures, not described by the originally proposed motif, are increasingly reported. Recent observations thus call for a comprehensive tool, capable of consolidating the expanding information on tested G4s, in order to conduct systematic comparative analyses of G4-promoting sequences. The G4RNA Database we propose was designed to help meet the need for easily-retrievable data on known RNA G4s. A user-friendly, flexible query system allows for data retrieval on experimentally tested sequences, from many separate genes, to assess G4-folding potential. Query output sorts data according to sequence position, G4 likelihood, experimental outcomes and associated bibliographical references. G4RNA also provides an ideal foundation to collect and store additional sequence and experimental data, considering the growing interest G4s currently generate.

Database URL:scottgroup.med.usherbrooke.ca/G4RNA

Graph-distance distribution of the Boltzmann ensemble of RNA secondary structures

Background

Large RNA molecules are often composed of multiple functional domains whose spatial arrangement strongly influences their function. Pre-mRNA splicing, for instance, relies on the spatial proximity of the splice junctions that can be separated by very long introns. Similar effects appear in the processing of RNA virus genomes. Albeit a crude measure, the distribution of spatial distances in thermodynamic equilibrium harbors useful information on the shape of the molecule that in turn can give insights into the interplay of its functional domains.

Result

Spatial distance can be approximated by the graph-distance in RNA secondary structure. We show here that the equilibrium distribution of graph-distances between a fixed pair of nucleotides can be computed in polynomial time by means of dynamic programming. While a naïve implementation would yield recursions with a very high time complexity of O(n⁶D⁵) for sequence length n and D distinct distance values, it is possible to reduce this to O(n⁴) for practical applications in which predominantly small distances are of of interest. Further reductions, however, seem to be difficult. Therefore, we introduced sampling approaches that are much easier to implement. They are also theoretically favorable for several real-life applications, in particular since these primarily concern long-range interactions in very large RNA molecules.

Conclusions

The graph-distance distribution can be computed using a dynamic programming approach. Although a crude approximation of reality, our initial results indicate that the graph-distance can be related to the smFRET data. The additional file and the software of our paper are available from http://www.rna.uni-jena.de/RNAgraphdist.html.

Optimal antisense target reducing INS intron 1 retention is adjacent to a parallel G quadruplex

Splice-switching oligonucleotides (SSOs) have been widely used to inhibit exon usage but antisense strategies that promote removal of entire introns to increase splicing-mediated gene expression have not been developed. Here we show reduction of INS intron 1 retention by SSOs that bind transcripts derived from a human haplotype expressing low levels of proinsulin. This haplotype is tagged by a polypyrimidine tract variant rs689 that decreases the efficiency of intron 1 splicing and increases the relative abundance of mRNAs with extended 5' untranslated region (5' UTR), which curtails translation. Co-expression of haplotype-specific reporter constructs with SSOs bound to splicing regulatory motifs and decoy splice sites in primary transcripts revealed a motif that significantly reduced intron 1-containing mRNAs. Using an antisense microwalk at a single nucleotide resolution, the optimal target was mapped to a splicing silencer containing two pseudoacceptor sites sandwiched between predicted RNA guanine (G) quadruplex structures. Circular dichroism spectroscopy and nuclear magnetic resonance of synthetic G-rich oligoribonucleotide tracts derived from this region showed formation of a stable parallel 2-quartet G-quadruplex on the 3' side of the antisense retention target and an equilibrium between quadruplexes and stable hairpin-loop structures bound by optimal SSOs. This region interacts with heterogeneous nuclear ribonucleoproteins F and H that may interfere with conformational transitions involving the antisense target. The SSO-assisted promotion of weak intron removal from the 5' UTR through competing noncanonical and canonical RNA structures may facilitate development of novel strategies to enhance gene expression.

Automated identification of RNA 3D modules with discriminative power in RNA structural alignments

Recent progress in predicting RNA structure is moving towards filling the ‘gap’ in 2D RNA structure prediction where, for example, predicted internal loops often form non-canonical base pairs. This is increasingly recognized with the steady increase of known RNA 3D modules. There is a general interest in matching structural modules known from one molecule to other molecules for which the 3D structure is not known yet. We have created a pipeline, metaRNAmodules, which completely automates extracting putative modules from the FR3D database and mapping of such modules to Rfam alignments to obtain comparative evidence. Subsequently, the modules, initially represented by a graph, are turned into models for the RMDetect program, which allows to test their discriminative power using real and randomized Rfam alignments. An initial extraction of 22 495 3D modules in all PDB files results in 977 internal loop and 17 hairpin modules with clear discriminatory power. Many of these modules describe only minor variants of each other. Indeed, mapping of the modules onto Rfam families results in 35 unique locations in 11 different families. The metaRNAmodules pipeline source for the internal loop modules is available at http://rth.dk/resources/mrm.