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

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.

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.

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.

A Short Chemically Modified dsRNA-Binding PNA (dbPNA) Inhibits Influenza Viral Replication by Targeting Viral RNA Panhandle Structure.

RNAs play critical roles in diverse catalytic and regulatory biological processes and are emerging as important disease biomarkers and therapeutic targets. Thus, developing chemical compounds for targeting any desired RNA structures has great potential in biomedical applications. The viral and cellular RNA sequence and structure databases lay the groundwork for developing RNA-binding chemical ligands through the recognition of both RNA sequence and RNA structure. Influenza A virion consists of eight segments of negative-strand viral RNA (vRNA), all of which contain a highly conserved panhandle duplex structure formed between the first 13 nucleotides at the 5' end and the last 12 nucleotides at the 3' end. Here, we report our binding and cell culture anti-influenza assays of a short 10-mer chemically modified double-stranded RNA (dsRNA)-binding peptide nucleic acid (PNA) designed to bind to the panhandle duplex structure through novel major-groove PNA·RNA2 triplex formation. We demonstrated that incorporation of chemically modified PNA residues thio-pseudoisocytosine (L) and guanidine-modified 5-methyl cytosine (Q) previously developed by us facilitates the sequence-specific recognition of Watson-Crick G-C and C-G pairs, respectively, at physiologically relevant conditions. Significantly, the chemically modified dsRNA-binding PNA (dbPNA) shows selective binding to the dsRNA region in panhandle structure over a single-stranded RNA (ssRNA) and a dsDNA containing the same sequence. The panhandle structure is not accessible to traditional antisense DNA or RNA with a similar length. Conjugation of the dbPNA with an aminosugar neamine enhances the cellular uptake. We observed that 2-5 µM dbPNA-neamine conjugate results in a significant reduction of viral replication. In addition, the 10-mer dbPNA inhibits innate immune receptor RIG-I binding to panhandle structure and thus RIG-I ATPase activity. These findings would provide the foundation for developing novel dbPNAs for the detection of influenza viral RNAs and therapeutics with optimal antiviral and immunomodulatory activities.

Secondary structure model of the naked segment 7 influenza A virus genomic RNA.

The influenza A virus (IAV) genome comprises eight negative-sense viral (v)RNA segments. The seventh segment of the genome encodes two essential viral proteins and is specifically packaged alongside the other seven vRNAs. To gain insights into the possible roles of RNA structure both within and without virions, a secondary structure model of a naked (protein-free) segment 7 vRNA (vRNA7) has been determined using chemical mapping and thermodynamic energy minimization. The proposed structure model was validated using microarray mapping, RNase H cleavage and comparative sequence analysis. Additionally, the detailed structures of three vRNA7 fragment constructs - comprising independently folded subdomains - were determined. Much of the proposed vRNA7 structure is preserved between IAV strains, suggesting their importance in the influenza replication cycle. Possible structure rearrangements, which allow or preclude long-range RNA interactions, are also proposed.

Microarrays for identifying binding sites and probing structure of RNAs.

Oligonucleotide microarrays are widely used in various biological studies. In this review, application of oligonucleotide microarrays for identifying binding sites and probing structure of RNAs is described. Deep sequencing allows fast determination of DNA and RNA sequence. High-throughput methods for determination of secondary structures of RNAs have also been developed. Those methods, however, do not reveal binding sites for oligonucleotides. In contrast, microarrays directly determine binding sites while also providing structural insights. Microarray mapping can be used over a wide range of experimental conditions, including temperature, pH, various cations at different concentrations and the presence of other molecules. Moreover, it is possible to make universal microarrays suitable for investigations of many different RNAs, and readout of results is rapid. Thus, microarrays are used to provide insight into oligonucleotide sequences potentially able to interfere with biological function. Better understanding of structure-function relationships of RNA can be facilitated by using microarrays to find RNA regions capable to bind oligonucleotides. That information is extremely important to design optimal sequences for antisense oligonucleotides and siRNA because both bind to single-stranded regions of target RNAs.

The contribution of pseudouridine to stabilities and structure of RNAs.

Thermodynamic data are reported revealing that pseudouridine (Ψ) can stabilize RNA duplexes when replacing U and forming Ψ-A, Ψ-G, Ψ-U and Ψ-C pairs. Stabilization is dependent on type of base pair, position of Ψ within the RNA duplex, and type and orientation of adjacent Watson-Crick pairs. NMR spectra demonstrate that for internal Ψ-A, Ψ-G and Ψ-U pairs, the N3 imino proton is hydrogen bonded to the opposite strand nucleotide and the N1 imino proton may also be hydrogen bonded. CD spectra show that general A-helix structure is preserved, but there is some shifting of peaks and changing of intensities. Ψ has two hydrogen donors (N1 and N3 imino protons) and two hydrogen bond acceptors because the glycosidic bond is C-C rather than C-N as in uridine. This greater structural potential may allow Ψ to behave as a kind of structurally driven universal base because it can enhance stability relative to U when paired with A, G, U or C inside a double helix. These structural and thermodynamic properties may contribute to the biological functions of Ψ.

Structural determinants for alternative splicing regulation of the MAPT pre-mRNA.

Alternative splicing at the MAPT gene exon 10 yields similar levels of the 3R and 4R tau protein isoforms. (1) The presence of mutations, particularly in exon 10 and intron 10-11, changes the quantity of tau isoforms. Domination each of the isoform yields tau protein aggregation and frontotemporal dementia and Parkinsonism linked to chromosome 17 (FTDP-17). Here, we report for the first time the secondary structure of the 194/195 nucleotide region for the wild type (WT) and 10 mutants of the MAPT gene pre-mRNA determined using both chemical and microarray mapping. Thermodynamic analyses indicate that single nucleotide mutations in the splicing regulatory element (SRE) that form a hairpin affect its stability by up to 4 and 7 kcal/mol. Moreover, binding the regulatory hairpin of small molecule ligands (neomycin, kanamycin, tobramycin and mitoxantrone) enhance its stability depending on the nature of the ligands and the RNA mutations. Experiments using the cos-7 cell line indicate that the presence of ligands and modified antisense oligonucleotides affect the quantity of 3R and 4R isoforms. This finding correlates with the thermodynamic stability of the regulatory hairpin. An alternative splicing regulation mechanism for exon 10 is postulated based on our experimental data and on published data.

Recognition of RNA duplexes by chemically modified triplex-forming oligonucleotides.

Triplex is emerging as an important RNA tertiary structure motif, in which consecutive non-canonical base pairs form between a duplex and a third strand. RNA duplex region is also often functionally important site for protein binding. Thus, triplex-forming oligonucleotides (TFOs) may be developed to regulate various biological functions involving RNA, such as viral ribosomal frameshifting and reverse transcription. How chemical modification in TFOs affects RNA triplex stability, however, is not well understood. Here, we incorporated locked nucleic acid, 2-thio U- and 2'-O methyl-modified residues in a series of all pyrimidine RNA TFOs, and we studied the binding to two RNA hairpin structures. The 12-base-triple major-groove pyrimidine-purine-pyrimidine triplex structures form between the duplex regions of RNA/DNA hairpins and the complementary RNA TFOs. Ultraviolet-absorbance-detected thermal melting studies reveal that the locked nucleic acid and 2-thio U modifications in TFOs strongly enhance triplex formation with both parental RNA and DNA duplex regions. In addition, we found that incorporation of 2'-O methyl-modified residues in a TFO destabilizes and stabilizes triplex formation with RNA and DNA duplex regions, respectively. The (de)stabilization of RNA triplex formation may be facilitated through modulation of van der Waals contact, base stacking, hydrogen bonding, backbone pre-organization, geometric compatibility and/or dehydration energy. Better understanding of the molecular determinants of RNA triplex structure stability lays the foundation for designing and discovering novel sequence-specific duplex-binding ligands as diagnostic and therapeutic agents targeting RNA.

Secondary structure of a conserved domain in an intron of influenza A M1 mRNA.

Influenza A virus utilizes RNA throughout infection. Little is known, however, about the roles of RNA structures. A previous bioinformatics survey predicted multiple regions of influenza A virus that are likely to generate evolutionarily conserved and stable RNA structures. One predicted conserved structure is in the pre-mRNA coding for essential proteins, M1 and M2. This structure starts 79 nucleotides downstream of the M2 mRNA 5' splice site. Here, a combination of biochemical structural mapping, mutagenesis, and NMR confirms the predicted three-way multibranch structure of this RNA. Imino proton NMR spectra reveal no change in secondary structure when 80 mM KCl is supplemented with 4 mM MgCl2. Optical melting curves in 1 M NaCl and in 100 mM KCl with 10 mM MgCl2 are very similar, with melting temperatures ∼14 °C higher than that for 100 mM KCl alone. These results provide a firm basis for designing experiments and potential therapeutics to test for function in cell culture.

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.

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

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.

Isoenergetic microarrays to study the structure and interactions of DsrA and OxyS RNAs in two- and three-component complexes.

Information on the secondary structure and interactions of RNA is important to understand the biological function of RNA as well as in applying RNA as a tool for therapeutic purposes. Recently, the isoenergetic microarray mapping method was developed to improve the prediction of RNA secondary structure. Herein, for the first time, isoenergetic microarrays were used to study the binding of RNA to protein or other RNAs as well as the interactions of two different RNAs and protein in a three-component complex. The RNAs used as models were the regulatory DsrA and OxyS RNAs from Escherichia coli, the fragments of their target mRNAs (fhlA and rpoS), and their complexes with Hfq protein. The collected results showed the advantages and some limitations of microarray mapping.

Structural diversity of triplet repeat RNAs.

Tandem repeats of various trinucleotide motifs are present in the human transcriptome, but the functions of these regular sequences, which likely depend on the structures they form, are still poorly understood. To gain new insight into the structural and functional properties of triplet repeats in RNA, we have performed a biochemical structural analysis of the complete set of triplet repeat transcripts, each composed of a single sequence repeated 17 times. We show that these transcripts fall into four structural classes. The repeated CAA, UUG, AAG, CUU, CCU, CCA, and UAA motifs did not form any higher order structure under any analyzed conditions. The CAU, CUA, UUA, AUG, and UAG repeats are ordered according to their increasing tendency to form semistable hairpins. The repeated CGA, CGU, and all CNG motifs form fairly stable hairpins, whereas AGG and UGG repeats fold into stable G-quadruplexes. The triplet repeats that formed the most stable structures were characterized further by biophysical methods. UV-monitored structure melting revealed that CGG and CCG repeats form, respectively, the most and least stable hairpins of all CNG repeats. Circular dichroism spectra showed that the AGG and UGG repeat quadruplexes are formed by parallel RNA strands. Furthermore, we demonstrated that the different susceptibility of various triplet repeat transcripts to serum nucleases can be explained by the sequence and structural features of the tested RNAs. The results of this study provide a comprehensive structural foundation for the functional analysis of triplet repeats in transcripts.