RNA structure: the long and the short of it

FebRNA: An automated fragment-ensemble-based model for building RNA 3D structures

Knowledge of RNA three-dimensional (3D) structures is critical to understanding the important biological functions of RNAs. Although various structure prediction models have been developed, the high-accuracy predictions of RNA 3D structures are still limited to the RNAs with short lengths or with simple topology. In this work, we proposed a new model, namely FebRNA, for building RNA 3D structures through fragment assembly based on coarse-grained (CG) fragment ensembles. Specifically, FebRNA is composed of four processes: establishing the library of different types of non-redundant CG fragment ensembles regardless of the sequences, building CG 3D structure ensemble through fragment assembly, identifying top-scored CG structures through a specific CG scoring function, and rebuilding the all-atom structures from the top-scored CG ones. Extensive examination against different types of RNA structures indicates that FebRNA consistently gives the reliable predictions on RNA 3D structures, including pseudoknots, three-way junctions, four-way and five-way junctions, and RNAs in the RNA-Puzzles. FebRNA is available on the Web site: https://github.com/Tan-group/FebRNA.

Kinetics of protein-assisted nucleic acid interconversion monitored by transient time resolved fluorescence in microfluidic droplets

Interconversions between nucleic acid structures play an important role in transcriptional and translational regulation and also in repair and recombination. These interconversions are frequently promoted by nucleic acid chaperone proteins. To monitor their kinetics, Förster resonance energy transfer (FRET) is widely exploited using ensemble fluorescence intensity measurements in pre-steady-state stopped-flow experiments. Such experiments only provide a weighted average of the emission of all species in solution and consume large quantities of materials. Herein, we lift these limitations by combining time-resolved fluorescence (TRF) with droplet microfluidics (DmF). We validate the innovative TRF-DmF approach by investigating the well characterized annealing of the HIV-1 (+)/(–) Primer Binding Sequences (PBS) promoted by a HIV-1 nucleocapsid peptide. Upon rapid mixing of the FRET-labelled (–)PBS with its complementary (+)PBS sequence inside microdroplets, the TRF-DmF set-up enables resolving the time evolution of sub-populations of reacting species and reveals an early intermediate with a ∼50 ps donor fluorescence lifetime never identified so far. TRF-DmF also favorably compares with single molecule experiments, as it offers an accurate control of concentrations with no upper limit, no need to graft one partner on a surface and no photobleaching issues.

Evaluating RNA Structural Flexibility: Viruses Lead the Way

Our understanding of RNA structure has lagged behind that of proteins and most other biological polymers, largely because of its ability to adopt multiple, and often very different, functional conformations within a single molecule. Flexibility and multifunctionality appear to be its hallmarks. Conventional biochemical and biophysical techniques all have limitations in solving RNA structure and to address this in recent years we have seen the emergence of a wide diversity of techniques applied to RNA structural analysis and an accompanying appreciation of its ubiquity and versatility. Viral RNA is a particularly productive area to study in that this economy of function within a single molecule admirably suits the minimalist lifestyle of viruses. Here, we review the major techniques that are being used to elucidate RNA conformational flexibility and exemplify how the structure and function are, as in all biology, tightly linked.

Engineering Crystal Packing in RNA Structures I: Past and Future Strategies for Engineering RNA Packing in Crystals

X-ray crystallography remains a powerful method to gain atomistic insights into the catalytic and regulatory functions of RNA molecules. However, the technique requires the preparation of diffraction-quality crystals. This is often a resource- and time-consuming venture because RNA crystallization is hindered by the conformational heterogeneity of RNA, as well as the limited opportunities for stereospecific intermolecular interactions between RNA molecules. The limited success at crystallization explains in part the smaller number of RNA-only structures in the Protein Data Bank. Several approaches have been developed to aid the formation of well-ordered RNA crystals. The majority of these are construct-engineering techniques that aim to introduce crystal contacts to favor the formation of well-diffracting crystals. A typical example is the insertion of tetraloop-tetraloop receptor pairs into non-essential RNA segments to promote intermolecular association. Other methods of promoting crystallization involve chaperones and crystallization-friendly molecules that increase RNA stability and improve crystal packing. In this review, we discuss the various techniques that have been successfully used to facilitate crystal packing of RNA molecules, recent advances in construct engineering, and directions for future research in this vital aspect of RNA crystallography.

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.

The effects of long non-coding ribonucleic acids on various cellular components in rheumatoid arthritis

RA is a chronic, autoimmune-mediated inflammatory pathology. Long non-coding RNAs (lncRNAs) are a novel group of non-coding RNAs with a length of >200 nucleotides. There are reports emerging that suggest that lncRNAs participate in establishing and sustaining autoimmune diseases, including RA. In this review article, we highlight the functions of lncRNAs in different cell types in RA. Our review indicates that lncRNAs affect various cellular components and are novel candidates that could constitute promising targets for the diagnosis and treatment of RA.

The interplay between molecular flexibility and RNA chemical probing reactivities analyzed at the nucleotide level via an extensive molecular dynamics study

Determination of the tridimensional structure of ribonucleic acid molecules is fundamental for understanding their function in the cell. A common method to investigate RNA structures of large molecules is the use of chemical probes such as SHAPE (2^'-hydroxyl acylation analyzed by primer extension) reagents, DMS (dimethyl sulfate) and CMCT (1-cyclohexyl-3-(2-morpholinoethyl) carbodiimide metho-p-toluene sulfate), the reaction of which is dependent on the local structural properties of each nucleotide. In order to understand the interplay between local flexibility, sugar pucker, canonical pairing and chemical reactivity of the probes, we performed all-atom molecular dynamics simulations on a set of RNA molecules for which both tridimensional structure and chemical probing data are available and we analyzed the correlations between geometrical parameters and the chemical reactivity. Our study confirms that SHAPE reactivity is guided by the local flexibility of the different chemical moieties but suggests that a combination of multiple parameters is needed to better understand the implications of the reactivity at the molecular level. This is also the case for DMS and CMCT for which the reactivity appears to be more complex than commonly accepted.

Tapping the RNA world for therapeutics

A recent revolution in RNA biology has led to the identification of new RNA classes with unanticipated functions, new types of RNA modifications, an unexpected multiplicity of alternative transcripts and widespread transcription of extragenic regions. This development in basic RNA biology has spawned a corresponding revolution in RNA-based strategies to generate new types of therapeutics. Here, I review RNA-based drug design and discuss barriers to broader applications and possible ways to overcome them. Because they target nucleic acids rather than proteins, RNA-based drugs promise to greatly extend the domain of 'druggable' targets beyond what can be achieved with small molecules and biologics.

IsRNA: An Iterative Simulated Reference State Approach to Modeling Correlated Interactions in RNA Folding

Coarse-grained RNA folding models promise great potential for RNA structure prediction. A key component in a coarse-grained folding model is the force field. One of the challenges in the coarse-grained force field calculation is how to treat the correlation between the different degrees of freedoms. Here, we describe a new approach (IsRNA) to extract the correlated energy functions from the known structures. Through iterative molecular dynamics simulations, we build the correlation effects into the reference states, from which we extract the energy functions. The validity of IsRNA is supported by the close agreement between the simulated Boltzmann-like probability distributions for all the structure parameters and those observed from the experimentally determined structures. The correlated energy functions derived here may provide a new tool for RNA 3D structure prediction.

Functional Studies of DNA-Protein Interactions Using FRET Techniques

Protein-DNA interactions underpin life and play key roles in all cellular processes and functions including DNA transcription, packaging, replication, and repair. Identifying and examining the nature of these interactions is therefore a crucial prerequisite to understand the molecular basis of how these fundamental processes take place. The application of fluorescence techniques and in particular fluorescence resonance energy transfer (FRET) to provide structural and kinetic information has experienced a stunning growth during the past decade. This has been mostly promoted by new advances in the preparation of dye-labeled nucleic acids and proteins and in optical sensitivity, where its implementation at the level of individual molecules has opened a new biophysical frontier. Nowadays, the application of FRET-based techniques to the analysis of protein-DNA interactions spans from the classical steady-state and time-resolved methods averaging over large ensembles to the analysis of distances, conformational changes, and enzymatic reactions in individual protein-DNA complexes. This chapter introduces the practical aspects of applying these methods for the study of protein-DNA interactions.

Embracing proteins: structural themes in aptamer-protein complexes

Understanding the structural rules that govern specific, high-affinity binding characteristic of aptamer-protein interactions is important in view of the increasing use of aptamers across many applications. From the modest number of 16 aptamer-protein structures currently available, trends are emerging. The flexible phosphodiester backbone allows folding into precise three-dimensional structures using known nucleic acid motifs as scaffolds that orient specific functional groups for target recognition. Still, completely novel motifs essential for structure and function are found in modified aptamers with diversity-enhancing side chains. Aptamers and antibodies, two classes of macromolecules used as affinity reagents with entirely different backbones and composition, recognize protein epitopes of similar size and with comparably high shape complementarity.

Phylogenetic Reconstruction of the Calosphaeriales and Togniniales Using Five Genes and Predicted RNA Secondary Structures of ITS, and Flabellascus tenuirostris gen. et sp. nov

The Calosphaeriales is revisited with new collection data, living cultures, morphological studies of ascoma centrum, secondary structures of the internal transcribed spacer (ITS) rDNA and phylogeny based on novel DNA sequences of five nuclear ribosomal and protein-coding loci. Morphological features, molecular evidence and information from predicted RNA secondary structures of ITS converged upon robust phylogenies of the Calosphaeriales and Togniniales. The current concept of the Calosphaeriales includes the Calosphaeriaceae and Pleurostomataceae encompassing five monophyletic genera, Calosphaeria, Flabellascus gen. nov., Jattaea, Pleurostoma and Togniniella, strongly supported by Bayesian and Maximum Likelihood methods. The structural elements of ITS1 form characteristic patterns that are phylogenetically conserved, corroborate observations based on morphology and have a high predictive value at the generic level. Three major clades containing 44 species of Phaeoacremonium were recovered in the closely related Togniniales based on ITS, actin and β-tubulin sequences. They are newly characterized by sexual and RNA structural characters and ecology. This approach is a first step towards understanding of the molecular systematics of Phaeoacremonium and possibly its new classification. In the Calosphaeriales, Jattaea aphanospora sp. nov. and J. ribicola sp. nov. are introduced, Calosphaeria taediosa is combined in Jattaea and epitypified. The sexual morph of Phaeoacremonium cinereum was encountered for the first time on decaying wood and obtained in vitro. In order to achieve a single nomenclature, the genera of asexual morphs linked with the Calosphaeriales are transferred to synonymy of their sexual morphs following the principle of priority, i.e. Calosphaeriophora to Calosphaeria, Phaeocrella to Togniniella and Pleurostomophora to Pleurostoma. Three new combinations are proposed, i.e. Pleurostoma ochraceum comb. nov., P. repens comb. nov. and P. richardsiae comb. nov. The morphology-based key is provided to facilitate identification of genera accepted in the Calosphaeriales.

Dwell-Time Distribution, Long Pausing and Arrest of Single-Ribosome Translation through the mRNA Duplex

Proteins in the cell are synthesized by a ribosome translating the genetic information encoded on the single-stranded messenger RNA (mRNA). It has been shown that the ribosome can also translate through the duplex region of the mRNA by unwinding the duplex. Here, based on our proposed model of the ribosome translation through the mRNA duplex we study theoretically the distribution of dwell times of the ribosome translation through the mRNA duplex under the effect of a pulling force externally applied to the ends of the mRNA to unzip the duplex. We provide quantitative explanations of the available single molecule experimental data on the distribution of dwell times with both short and long durations, on rescuing of the long paused ribosomes by raising the pulling force to unzip the duplex, on translational arrests induced by the mRNA duplex and Shine-Dalgarno(SD)-like sequence in the mRNA. The functional consequences of the pauses or arrests caused by the mRNA duplex and the SD sequence are discussed and compared with those obtained from other types of pausing, such as those induced by "hungry" codons or interactions of specific sequences in the nascent chain with the ribosomal exit tunnel.

An Efficient Minimum Free Energy Structure-Based Search Method for Riboswitch Identification Based on Inverse RNA Folding

Riboswitches are RNA genetic control elements that were originally discovered in bacteria and provide a unique mechanism of gene regulation. They work without the participation of proteins and are believed to represent ancient regulatory systems in the evolutionary timescale. One of the biggest challenges in riboswitch research is to find additional eukaryotic riboswitches since more than 20 riboswitch classes have been found in prokaryotes but only one class has been found in eukaryotes. Moreover, this single known class of eukaryotic riboswitch, namely the TPP riboswitch class, has been found in bacteria, archaea, fungi and plants but not in animals. The few examples of eukaryotic riboswitches were identified using sequence-based bioinformatics search methods such as a combination of BLAST and pattern matching techniques that incorporate base-pairing considerations. None of these approaches perform energy minimization structure predictions. There is a clear motivation to develop new bioinformatics methods, aside of the ongoing advances in covariance models, that will sample the sequence search space more flexibly using structural guidance while retaining the computational efficiency of sequence-based methods. We present a new energy minimization approach that transforms structure-based search into a sequence-based search, thereby enabling the utilization of well established sequence-based search utilities such as BLAST and FASTA. The transformation to sequence space is obtained by using an extended inverse RNA folding problem solver with sequence and structure constraints, available within RNAfbinv. Examples in applying the new method are presented for the purine and preQ1 riboswitches. The method is described in detail along with its findings in prokaryotes. Potential uses in finding novel eukaryotic riboswitches and optimizing pre-designed synthetic riboswitches based on ligand simulations are discussed. The method components are freely available for use.

Ab initio RNA folding

RNA molecules are essential cellular machines performing a wide variety of functions for which a specific three-dimensional structure is required. Over the last several years, the experimental determination of RNA structures through x-ray crystallography and NMR seems to have reached a plateau in the number of structures resolved each year, but as more and more RNA sequences are being discovered, the need for structure prediction tools to complement experimental data is strong. Theoretical approaches to RNA folding have been developed since the late nineties, when the first algorithms for secondary structure prediction appeared. Over the last 10 years a number of prediction methods for 3D structures have been developed, first based on bioinformatics and data-mining, and more recently based on a coarse-grained physical representation of the systems. In this review we are going to present the challenges of RNA structure prediction and the main ideas behind bioinformatic approaches and physics-based approaches. We will focus on the description of the more recent physics-based phenomenological models and on how they are built to include the specificity of the interactions of RNA bases, whose role is critical in folding. Through examples from different models, we will point out the strengths of physics-based approaches, which are able not only to predict equilibrium structures, but also to investigate dynamical and thermodynamical behavior, and the open challenges to include more key interactions ruling RNA folding.

Coarse-Grained HiRE-RNA Model for ab Initio RNA Folding beyond Simple Molecules, Including Noncanonical and Multiple Base Pairings

HiRE-RNA is a coarse-grained model for RNA structure prediction and the dynamical study of RNA folding. Using a reduced set of particles and detailed interactions accounting for base-pairing and stacking, we show that noncanonical and multiple base interactions are necessary to capture the full physical behavior of complex RNAs. In this paper, we give a full account of the model and present results on the folding, stability, and free energy surfaces of 16 systems with 12 to 76 nucleotides of increasingly complex architectures, ranging from monomers to dimers, using a total of 850 μs of simulation time.

RNA pseudo-knots simulated with a one-bead coarse-grained model

We present a revised version of a Monte Carlo simulation model for RNA molecules that was introduced in a previous communication [O. Taxilaga-Zetina, P. Pliego-Pastrana, and M. D. Carbajal-Tinoco, Phys. Rev. E 81, 041914 (2010)]. The basic model consists of a series of knowledge-based pair potentials that were obtained from the statistical analysis of large RNAs belonging to the Protein Data Bank. These effective interactions are then used to dress a polymeric chain that reproduces relatively simple secondary structures (e.g., small hairpins). In order to describe more complicated three-dimensional structures such as pseudo-knots, here we include orientational information for the interaction between nucleotides forming hydrogen bonds, as in the case of the Watson-Crick base pairs. As a result, the simulated molecules obtained through the modified model are now consistent with their corresponding experimental configurations.

Modeling of RNA nanotubes using molecular dynamics simulation

In this study, we construct novel RNA nanoclusters, RNA nanotubes made of several nanorings up to the size of 20 nm, utilizing the molecular dynamics simulation, and study their structural properties [i.e., the root mean square deviation, the radius of gyration and the radial distribution function (RDF)] in physiological solutions that can be used for drug delivery into the human body. The patterns of energy and temperature variations of the systems are also discussed. Furthermore, we study the concentration of ions around the tube as a function of time at a particular temperature. We have found that when the temperature increases, the number of ions increases within a certain distance of the tube. We report that the number of ions within this distance around the tubes decreases in quenched runs. This indicates that some ions evaporate with decrease in temperature, as has been observed in the case of the nanoring. RDF plots also demonstrate a similar trend with temperature, as was found in the case of RNA nanorings.

Sizing up long non-coding RNAs: do lncRNAs have secondary and tertiary structure?

Long noncoding RNAs (lncRNAs) play a key role in many important areas of epigenetics, stem cell biology, cancer, signaling and brain function. This emerging class of RNAs constitutes a large fraction of the transcriptome, with thousands of new lncRNAs reported each year. The molecular mechanisms of these RNAs are not well understood. Currently, very little structural data exist. We review the available lncRNA sequence and secondary structure data. Since almost no tertiary information is available for lncRNAs, we review crystallographic structures for other RNA systems and discuss the possibilities for lncRNAs in the context of existing constraints.

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.

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.

Voronoia4RNA--a database of atomic packing densities of RNA structures and their complexes

Voronoia4RNA (http://proteinformatics.charite.de/voronoia4rna/) is a structural database storing precalculated atomic volumes, atomic packing densities (PDs) and coordinates of internal cavities for currently 1869 RNAs and RNA-protein complexes. Atomic PDs are a measure for van der Waals interactions. Regions of low PD, containing water-sized internal cavities, refer to local structure flexibility or compressibility. RNA molecules build up the skeleton of large molecular machineries such as ribosomes or form smaller flexible structures such as riboswitches. The wealth of structural data on RNAs and their complexes allows setting up representative data sets and analysis of their structural features. We calculated atomic PDs from atomic volumes determined by the Voronoi cell method and internal cavities analytically by Delaunay triangulation. Reference internal PD values were derived from a non-redundant sub-data set of buried atoms. Comparison of internal PD values shows that RNA is more tightly packed than proteins. Finally, the relation between structure size, resolution and internal PD of the Voronoia4RNA entries is discussed. RNA, protein structures and their complexes can be visualized by the Jmol-based viewer Provi. Variations in PD are depicted by a color code. Internal cavities are represented by their molecular boundaries or schematically as balls.

Hydration sites of unpaired RNA bases: a statistical analysis of the PDB structures

BACKGROUND

Hydration is crucial for RNA structure and function. X-ray crystallography is the most commonly used method to determine RNA structures and hydration and, therefore, statistical surveys are based on crystallographic results, the number of which is quickly increasing.

RESULTS

A statistical analysis of the water molecule distribution in high-resolution X-ray structures of unpaired RNA nucleotides showed that: different bases have the same penchant to be surrounded by water molecules; clusters of water molecules indicate possible hydration sites, which, in some cases, match those of the major and minor grooves of RNA and DNA double helices; complex hydrogen bond networks characterize the solvation of the nucleotides, resulting in a significant rigidity of the base and its surrounding water molecules. Interestingly, the hydration sites around unpaired RNA bases do not match, in general, the positions that are occupied by the second nucleotide when the base-pair is formed.

CONCLUSIONS

The hydration sites around unpaired RNA bases were found. They do not replicate the atom positions of complementary bases in the Watson-Crick pairs.

Novel ligands for a purine riboswitch discovered by RNA-ligand docking

The increasing number of RNA crystal structures enables a structure-based approach to the discovery of new RNA-binding ligands. To develop the poorly explored area of RNA-ligand docking, we have conducted a virtual screening exercise for a purine riboswitch to probe the strengths and weaknesses of RNA-ligand docking. Using a standard protein-ligand docking program with only minor modifications, four new ligands with binding affinities in the micromolar range were identified, including two compounds based on molecular scaffolds not resembling known ligands. RNA-ligand docking performed comparably to protein-ligand docking indicating that this approach is a promising option to explore the wealth of RNA structures for structure-based ligand design.

A quantitative measure of chirality inside nucleic acid databank

We show the capability of a chirality index (Pietropaolo et al., Proteins 2008;70:667-677) to investigate nucleic acid structures because of its high sensitivity to helical conformations. By analyzing selected structures of DNA and RNA, we have found that sequences rich in cytosine and guanine have a tendency to left-handed chirality, in contrast to regions rich in adenine or thymine which show strong negative, right-handed, chirality values. We also analyze RNA structures, where specific loops and hairpin motifs are characterized by a well-defined chirality value. We find that in nucleosome the chirality is exalted, whereas in ribosome it is reduced. Our results illustrate the sensitivity of this descriptor for nucleic acid conformations.

Metal binding motif in the active site of the HDV ribozyme binds divalent and monovalent ions

The hepatitis delta virus (HDV) ribozyme uses both metal ion and nucleobase catalysis in its cleavage mechanism. A reverse G·U wobble was observed in a recent crystal structure of the precleaved state. This unusual base pair positions a Mg(2+) ion to participate in catalysis. Herein, we used molecular dynamics (MD) and X-ray crystallography to characterize the conformation and metal binding characteristics of this base pair in product and precleaved forms. Beginning with a crystal structure of the product form, we observed formation of the reverse G·U wobble during MD trajectories. We also demonstrated that this base pair is compatible with the diffraction data for the product-bound state. During MD trajectories of the product form, Na(+) ions interacted with the reverse G·U wobble in the RNA active site, and a Mg(2+) ion, introduced in certain trajectories, remained bound at this site. Beginning with a crystal structure of the precleaved form, the reverse G·U wobble with bound Mg(2+) remained intact during MD simulations. When we removed Mg(2+) from the starting precleaved structure, Na(+) ions interacted with the reverse G·U wobble. In support of the computational results, we observed competition between Na(+) and Mg(2+) in the precleaved ribozyme crystallographically. Nonlinear Poisson-Boltzmann calculations revealed a negatively charged patch near the reverse G·U wobble. This anionic pocket likely serves to bind metal ions and to help shift the pK(a) of the catalytic nucleobase, C75. Thus, the reverse G·U wobble motif serves to organize two catalytic elements, a metal ion and catalytic nucleobase, within the active site of the HDV ribozyme.

Identifying and searching for conserved RNA localisation signals

RNA localisation is an important mode of delivering proteins to their site of function. Cis-acting signals within the RNAs, which can be thought of as zip-codes, determine the site of localisation. There are few examples of fully characterised RNA signals, but the signals are thought to be defined through a combination of primary, secondary, and tertiary structures. In this chapter, we describe a selection of computational methods for predicting RNA secondary structure, identifying localisation signals, and searching for similar localisation signals on a genome-wide scale. The chapter is aimed at the biologist rather than presenting the details of each of the individual methods.

An RNA G-quadruplex is essential for cap-independent translation initiation in human VEGF IRES

RNA G-quadruplexes located within the 5'-UTR of mRNA are almost always known to be associated with repression of cap-dependent translation. However, in this report we present functional as well as structural evidence that sequence redundancy in a G-rich segment within the 5'-UTR of human VEGF mRNA supports a 'switchable' RNA G-quadruplex structure that is essential for IRES-mediated translation initiation. Additionally, utilization of a specific combination of G-tracts within this segment allows for the conformational switch that implies a tunable regulatory role of the quadruplex structure in translation initiation. A sequence engineered from a functionally handicapped mutant moderately rescued the activity, further indicating the importance of G-quadruplex structure for VEGF IRES-A function. This to our knowledge is the first report of a conformationally flexible RNA G-quadruplex which is essential for IRES-mediated translation initiation.

Role of unsatisfied hydrogen bond acceptors in RNA energetics and specificity

RNA plays essential roles in much of biology. These functions are dictated by structures mediated by hydrogen bonding, stacking, electrostatics, and steric interactions. Roles of unsatisfied hydrogen bond functionalities in these structures are less well understood. Herein, we evaluated the energetic contributions of unsatisfied hydrogen bonding groups by placing chemically modified substituents in select internal positions in RNA helices and conducting thermodynamic studies. We find that unsatisfied carbonyl groups make exceptional contributions to structure formation (approximately 3 kcal/mol in free energy), most likely due to a combination of strain and dehydration effects. Thus, unsatisfied hydrogen bonding groups are likely key determinants in the folding energetics and specificity of many RNA and DNA molecules and may be especially important in tertiary structure interactions.

Coarse-graining RNA nanostructures for molecular dynamics simulations

A series of coarse-grained models have been developed for study of the molecular dynamics of RNA nanostructures. The models in the series have one to three beads per nucleotide and include different amounts of detailed structural information. Such a treatment allows us to reach, for systems of thousands of nucleotides, a time scale of microseconds (i.e. by three orders of magnitude longer than in full atomistic modeling) and thus to enable simulations of large RNA polymers in the context of bionanotechnology. We find that the three-beads-per-nucleotide models, described by a set of just a few universal parameters, are able to describe different RNA conformations and are comparable in structural precision to the models where detailed values of the backbone P-C4' dihedrals taken from a reference structure are included. These findings are discussed in the context of RNA conformation classes.

FRASS: the web-server for RNA structural comparison

Background

The impressive increase of novel RNA structures, during the past few years, demands automated methods for structure comparison. While many algorithms handle only small motifs, few techniques, developed in recent years, (ARTS, DIAL, SARA, SARSA, and LaJolla) are available for the structural comparison of large and intact RNA molecules.

Results

The FRASS web-server represents a RNA chain with its Gauss integrals and allows one to compare structures of RNA chains and to find similar entries in a database derived from the Protein Data Bank. We observed that FRASS scores correlate well with the ARTS and LaJolla similarity scores. Moreover, the-web server can also reproduce satisfactorily the DARTS classification of RNA 3D structures and the classification of the SCOR functions that was obtained by the SARA method.

Conclusions

The FRASS web-server can be easily used to detect relationships among RNA molecules and to scan efficiently the rapidly enlarging structural databases.

The transDSL ligase ribozyme can utilize various forms of modules to clamp its substrate and enzyme units

TransDSL is an RNA ligase ribozyme whose enzyme unit joins two RNA fragments constituting a substrate. The enzyme unit recognizes the substrate by means of two clamp modules. We constructed active variants by replacing the original clamp module with various types of interactions. Such flexible modularity would be advantageous in the application of this ribozyme in nanobiotechnology.

RNA structural motifs that entail hydrogen bonds involving sugar-phosphate backbone atoms of RNA

The growing number of high-resolution crystal structures of large RNA molecules provides much information for understanding the principles of structural organization of these complex molecules. Several in-depth analyses of nucleobase-centered RNA structural motifs and backbone conformations have been published based on this information, including a systematic classification of base pairs by Leontis and Westhof. However, hydrogen bonds involving sugar-phosphate backbone atoms of RNA have not been analyzed systematically until recently, although such hydrogen bonds appear to be common both in local and tertiary interactions. Here we review some backbone structural motifs discussed in the literature and analyze a set of eight high-resolution multi-domain RNA structures. The analyzed RNAs are highly structured: among 5372 nucleotides in this set, 89% are involved in at least one "long-range" RNA-RNA hydrogen bond, i.e., hydrogen bonds between atoms in the same residue or sequential residues are ignored. These long-range hydrogen bonds frequently use backbone atoms as hydrogen bond acceptors, i.e., OP1, OP2, O2', O3', O4', or O5', or as a donor (2'OH). A surprisingly large number of such hydrogen bonds are found, considering that neither single-stranded nor double-stranded regions will contain such hydrogen bonds unless additional interactions with other residues exist. Among 8327 long-range hydrogen bonds found in this set of structures, 2811, or about one-third, are hydrogen bonds entailing RNA backbone atoms; they involve 39% of all nucleotides in the structures. The majority of them (2111) are hydrogen bonds entailing ribose hydroxyl groups, which can be used either as a donor or an acceptor; they constitute 25% of all hydrogen bonds and involve 31% of all nucleotides. The phosphate oxygens OP1 or OP2 are used as hydrogen bond acceptors in 12% of all nucleotides, and the ribose ring oxygen O4' and phosphodiester oxygens O3' and O5' are used in 4%, 4%, and 1% of all nucleotides, respectively. Distributions of geometric parameters and some examples of such hydrogen bonds are presented in this report. A novel motif involving backbone hydrogen bonds, the ribose-phosphate zipper, is also identified.

Genome packaging in influenza A virus

The negative-sense RNA genome of influenza A virus is composed of eight segments, which encode 12 proteins between them. At the final stage of viral assembly, these genomic virion (v)RNAs are incorporated into the virion as it buds from the apical plasma membrane of the cell. Genome segmentation confers evolutionary advantages on the virus, but also poses a problem during virion assembly as at least one copy of each of the eight segments is required to produce a fully infectious virus particle. Historically, arguments have been presented in favour of a specific packaging mechanism that ensures incorporation of a full genome complement, as well as for an alternative model in which segments are chosen at random but packaged in sufficient numbers to ensure that a reasonable proportion of virions are viable. The question has seen a resurgence of interest in recent years leading to a consensus that the vast majority of virions contain no more than eight segments and that a specific mechanism does indeed function to select one copy of each vRNA. This review summarizes work leading to this conclusion. In addition, we describe recent progress in identifying the specific packaging signals and discuss likely mechanisms by which these RNA elements might operate.

Raman crystallography of RNA

Raman crystallography is the application of Raman spectroscopy to single crystals. This technique has been applied to a variety of protein molecules where it has provided unique information about biopolymer folding, substrate binding, and catalysis. Here, we describe the application of Raman crystallography to functional RNA molecules. RNA represents unique opportunities and challenges for Raman crystallography. One issue that confounds studies of RNA is its tendency to adopt multiple non-functional folds. Raman crystallography has the advantage that it isolates a single state of the RNA within the crystal and can evaluate its fold, metal ion binding properties (ligand identity, stoichiometry, and affinity), proton binding properties (identity, stoichiometry, and affinity), and catalytic potential. In particular, base-specific stretches can be identified and then associated with the binding of metal ions and protons. Because measurements are carried out in the hanging drop at ambient, rather than cryo, conditions and because RNA crystals tend to be approximately 70% solvent, RNA dynamics and conformational changes become experimentally accessible. This review focuses on experimental setup and procedures, acquisition and interpretation of Raman data, and determination of physicochemical properties of the RNA. Raman crystallographic and solution biochemical experiments on the HDV RNA enzyme are summarized and found to be in excellent agreement. Remarkably, characterization of the crystalline state has proven to help rather than hinder functional characterization of functional RNA, most likely because the tendency of RNA to fold heterogeneously is limited in a crystalline environment. Future applications of Raman crystallography to RNA are briefly discussed.

Molecular dynamics study of the RNA ring nanostructure: a phenomenon of self-stabilization

We study mechanical and thermodynamic properties of RNA nanostructures focusing on a hexagonal nanoring discussed in Yingling and Shapiro (2007 Nano Lett. 7 2328). We are concerned with the following main issues: (i) the stability of the nanoring versus temperature; (ii) the effect of the environment (solvent, counterions) on its stability; (iii) conformations and dynamics under external force. The process of evaporation of the ions from the ring upon temperature drop has been found, demonstrating a surprising feature-the uptake of ions by the nanoring increases with the temperature. The connection of this behavior to the dielectric constant of water, hydration and structural changes in the nanoring is discussed. Several properties of the nanoring, such as elastic and transport coefficients, have been determined. A measure of the tensile elasticity of the ring against its uniform 2D in-plane compression has been given, as K(eff) < or = 0.01 GPa, which is a much lower value compared to typical values found for soft matter other than RNA.

Viroid intercellular trafficking: RNA motifs, cellular factors and broad impacts

Viroids are noncoding RNAs that infect plants. In order to establish systemic infection, these RNAs must traffic from an initially infected host cell into neighboring cells and ultimately throughout a whole plant. Recent studies have identified structural motifs in a viroid that are required for trafficking, enabling further studies on the mechanisms of their function. Some cellular proteins interact with viroids in vivo and may play a role in viroid trafficking, which can now be directly tested by using a virus-induced gene silencing system that functions efficiently in plant species from which these factors were identified. This review discusses these recent advances, unanswered questions and the use of viroid infection as an highly productive model to elucidate mechanisms of RNA trafficking that is of broad biological significance.

Tertiary motifs revealed in analyses of higher-order RNA junctions

RNA junctions are secondary-structure elements formed when three or more helices come together. They are present in diverse RNA molecules with various fundamental functions in the cell. To better understand the intricate architecture of three-dimensional (3D) RNAs, we analyze currently solved 3D RNA junctions in terms of base-pair interactions and 3D configurations. First, we study base-pair interaction diagrams for solved RNA junctions with 5 to 10 helices and discuss common features. Second, we compare these higher-order junctions to those containing 3 or 4 helices and identify global motif patterns such as coaxial stacking and parallel and perpendicular helical configurations. These analyses show that higher-order junctions organize their helical components in parallel and helical configurations similar to lower-order junctions. Their sub-junctions also resemble local helical configurations found in three- and four-way junctions and are stabilized by similar long-range interaction preferences such as A-minor interactions. Furthermore, loop regions within junctions are high in adenine but low in cytosine, and in agreement with previous studies, we suggest that coaxial stacking between helices likely forms when the common single-stranded loop is small in size; however, other factors such as stacking interactions involving noncanonical base pairs and proteins can greatly determine or disrupt coaxial stacking. Finally, we introduce the ribo-base interactions: when combined with the along-groove packing motif, these ribo-base interactions form novel motifs involved in perpendicular helix-helix interactions. Overall, these analyses suggest recurrent tertiary motifs that stabilize junction architecture, pack helices, and help form helical configurations that occur as sub-elements of larger junction networks. The frequent occurrence of similar helical motifs suggest nature's finite and perhaps limited repertoire of RNA helical conformation preferences. More generally, studies of RNA junctions and tertiary building blocks can ultimately help in the difficult task of RNA 3D structure prediction.

Hydrogen bonding and packing density are factors most strongly connected to limiting sites of high flexibility in the 16S rRNA in the 30S ribosome

BACKGROUND

Conformational flexibility in structured RNA frequently is critical to function. The 30S ribosomal subunit exists in different conformations in different functional states due to changes in the central part of the 16S rRNA. We are interested in evaluating the factors that might be responsible for restricting flexibility to specific parts of the 16S rRNA using biochemical data obtained from the 30S subunit in solution. This problem was approached taking advantage of the observation that there must be a high degree of conformational flexibility at sites where UV photocrosslinking occurs and a lack of flexibility inhibits photoreactivity at many other sites that are otherwise suitable for reaction.

RESULTS

We used 30S x-ray structures to quantify the properties of the nucleotide pairs at UV- and UVA-s4U-induced photocrosslinking sites in 16S rRNA and compared these to the properties of many hundreds of additional sites that have suitable geometry but do not undergo photocrosslinking. Five factors that might affect RNA flexibility were investigated - RNA interactions with ribosomal proteins, interactions with Mg2+ ions, the presence of long-range A minor motif interactions, hydrogen bonding and the count of neighboring heavy atoms around the center of each nucleobase to estimate the neighbor packing density. The two factors that are very different in the unreactive inflexible pairs compared to the reactive ones are the average number of hydrogen bonds and the average value for the number of neighboring atoms. In both cases, these factors are greater for the unreactive nucleotide pairs at a statistically very significant level.

CONCLUSION

The greater extent of hydrogen bonding and neighbor atom density in the unreactive nucleotide pairs is consistent with reduced flexibility at a majority of the unreactive sites. The reactive photocrosslinking sites are clustered in the 30S subunit and this indicates nonuniform patterns of hydrogen bonding and packing density in the 16S rRNA tertiary structure. Because this analysis addresses inter-nucleotide distances and geometry between nucleotides distant in the primary sequence, the results indicate regional and global flexibility of the rRNA.

Analysis of four-way junctions in RNA structures

RNA secondary structures can be divided into helical regions composed of canonical Watson-Crick and related base pairs, as well as single-stranded regions such as hairpin loops, internal loops, and junctions. These elements function as building blocks in the design of diverse RNA molecules with various fundamental functions in the cell. To better understand the intricate architecture of three-dimensional (3D) RNAs, we analyze existing RNA four-way junctions in terms of base-pair interactions and 3D configurations. Specifically, we identify nine broad junction families according to coaxial stacking patterns and helical configurations. We find that helices within junctions tend to arrange in roughly parallel and perpendicular patterns and stabilize their conformations using common tertiary motifs such as coaxial stacking, loop-helix interaction, and helix packing interaction. Our analysis also reveals a number of highly conserved base-pair interaction patterns and novel tertiary motifs such as A-minor-coaxial stacking combinations and sarcin/ricin motif variants. Such analyses of RNA building blocks can ultimately help in the difficult task of RNA 3D structure prediction.

T box transcription antitermination riboswitch: influence of nucleotide sequence and orientation on tRNA binding by the antiterminator element

Many bacteria utilize riboswitch transcription regulation to monitor and appropriately respond to cellular levels of important metabolites or effector molecules. The T box transcription antitermination riboswitch responds to cognate uncharged tRNA by specifically stabilizing an antiterminator element in the 5'-untranslated mRNA leader region and precluding formation of a thermodynamically more stable terminator element. Stabilization occurs when the tRNA acceptor end base pairs with the first four nucleotides in the seven nucleotide bulge of the highly conserved antiterminator element. The significance of the conservation of the antiterminator bulge nucleotides that do not base pair with the tRNA is unknown, but they are required for optimal function. In vitro selection was used to determine if the isolated antiterminator bulge context alone dictates the mode in which the tRNA acceptor end binds the bulge nucleotides. No sequence conservation beyond complementarity was observed and the location was not constrained to the first four bases of the bulge. The results indicate that formation of a structure that recognizes the tRNA acceptor end in isolation is not the determinant driving force for the high phylogenetic sequence conservation observed within the antiterminator bulge. Additional factors or T box leader features more likely influenced the phylogenetic sequence conservation.

Functional studies of DNA-protein interactions using FRET techniques

Protein-DNA interactions underpin life and play key roles in all cellular processes and functions including DNA transcription, packaging, replication, and repair. Identifying and examining the nature of these interactions is therefore a crucial prerequisite to understand the molecular basis of how these fundamental processes take place. The application of fluorescence techniques and in particular fluorescence resonance energy transfer (FRET) to provide structural and kinetic information has experienced a stunning growth during the past decade. This has been mostly promoted by new advances in the preparation of dye-labeled nucleic acids and proteins and in optical sensitivity, where its implementation at the level of individual molecules has opened a new biophysical frontier. Nowadays, the application of FRET-based techniques to the analysis of protein-DNA interactions spans from the classical steady-state and time-resolved methods averaging over large ensembles to the analysis of distances, conformational changes, and enzymatic reactions in individual Protein-DNA complexes. This chapter introduces the practical aspects of applying these methods for the study of Protein-DNA interactions.

Automated motif extraction and classification in RNA tertiary structures

We used a novel graph-based approach to extract RNA tertiary motifs. We cataloged them all and clustered them using an innovative graph similarity measure. We applied our method to three widely studied structures: Haloarcula marismortui 50S (H.m 50S), Escherichia coli 50S (E. coli 50S), and Thermus thermophilus 16S (T.th 16S) RNAs. We identified 10 known motifs without any prior knowledge of their shapes or positions. We additionally identified four putative new motifs.

The UA_handle: a versatile submotif in stable RNA architectures

Stable RNAs are modular and hierarchical 3D architectures taking advantage of recurrent structural motifs to form extensive non-covalent tertiary interactions. Sequence and atomic structure analysis has revealed a novel submotif involving a minimal set of five nucleotides, termed the UA_handle motif (5'XU/AN(n)X3'). It consists of a U:A Watson-Crick: Hoogsteen trans base pair stacked over a classic Watson-Crick base pair, and a bulge of one or more nucleotides that can act as a handle for making different types of long-range interactions. This motif is one of the most versatile building blocks identified in stable RNAs. It enters into the composition of numerous recurrent motifs of greater structural complexity such as the T-loop, the 11-nt receptor, the UAA/GAN and the G-ribo motifs. Several structural principles pertaining to RNA motifs are derived from our analysis. A limited set of basic submotifs can account for the formation of most structural motifs uncovered in ribosomal and stable RNAs. Structural motifs can act as structural scaffoldings and be functionally and topologically equivalent despite sequence and structural differences. The sequence network resulting from the structural relationships shared by these RNA motifs can be used as a proto-language for assisting prediction and rational design of RNA tertiary structures.