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

Double-stranded RNA drives SARS-CoV-2 nucleocapsid protein to undergo phase separation at specific temperatures

Nucleocapsid protein (N-protein) is required for multiple steps in betacoronaviruses replication. SARS-CoV-2-N-protein condenses with specific viral RNAs at particular temperatures making it a powerful model for deciphering RNA sequence specificity in condensates. We identify two separate and distinct double-stranded, RNA motifs (dsRNA stickers) that promote N-protein condensation. These dsRNA stickers are separately recognized by N-protein's two RNA binding domains (RBDs). RBD1 prefers structured RNA with sequences like the transcription-regulatory sequence (TRS). RBD2 prefers long stretches of dsRNA, independent of sequence. Thus, the two N-protein RBDs interact with distinct dsRNA stickers, and these interactions impart specific droplet physical properties that could support varied viral functions. Specifically, we find that addition of dsRNA lowers the condensation temperature dependent on RBD2 interactions and tunes translational repression. In contrast RBD1 sites are sequences critical for sub-genomic (sg) RNA generation and promote gRNA compression. The density of RBD1 binding motifs in proximity to TRS-L/B sequences is associated with levels of sub-genomic RNA generation. The switch to packaging is likely mediated by RBD1 interactions which generate particles that recapitulate the packaging unit of the virion. Thus, SARS-CoV-2 can achieve biochemical complexity, performing multiple functions in the same cytoplasm, with minimal protein components based on utilizing multiple distinct RNA motifs that control N-protein interactions.

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.

Geneticin shows selective antiviral activity against SARS-CoV-2 by interfering with programmed -1 ribosomal frameshifting

SARS-CoV-2 is currently causing an unprecedented pandemic. While vaccines are massively deployed, we still lack effective large-scale antiviral therapies. In the quest for antivirals targeting conserved structures, we focused on molecules able to bind viral RNA secondary structures. Aminoglycosides are a class of antibiotics known to interact with the ribosomal RNA of both prokaryotes and eukaryotes and have previously been shown to exert antiviral activities by interacting with viral RNA. Here we show that the aminoglycoside geneticin is endowed with antiviral activity against all tested variants of SARS-CoV-2, in different cell lines and in a respiratory tissue model at non-toxic concentrations. The mechanism of action is an early inhibition of RNA replication and protein expression related to a decrease in the efficiency of the -1 programmed ribosomal frameshift (PRF) signal of SARS-CoV-2. Using in silico modelling, we have identified a potential binding site of geneticin in the pseudoknot of frameshift RNA motif. Moreover, we have selected, through virtual screening, additional RNA binding compounds, interacting with the same site with increased potency.

Lessons Learned and Yet-to-Be Learned on the Importance of RNA Structure in SARS-CoV-2 Replication

SARS-CoV-2, the etiological agent responsible for the COVID-19 pandemic, is a member of the virus family Coronaviridae, known for relatively extensive (~30-kb) RNA genomes that not only encode for numerous proteins but are also capable of forming elaborate structures. As highlighted in this review, these structures perform critical functions in various steps of the viral life cycle, ultimately impacting pathogenesis and transmissibility. We examine these elements in the context of coronavirus evolutionary history and future directions for curbing the spread of SARS-CoV-2 and other potential human coronaviruses. While we focus on structures supported by a variety of biochemical, biophysical, and/or computational methods, we also touch here on recent evidence for novel structures in both protein-coding and noncoding regions of the genome, including an assessment of the potential role for RNA structure in the controversial finding of SARS-CoV-2 integration in “long COVID” patients. This review aims to serve as a consolidation of previous works on coronavirus and more recent investigation of SARS-CoV-2, emphasizing the need for improved understanding of the role of RNA structure in the evolution and adaptation of these human viruses.

Antisense oligonucleotide: A promising therapeutic option to beat COVID-19

The COVID-19 crisis and the development of the first approved mRNA vaccine have highlighted the power of RNA-based therapeutic strategies for the development of new medicines. Aside from RNA-vaccines, antisense oligonucleotides (ASOs) represent a new and very promising class of RNA-targeted therapy. Few drugs have already received approval from the Food and Drug Administration. Here, we underscored why and how ASOs hold the potential to change the therapeutic landscape to beat SARS-CoV-2 viral infections. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Small Molecule-RNA Interactions.

Putative Host-Derived Insertions in the Genomes of Circulating SARS-CoV-2 Variants

Insertions in the SARS-CoV-2 genome have the potential to drive viral evolution, but the source of the insertions is often unknown. Recent proposals have suggested that human RNAs could be a source of some insertions, but the small size of many insertions makes this difficult to confirm. Through an analysis of available direct RNA sequencing data from SARS-CoV-2-infected cells, we show that viral-host chimeric RNAs are formed through what are likely stochastic RNA-dependent RNA polymerase template-switching events. Through an analysis of the publicly available GISAID SARS-CoV-2 genome collection, we identified two genomic insertions in circulating SARS-CoV-2 variants that are identical to regions of the human 18S and 28S rRNAs. These results provide direct evidence of the formation of viral-host chimeric sequences and the integration of host genetic material into the SARS-CoV-2 genome, highlighting the potential importance of host-derived insertions in viral evolution. IMPORTANCE Throughout the COVID-19 pandemic, the sequencing of SARS-CoV-2 genomes has revealed the presence of insertions in multiple globally circulating lineages of SARS-CoV-2, including the Omicron variant. The human genome has been suggested to be the source of some of the larger insertions, but evidence for this kind of event occurring is still lacking. Here, we leverage direct RNA sequencing data and SARS-CoV-2 genomes to show that host-viral chimeric RNAs are generated in infected cells and two large genomic insertions have likely been formed through the incorporation of host rRNA fragments into the SARS-CoV-2 genome. These host-derived insertions may increase the genetic diversity of SARS-CoV-2 and expand its strategies to acquire genetic material, potentially enhancing its adaptability, virulence, and spread.

RNAvigator: A Pipeline to Identify Candidates for Functional RNA Structure Elements

Identifying structural elements in long and complex RNAs, such as long non-coding and RNA viruses, can shed light on the functionality and mechanisms of such RNAs. Here we present RNAvigator, a tool able to identify elements of structural importance by using experimental SHAPE data or SHAPE-like predictions in conjunction with stability and entropy assessments. RNAvigator recognizes regions that are the most stable, unambiguous, and structured on RNA molecules, and thus potentially functional. When relying on predictions, RNAvigator uses the CROSS algorithm, a neural network trained on experimental data that achieved an AUC of 0.74 on hepatitis C virus SHAPE-MaP data and which was able to improve the predictive power of Superfold. By using RNAvigator, we can identify known functional regions on the complete hepatitis C virus genome, including the regulatory regions CRE and IRES, and the 3’ UTR of dengue virus, a region known for the presence of structural elements essential for its replication, and functional regions of long non-coding RNAs such as XIST and HOTAIR. We envision that RNAvigator will be a useful tool for studying long and complex RNA molecules using known chemical probing data or, if they are not available, by employing predicted profiles.

Molecular Virology of SARS-CoV-2 and Related Coronaviruses

Coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The global COVID-19 pandemic continues to threaten the lives of hundreds of millions of people, with a severe negative impact on the global economy. Although several COVID-19 vaccines are currently being administered, none of them is 100% effective. Moreover, SARS-CoV-2 variants remain an important worldwide public health issue. Hence, the accelerated development of efficacious antiviral agents is urgently needed. Coronavirus depends on various host cell factors for replication. An ongoing research objective is the identification of host factors that could be exploited as targets for drugs and compounds effective against SARS-CoV-2. In the present review, we discuss the molecular mechanisms of SARS-CoV-2 and related coronaviruses, focusing on the host factors or pathways involved in SARS-CoV-2 replication that have been identified by genome-wide CRISPR screening.

Insertion-and-Deletion Mutations between the Genomes of SARS-CoV, SARS-CoV-2, and Bat Coronavirus RaTG13

The evolutional process of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) development remains inconclusive. This study compared the genome sequences of severe acute respiratory syndrome coronavirus (SARS-CoV), bat coronavirus RaTG13, and SARS-CoV-2. In total, the genomes of SARS-CoV-2 and RaTG13 were 77.9% and 77.7% identical to the genome of SARS-CoV, respectively. A total of 3.6% (1,068 bases) of the SARS-CoV-2 genome was derived from insertion and/or deletion (indel) mutations, and 18.6% (5,548 bases) was from point mutations from the genome of SARS-CoV. At least 35 indel sites were confirmed in the genome of SARS-CoV-2, in which 17 were with ≥10 consecutive bases long. Ten of these relatively long indels were located in the spike (S) gene, five in nonstructural protein 3 (Nsp3) gene of open reading frame (ORF) 1a, and one in ORF8 and noncoding region. Seventeen (48.6%) of the 35 indels were based on insertion-and-deletion mutations with exchanged gene sequences of 7–325 consecutive bases. Almost the complete ORF8 gene was replaced by a single 325 consecutive base-long indel. The distribution of these indels was roughly in accordance with the distribution of the rate of point mutation rate around the indels. The genome sequence of SARS-CoV-2 was 96.0% identical to that of RaTG13. There was no long insertion-and-deletion mutation between the genomes of RaTG13 and SARS-CoV-2. The findings of the uneven distribution of multiple indels and the presence of multiple long insertion-and-deletion mutations with exchanged consecutive base sequences in the viral genome may provide insights into SARS-CoV-2 development.

IMPORTANCE The developmental mechanism of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) remains inconclusive. This study compared the base sequence one-by-one between severe acute respiratory syndrome coronavirus (SARS-CoV) or bat coronavirus RaTG13 and SARS-CoV-2. The genomes of SARS-CoV-2 and RaTG13 were 77.9% and 77.7% identical to the genome of SARS-CoV, respectively. Seventeen of the 35 sites with insertion and/or deletion mutations between SARS-CoV-2 and SARS-CoV were based on insertion-and-deletion mutations with the replacement of 7–325 consecutive bases. Most of these long insertion-and-deletion sites were concentrated in the nonstructural protein 3 (Nsp3) gene of open reading frame (ORF) 1a, S1 domain of the spike protein, and ORF8 genes. Such long insertion-and-deletion mutations were not observed between the genomes of RaTG13 and SARS-CoV-2. The presence of multiple long insertion-and-deletion mutations in the genome of SARS-CoV-2 and their uneven distributions may provide further insights into the development of the virus.

Reconstitution of the SARS-CoV-2 ribonucleosome provides insights into genomic RNA packaging and regulation by phosphorylation

The nucleocapsid (N) protein of coronaviruses is responsible for compaction of the ∼30-kb RNA genome in the ∼100-nm virion. Cryo-electron tomography suggests that each virion contains 35-40 viral ribonucleoprotein (vRNP) complexes, or ribonucleosomes, arrayed along the genome. There is, however, little mechanistic understanding of the vRNP complex. Here, we show that N protein, when combined with viral RNA fragments in vitro, forms cylindrical 15-nm particles similar to the vRNP structures observed within coronavirus virions. These vRNPs form in the presence of stem-loop-containing RNA and depend on regions of N protein that promote protein-RNA and protein-protein interactions. Phosphorylation of N protein in its disordered serine/arginine (SR) region weakens these interactions and disrupts vRNP assembly. We propose that unmodified N binds stem-loop-rich regions in genomic RNA to form compact vRNP complexes within the nucleocapsid, while phosphorylated N maintains uncompacted viral RNA to promote the protein's transcriptional function.

The Emergence of SARS-CoV-2 Variants With a Lower Antibody Response: A Genomic and Clinical Perspective

The emergence of several novel SARS-CoV-2 variants regarded as variants of concern (VOCs) has exacerbated pathogenic and immunologic prominences, as well as reduced diagnostic sensitivity due to phenotype modification-capable mutations. Furthermore, latent and more virulent strains that have arisen as a result of unique mutations with increased evolutionary potential represent a threat to vaccine effectiveness in terms of incoming and existing variants. As a result, resisting natural immunity, which leads to higher reinfection rates, and avoiding vaccination-induced immunization, which leads to a lack of vaccine effectiveness, has become a crucial problem for public health around the world. This study attempts to review the genomic variation and pandemic impact of emerging variations of concern based on clinical characteristics management and immunization effectiveness. The goal of this study is to gain a better understanding of the link between genome level polymorphism, clinical symptom manifestation, and current vaccination in the instance of VOCs.

The highly conserved stem-loop II motif is important for the lifecycle of astroviruses but dispensable for SARS-CoV-2

The stem-loop II motif (s2m) is an RNA element present in viruses from divergent viral families, including astroviruses and coronaviruses, but its functional significance is unknown. We created deletions or substitutions of the s2m in astrovirus VA1 (VA1), classic human astrovirus 1 (HAstV1) and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). For VA1, recombinant virus could not be rescued upon partial deletion of the s2m or substitutions of G-C base pairs. Compensatory substitutions that restored the G-C base-pair enabled recovery of VA1. For HAstV1, a partial deletion of the s2m resulted in decreased viral titers compared to wild-type virus, and reduced activity in a replicon system. In contrast, deletion or mutation of the SARS-CoV-2 s2m had no effect on the ability to rescue the virus, growth in vitro , or growth in Syrian hamsters. Our study demonstrates the importance of the s2m is virus-dependent.

Putative host-derived insertions in the genomes of circulating SARS-CoV-2 variants

Insertions in the SARS-CoV-2 genome have the potential to drive viral evolution, but the source of the insertions is often unknown. Recent proposals have suggested that human RNAs could be a source of some insertions, but the small size of many insertions makes this difficult to confirm. Through an analysis of available direct RNA sequencing data from SARS-CoV-2 infected cells, we show that viral-host chimeric RNAs are formed through what are likely stochastic RNA-dependent RNA polymerase template switching events. Through an analysis of the publicly available GISAID SARS-CoV-2 genome collection, we identified two genomic insertions in circulating SARS-CoV-2 variants that are identical to regions of the human 18S and 28S rRNAs. These results provide direct evidence of the formation of viral-host chimeric sequences and the integration of host genetic material into the SARS-CoV-2 genome, highlighting the potential importance of host-derived insertions in viral evolution.

IMPORTANCE

Throughout the COVID-19 pandemic, the sequencing of SARS-CoV-2 genomes has revealed the presence of insertions in multiple globally circulating lineages of SARS-CoV-2, including the Omicron variant. The human genome has been suggested to be the source of some of the larger insertions, but evidence for this kind of event occurring is still lacking. Here, we leverage direct RNA sequencing data and SARS-CoV-2 genomes to show host-viral chimeric RNAs are generated in infected cells and two large genomic insertions have likely been formed through the incorporation of host rRNA fragments into the SARS-CoV-2 genome. These host-derived insertions may increase the genetic diversity of SARS-CoV-2 and expand its strategies to acquire genetic materials, potentially enhancing its adaptability, virulence, and spread.

The In Vivo and In Vitro Architecture of the Hepatitis C Virus RNA Genome Uncovers Functional RNA Secondary and Tertiary Structures

Hepatitis C virus (HCV) is a positive-strand RNA virus that remains one of the main contributors to chronic liver disease worldwide. Studies over the last 30 years have demonstrated that HCV contains a highly structured RNA genome and many of these structures play essential roles in the HCV life cycle. Despite the importance of riboregulation in this virus, most of the HCV RNA genome remains functionally unstudied. Here, we report a complete secondary structure map of the HCV RNA genome in vivo, which was studied in parallel with the secondary structure of the same RNA obtained in vitro. Our results show that HCV is folded extensively in the cellular context. By performing comprehensive structural analyses on both in vivo data and in vitro data, we identify compact and conserved secondary and tertiary structures throughout the genome. Genetic and evolutionary functional analyses demonstrate that many of these elements play important roles in the virus life cycle. In addition to providing a comprehensive map of RNA structures and riboregulatory elements in HCV, this work provides a resource for future studies aimed at identifying therapeutic targets and conducting further mechanistic studies on this important human pathogen. IMPORTANCE HCV has one of the most highly structured RNA genomes studied to date, and it is a valuable model system for studying the role of RNA structure in protein-coding genes. While previous studies have identified individual cases of regulatory RNA structures within the HCV genome, the full-length structure of the HCV genome has not been determined in vivo. Here, we present the complete secondary structure map of HCV determined both in cells and from corresponding transcripts generated in vitro. In addition to providing a comprehensive atlas of functional secondary structural elements throughout the genomic RNA, we identified a novel set of tertiary interactions and demonstrated their functional importance. In terms of broader implications, the pipeline developed in this study can be applied to other long RNAs, such as long noncoding RNAs. In addition, the RNA structural motifs characterized in this study broaden the repertoire of known riboregulatory elements.

Nucleic Acids as Biotools at the Interface between Chemistry and Nanomedicine in the COVID-19 Era

The recent development of mRNA vaccines against the SARS-CoV-2 infection has turned the spotlight on the potential of nucleic acids as innovative prophylactic agents and as diagnostic and therapeutic tools. Until now, their use has been severely limited by their reduced half-life in the biological environment and the difficulties related to their transport to target cells. These limiting aspects can now be overcome by resorting to chemical modifications in the drug and using appropriate nanocarriers, respectively. Oligonucleotides can interact with complementary sequences of nucleic acid targets, forming stable complexes and determining their loss of function. An alternative strategy uses nucleic acid aptamers that, like the antibodies, bind to specific proteins to modulate their activity. In this review, the authors will examine the recent literature on nucleic acids-based strategies in the COVID-19 era, focusing the attention on their applications for the prophylaxis of COVID-19, but also on antisense- and aptamer-based strategies directed to the diagnosis and therapy of the coronavirus pandemic.

1H, 13C and 15N chemical shift assignment of the stem-loops 5b + c from the 5'-UTR of SARS-CoV-2

The ongoing pandemic of the respiratory disease COVID-19 is caused by the SARS-CoV-2 (SCoV2) virus. SCoV2 is a member of the Betacoronavirus genus. The 30 kb positive sense, single stranded RNA genome of SCoV2 features 5'- and 3'-genomic ends that are highly conserved among Betacoronaviruses. These genomic ends contain structured cis-acting RNA elements, which are involved in the regulation of viral replication and translation. Structural information about these potential antiviral drug targets supports the development of novel classes of therapeutics against COVID-19. The highly conserved branched stem-loop 5 (SL5) found within the 5'-untranslated region (5'-UTR) consists of a basal stem and three stem-loops, namely SL5a, SL5b and SL5c. Both, SL5a and SL5b feature a 5'-UUUCGU-3' hexaloop that is also found among Alphacoronaviruses. Here, we report the extensive ¹H, ¹³C and ¹⁵N resonance assignment of the 37 nucleotides (nts) long sequence spanning SL5b and SL5c (SL5b + c), as basis for further in-depth structural studies by solution NMR spectroscopy.

RNA editing increases the nucleotide diversity of SARS-CoV-2 in human host cells

SARS-CoV-2 is a positive-sense, single-stranded RNA virus responsible for the COVID-19 pandemic. It remains unclear whether and to what extent the virus in human host cells undergoes RNA editing, a major RNA modification mechanism. Here we perform a robust bioinformatic analysis of metatranscriptomic data from multiple bronchoalveolar lavage fluid samples of COVID-19 patients, revealing an appreciable number of A-to-I RNA editing candidate sites in SARS-CoV-2. We confirm the enrichment of A-to-I RNA editing signals at these candidate sites through evaluating four characteristics specific to RNA editing: the inferred RNA editing sites exhibit (i) stronger ADAR1 binding affinity predicted by a deep-learning model built from ADAR1 CLIP-seq data, (ii) decreased editing levels in ADAR1-inhibited human lung cells, (iii) local clustering patterns, and (iv) higher RNA secondary structure propensity. Our results have critical implications in understanding the evolution of SARS-CoV-2 as well as in COVID-19 research, such as phylogenetic analysis and vaccine development.

The COVID-19 pandemic is caused by SARS-CoV-2, an RNA virus. In the cells of COVID-19 patients, SARS-CoV-2 interacts with human proteins and is potentially subjected to their enzymatic activities. Here we investigated whether human protein enzymes can change the nucleotide sequence of SARS-CoV-2, thereby leaving a unique molecular footprint. We developed a robust computational algorithm to analyze the sequence data of SARS-CoV-2 obtained from lung fluid samples of COVID-19 patients and found that the virus contains new nucleotide changes that are likely induced by ADAR1, a powerful human protein that can modify specific nucleotide positions in many human transcripts. We further confirmed that the characteristics of the nucleotide changes detected in SARS-CoV-2 are similar to those observed in the human genes. Thus, these ADAR1-induced nucleotide changes may represent an under-appreciated force that can affect the evolution of SARS-CoV-2. Our study helps researchers better understand the evolutionary trajectory of SARS-CoV-2.

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 an 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 highly structured "superfolder" mRNAs can be designed to improve both stability and expression with further enhancement 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.

Secondary structural ensembles of the SARS-CoV-2 RNA genome in infected cells

SARS-CoV-2 is a betacoronavirus with a single-stranded, positive-sense, 30-kilobase RNA genome responsible for the ongoing COVID-19 pandemic. Although population average structure models of the genome were recently reported, there is little experimental data on native structural ensembles, and most structures lack functional characterization. Here we report secondary structure heterogeneity of the entire SARS-CoV-2 genome in two lines of infected cells at single nucleotide resolution. Our results reveal alternative RNA conformations across the genome and at the critical frameshifting stimulation element (FSE) that are drastically different from prevailing population average models. Importantly, we find that this structural ensemble promotes frameshifting rates much higher than the canonical minimal FSE and similar to ribosome profiling studies. Our results highlight the value of studying RNA in its full length and cellular context. The genomic structures detailed here lay groundwork for coronavirus RNA biology and will guide the design of SARS-CoV-2 RNA-based therapeutics.

Coordinated regulation of interferon and inflammasome signaling pathways by SARS-CoV-2 proteins

Type I and III interferons (IFNs) and the nucleotide-binding domain (NBD) leucine-rich repeat (LRR)-containing receptor (NLR) family pyrin domain containing 3 (NLRP3) inflammasome play pivotal roles in the pathogenesis of SARS-CoV-2. While optimal IFN and inflammasome responses are essential for limiting SARS-CoV-2 infection, aberrant activation of these innate immune responses is associated with COVID-19 pathogenesis. In this review, we focus our discussion on recent findings on SARS-CoV-2-induced type I and III IFNs and NLRP3 inflammasome responses and the viral proteins regulating these mechanisms.

Ensemble cryo-EM reveals conformational states of the nsp13 helicase in the SARS-CoV-2 helicase replication-transcription complex

The SARS-CoV-2 nonstructural proteins coordinate genome replication and gene expression. Structural analyses revealed the basis for coupling of the essential nsp13 helicase with the RNA-dependent RNA polymerase (RdRp) where the holo-RdRp and RNA substrate (the replication-transcription complex or RTC) associated with two copies of nsp13 (nsp13₂-RTC). One copy of nsp13 interacts with the template-RNA in an opposing polarity to the RdRp and is envisaged to drive the RdRp backward on the RNA template (backtracking), prompting questions as to how the RdRp can efficiently synthesize RNA in the presence of nsp13. Here we use cryogenic-electron microscopy and molecular dynamics simulations to analyze the nsp13₂-RTC, revealing four distinct conformational states of the helicases. The results indicate a mechanism for the nsp13₂-RTC to turn backtracking on and off, using an allosteric mechanism to switch between RNA synthesis or backtracking in response to stimuli at the RdRp active site.

Dominant clade-featured SARS-CoV-2 co-occurring mutations reveal plausible epistasis: An in silico based hypothetical model

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has evolved into eight fundamental clades with four of these clades (G, GH, GR, and GV) globally prevalent in 2020. To explain plausible epistatic effects of the signature co-occurring mutations of these circulating clades on viral replication and transmission fitness, we proposed a hypothetical model using in silico approach. Molecular docking and dynamics analyses showed the higher infectiousness of a spike mutant through more favorable binding of G614 with the elastase-2. RdRp mutation p.P323L significantly increased genome-wide mutations (p < 0.0001), allowing for more flexible RdRp (mutated)-NSP8 interaction that may accelerate replication. Superior RNA stability and structural variation at NSP3:C241T might impact protein, RNA interactions, or both. Another silent 5'-UTR:C241T mutation might affect translational efficiency and viral packaging. These four G-clade-featured co-occurring mutations might increase viral replication. Sentinel GH-clade ORF3a:p.Q57H variants constricted the ion-channel through intertransmembrane-domain interaction of cysteine(C81)-histidine(H57). The GR-clade N:p.RG203-204KR would stabilize RNA interaction by a more flexible and hypo-phosphorylated SR-rich region. GV-clade viruses seemingly gained the evolutionary advantage of the confounding factors; nevertheless, N:p.A220V might modulate RNA binding with no phenotypic effect. Our hypothetical model needs further retrospective and prospective studies to understand detailed molecular events and their relationship to the fitness of SARS-CoV-2.

Secondary Structure of Subgenomic RNA M of SARS-CoV-2

SARS-CoV-2 belongs to the Coronavirinae family. Like other coronaviruses, SARS-CoV-2 is enveloped and possesses a positive-sense, single-stranded RNA genome of ~30 kb. Genomic RNA is used as the template for replication and transcription. During these processes, positive-sense genomic RNA (gRNA) and subgenomic RNAs (sgRNAs) are created. Several studies presented the importance of the genomic RNA secondary structure in SARS-CoV-2 replication. However, the structure of sgRNAs has remained largely unsolved so far. In this study, we probed the sgRNA M model of SARS-CoV-2 in vitro. The presented model molecule includes 5'UTR and a coding sequence of gene M. This is the first experimentally informed secondary structure model of sgRNA M, which presents features likely to be important in sgRNA M function. The knowledge of sgRNA M structure provides insights to better understand virus biology and could be used for designing new therapeutics.

Host cell-intrinsic innate immune recognition of SARS-CoV-2

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) emerged at the end of 2019 and caused the pandemic of coronavirus disease 2019 (COVID-19). Basic and clinical investigations indicate that severe forms of COVID-19 are due in part to dysregulated immune responses to virus infection. The innate immune system is the first line of host defense against most virus infections, with pathogen recognition receptors detecting SARS-CoV-2 RNA and protein components and initiating pro-inflammatory and antiviral responses. Notwithstanding this response, SARS-CoV-2 proteins evade, inhibit, and skew innate immune signaling early in infection. In this review, we highlight the components of cell-based recognition of SARS-CoV-2 infection and the mechanisms employed by the virus to modulate these innate immune host defense pathways.

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

SARS-CoV-2 produces two long viral protein precursors from one open reading frame using a highly conserved RNA pseudoknot that enhances programmed -1 ribosomal frameshifting. The 1.3 Å-resolution X-ray structure of the pseudoknot reveals three coaxially stacked helices buttressed by idiosyncratic base triples from loop residues. This structure represents a frameshift-stimulating state that must be deformed by the ribosome and exhibits base-triple-adjacent pockets that could be targeted by future small-molecule therapeutics.

Small Drugs, Huge Impact: The Extraordinary Impact of Antisense Oligonucleotides in Research and Drug Development

Antisense oligonucleotides (ASOs) are an increasingly represented class of drugs. These small sequences of nucleotides are designed to precisely target other oligonucleotides, usually RNA species, and are modified to protect them from degradation by nucleases. Their specificity is due to their sequence, so it is possible to target any RNA sequence that is already known. These molecules are very versatile and adaptable given that their sequence and chemistry can be custom manufactured. Based on the chemistry being used, their activity may significantly change and their effects on cell function and phenotypes can differ dramatically. While some will cause the target RNA to decay, others will only bind to the target and act as a steric blocker. Their incredible versatility is the key to manipulating several aspects of nucleic acid function as well as their process, and alter the transcriptome profile of a specific cell type or tissue. For example, they can be used to modify splicing or mask specific sites on a target. The entire design rather than just the sequence is essential to ensuring the specificity of the ASO to its target. Thus, it is vitally important to ensure that the complete process of drug design and testing is taken into account. ASOs' adaptability is a considerable advantage, and over the past decades has allowed multiple new drugs to be approved. This, in turn, has had a significant and positive impact on patient lives. Given current challenges presented by the COVID-19 pandemic, it is necessary to find new therapeutic strategies that would complement the vaccination efforts being used across the globe. ASOs may be a very powerful tool that can be used to target the virus RNA and provide a therapeutic paradigm. The proof of the efficacy of ASOs as an anti-viral agent is long-standing, yet no molecule currently has FDA approval. The emergence and widespread use of RNA vaccines during this health crisis might provide an ideal opportunity to develop the first anti-viral ASOs on the market. In this review, we describe the story of ASOs, the different characteristics of their chemistry, and how their characteristics translate into research and as a clinical tool.

Nuclear Magnetic Resonance Spectra and AMBER OL3 and ROC-RNA Simulations of UCUCGU Reveal Force Field Strengths and Weaknesses for Single-Stranded RNA

Single-stranded regions of RNA are important for folding of sequences into 3D structures and for design of therapeutics targeting RNA. Prediction of ensembles of 3D structures for single-stranded regions often involves classical mechanical approximations of interactions defined by quantum mechanical calculations on small model systems. Nuclear magnetic resonance (NMR) spectra and molecular dynamics (MD) simulations of short single strands provide tests for how well the approximations model many of the interactions. Here, the NMR spectra for UCUCGU at 2, 15, and 30 °C are compared to simulations with the AMBER force fields, OL3 and ROC-RNA. This is the first such comparison to an oligoribonucleotide containing an internal guanosine nucleotide (G). G is particularly interesting because of its many H-bonding groups, large dipole moment, and proclivity for both syn and anti conformations. Results reveal formation of a G amino to phosphate non-bridging oxygen H-bond. The results also demonstrate dramatic differences in details of the predicted structures. The variations emphasize the dependence of predictions on individual parameters and their balance with the rest of the force field. The NMR data can serve as a benchmark for future force fields.

Targeting a conserved structural element from the SARS-CoV-2 genome using l-DNA aptamers

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19, has resulted in an ongoing global pandemic that has already claimed more than 4 million lives. While most antiviral efforts have focused on essential SARS-CoV-2 proteins, RNA structural elements within the viral genome are also compelling targets. In this study, we identified high-affinity l-DNA aptamers against a SARS-CoV-2 stem-loop II-like motif (s2m), a highly conserved RNA structure with promising diagnostic and therapeutic potential. Optimized l-C1t and l-C3t aptamers bind selectively to s2m RNA with K d values in the nanomolar range, and are capable of differentiating the monomeric s2m stem-loop from the proposed homodimer duplex. This structure-specific mode of recognition also allows l-C1t and l-C3t to discriminate between s2m RNAs from SARS-CoV-2 and SARS-CoV-1, which differ by just two nucleotides. Finally, we show that l-C1t and l-C3t induce dramatic conformational changes in s2m structure upon binding, and thus, have the potential to block protein-s2m interactions. Overall, these results demonstrate the feasibility of targeting SARS-CoV-2 RNA using l-aptamers, which has important implications in the diagnosis and treatment of COVID-19. Moreover, the high affinity and selectivity of l-C1t and l-C3t, coupled with the intrinsic nuclease resistance of l-DNA, present an opportunity for generating new tools and probes for interrogating s2m function in SARS-CoV-2 and related viruses.

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

Conserved SARS-CoV-2 RNA regions of critical biological functions define excellent targets for anti-viral therapeutics against Covid-19 variants. One such region is the frameshifting element (FSE), responsible for correct translation of viral polyproteins. Here, we analyze molecular-dynamics motions of three FSE conformations, discovered by graph-theory analysis, and associated mutants designed by graph-based inverse folding: two distinct 3-stem H-type pseudoknots and a 3-way junction. We find that the prevalent H-type pseudoknot in literature adopts ring-like conformations, which in combination with 5' end threading could promote ribosomal pausing. An inherent shape switch from "L" to linear that may help trigger the frameshifting is suppressed in our designed mutant. The alternative conformation trajectories suggest a stable intermediate structure with mixed stem interactions of all three conformations, pointing to a possible transition pathway during ribosomal translation. These observations provide new insights into anti-viral strategies and frameshifting mechanisms.

Structures and functions of coronavirus replication-transcription complexes and their relevance for SARS-CoV-2 drug design

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has killed millions of people and continues to cause massive global upheaval. Coronaviruses are positive-strand RNA viruses with an unusually large genome of ~30 kb. They express an RNA-dependent RNA polymerase and a cohort of other replication enzymes and supporting factors to transcribe and replicate their genomes. The proteins performing these essential processes are prime antiviral drug targets, but drug discovery is hindered by our incomplete understanding of coronavirus RNA synthesis and processing. In infected cells, the RNA-dependent RNA polymerase must coordinate with other viral and host factors to produce both viral mRNAs and new genomes. Recent research aiming to decipher and contextualize the structures, functions and interplay of the subunits of the SARS-CoV-2 replication and transcription complex proteins has burgeoned. In this Review, we discuss recent advancements in our understanding of the molecular basis and complexity of the coronavirus RNA-synthesizing machinery. Specifically, we outline the mechanisms and regulation of RNA translation, replication and transcription. We also discuss the composition of the replication and transcription complexes and their suitability as targets for antiviral therapy.

Absolute quantitation of individual SARS-CoV-2 RNA molecules provides a new paradigm for infection dynamics and variant differences

Despite an unprecedented global research effort on SARS-CoV-2, early replication events remain poorly understood. Given the clinical importance of emergent viral variants with increased transmission, there is an urgent need to understand the early stages of viral replication and transcription. We used single-molecule fluorescence in situ hybridisation (smFISH) to quantify positive sense RNA genomes with 95% detection efficiency, while simultaneously visualising negative sense genomes, subgenomic RNAs, and viral proteins. Our absolute quantification of viral RNAs and replication factories revealed that SARS-CoV-2 genomic RNA is long-lived after entry, suggesting that it avoids degradation by cellular nucleases. Moreover, we observed that SARS-CoV-2 replication is highly variable between cells, with only a small cell population displaying high burden of viral RNA. Unexpectedly, the B.1.1.7 variant, first identified in the UK, exhibits significantly slower replication kinetics than the Victoria strain, suggesting a novel mechanism contributing to its higher transmissibility with important clinical implications.

LinearTurboFold: Linear-time global prediction of conserved structures for RNA homologs with applications to SARS-CoV-2

The constant emergence of COVID-19 variants reduces the effectiveness of existing vaccines and test kits. Therefore, it is critical to identify conserved structures in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes as potential targets for variant-proof diagnostics and therapeutics. However, the algorithms to predict these conserved structures, which simultaneously fold and align multiple RNA homologs, scale at best cubically with sequence length and are thus infeasible for coronaviruses, which possess the longest genomes (∼30,000 nt) among RNA viruses. As a result, existing efforts on modeling SARS-CoV-2 structures resort to single-sequence folding as well as local folding methods with short window sizes, which inevitably neglect long-range interactions that are crucial in RNA functions. Here we present LinearTurboFold, an efficient algorithm for folding RNA homologs that scales linearly with sequence length, enabling unprecedented global structural analysis on SARS-CoV-2. Surprisingly, on a group of SARS-CoV-2 and SARS-related genomes, LinearTurboFold's purely in silico prediction not only is close to experimentally guided models for local structures, but also goes far beyond them by capturing the end-to-end pairs between 5' and 3' untranslated regions (UTRs) (∼29,800 nt apart) that match perfectly with a purely experimental work. Furthermore, LinearTurboFold identifies undiscovered conserved structures and conserved accessible regions as potential targets for designing efficient and mutation-insensitive small-molecule drugs, antisense oligonucleotides, small interfering RNAs (siRNAs), CRISPR-Cas13 guide RNAs, and RT-PCR primers. LinearTurboFold is a general technique that can also be applied to other RNA viruses and full-length genome studies and will be a useful tool in fighting the current and future pandemics.

SARS-CoV-2 expresses a microRNA-like small RNA able to selectively repress host genes

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of coronavirus disease (COVID-19), continues to be a pressing health concern. In this study, we investigated the impact of SARS-CoV-2 infection on host microRNA (miRNA) populations in three human lung-derived cell lines, as well as in nasopharyngeal swabs from SARS-CoV-2-infected individuals. We did not detect any major and consistent differences in host miRNA levels after SARS-CoV-2 infection. However, we unexpectedly discovered a viral miRNA-like small RNA, named CoV2-miR-O7a (for SARS-CoV-2 miRNA-like ORF7a-derived small RNA). Its abundance ranges from low to moderate as compared to host miRNAs and it associates with Argonaute proteins-core components of the RNA interference pathway. We identify putative targets for CoV2-miR-O7a, including Basic Leucine Zipper ATF-Like Transcription Factor 2 (BATF2), which participates in interferon signaling. We demonstrate that CoV2-miR-O7a production relies on cellular machinery, yet is independent of Drosha protein, and is enhanced by the presence of a strong and evolutionarily conserved hairpin formed within the ORF7a sequence.

Inhibition of SARS-CoV-2 by Targeting Conserved Viral RNA Structures and Sequences

The ongoing COVID-19/Severe Acute Respiratory Syndrome CoV-2 (SARS-CoV-2) pandemic has become a significant threat to public health and has hugely impacted societies globally. Targeting conserved SARS-CoV-2 RNA structures and sequences essential for viral genome translation is a novel approach to inhibit viral infection and progression. This new pharmacological modality compasses two classes of RNA-targeting molecules: 1) synthetic small molecules that recognize secondary or tertiary RNA structures and 2) antisense oligonucleotides (ASOs) that recognize the RNA primary sequence. These molecules can also serve as a "bait" fragment in RNA degrading chimeras to eliminate the viral RNA genome. This new type of chimeric RNA degrader is recently named ribonuclease targeting chimera or RIBOTAC. This review paper summarizes the sequence conservation in SARS-CoV-2 and the current development of RNA-targeting molecules to combat this virus. These RNA-binding molecules will also serve as an emerging class of antiviral drug candidates that might pivot to address future viral outbreaks.

Highly conserved s2m element of SARS-CoV-2 dimerizes via a kissing complex and interacts with host miRNA-1307-3p

The ongoing COVID-19 pandemic highlights the necessity for a more fundamental understanding of the coronavirus life cycle. The causative agent of the disease, SARS-CoV-2, is being studied extensively from a structural standpoint in order to gain insight into key molecular mechanisms required for its survival. Contained within the untranslated regions of the SARS-CoV-2 genome are various conserved stem-loop elements that are believed to function in RNA replication, viral protein translation, and discontinuous transcription. While the majority of these regions are variable in sequence, a 41-nucleotide s2m element within the genome 3' untranslated region is highly conserved among coronaviruses and three other viral families. In this study, we demonstrate that the SARS-CoV-2 s2m element dimerizes by forming an intermediate homodimeric kissing complex structure that is subsequently converted to a thermodynamically stable duplex conformation. This process is aided by the viral nucleocapsid protein, potentially indicating a role in mediating genome dimerization. Furthermore, we demonstrate that the s2m element interacts with multiple copies of host cellular microRNA (miRNA) 1307-3p. Taken together, our results highlight the potential significance of the dimer structures formed by the s2m element in key biological processes and implicate the motif as a possible therapeutic drug target for COVID-19 and other coronavirus-related diseases.

Stepwise Evolution and Exceptional Conservation of ORF1a/b Overlap in Coronaviruses

The programmed frameshift element (PFE) rerouting translation from ORF1a to ORF1b is essential for the propagation of coronaviruses. The combination of genomic features that make up PFE-the overlap between the two reading frames, a slippery sequence, as well as an ensemble of complex secondary structure elements-places severe constraints on this region as most possible nucleotide substitution may disrupt one or more of these elements. The vast amount of SARS-CoV-2 sequencing data generated within the past year provides an opportunity to assess the evolutionary dynamics of PFE in great detail. Here, we performed a comparative analysis of all available coronaviral genomic data available to date. We show that the overlap between ORF1a and ORF1b evolved as a set of discrete 7, 16, 22, 25, and 31 nucleotide stretches with a well-defined phylogenetic specificity. We further examined sequencing data from over 1,500,000 complete genomes and 55,000 raw read data sets to demonstrate exceptional conservation and detect signatures of selection within the PFE region.

RNA helicases required for viral propagation in humans

RNA viruses cause many routine illnesses, such as the common cold and the flu. Recently, more deadly diseases have emerged from this family of viruses. The hepatitis C virus has had a devastating impact worldwide. Despite the cures developed in the U.S. and Europe, economically disadvantaged countries remain afflicted by HCV infection due to the high cost of these medications. More recently, COVID-19 has swept across the world, killing millions and disrupting economies and lifestyles; the virus responsible for this pandemic is a coronavirus. Our understanding of HCV and SARS CoV-2 replication is still in its infancy. Helicases play a critical role in the replication, transcription and translation of viruses. These key enzymes need extensive study not only as an essential player in the viral lifecycle, but also as targets for antiviral therapeutics. In this review, we highlight the current knowledge for RNA helicases of high importance to human health.

Template switching and duplications in SARS-CoV-2 genomes give rise to insertion variants that merit monitoring

The appearance of multiple new SARS-CoV-2 variants during the COVID-19 pandemic is a matter of grave concern. Some of these variants, such as B.1.617.2, B.1.1.7, and B.1.351, manifest higher infectivity and virulence than the earlier SARS-CoV-2 variants, with potential dramatic effects on the course of the pandemic. So far, analysis of new SARS-CoV-2 variants focused primarily on nucleotide substitutions and short deletions that are readily identifiable by comparison to consensus genome sequences. In contrast, insertions have largely escaped the attention of researchers although the furin site insert in the Spike (S) protein is thought to be a determinant of SARS-CoV-2 virulence. Here, we identify 346 unique inserts of different lengths in SARS-CoV-2 genomes and present evidence that these inserts reflect actual virus variance rather than sequencing artifacts. Two principal mechanisms appear to account for the inserts in the SARS-CoV-2 genomes, polymerase slippage and template switch that might be associated with the synthesis of subgenomic RNAs. At least three inserts in the N-terminal domain of the S protein are predicted to lead to escape from neutralizing antibodies, whereas other inserts might result in escape from T-cell immunity. Thus, inserts in the S protein can affect its antigenic properties and merit monitoring.

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

The SARS-CoV-2 pandemic, and the likelihood of future coronavirus pandemics, emphasized the urgent need for development of novel antivirals. Small-molecule chemical probes offer both to reveal aspects of virus replication and to serve as leads for antiviral therapeutic development. Here, we report on the identification of amiloride-based small molecules that potently inhibit OC43 and SARS-CoV-2 replication through targeting of conserved structured elements within the viral 5′-end. Nuclear magnetic resonance–based structural studies revealed specific amiloride interactions with stem loops containing bulge like structures and were predicted to be strongly bound by the lead amilorides in retrospective docking studies. Amilorides represent the first antiviral small molecules that target RNA structures within the 5′ untranslated regions and proximal region of the CoV genomes. 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.

LazySampling and LinearSampling: Fast Stochastic Sampling of RNA Secondary Structure with Applications to SARS-CoV-2

Many RNAs fold into multiple structures at equilibrium. The classical stochastic sampling algorithm can sample secondary structures according to their probabilities in the Boltzmann ensemble, and is widely used. However, this algorithm, consisting of a bottom-up partition function phase followed by a top-down sampling phase, suffers from three limitations: (a) the formulation and implementation of the sampling phase are unnecessarily complicated; (b) the sampling phase repeatedly recalculates many redundant recursions already done during the partition function phase; (c) the partition function runtime scales cubically with the sequence length. These issues prevent stochastic sampling from being used for very long RNAs such as the full genomes of SARS-CoV-2. To address these problems, we first adopt a hypergraph framework under which the sampling algorithm can be greatly simplified. We then present three sampling algorithms under this framework, among which the LazySampling algorithm is the fastest by eliminating redundant work in the sampling phase via on-demand caching. Based on LazySampling, we further replace the cubic-time partition function by a linear-time approximate one, and derive LinearSampling, an end-to-end linear-time sampling algorithm that is orders of magnitude faster than the standard one. For instance, LinearSampling is 176Ã- faster (38.9s vs. 1.9h) than Vienna RNAsubopt on the full genome of Ebola virus (18,959 nt ). More importantly, LinearSampling is the first RNA structure sampling algorithm to scale up to the full-genome of SARS-CoV-2 without local window constraints, taking only 69.2 seconds on its reference sequence (29,903 nt ). The resulting sample correlates well with the experimentally-guided structures. On the SARS-CoV-2 genome, LinearSampling finds 23 regions of 15 nt with high accessibilities, which are potential targets for COVID-19 diagnostics and drug design. See code: https://github.com/LinearFold/LinearSampling.

LinearTurboFold: Linear-Time Global Prediction of Conserved Structures for RNA Homologs with Applications to SARS-CoV-2

The constant emergence of COVID-19 variants reduces the effectiveness of existing vaccines and test kits. Therefore, it is critical to identify conserved structures in SARS-CoV-2 genomes as potential targets for variant-proof diagnostics and therapeutics. However, the algorithms to predict these conserved structures, which simultaneously fold and align multiple RNA homologs, scale at best cubically with sequence length, and are thus infeasible for coronaviruses, which possess the longest genomes (∼30,000 nt ) among RNA viruses. As a result, existing efforts on modeling SARS-CoV-2 structures resort to single sequence folding as well as local folding methods with short window sizes, which inevitably neglect long-range interactions that are crucial in RNA functions. Here we present LinearTurboFold, an efficient algorithm for folding RNA homologs that scales linearly with sequence length, enabling unprecedented global structural analysis on SARS-CoV-2. Surprisingly, on a group of SARS-CoV-2 and SARS-related genomes, LinearTurbo-Fold's purely in silico prediction not only is close to experimentally-guided models for local structures, but also goes far beyond them by capturing the end-to-end pairs between 5' and 3' UTRs (∼29,800 nt apart) that match perfectly with a purely experimental work. Furthermore, LinearTurboFold identifies novel conserved structures and conserved accessible regions as potential targets for designing efficient and mutation-insensitive small-molecule drugs, antisense oligonucleotides, siRNAs, CRISPR-Cas13 guide RNAs and RT-PCR primers. LinearTurboFold is a general technique that can also be applied to other RNA viruses and full-length genome studies, and will be a useful tool in fighting the current and future pandemics.

SIGNIFICANCE STATEMENT

Conserved RNA structures are critical for designing diagnostic and therapeutic tools for many diseases including COVID-19. However, existing algorithms are much too slow to model the global structures of full-length RNA viral genomes. We present LinearTurboFold, a linear-time algorithm that is orders of magnitude faster, making it the first method to simultaneously fold and align whole genomes of SARS-CoV-2 variants, the longest known RNA virus (∼30 kilobases). Our work enables unprecedented global structural analysis and captures long-range interactions that are out of reach for existing algorithms but crucial for RNA functions. LinearTurboFold is a general technique for full-length genome studies and can help fight the current and future pandemics.

CoV-er all the bases: Structural perspectives of SARS-CoV-2 RNA synthesis

The ongoing Covid-19 pandemic has spurred research in the biology of the nidovirus severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Much focus has been on the viral RNA synthesis machinery due to its fundamental role in viral propagation. The central and essential enzyme of the RNA synthesis process, the RNA-dependent RNA polymerase (RdRp), functions in conjunction with a coterie of viral-encoded enzymes that mediate crucial nucleic acid transactions. Some of these enzymes share common features with other RNA viruses, while others play roles unique to nidoviruses or CoVs. The RdRps are proven targets for viral pathogens, and many of the other nucleic acid processing enzymes are promising targets. The purpose of this review is to summarize recent advances in our understanding of the mechanisms of RNA synthesis in CoVs. By reflecting on these studies, we hope to emphasize the remaining gaps in our knowledge. The recent onslaught of structural information related to SARS-CoV-2 RNA synthesis, in combination with previous structural, genetic and biochemical studies, have vastly improved our understanding of how CoVs replicate and process their genomic RNA. Structural biology not only provides a blueprint for understanding the function of the enzymes and cofactors in molecular detail, but also provides a basis for drug design and optimization. The concerted efforts of researchers around the world, in combination with the renewed urgency toward understanding this deadly family of viruses, may eventually yield new and improved antivirals that provide relief to the current global devastation.

Evaluating RNA Structural Flexibility: Viruses Lead the Way

Our understanding of RNA structure has lagged behind that of proteins and most other biological polymers, largely because of its ability to adopt multiple, and often very different, functional conformations within a single molecule. Flexibility and multifunctionality appear to be its hallmarks. Conventional biochemical and biophysical techniques all have limitations in solving RNA structure and to address this in recent years we have seen the emergence of a wide diversity of techniques applied to RNA structural analysis and an accompanying appreciation of its ubiquity and versatility. Viral RNA is a particularly productive area to study in that this economy of function within a single molecule admirably suits the minimalist lifestyle of viruses. Here, we review the major techniques that are being used to elucidate RNA conformational flexibility and exemplify how the structure and function are, as in all biology, tightly linked.

Targeting SARS-CoV-2 Variants with Nucleic Acid Therapeutic Nanoparticle Conjugates

The emergence of SARS-CoV-2 variants is cause for concern, because these may become resistant to current vaccines and antiviral drugs in development. Current drugs target viral proteins, resulting in a critical need for RNA-targeted nanomedicines. To address this, a comparative analysis of SARS-CoV-2 variants was performed. Several highly conserved sites were identified, of which the most noteworthy is a partial homopurine palindrome site with >99% conservation within the coding region. This sequence was compared among recently emerged, highly infectious SARS-CoV-2 variants. Conservation of the site was maintained among these emerging variants, further contributing to its potential as a regulatory target site for SARS-CoV-2. RNAfold was used to predict the structures of the highly conserved sites, with some resulting structures being common among coronaviridae. An RNA-level regulatory map of the conserved regions of SARS-CoV-2 was produced based on the predicted structures, with each representing potential target sites for antisense oligonucleotides, triplex-forming oligomers, and aptamers. Additionally, homopurine/homopyrimidine sequences within the viral genome were identified. These sequences also demonstrate appropriate target sites for antisense oligonucleotides and triplex-forming oligonucleotides. An experimental strategy to investigate these is summarized along with potential nanoparticle types for delivery, and the advantages and disadvantages of each are discussed.

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.

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.

Comprehensive mapping of SARS-CoV-2 interactions in vivo reveals functional virus-host interactions

SARS-CoV-2 is a major threat to global health. Here, we investigate the RNA structure and RNA-RNA interactions of wildtype (WT) and a mutant (Δ382) SARS-CoV-2 in cells using Illumina and Nanopore platforms. We identify twelve potentially functional structural elements within the SARS-CoV-2 genome, observe that subgenomic RNAs can form different structures, and that WT and Δ382 virus genomes fold differently. Proximity ligation sequencing identify hundreds of RNA-RNA interactions within the virus genome and between the virus and host RNAs. SARS-CoV-2 genome binds strongly to mitochondrial and small nucleolar RNAs and is extensively 2'-O-methylated. 2'-O-methylation sites are enriched in viral untranslated regions, associated with increased virus pair-wise interactions, and are decreased in host mRNAs upon virus infection, suggesting that the virus sequesters methylation machinery from host RNAs towards its genome. These studies deepen our understanding of the molecular and cellular basis of SARS-CoV-2 pathogenicity and provide a platform for targeted therapy.

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.