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

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

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

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

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

Phase Space Invaders' podcast episode with Tamar Schlick: a trajectory from mathematics to biology

We present a transcript of the Phase Space Invaders podcast interview, with Tamar Schlick interviewed by Miłosz Wieczór. The conversation covers topics in computational biophysics and beyond: DNA and RNA research from genome organization to viral RNA frameshifting, transitioning from applied math to biology, developing algorithms and their utility in molecular dynamics and complex multiscale systems, the role of computers in biophysical research, writing reviews and books, collaborating in science, and using long-distance running as a template for building supportive communities.

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

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

High-throughput sequencing: a breakthrough in molecular diagnosis for precision medicine

High-resolution insights into the nucleotide arrangement within an organism's genome are pivotal for deciphering its genetic composition, function, and evolutionary trajectory. Over the years, nucleic acid sequencing has been instrumental in driving significant advancements in genomics and molecular biology. The advent of high-throughput or next-generation sequencing (NGS) technologies has revolutionized whole genome sequencing, revealing novel and intriguing features of genomes, such as single nucleotide polymorphisms and lethal mutations in both coding and non-coding regions. These platforms provide a practical approach to comprehensively identifying and analyzing whole genomes with remarkable throughput, accuracy, and scalability within a short time frame. The resulting data holds immense potential for enhancing healthcare systems, developing novel and personalized therapies, and preparing for future pandemics and outbreaks. Given the wide array of available high-throughput sequencing platforms, selecting the appropriate technology based on specific needs is crucial. However, there is limited information regarding sample preparation, sequencing principles, and output data to facilitate a comparative evaluation of these platforms. This review details various NGS technologies and approaches, examining their advantages, limitations, and future potential. Despite being in their early stages and facing challenges, ongoing advancements in NGS are expected to yield significant future benefits.

How Large is the Universe of RNA-Like Motifs? A Clustering Analysis of RNA Graph Motifs Using Topological Descriptors

Identifying novel and functional RNA structures remains a significant challenge in RNA motif design and is crucial for developing RNA-based therapeutics. Here we introduce a computational topology-based approach with unsupervised machine-learning algorithms to estimate the database size and content of RNA-like graph topologies. Specifically, we apply graph theory enumeration to generate all 110,667 possible 2D dual graphs for vertex numbers ranging from 2 to 9. Among them, only 0.11% (121 dual graphs) correspond to approximately 200,000 known RNA atomic fragments/substructures (collected in 2021) using the RNA-as-Graphs (RAG) mapping method. The remaining 99.89% of the dual graphs may be RNA-like or non-RNA-like. To determine which dual graphs in the 99.89% hypothetical set are more likely to be associated with RNA structures, we apply computational topology descriptors using the Persistent Spectral Graphs (PSG) method to characterize each graph using 19 PSG-based features and use clustering algorithms that partition all possible dual graphs into two clusters. The cluster with the higher percentage of known dual graphs for RNA is defined as the "RNA-like" cluster, while the other is considered as "non-RNA-like". The distance of each dual graph to the center of the RNA-like cluster represents the likelihood of it belonging to RNA structures. From validation, our PSG-based RNA-like cluster includes 97.3% of the 121 known RNA dual graphs, suggesting good performance. Furthermore, 46.017% of the hypothetical RNAs are predicted to be RNA-like. Among the top 15 graphs identified as high-likelihood candidates for novel RNA motifs, 4 were confirmed from the RNA dataset collected in 2022. Significantly, we observe that all the top 15 RNA-like dual graphs can be separated into multiple subgraphs, whereas the top 15 non-RNA-like dual graphs tend not to have any subgraphs (subgraphs preserve pseudoknots and junctions). Moreover, a significant topological difference between top RNA-like and non-RNA-like graphs is evident when comparing their topological features (e.g. Betti-0 and Betti-1 numbers). These findings provide valuable insights into the size of the RNA motif universe and RNA design strategies, offering a novel framework for predicting RNA graph topologies and guiding the discovery of novel RNA motifs, perhaps anti-viral therapeutics by subgraph assembly.

CParty: hierarchically constrained partition function of RNA pseudoknots

MOTIVATION

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

RESULTS

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

AVAILABILITY AND IMPLEMENTATION

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

Tree polynomials identify a link between co-transcriptional R-loops and nascent RNA folding

R-loops are a class of non-canonical nucleic acid structures that typically form during transcription when the nascent RNA hybridizes the DNA template strand, leaving the non-template DNA strand unpaired. These structures are abundant in nature and play important physiological and pathological roles. Recent research shows that DNA sequence and topology affect R-loops, yet it remains unclear how these and other factors contribute to R-loop formation. In this work, we investigate the link between nascent RNA folding and the formation of R-loops. We introduce tree-polynomials, a new class of representations of RNA secondary structures. A tree-polynomial representation consists of a rooted tree associated with an RNA secondary structure together with a polynomial that is uniquely identified with the rooted tree. Tree-polynomials enable accurate, interpretable and efficient data analysis of RNA secondary structures without pseudoknots. We develop a computational pipeline for investigating and predicting R-loop formation from a genomic sequence. The pipeline obtains nascent RNA secondary structures from a co-transcriptional RNA folding software, and computes the tree-polynomial representations of the structures. By applying this pipeline to plasmid sequences that contain R-loop forming genes, we establish a strong correlation between the coefficient sums of tree-polynomials and the experimental probability of R-loop formation. Such strong correlation indicates that the pipeline can be used for accurate R-loop prediction. Furthermore, the interpretability of tree-polynomials allows us to characterize the features of RNA secondary structure associated with R-loop formation. In particular, we identify that branches with short stems separated by bulges and interior loops are associated with R-loops.

Multisacle Jones Polynomial and Persistent Jones Polynomial for Knot Data Analysis

Many structures in science, engineering, and art can be viewed as curves in 3-space. The entanglement of these curves plays a crucial role in determining the functionality and physical properties of materials. Many concepts in knot theory provide theoretical tools to explore the complexity and entanglement of curves in 3-space. However, classical knot theory primarily focuses on global topological properties and lacks the consideration of local structural information, which is critical in practical applications. In this work, two localized models based on the Jones polynomial, namely the multiscale Jones polynomial and the persistent Jones polynomial, are proposed. The stability of these models, especially the insensitivity of the multiscale and persistent Jones polynomial models to small perturbations in curve collections, is analyzed, thus ensuring their robustness for real-world applications.

The role of structure in regulatory RNA elements

Regulatory RNA elements fulfill functions such as translational regulation, control of transcript levels, and regulation of viral genome replication. Trans-acting factors (i.e., RNA-binding proteins) bind the so-called cis elements and confer functionality to the complex. The specificity during protein-RNA complex (RNP) formation often exploits the structural plasticity of RNA. Functional integrity of cis-trans pairs depends on the availability of properly folded RNA elements, and RNA conformational transitions can cause diseases. Knowledge of RNA structure and the conformational space is needed for understanding complex formation and deducing functional effects. However, structure determination of RNAs under in vivo conditions remains challenging. This review provides an overview of structured eukaryotic and viral RNA cis elements and discusses the effect of RNA structural equilibria on RNP formation. We showcase implications of RNA structural changes for diseases, outline strategies for RNA structure-based drug targeting, and summarize the methodological toolbox for deciphering RNA structures.

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.

Knot data analysis using multiscale Gauss link integral

In the past decade, topological data analysis has emerged as a powerful algebraic topology approach in data science. Although knot theory and related subjects are a focus of study in mathematics, their success in practical applications is quite limited due to the lack of localization and quantization. We address these challenges by introducing knot data analysis (KDA), a paradigm that incorporates curve segmentation and multiscale analysis into the Gauss link integral. The resulting multiscale Gauss link integral (mGLI) recovers the global topological properties of knots and links at an appropriate scale and offers a multiscale geometric topology approach to capture the local structures and connectivities in data. By integration with machine learning or deep learning, the proposed mGLI significantly outperforms other state-of-the-art methods across various benchmark problems in 13 intricately complex biological datasets, including protein flexibility analysis, protein-ligand interactions, human Ether-à-go-go-Related Gene potassium channel blockade screening, and quantitative toxicity assessment. Our KDA opens a research area-knot deep learning-in data science.

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.

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.

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 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.

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.

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.

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.

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

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

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

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

SARS-CoV-2 external structures interacting with nanospheres using docking and molecular dynamics

Coronavirus is caused by the SARS-CoV-2 virus has shown rapid proliferation and scarcity of treatments with proven effectiveness. In this way, we simulated the hospitalization of carbon nanospheres, with external active sites of the SARS-CoV-2 virus (M-Pro, S-Gly and E-Pro), which can be adsorbed or inactivated when interacting with the nanospheres. The computational procedures performed in this work were developed with the SwissDock server for molecular docking and the GROMACS software for molecular dynamics, making it possible to extract relevant data on affinity energy, distance between molecules, free Gibbs energy and mean square deviation of atomic positions, surface area accessible to solvents. Molecular docking indicates that all ligands have an affinity for the receptor's active sites. The nanospheres interact favorably with all proteins, showing promising results, especially C60, which presented the best affinity energy and RMSD values ​​for all protein macromolecules investigated. The C60 with E-Pro exhibited the highest affinity energy of -9.361 kcal/mol, demonstrating stability in both molecular docking and molecular dynamics simulations. Our RMSD calculations indicated that the nanospheres remained predominantly stable, fluctuating within a range of 2 to 3 Å. Additionally, the analysis of other structures yielded promising results that hold potential for application in other proteases.Communicated by Ramaswamy H. Sarma.

Interactions between carbon nanotubes and external structures of SARS-CoV-2 using molecular docking and molecular dynamics

Molecular modeling techniques are used to describe the process of interaction between nanotubes and the main structures of the Covid-19 virus: the envelope protein, the main protease, and the Spike glycoprotein. Molecular docking studies show that the ligands have interaction characteristics capable of adsorbing the structures. Molecular dynamics simulations provide information on the mean squared deviation of atomic positions ​​between 0.5 and 3.0 Å. The Gibbs free energy model and solvent accessible surface area approaches are used. Through the results obtained through molecular dynamics simulations, it is noted that the zig-zag nanotube prefers to interact with E-pro, M-pro, and S-gly, respectively. Molecular couplings and free energy showed that the S-gly active site residues strongly interact with zigzag, chiral, and armchair nanotubes, in this order. The interactions demonstrated in this manuscript may predict some promising candidates for virus antagonists, which may be confirmed through experimental approaches.

Structural Modifications and Novel Protein-Binding Sites in Pre-miR-675-Explaining Its Regulatory Mechanism in Carcinogenesis

Pre-miR-675 is a microRNA expressed from the exon 1 of H19 long noncoding RNA, and the atypical expression of pre-miR-675 has been linked with several diseases and disorders including cancer. To execute its function inside the cell, pre-miR-675 is folded into a particular conformation, which aids in its interaction with several other biological molecules. However, the exact folding dynamics of pre-miR-675 and its protein-binding motifs are currently unknown. Moreover, how H19 lncRNA and pre-miR-675 crosstalk and modulate each other's activities is also unclear. The detailed structural analysis of pre-miR-675 in this study determines its earlier unknown conformation and identifies novel protein-binding sites on pre-miR-675, thus making it an excellent therapeutic target against cancer. Co-folding analysis between H19 lncRNA and pre-miR-675 determine structural transformations in pre-miR-675, thus describing the earlier unknown mechanism of interaction between these two molecules. Comprehensively, this study details the conformation of pre-miR-675 and its protein-binding sites and explains its relationship with H19 lncRNA, which can be interpreted to understand the role of pre-miR-675 in the development and progression of tumorigenesis and designing new therapeutics against cancers.

RNA folding pathways from all-atom simulations with a variationally improved history-dependent bias

Atomically detailed simulations of RNA folding have proven very challenging in view of the difficulties of developing realistic force fields and the intrinsic computational complexity of sampling rare conformational transitions. As a step forward in tackling these issues, we extend to RNA an enhanced path-sampling method previously successfully applied to proteins. In this scheme, the information about the RNA's native structure is harnessed by a soft history-dependent biasing force promoting the generation of productive folding trajectories in an all-atom force field with explicit solvent. A rigorous variational principle is then applied to minimize the effect of the bias. Here, we report on an application of this method to RNA molecules from 20 to 47 nucleotides long and increasing topological complexity. By comparison with analog simulations performed on small proteins with similar size and architecture, we show that the RNA folding landscape is significantly more frustrated, even for relatively small chains with a simple topology. The predicted RNA folding mechanisms are found to be consistent with the available experiments and some of the existing coarse-grained models. Due to its computational performance, this scheme provides a promising platform to efficiently gather atomistic RNA folding trajectories, thus retain the information about the chemical composition of the sequence.

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

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

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

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

The Functional RNA Identification (FRID) Pipeline: Identification of Potential Pseudoknot-Containing RNA Elements as Therapeutic Targets for SARS-CoV-2

The COVID-19 pandemic persists despite the development of effective vaccines. As such, it remains crucial to identify new targets for antiviral therapies. The causative virus of COVID-19, SARS-CoV-2, is a positive-sense RNA virus with RNA structures that could serve as therapeutic targets. One such RNA with established function is the frameshift stimulatory element (FSE), which promotes programmed ribosomal frameshifting. To accelerate identification of additional functional RNA elements, we introduce a novel computational approach termed the Functional RNA Identification (FRID) pipeline. The guiding principle of our pipeline, which uses established component programs as well as customized component programs, is that functional RNA elements have conserved secondary and pseudoknot structures that facilitate function. To assess the presence and conservation of putative functional RNA elements in SARS-CoV-2, we compared over 6,000 SARS-CoV-2 genomic isolates. We identified 22 functional RNA elements from the SARS-CoV-2 genome, 14 of which have conserved pseudoknots and serve as potential targets for small molecule or antisense oligonucleotide therapeutics. The FRID pipeline is general and can be applied to identify pseudoknotted RNAs for targeted therapeutics in genomes or transcriptomes from any virus or organism.

Probing the dynamic RNA structurome and its functions

RNA is a key regulator of almost every cellular process, and the structures adopted by RNA molecules are thought to be central to their functions. The recent fast-paced evolution of high-throughput sequencing-based RNA structure mapping methods has enabled the rapid in vivo structural interrogation of entire cellular transcriptomes. Collectively, these studies are shedding new light on the long underestimated complexity of the structural organization of the transcriptome - the RNA structurome. Moreover, recent analyses are challenging the view that the RNA structurome is a static entity by revealing how RNA molecules establish intricate networks of alternative intramolecular and intermolecular interactions and that these ensembles of RNA structures are dynamically regulated to finely tune RNA functions in living cells. This new understanding of how RNA can shape cell phenotypes has important implications for the development of RNA-targeted therapeutic strategies.

Shapify: Paths to SARS-CoV-2 frameshifting pseudoknot

Multiple coronaviruses including MERS-CoV causing Middle East Respiratory Syndrome, SARS-CoV causing SARS, and SARS-CoV-2 causing COVID-19, use a mechanism known as -1 programmed ribosomal frameshifting (-1 PRF) to replicate. SARS-CoV-2 possesses a unique RNA pseudoknotted structure that stimulates -1 PRF. Targeting -1 PRF in SARS-CoV-2 to impair viral replication can improve patients' prognoses. Crucial to developing these therapies is understanding the structure of the SARS-CoV-2 -1 PRF pseudoknot. Our goal is to expand knowledge of -1 PRF structural conformations. Following a structural alignment approach, we identify similarities in -1 PRF pseudoknots of SARS-CoV-2, SARS-CoV, and MERS-CoV. We provide in-depth analysis of the SARS-CoV-2 and MERS-CoV -1 PRF pseudoknots, including reference and noteworthy mutated sequences. To better understand the impact of mutations, we provide insight on -1 PRF pseudoknot sequence mutations and their effect on resulting structures. We introduce Shapify, a novel algorithm that given an RNA sequence incorporates structural reactivity (SHAPE) data and partial structure information to output an RNA secondary structure prediction within a biologically sound hierarchical folding approach. Shapify enhances our understanding of SARS-CoV-2 -1 PRF pseudoknot conformations by providing energetically favourable predictions that are relevant to structure-function and may correlate with -1 PRF efficiency. Applied to the SARS-CoV-2 -1 PRF pseudoknot, Shapify unveils previously unknown paths from initial stems to pseudoknotted structures. By contextualizing our work with available experimental data, our structure predictions motivate future RNA structure-function research and can aid 3-D modeling of pseudoknots.

RNA structure-altering mutations underlying positive selection on Spike protein reveal novel putative signatures to trace crossing host-species barriers in Betacoronavirus

Similar to other RNA viruses, the emergence of Betacoronavirus relies on cross-species viral transmission, which requires careful health surveillance monitoring of protein-coding information as well as genome-wide analysis. Although the evolutionary jump from natural reservoirs to humans may be mainly traced-back by studying the effect that hotspot mutations have on viral proteins, it is largely unexplored if other impacts might emerge on the structured RNA genome of Betacoronavirus. In this survey, the protein-coding and viral genome architecture were simultaneously studied to uncover novel insights into cross-species horizontal transmission events. We analysed 1,252,952 viral genomes of SARS-CoV, MERS-CoV, and SARS-CoV-2 distributed across the world in bats, intermediate animals, and humans to build a new landscape of changes in the RNA viral genome. Phylogenetic analyses suggest that bat viruses are the most closely related to the time of most recent common ancestor of Betacoronavirus, and missense mutations in viral proteins, mainly in the S protein S1 subunit: SARS-CoV (G > T; A577S); MERS-CoV (C > T; S746R and C > T; N762A); and SARS-CoV-2 (A > G; D614G) appear to have driven viral diversification. We also found that codon sites under positive selection on S protein overlap with non-compensatory mutations that disrupt secondary RNA structures in the RNA genome complement. These findings provide pivotal factors that might be underlying the eventual jumping the species barrier from bats to intermediate hosts. Lastly, we discovered that nearly half of the Betacoronavirus genomes carry highly conserved RNA structures, and more than 90% of these RNA structures show negative selection signals, suggesting essential functions in the biology of Betacoronavirus that have not been investigated to date. Further research is needed on negatively selected RNA structures to scan for emerging functions like the potential of coding virus-derived small RNAs and to develop new candidate antiviral therapeutic strategies.

Brief considerations on targeting RNA with small molecules

For more than three decades, RNA has been known to be a relevant and attractive macromolecule to target but figuring out which RNA should be targeted and how remains challenging. Recent years have seen the confluence of approaches for screening, drug optimization, and target validation that have led to the approval of a few RNA-targeting therapeutics for clinical applications. This focused perspective aims to highlight - but not exhaustively review - key factors accounting for these successes while pointing at crucial aspects worth considering for further breakthroughs.

Cis-mediated interactions of the SARS-CoV-2 frameshift RNA alter its conformations and affect function

The RNA genome of SARS-CoV-2 contains a frameshift stimulatory element (FSE) that allows access to an alternative reading frame through -1 programmed ribosomal frameshifting (PRF). -1PRF in the 1a/1b gene is essential for efficient viral replication and transcription of the viral genome. -1PRF efficiency relies on the presence of conserved RNA elements within the FSE. One of these elements is a three-stemmed pseudoknot, although alternative folds of the frameshift site might have functional roles as well. Here, by complementing ensemble and single-molecule structural analysis of SARS-CoV-2 frameshift RNA variants with functional data, we reveal a conformational interplay of the 5' and 3' immediate regions with the FSE and show that the extended FSE exists in multiple conformations. Furthermore, limiting the base pairing of the FSE with neighboring nucleotides can favor or impair the formation of the alternative folds, including the pseudoknot. Our results demonstrate that co-existing RNA structures can function together to fine-tune SARS-CoV-2 gene expression, which will aid efforts to design specific inhibitors of viral frameshifting.

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

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

Biomotors, viral assembly, and RNA nanobiotechnology: Current achievements and future directions

The International Society of RNA Nanotechnology and Nanomedicine (ISRNN) serves to further the development of a wide variety of functional nucleic acids and other related nanotechnology platforms. To aid in the dissemination of the most recent advancements, a biennial discussion focused on biomotors, viral assembly, and RNA nanobiotechnology has been established where international experts in interdisciplinary fields such as structural biology, biophysical chemistry, nanotechnology, cell and cancer biology, and pharmacology share their latest accomplishments and future perspectives. The results summarized here highlight advancements in our understanding of viral biology and the structure-function relationship of frame-shifting elements in genomic viral RNA, improvements in the predictions of SHAPE analysis of 3D RNA structures, and the understanding of dynamic RNA structures through a variety of experimental and computational means. Additionally, recent advances in the drug delivery, vaccine design, nanopore technologies, biomotor and biomachine development, DNA packaging, RNA nanotechnology, and drug delivery are included in this critical review. We emphasize some of the novel accomplishments, major discussion topics, and present current challenges and perspectives of these emerging fields.

Rotavirus RNA chaperone mediates global transcriptome-wide increase in RNA backbone flexibility

Due to genome segmentation, rotaviruses must co-package eleven distinct genomic RNAs. The packaging is mediated by virus-encoded RNA chaperones, such as the rotavirus NSP2 protein. While the activities of distinct RNA chaperones are well studied on smaller RNAs, little is known about their global effect on the entire viral transcriptome. Here, we used Selective 2'-hydroxyl Acylation Analyzed by Primer Extension and Mutational Profiling (SHAPE-MaP) to examine the secondary structure of the rotavirus transcriptome in the presence of increasing amounts of NSP2. SHAPE-MaP data reveals that despite the well-documented helix-unwinding activity of NSP2 in vitro, its incubation with cognate rotavirus transcripts does not induce a significant change in the SHAPE reactivities. However, a quantitative analysis of mutation rates measured by mutational profiling reveals a global 5-fold rate increase in the presence of NSP2. We demonstrate that the normalization procedure used in deriving SHAPE reactivities from mutation rates can mask an important global effect of an RNA chaperone. Analysis of the mutation rates reveals a larger effect on stems rather than loops. Together, these data provide the first experimentally derived secondary structure model of the rotavirus transcriptome and reveal that NSP2 acts by globally increasing RNA backbone flexibility in a concentration-dependent manner.

RNA-As-Graphs Motif Atlas—Dual Graph Library of RNA Modules and Viral Frameshifting-Element Applications

RNA motif classification is important for understanding structure/function connections and building phylogenetic relationships. Using our coarse-grained RNA-As-Graphs (RAG) representations, we identify recurrent dual graph motifs in experimentally solved RNA structures based on an improved search algorithm that finds and ranks independent RNA substructures. Our expanded list of 183 existing dual graph motifs reveals five common motifs found in transfer RNA, riboswitch, and ribosomal 5S RNA components. Moreover, we identify three motifs for available viral frameshifting RNA elements, suggesting a correlation between viral structural complexity and frameshifting efficiency. We further partition the RNA substructures into 1844 distinct submotifs, with pseudoknots and junctions retained intact. Common modules are internal loops and three-way junctions, and three submotifs are associated with riboswitches that bind nucleotides, ions, and signaling molecules. Together, our library of existing RNA motifs and submotifs adds to the growing universe of RNA modules, and provides a resource of structures and substructures for novel RNA design.

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

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

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

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

RNAspider: a webserver to analyze entanglements in RNA 3D structures

Advances in experimental and computational techniques enable the exploration of large and complex RNA 3D structures. These, in turn, reveal previously unstudied properties and motifs not characteristic for small molecules with simple architectures. Examples include entanglements of structural elements in RNA molecules and knot-like folds discovered, among others, in the genomes of RNA viruses. Recently, we presented the first classification of entanglements, determined by their topology and the type of entangled structural elements. Here, we introduce RNAspider - a web server to automatically identify, classify, and visualize primary and higher-order entanglements in RNA tertiary structures. The program applies to evaluate RNA 3D models obtained experimentally or by computational prediction. It supports the analysis of uncommon topologies in the pseudoknotted RNA structures. RNAspider is implemented as a publicly available tool with a user-friendly interface and can be freely accessed at https://rnaspider.cs.put.poznan.pl/.

Translational Control of COVID-19 and Its Therapeutic Implication

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

Thinking Outside the Frame: Impacting Genomes Capacity by Programmed Ribosomal Frameshifting

Translation facilitates the transfer of the genetic information stored in the genome via messenger RNAs to a functional protein and is therefore one of the most fundamental cellular processes. Programmed ribosomal frameshifting is a ubiquitous alternative translation event that is extensively used by viruses to regulate gene expression from overlapping open reading frames in a controlled manner. Recent technical advances in the translation field enabled the identification of precise mechanisms as to how and when ribosomes change the reading frame on mRNAs containing cis-acting signals. Several studies began also to illustrate that trans-acting RNA modulators can adjust the timing and efficiency of frameshifting illuminating that frameshifting can be a dynamically regulated process in cells. Here, we intend to summarize these new findings and emphasize how it fits in our current understanding of PRF mechanisms as previously described.

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

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

Expanding the known structure space for RNA binding: a test of 2,5-diketopiperazine

As the importance of RNA as a therapeutic target has become increasingly recognized, the need for new chemotypes able to bind RNA has grown in significance. We hypothesized that diketopiperazines (DKPs), common substructures in natural products and protein-targeting therapeutic agents, could serve as effective scaffolds for targeting RNA. To confirm this hypothesis, we designed and synthesized two analogs, one incorporating a DKP and one not, of compounds previously demonstrated to bind an RNA critical to the life cycle of HIV-1 with high affinity and specificity. Prior to compound synthesis, calculations employing density functional methods and molecular mechanics conformational searches were used to confirm that the DKP could present functionality in a similar (albeit not identical) orientation to the non DKP-containing compound. We found that both the DKP-containing and parent compound had similar affinities to the target RNA as measured by surface plasmon resonance (SPR). Both compounds were found to have modest but equal anti-HIV activity. These results establish the feasibility of using DKPs to target RNA.

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.

Targeting Ribosomal Frameshifting as an Antiviral Strategy: From HIV-1 to SARS-CoV-2

Treatment of HIV-1 has largely involved targeting viral enzymes using a cocktail of inhibitors. However, resistance to these inhibitors and toxicity in the long term have pushed the field to identify new therapeutic targets. To that end, -1 programmed ribosomal frameshifting (-1 PRF) has gained attention as a potential node for therapeutic intervention. In this process, a ribosome moves one nucleotide backward in the course of translating a mRNA, revealing a new reading frame for protein synthesis. In HIV-1, -1 PRF allows the virus to regulate the ratios of enzymatic and structural proteins as needed for correct viral particle assembly. Two RNA structural elements are central to -1 PRF in HIV: a slippery sequence and a highly conserved stable hairpin called the HIV-1 frameshifting stimulatory signal (FSS). Dysregulation of -1 PRF is deleterious for the virus. Thus, -1 PRF is an attractive target for new antiviral development. It is important to note that HIV-1 is not the only virus exploiting -1 PRF for regulating production of its proteins. Coronaviruses, including the COVID-19 pandemic virus SARS-CoV-2, also rely on -1 PRF. In SARS-CoV-2 and other coronaviruses, -1 PRF is required for synthesis of RNA-dependent RNA polymerase and several other nonstructural proteins. Coronaviruses employ a more complex RNA structural element for regulating -1 PRF called a pseudoknot. The purpose of this Account is primarily to review the development of molecules targeting HIV-1 -1 PRF. These approaches are case studies illustrating how the entire pipeline from screening to the generation of high-affinity leads might be implemented. We consider both target-based and function-based screening, with a particular focus on our group's approach beginning with a resin-bound dynamic combinatorial library (RBDCL) screen. We then used rational design approaches to optimize binding affinity, selectivity, and cellular bioavailability. Our tactic is, to the best of our knowledge, the only study resulting in compounds that bind specifically to the HIV-1 FSS RNA and reduce infectivity of laboratory and drug-resistant strains of HIV-1 in human cells. Lessons learned from strategies targeting -1 PRF HIV-1 might provide solutions in the development of antivirals in areas of unmet medical need. This includes the development of new frameshift-altering therapies for SARS-CoV-2, approaches to which are very recently beginning to appear.

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

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

A MULTISCALE VISION-ILLUSTRATIVE APPLICATIONS FROM BIOLOGY TO ENGINEERING

Modeling and simulation have quickly become equivalent pillars of research along with traditional theory and experimentation. The growing realization that most complex phenomena of interest span many orders of spatial and temporal scales has led to an exponential rise in the development and application of multiscale modeling and simulation over the past two decades. In this perspective, the associate editors of the International Journal for Multiscale Computational Engineering and their co-workers illustrate current applications in their respective fields spanning biomolecular structure and dynamics, civil engineering and materials science, computational mechanics, aerospace and mechanical engineering, and more. Such applications are highly tailored, exploit the latest and ever-evolving advances in both computer hardware and software, and contribute significantly to science, technology, and medical challenges in the 21st century.