Exploiting functional regions in the viral RNA genome as druggable entities
eLife Assessment
The study presents a fundamental advance in antiviral RNA research by adapting SHAPE-Map to chart the secondary structure of the porcine epidemic diarrhea virus (PEDV) genome in infected cells and pinpointing structurally conserved, accessible RNA elements as therapeutic targets. A broad, well-documented integration of biochemical probing, computational analysis, and functional validation provides convincing evidence that these regions are both biologically relevant and druggable. Beyond PEDV, the work offers a generalizable framework for RNA-guided antiviral discovery that will interest researchers in RNA therapeutics and viral genome biology.
https://doi.org/10.7554/eLife.103923.3.sa0Fundamental: Findings that substantially advance our understanding of major research questions
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Convincing: Appropriate and validated methodology in line with current state-of-the-art
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Abstract
RNA-targeting compounds or small interfering RNAs (siRNAs) offer a potent means to control viral infections. An essential prerequisite for their design depends on identifying conserved and functional viral RNA structures in cells. Techniques that probe RNA structures in situ have recently been developed including SHAPE-MaP, which has been helpful in the analysis of secondary structures of RNA. In this study, we report on the application of SHAPE-MaP to the porcine epidemic diarrhea virus RNA genome to categorize different functional regions, including potential quadruplex-forming sequence and target sites of siRNA. Our results show that these structures can be exploited to inhibit viral proliferation and that SHAPE-MaP is an effective method to identify secondary structures in RNA genomes.
Introduction
RNA-targeting drugs are expected to address the challenges posed by ‘undruggable’ protein targets and significantly expand the landscape of targetable macromolecules (Disney, 2019, Kovachka et al., 2024, Li and Wang, 2019, Ursu et al., 2020, Warner et al., 2018). The US Food and Drug Administration (FDA) has approved six siRNAs (Tang and Khvorova, 2024), nine antisense oligonucleotides (ASOs) (Bennett, 2019, Chen et al., 2024, Crooke et al., 2021), and the RNA-targeting compound Risdiplam (Ratni et al., 2021, Sheridan, 2021) for treating complicated diseases such as spinal muscular atrophies. The selection of RNAs that are functional, conserved, and accessible to be targeted in vivo is fundamental to the development of RNA-targeting drugs. Unlike RNAs in vitro, the folding of RNAs in cells is affected by the intricate intracellular environment (Manfredonia et al., 2020, Zhang et al., 2021b). Thus, in situ probing for the accessibility of viral RNAs by antivirals is essential for the identification of targetable RNAs.
In recent years, the integration of chemical probing methods and high-throughput sequencing technologies has been used to analyze the structure of the viral RNA genome in infected cells (Boerneke et al., 2019, Cao et al., 2021, Flynn et al., 2016, Huber et al., 2019, Huston et al., 2021, Lan et al., 2022, Li et al., 2018, Manfredonia et al., 2020, Mauger et al., 2015, Pollom et al., 2013, Siegfried et al., 2014, Smola et al., 2016, Yang et al., 2021, Zhang et al., 2021a, Ziv et al., 2020). These methods have confirmed the widespread presence of structurally organized RNA elements within the viral genome, which are composed of various types of RNA structures such as stem loop structures (with double-helical stems), RNA single-strands, RNA pseudoknots, and RNA G-quadruplexes (G4s), collectively forming a complex and dynamic global genome architecture (Boerneke et al., 2019, Kwok et al., 2016). Of these, RNA G4 has shown considerable promise due to its high stability and ability for modulation by small molecules (Zhao et al., 2021, Fang et al., 2023, Fleming et al., 2016, Kwok et al., 2016, Ruggiero and Richter, 2018, Wang et al., 2016, Zhang et al., 2021a). Furthermore, RNA interference (RNAi) has been shown to inhibit gene expression and viral replication (Bennett, 2019, Bowden-Reid et al., 2023, Piasecka et al., 2020, Qureshi et al., 2018, Sagan et al., 2010, Tang and Khvorova, 2024, Westerhout and Berkhout, 2007). Despite advances, it remains a challenge to identify RNA elements or regions that serve as targets within complex genomic structural contexts.
Selective 2’-hydroxyl acylation and primer extension mutational profile (SHAPE-MaP) stands out in ex vivo RNA structure probe technologies for its ability to quantify nucleotide flexibility and solvent accessibility at single nucleotide resolution (Lorenz et al., 2016, Siegfried et al., 2014, Smola et al., 2016, Lucks et al., 2011). In short, SHAPE reagents like NAI selectively modify flexible unpaired 2’-OH groups in RNA, and these modifications are detected as mutations during reverse transcription, allowing precise mapping of RNA secondary structures by sequencing (Siegfried et al., 2014). In the viral long-stranded RNA genome, SHAPE reactivities are used as constraints in the prediction of RNA structures based on nearest-neighbor rules (Deigan et al., 2009, Mathews, 2004a), which can result in thousands of possible conformations. The pairing probability of each nucleotide derived from the SHAPE reactivities was subsequently used to calculate the Shannon entropy. Regions with high Shannon entropy can adopt alternative conformations, while those with low Shannon entropy correspond to well-defined RNA structures or persistently single-stranded regions (Mathews et al., 2004b, Siegfried et al., 2014). Therefore, SHAPE reactivity provides information about the pairing probability and accessibility of nucleotides, whereas Shannon entropy reflects the likelihood of alternative conformations in a given region. RNA motifs with low SHAPE reactivity and low Shannon entropy are assumed to be well-folded structures and potential targets (Boerneke et al., 2019). In this study, we categorized RNAs with distinct features in the viral genome according to SHAPE reactivity and Shannon entropy and investigate their potential as targets (Figure 1).

Four regions with different characteristics in the porcine epidemic diarrhea virus (PEDV) genome.
(A) Schematic of SHAPE-MaP for probing the RNA structure of the PEDV genome in situ. (B) Well-folded regions (low SHAPE reactivity and low Shannon entropy; 26.40% of genome). These regions represent stably folded RNA structures with minimal conformational flexibility, serving as structural scaffolds or functional elements in viral replication. (C) Dynamic structured regions (low SHAPE reactivity and high Shannon entropy; 11.70% of genome). These conformationally plastic domains mediate regulatory switches between alternative secondary structures during infection. (D) Dynamic unpaired regions (high SHAPE reactivity and high Shannon entropy; 26.90% of genome). These regions are prone to form non-canonical nucleic acid structures (e.g., G-quadruplexes), which can be stabilized by small-molecule ligands to inhibit viral replication. (E) Persistent unpaired regions (high SHAPE reactivity and low Shannon entropy; 9.67% of genome). These regions are more accessible for siRNA binding, facilitating recruitment of Argonaute proteins and Dicer to form the RNA-induced silencing complex (RISC) for targeted cleavage.
The porcine epidemic diarrhea virus (PEDV) is an alpha-coronavirus with a positive single-stranded RNA genome of 28 kb (Jung et al., 2020, Lee, 2015). It results in severe diarrhea and vomiting, leading to the death of 90% of infected piglets. The spread of PEDV is worldwide and its infection is of primary concern in pig husbandry (Lee, 2015). The epizootic reemergence of PEDV has also been reported to infect human intestinal cells and shows a potential susceptibility to infection by different species (Niu et al., 2023).
In this study, we adapted the SHAPE-MaP method to detect the in situ structure of the PEDV RNA genome in infected cells. More specifically, the combination of SHAPE reactivity and Shannon entropy categorized the PEDV RNA genome into different functional types (Figure 1), based on distinctive profiles, including the 5’ untranslated region (5’ UTR), frameshifting stimulatory element (FSE), and 3’ untranslated region (3’ UTR). We particularly focused on guanine-rich regions, which could potentially fold into a stable G4 structure and sites suitable for binding small interfering RNA (siRNA). Our results show that G4 folding in a dynamic single-stranded region with high SHAPE reactivity and high Shannon entropy inhibited viral proliferation, while siRNA targeting persistent single-stranded regions with high SHAPE reactivity and low Shannon entropy exhibited greater success rates in inhibiting viral replication. Finally, this work proposes a strategy for selecting targetable RNAs based on the in situ folding characteristics of the viral genome.
Results
Mapping the RNA structures in the PEDV genome
SHAPE-MaP was used to characterize the structural features of the PEDV RNA genome in infected cells (Figure 1A). The high correlations of mutational signatures for biological replicates (ex vivo SHAPE R2=0.955; Figure 1—figure supplement 1A) and the high mutation rate with NAI treatment (Figure 1—figure supplement 1B) indicate the high quality of the sequencing data and the experimental workflow. Hence, two replicates were combined for downstream analyses.
As 5’UTR with conserved functions across coronaviruses has been extensively reported (Chen and Olsthoorn, 2010, Madhugiri et al., 2018, Yang and Leibowitz, 2015), we first examined the structure of this region in the PEDV genome for its stability ex vivo. The PEDV 5’ UTR structure modeled with SHAPE reactivity constraints contains four stem loops (SL1, SL2, SL4, and SL5; Figure 1—figure supplement 1C) similar to other alpha-coronaviruses (Chen and Olsthoorn, 2010, Madhugiri et al., 2018, Yang and Leibowitz, 2015). Our SHAPE dataset showed that 95.6% of highly reactive bases (SHAPE reactivity ≥0.7) are not involved in canonical base pairing or are located adjacent to helical termini or bulges/loops (Figure 1—figure supplement 1C). This result suggested that our dataset is reliable as a constraint to predict the intracellular genomic structure of PEDV.
To characterize the intracellular folding features of the PEDV genome, the median SHAPE reactivities and Shannon entropies were calculated in sliding, centered windows (window = 50 nt, step = 1 nt) and a full genome SHAPE structural map was constructed accordingly (Figure 2). The global median SHAPE reactivity of the genome is 0.233. Specifically, within the 1–8000 nt region at the 5’ end, the median reactivity is 0.173; it increases to 0.268 in the middle segment of the genome (8000–24,000 nt) and then decreases to 0.194 at the 3’ end (24,000–28,000 nt). These results indicate that the folding of PEDV genomic RNA within infected cells is nonuniform, with the 5’ and 3’ ends more compactly structured, while the central region is more flexible (Figure 2).

SHAPE structural map of the porcine epidemic diarrhea virus (PEDV) genome in infected cells.
From top to bottom: SHAPE reactivity, Shannon entropy, base-pairing probabilities, translation reading frames are indicated by arrows, and genomic coordinates are marked at the bottom. Well-folded regions with low SHAPE reactivity and low Shannon entropy are shaded in blue. Potential G4 forming sequences (PQSs) are marked with a black arrow.
Within the coding regions of the ORF1a, ORF3, M, and N genes, both SHAPE reactivity and Shannon entropy are below the global median, suggesting that these regions are stably folded in the genome (Figure 2—figure supplement 1). In contrast, the coding regions of the S and E genes exhibit more relaxed and variable conformations, with both SHAPE reactivity and Shannon entropy exceeding the global median (Figure 2—figure supplement 1). Among the known functional regions, the 5’ UTR (position: 1–310 nt), which is involved in translation initiation and genome packaging, exhibited lower SHAPE reactivity and lower Shannon entropy compared to the global median, indicating a well-folded overall structure. The FSE (position 12605–12690), which is involved in mediating the ribosomal frame shift, exhibited high SHAPE reactivity and low Shannon entropy values. The 3’ UTR (position 27,702–28,044 nt), which is responsible for genome cyclization, mRNA stability, and the initiation of negative-strand synthesis, showed low SHAPE reactivity and high Shannon entropy values. These functional regions in the PEDV genome exhibit distinct SHAPE reactivity, and Shannon entropy profiles suggest that these regions may adopt different structural folding patterns to perform specific functions (Figure 1).
Functional categorization of PEDV genome reveals potentially targetable RNA
To further analyze other functional regions through SHAPE reactivity and Shannon entropy, we combined these two parameters for the classification of structural features (see ‘Materials and methods’). Low SHAPE regions were defined as windows in which >75% of the bases exhibited SHAPE reactivities below the global median (calculated in the full PEDV genome), and low Shannon entropy regions were defined according to the same rule, following Manfredonia’s method (Manfredonia et al., 2020). Windows in which >75% of their bases exhibiting values above the global median were picked as high SHAPE or high Shannon regions. Finally, four categories of RNAs in the PEDV genome were defined according to the relative values of the SHAPE reactivity and Shannon entropy (Figure 1, Supplementary files 2 and 3).
Regions with high base-pairing rates (low SHAPE) accounted for 38.10% of the PEDV genome, including the 5’ UTR and 3’ UTR (Figure 1B and C). Of these, the low SHAPE-low Shannon regions (26.40%) maintain stable conformations throughout the viral life cycle and are assumed to represent the well-folded skeleton and the potentially functional regions of the viral genome (Figure 1B). Sixty well-folded structures were identified, which exhibited consistent conformations in two replicates (Figure 2—figure supplements 2–11). The low SHAPE-high Shannon regions (11.65%), which exhibit variable conformations, might perform different functions through conformational transitions (Figure 1C). The design of drugs targeting these regions would require the identification of the specific conformations that are responsible for these functions.
Regions with low basepairing rates (high SHAPE) represented 36.57% of the PEDV genome, including the TRS-Ls of the 5’ UTRs (positions 57–72 nt) and FSE regions (positions 12,605–12,690 nt) (Figure 1D and E). RNAs in the high SHAPE-high Shannon regions (26.90%) are prone to conformational changes, and their high accessibility suggests that they are more susceptible to interaction with small molecule drugs (Figure 1D). In contrast, RNAs in the high SHAPE-low Shannon regions (9.67% of the PEDV genome) maintain a stable single-stranded state within the cell, facilitating pairing with exogenous complementary sequences (Figure 1E). Consequently, these regions are considered ideal targets in the design of antiviral oligonucleotides (Manfredonia et al., 2020, Piasecka et al., 2020, Sagan et al., 2010).
Not all folded functional structures are potential antiviral targets
Although well-folded structures are thought to be potential targets for antiviral molecules (Childs-Disney et al., 2022, Sreeramulu et al., 2021, Huston et al., 2021, Manfredonia et al., 2020), there have been relatively limited experimental validations. The 5’UTR-SL5 of coronaviruses is required for defective interfering (DI) RNA replication (Brown et al., 2007) and is important for viral packaging (Masters, 2019). In the PEDV genome, the 5’UTR-SL5 with a four-way junction structure exhibits well-folded characteristics, similar to those of SARS-CoV-2 (Figure 2—figure supplement 12). Seven compounds (Supplementary file 4, compounds 1–7), previously reported to interact with 5’ UTR-SL5 of SARS-CoV-2 (Sreeramulu et al., 2021), were found to bind with the 5’ UTR-SL5 of PEDV in vitro, and two compounds (compounds 1 and 4) exhibited dissociation constants (KD) at the µM-level (Figure 2—figure supplement 13). However, none of these compounds inhibited the proliferation of PEDV at 50 µM in cells (Figure 2—figure supplement 14). This result suggests that compounds bound to well-folded functional structures in vitro do not necessarily exhibit antiviral activity in cells. The accessibility of compounds to potential RNA targets could be an important issue, suggesting that the targetability of the 60 well-folded RNA structures in the PEDV genome requires further validation.
Identification of PQSs in dynamic single-stranded regions
Dynamic RNAs in the high SHAPE-high Shannon regions can be induced by compounds to form stable structures that perform specific biological functions (Ratni et al., 2021). We identified 49 regions with high SHAPE reactivity and high Shannon entropy with at least 25 nucleotides in length (Supplementary file 3), spanning a total of 1748 bases (6.2% of the PEDV genome).
To investigate whether RNA motifs in high SHAPE-high Shannon regions can be stabilized to inhibit viral proliferation, we focused on the G4 structure in these regions. Nine putative PQSs were identified by at least three of the G4 prediction tools (Figure 3A, Figure 3—figure supplement 1). Given the abundance of PEDV mutant strains in nature, further conservation analyses identified seven PQSs that are highly conserved (>90%) among 761 PEDV strains as potential functional G4 candidates (PQS1–PQS7, Supplementary file 5). PQS2, PQS3, PQS4, and PQS6 form stable double-stranded structures with low SHAPE reactivity and low Shannon entropy, making them less likely to fold into stable G4 structures (Figure 3—figure supplements 2–5). PQS5 and PQS7 did not fall into our four defined regions and were also omitted from further investigation. In contrast, PQS1 is located in the dynamic single-stranded region with high SHAPE reactivity and high Shannon entropy (Figure 3B and C), with little competition from local hairpin structure folding.

The PQS1 within high SHAPE-high Shannon regions has the potential to form G4 structure in the porcine epidemic diarrhea virus (PEDV) genome.
(A) Distribution of PQSs in the PEDV genome. (B) SHAPE reactivity and Shannon entropy for the regions containing PQS1. (C) Local secondary structure of the region containing PQS1 predicted with SHAPE reactivity constraints. PQS1 is marked with a blue dashed box. (D) Fluorescence turn-on assays of ThT (1 µM) in the presence of PQS1 (0.5 µM), PQS1-long chain (0.5 µM), and PQS1mut-long chain (0.5 µM). (E) Fluorescence turn-on assays of ThT (1 µM) in the presence of PQS3 (0.5 µM) and PQS3-long chain (0.5 µM). Excitation wavelength was 442 nm. PRRSV-G4 RNA was used as control.
To validate whether SHAPE analysis could reflect the competitive conformational folding of PQSs in the PEDV genome, we performed in vitro transcription to obtain local intact structures containing PQSs within dynamic single-stranded regions and stable double-stranded regions (Supplementary file 6). Thioflavin T (ThT) fluorescence turn-on assays were performed under physiological K+ conditions (100 mM), with the G4 sequence of porcine reproductive and respiratory syndrome virus (PRRSV) serving as a positive control (Control-G4) (Fang et al., 2023). The results demonstrated that for short PQSs sequences containing only G4-forming motifs (Supplementary file 7), PQS1, PQS3, PQS4, and PQS6 all induced significant enhancement of ThT fluorescence (Figure 3D and E, Figure 3—figure supplement 6), confirming their ability to form G4 structures. However, in long RNA fragments that encompass PQSs and their flanking sequences (PQSs-long chain), only PQS1 and PQS4 exhibited pronounced ThT fluorescence responses (Figure 3D and E), while PQS2, PQS3, and PQS6 showed negligible signals (Figure 3E, Figure 3—figure supplement 6). Notably, the PQS1-long chain displayed the strongest fluorescence signal, while its mutant counterpart (PQS1mut-long chain) exhibited the lowest background fluorescence (Figure 3D). These findings indicate that although most PQSs can form G4 structures in vitro, PQS1—located in the high SHAPE-high Shannon entropy region—demonstrates the most robust G4-forming capability when competing with local secondary structures in the genomic context. Therefore, PQS1 was selected for further structural and functional validation.
PQS1 is located in the nsp3 gene of the PEDV genome (Figure 3F) and is highly conserved in the 761 PEDV strains (98.37%, Figure 3—figure supplement 7A). Native polyacrylamide gel electrophoresis (PAGE) and circular dichroism (CD) measurements were performed to further verify the formation of the PQS1 G4 structure. The migration rate of PQS1 accelerated noticeably with increasing K+ concentration (Figure 3—figure supplement 7B), while the migration rate of mutant sequences (PQS1m) that did not form G4 remained constant. CD spectra revealed that PQS1 adopted a parallel G4 structure in the presence of K+, with the corresponding typical features of a positive CD band at 265 nm and a negative peak near 240 nm (Figure 3—figure supplement 7C). 1H nuclear magnetic resonance (NMR) was used to further characterize the G4 structure of PQS1. Chemical shifts in the 10–12 ppm range are considered hallmarks of G4 structures (Adrian et al., 2012). The 1H NMR spectrum of PQS1 showed prominent imino proton peaks in this region, confirming the formation of the PQS1 G4 structure (Figure 4A).

The G-quadruplex structure, biological functions of PQS1, and antiviral effects of Braco-19.
(A) 1H nuclear magnetic resonance (NMR) analysis of PQS1. (B) Circular dichroism (CD) melting profiles of PQS1. (C) Quantitative fluorescence signal using the corrected total cell fluorescence method for EGFP in cells transfected with plasmids containing the empty vector, PQS1 and PQS1mut. (D) Proliferation curve of the PEDV wild type (WT) strain and PQS1 mutant strain. (E) The relative inhibition rates of Braco-19 against AJ1102-WT and AJ1102-PQS1mut. (F) Western blot analysis of the effects of Braco-19 on the viral N protein expression of AJ1102-WT and AJ1102-PQS1mut.
-
Figure 4—source data 1
Original western blots for Figure 4F, indicating the relevant bands and treatments.
- https://cdn.elifesciences.org/articles/103923/elife-103923-fig4-data1-v1.zip
-
Figure 4—source data 2
Original files for western blot analysis displayed in Figure 4F.
- https://cdn.elifesciences.org/articles/103923/elife-103923-fig4-data2-v1.zip
Stabilization of PQS1 as a G-quadruplex inhibits RNA synthesis, gene expression, and viral proliferation
The addition of the G4-specific ligand Braco-19 significantly increased the melting temperature (Tm) of PQS1 by 10.8℃ (Figure 4B), indicating that the RNA G4 structure of PQS1 is stabilized by Braco-19. Furthermore, RNA stop assays and EGFP expression inhibition experiments highlighted the effect of a stable G4 structure of PQS1 in inhibiting RNA replication (Figure 4—figure supplement 1) and protein expression (Figure 4C, Figure 4—figure supplement 2).
To explore whether PQS1 could be a potential drug target, we constructed a recombinant virus without PQS1 (AJ1102-PQS1mut) through a synonymous mutation (G3109A) that does not change the amino acid code (Figure 4—figure supplement 3A). The level of AJ1102-PQS1mut proliferation was significantly faster than that of the wild-type strain rAJ1102, suggesting that the disruption of the PQS1 folding facilitates viral proliferation (Figure 4D). Another potential G4 sequence, PQS3, located in a stable double-stranded region with low SHAPE reactivity and low Shannon entropy, was used as a control. The PQS3 mutant recombinant virus (AJ1102-PQS3mut) proliferated at the same level as the wild type after synonymous mutation (G12322A) of its G4 sequence (Figure 4—figure supplement 3B).
Braco-19 inhibited PEDV proliferation with a median effective inhibitory concentration (EC50) of 3.60 µM (Figure 4—figure supplement 4A), which is well below the median cytotoxic concentration (CC50 > 100 µM) of Braco-19 (Figure 4—figure supplement 4B). In contrast, no inhibition was observed for Braco-19 against the mutant strain AJ1102-PQS1mut (Figure 4E). Western blot analysis was used to determine the effects of Braco-19 on viral protein expression in infected cells using the PEDV N protein antibody. The production of the AJ1102-WT N protein was significantly reduced with increasing concentrations of Braco-19, while no change was observed in the production of the AJ1102-PQS1mut N protein (Figure 4F), indicating that Braco-19 inhibited the PEDV protein expression in cells by targeting PQS1. Furthermore, as a control, we observed nearly identical inhibitory activity of Braco-19 against both the PQS3 mutant strain (AJ1102-PQS3mut) and the wild-type virus (Figure 4—figure supplement 3C), demonstrating the specificity of Braco-19’s action on PQS1.
Efficient antiviral siRNAs design targeting stable single-stranded regions
High SHAPE-low Shannon regions represent stable, accessible RNAs in cells, exposing long Watson–Crick edges of multiple residues available for base-pairing with complementary sequences (Figure 1D). Our analysis of fold characterization (see ‘Materials and methods’) identified 73 such regions (Supplementary file 3) with at least 25 nucleotides in length, spanning a total of 2696 bases (9.6% of the PEDV genome). To test the suitability of these regions as targets for antiviral siRNAs, we designed four siRNAs targeting these regions that (i) have at least 10 consecutive unpaired nucleotides and (ii) were highly conserved in the different PEDV strains (Supplementary file 7, Figure 5A–D). As a control, we designed four additional highly conserved siRNAs targeting duplex regions (Figure 5—figure supplement 1) and one non-targeting siRNA that does not target any sequence of the viral and host genomes (Supplementary file 7).

Secondary structure of target regions and antiviral effects of siRNAs.
(A–D) Local secondary structures of the stable single-stranded regions targeted by siRNAs predicted by SHAPE reactivity as a constraint; (A) ss-1; (B) ss-2; (C) ss-3; (D) ss-4. (E) qPCR showing the relative abundance of the porcine epidemic diarrhea virus (PEDV) RNA genome in infected Vero cells. The four siRNAs targeting the high SHAPE-low Shannon regions and the four siRNAs targeting the duplex regions are labeled ss-1 to ss-4 and ds-1 to ds-4, respectively. ss (single-stranded targeting siRNAs); ds (dual-stranded targeting siRNAs). si-NC was a control siRNA that did not target any viral or host sequences, and the mock group was not inoculated with virus. (F) TCID50 assays for detecting virus titers. The presented results represent the means and standard deviations of data from three independent experiments. ns: no significant difference. *p<0.05; **p<0.01; ***p<0.001, Duncan’s multiple comparison test.
The results of RT–qPCR of infected cells 16 h after viral exposure revealed that four siRNAs targeting high SHAPE-low Shannon regions significantly inhibited viral proliferation in cells (Figure 5E) compared to the non-targeting siRNA control (siRNA-NC). Consistent with the RT-qPCR results, compared to siRNA-NC, we observed that these four siRNAs had minimal cytopathic effects (CPEs) in Vero cells (Figure 5—figure supplement 2). Further TCID50 analysis demonstrated that three of the four tested siRNAs targeting single-stranded regions reduced the TCID50 by 1.5 log10 titer compared to siRNA-NC (ss-1, ss-2, and ss-3; Figure 5F). In contrast, only one of the four siRNAs targeting duplexes could slightly prevent PEDV proliferation in cells (ds-4; Figure 5E and F, Figure 5—figure supplement 2). The target sequences of four single-stranded targeting siRNAs are more than 85% conserved in 761 known PEDV strains. The combination of two or more of these siRNAs is expected to be effective in inhibiting nearly all known PEDV strains. These results showed that siRNAs targeting single-stranded regions with high SHAPE reactivity and low Shannon entropy can be more effective in exerting their antiviral effects.
Discussion
Several RNA-targeting antiviral strategies have highlighted the potential for therapeutic application (Dai et al., 2024). However, it is challenging to identify these RNAs from the viral genome in situ. In SHAPE-MaP analysis, regions with low reactivity and low entropy values are more likely to form a single stable structure in cells and are therefore considered candidates for potential functional RNAs and targets for antiviral therapy. Novel well-defined RNA structures that are critical for viral fitness in HIV, HCV, DENV, and SINV have been identified by the analysis of low SHAPE and low entropy regions (Dethoff et al., 2018, Kutchko et al., 2018, Siegfried et al., 2014, Mauger et al., 2015). However, there are still no clinically approved small molecule drugs or oligonucleotides targeting stable folded RNA structures. Our data show that two compounds interacted with the well-folded 5’UTR-SL5 of PEDV in vitro at µM-level KD, but did not show antiviral activity at 50 µM in cells (Figure 2—figure supplements 13 and 14). This result suggests that, on the one hand, the low accessibility of the regions with low SHAPE reactivity makes them difficult to access by compounds in cells and may require compounds to have a better KD with RNA. However, the binding of compounds to the well-folded structure does not necessarily alter its function. In contrast, the FDA-approved compound Risdiplam exerts its biological activity by stabilizing transient RNA structures (Campagne et al., 2019, Ratni et al., 2021, Sheridan, 2021), suggesting that it may be possible to inhibit viral proliferation by stabilizing variable conformations. We therefore focused our attention on highly accessible regions with high SHAPE reactivity in cells.
We identified 42 low SHAPE-high Shannon regions and 49 high SHAPE-high Shannon regions in the PEDV genome by SHAPE-MaP. These regions represent conformationally variable RNAs in the viral genome. As the low SHAPE regions are highly structured, it is difficult for the G4 in these regions to outcompete the folding of long base-paired structures (Figure 6A). Therefore, we selected a potential G4 sequence (PQS1) located in the high SHAPE-high Shannon region of nsp3 (PLpro), where a G4-disrupting mutation promotes viral proliferation. We found that Braco-19 significantly stabilized the G4 structure of PQS1 in vitro and inhibited PEDV proliferation in cells. Crucially, Braco-19 showed no inhibitory activity against the PQS1-mutant strain while maintaining potent activity against the PQS3-mutant strain (Figure 4E, Figure 4—figure supplement 3C). This suggests that the compound can selectively target the PQS1 of the high SHAPE-high Shannon region in cells. This approach can be extended to other RNA viruses and potentially to DNA viruses where RNA intermediates play critical roles in the viral life cycle.

Single-stranded RNA regions in viral genomes as the targets for antiviral therapy.
(A) PQSs in the single strands are easier to be induced into G4s by ligands than those in the paired regions. (B) Influence of the secondary structures on the binding efficiency of siRNAs.
Recently, the FDA has approved several siRNAs for therapeutic purposes (Ali Zaidi et al., 2023). siRNA-mediated RNAi is particularly suited to inhibit viruses with single-stranded RNA genomes due to its potential to specifically silence virtually any therapeutic target. When screening siRNA target sequences, sequence conservation and functional importance of the target region are usually given priority. However, the RNA structure of the target region also has a significant impact on the efficiency of RNAi (Figure 6B; Cao et al., 2021, Manfredonia et al., 2020, Piasecka et al., 2020, Patzel et al., 2005, Sagan et al., 2010, Szabat et al., 2020). Our approach uniquely integrates in situ RNA structural data (SHAPE reactivity and Shannon entropy) to prioritize siRNA targets within stable single-stranded regions (high SHAPE reactivity, low Shannon entropy), which are experimentally validated as accessible in infected cells. This represents a significant departure from traditional siRNA design methods that rely primarily on sequence conservation, thermodynamic rules (e.g., Tuschl rules), or in vitro structural predictions (Ali Zaidi et al., 2023, Qureshi et al., 2018, Tang and Khvorova, 2024), which may not accurately reflect intracellular RNA accessibility. Bowden-Reid et al. designed 39 antiviral siRNAs against various SARS-CoV-2 variants based on sequence conservation, ultimately identifying 8 highly effective sequences (Bowden-Reid et al., 2023). Notably, five of these effective sequences targeted regions that were located in high SHAPE-high Shannon regions according to SARS-CoV-2 SHAPE datasets (Supplementary file 8; Manfredonia et al., 2020). This independent finding aligns perfectly with our conclusions and demonstrates that SHAPE-based siRNA design outperforms sequence/structure-agnostic approaches, at least in terms of significantly improving antiviral siRNA screening efficiency. Given the growing availability of SHAPE datasets for numerous viruses, we are confident that our methodology will facilitate a more precise design of antiviral siRNAs.
Materials and methods
Cell culture and PEDV infection
Request a detailed protocolVero E6 cells were cultured in 10 cm culture plates in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% fetal bovine serum at 37°C with 5% CO2. Cells were infected with PEDV strain AJ1102 (GenBank: MK584552.1) at an MOI of 0.1 when the cell confluence reached 95%.
In situ RNA probing and purification
Request a detailed protocolDuring the logarithmic phase of viral genome replication after infection (12 h), infected cells were washed three times with 10 mL cold PBS and then resuspended in PBS containing 2-methylnicotinic acid imidazolide (NAI) (Spitale et al., 2013) at a final concentration of 100 mM. Two independent replicates were performed for the ex vivo experiments. An equal amount of DMSO was added to another sample as a control. The samples were then incubated at room temperature with continuous shaking for 15 min, after which the reaction was terminated by adding DTT at a final concentration of 0.5 M. The supernatant was then removed and the cells were lysed by adding 6 mL TRIzol reagent. RNA was extracted by adding 1.2 mL of chloroform:isoamyl alcohol (24:1). 3 mL of the upper aqueous phase was transferred to a new tube, 6 mL of anhydrous ethanol was added, and the mixture was incubated overnight at –20°C. RNA was pelleted at 18,000 × g for 20 min at 4°C, washed twice with 80% EtOH, and then the pellet was air-dried. Finally, RNA was dissolved in RNase-free water and then stored at –80°C until further use.
Isolation and fragmentation of poly(A) PEDV genomic RNA
Request a detailed protocolWe isolated poly(A) RNA with VAHTS mRNA capture beads (Vazyme #N401) according to the manufacturer’s instructions. 5 µg of total RNA was diluted to 50 µL with RNase-free water and purified in two rounds with mRNA capture beads. Poly(A) RNA was fragmented to a median size of 150–200 nt by incubation at 94°C for 8 min in fragmentation buffer (65 mM Tris–HCl; 4 mM MgCl2; 95 mM KCl; pH 8.0) and then purified with an RNA Clean & Concentrator-5 kit (Zymo Research, R1015) following the manufacturer’s instructions.
Reverse transcription of modified PEDV genomic RNAs
Request a detailed protocolAliquots of 20 μL poly(A) RNA (~500 ng) were distributed into 0.65 mL RNase-free tubes, followed by the addition of 2 μL random primer (200 ng/μL) to each tube. The samples were incubated at 65°C for 5 min. Sixteen microliters of 2.5×MaP buffer (125 mM Tris, pH 8.0; 187.5 mM KCl; 15 mM MnCl2; 25 mM DTT and 1.25 mM dNTPs) were added to each tube, and the mixture was incubated at 25°C for 2 min. Then, 2 µL SuperScript II reverse transcriptase (200 U/µL; Life Technologies, cat. no. 18064-014) was added and the mixture was incubated at 42°C for 3 h. The reactions were heat-inactivated by incubation at 75°C for 15 min, after which the mixture was placed on ice. cDNA was purified from reactions using G-25 micro spin columns (GE Healthcare, cat. no. 27-5325-01) following the manufacturer’s instructions. The second-strand cDNA synthesis reaction was then immediately performed. The volume of first-strand cDNA was adjusted to 56 µL with nuclease-free H2O, 6.5 µL NEBNext Second Strand Synthesis Reaction Buffer (10×; New England Biolabs, cat. no. E6111S) and 3.5 µL NEBNext Second Strand Synthesis Enzyme Mix (New England Biolabs, cat. no. E6111S) were added and the mixture was then incubated at 16°C for 2.5 h.
SHAPE-MaP sequencing library construction and quality control
Request a detailed protocolSequencing libraries were generated using a VAHTS Universal V8 RNA-seq Library Prep Kit for MGI (Vazyme #NRM605) according to the manufacturer’s protocol. Libraries concentrations were measured using a Qubit dsDNA HS Assay Kit (Thermo Fisher, cat. no. Q32851). 1 µL of the final library product was removed and the library quality was analyzed using an Agilent DNA1000 Kit (Agilent, cat. no. 5067-1504).
SHAPE-MaP data analysis
Request a detailed protocolSHAPE-MaP data analysis was adapted from Manfredonia et al., 2020, and performed using the RNA Framework v2.6.9 (Incarnato et al., 2018; Incarnato and Brewer, 2025), employing the following steps:
(a) Data was preprocessed and aligned to the PEDV reference genome using the rf-map module with the parameters: -ctn -cmn 0 -cqo -cq5 20 -b2. The reads with internal Ns were discarded, and low-quality bases (Phred < 20) were excluded before alignment.
(b) The PEDV genome index was built using Bowtie2 software (Langmead and Salzberg, 2012). Individual nucleotide mutation signals were assessed with the rf-count module with parameters: -m -rd, generating a binary RC file with sequence data, mutation counts, coverage, and total aligned fragments.
(c) The rf-norm module was applied for data normalization with parameters: -sm 3 -nm 3 -rb AC -n 500 -mm 1, utilizing the scoring method developed by Siegfried et al., 2014.
PEDV RNA secondary structure modeling guided by SHAPE data
Request a detailed protocolThe rf-fold module, employing a previously established windowed approach, (Siegfried et al., 2014), predicts the RNA secondary structure of the entire PEDV genome using the normalization step’s XML file, with a default pseudo-free energy intercept and slope values of –0.6 (kcal/mol) and 1.8 (kcal/mol), respectively, with the following parameters: -fw 3000 -fo 300 -wt 200 -pw 1000 -po 250 -dp -sh -nlp -md 600. This means that RNA structure will generate a 3000 nt window and slide in 300 nt increments over the genome in 1000 nt partition windows, predicting the most reasonable secondary structure of the genome under the condition that the maximum base pairing distance does not exceed 600 nt. SHAPE reactivities were incorporated in the form of soft constraints using the Vienna RNA package (Lorenz et al., 2016). For each window, the first and last 100 nt were ignored to avoid terminal biases. Additionally, RF Fold can provide pairing probabilities between bases, Shannon entropy of individual bases, and potential pseudoknot structures. The Shannon entropy is commonly used to measure the diversity of RNA secondary structures:
where P represents the pairing probability between bases (Equation 1).
Identification of characterized regions in viral genomes
Request a detailed protocolBased on previous research methods to identify well-folded regions (Siegfried et al., 2014), modifications have been made as follows:
(a) For the identification of low SHAPE-low Shannon regions, we first calculated the local medians of SHAPE reactivities and Shannon entropies in sliding, centered, 50 nt windows, and subtracted the global median for visualization. Then, a 50 nt window was slid along the genome with a 1 nt step. Windows with >75% of the bases below both the global SHAPE and the Shannon median (calculated on the full PEDV genome) were selected, and windows less than 10 nt apart were merged. Following the initial identification of well-folded regions, entire structures exhibiting consistent conformations in both datasets and having at least 50% of their bases in well-folded regions were further selected. Subsequently, the well-folded RNA structures within cells were determined.
(b) For the identification of low SHAPE-high Shannon regions, a 50 nt window was slid along the PEDV genome. Windows with >75% of the bases below the global SHAPE median and above the global Shannon median (based on the full PEDV genome) were selected, and windows less than 10 nt apart were merged.
(c) For the identification of high SHAPE-low Shannon regions, instead, a window of 25 nt was slid along the PEDV genome. Windows with >75% of the bases were above the global SHAPE median and below the global Shannon median (calculated on the full PEDV genome), and >50% of the bases were predicted to be single-stranded in the MEA structure; windows less than 10 nt apart were selected and merged.
(d) For the identification of high SHAPE-high Shannon regions, a 25 nt window was slid along the PEDV genome. Windows with >75% of the bases above both the global SHAPE and Shannon median (based on the complete PEDV genome) were selected, and windows less than 10 nt apart were merged.
G-quadruplex forming sequence prediction and conservativeness analysis
Request a detailed protocolThe G4 sequence was defined as G≥2 N1–6 G≥2 N1–6 G≥2 N1–6 G≥2, where G refers to guanine and N refers to any nucleotide, including guanine. Four independent G4 prediction tools (QGRS-mapper, G4Catchall, ImGQfinder, and pqsfinder) were used to analyse the potential G4 forming sequences (PQSs) in the PEDV genome (Doluca, 2019, Hon et al., 2017, Kikin et al., 2006, Varizhuk et al., 2017). Conservation analysis was performed using Molecular Evolutionary Genetics Analysis (MEGA) software (version 6.0; available at https://www.megasoftware.net/). 761 PEDV genome sequences were retrieved from the National Center for Biotechnology Information website (NCBI, https://www.ncbi.nlm.nih.gov/). The graphical representation of the alignment sequence was generated by WebLogo 2 (https://weblogo.berkeley.edu/) (Crooks et al., 2004).
Verification of G-quadruplex structure formation by native PAGE
Request a detailed protocolCy2-tagged RNA was annealed in buffers (100 mM potassium arsenate, pH 7.0) containing different concentrations of KCl (0, 50 mM, 100 mM). The marker used in this experiment was a mixture of P15 (UAAUACGACUCACUA), M17 (UAAUACGACUCACUAUAUA), and M20 (UAAUACGACUCACUAUACGA). Acrylamide at a 20% concentration was used to separate the gel. The concentration of the RNA sample was 100 nM and electrophoresis was performed at 100 V for 2 h. Native PAGE was carried out under a controlled temperature (4.0°C) with a vertical electrophoresis instrument (DYCZ-22A; Bio-Rad, Beijing, China). After electrophoresis, RNA oligomers were scanned using a Pharos FX molecular imager for visualization.
CD spectroscopy
Request a detailed protocolRNAs (20.0 µM) were dissolved in 10 mM Tris–HCl buffer (pH 7.0) containing varying concentrations of KCl (0, 50, 100 mM). The CD experiment was carried out at 25°C using a JASCO-810 spectrophotometer (Jasco, Easton, MD, USA) with a quartz cell path length of 1.0 mm. CD spectra were collected from 320 to 200 nm with a scanning speed of 50 nm/min. All CD spectra were baseline corrected for the signal contribution of the buffer and represent the mean of at least two repeats.
NMR spectroscopy
Request a detailed protocol1H NMR spectra were recorded at 298 K using an 800-MHz Brock Avance DRX-800 spectrometer equipped with a low-temperature triple-resonance reverse autotuning and matching probe. The presumption of water is achieved by excitation engraving. RNA samples were dissolved in DEPC solution containing 200 mM KCl and 10% D2O at a final concentration of 0.5 mM. The sample was heated at 95°C for 5.0 min and slowly cooled to 25°C. The facility was set for 100 ms diffusion time (TD), 50 ms eddy recovery time (TE), and 2s relaxation delay (TR).
CD melting studies
Request a detailed protocolBraco-19 was dissolved in DMSO at a concentration of 20 mM. RNA samples were dissolved in buffer containing 10 mM Tris–HCl (pH 7.0) and 100 mM KCl at a final concentration of 20 µM. The samples were heated at 95°C for 10 min and then slowly annealed to room temperature. CD melting was measured on a Jasco-810 spectrophotometer equipped with a water bath temperature control attachment. CD melting curves were recorded using a heating rate of 0.2°C/min and absorbance values were collected at 1°C intervals. The thermal stability of G4 RNA was measured by recording the CD value at 264 nm in relation to temperature. The required temperature range was set from 4.0°C to 98°C.
RNA stop assay
Request a detailed protocolP15 (300 nM) and template RNA (600 nM) were annealed in buffer at 25°C (50 mM HEPES, 20 mM NaCl, 1.0 mM KCl, 5 mM MgCl2, and 4 mM DTT [pH 7.0]), heated at 95°C for 5.0 min, and slowly cooled to 4°C. Subsequently, NTPs (final concentration, 200 µM) and 3Dpol (recombinant RNA-dependent RNA polymerase; Abcam, ab277617, 0.02 mg/reaction) were added to the system, and the reaction was carried out at 33°C for 20 min. An equal amount of stop buffer (95% formamide) was added, and the reaction was stopped after heating at 90°C for 4.0 min. The products were loaded and separated on a 20% denatured polyacrylamide gel. Finally, the gel was scanned using a Pharos FX molecular Imager (Bio-Rad) operated in fluorescence mode.
EGFP expression repression by RNA G-quadruplex stabilization
Request a detailed protocolWe constructed pEGFP-C1 vector with wild-type PQS1 sequences or G4-mutated PQS1mut sequences fused to the N-terminus of an EGFP reporter gene. Vero cells were transfected with Lipofectamine 3000 (Invitrogen, USA) according to the manufacturer’s protocol and incubated in media supplemented with various concentrations of Braco-19 (Gowan et al., 2002) for 48 h. The cells were then washed with PBS and fixed at 25°C for 15 min with 4% paraformaldehyde. The samples were washed three times with PBS for 5 min each and then fixed with precooled methanol at –20°C for 10 min. Observations were conducted using a two-photon laser scanning confocal microscope (Olympus, fv1000mp).
Construction of the RNA G4-mutated recombinant virus
Request a detailed protocolThe pBAC-CMV-PEDV template was obtained from Prof. Xiao Shaobo (College of Veterinary Medicine, Huazhong Agricultural University). The primers used were chemically synthesized (Supplementary file 1). The G-rich region sequence (3109–3125 nt) was mutated (AAGGTTACAGGTGGTTGG to AAAGTTACAGGTGGTTGG). Two pairs of primers (PQS1Mut-F1 and PQS1Mut-R1, PQS1Mut-F2 and PQS1Mut-R2) were used to amplify the fragment containing the mutant site. The PQS1mut fragment containing the mutation sites was amplified by overlap PCR. A scaffold oligo was used as a template, PQS1m-sgRNA-a and PQS1m-sgRNA-a were used as primers to synthesize guide DNA, and then guide RNA (sgRNA) was obtained by reverse transcription. pBAC-CMV-PEDV was cut with the Cas9 enzyme and sgRNA. The PQS1mut fragment was cloned and inserted into a linearized pBAC-CMV-PEDV vector to generate the pBAC-PEDV-PQS1mut mutant plasmid. Vero cells were placed on 6-well plates and plasmid transfection was performed when the cells grew to 60% to 80% confluence according to the manufacturer’s protocol (Lipofectamine 3000, Thermo, L3000015). If no lesions appeared, the freeze–thaw solution was added to the six-well plate with 80% confluence after three repeated freeze–thaw cycles for further culture. When lesions appeared, the successful mutation of the mutation site was confirmed by reverse transcription sequencing and the PEDV-PQS1mut mutant strain was obtained.
Determination of the antiviral activity of Braco-19 against PEDV by RT-qPCR
Request a detailed protocolBraco-19 at different concentrations (0.5, 1, 2, and 5 µM) were inoculated into monolayer Vero cells and incubated for 2 h. PEDV-WT and PEDV-PQS1mut were inoculated at an MOI (multiplicity of infection) of 0.1 and incubated for 2 h. DMEM was replaced with DMEM containing Braco-19, and the mixture was incubated for 16 h. cDNA was synthesized from 100 ng of total RNA per sample using the Transcriptor First Strand cDNA Synthesis Kit (Roche, Switzerland) according to the manufacturer’s instructions. Quantitative PCR (qPCR) was carried out by following the Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE) guidelines (Bustin et al., 2009) using Taq Pro Universal SYBR qPCR Master Mix (Vazyme, Q712) according to the manufacturer’s protocol. mRNA expression was normalized to beta-actin levels with primers beta-actin-F/R (Supplementary file 1). PEDV genomic RNA was quantified at the nonstructural protein 12 (nsp12) gene using primers qPCR-nsp12-F/R (Supplementary file 1).
Western blot analysis to detect PEDV N gene expression
Request a detailed protocolProtein samples were separated by SDS–PAGE using a 12.5% polyacrylamide gel and then transferred to a polyvinylidene difluoride (PVDF) membrane (Millipore, USA). The membrane was blocked with 5% BSA in Tween-PBS buffer and incubated overnight at 4°C with primary antibodies. After three washes, it was incubated with HRP-conjugated secondary antibodies (Beyotime, China). Following three more washes, protein bands were detected using enhanced chemiluminescence (Bio-Rad, USA).
CCK-8 assays
Request a detailed protocolVero cells were treated with graded compound concentrations for 12 h in 96-well plates, then incubated with 10 μL CCK-8 reagent (Beyotime) for 2 h at 37°C. Absorbance at 450 nm was measured using a microplate reader to assess proliferation.
Indirect immunofluorescence assay (IFA)
Request a detailed protocolCells in 24-well plates were pretreated with test compounds, infected with PEDV (MOI=0.1, 2 h, 37°C), then maintained in compound-containing DMEM. After 24 hpi, cells were fixed, permeabilized, and washed before blocking and incubating with primary (mouse anti-PEDV N protein antibody) and Alexa Fluor 488-secondary antibody. Nuclei were DAPI-stained, and fluorescence images were captured using an Olympus IX73 microscope.
Isothermal titration calorimetry (ITC)
Request a detailed protocolITC experiments were carried out on an AutoiTC100 titration calorimeter (MicroCal) at 25°C. All solutions (buffer, RNA, and ligand) were degassed prior to experimentation. A 10 µM solution of RNA 5’ UTR-SL5 was introduced into the cell and a compound solution (250 μM) was added to the rotating syringe (750 rpm). The RNA solution received 40 μL of compound through 20 sequential injections (1.5 min intervals), preceded by a 0.4 μL priming injection to compensate for syringe diffusion during system equilibration. This initial injection was excluded from the data fitting. Furthermore, dilution experiments were performed with buffer (100 mM HEPES [pH 8.0], 100 mM NaCl, and 10 mM MgCl2) in the cell, and the compound remained in the syringe. We utilized Microcal ITC analysis software to analyze the raw ITC data using the two-site binding model. The dilution data were subtracted from the raw interaction data of the compound with RNA prior to analysis. Following fitting of these experimental data, the enthalpy change during the process was determined.
siRNA antiviral assay
Request a detailed protocolsiRNAs were synthesized in GENERAL BIOL and transfected into 8×104 Vero E6 cells in 24-well plates using 10 pmol of each siRNA and Lipofectamine RNAiMAX Transfection Reagent (Invitrogen, 13778075) according to the manufacturer’s protocol. Transfected cells were infected with PEDV (MOI = 0.1) after 24 h. After 16 h of PEDV infection, minimal CPEs on Vero cells were observed under the microscope. To quantify PEDV RNA levels in infected cells, total RNAs were extracted with TRIzol Reagent (Invitrogen, 15596026). For RT-qPCR, 100 ng of total RNAs were converted into cDNA using HiScript II 1st Strand cDNA Synthesis Kit (Vazyme, R212). qPCR was performed with ChamQ Universal SYBR qPCR Master Mix (Vazyme, Q711) on a StepOnePlus real-time PCR machine (Applied Biosystems). Primers targeting the nucleocapsid (N) gene of PEDV and the beta-actin gene of Vero cells were used for qPCR. For the 50% tissue culture infectious dose (TCID50) assay, monolayers of Vero cells seeded in 96-well plates were washed twice with DMEM. The supernatants of infected cells were serially diluted tenfold. Cells were then inoculated with eight replicates of the appropriate dilutions of the virus suspension. TCID50 was calculated using the Reed-Muench method after culturing the cells for 48 h at 37°C in 5% CO2.
Statistical analysis
Request a detailed protocolThe experimental data was analyzed by Student’s t-test or one-way or two-way ANOVA using GraphPad Prism 8 software. *p<0.05 was considered to indicate statistical significance.
Data availability
SHAPE-MaP datasets have been deposited in the NCBI Gene Expression Omnibus (GEO) under accession number GSE271098. Processed structural data (including the following files) are available at Zenodo (https://zenodo.org/doi/10.5281/zenodo.13743827): XML files with normalized SHAPE reactivities (SHAPE_react_rep1/2.xml); WIG files with Shannon entropies (GSE271098_PEDV_SHAPE_Shannon_incell_rep1/2.wig.gz); CT file with the full genome secondary structure of PEDV (PEDV_incell_secondary structure.ct).
-
ZenodoExploiting functional regions in the viral RNA genome as druggable entities.https://doi.org/10.5281/zenodo.13743827
-
NCBI Gene Expression OmnibusID GSE271098. Exploiting functional regions in the viral RNA genome as druggable entities.
References
-
Engineering siRNA therapeutics: challenges and strategiesJournal of Nanobiotechnology 21:381.https://doi.org/10.1186/s12951-023-02147-z
-
Therapeutic antisense oligonucleotides are coming of ageAnnual Review of Medicine 70:307–321.https://doi.org/10.1146/annurev-med-041217-010829
-
Physical and functional analysis of viral rna genomes by SHAPEAnnual Review of Virology 6:93–117.https://doi.org/10.1146/annurev-virology-092917-043315
-
Structural basis of a small molecule targeting RNA for a specific splicing correctionNature Chemical Biology 15:1191–1198.https://doi.org/10.1038/s41589-019-0384-5
-
The architecture of the SARS-CoV-2 RNA genome inside virionNature Communications 12:3917.https://doi.org/10.1038/s41467-021-22785-x
-
Targeting RNA structures with small moleculesNature Reviews. Drug Discovery 21:736–762.https://doi.org/10.1038/s41573-022-00521-4
-
Antisense technology: an overview and prospectusNature Reviews. Drug Discovery 20:427–453.https://doi.org/10.1038/s41573-021-00162-z
-
Targeting RNA with small molecules to capture opportunities at the intersection of chemistry, biology, and medicineJournal of the American Chemical Society 141:6776–6790.https://doi.org/10.1021/jacs.8b13419
-
G4Catchall: A G-quadruplex prediction approach considering atypical featuresJournal of Theoretical Biology 463:92–98.https://doi.org/10.1016/j.jtbi.2018.12.007
-
QGRS Mapper: a web-based server for predicting G-quadruplexes in nucleotide sequencesNucleic Acids Research 34:W676–W82.https://doi.org/10.1093/nar/gkl253
-
Small molecule approaches to targeting RNANature Reviews Chemistry 8:120–135.https://doi.org/10.1038/s41570-023-00569-9
-
Structural divergence creates new functional features in alphavirus genomesNucleic Acids Research 46:3657–3670.https://doi.org/10.1093/nar/gky012
-
Fast gapped-read alignment with Bowtie 2Nature Methods 9:357–359.https://doi.org/10.1038/nmeth.1923
-
RNA-targeted antiviral immunity: more than just RNA silencingTrends in Microbiology 27:792–805.https://doi.org/10.1016/j.tim.2019.05.007
-
Genome-wide mapping of SARS-CoV-2 RNA structures identifies therapeutically-relevant elementsNucleic Acids Research 48:12436–12452.https://doi.org/10.1093/nar/gkaa1053
-
Design of siRNAs producing unstructured guide-RNAs results in improved RNA interference efficiencyNature Biotechnology 23:1440–1444.https://doi.org/10.1038/nbt1151
-
RNA Secondary structure motifs of the influenza A Virus as targets for siRNA-mediated RNA interferenceMolecular Therapy. Nucleic Acids 19:627–642.https://doi.org/10.1016/j.omtn.2019.12.018
-
A review on current status of antiviral siRNAReviews in Medical Virology 28:e1976.https://doi.org/10.1002/rmv.1976
-
Risdiplam, the first approved small molecule splicing modifier drug as a blueprint for future transformative medicinesACS Medicinal Chemistry Letters 12:874–877.https://doi.org/10.1021/acsmedchemlett.0c00659
-
G-quadruplexes and G-quadruplex ligands: targets and tools in antiviral therapyNucleic Acids Research 46:3270–3283.https://doi.org/10.1093/nar/gky187
-
First small-molecule drug targeting RNA gains momentumNature Biotechnology 39:6–8.https://doi.org/10.1038/s41587-020-00788-1
-
RNA motif discovery by SHAPE and mutational profiling (SHAPE-MaP)Nature Methods 11:959–965.https://doi.org/10.1038/nmeth.3029
-
RNA SHAPE analysis in living cellsNature Chemical Biology 9:18–20.https://doi.org/10.1038/nchembio.1131
-
Exploring the Druggability of conserved RNA regulatory elements in the SARS-CoV-2 GenomeAngewandte Chemie 60:19191–19200.https://doi.org/10.1002/anie.202103693
-
RNAi-based drug design: considerations and future directionsNature Reviews. Drug Discovery 23:341–364.https://doi.org/10.1038/s41573-024-00912-9
-
Design of small molecules targeting RNA structure from sequenceChemical Society Reviews 49:7252–7270.https://doi.org/10.1039/d0cs00455c
-
Principles for targeting RNA with drug-like small moleculesNature Reviews. Drug Discovery 17:547–558.https://doi.org/10.1038/nrd.2018.93
-
A systematic analysis of the effect of target RNA structure on RNA interferenceNucleic Acids Research 35:4322–4330.https://doi.org/10.1093/nar/gkm437
-
In vivo structure and dynamics of the SARS-CoV-2 RNA genomeNature Communications 12:5695.https://doi.org/10.1038/s41467-021-25999-1
-
The short- and long-range RNA-RNA interactome of SARS-CoV-2Molecular Cell 80:1067–1077.https://doi.org/10.1016/j.molcel.2020.11.004
Article and author information
Author details
Funding
National Key Research and Development Program of China (2021YFD1801102)
- Dengguo Wei
Fundamental Research Funds for the Central Universities (2662023PY005)
- Dengguo Wei
National Natural Science Foundation of China (32473010)
- Dengguo Wei
China Postdoctoral Science Foundation (2024M761080)
- Dehua Luo
Natural Science Foundation of Hubei Province (2025AFB288)
- Dehua Luo
HZAU-AGIS Cooperation Fund (SZYJY2022016)
- Dengguo Wei
National Natural Science Foundation of China (22077043)
- Dengguo Wei
Natural Science Foundation of Hubei Province (2021CFA061)
- Dengguo Wei
State Key Laboratory of Agricultural Microbiology (AML2023B05)
- Dengguo Wei
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was supported by the National Key R&D Plan of China (2021YFD1801102), the Fundamental Research Funds for the Central Universities (2662023PY005), the National Natural Science Foundation of China (32473010), the China Postdoctoral Science Foundation (2024M761080), the Natural Science Foundation of Hubei Province (2025AFB288), the HZAU-AGIS Cooperation Fund (SZYJY2022016 to DW), the National Natural Science Foundation of China (22077043), the Natural Science Foundation of Hubei Province (2021CFA061 to DW), and funding from the State Key Laboratory of Agricultural Microbiology (AML2023B05 to DW). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. We thank Professor Shaobo Xiao (Huazhong Agricultural University) for providing the PEDV AJ1102 strain and for his guidance in virological experiments. The high-throughput sequencing, ITC, and confocal microscopy data were acquired at the State Key Laboratory of Agricultural Microbiology Core Facility.
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
- Version of Record updated:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.103923. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2024, Luo, Zheng, Huang et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 357
- views
-
- 18
- downloads
-
- 1
- citation
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Citations by DOI
-
- 1
- citation for umbrella DOI https://doi.org/10.7554/eLife.103923