key: cord-0790273-pjtnta8n authors: Roy, Arpita; Menon, Tarunya title: Evaluation of bioactive compounds from Boswellia serrata against SARS-CoV-2 date: 2021-11-17 journal: Vegetos DOI: 10.1007/s42535-021-00318-7 sha: d580f2d83de6c14e2398528263fb344ec95ed61d doc_id: 790273 cord_uid: pjtnta8n With the COVID-19 pandemic still wreaking havoc worldwide, new variants being discovered every month in some parts of the globe due to the mutating nature of the virus. There is no specific solution for this highly transmissible disease. In search of a lead molecule for the discovery and development of drug, extensive research is being conducted throughout the world. Many synthetic drugs are already in clinical trials and some are utilized for the treatment of this viral infection. Apart from synthetic drugs, phytocompounds from plants act as a potential drug candidate which can inhibit the growth of virus and thus able to prevent the viral infection. In this study, 26 ligands (bioactive compounds) from Boswellia serrata (an important medicinal plant) were tested against SARS-CoV-2 by using computational method. Selected ligands were shortlisted using Lipinski’s rule and then subjected to molecular docking against one of the main proteins of SARS-CoV-2, i.e., Mpro. Out of these compounds, Euphane, Ursane, α-Amyrin, Phytosterols, and 2,3-Dihydroxyurs-12-en-28-oic acid were potential to inhibit the Mpro activity with binding energies of − 10.47 kcal/mol, − 10.41 kcal/mol, − 9.99 kcal/mol, − 9.94 kcal/mol and − 9.72 kcal/mol respectively. A comparative study was performed using the best five ligands against four possible drug targets of SARS-CoV-2. It was found that Euphane showed highest negative binding energy against all the four crucial targets of SARS-CoV-2. Further, in-vitro experimentation is required to validate the use of Euphane as a potent drug against SARS-CoV-2. COVID-19 is one of the most transmissible and deadly diseases in recent years, resulting in the millions of deaths worldwide (Dong et al. 2020) . Still, there is no targeted drug or effective treatment developed for this virus because of its pathogenesis and mutation efficiency. Therefore, it is of utmost priority to identify and utilize effective antiviral agents that can be approved as drug candidates against SARS-CoV-2 in order to curb the spread of this pandemic (Jin et al. 2020; . SARS-CoV-2 is an enclosed ss-RNA virus with a nonsegmented genome of almost 30 kb in size (Andersen et al. 2020 ). The structure of SARS-CoV-2 is spherical and belongs to genus betacoronavirus . SARS-CoV and SARS-CoV-2 share the same spike protein's structure attached to the viral envelope. It recognizes and bind to Angiotensin-Converting Enzyme-2 (ACE2) and host cell's protease enzyme . SARS-CoV-2 has fourteen Open Reading Frames (ORFs) which codes for four structural proteins: spike (S), nucleocapsid (N), envelope (E), and membrane (M) proteins . It also has 5′ and 3′ untranslated regions (UTR), a poly (A) tail and many unknown non-structural ORFs. The S protein, around 150 kDa in size, is responsible for binding of virus to the receptors present on the host cells, which then allows the virus to enter into the host (Kirchdoerfer et al. 2016) . N protein is the sole protein that fuses with the RNA genome. It aids in the complete virion transformation by assembling the virus via the process of budding. The confirmation of virus' envelope is due to the most profuse structural protein (M protein) which is about 25-30 kDa in size and consists of three transmembrane domains . is also involving in the assembling and budding of virus (Venkatagopalan et al. 2015) . In addition, ORF codes for two polyproteins that build the complex involved in viral replication. These then proliferate into new viral particles after being processed further proteolytically by coronavirus main protease (Mpro) enzyme and papain-like protease (PLpro) enzyme (Hilgenfeld 2014). Mpro is one of the main enzymes of coronaviruses essential to the viral replication and transcription process. SARS-CoV-2 Mpro is similar to SARS-CoV Mpro with about 96% structural similarity (Ullrich and Nitsche 2020) . Thus, drug development studies aim to find therapeutic candidates that target this enzyme in order to inhibit viral growth. (Jin et al. 2020) . In-silico drug design technologies were used to develop an inhibitor against Mpro of SARS-CoV-2. The active site of SARS-CoV-2 Mpro is determined by the catalytic dyad which is responsible for the maturation of virus (Mirza and Froeyen 2020) . Plants are the most abundant source of bioactive compounds which are used for the potential drug development (Roy and Bharadvaja 2017; Roy 2018; Roy et al. 2018; Roy and Bhatia 2021) . Ayurveda has proved to be a fast-growing branch of medicine for treatment of various diseases due to its advantages, such as inexpensiveness, efficiency with no adverse effects, and ability to cater exponentially growing population (Mirza and Froeyen 2020) . Boswellia serrata is an Ayurvedic herb belonging to Burseraceae family and commonly known as 'Indian Frankincense' or 'Indian Olibanum.' This plant is native to the dry states of Peninsular India, such as Gujarat, Andhra Pradesh, Jharkhand, Madhya Pradesh and, Chhattisgarh. This traditional plant found in the Eastern part of the globe is proven to have anti-inflammatory effects, helps in relieving congestion due to colds and flu, act as a disinfectant, and even reduce stress and anxiety (Al-Yasiry and Kiczorowska 2016) . Biological activities of this plant are due to the presence of various bioactive compounds which includes monoterpenes, diterpenes, sterols, triterpenes, tetracyclic triterpenes like euphane, and pentacyclic triterpenes (Sultana et al. 2013) . In-silico analysis using molecular docking contributes immensely contributes to the drug development process by first screening the vast database of compounds and then reducing the burden of in vitro experiments Roy and Bhatia 2021; Bhatia et al. 2021 ). These tools apply inbuilt scientific calculation algorithms to the interactions of docked complex between receptor and ligand/drug molecule. The present study aims to find the potential of bioactive compounds of Boswellia serrata against Mpro enzyme of SARS-CoV-2 using computational approach. To determine the interactions of selected molecules against Mpro, computational techniques such as ADME analysis and molecular docking analysis were performed. The 3D structure of Mpro enzyme of SARS-CoV-2 from RCSB Protein Data Bank (PDB id: 6LU7) was extracted. Structure consists of a homodimer composed of two unique protein chains-A consisting of 306 amino acids complexed with an inhibitor N3 molecule. Other targets including prefusion SARS-CoV-2 spike protein (PDB id: 6VSB), SARS-CoV-2 spike glycoprotein (PDB id: 6VXX) and SARS-CoV-2 chimeric receptor-binding domain (PDB id: 6VW1) were also taken into consideration for a comparative study ( Table 1) . A total of 26 bioactive compounds from Boswellia serrata were selected. 3D structures for docking purposes were obtained in '.pdb' format. All the information on phytocompounds was collected from the IMPPAT: Indian Medicinal Plants, Phytochemistry and Therapeutics Database (Mohanraj et al. 2018). For evaluating compounds as potential ligands, their Absorption, Distribution, Metabolism, and Excretion properties (ADME) were determined. Lipinski's rule was applied using a web-based tool named SwissADME (https:// www. swiss adme. ch/). Any ligands showing violations of Lipinski's rule were eliminated from further studies. To predict the active binding site for most stable ligand-protein interaction, Computed Atlas for Surface Topography of proteins (CASTp) was used (Tian et al. 2018) . 10 residues were found in the active binding pocket of Mpro which were identified as follows: THR24, THR26, PHE140, ASN142, GLY143, CYS145, HIS163, HIS164, GLU166, HIS172. Molecular docking is a vital computational tool that is utilized for prescreening of compounds in the drug discovery process. It helps in the selection and filtration of potential inhibitors for targeted drug discovery by assessing the binding affinity of protein-ligand complex. Before the process of docking, preparation of protein and ligands was done. Autodock 4.2 was used for docking analysis. Preparing the optimized protein structure involved removing of inhibitor N3 and H 2 O molecules, adding polar hydrogens, and calculating Kollman charges and Gasteiger charges. Energy minimization of protein was done by using Swiss PDB Viewer. A grid box of 54 × 53 × 38 was calibrated as per the active site prediction of protein with 0.375 Å spacing with the center X: − 11.164, Y: 16.901, Z: 65.46. The Docking Log Files (DLG) were produced for checking the binding energies further search parameter was assigned as "Genetic algorithm" and Lamarckian GA to run the output. The docked complex with maximum negative binding energy was then chosen to visualized the 2D structure to evaluate the protein-ligand interactions. Then for comparative study, three other receptors (6VSB, 6VXX, and 6VW1) of SARS-CoV-2 were docked against selected ligands. Grid box dimensions of 60 × 60 × 60 and coordinates X: 204.457, Y: 199.799, Z: 246.898, 80 × 80 × 80, coordinates X: 223.221, Y: 190.083, Z: 263.907 and 75 × 75 × 75, coordinates X: 85, Y: − 10, Z: 180 respectively were selected with each having grid point spacing of 0.375 Å. All ligands were evaluated based on their physicochemical properties, i.e., ADME analysis before docking analysis. Compounds which did not violate two or more parameters of Lipinski's rule were then selected for further analysis. Compound must possess a molecular weight < 500 Da, Log P < 5, H-bond donors < 5, and hydrogen bond acceptors < 10 in order to follow Lipinski's rule (Table 2) . A molecular docking process was then performed with the selected ligands after ADME analysis. Ligands were ranked according to their binding energies, i.e., one with the least binding energy was ranked at the top, which reflected the maximum stability of the docked complex (Table 3) . It was then compared with the control N3 (native inhibitor) which revealed that ten compounds had more binding affinity to the protein because of its higher binding energy as compared to N3 (− 8.15 kcal/mol). Euphane showed binding energy of − 10.47 kcal/mol, and four different types of bonding interactions were observed, including van der Waals, pi-sigma, alkyl, and pi-alkyl. HIS163, MET165, CYS155, and LEU27 were bonded through pi-alkyl and alkyl bonds, and the remaining residues were weakly bound by van der Waals interactions (Fig. 1) . Ursane showed binding energy of − 10.41 kcal/mol. Only three types of hydrophobic interactions were identified. HIS163, CYS145, HIS41, and MET165 are involved in Pi-alkyl and alkyl bonding, while the remaining twelve residues interact with the ligand via van der Waal interactions (Fig. 1) . α-Amyrin-Mpro docked complex showed binding energy of − 9.99 kcal/mol. HIS163 was involved in the H-bonding and CYS145, and MET49 were interacting via the Alkyl bond. Thirteen other residues were interested in maintaining the stability of complex via van der Waals forces (Fig. 1) . Phytosterol indicated the lowest binding energy of − 9.94 kcal/mol. Five different types of bonding between amino acid residues and ligand atoms were observed. THR26 was involved in H-bonding. HIS41 was bonded with the ligand via Pi-sigma and Pi-alkyl bonds. MET49, LEU27, CYS145, and MET165 were also engaged in Pi-alkyl and alkyl interactions. The other thirteen residues in the vicinity stabilize the complex via van der Waal forces (Fig. 1) . 2,3-Dihydroxyurs-12-en-28-oic acid showed binding energy of − 9.72 kcal/mol. Five kinds of bonding interactions were seen. PHE140 and HIS41 were associated with hydrogen-bond formations. HIS163, CYS145, and LUE27 formed alkyl and pi-alkyl bonds with the ligand and helped in stabilizing the complex. ASN142 and HIS41 form carbon-hydrogen bonds with the ligand. Twelve other residues could be visualized interacting via van der Waals forces around the ligand (Fig. 1) . The antiviral drug ritonavir, which has already surpassed in-vitro screening and undergoes clinical trials reported to have amino acid residue interactions with Mpro where (thiazoly-5-yl) methylcarbamate oxygen was involved in hydrogen bonding with GLY143 and CYS145 and interacted with THR25, THR26, and LEU27 residues. The other atoms were stabilized at the active site by Pi alkyl and Pi sigma bonds with HIS41 (Kanhed et al. 2021) . One study reported molecular docking analysis of phytocompounds present in Caesalpinia minax. Docking analysis revealed that 5,7-dimethoxyflavanone-4′-O-dglucopyranoside has binding energy of − 9.23 kcal/mol. It was stabilized by hydrophobic interactions with other residues such as PHE140, LEU141, ASN142. The inhibitor α-ketoamide is a standard inhibitor with an IC50 value of 0.67 ± 0.18 μM against Mpro enzyme of SARS-CoV-2 and binding energy of − 8.24 kcal/mol ). In another study, similar amino acids were involved in the interaction where residues ASN142, CYS145, and GLU166 were bonded via hydrogen bonding and hydrophobic interactions via THR26, HIS41, PHE140, LEU141, GLY143, SER144, HIS163, HIS164, MET165, PRO168, ASP187, GLN189, THR190 and GLN192 (Gurung et al. 2020) . The top five compounds which showed maximum binding affinities i.e. Euphane, Ursane, α-Amyrin, Phytosterol, and 2,3-Dihydroxyurs-12-en-28-oic were further selected for comparative analysis. After obtaining binding energy of top five compounds docked with Mpro, the compounds with maximum binding energies were then docked against three other protein receptors of SARS-CoV-2, namely 6VSB, 6VXX, and 6VW1, respectively, using Autodock (Table 4 ). The results were then verified using PatchDock server, which is an open-access tool with a geometry complementarity of protein-ligand molecule-based algorithm online. (http:// bioin fo3d. cs. tau. ac. il/ Patch Dock/). Its principle is built around the geometry of complementarity criteria for docking (Table 5) . Patchdock evaluates the features of docked complexes which suggests the optimum molecular interaction between protein and ligand that reports the maximum interaction area covered by the ligand with protein having least steric hindrance. Atomic Contact Energy (ACE) is also evaluated which is defined as the amount of desolvation-free energy required to transfer the ligand from exterior (water) to the protein's interior. PatchDock analysis also validated better efficacy of Euphane amongst the other compounds of Boswellia serrata. Euphane has demonstrated the most potential ligand with highest average patch dock score compared to other phytocompounds. The Patchdock software was used in various docking analysis. In a recent study on in-silico docking analysis of Nigellidine against spike glycoprotein, Patchdock software was used, and results were analyzed and compared to other phytocompounds of Nigella sativa (Maiti et al. 2020 ). The ACE value reported for nigellidine with spike glycoprotein was − 340.50. Different values included − 181.08 and − 269.93 for thymoquinone and dithymoquinone, respectively. Nigellicine was also compared but showed lower affinity than nigellidine. Maiti et al. (2020) concluded that nigellidine was a more potent drug against SARS-CoV-2 than thymoquinone, dithymoquinone, and nigellicine. Another study reported inhibition of SARS-CoV-2 spike glycoprotein by tea flavonoids. The top 20 docked results obtained by Patchdock for each complex were selected based on their geometric complementarity analysis and those with the most negligible ACE value. The ligands chosen for the study showed high binding affinities against spike protein (Maiti and Banerjee 2021) . The highest ACE value was reported for theaflavin digallate (− 465.17) and theaflavin monogallate (− 434.42), followed by Epigallocatechin 3-gallate (− 407.58). On the other hand, a drug already in use to treat COVID-19, hydroxychloroquine, showed an ACE value of only − 293.32 (Maiti and Banerjee 2021) . In a study, Tucatinib, a drug developed for HER2-positive breast cancer, was repurposed against SARS-CoV-2. The docking results with the Mpro enzyme using PatchDock reported docking score of 5640 and ACE value of − 348.62 (Alsalme et al. 2020 ). Studies have been conducted in order to find potential drug candidate against SARS-CoV-2. Most of the studies are targeted on finding a cure by using chemically synthesized or previously used drugs (drug repurposing). Since time immemorial, bioactive compounds from plants have also shown great efficacy in controlling diseases, thus utilizing phytocompounds for drug development can help combat COVID-19. From the present study, it can be concluded that the top five ligands, except α-Amyrin, bind to the catalytic dyad amino acid residues of Mpro by different bonding interactions. A comparative study using AutoDock and PATCHDOCK revealed the affinity of top five ligands. It was found that Euphane possesses most significant inhibitory potential against all of four receptors. Further in-vitro experiments are required to validate potential of Euphane as a clinical drug for the treatment of COVID-19. Funding Dr. Arpita Roy is thankful to Sharda University for providing seed fund (Seed fund-4 2001 (SUSF2001/12)). Modelling the structural and reactivity landscapes of tucatinib with special reference to its wavefunction-dependent properties and screening for potential antiviral activity Frankincense-therapeutic properties The proximal origin of SARS-CoV-2 Evaluation of different phytochemicals against BRCA2 receptor An interactive web-based dashboard to track COVID-19 in real time In-silico analysis of selected alkaloids against main protease (Mpro) of SARS-CoV-2 Molecular docking analysis of selected phytochemicals against SARS-CoV-2 M pro receptor Unravelling lead antiviral phytochemicals for the inhibition of SARS-CoV-2 Mpro enzyme through in-silico approach From SARS to MERS: crystallographic studies on coronaviral proteases enable antiviral drug design Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors Identification of potential Mpro inhibitors for the treatment of COVID-19 by using systematic virtual screening approach Pre-fusion structure of a human coronavirus spike protein Molecular immune pathogenesis and diagnosis of COVID-19 Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding Epigallocatechin-gallate and theaflavingallate interaction in SARS CoV-2 spike-protein central-channel with reference to the hydroxychloroquine interaction: bioinformatics and molecular docking study Active-site molecular docking of Nigellidine with nucleocapsid-NSP2-MPro of COVID-19 and to human IL1R-IL6R and strong antioxidant role of Nigella-sativa in experimental rats Structural elucidation of SARS-CoV-2 vital proteins: computational methods reveal potential drug candidates against main protease, Nsp12 polymerase and Nsp13 helicase IMPPAT: A curated database of Indian medicinal plants Role of medicinal plants against Alzheimer's disease Establishment of the shoot and callus culture of an important medicinal plant Plumbago zeylanica In silico analysis of plumbagin against cyclindependent kinases receptor Qualitative and quantitative phytochemical analysis of Centella asiatica Structural basis of receptor recognition by SARS-CoV-2 Boswellia serrata Roxb. a traditional herb with versatile pharmacological activity: a review CASTp 3.0: computed atlas of surface topography of proteins The SARS-CoV-2 main protease as drug target Coronavirus envelope protein remains at the site of assembly Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved α-ketoamide inhibitors Network-based drug repurposing for novel coronavirus 2019-nCoV/SARS-CoV-2 The authors declare no conflict of interest.