Viral RNA pseudoknots: versatile motifs in gene expression and replication

RNase P cleavage of pseudoknot substrates reveals differences in active site architecture that depend on residue N-1 in the 5' leader

We show that a small biotin-binding RNA aptamer that folds into a pseudoknot structure acts as a substrate for bacterial RNase P RNA (RPR) with and without the RNase P C5 protein. Cleavage in the single-stranded region in loop 1 was shown to depend on the presence of a RCCA-motif at the 3' end of the substrate. The nucleobase and the 2'hydroxyl at the position immediately 5' of the cleavage site contribute to both cleavage efficiency and site selection, where C at this position induces significant cleavage at an alternative site, one base upstream of the main cleavage site. The frequencies of cleavage at these two sites and Mg2+ binding change upon altering the structural topology in the vicinity of the cleavage site as well as by replacing Mg2+ with other divalent metal ions. Modelling studies of RPR in complex with the pseudoknot substrates suggest alternative structural topologies for cleavage at the main and the alternative site and a shift in positioning of Mg2+ that activates the H2O nucleophile. Together, our data are consistent with a model where the organization of the active site structure and positioning of Mg2+ is influenced by the identities of residues at and in the vicinity of the site of cleavage.

Optimizing On-the-Fly Probability Enhanced Sampling for Complex RNA Systems: Sampling Free Energy Surfaces of an H-Type Pseudoknot

All-atom molecular dynamics (MD) simulations offer crucial insights into biomolecular dynamics, but inherent time scale constraints often limit their effectiveness. Advanced sampling techniques help overcome these limitations, enabling predictions of deeply rugged folding free energy surfaces (FES) of RNA at atomistic resolution. The Multithermal-Multiumbrella On-the-Fly Probability Enhanced Sampling (MM-OPES) method, which combines temperature and collective variables (CVs) to accelerate sampling, has shown promise and cost-effectiveness. However, the applications have so far been limited to simpler RNA systems, such as stem-loops. In this study, we optimized the MM-OPES method to explore the FES of an H-type RNA pseudoknot, a more complex fundamental RNA folding unit. Through systematic exploration of CV combinations and temperature ranges, we identified an optimal strategy for both sampling and analysis. Our findings demonstrate that treating the native-like contacts in two stems as independent CVs and using a temperature range of 300-480 K provides the most effective sampling, while projections onto native Watson-Crick-type hydrogen bond CVs yield the best resolution FES prediction. Additionally, our sampling scheme also revealed various folding/unfolding pathways. This study provides practical insights and detailed decision-making strategies for adopting the MM-OPES method, facilitating its application to complex RNA structures at atomistic resolution.

Live-cell RNA imaging with the inactivated endonuclease Csy4 enables new insights into plant virus transport through plasmodesmata

Plant-infecting viruses spread through their hosts by transporting their infectious genomes through intercellular nano-channels called plasmodesmata. This process is mediated by virus-encoded movement proteins. Whilst the sub-cellular localisations of movement proteins have been intensively studied, live-cell RNA imaging systems have so far not been able to detect viral genomes inside the plasmodesmata. Here, we describe a highly sensitive RNA live-cell reporter based on an enzymatically inactive form of the small bacterial endonuclease Csy4, which binds to its cognate stem-loop with picomolar affinity. This system allows imaging of plant viral RNA genomes inside plasmodesmata and shows that potato virus X RNA remains accessible within the channels and is therefore not fully encapsidated during movement. We also combine Csy4-based RNA-imaging with interspecies movement complementation to show that an unrelated movement protein from tobacco mosaic virus can recruit potato virus X replication complexes adjacent to plasmodesmata. Therefore, recruitment of potato virus X replicase is mediated non-specifically, likely by indirect coupling of movement proteins and viral replicase via the viral RNA or co-compartmentalisation, potentially contributing to transport specificity. Lastly, we show that a 'self-tracking' virus can express the Csy4-based reporter during the progress of infection. However, expression of the RNA-binding protein in cis interferes with viral movement by an unidentified mechanism when cognate stem-loops are present in the viral RNA.

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

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

Molecular Dynamics Characterization of the Free and Encapsidated RNA2 of CCMV with the oxRNA Model

The cowpea chlorotic mottle virus (CCMV) has emerged as a model system to assess the balance between electrostatic and topological features of single-stranded RNA viruses, specifically in the context of the viral self-assembly. Yet, despite its biophysical significance, little structural data on the RNA content of the CCMV virion is available. Here, the conformational dynamics of the RNA2 fragment of CCMV was assessed via coarse-grained molecular dynamics simulations, employing the oxRNA2 force field. The behavior of RNA2 was characterized both as a freely-folding molecule and within a mean-field depiction of the capsid. For the former, the role of the salt concentration, the temperature and of ad hoc constraints on the RNA termini was verified on the equilibrium properties of RNA2. For the latter, a multi-scale approach was employed to derive a potential profile of the viral cavity from atomistic structures of the CCMV capsid in solution. The conformational ensembles of the encapsidated RNA2 were significantly altered with respect to the freely-folding counterparts, as shown by the emergence of long-range motifs and pseudoknots. Finally, the role of the N-terminal tails of the CCMV subunits is highlighted as a critical feature in the construction of a proper electrostatic model of the CCMV capsid.

Magnesium Regulates RNA Ring Dynamics and Folding in Subgenomic Flaviviral RNA

Mosquito-borne flaviviruses including dengue, Zika, yellow fever, and regional encephalitis produce a large amount of short subgenomic flaviviral RNAs during infection. A segment of these RNAs named as xrRNA1 features a multi-pseudoknot (PK)-associated structure, which resists the host cell enzyme (XRN1) from degrading the viral RNA. We investigate how this long-range RNA PK folds in the presence of counterions, specifically in a mix of monovalent (K+) and divalent (Mg2+) salts at physiological concentrations. In this study, we use extensive explicit solvent molecular dynamics (MD) simulations to characterize the RNA ion environment of the folded RNA conformation, as determined by the crystal structure. This allowed us to identify the precise locations of various coordinated RNA-Mg2+ interactions, including inner-sphere/chelated and outer-sphere coordinated Mg2+. Given that RNA folding involves large-scale conformational changes, making it challenging to explore through classical MD simulations, we investigate the folding mechanism of xrRNA1 using an all-atom structure-based RNA model with a hybrid implicit-explicit treatment of the ion environment via the dynamic counterion condensation model, both with and without physiological Mg2+ concentration. The study reveals potential folding pathways for this xrRNA1, which is consistent with the results obtained from optical tweezer experiments. The equilibrium and free energy simulations both capture a dynamic equilibrium between the ring-open and ring-close states of the RNA, driven by a long-range PK interaction. Free energy calculations reveal that with the addition of Mg2+ ions, the equilibrium shifts more toward the ring-close state. A detailed analysis of the free energy pathways and ion-mediated contact probability map highlights the critical role of Mg2+ in bridging G50 and A33. This Mg2+-mediated connection helps form the long-range PK which in turn controls the transition between the ring-open and ring-close states. The study underscores the critical role of Mg2+ in the RNA folding transition, highlighting specific locations of Mg2+ contributing to the stabilization of long-range PK connections likely to enhance the robustness of Xrn1 resistance of flaviviral xrRNAs.

RNA-Binding Small Molecules in Drug Discovery and Delivery: An Overview from Fundamentals

RNA molecules, similar to proteins, fold into complex structures to confer diverse functions in cells. The intertwining of functions with RNA structures offers a new therapeutic opportunity for small molecules to bind and manipulate disease-relevant RNA pathways, thus creating a therapeutic realm of RNA-binding small molecules. The ongoing interest in RNA targeting and subsequent screening campaigns have led to the identification of numerous compounds that can regulate RNAs from splicing, degradation to malfunctions, with therapeutic benefits for a variety of diseases. Moreover, along with the rise of RNA-based therapeutics, RNA-binding small molecules have expanded their application to the modification, regulation, and delivery of RNA drugs, leading to the burgeoning interest in this field. This Perspective overviews the emerging roles of RNA-binding small molecules in drug discovery and delivery, covering aspects from their action fundamentals to therapeutic applications, which may inspire researchers to advance the field.

Mapping the structural landscape of the yeast Ty3 retrotransposon RNA genome

Long terminal repeat (LTR)-retrotransposons are significant contributors to the evolution and diversity of eukaryotic genomes. Their RNA genomes (gRNA) serve as a template for protein synthesis and reverse transcription to a DNA copy, which can integrate into the host genome. Here, we used the SHAPE-MaP strategy to explore Ty3 retrotransposon gRNA structure in yeast and under cell-free conditions. Our study reveals the structural dynamics of Ty3 gRNA and the well-folded core, formed independently of the cellular environment. Based on the detailed map of Ty3 gRNA structure, we characterized the structural context of cis-acting sequences involved in reverse transcription and frameshifting. We also identified a novel functional sequence as a potential initiator for Ty3 gRNA dimerization. Our data indicate that the dimer is maintained by direct interaction between short palindromic sequences at the 5' ends of the two Ty3 gRNAs, resembling the model characteristic for other retroelements like HIV-1 and Ty1. This work points out a range of cell-dependent and -independent Ty3 gRNA structural changes that provide a solid background for studies on RNA structure-function relationships important for retroelement biology.

Help or Hinder: Protein Host Factors That Impact HIV-1 Replication

Interactions between human immunodeficiency virus type 1 (HIV-1) and the host factors or restriction factors of its target cells determine the cell's susceptibility to, and outcome of, infection. Factors intrinsic to the cell are involved at every step of the HIV-1 replication cycle, contributing to productive infection and replication, or severely attenuating the chances of success. Furthermore, factors unique to certain cell types contribute to the differences in infection between these cell types. Understanding the involvement of these factors in HIV-1 infection is a key requirement for the development of anti-HIV-1 therapies. As the list of factors grows, and the dynamic interactions between these factors and the virus are elucidated, comprehensive and up-to-date summaries that recount the knowledge gathered after decades of research are beneficial to the field, displaying what is known so that researchers can build off the groundwork of others to investigate what is unknown. Herein, we aim to provide a review focusing on protein host factors, both well-known and relatively new, that impact HIV-1 replication in a positive or negative manner at each stage of the replication cycle, highlighting factors unique to the various HIV-1 target cell types where appropriate.

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

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

Recruitment of the 40S ribosomal subunit by the West Nile virus 3' UTR promotes the cross-talk between the viral genomic ends for translation regulation

Flaviviral RNA genomes are composed of discrete RNA structural units arranged in an ordered fashion and grouped into complex folded domains that regulate essential viral functions, e.g. replication and translation. This is achieved by adjusting the overall structure of the RNA genome via the establishment of inter- and intramolecular interactions. Translation regulation is likely the main process controlling flaviviral gene expression. Although the genomic 3' UTR is a key player in this regulation, little is known about the molecular mechanisms underlying this role. The present work provides evidence for the specific recruitment of the 40S ribosomal subunit by the 3' UTR of the West Nile virus RNA genome, showing that the joint action of both genomic ends contributes the positioning of the 40S subunit at the 5' end. The combination of structural mapping techniques revealed specific conformational requirements at the 3' UTR for 40S binding, involving the highly conserved SL-III, 5'DB, 3'DB and 3'SL elements, all involved in the translation regulation. These results point to the 40S subunit as a bridge to ensure cross-talk between both genomic ends during viral translation and support a link between 40S recruitment by the 3' UTR and translation control.

Unveiling hidden structural patterns in the SARS-CoV-2 genome: Computational insights and comparative analysis

SARS-CoV-2, the causative agent of COVID-19, is known to exhibit secondary structures in its 5' and 3' untranslated regions, along with the frameshifting stimulatory element situated between ORF1a and 1b. To identify additional regions containing conserved structures, we utilized a multiple sequence alignment with related coronaviruses as a starting point. We applied a computational pipeline developed for identifying non-coding RNA elements. Our pipeline employed three different RNA structural prediction approaches. We identified forty genomic regions likely to harbor structures, with ten of them showing three-way consensus substructure predictions among our predictive utilities. We conducted intracomparisons of the predictive utilities within the pipeline and intercomparisons with four previously published SARS-CoV-2 structural datasets. While there was limited agreement on the precise structure, different approaches seemed to converge on regions likely to contain structures in the viral genome. By comparing and combining various computational approaches, we can predict regions most likely to form structures, as well as a probable structure or ensemble of structures. These predictions can be used to guide surveillance, prophylactic measures, or therapeutic efforts. Data and scripts employed in this study may be found at https://doi.org/10.5281/zenodo.8298680.

The non-template functions of helper virus RNAs create optimal replication conditions to enhance the proliferation of satellite RNAs

As a type of parasitic agent, satellite RNAs (satRNAs) rely on cognate helper viruses to achieve their replication and transmission. During the infection of satRNAs, helper virus RNAs serve as templates for synthesizing viral proteins, including the replication proteins essential for satRNA replication. However, the role of non-template functions of helper virus RNAs in satRNA replication remains unexploited. Here we employed the well-studied model that is composed of cucumber mosaic virus (CMV) and its associated satRNA. In the experiments employing the CMV trans-replication system, we observed an unexpected phenomenon the replication proteins of the mild strain LS-CMV exhibited defective in supporting satRNA replication, unlike those of the severe strain Fny-CMV. Independent of translation products, all CMV genomic RNAs could enhance satRNA replication, when combined with the replication proteins of CMV. This enhancement is contingent upon the recruitment and complete replication of helper virus RNAs. Using the method developed for analyzing the satRNA recruitment, we observed a markedly distinct ability of the replication proteins from both CMV strains to recruit the positive-sense satRNA-harboring RNA3 mutant for replication. This is in agreement with the differential ability of both 1a proteins in binding satRNAs in plants. The discrepancies provide a convincing explanation for the variation of the replication proteins of both CMV strains in replicating satRNAs. Taken together, our work provides compelling evidence that the non-template functions of helper virus RNAs create an optimal replication environment to enhance satRNA proliferation.

RNA-Small-Molecule Interaction: Challenging the "Undruggable" Tag

RNA targeting, specifically with small molecules, is a relatively new and rapidly emerging avenue with the promise to expand the target space in the drug discovery field. From being "disregarded" as an "undruggable" messenger molecule to FDA approval of an RNA-targeting small-molecule drug Risdiplam, a radical change in perspective toward RNA has been observed in the past decade. RNAs serve important regulatory functions beyond canonical protein synthesis, and their dysregulation has been reported in many diseases. A deeper understanding of RNA biology reveals that RNA molecules can adopt a variety of structures, carrying defined binding pockets that can accommodate small-molecule drugs. Due to its functional diversity and structural complexity, RNA can be perceived as a prospective target for therapeutic intervention. This perspective highlights the proof of concept of RNA-small-molecule interactions, exemplified by targeting of various transcripts with functional modulators. The advent of RNA-oriented knowledge would help expedite drug discovery.

Accurate prediction of RNA secondary structure including pseudoknots through solving minimum-cost flow with learned potentials

Pseudoknots are key structure motifs of RNA and pseudoknotted RNAs play important roles in a variety of biological processes. Here, we present KnotFold, an accurate approach to the prediction of RNA secondary structure including pseudoknots. The key elements of KnotFold include a learned potential function and a minimum-cost flow algorithm to find the secondary structure with the lowest potential. KnotFold learns the potential from the RNAs with known structures using an attention-based neural network, thus avoiding the inaccuracy of hand-crafted energy functions. The specially designed minimum-cost flow algorithm used by KnotFold considers all possible combinations of base pairs and selects from them the optimal combination. The algorithm breaks the restriction of nested base pairs required by the widely used dynamic programming algorithms, thus enabling the identification of pseudoknots. Using 1,009 pseudoknotted RNAs as representatives, we demonstrate the successful application of KnotFold in predicting RNA secondary structures including pseudoknots with accuracy higher than the state-of-the-art approaches. We anticipate that KnotFold, with its superior accuracy, will greatly facilitate the understanding of RNA structures and functionalities.

A Hitchhiker's guide to RNA-RNA structure and interaction prediction tools

RNA biology has risen to prominence after a remarkable discovery of diverse functions of noncoding RNA (ncRNA). Most untranslated transcripts often exert their regulatory functions into RNA-RNA complexes via base pairing with complementary sequences in other RNAs. An interplay between RNAs is essential, as it possesses various functional roles in human cells, including genetic translation, RNA splicing, editing, ribosomal RNA maturation, RNA degradation and the regulation of metabolic pathways/riboswitches. Moreover, the pervasive transcription of the human genome allows for the discovery of novel genomic functions via RNA interactome investigation. The advancement of experimental procedures has resulted in an explosion of documented data, necessitating the development of efficient and precise computational tools and algorithms. This review provides an extensive update on RNA-RNA interaction (RRI) analysis via thermodynamic- and comparative-based RNA secondary structure prediction (RSP) and RNA-RNA interaction prediction (RIP) tools and their general functions. We also highlighted the current knowledge of RRIs and the limitations of RNA interactome mapping via experimental data. Then, the gap between RSP and RIP, the importance of RNA homologues, the relationship between pseudoknots, and RNA folding thermodynamics are discussed. It is hoped that these emerging prediction tools will deepen the understanding of RNA-associated interactions in human diseases and hasten treatment processes.

Identification of over ten thousand candidate structured RNAs in viruses and phages

Structured RNAs play crucial roles in viruses, exerting influence over both viral and host gene expression. However, the extensive diversity of structured RNAs and their ability to act in cis or trans positions pose challenges for predicting and assigning their functions. While comparative genomics approaches have successfully predicted candidate structured RNAs in microbes on a large scale, similar efforts for viruses have been lacking. In this study, we screened over 5 million DNA and RNA viral sequences, resulting in the prediction of 10,006 novel candidate structured RNAs. These predictions are widely distributed across taxonomy and ecosystem. We found transcriptional evidence for 206 of these candidate structured RNAs in the human fecal microbiome. These candidate RNAs exhibited evidence of nucleotide covariation, indicative of selective pressure maintaining the predicted secondary structures. Our analysis revealed a diverse repertoire of candidate structured RNAs, encompassing a substantial number of putative tRNAs or tRNA-like structures, Rho-independent transcription terminators, and potentially cis-regulatory structures consistently positioned upstream of genes. In summary, our findings shed light on the extensive diversity of structured RNAs in viruses, offering a valuable resource for further investigations into their functional roles and implications in viral gene expression and pave the way for a deeper understanding of the intricate interplay between viruses and their hosts at the molecular level.

The plant noncoding transcriptome: a versatile environmental sensor

Plant noncoding RNA transcripts have gained increasing attention in recent years due to growing evidence that they can regulate developmental plasticity. In this review article, we comprehensively analyze the relationship between noncoding RNA transcripts in plants and their response to environmental cues. We first provide an overview of the various noncoding transcript types, including long and small RNAs, and how the environment modulates their performance. We then highlight the importance of noncoding RNA secondary structure for their molecular and biological functions. Finally, we discuss recent studies that have unveiled the functional significance of specific long noncoding transcripts and their molecular partners within ribonucleoprotein complexes during development and in response to biotic and abiotic stress. Overall, this review sheds light on the fascinating and complex relationship between dynamic noncoding transcription and plant environmental responses, and highlights the need for further research to uncover the underlying molecular mechanisms and exploit the potential of noncoding transcripts for crop resilience in the context of global warming.

Improved prime editing allows for routine predictable gene editing in Physcomitrium patens

Efficient and precise gene editing is the gold standard of any reverse genetic study. The recently developed prime editing approach, a modified CRISPR/Cas9 [clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein] editing method, has reached the precision goal but its editing rate can be improved. We present an improved methodology that allows for routine prime editing in the model plant Physcomitrium patens, whilst exploring potential new prime editing improvements. Using a standardized protoplast transfection procedure, multiple prime editing guide RNA (pegRNA) structural and prime editor variants were evaluated targeting the APT reporter gene through direct plant selection. Together, enhancements of expression of the prime editor, modifications of the 3' extension of the pegRNA, and the addition of synonymous mutation in the reverse transcriptase template sequence of the pegRNA dramatically improve the editing rate without affecting the quality of the edits. Furthermore, we show that prime editing is amenable to edit a gene of interest through indirect selection, as demonstrated by the generation of a Ppdek10 mutant. Additionally, we determine that a plant retrotransposon reverse transcriptase enables prime editing. Finally, we show for the first time the possibility of performing prime editing with two independently coded peptides.

Structural and Functional Insights into Viral Programmed Ribosomal Frameshifting

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

Small Molecules Targeting Viral RNA

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

Machine learning modeling of RNA structures: methods, challenges and future perspectives

The three-dimensional structure of RNA molecules plays a critical role in a wide range of cellular processes encompassing functions from riboswitches to epigenetic regulation. These RNA structures are incredibly dynamic and can indeed be described aptly as an ensemble of structures that shifts in distribution depending on different cellular conditions. Thus, the computational prediction of RNA structure poses a unique challenge, even as computational protein folding has seen great advances. In this review, we focus on a variety of machine learning-based methods that have been developed to predict RNA molecules' secondary structure, as well as more complex tertiary structures. We survey commonly used modeling strategies, and how many are inspired by or incorporate thermodynamic principles. We discuss the shortcomings that various design decisions entail and propose future directions that could build off these methods to yield more robust, accurate RNA structure predictions.

Toxin-antitoxin systems in bacterial pathogenesis

Toxin-Antitoxin (TA) systems are abundant in prokaryotes and play an important role in various biological processes such as plasmid maintenance, phage inhibition, stress response, biofilm formation, and dormant persister cell generation. TA loci are abundant in pathogenic intracellular micro-organisms and help in their adaptation to the harsh host environment such as nutrient deprivation, oxidation, immune response, and antimicrobials. Several studies have reported the involvement of TA loci in establishing successful infection, intracellular survival, better colonization, adaptation to host stresses, and chronic infection. Overall, the TA loci play a crucial role in bacterial virulence and pathogenesis. Nonetheless, there are some controversies about the role of TA system in stress response, biofilm and persister formation. In this review, we describe the role of the TA systems in bacterial virulence. We discuss the important features of each type of TA system and the recent discoveries identifying key contributions of TA loci in bacterial pathogenesis.

Looking at Biomolecular Interactions through the Lens of Correlated Fluorescence Microscopy and Optical Tweezers

Understanding complex biological events at the molecular level paves the path to determine mechanistic processes across the timescale necessary for breakthrough discoveries. While various conventional biophysical methods provide some information for understanding biological systems, they often lack a complete picture of the molecular-level details of such dynamic processes. Studies at the single-molecule level have emerged to provide crucial missing links to understanding complex and dynamic pathways in biological systems, which are often superseded by bulk biophysical and biochemical studies. Latest developments in techniques combining single-molecule manipulation tools such as optical tweezers and visualization tools such as fluorescence or label-free microscopy have enabled the investigation of complex and dynamic biomolecular interactions at the single-molecule level. In this review, we present recent advances using correlated single-molecule manipulation and visualization-based approaches to obtain a more advanced understanding of the pathways for fundamental biological processes, and how this combination technique is facilitating research in the dynamic single-molecule (DSM), cell biology, and nanomaterials fields.

Genomic Analysis of Non-B Nucleic Acids Structures in SARS-CoV-2: Potential Key Roles for These Structures in Mutability, Translation, and Replication?

Non-B nucleic acids structures have arisen as key contributors to genetic variation in SARS-CoV-2. Herein, we investigated the presence of defining spike protein mutations falling within inverted repeats (IRs) for 18 SARS-CoV-2 variants, discussed the potential roles of G-quadruplexes (G4s) in SARS-CoV-2 biology, and identified potential pseudoknots within the SARS-CoV-2 genome. Surprisingly, there was a large variation in the number of defining spike protein mutations arising within IRs between variants and these were more likely to occur in the stem region of the predicted hairpin stem-loop secondary structure. Notably, mutations implicated in ACE2 binding and propagation (e.g., ΔH69/V70, N501Y, and D614G) were likely to occur within IRs, whilst mutations involved in antibody neutralization and reduced vaccine efficacy (e.g., T19R, ΔE156, ΔF157, R158G, and G446S) were rarely found within IRs. We also predicted that RNA pseudoknots could predominantly be found within, or next to, 29 mutations found in the SARS-CoV-2 spike protein. Finally, the Omicron variants BA.2, BA.4, BA.5, BA.2.12.1, and BA.2.75 appear to have lost two of the predicted G4-forming sequences found in other variants. These were found in nsp2 and the sequence complementary to the conserved stem-loop II-like motif (S2M) in the 3' untranslated region (UTR). Taken together, non-B nucleic acids structures likely play an integral role in SARS-CoV-2 evolution and genetic diversity.

Computer-Assisted Design and Characterization of RNA Nanostructures

Molecular dynamics (MD) simulations can aid in the design and characterization of RNA nanomaterials, providing details about structural and dynamical properties as a function of sequence and environment. Here, we describe how to perform explicit and implicit solvent all-atom MD simulations for RNA nanoring systems.

Examples of Structural Motifs in Viral Genomes and Approaches for RNA Structure Characterization

The relationship between conserved structural motifs and their biological function in the virus replication cycle is the interest of many researchers around the world. RNA structure is closely related to RNA function. Therefore, technological progress in high-throughput approaches for RNA structure analysis and the development of new ones are very important. In this mini review, we discuss a few perspectives on the structural elements of viral genomes and some methods used for RNA structure prediction and characterization. Based on the recent literature, we describe several examples of studies concerning the viral genomes, especially severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and influenza A virus (IAV). Herein, we emphasize that a better understanding of viral genome architecture allows for the discovery of the structure-function relationship, and as a result, the discovery of new potential antiviral therapeutics.

Dynamic effects of the spine of hydrated magnesium on viral RNA pseudoknot structure

In the cellular environment, a viral RNA Pseudoknot (PK) structure is responsive to its prevailing ion atmosphere created by a mixture of monovalent and divalent cations. We investigate the influence of such a mixed-salt environment on RNA-PK structure at an atomic resolution through three sets of 1.5 μs explicit solvent molecular dynamics (MD) simulations and also by building a dynamic counterion-condensation (DCC) model at varying divalent Mg²⁺ concentrations. The DCC model includes explicit interaction of the ligand and adjacent chelated Mg²⁺ by extending the recently developed generalized Manning condensation model. Its implementation within an all-atom structure-based molecular dynamics framework bolsters its opportunity to explore large-length scale and long-timescale phenomena associated with complex RNA systems immersed in its dynamic ion environment. In the present case of RNA-PK, both explicit MD and DCC simulations reveal a spine of hydrated-Mg²⁺ to induce stem-I and stem-II closure where the minor groove between these stems is akin to breathing. Mg²⁺ mediated minor-groove narrowing is coupled with local base-flipping dynamics of a base triple and quadruple, changing the stem structure of such RNA. Contrary to the conversational view of the indispensable need for Mg²⁺ for the tertiary structure of RNA, the study warns about the plausible detrimental effect of specific Mg²⁺-phosphate interactions on the RNA-PK structure beyond a certain concentration of Mg²⁺.

The 5′-end motif of Senecavirus A cDNA clone is genetically modified in 36 different ways for uncovering profiles of virus recovery

Senecavirus A (SVA) is an emerging picornavirus. Its genome is one positive-sense, single-stranded RNA. The viral protein (VPg) is covalently linked to the extreme 5′ end of the SVA genome. A complex hairpin-pseudoknot-hairpin (HPH) RNA structure was computationally predicted to form at the 5′ end of the SVA genome. A total of three extra “U” residues (UUU) served as a linker between the HPH structure and the VPg, causing putative UUU–HPH formation at the extreme 5′ end of the SVA genome. It is unclear how the UUU–HPH structure functions. One SVA cDNA clone (N0) was constructed previously in our laboratory. Here, the N0 was genetically tailored for reconstructing a set of 36 modified cDNA clones (N1 to N36) in an attempt to rescue replication-competent SVAs using reverse genetics. The results showed that a total of nine viruses were successfully recovered. Out of them, five were independently rescued from the N1 to N5, reconstructed by deleting the first five nucleotides (TTTGA) one by one from the extreme 5′ end of N0. Interestingly, these five viral progenies reverted to the wild-type or/and wild-type-like genotype, suggesting that SVA with an ability to repair nucleotide defects in its extreme 5′ end. The other four were independently rescued from the N26 to N29, containing different loop-modifying motifs in the first hairpin of the HPH structure. These four loop-modifying motifs were genetically stable after serial passages, implying the wild-type loop motif was not a high-fidelity element in the first hairpin during SVA replication. The other genetically modified sequences were demonstrated to be lethal elements in the HPH structure for SVA recovery, suggesting that the putative HPH formation was a crucial cis-acting replication element for SVA propagation.

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.

Highly sensitive EGFRvIII detection in circulating extracellular vesicle RNA of glioma patients

PURPOSE

Liquid biopsy offers an attractive platform for non-invasive tumor diagnosis, prognostication and prediction of glioblastoma clinical outcomes. Prior studies report that 30-50% of GBM lesions characterized by EGFR amplification also harbor the EGFRvIII mutation.

EXPERIMENTAL DESIGN

A novel digital droplet PCR (ddPCR) assay for high GC content amplicons was developed and optimized for sensitive detection of EGFRvIII in tumor tissue and circulating extracellular vesicle RNA (EV RNA) isolated from the plasma of glioma patients.

RESULTS

Our optimized qPCR assay detected EGFRvIII mRNA in 81% (95% CI, 68% - 94%) of EGFR amplified glioma tumor tissue, indicating a higher than previously reported prevalence of EGFRvIII in glioma. Using the optimized ddPCR assay in discovery and blinded validation cohorts, we detected EGFRvIII mutation in 73% (95% CI, 64% - 82%) of patients with a specificity of 98% (95% CI, 87% - 100%), compared with qPCR tumor tissue analysis. Additionally, upon longitudinal monitoring in 4 patients, we report detection of EGFRvIII in the plasma of patients with different clinical outcomes, rising with tumor progression, and decreasing in response to treatment.

CONCLUSION

This study demonstrates the feasibility of detecting EGFRvIII mutation in plasma using a highly sensitive and specific ddPCR assay. We also show a higher than previously reported EGFRvIII prevalence in glioma tumor tissue. Several features of the assay are favorable for clinical implementation for detection and monitoring of EGFRvIII positive tumors.

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.

Applicability of the nearest-neighbour model for pseudoknot RNAs

The validity of the nearest-neighbour (NN) model was verified in an RNA pseudoknot (PK) structure. The thermodynamic parameters of the second hairpin stem (S2) region, which separates the PK from a hairpin structure, were monitored using CD and UV melting. Different PKs with identical NN base pairs in the S2 region exhibited similar thermodynamic parameters, highlighting the validity of the NN model in this RNA tertiary structure motif.

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.

Development of RNA G-quadruplex (rG4)-targeting L-RNA aptamers by rG4-SELEX

RNA G-quadruplex (rG4)-SELEX is a method that generates L-RNA aptamers to target an rG4 structure of interest, which can be applied to inhibit G-quadruplex-mediated interactions that have important roles in gene regulation and function. Here we present a Protocol Extension substantially modifying an existing SELEX protocol to describe in detail the procedures involved in performing rG4-SELEX to identify rG4-specific binders that can effectively suppress rG4-peptide and rG4-protein associations. This Protocol Extension improves the speed of aptamer discovery and identification, offering a suite of techniques to characterize the aptamer secondary structure and monitor binding affinity and specificity, and demonstrating the utility of the L-RNA aptamer. The previous protocol mainly describes the identification of RNA aptamers against proteins of interest, whereas in this Protocol Extension we present the development of an unnatural RNA aptamer against an RNA structure of interest, with the potential to be applicable to other nucleic acid motifs or biomolecules. rG4-SELEX starts with a random D-RNA library incubated with the L-rG4 target of interest, followed by binding, washing and elution of the library. Enriched D-aptamer candidates are sequenced and structurally characterized. Then, the L-aptamer is synthesized and used for different applications. rG4-SELEX can be carried out by an experienced molecular biologist with a basic understanding of nucleic acids. The development of rG4-targeting L-RNA aptamers expands the current rG4 toolkit to explore innovative rG4-related applications, and opens new doors to discovering novel rG4 biology in the near future. The duration of each selection cycle as outlined in the protocol is ~2 d.

Insights into the secondary and tertiary structure of the Bovine Viral Diarrhea Virus Internal Ribosome Entry Site

The internal ribosome entry site (IRES) RNA of bovine viral diarrhoea virus (BVDV), an economically significant Pestivirus, is required for the cap-independent translation of viral genomic RNA. Thus, it is essential for viral replication and pathogenesis. We applied a combination of high-throughput biochemical RNA structure probing (SHAPE-MaP) and in silico modelling approaches to gain insight into the secondary and tertiary structures of BVDV IRES RNA. Our study demonstrated that BVDV IRES RNA in solution forms a modular architecture composed of three distinct structural domains (I-III). Two regions within domain III are represented in tertiary interactions to form an H-type pseudoknot. Computational modelling of the pseudoknot motif provided a fine-grained picture of the tertiary structure and local arrangement of helices in the BVDV IRES. Furthermore, comparative genomics and consensus structure predictions revealed that the pseudoknot is evolutionarily conserved among many Pestivirus species. These studies provide detailed insight into the structural arrangement of BVDV IRES RNA H-type pseudoknot and encompassing motifs that likely contribute to the optimal functionality of viral cap-independent translation element.

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.

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

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

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

Pseudoknot length modulates the folding, conformational dynamics, and robustness of Xrn1 resistance of flaviviral xrRNAs

To understand how RNA dynamics is regulated and connected to its function, we investigate the folding, conformational dynamics and robustness of Xrn1 resistance of a set of flaviviral xrRNAs using SAXS, smFRET and in vitro enzymatic assays. Flaviviral xrRNAs form discrete ring-like 3D structures, in which the length of a conserved long-range pseudoknot (PK2) ranges from 2 bp to 7 bp. We find that xrRNAs' folding, conformational dynamics and Xrn1 resistance are strongly correlated and highly Mg²⁺-dependent, furthermore, the Mg²⁺-dependence is modulated by PK2 length variations. xrRNAs with long PK2 require less Mg²⁺ to stabilize their folding, exhibit reduced conformational dynamics and strong Xrn1 resistance even at low Mg²⁺, and tolerate mutations at key tertiary motifs at high Mg²⁺, which generally are destructive to xrRNAs with short PK2. These results demonstrate an unusual regulatory mechanism of RNA dynamics providing insights into the functions and future biomedical applications of xrRNAs.

Prediction of RNA secondary structure including pseudoknots for long sequences

RNA structural elements called pseudoknots are involved in various biological phenomena including ribosomal frameshifts. Because it is infeasible to construct an efficiently computable secondary structure model including pseudoknots, secondary structure prediction methods considering pseudoknots are not yet widely available. We developed IPknot, which uses heuristics to speed up computations, but it has remained difficult to apply it to long sequences, such as messenger RNA and viral RNA, because it requires cubic computational time with respect to sequence length and has threshold parameters that need to be manually adjusted. Here, we propose an improvement of IPknot that enables calculation in linear time by employing the LinearPartition model and automatically selects the optimal threshold parameters based on the pseudo-expected accuracy. In addition, IPknot showed favorable prediction accuracy across a wide range of conditions in our exhaustive benchmarking, not only for single sequences but also for multiple alignments.

Mass Spectrometry of Nucleic Acid Noncovalent Complexes

Nucleic acids have been among the first targets for antitumor drugs and antibiotics. With the unveiling of new biological roles in regulation of gene expression, specific DNA and RNA structures have become very attractive targets, especially when the corresponding proteins are undruggable. Biophysical assays to assess target structure as well as ligand binding stoichiometry, affinity, specificity, and binding modes are part of the drug development process. Mass spectrometry offers unique advantages as a biophysical method owing to its ability to distinguish each stoichiometry present in a mixture. In addition, advanced mass spectrometry approaches (reactive probing, fragmentation techniques, ion mobility spectrometry, ion spectroscopy) provide more detailed information on the complexes. Here, we review the fundamentals of mass spectrometry and all its particularities when studying noncovalent nucleic acid structures, and then review what has been learned thanks to mass spectrometry on nucleic acid structures, self-assemblies (e.g., duplexes or G-quadruplexes), and their complexes with ligands.

Motif mutations in pseudoknot stem I upstream of start codon in Senecavirus A genome: Impacts on activity of viral IRES and on rescue of recombinant virus

Senecavirus A (SVA), formerly known as Seneca Valley virus, is classified into the genus Senecavirus in the family Picornaviridae. Mature virion harbors an approximately 7 300-nt-long, positive-sense, and single-stranded RNA genome, which contains 5' and 3' untranslated regions (UTRs). Internal ribosome entry site (IRES) is identified in the SVA 5' UTR, and includes a RNA pseudoknot upstream of the start codon. This pseudoknot contains two stem structures, pseudoknot stem I and II (PKS-I and -II). The PKS-I is composed of two base-paired motifs (PKS-Ia and -Ib), between which there is an unpaired spacing (UpS). We reported previously that motif mutation in the PKS-II did not abolish the IRES activity, but interfered with SVA recovery from cDNA clone. In this study, we constructed five SVA minigenomes with point mutations in the PKS-I motif. Dual-luciferase reporter assay showed that motif mutations in PKS-I did not significantly interfere with the IRES activity to initiate protein expression. Correspondingly, we constructed five SVA cDNA clones with point mutations in the PKS-I motif. These genetically modified cDNA clones were separately transfected into BSR-T7/5 cells in attempting to rescue competent SVAs. However, only two viruses, namely PKS-Ia- and UpS-mutated recombinants, could be recovered from their individual cDNA clones. It can be concluded that the PKS-Ib is indispensable for viral growth.

Stem II-disrupting pseudoknot does not abolish ability of Senecavirus A IRES to initiate protein expression, but inhibits recovery of virus from cDNA clone

Senecavirus A (SVA) is classified into the genus Senecavirus in the family Picornaviridae. Its genome is a positive-sense, single-stranded and nonsegmented RNA, approximately 7300 nucleotides in length. A picornaviral genome is essentially an mRNA, which, albeit unmodified with 5' cap structure, can still initiate protein expression by the internal ribosome entry site (IRES). The SVA genome contains a hepatitis C virus-like IRES, in which a pseudoknot structure plays an important role in initiating protein expression. In this study, we constructed a set of SVA (CH-LX-01-2016 strain) minigenomes with all combinations of point mutations in its pseudoknot stem II (PKS-II). The results showed that any combination of point mutations could not significantly interfere with the IRES to initiate protein expression. Further, we constructed a full-length SVA cDNA clone, in which the PKS-II-forming cDNA motif was subjected to site-directed mutagenesis for totally disrupting the PKS-II formation in IRES. Such a modified SVA cDNA clone was transfected into BSR-T7/5 cells, consequently demonstrating that the PKS-II-disrupting IRES interfered neither with protein expression nor with antigenome replication, whereas a competent SVA could not be rescued from the cDNA clone. It was speculated that the mutated motif possibly disrupted a packaging signal for virion assembly, therefore causing the failure of SVA rescue.

Molecular motor traffic with a slow binding site

We discuss how the presence of a slow binding site in molecular motor traffic gives rise to defect-induced "traffic jams" that have properties different from those of the well-studied boundary-induced jams that originate from an imbalance between initiation and termination. To this end we analyze in detail the stationary distribution of a lattice gas model for traffic of molecular motors with a defect. In particular, we obtain analytically the exact spatial distribution of motors, the probability distribution of the random position of the molecular traffic jam and we report unexpected spatial anticorrelations between local molecular motor densities near the defect.

Programmed -1 Ribosomal Frameshifting in coronaviruses: A therapeutic target

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

ATTfold: RNA Secondary Structure Prediction With Pseudoknots Based on Attention Mechanism

Accurate RNA secondary structure information is the cornerstone of gene function research and RNA tertiary structure prediction. However, most traditional RNA secondary structure prediction algorithms are based on the dynamic programming (DP) algorithm, according to the minimum free energy theory, with both hard and soft constraints. The accuracy is particularly dependent on the accuracy of soft constraints (from experimental data like chemical and enzyme detection). With the elongation of the RNA sequence, the time complexity of DP-based algorithms will increase geometrically, as a result, they are not good at coping with relatively long sequences. Furthermore, due to the complexity of the pseudoknots structure, the secondary structure prediction method, based on traditional algorithms, has great defects which cannot predict the secondary structure with pseudoknots well. Therefore, few algorithms have been available for pseudoknots prediction in the past. The ATTfold algorithm proposed in this article is a deep learning algorithm based on an attention mechanism. It analyzes the global information of the RNA sequence via the characteristics of the attention mechanism, focuses on the correlation between paired bases, and solves the problem of long sequence prediction. Moreover, this algorithm also extracts the effective multi-dimensional features from a great number of RNA sequences and structure information, by combining the exclusive hard constraints of RNA secondary structure. Hence, it accurately determines the pairing position of each base, and obtains the real and effective RNA secondary structure, including pseudoknots. Finally, after training the ATTfold algorithm model through tens of thousands of RNA sequences and their real secondary structures, this algorithm was compared with four classic RNA secondary structure prediction algorithms. The results show that our algorithm significantly outperforms others and more accurately showed the secondary structure of RNA. As the data in RNA sequence databases increase, our deep learning-based algorithm will have superior performance. In the future, this kind of algorithm will be more indispensable.

Tertiary Base Triple Formation in the SRV-1 Frameshifting Pseudoknot Stabilizes Secondary Structure Components

Minor-groove base triples formed between stem 1 and loop 2 of the simian retrovirus type 1 (SRV-1) mRNA frameshifting pseudoknot are essential in stimulating -1 ribosomal frameshifting. How tertiary base triple formation affects the local stabilities of secondary structures (stem 1 and stem 2) and thus ribosomal frameshifting efficiency is not well understood. We made a short peptide nucleic acid (PNA) that is expected to invade stem 1 of the SRV-1 pseudoknot by PNA-RNA duplex formation to mimic the stem 1 unwinding process by a translating ribosome. In addition, we used a PNA for invading stem 2 in the SRV-1 pseudoknot. Our nondenaturing polyacrylamide gel electrophoresis data for the binding of PNA to the SRV-1 pseudoknot and mutants reveal that mutations in loop 2 disrupting base triple formation between loop 2 and stem 1 in the SRV-1 pseudoknot result in enhanced invasion by both PNAs. Our data suggest that tertiary stem 1-loop 2 base triple interactions in the SRV-1 pseudoknot can stabilize both of the secondary structural components, stem 1 and stem 2. Stem 2 stability is thus coupled to the structural stability of stem 1-loop 2 base triples, mediated through a long-range effect. The apparent dissociation constants of both PNAs are positively correlated with the pseudoknot mechanical stabilities and frameshifting efficiencies. The relatively simple PNA local invasion experiment may be used to characterize the energetic contribution of tertiary interactions and ligand binding in many other RNA and DNA structures.

Co-infection of a novel fusagravirus and a partitivirus in a Korean isolate of Rosellinia necatrix KACC40168

We report here the presence of dsRNA mycoviruses in a Korean isolate of Rosellinia necatrix. A multiple band pattern of double-stranded RNA (dsRNA) from R. necatrix suggested mixed mycovirus infection. Next-generation sequencing analysis of purified dsRNAs indicated the presence of two dsRNA mycoviruses related to the members of families "Fusagraviridae" (proposed) and Partitiviridae. The first dsRNA virus revealed that the complete genome sequence was 8868 bp in size and contained two large open reading frames (ORFs 1 and 2), overlapped by 22 bp containing a canonical (- 1) slippery heptanucelotide sequence of UUUAAAC. The deduced amino acid sequence of ORF1 and ORF2 showed highest similarity to the hypothetical protein and RNA-dependent RNA polymerase (RdRp) of Rosellinia necatrix fusagravirus 3 (RnFGV3). Phylogenetic analysis showed that this dsRNA virus clustered with RnFGV3 and other fusagraviruses. Gene organization, sequence similarity, and phylogenetic analysis indicate that this virus seems to belong to a novel species of "Fusagraviridae", which we have named Rosellinia necatrix fusagravirus 4. The second virus has two dsRNA segments with sizes of 1907 bp and 1918 bp, each of which encoded a single ORF showing highest similarity to the RdRp and capsid protein of known members of Partitiviridae. Evaluation of genome structure, sequence similarity, and phylogeny indicate this to be a new member of the genus Alphapartitivirus in the family Partitiviridae, hereafter designated as Rosellinia necatrix partitivirus 26. This is the first report of the presence of a fusagravirus in an Asian R. necatrix isolate and of its mixed infection with a partitivirus.

Characterization of a novel dsRNA mycovirus of Trichoderma atroviride NFCF377 reveals a member of "Fusagraviridae" with changes in antifungal activity of the host fungus

Trichoderma atroviride is a common fungus found in various ecosystems that shows mycoparasitic ability on other fungi. A novel dsRNA virus was isolated from T. atroviride NFCF377 strain and its molecular features were analyzed. The viral genome consists of a single segmented double-stranded RNA and is 9,584 bp in length, with two discontinuous open reading frames (ORF1 and ORF2). A mycoviral structural protein and an RNA-dependent RNA polymerase (RdRp) are encoded by ORF1 and ORF2, respectively, between which is found a canonical shifty heptameric signal motif (AAAAAAC) followed by an RNA pseudoknot. Analysis of sequence similarity and phylogeny showed that it is closely related to members of the proposed family "Fusagraviridae", with a highest similarity to the Trichoderma atroviride mycovirus 1 (TaMV1). Although the sequence similarity of deduced amino acid to TaMV1 was evident, sequence deviations were distinctive at untranslated regions (UTRs) due to the extended size. Thus, we inferred this dsRNA to be a different strain of Trichoderma atroviride mycovirus 1 (TaMV1-NFCF377). Electron microscopy image exhibited an icosahedral viral particle of 40 nm diameter. Virus-cured isogenic isolates were generated and no differences in growth rate, colony morphology, or conidia production were observed between virus-infected and virus-cured strains. However, culture filtrates of TaMV1-NFCF377-infected strain showed enhanced antifungal activity against the plant pathogen Rhizoctonia solani but not to edible mushroom Pleurotus ostreatus. These results suggested that TaMV1-NFCF377 affected the metabolism of the fungal host to potentiate antifungal compounds against a plant pahogen, but this enhanced antifungal activity appeared to be species-specific.