Visualization of distinct substrate-recruitment pathways in the yeast exosome by EM.

The eukaryotic exosome is a multisubunit complex typically composed of a catalytically inactive core and the Rrp44 protein, which contains 3'-to-5' exo- and endo-RNase activities. RNA substrates have been shown to be recruited through the core to reach Rrp44's exo-RNase (EXO) site. Using single-particle EM and biochemical analysis, we provide visual evidence that two distinct substrate-recruitment pathways exist. In the through-core route, channeling of the single-stranded substrates from the core to Rrp44 induces a characteristic conformational change in Rrp44. In the alternative direct-access route, this conformational change does not take place, and the RNA substrate is visualized to avoid the core and enter Rrp44's EXO site directly. Our results provide mechanistic explanations for several RNA processing scenarios by the eukaryotic exosome and indicate substrate-specific modes of degradation by this complex.

Introduction

The eukaryotic exosome is a multi-subunit protein complex crucial for RNA maturation, surveillance, and turnover[1-7]. The exosome core is composed of six RNase PH like subunits (Rrp41, Rrp45, Rrp42, Rrp43, Mtr3, and Rrp46), and three capping subunits containing RNA-binding domains (Rrp4, Rrp40, and Csl4)[8], in a similar architecture to the archaeal exosome[9-11]. However, the exonuclease activity of the core has been lost in yeast and humans[8,12] but gains its RNase activity by binding a tenth subunit, Rrp44 [12,13].

Rrp44 contains multiple functional domains (). Its PilT N-terminus (PIN) domain exhibits manganese-dependent endonuclease (ENDO) activity[14-16]. The C-terminal EXO region consists of tandem cold shock domains (CSD1 and CSD2), a magnesium-dependent 3’- to 5’-exoribonuclease RNB domain, and an S1 domain. The overall architecture of Rrp44's EXO region is similar to that of bacterial RNase II, but the recruitment of RNA to the two enzymes differs[17-20]. Also, while RNase II can only process ssRNA, Rrp44 can unwind and degrade duplex RNA[17], likely via elastic based helicase-like activity[21].

Rrp44 associates with the core complex by binding of the PIN domain to subunits Rrp41 and Rrp45 [18,22,23]. In tune with a similar function of the archaeal exosome, RNA substrates with long single-stranded (ss) 3’-overhangs are first channeled through the eukaryotic exosome core before being degraded by Rrp44 [22,24]. RNA degradation and protection assays indicate that ~31-33 stretched nucleotides are required to reach Rrp44's EXO site from the top of exosome core's RNA binding subunits[22]. This through-core route is clearly demonstrated by the crystal structure of the yeast Rrp44-exosome (RE) in complex with an RNA substrate with a 5’-hairpin and long 3’-ss-overhang[18]. However, degradation of hypomethylated yeast initiator methionine tRNA (tRNA Meti)[6] and tRNAs with a double CCA motif at the 3’-end[25] by Rrp44 implies the presence of an RNA processing pathway that does not necessarily involve channeling through the core. Recent transcriptome data also suggests the presence of alternative routes bypassing the core for RNA substrates with shorter 3’-ss-overhangs to be processed by the exosome[26]. The multi-porous structure of the apo-RE complex provides the potential for RNA substrates to take multiple routes including the through-core and direct access pathways to reach Rrp44's EXO site[22].

To further dissect the mechanism of RNA recruitment to the exosome, we performed biochemical and single particle electron microscopy (EM) analysis on the Saccharomyces cerevisiae RE in concert with different RNA substrates. Our single particle analysis revealed a substrate induced conformational change of the complex upon RNA binding and substrate-specific alternative routes of RNA recruitment by the exosome complex.

Results

RNase assays for both through-core and direct-access routes

The through-core route and direct-access route predict distinct outcomes in the processing of the 3’-ss-overhang of a highly structured RNA substrate. While the through-core route stalls when the long ss-overhang is trimmed down to ~30 nt, the direct-access route bypasses this restriction and predicts a processed ss-overhang of likely less than 10 nt. We used a set of molecular ruler experiments to distinguish these two pathways, by fusing unstructured AU-rich sequences of various lengths at the 3’-end of the highly structured Hepatitis Delta Virus (HDV) ribozyme ( and ).

With a C75U mutation introduced to prevent the HDV ribozyme cleavage, 5’-end 32P-labeled substrates with a stable tertiary structure[27-29] were properly refolded for processing assays (). We carried out RNA processing assays of these substrates by Rrp44 alone and RE complex under single turnover conditions in low Mg2+ (100 μM), where RE was more active[12]. The low Mg2+ condition, however, inadvertently resulted in the partial destabilization of the 3’-portion of the HDV structure beyond the tightly folded pseudo-knot (after nucleotide G40), increasing the length of the 3’-ss-overhang by 33-nt. All HDV RNAs were efficiently processed by both Rrp44 and RE, producing ~5-nt completely digested products as well as one or two predominant processing intermediates (). The EXO active site mutant (D551N, EX−) but not the ENDO active site mutant (D171N, EN−) resulted in the loss of degradation, similar to the behavior of the double mutant (D551N and D171N, DM) (), suggesting predominant exonuclease activity from Rrp44 and RE in our assays. Mapping of the HDV ribozyme with RNase T1 and OH ladders ( and ) indicated that the single processing intermediate in Rrp44-alone reactions represents a stalling event near residue G40 of HDV RNA, presumably due to the blockage of the direct-access pathway by the tightly folded 5’-pseudoknot. Interestingly, the same processing intermediate was also present in RE reactions, suggesting that a similar, direct-access route is present in the complex. The short processing intermediate in RE reactions is unlikely due to excess free Rrp44s because our purification procedure ensures a 1:1 molar ratio between Rrp44 and the exosome core (). The complex's integrity was further supported by EM (). In the RE reactions, an additional processing intermediate ~26 nt longer than the Rrp44 reaction intermediate was evident for all HDV RNAs, consistent with a processing event by the through-core route (). The above results suggest that while the through-core route is predominant (~70%), a direct-access route is simultaneously present for the HDV RNA degradation.

Two distinct RE conformations with different RNA substrates

We further performed 3D reconstruction of RE incubated with ssRNA substrates with different lengths () under non-digesting conditions (with 5 mM EDTA present to sequester the catalytic metal ions), where all RNA substrates were capable to bind RE without detectable degradation (). Interestingly, the series of 3D reconstructions revealed two distinct conformations of RE, depending on the ssRNA length ( and ). When incubated with RNA substrates shorter than 12 nt, the resulting EM reconstruction of the complexes revealed a conformation (RE-short) very similar to that of apo-RE (). By contrast, when the RNA substrates were longer than 14 nt, RE adopted a distinct conformation (RE-long) (). The major difference occurred within the Rrp44 protein, which is more compact in RE-long than in RE-short conformation, but the structure of the core complex is almost unchanged between the two conformers (, and ). This RNA-induced conformational change is not necessarily restricted to purely ssRNAs as the same conformation of RE-long was also observed in the RE-HDV40 complex ().

We docked atomic models into the two conformations to rationalize the observed conformational changes. For the RE-short 3D reconstruction, the crystal structure of the human core exosome (PDB 2NN6)[8] fits nicely in the map as a rigid body, while the crystal structures of Rrp44's PIN domain and C-terminal EXO region from the Rrp44-Rrp41-Rrp45 subcomplex structure (PDB 2WP8)[22] fit well in the map as separate rigid bodies ( and ). In the docked model, the overall spatial relationship among Rrp44, Rrp41, and Rrp45 is very similar as in the crystal structure of the sub-complex except that the PIN domain moves ~10 Å away from Rrp44's EXO region. For the RE-long 3D reconstruction, we took the recently published crystal structure of the yeast RE in complex with a partial hairpin RNA substrate (PDB 4IFD)[18] and docked it in the map as a rigid body. The crystal structure fit with good agreement into the 3D envelope of the RE-long structure, except that the PIN domain of Rrp44 needed to be shifted by about 10 Å closer to the EXO region ( and ). The high level of agreement between the two independently determined structures indicates the strong biological relevance of these two conformations. Compared to the crystal structure, the EM reconstructions all had a substantially weaker density for Csl4 on top of the core, in consistence with previous observations that Csl4 is prone to dissociation from the complex during purification[23,26].

Close examination of the two models revealed that the major structural change occurring in Rrp44 concerns a rigid body rotation (~120 degrees) of its EXO region from the RE-short to RE-long state (RMSD of 4.3 Å between the two atomic models using the ProFit program, Martin, A.C.R., http://www.bioinf.org.uk/software/profit/) (). In the RE-short conformation, there are two solvent exposed openings in Rrp44: the opening between the PIN and the CSD domains, and a major cleft between the CSD1 and the RNB domain that makes the EXO site exposed to solvent (). This cleft has been shown to bind ssRNA substrates that enter the EXO site for degradation[17,18]. The two openings are separated by the CSD domains in the center, with the CSD2 domain blocking in way from the core channel's bottom exit to Rrp44's EXO site ( and ). In contrast, in the RE-long conformation, the Rrp44's region undergoes a dramatic rotation as it rolls against the bottom of the core and pivots around the CSD1 domain, so that the CSD2 domain moves away from the core channel while the RNB domain moves closer to the core channel ( and ). As a result, the cleft between the CSD1 and the RNB domain moves just to the core channel's bottom exit, thus forming a continuous route for ssRNA substrates to thread from the core to Rrp44's EXO site. In agreement with this notion, the RNA molecule from the crystal structure fits nicely within the central channel ().

Our series of 3D reconstructions revealed that only ssRNAs longer than 14 nt induce the dramatic conformational change in Rrp44 within RE, while ssRNAs shorter than 12 nt maintain the same apo-RE conformation. Such a dramatic structural transition in response to the ssRNA's length agrees with the crystal structure of the RNA-RE complex, in which interpretable electron density for the RNA substrate's very 3’-end extends from the bottom exit of the core to Rrp44's EXO site[18]. In the crystal structure, nucleotides numbered −13 ~ −15 of the RNA interact with the Rrp45 and Rrp41 subunits, while −1 ~ −12 interact entirely with Rrp44. This indicates that ssRNA substrates with enough length to connect the core to Rrp44 are ready to induce Rrp44's dramatic rearrangement.

EM of RE in complex with RNAs with long 3’-ss-overhangs

In order to visualize the RNA recruitment pathways to the exosome, we pre-formed a complex between streptavidin (SA) and 5’-biotinylated ss-RNA oligos of various lengths (SA-RNA, and ), incubated the SA-RNAs with RE under non-digesting condition, and sought to locate the SA density in the resulting negative stain EM reconstructions as an indicator of the ss-RNA substrates’ 5’-end. Single particle 2D classification of such samples allowed us to detect additional globular densities attached to the RE (), the size of which agrees well with a SA tetramer. There were substantially more class averages showing the additional SA densities for the ssRNA substrates longer than 36 nt, but very few with SA densities for the ssRNA substrates shorter than 24 nt (). This result was in agreement with pull-down assays using RE as bait to detect horseradish peroxidase (HRP)-conjugated SA-RNA binding (). The observation reflects the fact that SA prevents RNAs with short 3’-ss-overhangs (<24 nt) from forming stable complex with RE but does not affect those with long enough 3’-ss-overhangs (>36 nt) to thread through the core's entry site to Rrp44's EXO site to form a stable complex. We performed 2D classification and 3D reconstruction of the SA-RNA47-RE and SA-RNA50-RE and revealed clear shape of RE with a globular density attached to the entry side of core complex in close proximity to the Rrp4 protein ( and ). The dimension of the globular density agrees very well with the size of a SA tetramer so that we can unambiguously dock a SA atomic model within it ( and ). Interestingly, the 3D model shows the same conformational state of Rrp44 seen in the RE-long conformation ( and ), as verified by comparing the 2D class averages with reprojections from the 3D models of RE-short and RE-long ( and ). This further proves that RNAs going through the core channel and reaching Rrp44 are responsible for RE's conformational change. The fact that RNAs with very different sequences as the 50-mer and 47-mer behave similarly in their interactions with RE underscores the importance of the single stranded nature of RNA substrates over their specific sequences.

EM of RE in complex with RNAs with short 3’-ss-overhangs

As we have revealed, in contrast to the longer RNA substrates, SA labeled RNA substrates shorter than 24 nt do not seem to form a very stable complex with RE in the non-digestive condition in the presence of 5 mM EDTA (), because SA prevents the relatively short ss-overhang to reach Rrp44's EXO site through the core. However, in the presence of 2 mM Mg2+, when the RE is active for RNA processing, we found that SA-RNA substrates with a 3’-ss-overhang shorter than 24 nt can also form reasonably stable complexes with RE that can be characterized by EM (). We were able to detect additional densities corresponding to the dimension of a SA in a portion of class averages from the SA-RNA24-RE sample, albeit fewer than those in the SA-RNA47-RE. The RE portion of 2D class averages of SA-RNA24-RE compared reasonably well with reprojections of the 3D model of apo-RE, allowing us to assign each class average to a specific Euler angle in 3D space. The comparison between the 2D class averages and corresponding reprojections showed almost no SA additional density appearing at the top of the core complex, in contrast to what was seen for the SA-RNA47-RE sample, but more occasional SA densities attached to other parts of RE. Unfortunately, our efforts to get a 3D reconstruction from these 2D class averages were not successful, presumably due to two reasons: 1) The 2D class averages with clear SA densities were not enough to cover the full set of views necessary for 3D reconstruction; 2) There may exist more than one SA labeling sites on RE. Based on the Euler angle assigned to each 2D class average, we back-projected the 2D class averages into 3D space with their SA density as a rod (). This allowed us to capture two major cross points of SA back-projection rods around the RE 3D model, thus allowing us to locate the most probable SA locations attached to RE in 3D space ().

Compared with the SA-RNA47-RE sample, the SA density appeared mostly close to Rrp44 in the SA-RNA24-RE sample. The locations identified are around the RNB domain, facing its front and back sides, respectively. Both locations are close to the cleft between Rrp44's CSD and RNB domains which has been shown by crystallography to be able to recruit short ssRNAs[17,18,22] (). Therefore, the SA densities directly attached to these two sites provide strong evidence of RNA substrates’ direct access to Rrp44 bypassing the core. In correlation with the RNase assays of the HDV ribozyme substrates by Rrp44 and RE, this observation suggests that Rrp44 maintains the capability to directly recruit RNAs when in complex with the core. We also detected a minor population of 2D class averages with the SA density near Rrp44's PIN domain, implying that PIN domain's endonuclease site may also recruit RNA substrates directly.

Time-resolved EM of RNA degradation by RE

To directly monitor exosome's RNA processing reaction, we mixed RE with SA-labeled RNA47 under the processing condition (2 mM Mg2+ in buffer) and took samples at different reaction time points for single particle EM analysis as done above. This revealed 10~30% of 2D class averages with clear additional densities corresponding to RE-bound SAs. Based on the locations of SA in the 2D class averages, we grouped the actively processing SA-RNA-RE complex into three major families (): 1) Particles with a single SA density at the top of RE, corresponding to the complex recruiting RNA substrates via the through-core pathway (); 2) Particles with one single SA density around the bottom part of RE, corresponding to the complex recruiting RNA substrates by direct-access pathway (); and 3) Particles with two SA densities attached to RE at the top and bottom, respectively, corresponding to the complex recruiting RNA substrates via both pathways simultaneously ().

We examined the distribution of particles within each family during the reaction (). Because RNA47 has a long ss-overhang, it is expected to be mainly recruited to the EXO site of Rrp44 via the through-core route. Indeed, in all the reaction time points, family 1 was the most populated among the three families. In agreement with the processive degradation, the population of family 1 particles decreased as the reaction time went on, indicative of the release of shortened SA-RNA substrates from the core channel. In contrast, family 2, albeit less populated than family 1, existed and kept at roughly the same percentage at all reaction time points. Together with a similar assay of SA-RNA24 substrates (), this suggests that the direct-access route of RNA substrate recruitment is less competitive than the through-core route on RNAs with long 3’-ss-overhangs but works as a distinct route for RNA substrates with short 3’-ss-overhangs towards deeper processing. The very minor population of family 3 particles at the beginning of the reaction may correspond to some RE complexes competing the through-core and direct-access routes for RNA recruitment.

The results above suggest that RE may adopt a mixture of states during the processing of RNA substrates. Agreeing with this, 3D reconstruction from a full set of particle images taken of RE in the presence of both RNA and Mg2+ showed a much weaker density of the Rrp44 region (), indicative of a highly mobile nature of the Rrp44 during RNA processing. Supervised classification of the above dataset revealed the coexistence of roughly 30% RE-short and 70% of RE-long conformers (), verifying the coexistence and possible dynamic transition between the two recruitment pathways.

Cryo-EM reconstruction of RE in complex with tRNA substrates

tRNAs are major substrates under the quality control by the exosome and related complexes such as the TRAMP complex, which adds short poly A tail (6 to 15 nt) to the 3’-end of RNAs before processing[30-32]. To examine RE's transactions with tRNA substrates, we produced various tRNAs with different 3’-end polyA tail from 0 nt to 20 nt using in vitro transcription (lacking post-transcriptional modifications), and performed degradation assays with Rrp44 and RE. tRNAs with tails of at least 10 nt were efficiently processed by Rrp44, while the tail-less tRNA only produced faint degradation products (). Note that these results differ slightly from another published report[6] and this is likely because our assays used a four-fold higher concentration of Mg2+ to prevent unfolding of the tRNA. The removal of the 3’-ss-overhang from the tRNA by RE occurs only when the overhang is 20 nt or longer (). Because an intermediate overhang length of 20 nt is too short for channeling through the core but longer than required for degradation by Rrp44 alone, it suggests that direct-access products can be generated by RE, and that the core slightly modifies the pathway of the RNA to Rrp44.

In order to directly visualize tRNA's interaction with RE, we performed cryo-EM 3D analysis of a frozen-hydrated specimen of RE in complex with unmodified yeast tRNAMet. Using single particle 3D maximum-likelihood classification method[33], we classified the dataset and reconstructed two major conformers, one as a model of apo-RE and the other as a model of RE with clear additional density attached to the pocket formed by Rrp44, Rrp45 and Rrp43 at the predicted direct-access route (). The additional density of the latter model can be fitted with the atomic model of a tRNA (PDB ID 4TNA)[34] with its 3’-end facing the tunnel to the EXO site of Rrp44 (). Compared to the apo-RE, the tRNA-bound RE has its C-terminal EXO domain shifting down with about 20 Angstroms likely due to the tRNA insertion (). These results further support the direct-access route for biologically relevant structured RNA substrates within RE.

Discussion

Compared to its bacterial and archaeal homologues, the eukaryotic exosome core complex has lost the phosphorolytic activity but evolves a function to channel the 3’-end of long ssRNA substrates into the hydrolytic exonuclease Rrp44 for processive degradation. Such a through-core route has been supported by biochemical data, EM analysis, and, most recently, X-ray crystallography[18,22,24]. To date, a few crystal structures of the RE complex, its sub-complexes, or individual proteins are available but don't fully agree: Rrp44's EXO region adopts very different orientations in the 5’-hairpin-RNA-bound RE complex[26] and the Rrp44-Rrp41-Rrp45 sub-complex[22]. Our work shows that, depending on the RNA substrates, RE adopts at least two conformations, in which Rrp44's orientations agree to the two crystal structures, respectively. Furthermore, the 3D reconstructions of RE bound with a series of ssRNA oligos illustrate that Rrp44's conformational change observed in the crystal structure of the holo-RE complex with RNA substrate is solely induced by the 3’-ss portion of the RNA substrate rather than the presence of the 5’-hairpin or crystallographic packing.

Another agreement between our single particle EM reconstruction analysis and the crystal structure of RE with RNA substrates is the identification of the minimum length of 3’-ss portion as 14 nt to induce the conformational change of Rrp44. RNAs longer than 13 nt would form a continuous linker bridging Rrp45 and Rrp40 of the core and Rrp44. Thus, the continuity of the long ssRNA chain and its multiple interaction sites with RE is essential in holding Rrp44's EXO region in a new orientation so that the cleft between the CSD1 and RNB domains aligns to the exit of core channel (). The fact that such a conformation exists only in the presence of long 3’-ss-overhangs indicates that this conformation is unstable upon the loss of RNA bridging between the core and Rrp44. Rrp44 has the tendency to resume its apo-state, when the ssRNA becomes too short to hold Rrp44's EXO region with the core. We therefore hypothesize that RE breathes between the apo-RE and RE-long conformations so that it can generate a Brownian Ratchet-like pulling force on the 3’-ss-overhang inside the core towards Rrp44's EXO site (). This force, together with the intrinsic elasticity of Rrp44 [21], probably contributes to the processivity and weak RNA duplex unwinding activity of the entire complex. Furthermore, this suggests a mechanism of how RE releases RNA substrates with stable secondary or ternary structure at the 5’-end that cannot be unwound by the complex. When the 3’-end of such RNA substrates is processed to a length that can barely hold Rrp44 in the active processing conformation, the tension caused by resuming the apo-conformation could release the RNA's 3’-end from the protein and eventually release those stable RNA substrates from RE for downstream processing and vacate the complex for recruiting new RNA substrates ().

The physiological role of the through-core route for RNA degradation has been firmly established from both in vitro and in vivo assays. The most recent yeast genetics analysis revealed that the central channel of the exosome core is essential for function and viability of budding yeast[35,36]. But this does not preclude the presence of an alternative RNA recruitment route, which could be very important in processing a subset of RNA substrates. The existence of a possible direct-access pathway for RNA recruitment and processing by RE was first hypothesized based on our EM reconstruction of the RE complex[23], as well as implied from genetic and biochemical work from others[6,25], but no direct observation of such a route was available. Our data here provides the first direct visual evidence of such an RNA recruitment pathway that bypasses the core channel. Compared to the through-core route, the direct-access route has a weaker affinity for RNA, yet contributes significantly to the RNA processing even at the beginning of a reaction. The direct-access route likely plays a more dominant role in the processing or degradation of RNA substrates with 3’-ss-overhangs shorter than 30 nt, for both single stranded and structured RNAs. This agrees with the transcriptome analysis that reported structured RNAs with short 3’-ss tails contacting directly with Rrp6 and Rrp44 while bypassing the core complex[26]. The same work also revealed significant substrate spectrum overlap between Rrp6 and Rrp44 with each enzyme being responsible for the processing of different regions of the same substrate. It is plausible that evolutionary pressure preserves the two conformations of RE observed by us for different purposes of RNA processing. While the through-core pathway may be mainly for RNA degradation, the direct-access pathway might be for RNA processing, which removes the 3’-ss-overhangs from structured RNAs. It is also possible that a structured RNA with a long 3’-ss-overhang may first be processed via the through-core pathway, then further cleaved by the direct-access pathway after retracting from the core. In vivo, activities of the two recruitment pathways may be regulated by multiple factors such as the RNA substrate themselves, as for tRNA in this work, and cofactors such as the Rrp6, TRAMP, and SKI complexes. Indeed, Rrp6 has been shown to deviate certain RNA substrates from the core directly to Rrp44 [36]. Future research needs to be conducted to elucidate the mode of RNA recruitment to the exosome in the context of multiple cofactors to more holistically define the role of the exosome in RNA metabolism in the eukaryotic cell.

Online Methods

Tandem affinity purification (TAP) and reconstitution of S. cerevisiae exosomes

Purification of the core and Rrp44-exosome was done as described with minor modification[23]. Briefly, ~500 gram of TAP-Rrp46ΔRrp6 yeast cells were lysed with bead-beaters, and after ultracentrifugation the supernatant was loaded onto 3 ml of Rabbit IgG resin (Sigma-Aldrich). The resin was extensively washed with high salt buffer (500 M NaCl), and the Rrp44-exosome was eluted by on-column TEV cleavage to remove the protein A tag. For purifying the core exosome, the Rrp44 was removed from the complex by washing with 40 ml of 1.2 M MgCl2. The mutant Rrp44-exosome complexes were formed by incubating a 50% excess amount of purified recombinant mutant Rrp44 (EX−: D551N, EN−: D117N, or DM: D551N D117N) with the core exosome at 4 °C overnight. All Rrp44-exosome complexes used in the biochemical assays, wild-type or mutant, were further purified on an ion-exchange Mono Q 5/50 GL column (GE Healthcare), where any excess Rrp44 and endogenous nucleic acids were well-separated from the Rrp44-exosome complex (). Samples were flash frozen in liquid nitrogen and stored at −80 °C until later use.

RNA degradation assays

HDV and tRNA substrates were produced by RNA polymerase T7 based in vitro transcription. RNA was separated on a 10% polyacrylamide/8 M urea gel and RNA visualized by UV shadowing. The RNA band was excised from the gel, crushed, and eluted in nuclease free deionized water at 4 °C overnight.

The purified RNAs were dephosphorylated with FastAP Alkaline Phosphotase (Fermentas), phenol extracted, ethanol precipitated, and dissolved in nuclease free water. RNA was 5’ end labeled with Polynucleotide Kinase (New England Biolabs) and 2 μl of 3000 Ci/mmol 32γP-ATP (Perkin Elmer), purified on 10% polyacrylamide/8 M gels, eluted and ethanol precipitated. RNA pellets were washed with 80% ethanol, air dried, dissolved in a final volume of 30 μl nuclease free water, and stored at −20°C until use.

Immediately before use, RNA substrates were refolded in 2X reaction buffer (20 mM HEPES pH 7.5, 100 mM NaCl, and twice the concentration of MgCl2 used in the assay) by heating to 80 °C for 5 minutes followed by snap cooling at 4 °C. All reactions were performed in 1X reaction buffer, where the molar ratio of protein to RNA was 2:1, and incubated at 30 °C for indicated time points. The reactions were terminated by adding an equal volume of formamide plus 10 mM EDTA and heating the solution at 65°C for 10 minutes. The reaction products were analyzed on 16% acrylamide/8M Urea sequencing gels in 0.5X TBE and quantified by a Typhoon Phosphorimager (GE Healthcare).

Preparation of exosome coupled IgG resin beads

An amount of 100 ml supernatant of TAP-Rrp46-ΔRrp6 yeast lysate was incubated with rabbit IgG affinity resin beads (Sigma-Aldrich, St. Louis, MO) at 4 °C for 2 hours with constant stirring. The resin was then centrifuged at 500 rpm for 2 minutes and washed three times each with 10 ml washing buffer (50 mM Tris-HCl pH 8.0, 150 mM NaCl, 10% glycerol with 2mM DTT, 2mM EGTA and 1 tablet of complete EDTA-free protease inhibitor mixture from Roche Applied Science). The resin was then washed three times each with 10 ml of high-salt buffer (washing buffer plus 250 mM NaCl) and the beads were re-equilibrated three times with 5 ml washing buffer.

Gel shifting and pull down assays

Streptavidin (SA) was purchased from Sigma and HRP-coupled streptavidin (HRP-SA) was from Abicom. For gel shifting assays, proteins and RNAs were mixed together in a 4:1 molar ratio in nuclease free water and incubated on ice for 5 minutes, then run on a 6% acrylamide native gel in 0.5X TBE buffer. The gels were subsequently stained with ethidium bromide and imaged in a regular UV illuminated gel scanner.

For pull down assays, HRP-SA and RNAs were mixed together in a 4:1 molar ratio and incubated on ice for 5 minutes. 1μl of HRP-SA-RNA complex (2.5uM) was incubated with 30μl of exosome coupled IgG resin in 500 μl of reaction buffer (50 mM Tris·HCl pH 8.0, 150 mM NaCl, 5mM EDTA, 1mM DTT) at room temperature (~22 °C) for 15 minutes with occasional stirring. The mixture was centrifuged at 500 rpm for 2 minutes, and the resin was washed three times each with 1.5 ml reaction buffer. 2 μl supernatant or beads were loaded on nitrocellose film. The film was subsequently incubated in ECA agent for HRP signal development as in normal Western Blotting.

EM sample preparation

The RNA-RE complex for EM analysis was prepared by mixing RNA and RE in a 1.2:1 molar in a non-digestive buffer containing 50 mM Tris-HCl pH 8.0, 150 mM NaCl, 5mM EDTA, and 1mM DTT at room temperature for 15 minutes. For preparing SA-RNA-RE sample in non-digestive condition, SA and biotinylated RNAs were first mixed together in nuclease free water and incubated on ice for 15 minutes in a 2:1 molar ratio and diluted to a final concentration of 6 μM of RNA. Subsequently, the SA-RNA mixture was incubated with RE complex in non-digestive buffer as described above in a molar ratio of 1.2:1. All samples were diluted at a final concentration of ~80 nM of the exosome in the non-digestive buffer and negatively stained in 2% (w/v) uranyl acetate solution following the standard deep stain procedure on holey-carbon coated EM copper grids covered with a thin layer of continuous carbon [38]. For SA-RNA-RE samples in magnesium condition, everything followed the same procedure except for that the non-digestive reaction buffer was replaced by a magnesium buffer (50 mM Tirs-HCl pH 8.0, 150 mM NaCl, and 2 mM MgCl2) and the incubation time of SA-RNA and RE was 5 minutes.

Time-resolved EM sample preparation

To monitor the structure of SA-RNA being processed by RE in a degradation reaction, we first mixed 4 μl of SA-RNAs produced as described above with 2 μl of 3X reaction buffer (150 mM Tirs-HCl pH 8.0, 450 mM NaCl, 6 mM MgCl2) at 30 °C. Subsequently, we added 4 μl of RE (2.5 μM dissolved in 50 mM Tris-HCl pH 8.0, 150 mM NaCl, and 2 mM MgCl2) to the reaction system to initiate the reaction. At given reaction time points, the reaction was terminated with the addition of 10 μl termination buffer (50 mM Tris-HCl pH 8.0, 150 mM NaCl, 20 mM EDTA) and immediate cooled on ice. 4 μl of the final solution was then negatively stained.

Electron microscopy and image processing

We examined all the negative stained specimens on an FEI Tecnai-F20 electron microscope operated at 200 kV acceleration voltage at 50,000 nominal magnification with a range of defocus from 0.8 to 1.2 μm. The electron micrographs were taken on a Gatan Ultrascan4000 4k X 4k CCD camera.

The digital micrographs taken from CCD camera was then processed using EMAN2 package [39] to semi-automatically pick the individual particle images of the complexes and using IMAGIC-4D package[40] to determine and correct the contrast transfer functions using phase-flip method. For each sample, 20,000-50,000 negatively stained particle images were analyzed two-dimensionally using IMAGIC-4D. Briefly, we performed iterative multi-variant statistical analysis (MSA) and multi-reference alignment (MRA) cycles on the particle stacks using IMAGIC-4D to obtain 2D class averages as described previously[38]. We performed 3D reconstruction of all the samples using projection-matching refinement of particle stacks against the same initial model of apo-RE (EMD1439 [23]) low-pass filtered to 70 Å in SPIDER[41] as described previously[38]. The resolution of final reconstruction was estimated using Fourier Shell Correlation (FSC) algorithm at a criteria of 0.5. The number of particles contributing to and the resolution of each 3D reconstruction (measured by FSC 0.5 criteria) are summarized in . The distinct differences among the final reconstructions indicated that our results are free of model bias. This was also verified by using a totally different initial model of an elongated bar to converge to the same structures of RE sample after about 30 rounds of projection-matching refining iterations. The model-bias free nature of the RE sample is very likely due to its strong intrinsic structural features in the high contrast negative stained images.

The docking of atomic models were performed by UCSF-Chimera[42] using fit-model-in-map functions. All the 3D models in the figures are rendered in UCSF-Chimera.

Cryo-EM of RE in complex with tRNA

We incubated 250 nM RE with 1 μM unmodified yeast tRNAMet in reaction buffer (50mM Tris-Hcl, 150mM NaCl, 2mM MgCl2) on ice for 10 min. 4 μl of the reaction solution was then applied to glow-discharged C-flat grids (1.2/1.3) covered with a layer of continuous carbon with a thickness of ~4 nm. The grids were then blotted and plunged into liquid ethane in a FEI Vitrobot Mark IV. The frozen-hydrated grids were subsequently examined under a FEI Titan Krios electron microscope operated at 120 kV acceleration voltage. Micrographs of RE were collected using the AutoEMation software [43] installed on the microscope at low-dose condition with a dose of ~20 electron/Å2 and a defocus value ranging from −1.2 to −4 μm. The micrographs were collected on a FEI Eagle CCD camera with a pixel size of 1.5 Å. We used EMAN2 package to perform raw micrograph screening, CTF determination, and semi-automatic particle picking[39] and obtained a dataset of about 61,000 particle images. We then performed 3D refinement and classification of the particle images using the Relion package[44]. After the first 25 iterations of classification, we identified two major classes of particles (accounting for more than 56,000 images) with good quality of reconstruction out of six classes. We further performed 3D refinement of the two classes and got a final reconstruction of the RE-tRNA complex from about 28,800 particles and of the apo-RE complex from about 27,200 particles, respectively. The number of particles contributing to and the resolution of each 3D reconstruction (measured by FSC 0.5 criteria) are summarized in .

Back-projection analysis of SA-RNA24-RE complex

For back-projection analysis of the SA-RNA24-RE complex, we calculated the Euler angle of the 2D class averages by firstly masking off the SA densities in the 2D class averages and then matching them to reprojections of the 3D model of apo-RE. The corresponding view of the 3D model of apo-RE matching to a 2D class average was then displayed in UCSF-Chimera. A 3D model of rod perpendicular to the view was added in UCSF-Chimera to represent the back-projection of the corresponding SA density in the 2D class average. By applying this procedure to all the views of the 3D model matching to the 2D class averages with SA density, we finally visualize the cross-point of all the rods, assuming that the most crossed points have a higher probability to represent the location of SA in the 3D model.

Table 1

Cross-correlation coefficients among different 3D models

	Apo-RE	RNA08-RE	RNA10-RE	RNA12-RE	SA-RNA47-RE
RNA47-RE Core	0.9323	0.9265	0.9287	0.9135	0.8892[*]
RNA47-RE Rrp44	0.6822	0.6443	0.6573	0.6573	0.6849

The upper row is the comparison among different models at the core; the bottom row is the comparison among different models at the Rrp44 portion.

The comparison is calculated at the core region only, with SA density removed.