key: cord-0788832-17rj71x4
authors: Dutta, Mycal; Nezam, Mohammad; Chowdhury, Subrata; Rakib, Ahmed; Paul, Arkajyoti; Sami, Saad Ahmed; Uddin, Md. Zia; Rana, Md. Sohel; Hossain, Shahadat; Effendi, Yunus; Idroes, Rinaldi; Tallei, Trina; Alqahtani, Ali M.; Emran, Talha Bin
title: Appraisals of the Bangladeshi Medicinal Plant Calotropis gigantea Used by Folk Medicine Practitioners in the Management of COVID-19: A Biochemical and Computational Approach
date: 2021-05-26
journal: Front Mol Biosci
DOI: 10.3389/fmolb.2021.625391
sha: a154cc216405525f4f7f16d6486ea7655a8353e0
doc_id: 788832
cord_uid: 17rj71x4

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was first recognized in Wuhan in late 2019 and, since then, had spread globally, eventually culminating in the ongoing pandemic. As there is a lack of targeted therapeutics, there is certain opportunity for the scientific community to develop new drugs or vaccines against COVID-19 and so many synthetic bioactive compounds are undergoing clinical trials. In most of the countries, due to the broad therapeutic spectrum and minimal side effects, medicinal plants have been used widely throughout history as traditional healing remedy. Because of the unavailability of synthetic bioactive antiviral drugs, hence all possible efforts have been focused on the search for new drugs and alternative medicines from different herbal formulations. In recent times, it has been assured that the M(pro), also called 3CL(pro), is the SARS-CoV-2 main protease enzyme responsible for viral reproduction and thereby impeding the host’s immune response. As such, M(pro) represents a highly specified target for drugs capable of inhibitory action against coronavirus disease 2019 (COVID-19). As there continue to be no clear options for the treatment of COVID-19, the identification of potential candidates has become a necessity. The present investigation focuses on the in silico pharmacological activity of Calotropis gigantea, a large shrub, as a potential option for COVID-19 M(pro) inhibition and includes an ADME/T profile analysis of that ligand. For this study, with the help of gas chromatography–mass spectrometry analysis of C. gigantea methanolic leaf extract, a total of 30 bioactive compounds were selected. Our analyses unveiled the top four options that might turn out to be prospective anti–SARS-CoV-2 lead molecules; these warrant further exploration as well as possible application in processes of drug development to combat COVID-19.

Coronavirus disease 2019 is an infectious disease triggered by infection with the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus. In late 2019, cases of a pneumonia of a then-unknown origin began to appear in the capital of China's Hubei province, Wuhan, later progressively spreading throughout the globe, resulting in the current COVID-19 pandemic (Hui et al., 2020) . COVID-19 is not the first respiratory disease that has caused an epidemic in the last two decades; previously, outbreaks of Middle East respiratory syndrome (MERS) and severe acute respiratory syndrome (SARS), both caused by other coronaviruses, have occurred (de Wit et al., 2016) . Coronaviruses are enveloped positivesense single-stranded RNA genome (26-32 kb) viruses (Su et al., 2016) . To date, four genera of coronaviruses have been identified: coronavirus genera α (alpha), coronavirus genera β (beta), coronavirus genera γ (gamma), and coronavirus genera δ (delta). Human coronaviruses (HCoVs) detected thus far include the HCoV-229E and NL63 in the α-coronavirus genus and MERS-CoV, SARS-CoV, HCoV-OC43, and HCoV-HKU1 in the β-coronavirus genus (Perlman and Netland, 2009) . Following virus genome sequencing in December 2019 of five hospitalized patients in Wuhan with the mysterious pneumonia, it was divulged the existence of a formerly unknown β-CoV strain existed in all of them. This isolated neoteric β-CoV manifests an 88% degree of specification to the sequence of two bat-derived SARS-like coronaviruses, bat-SL-CoVZC45 and bat-SL-CoVZXC21, and an almost 50% degree of specification to the sequence of MERS-CoV (Lu et al., 2020) . Liu et al., recently determined that there is a presence of the main protease enzyme in SARS-CoV-2 responsible for inhibiting SARS-CoV-2 replication. Therefore, targeting this specific enzyme may lead to potential therapeutic advantages in drug discovery and efficacy for managing this dangerous virus. Elsewhere, experiments with four drugs approved by the United States Food and Drug Administration, including praziquantel, perampanel, nelfinavir, and pitavastatin as potent inhibitors for SARS-CoV-2 M pro enzyme, are ongoing (Yusuf et al., 1994) . Recently, researchers are also reviewing medicinal plants for possible efficacy in the context of SARS-CoV-2 infection (Khaerunnisa et al., 2020) .

Having relatively less non-toxic property, synthetic bioactive natural compounds have diverse bioavailability. Various phytoconstituents like phenols, steroids, flavonoid molecules are present in the natural bioactive compounds and also they have been reported and screened for their possible therapeutic effects against various viral diseases. Certain natural compounds such as digitoxigenin, crocin, withanone, β-eudesmol, caffeic acid phenethyl ester and other various bioactive molecules have also been reported to interact with SARS-CoV-2 main protease (Tripathi et al., 2020) . Traditional medicine employs a holistic approach to the prevention of disease and act as a source of many components with a high therapeutic value. Globally, many studies have been conducted on these herbs and revealing unique active constituents that might activate the innate immune system through the stimulation of macrophages and lymphocytes and modulation of the cytokine profile also . The potent therapeutics that might be effective to combat SARS-CoV-2 infection include virus binding molecules, helicase inhibitors molecules or inhibitors targeting particular enzymes implicated in replication and transcription process of the virus, protease inhibitors of host cells, endocytosis inhibitors, short interfering RNA (siRNA), neutralizing antibodies and natural drugs or medicines (Dhama et al., 2020) .

Various kinds of secondary metabolites such glycosides, alkaloids, tannins, lignins and terpenes are present abundantly in various folklore medicinal plants and, thus, many studies have already explored the use of these plants in the treatment of various infections and diseases with the help of these secondary metabolites. Different extractable biochemical and bioactive compounds are contained in traditional medicinal plants might cure or prevent several viral diseases and infections by targeting certain specific sites of viruses. Despite the healing histories or scientific evidences, the use of medicinal plants and natural products as prophylactics, therapeutics and their multiple beneficial applications have gained strong momentum over the last few decades Brendler et al., 2020) .

Henceforth, compound isolation and subsequent further investigation for several pharmacological activities has been established as a crucial process in the areas of drug design and discovery of novel drugs as well. The present study considered some isolated phytochemicals from a medicinal plant Calotropis gigantea, a member of the Apocynaceae family, for the treatment of COVID-19 using computational biology techniques. C. gigantea, frequently called giant milkweed or crown flower, is a large shrub or small tree native to Bangladesh, India, Malaysia, Southern China, and Indonesia. The plant is also locally known as Gulancha. These large shrubs have also been extensively cultivated in tropical areas around the world and typically grow to eight to 15 feet tall.

As calotropin offers potent bioactivity, C. gigantea has been used by the folk practitioners in India for many years for a variety of purposes. In Ayurveda, the root and leaf of Calotropis procera have been being used in managing asthma and breathing problems such as shortness of breath. The bark of C. procera is also being applied in other types of diseases such as spleen and liver conditions. This plant has been reported to be effective as well in treating skin, circulatory, digestive, and neurological issues and was also used previously to treat fevers, nausea, vomiting, and diarrhea. However, these reports are usually based on folklore uses and, therefore, more scientific investigations are necessary to affirm the rational and clinical usefulness of the leaves, latex, and bark of this family. Moreover, recent research has suggested the use of calotropin as a contraceptive and as a promising medication for cancer patients (Wang et al., 2009) . In one study, the extracts of C. gigantea appeared to be strongly cytotoxic against non-small-cell lung carcinoma (A549), hepatocellular carcinoma (Hep G2), and colon carcinoma (HCT 116) . In other research, leaf extracts of this plant have also shown promising activities as cancer medications (Jacinto et al., 2011) . One species of the Calotropis genus, C. procera, displayed potent anti-human immunodeficiency virus 1 activity when the aqueous crude extract of the leaves was analyzed (Rakshit et al., 2010 ). An ethanolic extract of C. gigantea exhibited activity against herpes simplex type 1 viruses and vesicular stomatitis viruses (Ali et al., 1996) . A novel lignan glycoside, (+)-pinoresinol 4-O-[6″-O-vanilloyl]-β-dglucopyranoside, isolated from the latex of C. gigantea exhibited anti-influenza activity in an in vitro study (Parhira et al., 2014) . Meanwhile, several other research works have suggested the protease inhibitory effect of C. gigantea.

In alignment with these previous findings, we speculate that the phytochemicals obtained from C. gigantea might be able to impede the main protease enzyme of SARS-CoV-2 (Kaladhar et al., 2013; Gupta and Chaphalkar, 2016) . Through the process of gas chromatography-mass spectrophotometry (GC-MS), several compounds have been isolated for molecular docking. Molecular docking studies are of terrific significance in establishing plans for the layout of new drug substances; specifically, their strategic aim is to explore the binding affinity and the mode of binding of a small molecule within the specific binding sites of the receptor target of interest (Chikhale et al., 2020) . This methodology is based on the prediction of native ligands that pose or bind with the proper receptor binding site, i.e., to find the experimental ligand geometry associated with physicochemical molecular interactions. The aim of this study was to screen out potent bioactive ligands from the GC-MS compound list of C. gigantea that could be probable inhibitors of M pro .

The methanolic leaf extract of C. gigantea was inspected in a mass spectrometer (TQ8040; Shimadzu Corporation, Kyoto, Japan) with assistance from of electron-impact ionization method and a gas chromatograph (GC-17A; Shimadzu Corporation) fused with a silica capillary column made by silica with the specifications Rxi-5 ms, 0.25 m film, 30 m long, and internal diameter 0.32 mm coated with DB-1 (J&W). The oven temperature was set at 70°C (0 min); 10°C, 150°C (5 min); 12 and 200°C (15 min); 12 and 220°C (5 min), with a hold time of 10 min for each. A temperature of 260°C was set as the inlet temperature. The flow rate was 0.6 ml/min using helium gas at a constant pressure of 90 kPa. There was an interface temperature from GC to MS, which was adjusted at a constant 280°C. The scanning range of the MS was 40-350 amu and the ionization mode (EI) was established with the mass range of 50-550 m/z. The injection volume of the sample was 1 μl. During 50 min, the total GC-MS procedure was conducted. A comparison with the National Institute of Standards and Technology (NIST) GC-MS library version 08-S was performed for the identification of compounds in the peak areas .

The GC helped the blend of components to be separated into individual ones and, based on their masses, MS identified the molecules. As such, GC-MS is an analytical method combining two instruments to yield a powerful separation and identification process. From the GC-MS process, 74 compounds were isolated and 30 bioactive compounds were finally enlisted and further assessed regarding their role as ligands. Their lists as well as details of the structure and PubChem ID number are given in Table 1 .

Protein crystal structures are crucial for performing molecular docking simulation. In this study, we used the crystal structure of SARS-CoV-2 in concert with X77 (N-(4-tert-butylphenyl)-N-[(1R)-2-(cyclohexylamino)-2-oxo-1-(pyridin-3-yl)ethyl]-1Himidazole-4-carboxamide) (PDB ID: 6W63) with reference to a recent study by Dubey et al., who adopted the structure for performing docking studies by using some flavonoid compounds (Dubey and Dubey, 2020) . The three-dimensional crystal structure of SARS-CoV-2 M pro was downloaded from the Protein Data Bank in protein data bank (PDB) format having the PDB ID: 6W63 and was refined further using the protein preparation wizard of Maestro version 10.1 (Schrödinger, New York. NY, United States). Furthermore, the bond orders and charges together with the addition of hydrogen to the heavy metals; conversion of selenomethionines to the methionines; and, finally, the deletion of water were performed by using the force field OPLS_2005. Then, by setting the maximum heavy atom root mean-square-deviation to 0.30 Å, the minimization was carried out (Chowdhury et al., 2020) . In addition, the structures of the chemical compounds were obtained from the PubChem database (https://pubchem.ncbi.nlm.nih.gov) and then ligands were prepared by neutralizing the pH to 7.0 ± 2.0 and minimized by the application of force field OPLS_2005 embedded in Maestro version 10.1 (Dash et al., 2015) .

For prepared proteins, the receptor grids were calculated in such a way that various ligands predicted during docking were bound within the active site. In Glide, grids were brought out by keeping the default parameters of the Van der Waals scaling factor at 1.00 and the charge cutoff at 0.25, which were exposed to OPLS_2005 force field. On the centroid, a cubic box of specific dimensions was centered at the active site residues and the bounding box was set to 17 × 17 × 17 Å for docking experiments.

With the help of Maestro v 10.1 (Schrödinger), ligand docking was carried out, which were furthermore enforced to non-cis/ trans amide bonds and the Van der Waals scaling factor; thereafter, the partial charge and cutoff values were selected (Friesner et al., 2004; Friesner et al., 2006) . The ultimate scoring was executed on energy-minimized poses and represented as the glide score. For each ligand, the bestdocked pose with the lowest glide score value was listed.

The amino acids from the receptors that bind with the ligand are examined using PDBsum (http://www.ebi.ac.uk/pdbsum), a server that provides various pictorial analyses. This server integrates LIGPLOT to show the interactions of the analyzed complexes. A LIGPLOT diagram depicts the hydrogen bonds and non-bonded interactions between each ligand molecule and the protein residues in which it interacts with each ligand molecule.

To assess the bioactivity of the lead compounds, the calculation of their activity scores as kinase inhibitors and enzyme inhibitors was completed. With the help of the software Molinspiration (www.molinspiration.com; Molinspiration, Slovensky Grob, Slovak Republic), all parameters were scrutinized properly (Rakib et al., 2020a; Rakib et al., 2020b; Tareq et al., 2020) and calculated drug-likeness scores of each compound were compared with the specific activity of each compound as well as those of standard drugs (i.e., nelfinavir and lopinavir).

The QikProp module of the Maestro version 10.1 (Schrödinger) is an easily assessable absorption, distribution, metabolism, and excretion (ADME) prediction program that has been designed specifically to obtain certain descriptors linked to ADME. It is helpful to anticipate both the physicochemical and pharmacokinetic descriptors' relevant properties. QikProp follows Lipinski's rule of five for the evaluation of ADME properties, ascertaining drugsimilar activities of ligand molecules. ADME/T properties of the compounds (i.e., vitamin E acetate; methyl gamma-linolenate; clionasterol; juniper camphor; ethyl 4-fluoro-1-methyl-1Himidazole-5-carboxylate; bicyclo[3.3.1]nonane-2,4-dione; 9,9-dimethoxy-; alpha.-amyrin; moretenone; 6-methoxy-2,5,8trimethyl 2-(4,8,12trimethyltridecyl) 3,4-dihydrochromene; betulinaldehyde; bicyclo[4.3.0]nonane; 1 isopropenyl-4,5-dimethyl-5phenylsulfonylmethyl; tert-butyl(5-isopropyl-2-methylphenoxy) dimethylsilane; behenyl behenate; ergost-5-en-3-ol, (3.beta.,24R)-; erucic acid; glyceryl palmitate; methyl heneicosanoate; methyl palmitate; ethyl linoleate; 1,3-dihydroxypropan-2-yl (9E,12E, 15E)-octadeca 9,12,15-trienoate; palmitic acid; stearic acid; glyceryl monostearate; olean-12-en-3-ol; acetate, (3.beta.)-; methyl pentadecanoate; phytol, tridecanedial; trilinolein; methyl n-undecanoate; and vitamin E) were analyzed using the QikProp 3.2 module (Rahman et al., 2020) .

The inhibition constant in this study was determined by AutoDock 4.0 tools. AutoDock software calculates and predicts the interactions between a ligand molecule and a protein molecule, based on predefined parameters. The interactions between the molecules can be calculated at a specified user-defined region of the protein. This region must be defined by the user through the Grip map option or GridBox. Thus, the use of the GridBox at the binding site or active site, or at other essential protein regions, is essential to perform AutoDock analysis. Before executing the AutoDock, the ".pdb" files of the protein and ligand have to be moved into one folder. Analysis in AutoDock can be divided into the following categories: a) Initialising molecules; and b) Running AutoGrid.

Initialising the molecule involved the addition of hydrogen atoms. The ligand molecule required the addition of the Gasteiger charge, the identification of aromatic carbons, the detection of rotatable bonds, and the setting of the torsional degrees of freedom (TORSDOF) value. The protein must be initialised manually, whereas the ligand is automatically initialised when opened in the tool. The receptor and ligand files were saved in.pdbqt format. Next, the receptor was opened again using the "Grid" menu and the "Macromolecules" submenu, and "open" was selected. The ligand opens using the "Set Map Types" sub-menu and "Open Ligand." Then, the GridBox is set in AutoDock to cover the identified binding sites. AutoDock only analyses the ligand molecule's interactions and the amino acids that are present within the GridBox. The GridBox's size can be increased or decreased using the number of points in the X/Y/ Z dimension. The position of the GridBox can be adjusted to cover the binding site or binding residues, using the "Center GridBox" field that moves the GridBox in the X, Y, and Z-axis.

After setting up the GridBox, the AutoGrid file was saved. When in the "Grid" menu, select the "Output" and "Save GPF" submenu. The file saves in.gpf format. Meanwhile, the AutoDock file saving procedure is to choose the "Docking" menu and select "Lamarckian GA" in the "Output" menu. The file saves in.dpf format. The AutoGrid and AutoDock programs was run through a command prompt directed to the folder to be docked using the command: autogrid4 -p control.gpf -l control.glg. After the AutoGrid calculation was completed successfully, AutoDock was run using the command: autodock4 -p control.dpf -l control.dlg. The results of the docking calculations are obtained in a notepad file format, reporting the values for binding energy and the estimated inhibition constant.

From the methanolic leaf extract of C. gigantea, 52 compounds were eluted in between the retention time of 5 and 32 min ( Figure 1 ). Thirty bioactive compounds were selected for this study ( Table 1 ). The compound identification was carried out in comparison with the NIST GC-MS library version 08-S. The major phytochemicals included phenols, sterols, terpenoids, esters, and other organic compounds. The major compounds identified are trilinolein (18.79%), glyceryl monostearate (12.23%), erucic acid (11.64%), betulinaldehyde (10.38%), ergost-5-en-3-ol, (3.beta.,24R)-(7.03%), tridecanedial (6.90%), palmitic acid (6.78%), and methyl n-undecanoate (5.03%).

A virtual screening process, including molecular docking, is a feasible approach for identifying antiviral therapeutics among diverse peptide sequences available in the online-based database. Notably, the interaction type, and the binding energy in conjunction with the bond distance of interaction between a protein substrate and a ligand, can be evaluated through the utility of molecular docking. Therefore, screening a ligand among numerous plant sources based on binding energy assessment is possible within minimum time by running molecular docking to evaluate the best ligand candidates according to the docking interaction.

In this study, we determined the interaction between several isolated compounds of C. gigantea and the M pro enzyme. The docking results of the compounds that docked with the M pro enzyme of SARS-CoV-2 are presented in Tables 2, 3; Figures 2 through 5. Meanwhile, the two-and three-dimensional interactions of nelfinavir and lopinavir with the active site of the SARS-CoV-2 M pro are shown in Figures 6 and 7 .

Among the molecular docking results, nelfinavir and lopinavir were confirmed to have docked with the main protease and had docking scores of −7.596 and −8.251 kcal/mol, respectively. According to the docking scores, juniper camphor was better in comparison with the standard drugs. Bicyclo[4.3.0]nonane, 1 isopropenyl-4,5-dimethyl-5-phenylsulfonylmethyl, olean-12-en-3-ol, acetate, (3.beta.)-, and ethyl 4-fluoro-1-methyl-1H-imidazole-5-carboxylate also exhibited better scores in a moderate fashion. Moreover, methyl pentadecanoate shows lowest docking score (0.806 kcal/mol). If arranging these compounds according to their scores, then the arrangement would be, from highest to lowest, juniper camphor (−6.06 kcal/mol); bicyclo [4.3.0]nonane, 1 isopropenyl-4,5-dimethyl-5-phenylsulfonylmethyl (−5.808 kcal/mol); olean-12-en-3-ol, acetate, (3.beta.)-(−5.588 kcal/mol); and ethyl 4-fluoro-1-methyl-1H-imidazole-5carboxylate (−5.558 kcal/mol).

Again As seen in Table 3 , we also found that ethyl 4-fluoro-1methyl-1H-imidazole-5-carboxylate interacted with the protein through a single H-bond and six hydrophobic bonds as follows: HIS 41 (carbon), ASP 187 (fluorine), MET 165 (pi-sulfur), MET 49 (alkyl), HIS 41 (pi-pi-T-shaped), and ARG (amide-pistacked).

Specific parameters, protease inhibitors, and enzyme inhibitors of all selected compounds are shown in Table 4 .

Pharmacokinetic and physicochemical properties of the 30 selected bioactive compounds were screened for drugcandidacy using the QikProp module of the Schrödinger suite-Maestro version 10.1. The ligand molecule's similar drug attitude was determined using ADME properties by the QikProp module of Schrodinger Maestro v10.1. The ADME/T properties of the selected compounds were determined with the QikProp module of Schrodinger, shown in Table 5 .

The docking conformation of the ligands and receptors were determined by selecting the pose with the highest affinity (most negative Gibbs free energy). The results from each affinity value of the conformation were obtained in the notepad file. The visualization of three-dimensional (3D) and two-dimensional (2D) structures from the docking calculations were carried out. The visualization and potential interaction between each ligand molecule with protein residues results are presented in Supplementary Table S1; Supplementary Figure S1 . The interaction of the receptor's amino acid with the ligand in 3D and 2D were saved as image files. 

The current study determined the inhibition constant of major bioactive compounds from methanol extract of C. gigantea. toward the modeled SARS-CoV-2 M pro . The theoretical halfmaximal inhibitory concentration (IC 50 ) was calculated using AutoDock 4.0 tools, and the results of the analysis are shown in Table 6 .

Currently, SARS-CoV-2 infection has emerged as a global pandemic and requires an international public health emergence (Matveeva et al., 2020; Remali and Aizat, 2021) . However, different research organizations and pharmaceutical companies are working relentlessly for the management of the infection. By far, several vaccines have been made available as immunomodulation is characterized as one of the pivotal characteristics in the maintenance of human health and productivity Adhikari et al., 2020; Tahir Ul Qamar et al., 2020; Tallei et al., 2020 ) .

Researchers are learning new things about this virus every day about the mutation of the virus (Sinha et al., 2021) . Surprisingly, according to Centers for Disease Control and Prevention (CDC), multiple SARS-CoV-2 variants are available globally. In United Kingdom, the variant of SARS-CoV-2 showed a large number of mutations In South Africa, another variant emerged independently and share several characteristics with other strain. As the mutation of the virus is rising rapidly, it is crucial to identify potential therapeutic agents to overcome the pathetic situation. So far, we know that initially, some patients having attacked by COVID-19 may not figure out any symptoms, and hence they may carry the virus for two days or up to two weeks without noticing any symptoms (Tahmasbi et al., 2021) . Some common symptoms have explicitly been notified to COVID-19. These include shortness of breath, a low-grade fever that gradually increases in temperature, having a cough that gets more severe over time, and fatigue. Less common symptoms include chills, repeated shaking with colds, sore throat, headache, muscle aches and pains, loss of taste, and smell loss. Some people may need to seek emergency medical services to have any of the symptoms more severe, like trouble breathing with persistent pain or pressure in the chest and excessive drowsiness (--).

Although several research works are being carried out in renowned laboratories in various countries, the vaccine development for a disease generally requires a long time.

Despite plenty of experimental trials related to COVID-19 treatment and medications, scientists are still seeking specific therapeutic drugs. Several previous studies have already performed for the identification of potential therapeutic agents against SARS-CoV-2 from medicinal plants (Rabaan et al., 2020; Shoaib et al., 2020; Emirik, 2020) . Medicinal plants are abundant with plenty of phytocompounds and have been used to treat several kinds of diseases (Shree et al., 2020) . Different secondary metabolites, including lignin's, glycosides, alkaloids, terpenoids, and amino acids that are fundamental for the functioning and growth of plants as well as having numerous pharmacological attributes to fight against several abnormal conditions (Rahman et al., 2013; Emran et al., 2018; Diniz et al., 2021 ) . We had selected the plant C. gigantea based on their traditional use and ethnomedicinal value and their contribution as an antiviral drug for certain viruses that were analyzed before. For this study, with the help of the GC-MS technique, seventy-four compounds had been isolated, among which thirty bioactive compounds were enlisted and subjected to play their role as ligands. In recent times, researchers are working to find out some potential lead compounds from medicinal plants that are active against several enzymes and other proteins responsible for viral replication and growth (Gupta et al., 2020b; Krupanidhi et al., 2020; Pendyala and Patras, 2020 ) . Several studies reported the inhibitory potential of several herb derived phytochemicals against SARS-CoV-2 Gupta et al., 2020b; Jalali et al., 2020) . The advancement of numerous computational technology has affected vastly in biology and computational biology, incorporating both computational algorithms and biological data, which has emerged as a breakthrough in modern drug discovery. High-throughput screening (HTS) often requires the screening of the plenty of chemical compounds specific for a targeted receptor, utilizes molecular docking as essential tools (Brooijmans and Kuntz, 2013) . We predicted the interaction between the SARS-CoV-2 M pro and the targeted phytochemical through molecular docking simulation in the current study. The molecular docking analysis findings revealed that the selected phytocompounds were able to interact with the binding pockets of the SARS-CoV-2 M pro . The compound Juniper camphor exhibited the most significant interaction compared with other compounds. Like two standard drugs, nelfinavir, and lopinavir, it yielded hydrogen bonding with Met165 residue. Besides, our selected compounds were able to interact with the SARS-CoV-2 M pro , and previous other studies also suggested the residues were responsible for interacting with SARS-CoV-2 M pro (Joshi et al., 2020; Khan et al., 2020; Wu et al., 2020) . Molinspiration is used to predict the biological activity parameters like protease inhibitory activity and enzyme inhibitory action. The increase scores specify better inhibition properties (Yamamoto et al., 2004; Chang et al., 2019; Das et al., 2020; Gentile et al., 2020; Li et al., 2020) ; most of the targeted compounds were shown greater inhibition properties. QikProp module of the Schrödinger suite predicts ADMET properties related to molecular weight, partition coefficient (QPlogPo/w), stars, % oral absorption values, and Lipinski's properties (Rule of five). Different ADME parameters were predicted, as the molecular weight ranged between 200.32 and 649.178, where only one phytochemical was outside the specified range. The estimated number of hydrogen bond donors (HB donors) by the solute to water molecules in an aqueous solution lay in the range between 1 and 2. In contrast, in the case of hydrogen bond acceptor (HB acceptors), the number was found between 0.75 and 5.5. The value of total solvent accessible surface area (SASA) ranged from 390.235 to 1729.5, where only Behenyl behenate was recorded outside the acceptable range. The value of the octanol/water partition coefficient (QPlog Po/w) lay in the range between 0.545 and 15.728. All the bioactive compounds have shown good Brain/blood partition coefficient (QPlogBB), the drug's ability to pass through the blood-brain barrier (BBB). The QPlogBB value ranged from −1.633 to 0.986 suggested all the compounds are within the specified limits; aqueous solubility, log S (QPlogS) ranged between −27.42 and −0.732. Interestingly, all the compounds except Behenyl behenate fulfilled the Lipinski's ROF as it was found that Lipinski violations were ≤1. The percent of human oral absorption was high for all the selected compounds. All the phytochemicals exceeded 80% and showed more oral absorption than both the standard drugs. Notably, some phytochemicals exhibited 100% oral bioavailability, which enhanced the drug candidacy of those compounds.

Enzyme activation or inhibition is the result of the binding interaction between a small molecule ligand and an enzyme protein. The protein is supposed to be a receptor, and thus binding with the ligand may result in agonism or antagonism. As a form of bioinformatics analysis, Molecular docking simulation reveals the binding affinities of ligand molecules with a specific receptor, where the lower binding energy predicts the greater binding affinity (Brooijmans and Kuntz, 2003; Lalitha and Sivakamasundari, 2010; Ahmed et al., 2019; Rakib et al., 2020c; Rakib et al., 2020d; Jahan et al., 2020; Dutta et al., 2019) . Despite showing lower binding affinities than nelfinavir and lopinavir, our selected compounds interacted with the active pockets of the SARS-CoV-2 M pro enzyme like the standard compounds. Among the selected 30 ligands, 16 ligands docked with active pockets of SARS-CoV-2 M pro enzymes. Four of them, namely Juniper camphor, Ethyl 4fluoro- Bicyclo[4.3.0] nonane, 1 isopropenyl-4,5-dimethyl-5-phenylsulfonylmethyl and Olean-12-en-3-ol, acetate, (3.beta.) possessed docking scores more likely to nelfinavir. Moreover, those compounds also followed the Lipinski's rule of five for drug-likeness properties. Hence, these four may be potential inhibitors of COVID-19, which warrants further investigation. While so many synthetic bioactive compounds are undergoing clinical trials, including medicinal plants a Molecular weight (acceptable range: < 500). b Hydrogen bond donor (acceptable range: ≤5). c Hydrogen bond acceptor (acceptable range: ≤10). d Total Solvent Accessible Surface Area in using a probe with a 1.4 radius (acceptable range: 300-1,000). e Predicted octanol/water partition coefficient (acceptable range: − 2-6.5).

f Predicted blood-brain partition co-efficient (acceptable range: −3-1.2). g Predicted aqueous solubility, S in mol/dm −3 (acceptable range: −6.5-0.5). h Predicted human oral absorption on 0-100% scale (<25% is poor and >80% is high). derived phytocompounds, this purely bioinformatic work was limited by a lack of direct experimental evidence. We were unable to conduct follow-up experiments using other techniques, such as Western blotting or qRT-PCR, to confirm our bioactive compounds' antiviral properties generated only from computational prediction. This study could be expanded to examine the antiviral potentiality of plant-derived bioactive compounds, specifically from several traditional medicinal plants available in the Bangladesh region. Based on results from this study, the top compounds that might turn out to be prospective anti-SARS-CoV-2 lead molecules warrant further investigation and possible application in drug development processes against COVID-19.

Increasing use of HTS and molecular docking simulation is emerged as a new era in the field of drug discovery. Medicinal plant derived phytochemicals are used for the treatment from since the dawn of civilization. Recently, the enriched extraction procedures and sophisticated computer-aided program have already added a new dimension in the field of drug discovery. To sum up, this study provides significant data for the use of isolated phytocompounds from C. gigantea for the treatment of COVID-19 by targeting one of the crucial enzymes that is essential replication of SARS-CoV-2. Further experimental analysis might be needed to justify the activity of the phytochemicals at the molecular level. In addition, many more such bioactive ingredients from medicinal plants existing in the rich Bangladeshi biodiversity need to be further explored. The study initiated the window toward plant-based therapy against COVID-19, though extensive research and wet lab validation need to make it usable for the patient.

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material. 

Potential Roles of Medicinal Plants for the Treatment of Viral Diseases Focusing on COVID -19: A Review

In Vivo and In Vitro Pharmacological Activities of Tacca Integrifolia Rhizome and Investigation of Possible lead Compounds against Breast Cancer through In Silico Approaches

Antiviral and Cytotoxic Activities of Some Plants Used in Malaysian Indigenous Medicine

Botanical Drugs and Supplements Affecting the Immune Response in the Time of COVID-19: Implications for Research and Clinical Practice

Molecular Recognition and Docking Algorithms

Molecular Recognition and Docking Algorithms

Antiviral Drug Discovery: Norovirus Proteases and Development of Inhibitors

Sars-cov-2 Host Entry and Replication Inhibitors from Indian Ginseng: Anin-Silicoapproach

Drug Repurposing Approach against Novel Coronavirus Disease (COVID-19) through Virtual Screening Targeting SARS-CoV-2 Main Protease

An Investigation into the Identification of Potential Inhibitors of SARS-CoV-2 Main Protease Using Molecular Docking Study

Evolution of Selective COX-2 Inhibitor from Alangium Salvifolium: an In Silico Approach

SARS and MERS: Recent Insights into Emerging Coronaviruses

Medicinal and Therapeutic Potential of Herbs and Plant Metabolites/Extracts Countering Viral Pathogens -Current Knowledge and Future Prospects

Coronavirus Disease 2019-COVID-19

Bioactive Terpenes and Their Derivatives as Potential SARS-CoV-2 Proteases Inhibitors from Molecular Modeling Studies

Computation Screening of Narcissoside a Glycosyloxyflavone for Potential Novel Coronavirus

Potential Therapeutic Effect of Turmeric Contents against SARS-CoV-2 Compared with Experimental COVID-19 Therapies: In Silico Study

Antipyretic Effects of a Chloroform Extract from the Leaves of Urena Sinuata in Rodents

Pharmacological evidence for the use of Cissus assamica as a medicinal plant in the management of pain and pyrexia

Glide: a New Approach for Rapid, Accurate Docking and Scoring

Method and Assessment of Docking Accuracy

Extra Precision Glide: Docking and Scoring Incorporating a Model of Hydrophobic Enclosure for Protein−Ligand Complexes

Putative Inhibitors of SARS-CoV-2 Main Protease from a Library of marine Natural Products: A Virtual Screening and Molecular Modeling Study

Anti-diabetic Activity of Calotropis Gigantea in Human Whole Blood

Secondary Metabolites from Spice and Herbs as Potential Multitarget Inhibitors of SARS-CoV-2 Proteins

The Continuing 2019-nCoV Epidemic Threat of Novel Coronaviruses to Global Health -the Latest 2019 Novel Coronavirus Outbreak in Wuhan

Cytotoxic Cardenolide and Sterols from Calotropis Gigantea

Pharmacological Properties, and In Silico Studies of Chukrasia Velutina Leaves: A Novel Source for Bioactive Agents

A Pharmacology-based Comprehensive Review on Medicinal Plants and Phytoactive Constituents Possibly Effective in the Management of COVID -19

Discovery of Potential Multi-Target-Directed Ligands by Targeting Host-specific SARS-CoV-2 Structurally Conserved Main Protease

In Vitro protease Inhibition, Modulation of PLA2 Activity and Protein Interaction Studies of Calotropis Gigantea

Potential Inhibitor of COVID-19 Main Protease (M pro ) from Several Medicinal Plant Compounds by Molecular Docking Study

Identification of Chymotrypsin-like Protease Inhibitors of SARS-CoV-2 via Integrated Computational Approach

Screening of Phytochemical Compounds ofTinospora Cordifoliafor Their Inhibitory Activity on SARS-CoV-2: Anin Silicostudy

Calculation of Molecular Lipophilicity and Drug Likeness for Few Heterocycles

The Epidemic of 2019-Novel-Coronavirus (2019-nCoV) Pneumonia and Insights for Emerging Infectious Diseases in the Future

The Crytal Structure of 2019-NCoV Main Protease in Complex with an Inhibitor N3. RCSB Protein Data Bank

Genomic Characterisation and Epidemiology of 2019 Novel Coronavirus: Implications for Virus Origins and Receptor Binding

Search of Herbal Anti-SARS-cov2 Compounds

Zhi-Hong Jiang. In Vitro Antiinfluenza Virus Activities of a New Lignan Glycoside from the Latex of Calotropis Gigantea

In Silico Screening of Food Bioactive Compounds to Predict Potential Inhibitors of COVID-19 Main Protease (Mpro) and RNA-dependent RNA Polymerase (RdRp)

Coronaviruses post-SARS: Update on Replication and Pathogenesis

SARS-CoV-2/COVID-19 and Advances in Developing Potential Therapeutics and Vaccines to Counter This Emerging Pandemic

Biological Evaluation, DFT Calculations and Molecular Docking Studies on the Antidepressant and Cytotoxicity Activities of Cycas Pectinata Buch

Antioxidative, Antimicrobial and Cytotoxic Effects of the Phenolics of Leea Indica Leaf Extract

Biochemical and Computational Approach of Phytocompounds from Tinospora Crispa in the Management of COVID-19

Antipyretic and Hepatoprotective Potential of Tinospora Crispa and Investigation of Possible lead Compounds through

Epitope-Based Immunoinformatics Approach on Nucleocapsid Protein of Severe Acute Respiratory Syndrome-Coronavirus-2

Frontiers in Molecular Biosciences | www.frontiersin.org

Immunoinformatics-guided Design of an Epitope-Based Vaccine against Severe Acute Respiratory Syndrome Coronavirus 2 Spike Glycoprotein

Anti-HIV-1 and Antimicrobial Activity of the Leaf Extracts of Calotropis Procera

A Review on Plant Bioactive Compounds and Their Modes of Action against Coronavirus Infection

Properties of Ethno-Medicinal Plants and Their Bioactives-Possible Use for Prevention and Treatment of COVID-19: A Review

Targeting COVID-19 (SARS-CoV-2) Main Protease through Active Phytochemicals of Ayurvedic Medicinal Plants -Withania somnifera(Ashwagandha),Tinospora cordifolia(Giloy) andOcimum sanctum(Tulsi) -a Molecular Docking Study

Potential Leads from Liquorice against SARS-CoV-2 Main Protease Using Molecular Docking Simulation Studies

Epidemiology, Genetic Recombination, and Pathogenesis of Coronaviruses

Structural Basis of SARS-CoV-2 3CLpro and Anti-COVID-19 Drug Discovery from Medicinal Plants

Herbal Medication to Enhance or Modulate Viral Infections

Potential of Plant Bioactive Compounds as SARS-CoV-2 Main Protease (Mpro) and Spike (S) Glycoprotein Inhibitors: A Molecular Docking Study

Chemical Profiles, Pharmacological Properties, and In Silico Studies Provide New Insights on Cycas Pectinata

Herbal Immunomodulators -A Remedial Panacea for Designing and Developing Effective Drugs and Medicines: Current Scenario and Future Prospects

Identification of Bioactive Molecule fromWithania somnifera(Ashwagandha) as SARS-CoV-2 Main Protease Inhibitor

Cytotoxicity of Calotropin Is through Caspase Activation and Downregulation of Anti-apoptotic Proteins in K562 Cells

Analysis of Therapeutic Targets for SARS-CoV-2 and Discovery of Potential Drugs by Computational Methods

HIV Protease Inhibitor Nelfinavir Inhibits Replication of SARS-Associated Coronavirus

Medicinal Plants of Bangladesh

The authors would like to thank Cambridge Proofreading ® & Editing LLC. (https://proofreading.org/) for editing a draft of this manuscript.

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2021.625391/ full#supplementary-material Conflict of Interest: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.Copyright © 2021 Dutta, Nezam, Chowdhury, Rakib, Paul, Sami, Uddin, Rana, Hossain, Effendi, Idroes, Tallei, Alqahtani and Emran. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.