key: cord-0866370-8dl3mbr9
authors: Zothantluanga, James H.; Gogoi, Neelutpal; Shakya, Anshul; Chetia, Dipak; Lalthanzara, H.
title: Computational guided identification of potential leads from Acacia pennata (L.) Willd. as inhibitors for cellular entry and viral replication of SARS-CoV-2
date: 2021-10-09
journal: Futur J Pharm Sci
DOI: 10.1186/s43094-021-00348-7
sha: 4bf61d3edaa1aca8e078989a17f87351d7bf7f90
doc_id: 866370
cord_uid: 8dl3mbr9

BACKGROUND: Coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) started in 2019 and is still an on-going pandemic. SARS-CoV-2 uses a human protease called furin to aid in cellular entry and its main protease (M(pro)) to achieve viral replication. By targeting these proteins, scientists are trying to identify phytoconstituents of medicinal plants as potential therapeutics for COVID-19. Therefore, our study was aimed to identify promising leads as potential inhibitors of SARS-CoV-2 M(pro) and furin using the phytocompounds reported to be isolated from Acacia pennata (L.) Willd. RESULTS: A total of 29 phytocompounds were reported to be isolated from A. pennata. Molecular docking simulation studies revealed 9 phytocompounds as having the top 5 binding affinities towards SARS-CoV-2 M(pro) and furin. Among these phytocompounds, quercetin-3-O-α-L-rhamnopyranoside (C_18), kaempferol 3-O-α-L-rhamnopyranosyl-(1 → 4)-β-D-glucopyranoside (C_4), and isovitexin (C_5) have the highest drug score. However, C_18 and C_4 were not selected for further studies due to bioavailability issues and low synthetic accessibility. Based on binding affinity, molecular properties, drug-likeness, toxicity parameters, ligand interactions, bioavailability, synthetic accessibility, structure–activity relationship, and comparative analysis of our experimental findings with other studies, C_5 was identified as the most promising phytocompound. C_5 interacted with the active site residues of SARS-CoV-2 M(pro) (GLU166, ARG188, GLN189) and furin (ASN295, ARG298, HIS364, THR365). Many phytocompounds that interacted with these amino acid residues were reported by other studies as potential inhibitors of SARS-CoV-2 M(pro) and furin. The oxygen atom at position 18, the –OH group at position 19, and the 6-C-glucoside were identified as the pharmacophores in isovitexin (also known as apigenin-6-C-glucoside). Other in-silico studies reported apigenin as a potential inhibitor of SARS-CoV-2 M(pro) and apigenin-o-7-glucuronide was reported to show stable conformation during MD simulations with SARS-CoV-2 M(pro). CONCLUSION: The present study found isovitexin as the most promising phytocompound to potentially inhibit the cellular entry and viral replication of SARS-CoV-2. We also conclude that compounds having oxygen atom at position 18 (C-ring), –OH group at position 19 (A-ring), and 6-C-glucoside attached to the A-ring at position 3 on a C(6)–C(3)–C(6) flavonoid scaffold could offer the best alternative to develop new leads against SARS-CoV-2.

The coronavirus disease 2019 pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) started in December 2019 from Wuhan, Hubei Province of China [1] . To date, COVID-19 is still an on-going pandemic and as of 2nd February, 2021, it has affected more than 102 million people globally and has killed more than 2 million people [2] . The clinical manifestation of COVID-19 includes shortness of breath or difficulty in breathing, fever, cough, headache, body ache, fatigue, sore throat, loss of taste or smell, nausea, vomiting, congestion or runny nose, and diarrhea [3] . Since the initial outbreak, scientists and researchers are working to develop therapies against SARS-CoV-2 [1] . Although a specific drug for COVID-19 is still not available, studies to repurpose existing drugs for COVID-19 is underway [4] .

With time, the pharmacological treatment options for COVID-19 had increased significantly [5] . However, problems like adverse effects, toxicity, or drug interactions seem to hinder the clinical utility of repurposed drugs for COVID-19 [4] . In the latest development, few vaccines have been authorized or approved for use against SARS-CoV-2 infection in humans [6] . Despite the breakthrough in vaccines, several deaths were reported among the recipients of those that receive the COVID-19 vaccinations [7] . Amidst the dilemma surrounding the available therapeutic options for COVID-19, phytotherapy may offer a safe and effective treatment against SARS-CoV-2 infection [8, 9] . Studies reveal that viral diseases were successfully treated with the bioactive compounds from medicinal plants [1, 10, 11] . Therefore, phytomedicines may be a promising prospect for COVID-19 therapy [1] .

To develop an effective antiviral therapy, it is logical to validate the possible drug targets by identifying the primary proteins involved in a viral replication process. Proteolytic activation occurs when a human protease 'furin' cleaves the spike protein (SP) of SARS-CoV-2 [12, 13] . After the SP of SARS-CoV-2 is cleaved by furin, SARS-CoV-2 uses the S1 subunit of the spike protein (SP) to bind to the angiotensin-converting enzyme 2 and then the S2 subunit of the SP to fuse with the host cell to release its viral RNA [14] . After the viral RNA is released into the host cell, SARS-CoV-2 uses the host cell machinery to translate polyproteins from the RNA genome. Finally, cleavage of the polyproteins by the main protease (M pro ) of SARS-CoV-2 results in replication and transcription of the viral genome [14] [15] [16] . The M pro of SARS-CoV-2 is crucial for the replication of the virus [15] . Therefore, the M pro of SARS-CoV-2 ( Fig. 1 ) is recognized as the most favorable drug target SARS-CoV-2 [15, 17, 18] . Cleavage of SARS-CoV-2 SP by furin is necessary for proteolytic activation, viral fusion, and viral entry into the host cell [12, 13] . As inhibition of furin would prevent the entry of SARS-CoV-2 into the host cell, furin ( Fig. 1) is also emerging as a favorable drug target in SARS-CoV-2 [12] .

In the quest to find a safe and effective antiviral therapy for COVID-19, bioactive compounds of medicinal plants have gained the attention of researchers and scientists [8, 9, 19] . The endeavor to find possible leads from medicinal plants against SARS-CoV-2 has already started [1, 20, 21] . Interestingly, several phytochemicals of medicinal plants were reported to elicit an anti-viral activity [10, 11] . The potential of flavonoids as an antiviral agent against SARS-CoV-2 including other respiratory diseases has been reported [22, 23] . Moreover, the anti-inflammatory and antioxidant properties of flavonoids may also be beneficial to alleviate the severity of SARS-CoV-2 infection [24] . Therefore, in the search for an alternative anti-SARS-CoV-2 therapy, exploring a medicinal plant that is rich in flavonoids might be helpful.

Acacia pennata (L.) Willd. (Family: Mimosaceae) is an important Southeast Asian medicinal plant that is rich in flavonoids [25] [26] [27] . A. pennata is distributed in China, Thailand, India, Sri Lanka, Myanmar, Bhutan, and Bangladesh [25] . The plant parts of Acacia species like the leaves, barks, roots, pods, twigs, gum, and flowers are traditionally used to treat various health ailments [28] . The bark and the root bark of A. pennata are traditionally used to treat respiratory ailments like bronchitis and asthma [25, 28] . Headaches and fevers are also treated with A. pennata [25] . The stem bark of A. pennata is used as a traditional anti-inflammatory agent [29] . The bark of A. pennata is also traditionally used to treat gastrointestinal ailments like cholera and dysentery [25, 29] . Pharmacological activities like antiviral [30, 31] , anti-inflammatory, antinociceptive [32] , antioxidant [25] , anti-parasitic [33] , antimicrobial [34] , and antidiabetic [35] had been scientifically investigated for A. pennata. The traditional utility of A. pennata to treat various health ailments may be attributed to the presence of flavonoids [26] , terpenoids [27] , phenols [28] , and saponins [30] in the plant. Also, the young shoot tips of A. pennata is consumed as a vegetable as it is rich in nutrients and minerals [36] [37] [38] .

A. pennata is widely available and is traditionally used to treat respiratory ailments. It is also used to treat other health ailments that are associated with COVID-19 such as fever, headache, inflammation, diarrhea like gastrointestinal issues such as dysentery, and cholera. The plant is also reported to possess anti-infective properties such as antiparasitic and antiviral activities. The facts provided above justifies the rationale to explore the phytocompounds of A. pennata for potential antiviral agents against SARS-CoV-2. Under urgent circumstances, computational approaches provide a great opportunity to identify natural compounds as potential inhibitors for cellular entry and viral replication of SARS-CoV-2. Therefore, the present study aims to identify potential leads against SARS-CoV-2 M pro and furin using computational studies on the phytocompounds isolated from A. pennata.

The research methodology of the current study is summarized in Fig. 2 .

An exhaustive online literature survey was conducted between November 2020 up to February 2021 to identify the phytocompounds that had been isolated from A. pennata. Academic databases like 'PubMed' , 'Science Direct' , and 'Google Scholar' were utilized for the literature survey. To search relevant papers, keywords such as 'phytochemistry' , 'phytoconstituents' , and 'phytocompounds' were used in combination with 'Acacia pennata' . The phytocompounds that had been isolated from A. pennata were identified from the retrieved literature.

The ligands used in the study were prepared manually. The chemical structures of the isolated phytocompounds 

The X-ray crystal structure of SARS-CoV-2 M pro (PDB ID: 6M2N) was downloaded in '.pdb' format from the RCSB-PDB website [40] . The co-crystal inhibitor (3WL) of SARS-CoV-2 M pro was downloaded in '.sdf ' format from the RCSB-PDB website [41] . The X-ray crystal structure of human furin (PDB ID: 4RYD) was downloaded in '.pdb' format from the RCSB-PDB website [42] . The co-crystal inhibitor (para-guanidinomethyl-Phac-R-Tle-R-Amba) of human furin was also downloaded in '.sdf ' format from the RCSB-PDB website [43] . The protein structures were pre-processed to remove water. The co-crystal inhibitors obtained along with the target proteins will be used as standard references.

An anticancer drug 'carmofur' was reported to inhibit SARS-CoV-2 M pro in Vero E6 cells (EC 50 = 24.30 ± 3.61 µM) [44] . Another study reported that naphthofluorescein inhibited furin that is responsible for the cleavage of SARS-CoV-2 SP [12] . Owing to their experimental evidence against SARS-CoV-2 M pro and furin, carmofur and naphthofluorescein will be used as a standard reference in addition to the co-crystal inhibitors. The chemical structures of carmofur and naphthofluorescein were prepared with MarvinSketch 20.10 software [39] . The correctness of their structures was checked using the 'Structure checker' add-on available on the MarvinSketch 20.10 software. Their structures were then saved in MDL SDfile ' * .sdf * .sd' format.

To predict the binding affinity between the target proteins and the phytocompounds, molecular docking simulation studies were carried out with Autodock Vina 1.1.2 on PyRx 0.8 virtual screening platform [45, 46] . In this study, chain A of SARS-CoV-2 M pro and chain A of furin were used as the drug targets [40, 42] . According to the standard protocols, the three-dimensional affinity grid box in the Vina search space of PyRx 0.8 tool should enclose the entire protein for blind docking [47, 48] . However, the protein-co-crystal inhibitor complex was used to manually identify the three-dimensional coordinates of the active binding sites for SARS-CoV-2 M pro and furin on the PyRx 0.8 tool. After the target protein was loaded onto the PyRx virtual screening platform, the target protein was expanded to reveal different chains that made up the protein. All the chains except chain A were removed from the scene for each protein. The protein-data bank format of the chain A of each protein was converted to 'pdbqt. ' format wherein hydrogens were added during the conversion process. The sequence of the amino acids including the co-crystal inhibitor were revealed by expanding chain A. The atoms of the cocrystal inhibitor were labelled to identify the accurate location of the co-crystal inhibitor that is present at the active binding site of the protein. The centre of the 3D affinity grid box (x = 25.0 Å; y = 25.0 Å; z = 25.0 Å) and the co-crystal inhibitor were adjusted to align so that the affinity grid box was able to cover all the amino acid residues at the active binding site of the protein. Then, the active binding site coordinates of SARS-CoV-2 M pro (x = − 33.1066; y = − 64.6263; z = 41.2995) and furin (x = 28.0166; y = 28.3695; z = − 6.0437) were manually adjusted for each docking process to simulate the active binding sites of the co-crystal ligands. The rest of the parameters such as energy minimization for protein and ligands were kept default. The molecular docking simulation studies were carried out as per the standard protocols of PyRx software for the phytocompounds and the standard references [47] . The phytocompounds with the top 5 binding affinities towards the active binding site of SARS-CoV-2 M pro and furin were selected for further analysis.

The drug score of the standard references and the phytocompounds with the top 5 binding affinities for SARS-CoV-2 M pro and furin were calculated using ORISIS Data Warrior v5.2.1 software [49] . The Data Warrior v.5.2.1 software takes several parameters (molecular properties, drug-likeness, and toxicity) into consideration to calculate the drug score. To select a phytocompound for further studies, it is important to make sure that the compound is safe and has good molecular properties. The Data Warrior v.5.2.1 software provides a low drug score for toxic compounds with undesirable molecular properties. Therefore, three phytocompounds with the highest drug scores were selected for further analysis.

Even though a phytocompound is found to have a high binding affinity towards a target protein, it is important to make sure that the phytocompound interacts with the amino acid residues at the active binding site of the protein. Protein-ligand interactions like hydrogen bonding, hydrophobic interactions, electrostatic interactions, and the interacting active site residues of the proteins were visualized for both the standard references and the three phytocompounds with the highest drug scores. Discovery Studio Visualizer v20.1.0.19295 software was used to visualize the 2-dimensional ligand interactions. The 3-dimensional binding pose of the phytocompounds and the standard references towards the target proteins were visualized using PyMOL molecular graphics system, Version 2.4.1 Schrodinger, LLC. [50] . The ligand interactions of the standard references and the potential leads were comparatively analyzed.

A compound can be effective as a drug only if it is bioavailable [51] . Since bioavailability issues can slow down the process of drug development [52] , preliminary investigation on the bioavailability of the phytocompounds is important. Synthetic accessibility is a fingerprint-based computational approach to determine the level of difficulty for synthesizing a compound [51] . As we aim to identify a potential lead from phytocompounds, it will be favourable for the identified lead to have high synthetic accessibility so that potent antiviral agents can be synthesized using the scaffold and pharmacophores of the lead compound. Therefore, after it was confirmed that multiple interactions occurred between the phytocompounds and the target protein, the bioavailability score and the synthetic accessibility of the phytocompounds were predicted with the SwissADME web tool [51] .

The structure-activity relationship of the phytocompound with the highest bioavailability score and the best synthetic accessibility was analyzed. A mapped structure of the phytocompound was prepared with MarvinSketch 20.10 software. The functional groups, chains, or atoms of the structure of the phytocompound that interacted with the amino acid residues at the active binding site were manually identified. The observations made in the study were compared with several other published data.

The list of phytocompounds that had been reported to be isolated from A. pennata is given in Table 1 . The chemical structures of all the isolated phytocompounds are also presented in Fig. 3 .

The binding affinities of the phytocompounds and the standard references (co-crystal inhibitors and additional standard references) towards SARS-CoV-2 M pro and furin are given in Table 2 .

The molecular properties (molecular weight and lipophilicity), drug-likeness, toxicity (mutagenicity, tumorigenic, reproductive effective, and irritant), and the overall drug score of the standard references and the phytocompounds with the top 5 binding affinities towards SARS-CoV-2 M pro and furin are given in Tables 3 and 4 respectively.

Ligand interactions of the three phytocompounds with the highest drug scores (C_18, C_4, and C_5) and the standard references (3WL and carmofur) with the amino acid residues at the active binding site of SARS-CoV-2 M pro in 2-dimensional view is given in Fig. 4 . The 3-dimensional binding pose of the phytocompounds and the standard references at the active binding pockets of SARS-CoV-2 M pro is given in Fig. 5 . The binding parameters of the standard references and the phytocompounds with the amino acid residues at the active binding site of SARS-CoV-2 M pro are given in Table 5 .

Ligand interactions of the three phytocompounds with the highest drug scores (C_18, C_4, and C_5) and the standard references (PGPRTRA and naphthofluorescein) with the amino acid residues at the active binding site of furin in 2-dimensional view is given in Fig. 6 . The 3-dimensional binding pose of the phytocompounds and the standard references at the active binding pockets of furin is given in Fig. 7 . The binding parameters of the standard references and the phytocompounds with the active binding site amino acid residues of furin are given in Table 6 .

The bioavailability score and synthetic accessibility of the three phytocompounds with the highest drug scores are given in Table 7 . A compound should have a high bioavailability score with a low numerical value for its synthetic accessibility (1 = easy to synthesize; 10 = difficult to synthesize) [51] .

The mapped structure of the phytocompound (C_5) with the highest bioactivity score and the best synthetic accessibility is given in Fig. 8 .

The hydroxyl (-OH) group at position 19 attached to the A-ring of the C 6 -C 3 -C 6 scaffold interacts with SARS-COV-2 M pro by forming two conventional hydrogen bonds with GLU166 at the active binding site of the protein (Figs. 4 and 8) . The -OH group at positions 28 and 31 of the 6-C glucoside attached to the A-ring of the C 6 -C 3 -C 6 scaffold interacts with SARS-CoV-2 M pro 

Literature survey reveals that 29 phytocompounds belonging to different classes of secondary metabolites like flavonoids, alkaloids, terpenoids, and saponins had been isolated from the stem, aerial parts, twigs, and leaves A. pennata (Table 1, Fig. 3 ). Recently, a team of researchers found that the hydromethanolic extract of A. pennata offers protection against DNA damage induced by acetaminophen in rat liver and they had also detected 5 new flavonoid-glycosides (Apigenin 6,8-di-C-hexoside, Apigenin-6-C-pentoside-8-C-hexoside, Apigenin-6-Chexoside-8-C-pentoside, Kaempferol 3,7-di-O-hexoside, and Luteolin-6-C-pentoside-8-C-hexoside) in the hydromethanolic extract [54] . However, these newly detected phytocompounds were not included in the study as there is still no report on their isolation. Another recent study reported that phosphodiesterase-5 was inhibited by the ethanolic extract of the leaves of A. pennata [55] . The same study highlighted four phytocompounds (quercetin-3-O-β-D-glucopyranoside, chrysin-7-O-β-D-glucopyranoside, koaburanin, and pinocembrin-7-O-β-D-glucopyranoside) as identified phytochemicals in the ethanolic extract [55] . The isolation of these four phytocompounds had been reported and are included in the present study. The majority of the isolated phytocompounds (n = 17) from A. pennata belonged to the class of flavonoids. The antiviral activity of flavonoids is well explored [56, 57] and their potential role against coronavirus infection had also been thoroughly discussed [24, 58] . Since A. pennata is rich in flavonoids, the chances of finding a potential lead from A. pennata as a SARS-CoV-2 inhibitor also increases. The study was initially proposed to limit the virtual screening to isolated flavonoids only. However, a literature survey revealed that a saponin '21β-O-[(2E)-6-hydroxyl-2,6dimethyl-2,7-octadienoyl] pitheduloside G' isolated from A. pennata exhibited an antiviral activity [30] . Therefore, in addition to flavonoids, the chances to find a potential lead from other classes of secondary metabolites also exist. Thus, irrespective of their chemical class, it was decided to screen all the reported phytocompounds that had been isolated from A. pennata.

The molecular docking simulation experiment revealed the binding affinities of the phytocompounds and the standard references (co-crystal inhibitors and additional standard references) towards SARS-CoV-2 M pro and furin (Table 2) . A high binding affinity signifies that a ligand will bind to the target protein with minimum energy [1] . Thus, a ligand with a high binding affinity towards the target protein is often desirable to be selected for further studies. In general, the binding affinities of the phytocompounds towards SARS-CoV-2 M pro ranged from − 8.0 to − 1.7 kcal/mol while the binding affinities of the phytocompounds towards furin ranged from − 9.0 to − 2.0 kcal/mol. When a molecular docking simulation study was conducted on a PyRx 0.8 software, a total of 9 poses were generated for each ligand [47] . Table 2 ). The co-crystal inhibitor '3WL' and the additional standard reference 'Carmofur' showed a binding affinity of − 6.7 kcal/mol and − 5.6 kcal/ mol respectively ( Table 2) Table 2 ). The co-crystal inhibitor 'Para-guanidinomethyl-Phac-R-Tle-R-Amba' (PGPRTRA) and the additional standard reference 'Naphthofluorescein' showed a binding affinity of − 7.0 kcal/ mol and − 10.0 kcal/mol respectively ( Table 2 ). The phytocompounds (C_29, C_28, C_1, C_26, C_14, C_23, C_4, C_5, and C_18) showed a better binding affinity towards furin than the co-crystal inhibitor (PGPRTRA) but not a better binding affinity than the additional standard reference (Naphthofluorescein).

The molecular properties, drug-likeness, toxicity parameters, and the overall drug score of the standard references and the phytocompounds with the top 5 binding affinities towards SARS-CoV-2 M pro and furin are given in Tables 3 and 4 respectively. The ORISIS Data Warrior v5.2.1 software uses molecular properties like molecular weight, cLogP, hydrogen acceptors, hydrogen donors, drug-likeness, and toxicity parameters like mutagenicity, tumorigenic, reproductive effective, and irritant to calculate the overall drug score of a compound. Toxicity is one of the reasons why many drugs had to be withdrawn from the market [59] . In addition to the binding affinity, the drug score can be used as a filter to narrow down and select suitable compounds for further studies. Thus, the drug score provides a referential value that can be used to judge the overall potential of a compound to qualify as a drug [60] . The reference standards against SARS-CoV-2 M pro had a drug score of 0.6981475 (3WL) and 0.16024 Table 3 Molecular properties, drug-likeness, toxicity parameters, and the overall drug score of phytocompounds with the top 5 binding affinities and the standard references towards SARS-CoV-2 M pro (Tables 3 and 4 ).

The drug score of all the three phytocompounds is higher than all the standard references except for the co-crystal inhibitor of SARS-CoV-2 M pro (3WL). As C_18, C_4 and C_5 were free from all possible signs of toxicity, they were subjected to further analysis to check whether they had interacted with the amino acid residues at the active binding site of the target proteins.

For SARS-CoV-2 M pro , the 2-dimensional interactions of the standard references (3WL and carmofur) and the phytocompounds (C_18, C_4, and C_5) with the amino acid residues at the active binding site is given in Fig. 4 . The 3-dimensional biding pose of the phytocompounds and the standard references at the active binding pockets of SARS-CoV-2 M pro is also given in Fig. 5 . In general, the interacting amino acid residues of SARS-CoV-2 M pro are GLN189, GLN192, MET165, PRO168, THR190, GLU166, ARG188, MET49, ASN142, HIS41, THR25, CYS44, GLY143 and THR26 (Fig. 4) . The binding parameters of the standard references and the phytocompounds with the amino acid residues at the active binding site of SARS-CoV-2 M pro are given in Table 5 . The number of conventional hydrogen bonds formed by C_18 (n = 3), C_4 (n = 5), and C_5 (n = 4) with the amino acid residues at the active binding site of SARS-CoV-2 M pro is more than the number of hydrogen bonds formed by the standard references '3WL' (n = 2) and 'carmofur' (n = 2) (Fig. 4) . Overall, C_18, C_4, and C_5 showed different types of interactions with 6, 7, and 4 amino acid residues at the active binding site of SARS-CoV-2 M pro respectively. Among the interacting amino acid residues of the active binding site of SARS-CoV-2 M pro , GLN189 has the highest occurrence (n = 3) followed by ARG188 (n = 2), ASN142 (n = 2), MET49 (n = 2), GLU166 (n = 2), GLN192 (n = 2), THR190 (n = 1), THR25 (n = 1), HIS41 (n = 1), CYS44 (n = 1), THR26 (n = 1), GLY143 (n = 1), MET165 (n = 1), and PRO168 (n = 1). For furin, the 2-dimensional interactions of the standard references (PGPRTRA and naphthofluorescein) and the phytocompounds (C_18, C_4, and C_5) with the amino acid residues at the active binding site is given in Fig. 6 . The 3-dimensional biding pose of the phytocompounds and the standard references at the active binding pockets of furin is also given in Fig. 7 . In general, the interacting amino acid residues of furin are ASP258, GLU299, ASN295, FMT601, SER368, HIS194, GLY229, LEU227, ASP191, THR365, ARG193, ARG197, ARG298 and HIS364 (Fig. 6) . The binding parameters of the standard references and the phytocompounds with the active binding site amino acid residues of furin are given in Table 6 . The number of conventional hydrogen bonds formed by C_18, C_4, C_5, PGPRTRA, and naphthofluorescein are 8, 3, 5, 7, and 2 respectively (Fig. 6) . Overall, C_18, C_4, and C_5 showed different types of interactions with 6, 5, and 5 amino acid residues at the active binding site of furin respectively. Among the interacting amino acid residues of the active binding site of furin, FMT601 has the highest occurrence (n = 4) followed by ASN295 (n = 3), SER368 (n = 3), THR365 (n = 3), HIS194 (n = 2), LEU227 (n = 2), ARG298 (n = 2), ASP258 (n = 1), GLU299 (n = 1), GLY229 (n = 1), ASP191 (n = 1), ARG193 (n = 1), ARG197 (n = 1), and HIS364 (n = 1). Table 5 Binding parameters of the phytocompounds and the standard references with the amino acid residues at the active binding site of SARS-CoV-2 M pro

Electrostatic interaction   1  C_18  THR190, ARG188, ASN142  MET49, GLN189  GLU166  -2  C_4  THR25, HIS41, CYS44, ASN142, THR26  MET49  GLY143  -3  C_5  ARG188, GLN189, GLU166  ---4  3WL  GLN189, GLN192  --MET165   5  Carmofur  GLN192 PRO168 -- A hydrogen bond is an important facilitator for proper binding between a protein and a ligand [61] . In a similar fashion with the reference standards, the phytocompounds (C_18, C_4, and C_5) showed a good number of conventional hydrogen bonding with the amino acid residues at the active binding site of the target proteins (SARS-CoV-2 M pro and furin) (Figs. 4 and 6). Hydrophobic interactions are noncovalent bonding interactions that are considered to be crucial for protein folding and protein-ligand interactions [62, 63] . Two phytocompounds (C_18 and C_5) showed hydrophobic interactions with the amino acid residues at the active binding site of SARS-CoV-2 M pro (Fig. 4) . Electrostatic interactions are also reported to be important for protein stability, function, flexibility, and folding [63] . Two phytocompounds (C_4 and C_5) showed electrostatic interactions with the amino acid residues at the active binding site of furin (Fig. 6) . Analysis of ligand interactions revealed that all the three phytocompounds (C_18, C_4, and C_5) interacted with different amino acid residues at the active binding site of the target proteins (SARS-CoV-2 M pro and furin).

Quercetin-3-O-α-L-rhamnopyranoside (C_18) has been reported to exhibit antiviral activity against the influenza A virus (H1N1) [64] . SARS-CoV-2 is a viral pathogen that primarily affects the respiratory system [65] . Pharmacological activities of C_18 like antioxidant, anti-obesity, and other activities have also been reported [66] [67] [68] . The antioxidant activity of C_18 might prove beneficial in SARS-CoV-2 infection as COVID-19 is associated with oxidative stress [65] . In addition to A. pennata, Quercetin-3-O-α-L-rhamnopyranoside is also found in other plants [64, [66] [67] [68] [69] [70] [71] . Kaempferol 3-O-α-Lrhamnopyranosyl-(1 → 4)-β-D-glucopyranoside (C_4) and isovitexin (C_5) were reported to inhibit cyclooxygenase (COX)-1 and COX-2 enzyme in a COX-1/COX-2-catalyzed prostaglandin biosynthesis assay [32] . COX enzymes are of clinical relevance as they are pro-inflammatory agents that are inhibited by anti-inflammatory drugs [72] . The anti-inflammatory activity of C_5 against lipopolysaccharide-induced neuroinflammation has also been reported [73] . COVID-19 is associated with inflammation [74] and some studies suggest the treatment of COVID-19 by inhibition of COX-2 [75] . Moreover, the antioxidant activity of C_5 has also been reported [76] .

In addition to A. pennata, isovitexin (C_5) had also been isolated from other plants [77] [78] [79] [80] [81] . The existing data of the phytocompounds suggests that they could be considered as potential leads. However, the bioavailability and the degree of difficulty to synthesize a compound must be evaluated for the promising phytocompounds before they could be considered as potential leads.

The process of drug development can be slowed down by bioavailability issues associated with a compound [52] . For a compound to be effective as a drug, it is important that a sufficient concentration of a compound is available within the systemic circulation for a specified period so that a compound can exert its pharmacological action on the body [51] . Just like toxicity, it is important to assess the bioavailability of a compound at an earlier stage in the drug development process to avoid unfruitful outcomes in the future. Therefore, the bioavailability of the phytocompounds (C_18, C_4, and C_5) was determined with the SwissADME web tool. The SwissADME web tool determines the bioavailability based on the molecular properties and lipophilicity of a compound by applying different principles such as Lipinski's rule of 5, Ghose filter, Veber filter, Egan filter and Muegge filter [51] . Among the phytocompounds, C_5 has the highest bioavailability score (0.55) while C_4 and C_18 each have a low bioavailability score (0.17) ( Table 7) . The bioavailability study revealed that C_5 has better bioavailability than the other phytocompounds.

During the process of virtual screening to identify promising lead compounds, it is preferable to filter out a non-toxic, biologically active compound having good bioavailability [51] . Along with this, the degree of difficulty to synthesize a compound is also a factor that needs to be taken into consideration while selecting the most promising compound [51] . In our study, we have assessed the degree of difficulty to synthesize a compound with the SwissADME web tool. Synthetic accessibility is a fingerprint-based computational approach that can be used to determine the level of difficulty for synthesizing a compound [51] . A compound with a synthetic accessibility score of 1 indicates that the compound can be easily synthesized while a synthetic accessibility score of 10 suggests that the compound will be very difficult to synthesize [51] . Among the phytocompounds, C_5 (4.99) has the best synthetic accessibility score followed by C_18 (5.19 ) and C_4 (6.35) ( Table 7) . The synthetic accessibility study revealed that C_5 will be easier to be synthesized in comparison to the other phytocompounds.

Based on the binding affinity towards the target proteins, molecular properties, drug-likeness, toxicity, ligand interactions, bioavailability and synthetic accessibility, C_5 (Isovitexin) has been found as the most promising phytocompound that could act as a potential lead. Although C_5 showed good ligand interactions with SARS-CoV-2 M pro and furin, it is important to make sure that C_5 interacted specifically with the active amino acid residues at the active binding site of the target proteins. Therefore, we decided to compare our experimental findings with several studies published by other researchers. From Figs. 4, 6, and 8, it can be observed that the same functional group located at different positions on C_5 interacted with different amino acid residues at the active binding site of SARS-CoV-2 M pro and furin. The structure-activity relationship of C_5 revealed that the -OH group at position 19 of the A-ring interacted with GLU166 of SARS-CoV-2 M pro (Figs. 4 and 8) and also interacted with ASN295 of furin (Figs. 6 and 8). Several in-silico studies reported that many phytocompounds either formed conventional hydrogen bonds or interacted with GLU166 of SARS-CoV-2 M pro [82] [83] [84] . An in-silico study also reported that several phytocompounds either formed conventional hydrogen bonds or interacted with ASN295 of furin [85] . The -OH group at position 28 of the 6-C glucoside attached to the A-ring of C_5 interacted with ARG188 of SARS-CoV-2 M pro (Figs. 4 and 8) and also interacted with HIS364 and THR365 of furin (Figs. 6 and 8). Interestingly, other phytocompounds have also been reported to interact with ARG188 of SARS-CoV-2 M pro [82, 86] . Also, several phytocompounds were reported to either form conventional hydrogen bonds or interact with HIS364 and THR365 of furin [85] .

The -OH group at position 31 of the 6-C glucoside attached to the A-ring of C_5 interacted with GLN189 of SARS-CoV-2 M pro (Figs. 4 and 8) . Other compounds were also reported to either form conventional hydrogen bonds or interacted with GLN189 of SARS-CoV-2 M pro [86, 87] . The oxygen atom at position 18 of the C-ring formed a conventional hydrogen bond with ASN295 of furin ( Figs. 6 and 8) . Several phytocompounds were also reported to either form a conventional hydrogen bond or interacted with ASN295 of furin [85] . The -OH group at position 27 of the 6-C glucoside attached to the A-ring formed a conventional hydrogen bond with THR365 [85] . The B-ring of the C 6 -C 3 -C 6 scaffold of C_5 formed electrostatic interaction with ARG298 of furin ( Figs. 6 and 8) . Several phytocompounds were also reported to interact with ARG298 of furin [85] . It was observed that phytocompounds that interacted with the same active site residues of the target proteins similar to what we reported for isovitexin were published as potential inhibitors of SARS-CoV-2 M pro and furin. After comparison with other studies, it was observed that C_5 has the potential to interact specifically with the active amino acid residues at the active binding site of SARS-CoV-2 M pro and furin. Isovitexin (C_5) is apigenin with a 6-C glucoside attached to position 3 (Fig. 8) of the A-ring [88] . Recently, a study showed that the methanol-trifluoroacetic acid leaf extract of Vitis vinifera was able to effectively inhibit the replication of SARS-CoV-2 in-vitro. From about 40 phenolic compounds present in the leaf extract of V. vinifera, isovitexin (C_5, apigenin-6-C glucoside) was among the most abundant phenolic compounds present in the leaf extract [89] . Although the study did not identify the phenolic compound responsible for the inhibition of SARS-CoV-2, the fact that isovitexin was among the major phytoconstituents present in the plant extract suggests that isovitexin might play an important role in the inhibition of SARS-CoV-2.

Interestingly, an in-silico study found apigenin as a potential inhibitor of SARS-CoV-2 M pro [83] . A molecular docking simulation-based in-silico study found apigenino-7-glucuronide as a potential inhibitor of SARS-CoV-2 M pro . Molecular dynamics simulation reveal stable conformation between the SARS-CoV-2 M pro and apigenino-7-glucuronide complex [84] . The in-silico findings reported by other researchers suggest a stable proteinligand interaction between apigenin (with or without its sugar moiety) with SARS-CoV-2 M pro .

Pharmacophores could be any part of the structure of a compound that interacts with a target protein to exert a pharmacological activity [90] . The oxygen atom at position 18 of the C-ring, the -OH group at position 19 of the A-ring, the -OH group at position 27, 28, and 31 of the 6-C-glucoside attached to the A-ring, and the B-ring could be considered as the pharmacophores in the structure of isovitexin as they had interacted with the active amino acid residues at the active binding site of SARS-CoV-2 M pro and furin (Figs. 4, 6, and 8 ). All these pharmacophores are attached to the basic flavonoid (C 6 -C 3 -C 6 ) scaffold. Therefore, we conclude that compounds having oxygen atom at position 18 attached to the C-ring, -OH group at position 19 attached to the A-ring, and the presence of a 6-C-glucoside that is attached to the A-ring of a C 6 -C 3 -C 6 scaffold (Fig. 8 ) could offer the best alternative to develop new leads against SARS-CoV-2.

Since SARS-CoV-2 utilizes its main viral protease 'M pro ' and human protease 'furin' to achieve cellular entry and viral replication, we have screened 29 phytocompounds reported to be isolated from A. pennata for potential leads using computational studies. Initially, our computational guided study revealed three flavonoids viz. quercetin-3-O-α-L-rhamnopyranoside, kaempferol 3-O-α-Lrhamnopyranosyl-(1 → 4)-β-D-glucopyranoside, and isovitexin as promising phytocompounds against SARS-CoV-2. However, quercetin-3-O-α-L-rhamnopyranoside and kaempferol 3-O-α-L-rhamnopyranosyl-(1 → 4)-β-Dglucopyranoside were not selected for further studies as they had bioavailability issues and were also found as difficult to synthesize. Based on binding affinity, molecular properties, drug-likeness, toxicity parameters, ligand interactions, bioavailability score, synthetic accessibility, structure-activity relationship and other supporting literary evidence, we found that isovitexin (apigenin-6-Cglucoside) was the most promising phytocompound to potentially inhibit the cellular entry and viral replication of SARS-CoV-2. Also, based on the structure-activity relationship, we conclude that compounds having oxygen atom at position 18 attached to the C-ring, -OH group at position 19 attached to the A-ring, and the presence of a 6-C-glucoside that is attached to the A-ring of a C 6 -C 3 -C 6 scaffold could offer the best alternative to develop new leads against SARS-CoV-2. As the evidence provided in our study is limited to in-silico data only, further investigations such as in-vitro studies are suggested to understand the complete inhibitory potential of the phytocompound against SARS-CoV-2. 

Recognition of natural products as potential inhibitors of COVID-19 main protease (Mpro): In-Silico evidences

World Health Organization (2021) WHO Coronavirus disease (COVID-19) Dashboard

Centers for Disease Control and Prevention (2020) CDC COVID-19 Symptoms of Coronavirus

Drug repurposing: new strategies for addressing COVID-19 outbreak

Treatment for COVID-19: an overview

Regulatory Affairs Professional Society (2021) COVID-19 vaccine tracker

COVID-19: Norway investigates 23 deaths in frail elderly patients after vaccination

Natural products' role against COVID-19

Natural products and their derivative against coronavirus: a review of the non-clinical and pre-clinical data

Phenolic phytochemical displaying SARS-CoV papain-like protease inhibition from the seeds of Psoralea corylifolia

Computational screening of medicinal plant phytochemicals to discover potent pan-serotype inhibitors against dengue virus

Furin inhibitors block SARS-CoV-2 spike protein cleavage to suppress virus production and cytopathic effects

Furin-mediated protein processing in infectious diseases and cancer

Structural and functional properties of SARS-CoV-2 spike protein: potential antivirus drug development for COVID-19

Structure of M pro from SARS-CoV-2 and discovery of its inhibitors

Expanding our understanding of the role polyprotein conformation plays in the coronavirus life cycle

Drug targets for corona virus: a systematic review

Preliminary bioinformatics studies on the design of a synthetic vaccine and a preventative peptidomimetic antagonist against the SARS-CoV-2 (2019-nCoV, COVID-19) coronavirus

Turning the tide: Natural products and naturalproduct-inspired chemicals as potential counters to SARS-CoV-2 infection

In silico screening of natural compounds against COVID-19 by targeting Mpro and ACE2 using molecular docking

Anthraquinone derivatives as an immune booster and their therapeutic option against COVID-19

Inhibition of SARS-CoV 3CL protease by flavonoids

Review of evidence available on hesperidin-rich products as potential tools against COVID-19 and hydrodynamic cavitation-based extraction as a method of increasing their production

Roles of flavonoids against coronavirus infection

A Review on the phytochemistry, medicinal properties and pharmacological activities of 15 selected Myanmar medicinal plants

Flavonoid glycosides from the aerial parts of Acacia pennata in Myanmar

Terpenoids and a flavonoid glycoside from Acacia pennata leaves as hedgehog/GLImediated transcriptional inhibitors

Antioxidant and free radical scavenging activities of Indian Acacias: Acacia leucophloea (Roxb.) Willd., Acacia ferruginea Dc., Acacia dealbata Link. and Acacia pennata (L.) Willd

Ethnomedicinal value of various plants used in the preparation of traditional rice beer by different tribes of Assam

A new saponin with anti-HIV-1 protease activity from Acacia pennata

Efficacy of Thai medicinal plant extracts against herpes simplex virus type 1 infection in vitro and in vivo

Flavonoids from Acacia pennata and their cyclooxygenase (COX-1 and COX-2) inhibitory activities

Efficacy and structural effects of Acacia pennata root bark upon the avian parasitic helminth, Raillietina echinobothrida

Antimicrobial, antioxidant and anticancer activities of Thai local vegetables

In vitro screening of phenolic compounds, potential inhibition against α-amylase and α-glucosidase of culinary herbs in Thailand

Total and soluble oxalate contents in Thai vegetables, cereal grains and legume seeds and their changes after cooking

Larvicidal and pupicidal activities of crude and fractionated extracts of Acacia pennata (L.) Willd

Subsp Insuavis shoot tips against Aedes aegypti (L.) (Diptera: Culicidae)

Folk medicinal plants of the Nagas in India

RCSB Protein Data Bank (2021) 6M2N SARS-CoV-2 3CL protease (3CL pro) in complex with a novel inhibitor

RCSB Protein Data Bank (2021) 3WL 5,6,7-trihydorxy-2-phenyl-4H-chromen-4-one

RCSB Protein Data Bank (2021) 4RYD X-ray structure of human furin in complex with the competitive inhibitor para-guanidinomethyl-Phac-R-Tle-R-Amba

RCSB Protein Data Bank (2021) PRD_001257 para-guanidinomethyl-phenylacetyl-Arg-(3-methylvaline)-Arg-(amidomethyl)benzamidine

Structural basis for the inhibition of SARS-CoV-2 main protease by antineoplastic drug carmofur

PyRx website (2021) Python Prescription Virtual Screening Tool

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

Small-molecule library screening by docking with PyRx

Protein-ligand blind docking using QuickVina-W with inter-process spatio-temporal integration

DataWarrior: an opensource program for chemistry aware data visualization and analysis

PyMOL (2021) The PyMOL Molecular Graphics System

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

Oral bioavailability: issues and solutions via nanoformulations

Bioactivity-guided separation of the active compounds in Acacia pennata responsible for the prevention of Alzheimer's disease

Acacia pennata L. leaves: chemical profiling and impact on DNA damage, alteration of genotoxicity-related genes expression and ROS generation in hepatic tissues of acetaminophen treated male rats

Medicinal plants as a potential source of Phosphodiesterase-5 inhibitors: a review

Antiviral activity of quercetin-3-β-O-D-glucoside against Zika virus infection

Antiviral activity of selected flavonoids against Chikungunya virus

Flavonoids with inhibitory activity against SARS-CoV-2 3CLpro

Mechanisms of drug toxicity and relevance to pharmaceutical development

html#: ~: text= The% 20drug% 20sco re% 20com bines% 20dru glike ness

Regulation of protein-ligand binding affinity by hydrogen bond pairing

A systematic analysis of atomic proteinligand interactions in the PDB

Close-range electrostatic interactions in proteins

Experimental validation and computational modeling of anti-influenza effects of quercetin-3-O-α-L-rhamnopyranoside from indigenous south African medicinal plant Rapanea melanophloeos

COVID-19 and oxidative stress

Inhibitory effect of quercetin-3-O-αα-L-rhamnopyranoside from Chamaecyparis obtuse on aldose reductase and sorbitol accumulation

Inhibitory characteristics of flavonol-3-O-glycosides from Polygonum aviculare L. (common knotgrass) against porcine pancreatic lipase

Quercetin-3-O-α-L-rhamnopyranoside derived from the leaves of Lindera aggregata (Sims) Kosterm. evokes the autophagy-induced nuclear factor erythroid 2-related factor 2 antioxidant pathway in human umbilical vein endothelial cells

Isolation and characterization of polyphenols from Euphorbia heterophylla L. (Euphorbiaceae) leaves

Quantitative and qualitative determination of flavonoids and phenolic acid derivatives from pericarp of hot pepper fruit cv. Bronowicka Ostra

Two new acylated flavonol glycosides from Mimosa pigra L. leaves sub-family Mimosoideae

Cyclooxygenase enzymes: regulation and function

Isovitexinmediated regulation of microglial polarization in lipopolysaccharideinduced neuroinflammation via activation of the CaMKKβ/AMPK-PGC-1α signaling axis

The trinity of COVID-19: immunity, inflammation and intervention

Yavuz E (2020) COX2 inhibition in the treatment of COVID-19: Review of literature to propose repositioning of celecoxib for randomized controlled studies

Isovitexin exerts anti-inflammatory and anti-oxidant activities on lipopolysaccharideinduced acute lung injury by inhibiting MAPK and NF-κB and activating HO-1/Nrf2 pathways

The glycosylated flavonoids vitexin, isovitexin, and quercetrin isolated from Serjania erecta Radlk (Sapindaceae) leaves protect PC12 cells against amyloid-β25-35 peptideinduced toxicity

Preparative separation of vitexin and isovitexin from pigeonpea extracts with macroporous resins

Isolation and purification of orientin and isovitexin from Thlaspi arvense Linn

Fast identification of radical scavengers from Securigera varia by combining 13 C-NMR-Based dereplication to bioactivity-guided fractionation

Separation and purification of four tannins from Potentilla parvifolia Fisch. (Rosaceae) using high-speed counter-current chromatography

Potential bioactive glycosylated flavonoids as SARS-CoV-2 main protease inhibitors: A molecular docking and simulation studies

Virtual screening, ADMET profiling, PASS prediction, and bioactivity studies of potential inhibitory roles of alkaloids, phytosterols, and flavonoids against COVID-19 main protease (M pro )

Promising phytochemicals of traditional Indian herbal steam inhalation therapy to combat COVID-19-an in silico study

Virtual screening by targeting proteolytic sites of furin and TMPRSS2 to propose potential compounds obstructing the entry of SARS-CoV-2 virus into human host cells

Unravelling lead antiviral phytochemicals for the inhibition of SARS-CoV-2 M pro enzyme through in silico approach

In silico screening of FDA approved drugs reveals ergotamine and dihydroergotamine as potential coronavirus main protease enzyme inhibitors

A review on the pharmacological effects of vitexin and isovitexin

Antiviral Activity of Vitis vinifera Leaf Extract against SARS-CoV-2 and HSV-1

The pharmacophore concept and its applications in computer-aided drug design

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

The authors are thankful to the Department of Pharmaceutical Sciences, Dibrugarh University for providing the necessary facilities to carry out the work.

Authors' contributions JHZ and HL developed the idea of the work. JHZ, NG and AS came up with the study design. JHZ carried out the works. NG, AS and DC analyzed and interpreted the results. JHZ and HL drafted the manuscript. NG, AS and DC provided additional inputs. HL carried out critical reviewing of the manuscript. All five authors wrote, read and approve the manuscript.

Not applicable.