RNA Secondary Structure as a First Step for Rational Design of the Oligonucleotides towards Inhibition of Influenza A Virus Replication

In vivo secondary structural analysis of Influenza A virus genomic RNA

Influenza A virus (IAV) is a respiratory virus that causes epidemics and pandemics. Knowledge of IAV RNA secondary structure in vivo is crucial for a better understanding of virus biology. Moreover, it is a fundament for the development of new RNA-targeting antivirals. Chemical RNA mapping using selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) coupled with Mutational Profiling (MaP) allows for the thorough examination of secondary structures in low-abundance RNAs in their biological context. So far, the method has been used for analyzing the RNA secondary structures of several viruses including SARS-CoV-2 in virio and in cellulo. Here, we used SHAPE-MaP and dimethyl sulfate mutational profiling with sequencing (DMS-MaPseq) for genome-wide secondary structure analysis of viral RNA (vRNA) of the pandemic influenza A/California/04/2009 (H1N1) strain in both in virio and in cellulo environments. Experimental data allowed the prediction of the secondary structures of all eight vRNA segments in virio and, for the first time, the structures of vRNA5, 7, and 8 in cellulo. We conducted a comprehensive structural analysis of the proposed vRNA structures to reveal the motifs predicted with the highest accuracy. We also performed a base-pairs conservation analysis of the predicted vRNA structures and revealed many highly conserved vRNA motifs among the IAVs. The structural motifs presented herein are potential candidates for new IAV antiviral strategies.

Fighting the flu: a brief review on anti-influenza agents

The influenza virus causes one of the most prevalent and lethal infectious viral diseases of the respiratory system; the disease progression varies from acute self-limiting mild fever to disease chronicity and death. Although both the preventive and treatment measures have been vital in protecting humans against seasonal epidemics or sporadic pandemics, there are several challenges to curb the influenza virus such as limited or poor cross-protection against circulating virus strains, moderate protection in immune-compromised patients, and rapid emergence of resistance. Currently, there are four US-FDA-approved anti-influenza drugs to treat flu infection, viz. Rapivab, Relenza, Tamiflu, and Xofluza. These drugs are classified based on their mode of action against the viral replication cycle with the first three being Neuraminidase inhibitors, and the fourth one targeting the viral polymerase. The emergence of the drug-resistant strains of influenza, however, underscores the need for continuous innovation towards development and discovery of new anti-influenza agents with enhanced antiviral effects, greater safety, and improved tolerability. Here in this review, we highlighted commercially available antiviral agents besides those that are at different stages of development including under clinical trials, with a brief account of their antiviral mechanisms.

Structural and Functional RNA Motifs of SARS-CoV-2 and Influenza A Virus as a Target of Viral Inhibitors

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the COVID-19 pandemic, whereas the influenza A virus (IAV) causes seasonal epidemics and occasional pandemics. Both viruses lead to widespread infection and death. SARS-CoV-2 and the influenza virus are RNA viruses. The SARS-CoV-2 genome is an approximately 30 kb, positive sense, 5' capped single-stranded RNA molecule. The influenza A virus genome possesses eight single-stranded negative-sense segments. The RNA secondary structure in the untranslated and coding regions is crucial in the viral replication cycle. The secondary structure within the RNA of SARS-CoV-2 and the influenza virus has been intensively studied. Because the whole of the SARS-CoV-2 and influenza virus replication cycles are dependent on RNA with no DNA intermediate, the RNA is a natural and promising target for the development of inhibitors. There are a lot of RNA-targeting strategies for regulating pathogenic RNA, such as small interfering RNA for RNA interference, antisense oligonucleotides, catalytic nucleic acids, and small molecules. In this review, we summarized the knowledge about the inhibition of SARS-CoV-2 and influenza A virus propagation by targeting their RNA secondary structure.

ScanFold 2.0: a rapid approach for identifying potential structured RNA targets in genomes and transcriptomes

A major limiting factor in target discovery for both basic research and therapeutic intervention is the identification of structural and/or functional RNA elements in genomes and transcriptomes. This was the impetus for the original ScanFold algorithm, which provides maps of local RNA structural stability, evidence of sequence-ordered (potentially evolved) structure, and unique model structures comprised of recurring base pairs with the greatest structural bias. A key step in quantifying this propensity for ordered structure is the prediction of secondary structural stability for randomized sequences which, in the original implementation of ScanFold, is explicitly evaluated. This slow process has limited the rapid identification of ordered structures in large genomes/transcriptomes, which we seek to overcome in this current work introducing ScanFold 2.0. In this revised version of ScanFold, we no longer explicitly evaluate randomized sequence folding energy, but rather estimate it using a machine learning approach. For high randomization numbers, this can increase prediction speeds over 100-fold compared to ScanFold 1.0, allowing for the analysis of large sequences, as well as the use of additional folding algorithms that may be computationally expensive. In the testing of ScanFold 2.0, we re-evaluate the Zika, HIV, and SARS-CoV-2 genomes and compare both the consistency of results and the time of each run to ScanFold 1.0. We also re-evaluate the SARS-CoV-2 genome to assess the quality of ScanFold 2.0 predictions vs several biochemical structure probing datasets and compare the results to those of the original ScanFold program.

Antiviral Activity of Oligonucleotides Targeting the SARS-CoV-2 Genomic RNA Stem-Loop Sequences within the 3'-End of the ORF1b

Increased evidence shows vaccines against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) exhibited no long-term efficacy and limited worldwide availability, while existing antivirals and treatment options have only limited efficacy. In this study, the main objective was the development of antiviral strategies using nucleic acid-based molecules. To this purpose, partially overlapped 6-19-mer phosphorothioate deoxyoligonucleotides (S-ONs) designed on the SARS-CoV-2 genomic RNA stem-loop packaging sequences within the 3' end of the ORF1b were synthetized using the direct and complementary sequence. Among the S-ONs tested, several oligonucleotides exhibited a fifty percent inhibitory concentration antiviral activity ranging from 0.27 to 34 μM, in the absence of cytotoxicity. The S-ON with a scrambled sequence used in the same conditions was not active. Moreover, selected 10-mer S-ONs were tested using different infectious doses and against different SARS-CoV-2 variants, showing comparable antiviral activity that was abrogated when the central sequence was mutated. Experiments to evaluate the intracellular functional target localization of the S-ON inhibitory activity were also performed. Collectively the data indicate that the SARS-CoV-2 packaging region in the 3' end of the ORF1b may be a promising target candidate for further investigation to develop innovative nucleic-acid-based antiviral therapy.

Structural variants and modifications of hammerhead ribozymes targeting influenza A virus conserved structural motifs

The naturally occurring structure and biological functions of RNA are correlated, which includes hammerhead ribozymes. We proposed new variants of hammerhead ribozymes targeting conserved structural motifs of segment 5 of influenza A virus (IAV) (+)RNA. The variants carry structural and chemical modifications aiming to improve the RNA cleavage activity of ribozymes. We introduced an additional hairpin motif and attempted to select ribozyme-target pairs with sequence features that enable the potential formation of the trans-Hoogsteen interactions that are present in full-length, highly active hammerhead ribozymes. We placed structurally defined guanosine analogs into the ribozyme catalytic core. Herein, the significantly improved synthesis of 2′-deoxy-2′-fluoroarabinoguanosine derivatives is described. The most potent hammerhead ribozymes were applied to chimeric short hairpin RNA (shRNA)-ribozyme plasmid constructs to improve the antiviral activity of the two components. The modified hammerhead ribozymes showed moderate cleavage activity. Treatment of IAV-infected Madin-Darby canine kidney (MDCK) cells with the plasmid constructs resulted in significant inhibition of virus replication. Real-time PCR analysis revealed a significant (80%–88%) reduction in viral RNA when plasmids carriers were used. A focus formation assay (FFA) for chimeric plasmids showed inhibition of virus replication by 1.6–1.7 log₁₀ units, whereas the use of plasmids carrying ribozymes or shRNAs alone resulted in lower inhibition.

Expansion of the RNAStructuromeDB to include secondary structural data spanning the human protein-coding transcriptome

RNA plays vital functional roles in almost every component of biology, and these functional roles are often influenced by its folding into secondary and tertiary structures. An important role of RNA secondary structure is in maintaining proper gene regulation; therefore, making accurate predictions of the structures involved in these processes is important. In this study, we have expanded on our previous work that led to the creation of the RNAStructuromeDB. Unlike this previous study that analyzed the human genome at low resolution, we have now scanned the protein-coding human transcriptome at high (single nt) resolution. This provides more robust structure predictions for over 100,000 isoforms of known protein-coding genes. Notably, we also utilize the motif identification tool, ScanFold, to model structures with high propensity for ordered/evolved stability. All data have been uploaded to the RNAStructuromeDB, allowing for easy searching of transcripts, visualization of data tracks (via the Integrative Genomics Viewer or IGV), and download of ScanFold data—including unique highly-ordered motifs. Herein, we provide an example analysis of MAT2A to demonstrate the utility of ScanFold at finding known and novel secondary structures, highlighting regions of potential functionality, and guiding generation of functional hypotheses through use of the data.

Nearest neighbor rules for RNA helix folding thermodynamics: improved end effects

Nearest neighbor parameters for estimating the folding stability of RNA secondary structures are in widespread use. For helices, current parameters penalize terminal AU base pairs relative to terminal GC base pairs. We curated an expanded database of helix stabilities determined by optical melting experiments. Analysis of the updated database shows that terminal penalties depend on the sequence identity of the adjacent penultimate base pair. New nearest neighbor parameters that include this additional sequence dependence accurately predict the measured values of 271 helices in an updated database with a correlation coefficient of 0.982. This refined understanding of helix ends facilitates fitting terms for base pair stacks with GU pairs. Prior parameter sets treated 5′GGUC3′ paired to 3′CUGG5′ separately from other 5′GU3′/3′UG5′ stacks. The improved understanding of helix end stability, however, makes the separate treatment unnecessary. Introduction of the additional terms was tested with three optical melting experiments. The average absolute difference between measured and predicted free energy changes at 37°C for these three duplexes containing terminal adjacent AU and GU pairs improved from 1.38 to 0.27 kcal/mol. This confirms the need for the additional sequence dependence in the model.

Nuclear magnetic resonance reveals a two hairpin equilibrium near the 3'-splice site of influenza A segment 7 mRNA that can be shifted by oligonucleotides

Influenza A kills hundreds of thousands of people globally every year and has the potential to generate more severe pandemics. Influenza A's RNA genome and transcriptome provide many potential therapeutic targets. Here, nuclear magnetic resonance (NMR) experiments suggest that one such target could be a hairpin loop of 8 nucleotides in a pseudoknot that sequesters a 3' splice site in canonical pairs until a conformational change releases it into a dynamic 2 × 2-nt internal loop. NMR experiments reveal that the hairpin loop is dynamic and able to bind oligonucleotides as short as pentamers. A 3D NMR structure of the complex contains 4 and likely 5 bp between pentamer and loop. Moreover, a hairpin sequence was discovered that mimics the equilibrium of the influenza hairpin between its structure in the pseudoknot and upon release of the splice site. Oligonucleotide binding shifts the equilibrium completely to the hairpin secondary structure required for pseudoknot folding. The results suggest this hairpin can be used to screen for compounds that stabilize the pseudoknot and potentially reduce splicing.

In Vitro Methods to Decipher the Structure of Viral RNA Genomes

RNA viruses encode essential information in their genomes as conserved structural elements that are involved in efficient viral protein synthesis, replication, and encapsidation. These elements can also establish complex networks of RNA-RNA interactions, the so-called RNA interactome, to shape the viral genome and control different events during intracellular infection. In recent years, targeting these conserved structural elements has become a promising strategy for the development of new antiviral tools due to their sequence and structural conservation. In this context, RNA-based specific therapeutic strategies, such as the use of siRNAs have been extensively pursued to target the genome of different viruses. Importantly, siRNA-mediated targeting is not a straightforward approach and its efficiency is highly dependent on the structure of the target region. Therefore, the knowledge of the viral structure is critical for the identification of potentially good target sites. Here, we describe detailed protocols used in our laboratory for the in vitro study of the structure of viral RNA genomes. These protocols include DMS (dimethylsulfate) probing, SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) analysis, and HMX (2'-hydroxyl molecular interference). These methodologies involve the use of high-throughput analysis techniques that provide extensive information about the 3D folding of the RNA under study and the structural tuning derived from the interactome activity. They are therefore a good tool for the development of new RNA-based antiviral compounds.

Effective Inhibition of Newly Emerged A/H7N9 Virus with Oligonucleotides Targeted to Conserved Regions of the Virus Genome

Newly emerged highly pathogenic A/H7N9 viruses with pandemic potential are effectively transmitted from birds to humans and require the development of novel antiviral drugs. For the first time, we studied the in vitro and in vivo antiviral activity against A/H7N9 of oligodeoxyribonucleotides (ODNs), which were delivered into the cells in the proposed TiO₂-based nanocomposites (TiO₂∼ODN). The highest inhibition of A/H7N9 in vitro (∼400-fold) and efficient, sequence-specific, and dose-dependent protection (up to 100%) of A/H7N9-infected mice was revealed when ODN was targeted to the conserved terminal 3'-noncoding region of viral (-)RNA. After the treatment with ODN, the virus titer values in the lungs of mice decreased by several orders of magnitude. The TiO₂∼ODN nanocomposite did not show toxicity in mice under the treatment conditions. The proposed approach for effective inhibition of the A/H7N9 can be tested against other viruses, for example, new emerging influenza viruses and coronaviruses with pandemic potential.

Identification and Structural Aspects of G-Quadruplex-Forming Sequences from the Influenza A Virus Genome

Influenza A virus (IAV) causes seasonal epidemics and sporadic pandemics, therefore is an important research subject for scientists around the world. Despite the high variability of its genome, the structure of viral RNA (vRNA) possesses features that remain constant between strains and are biologically important for virus replication. Therefore, conserved structural motifs of vRNA can represent a novel therapeutic target. Here, we focused on the presence of G-rich sequences within the influenza A/California/07/2009(H1N1) genome and their ability to form RNA G-quadruplex structures (G4s). We identified 12 potential quadruplex-forming sequences (PQS) and determined their conservation among the IAV strains using bioinformatics tools. Then we examined the propensity of PQS to fold into G4s by various biophysical methods. Our results revealed that six PQS oligomers could form RNA G-quadruplexes. However, three of them were confirmed to adopt G4 structures by all utilized methods. Moreover, we showed that these PQS motifs are present within segments encoding polymerase complex proteins indicating their possible role in the virus biology.

In Vitro Inhibitory Analysis of Rationally Designed siRNAs against MERS-CoV Replication in Huh7 Cells

MERS-CoV was identified for the first time in Jeddah, Saudi Arabia in 2012 in a hospitalized patient. This virus subsequently spread to 27 countries with a total of 939 deaths and 2586 confirmed cases and now has become a serious concern globally. Camels are well known for the transmission of the virus to the human population. In this report, we have discussed the prediction, designing, and evaluation of potential siRNA targeting the ORF1ab gene for the inhibition of MERS-CoV replication. The online software, siDirect 2.0 was used to predict and design the siRNAs, their secondary structure and their target accessibility. ORF1ab gene folding was performed by RNAxs and RNAfold software. A total of twenty-one siRNAs were selected from 462 siRNAs according to their scoring and specificity. siRNAs were evaluated in vitro for their cytotoxicity and antiviral efficacy in Huh7 cell line. No significant cytotoxicity was observed for all siRNAs in Huh7 cells. The in vitro study showed the inhibition of viral replication by three siRNAs. The data generated in this study provide preliminary and encouraging information to evaluate the siRNAs separately as well as in combination against MERS-CoV replication in other cell lines. The prediction of siRNAs using online software resulted in the filtration and selection of potential siRNAs with high accuracy and strength. This computational approach resulted in three effective siRNAs that can be taken further to in vivo animal studies and can be used to develop safe and effective antiviral therapies for other prevalent disease-causing viruses.

Developing therapeutic approaches for twenty-first-century emerging infectious viral diseases

The twenty-first century has already recorded more than ten major epidemic or pandemic virus emergence events, including the ongoing and devastating coronavirus disease 2019 (COVID-19) pandemic. As viral disease emergence is expected to accelerate, these data dictate a need for proactive approaches to develop broadly active family-specific and cross-family therapeutics for use in future disease outbreaks. Emphasis should focus not only on the development of broad-spectrum small-molecule and antibody direct-acting antivirals, but also on host-factor therapeutics, including repurposing previously approved or in-pipeline drugs. Another new class of therapeutics with great antiviral therapeutic potential is RNA-based therapeutics. Rather than only focusing on known risks, dedicated efforts must be made toward pre-emptive research focused on outbreak-prone virus families, ultimately offering a strategy to shorten the gap between outbreak and response. Emphasis should also focus on orally available drugs for outpatient use, if possible, and on identifying combination therapies that combat viral and immune-mediated pathologies, extend the effectiveness of therapeutic windows and reduce drug resistance. While such an undertaking will require new vision, dedicated funding and private, federal and academic partnerships, this approach offers hope that global populations need never experience future pandemics such as COVID-19.

Targeting the RdRp of Emerging RNA Viruses: The Structure-Based Drug Design Challenge

The RNA-dependent RNA polymerase (RdRp) is an essential enzyme for the viral replication process, catalyzing the viral RNA synthesis using a metal ion-dependent mechanism. In recent years, RdRp has emerged as an optimal target for the development of antiviral drugs, as demonstrated by recent approvals of sofosbuvir and remdesivir against Hepatitis C virus (HCV) and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), respectively. In this work, we overview the main sequence and structural features of the RdRp of emerging RNA viruses such as Coronaviruses, Flaviviruses, and HCV, as well as inhibition strategies implemented so far. While analyzing the structural information available on the RdRp of emerging RNA viruses, we provide examples of success stories such as for HCV and SARS-CoV-2. In contrast, Flaviviruses' story has raised attention about how the lack of structural details on catalytically-competent or ligand-bound RdRp strongly hampers the application of structure-based drug design, either in repurposing and conventional approaches.