Non-Watson-Crick basepairing and hydration in RNA motifs: molecular dynamics of 5S rRNA loop E

Explicit solvent and counterion molecular dynamics simulations have been carried out for a total of >80 ns on the bacterial and spinach chloroplast 5S rRNA Loop E motifs. The Loop E sequences form unique duplex architectures composed of seven consecutive non-Watson-Crick basepairs. The starting structure of spinach chloroplast Loop E was modeled using isostericity principles, and the simulations refined the geometries of the three non-Watson-Crick basepairs that differ from the consensus bacterial sequence. The deep groove of Loop E motifs provides unique sites for cation binding. Binding of Mg(2+) rigidifies Loop E and stabilizes its major groove at an intermediate width. In the absence of Mg(2+), the Loop E motifs show an unprecedented degree of inner-shell binding of monovalent cations that, in contrast to Mg(2+), penetrate into the most negative regions inside the deep groove. The spinach chloroplast Loop E shows a marked tendency to compress its deep groove compared with the bacterial consensus. Structures with a narrow deep groove essentially collapse around a string of Na(+) cations with long coordination times. The Loop E non-Watson-Crick basepairing is complemented by highly specific hydration sites ranging from water bridges to hydration pockets hosting 2 to 3 long-residing waters. The ordered hydration is intimately connected with RNA local conformational variations.     

G.U base pairing motifs in ribosomal RNA.

An increasing number of recognition mechanisms in RNA are found to involve G.U base pairs. In order to detect new functional sites of this type, we exhaustively analyzed the sequence alignments and secondary structures of eubacterial and chloroplast 16S and 23S rRNA, seeking positions with high levels of G.U pairs. Approximately 120 such sites were identified and classified according to their secondary structure and sequence environment. Overall biases in the distribution of G.U pairs are consistent with previously proposed structural rules: the side of the wobble pair that is subject to a loss of stacking is preferentially exposed to a secondary structure loop, where stacking is not as essential as in helical regions. However, multiple sites violate these rules and display highly conserved G.U pairs in orientations that could cause severe stacking problems. In addition, three motifs displaying a conserved G.U pair in a specific sequence/structure environment occur at an unusually high frequency. These motifs, of which two had not been reported before, involve sequences 5''UG3'' 3''GA5'' and 5''UG3'' 3''GU5'', as well as G.U pairs flanked by a bulge loop 3'' of U. The possible structures and functions of these recurrent motifs are discussed.     

The 5S rRNA loop E: chemical probing and phylogenetic data versus crystal structure.

A significant fraction of the bases in a folded, structured RNA molecule participate in noncanonical base pairing interactions, often in the context of internal loops or multi-helix junction loops. The appearance of each new high-resolution RNA structure provides welcome data to guide efforts to understand and predict RNA 3D structure, especially when the RNA in question is a functionally conserved molecule. The recent publication of the crystal structure of the "Loop E" region of bacterial 5S ribosomal RNA is such an event [Correll CC, Freeborn B, Moore PB, Steitz TA, 1997, Cell 91:705-712]. In addition to providing more examples of already established noncanonical base pairs, such as purine-purine sheared pairings, trans-Hoogsteen UA, and GU wobble pairs, the structure provides the first high-resolution views of two new purine-purine pairings and a new GU pairing. The goal of the present analysis is to expand the capabilities of both chemical probing and phylogenetic analysis to predict with greater accuracy the structures of RNA molecules. First, in light of existing chemical probing data, we investigate what lessons could be learned regarding the interpretation of this widely used method of RNA structure probing. Then we analyze the 3D structure with reference to molecular phylogeny data (assuming conservation of function) to discover what alternative base pairings are geometrically compatible with the structure. The comparisons between previous modeling efforts and crystal structures show that the intricate involvements of ions and water molecules in the maintenance of non-Watson-Crick pairs render the process of correctly identifying the interacting sites in such pairs treacherous, except in cases of trans-Hoogsteen A/U or sheared A/G pairs for the adenine N1 site. The phylogenetic analysis identifies A/A, A/C, A/U and C/A, C/C, and C/U pairings isosteric with sheared A/G, as well as A/A and A/C pairings isosteric with both G/U and G/G bifurcated pairings. Thus, each non-Watson-Crick pair could be characterized by a phylogenetic signature of variations between isosteric-like pairings. In addition to the conservative changes, which form a dictionary of pairings isosterically compatible with those observed in the crystal structure, concerted changes involving several base pairs also occur. The latter covariations may indicate transitions between related but distinctive motifs within the loop E of 5S ribosomal RNA.     

The lonepair triloop: a new motif in RNA structure

The lonepair triloop (LPTL) is an RNA structural motif that contains a single ("lone") base-pair capped by a hairpin loop containing three nucleotides. The two nucleotides immediately outside of this motif (5' and 3' to the lonepair) are not base-paired to one another, restricting the length of this helix to a single base-pair. Four examples of this motif, along with three tentative examples, were initially identified in the 16S and 23S rRNAs with covariation analysis. An evaluation of the recently determined crystal structures of the Thermus thermophilus 30S and Haloarcula marismortui 50S ribosomal subunits revealed the authenticity for all of these proposed interactions and identified 16 more LPTLs in the 5S, 16S and 23S rRNAs. This motif is found in the T loop in the tRNA crystal structures. The lonepairs are positioned, in nearly all examples, immediately 3' to a regular secondary structure helix and are stabilized by coaxial stacking onto this flanking helix. In all but two cases, the nucleotides in the triloop are involved in a tertiary interaction with another section of the rRNA, establishing an overall three-dimensional function for this motif. Of these 24 examples, 14 occur in multi-stem loops, seven in hairpin loops and three in internal loops. While the most common lonepair, U:A, occurs in ten of the 24 LPTLs, the remaining 14 LPTLs contain seven different base-pair types. Only a few of these lonepairs adopt the standard Watson-Crick base-pair conformations, while the majority of the base-pairs have non-standard conformations. While the general three-dimensional conformation is similar for all examples of this motif, characteristic differences lead to several subtypes present in different structural environments. At least one triloop nucleotide in 22 of the 24 LPTLs in the rRNAs and tRNAs forms a tertiary interaction with another part of the RNA. When a LPTL containing the GNR or UYR triloop sequence forms a tertiary interaction with the first (and second) triloop nucleotide, it recruits a fourth nucleotide to mediate stacking and mimic the tetraloop conformation. Approximately half of the LPTL motifs are in close association with proteins. The majority of these LPTLs are positioned at sites in rRNAs that are conserved in the three phylogenetic domains; a few of these occur in regions of the rRNA associated with ribosomal function, including the presumed site of peptidyl transferase activity in the 23S rRNA.     

FR3D: finding local and composite recurrent structural motifs in RNA 3D structures

New methods are described for finding recurrent three-dimensional (3D) motifs in RNA atomic-resolution structures. Recurrent RNA 3D motifs are sets of RNA nucleotides with similar spatial arrangements. They can be local or composite. Local motifs comprise nucleotides that occur in the same hairpin or internal loop. Composite motifs comprise nucleotides belonging to three or more different RNA strand segments or molecules. We use a base-centered approach to construct efficient, yet exhaustive search procedures using geometric, symbolic, or mixed representations of RNA structure that we implement in a suite of MATLAB programs, “Find RNA 3D” (FR3D). The first modules of FR3D preprocess structure files to classify base-pair and -stacking interactions. Each base is represented geometrically by the position of its glycosidic nitrogen in 3D space and by the rotation matrix that describes its orientation with respect to a common frame. Base-pairing and base-stacking interactions are calculated from the base geometries and are represented symbolically according to the Leontis/Westhof basepairing classification, extended to include base-stacking. These data are stored and used to organize motif searches. For geometric searches, the user supplies the 3D structure of a query motif which FR3D uses to find and score geometrically similar candidate motifs, without regard to the sequential position of their nucleotides in the RNA chain or the identity of their bases. To score and rank candidate motifs, FR3D calculates a geometric discrepancy by rigidly rotating candidates to align optimally with the query motif and then comparing the relative orientations of the corresponding bases in the query and candidate motifs. Given the growing size of the RNA structure database, it is impossible to explicitly compute the discrepancy for all conceivable candidate motifs, even for motifs with less than ten nucleotides. The screening algorithm that we describe finds all candidate motifs whose geometric discrepancy with respect to the query motif falls below a user-specified cutoff discrepancy. This technique can be applied to RMSD searches. Candidate motifs identified geometrically may be further screened symbolically to identify those that contain particular basepair types or base-stacking arrangements or that conform to sequence continuity or nucleotide identity constraints. Purely symbolic searches for motifs containing user-defined sequence, continuity and interaction constraints have also been implemented. We demonstrate that FR3D finds all occurrences, both local and composite and with nucleotide substitutions, of sarcin/ricin and kink-turn motifs in the 23S and 5S ribosomal RNA 3D structures of the H. marismortui 50S ribosomal subunit and assigns the lowest discrepancy scores to bona fide examples of these motifs. The search algorithms have been optimized for speed to allow users to search the non-redundant RNA 3D structure database on a personal computer in a matter of minutes.     

A three-dimensional RNA motif in Potato spindle tuber viroid mediates trafficking from palisade mesophyll to spongy mesophyll in Nicotiana benthamiana

Cell-to-cell trafficking of RNA is an emerging biological principle that integrates systemic gene regulation, viral infection, antiviral response, and cell-to-cell communication. A key mechanistic question is how an RNA is specifically selected for trafficking from one type of cell into another type. Here, we report the identification of an RNA motif in Potato spindle tuber viroid (PSTVd) required for trafficking from palisade mesophyll to spongy mesophyll in Nicotiana benthamiana leaves. This motif, called loop 6, has the sequence 5'-CGA-3'...5'-GAC-3' flanked on both sides by cis Watson-Crick G/C and G/U wobble base pairs. We present a three-dimensional (3D) structural model of loop 6 that specifies all non-Watson-Crick base pair interactions, derived by isostericity-based sequence comparisons with 3D RNA motifs from the RNA x-ray crystal structure database. The model is supported by available chemical modification patterns, natural sequence conservation/variations in PSTVd isolates and related species, and functional characterization of all possible mutants for each of the loop 6 base pairs. Our findings and approaches have broad implications for studying the 3D RNA structural motifs mediating trafficking of diverse RNA species across specific cellular boundaries and for studying the structure-function relationships of RNA motifs in other biological processes.     

The annotation of RNA motifs

The recent deluge of new RNA structures, including complete atomic-resolution views of both subunits of the ribosome, has on the one hand literally overwhelmed our individual abilities to comprehend the diversity of RNA structure, and on the other hand presented us with new opportunities for comprehensive use of RNA sequences for comparative genetic, evolutionary and phylogenetic studies. Two concepts are key to understanding RNA structure: hierarchical organization of global structure and isostericity of local interactions. Global structure changes extremely slowly, as it relies on conserved long-range tertiary interactions. Tertiary RNA-RNA and quaternary RNA-protein interactions are mediated by RNA motifs, defined as recurrent and ordered arrays of non-Watson-Crick base-pairs. A single RNA motif comprises a family of sequences, all of which can fold into the same three-dimensional structure and can mediate the same interaction(s). The chemistry and geometry of base pairing constrain the evolution of motifs in such a way that random mutations that occur within motifs are accepted or rejected insofar as they can mediate a similar ordered array of interactions. The steps involved in the analysis and annotation of RNA motifs in 3D structures are: (a) decomposition of each motif into non-Watson-Crick base-pairs; (b) geometric classification of each basepair; (c) identification of isosteric substitutions for each basepair by comparison to isostericity matrices; (d) alignment of homologous sequences using the isostericity matrices to identify corresponding positions in the crystal structure; (e) acceptance or rejection of the null hypothesis that the motif is conserved.     

Stacking geometry between two sheared Watson-Crick basepairs: Computational chemistry and bioinformatics based prediction

BackgroundMolecular modeling of RNA double helices is possible using most probable values of basepair parameters obtained from crystal structure database. The A:A w:wC non-canonical basepair, involving Watson-Crick edges of two Adenines in cis orientation, appears quite frequently in database. Bimodal distribution of its Shear, due to two different H-bonding schemes, introduces the confusion in assigning most the probable value. Its effect is pronounced when the A:A w:wC basepair stacks on Sheared wobble G:U W:WC basepairs.MethodsWe employed molecular dynamics simulations of three possible double helices with GAG, UAG and GAU sequence motifs at their centers and quantum chemical calculation for non-canonical A:A w:wC basepair stacked on G:U W:WC basepair.ResultsWe noticed stable structures of GAG motif with specifically negative Shear of the A:A basepair but stabilities of the other motifs were not found with A:A w:wC basepairing. Hybrid DFT-D and MP2 stacking energy analyses on dinucleotide step sequences, A:A w:wC::G:U W:WC and A:A w:wC::U:G W:WC reveal that viable orientation of A:A::G:U prefers one of the H-bonding modes with negative Shear, supported by crystal structure database. The A:A::U:G dinucleotide, however, prefers structure with only positive Shear.ConclusionsThe quantum chemical calculations explain why MD simulations of GAG sequence motif only appear stable. In the cases of the GAU and UAG motifs “tug of war” situation between positive and negative Shears of A:A w:wC basepair induces conformational plasticity.General significanceWe have projected comprehensive reason behind the promiscuous nature of A:A w:wC basepair which brings occasional structural plasticity.     

Attachment sites of primary binding proteins L1, L2 and L23 on 23 S ribosomal RNA of Escherichia coli

The attachment sites of the primary binding proteins L1, L2 and L23 on 23 S ribosomal RNA of Escherichia coli were examined by a chemical and ribonuclease footprinting method using several probes with different specificities. The results show that the sites are confined to localized RNA regions within the large ribonuclease-protected ribonucleoprotein fragments that were characterized earlier. They are as follows: (1) L1 recognizes a tertiary structural motif in domain V centred on two interacting internal loops; the main protein interaction sites occur at the internal loop/helix junctions. (2) The L2 site constitutes a single irregular stem/loop structure in the centre of domain IV where non-Watson-Crick pairing is likely to occur. (3) L23 recognizes a tertiary structural motif involving a single terminal loop structure and part of an adjacent internal loop at the centre of domain III. Each of the three primary binding proteins, whose presence is essential for ribosomal assembly, has been associated with important ribosomal functions: L1 lies in the E-site for deacylated tRNA binding while L2 and L23 have been implicated in the P and A substrate sites, respectively, of the peptidyl transferase centre. Moreover, each of the protein sites, but particularly those of L2 and L23, lies at the centre of RNA domains where they can maximally influence both the assembly of secondary binding proteins and the function of the RNA region.     

Architecture of the U5 small nuclear RNA.

We have used comparative sequence analysis and deletion analysis to examine the secondary structure of the U5 small nuclear RNA (snRNA), an essential component of the pre-mRNA splicing apparatus. The secondary structure of Saccharomyces cerevisiae U5 snRNA was studied in detail, while sequences from six other fungal species were included in the phylogenetic analysis. Our results indicate that fungal U5 snRNAs, like their counterparts from other taxa, can be folded into a secondary structure characterized by a highly conserved stem-loop (stem-loop 1) that is flanked by a moderately conserved internal loop (internal loop 1). In addition, several of the fungal U5 snRNAs include a novel stem-loop structure (ca. 30 nucleotides) that is adjacent to stem-loop 1. By deletion analysis of the S. cerevisiae snRNA, we have demonstrated that the minimal U5 snRNA that can complement the lethal phenotype of a U5 gene disruption consists of (i) stem-loop 1, (ii) internal loop 1, (iii) a stem-closing internal loop 1, and (iv) the conserved Sm protein binding site. Remarkably, all essential, U5-specific primary sequence elements are encoded by a 39-nucleotide domain consisting of stem-loop 1 and internal loop 1. This domain must, therefore, contain all U5-specific sequences that are essential for splicing activity, including binding sites for U5-specific proteins.     

Measuring the thermodynamics of RNA secondary structure formation

The thermodynamics of RNA secondary structure formation in small model systems provides a database for predicting RNA structure from sequence. Methods for making these measurements are reviewed with emphasis on optical methods and treatment of experimental errors. Analysis of experimental results in terms of simple nearest-neighbor models is presented. Some measured sequence dependences of non-Watson-Crick motifs are discussed. © 1998 John Wiley & Sons, Inc. Biopoly 44: 309–319, 1997     

Is there a special function for U.G basepairs in ribosomal RNA?

U.G basepairs are well-established elements of RNA structure. The geometry of this pair is different, however, from classical Watson-Crick basepairs. This leads to an unusual stacking of the basepair: overlap with the basepair at the 5' side of the U (and the 3' side of the G) is strong (stacked) while it is weak with the basepair on the other side (destacked). The closure of an RNA helix by a U.G pair will be energetically unfavourable when the U residue occupies the 5' end. In transfer RNA there is a strong selection against a 'destacked' U.G pair at helix ends. In the 16S rRNA model of Escherichia coli there are 72 U.G pairs of which 36 or 22 occupy a helix end, depending on how such an end is defined. There is a slight preference for 'stacked' U.G's in these positions. It is remarkable, however, that of 13 very conserved U.G pairs in the 16S (-like) rRNA, 7 occur at helix ends and that 5 of these have the 'destacked' configuration. It is suggested that these pairs, if they exist at all in a hydrogen-bounded form, are stabilized by co-axial stacking with other helices or by interaction with a protein.     

Analysis of RNA motifs

RNA motifs are directed and ordered stacked arrays of non-Watson-Crick base pairs forming distinctive foldings of the phosphodiester backbones of the interacting RNA strands. They correspond to the 'loops' - hairpin, internal and junction - that intersperse the Watson-Crick two-dimensional helices as seen in two-dimensional representations of RNA structure. RNA motifs mediate the specific interactions that induce the compact folding of complex RNAs. RNA motifs also constitute specific protein or ligand binding sites. A given motif is characterized by all the sequences that fold into essentially identical three-dimensional structures with the same ordered array of isosteric non-Watson-Crick base pairs. It is therefore crucial, when analyzing a three-dimensional RNA structure in order to identify and compare motifs, to first classify its non-Watson-Crick base pairs geometrically.     

Two distinct SECIS structures capable of directing selenocysteine incorporation in eukaryotes.

Translation of UGA as selenocysteine requires specific RNA secondary structures in the mRNAs of selenoproteins. These elements differ in sequence, structure, and location in the mRNA, that is, coding versus 3' untranslated region, in prokaryotes, eukaryotes, and archaea. Analyses of eukaryotic selenocysteine insertion sequence (SECIS) elements via computer folding programs, mutagenesis studies, and chemical and enzymatic probing has led to the derivation of a predicted consensus structural model for these elements. This model consists of a stem-loop or hairpin, with conserved nucleotides in the loop and in a non-Watson-Crick motif at the base of the stem. However, the sequences of a number of SECIS elements predict that they would diverge from the consensus structure in the loop region. Using site-directed mutagenesis to introduce mutations predicted to either disrupt or restore structure, or to manipulate loop size or stem length, we show that eukaryotic SECIS elements fall into two distinct classes, termed forms 1 and 2. Form 2 elements have additional secondary structures not present in form 1 elements. By either insertion or deletion of the sequences and structures distinguishing the two classes of elements while maintaining appropriate loop size, conversion of a form 1 element to a functional form 2-like element and of a form 2 to a functional form 1-like element was achieved. These results suggest commonality of function of the two classes. The information obtained regarding the existence of two classes of SECIS elements and the tolerances for manipulations of stem length and loop size should facilitate designing RNA molecules for obtaining high-resolution structural information about these elements.     

Solution structure of the pseudo-5' splice site of a retroviral splicing suppressor

Control of Rous sarcoma virus RNA splicing depends in part on the interaction of U1 and U11 snRNPs with an intronic RNA element called the negative regulator of splicing (NRS). A 23mer RNA hairpin (NRS23) of the NRS directly binds U1 and U11 snRNPs. Mutations that disrupt base-pairing between the loop of NRS23 and U1 snRNA abolish its negative control of splicing. We have determined the solution structure of NRS23 using NOEs, torsion angles, and residual dipolar couplings that were extracted from multidimensional heteronuclear NMR spectra. Our structure showed that the 6-bp stem of NRS23 adopts a nearly A-form duplex conformation. The loop, which consists of 11 residues according to secondary structure probing, was in a closed conformation. U913, the first residue in the loop, was bulged out or dynamic, and loop residues G914-C923, G915-U922, and U916-A921 were base-paired. The remaining UUGU tetraloop sequence did not adopt a stable structure and appears flexible in solution. This tetraloop differs from the well-known classes of tetraloops (GNRA, CUYG, UNCG) in terms of its stability, structure, and function. Deletion of the bulged U913, which is not complementary to U1 snRNA, increased the melting temperature of the RNA hairpin. This hyperstable hairpin exhibited a significant decrease in binding to U1 snRNP. Thus, the structure of the NRS RNA, as well as its sequence, is important for interaction with U1 snRNP and for splicing suppression.     

Identifying novel sequence variants of RNA 3D motifs

Predicting RNA 3D structure from sequence is a major challenge in biophysics. An important sub-goal is accurately identifying recurrent 3D motifs from RNA internal and hairpin loop sequences extracted from secondary structure (2D) diagrams. We have developed and validated new probabilistic models for 3D motif sequences based on hybrid Stochastic Context-Free Grammars and Markov Random Fields (SCFG/MRF). The SCFG/MRF models are constructed using atomic-resolution RNA 3D structures. To parameterize each model, we use all instances of each motif found in the RNA 3D Motif Atlas and annotations of pairwise nucleotide interactions generated by the FR3D software. Isostericity relations between non-Watson–Crick basepairs are used in scoring sequence variants. SCFG techniques model nested pairs and insertions, while MRF ideas handle crossing interactions and base triples. We use test sets of randomly-generated sequences to set acceptance and rejection thresholds for each motif group and thus control the false positive rate. Validation was carried out by comparing results for four motif groups to RMDetect. The software developed for sequence scoring (JAR3D) is structured to automatically incorporate new motifs as they accumulate in the RNA 3D Motif Atlas when new structures are solved and is available free for download.     

The selenocysteine incorporation machinery: interactions between the SECIS RNA and the SECIS-binding protein SBP2.

The decoding of UGA as a selenocysteine (Sec) codon in mammalian selenoprotein mRNAs requires a selenocysteine insertion sequence (SECIS) element in the 3' untranslated region. The SECIS is a hairpin structure that contains a non-Watson-Crick base-pair quartet with a conserved G.A/A.G tandem in the core of the upper helix. Another essential component of the Sec insertion machinery is SECIS-binding protein 2 (SBP2). In this study, we define the binding site of SBP2 on six different SECIS RNAs using enzymatic and hydroxyl radical footprinting, gel mobility shift analysis, and phosphate-ethylation binding interference. We show that SBP2 binds to a variety of mammalian SECIS elements with similar affinity and that the SBP2 binding site is conserved across species. Based on footprinting studies, SBP2 protects the proximal part of the hairpin and both strands of the lower half of the upper helix that contains the non-Watson-Crick base pair quartet. Gel mobility shift assays showed that the G.A/A.G tandem and internal loop are critical for the binding of SBP2. Modification of phosphates by ethylnitrosourea along both strands of the non-Watson-Crick base pair quartet, on the 5' strand of the lower helix and part of the 5' strand of the internal loop, prevented binding of SBP2. We propose a model in which SBP2 covers the central part of the SECIS RNA, binding to the non-Watson-Crick base pair quartet and to the 5' strands of the lower helix and internal loop. Our results suggest that the affinity of SBP2 for different SECIS elements is not responsible for the hierarchy of selenoprotein expression that is observed in vivo.     

A novel RNA structural motif in the selenocysteine insertion element of eukaryotic selenoprotein mRNAs.

In eukaryotes, co-translational insertion of selenocysteine into selenoproteins necessitates the participation of the selenocysteine insertion sequence (SECIS), an element lying in the 3'-untranslated region of selenoprotein mRNAs. We report a detailed experimental study of the secondary structures of the SECIS elements of three selenoprotein mRNAs, the rat and human type I iodothyronine deiodinase (5'DI) and rat glutathione peroxidase (GPx). Based on RNase and chemical probing, a new secondary structure model is established. It is characterized by a stem-loop structure, comprising two helices (I and II) separated by an internal loop, with an apical loop surmounting helix II. Sequence comparisons of 20 SECIS elements, arising from 2 5'DI, 13 GPx, 2 selenoprotein P, and 1 selenoprotein W mRNAs, confirm the secondary structure model. The most striking finding of the experimental study concerns a set of conserved sequences in helix II that interact to form a novel RNA structural motif consisting of a quartet composed of non-Watson-Crick base pairs 5'UGAY3': 5'UGAU3'. The potential for forming the quartet is preserved in 15 SECIS elements, but three consecutive non-Watson-Crick base pairs can nevertheless form in the other five SECIS, the central G.A tandem being invariant in all cases. A 3D model, derived by computer modeling with the use of the solution data, suggests that the base pairing interactions in the G.A tandem are of the type found in GNRA loops. The 3D model displays the quartet lying in an accessible position at the foot of helix II, which is bent at the internal loop, suggesting that the non-Watson-Crick base pair arrangement provides an unusual pattern of chemical groups for putative ligand interaction.     

RNA structural motifs: building blocks of a modular biomolecule

RNAs are modular biomolecules, composed largely of conserved structural subunits, or motifs. These structural motifs comprise the secondary structure of RNA and are knit together via tertiary interactions into a compact, functional, three-dimensional structure and are to be distinguished from motifs defined by sequence or function. A relatively small number of structural motifs are found repeatedly in RNA hairpin and internal loops, and are observed to be composed of a limited number of common 'structural elements'. In addition to secondary and tertiary structure motifs, there are functional motifs specific for certain biological roles and binding motifs that serve to complex metals or other ligands. Research is continuing into the identification and classification of RNA structural motifs and is being initiated to predict motifs from sequence, to trace their phylogenetic relationships and to use them as building blocks in RNA engineering.