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

Wastewater sample storage for physicochemical and microbiological analysis

Wastewater-based epidemiology (WBE) is a crucial tool for health and environmental monitoring, providing real-time data on public health indicators by analysis of sewage samples. Ensuring the integrity of these samples from collection to analysis is paramount. This study investigates the effects of different cold-storage conditions on the integrity of wastewater samples, focusing on both microbiological markers (such as extractable nucleic acids, SARS-CoV-2, and crAssphage) and physicochemical parameters (including ammonium, orthophosphate, pH, conductivity, and turbidity). Composite samples from the combined raw wastewater influent from five wastewater treatment works in South Wales, UK, were stored at 4°C, -20°C, and -80°C, and subjected to up to six freeze-thaw cycles over one year. The study found significant effects of storage temperature on the preservation of certain WBE markers, with the best yield most frequently seen in samples stored at -80°C. However, the majority of WBE markers showed no significant difference between storage at -80°C or at 4°C, demonstrating that it may not always be necessary to archive wastewater samples at ultra-low temperatures, thus reducing CO2 emissions and laboratory energy costs. These findings underscore the importance of optimized storage conditions to maintain sample integrity, while ensuring accurate and reliable WBE data for public health and environmental monitoring.

Bat RNA viruses employ viral RHIMs orchestrating species-specific cell death programs linked to Z-RNA sensing and ZBP1-RIPK3 signaling

RHIM is a protein motif facilitating the assembly of large signaling complexes triggering regulated cell death. A few DNA viruses employ viral RHIMs mimicking host RHIMs and counteract cell death by interacting with host RHIM-proteins to alleviate antiviral defenses. Whether RNA viruses operate such viral RHIMs remains unknown. Here, we identified viral RHIMs in Nsp13 of SARS-CoV-2 and other bat RNA viruses, providing the basis for bats as the hosts for their evolution. Nsp13 promoted viral RHIM and RNA-binding channel-dependent cell death. However, Nsp13 viral RHIM is more critical for human cell death than in bat-derived Tb1 Lu cells, suggesting species-specific regulation. Nsp13 showed RHIM-dependent interactions with ZBP1 and RIPK3, forming large complexes and promoting ZBP1-RIPK3 signaling-mediated cell death. Intriguingly, the SARS-CoV-2 genome consisted of Z-RNA-forming segments promoting Nsp13-dependent cell death. Our findings reveal the functional viral RHIMs of bat-originated RNA viruses regulating host cell death associated with ZBP1-RIPK3 signaling, indicating possible mechanisms of cellular damage and cytokine storm in bat-originated RNA virus infections.

Analysis of natural structures and chemical mapping data reveals local stability compensation in RNA

RNA molecules adopt complex structures that perform essential biological functions across all forms of life, making them promising candidates for therapeutic applications. However, our ability to design new RNA structures remains limited by an incomplete understanding of their folding principles. While global metrics such as the minimum free energy are widely used, they are at odds with naturally occurring structures and incompatible with established design rules. Here, we introduce local stability compensation (LSC), a principle that RNA folding is governed by the local balance between destabilizing loops and their stabilizing adjacent stems, challenging the focus on global energetic optimization. Analysis of over 100,000 RNA structures revealed that LSC signatures are particularly pronounced in bulges and their adjacent stems, with distinct patterns across different RNA families that align with their biological functions. To validate LSC experimentally, we systematically analyzed thousands of RNA variants using DMS chemical mapping. Our results demonstrate that stem reactivity correlates strongly with LSC (R2 = 0.458 for hairpin loops) and that structural perturbations affect folding primarily within ~6 nucleotides from the loop. These findings establish LSC as a fundamental principle that could enhance the rational design of functional RNAs.

C→U transition biases in SARS-CoV-2: still rampant 4 years from the start of the COVID-19 pandemic

The evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in the pandemic and post-pandemic periods has been characterized by rapid adaptive changes that confer immune escape and enhanced human-to-human transmissibility. Sequence change is additionally marked by an excess number of C→U transitions suggested as being due to host-mediated genome editing. To investigate how these influence the evolutionary trajectory of SARS-CoV-2, 2,000 high-quality, coding complete genome sequences of SARS-CoV-2 variants collected pre-September 2020 and from each subsequently appearing alpha, delta, BA.1, BA.2, BA.5, XBB, EG, HK, and JN.1 lineages were downloaded from NCBI Virus in April 2024. C→U transitions were the most common substitution during the diversification of SARS-CoV-2 lineages over the 4-year observation period. A net loss of C bases and accumulation of U's occurred at a constant rate of approximately 0.2%-0.25%/decade. C→U transitions occurred in over a quarter of all sites with a C (26.5%; range 20.0%-37.2%) around five times more than observed for the other transitions (5.3%-6.8%). In contrast to an approximately random distribution of other transitions across the genome, most C→U substitutions occurred at statistically preferred sites in each lineage. However, only the most C→U polymorphic sites showed evidence for a preferred 5'U context previously associated with APOBEC 3A editing. There was a similarly weak preference for unpaired bases suggesting much less stringent targeting of RNA than mediated by A3 deaminases in DNA editing. Future functional studies are required to determine editing preferences, impacts on replication fitness in vivo of SARS-CoV-2 and other RNA viruses, and impact on host tropism.

IMPORTANCE

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in the pandemic and post-pandemic periods has shown a remarkable capacity to adapt and evade human immune responses and increase its human-to-human transmissibility. The genome of SARS-CoV-2 is also increasingly scarred by the effects of multiple C→U mutations from host genome editing as a cellular defense mechanism akin to restriction factors for retroviruses. Through the analysis of large data sets of SARS-CoV-2 isolate sequences collected throughout the pandemic period and beyond, we show that C→U transitions have driven a base compositional change over time amounting to a net loss of C bases and accumulation of U's at a rate of approximately 0.2%-0.25%/decade. Most C→U substitutions occurred in the absence of the preferred upstream-base context or targeting of unpaired RNA bases previously associated with the host RNA editing protein, APOBEC 3A. The analyses provide a series of testable hypotheses that can be experimentally investigated in the future.

Conformational penalties: New insights into nucleic acid recognition

The energy cost accompanying changes in the structures of nucleic acids when they bind partner molecules is a significant but underappreciated thermodynamic contribution to binding affinity and specificity. This review highlights recent advances in measuring conformational penalties and determining their contribution to the recognition, folding, and regulatory activities of nucleic acids. Notable progress includes methods for measuring and structurally characterizing lowly populated conformational states, obtaining ensemble information in high throughput, for large macromolecular assemblies, and in complex cellular environments. Additionally, quantitative and predictive thermodynamic models have been developed that relate conformational penalties to nucleic acid-protein association and cellular activity. These studies underscore the crucial role of conformational penalties in nucleic acid recognition.

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

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

A quantitative framework for structural interpretation of DMS reactivity

Dimethyl sulfate (DMS) chemical mapping is widely used for probing RNA structure, with low reactivity interpreted as Watson-Crick (WC) base pairs and high reactivity as unpaired nucleotides. Despite its widespread use, a quantitative understanding of how DMS reactivity relates to specific RNA 3D structural features remains incomplete. To address this gap, we systematically analyzed DMS reactivity patterns with a massive library of 7,500 RNA constructs containing two-way junctions with known 3D structures. Our results reveal that DMS reactivity exists on a continuous spectrum rather than discrete high and low bins. Approximately 10% overlap in reactivity between WC and non-WC nucleotides demonstrates that simple thresholds cannot accurately determine base-pairing status. In flanking WC pairs, DMS reactivity correlates with base stacking strength and junction dynamics. For non-WC nucleotides, increased hydrogen bonding and decreased solvent accessibility led to WC-like DMS protection. Most significantly, we discover that DMS reactivity in non-canonical pairs strongly correlates with atomic distances and base pair geometry, enabling discrimination between different 3D conformations. These quantitative relationships establish novel metrics for evaluating RNA structural models and provide a new framework for incorporating DMS reactivity patterns into structure prediction algorithms.

Analysis of the abundance and diversity of RNA secondary structure elements in RNA viruses using the RNAsselem Python package

Recent advancements in experimental and computational methods for RNA secondary structure detection have revealed the crucial role of RNA structural elements in diverse molecular processes within living cells. It has been demonstrated that the secondary structure of the entire viral genome is often responsible for performing crucial functions in the viral life cycle and also influences virus evolution. To investigate the role of viral RNA secondary structure, alongside experimental techniques, the use of bioinformatics tools is important for analyzing various secondary structure patterns, including hairpin loops, internal loops, multifurcations, external loops, bulges, stems, and pseudoknots. Here, we have introduced a Python package for analyzing RNA secondary structure elements in viral genomes, which includes the recognition of common secondary structure patterns, the generation of descriptive statistics for these structural elements, and the provision of their basic properties. We applied the developed package to analyze the secondary structures of complete viral genomes collected from the literature, aiming to gain insights into viral function and evolution. Both the package and the collection of secondary structures of viral genomes are available at http://github.com/KazanovLab/RNAsselem .

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

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

Direct Zeptomole Detection of RNA Biomarkers by Ultrabright Fluorescent Nanoparticles on Magnetic Beads

Nucleic acids are important biomarkers in cancer and viral diseases. However, their ultralow concentration in biological/clinical samples makes direct target detection challenging, because it leads to slow hybridization kinetics with the probe and its insufficient signal-to-noise ratio. Therefore, RNA target detection is done by molecular (target) amplification, notably by RT-PCR, which is a tedious multistep method that includes nucleic acid extraction and reverse transcription. Here, a direct method based on ultrabright dye-loaded polymeric nanoparticles in a sandwich-like hybridization assay with magnetic beads is reported. The ultrabright DNA-functionalized nanoparticle, equivalent to ≈10 000 strongly emissive rhodamine dyes, is hybridized with the magnetic bead to the RNA target, providing the signal amplification for the detection. This concept (magneto-fluorescent sandwich) enables high-throughput detection of DNA and RNA sequences of varied lengths from 48 to 1362 nt with the limit of detection down to 0.3 fm using a plate reader (15 zeptomoles), among the best reported for optical sandwich assays. Moreover, it allows semi-quantitative detection of SARS-CoV-2 viral RNA directly in clinical samples without a dedicated RNA extraction step. The developed technology, combining ultrabright nanoparticles with magnetic beads, addresses fundamental challenges in RNA detection; it is expected to accelerate molecular diagnostics of diseases.

Structural basis for the participation of the SARS-CoV-2 nucleocapsid protein in the template switch mechanism and genomic RNA reorganization

The COVID-19 pandemic has resulted in a significant toll of deaths worldwide, exceeding seven million individuals, prompting intensive research efforts aimed at elucidating the molecular mechanisms underlying the pathogenesis of SARS-CoV-2 infection. Despite the rapid development of effective vaccines and therapeutic interventions, COVID-19 remains a threat to humans due to the emergence of novel variants and largely unknown long-term consequences. Among the viral proteins, the nucleocapsid protein (N) stands out as the most conserved and abundant, playing the primary role in nucleocapsid assembly and genome packaging. The N protein is promiscuous for the recognition of RNA, yet it can perform specific functions. Here, we discuss the structural basis of specificity, which is directly linked to its regulatory role. Notably, the RNA chaperone activity of N is central to its multiple roles throughout the viral life cycle. This activity encompasses double-stranded RNA (dsRNA) annealing and melting and facilitates template switching, enabling discontinuous transcription. N also promotes the formation of membrane-less compartments through liquid-liquid phase separation, thereby facilitating the congregation of the replication and transcription complex. Considering the information available regarding the catalytic activities and binding signatures of the N protein-RNA interaction, this review focuses on the regulatory role of the SARS-CoV-2 N protein. We emphasize the participation of the N protein in discontinuous transcription, template switching, and RNA chaperone activity, including double-stranded RNA melting and annealing activities.

Viral RNA Interactome: The Ultimate Researcher's Guide to RNA-Protein Interactions

RNA molecules in the cell are bound by a multitude of RNA-binding proteins (RBPs) with a variety of regulatory consequences. Often, interactions with these RNA-binding proteins are facilitated by the complex secondary and tertiary structures of RNA molecules. Viral RNAs especially are known to be heavily structured and interact with many RBPs, with roles including genome packaging, immune evasion, enhancing replication and transcription, and increasing translation efficiency. As such, the RNA-protein interactome represents a critical facet of the viral replication cycle. Characterization of these interactions is necessary for the development of novel therapeutics targeted at the disruption of essential replication cycle events. In this review, we aim to summarize the various roles of RNA structures in shaping the RNA-protein interactome, the regulatory roles of these interactions, as well as up-to-date methods developed for the characterization of the interactome and directions for novel, RNA-directed therapeutics.

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

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

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

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

Making target sites in large structured RNAs accessible to RNA-cleaving DNAzymes through hybridization with synthetic DNA oligonucleotides

The 10-23 DNAzyme is one of the most active DNA-based enzymes, and in theory, can be designed to target any purine-pyrimidine junction within an RNA sequence for cleavage. However, purine-pyrimidine junctions within a large, structured RNA (lsRNA) molecule of biological origin are not always accessible to 10-23, negating its general utility as an RNA-cutting molecular scissor. Herein, we report a generalizable strategy that allows 10-23 to access any purine-pyrimidine junction within an lsRNA. Using three large SARS-CoV-2 mRNA sequences of 566, 584 and 831 nucleotides in length as model systems, we show that the use of antisense DNA oligonucleotides (ASOs) that target the upstream and downstream regions flanking the cleavage site can restore the activity (kobs) of previously poorly active 10-23 DNAzyme systems by up to 2000-fold. We corroborated these findings mechanistically using in-line probing to demonstrate that ASOs reduced 10-23 DNAzyme target site structure within the lsRNA substrates. This approach represents a simple, efficient, cost-effective, and generalizable way to improve the accessibility of 10-23 to a chosen target site within an lsRNA molecule, especially where direct access to the genomic RNA target is necessary.

Untangling the pseudoknots of SARS-CoV-2: Insights into structural heterogeneity and plasticity

Since the onset of the COVID-19 pandemic, one productive area of research has focused on the intricate two- and three-dimensional structures taken on by SARS-CoV-2's RNA genome. These structures control essential viral processes, making them tempting targets for therapeutic intervention. This review focuses on two such structured regions, the frameshift stimulation element (FSE), which controls the translation of viral protein, and the 3' untranslated region (3' UTR), which is thought to regulate genome replication. For the FSE, we discuss its canonical pseudoknot's threaded and unthreaded topologies, as well as the diversity of competing two-dimensional structures formed by local and long-distance base pairing. For the 3' UTR, we review the evidence both for and against the formation of its replication-enabling pseudoknot.

Structural and biophysical dissection of RNA conformational ensembles

RNA's ability to form and interconvert between multiple secondary and tertiary structures is critical to its functional versatility and the traditional view of RNA structures as static entities has shifted towards understanding them as dynamic conformational ensembles. In this review we discuss RNA structural ensembles and their dynamics, highlighting the concept of conformational energy landscapes as a unifying framework for understanding RNA processes such as folding, misfolding, conformational changes, and complex formation. Ongoing advancements in cryo-electron microscopy and chemical probing techniques are significantly enhancing our ability to investigate multiple structures adopted by conformationally dynamic RNAs, while traditional methods such as nuclear magnetic resonance spectroscopy continue to play a crucial role in providing high-resolution, quantitative spatial and temporal information. We discuss how these methods, when used synergistically, can provide a comprehensive understanding of RNA conformational ensembles, offering new insights into their regulatory functions.

Context-dependent structure formation of hairpin motifs in bacteriophage MS2 genomic RNA

Many functions of ribonucleic acid (RNA) rely on its ability to assume specific sequence-structure motifs. Packaging signals found in certain RNA viruses are one such prominent example of functional RNA motifs. These signals are short hairpin loops that interact with coat proteins and drive viral self-assembly. As they are found in different positions along the much longer genomic RNA, the formation of their correct structure occurs as a part of a larger context. Any changes to this context can consequently lead to changes in the structure of the motifs themselves. In fact, previous studies have shown that structure and function of RNA motifs can be highly context sensitive to the flanking sequence surrounding them. However, in what ways different flanking sequences influence the structure of an RNA motif they surround has yet to be studied in detail. We focus on a hairpin-rich region of the RNA genome of bacteriophage MS2-a well-studied RNA virus with a wide potential for use in biotechnology-and systematically examine context-dependent structural stability of 14 previously identified hairpin motifs, which include putative and confirmed packaging signals. Combining secondary and tertiary RNA structure prediction of the hairpin motifs placed in different contexts, ranging from the native genomic sequence to random RNA sequences and unstructured poly-U sequences, we determine different measures of motif structural stability. In this way, we show that while some motif structures can be stable in any context, others require specific context provided by the genome. Our results demonstrate the importance of context in RNA structure formation and how changes in the flanking sequence of an RNA motif sometimes lead to drastic changes in its structure. Structural stability of a motif in different contexts could provide additional insights into its functionality as well as assist in determining whether it remains functional when intentionally placed in other contexts.

Visualizing RNA structure ensembles by single-molecule correlated chemical probing

RNA molecules fold to form complex internal structures. Many of these RNA structures populate ensembles with rheostat-like properties, with each state having a distinct function. Until recently, analysis of RNA structures, especially within cells, was limited to modeling either a single averaged structure or computationally-modeled ensembles. These approaches obscure the intrinsic heterogeneity of many structured RNAs. Single-molecule correlated chemical probing (smCCP) strategies are now making it possible to measure and deconvolute RNA structure ensembles based on efficiently executed chemical probing experiments. Here, we provide an overview of fundamental single-molecule probing principles, review current ensemble deconvolution strategies, and discuss recent applications to diverse biological systems. smCCP is enabling a revolution in understanding how the plasticity of RNA structure is exploited in biological systems to respond to stimuli and alter gene function. The energetics of RNA ensembles are often subtle and a subset can likely be targeted to modulate disease-associated biological processes.

Prediction of the effects of the top 10 synonymous mutations from 26645 SARS-CoV-2 genomes of early pandemic phase

Background

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) had led to a global pandemic since December 2019. SARS-CoV-2 is a single-stranded RNA virus, which mutates at a higher rate. Multiple works had been done to study nonsynonymous mutations, which change protein sequences. However, there is little study on the effects of SARS-CoV-2 synonymous mutations, which may affect viral fitness. This study aims to predict the effect of synonymous mutations on the SARS-CoV-2 genome.

Methods

A total of 26645 SARS-CoV-2 genomic sequences retrieved from Global Initiative on Sharing all Influenza Data (GISAID) database were aligned using MAFFT. Then, the mutations and their respective frequency were identified. Multiple RNA secondary structures prediction tools, namely RNAfold, IPknot++ and MXfold2 were applied to predict the effect of the mutations on RNA secondary structure and their base pair probabilities was estimated using MutaRNA. Relative synonymous codon usage (RSCU) analysis was also performed to measure the codon usage bias (CUB) of SARS-CoV-2.

Results

A total of 150 synonymous mutations were identified. The synonymous mutation identified with the highest frequency is C3037U mutation in the nsp3 of ORF1a. Of these top 10 highest frequency synonymous mutations, C913U, C3037U, U16176C and C18877U mutants show pronounced changes between wild type and mutant in all 3 RNA secondary structure prediction tools, suggesting these mutations may have some biological impact on viral fitness. These four mutations show changes in base pair probabilities. All mutations except U16176C change the codon to a more preferred codon, which may result in higher translation efficiency.

Conclusion

Synonymous mutations in SARS-CoV-2 genome may affect RNA secondary structure, changing base pair probabilities and possibly resulting in a higher translation rate. However, lab experiments are required to validate the results obtained from prediction analysis.

Requirement of the N-terminal region of nonstructural protein 1 in cis for SARS-CoV-2 defective RNA replication

SARS-CoV-2 belongs to the family Coronaviridae and carries a single-stranded positive-sense RNA genome. During coronavirus (CoV) replication, defective or defective interfering RNAs that lack a large portion of the genome often emerge. These defective RNAs typically carry the necessary RNA elements that are required for replication and packaging. We identified the minimum requirement of the 5' proximal region necessary for viral RNA replication by using artificially generated SARS-CoV-2 minigenomes. The minigenomes consist of the 5'-proximal region, an open reading frame (ORF) that encodes a fusion protein consisting of the N-terminal of viral NSP1 and a reporter gene, and the 3' untranslated region of the SARS-CoV-2 genome. We used a modified SARS-CoV-2 variant to support replication of the minigenomes. A minigenome carrying the 5' proximal 634 nucleotides replicated, whereas those carrying shorter than 634 nucleotides did not, demonstrating that the entire 265 nt-long 5' untranslated region and N-terminal portion of the NSP1 coding region are required for the minigenome replication. Minigenome RNAs carrying a specific amino acid substitution or frame shift insertions in the partial NSP1 coding sequence abrogated minigenome replication. Introduction of synonymous mutations in the minigenome RNAs also affected the replication efficiency of the minigenomes. These data suggest that the expression of the N-terminal portion of NSP1 and the primary sequence of the 5' proximal 634 nucleotides are important for minigenome replication.IMPORTANCESARS-CoV-2, the causative agent of COVID-19, is highly transmissible and continues to have a significant impact on public health and the global economy. While several vaccines mitigate the severe consequences of SARS-CoV-2 infection, mutant viruses with reduced reactivity to current vaccines continue to emerge and circulate. This study aimed to identify the minimal 5' proximal region of SARS-CoV-2 genomic RNA required for SARS-CoV-2 defective RNA replication and investigate the importance of an ORF encoded in these defective RNAs. Identifying cis-acting replication signals of SARS-CoV-2 genomic RNA is critical for the development of antivirals that target these signals. Additionally, replication-competent defective RNAs can serve as therapeutic reagents to interfere with SARS-CoV-2 replication. Our findings provide valuable insights into the mechanisms of SARS-CoV-2 RNA replication and the development of reagents that suppress SARS-CoV-2 replication.

High-throughput determination of RNA tertiary contact thermodynamics by quantitative DMS chemical mapping

Structured RNAs often contain long-range tertiary contacts that are critical to their function. Despite the importance of tertiary contacts, methods to measure their thermodynamics are low throughput or require specialized instruments. Here, we introduce a new quantitative chemical mapping method (qMaPseq) to measure Mg2+-induced formation of tertiary contact thermodynamics in a high-throughput manner using standard biochemistry equipment. With qMaPseq, we measured the ΔG of 98 unique tetraloop/tetraloop receptor (TL/TLR) variants in a one-pot reaction. These results agree well with measurements from specialized instruments (R2= 0.64). Furthermore, the DMS reactivity of the TL directly correlates to the stability of the contact (R2= 0.68), the first direct evidence that a single DMS reactivity measurement reports on thermodynamics. Combined with structure prediction, DMS reactivity allowed the development of experimentally accurate 3D models of TLR mutants. These results demonstrate that qMaPseq is broadly accessible, high-throughput and directly links DMS reactivity to thermodynamics.

Telomerase RNA structural heterogeneity in living human cells detected by DMS-MaPseq

Telomerase is a specialized reverse transcriptase that uses an intrinsic RNA subunit as the template for telomeric DNA synthesis. Biogenesis of human telomerase requires its RNA subunit (hTR) to fold into a multi-domain architecture that includes the template-containing pseudoknot (t/PK) and the three-way junction (CR4/5). These two hTR domains bind the telomerase reverse transcriptase (hTERT) protein and are thus essential for telomerase catalytic activity. Here, we probe the structure of hTR in living cells using dimethyl sulfate mutational profiling with sequencing (DMS-MaPseq) and ensemble deconvolution analysis. Unexpectedly, approximately 15% of the steady state population of hTR has a CR4/5 conformation lacking features thought to be required for hTERT binding. The proportion of hTR CR4/5 that is folded into the primary functional conformation does not require hTERT expression and the fraction of hTR that assumes a misfolded CR4/5 domain is not refolded by overexpression of its hTERT binding partner. This result suggests a functional role for an RNA folding cofactor other than hTERT during telomerase biogenesis. Mutagenesis demonstrates that stabilization of the alternative CR4/5 conformation is detrimental to telomerase assembly and activity. Moreover, the alternative CR4/5 conformation is not found in telomerase RNP complexes purified from cells via an epitope tag on hTERT, supporting the hypothesis that only the major CR4/5 conformer is active. We propose that this misfolded portion of the cellular hTR pool is either slowly refolded or degraded. Thus, kinetic traps for RNA folding that have been so well-studied in vitro may also present barriers for assembly of ribonucleoprotein complexes in vivo.

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

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

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

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

Molecular Bases and Specificity behind the Activation of the Immune System OAS/RNAse L Pathway by Viral RNA

The first line of defense against invading pathogens usually relies on innate immune systems. In this context, the recognition of exogenous RNA structures is primordial to fight, notably, against RNA viruses. One of the most efficient immune response pathways is based on the sensing of RNA double helical motifs by the oligoadenylate synthase (OAS) proteins, which in turn triggers the activity of RNase L and, thus, cleaves cellular and viral RNA. In this contribution, by using long-range molecular dynamics simulations, complemented with enhanced sampling techniques, we elucidate the structural features leading to the activation of OAS by interaction with a model double-strand RNA oligomer mimicking a viral RNA. We characterize the allosteric regulation induced by the nucleic acid leading to the population of the active form of the protein. Furthermore, we also identify the free energy profile connected to the active vs. inactive conformational transitions in the presence and absence of RNA. Finally, the role of two RNA mutations, identified as able to downregulate OAS activation, in shaping the protein/nucleic acid interface and the conformational landscape of OAS is also analyzed.

Targeting sgRNA N secondary structure as a way of inhibiting SARS-CoV-2 replication

SARS-CoV-2 is a betacoronavirus that causes COVID-19, a global pandemic that has resulted in many infections, deaths, and socio-economic challenges. The virus has a large positive-sense, single-stranded RNA genome of ∼30 kb, which produces subgenomic RNAs (sgRNAs) through discontinuous transcription. The most abundant sgRNA is sgRNA N, which encodes the nucleocapsid (N) protein. In this study, we probed the secondary structure of sgRNA N and a shorter model without a 3' UTR in vitro, using the SHAPE (selective 2'-hydroxyl acylation analyzed by a primer extension) method and chemical mapping with dimethyl sulfate and 1-cyclohexyl-(2-morpholinoethyl) carbodiimide metho-p-toluene sulfonate. We revealed the secondary structure of sgRNA N and its shorter variant for the first time and compared them with the genomic RNA N structure. Based on the structural information, we designed gapmers, siRNAs and antisense oligonucleotides (ASOs) to target the N protein coding region of sgRNA N. We also generated eukaryotic expression vectors containing the complete sequence of sgRNA N and used them to screen for new SARS-CoV-2 gene N expression inhibitors. Our study provides novel insights into the structure and function of sgRNA N and potential therapeutic tools against SARS-CoV-2.

Examining the stability of viral RNA and DNA in wastewater: Effects of storage time, temperature, and freeze-thaw cycles

Wastewater-based epidemiology (WBE) has been demonstrably successful as a relatively unbiased tool for monitoring levels of SARS-CoV-2 virus circulating in communities during the COVID-19 pandemic. Accumulated biobanks of wastewater samples allow retrospective exploration of spatial and temporal trends for public health indicators such as chemicals, viruses, antimicrobial resistance genes, and the possible emergence of novel human or zoonotic pathogens. We investigated virus resilience to time, temperature, and freeze-thaw cycles, plus the optimal storage conditions to maintain the stability of genetic material (RNA/DNA) of viral +ssRNA (Envelope - E, Nucleocapsid - N and Spike protein - S genes of SARS-CoV-2), dsRNA (Phi6 phage) and circular dsDNA (crAssphage) in wastewater. Samples consisted of (i) processed and extracted wastewater samples, (ii) processed and extracted distilled water samples, and (iii) raw, unprocessed wastewater samples. Samples were stored at -80 °C, -20 °C, 4 °C, or 20 °C for 10 days, going through up to 10 freeze-thaw cycles (once per day). Sample stability was measured using reverse transcription quantitative PCR, quantitative PCR, automated electrophoresis, and short-read whole genome sequencing. Exploring different areas of the SARS-CoV-2 genome demonstrated that the S gene in processed and extracted samples showed greater sensitivity to freeze-thaw cycles than the E or N genes. Investigating surrogate and normalisation viruses showed that Phi6 remains a stable comparison for SARS-CoV-2 in a laboratory setting and crAssphage was relatively resilient to temperature variation. Recovery of SARS-CoV-2 in raw unprocessed samples was significantly greater when stored at 4 °C, which was supported by the sequencing data for all viruses - both time and freeze-thaw cycles negatively impacted sequencing metrics. Historical extracts stored at -80 °C that were re-quantified 12, 14 and 16 months after original quantification showed no major changes. This study highlights the importance of the fast processing and extraction of wastewater samples, following which viruses are relatively robust to storage at a range of temperatures.

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

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

The 5'-terminal stem-loop RNA element of SARS-CoV-2 features highly dynamic structural elements that are sensitive to differences in cellular pH

We present the nuclear magnetic resonance spectroscopy (NMR) solution structure of the 5'-terminal stem loop 5_SL1 (SL1) of the SARS-CoV-2 genome. SL1 contains two A-form helical elements and two regions with non-canonical structure, namely an apical pyrimidine-rich loop and an asymmetric internal loop with one and two nucleotides at the 5'- and 3'-terminal part of the sequence, respectively. The conformational ensemble representing the averaged solution structure of SL1 was validated using NMR residual dipolar coupling (RDC) and small-angle X-ray scattering (SAXS) data. We show that the internal loop is the major binding site for fragments of low molecular weight. This internal loop of SL1 can be stabilized by an A12-C28 interaction that promotes the transient formation of an A+•C base pair. As a consequence, the pKa of the internal loop adenosine A12 is shifted to 5.8, compared to a pKa of 3.63 of free adenosine. Furthermore, applying a recently developed pH-differential mutational profiling (PD-MaP) approach, we not only recapitulated our NMR findings of SL1 but also unveiled multiple sites potentially sensitive to pH across the 5'-UTR of SARS-CoV-2.

The human mitochondrial mRNA structurome reveals mechanisms of gene expression

The human mitochondrial genome encodes crucial oxidative phosphorylation system proteins, pivotal for aerobic energy transduction. They are translated from nine monocistronic and two bicistronic transcripts whose native structures remain unexplored, posing a gap in understanding mitochondrial gene expression. In this work, we devised the mitochondrial dimethyl sulfate mutational profiling with sequencing (mitoDMS-MaPseq) method and applied detection of RNA folding ensembles using expectation-maximization (DREEM) clustering to unravel the native mitochondrial messenger RNA (mt-mRNA) structurome in wild-type (WT) and leucine-rich pentatricopeptide repeat-containing protein (LRPPRC)-deficient cells. Our findings elucidate LRPPRC's role as a holdase contributing to maintaining mt-mRNA folding and efficient translation. mt-mRNA structural insights in WT mitochondria, coupled with metabolic labeling, unveil potential mRNA-programmed translational pausing and a distinct programmed ribosomal frameshifting mechanism. Our data define a critical layer of mitochondrial gene expression regulation. These mt-mRNA folding maps provide a reference for studying mt-mRNA structures in diverse physiological and pathological contexts.

Peptide nucleic acids (PNAs) control function of SARS-CoV-2 frameshifting stimulatory element trough PNA-RNA-PNA triplex formation

The highly structured nature of the SARS-CoV-2 genome provides many promising antiviral drug targets. One particularly promising target is a cis-acting RNA pseudoknot found within a critical region called the frameshifting stimulatory element (FSE). In this study, peptide nucleic acids (PNAs) binding to stem 2 of FSE RNA inhibited protein translation and frameshifting, as measured by a cell-free dual luciferase assay, more effectively than PNAs binding to stem 1, stem 3, or the slippery site. Surprisingly, simple antisense PNAs were stronger disruptors of frameshifting than PNA tail-clamps, despite higher thermal stability of the PNA-RNA-PNA triplexes formed by the latter. Another unexpected result was a strong and sequence non-specific enhancement of frameshifting inhibition when using a cationic triplex-forming PNA in conjunction with an antisense PNA targeting key regions of the frameshifting element. Our results illustrate both the potential and the challenges of using antisense PNAs to target highly structured RNAs, such as SARS-CoV-2 pseudoknots. While triplex forming PNAs, including PNA tail-clamps, are emerging as promising ligands for RNA recognition, the binding affinity enhancements when using cationic modifications in triplex-forming PNAs must be carefully balanced to avoid loss of sequence specificity in complex biological systems.

An intricate balancing act: Upstream and downstream frameshift co-regulatory elements

Targeting ribosomal frameshifting has emerged as a potential therapeutic intervention strategy against Covid-19. During ribosomal translation, a fraction of elongating ribosomes slips by one base in the 5' direction and enters a new reading frame for viral protein synthesis. Any interference with this process profoundly affects viral replication and propagation. For Covid-19, two RNA sites associated with ribosomal frameshifting for SARS-CoV-2 are positioned on the 5' and 3' of the frameshifting residues. Although much attention has been on the 3' frameshift element (FSE), the 5' stem-loop (attenuator hairpin, AH) can play a role. The formation of AH has been suggested to occur as refolding of the 3' RNA structure is triggered by ribosomal unwinding. However, the attenuation activity and the relationship between the two regions are unknown. To gain more insight into these two related viral RNAs and to further enrich our understanding of ribosomal frameshifting for SARS-CoV-2, we explore the RNA folding of both 5' and 3' regions associated with frameshifting. Using our graph-theory-based modeling tools to represent RNA secondary structures, "RAG" (RNA- As-Graphs), and conformational landscapes to analyze length-dependent conformational distributions, we show that AH coexists with the 3-stem pseudoknot of the 3' FSE (graph 3_6 in our dual graph notation) and alternative pseudoknot (graph 3_3) but less likely with other 3' FSE alternative folds (such as 3-way junction 3_5). This is because an alternative length-dependent Stem 1 (AS1) can disrupt the FSE pseudoknots and trigger other folds. In addition, we design four mutants for long lengths that stabilize or disrupt AH, AS1 or FSE pseudoknot to illustrate the deduced AH/AS1 roles and favor the 3_5, 3_6 or stem-loop. These mutants further show how a strengthened pseudoknot can result from a weakened AS1, while a dominant stem-loop occurs with a strengthened AS1. These structural and mutational insights into both ends of the FSE in SARS-CoV-2 advance our understanding of the SARS-CoV-2 frameshifting mechanism by suggesting a sequence of length-dependent folds, which in turn define potential therapeutic intervention techniques involving both elements. Our work also highlights the complexity of viral landscapes with length-dependent folds, and challenges in analyzing these multiple conformations.

Structure of an internal loop motif with three consecutive U•U mismatches from stem-loop 1 in the 3'-UTR of the SARS-CoV-2 genomic RNA

The single-stranded RNA genome of SARS-CoV-2 is highly structured. Numerous helical stem-loop structures interrupted by mismatch motifs are present in the functionally important 5'- and 3'-UTRs. These mismatches modulate local helical geometries and feature unusual arrays of hydrogen bonding donor and acceptor groups. However, their conformational and dynamical properties cannot be directly inferred from chemical probing and are difficult to predict theoretically. A mismatch motif (SL1-motif) consisting of three consecutive U•U base pairs is located in stem-loop 1 of the 3'-UTR. We combined NMR-spectroscopy and MD-simulations to investigate its structure and dynamics. All three U•U base pairs feature two direct hydrogen bonds and are as stable as Watson-Crick A:U base pairs. Plasmodium falciparum 25S rRNA contains a triple U•U mismatch motif (Pf-motif) differing from SL1-motif only with respect to the orientation of the two closing base pairs. Interestingly, while the geometry of the outer two U•U mismatches was identical in both motifs the preferred orientation of the central U•U mismatch was different. MD simulations and potassium ion titrations revealed that the potassium ion-binding mode to the major groove is connected to the different preferred geometries of the central base pair in the two motifs.

Isoform-specific RNA structure determination using Nano-DMS-MaP

RNA structure determination is essential to understand how RNA carries out its diverse biological functions. In cells, RNA isoforms are readily expressed with partial variations within their sequences due, for example, to alternative splicing, heterogeneity in the transcription start site, RNA processing or differential termination/polyadenylation. Nanopore dimethyl sulfate mutational profiling (Nano-DMS-MaP) is a method for in situ isoform-specific RNA structure determination. Unlike similar methods that rely on short sequencing reads, Nano-DMS-MaP employs nanopore sequencing to resolve the structures of long and highly similar RNA molecules to reveal their previously hidden structural differences. This Protocol describes the development and applications of Nano-DMS-MaP and outlines the main considerations for designing and implementing a successful experiment: from bench to data analysis. In cell probing experiments can be carried out by an experienced molecular biologist in 3-4 d. Data analysis requires good knowledge of command line tools and Python scripts and requires a further 3-5 d.

Nanoscale cellular organization of viral RNA and proteins in SARS-CoV-2 replication organelles

The SARS-CoV-2 viral infection transforms host cells and produces special organelles in many ways, and we focus on the replication organelles, the sites of replication of viral genomic RNA (vgRNA). To date, the precise cellular localization of key RNA molecules and replication intermediates has been elusive in electron microscopy studies. We use super-resolution fluorescence microscopy and specific labeling to reveal the nanoscopic organization of replication organelles that contain numerous vgRNA molecules along with the replication enzymes and clusters of viral double-stranded RNA (dsRNA). We show that the replication organelles are organized differently at early and late stages of infection. Surprisingly, vgRNA accumulates into distinct globular clusters in the cytoplasmic perinuclear region, which grow and accommodate more vgRNA molecules as infection time increases. The localization of endoplasmic reticulum (ER) markers and nsp3 (a component of the double-membrane vesicle, DMV) at the periphery of the vgRNA clusters suggests that replication organelles are encapsulated into DMVs, which have membranes derived from the host ER. These organelles merge into larger vesicle packets as infection advances. Precise co-imaging of the nanoscale cellular organization of vgRNA, dsRNA, and viral proteins in replication organelles of SARS-CoV-2 may inform therapeutic approaches that target viral replication and associated processes.

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

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

Tying the knot: Unraveling the intricacies of the coronavirus frameshift pseudoknot

Understanding and targeting functional RNA structures towards treatment of coronavirus infection can help us to prepare for novel variants of SARS-CoV-2 (the virus causing COVID-19), and any other coronaviruses that could emerge via human-to-human transmission or potential zoonotic (inter-species) events. Leveraging the fact that all coronaviruses use a mechanism known as -1 programmed ribosomal frameshifting (-1 PRF) to replicate, we apply algorithms to predict the most energetically favourable secondary structures (each nucleotide involved in at most one pairing) that may be involved in regulating the -1 PRF event in coronaviruses, especially SARS-CoV-2. We compute previously unknown most stable structure predictions for the frameshift site of coronaviruses via hierarchical folding, a biologically motivated framework where initial non-crossing structure folds first, followed by subsequent, possibly crossing (pseudoknotted), structures. Using mutual information from 181 coronavirus sequences, in conjunction with the algorithm KnotAli, we compute secondary structure predictions for the frameshift site of different coronaviruses. We then utilize the Shapify algorithm to obtain most stable SARS-CoV-2 secondary structure predictions guided by frameshift sequence-specific and genome-wide experimental data. We build on our previous secondary structure investigation of the singular SARS-CoV-2 68 nt frameshift element sequence, by using Shapify to obtain predictions for 132 extended sequences and including covariation information. Previous investigations have not applied hierarchical folding to extended length SARS-CoV-2 frameshift sequences. By doing so, we simulate the effects of ribosome interaction with the frameshift site, providing insight to biological function. We contribute in-depth discussion to contextualize secondary structure dual-graph motifs for SARS-CoV-2, highlighting the energetic stability of the previously identified 3_8 motif alongside the known dominant 3_3 and 3_6 (native-type) -1 PRF structures. Using a combination of thermodynamic methods and sequence covariation, our novel predictions suggest function of the attenuator hairpin via previously unknown pseudoknotted base pairing. While certain initial RNA folding is consistent, other pseudoknotted base pairs form which indicate potential conformational switching between the two structures.

NMR Studies of Genomic RNA in 3' Untranslated Region Unveil Pseudoknot Structure that Initiates Viral RNA Replication in SARS-CoV-2

In the 3' untranslated region of the SARS-CoV-2 virus RNA genome, genomic RNA replication is initiated in the highly conserved region called 3'PK, containing three stem structures (P1pk, P2, and P5). According to one proposed mechanism, P1pk and distal P2 stems switch their structure to a pseudoknot through base-pairing, thereby initiating transcription by recruiting RNA-dependent RNA polymerase complexed with nonstructural proteins (nsp)7 and nsp8. However, experimental evidence of pseudoknot formation or structural switching is unavailable. Using SARS-CoV-2 3'PK fragments, we show that 3'PK adopted stem-loop and pseudoknot forms in a mutually exclusive manner. When P1pk and P2 formed a pseudoknot, the P5 stem, which includes a sequence at the 3' end, exited from the stem-loop structure and opened up. Interaction with the nsp7/nsp8 complex destabilized the stem-loop form but did not alter the pseudoknot form. These results suggest that the interaction between the pseudoknot and nsp7/nsp8 complex transformed the 3' end of viral genomic RNA into single-stranded RNA ready for synthesis, presenting the unique pseudoknot structure as a potential pharmacological target.

Identification of Hammerhead-variant ribozyme sequences in SARS-CoV-2

The SARS-CoV-2 RNA virus and variants, responsible for the COVID-19 pandemic has become endemic, raised a need for further understanding of the viral genome and biology. Despite vast research on SARS-CoV-2, no ribozymes have been found in the virus genome. Here we report the identification of 39 Hammerhead-variant ribozyme sequences (CoV-HHRz) in SARS-CoV-2. These sequences are highly conserved within SARS-CoV-2 variants but show large diversity among other coronaviruses. In vitro CoV-HHRz sequences possess the characteristics of typical ribozymes; cleavage is pH and ion dependent, although their activity is relatively low and Mn2+ is required for cleavage. The cleavage sites of four CoV-HHRz coincide with the breakpoint of expressed subgenomic RNA (sgRNAs) in SARS-CoV-2 transcriptome data suggesting in vivo activity. The CoV-HHRz are involved in processing sgRNAs for ORF7b, ORF 10 and ORF1ab nsp13 which are essential for viral packaging and life cycle.

Conserved structures and dynamics in 5'-proximal regions of Betacoronavirus RNA genomes

Betacoronaviruses are a genus within the Coronaviridae family of RNA viruses. They are capable of infecting vertebrates and causing epidemics as well as global pandemics in humans. Mitigating the threat posed by Betacoronaviruses requires an understanding of their molecular diversity. The development of novel antivirals hinges on understanding the key regulatory elements within the viral RNA genomes, in particular the 5'-proximal region, which is pivotal for viral protein synthesis. Using a combination of cryo-electron microscopy, atomic force microscopy, chemical probing, and computational modeling, we determined the structures of 5'-proximal regions in RNA genomes of Betacoronaviruses from four subgenera: OC43-CoV, SARS-CoV-2, MERS-CoV, and Rousettus bat-CoV. We obtained cryo-electron microscopy maps and determined atomic-resolution models for the stem-loop-5 (SL5) region at the translation start site and found that despite low sequence similarity and variable length of the helical elements it exhibits a remarkable structural conservation. Atomic force microscopy imaging revealed a common domain organization and a dynamic arrangement of structural elements connected with flexible linkers across all four Betacoronavirus subgenera. Together, these results reveal common features of a critical regulatory region shared between different Betacoronavirus RNA genomes, which may allow targeting of these RNAs by broad-spectrum antiviral therapeutics.

SARS-CoV-2 biology and host interactions

The zoonotic emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and the ensuing coronavirus disease 2019 (COVID-19) pandemic have profoundly affected our society. The rapid spread and continuous evolution of new SARS-CoV-2 variants continue to threaten global public health. Recent scientific advances have dissected many of the molecular and cellular mechanisms involved in coronavirus infections, and large-scale screens have uncovered novel host-cell factors that are vitally important for the virus life cycle. In this Review, we provide an updated summary of the SARS-CoV-2 life cycle, gene function and virus-host interactions, including recent landmark findings on general aspects of coronavirus biology and newly discovered host factors necessary for virus replication.

Nanoscale cellular organization of viral RNA and proteins in SARS-CoV-2 replication organelles

The SARS-CoV-2 viral infection transforms host cells and produces special organelles in many ways, and we focus on the replication organelle where the replication of viral genomic RNA (vgRNA) occurs. To date, the precise cellular localization of key RNA molecules and replication intermediates has been elusive in electron microscopy studies. We use super-resolution fluorescence microscopy and specific labeling to reveal the nanoscopic organization of replication organelles that contain vgRNA clusters along with viral double-stranded RNA (dsRNA) clusters and the replication enzyme, encapsulated by membranes derived from the host endoplasmic reticulum (ER). We show that the replication organelles are organized differently at early and late stages of infection. Surprisingly, vgRNA accumulates into distinct globular clusters in the cytoplasmic perinuclear region, which grow and accommodate more vgRNA molecules as infection time increases. The localization of ER labels and nsp3 (a component of the double-membrane vesicle, DMV) at the periphery of the vgRNA clusters suggests that replication organelles are enclosed by DMVs at early infection stages which then merge into vesicle packets as infection progresses. Precise co-imaging of the nanoscale cellular organization of vgRNA, dsRNA, and viral proteins in replication organelles of SARS-CoV-2 may inform therapeutic approaches that target viral replication and associated processes.

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

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

RNAvigate: efficient exploration of RNA chemical probing datasets

Chemical probing technologies enable high-throughput examination of diverse structural features of RNA, including local nucleotide flexibility, RNA secondary structure, protein and ligand binding, through-space interaction networks, and multistate structural ensembles. Deep understanding of RNA structure-function relationships typically requires evaluating a system under structure- and function-altering conditions, linking these data with additional information, and visualizing multilayered relationships. Current platforms lack the broad accessibility, flexibility and efficiency needed to iterate on integrative analyses of these diverse, complex data. Here, we share the RNA visualization and graphical analysis toolset RNAvigate, a straightforward and flexible Python library that automatically parses 21 standard file formats (primary sequence annotations, per- and internucleotide data, and secondary and tertiary structures) and outputs 18 plot types. RNAvigate enables efficient exploration of nuanced relationships between multiple layers of RNA structure information and across multiple experimental conditions. Compatibility with Jupyter notebooks enables nonburdensome, reproducible, transparent and organized sharing of multistep analyses and data visualization strategies. RNAvigate simplifies and accelerates discovery and characterization of RNA-centric functions in biology.

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

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

Graph Convolutional Network for predicting secondary structure of RNA

The prediction of RNA secondary structures is essential for understanding its underlying principles and applications in diverse fields, including molecular diagnostics and RNA-based therapeutic strategies. However, the complexity of the search space presents a challenge. This work proposes a Graph Convolutional Network (GCNfold) for predicting the RNA secondary structure. GCNfold considers an RNA sequence as graph-structured data and predicts posterior base-pairing probabilities given the prior base-pairing probabilities, calculated using McCaskill's partition function. The performance of GCNfold surpasses that of the state-of-the-art folding algorithms, as we have incorporated minimum free energy information into the richly parameterized network, enhancing its robustness in predicting non-homologous RNA secondary structures. A Symmetric Argmax Post-processing algorithm ensures that GCNfold formulates valid structures. To validate our algorithm, we applied it to the SARS-CoV-2 E gene and determined the secondary structure of the E-gene across the Betacoronavirus subgenera.

Causes, functions, and therapeutic possibilities of RNA secondary structure ensembles and alternative states

RNA secondary structure plays essential roles in encoding RNA regulatory fate and function. Most RNAs populate ensembles of alternatively paired states and are continually unfolded and refolded by cellular processes. Measuring these structural ensembles and their contributions to cellular function has traditionally posed major challenges, but new methods and conceptual frameworks are beginning to fill this void. In this review, we provide a mechanism- and function-centric compendium of the roles of RNA secondary structural ensembles and minority states in regulating the RNA life cycle, from transcription to degradation. We further explore how dysregulation of RNA structural ensembles contributes to human disease and discuss the potential of drugging alternative RNA states to therapeutically modulate RNA activity. The emerging paradigm of RNA structural ensembles as central to RNA function provides a foundation for a deeper understanding of RNA biology and new therapeutic possibilities.

Specific binding of G-quadruplex in SARS-CoV-2 RNA by RHAU peptide

G-quadruplexes (G4s) are reported to present on the SARS-CoV-2 RNA genome and control various viral activities. Specific ligands targeting those viral nucleic acid structures could be investigated as promising detection methods or antiviral reagents to suppress this menacing virus. Herein, we demonstrate the binding between a G4 structure in the RNA of SARS-CoV-2 and a fluorescent probe created by fusing a parallel-G4 specific RHAU53 and a cyan fluorescent protein. The specific binding of G4 in SARS-CoV-2 by RHAU peptide was easily detected under the fluorescence spectrometer. The drawbacks of this approach and potential solutions are also discussed.

Decoding the genome of SARS-CoV-2: a pathway to drug development through translation inhibition

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes the coronavirus disease 2019 (COVID-19) pandemic and is continuously spreading globally. The continuous emergence of new SARS-CoV-2 variants keeps posing threats, highlighting the need for fast-acting, mutation-resistant broad-spectrum therapeutics. Protein translation is vital for SARS-CoV-2 replication, producing early non-structural proteins for RNA replication and transcription, and late structural proteins for virion assembly. Targeted blocking of viral protein translation is thus a potential approach to developing effective anti-SARS-CoV-2 drugs. SARS-CoV-2, as an obligate parasite, utilizes the host's translation machinery. Translation-blocking strategies that target the SARS-CoV-2 mRNA, especially those that target its conserved elements are generally preferred. In this review, we discuss the current understanding of SARS-CoV-2 translation, highlighting the important conserved motifs and structures involved in its regulation. We also discuss the current strategies for blocking SARS-CoV-2 translation through viral RNA degradation or RNA element dysfunction.