Next Article in Journal
The Bacterial Nucleoid: From Electron Microscopy to Polymer Physics—A Personal Recollection
Next Article in Special Issue
Preparation and Evaluation of Nanofibrous and Film-Structured Ciprofloxacin Hydrochloride Inserts for Sustained Ocular Delivery: Pharmacokinetic Study in Rabbit’s Eye
Previous Article in Journal
Endoscopic Submucosal Dissection for Esophageal Cancer: Current and Future
Previous Article in Special Issue
The Role of Mu Opioid Receptors in High Fat Diet-Induced Reward and Potentiation of the Rewarding Effect of Oxycodone
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Life
Volume 13
Issue 4
10.3390/life13040893
life-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bioactive Compounds and Signaling Pathways of Wolfiporia extensa in Suppressing Inflammatory Response by Network Pharmacology

by
Juri Jin
1,†,
Md. Helal Uddin Chowdhury
2,†,
Md. Hafizur Rahman
3,
Ki-Young Choi
1,* and
Md. Adnan
3,*
1
Division of Future Agriculture Convergence, College of Agriculture and Life Science, Kangwon National University, Chuncheon 24341, Republic of Korea
2
Ethnobotany and Pharmacognosy Lab, Department of Botany, University of Chittagong, Chattogram 4331, Bangladesh
3
Department of Bio-Health Convergence, College of Biomedical Science, Kangwon National University, Chuncheon 24341, Republic of Korea
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Life 2023, 13(4), 893; https://doi.org/10.3390/life13040893
Submission received: 8 March 2023 / Revised: 18 March 2023 / Accepted: 21 March 2023 / Published: 27 March 2023
(This article belongs to the Special Issue New Trends in Pharmaceutical Science: 2nd Edition)
Browse Figures
Versions Notes

Abstract

:
Wolfiporia extensa (WE) is a medicinal mushroom and an excellent source of naturally occurring anti-inflammatory substances. However, the particular bioactive compound(s) and mechanism(s) of action against inflammation have yet to be determined. Here, we studied anti-inflammatory bioactive compounds and their molecular mechanisms through network pharmacology. Methanol (ME) extract of WE (MEWE) was used for GC-MS analysis to identify the bioactives, which were screened by following Lipinski’s rules. Public databases were used to extract selected bioactives and inflammation-related targets, and Venn diagrams exposed the common targets. Then, STRING and Cytoscape tools were used to construct protein-protein (PPI) network and mushroom-bioactives-target (M-C-T) networks. Gene Ontology and KEGG pathway analysis were performed by accessing the DAVID database and molecular docking was conducted to validate the findings. The chemical reactivity of key compounds and standard drugs was explored by the computational quantum mechanical modelling method (DFT study). Results from GC-MS revealed 27 bioactives, and all obeyed Lipinski’s rules. The public databases uncovered 284 compound-related targets and 7283 inflammation targets. A Venn diagram pointed to 42 common targets which were manifested in the PPI and M-C-T networks. KEGG analysis pointed to the HIF-1 signaling pathway and, hence, the suggested strategy for preventing the onset of inflammatory response was inhibition of downstream NFKB, MAPK, mTOR, and PI3K-Akt signaling cascades. Molecular docking revealed the strongest binding affinity for “N-(3-chlorophenyl) naphthyl carboxamide” on five target proteins associated with the HIF-1 signaling pathway. Compared to the standard drug utilized in the DFT (Density Functional Theory) analysis, the proposed bioactive showed a good electron donor component and a reduced chemical hardness energy. Our research pinpoints the therapeutic efficiency of MEWE and this work suggests a key bioactive compound and its action mechanism against inflammation.

1. Introduction

Inflammation is a crucial component of innate immunity’s protection against infections or tissue damage brought on by pathogen invasion, non-microbial stimuli, chemical stimulants, and toxicants, as well as by improper autoimmune reactions [1,2,3]. It may be distinguished by the presence of fluid accumulation and active cells, which are often symptoms of tissue degeneration [4]. Thus, the primary functions of inflammation are the elimination of pathogens and cell repair, but effective inflammation resolution and the restoration of homeostasis are essential for preventing inflammatory disorders [5,6]. Numerous variables, such as a shortage of antioxidants, vitamins, and anti-inflammatory substances such as zinc and selenium, might contribute to an insufficient inflammatory resolution. Reactive oxygen species (ROS) are formed during the inflammatory process as a consequence of the synthesis of cytokines and other pro-inflammatory mediators [7,8,9,10]. Hence, nonsteroidal anti-inflammatory drugs (NSAIDs) have been utilized as successful treatments for inflammatory disorders. However, these medications might cause serious gastrointestinal toxicity that is linked to thrombus formation, elevated blood pressure, and congestive heart failure [11,12,13]. Because of these drawbacks, there is a pressing need for innovative anti-inflammatory medicines originating from natural sources, notably plants and mushrooms, that have the ability to decrease the release of inflammatory mediators while also minimizing adverse side effects.
Mushrooms are valued functional foods used for their culinary values (taste and flavor) and medicinal (disease prevention) properties [14,15,16,17], notably hypocholesterolemic [18] and antiproliferative [19] effects. Wolfiporia extensa (Peck) Ginns (syn. Poria cocos F.A.Wolf) is a nutritious pharmaceutical mushroom, traditionally used to treat various ailments such as inflammation, acute gastroenteric catarrh, chronic gastritis, dizziness, edema, nausea, nephrosis, and gastric atony [20]. The primary metabolite from W. extensa, “Polysaccharide-II or 1,6-branched 1,3-alpha-D-galactan”, successfully delayed the production of the IFNγ activated IP-10 inflammation marker through regulating the expression of the IP-10 gene at the translational level while assuring human vascular endothelial cell (EC) toxicological safety, and also showed its efficiency as a new anti-inflammatory agent [21]. In mice, the colitis induced by TNBS was successfully alleviated by mitigating pro-inflammatory cytokines and boosting anti-inflammatory mediators in blood and colon-tissue by employing CMP33 from W. extensa [22]. Exploratory research concerning secondary metabolites revealed an anti-inflammatory action shown by ethanolic extract of W. extensa sclerotium on paw edema in acute and chronic stages [23]. Lee et al. deduced that six triterpenoids extracted from W. extensa sclerotia regulated NO and PGE2 levels in LPS-induced Raw 264.7 cells, downregulating iNOS and COX-2 expression [24,25]. Despite such evidence of W. extensa’s anti-inflammatory potential, there is no information on the cellular signaling mechanism of its therapeutic action or how it might exert anti-inflammatory effects to restore tissue damage or injury in individuals.
Hopkins first introduced Network Pharmacology in 2007 [26], an integrated drug discovery endeavor that includes system biology, bioinformatics and multidirectional pharmacology of the system [27]. This newly emerged process can pinpoint the molecular system of disease and the pathways of the biochemical framework of drugs at the molecular level [28,29]. It has changed the research paradigm from the “one target, one medication” model, to the newly postulated “multiple targets, multi-components” approach, by offering comprehensive composite target-compound and target-pathway networks to facilitate rationality and compatibility of medicines [30].
In this study, a comprehensive network pharmacology approach was employed to explore the potential molecular mechanisms of the bioactive compounds of WE supporting its anti-inflammatory action. To identify bioactives from WE, gas chromatography and mass spectroscopy (GC-MS) analyses were conducted. Drug-likeness parameters were additionally filtered; the targets related to filtered bioactives and inflammation were acquired by accessing public databases. Common targets among them were culled, generating a substantial network to find core targets and bioactives. Gene set enrichment analysis was employed to unearth the targets linked to the hub signaling pathway, biological process, cellular component, and molecular function prediction. Finally, the selected targets were subjected to molecular docking simulation and quantum parameter analysis (DFT—Density Functional Theory) for extensive affirmation of WE’s chemical reactivity to determine the most effective anti-inflammation candidate (Scheme 1).

2. Materials and Methods

2.1. Mushroom Collection, Identification, and Extraction

The mushrooms (Wolfiporia extensa) were collected from (Latitude: 36. 666700, Longitude: 128. 510729, Gyeongsangbuk-do, Republic of Korea) in December 2020. The collected mushrooms (200 g) were dried at ambient temperature (20–22 °C; for 7 days) and ground into a coarse powder using an automated grinder. The refined powder (30 g) was soaked in 300 mL of methanol (Daejung, Siheung, Republic of Korea). The extraction (repeated 3 times at room temperature) was carried out in a sealed bottle, with continuous shaking and stirring (for 5 days) using an electric shaker machine in order to increase the yield rate. The mixture was filtered (Whatman qualitative filter paper Grade 1) and evaporated using a vacuum evaporator (IKA, Staufen, Germany). The evaporated sample (MEWE) was dried using a hot water bath (IKA, Staufen, Germany) at 40 °C and preserved in a refrigerator (−4 °C) for GC-MS analysis.

2.2. GC-MS Analysis

MEWE was analyzed by GC-MS technique using the GC-MS system (Agilent 7890A, 5975C Agilent Technologies Inc., Santa Rosa, CA, USA) and a DB-5MS capillary column (30 m × 0.25 µm × 0.25 mm). The detailed protocol was described in our previous study [31].

2.3. Bioactive Compounds Filtration

The bioactive compounds (by GC-MS) were detected by the drug-likeness protocol “Lipinski’s rule of five” to obtain potentially active compounds with a drug-like character. Each compound was considered in terms of the absorption, distribution, metabolism and excretion (ADME) framework and the requirement for an oral bioavailability score > 0.50. Here, the essential pharmacokinetic factor is the oral bioavailability (OB) aspect of the ADME processes [32]. This analysis was performed through accessing an online tool, Swiss ADME [33]. SMILES notations for the compounds were obtained from PubChem (https://pubchem.ncbi.nlm.nih.gov/, accessed on 23 December 2021) database.

2.4. Extraction of Compound Associated Targets and Inflammatory Targets

Using ‘Homo sapiens’ mode, we compiled the target genes associated with filtered bioactives by entering their SMILES into the SEA (Similarity Ensemble Approach) (http://sea.bkslab.org/, accessed on 7 January 2022) and STP (Swiss Target Prediction) (http://www.swisstargetprediction.ch/, accessed on 9 January 2022) databases. In contrast, inflammation-related genes were gathered accessing DisGeNeT (https://www.disgenet.org/search, accessed on 18 January 2022) [34], Malacards (https://www.malacards.org/, accessed on 24 January 2022) [35] and OMIM (https://www.ncbi.nlm.nih.gov/omim, accessed on 5 February 2022) [36] databases. VENNY 2.1 (https://bioinfogp.cnb.csic.es/tools/venny/, accessed on 13 February 2022) displayed the common targets among MEWE bioactives and inflammatory target genes.

2.5. Common Targets Network Construction

To assess possible protein interactions, the common targets were entered into the String Database (https://string-db.org/, accessed on 25 February 2022) in ‘Homo sapiens’ mode, with a medium confidence level of 0.400. For better visualization, we imported the network into Cytoscape 3.8.2 software [37], and the whole network was analyzed in Cytoscape using the CytoHubba plugin and the following three algorithms: Maximal Clique Centrality (MCC) [38], Maximum Neighborhood Component (MNC), and Degree value [39]. Formulas are given below:
MCC(v) = ∑C∈S(v) (|C| − 1)!
where S(v) is the set of maximum cliques containing v and (|C| − 1)! is the sum of all positive integers that are smaller than |C|.
MNC(v) = |V(MC(v))|
where MC(v) is the G[N(v)]’s mostly linked component and G[N(v)] is the G’s induced subgraph by N(v).
Deg(v) = |N(v)|
where N(v) represents the connections of v’s neighbors, and v is their respective node.

2.6. Mushroom-Bioactives-Targets Network Construction

The bioactive compounds of MEWE, common inflammatory targets, and mushroom were loaded into Cytoscape 3.8.2 software to generate a graphical representation of the Mushroom-Compound-Target network. The merging function plugin Cytoscape was used to create this network. The Network Analyzer was used to evaluate network topology parameters. Nodes represent bioactives, targets, and mushrooms, and edges indicate their interaction. The frequency of connective neighbors of a node is referred to as its degree. The greater the number of linked edges in a node, the greater the impact [40].

2.7. Analysis of GO and KEGG Pathway Involvement within Common Targets

All common targets were submitted to the Annotation, Visualization, and Integrated Discovery (DAVID, https://david.ncifcrf.gov/tools.jsp, accessed on 10 March 2022) database for molecular functional annotation and KEGG pathway analysis to uncover their involvement in signal transduction. We chose OFFICIAL GENE SYMBOL as the identifier and Homo sapiens as the species for this enrichment analysis. In network pharmacology, the KEGG database has immense significance in illustrating targets involved in a disease’s molecular mechanism. The GO database depicts the biological descriptors of those targets, notably, Biological Process (BP), Cellular Component (CC) and Molecular Function (MF) [41]. A threshold value for GO terms and pathways enrichment was selected as p-value < 0.05. The FDR error control approach was used to correct the p-value, and the outcome was referred to as the Q value. The bubble plot map of KEGG pathways was graphically displayed utilizing Origin Pro 2021 to analyze the pathways.

2.8. Preparation of Ligand and Receptor Protein

Preferential ligands revealed from compound-target network research included conventional drugs such as Aspirin and Indomethacin; co-crystallized protein ligands were acquired in .sdf format using the PubChem chemical library. Furthermore, the metabolites were processed for molecular docking assays using the LigPrep program included in the Schrödinger suite-Maestro v12.5, adopting the previously reported techniques [42]. In addition, 5 receptor proteins of the hub signaling pathways targets with crystal structures had been selected for accessing in the RCSB Protein Data Bank (https://www.rcsb.org/, accessed on 19 March 2022) and the UniProt database (https://www.uniprot.org/, accessed on 19 March 2022), within which each of them have been accessible, i.e., TLR4 (PDB ID: 3UL7), EGFR (PDB: 5WB7), FLT1 (PDB: 3HNG), NOS3 (PDB: 1M9J) and NOS2 (PDB: 1NSI). The Schrödinger Suite-Maestro v12,5 embedded Protein Preparation Wizard tools have been configured once the 3D crystal structure have been found from the RCSB database following our previously described protocols [43,44].

2.9. Molecular Docking Simulation Using Glide

We utilized Glide tools plugged-in Schrödinger Suite-Maestro version 12.5 software for molecular docking investigations to identify receptor grids’ active site molecules (co-crystallized ligand site) [42]. During grid preparation, we employed default topological parameters such as the vdw (van der Waals) scaling factor 1.00, the OPLS3 force field and charge cut-off value 0.25 for individual 3D protein structures. At the plausible docking site, a cubic box of specific dimensions was set on kernel active site residues of macromolecules and the box size was given 14 Å × 14 Å × 14 Å co-ordinates. The docking experiments were subsequently implemented deploying the standard precision (SP) scoring algorithm of Glide, with each ligand noted individually in terms of the highest rating posture and docking score.

2.10. Quantum Chemistry of Frontier Molecular Orbitals

By following the Lee Yang Parr (B3LYP-D3) correlation functional approach at the 6–31G++ (d,p) level basis set, all compounds’ (hub and standard drugs) structural coordinates were thoroughly optimized using Jaguar panel of Maestro 12.5 software [45]. This optimized geometry also yielded frontier molecular orbital energies of HOMO (highest occupied molecular orbital) and LUMO (lowest occupied molecular orbital). In order to compute the HOMO-LUMO gaps of each chemical, LUMO energy was subtracted from the appropriate HOMO energy value. Depending on the energies of frontier HOMO and LUMO, the following formulae computed each compound’s hardness (ɳ) and softness (S) to hypothesize their chemical reactivity.
ɳ = (HOMO − LUMO)/2
S = 1/ɳ

3. Results

3.1. Chemical Composition of MEWE

Gas Chromatography-Mass Spectrometry (GC-MS) analysis revealed the presence of 27 bioactives in MEWE (Figure 1). The spectrometric data of identified bioactives, including RT (Retention Time), area (%), chemical formula and the names of bioactives, are displayed in Table 1. Several diversified chemical classes including Organooxygen compounds, Pyrans, Acyl halides, Isothiocyanates, Azolines, Boronic acid derivatives, Carboxylic acids and derivatives, Fatty acyls, Naphthalenes, Benzene and substituted derivatives, Pyrimidine nucleosides, Glycerolipids, Organonitrogen compounds, and Steroids and steroid derivatives were found in MEWE.

3.2. Bioactive Compounds Filtration

Following Lipinski’s rule, the bioactives which have a molecular weight not more than 500, hydrogen bond acceptor (HBA) not more than 10, hydrogen bond donor (HBD) not exceeding 5, Moriguchi octanol-water partition coefficient score not exceeding 4.15 and maintaining a standard ‘Abott Bioavailability Score’ of not more than 0.1 were considered as biologically active constituents. Following these filtration criteria, all bioactives (27) showed drug likeliness while not violating more than one of Lipinski’s assertions (Supplementary File S1: Table S1).

3.3. Common Targets of Bioactives Intersected from SEA and STP Databases

After filtration, bioactives-associated targets were extracted from the SEA and STP databases. Inputting bioactives-specific SMILES in the respective databases, a total of 608 targets from SEA and 730 targets from STP were gathered after removing duplication (File S2). Finally, Venn analysis revealed 284 common targets (File S3) between two public databases (Figure 2A).

3.4. Identification of Inflammation Targets within Disease Targets and 284 Common Targets

Three public databases, namely, DisGeNET, OMIM, and Malacards, yielded a total of 7283 targets associated with inflammation (File S4). The Venn diagram showed 42 common targets (File S5) when comparing the inflammatory targets (7283) and 284 bioactives-related targets (Figure 2B). However, the identified 42 common targets were connected to 21 out of the 27 bioactives. No targets were found in either database that connected with the remaining 6 bioactives, namely, 2,4-Dimethyl-1,3-cyclopentanedione; 3-Hydroxy-2,3-dihydromaltol; 2-Butynol; 2-Butenoyl chloride; alpha-Hydroxyisobutyric acid cyclic butaneboronate and N4 hydroxycytosine.

3.5. Network Construction of 42 Common Targets

The extracted 42 common targets were added to the STRING database to generate a PPI network to uncover possible mechanisms of MEWE anti-inflammatory action. The STRING network analysis indicated that each target in the network was related to the others via 42 nodes, 177 edges, and an average number of neighbors of 8.429, where network diameter and radius were 4 and 2, respectively. The Cytoscape network was used to explore key targets in the network, by means of three algorithms (MCC, MNC and Degree value) integrated with cytoHUBBA to increase network node precision and accuracy. Interestingly, TLR4, EGFR, and NOS3 were the top three essential targets in all algorithm-based networks studied (File S1: Figure S1). Table 2 shows the ranking of each gene present in the network according to the three algorithms.

3.6. Analysis of Mushroom-Compound-Targets Network

The mushroom-compound-targets network was constructed to evaluate the interconnection between mushrooms components and common inflammatory targets. The Network Analyzer applications in Cytoscape ascertained that 100 edges were bound to 64 nodes in the network; 42 nodes were common disease targets, 21 nodes referred to compounds and one node referred to the mushroom, viz. W.extensa (Figure 3). To detect key chemicals within the network, compounds that interacted with targets and mushrooms were assessed based on their relationship with adjacent targets, referred to as the ‘Degree Value’. Finally, N-(3-chlorophenyl)naphthylcarboxamide was exposed as the most active and potent metabolite in the network, which might protect against inflammation [File S1: Table S2].

3.7. Gene Ontology (GO) Analysis of 42 Common Targets

Gene ontology (GO) enrichment assessment of potential targets was performed to elucidate the crucial pharmacological processes as well as to corroborate the biological features (biological processes, chemical contents and molecular function) of the 42 inflammation targets. Here, the first ten functional keywords were picked based on gene percentages. File S1: Figure S2A demonstrates that the top 10 biological processes (BP) were the release of sequestered calcium ions into the cytosol, response to lipopolysaccharide, positive regulation of phosphatidylinositol 3-kinase signaling, negative regulation of blood pressure, inflammatory response, positive regulation of ERK1 and ERK2 cascade, activation of MAPKK activity, negative regulation of the apoptotic process, positive regulation of cell proliferation and positive regulation of cytosolic calcium ion concentration. The top 10 chemical contents (CC) were the integral component of the plasma membrane, plasma membrane, membrane raft, cell surface, an intrinsic component of the plasma membrane, nuclear envelope, extracellular space, caveola, and the extrinsic component of the cytoplasmic side of plasma membrane focal adhesion (File S1: Figure S2B). The molecular function (MF) output also indicated the presence of protein tyrosine kinase activity, enzyme binding, receptor binding, phosphatidylinositol-4,5-bisphosphate 3-kinase activity, protein phosphatase binding, drug binding, phosphatidylinositol 3-kinase binding, protein kinase binding, Ras guanyl-nucleotide exchange factor activity and RNA polymerase II transcription factor activity and ligand-activated sequence-specific DNA binding characteristics, which may contribute to the anti-inflammatory activity of MEWE (File S1: Figure S2C).

3.8. KEGG Pathway Enrichment Analysis of Identified 42 Common Targets

Data in File S1: Table S3 show that 21 of the 42 common targets were actively engaged in inflammation progression connected to 10 KEGG (Kyoto Encyclopedia of Genes and Genomes) annotation pathways (filtered at less than 0.05 p-value). In order to undertake the pathway enrichment analysis, the rich factor and corrected False Discovery Rate (FDR) value (Q-value) were considered. The rich factor shows the extent of the pathway enrichment with a significantly lower Q value. Figure 4 elucidates that the HIF-1 signaling pathway was the most enriched within the selected targets in conformity with prior filtration. Importantly, the Target-Pathway Network in Figure 5 provides intrinsic insight into targets and the pathway relationships. Details of the core pathway, namely the HIF-1 signaling pathway, are shown in Figure 6.

3.9. Docking Score Assessment of Key Bioactives and 5 HIF-1 Signaling Pathway Targets

The compound-target network analysis identified “N-(3-chlorophenyl)naphthylcarboxamide” as the best bioactive component. It was therefore investigated to evaluate its binding affinity with five targets of the HIF-1 signaling pathway linked to inflammatory function (Table 3). The co-crystallized ligands of macromolecules were re-docked to assess complex binding stability energy and compare disease reaction to a typical medication (Aspirin, Indomethacin) for further validation.
In terms of docking affinity, N-(3-chlorophenyl)naphthylcarboxamide connecting to TRP-356, VAL-336, PHE-353, and PRO-334 residues (forming one hydrogen and three hydrophobic bonds) of 1M9J exposed a substantial impact on the linkage with NOS3 (PDB ID: 1M9J) (Figure 7A). Compared to conventional medicines such as Aspirin (−5.262 Kcal/mol) and Indomethacin (−7.451 kcal/mol), the affinity of the interaction complex (bioactive-NOS3) was −7.236 kcal/mol.
An interface complex of N-(3-chlorophenyl)naphthylcarboxamide with NOS2 (PDB ID: 1NSI) was stabilized by the 2-H bond (CYS-200 and ARG-700) and 10 hydrophobic (TRP-194, PHE-369, MET-434, TRP-372, and CYS-200) and one pi-sulfur (MET-434) interactions. This interaction revealed an optimal docking energy of −8.85 kcal/mol as compared to the reference medicine Aspirin (−6.467 kcal/mol) and Indomethacin (−7.491 kcal/mol) (Figure 7B).
The binding between the bioactive and the active site of FLT1 (PDB ID: 3HNG) target demonstrated −8.88 kcal/mol energy by three H-bonds (ASP-1040, GLU-878, and CYS-1039), ten hydrophobic bonds (VAL-891, LEU-1013, ILE-1038, VAL-841, ALA-859, LYS-861, VAL-892, VAL-909, and LEU-882), two electrostatic bonds (LYS-861) and one pi-sulfur bond (LYS-1039). This energy was substantially lower than for the conventional medicine Aspirin and Indomethacin, with energy values of −6.794 kcal/mol and −5.538 kcal/mol, respectively (Figure 7C).
The TLR4 gene (PDB ID: 3UL7) exhibited the comparable binding affinity with the bioactive (−4.027 kcal/mol), through one hydrogen bond of LEU-119 residues, three hydrophobic bonds of HIS-148, PRO-145, and LEU-119 residues (Figure 7D), which was a contrast to standard medications Aspirin (−4.178 kcal/mol) and Indomethacin (−4.168 kcal/mol).
In the 5WB7 (EGFR target) scenario, interaction with the bioactive revealed that the active pocket included one H-bonding with HIS-409 and two H-phobic bondings with PHE-45 and ARG-29 residues (Figure 7E). The resulted docking score (−4.123 kcal/mol) was almost parallel to the competitive standards Aspirin (−4.182 kcal/mol) and Indomethacin (−4.707 kcal/mol). Details of binding information are given in File S1: Table S4.
However, this best bioactive “N-(3-chlorophenyl)naphthylcarboxamide” interaction with TLR4 and EGFR targets displayed almost equivalent binding affinity relative to its co-crystallized complex energy. Aside from that, the bioactive complex with FLT1 and NOS2 showed relatively low docking energy, unlike the respective co-crystallized ligands; values <−8 kcal/mol indicate significantly identical stability energy for any complex. However, the complex with the NOS3 target had greater affinity than the co-crystallized ligands.

3.10. DFT Evaluation of Key Compound and Standard Drugs

The DFT (Density Functional Theory) of key compounds and standard drugs was checked to explore their chemical reactivity with other species or targets. Generally, how a chemical molecule donates or accept its valence electrons to the ligand can be confirmed by the HOMO and LUMO levels. Here, the key bioactive “N-(3-chlorophenyl)naphthyl carboxamide” showed a considerable HOMO energy (−0.224 Kcal/mol) and appears to be a good electron donor constituent compared to the standard drugs, Indomethacin and Aspirin. The stability of a compound depends on the HOMO-LUMO energy gap, and lower energy is linked to a soft molecule. The key bioactive manifested a lower chemical hardness energy (0.08255 Kcal/mol) than standard drug Aspirin (File S1: Table S5). Figure 8 displays the frontier molecular orbitals localization pattern of the ground state (HOMO) and the first excited state (LUMO) of corresponding compounds.

4. Discussion

From the pathophysiological point of view, inflammatory disorders are complex, including a cascade of events that may result in severe sickness and include several proteins and pathways. Herbal treatments have long been used as an integral part of standard medicine practice due to their rich chemical components [46]. However, in most situations, the pharmacological mechanism of action of traditional medicines is still unknown. Network pharmacology, in this regard, provides a fresh perspective on the search for effective herbal substances against various diseases [47].
Through utilizing the GC-MS technique, a total of 27 bioactives from MEWE were screened. Of them, 21 out of 27 bioactives were potentially directly involved in the therapeutic efficiency of MEWE against inflammation. In addition, the compound-target network exposed 42 inflammatory targets intimately associated with 21 bioactives in the mechanism of inflammation. Among them, N-(3-chlorophenyl)naphthylcarboxamide was classified as a core essential bioactive in the network. In addition, the KEGG pathway enrichment analysis of 42 common targets disclosed that the BP, CC and MF activities interacted with 10 signaling pathways, among which 8 signaling pathways were closely related to the inflammation process and development. The causality of the 10 signaling pathways in inflammation are outlined here.
VEGF signaling pathway: VEGF activation induces angiogenesis in rheumatoid arthritis, which increases synovium nutrient flow, leukocyte motility, and cytokine release. Angiogenic factors make tumors more vascular, leading to a faster spread [48,49]. Platelet activation pathway: activated platelets produce IL-6, IL-8, IL-1, and TNF- α (insoluble versions) and regulate pro-inflammatory actions such as phagocytosis, leukocyte migration, and ROS generation. These mediators affect vascular inflammation, asthma, atherosclerosis, and rheumatoid arthritis [50]. Regulation of actin cytoskeleton: immunodeficiency or autoinflammatory diseases are linked to protein scrappiness (actin severing proteins, nucleation promoting factors, stabilizing protein of actin, de-polymerizing protein of actin and actin nucleators). These proteins subsequently intertwine with the actin cytoskeleton [51]. Ras signaling pathway: Ras activation produces pro-inflammatory cytokines, contributing to rheumatoid arthritis and vascular inflammation [52]. Pathways in cancer: intrinsic and extrinsic processes connect cancer with inflammation, activating transcription factors including NFKB, STAT-3, and HIF-1. These variables cause tumor cell growth. Thus, tumor-associated inflammation rises, halting immunological defenses [53]. Toxoplasmosis: parasite effectors may slow IFN-triggered toxoplasmacidal processes, causing toxoplasma-induced inflammation. Such effectors affect STAT3/6 (upregulated by ROP16), NFKB (upregulated by GRA15), and MAPK (induction by ROP38) signaling pathways that impact cytokine production [54]. PI3K-Akt signaling pathway: cytokine TNF-α induced phosphatidylinositol-3-kinase and its downstream target Akt stimulation lead to the phosphorylation of IKK, which activates NFKB, and subsequently induces vascular diseases [55]. HIF-1 signaling pathway: hypoxia-activated NFkB stimulates the synthesis of pro-inflammatory cytokines and growth factors via HIF signaling pathways, causing hyperglycemia, cancer, atherosclerosis, and rheumatoid arthritis [56]. Rap1 signaling pathway: macrophage Rap1 enhanced IKB and p65 phosphorylation, allowing NFkB binding to DNA kB sites and influencing pro-inflammatory gene transcription. Rap1 in inflammatory macrophages may promote atherosclerosis [57]. Calcium signaling pathway: channels in endothelial cells are opened when inflammatory mediators such as vasoactive amines, peptides, protease thrombin, and eicosanoids interact with receptors on these cells [58].
These 10 signaling pathways were shown to have direct involvement in the initiation of inflammation. However, there was a much higher enrichment level for the HIF-1 signaling pathway than for other similar rich factors. Significantly, the largest enrichment occurs at the highest rich factor [59]. In a mechanistic sense, the hypoxia caused a significant reduction in the activity of the HIF hydroxylase enzyme. In addition, HIF activation by TLRs and the EGFR downstream channel stabilized HIF, allowing dimerization with HIF-1 and attachment to p300 co-activators following nucleus accumulation. In addition, hypoxia-induced downregulation of IKK2 results in phosphorylation and degradation of IkB and activation of NFkB. The sensitivity of several variables involved in angiogenesis, nitric oxide synthase, and inflammation is influenced by HIF signaling because of this enriched and downstream route [60]. Ultimately, “TLR4” controls the synthesis of inflammatory cytokines via regulating the activation of the NFkB pathway’s RelA/p50 transcription factor complexes [6]. Consistently activated NFkB perpetuates an invasive phenotype by upregulating cell cycle regulators, anti-apoptotic, proteolytic factors, and pro-inflammatory cytokines [53]. The EGFR gene also activates downstream pathways, including the MAPK and PI3K-Akt signaling pathways that regulate HIF-1 upon hypoxia. Nitric oxide (NO) production caused by the activation of nitric oxide synthases (NOS2) triggers several pre-cancerous and malignant lesions, such as Barrett’s mucosa [61]. In addition, the aberrant expression of the iNOS/eNOS enzymes can induce inflammation-related cardiomyocyte mortality and protein nitration disorder [62]. The elevated expression of FLT1 may also pertain to the development of rheumatoid arthritis inflammation [63]. These findings support the hypothesis that the HIF-1 signaling pathway is associated with inflammatory illnesses such as rheumatoid arthritis, inflammatory bowel disease, chronic renal disease, atherosclerosis, and asthma [64]. As a result, inhibiting NFKB activation and HIF-1’s downstream pathways may be a viable therapeutic strategy for treating inflammation (Figure 9). Five targets in the HIF-1 signaling pathway were rigorously docked with the essential molecule discussed here to ascertain the sensitivity of this pathway.
The docking assay concluded that FLT1, NOS3 and NOS2 demonstrated superior binding interaction to the key essential component N-(3-chlorophenyl)naphthyl carboxamide compared to either reference medicine (Aspirin, Indomethacin). TLR4 and EGFR transcribed protein’s binding energy and stability were almost identical to Aspirin, Indomethacin and co-crystallized ligands. Conversely, FLT1 and NOS2 had indistinguishable compound docking complex stability in response to corresponding co-crystallized docking affinity, but NOS3 had the best stability. Quantum chemical analysis at the DFT (Density Functional Theory) level has confirmed the compound’s chemical reactivity. The kinetic stability and chemical reactivity of a molecule are largely determined by its HOMO and LUMO energy gaps [65]. The high kinetic stability and low chemical reactivity may be attributed to a sizable HOMO-LUMO gap which helps to explain the chemical function descriptors like hardness and softness [66,67,68]. The softness of a chemical boosts the reactivity and our explored hub compound “N-(3-chlorophenyl)naphthyl carboxamide” demonstrated a greater degree of softness energy than Aspirin, suggesting that the compound had strong binding affinity to the targets. Overall, these findings imply that N-(3-chlorophenyl)naphthyl carboxamide could block key molecular targets that support inflammatory escalation.

5. Conclusions

In conclusion, we used a computer-assisted network pharmacology prediction to explore the molecular pathway process by which MEWE can act against inflammation, in order to learn more about how the HIF-1 signaling pathway and other important pathways influence inflammation. The docking simulation showed that “N-(3-chlorophenyl)naphthyl carboxamide” successfully inhibited hypoxia-induced HIF-1α activation through downregulation of downstream NFKB, MAPK, mTOR, and PI3K-Akt signaling channels related to inflammation in the HIF-1 signaling pathway. The low HOMO-LUMO energy gap of the compound confirmed its robust chemical reactivity behavior with potential targets. In spite of this, more pharmacodynamic and mechanistic research studies are required in order to obtain a complete understanding of the intricate synergistic activities that unite the pharmacological effectiveness of WE on inflammation, as described in this work.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/life13040893/s1, File S1: Table S1: Drug-likeness properties of 27 bioactives from methanolic extract of W. extensa. Table S2: Bioactives ranking order in Bioactives-Targets network based on degree value. Table S3: Inflammation related targets enrichment in 10 signaling pathways. Table S4: Details binding interaction of key compounds with HIF-1 signaling pathways genes. Table S5: Different quantum parameters of key compound of WE and standard drugs. File S1: Figure S1: Protein-Protein Interaction using different algorithm. File S1: Figure S2: Gene Ontology (GO) analysis of common targets between bioactives and inflammation. File S2: Compound related gene. File S3: Common genes between SEA and STP. File S4: Inflammatory disease related genes. File S5: Common genes between inflammation related targets and compound related overlapping genes.

Author Contributions

J.J. and M.H.U.C.: Methodology, Formal analysis, Investigation, Software, Data Curation, Writing—Original Draft; M.H.R.: Validation, Software, Data Curation, Writing—Review and Editing; K.-Y.C.: Writing—Review and Editing, Visualization, Supervision; M.A.: Conceptualization, Resources, Visualization, Supervision, Project administration. All authors have read and agreed to the published version of the manuscript.

Funding

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

All data are provided in this paper and supplementary materials.

Acknowledgments

All authors are thankful to the Department of Bio-Health Convergence, Kangwon National University, Chuncheon 24341, Republic of Korea.

Conflicts of Interest

The authors have declared no conflict of interest. They have no known competing financial interests or personal relationships that could have appeared to influence the research reported in this publication.

Abbreviations

SEASimilarity Ensemble Approach
STPSwiss Target Prediction
OMIMOnline Mendelian Inheritance in Man
MEWEMethanolic Extract of Wolfiporia extensa
HIFHypoxia Inducible Factor
MNCMaximum Neighborhood Component
MCCMaximal Clique Centrality
KEGGKyoto Encyclopedia of Genes and Genomes
FDRFalse discovery rate
PPIProtein-protein interaction
SMILESSimplified molecular input line entry system
NFkBNuclear factor kappa B
EGFREpidermal Growth Factor Receptor
TLR4Toll Like Receptor 4
FLT1fms related tyrosine kinase 1
NOS2Nitric Oxide Synthase 2
NOS3Nitric Oxide Synthase 23

References

  1. Ilic, N.M.; Dey, M.; Poulev, A.A.; Logendra, S.; Kuhn, P.E.; Raskin, I. Anti-inflammatory activity of grains of paradise (Aframomum Melegueta Schum) extract. J. Agric. Food Chem. 2014, 62, 10452–10457. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Xu, Q.; Qiao, Y.; Zhang, Z.; Deng, Y.; Chen, T.; Tao, L.; Xu, Q.; Liu, J.; Sun, W.; Ye, Y.; et al. New Polyketides with Anti-Inflammatory Activity from the Fungus Aspergillus rugulosa. Front. Pharmacol. 2021, 12, 1–13. [Google Scholar] [CrossRef] [PubMed]
  3. Medzhitov, R. Inflammation 2010: New Adventures of an Old Flame. Cell 2010, 140, 771–776. [Google Scholar] [CrossRef] [Green Version]
  4. Chauhan, A.S.; Diwan, P.V.; Jain, N.K.; Tomalia, D.A. Unexpected in vivo anti-inflammatory activity observed for simple, surface functionalized poly(amidoamine) dendrimers. Biomacromolecules 2009, 10, 1195–1202. [Google Scholar] [CrossRef] [PubMed]
  5. Wu, B.C.; Skovbakke, S.L.; Masoudi, H.; Hancock, R.E.W.; Franzyk, H. In Vivo Anti-inflammatory Activity of Lipidated Peptidomimetics Pam-(Lys-βNspe)6-NH2 and Lau-(Lys-βNspe)6-NH2 against PMA-Induced Acute Inflammation. Front. Immunol. 2020, 11, 1–11. [Google Scholar] [CrossRef]
  6. Chen, L.; Deng, H.; Cui, H.; Fang, J.; Zuo, Z.; Deng, J.; Li, Y.; Wang, X.; Zhao, L.; Chen, L.; et al. Inflammatory responses and inflammation-associated diseases in organs. Oncotarget 2017, 9, 7204–7218. [Google Scholar] [CrossRef] [Green Version]
  7. Taofiq, O.; Calhelha, R.C.; Heleno, S.; Barros, L.; Martins, A.; Santos-Buelga, C.; Queiroz, M.J.R.P.; Ferreira, I.C.F.R. The contribution of phenolic acids to the anti-inflammatory activity of mushrooms: Screening in phenolic extracts, individual parent molecules and synthesized glucuronated and methylated derivatives. Food Res. Int. 2015, 76, 821–827. [Google Scholar] [CrossRef] [Green Version]
  8. Kanwar, J.R.; Kanwar, R.K.; Burrow, H.; Baratchi, S. Recent advances on the roles of NO in cancer and chronic inflammatory disorders. Curr. Med. Chem. 2009, 16, 2373–2394. [Google Scholar] [CrossRef]
  9. Moro, C.; Palacios, I.; Lozano, M.; D’Arrigo, M.; Guillamón, E.; Villares, A.; Martínez, J.A.; García-Lafuente, A. Anti-inflammatory activity of methanolic extracts from edible mushrooms in LPS activated RAW 264.7 macrophages. Food Chem. 2012, 130, 350–355. [Google Scholar] [CrossRef]
  10. Taofiq, O.; Martins, A.; Barreiro, M.F.; Ferreira, I.C.F.R. Anti-inflammatory potential of mushroom extracts and isolated metabolites. Trends Food Sci. Technol. 2016, 50, 193–210. [Google Scholar] [CrossRef] [Green Version]
  11. Gunawardena, D.; Bennett, L.; Shanmugam, K.; King, K.; Williams, R.; Zabaras, D.; Head, R.; Ooi, L.; Gyengesi, E.; Münch, G. Anti-inflammatory effects of five commercially available mushroom species determined in lipopolysaccharide and interferon-γ activated murine macrophages. Food Chem. 2014, 148, 92–96. [Google Scholar] [CrossRef] [PubMed]
  12. Jayasuriya, W.J.A.B.N.; Handunnetti, S.M.; Wanigatunge, C.A.; Fernando, G.H.; Abeytunga, D.T.U.; Suresh, T.S. Anti-Inflammatory Activity of Pleurotus ostreatus, a Culinary Medicinal Mushroom, in Wistar Rats. Evid. Based Complement. Altern. Med. 2020, 2020, 6845383. [Google Scholar] [CrossRef] [Green Version]
  13. Lee, S.S.; Tan, N.H.; Fung, S.Y.; Sim, S.M.; Tan, C.S.; Ng, S.T. Anti-inflammatory effect of the sclerotium of Lignosus rhinocerotis (Cooke) Ryvarden, the Tiger Milk mushroom. BMC Complement. Altern. Med. 2014, 14, 1–8. [Google Scholar] [CrossRef] [Green Version]
  14. Kalač, P. A review of chemical composition and nutritional value of wild-growing and cultivated mushrooms. J. Sci. Food Agric. 2013, 93, 209–218. [Google Scholar] [CrossRef] [PubMed]
  15. Dore, C.M.P.G.; Alves, M.G.d.C.F.; Santos, M.D.G.L.; Souza, L.A.R.D.; Baseia, I.G.; Leite, E.L. Antioxidant and Anti-Inflammatory Properties of an Extract Rich in Polysaccharides of the Mushroom Polyporus dermoporus. Antioxidants 2014, 3, 730–744. [Google Scholar] [CrossRef]
  16. Yuan, B.; Zhao, L.; Rakariyatham, K.; Han, Y.; Gao, Z.; Kimatu, B.M.; Hu, Q.; Xiao, H. Isolation of a novel bioactive protein from an edible mushroom Pleurotus eryngii and its anti-inflammatory potential. Food Funct. 2017, 8, 2175–2183. [Google Scholar] [CrossRef] [PubMed]
  17. Souilem, F.; Fernandes, Â.; Calhelha, R.C.; Barreira, J.C.M.; Barros, L.; Skhiri, F.; Martins, A.; Ferreira, I.C.F.R. Wild mushrooms and their mycelia as sources of bioactive compounds: Antioxidant, anti-inflammatory and cytotoxic properties. Food Chem. 2017, 230, 40–48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Morales, D.; Tabernero, M.; Largo, C.; Polo, G.; Piris, A.J.; Soler-Rivas, C. Effect of traditional and modern culinary processing, bioaccessibility, biosafety and bioavailability of eritadenine, a hypocholesterolemic compound from edible mushrooms. Food Funct. 2018, 9, 6360–6368. [Google Scholar] [CrossRef] [Green Version]
  19. Haque, M.A.; Reza, A.A.; Nasrin, M.S.; Rahman, M.A. Pleurotus highking mushrooms potentiate antiproliferative and antimigratory activity against triple-negative breast cancer cells by suppressing Akt signaling. Integr. Cancer Ther. 2020, 19, 1534735420969809. [Google Scholar] [CrossRef]
  20. Park, W.; Joo, S.; Park, K.; Chang, Y.; Kim, C. Effects of the Geiji-Bokryung-Hwan on carrageenan-induced inflammation in mice and cyclooxygenase-2 in hepatoma cells of HepG2 and Hep3B. Immunopharmacol. Immunotoxicol. 2004, 26, 103–112. [Google Scholar] [CrossRef]
  21. Lu, M.K.; Cheng, J.J.; Lin, C.Y.; Chang, C.C. Purification, structural elucidation, and anti-inflammatory effect of a water-soluble 1,6-branched 1,3-α-d-galactan from cultured mycelia of Poria cocos. Food Chem. 2010, 118, 349–356. [Google Scholar] [CrossRef]
  22. Liu, X.; Yu, X.; Xu, X.; Zhang, X.; Zhang, X. The protective effects of Poria cocos-derived polysaccharide CMP33 against IBD in mice and its molecular mechanism. Food Funct. 2018, 9, 5936–5949. [Google Scholar] [CrossRef]
  23. Park, W.-H.; Lee, S.-K.; Oh, H.-K.; Bae, J.-Y.; Kim, C.-H. Tumor initiation inhibition through inhibition COX-1 activity of a traditional Korean herbal prescription, Geiji-Bokryung-Hwan, in human hepatocarcinoma cells. Immunopharmacol. Immunotoxicol. 2005, 27, 473–483. [Google Scholar] [CrossRef]
  24. Nie, A.; Chao, Y.; Zhang, X.; Jia, W.; Zhou, Z.; Zhu, C. Phytochemistry and Pharmacological Activities of Wolfiporia cocos (F.A. Wolf) Ryvarden & Gilb. Front. Pharmacol. 2020, 11, 1–17. [Google Scholar] [CrossRef]
  25. Lee, S.R.; Lee, S.; Moon, E.; Park, H.J.; Park, H.B.; Kim, K.H. Bioactivity-guided isolation of anti-inflammatory triterpenoids from the sclerotia of Poria cocos using LPS-stimulated Raw264.7 cells. Bioorg. Chem. 2017, 70, 94–99. [Google Scholar] [CrossRef] [PubMed]
  26. Hopkins, A.L. Network pharmacology. Nat. Biotechnol. 2007, 25, 1110–1111. [Google Scholar] [CrossRef]
  27. Chen, W.; Li, J.; Li, C.; Fan, H.N.; Zhang, J.; Zhu, J.S. Network pharmacology-based identification of the antitumor effects of taraxasterol in gastric cancer. Int. J. Immunopathol. Pharmacol. 2020, 34, 2058738420933107. [Google Scholar] [CrossRef]
  28. Jiang, L.; Lu, J.; Qin, Y.; Jiang, W.; Wang, Y. Antitumor effect of guava leaves on lung cancer: A network pharmacology study. Arab. J. Chem. 2020, 13, 7773–7797. [Google Scholar] [CrossRef]
  29. Chen, Y.; Kern, T.S.; Kiser, P.D.; Palczewski, K. Eyes on systems pharmacology. Pharmacol. Res. 2016, 114, 39–41. [Google Scholar] [CrossRef] [Green Version]
  30. Guo, W.; Huang, J.; Wang, N.; Tan, H.Y.; Cheung, F.; Chen, F.; Feng, Y. Integrating network pharmacology and pharmacological evaluation for deciphering the action mechanism of herbal formula Zuojin Pill in suppressing hepatocellular carcinoma. Front. Pharmacol. 2019, 10, 1–21. [Google Scholar] [CrossRef]
  31. Oh, K.; Adnan, M.; Cho, D. Uncovering Mechanisms of Zanthoxylum piperitum Fruits for the Alleviation of Rheumatoid Arthritis Based on Network Pharmacology. Biology 2021, 10, 703. [Google Scholar] [CrossRef]
  32. Xu, X.; Zhang, W.; Huang, C.; Li, Y.; Yu, H.; Wang, Y.; Duan, J.; Ling, Y. A Novel Chemometric Method for the Prediction of Human Oral Bioavailability. Int. J. Mol. Sci. 2012, 13, 6964–6982. [Google Scholar] [CrossRef]
  33. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 2017, 7, 42717. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Piñero, J.; Ramírez-Anguita, J.M.; Saüch-Pitarch, J.; Ronzano, F.; Centeno, E.; Sanz, F.; Furlong, L.I. The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. 2020, 48, D845–D855. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Rappaport, N.; Twik, M.; Plaschkes, I.; Nudel, R.; Iny Stein, T.; Levitt, J.; Gershoni, M.; Morrey, C.P.; Safran, M.; Lancet, D. MalaCards: An amalgamated human disease compendium with diverse clinical and genetic annotation and structured search. Nucleic Acids Res. 2017, 45, D877–D887. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Hamosh, A.; Scott, A.F.; Amberger, J.S.; Bocchini, C.A.; McKusick, V.A. Online Mendelian Inheritance in Man (OMIM), a knowledgebase of human genes and genetic disorders. Nucleic Acids Res. 2005, 33, D514. [Google Scholar] [CrossRef] [PubMed]
  37. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software Environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  38. Lu, C.; Hu, X.; Wang, G.; Leach, L.J.; Yang, S.; Kearsey, M.J.; Luo, Z.W. Why do essential proteins tend to be clustered in the yeast interactome network? Mol. Biosyst. 2010, 6, 871–877. [Google Scholar] [CrossRef]
  39. Chin, C.; Chen, S.; Wu, H.; Ho, C.; Ko, M.; Lin, C. cytoHubba: Identifying hub objects and sub- networks from complex interactome. BMC Syst. Biol. 2014, 8, S11. [Google Scholar] [CrossRef] [Green Version]
  40. Huang, Z.; Shi, X.; Li, X.; Zhang, L.; Wu, P.; Mao, J.; Xing, R.; Zhang, N.; Wang, P. Network Pharmacology Approach to Uncover the Mechanism Governing the Effect of Simiao Powder on Knee Osteoarthritis. Biomed Res. Int. 2020, 2020, 6971503. [Google Scholar] [CrossRef]
  41. Xiong, Y.; Yang, Y.; Xiong, W.; Yao, Y.; Wu, H.; Zhang, M. Network pharmacology-based research on the active component and mechanism of the antihepatoma effect of Rubia cordifolia L. J. Cell. Biochem. 2019, 120, 12461–12472. [Google Scholar] [CrossRef] [PubMed]
  42. Rudra, S.; Tahamina, A.; Emon, N.U.; Adnan, M.; Shakil, M.; Chowdhury, M.H.U.; Barlow, J.W.; Alwahibi, M.S.; Soliman Elshikh, M.; Faruque, M.O.; et al. Evaluation of Various Solvent Extracts of Tetrastigma leucostaphylum (Dennst.) Alston Leaves, a Bangladeshi Traditional Medicine Used for the Treatment of Diarrhea. Molecules 2020, 25, 4994. [Google Scholar] [CrossRef] [PubMed]
  43. Adnan, M.; Nazim Uddin Chy, M.; Mostafa Kamal, A.T.M.; Azad, M.O.K.; Paul, A.; Uddin, S.B.; Barlow, J.W.; Faruque, M.O.; Park, C.H.; Cho, D.H. Investigation of the biological activities and characterization of bioactive constituents of ophiorrhiza rugosa var. prostrata (D.Don) & Mondal leaves through in vivo, in vitro, and in silico approaches. Molecules 2019, 24, 1367. [Google Scholar] [CrossRef] [Green Version]
  44. Shovo, M.A.R.B.; Tona, M.R.; Mouah, J.; Islam, F.; Chowdhury, M.H.U.; Das, T.; Paul, A.; Ağagündüz, D.; Rahman, M.M.; Emran, T.B.; et al. Computational and Pharmacological Studies on the Antioxidant, Thrombolytic, Anti-Inflammatory, and Analgesic Activity of Molineria capitulata. Curr. Issues Mol. Biol. 2021, 43, 434–456. [Google Scholar] [CrossRef]
  45. Bochevarov, A.D.; Harder, E.; Hughes, T.F.; Greenwood, J.R.; Braden, D.A.; Philipp, D.M.; Rinaldo, D.; Halls, M.D.; Zhang, J.; Friesner, R.A. Jaguar: A high-performance quantum chemistry software program with strengths in life and materials sciences. Int. J. Quantum Chem. 2013, 113, 2110–2142. [Google Scholar] [CrossRef]
  46. Ko, S.G.; Yin, C.S.; Du, B.; Kim, K. Herbal medicines for inflammatory diseases. Mediators Inflamm. 2014, 2014, 982635. [Google Scholar] [CrossRef]
  47. Qi, Q.; Li, R.; Li, H.Y.; Cao, Y.B.; Bai, M.; Fan, X.J.; Wang, S.Y.; Zhang, B.; Li, S. Identification of the anti-tumor activity and mechanisms of nuciferine through a network pharmacology approach. Acta Pharmacol. Sin. 2016, 37, 963–972. [Google Scholar] [CrossRef] [Green Version]
  48. Ramakrishnan, S.; Anand, V.; Roy, S. Vascular Endothelial growth factor signaling in hypoxia and Inflammation. J. Neuroimmune Pharmacol. 2014, 9, 142. [Google Scholar] [CrossRef] [Green Version]
  49. Yoo, S.-A.; Bae, D.-G.; Ryoo, J.-W.; Kim, H.-R.; Park, G.-S.; Cho, C.-S.; Chae, C.-B.; Kim, W.-U. Arginine-Rich Anti-Vascular Endothelial Growth Factor (Anti-VEGF) Hexapeptide Inhibits Collagen-Induced Arthritis and VEGF-Stimulated Productions of TNF-α and IL-6 by Human Monocytes. J. Immunol. 2005, 174, 5846–5855. [Google Scholar] [CrossRef] [Green Version]
  50. Arman, M.; Payne, H.; Ponomaryov, T.; Brill, A. Role of Platelets in Inflammation. In The Non-Thrombotic Role of Platelets in Health and Disease; IntechOpen: London, UK, 2015. [Google Scholar] [CrossRef]
  51. Kopecki, Z.; Ludwig, R.J.; Cowin, A.J. Cytoskeletal Regulation of Inflammation and Its Impact on Skin Blistering Disease Epidermolysis Bullosa Acquisita. Int. J. Mol. Sci. 2016, 17, 1116. [Google Scholar] [CrossRef] [Green Version]
  52. Catanzaro, J.M.; Sheshadri, N.; Pan, J.-A.; Sun, Y.; Shi, C.; Li, J.; Powers, R.S.; Crawford, H.C.; Zong, W.-X. Oncogenic Ras induces inflammatory cytokine production by upregulating the squamous cell carcinoma antigens SerpinB3/B4. Nat. Commun. 2014, 5, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Multhoff, G.; Molls, M.; Radons, J. Chronic Inflammation in Cancer Development. Front. Immunol. 2012, 2, 98. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Melo, M.B.; Jensen, K.D.C.; Saeij, J.P.J. Toxoplasma gondii effectors are master regulators of the inflammatory response. Trends Parasitol. 2011, 27, 487–495. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Nidai Ozes, O.; Mayo, L.D.; Gustin, J.A.; Pfeffer, S.R.; Pfeffer, L.M.; Donner, D.B. NF-κB activation by tumour necrosis factor requires the Akt serine–threonine kinase. Nature 1999, 401, 82–85. [Google Scholar] [CrossRef]
  56. Taylor, C.T. Interdependent roles for hypoxia inducible factor and nuclear factor-κB in hypoxic inflammation. J. Physiol. 2008, 586, 4055–4059. [Google Scholar] [CrossRef]
  57. Cai, Y.; Sukhova, G.K.; Wong, H.K.; Xu, A.; Tergaonkar, V.; Vanhoutte, P.M.; Tang, E.H.C. Rap1 induces cytokine production in pro-inflammatory macrophages through NFκB signaling and is highly expressed in human atherosclerotic lesions. Cell Cycle 2015, 14, 3580. [Google Scholar] [CrossRef] [Green Version]
  58. Dalal, P.J.; Muller, W.A.; Sullivan, D.P. Endothelial Cell Calcium Signaling during Barrier Function and Inflammation. Am. J. Pathol. 2020, 190, 535–542. [Google Scholar] [CrossRef] [Green Version]
  59. Chen, Q.; He, G.; Zhang, W.; Xu, T.; Qi, H.; Li, J.; Zhang, Y.; Gao, M.Q. Stromal fibroblasts derived from mammary gland of bovine with mastitis display inflammation-specific changes. Sci. Rep. 2016, 6, 1–13. [Google Scholar] [CrossRef] [Green Version]
  60. Konisti, S.; Kiriakidis, S.; Paleolog, E.M. Hypoxia-a key regulator of angiogenesis and inflammation in rheumatoid arthritis. Nat. Rev. Rheumatol. 2012, 8, 153–162. [Google Scholar] [CrossRef]
  61. Wilson, K.T.; Fu, S.; Ramanujam, K.S.; Meltzer, S.J. Increased Expression of Inducible Nitric Oxide Synthase and Cyclooxygenase-2 in Barrett’s Esophagus and Associated Adenocarcinomas. Cancer Res. 1998, 58, 2929–2934. [Google Scholar]
  62. Ying, L.; Hofseth, L.J. An Emerging Role for Endothelial Nitric Oxide Synthase in Chronic Inflammation and Cancer. Cancer Res. 2007, 67, 1407–1410. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Yoo, S.-A.; Yoon, H.-J.; Kim, H.-S.; Chae, C.-B.; De Falco, S.; Cho, C.-S.; Kim, W.-U. Role of placenta growth factor and its receptor flt-1 in rheumatoid inflammation: A link between angiogenesis and inflammation. Arthritis Rheum. 2009, 60, 345–354. [Google Scholar] [CrossRef] [PubMed]
  64. Olson, N.; Vliet, A. van der Interactions between Nitric Oxide and Hypoxia-Inducible Factor Signaling Pathways in Inflammatory Disease. Nitric Oxide 2011, 25, 125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Khan, A.M.; Shawon, J.; Halim, M.A. Multiple receptor conformers based molecular docking study of fluorine enhanced ethionamide with mycobacterium enoyl ACP reductase (InhA). J. Mol. Graph. Model. 2017, 77, 386–398. [Google Scholar] [CrossRef] [PubMed]
  66. Shahinozzaman, M.; Taira, N.; Ishii, T.; Halim, M.; Hossain, M.; Tawata, S. Anti-Inflammatory, Anti-Diabetic, and Anti-Alzheimer’s Effects of Prenylated Flavonoids from Okinawa Propolis: An Investigation by Experimental and Computational Studies. Molecules 2018, 23, 2479. [Google Scholar] [CrossRef] [Green Version]
  67. Pearson, R.G. Absolute electronegativity and hardness correlated with molecular orbital theory. Proc. Natl. Acad. Sci. USA 1986, 83, 8440–8441. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Lucido, M.J.; Orlando, B.J.; Vecchio, A.J.; Malkowski, M.G. Crystal Structure of Aspirin-Acetylated Human Cyclooxygenase-2: Insight into the Formation of Products with Reversed Stereochemistry. Biochemistry 2016, 55, 1226–1238. [Google Scholar] [CrossRef] [Green Version]
Life 13 00893 sch001 550
Scheme 1. Schematic representation of network pharmacology approach for W. extensa.
Scheme 1. Schematic representation of network pharmacology approach for W. extensa.
Life 13 00893 sch001
Life 13 00893 g001 550
Figure 1. GC-MS peak of the methanolic extract of W. extensa and an indication of the key bioactive.
Figure 1. GC-MS peak of the methanolic extract of W. extensa and an indication of the key bioactive.
Life 13 00893 g001
Life 13 00893 g002 550
Figure 2. Intersecting targets. (A) Bioactives related targets from SEA and STP; (B) Inflammation targets from DisGeNet, Malacards and OMIM, and common bioactives-related targets.
Figure 2. Intersecting targets. (A) Bioactives related targets from SEA and STP; (B) Inflammation targets from DisGeNet, Malacards and OMIM, and common bioactives-related targets.
Life 13 00893 g002
Life 13 00893 g003 550
Figure 3. Bioactives—Target network interaction between WE bioactives and common inflammation targets.
Figure 3. Bioactives—Target network interaction between WE bioactives and common inflammation targets.
Life 13 00893 g003
Life 13 00893 g004 550
Figure 4. Bubble plot of significantly enriched pathways of common targets in KEGG evaluation.
Figure 4. Bubble plot of significantly enriched pathways of common targets in KEGG evaluation.
Life 13 00893 g004
Life 13 00893 g005 550
Figure 5. Pathway-Gene network of common targets.
Figure 5. Pathway-Gene network of common targets.
Life 13 00893 g005
Life 13 00893 g006 550
Figure 6. Illustration of HIF-1 signaling pathway. The red node represents the enriched targets in the present study.
Figure 6. Illustration of HIF-1 signaling pathway. The red node represents the enriched targets in the present study.
Life 13 00893 g006
Life 13 00893 g007 550
Figure 7. Binding relationship between targets enriched in the HIF-1 signaling pathway and the key bioactive “N-(3-chlorophenyl)naphthylcarboxamide” from MEWE. (A) 1M9J (B) 1NSI (C) 3HNG (D) 3UL7 and (E) 5WB7.
Figure 7. Binding relationship between targets enriched in the HIF-1 signaling pathway and the key bioactive “N-(3-chlorophenyl)naphthylcarboxamide” from MEWE. (A) 1M9J (B) 1NSI (C) 3HNG (D) 3UL7 and (E) 5WB7.
Life 13 00893 g007
Life 13 00893 g008 550
Figure 8. Frontier Molecular Orbital plots and their energy for key compounds from MEWE and standard medicines.
Figure 8. Frontier Molecular Orbital plots and their energy for key compounds from MEWE and standard medicines.
Life 13 00893 g008
Life 13 00893 g009 550
Figure 9. Mechanistic mode of action of W. extensa against inflammation.
Figure 9. Mechanistic mode of action of W. extensa against inflammation.
Life 13 00893 g009
Table
Table 1. A list of 27 GC-MS-detected chemical components in methanolic extract of W. extensa.
Table 1. A list of 27 GC-MS-detected chemical components in methanolic extract of W. extensa.
S.N.R.T. (min)Area (%)PubChem
CID		Chemical FormulaBioactives ClassBioactives Name
14.380.4354544338C5H6N2O2Organooxygen compoundsN-Cyano-3-oxobutanamide
24.430.8580975C7H10O2Organooxygen compounds1,3-Cyclopentanedione, 2,4-dimethyl- or 2,4-Dimethyl-1,3-cyclopentanedione
34.831.76119838C6H8O4Pyrans2,3-Dihydro-3,5-dihydroxy-6-methyl-4h-pyran-4-one or 3-Hydroxy-2,3-dihydromaltol
45.454.5812991C4H6OOrganooxygen compounds2-Butyn-1-ol or, 2-Butynol
55.680.74643131C4H5ClOAcyl halides(2E)-2-Butenoyl chloride or 2-Butenoyl chloride
65.750.82123411C4H3NSIsothiocyanatesThiocyanic acid, 2-propynyl ester or Propargyl isothiocyanate
75.990.9976665C6H10N2OAzolines3H-Pyrazol-3-one, 2,4-dihydro-4,4,5-trimethyl- or 3,4-Trimethyl-5-pyrazolone
86.080.685362763C8H14OOrganooxygen compounds2-Heptenal, 2-methyl- or 2-methyl-2-heptenal
96.21.968102C6H15NOrganooxygen compounds1-Hexanamine or Hexylamine
106.492.41538272C8H15BO3Boronic acid derivativesLactic acid, 2-methyl-, monoanhydride with 1-butaneboronic acid, cyclic ester or alpha-Hydroxyisobutyric acid cyclic butaneboronate
116.872.8811850C6H14O6Organooxygen compoundsGalactitol
127.087.145951C3H7NO3Carboxylic acids and derivativesL-Serine or Serine
137.310.85366263C19H38OOrganooxygen compoundsEther, methyl 1-octadecenyl or 1-Methoxy-1-octadecene
147.350.4111005C14H28O2Fatty acylsTetradecanoic acid or Myristic acid
157.6424.35780C6H14O6Organooxygen compoundsSorbitol
168.740.41554151C17H34O2Fatty acylsPentadecanoic acid, 13-methyl-, methyl ester or 13-Methylpentadecanoic acid methyl ester
178.91.21985C16H32O2Fatty acylsHexadecanoic acid or Palmitic acid
189.391.275284421C19H34O2Fatty acyls9,12-Octadecadienoic acid (Z,Z)-, methyl ester or Methyl linoleate
199.571.335280450C18H32O2Fatty acyls9,12-Octadecadienoic acid (Z,Z)- or, Linoleic acid
2010.21.12610075C10H6ClNO2NaphthalenesNaphthalene, 6-chloro-1-nitro- or 6-Chloro-1-nitronaphthalene
2110.31.07581332C15H23Cl2N3Benzene and substituted derivativesHexahydropyrazin-1-propylamine, 4-[2-[3,4-dichlorophenyl]ethyl]-
22113.295789C10H14N2O5Pyrimidine nucleosidesThymidine
2312.94.9610850C27H50O6GlycerolipidsGlycerol tricaprylate or Tricaprylin
2414.36.43159892C4H5N3O2Organonitrogen compounds2-Hydroxy-4-hydroxylaminopyrimidine or, N4-hydroxycytosine
2515.65.756432563C28H44OSteroids and steroid derivativesErgosta-5,7,22-trien-3-ol, (3.beta.,22E)- or (3.beta.)-Ergosta-5,7,22- trien-3-ol
2616.34.24532200C17H12ClNONaphthalenes1-Naphthalenecarboxamide, N-(3-chlorophenyl)- or N-(3-chlorophenyl)naphthylcarboxamide
2716.61.8312796C28H48OSteroids and steroid derivatives.alpha.-Ergostenol
Notes: S.N. = Serial number; R.T. = Retention time.
Table
Table 2. Classified orders of PPI network of 42 common targets based on different algorithms.
Table 2. Classified orders of PPI network of 42 common targets based on different algorithms.
MCCMNCDegree
RankGenesScoreRankGenesScoreRankGenesDegree
1TLR412981TLR4181TLR418
2EGFR12321EGFR182EGFR18
3NOS311803NOS3163NOS316
4PTPRC11583PTPRC164PTPRC16
5MAPK1410855F2155MAPK1415
6PPARG9696MAPK14146F215
7FGF29207TRPV1137TRPV113
8JAK22178FGF2128FGF212
9NOS21489ALOX5119ALOX511
10CNR114610CNR11010PPARG11
11SIRT114210PPARG1011CNR110
12CNR213410PTGER41012PTGER410
13ADORA113210SIRT11013SIRT110
14PTGER313014PTGER3914PTGER39
15PTGDR212214PTGES915PTGES9
16NPY5R12116CNR2816JAK29
17LCK9616NOS2817CNR28
18PTGER48416JAK2818PLA2G4A8
19TRPV18216PLA2G4A819LTB4R8
20F27220ADORA1720F2R8
21PDGFRB5420HRH1721NOS28
21ALOX55420LCK722ADORA17
23PTGES3820TRPA1723PTGDR27
24FLT13620F2R724LCK7
25F2R3325PDGFRB625TRPA17
26PLA2G4A3225FLT1626PPARA7
26HRH13225PLAU627HRH17
26CD383225CD38628NPY5R6
29LTB4R3029PTGDR2529CD386
30PTGFR2629NPY5R530PDGFRB6
31PLAU2429PARP1531FLT16
32TRPA11629PLAT532PLAU6
33PLAT1429PTGFR533NR1H45
34PLA2G2A1229PPARA534PARP15
34PPARA1235NR1H4435PLAT5
36PARP11035LTB4R436PTGFR5
37NR1H4735PLA2G2A437CYP2C194
37CYP2C19738CYP2C19338PLA2G2A4
39KCNK3638KCNK3339FABP14
39KCNK9638KCNK9340KCNK33
41FABP1538FABP1341KCNK93
42GBA242GBA242GBA2
Table
Table 3. Docking score of highly enriched pathway targets of inflammation with key bioactives, standard medicines and theirs co-crystallized ligands.
Table 3. Docking score of highly enriched pathway targets of inflammation with key bioactives, standard medicines and theirs co-crystallized ligands.
GenePDB IDBioactivesDocking Score
Kcal/mol
TLR43UL7N-(3-chlorophenyl)naphthylcarboxamide−4.027
Aspirin *−4.178
Indomethacin *−4.168
2-acetamido−2-deoxy-beta-d-glucopyranose−4.797
EGFR5WB7N-(3-chlorophenyl)naphthylcarboxamide−4.123
Aspirin *−4.182
Indomethacin *−4.707
2-acetamido-2-deoxy-beta-d-glucopyranose−4.205
FLT13HNGN-(3-chlorophenyl)naphthylcarboxamide−8.88
Aspirin * −6.794
Indomethacin *−5.538
N-(4-chlorophenyl)-2-[(pyridin-4-ylmethyl)amino]benzamide−11.044
NOS31M9JN-(3-chlorophenyl)naphthylcarboxamide−7.236
Aspirin *−5.262
Indomethacin *−7.451
Chlorzoxazone−6.118
NOS21NSIN-(3-chlorophenyl)naphthylcarboxamide−8.85
Aspirin *−6.467
Indomethacin *−7.491
Protoporphyrin Ix Containing Fe−14.755
Notes: * denotes standard medicine.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jin, J.; Chowdhury, M.H.U.; Hafizur Rahman, M.; Choi, K.-Y.; Adnan, M. Bioactive Compounds and Signaling Pathways of Wolfiporia extensa in Suppressing Inflammatory Response by Network Pharmacology. Life 2023, 13, 893. https://doi.org/10.3390/life13040893

AMA Style

Jin J, Chowdhury MHU, Hafizur Rahman M, Choi K-Y, Adnan M. Bioactive Compounds and Signaling Pathways of Wolfiporia extensa in Suppressing Inflammatory Response by Network Pharmacology. Life. 2023; 13(4):893. https://doi.org/10.3390/life13040893

Chicago/Turabian Style

Jin, Juri, Md. Helal Uddin Chowdhury, Md. Hafizur Rahman, Ki-Young Choi, and Md. Adnan. 2023. "Bioactive Compounds and Signaling Pathways of Wolfiporia extensa in Suppressing Inflammatory Response by Network Pharmacology" Life 13, no. 4: 893. https://doi.org/10.3390/life13040893

APA Style

Jin, J., Chowdhury, M. H. U., Hafizur Rahman, M., Choi, K. -Y., & Adnan, M. (2023). Bioactive Compounds and Signaling Pathways of Wolfiporia extensa in Suppressing Inflammatory Response by Network Pharmacology. Life, 13(4), 893. https://doi.org/10.3390/life13040893

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop