Genome-wide mapping of SARS-CoV-2 RNA structures identifies therapeutically-relevant elements

In vivo structure and dynamics of the SARS-CoV-2 RNA genome

The dynamics of SARS-CoV-2 RNA structure and their functional relevance are largely unknown. Here we develop a simplified SPLASH assay and comprehensively map the in vivo RNA-RNA interactome of SARS-CoV-2 genome across viral life cycle. We report canonical and alternative structures including 5'-UTR and 3'-UTR, frameshifting element (FSE) pseudoknot and genome cyclization in both cells and virions. We provide direct evidence of interactions between Transcription Regulating Sequences, which facilitate discontinuous transcription. In addition, we reveal alternative short and long distance arches around FSE. More importantly, we find that within virions, while SARS-CoV-2 genome RNA undergoes intensive compaction, genome domains remain stable but with strengthened demarcation of local domains and weakened global cyclization. Taken together, our analysis reveals the structural basis for the regulation of replication, discontinuous transcription and translational frameshifting, the alternative conformations and the maintenance of global genome organization during the whole life cycle of SARS-CoV-2, which we anticipate will help develop better antiviral strategies.

Chemical and Enzymatic Probing of Viral RNAs: From Infancy to Maturity and Beyond

RNA molecules are key players in a variety of biological events, and this is particularly true for viral RNAs. To better understand the replication of those pathogens and try to block them, special attention has been paid to the structure of their RNAs. Methods to probe RNA structures have been developed since the 1960s; even if they have evolved over the years, they are still in use today and provide useful information on the folding of RNA molecules, including viral RNAs. The aim of this review is to offer a historical perspective on the structural probing methods used to decipher RNA structures before the development of the selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) methodology and to show how they have influenced the current probing techniques. Actually, these technological breakthroughs, which involved advanced detection methods, were made possible thanks to the development of next-generation sequencing (NGS) but also to the previous works accumulated in the field of structural RNA biology. Finally, we will also discuss how high-throughput SHAPE (hSHAPE) paved the way for the development of sophisticated RNA structural techniques.

Characterisation of SARS-CoV-2 clades based on signature SNPs unveils continuous evolution

Since the emergence of SARS-CoV-2 in Wuhan, China more than a year ago, it has spread across the world in a very short span of time. Although, different forms of vaccines are being rolled out for vaccination programs around the globe, the mutation of the virus is still a cause of concern among the research communities. Hence, it is important to study the constantly evolving virus and its strains in order to provide a much more stable form of cure. This fact motivated us to conduct this research where we have initially carried out multiple sequence alignment of 15359 and 3033 global dataset without Indian and the dataset of exclusive Indian SARS-CoV-2 genomes respectively, using MAFFT. Subsequently, phylogenetic analyses are performed using Nextstrain to identify virus clades. Consequently, the virus strains are found to be distributed among 5 major clades or clusters viz. 19A, 19B, 20A, 20B and 20C. Thereafter, mutation points as SNPs are identified in each clade. Henceforth, from each clade top 10 signature SNPs are identified based on their frequency i.e. number of occurrences in the virus genome. As a result, 50 such signature SNPs are individually identified for global dataset without Indian and dataset of exclusive Indian SARS-CoV-2 genomes respectively. Out of each 50 signature SNPs, 39 and 41 unique SNPs are identified among which 25 non-synonymous signature SNPs (out of 39) resulted in 30 amino acid changes in protein while 27 changes in amino acid are identified from 22 non-synonymous signature SNPs (out of 41). These 30 and 27 amino acid changes for the non-synonymous signature SNPs are visualised in their respective protein structure as well. Finally, in order to judge the characteristics of the identified clades, the non-synonymous signature SNPs are considered to evaluate the changes in proteins as biological functions with the sequences using PROVEAN and PolyPhen-2 while I-Mutant 2.0 is used to evaluate their structural stability. As a consequence, for global dataset without Indian sequences, G251V in ORF3a in clade 19A, F308Y and G196V in NSP4 and ORF3a in 19B are the unique amino acid changes which are responsible for defining each clade as they are all deleterious and unstable. Such changes which are common for both global dataset without Indian and dataset of exclusive Indian sequences are R203M in Nucleocapsid for 20B, T85I and Q57H in NSP2 and ORF3a respectively for 20C while for exclusive Indian sequences such unique changes are A97V in RdRp, G339S and G339C in NSP2 in 19A and Q57H in ORF3a in 20A.

Potential Achilles heels of SARS-CoV-2 are best displayed by the base order-dependent component of RNA folding energy

The base order-dependent component of folding energy has revealed a highly conserved region in HIV-1 genomes that associates with RNA structure. This corresponds to a packaging signal that is recognized by the nucleocapsid domain of the Gag polyprotein. Long viewed as a potential HIV-1 "Achilles heel," the signal can be targeted by a new antiviral compound. Although SARS-CoV-2 differs in many respects from HIV-1, the same technology displays regions with a high base order-dependent folding energy component, which are also highly conserved. This indicates structural invariance (SI) sustained by natural selection. While the regions are often also protein-encoding (e. g. NSP3, ORF3a), we suggest that their nucleic acid level functions can be considered potential "Achilles heels" for SARS-CoV-2, perhaps susceptible to therapies like those envisaged for AIDS. The ribosomal frameshifting element scored well, but higher SI scores were obtained in other regions, including those encoding NSP13 and the nucleocapsid (N) protein.

Rapid structure-function insights via hairpin-centric analysis of big RNA structure probing datasets

The functions of RNA are often tied to its structure, hence analyzing structure is of significant interest when studying cellular processes. Recently, large-scale structure probing (SP) studies have enabled assessment of global structure-function relationships via standard data summarizations or local folding. Here, we approach structure quantification from a hairpin-centric perspective where putative hairpins are identified in SP datasets and used as a means to capture local structural effects. This has the advantage of rapid processing of big (e.g. transcriptome-wide) data as RNA folding is circumvented, yet it captures more information than simple data summarizations. We reformulate a statistical learning algorithm we previously developed to significantly improve precision of hairpin detection, then introduce a novel nucleotide-wise measure, termed the hairpin-derived structure level (HDSL), which captures local structuredness by accounting for the presence of likely hairpin elements. Applying HDSL to data from recent studies recapitulates, strengthens and expands on their findings which were obtained by more comprehensive folding algorithms, yet our analyses are orders of magnitude faster. These results demonstrate that hairpin detection is a promising avenue for global and rapid structure-function analysis, furthering our understanding of RNA biology and the principal features which drive biological insights from SP data.

Cryo-EM and antisense targeting of the 28-kDa frameshift stimulation element from the SARS-CoV-2 RNA genome

Drug discovery campaigns against COVID-19 are beginning to target the SARS-CoV-2 RNA genome. The highly conserved frameshift stimulation element (FSE), required for balanced expression of viral proteins, is a particularly attractive SARS-CoV-2 RNA target. Here we present a 6.9 Å resolution cryo-EM structure of the FSE (88 nucleotides, ~28 kDa), validated through an RNA nanostructure tagging method. The tertiary structure presents a topologically complex fold in which the 5' end is threaded through a ring formed inside a three-stem pseudoknot. Guided by this structure, we develop antisense oligonucleotides that impair FSE function in frameshifting assays and knock down SARS-CoV-2 virus replication in A549-ACE2 cells at 100 nM concentration.

Structure-based design of antisense oligonucleotides that inhibit SARS-CoV-2 replication

Antisense oligonucleotides (ASOs) are an emerging class of drugs that target RNAs. Current ASO designs strictly follow the rule of Watson-Crick base pairing along target sequences. However, RNAs often fold into structures that interfere with ASO hybridization. Here we developed a structure-based ASO design method and applied it to target severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Our method makes sure that ASO binding is compatible with target structures in three-dimensional (3D) space by employing structural design templates. These 3D-ASOs recognize the shapes and hydrogen bonding patterns of targets via tertiary interactions, achieving enhanced affinity and specificity. We designed 3D-ASOs that bind to the frameshift stimulation element and transcription regulatory sequence of SARS-CoV-2 and identified lead ASOs that strongly inhibit viral replication in human cells. We further optimized the lead sequences and characterized structure-activity relationship. The 3D-ASO technology helps fight coronavirus disease-2019 and is broadly applicable to ASO drug development.

Exploring the Druggability of Conserved RNA Regulatory Elements in the SARS-CoV-2 Genome

SARS-CoV-2 contains a positive single-stranded RNA genome of approximately 30 000 nucleotides. Within this genome, 15 RNA elements were identified as conserved between SARS-CoV and SARS-CoV-2. By nuclear magnetic resonance (NMR) spectroscopy, we previously determined that these elements fold independently, in line with data from in vivo and ex-vivo structural probing experiments. These elements contain non-base-paired regions that potentially harbor ligand-binding pockets. Here, we performed an NMR-based screening of a poised fragment library of 768 compounds for binding to these RNAs, employing three different 1 H-based 1D NMR binding assays. The screening identified common as well as RNA-element specific hits. The results allow selection of the most promising of the 15 RNA elements as putative drug targets. Based on the identified hits, we derive key functional units and groups in ligands for effective targeting of the RNA of SARS-CoV-2.

The SARS-CoV-2 Programmed -1 Ribosomal Frameshifting Element Crystal Structure Solved to 2.09 Å Using Chaperone-Assisted RNA Crystallography

The programmed -1 ribosomal frameshifting element (PFSE) of SARS-CoV-2 is a well conserved structured RNA found in all coronaviruses' genomes. By adopting a pseudoknot structure in the presence of the ribosome, the PFSE promotes a ribosomal frameshifting event near the stop codon of the first open reading frame Orf1a during translation of the polyprotein pp1a. Frameshifting results in continuation of pp1a via a new open reading frame, Orf1b, that produces the longer pp1ab polyprotein. Polyproteins pp1a and pp1ab produce nonstructural proteins NSPs 1-10 and NSPs 1-16, respectively, which contribute vital functions during the viral life cycle and must be present in the proper stoichiometry. Both drugs and sequence alterations that affect the stability of the -1 programmed ribosomal frameshifting element disrupt the stoichiometry of the NSPs produced, which compromise viral replication. For this reason, the -1 programmed frameshifting element is considered a promising drug target. Using chaperone assisted RNA crystallography, we successfully crystallized and solved the three-dimensional structure of the PFSE. We observe a three-stem H-type pseudoknot structure with the three stems stacked in a vertical orientation stabilized by two triple base pairs at the stem 1/stem 2 and stem 1/stem 3 junctions. This structure provides a new conformation of PFSE distinct from the bent conformations inferred from midresolution cryo-EM models and provides a high-resolution framework for mechanistic investigations and structure-based drug design.

Characterization of SARS2 Nsp15 nuclease activity reveals it's mad about U

Nsp15 is a uridine specific endoribonuclease that coronaviruses employ to cleave viral RNA and evade host immune defense systems. Previous structures of Nsp15 from across Coronaviridae revealed that Nsp15 assembles into a homo-hexamer and has a conserved active site similar to RNase A. Beyond a preference for cleaving RNA 3′ of uridines, it is unknown if Nsp15 has any additional substrate preferences. Here, we used cryo-EM to capture structures of Nsp15 bound to RNA in pre- and post-cleavage states. The structures along with molecular dynamics and biochemical assays revealed critical residues involved in substrate specificity, nuclease activity, and oligomerization. Moreover, we determined how the sequence of the RNA substrate dictates cleavage and found that outside of polyU tracts, Nsp15 has a strong preference for purines 3′ of the cleaved uridine. This work advances our understanding of how Nsp15 recognizes and processes viral RNA, and will aid in the development of new anti-viral therapeutics.

Supramolecular Cylinders Target Bulge Structures in the 5' UTR of the RNA Genome of SARS-CoV-2 and Inhibit Viral Replication

The untranslated regions (UTRs) of viral genomes contain a variety of conserved yet dynamic structures crucial for viral replication, providing drug targets for the development of broad spectrum anti-virals. We combine in vitro RNA analysis with molecular dynamics simulations to build the first 3D models of the structure and dynamics of key regions of the 5' UTR of the SARS-CoV-2 genome. Furthermore, we determine the binding of metallo-supramolecular helicates (cylinders) to this RNA structure. These nano-size agents are uniquely able to thread through RNA junctions and we identify their binding to a 3-base bulge and the central cross 4-way junction located in stem loop 5. Finally, we show these RNA-binding cylinders suppress SARS-CoV-2 replication, highlighting their potential as novel anti-viral agents.

No species-level losses of s2m suggests critical role in replication of SARS-related coronaviruses

The genetic element s2m has been acquired through horizontal transfer by many distantly related viruses, including the SARS-related coronaviruses. Here we show that s2m is evolutionarily conserved in these viruses. Though several lineages of severe acute respiratory syndrome coronavirus 2 (SARS‑CoV‑2) devoid of the element can be found, these variants seem to have been short lived, indicating that they were less evolutionary fit than their s2m-containing counterparts. On a species-level, however, there do not appear to be any losses and this pattern strongly suggests that the s2m element is essential to virus replication in SARS-CoV-2 and related viruses. Further experiments are needed to elucidate the function of s2m.

Supramolecular Cylinders Target Bulge Structures in the 5' UTR of the RNA Genome of SARS-CoV-2 and Inhibit Viral Replication*

The untranslated regions (UTRs) of viral genomes contain a variety of conserved yet dynamic structures crucial for viral replication, providing drug targets for the development of broad spectrum anti-virals. We combine in vitro RNA analysis with molecular dynamics simulations to build the first 3D models of the structure and dynamics of key regions of the 5' UTR of the SARS-CoV-2 genome. Furthermore, we determine the binding of metallo-supramolecular helicates (cylinders) to this RNA structure. These nano-size agents are uniquely able to thread through RNA junctions and we identify their binding to a 3-base bulge and the central cross 4-way junction located in stem loop 5. Finally, we show these RNA-binding cylinders suppress SARS-CoV-2 replication, highlighting their potential as novel anti-viral agents.

Structural dynamics of single SARS-CoV-2 pseudoknot molecules reveal topologically distinct conformers

The RNA pseudoknot that stimulates programmed ribosomal frameshifting in SARS-CoV-2 is a possible drug target. To understand how it responds to mechanical tension applied by ribosomes, thought to play a key role during frameshifting, we probe its structural dynamics using optical tweezers. We find that it forms multiple structures: two pseudoknotted conformers with different stability and barriers, and alternative stem-loop structures. The pseudoknotted conformers have distinct topologies, one threading the 5′ end through a 3-helix junction to create a knot-like fold, the other with unthreaded 5′ end, consistent with structures observed via cryo-EM and simulations. Refolding of the pseudoknotted conformers starts with stem 1, followed by stem 3 and lastly stem 2; Mg²⁺ ions are not required, but increase pseudoknot mechanical rigidity and favor formation of the knot-like conformer. These results resolve the SARS-CoV-2 frameshift signal folding mechanism and highlight its conformational heterogeneity, with important implications for structure-based drug-discovery efforts.

The RNA pseudoknot of SARS-CoV-2 promotes -1 programmed ribosomal frameshifting. Here the authors use single molecule force spectroscopy to study the folding of this pseudoknot, showing that it forms at least two different pseudoknot conformers with distinct fold topologies.

Molecular docking and dynamics study to explore phytochemical ligand molecules against the main protease of SARS-CoV-2 from extensive phytochemical datasets

BACKGROUND

The high transmission and pathogenicity of SARS-CoV-2 has led to a pandemic that has halted the world's economy and health. The newly evolved strains and scarcity of vaccines has worsened the situation. The main protease (M^pro) of SARS-CoV-2 can act as a potential target due to its role in viral replication and conservation level.

METHODS

In this study, we have enlisted more than 1100 phytochemicals from Asian plants based on deep literature mining. The compounds library was screened against the M^pro of SARS-CoV-2.

RESULTS

The selected three ligands, Flemichin, Delta-Oleanolic acid, and Emodin 1-O-beta-D-glucoside had a binding energy of -8.9, -8.9, -8.7 KJ/mol respectively. The compounds bind to the active groove of the main protease at; Cys145, Glu166, His41, Met49, Pro168, Met165, Gln189. The multiple descriptors from the simulation study; root mean square deviation, root mean square fluctuation, radius of gyration, hydrogen bond, solvent accessible surface area confirms the stable nature of the protein-ligand complexes. Furthermore, post-md analysis confirms the rigidness in the docked poses over the simulation trajectories.

CONCLUSIONS

Our combinatorial drug design approaches may help researchers to identify suitable drug candidates against SARS-CoV-2.

To Knot or Not to Knot: Multiple Conformations of the SARS-CoV-2 Frameshifting RNA Element

The SARS-CoV-2 frameshifting RNA element (FSE) is an excellent target for therapeutic intervention against Covid-19. This small gene element employs a shifting mechanism to pause and backtrack the ribosome during translation between Open Reading Frames 1a and 1b, which code for viral polyproteins. Any interference with this process has a profound effect on viral replication and propagation. Pinpointing the structures adapted by the FSE and associated structural transformations involved in frameshifting has been a challenge. Using our graph-theory-based modeling tools for representing RNA secondary structures, "RAG" (RNA-As-Graphs), and chemical structure probing experiments, we show that the 3-stem H-type pseudoknot (3_6 dual graph), long assumed to be the dominant structure, has a viable alternative, an HL-type 3-stem pseudoknot (3_3) for longer constructs. In addition, an unknotted 3-way junction RNA (3_5) emerges as a minor conformation. These three conformations share Stems 1 and 3, while the different Stem 2 may be involved in a conformational switch and possibly associations with the ribosome during translation. For full-length genomes, a stem-loop motif (2_2) may compete with these forms. These structural and mechanistic insights advance our understanding of the SARS-CoV-2 frameshifting process and concomitant virus life cycle, and point to three avenues of therapeutic intervention.

To knot or not to knot: Multiple conformations of the SARS-CoV-2 frameshifting RNA element

The SARS-CoV-2 frameshifting RNA element (FSE) is an excellent target for therapeutic intervention against Covid-19. This small gene element employs a shifting mechanism to pause and backtrack the ribosome during translation between Open Reading Frames 1a and 1b, which code for viral polyproteins. Any interference with this process has profound effect on viral replication and propagation. Pinpointing the structures adapted by the FSE and associated structural transformations involved in frameshifting has been a challenge. Using our graph-theory-based modeling tools for representing RNA secondary structures, "RAG" (RNA-As-Graphs), and chemical structure probing experiments, we show that the 3-stem H-type pseudoknot (3_6 dual graph), long assumed to be the dominant structure has a viable alternative, an HL-type 3-stem pseudoknot (3_3) for longer constructs. In addition, an unknotted 3-way junction RNA (3_5) emerges as a minor conformation. These three conformations share Stems 1 and 3, while the different Stem 2 may be involved in a conformational switch and possibly associations with the ribosome during translation. For full-length genomes, a stem-loop motif (2_2) may compete with these forms. These structural and mechanistic insights advance our understanding of the SARS-CoV-2 frameshifting process and concomitant virus life cycle, and point to three avenues of therapeutic intervention.

Targeting the Conserved Stem Loop 2 Motif in the SARS-CoV-2 Genome

RNA structural elements occur in numerous single-stranded positive-sense RNA viruses. The stem-loop 2 motif (s2m) is one such element with an unusually high degree of sequence conservation, being found in the 3' untranslated region (UTR) in the genomes of many astroviruses, some picornaviruses and noroviruses, and a variety of coronaviruses, including severe acute respiratory syndrome coronavirus (SARS-CoV) and SARS-CoV-2. The evolutionary conservation and its occurrence in all viral subgenomic transcripts imply a key role for s2m in the viral infection cycle. Our findings indicate that the element, while stably folded, can nonetheless be invaded and remodeled spontaneously by antisense oligonucleotides (ASOs) that initiate pairing in exposed loops and trigger efficient sequence-specific RNA cleavage in reporter assays. ASOs also act to inhibit replication in an astrovirus replicon model system in a sequence-specific, dose-dependent manner and inhibit SARS-CoV-2 replication in cell culture. Our results thus permit us to suggest that the s2m element is readily targeted by ASOs, which show promise as antiviral agents. IMPORTANCE The highly conserved stem-loop 2 motif (s2m) is found in the genomes of many RNA viruses, including SARS-CoV-2. Our findings indicate that the s2m element can be targeted by antisense oligonucleotides. The antiviral potential of this element represents a promising start for further research into targeting conserved elements in RNA viruses.

The architecture of the SARS-CoV-2 RNA genome inside virion

SARS-CoV-2 carries the largest single-stranded RNA genome and is the causal pathogen of the ongoing COVID-19 pandemic. How the SARS-CoV-2 RNA genome is folded in the virion remains unknown. To fill the knowledge gap and facilitate structure-based drug development, we develop a virion RNA in situ conformation sequencing technology, named vRIC-seq, for probing viral RNA genome structure unbiasedly. Using vRIC-seq data, we reconstruct the tertiary structure of the SARS-CoV-2 genome and reveal a surprisingly "unentangled globule" conformation. We uncover many long-range duplexes and higher-order junctions, both of which are under purifying selections and contribute to the sequential package of the SARS-CoV-2 genome. Unexpectedly, the D614G and the other two accompanying mutations may remodel duplexes into more stable forms. Lastly, the structure-guided design of potent small interfering RNAs can obliterate the SARS-CoV-2 in Vero cells. Overall, our work provides a framework for studying the genome structure, function, and dynamics of emerging deadly RNA viruses.

Discovery of 17 conserved structural RNAs in fungi

Many non-coding RNAs with known functions are structurally conserved: their intramolecular secondary and tertiary interactions are maintained across evolutionary time. Consequently, the presence of conserved structure in multiple sequence alignments can be used to identify candidate functional non-coding RNAs. Here, we present a bioinformatics method that couples iterative homology search with covariation analysis to assess whether a genomic region has evidence of conserved RNA structure. We used this method to examine all unannotated regions of five well-studied fungal genomes (Saccharomyces cerevisiae, Candida albicans, Neurospora crassa, Aspergillus fumigatus, and Schizosaccharomyces pombe). We identified 17 novel structurally conserved non-coding RNA candidates, which include four H/ACA box small nucleolar RNAs, four intergenic RNAs and nine RNA structures located within the introns and untranslated regions (UTRs) of mRNAs. For the two structures in the 3' UTRs of the metabolic genes GLY1 and MET13, we performed experiments that provide evidence against them being eukaryotic riboswitches.

Conformational Ensembles of Noncoding Elements in the SARS-CoV-2 Genome from Molecular Dynamics Simulations

The 5' untranslated region (UTR) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome is a conserved, functional and structured genomic region consisting of several RNA stem-loop elements. While the secondary structure of such elements has been determined experimentally, their three-dimensional structures are not known yet. Here, we predict structure and dynamics of five RNA stem loops in the 5'-UTR of SARS-CoV-2 by extensive atomistic molecular dynamics simulations, more than 0.5 ms of aggregate simulation time, in combination with enhanced sampling techniques. We compare simulations with available experimental data, describe the resulting conformational ensembles, and identify the presence of specific structural rearrangements in apical and internal loops that may be functionally relevant. Our atomic-detailed structural predictions reveal a rich dynamics in these RNA molecules, could help the experimental characterization of these systems, and provide putative three-dimensional models for structure-based drug design studies.

A map of the SARS-CoV-2 RNA structurome

SARS-CoV-2 has exploded throughout the human population. To facilitate efforts to gain insights into SARS-CoV-2 biology and to target the virus therapeutically, it is essential to have a roadmap of likely functional regions embedded in its RNA genome. In this report, we used a bioinformatics approach, ScanFold, to deduce the local RNA structural landscape of the SARS-CoV-2 genome with the highest likelihood of being functional. We recapitulate previously-known elements of RNA structure and provide a model for the folding of an essential frameshift signal. Our results find that SARS-CoV-2 is greatly enriched in unusually stable and likely evolutionarily ordered RNA structure, which provides a large reservoir of potential drug targets for RNA-binding small molecules. Results are enhanced via the re-analyses of publicly-available genome-wide biochemical structure probing datasets that are broadly in agreement with our models. Additionally, ScanFold was updated to incorporate experimental data as constraints in the analysis to facilitate comparisons between ScanFold and other RNA modelling approaches. Ultimately, ScanFold was able to identify eight highly structured/conserved motifs in SARS-CoV-2 that agree with experimental data, without explicitly using these data. All results are made available via a public database (the RNAStructuromeDB: https://structurome.bb.iastate.edu/sars-cov-2) and model comparisons are readily viewable at https://structurome.bb.iastate.edu/sars-cov-2-global-model-comparisons.

SHAPE Directed Discovery of New Functions in Large RNAs

RNA lies upstream of nearly all biology and functions as the central conduit of information exchange in all cells. RNA molecules encode information both in their primary sequences and in complex structures that form when an RNA folds back on itself. From the time of discovery of mRNA in the late 1950s until quite recently, we had only a rudimentary understanding of RNA structure across vast regions of most messenger and noncoding RNAs. This deficit is now rapidly being addressed, especially by selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry, mutational profiling (MaP), and closely related platform technologies that, collectively, create chemical microscopes for RNA. These technologies make it possible to interrogate RNA structure, quantitatively, at nucleotide resolution, and at large scales, for entire mRNAs, noncoding RNAs, and viral RNA genomes. By applying comprehensive structure probing to diverse problems, we and others are showing that control of biological function mediated by RNA structure is ubiquitous across prokaryotic and eukaryotic organisms.Work over the past decade using SHAPE-based analyses has clarified key principles. First, the method of RNA structure probing matters. SHAPE-MaP, with its direct and one-step readout that probes nearly every nucleotide by reaction at the 2'-hydroxyl, gives a more detailed and accurate readout than alternatives. Second, comprehensive chemical probing is essential. Focusing on fragments of large RNAs or using meta-gene or statistical analyses to compensate for sparse data sets misses critical features and often yields structure models with poor predictive power. Finally, every RNA has its own internal structural personality. There are myriad ways in which RNA structure modulates sequence accessibility, protein binding, translation, splice-site choice, phase separation, and other fundamental biological processes. In essentially every instance where we have applied rigorous and quantitative SHAPE technologies to study RNA structure-function interrelationships, new insights regarding biological regulatory mechanisms have emerged. RNA elements with more complex higher-order structures appear more likely to contain high-information-content clefts and pockets that bind small molecules, broadly informing a vigorous field of RNA-targeted drug discovery.The broad implications of this collective work are twofold. First, it is long past time to abandon depiction of large RNAs as simple noodle-like or gently flowing molecules. Instead, we need to emphasize that nearly all RNAs are punctuated with distinctive internal structures, a subset of which modulate function in profound ways. Second, structure probing should be an integral component of any effort that seeks to understand the functional nexuses and biological roles of large 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.

1H, 13C and 15N chemical shift assignment of the stem-loop 5a from the 5'-UTR of SARS-CoV-2

The SARS-CoV-2 (SCoV-2) virus is the causative agent of the ongoing COVID-19 pandemic. It contains a positive sense single-stranded RNA genome and belongs to the genus of Betacoronaviruses. The 5'- and 3'-genomic ends of the 30 kb SCoV-2 genome are potential antiviral drug targets. Major parts of these sequences are highly conserved among Betacoronaviruses and contain cis-acting RNA elements that affect RNA translation and replication. The 31 nucleotide (nt) long highly conserved stem-loop 5a (SL5a) is located within the 5'-untranslated region (5'-UTR) important for viral replication. SL5a features a U-rich asymmetric bulge and is capped with a 5'-UUUCGU-3' hexaloop, which is also found in stem-loop 5b (SL5b). We herein report the extensive ¹H, ¹³C and ¹⁵N resonance assignment of SL5a as basis for in-depth structural studies by solution NMR spectroscopy.

Combinatorial optimization of mRNA structure, stability, and translation for RNA-based therapeutics

Therapeutic mRNAs and vaccines are being developed for a broad range of human diseases, including COVID-19. However, their optimization is hindered by mRNA instability and inefficient protein expression. Here, we describe design principles that overcome these barriers. We develop a new RNA sequencing-based platform called PERSIST-seq to systematically delineate in-cell mRNA stability, ribosome load, as well as in-solution stability of a library of diverse mRNAs. We find that, surprisingly, in-cell stability is a greater driver of protein output than high ribosome load. We further introduce a method called In-line-seq, applied to thousands of diverse RNAs, that reveals sequence and structure-based rules for mitigating hydrolytic degradation. Our findings show that "superfolder" mRNAs can be designed to improve both stability and expression that are further enhanced through pseudouridine nucleoside modification. Together, our study demonstrates simultaneous improvement of mRNA stability and protein expression and provides a computational-experimental platform for the enhancement of mRNA medicines.

In vivo structure of the Ty1 retrotransposon RNA genome

Long terminal repeat (LTR)-retrotransposons constitute a significant part of eukaryotic genomes and influence their function and evolution. Like other RNA viruses, LTR-retrotransposons efficiently utilize their RNA genome to interact with host cell machinery during replication. Here, we provide the first genome-wide RNA secondary structure model for a LTR-retrotransposon in living cells. Using SHAPE probing, we explore the secondary structure of the yeast Ty1 retrotransposon RNA genome in its native in vivo state and under defined in vitro conditions. Comparative analyses reveal the strong impact of the cellular environment on folding of Ty1 RNA. In vivo, Ty1 genome RNA is significantly less structured and more dynamic but retains specific well-structured regions harboring functional cis-acting sequences. Ribosomes participate in the unfolding and remodeling of Ty1 RNA, and inhibition of translation initiation stabilizes Ty1 RNA structure. Together, our findings support the dual role of Ty1 genomic RNA as a template for protein synthesis and reverse transcription. This study also contributes to understanding how a complex multifunctional RNA genome folds in vivo, and strengthens the need for studying RNA structure in its natural cellular context.

Evolutionary conservation of RNA sequence and structure

An RNA structure prediction from a single‐sequence RNA folding program is not evidence for an RNA whose structure is important for function. Random sequences have plausible and complex predicted structures not easily distinguishable from those of structural RNAs. How to tell when an RNA has a conserved structure is a question that requires looking at the evolutionary signature left by the conserved RNA. This question is important not just for long noncoding RNAs which usually lack an identified function, but also for RNA binding protein motifs which can be single stranded RNAs or structures. Here we review recent advances using sequence and structural analysis to determine when RNA structure is conserved or not. Although covariation measures assess structural RNA conservation, one must distinguish covariation due to RNA structure from covariation due to independent phylogenetic substitutions. We review a statistical test to measure false positives expected under the null hypothesis of phylogenetic covariation alone (specificity). We also review a complementary test that measures power, that is, expected covariation derived from sequence variation alone (sensitivity). Power in the absence of covariation signals the absence of a conserved RNA structure. We analyze artifacts that falsely identify conserved RNA structure such as the misuse of programs that do not assess significance, the use of inappropriate statistics confounded by signals other than covariation, or misalignments that induce spurious covariation. Among artifacts that obscure the signal of a conserved RNA structure, we discuss the inclusion of pseudogenes in alignments which increase power but destroy covariation.

This article is categorized under:

RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry

RNA Evolution and Genomics > Computational Analyses of RNA

RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution

Genome composition and genetic characterization of SARS-CoV-2

SARS-CoV-2 is a type of Betacoronaviruses responsible for COVID-19 pandemic disease, with more than 1.745 million fatalities globally as of December-2020. Genetically, it is considered the second largest genome of all RNA viruses with a 5' cap and 3' poly-A tail. Phylogenetic analyses of coronaviruses reveal that SARS-CoV-2 is genetically closely related to the Bat-SARS Like-Corona virus (Bat-SL-Cov) with 96% whole-genome identity. SARS-CoV-2 genome consists of 15 ORFs coded into 29 proteins. At the 5' terminal of the genome, we have ORF1ab and ORF1a, which encode the 1ab and 1a polypeptides that are proteolytically cleaved into 16 different nonstructural proteins (NSPs). The 3' terminal of the genome represents four structural (spike, envelope, matrix, and nucleocapsid) and nine accessory (3a, 3b, 6, 7a, 7b, 8b, 9a, 9b, and orf10) proteins. As the number of COVID-19 patients increases dramatically worldwide, there is an urgent need to find a quick and sensitive diagnostic tool for controlling the outbreak of SARS-CoV-2 in the community. Today, molecular testing methods utilizing viral genetic material (e.g., PCR) represent the crucial diagnostic tool for the SARS-CoV-2 virus despite its low sensitivity in the early stage of viral infection. This review summarizes the genome composition and genetic characterization of the SARS-CoV-2.

Structure and regulation of coronavirus genomes: state-of-the-art and novel insights from SARS-CoV-2 studies

Coronaviruses (CoV) are positive-sense single-stranded RNA viruses, harboring the largest viral RNA genomes known to date. Apart from the primary sequence encoding for all the viral proteins needed for the generation of new viral particles, certain regions of CoV genomes are known to fold into stable structures, controlling several aspects of CoV life cycle, from the regulation of the discontinuous transcription of subgenomic mRNAs, to the packaging of the genome into new virions. Here we review the current knowledge on CoV RNA structures, discussing it in light of the most recent discoveries made possible by analyses of the SARS-CoV-2 genome.

Genome-scale deconvolution of RNA structure ensembles

RNA structure heterogeneity is a major challenge when querying RNA structures with chemical probing. We introduce DRACO, an algorithm for the deconvolution of coexisting RNA conformations from mutational profiling experiments. Analysis of the SARS-CoV-2 genome using dimethyl sulfate mutational profiling with sequencing (DMS-MaPseq) and DRACO, identifies multiple regions that fold into two mutually exclusive conformations, including a conserved structural switch in the 3' untranslated region. This work may open the way to dissecting the heterogeneity of the RNA structurome.

Comprehensive in vivo secondary structure of the SARS-CoV-2 genome reveals novel regulatory motifs and mechanisms

Severe-acute-respiratory-syndrome-related coronavirus 2 (SARS-CoV-2) is the positive-sense RNA virus that causes coronavirus disease 2019 (COVID-19). The genome of SARS-CoV-2 is unique among viral RNAs in its vast potential to form RNA structures, yet as much as 97% of its 30 kilobases have not been structurally explored. Here, we apply a novel long amplicon strategy to determine the secondary structure of the SARS-CoV-2 RNA genome at single-nucleotide resolution in infected cells. Our in-depth structural analysis reveals networks of well-folded RNA structures throughout Orf1ab and reveals aspects of SARS-CoV-2 genome architecture that distinguish it from other RNA viruses. Evolutionary analysis shows that several features of the SARS-CoV-2 genomic structure are conserved across β-coronaviruses, and we pinpoint regions of well-folded RNA structure that merit downstream functional analysis. The native, secondary structure of SARS-CoV-2 presented here is a roadmap that will facilitate focused studies on the viral life cycle, facilitate primer design, and guide the identification of RNA drug targets against COVID-19.

Programmed -1 Ribosomal Frameshifting in coronaviruses: A therapeutic target

Human population growth, climate change, and globalization are accelerating the emergence of novel pathogenic viruses. In the past two decades alone, three such members of the coronavirus family have posed serious threats, spurring intense efforts to understand their biology as a way to identify targetable vulnerabilities. Coronaviruses use a programmed -1 ribosomal frameshift (-1 PRF) mechanism to direct synthesis of their replicase proteins. This is a critical switch in their replication program that can be therapeutically targeted. Here, we discuss how nearly half a century of research into -1 PRF have provided insight into the virological importance of -1 PRF, the molecular mechanisms that drive it, and approaches that can be used to manipulate it towards therapeutic outcomes with particular emphasis on SARS-CoV-2.

Amplification-free detection of SARS-CoV-2 with CRISPR-Cas13a and mobile phone microscopy

The December 2019 outbreak of a novel respiratory virus, SARS-CoV-2, has become an ongoing global pandemic due in part to the challenge of identifying symptomatic, asymptomatic, and pre-symptomatic carriers of the virus. CRISPR diagnostics can augment gold-standard PCR-based testing if they can be made rapid, portable, and accurate. Here, we report the development of an amplification-free CRISPR-Cas13a assay for direct detection of SARS-CoV-2 from nasal swab RNA that can be read with a mobile phone microscope. The assay achieved ∼100 copies/μL sensitivity in under 30 min of measurement time and accurately detected pre-extracted RNA from a set of positive clinical samples in under 5 min. We combined crRNAs targeting SARS-CoV-2 RNA to improve sensitivity and specificity and directly quantified viral load using enzyme kinetics. Integrated with a reader device based on a mobile phone, this assay has the potential to enable rapid, low-cost, point-of-care screening for SARS-CoV-2.

The RNA Architecture of the SARS-CoV-2 3'-Untranslated Region

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the current COVID-19 pandemic. The 3' untranslated region (UTR) of this β-CoV contains essential cis-acting RNA elements for the viral genome transcription and replication. These elements include an equilibrium between an extended bulged stem-loop (BSL) and a pseudoknot. The existence of such an equilibrium is supported by reverse genetic studies and phylogenetic covariation analysis and is further proposed as a molecular switch essential for the control of the viral RNA polymerase binding. Here, we report the SARS-CoV-2 3' UTR structures in cells that transcribe the viral UTRs harbored in a minigene plasmid and isolated infectious virions using a chemical probing technique, namely dimethyl sulfate (DMS)-mutational profiling with sequencing (MaPseq). Interestingly, the putative pseudoknotted conformation was not observed, indicating that its abundance in our systems is low in the absence of the viral nonstructural proteins (nsps). Similarly, our results also suggest that another functional cis-acting element, the three-helix junction, cannot stably form. The overall architectures of the viral 3' UTRs in the infectious virions and the minigene-transfected cells are almost identical.

Amilorides inhibit SARS-CoV-2 replication in vitro by targeting RNA structures

The SARS-CoV-2 pandemic, and the likelihood of future coronavirus pandemics, has rendered our understanding of coronavirus biology more essential than ever. Small molecule chemical probes offer to both reveal novel aspects of virus replication and to serve as leads for antiviral therapeutic development. The RNA-biased amiloride scaffold was recently tuned to target a viral RNA structure critical for translation in enterovirus 71, ultimately uncovering a novel mechanism to modulate positive-sense RNA viral translation and replication. Analysis of CoV RNA genomes reveal many conserved RNA structures in the 5'-UTR and proximal region critical for viral translation and replication, including several containing bulge-like secondary structures suitable for small molecule targeting. Following phylogenetic conservation analysis of this region, we screened an amiloride-based small molecule library against a less virulent human coronavirus, OC43, to identify lead ligands. Amilorides inhibited OC43 replication as seen in viral plaque assays. Select amilorides also potently inhibited replication competent SARS-CoV-2 as evident in the decreased levels of cell free virions in cell culture supernatants of treated cells. Reporter screens confirmed the importance of RNA structures in the 5'-end of the viral genome for small molecule activity. Finally, NMR chemical shift perturbation studies of the first six stem loops of the 5'-end revealed specific amiloride interactions with stem loops 4, 5a, and 6, all of which contain bulge like structures and were predicted to be strongly bound by the lead amilorides in retrospective docking studies. Taken together, the use of multiple orthogonal approaches allowed us to identify the first small molecules aimed at targeting RNA structures within the 5'-UTR and proximal region of the CoV genome. These molecules will serve as chemical probes to further understand CoV RNA biology and can pave the way for the development of specific CoV RNA-targeted antivirals.

Amplification-free detection of SARS-CoV-2 with CRISPR-Cas13a and mobile phone microscopy

The December 2019 outbreak of a novel respiratory virus, SARS-CoV-2, has become an ongoing global pandemic due in part to the challenge of identifying symptomatic, asymptomatic, and pre-symptomatic carriers of the virus. CRISPR diagnostics can augment gold-standard PCR-based testing if they can be made rapid, portable, and accurate. Here, we report the development of an amplification-free CRISPR-Cas13a assay for direct detection of SARS-CoV-2 from nasal swab RNA that can be read with a mobile phone microscope. The assay achieved ∼100 copies/μL sensitivity in under 30 min of measurement time and accurately detected pre-extracted RNA from a set of positive clinical samples in under 5 min. We combined crRNAs targeting SARS-CoV-2 RNA to improve sensitivity and specificity and directly quantified viral load using enzyme kinetics. Integrated with a reader device based on a mobile phone, this assay has the potential to enable rapid, low-cost, point-of-care screening for SARS-CoV-2.

Computational search for potential COVID-19 drugs from FDAapproved drugs and small molecules of natural origin identifies several anti-virals and plant products

The world is currently facing the COVID-19 pandemic, for which mild symptoms include fever and dry cough. In severe cases, it could lead to pneumonia and ultimately death in some instances. Moreover, the causative pathogen is highly contagious and there are no drugs or vaccines for it yet. The pathogen, SARS-CoV-2, is one of the human coronaviruses which was identified to infect humans first in December 2019. SARS-CoV-2 shares evolutionary relationship to other highly pathogenic viruses such as Severe Acute Respiratory Syndrome (SARS) and Middle East respiratory syndrome (MERS). We have exploited this similarity to model a target non-structural protein, NSP1, since it is implicated in the regulation of host gene expression by the virus and hijacking of host machinery. We next interrogated the capacity to repurpose around 2300 FDA-approved drugs and more than 3,00,000 small molecules of natural origin towards drug identification through virtual screening and molecular dynamics. Interestingly, we observed simple molecules like lactose, previously known anti-virals and few secondary metabolites of plants as promising hits. These herbal plants are already practiced in Ayurveda over centuries to treat respiratory problems and inflammation. Disclaimer: we would not like to recommend uptake of these small molecules for suspect COVID patients until it is approved by competent national or international authorities.