Next Article in Journal
A Novel N-Tert-Butyl Derivatives of Pseudothiohydantoin as Potential Target in Anti-Cancer Therapy
Next Article in Special Issue
In Vitro and In Silico Anti-Arboviral Activities of Dihalogenated Phenolic Derivates of L-Tyrosine
Previous Article in Journal
Seaweed Polysaccharide Based Products and Materials: An Assessment on Their Production from a Sustainability Point of View
Previous Article in Special Issue
Synergistic Effect by Combining a gp120-Binding Protein and a gp41-Binding Antibody to Inactivate HIV-1 Virions and Inhibit HIV-1 Infection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Sherif Aly El-Kafrawy
1,2,
Sayed Sartaj Sohrab
1,2,*,
Zeenat Mirza
2,3,
Ahmed M. Hassan
1,
Fatima Alsaqaf
1 and
Esam Ibraheem Azhar
1,2
1
Special Infectious Agents Unit, King Fahd Medical Research Center, King Abdulaziz University, Jeddah 21589, Saudi Arabia
2
Department of Medical Laboratory Technology, Faculty of Applied Medical Sciences, King Abdulaziz University, Jeddah 21589, Saudi Arabia
3
King Fahd Medical Research Center, King Abdulaziz University, Jeddah 21589, Saudi Arabia
*
Author to whom correspondence should be addressed.
Molecules 2021, 26(9), 2610; https://doi.org/10.3390/molecules26092610
Submission received: 18 March 2021 / Revised: 23 April 2021 / Accepted: 27 April 2021 / Published: 29 April 2021
(This article belongs to the Special Issue Antiviral Agents for RNA-Virus Infection)

Abstract

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

1. Introduction

Human coronaviruses are well known to cause both upper and lower respiratory tract infections in humans. A new human coronavirus, isolated in June 2012, was causing respiratory illness in a 60-year-old patient who died 11 days after hospitalization in Jeddah, Saudi Arabia. The etiological agent was designated as Middle East Respiratory Syndrome Coronavirus (MERS-CoV). MERS-CoV is known as the sixth coronavirus to cause respiratory illness and has now become a serious issue for the human as well as the camel population globally. The highest number of MERS-CoV cases was reported from Saudi Arabia, 2121 cases with 788 related deaths, and the case fatality rate (CFR) peaked at 37.1% [1]. This was the first highly pathogenic coronavirus to emerge after SARS-CoV. MERS-CoV-infected individuals developed symptoms like fever, shortness of breath, and even multiorgan failure in patients with comorbidities [2,3]. The large number of pilgrims performing Hajj and Umrah in Makkah, Saudi Arabia have the potential to increase the spread of MERS-CoV infections globally. The government of Saudi Arabia has taken several control measures to avoid this spread and no cases were reported to be transmitted in the Hajj seasons since the discovery of the virus. Currently, this virus has caused 940 deaths with 2589 confirmed cases and has spread to 27 countries (https://www.ecdc.europa.eu/en/middle-east-respiratory-syndrome-coronavirus-mers-cov-situation-update accessed on 10 April 2021) including a major outbreak in the Republic of Korea in 2015 [4]. Dromedary camels were identified as the animal reservoir for the transmission of MERS-CoV to the human population [4,5,6,7,8]. Additionally, the camel workers played a possible intermediary role for the spread of virus to the human population [8], but several MERS-CoV cases were reported that had no camel contact history [9].
The MERS-CoV belongs to the lineage C of Betacoronavirus (βCoVs), of the family Coronaviridae which is a large group of viruses known to infect humans and other species [10]. Currently, coronaviruses have been divided into four main groups—Alpha, Beta, Gamma, and Delta coronaviruses. The coronaviruses are enveloped, single-stranded, positive-sense RNA viruses with a genome size of 25 to 32 kilobases with high genetic diversity which favors a high rate of recombination and the emergence of new viral strains with novel characteristics and an extended host range [11]. Studies on the phylogeny of lineage C betacoronaviruses suggest the evolution of MERS-CoV in camels occurred prior to that in humans and with a possible exchange of genetic elements among ancestral viruses either in bats, or within the camel genetic ‘mixing vessel’, leading to MERS-CoV emergence [12].
Currently, there are no approved vaccines or antiviral therapy available against MERS-CoV, but many have reached advanced stages of development [13,14,15]. The technological progression is reflected upon a series of genome-wide molecular screening platforms and computational biology approaches that provide novel insights to prompt response against emerging viral diseases. The RNA interference (RNAi) using short interfering RNAs (siRNAs) and micro-RNA (miRNAs) technology plays a significant role in antiviral therapy against human virus infection [16,17]. There are several classes of RNAi mediators including siRNAs to facilitate antiviral immunity in many mammals including humans [17]. The status and future directions for RNAi-based drugs has been significantly reviewed recently [18,19]. In 2018, an siRNAs-based drug known as patisiran (Onpattro; Alnylam Pharmaceuticals) was approved by USFDA for the treatment of hereditary transthyretin amyloidosis (hATTR) with polyneuropathy [19].
siRNAs are short sequences of RNA ranging from 21 to 23 base pairs associated with a 5′ phosphate group and a 3′ hydroxyl group. Currently, therapeutic siRNAs and miRNAs are found to be some of the most promising biopharmaceuticals in commercial space as oligonucleotide-based next-generation medicines [16]. Various miRNA and siRNA-based candidate drugs are being evaluated in more than 20 clinical trials and many of them have been evaluated against different viruses including HIV [20], Flock house virus (FHV) [21], Dengue virus [22], Hepatitis C virus (HCV) [23], Influenza virus [24,25,26], Hepatitis B virus (HBV) [27], Human Papillomavirus (HPV) [28], SARS coronavirus (SARS-CoV) [29], MERS-CoV [30], and SARS-CoV-2 [31,32]. The potential siRNAs can be predicted and designed in silico, using a computational approach to provide filtration of potential siRNA candidates with high accuracy, target specificity, and reduced off-target effects.
A detailed analysis of the viral RNA structure and information on RNA interactions is vital for the virus replication cycle and to allow an understanding of the structure−function relationships and can help in the prediction of the genomic regions that need to be considered during the designing and development of potential antiviral agents against several viruses such as influenza [24]. Other structure-based approaches use the ability of small molecules to target specific RNA structures or recognize and bind to RNA targets based on their secondary or tertiary structures [24,33,34]. Recently, in silico designing of siRNAs against respiratory viruses has been reported [18,35,36] including evaluation of siRNAs against MERS-CoV. We have previously performed in silico prediction alongside a pilot study for investigating the antiviral activity for 10 of the described siRNAs in a different cell culture system [30,37].
In continuation of our previous investigation on the potential use of siRNA to inhibit MERS-CoV replication, we report here the in silico prediction, designing as well as cytotoxicity and antiviral activity of MERS-CoV-siRNAs in Huh7 cell lines. The purpose of this study is to evaluate the antiviral and cytotoxic effect of in silico designed siRNA molecules against MERS-CoV in Huh7 which are well differentiated human hepatocarcinoma cell line pre-exhibited for good in vitro growth and increased viral titer of MERS-CoV [38]. Huh 7 cell lines were selected as a proof of concept to evaluate the cytotoxicity and antiviral activity of siRNAs in a human cell line in order to show their safety and efficacy for future use in preclinical and clinical trials.

2. Results

2.1. Sequence Isolation and Multiple Sequence Analysis

The full-genome sequence of MERS-CoV was retrieved from NCBI-PubMed and multiple sequence alignment was performed using BioEdit software (V.7.02). The resulting sequence alignment of the MERS-CoV isolates showed high conservation amongst all the analyzed sequences in the ORF1ab gene. Based on the high sequence homology and role in virus replication, the ORF1ab gene was identified and used as a target for siRNA prediction and designing. The multiple sequence alignment result has been presented in Figure 1. The size of MERS-CoV genome is around 25–32 kb and Figure 1 shows the high sequence similarity in the ORF1ab gene among different strains.

2.2. Prediction and Selection of Potential siRNAs

In silico analysis, prediction, designing and scoring of potential siRNAs targeting the ORF1ab gene of the MERS-CoV was performed using online software. Potential siRNAs with no off-target matches effect on any human mRNA sequences were filtered and selected for further evaluation. During in silico prediction, many siRNAs were found to fulfill the less favorable criteria targeting a region in the ORF1ab gene of MERS-CoV. By using strong selection criteria and a stringent strategy, we have selected and filtered out a total of twenty-one functional, off-target siRNAs that were shortlisted from four hundred and sixty-two siRNAs as per the guidelines and basic rules of filtration [18,35,36,39,40]. The designed siRNAs’ length, nucleotide content and specificity and absence of off-target effects and secondary structures in the target site were taken into consideration by the used software during the design and filtering of the siRNAs as this would influence the efficiency, precision, and functionality of siRNAs, hence leading to a better silencing outcome [35].
The list of predicted siRNAs and their characteristics are presented in Table 1. The predicted siRNAs were expected to be highly specific and potent against MERS-CoV replication. Actually, 10 of the 21 siRNAs have been previously evaluated and described [30,37], and only 4 showed an antiviral effect against MERS-CoV replication. The predicted and designed 21 siRNAs were further used to evaluate the cytotoxicity and inhibition of virus replication in Huh7 cells.

2.3. Target Accessibility and Secondary Structure Prediction

The target accessibility plot and secondary structure prediction were performed for each siRNA. Based on the results obtained from the analysis (RNAxs software), it was observed that each siRNA had strong target accessibility and binding at a specific position with minimum free energy and frequency of thermodynamics in the viral genome. In our study, the predicted siRNAs showed minimum free energy (MFE–kcal/mol) ranging from −26.53 to −17.65 kcal/mol and frequency of thermodynamic ensemble ranging from 83.64 to 35.22% (Table 1). These values are in the range of optimum values for efficient accessibility and binding to the target RNA (RNAxs web server http://rna.tbi.univie.ac.at accessed on 10 April 2021).
The RNA secondary structures of the MERS-CoV-orf1ab were predicted using RNAfold. As an alternative to the MFE structure, we proposed the centroid structure. The centroid secondary structure in dot-bracket notation with a minimum free energy (MFE) of −2014.10 kcal/mol and the ensemble diversity is 1566.99. The free energy of the thermodynamic ensemble is −2558.78 kcal/mol. In order to predict using RNAfold, the genomic regions were broken into fragments consisting of the siRNA elements, and the secondary structure topology was identified from the larger genome sequence. The ensemble centroid, the centroid of the largest cluster, and the best centroid are closer in base-pair distance to the structure determined by comparative sequence analysis than is the MFE structure.
The target accessibility and the local alteration of its structure may have an inhibitory effect on target recognition and siRNA binding. Therefore, it is reasonable to consider it in the design of siRNA tools. Applying knowledge about viral RNA secondary structures as additional criteria in selecting siRNA molecules seems to be a promising approach. The binding of oligonucleotides depends on the thermodynamic properties and stability of the desired sites for antiviral agents and the secondary structure plays a crucial role in developing successful antiviral strategies. RNA structural motifs with long unpaired nucleotide tracks such as internal loops and bulges, long hairpin loops and long terminal single stranded regions are favorable targets. The target stability plot indicates the binding position of siRNA at a specific location in the viral genome and indicates the efficient binding and stability of siRNAs in the genome. The secondary structure of each siRNA and centroid secondary structure of ORF1ab region are presented in Figure 2 and Figure 3. The results of the target accessibility plot for each siRNA are presented in Figure 4 and Table 1.
We analyzed the secondary structure of siRNAs in MERS-CoV-orf1ab gene using an online tool RNAfold. The software predicted a total of 21 siRNAs in the MERS-CoV-orf1ab genome. All the siRNAs were designed and downloaded in forna format; siRNAs were potentially different with specific features. Detailed information on these siRNAs, including location and other features are provided in Table 1. The stem-loop RNA structure of all siRNAs and its conservation are presented across MERS-CoV-orf1ab genome.
The target accessibility plot of each siRNA (Figure 4) was predicted and plotted by using the RNAxs software. Each plot indicated the significant binding of each siRNA at their defined positions. The target accessibility and effective binding position of each siRNA have been indicated by using the dotted line in all the figures. The binding position of each siRNA has been shown separately and designated as siRNA-1 to siRA-21.

2.4. siRNA Transfection and Cytotoxicity Assay

The reverse transfection method was applied, and siRNAs were delivered by using Lipofectamine 2000 (ThermoFisher Scientific, Waltham, MA, USA) in grown and 60–80% confluent Huh7 cells. The cytotoxicity of siRNAs in Huh7 cells was found to be concentration dependent. None of the evaluated siRNAs showed significant cytotoxicity to Huh7 cells and the CC50 for all siRNAs were found to be >100 nM (Table 2).

2.5. Evaluation of Virus Replication Inhibition

In vitro evaluation of the inhibitory activity of siRNAs was performed through relative quantification of the viral RNA by real-time RT-PCR in both cell lysate and cell culture supernatant using positive virus control as a reference. The inhibition of MERS-CoV replication was found to be concentration dependent for 3 of the investigated 21 siRNAs (siRNA 6,16 and 19) both in the cell lysate and in the culture supernatant. The siRNAs showed IC50s of 26.22, 13.87 and 11.12nM for siRNA 6,16 and 19; respectively in supernatants and, 17.56, 16.25 and 5.17nM for siRNA 6,16 and 19; respectively, in cell lysate (Table 2). The inhibition curves for the three siRNAs are shown in Figure 5a,b.

3. Discussion

The emergence of MERS-CoV was reported for the first time in June 2012 from Jeddah, Saudi Arabia and this virus has reportedly spread to 27 countries with 2589 confirmed cases and 940 deaths globally (https://www.ecdc.europa.eu/en/middle-east-respiratory-syndrome-coronavirus-mers-cov-situation-update accessed on 10 April 2021) [1]. The highest fatality rate was reported in the Arabian Peninsula. Several research studies have been conducted and published globally, significantly contributing to the design and development of the measures for control of virus spread and disease management. Currently, there are no approved vaccine or antivirals available, but many studies have reached advanced stages of research and development focusing on antivirals and vaccines against MERS-CoV globally [13,41].
The nucleic acid-based siRNAs/miRNAs therapeutics have emerged as an alternative antiviral therapeutic option and contributed significantly to the development of antivirals against MERS-CoV [18,30,37,40]. Currently, many siRNA/miRNA-based therapeutics are being evaluated against various diseases, including those viral-borne, and RNAi technology is being applied to silence the expression of desired/undesired genes [16,42]. The first siRNA with documented effect in humans, ALNRSV01, is a 19 bp RNA duplex with two (2′-deoxy) thymidine overhangs on both 3′ ends to prevent its nuclease degradation [17]. The replication of MERS-CoV is mediated by the ORF1ab gene while the spike protein (S) gene enables the virus attachment to the host cells and is subject to immune pressure which might lead to immune escape mutations. Multiple advanced platforms and innovative strategies are being used to develop effective vaccines and antivirals against MERS-CoV. The MERS-CoV- S gene and receptor binding domain (RBD) have been used extensively for development of MERS-CoV vaccines [13,14,15,41,43,44,45]. The ORF1ab region includes two-thirds of the coronavirus genome and encodes for nonstructural proteins [39]. The benefits of using RNAi technology in therapeutics was hampered by many obstacles, such as: off-target effects, delivery, stability, and stimulation of immune responses. The continuous research effort was successful in reducing most of the obstacles to the minimum and many siRNA-based therapeutics have reached an advanced stage of research against various diseases [46,47,48]. Recently, a few siRNAs against HCV and MERS-CoV have been evaluated and found to be very efficient to inhibit the viral replication [30,49]. Only two studies have investigated the potential use of siRNAs against MERS-CoV but they only reported the in silico design using online software but none of them have evaluated the predicted siRNA for their cytotoxicity or antiviral activity [39,40]. The in silico prediction and designing of siRNAs, provided an opportunity to filter, screen and select the potential siRNAs against MERS-CoV. The effective use of online software also helps to minimize the siRNAs with undesired properties. The use of ORF1ab gene as a potential target for the design of siRNAs has been selected and published by other researchers with no in vitro validation to investigate the efficacy and potential of the designed siRNAs [39,40].
In this study, we have predicted, designed, selected, screened, and validated the cytotoxicity and antiviral efficacy of potential siRNAs in the Huh 7 cell line. Using the computational approach and bioinformatic analysis, this resulted in several siRNAs against MERS-CoV ORF1ab gene. We have designed and evaluated the secondary structure, target accessibility plot and thermodynamic properties of all 21 siRNAs and centroid secondary structure of ORF1ab gene using online software. We have evaluated all the twenty-one siRNAs because of their better target accessibility. The predicted siRNAs were chemically synthesized and evaluated by using Lipofectamine 2000 mediated delivery in Huh 7 cells. The lower binding energy of the siRNAs indicates better interaction therefore a better chance of target inhibition. The free binding energy of our predicted siRNAs ranged from −26.53 to −17.65 kcal/mol (Table 1). This low binding energy gives a better chance for the siRNA to bind to its target RNA for increased potential of antiviral activity.
The secondary structure of siRNAs provided useful information about the locations in the viral genome and their minimum free energy (MFE) with free energy of the thermodynamic ensemble (TE). It has been observed that the MFE and TE of each siRNA varied significantly, and these variations are significantly responsible for the effective binding of each siRNA in the viral genome and their potential effects on the degradation of the viral genome with the induction of RSIC and RNAi process. Reportedly, the secondary structure of siRNA provides vital information about the features for efficient binding and significantly regulating the translation and replication processes of the viral genome [24].
In our analysis, we have designed and predicted the secondary structure of 21 siRNAs. All the siRNAs formed the secondary structure and showed the variable binding free energy and thermodynamic properties. The secondary structure of a given viral RNA can be predicted and designed rationally using multiple techniques. Many bioinformatics tools for prediction and designing have been designed, developed and are being successfully used to identify new biological insights pertaining to RNA structures [50]. siRNAs can silence an unknown number of unintended genes due to the promiscuous entry of siRNAs into endogenous miRNA machinery and the recognition of complementarity to seed regions [51]. The bioinformatics-based methods become indispensable tools for designing siRNA-based therapeutics with high efficacy and minimal off-target effects. The in silico prediction of the RNA secondary structure helps to understand the structure–function relationships of target RNA to design and identify novel bioactive therapeutic compounds, effective diagnostics, and successful antisense strategy, especially against undruggable targets. Online tools use different algorithms, search databases of different species and provide comprehensive filtering, and evaluate seed duplex Tm for reduced off-target effects. The thermodynamics approaches apply multiple algorithms, such as minimizing free energy (MFE), high and defined accuracy, and sampling-based models [52]. The binding of siRNA efficacy strongly depends on thermodynamic properties and accessibility of the 3′-end of their binding sites and significantly reflects the underlying mechanisms of siRNA function for therapeutic agents [24,53]. For the MFE structure, the ensemble centroid (reflective of the high-frequency base pairs in the structure sample) and the centroid of the largest cluster, the results are comparable with marginal overall improvements by the centroids. Although the best centroids are the best predictors, these centroids cannot be defined when a reference structure is unavailable. However, it is an appealing feature that the best centroid predictions are based on only three to four clusters, on average. The centroid for a given set of structures is the structure in the entire structure ensemble that has the minimum total base-pair distance to the structures in the set. Thus, the centroid structure can be considered as the single structure that best represents the central tendency of the set. A centroid is referred to as the ensemble centroid when the set is the entire collection of structures sampled from the ensemble [54]. Despite the potential of siRNA-based therapeutics, challenges like rapid degradation, poor cellular uptake and off-target effects remain.
The centroid secondary structure of MERS-CoV-orf1ab genome was designed by online software and showed the binding position of siRNAs at different nucleotide sequences. The centroid secondary structure showed the binding of high frequency base pairs with the siRNAs. Target accessibility plays an important role in the efficiency of an siRNA. Currently, a few programs are available for computing the accessibility using the energy model. Sfold computes the accessibilities using posterior sampling techniques. Recently, a few reports have been published discussing the importance of accessibility around the target regions of functional RNAs. By using the computational predictions and experimental validation it has been demonstrated that the secondary structure of RNA and the binding site of siRNA and their effective accessibility around the target sites influence the siRNAs’ efficacy [53,55]. Additionally, the silencing activity of structurally similar duplexes with different sequences varies significantly. The thermodynamic properties of siRNA play a key role on the stages of RISC activation and mRNA target cleavage, determining the efficiency of strand dissociation, strand selection and mRNA cleavage. Any changes in the nucleotide affecting the thermodynamic stability of the duplex could change the silencing activity of siRNA [56]. Based on experimental data it has been suggested that highly active siRNAs are likely to have lower internal stabilities than less active siRNAs [57]. An siRNA or shRNA targeted to an accessible region will not necessarily be functional if the guide strand cannot successfully assemble into the RISC. The asymmetry of siRNA duplex ends is important for RISC assembly, whereas target accessibility is important for the downstream step of target recognition in the RNAi pathway. An siRNA designed for an accessible target site will not necessarily be functional, if it does not have the favorable differential stability of siRNA duplex ends for effective assembly of the guide strand into RISC [58].
In our study, the target accessibility plot of all 21 siRNAs generated from online software showed the binding position of siRNAs in the MERS-CoV-orf1ab genome. The efficient binding of siRNAs in the viral genome reflects the potential activity in the viral genome. As it is well reported that the efficient binding of siRNAs plays an important role in the RNAi pathway, the information about the secondary structure and target accessibility of siRNAs is very important for their efficient binding to the target. The asymmetry of the siRNA duplex plays a significant role to facilitate the assembly of the RNA-induced silencing complex (RISC). The estimation of target accessibility and duplex asymmetry can improve the target knockdown level significantly by nearly 40 and 26%, respectively. The duplex asymmetry has a significant upstream effect on RISC assembly and target accessibility has a strong downstream effect on target recognition [58]. The secondary structure of viral RNA is an important regulator of virus biology because it plays an important role in virion ribonucleoproteins and virus assembly and release of infectious virus particles. The terminal (5′ and 3′) ends of viral RNA fragments are highly conserved and known as virus promotors for virus replication and termination [24]. These regions form a base pairing resulting in a double stranded structure (Panhandle Structure, Figure 2) which is finally recognized and bound by viral polymerases. The fork or corkscrew model has been proposed to explain the initiation of viral RNA transcription.
The prediction of RNA secondary structure by using bioinformatics tools and in combination with experimental data, enables the design and development of novel, efficient and highly bioactive siRNA molecules with high potency. The RNA structure-based target-site selecting method provides novel prospects for other antiviral therapeutic agents, such as small molecules. To control and silence gene expression, synthetic oligonucleotides, including antisense oligonucleotides, small interfering RNAs (siRNAs) and microRNAs (miRNAs) as well as catalytic nucleic acids, have been broadly investigated and their mechanism of action is based on selective binding to complementary sequences of target viral RNAs [59].
The cytotoxicity of each siRNA was evaluated, and the investigated siRNAs were found to have minimum cytotoxicity on the Huh7 cells with CC50s > 100 nM. Three of the investigated siRNAs were found to have concentration-dependent antiviral activities in both supernatants and cell lysates as shown in Figure 5a,b. Based on their IC50, three of the designed siRNAs exhibited better concentration-dependent viral inhibition as compared to others in both cell lysate and supernatant. The results obtained from this work encourage us to evaluate the combination of siRNAs in other multiple cell lines by using multiple siRNAs transfection and delivery methods in future studies. This will provide better and clearer understanding of the effectiveness and potential use of siRNA-based antiviral therapeutics against MERS-CoV.
In previous pilot studies [30,37], we reported the experimental evaluation of 10 of the designed siRNAs in HEK 293 and Vero cells and found that siRNAs 1 and 4 were the most effective in Vero cells while siRNAs 1, 2, 4, 6 and 9 were most effective in HEK 293 cells. In this study, siRNAs 2, 6, 16 and 19 were found to be the most effective indicating that the cell type affects the inhibitory activity of the siRNA. This change in activity with the cell type needs to be taken into account when performing the in vivo evaluation of the siRNAs. Dromedary camels are known to be the animal reservoir for MERS-CoV transmission to human. A search for a camel cell line is needed to evaluate the effect of the designed siRNAs in camels in order to further perform an in vivo evaluation of the candidate siRNAs.

4. Materials and Methods

4.1. Sequence Selection and Multiple Sequence Analysis

To perform multiple sequence analysis and identification of targets for siRNA design, filtration and selection, total thirty-four full-genome sequences of MERS-CoV from human and camel were retrieved from NCBI-GenBank. The multiple sequence alignment was performed to select the highly conserved region in the viral genome using BioEdit software (Version 7.2). The accession numbers of each of the isolates with the aligned homologous region have been presented in Figure 1. Due to lack of space and the large size of MERS-CoV genome, Figure 1 is showing only the start and end regions of the genome.

4.2. In Silico Prediction, Selection and Synthesis of siRNAs

The full-genome of MERS-CoV was retrieved from NCBI-GenBank and used for prediction, filtration and designing of siRNAs. The ORF1ab region was selected as the target gene for siRNAs prediction and designing. Online freely available softwares were used as an integrated bioinformatics approach to predict, design, and for the final selection of siRNAs [60,61,62,63,64,65,66]. Based on the strict selection and filtration criteria, with different scoring tools and i-SCORE, a total of 21 different siRNAs were finally selected for further study [18]. Smaller siRNAs are preferable to use for mammalian cells as longer siRNAs can induce mammalian immune response. The RNA sequences and the locations of each siRNA in the viral genome with minimum free energy are presented in Table 1. The predicted and selected siRNAs were used for custom synthesis (Integrated DNA Technologies (IDT-USA).

4.3. Secondary Structure and Target Accessibility Prediction

The secondary structure and target accessibility plot of each siRNA was predicted by using online software known as RNAxs tool (http://rna.tbi.univie.ac.at/cgi-bin/RNAxs/RNAxs.cgi accessed on 10 April 2021). The sequences of siRNAs were separately used for secondary structure prediction. The secondary structure of each siRNA is presented in Figure 2 and centroid secondary structure of ORF1ab gene in Figure 3. The target accessibility plot of each siRNA was performed by (http://rna.tbi.univie.ac.at, http://rna.tbi.univie.ac.at/cgi-bin/RNAWebSuite/RNAfold.cgi accessed on 10 April 2021) and is presented in Figure 4. The secondary structure of the siRNAs was visualized in forna format by using online software (http://rna.tbi.univie.ac.at/forna/documentation.html accessed on 10 April 2021). The color of the secondary structure indicated green: stems (canonical helices), red: multiloops (junctions), yellow: interior loops, blue: hairpin loops.

4.4. Transfection of siRNA to Huh7 Cells

The chemically synthesized siRNAs were used for Huh7 cells transfection. The cells were purchased from Creative Bioarray (Shirley, NY, USA) and grown in standard DMEM medium at 37 °C in 96-well plates with 60–80% confluency (1 × 104). A reverse transfection method was used by applying the Lipofectamine 2000 (ThermoFisher Scientific, Waltham, MA, USA) as transfection reagent following the manufacturer’s instructions. All the experiments were performed in triplicate. A defined serial dilution (0.1 to 50 nM) of each siRNA was prepared from 50 μM stocks with the addition of 100 μL Opti-MEM and Lipofectamine 2000 following the 30 min incubation at room temperature. The lipid complex with siRNAs was added gently to the grown Huh7 cells and incubated at 37 °C for 24 h.

4.5. Cytotoxicity Assay

The cytotoxicity of each siRNA in Huh7 cells was evaluated by using the MTT ([3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide]) assay kit (Invitrogen) following the manufacturer’s instructions. After transfection of cells, the media containing transfection reagents were removed and 100 μL of fresh DMEM media was added to each well and incubated for 24 h at 37 °C and 5% CO2. After incubation, 10 μL of MTT (12 mM) was added to each well followed by further incubation at 37 °C for 4 h. The precipitated formazan crystals were dissolved by mixing the 10% SDS-HCL solution (100 μL/well). The cells were further incubated for 4 h at 37 °C. The absorbance was measured at 570 nm using SpectraMax i3x imaging cytometer and the mean OD value was used for calculation of cytotoxicity of each siRNA applying the standard formula.

4.6. Analysis of MERS-CoV Replication Inhibition in Huh7 Cells

To validate the dose dependent efficacy of the synthesized siRNAs, we performed the experimental evaluation of siRNAs transfection at variable doses (50, 25, 10, 5.0, 1.0, 0.5, 0.25, and 0.1 nm) as previously described [30,37] except Huh7 cells were used for virus inoculation and real-time RT-PCR for MERS-CoV in the cell lysate and cell culture supernatant.

4.7. Inoculation of MERS-CoV to Huh7 Cells

The inoculation of MERS-CoV to the siRNA transfected Huh7 cells was performed following previously published protocols [30,37]. The transfected cells were virus inoculated and cells were further grown at 37 °C in 5% CO2 atmosphere and cytopathic effect (CPE) in Huh7 cells was observed daily for 3 days or until full CPE was developed in the positive virus control cells. All the experiments were performed in triplicate, the triplicate wells of each concentration were pooled together to provide enough volume for subsequent RNA extraction. Each experiment included positive virus control wells with no siRNA added and negative control wells (with no virus added). The cell supernatant and lysate were isolated and after adding lysis buffer, they were further used for viral RNA isolation by using QIAamp Viral RNA Mini Kit (Qiagen, Hilden, Germany) as per manufacturer’s instructions.

4.8. Real-Time PCR Assay

The purified viral RNA was used to perform quantitative real-time RT-PCR for evaluation of virus replication inhibition using MERS-CoV specific primers and probes (ORF1a and ORF1b) [5]. The Ct values of real-time PCR generated variable and comparable slopes of inhibition for each siRNA at various concentrations. each experiment included the positive and negative controls of the antiviral assays and an external positive and negative control. The variation of the Ct values for positive controls were found to have less than 10% variability between runs.

5. Conclusions

Based on the bioinformatics analysis, prediction and design of siRNAs, the secondary structure of siRNAs, centroid secondary structure of MERS-CoV-orf1ab gene, target accessibility plot of each siRNA and low cytotoxicity of all 21 siRNAs and the concentration dependent and low IC50 of the three effective siRNAs in this study, the future investigation will involve the in vivo evaluation of the selected siRNAs in a preclinical in vivo trial in small animal models. The use of the computational and bioinformatics approach for siRNA design, analysis, and real-time PCR and in vitro evaluation of siRNAs activity against MERS-CoV can be effectively used as an alternative approach for screening oligonucleotide-based antiviral therapeutics. By using this technique, a novel siRNA targeting other MERS-CoV genetic regions can be predicted, designed and evaluated as oligonucleotide-based antiviral therapeutics against MERS-CoV as well as other emerging viruses, including the culprit responsible for the currently prevalent pandemic, COVID-19.

Author Contributions

S.S.S., S.A.E.-K. conceptualized, designed, and executed the experiments, Z.M., S.S.S. performed bioinformatics study. S.S.S., S.A.E.-K., Z.M. wrote the manuscript. E.I.A.: contributed to the design and execution of experiments and critically reviewed the manuscript. S.S.S., A.M.H., F.A. performed the experiments. All authors provided critical feedback and analysis of manuscript. All authors reviewed and approved the manuscript.

Funding

This study was financially supported by King Abdulaziz City for Science and Technology (KACST), Riyadh, Saudi Arabia by providing the special grant on MERS-CoV (Project number 39-2). Therefore, all the authors are grateful for financial support from the KACST, Riyadh, Saudi Arabia.

Data Availability Statement

Data is contained within the article.

Acknowledgments

The authors would also like to gratefully acknowledge the research facility provided by Special Infectious Agents Unit (SIAU), King Fahd Medical Research Center (KFMRC), King Abdulaziz University, Jeddah, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. WHO. Middle East Respiratory Syndrome Coronavirus (MERS-CoV). Available online: http://www.emro.who.int/pandemic-epidemic-diseases/mers-cov/mers-situation-update-january-2020.html (accessed on 25 February 2021).
  2. Assiri, A.; Al-Tawfiq, J.A.; Al-Rabeeah, A.A.; Al-Rabiah, F.A.; Al-Hajjar, S.; Al-Barrak, A.; Flemban, H.; Al-Nassir, W.N.; Balkhy, H.H.; Al-Hakeem, R.F. Epidemiological, demographic, and clinical characteristics of 47 cases of Middle East respira-tory syndrome coronavirus disease from Saudi Arabia: A descriptive study. Lancet Infect. Dis. 2013, 13, 752–761. [Google Scholar] [CrossRef] [Green Version]
  3. Yin, Y.; Wunderink, R.G. MERS, SARS and other coronaviruses as causes of pneumonia. Respirology 2017, 23, 130–137. [Google Scholar] [CrossRef] [Green Version]
  4. Lee, S.S.; Wong, N.S. Probable transmission chains of Middle East respiratory syndrome coronavirus and the multiple generations of secondary infection in South Korea. Int. J. Infect. Dis. 2015, 38, 65–67. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Azhar, E.I.; El-Kafrawy, S.A.; Farraj, S.A.; Hassan, A.M.; Al-Saeed, M.S.; Hashem, A.M.; Madani, T.A. Evidence for cam-el-to-human transmission of MERS coronavirus. N. Engl. J. Med. 2014, 370, 2499–2505. [Google Scholar] [CrossRef] [PubMed]
  6. Oboho, I.K.; Tomczyk, S.M.; Al-Asmari, A.M.; Banjar, A.A.; Al-Mugti, H.; Aloraini, M.S.; Alkhaldi, K.Z.; Almohammadi, E.L.; Alraddadi, B.M.; Gerber, S.I.; et al. 2014 MERS-CoV outbreak in Jeddah--a link to health care facilities. N. Engl. J. Med. 2015, 372, 846–854. [Google Scholar] [CrossRef]
  7. Memish, Z.A.; Cotten, M.; Watson, S.J.; Kellam, P.; Zumla, A.; Alhakeem, R.F.; Assiri, A.; Al Rabeeah, A.A.; Al-Tawfiq, J.A. Community Case Clusters of Middle East Respiratory Syndrome Coronavirus in Hafr Al-Batin, Kingdom of Saudi Arabia: A Descriptive Genomic study. Int. J. Infect. Dis. 2014, 23, 63–68. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Kiyong’a, A.N.; Cook, E.A.; Okba, N.; Kivali, V.; Reusken, C.; Haagmans, B.L.; Fèvre, E.M. Middle East respiratory syndrome coronavirus (MERS-CoV) seropositive camel handlers in Kenya. Viruses 2020, 12, 396. [Google Scholar] [CrossRef] [Green Version]
  9. Al Awaidy, S.T.; Khamis, F. Middle east respiratory syndrome coronavirus (Mers-cov) in Oman: Current situation and going forward. Oman Med. J. 2019, 34, 181. [Google Scholar] [CrossRef] [PubMed]
  10. Chafekar, A.; Fielding, B.C. MERS-CoV: Understanding the latest human coronavirus threat. Viruses 2018, 10, 93. [Google Scholar] [CrossRef]
  11. Al-Omari, A.; Rabaan, A.A.; Salih, S.; Al-Tawfiq, J.A.; Memish, Z.A. MERS coronavirus outbreak: Implications for emerging viral infections. Diagn. Microbiol. Infect. Dis. 2019, 93, 265–285. [Google Scholar] [CrossRef]
  12. Corman, V.M.; Ithete, N.L.; Richards, L.R.; Schoeman, M.C.; Preiser, W.; Drosten, C.; Drexler, J.F.; Raposo, R.A.S.; Abdel-Mohsen, M.; Deng, X.; et al. Rooting the Phylogenetic Tree of Middle East Respiratory Syndrome Coronavirus by Characterization of a Conspecific Virus from an African Bat. J. Virol. 2014, 88, 11297–11303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Mubarak, A.; Alturaiki, W.; Hemida, M.G. Middle East Respiratory Syndrome Coronavirus (MERS-CoV): Infection, Immunological Response, and Vaccine Development. J. Immunol. Res. 2019, 2019, 1–11. [Google Scholar] [CrossRef] [PubMed]
  14. Yong, C.Y.; Ong, H.K.; Yeap, S.K.; Ho, K.L.; Tan, W.S. Recent Advances in the Vaccine Development Against Middle East Respiratory Syndrome-Coronavirus. Front. Microbiol. 2019, 10, 1781. [Google Scholar] [CrossRef] [Green Version]
  15. Li, Y.-D.; Chi, W.-Y.; Su, J.-H.; Ferrall, L.; Hung, C.-F.; Wu, T.-C. Coronavirus vaccine development: From SARS and MERS to COVID-19. J. Biomed. Sci. 2020, 27, 1–23. [Google Scholar] [CrossRef]
  16. Chakraborty, C.; Sharma, A.R.; Sharma, G.; Doss, C.G.P.; Lee, S.-S. Therapeutic miRNA and siRNA: Moving from Bench to Clinic as Next Generation Medicine. Mol. Ther. Nucleic Acids 2017, 8, 132–143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Levanova, A.; Poranen, M.M. RNA Interference as a Prospective Tool for the Control of Human Viral Infections. Front. Microbiol. 2018, 9, 2151. [Google Scholar] [CrossRef] [PubMed]
  18. Sohrab, S.S.; El-Kafrawy, S.A.; Mirza, Z.; Kamal, M.A.; Azhar, E.I. Design and Delivery of Therapeutic siRNAs: Application to MERS-Coronavirus. Curr. Pharm. Des. 2018, 24, 62–77. [Google Scholar] [CrossRef] [PubMed]
  19. Setten, R.L.; Rossi, J.J.; Han, S.-P. The current state and future directions of RNAi-based therapeutics. Nat. Rev. Drug Discov. 2019, 18, 421–446. [Google Scholar] [CrossRef]
  20. Liu, C.; Liang, Z.; Kong, X. Efficacy Analysis of Combinatorial siRNAs against HIV Derived from One Double Hairpin RNA Precursor. Front. Microbiol. 2017, 8, 1651. [Google Scholar] [CrossRef]
  21. Taning, C.N.T.; Christiaens, O.; Li, X.; Swevers, L.; Casteels, H.; Maes, M.; Smagghe, G. Engineered Flock House Virus for Targeted Gene Suppression Through RNAi in Fruit Flies (Drosophila melanogaster) in Vitro and in Vivo. Front. Physiol. 2018, 9, 805. [Google Scholar] [CrossRef] [Green Version]
  22. Idrees, S.; A Ashfaq, U. RNAi: Antiviral therapy against dengue virus. Asian Pac. J. Trop. Biomed. 2013, 3, 232–236. [Google Scholar] [CrossRef] [Green Version]
  23. Shahid, I.; AlRabia, M.W.; Mukhtar, M.H.; Almalki, S.S.R.; Alkahtani, S.A.; Ashgar, S.S.; Faidah, H.S.; Hafeez, M.H.; Almalki, W.H. In vitro inhibitory analysis of consensus siRNAs against NS3 gene of hepatitis C virus 1a genotype. Asian Pac. J. Trop. Med. 2017, 10, 701–709. [Google Scholar] [CrossRef] [PubMed]
  24. Szabat, M.; Lorent, D.; Czapik, T.; Tomaszewska, M.; Kierzek, E.; Kierzek, R. RNA Secondary Structure as a First Step for Rational Design of the Oligonucleotides towards Inhibition of Influenza A Virus Replication. Pathogens 2020, 9, 925. [Google Scholar] [CrossRef]
  25. Huang, D.T.-N.; Lu, C.-Y.; Shao, P.-L.; Chang, L.-Y.; Wang, J.-Y.; Chang, Y.-H.; Lai, M.-J.; Chi, Y.-H.; Huang, L.-M. In vivo inhibition of influenza A virus replication by RNA interference targeting the PB2 subunit via intratracheal delivery. PLoS ONE 2017, 12, e0174523. [Google Scholar] [CrossRef]
  26. Tsai, K.; Courtney, D.G.; Kennedy, E.M.; Cullen, B.R. Influenza A virus-derived siRNAs increase in the absence of NS1 yet fail to inhibit virus replication. RNA 2018, 24, 1172–1182. [Google Scholar] [CrossRef] [PubMed]
  27. Wang, Y.-J.; Yang, L.; Zuo, J.-P. Recent developments in antivirals against hepatitis B virus. Virus Res. 2016, 213, 205–213. [Google Scholar] [CrossRef]
  28. Zeng, J.; Yang, S.; Wang, X.; Gao, Y.; Zhang, M. Effects of siRNA-mediated suppression of HPV-11 L1 expression on the proliferation and apoptosis of vaginal epithelial cells. Exp. Ther. Med. 2017, 13, 1561–1565. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Kumar, V.; Jung, Y.-S.; Liang, P.-H. Anti-SARS coronavirus agents: A patent review (2008–present). Expert Opin. Ther. Patents 2013, 23, 1337–1348. [Google Scholar] [CrossRef] [Green Version]
  30. Sohrab, S.S.; El-Kafrawy, S.A.; Mirza, Z.; Hassan, A.M.; Alsaqaf, F.; Azhar, E.I. Designing and evaluation of MERS-CoV siRNAs in HEK-293 Cell line. J. Infect. Public Health 2020, 4, 238–243. [Google Scholar] [CrossRef]
  31. Chowdhury, U.F.; Shohan, M.U.S.; Hoque, K.I.; Beg, M.A.; Siam, M.K.S.; Moni, M.A. A computational approach to design potential siRNA molecules as a prospective tool for silencing nucleocapsid phosphoprotein and surface glycoprotein gene of SARS-CoV-2. Genomics 2020, 113, 331–343. [Google Scholar] [CrossRef]
  32. Sohrab, S.S.; El-Kafrawy, S.A.; Abbas, A.T.; Bajrai, L.H.; Azhar, E.I. In silico Prediction and Designing of Potential siRNAs to be Used as Antivirals Against SARS-CoV-2. Curr. Pharm. Des. 2021, 27, 1–12. [Google Scholar] [CrossRef]
  33. McKimm-Breschkin, J.L.; Jiang, S.; Hui, D.S.; Beigel, J.H.; Govorkova, E.A.; Lee, N. Prevention and treatment of respiratory viral infections: Presentations on antivirals, traditional therapies and host-directed interventions at the 5th ISIRV Antiviral Group conference. Antivir. Res. 2018, 149, 118–142. [Google Scholar] [CrossRef] [PubMed]
  34. Roberts, T.C.; Langer, R.; Wood, M.J.A. Advances in oligonucleotide drug delivery. Nat. Rev. Drug Discov. 2020, 19, 673–694. [Google Scholar] [CrossRef] [PubMed]
  35. Fakhr, E.; Zare, F.; Teimoori-Toolabi, L. Precise and efficient siRNA design: A key point in competent gene silencing. Cancer Gene Ther. 2016, 23, 73–82. [Google Scholar] [CrossRef] [PubMed]
  36. Zhang, D.; Lu, J. In Silico Design of siRNAs Targeting Existing and Future Respiratory Viruses with VirusSi. bioRxiv 2020. [Google Scholar] [CrossRef]
  37. Sohrab, S.S.; Aly El-Kafrawy, S.; Mirza, Z.; Hassan, A.M.; Alsaqaf, F.; Azhar, E.I. In silico prediction and experimental vali-dation of siRNAs targeting ORF1ab of MERS-CoV in Vero cell line. Saudi J. Biol. Sci. 2020, 28, 1348–1355. [Google Scholar] [CrossRef]
  38. Eckerle, I.; Corman, V.M.; Müller, M.A.; Lenk, M.; Ulrich, R.G.; Drosten, C. Replicative Capacity of MERS Coronavirus in Livestock Cell Lines. Emerg. Infect. Dis. 2014, 20, 276. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Hasan, M.M.; Akter, R.; Ullah, M.; Abedin, M.; Ullah, G.M.; Hossain, M. A computational approach for predicting role of human microRNAs in MERS-CoV genome. Adv. Bioinform. 2014, 2014. [Google Scholar] [CrossRef] [Green Version]
  40. Nur, S.M.; Hasan, M.A.; Amin, M.A.; Hossain, M.; Sharmin, T. Design of Potential RNAi (miRNA and siRNA) Molecules for Middle East Respiratory Syndrome Coronavirus (MERS-CoV) Gene Silencing by Computational Method. Interdiscip. Sci. Comput. Life Sci. 2015, 7, 257–265. [Google Scholar] [CrossRef]
  41. Hashem, A.M.; Algaissi, A.; Agrawal, A.S.; Al-Amri, S.S.; Alhabbab, R.Y.; Sohrab, S.S.; Almasoud, A.S.; Alharbi, N.K.; Peng, B.-H.; Russell, M.; et al. A Highly Immunogenic, Protective, and Safe Adenovirus-Based Vaccine Expressing Middle East Respiratory Syndrome Coronavirus S1-CD40L Fusion Protein in a Transgenic Human Dipeptidyl Peptidase 4 Mouse Model. J. Infect. Dis. 2019, 220, 1558–1567. [Google Scholar] [CrossRef]
  42. Jeang, K.-T. RNAi in the regulation of mammalian viral infections. BMC Biol. 2012, 10, 1–6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Zhou, Y.; Yang, Y.; Huang, J.; Jiang, S.; Du, L. Advances in MERS-CoV vaccines and therapeutics based on the recep-tor-binding domain. Viruses 2019, 11, 60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Xu, J.; Jia, W.; Wang, P.; Zhang, S.; Shi, X.; Wang, X.; Zhang, L. Antibodies and vaccines against Middle East respiratory syndrome coronavirus. Emerg. Microbes Infect. 2019, 8, 841–856. [Google Scholar] [CrossRef] [Green Version]
  45. Folegatti, P.M.; Bittaye, M.; Flaxman, A.; Lopez, F.R.; Bellamy, D.; Kupke, A.; Mair, C.; Makinson, R.; Sheridan, J.; Rohde, C.; et al. Safety and immunogenicity of a candidate Middle East respiratory syndrome coronavirus viral-vectored vaccine: A dose-escalation, open-label, non-randomised, uncontrolled, phase 1 trial. Lancet Infect. Dis. 2020, 20, 816–826. [Google Scholar] [CrossRef]
  46. Stein, D.A.; Perry, S.T.; Buck, M.D.; Oehmen, C.S.; Fischer, M.A.; Poore, E.; Smith, J.L.; Lancaster, A.M.; Hirsch, A.J.; Slifka, M.K.; et al. Inhibition of Dengue Virus Infections in Cell Cultures and in AG129 Mice by a Small Interfering RNA Targeting a Highly Conserved Sequence. J. Virol. 2011, 85, 10154–10166. [Google Scholar] [CrossRef] [Green Version]
  47. Carneiro, B.; Braga, A.C.S.; Batista, M.N.; Harris, M.; Rahal, P. Evaluation of Canonical siRNA and Dicer Substrate RNA for Inhibition of Hepatitis C Virus Genome Replication—A Comparative Study. PLoS ONE 2015, 10, e0117742. [Google Scholar] [CrossRef] [Green Version]
  48. Moon, J.-S.; Lee, S.-H.; Kim, E.-J.; Cho, H.; Lee, W.; Kim, G.-W.; Park, H.-J.; Cho, S.-W.; Lee, C.; Oh, J.-W. Inhibition of Hepatitis C Virus in Mice by a Small Interfering RNA Targeting a Highly Conserved Sequence in Viral IRES Pseudoknot. PLoS ONE 2016, 11, e0146710. [Google Scholar] [CrossRef] [Green Version]
  49. Elhefnawi, M.; Kim, T.; Kamar, M.A.; Min, S.; Hassan, N.M.; El-Ahwany, E.; Kim, H.; Zada, S.; Amer, M.; Windisch, M.P. In Silico Design and Experimental Validation of siRNAs Targeting Conserved Regions of Multiple Hepatitis C Virus Genotypes. PLoS ONE 2016, 11, e0159211. [Google Scholar] [CrossRef]
  50. Ahmed, F.; Sharma, M.; Al-Ghamdi, A.A.; Al-Yami, S.M.; Al-Salami, A.M.; Refai, M.Y.; Warsi, M.K.; Howladar, S.M.; Baeshen, M.N. A Comprehensive Analysis of cis-Acting RNA Elements in the SARS-CoV-2 Genome by a Bioinformatics Approach. Front. Genet. 2020, 11, 572702. [Google Scholar] [CrossRef]
  51. Doench, J.G.; Petersen, C.P.; Sharp, P.A. siRNAs can function as miRNAs. Genes Dev. 2003, 17, 438–442. [Google Scholar] [CrossRef] [Green Version]
  52. Wu, Y.; Shi, B.; Ding, X.; Liu, T.; Hu, X.; Yip, K.Y.; Yang, Z.R.; Mathews, D.H.; Lu, Z.J. Improved prediction of RNA secondary structure by integrating the free energy model with restraints derived from experimental probing data. Nucleic Acids Res. 2015, 43, 7247–7259. [Google Scholar] [CrossRef] [Green Version]
  53. Kiryu, H.; Terai, G.; Imamura, O.; Yoneyama, H.; Suzuki, K.; Asai, K. A detailed investigation of accessibilities around target sites of siRNAs and miRNAs. Bioinformatics 2011, 27, 1788–1797. [Google Scholar] [CrossRef]
  54. Ding, Y.; Chan, C.Y.; Lawrence, C.E. RNA secondary structure prediction by centroids in a Boltzmann weighted ensemble. RNA 2005, 11, 1157–1166. [Google Scholar] [CrossRef] [Green Version]
  55. Lu, Z.J.; Mathews, D.H. Efficient siRNA selection using hybridization thermodynamics. Nucleic Acids Res. 2007, 36, 640–647. [Google Scholar] [CrossRef] [Green Version]
  56. Petrova, N.S.; Zenkova, M.A.; Chernolovskaya, E.L. Structure—Functions Relations in Small Interfering RNAs. Pract. Appl. Biomed. Eng. 2013. [Google Scholar] [CrossRef] [Green Version]
  57. Walton, S.P.; Wu, M.; Gredell, J.A.; Chan, C. Designing highly active siRNAs for therapeutic applications. FEBS J. 2010, 277, 4806–4813. [Google Scholar] [CrossRef] [PubMed]
  58. Shao, Y.; Chan, C.Y.; Maliyekkel, A.; Lawrence, C.E.; Roninson, I.B.; Ding, Y. Effect of target secondary structure on RNAi efficiency. RNA 2007, 13, 1631–1640. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Bajan, S.; Hutvagner, G. RNA-Based Therapeutics: From Antisense Oligonucleotides to miRNAs. Cells 2020, 9, 137. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Naito, Y.; Ui-Tei, K. siRNA Design Software for a Target Gene-Specific RNA Interference. Front. Genet. 2012, 3, 102. [Google Scholar] [CrossRef] [Green Version]
  61. Naito, Y.; Ui-Tei, K.; Nishikawa, T.; Takebe, Y.; Saigo, K. siVirus: Web-based antiviral siRNA design software for highly divergent viral sequences. Nucleic Acids Res. 2006, 34, W448–W450. [Google Scholar] [CrossRef] [Green Version]
  62. Naito, Y.; Yoshimura, J.; Morishita, S.; Ui-Tei, K. siDirect 2.0: Updated software for designing functional siRNA with reduced seed-dependent off-target effect. BMC Bioinform. 2009, 10, 392. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Reynolds, A.; Leake, D.; Boese, Q.; Scaringe, S.; Marshall, W.S.; Khvorova, A. Rational siRNA design for RNA interference. Nat. Biotechnol. 2004, 22, 326–330. [Google Scholar] [CrossRef] [PubMed]
  64. Amarzguioui, M.; Prydz, H. An algorithm for selection of functional siRNA sequences. Biochem. Biophys. Res. Commun. 2004, 316, 1050–1058. [Google Scholar] [CrossRef] [PubMed]
  65. Ichihara, M.; Murakumo, Y.; Masuda, A.; Matsuura, T.; Asai, N.; Jijiwa, M.; Ishida, M.; Shinmi, J.; Yatsuya, H.; Qiao, S.; et al. Thermodynamic instability of siRNA duplex is a prerequisite for dependable prediction of siRNA activities. Nucleic Acids Res. 2007, 35, e123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Lee, H.S.; Ahn, J.; Jun, E.J.; Yang, S.; Joo, C.H.; Kim, Y.K.; Lee, H. A novel program to design siRNAs simultaneously effective to highly variable virus genomes. Biochem. Biophys. Res. Commun. 2009, 384, 431–435. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Multiple sequence alignment of selected MERS-CoVs (Nucleotide sequences start. 1……300,000 bp. Due to space limitations, we have selected the conserved regions only to display in the figure).
Figure 1. Multiple sequence alignment of selected MERS-CoVs (Nucleotide sequences start. 1……300,000 bp. Due to space limitations, we have selected the conserved regions only to display in the figure).
Molecules 26 02610 g001
Figure 2. Secondary structure of siRNA MERS-CoV genome predicted by using RNAfold software and were visualized in forna Figure 3. The proposed centroid secondary structure of the MERS-CoV-ORF1ab gene with nucleotide positions was predicted using RNAfold software and downloaded in PNG format. The phosphodiester backbone is shown in orange color and the nucleotides are shown in the sticks and filled rings with elemental coloring; Green, red, yellow and blue. The color of the secondary structure indicates; green: C (canonical helices), red: G (junctions), yellow: A, blue: U hairpin loops. The structure shows all the nucleotides and their bindings with each other and forming stem, hairpin loops, canonical helices, and junctions. Chemically modified nucleotides are unpaired, in A-U or G-C pairs at helix ends, in G-U pairs anywhere, or adjacent to G-U pairs.
Figure 2. Secondary structure of siRNA MERS-CoV genome predicted by using RNAfold software and were visualized in forna Figure 3. The proposed centroid secondary structure of the MERS-CoV-ORF1ab gene with nucleotide positions was predicted using RNAfold software and downloaded in PNG format. The phosphodiester backbone is shown in orange color and the nucleotides are shown in the sticks and filled rings with elemental coloring; Green, red, yellow and blue. The color of the secondary structure indicates; green: C (canonical helices), red: G (junctions), yellow: A, blue: U hairpin loops. The structure shows all the nucleotides and their bindings with each other and forming stem, hairpin loops, canonical helices, and junctions. Chemically modified nucleotides are unpaired, in A-U or G-C pairs at helix ends, in G-U pairs anywhere, or adjacent to G-U pairs.
Molecules 26 02610 g002
Figure 3. The proposed centroid secondary structure of the MERS-CoV-ORF1ab gene with nucleotide positions was predicted using RNAfold software and downloaded in PNG format. The phosphodiester backbone is shown in orange color and the nucleotides are shown in the sticks and filled rings with elemental coloring; Green, red, yellow and blue. The color of the secondary structure indicates; green: C (canonical helices), red: G (junctions), yellow: A, blue: U hairpin loops. The structure shows all the nucleotides and their bindings with each other and forming stem, hairpin loops, canonical helices, and junctions. Chemically modified nucleotides are unpaired, in A-U or G-C pairs at helix ends, in G-U pairs anywhere, or adjacent to G-U pairs.
Figure 3. The proposed centroid secondary structure of the MERS-CoV-ORF1ab gene with nucleotide positions was predicted using RNAfold software and downloaded in PNG format. The phosphodiester backbone is shown in orange color and the nucleotides are shown in the sticks and filled rings with elemental coloring; Green, red, yellow and blue. The color of the secondary structure indicates; green: C (canonical helices), red: G (junctions), yellow: A, blue: U hairpin loops. The structure shows all the nucleotides and their bindings with each other and forming stem, hairpin loops, canonical helices, and junctions. Chemically modified nucleotides are unpaired, in A-U or G-C pairs at helix ends, in G-U pairs anywhere, or adjacent to G-U pairs.
Molecules 26 02610 g003
Figure 4. The accessibility to the target site for each nucleotide in the sequence. The proposed target accessibility plot of different siRNAs (siRNAs 1–21) was determined by RNAxs online software. Each full-length sequence target was submitted to the RNAxs online server. The target accessibility probability profile for each site targeted by the siRNA is displayed. Accessibility and targeted binding position of each siRNA have been shown as a dotted line. Each siRNA has different binding positions in the viral genome. X axis represents the efficiency and Y axis represents the binding position of siRNAs in the viral genome.
Figure 4. The accessibility to the target site for each nucleotide in the sequence. The proposed target accessibility plot of different siRNAs (siRNAs 1–21) was determined by RNAxs online software. Each full-length sequence target was submitted to the RNAxs online server. The target accessibility probability profile for each site targeted by the siRNA is displayed. Accessibility and targeted binding position of each siRNA have been shown as a dotted line. Each siRNA has different binding positions in the viral genome. X axis represents the efficiency and Y axis represents the binding position of siRNAs in the viral genome.
Molecules 26 02610 g004aMolecules 26 02610 g004bMolecules 26 02610 g004c
Figure 5. Graphical representation of Ct value of (a) real time-PCR result-cell lysate and (b) real time-PCR result-cell supernatant.
Figure 5. Graphical representation of Ct value of (a) real time-PCR result-cell lysate and (b) real time-PCR result-cell supernatant.
Molecules 26 02610 g005
Table 1. List of predicted siRNAs from MERS-CoV ORF1ab gene (KF958702).
Table 1. List of predicted siRNAs from MERS-CoV ORF1ab gene (KF958702).
S.N.Position of siRNA in the Genome
(Start–End)
Target SequencePredicted RNA Oligo Sequences (5′→3′)Minimum Free Energy (MFE- kcal/mol) and Frequency of Thermodynamic Ensemble (%)
1791–813agcaatctatttttactattaat UAAUAGUAAAAAUAGAUUGCU
CAAUCUAUUUUUACUAUUAAU
−17.96, 65.21
21615–1637atggataatgctattaatgttgg AACAUUAAUAGCAUUAUCCAU
GGAUAAUGCUAUUAAUGUUGG
−21.80, 77.37
31910–1932gcgactttatgtctacaattatt UAAUUGUAGACAUAAAGUCGC
GACUUUAUGUCUACAAUUAUU
−22.02,69.74
44018–4040gacactttagatgatatcttaca UAAGAUAUCAUCUAAAGUGUC
CACUUUAGAUGAUAUCUUACA
−22.62,69.54
55597–5619atgctattagtttgagttttaat UAAAACUCAAACUAAUAGCAU
GCUAUUAGUUUGAGUUUUAAU
−21.51, 83.64
65598–5620tgctattagtttgagttttaata UUAAAACUCAAACUAAUAGCA
CUAUUAGUUUGAGUUUUAAUA
−19.74, 57.91
75819–5841gagctagtttgcgtcaaattttt AAAUUUGACGCAAACUAGCUC
GCUAGUUUGCGUCAAAUUUUU
−24.28, 53.63
89495–9517ctctaatatctttgttattaaca UUAAUAACAAAGAUAUUAGAG
CUAAUAUCUUUGUUAUUAACA
−17.97, 54.45
99533–9555ctcttagaaactctttaactaat UAGUUAAAGAGUUUCUAAGAG
CUUAGAAACUCUUUAACUAAU
−22.37, 64.54
1013,605–13,627tggtttgattttgttgaaaatcc AUUUUCAACAAAAUCAAACCA
GUUUGAUUUUGUUGAAAAUCC
−18.34, 35.22
1114,005–14,027acgcaaattgcgttaattgtact UACAAUUAACGCAAUUUGCGU
GCAAAUUGCGUUAAUUGUACU
−22.34, 79.46
1214,389–14,411tggtatctaaaggtttctttaag UAAAGAAACCUUUAGAUACCA
GUAUCUAAAGGUUUCUUUAAG
−22.04, 67.95
1316,177–16,199gtcttgtattcggcttatacaag UGUAUAAGCCGAAUACAAGAC
CUUGUAUUCGGCUUAUACAAG
−26.53, 58.68
1416,217–16,239tccttctatagttgaatttaata UUAAAUUCAACUAUAGAAGGA
CUUCUAUAGUUGAAUUUAAUA
−20.24, 48.81
1517,283–17,305gtctacaataataaattgttagc UAACAAUUUAUUAUUGUAGAC
CUACAAUAAUAAAUUGUUAGC
−17.87, 75.42
1617,583–17,605aacaacattaacagatttaatgt AUUAAAUCUGUUAAUGUUGUU
CAACAUUAACAGAUUUAAUGU
−19.59, 62.23
1718,028–18,050ctctacaattaggattttcaact UUGAAAAUCCUAAUUGUAGAG
CUACAAUUAGGAUUUUCAACU
−22.08, 53.94
1819,806–19,828ttgtataagaaagtcaataatga AUUAUUGACUUUCUUAUACAA
GUAUAAGAAAGUCAAUAAUGA
−19.97, 64.53
1920,090–20,112ctcaactattcataactatttta AAAUAGUUAUGAAUAGUUGAG
CAACUAUUCAUAACUAUUUUA
−19.63, 42.01
2020,498–20,520tgccaatatgcgtgttatacatt UGUAUAACACGCAUAUUGGCA
CCAAUAUGCGUGUUAUACAUU
−25.98, 74.21
2120,948–20,970gggtactattaaagaaaatatag AUAUUUUCUUUAAUAGUACCC
GUACUAUUAAAGAAAAUAUAG
−17.65, 66.76
S.N.: serial number
Table 2. Cytotoxicity 50% (CC50) and inhibitory 50% (IC50) concentrations of the different siRNAs.
Table 2. Cytotoxicity 50% (CC50) and inhibitory 50% (IC50) concentrations of the different siRNAs.
siRNAsCC50 (nM)IC50 (nM)
SupernatantLysate
siRNA-1>100NDND
siRNA-2>100NDND
siRNA-3>100NDND
siRNA-4>100NDND
siRNA-5>100NDND
siRNA-6>10026.2217.56
siRNA-7>100NDND
siRNA-8>100NDND
siRNA-9>100NDND
siRNA-10>100NDND
siRNA-11>100NDND
siRNA-12>100NDND
siRNA-13>100NDND
siRNA-14>100NDND
siRNA-15>100NDND
siRNA-16>10013.8716.25
siRNA-17>100NDND
siRNA-18>100NDND
siRNA-19>10011.125.17
siRNA-20>100NDND
siRNA-21>100NDND
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

El-Kafrawy, S.A.; Sohrab, S.S.; Mirza, Z.; Hassan, A.M.; Alsaqaf, F.; Azhar, E.I. In Vitro Inhibitory Analysis of Rationally Designed siRNAs against MERS-CoV Replication in Huh7 Cells. Molecules 2021, 26, 2610. https://doi.org/10.3390/molecules26092610

AMA Style

El-Kafrawy SA, Sohrab SS, Mirza Z, Hassan AM, Alsaqaf F, Azhar EI. In Vitro Inhibitory Analysis of Rationally Designed siRNAs against MERS-CoV Replication in Huh7 Cells. Molecules. 2021; 26(9):2610. https://doi.org/10.3390/molecules26092610

Chicago/Turabian Style

El-Kafrawy, Sherif Aly, Sayed Sartaj Sohrab, Zeenat Mirza, Ahmed M. Hassan, Fatima Alsaqaf, and Esam Ibraheem Azhar. 2021. "In Vitro Inhibitory Analysis of Rationally Designed siRNAs against MERS-CoV Replication in Huh7 Cells" Molecules 26, no. 9: 2610. https://doi.org/10.3390/molecules26092610

APA Style

El-Kafrawy, S. A., Sohrab, S. S., Mirza, Z., Hassan, A. M., Alsaqaf, F., & Azhar, E. I. (2021). In Vitro Inhibitory Analysis of Rationally Designed siRNAs against MERS-CoV Replication in Huh7 Cells. Molecules, 26(9), 2610. https://doi.org/10.3390/molecules26092610

Article Metrics

Back to TopTop