Abstract
Background: The current pandemic of the coronavirus disease-2019 (COVID-19) has badly affected our life during the year 2020. SARS-CoV-2 is the primary causative agent of the newly emerged pandemic. Natural flavonoids, Terpenoid and Thymoquinone are tested against different viral and host-cell protein targets. These natural compounds have a good history in treating Hepatitis C Virus (HCV) and Human Immunodeficiency Virus (HIV).
Methods: Molecular docking combined with cytotoxicity and plaque reduction assay is used to test the natural compounds against different viral (Spike, RdRp, and Mpro) and host-cell (TMPRSS II, keap 1, and ACE2) targets.
Results: The results demonstrate the binding possibility of the natural compounds (Thymol, Carvacrol, Hesperidine, and Thymoquinone) to the viral main protease (Mpro). Some of these natural compounds were approved to start clinical trail from Egypt Center for Research and Regenerative Medicine ECRRM IRB (Certificate No.IRB00012517)
Conclusion: Development of an effective anti-viral for SARS-CoV-2 could help to limit the viral load. Benchmarking testing of those natural compounds against other potential antivirals for SARS-CoV-2 with alternative mechanisms of action would thus be important as soon as practicable.
Keywords: COVID-19; SARS-CoV-2; Natural Products; Antivirals; Docking; Plaque Reduction Assay
Abbreviations: COVID-19: Coronavirus Disease 2019; SARS-CoV-2: Severe Acute Respiratory Syndrome Coronavirus 2; S Protein: Spike Protein; RBD: Receptor Binding Domain; RdRp : RNA Dependent RNA Polymerase; IC50: Inhibition Concentration; NS: Nigella Sativa; TQ: Thymoquinone; HCV: Hepatitis C Virus; ACE2: Angiotensin Converting Enzyme 2; HIV: Human Immunodeficiency Virus; TMPRSS: Transmembrane Protease Serine
Introduction
By the end of 2019, an outbreak of a novel coronavirus (SARSCoV- 2) in Wuhan city in China was detected and spread all over the world [1]. On 10th October 2020, the number of confirmed cases of coronavirus disease (COVID-19) reached more than 37 M worldwide with +1M total death, as reported in the World Health Organization (WHO). The associated pneumonia with the novel viral infection, COVID-19, is divided into three phases that correspond to different clinical stages of the disease [2]. Stage 1 is the asymptomatic stage, where the inhaled virus binds to nasal epithelial cells in the nasal cavity and starts replicating. Stage 2 is the upper airway stage, where the virus propagates, migrates down the respiratory tract along the conducting airways, and a more robust innate immune response is triggered. About 20% of the infected patients will progress to stage 3 disease and develop pulmonary infiltrates. Some of these patients will develop a very severe disease as the virus reaches alveoli in the lung and infects alveolar type II cells in peripheral and sub-pleural areas of the lung [3]. SARS-CoV-2 propagates within type II cells, large numbers of viral particles are released, and the cells undergo apoptosis and die. Therefore, the spectrum of symptomatic COVID-19 ranges from mild respiratory tract infection to severe pneumonia that may progress to fatal respiratory syndrome and multi-organ malfunctions [2].
Thymol, known as 2-isopropyl-5-methylphenol, is a natural mono- terpenoid phenol derivative of Cymene. It is found in thyme oil and extracted from Thymus vulgaris [4]. Thymol is a white crystalline substance that has a pleasant aromatic odor. Thymol also provides the distinctive and robust flavor of the culinary herb thyme. Carvacrol is known as monoterpenoid phenol and extracted from Oregano. It has a characteristic pungent and warm odor [5]. Hesperidine is a common flavone glycoside found in citrus fruit such as lemons and sweet oranges [6,7]. It has several pharmacological activities such as antihyperlipidemic, anti-atherogenic, venotonic, antidiabetic, cardioprotective, anti- antihypertensive, and inflammatory actions [6,7]. The anti-inflammatory activity of hesperidin was mainly attributed to its antioxidant defense mechanism and suppression of pro-inflammatory cytokine production [6]. Hesperidin exhibited antiviral activity against the influenza virus through a significant reduction of viral replication.
Nigella sativa (NS) contains many active molecules, such as Thymoquinone (TQ), two forms of alkaloids: isoquinoline alkaloid that includes nigellicimine, nigellicimine n-oxide and pyrazol alkaloid that includes nigellidine and nigellicine [8,9]. TQ is the most abundant constituent in the volatile oil of Nigella sativa seeds, and most of the herb’s properties are attributed to it [10,11]. It has been reported that NS oil can decrease the viral count of HCV in patients received capsules of NS oil (450 mg) three times a day over a 3-month period [12]. Moreover, two clinical studies documented to sustained sero-reversion of the HIV virus over treatment period of 6 to 12 months [13-15]. Molecular docking represents a promising in silico method used to predict the binding affinities of small molecules to proteins as a first step in structure-based drug design [16-21]. This study investigated many active ingredients that showed antiviral activities against SARS-CoV-2, such as Thymol, Carvacrol, Hesperidine, and Thymoquinone. Molecular docking is used to test the binding affinities of these natural product derived compounds against different viral and host cell proteins. Additionally cytotoxicity assay and plaque reduction assay are used to verify their antiviral activity against SARS-CoV-2 collected from Egyptian patients.
Materials and Methods
In silico Testing
Before performing the docking studies, the tested compounds are retrieved from the PubChem database and then prepared using PyMOL software [22,23]. The structures of the proteins are downloaded from the protein data bank [24]. Autodock Tools are used to prepare the docking input files after adding charges [25]. The docking experiments are performed utilizing Autodock Vina software in triplicates to have a clear depiction of the mode of action of compounds against their protein targets [26]. SARS-CoV-2 Spike protein (S), The RNA dependent RNA polymerase (RdRp), the main protease (Mpro) are used as protein targets due to its critical role in maintaining the viral infection [27-31]. Additionally, the hostcell receptors, Transmembrane protease, serine 2 (TMPRSS II), Kelch-like ECH-associated protein 1 (KEAP 1) and the Angiotensin Converting Enzyme 2 (ACE2) are targeted due to its fundamental in viral recognition and maintain infectivity for SARS- CoV-2 [32- 35]. For each compound, ten interactions were generated and one with best binding affinity was selected. PyMOL software was used to represent and analyze the docking complexes.
Experimental Section
All the chemical compounds are purchased from different sources as follow; Thymol purchased from upnature as THYME 100% pure and natural 118 ml, Carvacrol purchased from Zane Hellas as oil of Oregano 30 ml, Hesperidin purchased from sciencebased nutrition as hesperidin methyl chalcone 500 mg -60 veggie caps, and Thymoquinone purchased from prime natural black seed USDA organic 240 ml.
Cytotoxicity Assay: Samples were diluted with Dulbecco’s Modified Eagle’s Medium (DMEM). Stock solutions of the test compounds were prepared in 10 % DMSO in ddH2O. The cytotoxic activity of the extracts were tested in Vero E6 cells by using the 3-(4, 5-dimethylthiazol -2-yl)-2, 5- diphenyltetrazolium bromide (MTT) method [36] with minor modification. Briefly, the cells were seeded in 96 well-plates (100 μl/well at a density of 3×105 cells/ ml) and incubated for 24 hrs at 37 °C in 5%CO2. After 24 hrs, cells were treated with various concentrations of the tested compounds in triplicates. After further 24 hrs, the supernatant was discarded, and cell monolayers were washed with sterile phosphate buffer saline (PBS) 3 times and MTT solution (20 μl of 5 mg/ml stock solution) was added to each well and incubated at 37 °C for 4 hrs followed by medium aspiration. In each well, the formed formazan crystals were dissolved with 200 μl of acidified isopropanol (0.04 M HCl in absolute isopropanol = 0.073 ml HCL in 50 ml isopropanol). Absorbance of formazan solutions were measured at λmax 540 nm with 620 nm as a reference wavelength using a multi-well plate reader. The percentage of cytotoxicity compared to the untreated cells was determined with the following equation.
The plot of % cytotoxicity versus sample concentration was used to calculate the concentration which exhibited 50% cytotoxicity (IC50).
Plaque Reduction Assay
Assay was carried out according to the method of [37] in a six well plate where Vero E6 cells (105 cells / ml) were cultivated for 24 hrs at 37 °C. Sever Acute Respiratory Syndrome Coronavirus (SARSCoV2) virus was diluted to give 103 PFU /well and mixed with the safe concentration of the tested compounds and incubated for 1 hour at 37 °C before being added to the cells. Growth medium was removed from the cell culture plates and the cells were inoculated with (100 μl /well) virus with the tested compounds, After 1 hour contact time for virus adsorption, 3 ml of DMEM supplemented with 2% agarose and the tested compounds was added onto the cell monolayer, plates were left to solidify and incubated at 37 °C till formation of viral plaques (3 to 4 days). Formalin (10%) was added for two hours then plates were stained with 0.1 %crystal violet in distilled water. Control wells were included where untreated virus was incubated with Vero E6 cells and finally plaques were counted and percentage reduction in plaques formation in comparison to control wells was recorded as following:
% Inhibition = [viral count (untreated) - viral count (treated)/viral count (untreated)] x 100
Statistical Analysis
Analysis was performed using Graphpad Prism 8.0.2. Data are represented as mean ± SD and statistical significance was evaluated using one-way ANOVA followed by tukey multiple comparison tests.
Results
Many reports showed that different natural product derived compounds have promising results against inflammation and viral infections including the newly emerged viral infection (SARS-CoV-2) [15,38-40]. Different natural compounds and their derivatives proved its binding affinity against different SARS-CoV-2 and host-cell targets [38,39,41]. Figure 1 shows the 2D structures of the selected natural compounds used in the current study against viral and host-cell protein targets. Three SARS-CoV-2 proteins are targeted in this study including the viral, host-cell recognizing critical element, spike protein (Figure 2A), the vital viral enzyme RNA dependent RNA polymerase (RdRp), responsible for the polymerization of the complement RNA copy, and the Main protease (Mpro) of SARS-CoV-2, which is critical for polyprotein processing (Figure 2B) [42-47]. The spike protein trimers over the virions take different conformations during infection such as the prefusion and postfusion. The prefusion has open and closed conformations in which one or two of the receptor binding domains (RBD) are exposed (open) or immersed into the trimer (close) [30,48]. Table 1 shows the binding affinity calculated using AutoDock Vina software for the docking of the natural compounds (Thymol, Carvacrol, Hesperidine, and Thymoquinone) against the SARS-CoV-2 Mpro as a protein target. The standard compound Chloroquine is used to assess the binding affinity of the natural compounds against the Mpro. As reflected from the values, Carvacrol, Hesperidine, and Thymoquinone show comparable binding affinities (-7.0, -6.9, and - 6.9 kcal/mol, respectively) to SARS-CoV-2 Mpro compared to that of the standard compound (-7.2 kcal/mol). Thymol show slightly higher (worse) binding affinity value (-5.8 kcal/mol) compared to Chloroquine but still able to bind the SARS-CoV-2 Mpro tightly. Figures 3A & 3B show the 3D poses for the docking complexes. The four compounds (Thymol, Carvacrol, Hesperidine, and Thymoquinone) are able to bind to the active site of the Mpro (His41 and Cys145).
The Half-Maximal Inhibitory Concentration (Ic50) of Thymol, Carvacrol, Hesperidin, and Thymoquinone
The effect of different concentrations of the compounds on the cellular proliferation of Vero E6 cell line following 24 h of treatment was determined using MTT assay. Results of the present work revealed a concentration-dependent cytotoxic effect of Thymol (617 μM), Carvacrol (464 μM), Hesperidine (5.58 μM) and Thymoquinone (3.9 μM) on Vero E6 cell line as shown in Figure 4.
Discussion
Our results revealed that Hesperidine was the most affective natural product against SARS-CoV- 2 in Egyptians patients, as it caused inhibition to approximately 100% for infected Vero E6 cell line, Carvacrol, Thymoquinone caused inhibition to about 98%, Chloroquine as a positive control caused inhibition to 98% and Thymol show the lowest inhibition percentage 96% for infected Vero E6 cell line. The results showed significantly difference between natural products effect with (P-value < 0.0001) as shown in Figure 5. These results indicate that Hesperidine and other natural compounds as Carvacrol and Thymoquinone could be therapeutic agent against SARS-CoV-2 in Egypt as treatment resulted in the effective loss of essentially all viral material by time.Eventually, development of an effective anti-viral for SARS-CoV-2, if given to patients early in infection, could help to limit the viral load, prevent severe disease progression and limit person- person transmission. Benchmarking testing of those natural compounds against other potential antivirals for SARS-CoV-2 with alternative mechanisms of action would thus be important as soon as practicable.
Conclusion
These results indicate that Hesperidine and other natural compounds as Carvacrol and Thymoquinone could be therapeutic agent against SARS-CoV-2 in Egypt as treatment resulted in the effective loss of essentially all viral material by time. Eventually, development of an effective anti-viral for SARS-CoV-2, if given to patients early in infection, could help to limit theviral load, prevent severe disease progression and limit person- person transmission. Benchmarking testing of those natural compounds against other potential antivirals for SARS-CoV-2 with alternative mechanisms of action would thus be important as soon as practicable.
Declaration
Ethics Approval and Consent to Participate
Not applicable.
Consent to Publication
Not applicable.
Availability of Data and Material
The docking structures are available upon request from the corresponding author.
Competing Interests
All the authors declare that there is no competing interest in this work.
Funding
Not applicable.
Authors’ Contributions
A. Mohamed Gomaa Seadawy
In Vitro work, Writing, Editing and Molecular Docking
B. Ahmed Fawzy Gad
Molecular docking
C. Mohamed Shamel Eldesoky
Data analysis
D. Bassem Elsayed Elharty
Molecular docking
E. Mostfa Fetooh Mohamed
In vitro study
F. Abdo Elfiky
Molecular docking
G. Aya Ahmed
Molecular docking and Editing
H. Abdel Rahman Nabwi Zekri
Reviewing
Acknowledgements
Not applicable.
References
- Chan JFW, Yuan S, Kok KH, Kelvin KWT, Hin C, et al. (2020) A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster. The Lancet 395(10223): 514-523.
- Wu Z, McGoogan JM (2020) Characteristics of and Important Lessons from the Coronavirus Disease 2019 (COVID-19) Outbreak in China: Summary of a Report of 72 314 Cases From the Chinese Center for Disease Control and Prevention. JAMA 323(13): 1239-1242.
- Wu J, Wu X, Zeng W, Dajing G, Zheng F, et al. (2020) Chest CT Findings in Patients with Coronavirus Disease 2019 and Its Relationship With Clinical Features. Investigative Radiology 55(5): 257-261.
- Ultee A, Slump RA, Steging G, Smid EJ (2000) Antimicrobial activity of carvacrol toward Bacillus cereus on rice. Journal of food protection 63(5): 620-624.
- ULTEE A, SLUMP RA, STEGING G, SMID EJ (2000) Antimicrobial Activity of Carvacrol toward Bacillus cereus on Rice. Journal of food protection 63(5): 620-624.
- Zanwar AA, Badole SL, Shende PS, Hegde MV, Bodhankar SL (2014) Cardiovascular effects of hesperidin: A flavanone glycoside. Polyphenols in human health and disease: Elsevier pp. 989-992.
- Jadeja RN, Devkar RV (2014) Polyphenols and flavonoids in controlling non- alcoholic steatohepatitis. Polyphenols in Human Health and Disease: Elsevier pp. 615-623.
- Mollazadeh H, Hosseinzadeh H (2014) The protective effect of Nigella sativa against liver injury: a review. Iranian journal of basic medical sciences 17(12): 958-966.
- Rahmani AH, Aly SM (2015) Nigella sativa and its active constituents thymoquinone shows pivotal role in the disease’s prevention and treatment. Asian J Pharm Clin Res 8(1): 48-53.
- Smith RD (1980) Avicenna and the Canon of Medicine: a millennial tribute. Western Journal of Medicine 133(4): 367.
- Ahmad A, Husain A, Mujeeb M, Shah Alam Khan, Abdul Kalam N, et al. (2013) A review on therapeutic potential of Nigella sativa: A miracle herb. Asian Pacific journal of tropical biomedicine 3(5): 337-352.
- Ulasli M, Gurses SA, Bayraktar R, Yumrutas O, Oztuzcu S, et al. (2014) The effects of Nigella sativa (Ns), Anthemis hyalina (Ah) and Citrus sinensis (Cs) extracts on the replication of coronavirus and the expression of TRP genes family. Mol Biol Rep 41(3): 1703-1711.
- Salem AM, Bamosa AO, Qutub HO, Rakesh K, Ahmed B, et al. (2017) Effect of Nigella sativa supplementation on lung function and inflammatory mediatorsin partly controlled asthma: a randomized controlled trial. Annals of Saudi medicine 37(1): 64-71.
- Boskabady MH, Javan H, Sajady M, Rakhshandeh H (2007) The possible prophylactic effect of Nigella sativa seed extract in asthmatic patients. Fundamental & clinical pharmacology 21(5): 559-566.
- Mahboubi M (2018) Natural therapeutic approach of Nigella sativa (Black seed) fixed oil in management of Sinusitis. Integrative medicine research 7(1): 27-32.
- Elfiky AA, Ismail AM, Elshemey WM (2020) Recognition of gluconeogenic enzymes; Icl1, Fbp1, and Mdh2 by Gid4 ligase: A molecular docking study. Journal of Molecular Recognition 33(5): e2831.
- Leach A (2001) Molecular Modelling: Principles and Applications (2nd,). Prentice Hall, England.
- Elfiky AA (2020) Novel guanosine derivatives against Zika virus polymerase in silico. Journal of Medical Virology 92(1): 11-16.
- Ganesan A, Barakat K (2017) Applications of Computer-Aided Approaches in The Development of Hepatitis C Antiviral Agents. Expert Opinion on Drug Discovery 12(4): 407-425.
- Hasani HJ, Barakat KH (2016) Protein-Protein Docking: Are We There Yet? Methods and Algorithms for Molecular Docking-Based Drug Design and Discovery: IGI Global pp. 173-195.
- Elfiky AA, Ismail A (2019) Molecular dynamics and docking reveal the potency of novel GTP derivatives against RNA dependent RNA polymerase of genotype 4a HCV. Life Sciences 238: 116958.
- 7.6 V. The PyMOL Molecular Graphics System, Version 1.7.6 Schrödinger, LLC.
- Kim S, Thiessen PA, Bolton EE, Jie C, Gang Fu, et al. (2016) PubChem substance and compound databases. Nucleic acids research 44(D1): D1202-D1213.
- Berman H, Henrick K, Nakamura H (2003) Announcing the worldwide Protein Data Bank. Nat Struct Biol 10(12): 980.
- Morris GM, Huey R, Lindstrom W, Miche FS, Richard KB, et al. (2009) AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J Comput Chem 30(16): 2785-2791.
- Trott O, Olson AJ (2010) AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem 31(2): 455-461.
- Elfiky AA (2020) Reply to a letter to the editor. Life sciences 252: 117715-.
- Abu Melha S, Edrees MM, Riyadh SM, Abdelaziz MR, Elfiky AA, et al. (2020) Clean Grinding Technique: A Facile Synthesis and In Silico Antiviral Activity of Hydrazones, Pyrazoles, and Pyrazines Bearing Thiazole Moiety against SARS- CoV-2 Main Protease (M(pro)). Molecules 25(19): 4565.
- Elfiky AA, Ibrahim IM, Ismail AM, Elshemey WM (2021) A possible role for GRP78 in cross vaccination against COVID-19. Journal of Infection 82(2): 282-327.
- Alaa Ismail AE, Abdo A (2020) SARS-CoV-2 Spike behavior in situ: A Cryo-EM Images for a better understanding of the COVID-19 pandemic. Signal transduction and targeted therapy.
- Elfiky A (2020) SARS-CoV-2 Spike-Heat Shock Protein A5 (GRP78) recognition may be related to the immersed human coronaviruses. Frontiers in Pharmacology.
- Yan R, Zhang Y, Li Y, Xia L, Guo Y, et al. (2020) Structural basis for the recognition of SARS-CoV-2 by full-length human ACE2. Science 367(6485): 1444-1448.
- Hoffmann M, Kleine Weber H, Schroeder S, Nadine K, et al. (2020) SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor. Cell 181(2): 271-280.
- Mendonca P, Soliman KF (2020) Flavonoids Activation of the Transcription Factor Nrf2 as a Hypothesis Approach for the Prevention and Modulation of SARS-CoV-2 Infection Severity. Antioxidants 9(8): 659.
- Mariappan V, SR Rao, Balakrishna Pillai A (2020) Angiotensin-converting enzyme 2: A protective factor in regulating disease virulence of SARS-COV-2. IUBMB Life 72(12): 2533-2545.
- Mosmann T (1983) Rapid colorimetric assay for cellular growth and survival: Application to proliferation and cytotoxicity assays. Journal of Immunological Methods 65(1): 55-63.
- Hayden FG, Cote KM, Douglas RG (1980) Plaque inhibition assay for drug susceptibility testing of influenza viruses. Antimicrobial agents and chemotherapy 17(5): 865-870.
- Adem Ş, Eyupoglu V, Sarfraz I, Rasul A, Amer FZ, et al. (2020) Caffeic acid derivatives (CAFDs) as inhibitors of SARS-CoV-2: CAFDs-based functional foods as a potential alternative approach to combat COVID-19. Phytomedicine pp. 153310.
- Elfiky AA (2020) Natural products may interfere with SARS-CoV-2 attachment to the host cell. Journal of Biomolecular Structure and Dynamics p. 1-10.
- Noor NA, Fahmy HM, Mohammed FF, Elsayed AA, Radwan NM (2015) Nigella sativa amliorates inflammation and demyelination in the experimental autoimmune encephalomyelitis induced Wistar rats. Int J Clin Exp Pathol 8(6): 6269-6286.
- Elfiky AA, Baghdady AM, Ali SA, Ahmed MI (2020) GRP78 targeting: Hitting two birds with a stone. Life Sciences 260: 118317.
- Elfiky AA (2020) Anti-HCV, nucleotide inhibitors, repurposing against COVID-19. Life Science 248: 117477.
- Elfiky AA (2020) Ribavirin, Remdesivir, Sofosbuvir, Galidesivir, and Tenofovir against SARS-CoV-2 RNA dependent RNA polymerase (RdRp): A molecular docking study. Life Sciences 253: 117592.
- Elfiky AA (2020) SARS-CoV-2 RNA dependent RNA polymerase (RdRp) targeting: an in silico perspective. Journal of Biomolecular Structure and Dynamics p. 1-9.
- Jin Z, Du X, Xu Y, Yongjuaing D, Meiqiin L, et al. (2020) Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors. Nature 582(7811): 289-293.
- Ibrahim IM, Abdelmalek DH, Elshahat ME, Elfiky AA (2020) COVID-19 spike- host cell receptor GRP78 binding site prediction. Journal of Infection 80(5): 554-562.
- Korber B, Fischer WM, Gnanakaran S, Yoon H, Theiler H, et al. (2020) Tracking Changes in SARS- CoV-2 Spike: Evidence that D614G Increases Infectivity of the COVID-19 Virus. Cell 182(4): 812-827.e19.
- Ke Z, Oton J, Qu K, Mirko C, Zilla V, et al. (2020) Structures and distributions of SARS-CoV-2 spike proteins on intact virions. Nature.