key: cord-0808107-pf0k4px7 authors: Ahmad Mir, Shabir; Firoz, Ahmad; Alaidarous, Mohammed; Alshehri, Bader; Aziz Bin Dukhyil, Abdul; Banawas, Saeed; Alsagaby, Suliman A; Alturaiki, Wael; Ahmad Bhat, Gulzar; Kashoo, Faizan; Abdel-Hadi, Ahmad M title: Identification of SARS-CoV-2 RNA-dependent RNA polymerase inhibitors from the major phytochemicals of Nigella sativa: An in silico approach date: 2021-09-08 journal: Saudi J Biol Sci DOI: 10.1016/j.sjbs.2021.09.002 sha: 11eef1648d1c093fac2d3b1baac32d7bdddeb39c doc_id: 808107 cord_uid: pf0k4px7 The coronavirus disease 2019 (COVID-19), which emerged in December 2019, continues to be a serious health concern worldwide. There is an urgent need to develop effective drugs and vaccines to control the spread of this disease. In the current study, the main phytochemical compounds of Nigella sativa were screened for their binding affinity for the active site of the RNA-dependent RNA polymerase (RdRp) enzyme of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). The binding affinity was investigated using molecular docking methods, and the interaction of phytochemicals with the RdRp active site was analyzed and visualized using suitable software. Out of the nine phytochemicals of N. sativa screened in this study, a significant docking score was observed for four compounds, namely α-hederin, dithymoquinone, nigellicine, and nigellidine. Based on the findings of our study, we report that α-hederin, which was found to possess the lowest binding energy (–8.6 kcal/mol) and hence the best binding affinity, is the best inhibitor of RdRp of SARS-CoV-2, among all the compounds screened here. Our results prove that the top four potential phytochemical molecules of N. sativa, especially α-hederin, could be considered for ongoing drug development strategies against SARS-CoV-2. However, further in vitro and in vivo testing are required to confirm the findings of this study. Bioinformatics Protein Data Bank; RdRp, RNA-dependent RNA polymerase; RMSD, root mean square deviation; RMSF, root mean square fluctuations; SARS-CoV-2, severe acute respiratory syndrome-coronavirus-2 The coronavirus disease 2019 , caused by the severe acute respiratory syndromecoronavirus (SARS-CoV)-2, has adversely affected the health of people across the globe (Mali et al., 2020) . As of January 16, 2021, there are 92,262,621 confirmed cases of COVID-19 worldwide, including 1,995,037 deaths (https://www.who.int/emergencies/diseases/novelcoronavirus-2019). SARS-CoV-2 belongs to the Coronaviridae family and is closely related to the other members of the family, especially SARS-CoV and the Middle East respiratory syndrome coronavirus (MERS-CoV) (Hui et al., 2020; Mali et al., 2020) . The positive stranded RNA genome of SARS-CoV-2 is surrounded by a lipid envelope which contains the spike proteins as well as the membrane proteins. The spike proteins of SARS-CoV-2 bind to the host cell receptors and help the virus release the viral genome into the host cell where it is translated into two polyproteins and the structural proteins, and replication of the viral genome is initiated (V'kovski et al., 2020) . The two-third of the RNA genome of SARS-CoV-2 encodes viral RNAdependent RNA polymerase (RdRp), the associated accessory proteins, and the two large nonstructural polyproteins. The remaining one-third of the genome codes for the four structural proteins (spike, envelope, membrane and nucleocapsid), and the other helper proteins (Luk et al., 2019; Shankar et al., 2020) . RdRp is very important for replication and transcription of the viral genome and is highly conserved among different RNA viruses. The core protein of RdRp consisting of a single chain of approximately 900 amino acid residues, shows minimal activity. However, enhanced activity is achieved by the attachment of additional key subunits (Ahn et al., 2012; Kirchdoerfer and Ward 2019; Subissi et al., 2014) . Being a very essential enzyme in the life cycle of RNA viruses including SARS-CoV-2, the RdRp has already been targeted in various viral infections, including the Zika virus, hepatitis C virus, and coronaviruses (Elfiky et al., 2013; Elfiky 2016; Elfiky 2017; Elfiky 2019; Elfiky and Elshemey 2016; Elfiky and Elshemey 2018; Elfiky and Ismail 2017; Elfiky and Ismail 2019; Ganesan and Barakat 2017) . The active site of RdRp has has been reported to be highly conserved among the different corona viruses (Doublie and Ellenberger 1998; Elfiky 2020a; Elfiky and Ismail 2018) . Therefore, RdRp is the main drug target for SARS-Cov-2 and other coronaviruses (Aftab et al., 2020; Elfiky 2021b; Elfiky 2020b; Yin et al., 2020) . The virus generally spreads from the infected person through close contact along with the droplets spilled during talking, coughing and sneezing (Chan et al., 2020) . The lack of approved drugs or vaccines for COVID-19 is the main concern of the ongoing pandemic. Therefore, developing a promising vaccine or drug intervention is of prime interest and importance in combating this viral disease (Vardhan and Sahoo 2020) . Various efforts are being carried out by the researchers throughout the world for developing the effective treatment strategy for this infectious disease. As a result of these continued efforts various clinical interventions are on trial worldwide. In addition to drug repurposing and synthesis of new drug molecules, the herbal medicine system has gained global emphasis in providing favorable interventions to combat this pandemic viral disease (Patwardhan et al., 2020; Rastogi et al., 2020; Elfiky 2021a) . One of the herbs, Nigella sativa (also known as black seed or Prophetic medicine) has been reported to possess several pharmacological properties including anti-microbial, anti-inflammatory and immunostimulatory activities (Ahmad et al., 2013; Molla et al., 2019) . N. sativa belongs to the plant family 'Ranunculaceae' and its seeds in have been consumed to treat various diseases and infirmities (Mazaheri et al., 2019; Majeed et al., 2020; Yimer et al., 2019) . Besides its use as a spice and a food preservative, it has been traditionally used as a protective and curative remedy for several disorders (Majeed et al., 2020) . Moreover, the beneficial effects and safety of N. Koshak et al., 2017; Namazi et al., 2018; Sahebkar et al., 2016a; Sahebkar et al., 2016b; Tavakkoli et al., 2017) . Recently, various studies have highlighted the efficacy of N. sativa in the treatment of viral diseases Barakat et al., 2013; Bouchentouf et al., 2020; Dorra et al., 2019; Elfiky, 2021a; Onifade et al., 2013a; Onifade et al., 2013b; Sekiou 2020) . N. sativa has been reported to inhibit cytomegalovirus, human immunodeficiency virus, herpes simplex virus and hepatitis A virus infections (Salem and Hossain, 2000; Barakat et al., 2010) . Some studies in literature have revealed that N. sativa inhibited growth of influenza virus H5N1 and replication of Hepatitis C virus (Oyero et al., 2016; Dorra et al., 2019) , in addition to decreasing the coronavirus load in infected HeLa cells (Ulasli et al., 2014) . Few studies have also focused on the in silico screening of its main phytochemicals against some drug targets of SARS-CoV-2 Bouchentouf et al., Molla et al., 2019; 2020; Maiti et al., 2020) . However, the inhibitory potential of these phytochemicals against the prominent drug target, the RNA-dependent RNA polymerase of SARS-CoV-2 has not yet been reported. Therefore, based on some previous studies in literature (Akram Khan and Afzal, 2016; Ahmad et al., 2020; Bouchentouf et al., Molla et al., 2019; 2020; Maiti et al., 2020) , here we selected the major phytochemicals of N. sativa and screened them for their potential to inhibit the RdRp of SARS-CoV-2 using computer aided molecular docking methods. The crystallographic three-dimensional (3D) structure of the RdRp (PDB ID: 6M71) of SARS- (Goodsell and Olson 1990) . The water molecules were removed, polar hydrogen atoms and the Kollman charges were added to the protein, and saved as a PDBQT file (.pdbqt). A grid box covering the active site pocket of 6M71 was generated and the grid parameters were saved as a TXT file (.txt) for input during docking (Iheagwam and Rotimi 2020; Vardhan and Sahoo 2020). The 3D We did not investigate Lipinski's physicochemical parameters of the ligands here because this information is already available in the literature Bouchentouf and Noureddine 2020). The molecular docking of individual ligands into the target protein was performed using the AutoDock Vina software (Trott O and Olson 2009) RdRp is crucial for the replication of SARS-CoV-2 and is therefore considered a key target for the development of antiviral drugs ( Therefore, to fill the gap, the same strategy has been applied in the current study to screen the N. sativa phytochemicals for their potential to inhibit RdRp of SARS-CoV-2. Our results revealed that out of the total nine compounds of N. sativa that we screened by in silico analysis, four (α-hederin, dithymoquinone, nigellicine and nigellidine) had significant binding affinity towards the active site of RdRp as indicated by their considerably lower binding energies (ΔG values less than the cut-off value of -6 kcal/mol) (Shityakov and Foerster, 2014) . α-hederin possessed the lowest binding energy (ΔG = -8.6 Kcal/mol) in complex with 6M71, followed by dithymoquinone (ΔG = -6.1 kcal/mol), nigellicine (ΔG = -6.1 kcal/mol) and nigellidine (ΔG = -6.0 kcal/mol). The binding affinity (ΔG = -8.6 Kcal/mol) of α-hederin was found to be higher than that of remdesivir (-7.6 kcal/mol), the already known antiviral compound used as a positive control in this study. The binding energies of the other three compounds were lower, but close to that of remdesivir (Table 1) . In the current study we have used the recently determined X-Ray crystal structure of RdRp (PDB ID: 6M71) for screening the selected phytochemicals. According to the available literature, the key residues of RdRp involved in the interaction include; Y618, C622, N691, N695, M755, I756, L757, L758, S759, D760, D761, A762, V763, E811, F812, C813 and S814 (numbering for PDB ID: 6M71). The key residues in the active site are D761 and D762. These are involved in the actual reaction of the RdRp enzyme (Elfiky 2021b; Elfiky 2020b; Yin et al., 2020) . We used the RdRp holoenzyme for the initial screening of the selected phytochemicals. The structure of RdRp (PDB ID: 6M71) and the binding of high-affinity phytochemicals/ligands to its active site is shown in Fig. 1 . The 3D and 2D visualization of the interaction of the top four phytochemicals of N. sativa with the active site residues of RdRp is shown in Fig. 2 . The phytochemicals α-hederin and negillicine interacted with the active site residues through four and three conventional hydrogen bonds, respectively, whereas the other two phytochemicals, dithymoquinone and negillidine interacted with the active site residues through two and one hydrogen bonds, respectively. There are also some non-bonded interactions, such as van der Waals forces, pi-alkyl, pi-cation, etc., between the phytochemicals and active site amino acid residues of RdRp as depicted in the 2D illustration (Fig. 2) . The control drug, remdesivir is a nucleotide analog that binds to RdRp in a similar manner as a nucleotide does, thereby inhibiting viral RdRp activity through RNA chain termination. The interaction of various active site amino acid residues of RdRp with remdesivir is available in the literature (Elfiky 2020b; Yin et al., 2020) . In the present study, we observed that upon in silico screening of the selected phytochemicals, the top four compounds (α-hederin, dithymoquinone, negillicine and negillidine), especially α-hederin efficiently bound to the binding pocket of RdRp and interacted with one or more interacting residues present in the RNA binding tunnel of RdRP (Fig. 2) which might lead to the inhibition of its activity. Our results support those of earlier studies and suggest that N. sativa and its phytochemicals are worth studying further and could be recommended as an antiviral herbal supplement against COVID-19. Based on its lowest docking energy on docking with RdRp (PDB ID: 6M71) and its efficient interaction with the active site residues of 6M71, we report here α-hederin as the most potent inhibitor of SARS-CoV-2 RdRp, among the nine compounds screened in this study. For further confirmation, we docked the top-ranked phytochemicals (α-hederin, dithymoquinone, negillicine and negillidine) with a single chain of RdRp (Chain A, the core enzyme). We observed that α-hederin, dithymoquinone, negillidine and negillicine bind into the active site of RdRp (Chain A) with a docking score of -8.6, -6.7, -6.6 and -6.1 kcal/mol, respectively. The stability of the docking complex of the top-ranked compound (α-hederin) was confirmed by molecular dynamic simulation using the Molecular Dynamics on Web (MDWEB) server (https://mmb.irbbarcelona.org/MDWeb/) and CABS flex server (http://212.87.3.12/CABSflex2/index) as shown in Fig. 3 . In this in silico study, we identified four phytochemicals from N. sativa that have the potential The authors declare that they have no conflict of interest regarding the work and content reported in this manuscript. demonstrating the binding of the high affinity phytochemicals to its active site (binding pocket). Chains A, B, C, and D of RdRp are shown in red, green, tv-green, and lime-green color, respectively, whereas, the ligands, α-hederin, dithymoquinone, negillicine, negillidine, and remdesivir, are shown in magenta, cyan, yellow, orange, and blue color, respectively. Analysis of SARS-CoV-2 RNA-dependent RNA polymerase as a potential therapeutic drug target using a computational approach A review on therapeutic potential of Nigella sativa: A miracle herb Molecular docking, simulation and MM-PBSA studies of Nigella sativa compounds: A computational quest to identify potential natural antiviral for COVID-19 treatment Biochemical characterization of a recombinant SARS coronavirus nsp12 RNA-dependent RNA polymerase capable of copying viral RNA templates Chemical composition of Nigella sativa Linn: Part 2 Effects of Nigella sativa on outcome of hepatitis C in Egypt Identification of compounds from Nigella sativa as new potential inhibitors of 2019 novel coronavirus (COVID-19): Molecular docking study A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: A study of a family cluster Nigella sativa improves glucose homeostasis and serum lipids in type 2 diabetes: A systematic review and meta-analysis Discovery Studio Visualizer, v17.2.0.16349; Dassault Systèmes Pymol: An open-source molecular graphics tool. CCP4 Newsl on Protein Crystallogr Evaluation of antiviral and antioxidant activity of selected herbal extracts The mechanism of action of T7 DNA polymerase Zika viral polymerase inhibition using anti-HCV drugs both in market and under clinical trials Zika virus: Novel guanosine derivatives revealed strong binding and possible inhibition of the polymerase Novel guanosine derivatives as anti-HCV NS5b polymerase: A QSAR and molecular docking study Anti-HCV, nucleotide inhibitors, repurposing against COVID-19 Ribavirin, remdesivir, sofosbuvir, galidesivir, and tenofovir against SARS-CoV-2 RNA-dependent RNA polymerase (RdRp): A molecular docking study Natural products may interfere with SARS-CoV-2 attachment to the host cell SARS-CoV-2 RNA dependent RNA polymerase (RdRp) targeting: an in silico perspective Molecular modeling comparison of the performance of NS5b polymerase inhibitor (PSI -7977) on prevalent HCV genotypes IDX -184 is a superior HCV direct-acting antiviral drug: A QSAR study Molecular dynamics simulation revealed binding of nucleotide inhibitors to ZIKV polymerase over 444 nanoseconds Molecular modeling and docking revealed superiority of IDX-184 as HCV polymerase inhibitor Molecular docking revealed the binding of nucleotide/side inhibitors to Zika viral polymerase solved structures Molecular dynamics and docking reveal the potency of novel GTP derivatives against RNA-dependent RNA polymerase of genotype 4a HCV Applications of computer-aided approaches in the development of hepatitis C antiviral agents Automated docking of substrates to proteins by simulated annealing The influency of Nigella sativa for asthma control: A meta-analysis The continuing 2019-nCoV epidemic threat of novel coronaviruses to global health -The latest 2019 novel coronavirus outbreak in Wuhan Computer-aided analysis of multiple SARS-CoV-2 therapeutic targets: Identification of potent molecules from African medicinal plants Structure of the SARS-CoV nsp12 polymerase bound to nsp7 and nsp8 co-factors Nigella sativa supplementation improves asthma control and biomarkers: A randomized, double-blind, placebo-controlled trial Molecular epidemiology, evolution and phylogeny of SARS coronavirus Active-site Molecular docking of Nigellidine with nucleocapsid-NSP2-MPro of COVID-19 and to human IL1R-IL6R and strong antioxidant role of Nigella-sativa in experimental rats Nigella sativa L.: Uses in traditional and contemporary medicines-An overview The rise of new coronavirus infection-(COVID-19): A recent update A comprehensive review of the physicochemical, quality, and nutritional properties of Nigella sativa oil A review on antiviral effects of nigella sativa L The effects of Nigella sativa L. on obesity: A systematic review and meta-analysis Open Babel: An Open chemical toolbox Nigella sativa concoction induced sustained seroreversion in HIV patient Effectiveness of a herbal remedy in six HIV patients in Nigeria Selective inhibition of hepatitis C virus replication by alpha-zam, a nigella sativa seed formulation Ayurveda rasayana in prophylaxis of COVID-19 Potential of inhibiting the receptor binding mechanism of sarscov-2 using phytochemical extracts of medicinal herb; moleculer docking study COVID-19 pandemic: A pragmatic plan for ayurveda intervention Nigella sativa (black seed) effects on plasma lipid concentrations in humans: A systematic review and metaanalysis of randomized placebo-controlled trials A systematic review and meta-analysis of randomized controlled trials investigating the effects of supplementation with Nigella sativa (black seed) on blood pressure Protective effect of black seed oil from Nigella sativa against murine cytomegalovirus infection Molecular docking analysis of selected natural products from plants for inhibition of SARS-CoV-2 main protease In-silico identification of potent inhibitors of COVID-19 main protease (Mpro) and Angiotensin Converting Enzyme 2 (ACE2) from Natural Products: Quercetin, Hispidulin, and Cirsimaritin Exhibited Better Potential Inhibition than Hydroxy-Chloroquine Against COVID-19 Main Protease Active Site and ACE2. ChemRxiv. Cambridge: Cambridge Open Engage An update on COVID-19: SARS-CoV-2 life cycle, immunopathology, and BCG vaccination In silico predictive model to determine vector-mediated transport properties for the blood-brain barrier choline transporter Essential oils as antiviral agents. Potential of essential oils to treat SARS-CoV-2 infection: An in-silico investigation One severe acute respiratory syndrome coronavirus protein complex integrates processive RNA polymerase and exonuclease activities Review on clinical trials of black seed (Nigella sativa) and its active constituent, thymoquinone AutoDockVina: Improving the speed and accuracy of docking with a new scoring function, ecient optimization, and multithreading The effects of Nigella sativa (Ns), Anthemishyalina (Ah) and Citrus sinensis (Cs) extracts on the replication of coronavirus and the expression of TRP genes family Coronavirus biology and replication: implications for SARS-CoV-2 Searching inhibitors for three important proteins of COVID-19 through molecular docking studies Nigella sativa L. (Black Cumin): A Promising Natural Remedy for Wide Range of Illnesses Structural basis for inhibition of the RNA-dependent RNA polymerase from SARS-CoV-2 by remdesivir The authors are highly thankful to the Deanship of Scientific Research (DSR) at Majmaah University for their support in proofreading the manuscript through professional editors at Editage, a division of Cactus Communications. ☒ The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.☐The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: