The lonepair triloop: a new motif in RNA structure

Features and Functions of the A-Minor Motif, the Most Common Motif in RNA Structure

A-minor motifs are RNA tertiary structure motifs that generally involve a canonical base pair and an adenine base forming hydrogen bonds with the minor groove of the base pair. Such motifs are among the most common tertiary interactions in known RNA structures, comparable in number with the non-canonical base pairs. They are often found in functionally important regions of non-coding RNAs and, in particular, play a central role in protein synthesis. Here, we review local variations of the A-minor geometry and discuss difficulties associated with their annotation, as well as various structural contexts and common A-minor co-motifs, and diverse functions of A-minors in various processes in a living cell.

Conservation in the Iron Responsive Element Family

Iron responsive elements (IREs) are mRNA stem-loop targets for translational control by the two iron regulatory proteins IRP1 and IRP2. They are found in the untranslated regions (UTRs) of genes that code for proteins involved in iron metabolism. There are ten “classic” IRE types that define the conserved secondary and tertiary structure elements necessary for proper IRP binding, and there are 83 published “IRE-like” sequences, most of which depart from the established IRE model. Here are structurally-guided discussions regarding the essential features of an IRE and what is important for IRE family membership.

A functional genetic screen reveals sequence preferences within a key tertiary interaction in cobalamin riboswitches required for ligand selectivity

Riboswitches are structured mRNA sequences that regulate gene expression by directly binding intracellular metabolites. Generating the appropriate regulatory response requires the RNA rapidly and stably acquire higher-order structure to form the binding pocket, bind the appropriate effector molecule and undergo a structural transition to inform the expression machinery. These requirements place riboswitches under strong kinetic constraints, likely restricting the sequence space accessible by recurrent structural modules such as the kink turn and the T-loop. Class-II cobalamin riboswitches contain two T-loop modules: one directing global folding of the RNA and another buttressing the ligand binding pocket. While the T-loop module directing folding is highly conserved, the T-loop associated with binding is substantially less so, with no clear consensus sequence. To further understand the functional role of the binding-associated module, a functional genetic screen of a library of riboswitches with the T-loop and its interacting nucleotides was used to build an experimental phylogeny comprised of sequences that possess a wide range of cobalamin-dependent regulatory activity. Our results reveal conservation patterns of the T-loop and its interaction with the binding core that allow for rapid tertiary structure formation and demonstrate its importance for generating strong ligand-dependent repression of mRNA expression.

A conserved RNA structural motif for organizing topology within picornaviral internal ribosome entry sites

Picornaviral IRES elements are essential for initiating the cap-independent viral translation. However, three-dimensional structures of these elements remain elusive. Here, we report a 2.84-Å resolution crystal structure of hepatitis A virus IRES domain V (dV) in complex with a synthetic antibody fragment-a crystallization chaperone. The RNA adopts a three-way junction structure, topologically organized by an adenine-rich stem-loop motif. Despite no obvious sequence homology, the dV architecture shows a striking similarity to a circularly permuted form of encephalomyocarditis virus J-K domain, suggesting a conserved strategy for organizing the domain architecture. Recurrence of the motif led us to use homology modeling tools to compute a 3-dimensional structure of the corresponding domain of foot-and-mouth disease virus, revealing an analogous domain organizing motif. The topological conservation observed among these IRESs and other viral domains implicates a structured three-way junction as an architectural scaffold to pre-organize helical domains for recruiting the translation initiation machinery.

Secondary structure of a conserved domain in an intron of influenza A M1 mRNA

Influenza A virus utilizes RNA throughout infection. Little is known, however, about the roles of RNA structures. A previous bioinformatics survey predicted multiple regions of influenza A virus that are likely to generate evolutionarily conserved and stable RNA structures. One predicted conserved structure is in the pre-mRNA coding for essential proteins, M1 and M2. This structure starts 79 nucleotides downstream of the M2 mRNA 5' splice site. Here, a combination of biochemical structural mapping, mutagenesis, and NMR confirms the predicted three-way multibranch structure of this RNA. Imino proton NMR spectra reveal no change in secondary structure when 80 mM KCl is supplemented with 4 mM MgCl2. Optical melting curves in 1 M NaCl and in 100 mM KCl with 10 mM MgCl2 are very similar, with melting temperatures ∼14 °C higher than that for 100 mM KCl alone. These results provide a firm basis for designing experiments and potential therapeutics to test for function in cell culture.

Ten lessons with Carl Woese about RNA and comparative analysis

A few years before I started my graduate studies, Carl Woese was establishing a collaboration with his friend, colleague, and my PhD advisor, Harry Noller. Carl was introducing comparative methods to Harry's lab to determine the secondary structure for the 16S and 23S rRNAs. In addition to an experimental project that had minimal to no success, I was attempting to predict an RNA secondary structure from a single sequence. I determined after a few months that the complexity of RNA folding was much greater than ever anticipated. Ten lessons were learned about the dynamics of RNA folding, the comparative methods used to accurately predict the RNAs secondary structure and the beginnings of its tertiary structure, the use of comparative methods to reveal much more than ever anticipated about RNA structure, other applications beyond RNA structure, and the lessons about the process of scientific discovery.

Helix capping in RNA structure

Helices are an essential element in defining the three-dimensional architecture of structured RNAs. While internal basepairs in a canonical helix stack on both sides, the ends of the helix stack on only one side and are exposed to the loop side, thus susceptible to fraying unless they are protected. While coaxial stacking has long been known to stabilize helix ends by directly stacking two canonical helices coaxially, based on analysis of helix-loop junctions in RNA crystal structures, herein we describe helix capping, topological stacking of a helix end with a basepair or an unpaired nucleotide from the loop side, which in turn protects helix ends. Beyond the topological protection of helix ends against fraying, helix capping should confer greater stability onto the resulting composite helices. Our analysis also reveals that this general motif is associated with the formation of tertiary structure interactions. Greater knowledge about the dynamics at the helix-junctions in the secondary structure should enhance the prediction of RNA secondary structure with a richer set of energetic rules and help better understand the folding of a secondary structure into its three-dimensional structure. These together suggest that helix capping likely play a fundamental role in driving RNA folding.

Structural and thermodynamic signatures that define pseudotriloop RNA hairpins

Pseudotriloop (PTL) structures in RNAs have been recognized as essential elements in RNA folding and recognition of proteins. PTL structures are derived from hexaloops by formation of a cross-loop base pair leaving a triloop and 3' bulged out residue. Despite their common presence and functional importance, insufficient structural and thermodynamic data are available that can be used to predict formation of PTLs from sequence alone. Using NMR spectroscopy and UV-melting data we established factors that contribute to the formation and stability of PTL structures derived from hepatitis B virus and human foamy virus. The NMR data show that, besides the cross-loop base pair, also a 3' pyrimidine bulge and a G-C loop-closing base pair are primary determinants of PTL formation. By changing the G-C closing base pair into C-G, the PTL switches into a hexaloop. Comparison of these rules with regular triloop hairpins and PTLs from other sources is discussed as well as the conservation of a PTL in human foamy virus and other spumaretroviruses.

RNA modularity for synthetic biology

RNA molecules are highly modular components that can be used in a variety of contexts for building new metabolic, regulatory and genetic circuits in cells. The majority of synthetic RNA systems to date predominately rely on two-dimensional modularity. However, a better understanding and integration of three-dimensional RNA modularity at structural and functional levels is critical to the development of more complex, functional bio-systems and molecular machines for synthetic biology applications.

Secondary structure of a conserved domain in the intron of influenza A NS1 mRNA

Influenza A virus is a segmented single-stranded (−)RNA virus that causes substantial annual morbidity and mortality. The transcriptome of influenza A is predicted to have extensive RNA secondary structure. The smallest genome segment, segment 8, encodes two proteins, NS1 and NEP, via alternative splicing. A conserved RNA domain in the intron of segment 8 may be important for regulating production of NS1. Two different multi-branch loop structures have been proposed for this region. A combination of in vitro chemical mapping and isoenergetic microarray techniques demonstrate that the consensus sequence for this region folds into a hairpin conformation. These results provide an alternative folding for this region and a foundation for designing experiments to probe its functional role in the influenza life cycle.

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.

Twist-joints and double twist-joints in RNA structure

Analysis of available RNA crystal structures has allowed us to identify a new family of RNA arrangements that we call double twist-joints, or DTJs. Each DTJ is composed of a double helix that contains two bulges incorporated into different strands and separated from each other by 2 or 3 bp. At each bulge, the double helix is over-twisted, while the unpaired nucleotides of both bulges form a complex network of stacking and hydrogen-bonding with nucleotides of helical regions. In total, we identified 14 DTJ cases, which can be combined in three groups based on common structural characteristics. One DTJ is found in a functional center of the ribosome, another DTJ mediates binding of the pre-tRNA to the RNase P, and two more DTJs form the sensing domains in the glycine riboswitch.

Determining RNA three-dimensional structures using low-resolution data

Knowing the 3-D structure of an RNA is fundamental to understand its biological function. Nowadays X-ray crystallography and NMR spectroscopy are systematically applied to newly discovered RNAs. However, the application of these high-resolution techniques is not always possible, and thus scientists must turn to lower resolution alternatives. Here, we introduce a pipeline to systematically generate atomic resolution 3-D structures that are consistent with low-resolution data sets. We compare and evaluate the discriminative power of a number of low-resolution experimental techniques to reproduce the structure of the Escherichia coli tRNA(VAL) and P4-P6 domain of the Tetrahymena thermophila group I intron. We test single and combinations of the most accessible low-resolution techniques, i.e. hydroxyl radical footprinting (OH), methidiumpropyl-EDTA (MPE), multiplexed hydroxyl radical cleavage (MOHCA), and small-angle X-ray scattering (SAXS). We show that OH-derived constraints are accurate to discriminate structures at the atomic level, whereas EDTA-based constraints apply to global shape determination. We provide a guide for choosing which experimental techniques or combination of thereof is best in which context. The pipeline represents an important step towards high-throughput low-resolution RNA structure determination.

RNA2DMap: A Visual Exploration Tool of the Information in RNA's Higher-Order Structure

A new and emerging paradigm in molecular biology is revealing that RNA is implicated in nearly every aspect of the metabolism in the cell. To enhance our understanding of the function of these RNA molecules in the cell, it is essential that we have a complete understanding of their higher-order structures. While many computational tools have been developed to predict and analyse these higher-order RNA structures, few are able to visualize them for analytical purposes. In this paper, we present an interactive visualization tool of the secondary structure of RNA, named RNA2DMap. This program enables multiple-dimensions of information about RNA structure to be selected, customized and displayed to visually identify patterns and relationships. RNA2DMap facilitates the comparative analysis and understanding of RNAs that cannot be readily obtained with other graphical or text output from computer programs. Three use cases are presented to illustrate how RNA2DMap aids structural analysis.

RNA structure determination using SAXS data

Exploiting the experimental information from small-angle X-ray solution scattering (SAXS) in conjunction with structure prediction algorithms can be advantageous in the case of ribonucleic acids (RNA), where global restraints on the 3D fold are often lacking. Traditional usage of SAXS data often starts by attempting to reconstruct the molecular shape ab initio, which is subsequently used to assess the quality of a model. Here, an alternative strategy is explored whereby the models from a very large decoy set are directly sorted according to their fit to the SAXS data. For rapid computation of SAXS patterns, the method developed here makes use of a coarse-grained representation of RNA. It also accounts for the explicit treatment of the contribution to the scattering of water molecules and ions surrounding the RNA. The method, called Fast-SAXS-RNA, is first calibrated using a tRNA (tRNA-val) and then tested on the P4-P6 fragment of group I intron (P4-P6). Fast-SAXS-RNA is then used as a filter for decoy models generated by the MC-Fold and MC-Sym pipeline, a suite of RNA 3D all-atom structure algorithms that encode and exploit RNA 3D architectural principles. The ability of Fast-SAXS-RNA to discriminate native folds is tested against three widely used RNA molecules in molecular modeling benchmarks: the tRNA, the P4-P6, and a synthetic hairpin suspected to assemble into a homodimer. For each molecule, a large pool of decoys are generated, scored, and ranked using Fast-SAXS-RNA. The method is able to identify low-rmsd models among top ranking structures, for both tRNA and P4-P6. For the hairpin, the approach correctly identifies the dimeric state as the solution structure over the monomeric state and alternative secondary structures. The method offers a powerful strategy for recognizing native RNA conformations as well as multimeric assemblies and alternative secondary structures, thus enabling high-throughput RNA structure determination using SAXS data.

Promoting RNA helical stacking via A-minor junctions

RNA molecules take advantage of prevalent structural motifs to fold and assemble into well-defined 3D architectures. The A-minor junction is a class of RNA motifs that specifically controls coaxial stacking of helices in natural RNAs. A sensitive self-assembling supra-molecular system was used as an assay to compare several natural and previously unidentified A-minor junctions by native polyacrylamide gel electrophoresis and atomic force microscopy. This class of modular motifs follows a topological rule that can accommodate a variety of interchangeable A-minor interactions with distinct local structural motifs. Overall, two different types of A-minor junctions can be distinguished based on their functional self-assembling behavior: one group makes use of triloops or GNRA and GNRA-like loops assembling with helices, while the other takes advantage of more complex tertiary receptors specific for the loop to gain higher stability. This study demonstrates how different structural motifs of RNA can contribute to the formation of topologically equivalent helical stacks. It also exemplifies the need of classifying RNA motifs based on their tertiary structural features rather than secondary structural features. The A-minor junction rule can be used to facilitate tertiary structure prediction of RNAs and rational design of RNA parts for nanobiotechnology and synthetic biology.

Arrangement of 3D structural motifs in ribosomal RNA

Structural 3D motifs in RNA play an important role in the RNA stability and function. Previous studies have focused on the characterization and discovery of 3D motifs in RNA secondary and tertiary structures. However, statistical analyses of the distribution of 3D motifs along the RNA appear to be lacking. Herein, we present a novel strategy for evaluating the distribution of 3D motifs along the RNA chain and those motifs whose distributions are significantly non-random are identified. By applying it to the X-ray structure of the large ribosomal subunit from Haloarcula marismortui, helical motifs were found to cluster together along the chain and in the 3D structure, whereas the known tetraloops tend to be sequentially and spatially dispersed. That the distribution of key structural motifs such as tetraloops differ significantly from a random one suggests that our method could also be used to detect novel 3D motifs of any size in sufficiently long/large RNA structures. The motif distribution type can help in the prediction and design of 3D structures of large RNA molecules.

Tertiary architecture of the Oceanobacillus iheyensis group II intron

Group II introns are large ribozymes that act as self-splicing and retrotransposable RNA molecules. They are of great interest because of their potential evolutionary relationship to the eukaryotic spliceosome, their continued influence on the organization of many genomes in bacteria and eukaryotes, and their potential utility as tools for gene therapy and biotechnology. One of the most interesting features of group II introns is their relative lack of nucleobase conservation and covariation, which has long suggested that group II intron structures are stabilized by numerous unusual tertiary interactions and backbone-mediated contacts. Here, we provide a detailed description of the tertiary interaction networks within the Oceanobacillus iheyensis group IIC intron, for which a crystal structure was recently solved to 3.1 A resolution. The structure can be described as a set of several intricately constructed tertiary interaction nodes, each of which contains a core of extended stacking networks and elaborate motifs. Many of these nodes are surrounded by a web of ribose zippers, which appear to further stabilize local structure. As predicted from biochemical and genetic studies, the group II intron provides a wealth of new information on strategies for RNA folding and tertiary structural organization.

Three-way RNA junctions with remote tertiary contacts: a recurrent and highly versatile fold

Three-way junction RNAs adopt a recurrent Y shape when two of the helices form a coaxial stack and the third helix establishes one or more tertiary contacts several base pairs away from the junction. In this review, the structure, distribution, and functional relevance of these motifs are examined. Structurally, the folds exhibit conserved junction topologies, and the distal tertiary interactions play a crucial role in determining the final shape of the structures. The junctions and remote tertiary contacts behave as flexible hinge motifs that respond to changes in the other region, providing these folds with switching mechanisms that have been shown to be functionally useful in a variety of contexts. In addition, the juxtaposition of RNA domains at the junction and at the distal tertiary complexes enables the RNA helices to adopt unusual conformations that are frequently used by proteins, RNA molecules, and antibiotics as platforms for specific binding. As a consequence of these properties, Y-shaped junctions are widely distributed in all kingdoms of life, having been observed in small naked RNAs such as riboswitches and ribozymes or embedded in complex ribonucleoprotein systems like ribosomal RNAs, RNase P, or the signal recognition particle. In all cases, the folds were found to play an essential role for the functioning or assembly of the RNA or ribonucleoprotein systems that contain them.

Correlation of RNA secondary structure statistics with thermodynamic stability and applications to folding

The accurate prediction of the secondary and tertiary structure of an RNA with different folding algorithms is dependent on several factors, including the energy functions. However, an RNA higher-order structure cannot be predicted accurately from its sequence based on a limited set of energy parameters. The inter- and intramolecular forces between this RNA and other small molecules and macromolecules, in addition to other factors in the cell such as pH, ionic strength, and temperature, influence the complex dynamics associated with transition of a single stranded RNA to its secondary and tertiary structure. Since all of the factors that affect the formation of an RNAs 3D structure cannot be determined experimentally, statistically derived potential energy has been used in the prediction of protein structure. In the current work, we evaluate the statistical free energy of various secondary structure motifs, including base-pair stacks, hairpin loops, and internal loops, using their statistical frequency obtained from the comparative analysis of more than 50,000 RNA sequences stored in the RNA Comparative Analysis Database (rCAD) at the Comparative RNA Web (CRW) Site. Statistical energy was computed from the structural statistics for several datasets. While the statistical energy for a base-pair stack correlates with experimentally derived free energy values, suggesting a Boltzmann-like distribution, variation is observed between different molecules and their location on the phylogenetic tree of life. Our statistical energy values calculated for several structural elements were utilized in the Mfold RNA-folding algorithm. The combined statistical energy values for base-pair stacks, hairpins and internal loop flanks result in a significant improvement in the accuracy of secondary structure prediction; the hairpin flanks contribute the most.

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.

The distributions, mechanisms, and structures of metabolite-binding riboswitches

BACKGROUND

Riboswitches are noncoding RNA structures that appropriately regulate genes in response to changing cellular conditions. The expression of many proteins involved in fundamental metabolic processes is controlled by riboswitches that sense relevant small molecule ligands. Metabolite-binding riboswitches that recognize adenosylcobalamin (AdoCbl), thiamin pyrophosphate (TPP), lysine, glycine, flavin mononucleotide (FMN), guanine, adenine, glucosamine-6-phosphate (GlcN6P), 7-aminoethyl 7-deazaguanine (preQ1), and S-adenosylmethionine (SAM) have been reported.

RESULTS

We have used covariance model searches to identify examples of ten widespread riboswitch classes in the genomes of organisms from all three domains of life. This data set rigorously defines the phylogenetic distributions of these riboswitch classes and reveals how their gene control mechanisms vary across different microbial groups. By examining the expanded aptamer sequence alignments resulting from these searches, we have also re-evaluated and refined their consensus secondary structures. Updated riboswitch structure models highlight additional RNA structure motifs, including an unusual double T-loop arrangement common to AdoCbl and FMN riboswitch aptamers, and incorporate new, sometimes noncanonical, base-base interactions predicted by a mutual information analysis.

CONCLUSION

Riboswitches are vital components of many genomes. The additional riboswitch variants and updated aptamer structure models reported here will improve future efforts to annotate these widespread regulatory RNAs in genomic sequences and inform ongoing structural biology efforts. There remain significant questions about what physiological and evolutionary forces influence the distributions and mechanisms of riboswitches and about what forms of regulation substitute for riboswitches that appear to be missing in certain lineages.

Structural principles from large RNAs

Since the year 2000 a number of large RNA three-dimensional structures have been determined by X-ray crystallography. Structures composed of more than 100 nucleotide residues include the signal recognition particle RNA, group I intron, the GlmS ribozyme, RNAseP RNA, and ribosomal RNAs from Haloarcula morismortui, Escherichia coli, Thermus thermophilus, and Deinococcus radiodurans. These large RNAs are constructed from the same secondary and tertiary structural motifs identified in smaller RNAs but appear to have a larger organizational architecture. They are dominated by long continuous interhelical base stacking, tend to segregate into domains, and are planar in overall shape as opposed to their globular protein counterparts. These findings have consequences in RNA folding, intermolecular interaction, and packing, in addition to studies of design and engineering and structure prediction.

Structures, kinetics, thermodynamics, and biological functions of RNA hairpins

Most RNA comprises one strand and therefore can fold back on itself to form complex structures. At the heart of these structures is the hairpin, which is composed of a stem having Watson-Crick base pairing and a loop wherein the backbone changes directionality. First, we review the structure of hairpins including diversity in the stem, loop, and closing base pair. The function of RNA hairpins in biology is discussed next, including roles for isolated hairpins, as well as hairpins in the context of complex tertiary structures. We describe the kinetics and thermodynamics of hairpin folding including models for hairpin folding, folding transition states, and the cooperativity of folding. Lastly, we discuss some ways in which hairpins can influence the folding and function of tertiary structures, both directly and indirectly. RNA hairpins provide a simple means of controlling gene expression that can be understood in the language of physical chemistry.

Cartilage-hair hypoplasia-associated mutations in the RNase MRP P3 domain affect RNA folding and ribonucleoprotein assembly

Cartilage-hair hypoplasia (CHH) is caused by mutations in the gene encoding the RNA component of RNase MRP. Currently it is unknown how these mutations affect the function of this endoribonuclease. In this study we investigated the effect of mutations in the P3 domain on protein binding and RNA folding. Our data demonstrate that a number of P3 nucleotide substitutions reduced the efficiency of its interaction with Rpp25 and Rpp20, two protein subunits binding as a heterodimer to this domain. The CHH-associated 40G>A substitution, as well as the replacement of residue 47, almost completely abrogated Rpp25 and Rpp20 binding in different assays. Also other CHH-associated P3 mutations reduced the efficiency by which the RNase MRP RNA is bound by Rpp25-Rpp20. These data demonstrate that the most important residues for binding of the Rpp25-Rpp20 dimer reside in the apical stem-loop of the P3 domain. Structural analyses by NMR not only showed that this loop may adopt a pseudo-triloop structure, but also demonstrated that the 40G>A substitution alters the folding of this part of the P3 domain. Our data are the first to provide insight into the molecular mechanism by which CHH-associated mutations affect the function of RNase MRP.

A comparative analysis of the triloops in all high-resolution RNA structures reveals sequence structure relationships

Despite an increasing number of experimentally determined RNA structures, the gap between the number of structures and that of RNA families is still growing. To overcome this limitation, efficient and reliable RNA modeling methodologies must be developed. In order to reach this goal, here, we show how triloop sequence-structure relationships have been inferred through a systematic analysis of all triloops found in available high-resolution structures. The structural annotation of all triloops allowed us to define discrete states of the triloop's conformational space, and therefore an explicit sequence-to-structure relation. The sequence-structure relationships inferred from this explicit relation are presented in a convenient modeling table that provides a limited set of possible three-dimensional structures given any triloop sequence. The table is indexed by the two nucleotides that form the triloop's flanking base pair, since they are shown to provide the most information about the triloop three-dimensional structures. We also report the observations in the X-ray crystallographic structures of important conformational variations, which we believe might be the result of RNA dynamic.

Active suppressor tRNAs with a double helix between the D- and T-loops

One of the most conserved elements of the tRNA structure is the reverse-Hoogsteen base-pair T54--A58 in the T-loop, which plays a major role in the maintenance of the standard L-shape conformation. Here, we present the results of in vivo selection of 51 active suppressor tRNA clones, none of which contains base-pair T54--A58. In 49 clones, we found two regions in the D and T-loops that are complementary to each other. This finding suggests the existence of an inter-loop double helix consisting of three base-pairs, which could have the same role as base-pair T54--A58 in the fixation of the juxtaposition of the two helical domains within the L-shape. From this point of view, the appearance of the inter-loop double helix represents a compensatory effect for the absence of base-pair T54--A58. The results shed new light on the role of different elements of the tRNA structure in the formation of the standard L-shape conformation and on the possibility of synonymous replacements of one arrangement by another in functional RNA molecules.

Thermodynamics and NMR studies on Duck, Heron and Human HBV encapsidation signals

Hepatitis B virus (HBV) replication is initiated by binding of its reverse transcriptase (P) to the apical stem-loop (AL) and primer loop (PL) of epsilon, a highly conserved RNA element at the 5'-end of the RNA pregenome. Mutation studies on duck/heron and human in vitro systems have shown similarities but also differences between their P-epsilon interaction. Here, NMR and UV thermodynamic data on AL (and PL) from these three species are presented. The stabilities of the duck and heron ALs were found to be similar, and much lower than that of human. NMR data show that this low stability stems from an 11-nt internal bulge destabilizing the stem of heron AL. In duck, although structured at low temperature, this region also forms a weak point as its imino resonances broaden to disappearance between 30 and 35 degrees C well below the overall AL melting temperature. Surprisingly, the duck- and heron ALs were both found to be capped by a stable well-structured UGUU tetraloop. All avian ALs are expected to adhere to this because of their conserved sequence. Duck PL is stable and structured and, in view of sequence similarities, the same is expected for heron - and human PL.

The interaction networks of structured RNAs

All pairwise interactions occurring between bases which could be detected in three-dimensional structures of crystallized RNA molecules are annotated on new planar diagrams. The diagrams attempt to map the underlying complex networks of base–base interactions and, especially, they aim at conveying key relationships between helical domains: co-axial stacking, bending and all Watson–Crick as well as non-Watson–Crick base pairs. Although such wiring diagrams cannot replace full stereographic images for correct spatial understanding and representation, they reveal structural similarities as well as the conserved patterns and distances between motifs which are present within the interaction networks of folded RNAs of similar or unrelated functions. Finally, the diagrams could help devising methods for meaningfully transforming RNA structures into graphs amenable to network analysis.

Characteristics of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) rRNA genes of Apis mellifera (Insecta: Hymenoptera): structure, organization, and retrotransposable elements

As an accompanying manuscript to the release of the honey bee genome, we report the entire sequence of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) ribosomal RNA (rRNA)-encoding gene sequences (rDNA) and related internally and externally transcribed spacer regions of Apis mellifera (Insecta: Hymenoptera: Apocrita). Additionally, we predict secondary structures for the mature rRNA molecules based on comparative sequence analyses with other arthropod taxa and reference to recently published crystal structures of the ribosome. In general, the structures of honey bee rRNAs are in agreement with previously predicted rRNA models from other arthropods in core regions of the rRNA, with little additional expansion in non-conserved regions. Our multiple sequence alignments are made available on several public databases and provide a preliminary establishment of a global structural model of all rRNAs from the insects. Additionally, we provide conserved stretches of sequences flanking the rDNA cistrons that comprise the externally transcribed spacer regions (ETS) and part of the intergenic spacer region (IGS), including several repetitive motifs. Finally, we report the occurrence of retrotransposition in the nuclear large subunit rDNA, as R2 elements are present in the usual insertion points found in other arthropods. Interestingly, functional R1 elements usually present in the genomes of insects were not detected in the honey bee rRNA genes. The reverse transcriptase products of the R2 elements are deduced from their putative open reading frames and structurally aligned with those from another hymenopteran insect, the jewel wasp Nasonia (Pteromalidae). Stretches of conserved amino acids shared between Apis and Nasonia are illustrated and serve as potential sites for primer design, as target amplicons within these R2 elements may serve as novel phylogenetic markers for Hymenoptera. Given the impending completion of the sequencing of the Nasonia genome, we expect our report eventually to shed light on the evolution of the hymenopteran genome within higher insects, particularly regarding the relative maintenance of conserved rDNA genes, related variable spacer regions and retrotransposable elements.

Identification of conserved regions in the plastid genome: implications for DNA barcoding and biological function

All oligonucleotides of the sugarcane chloroplast genome that are conserved in one or more of 36 other completed plastid genomes have been identified by computer-assisted sequence comparison. These regions are of interest because they (i) are indicative of strong selection pressures to maintain specific nucleotide sequences that may yield insights into plastid biology and (ii) can be used as priming sites for amplifying intervening sequences that represent potential DNA barcodes for species identification. The majority of conserved sites are located in the inverted repeat (IR) region, but several sites in the single copy region (predominantly in tRNA and psa/psb genes) are conserved among chloroplasts of all higher plants examined here. Of particular interest are protein coding regions that have been conserved at the nucleotide level, as these may be involved in transcript regulation. This analysis also provides the basis for rational design of a DNA barcode for plastids, and several potential barcode regions have been identified. In particular, two oligonucleotides of length 33 and 25, and separated by approximately 362 nucleotides, are found in all cyanobacteria, red, brown and green algae, as well as diatoms, euglenids, apicomplexans and land plants that have been examined to date. Their widespread occurrence makes the intervening sequence a universal marker for all photosynthetic lineages. Analysis of 160 GenBank accessions illustrates that this region discriminates many algae at the species level, but lacks sufficient variation among the more recently diverged land plants to serve as a single DNA barcode for this taxon. However, this marker should be particularly useful for the DNA barcoding of algal lineages and lichens, as well as for environmental sampling. More rapidly evolving regions of the plastid genome also identified here serve as a starting point to design and test barcodes for more narrowly defined lineages, including the more recently diverged angiosperms.

The UAA/GAN internal loop motif: a new RNA structural element that forms a cross-strand AAA stack and long-range tertiary interactions

Analysis of aligned RNA sequences and high-resolution crystal structures has revealed a new RNA structural element, termed the UAA/GAN motif. Found in internal loops of the 23 S rRNA, as well as in RNase P RNA and group I and II introns, this six-nucleotide motif adopts a distinctive local structure that includes two base-pairs with non-canonical conformations and three conserved adenine bases, which form a cross-strand AAA stack in the minor groove. Most importantly, the motif invariably forms long-range tertiary contacts, as the AAA stack typically forms A-minor interactions and the flipped-out N nucleotide forms additional contacts that are specific to the structural context of each loop. The widespread presence of this motif and its propensity to form long-range contacts suggest that it plays a critical role in defining the architectures of structured RNAs.

Single nucleotide RNA choreography

New structural analysis methods, and a tree formalism re-define and expand the RNA motif concept, unifying what previously appeared to be disparate groups of structures. We find RNA tetraloops at high frequencies, in new contexts, with unexpected lengths, and in novel topologies. The results, with broad implications for RNA structure in general, show that even at this most elementary level of organization, RNA tolerates astounding variation in conformation, length, sequence and context. However the variation is not random; it is well-described by four distinct modes, which are 3-2 switches (backbone topology variations), insertions, deletions and strand clips.

A structural model for the large subunit of the mammalian mitochondrial ribosome

Protein translation is essential for all forms of life and is conducted by a macromolecular complex, the ribosome. Evolutionary changes in protein and RNA sequences can affect the 3D organization of structural features in ribosomes in different species. The most dramatic changes occur in animal mitochondria, whose genomes have been reduced and altered significantly. The RNA component of the mitochondrial ribosome (mitoribosome) is reduced in size, with a compensatory increase in protein content. Until recently, it was unclear how these changes affect the 3D structure of the mitoribosome. Here, we present a structural model of the large subunit of the mammalian mitoribosome developed by combining molecular modeling techniques with cryo-electron microscopic data at 12.1A resolution. The model contains 93% of the mitochondrial rRNA sequence and 16 mitochondrial ribosomal proteins in the large subunit of the mitoribosome. Despite the smaller mitochondrial rRNA, the spatial positions of RNA domains known to be involved directly in protein synthesis are essentially the same as in bacterial and archaeal ribosomes. However, the dramatic reduction in rRNA content necessitates evolution of unique structural features to maintain connectivity between RNA domains. The smaller rRNA sequence also limits the likelihood of tRNA binding at the E-site of the mitoribosome, and correlates with the reduced size of D-loops and T-loops in some animal mitochondrial tRNAs, suggesting co-evolution of mitochondrial rRNA and tRNA structures.

Accurate energies of hydrogen bonded nucleic acid base pairs and triplets in tRNA tertiary interactions

Tertiary interactions are crucial in maintaining the tRNA structure and functionality. We used a combined sequence analysis and quantum mechanics approach to calculate accurate energies of the most frequent tRNA tertiary base pairing interactions. Our analysis indicates that six out of the nine classical tertiary interactions are held in place mainly by H-bonds between the bases. In the remaining three cases other effects have to be considered. Tertiary base pairing interaction energies range from -8 to -38 kcal/mol in yeast tRNA(Phe) and are estimated to contribute roughly 25% of the overall tRNA base pairing interaction energy. Six analyzed posttranslational chemical modifications were shown to have minor effect on the geometry of the tertiary interactions. Modifications that introduce a positive charge strongly stabilize the corresponding tertiary interactions. Non-additive effects contribute to the stability of base triplets.

Topology of three-way junctions in folded RNAs

The three-way junctions contained in X-ray structures of folded RNAs have been compiled and analyzed. Three-way junctions with two helices approximately coaxially stacked can be divided into three main families depending on the relative lengths of the segments linking the three Watson-Crick helices. Each family has topological characteristics with some conservation in the non-Watson-Crick pairs within the linking segments as well as in the types of contacts between the segments and the helices. The most populated family presents tertiary interactions between two helices as well as extensive shallow/minor groove contacts between a linking segment and the third helix. On the basis of the lengths of the linking segments, some guidelines could be deduced for choosing a topology for a three-way junction on the basis of a secondary structure. Examples and prediction bas'ed on those rules are discussed.

Assessment of adenyl residue reactivity within model nucleic acids by surface enhanced Raman spectroscopy

We rank the reactivity of the adenyl residues (A) of model DNA and RNA molecules with electropositive subnano size [Ag]n+ sites as a function of nucleic acid primary sequences and secondary structures and in the presence of biological amounts of Cl- and Na+ or Mg2+ ions. In these conditions A is markedly more reactive than any other nucleic acid bases. A reactivity is higher in ribo (r) than in deoxyribo (d) species [pA>pdA and (pA)n>>(pdA)n]. Base pairing decreases A reactivity in corresponding duplexes but much less in r than in d. In linear single and paired dCAG or dGAC loci, base stacking inhibits A reactivity even if A is bulged or mispaired (A.A). dA tracts are highly reactive only when dilution prevents self-association and duplex structures. In d hairpins the solvent-exposed A residues are reactive in CAG and GAC triloops and even more in ATC loops. Among the eight rG1N2R3A4 loops, those bearing a single A (A4) are the least reactive. The solvent-exposed A2 is reactive, but synergistic structural transitions make the initially stacked A residues of any rGNAA loop much more reactive. Mg2+ cross-bridging single strands via phosphates may screen A reactivity. In contrast d duplexes cross-bridging enables "A flipping" much more in rA.U pairs than in dA.T. Mg2+ promotes A reactivity in unpaired strands. For hairpins Mg2+ binding stabilizes the stems, but according to A position in the loops, A reactivity may be abolished, reduced, or enhanced. It is emphasized that not only accessibility but also local flexibility, concerted docking, and cation and anion binding control A reactivity.

On the structural features of hairpin triloops in rRNA: from nucleotide to global conformational change upon ligand binding

RNA structure can be viewed as both a construct composed of various structural motifs and a flexible polymer that is substantially influenced by its environment. In this light, the present paper represents an attempt to reconcile the two standpoints. By using the 3D structures both of four (16S and 23S) portions of unbound 50S, H50S, and T30S ribosomal subunits and of 38 large ribonucleoligand complexes as the starting point, the behavior, which is induced by ligand binding, of 73 hairpin triloops with closing g-c and c-g base pairs was investigated using root-mean-square deviation (RMSD) approach and pseudotorsional (eta,theta) convention at the nucleotide-by-nucleotide level. Triloops were annotated in accordance with a recent proposal of geometric nomenclature. A simple measure for the determination of the strain of a triloop is introduced. It is believed that a possible classification of the interior triloops, based on the 2D eta-theta unique path, will aid to conceive their local behavior upon ligand binding. All rRNA residues in contact with ligands as well as regions of considerable conformational changes upon complex formation were identified. The analysis offers the answer to: how proximal to and how far from the actual ligand-binding sites the structural changes occur?

Functional hammerhead ribozymes naturally encoded in the genome of Arabidopsis thaliana

The hammerhead ribozyme (HHRz) is an autocatalytic RNA motif found in subviral plant pathogens and transcripts of repetitive DNA sequences in animals. Here, we report the discovery and characterization of unique HHRzs encoded in a plant genome. Two novel sequences were identified on chromosome IV of Arabidopsis thaliana in a database search, which took into account recently defined structural requirements. The HHRzs are expressed in several tissues and coexist in vivo as both cleaved and noncleaved species. In vitro, both sequences cleave efficiently at physiological Mg(2+) concentrations, indicative of functional loop-loop interactions. Kinetic analysis of loop nucleotide variants was used to determine a three-dimensional model of these tertiary interactions. Based on these results, on the lack of infectivity of hammerhead-carrying viroids in Arabidopsis, and on extensive sequence comparisons, we propose that the ribozyme sequences did not invade this plant by horizontal transfer but have evolved independently to perform a specific, yet unidentified, biological function.

Nucleic acid crystallography: current progress

Fifty years after the publication of the DNA double helix model by Watson and Crick, new nucleic acid structures keep emerging at an ever-increasing rate. The past three years have brought a flurry of new oligonucleotide structures, including those of a Hoogsteen-paired DNA duplex, Holliday junctions, DNA-drug complexes, quadruplexes, a host of RNA motifs and various nucleic acid analogues. Major advances were also made in terms of the structure and function of catalytic RNAs. These range from improved models of the phosphodiester cleavage reactions catalyzed by the hairpin and hepatitis delta virus ribozymes to the visualization of a complete active site of a group I self-splicing intron with bound 5'- and 3'-exons. These triumphs are complemented by a refined understanding of cation-nucleic-acid interactions and new routes to the generation of derivatives for phasing of DNA and RNA structures.

Structural transitions induced by the interaction between tRNA(Gly) and the Bacillus subtilis glyQS T box leader RNA

The T box system regulates expression of amino acid-related genes in Gram-positive bacteria through premature termination of transcription. Synthesis of the full-length mRNA requires stabilization of an antiterminator element in the 5' untranslated leader RNA by the cognate uncharged tRNA. tRNA(Gly)-dependent antitermination of the Bacillus subtilis glyQS gene (encoding glycyl-tRNA synthetase) can be reproduced in a purified in vitro transcription system, indicating that the nascent transcript is sufficient for interaction with the tRNA. Genetic analyses previously demonstrated base pairing of a single codon in the leader RNA with the tRNA anticodon, and between the antiterminator and the tRNA acceptor end. In this study, we established conditions for specific binding of tRNA(Gly) to glyQS leader RNA generated by phage T7 RNA polymerase. Structural mapping studies revealed tRNA(Gly)-induced protection in the glyQS leader RNA at the two known sites of interaction with the tRNA, as well as at other regions between these sites. The proposed tRNA-dependent structural switch between the competing terminator and antiterminator forms of the leader RNA was demonstrated directly. Changes in tRNA(Gly) upon binding to glyQS leader RNA were detected in the anticodon loop, consistent with pairing with the specifier sequence, and in the highly conserved G19 in the D-loop, similar to effects induced by codon-anticodon interaction in the ribosome. This study provides biochemical evidence for direct interaction of tRNA(Gly) with full-length in vitro transcribed glyQS leader RNA, and an initial view of structural modulations of both RNA partners within the complex.

The application of cluster analysis in the intercomparison of loop structures in RNA

We have developed a computational approach for the comparison and classification of RNA loop structures. Hairpin or interior loops identified in atomic resolution RNA structures were intercompared by conformational matching. The root-mean-square deviation (RMSD) values between all pairs of RNA fragments of interest, even if from different molecules, are calculated. Subsequently, cluster analysis is performed on the resulting matrix of RMSD distances using the unweighted pair group method with arithmetic mean (UPGMA). The cluster analysis objectively reveals groups of folds that resemble one another. To demonstrate the utility of the approach, a comprehensive analysis of all the terminal hairpin tetraloops that have been observed in 15 RNA structures that have been determined by X-ray crystallography was undertaken. The method found major clusters corresponding to the well-known GNRA and UNCG types. In addition, two tetraloops with the unusual primary sequence UMAC (M is A or C) were successfully assigned to the GNRA cluster. Larger loop structures were also examined and the clustering results confirmed the occurrence of variations of the GNRA and UNCG tetraloops in these loops and provided a systematic means for locating them. Nineteen examples of larger loops that closely resemble either the GNRA or UNCG tetraloop were found in the large ribosomal RNAs. When the clustering approach was extended to include all structures in the SCOR database, novel relationships were detected including one between the ANYA motif and a less common folding of the GAAA tetraloop sequence.

Selection on the structural stability of a ribosomal RNA expansion segment in Daphnia obtusa

The high rate of sequence divergence in nuclear ribosomal RNA (rRNA) expansion segments offers a unique opportunity to study the importance of natural selection in their evolution. To this end, we polymerase chain reaction amplified and cloned a 589-nt fragment of the 18S rRNA gene containing expansion segments 43/e1 and 43/e4 from six individual Daphnia obtusa from four populations. We screened 2,588 clones using single-stranded conformation polymorphism analysis and identified 103 unique haplotype sequences. We detected two pairs of indel sites in segment 43/e4 that complement each other when the secondary structure of the linear sequence is formed. Seven of the 12 observed combinations of length variants at these four sites (haplotypes) are shared between individuals from different populations, which may suggest that some of the length variation was present in their common ancestor. Haplotypes with uncompensated indels were only observed at low frequencies, while compensated indel haplotypes were found at a wide range of frequencies, supporting the hypothesis that the energetic stability of expansion segments is a trait under natural selection. In addition, there was strong linkage disequilibrium between the four complementary indel sites, particularly those that pair with one another in the secondary structure. Despite selection against unpaired bulges at these four indel sites, some nucleotides that form unpaired bulges are highly conserved in segment 43/e4, indicating that they are under a different selective constraint, possibly due to their role in higher level structural interactions.