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

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

Short 2'-O-methyl/LNA oligomers as highly-selective inhibitors of miRNA production in vitro and in vivo.

MicroRNAs (miRNAs) that share identical or near-identical sequences constitute miRNA families and are predicted to act redundantly. Yet recent evidence suggests that members of the same miRNA family with high sequence similarity might have different roles and that this functional divergence might be rooted in their precursors' sequence. Current knock-down strategies such as antisense oligonucleotides (ASOs) or miRNA sponges cannot distinguish between identical or near identical miRNAs originating from different precursors to allow exploring unique functions of these miRNAs. We here develop a novel strategy based on short 2'-OMe/LNA-modified oligonucleotides to selectively target specific precursor molecules and ablate the production of individual members of miRNA families in vitro and in vivo. Leveraging the highly conserved Xenopus miR-181a family as proof-of-concept, we demonstrate that 2'-OMe/LNA-ASOs targeting the apical region of pre-miRNAs achieve precursor-selective inhibition of mature miRNA-5p production. Furthermore, we extend the applicability of our approach to the human miR-16 family, illustrating its universality in targeting precursors generating identical miRNAs. Overall, our strategy enables efficient manipulation of miRNA expression, offering a powerful tool to dissect the functions of identical or highly similar miRNAs derived from different precursors within miRNA families.

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.

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 log10 units, whereas the use of plasmids carrying ribozymes or shRNAs alone resulted in lower inhibition.

A Structural Potential of Rare Trinucleotide Repeat Tracts in RNA.

Among types of trinucleotide repeats, there is some disproportion in the frequency of their occurrence in the human exome. This research presents new data describing the folding and thermodynamic stability of short, tandem RNA repeats of 23 types, focusing on the rare, yet poorly analyzed ones. UV-melting experiments included the presence of PEG or potassium and magnesium ions to determine their effect on the stability of RNA repeats structures. Rare repeats predominantly stayed single-stranded but had the potential for base pairing with other partially complementary repeat tracts. A coexistence of suitably complementary repeat types in a single RNA creates opportunities for interaction in the context of the secondary structure of RNA. We searched the human transcriptome for model RNAs in which different, particularly rare trinucleotide repeats coexist and selected the GABRA4 and CHIC1 RNAs to study intramolecular interactions between the repeat tracts that they contain. In vitro secondary structure probing results showed that the UAA and UUG repeat tracts, present in GABRA4 3' UTR, form a double helix, which separates one of its structural domains. For the RNA CHIC1 ORF fragment containing four short AGG repeat tracts and the CGU tract, we proved the formation of quadruplexes that blocked reverse transcription.

A Test and Refinement of Folding Free Energy Nearest Neighbor Parameters for RNA Including N6-Methyladenosine.

RNA folding free energy change parameters are widely used to predict RNA secondary structure and to design RNA sequences. These parameters include terms for the folding free energies of helices and loops. Although the full set of parameters has only been traditionally available for the four common bases and backbone, it is well known that covalent modifications of nucleotides are widespread in natural RNAs. Covalent modifications are also widely used in engineered sequences. We recently derived a full set of nearest neighbor terms for RNA that includes N6-methyladenosine (m6A). In this work, we test the model using 98 optical melting experiments, matching duplexes with or without N6-methylation of A. Most experiments place RRACH, the consensus site of N6-methylation, in a variety of contexts, including helices, bulge loops, internal loops, dangling ends, and terminal mismatches. For matched sets of experiments that include either A or m6A in the same context, we find that the parameters for m6A are as accurate as those for A. Across all experiments, the root mean squared deviation between estimated and experimental free energy changes is 0.67 kcal/mol. We used the new experimental data to refine the set of nearest neighbor parameter terms for m6A. These parameters enable prediction of RNA secondary structures including m6A, which can be used to model how N6-methylation of A affects RNA structure.

Secondary structure prediction for RNA sequences including N6-methyladenosine.

There is increasing interest in the roles of covalently modified nucleotides in RNA. There has been, however, an inability to account for modifications in secondary structure prediction because of a lack of software and thermodynamic parameters. We report the solution for these issues for N6-methyladenosine (m6A), allowing secondary structure prediction for an alphabet of A, C, G, U, and m6A. The RNAstructure software now works with user-defined nucleotide alphabets of any size. We also report a set of nearest neighbor parameters for helices and loops containing m6A, using experiments. Interestingly, N6-methylation decreases folding stability for adenosines in the middle of a helix, has little effect on folding stability for adenosines at the ends of helices, and increases folding stability for unpaired adenosines stacked on a helix. We demonstrate predictions for an N6-methylation-activated protein recognition site from MALAT1 and human transcriptome-wide effects of N6-methylation on the probability of adenosine being buried in a helix.

Secondary Structure of Influenza A Virus Genomic Segment 8 RNA Folded in a Cellular Environment.

Influenza A virus (IAV) is a member of the single-stranded RNA (ssRNA) family of viruses. The most recent global pandemic caused by the SARS-CoV-2 virus has shown the major threat that RNA viruses can pose to humanity. In comparison, influenza has an even higher pandemic potential as a result of its high rate of mutations within its relatively short (<13 kbp) genome, as well as its capability to undergo genetic reassortment. In light of this threat, and the fact that RNA structure is connected to a broad range of known biological functions, deeper investigation of viral RNA (vRNA) structures is of high interest. Here, for the first time, we propose a secondary structure for segment 8 vRNA (vRNA8) of A/California/04/2009 (H1N1) formed in the presence of cellular and viral components. This structure shows similarities with prior in vitro experiments. Additionally, we determined the location of several well-defined, conserved structural motifs of vRNA8 within IAV strains with possible functionality. These RNA motifs appear to fold independently of regional nucleoprotein (NP)-binding affinity, but a low or uneven distribution of NP in each motif region is noted. This research also highlights several accessible sites for oligonucleotide tools and small molecules in vRNA8 in a cellular environment that might be a target for influenza A virus inhibition on the RNA level.

Secondary Structure of Subgenomic RNA M of SARS-CoV-2.

SARS-CoV-2 belongs to the Coronavirinae family. Like other coronaviruses, SARS-CoV-2 is enveloped and possesses a positive-sense, single-stranded RNA genome of ~30 kb. Genomic RNA is used as the template for replication and transcription. During these processes, positive-sense genomic RNA (gRNA) and subgenomic RNAs (sgRNAs) are created. Several studies presented the importance of the genomic RNA secondary structure in SARS-CoV-2 replication. However, the structure of sgRNAs has remained largely unsolved so far. In this study, we probed the sgRNA M model of SARS-CoV-2 in vitro. The presented model molecule includes 5'UTR and a coding sequence of gene M. This is the first experimentally informed secondary structure model of sgRNA M, which presents features likely to be important in sgRNA M function. The knowledge of sgRNA M structure provides insights to better understand virus biology and could be used for designing new therapeutics.

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.

Universal and strain specific structure features of segment 8 genomic RNA of influenza A virus-application of 4-thiouridine photocrosslinking.

RNA structure in the influenza A virus (IAV) has been the focus of several studies that have shown connections between conserved secondary structure motifs and their biological function in the virus replication cycle. Questions have arisen on how to best recognize and understand the pandemic properties of IAV strains from an RNA perspective, but determination of the RNA secondary structure has been challenging. Herein, we used chemical mapping to determine the secondary structure of segment 8 viral RNA (vRNA) of the pandemic A/California/04/2009 (H1N1) strain of IAV. Additionally, this long, naturally occurring RNA served as a model to evaluate RNA mapping with 4-thiouridine (4sU) crosslinking. We explored 4-thiouridine as a probe of nucleotides in close proximity, through its incorporation into newly transcribed RNA and subsequent photoactivation. RNA secondary structural features both universal to type A strains and unique to the A/California/04/2009 (H1N1) strain were recognized. 4sU mapping confirmed and facilitated RNA structure prediction, according to several rules: 4sU photocross-linking forms efficiently in the double-stranded region of RNA with some flexibility, in the ends of helices, and across bulges and loops when their structural mobility is permitted. This method highlighted three-dimensional properties of segment 8 vRNA secondary structure motifs and allowed to propose several long-range three-dimensional interactions. 4sU mapping combined with chemical mapping and bioinformatic analysis could be used to enhance the RNA structure determination as well as recognition of target regions for antisense strategies or viral RNA detection.

Impact of a Single Nucleotide Change or Non-Nucleoside Modifications in G-Rich Region on the Quadruplex-Duplex Hybrid Formation.

In this paper, a method to discriminate between two target RNA sequences that differ by one nucleotide only is presented. The method relies on the formation of alternative structures, i.e., quadruplex-duplex hybrid (QDH) and duplex with dangling ends (Dss), after hybridization of DNA or RNA G-rich oligonucleotides with target sequences containing 5'-GGGCUGG-3' or 5'-GGGCGGG-3' fragments. Using biophysical methods, we studied the effect of oligonucleotide types (DNA, RNA), non-nucleotide modifications (aliphatic linkers or abasic), and covalently attached G4 ligand on the ability of G-rich oligonucleotides to assemble a G-quadruplex motif. We demonstrated that all examined non-nucleotide modifications could mimic the external loops in the G-quadruplex domain of QDH structures without affecting their stability. Additionally, some modifications, in particular the presence of two abasic residues in the G-rich oligonucleotide, can induce the formation of non-canonical QDH instead of the Dss structure upon hybridization to a target sequence containing the GGGCUGG motif. Our results offer new insight into the sequential requirements for the formation of G-quadruplexes and provide important data on the effects of non-nucleotide modifications on G-quadruplex formation.

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.

Prevalence of Anti-SARS-CoV-2 Antibodies in Poznan, Poland, after the First Wave of the COVID-19 Pandemic.

In comparison to other European countries, during the first months of the COVID-19 pandemic, Poland reported a relatively low number of confirmed cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections. To estimate the scale of the pandemic in Poland, a serosurvey of antibodies against SARS-CoV-2 was performed after the first wave of COVID-19 in Europe (March-May 2020). Within this study, we collected samples from 28 July to 24 September 2020 and, based on the ELISA results, we found that 1.67% (25/1500, 95% CI 1.13-2.45) of the Poznan (Poland) metropolitan area's population had antibodies against SARS-CoV-2 after the first wave of COVID-19. However, the presence of anti-SARS-CoV-2 IgG antibodies was confirmed with immunoblotting in 56% (14/25) samples, which finally resulted in a decrease in seroprevalence, i.e., 0.93% (14/1500, 95% CI 0.56-1.56). The positive anti-SARS-CoV-2 IgG results were associated with age, occupation involving constant contact with people, travelling abroad, non-compliance with epidemiological recommendations and direct contact with the novel coronavirus. Our findings confirm the low SARS-CoV-2 incidence in Poland and imply that the population had little herd immunity heading into the second and third wave of the pandemic, and therefore, that herd immunity contributed little to preventing the high numbers of SARS-CoV-2 infections and COVID-19-related deaths in Poland during these subsequent waves.

Conserved Structural Motifs of Two Distant IAV Subtypes in Genomic Segment 5 RNA.

The functionality of RNA is fully dependent on its structure. For the influenza A virus (IAV), there are confirmed structural motifs mediating processes which are important for the viral replication cycle, including genome assembly and viral packaging. Although the RNA of strains originating from distant IAV subtypes might fold differently, some structural motifs are conserved, and thus, are functionally important. Nowadays, NGS-based structure modeling is a source of new in vivo data helping to understand RNA biology. However, for accurate modeling of in vivo RNA structures, these high-throughput methods should be supported with other analyses facilitating data interpretation. In vitro RNA structural models complement such approaches and offer RNA structures based on experimental data obtained in a simplified environment, which are needed for proper optimization and analysis. Herein, we present the secondary structure of the influenza A virus segment 5 vRNA of A/California/04/2009 (H1N1) strain, based on experimental data from DMS chemical mapping and SHAPE using NMIA, supported by base-pairing probability calculations and bioinformatic analyses. A comparison of the available vRNA5 structures among distant IAV strains revealed that a number of motifs present in the A/California/04/2009 (H1N1) vRNA5 model are highly conserved despite sequence differences, located within previously identified packaging signals, and the formation of which in in virio conditions has been confirmed. These results support functional roles of the RNA secondary structure motifs, which may serve as candidates for universal RNA-targeting inhibitory methods.

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

Influenza is an important research subject around the world because of its threat to humanity. Influenza A virus (IAV) causes seasonal epidemics and sporadic, but dangerous pandemics. A rapid antigen changes and recombination of the viral RNA genome contribute to the reduced effectiveness of vaccination and anti-influenza drugs. Hence, there is a necessity to develop new antiviral drugs and strategies to limit the influenza spread. IAV is a single-stranded negative sense RNA virus with a genome (viral RNA-vRNA) consisting of eight segments. Segments within influenza virion are assembled into viral ribonucleoprotein (vRNP) complexes that are independent transcription-replication units. Each step in the influenza life cycle is regulated by the RNA and is dependent on its interplay and dynamics. Therefore, viral RNA can be a proper target to design novel therapeutics. Here, we briefly described examples of anti-influenza strategies based on the antisense oligonucleotide (ASO), small interfering RNA (siRNA), microRNA (miRNA) and catalytic nucleic acids. In particular we focused on the vRNA structure-function relationship as well as presented the advantages of using secondary structure information in predicting therapeutic targets and the potential future of this field.

The origin of the high stability of 3'-terminal uridine tetrads: contributions of hydrogen bonding, stacking interactions, and steric factors evaluated using modified oligonucleotide analogs.

RNA G-quadruplexes fold almost exclusively into parallel-stranded structures and thus display much less structural diversity than their DNA counterparts. However, also among RNA G-quadruplexes peculiar structural elements can be found which are capable of reshaping the physico-chemical properties of the folded structure. A striking example is provided by a uridine tetrad (U-tetrad) placed on the 3'-terminus of the tetramolecular G-quadruplex. In this context, the U-tetrad adopts a unique conformation involving chain reversal and is responsible for a tremendous stabilization of the G-quadruplex (ΔTm up to 30°C). In this report, we attempt to rationalize the origin of this stabilizing effect by concurrent structural, thermal stability, and molecular dynamics studies of a series of G-quadruplexes with subtle chemical modifications at their 3'-termini. Our results provide detailed insights into the energetics of the "reversed" U-tetrad motif and the requirements for its formation. They point to the importance of the 2'OH to phosphate hydrogen bond and preferential stacking interactions for the formation propensity and stability of the motif.

Conscious uncoupling of riboswitch functions.

Riboswitches alter gene expression in response to ligand binding, coupling sensing and regulatory functions to help bacteria respond to their environment. The structural determinants of ligand binding in the prequeuosine (7-aminomethyl-7-deazaguanine, preQ1) bacterial riboswitches have been studied, but the functional consequences of structural perturbations are less known. A new article combining biophysical and cell-based readouts of 15 mutants of the preQ1-II riboswitch from Lactobacillus rhamnosus demonstrates that ligand binding does not ensure successful gene regulation, providing new insights into these shapeshifting sequences.

RNA Secondary Structure Motifs of the Influenza A Virus as Targets for siRNA-Mediated RNA Interference.

The influenza A virus is a human pathogen that poses a serious public health threat due to rapid antigen changes and emergence of new, highly pathogenic strains with the potential to become easily transmitted in the human population. The viral genome is encoded by eight RNA segments, and all stages of the replication cycle are dependent on RNA. In this study, we designed small interfering RNA (siRNA) targeting influenza segment 5 nucleoprotein (NP) mRNA structural motifs that encode important functions. The new criterion for choosing the siRNA target was the prediction of accessible regions based on the secondary structure of segment 5 (+)RNA. This design led to siRNAs that significantly inhibit influenza virus type A replication in Madin-Darby canine kidney (MDCK) cells. Additionally, chemical modifications with the potential to improve siRNA properties were introduced and systematically validated in MDCK cells against the virus. A substantial and maximum inhibitory effect was achieved at concentrations as low as 8 nM. The inhibition of viral replication reached approximately 90% for the best siRNA variants. Additionally, selected siRNAs were compared with antisense oligonucleotides targeting the same regions; this revealed that effectiveness depends on both the target accessibility and oligonucleotide antiviral strategy. Our new approach of target-site preselection based on segment 5 (+)RNA secondary structure led to effective viral inhibition and a better understanding of the impact of RNA structural motifs on the influenza replication cycle.

Computational and NMR studies of RNA duplexes with an internal pseudouridine-adenosine base pair.

Pseudouridine (Ψ) is the most common chemical modification present in RNA. In general, Ψ increases the thermodynamic stability of RNA. However, the degree of stabilization depends on the sequence and structural context. To explain experimentally observed sequence dependence of the effect of Ψ on the thermodynamic stability of RNA duplexes, we investigated the structure, dynamics and hydration of RNA duplexes with an internal Ψ-A base pair in different nearest-neighbor sequence contexts. The structures of two RNA duplexes containing 5'-GΨC/3'-CAG and 5'-CΨG/3'-GAC motifs were determined using NMR spectroscopy. To gain insight into the effect of Ψ on duplex dynamics and hydration, we performed molecular dynamics (MD) simulations of RNA duplexes with 5'-GΨC/3'-CAG, 5'-CΨG/3'-GAC, 5'-AΨU/3'-UAA and 5'-UΨA/3'-AAU motifs and their unmodified counterparts. Our results showed a subtle impact from Ψ modification on the structure and dynamics of the RNA duplexes studied. The MD simulations confirmed the change in hydration pattern when U is replaced with Ψ. Quantum chemical calculations showed that the replacement of U with Ψ affected the intrinsic stacking energies at the base pair steps depending on the sequence context. The calculated intrinsic stacking energies help to explain the experimentally observed sequence dependent changes in the duplex stability from Ψ modification.