Structural principles from large RNAs

Temperature-dependent elasticity of DNA, RNA, and hybrid double helices

Nucleic acid double helices in their DNA, RNA, and DNA-RNA hybrid form play a fundamental role in biology and are main building blocks of artificial nanostructures, but how their properties depend on temperature remains poorly understood. Here, we report thermal dependence of dynamic bending persistence length, twist rigidity, stretch modulus, and twist-stretch coupling for DNA, RNA, and hybrid duplexes between 7°C and 47°C. The results are based on all-atom molecular dynamics simulations using different force field parameterizations. We first demonstrate that unrestrained molecular dynamics can reproduce experimentally known mechanical properties of the duplexes at room temperature. Beyond experimentally known features, we also infer the twist rigidity and twist-stretch coupling of the hybrid duplex. As for the temperature dependence, we found that increasing temperature softens all the duplexes with respect to bending, twisting, and stretching. The relative decrease of the stretch moduli is 0.003-0.004/°C, similar for all the duplex variants despite their very different stretching stiffness, whereas RNA twist stiffness decreases by 0.003/°C, and smaller values are found for the other elastic moduli. The twist-stretch couplings are nearly unaffected by temperature. The stretching, bending, and twisting stiffness all include an important entropic component. Relation of our results to the two-state model of DNA flexibility is discussed. Our work provides temperature-dependent elasticity of nucleic acid duplexes at the microsecond scale relevant for initial stages of protein binding.

Temperature-Dependent Twist of Double-Stranded RNA Probed by Magnetic Tweezer Experiments and Molecular Dynamics Simulations

RNA plays critical roles in the transmission and regulation of genetic information and is increasingly used in biomedical and biotechnological applications. Functional RNAs contain extended double-stranded regions, and the structure of double-stranded RNA (dsRNA) has been revealed at high resolution. However, the dependence of the properties of the RNA double helix on environmental effects, notably temperature, is still poorly understood. Here, we use single-molecule magnetic tweezer measurements to determine the dependence of the dsRNA twist on temperature. We find that dsRNA unwinds with increasing temperature, even more than DNA, with ΔTwRNA = -14.4 ± 0.7°/(°C·kbp), compared to ΔTwDNA = -11.0 ± 1.2°/(°C·kbp). All-atom molecular dynamics (MD) simulations using a range of nucleic acid force fields, ion parameters, and water models correctly predict that dsRNA unwinds with rising temperature but significantly underestimate the magnitude of the effect. These MD data, together with additional MD simulations involving DNA and DNA-RNA hybrid duplexes, reveal a linear correlation between the twist temperature decrease and the helical rise, in line with DNA but at variance with RNA experimental data. We speculate that this discrepancy might be caused by some unknown bias in the RNA force fields tested or by as yet undiscovered transient alternative structures in the RNA duplex. Our results provide a baseline to model more complex RNA assemblies and to test and develop new parametrizations for RNA simulations. They may also inspire physical models of the temperature-dependent dsRNA structure.

Context-dependence of T-loop Mediated Long-range RNA Tertiary Interactions

The architecture and folding of complex RNAs is governed by a limited set of highly recurrent structural motifs that form long-range tertiary interactions. One of these motifs is the T-loop, which was first identified in tRNA but is broadly distributed across biological RNAs. While the T-loop has been examined in detail in different biological contexts, the various receptors that it interacts with are not as well defined. In this study, we use a cell-based genetic screen in concert with bioinformatic analysis to examine three different, but related, T-loop receptor motifs found in the flavin mononucleotide (FMN) and cobalamin (Cbl) riboswitches. As a host for different T-loop receptors, we employed the env8 class-II Cbl riboswitch, an RNA that uses two T-loop motifs for both folding and supporting the ligand binding pocket. A set of libraries was created in which select nucleotides that participate in the T-loop/T-loop receptor (TL/TLR) interaction were fully randomized. Library members were screened for their ability to support Cbl-dependent expression of a reporter gene. While T-loops appear to be variable in sequence, we find that the functional sequence space is more restricted in the Cbl riboswitch, suggesting that TL/TLR interactions are context dependent. Our data reveal clear sequence signatures for the different types of receptor motifs that align with phylogenic analysis of these motifs in the FMN and Cbl riboswitches. Finally, our data suggest the functional contribution of various nucleobase-mediated long-range interactions within the riboswitch subclass of TL/TLR interactions that are distinct from those found in other RNAs.

Riboswitch-mediated regulation of riboflavin biosynthesis genes in prokaryotes

Prokaryotic organisms frequently use riboswitches to quantify intracellular metabolite concentration via high-affinity metabolite receptors. Riboswitches possess a metabolite-sensing system that controls gene regulation in a cis-acting fashion at the initiation of transcriptional/translational level by binding with a specific metabolite and controlling various biochemical pathways. Riboswitch binds with flavin mononucleotide (FMN), a phosphorylated form of riboflavin and controls gene expression involved in riboflavin biosynthesis and transport pathway. The first step of the riboflavin biosynthesis pathway is initiated by the conversion of guanine nucleotide triphosphate (GTP), which is an intermediate of the purine biosynthesis pathway. An alternative pentose phosphate pathway of riboflavin biosynthesis includes the enzymatic conversion of ribulose-5-phosphate into 3, 4 dihydroxy-2-butanone-4-phosphates by DHBP synthase. The product of ribAB interferes with both GTP cyclohydrolase II as well as DHBP synthase activities, which catalyze the cleavage of GTP and converts DHBP Ribu5P in the initial steps of both riboflavin biosynthesis branches. Riboswitches are located in the 5' untranslated region (5' UTR) of messenger RNAs and contain an aptamer domain (highly conserved in sequence) where metabolite binding leads to a conformational change in an aptamer domain, which modulate the regulation of gene expression located on bacterial mRNA. In this review, we focus on how riboswitch regulates the riboflavin biosynthesis pathway in Bacillus subtilis and Lactobacillus plantarum.

MetBP: A Software Tool for Detection of Interaction between Metal Ion-RNA Base Pairs

MOTIVATION

The role of metals in shaping and functioning of RNA is a well established fact and the understanding of that through the analysis of structural data has biological relevance. Often metal ions bind to one or more atoms of the nucleobase of an RNA. This fact becomes more interesting when such bases form a base pair with any other base. Furthermore, when metal ions bind to any residue of an RNA, the secondary structural features of the residue (helix, loop, unpaired etc) are also biologically important. The available metal binding related software tools cannot address such type specific queries.

RESULTS

To fill this limitation, we have designed a software tool, called MetBP, that meets the goal. This tool is a stand-alone command line based tool and has no dependency on the other existing software. It accepts a structure file in mmCIF or PDB format and computes the base pairs and thereafter reports all metals that bind to one or more nucleotides that form pairs with another. It reports binding distance, angles along with base pair stability. It also reports several other important aspects, e.g. secondary structure of the residue in the RNA. MetBP can be used as a generalized metal binding site detection tool for Proteins and DNA as well.

AVAILABILITY

https://github.com/computational-biology/metbp.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Phase separation of DNA: From past to present

Phase separation of biological molecules, such as nucleic acids and proteins, has garnered widespread attention across many fields in recent years. For instance, liquid-liquid phase separation has been implicated not only in membraneless intracellular organization but also in many biochemical processes, including transcription, translation, and cellular signaling. Here, we present a historical background of biological phase separation and survey current work on nuclear organization and its connection to DNA phase separation from the perspective of DNA sequence, structure, and genomic context.

Nanosized Particles Assembled by a Recombinant Virus Protein Are Able to Encapsulate Negatively Charged Molecules and Structured RNA

RNA-based molecules have recently become hot candidates to be developed into therapeutic agents. However, successful applications of RNA-based therapeutics might require suitable carriers to protect the RNA from enzymatic degradation by ubiquitous RNases in vivo. Because of their better biocompatibility and biodegradability, protein-based nanoparticles are considered to be alternatives to their synthetic polymer-based counterparts for drug delivery. Hepatitis C virus (HCV) core protein has been suggested to be able to self-assemble into nucleocapsid-like particles in vitro. In this study, the genomic RNA-binding domain of HCV core protein consisting of 116 amino acids (p116) was overexpressed with E. coli for investigation. The recombinant p116 was able to assemble into particles with an average diameter of approximately 27 nm, as visualized by electron microscopy and atomic force microscopy. Measurements with fluorescence spectroscopy, flow cytometry, and fluorescence quenching indicated that the p116-assembled nanoparticles were able to encapsulate small anionic molecules and structured RNA. This study demonstrates methods that exploit the self-assembly nature of a virus-derived protein for nanoparticle production. This study also suggests that the virus-derived protein-assembled particles could possibly be developed into potential carriers for anionic molecular drugs and structured RNA-based therapeutics.

Lighting Up DNA with the Environment-Sensitive Bright Adenine Analogue qAN4

The fluorescent adenine analogue qAN4 was recently shown to possess promising photophysical properties, including a high brightness as a monomer. Here we report the synthesis of the phosphoramidite of qAN4 and its successful incorporation into DNA oligonucleotides using standard solid-phase synthesis. Circular dichroism and thermal melting studies indicate that the qAN4-modification has a stabilizing effect on the B-form of DNA. Moreover, qAN4 base-pairs selectively with thymine with mismatch penalties similar to those of mismatches of adenine. The low energy absorption band of qAN4 inside DNA has its peak around 358 nm and the emission in duplex DNA is partly quenched and blue-shifted (ca. 410 nm), compared to the monomeric form. The spectral properties of the fluorophore also show sensitivity to pH; a property that may find biological applications. Quantum yields in single-stranded DNA range from 1-29 % and in duplex DNA from 1-7 %. In combination with the absorptive properties, this gives an average brightness inside duplex DNA of 275 M^-1  cm^-1 , more than five times higher than the most used environment-sensitive fluorescent base analogue, 2-aminopurine. Finally, we show that qAN4 can be used to advantage as a donor for interbase FRET applications in combination with adenine analogue qA_nitro as an acceptor.

Probing G-quadruplex topologies and recognition concurrently in real time and 3D using a dual-app nucleoside probe

Comprehensive understanding of structure and recognition properties of regulatory nucleic acid elements in real time and atomic level is highly important to devise efficient therapeutic strategies. Here, we report the establishment of an innovative biophysical platform using a dual-app nucleoside analog, which serves as a common probe to detect and correlate different GQ structures and ligand binding under equilibrium conditions and in 3D by fluorescence and X-ray crystallography techniques. The probe (SedU) is composed of a microenvironment-sensitive fluorophore and an excellent anomalous X-ray scatterer (Se), which is assembled by attaching a selenophene ring at 5-position of 2'-deoxyuridine. SedU incorporated into the loop region of human telomeric DNA repeat fluorescently distinguished subtle differences in GQ topologies and enabled quantify ligand binding to different topologies. Importantly, anomalous X-ray dispersion signal from Se could be used to determine the structure of GQs. As the probe is minimally perturbing, a direct comparison of fluorescence data and crystal structures provided structural insights on how the probe senses different GQ conformations without affecting the native fold. Taken together, our dual-app probe represents a new class of tool that opens up new experimental strategies to concurrently investigate nucleic acid structure and recognition in real time and 3D.

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.

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

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

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

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

Large-scale frequent stem pattern mining in RNA families

Functionally similar non-coding RNAs are expected to be similar in certain regions of their secondary structures. These similar regions are called common structure motifs, and are structurally conserved throughout evolution to maintain their functional roles. Common structure motif identification is one of the critical tasks in RNA secondary structure analysis. Nevertheless, current approaches suffer several limitations, and/or do not scale with both structure size and the number of input secondary structures. In this work, we present a method to transform the conserved base pair stems into transaction items and apply frequent itemset mining to identify common structure motifs existing in a majority of input structures. Our experimental results on telomerase and ribosomal RNA secondary structures report frequent stem patterns that are of biological significance. Moreover, the algorithms utilized in our method are scalable and frequent stem patterns can be identified efficiently among many large structures.

Fluorescent nucleobase analogues for base-base FRET in nucleic acids: synthesis, photophysics and applications

Förster resonance energy transfer (FRET) between a donor nucleobase analogue and an acceptor nucleobase analogue, base–base FRET, works as a spectroscopic ruler and protractor. With their firm stacking and ability to replace the natural nucleic acid bases inside the base-stack, base analogue donor and acceptor molecules complement external fluorophores like the Cy-, Alexa- and ATTO-dyes and enable detailed investigations of structure and dynamics of nucleic acid containing systems. The first base–base FRET pair, tC^O–tC_nitro, has recently been complemented with among others the adenine analogue FRET pair, qAN1–qA_nitro, increasing the flexibility of the methodology. Here we present the design, synthesis, photophysical characterization and use of such base analogues. They enable a higher control of the FRET orientation factor, κ², have a different distance window of opportunity than external fluorophores, and, thus, have the potential to facilitate better structure resolution. Netropsin DNA binding and the B-to-Z-DNA transition are examples of structure investigations that recently have been performed using base–base FRET and that are described here. Base–base FRET has been around for less than a decade, only in 2017 expanded beyond one FRET pair, and represents a highly promising structure and dynamics methodology for the field of nucleic acids. Here we bring up its advantages as well as disadvantages and touch upon potential future applications.

Fluorescent RNA cytosine analogue - an internal probe for detailed structure and dynamics investigations

The bright fluorescent cytosine analogue tCO stands out among fluorescent bases due to its virtually unquenched fluorescence emission in duplex DNA. However, like most reported base analogues, it has not been thoroughly characterized in RNA. We here report on the first synthesis and RNA-incorporation of tCO, and characterize its base-mimicking and fluorescence properties in RNA. As in DNA, we find a high quantum yield inside RNA duplexes (<ΦF> = 0.22) that is virtually unaffected by the neighbouring bases (ΦF = 0.20-0.25), resulting in an average brightness of 1900 M-1 cm-1. The average fluorescence lifetime in RNA duplexes is 4.3 ns and generally two lifetimes are required to fit the exponential decays. Fluorescence properties in ssRNA are defined by a small increase in average quantum yield (<ΦF > = 0.24) compared to dsRNA, with a broader distribution (ΦF = 0.17-0.34) and slightly shorter average lifetimes. Using circular dichroism, we find that the tCO-modified RNA duplexes form regular A-form helices and in UV-melting experiments the stability of the duplexes is only slightly higher than that of the corresponding natural RNA (<ΔT m> = + 2.3 °C). These properties make tCO a highly interesting fluorescent RNA base analogue for detailed FRET-based structural measurements, as a bright internal label in microscopy, and for fluorescence anisotropy measurements of RNA dynamics.

CHSalign: A Web Server That Builds upon Junction-Explorer and RNAJAG for Pairwise Alignment of RNA Secondary Structures with Coaxial Helical Stacking

RNA junctions are important structural elements of RNA molecules. They are formed when three or more helices come together in three-dimensional space. Recent studies have focused on the annotation and prediction of coaxial helical stacking (CHS) motifs within junctions. Here we exploit such predictions to develop an efficient alignment tool to handle RNA secondary structures with CHS motifs. Specifically, we build upon our Junction-Explorer software for predicting coaxial stacking and RNAJAG for modelling junction topologies as tree graphs to incorporate constrained tree matching and dynamic programming algorithms into a new method, called CHSalign, for aligning the secondary structures of RNA molecules containing CHS motifs. Thus, CHSalign is intended to be an efficient alignment tool for RNAs containing similar junctions. Experimental results based on thousands of alignments demonstrate that CHSalign can align two RNA secondary structures containing CHS motifs more accurately than other RNA secondary structure alignment tools. CHSalign yields a high score when aligning two RNA secondary structures with similar CHS motifs or helical arrangement patterns, and a low score otherwise. This new method has been implemented in a web server, and the program is also made freely available, at http://bioinformatics.njit.edu/CHSalign/.

Rsite2: an efficient computational method to predict the functional sites of noncoding RNAs

Noncoding RNAs (ncRNAs) represent a big class of important RNA molecules. Given the large number of ncRNAs, identifying their functional sites is becoming one of the most important topics in the post-genomic era, but available computational methods are limited. For the above purpose, we previously presented a tertiary structure based method, Rsite, which first calculates the distance metrics defined in Methods with the tertiary structure of an ncRNA and then identifies the nucleotides located within the extreme points in the distance curve as the functional sites of the given ncRNA. However, the application of Rsite is largely limited because of limited RNA tertiary structures. Here we present a secondary structure based computational method, Rsite2, based on the observation that the secondary structure based nucleotide distance is strongly positively correlated with that derived from tertiary structure. This makes it reasonable to replace tertiary structure with secondary structure, which is much easier to obtain and process. Moreover, we applied Rsite2 to three ncRNAs (tRNA (Lys), Diels-Alder ribozyme, and RNase P) and a list of human mitochondria transcripts. The results show that Rsite2 works well with nearly equivalent accuracy as Rsite but is much more feasible and efficient. Finally, a web-server, the source codes, and the dataset of Rsite2 are available at http://www.cuialb.cn/rsite2.

Computational study of stability of an H-H-type pseudoknot motif

Motifs in RNA tertiary structures are important to their structural organizations and biological functions. Here we consider an H-H-type pseudoknot (HHpk) motif that consists of two hairpins connected by a junction loop and with kissing interactions between the two hairpin loops. Such a tertiary structural motif is recurrently found in RNA tertiary structures, but is difficult to predict computationally. So it is important to understand the mechanism of its formation and stability. Here we investigate the stability of the HHpk tertiary structure by using an all-atom molecular dynamics simulation. The results indicate that the HHpk tertiary structure is stable. However, it is found that this stability is not due to the helix-helix packing, as is usually expected, but is maintained by the combined action of the kissing hairpin loops and junctions, although the former plays the main role. Stable HHpk motifs may form structural platforms for the molecules to realize their biological functions. These results are useful for understanding the construction principle of RNA tertiary structures and structure prediction.

Synthesis and Substrate Evaluation of (E)-5-[(3-Selenophene-2-Carboxamido)Prop-1-en-1-yl]-Uridine-5'-O-Triphosphate for RNA Polymerase

Design, synthesis and T7 RNA polymerase substrate evaluation of (E)-5-[(3-selenophene-2-carboxamido)prop-1-en-1-yl]-uridine-5'-O-triphosphate is reported. The title compound is shown to be a good substrate for RNA polymerase by RNA labeling through in vitro transcription. pTRI-plasmid DNA with β-actin gene sequence (∼300 base pairs) with T7 promoter was used as a template for the in vitro transcription. Transcribed product is characterized for incorporation by gel assay and for integrity, full length and size by bioanalyzer. The title compound will be very useful in biophysical techniques to obtain information on dynamics and recognition properties in real time as well as 3D structure of nucleic acids.

MultiSETTER: web server for multiple RNA structure comparison

Background

Understanding the architecture and function of RNA molecules requires methods for comparing and analyzing their tertiary and quaternary structures. While structural superposition of short RNAs is achievable in a reasonable time, large structures represent much bigger challenge. Therefore, we have developed a fast and accurate algorithm for RNA pairwise structure superposition called SETTER and implemented it in the SETTER web server. However, though biological relationships can be inferred by a pairwise structure alignment, key features preserved by evolution can be identified only from a multiple structure alignment. Thus, we extended the SETTER algorithm to the alignment of multiple RNA structures and developed the MultiSETTER algorithm.

Results

In this paper, we present the updated version of the SETTER web server that implements a user friendly interface to the MultiSETTER algorithm. The server accepts RNA structures either as the list of PDB IDs or as user-defined PDB files. After the superposition is computed, structures are visualized in 3D and several reports and statistics are generated.

Conclusion

To the best of our knowledge, the MultiSETTER web server is the first publicly available tool for a multiple RNA structure alignment. The MultiSETTER server offers the visual inspection of an alignment in 3D space which may reveal structural and functional relationships not captured by other multiple alignment methods based either on a sequence or on secondary structure motifs.

Multiple 3D RNA Structure Superposition Using Neighbor Joining

Recent advances in RNA research and the steady growth of available RNA structures call for bioinformatics methods for handling and analyzing RNA structural data. Recently, we introduced SETTER-a fast and accurate method for RNA pairwise structure alignment. In this paper, we describe MultiSETTER, SETTER extension for multiple RNA structure alignment. MultiSETTER combines SETTER's decomposition of RNA structures into non-overlapping structural subunits with the multiple sequence alignment algorithm ClustalW adapted for the structure alignment. The accuracy of MultiSETTER was assessed by the automatic classification of RNA structures and its comparison to SCOR annotations. In addition, MultiSETTER classification was also compared to multiple sequence alignment-based and secondary structure alignment-based classifications provided by LocARNA and RNADistance tools, respectively. MultiSETTER precompiled Windows libraries, as well as the C++ source code, are freely available from http://siret.cz/multisetter.

Small-angle X-ray scattering: a bridge between RNA secondary structures and three-dimensional topological structures

Whereas the structures of small to medium-sized well folded RNA molecules often can be determined by either X-ray crystallography or NMR spectroscopy, obtaining structural information for large RNAs using experimental, computational, or combined approaches remains a major interest and challenge. RNA is very sensitive to small-angle X-ray scattering (SAXS) due to high electron density along phosphate-sugar backbones, whose scattering contribution dominates SAXS intensity. For this reason, SAXS is particularly useful in obtaining global RNA structural information that outlines backbone topologies and, therefore, molecular envelopes. Such information is extremely valuable in bridging the gap between the secondary structures and three-dimensional topological structures of RNA molecules, particularly those that have proven difficult to study using other structure-determination methods. Here we review published results of RNA topological structures derived from SAXS data or in combination with other experimental data, as well as details on RNA sample preparation for SAXS experiments.

Structure-specific nucleic acid recognition by L-motifs and their diverse roles in expression and regulation of the genome

The high-mobility group (HMG) domain containing proteins regulate transcription, DNA replication and recombination. They adopt L-shaped folds and are structure-specific DNA binding motifs. Here, I define the L-motif super-family that consists of DNA-binding HMG-box proteins and the L-motif of the histone mRNA binding domain of stem-loop binding protein (SLBP). The SLBP L-motif and HMG-box domains adopt similar L-shaped folds with three α-helices and two or three small hydrophobic cores that stabilize the overall fold, but have very different and distinct modes of nucleic acid recognition. A comparison of the structure, dynamics, protein-protein and nucleic acid interactions, and regulation by PTMs of the SLBP and the HMG-box L-motifs reveals the versatile and diverse modes by which L-motifs utilize their surfaces for structure-specific recognition of nucleic acids to regulate gene expression.

3dRNAscore: a distance and torsion angle dependent evaluation function of 3D RNA structures

Model evaluation is a necessary step for better prediction and design of 3D RNA structures. For proteins, this has been widely studied and the knowledge-based statistical potential has been proved to be one of effective ways to solve this problem. Currently, a few knowledge-based statistical potentials have also been proposed to evaluate predicted models of RNA tertiary structures. The benchmark tests showed that they can identify the native structures effectively but further improvements are needed to identify near-native structures and those with non-canonical base pairs. Here, we present a novel knowledge-based potential, 3dRNAscore, which combines distance-dependent and dihedral-dependent energies. The benchmarks on different testing datasets all show that 3dRNAscore are more efficient than existing evaluation methods in recognizing native state from a pool of near-native states of RNAs as well as in ranking near-native states of RNA models.

Paradigms of ribosome synthesis: Lessons learned from ribosomal proteins

The proteome in all cells is manufactured via the intricate process of translation by multimolecular factories called ribosomes. Nevertheless, these ribonucleoprotein particles, the largest of their kind, also have an elaborate assembly line of their own. Groundbreaking discoveries that bacterial ribosomal subunits can be self-assembled in vitro jumpstarted studies on how ribosomes are constructed. Until recently, ribosome assembly has been investigated almost entirely in vitro with bacterial small subunits under equilibrium conditions. In light of high-resolution ribosome structures and a more sophisticated toolkit, the past decade has been defined by a burst of kinetic studies in vitro and, importantly, also a shift to examining ribosome maturation in living cells, especially in eukaryotes. In this review, we summarize the principles governing ribosome assembly that emerged from studies focusing on ribosomal proteins and their interactions with rRNA. Understanding these paradigms has taken center stage, given the linkage between anomalous ribosome biogenesis and proliferative disorders.

Molecular crowding overcomes the destabilizing effects of mutations in a bacterial ribozyme

The native structure of the Azoarcus group I ribozyme is stabilized by the cooperative formation of tertiary interactions between double helical domains. Thus, even single mutations that break this network of tertiary interactions reduce ribozyme activity in physiological Mg(2+) concentrations. Here, we report that molecular crowding comparable to that in the cell compensates for destabilizing mutations in the Azoarcus ribozyme. Small angle X-ray scattering, native polyacrylamide gel electrophoresis and activity assays were used to compare folding free energies in dilute and crowded solutions containing 18% PEG1000. Crowder molecules allowed the wild-type and mutant ribozymes to fold at similarly low Mg(2+) concentrations and stabilized the active structure of the mutant ribozymes under physiological conditions. This compensation helps explains why ribozyme mutations are often less deleterious in the cell than in the test tube. Nevertheless, crowding did not rescue the high fraction of folded but less active structures formed by double and triple mutants. We conclude that crowding broadens the fitness landscape by stabilizing compact RNA structures without improving the specificity of self-assembly.

Heavy atom containing fluorescent ribonucleoside analog probe for the fluorescence detection of RNA-ligand binding

Although numerous biophysical tools have provided effective systems to study nucleic acids, our current knowledge on how RNA structure complements its function is limited. Therefore, development of robust tools to study the structure–function relationship of RNA is highly desired. Toward this endeavor, we have developed a new ribonucleoside analog, based on a (selenophen-2-yl)pyrimidine core, which could serve as a fluorescence probe to study the function of RNA in real time and as an anomalous scattering label (selenium atom) for the phase determination in X-ray crystallography. The fluorescent selenophene-modified uridine analog is minimally perturbing and exhibits probe-like properties such as sensitivity to microenvironment and conformation changes. Utilizing these properties and amicability of the corresponding ribonucleotide analog to enzymatic incorporation, we have synthesized a fluorescent bacterial ribosomal decoding site (A-site) RNA construct and have developed a fluorescence binding assay to effectively monitor the binding of aminoglycoside antibiotics to the A-site. Our results demonstrate that this simple approach of building a dual probe could provide new avenues to study the structure–function relationship of not only nucleic acids, but also other biomacromolecules.

Topological constraints are major determinants of tRNA tertiary structure and dynamics and provide basis for tertiary folding cooperativity

Recent studies have shown that basic steric and connectivity constraints encoded at the secondary structure level are key determinants of 3D structure and dynamics in simple two-way RNA junctions. However, the role of these topological constraints in higher order RNA junctions remains poorly understood. Here, we use a specialized coarse-grained molecular dynamics model to directly probe the thermodynamic contributions of topological constraints in defining the 3D architecture and dynamics of transfer RNA (tRNA). Topological constraints alone restrict tRNA's allowed conformational space by over an order of magnitude and strongly discriminate against formation of non-native tertiary contacts, providing a sequence independent source of folding specificity. Topological constraints also give rise to long-range correlations between the relative orientation of tRNA's helices, which in turn provides a mechanism for encoding thermodynamic cooperativity between distinct tertiary interactions. These aspects of topological constraints make it such that only several tertiary interactions are needed to confine tRNA to its native global structure and specify functionally important 3D dynamics. We further show that topological constraints are conserved across tRNA's different naturally occurring secondary structures. Taken together, our results emphasize the central role of secondary-structure-encoded topological constraints in defining RNA 3D structure, dynamics and folding.

Viral RNAs are unusually compact

A majority of viruses are composed of long single-stranded genomic RNA molecules encapsulated by protein shells with diameters of just a few tens of nanometers. We examine the extent to which these viral RNAs have evolved to be physically compact molecules to facilitate encapsulation. Measurements of equal-length viral, non-viral, coding and non-coding RNAs show viral RNAs to have among the smallest sizes in solution, i.e., the highest gel-electrophoretic mobilities and the smallest hydrodynamic radii. Using graph-theoretical analyses we demonstrate that their sizes correlate with the compactness of branching patterns in predicted secondary structure ensembles. The density of branching is determined by the number and relative positions of 3-helix junctions, and is highly sensitive to the presence of rare higher-order junctions with 4 or more helices. Compact branching arises from a preponderance of base pairing between nucleotides close to each other in the primary sequence. The density of branching represents a degree of freedom optimized by viral RNA genomes in response to the evolutionary pressure to be packaged reliably. Several families of viruses are analyzed to delineate the effects of capsid geometry, size and charge stabilization on the selective pressure for RNA compactness. Compact branching has important implications for RNA folding and viral assembly.

A hierarchical model for assembly of eukaryotic 60S ribosomal subunit domains

Despite having structures for eukaryotic large ribosomal subunits, it has been unclear how these ribonucleoprotein complexes are constructed in living cells. Here, Gamalinda et al. used proteomics to assay changes in the levels of ribosomal protein and assembly factors in preribosomes upon depletion of ribosomal proteins that function in 60S assembly. The results show that the eukaryotic 60S ribosomal subunits are assembled in a hierarchical fashion. This study further reveals striking differences and similarities between bacterial and eukaryotic large ribosomal subunits, providing insights into the evolution of RNA–protein particles.

Despite having high-resolution structures for eukaryotic large ribosomal subunits, it remained unclear how these ribonucleoprotein complexes are constructed in living cells. Nevertheless, knowing where ribosomal proteins interact with ribosomal RNA (rRNA) provides a strategic platform to investigate the connection between spatial and temporal aspects of 60S subunit biogenesis. We previously found that the function of individual yeast large subunit ribosomal proteins (RPLs) in precursor rRNA (pre-rRNA) processing correlates with their location in the structure of mature 60S subunits. This observation suggested that there is an order by which 60S subunits are formed. To test this model, we used proteomic approaches to assay changes in the levels of ribosomal proteins and assembly factors in preribosomes when RPLs functioning in early, middle, and late steps of pre-60S assembly are depleted. Our results demonstrate that structural domains of eukaryotic 60S ribosomal subunits are formed in a hierarchical fashion. Assembly begins at the convex solvent side, followed by the polypeptide exit tunnel, the intersubunit side, and finally the central protuberance. This model provides an initial paradigm for the sequential assembly of eukaryotic 60S subunits. Our results reveal striking differences and similarities between assembly of bacterial and eukaryotic large ribosomal subunits, providing insights into how these RNA–protein particles evolved.

An unusual topological structure of the HIV-1 Rev response element

Nuclear export of unspliced and singly spliced viral mRNA is a critical step in the HIV life cycle. The structural basis by which the virus selects its own mRNA among more abundant host cellular RNAs for export has been a mystery for more than 25 years. Here, we describe an unusual topological structure that the virus uses to recognize its own mRNA. The viral Rev response element (RRE) adopts an "A"-like structure in which the two legs constitute two tracks of binding sites for the viral Rev protein and position the two primary known Rev-binding sites ~55 Å apart, matching the distance between the two RNA-binding motifs in the Rev dimer. Both the legs of the "A" and the separation between them are required for optimal RRE function. This structure accounts for the specificity of Rev for the RRE and thus the specific recognition of the viral RNA.

Coarse-grained simulations of RNA and DNA duplexes

Although RNAs play many cellular functions, little is known about the dynamics and thermodynamics of these molecules. In principle, all-atom molecular dynamics simulations can investigate these issues, but with current computer facilities, these simulations have been limited to small RNAs and to short times. HiRe-RNA, a recently proposed high-resolution coarse-grained RNA that captures many geometric details such as base pairing and stacking, is able to fold RNA molecules to near-native structures in a short computational time. So far, it had been applied to simple hairpins, and here we present its application to duplexes of a couple dozen nucleotides and show how with replica exchange molecular dynamics (REMD) we can easily predict the correct double helix from a completely random configuration and study the dissociation curve. To show the versatility of our model, we present an application to a double stranded DNA molecule as well. A reconstruction algorithm allows us to obtain full atom structures from the coarse-grained model. Through atomistic molecular dynamics (MD), we can compare the dynamics starting from a representative structure of a low temperature replica or from the experimental structure, and show how the two are statistically identical, highlighting the validity of a coarse-grained approach for structured RNAs and DNAs.

Structure and function of the T-loop structural motif in noncoding RNAs

The T-loop is a frequently occurring five-nucleotide motif found in the structure of noncoding RNAs where it is commonly assumed to play an important role in stabilizing the tertiary RNA structure by facilitating long-range interactions between different regions of the molecule. T-loops were first identified in tRNA(Phe) and a formal consensus sequence for this motif was formulated and later revised based on analyses of the crystal structures of prokaryotic ribosomal RNAs and RNase P and the corresponding primary sequence of their orthologues. In the past decade, several new structures of large RNA molecules have been added to the RCSB Protein Data Bank, including the eukaryotic ribosome, a self-splicing group II intron, numerous synthetases in complex with their cognate transfer RNAs (tRNAs), transfer-messenger RNA (tmRNA) in complex with SmpB, several riboswitches, and a complex of bacterial RNase P bound to its tRNA substrate. In this review, the search for T-loops is extended to these new RNA molecules based on the previously established structure-based criteria. The review highlights and discusses the function and additional roles of T-loops in four broad categories of RNA molecules, namely tRNAs, ribosomal RNAs (rRNAs), P RNAs, and RNA genetic elements. Additionally, the potential application for T-loops as interaction modules is also discussed.

RNA-DNA hybrid origami: folding of a long RNA single strand into complex nanostructures using short DNA helper strands

Quick folding of a long RNA strand using short DNA staple strands (at a 1 : 1 ratio) into various pre-designed nanostructures in high yields has been demonstrated.

Novel cis-acting element within the capsid-coding region enhances flavivirus viral-RNA replication by regulating genome cyclization

cis-Acting elements in the viral genome RNA (vRNA) are essential for the translation, replication, and/or encapsidation of RNA viruses. In this study, a novel conserved cis-acting element was identified in the capsid-coding region of mosquito-borne flavivirus. The downstream of 5' cyclization sequence (5'CS) pseudoknot (DCS-PK) element has a three-stem pseudoknot structure, as demonstrated by structure prediction and biochemical analysis. Using dengue virus as a model, we show that DCS-PK enhances vRNA replication and that its function depends on its secondary structure and specific primary sequence. Mutagenesis revealed that the highly conserved stem 1 and loop 2, which are involved in potential loop-helix interactions, are crucial for DCS-PK function. A predicted loop 1-stem 3 base triple interaction is important for the structural stability and function of DCS-PK. Moreover, the function of DCS-PK depends on its position relative to the 5'CS, and the presence of DCS-PK facilitates the formation of 5'-3' RNA complexes. Taken together, our results reveal that the cis-acting element DCS-PK enhances vRNA replication by regulating genome cyclization, and DCS-PK might interplay with other cis-acting elements to form a functional vRNA cyclization domain, thus playing critical roles during the flavivirus life cycle and evolution.

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.

The GA-minor submotif as a case study of RNA modularity, prediction, and design

Complex natural RNAs such as the ribosome, group I and group II introns, and RNase P exemplify the fact that three-dimensional (3D) RNA structures are highly modular and hierarchical in nature. Tertiary RNA folding typically takes advantage of a rather limited set of recurrent structural motifs that are responsible for controlling bends or stacks between adjacent helices. Herein, the GA minor and related structural motifs are presented as a case study to highlight several structural and folding principles, to gain further insight into the structural evolution of naturally occurring RNAs, as well as to assist the rational design of artificial RNAs.

FRETmatrix: a general methodology for the simulation and analysis of FRET in nucleic acids

Förster resonance energy transfer (FRET) is a technique commonly used to unravel the structure and conformational changes of biomolecules being vital for all living organisms. Typically, FRET is performed using dyes attached externally to nucleic acids through a linker that complicates quantitative interpretation of experiments because of dye diffusion and reorientation. Here, we report a versatile, general methodology for the simulation and analysis of FRET in nucleic acids, and demonstrate its particular power for modelling FRET between probes possessing limited diffusional and rotational freedom, such as our recently developed nucleobase analogue FRET pairs (base–base FRET). These probes are positioned inside the DNA/RNA structures as a replacement for one of the natural bases, thus, providing unique control of their position and orientation and the advantage of reporting from inside sites of interest. In demonstration studies, not requiring molecular dynamics modelling, we obtain previously inaccessible insight into the orientation and nanosecond dynamics of the bases inside double-stranded DNA, and we reconstruct high resolution 3D structures of kinked DNA. The reported methodology is accompanied by a freely available software package, FRETmatrix, for the design and analysis of FRET in nucleic acid containing systems.

Advances in quantitative FRET-based methods for studying nucleic acids

Förster resonance energy transfer (FRET) is a powerful tool for monitoring molecular distances and interactions at the nanoscale level. The strong dependence of transfer efficiency on probe separation makes FRET perfectly suited for "on/off" experiments. To use FRET to obtain quantitative distances and three-dimensional structures, however, is more challenging. This review summarises recent studies and technological advances that have improved FRET as a quantitative molecular ruler in nucleic acid systems, both at the ensemble and at the single-molecule levels.

Fingerprinting noncanonical and tertiary RNA structures by differential SHAPE reactivity

Many RNA structures are composed of simple secondary structure elements linked by a few critical tertiary interactions. SHAPE chemistry has made interrogation of RNA dynamics at single-nucleotide resolution straightforward. However, de novo identification of nucleotides involved in tertiary interactions remains a challenge. Here we show that nucleotides that form noncanonical or tertiary contacts can be detected by comparing information obtained using two SHAPE reagents, N-methylisatoic anhydride (NMIA) and 1-methyl-6-nitroisatoic anhydride (1M6). Nucleotides that react preferentially with NMIA exhibit slow local nucleotide dynamics and usually adopt the less common C2'-endo ribose conformation. Experiments and first-principles calculations show that 1M6 reacts preferentially with nucleotides in which one face of the nucleobase allows an unhindered stacking interaction with the reagent. Differential SHAPE reactivities were used to detect noncanonical and tertiary interactions in four RNAs with diverse structures and to identify preformed noncanonical interactions in partially folded RNAs. Differential SHAPE reactivity analysis will enable experimentally concise, large-scale identification of tertiary structure elements and ligand binding sites in complex RNAs and in diverse biological environments.

Structure and mechanism of purine-binding riboswitches

A riboswitch is a non-protein coding sequence capable of directly binding a small molecule effector without the assistance of accessory proteins to regulate expression of the mRNA in which it is embedded. Currently, over 20 different classes of riboswitches have been validated in bacteria with the promise of many more to come, making them an important means of regulating the genome in the bacterial kingdom. Strikingly, half of the known riboswitches recognize effector compounds that contain a purine or related moiety. In the last decade, significant progress has been made to determine how riboswitches specifically recognize these compounds against the background of many other similar cellular metabolites and transduce this signal into a regulatory response. Of the known riboswitches, the purine family containing guanine, adenine and 2'-deoxyguanosine-binding classes are the most extensively studied, serving as a simple and useful paradigm for understanding how these regulatory RNAs function. This review provides a comprehensive summary of the current state of knowledge regarding the structure and mechanism of these riboswitches, as well as insights into how they might be exploited as therapeutic targets and novel biosensors.

Conformational features of topologically classified RNA secondary structures

Background

Current RNA secondary structure prediction approaches predict prevalent pseudoknots such as the H-pseudoknot and kissing hairpin. The number of possible structures increases drastically when more complex pseudoknots are considered, thus leading to computational limitations. On the other hand, the enormous population of possible structures means not all of them appear in real RNA molecules. Therefore, it is of interest to understand how many of them really exist and the reasons for their preferred existence over the others, as any new findings revealed by this study might enhance the capability of future structure prediction algorithms for more accurate prediction of complex pseudoknots.

Methodology/Principal Findings

A novel algorithm was devised to estimate the exact number of structural possibilities for a pseudoknot constructed with a specified number of base pair stems. Then, topological classification was applied to classify RNA pseudoknotted structures from data in the RNA STRAND database. By showing the vast possibilities and the real population, it is clear that most of these plausible complex pseudoknots are not observed. Moreover, from these classified motifs that exist in nature, some features were identified for further investigation. It was found that some features are related to helical stacking. Other features are still left open to discover underlying tertiary interactions.

Conclusions

Results from topological classification suggest that complex pseudoknots are usually some well-known motifs that are themselves complex or the interaction results of some special motifs. Heuristics can be proposed to predict the essential parts of these complex motifs, even if the required thermodynamic parameters are currently unknown.

SETTER: web server for RNA structure comparison

The recent discoveries of regulatory non-coding RNAs changed our view of RNA as a simple information transfer molecule. Understanding the architecture and function of active RNA molecules requires methods for comparing and analyzing their 3D structures. While structural alignment of short RNAs is achievable in a reasonable amount of time, large structures represent much bigger challenge. Here, we present the SETTER web server for the RNA structure pairwise comparison utilizing the SETTER (SEcondary sTructure-based TERtiary Structure Similarity Algorithm) algorithm. The SETTER method divides an RNA structure into the set of non-overlapping structural elements called generalized secondary structure units (GSSUs). The SETTER algorithm scales as O(n²) with the size of a GSSUs and as O(n) with the number of GSSUs in the structure. This scaling gives SETTER its high speed as the average size of the GSSU remains constant irrespective of the size of the structure. However, the favorable speed of the algorithm does not compromise its accuracy. The SETTER web server together with the stand-alone implementation of the SETTER algorithm are freely accessible at http://siret.cz/setter.

Cooperative tertiary interaction network guides RNA folding

Noncoding RNAs form unique 3D structures, which perform many regulatory functions. To understand how RNAs fold uniquely despite a small number of tertiary interaction motifs, we mutated the major tertiary interactions in a group I ribozyme by single-base substitutions. The resulting perturbations to the folding energy landscape were measured using SAXS, ribozyme activity, hydroxyl radical footprinting, and native PAGE. Double- and triple-mutant cycles show that most tertiary interactions have a small effect on the stability of the native state. Instead, the formation of core and peripheral structural motifs is cooperatively linked in near-native folding intermediates, and this cooperativity depends on the native helix orientation. The emergence of a cooperative interaction network at an early stage of folding suppresses nonnative structures and guides the search for the native state. We suggest that cooperativity in noncoding RNAs arose from natural selection of architectures conducive to forming a unique, stable fold.