Quercetin attenuates viral infections by interacting with target proteins and linked genes in chemicobiological models.

Medicinally active compounds in the flavonoid class of phytochemicals are being studied for antiviral action against various DNA and RNA viruses. Quercetin is a flavonoid present in a wide range of foods, including fruits and vegetables. It is said to be efficient against a wide range of viruses. This research investigated the usefulness of Quercetin against Hepatitis C virus, Dengue type 2 virus, Ebola virus, and Influenza A using computational models. A molecular docking study using the online tool PockDrug was accomplished to identify the best binding sites between Quercetin and PubChem-based receptors. Network-pharmacological assay to opt to verify function-specific gene-compound interactions using STITCH, STRING, GSEA, Cytoscape plugin cytoHubba. Quercetin explored tremendous binding affinity against NS5A protein for HCV with a docking score of - 6.268 kcal/mol, NS5 for DENV-2 with a docking score of - 5.393 kcal/mol, VP35 protein for EBOV with a docking score of - 4.524 kcal/mol, and NP protein for IAV with a docking score of - 6.954 kcal/mol. In the network-pharmacology study, out of 39 hub genes, 38 genes have been found to interact with Quercetin and the top interconnected nodes in the protein-protein network were (based on the degree of interaction with other nodes) AKT1, EGFR, SRC, MMP9, MMP2, KDR, IGF1R, PTK2, ABCG2, and MET. Negative binding energies were noticed in Quercetin-receptor interaction. Results demonstrate that Quercetin could be a potential antiviral agent against these viral diseases with further study in in-vivo models. Supplementary Information: The online version contains supplementary material available at 10.1007/s40203-022-00132-2.

Introduction

Throughout history, virus-borne infectious diseases have caused millions of deaths and continue to be a hazard to public health. Scientists all over the world have made tremendous efforts to identify, treat, and prevent viral infections. However, discovering novel medicines may be difficult because of the complex life cycle of viruses and their rapid genetic alterations. Despite the availability of vaccines for many viral infectious diseases, some of the deadliest viruses, such as influenza A, Hepatitis C, Dengue type 2, and Ebola virus, continue to cause outbreaks worldwide and have resulted in significant morbidity and mortality in humans through epidemics or pandemics (Mourya et al. 2019). Highly contagious viruses are increasing the demand for vaccine research, but vaccine production is a long and complicated process. However, vaccines and antiviral medicines are often expensive, making them unavailable to most of the world's population. Both have a significant negative impact on the human body.

Natural compounds may be the best alternative for developing novel drugs. Many researchers have been working to synthesize new antiviral drugs over the last few years (Badshah et al. 2021). Natural substances have been more popular in recent decades because of their wide range of biological functions and direct usage as medications (Badshah et al. 2021; Ninfali et al. 2020). Quercetin is a unique flavonoid component that is well known for its antiviral properties. It has little or no side effects compared to synthetic medications, which could be an effective treatment option for different viral infections (Ninfali et al. 2020).

The name Quercetin comes from the Latin word ''Quercetum'' which means ‘Oak Forest’ (Batiha et al. 2020). The chemical formula of quercetin is 2-(3,4-dihydroxy phenyl)-3,5,7-trihydroxychromen-4-one and C15H10O7 (Batiha et al. 2020). Quercetin belongs to the flavonol class of flavonoids found widely among vegetables and fruits, including berries, lovage, capers, cilantro, dill, apples, and onions (Anand David et al. 2016). The color of quercetin is yellow and entirely soluble in lipids and alcohol, slightly soluble in hot water but insoluble in cold water (Batiha et al. 2020).

Quercetin (Fig. 1) is an active therapeutic agent, which shows potential antiviral activity and inhibits viruses by targeting viral infections at multiple stages. The main molecular mechanisms of quercetin's antiviral activities are the suppression of viral neuraminidase, proteases, and DNA/RNA polymerases and the alteration of many viral proteins. Quercetin has shown antiviral activity towards a wide range of viruses, including hepatitis C virus (HCV) (Khachatoorian et al. 2012), Mayaro virus (dos Santos et al. 2014), influenza A virus (IAV) (Kim et al. 2010), Chikungunya virus (CHIKV) (Lani et al. 2015), and Epstein-Barr virus (EBV) (Lee et al. 2015). Quercetin has been shown to inhibit HCV through binding and inactivating the viral NS3 protease (Bachmetov et al. 2012). Quercetin prevents DENV-2 replication by inhibiting viral RNA polymerase and blocking Ebola viral infection via VP24 interferon-inhibitory function (Fanunza et al. 2020). This manuscript, the mechanism of quercetin against four viruses such as influenza A, Hepatitis C, Dengue type 2, and Ebola virus will be discussed.

Fig. 1.

2D structure of quercetin (Ashwini et al. 2017)

Materials and methods

Antiviral drug design of Quercetin by computer-aided approach

Molecular docking study

Docking tools

The molecular docking was performed using Schrodinger suites-Maestro 2021 (Maestro 2021). Pockdrug online server was used to predict the best binding pocket and probable drug ability. Discovery studio (v 4.1) was used for the visualization.

Ligand preparation

The chemical structure of quercetin was extracted from the PubChem repository (https://pubchem.ncbi.nlm.nih.gov/). The ligand was prepared by using the LigPrep tool, which was embedded in Schrödinger suite-Maestro v 11.1, where the following parameters were used as follows: neutralized at pH 7.0 ± 2.0 using Epik 2.2 and the OPLS3 force field were used for minimization (Md. Amjad Hossen et al. 2021).

Protein preparation

To unravel quercetin activity against different viruses, four major proteins following NS5A (PDB ID: 4CL1) for HCV, NS5 (PDB ID: 3EVG) for DENV-2, VP35 (PDB ID: 4IBK) for EBOV, and NP-specific inhibitor (PDB ID: 6J1U) for IAV, were selected for molecular docking study based on literature review (Yang et al. 2021). The crystal structure of these proteins has been collected from the RCSB Protein Data Bank (PDB) an online database (https://www.rcsb.org/) and their active binding sites were selected based on drug probability score by using an online tool PockDrug (http://pockdrug.rpbs.univ-paris-diderot.fr/cgi-bin/index.py?page=home). Preprocessing, optimization and minimization processes were done by using Protein Preparation Wizard. This process is included in Schrodinger suit-maestro (v11.1). The structures were optimized at pH 7.0, water molecules fewer than 3 H-bonds to non-waters were removed. Restrained minimization was done where heavy atoms are converged to RMSD of 0.30 Å on the implemented OPLS3 force field. Then the receptor grids were generated after selecting the best binding sites by using an online tool PockDrug (Hossen et al. 2021).

Glide ligand molecular docking

To evaluate antiviral effects, prepared ligand and proteins were interacted by using Schrodinger suite-maestro (v11.1). Consequently, the spreadsheets and structures were exported for further research. Discovery Studio (v. 4.1) software was used to visualize receptor-ligand binding interactions in three dimensions (Md. Amjad Hossen et al. 2021).

Identifying the quercetin-interacting genes associated with functional enrichment in human

Quercetin–target protein network construction

We used STITCH 5 (http://stitch.embl.de/, ver. 5.0) to identify target proteins related to quercetin (Szklarczyk et al. 2016). It calculates a score for each pair of protein-chemical interactions.. We used the “single item by name” module of STITCH 5 tool to enter the term “quercetin” into the STITCH 5 tool. Then, we selected the organism as “Homo sapiens” and the “medium required interaction score” ⩾0.4 was set in the “settings” module. We used the “maximum number of interaction” as 500 by using “custom value” for getting the more interacted target genes with quercetin. We also utilized the Swiss Target Prediction tool for performing a combination of similarity measurements based on known 2D and 3D chemical structures to predict the corresponding potential bioactive targets (probability > 0.2, http://www.swisstargetprediction.ch/) (Gfeller et al. 2014). The online tool “calculate and draw custom Venn diagrams” (http://bioinformatics.psb.ugent.be/webtools/Venn/) was used to identify commonly interacted genes. We entered the total genes found in STITCH tool and the Swiss Target Prediction tool into the Venn diagrams tool (http://bioinformatics.psb.ugent.be/webtools/Venn/) to identify the commonly interacted genes with quercetin. Finally, we identified the common genes that are interacted with quercetin not only in STITCH tool but also in the Swiss Target Prediction tool, which indicated the interaction of quercetin with target genes in both tools.

Construction of protein–protein interaction (PPI) network and top-ranked gene identification of the predicted genes

We constructed a PPI network of the predicted 39 common genes by search tool for the retrieval of interacting genes (STRING) database (https://string-db.org/cgi/input.pl; STRING-DB v11.0) (Szklarczyk et al. 2019). We used the “multiple proteins” module of the STRING tool to enter the 39 commonly found genes into the STRING tool. Then, we selected the organism as “Homo sapiens” and the “medium required interaction score” ⩾0.4 was set in the “settings” module. We exported all interacted information of PPI network by using “exports” module. Then, we inputted the STRING-based information into the Cytoscape tool 3.8.2 (Shannon et al, 2003) to construct the PPI network. The cytoHubba (Chin et al. 2014) a user-friendly interface in Cytoscape tool to explore the nodes in the PPI networks. We installed the cytoHubba in the Cytoscape tool by using “app manager” module. Then calculated the degree of interaction in the original PPI network by using the “target network” module in the cytoHubba. Based on the degree of interaction, the rank of the target proteins in the PPI network was identified using the “Hubba nodes” module of cytoHubba. The obtained protein interaction data of 39 target proteins were imported into Cytoscape 3.8.2 software for visualizing the PPI network (Shannon et al, 2003).

Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analyses of the target proteins

To identify the role of target proteins that interact with quercetin in gene functions and cellular signaling pathway, gene ontology and gene set enrichment analysis (GSEA) was employed (Subramanian et al. 2005). We inputted the commonly found 39 interacted target proteins into the GSEA by using the “Investigate Gene Sets” module. Then, the KEGG (Kanehisa et al. 2019) pathways significantly associated with the predicted genes were identified by selecting “KEGG: KEGG gene sets” in the GSEA. Moreover, we analyzed the gene ontology (GO) function enrichment of proteins (39 interacted target proteins with quercetin, Table 2) involved in the PPI network by selecting the “GO: Gene Ontology” module in the GSEA. The target proteins involved with the biological process (BP), cellular components (CC), molecular function (MF), and the KEGG pathways were also described. The FDR value < 0.05, calculated by the Benjamini–Hochberg method (Benjamini and Hochberg 1995) was considered statistically significant.

Table 2

Binding affinity and ligand efficiencies estimation of quercetin with different viruses including HCV, DENV-2, EBOV, and IAV via using MM-GBSA

Name of proteins with PDB ID	MM-GBSA (kcal/mol)	Ligand efficiency
NS5A (PDB ID: 4CL1) for HCV	– 50.31	– 0.285
NS5 (PDB ID: 3EVG) for DENV-2	– 33.36	– 0.245
VP35 (PDB ID: 4IBK) for EBOV	– 35.31	– 0.206
NP specific inhibitor (PDB ID: 6J1U) for IAV	– 41.69	– 0.316

Results

Molecular docking analysis

Quercetin demonstrated promising binding interaction against different viruses including HCV, DENV-2, EBOV, and IAV while the results have been shown in (Table 1) and (Fig. 2). Binding affinity and ligand efficiencies estimation of quercetin with different viruses including HCV, DENV-2, EBOV, and IAV via using MM-GBSA are summarized in Table 2. Firstly, the compound quercetin explored tremendous binding affinity against NS5A protein for HCV with a docking score of -6.268 kcal/mol (Fig. 2A). Besides, it also interacted with 5 conventional H-bond such as PRO 154 (2.78 Å), SER 157 (2.46 Å), GLN 158 (2.76 Å), CYS 111 (1.99 Å) along with single C-H bond PRO 160 (2.54 Å). At the same time, it interacted with two pi-alkyl including VAL 101 (4.20 Å), LYS 110 (5.05 Å), and one pi-pi stacked bond TYR 153 (5.55 Å) interactions. In addition, NS5 for DENV-2 was interacted (Fig. 2B) with quercetin which elucidated a docking score of -5.393 kcal/mol with a single water hydrogen bond HOH 546 (2.63 Å) and three conventional H-bond like LYS 110 (2.19 Å), as well as SO4904, showed two bands with their distance (1.62 and 1.89 Å). Conversely, the compound showed an unfavorable donor group (2.48 Å) but two pi anion bonds with ASP 146 (3.83 and 3.24 Å) and one conventional hydrogen bond GLY 83 (3.07 Å). Furthermore, quercetin was interacted with VP35 for EBOV remarkably (Fig. 2C) by four conventional H-bond following ALA 221 (2.45 Å), ARG 225 (2.82 Å), DMS 403 (3.06 Å), and VAL 294 (3.60 Å) whereas, docking score was 4.524 kcal/mol. The compound also demonstrated a single pi-alkyl bond ILE 295 (4.60 Å), two pi- cation bonds (4.64 and 4.97 Å) with LYS, and a pi-lone pair bond GLN 254 (3.00 Å). Moreover, the compound also showed a favorable binding score of -6.954 kcal/mol against NP for IAV while the compound was interacted with (Fig. 2D) by five H-bond such as ASP 51 (1.64 Å), SER 283 (1.68 and 2.61 Å), CYS 279 (2.68 Å) and THR 45 (1.83 Å) three pi alkyl bond following ALA 286 (4.41 and 4.16 Å), CYS 44 (4.86 Å)) and single amide pi stacked bond GLY 282 (4.22 Å). Hence, quercetin was unveiled as a potential agent based on the activity.

Table 1

Impacts of Quercetin against different viruses including HCV, DENV-2, EBOV, and IAV via molecular docking analysis

Name of proteins with PDB ID	Docking score	Glide emodel	Glide energy	Bond interaction types	Involved amino acid residues	Distances
NS5A (PDB ID: 4CL1) for HCV	– 6.268 kcal/mol	– 60.248 kcal/mol	–44.298 kcal/mol	Hydrogen bond	PRO 154	2.78 Å
					SER 157	2.46 Å
					GLN 158	2.76 Å
					CYS 111	1.99 Å
					PRO 160	2.54 Å
				Pi-pi stacked	TYR 153	5.55 Å
				Pi-alkyl	VAL 101	4.20 Å
					CYS 111	5.14 Å
					CYS 111	4.72 Å
					LYS 110	5.05 Å
NS5 (PDB ID: 3EVG) for DENV-2	– 5.393 kcal/mol	– 62.617 kcal/mol	– 47.496 kcal/mol	Water H-bond	HOH 546	2.63 Å
				Conventional H-bond	LYS 110	2.19 Å
					SO 4904	1.62 Å
				Pi anion bond	SO 490	1.89 Å
					ASP 146	3.83 Å
					ASP 146	3.24 Å
VP35 (PDB ID: 4IBK) for EBOV	4.524 kcal/mol	– 40.324	– 31.814	Conventional H-bond	ALA 221	2.45 Å
					ARG 225	2.82 Å
					DMS 403	3.06 Å
					VAL 294	3.60 Å
				Pi-alkyl bond	ILE 295	4.60 Å
				Pi- cation bond	LYS 248	4.64 Å
					LYS 248	4.97 Å
NP specific inhibitor (PDB ID: 6J1U) for IAV	– 6.954 kcal/mol	– 56.237	– 42.26	H-bond	ASP 51	1.64 Å
					SER 283	1.68 Å
					SER 283	2.61 Å
					CYS 279	2.68 Å
					THR 45	1.83 Å
				Pi alkyl bond	ALA 286	4.41 Å
					ALA 286	4.16 Å
					CYS 44	4.86 Å
				Amide pi stacked bond	GLY 282	4.22 Å

Fig. 2

Illustration of Quercetin binding mode 3D (left) and 2D (right) against different viruses including HCV (A), DENV-2 (B), EBOV (C), and IAV (D) via molecular docking analysis

Quercetin significantly interacted with various genes

By using the STITCH tool, we found that the quercetin interacted with the 299 genes (Fig. 3 and Supplementary Table S1). Besides, we identified the 100 genes that interacted with quercetin with a minimum probability of 0.20 (Fig. 3 and Supplementary Table S2). After doing the Venn diagram, we found that 39 genes that are identified by both tools interact with quercetin (Fig. 3 and Supplementary Table S3). The interaction of quercetin with 39 genes is shown in Fig. 4 (supplementary Tables S3).

Fig. 3

Identifying common 39 genes that are interacted with quercetin. Interaction of common 39 genes with quercetin

Fig. 4

Interaction of common 39 genes with quercetin. The Cytoscape tool was utilized to illustrate the interactions

Quercetin-associated genes are involved in protein–protein interaction

We inputted all 39 common genes into the STRING to investigate their involvement in the PPI network. Interestingly, we found that 38 genes, except PYGL, are involved in the PPI network with PPI enrichment p-value < 1.0e-16, number of edges 153, average node degree 7.85, and average local clustering coefficient is 0.589. The STRING-based information was inputted into the Cytoscape tool for constructing the PPI network. The PPI network of 38 genes is displayed in Fig. 5. The redder node indicated the highest number of degrees of interaction in the network and yellow node indicated the smaller number of interactions with other genes. Then, we used the cytoHubba (Chin et al. 2014) a user-friendly interface in Cytoscape tool for exploring the highly interconnected nodes in the PPI network. We identified the top interconnected nodes in the PPI network from the original PPI network by using the cytoHubba. The top 10 interconnected nodes (based on the degree of interaction with other nodes) are AKT1, EGFR, SRC, MMP9, MMP2, KDR, IGF1R, PTK2, ABCG2, and MET (Table 3), indicating that the cellular signaling controlled by these genes may be regulated by the quercetin.

Fig. 5

The protein–protein interaction (PPI) network of common 39 genes that are interacted with quercetin in STITCH and the Swiss Target Prediction tool

Table 3

Top ten hub nodes in the PPI

Rank	Symbol	Name	Degree of interaction
1	AKT1	AKT serine/threonine kinase 1	27
2	EGFR	epidermal growth factor receptor	23
3	SRC	"SRC proto-oncogene, non-receptor tyrosine kinase	22
4	MMP9	matrix metallopeptidase 9	16
5	MMP2	matrix metallopeptidase 2	14
5	KDR	kinase insert domain receptor	14
7	IGF1R	Insulin-like growth factor 1 receptor	12
7	PTK2	protein tyrosine kinase 2	12
9	ABCG2	ATP binding cassette subfamily G member 2	10
10	MET	"MET proto-oncogene, receptor tyrosine kinase	9

Functional enrichment analysis of quercetin associated genes

The functional analysis of gene ontology included the biological processes (Fig. 6, Supplementary Table S5), cellular components (Supplementary Table S6), and molecular functions (Supplementary Table S7). The significantly enriched top biological processes are a response to oxygen-containing compounds, cellular response to oxygen-containing compounds, regulation of cell death, cellular response to chemical stress, positive regulation of signaling, response to oxidative stress, apoptotic process, cell population proliferation, cellular response to reactive oxygen species, and response to nitrogen compound. The significant cellular components are mitochondrion, receptor complex, an extrinsic component of the cytoplasmic side of the plasma membrane, external side of the apical plasma membrane, an extrinsic component of membrane, plasma membrane region, an extrinsic component of the plasma membrane, apical plasma membrane, side of the membrane, and intrinsic component of the plasma membrane. Moreover, the top molecular functions are adenyl nucleotide binding, protein kinase activity, ribonucleotide binding, kinase activity, transferase activity transferring phosphorus-containing groups, identical protein binding, protein tyrosine kinase activity, oxidoreductase activity, transition metal ion binding, and transmembrane receptor protein tyrosine kinase activity.

Fig. 6

The top ten biological processes (BP), cellular components (CC), and molecular functions (MF) are significantly associated with quercetin-targeted genes

Furthermore, KEGG pathways enrichment analysis identifies the 50 significant terms associated with quercetin interacted genes (Table 4). The pathways mainly involved cancer, metabolism, cellular signaling cascades, and immune regulations. The cancer-associated pathways are pathways in cancer, prostate cancer, melanoma, endometrial cancer, acute myeloid leukemia, glioma, small cell lung cancer, bladder cancer, non-small cell lung cancer, colorectal cancer, pancreatic cancer, renal cell carcinoma, and chronic myeloid leukemia. The metabolic pathways are arachidonic acid metabolism, tryptophan metabolism, tyrosine metabolism, and metabolism of xenobiotics by cytochrome P450. The cellular signaling and developmental pathways are focal adhesion, ErbB signaling pathway, VEGF signaling pathway, progesterone-mediated oocyte maturation, adherens junction, Endocytosis Insulin signaling pathway, ABC transporters, steroid hormone biosynthesis, gap junction, GnRH signaling pathway, oocyte meiosis, cell cycle, neurotrophin signaling pathway, axon guidance, Jak-STAT signaling pathway, regulation of actin cytoskeleton, mTOR signaling pathway, p53 signaling pathway, apoptosis, melanogenesis, tight junction, calcium signaling pathway. The enriched immune regulatory pathways are leukocyte transendothelial migration, chemokine signaling pathway, epithelial cell signaling in Helicobacter pylori infection B cell receptor signaling pathway, cytokine-cytokine receptor interaction, T cell receptor signaling pathway, Fc epsilon RI signaling pathway, hematopoietic cell lineage, Fc gamma R-mediated phagocytosis, and Toll-like receptor signaling pathway. These findings indicated that the quercetin-associated interactions of genes are regulating the key biological functions.

Table 4

The significantly enriched KEGG pathways associated with genes that are interacted with quercetin

KEGG pathways	Genes in overlap	p-value	FDR
Pathways in cancer	11	1.08E-14	2.01E-12
Focal adhesion	9	2.71E-13	2.52E-11
ErbB signaling pathway	6	2.61E-10	1.40E-08
Prostate cancer	6	3.00E-10	1.40E-08
Melanoma	5	8.07E-09	3.00E-07
VEGF signaling pathway	5	1.14E-08	3.54E-07
Progesterone-mediated oocyte maturation	5	2.01E-08	5.35E-07
Endometrial cancer	4	1.96E-07	4.55E-06
Acute myeloid leukemia	4	2.85E-07	5.68E-06
Arachidonic acid metabolism	4	3.06E-07	5.68E-06
Glioma	4	4.85E-07	8.20E-06
Adherens junction	4	7.76E-07	1.20E-05
Endocytosis	5	8.78E-07	1.26E-05
Small cell lung cancer	4	1.36E-06	1.81E-05
Leukocyte transendothelial migration	4	4.95E-06	6.14E-05
Bladder cancer	3	9.36E-06	1.05E-04
Insulin signaling pathway	4	9.57E-06	1.05E-04
ABC transporters	3	1.08E-05	1.11E-04
Non-small cell lung cancer	3	2.01E-05	1.96E-04
Steroid hormone biosynthesis	3	2.12E-05	1.97E-04
Colorectal cancer	3	3.04E-05	2.70E-04
Chemokine signaling pathway	4	3.39E-05	2.86E-04
Epithelial cell signaling in Helicobacter pylori infection	3	4.02E-05	3.25E-04
Pancreatic cancer	3	4.38E-05	3.26E-04
Renal cell carcinoma	3	4.38E-05	3.26E-04
B cell receptor signaling pathway	3	5.39E-05	3.85E-04
Gap junction	3	9.28E-05	6.39E-04
Cytokine-cytokine receptor interaction	4	1.25E-04	8.32E-04
GnRH signaling pathway	3	1.31E-04	8.38E-04
T cell receptor signaling pathway	3	1.59E-04	9.88E-04
Oocyte meiosis	3	1.82E-04	1.09E-03
Cell cycle	3	2.45E-04	1.41E-03
Neurotrophin signaling pathway	3	2.51E-04	1.41E-03
Axon guidance	3	2.69E-04	1.47E-03
Jak-STAT signaling pathway	3	4.60E-04	2.45E-03
Tryptophan metabolism	2	6.95E-04	3.59E-03
Tyrosine metabolism	2	7.66E-04	3.85E-03
Regulation of actin cytoskeleton	3	1.16E-03	5.59E-03
mTOR signaling pathway	2	1.17E-03	5.59E-03
p53 signaling pathway	2	2.00E-03	9.28E-03
Metabolism of xenobiotics by cytochrome P450	2	2.11E-03	9.59E-03
Chronic myeloid leukemia	2	2.29E-03	1.02E-02
Fc epsilon RI signaling pathway	2	2.68E-03	1.16E-02
Apoptosis	2	3.24E-03	1.34E-02
Hematopoietic cell lineage	2	3.24E-03	1.34E-02
Fc gamma R-mediated phagocytosis	2	3.93E-03	1.59E-02
Melanogenesis	2	4.34E-03	1.71E-02
Toll-like receptor signaling pathway	2	4.42E-03	1.71E-02
Tight junction	2	7.28E-03	2.76E-02
Calcium signaling pathway	2	1.29E-02	4.80E-02

Discussion

The in vitro and in vivo findings are usually verified using the computational model in a drug discovery process. Molecular docking, a sort of computer study, is often used to predict ligand-target interactions and acquire a better understanding of natural product biological activity. It also outlines the likely mechanism(s) of action as well as the binding modalities within enzyme binding pockets (Rodenhuis-Zybert et al. 2010). Our anticipated compound Quercetin, a flavonoid, was found to be highly interacted with the receptor of four proteins such as NS5A (PDB ID: 4CL1) for HCV, NS5 (PDB ID: 3EVG) for DENV-2, VP35 (PDB ID: 4IBK) for EBOV, NP-specific inhibitor (PDB ID: 6J1U) for IAV while the bond-energy of quercetin with hepatitis C and influenza viruses are more negative than two others implying the higher binding affinity.

In the Big data era, molecular docking and network pharmacological studies are offering a great development opportunity in terms of algorithm development, theoretical analyses, and applications. However, researchers are concerned about the difficulty of combining massive clinical and experimental data with scientific verification to uncover the regulatory mechanisms of network pharmacology and conduct research more effectively. To avoid ambiguous usage of network pharmacology, data validation before the clinical application is an inevitable step (Hoofnagle, 2002). As a result, experimental validation is required before applying these findings to clinical practice. The PPI network is involved in numerous biological processes, and faulty PPI is at the root of many clinical disorders (Falck-Ytter et al. 2002). We discovered top 10 interconnected nodes (based on the degree of interaction with other nodes) are AKT1, EGFR, SRC, MMP9, MMP2, KDR, IGF1R, PTK2, ABCG2, and MET which are strongly linked to antivirals using the PPI of a targeted protein. The epidermal growth factor receptor (EGFR) is involved in the regulation of essential epithelial processes as well as transcription factor upregulation. IRF1 is first implicated in increased interferon-kappa production, which leads to STAT1 activation and further amplification of downstream interferon-induced genes, such as antiviral effectors and chemokines (Lulli et al. 2016). Matrix Metalloproteinase 9 exerts antiviral activity against respiratory syncytial virus. Persistent MMP9 signaling is associated with tissue remodeling (Dabo et al. 2015) and the progression of many lung diseases including chronic obstructive pulmonary disease (COPD), and asthma (Ventura et al. 2014), and airway infections (Lee et al. 2013).

In addition, we identified GO pathways enrichment analysis of target protein. The GO analysis indicates that the target proteins may bind with the plasma membrane, chromosome, chromatin, regulatory region nucleic acid, or/and cellular receptor of cells for mediating the process of metabolic, immunological, signaling, and/or other activities, to exert signaling and antiviral potential. Pathway analysis revealed some pathways which are involved in the attenuation of oxidative stress and cytokine production. Oxidative stress is found to upregulate the biological enzymes NADPH oxidase, xanthine oxidase, and dual oxidase which altogether defend the antioxidative enzymes (superoxide dismutase, glutathione peroxidase, catalase, etc.), and cause several metabolic and physiological disorders in viral infections while Quercetin reverses the phenomenon and enhances the healing process (Borchi et al. 2010). Moreover, antioxidative enzymes are controlling the PPI network of target proteins which are significantly associated with enriched pathways, suggesting that they are responsible for controlling these signaling networks and enriched pathways to control viral infections.

KEGG pathways are primarily involved in metabolism (cytochrome P450 xenobiotic metabolism, cytochrome metabolism of drugs, drug metabolism, other enzymes, biosynthesis of steroid hormones, ascorbate, and alternate metabolism, interconversion of pentose and glucuronate, porphyrin and chlorophyll metabolism, and starch and sucrose metabolism), immune regulation (Cytokine-cytokine receptor interaction, T cell receptor signaling pathway, natural killer cell-mediated cytotoxicity, adipocytokine signaling pathway, Toll-like receptor signaling pathway, allograft rejection, and systemic lupus erythematosus), neurological regulation (amyotrophic lateral sclerosis), cellular development (apoptosis and p53 signaling pathway), cellular signaling (NOD-like receptor signaling pathway, RIG-I-like receptor signaling pathway, and PPAR signaling pathway), molecular transportation (ABC transporters) and interconversion of pentose and glucuronate metabolism and PPAR signaling pathway (Shen et al. 2021). Our anticipated Quercetin molecule is found to be highly interacted with the JAK-STAT Pathway which is considered a novel target to tackle viral infections because the Janus kinase–signal transducer and activator of transcription (JAK-STAT) pathway is especially relevant due to its essential role in the regulation of local and systemic inflammation in response to viral infections, being, therefore, a putative therapeutic target. Therefore, the genes that interacted with Quercetin are responsible for controlling these signaling networks and enriched pathways (Ezeonwumelu et al. 2021). Based on the compounds–protein targets relationships, it is evident that Quercetin is acting on the multi-target genes and proteins which are substantial contributions of Quercetin to biological and physiological activities to manage viral infections.

Conclusion

Antiviral effects of Quercetin have been investigated in this research through computational model. Molecular docking study showed that Quercetin had a strong binding affinity against HCV and IAV among other studied viruses. Highest binding energy – 6.954 kcal/mol was noticed for IAV. In the network pharmacological study, 38 hub genes showed a top ten interconnected nodes were triggered by Quercetin leading to expedite its antiviral action. Results demonstrate that the Quercetin could be a potential antiviral agent through a comprehensive study in experimental antiviral models.

Below is the link to the electronic supplementary material.

Supplementary file1 (XLSX 48 kb)