In silico identification of compounds from Nigella sativa seed oil as potential inhibitors of SARS-CoV-2 targets

Background The growing number of cases, severity and fatality of the COVID-19 pandemic, coupled with the fact that no cure has been found has made infected individuals especially in Africa, to resort to the consumption of different natural products to alleviate their condition. One of such plant materials that have been consumed to remedy the severity of this viral infection is the oil of Nigella sativa seed commonly called black seed oil. In this study, we extracted and characterized the oil from this seed using gas chromatography coupled to a mass selective detector to identify the component phytochemicals. Site-directed multiligand docking of the identified compounds was performed on SARS-CoV-2 molecular targets- Replicase polyprotein 1a, RNA binding protein of NSP9, ADP ribose phosphatase of NSP3, 3-chymotrypsin-like protease 3CLpro, and RNA-dependent RNA polymerase RDRP, and ACE2–angiotensin-converting enzyme from the Homo sapiens. Results The binding affinity of caryophyllene oxide was the highest on 3CLpro (− 6.0 kcal/mol), NSP3 (− 6.3 kcal/mol), NSP9 (− 6.3 kcal/mol), and RDRP (− 6.9 kcal/mol) targets, while α-bergamotene gave the best binding affinity on RPIA (5.7 kcal/mol) target. The binding affinity of β-bisabolene on the ACE2 target (− 8.0 kcal/mol) was almost the same as Remdesivir (− 8.1 kcal/mol). The ADMET properties of these three phytochemicals showed that they are good drug leads for these SARS-CoV-2 receptors. Conclusion The findings from this study strongly indicate that the reported recovery from COVID-19 infection claimed by patients who consumed black seed oil could be linked to the presence of caryophyllene oxide, α-bergamotene, and β-bisabolene in this natural product.


Background
Black cumin, scientifically known as Nigella sativa L., is an annual dicotyledonous herb that belongs to the family Ranunculaceae (Datta et al. 2012). It is a small shrub ranging between 20 and 90 cm in height, with tapering green leaves and rosaceous pink, purple, white, yellow, or pale blue flowers with 5-10 petals. Flowering starts in July and the seeds ripen in September. The ripe fruit is a capsule with 3-7 united follicles containing numerous tiny, pitch-black seeds with an aromatic and slightly bitter taste (Ahmad et al. 2013). It is native to northeastern Africa, the eastern Mediterranean, and southwestern Asia (Vavilov and Dorofeev 1992). The species is grown in many countries like India, Turkey, Middle East, Pakistan, Syria, and Saudi Arabia (Naz 2011). The seeds of N. sativa have been widely used as diuretics, anti-diarrheal, analgesic, antibacterial, antidiabetic, anticancer, immunomodulator, anti-inflammatory, spasmolytic, bronchodilator, hepatoprotective, renal protective, gastro-protective, antioxidant, etc. (Ahmed et al. 2013;Sultan et al. 2014).
Page 2 of 13 Duru et al. Bull Natl Res Cent (2021) 45:57 The chemical composition of seed extracts of N. sativa seed varies from one country to another and also from one region in a country to another (Shariq et al. 2015). The major components in N. sativa seed are fatty oils and aliphatic hydrocarbons, arachidonic acid, tocopherols, γ-linolenic acid, and essential oil, which contains carvacrol, ρ-cymene, thymoquinone, α-pinene, and indazole-type alkaloids, isoquinoline alkaloids, as well as dolabellane-type diterpene alkaloids. The food and therapeutic applications of N. sativa oil and seed have a long history in many cultures. Islamic medicine also cataloged the numerous benefits of this plant (Sharma et al. 2009). Its compounds have been shown to have various pharmacological effects on different body parts. Empirical studies have reported that the seed oil and seed extracts have antimicrobial activities, which some workers have attributed to the presence of compounds such as thymohydroquinone and melanin (Kalim 2013;Al Yahya 1986). The seed extract contains numerous antioxidants like carvacrol, t-anethole, and 4-terpineol, and their activities have been reported in different studies (Awadalla 2012). Also, available reports have shown the anti-inflammatory, anticancer, anti-diabetic, and hepatoprotective activities of its oil and extracts (Umar et al. 2012;Peng et al. 2013;Abdelmeguid et al. 2010).
The Coronavirus disease 2019 (COVID-19) caused by the severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2) began in Wuhan, China, and was reported to the World Health Organization on December 31, 2019 (WHO 2020). The pandemic ravaged the entire globe, claiming hundreds of thousands of lives. Community to community transmission is now evident in many countries, with a daily rise in the number of cases. There are potential drug targets essential for the survival of the coronavirus, which include a 3-chymotrypsin-like protease (3CLpro), Angiotensin-converting enzyme (ACE2), ADP ribose phosphatase of Nonstructured Protein 3 (NSP3), Nonstructured Protein 9 RNA binding protein (NSP9), RNA-dependent RNA polymerase (RdRp), and Replicase Polyprotein 1a (RP1a). To infect the host, the structural spike glycoprotein interacts with the transmembrane protein of the human host cell receptor ACE2 (Guy et al. 2005). The main protease 3CLpro acts in the viral polyprotein large Replicase Polyprotein 1a (RP1a) proteolytic cleavage, cleaving it into functional units to produce nonstructural viral constituents like RdRp, NSP3, and NSP9. This makes 3CLpro and ACE2 important drug targets to inhibit Nonstructured Proteins production and terminate viral action. The RP1a and the NSPs (RdRp, NSP3, and NSP9) are important targets involved in the replication and translation process, which lead to a viral proliferation in host cells (Kumar et al. 2020). In vitro and in silico studies of possible cures for this dreaded disease have become hot topics in current research. Although no drug is currently available to cure this disease, Remdesivir, an antiviral drug has been clinically tested and confirmed useable in the interim (Williamson et al. 2020).
Recently, the governor of Oyo State in Nigeria revealed that the oil from N. sativa seed was one of the remedies he consumed to combat the disease, a week after he tested positive for COVID-19 (Punch Newspaper, April 7, 2020). So far, no attempt has been made to investigate the veracity of this claim. In this study, the oil from the seed of N. sativa was extracted, and the identified phytochemical components docked on some selected SARS-CoV-2 proteins. The dock scores of these compounds were compared to the values obtained with a standard COVID-19 prescription drug.

Sample collection and preparation
Mature seeds of N. sativa were collected toward the end of rainy season in October 2020 (average temperature 29 °C), from a fallow farmland in Northern Nigeria (N 10° 9′ 32.2524", E 8° 8′ 1.8816"), and no permissions or licenses were required for the collection and study. They were identified by a professional taxonomist Prof. F.N. Mbagwu of the Department of Plant Science and Biotechnology, Imo State University, Owerri, Nigeria, and the voucher specimen was deposited at the Imo State University Herbarium, with herbarium number IMSUH-468. They were washed with tap water and left to dry on a laminated paper for one week under room conditions. The dried seeds were subjected to preliminary grinding using mortar and pestle, and the particle size was further reduced using an electric blender. The seed powder was transferred in an airtight plastic container and stored in a refrigerator at 4 °C.

Extraction of oil
Oil was extracted from the ground seed powder using a Soxhlet apparatus mounted on a heating mantle (Duru 2020). The seed powder of weight 100 g was encapsulated in a clean white cotton cloth previously degreased in hexane (≥ 98.5%, Sigma-Aldrich, Darmstadt, Germany). The sample was held in the thimble-holder of the extractor which was gradually filled with condensed redistilled hexane from a distillation flask placed on a heating mantle. When the solvent reaches the overflow level in the thimble-holder, it is siphoned back into the distillation flask, carrying the extracted oil with it. The operation was repeated until complete extraction was achieved. The extracted oil was then recovered from the solvent by rotary evaporation.
Page 3 of 13 Duru et al. Bull Natl Res Cent (2021) 45:57 Characterization of oil The oil sample was characterized using gas chromatography-mass Spectrometer (GC-MS) instrument (Model: 7890 GC and 5977B MSD, Agilent Technologies, USA). An HP-5 MS capillary standard nonpolar column L × I.D., 30 m × 0.25 mm and film thickness 0.25 µm was used. The flow rate of mobile phase (carrier gas: He) was set at 1.0 mL/min. In the gas chromatography part, the temperature program (oven temperature) was set at 40 °C and raised to 250 °C at 5 °C/min and injection volume 1 µL. The oil sample was dissolved in methanol (≥ 99.8%, Sigma-Aldrich, Darmstadt, Germany) and fully scanned at a range of 40-650 m/z, and the results were compared using NIST mass spectral library search program (El-Sawi et al. 2020).

Identification and preparation of ligands
The 3D structure-data files (SDF) of the compounds in the crude oil sample were identified and downloaded from the PubChem database. They were minimized in PyRx virtual screening tool, using Universal Force Field at 200 steps. They were then converted to AutoDock ligands (pdbqt) and used for the docking analysis.

Docking procedure and analysis of results
The screening of the phytochemical compounds from the oil was performed by docking them on selected binding pockets of SARS-CoV-2 proteins and ranked based on their binding energies. The multiple docking of the ligands and proteins was done with Autodock Vina in PyRx software (Tsao et al. 2020;Johnson et al. 2020). The molecular docking results were organized on an Excel spreadsheet, and the Heat Map of the data was viewed using the Conditional Formatting function.

Absorption, distribution, metabolism and elimination and toxicity (ADMET) analysis
The compounds with the lowest binding energy for each protein were selected and submitted to the ADMETsar 2 server to examine their drug-like properties, pharmacokinetics, and pharmacodynamics parameters (Yang et al. 2019).

Analysis of protein-ligand interactions
The interactions between the most potent compounds in the oil and amino acid residues of the proteins were visualized using Biovia Discovery studio client 20.1 (BIOVIA 2020) and UCSF Chimera software.

Chemical composition of oil extract
The GC-MS analysis of N. sativa seed oil gave 27 peaks for the compounds detected, and the chromatogram is shown in Fig. 2.
The elution information of the compounds and their therapeutic uses are summarized in Table 1.
The binding affinities of the compounds from the oil extract to the SARS-CoV-2 receptors are shown in Table 2.

Key: Good Fair Poor
Page 7 of 13 Duru et al. Bull Natl Res Cent (2021) 45:57 characteristics, molecule size, and flexibility, as well as the presence of various pharmacophores that influence the behavior of a molecule in a living organism, were used to predict the drug likeliness of the compounds. A good drug candidate should not violate more than one of the rules (Lipinski 2016). The molecular weights of the compounds were < 500, and their hydrophobicity (log P) did not exceed 5. The hydrogen bond donor ( ≤ 5 hydrogens) and hydrogen bond acceptor (not more than 10 hydrogens) of the compounds were in line with the rule. The rotatable bonds (not more than 3) were in line with the rule of three, and their acute oral toxicity values were well below 5 mg/kg. Caryophyllene oxide is soluble in water, while α-bergamotene and β-bisabolene are moderately soluble. The human intestinal adsorption and blood-brain barrier of the compounds were very high, and the compounds were not carcinogenic. These results are indications that α-bergamotene, β-bisabolene, and caryophyllene oxide are excellent drug candidates (Zhong 2017).
The 3D and 2D protein-ligand interaction images for α-bergamotene, β-bisabolene, and caryophyllene oxide are shown in Fig. 5. The amino acid residues involved in the interactions with the compounds and the control drug Remdesivir are depicted in Table 4.
The hydrogen bond, pi-sigma, pi-alkyl, and alkyl type interactions were the main forces that held the compounds in the protein pockets. At least two hydrogen bonds were observed in all the protein-Remdesivir interactions, and one was found only in the RDRP-caryophyllene oxide interaction. The protein-Remdesivir interactions in all the protein pockets had at least two hydrogen bonds and the other interacting forces, in the pockets of the studied proteins. Caryophyllene oxide and Remdesivir interacted with ASN 209 by hydrogen bonding in RDRP, while α-bergamotene and Remdesivir interacted by alkyl forces with ILE 120, ALA 150, and ALA 102. These suggested that the mechanism of inhibition of the proteins at these sites by caryophyllene oxide and α-bergamotene was similar to the action of Remdesivir at those sites.

Conclusions
The phytochemical composition and binding affinity of compounds in the oil extract from Nigella sativa seed on SARS-CoV-2 molecular targets were investigated. The  GC-MS analysis of the oil showed that it contained many compounds with known biological and therapeutic activities. The docking studies of these compounds against some major targets on the SARS-CoV-2 revealed that caryophyllene oxide, α-bergamotene, and β-bisabolene were promising candidates in the inhibition of the activity of these proteins. Their ADMET results further indicated that they would have good pharmacological and pharmacodynamic properties in the human system. There is therefore need for further in vitro, in vivo, and clinical trials to validate the potentials of these compounds in the eradication of the coronavirus pandemic.