RNA helical packing in solution: NMR structure of a 30 kDa GAAA tetraloop-receptor complex

Three Biopolymers and Origin of Life Scenarios

To track down the possible roots of life, various models for the initial living system composed of different combinations of the three extant biopolymers, RNA, DNA, and proteins, are presented. The suitability of each molecular set is assessed according to its ability to emerge autonomously, sustain, and evolve continuously towards life as we know it. The analysis incorporates current biological knowledge gained from high-resolution structural data and large sequence datasets, together with experimental results concerned with RNA replication and with the activity demonstrated by standalone constructs of the ribosomal Peptidyl Transferase Center region. The scrutiny excludes the DNA-protein combination and assigns negligible likelihood to the existence of an RNA-DNA world, as well as to an RNA world that contained a replicase made of RNA. It points to the precedence of an RNA-protein system, whose model of emergence suggests specific processes whereby a coded proto-ribosome ribozyme, specifically aminoacylated proto-tRNAs and a proto-polymerase enzyme, could have autonomously emerged, cross-catalyzing the formation of each other. This molecular set constitutes a feasible starting point for a continuous evolutionary path, proceeding via natural processes from the inanimate matter towards life as we know it.

Deep Boosted Molecular Dynamics (DBMD): Accelerating molecular simulations with Gaussian boost potentials generated using probabilistic Bayesian deep neural network

We have developed a new Deep Boosted Molecular Dynamics (DBMD) method. Probabilistic Bayesian neural network models were implemented to construct boost potentials that exhibit Gaussian distribution with minimized anharmonicity, thereby allowing for accurate energetic reweighting and enhanced sampling of molecular simulations. DBMD was demonstrated on model systems of alanine dipeptide and the fast-folding protein and RNA structures. For alanine dipeptide, 30ns DBMD simulations captured up to 83-125 times more backbone dihedral transitions than 1μs conventional molecular dynamics (cMD) simulations and were able to accurately reproduce the original free energy profiles. Moreover, DBMD sampled multiple folding and unfolding events within 300ns simulations of the chignolin model protein and identified low-energy conformational states comparable to previous simulation findings. Finally, DBMD captured a general folding pathway of three hairpin RNAs with the GCAA, GAAA, and UUCG tetraloops. Based on Deep Learning neural network, DBMD provides a powerful and generally applicable approach to boosting biomolecular simulations. DBMD is available with open source in OpenMM at https://github.com/MiaoLab20/DBMD/.

A drug discovery toolbox for Nuclear Magnetic Resonance (NMR) characterization of ligands and their targets

Information about the structure, dynamics, and ligand-binding properties of biomolecules can be derived from Nuclear Magnetic Resonance (NMR) spectroscopy and provides valuable information for drug discovery. A multitude of experimental approaches provides a wealth of information that can be tailored to the system of interest. Methods to study the behavior of ligands upon target binding enable the identification of weak binders in a robust manner that is critical for the identification of truly novel binding interactions. This is particularly important for challenging targets. Observing the solution behavior of biomolecules yields information about their structure, dynamics, and interactions. This review describes the breadth of approaches that are available, many of which are under-utilized in a drug-discovery environment, and focuses on recent advances that continue to emerge.

Strategies for RNA Resonance Assignment by 13C/15N- and 1H-Detected Solid-State NMR Spectroscopy

Magic angle spinning (MAS) solid-state NMR (ssNMR) is an established tool that can be applied to non-soluble or non-crystalline biomolecules of any size or complexity. The ssNMR method advances rapidly due to technical improvements and the development of advanced isotope labeling schemes. While ssNMR has shown significant progress in structural studies of proteins, the number of RNA studies remains limited due to ssNMR methodology that is still underdeveloped. Resonance assignment is the most critical and limiting step in the structure determination protocol that defines the feasibility of NMR studies. In this review, we summarize the recent progress in RNA resonance assignment methods and approaches for secondary structure determination by ssNMR. We critically discuss advantages and limitations of conventional 13C- and 15N-detected experiments and novel 1H-detected methods, identify optimal regimes for RNA studies by ssNMR, and provide our view on future ssNMR studies of RNA in large RNP complexes.

Structures of artificially designed discrete RNA nanoarchitectures at near-atomic resolution

Although advances in nanotechnology have enabled the construction of complex and functional synthetic nucleic acid–based nanoarchitectures, high-resolution discrete structures are lacking because of the difficulty in obtaining good diffracting crystals. Here, we report the design and construction of RNA nanostructures based on homooligomerizable one-stranded tiles for x-ray crystallographic determination. We solved three structures to near-atomic resolution: a 2D parallelogram, a 3D nanobracelet unexpectedly formed from an RNA designed for a nanocage, and, eventually, a bona fide 3D nanocage designed with the guidance of the two previous structures. Structural details of their constituent motifs, such as kissing loops, branched kissing loops, and T-junctions, that resemble natural RNA motifs and resisted x-ray determination are revealed, providing insights into those natural motifs. This work unveils the largely unexplored potential of crystallography in gaining high-resolution feedback for nanoarchitectural design and suggests a route to investigate RNA motif structures by configuring them into nanoarchitectures.

High-throughput dissection of the thermodynamic and conformational properties of a ubiquitous class of RNA tertiary contact motifs

Despite RNA's diverse secondary and tertiary structures and its complex conformational changes, nature utilizes a limited set of structural "motifs"-helices, junctions, and tertiary contact modules-to build diverse functional RNAs. Thus, in-depth descriptions of a relatively small universe of RNA motifs may lead to predictive models of RNA tertiary conformational landscapes. Motifs may have different properties depending on sequence and secondary structure, giving rise to subclasses that expand the universe of RNA building blocks. Yet we know very little about motif subclasses, given the challenges in mapping conformational properties in high throughput. Previously, we used "RNA on a massively parallel array" (RNA-MaP), a quantitative, high-throughput technique, to study thousands of helices and two-way junctions. Here, we adapt RNA-MaP to study the thermodynamic and conformational properties of tetraloop/tetraloop receptor (TL/TLR) tertiary contact motifs, analyzing 1,493 TLR sequences from different classes. Clustering analyses revealed variability in TL specificity, stability, and conformational behavior. Nevertheless, natural GAAA/11ntR TL/TLRs, while varying in tertiary stability by ∼2.5 kcal/mol, exhibited conserved TL specificity and conformational properties. Thus, RNAs may tune stability without altering the overall structure of these TL/TLRs. Furthermore, their stability correlated with natural frequency, suggesting thermodynamics as the dominant selection pressure. In contrast, other TL/TLRs displayed heterogenous conformational behavior and appear to not be under strong thermodynamic selection. Our results build toward a generalizable model of RNA-folding thermodynamics based on the properties of isolated motifs, and our characterized TL/TLR library can be used to engineer RNAs with predictable thermodynamic and conformational behavior.

In vitro selected GUAA tetraloop-binding receptors with structural plasticity and evolvability towards natural RNA structural modules

GNRA tetraloop-binding receptor interactions are key components in the macromolecular assembly of a variety of functional RNAs. In nature, there is an apparent bias for GAAA/11nt receptor and GYRA/helix interactions, with the former interaction being thermodynamically more stable than the latter. While past in vitro selections allowed isolation of novel GGAA and GUGA receptors, we report herein an in vitro selection that revealed several novel classes of specific GUAA receptors with binding affinities comparable to those from natural GAAA/11nt interactions. These GUAA receptors have structural homology with double-locked bulge RNA modules naturally occurring in ribosomal RNAs. They display mutational robustness that enables exploration of the sequence/phenotypic space associated to GNRA/receptor interactions through epistasis. Their thermodynamic self-assembly fitness landscape is characterized by a rugged neutral network with possible evolutionary trajectories toward natural GNRA/receptor interactions. High throughput sequencing analysis revealed synergetic mutations located away from the tertiary interactions that positively contribute to assembly fitness. Our study suggests that the repertoire of GNRA/receptor interactions is much larger than initially thought from the analysis of natural stable RNA molecules and also provides clues for their evolution towards natural GNRA/receptors.

NMR Spectroscopy of Large Functional RNAs: From Sample Preparation to Low-Gamma Detection

NMR spectroscopy is a potent method for the structural and biophysical characterization of RNAs. The application of NMR spectroscopy is restricted in RNA size and most often requires isotope-labeled or even selectively labeled RNAs. Additionally, new NMR pulse sequences, such as the heteronuclear-detected NMR experiments, are introduced. We herein provide detailed protocols for the preparation of isotope-labeled RNA for NMR spectroscopy via in vitro transcription. This protocol covers all steps, from the preparation of DNA template to the transcription of milligram RNA quantities. Moreover, we present a protocol for a chemo-enzymatic approach to introduce a single modified nucleotide at any position of any RNA. Regarding NMR methodology, we share protocols for the implementation of a suite of heteronuclear-detected NMR experiments including ¹³ C-detected experiments for ribose assignment and amino groups, the CN-spin filter heteronuclear single quantum coherence (HSQC) for imino groups and the ¹⁵ N-detected band-selective excitation short transient transverse-relaxation-optimized spectroscopy (BEST-TROSY) experiment. © 2020 The Authors. Basic Protocol 1: Preparation of isotope-labeled RNA samples with in vitro transcription using T7 RNAP, DEAE chromatography, and RP-HPLC purification Alternate Protocol 1: Purification of isotope-labeled RNA from in vitro transcription with preparative PAGE Alternate Protocol 2: Purification of isotope-labeled RNA samples from in vitro transcription via centrifugal concentration Support Protocol 1: Preparation of DNA template from plasmid Support Protocol 2: Preparation of PCR DNA as template Support Protocol 3: Preparation of T7 RNA Polymerase (T7 RNAP) Support Protocol 4: Preparation of yeast inorganic pyrophosphatase (YIPP) Basic Protocol 2: Preparation of site-specific labeled RNAs using a chemo-enzymatic synthesis Support Protocol 5: Synthesis of modified nucleoside 3',5'-bisphosphates Support Protocol 6: Preparation of T4 RNA Ligase 2 Support Protocol 7: Setup of NMR spectrometer for heteronuclear-detected NMR experiments Support Protocol 8: IPAP and DIPAP for homonuclear decoupling Basic Protocol 3: ¹³ C-detected 3D (H)CC-TOCSY, (H)CPC, and (H)CPC-CCH-TOCSY experiments for ribose assignment Basic Protocol 4: ¹³ C-detected 2D CN-spin filter HSQC experiment Basic Protocol 5: ¹³ C-detected C(N)H-HDQC experiment for the detection of amino groups Support Protocol 9: ¹³ C-detected CN-HSQC experiment for amino groups Basic Protocol 6: ¹³ C-detected "amino"-NOESY experiment Basic Protocol 7: ¹⁵ N-detected BEST-TROSY experiment.

Structural Insights into RNA Dimerization: Motifs, Interfaces and Functions

In comparison with the pervasive use of protein dimers and multimers in all domains of life, functional RNA oligomers have so far rarely been observed in nature. Their diminished occurrence contrasts starkly with the robust intrinsic potential of RNA to multimerize through long-range base-pairing ("kissing") interactions, self-annealing of palindromic or complementary sequences, and stable tertiary contact motifs, such as the GNRA tetraloop-receptors. To explore the general mechanics of RNA dimerization, we performed a meta-analysis of a collection of exemplary RNA homodimer structures consisting of viral genomic elements, ribozymes, riboswitches, etc., encompassing both functional and fortuitous dimers. Globally, we found that domain-swapped dimers and antiparallel, head-to-tail arrangements are predominant architectural themes. Locally, we observed that the same structural motifs, interfaces and forces that enable tertiary RNA folding also drive their higher-order assemblies. These feature prominently long-range kissing loops, pseudoknots, reciprocal base intercalations and A-minor interactions. We postulate that the scarcity of functional RNA multimers and limited diversity in multimerization motifs may reflect evolutionary constraints imposed by host antiviral immune surveillance and stress sensing. A deepening mechanistic understanding of RNA multimerization is expected to facilitate investigations into RNA and RNP assemblies, condensates, and granules and enable their potential therapeutical targeting.

Responsive self-assembly of tectoRNAs with loop-receptor interactions from the tetrahydrofolate (THF) riboswitch

Naturally occurring RNAs are known to exhibit a high degree of modularity, whereby specific structural modules (or motifs) can be mixed and matched to create new molecular architectures. The modular nature of RNA also affords researchers the ability to characterize individual structural elements in controlled synthetic contexts in order to gain new and critical insights into their particular structural features and overall performance. Here, we characterized the binding affinity of a unique loop–receptor interaction found in the tetrahydrofolate (THF) riboswitch using rationally designed self-assembling tectoRNAs. Our work suggests that the THF loop–receptor interaction has been fine-tuned for its particular role as a riboswitch component. We also demonstrate that the thermodynamic stability of this interaction can be modulated by the presence of folinic acid, which induces a local structural change at the level of the loop–receptor. This corroborates the existence of a THF binding site within this tertiary module and paves the way for its potential use as a THF responsive module for RNA nanotechnology and synthetic biology.

Deducing putative ancestral forms of GNRA/receptor interactions from the ribosome

Stable RNAs rely on a vast repertoire of long-range interactions to assist in the folding of complex cellular machineries such as the ribosome. The universally conserved L39/H89 interaction is a long-range GNRA-like/receptor interaction localized in proximity to the peptidyl transferase center of the large subunit of the ribosome. Because of its central location, L39/H89 likely originated at an early evolutionary stage of the ribosome and played a significant role in its early function. However, L39/H89 self-assembly is impaired outside the ribosomal context. Herein, we demonstrate that structural modularity principles can be used to re-engineer L39/H89 to self-assemble in vitro. The new versions of L39/H89 improve affinity and loop selectivity by several orders of magnitude and retain the structural and functional features of their natural counterparts. These versions of L39/H89 are proposed to be ancestral forms of L39/H89 that were capable of assembling and folding independently from proteins and post-transcriptional modifications. This work demonstrates that novel RNA modules can be rationally designed by taking advantage of the modular syntax of RNA. It offers the prospect of creating new biochemical models of the ancestral ribosome and increases the tool kit for RNA nanotechnology and synthetic biology.

Small-Angle Scattering as a Structural Probe for Nucleic Acid Nanoparticles (NANPs) in a Dynamic Solution Environment

Nucleic acid-based technologies are an emerging research focus area for pharmacological and biological studies because they are biocompatible and can be designed to produce a variety of scaffolds at the nanometer scale. The use of nucleic acids (ribonucleic acid (RNA) and/or deoxyribonucleic acid (DNA)) as building materials in programming the assemblies and their further functionalization has recently established a new exciting field of RNA and DNA nanotechnology, which have both already produced a variety of different functional nanostructures and nanodevices. It is evident that the resultant architectures require detailed structural and functional characterization and that a variety of technical approaches must be employed to promote the development of the emerging fields. Small-angle X-ray and neutron scattering (SAS) are structural characterization techniques that are well placed to determine the conformation of nucleic acid nanoparticles (NANPs) under varying solution conditions, thus allowing for the optimization of their design. SAS experiments provide information on the overall shapes and particle dimensions of macromolecules and are ideal for following conformational changes of the molecular ensemble as it behaves in solution. In addition, the inherent differences in the neutron scattering of nucleic acids, lipids, and proteins, as well as the different neutron scattering properties of the isotopes of hydrogen, combined with the ability to uniformly label biological macromolecules with deuterium, allow one to characterize the conformations and relative dispositions of the individual components within an assembly of biomolecules. This article will review the application of SAS methods and provide a summary of their successful utilization in the emerging field of NANP technology to date, as well as share our vision on its use in complementing a broad suite of structural characterization tools with some simulated results that have never been shared before.

Advances that facilitate the study of large RNA structure and dynamics by nuclear magnetic resonance spectroscopy

The characterization of functional yet nonprotein coding (nc) RNAs has expanded the role of RNA in the cell from a passive player in the central dogma of molecular biology to an active regulator of gene expression. The misregulation of ncRNA function has been linked with a variety of diseases and disorders ranging from cancers to neurodegeneration. However, a detailed molecular understanding of how ncRNAs function has been limited; due, in part, to the difficulties associated with obtaining high-resolution structures of large RNAs. Tertiary structure determination of RNA as a whole is hampered by various technical challenges, all of which are exacerbated as the size of the RNA increases. Namely, RNAs tend to be highly flexible and dynamic molecules, which are difficult to crystallize. Biomolecular nuclear magnetic resonance (NMR) spectroscopy offers a viable alternative to determining the structure of large RNA molecules that do not readily crystallize, but is itself hindered by some technical limitations. Recently, a series of advancements have allowed the biomolecular NMR field to overcome, at least in part, some of these limitations. These advances include improvements in sample preparation strategies as well as methodological improvements. Together, these innovations pave the way for the study of ever larger RNA molecules that have important biological function. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

Computational modeling of RNA 3D structure based on experimental data

RNA molecules are master regulators of cells. They are involved in a variety of molecular processes: they transmit genetic information, sense cellular signals and communicate responses, and even catalyze chemical reactions. As in the case of proteins, RNA function is dictated by its structure and by its ability to adopt different conformations, which in turn is encoded in the sequence. Experimental determination of high-resolution RNA structures is both laborious and difficult, and therefore the majority of known RNAs remain structurally uncharacterized. To address this problem, predictive computational methods were developed based on the accumulated knowledge of RNA structures determined so far, the physical basis of the RNA folding, and taking into account evolutionary considerations, such as conservation of functionally important motifs. However, all theoretical methods suffer from various limitations, and they are generally unable to accurately predict structures for RNA sequences longer than 100-nt residues unless aided by additional experimental data. In this article, we review experimental methods that can generate data usable by computational methods, as well as computational approaches for RNA structure prediction that can utilize data from experimental analyses. We outline methods and data types that can be potentially useful for RNA 3D structure modeling but are not commonly used by the existing software, suggesting directions for future development.

The Story of RNA Folding, as Told in Epochs

The past decades have witnessed tremendous developments in our understanding of RNA biology. At the core of these advances have been studies aimed at discerning RNA structure and at understanding the forces that influence the RNA folding process. It is easy to take the present state of understanding for granted, but there is much to be learned by considering the path to our current understanding, which has been tortuous, with the birth and death of models, the adaptation of experimental tools originally developed for characterization of protein structure and catalysis, and the development of novel tools for probing RNA. In this review we tour the stages of RNA folding studies, considering them as "epochs" that can be generalized across scientific disciplines. These epochs span from the discovery of catalytic RNA, through biophysical insights into the putative primordial RNA World, to characterization of structured RNAs, the building and testing of models, and, finally, to the development of models with the potential to yield generalizable predictive and quantitative models for RNA conformational, thermodynamic, and kinetic behavior. We hope that this accounting will aid others as they navigate the many fascinating questions about RNA and its roles in biology, in the past, present, and future.

Isotope labeling for studying RNA by solid-state NMR spectroscopy

Nucleic acids play key roles in most biological processes, either in isolation or in complex with proteins. Often they are difficult targets for structural studies, due to their dynamic behavior and high molecular weight. Solid-state nuclear magnetic resonance spectroscopy (ssNMR) provides a unique opportunity to study large biomolecules in a non-crystalline state at atomic resolution. Application of ssNMR to RNA, however, is still at an early stage of development and presents considerable challenges due to broad resonances and poor dispersion. Isotope labeling, either as nucleotide-specific, atom-specific or segmental labeling, can resolve resonance overlaps and reduce the line width, thus allowing ssNMR studies of RNA domains as part of large biomolecules or complexes. In this review we discuss the methods for RNA production and purification as well as numerous approaches for isotope labeling of RNA. Furthermore, we give a few examples that emphasize the instrumental role of isotope labeling and ssNMR for studying RNA as part of large ribonucleoprotein complexes.

Accurate Measurement of Residual Dipolar Couplings in Large RNAs by Variable Flip Angle NMR

NMR approaches using nucleotide-specific deuterium labeling schemes have enabled structural studies of biologically relevant RNAs of increasing size and complexity. Although local structure is well-determined using these methods, definition of global structural features, including relative orientations of independent helices, remains a challenge. Residual dipolar couplings, a potential source of orientation information, have not been obtainable for large RNAs due to poor sensitivity resulting from rapid heteronuclear signal decay. Here we report a novel multiple quantum NMR method for RDC determination that employs flip angle variation rather than a coupling evolution period. The accuracy of the method and its utility for establishing interhelical orientations are demonstrated for a 36-nucleotide RNA, for which comparative data could be obtained. Applied to a 78 kDa Rev response element from the HIV-1 virus, which has an effective rotational correlation time of ca. 160 ns, the method yields sensitivity gains of an order of magnitude or greater over existing approaches. Solution-state access to structural organization in RNAs of at least 230 nucleotides is now possible.

Influence of Na+ and Mg2+ ions on RNA structures studied with molecular dynamics simulations

The structure of ribonucleic acid (RNA) polymers is strongly dependent on the presence of, in particular Mg2+ cations to stabilize structural features. Only in high-resolution X-ray crystallography structures can ions be identified reliably. Here, we perform molecular dynamics simulations of 24 RNA structures with varying ion concentrations. Twelve of the structures were helical and the others complex folded. The aim of the study is to predict ion positions but also to evaluate the impact of different types of ions (Na+ or Mg2+) and the ionic strength on structural stability and variations of RNA. As a general conclusion Mg2+ is found to conserve the experimental structure better than Na+ and, where experimental ion positions are available, they can be reproduced with reasonable accuracy. If a large surplus of ions is present the added electrostatic screening makes prediction of binding-sites less reproducible. Distinct differences in ion-binding between helical and complex folded structures are found. The strength of binding (ΔG‡ for breaking RNA atom-ion interactions) is found to differ between roughly 10 and 26 kJ/mol for the different RNA atoms. Differences in stability between helical and complex folded structures and of the influence of metal ions on either are discussed.

Single-Molecule Fluorescence Reveals Commonalities and Distinctions among Natural and in Vitro-Selected RNA Tertiary Motifs in a Multistep Folding Pathway

Decades of study of the RNA folding problem have revealed that diverse and complex structured RNAs are built from a common set of recurring structural motifs, leading to the perspective that a generalizable model of RNA folding may be developed from understanding of the folding properties of individual structural motifs. We used single-molecule fluorescence to dissect the kinetic and thermodynamic properties of a set of variants of a common tertiary structural motif, the tetraloop/tetraloop-receptor (TL/TLR). Our results revealed a multistep TL/TLR folding pathway in which preorganization of the ubiquitous AA-platform submotif precedes the formation of the docking transition state and tertiary A-minor hydrogen bond interactions form after the docking transition state. Differences in ion dependences between TL/TLR variants indicated the occurrence of sequence-dependent conformational rearrangements prior to and after the formation of the docking transition state. Nevertheless, varying the junction connecting the TL/TLR produced a common kinetic and ionic effect for all variants, suggesting that the global conformational search and compaction electrostatics are energetically independent from the formation of the tertiary motif contacts. We also found that in vitro-selected variants, despite their similar stability at high Mg²⁺ concentrations, are considerably less stable than natural variants under near-physiological ionic conditions, and the occurrence of the TL/TLR sequence variants in Nature correlates with their thermodynamic stability in isolation. Overall, our findings are consistent with modular but complex energetic properties of RNA structural motifs and will aid in the eventual quantitative description of RNA folding from its secondary and tertiary structural elements.

A homodimer interface without base pairs in an RNA mimic of red fluorescent protein

Corn, a 28-nucleotide RNA, increases yellow fluorescence of its cognate ligand 3,5-difluoro-4-hydroxybenzylidene-imidazolinone-2-oxime (DFHO) by >400-fold. Corn was selected in vitro to overcome limitations of other fluorogenic RNAs, particularly rapid photobleaching. We now report the Corn-DFHO co-crystal structure, discovering that the functional species is a quasisymmetric homodimer. Unusually, the dimer interface, in which six unpaired adenosines break overall two-fold symmetry, lacks any intermolecular base pairs. The homodimer encapsulates one DFHO at its interprotomer interface, sandwiching it with a G-quadruplex from each protomer. Corn and the green-fluorescent Spinach RNA are structurally unrelated. Their convergent use of G-quadruplexes underscores the usefulness of this motif for RNA-induced small-molecule fluorescence. The asymmetric dimer interface of Corn could provide a basis for the development of mutants that only fluoresce as heterodimers. Such variants would be analogous to Split GFP, and may be useful for analyzing RNA co-expression or association, or for designing self-assembling RNA nanostructures.

Applications of NMR to structure determination of RNAs large and small

Nuclear magnetic resonance (NMR) spectroscopy is a powerful tool to investigate the structure and dynamics of RNA, because many biologically important RNAs have conformationally flexible structures, which makes them difficult to crystallize. Functional, independently folded RNA domains, range in size between simple stem-loops of as few as 10-20 nucleotides, to 50-70 nucleotides, the size of tRNA and many small ribozymes, to a few hundred nucleotides, the size of more complex RNA enzymes and of the functional domains of non-coding transcripts. In this review, we discuss new methods for sample preparation, assignment strategies and structure determination for independently folded RNA domains of up to 100 kDa in molecular weight.

A New Method to Predict Ion Effects in RNA Folding

The strong interaction between metal ions in solution and highly charged RNA molecules is critical for RNA structure formation and stabilization. Metal ions binding to RNA can induce RNA structural changes that are important for RNA cellular functions. Therefore, quantitative modeling of the ion effects is essential for RNA structure prediction and RNA-based molecular design. Recently, inspired by the increasing experimental evidence that supports the importance of ion correlation and fluctuation in ion-RNA interactions, we developed a new computational model, Monte Carlo Tightly Bound Ion (MCTBI) model. The validity of the model is shown by the improved accuracy in the predictions for ion binding properties and ion-dependent free energies for RNA structures. In this chapter, using homodimeric tetraloop-receptor docking as an illustrative example, we showcase the MCTBI method for the computational prediction of the ion effects in RNA folding.

Could a Proto-Ribosome Emerge Spontaneously in the Prebiotic World?

An indispensable prerequisite for establishing a scenario of life emerging by natural processes is the requirement that the first simple proto-molecules could have had a realistic probability of self-assembly from random molecular polymers in the prebiotic world. The vestige of the proto-ribosome, which is believed to be still embedded in the contemporary ribosome, is used to assess the feasibility of such spontaneous emergence. Three concentric structural elements of different magnitudes, having a dimeric nature derived from the symmetrical region of the ribosomal large subunit, were suggested to constitute the vestige of the proto-ribosome. It is assumed to have materialized spontaneously in the prebiotic world, catalyzing non-coded peptide bond formation and simple elongation. Probabilistic and energetic considerations are applied in order to evaluate the suitability of the three contenders for being the initial proto-ribosome. The analysis points to the simplest proto-ribosome, comprised of a dimer of tRNA-like molecules presently embedded in the core of the symmetrical region, as the only one having a realistic statistical likelihood of spontaneous emergence from random RNA chains. Hence it offers a feasible starting point for a continuous evolutionary path from the prebiotic matter, through natural processes, into the intricate modern translation system.

NMR solution structure determination of large RNA-protein complexes

Structure determination of RNA-protein complexes is essential for our understanding of the multiple layers of RNA-mediated posttranscriptional regulation of gene expression. Over the past 20years, NMR spectroscopy became a key tool for structural studies of RNA-protein interactions. Here, we review the progress being made in NMR structure determination of large ribonucleoprotein assemblies. We discuss approaches for the design of RNA-protein complexes for NMR structural studies, established and emerging isotope and segmental labeling schemes suitable for large RNPs and how to gain distance restraints from NOEs, PREs and EPR and orientational information from RDCs and SAXS/SANS in such systems. The new combination of NMR measurements with MD simulations and its potential will also be discussed. Application and combination of these various methods for structure determination of large RNPs will be illustrated with three large RNA-protein complexes (>40kDa) and other interesting complexes determined in the past six and a half years.

Folding of RNA tertiary structure: Linkages between backbone phosphates, ions, and water

The functional forms of many RNAs have compact architectures. The placement of phosphates within such structures must be influenced not only by the strong electrostatic repulsion between phosphates, but also by networks of interactions between phosphates, water, and mobile ions. This review first explores what has been learned of the basic thermodynamic constraints on these arrangements from studies of hydration and ions in simple DNA molecules, and then gives an overview of what is known about ion and water interactions with RNA structures. A brief survey of RNA crystal structures identifies several interesting architectures in which closely spaced phosphates share hydration shells or phosphates are buried in environments that provide intramolecular hydrogen bonds or site-bound cations. Formation of these structures must require strong coupling between the uptake of ions and release of water.

Chemo-enzymatic labeling for rapid assignment of RNA molecules

Even though Nuclear Magnetic Resonance (NMR) spectroscopy is one of the few techniques capable of determining atomic resolution structures of RNA, it is constrained by two major problems of chemical shift overlap of resonances and rapid signal loss due to line broadening. Emerging tools to tackle these problems include synthesis of atom specifically labeled or chemically modified nucleotides. Herein we review the synthesis of these nucleotides, the design and production of appropriate RNA samples, and the application and analysis of the NMR experiments that take advantage of these labels.

Studying metal ion binding properties of a three-way junction RNA by heteronuclear NMR

Self-splicing group II introns are highly structured RNA molecules, containing a characteristic secondary and catalytically active tertiary structure, which is formed only in the presence of Mg(II). Mg(II) initiates the first folding step governed by the κζ element within domain 1 (D1κζ). We recently solved the NMR structure of D1κζ derived from the mitochondrial group II intron ribozyme Sc.ai5γ and demonstrated that Mg(II) is essential for its stabilization. Here, we performed a detailed multinuclear NMR study of metal ion interactions with D1κζ, using Cd(II) and cobalt(III)hexammine to probe inner- and outer-sphere coordination of Mg(II) and thus to better characterize its binding sites. Accordingly, we mapped (1)H, (15)N, (13)C, and (31)P spectral changes upon addition of different amounts of the metal ions. Our NMR data reveal a Cd(II)-assisted macrochelate formation at the 5'-end triphosphate, a preferential Cd(II) binding to guanines in a helical context, an electrostatic interaction in the ζ tetraloop receptor and various metal ion interactions in the GAAA tetraloop and κ element. These results together with our recently published data on Mg(II) interaction provide a much better understanding of Mg(II) binding to D1κζ, and reveal how intricate and complex metal ion interactions can be.

(19)F-labeling of the adenine H2-site to study large RNAs by NMR spectroscopy

In comparison to proteins and protein complexes, the size of RNA amenable to NMR studies is limited despite the development of new isotopic labeling strategies including deuteration and ligation of differentially labeled RNAs. Due to the restricted chemical shift dispersion in only four different nucleotides spectral resolution remains limited in larger RNAs. Labeling RNAs with the NMR-active nucleus (19)F has previously been introduced for small RNAs up to 40 nucleotides (nt). In the presented work, we study the natural occurring RNA aptamer domain of the guanine-sensing riboswitch comprising 73 nucleotides from Bacillus subtilis. The work includes protocols for improved in vitro transcription of 2-fluoroadenosine-5'-triphosphat (2F-ATP) using the mutant P266L of the T7 RNA polymerase. Our NMR analysis shows that the secondary and tertiary structure of the riboswitch is fully maintained and that the specific binding of the cognate ligand hypoxanthine is not impaired by the introduction of the (19)F isotope. The thermal stability of the (19)F-labeled riboswitch is not altered compared to the unmodified sequence, but local base pair stabilities, as measured by hydrogen exchange experiments, are modulated. The characteristic change in the chemical shift of the imino resonances detected in a (1)H,(15)N-HSQC allow the identification of Watson-Crick base paired uridine signals and the (19)F resonances can be used as reporters for tertiary and secondary structure transitions, confirming the potential of (19)F-labeling even for sizeable RNAs in the range of 70 nucleotides.

Application of Nuclear Magnetic Resonance and Hybrid Methods to Structure Determination of Complex Systems

The current main challenge of Structural Biology is to undertake the structure determination of increasingly complex systems in the attempt to better understand their biological function. As systems become more challenging, however, there is an increasing demand for the parallel use of more than one independent technique to allow pushing the frontiers of structure determination and, at the same time, obtaining independent structural validation. The combination of different Structural Biology methods has been named hybrid approaches. The aim of this review is to critically discuss the most recent examples and new developments that have allowed structure determination or experimentally-based modelling of various molecular complexes selecting them among those that combine the use of nuclear magnetic resonance and small angle scattering techniques. We provide a selective but focused account of some of the most exciting recent approaches and discuss their possible further developments.

NMR Methods for Characterization of RNA Secondary Structure

Knowledge of RNA secondary structure is often sufficient to identify relationships between the structure of RNA and processing pathways, and the design of therapeutics. Nuclear magnetic resonance (NMR) can identify types of nucleotide base pairs and the sequence, thus limiting possible secondary structures. Because NMR experiments, like chemical mapping, are performed in solution, not in single crystals, experiments can be initiated as soon as the biomolecule is expressed and purified. This chapter summarizes NMR methods that permit rapid identification of RNA secondary structure, information that can be used as supplements to chemical mapping, and/or as preliminary steps required for 3D structure determination. The primary aim is to provide guidelines to enable a researcher with minimal knowledge of NMR to quickly extract secondary structure information from basic datasets. Instrumental and sample considerations that can maximize data quality are discussed along with some details for optimal data acquisition and processing parameters. Approaches for identifying base pair types in both unlabeled and isotopically labeled RNA are covered. Common problems, such as missing signals and overlaps, and approaches to address them are considered. Programs under development for merging NMR data with structure prediction algorithms are briefly discussed.

Metal Ion-Mediated Nucleobase Recognition by the ZTP Riboswitch

The ZTP riboswitch is a widespread family of regulatory RNAs that upregulate de novo purine synthesis in response to increased intracellular levels of ZTP or ZMP. As an important intermediate in purine biosynthesis, ZMP also serves as a proxy for the concentration of N10-formyl-tetrahydrofolate, a key component of one-carbon metabolism. Here, we report the structure of the ZTP riboswitch bound to ZMP at a resolution of 1.80 Å. The RNA contains two subdomains brought together through a long-range pseudoknot further stabilized through helix-helix packing. ZMP is bound at the subdomain interface of the RNA through a set of interactions with the base, ribose sugar, and phosphate moieties of the ligand. Unique to nucleobase recognition by RNAs, the Z base is inner-sphere coordinated to a magnesium cation bound by two backbone phosphates. This interaction, along with steric hindrance by the backbone, imparts specificity over chemically similar compounds such as ATP/AMP.

Structural features of a 3' splice site in influenza a

Influenza A is an RNA virus with a genome of eight negative sense segments. Segment 7 mRNA contains a 3′ splice site for alternative splicing to encode the essential M2 protein. On the basis of sequence alignment and chemical mapping experiments, the secondary structure surrounding the 3′ splice site has an internal loop, adenine bulge, and hairpin loop when it is in the hairpin conformation that exposes the 3′ splice site. We report structural features of a three-dimensional model of the hairpin derived from nuclear magnetic resonance spectra and simulated annealing with restrained molecular dynamics. Additional insight was provided by modeling based on ¹H chemical shifts. The internal loop containing the 3′ splice site has a dynamic guanosine and a stable imino (cis Watson–Crick/Watson–Crick) GA pair. The adenine bulge also appears to be dynamic with the A either stacked in the stem or forming a base triple with a Watson–Crick GC pair. The hairpin loop is a GAAA tetraloop closed by an AC pair.

Silver nanoclusters for RNA nanotechnology: steps towards visualization and tracking of RNA nanoparticle assemblies

The growing interest in designing functionalized, RNA-based nanoparticles (NPs) for applications such as cancer therapeutics requires simple, efficient assembly assays. Common methods for tracking RNA assemblies such as native polyacrylamide gels and atomic force microscopy are often time-intensive and, therefore, undesirable. Here we describe a technique for rapid analysis of RNA NP assembly stages using the formation of fluorescent silver nanoclusters (Ag NCs). This method exploits the single-stranded specificity and sequence dependence of Ag NC formation to produce unique optical readouts for each stage of RNA NP assembly, obtained readily after synthesis.

Progress and outlook in structural biology of large viral RNAs

The field of viral molecular biology has reached a precipice for which pioneering studies on the structure of viral RNAs are beginning to bridge the gap. It has become clear that viral genomic RNAs are not simply carriers of hereditary information, but rather are active players in many critical stages during replication. Indeed, functions such as cap-independent translation initiation mechanisms are, in some cases, primarily driven by RNA structural determinants. Other stages including reverse transcription initiation in retroviruses, nuclear export and viral packaging are specifically dependent on the proper 3-dimensional folding of multiple RNA domains to recruit necessary viral and host factors required for activity. Furthermore, a large-scale conformational change within the 5'-untranslated region of HIV-1 has been proposed to regulate the temporal switch between viral protein synthesis and packaging. These RNA-dependent functions are necessary for replication of many human disease-causing viruses such as severe acute respiratory syndrome (SARS)-associated coronavirus, West Nile virus, and HIV-1. The potential for antiviral development is currently hindered by a poor understanding of RNA-driven molecular mechanisms, resulting from a lack of structural information on large RNAs and ribonucleoprotein complexes. Herein, we describe the recent progress that has been made on characterizing these large RNAs and provide brief descriptions of the techniques that will be at the forefront of future advances. Ongoing and future work will contribute to a more complete understanding of the lifecycles of retroviruses and RNA viruses and potentially lead to novel antiviral strategies.

Structure determination of noncanonical RNA motifs guided by ¹H NMR chemical shifts

Structured noncoding RNAs underlie fundamental cellular processes, but determining their three-dimensional structures remains challenging. We demonstrate that integrating ¹H NMR chemical shift data with Rosetta de novo modeling can be used to consistently determine high-resolution RNA structures. On a benchmark set of 23 noncanonical RNA motifs, including 11 'blind' targets, chemical-shift Rosetta for RNA (CS-Rosetta-RNA) recovered experimental structures with high accuracy (0.6-2.0 Å all-heavy-atom r.m.s. deviation) in 18 cases.

Exploring the electrostatic energy landscape for tetraloop-receptor docking

It has long been appreciated that Mg(2+) is essential for the stabilization of RNA tertiary structure. However, the problem of quantitative prediction for the ion effect in tertiary structure folding remains. By using the virtual bond RNA folding model (Vfold) to generate RNA conformations and the newly improved tightly bound ion model (TBI) to treat ion-RNA interactions, we investigate Mg(2+)-facilitated tetraloop-receptor docking. For the specific construct of the tetraloop-receptor system, the theoretical analysis shows that the Mg(2+)-induced stabilizing force for the docked state is predominantly entropic and the major contribution comes from the entropy of the diffusive ions. Furthermore, our results show that Mg(2+) ions promote tetraloop-receptor docking mainly through the entropy of the diffusive ions. The theoretical prediction agrees with experimental analysis. The method developed in this paper, which combines the theory for the (Mg(2+)) ion effects in RNA folding and RNA conformational sampling, may provide a useful framework for studying the ion effect in the folding of more complex RNA structures.

Ebola virus RNA editing depends on the primary editing site sequence and an upstream secondary structure

Ebolavirus (EBOV), the causative agent of a severe hemorrhagic fever and a biosafety level 4 pathogen, increases its genome coding capacity by producing multiple transcripts encoding for structural and nonstructural glycoproteins from a single gene. This is achieved through RNA editing, during which non-template adenosine residues are incorporated into the EBOV mRNAs at an editing site encoding for 7 adenosine residues. However, the mechanism of EBOV RNA editing is currently not understood. In this study, we report for the first time that minigenomes containing the glycoprotein gene editing site can undergo RNA editing, thereby eliminating the requirement for a biosafety level 4 laboratory to study EBOV RNA editing. Using a newly developed dual-reporter minigenome, we have characterized the mechanism of EBOV RNA editing, and have identified cis-acting sequences that are required for editing, located between 9 nt upstream and 9 nt downstream of the editing site. Moreover, we show that a secondary structure in the upstream cis-acting sequence plays an important role in RNA editing. EBOV RNA editing is glycoprotein gene-specific, as a stretch encoding for 7 adenosine residues located in the viral polymerase gene did not serve as an editing site, most likely due to an absence of the necessary cis-acting sequences. Finally, the EBOV protein VP30 was identified as a trans-acting factor for RNA editing, constituting a novel function for this protein. Overall, our results provide novel insights into the RNA editing mechanism of EBOV, further understanding of which might result in novel intervention strategies against this viral pathogen.

Comparison of interactions of diamine and Mg²⁺ with RNA tertiary structures: similar versus differential effects on the stabilities of diverse RNA folds

Cations play a large role in stabilizing the native state of RNA in vivo. In addition to Mg²⁺, putrescine²⁺ is an abundant divalent cation in bacterial cells, but its effect on the folding of RNA tertiary structure has not been widely explored. In this study, we look at how the stabilities of four structured RNAs, each with a different degree of dependence on K⁺ and Mg²⁺, are affected by putrescine²⁺ relative to Mg²⁺. Through the use of thermal melts, we observe that (i) at a given concentration, putrescine²⁺ is less effective than Mg²⁺ at stabilizing RNA, (ii) the stability imparted to RNA by various diamines is a function of charge density (average separation distance between charges) as well as the flexibility of the counterion, and (iii) when Mg²⁺ is already present in a buffer, further addition of putrescine²⁺ may either destabilize or stabilize RNA structure, depending on whether the native RNA does or does not chelate Mg²⁺ ion, respectively. At ion concentrations likely to be found in vivo, the effect of putrescine²⁺ on the free energy of folding of an RNA tertiary structure is probably quite small compared to that of Mg²⁺, but the ability of mixed Mg²⁺/putrescine²⁺ environments to (in effect) discriminate between different RNA architectures suggests that, in some cells, the evolution of functional RNA structures may have been influenced by the presence of putrescine²⁺.

An RNA folding motif: GNRA tetraloop-receptor interactions

Nearly two decades after Westhof and Michel first proposed that RNA tetraloops may interact with distal helices, tetraloop–receptor interactions have been recognized as ubiquitous elements of RNA tertiary structure. The unique architecture of GNRA tetraloops (N=any nucleotide, R=purine) enables interaction with a variety of receptors, e.g., helical minor grooves and asymmetric internal loops. The most common example of the latter is the GAAA tetraloop–11 nt tetraloop receptor motif. Biophysical characterization of this motif provided evidence for the modularity of RNA structure, with applications spanning improved crystallization methods to RNA tectonics. In this review, we identify and compare types of GNRA tetraloop–receptor interactions. Then we explore the abundance of structural, kinetic, and thermodynamic information on the frequently occurring and most widely studied GAAA tetraloop–11 nt receptor motif. Studies of this interaction have revealed powerful paradigms for structural assembly of RNA, as well as providing new insights into the roles of cations, transition states and protein chaperones in RNA folding pathways. However, further research will clearly be necessary to characterize other tetraloop–receptor and long-range tertiary binding interactions in detail – an important milestone in the quantitative prediction of free energy landscapes for RNA folding.

An in vitro-selected RNA receptor for the GAAC loop: modular receptor for non-GNRA-type tetraloop

Although artificial RNA motifs that can functionally replace the GNRA/receptor interaction, a class of RNA-RNA interacting motifs, were isolated from RNA libraries and used to generate designer RNA structures, receptors for non-GNRA tetraloops have not been found in nature or selected from RNA libraries. In this study, we report successful isolation of a receptor motif interacting with GAAC, a non-GNRA tetraloop, from randomized sequences embedded in a catalytic RNA. Biochemical characterization of the GAAC/receptor interacting motif within three structural contexts showed its binding affinity, selectivity and structural autonomy. The motif has binding affinity comparable with that of a GNRA/receptor, selectivity orthogonal to GNRA/receptors and structural autonomy even in a large RNA context. These features would be advantageous for usage of the motif as a building block for designer RNAs. The isolated motif can also be used as a query sequence to search for unidentified naturally occurring GANC receptor motifs.

Database proton NMR chemical shifts for RNA signal assignment and validation

The Biological Magnetic Resonance Data Bank contains NMR chemical shift depositions for 132 RNAs and RNA-containing complexes. We have analyzed the (1)H NMR chemical shifts reported for non-exchangeable protons of residues that reside within A-form helical regions of these RNAs. The analysis focused on the central base pair within a stretch of three adjacent base pairs (BP triplets), and included both Watson-Crick (WC; G:C, A:U) and G:U wobble pairs. Chemical shift values were included for all 4(3) possible WC-BP triplets, as well as 137 additional triplets that contain one or more G:U wobbles. Sequence-dependent chemical shift correlations were identified, including correlations involving terminating base pairs within the triplets and canonical and non-canonical structures adjacent to the BP triplets (i.e. bulges, loops, WC and non-WC BPs), despite the fact that the NMR data were obtained under different conditions of pH, buffer, ionic strength, and temperature. A computer program (RNAShifts) was developed that enables convenient comparison of RNA (1)H NMR assignments with database predictions, which should facilitate future signal assignment/validation efforts and enable rapid identification of non-canonical RNA structures and RNA-ligand/protein interaction sites.

The structural stabilization of the κ three-way junction by Mg(II) represents the first step in the folding of a group II intron

Folding of group II introns is characterized by a first slow compaction of domain 1 (D1) followed by the rapid docking of other domains to this scaffold. D1 compaction initiates in a small subregion encompassing the κ and ζ elements. These two tertiary elements are also the major interaction sites with domain 5 to form the catalytic core. Here, we provide the first characterization of the structure adopted at an early folding step and show that the folding control element can be narrowed down to the three-way junction with the κ motif. In our nuclear magnetic resonance studies of this substructure derived from the yeast mitochondrial group II intron Sc.ai5γ, we show that a high affinity Mg(II) ion stabilizes the κ element and enables coaxial stacking between helices d′ and d′′, favoring a rigid duplex across the three-way junction. The κ-element folds into a stable GAAA-tetraloop motif and engages in A-minor interactions with helix d′. The addition of cobalt(III)hexammine reveals three distinct binding sites. The Mg(II)-promoted structural rearrangement and rigidification of the D1 core can be identified as the first micro-step of D1 folding.

Spectroscopic probes of RNA structure and dynamics

A single technique is insufficient to characterize the properties of an RNA molecule so this Guide provides advice and suggestions for use of several spectroscopic methods applied to RNA. It begins with a discussion of design features to synthesize a suitable molecule for study, assuming that the reader is familiar with in vitro enzymatic and chemical synthesis of RNA. With the RNA in hand, the application of UV melting studies to characterize its folding/unfolding transitions is presented, followed by instructions for fluorescence experiments to augment the UV data. Since RNAs are flexible molecules, it is important to characterize their mobility, and so time-resolved fluorescence data are analyzed. A brief exposition of the power of NMR spectroscopy to identify ion-binding sites is provided. RNA examples are described for each method to give a perspective on what can and cannot be learned.

The role of counterion valence and size in GAAA tetraloop-receptor docking/undocking kinetics

For RNA to fold into compact, ordered structures, it must overcome electrostatic repulsion between negatively charged phosphate groups by counterion recruitment. A physical understanding of the counterion-assisted folding process requires addressing how cations kinetically and thermodynamically control the folding equilibrium for each tertiary interaction in a full-length RNA. In this work, single-molecule FRET (fluorescence resonance energy transfer) techniques are exploited to isolate and explore the cation-concentration-dependent kinetics for formation of a ubiquitous RNA tertiary interaction, that is, the docking/undocking of a GAAA tetraloop with its 11-nt receptor. Rate constants for docking (k(dock)) and undocking (k(undock)) are obtained as a function of cation concentration, size, and valence, specifically for the series Na(+), K(+), Mg(2+), Ca(2+), Co(NH(3))(6)(3+), and spermidine(3+). Increasing cation concentration acceleratesk(dock)dramatically but achieves only a slight decrease in k(undock). These results can be kinetically modeled using parallel cation-dependent and cation-independent docking pathways, which allows for isolation of the folding kinetics from the interaction energetics of the cations with the undocked and docked states, respectively. This analysis reveals a preferential interaction of the cations with the transition state and docked state as compared to the undocked RNA, with the ion-RNA interaction strength growing with cation valence. However, the corresponding number of cations that are taken up by the RNA upon folding decreases with charge density of the cation. The only exception to these behaviors is spermidine(3+), whose weaker influence on the docking equilibria with respect to Co(NH(3))(6)(3+) can be ascribed to steric effects preventing complete neutralization of the RNA phosphate groups.