Research Article

Structural basis of ribosomal frameshifting during translation of the SARS-CoV-2 RNA genome

See allHide authors and affiliations

Science  13 May 2021:
eabf3546
DOI: 10.1126/science.abf3546

Abstract

Programmed ribosomal frameshifting is a key event during translation of the SARS-CoV-2 RNA genome allowing synthesis of the viral RNA-dependent RNA polymerase and downstream proteins. Here we present the cryo-electron microscopy structure of a translating mammalian ribosome primed for frameshifting on the viral RNA. The viral RNA adopts a pseudoknot structure that lodges at the entry to the ribosomal mRNA channel to generate tension in the mRNA and promote frameshifting, whereas the nascent viral polyprotein forms distinct interactions with the ribosomal tunnel. Biochemical experiments validate the structural observations and reveal mechanistic and regulatory features that influence frameshifting efficiency. Finally, we compare compounds previously shown to reduce frameshifting with respect to their ability to inhibit SARS-CoV-2 replication, establishing coronavirus frameshifting as a target for antiviral intervention.

Ribosomal frameshifting, a process during which the reading frame of translation is changed at the junction between open reading frames 1a and 1b, is one of the key events during translation of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) positive sense single-stranded RNA genome. This programmed -1 translational frameshifting is conserved in all coronaviruses and is necessary for synthesis of viral RNA-dependent RNA polymerase (RdRp or Nsp12) and downstream viral non-structural proteins encoding core enzymatic functions involved in capping of viral RNA, RNA modification and processing, and RNA proof-reading (1). Although the translational machinery typically prevents frameshifting as a potential source of one of the most disruptive errors in translation (2, 3), many viruses rely on programmed ribosomal frameshifting to expand and fine-tune the repertoire and stoichiometry of expressed proteins (4).

Programmed -1 frameshifting in SARS-related coronaviruses occurs at the slippery sequence U_UUA_AAC in the context of a 3′ stimulatory RNA sequence that was predicted to form a 3-stemmed pseudoknot structure (5), and in parallel was independently tested by our lab and others (68). The frameshifting occurs with high efficiency (25-75%) depending on the system used (6, 7, 911) and changes the reading frame to UUU_AAA_C (12) (Fig. 1A). Consequently, two viral polyproteins are synthesized, one encoded by the ORF1a when frameshifting does not take place, whereas ORF1ab is expressed as a result of frameshifting. Translation of ORF1a produces polyprotein 1a ending with Nsp10 followed by the short Nsp11. On the other hand, when the frameshift occurs, the polyprotein 1ab is generated, which contains almost 2700 additional amino acids and in which the viral RdRp, Nsp12, is produced after Nsp10 as a consequence of translation in the -1 frame. A putative secondary structure element in the viral RNA that forms a loop upstream of the shift site has been proposed to play an attenuating role in frameshifting and is referred to as the 5′ attenuator loop (8). Maintaining the precise level of coronavirus frameshifting efficiency is crucial for viral infectivity, evidenced by the remarkable fact that mutation of a single nucleotide in the frameshifting region of the SARS-CoV-1 RNA results in a concomitant abrogation of viral replication (13). Therefore, the importance of 3-stemmed pseudoknot-dependent -1 ribosomal frameshifting for the propagation of SARS-related coronaviruses, a process that has not been seen to occur on any endogenous human transcript in human cells, presents itself as an opportune drug-target with minimal tolerance for drug-resistant mutations.

Fig. 1 The SARS-CoV-2 pseudoknot interacts with the ribosome and pauses translation upstream of the slippery site.

(A) Schematic of the SARS-CoV-2 main ORF. In the close up view of the frameshift event, codons and corresponding amino acids are shown. During -1 frameshifting, the ‘slippery site’ codons UUA (Leu) and AAC (Asn) are the last codons decoded in the 0 frame. Upon -1 frameshifting of the AAC codon to AAA, translation resumes at the CGG (Arg) triplet, where elongation proceeds uninterrupted to produce full-length Nsp12. (B) In vitro translation reaction depicting pausing at the frameshift site. Efficient frameshifting is observed for the WT template, consistent with our dual luciferase assays (see methods). Samples for cryo-EM originally intended to be trapped by dominant negative eRF1 (AAQ) show a tRNA-bound pause in proximity of the frameshift site. The tRNA-associated band is lost upon RNase treatment. Reactions without added eRF1 (AAQ) produce a similarly paused product. (C) Overview of the density low pass filtered to 6Å with the pseudoknot found close to the entry of the mRNA channel on the small subunit (SSU). The SSU proteins are colored in yellow, the large subunit (LSU) proteins in blue and the rRNA in grey. The pseudoknot is colored according to its secondary structure as in (F), and the P-site tRNA is colored in green. (D) Close-up view of the pseudoknot from the solvent-exposed side of the SSU. Helix h16 of the 18S rRNA interacts with the base of Stem 1. Unpaired loop-forming nucleotides are colored in cyan. (E) P-site codon-anticodon interactions reveal a Phe (UUU) codon interacting with tRNA(Phe). (F) Schematic of the revised secondary structure elements in the pseudoknot necessary for -1 PRF with different functional regions labeled and colored accordingly.

Due to its importance in the life cycle of many important viruses and coronaviruses in particular, programmed frameshifting has been extensively studied using a range of structural and functional approaches (4). The structure of a 3′ stimulatory pseudoknot in isolation or in context of the viral genome has been proposed recently by various groups using techniques that include molecular dynamics, nuclease mapping, in vivo selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE), nuclear magnetic resonance (NMR) and cryo-electron microscopy (cryo-EM) (7, 1417). Furthermore, a ribosomal complex with a frameshift stimulatory pseudoknot from the avian infectious bronchitis virus was reported at low resolution (18). Here, to provide a structural and mechanistic description of the events during ribosomal frameshifting, we investigated mammalian ribosomes captured in distinct functional states during translation of a region of SARS-CoV-2 genomic RNA where -1 programmed frameshifting occurs.

Structure determination of a frameshifting-primed ribosomal complex

We captured a 0 frame, pre-frameshift ribosomal complex by introducing a stop codon in place of the second codon of the slippery site (U_UUA_AAC to U_UUA_UAA) (Fig. 1A) and adding mutant eukaryotic Release Factor 1 [eRF1 (AAQ)] that is unable to release the nascent polypeptide. Translating complexes were prepared in an in vitro translation reaction using an in-house generated rabbit reticulocyte lysate (RRL) system that supported efficient frameshifting in the previously reported range of around 50% (19) according to dual luciferase experiments (see methods). The ribosomes were programmed with mRNA encoding an affinity tag and harboring a region of the SARS-CoV-2 genome that encodes proteins Nsp10 (C terminus), Nsp11 and the majority of Nsp12. Western blotting showed that when using the WT RNA template, frameshifting was efficient, while the stop codon mutation prevented frameshifting and led to ribosome pausing. This effect was further enhanced when eRF1 (AAQ) was present in excess over endogenous wild type eRF1 (Fig. 1B).

The cryo-EM 3D reconstruction of ribosome-nascent chain complexes (RNCs) affinity purified from the reactions supplemented with eRF1 (AAQ) revealed two distinct ribosomal complexes captured in the process of translating the slippery sequence (figs. S1 and S2). One represented a termination complex that contained the ATP-binding cassette transporter 1 (ABCE1) known to be involved in termination and recycling together with mutant eRF1 interacting with the stop codon (fig. S3). The second reconstruction resolved translating 80S ribosomes containing P- and E-site tRNAs bound (fig. S2). This reconstruction at 2.2 Å resolution allowed us to build the most accurate structure of a mammalian 80S ribosome so far and directly visualize many protein and virtually all rRNA modifications identified for the human ribosome based on quantitative mass spectrometry and as interpreted in a recent human ribosome structure (20, 21), consistent with the complete conservation of all modified residues between rabbit and human rRNAs (figs. S4 and S5; and tables S1 to S3). Importantly, this reconstruction also featured additional density at the entrance to the mRNA channel suggestive of a structured RNA, which after focused classification revealed a prominent density for a complete 3′ frameshifting stimulatory pseudoknot at the entry of the mRNA channel on the 40S subunit (Fig. 1, C and D). The resolution of this reconstruction ranged from 2.4 Å at the core of the ribosome to ~7 Å at the periphery, where the most flexible regions of the pseudoknot are located (figs. S2 and S6). Based on the high-resolution maps that allowed visualization of the codon-anticodon interactions and modifications in the tRNA (Fig. 1E and fig. S6, A and B), we could unequivocally determine that a Phe-tRNA(Phe) was bound at the P-site (22). The mRNA does not adopt any unusual structure in the A-site of the ribosome as was observed for the HIV-1 frameshifting sequence visualized on the bacterial ribosome (23). This implied that the ribosome is paused by the downstream pseudoknot located at the entrance to the mRNA channel such that the P-site tRNA interacts with the UUU codon just prior to the first codon, UUA, of the slippery site (Fig. 2A).

Fig. 2 Critical features of the ribosome-bound pseudoknot.

(A) Overview of the frameshift-primed state. The stimulatory pseudoknot pauses the ribosome at the penultimate codon (UUU) of the slippery site (red), with P- (green) and E- (pink) sites occupied by tRNAs, and an empty A-site awaiting decoding in the non-rotated state. The length of the spacer region (grey) is critical for exact positioning of the pseudoknot as the spacer exerts tension at the entry of the mRNA channel (fig. S6C). (B) The backbone of Loop 1 (UGC) (cyan) of the pseudoknot interacts with the N-terminal domain of uS3 (red) and the C-terminal tail of eS10 (orange). mRNA residue G13486 is flipped out and interacts with uS3 (fig. S6D). (C) Mutagenesis experiments using dual luciferase assays in HEK293T cells indicate that the G13486 interaction is specific. Mutation of G13486 to other residues leads to a marked reduction in frameshifting efficiency, and deletion of Loop 1 (Δ L1) completely abolishes frameshifting. Similarly, deletion of a single nucleotide (A13537) in Loop 2 reduces frameshifting, while deletion of the entire loop (ΔL2) abolishes frameshifting. Normalized (Firefly/Renilla) luciferase activities were calculated for each construct as a percentage of their individual normalized in-frame controls. Data are presented as mean values ± standard deviations of three biological replicates (sets of translation reactions) averaged after three measurements, with error bars representing standard deviations. (D) Mutagenesis experiments using dual luciferase reporter assays in HEK293T cells show that the position of the 0 frame stop codon influences frameshifting. Leaving the pseudoknot unaltered, incremental increase in the distance of the 0 frame stop codon from the frameshift site leads to a concomitant decrease in frameshifting levels. Loss of the stop codon in 0 frame leads to a sharp decline in frameshifting levels. This reduction is rescued by ~45% upon decreasing ribosome loading levels by implementing weaker initiation codons. The graph is normalized relative to the WT frameshifting of 25%. Mutations and complementary mutations are shown in fig. S8.

The pseudoknot causes ribosomal pausing prior to -1 frameshifting

The observation that the pseudoknot acts as an obstacle to slow down translation as the ribosome approaches the slippery site is mechanistically reasonable. Since the pseudoknot is a stable structural element in the mRNA, it will resist unfolding and consequently generate a back-pull on the viral RNA, resulting in an increased chance of -1 frameshifting as the tRNAs are translocated. A pause in translocation at a codon that precedes the slippery site, characterized by a >10 times longer occupancy prior to the slippage event, was observed in an analogous case of heptanucleotide -1 frameshifting on the bacterial dnaX gene using single molecule experiments (24). According to this model, it would be anticipated that a further round of translocation results in unwinding of Stem 1 of the downstream stimulatory pseudoknot structure. Consistently, in our structure of the eRF1 (AAQ)-bound ribosome that advanced one codon further along the mRNA, no clear secondary structure is visible at the entrance to the mRNA channel as the mRNA now becomes disordered at this position (figs. S1 and S3, A and B).

In order to investigate the slowdown of translation on the wild type slippery sequence, we performed disome footprint profiling, a method to identify translational pause sites through the analysis of transitory ribosome collisions (2527) (see methods). Notably, recent studies using conventional ribosome profiling methodology reported a lack in monosome footprint coverage across the frameshifting region on the SARS-CoV-2 RNA (11, 28) – possibly because ribosomes in this area became trapped in temporary collisions. Moreover, the highly structured pseudoknot at the entry to the mRNA channel would likely preclude efficient trimming by RNase I, the enzyme used for footprint generation, further reducing efficient monosome footprint capture. Using a modified nuclease treatment protocol (see methods) that recovered monosome footprints from the frameshift region (Fig. 3, A and C), our experiments revealed that ribosome collisions occur as a result of ribosomal pausing at the same position that is observed in the structure of the pseudoknot-engaged ribosome (Fig. 3, B and D). Apparently, although the base substitutions creating a stop codon in the 3′ adjacent slippery site did not change the features of pausing, it increased the dwell time of the ribosomes at the pause site sufficiently to allow visualization in the cryo-EM experiment.

Fig. 3 Pseudoknot-mediated pause occurs prior to the -1 frameshifting event.

(A and B) Footprint coverage for WT and mutant constructs determined by monosome- and disome-selective ribosome profiling. Pileup of reads from the indicated areas are plotted separately for reads that overlap (pink) vs. do not overlap (grey) the frameshift site (indicated by red bar below x-axis). The predicted A-sites of the ribosomes giving rise to the footprints are depicted as blue peaks. A-site predictions were carried out as described in Supplementary Material. (C) Zoom into the frameshift region from (A) and (B) reveals that monosome profiles show transient occupancy in the vicinity of the frameshift site, while disome profiles – indicative of strong pause sites – show a similarly enhanced occupancy at the first codon (UUA) of the frameshift site in both WT and mutant constructs. A-site codons of the leading and trailing ribosome are highlighted with a translucent bar and correspond to those seen in the disome structure in (D). (D) In high resolution cryo-EM reconstructions of disomes at the frameshift site, the P- and A-sites of the trailing ribosome show occupancy of CCC and AUG codons, respectively, corresponding to the positions estimated by disome profiling. Disome maps were calculated by separately refining the orientational parameters for each ribosome.

The results of our disome profiling experiments prompted us to structurally investigate disomes by cryo-EM. We were able to visualize the pseudoknot-paused ribosome followed by a closely trailing ribosome. Upon focused refinement, we obtained a high-resolution (3.1 Å) structure of the trailing ribosome in a rotated state (fig. S1). In congruence with our estimated positioning of the ribosomes in disome profiling (Fig. 3D), the purine-pyrimidine pattern of codon-anticodon pairs in the structure of the colliding ribosome revealed that the pause occurs with CCC and AUG triplets in the P- and A-sites, respectively (Fig. 3C).

The SARS-CoV-2 RNA pseudoknot specifically interacts with ribosomal proteins and 18S rRNA

The intermediate local resolution (5-7 Å) of the cryo-EM map in the area of the pseudoknot allowed us to visualize the overall fold of the RNA and readjust its previously predicted secondary structure (1417, 19) (Fig. 1, C, D, and F). The stimulatory pseudoknot forms an H-type pseudoknot with Stem 1 and Stem 2 coaxially stacked on top of each other to form a quasi-continuous helix, while Stem 3 stands out almost perpendicular to this plane (Figs. 1D and 2B). This corkscrew-like formation provides a bulky and well-structured obstacle wedged at the mRNA entry channel, having the potential to resist unwinding by the helicase activity of the ribosome and generating tension on the upstream mRNA up to the decoding center. Stem 1 of the pseudoknot forms a 9 bp helix which is GC rich at the bottom (Fig. 1F). The penultimate nucleotides of the ‘spacer region’ prior to Stem 1 are located at the mRNA entry tunnel, where they interact with several basic residues in the C-terminal domain of uS3 on one side and are supported by uS5 from the other, with an additional weak contact contributed by the C-terminal end of eS30. uS3 and eS30 are primary components of the ribosome helicase and uS5 has been proposed to be a component of the ribosomal helicase processivity clamp at the mRNA entry site (29, 30). The observed distance between the P-site UUU codon and Stem 1 of the pseudoknot underscores the critical dependence of the frameshifting efficiency on the length of the spacer region (31). Translocation to the next codon would place the frameshifting codon UUA into the P-site, with a simultaneous increase in the tension of the mRNA and unwinding of the GC-rich base of Stem 1 upon entering into the mRNA entry channel, comparable to the situation when the ribosome proceeds to the engineered stop codon as observed in our eRF1 (AAQ)-stalled structure (fig. S3).

The pseudoknot structure also reveals a hitherto unobserved and possibly unappreciated role for the distal site of the mRNA entrance channel in helicase activity. While mRNA unwinding studies outside the mRNA entrance channel have so far implicated only a helix in the C-terminal domain of uS3 (32), we notice that Loop 1 of the pseudoknot contacts the N-terminal domain of uS3 as well as the C-terminal tail of eS10 (Fig. 2B and fig. S6D), whereas the flipped-out base G13486 in this loop forms specific interactions (Fig. 2B). Furthermore, as the pseudoknot is located at the entry to the mRNA channel, helix h16 of the 18S rRNA is noticeably pushed outwards due to a direct contact with the minor groove of Stem 1 (Fig. 2B and fig. S7A). Since the pseudoknot wedges between the head and the body of the small ribosomal subunit, it would restrict their relative motions that need to take place during translocation. This is consistent with the studies on dynamics of coronavirus frameshifting, which revealed that the mechanism of -1 frameshifting involves restriction of small subunit head motion (33).

The structure also reveals another key aspect of the architecture of the pseudoknot as the ribosome encounters it. The start of the pseudoknot is shifted relative to the predicted secondary structure (1417, 19) by two nucleotides. The two opposed nucleotides, which were assumed to base pair with Stem 1, are actually forming the start of Stem 3 by pairing with bases predicted to be in the single-stranded linker 2 (Fig. 1F and fig. S7, B and C). Our cryo-EM density reveals that Loop 3 accommodates a total of 4 nucleotides, three of which were originally attributed to Stem 2. Thus, we observe that Loop 3 is shifted and expanded relative to the initially predicted secondary structures (1417, 19).

To functionally support our structural findings and confirm the nature and specificity of the pseudoknot interactions, we performed structure-guided mutagenesis experiments using dual luciferase reporter assays in HEK293T cells (see methods) and monitored the frameshifting efficiency relative to the WT (Fig. 2C). Mutation of G13486 of Loop 1 to another purine reduced the frameshifting efficiency to 30% of the WT level, and mutation of this base to a pyrimidine further reduced frameshifting to 15%. As expected from our structural data, deletions of the nucleotides of the spacer regions also had a deteriorating effect on frameshifting. Loss of Loop 1 entirely abolished frameshifting. Deletion of a single nucleotide of Loop 3 in agreement with its proposed role in forming the base pairing interactions diminished the frameshifting rate to 25% of the WT level. Loss of the entire Loop 3 reduced frameshifting to 10% of WT levels.

Frameshifting efficiency depends on the position of the “0” frame stop codon

In SARS-CoV-2, the 0 frame stop codon is located 5 codons downstream of the frameshift site and is a constituent of Stem 1. The placement of the stop codon in such proximity to the frameshift site is a common feature in coronaviruses, and its presence in a critical region of the stimulatory pseudoknot prompted us to probe the effect of the distance of the 0 frame stop codon on frameshifting. To this end, knowledge of the 3D structure of the pseudoknot helped us to confidently manipulate the stop codon without hampering pseudoknot formation. We introduced mutations to incrementally extend the stop codon from the WT position and to completely remove the occurrence of a stop codon in the 0 frame (Fig. 2D and fig. S8). While introducing a stop codon 6 nucleotides downstream of the WT position only marginally decreased the frameshifting rate (98% of WT), a stronger attenuation was observed when the distance of the stop codon was increased to 15 nucleotides from the WT stop (80% of WT). Finally, removal of the stop codon by two different point mutations led to a reduction of frameshifting efficiency to 50% of WT levels. To test whether reduced ribosomal loading rescues the effect of stop codon removal, we analyzed the frameshifting efficiency in the context of weaker initiation codons such as CUG and AUU (Fig. 2D). These constructs led to a 45% rescue of the reduction in frameshifting compared to stop codon mutants initiating at an AUG start.

Taken together, these observations suggest that the stop codon position plays an important role in maintaining optimum frameshift efficiency. We propose that the stop codon serves to prevent the closely trailing ribosome from encountering a viral RNA that was unfolded by the leading ribosome. In this case, upon encountering a stop codon, termination and subunit disassembly will occur, which will provide an opportunity for the pseudoknot to refold without the constraints of the mRNA channel (see Conclusions). According to this model, although the wild type stop codon will make the frameshifting efficiency less sensitive to ribosome loading in the “no-frameshifting” scenario, the frameshifting events that occur following a -1 frameshift will nevertheless be more likely when the ribosomes are spaced further apart. Our measurements of the efficiency of frameshifting for the wild type sequence in the context of different rates of translation initiation are in agreement with this hypothesis (fig. S9). This mechanism, consistent with our biochemical data, increases the efficiency of frameshifting to the levels required by SARS-CoV-2 and may be used by viruses in general when high-efficiency frameshifting is required.

Nascent chain forms specific interactions with the ribosomal exit tunnel

Strikingly, in the reconstruction of the paused translating ribosome, the nascent chain that corresponds to the viral polyprotein was visible along the entire length of the ribosomal exit tunnel (Fig. 4A). The density corresponded to the C-terminal region of Nsp10, which is the activator of the viral proofreading exonuclease and N7-methyltransferase Nsp14 (34, 35), and then (depending on the frameshifting event) continues as either the viral RNA-dependent RNA polymerase Nsp12 (6) or as protein Nsp11, whose function is yet unknown (Figs. 1A and 4B). The nascent chain makes several specific interactions with the ribosomal tunnel, one of which is at the constriction site where arginine 4387 of Nsp10 interacts with A1555 of the 28S rRNA [corresponding to A1600 in humans, numbering according to PDB 6EK0 (36)] and is stabilized by the preceding leucine 4386 (Fig. 4C). Notably, these two amino acids are highly conserved across multiple coronaviruses (Fig. 4G), although they are located in the unstructured C-terminal region of Nsp10 and therefore considered not to be important for the fold of the protein (37).

Fig. 4 The nascent viral polypeptide co-translationally folds and specifically interacts with the ribosomal tunnel.

(A) Cross-section of the pseudoknot-paused ribosome structure showing the exit tunnel. The nascent C terminus of Nsp10 (orange) and the N terminus of Nsp11/12 (purple) are visible from the PTC to the periphery of the ribosome exit tunnel (LSU in blue). (B) Schematic representation of the path of nascent peptide along the exit tunnel. Arg 4387 stacks with 28S rRNA residue A1555 at the constriction site. Further down, where the tunnel widens, the C-terminal zinc finger domain of Nsp10 folds co-translationally, with Trp 4376 stacking on A2261 of 28S rRNA. (C) Well-ordered density is visible for Arg 4387 of Nsp10 as it stacks onto A1555 of 28S rRNA at the constriction site and is stabilized by Leu 4386. The structure is shown within the cryo-EM map contoured at two different levels (grey and red). (D and E) The overlay of the co-translationally folded zinc finger domain with the crystal structure of Nsp10 [green, PDB 2FYG (37)] reveals the structural similarity. (F) Probing the role of nascent chain interactions with the ribosome exit tunnel using an RRL in vitro system. Mutations of the interacting residues were tested for their effect on frameshifting shown in comparison to the wild type frameshifting (41% frameshifting was normalized to 100%). Replacement of the entire nascent chain with an unrelated sequence leads to a 35% relative increase in frameshifting, which is only in part due to loss of the 5′ attenuator loop. Interactions around the constriction site likely serve to attenuate frameshifting, as replacement of the interacting Arg 4387 and stabilizing Leu 4386 with Ala increases frameshifting by 30%. (G) Alignment of SARS2 with closely related sequences of other coronaviruses highlighting the conservation of the mutated residues [colored as in (F)]. The shown sequence stretch encompasses the C-terminal zinc finger domain of Nsp10 (orange) and parts of Nsp11/Nsp12 (purple) visible in our reconstruction. Nascent chain residues Leu 4386 and Arg 4387 interacting with the ribosomal exit tunnel are strictly conserved, while the conservation of neighboring residues is lower. Stars represent the four cysteines of the Nsp10 zinc finger.

Further down the tunnel, the C-terminal end of Nsp10 adopts a partially folded zinc finger motif (Fig. 4, D and E), which upon superposition reveals similarity with the corresponding fully folded C-terminal domain previously observed in the crystal structure of SARS-CoV-1 Nsp10 (37). Tryptophan 4376 located between the two pairs of cysteines that form the zinc finger stacks with A2261 (A2418), an interaction that might serve to promote the change of nascent chain direction and facilitate folding of the zinc finger at the end of the exit tunnel. Co-translational events, such as insertion of a transmembrane domain at the exit of the ribosomal tunnel, was shown to promote -1 ribosomal frameshifting in alphaviruses (38).

To investigate whether the observed contacts between the nascent chain and the ribosomal tunnel are specific, and whether these interactions and co-translational folding of Nsp10 might play a role in modulating the frameshifting process, we employed our dual luciferase reporter assay to measure the frameshifting efficiency of WT and mutant nascent chain sequence constructs. As our measurements in HEK293T cells did not reveal an appreciable change of frameshift efficiency, we carried out the same experiments in vitro using RRL to monitor the effects in a single mRNA setup. Replacement of the entire nascent chain with an unrelated sequence leads to a 35% increase in frameshifting (Fig. 4F). Importantly, this effect was provoked by the change in peptide sequence and not simply by the loss of the 5′ attenuator loop, given that a reporter containing silent attenuator loop mutations resulted in only a slight increase in frameshifting (Fig. 4F). Mutation of the leucine 4386 and arginine 4387 to alanine led to a considerable (30%) increase in frameshifting (Fig. 4, F and G), implying that these nascent chain interactions with the ribosomal exit tunnel play an important role in regulating frameshifting levels, possibly mechanistically akin to the well-studied SecM stalling system in bacteria (39), where it was shown that co-translational folding and the translocon-induced mechanical force can rescue the stall induced by interactions between the nascent chain and the ribosomal tunnel (40). These observations also suggest that any cellular nascent-chain factors (41, 42) might influence the rate of frameshifting.

Inhibition of viral replication by a compound that targets the SARS-CoV-2 pseudoknot

The sensitivity of the coronavirus to the finely controlled frameshifting levels (13) may present an opportunity to develop compounds that interfere with the frameshifting process and thus inhibit replication of the virus. Using computational modeling and reporter assays, compounds that have been predicted to bind the pseudoknot and inhibit SARS-CoV-2 frameshifting were described (19, 43), but never tested with respect to their ability to inhibit viral replication. Furthermore, the fluoroquinolone compound merafloxacin was recently reported to also inhibit -1 frameshifting efficiency of SARS-CoV-2 and other betacoronaviruses (44). To demonstrate that the inhibition of frameshifting is a plausible strategy for drug development, we compared two of the previously described compounds with respect to their ability to reduce viral levels in infected African green monkey VeroE6 cells (fig. S10 and methods). Our experiments show that merafloxacin is a better candidate compound as it showed a concentration dependent inhibition of frameshifting, whereas, contrary to earlier reports (19, 43) MTDB under our experimental conditions did not specifically inhibit frameshifting (fig. S10). The two compounds showed no cellular toxicity and resulted in a 3 to 4 orders of magnitude reduction of SARS-CoV-2 titer, with the half maximal inhibitory concentration (IC50) of 48 μΜ for MTDB and an order of magnitude higher efficacy of merafloxacin with an IC50 of 4.3 μΜ (fig. S10). Since MTDB did not appear to affect frameshifting in our reporter construct experiments in vitro and in vivo, it is possible that it inhibits SARS-CoV-2 replication by a different mechanism. Although the potency range for these compounds is not what would be expected from potential drug candidates, it nevertheless provides a starting point for high-throughput screening and establishes that frameshifting is a viable target for therapeutic intervention against SARS-CoV-2.

Conclusions

Our results provide a mechanistic description of frameshifting that occurs during translation of the SARS-CoV-2 genome and reveal the features that may be exploited by the virus to finely control the stoichiometry of viral proteins at different stages of infection (Fig. 5). Interfering with the frameshifting process at the level of nascent chain interactions with the ribosomal tunnel, at the level of RNA folding that leads to the formation of the frameshift stimulatory pseudoknot, or to perturb the interactions between the pseudoknot and the mRNA channel, represent a viable strategy in our search for new drugs against SARS-CoV-2, the virus that is currently causing the global COVID-19 pandemic. Our results will also be useful for understanding the mechanism of programmed ribosomal “-1” frameshifting (4) including that employed by many other medically important viruses.

Fig. 5 Structure based model for -1 programmed frameshifting in coronaviruses and its regulation.

The observed interactions between the pseudoknot and the ribosome prime the system for frameshifting. The features of the pseudoknot and the interactions between the nascent chain and the ribosomal tunnel play a role in the efficiency of frameshifting. The efficiency of frameshifting is increased by the presence of a stop codon near the frameshifting site. Ribosomes that progress beyond the frameshifting site in the 0 frame quickly terminate and disassemble, thereby increasing the chances that the pseudoknot will refold before it is encountered by the closely trailing ribosome. The trailing ribosome in turn encounters the pseudoknot, which increases the possibility of undergoing -1 frameshifting.

Supplementary Materials

science.sciencemag.org/cgi/content/full/science.abf3546/DC1

Materials and Methods

Figs. S1 to S11

Tables S1 to S4

References (4478)

MDAR Reproducibility Checklist

https://creativecommons.org/licenses/by/4.0/

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 work is properly cited.

References and Notes

Acknowledgments: We thank A. Jomaa for advice on cryo-EM data processing, Noreen Casey for support with technical resources and A. Picenoni for help with the grid preparation. We are indebted to the ETH scientific center for optical and electron microscopy (ScopeM) for access to electron microscopes, and in particular to M. Peterek and D. Boehringer. We thank Junjie Guo for helpful discussions and for providing merafloxacin compound. Funding: This work was supported by the Swiss National Science Foundation (SNSF) to NB and VT (grant 31003A_182341 and 310030_173085) via the National Center of Excellence in RNA and Disease (project funding 182880) to N. Ban, V. Thiel and D. Gatfield; The Federal Ministry of Education and Research, Germany (BMBF; grant RAPID, #01KI1723A) to AK and VT; The Irish Research Council Advanced Laureate (IRCLA/2019/74) to J.F.A., who thanks Prof. G.F. Fitzgerald for funds from Carbery Group Ltd. Author contributions: PRB and NB initiated the project, and together with JFA designed the experiments. PRB carried out biochemical experiments and sample preparation. AS and PRB prepared grids. AS carried out data collection and calculated EM maps. AS and ML performed molecular model building and refinement. AS, ML, PRB and NB interpreted the structure. PRB, AS, ML and NB drafted the manuscript. AS prepared the figures. GL and KMOC performed dual-luciferase reporter assays. AM prepared and purified the drug-like ligand MTDB. AK and VT formulated experiments to test viral infected cells. AK carried out testing of ligand on viral-infected cells. RM, RD and DG designed, carried out, analyzed and interpreted ribosome profiling experiments. All authors contributed to the final version of the manuscript. Competing interests: The authors declare no competing interest. Data and materials availability: The structure and cryo-EM map for the high-resolution reconstruction is available in the Protein Data Bank (PDB) as PDB ID 7O7Y and in the Electron Microscopy Data Bank (EMDB) as EMD-12756, respectively. The structure refined into the further classified particle set reconstruction and the corresponding maps are available as PDB-7O7Z and EMD-12757, respectively. The structure of the ribosome bound with eRF1 (AAQ) and ABCE1 reconstruction and the corresponding maps are available as PDB-7O80 and EMD-12758, respectively. The structure of the colliding ribosome reconstruction and the corresponding maps are available as PDB-7O81 and EMD-12759, respectively. Structure of the composite map of the disome is available in the EMDB as EMD-12760. Ribosome profiling data have been deposited at GSE167421. This work is licensed under a Creative Commons Attribution 4.0 International (CC BY 4.0) license, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. To view a copy of this license, visit https://creativecommons.org/licenses/by/4.0/. This license does not apply to figures/photos/artwork or other content included in the article that is credited to a third party; obtain authorization from the rights holder before using such material.

Stay Connected to Science

Subjects

Navigate This Article