Phage anti-CRISPR control by an RNA- and DNA-binding helix-turn-helix protein.

In all organisms, regulation of gene expression must be adjusted to meet cellular requirements and frequently involves helix-turn-helix (HTH) domain proteins1. For instance, in the arms race between bacteria and bacteriophages, rapid expression of phage anti-CRISPR (acr) genes upon infection enables evasion from CRISPR-Cas defence; transcription is then repressed by an HTH-domain-containing anti-CRISPR-associated (Aca) protein, probably to reduce fitness costs from excessive expression2-5. However, how a single HTH regulator adjusts anti-CRISPR production to cope with increasing phage genome copies and accumulating acr mRNA is unknown. Here we show that the HTH domain of the regulator Aca2, in addition to repressing Acr synthesis transcriptionally through DNA binding, inhibits translation of mRNAs by binding conserved RNA stem-loops and blocking ribosome access. The cryo-electron microscopy structure of the approximately 40 kDa Aca2-RNA complex demonstrates how the versatile HTH domain specifically discriminates RNA from DNA binding sites. These combined regulatory modes are widespread in the Aca2 family and facilitate CRISPR-Cas inhibition in the face of rapid phage DNA replication without toxic acr overexpression. Given the ubiquity of HTH-domain-containing proteins, it is anticipated that many more of them elicit regulatory control by dual DNA and RNA binding.

Structural basis of exoribonuclease-mediated mRNA transcription termination.

Efficient termination is required for robust gene transcription. Eukaryotic organisms use a conserved exoribonuclease-mediated mechanism to terminate the mRNA transcription by RNA polymerase II (Pol II)1-5. Here we report two cryogenic electron microscopy structures of Saccharomyces cerevisiae Pol II pre-termination transcription complexes bound to the 5'-to-3' exoribonuclease Rat1 and its partner Rai1. Our structures show that Rat1 displaces the elongation factor Spt5 to dock at the Pol II stalk domain. Rat1 shields the RNA exit channel of Pol II, guides the nascent RNA towards its active centre and stacks three nucleotides at the 5' terminus of the nascent RNA. The structures further show that Rat1 rotates towards Pol II as it shortens RNA. Our results provide the structural mechanism for the Rat1-mediated termination of mRNA transcription by Pol II in yeast and the exoribonuclease-mediated termination of mRNA transcription in other eukaryotes.

Cryo-EM structure of SRP68/72 reveals an extended dimerization domain with RNA-binding activity.

The signal recognition particle (SRP) is a critical component in protein sorting pathways in all domains of life. Human SRP contains six proteins bound to the 7S RNA and their structures and functions have been mostly elucidated. The SRP68/72 dimer is the largest SRP component and is essential for SRP function. Although the structures of the SRP68/72 RNA binding and dimerization domains have been previously reported, the structure and function of large portions of the SRP68/72 dimer remain unknown. Here, we analyse full-length SRP68/72 using cryo-EM and report that SRP68/72 depend on each other for stability and form an extended dimerization domain. This newly observed dimerization domain is both a protein- and RNA-binding domain. Comparative analysis with current structural models suggests that this dimerization domain undergoes dramatic translocation upon SRP docking onto SRP receptor and eventually comes close to the Alu domain. We propose that the SRP68/72 dimerization domain functions by binding and detaching the Alu domain and SRP9/14 from the ribosomal surface, thus releasing elongation arrest upon docking onto the ER membrane.

Structure of the Dicer-2-R2D2 heterodimer bound to a small RNA duplex.

In flies, Argonaute2 (Ago2) and small interfering RNA (siRNA) form an RNA-induced silencing complex to repress viral transcripts1. The RNase III enzyme Dicer-2 associates with its partner protein R2D2 and cleaves long double-stranded RNAs to produce 21-nucleotide siRNA duplexes, which are then loaded into Ago2 in a defined orientation2-5. Here we report cryo-electron microscopy structures of the Dicer-2-R2D2 and Dicer-2-R2D2-siRNA complexes. R2D2 interacts with the helicase domain and the central linker of Dicer-2 to inhibit the promiscuous processing of microRNA precursors by Dicer-2. Notably, our structure represents the strand-selection state in the siRNA-loading process, and reveals that R2D2 asymmetrically recognizes the end of the siRNA duplex with the higher base-pairing stability, and the other end is exposed to the solvent and is accessible by Ago2. Our findings explain how R2D2 senses the thermodynamic asymmetry of the siRNA and facilitates the siRNA loading into Ago2 in a defined orientation, thereby determining which strand of the siRNA duplex is used by Ago2 as the guide strand for target silencing.

Structural insights into dsRNA processing by Drosophila Dicer-2-Loqs-PD.

Small interfering RNAs (siRNAs) are the key components for RNA interference (RNAi), a conserved RNA-silencing mechanism in many eukaryotes1,2. In Drosophila, an RNase III enzyme Dicer-2 (Dcr-2), aided by its cofactor Loquacious-PD (Loqs-PD), has an important role in generating 21 bp siRNA duplexes from long double-stranded RNAs (dsRNAs)3,4. ATP hydrolysis by the helicase domain of Dcr-2 is critical to the successful processing of a long dsRNA into consecutive siRNA duplexes5,6. Here we report the cryo-electron microscopy structures of Dcr-2-Loqs-PD in the apo state and in multiple states in which it is processing a 50 bp dsRNA substrate. The structures elucidated interactions between Dcr-2 and Loqs-PD, and substantial conformational changes of Dcr-2 during a dsRNA-processing cycle. The N-terminal helicase and domain of unknown function 283 (DUF283) domains undergo conformational changes after initial dsRNA binding, forming an ATP-binding pocket and a 5'-phosphate-binding pocket. The overall conformation of Dcr-2-Loqs-PD is relatively rigid during translocating along the dsRNA in the presence of ATP, whereas the interactions between the DUF283 and RIIIDb domains prevent non-specific cleavage during translocation by blocking the access of dsRNA to the RNase active centre. Additional ATP-dependent conformational changes are required to form an active dicing state and precisely cleave the dsRNA into a 21 bp siRNA duplex as confirmed by the structure in the post-dicing state. Collectively, this study revealed the molecular mechanism for the full cycle of ATP-dependent dsRNA processing by Dcr-2-Loqs-PD.

Structure and dynamics of the quaternary hunchback mRNA translation repression complex.

A key regulatory process during Drosophila development is the localized suppression of the hunchback mRNA translation at the posterior, which gives rise to a hunchback gradient governing the formation of the anterior-posterior body axis. This suppression is achieved by a concerted action of Brain Tumour (Brat), Pumilio (Pum) and Nanos. Each protein is necessary for proper Drosophila development. The RNA contacts have been elucidated for the proteins individually in several atomic-resolution structures. However, the interplay of all three proteins during RNA suppression remains a long-standing open question. Here, we characterize the quaternary complex of the RNA-binding domains of Brat, Pum and Nanos with hunchback mRNA by combining NMR spectroscopy, SANS/SAXS, XL/MS with MD simulations and ITC assays. The quaternary hunchback mRNA suppression complex comprising the RNA binding domains is flexible with unoccupied nucleotides functioning as a flexible linker between the Brat and Pum-Nanos moieties of the complex. Moreover, the presence of the Pum-HD/Nanos-ZnF complex has no effect on the equilibrium RNA binding affinity of the Brat RNA binding domain. This is in accordance with previous studies, which showed that Brat can suppress mRNA independently and is distributed uniformly throughout the embryo.

Opposing roles for Egalitarian and Staufen in transport, anchoring and localization of oskar mRNA in the Drosophila oocyte.

Localization of oskar mRNA includes two distinct phases: transport from nurse cells to the oocyte, a process typically accompanied by cortical anchoring in the oocyte, followed by posterior localization within the oocyte. Signals within the oskar 3' UTR directing transport are individually weak, a feature previously hypothesized to facilitate exchange between the different localization machineries. We show that alteration of the SL2a stem-loop structure containing the oskar transport and anchoring signal (TAS) removes an inhibitory effect such that in vitro binding by the RNA transport factor, Egalitarian, is elevated as is in vivo transport from the nurse cells into the oocyte. Cortical anchoring within the oocyte is also enhanced, interfering with posterior localization. We also show that mutation of Staufen recognized structures (SRSs), predicted binding sites for Staufen, disrupts posterior localization of oskar mRNA just as in staufen mutants. Two SRSs in SL2a, one overlapping the Egalitarian binding site, are inferred to mediate Staufen-dependent inhibition of TAS anchoring activity, thereby promoting posterior localization. The other three SRSs in the oskar 3' UTR are also required for posterior localization, including two located distant from any known transport signal. Staufen, thus, plays multiple roles in localization of oskar mRNA.

Cryo-EM structure of the yeast TREX complex and coordination with the SR-like protein Gbp2.

The evolutionarily conserved TRanscript-EXport (TREX) complex plays central roles during mRNP (messenger ribonucleoprotein) maturation and export from the nucleus to the cytoplasm. In yeast, TREX is composed of the THO sub-complex (Tho2, Hpr1, Tex1, Mft1, and Thp2), the DEAD box ATPase Sub2, and Yra1. Here we present a 3.7 Šcryo-EM structure of the yeast THO•Sub2 complex. The structure reveals the intimate assembly of THO revolving around its largest subunit Tho2. THO stabilizes a semi-open conformation of the Sub2 ATPase via interactions with Tho2. We show that THO interacts with the serine-arginine (SR)-like protein Gbp2 through both the RS domain and RRM domains of Gbp2. Cross-linking mass spectrometry analysis supports the extensive interactions between THO and Gbp2, further revealing that RRM domains of Gbp2 are in close proximity to the C-terminal domain of Tho2. We propose that THO serves as a landing pad to configure Gbp2 to facilitate its loading onto mRNP.

Structural design principles for specific ultra-high affinity interactions between colicins/pyocins and immunity proteins.

The interactions of the antibiotic proteins colicins/pyocins with immunity proteins is a seminal model system for studying protein-protein interactions and specificity. Yet, a precise and quantitative determination of which structural elements and residues determine their binding affinity and specificity is still lacking. Here, we used comparative structure-based energy calculations to map residues that substantially contribute to interactions across native and engineered complexes of colicins/pyocins and immunity proteins. We show that the immunity protein α1-α2 motif is a unique structurally-dissimilar element that restricts interaction specificity towards all colicins/pyocins, in both engineered and native complexes. This motif combines with a diverse and extensive array of electrostatic/polar interactions that enable the exquisite specificity that characterizes these interactions while achieving ultra-high affinity. Surprisingly, the divergence of these contributing colicin residues is reciprocal to residue conservation in immunity proteins. The structurally-dissimilar immunity protein α1-α2 motif is recognized by divergent colicins similarly, while the conserved immunity protein α3 helix interacts with diverse colicin residues. Electrostatics thus plays a key role in setting interaction specificity across all colicins and immunity proteins. Our analysis and resulting residue-level maps illuminate the molecular basis for these protein-protein interactions, with implications for drug development and rational engineering of these interfaces.

Recognition of non-CpG repeats in Alu and ribosomal RNAs by the Z-RNA binding domain of ADAR1 induces A-Z junctions.

Adenosine-to-inosine (A-to-I) editing of eukaryotic cellular RNAs is essential for protection against auto-immune disorders. Editing is carried out by ADAR1, whose innate immune response-specific cytoplasmic isoform possesses a Z-DNA binding domain (Zα) of unknown function. Zα also binds to CpG repeats in RNA, which are a hallmark of Z-RNA formation. Unexpectedly, Zα has been predicted - and in some cases even shown - to bind to specific regions within mRNA and rRNA devoid of such repeats. Here, we use NMR, circular dichroism, and other biophysical approaches to demonstrate and characterize the binding of Zα to mRNA and rRNA fragments. Our results reveal a broad range of RNA sequences that bind to Zα and adopt Z-RNA conformations. Binding is accompanied by destabilization of neighboring A-form regions which is similar in character to what has been observed for B-Z-DNA junctions. The binding of Zα to non-CpG sequences is specific, cooperative and occurs with an affinity in the low micromolar range. This work allows us to propose a model for how Zα could influence the RNA binding specificity of ADAR1.

Structural basis for the final steps of human 40S ribosome maturation.

Eukaryotic ribosomes consist of a small 40S and a large 60S subunit that are assembled in a highly coordinated manner. More than 200 factors ensure correct modification, processing and folding of ribosomal RNA and the timely incorporation of ribosomal proteins1,2. Small subunit maturation ends in the cytosol, when the final rRNA precursor, 18S-E, is cleaved at site 3 by the endonuclease NOB13. Previous structures of human 40S precursors have shown that NOB1 is kept in an inactive state by its partner PNO14. The final maturation events, including the activation of NOB1 for the decisive rRNA-cleavage step and the mechanisms driving the dissociation of the last biogenesis factors have, however, remained unresolved. Here we report five cryo-electron microscopy structures of human 40S subunit precursors, which describe the compositional and conformational progression during the final steps of 40S assembly. Our structures explain the central role of RIOK1 in the displacement and dissociation of PNO1, which in turn allows conformational changes and activation of the endonuclease NOB1. In addition, we observe two factors, eukaryotic translation initiation factor 1A domain-containing protein (EIF1AD) and leucine-rich repeat-containing protein 47 (LRRC47), which bind to late pre-40S particles near RIOK1 and the central rRNA helix 44. Finally, functional data shows that EIF1AD is required for efficient assembly factor recycling and 18S-E processing. Our results thus enable a detailed understanding of the last steps in 40S formation in human cells and, in addition, provide evidence for principal differences in small ribosomal subunit formation between humans and the model organism Saccharomyces cerevisiae.

Host ANP32A mediates the assembly of the influenza virus replicase.

Aquatic birds represent a vast reservoir from which new pandemic influenza A viruses can emerge1. Influenza viruses contain a negative-sense segmented RNA genome that is transcribed and replicated by the viral heterotrimeric RNA polymerase (FluPol) in the context of viral ribonucleoprotein complexes2,3. RNA polymerases of avian influenza A viruses (FluPolA) replicate viral RNA inefficiently in human cells because of species-specific differences in acidic nuclear phosphoprotein 32 (ANP32), a family of essential host proteins for FluPol activity4. Host-adaptive mutations, particularly a glutamic-acid-to-lysine mutation at amino acid residue 627 (E627K) in the 627 domain of the PB2 subunit, enable avian FluPolA to overcome this restriction and efficiently replicate viral RNA in the presence of human ANP32 proteins. However, the molecular mechanisms of genome replication and the interplay with ANP32 proteins remain largely unknown. Here we report cryo-electron microscopy structures of influenza C virus polymerase (FluPolC) in complex with human and chicken ANP32A. In both structures, two FluPolC molecules form an asymmetric dimer bridged by the N-terminal leucine-rich repeat domain of ANP32A. The C-terminal low-complexity acidic region of ANP32A inserts between the two juxtaposed PB2 627 domains of the asymmetric FluPolA dimer, suggesting a mechanism for how the adaptive PB2(E627K) mutation enables the replication of viral RNA in mammalian hosts. We propose that this complex represents a replication platform for the viral RNA genome, in which one of the FluPol molecules acts as a replicase while the other initiates the assembly of the nascent replication product into a viral ribonucleoprotein complex.

Structure of the human core transcription-export complex reveals a hub for multivalent interactions.

The export of mRNA from nucleus to cytoplasm requires the conserved and essential transcription and export (TREX) complex (THO-UAP56/DDX39B-ALYREF). TREX selectively binds mRNA maturation marks and licenses mRNA for nuclear export by loading the export factor NXF1-NXT1. How TREX integrates these marks and achieves high selectivity for mature mRNA is poorly understood. Here, we report the cryo-electron microscopy structure of the human THO-UAP56/DDX39B complex at 3.3 Å resolution. The seven-subunit THO-UAP56/DDX39B complex multimerizes into a 28-subunit tetrameric assembly, suggesting that selective recognition of mature mRNA is facilitated by the simultaneous sensing of multiple, spatially distant mRNA regions and maturation marks. Two UAP56/DDX39B RNA helicases are juxtaposed at each end of the tetramer, which would allow one bivalent ALYREF protein to bridge adjacent helicases and regulate the TREX-mRNA interaction. Our structural and biochemical results suggest a conserved model for TREX complex function that depends on multivalent interactions between proteins and mRNA.

The DNA of human and other eukaryotic cells is stored inside a compartment called the nucleus. DNA carries the genetic code and provides a blueprint for all of the cell's proteins. However, protein production occurs outside the nucleus, in the main body of the cell. To transmit genetic information from one compartment to the other, the DNA sequences are first transcribed into another molecule called messenger RNA, or mRNA for short. Once made, mRNA exits the nucleus and enters the cell's main body to encounter the machinery that translates its sequence into a protein. Before mRNA can exit the nucleus, it must first undergo a series of modifications, which result in the mRNA molecule being successively bound to specific proteins. Once mRNA has passed through these steps, it is recognized by the transcription-and-export complex, or TREX for short, which is comprised of several proteins. When TREX binds to mRNA, it adds on a final protein which allows the mRNA molecule to be transported out of the nucleus. However, it remained unclear how TREX selects the completed mRNA-protein complexes that are ready for export while at the same time recognizing the wide variety of mRNA molecules produced by cells. Now, Pühringer and Hohmann et al. have identified the first three-dimensional structure of the core of the human TREX complex using a technique called cryo-electron microscopy. This revealed that the seven proteins of the TREX core assemble into a large complex that has four copies of each protein. The structure suggests that TREX can bind to mRNA and its attached proteins in various ways. These different binding arrangements may help the complex select which mRNA molecules are fully modified and ready to be exported. The structure also sheds light on how mutations in this complex can lead to diseases such as Beaulieu­Boycott­Innes syndrome (BBIS). This work will help guide future research into the activity of TREX, including how its structure changes when it binds to mRNA and deposits the final transport protein. Identifying these structures will make it easier to design experiments that target specific aspects of TREX activity and provide new insights into how these complexes work.

The TFIIH subunits p44/p62 act as a damage sensor during nucleotide excision repair.

Nucleotide excision repair (NER) in eukaryotes is orchestrated by the core form of the general transcription factor TFIIH, containing the helicases XPB, XPD and five 'structural' subunits, p62, p44, p34, p52 and p8. Recent cryo-EM structures show that p62 makes extensive contacts with p44 and in part occupies XPD's DNA binding site. While p44 is known to regulate the helicase activity of XPD during NER, p62 is thought to be purely structural. Here, using helicase and adenosine triphosphatase assays we show that a complex containing p44 and p62 enhances XPD's affinity for dsDNA 3-fold over p44 alone. Remarkably, the relative affinity is further increased to 60-fold by dsDNA damage. Direct binding studies show this preference derives from p44/p62's high affinity (20 nM) for damaged ssDNA. Single molecule imaging of p44/p62 complexes without XPD reveals they bind to and randomly diffuse on DNA, however, in the presence of UV-induced DNA lesions these complexes stall. Combined with the analysis of a recent cryo-EM structure, we suggest that p44/p62 acts as a novel DNA-binding entity that enhances damage recognition in TFIIH. This revises our understanding of TFIIH and prompts investigation into the core subunits for an active role during DNA repair and/or transcription.

Dual conformational recognition by Z-DNA binding protein is important for the B-Z transition process.

Left-handed Z-DNA is radically different from the most common right-handed B-DNA and can be stabilized by interactions with the Zα domain, which is found in a group of proteins, such as human ADAR1 and viral E3L proteins. It is well-known that most Zα domains bind to Z-DNA in a conformation-specific manner and induce rapid B-Z transition in physiological conditions. Although many structural and biochemical studies have identified the detailed interactions between the Zα domain and Z-DNA, little is known about the molecular basis of the B-Z transition process. In this study, we successfully converted the B-Z transition-defective Zα domain, vvZαE3L, into a B-Z converter by improving B-DNA binding ability, suggesting that B-DNA binding is involved in the B-Z transition. In addition, we engineered the canonical B-DNA binding protein GH5 into a Zα-like protein having both Z-DNA binding and B-Z transition activities by introducing Z-DNA interacting residues. Crystal structures of these mutants of vvZαE3L and GH5 complexed with Z-DNA confirmed the significance of conserved Z-DNA binding interactions. Altogether, our results provide molecular insight into how Zα domains obtain unusual conformational specificity and induce the B-Z transition.

Periphilin self-association underpins epigenetic silencing by the HUSH complex.

Transcription of integrated DNA from viruses or transposable elements is tightly regulated to prevent pathogenesis. The Human Silencing Hub (HUSH), composed of Periphilin, TASOR and MPP8, silences transcriptionally active viral and endogenous transgenes. HUSH recruits effectors that alter the epigenetic landscape and chromatin structure, but how HUSH recognizes target loci and represses their expression remains unclear. We identify the physicochemical properties of Periphilin necessary for HUSH assembly and silencing. A disordered N-terminal domain (NTD) and structured C-terminal domain are essential for silencing. A crystal structure of the Periphilin-TASOR minimal core complex shows Periphilin forms an α-helical homodimer, bound by a single TASOR molecule. The NTD forms insoluble aggregates through an arginine/tyrosine-rich sequence reminiscent of low-complexity regions from self-associating RNA-binding proteins. Residues required for TASOR binding and aggregation were required for HUSH-dependent silencing and genome-wide deposition of repressive mark H3K9me3. The NTD was functionally complemented by low-complexity regions from certain RNA-binding proteins and proteins that form condensates or fibrils. Our work suggests the associative properties of Periphilin promote HUSH aggregation at target loci.

Mechanism of ribosome rescue by alternative ribosome-rescue factor B.

Alternative ribosome-rescue factor B (ArfB) rescues ribosomes stalled on non-stop mRNAs by releasing the nascent polypeptide from the peptidyl-tRNA. By rapid kinetics we show that ArfB selects ribosomes stalled on short truncated mRNAs, rather than on longer mRNAs mimicking pausing on rare codon clusters. In combination with cryo-electron microscopy we dissect the multistep rescue pathway of ArfB, which first binds to ribosomes very rapidly regardless of the mRNA length. The selectivity for shorter mRNAs arises from the subsequent slow engagement step, as it requires longer mRNA to shift to enable ArfB binding. Engagement results in specific interactions of the ArfB C-terminal domain with the mRNA entry channel, which activates peptidyl-tRNA hydrolysis by the N-terminal domain. These data reveal how protein dynamics translate into specificity of substrate recognition and provide insights into the action of a putative rescue factor in mitochondria.

CryoEM structure of the low-complexity domain of hnRNPA2 and its conversion to pathogenic amyloid.

hnRNPA2 is a human ribonucleoprotein (RNP) involved in RNA metabolism. It forms fibrils both under cellular stress and in mutated form in neurodegenerative conditions. Previous work established that the C-terminal low-complexity domain (LCD) of hnRNPA2 fibrillizes under stress, and missense mutations in this domain are found in the disease multisystem proteinopathy (MSP). However, little is known at the atomic level about the hnRNPA2 LCD structure that is involved in those processes and how disease mutations cause structural change. Here we present the cryo-electron microscopy (cryoEM) structure of the hnRNPA2 LCD fibril core and demonstrate its capability to form a reversible hydrogel in vitro containing amyloid-like fibrils. Whereas these fibrils, like pathogenic amyloid, are formed from protein chains stacked into ß-sheets by backbone hydrogen bonds, they display distinct structural differences: the chains are kinked, enabling non-covalent cross-linking of fibrils and disfavoring formation of pathogenic steric zippers. Both reversibility and energetic calculations suggest these fibrils are less stable than pathogenic amyloid. Moreover, the crystal structure of the disease-mutation-containing segment (D290V) of hnRNPA2 suggests that the replacement fundamentally alters the fibril structure to a more stable energetic state. These findings illuminate how molecular interactions promote protein fibril networks and how mutation can transform fibril structure from functional to a pathogenic form.

Structural snapshots of human pre-60S ribosomal particles before and after nuclear export.

Ribosome biogenesis is an elaborate and energetically expensive program that involve two hundred protein factors in eukaryotes. Nuclear export of pre-ribosomal particles is one central step which also serves as an internal structural checkpoint to ensure the proper completion of nuclear assembly events. Here we present four structures of human pre-60S particles isolated through a nuclear export factor NMD3, representing assembly stages immediately before and after nuclear export. These structures reveal locations of a dozen of human factors, including an uncharacterized factor TMA16 localized between the 5S RNA and the P0 stalk. Comparison of these structures shows a progressive maturation for the functional regions, such as peptidyl transferase centre and peptide exit tunnel, and illustrate a sequence of factor-assisted rRNA maturation events. These data facilitate our understanding of the global conservation of ribosome assembly in eukaryotes and species-specific features of human assembly factors.

Base-intercalated and base-wedged stacking elements in 3D-structure of RNA and RNA-protein complexes.

Along with nucleobase pairing, base-base stacking interactions are one of the two main types of strong non-covalent interactions that define the unique secondary and tertiary structure of RNA. In this paper we studied two subfamilies of nucleobase-inserted stacking structures: (i) with any base intercalated between neighboring nucleotide residues (base-intercalated element, BIE, i + 1); (ii) with any base wedged into a hydrophobic cavity formed by heterocyclic bases of two nucleotides which are one nucleotide apart in sequence (base-wedged element, BWE, i + 2). We have exploited the growing database of natively folded RNA structures in Protein Data Bank to analyze the distribution and structural role of these motifs in RNA. We found that these structural elements initially found in yeast tRNAPhe are quite widespread among the tertiary structures of various RNAs. These motifs perform diverse roles in RNA 3D structure formation and its maintenance. They contribute to the folding of RNA bulges and loops and participate in long-range interactions of single-stranded stretches within RNA macromolecules. Furthermore, both base-intercalated and base-wedged motifs participate directly or indirectly in the formation of RNA functional centers, which interact with various ligands, antibiotics and proteins.