The RNA Base-Pairing Problem and Base-Pairing Solutions

RNA molecules are folded into structures and complexes to perform a wide variety of functions. Determination of RNA structures and their interactions is a fundamental problem in RNA biology. Most RNA molecules in living cells are large and dynamic, posing unique challenges to structure analysis. Here we review progress in RNA structure analysis, focusing on methods that use the "cross-link, proximally ligate, and sequence" principle for high-throughput detection of base-pairing interactions in living cells. Beginning with a comparison of commonly used methods in structure determination and a brief historical account of psoralen cross-linking studies, we highlight the important features of cross-linking methods and new biological insights into RNA structures and interactions from recent studies. Further improvement of these cross-linking methods and application to previously intractable problems will shed new light on the mechanisms of the "modern RNA world."

Functional 5' UTR mRNA structures in eukaryotic translation regulation and how to find them

RNA molecules can fold into intricate shapes that can provide an additional layer of control of gene expression beyond that of their sequence. In this Review, we discuss the current mechanistic understanding of structures in 5' untranslated regions (UTRs) of eukaryotic mRNAs and the emerging methodologies used to explore them. These structures may regulate cap-dependent translation initiation through helicase-mediated remodelling of RNA structures and higher-order RNA interactions, as well as cap-independent translation initiation through internal ribosome entry sites (IRESs), mRNA modifications and other specialized translation pathways. We discuss known 5' UTR RNA structures and how new structure probing technologies coupled with prospective validation, particularly compensatory mutagenesis, are likely to identify classes of structured RNA elements that shape post-transcriptional control of gene expression and the development of multicellular organisms.

RNA structure inference through chemical mapping after accidental or intentional mutations

Despite the critical roles RNA structures play in regulating gene expression, sequencing-based methods for experimentally determining RNA base pairs have remained inaccurate. Here, we describe a multidimensional chemical-mapping method called "mutate-and-map read out through next-generation sequencing" (M2-seq) that takes advantage of sparsely mutated nucleotides to induce structural perturbations at partner nucleotides and then detects these events through dimethyl sulfate (DMS) probing and mutational profiling. In special cases, fortuitous errors introduced during DNA template preparation and RNA transcription are sufficient to give M2-seq helix signatures; these signals were previously overlooked or mistaken for correlated double-DMS events. When mutations are enhanced through error-prone PCR, in vitro M2-seq experimentally resolves 33 of 68 helices in diverse structured RNAs including ribozyme domains, riboswitch aptamers, and viral RNA domains with a single false positive. These inferences do not require energy minimization algorithms and can be made by either direct visual inspection or by a neural-network-inspired algorithm called M2-net. Measurements on the P4-P6 domain of the Tetrahymena group I ribozyme embedded in Xenopus egg extract demonstrate the ability of M2-seq to detect RNA helices in a complex biological environment.

RNA structure through multidimensional chemical mapping

The discoveries of myriad non-coding RNA molecules, each transiting through multiple flexible states in cells or virions, present major challenges for structure determination. Advances in high-throughput chemical mapping give new routes for characterizing entire transcriptomes in vivo, but the resulting one-dimensional data generally remain too information-poor to allow accurate de novo structure determination. Multidimensional chemical mapping (MCM) methods seek to address this challenge. Mutate-and-map (M2), RNA interaction groups by mutational profiling (RING-MaP and MaP-2D analysis) and multiplexed •OH cleavage analysis (MOHCA) measure how the chemical reactivities of every nucleotide in an RNA molecule change in response to modifications at every other nucleotide. A growing body of in vitro blind tests and compensatory mutation/rescue experiments indicate that MCM methods give consistently accurate secondary structures and global tertiary structures for ribozymes, ribosomal domains and ligand-bound riboswitch aptamers up to 200 nucleotides in length. Importantly, MCM analyses provide detailed information on structurally heterogeneous RNA states, such as ligand-free riboswitches that are functionally important but difficult to resolve with other approaches. The sequencing requirements of currently available MCM protocols scale at least quadratically with RNA length, precluding general application to transcriptomes or viral genomes at present. We propose a modify-cross-link-map (MXM) expansion to overcome this and other current limitations to resolving the in vivo 'RNA structurome'.

RNAComposer and RNA 3D structure prediction for nanotechnology

RNAs adopt specific, stable tertiary architectures to perform their activities. Knowledge of RNA tertiary structure is fundamental to understand RNA functions beginning with transcription and ending with turnover. Contrary to advanced RNA secondary structure prediction algorithms, which allow good accuracy when experimental data are integrated into the prediction, tertiary structure prediction of large RNAs still remains a significant challenge. However, the field of RNA tertiary structure prediction is rapidly developing and new computational methods based on different strategies are emerging. RNAComposer is a user-friendly and freely available server for 3D structure prediction of RNA up to 500 nucleotide residues. RNAComposer employs fully automated fragment assembly based on RNA secondary structure specified by the user. Importantly, this method allows incorporation of distance restraints derived from the experimental data to strengthen the 3D predictions. The potential and limitations of RNAComposer are discussed and an application to RNA design for nanotechnology is presented.

Automated RNA 3D Structure Prediction with RNAComposer

RNAs adopt specific structures to perform their activities and these are critical to virtually all RNA-mediated processes. Because of difficulties in experimentally assessing structures of large RNAs using NMR, X-ray crystallography, or cryo-microscopy, there is currently great demand for new high-resolution 3D structure prediction methods. Recently we reported on RNAComposer, a knowledge-based method for the fully automated RNA 3D structure prediction from a user-defined secondary structure. RNAComposer method is especially suited for structural biology users. Since our initial report in 2012, both servers, freely available at http://rnacomposer.ibch.poznan.pl and http://rnacomposer.cs.put.poznan.pl have been often visited. Therefore this chapter provides guidance for using RNAComposer and discusses points that should be considered when predicting 3D RNA structure. An application example presents current scope and limitations of RNAComposer.

Elucidating the higher-order structure of biopolymers by structural probing and mass spectrometry: MS3D

Chemical probing represents a very versatile alternative for studying the structure and dynamics of substrates that are intractable by established high-resolution techniques. The implementation of MS-based strategies for the characterization of probing products has not only extended the range of applicability to virtually all types of biopolymers but has also paved the way for the introduction of new reagents that would not have been viable with traditional analytical platforms. As the availability of probing data is steadily increasing on the wings of the development of dedicated interpretation aids, powerful computational approaches have been explored to enable the effective utilization of such information to generate valid molecular models. This combination of factors has contributed to making the possibility of obtaining actual 3D structures by MS-based technologies (MS3D) a reality. Although approaches for achieving structure determination of unknown targets or assessing the dynamics of known structures may share similar reagents and development trajectories, they clearly involve distinctive experimental strategies, analytical concerns and interpretation paradigms. This Perspective offers a commentary on methods aimed at obtaining distance constraints for the modeling of full-fledged structures while highlighting common elements, salient distinctions and complementary capabilities exhibited by methods used in dynamics studies. We discuss critical factors to be addressed for completing effective structural determinations and expose possible pitfalls of chemical methods. We survey programs developed for facilitating the interpretation of experimental data and discuss possible computational strategies for translating sparse spatial constraints into all-atom models. Examples are provided to illustrate how the concerted application of very diverse probing techniques can lead to the solution of actual biological systems.

Arrangement of 3'-terminus of tRNA on the human ribosome as revealed from cross-linking data

This study is directed towards an important problem concerning the organization of the peptidyl transferase center (PTC) on the mammalian ribosome that cannot be studied by X-ray analysis since crystals of 80S ribosomes are still unavailable. Here, we investigated the arrangement of the 3'-end of tRNA in the 80S ribosomal A and P sites using a tRNA(Asp) analogue that bears a 4-thiouridine (s(4)U) attached to the 3'-terminal adenosine. It was shown that an additional nucleotide s(4)U77 on the 3'-end does not impede codon-dependent binding of the tRNA to the A and P sites of 80S ribosome. Mild UV-irradiation of the ribosomal complexes containing a short appropriately designed mRNA and the tRNA analogue resulted in cross-linking of the analogue exclusively to 28S rRNA. The cross-linking site was detected in the 4302-4540 fragment of the 28S rRNA which belongs to the highly conserved domain V that in prokaryotic ribosomes is involved in the formation of the PTC. Nucleotides cross-linked to the tRNA analogue were determined by means of reverse transcription. A comparison of the results obtained with a dynamic model of mutual arrangement of s(4)U77 of the A site tRNA and nucleotides of 23S rRNA built on the basis of an atomic model for the prokaryotic PTC led to the conclusion that environments of the tRNA 3'-terminus in prokaryotic and eukaryotic ribosomes share a significant extent of similarity, although pronounced differences are also detectable.

Ribosomal position and contacts of mRNA in eukaryotic translation initiation complexes

The position of mRNA on 40S ribosomal subunits in eukaryotic initiation complexes was determined by UV crosslinking using mRNAs containing uniquely positioned 4-thiouridines. Crosslinking of mRNA positions (+)11 to ribosomal protein (rp) rpS2(S5p) and rpS3(S3p), and (+)9-(+)11 and (+)8-(+)9 to h18 and h34 of 18S rRNA, respectively, indicated that mRNA enters the mRNA-binding channel through the same layers of rRNA and proteins as in prokaryotes. Upstream of the P-site, the proximity of positions (-)3/(-)4 to rpS5(S7p) and h23b, (-)6/(-)7 to rpS14(S11p), and (-)8-(-)11 to the 3'-terminus of 18S rRNA (mRNA/rRNA elements forming the bacterial Shine-Dalgarno duplex) also resembles elements of the bacterial mRNA path. In addition to these striking parallels, differences between mRNA paths included the proximity in eukaryotic initiation complexes of positions (+)7/(+)8 to the central region of h28, (+)4/(+)5 to rpS15(S19p), and (-)6 and (-)7/(-)10 to eukaryote-specific rpS26 and rpS28, respectively. Moreover, we previously determined that eukaryotic initiation factor2alpha (eIF2alpha) contacts position (-)3, and now report that eIF3 interacts with positions (-)8-(-)17, forming an extension of the mRNA-binding channel that likely contributes to unique aspects of eukaryotic initiation.

Developing limited proteolysis and mass spectrometry for the characterization of ribosome topography

An approach that combines limited proteolysis and matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) has been developed to probe protease-accessible sites of ribosomal proteins from intact ribosomes. Escherichia coli and Thermus thermophilus 70S ribosomes were subjected to limited proteolysis using different proteases under strictly controlled conditions. Intact ribosomal proteins and large proteolytic peptides were recovered and directly analyzed by MALDI-MS, which allows for the determination of proteins that are resistant to proteolytic digestion by accurate measurement of molecular weights. Larger proteolytic peptides can be directly identified by the combination of measured mass, enzyme specificity, and protein database searching. Sucrose density gradient centrifugation revealed that the majority of the 70S ribosome dissociates into intact 30S and 50S subunits after 120 min of limited proteolysis. Thus, examination of ribosome populations within the first 30 to 60 min of incubation provides insight into 70S structural features. Results from E. coli and T. thermophilus revealed that a significantly larger fraction of 50S ribosomal proteins have similar limited proteolysis behavior than the 30S ribosomal proteins of these two organisms. The data obtained by this approach correlate with information available from the high-resolution crystal structures of both organisms. This new approach will be applicable to investigations of other large ribonucleoprotein complexes, is readily extendable to ribosomes from other organisms, and can facilitate additional structural studies on ribosome assembly intermediates.

Internucleotide movements during formation of 16 S rRNA-rRNA photocrosslinks and their connection to the 30 S subunit conformational dynamics

UV light-induced RNA photocrosslinks are formed at a limited number of specific sites in the Escherichia coli and in other eubacterial 16 S rRNAs. To determine if unusually favorable internucleotide geometries could explain the restricted crosslinking patterns, parameters describing the internucleotide geometries were calculated from the Thermus thermophilus 30 S subunit X-ray structure and compared to crosslinking frequencies. Significant structural adjustments between the nucleotide pairs usually are needed for crosslinking. Correlations between the crosslinking frequencies and the geometrical parameters indicate that nucleotide pairs closer to the orientation needed for photoreaction have higher crosslinking frequencies. These data are consistent with transient conformational changes during crosslink formation in which the arrangements needed for photochemical reaction are attained during the electronic excitation times. The average structural rearrangement for UVA-4-thiouridine (s4U)-induced crosslinking is larger than that for UVB or UVC-induced crosslinking; this is associated with the longer excitation time for s4U and is also consistent with transient conformational changes. The geometrical parameters do not completely predict the crosslinking frequencies, implicating other aspects of the tertiary structure or conformational flexibility in determining the frequencies and the locations of the crosslinking sites. The majority of the UVB/C and UVA-s4U-induced crosslinks are located in four regions in the 30 S subunit, within or at the ends of RNA helix 34, in the tRNA P-site, in the distal end of helix 28 and in the helix 19/helix 27 region. These regions are implicated in different aspects of tRNA accommodation, translocation and in the termination reaction. These results show that photocrosslinking is an indicator for sites where there is internucleotide conformational flexibility and these sites are largely restricted to parts of the 30 S subunit associated with ribosome function.

Untying the FIV frameshifting pseudoknot structure by MS3D

The structure of the putative feline immunodeficiency virus (FIV) ribosomal frameshifting pseudoknot (PK) has been investigated by a mass spectrometric three-dimensional (MS3D) approach, which involves the application of established solvent-accessibility probes and chemical crosslinkers with detection by electrospray ionization (ESI) Fourier transform mass spectrometry (FTMS). Regardless of their size, probed substrates can be treated with ribonucleases and analyzed by ESI-FTMS to obtain the correct position of chemically modified nucleotides. Protection maps and distance information can be utilized to generate 3D models using the constraint satisfaction algorithm provided by MC-SYM and the energy minimization modules included in CNS. Control experiments were performed on a mutant of mouse mammary tumor virus pseudoknot (VPK), for which an NMR structure is available. Comparison between the MS3D model and the high-resolution structure provided a approximately 3A root-mean-square deviation calculated from all the atoms present in double-stranded regions. Applied to FIV-PK, the MS3D approach confirmed that the selected sequence could fold into an actual pseudoknot, supporting the sequence alignment predictions. Characteristic features of H-type pseudoknots were recognized immediately, but a putative A13-U30 pair was not observed at the stem junction, making FIV-PK resemble VPK more closely than the initially suggested simian retrovirus type-1 pseudoknot. In our model, the unpaired U30 protrudes into the medium, while the hinging A13 assumes a stacked conformation that enables the stems to form a approximately 60 degrees bend and relieve the strain caused by a short loop 1. The model provided the basis to explain the different alkylation patterns observed in the absence and presence of Mg(2+), suggesting the possible formation of a specific metal-binding site between loop 1 and stem 2. This instance illustrates how the MS3D model of FIV-PK can be utilized effectively to generate hypotheses and support functional observations in the absence of a high-resolution structure.

Positioning of the mRNA stop signal with respect to polypeptide chain release factors and ribosomal proteins in 80S ribosomes

To study positioning of the mRNA stop signal with respect to polypeptide chain release factors (RFs) and ribosomal components within human 80S ribosomes, photoreactive mRNA analogs were applied. Derivatives of the UUCUAAA heptaribonucleotide containing the UUC codon for Phe and the stop signal UAAA, which bore a perfluoroaryl azido group at either the fourth nucleotide or the 3'-terminal phosphate, were synthesized. The UUC codon was directed to the ribosomal P site by the cognate tRNA(Phe), targeting the UAA stop codon to the A site. Mild UV irradiation of the ternary complexes consisting of the 80S ribosome, the mRNA analog and tRNA resulted in tRNA-dependent crosslinking of the mRNA analogs to the 40S ribosomal proteins and the 18S rRNA. mRNA analogs with the photoreactive group at the fourth uridine (the first base of the stop codon) crosslinked mainly to protein S15 (and much less to S2). For the 3'-modified mRNA analog, the major crosslinking target was protein S2, while protein S15 was much less crosslinked. Crosslinking of eukaryotic (e) RF1 was entirely dependent on the presence of a stop signal in the mRNA analog. eRF3 in the presence of eRF1 did not crosslink, but decreased the yield of eRF1 crosslinking. We conclude that (i) proteins S15 and S2 of the 40S ribosomal subunit are located near the A site-bound codon; (ii) eRF1 can induce spatial rearrangement of the 80S ribosome leading to movement of protein L4 of the 60S ribosomal subunit closer to the codon located at the A site; (iii) within the 80S ribosome, eRF3 in the presence of eRF1 does not contact the stop codon at the A site and is probably located mostly (if not entirely) on the 60S subunit.

Complex of transfer-messenger RNA and elongation factor Tu. Unexpected modes of interaction

Transfer-messenger RNA (tmRNA) is a stable RNA in bacteria of 360 +/- 40 nucleotides that can be charged with alanine and can function as both tRNA and mRNA. Ribosomes that are stalled either in a coding region of mRNA or at the 3' end of an mRNA fragment lacking a stop codon are rescued by replacing their mRNA for tmRNA. Here we demonstrate that the interaction of tmRNA with the elongation factor Tu shows unexpected features. Deacylated tmRNA can form a complex with either EF-Tu.GDP or EF-Tu.GTP, the association constants are about one order of magnitude smaller than that of an Ala-tRNA.EF-Tu.GTP complex. tmRNA as well as Ala-tmRNA can be efficiently cross-linked with EF-Tu.GDP using a zero-length cross-link. The efficiency of cross-linking in the case of deacylated tmRNA does not depend on an intact CCA-3' end and is about the same, regardless whether protein mixtures such as the post-ribosomal supernatant (S100 enzymes) or purified EF-Tu are present. Two cross-linking sites with EF-Tu.GDP have been identified that are located outside the tRNA part of tmRNA, indicating an unusual interaction of tmRNA with EF-Tu.GDP.