key: cord-1042165-ekxm1sem
authors: Sohrab, Sayed Sartaj; El-Kafrawy, Sherif Aly; Azhar, Esam Ibraheem
title: Effect of insilico predicted and designed potential siRNAs on inhibition of SARS-CoV-2 in HEK-293 cells
date: 2022-03-16
journal: J King Saud Univ Sci
DOI: 10.1016/j.jksus.2022.101965
sha: 5de7bcc8ec42c43e4dfb67b49585f54ca8876769
doc_id: 1042165
cord_uid: ekxm1sem

OBJECTIVES: The COVID-19 was identified for the first time from the sea food market, Wuhan city, China in 2019 and the pathogenic organism was identified as SARS-CoV-2. Currently, this virus has spread to 223 countries and territories and known as a serious issue for the global human community. Many vaccines have been developed and used for immunization METHODS: We have reported the insilico prediction, designing, secondary structure prediction, molecular docking analysis, and in vitro assessment of siRNAs against SARS-CoV-2. The online bioinformatic approach was used for siRNAs selection and designing. The selected siRNAs were evaluated for antiviral efficacy by using Lipofectamine 2000 as delivery agent to HEK-293 cells. The MTT assay was used for cytotoxicity determination. The antiviral efficacy of potential siRNAs was determined based on the Ct value of q-RT-PCR and the data analysis was done by Prism-GraphPad software. RESULTS: The analyzed data resulted in the selection of only three siRNAs out of twenty-six siRNAs generated by online software. The secondary structure prediction and molecular docking analysis of siRNAs revealed the efficient binding to the target. There was no cellular toxicity observed in the HEK-293 cells at any tested concentrations of siRNAs. The purification of RNA was completed from inoculated cells and subjected to q-RT-PCR. The highest Ct value was observed in siRNA 3 than the others. The results offered valuable evidence and invigorated us to assess the potency of siRNAs by using alone or in combination in other human cells. CONCLUSION: The data generated from this study indicates the significance of in silico prediction and narrow down the potential siRNA' against SARS-CoV-2, and molecular docking investigation offered the effective siRNAs binding with the target. Finally, it is concluded that the online bioinformatics approach provided the prediction and selection of siRNAs with better antiviral efficacy. The siRNA-3 was observed to be the best for reduction of viral RNA in cells.

A new case of unusual pneumonia was observed in infected patients in Wuhan, China in late 2019. The infected patients developed the most common symptoms like, breathing problems, fever, and failure of the multi organ in severe cases. The causative agent was identified and sequenced and, based on the high sequence similarity and phylogenetic relationships with bat coronavirus (SARS-bat virus), it was finally designated as Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2), and this disease was designated as COVID-19 in February 2020 and pandemic was declared on March 11, 2020. Currently, the SARS-CoV-2 has been reported from 223 countries with 379,303,864 cases, 5,693,527 deaths, and 299,352,192 recovery (last accessed on 1.2.2022-10.44 GMT) and represented a serious threat to human health globally with approximately 2-5% mortality rate of infected individuals (Azoulay et al., 2020; Zhou et al., 2020) . The virus is attached to the host cell by Spike (S) protein gene and enters in cells by angiotensin−converting enzyme 2 (ACE2) and affects the critical organs resulting in malfunction of airways, lungs, kidneys, gastrointestinal tract as well as the associated cardiovascular system which finally leads to loss of life. This virus falls under Coronaviridae family, with +ve sense single-stranded RNA with an approximately 26-32 kb genome. Currently, many genetic variants have emerged and spread from circulating virions which are found to be more deleterious in fast spread among human population. The first mutation of SARS-CoV-2 was detected in the S protein gene within, a few months of the evolution of the virus. The small mutations tend to impact the virus phenotypes, biology, pathogenicity, more virulence, infectivity, host adaptability, as well as transmissibility and/or antigenicity. Due to the dynamic nature of SARS-CoV-2, frequent mutations, and new variants, it is difficult to predict and develop effective antiviral therapeutics in near future. An extreme measure to be taken to minimize the patient mortality and loss of life.

Recently, a new variant known as Omicron has been identified from South Africa, but many cases have been reported globally. Currently, five variants of concern have been reported from different locations and they are designated as Alpha (B.1.1.7), Beta (B.1.351), Delta (B.1.617.2), Gamma (P.1) and Omicron (B.1.1.529) (WHO; 2021; Khan et al., 2022; Aleem et al.,2022) . Currently, many vaccines have been developed and approved by WHO by using different platforms against multiple targets of SARS-CoV-2 and being used for vaccination of the human population globally, but still an effective broad spectrum antiviral therapeutics are urgently needed to protect the global human population.

Many antiviral compounds and therapeutic molecules have been identified, designed, developed, and are being investigated at various stages. Many compounds have reached an advanced stage of investigation. In continuation of the antiviral therapy development, a new and advanced technology known as RNA interference (RNAi) can also be applied. This technology has been adopted for developing antiviral therapy against many viral and deadly diseases by regulating and silencing the expression of desired genes using siRNAs/shRNAs and degradation of corresponding mRNAs (Hasan et al., 2014; Sohrab et al., 2020 , Sohrab et al., 2021 El-Kafrawy et al., 2021) . The basic criteria for computational prediction, filtration and selection of siRNAs have been discussed in many reports (Naito et al., 2009; Naito et al., 2012; Sohrab et a., 2018 , Uludag et al., 2020 Donia and Bokhari 2021) . The efficacy of miRNA and siRNA are under the stage of investigation in more than 20 clinal trials. Some of them are found to be the most promising and can be considered against respiratory viruses (Chakraborty et al., 2017; Alyan pharmaceuticals 2020) .

For SARS-CoV-2, the computational identification of siRNAs, aptamers/siRNAs chimera, siRNA-Nanoparticle, siRNA-peptide dendrimer formulation utilizing multiple targets have been reported in many published papers (Panda et al., 2020; Chen, et al., 2020; Bappy et al., 2021; Khanali et al., 2020; Shawan et al., 2021; Chowdhury et al., 2021; Idris et al., 2021; Pandey and Verma, 2021; Wu and Luo, 2021; Rohani et al., 2021; Khaitov et al., 2021; Sohrab et al., 2021a) . But none of them were experimentally evaluated except one tested in Vero cells and found some promising results (Niktab et al., 2021) . Recently, we also have in silico predicted, designed, and experimentally evaluated the siRNAs against MERS-CoV, and promising results were observed to inhibit the MERS-CoV proliferation and replication in different cell lines (Sohrab et al., 2020 , Sohrab et al., 2021b , Sohrab et al., 2021c El-Kafrawy et al., 2021) .

Recently, for SARS-CoV-2, many papers have reported about the designing and selection of siRNAs using many target genes. In another study, the leader protein sequences were used and total of four siRNAs were finally designed, but none of them were experimentally validated (Pandey and Verma 2021; Tolksdorf et al., 2021) . So, by using the similar leader protein sequences of SARS-CoV-2/human/SAU/85791C/2020, we tried to predict, design, and evaluate the potential siRNAs in HEK-293 cells. We have used different concentrations of siRNAs and delivered them to grown HEK-293 cells using Lipofectamine 2000 followed by virus inoculation. The viral RNA was purified and used for q-RT-PCR to determine the reduction of viral RNA based on the Ct value in the HEK-293 cells. The results obtained from this work concluded that only one siRNA-3 has the highest antiviral potency to inhibit and reduce the viral RNA in HEK-293 cells.

The leader protein gene sequences were selected from SIAU isolate (GenBank:MT630432). for target-specific prediction, design, and filtration of siRNAs. This tool applies combined utilization Ui-Tei, Amarzguioui, and Reynolds rules. Generally, the strict criteria should be followed for siRNAs designing and filtration as described in the published paper (Naito et al., 2009; Naito et al., 2012 Sohrab et al., 2018 . The predicted siRNAs were further screened by using NCBI Human genomic and transcript database to avoid any similarities with human genome sequences as well as to filter any off-target sequences.

For the analysis of interaction patterns, the computer-assisted molecular docking and molecular dynamics simulation was performed using a guide strand of predicted siRNAs with leader protein sequences of SARS-CoV-2, to determine the effective binding to the target, using online software HNADOCK (He et al., 2019 ). The 3D model structure was generated for guide strand and leader protein sequences and visualized in PyMOL 2.5.2. Based on the lowest Tm value for guide strand and the highest binding score of leader protein sequences as target were considered for best docking analysis. We have selected the best siRNAs model from the top ten docked models generated and used for further study.

To analyze the interaction between predicted siRNA and target leader protein gene sequences based on the thermodynamics as well as the secondary structure of siRNA, the online software like RNAfold 2.4.18 (http://rna.tbi.univie.ac.at) and Mfold server (http://www.mfold.-rna.albany.edu/) was used. The guide and passenger strand sequences were exported to the given locations of program and secondary structure was predicted and generated.

The predicted secondary structure was visualized and downloaded in Forna format (Kerpedjiev et al.,2015) . The secondary structure of siRNAs provides vital information about the effective binding to the target as well as the possibility of any undesired formation of hairpin structure that affect the binding of siRNAs to the selected target.

The predcited siRNAs were filtered for chemical synthesis as per the results generted by online software; siRNA Direct 2.0, HNADOCK, and RNAfold 2.4.18. The synthesis of potential siRNAs was done by Iintegrated DNA Technologies and used for further analysis for antiviral potency against SARS-CoV-2 in HEK-293 cells.

Before testing the effect of siRNAs on the reduction of virus multiplication and proliferation, we tested them on normal and non-infected HEK-293 cells and the cell viability was compared with normal cells. To perform the siRNAs transfection, cells were initially grown in 96 well plates using standard DMEM media at 37°C with 5 % CO 2 . The The cytotoxicity of each siRNAs was determined by using MTT Assay Kit (ThermoFisher Scientific-USA) as per the manufacture's instruction. In brief, the transfected cells were replenished with freshly prepared DMEM media (100 l/well) as well as 10 l MTT (12 mM) and further incubated for 4 hours at 37°C. The absorbance was measured at 570nm by using SpectraMax i3x imaging cytometer and cytotoxicity was calculated based on the mean OD value recorded from SpectraMax i3x using the standard formula for each siRNAs separately. We used Lipofectamine 2000 alone Lipofectamine and Opti-MEM, only Opti-MEM, and only cells as negative control in this study to evaluate their effect on grown cells.

The antiviral activity of predicted, designed, and selected siRNAs were evaluated in grown HEK-293 cells at different concentrations (0.1-50 nM). The 60-80% confluent (1×10 -4 ) HEK-293 cells were transfected with siRNAs and virus inoculation was performed. The cells were further incubated and observed till full CPE in positive control as described (Sohrab et al., 2020) . The cells were harvested after full CPE and the commercial QIAamp Viral RNA Mini Kit was used to purify the viral RNA from lysed cells.

The QIAamp Viral RNA Mini Kit (Qiagen) was used to extract and purify the RNA from harvested and lysed cells. The q-RT-PCR was used to determine the relative expression of the viral genes in HEK-293 cells. The purified RNA was for q-RT-PCR and the level of viral RNA was determined by using PowerCheck SARS-CoV-2 real-time-PCR kit (Kogenebiotech-Korea) following the manufacturer instructions. The data were statistically analyzed, and figures and graphs were generated by using GraphPad Prism-9.3.0.

The online software siDirect 2.0 generated twenty-six siRNAs from leader protein sequence of SARS-CoV-2/human/SAU/85791C/2020 but, we selected only three siRNAs as per the strict selection criteria followed with different scoring tools, functional, Tm, off-targets reduced threshold, and algorithm scores (Naito et al., 2009; Naito et al., 2012 , Sohrab et al., 2018 . The sequences of siRNAs (guide/passenger) and their target location in the viral genome with seedduplex stability (Tm/°C) have been presented in Table 1 .

The secondary structure of each siRNAs and leader protein sequences was predicted by using online software RNAfold 2.4.18. The results of the secondary structure prediction have been computed and visualized in Forna format with variable MFE structure drawing encoding base-pair probabilities and presented in Fig. 1A , B. A variable range of MFE of each siRNAs was observed. The colouring of base pairs denotes the probability of their binding efficiency. There was no probability and structural complexity observed in any of the predicted siRNAs. The structure showed the improved binding efficacy to the target region which can regulate and silence the function of the desired gene.

The molecular docking analysis was performed by HNADOCK online software to predict and reveal the efficient binding capacity to the target region. The software resulted in ten models and the best model was selected based on the docking score of siRNAs/target and visualized in PyMOL 2.5.2. software (Fig. 2) . The docking root-mean-square deviation (RMSD) score ranged from -285.03/153.01 (Å)-siRNA-1, -261.10/350.30(Å)-siRNA-2, -328.83/360.59(Å)-siRNA-3 (Table 1 ). The root-mean-square deviation (RMSD) is used to define the average distance among the atoms.

The highest binding score denotes the best binding efficiency with the target. The siRNA-3 showed the highest RMSD score -328.83/360.59(Å) with SARS-CoV-2-leader protein sequence, which confirms the best binding efficiency and is thereby expected to be the best siRNA for antiviral activity. The antiviral activity was confirmed by in vitro analysis in HEK-293 cells.

The selected siRNAs were used for evaluation of cytotoxicity in HEK-293 cells. The results showed the dosedependent cytotoxicity of siRNAs in tested cells. The OD and CC50 value of each siRNA was variable at different concentrations of siRNAs (Table 2, Fig. 3 ). There was no significant cytotoxicity in HEK-293 cells observed in any siRNAs tested up to 50nM concentration.

inoculation. The inoculated cells were harvested, and viral RNA was purified after 72 hours of post-inoculation with the full CPE observed in positive control. The Ct value of quantitative RT-PCR indicated the reduction of viral RNA at various concentrations of siRNAs (Table 3 ). The dose-dependent reduction of the level of viral RNA was observed in tested siRNAs. Interestingly, a significant reduction of viral RNA was reported in siRNA-3 than others at 5.0 nM and

A new disease known as COVID-19 was reported in December 2019 from China. The etiological agent was identified as SARS-CoV-2. As of 1.2.2022, this virus has spread to 223 countries and territories with 5,693,527 deaths.

Currently, the pathogenic virus has been fully characterized and sequenced. A significant amount of information generated, and various vaccines have been developed by using different platforms and multiple targets including Spike (S) protein gene and being used for immunization. The S gene plays an important role in host cell attachment. Currently, many mutations have been identified especially in the S gene which favoured the introduction of many variants with different properties which attract the researchers to design and develop antivirals, therapeutic molecules, and targeted drugs against this virus for broad spectrum protection of the human population (Harvey et al., 2021; Greaneyet al., 2021; Khateeb et al., 2021; Starr et al., 2021) .

The RNAi technology has been used against many viruses and significant results have been achieved. The nucleic acid-based treatments can be offered against SARS-CoV-2 as co-treatment due to the unique and high range of target sequence specificity and to achieve broad-spectrum resistance against multiple variants. This technology is being used to regulate by silencing of a particular gene/target using siRNA/shRNA molecules ( Chakraborty et al., 2017) . The predesigned siRNAs against unique target sequences of the viral genome will regulate and terminate the virus replication and viral load in the infected cells. Currently, many siRNAs have been predicted, designed and their efficacy was evaluated against many viruses including MERS-CoV, SARS-CoV-1, SARS-CoV-2 and showed a significant inhibition of virus replication (Sohrab et al., 2018; Sohrab et al., 2020; Uludag et al., 2020; Panda et al., 2020; Chen et al., 2020; Sohrab et al., 2021a; Sohrab et al., 2021b; El-Kafrawy et al., 2021; Khaitov et al., 2021; Idris et al., 2021; Pandey and Verma, 2021; Donia and Bokhari 2021; Bappy et al., 2021; Niktab et al., 2021; Khanali et al., 2021; Shawan et al., 2021; Chowdhury et al., 2021; Wu and KQ, 2021; Rohani, et al., 2021; Khaitov et al., 2021) . For designing of potential siRNAs molecule against a particular target gene, several challenges arise and need to be considered. These include siRNA delivery, stability, immunogenic response, off-target effect, possible cross-reactions with other sequences. To overcome these challenges, scientists have improved the model for selection and filtration of potential siRNAs in many published reports (Naito et al., 2012; Sohrab et al., 2018) . Additionally, the secondary structure prediction as well as thermodynamics of RNA-RNA interaction of predicted siRNAs, play an important role in efficient binding to the target.

In this study, we have followed the strict criteria for prediction and filtration of siRNAs against the leader protein gene of SARS-CoV-2 by applying automated online software and evaluated in the HEK-293 cells. (Naito et al., 2009; Naito et al., 2012; Sohrab et al.,2018) . The online software resulted in twenty-six siRNAs against the target, but only three siRNAs were selected for the evaluation of their efficacy in HEK-293 cells. The selected siRNAs were analyzed by secondary structure prediction and molecular docking simulation for their effective binding to the target. The secondary structure analysis results provide vital information about the probability efficient binding to the specific target. The secondary structure formation of a siRNA may inhibit the regulation and expression and cleavage of a particular target mediated by RNA-induced silencing complex (RSIC). Therefore, the secondary structure prediction and determination of minimum free energy are very critical for the efficacy of siRNAs. The most common algorithms were used by a web server for RNA secondary structure prediction by calculating the MFE (Bellaousov et al.,2013) . Additionally, three different siRNAs were designed against E, M and N protein genes of SARS-CoV-1, and an 80% reduction of target gene expression was observed in treated cells (Shi et al., 2005) . The results of our study should be evaluated further for the efficacy of selected siRNA in other cell lines either by using in combinations or alone as well as by using different delivery methods and transfection reagents soon study. Based on the overall results; it is obvious that the in silico designing, selection, filtration and experimental evaluation of selected siRNAs could be very useful for development of next-generation oligonucleotide-based therapeutics molecules against SARS-CoV-2 to protect the human population globally.

Based on the findings, it is concluded that the siRNAs could be considered as an alternative antiviral molecule against SARS-CoV-2. Additionally, more siRNAs can be designed and evaluated in multiple cell lines to get detail information against SARS-CoV. 

Vir and Alnylam Expand Collaboration to Advance RNAi Therapeutics for the Treatment of Coronavirus Infection, Including Covid-19

Emerging Variants of SARS-CoV-2 And Novel Therapeutics Against Coronavirus (COVID-19)

Increased mortality in patients with severe SARS-CoV-2 infection admitted within seven days of disease onset

Designing potential siRNA molecule for the nucleocapsid(N) gene silencing of different SARS-CoV-2 strains of Bangladesh: Computational approach

RNA structure: Web servers for RNA secondary structure prediction and analysis

Therapeutic miRNA and siRNA: Moving from Bench to Clinic as Next Generation Meicine

Computational Identification of Small Interfering RNA Targets in SARS-CoV-2

A computational approach to design potential siRNA molecules as a prospective tool for silencing nucleocapsid phosphoprotein and surface glycoprotein gene of SARS-CoV-2

RNA interference as a promising treatment against SARS-CoV-2

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

Mapping mutations to the SARS-CoV-2 RBD that escape binding by different classes of

SARS-CoV-2 variants, spike mutations and immune escape

HNADOCK: a nucleic acid docking server for modelling RNA/DNA-RNA/DNA 3D complex structures

A SARS-CoV-2 targeted siRNA-nanoparticle therapy for COVID-19

An aptamer/ siRNA Chimera against the SARS-CoV-2: A Dual Therapeutic Strategy for The Virus Neutralizing and RNA Interfering

Forna (force-directed RNA): Simple and effective online RNA secondary structure diagrams

Silencing of SARS-CoV-2 with modified siRNA-peptide dendrimer formulation

The Omicron (B.1.1.529) variant of SARS-CoV-2 binds to the hACE2 receptor more strongly and escapes the antibody response: Insights from structural and simulation data

Emerging SARS-CoV-2 variants of concern and potential intervention approaches

siRNA Design Software for a Target Gene-Specific RNA Interference. Frontiers in genetics. 3, 102

siDirect 2.0: updated software for designing functional siRNA with reduced seed-dependent off-target effect

Design of advanced siRNA therapeutics for the treatment of COVID-19. Meta gene. 29, 100910. Panda et al., 2020. Prediction of potential small interfering RNA molecules for silencing of the spike gene of SARS-CoV-2. The Indian Journal of Medical Research

An in-silico analysis of effective siRNAs against COVID-19 by targeting the leader sequence of SARS-CoV-2. Advances in cell and gene therapy

DisCoVering potential candidates of RNAi-based therapy for COVID-19 using computational methods

Designing an effective therapeutic siRNA to silence RdRp gene of SARS-CoV-2. Infection, genetics and evolution: journal of molecular epidemiology and evolutionary genetics in infectious diseases

Inhibition of genes expression of SARS coronavirus by synthetic small interfering RNAs

In silico prediction and experimental validation of siRNAs targeting ORF1ab of MERS-CoV in Vero cell line

In Silico Prediction and Designing of Potential siRNAs to be used as Antivirals Against SARS-CoV-2

Designing and evaluation of MERS-CoV siRNAs in HEK-293 cell line

Design and Delivery of Therapeutic siRNAs: Application to MERS-Coronavirus

SARS-CoV-2 RBD antibodies that maximize breadth and resistance to escape

Inhibition of SARS-CoV-2 Replication by a Small Interfering RNA Targeting the Leader Sequence

Prospects for RNAi Therapy of COVID-19

Coronavirus disease 2019 (COVID-19): a clinical update