key: cord-1046817-k1jqkeoj authors: Chen, Wei; Feng, Pengmian; Liu, Kewei; Wu, Meng; Lin, Hao title: Computational Identification of Small Interfering RNA Targets in SARS-CoV-2 date: 2020-04-15 journal: Virol Sin DOI: 10.1007/s12250-020-00221-6 sha: 4928fa8e83ecacd1284752353a8d5b9c0614a542 doc_id: 1046817 cord_uid: k1jqkeoj nan dsRNA is a phenomenon of homology-dependent gene silencing and may play certain roles in affecting the process of virus expression and proliferation. Recently, several reports have demonstrated the use of RNAi in blocking virus infection and replication in animal cells (Ge et al. 2003) , suggesting that the small interfering RNA (siRNA, 21-25 nt long) plays an important role in RNAi-related gene silencing pathways (Elbashir et al. 2001) . Progress has been made in anti-HIV and anti-HCV drug design by applying the method of RNA interference (Wilson et al. 2003) . The effectiveness of siRNA for inhibiting SARS coronavirus genes expression was also demonstrated by Shi et al. (2005) . Besides silencing the targeted genes, the siRNAs can also inhibit the replication of the virus. For example, it has been demonstrated that, by targeting the Leader sequence of SARS-CoV, the siRNA demonstrate a strong inhibitory effect on SARS-CoV replication (Li et al. 2005) . More recently, a CRISPR/Cas13d system was proposed for the treatment of SARS-COV-2 (Nguyen et al. 2020) . These results indicate that both RNAi and CRISPR/ Cas technology might become potential therapeutic approaches for treating viral diseases. Accordingly, as complementary to the CRISPR/Cas13d system, we proposed an RNAi based strategy that might interfere the gene expression and block the replication of SARS-COV-2. The main idea of this strategy is to search for siRNA targets in the virus genome, which will be recognized and cleaved by the RNA-induced silencing complex (RISC). In this work, we performed theoretical predictions of the potential siRNA targets in the virus genome. We firstly collected the representative SARS-COV-2 genome (MN908947, https://www.ncbi.nlm.nih.gov/nuccore/ MN908947) and the mutation information of the SARS-COV-2 genomes from the 2019nCoVR database (Zhao et al. 2020) , which is available at https://bigd.big.ac.cn/ ncov/. The 2019nCoVR database not only integrates genomic and proteomic sequences of SARS-COV-2 from different resources, but also provides a series of scientific services, such as variation visualization, variation annotations, AI diagnosis, etc. Next, we folded the SARS-COV-2 genome (MN908947) in a window of 3000 nucleotides with the step of 1500 nucleotides by using RNAstructure (version 4.5) program (Bellaousov et al. 2013) . Only those 21-25 nt long non-base-paired regions can be served as the potential targets of siRNA (Huang et al. 2008) , which is called free segments. The long non-base-paired region containing one or several short stems (total length of stems 1-3 base pairs), called quasi-free segments (Ji and Luo 2004) , was also considered in the present work. A given RNA sequence segment may have different configurations of secondary structure with lower free energy. The total frequency of a segment occurring in nonbase-paired region of different folds (20 folds are selected for each segment) is called appearance rate (AR). If each quasi-free case is multiplied by a reduced factor in numeration, namely, by 0.9 for 1 base pair, 0.8 for 2 base pair, and 0.7 for 3 base pairs (base pairs may be continuous in structure or disconnected) then the total number of folds is called reduced appearance rate (RAR) (Ji and Luo 2004) . To guarantee the safety of the designed drug, we further performed alignment of the free and quasi-free segments with human genome (hg 38) by using BLAST and deleted the matching ones in siRNA target candidates. Finally, we obtained nine potential siRNA targets in the SARS-COV-2 genome (MN908947). The information about their position and region in the virus genome, length, AR and RAR was provided in Table 1 . In addition, we also analyzed the mutations of the target sequences by comparing all the 143 high quality strains in the 2019nCoVR database (as of March 15, 2020). SNP were found in two of the nine target sequences (indicated by bold character in Table 1 ). For the potential target 'AAUAGUUUAAAAAUUACAGAAGA', only one SNP was found in the strain BetaCoV/Wuhan/HBCDC-HB-05/ 2020, which is a coding_sequence_variant that changes the coding sequence. For 'CAACUAUAAAUUAAACA-CAGA', the SNP was found in the strain BetaCoV/Singapore/6/2020 and BetaCoV/Singapore/2/2020, respectively, which is a missense_variant that changes G to A resulting in a different amino acid sequence. These results indicate that the selected targets are conserved among the existing SARS-COV-2 genomes. Although there are still some challenges that needed to be overcome for the clinic applications of siRNA, progresses have been made to solve the fundamental problems, such as off-target effects and effective delivery. For example, the position-specific chemical modification of siRNAs could can significantly reduce off targeting; safe and effective in vivo delivery systems have also been developed, such as nanoparticles, cationic lipids, antibodies, cholesterol, aptamers delivery strategies. Therefore, we hope that the above results would be useful in drug design and treatments against SARS-COV-2. Animal and Human Rights Statement This article does not contain any studies with human or animal subjects performed by any of the authors. The bold and underlined characters indicate the SNP found in different strains. Virologica Sinica Rnastructure: web servers for RNA secondary structure prediction and analysis The 2019-new coronavirus epidemic: evidence for virus evolution Rna interference (rnai)-based therapeutics: delivering on the promise? Computational antisense oligo prediction with a neural network model Duplexes of 21-nucleotide rnas mediate rna interference in cultured mammalian cells Rna interference of influenza virus production by directly targeting mrna for degradation and indirectly inhibiting all viral RNA transcription Small interfering rna therapy in cancer: mechanism, potential targets, and clinical applications Prediction for target sites of small interfering RNA duplexes in sars coronavirus Sirna targeting the leader sequence of sars-cov inhibits virus replication Virus against virus: a potential treatment for 2019-ncov (sars-cov-2) and other RNA viruses Novel coronavirus: from discovery to clinical diagnostics Inhibition of genes expression of sars coronavirus by synthetic small interfering rnas Rna interference blocks gene expression and RNA synthesis from hepatitis c replicons propagated in human liver cells The 2019 novel coronavirus resource China Novel Coronavirus I, Research T (2020) A novel coronavirus from patients with pneumonia in china Acknowledgements The authors would like to express gratitude to the Editor and anonymous reviewers for their constructive comments. This work was supported by the National Nature Scientific Foundation of China (31771471, 61772119) and the Natural Science Foundation for Distinguished Young Scholar of Hebei Province (No. C2017209244).