Natural Compounds from Hatikana Extract Potentiate Antidiabetic Actions as Displayed by In Vivo Assays and Verified by Network Pharmacological Tools

Rahman, Md. Atiar; Uddin, Md. Nazim; Babteen, Nouf Abubakr; Alnajeebi, Afnan M.; Zakaria, Zainul Amiruddin; Aboelenin, Salama Mostafa

doi:https://doi.org/10.1155/2021/6978450

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Data Availability Ethical Approval Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 6978450 | https://doi.org/10.1155/2021/6978450

Natural Compounds from Hatikana Extract Potentiate Antidiabetic Actions as Displayed by In Vivo Assays and Verified by Network Pharmacological Tools

Md. Atiar Rahman,¹Md. Nazim Uddin,²Nouf Abubakr Babteen,³Afnan M. Alnajeebi,³Zainul Amiruddin Zakaria,^4,5and Salama Mostafa Aboelenin⁶

Academic Editor: Kazim Husain

Received18 May 2021

Accepted29 Sept 2021

Published23 Oct 2021

Abstract

Background. Hatikana is a traditional medicinal plant used to treat inflammation, urolithiasis, goiter, cancer, wounds and sores, gastrointestinal, tumor, tetanus, arthritis, hepatic damage, neurodegeneration, and other ailments. The goal of this study is to investigate the antidiabetic properties of Hatikana extract (HKEx) and to construct the effects of its natural constituents on the genes and biochemical indices that are connected with them. Methods. HKEx was evaluated using GC-MS and undertaken for a three-week intervention in fructose-fed STZ-induced Wistar albino rats at the doses of HKEx50, HKEx100, and HKEx200 mg/kg bw. Following intervention, blood serum was examined for biochemical markers, and liver tissue was investigated for the mRNA expression of catalase (CAT), glutathione peroxidase (GPx), and superoxide dismutase (SOD1) by RTPCR analysis. Most abundant compounds (oleanolic acid, 7α, 28-olean diol, and stigmasterol) from GC-MS were chosen for the network pharmacological assay to verify function-specific gene-compound interactions using STITCH, STRING, GSEA, and Cytoscape plugin cytoHubba. Results. In vivo results showed a significant () decrease of blood sugar, aspartate aminotransferase (AST), alanine aminotransferase (ALT), creatinine kinase (CK-MB), and lactate dehydrogenase (LDH) and increase of liver glycogen, glucose load, and serum insulin. Out of three antioxidative genes, catalase (CAT) and superoxide dismutase (SOD1) were found to be few fold increased. Oleanolic acid and stigmasterol were noticed to strongly interact with 27 target proteins. Oleanolic acid interacted with the proteins AKR1B10, CASP3, CASP8, CYP1A2, CYP1A2, HMGB1, NAMPT, NFE2L2, NQO1, PPARA, PTGIR, TOP1, TOP2A, UGT2B10, and UGT2B11 and stigmasterol with ABCA1, ABCG5, ABCG8, CTSE, HMGCR, IL10, CXCL8, NR1H2, NR1H3, SLCO1B1, SREBF2, and TNF. Protein-protein interaction (PPI) analysis revealed the involvement of 25 target proteins out of twenty seven. Cytoscape plugin cytoHubba identified TNF, CXCL8, CASP3, PPARA, SREBF2, and IL10 as top hub genes. Pathway analysis identified 31 KEGG metabolic, signaling, and immunogenic pathways associated with diabetes. Notable degree of PPI enrichment showed that SOD1 and CAT are responsible for controlling signaling networks and enriched pathways. Conclusion. The findings show that antioxidative genes have regulatory potential, allowing the HKEx to be employed as a possible antidiabetic source pending further validation.

1. Introduction

Diabetes mellitus (DM), a chronic metabolic disorder characterized by hyperglycemia resulting from increased hepatic glucose production, impaired insulin production by pancreatic β cells and insulin resistance [1]. Long-term impact of DM contributes to the damage of the renal, cardiovascular, endocrine and nervous tissues, and organs [2]. The use of plant-based medicines has recently been attracted much due to the unaffordability, noncompliance, and adverse effects of synthetic drugs. Network pharmacology has been used as a very powerful instrument to analyze the molecular pathways in order to fully determine the pharmacological effects of plant-derived drugs [3, 4]. This tool has drawn considerable interest to conduct systematic pharmacological studies on the relationship between biological processes and the treatment of herbal medicine. Indeed, network pharmacology methods were developed with the goal of prioritizing disease-associated genes, predicting the target profiles and pharmacological actions of herbal compounds, revealing drug-gene-disease comodule associations, screening synergistic multicompounds from herbal formulae in a high-throughput manner, and interpreting combinatorial rules and network regulations [5]. Studying the basic biological knowledge of plant-based drugs on the basis of network interaction will provide indepth insight or scientific evidence for the discovery of alternative medicine and enable us to understand their pharmacological function at the biomolecule stage [6]. The usefulness of network pharmacology has made it a dynamic tool for investigating the interactions between drugs and biological systems, including human organs, diseases, metabolic processes, and target proteins [7, 8]. Natural compounds from numerous plants have been approached for network pharmacology to be explored as new drugs in the current drug discovery method.

Leea macrophylla, locally named as Hatikana or Hostikorna, is an extraordinary medicinal plant widely distibuted in Bangladesh, India, Srilanka, Nepal, and China. Leaves and roots of Leea macrophy have been studied and so far evident as antioxidant [9], anti-inflammatory [10], antithrombotic [11] antilithiatic [12], and protective for goiter, gastric tumor, lipoma, and tetanus [13, 14]. In this study, the antidiabetic effect of Leea macrophylla extract (HKEx) has been evaluated using the animal model and the combination of pharmacology, and pharmacodynamics has been successfully used to clarify the in vivo effect at the stage of the molecular network [8].

2. Materials and Methods

2.1. Plant Material and Extraction

Hatikana samples (Leea macrophylla root samples) were collected from the cultivation plots of the Bangladesh Council of Science and Industrial Research (BCSIR), Rajshahi, in July 2015. Dr. Sheikh Bokhtear Uddin, a taxonomist and professor at the Department of Botany, University of Chittagong, identified the plant. The nomenclature of the plant has been confirmed from International Plant Name Index (IPNI) using http://powo.science.kew.org/. A sample specimen of the plant has been preserved in the institutional herbarium as an accession number SBU1109CTGUH for future reference. Hatikana samples were thoroughly washed to remove debris and earth remains followed by chopping into small pieces and shade-dried for about 15 days. Dried pieces were subsequently ground into powder using mechanical grinder (Miyako, Model No: DL-718, China). Resulting powder (800 g) was defatted by n-hexane (2.0 L) for two days of maceration with occasional stirring. The fat-free residue was further extracted with maceration in ethanol (99.9%) for 3 days, and collected supernatant was evaporated by using a rotary evaporator (RE 200, Bibby Sterling Ltd., UK) set at reduced pressure within 50-55°C. The concentrated crude extract of Hatikana (HKEx), semisolid dark-brown extract (w/W, 4.7%), was left in the air in petridish and allowed for a complete solvent evaporation and stored at 4°C for further use.

2.2. Phytochemical Analyses of the Extract

HKEx was assessed with GC-MS using the electron impact ionization (EI) process on a gas chromatograph (GC-17A, Shimadzu Corporation, Kyoto, Japan) coupled to a mass spectrometer (GC-MS TQ 8040, Shimadzu Corporation, Kyoto, Japan). A fused capillary column (Rxi-5 ms; 0.25 m film thickness, 0.32 mm internal diameter) coated with silica was used for the separation. The temperature of the inlet was set to 260°C, and the temperature of the oven was programmed to be 70°C (0 min); 10°C, 150°C (5 min); 12°C, 200°C (15 min); and 12°C, 220°C (5 min). At a constant pressure of 90 KPa, the column flow rate was 0.6 mL/min helium gas. The temperature of the auxiliary gas flow rate (GC to MS interface) was set at 280°C. With a scanning range of 40-350 amu, the MS was set in scan mode; the ionization mode was type EI (electron ionization). The mass range was set at the range of 50-550 m/z. For GC/MS analysis, the prepared one microliter (1 μL) sample was injected and run in spiltless modes. The overall runtime of the GC-MS was 35 minutes. In the GC-MS library version NIST 08-S, all peak areas were compared to the database.

2.3. Experiments in the Animal Model

In an animal experiment, thirty (age: six-seven weeks, average body weight: ) fructose-fed (10% fructose solution, an appropriate concentration to induce partial pancreatic dysfunction by STZ) Wistar albino rats were randomly classified into normal control (NC, received distilled water), diabetes control (DBC, received STZ), and three treatment groups (HKEx at 50, 100, and 200 mg/kg bw). Following a week-long acclimatization in a confined setting (12-hour day-light cycle, room temperature ; humidity 60-65%), the animals were injected by a single dose of freshly prepared streptozotocin (60 mg/kg bw, dissolved in 0.1 M citrate buffer). After a week, animals with a nonfasting blood glucose level of around 300 mg/dL were considered to be diabetic and chosen for the experiment. Under the institutional ethical guidelines of the University of Chittagong (Faculty of Biological Sciences, AERB/CUBS/2019-05), the handling and treatment of the experimental animals were guaranteed.

2.4. Biochemical and Molecular Analyses of Collected Blood

At the third week of animal intervention, a single dose of glucose solution (2 g/kg body weight) was administered to assess the glucose tolerance capacity by measuring the subsequent blood glucose levels at 0, 30, 60, 90, and 120 min after administration by the tail prick method. At the end of the intervention, halothane anesthesia was used to sacrifice the animals, and blood and organs were collected. The entire blood, collected by cardiac puncture, was centrifuged for 15 min at 1100 g, and the separated serum was retained at -30°C for further analysis of lipid profile, aspartate transaminase (AST), alanine transaminase (ALT), creatinine kinase (CK-MB), lactate dehydrogenase (LDH), creatinine, and uric acids. Liver glycogen concentrations were measured by a phenol-sulfuric acid method as described by Lo et al. [15]. Total RNA extraction and mRNA expression of liver antioxidative enzymes: glutathione peroxidase, CAT, and SOD1 were accomplished using the protocol described by Mawa et al. [16]. Briefly, 500 mg of each liver sample was added shortly after dissection to 1 mL of TRIzol Reagent (Invitrogen, Ambion, USA) and homogenized with a tissue homogenizer (IKA, UltraTurrax, Germany). One mL of tissue homogenate was transferred to a microfuge tube for total RNA extraction by adding 0.2 mL chloroform and vigorously vortexed for 15 sec before incubating for 3 min at room temperature. The aqueous phase containing RNA was transferred to new tubes after centrifugation at 12,000 g for 15 min at 4°C, and the RNA was precipitated by mixing the aqueous phase with 0.5 mL isopropyl alcohol and incubating at room temperature for 10 min. After centrifuging the RNA pellets at 12,000 g for 10 min at 4°C, they were washed by mixing and vortexing with 1 mL 75 percent ethanol and recentrifuged at 7500 g for 5 min at 4°C. A Reverse Transcription System Kit (Promega) and aSimpliAMP Thermal Cycler were used to generate cDNA from 2 g of total RNA (Life Technologies, Applied Biosystem, USA). A Reverse Transcription System Kit (Promega) and aSimpliAMP Thermal Cycler were used to generate cDNA from 2 g of total RNA (Life Technologies, Applied Biosystem, USA). Briefly, total RNA was activated at 70°C for 10 min, and a 20 μL reaction mix was made with 2.5 μL MgCl₂, 4 μL of reverse transcription 10× buffer, 1 μL of 10 mM dNTTP mixture, ribonuclease inhibitor (0.5 μL), 1 μL oligo DT, 1 μL of GoScript reverse transcriptase enzyme, 1 ng RNA, and nuclease-free water to a final volume of 20 μL. The reaction mixture was then incubated at 25°C for 5 min and 42°C for 60 min followed by incubation at 70°C for 15 min. For PCR amplification, the cDNA was diluted up to 50 l with nuclease-free water. The RT-qPCR was performed with a SYBR Green PCR Master Mix kit (Promega) and specific primers for the antioxidant enzymes-related genes: CAT and SOD1. Each 18 μL reaction contained 9 μL of master mix, 1.5 μL of forward primer (10 μM/liter), 1.5 μL of reverse primer (10 μM/liter), and 6 μL of cDNA. The cycling parameters were as follows: initial incubation at 95°C for 3 min, 40 cycles of 95°C for 30 s, 51°C for 15 s, and 72°C for 30s, with a final extension at 72°C for 10 min. The specificity of the obtained products was confirmed by analysis of the amplified product dissociation curves. The data obtained were analyzed using the 2^-∆∆CT method. Within each sample, the target genes were then normalized to GAPDH.

2.5. Target Identification for the Computational Model

Three compounds oleanolic acid, stigmasterol, and 7 α, 28-olean diol from the list of GC-MS characterized compounds were selected for bioinformatic-based network pharmacological analyses. These compounds were preferred because they were also identified by other authors via vacuum liquid chromatography (VLC) analysis [17] from Hatikana as well as their abundance of occurrence in GC-MS analysis.

2.6. Bioactive Compound-Target Protein Network Construction

Target compounds for network pharmacology were selected from the GC-MS analysis data and the bioactive compounds isolated through vacuum liquid chromatography (VLC) for fractionation over silica gel (kieselgel 60H, 70-230 mesh) with a mobile phase of n-hexane: ethyl acetate [17]. STITCH 5 (http://stitch.embl.de/, ver.5.0) was used on the basis of network pharmacology-based prediction to classify target proteins associated with bioactive phytochemicals that were found in HKEx [18]. It calculates a score for each pair of protein-chemical interactions. Chemical names of bioactive compounds (oleanolic acid, NZ-15; 7 α, 28-olean diol, NZ-38; and stigmasterol, NZ-14) were put into STITCH 5 together to match their potential targets, with the organism selected as “Homo sapiens” and medium required interaction score being 0.4 (threshold level). Twenty-seven medium confidence target proteins were predicted, 15 for oleanolic acid (NZ-15), and 12 for stigmasterol (NZ-14), but no target was identified for the third compound, 7 alpha, 28-olean diol. The compound targets with no association with the interactions of compound proteins were not taken into account for further analysis. Compound-protein interaction data were imported into Cytoscape 3.6.1 software to compare and construct a compound-protein interaction network [19].

2.7. Construction of Protein-Protein Interaction (PPI) Network of the Predicted Genes

A PPI network of predicted genes was created through a search tool for the recovery of the interacting gene (STRING) database (https://string-db.org/cgi/input.pl; STRING-DB v11.0) [20]. The rank of the target proteins based on the degree of interactions in the PPI network was defined using the cytoHubba plugin Cytoscape [21]. The obtained protein interaction data of 27 target proteins were imported into Cytoscape 3.6.1 software to construct a PPI protein interaction network [22].

2.8. Gene Ontology (GO) and Kyoto Encyclopedia of Genomes and Genes (KEGG) Pathway Target Protein’s Enrichment Analyses

Gene set enrichment analysis (GSEA) was used to identify the function of target proteins that interact with the active ingredients of LMR in gene function and signaling pathway [21]. The KEGG [22] pathways significantly associated with the predicted genes were identified. We analyzed the Gene Ontology (GO) function and KEGG pathway enrichment of proteins (25 interacted target proteins) involved in the PPI network. The target proteins involved in the molecular function (MF), cellular components (CC), biological process (BP), and the KEGG pathways were also described. The , calculated by the Benjamini–Hochberg method [22], was considered as significant.

2.9. Statistical Analysis

Data are presented as a . Data analysis was made by SPSS (Statistical package for Social Science, Version 22.0) using one-way analysis of variance (ANOVA) followed by Tukey’s multiple posthoc tests. At , the values were found to be significantly different.

3. Results

3.1. Primary Phytochemical Screening and GC-MS Analysis

Figure 1 provides an array of compounds (Table 1) obtained from the GC-MS analysis of HKEx. The constituents were identified by comparing their mass spectrum with those in the machine library and with genuine compounds. Oleanolic acid, stigmasterol, NZ-14 N-Hexadecanoic acid (37.15%), 7 α, 28-olean diol, 9-octadecenoic acid, octadecanoic acid (12.56%), γ-sitostenone (5.88%), 1,2,3-propanetriyl ester, (E, E, E)-(18.87%), and γ-Sitosterol (4.13%) were the main components. Among these compounds, three most abundant compounds (oleanolic acid, stigmasterol, and 7 α, 28-olean diol in Figure 2) were also identified by other authors via vacuum liquid chromatography (VLC) analysis [17].

(a)

(b)

(c)

3.2. Changes of Blood and Biochemical Parameters

Figure 3(a) demonstrates the oral glucose tolerance level of the experimental animals over the third week of intervention. The glucose tolerance ability of the lower dose (HKEx50) was significantly () higher than the other groups including the diabetic control (DBC) group. The amount of liver glycogen present in experimental groups is shown in Figure 3(b). Liver glycogen in the HKEx-treated groups was significantly higher than in the DBC group, and even HKEx50 seemed to be more effective. Serum creatinine in the HKEx100 and HKEx200 was found to be significantly lower than in the DBC group as shown in Figure 3(c), but serum uric acid levels were increased in all the treatment groups (Figure 3(c)). The HKEx treatments resulted in a significant increase in insulin production compared to DBC groups (Figure 3(d)).

(a)

(b)

(c)

(d)

Changes in serum alanine aminotransferase (ALT) and aspartate aminotransferase (AST) enzymes were summarized in Figure 4(a). All doses of HKEx were found to significantly reduce serum ALT levels while only HKEx50 and HKEx200 were able to decrease AST level significantly compared to DBC. Figure 4(b) presents lactate dehydrogenase (LDH) and the heart-specific isozyme creatinine kinase (CK-MB) values. The serum LDH and CK-MB concentrations of the treatment groups HKEx50 and HKEx100 were significantly () reduced compared to diabetic control animal. Figure 4(c) displays the data for the serum lipid profile. The HDL (high-density lipoprotein) and LDL (low-density lipoprotein) of the treatment groups were significantly () normalized by HKEx100. Triglyceride concentration was also found to be lowered in the treatment group (HKEx100). Figure 4(d) illustrates the mRNA expression of genes encoded for biological enzymes SOD1 and CAT. In RT-qPCR, the mRNA expression of superoxide dismutase 1 (SOD1) was increased with all doses of HKEx. HKEx50 had a substantial () elevation in the SOD1 mRNA expression among the treatment groups. HKEx100 was only able to upregulate the mRNA expression for CAT.

(a)

(b)

(c)

(d)

3.3. Analysis of Interactions with Active Ingredients and Target Proteins

Antidiabetic mechanism of HKEx was evaluated through interactions of the target proteins with the functioning ingredients. We found that oleanolic acid (NZ-15) and stigmasterol (NZ-14) are significantly (PPI enrichment value is 0.0001) interacted with 27 target proteins (Figure 5). Oleanolic acid interacted with AKR1B10, CASP3, CASP8, CYP1A2, CYP1A2, HMGB1, NAMPT, NFE2L2, NQO1, PPARA, PTGIR, TOP1, TOP2A, UGT2B10, and UGT2B11 and stigmasterol interacted with ABCA1, ABCG5, ABCG8, CTSE, HMGCR, IL10, CXCL8, NR1H2, NR1H3, SLCO1B1, SREBF2, and TNF. Interestingly, we found that oleanolic acid and stigmasterol also interact with each other in the compound targets interactions. Compounds and target protein interactions are displayed in Figure 5.

3.4. Construction and Analysis of Target Protein PPI Network

PPI network plays substantial roles in the molecular processes, and abnormal PPI is the basis of many pathological conditions. Using the STRING database and Cytoscape plugin cytoHubba, all target proteins were mapped into the PPI network. Interestingly, we found that 25 target proteins are involved in PPI which have 77 edges, and average node degree 5.7 with the PPI enrichment value was less than . This PPI network indicates that the greater the node degree, the better the relationship between the proteins which signify their role in the entire interaction network. Only CTSE and FTGIR proteins are not included in PPI. We got only one subnetwork in PPI, which included 25 target proteins (listed with degree of interaction in supplementary Table S1). Twenty five target proteins with top degree of interactions with other proteins illustrated in Figure 6 and in supplementary Table S1. ABCA1, TNF, CXCL8, CASP3, HMGCR, PPARA, SREBF2, ABCG8, ABCG8, ABCG5, and IL10 are the top ten hub proteins with maximum degree of interaction (Figure 6 and Supplementary Table S1).

3.5. GO Analysis of Interacted Target Proteins

GO enrichment analysis of interacted target proteins (total 25) that interact with compounds was performed by GSEA. The top 20 significantly enriched terms in biological process (BP), 20 molecular functions (MF), and total 14 cellular components (CC) were selected, according to . Significant BP was listed in supplementary Table S2, and top 20 BP are presented in Figure 7(a). Fourteen top identified CCs are illustrated in Figure 7(b) and supplementary Table S4. In molecular function (MF) category, top twenty molecular functions such as cytokine activity, cytokine receptor binding, and tumor necrosis factor receptor superfamily binding, lipid transporter activity, lipid transfer activity, and some other receptor activities are presented in Figure 7(c) and supplementary Table S3.

(a)

(b)

(c)

3.6. Target-Protein Set Enrichment Analysis of KEGG Pathways

Link between the target proteins and the routes was further evaluated through identification of 31 KEGG pathways which were significantly associated with the target proteins (Figure 8 and Supplementary Table S5). These pathways were primarily involved in metabolism, immune regulation, neurological regulation, cellular development (apoptosis and p53 signaling pathway), cellular signaling, and PPAR signaling pathway.

3.7. Pathway-Associated Target Proteins Are Interacting with SOD1 and CAT

We found that 15 target proteins (TNF, CXCL8, IL10, PPARA, CASP3, CASP8, AKR1B10, ABCA1, ABCG5, ABCG8, NR1H3, CYP3A4, CYP1A2, UGT2B11, and UGT2B10) are associated with the enrichment of pathways (Figure 9). The mRNAs of antioxidative enzymes SOD1 and CAT were found to increase multifold in animal intervention study. PPI interaction of SOD1, CAT, and pathway-associated target proteins was further investigated by using STRING tool and interestingly, 17 proteins were found to be associated with a single network with the PPI enrichment value 6.6-16 (Figure 9). CAT interaction degree was 10, and SOD1 interaction degree was 5.

4. Discussion

Plant-based therapies have long been used to treat T2DM due to their efficacy, affordability, and lack of adverse effects. This study [23] demonstrates the antidiabetic mechanism of a common antihyperglycemic plant combining both animal models and bioinformatic-based network pharmacology. Type 2 diabetes mellitus is a metabolic derangement with abnormal reaction to insulin and β-cell dysfunction of the pancreas leading to unexplained weight loss, energy balance changes, retinopathy, neuropathy, nephropathy, and vascular complications. This study evaluated the effects of HKEx on the involvement of SOD1 and CAT antioxidative enzymes and associated genes via a network pharmacological evaluation.

In this research, streptozotocin was chosen to induce diabetes, because its alkylating properties directly target pancreatic β-cells to damage [21]. Because of the structural similarity with glucose, selectivity for β-cells is correlated with preferential accumulation of streptozotocin in β-cells after entry through GLUT-2 glucose transporter receptors [24]. As part of drug discovery, plants have always been good source while many of the currently available drugs have been directly or indirectly derived from them. Antidiabetic action has been found in a wide range of plant-derived active principles. In that light, HKEx containing number of phytoconstituents could attribute to the decrease of blood glucose levels in animal model. In particular, the glucose tolerance capacity of the extract was noted at the lower dose of HKEx50 mg/kg during the third week of intervention. Moreover, liver glycogen levels were found substantially higher in the treatment groups than in the DBC group. The liver is essential for maintaining appropriate glucose homeostasis, as it delivers glucose during fasting and stores glucose afterward. However, these hepatic processes are dysregulated in diabetes mellitus, and this imbalance contributes to hyperglycemia in the fasted and postprandial states which is consistent with the decreased glycogen level in the liver of the DBC group [25]. The significant decrease of creatinine levels is an indication of protective effect of HKEx in diabetic nephropathy. Stimulation of insulin secretion by plant extract is documented by many researchers [26] who have noticed the potentiation of glucose-induced insulin release from isolate mouse islets. Our results are consistent to increase the insulin secretion by two of the doses.

Insulin resistance, on the other hand, results in hyperinsulinemia, which is responsible for dyslipidemia, increased sensitivity of LDL to oxidation, and decreased tolerance to glucose [27]. In several recent studies, serum ALT and AST have been found to be increased due to liver damage, and serum LDH and CK-MB have been increased due to cardiac damage in diabetic conditions, called diabetic cardiomyopathy [26, 28, 29]. Hyperglycemia is also associated with insulin resistance in T2DM and is a highly independent cardiovascular disease indicator. In treatment groups, the decrease in the CK-MB cardiac marker suggests that the extract could be showing a successful activity to boost cardiac conditions [30]. The information suggests that HKEx has a protective effect on diabetes-related heart damage. Serum ALT, AST, LDH, and CK-MB elevation in DBC appear to be due to STZ-induced cellular damage while ALT and AST were found to be substantially recovered by HKEx treatments implying the increased activity of liver and liver aminotransferase. Lipid homeostasis in this research is mainly achieved through the increase of HDL and decrease of LDL. This effect of HKEx might be attributable to the reduced synthesis of fatty acids and other cholesterol derivatives [31].

The major phytochemical constituents in addition to the three target compounds oleanolic acid, 7 α, 28-olean diol, and stigmasterol for network pharmacology, some fatty acids appear to be major more prevalent compounds in GC-MS analyses. The role of saturated fatty acids in the control of anti-inflammatory pathways has been justified by previous investigations. N-Hexadecanoic acid is literally displayed as antioxidant, anti-inflammatory, and hypocholesterolemic. The enzyme kinetic analysis showed that n-hexadecanoic acid acted as a competitive phospholipase A2 inhibitor and controlled inflammatory cytokine synthesis by controlling the release of arachidonic acid [32], octadecanoic acid and 9-octadecenoic acid, 1, 2, 3-propanetriyl ester, (E, E, E)-exercise antioxidants, immunomodulators, and immunostimulants [33]. Again, by growing the activity of Cu/Zn SOD and CAT in cultured cortical neurons, octadecanoic acid (stearic acid) exhibits defense against oxidative stress by substantially decreasing H2O2 treatment by 24 h and thus exerting neuroprotective effects [34]. The levels of glutathione, SOD1, and CAT are significantly increased by phytosterol and stigmasterol [35] which is consistent with our observation.

Antioxidant effects are generally known to be the most effective way to protect the number of diseases including diabetes. This plant was evident in promising in vitro antioxidant potentials that were confirmed with certain antioxidant enzymes’ gene expression profile. Interestingly, however, the administration of HKEx to the STZ-induced animals upregulated the expression profile of these genes, implying that an antioxidant-dependent event mediated the pancreatic protective effects of HKEx. In fact, the twofold increase in the SOD1 expression observed was remarkable because it could imply a very rapid HKEx-induced antioxidative response. Similar findings have been recorded with other extracts attributed to the potent phytochemicals of antioxidants, as observed in HKEx.

The in vitro and in vivo findings have been verified using network pharmacological study which is currently facing a great development opportunities and challenges in terms of theoretical analyses, algorithm development, and applications in the Big Data era, but the challenge of merging enormous clinical and experimental data with scientific verification to expose the regulatory mechanisms of network pharmacology in order to conduct research more effectively has become a major concern for researchers. Therefore, data validation before clinical application is an inevitable step to be followed to avoid ambiguous use of network pharmacology [36]. Therefore, an experimental validation needs to be warranted before applying these findings into clinical application. PPI network plays substantial roles in molecular processes, and abnormal PPI is the basis of many pathological conditions [37]. From the PPI of targeted protein, we identified hub nodes which are markedly associated with diabetes. For example, the 10-degree interaction protein ABCA1 plays an important role in controlling lipid metabolism, and a defect in this gene interferes with the lipid transport of HDL cholesterol associated with the production of T2DM [38]. Peripheral insulin resistance and elevated plasma glucose and insulin levels are correlated with TNF-alpha concentration prior to the onset of T2DM [39]. In male mice with diabetes, the CXCL8 antagonist is evident to enhance diabetic nephropathy and attenuate high glucose-induced mesangial injury [40]. Caspase-3 promotes diabetic kidney disease [41]. It was stated that the HMGCR gene in population studies was associated with higher risk of T2DM [42]. It is proposed that interleukin-10 (IL10), an anti-inflammatory cytokine, has a protective role in T2D [43]. Our results shows that the immunological target proteins including IL10, CXCL8, and TNF are located in the PPI networks with top degree of interaction, indicating that this PPI network is associated with immunological activities. Human chemokines have been shown to be associated with or involved in the pathogenesis of type 1 diabetes [44].

In addition, we identified GO pathway enrichment analysis of target protein. The GO analysis indicating that the target proteins may bind with plasma membrane, chromosome, chromatin, regulatory region nucleic acid, or/and cellular receptor of cells for mediating the process of metabolic, immunological, signaling, and/or other activities, so as to exert signaling and antidiabetic potential of compounds. Pathway analysis revealed some pathways which are also associated with diabetes. Retinol metabolism is the most significantly () enriched pathway. Retinoids and retinoid-related proteins involved with signaling molecules linking obesity with the production of type II diabetes and in pancreatic β-cell biology/insulin secretion [45]. Linoleic acid metabolism is metabolism is another enriched pathway, and it was indicated that the prevention of T2DM should include increased linoleic acid consumption [46]. Metabolic pathway research showed that possible metabolic pathways are primarily involved in the pathway of interconversion of pentose and glucoronate and the metabolism of starch and sucrose, which changed specifically in the T2D environment [47]. PPAR signaling pathway and adipocytokine signaling pathway are also dysregulated in T2DM [21, 48, 49]. Moreover, we found that CAT and SOD1 are controlling the PPI network of target proteins which are significantly associated with enriched pathways, suggesting that these two genes are responsible for controlling this signaling networks and enriched pathways. Altogether, it indicates that the target proteins are associated with the alteration or dysregulation of diabetic pathways and the phytoconstituents of L. macrophylla are highly prospective to be processed further for befitted uses in diabetes.

KEGG pathways are primarily involved in metabolism (retinol metabolism, cytochrome P450 xenobiotic metabolism, cytochrome metabolism of drugs, linoleic acid metabolism, drug metabolism, other enzymes, biosynthesis of steroid hormones, ascorbate and aldarate metabolism, interconversion of pentose and glucuronate, porphyrin and chlorophyll metabolism, and starch and sucrose metabolism), immune regulation (cytokine-cytokine receptor interaction, T cell receptor signaling pathway, natural killer cell-mediated cytotoxicity, adipocytokine signaling pathway, Toll-like receptor signaling pathway, allograft rejection, and systemic lupus erythematosus), neurological regulation (amyotrophic lateral sclerosis), cellular development (apoptosis and p53 signaling pathway), cellular signaling (NOD-like receptor signaling pathway, RIG-I-like receptor signaling pathway, and PPAR signaling pathway), molecular transportation (ABC transporters), and retinol metabolism linoleic acid consumption, interconversion of pentose and glucuronate metabolism [50] and PPAR signaling pathway [51–53]. Our anticipated SOD1 and CAT genes are very likely to influence the networks of pathways because all 17 proteins for PPI interaction of SOD1, CAT, and pathway-associated target proteins used in STRING tool were found to be associated with a single network with the PPI enrichment value 6.6-16. Importantly, CAT interaction degree 10 and SOD1 interaction degree 5 in the PPI network suggest that these two genes are responsible for controlling this signaling networks and enriched pathways. Based on the compound–protein target relationships, it is evident that compounds are acting on the multitarget genes and proteins which are substantial contributions of HKEx on biological and physiological activities to manage diabetes.

5. Conclusions

This research has meaningfully justified the use of L. macrophylla in T2DM which is a multigene-metabolic complication pivotalizing the oxidative stress. Network pharmacological tools successfully bridged the role of bioactive compounds from L. macrophylla in controlling different metabolic and molecular pathways including the upregulation of most prominent biological antioxidative enzymes SOD1 and CAT. This worthy piece of work can lead to unveil the prospects of this plant material through justifying other formulatory necessities.

Data Availability

Data used in this sutdy are included in the article and rest of them are reposited as Supplementary Materials.

Ethical Approval

The experimental protocol was developed in compliance with ethical guidelines and approved by the Ethical Committee of the University of Chittagong.

Conflicts of Interest

The authors declare no conflict of interest.

Authors’ Contributions

Md. Atiar Rahman (MAR) made the research question, stated the problem, and hypothesized the work plan. He made the research design, carried out the bench work, and continued data acquisition with the aid of Md Nazim Uddin (MNU). MNU particularly focused on the network pharmacological part of the work. Nouf Abubakr Babteen (NAB) and Afnan M. Alnajeebi (AMA) assisted in manuscript reviewing and fund acquisition. Zainul Amiruddin Zakaria (ZAZ) accomplished the network pharmacological work with MNU and provided the grant support for the study. MAR finally revised the manuscript. SMA assisted in research funding. All authors have gone through the manuscript and agreed to submit.

Acknowledgments

This work was supported by the Ministry of Higher Education, Malaysia awarded the Fundamental Research Grant Scheme (FRGS, Reference Number: FRGS/1/2018/SKK10/UPM/02/4). We appreciate and thank Taif University for the financial support for Taif University Researchers Supporting Project (TURSP-2020/105), Taif University, Taif, Saudi Arabia.

Supplementary Materials

The following materials are available in the journal’s online system: table S1: list of target proteins with degree of interactions in PPI network. Table S2: list of biological process. Table S3: list of molecular functions. Table S4: list of cellular components. Table S5: list of KEGG pathways. (Supplementary Materials)

References

World Health Organization, Global status report on non-communicable diseases 2010, WHO, Geneva, 2011.
American Diabetes Association, “Diagnosis and classification of diabetes mellitus,” Diabetes Care, vol. 33, Supplement_1, pp. S62–S69, 2010.
View at: Publisher Site | Google Scholar
R. H. Eckel, “Preventive cardiology by lifestyle intervention: opportunity and/or challenge? Presidential address at the 2005 American Heart Association Scientific Sessions,” Circulation, vol. 113, pp. 2657–2661, 2016.
View at: Google Scholar
S. H. Pooya, M. D. Jalali, A. D. Jazayery, A. Saedisomeolia, M. R. Eshraghian, and F. Toorang, “The efficacy of omega-3 fatty acid supplementation on plasma homocysteine and malondialdehyde levels of type 2 diabetic patients,” Nutrition Metabolism & Cardiovascular Disease, vol. 20, no. 5, pp. 326–331, 2010.
View at: Publisher Site | Google Scholar
S. Li and B. Zhang, “Traditional Chinese medicine network pharmacology: theory, methodology and application,” Chinical Journal of Natural Medicine, vol. 11, no. 2, pp. 110–120, 2013.
View at: Publisher Site | Google Scholar
L. Xu, Y. Li, Y. Dai, and J. Peng, “Natural products for the treatment of type 2 diabetes mellitus: pharmacology and mechanisms,” Pharmacological Research, vol. S1043-6618, p. 31519, 2018.
View at: Google Scholar
Q. Liu, S. Li, H. Quan, and J. Li, “Vitamin B12 status in metformin treated patients: systematic review,” PLoS One, vol. 9, no. 6, article e100379, 2014.
View at: Publisher Site | Google Scholar
M. K. Moon, K. Y. Hur, S. H. Ko et al., “Combination therapy of oral hypoglycemic agents in patients with type 2 diabetes mellitus,” The Korean Journal of Internal Medicine, vol. 32, no. 6, pp. 974–983, 2017.
View at: Publisher Site | Google Scholar
S. Akhter, M. A. Rahman, J. Aklima, M. R. Hasan, and J. M. K. Hasan Chowdhury, “Antioxidative Role of Hatikana (Leea macrophylla Roxb.) Partially Improves the Hepatic Damage Induced by CCl4 in Wistar Albino Rats,” BioMed Research International, vol. 2015, 356712 pages, 2015.
View at: Publisher Site | Google Scholar
K. H. Hossain, M. A. Rahman, M. Taher et al., “Hot methanol extract of _Leea macrophylla_ (Roxb.) manages chemical-induced inflammation in rodent model,” Journal of King Saud University-Science, vol. 32, no. 6, pp. 2892–2899, 2020.
View at: Publisher Site | Google Scholar
M. A. Rahman, J. M. K. H. Chowdhury, J. Aklima, and M. A. Azadi, “Leea macrophylla Roxb. leaf extract potentially helps normalize islet of β-cells damaged in STZ-induced albino rats,” Food Science and Nutrition, vol. 6, no. 4, pp. 943–952, 2018.
View at: Publisher Site | Google Scholar
A. N. Nizami, M. A. Rahman, N. U. Ahmed, and M. S. Islam, “Whole Leea macrophylla ethanolic extract normalizes kidney deposits and recovers renal impairments in an ethylene glycol -induced urolithiasis model of rats,” Asian Pacific Journal of Tropical Medicine, vol. 5, no. 7, pp. 533–538, 2012.
View at: Publisher Site | Google Scholar
R. S. Singh and A. N. Singh, “On the identity and economic-medical uses of Hastikarnapalsa (Leea macrophylla Roxb., family Ampelidaceae) as evinced in the ancient texts and traditions,” Indian Journal of History Science, vol. 16, pp. 219–222, 1981.
View at: Google Scholar
M. B. Islam, M. M. H. Sarkar, M. Z. Shafique, M. A. Jalil, M. Z. Haque, and R. Amin, “Phytochemical screening and anti-microbial activity studies on Leea macrophylla seed extracts,” Journal of Scientific Research, vol. 5, no. 2, pp. 399–405, 2013.
View at: Publisher Site | Google Scholar
S. Lo, J. C. Russell, and A. W. Taylor, “Determination of glycogen in small tissue samples,” Journal of Applied Physiology, vol. 28, no. 2, pp. 234–236, 1970.
View at: Publisher Site | Google Scholar
J. Mawa, M. A. Rahman, M. A. Hashem, and M. J. Hosen, “Leea macrophylla root extract upregulates the mRNA expression for antioxidative enzymes and repairs the necrosis of pancreatic β-cell and kidney tissues in fructose-fed Type 2 diabetic rats,” Biomedicine and Pharmacotherapy, vol. 110, pp. 74–84, 2019.
View at: Publisher Site | Google Scholar
Z. A. Mahmud, S. C. Bachar, C. M. Hasan, T. B. Emran, N. Qais, and M. M. N. Uddin, “Phytochemical investigations and antioxidant potential of roots of Leea macrophylla (Roxb.),” BMC Research Notes, vol. 10, no. 1, p. 245, 2017.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. Santos, C. von Mering, J. L. Jensen, P. Bork, and M. Kuhn, “STITCH 5: augmenting protein-chemical interaction networks with tissue and affinity data,” Nucleic Acids Research, vol. 44, no. D1, pp. D380–D384, 2016.
View at: Publisher Site | Google Scholar
M. Kanehisa, M. Furumichi, M. Tanabe, Y. Sato, and K. Morishima, “KEGG: new perspectives on genomes, pathways, diseases and drugs,” Nucleic Acids Research, vol. 45, no. D1, pp. D353–D361, 2017.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. L. Gable, D. Lyon et al., “STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets,” Nucleic Acids Research, vol. 47, no. D1, pp. D607–D613, 2019.
View at: Publisher Site | Google Scholar
C. H. Chin, S. H. Chen, H. H. Wu, C. W. Ho, M. T. Ko, and C. Y. Lin, “CytoHubba: identifying hub objects and sub-networks from complex interactome,” BMC Systems Biology, vol. 8, no. S4, p. S11, 2014.
View at: Publisher Site | Google Scholar
Y. Benjamini and Y. Hochberg, “Controlling the false discovery rate: a practical and powerful approach to multiple testing,” Journal of the Royal Statistical Society Series B, Statistical Methodology, vol. 57, no. 1, pp. 289–300, 1995.
View at: Publisher Site | Google Scholar
H. Li, H. Li, J. Ran, and X. Mao, “Dysregulated pathways in type 2 diabetes mellitus,” Biomedical Research-India, vol. 28, pp. 1426–1443, 2017.
View at: Google Scholar
M. L. Graham, J. L. Janecek, J. A. Kittredge, B. J. Hering, and H. J. Schuurman, “The Streptozotocin-induced diabetic nude mouse model: differences between animals from different sources,” Comparative Medicine, vol. 61, no. 4, pp. 356–360, 2011.
View at: Google Scholar
P. Ansari, P. R. Flatt, P. Harriott, and Y. H. A. Abdel-Wahab, “Evaluation of the antidiabetic and insulin releasing. effects of A. squamosa, including isolation and characterization of active phytochemicals,” Plants (Basel), vol. 9, no. 10, article 1348, 2020.
View at: Google Scholar
D. K. Patel, S. K. Prasad, R. Kumar, and S. Hemalatha, “An overview on antidiabetic medicinal plants having insulin mimetic property,” Asian Pacific Journal of Tropical Biomedicine, vol. 2, no. 4, pp. 320–330, 2012.
View at: Publisher Site | Google Scholar
K. Sharabi, D. Clint, J. Tavares, A. K. Rines, and P. Puigserver, “Molecular pathophysiology of hepatic glucose production,” Molecular Aspects of Medicine, vol. 46, pp. 21–33, 2015.
View at: Publisher Site | Google Scholar
V. V. Kuzhandai, R. Ramesh, and A. R. Srinivasan, “Circulating microRNAs as biomarkers in health and disease,” Journal of Clinical Diagnostics and Research, vol. 6, pp. 1791–1795, 2012.
View at: Google Scholar
M. Sankaran and A. Vadivel, “Antioxidant and antidiabetic effect of hibiscus rosasinensis flower extract on streptozotocin induced experimental rats-a dose response study,” Notulae Scientia Biologia, vol. 3, no. 4, pp. 13–21, 2011.
View at: Publisher Site | Google Scholar
N. Maghamiour and N. Safaie, “High creatine kinase (CK)-MB and lactate dehydrogenase in the absence of myocardial injury or infarction: a case report,” Journal of Cardiovascular and Thoracic Research, vol. 6, no. 1, pp. 69-70, 2014.
View at: Publisher Site | Google Scholar
A. J. Afolayan and E. A. Adewusi, “A review of natural products with hepatoprotective activity,” Journal of Medicinal Plants Research, vol. 4, pp. 1318–1334, 2010.
View at: Google Scholar
B. B. Patel, M. Di Iorio, and L. E. Chalifour, “Metabolic response to chronic bisphenol A exposure in C57bl/6n mice,” Toxicology Report, vol. 7, pp. 522–532, 2014.
View at: Google Scholar
R. Ruiz, V. Jideonwo, M. Ahn et al., “Sterol regulatory element-binding protein-1 (SREBP-1) is required to regulate glycogen synthesis and gluconeogenic gene expression in mouse liver,” Journal of Biological Chemistry, vol. 289, no. 9, pp. 5510–5517, 2014.
View at: Publisher Site | Google Scholar
M. Abdollahi, A. B. Z. Zuki, Y. M. Goh, A. Rezaeizadeh, and M. M. Nooedin, “Effects of Momordica charantia on pancreatic histopathological changes associated with streptozotocin-induced diabetes in neonatal rats,” Histology and Histopathology, vol. 26, no. 1, pp. 13–21, 2011.
View at: Publisher Site | Google Scholar
P. H. Pan, S. Y. Lin, Y. C. Ou et al., “Stearic acid attenuates cholestasis-induced liver injury,” Biochemistry Biophysics Research Communication, vol. 391, no. 3, pp. 1537–1542, 2010.
View at: Publisher Site | Google Scholar
S. Li, “Network pharmacology evaluation method guidance-draft,” World Journal of Traditional Chinese Medicine, vol. 7, no. 1, pp. 148–154, 2021.
View at: Publisher Site | Google Scholar
Z. J. Wang, C. L. Liang, G. M. Li, C. Y. Yu, and M. Yin, “Stearic acid protects primary cultured cortical neurons against oxidative stress,” Acta Pharmacologica Sinica, vol. 28, no. 3, pp. 315–326, 2007.
View at: Publisher Site | Google Scholar
H. Ali, S. Dixit, D. Ali, S. M. Alqahtani, S. Alkahtani, and S. Alarifi, “Isolation and evaluation of anticancer efficacy of stigmasterol in a mouse model of DMBA-induced skin carcinoma,” Drug Design Development and Therapy, vol. 9, pp. 2793–2800, 2015.
View at: Google Scholar
G. Wilcox, “Insulin and insulin resistance,” Clinical Biochemistry Reviews, vol. 26, no. 2, pp. 19–39, 2005.
View at: Google Scholar
B. S. Haerian, M. S. Haerian, A. Roohi, and H. Mehrad-Majd, “ABCA1 genetic polymorphisms and type 2 diabetes mellitus and its complications,” Meta Gene, vol. 13, pp. 104–114, 2017.
View at: Publisher Site | Google Scholar
Y. Miyazaki, R. Pipek, L. J. Mandarino, and R. A. DeFronzo, “Tumor necrosis factor α and insulin resistance in obese type 2 diabetic patients,” International Journal of Obesity Related Metabolic Disorder, vol. 27, no. 1, pp. 88–94, 2003.
View at: Publisher Site | Google Scholar
S. Wen, Z. H. Wang, C. X. Zhang, Y. Yang, and Q. L. Fan, “Caspase-3 promotes diabetic kidney disease through gasdermin E-mediated progression to secondary necrosis during apoptosis,” Metabolic Syndrome and Obesity, vol. 13, pp. 313–323, 2020.
View at: Publisher Site | Google Scholar
D. I. Swerdlow, D. Preiss, K. B. Kuchenbaecker et al., “HMG-coenzyme A reductase inhibition, type 2 diabetes, and bodyweight: evidence from genetic analysis and randomised trials,” The Lancet, vol. 385, no. 9965, pp. 351–361, 2015.
View at: Publisher Site | Google Scholar
J. C. Barry, S. Shakibakho, C. Durrer et al., “Hyporesponsiveness to the anti-inflammatory action of interleukin-10 in type 2 diabetes,” Scientific Reports, vol. 6, no. 1, article 21244, 2016.
View at: Publisher Site | Google Scholar
S. A. Sarkar, C. E. Lee, F. Victorino et al., “Expression and regulation of chemokines in murine and human type 1 diabetes,” Diabetes, vol. 61, no. 2, pp. 436–446, 2012.
View at: Publisher Site | Google Scholar
P. J. Brun, K. J. Z. Yang, S. A. Lee, J. J. Yuen, and W. S. Blaner, “Retinoids: potent regulators of metabolism,” BioFactors, vol. 39, no. 2, pp. 151–163, 2013.
View at: Publisher Site | Google Scholar
G. Henderson, C. Crofts, and G. Schofield, “Linoleic acid and diabetes prevention,” The Lancet Diabetes & Endocrinol, vol. 6, no. 1, pp. 12-13, 2018.
View at: Publisher Site | Google Scholar
H. Sun, S. Zhang, A. Zhang et al., “Metabolomic analysis of diet-induced type 2 diabetes using UPLC/MS integrated with pattern recognition approach,” PLoS One, vol. 9, no. 3, article e93384, 2014.
View at: Publisher Site | Google Scholar
F. H. Wong, C. Y. Huang, L. J. Su et al., “Combination of microarray profiling and protein-protein interaction databases delineates the minimal discriminators as a metastasis network for esophageal squamous cell carcinoma,” International Journal of Oncology, vol. 34, pp. 117–128, 2009.
View at: Publisher Site | Google Scholar
S. Cui, Y. Zhu, J. Du et al., “CXCL8 antagonist improves diabetic nephropathy in male mice with diabetes and attenuates high glucose–induced mesangial injury,” Endocrinology, vol. 158, no. 6, pp. 1671–1684, 2017.
View at: Publisher Site | Google Scholar
M. Shen, M. Xu, F. Zhong et al., “A multi-omics study revealing the metabolic effects of estrogen in liver cancer cells HepG2,” Cells, vol. 10, no. 2, p. 455, 2021.
View at: Publisher Site | Google Scholar
P. Subramanian, V. K. Tamayo, S. Mootha et al., “Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles,” Proceedings National Academy Sciences, vol. 102, no. 43, pp. 15545–15550, 2005.
View at: Publisher Site | Google Scholar
P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Md. Atiar Rahman et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

505

Downloads

1056

Citations