Specificity of mRNA Folding and Its Association with Evolutionarily Adaptive mRNA Secondary Structures

The secondary structure is a fundamental feature of both non-coding RNAs (ncRNAs) and messenger RNAs (mRNAs). However, our understanding of the secondary structures of mRNAs, especially those of the coding regions, remains elusive, likely due to translation and the lack of RNA-binding proteins that sustain the consensus structure like those binding to ncRNAs. Indeed, mRNAs have recently been found to adopt diverse alternative structures, but the overall functional significance remains untested. We hereby approach this problem by estimating the folding specificity, i.e., the probability that a fragment of an mRNA folds back to the same partner once refolded. We show that the folding specificity of mRNAs is lower than that of ncRNAs and exhibits moderate evolutionary conservation. Notably, we find that specific rather than alternative folding is likely evolutionarily adaptive since specific folding is frequently associated with functionally important genes or sites within a gene. Additional analysis in combination with ribosome density suggests the ability to modulate ribosome movement as one potential functional advantage provided by specific folding. Our findings reveal a novel facet of the RNA structurome with important functional and evolutionary implications and indicate a potential method for distinguishing the mRNA secondary structures maintained by natural selection from molecular noise.     

Local absence of secondary structure permits translation of mRNAs that lack ribosome-binding sites

The initiation of translation is a fundamental and highly regulated process in gene expression. Translation initiation in prokaryotic systems usually requires interaction between the ribosome and an mRNA sequence upstream of the initiation codon, the so-called ribosome-binding site (Shine-Dalgarno sequence). However, a large number of genes do not possess Shine-Dalgarno sequences, and it is unknown how start codon recognition occurs in these mRNAs. We have performed genome-wide searches in various groups of prokaryotes in order to identify sequence elements and/or RNA secondary structural motifs that could mediate translation initiation in mRNAs lacking Shine-Dalgarno sequences. We find that mRNAs without a Shine-Dalgarno sequence are generally less structured in their translation initiation region and show a minimum of mRNA folding at the start codon. Using reporter gene constructs in bacteria, we also provide experimental support for local RNA unfoldedness determining start codon recognition in Shine-Dalgarno--independent translation. Consistent with this, we show that AUG start codons reside in single-stranded regions, whereas internal AUG codons are usually in structured regions of the mRNA. Taken together, our bioinformatics analyses and experimental data suggest that local absence of RNA secondary structure is necessary and sufficient to initiate Shine-Dalgarno--independent translation. Thus, our results provide a plausible mechanism for how the correct translation initiation site is recognized in the absence of a ribosome-binding site.     

Deciphering the rules by which dynamics of mRNA secondary structure affect translation efficiency in Saccharomyces cerevisiae

Messenger RNA (mRNA) secondary structure decreases the elongation rate, as ribosomes must unwind every structure they encounter during translation. Therefore, the strength of mRNA secondary structure is assumed to be reduced in highly translated mRNAs. However, previous studies in vitro reported a positive correlation between mRNA folding strength and protein abundance. The counterintuitive finding suggests that mRNA secondary structure affects translation efficiency in an undetermined manner. Here, we analyzed the folding behavior of mRNA during translation and its effect on translation efficiency. We simulated translation process based on a novel computational model, taking into account the interactions among ribosomes, codon usage and mRNA secondary structures. We showed that mRNA secondary structure shortens ribosomal distance through the dynamics of folding strength. Notably, when adjacent ribosomes are close, mRNA secondary structures between them disappear, and codon usage determines the elongation rate. More importantly, our results showed that the combined effect of mRNA secondary structure and codon usage in highly translated mRNAs causes a short ribosomal distance in structural regions, which in turn eliminates the structures during translation, leading to a high elongation rate. Together, these findings reveal how the dynamics of mRNA secondary structure coupling with codon usage affect translation efficiency.     

The separation between the 5′-3′ ends in long RNA molecules is short and nearly constant

RNA molecules play different roles in coding, decoding and gene expression regulation. Such roles are often associated to the RNA secondary or tertiary structures. The folding dynamics lead to multiple secondary structures of long RNA molecules, since an RNA molecule might fold into multiple distinct native states. Despite an ensemble of different structures, it has been theoretically proposed that the separation between the 5′ and 3′ ends of long single-stranded RNA molecules (ssRNA) remains constant, independent of their base content and length. Here, we present the first experimental measurements of the end-to-end separation in long ssRNA molecules. To determine this separation, we use single molecule Fluorescence Resonance Energy Transfer of fluorescently end-labeled ssRNA molecules ranging from 500 to 5500 nucleotides in length, obtained from two viruses and a fungus. We found that the end-to-end separation is indeed short, within 5–9 nm. It is remarkable that the separation of the ends of all RNA molecules studied remains small and similar, despite the origin, length and differences in their secondary structure. This implies that the ssRNA molecules are ‘effectively circularized’ something that might be a general feature of RNAs, and could result in fine-tuning for translation and gene expression regulation.     

From structure prediction to genomic screens for novel non-coding RNAs

Non-coding RNAs (ncRNAs) are receiving more and more attention not only as an abundant class of genes, but also as regulatory structural elements (some located in mRNAs). A key feature of RNA function is its structure. Computational methods were developed early for folding and prediction of RNA structure with the aim of assisting in functional analysis. With the discovery of more and more ncRNAs, it has become clear that a large fraction of these are highly structured. Interestingly, a large part of the structure is comprised of regular Watson-Crick and GU wobble base pairs. This and the increased amount of available genomes have made it possible to employ structure-based methods for genomic screens. The field has moved from folding prediction of single sequences to computational screens for ncRNAs in genomic sequence using the RNA structure as the main characteristic feature. Whereas early methods focused on energy-directed folding of single sequences, comparative analysis based on structure preserving changes of base pairs has been efficient in improving accuracy, and today this constitutes a key component in genomic screens. Here, we cover the basic principles of RNA folding and touch upon some of the concepts in current methods that have been applied in genomic screens for de novo RNA structures in searches for novel ncRNA genes and regulatory RNA structure on mRNAs. We discuss the strengths and weaknesses of the different strategies and how they can complement each other.     

Genome-wide analysis demonstrates conserved localization of messenger RNAs to mitotic microtubules

RNA localization is of critical importance in many fundamental cell biological and developmental processes by regulating the spatial control of gene expression. To investigate how spindle-localized RNAs might influence mitosis, we comprehensively surveyed all messenger RNAs (mRNAs) that bound to microtubules during metaphase in both Xenopus laevis egg extracts and mitotic human cell extracts. We identify conserved classes of mRNAs that are enriched on microtubules in both human and X. laevis. Active mitotic translation occurs on X. laevis meiotic spindles, and a subset of microtubule-bound mRNAs (MT-mRNAs) associate with polyribosomes. Although many MT-mRNAs associate with polyribosomes, we find that active translation is not required for mRNA localization to mitotic microtubules. Our results represent the first genome-wide survey of mRNAs localized to a specific cytoskeletal component and suggest that microtubule localization of specific mRNAs is likely to function in mitotic regulation and mRNA segregation during cell division.     

Major role for mRNA binding and restructuring in sRNA recruitment by Hfq

Bacterial small RNAs (sRNAs) modulate gene expression by base-pairing with target mRNAs. Many sRNAs require the Sm-like RNA binding protein Hfq as a cofactor. Well-characterized interactions between DsrA sRNA and the rpoS mRNA leader were used to understand how Hfq stimulates sRNA pairing with target mRNAs. DsrA annealing stimulates expression of rpoS by disrupting a secondary structure in the rpoS leader, which otherwise prevents translation. Both RNAs bind Hfq with similar affinity but interact with opposite faces of the Hfq hexamer. Using mutations that block interactions between two of the three components, we demonstrate that Hfq binding to a functionally critical (AAN)(4) motif in rpoS mRNA rescues DsrA binding to a hyperstable rpoS mutant. We also show that Hfq cannot stably bridge the RNAs. Persistent ternary complexes only form when the two RNAs are complementary. Thus, Hfq mainly acts by binding and restructuring the rpoS mRNA. However, Hfq binding to DsrA is needed for maximum annealing in vitro, indicating that transient interactions with both RNAs contribute to the regulatory mechanism.     

Small noncoding RNAs in the germline

Small noncoding RNAs have emerged as potent regulators of gene expression, especially in the germline. We review the biogenesis and regulatory function of three major small noncoding RNA pathways in the germline: The small interfering RNA (siRNA) pathway that leads to the degradation of target mRNAs, the microRNA (miRNA) pathway that mostly represses the translation of target mRNAs, and the newly discovered Piwi-interacting RNA (piRNA) pathway that appears to have diverse functions in epigenetic programming, transposon silencing, and the regulation of mRNA translation and stability. The siRNA and miRNA pathways are present in the germline as well as many somatic tissues, whereas the piRNA pathway is predominantly confined to the germline. Investigation of the three small RNA pathways has started to reveal a new dimension of gene regulation with defining roles in germline specification and development.     

Crystal structure of the lysine riboswitch regulatory mRNA element

Riboswitches are metabolite-sensitive elements found in mRNAs that control gene expression through a regulatory secondary structural switch. Along with regulation of lysine biosynthetic genes, mutations within the lysine-responsive riboswitch (L-box) play a role in the acquisition of resistance to antimicrobial lysine analogs. To understand the structural basis for lysine binding, we have determined the 2.8 angstroms resolution crystal structure of lysine bound to the Thermotoga maritima asd lysine riboswitch ligand-binding domain. The structure reveals a complex architecture scaffolding a binding pocket completely enveloping lysine. Mutations conferring antimicrobial resistance cluster around this site as well as highly conserved long range interactions, indicating that they disrupt lysine binding or proper folding of the RNA. Comparison of the free and bound forms by x-ray crystallography, small angle x-ray scattering, and chemical probing reveals almost identical structures, indicating that lysine induces only limited and local conformational changes upon binding.