Structure-altering mutations of the SARS-CoV-2 frameshifting RNA element

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.

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.

Unwinding the SARS-CoV-2 Ribosomal Frameshifting Pseudoknot with LNA and G-Clamp-Modified Phosphorothioate Oligonucleotides Inhibits Viral Replication

Ribosomal frameshifting (RFS) at the slippery site of SARS-CoV-2 RNA is essential for the biosynthesis of the viral replication machinery. It requires the formation of a pseudoknot (PK) structure near the slippery site and can be inhibited by PK-disrupting oligonucleotide-based antivirals. We obtained and compared three types of such antiviral candidates, namely locked nucleic acids (LNA), LNA-DNA gapmers, and G-clamp-containing phosphorothioates (CPSs) complementary to PK stems. Using optical and electrophoretic methods, we showed that stem 2-targeting oligonucleotide analogs induced PK unfolding at nanomolar concentrations, and this effect was particularly pronounced in the case of LNA. For the leading PK-unfolding LNA and CPS oligonucleotide analogs, we also demonstrated dose-dependent RSF inhibition in dual luciferase assays (DLAs). Finally, we showed that the leading oligonucleotide analogs reduced SARS-CoV-2 replication at subtoxic concentrations in the nanomolar range in two human cell lines. Our findings highlight the promise of PK targeting, illustrate the advantages and limitations of various types of DNA modifications and may promote the future development of oligonucleotide-based antivirals.

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.

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.

Evolution of Sequence and Structure of SARS-CoV-2 Spike Protein: A Dynamic Perspective

Novel coronavirus (SARS-CoV-2) enters its host cell through a surface spike protein. The viral spike protein has undergone several modifications/mutations at the genomic level, through which it modulated its structure-function and passed through several variants of concern. Recent advances in high-resolution structure determination and multiscale imaging techniques, cost-effective next-generation sequencing, and development of new computational methods (including information theory, statistical methods, machine learning, and many other artificial intelligence-based techniques) have hugely contributed to the characterization of sequence, structure, function of spike proteins, and its different variants to understand viral pathogenesis, evolutions, and transmission. Laying on the foundation of the sequence-structure-function paradigm, this review summarizes not only the important findings on structure/function but also the structural dynamics of different spike components, highlighting the effects of mutations on them. As dynamic fluctuations of three-dimensional spike structure often provide important clues for functional modulation, quantifying time-dependent fluctuations of mutational events over spike structure and its genetic/amino acidic sequence helps identify alarming functional transitions having implications for enhanced fusogenicity and pathogenicity of the virus. Although these dynamic events are more difficult to capture than quantifying a static, average property, this review encompasses those challenging aspects of characterizing the evolutionary dynamics of spike sequence and structure and their implications for functions.

Predicting 3D structures and stabilities for complex RNA pseudoknots in ion solutions

RNA pseudoknots are a kind of important tertiary motif, and the structures and stabilities of pseudoknots are generally critical to the biological functions of RNAs with the motifs. In this work, we have carefully refined our previously developed coarse-grained model with salt effect through involving a new coarse-grained force field and a replica-exchange Monte Carlo algorithm, and employed the model to predict structures and stabilities of complex RNA pseudoknots in ion solutions beyond minimal H-type pseudoknots. Compared with available experimental data, the newly refined model can successfully predict 3D structures from sequences for the complex RNA pseudoknots including SARS-CoV-2 programming-1 ribosomal frameshifting element and Zika virus xrRNA, and can reliably predict the thermal stabilities of RNA pseudoknots with various sequences and lengths over broad ranges of monovalent/divalent salts. In addition, for complex pseudoknots including SARS-CoV-2 frameshifting element, our analyses show that their thermally unfolding pathways are mainly dependent on the relative stabilities of unfolded intermediate states, in analogy to those of minimal H-type pseudoknots.

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.

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.

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.

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.

Variant-Specific Analysis Reveals a Novel Long-Range RNA-RNA Interaction in SARS-CoV-2 Orf1a

Since the start of the COVID-19 pandemic, understanding the pathology of the SARS-CoV-2 RNA virus and its life cycle has been the priority of many researchers. Currently, new variants of the virus have emerged with various levels of pathogenicity and abundance within the human-host population. Although much of viral pathogenicity is attributed to the viral Spike protein’s binding affinity to human lung cells’ ACE2 receptor, comprehensive knowledge on the distinctive features of viral variants that might affect their life cycle and pathogenicity is yet to be attained. Recent in vivo studies into the RNA structure of the SARS-CoV-2 genome have revealed certain long-range RNA-RNA interactions. Using in silico predictions and a large population of SARS-CoV-2 sequences, we observed variant-specific evolutionary changes for certain long-range RRIs. We also found statistical evidence for the existence of one of the thermodynamic-based RRI predictions, namely Comp1, in the Beta variant sequences. A similar test that disregarded sequence variant information did not, however, lead to significant results. When performing population-based analyses, aggregate tests may fail to identify novel interactions due to variant-specific changes. Variant-specific analyses can result in de novo RRI identification.

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.

Continent-wide evolutionary trends of emerging SARS-CoV-2 variants: dynamic profiles from Alpha to Omicron

The ongoing SARS-CoV-2 evolution process has generated several variants due to its continuous mutations, making pandemics more critical. The present study illustrates SARS-CoV-2 evolution and its emerging mutations in five directions. First, the significant mutations in the genome and S-glycoprotein were analyzed in different variants. Three linear models were developed with the regression line to depict the mutational load for S-glycoprotein, total genome excluding S-glycoprotein, and whole genome. Second, the continent-wide evolution of SARS-CoV-2 and its variants with their clades and divergence were evaluated. It showed the region-wise evolution of the SARS-CoV-2 variants and their clustering event. The major clades for each variant were identified. One example is clade 21K, a major clade of the Omicron variant. Third, lineage dynamics and comparison between SARS-CoV-2 lineages across different countries are also illustrated, demonstrating dominant variants in various countries over time. Fourth, gene-wise mutation patterns and genetic variability of SARS-CoV-2 variants across various countries are illustrated. High mutation patterns were found in the ORF10, ORF6, S, and low mutation pattern E genes. Finally, emerging AA point mutations (T478K, L452R, N501Y, S477N, E484A, Q498R, and Y505H), their frequencies, and country-wise occurrence were identified, and the highest event of two mutations (T478K and L452R) was observed.

Finite element modeling of alpha-helices and tropocollagen molecules with reference to the spike of SARS-CoV-2

The newly developed finite element modeling at the atomic scale was used to predict the static and dynamic response of the alpha-helix (AH) and tropocollagen (TC) protein fragments, the main building blocks of the spike of the SARS-CoV-2. The geometry and morphology of the spike's stalk and its connection to the viral envelope were determined from the combination of most recent Molecular Dynamics simulation and images of Cryo-Electron microscope. The stiffness parameters of the covalent bonds in the main chain of the helix were taken from the literature. The AH and TC were modeled using both beam elements (wire model) and shell elements (ribbon model) in finite element analysis to predict their mechanical properties under tension. The asymptotic stiffening features of AH and TC under tensile loading were revealed and compared with a new analytical solution. The mechanical stiffnesses under other loading conditions, including compression, torsion and bending were also predicted numerically and correlated with the results of the existing MD simulations and tests. The mode shapes and natural frequencies of the spike were predicted using the built FE model. The frequencies were shown to be within the safe range of 1-20 MHz routinely used for medical imaging and diagnosis by means of ultrasound. These results provide a solid theoretical basis for using ultrasound to study damaging coronavirus through transient and resonant vibration at large deformations.

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.

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.

Computer-aided comprehensive explorations of RNA structural polymorphism through complementary simulation methods

While RNA folding was originally seen as a simple problem to solve, it has been shown that the promiscuous interactions of the nucleobases result in structural polymorphism, with several competing structures generally observed for non-coding RNA. This inherent complexity limits our understanding of these molecules from experiments alone, and computational methods are commonly used to study RNA. Here, we discuss three advanced sampling schemes, namely Hamiltonian-replica exchange molecular dynamics (MD), ratchet-and-pawl MD and discrete path sampling, as well as the HiRE-RNA coarse-graining scheme, and highlight how these approaches are complementary with reference to recent case studies. While all computational methods have their shortcomings, the plurality of simulation methods leads to a better understanding of experimental findings and can inform and guide experimental work on RNA polymorphism.

The short isoform of the host antiviral protein ZAP acts as an inhibitor of SARS-CoV-2 programmed ribosomal frameshifting

Programmed ribosomal frameshifting (PRF) is a fundamental gene expression event in many viruses, including SARS-CoV-2. It allows production of essential viral, structural and replicative enzymes that are encoded in an alternative reading frame. Despite the importance of PRF for the viral life cycle, it is still largely unknown how and to what extent cellular factors alter mechanical properties of frameshift elements and thereby impact virulence. This prompted us to comprehensively dissect the interplay between the SARS-CoV-2 frameshift element and the host proteome. We reveal that the short isoform of the zinc-finger antiviral protein (ZAP-S) is a direct regulator of PRF in SARS-CoV-2 infected cells. ZAP-S overexpression strongly impairs frameshifting and inhibits viral replication. Using in vitro ensemble and single-molecule techniques, we further demonstrate that ZAP-S directly interacts with the SARS-CoV-2 RNA and interferes with the folding of the frameshift RNA element. Together, these data identify ZAP-S as a host-encoded inhibitor of SARS-CoV-2 frameshifting and expand our understanding of RNA-based gene regulation.

Sterically confined rearrangements of SARS-CoV-2 Spike protein control cell invasion

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is highly contagious, and transmission involves a series of processes that may be targeted by vaccines and therapeutics. During transmission, host cell invasion is controlled by a large-scale (200–300 Å) conformational change of the Spike protein. This conformational rearrangement leads to membrane fusion, which creates transmembrane pores through which the viral genome is passed to the host. During Spike-protein-mediated fusion, the fusion peptides must be released from the core of the protein and associate with the host membrane. While infection relies on this transition between the prefusion and postfusion conformations, there has yet to be a biophysical characterization reported for this rearrangement. That is, structures are available for the endpoints, though the intermediate conformational processes have not been described. Interestingly, the Spike protein possesses many post-translational modifications, in the form of branched glycans that flank the surface of the assembly. With the current lack of data on the pre-to-post transition, the precise role of glycans during cell invasion has also remained unclear. To provide an initial mechanistic description of the pre-to-post rearrangement, an all-atom model with simplified energetics was used to perform thousands of simulations in which the protein transitions between the prefusion and postfusion conformations. These simulations indicate that the steric composition of the glycans can induce a pause during the Spike protein conformational change. We additionally show that this glycan-induced delay provides a critical opportunity for the fusion peptides to capture the host cell. In contrast, in the absence of glycans, the viral particle would likely fail to enter the host. This analysis reveals how the glycosylation state can regulate infectivity, while providing a much-needed structural framework for studying the dynamics of this pervasive pathogen.

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.

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

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

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

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

Biomolecular Modeling and Simulation: A Prospering Multidisciplinary Field

We reassess progress in the field of biomolecular modeling and simulation, following up on our perspective published in 2011. By reviewing metrics for the field's productivity and providing examples of success, we underscore the productive phase of the field, whose short-term expectations were overestimated and long-term effects underestimated. Such successes include prediction of structures and mechanisms; generation of new insights into biomolecular activity; and thriving collaborations between modeling and experimentation, including experiments driven by modeling. We also discuss the impact of field exercises and web games on the field's progress. Overall, we note tremendous success by the biomolecular modeling community in utilization of computer power; improvement in force fields; and development and application of new algorithms, notably machine learning and artificial intelligence. The combined advances are enhancing the accuracy andscope of modeling and simulation, establishing an exemplary discipline where experiment and theory or simulations are full partners.

Biomolecular modeling thrives in the age of technology

The biomolecular modeling field has flourished since its early days in the 1970s due to the rapid adaptation and tailoring of state-of-the-art technology. The resulting dramatic increase in size and timespan of biomolecular simulations has outpaced Moore's law. Here, we discuss the role of knowledge-based versus physics-based methods and hardware versus software advances in propelling the field forward. This rapid adaptation and outreach suggests a bright future for modeling, where theory, experimentation and simulation define three pillars needed to address future scientific and biomedical challenges.

A Fiedler Vector Scoring Approach for Novel RNA Motif Selection

Novel RNA motif design is of great practical importance for technology and medicine. Increasingly, computational design plays an important role in such efforts. Our coarse-grained RAG (RNA-As-Graphs) framework offers strategies for enumerating the universe of RNA 2D folds, selecting "RNA-like" candidates for design, and determining sequences that fold onto these candidates. In RAG, RNA secondary structures are represented as tree or dual graphs. Graphs with known RNA structures are called "existing", and the others are labeled "hypothetical". By using simplified features for RNA graphs, we have clustered the hypothetical graphs into "RNA-like" and "non-RNA-like" groups and proposed RNA-like graphs as candidates for design. Here, we propose a new way of designing graph features by using Fiedler vectors. The new features reflect graph shapes better, and they lead to a more clustered organization of existing graphs. We show significant increases in K-means clustering accuracy by using the new features (e.g., up to 95% and 98% accuracy for tree and dual graphs, respectively). In addition, we propose a scoring model for top graph candidate selection. This scoring model allows users to set a threshold for candidates, and it incorporates weighing of existing graphs based on their corresponding number of known RNAs. We include a list of top scored RNA-like candidates, which we hope will stimulate future novel RNA design.

A short review on nanotechnology interventions against COVID-19

The COVID-19 has affected all major aspects of the society in a global perspective. The role of nanotechnology is much sought after in fighting this pandemic. Advanced materials based on nanotechnology are the basis of several technologies starting from masks and personal protection equipment to specific diagnostic tools that could diminish the impact of COVID-19. Development of nanotechnology-based products is therefore an absolute necessity for fight against COVID-19. We examine the fundamental concepts related to virology, histopathologic findings and how nanotechnology can help in fighting the disease. In this review we discuss the state of the art and ongoing nanotechnology-based strategies like antiviral coatings, 3D printing and therapeutics to fight against this deadly disease. The importance of using nanoparticles in point of care tests and biosensors is also highlighted.

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.