Modeling the structure of the frameshift-stimulatory pseudoknot in SARS-CoV-2 reveals multiple possible conformers

The coronavirus causing the COVID-19 pandemic, SARS-CoV-2, uses −1 programmed ribosomal frameshifting (−1 PRF) to control the relative expression of viral proteins. As modulating −1 PRF can inhibit viral replication, the RNA pseudoknot stimulating −1 PRF may be a fruitful target for therapeutics treating COVID-19. We modeled the unusual 3-stem structure of the stimulatory pseudoknot of SARS-CoV-2 computationally, using multiple blind structural prediction tools followed by μs-long molecular dynamics simulations. The results were compared for consistency with nuclease-protection assays and single-molecule force spectroscopy measurements of the SARS-CoV-1 pseudoknot, to determine the most likely conformations. We found several possible conformations for the SARS-CoV-2 pseudoknot, all having an extended stem 3 but with different packing of stems 1 and 2. Several conformations featured rarely-seen threading of a single strand through junctions formed between two helices. These structural models may help interpret future experiments and support efforts to discover ligands inhibiting −1 PRF in SARS-CoV-2.     

An atypical RNA pseudoknot stimulator and an upstream attenuation signal for -1 ribosomal frameshifting of SARS coronavirus

The −1 ribosomal frameshifting requires the existence of an in cis RNA slippery sequence and is promoted by a downstream stimulator RNA. An atypical RNA pseudoknot with an extra stem formed by complementary sequences within loop 2 of an H-type pseudoknot is characterized in the severe acute respiratory syndrome coronavirus (SARS CoV) genome. This pseudoknot can serve as an efficient stimulator for −1 frameshifting in vitro. Mutational analysis of the extra stem suggests frameshift efficiency can be modulated via manipulation of the secondary structure within the loop 2 of an infectious bronchitis virus-type pseudoknot. More importantly, an upstream RNA sequence separated by a linker 5′ to the slippery site is also identified to be capable of modulating the −1 frameshift efficiency. RNA sequence containing this attenuation element can downregulate −1 frameshifting promoted by an atypical pseudoknot of SARS CoV and two other pseudoknot stimulators. Furthermore, frameshift efficiency can be reduced to half in the presence of the attenuation signal in vivo. Therefore, this in cis RNA attenuator represents a novel negative determinant of general importance for the regulation of −1 frameshift efficiency, and is thus a potential antiviral target.     

Rational design of a synthetic mammalian riboswitch as a ligand-responsive -1 ribosomal frame-shifting stimulator

Metabolite-responsive RNA pseudoknots derived from prokaryotic riboswitches have been shown to stimulate −1 programmed ribosomal frameshifting (PRF), suggesting −1 PRF as a promising gene expression platform to extend riboswitch applications in higher eukaryotes. However, its general application has been hampered by difficulty in identifying a specific ligand-responsive pseudoknot that also functions as a ligand-dependent -1 PRF stimulator. We addressed this problem by using the −1 PRF stimulation pseudoknot of SARS-CoV (SARS-PK) to build a ligand-dependent −1 PRF stimulator. In particular, the extra stem of SARS-PK was replaced by an RNA aptamer of theophylline and designed to couple theophylline binding with the stimulation of −1 PRF. Conformational and functional analyses indicate that the engineered theophylline-responsive RNA functions as a mammalian riboswitch with robust theophylline-dependent −1 PRF stimulation activity in a stable human 293T cell-line. Thus, RNA–ligand interaction repertoire provided by in vitro selection becomes accessible to ligand-specific −1 PRF stimulator engineering using SARS-PK as the scaffold for synthetic biology application.     

An intermolecular RNA triplex provides insight into structural determinants for the pseudoknot stimulator of ?1 ribosomal frameshifting

An efficient −1 programmed ribosomal frameshifting (PRF) signal requires an RNA slippery sequence and a downstream RNA stimulator, and the hairpin-type pseudoknot is the most common stimulator. However, a pseudoknot is not sufficient to promote −1 PRF. hTPK-DU177, a pseudoknot derived from human telomerase RNA, shares structural similarities with several −1 PRF pseudoknots and is used to dissect the roles of distinct structural features in the stimulator of −1 PRF. Structure-based mutagenesis on hTPK-DU177 reveals that the −1 PRF efficiency of this stimulator can be modulated by sequential removal of base–triple interactions surrounding the helical junction. Further analysis of the junction-flanking base triples indicates that specific stem–loop interactions and their relative positions to the helical junction play crucial roles for the −1 PRF activity of this pseudoknot. Intriguingly, a bimolecular pseudoknot approach based on hTPK-DU177 reveals that continuing triplex structure spanning the helical junction, lacking one of the loop-closure features embedded in pseudoknot topology, can stimulate −1 PRF. Therefore, the triplex structure is an essential determinant for the DU177 pseudoknot to stimulate −1 PRF. Furthermore, it suggests that −1 PRF, induced by an in-trans RNA via specific base–triple interactions with messenger RNAs, can be a plausible regulatory function for non-coding RNAs.     

De novo 3D models of SARS-CoV-2 RNA elements from consensus experimental secondary structures

The rapid spread of COVID-19 is motivating development of antivirals targeting conserved SARS-CoV-2 molecular machinery. The SARS-CoV-2 genome includes conserved RNA elements that offer potential small-molecule drug targets, but most of their 3D structures have not been experimentally characterized. Here, we provide a compilation of chemical mapping data from our and other labs, secondary structure models, and 3D model ensembles based on Rosetta''s FARFAR2 algorithm for SARS-CoV-2 RNA regions including the individual stems SL1-8 in the extended 5′ UTR; the reverse complement of the 5′ UTR SL1-4; the frameshift stimulating element (FSE); and the extended pseudoknot, hypervariable region, and s2m of the 3′ UTR. For eleven of these elements (the stems in SL1–8, reverse complement of SL1–4, FSE, s2m and 3′ UTR pseudoknot), modeling convergence supports the accuracy of predicted low energy states; subsequent cryo-EM characterization of the FSE confirms modeling accuracy. To aid efforts to discover small molecule RNA binders guided by computational models, we provide a second set of similarly prepared models for RNA riboswitches that bind small molecules. Both datasets (‘FARFAR2-SARS-CoV-2’, https://github.com/DasLab/FARFAR2-SARS-CoV-2; and ‘FARFAR2-Apo-Riboswitch’, at https://github.com/DasLab/FARFAR2-Apo-Riboswitch’) include up to 400 models for each RNA element, which may facilitate drug discovery approaches targeting dynamic ensembles of RNA molecules.     

Structure-altering mutations of the SARS-CoV-2 frameshifting RNA element

With the rapid rate of COVID-19 infections and deaths, treatments and cures besides hand washing, social distancing, masks, isolation, and quarantines are urgently needed. The treatments and vaccines rely on the basic biophysics of the complex viral apparatus. Although proteins are serving as main drug and vaccine targets, therapeutic approaches targeting the 30,000 nucleotide RNA viral genome form important complementary approaches. Indeed, the high conservation of the viral genome, its close evolutionary relationship to other viruses, and the rise of gene editing and RNA-based vaccines all argue for a focus on the RNA agent itself. One of the key steps in the viral replication cycle inside host cells is the ribosomal frameshifting required for translation of overlapping open reading frames. The RNA frameshifting element (FSE), one of three highly conserved regions of coronaviruses, is believed to include a pseudoknot considered essential for this ribosomal switching. In this work, we apply our graph-theory-based framework for representing RNA secondary structures, “RAG (or RNA-As-Graphs),” to alter key structural features of the FSE of the SARS-CoV-2 virus. Specifically, using RAG machinery of genetic algorithms for inverse folding adapted for RNA structures with pseudoknots, we computationally predict minimal mutations that destroy a structurally important stem and/or the pseudoknot of the FSE, potentially dismantling the virus against translation of the polyproteins. Our microsecond molecular dynamics simulations of mutant structures indicate relatively stable secondary structures. These findings not only advance our computational design of RNAs containing pseudoknots, they pinpoint key residues of the SARS-CoV-2 virus as targets for antiviral drugs and gene editing approaches.     

Highly conserved RNA pseudoknots at the gag-pol junction of HIV-1 suggest a novel mechanism of ?1 ribosomal frameshifting

−1 programmed ribosomal frameshifting (PRF) is utilized by many viruses to synthesize their enzymatic (Pol) and structural (Gag) proteins at a defined ratio. For efficient −1 PRF, two cis-acting elements are required: a heptanucleotide frameshift site and a downstream stimulator such as a pseudoknot. We have analyzed the gag-pol junction sequences from 4254 HIV-1 strains. Approximately ninety-five percent of the sequences can form four pseudoknots PK1–PK4 (∼97% contain PK1, PK3, and PK4), covering ∼72 nt including the frameshift site. Some pseudoknots are mutually excluded due to sequence overlap. PK1 and PK3 arrange tandemly. Their stems form a quasi-continuous helix of ∼22 bp. We propose a novel mechanism for possible roles of these pseudoknots. Multiple alternative structures may exist at the gag-pol junction. In most strains, the PK1–PK3 tandem pseudoknots may dominate the structurally heterogeneous pool of RNA due to their greater overall stability. The tandem pseudoknots may function as a breaking system to slow down the ribosome. The ribosome unwinds PK1 and stem 1 of PK3 before it can reach the frameshift site. Then, PK4 can form rapidly because the intact stem 2 of PK3 makes up a large part of the stem 1 of PK4. The newly formed PK4 jams the entrance of the mRNA tunnel. The process then proceeds as in a typical case of −1 PRF. This mechanism incorporates several exquisite new features while still being consistent with the current paradigm of pseudoknot-dependent −1 PRF.     

Single-Molecule Measurements of the CCR5 mRNA Unfolding Pathways

Secondary or tertiary structure in an mRNA, such as a pseudoknot, can create a physical barrier that requires the ribosome to generate additional force to translocate. The presence of such a barrier can dramatically increase the probability that the ribosome will shift into an alternate reading frame, in which a different set of codons is recognized. The detailed biophysical mechanism by which frameshifting is induced remains unknown. Here we employ optical trapping techniques to investigate the structure of a −1 programmed ribosomal frameshift (−1 PRF) sequence element located in the CCR5 mRNA, which encodes a coreceptor for HIV-1 and is, to our knowledge, the first known human −1 PRF signal of nonviral origin. We begin by presenting a set of computationally predicted structures that include pseudoknots. We then employ what we believe to be new analytical techniques for measuring the effective free energy landscapes of biomolecules. We find that the −1 PRF element manifests several distinct unfolding pathways when subject to end-to-end force, one of which is consistent with a proposed pseudoknot conformation, and another of which we have identified as a folding intermediate. The dynamic ensemble of conformations that CCR5 mRNA exhibits in the single-molecule experiments may be a significant feature of the frameshifting mechanism.     

A −1 frameshift in the HIV-1 env gene is enhanced by arginine deficiency via a hungry codon mechanism

Ribosomal frameshifting is used by various organisms to maximize protein coding potential of genomic sequences. It is commonly exploited by RNA viruses to overcome the constraint of their limited genome size. Frameshifting requires specific RNA structural features, such as a suitable heptanucleotide “slippery” sequence and an RNA pseudoknot. Previous genomic analysis of HIV-1 indicated the potential for several hidden genes encoded through frameshifting; one of these, overlapping the envelope gene, has an RNA pseudoknot just downstream from a slippery sequence, AAAAAGA that features an adenine quadruplet prior to a potential hungry arginine codon (AGA). This env-frameshift (env-fs) gene has been shown to encode a truncated glutathione peroxidase homologue, with both antioxidant and anti-apoptotic activities in transfected cells. Using a dual reporter cell-based frameshift assay, we demonstrate that the env-fs frameshift sequence is active in vitro. Furthermore, in arginine deficient media, env-fs frameshifting increased over 100% (p < 0.005), consistent with the hypothesized hungry codon mechanism. As a response to arginine deficiency, increased expression of the antioxidant viral GPx gene (env-fs) by upregulation of frameshifting could be protective to HIV-infected cells, as a countermeasure to the increased oxidative stress induced by arginine deficiency (because NO is a known scavenger of hydroxyl radical).     

A three-stemmed mRNA pseudoknot in the SARS coronavirus frameshift signal

A wide range of RNA viruses use programmed −1 ribosomal frameshifting for the production of viral fusion proteins. Inspection of the overlap regions between ORF1a and ORF1b of the SARS-CoV genome revealed that, similar to all coronaviruses, a programmed −1 ribosomal frameshift could be used by the virus to produce a fusion protein. Computational analyses of the frameshift signal predicted the presence of an mRNA pseudoknot containing three double-stranded RNA stem structures rather than two. Phylogenetic analyses showed the conservation of potential three-stemmed pseudoknots in the frameshift signals of all other coronaviruses in the GenBank database. Though the presence of the three-stemmed structure is supported by nuclease mapping and two-dimensional nuclear magnetic resonance studies, our findings suggest that interactions between the stem structures may result in local distortions in the A-form RNA. These distortions are particularly evident in the vicinity of predicted A-bulges in stems 2 and 3. In vitro and in vivo frameshifting assays showed that the SARS-CoV frameshift signal is functionally similar to other viral frameshift signals: it promotes efficient frameshifting in all of the standard assay systems, and it is sensitive to a drug and a genetic mutation that are known to affect frameshifting efficiency of a yeast virus. Mutagenesis studies reveal that both the specific sequences and structures of stems 2 and 3 are important for efficient frameshifting. We have identified a new RNA structural motif that is capable of promoting efficient programmed ribosomal frameshifting. The high degree of conservation of three-stemmed mRNA pseudoknot structures among the coronaviruses suggests that this presents a novel target for antiviral therapeutics.     

Secondary structure determination of conserved SARS-CoV-2 RNA elements by NMR spectroscopy

The current pandemic situation caused by the Betacoronavirus SARS-CoV-2 (SCoV2) highlights the need for coordinated research to combat COVID-19. A particularly important aspect is the development of medication. In addition to viral proteins, structured RNA elements represent a potent alternative as drug targets. The search for drugs that target RNA requires their high-resolution structural characterization. Using nuclear magnetic resonance (NMR) spectroscopy, a worldwide consortium of NMR researchers aims to characterize potential RNA drug targets of SCoV2. Here, we report the characterization of 15 conserved RNA elements located at the 5′ end, the ribosomal frameshift segment and the 3′-untranslated region (3′-UTR) of the SCoV2 genome, their large-scale production and NMR-based secondary structure determination. The NMR data are corroborated with secondary structure probing by DMS footprinting experiments. The close agreement of NMR secondary structure determination of isolated RNA elements with DMS footprinting and NMR performed on larger RNA regions shows that the secondary structure elements fold independently. The NMR data reported here provide the basis for NMR investigations of RNA function, RNA interactions with viral and host proteins and screening campaigns to identify potential RNA binders for pharmaceutical intervention.     

FSscan: a mechanism-based program to identify +1 ribosomal frameshift hotspots

In +1 programmed ribosomal frameshifting (PRF), ribosomes skip one nucleotide toward the 3′-end during translation. Most of the genes known to demonstrate +1 PRF have been discovered by chance or by searching homologous genes. Here, a bioinformatic framework called FSscan is developed to perform a systematic search for potential +1 frameshift sites in the Escherichia coli genome. Based on a current state of the art understanding of the mechanism of +1 PRF, FSscan calculates scores for a 16-nt window along a gene sequence according to different effects of the stimulatory signals, and ribosome E-, P- and A-site interactions. FSscan successfully identified the +1 PRF site in prfB and predicted yehP, pepP, nuoE and cheA as +1 frameshift candidates in the E. coli genome. Empirical results demonstrated that potential +1 frameshift sequences identified promoted significant levels of +1 frameshifting in vivo. Mass spectrometry analysis confirmed the presence of the frameshifted proteins expressed from a yehP-egfp fusion construct. FSscan allows a genome-wide and systematic search for +1 frameshift sites in E. coli. The results have implications for bioinformatic identification of novel frameshift proteins, ribosomal frameshifting, coding sequence detection and the application of mass spectrometry on studying frameshift proteins.     

Possible involvement of coaxially stacked double pseudoknots in the regulation of −1 programmed ribosomal frameshifting in RNA viruses

?1 programmed ribosomal frameshifting (PRF) in viruses is often stimulated by a pseudoknot downstream from the slippery sequence. At the PRF junction of HIV-1, transmissible gastroenteritis virus (TGEV), Barmah Forest virus (BFV), Fort Morgan virus (FMV), and Equine arteritis virus (EAV), we identified potential double pseudoknots in either a tandem mode or embedded mode. In viruses with tandem pseudoknots (5′PK & 3′PK), the slippery sequence is encompassed in the 5′PK. The ribosome needs to unwind the 5′PK to get to the slippery sequence. In HIV-1, the 3′PK and several alternative structures are mutually exclusive. Disruption of the tandem pseudoknots may enable one of the alternative structures to form as the effective frameshift stimulator. In TGEV/BFV/FMV, the 3′PK is a conventional frameshift stimulator. In all cases, the tandem pseudoknots may slow down the ribosome before it reaches the conventional PRF signals. In EAV, a compact pseudoknot is embedded within loop2 of the otherwise conventional frameshift-stimulating pseudoknot. All double pseudoknots have the potential to stack their stems coaxially. We built structural models of the HIV-1 and EAV double pseudoknots to show that both the tandem and embedded modes are feasible and reasonable. We hypothesize that the fundamental reason for the viruses to utilize coaxially stacked double pseudoknots is to increase the overall stability of the frameshift regulating structure, and avoid an ultra-stable single pseudoknot which may become a ribosomal roadblock. Our results significantly expand the repertoire of RNA structures and dynamics that may potentially involve in ?1 PRF regulation.     

Three-dimensional structure of a flavivirus dumbbell RNA reveals molecular details of an RNA regulator of replication

Mosquito-borne flaviviruses (MBFVs) including dengue, West Nile, yellow fever, and Zika viruses have an RNA genome encoding one open reading frame flanked by 5′ and 3′ untranslated regions (UTRs). The 3′ UTRs of MBFVs contain regions of high sequence conservation in structured RNA elements known as dumbbells (DBs). DBs regulate translation and replication of the viral RNA genome, functions proposed to depend on the formation of an RNA pseudoknot. To understand how DB structure provides this function, we solved the x-ray crystal structure of the Donggang virus DB to 2.1Å resolution and used structural modeling to reveal the details of its three-dimensional fold. The structure confirmed the predicted pseudoknot and molecular modeling revealed how conserved sequences form a four-way junction that appears to stabilize the pseudoknot. Single-molecule FRET suggests that the DB pseudoknot is a stable element that can regulate the switch between translation and replication during the viral lifecycle by modulating long-range RNA conformational changes.     

A Nascent Peptide Signal Responsive to Endogenous Levels of Polyamines Acts to Stimulate Regulatory Frameshifting on Antizyme mRNA

The protein antizyme is a negative regulator of cellular polyamine concentrations from yeast to mammals. Synthesis of functional antizyme requires programmed +1 ribosomal frameshifting at the 3′ end of the first of two partially overlapping ORFs. The frameshift is the sensor and effector in an autoregulatory circuit. Except for Saccharomyces cerevisiae antizyme mRNA, the frameshift site alone only supports low levels of frameshifting. The high levels usually observed depend on the presence of cis-acting stimulatory elements located 5′ and 3′ of the frameshift site. Antizyme genes from different evolutionary branches have evolved different stimulatory elements. Prior and new multiple alignments of fungal antizyme mRNA sequences from the Agaricomycetes class of Basidiomycota show a distinct pattern of conservation 5′ of the frameshift site consistent with a function at the amino acid level. As shown here when tested in Schizosaccharomyces pombe and mammalian HEK293T cells, the 5′ part of this conserved sequence acts at the nascent peptide level to stimulate the frameshifting, without involving stalling detectable by toe-printing. However, the peptide is only part of the signal. The 3′ part of the stimulator functions largely independently and acts at least mostly at the nucleotide level. When polyamine levels were varied, the stimulatory effect was seen to be especially responsive in the endogenous polyamine concentration range, and this effect may be more general. A conserved RNA secondary structure 3′ of the frameshift site has weaker stimulatory and polyamine sensitizing effects on frameshifting.     

Ribosomal Protein L3 Mutants Alter Translational Fidelity and Promote Rapid Loss of the Yeast Killer Virus

Programmed −1 ribosomal frameshifting is utilized by a number of RNA viruses as a means of ensuring the correct ratio of viral structural to enzymatic proteins available for viral particle assembly. Altering frameshifting efficiencies upsets this ratio, interfering with virus propagation. We have previously demonstrated that compounds that alter the kinetics of the peptidyl-transfer reaction affect programmed −1 ribosomal frameshift efficiencies and interfere with viral propagation in yeast. Here, the use of a genetic approach lends further support to the hypothesis that alterations affecting the ribosome’s peptidyltransferase activity lead to changes in frameshifting efficiency and virus loss. Mutations in the RPL3 gene, which encodes a ribosomal protein located at the peptidyltransferase center, promote approximately three- to fourfold increases in programmed −1 ribosomal frameshift efficiencies and loss of the M₁ killer virus of yeast. The mak8-1 allele of RPL3 contains two adjacent missense mutations which are predicted to structurally alter the Mak8-1p. Furthermore, a second allele that encodes the N-terminal 100 amino acids of L3 (called L3Δ) exerts a trans-dominant effect on programmed −1 ribosomal frameshifting and killer virus maintenance. Taken together, these results support the hypothesis that alterations in the peptidyltransferase center affect programmed −1 ribosomal frameshifting.     

KnotInFrame: prediction of -1 ribosomal frameshift events

Programmed -1 ribosomal frameshift (-1 PRF) allows for alternative reading frames within one mRNA. First found in several viruses, it is now believed to exist in all kingdoms of life. Strong stimulators for -1 PRF are a heptameric slippery site and an RNA pseudoknot. Here, we present a new algorithm KnotInFrame, for the automatic detection of -1 PRF signals from genomic sequences. It finds the frameshifting stimulators by means of a specialized RNA-pseudoknot folding program, fast enough for genome-wide analyses. Evaluations on known -1 PRF signals demonstrate a high sensitivity.     

Regulation of Programmed Ribosomal Frameshifting by Co-Translational Refolding RNA Hairpins

RNA structures are unwound for decoding. In the process, they can pause the elongating ribosome for regulation. An example is the stimulation of -1 programmed ribosomal frameshifting, leading to 3′ direction slippage of the reading-frame during elongation, by specific pseudoknot stimulators downstream of the frameshifting site. By investigating a recently identified regulatory element upstream of the SARS coronavirus (SARS-CoV) −1 frameshifting site, it is shown that a minimal functional element with hairpin forming potential is sufficient to down-regulate−1 frameshifting activity. Mutagenesis to disrupt or restore base pairs in the potential hairpin stem reveals that base-pair formation is required for−1 frameshifting attenuation in vitro and in 293T cells. The attenuation efficiency of a hairpin is determined by its stability and proximity to the frameshifting site; however, it is insensitive to E site sequence variation. Additionally, using a dual luciferase assay, it can be shown that a hairpin stimulated +1 frameshifting when placed upstream of a +1 shifty site in yeast. The investigations indicate that the hairpin is indeed a cis-acting programmed reading-frame switch modulator. This result provides insight into mechanisms governing−1 frameshifting stimulation and attenuation. Since the upstream hairpin is unwound (by a marching ribosome) before the downstream stimulator, this study’s findings suggest a new mode of translational regulation that is mediated by the reformed stem of a ribosomal unwound RNA hairpin during elongation.