
==== Front
Nucleic Acids ResNucleic Acids ResearchNucleic Acids Research0305-10481362-4962Oxford University Press 10.1093/nar/gkl96317135184Survey and SummaryThe interaction networks of structured RNAs Lescoute A. Westhof E. *Architecture et Réactivité de l'ARN, Université Louis PasteurIBMC, CNRS, 15 rue R.Descartes, F-67084 Strasbourg, France*To whom correspondence should be addressed. Tel/Fax: +33 388 41 70 46; Email : E.Westhof@ibmc.u-strasbg.fr12 2006 12 2006 28 11 2006 34 22 6587 6604 25 8 2006 21 10 2006 25 10 2006 © 2006 The Author(s).2006This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/2.0/uk/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

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.
==== Body
INTRODUCTION
RNA architecture is now visualized as the hierarchical assembly of preorganized double-stranded helices, formed by Watson–Crick base pairs, and which are connected and interlinked by RNA modules maintained by non-Watson–Crick base pairs. The folded and native architecture of a structured RNA molecule is so complex and intimate that each nucleotide forms several types of non-bonded interactions with the neighbouring nucleotides which are brought into contact by the folding process. Among the interactions, one can find: (i) phosphate–phosphate contacts mediated by water molecules or positively charged cations; (ii) phosphate–sugar H-bonding usually mediated by the 2′-hydroxyl group; (iii) sugar–sugar H-bonding interactions; (iv) base to phosphate or base to sugar H-bonding contacts; (v) base to phosphate or base to sugar stacking interactions; (vi) base to base H-bonding and (vii) base to base stacking interactions. Interactions (ii) to (vii) can occur either directly or indirectly via water molecules or ions (sometimes both modes occur simultaneously). The number of interacting strands generally varies between two and four. Since the polynucleotide sugar–phosphate backbone is monotonous, the specificity of the architectural fold and of the molecular recognition properties resides essentially in the base sequence, although as stated above several interactions either do not involve the bases or are sequence independent. In several instances, such contacts can be considered opportunistic or contingent to the presence of correct neighbouring base–base interactions. These molecular interactions do, however, contribute to the stability of the overall fold. Furthermore, the recent crystal structures clearly reveal the dominant role continuous base stacking plays in the maintenance of RNA architecture. It is, therefore, of central importance to represent and annotate base–base interactions as systematically and concisely as possible. Here, we present maps of the complex networks of interactions between nucleotide bases in known three-dimensional structures of RNAs. One major aim of the present analysis and survey is to facilitate the extraction and the recognition of the sequence constraints underlying a fold or three-dimensional motif. The hope is that the knowledge of such sequence constraints will help predictions of tertiary structure of RNAs on the sole basis of a genomic sequence.

The mapping is based on a nomenclature for base–base contacts that has been previously proposed (1). In such a scheme, any interacting assembly involving base–base interactions can be decomposed as a set of one-to-one base–base contacts. Using the proposed annotation symbols, one can therefore represent on a diagram the ensemble of all base–base contacts in a folded RNA, the canonical Watson–Crick base pairs of the secondary structure and the non-Watson–Crick pairs of the tertiary interactions. Here, we present planar diagrams for several types of structured RNAs. Whengenerating such diagrams, attention is paid to the coaxial stacking or parallel orientation of helices. Thus, such diagrams represent the network of base–base interactions occurring within a RNA architecture. Such diagrams cannot include all the overwhelming richness of any RNA architecture and, thus, cannot compensate for three-dimensional views of the structured RNAs. However, recurring motifs andthe networks they form are easily recognizable and local comparisons between structures facilitated (2). Such diagrams contain more information than that conveyed, with deceptive simplicity, by secondary structure drawings. The improvement in information content is achieved at a cost: (i) the need to know the symbolic nomenclature for the interactions; (ii) the transformation of the conventional secondary structure drawings into new ones displaying co-axial stacking of helices. Finally, whatever the extra information, the diagrams cannot replace three-dimensional views of the folded RNAs. On the contrary, they should be used as a complement to stereo views helping to organize and order the striking and amazing complexity of three-dimensional structures.

MATERIALS AND METHODS
RNA purine and pyrimidine bases present three edges for hydrogen bonding interactions (Figure 1): the Watson-Crick edge, the Hoogsteen edge (for purine) or the C–H edge (for pyrimidine), and the Sugar edge (which includes the sugar 2′-hydroxyl group) (1). A given edge of one base can potentially interact in a plane with any one of the three edges of a second base. The base–base interactions can occur in either the cis or trans orientation of the glycosyl bonds, i.e. with the attached sugars on the same side or on opposite sides of a line linking the interacting edges. This scheme leads to twelve possible, distinct edge-to-edge base-pairing geometries (or families). Each pairing geometry is designated by stating the interacting edges of the two bases (Watson–Crick represented by a circle, Hoogsteen represented by a square, or Sugar edge represented by a triangle) and the relative orientation of the glycosidic bonds, cis (filled symbols) or trans (empty symbols) (Figure 1). A historically based priority rule is invoked for listing the bases in a pair: Watson–Crick edge > Hoogsteen edge > Sugar edge.

This nomenclature allows for computer-based automatic search of base pairings (3–5) with planar drawings of the connected pairs (6). Here all contacts have been visually verified and all drawings done anew in Adobe Illustrator in order to emphasize their biological relevance as described above. Table 1 gives the list of the structures for which we present the diagram of the network of interactions. The selection spans the sizes of the small ribozymes to the large 16S rRNA of the 30S ribosomal subunit. An editable version of all drawings can be retrieved from the laboratory website ().

The evolving view of a structured RNA
The secondary structures of most structured RNAs with a defined and precise function occur in many flavours. For example, group I introns are classified in at least four sub-groups characterized by the number and positions of helices attached to an invariant secondary structure core (7). Despite this variability, the overall architectures of the ribozymes promote the stabilization of the helical stems building up the core and the correct positioning of the helical substrates [e.g. in group I introns, see Refs. (8–13)]. This is mainly achieved by the properties of the RNA anchoring motifs which allow for the formation of different and often mutually exclusive long–range contacts between non-homologous peripheral elements. Recurrent and systematic use of essentially two main types of long-range RNA–RNA anchors are observed in domain assembly: GNRA tetraloops with their receptors (14,15) and loop–loop Watson–Crick or non-Watson–Crick base pairings (8,15,16). In order to promote such long-range contacts, sub-domains, which are usually subtended by complex and diverse sets of molecular interactions, have to be assembled. Among them, three-way junctions constitute frequent and critical structured sub-domains necessary to promote further long-range RNA–RNA contacts. For example, the three-way junction that forms the catalytic core of the hammerhead ribozyme is constrained by tertiary interactions between peripheral elements. These constraints accelerate the folding of the ribozyme which is hundred-times more efficient than the minimally reduced ribozyme (17–19).

Group I introns
Figure 2 illustrates the evolving process of our structural understanding of a large RNA molecule, the prototypical group I introns. First, the secondary structure elements were established using sequence alignments (20) and a planar structure diagram agreed upon (21). After modelling of the three-dimensional structure (7), a new secondary structure diagram was proposed (22). The latter representation (22) included numerous crossings of strands but allowed to visualize the co-axial stacking of helices and the relative position of the co-axial stacks with respect to each other. On the basis of such an architecture of the core, long-range and unsuspected contacts between peripheral elements could be established and modelled (8). Finally, upon completion of crystal structures at sufficient resolution, full interaction diagrams could be produced (9–11).

Supplementary Figure S1 shows the interaction network diagrams for the three crystal structures of group I introns published to date. A comparison between those diagrams shows quickly several points. The numerous A-minor contacts [e.g. J3/4-P6 (23)] and especially those made by the GNRA tetraloops with the 11nt-receptor or two stacked helical base pairs are conspicuous. The 11nt-receptor displays variability in the conformation of the bulging U (compare J8/8a and J5/5a in the Azoarcus intron). Differences between the various subgroups appear also. L9 always contact the end of P5, but P9 folds back on P5 through distinct local topology (a three-way junction in Twort; a ‘reverse Kink-turn’ in Azoarcus (24); an undefined bend in Tetrahymena). Depending on the P1/P10 substrate, J8/7 displays different conformations; some nucleotides are in the syn conformation in the absence of substrate and new A-minor contacts occur in presence of substrate.

The RNase P structures
A similar evolving process is apparent in the RNase P structures (Figure 3). First, phylogenetic analysis led to the establishment of the secondary structures (25) for the two families (types A and B) (26). Later, modelling studies (27,28) led to the predictions of co-axial stacking between definite helices in coherence with long-range contacts. This formed the basis for a second type of drawings. Recently, crystal structures of the specificity domains (29,30) and of the catalytic domains (31,32) were published.

A striking result of the two structures of the specificity domain was the fold of the large L11/12 internal loop which displayed two T-like loops (33) characterized by a 5-membered loop closed by a trans Watson–Crick/Hoogsteen base pair. Despite the variability between the two structures, in part due to the different junctions, some unexpected contacts are maintained, e.g. the last residue of L11/12, C175 in type A, stacks with a purine residue following the first T-like loop while it is G168, part of a bacterial loop E motif in P10.1 of type B (28), which forms the same stack. This is not the only example (see Figure 4): P11 contains a stack of two base triples made of a C=G and A–U Watson–Crick pair to which, respectively, a A forms a Sugar/Hoogsteen cis and an amino base a Watson–Crick/Hoogsteen trans. But, in type A, the amino base is a A which comes from an internal loop 5′ of P11 and, in type B, it is a C which is a single-stranded residue linking the 3′end of P10 to P10.1. However, the A forming the triple with the C=G base pair is structurally identical in both types. Further, in both types a bulging A from stem P9 stacks 5′ of the A of the A–U pair. These structures stress the role of bulging residues for linking two helical domains (e.g. bulged A in P9 linking P9 and P11 or the last residue of L13 which stacks 5 base pairs below on a bulged G coming from the co-axially stacked P14) (34). An almost identical situation occurs in helix h11 of both 16S rRNA structures (35,36) (see below). For completeness, the catalytic domain (32) and the full ribozyme (31) have been also drawn but a somewhat lower resolution precludes comparisons between these two structures (see Supplementary Figure S2). The multibranch junction around P1–P4 is particularly intriguing and present some analogies with the active site of the 23S rRNA (see below).

The Diels-Adler and nucleolytic ribozymes
In Figure 5 are shown the interaction diagrams for the published crystal structures of small ribozymes. The Diels-Adler ribozyme (Figure 5a) (37) displays three non-Watson–Crick pairs and a two base pair pseudoknot at the 5′ end (38). The hepatitis delta virus ribozyme presents a pseudoknotted helix at its 3′ end (39) and an additional two base pair pseudoknot below the cleavage site (Figure 5b) (40,41). This latter pseudoknot is stabilized by two base triples. The single-strand leading to the 3′ end pseudoknot forms A-minor contacts with a co-axially stacked helix (40,42). It is striking to see how four nucleotides of the 6-membered hairpin loop L3 form the stable UNCG tetraloop characterized by a G in the syn conformation (43). The hairpin ribozyme (Figure 5c) has a four-way junction with interactions mediated by structured internal loops between two of the helices (44). The ribose zipper (between nucleotides 10/11 and 24/25) critical for the folding (45,46) is clearly apparent. The loop E like structure of loop B (46,47) displays the characteristic set of non-Watson–Crick pairs. It is interesting to notice, however, that the bulging U42 forms a cis Watson–Crick pair with A23.

The two main crystal structures of the hammerhead ribozymes are shown on Figure 5d and 5e (48,49). Recently, it has been clear that the structure shown on Figure 5d does not correspond to the active structure of the hammerhead ribozyme. The new crystal structure of the hammerhead ribozyme from Schistosome mansoni (49) is stunning in several respects. Interestingly, the loop–loop interactions are dominated by trans Watson–Crick/Hoogsteen pairs. A comparison between the structures of the hairpin and hammerhead ribozymes reveals an intriguing pattern (Figure 5f). In both cases, cleavage occurs between two trans Hoogsteen/Sugar–edge pairs (tandem sheared base pairs) with an intermediate residue forming a canonical Watson–Crick G=C pair.

The riboswitches
Five riboswitch structures have now been determined, all complexed with activating ligand (Figure 6): two structures of the guanine riboswitch (50,51), the SAM riboswitch (52) and two structures of the thiamine pyrophosphate (TPP) riboswitch (53,54). Each of those structures contains at leastone motif found recurrently in other RNAs: a common three-way junction in the guanine riboswitch (Figure 6a) (55,56); a kink-turn motif (57,58) in the SAM riboswitch (Figure 6b); a T-like loop (33) in the TPP riboswitch (Figure 6c). As in the Schistosome hammerhead, the trans Watson–Crick/Hoogsteen pairs are frequent in the loop–loop interactions of the guanine riboswitch. Interestingly, the latter type of base pairs leads to a local antiparallel orientation of the strands. Although the SAM riboswitch displays a standard A-minor type of contact with alternating type I and type II pairs (59), the TPP riboswitch presents an example of the unfrequent double type I/type I contact (Figure 7) (60).

The transpeptidylation center of the 23S rRNA structure
The transpeptidylation center is a multibranch loop in domain V of the 23S rRNA (61). Two interaction diagrams are shown on Figures 8 and 9. Figure 8 corresponds to the structure in the Haloarcula marismortui 23S rRNA and shows also the contacts with the -CCA end of a tRNA in the A site and P site (62). Figure 9 corresponds to the E.coli structure (36). To help the understanding, the usual secondary structure and a 3D representation are shown in Figure 10. The compactness of the PTC multibranch loop leads to a T-like shape with the parallel stacks of H89 and H90/H91 forming the vertical branch of the ‘T’, helices H73 and H93 forming the horizontal branch of the ‘T’ and the active nucleotides at the center of the ‘T’. Helices H90/H91/H92 adopt the frequent family C fold of three-way junction (56), which brings the apical loop of H92 close to the important residues U2506 and U2585 thereby forming part of the binding site for the A site tRNA on the 23S rRNA. The center contains three occurrences (helices H73, H74 and H91) of a motif described above (Figure 4) which is present in P10/P11 of RNase P RNAs. One of those motifs includes A2451 of the active site (H74) and forms part of the binding site of P site tRNA on the 23S rRNA.

The 16S rRNA structures
The interaction network diagrams for the 16S rRNA of T.thermophilus and E.coli are shown in Figures 11 and 12. The characteristic feature of the 16S rRNA is the clear decomposition into large domains (35): the 5′, the central, the 3′major and the 3′minor domains. The 5′domain forms the ‘body’, the 3′major domain forms the ‘head’ and the central domain forms the ‘platform’ (63,64). In the interaction network diagram, the 5′pseudoknot (formed between the loop of h1 and the junction between helices h27 and h28) is set at the center because it interconnects the 5′domain, the 3′major and the central domain. The 3′major domain has pronounced intra-domain connections and only two inter-domain connections (h36 with the 5′ pseudoknot and h28, which is stacked with the 5′pseudoknot, interacts with the base of helix h44 in the 3′minor domain). A single primary binding protein (65–67), S7, interacts in the multibranch junctions which follow helix h28. On the contrary, the central domain has less intra-domain contacts but several inter-domain contacts: one with the 3′minor (h24 with h45); two with the 5′pseudoknot (h26a and h27), and four with the 5′domain (h19 with h12, h27 with h11 and the long helix h21 with h4 and h12). Two primary binding proteins bind to the central domain: proteins S8 and S15. The 5′domain presents rather extensive intra-domain contacts, has two primary binding proteins S4 and S17 binding to it, and besides the four inter-domain contacts with the central domain presents additional links with the long decoding helix h44 (via h13 in both structures and also via h8 in the E.coli structure). In short, the head is rather free to move with respect to the other domains, while the platform and body could move in concert with the movements of the body transmitted to the penultimate helix h44. These long-range contacts therefore reflect the known conformational changes occurring in the ribosome (9,36,68,69).

The E.coli structure has revealed the molecular details of three important RNA–RNA inter-subunit contacts (36,69–71). The inter-subunit bridge, B7a, occurs between h23 and H68 (see Figure 13). The inter-subunit bridge, B2a, occurs between the A site of h44 and the loop of H69, while the inter-subunit bridge, B3, occurs below the A site of h44 and with helix H71 (see Figure 14). The B7a bridge is essentially a cis sugar-edge/sugar-edge base pair which interestingly involves an unusual K-turn in the 16S [KT-23, see (58)]. The B2a bridge involves the A site of the 16S and especially an invariant UoU pair. Although the conformation of that UoU pair in the free state is not precisely known, in presence of antibiotics of the aminoglycoside family, the conformation of the UoU pair is always bifurcated with the O4 of U1406 pointing to the N3-H of U1495 (72,73). Thus, the effects of antibiotic binding could be transmitted to the large subunit via the bridge B2a. The B3 bridge exploits the less frequent A-minor motif with two type I contacts (see Figure 7) between a tandem trans Hoogsteen/sugar-edge GoA pairs (sheared) and two stacked G=C pairs of H71.

Structural consequences of modularity
The modular architecture of folded RNAs implies that distances between interacting parts are conserved in functionally homologous molecules. Thus, elongator tRNAs have generally four base pairs in the dihydrouridine helix and five in the thymine helix and the maintenance of those distances allows for specific base pairing contacts between a conserved set of two Gs in the dihydrouridine loop and the conserved TΨCG segment of the thymine loop. Also, in group I introns (see Figures 2 and Supplementary Figure S1), the distance between the GoU pair at the cleavage site of helix P1 and the GAAA capping loop of the co-axially stacked P2 helix is always of 12 bp. This invariant distance guarantees formation of the contact between GAAA and its receptor in P8 (3bp below the interface with helix P3) (7).

The 16S rRNA presents similar examples of conserved numbered of base pairs between interacting modules. The distance between the two sets of interacting As in helix h21 is 7 bp (with two unpaired nucleotides); the two sets of As contact identical positions of base pairs in h4 and h12. In helix h44, the distance between the motif involved in bridge B3 (tandem trans Hoogsteen/sugar-edge GoA pairs followed by a GoU pair) and the motif forming an inter-domain contact with h13 is of 12 bp (despite the fact that the motif of h44 interacting with h13 is not conserved). In helix h17, the distance between the interface with h16 and the A contacting a G=C pair in h15 is 13 bp (not counting bulged residues); at the same time, the distance between the same G=C pair in h15 and the end of the co-axially stacked h4 is 11 bp. Both helices h22 and h26 start with a trans Hoogsteen/sugar-edge base pair and the sixth base pair above it in h22 forms a sugar-edge/sugar-edge base pair with the eighth base pair in h26 of both 16S rRNA structures (the overall length of h26 is variable).

The capping GNRA loop of h8 forms an identical contact with the last base pair of h14 and 6 bp from the capping GNRA a striking motif in h8 forms identical contacts with an identical motif at the base of h6. However, beyond that conserved motif, the length of h6 is variable. In Thermus thermophilus 16S rRNA, helix h9 is longer than in E.coli 16S rRNA, but there is an additional residue in h7 which forms a trans Watson–Crick/sugar-edge within the capping loop of h9. At the same time, helix h10 is shorter in T.thermophilus, which allows an unpaired U of h8 to form a trans Hoogsteen/sugar-edge with the first base pair of h10. There is an interdomain contact between the tip of h44 and the beginning of h8 in E.coli which is absent in T.thermophilus where h44 is too short to form such contacts. Such conservations in distances and contacts are significant considering the distance between the two bacterial species to which these 16S rRNAs belong to.

Similarities in motifs and contacts made are apparent from such diagrams. For example, the striking motif in h8, forming a tight contact with h6 (74), is found also in h41 where it forms also a tight contact with h43. Helices h20 and h24 both contain an identical motif of three consecutive non-Watson–Crick pairs in trans (Hoogsteen/sugar-edge, Watson–Crick/sugar-edge). Surprisingly, they form no further contacts but, two base pairs above, two Watson–Crick pairs are engaged in A-minor interactions (in which the As originate from different types of motifs in h23a and h23b, respectively).

CONCLUSIONS
The wiring diagrams representing the networks of all pairwise base–base interactions in architectures of structured RNAs illustrate several principles of RNA folding. A striking observation is the ubiquitous presence of RNA motifs with a defined order of non-Watson–Crick base pairs. They are found recurrently in structured RNAs with different functions and stemming from the three kingdoms of life. Some motifs are more frequent and more pervasive than others (e.g. loop E or K-turn motifs). For those motifs, we have progressed in the establishment of their sequence constraints (58,75). Although a complete survey of RNA motifs is not yet achieved, their number is clearly limited (76). The task of enumerating and cataloguing RNA motifs is problematic because of the lack of accepted definition on the nature and size of RNA motifs (2,77). This is compounded by the observation that most RNA motifs can be assembled like Russian dolls such that only portions of RNA motifs can occur (78,79). For example, the variants of the T-loop motifs (33,80,81) or those of the loop E family (79,82).

In statistical physics, complex systems are analysed as networks of interactions between the components of the system. Recently, network theory has permeated all realms of science (83,84). Some major real-life networks have striking properties. They are said ‘scale-free’, meaning that their properties are controlled by a small number of highly connected nodes. In such scale-free networks, any two nodes can be connected by a very small number of intermediate nodes, the small-world effect. The main attribute of such networks is the presence of clusters of local interactions with long-range interactions between the clusters (85,86). Visually, in the wiring diagrams of the 16S rRNAs, the clusters are apparent as well as the connections between them. Such diagrams appear like hierarchical networks which have the properties of scale-free networks with embedded modularity (84,87). However, as shown for protein structures (88–90), in the process of conceptualising macromolecules as networks several simplifications and hypotheses have to be formulated and evaluated. How should the number and types of contacts be considered? Should we consider all atomic contacts or only hydrogen bonds? Such analysis has not been performed yet for RNA systems. The hypothesis that RNA three-dimensional structures form scale-free networks could lead to valuable knowledge on the evolution pathways of RNA molecules. The origin of scale-free topology is rooted in the self-organization of networks because networks grow by preferential linking and attachment of links to nodes with a high number of connections (86). The reverse process could therefore lead to the set of nucleotides constituting the common ancestor of 16S rRNA.

SUPPLEMENTARY DATA
Supplementary Data are available at NAR Online.

Supplementary Material
[Supplementary Data]
 We are grateful to Neocles Leontis for numerous discussions and sharing of data. Funding to pay the Open Access publication charges for this article was provided by The Human Frontier Science Program (grant RGP0032/2005-C to E.W.).

Conflict of interest statement. None declared.

Figures and Tables
Figure 1 The three edges of the purine and pyrimidine nucleotides (Watson–Crick, Hoogsteen and Sugar) with the two possible orientations of the attached glycoside bongs. They can be cis or trans depending on whether they are on the same side (cis) or on opposite sides (trans) of a line median to the base-base H-bonds. Redrawn after (1).The symbols used for annotating the diagrams are shown below. A nucleotide with a base in the syn orientation with respect to the sugar is represented either as a hollowed or as a bold uppercase letter. The annotation for stacking (empty rectangle used as a placeholder) was introduced by (9). In some instances where a single base–base H-bond is present, the connected atoms are given and separated by a dashed line or the symbols are separated by a dashed line, in order to help visualize the facing edges. Within the text, the non-Watson–Crick base pairs is indicated by an ‘o’ between the interacting nucleotides instead of the symbols and when necessary the full name is given. There is no special symbol for the backbone connectivity; the line or the arrow is chosen so as to prevent confusion with other base-base symbols.

Figure 2 Evolution of the structural representations of group I introns. (A) Drawing based on the analysis of sequences (21). (B) Drawing with the inclusion of the co-axial stacking based on three-dimensional modelling (7,22). In (A) and (B), the numbers in circles correspond to the domain numbers. (C) Drawing with annotated tertiary contacts based on the crystal structure of the Tetrahymena thermophyla group I intron (11).

Figure 3 Evolution of the structural representations of the RNA of the RNase P ribozyme. On the left for the type A subgroup and on the right for the type B subgroup. At the top, the representation based on the phylogeny with some non-Watson–Crick base pairs shown as filled circles. (26). In the middle, the representation following the 3D modelling (28). The double ended red arrows correspond to the long-range contacts between secondary structure elements. Below, the annotated tertiary contacts based on the crystal structures (29,30).

Figure 4 Consensus diagram of a motif observed in P RNA in stem P11 and in the transpeptidylation center of the 23S rRNA.

Figure 5 Interaction network diagrams for (A) the Diels-Adler ribozyme (37); (B) the hepatitis delta virus ribozyme (40); (C) the hairpin ribozyme (44); (D) the minimal hammerhead ribozyme (48); (e) the Schistosome hammerhead ribozyme (49); (F) comparison of the active sites of the hairpin and hammerhead ribozymes. The single ended red arrows indicate the cleavage site. The colored nucleotides indicate those that change state between the minimal and complete hammerhead ribozymes.

Figure 6 Interaction network diagrams for (A) the guanine riboswitch in presence of hypoxanthine (HX) (51) and on the insert at the right the variations in presnce of guanine and adenine (50); (B) the SAM riboswitch (52); (C) the thiamine pyrophosphate riboswitch (54). The filled red circles indicate those bases contacted by the ligand in each case and the dashed red circles indicate van der Waals contacts (52).

Figure 7 Example of the type I and type II A-minor contacts (59) and their observed associations into motifs. The most frequently observed one is typeII/typeI from the 5′ end to the 3′ end. The double combination typeI/typeI does occur also but it is less frequent.

Figure 8 Interaction network diagrams for the transpeptidylation center of the 23S rRNA of H. marismortui (62) (69,94) showing the contacts with the –CCA ends of A- and P-site tRNAs. The circles with a central dot indicates that the helix points towards the reader and the circle with a cross the reverse. Known modifications are indicated by blue stars (95).

Figure 9 Interaction network diagrams for the transpeptidylation center of the 23S rRNA of E.coli (36). Known modifications are indicated by blue stars (96).

Figure 10 The conventional secondary structure and the tertiary fold of the transpeptidylation center of H.marismortui (62). The lowercase letters indicate nucleotides which are modified.

Figure 11 Interaction network diagrams for the 16S rRNA of T.thermophilus in the 30S subunit (35). The contacts of the primary proteins (S4, S7, S8, S15, S17) are colored on the contacted nucleotides following the indicated color code (97).

Figure 12 Interaction network diagrams for the for the 16S rRNA of E.coli in the 70S ribosome (36). The three inter-subunit bridges RNA mediated are shown in blue (36). The modifications are shown by stars (in orange: 2′-O-methylations; green: pseudouridines) (98).

Figure 13 The central domain of the 16S rRNA in the E.coli ribosome (36) showing the tertiary fold, the conventional secondary structure, the view emphasizing co-axial stacking and the network diagram.

Figure 14 The 3′minor domain of the 16S rRNA in the E.coli ribosome (36) showing the tertiary fold, the conventional secondary structure, the view emphasizing co-axial stacking and the network diagram. For comparison, the 3′minor domain of the 16S rRNA in the T.thermophilus 30S particle (35) is also shown.

Table 1 List of structures with the Protein Data Bank (93) identification code and the crystallographic resolution

Structure	Organism	PDB id	Resolution (Å)	References	
Group I intron	with exons	Azoarcus	1U6B	3.10	(9)	
	ribozyme	Tetrahymena	1X8W	3.80	(11)	
	ribozyme/product	Twort	1Y0Q	3.60	(10)	
RNaseP Type A RNA	Spec. Domain	T.thermophilus	1US9	2.90	(33)	
	Full length	T.maritima	2A2E	3.85	(32)	
RNaseP Type B RNA	Spec. Domain	B. subtilis	1NBS	3.15	(29)	
	Full length	B.stearothermophilus	2A64	3.30	(31)	
16S rRNA		T.thermophilus	1J5E	3.05	(35)	
		E.coli	2AVY	3.46	(36)	
23S rRNA	Transpeptidation site	H.marismortui	1S72	2.40	(91)	
	TPC + CCA	H.marismortui	1Q86	3.00	(92)	
	TPC	E.coli	2AW4	3.46	(36)	
Hepatitis delta virus			1DRZ	2.30	(40)	
Diels-Alder ribozyme		Aptamer	1YLS	3.00	(37)	
Hairpin ribozyme		Viroïd RNA	1HP6	2.40	(44)	
Hammerhead ribozyme	All RNA		1MME	3.10	(48)	
	Whole	S.mansoni	2G0Z	2.20	(49)	
Guanine riboswitch	With hypoxanthine	B.subtilis	1U8D	1.95	(51)	
	With guanine		1Y27	2.40	(50)	
	With adenine		1Y26	2.10	(50)	
TPP riboswitch		E.coli	2GDI	2.05	(54)	
SAM riboswitch		T.tengcongensis	2GIS	2.90	(52)
==== Refs
REFERENCES
1 Leontis N.B. Westhof E.  Geometric nomenclature and classification of RNA base pairs RNA 2001 7 499 512 11345429 
2 Leontis N.B. Lescoute A. Westhof E.  The building blocks and motifs of RNA architecture Curr. Opin. Struct. Biol. 2006 16 279 287 16713707 
3 Yang H. Jossinet F. Leontis N. Chen L. Westbrook J. Berman H. Westhof E.  Tools for the automatic identification and classification of RNA base pairs Nucleic Acids Res. 2003 31 3450 3460 12824344 
4 Gendron P. Lemieux S. Major F.  Quantitative analysis of nucleic acid three-dimensional structures J. Mol. Biol. 2001 308 919 936 11352582 
5 Lemieux S. Major F.  RNA canonical and non-canonical base pairing types: a recognition method and complete repertoire Nucleic Acids Res. 2002 30 4250 4263 12364604 
6 Jossinet F. Westhof E.  Sequence to Structure (S2S): display, manipulate and interconnect RNA data from sequence to structure Bioinformatics 2005 21 3320 3321 15905274 
7 Michel F. Westhof E.  Modelling of the three-dimensional architecture of group I catalytic introns based on comparative sequence analysis J. Mol. Biol. 1990 216 585 610 2258934 
8 Lehnert V. Jaeger L. Michel F. Westhof E.  New loop-loop tertiary interactions in self-splicing introns of subgroup IC and ID: a complete 3D model of the Tetrahymena thermophila ribozyme Chem. Biol. 1996 3 993 1009 9000010 
9 Adams P.L. Stahley M.R. Gill M.L. Kosek A.B. Wang J. Strobel S.A.  Crystal structure of a group I intron splicing intermediate RNA 2004 10 1867 1887 15547134 
10 Golden B.L. Kim H. Chase E.  Crystal structure of a phage Twort group I ribozyme-product complex Nature Struct. Mol. Biol. 2005 12 82 89 15580277 
11 Guo F. Gooding A.R. Cech T.R.  Structure of the Tetrahymena ribozyme: base triple sandwich and metal ion at the active site Mol. Cell 2004 16 351 362 15525509 
12 Woodson S.A.  Structure and assembly of group I introns Curr. Opin. Struct. Biol. 2005 15 324 330 15922592 
13 Vicens Q. Cech T.R.  Atomic level architecture of group I introns revealed Trends Biochem. Sci. 2006 31 41 51 16356725 
14 Costa M. Michel F.  Frequent use of the same tertiary motif by self-folding RNAs EMBO J. 1995 14 1276 1285 7720718 
15 Costa M. Michel F.  Rules for RNA recognition of GNRA tetraloops deduced by in vitro selection: comparison with in vivo evolution EMBO J. 1997 16 3289 3302 9214644 
16 Costa M. Michel F. Westhof E.  A three-dimensional perspective on exon binding by a group II self-splicing intron EMBO J. 2000 19 5007 5018 10990464 
17 Khvorova A. Lescoute A. Westhof E. Jayasena S.D.  Sequence elements outside the hammerhead ribozyme catalytic core enable intracellular activity Nature Struct. Biol. 2003 10 708 712 12881719 
18 Penedo J.C. Wilson T.J. Jayasena S.D. Khvorova A. Lilley D.M.  Folding of the natural hammerhead ribozyme is enhanced by interaction of auxiliary elements RNA 2004 10 880 888 15100442 
19 Canny M.D. Jucker F.M. Kellogg E. Khvorova A. Jayasena S.D. Pardi A.  Fast cleavage kinetics of a natural hammerhead ribozyme J. Am. Chem. Soc. 2004 126 10848 10849 15339162 
20 Michel F. Jacquier A. Dujon B.  Comparison of fungal mitochondrial introns reveals extensive homologies in RNA secondary structure Biochimie 1982 64 867 881 6817818 
21 Burke J.M. Belfort M. Cech T.R. Davies R.W. Schweyen R.J. Shub D.A. Szostak J.W. Tabak H.F.  Structural conventions for group I introns Nucleic Acids Res. 1987 15 7217 7221 3658691 
22 Cech T.R. Damberger S.H. Gutell R.R.  Representation of the secondary and tertiary structure of group I introns Nature Struct. Biol. 1994 1 273 280 7545072 
23 Michel F. Ellington A.D. Couture S. Szostak J.W.  Phylogenetic and genetic evidence for base-triples in the catalytic domain of group I introns Nature 1990 347 578 580 2215683 
24 Strobel S.A. Adams P.L. Stahley M.R. Wang J.  RNA kink turns to the left and to the right RNA 2004 10 1852 1854 15547133 
25 James B.D. Olsen G.J. Liu J. Pace N.R.  The secondary structure of ribonuclease P RNA, the catalytic element of a ribonucleoprotein enzyme Cell 1988 52 19 26 2449969 
26 Brown J.W.  The Ribonuclease P Database Nucleic Acids Res. 1999 27 314 9847214 
27 Harris M.E. Kazantchev A.V. Chen J.L. Pace N.R.  Analysis of the tertiary structure of the ribonuclease P ribozyme-substrate complex by site-specific photoaffinity crosslinking RNA 1997 3 561 574 9174092 
28 Massire C. Jaeger L. Westhof E.  Derivation of the three-dimensional architecture of bacterial ribonuclease P RNAs from comparative sequence analysis J. Mol. Biol. 1998 279 773 793 9642060 
29 Krasilnikov A.S. Xiao Y. Pan T. Mondragon A.  Basis for structural diversity in homologous RNAs Science 2004 306 104 107 15459389 
30 Krasilnikov A.S. Yang X. Pan T. Mondragon A.  Crystal structure of the specificity domain of ribonuclease P Nature 2003 421 760 764 12610630 
31 Kazantsev A.V. Krivenko A.A. Harrington D.J. Holbrook S.R. Adams P.D. Pace N.R.  Crystal structure of a bacterial ribonuclease P RNA Proc. Natl Acad. Sci. USA 2005 102 13392 13397 16157868 
32 Torres-Larios A. Swinger K.K. Krasilnikov A.S. Pan T. Mondragon A.  Crystal structure of the RNA component of bacterial ribonuclease P Nature 2005 437 584 587 16113684 
33 Krasilnikov A.S. Mondragon A.  On the occurrence of the T-loop RNA folding motif in large RNA molecules RNA 2003 9 640 643 12756321 
34 Hermann T. Patel D.J.  RNA bulges as architectural and recognition motifs Structure 2000 8 R47 R54 10745015 
35 Wimberly B.T. Brodersen D.E. Clemons W.M. JrMorgan-Warren R.J. Carter A.P. Vonrhein C. Hartsch T. Ramakrishnan V.  Structure of the 30S ribosomal subunit Nature 2000 407 327 339 11014182 
36 Schuwirth B.S. Borovinskaya M.A. Hau C.W. Zhang W. Vila-Sanjurjo A. Holton J.M. Cate J.H.  Structures of the bacterial ribosome at 3.5 Å resolution Science 2005 310 827 834 16272117 
37 Serganov A. Keiper S. Malinina L. Tereshko V. Skripkin E. Hobartner C. Polonskaia A. Phan A.T. Wombacher R. Micura R.  Structural basis for Diels-Alder ribozyme-catalyzed carbon-carbon bond formation Nature Struct. Mol. Biol. 2005 12 218 224 15723077 
38 Keiper S. Bebenroth D. Seelig B. Westhof E. Jaschke A.  Architecture of a Diels-Alderase ribozyme with a preformed catalytic pocket Chem. Biol. 2004 11 1217 1227 15380182 
39 Perrotta A.T. Been M.D.  A pseudoknot-like structure required for efficient self-cleavage of hepatitis delta virus RNA Nature 1991 350 434 436 2011192 
40 Ferre-D'Amare A.R. Zhou K. Doudna J.A.  Crystal structure of a hepatitis delta virus ribozyme Nature 1998 395 567 574 9783582 
41 Wadkins T.S. Perrotta A.T. Ferre-D'Amare A.R. Doudna J.A. Been M.D.  A nested double pseudoknot is required for self-cleavage activity of both the genomic and antigenomic hepatitis delta virus ribozymes RNA 1999 5 720 727 10376872 
42 Tanner N.K. Thill G. Petit-Kostas E. Crain-Denoyelle A.M. Westhof E.  A three-dimensional model of hepatitis delta virus ribozyme based on biochemical and mutational analyses Curr. Biol. 1994 4 488 498 7922369 
43 Ennifar E. Nikulin A. Tishchenko S. Serganov A. Nevskaya N. Garber M. Ehresmann B. Ehresmann C. Nikonov S. Dumas P.  The crystal structure of UUCG tetraloop J. Mol. Biol. 2000 304 35 42 11071808 
44 Rupert P.B. Ferre-D'Amare A.R.  Crystal structure of a hairpin ribozyme-inhibitor complex with implications for catalysis Nature 2001 410 780 786 11298439 
45 Chowrira B.M. Berzal-Herranz A. Keller C.F. Burke J.M.  Four ribose 2′-hydroxyl groups essential for catalytic function of the hairpin ribozyme J. Biol. Chem. 1993 268 19458 19462 7690032 
46 Earnshaw D.J. Masquida B. Müller S. Sigurdsson S.T. Eckstein F. Westhof E. Gait M.J.  Inter-domain cross-linking and molecular modelling of the hairpin ribozyme J. Mol. Biol. 1997 274 197 212 9398527 
47 Butcher S.E. Allain F.H.-T. Feigon J.  Solution structure of the loop B from the hairpin ribozyme Nature Struct. Biol. 1999 6 212 216 10074938 
48 Scott W.G. Finch J.T. Klug A.  The crystal structure of an all-RNA hammerhead ribozyme Nucleic Acids Symp. Ser. 1995 34 214 216 8841628 
49 Martick M. Scott W.G.  Tertiary contacts distant from the active site prime a ribozyme for catalysis Cell 2006 126 309 320 16859740 
50 Serganov A. Yuan Y.R. Pikovskaya O. Polonskaia A. Malinina L. Phan A.T. Hobartner C. Micura R. Breaker R.R. Patel D.J.  Structural basis for discriminative regulation of gene expression by adenine- and guanine-sensing mRNAs Chem. Biol. 2004 11 1729 1741 15610857 
51 Batey R.T. Gilbert S.D. Montange R.K.  Structure of a natural guanine-responsive riboswitch complexed with the metabolite hypoxanthine Nature 2004 432 411 415 15549109 
52 Montange R.K. Batey R.T.  Structure of the S-adenosylmethionine riboswitch regulatory mRNA element Nature 2006 441 1172 1175 16810258 
53 Thore S. Leibundgut M. Ban N.  Structure of the eukaryotic thiamine pyrophosphate riboswitch with its regulatory ligand Science 2006 312 1208 1211 16675665 
54 Serganov A. Polonskaia A. Phan A.T. Breaker R.R. Patel D.J.  Structural basis for gene regulation by a thiamine pyrophosphate-sensing riboswitch Nature 2006 441 1167 1171 16728979 
55 Lescoute A. Westhof E.  Riboswitch structures: purine ligands replace tertiary contacts Chem. Biol. 2005 12 10 13 15664510 
56 Lescoute A. Westhof E.  Topology of three-way junctions in folded RNAs RNA 2006 12 83 93 16373494 
57 Klein D.J. Schmeing T.M. Moore P.B. Steitz T.A.  The kink-turn: a new RNA secondary structure motif EMBO J. 2001 20 4214 4221 11483524 
58 Lescoute A. Leontis N.B. Massire C. Westhof E.  Recurrent structural RNA motifs, Isostericity Matrices and sequence alignments Nucleic Acids Res. 2005 33 2395 2409 15860776 
59 Nissen P. Ippolito J.A. Ban N. Moore P.B. Steitz T.A.  RNA tertiary interactions in the large ribosomal subunit: the A-minor motif Proc. Natl Acad. Sci. USA 2001 98 4899 4903 11296253 
60 Lescoute A. Westhof E.  The A-minor motifs in the decoding recognition process Biochimie 2006 88 993 999 16889885 
61 Barta A. Steiner G. Brosius J. Noller H.F. Kuechler E.  Identification of a site on 23S ribosomal RNA located at the peptidyl transferase center Proc. Natl Acad. Sci. USA 1984 81 3607 3611 6374660 
62 Ban N. Nissen P. Hansen J. Moore P.B. Steitz T.A.  The complete atomic structure of the large ribosomal subunit at 2.4 Å resolution Science 2000 289 905 920 10937989 
63 Lata K.R. Agrawal R.K. Penczek P. Grassucci R. Zhu J. Frank J.  Three-dimensional reconstruction of the Escherichia coli 30S ribosomal subunit in ice J. Mol. Biol. 1996 262 43 52 8809178 
64 Clemons W.M. JrMay J.L. Wimberly B.T. McCutcheon J.P. Capel M.S. Ramakrishnan V.  Structure of a bacterial 30S ribosomal subunit at 5.5 Å resolution Nature 1999 400 833 840 10476960 
65 Mizushima S. Nomura M.  Assembly mapping of 30S ribosomal proteins from E.coli Nature 1970 226 1214 4912319 
66 Nomura M. Erdmann V.A.  Reconstitution of 50S ribosomal subunits from dissociated molecular components Nature 1970 228 744 748 5472962 
67 Grondek J.F. Culver G.M.  Assembly of the 30S ribosomal subunit: positioning ribosomal protein S13 in the S7 assembly branch RNA 2004 10 1861 1866 15525707 
68 Serdyuk I. Baranov V. Tsalkova T. Gulyamova D. Pavlov M. Spirin A. May R.  Structural dynamics of translating ribosomes Biochimie 1992 74 299 306 1379074 
69 Yusupov M.M. Yusupova G.Z. Baucom A. Lieberman K. Earnest T.N. Cate J.H. Noller H.F.  Crystal structure of the ribosome at 5.5 A resolution Science 2001 292 883 896 11283358 
70 Frank J. Verschoor A. Li Y. Zhu J. Lata R.K. Radermacher M. Penczek P. Grassucci R. Agrawal R.K. Srivastava S.  A model of the translational apparatus based on a three-dimensional reconstruction of the Escherichia coli ribosome Biochem. Cell Biol. 1995 73 757 765 8721992 
71 Cate J.H. Yusupov M.M. Yusupova G.Z. Earnest T.N. Noller H.F.  X-ray crystal structures of 70S ribosome functional complexes Science 1999 285 2095 2104 10497122 
72 Vicens Q. Westhof E.  Crystal structure of paromomycin docked into the eubacterial ribosomal decoding A site Structure 2001 9 647 658 11587639 
73 Francois B. Russell R.J. Murray J.B. Aboul-ela F. Masquida B. Vicens Q. Westhof E.  Crystal structures of complexes between aminoglycosides and decoding A site oligonucleotides: role of the number of rings and positive charges in the specific binding leading to miscoding Nucleic Acids Res. 2005 33 5677 5690 16214802 
74 Noller H.F.  RNA structure: reading the ribosome Science 2005 309 1508 1514 16141058 
75 Leontis N.B. Westhof E.  A common motif organizes the structure of multi-helix loops in 16S and 23S ribosomal RNAs J. Mol. Biol. 1998 283 571 583 9784367 
76 Moore P.B.  Structural motifs in RNA Annu. Rev. Biochem. 1999 68 287 300 10872451 
77 Leontis N.B. Altman R.B. Berman H.M. Brenner S.E. Brown J.W. Engelke D.R. Harvey S.C. Holbrook S.R. Jossinet F. Lewis S.E.  The RNA Ontology Consortium: an open invitation to the RNA community RNA 2006 12 533 541 16484377 
78 Westhof E. Fritsch V.  RNA folding: beyond Watson–Crick pairs Structure 2000 8 R55 R65 10745012 
79 Leontis N.B. Stombaugh J. Westhof E.  Motif prediction in ribosomal RNAs lessons and prospects for automated motif prediction in homologous RNA molecules Biochimie 2002 84 961 973 12458088 
80 Leontis N.B. Westhof E.  Analysis of RNA motifs Curr. Opin. Struct. Biol. 2003 13 300 308 12831880 
81 Lee J.C. Cannone J.J. Gutell R.R.  The lonepair triloop: a new motif in RNA structure J. Mol. Biol. 2003 325 65 83 12473452 
82 Lee J.C. Gutell R.R. Russell R.  The UAA/GAN internal loop motif: a new RNA structural element that forms a cross-strand AAA stack and long-range tertiary interactions J. Mol. Biol. 2006 360 978 988 16828489 
83 Strogatz S.H.  Exploring complex networks Nature 2001 410 268 276 11258382 
84 Barabasi A.L. Oltvai Z.N.  Network biology: understanding the cell's functional organization Nature Rev. Genet. 2004 5 101 113 14735121 
85 Watts D.J. Strogatz S.H.  Collective dynamics of ‘small-world’ networks Nature 1998 393 440 442 9623998 
86 Barabasi A.L. Albert R.  Emergence of scaling in random networks Science 1999 286 509 512 10521342 
87 Ravasz E. Somera A.L. Mongru D.A. Oltvai Z.N. Barabasi A.L.  Hierarchical organization of modularity in metabolic networks Science 2002 297 1551 1555 12202830 
88 Greene L.H. Higman V.A.  Uncovering network systems within protein structures J. Mol. Biol. 2003 334 781 791 14636602 
89 Amitai G. Shemesh A. Sitbon E. Shklar M. Netanely D. Venger I. Pietrokovski S.  Network analysis of protein structures identifies functional residues J. Mol. Biol. 2004 344 1135 1146 15544817 
90 del Sol A. Fujihashi H. Amoros D. Nussinov R.  Residues crucial for maintaining short paths in network communication mediate signaling in proteins Mol. Syst. Biol. 2006 2 2006. 0019 
91 Klein D.J. Moore P.B. Steitz T.A.  The contribution of metal ions to the structural stability of the large ribosomal subunit RNA 2004 10 1366 1379 15317974 
92 Hansen J.L. Schmeing T.M. Moore P.B. Steitz T.A.  Structural insights into peptide bond formation Proc. Natl Acad. Sci. USA 2002 99 11670 11675 12185246 
93 Berman H.M. Westbrook J. Feng Z. Gilliland G. Bhat T.N. Weissig H. Shindyalov I.N. Bourne P.E.  The Protein Data Bank Nucleic Acids Res. 2000 28 235 242 10592235 
94 Nissen P. Hansen J. Ban N. Moore P.B. Steitz T.A.  The structural basis of ribosome activity in peptide bond synthesis Science 2000 289 920 930 10937990 
95 Mengel-Jorgensen J. Jensen S.S. Rasmussen A. Poehlsgaard J. Iversen J.J. Kirpekar F.  Modifications in Thermus thermophilus 23S ribosomal RNA are centered in regions of RNA-RNA contact J. Biol. Chem. 2006 281 22108 22117 16731530 
96 Kowalak J.A. Bruenger E. McCloskey J.A.  Posttranscriptional modification of the central loop of domain V in Escherichia coli 23S ribosomal RNA J. Biol. Chem. 1995 270 17758 17764 7629075 
97 Brodersen D.E. Clemons W.M. JrCarter A.P. Wimberly B.T. Ramakrishnan V.  Crystal structure of the 30S ribosomal subunit from Thermus thermophilus: structure of the proteins and their interactions with 16S RNA J. Mol. Biol. 2002 316 725 768 11866529 
98 McCloskey J.A. Rozenski J.  The small subunit rRNA modification database Nucleic Acids Res. 2005 33 D135 D138 15608163

