Next Article in Journal
Objective Perfusion Assessment in Gracilis Muscle Interposition—A Novel Software-Based Approach to Indocyanine Green Derived Near-Infrared Fluorescence in Reconstructive Surgery
Next Article in Special Issue
The Case for Clinical Trials with Novel GABAergic Drugs in Diabetes Mellitus and Obesity
Previous Article in Journal
Combining Direct PCR Technology and Capillary Electrophoresis for an Easy-to-Operate and Highly Sensitive Infectious Disease Detection System for Shrimp
Previous Article in Special Issue
Pharmacological Effect of GABA Analogues on GABA-ϱ2 Receptors and Their Subtype Selectivity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Network Pharmacology Study to Reveal the Potentiality of a Methanol Extract of Caesalpinia sappan L. Wood against Type-2 Diabetes Mellitus

1
Department of Bio-Health Convergence, College of Biomedical Science, Kangwon National University, Chuncheon 24341, Korea
2
Ethnobotany and Pharmacognosy Lab, Department of Botany, University of Chittagong, Chattogram 4331, Bangladesh
3
Department of Microbiology, University of Chittagong, Chattogram 4331, Bangladesh
4
Department of Pharmacy, International Islamic University Chittagong, Chittagong 4318, Bangladesh
*
Author to whom correspondence should be addressed.
Life 2022, 12(2), 277; https://doi.org/10.3390/life12020277
Submission received: 2 January 2022 / Revised: 8 February 2022 / Accepted: 10 February 2022 / Published: 13 February 2022
(This article belongs to the Special Issue New Trends in Pharmaceutical Science)

Abstract

:
Caesalpinia sappan L. (CS) is widely used to treat diabetic complications in south-east Asia, specifically in traditional Chinese medicine. This study intends to explain the molecular mechanism of how chemical constituents of CS interrelate with different signaling pathways and receptors involved in T2DM. GC-MS was employed to identify the chemical compounds from the methanol extract of CS wood (MECSW). Lipinski’s rule of five was applied, and 33 bioactive constituents have been screened from the CS extract. After that, 124 common targets and 26 compounds associated with T2DM were identified by mining several public databases. Protein–protein interactions and compound-target network were constructed using the STRING database and Cytoscape tool. Protein–protein interactions were identified in 121 interconnected nodes active in T2DM and peroxisome proliferator-activated receptor gamma (PPARG) as key target receptors. Furthermore, pathway compound target (PCT) analysis using the merger algorithm plugin of Cytoscape revealed 121 nodes from common T2DM targets, 33 nodes from MECSW compounds and 9 nodes of the KEGG pathway. Moreover, network topology analysis determined “Fisetin tetramethyl ether” as the key chemical compound. The DAVID online tool determined seven signaling receptors, among which PPARG was found most significant in T2DM progression. Gene ontology and KEGG pathway analysis implied the involvement of nine pathways, and the peroxisome proliferator-activated receptor (PPAR) pathway was selected as the hub signaling pathway. Finally, molecular docking and quantum chemistry analysis confirmed the strong binding affinity and reactive chemical nature of fisetin tetramethyl ether with target receptors exceeding that of the conventional drug (metformin), PPARs agonist (rosiglitazone) and co-crystallized ligands, indicating that fisetin could be a potential drug of choice in T2DM management. This study depicts the interrelationship of the bioactive compounds of MECSW with the T2DM-associated signaling pathways and target receptors. It also proposes a more pharmaceutically effective substance, fisetin tetramethyl ether, over the standard drug that activates PPARG protein in the PPAR signaling pathway of T2DM.

1. Introduction

Type-2 diabetes mellitus (T2DM) is a global epidemic attributed to the dysregulation of carbohydrate, lipid, and protein metabolism resulting from impaired insulin secretion, insulin resistance, or a combination of both. According to the International Diabetes Federation (IDF) report in 2021, one in ten adults (20 to 79 years) has diabetes mellitus. The IDF also estimated this rate to rise to 643 million by 2030 and 783 million by 2045. Such a significant increase will come from the economic transition from low-income to middle-income levels [1].
T2DM is a complex chronic disease involving genetic predisposition, environmental and behavioral risk factors [2,3,4]. Upon the dysfunction of feedback loops between insulin action and insulin secretion, insulin-sensitive tissues, such as liver, muscle, and the adipose tissue, are affected, resulting in abnormal insulin secretion by pancreatic islet β cells and abnormal glucose levels in the blood; this impaired insulin secretion comes forth as rooted sake of T2DM [5]. Insulin resistance causes type-2 diabetes, which mediates increased glucose synthesis in the liver and decreases glucose intake in muscle and adipose tissue at a certain insulin level. Moreover, chronic inflammation leads to impaired control of blood glucose levels, resulting in vascular complications [6]. Likewise, β-cell dysfunction causes reduced insulin release resulting in blood glucose homeostasis failure [7] and the accumulation of reactive oxygen species (ROS) in β cells, also accountable for insulin deficiency through immutable damage in mitochondria [8]. Individuals who suffer from this noteworthy metabolic condition are at a greater risk of mortality from myocardial disease, renal disease, virus-related sickness, and respiratory diseases, among other causes [9,10,11].
Pharmacologically significant agents, such as biguanides, sulfonylureas, meglitinides, thiazolidinediones, alpha-glucosidase inhibitors, incretin-based therapies, dipeptidyl-peptidase IV inhibitors, insulin analogs, and bromocriptine, have been reported as potential for T2DM patients in clinical scenarios [9,12,13,14,15]. Several oral anti-diabetic drugs, including metformin, glimepiride, repaglinide, pioglitazone, sitagliptin, and acarbose, have been used frequently as therapeutic agents. However, these drugs often create side effects, for instance, lactic acidosis, hypoglycemia, anorexia, nausea, dyspeptic episodes, impaired renal function, and other gastrointestinal issues [16,17,18]. While synthetic drugs pose some adverse effects in the patient, a wide range of plant-derived compounds are showing promising outcomes in managing T2DM [19]. That is why plant-derived metabolites have been way ahead of synthetic ones in the quest for a safer therapeutic agent of T2DM. Since different medicinal plants are reported and widely used to manage diabetes as per traditional practice in the south Asian region, this study aims to draw some insight into the selected plant’s anti-diabetic activity, Caesalpinia sappan, which has been reported as having anti-diabetic activity in several recent studies [20,21,22,23,24].
Caesalpinia sappan L. is also known as Biancaea sappan (L.) Tod and belongs to the Caesalpiniaceae family. It is commonly recognized as sappan wood and Indian redwood. It has traditionally been utilized in Ayurvedic diabetes therapy [25]. Water extracts of C. Sappan wood (CSW) have been widely used to treat diabetic complications in south-east Asia, specifically in Chinese traditional medicine [26]. CSW exhibited inhibitory activities towards LPS-induced NO production in macrophage, cytotoxicity against human pancreatic cancer cells, breast cancer cell, and colon cancer cell lines [27,28,29,30]. It also showed significant anti-inflammatory activity in many studies [31,32,33]. Moreover, a methanolic extract of C. sappan showed significant anti-diabetic activity in alloxan-induced diabetes mellitus in rats and alpha-glucosidase inhibitory activity [34,35]. Brazilin, extracted from CSW, increased the synthesis of fructose-2,6-bisphosphate [36]. Fructose-2,6-bisphosphate synthesized from CS stimulated the process of glycolysis by activating phosphofructokinase-1 resulted in reduced blood glucose levels [37]. Furthermore, some important compounds—including brazilin, sappanol, and episappanol [38]; dibenz[b,d]oxocins protosappanin B, C, and isoprotosappanin B [39]; and homoisoflavonoids 3′-deoxysappanol, 3-deoxysappanone B, and 4-O-methylsappanol [40]—were extracted from CSW by using high-performance liquid chromatography. Saponin, another compound from CSW, exposed potential anti-diabetic activity by inhibiting glucose transport and stimulating insulin secretion in pancreatic β cells [41,42]. CSW also contains tannins capable of lowering blood glucose by inhibiting α-amylase and α-glucosidase, leading to carbohydrate digestion and glucose absorption during high blood glucose levels after meals [43]. Hence, due to the promising anti-diabetic activity of C. sappan wood (CSW), we intend to predict its interaction mechanism with metabolic pathways and genes that express proteins involved in the progression of T2DM.
A more efficient and safer treatment for T2DM demands a proper explanation of its molecular etiology. Network pharmacology is a well-defined method to investigate the interactions between target compounds, genes, and proteins associated with the disease [44,45]. It can reveal the molecular mechanism of the compounds from a multifunctional point of view, highlighting the interaction of diverse factors involved in the disease progression and elucidating the synergistic impact of the bioactive compounds in a living cell [43,46]. With the advent of modern bioinformatics and systems biology, poly-pharmacological methods contribute suggestively to network-based drug discovery as a cost-effective drug developing tool [47]. Network pharmacology has proven beneficial in explaining the underlying complex relationship between the pharmacological properties of given compounds and the whole biochemical pathway involved in a particular disease [48,49]. In our study, network pharmacology is employed to evaluate the bioactive compounds from CSW and their mechanism against T2DM. Bioactive compounds from CSW are identified using GC-MS analysis. Following that, overlapping genes related to the selected compounds and T2DM are identified using the public databases. Then genes involved in hub signaling are chosen by analyzing gene set analysis. Finally, the most potent candidates of CSW against T2DM are determined by implementing molecular docking analysis for the selected genes. The detailed process is depicted in Scheme 1.

2. Materials and methods

2.1. Plant Collection and Extraction

Caesalpinia sappan wood (CSW) was collected from Chuncheon local market and was authenticated by Dr. Dong Ha Cho, plant biologist and Professor, Department of Bio-Health Convergence, College of Biomedical Science, Kangwon National University. A voucher number (CRT 215) was stored at Kenaf Corporation in the Department of Bio-Health Convergence, and the collected material can be used only for research purposes. The collected dried wood (500 g) was grinded into a coarse powder using an automated grinder. The fined powder (100 g) was soaked in 500 mL of methanol (Daejung, Siheung City, Korea). The mixture was put in a sealed bottle and placed for continuous shaking and stirring (for 3 days) through an automated shaker machine (repeated 3 times for the highest yield). The mixture was filtered (Whatman qualitative filter paper Grade 1) and evaporated using a vacuum evaporator. The evaporated sample (MECSW) was dried under a hot water bath (IKA, Staufen city, Germany) at 40 °C. The yield was found 3.4 g, which was preserved in a refrigerator (−4 °C) for further GC-MS analysis.

2.2. GC-MS Analysis

In order to perform the GC-MS study on MECSW, we adopted the GC-MS equipment together with an analytical capillary column DB-5MS (30 m × 0.25 µm × 0.25 mm). Our prior research outlined the whole methodology in precise detail [50].

2.3. Filtration of Bioactive Constituents from MECSW

The drug-likeness approach (Lipinski’s rule of five) was used to screen MECSW’s bioactive compounds (identified by GC-MS), which overcome ADME (Absorption, Distribution, metabolism and Excretion) restrictions while securing oral bioavailability score > 0.50. An online program called Swiss ADME [51] was used to explore the drug-likeness properties of identified compounds. In order to do this calculation, we utilized compound’s SMILES from PubChem (https://pubchem.ncbi.nlm.nih.gov/, accessed on 27 July 2021) database.

2.4. Acquisition of Compound and T2DM Associated Targets

By using the Homo sapiens mode, we collected targets linked to the filtered bioactive compounds by putting their SMILES into the SEA (Similarity Ensemble Approach) (http://sea.bkslab.org/, accessed on 2 August 2021) and STP (Swiss Target Prediction) (http://www.swisstargetprediction.ch/, accessed on 3 August 2021) databases, respectively. The DisGeNeT (https://www.disgenet.org/search, accessed on 5 August 2021) [52], Malacards (https://www.malacards.org/, accessed on 6 August 2021) [53], and the OMIM (https://www.ncbi.nlm.nih.gov/omim, accessed on 6 August 2021) [54] databases were used to collect targets interacted with T2DM. On the other hand, VENNY 2.1 (https://bioinfogp.cnb.csic.es/tools/venny/, accessed on 15 August 2021) identified and exhibited the common overlapping targets between MECSW bioactive compounds and T2DM targets.

2.5. Creating a Network Involving Intersecting Targets

Homo sapiens with a confidence level of >0.4 in protein interactions, the intersected targets between the compound-related gene and T2DM target genes were included into the STRING database Version 11.0b (https://string-db.org/, accessed on 24 August 2021) for protein queries. The Cytoscape 3.8.2 software program [55] was used to subsequently classify the entire network employing the CytoHubba module contained in Cytoscape and following the degree algorithm to identify its key essential genes:
Deg(v) = |N(v)|
where, a node’s neighbors are represented by N(v), while each node’s neighbors are represented by v.

2.6. Network Layout for the Pathway Compound Target (PCT)

A graphical depiction of the pathway compound target (PCT) network was created using the preprocessing output of bioactive chemicals from MECSW and frequent T2DM targets that intersected with the MECSW. This network was built using Cytoscape’s merger algorithm plugin (Cytoscape 3.8.2). An analysis of network topology parameters was performed with the help of the network analyzer. Nodes represented bioactive compounds, targets, and pathways, and interactions between these components were shown along the edges. The degree also referred to the prevalence of a node’s interconnected neighbors. The greater the percentage of a node’s that are directly linked to each one, the greater the impact [56].

2.7. Investigation of the Role of GO and KEGG Pathways in Common Intersected Targets

The KEGG pathway interpretation and molecular functional annotation of all intersected targets were undertaken on the DAVID (https://david.ncifcrf.gov/tools.jsp, accessed on 19 September 2021) (Database for Annotation, Visualization and Integrated Discovery) database to determine their role in signaling pathways. For this enrichment study, the identifier was OFFICIAL GENE SYMBOL, and the species was Homo sapiens. This database is exceptionally crucial in network pharmacology because it shows targets implicated in disease underlying molecular mechanisms, and the GO database exhibits the descriptive biological terms of those targets, such as biological process (BP), cellular component (CC), and molecular function (MF) [57]. A p-value <0.05 was used as a cutoff for GO and pathway enrichment. The p-value was corrected using a false discovery rate (FDR) error control technique, and the result was known as the Q value. It was possible to visualize the KEGG pathway bubble plot map visually by using Origin Pro 2021 to examine the routes.

2.8. Formulation and Purification of the Ligand and Receptor Protein

The PubChem chemical library was used to obtain the .sdf files of the preferred ligands revealed via compound-target network research along with conventional pharmaceuticals, such as metformin, and co-crystallized protein ligands. As previously reported, the constituents were indeed prepared for molecular docking investigations incorporating the LigPrep tool in Schrödinger Suite-Maestro v12.5 [58]. PPARG (PDB ID: 3E00), PPARA (PDB ID: 1K7L), PPARD (PDB ID: 5U3Q), FABP3 (PDB ID: 5HZ9), FABP4 (PDB ID: 3P6D), MMP1 (PDB ID: 1SU3) and NR1H3 (PDB ID:1UHL), seven receptor proteins essential in hub signaling pathways, were chosen because their crystal structures were available in the RCSB Protein Data Bank (https://www.rcsb.org/, accessed on 11 October 2021) and UniProt database (https://www.uniprot.org/, accessed on 11 October 2021), respectively. Following our previously reported processes, we set the Schrödinger Suite-Maestro v1. 5 integrated Protein Preparation Wizard tools after the 3D crystal structure was located in the RCSB database [59,60].

2.9. Glide Directed Molecular Docking Assay

In order to hatch receptor grids overactive molecules (co-crystallized ligand site), we used Glide tools integrated in Schrödinger Suite-Maestro version 12.5 software [61]. Their default topological settings used during grid creation were a scaling factor 1.00, the OPLS3 force field, and a cut-off value of 0.25% for three-dimensional protein assemblies. On the macromolecules’ kernel active site residues, a cubic box of precise facets was placed in 14 Å × 14 Å × 14 Å grid points, at a feasible docking site. Later, docking studies used Glide’s standard precision (SP) scoring system, with each ligand’s best rating conformation and binding result recorded separately.

2.10. Quantum Chemistry of Key Ingredients

Utilizing the Jaguar panel of Maestro 12.5 software, the key compound (fisetin tetramethyl ether) and standard drugs’ (metformin) structural variables were absolutely optimized using the Lee–Yang–Parr (B3LYP-D3) correlation functional technique at the 6–31G++ (d,p) level basis set [62]. This optimized configuration also generated enthalpy, Gibbs free energy, highest occupied molecular orbital (HOMO) and lowest occupied molecular orbital (LUMO) border molecular orbital energies using the same level of theory. LUMO energy was deducted from the relevant HOMO energy value to compute each chemical’s HOMO–LUMO gaps (Eg). The following equations estimated the hardness (η) and softness (S) based on border HOMO and LUMO energies.
η = (HOMO − LUMO)/2
S = 1/η

3. Results

3.1. Exploration of MECSW Ingredients Using GC-MS

The data obtained from the gas chromatography-mass spectrometry (GC-MS) study revealed 33 significant bioactive constituents in the MECSW (Figure 1). Table 1 shows the retention time (RT), peak area (percent), chemical formula, and tentative identities of the bioactive compounds. All compounds were categorized (belong to the organic kingdom): namely lipids and lipid-like molecules, benzenoids, organic acids and derivatives, organic oxygen compounds, organic 1,3-dipolar compounds, organoheterocyclic compounds, phenylpropanoids and polyketides, and alkaloids and derivatives.

3.2. Drug Candidates Filtering

The primary bioactive ingredients from MECSW were screened based on the drug-likeness properties of recorded small molecules. Those properties were stated as molecular weight must not be over 500, a number of H-bond donor definitely below 10, contrarily H-bond acceptor must not exceed 5, moriguchi octanol–water partition coefficient value must be below or equal to 4.15, and ‘Abott Bioavailability Score’ should be under 0.1 standard value. Remarkably, all of the identified small molecules (33) occupied the above-mentioned criteria and were classified as significant bioactive substances without infringing more than one of the features mentioned earlier (Table 2).

3.3. Common Intersected Targets of Compounds within SEA and STP Database

The screened bioactive compounds were subjected to acquiring compound pertaining targets from public databases. In parallel, each component’s SMILES code was acquired from the PubChem chemical library and entered into the SEA and STP database queries. The removal of duplication targets conveyed the presence of 844 compounds linked to targets from the STP database and 489 targets from the SEA database (Table S1). The result of the Venn diagram analysis disclosed the presence of 21.7% (238) prevalent targets between those 2 databases (Figure 2A).

3.4. Potential Overlapping Targets between T2DM Targets and Compound Linked 238 Intersecting Targets

T2DM targets were gleaned from accessing three disease-related public databases, including DisGeNeT, OMIM and Malacards, which resulted in the procurement of 3081 disease-related targets (Table S2). After that, the culled targets were examined with 238 compound-related targets to determine the T2DM targets that were directly connected with the various MECSW compounds (Table S3). Consequently, 124 common targets (Table S4) directly related to T2DM and MECSW compounds were outlined (Figure 2B), in which 26 molecules were closely connected to typical T2DM targets, while the remaining seven compounds, including (3E)-5-Hydroxy-2-methyl-3-hexenoic acid, Disalicylalpropylenediimine, 2-[(4E)-4-Hexenyl]-6-nitrocyclohexanone dimethylhydrazone, 2-(4-methylphenyl) Indolizine, Aspidodispermine, O-methyl-, Nadolol di-methylboronic acid, and 1,3-Dimethyl-5,6-dicarbethoxy-5,6,7,8-tetrahydro-6,7-diazalumazine, were not affiliated with any targets in neither database (SEA and STP). Notably, those compounds were associated with the genes involved in T2DM progression. These compounds were resulting from the intersection of T2DM related genes and target genes related to MECSW compounds

3.5. PPI Network Analysis of 124 Common Targets

To unearth the potential mechanistic insight of MECSW to treat T2DM, we introduced those 124 common targets into the STRING database to build up a network within them. Meanwhile, the STRING algorithm expressed those 121 nodes connected by forming 596 edges (Figure 3). The average number of neighbors was six, whilst the network diameter was three. However, three targets, namely NEK6, ST6GAL1 and PDE4D, did not interact with any other nodes and were thus omitted from this analysis. To scrutinize the essential key target in the network of T2DM, we further visualized this network in Cytoscape, where we used a degree value algorithm by employing cytoHubba apps. The number of edges connecting to the corresponding target nodes was defined as the number of degrees for each target. Notably, a higher degree value pinpoints the best target in the network. With such conformity, PPARG was designated as a key target (36-degree value) in the network for T2DM progression (Table S5).

3.6. Pathway Compound Target (PCT) Network Analysis

To ascertain the most significant key components among the detected ingredients from MECSW, we constructed a pathway compound target network utilizing Cytoscape. This network displayed how pathway, compound and T2DM targets correlated among them. The topological parameter analysis of PCT network via network analyzer apps integrated into Cytoscape revealed that this network consists of a total of 166 (121 nodes from common T2DM targets, 33 nodes from MECSW compounds and 9 nodes of KEGG pathway) nodes along with interacting within them through 319 edges (Figure 4), whereas we classified the constituents in the network based on their degree value, referred to as their number of connections to relevant targets. Based on this, the fisetin tetramethyl ether component was picked as a key hub substance (27 degree) in the network that might have an influential therapeutic impact on T2DM (Table S6).

3.7. Gene-Ontology (GO) and KEGG Pathway Enrichment Analysis of 124 Common Targets

The GO and KEGG pathway appraisal truly reflects the intersected 124 common targets engaged in the T2DM functional process through relevant targets molecular function (MF), the biological process (BP) in which it takes part, and its cellular localization (CL). We utilized a web-based tool, “DAVID”, to analyze GO and KEGG pathways. With respect to GO, we determined the top 10 MF, BP, and chemical contents based on the proportion of targets that were enriched in those categories (Figure 5), where MF is mostly involved in binding heme, DNA, ligand activation and sequence-specific DNA, enzyme, protein, zinc ion and protein, and steroid hormone receptor activity. The BP in which those targets significantly participate were apoptotic process, cell proliferation, inflammatory response, positive–negative regulation of transcription, transcription initiation from RNA polymerase II promoter, drug response, and oxidation–reduction process. However, the above activities took place in the following top 10 CLs: mitochondrion, extracellular space, endoplasmic reticulum membrane, an integral component of the plasma membrane, extracellular region, nucleoplasm, cytosol, extracellular exosome, plasma membrane and integral component of membrane.
Similarly, the KEGG pathway enrichment analysis was conducted using DAVID online tool on the 124 putative therapeutic targets for T2DM and discovered 9 KEGG pathways at a threshold level of p-value < 0.05. Further, we took into account the rich factor and adjusted the FDR value (Q-value) to define the pathway enrichment analysis, in which the rich factor unveils the amount of pathway enrichment with a substantially lower Q value (Table S7). According to preceding filtering, the PPAR signaling pathway was the most abundantly enriched within the assigned targets, as shown in Figure 6.

3.8. Docking Interaction of a Key Substance with PPAR Signaling Pathway Enriched Targets

Seven targets of the PPAR signaling pathway involved in the T2DM process were subjected to bind with the selected key compound “fisetin tetramethyl ether” using the Glide tools of Schrodinger. We re-docked each co-crystallized ligand to compare each corresponding complex with their synthesized complex pattern. For additional comparison, we also docked metformin with target transcribed protein molecules. Remarkably, we docked thiazolidine derivative and fibrate medicines as of PPARs agonist to determine their activation energy, whereas PPARG (PDB ID: 3EOO) demonstrated a greater binding affinity of −6.092 Kcal/mol than thiazolidine derivatives (−4.849 Kcal/mol), co-crystallized ligand (−4.381 Kcal/mol) and metformin (−3.761 Kcal/mol) via interacting eleven H-bond of ARG-288, SER-342, ILE-262, GLY-258, LEU-340, and GLU-259 residues; and twelve H-phobic bonds of ARG-288, ILE-249, LEU-333, LEU-255, ARG-280, ILE-281, CYS-285, and ILE-341 residues with fisetin tetramethyl ether (Figure 7A).
The docking feature of PPARA (PDB ID: 1K7L) also reflected the advancement stabilization of 1K7L—fisetin tetramethyl ether interface complex having −5.563 Kcal/mol of seven hydrogen bonds (ASN-219, GLU-286, TYR-334, CYS-278, ILE-317, and MET-220); seven hydrophobic interactions (TYR-334, ALA-333, CYS-278, ILE-317, and LEU-321), and one pi-sulfur bond (MET-320) in contrast to co-crystallized ligand (−5.279 Kcal/mol), fibrate drug (−5.204 Kcal/mol), rosiglitazone (−4.156 Kcal/mol) and the standard medicine metformin (−3.145 Kcal/mol) (Figure 7B).
In the same way, PPARD (PDB ID: 5U3Q) had the most apparent binding affinity towards fisetin tetramethyl ether aiming at −5.58 Kcal/mol, which was markedly different from rosiglitazone (−3.372 Kcal/mol), metformin (−3.319 Kcal/mol) and co-crystallized ligand (−2.845 Kcal/mol) molecules binding energy. Of that prominent interaction was worked out forming eight hydrogen bonds of MET-192, ALA-306, ASN-307, PHE-190, GLU-259, LEU-304, and ILE-290; five hydrophobic bonding of LEU-304, ILE-290, LEU-294, and LYS-229; one electrostatic bond of GLU-259; and one pi-sulfur bond of MET-293 residues (Figure 7C).
Following that, fisetin tetramethyl ether formed nine hydrogens and two electrostatic bonds with GLU-62, THR-61, THR-74, THR-75, ASN-60, ASP-72, ASP-78, LYS-59 and ASP-78 residues of 5HZ9, which all have a significant influence on the FABP3 association (PDB ID: 5HZ9) (Figure 7D). The essence of this interaction complex was −3.924 kcal/mol, opposed to rosiglitazone of −3.339 Kcal/mol, co-crystallized ligand of −3.329 Kcal/mol and a standard drug of −2.908 kcal/mol (metformin).
Fisetin tetramethyl ether binds to the FABP4 (PDB ID: 3P6D) receptor, resulting in −3.816 kcal/mol across nine H-bonds (GLU-61, THR-60, ASN-59, VAL-73, and ASP-71) and three H-phobic bonds (VAL-73), which was considerably lower than co-crystallized ligand −3.74 Kcal/mol, rosiglitazone −3.196 Kcal/mol and the conventional medicine metformin −2.81 kcal/mol (Figure 7E).
The active pocket of MMP1 (PDB ID: 1SU3) association with fisetin tetramethyl ether revealed eight H-bonds with HOH-951, GLN-50, SER-172, PRO-173, GLU-39, LYS-36, ASN-43 and six H-phobic bonds with LYS-36, LYS-40, and PRO-95 residues, ultimately exposing a docking score of −4.043 kcal/mol, which was of much lower energy than the competitive agonist rosiglitazone (−3.79 Kcal/mol) and conventional medicine metformin (−2.74 kcal/mol), but nearly identical to the co-crystallized ligand (−4.365 kcal/mol) (Figure 7F). Meanwhile, the NR1H3 (PDB ID: 1UHL)– fisetin tetramethyl ether complex showed more remarkable affinities of −2.967 Kcal/mol than the co-crystallized ligand (−2.22 Kcal/mol), but lower affinity in the case of standard medication metformin (−4.431 Kcal/mol) and rosiglitazone (−4.983 Kcal/mol) (Figure 7G). Table S8 lists the binding pocket residues in detail.
Nevertheless, compared to the corresponding agonist, the co-crystallized ligand and metformin, the vital essential component of MECSW, showed excellent binding affinity, implying that those binding complexes had greater binding stability than their synthetic counterparts. In contrast, this remarkable finding was not lined with the NR1H3 target in terms of metformin and rosiglitazone, which had a comparatively lower binding affinity. Table 3 provides information on the docking score, co-crystallized ligand, thiazolidine derivatives, and fibrate drug nomenclature.

3.9. Quantum Chemistry of Key Ingredients

Density functional theory (DFT) was used to figure out molecular descriptors scoring functions, such as enthalpy, Gibbs free energy, HOMO, LUMO, hardness, and softness energy of key compound and standard medicine, to clearly define their reactive chemical nature, structural properties and regions of the molecules. Whilst the negative values of thermodynamic characteristics of fisetin tetramethyl ether are higher than metformin, the thermodynamic attributes of the key compounds are preferably excellent. The compound’s electron-donating and accepting aptitude are symbolized by the HOMO and LUMO frontiers molecular orbitals, respectively. Similarly, fisetin tetramethyl ether also has the highest negative energy of −0.21123 Kcal/mol, signifying that it is the ideal electron donor ingredient than metformin. Hardness and softness are defined as an energy gap, which are determinant hallmarks of compounds’ chemical reactivity [63]. Appropriately, soft molecules have a smaller energy gap over complex molecules, while complicated molecules have an enormous energy difference. Fisetin tetramethyl ether has a relatively low hardness energy of 0.0759 Kcal/mol than metformin. As the responsiveness of medications accelerates with their softness, this prior articulated pattern is equivalent to molecular softness (Table 4). Figure 8 compares the ground state (HOMO) to the first excited state (LUMO) localization pattern of comparable compounds’ frontier molecular orbitals.

4. Discussion

A synergistic attributable pattern of multi-pathways, multi-components and multi-targets has constantly been observed in medicinal plants. It is more complex to evaluate the herbal treatment efficiency in one drug-one target model due to the presence of multi-diverse chemical ingredients in plants [64]. To metaphrase the traditional potentiality of medicinal plants, network pharmacology research is an excellent technique that prophesied the therapeutic effects of herbal compounds. Such a strategy unearths their target profiles using a computation-based action plan to generate co-module associations of gene–compound–disease that hypothesize the synergetic rules and network modulatory actions of botanical medicinal formula [65]. In this study, we explored the bioactive present in MECSW employing GC-MS techniques and portended the relevant compound’s target profiles by the sequential intersection of their association with disease targets to dig up the mechanistic molecular process of MECSW in the treatment of T2DM. Pertinently, CSW has been widely used in traditional medicine in the south Asian region. However, the exact mechanism in which it exerts antidiabetic activity has not been clearly elucidated yet. Our current study predicted the mechanistic action of compounds extracted from MECSW at the molecular level. It predicted the interaction of target compounds with the genes and pathways involved in T2DM. Yet, further sophisticated experiments still warrant the establishment of this finding.
Pathway compound target (PCT) network revealed that 124 T2DM genes were linked to 26 chemicals (out of 33 chemicals) and 9 pathways in the diagnostic mechanism of T2DM. Similarly, the KEGG pathway also suggested that the nine pathways (four were signaling pathways) influenced the progression and diagnostic process of T2DM, where the gene ontology (GO) resolution of MF, BP and CC was used to perform the functional process of the common targets between T2DM and compound-related genes. The biological process analysis found that common targets were mainly enriched in oxidation–reduction processes backed by the positive–negative regulation of transcription, transcription initiation, signal transduction, etc., which were possibly located in the integral component of the membrane followed by plasma membrane extracellular exosome, cytosol and others; at a molecular level, the genes participated in zinc ion binding, DNA binding, enzyme binding, and heme binding with notably ligand-activated sequence-specific binding potentiality. The relationship between the nine different types of T2DM related pathways is discussed below.
Prostate cancer: insulin-deprived environment in long-term diabetes results in a lower amount of insulin-like growth factor-1 receptor (IGF-1R), and higher plasma IGF-1 is associated with a higher risk of prostate cancer [66,67,68]. Transcriptional misregulation in cancer: a broad range of human cancer had been correlated with the overexpression of insulin receptors (IR) and IGF-1R [69,70,71,72,73]. Moreover, IR-A activated by IGF-2 promotes metastasis and tumor progression, which is highly correlated with hormone-resistant breast cancer [74]. AMPK signaling pathway: AMPK stimulates glucose uptake and regulates the concentration of AMP in pancreatic beta cells in response to blood glucose levels, suggesting the role of AMPK as a key molecular sensor in regulating glucose and lipid metabolism [75,76,77,78]. Thyroid hormone signaling pathway: the thyroid hormone also regulates glucose homeostasis by increasing hepatic glucose output, altering glucose metabolism, decreasing active insulin secretion, and increasing renal insulin clearance [79]. Hyperthyroidism is associated with poor glycemic control, such as hyperglycemia and insulinopenia, responsible for overt diabetes [80]. Metabolic pathways: T2DM-associated metabolic symptoms resulting from chronic inflammation often lead to cell stress and insulin resistance. Insulin deficiency also interferes with some metabolic pathways resulting in hyperglycemia [81]. Calcium signaling pathway: increased Ca2+ influx was found correlated with increased insulin-mediated glucose uptake, indicating the significance of the calcium signaling pathway in T2DM progression [82]. Furthermore, insulin resistance in humans has been found to be positively correlated with increased serum Ca2+ levels [83,84]. Pathways in cancer: hyperglycemia-induced oxidative stress and DNA damage trigger the primary phases of tumor formation [85,86]. Even high glucose concentrations have been reported to modify the gene expression related to proliferation, migration and cell adhesion, thus leading to cancer [87]. Bile secretion: two of the most prominent signaling molecules in the bile secretion pathway are the farnesoid X receptor (FXR) and the G-protein-coupled membrane receptor (TGR5), which play a significant role in the regulation of lipid, glucose and energy metabolism [88]. The alteration of this pathway may lead to T2DM and relevant metabolic disorders [89]. PPAR signaling pathway: direct action of peroxisome proliferator-activated receptors (PPAR)-γ often results in improved insulin sensitivity, which also acts as major drug target for thiazolidinedione (TZD) in the treatment of diabetes mellitus. TZD activated PPAR-γ inhibits the transcription of genes associated with glucose and lipid metabolism [90,91]. The above mentioned nine pathways are integrally linked to the emergence of T2DM and their rich factors indicate that the PPAR signaling pathway has a higher enrichment level. It was reported that the greater rich factor is particularly linked to the greatest extent of enrichment [92]. Farther, this study only focused on the mechanistic action of one specific metabolic pathway involving the hub gene (PPARG) with highest degree value. Thus, other targets with lower value are not considered in this study. However, they may also possess potential to be a good target, which should be further investigated.
The upregulation of PPARG in the PPAR signaling pathway regulates glucose homeostasis through modulating glucose transporter type 4 (GLUT4), beta-glucokinase, and c-Cbl–associated protein (CAP) in the adipose tissue [93]. In addition, PPARG regulates several adipose-tissue-mediated factors, such as TNF-α, adiponectin, leptin, and resistin, which are essential for insulin sensitivity. Hence, the PPARG agonist is important for increasing glucose tolerance by boosting insulin sensitivity and the functionality of beta cells in diabetics [91]. In contrast, insulin sensitivity may be improved by activating PPARG and PPARA [94]. Severe hyperglycemia has been discerned in individuals with the dominant-negative PPARG genetic defect, suggesting a biological link between type-2 diabetes and the PPARG gene [95]. Furthermore, one of the clock genes, “PPARD”, expressed abnormally to the patients with gestational diabetes mellitus (GDM) and T2DM pregnant women [96]. In that way, activation of clock genes, namely PPARA, PPARD and PPARG, would be a better therapeutic approach for T2DM.
In the context of clockwise gene (PPARA, PPARD, and PPARG) activation, we conducted docking assays using the key compound “fisetin tetramethyl ether”, co-crystallized ligand, corresponding gene activation agonist (rosiglitazone and bezfibrate) and the standard medication metformin. As an equivalent to those gene agonists, the key ingredient had excellent activation energy over co-crystallized ligand affinity and metformin binding affinity. This indicates that fisetin tetramethyl ether has the potential to activate PPARs. Several pre-clinical studies also reported fisetin for improved antihyperlipidemic effect and antihyperglycemic effect in streptozotocin-induced diabetic rats [97,98,99]. However, FABP3 and FABP4 also exposed efficient binding affinity and stability with fisetin tetramethyl ether than either metformin or co-crystallized molecules, and MMP1 genes possessed the equivalent or analogous binding interactions and stability positions with respect to the key component and co-crystallized ligand. Compared to the co-crystallized ligand, the NR1H3 gene had more signatory stability, almost equivalent to metformin complex energy. Figure 9 depicts a detailed insight of the PPAR signaling pathway. Due to rosiglitazone’s side effect profile, another contemporary drug, metformin, was used as a control in this study [100]. However, fisetin showed a higher affinity towards PPARG activation than that of both metformin and thiazolidine derivatives.
Additionally, the quantum chemical assessment of the key compound and metformin at the DFT (density functional theory) level was used to verify their chemical reactivity with targets or other species. The HOMO and LUMO energy gaps are essential in influencing medicines’ kinetic stability and chemical reactivity [101]. A large HOMO–LUMO gap can be linked to excellent kinetic stability and inadequate chemical reactivity. Chemical function descriptors, such as hardness and softness, may account for the disparity. Chemical reactivity is exacerbated by their softness [102,103,104]. Simultaneously, our targeted hub component fisetin tetramethyl ether outperformed metformin in robust softness energy. Hence, this progressed as an effective bioactive in treating T2DM. Therefore, fisetin tetramethyl ether might be able to control T2DM. Other studies had also corroborated the antihyperglycemic efficacy of this compound [105,106,107]. However, this research anticipates that fisetin tetramethyl ether could be a multi-target antagonist of T2DM in the network. Even though this research founded the pharmacological effectiveness of MECSW in the treatment of T2DM, more pharmacodynamic and molecular investigations are essential to precisely comprehend the complex synergistic action that underlies the success of MECSW.

5. Conclusions

This study hypothesized the molecular insight of the bioactive presence in MECSW to treat T2DM employing network pharmacology. Importantly, this study’s findings show that MECSW might work through targeting the PPAR signaling pathway, where the key bioactive “fisetin tetramethyl ether” exhibited tremendous impact on this pathway in the regulation towards T2DM. These speculative mechanistic actions were further affirmed by molecular docking simulations, in which this key constituent showed excellent binding affinity with each protein agonist involved in the PPAR signaling cascade and conventional pharmaceuticals. In addition, this compound’s quantum chemistry was also characterized by its superior chemical reactivity over the usual medicine. Ultimately, this study provides empirical testimony to substantiate the treatment effectiveness of MECSW on T2DM and outlines profound insights into the bioactive, interacting potential target, and modes of action of MECSW against T2DM, thus adding to the existing knowledge.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/life12020277/s1, Table S1: Compound related gene, Table S2: T2DM disease related genes, Table S3: Common genes between SEA and STP, Table S4: Common genes between T2DM related targets and compound related overlapping genes, Table S5: Classification of PPI network ordering based on Degree algorithm for 124 common targets.

Author Contributions

Conceptualization, methodology, software, validation, formal analysis, writing—original draft preparation, M.A., B.-B.J. and M.H.U.C.; investigation, data curation, writing—review and editing K.-K.O., T.D. and M.N.U.C.; supervision, project administration, funding acquisition, D.-H.C. All authors have read and agreed to the published version of the manuscript.

Funding

This study was performed with the support of a research grant from Kangwon National University in 2022.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data and materials used are all available in the manuscript.

Acknowledgments

This research was supported by the Department of Bio-Health Convergence, Kangwon National University, Chuncheon 24341, Republic of Korea.

Conflicts of Interest

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

Abbreviations

T2DMType-2 diabetes mellitus
PPIProtein–protein interaction
SEASimilarity ensemble approach
SMILESSimplified molecular input line entry system
STPSwiss target prediction
NR1H3Nuclear receptor subfamily 1 group H member 3
OMIMOnline mendelian inheritance in man
PPARGPeroxisome proliferator activated receptor gamma
MECSWMethanolic extract of Caesalpinia sappan wood
FABP3Fatty acid binding protein 3
PPARAPeroxisome proliferator activated receptor alpha
FABP4Fatty acid binding protein 4
PPARDPeroxisome proliferator activated receptor delta
MMP1Matrix metallopeptidase 1
KEGGKyoto Encyclopedia of Genes and Genomes
FDRFalse discovery rate

References

  1. IDF. Atlas IDF Diabetes Atlas, 10th ed.; 2021; Available online: https://diabetesatlas.org/atlas/tenth-edition/ (accessed on 3 August 2020).
  2. Ling, C.; Groop, L. Epigenetics: A molecular link between environmental factors and type 2 diabetes. Diabetes 2009, 58, 2718–2725. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Annis, A.M.; Caulder, M.S.; Cook, M.L.; Duquette, D. Family history, diabetes, and other demographic and risk factors among participants of the national health and nutrition examination survey 1999–2002. Prev. Chronic Dis. 2005, 2, A19. [Google Scholar] [PubMed]
  4. Sirdah, M.M.; Reading, N.S. Genetic predisposition in type 2 diabetes: A promising approach toward a personalized management of diabetes. Clin. Genet. 2020, 98, 525–547. [Google Scholar] [CrossRef] [PubMed]
  5. Stumvoll, M.; Goldstein, B.J.; Van Haeften, T.W. Type 2 diabetes: Principles of pathogenesis and therapy. Lancet 2005, 365, 1333–1346. [Google Scholar] [CrossRef]
  6. Defronzo, R.A. From the triumvirate to the ominous octet: A new paradigm for the treatment of type 2 diabetes mellitus. Diabetes 2009, 58, 773–795. [Google Scholar] [CrossRef] [Green Version]
  7. Reaven, G.M. Role of insulin resistance in human disease. Diabetes 1988, 37, 1595–1607. [Google Scholar] [CrossRef]
  8. Pinti, M.V.; Fink, G.K.; Hathaway, Q.A.; Durr, A.J.; Kunovac, A.; Hollander, J.M. Mitochondrial dysfunction in type 2 diabetes mellitus: An organ-based analysis. Am. J. Physiol. Metab. 2019, 316, E268–E285. [Google Scholar] [CrossRef]
  9. Wu, P.P.; Kor, C.T.; Hsieh, M.C.; Hsieh, Y.P. Association between End-Stage Renal Disease and Incident Diabetes Mellitus—A Nationwide Population-Based Cohort Study. J. Clin. Med. 2018, 7, 343. [Google Scholar] [CrossRef] [Green Version]
  10. Hsu, Y.C.; Lin, J.T.; Ho, H.J.; Kao, Y.H.; Huang, Y.T.; Hsiao, N.W.; Wu, M.S.; Liu, Y.Y.; Wu, C.Y. Antiviral treatment for hepatitis C virus infection is associated with improved renal and cardiovascular outcomes in diabetic patients. Hepatology 2014, 59, 1293–1302. [Google Scholar] [CrossRef] [Green Version]
  11. Gordon-Dseagu, V.L.Z.; Shelton, N.; Mindell, J. Diabetes mellitus and mortality from all-causes, cancer, cardiovascular and respiratory disease: Evidence from the Health Survey for England and Scottish Health Survey cohorts. J. Diabetes Complicat. 2014, 28, 791–797. [Google Scholar] [CrossRef] [Green Version]
  12. Olokoba, A.B.; Obateru, O.A.; Olokoba, L.B. Type 2 diabetes mellitus: A review of current trends. Oman Med. J. 2012, 27, 269–273. [Google Scholar] [CrossRef] [PubMed]
  13. Mazzola, N. Review of current and emerging therapies in type 2 diabetes mellitus. Am. J. Manag. Care 2012, 18, S17–S26. [Google Scholar] [PubMed]
  14. Zintzaras, E.; Miligkos, M.; Ziakas, P.; Balk, E.M.; Mademtzoglou, D.; Doxani, C.; Mprotsis, T.; Gowri, R.; Xanthopoulou, P.; Mpoulimari, I.; et al. Assessment of the relative effectiveness and tolerability of treatments of type 2 diabetes mellitus: A network meta-analysis. Clin. Ther. 2014, 36, 1443–1453.e9. [Google Scholar] [CrossRef] [PubMed]
  15. El-Kaissi, S.; Sherbeeni, S. Pharmacological Management of Type 2 Diabetes Mellitus: An Update. Curr. Diabetes Rev. 2011, 7, 392–405. [Google Scholar] [CrossRef] [PubMed]
  16. Tsang, M.-W. The Management of Type 2 Diabetic Patients with Hypoglycaemic Agents. ISRN Endocrinol. 2012, 2012, 478120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Palmer, S.C.; Mavridis, D.; Nicolucci, A.; Johnson, D.W.; Tonelli, M.; Craig, J.C.; Maggo, J.; Gray, V.; De Berardis, G.; Ruospo, M.; et al. Comparison of clinical outcomes and adverse events associated with glucose-lowering drugs in patients with type 2 diabetes a meta-analysis. J. Am. Med. Assoc. 2016, 316, 313–324. [Google Scholar] [CrossRef] [PubMed]
  18. Bolen, S.; Feldman, L.; Vassy, J.; Wilson, L.; Yeh, H.C.; Marinopoulos, S.; Wiley, C.; Selvin, E.; Wilson, R.; Bass, E.B.; et al. Systematic review: Comparative effectiveness and safety of oral medications for type 2 diabetes mellitus. Ann. Intern. Med. 2007, 147, 386–399. [Google Scholar] [CrossRef]
  19. Hung, H.Y.; Qian, K.; Morris-Natschke, S.L.; Hsu, C.S.; Lee, K.H. Recent discovery of plant-derived anti-diabetic natural products. Nat. Prod. Rep. 2012, 29, 580–606. [Google Scholar] [CrossRef]
  20. Covington, M.B. Traditional Chinese Medicine in the Treatment of Diabetes. Diabetes Spectr. 2001, 14, 154–159. [Google Scholar] [CrossRef] [Green Version]
  21. Ye, X.P.; Song, C.Q.; Yuan, P.; Mao, R.G. α-Glucosidase and α-Amylase Inhibitory Activity of Common Constituents from Traditional Chinese Medicine Used for Diabetes Mellitus. Chin. J. Nat. Med. 2010, 8, 349–352. [Google Scholar] [CrossRef]
  22. Li, J.; Bai, L.; Li, X.; He, L.; Zheng, Y.; Lu, H.; Li, J.; Zhong, L.; Tong, R.; Jiang, Z.; et al. Antidiabetic potential of flavonoids from traditional Chinese medicine: A review. Am. J. Chin. Med. 2019, 47, 933–957. [Google Scholar] [CrossRef]
  23. Masaenah, E.; Elya, B.; Setiawan, H.; Fadhilah, Z.; Wediasari, F.; Nugroho, G.A.; Elfahmi; Mozef, T. Antidiabetic activity and acute toxicity of combined extract of Andrographis paniculata, Syzygium cumini, and Caesalpinia sappan. Heliyon 2021, 7, e08561. [Google Scholar] [CrossRef] [PubMed]
  24. Muti, A.F.; Pradana, D.L.C.; Rahmi, E.P. Extract of Caesalpinia sappan L. heartwood as food treatment anti-diabetic: A narrative review. IOP Conf. Ser. Earth Environ. Sci. 2021, 755, 012042. [Google Scholar] [CrossRef]
  25. Srilakshmi, V.S.; Vijayan, P.; Raj, P.V.; Dhanaraj, S.A.; Chandrashekhar, H.R. Hepatoprotective properties of Caesalpinia sappan Linn. heartwood on carbon tetrachloride induced toxicity. Indian J. Exp. Biol. 2010, 48, 905–910. [Google Scholar]
  26. Nirmal, N.P.; Rajput, M.S.; Prasad, R.G.S.V.; Ahmad, M. Brazilin from Caesalpinia sappan heartwood and its pharmacological activities: A review. Asian Pac. J. Trop. Med. 2015, 8, 421–430. [Google Scholar] [CrossRef] [Green Version]
  27. Nguyen, H.X.; Nguyen, M.T.T.; Nguyen, T.A.; Nguyen, N.Y.T.; Phan, D.A.T.; Thi, P.H.; Nguyen, T.H.P.; Dang, P.H.; Nguyen, N.T.; Ueda, J.Y.; et al. Cleistanthane diterpenes from the seed of Caesalpinia sappan and their antiausterity activity against PANC-1 human pancreatic cancer cell line. Fitoterapia 2013, 91, 148–153. [Google Scholar] [CrossRef]
  28. Zhu, N.L.; Xu, X.D.; Hu, M.G.; Sun, Z.H.; Wu, T.Y.; Yuan, J.Q.; Wu, H.F.; Tian, Y.; Li, P.F.; Yang, J.S.; et al. Bioactive cassane diterpenoids from the seeds of Caesalpinia sappan. Phytochem. Lett. 2017, 22, 113–116. [Google Scholar] [CrossRef]
  29. Min, B.S.; Cuong, T.D.; Hung, T.M.; Min, B.K.; Shin, B.S.; Woo, M.H. Compounds from the heartwood of Caesalpinia sappan and their anti-inflammatory activity. Bioorg. Med. Chem. Lett. 2012, 22, 7436–7439. [Google Scholar] [CrossRef]
  30. Nguyen, H.X.; Nguyen, N.T.; Dang, P.H.; Thi Ho, P.; Nguyen, M.T.T.; Van Can, M.; Dibwe, D.F.; Ueda, J.Y.; Awale, S. Cassane diterpenes from the seed kernels of Caesalpinia sappan. Phytochemistry 2016, 122, 286–293. [Google Scholar] [CrossRef]
  31. Bae, I.K.; Min, H.Y.; Han, A.R.; Seo, E.K.; Sang, K.L. Suppression of lipopolysaccharide-induced expression of inducible nitric oxide synthase by brazilin in RAW 264. 7 macrophage cells. Eur. J. Pharmacol. 2005, 513, 237–242. [Google Scholar] [CrossRef]
  32. Zeng, K.W.; Yu, Q.; Song, F.J.; Liao, L.X.; Zhao, M.B.; Dong, X.; Jiang, Y.; Tu, P.F. Deoxysappanone B, a homoisoflavone from the Chinese medicinal plant Caesalpinia sappan L., protects neurons from microglia-mediated inflammatory injuries via inhibition of IκB kinase (IKK)-NF-κB and p38/ERK MAPK pathways. Eur. J. Pharmacol. 2015, 748, 18–29. [Google Scholar] [CrossRef]
  33. Shen, J.; Zhang, H.; Lin, H.; Su, H.; Xing, D.; Du, L. Brazilein protects the brain against focal cerebral ischemia reperfusion injury correlating to inflammatory response suppression. Eur. J. Pharmacol. 2007, 558, 88–95. [Google Scholar] [CrossRef] [PubMed]
  34. Hong, Y.J.; Jeong, G.H.; Jeong, Y.H.; Kim, T.H. Evaluation of Antioxidant and Anti-diabetic Effects of Sappan Lignum by Extraction Method. Korea J. Herbol. 2017, 32, 1–7. [Google Scholar]
  35. Chinnala, K.M.; Elsani, M.M.; Nalla, M.K.; Chinnala, K.M. Anti diabetic activity of methanolic extract of Caesalpinia sappan Linn. on alloxan induced diabetes mellitus in rats. Int. J. Exp. Pharmacol. 2015, 5, 65–69. [Google Scholar]
  36. You, E.J.; Khil, L.Y.; Kwak, W.J.; Won, H.S.; Chae, S.H.; Lee, B.H.; Moon, C.K. Effects of brazilin on the production of fructose-2,6-bisphosphate in rat hepatocytes. J. Ethnopharmacol. 2005, 102, 53–57. [Google Scholar] [CrossRef]
  37. Holidah, D.; Dewi, I.P.; Christianty, F.M.; Muhammadiy, N.S.; Huda, N. Antidiabetic and antidyslipidemic activity of secang (Caesalpinia sappan l.) wood extract on diabetic rat. Res. J. Pharm. Technol. 2021, 14, 2801–2806. [Google Scholar] [CrossRef]
  38. Jia, Y.; Wang, H.; Song, Y.; Liu, K.; Dou, F.; Lu, C.; Ge, J.; Chi, N.; Ding, Y.; Hai, W.; et al. Application of a liquid chromatography-tandem mass spectrometry method to the pharmacokinetics, tissue distribution and excretion studies of brazilin in rats. J. Chromatogr. B Anal. Technol. Biomed. Life Sci. 2013, 931, 61–67. [Google Scholar] [CrossRef]
  39. Mueller, M.; Weinmann, D.; Toegel, S.; Holzer, W.; Unger, F.M.; Viernstein, H. Compounds from Caesalpinia sappan with anti-inflammatory properties in macrophages and chondrocytes. Food Funct. 2016, 7, 1671–1679. [Google Scholar] [CrossRef] [Green Version]
  40. Xu, P.; Guan, S.; Feng, R.; Tang, R.; Guo, D. Separation of four homoisoflavonoids from Caesalpinia sappan by high-speed counter-current chromatography. Phytochem. Anal. 2012, 23, 228–231. [Google Scholar] [CrossRef]
  41. Atangwho, I.J.; Ebong, P.E.; Egbung, G.E.; Obi, A.U. Extract of Vernonia amygdalina Del. (African bitter leaf) can reverse pancreatic cellular lesion after alloxan damage in the rat. Aust. J. Basic Appl. Sci. 2010, 4, 711–716. [Google Scholar]
  42. Meliani, N.; Dib, M.E.A.; Allali, H.; Tabti, B. Hypoglycaemic effect of Berberis vulgaris L. in normal and streptozotocin-induced diabetic rats. Asian Pac. J. Trop. Biomed. 2011, 1, 468–471. [Google Scholar] [CrossRef] [Green Version]
  43. Wang, Y.; Hu, B.; Feng, S.; Wang, J.; Zhang, F. Target recognition and network pharmacology for revealing anti-diabetes mechanisms of natural product. J. Comput. Sci. 2020, 45, 101186. [Google Scholar] [CrossRef]
  44. 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, 1185. [Google Scholar] [CrossRef] [PubMed]
  45. Hopkins, A.L. Network pharmacology. Nat. Biotechnol. 2007, 25, 1110–1111. [Google Scholar] [CrossRef] [PubMed]
  46. Hopkins, A.L. Network pharmacology: The next paradigm in drug discovery. Nat. Chem. Biol. 2008, 4, 682–690. [Google Scholar] [CrossRef]
  47. Li, S.; Zhang, B. Traditional Chinese medicine network pharmacology: Theory, methodology and application. Chin. J. Nat. Med. 2013, 11, 110–120. [Google Scholar] [CrossRef]
  48. Zhang, B.; Wang, X.; Li, S. An integrative platform of TCM network pharmacology and its application on a herbal formula, Qing-Luo-Yin. Evid.-Based Complement. Altern. Med. 2013, 2013, 456747. [Google Scholar] [CrossRef] [Green Version]
  49. Berger, S.I.; Iyengar, R. Network analyses in systems pharmacology. Bioinformatics 2009, 25, 2466–2472. [Google Scholar] [CrossRef]
  50. 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]
  51. 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] [Green Version]
  52. 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]
  53. 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]
  54. 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]
  55. 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] [PubMed]
  56. 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]
  57. 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]
  58. Auniq, R.; Chy, M.; Adnan, M.; Roy, A.; Islam, M.; Khan, T.; Hasan, M.; Ahmed, S.; Khan, M.; Islam, N.; et al. Assessment of anti-nociceptive and anthelmintic activities of Vitex Peduncularis Wall. leaves and in silico molecular docking, ADME/T, and PASS prediction studies of its isolated compounds. J. Complement. Med. Res. 2019, 10, 170. [Google Scholar] [CrossRef]
  59. 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]
  60. 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, 35. [Google Scholar] [CrossRef]
  61. Adnan, M.; Chy, M.; Uddin, N.; Kama, A.T.M.; Azad, M.; Kalam, O.; Chowdhury, K.A.A.; Kabir, M.S.H.; Gupta, S.D.; Chowdhury, M. Comparative Study of Piper sylvaticum Roxb. Leaves and Stems for Anxiolytic and Antioxidant Properties Through in vivo, in vitro, and in silico Approaches. Biomedicines 2020, 8, 68. [Google Scholar] [CrossRef] [Green Version]
  62. 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]
  63. Azam, F.; Alabdullah, N.H.; Ehmedat, H.M.; Abulifa, A.R.; Taban, I.; Upadhyayula, S. NSAIDs as potential treatment option for preventing amyloid β toxicity in Alzheimer’s disease: An investigation by docking, molecular dynamics, and DFT studies. J. Biomol. Struct. Dyn. 2017, 36, 2099–2117. [Google Scholar] [CrossRef] [PubMed]
  64. Li, L.; Dai, W.; Li, W.; Zhang, Y.; Wu, Y.; Guan, C.; Zhang, A.; Huang, H.; Li, Y. Integrated Network Pharmacology and Metabonomics to Reveal the Myocardial Protection Effect of Huang-Lian-Jie-Du-Tang on Myocardial Ischemia. Front. Pharmacol. 2020, 11, 2246. [Google Scholar] [CrossRef] [PubMed]
  65. Dong, Y.; Qiu, P.; Zhu, R.; Zhao, L.; Zhang, P.; Wang, Y.; Li, C.; Chai, K.; Shou, D.; Zhao, H. A Combined Phytochemistry and Network Pharmacology Approach to Reveal the Potential Antitumor Effective Substances and Mechanism of Phellinus igniarius. Front. Pharmacol. 2019, 10, 266. [Google Scholar] [CrossRef] [PubMed]
  66. Kasper, J.S.; Liu, Y.; Pollak, M.N.; Rifai, N.; Giovannucci, E. Hormonal profile of diabetic men and the potential link to prostate cancer. Cancer Causes Control 2008, 19, 703–710. [Google Scholar] [CrossRef] [PubMed]
  67. Travis, R.C.; Appleby, P.N.; Martin, R.M.; Holly, J.M.P.; Albanes, D.; Black, A.; Bueno-de-Mesquita, H.B.A.; Chan, J.M.; Chen, C.; Chirlaque, M.-D.; et al. A Meta-analysis of Individual Participant Data Reveals an Association between Circulating Levels of IGF-I and Prostate Cancer Risk. Cancer Res. 2016, 76, 2288–2300. [Google Scholar] [CrossRef] [Green Version]
  68. Feng, X.; Song, M.; Preston, M.A.; Ma, W.; Hu, Y.; Pernar, C.H.; Stopsack, K.H.; Ebot, E.M.; Fu, B.C.; Zhang, Y.; et al. The association of diabetes with risk of prostate cancer defined by clinical and molecular features. Br. J. Cancer 2020, 123, 657–665. [Google Scholar] [CrossRef]
  69. Frasca, F.; Pandini, G.; Scalia, P.; Sciacca, L.; Mineo, R.; Costantino, A.; Goldfine, I.D.; Belfiore, A.; Vigneri, R. Insulin Receptor Isoform A, a Newly Recognized, High-Affinity Insulin-Like Growth Factor II Receptor in Fetal and Cancer Cells. Mol. Cell. Biol. 1999, 19, 3278–3288. [Google Scholar] [CrossRef] [Green Version]
  70. Sciacca, L.; Mineo, R.; Pandini, G.; Murabito, A.; Vigneri, R.; Belfiore, A. In IGF-I receptor-deficient leiomyosarcoma cells autocrine IGF-II induces cell invasion and protection from apoptosis via the insulin receptor isoform A. Oncogene 2002, 21, 8240–8250. [Google Scholar] [CrossRef] [Green Version]
  71. Vella, V.; Pandini, G.; Sciacca, L.; Mineo, R.; Vigneri, R.; Pezzino, V.; Belfiore, A. A novel autocrine loop involving IGF-II and the insulin receptor isoform-A stimulates growth of thyroid cancer. J. Clin. Endocrinol. Metab. 2002, 87, 245–254. [Google Scholar] [CrossRef]
  72. Kalli, K.R.; Falowo, O.I.; Bale, L.K.; Zschunke, M.A.; Roche, P.C.; Conover, C.A. Functional insulin receptors on human epithelial ovarian carcinoma cells: Implications for IGF-II mitogenic signaling. Endocrinology 2002, 143, 3259–3267. [Google Scholar] [CrossRef] [PubMed]
  73. Cox, M.E.; Gleave, M.E.; Zakikhani, M.; Bell, R.H.; Piura, E.; Vickers, E.; Cunningham, M.; Larsson, O.; Fazli, L.; Pollak, M. Insulin receptor expression by human prostate cancers. Prostate 2009, 69, 33–40. [Google Scholar] [CrossRef] [PubMed]
  74. Zelenko, Z.; Gallagher, E.J.; Antoniou, I.M.; Sachdev, D.; Nayak, A.; Yee, D.; LeRoith, D. EMT reversal in human cancer cells after IR knockdown in hyperinsulinemic mice. Endocr. Relat. Cancer 2016, 23, 747–758. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Detimary, P.; Gilon, P.; Henquin, J.C. Interplay between cytoplasmic Ca2+ and the ATP/ADP ratio: A feedback control mechanism in mouse pancreatic islets. Biochem. J. 1998, 333, 269–274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Salt, I.P.; Johnson, G.; Ashcroft, S.J.H.; Hardie, D.G. AMP-activated protein kinase is activated by low glucose in cell lines derived from pancreatic β cells, and may regulate insulin release. Biochem. J. 1998, 335, 533–539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Misra, P.; Chakrabarti, R. The role of AMP kinase in diabetes. Indian J. Med. Res. 2007, 125, 389–398. [Google Scholar] [PubMed]
  78. Jeon, S.M. Regulation and function of AMPK in physiology and diseases. Exp. Mol. Med. 2016, 48, e245. [Google Scholar] [CrossRef]
  79. Potenza, M.; Via, M.A.; Yanagisawa, R.T. Excess thyroid hormone and carbohydrate meta bolism. Endocr. Pract. 2009, 15, 254–262. [Google Scholar] [CrossRef]
  80. Ray, S.; Ghosh, S. Thyroid Disorders and Diabetes Mellitus: Double Trouble. J. Diabetes Res. Ther. 2016, 2, 1–7. [Google Scholar] [CrossRef]
  81. Hameed, I.; Masoodi, S.R.; Mir, S.A.; Nabi, M.; Ghazanfar, K.; Ganai, B.A. Type 2 diabetes mellitus: From a metabolic disorder to an inflammatory condition. World J. Diabetes 2015, 6, 598. [Google Scholar] [CrossRef]
  82. Lanner, J.T.; Katz, A.; Tavi, P.; Sandström, M.E.; Zhang, S.J.; Wretman, C.; James, S.; Fauconnier, J.; Lännergren, J.; Bruton, J.D.; et al. The role of Ca2+ influx for insulin-mediated glucose uptake in skeletal muscle. Diabetes 2006, 55, 2077–2083. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Wang, C.H.; Wei, Y.H. Role of mitochondrial dysfunction and dysregulation of Ca2+ homeostasis in the pathophysiology of insulin resistance and type 2 diabetes. J. Biomed. Sci. 2017, 24, 70. [Google Scholar] [CrossRef] [PubMed]
  84. Yamaguchi, T.; Kanazawa, I.; Takaoka, S.; Sugimoto, T. Serum calcium is positively correlated with fasting plasma glucose and insulin resistance, independent of parathyroid hormone, in male patients with type 2 diabetes mellitus. Metabolism 2011, 60, 1334–1339. [Google Scholar] [CrossRef] [PubMed]
  85. Calle, E.E.; Kaaks, R. Overweight, obesity and cancer: Epidemiological evidence and proposed mechanisms. Nat. Rev. Cancer 2004, 4, 579–591. [Google Scholar] [CrossRef]
  86. Ricart, W.; Fernández-Real, J.M. No decrease in free IGF-I with increasing insulin in obesity-related insulin resistance. Obes. Res. 2001, 9, 631–636. [Google Scholar] [CrossRef]
  87. Masur, K.; Vetter, C.; Hinz, A.; Tomas, N.; Henrich, H.; Niggemann, B.; Zänker, K.S. Diabetogenic glucose and insulin concentrations modulate transcriptom and protein levels involved in tumour cell migration, adhesion and proliferation. Br. J. Cancer 2011, 104, 345–352. [Google Scholar] [CrossRef] [Green Version]
  88. Lefebvre, P.; Cariou, B.; Lien, F.; Kuipers, F.; Staels, B. Role of bile acids and bile acid receptors in metabolic regulation. Physiol. Rev. 2009, 89, 147–191. [Google Scholar] [CrossRef] [Green Version]
  89. Prawitt, J.; Caron, S.; Staels, B. Bile acid metabolism and the pathogenesis of type 2 diabetes. Curr. Diabetes Rep. 2011, 11, 160–166. [Google Scholar] [CrossRef] [Green Version]
  90. Hauner, H. The mode of action of thiazolidinediones. Diabetes Metab. Res. Rev. 2002, 18, S10–S15. [Google Scholar] [CrossRef]
  91. Kim, H.I.; Ahn, Y.H. Role of Peroxisome Proliferator-Activated Receptor-γ in the Glucose-Sensing Apparatus of Liver and β-Cells. Diabetes 2004, 53, 1–6. [Google Scholar] [CrossRef] [Green Version]
  92. Chowdhury, H.U.; Adnan, M.; Oh, K.K.; Cho, D.H. Decrypting molecular mechanism insight of Phyllanthus emblica L. fruit in the treatment of type 2 diabetes mellitus by network pharmacology. Phytomedicine Plus 2021, 1, 100144. [Google Scholar] [CrossRef]
  93. Ahmadian, M.; Suh, J.M.; Hah, N.; Liddle, C.; Atkins, A.R.; Downes, M.; Evans, R.M. PPARγ signaling and metabolism: The good, the bad and the future. Nat. Med. 2013, 19, 557–566. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Wargent, E.; Sennitt, M.V.; Stocker, C.; Mayes, A.E.; Brown, L.; O’Dowd, J.; Wang, S.; Einerhand, A.W.C.; Mohede, I.; Arch, J.R.S.; et al. Prolonged treatment of genetically obese mice with conjugated linoleic acid improves glucose tolerance and lowers plasma insulin concentration: Possible involvement of PPAR activation. Lipids Health Dis. 2005, 4, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Barroso, I.; Gurnell, M.; Crowley, V.E.F.; Agostini, M.; Schwabe, J.W.; Soos, M.A.; Maslen, G.L.; Williams, T.D.M.; Lewis, H.; Schafer, A.J.; et al. Dominant negative mutations in human PPARgamma associated with severe insulin resistance, diabetes mellitus and hypertension. Nature 1999, 402, 880–883. [Google Scholar] [CrossRef]
  96. Pappa, K.I.; Gazouli, M.; Anastasiou, E.; Iliodromiti, Z.; Antsaklis, A.; Anagnou, N.P. Circadian clock gene expression is impaired in gestational diabetes mellitus. Gynecol. Endocrinol. 2013, 29, 331–335. [Google Scholar] [CrossRef]
  97. Prasath, G.S.; Pillai, S.I.; Subramanian, S.P. Fisetin improves glucose homeostasis through the inhibition of gluconeogenic enzymes in hepatic tissues of streptozotocin induced diabetic rats. Eur. J. Pharmacol. 2014, 740, 248–254. [Google Scholar] [CrossRef]
  98. Prasath, G.S.; Subramanian, S.P. Antihyperlipidemic Effect of Fisetin, a Bioflavonoid of Strawberries, Studied in Streptozotocin-Induced Diabetic Rats. J. Biochem. Mol. Toxicol. 2014, 28, 442–449. [Google Scholar] [CrossRef]
  99. Prasath, G.S.; Subramanian, S.P. Fisetin, a tetra hydroxy flavone recuperates antioxidant status and protects hepatocellular ultrastructure from hyperglycemia mediated oxidative stress in streptozotocin induced experimental diabetes in rats. Food Chem. Toxicol. 2013, 59, 249–255. [Google Scholar] [CrossRef]
  100. Chiang, C.-K.; Ho, T.-I.; Peng, Y.-S.; Hsu, S.-P.; Pai, M.-F.; Yang, S.-Y.; Hung, K.-Y.; Wu, K.-D. Rosiglitazone in Diabetes Control in Hemodialysis Patients With and Without Viral Hepatitis Infection. Diabetes Care 2007, 30, 3–7. [Google Scholar] [CrossRef] [Green Version]
  101. 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]
  102. 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] [PubMed] [Green Version]
  103. 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]
  104. 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] [PubMed] [Green Version]
  105. Ranga Rao, R.; Tiwari, A.K.; Prabhakar Reddy, P.; Suresh Babu, K.; Ali, A.Z.; Madhusudana, K.; Madhusudana Rao, J. New furanoflavanoids, intestinal α-glucosidase inhibitory and free-radical (DPPH) scavenging, activity from antihyperglycemic root extract of Derris indica (Lam.). Bioorg. Med. Chem. 2009, 17, 5170–5175. [Google Scholar] [CrossRef]
  106. Maher, P.; Dargusch, R.; Ehren, J.L.; Okada, S.; Sharma, K.; Schubert, D. Fisetin Lowers Methylglyoxal Dependent Protein Glycation and Limits the Complications of Diabetes. PLoS ONE 2011, 6, e21226. [Google Scholar] [CrossRef]
  107. Miyata, T.; Van Ypersele de Strihou, C.; Imasawa, T.; Yoshino, A.; Ueda, Y.; Ogura, H.; Kominami, K.; Onogi, H.; Inagi, R.; Nangaku, M.; et al. Glyoxalase I deficiency is associated with an unusual level of advanced glycation end products in a hemodialysis patient. Kidney Int. 2001, 60, 2351–2359. [Google Scholar] [CrossRef] [Green Version]
Scheme 1. The current framework implies how a comprehensive network pharmacology protocol was used to uncover key targets, possible functions, and pharmacological pathways in T2DM patients treated with MECSW.
Scheme 1. The current framework implies how a comprehensive network pharmacology protocol was used to uncover key targets, possible functions, and pharmacological pathways in T2DM patients treated with MECSW.
Life 12 00277 sch001
Figure 1. GC-MS peak of the methanolic extract of C. sappan and an indication of the key bioactive (fisetin tetramethyl ether).
Figure 1. GC-MS peak of the methanolic extract of C. sappan and an indication of the key bioactive (fisetin tetramethyl ether).
Life 12 00277 g001
Figure 2. Overlapping genes: (A) compounds related genes from SEA and STP; (B) T2DM targets from DisGeNet, Malacards, and OMIM.
Figure 2. Overlapping genes: (A) compounds related genes from SEA and STP; (B) T2DM targets from DisGeNet, Malacards, and OMIM.
Life 12 00277 g002
Figure 3. Interactions between proteins (PPI) of common interacting genes utilizing the degree algorithm.
Figure 3. Interactions between proteins (PPI) of common interacting genes utilizing the degree algorithm.
Life 12 00277 g003
Figure 4. Interactions between MECSW components and common T2DM targets along with their pathways through the pathway compound target (PCT) network.
Figure 4. Interactions between MECSW components and common T2DM targets along with their pathways through the pathway compound target (PCT) network.
Life 12 00277 g004
Figure 5. Gene ontology (GO) exploration of compounds’ and T2DM’s shared targets.
Figure 5. Gene ontology (GO) exploration of compounds’ and T2DM’s shared targets.
Life 12 00277 g005
Figure 6. KEGG assessment bubble map showing strongly enriched pathways for common targets.
Figure 6. KEGG assessment bubble map showing strongly enriched pathways for common targets.
Life 12 00277 g006
Figure 7. PPAR signaling pathway genes and MECSW key molecules binding interaction. (A) 3E00—fisetin tetramethyl ether, (B) 1K7L—fisetin tetramethyl ether, (C) 5U3Q—fisetin tetramethyl ether, (D) 5HZ9—fisetin tetramethyl ether, (E) 3P6D—fisetin tetramethyl ether, (F) 1SU3—fisetin tetramethyl ether, and (G) 1UHL—fisetin tetramethyl ether.
Figure 7. PPAR signaling pathway genes and MECSW key molecules binding interaction. (A) 3E00—fisetin tetramethyl ether, (B) 1K7L—fisetin tetramethyl ether, (C) 5U3Q—fisetin tetramethyl ether, (D) 5HZ9—fisetin tetramethyl ether, (E) 3P6D—fisetin tetramethyl ether, (F) 1SU3—fisetin tetramethyl ether, and (G) 1UHL—fisetin tetramethyl ether.
Life 12 00277 g007aLife 12 00277 g007b
Figure 8. Plotting the energy of key compound and conventional drug on the frontier molecular orbital.
Figure 8. Plotting the energy of key compound and conventional drug on the frontier molecular orbital.
Life 12 00277 g008
Figure 9. Depiction of the PPAR signaling pathway with the enriched genes in this study are shown by the red node.
Figure 9. Depiction of the PPAR signaling pathway with the enriched genes in this study are shown by the red node.
Life 12 00277 g009
Table 1. A list of 33 chemical components noted by GC-MS revealed from methanolic extract of C. sappan wood (MECSW).
Table 1. A list of 33 chemical components noted by GC-MS revealed from methanolic extract of C. sappan wood (MECSW).
SL. No.RT Time (min)Area (%)PubChem CIDChemical FormulaCompound ClassCompounds Name
15.3660.225363448C7H12O3Fatty acyls(3E)-5-Hydroxy-2-methyl-3-hexenoic acid
25.6831.35055C8H8O3Phenol esters1,3-Benzenediol, Monoacetate
36.6060.537342C6H12O2Carboxylic acids and derivativesEthyl 2-Methylpropanoate
46.8560.75054C6H6O2Phenols1,3-Benzenediol
57.5292.1462378C11H18OOrganooxygen compounds3-Methyl-2-Pentyl-2-Cyclopenten-1-One
68.1160.5281750C9H12O3Benzene and substituted derivatives2,4-Dimethoxybenzyl alcohol
78.2310.9291477C27H44OSteroids and steroid derivativesCholest-4-En-3-One
88.511.26785C6H6O2PhenolsHydroquinone
98.7310.965284421C19H34O2Fatty acyls9,12-Octadecadienoic acid, methyl ester
108.8180.1521206C16H32O2Fatty acylsTetradecanoic acid, 12-Methyl-, Methyl Ester
118.9143.955280450C18H32O2Fatty acyls9,12-Octadecanoic Acid(Z,Z)-
128.9813.345363372C18H32O2Fatty acylsZ,Z-10,12-Hexadecadien-1-ol acetate
139.3750.1841133C14H24OOrganooxygen compounds2(1H)-Naphthalenone, octahydro-4a-methyl-7-(1-methylethyl)-, (4a.alpha.,7.beta.,8a.beta.)-
1410.010.227210C17H18N2O2Phenols2-((E)-[((E)-2-([(E)-(2-Hydroxyphenyl)Methylidene]Amino)Propyl)Imino]Phenol or, Disalicylalpropylenediimine
1510.190.475364132C14H25N3O2Allyl-type 1,3-dipolar organic compounds2-[(4E)-4-Hexenyl]-6-nitrocyclohexanone dimethylhydrazone
1610.650.664506C15H11N3O3Benzodiazepines7-Nitro-1,3-dihydro-5-phenyl-2H-1,4-benzodiazepin-2-one or, Neozepam
1710.730.52610113C11H13NO3Indoles and derivativesIndole-2-one, 2,3-dihydro-N-hydroxy-4-methoxy-3,3-dimethyl-
1810.883.94610177C13H9N3DiazinesPyrido[2,3-d]pyrimidine, 4-phenyl-
1911.192.935116643C24H37NOCarboxylic acids and derivativesAcetamide, 2-(adamantan-1-yl)-N-(1-adamantan-1-ylethyl)- or, 2-adamantanyl-N-(adamantanylethyl)acetamide
2011.281.225379033C20H22O42-arylbenzofuran flavonoidsPhenol, 4-[2,3-dihydro-7-methoxy-3-methyl-5-(1-propenyl)-2-benzofuranyl]-2-methoxy- or, Dehydrodiisoeugenol
2111.6513.935379034C19H18O3SPhenol ethersPropenone, 1-[5-(3-hydroxy-3-methyl-1-butynyl)-2-thienyl]-3-(4-methoxyphenyl)-
2211.91.78346948C15H13N PyrrolesIndolizine, 2-(4-methylphenyl)-
2311.996.15610018C17H18N4O3Pyrans6-Amino-5-cyano-4-(5-cyano-2,4-dimethyl-1H-pyrrol-3-yl)-2-methyl-4H-pyran-3-carboxylic acid ethyl ester
2412.5111.82631121C20H26N2O3Plumeran-type alkaloidsAspidodispermine, O-methyl-
2512.633.08631095C19H18O6Flavonoids4H-1-Benzopyran-4-one, 3,5,7-trimethoxy-2-(4-methoxyphenyl)- or, Tetramethylkaempferol
2612.733.38631112C19H29B2NO4TetralinsNadolol di-methylboronic acid
2712.855.34631171C19H18O6Flavonoids4H-1-Benzopyran-4-one, 2-(3,4-dimethoxyphenyl)-3,7-dimethoxy- or, Fisetin tetramethyl ether
2812.9512.81284060C12H18N6O6Diazines1,3-Dimethyl-5,6-dicarbethoxy-5,6,7,8-tetrahydro-6,7-diazalumazine
2913.891.75541560C28H48OSteroids and steroid derivativesCholestane, 3,4-epoxy-2-methyl-, (2.alpha.,3.alpha.,4.alpha.,5.alpha.)-
3014.091.6591733922C12H17NO2Isoindoles and derivatives4-Cyclohexene-1,2-dicarboximide, N-butyl-, cis-
3114.655.69457801C29H50OSteroids and steroid derivativesStigmast-5-en-3-ol, (3beta,24S)- or, Clionasterol
3215.270.483247C15H13NIndoles and derivatives1H-Indole, 5-methyl-2-phenyl-
3315.660.68610182C15H13NQuinolines and derivativesBenzo[h]quinoline, 2,4-dimethyl-
Table 2. Documentation of drug-likeness character of 33 compounds from MECSW.
Table 2. Documentation of drug-likeness character of 33 compounds from MECSW.
Compounds NameMWHBAHBDMLogPNumber of ViolationsBioavailability
<500<10≤5≤4.15≤1>0.1
(3E)-5-Hydroxy-2-methyl-3-hexenoic acid144.17320.6400.85
1,3-Benzenediol, Monoacetate152.15311.3200.55
Ethyl 2-Methylpropanoate116.16201.2700.55
1,3-Benzenediol110.11220.7900.55
3-Methyl-2-Pentyl-2-Cyclopenten-1-One166.26102.4900.55
2,4-Dimethoxybenzyl alcohol168.19310.9200.55
Cholest-4-En-3-One384.64106.2310.55
Hydroquinone110.11220.7900.55
9,12-Octadecadienoic acid, methyl ester294.47204.710.55
Tetradecanoic acid, 12-Methyl-, Methyl Ester256.42204.1910.55
9,12-Octadecanoic Acid (Z,Z)-280.45214.4710.85
Z,Z-10,12-Hexadecadien-1-ol acetate280.45204.4710.55
2(1H)-Naphthalenone, octahydro-4a-methyl-7-(1-methylethyl)-, (4a.alpha.,7.beta.,8a.beta.)-208.34103.4100.55
Disalicylalpropylenediimine282.34421.8100.55
2-[(4E)-4-Hexenyl]-6-nitrocyclohexanone dimethylhydrazone267.37301.8200.55
Neozepam281.27410.900.55
Indole-2-one, 2,3-dihydro-N-hydroxy-4-methoxy-3,3-dimethyl-207.23311.400.55
4-phenylpyrido[2,3-d]pyrimidine207.23302.0300.55
2-adamantanyl-N-(adamantanylethyl)acetamide355.56115.210.55
Dehydrodiisoeugenol326.39412.7300.55
Propenone, 1-[5-(3-hydroxy-3-methyl-1-butynyl)-2-thienyl]-3-(4-methoxyphenyl)-326.41312.4300.55
2-(4-methylphenyl) Indolizine207.27003.3200.55
6-Amino-5-cyano-4-(5-cyano-2,4-dimethyl-1H-pyrrol-3-yl)-2-methyl-4H-pyran-3-carboxylic acid ethyl ester326.3552-0.300.56
Aspidodispermine, O-methyl-342.43412.0200.55
Tetramethylkaempferol342.34600.9400.55
Nadolol di-methylboronic acid357.06501.2100.55
Fisetin tetramethyl ether342.34600.9400.55
1,3-Dimethyl-5,6-dicarbethoxy-5,6,7,8-tetrahydro-6,7-diazalumazine342.31721.2210.55
Cholestane, 3,4-epoxy-2-methyl-, (2.alpha.,3.alpha.,4.alpha.,5.alpha.)-400.68106.6810.55
4-Cyclohexene-1,2-dicarboximide, N-butyl-, cis-207.27201.6800.55
Clionasterol414.71116.7310.55
5-methyl-2-phenyl-1H-Indole207.27013.3200.55
Benzo[h]quinoline, 2,4-dimethyl-207.27103.3200.55
MW = molecular weight: HBA = hydrogen bond acceptor; HBD = hydrogen bond donor.
Table 3. Scores for bioactive compound, conventional medication, and their co-crystallized ligands docking with substantially enriched pathway genes of T2DM.
Table 3. Scores for bioactive compound, conventional medication, and their co-crystallized ligands docking with substantially enriched pathway genes of T2DM.
GenePDB IDCompoundDocking Score
PPARG3E00Fisetin Tetramethyl Ether−6.092
Rosiglitazone−4.73
Metformin *−3.761
2-Chloro-5-Nitro-N-Phenylbenzamide−4.381
PPARA1K7LFisetin Tetramethyl Ether−5.563
Bezfibrate−5.204
Rosiglitazone−4.156
Metformin *−3.145
2-(1-Methyl-3-Oxo-3-Phenyl-Propylamino)-3-{4-[2-(5-Methyl-2-Phenyl-Oxazol-4-Yl)-Ethoxy]-Phenyl}-Propionic acid−5.279
PPARD5U3QFisetin Tetramethyl Ether−5.58
Bezfibrate−4.06
Rosiglitazone−3.372
Metformin *−3.319
6-(2-{[([1,1′-Biphenyl]-4-Carbonyl)(propan-2-Yl) amino]methyl}phenoxy)hexanoic acid−2.845
FABP35HZ9Fisetin Tetramethyl Ether−3.924
Rosiglitazone−3.339
Metformin *−2.908
6-Chloranyl-2-Methyl-4-Phenyl-Quinoline-3-Carboxylic acid−3.329
FABP43P6DFisetin Tetramethyl Ether−3.816
Rosiglitazone−3.196
Metformin *−2.81
3-(4-Methoxy-3-Methylphenyl)propanoic acid−3.75
MMP11SU3Fisetin Tetramethyl Ether−4.043
Rosiglitazone−3.79
Metformin *−2.74
4-(2-Hydroxyethyl)-1-Piperazine ethanesulfonic acid−4.365
NR1H31UHLFisetin Tetramethyl Ether−2.967
Rosiglitazone−4.983
Metformin *−4.431
N-(2,2,2-Trifluoroethyl)-N-{4-[2,2,2-Trifluoro-1-Hydroxy-1-(Trifluoromethyl) ethyl]phenyl}benzenesulfonamide−2.22
* = Standard medicine.
Table 4. Key chemical and conventional medication quantum properties.
Table 4. Key chemical and conventional medication quantum properties.
CompoundEnthalpyGibbs Free EnergyHOMOLUMOEgηS
Fisetin Tetramethyl Ether−1186.191−1186.19096−0.21123−0.05941−0.1518−0.0759−13.173
Metformin−432.845−432.844572−0.202150.02114−0.2233−0.1116−8.957
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Adnan, M.; Jeon, B.-B.; Chowdhury, M.H.U.; Oh, K.-K.; Das, T.; Chy, M.N.U.; Cho, D.-H. Network Pharmacology Study to Reveal the Potentiality of a Methanol Extract of Caesalpinia sappan L. Wood against Type-2 Diabetes Mellitus. Life 2022, 12, 277. https://doi.org/10.3390/life12020277

AMA Style

Adnan M, Jeon B-B, Chowdhury MHU, Oh K-K, Das T, Chy MNU, Cho D-H. Network Pharmacology Study to Reveal the Potentiality of a Methanol Extract of Caesalpinia sappan L. Wood against Type-2 Diabetes Mellitus. Life. 2022; 12(2):277. https://doi.org/10.3390/life12020277

Chicago/Turabian Style

Adnan, Md., Byeong-Bae Jeon, Md. Helal Uddin Chowdhury, Ki-Kwang Oh, Tuhin Das, Md. Nazim Uddin Chy, and Dong-Ha Cho. 2022. "Network Pharmacology Study to Reveal the Potentiality of a Methanol Extract of Caesalpinia sappan L. Wood against Type-2 Diabetes Mellitus" Life 12, no. 2: 277. https://doi.org/10.3390/life12020277

APA Style

Adnan, M., Jeon, B. -B., Chowdhury, M. H. U., Oh, K. -K., Das, T., Chy, M. N. U., & Cho, D. -H. (2022). Network Pharmacology Study to Reveal the Potentiality of a Methanol Extract of Caesalpinia sappan L. Wood against Type-2 Diabetes Mellitus. Life, 12(2), 277. https://doi.org/10.3390/life12020277

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