A Conserved Secondary Structural Element in the Coding Region of the Influenza A Virus Nucleoprotein (NP) mRNA Is Important for the Regulation of Viral Proliferation (original) (raw)

Open Access

Peer-reviewed

Research Article

A Conserved Secondary Structural Element in the Coding Region of the Influenza A Virus Nucleoprotein (NP) mRNA Is Important for the Regulation of Viral Proliferation

PLOS

x

Figures

Abstract

Influenza A virus is a threat to humans due to seasonal epidemics and infrequent, but dangerous, pandemics that lead to widespread infection and death. Eight segments of RNA constitute the genome of this virus and they encode greater than eight proteins via alternative splicing of coding (+)RNAs generated from the genomic (-)RNA template strand. RNA is essential in its life cycle. A bioinformatics analysis of segment 5, which encodes nucleoprotein, revealed a conserved structural motif in the (+)RNA. The secondary structure proposed by energy minimization and comparative analysis agrees with structure predicted based on experimental data using a 121 nucleotide in vitro RNA construct comprising an influenza A virus consensus sequence and also an entire segment 5 (+)RNA (strain A/VietNam/1203/2004 (H5N1)). The conserved motif consists of three hairpins with one being especially thermodynamically stable. The biological importance of this conserved secondary structure is supported in experiments using antisense oligonucleotides in cell line, which found that disruption of this motif led to inhibition of viral fitness. These results suggest that this conserved motif in the segment 5 (+)RNA might be a candidate for oligonucleotide-based antiviral therapy.

Citation: Soszynska-Jozwiak M, Michalak P, Moss WN, Kierzek R, Kierzek E (2015) A Conserved Secondary Structural Element in the Coding Region of the Influenza A Virus Nucleoprotein (NP) mRNA Is Important for the Regulation of Viral Proliferation. PLoS ONE 10(10): e0141132. https://doi.org/10.1371/journal.pone.0141132

Editor: Eric Jan, University of British Columbia, CANADA

Received: August 11, 2015; Accepted: October 5, 2015; Published: October 21, 2015

Copyright: © 2015 Soszynska-Jozwiak et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Data Availability: All relevant data are within the paper.

Funding: This work was supported by National Science Center grants: UMO-2013/08/M/NZ1/01062 to EK and UMO-2013/08/A/ST5/00295 to RK, American Cancer Society grant 126117-PF-14-139-01-MPC and National Institutes of Health grant K99GM112877-01 to WNM, and by NIH grant R03TW008739 to EK. This publication was also supported by the Polish Ministry of Science and Higher Education, under the KNOW program.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Influenza A virus is a grave threat to human health. In 1918 the Spanish flu (H1N1strain) caused the deaths of over 50 million people [1, 2]. In 1997, migration of the H5N1 strain from bird to human was verified by molecular analysis [3]. In 2007, a new pandemic strain of influenza virus H1N1 was observed. This is particularly virulent and can rapidly disseminate itself through humans [4]. New pandemic strains of influenza virus are likely to arise in the 21st century, which makes the investigation of novel therapeutic targets in influenza especially relevant.

Influenza A virus belongs to the family Orthomyxoviridae and possesses a segmented, negative-sense RNA genomic viral (v)RNA. RNA is used throughout infection and plays roles in every process of viral life cycle. Replication begins with vRNA acting as a template to produce two plus-sense (+)RNAs: a complementary (c)RNA intermediate, which in turn becomes the template for producing more vRNA strands, as well as protein coding mRNA [5, 6]. A number of publications have suggested that RNA secondary structure plays important roles in influenza A infection [7–15]. A bioinformatics analysis based on identifying regions of unusual thermodynamic stability and structural conservation revealed that the (+) sense influenza RNA contains at least twelve structured motifs with likely function [12, 15]. Several of these regions overlapped areas of suppressed synonymous codon usage [12, 16], which suggests that RNA structure is exerting an evolutionary constraint on influenza A virus codon evolution.

The best-studied structured regions occur in the mRNAs of segments 7 and 8, where structures are modeled at or near splice sites. Indeed, predicted structural regions appear at or near splice sites in influenza B and C, suggesting common roles for structure in the regulation of influenza alternative splicing [16]. In influenza A, for example, a 63 nucleotide (nt) conserved region was identified in the segment 7 mRNA. This region can fold in two conformations: a hairpin and pseudoknot [11]. Changing between these two conformations places splicing regulatory elements into varying structural contexts, which has likely implications on segment 7 splicing. Another conserved structure occurs in the intron of segment 8 mRNA [14]. The function of this domain is unclear, but its proximity to the 5’ splice site makes it a possible intronic splicing enhancer/inhibitor. Also the defined structure could be a tag for recognition of proteins to distinguish unspliced NS1 mRNA from spliced NEP mRNA. Another conserved motif, in influenza mRNA of segment 7 was determined biochemically [13]. A multibranch loop structure is proposed to modulate alternative splicing of segment 7; as deduced by comparing previous point mutations studies [17] to the model [13]. Understanding the roles of influenza RNA thus plays an important role in gaining mechanistic insights into influenza virology and, significantly, in designing new drugs that can target viral RNA/RNA structure.

A previous bionformatic analysis of six (+)RNA sequences from segment 5 predicted a high probability of structure in the region spanning nts 1031–1250 [12]. This regions has unusually stable predicted thermodynamic stability and the model base pairs are conserved between homologous sequences. Additionally, the structure in this region appears to be influencing amino acid codon evolution, where synonymous codon use is highly suppressed (vs. other regions of influenza); presumably due to the need of maintaining structure in base-paired third codon (wobble) sites. In this study, we focused on a structured region of influenza A segment 5. This segment encodes nucleoprotein (NP), which is a structural protein and also takes part in regulation of transcription and replication. NP with viral polymerase and vRNAs create eight ribonucleoprotein complexes (RNP) which determine virus particle structure. Here, a new predicted motif in (+) RNA of segment 5 (nts 1051 to 1171) was confirmed in vitro and in cell line experiments. Chemical and enzymatic probing, Pb2+ cleavage and isoenergetic microarray mapping were used to probe the secondary structure a truncated in vitro construct comprising the consensus sequence of this NP mRNA motif. This isolated structure was also confirmed by probing the entire segment 5 (+)RNA of strain A/VietNam/1203/2004 (H5N1). Antisense oligonucleotides targeting a hairpin in this conserved motif are able inhibit virus proliferation in cell line, suggesting functional significance for this structure.

Materials and Methods

Oligonucleotides synthesis

Primers for PCR and reverse transcription and microarrays probes and antisense oligonucleotides were synthesized on MerMade12 (BioAutomation) synthesizer, using β-cyanoethyl phosphoramidite chemistry on solid support. Oligonucleotides were deprotected and purified according to published protocol [18, 19]. Concentrations of all oligonucleotides were measured using a Spectrophotometer UV (Picodrop-Syngen).

Experimental Constructs

DNA templates for synthesis of M121 was obtained from PCRs reactions using appropriate primers (Table 1). Two partially complementary DNA oligonucleotides S5-1 and S5-2 and were used in first PCR reaction. The product from the PCR was temple for second PCR reaction with primers contains EcoRI and PstI restriction sites: M-F and M-R. The PCR final product was ligated into a digested pUC19 vector. The resulting plasmid was cloned into DH5-alpha cells. Plasmid was isolated from cloned cells using a Qiagen mini-prep kit and the sequence was confirmed by sequencing.

M121 RNA and (+)RNA5 synthesis

PCR reaction with specific primers T1-F and T2-R on above digested plasmid was perform to produce templates for transcription. In vitro transcription reactions was performed using an Ampliscribe T7 Flash kit (Epicenter). Product was purified by denaturing 8% PAGE gel and eluted. M121 RNA for folding experiment, enzymatic mapping and lead ion cleavage was 5’ end labeled with γ-32P ATP (Hartmann Analityc GmbH), then re-purified by denaturing PAGE.

DNA template for in vitro transcriptions of (+)RNA5 (1565 nt) was obtained by PCR from vector pPol1 using primers For_c5 and Rev_c5 (Table 1). The pPol1 vector containing DNA of segment 5 influenza strain A/Viet Nam/1203/2004 (H5N1) was received from Prof. Baek Kim, University of Rochester. (+)RNA5 was purified using RNeasy MiniElute Cleanup kit (Qiagen) and quality of RNA was checked on agarose gel.

M121 RNA and (+)RNA5 Folding

For each sample, about 20,000 cpm of 5' labeled M121 RNA in folding buffer (100 mM Tris-HCl pH 7, 100 mM KCl and different concentration of MgCl2) was heated in 65°C for 5 min and slowly cooled to room temperature for 50 min. To study the multivalent cation dependent folding, MgCl2 was used to get a range of final concentration from 2.5 mM to 25 mM. Also Co3+ dependence folding was studied using Co(NH3)6Cl3 added to span 0.002 to 10 mM. Folding was analyzed by native gel electrophoresis using 6% polyacrylamide gel with 1x THEM (34 mM Tris Base, 57 mM HEPES, 0.1 mM EDTA, 10 mM MgCl2) buffer running at 4°C with low voltage. In all conditions one band was achieved. For next experiments the folding buffer contained 100 mM KCl and 5 mM MgCl2 was chosen. Entire (+)RNA5 was folded in the same condition and obtaining of one structure was checked on agarose gel in native conditions (running at 4°C).

Chemical and Enzymatic Mapping and lead ion cleavage

M121 and (+)RNA5 for SHAPE (Selective 2′-Hydroxyl Acylation analyzed by Primer Extension) and chemical mapping experiments were folded in buffer 100 mM KCl, 5 mM MgCl2, 50 mM HEPES pH 7.0. Enzymatic mapping and lead cleavage used 10 mM Tris-HCl pH 7.0 instead of 50 mM HEPES pH 7.0 and were done only for M121 RNA. Enzymatic and chemical mapping was carried out at room temperature.

Conditions for RNase S1, V1, T1 cleavages and RNase T1 ladder were adapted from manufacturer's protocol (Ambion). Optimal enzymatic concentrations were determined with enzymes titration. Lead ion cleavage reactions were carried out by incubating 5'-end labeled M121 RNA with 1 mM Pb(OAc)2. Aliquots with RNA (20,000 cpm per aliquot) were removed at 0, 1, 5, 15, 30, 60 min. The enzymatic or hydrolysis reactions were stopped by placing the aliquots in the gel loading buffer and freezing at -80°C until they were fractionated on a denaturing, 12% polyacrylamide gel. All gel were dried, exposed to phosphorimager screen, imaged and analyzed in MultiGauge program.

M121 RNA and (+)RNA5 were folded in folding buffer (100 mM KCl, 5 mM MgCl2, 50 mM HEPES, pH 7.0) as described above. Next, RNA was modified with optimized concentrations of DMS (dimethyl sulfate), CMCT (1-cyclohexyl-(2-morpholinoethyl) carbodiimide metho-p-toluene sulfonate), kethoxal and NMIA (N-methylisatoic anhydride) using slightly modified published protocols [20, 21]. For single reaction 1 pmol of RNA was used. Chemical mapping was carried out at room temperature with DMS, CMCT, kethoxal or NMIA for 15, 20, 30 or 40 min, respectively. Control reaction was done in the same conditions but without reagents. Modifications were read out by primer extension using primers listed in Table 2, with reverse transcriptase SuperScript III and Invitrogen protocols. RNA fragments and ddNTP ladders were separated by capillary electrophoresis.

Data were analyzed by PeakScanner. Quantitative NMIA, DMS, CMCT and kethoxal reactivaties for individual datasets were normalized to a scale in which 0 indicates an unreactive site and the average intensity at highly reactive sites is set to 1.0. The normalization factor for each datasets was determined by first excluding the most reactive 2% of peak intensities and then calculating the average for the next 8% of peak intensities. All reactivities were then divided by this average.

Hybridization to Isoenergetic Microarrays

Isoenergetic microarrays were prepared and microarrays mapping was performed similarly to published procedures [22–24]. About 200,000 cpm of labeled RNA was folded as described above before hybridization to isoenergetic microarrays. Probes of universal microarray was spotted in triplicate. Spotting buffer, monomer U, and pentamer UUUUU (2’-O-methyl-RNA), which should show no binding to RNA, were also printed on the microarray as internal negative controls. Hybridizations were carried out in folding buffer for 18 h at 4, 23 and 37°C. Microarrays were washed for 1 min at the same temperature as hybridization occurred and then dried by centrifugation. Hybridization was visualized by exposure to a phosphorimager screen and quantitative analysis was performed with ArrayGaugeV2.1. Binding was considered strong, medium and weak, when the integrated spot intensity was ≥ 1/3, ≥ 1/9, or ≥ 1/27 of the strongest integrated spot intensity, respectively. Alternative binding sites were predicted using RNA/RNA thermodynamics.

Sequence Alignment and Analysis

All available influenza A sequences for segment 5 were acquired from the NCBI influenza virus resource [25]. Sequences were filtered to remove identical sequences and leave a set of unique RNAs. From this set of 2515 sequences subsets of 100 were randomly selected for analysis with the program RNAz [26, 27] which scans for RNA structure using overlapping prediction windows of 120 nt in 10 nt steps. Overlapping windows with high propensity of forming functional RNA structure were concatenated and that region of the alignment was extracted for initial in silico modeling with the program RNAalifold [28] that predicts consensus folds for fixed alignments. Model structures were then evaluated for conservation with respect to the full alignment of all unique influenza A sequences.

Cells experiments

Cell line experiments were performed using single-cycle infectious influenza A virus (sciIAV) of strain A/California/04_NYICE_E3/2009 (H1N1) in MDCK-HA (Madin–Darby canine kidney epithelial cells-Hemagglutinin) cells, both were gift from Prof. Luis Martinez-Sobrido, University of Rochester [29, 30]. Cells were grown at 37°C with 5% CO2 in Dulbecco's modified Eagle's medium (DMEM) supplemented with 10% fetal bovine serum (FBS) and penicillin-streptomycin-glutamine. Every third passage hygromycin B was added to constitutive express hemagglutinin. 96-Well plates were seeded with MDCK-HA cells (10,000 cells per well) 24 h before transfection. Lipofectamine 2000 was used to transfect 4 uM antisense oligonucleotides and 12 h after transfection, medium was changed and cells were incubated for 6 h at 37°C. Next, infection with low MOI (0.001) was performed. Post infection takes 40 h at 33°C (two viral life cycle) and all supernatants were collected for analysis. Next, serial dilutions of cell culture supernatant were used for infections of MDCK-HA monolayers in 96-well plates. Post infection time was 20 h at 33°C (one viral life cycle); this unable counting GFP foci and calculation of FFU/ml. Controls were infected cells in medium, infected cells treated with lipofectamine 2000 without any antisense oligonucleotide and infected cells transfected with 2’-O-methyl-RNA oligonucleotide AGACCUCUAUAGCAGCU (N) as negative controls [31]. Ribavirin at concentration 40 uM added in postinfection media served as positive control.

Results and Discussion

Identification of a conserved structural motif in segment 5 (+)RNA

Here, we reanalyzed the segment 5 RNA using a larger set of 100 sequences with the ncRNA discovery program RNAz [26, 32]. This new analysis revealed one predicted region with a high propensity of forming structured RNA: support vector machine class probability (p-class) 0.90 and z-score of -2.49. The high p-class indicates that this region has similar calculated metrics (e.g. stability and structure conservation) as known structured ncRNAs [26, 32], and the z-score indicates that, compared to random sequence, this region is ~2.5 standard deviations more thermodynamically stable. Furthermore, compared to previous work on segment 5, the identified structured region is narrowed to nts 1051–1171, whereas the previously identified region spanned nt 1031–1250 [12]. The program RNAalifold [28] was used to generate a preliminary consensus sequence/structure for nt 1031–1250 based on energy minimization and structure conservation. The structure predicted without any experimental data is a multibranch loop with three hairpins, including one very thermodynamically stable long stem capped by a tetraloop (Fig 1A). When this model is compared to an alignment of 2515 unique influenza A sequences, base pairs are 91% conserved. When mutations occur, they preserve pairing via consistent (single point) and compensatory (double point) mutations (Fig 1).

thumbnail

Fig 1. Secondary structure of M121 predicted by RNAalifold and bioinformatics calculations for possibility of each base pair.

A. Base pairs in secondary structure are colored according to % of canonical base pairs calculated for type A influenza (last column of table): red >98%; 95 ≤ orange < 98%; 90 ≤ green < 95%; 85 ≤ blue < 90%. B. In table are marked preserving mutations: indicated by green shadow are evidence of compensatory mutation.

https://doi.org/10.1371/journal.pone.0141132.g001

Secondary structure mapping of segment 5 1031–1250 (+)RNA structural motif

The consensus sequence of nts 1051–1171 (+)RNA of segment 5 was synthesized (M121) for in vitro structure mapping. The possibility of alternative folding of M121 was checked by using different annealing buffers. Native gel electrophoresis revealed that M121 folds into one conformation in wide range of Mg2+ concentration, from 5 mM to 25 mM, and also in [(Co(NH3)6]3+ at concentrations from 5 to 25 mM.

RNA secondary structure was mapped in a buffer containing 100 mM KCl, 5 mM MgCl2 at 23°C and pH 7. For enzymatic mapping RNase T1 (cleaves after unpaired G), RNase S1 (cleaves after unpaired nucleotide), and RNase V1 (cleaves after paired nucleotide) were used. Small molecules CMCT (modifies N3 of U and N1 of G when unpaired), DMS (methylates N1 of A and N3 of C when unpaired), and kethoxal (modifies G when unpaired) were used to probe the Watson–Crick face of bases. Lead ion cleavage and SHAPE mapping were used to identify flexible backbone regions.

The results of chemical mapping were used to constrain minimum free energy secondary structure prediction using RNAstructure v. 5.5 [33] (Fig 2A). The experimentally constrained structure is almost identical to the RNAalifold model; only one base pair (C1088-G1093; Fig 2A) is forced to be unpaired because of a strong CMCT hit at U1094. Otherwise, all chemical mapping data agree with the RNAalifold consensus model. Furthermore, incorporating experimental data into the calculation of the base pair partition function predicts the model pairs with a high calculated probability (Fig 3). Higher probability pairs are better-defined and correlate with an increased quality of prediction [34]. In particular nts 1102–1149, 1051–1063 and 1158–1171 are predicted with a very high level of confidence.

thumbnail

Fig 2. Secondary structure of influenza RNA motif.

A. M121, isolated RNA with consensus sequence; secondary structure was predicted by RNAstructure 5.5 using chemical mapping results from SHAPE. B. Chemical mapping results of M121 in entire segment 5 (+)RNA (A/VietNam/1203/2004 (H5N1)); in blue is marked difference in sequence comparing to M121 on panel A.

https://doi.org/10.1371/journal.pone.0141132.g002

Chemical and enzymatic mapping are complementary to each other and show two the most flexible regions (Figs 2A, 4 and 5). Nucleotides 1070–1075 were the most reactive: modified by kethoxal, CMCT and NMIA and cut strongly by RNase T1 and S1—all in a hairpin loop (Fig 2A). The next, most reactive sites were nts 1124–1128, where DMS, kethoxal and NMIA hits were observed. Furthermore, Pb2+ cleavages were occurred after 1124 and 1126 nucleotides. There were RNase T1 cuts after G1124 and G1126 and S1 cuts after nucleotides 1124–1128 and no cleavage by RNase V1, which offers good support for this region being unpaired.

thumbnail

Fig 4. Enzymes mapping of M121.

All reactions were conducted in 23°C for 30 min. Lane 1 –control reaction: M121 incubated in 100 mM KCl and 5 mM MgCl2, 10 mM Tris-HCl, pH 7, for 30 min in 23°C. Lanes 2–4—RNase V1 cuts in increasing concentration of enzyme: 0.5x10-3 U/μl, 1x10-3 U/μl and 3x10-3 U/μl, respectively. Lanes 5–7—RNase T1 cuts in increasing concentration of enzyme: 0.15 U/μl, 0.25 U/μl and 0.75 U/μl, respectively. Lanes 8–10—RNase S1 cuts in increasing concentration of enzyme: 0.05 U/μl, 0.3 U/μl, 1 U/μl, respectively, Lane 11—formamide ladder, Lane 12—RNase T1 ladder.

https://doi.org/10.1371/journal.pone.0141132.g004

thumbnail

Fig 5. Lead ion cleavage of M121.

RNA was incubated with 1 mM Pb(OAc)2, 100 mM KCl and 5 mM MgCl2, 10 mM Tris-HCl pH 7 in time course: lanes 1–6 - 0, 1, 5, 15, 30 and 60 min, respectively. Lane 7—control reaction: M121 incubated in 100 mM KCl and 5 mM MgCl2, 10 mM Tris-HCl pH 7 for 60 min. Lane 8—RNase T1 ladder. Lane 9—formamide ladder.

https://doi.org/10.1371/journal.pone.0141132.g005

The predicted hairpin with looped nts 1087–1093 is slightly modified by chemicals: only A1089 and A1090 were reactive to DMS and U1094 to CMCT. Strong Pb2+ cleavage was observed after A1090 and there were no enzymes cleavages. Single stranded regions between hairpins undergo strong modification in 1081–1084 by NMIA, DMS and kethoxal. In this region strong Pb2+ cleavages were observed in region 1080–1082. This fragment is more accessible than hairpin 1085–1096. Lead ion cleavage also indicates additional flexible nucleotides: 1096–1098 and 1150, 1151 1154, 1155. These predicted single-stranded regions have only two identified chemical mapping hits (1101, 1102) as nts 1171–1132 were not able to be read-out by reverse transcription. Hairpins 1064–1080 and 1103–1147 and helix 1053-1062/1168-1160 are well confirmed in structural analysis. Additionally, mapped by NMIA, DMS and kethoxal nucleotides confirmed bulges: A1111-A1112, C1116, A1120-A1121, G1132, and A1058. After 1110 and 1116 nucleotides strong Pb2+ cleavages were observed, where 1110 is the strongest from all. (Figs 2A and 5). RNase V1 cuts indicate presence of helixes: 1053-1057/1168-1164, 1064-1068/1076-1080, 1103-1110/1140-1147 and 1117-1119/1133-1135. Small hairpin 1185–1096 and single-stranded multibranch-loop linker region between hairpins could have additional or alternative intramolecular interactions.

Hybridization to Oligonucleotide Microarrays

Accessibility of M121 to oligonucleotide binding was studied using isoenergetic microarrays. Isoenergetic microarrays are built with short penta- and hexanucleotide probes which are 2’-O-methyl-RNA with locked nucleic acid (LNA) and 2,6-diaminopurine modifications at selected positions [22, 35]. In designing probes, the roughly-equal free energies of hybridization of probe/RNA duplexes was the goal [36–38]. Additionally, all chosen modifications increase the thermodynamic stability of hybridization duplex (probe/target RNA duplex), which allows the use of short probes. For example: 2,6-diaminopurine riboside (D) substituted for adenosine allows an additional hydrogen bond to form, making D-U pairs similar in stability to G-C [36]. Based on previous studies, LNAs in selected positions also increase the thermodynamic stability of duplexes. Careful design makes all probes used close to isoenergetic (in free energies of binding) [39]. In previous microarray mapping studies, primarily RNA secondary structure dictated the binding affinity of probes [22–24, 40].

The results of hybridizations of M121 to isoenergetic microarrays for three temperatures (4, 23 and 37°C) are collected in Table 2. The binding of modified probes are almost identical in 4°C and 23°C. Higher temperature eliminated binding of several probes, leaving only the most thermodynamically stable and specific. Probes which bind to target RNA were analyzed in terms of possible alternative binding sites. Predicted, other than target site, the possible binding sites (thermodynamically stable but not fully complementary; Table 2, last column) were used for identification of the unambiguous binding sites and possible accessible regions. All binding sites (the middle nt of the hybridized region for a probe) at 23°C were marked on the M121 secondary structure (Fig 2A). Accessible sites are in single stranded regions: 1073, 1154, 1155 and 1156. Binding site 1059 is in a weakly structured regions near bulge. Probe DCLUCLCGL (symbol NL indicate specific locked nucleic acid (LNA) nucleotide) has binding site at 1105 which is helix region. It is possible, because of G and A rich contest of target, that this probe form triplex instead of disrupting the double helix. Additionally, possible binding sites are in single stranded regions: 1084–1085, 1096, 1099, 1101–1103, 1126. Binding sites at 37°C provide target regions for potentially clinically useful antisense oligonucleotides.

M121 motif in segment 5 (+)RNA of A/VietNam/1203/2004(H5N1)

The structure of M121 was checked in context of entire segment 5 (+)RNA sequence [A/VietNam/1203/2004(H5N1)]. The results of chemical mapping of nts 1051–1171 using NMIA, DMS, CMCT and kethoxal is presented on Fig 2B, where differences between this H5N1 sequence and the consensus sequence used for the in vitro construct are shown by blue letters. The binding site for the reverse transcription primer was chosen to be an appropriate distance to motif and enable to read out nts 1051–1171. In general, reactivity is in agreement with proposed secondary structure: e.g. hairpins 1064–1080 and 1103–1147 (Fig 2B). Also, medium RNase H cleavage after A1127 and A1128 in presence of DNA oligonucleotide TGTTCTC complementary to 1124–1130 show the accessibility of the hairpin loop (Fig 2B). The helix formed from nts 1053-1062/1168-1160 has several hits in the stem region: strong CMCT and medium NMIA modifications at 1056 and also reactivity at 1061 (NMIA) and 1163 (NMIA, DMS). The most noticeable difference between the truncated in vitro construct and the analogical region in the whole (+)RNA is at nts 1150–1156, which are unreactive in the (+)RNA. No alternative folds of the studied region (including pseudoknot) could explain these differences; they may be due to a tertiary interaction with another region in the (+)RNA, for example, the nearby region (nts 991–997) is complementary to 1150–1156.

Inhibitory effect of oligonucleotides targeting the M121 motif in cell culture

Influenza RNAs were previously targeted with oligonucleotides as potential antiviral agents [41–43]. Several approaches using antisense oligonucleotides, siRNAs or miRNAs were published, where target regions were generally unique, well-conserved sequences, including genome packaging signals or start codons [31, 44–46]. Common used and good target were NP (segment 5) and the three polymerase RNAs (segment 1, 2, 3) [31, 47–50]. Oligonucleotide tools blocking or causing degradation of these segments not only disturb specific mRNAs but also block the accumulation of all viral RNAs, making them especially potent as influenza inhibitors. The conserved RNA structural motif presented in this current work contains possible new targets for oligonucleotides, that may complement those already found for the NP (+)RNA. Furthermore, here we use principles of oligonucleotide design that differ from previous studies.

The importance of the M121 structural motif was investigated using an antisense oligonucleotide inhibition assay with a single-cycle infectious influenza A virus (sciIAV) [strain A/California/04_NYICE_E3/2009 (H1N1)] in MDCK-HA cells. The chosen strain differs slightly (82% sequence conservation, Fig 6) from the consensus sequence of M121: this was taken into account during oligonucleotides design to make them fully complementary. Oligonucleotide target sites were the two well-characterized hairpins at nts 1064–1080 and 1103–1147 and the single-stranded linker that possibly is involved in a tertiary interaction.

thumbnail

Fig 6.

A. (+) RNA5 motif with marked complementary region to anisense oligonucleotides. In green were marked differences in sequence between A/California/04_NYICE_E3/2009 strain and consensus sequence of M121. B. Effect of antisense oligonucleotides targeting M121 motif of scIAV A/California/04_NYICE_E3/2009 in cell line MDCK-HA. C—control; L—control with lipofectamine, R—positive control with ribavirin; N—negative control with oligonucleotide N.

https://doi.org/10.1371/journal.pone.0141132.g006

2’-O-Methylated oligonucleotides with LNA modifications: CAULUCULCAULUU (2A) and CAULUCULCAULUUGAAGC (3A) targeted the long hairpin and internal loop of M121 (region 1120–1130 and 1115–1130, respectively) (see the complementary regions, and strain sequence differences in Fig 6). These oligonucleotides were able to inhibit virus proliferation by a factor of 8.40 (0.93 LOG10 FFU/ml difference) (Fig 6). This hairpin thus appears to serve an important function in influenza A proliferation. Unmodified, 2’-O-methyl oligonucleotide CAUUCUCAUUU (1A) slightly inhibit virus: FFU/ml decrease by 1.64 times, which gives 0.16 LOG10 FFU/ml difference.

For (+)RNA5 of A/California/04_NYICE_E3/2009 the long distance interactions between regions: 1151-1157/991-997 is possible, but weaker: there are two mismatches (5’ACCCUGG/3’UGAGAAC). Targeting nts 1148–1158 with 2’-O-methyl oligonucleotides: 1B or 2B (UCCAGGGUAUU, UCCLAGGLGUALUU), respectively, did not yield significant effects and only moderate inhibition for both 1B and 2B were observed. This is in agreement with in vitro mapping results, which show partial protection of the region from chemicals and lack of accessibility by microarray. Additionally, RNase H cleavage in presence of DNA oligonucleotide TCAAGAGTG complementary to 1150–1158 resulted in only weak cutting after U1155. It is likely that nts 1151–1157 is involved in weak tertiary interactions.

2’-O-Methyl oligonucleotide AUCACUUUCUU (1C) complementary to region 1067–1077 did not influence the virus replication, suggesting that this hairpin is less important to the virus or that it is less accessible to the oligonucleotide.

Presenting here approach leading from finding conserved secondary structure motifs in virus RNA to application of antisense oligonucleotides disturbing the motif possible function would be useful in therapeutics development.

Conclusion

A bioinformatics analysis of the viral genome segment 5, which encodes nucleoprotein, revealed a conserved structural motif in the (+)RNA. The secondary structure proposed by energy minimization and comparative analysis agrees with structure predicted based on experimental data: using a 121 nucleotide in vitro RNA construct comprising an influenza A virus consensus sequence and also an entire segment 5 (+)RNA (strain A/VietNam/1203/2004 (H5N1)). The conserved motif consists of three hairpins with one being especially thermodynamically stable. The biological importance of this conserved secondary structure is supported in experiments using antisense oligonucleotides in cell line, which found that disruption of this motif led to inhibition of viral fitness. These results suggest that this conserved motif in the segment 5 (+)RNA might be a candidate for oligonucleotide-based antiviral therapy.

Author Contributions

Conceived and designed the experiments: EK MSJ. Performed the experiments: EK MSJ PM RK WNM. Analyzed the data: EK MSJ WNM. Contributed reagents/materials/analysis tools: EK RK WNM. Wrote the paper: EK MSJ WNM.

References

  1. 1.Shapshak P, Chiappelli F, Somboonwit C, Sinnott J. The Influenza Pandemic of 2009 Lessons and Implications. Molecular Diagnosis & Therapy. 2011 2011;15(2):63–81.
  2. 2.Franco-Paredes C, Hernandez-Ramos I, Del Rio C, Alexander KT, Tapia-Conyer R, Santos-Preciado JI. H1N1 Influenza Pandemics: Comparing the Events of 2009 in Mexico with those of 1976 and 1918–1919. Archives of Medical Research. 2009 Nov;40(8):669–72. pmid:20304254
  3. 3.Van Kerkhove MD, Asikainen T, Becker NG, Bjorge S, Desenclos J-C, dos Santos T, et al. Studies Needed to Address Public Health Challenges of the 2009 H1N1 Influenza Pandemic: Insights from Modeling. Plos Medicine. 2010 Jun;7(6).
  4. 4.Neumann G, Noda T, Kawaoka Y. Emergence and pandemic potential of swine-origin H1N1 influenza virus. Nature. 2009 Jun 18;459(7249):931–9. pmid:19525932
  5. 5.York A, Fodor E. Biogenesis, assembly, and export of viral messenger ribonucleoproteins in the influenza A virus infected cell. Rna Biology. 2013 Aug 1;10(8):1274–82. pmid:23807439
  6. 6.Bouvier NM, Palese P. The biology of influenza viruses. Vaccine. 2008 Sep;26:D49–D53. pmid:19230160
  7. 7.Gultyaev AP, Fouchier RAM, Olsthoorn RCL. Influenza Virus RNA Structure: Unique and Common Features. International Reviews of Immunology. 2010 Dec;29(6):533–56. pmid:20923332
  8. 8.Gultyaev AP, Olsthoorn RCL. A family of non-classical pseudoknots in influenza A and B viruses. RNA Biology. 2010 Mar-Apr;7(2):125–9. pmid:20200490
  9. 9.Gultyaev A, Tsyganov-Bodounov A, Spronken M, van der Kooij S, Fouchier R, Olsthoorn R. RNA structural constraints in the evolution of the influenza A virus genome NP segment. RNA Biology 2014;11:942–52. pmid:25180940
  10. 10.Bae SH, Cheong HK, Lee JH, Cheong C, Kainosho M, Choi BS. Structural features of an influenza virus promoter and their implications for viral RNA synthesis. Proc Natl Acad Sci U S A. 2001 Sep 11;98(19):10602–7. pmid:11553808
  1. 11.Moss WN, Dela-Moss LI, Kierzek E, Kierzek R, Priore SF, Turner DH. The 3' splice site of influenza A segment 7 mRNA can exist in two conformations: a pseudoknot and a hairpin. PloS One. 2012;7(6):e38323–e. pmid:22685560
  1. 12.Moss WN, Priore SF, Turner DH. Identification of potential conserved RNA secondary structure throughout influenza A coding regions. RNA. 2011 Jun;17(6):991–1011. pmid:21536710
  1. 13.Jiang T, Kennedy SD, Moss WN, Kierzek E, Turner DH. Secondary structure of a conserved domain in an intron of influenza A M1 mRNA. Biochemistry. 2014;53:5236–48. pmid:25026548
  1. 14.Priore SF, Kierzek E, Kierzek R, Baman J, Moss WN, Dela-Moss LI, et al. Secondary Structure of a Conserved Domain in the Intron of Influenza A NS1 mRNA. Plos One. 2013;8(9):e70615. pmid:24023714
  1. 15.Priore SF, Moss WN, Turner DH. Influenza A Virus Coding Regions Exhibit Host-Specific Global Ordered RNA Structure. Plos One. 2012 Apr;7(4):e35989. pmid:22558296
  1. 16.Dela-Moss LI, Moss WN, Turner DH. Identification of conserved RNA secondary structures at influenza B and C splice sites reveals similarities and differences between influenza A, B, and C. BMC Research Notes. 2014 2014 Jan;7:22-. pmid:24405943
  1. 17.Chiang C, Chen G-W, Shih S-R. Mutations at Alternative 5 ' Splice Sites of M1 mRNA Negatively Affect Influenza A Virus Viability and Growth Rate. J Virol. 2008 Nov;82(21):10873–86. pmid:18768984
  1. 18.Xia TB, SantaLucia J, Burkard ME, Kierzek R, Schroeder SJ, Jiao XQ, et al. Thermodynamic parameters for an expanded nearest-neighbor model for formation of RNA duplexes with Watson-Crick base pairs. Biochemistry. 1998 Oct;37(42):14719–35. pmid:9778347
  1. 19.Kierzek E, Kierzek R. The thermodynamic stability of RNA duplexes and hairpins containing N-6-alkyladenosines and 2-methylthio-N-6-alkyladenosines. Nucleic Acids Res. 2003 Aug;31(15):4472–80. pmid:12888507
  1. 20.Ehresmann C, Baudin F, Mougel M, Romby P, Ebel JP, Ehresmann B. Probing the structure of RNAs in solution. Nucleic Acids Res. 1987;15:9109–28. pmid:2446263
  1. 21.Merino EJ, Wilkinson KA, Coughlan JL, Weeks KM. RNA structure analysis at single nucleotide resolution by selective 2 '-hydroxyl acylation and primer extension (SHAPE). J Am Chem Soc. 2005 Mar;127(12):4223–31. pmid:15783204
  1. 22.Kierzek E, Kierzek R, Moss WN, Christensen SM, Eickbush TH, Turner DH. Isoenergetic penta- and hexanucleotide microarray probing and chemical mapping provide a secondary structure model for an RNA element orchestrating R2 retrotransposon protein function. Nucleic Acids Res. 2008 Apr;36(6):1770–82. pmid:18252773
  1. 23.Kierzek E. Binding of short oligonucleotides to RNA: studies of the binding of common RNA structural motifs to isoenergetic microarrays. Biochemistry. 2009 Dec;48(48):11344–56. pmid:19835418
  1. 24.Fratczak A, Kierzek R, Kierzek E. Isoenergetic Microarrays To Study the Structure and Interactions of DsrA and OxyS RNAs in Two- and Three-Component Complexes. Biochemistry. 2011 Sep;50(35):7647–65. pmid:21793590
  1. 25.Bao Y, Bolotov P, Dernovoy D, Kiryutin B, Zaslavsky L, Tatusova T, et al. The influenza virus resource at the national center for biotechnology information. J Virol. 2008 Jan;82(2):596–601. pmid:17942553
  1. 26.Gruber AR, Findeiss S, Washietl S, Hofacker IL, Stadler PF. RNAz 2.0: improved noncoding RNA detection. Pac Symp Biocomp. 2010;15:69–79.
  1. 27.Gruber AR, Neuboeck R, Hofacker IL, Washietl S. The RNAz web server: prediction of thermodynamically stable and evolutionarily conserved RNA structures. Nucleic Acids Res. 2007 Jul;35:W335–W8. pmid:17452347
  1. 28.Bernhart SH, Hofacker IL, Will S, Gruber AR, Stadler PF. RNAalifold: improved consensus structure prediction for RNA alignments. Bmc Bioinformatics. 2008 Nov 11;9:474. pmid:19014431
  1. 29.Martinez-Sobrido L, Cadagan R, Steel J, Basler CF, Palese P, Moran TM, et al. Hemagglutinin-Pseudotyped Green Fluorescent Protein-Expressing Influenza Viruses for the Detection of Influenza Virus Neutralizing Antibodies. J Virol. 2010 Feb 15;84(4):2157–63. pmid:19939917
  1. 30.Baker SF, Guo H, Albrecht RA, Garcia-Sastre A, Topham DJ, Martinez-Sobrido L. Protection against Lethal Influenza with a Viral Mimic. J Virol. 2013 Aug;87(15):8591–605. pmid:23720727
  1. 31.Zhang T, Wang T-C, Zhao P-S, Liang M, Gao Y-W, Yang S-T, et al. Antisense oligonucleotides targeting the RNA binding region of the NP gene inhibit replication of highly pathogenic avian influenza virus H5N1. International Immunopharmacology. 2011 Dec;11(12):2057–61. pmid:21933722
  1. 32.Washietl S, Hofacker IL, Stadler PF. Fast and reliable prediction of noncoding RNAs. Proc Natl Acad Sci U S A. 2005 Feb 15;102(7):2454–9. pmid:15665081
  1. 33.Reuter JS, Mathews DH. RNAstructure: software for RNA secondary structure prediction and analysis. Bmc Bioinformatics. 2010 Mar;11.
  1. 34.Mathews DH, Disney MD, Childs JL, Schroeder SJ, Zuker M, Turner DH. Incorporating chemical modification constraints into a dynamic programming algorithm for prediction of RNA secondary structure. Proc Natl Acad Sci U S A. 2004 May;101(19):7287–92. pmid:15123812
  1. 35.Kierzek E, Christensen SM, Eickbush TH, Kierzek R, Turner DH, Moss WN. Secondary structures for 5' regions of R2 retrotransposon RNAs reveal a novel conserved pseudoknot and regions that evolve under different constraints. J Mol Biol. 2009;390:428–42. pmid:19397915
  1. 36.Pasternak A, Kierzek E, Pasternak K, Turner DH, Kierzek R. A chemical synthesis of LNA-2,6-diaminopurine riboside, and the influence of 2'-O-methyl-2,6-diaminopurine and LNA-2,6-diaminopurine ribosides on the thermodynamic properties of 2'-O-methyl RNA/RNA heteroduplexes. Nucleic Acids Res. 2007;35(12):4055–63. pmid:17567607
  1. 37.Pasternak A, Kierzek E, Pasternak K, Fratczak A, Turner DH, Kierzek R. The thermodynamics of 3'-terminal pyrene and guanosine for the design of isoenergetic 2'-O-methyl-RNA-LNA chimeric oligonucleotide probes of RNA structure. Biochemistry. 2008 Feb;47(5):1249–58. pmid:18177022
  1. 38.Kierzek E, Ciesielska A, Pasternak K, Mathews DH, Turner DH, Kierzek R. The influence of locked nucleic acid residues on the thermodynamic properties of 2'-O-methyl RNA/RNA heteroduplexes. Nucleic Acids Res. 2005;33(16):5082–93. pmid:16155181
  1. 39.Kierzek E, Fratczak A, Pasternak A, Turner DH, Kierzek R. Isoenergetic RNA microarrays, a new method to study the structure and interactions of RNA. International Proceedings Division. 2007:25–30.
  1. 40.Kierzek E, Kierzek R, Turner DH, Catrina IE. Facilitating RNA structure prediction with microarrays. Biochemistry. 2006 Jan;45(2):581–93. pmid:16401087
  1. 41.Mizuta T, Fujiwara M, Hatta T, Abe T, Miyano-Kurosaki N, Shigeta S, et al. Antisense oligonucleotides directed against the viral RNA polymerase gene enhance survival of mice infected with influenza A. Nature Biotechnology. 1999 Jun;17(6):583–7. pmid:10385324
  1. 42.Ge Q, Pastey M, Kobasa D, Puthavathana P, Lupfer C, Bestwick RK, et al. Inhibition of multiple subtypes of influenza A virus in cell cultures with morpholino oligomers. Antimicrob Agents Chemother. 2006 Nov;50(11):3724–33. pmid:16966399
  1. 43.Wu Y, Zhang G, Li Y, Jin Y, Dale R, Sun L-Q, et al. Inhibition of highly pathogenic avian H5N1 influenza virus replication by RNA oligonucleotides targeting NS1 gene. Biochem Biophys Res Commun. 2008;365:369–74. pmid:17996729
  1. 44.Kumar B, Khanna M, Kumar P, Sood V, Vyas R, Banerjea AC. Nucleic Acid-Mediated Cleavage of M1 Gene of Influenza A Virus Is Significantly Augmented by Antisense Molecules Targeted to Hybridize Close to the Cleavage Site. Molecular Biotechnology. 2012 May;51(1):27–36. pmid:21744034
  1. 45.Li W, Yang X, Jiang Y, Wang B, Yang Y, Jiang Z, et al. Inhibition of influenza A virus replication by RNA interference targeted against the PB1 subunit of the RNA polymerase gene. Archives of Virology. 2011 Nov;156(11):1979–87. pmid:21837415
  1. 46.Ge Q, Filip L, Bai AL, Nguyen T, Eisen HN, Chen J. Inhibition of influenza virus production in virus-infected mice by RNA interference. Proc Natl Acad Sci U S A. 2004 Jun 8;101(23):8676–81. pmid:15173599
  1. 47.Ge Q, McManus MT, Nguyen T, Shen CH, Sharp PA, Eisen HN, et al. RNA interference of influenza virus production by directly targeting rnRNA for degradation and indirectly inhibiting all viral RNA transcription. Proc Natl Acad Sci U S A. 2003 Mar;100(5):2718–23. pmid:12594334
  1. 48.Martelli F, Salata C, Calistri A, Parolin C, Azzi A, Palù G, et al. Small RNAs targeting the 5' end of the viral polymerase gene segments specifically interfere with influenza type A virus replication. J Biotechnol 2015;210:85–90. pmid:26091771
  1. 49.Zhang T, Zhao P-S, Zhang W, Liang M, Gao Y-W, Yang S-T, et al. Antisense oligonucleotide inhibits avian influenza virus H5N1 replication by single chain antibody delivery system. Vaccine. 2011 Feb 11;29(8):1558–64. pmid:21215341
  1. 50.Giannecchini S, Wise HM, Digard P, Clausi V, Del Poggetto E, Vesco L, et al. Packaging signals in the 5 '-ends of influenza virus PA, PB1, and PB2 genes as potential targets to develop nucleic-acid based antiviral molecules. Antiviral Res. 2011 Oct;92(1):64–72. pmid:21741410