key: cord-0969218-eil3gpx4
authors: Falade, Victoria Adeola; Adelusi, Temitope Isaac; Adedotun, Ibrahim Olaide; Abdul-Hammed, Misbaudeen; Lawal, Teslim Alabi; Agboluaje, Saheed Alabi
title: In silico investigation of saponins and tannins as potential inhibitors of SARS-CoV-2 main protease (M(pro))
date: 2021-01-06
journal: In Silico Pharmacol
DOI: 10.1007/s40203-020-00071-w
sha: a10ff4a048e063c82a8ad229a2242c6af719756b
doc_id: 969218
cord_uid: eil3gpx4

It is no longer news that a novel strain of coronavirus named SARS-CoV-2 is ravaging the health sector worldwide, several attempts have been made to curtail this pandemic via repurposing of old drugs but at the present, available drugs are not adequately effective. Over the years, plant phytochemicals are increasingly becoming alternative sources of antimicrobial agents with novel mechanisms of action and limited side effects compared to synthetic drugs. Isolated saponins and tannins were evaluated for antiviral activity against SARS-CoV-2 (M(pro)) via Molecular Docking and it was observed that a handsome number of the phytochemicals had binding affinities much better than Remdesivir, Dexamethasone, and N3 inhibitor which were used as the standards in this study. Further investigation of drug-likeness, ADMET profile, PASS profile, oral bioavailability, bioactivity, binding mode, and molecular interactions of these phytochemicals revealed that binding affinity alone is not enough to justify the potency of a molecule in the drug discovery process, as only 4 among the screened compounds passed all the analyses and are identified as potential inhibitors of SARS-CoV-2 (M(pro)). This preliminary study thereby recommends Ellagic acid (− 8.4 kcal/mol), Arjunic Acid (− 8.1 kcal/mol), Theasapogenol B (− 8.1 kcal/mol), and Euscaphic Acid (− 8.0 kcal/mol) as potential inhibitors of SARS-CoV-2 (M(pro)) with better pharmacokinetics and bioavailability compared to Remdesivir which is currently used compassionately. [Image: see text]

The outbreak of a novel strain of coronavirus has been an unprecedented challenge in the healthcare system globally. This was initially identified in a cluster of Pneumonia patients in Wuhan, Hubei province, China in December 2019 and had since spread to other continents of the world Chang et al. 2020 ). The causative virus was found to be closely related to a previously known Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV), hence the name SARS- CoV-2 (Wand et al. 2020) . Clinical symptoms include pneumonia, cold, dry cough, sore throat, fever, shortness of breath, respiratory tract infection, headache, diarrhea, and abdominal pain (Chen et al. 2020) . Over twenty (20) million symptomatic and asymptomatic cases had been reported, with more than a million deaths (https :// coron aviru s.jhu.edu/map.html). The infection was reported to be severe in elderly patients and more critical in individuals with underlying health conditions like diabetics, cardiovascular disease, obesity, pulmonary fibrosis, asthma, cancer, etc (Wu et al. 2020) .

As of the time of this report, no drug/vaccine has been proven to effectively cure/manage this infection. The current strategy in the management includes early diagnosis, isolation of infected individuals, and maintenance of personal hygiene, physical distance, and use of personal protective equipment. Attempts have also been made to reposition/ repurpose old drugs of promising antiviral potential in a bid to tackle this fast-spreading pandemic. Remdesivir, lopinavir-ritonavir, favipiravir, etc. are among the drugs used in the management of this pandemic (Wu et al. 2020 ). More recently dexamethasone has been found to reduce mortality in severe cases of COVID-19 disease (Johnson and Vinetz 2020) .

Phytochemicals are increasingly becoming the alternative therapeutic and pharmacological agent of great importance in drug discovery and development. They have limited or no side effects on administration, possess a novel mechanism of action, and great chemical diversity which enhances their therapeutic interaction with varied biological targets compared to synthetic drugs. Saponins are naturally occurring non-volatile and surface-active glycosides of triterpenes and steroids with a wide range of pharmacological activities, including anti-inflammatory, anticancer, antiviral, antitumor, antifungal, hypoglycemic, and cytotoxic activity. They are potential vaccine adjuvant due to their ability to stimulate and activate the immune system in mammals (Sun et al. 2009; Skene and Sutton 2006) . The sugar side chain in saponins could be responsible for their adjuvant. Similarly, saponins have been found to possess a novel mechanism of action on the virus which involves the destruction of viral envelop, loss of binding sites, and subsequent prevention of binding of the virus to host cell (Roner et al. 2007) . They have also found application in food, confectionery, beverages, cosmetics, preparation of pharmaceutical products, as well as industrial processes as surface-active and foaming agents. Rich sources of saponins include soya beans, capsicum peppers, spinach, horse chestnut, fenugreek, ginseng, etc. (Francis et al. 2002; Vincken 2007) . Moreover, much toxicity is observed when saponins are administered intravenously while lower toxicity is observed when is administered orally due to low absorption in the alimentary canal but it is rather hydrolyzed to sapogenins via enzymatic action (Guclu-Ustundag and Mazza 2007) .

Tannins on the other hand are phenolic compounds of great antioxidant potential and antiviral properties, they are high molecular weight compounds (ranging between 500 and 3000 gmol-1) found in different parts of plants including leaf, stem, root, or bark (Serrano et al. 2009 ). This study aims at investigating the antiviral potential of saponins and tannins against SARS-CoV-2 M pro through molecular docking coupled with ADMET-studies, pharmacokinetic evaluation, drug likeliness among other analyses at a therapeutic dose.

In this research, the inhibitory activities of twenty-seven (27) phytochemicals from Saponins and Tannins against SARS-CoV-2 main protease (M pro ) PDB ID: 6LU7 were studied. The selected isolated phytochemicals from different species of plants have been reported to possess broad-spectrum biological activities including antiviral, anti-inflammatory, cytotoxic activities, anti-diabetic, neuroprotective, immunomodulatory, antioxidant, and anticancer, to mention a few. Their reported biological activities arouse our interest to investigate them as potential inhibitors of the novel SARS-CoV-2 main protease (M pro ) towards designing a new therapeutic agent. The sources, species name, and biological activities of the studied compounds are shown in Table 1 . Also, three standards (Remdesivir, Dexamethasone, and N3 inhibitor) whose completed randomized trials and inhibitory activities against the novel SARS-CoV-2 (M pro ) have been reported (Kaddoura et al. 2020; Jin et al. 2020) were used for comparison. The Saponins/Tannins used are Priverogenin A, Arjunic acid, Theasapogenol B, Euscaphic Acid, Camelliagenin C, Medicagenic Acid, Protoescigenin, Arjunolic acid, Asiatic acid, Protobassic acid, Arjugenin, Polygalacic acid, Primulagenin A, Soyasapogenol B, Tomentosic acid, Presenegenin, Punicalagin, Punicalin, Ellagic acid, Corilagin, Gallagic acid, Terflavin B. Catechin, Chebulinic acid, Hexahydroxydiphenic acid, Gallic acid, and Catechol respectively.

Using conformer distribution with Molecular mechanics/ MMFF in Spartan 14 version 1.1.4, the most stable conformers of all the studied compounds were obtained. The conformers were optimized using equilibrium geometry density functional theory method (DFT) at B3LYP and 6 − 31 + G* as the basis set on HP Desktop computer, 2 terabytes hard disk, 64-gigabyte random accessed memory (RAM), Intel ® Core™i7-2600 CPU, 3.40 GHz processor, and 4 gigabytes dedicated video memory to generate important molecular properties and well optimized structures for molecular docking simulation.

The crystal structure of SARS-CoV-2 M pro (PDB ID: 6LU7) ( Fig. 1 ) was used as a target receptor in this study. The structure was retrieved from the protein data bank (RCSB) (http:// www.rcsb.org/pdb). SARS-CoV-2 M pro is the main protease of the novel strain of the 2019 coronavirus disease. It is a protein that mediates replication and transcription of the virus (Jin et al. 2020) , therefore, the protease is often the target of the potential inhibitors of the virus in the drug discovery and development processes. The resolution of 6LU7 (2.16 Å) approximately falls within the 2.0 Å recommended resolution for a protein of good quality (Hajduk and Tse 2005) . Water molecules and other unwanted complexes were removed from the downloaded protease to avoid undesired molecular interactions and to ensure that no molecule interfered with the potential binding site of the target protease during the docking simulation. The Ramachandran plot ( Fig. 2) which revealed the quality of the receptor under study was obtained using the Volume, Area, Dihedral Angle Reporter (VADAR) webserver. The binding pocket X, Y, and Z coordinates were defined as − 26.284, 12.603, and 58.96 respectively for the screening exercise using the grid box of the native ligand inhibitor (N3 inhibitor) complexed with the target receptor as a basis.

Binding pocket, ligand interactions, and all amino acids in the active site of SARS-CoV-2 M pro were established using the Computed Atlas for Surface Topography of Proteins (CASTp) (http://sts.bioe.uic.edu/castp /index .html?2011) (Tian et al. 2019) and Biovia Discovery Studio (2019). The obtained data were compared and validated with the previously reported experimental data for the SARS-CoV-2 M pro active site complexed with N3 native ligand (Jin et al. 2020) . 

AutoDock/Vina tool (Trott and Olson 2010) was used for the docking simulation. It is a reliable protein-ligand docking tool that uses the Broyden-Fletcher-Goldfarb-Shanno algorithm which significantly improves the average accuracy of the binding mode prediction (Azam and Abbasi 2013). AutoDock tool was used to generate the (.pdbqt) files using the grid size 58, 62, and 40 for x, y, and z axes, grid center − 26.284 × 12.603 × 58.960 respectively, and 1.000 Å spacing. The docking simulation was performed using AutoDock

Vina and the binding energies (∆G) (kcal/mol) of the docked ligands were obtained. The inhibition constants and inhibitory efficiencies of the docked ligands were calculated using (Eqs. 1 and 2). PyMol and Biovia-2019 discovery studio were used to view and analyse the docking results.

where R = Gas constant (1.987 × 10 -3 kcal/mol); T = 298.15 K (absolute temperature); K i = Inhibition constant; ΔG = Binding energy.

Prediction of Activity Spectra for Substances (PASS) web server http://www.pharm aexpe rt.ru/passo nline /) (Lagunin 2000) was used to predict the biological activities of the studied ligands.

Absorption, distribution, metabolism, excretion and toxicity (ADMET) properties of the selected compounds were predicted using ADMET SAR2 web-server, while drug-like 

The X-ray crystallographic structure of SARS-CoV-2 M pro (PDB ID: 6LU7) (Fig. 1) contains 306 amino acid residues complexed with an inhibitor (N3-(N-

. It consists of 23%, 31%, 45% and 28% α-helix, β-sheets, Coil and Turns respectively. The resolution of the protease as revealed by X-ray diffraction was 2.16 Å, crystal dimension is a = 97.93 Å, b = 79.48 Å, and c = 51.08 Å with angles α (900), β (114.550), and γ (900) SARS-CoV-2 M pro active site is located in the cleft between Domain I and II and contains a Cys-His catalytic dyad. Amino acid residue at the active site is as follows Thr24, Thr25, Thr26, His41, Met49, Tyr54, Phe140, Leu141, Asn142, Gly143, Ser144, Cys145, His163, His164, Met165, Glu166, Leu176, Pro168, His172, Asp187, Arg188, Gln189, Thr190, Ala191, and Gln192 (Jin et al. 2020 ).

Recent developments in drug discovery have led to a renewed interest in the computational study which involves the use of algorithms and programs for predictions of therapeutic interventions in biological processes (Shoichet et al. 2002) . Molecular Docking is a structure-based drug design approach that predicts binding interactions between ligand and target receptor at the binding site (Ferreira et al. 2015) . It is an important virtual screening technique that could screen several thousands of ligands against the target as well as identify potential inhibitors of the target receptor with speed and accuracy (Dias et al. 2008) .

To investigate potential inhibitors of SARS-CoV-2 M pro , AutoDock/Vina (MGL tools-1.5.6), PyMOL Console Edu, and Biovia Discovery studio 4.5 were used. The results obtained from the docking of selected saponins and tannins against SARS-CoV-2 M pro were as shown in Table 2 . Binding affinity is a reflection of the inhibitory activity of the plant extract against SARS-CoV-2 M pro . It is apparent from this table that most of the selected saponins and tannins had better binding affinities and inhibitory activities against SARS-CoV-2 M pro compared to Remdesivir and Dexamethasone which are standard drugs used in the management of this pandemic disease. Binding affinity (BA) for the selected saponins range between − 8.3 k cal/mol and − 7.1 kcal/mol, while those of tannins range between − 9.0 kcal/mol and − 4.7 kcal/mol. It can be seen that selected tannins had far greater binding affinities compared to saponins, with Punicalagin (− 9.0 kcal/ mol) having the outstanding inhibitory activity. To validate the inhibitory effects observed through docking scores, the hit compounds among the docked ligands were subjected for further analysis.

Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) profile of a molecule is an important assay in the early stage of drug discovery. ADMET data enhances the selection and identification of molecules with optimum safety profile at a therapeutic dose along the discovery process rather than at the final stage, as this helps in avoiding waste of time and precious resources on drug molecules that may eventually be discarded (Tsaioun and Kates 2010) .

Considering the binding affinities and inhibition constants ( Table 2) which are expected to be within (0.1 µM and 1.0 µM), only nine (9) of the docked compounds qualified as hits and were subjected to ADMET analysis using ADMET SAR-2 web-server (Cheng et al. 2012) . Notably, only Ellagic acid, Arjunic Acid, Thesapogenol B, and Euscaphic Acid (coded C-1, C-2, C-3, and C-4 respectively) showed excellent ADMET profile, their structures were shown in Fig. 3 , their ADMET properties were discussed in Table 3 and they were selected for further analyses. Although Punicalagin, Punicalin, and Priverogenin A had better inhibitory activities and binding affinities than ellagic acid, they were shunned on the basis of toxicity and the possibility of poor absorption and permeability across a biological membrane. As part of the drug ADMET profile, a drug molecule should have good human intestinal absorption (HIA), solubility (Log S) range between -1 and -5, should be a noninhibitor of cytochrome P450 enzymes, and should be non-Ames toxic. Others include non-carcinogenicity, noninhibition of hERG, and no or low level of toxicity (Tsaioun and Kates 2010). C-1, C-2, C-3, and C-4 (Table 3) were well absorbed in the human intestine, only Thesapogenol B (C-3) was found to cross the blood-brain barrier, although an oral drug does not necessarily need to cross the blood-brain barrier, only central nervous system target drug need to (Hughes et al. 2011) .

It was also found in the prediction, that all the four (4) selected hit compounds were non-inhibitor of microsomal enzymes (Cytochrome P450), which is an indication of good metabolism of the drug in the liver (Stevens 2014). The potential of a drug molecule to cause mutation in DNA is revealed by Ames toxicity value and could be a major reason for excluding a drug molecule along the discovery process (Mccarren et al. 2011) . As shown in Table 3 , all the selected hit compounds are non-genotoxic and non-carcinogenic. Similarly, the hits possess type III acute oral toxicity values (slightly toxic) which could easily be converted to type IV (nontoxic) during hit-lead optimization. The human ether-a-go-go related gene (hERG) potassium ion channel plays important role in cardiac repolarization, blockage of which may be caused by inherited mutation or some drug molecules, leading to long QT syndrome and eventual death (Sanguinetti and Tristani-firouzi 2006) . Interestingly all the hit compounds are non-blockers of the hERG potassium channel.

Comparison of ADMET profile of the hits with the standard drugs revealed that Remdesivir, Dexamethasone, and N3 inhibitor had good ADMET profiles, but both dexamethasone and N3 inhibitor were found to be blockers of hERG potassium ion channel while Remdesivir is not. This result, therefore, suggests that the hit compounds have a similar ADMET profile just like Remdesivir, and could be developed further in the quest of finding a new therapeutic agent in COVID-19 management. It is important to note that the ADMET profile of dexamethasone revealed that it is a blocker of hERG potassium ion channel, therefore caution should be taken in its administration to COVID-19 patients. 

Evaluation of physicochemical and drug-likeness of the potential active compound is an important step in drug discovery. As proposed by Lipinski, an effective oral therapeutic drug must obey the 'rule of five' with not more than one (1) violation, this is due to the fact that an orally bioavailable drug must possess molecular weight (MW) ≤ 500 Da, hydrogen bond donor (HBDs) ≤ 5, hydrogen bond acceptor (HBAs) ≤ 10 and log P (octanol-water partition coefficient) ≤ 5 (Lipinski 2004) . These descriptors of oral bioavailability are important as they predict the permeability and absorption of such drug across a biological membrane such as an epithelium cell, partition coefficient value (log p) is especially important in predicting intestinal absorption of such drug (Aucamp et al. 2015) . Drug-likeness of the selected hits were evaluated with Molinspiration online (http://www.molin spira tion.com/) as shown in Table 4 , it is apparent from the table that none of the selected hits had more than one violation of the 'rule of five' which is an indication of good oral bioavailability and permeability. It is also interesting to note that all the selected hits and Dexamethasone (SD-2) had better drug-like properties compared to Remdesivir, and N3-Inhibitor (SD-1 and SD-3) with more than 1 violation. This result revealed that the selected hits (C-1 to C-4) possessed excellent drug-like properties and could be developed further as oral drugs.

Oral-bioavailability of the selected hit compounds and standards Table 5 shows the oral-bioavailability profile of the selected hit ligands and standards obtained using the swissADME tool (Daina et al. 2017 ). The bioavailability radar (Fig. 4) revealed Table 4 Drug-likeness of selected hit compounds C-1 = Ellagic Acid, C-2 = Arjunic Acid, C-3 = Thesapogenol B, C-4 = Euscaphic Acid, SD-1 = Remdesivir, SD-2 = Dexamethasone, SD-3 = N3 Inhibitor Table 5 Oral bioavailability of the selected hit compounds C-1 = Ellagic, C-2 = Arjunic acid, C-3 = Theasapogenol B, C-4 = Euscaphic acid SD-1 = Remedesivir, SD-2 = Dexamethasone, SD-3 = N3 Inhibitor LIGAND C-1 C-2 C-3 C-4 SD-1 SD-2 SD-3 the oral-bioavailability profile of the selected hits and Standard drug in a catch sight. The pink area of the radar shows the optimum zone for each of the properties i.e., (POLAR, FLEX, SIZE LIPO, INSATU, and INSOLU). As shown in Table 5 , all the selected hits obeyed the recommended size of 500 gmol-1 for effective drug candidates as reported by Lipinski, compared to 602.6 gmol-1 and 680.79 gmol-1 reported for SD-1 and SD-3 respectively. The polarity (POLAR) of the selected hits were evaluated using their Total Polarity Surface Area (TPSA). A polar compound is expected to have TPSA value within the range of 20 and 130 Å2. The TPSA of C-2 and C-3 fall within the recommended range, while all the hit compounds have lower TPSA values compared to the value obtained for SD-1 and SD-3. However, SD-2 has the best TPSA value (94Å2), thus has the highest POLARITY property. The INSOLU (insolubility) requirement of the selected compounds and standards as depicted in their ESOL (Log S) and ESOL Class showed that C-1 and SD-2 are soluble; C-2 has a poor solubility profile while C-3, C-4, SD-1, and SD-3 are moderately soluble. Notably, C-1 has the most outstanding Table 6 shows the bioactivity properties of the four (4) selected hits and standards. The inverse relationship existing between binding energy and inhibition constant is as shown in Eqs. 1 and 2, indicating that the higher the binding energy the lower the inhibition constant. A potential hit compound is expected to have inhibition constant values ranging between 0.1 and 1.0 µM and not more than 10 nM for a drug (Stevens 2014) .

Inhibition constant values of the hit compounds range from 0.70 to 1.37 µM. This observation revealed that all the four (4) selected compounds are qualified as hit with Ellagic acid (0.70 µM) being the most potent of all. Ligand Efficiency (LE), Fit Quality (FQ), and Ligand-efficiency-dependent lipophilicity (LELP) were also calculated according to Eq. 3-6. C-1 has the ligand efficiency (LE) value within the recommended value of (≥ 0.3), while both C-1 and SD-2 have LELP values within the recommended range of (− 10 and 10). Notably, all the hit compounds had FQ values within the recommended range (≥ 0.8), which is higher than the FQ values reported for all the standards (Hopkins et al. 2014 ). Table 7 shows the prediction of the biological activity spectra of the selected hits and standards using the PASS online web tool (Filimonov et al. 2014) . It shows predictions based on the structure-activity relationship (SAR). Pa is the probability that shows the similarity in structure of specific molecules in comparison to the molecules of the most active compound of similar activity. For this activity spectrum, the probability of validating the biological activity of a molecule (Filimonov et al. 2014) . It can be observed from Table 7 that all the selected hits (Ellagic acid, Arjunic acid, Theasapogenol, and Euscaphic acid) and SD-1 (Remdesivir) possess antiviral activities while (SD-2) Dexamethasone showed no antiviral activity within the spectrum. Also, Table 7 revealed that SD-3 had higher activity spectra for severe acute respiratory syndrome treatment and Protease Inhibition potential with Pa values of 0.698 and 0.665 respectively, which is an indication of its tendency to inhibit the SARS CoV-2 main protease. Notably, The Pa values observed for the hit compounds were greater than Pi, thus validating their biological activities. Therefore, all the selected hits possessed excellent antiviral activities and could be explored in the design and development of new drugs/vaccines in the management of the novel SARS-CoV-2 ravaging the globe.

Identifying the binding site and evaluating the binding interactions of the ligands in the active pocket of the target receptor plays a crucial role in drug design. This eases the improvement of ligand affinity to the pocket during the lead optimization stage of drug discovery. The binding interactions of the selected hits and standards with SARS-CoV-2 main protease (M pro ) are shown in Tables 8 and 9 . The binding affinities obtained from the docking of C-1, C-2, C-3, C-4, SD-1, and SD-2 as shown in Table 8 are (− 8.4, − 8.1, − 8.1, − 8.0, − 7.6, and − 7.7 kcal/mol) respectively. C-1 (Table 9 ) formed a conventional hydrogen bond with Gln189, Leu141, Gly143, Ser144, His163, and Asn142, carbon-hydrogen bond with His172, pi-donor hydrogen bond with Met165, pi-alkyl interaction with Cys145. However, a close examination of the main active site of the target receptor (6LU7) revealed that its binding pocket (active site) was located in the cleft between domains I and II with 

This study evaluates twenty-seven (27) bioactive saponins/ tannins against SARS-CoV-2 main protease (M pro ) via in silico studies (structure-based drug design). The compounds were screened using the reliable AutoDock/Vina tool, ADMET SAR-2, swissADME, Molinspiration web server, PASS software, among others. CASTp and VADAR web tools were used to establish the active site and Ramachandran plot of the target receptor (6LU7) respectively. The results obtained flourish Ellagic acid (− 8.4 kcal/mol), Arjunic acid (− 8.1 kcal/mol), Theasapogenol B (− 8.1 kcal/mol), and Euscaphic Acid (− 8.0 kcal/mol) as probable inhibitors of SARS-CoV-2 Mpro due to their excellent binding energies, ADMET profile, drug-likeness, oral-bioavailability properties, PASS properties, Bioactivity, outstanding binding mode and molecular interactions with the target receptor. C-1 (Ellagic acid) share the same pocket and interact completely with the active site of the target receptor, while C-2, C-3, and C-4 (Arjunic acid, Theasapogenol B, Euscaphic Table 9 Binding mode and molecular interactions of the selected hits against 6LU7

C-1 = Ellagic acid; C-2 = Arjunic acid; C-3 = Theasapogenol B; C-4 = Euscaphic acid Acid) respectively share the same probable allosteric pocket and interactions with SD-1 (Remdesivir), and SD-2 (Dexamethasone) whose randomized clinical trials against the target receptor have been completed. Although, it is widely accepted that computational drug-design plays a pivotal role in the world of drug design and development today. It is used in hit/lead identification and optimization towards the design and development of new therapeutic agents, however, the need for experimental work to buttress our findings is duly acknowledged, but lack of funding among other factors limited our scope. We hereby recommend that the fourhit compounds identified in this study should be developed further towards finding a reliable and effective drug/vaccine to arrest the new disease (COVID-19) terrorizing the entire human race.

Structural and Active site analysis of SARS-CoV-2 Mpro complexed with N3 inhibitor (PDB ID: 6LU7)"). This implies that C-1 shares the same pocket and interacts effectively with the active site of the target receptor and could be its potential inhibitor. Furthermore, strong similarities were observed in the interactions of C-2, C-3, and C-4 with 6LU7 (Table 9) as compared with standards SD-1 and SD-2 (Table 8). C-2 formed a conventional hydrogen bond with Leu287, and alkyl interaction with Leu286

Structure of coronavirus main proteinase reveals combination of a chymotrypsin fold with an extra alpha-helical domain

Amorphous azithromycin with improved aqueous solubility and intestinal membrane permeability

Molecular docking studies for the identification of novel melatoninergic inhibitors for acetylserotonin-O-methyltransferase using different docking routines

Coronavirus disease 2019: coronaviruses and blood safety

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study

admetSAR: a comprehensive source and free tool for assessment of chemical ADMET properties

SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules

Molecular docking algorithms

Molecular docking and structure-based drug design strategies

Prediction of the biological activity spectra of organic compounds using the PASS online web resource. Chem Heterocyclic Compounds

The biological action of saponins in animal systems: a review

Saponins: properties, applications and processing

Predicting protein druggability

The role of ligand efficiency metrics in drug discovery

Principles of early drug

Structure of Mpro from COVID-19 virus and discovery of its inhibitors

Dexamethasone in the management of covid -19

COVID-19 therapeutic options under investigation

PASS: prediction of activity spectra for biologically active subtances

Lead profiling Lead-and drug-like compounds: the rule-of-five revolution

An investigation into pharmaceutically relevant mutagenicity data and the influence on Ames predictive potential

Antiviral activity obtained from aqueous extracts of the Chilean soapbark tree (Quillaja saponaria Molina)

hERG potassium channels and cardiac arrhythmia

Review Tannins: current knowledge of food sources, intake, bioavailability and biological effects

Lead discovery using molecular docking

Saponin-adjuvanted particulate vaccines for clinical use

Lead Discovery. In: Jaworski A (ed) Medicinal chemistry: modern drug discovery process

Advances in saponin-based adjuvants

CASTp 3.0: computed atlas of surface topography of proteins

Software news and update AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading

ADMET for medicinal chemists: a practical guide

Saponins, classification and occurrence in the plant kingdom

Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro

An update on current therapeutic drugs treating COVID-19

A pneumonia outbreak associated with a new coronavirus of probable bat origin

The authors wish to acknowledge members of the Computational and Biophysical Chemistry Research Group in the Department of Pure and Applied Chemistry, Ladoke Akintola University of Technology, LAUTECH, Ogbomoso, Oyo State, Nigeria.Author contributions MA-H, ATI: conceptualization, methodology, validation, supervision, writing-review and editing, visualization, data curation, resources. VAF and IOA: conceptualization, writing-original draft, writing-review and editing, methodology, investigation, data curation, visualization. TAL and SAA: conceptualization, methodology, data curation, investigation, visualization.Funding The author(s) received no financial aid or support in any form for this research work.

Conflict of interest The authors declare that there is no conflict of interest in this manuscript.