Crystal structure of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) frameshifting pseudoknot

A Cascade of Conformational Switches in SARS-CoV-2 Frameshifting: Coregulation by Upstream and Downstream Elements

Targeting ribosomal frameshifting has emerged as a potential therapeutic intervention strategy against COVID-19. In this process, a -1 shift in the ribosomal reading frame encodes alternative viral proteins. Any interference with this process profoundly affects viral replication and propagation. For SARS-CoV-2, two RNA sites associated with ribosomal frameshifting are positioned on the 5' and 3' of the frameshifting residues. Although much attention has been focused on the 3' frameshift element (FSE), the 5' stem-loop (attenuator hairpin, AH) can play a role. Yet the relationship between the two regions is unknown. In addition, multiple folds of the FSE and FSE-containing RNA regions have been discovered. To gain more insight into these RNA folds in the larger sequence context that includes AH, we apply our graph-theory-based modeling tools to represent RNA secondary structures, "RAG" (RNA-As-Graphs), to generate conformational landscapes that suggest length-dependent conformational distributions. We show that the AH region can coexist as a stem-loop with main and alternative 3-stem pseudoknots of the FSE (dual graphs 3_6 and 3_3 in our notation) but that an alternative stem 1 (AS1) can disrupt the FSE pseudoknots and trigger other folds. A critical length for AS1 of 10-bp regulates key folding transitions. Together with designed mutants and available experimental data, we present a sequential view of length-dependent folds during frameshifting and suggest their mechanistic roles. These structural and mutational insights into both ends of the FSE advance our understanding of the SARS-CoV-2 frameshifting mechanism by suggesting how alternative folds play a role in frameshifting and defining potential therapeutic intervention techniques that target specific folds.

RNA-Puzzles Round V: blind predictions of 23 RNA structures

RNA-Puzzles is a collective endeavor dedicated to the advancement and improvement of RNA three-dimensional structure prediction. With agreement from structural biologists, RNA structures are predicted by modeling groups before publication of the experimental structures. We report a large-scale set of predictions by 18 groups for 23 RNA-Puzzles: 4 RNA elements, 2 Aptamers, 4 Viral elements, 5 Ribozymes and 8 Riboswitches. We describe automatic assessment protocols for comparisons between prediction and experiment. Our analyses reveal some critical steps to be overcome to achieve good accuracy in modeling RNA structures: identification of helix-forming pairs and of non-Watson-Crick modules, correct coaxial stacking between helices and avoidance of entanglements. Three of the top four modeling groups in this round also ranked among the top four in the CASP15 contest.

Elaborated pseudoknots that stimulate -1 programmed ribosomal frameshifting or stop codon readthrough in RNA viruses

Pseudoknots assume various functions including stimulation of -1 programmed ribosomal frameshifting (PRF) or stop codon readthrough (SCR) in RNA viruses. These pseudoknots vary greatly in sizes and structural complexities. Recent biochemical and structural studies confirm the three-stemmed pseudoknots as the -1 PRF stimulators in the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and related coronaviruses. We reexamined previously reported -1 PRF or SCR stimulating pseudoknots, especially those containing a relatively long connecting loop between the two pseudoknot-forming stems, for their ability to form elaborated structures. Many potential elaborated pseudoknots were identified that contain one or more of the following extra structural elements: stem-loop, embedded pseudoknot, kissing hairpins, and additional loop-loop interactions. The elaborated pseudoknots are found in several different virus families that utilize either the -1 PRF or SCR recoding mechanisms. Model-building studies were performed to not only establish the structural feasibility of the elaborated pseudoknots but also reveal potential additional structural features that cannot be readily inferred from the predicted secondary structures. Some of the structures, such as embedded double pseudoknots and compact loop-loop pseudoknots mediated by the previously established common pseudoknot motif-1 (CPK-1), represent the first of its kind in the literatures. By advancing discovery of new functional RNA structures, we significantly expand the repertoire of known elaborated pseudoknots that could potentially play a role in -1 PRF and SCR regulation. These results contribute to a better understanding of RNA structures in general, facilitating the design of engineering RNA molecules with certain desired functions.Communicated by Ramaswamy H. Sarma.

Optimization of Structure-Guided Development of Chemical Probes for the Pseudoknot RNA of the Frameshift Element in SARS-CoV-2

Targeting the RNA genome of SARS-CoV-2 is a viable option for antiviral drug development. We explored three ligand binding sites of the core pseudoknot RNA of the SARS-CoV-2 frameshift element. We iteratively optimized ligands, based on improved affinities, targeting these binding sites and report on structural and dynamic properties of the three identified binding sites. Available experimental 3D structures of the pseudoknot element were compared to SAXS and NMR data to validate its dominant folding state in solution. In order to experimentally map in silico predicted binding sites, NMR assignments of the majority of nucleobases were achieved by segmental labeling of the pseudoknot RNA and isotope-filtered NMR experiments at 1.2 GHz, demonstrating the value of NMR spectroscopy to supplement modelling and docking data. Optimized ligands with enhanced affinity were shown to specifically inhibit frameshifting without affecting 0-frame translation in cell-free translation assays, establishing the frameshift element as target for drug-like ligands of low molecular weight.

Heterogeneous and multiple conformational transition pathways between pseudoknots of the SARS-CoV-2 frameshift element

Frameshifting is an essential mechanism employed by many viruses including coronaviruses to produce viral proteins from a compact RNA genome. It is facilitated by specific RNA folds in the frameshift element (FSE), which has emerged as an important therapeutic target. For SARS-CoV-2, a specific 3-stem pseudoknot has been identified to stimulate frameshifting. However, prior studies and our RNA-As-Graphs analysis coupled to chemical reactivity experiments revealed other folds, including a different pseudoknot. Although structural plasticity has been proposed to play a key role in frameshifting, paths between different FSE RNA folds have not been yet identified. Here, we capture atomic-level transition pathways between two key FSE pseudoknots by transition path sampling coupled to Markov State Modeling and our BOLAS free energy method. We reveal multiple transition paths within a heterogeneous, multihub conformational landscape. A shared folding mechanism involves RNA stem unpairing followed by a 5'-chain end release. Significantly, this pseudoknot transition critically tunes the tension through the RNA spacer region and places the viral RNA in the narrow ribosomal channel. Our work further explains the role of the alternative pseudoknot in ribosomal pausing and clarifies why the experimentally captured pseudoknot is preferred for frameshifting. Our capturing of this large-scale transition of RNA secondary and tertiary structure highlights the complex pathways of biomolecules and the inherent multifarious aspects that viruses developed to ensure virulence and survival. This enhanced understanding of viral frameshifting also provides insights to target key transitions for therapeutic applications. Our methods are generally applicable to other large-scale biomolecular transitions.

CParty: hierarchically constrained partition function of RNA pseudoknots

MOTIVATION

Biologically relevant RNA secondary structures are routinely predicted by efficient dynamic programming algorithms that minimize their free energy. Starting from such algorithms, one can devise partition function algorithms, which enable stochastic perspectives on RNA structure ensembles. As the most prominent example, McCaskill's partition function algorithm is derived from pseudoknot-free energy minimization. While this algorithm became hugely successful for the analysis of pseudoknot-free RNA structure ensembles, as of yet there exists only one pseudoknotted partition function implementation, which covers only simple pseudoknots and comes with a borderline-prohibitive complexity of O(n5) in the RNA length n.

RESULTS

Here, we develop a partition function algorithm corresponding to the hierarchical pseudoknot prediction of HFold, which performs exact optimization in a realistic pseudoknot energy model. In consequence, our algorithm CParty carries over HFold's advantages over classical pseudoknot prediction in characterizing the Boltzmann ensemble at equilibrium. Given an RNA sequence S and a pseudoknot-free structure G, CParty computes the partition function over all possibly pseudoknotted density-2 structures G∪G' of S that extend the fixed G by a disjoint pseudoknot-free structure G'. Thus, CParty follows the common hypothesis of hierarchical pseudoknot formation, where pseudoknots form as tertiary contacts only after a first pseudoknot-free "core" G and we call the computed partition function hierarchically constrained (by G). Like HFold, the dynamic programming algorithm CParty is very efficient, achieving the low complexity of the pseudoknot-free algorithm, i.e. cubic time and quadratic space. Finally, by computing pseudoknotted ensemble energies, we unveil kinetics features of a therapeutic target in SARS-CoV-2.

AVAILABILITY AND IMPLEMENTATION

CParty is available at https://github.com/HosnaJabbari/CParty.

Cyclic peptides targeting the SARS-CoV-2 programmed ribosomal frameshifting RNA from a multiplexed phage display library

RNA provides the genetic blueprint for many pathogenic viruses, including SARS-CoV-2. The propensity of RNA to fold into specific tertiary structures enables the biomolecular recognition of cavities and crevices suited for the binding of drug-like molecules. Despite increasing interest in RNA as a target for chemical biology and therapeutic applications, the development of molecules that recognize RNA with high affinity and specificity represents a significant challenge. Here, we report a strategy for the discovery and selection of RNA-targeted macrocyclic peptides derived from combinatorial libraries of peptide macrocycles displayed by bacteriophages. Specifically, a platform for phage display of macrocyclic organo-peptide hybrids (MOrPH-PhD) was combined with a diverse set of non-canonical amino acid-based cyclization modules to produce large libraries of 107 structurally diverse, genetically encoded peptide macrocycles. These libraries were panned against the -1 programmed ribosomal frameshifting stimulatory sequence (FSS) RNA pseudoknot of SARS-CoV-2, which revealed specific macrocyclic peptide sequences that bind this essential motif with high affinity and selectivity. Peptide binding localizes to the FSS dimerization loop based on chemical modification analysis and binding assays and the cyclic peptides show specificity toward the target RNA over unrelated RNA pseudoknots. This work introduces a novel system for the generation and high-throughput screening of topologically diverse cyclopeptide scaffolds (multiplexed MOrPH-PhD), and it provides a blueprint for the exploration and evolution of genetically encoded macrocyclic peptides that target specific RNAs.

Dissecting the Conformational Heterogeneity of Stem-Loop Substructures of the Fifth Element in the 5'-Untranslated Region of SARS-CoV-2

Throughout the family of coronaviruses, structured RNA elements within the 5' region of the genome are highly conserved. The fifth stem-loop element from SARS-CoV-2 (5_SL5) represents an example of an RNA structural element, repeatedly occurring in coronaviruses. It contains a conserved, repetitive fold within its substructures SL5a and SL5b. We herein report the detailed characterization of the structure and dynamics of elements SL5a and SL5b that are located immediately upstream of the SARS-CoV-2 ORF1a/b start codon. Exploiting the unique ability of solution NMR methods, we show that the structures of both apical loops are modulated by structural differences in the remote parts located in their stem regions. We further integrated our high-resolution models of SL5a/b into the context of full-length 5_SL5 structures by combining different structural biology methods. Finally, we evaluated the impact of the two most common VoC mutations within 5_SL5 with respect to individual base-pair stability.

Sarbecovirus programmed ribosome frameshift RNA element folding studied by NMR spectroscopy and comparative analyses

The programmed ribosomal frameshift (PRF) region is found in the RNA genome of all coronaviruses and shifts the ribosome reading frame through formation of a three-stem pseudoknot structure, allowing the translation of essential viral proteins. Using NMR spectroscopy, comparative sequence analyses and functional assays we show that, in the absence of the ribosome, a 123-nucleotide sequence encompassing the PRF element of SARS-CoV-2 adopts a well-defined two-stem loop structure that is conserved in all SARS-like coronaviruses. In this conformation, the attenuator hairpin and slippery site nucleotides are exposed in the first stem-loop and two pseudoknot stems are present in the second stem-loop, separated by an 8-nucleotide bulge. Formation of the third pseudoknot stem depends on pairing between bulge nucleotides and base-paired nucleotides of the upstream stem-loop, as shown by a PRF construct where residues of the upstream stem were removed, which formed the pseudoknot structure and had increased frameshifting activity in a dual-luciferase assay. The base-pair switch driving PRF pseudoknot folding was found to be conserved in several human non-SARS coronaviruses. The collective results suggest that the frameshifting pseudoknot structure of these viruses only forms transiently in the presence of the translating ribosome. These findings clarify the frameshifting mechanism in coronaviruses and can have a beneficial impact on antiviral drug discovery.

Abolished frameshifting for predicted structure-stabilizing SARS-CoV-2 mutants: implications to alternative conformations and their statistical structural analyses

The SARS-CoV-2 frameshifting element (FSE) has been intensely studied and explored as a therapeutic target for coronavirus diseases, including COVID-19. Besides the intriguing virology, this small RNA is known to adopt many length-dependent conformations, as verified by multiple experimental and computational approaches. However, the role these alternative conformations play in the frameshifting mechanism and how to quantify this structural abundance has been an ongoing challenge. Here, we show by DMS and dual-luciferase functional assays that previously predicted FSE mutants (using the RAG graph theory approach) suppress structural transitions and abolish frameshifting. Furthermore, correlated mutation analysis of DMS data by three programs (DREEM, DRACO, and DANCE-MaP) reveals important differences in their estimation of specific RNA conformations, suggesting caution in the interpretation of such complex conformational landscapes. Overall, the abolished frameshifting in three different mutants confirms that all alternative conformations play a role in the pathways of ribosomal transition.

Novel Inhibitors of SARS-CoV-2 RNA Identified through Virtual Screening

We currently lack antivirals for most human viruses. In a quest for new molecules, focusing on viral RNA, instead of viral proteins, can represent a promising strategy. In this study, new inhibitors were identified starting from a published crystal structure of the tertiary SARS-CoV-2 RNA involved in the -1 programmed ribosomal frameshift. The pseudoknot structure was refined, and a virtual screening was performed using the repository of binders to the nucleic acid library, taking into consideration RNA flexibility. Hit compounds were validated against the wild-type virus and with a dual-luciferase assay measuring the frameshift efficiency. Several active molecules were identified. Our study reveals new inhibitors of SARS-CoV-2 but also highlights the feasibility of targeting RNA starting from virtual screening, a strategy that could be broadly applied to drug development.

Discovery and Quantification of Long-Range RNA Base Pairs in Coronavirus Genomes with SEARCH-MaP and SEISMIC-RNA

RNA molecules perform a diversity of essential functions for which their linear sequences must fold into higher-order structures. Techniques including crystallography and cryogenic electron microscopy have revealed 3D structures of ribosomal, transfer, and other well-structured RNAs; while chemical probing with sequencing facilitates secondary structure modeling of any RNAs of interest, even within cells. Ongoing efforts continue increasing the accuracy, resolution, and ability to distinguish coexisting alternative structures. However, no method can discover and quantify alternative structures with base pairs spanning arbitrarily long distances - an obstacle for studying viral, messenger, and long noncoding RNAs, which may form long-range base pairs. Here, we introduce the method of Structure Ensemble Ablation by Reverse Complement Hybridization with Mutational Profiling (SEARCH-MaP) and software for Structure Ensemble Inference by Sequencing, Mutation Identification, and Clustering of RNA (SEISMIC-RNA). We use SEARCH-MaP and SEISMIC-RNA to discover that the frameshift stimulating element of SARS coronavirus 2 base-pairs with another element 1 kilobase downstream in nearly half of RNA molecules, and that this structure competes with a pseudoknot that stimulates ribosomal frameshifting. Moreover, we identify long-range base pairs involving the frameshift stimulating element in other coronaviruses including SARS coronavirus 1 and transmissible gastroenteritis virus, and model the full genomic secondary structure of the latter. These findings suggest that long-range base pairs are common in coronaviruses and may regulate ribosomal frameshifting, which is essential for viral RNA synthesis. We anticipate that SEARCH-MaP will enable solving many RNA structure ensembles that have eluded characterization, thereby enhancing our general understanding of RNA structures and their functions. SEISMIC-RNA, software for analyzing mutational profiling data at any scale, could power future studies on RNA structure and is available on GitHub and the Python Package Index.

Discovery and Quantification of Long-Range RNA Base Pairs in Coronavirus Genomes with SEARCH-MaP and SEISMIC-RNA

RNA molecules perform a diversity of essential functions for which their linear sequences must fold into higher-order structures. Techniques including crystallography and cryogenic electron microscopy have revealed 3D structures of ribosomal, transfer, and other well-structured RNAs; while chemical probing with sequencing facilitates secondary structure modeling of any RNAs of interest, even within cells. Ongoing efforts continue increasing the accuracy, resolution, and ability to distinguish coexisting alternative structures. However, no method can discover and quantify alternative structures with base pairs spanning arbitrarily long distances - an obstacle for studying viral, messenger, and long noncoding RNAs, which may form long-range base pairs. Here, we introduce the method of Structure Ensemble Ablation by Reverse Complement Hybridization with Mutational Profiling (SEARCH-MaP) and software for Structure Ensemble Inference by Sequencing, Mutation Identification, and Clustering of RNA (SEISMIC-RNA). We use SEARCH-MaP and SEISMIC-RNA to discover that the frameshift stimulating element of SARS coronavirus 2 base-pairs with another element 1 kilobase downstream in nearly half of RNA molecules, and that this structure competes with a pseudoknot that stimulates ribosomal frameshifting. Moreover, we identify long-range base pairs involving the frameshift stimulating element in other coronaviruses including SARS coronavirus 1 and transmissible gastroenteritis virus, and model the full genomic secondary structure of the latter. These findings suggest that long-range base pairs are common in coronaviruses and may regulate ribosomal frameshifting, which is essential for viral RNA synthesis. We anticipate that SEARCH-MaP will enable solving many RNA structure ensembles that have eluded characterization, thereby enhancing our general understanding of RNA structures and their functions. SEISMIC-RNA, software for analyzing mutational profiling data at any scale, could power future studies on RNA structure and is available on GitHub and the Python Package Index.

Possible involvement of three-stemmed pseudoknots in regulating translational initiation in human mRNAs

RNA pseudoknots play a crucial role in various cellular functions. Established pseudoknots show significant variation in both size and structural complexity. Specifically, three-stemmed pseudoknots are characterized by an additional stem-loop embedded in their structure. Recent findings highlight these pseudoknots as bacterial riboswitches and potent stimulators for programmed ribosomal frameshifting in RNA viruses like SARS-CoV2. To investigate the possible presence of functional three-stemmed pseudoknots in human mRNAs, we employed in-house developed computational methods to detect such structures within a dataset comprising 21,780 full-length human mRNA sequences. Numerous three-stemmed pseudoknots were identified. A selected set of 14 potential instances are presented, in which the start codon of the mRNA is found in close proximity either upstream, downstream, or within the identified three-stemmed pseudoknot. These pseudoknots likely play a role in translational initiation regulation. The probability of their existence gains support from their ranking as the most stable pseudoknot identified in the entire mRNA sequence, structural conservation across homologous mRNAs, stereochemical feasibility as demonstrated by structural modeling, and classification as members of the CPK-1 pseudoknot family, which includes many well-established pseudoknots. Furthermore, in four of the mRNAs, two or three closely spaced or tandem three-stemmed pseudoknots were identified. These findings suggest the frequent occurrence of three-stemmed pseudoknots in human mRNAs. A stepwise co-transcriptional folding mechanism is proposed for the formation of a three-stemmed pseudoknot structure. Our results not only provide fresh insights into the structures and functions of pseudoknots but also unveil the potential to target pseudoknots for treating human diseases.

Structure of the SARS-CoV-2 Frameshift Stimulatory Element with an Upstream Multibranch Loop

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) frameshift stimulatory element (FSE) is necessary for programmed -1 ribosomal frameshifting (-1 PRF) and optimized viral efficacy. The FSE has an abundance of context-dependent alternate conformations, but two of the structures most crucial to -1 PRF are an attenuator hairpin and a three-stem H-type pseudoknot structure. A crystal structure of the pseudoknot alone features three RNA stems in a helically stacked linear structure, whereas a 6.9 Å cryo-EM structure including the upstream heptameric slippery site resulted in a bend between two stems. Our previous research alluded to an extended upstream multibranch loop that includes both the attenuator hairpin and the slippery site-a conformation not previously modeled. We aim to provide further context to the SARS-CoV-2 FSE via computational and medium resolution cryo-EM approaches, by presenting a 6.1 Å cryo-EM structure featuring a linear pseudoknot structure and a dynamic upstream multibranch loop.

Abolished frameshifting for predicted structure-stabilizing SARS-CoV-2 mutants: Implications to alternative conformations and their statistical structural analyses

The SARS-CoV-2 frameshifting element (FSE) has been intensely studied and explored as a therapeutic target for coronavirus diseases including COVID-19. Besides the intriguing virology, this small RNA is known to adopt many length-dependent conformations, as verified by multiple experimental and computational approaches. However, the role these alternative conformations play in the frameshifting mechanism and how to quantify this structural abundance has been an ongoing challenge. Here, we show by DMS and dual-luciferase functional assays that previously predicted FSE mutants (using the RAG graph theory approach) suppress structural transitions and abolish frameshifting. Furthermore, correlated mutation analysis of DMS data by three programs (DREEM, DRACO, and DANCE-MaP) reveals important differences in their estimation of specific RNA conformations, suggesting caution in the interpretation of such complex conformational landscapes. Overall, the abolished frameshifting in three different mutants confirms that all alternative conformations play a role in the pathways of ribosomal transition.

Tertiary folds of the SL5 RNA from the 5' proximal region of SARS-CoV-2 and related coronaviruses

Coronavirus genomes sequester their start codons within stem-loop 5 (SL5), a structured, 5' genomic RNA element. In most alpha- and betacoronaviruses, the secondary structure of SL5 is predicted to contain a four-way junction of helical stems, some of which are capped with UUYYGU hexaloops. Here, using cryogenic electron microscopy (cryo-EM) and computational modeling with biochemically determined secondary structures, we present three-dimensional structures of SL5 from six coronaviruses. The SL5 domain of betacoronavirus severe-acute-respiratory-syndrome-related coronavirus 2 (SARS-CoV-2), resolved at 4.7 Å resolution, exhibits a T-shaped structure, with its UUYYGU hexaloops at opposing ends of a coaxial stack, the T's "arms." Further analysis of SL5 domains from SARS-CoV-1 and MERS (7.1 and 6.4 to 6.9 Å resolution, respectively) indicate that the junction geometry and inter-hexaloop distances are conserved features across these human-infecting betacoronaviruses. The MERS SL5 domain displays an additional tertiary interaction, which is also observed in the non-human-infecting betacoronavirus BtCoV-HKU5 (5.9 to 8.0 Å resolution). SL5s from human-infecting alphacoronaviruses, HCoV-229E and HCoV-NL63 (6.5 and 8.4 to 9.0 Å resolution, respectively), exhibit the same coaxial stacks, including the UUYYGU-capped arms, but with a phylogenetically distinct crossing angle, an X-shape. As such, all SL5 domains studied herein fold into stable tertiary structures with cross-genus similarities and notable differences, with implications for potential protein-binding modes and therapeutic targets.

Coronavirus takeover of host cell translation and intracellular antiviral response: a molecular perspective

Coronaviruses are a group of related RNA viruses that cause respiratory diseases in humans and animals. Understanding the mechanisms of translation regulation during coronaviral infections is critical for developing antiviral therapies and preventing viral spread. Translation of the viral single-stranded RNA genome in the host cell cytoplasm is an essential step in the life cycle of coronaviruses, which affects the cellular mRNA translation landscape in many ways. Here we discuss various viral strategies of translation control, including how members of the Betacoronavirus genus shut down host cell translation and suppress host innate immune functions, as well as the role of the viral non-structural protein 1 (Nsp1) in the process. We also outline the fate of viral RNA, considering stress response mechanisms triggered in infected cells, and describe how unique viral RNA features contribute to programmed ribosomal -1 frameshifting, RNA editing, and translation shutdown evasion.

Programmable modulation of ribosomal frameshifting by mRNA targeting CRISPR-Cas12a system

Minus 1 programmed ribosomal frameshifting (-1 PRF) is a conserved translational regulation event essential for critical biological processes, including the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) replication. Efficient trans-modulation of the structured RNA element crucial to -1 PRF will endow the therapeutic application. Here, we demonstrate that CRISPR RNA can stimulate efficient -1 PRF. Assembled CRISPR-Cas12a, but not CRISPR-Cas9, complex further enhances -1 PRF efficiency through its higher capacity to stall translating ribosomes. We additionally perform CRISPR-Cas12a targeting to impair the SARS-CoV-2 frameshifting pseudoknot structure via a focused screening. We demonstrate that targeting CRISPR-Cas12a results in more than 70% suppression of -1 PRF in vitro and about 50% suppression in mammalian cells. Our results show the expanded function of the CRISPR-Cas12 system in modulating -1 PRF efficiency through stalling ribosomes and deforming frameshifting stimulatory signals, which could serve as a new strategy for future coronavirus pandemics.

Tertiary folds of the SL5 RNA from the 5' proximal region of SARS-CoV-2 and related coronaviruses

Coronavirus genomes sequester their start codons within stem-loop 5 (SL5), a structured, 5' genomic RNA element. In most alpha- and betacoronaviruses, the secondary structure of SL5 is predicted to contain a four-way junction of helical stems, some of which are capped with UUYYGU hexaloops. Here, using cryogenic electron microscopy (cryo-EM) and computational modeling with biochemically-determined secondary structures, we present three-dimensional structures of SL5 from six coronaviruses. The SL5 domain of betacoronavirus SARS-CoV-2, resolved at 4.7 Å resolution, exhibits a T-shaped structure, with its UUYYGU hexaloops at opposing ends of a coaxial stack, the T's "arms." Further analysis of SL5 domains from SARS-CoV-1 and MERS (7.1 and 6.4-6.9 Å resolution, respectively) indicate that the junction geometry and inter-hexaloop distances are conserved features across the studied human-infecting betacoronaviruses. The MERS SL5 domain displays an additional tertiary interaction, which is also observed in the non-human-infecting betacoronavirus BtCoV-HKU5 (5.9-8.0 Å resolution). SL5s from human-infecting alphacoronaviruses, HCoV-229E and HCoV-NL63 (6.5 and 8.4-9.0 Å resolution, respectively), exhibit the same coaxial stacks, including the UUYYGU-capped arms, but with a phylogenetically distinct crossing angle, an X-shape. As such, all SL5 domains studied herein fold into stable tertiary structures with cross-genus similarities, with implications for potential protein-binding modes and therapeutic targets.

Dynamic Counterion Condensation Model Decodes Functional Dynamics of RNA Pseudoknot in SARS-CoV-2: Control of Ion-Mediated Pierced Lasso Topology

The programmed frameshifting stimulatory element, a promising drug target for COVID-19 treatment, involves a RNA pseudoknot (PK) structure. This RNA PK facilitates frameshifting, enabling RNA viruses to translate multiple proteins from a single mRNA, which is a key strategy for their rapid evolution. Overcoming the challenges of capturing large-scale structural changes of RNA under the influence of a dynamic counterion environment (K+ and Mg2+), the study extended the applications of a newly developed dynamic counterion condensation (DCC) model. DCC simulations reveal potential folding pathways of this RNA PK, supported by the experimental findings obtained using optical tweezers. The study elucidates the pivotal role of Mg2+ ions in crafting a lasso-like RNA topology, a novel RNA motif that governs dynamic transitions between the ring-opened and ring-closed states of the RNA. The pierced lasso component guided by Mg2+-mediated interactions orchestrates inward and outward motion fine-tuning tension on the slippery segment, a critical factor for optimizing frameshifting efficiency.

Unwinding the SARS-CoV-2 Ribosomal Frameshifting Pseudoknot with LNA and G-Clamp-Modified Phosphorothioate Oligonucleotides Inhibits Viral Replication

Ribosomal frameshifting (RFS) at the slippery site of SARS-CoV-2 RNA is essential for the biosynthesis of the viral replication machinery. It requires the formation of a pseudoknot (PK) structure near the slippery site and can be inhibited by PK-disrupting oligonucleotide-based antivirals. We obtained and compared three types of such antiviral candidates, namely locked nucleic acids (LNA), LNA-DNA gapmers, and G-clamp-containing phosphorothioates (CPSs) complementary to PK stems. Using optical and electrophoretic methods, we showed that stem 2-targeting oligonucleotide analogs induced PK unfolding at nanomolar concentrations, and this effect was particularly pronounced in the case of LNA. For the leading PK-unfolding LNA and CPS oligonucleotide analogs, we also demonstrated dose-dependent RSF inhibition in dual luciferase assays (DLAs). Finally, we showed that the leading oligonucleotide analogs reduced SARS-CoV-2 replication at subtoxic concentrations in the nanomolar range in two human cell lines. Our findings highlight the promise of PK targeting, illustrate the advantages and limitations of various types of DNA modifications and may promote the future development of oligonucleotide-based antivirals.

Atomistic structure of the SARS-CoV-2 pseudoknot in solution from SAXS-driven molecular dynamics

SARS-CoV-2 depends on -1 programmed ribosomal frameshifting (-1 PRF) to express proteins essential for its replication. The RNA pseudoknot stimulating -1 PRF is thus an attractive drug target. However, the structural models of this pseudoknot obtained from cryo-EM and crystallography differ in some important features, leaving the pseudoknot structure unclear. We measured the solution structure of the pseudoknot using small-angle X-ray scattering (SAXS). The measured profile did not agree with profiles computed from the previously solved structures. Beginning with each of these solved structures, we used the SAXS data to direct all atom molecular dynamics (MD) simulations to improve the agreement in profiles. In all cases, this refinement resulted in a bent conformation that more closely resembled the cryo-EM structures than the crystal structure. Applying the same approach to a point mutant abolishing -1 PRF revealed a notably more bent structure with reoriented helices. This work clarifies the dynamic structures of the SARS-CoV-2 pseudoknot in solution.

High-resolution structure of stem-loop 4 from the 5'-UTR of SARS-CoV-2 solved by solution state NMR

We present the high-resolution structure of stem-loop 4 of the 5'-untranslated region (5_SL4) of the severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) genome solved by solution state nuclear magnetic resonance spectroscopy. 5_SL4 adopts an extended rod-like structure with a single flexible looped-out nucleotide and two mixed tandem mismatches, each composed of a G•U wobble base pair and a pyrimidine•pyrimidine mismatch, which are incorporated into the stem-loop structure. Both the tandem mismatches and the looped-out residue destabilize the stem-loop structure locally. Their distribution along the 5_SL4 stem-loop suggests a role of these non-canonical elements in retaining functionally important structural plasticity in particular with regard to the accessibility of the start codon of an upstream open reading frame located in the RNA's apical loop. The apical loop-although mostly flexible-harbors residual structural features suggesting an additional role in molecular recognition processes. 5_SL4 is highly conserved among the different variants of SARS-CoV-2 and can be targeted by small molecule ligands, which it binds with intermediate affinity in the vicinity of the non-canonical elements within the stem-loop structure.

Structural and Functional Insights into Viral Programmed Ribosomal Frameshifting

Protein synthesis by the ribosome is the final stage of biological information transfer and represents an irreversible commitment to gene expression. Accurate translation of messenger RNA is therefore essential to all life, and spontaneous errors by the translational machinery are highly infrequent (∼1/100,000 codons). Programmed -1 ribosomal frameshifting (-1PRF) is a mechanism in which the elongating ribosome is induced at high frequency to slip backward by one nucleotide at a defined position and to continue translation in the new reading frame. This is exploited as a translational regulation strategy by hundreds of RNA viruses, which rely on -1PRF during genome translation to control the stoichiometry of viral proteins. While early investigations of -1PRF focused on virological and biochemical aspects, the application of X-ray crystallography and cryo-electron microscopy (cryo-EM), and the advent of deep sequencing and single-molecule approaches have revealed unexpected structural diversity and mechanistic complexity. Molecular players from several model systems have now been characterized in detail, both in isolation and, more recently, in the context of the elongating ribosome. Here we provide a summary of recent advances and discuss to what extent a general model for -1PRF remains a useful way of thinking.

Terahertz waves regulate the mechanical unfolding of tau pre-mRNA hairpins

Intermolecular interactions, including hydrogen bonds, dominate the pairing and unpairing of nucleic acid chains in the transfer process of genetic information. The energy of THz waves just matches with the weak interactions, so THz waves may interact with biomolecules. Here, the dynamic effects of THz electromagnetic (EM) waves on the mechanical unfolding process of RNA hairpins (WT-30nt and its mutants, rHP, SARS-CoV-2, and SRV-1 SF206) are investigated using steered molecular dynamics (SMD) simulations. The results show that THz waves can either promote the unfolding of the double helix of the RNA hairpin during the initial unfolding phase (4-21.8 THz) or significantly enhance (23.8 and 25.5 THz) or weaken (37.4 and 41.2 THz) its structural stability during unfolding. Our findings have important implications for applying THz waves to regulate dynamic deconvolution processes, such as gene replication, transcription, and translation.

Comparative analysis of RNA secondary structure accuracy on predicted RNA 3D models

RNA structure is conformationally dynamic, and accurate all-atom tertiary (3D) structure modeling of RNA remains challenging with the prevailing tools. Secondary structure (2D) information is the standard prerequisite for most RNA 3D modeling. Despite several 2D and 3D structure prediction tools proposed in recent years, one of the challenges is to choose the best combination for accurate RNA 3D structure prediction. Here, we benchmarked seven small RNA PDB structures (40 to 90 nucleotides) with different topologies to understand the effects of different 2D structure predictions on the accuracy of 3D modeling. The current study explores the blind challenge of 2D to 3D conversions and highlights the performances of de novo RNA 3D modeling from their predicted 2D structure constraints. Our results show that conformational sampling-based methods such as SimRNA and IsRNA1 depend less on 2D accuracy, whereas motif-based methods account for 2D evidence. Our observations illustrate the disparities in available 3D and 2D prediction methods and may further offer insights into developing topology-specific or family-specific RNA structure prediction pipelines.

RNA structure determination: From 2D to 3D

RNA molecules serve a wide range of functions that are closely linked to their structures. The basic structural units of RNA consist of single- and double-stranded regions. In order to carry out advanced functions such as catalysis and ligand binding, certain types of RNAs can adopt higher-order structures. The analysis of RNA structures has progressed alongside advancements in structural biology techniques, but it comes with its own set of challenges and corresponding solutions. In this review, we will discuss recent advances in RNA structure analysis techniques, including structural probing methods, X-ray crystallography, nuclear magnetic resonance, cryo-electron microscopy, and small-angle X-ray scattering. Often, a combination of multiple techniques is employed for the integrated analysis of RNA structures. We also survey important RNA structures that have been recently determined using various techniques.

Small Molecules Targeting Viral RNA

The majority of antivirals available target viral proteins; however, RNA is emerging as a new and promising antiviral target due to the presence of highly structured RNA in viral genomes fundamental for their replication cycle. Here, we discuss methods for the identification of RNA-targeting compounds, starting from the determination of RNA structures either from purified RNA or in living cells, followed by in silico screening on RNA and phenotypic assays to evaluate viral inhibition. Moreover, we review the small molecules known to target the programmed ribosomal frameshifting element of SARS-CoV-2, the internal ribosomal entry site of different viruses, and RNA elements of HIV.

Discovery of Small Molecules Targeting the Frameshifting Element RNA in SARS-CoV-2 Viral Genome

Targeting structured RNA elements in the SARS-CoV-2 viral genome with small molecules is an attractive strategy for pharmacological control over viral replication. In this work, we report the discovery of small molecules that target the frameshifting element (FSE) in the SARS-CoV-2 RNA genome using high-throughput small-molecule microarray (SMM) screening. A new class of aminoquinazoline ligands for the SARS-CoV-2 FSE are synthesized and characterized using multiple orthogonal biophysical assays and structure-activity relationship (SAR) studies. This work reveals compounds with mid-micromolar binding affinity (KD = 60 ± 6 μM) to the FSE RNA and supports a binding mode distinct from previously reported FSE binders MTDB and merafloxacin. In addition, compounds are active in in vitro dual-luciferase and in-cell dual-fluorescent-reporter frameshifting assays, highlighting the promise of targeting structured elements of RNAs with druglike compounds to alter expression of viral proteins.

Modelling the structures of frameshift-stimulatory pseudoknots from representative bat coronaviruses

Coronaviruses (CoVs) use -1 programmed ribosomal frameshifting stimulated by RNA pseudoknots in the viral genome to control expression of enzymes essential for replication, making CoV pseudoknots a promising target for anti-coronaviral drugs. Bats represent one of the largest reservoirs of CoVs and are the ultimate source of most CoVs infecting humans, including those causing SARS, MERS, and COVID-19. However, the structures of bat-CoV frameshift-stimulatory pseudoknots remain largely unexplored. Here we use a combination of blind structure prediction followed by all-atom molecular dynamics simulations to model the structures of eight pseudoknots that, together with the SARS-CoV-2 pseudoknot, are representative of the range of pseudoknot sequences in bat CoVs. We find that they all share some key qualitative features with the pseudoknot from SARS-CoV-2, notably the presence of conformers with two distinct fold topologies differing in whether or not the 5' end of the RNA is threaded through a junction, and similar conformations for stem 1. However, they differed in the number of helices present, with half sharing the 3-helix architecture of the SARS-CoV-2 pseudoknot but two containing 4 helices and two others only 2. These structure models should be helpful for future work studying bat-CoV pseudoknots as potential therapeutic targets.

A riboswitch separated from its ribosome-binding site still regulates translation

Riboswitches regulate downstream gene expression by binding cellular metabolites. Regulation of translation initiation by riboswitches is posited to occur by metabolite-mediated sequestration of the Shine-Dalgarno sequence (SDS), causing bypass by the ribosome. Recently, we solved a co-crystal structure of a prequeuosine1-sensing riboswitch from Carnobacterium antarcticum that binds two metabolites in a single pocket. The structure revealed that the second nucleotide within the gene-regulatory SDS, G34, engages in a crystal contact, obscuring the molecular basis of gene regulation. Here, we report a co-crystal structure wherein C10 pairs with G34. However, molecular dynamics simulations reveal quick dissolution of the pair, which fails to reform. Functional and chemical probing assays inside live bacterial cells corroborate the dispensability of the C10-G34 pair in gene regulation, leading to the hypothesis that the compact pseudoknot fold is sufficient for translation attenuation. Remarkably, the C. antarcticum aptamer retained significant gene-regulatory activity when uncoupled from the SDS using unstructured spacers up to 10 nucleotides away from the riboswitch-akin to steric-blocking employed by sRNAs. Accordingly, our work reveals that the RNA fold regulates translation without SDS sequestration, expanding known riboswitch-mediated gene-regulatory mechanisms. The results infer that riboswitches exist wherein the SDS is not embedded inside a stable fold.

Length-dependent motions of SARS-CoV-2 frameshifting RNA pseudoknot and alternative conformations suggest avenues for frameshifting suppression

The SARS-CoV-2 frameshifting element (FSE), a highly conserved mRNA region required for correct translation of viral polyproteins, defines an excellent therapeutic target against Covid-19. As discovered by our prior graph-theory analysis with SHAPE experiments, the FSE adopts a heterogeneous, length-dependent conformational landscape consisting of an assumed 3-stem H-type pseudoknot (graph motif 3_6), and two alternative motifs (3_3 and 3_5). Here, for the first time, we build and simulate, by microsecond molecular dynamics, 30 models for all three motifs plus motif-stabilizing mutants at different lengths. Our 3_6 pseudoknot systems, which agree with experimental structures, reveal interconvertible L and linear conformations likely related to ribosomal pausing and frameshifting. The 3_6 mutant inhibits this transformation and could hamper frameshifting. Our 3_3 systems exhibit length-dependent stem interactions that point to a potential transition pathway connecting the three motifs during ribosomal elongation. Together, our observations provide new insights into frameshifting mechanisms and anti-viral strategies.

Viral and cellular translation during SARS‐CoV‐2 infection

SARS‐CoV‐2 is a betacoronavirus that emerged in China in December 2019 and which is the causative agent of the Covid‐19 pandemic. This enveloped virus contains a large positive‐sense single‐stranded RNA genome. In this review, we summarize the current knowledge on the molecular mechanisms for the translation of both viral transcripts and cellular messenger RNAs. Non‐structural proteins are encoded by the genomic RNA and are produced in the early steps of infection. In contrast, the structural proteins are produced from subgenomic RNAs that are translated in the late phase of the infectious program. Non‐structural protein 1 (NSP1) is a key molecule that regulates both viral and cellular translation. In addition, NSP1 interferes with multiple steps of the interferon I pathway and thereby blocks host antiviral responses. Therefore, NSP1 is a drug target of choice for the development of antiviral therapies.

Translational Control of COVID-19 and Its Therapeutic Implication

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of COVID-19, which has broken out worldwide for more than two years. However, due to limited treatment, new cases of infection are still rising. Therefore, there is an urgent need to understand the basic molecular biology of SARS-CoV-2 to control this virus. SARS-CoV-2 replication and spread depend on the recruitment of host ribosomes to translate viral messenger RNA (mRNA). To ensure the translation of their own mRNAs, the SARS-CoV-2 has developed multiple strategies to globally inhibit the translation of host mRNAs and block the cellular innate immune response. This review provides a comprehensive picture of recent advancements in our understanding of the molecular basis and complexity of SARS-CoV-2 protein translation. Specifically, we summarize how this viral infection inhibits host mRNA translation to better utilize translation elements for translation of its own mRNA. Finally, we discuss the potential of translational components as targets for therapeutic interventions.

Short Oligonucleotides Facilitate Co-transcriptional Labeling of RNA at Specific Positions

Labeling RNA molecules at specific positions is critical for RNA research and applications. Such methods are in high demand but still a challenge, especially those that enable native co-synthesis rather than post-synthesis labeling of long RNAs. The method we developed in this work meets these requirements, in which a leader RNA is extended on the hybrid solid-liquid phase by an engineered transcriptional complex following the pause-restart mode. A custom-designed short oligonucleotide is used to functionalize the engineered complex. This remarkable co-transcriptional labeling method incorporates labels into RNAs in high yields with great flexibility. We demonstrate the method by successfully introducing natural modifications, a fluorescent nucleotide analogue and a donor-acceptor fluorophore pair to specific sites located at an internal loop, a pseudoknot, a junction, a helix, and the middle of consecutive identical nucleotides of various RNAs. This newly developed method overcomes efficiency and position-choosing constraints that have hampered routine strategies to label RNAs beyond 200 nucleotides (nt).

RT-qPCR as a screening platform for mutational and small molecule impacts on structural stability of RNA tertiary structures

The exponential increase in the discovery and characterization of RNA tertiary structures has highlighted their active role in a variety of human diseases, yet often their interactome and specific function remain unknown. Small molecules offer opportunities to both decode these cellular roles and develop therapeutics, however there are few examples of small molecules that target biologically relevant RNA tertiary structures. While RNA triple helices are a particularly attractive target, discovery of triple helix modulators has been hindered by the lack of correlation between small molecule affinity and effect on structural modulation, thereby limiting the utility of affinity-based screening as a primary filtering method. To address this challenge, we developed a high-throughput RT-qPCR screening platform that reports on the effect of mutations and additives, such as small molecules, on the stability of triple helices. Using the 3′-end of the oncogenic long non-coding RNA MALAT1 as a proof-of-concept, we demonstrated the applicability of both a two-step and a one-pot method to assess the impact of mutations and small molecules on the stability of the triple helix. We demonstrated the adaptability of the assay to diverse RNA tertiary structures by applying it to the SARS-CoV-2 pseudoknot, a key viral RNA structure recently identified as an attractive therapeutic target for the development of antivirals. Employment of a functional high-throughput assay as a primary screen will significantly expedite the discovery of probes that modulate the structural landscape of RNA structures and, consequently, help gain insight into the roles of these pervasive structures.

RT-qPCR can be harnessed as a small molecule screening platform to read out the effect of small molecules on the structural stability of a variety of RNA targets.