RNA silencing refers to a group of RNA-induced gene-silencing mechanisms that developed early in the eukaryotic lineage, probably for defence against pathogens and regulation of gene expression. In plants, protozoa, fungi, and nematodes, but apparently not insects and vertebrates, it involves a cell-encoded RNA-dependent RNA polymerase (cRdRP) that produces double-stranded RNA triggers from aberrant single-stranded RNA. We report the 2.3-Å resolution crystal structure of QDE-1, a cRdRP from Neurospora crassa, and find that it forms a relatively compact dimeric molecule, each subunit of which comprises several domains with, at its core, a catalytic apparatus and protein fold strikingly similar to the catalytic core of the DNA-dependent RNA polymerases responsible for transcription. This evolutionary link between the two enzyme types suggests that aspects of RNA silencing in some organisms may recapitulate transcription/replication pathways functioning in the ancient RNA-based world.
Citation: Salgado PS, Koivunen MRL, Makeyev EV, Bamford DH, Stuart DI, Grimes JM (2006) The Structure of an RNAi Polymerase Links RNA Silencing and Transcription. PLoS Biol 4(12): e434. doi:10.1371/journal.pbio.0040434
Academic Editor: Gregory Petsko, Brandeis University, United States of America
Received: July 20, 2006; Accepted: October 13, 2006; Published: December 5, 2006
Copyright: © 2006 Salgado et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The work was supported by the Human Frontier Science Program, the Academy of Finland (Center of Excellence Program 2006-20011, grants 1213467 and 1213992 to DHB) and the Medical Research Council, UK. MRLK is supported by Helsinki Graduate School in Biotechnology and Molecular Biology, EVM is a Leukemia and Lymphoma Society Fellow. JMG is supported by the Royal Society; DIS and the Oxford Protein Production Facility (OPPF) are supported by the UK MRC and European Commission grant number QLG2-CT-2002–00988 (SPINE).
Competing interests: The authors have declared that no competing interests exist.
Abbreviations: cRdRP, cell-encoded RNA-dependent RNA polymerase; DdRP, DNA-dependent RNA polymerase; DPBB, double-psi β-barrel; dsRNA, double-stranded RNA; MAD, multiwavelength anomalous diffraction; NTP, nucleotide; PTGS, posttranscriptional gene silencing; rmsd, root mean square deviation; RNAi, RNA interference; TGS, transcriptional gene silencing
RNA silencing, or RNA interference (RNAi), refers to a group of RNA-induced gene silencing mechanisms that developed early in the eukaryotic lineage and play essential roles in cellular immunity, modulation of chromatin structure, and development [1–3]. RNAi can induce transcriptional gene silencing (TGS) via chromatin repression or posttranscriptional gene silencing (PTGS) by degradation of target RNAs. RNA silencing pathways use double-stranded RNA (dsRNA) triggers, processed by Dicer, to form short interfering RNAs (siRNAs) of 21–25 nucleotides . One siRNA strand is recruited by an effector complex containing the Argonaute protein and used as a guide for sequence-specific degradation of target mRNAs (in PTGS) or directed silencing of cognate chromatin domains (in TGS) [5,6]. A cell-encoded RNA-dependent RNA polymerase (cRdRP) is also involved, in plants, fungi, protozoa, and certain animals, but apparently not insects and vertebrates, producing dsRNA triggers and hence amplifying the PTGS response [7–9]. Furthermore, cRdRPs may also interact with the cellular transcription apparatus and effect chromatin silencing .
Transgene silencing in Neurospora crassa (quelling) is one of the best studied models for RNAi. Three genes products were identified by screening for quelling defective phenotype, QDE-1 (cRdRP), QDE-2 (Argonaute), and QDE-3 (RecQ-like helicase) . Subsequent mechanism studies demonstrated the additional involvement of two Dicer-like proteins: DCL-1 and DCL-2 . We have previously isolated recombinant QDE-1 and an enzymatically active C-terminal portion of the molecule (QDE-1 ΔN, residues 376–1,402), and have shown that QDE-1 efficiently produces full-length copies or short 9–21-nt copies scattered throughout the input ssRNA templates, as well as extending, rather inefficiently, complementary primers .
In order to understand the functionality of the QDE-1 in more detail, we have determined the high-resolution crystal structure. This reveals the molecule to be dimeric and containing at its core two subdomains responsible for catalysis. These subdomains each have the topology of double-psi β-barrels (DPBBs), and are similar to (and disposed in a similar fashion to) two separate subunits in the DNA-dependent RNA polymerases (DdRPs) that perform transcription across the domains of life. The structure not only suggests how the molecule might efficiently produce dsRNA triggers, but may also add a piece to the jigsaw that relates the RNA world to the DNA world, and provide a model for all cellular RdRPs.
QDE-1 ΔN Is a Dimer with an Active Site Formed from DPBBs
The structure of QDE-1 ΔN was initially solved by multiwavelength anomalous diffraction (MAD) analysis of a crystal of selenomethionated protein, expressed in yeast, to a resolution of 3.2 Å . The structure was then refined against a higher resolution dataset extending to 2.3 Å. These crystals belong to space group P21, and the crystallographic asymmetric unit contains two subunits, which together form a compact pyramidal object with dimensions of 90 × 70 × 127 Å3. The two constituent subunits are related by an approximate 2-fold axis and are in such intimate contact (over 2,000 Å2 of contact area per subunit) that we would expect this to represent a functional dimer . This was confirmed by gel filtration and sedimentation assays, with a molecular weight (MW) approximately equal to 230 kDa (compared to a predicted ∼120 kDa/monomer). QDE-1 ΔN has 1,026 residues, and the final model contains 933 residues in subunit A and 930 residues in subunit B (residues are missing from the N- and C-termini and from some flexible loops, see Figure 1). The refined structure at 2.3 Å has an Rwork of 21.7% (Rfree = 26.4%), and the stereochemistry is good (root mean square deviation [rmsd] bond = 0.013 Å, rmsd angle = 1.9°, 2.3% of residues are in disallowed regions of the Ramachandran plot; Table 1). In addition, a lower resolution structure (3.5-Å resolution) was determined in space group C2 (see Table 1). In space group C2, the dimer is formed by two subunits related by exact crystallographic 2-fold symmetry. Both crystal forms use the same residues to stabilize the dimer, and these residues are primarily contributed by the head domains (head domain contacts account for 1,720/1,710 Å2 of the contact area of 2,215/2,100 Å2 per subunit in the P21/C2 space groups, respectively ).
(A) Cartoon representation of the QDE-1 ΔN P21 dimer–subunit A coloured according to domains: slab, blue; catalytic, purple; neck, pink; and head, orange; subunit B coloured grey. The approximate non-crystallographic 2-fold is represented as a grey line. Green spheres mark Mg2+ ions. Disordered regions correspond to ten residues at the N-terminus, 30 residues at the C-terminal, 45 residues in monomer A, and 48 in monomer B belonging to four loops (A: residues 590–603, 628–640, 1,241–1,251, and 1,271–1,281; B: 591–606, 627–640, 1,241–1,251, and 1,271–1,281).
(B) Topology of QDE-1 ΔN subunit A, coloured as in (A). The catalytic subdomains DPBB1, DPBB2, and flap are denoted by boxes. The non-crystallographic 2-fold is represented as in (A). Disordered loops are represented by dashed lines.
(C) View of presumed active site. The two DPBBs that form the active cleft, DPBB1 (residues 680–782) and DPBB2 (residues 916–1,018), are labelled. Sequence motifs conserved across other cRdRPs are highlighted: motif 1, red; 2, orange; 3, dark yellow, 4, purple; 5, dark pink; 6, bright pink; and 7, blue (see Figure S1). Invariant residues in cRdRPs are shown in ball-and-stick representation (O, red; and N, blue); Mg2+ ion shown as green sphere.
(D) Zoom of active site region, with conserved residues labelled. Representation is as for (C).
Each QDE-1 subunit contains 41 α-helices and 25 β-strands, creating a four domain fold previously undescribed for RdRPs (Figure 1A and 1B). Distal to the molecular twofold axis is a mixed α/β “slab” domain composed of the approximately 250 N-terminal residues (390–646). The polypeptide chain then leads into a “catalytic” domain (residues 647–807; 914–1,161) that houses the three proposed catalytic aspartic acid residues  (D1007, D1009, and D1011; Figure 1C and 1D) within an active site cleft formed between two DPBB  subdomains (DPBB1, residues 690–792, and DPBB2, residues 916–1,018). The catalytic domain also contains a separate, mainly α–helical, “flap” subdomain (residues 1,025–1,161), peripheral to the active site cleft. The “neck” domain comprises three long α-helices (residues 808–836, 817–913, and 1,162–1,195), which lie close to the molecular twofold axis, connecting the catalytic domain to the mainly α-helical “head” domain (residues 837–888 and 1,196-1372; Figure 1A and 1B).
Multiple alignment of cRdRP amino acid sequences reveals the presence of seven motifs containing invariant residues: motifs 1–3 map to DPBB1, 4–6 to DPBB2, and 7 corresponds to α-helices 29 and 30 at the inner face of DPBB2 (Figures 1B, 1D, and S1). The conserved motifs therefore cohere in three dimensions, with at their heart the proposed catalytic aspartates that reside in a loop (residues 1,007–1,011, motif 6) in DPBB2 at the interface of the two β-barrels. Analysis of electron density maps showed a Mg2+ ion coordinated by these aspartic acid side chains (the functional significance of this is underscored by the observation that the QDE-1 polymerase activity is dependent on divalent cations ). DPBB1 contributes several positively charged residues to the active cleft, which include three invariant residues: Q736, K743 (motif 2), and K767 (motif 3). Together these establish a network of hydrogen bonds with water molecules, linking the two DPBB subdomains (Figure 1D). Although we have been unable to determine structures for complexes of QDE-1 with RNA and/or nucleotides (NTPs), we are able to infer functional aspects from the architecture of the uncomplexed molecule and its unexpected similarity to other polymerases, as discussed below.
The QDE-1 ΔN molecule has several distinct channels and cavities: first, there is a channel formed between the slab and head of each subunit, which is highly positively charged and leads to the active site (Figure 2A); we propose that these channels accommodate dsRNA product. Second, there is a small, negatively charged tunnel at the bottom of each subunit (formed between the flap and DPBB subdomains), which communicates with the active site and may be a route of entry for NTPs (Figure 2A). Third, there is a single tunnel, for which we cannot propose a function, formed between the neck domains and the catalytic domains, bridging the two active sites. The proposed dsRNA product binding channels are not identical in the two subunits, because the disposition of the domains is different: subunit A has a “closed” conformation, with the head and slab clamped down on the active site cleft, whereas an 11° rotation of the head and 2° rotation of the slab render the B subunit more open and provide space for an RNA duplex (Figure 2B and 2C). In the lower resolution C2 crystal form, the molecule forms a symmetric dimer. Here, the subunits both assume a partially closed conformation (the head is rotated by 4° relative to A and 8° to B, and the slab is displaced upwards and outwards by 4° relative to A and 2° to B). Overall, QDE-1 exists as a compact, but flexible, dimeric enzyme with metal binding sites confirming positions of the catalytic sites but with extensive additional structure whose biological functions are less immediately clear.
(A) Surface charge representation of subunit A of QDE-1 ΔN (subunit B coloured grey). Blue arrow indicates proposed RNA product groove; red arrow indicates proposed NTP tunnel; and purple arrow indicates tunnel linking active sites in the dimer. The Mg2+ ions are shown as green spheres. Left panel: dimer (view similar to Figure 1A). Right panel: sliced view of left panel to reveal the tunnels.
(B) Different conformations of QDE-1 ΔN subunits modify polymerase active site accessibility. Left panel: cartoon representation of superposition of observed subunit conformations, coloured according to domain definition: subunit A as in Figure 1A; subunit B: slab, cyan; catalytic, magenta; neck, wheat; and head, yellow; C2 crystal form subunit (semi-transparent): slab, marine; catalytic, violet; neck, light salmon; and head, orange. In all panels, Mg2+ ions are represented as green spheres. The direction of the movement of the slab and head domains is indicated by grey arrows, and the non-crystallographic 2-fold is shown. Central and right panels: RNA-DNA duplex (from an elongation complex of yeast RNApolII) fitted (using the operators for superposition of yeast RNApolII onto QDE-1) into QDE-1 ΔN closed and open subunits (represented as molecular surfaces). Centre panel colours subunit A (closed conformation). Right panel shows a view rotated by 180° with subunit B (open conformation) coloured. RNA-DNA duplex model is show in ball-and-stick representation.
(C) Close-up of the duplex model from (B). Domains are coloured as for the central panel of (B). The slab domain of molecule A would lie across the front of the figure and has therefore been removed for clarity.
QDE-1 Active Site Is Closely Similar to Those of DdRPs
A search using the QDE-1 ΔN active site residues (ASSAM ) identified similarities with the yeast  and bacterial  DdRPs. Superposition revealed that the DPBB subdomains in QDE-1 and the DdRPs are structurally very similar and almost identically disposed (∼10° change in the relative positions of the subdomains; Figure 3). The DdRPs have a DPBB in each of two largest subunits (β′ and β subunits in the bacterial enzyme, and Rbp1 and Rbp2 subunits in yeast RNA pol II), with the first contributing the catalytic aspartates to the active site and the other a set of positively charged residues . In QDE-1 ΔN, the two subdomains have a similar segregation of chemical roles, but they are arranged sequentially on a single polypeptide chain. The similarity of the DPBB bearing the catalytic aspartates between the cRdRPs and DdRPs had been predicted from sequence analysis ; however, there is no sequence homology detectable in the second DPBB or elsewhere in the molecule (see Figure 4). Superposition matches 81 residues of DPBB2 in QDE-1 ΔN with the bacterial β′ DPBB and 85 residues with the yeast Rbp1 DPBB (2.2 and 2.1 Å rmsd in Cαs, respectively). DPBB1 is somewhat less similar to the homologous domain in the DdRPs (74 and 67 residues matched with rmsd of 3.0 and 3.1 Å for bacterial and yeast, respectively; Figure 3). The QDE-1 ΔN catalytic aspartates lie within 1.4 Å of the bacterial and yeast catalytic residues, with equivalent Mg2+ coordination to metal A in the yeast enzyme . Moreover, the bridge helices in DdRPs, proposed to be important for nucleic acid–protein interactions during translocation of the duplex [18,21], are structurally equivalent to helices 27 and 28 in QDE-1 ΔN, suggesting a similar role. Superposition onto QDE-1 of the structure of yeast DdRP complexed with an RNA–DNA duplex  maps the duplex into the putative RNA product groove in QDE-1 (Figure 2B) and the proposed NTP tunnel matches well with that proposed in yeast RNA pol II . Intriguingly, only about ten base pairs of duplex RNA can be modelled into this groove without severe steric clashes with the head domain, suggesting a steric basis for modulating the length of RNA synthesised (Figure 2C).
Stereo representation of the superposition of QDE-1 ΔN and yeast DdRP DPBBs (view as in Figure 1C). Structurally equivalent residues in QDE-1 ΔN are coloured dark purple (non-equivalent residues in light purple). Equivalent residues in yeast DPBBs are coloured green (non-equivalent residues in semi-transparent grey). QDE-1 ΔN and yeast (D481, D483, and D485) active site aspartates are coloured yellow and green respectively.
Alignment of QDE1 ΔN (top sequence), bacterial (middle sequence, orange), and yeast (bottom sequence, green) polymerases is based on structurally equivalent residues, as determined by SHP . Residues structurally equivalent in all polymerases are shaded purple, green (light for Rbp1, dark for Rp2) if equivalent in yeast DdRP and QDE1 ΔN, and orange if only equivalent in QDE1 ΔN and bacterial DdRP. Invariant residues are shaded in red. Conserved sequence motifs identified in cRdRPs are represented as in Figure S1, marked on QDE1 sequence. QDE1 secondary structure elements are shown on top, coloured according to domain definition (slab, blue; catalytic, deep purple; neck, pink; and head, red). DPBB1 and DPBB2 are outlined by deep purple boxes. The flap sub-domain and the potential “bridge helix” are also represented by boxes, coloured light purple and grey, respectively.
Overall, the DdRP structures are much more elaborate than QDE-1 ΔN, and we can detect no significant structural similarity beyond the vicinity of the active site. Nevertheless, there may be some functional relationship at the level of protein domains. The flexible head that forms part of the proposed QDE-1 product groove might be equivalent to the clamp in yeast Rbp1 : the head domain closing down on the slab to stabilize the RNA product during polymerisation (Figure 2B). The yeast Rbp2 protrusion-lobe is also displaced during transcription  to accommodate and stabilize the RNA–DNA product, and the somewhat flexible QDE-1 slab may play a similar role. These observations suggest a similar mode of action for the DdRP and cRdRPs, a view reinforced by the proposal that a human DdRP, Pol II, is involved in replication of hepatitis Δ RNA .
QDE-1 May Act as a Two-Stroke Motor for the Production of dsRNA Triggers
QDE-1 ΔN is active in solution , and we have shown here that this form is dimeric. What, if anything, might the functional significance of the dimer be? Three roles for cRdRPs have been discussed: (1) conversion of aberrant RNAs into dsRNA triggers for PTGS, (2) primer-dependent amplification of RNAi triggers, and (3) chromatin silencing [2,8–10]. Reaction (1) presumably involves initiation at the 3′ end of aberrant mRNA and processive RNA synthesis, whereas reactions (2) and (3) would require internal recruitment of the polymerase to its RNA targets. Both initiation modes function in vitro [9,24]. The more open conformation of subunit B would allow internal initiation, but conversely, the partial closure of the RNA product groove in the dimer might favour initiation at the 3′ end of an RNA template. It seems likely that, in the dimer, only one catalytic site is active at any given time, represented by the closed conformation of subunit A, with the inactive subunit being held open. This suggests a mechanism whereby binding to one active site primes the other. Thus, the molecular architecture might favour the molecule working as a “two-stroke motor,” with facilitated active site switching as the subunits cycle back and forth between conformations in response to RNA binding. This model has two attractive features: first, by tethering the RNA template to the dimer, re-initiation will be efficient, and second, initiation at active site B can be coupled to the activity at active site A, by molecular switching, driven by the steric clashes with the rather stiff dsRNA product suggested by our modelling studies. Overall, this mechanism might lead to the effective production of appropriate-length dsRNA triggers.
A Polymerase for the RNA World
From the structure of QDE-1, it appears that all cellular RdRPs as well as DdRPs are related but structurally distinct from viral RdRPs, which have the “right-hand” architecture [25–30]. The similarity revealed here between the active sites of cRdRPs and the DdRPs indicates that they share a common ancestor. We refer to this family of enzymes as the “double-barrel” polymerases. The RNAi polymerases present the entire active site on a single polypeptide chain, suggesting that they are closer to the ancestral protein than the present day transcription polymerases, which are composed of up to 12 separate polypeptide chains, and in general bear the two active site β-barrels on separate chains. Interestingly the DdRP of Helicobacter pylori also has its active site on a single chain, due to the fusion of the rpoB and rpoC genes, coding for the β and β′ subunits . Indeed, it has been shown in vitro that a fused Escherichia coli β-β′ protein can assemble into a functional polymerase, although there is presumably some advantage in vivo for separate β and β′ chains . Classifying the evolutionary relatedness of polymerases on the basis of their structural similarity, reveals that the “right-handed” and “double-barrel” polymerases are tightly grouped and quite separate from each other and also from Polβ (Figure 5A). The implied path for the evolution of the double-barrelled RNA polymerases is shown in Figure 5B. The original enzyme presumably possessed a single DPBB, bearing the essential catalytic apparatus. This molecule may have taken over the role of RNA polymerisation, presumably from an RNA molecule, very early in evolution. Gene duplication led to a molecule with two DPBBs on a single polypeptide chain, which differentiated to produce a QDE-1–like molecule capable of efficient RNA polymerization, a key development in the elaboration of the RNA-based world. During this process, DPBB2, containing the active site aspartate residues, was conserved relatively strongly, whereas DPBB1, which acts as an accessory domain, essentially lost any detectable sequence conservation (Figure 4). The ability of such molecules to act on DNA could then underpin the switch to DNA as the repository of genomic information. Splitting of the DPBBs onto separate subunits would facilitate the radical evolution of this first DdRP into the complex, highly regulated, transcription machines with particular functions delegated to specialist subunits, which we observe today. In summary, it seems that this polymerase part of the RNA silencing machinery may give us a glimpse far back in time, providing insight into the evolution of a protein-based mechanism for the transmission of RNA genomic information in an RNA-based world.
(A) Evolutionary phylogenetic tree for known polymerases based on structural similarity (for description of the method see ). The “right-handed” polymerases, “double-barrel” polymerases, and polβ appear to form three separate, unrelated families. Within the “double-barrel” family, the DPBB2 domains containing the active site aspartate residues are more structurally conserved than DPBB1. Branches are coloured according to structural fold: green, right hand (dark, cellular; and light, viral), Polβ, blue; and DPBB-containing fold, dark magenta. Since we do not believe that all polymerases originate from a common ancestor, the central node of the tree is shaded grey. The key to the additional structures is: Yeast DPBB1, yeast RNApolII DPBB from Rbp2; Yeast DPBB2, yeast RNApolII DPBB from Rbp1; Bact. DPBB1, bacterial β subunit DPBB; Bact. DPBB2, bacterial β′ subunit DPBB; Polβ, rat DNA polymerase β; T7pol, bacteriophage T7 DdRP; KF, Klenow fragment of DNA polymerase I; HIV1-RT, human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT); λ3 reo, reovirus RdRP; Φ6, Φ6 bacteriophage RdRP; HCV, hepatitis C virus RdRP; BVDV, bovine viral diarrhoea virus RdRP; PV, poliovirus RdRP; RHDV, rabbit hemorrhagic disease virus RdRP; NV, Norwalk virus RdRP; HRV, human rhinovirus RdRP; and FMDV, foot–and-mouth disease virus RdRP.
(B) Originally, in an all RNA world, RNA self-replicates until the advent of a protein-based, primeval RNA-dependent RNA polymerase. Initially this possesses a single DPBB domain on a single polypeptide chain. Gene duplication leads to a polypeptide chain containing two copies of the DPBB domain. Differentiation of the two DPBB domains then results in QDE-1–like RdRPs. Emergence of DNA and associated increase of complexity lead to segregation of the DPBB into different polypeptidic chains, giving rise to the complex multi-subunit DdRP machinery observed today.
Materials and Methods
The protein expression in yeast, purification, crystallization, and data collection of a cryo-cooled selenomethionated QDE-1 ΔN crystal (space group P21), and two native datasets (space group P21 and C2) have been described previously . Briefly, a three-wavelength MAD experiment was performed using a cryo-cooled selenomethionated crystal of QDE-1 ΔN, to a resolution of 3.2 Å, on the Medical Research Council (MRC) MAD beamline, BM14, at the European Synchrotron Radiation Facility (ESRF) (Grenoble) using a MarCCD detector (Mar USA, Evanston, Illinois, United States). Subsequently, a native dataset to 2.3-Å resolution was collected on BM14 using a MarMosaic 225 CCD detector (Mar USA). The datasets from the three-wavelength MAD experiment were scaled and merged together and used in SnB  to solve the selenium substructure. The top 54 Se sites from SnB were refined using SOLVE  to obtain initial phases, and anomalous difference Fourier maps allowed the identification of a subset of 46 correct Se atoms that were refined using SHARP . Phase improvement using density modification with RESOLVE  led to maps which, in combination with Se positions, revealed two articulated subunits, related by two slightly displaced non-crystallographic twofold axes, which precluded simple averaging. The initial RESOLVE model was completed using manual model building in O  and Coot , and automated model building and water placement with ARP/wARP . The model was refined with REFMAC5 , using TLS refinement and imposing non-crystallographic restraints for the core regions of subunits, against the native P21 data to 2.3-Å resolution to yield the final model described in the main text and Table 1. The C2 crystal form with one molecule in the asymmetric unit was solved by molecular replacement and refined, keeping the domains as rigid bodies (AMORE , Table 1).
The program ASSAM  was used to find structural similarities with the proposed three catalytic aspartates and the coordinating magnesium ion. Structures of DNA-directed RNA polymerase II largest subunit were amongst the best matches. Superimposition operators for yeast and bacterial DdRP models onto to QDE-1 ΔN were optimized using SHP .
Oligomeric state of QDE-1 ΔN.
QDE-1 ΔN size was determined by (1) gel filtration (Superdex 200 16/60, equilibrated with 25 mM Tris-HCl [pH 8.0], 200 mM NaCl) with blue dextran (2,000 kDa), catalase (232 kDa), aldolase (158 kDa), bovine serum albumin (68 kDa), ovalbumin (43 kDa), chymotrypsin A (25 kDa), and RNase A (13.7 kDa) as molecular mass markers; and (2) by sedimentation in a linear 10% to 40% sucrose gradient in 50 mM Tris-HCl (pH 8.0), 350 mM NaCl (Sorvall TH641 rotor [Thermo Electron Corporation, Waltham, Massachusetts, United States], 35,000 rpm, 42 h, 15 °C), with catalase (232 kDa), bovine serum albumin (68 kDa), phage PM2 protein P2 (30.2 kDa), and lysozyme (14 kDa) as molecular mass markers.
Figure S1. Conserved Sequence Motifs in Cellular RdRPs
Multiple sequence alignment of a representative subset of cRdRPs. Amino acid sequences of 30 cRdRPs from fungi from the groups of Ascomycota (Schizosaccharomyces pombe, Spo; Neurospora crassa, Ncr; and Gibberella zeae, Gze) and Basidiomycota (Cryptococcus neoformans, Cne), slime molds (Dictyostelium discoideum, Ddi), dicot plants (Arabidopsis thaliana, Ath; Solanum tuberosum, Stu; and Nicotiana tabacum, Ntu), monocot plants (Oryza sativa, Osa), protozoa (Entamoeba histolytica, Ehi), and nematodes (Caenorhabditis elegans, Cel) were aligned using standard settings of ClustalW algorithm. Local alignment was improved by manual editing. N. crassa QDE-1 protein sequence is shown on the top. N. crassa contains two additional non-allelic cRdRP genes—SAD-1 (essential for meiotic silencing by unpaired DNA ) and RdRP-3—that likely function in distinct cellular pathways. Invariant residues are shaded in black; other residues with 80% or more conservation are shaded in grey. Conserved sequence motifs comprising invariant residues are outlined: motif 1, red; motif 2, orange; motif 3, dark yellow; motif 4, purple; motif 5, violet; motif 6, light pink; and motif 7, blue. QDE-1 secondary structure elements are shown on top, coloured according to domain definition (slab, blue; catalytic, deep purple; neck, pink; and head, red). The identified double-psi β-barrels DPBB1 and DPBB2 are outlined by deep purple boxes. The flap sub-domain and the potential “bridge helices” are also represented by boxes, coloured light purple and grey, respectively.
(914 KB DOC)
Coordinates and structure factors have been deposited in the Protein Data Bank (PDB; http://www.rcsb.org/pdb) as accession numbers 2J7N and 2J7O. The GenPept accession numbers for the genes and gene products mentioned in this paper are N. crassa QDE-1 protein sequence (EAA29811); dRP-3 (EAA34169); and N. crassa SAD-1 (AAK31733).
E. Mancini, G. Sutton, and C. Siebold helped with synchrotron data collection, and R. Esnouf and J. Dong helped with computing. We thank the staff at the UK MAD station, BM14, ESRF, Grenoble, in particular, M. Walsh and G. Fox for assistance. We also thank P. Artymuik for structural searches with ASSAM.
PSS, MRLK, EVM, DHB, DIS, and JMG conceived and designed the experiments, analyzed the data, and wrote the paper. PSS, MRLK, EVM, and JMG performed the experiments.
- 1. Mello CC, Conte D Jr (2004) Revealing the world of RNA interference. Nature 431. pp. 338–342.
- 2. Baulcombe D (2004) RNA silencing in plants. Nature 431: 356–363.
- 3. Lippman Z, Martienssen R (2004) The role of RNA interference in heterochromatic silencing. Nature 431: 364–370.
- 4. Bernstein E, Caudy AA, Hammond SM, Hannon GJ (2001) Role for a bidentate ribonuclease in the initiation step of RNA interference. Nature 409: 363–366.
- 5. Hammond SM, Boettcher S, Caudy AA, Kobayashi R, Hannon GJ (2001) Argonaute2, a link between genetic and biochemical analyses of RNAi. Science 293: 1146–1150.
- 6. Verdel A, Jia S, Gerber S, Sugiyama T, Gygi S, et al. (2004) RNAi-mediated targeting of heterochromatin by the RITS complex. Science 303: 672–676.
- 7. Cogoni C, Macino G (1999) Gene silencing in Neurospora crassa requires a protein homologous to RNA-dependent RNA polymerase. Nature 399: 166–169.
- 8. Sijen T, Fleenor J, Simmer F, Thijssen KL, Parrish S, et al. (2001) On the role of RNA amplification in dsRNA-triggered gene silencing. Cell 107: 465–476.
- 9. Makeyev EV, Bamford DH (2002) Cellular RNA-dependent RNA polymerase involved in posttranscriptional gene silencing has two distinct activity modes. Mol Cell 10: 1417–1427.
- 10. Motamedi MR, Verdel A, Colmenares SU, Gerber SA, Gygi SP, et al. (2004) Two RNAi complexes, RITS and RDRC, physically interact and localize to noncoding centromeric RNAs. Cell 119: 789–802.
- 11. Pickford AS, Catalanotto C, Cogoni C, Macino G (2002) Quelling in Neurospora crassa. Adv Genet 46: 277–303.
- 12. Catalanotto C, Pallotta M, ReFalo P, Sachs MS, Vayssie L, et al. (2004) Redundancy of the two dicer genes in transgene-induced posttranscriptional gene silencing in Neurospora crassa. Mol Cell Biol 24: 2536–2545.
- 13. Laurila MR, Salgado PS, Makeyev EV, Nettelship J, Stuart DI, et al. (2005) Gene silencing pathway RNA-dependent RNA polymerase of Neurospora crassa: Yeast expression and crystallization of selenomethionated QDE-1 protein. J Struct Biol 149: 111–115.
- 14. Janin J, Miller S, Chothia C (1988) Surface, subunit interfaces and interior of oligomeric proteins. J Mol Biol 204: 155–164.
- 15. Krissinel E, Henrick K (2005) Detection of protein assemblies in crystals. In: Berthold MR, Glen R, Diederichs K, Kohlbacher O, Fischer I, editors. Computational life sciences: First international symposium. Berlin: Springer-Verlag. pp. 163–174.
- 16. Castillo RM, Mizuguchi K, Dhanaraj V, Albert A, Blundell TL, et al. (1999) A six-stranded double-psi b-barrel is shared by several protein superfamilies. Structure 7: 227–236.
- 17. Artymiuk PJ, Poirrette AR, Grindley HM, Rice DW, Willett P (1994) A graph-theoretic approach to the identification of three-dimensional patterns of amino acid side-chains in protein structures. J Mol Biol 243: 327–344.
- 18. Cramer P, Bushnell DA, Kornberg RD (2001) Structural basis of transcription: RNA polymerase II at 2.8 angstrom resolution. Science 292: 1863–1876.
- 19. Vassylyev DG, Sekine S, Laptenko O, Lee J, Vassylyeva MN, et al. (2002) Crystal structure of a bacterial RNA polymerase holoenzyme at 2.6 A resolution. Nature 417: 712–719.
- 20. Iyer LM, Koonin EV, Aravind L (2003) Evolutionary connection between the catalytic subunits of DNA-dependent RNA polymerases and eukaryotic RNA-dependent RNA polymerases and the origin of RNA polymerases. BMC Struct Biol 3: 1.
- 21. Gnatt AL, Cramer P, Fu J, Bushnell DA, Kornberg RD (2001) Structural basis of transcription: An RNA polymerase II elongation complex at 3.3 A resolution. Science 292: 1876–1882.
- 22. Westover KD, Bushnell DA, Kornberg RD (2004) Structural basis of transcription: Nucleotide selection by rotation in the RNA polymerase II active center. Cell 119: 481–489.
- 23. Lai MM (2005) RNA replication without RNA-dependent RNA polymerase: Surprises from hepatitis delta virus. J Virol 79: 7951–7958.
- 24. Schiebel W, Pelissier T, Riedel L, Thalmeir S, Schiebel R, et al. (1998) Isolation of an RNA-directed RNA polymerase-specific cDNA clone from tomato. Plant Cell 10: 2087–2101.
- 25. Bressanelli S, Tomei L, Roussel A, Incitti I, Vitale RL, et al. (1999) Crystal structure of the RNA-dependent RNA polymerase of hepatitis C virus. Proc Natl Acad Sci U S A 96: 13034–13039.
- 26. Butcher SJ, Grimes JM, Makeyev EV, Bamford DH, Stuart DI (2001) A mechanism for initiating RNA-dependent RNA polymerization. Nature 410: 235–240.
- 27. Choi KH, Groarke JM, Young DC, Kuhn RJ, Smith JL, et al. (2004) The structure of the RNA-dependent RNA polymerase from bovine viral diarrhea virus establishes the role of GTP in de novo initiation. Proc Natl Acad Sci U S A 101: 4425–4430.
- 28. Lesburg CA, Cable MB, Ferrari E, Hong Z, Mannarino AF, et al. (1999) Crystal structure of the RNA-dependent RNA polymerase from hepatitis C virus reveals a fully encircled active site. Nat Struct Biol 6: 937–943.
- 29. Ferrer-Orta C, Arias A, Perez-Luque R, Escarmis C, Domingo E, et al. (2004) Structure of foot-and-mouth disease virus RNA-dependent RNA polymerase and its complex with a template-primer RNA. J Biol Chem 279: 47212–47221.
- 30. Tao Y, Farsetta DL, Nibert ML, Harrison SC (2002) RNA synthesis in a cage—structural studies of reovirus polymerase lambda3. Cell 111: 733–745.
- 31. Zakharova N, Hoffman PS, Berg DE, Severinov K (1998) The largest subunits of RNA polymerase from gastric helicobacters are tethered. J Biol Chem 273: 19371–19374.
- 32. Severinov K, Mooney R, Darst SA, Landick R (1997) Tethering of the large subunits of Escherichia coli RNA polymerase. J Biol Chem 272: 24137–24140.
- 33. Weeks CM, Miller R (1999) The design and implementation of SnB 2.0. J Appl Cryst 32: 120–124.
- 34. Terwilliger TC (2000) Maximum-likelihood density modification. Acta Crystallogr D Biol Crystallogr 56(Pt 8): 965–972.
- 35. de la Fortelle E, Bricogne G (1997) Maximum-likelihood heavy-atom parameter refinement for multiple isomorphous replacement and multiwavelength anomalous diffraction methods. Methods Enzymol 276: 472–494.
- 36. Jones TA, Zou JY, Cowan SW, Kjeldgaard M (1991) Improved methods for building protein models in electron density maps and the location of errors in these models. Acta Crystallogr A 47 ( Pt 2). pp. 110–119.
- 37. Emsley P, Cowtan K (2004) Coot: Model-building tools for molecular graphics. Acta Crystallogr D Biol Crystallogr 60: 2126–2132.
- 38. Perrakis A, Harkiolaki M, Wilson KS, Lamzin VS (2001) ARP/wARP and molecular replacement. Acta Crystallogr D Biol Crystallogr 57: 1445–1450.
- 39. Murshudov GN, Vagin AA, Dodson EJ (1997) Refinement of macromolecular structures by the maximum-likelihood method. Acta Crystallogr D Biol Crystallogr 53: 240–255.
- 40. Navaza J (1994) AMoRe: An automated package for molecular replacement. Acta Crystallographica A 50: 157–163.
- 41. Stuart DI, Levine M, Muirhead H, Stammers DK (1979) Crystal structure of cat muscle pyruvate kinase at a resolution of 2.6 A. J Mol Biol 134: 109–142.
- 42. Gouet P, Courcelle E, Stuart DI, Metoz F (1999) ESPript: Analysis of multiple sequence alignments in PostScript. Bioinformatics 15: 305–308.
- 43. Bamford DH, Grimes JM, Stuart DI (2005) What does structure tell us about virus evolution? Curr Opin Struct Biol 15: 655–663.