Skip to main content
  • Research article
  • Open access
  • Published:

A systematic analysis of natural α-glucosidase inhibitors from flavonoids of Radix scutellariae using ultrafiltration UPLC-TripleTOF-MS/MS and network pharmacology

Abstract

Background

Flavonoids from plant medicines are supposed to be viable alternatives for the treatment of type 2 diabetes (T2D) as less toxicity and side effects. Radix scutellariae (RS) is a widely used traditional medicine in Asia. It has shown great potential in the research of T2D. However, the pharmacological actions remain obscured due to the complex chemical nature of plant medicines.

Methods

In the present study, a systematic method combining ultrafiltration UPLC-TripleTOF-MS/MS and network pharmacology was developed to screen α-glucosidase inhibitors from flavonoids of RS, and explore the underlying mechanism for the treatment of T2D.

Results

The n-butanol part of ethanol extract from RS showed a strong α-glucosidase inhibition activity (90.55%, IC50 0.551 mg/mL) against positive control acarbose (90.59%, IC50 1.079 mg/mL). A total of 32 kinds of flavonoids were identified from the extract, and their ESI-MS/MS behaviors were elucidated. Thirteen compounds were screened as α-glucosidase inhibitors, including viscidulin III, 2′,3,5,6′,7-pentahydroxyflavanone, and so on. A compound-target-pathway (CTP) network was constructed by integrating these α-glucosidase inhibitors, target proteins, and related pathways. This network exhibited an uneven distribution and approximate scale-free property. Chrysin (k = 87), 5,8,2′-trihydroxy-7-methoxyflavone (k = 21) and wogonin (k = 20) were selected as the main active constituents with much higher degree values. A protein-protein interaction (PPI) weighted network was built for target proteins of these α-glucosidase inhibitors and drug targets of T2D. PPARG (Cd = 0.165, Cb = 0.232, Cc = 0.401), ACACB (Cd = 0.155, Cb = 0.184, Cc = 0.318), NFKB1 (Cd = 0.233, Cb = 0.161, Cc = 0.431), and PGH2 (Cd = 0.194, Cb = 0.157, Cc = 0.427) exhibited as key targets with the highest scores of centrality indices. Furthermore, a core subnetwork was extracted from the CTP and PPI weighted network. Type II diabetes mellitus (hsa04930) and PPAR signaling pathway (hsa03320) were confirmed as the critical pathways.

Conclusions

These results improved current understanding of natural flavonoids on the treatment of T2D. The combination of ultrafiltration UPLC-TripleTOF-MS/MS and network pharmacology provides a novel strategy for the research of plant medicines and complex diseases.

Peer Review reports

Background

Type 2 diabetes (T2D) is one of the most serious chronic metabolic disorders characterized by persistent hyperglycemia. It accounts for more than 90% of all diabetes [1], and is directly linked to pathogenic consequences in the eyes, kidney, and cardiovascular diseases [2]. Natural products with less side effects have been used to treat diabetes for thousands of years [3]. Previously, various natural products were found to exhibit anti-diabetic effects, such as herbal formulas, fruits, vegetables, spices, or natural beverages [4]. These remedies are more accessible and affordable than modern pharmaceuticals [5].

Radix scutellariae (RS) is the dried root of Scutellaria baicalensis [6]. It is widely used as herbal medicine in Asia for thousands of years [7]. This medicine has various therapeutic functions, including antitumor, cardiovascular, neuroprotective and anti-inflammatory effects [6, 8, 9]. A growing body of exciting evidences also indicated an antidiabetic effect of RS. For instance, water extract of RS showed a reduction of body weight and blood triglyceride in type 2 diabetic db/db mice [10]. Methanol extract of RS had strong α-glucosidase inhibition [11]. Another traditional medicine coptis together with RS demonstrated potent anti-hyperglycemic effect on diabetic rats [12]. In addition, ethanolic extract of RS was found to enhance the antidiabetic effect of metformin, and increase pancreatic insulin content as well as improving the lipid profile in diabetic Wistar rats [13]. These reports suggested great potential of RS in the drug discovery of T2D.

The α-glucosidase is an exo-type carbohydrate enzyme that catalyzes the liberation of α-glucose from the non-reducing end of the carbohydrates. It locates in the brush border surface membrane of the small intestinal cells. This enzyme accelerates glucose reabsorption in the intestine [14]. Inhibition of α-glucosidase could delay the digestion, absorption of carbohydrates, and suppress postprandial hyperglycemia [15]. Natural α-glucosidase inhibitors have presented viable alternatives to the treatment of T2D as fewer toxicity and adverse effects [16]. More than one hundred herbal medicines have exhibited great potency in α-glucosidase inhibition and equivalent efficacies to synthetic drugs in managing diabetes [17].

Numerous compounds of natural origin are considered as models for drug discovery of T2D, such as flavonoids, polyphenols, terpenoids, alkaloids, saponins, quinones [18]. Flavonoids are a group of natural polyphenolic derivatives that widely exist in traditional medicines [19]. In recent years, considerable portions of natural flavonoids displayed anti-diabetic effects, including quercetin, rutin, naringin, baicalein [20,21,22]. Many flavonoids were also found in RS, such as baicalein, baicalin, wogonin, wogonoside, and so on [23]. A part of these compounds were reported to exhibit α-glucosidase inhibition [22]. Therefore, RS is considered as a source of natural α-glucosidase inhibitors.

The identification of the pharmacological profile of natural products is always a challenging task. Ultrafiltration method has attracted much attention in the screening and analysis of bioactive compounds from botanical extracts [24]. It has the advantage of high-speed and high-reliability, which facilitates the separation of ligand-receptor complexes for unbound compounds [25, 26]. For instance, Zhang et al. established an ultrafiltration LC-MS method for screening and characterizing thrombin inhibitors from Rhizoma. Wang et al. applied ultrafiltration LC-MS combined with reverse phase-medium pressure liquid chromatography for screening and isolation potential α-glucosidase inhibitors from RS [27]. However, these studies mainly focused on the screening of bioactive molecules. Further researches are urgently needed to elucidate the underlying mechanism.

Molecular mechanism of natural products is always difficult and confused as the complex chemical nature [28, 29]. Recently, network pharmacology approach, also known as system pharmacology, has emerged as a powerful tool to solve the problem [30, 31]. This methodology holds a significant potential for extracting biological information from large amounts of chemical data [32], and enables to predict the target profiles and pharmacological actions of herbal compounds [33, 34]. Chen et al. constructed a multi-parameter network model on the basis of three important parameters to tentatively explain the anti-fibrosis mechanism of herbal medicine Sophora flavescens [35]. Luo et al. used systems pharmacology strategies for anti-cancer drug discovery based on natural products [36]. Gogoi et al developed a network pharmacology-based virtual screening of natural products from Clerodendrum species for identification of novel anti-cancer therapeutics [37]. These studies demonstrated that network pharmacology approach had real potential in the mechanism research of natural products [38, 39].

In the present study, we developed a systematic method to screen α-glucosidase inhibitors from plant flavonoids and explore the underlying mechanism. Ultrafiltration UPLC-TripleTOF-MS/MS was used to identify flavonoids from the RS extract and screen potential α-glucosidase inhibitors. Network pharmacology was applied to investigate the interrelationships between these compounds, related target proteins and pathways. Several networks were constructed, and a series of topological characteristics were calculated to determine the main active constituents, key targets and critical pathways.

Methods

Materials and reagents

Crude Radix scutellariae was purchased from Baoji Medicinal Material Company (Shaanxi, China). The plant species was authenticated by Prof. Xiaomei Wang from Shaanxi Key Laboratory of Phytochemistry in Baoji University of Art and Sciences. Wogonin (HW158604), baicalin (HB158602), wogonoside (HW158601), oroxylin A (HB158728), chrysin (C110078), skullcapflavone II (HA062620), oroxylin A− 7-O-β-D-glucuronopyranoside (HO158605), baicalein-6-O-β-D-glucuronopyranoside (XB161661), and chrysin-7-O-β-D-glucuronopyranoside (HA061609) were obtained from Chenguang Biotech Co. Ltd. (Shaanxi, China) with a purity higher than 98%. The α-Glucosidase (from Saccharomyces cerevisiae) was purchased from Sigma-Aldrich (St. Louis, MO, USA). Acarbose and p-nitrophenyl-α-D-glucopyranoside (α-pNPG) were acquired from Aladdin Industrial Corporation (Shanghai, China). Methanol of HPLC grade was supplied by Merck (Darmstadt, Germany). Formic acid (HPLC grade) was purchased from Tedia (Fairfield, OH, USA). Ultrapure water was obtained using a Milli-Q purification system (Millipore Co., USA). All other reagents were analytically pure.

Standards and sample preparation

Reference substances were accurately weighed and dissolved in methanol (5 μg/mL). The solutions were stored at 4 °C until use. Dried Radix scutellariae powders (100 g) were passed through 100-mesh sieves, then orderly extracted with n-hexane, chloroform, 70% ethanol by heating reflux for 2–3 h, three times. Then, the 70% ethanol solution was leached with n-butanol (saturated by water). The solvents were removed by evaporation in vacuo, and the extracts were stored at − 20 °C until required, thawed at room temperature, dissolved in methanol (1 mg/mL). Finally, the solution was filtered with 0.22 μm Millipore filter membrane, and used directly for LC-MS.

α-Glucosidase inhibition assay

The α-glucosidase inhibitory activity was evaluated based on the slightly modified method of the literature [40]. The assay mixture (160 μL) contained 20 μL of phosphate buffer (0.1 M, pH 7.0) in 96-well plates, 20 μL of enzyme solution (0.1 U /mL α-glucosidase in phosphate buffer), 20 μL of sample in phosphate buffer with different concentrations to be mixed and incubated at 37 °C for 15 min. Then, the reaction was initiated by adding 20 μL of α-pNPG (2.5 mM in phosphate buffer). After 15 min at 37 °C, the reaction was stopped by adding 80 μL Na2CO3 (0.2 M) solution. Amount of released PNP (4-nitrophenol) was quantified by a microplate reader at the absorbance of 405 nm. The inhibitory rates (%) were calculated as follows:

$$ \mathrm{Inhibition}\%=1-\left[\left({\mathrm{A}}_{\mathrm{s}}-{\mathrm{A}}_{\mathrm{c}}\right)/{\mathrm{A}}_{\mathrm{b}}\right]\times 100 $$

Where the symbol ‘As’ is the absorbance of the test sample; ‘Ac’ is the absorbance of the sample contrast (without enzyme solution); and ‘Ab’ is the absorbance of the blank (without tested sample). All reactions were conducted in three replications and acarbose was used as positive control. The half maximal inhibitory concentration of the test sample (IC50) was calculated using the modified Karber’s method.

Screening of α-glucosidase inhibitors from RS

A 2 μL aliquot of n-butanol part of ethanol extract from RS (50 mg/mL) was incubated with 8 μL of α-glucosidase (100 μM, dissolved in 10 mM ammonium acetate buffer, pH 6.86) for 30 min at 37 °C. After incubation, each sample was filtered through a Vivaspin 2 concentrator (MWCO 10 kDa, Sartorius, Göttingen, Germany) at 10,000 g for 10 min. Then, the filter was washed three times with 200 μL ammonium acetate buffer (PH 6.86) to remove the unbound compounds. The bound ligands were released by adding 200 μL of methanol/water mixtures (1:1, v/v) adjusted with acetic acid to pH 3.30, followed by centrifugation at 10,000 g for 15 min. This procedure was repeated three times. A control experiment in which α-glucosidase omitted was also carried out before each screening experiment. The released ligands were then analyzed by LC-MS.

UPLC-TripleTOF-MS/MS analysis

Chromatographic separations were achieved on LC-20ADXR (Shimadzu, Tokyo, Japan) coupled with a Shim-pack XR-ODS column (100 mm × 2.0 mm, 2.2 μm, Shimadzu). The mobile phases consisted of eluent A (0.1% formic acid in water, v/v) and eluent B (0.1% formic acid in methanol, v/v). The gradient elution program was set as follows: 10 to 48% B from 0 to 8 min, holding for 6 min, 48 to 100% B from 14 to 20 min. After holding 100% B for next 5 min, the ratio was returned to its starting condition. The injection volume was 20 μL (20 μg/mL) at a flow rate of 0.3 mL/min. The column was maintained at 40 °C. MS analysis was performed on a TripleTOF 4600 mass analyzer (AB SCIEX, USA) equipped with electrospray ionization (ESI) source. The instrument was operated in positive ESI mode with a declustering potential voltage (DP) of 100 V and ionspray voltage of 5.5 KV. The nebulization temperature was 550 °C. GS1, GS2 and curtain gas were maintained at 55, 55 and 30 psi, respectively. Collision energy was 10 eV for MS and 50 eV for MS/MS. An automated calibration delivery system (CDS) was applied to regulate the MS and the MS/MS. The constituents were identified by the comparison with reference standards, the accurate molecular weights (with a mass tolerance of ±5 ppm), as well as the MS/MS fragment patterns. The operations, acquisition, and analysis of data were monitored by Analyst TF 1.7 (AB SCIEX, Concord, Canada) and PeakView 2.0 (AB SCIEX, Concord, Canada).

Collection of target proteins and pathways enrichment analysis

Target proteins of the α-glucosidase inhibitors from RS were collected using SuperPred (http://prediction.charite.de/) and DrugBank (https://www.drugbank.ca/). The target prediction is based on the similarity distribution among the targets’ ligands. The distributions are utilized for estimating individual thresholds and probabilities for a specific target. By means of these individual thresholds and probabilities, the input compound was screened against a database containing about 341,000 compounds, 1800 targets and 665,000 compound-target interactions [41]. Information of all the targets was uniformed by Uniprot (http://www.uniprot.org/). Pathway analysis was applied to these proteins by DAVID 6.8 (https://david.ncifcrf.gov/). The queried species was Homo sapiens. Raw P-values were adjusted using Benjamini & Hochberg procedure [42]. Pathways with adjusted P-values less than 0.05 were considered as significant.

Construction of networks for the α-glucosidase inhibitors from RS

A complex network analysis was performed on the collected data for further interpretation. First, a compound-target-pathway (CTP) network was constructed to screen main active components from RS. This network consisted of numerous nodes and edges. Nodes represented the α-glucosidase inhibitors from RS, corresponding target proteins and related pathways, respectively. Edges referred to interactions between them. If a protein was the hit target of particular inhibitor, or involved in any pathways, connections were made between these elements.

Subsequently, we collected therapeutic targets of T2D from TTD database (https://db.idrblab.org/ttd/) [43], and integrated with targets of the α-glucosidase inhibitors into a protein-protein interaction (PPI) weighted network. This network was designed to evaluate the closeness of interaction between RS and T2D. Interactions between the two groups of proteins were calculated by Search Tool for the Retrieval of Interacting Genes/Proteins (STRING, https://string-db.org/) [44]. STRING uses a scoring mechanism to give a comprehensive score to the results obtained by these different methods, including experimental data, data mined from PubMed abstract text, database data, and results predicted by bioinformatics methods. A weighted protein-protein interaction (PPI) network was constructed on the basis of these data. Nodes indicated the proteins, and that connections represented interactions between them with scores higher than 0.7.

Moreover, key nodes of the PPI network and their neighbor nodes were extracted, as well as the directly connected α-glucosidase inhibitors and related pathways. These elements were reconstructed as a core subnetwork to explore the underlying pharmaceutical mechanism of RS. Construction and visualization of all the networks were performed by Pajek ver. 2.00 (Batagelj and Mrvar, 2009).

Statistical and topological analysis of the network

To interpret the behavior of the α-glucosidase inhibitors from RS and T2D, several topological parameters of the network were analyzed (Table 1). The degree ki is the number of its connections attached to a given node i. The directly linked nodes are called neighbors of node i. Mean value of k of all nodes is defined as average degree k. Degree distribution is the proportion of randomly selected nodes with specific number of connections, and denoted as P(k). Average path length (L) refers to the mean distance between each pairs of nodes, which measures the overall navigability of a network. The diameter (D) is the maximum distance between any pair of nodes.

Table 1 Definitions of the topological parameters used in the network analysis

Centrality measures the relative influence of a node within the overall architecture of a network. In this study, three centrality metrics were comprehensively evaluated. Each of them focused on specific influence of a node on other nodes. Degree centrality (Cd) indicates the proportion of other nodes adjacent to a node, representing the immediate influence that the closest nodes produce on the corresponding vertex. Betweenness centrality (Cb) refers to the total number of shortest paths going through a node, which directly reflects the influence of a node has on the spread of information through the network. Closeness centrality (Cc) is the number of other nodes divided by the sum of distances between one node and the others, reflecting how close a node is to others. The statistical analysis was performed with MATLAB 2016a (The MathWorks Inc., Natick, MA, USA).

Results

Identification of flavonoids from Radix Scutellariae

The n-butanol part of ethanol extract from Radix Scutellariae were analyzed by UPLC-TripleTOF-MS/MS. A total of 32 kinds of flavonoids were identified within 5 ppm mass tolerance. Nine of them (compound 10, 11, 12, 14, 16, 24, 26, 28, 29) were confirmed by the reference standards, and the others by fragmentation analysis. The identification results were also compared with those from previous studies to ensure the accuracy. Detailed MS data of these compounds are listed in Table 2, and the MS/MS fragmentation patterns of typical flavonoids from RS are shown in Additional file 1.

Table 2 Characterization of flavonoids from Radix Scutellariae extracts by UPLC-TripleTOF-MS/MS

The identified flavonoids contained 20 aglycones (compound 1, 3, 5, 6, 8, 9, 17–21, 23–31) and 12 glycosides (compound 2, 4, 7, 10–16, 22, 32). Compound 1, 3, 8, 28 belonged to aglycones, substituted by several hydroxyl groups. Compound 28 showed [M + H]+ peak at m/z 255.0654 (Additional file 1a). This protonated molecular ion yielded the product ions at m/z 153.0192, 103.0547 in the MS/MS spectra. The two ions were attributed to the 1,3A+ and 1,3B+, indicating the occurrence of two OH groups in A-ring and none OH in B-ring. It was consistent with the report of chrysin by Luo et al [7]. Hence, compound 28 was tentatively identified as chrysin, which was further confirmed by the standard. Similarly, compound 1, 3, 8 were identified as 2′,3,5,6′,7-pentahydroxyflavanone, 2′,5,6′,7-tetrahydroxyflavane, and dihydrobaicalein, respectively [45, 46].

Compound 5, 6, 9, 17–21, 23–27, 29–31 belonged to the methoxylated flavonoid aglycones. They exhibited a characteristic ion (15*n Da) due to the loss of CH3 radicals. Protonated molecular ion of compound 26 was observed at m/z 285.0764 (Additional file 1b). The sole flavone aglycone easily gave a prominent ion [M + H-15]+ at m/z 270.2520, originated from losing one CH3 (15 Da). It also lost a CHO (29 Da) from C-ring, and produced the fragment at m/z 241.0494. In addition, the neutral loss of H2O (18 Da) from m/z 270.2520 produced the ions at m/z 252.0420 and 179.0462. Diagnostic fragment ions originated from Retro-Diels-Alder (RDA) reaction are often helpful to the structural determination of A- and B-ring substitution patterns [45]. Our data showed fragment ions 1,3A+ (m/z 168.0076) and 1,3B+ (m/z 103.0523), originated from the cleavage of the bond at position 1/3 of C-ring. It was also annotated by Ma et al. [47], and that compound 26 was finally identified as wogonin.

Compound 17 showed [M + H]+ peak at m/z 301.0703 (Additional file 1c). MS/MS spectra of this compound exhibited a methoxylated flavone characteristic loss of CH3 (15 Da), resulting in a product ion at m/z 286.0501. Besides, the parent ion m/z 301.0703 also yielded the ions at m/z 184.0006, 156.0054, 137.9944, and 119.0452. The product ion m/z 184.0006 was attributed to the 1,3A+, indicating that the substituent groups of two OH and an OCH3 were located in A-ring. The ion m/z 156.0054 was produced by the neutral loss of CO and H2O from the 1,3A+. This compound was finally identified as 4′-hydroxywogonin [46]. Compound 23 (Additional file 1d) exhibited a same characteristic loss of CH3 (15 Da) at m/z 286.0456. Moreover, other RDA fragments from the fragment ion, 1,3A+ at m/z 168.0050 and 1,4A+ at m/z 140.0096, could also be observed. It was identified as 5,8,2′-trihydroxy-7-methoxyflavone, the isomer of compound 17. Identification of other compounds in this group was also conducted by comparison with previous reports, including compound 5, 6, 9 [7, 45, 48], compound 18–21 [49, 50], compound 25 [51], compound 27 [52], compound 30–31 [8].

Compound 2, 4, 7, 10–16, 22, 32 belonged to flavonoid glycosides, glycosylated in different positions. Neutral loss of glucuronic acid (176 Da) or glucose (162 Da) is common in flavone glycoside, as the O-glucosylic bond is easily cleaved to generate aglycone. Reference standards of compound 10, 11, 12, 14, and 16 showed a fragment ion at [M + H-176]+ due to the loss of glucuronic acid. These compounds could be further distinguished by the fragment of residual aglycone. Compound 14 (Additional file 1e) and compound 16 (Additional file 1f) were a pair of isomers. They both went through the loss of glucuronic acid (176 Da), and produced the aglycone ions wogonin (8-OCH3) and oroxylin A (6-OCH3), respectively. Furthermore, they both exhibited a base peak [M + H-176-CH3] at m/z 270.0494. However, these two compounds could be distinguished by the relative abundances of the losses of CO and CHO from the parent ion [M + H-176-CH3]. In the MS/MS of compound 16, the relative abundance of ion [M + H-176-CH3-CO] at m/z 242.0598 was lower than [M + H-176-CH3-CHO] at m/z 241.0494, since the loss of CHO could produce a more stable p-quinoid skeleton. Nevertheless, it was contrary to compound 14. Finally, compound 14 was identified as oroxylin A-7-O-β-D-glucuronopyranoside, and compound 16 was wogonoside.

Compound 12 also exhibited the characteristic fragment ion [M + H-176] at m/z 255.0683. Moreover, other RDA fragments from the aglycone ion of chrysin (1,3A0+ at m/z 153.0186 and 1,3B0+ at m/z 103.0558) were observed. This compound was identified as chrysin-7-O-β-D-glucuronopyranoside, which was supported by the report by Luo et al [7]. Identification of other compounds in this group was according to previous studies, including compound 2 and 4 [7], compound 7 [53], compound 13 and 15 [8], compound 22 and 32 [54]. However, these results are only based on LC-MS/MS, which might be limited by various factors. More reference standards and analytical tools would be used to check the accuracy of identification in our next study.

Potential α-glucosidase inhibitors, target proteins and related pathways

An in vitro α-glucosidase inhibition assay was performed on the n-butanol part of ethanol extract from RS (Additional file 2). It showed higher α-glucosidase inhibitory activity (IC50 = 0.551 mg/mL) than the positive control (IC50 = 1.079 mg/mL). The inhibition rate reached 90.55% at the concentration of 2.34 mg/mL, whereas that of the positive control was 90.59% at 15 mg/mL. These data indicated that the crude extract of RS was much more potent than acarbose on α-glucosidase inhibition.

Ultrafiltration LC-MS/MS has been widely used to screen bioactive compounds from natural products [55]. In this study, the complexes of α-glucosidase and ligands from RS were retained by an ultrafiltration membrane, whereas the unbound, low molecular weight compounds were washed away from the chamber. Subsequently, the remainings were dissociated, and the ligands were identified by LC-MS/MS. Finally, a total of 13 peaks were detected, yet not presented in control group, indicating a specific binding to α-glucosidase. Chemical structures of the trapped ligands are shown in Fig. 1, including wogonin, chrysin, oroxylin A, 5,8,2′-trihydroxy-7-methoxyflavone, and so on. These compounds were considered as potential α-glucosidase inhibitors, and reorganized as a chemical ingredients database for the next analysis.

Fig. 1
figure 1

Chemical structures of the potential α-glucosidase inhibitors from flavonoids of Radix Scutellariae extract

Interactions between small molecules and proteins are often highly valued in biomedical and pharmaceutical sciences [56]. Numerous target proteins have been used for the treatment of T2D, such as insulin receptor, peroxisome proliferator activated receptor gamma (PPARG), and α-glucosidase [57]. In this study, target proteins of the 13 α-glucosidase inhibitors were collected using web tools. A total of 117 target proteins were obtained (Additional file 3). Some of them were therapeutic targets of T2D, such as bile acid receptor, histone deacetylase 1, prostaglandin G/H synthase 2, and so on [58]. It indicated that these α-glucosidase inhibitors might affect T2D through multi-targets.

A pathway contains a panel of cascade reactions among various biomolecules [59]. Pathway analysis demonstrated that the 117 targets of α-glucosidase inhibitors were involved in 86 pathways (Additional file 3), including metabolism of xenobiotics by cytochrome P450, steroid hormone biosynthesis, insulin resistance, retinol metabolism, and so on. These data were further analyzed by network pharmacology.

Compound-target-pathway (CTP) network and main active constituents of the α-glucosidase inhibitors from RS

Natural products have advantages of multi-components and multi-targets [60]. Their pharmacological effects are the consequence of a series of interactional biochemical reactions. On the other hand, T2D is a chronic degenerative disease involving various genetic and environmental factors [61]. These elements cause great difficulty in the related researches. Network methodology provides us a great opportunity to solve the problem from a systemic perspective. In the present study, a compound-target-pathway network was first built for the α-glucosidase inhibitors, target proteins and related pathways (Fig. 2). This network contained 216 nodes (N = 216) and 877 connections (E = 877). The nodes consisted of 13 α-glucosidase inhibitors (red nodes), 117 target proteins (yellow nodes), and 86 pathways (green nodes). The larger circles denote the nodes with the most connections, and that gray lines represent connections. Annotations of these nodes were listed in Additional file 3.

Fig. 2
figure 2

Compound-target-pathway (CTP) network of the potential α-glucosidase inhibitors from RS. The network consists of 13 red nodes (potential α-glucosidase inhibitors), 117 yellow nodes (target proteins), 86 green nodes (pathways), and 877 connections. The larger circles denote key nodes with the most connections. Node information is listed in Additional file 4. Gray lines represent connections

Several topological parameters were calculated to describe characteristics of the network. Degree is the most elementary characteristic for a node, which tells us how many directly connected neighbors a node holds. k of the CTP network was 8.12, demonstrating that an average of more than eight neighbors were connected with one node. On the other hand, degree distribution measures the diversity of a network. Obviously, a few nodes had numerous neighbors, whereas others only had a small number of connections. Figure 3a shows that the red and blue nodes had an uneven distribution, whereas green nodes exhibited approximate scale-free property. The difference of P(k) for α-glucosidase inhibitors, hit targets and related pathways might be due to the complexity of natural products. The data suggested that a few highly connected nodes existed in the CTP network.

Fig. 3
figure 3

a Degree distribution of the CTP network. k represents degree values, and that P(k) indicates degree distribution. b Degree values (k) of all nodes in the CTP network, ranked in a descending order

As interactions between molecules play a critical role in modulating the intrinsic biological processes, more attention should be paid to the highly-connected elements [62]. Within the framework of network science, the nodes with most connections are generally considered as hubs [63]. Although hubs are few in number, they are generally positioned to make strong contributions to global network function [63]. Disturbances on these hubs would spread rapidly throughout the entire network. Figure 3b shows a descending order of degree values for all nodes of the CTP network.

Among the α-glucosidase inhibitors, chrysin (k = 87) had the largest k, followed by 5,8,2′-trihydroxy-7-methoxyflavone (k = 21), and wogonin (k = 20). The three compounds exhibited much higher degree values than average (k = 8.12). Chrysin treatment has been used to improve diabetes in rats, which could attenuate diabetes-induced impairment [64]. Wogonin could enhance the intracellular level of adiponectin, a therapeutic target for insulin resistance, diabetes, and diabetes-related complications [65]. Although there were few reports about bioactivity of 5,8,2′-trihydroxy-7-methoxyflavone, the activities of flavonoids are structure dependent with the hydroxylated phenolic structure [66], which should be tested in the future. These compounds have significant impacts on the global network function, and were considered as main active ingredients. They might contribute the most to the pharmacological effects of α-glucosidase inhibitors from RS.

Multi-targets could provide superior therapeutic effect and less side effect compared to a single action, especially in the treatment of complex diseases [67]. In the target proteins, mitogen-activated protein kinase 3 (MAPK3, k = 60), mitogen-activated protein kinase 1 (MAPK1, k = 59) and phosphoinositide-3-kinase regulatory subunit 1 (PIK3R1, k = 49) had much higher degree values than others. The two MAPKs have been found to be increased in human and rodent adipose tissue in the diabetic state [68]. PIK3R1 plays a key role in insulin signaling and diabetes [69]. The three proteins were interconnected with the most compounds and pathways in the CTP network, and were also considered as hub nodes.

Normal pathways maintain balance between complex intracellular and intercellular networks. The most highly connected pathway in the CTP network was metabolic pathways (hsa01100, k = 29), followed by pathways in cancer (hsa05200, k = 21), metabolism of xenobiotics by cytochrome P450 (hsa00980, k = 15) and PI3K-Akt signaling pathway (hsa04151, k = 15). Hsa01100 is a mega pathway defined in Kyoto Encyclopedia of Genes and Genomes (KEGG), that encompasses several other pathways [70], and was excluded to avoid redundant data. Hsa05200 is related to many diseases. A cross talk is existed between diabetes and obesity, and that diabetes has been shown to increase cancer risk [71]. Some drugs appear to reduce cancer incidence and improve prognosis of patients with diabetes [72]. Hsa00980 takes part in biotransformation of medicines in vivo [73]. Hsa04151 is required for insulin-dependent regulation on cellular metabolism, which was directly associated with T2D [74]. These reports indicated that the α-glucosidase inhibitors from RS might be work through various pathways.

Protein-protein interaction (PPI) weighted network and key targets of the α-glucosidase inhibitors from RS

Many biological systems found in biology contain numerous components, and the interactions between individual agents cause the emergence of structures and functions [75]. T2D is a typical polygenic disease affected by various therapeutic targets [76]. Elucidation of the interactions between targets of T2D and ligands of α-glucosidase inhibitors from RS would help to understand the molecular mechanisms. In this study, a total of 64 targets of T2D were collected, including 34 successful and 30 clinical trial targets. Coincidentally, nine of them were also the targets of α-glucosidase inhibitors from RS, such as bile acid receptor (NR1H4), histone deacetylase 1 (HDAC1), prostaglandin G/H synthase 2 (PGH2), and so on. We analyzed the functional associations between the two groups of proteins using STRING database. After preliminary exclusion of isolated nodes, 69 ligands of α-glucosidase inhibitors from RS and 35 drug-targets of T2D were preserved. A protein-protein interaction (PPI) weighted network was then constructed (Fig. 4) for the two groups of proteins, containing 104 nodes and 228 connections.

Fig. 4
figure 4

Protein-protein interaction (PPI) weighted network for the ligands of α-glucosidase inhibitors from RS and targets of T2D, containing 104 nodes and 228 connections. The yellow nodes are the targets of potential α-glucosidase inhibitors from RS, and that blue nodes represent therapeutic targets of T2D

In this PPI weighted network, degree ki represents the number of other proteins interacted with node i. Node strength (si) describes the interactive intensities between the two groups of proteins. This parameter comprehensively reflects local information of node i by considering neighbor nodes and edge weights. Correlation between degree k and the average strength (k) for nodes with specific k values were investigated. If s(k)  kβ with an exponent β ≠ 1, the edge weight is correlated with the network topology. Figure 5a depicts the correlation between k and s(k) for the PPI network. The s(k) increased with k as s(k)  kβ with the exponent β ≈ 0.87, indicating that the node strengths were closely associated with degrees.

Fig. 5
figure 5

a Average strength s(k) as a function of degree k in logarithmic coordinates. The data points are fitted to a straight line, showing the relation s(k)  kβ.b Node strengths of the PPI network sorted in descending order. c Disparity Y(k) in the edge weight decays as a function of k. The data points are well approximated by the curve Y(k) = 1/k

Node strengths of the PPI network were then sorted in descending order (Fig. 5b). Larger strengths point to nodes with larger degrees. Nuclear factor kappa B subunit 1 (NFKB1) had the strongest interactions (s = 21.26) with the targets of the potential α-glucosidase inhibitors, followed by phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit gamma (PK3CG, s = 16.61), acetyl-CoA carboxylase beta (ACACB, s = 13.32), and insulin precursor (INS, s = 12.96). In the targets of the potential α-glucosidase inhibitors, PGH2 (s = 17.47) and PPARG (s = 14.66) interacted with the most nodes of T2D. Moreover, the two proteins were also successful targets of T2D. These nodes occupied important positions in the overall organization of PPI network.

System behaviors are dependent largely on the overall structure rather than the individual parts. Disparity Yi depicts the dispersion of weight distribution [77]. If the weights of all edges are approximately equal, Yi 1/ki. If one edge weight is important whereas the others negligible, Yi ≈ 1. It is obvious that Yi is associated with ki. Y(k) often attracts more attentions in the weighted network. It is the average of Yi for all nodes with degree k. If the weight distribution is relatively uniform, Y(k) 1/k, otherwise Y(k) ≈ 1. As shown in Fig. 5c, the average disparity Y(k) 1/k in the PPI network. It demonstrated that the edge weight distribution for nodes with the same degree k was approximately equal.

Changes in central positions of the network are always more important than those in marginal or relatively isolated positions [78]. To determine central nodes of the PPI network, three centrality indices, degree centrality, betweenness centrality, and closeness centrality were integrated into a three-dimensional diagram (Fig. 6). Distribution of these values seemed roughly uniform. However, a few nodes appeared as outliers. PPARG (Cd = 0.165, Cb = 0.232, Cc = 0.401), ACACB (Cd = 0.155, Cb = 0.184, Cc = 0.318), NFKB1 (Cd = 0.233, Cb = 0.161, Cc = 0.431), and PGH2 (Cd = 0.194, Cb = 0.157, Cc = 0.427) showed higher centrality scores than other nodes. A total of 54 neighbors were found to be connected with these central nodes, accounting for 51.9% of the total target proteins. The highly connectivity of the four proteins indicated that they could affect the PPI weighted network greatly.

Fig. 6
figure 6

Three-dimensional diagram of degree centrality (Cd), betweenness centrality (Cb) and closeness centrality (Cc) for the nodes in PPI network

PPARG and PGH2 are both important targets of the α-glucosidase inhibitors and T2D. Recent studies have demonstrated the association of PPARG with T2D. PPARG is a master transcriptional regulator of adipocyte differentiation. Variants in PPARG with decreased activity in adipocyte differentiation were found to be associated with increased risk of T2D [79]. A family-based study of Mexican Americans showed that variation in PPARG contributes to declining insulin resistance and concomitant deterioration in β-cell function at risk for T2D [80]. PGH2 generates prostaglandins and causes insulin insensitivity. PGH2 polymorphisms were found to be associated with T2D in Pima Indians comprising 1000 subjects [81]. Another variant of PGH2 had a protective role against T2D in two German cohorts [82].

ACACB and NFKB1 are therapeutic targets of T2D, and have a strong relationship with targets of α-glucosidase inhibitors from RS. ACACB is a regulator of fatty acid metabolism. It catalyzes the carboxylation of acetyl-CoA to malonyl-CoA. The problems in fatty acid metabolism can lead to insulin resistance, which is a precursor for T2D. Polymorphisms of ACACB are associated with T2D in postmenopausal women and Pakistani Punjabis [83, 84]. NFKB1 encodes a subunit of NF-kappa B. It is specifically involved in anti-inflammatory effects, and that inflammation is linked to insulin resistance and diabetes. Two common NFKB1 variants were found to be involved in T2D in an elderly cohort [85]. A transcriptome and proteome study demonstrated that NFKB1 were increased in expression in diabetic subjects [86]. These reports further confirmed the importance of PPARG, PGH2, ACACB, and NFKB1 to the α-glucosidase inhibitors from RS.

Core subnetwork and critical pathways of the α-glucosidase inhibitors from RS

In order to get further understanding of the key targets, the nodes connected to PPARG, ACACB, NFKB1, and PGH2 were extracted from PPI network. A total of 45 targets of the α-glucosidase inhibitors and 13 drug targets of T2D were selected. Pathway analysis indicated that these proteins were involved in 91 pathways. After querying KEGG DISEASE database, the type II diabetes mellitus pathway (hsa04930) and PPAR signaling pathway (hsa03320) showed a direct correlation with T2D. Therefore, the two processes might play significant roles in the pharmacological activities of the α-glucosidase inhibitors from RS. Moreover, we also extracted the α-glucosidase inhibitors connected to these key targets from CTP network, including chrysin, 5,8,2′-trihydroxy-7-methoxyflavone, and wogonin. All these elements were reorganized as a core subnetwork (Fig. 7), consisted of 63 nodes and 220 connections.

Fig. 7
figure 7

Core subnetwork for the potential α-glucosidase inhibitors and type 2 diabetes mellitus, consisted of 29 nodes and 47 connections. The yellow nodes are the targets of potential α-glucosidase inhibitors from RS, and that blue nodes indicate therapeutic targets of T2D. Red and green nodes represent the related α-glucosidase inhibitors and pathways, respectively

In this core subnetwork, chrysin, 5,8,2′-trihydroxy-7-methoxyflavone, and seven targets, MAPK1, MAPK3, PIK3R1, protein kinase C delta (PRKCD), INS, insulin receptor (INSR), solute carrier family 2 member 4 (GLUT4) belonged to hsa04930. As Fig. 8a shown, type II diabetes mellitus contains various kinases. MAPK1 and MAPK3 (also known as ERK2 and ERK1) play an important role in the MAPK/ERK cascade. Diabetogenic factors have been found to affect insulin signaling through activation of the ERK signaling pathway [87]. Previous research has revealed that targeting of the ERK pathway held promise for the treatment of T2D [88]. In the present study, MAPK1 and MAPK3 were the targets of chrysin and 5,8,2′-trihydroxy-7-methoxyflavone. In addition, PIK3R1 is necessary for the insulin-stimulated increase in glucose uptake and glycogen synthesis in insulin-sensitive tissues. Mutations of PIK3R1 could cause insulin resistance, which is strongly associated with insulin resistance [69]. Interestingly, chrysin, 5,8,2′-trihydroxy-7-methoxyflavone, MAPK1, MAPK3 and PIK3R1 were also hub nodes of the CTP network.

Fig. 8
figure 8

Critical pathways of the potential α-glucosidase inhibitors from RS. a Type 2 diabetes mellitus pathway. b PPAR signaling pathway. The yellow nodes are the targets of potential α-glucosidase inhibitors from RS, blue nodes indicate therapeutic targets of T2D, and that pink nodes denote targets belonged to both the two groups. Red nodes represent the related α-glucosidase inhibitors

Chrysin, 5,8,2′-trihydroxy-7-methoxyflavone, and another six targets, peroxisome proliferator activated receptor alpha (PPARA), peroxisome proliferator activated receptor delta (PPARD), PPARG, retinoid X receptor alpha (RXRA), 3-phosphoinositide dependent protein kinase 1 (PDPK1), stearoyl-CoA desaturase (ACOD) were involved in hsa03320 (Fig. 8b). Accumulating evidence highlighted the role of PPAR signaling in T2D [89]. PPARs (Peroxisome proliferator-activated receptors) are nuclear hormone receptors that are activated by fatty acids and their derivatives, containing three subtypes (PPAR alpha, beta/delta, and gamma). The three PPAR isoforms all appeared in hsa03320. They were both the targets of the α-glucosidase inhibitors from RS and T2D. Moreover, PPARG was the central node of PPI weight network. It could promote adipocyte differentiation to enhance blood glucose uptake. In hsa03320 process, chrysin and 5,8,2′-trihydroxy-7-methoxyflavone were directly connected with the PPARs. This further confirmed the importance of the two α-glucosidase inhibitors. These data supported the hypothesis that the α-glucosidase inhibitors from RS might influence T2D through hsa04930 and hsa03320 processes.

Discussion

Many previous studies have investigated the antidiabetic constituents of Radix scutellariae. For instance, Cui et al. used a tandem quadrupole mass spectrometer coupled with enzyme activity analysis to explore hypoglycemic effect of RS and Coptidis rhizome [90]. This method was accurate and sensitive enough for quantitative evaluation of seven major components and six enzymes. The results indicated that combined extract had stronger effects on T2D through multiple components against multiple targets. Tahtah et al. applied triple aldose reductase/α-glucosidase/radical scavenging high-resolution profiling combined with high-performance liquid chromatography-high-resolution mass spectrometry-solid-phase extraction-nuclear magnetic resonance spectroscopy to identify antidiabetic constituents from RS [91]. Baicalein was screened as α-glucosidase inhibitor. In another approach, a systematic study on metabolism and activity evaluation of Radix Scutellaria extract in rat plasma was conducted, using UHPLC with quadrupole time-of-flight mass spectrometry [92]. Wogonoside, norwogoin-7-O-Glu acid and oroxyloside exhibited better binding affinities with α-glucosidase. These results are partially consistent with those obtained in our study. While they demonstrate versatility and success of phytochemical analysis in the identification of novel ligands for therapeutic targets, these studies are labor-intensive and time consuming [55]. Moreover, the mechanism analysis of active constituents is often limited by existing knowledges and experiences. Our study constructs a network model of compounds, target proteins and pathways to explore mechanism of α-glucosidase inhibitors identified by ultrafiltration LC-MS from RS. This approach is more rapidly and extensive as the application of computational tools as well as systems biology. The selected main active components, key targets and critical pathways would provide more information for the interpretation of Radix Scutellaria and T2D. However, these results are mainly based on statistical analysis and prediction. Further studies on cell and animal models are required to give a distinct answer.

Conclusion

This study presents one of the first systematic analyses of α-glucosidase inhibitors from natural products using ultrafiltration LC-MS/MS and network pharmacology. Our findings suggested a possible application of flavonoids from Radix scutellariae in the treatment of T2D. The n-butanol part of ethanol extract from RS showed strong α-glucosidase inhibition activity. Thirty-two kinds of flavonoids were identified from the extract, and 13 of them were screened as α-glucosidase inhibitors, including viscidulin III, oroxylin A, 2′,3,5,6′,7-pentahydroxyflavanone, and so on. The results were strongly supported by previous reports about natural flavonoids and T2D [4, 5, 14, 22, 66]. These compounds, together with their target proteins and related pathways, were integrated into three complex networks. The underlying mechanism of these natural α-glucosidase inhibitors were revealed by network analyses. Chrysin, 5,8,2′-trihydroxy-7-methoxyflavone and wogonin were the main active constituents. PPARG, PGH2, ACACB, and NFKB1 were key targets. The α-glucosidase inhibitors from RS might influence T2D progression through the type II diabetes mellitus and PPAR signaling pathways. In all, the combination of ultrafiltration LC-MS and system pharmacology would enable to generate novel insight into the research of plant medicines.

Availability of data and materials

Datasets supporting the results of this article have been included in the additional files.

Abbreviations

ACACB:

Acetyl-CoA carboxylase beta

ACOD:

Stearoyl-CoA desaturase

CI:

Centrality indices

CTP:

Compound-target-pathway

GLUT4:

Solute carrier family 2 member 4

HDAC1:

Histone deacetylase 1

hsa00980:

Metabolism of xenobiotics by cytochrome P450

hsa01100:

Metabolic pathways

hsa03320:

PPAR signaling pathway

hsa04151:

PI3K-Akt signaling pathway

hsa04930:

Type II diabetes mellitus pathway

INS:

Insulin precursor

INSR:

Insulin receptor

KEGG:

Kyoto Encyclopedia of Genes and Genomes

MAPK1:

Mitogen-activated protein kinase 1

MAPK3:

Mitogen-activated protein kinase 3

NFKB1:

Nuclear factor kappa B subunit 1

NR1H4:

Bile acid receptor

PDPK1:

3-phosphoinositide dependent protein kinase 1

PGH2:

Prostaglandin G/H synthase 2

PIK3R1:

Phosphoinositide-3-kinase regulatory subunit 1

PK3CG:

Phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit gamma

PPARA:

Peroxisome proliferator activated receptor alpha

PPARD:

Peroxisome proliferator activated receptor delta

PPARG:

Peroxisome proliferator activated receptor gamma

PPI:

Protein-protein interaction

PRKCD:

Protein kinase C delta

RS :

Radix scutellariae

RXRA:

Retinoid X receptor alpha

STRING:

Search Tool for the Retrieval of Interacting Genes/Proteins

T2D:

Type 2 diabetes

Triple TOF-MS:

Triple time-of-flight mass spectrometer

UPLC:

Ultra high performance liquid chromatography

References

  1. Chatterjee S, Khunti K, Davies MJ. Type 2 diabetes. Lancet. 2017;58(1):63.

    Google Scholar 

  2. Hays NP, Galassetti PR, Coker RH. Prevention and treatment of type 2 diabetes: current role of lifestyle, natural product, and pharmacological interventions. Pharmacol Ther. 2008;118(2):181–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Shapiro K, Gong WC. Natural products used for diabetes. J Am Pharm Assoc. 2002;42(2):217–26.

    Google Scholar 

  4. Alam F, Islam MA, Kamal MA, Gan SH. Updates on managing type 2 diabetes mellitus with natural products: towards antidiabetic drug development. Curr Med Chem. 2016;25(39):5395–431.

    Article  CAS  Google Scholar 

  5. Salimifar M, Fatehi-Hassanabad Z, Fatehi M. A review on natural products for controlling type 2 diabetes with an emphasis on their mechanisms of actions. Curr Diabetes Rev. 2013;9(5):402–11.

    Article  CAS  PubMed  Google Scholar 

  6. Li C, Lin G, Zuo Z. Pharmacological effects and pharmacokinetics properties of Radix Scutellariae and its bioactive flavones. Biopharm Drug Dispos. 2011;32(8):427–45.

    Article  PubMed  CAS  Google Scholar 

  7. Luo J-L, Lu F-L, Liu Y-C, Lo C-F. Identification of Scutellaria Baicalensis in traditional Chinese medicine preparations by LC/MS/MS fingerprinting method. J Food Drug Anal. 2012;20(4):887–899+984.

    CAS  Google Scholar 

  8. Liu G, Ma J, Chen Y, Tian Q, Shen Y, Wang X, Chen B, Yao S. Investigation of flavonoid profile of Scutellaria bacalensis Georgi by high performance liquid chromatography with diode array detection and electrospray ion trap mass spectrometry. J Chromatogr A. 2009;1216(23):4809–14.

    Article  CAS  PubMed  Google Scholar 

  9. Zhao T, Tang H, Xie L, Zheng Y, Ma Z, Sun Q, Li X. Scutellaria baicalensis Georgi.(Lamiaceae): a review of its traditional uses, botany, phytochemistry, pharmacology and toxicology. J Pharm Pharmacol. 2019;71(9):1353–69.

    Article  CAS  PubMed  Google Scholar 

  10. Song KH, Lee SH, Kim BY, Park AY, Kim JY. Extracts of Scutellaria baicalensis reduced body weight and blood triglyceride in db/db mice. Phytother Res. 2013;27(2):244–50.

    Article  PubMed  Google Scholar 

  11. Park JH, Kim RY, Park E. Antioxidant and α-glucosidase inhibitory activities of different solvent extracts of skullcap (Scutellaria baicalensis). Food Sci Biotechnol. 2011;20(4):1107–12.

    Article  CAS  Google Scholar 

  12. Liu SZ, Deng YX, Chen B, Zhang XJ, Shi QZ, Qiu XM. Antihyperglycemic effect of the traditional Chinese scutellaria-coptis herb couple and its main components in streptozotocin-induced diabetic rats. J Ethnopharmacol. 2013;145(2):490–8.

    Article  CAS  PubMed  Google Scholar 

  13. Waisundara VY, Hsu A, Huang D, Tan BK. Scutellaria baicalensis enhances the anti-diabetic activity of metformin in streptozotocin-induced diabetic Wistar rats. Am J Chin Med. 2008;36(3):517.

    Article  CAS  PubMed  Google Scholar 

  14. Park J, Jang HJ. Anti-diabetic effects of natural products an overview of therapeutic strategies. Mol Cell Toxicol. 2017;13(1):1–20.

    Article  CAS  Google Scholar 

  15. Bischoff H. The mechanism of alpha-glucosidase inhibition in the management of diabetes. Clin Invest Med. 1995;18(4):303.

    CAS  PubMed  Google Scholar 

  16. Yang Y, Lian G, Yu B. Naturally occurring polyphenolic glucosidase inhibitors. Cheminform. 2015;55(3–4):268–84.

    CAS  Google Scholar 

  17. Chen J, Cheng YQ, Yamaki K, Li LT. Anti-α-glucosidase activity of Chinese traditionally fermented soybean (douchi). Food Chem. 2007;103(4):1091–6.

    Article  CAS  Google Scholar 

  18. Ríos JL, Francini F, Schinella GR. Natural products for the treatment of type 2 diabetes mellitus. Planta Med. 2015;81(12–13):975.

    PubMed  Google Scholar 

  19. Habtemariam S. Targeting intestinal digestive enzymes by natural products: synergistic effect of flavonoids. Planta Med. 2011;77(12):1404.

    Article  Google Scholar 

  20. Shibano M, Kakutani K, Taniguchi M, Yasuda M, Baba K. Antioxidant constituents in the dayflower (Commelina communis L.) and their alpha-glucosidase-inhibitory activity. J Nat Med. 2008;62(3):349–53.

    Article  CAS  PubMed  Google Scholar 

  21. Matsuda H, Morikawa T, Toguchida I, Yoshikawa M. Structural requirements of flavonoids and related compounds for aldose reductase inhibitory activity. Chem Pharm Bull (Tokyo). 2002;50(6):788.

    Article  CAS  Google Scholar 

  22. Xu L-N, Li Y, Dai Y, Peng J-Y. Natural products for the treatment of type 2 diabetes mellitus: pharmacology and mechanisms. Pharmacol Res. 2018;130:451.

    Article  CAS  PubMed  Google Scholar 

  23. Yang JR, Luo JG, Kong LY. Determination of α-glucosidase inhibitors from Scutellaria baicalensis using liquid chromatography with quadrupole time of flight tandem mass spectrometry coupled with centrifugal ultrafiltration. Chin J Nat Med. 2015;13(3):208–14.

    CAS  PubMed  Google Scholar 

  24. Kishida K. Simplified extraction of tetracycline antibiotics from milk using a centrifugal ultrafiltration device. Food Chem. 2011;126(2):687–90.

    Article  CAS  Google Scholar 

  25. Yang Z, Zhang Y, Sun L, Wang Y, Gao X, Cheng Y. An ultrafiltration high-performance liquid chromatography coupled with diode array detector and mass spectrometry approach for screening and characterising tyrosinase inhibitors from mulberry leaves. Anal Chim Acta. 2012;719(6):87.

    Article  CAS  PubMed  Google Scholar 

  26. Zhang Y, Peng M, Liu L, Shi S, Peng S. Screening, identification, and potential interaction of active compounds from Eucommia ulmodies leaves binding with bovine serum albumin. J Agric Food Chem. 2012;60(12):3119–25.

    Article  CAS  PubMed  Google Scholar 

  27. Wang J, Liu S, Li S, Song F, Zhang Y, Liu Z, Liu C-M. Ultrafiltration LC-PDA-ESI/MS combined with reverse phase-medium pressure liquid chromatography for screening and isolation potential α-glucosidase inhibitors from Scutellaria baicalensis Georgi. Anal Methods. 2014;6(15):5918.

    Article  CAS  Google Scholar 

  28. Beutler JA. Natural products as a foundation for drug discovery. Curr Protoc Pharmacol. 2009;46(1):1–21.

    Article  Google Scholar 

  29. Lam KS. New aspects of natural products in drug discovery. Trends Microbiol. 2007;15(6):279–89.

    Article  CAS  PubMed  Google Scholar 

  30. Hopkins AL. Network pharmacology: the next paradigm in drug discovery. Nat Chem Biol. 2008;4(11):682–90.

    Article  CAS  PubMed  Google Scholar 

  31. Zhang X-Z, Gu J-Y, Cao L, Li N, Ma Y-M, Su Z-Z, Ding G, Chen L-R, Xu X-J, Xiao W. Network pharmacology study on the mechanism of traditional Chinese medicine for upper respiratory tract infection. Mol BioSyst. 2014;10(10):2517.

    Article  CAS  PubMed  Google Scholar 

  32. Li B-H, Tao W-Y, Zheng C-L, Shar PA, Huang C, Fu Y-X, Wang Y-H. Systems pharmacology-based approach for dissecting the addition and subtraction theory of traditional Chinese medicine: an example using Xiao-Chaihu-Decoction and Da-Chaihu-Decoction. Comput Biol Med. 2014;53:19–29.

    Article  PubMed  Google Scholar 

  33. Hopkins AL. Network pharmacology. Nat Biotechnol. 2007;25(10):1110–1.

    Article  CAS  PubMed  Google Scholar 

  34. Wang L-L, Li Z, Shao Q, Li X, Ai N, Zhao X-P, Fan X-H. Dissecting active ingredients of Chinese medicine by content-weighted ingredient-target network. Mol BioSyst. 2014;10(7):1905.

    Article  CAS  PubMed  Google Scholar 

  35. Chen L, Lv D, Wang D, Chen X, Zhu Z, Cao Y, Chai Y. A novel strategy of profiling the mechanism of herbal medicines by combining network pharmacology with plasma concentration determination and affinity constant measurement. Mol BioSyst. 2016;12(11):3347.

    Article  CAS  PubMed  Google Scholar 

  36. Luo F, Gu J, Chen L, Xu X. Systems pharmacology strategies for anticancer drug discovery based on natural products. Mol BioSyst. 2014;10(7):1912.

    Article  CAS  PubMed  Google Scholar 

  37. Gogoi B, Gogoi D, Silla Y, Kakoti BB, Bhau BS. Network pharmacology-based virtual screening of natural products from Clerodendrum species for identification of novel anti-cancer therapeutics. Mol BioSyst. 2017;13(2):406–16.

    Article  CAS  PubMed  Google Scholar 

  38. Kibble M, Saarinen N, Tang J, Wennerberg K, Mäkelä S, Aittokallio T. Network pharmacology applications to map the unexplored target space and therapeutic potential of natural products. Nat Prod Rep. 2015;32(8):1249–66.

    Article  CAS  PubMed  Google Scholar 

  39. Lagunin AA, Goel RK, Gawande DY, Pahwa P, Gloriozova TA, Dmitriev AV, Ivanov SM, Rudik AV, Konova VI, Pogodin PV. Chemo- and bioinformatics resources for in silico drug discovery from medicinal plants beyond their traditional use: a critical review. Nat Prod Rep. 2014;31(11):1585–611.

    Article  CAS  PubMed  Google Scholar 

  40. Kang WY, Song YL, Zhang L. α-Glucosidase inhibitory and antioxidant properties and antidiabetic activity of Hypericum ascyron L. Med Chem Res. 2010;20(7):809–16.

    Article  CAS  Google Scholar 

  41. Nickel J, Gohlke BO, Erehman J, Banerjee P, Rong WW, Goede A, Dunkel M, Preissner R. SuperPred: update on drug classification and target prediction. Nucleic Acids Res. 2014;42(1):26–31.

    Article  CAS  Google Scholar 

  42. Benjamini Y, Hochberg Y. Controlling the false discovery rate - a practical and powerful approach to multiple testing. J R Stat Soc B. 1995;57(1):289–300.

    Google Scholar 

  43. Li YH, Yu CY, Li XX, Zhang P, Tang J, Yang Q, Fu T, Zhang X, Cui X, Tu G. Therapeutic target database update 2018: enriched resource for facilitating bench-to-clinic research of targeted therapeutics. Nucleic Acids Res. 2017;46(1):1121–7.

    Google Scholar 

  44. Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, Santos A, Doncheva NT, Roth A, Bork P. The STRING database in 2017: quality-controlled protein–protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45(Database issue):D362–8.

    Article  CAS  PubMed  Google Scholar 

  45. Zhang L, Zhang R-W, Li Q, Lian J-W, Liang J, Chen X-H, Bi K-S. Development of the fingerprints for the quality evaluation of Scutellariae Radix by HPLC-DAD and LC-MS-MS. Chromatographia. 2007;66(1–2):13–20.

    Article  CAS  Google Scholar 

  46. Liu G, Ma J, Chen Y, Tian Q, Shen Y, Wang X, Yao S. Investigation of flavonoid profile of Scutellaria bacalensis Georgi by high performance liquid chromatography with diode array detection and electrospray ion trap mass spectrometry. J Chromatogr. 2009;1216(23):4809–14.

    Article  CAS  Google Scholar 

  47. Ma YL, Li QM, Van den Heuvel H, Claeys M. Characterization of flavone and flavonol aglycones by collision-induced dissociation tandem mass spectrometry. Rapid Commun Mass Spectrom. 2015;11(12):1357–64.

    Article  Google Scholar 

  48. Reddy B, Reddy M, Gunasekar D, Murthy MM, Caux C, Bodo B. Two new flavonoids from Andrographis macrobotrys. Indian J Chem B. 2005;44(9):1966–9.

    Google Scholar 

  49. He L, Zhang Z, Lu L, Liu Y, Li S, Wang J, Song Z, Yan Z, Miao J. Rapid identification and quantitative analysis of the chemical constituents in Scutellaria indica L. by UHPLC–QTOF–MS and UHPLC–MS/MS. J Pharm Biomed Anal. 2016;117:125–39.

    Article  CAS  PubMed  Google Scholar 

  50. Azimova SS, Vinogradova VI. Natural compounds: flavonoids. New York: Springer; 2013. p. 104–18.

  51. Xiao L, Wang H, Song S, Zhang G, Song H, Xu S. Isolation and identification of the chemical constituents of roots of Scutellaria amoena CH Wright. Shenyang Yaoke Daxue Xuebao. 2003;20(3):181–3.

    CAS  Google Scholar 

  52. Barberan F, Ferreres F, Tomas F, Guirado A. Electron impact mass spectrometric differentiation of 5, 6-dihydroxy-7, 8-dimethoxy-and 5, 8-dihydroxy-6, 7-dimethoxyflavones. Phytochemistry. 1986;25(4):923–5.

    Article  CAS  Google Scholar 

  53. Han J, Ye M, Xu M, Sun J, Wang B, Guo D. Characterization of flavonoids in the traditional Chinese herbal medicine-Huangqin by liquid chromatography coupled with electrospray ionization mass spectrometry. J Chromatogr B Analyt Technol Biomed Life Sci. 2007;848(2):355–62.

    Article  CAS  PubMed  Google Scholar 

  54. Guo M, Zhang Y, Wang Y, Zhao X. Development of a rapid screening method for discovering neuroprotective components from traditional Chinese medicine. Zhongguo Zhong Yao Za Zhi. 2013;38(10):1581–4.

    CAS  PubMed  Google Scholar 

  55. Cieśla Ł, Moaddel R. Comparison of analytical techniques for the identification of bioactive compounds from natural products. Nat Prod Rep. 2016;33(10):1131–45.

    Article  PubMed  Google Scholar 

  56. Szklarczyk D, Santos A, Von Mering C, Jensen LJ, Bork P, Kuhn M. STITCH 5: augmenting protein–chemical interaction networks with tissue and affinity data. Nucleic Acids Res. 2016;44(1):380–4.

    Article  CAS  Google Scholar 

  57. Moller DE. New drug targets for type 2 diabetes and the metabolic syndrome. Nature. 2001;414(6865):821–7.

    Article  CAS  PubMed  Google Scholar 

  58. Kaur P, Mittal A, Nayak S, Vyas M, Mishra V, Khatik G. Current strategies and drug targets in the management of type 2 diabetes mellitus. Curr Drug Targets. 2018;19(15):1738–66.

    Article  CAS  PubMed  Google Scholar 

  59. Carmeli B, Erhan B, Koyama T, Rhrissorrakrai K, Royyuru AK, Utro F, Waks Z. Relevancy assessment and visualization of biological pathways. In: U.S. patent application, vol. 14/745; 2016. p. 616.

    Google Scholar 

  60. Wang X-J, Zhang A-H, Sun H, Han Y, Yan G-L. Discovery and development of innovative drug from traditional medicine by integrated chinmedomics strategies in the post-genomic era. TrAC Trends Anal Chem. 2016;76:86–94.

    Article  CAS  Google Scholar 

  61. Martínez Gómez LE, Cruz M, Martínez Nava GA, Madrid Marina V, Parra E, García Mena J, Espinoza Rojo M, Estrada Velasco BI, Piza Roman LF, Aguilera P. A replication study of the IRS1, CAPN10, TCF7L2, and PPARG gene polymorphisms associated with type 2 diabetes in two different populations of Mexico. Ann Hum Genet. 2011;75(5):612–20.

    Article  PubMed  CAS  Google Scholar 

  62. Wang L, Tan N, Hu J, Wang H, Duan D, Ma L, Xiao J, Wang X. Analysis of the main active ingredients and bioactivities of essential oil from Osmanthus fragrans Var. thunbergii using a complex network approach. BMC Syst Biol. 2017;11(1):144.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  63. Mp VDH, Sporns O. Network hubs in the human brain. Trends Cogn Sci. 2013;17(12):683–96.

    Article  Google Scholar 

  64. Samarghandian PS, Nezhad MA, Samini F, Farkhondeh T. Chrysin treatment improves diabetes and its complications in liver, brain, and pancreas in streptozotocin-induced diabetic rats. Can J Physiol Pharmacol. 2016;94(4):388–93.

    Article  CAS  PubMed  Google Scholar 

  65. Yang T, Liu H, Zhao B, Xia Z, Zhang Y, Zhang D, Li M, Cao Y, Zhang Z, Bi Y. Wogonin enhances intracellular adiponectin levels and suppresses adiponectin secretion in 3T3-L1 adipocytes. Endocr J. 2016;64(1):15–26.

    Article  PubMed  Google Scholar 

  66. Bobilya DJ. Flavonoid antioxidants: chemistry, metabolism and structure-activity relationships. J Nutr Biochem. 2002;13(10):572–84.

    Article  PubMed  Google Scholar 

  67. Li Y, Huang TH, Yamahara J. Salacia root, a unique Ayurvedic medicine, meets multiple targets in diabetes and obesity. Life Sci. 2008;82(21):1045–9.

    Article  CAS  PubMed  Google Scholar 

  68. Bashan N, Dorfman K, Tarnovscki T, Harman Boehm I, Liberty IF, Blüher M, Ovadia S, Maymon Zilberstein T, Potashnik R, Stumvoll M. Mitogen-activated protein kinases, inhibitory-κB kinase, and insulin signaling in human omental versus subcutaneous adipose tissue in obesity. Endocrinology. 2007;148(6):2955–62.

    Article  CAS  PubMed  Google Scholar 

  69. Thauvinrobinet C, Auclair M, Duplomb L, Carondebarle M, Avila M, Stonge J, Merrer ML, Luyer BL, Héron D, Mathieudramard M. PIK3R1 mutations cause syndromic insulin resistance with lipoatrophy. Am J Hum Genet. 2013;93(1):141–9.

    Article  CAS  Google Scholar 

  70. Liu C, Bousman C, Pantelis C, Skafidas E, Zhang D, Yue W, Everall I. Pathway-wide association study identifies five shared pathways associated with schizophrenia in three ancestral distinct populations. Transl Psychiatry. 2017;7(2):e1037.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. La LP, Listì A, Caruso S, Amodeo V, Passiglia F, Bazan V, Fanale D. Potential role of ANGPTL4 in the cross talk between metabolism and cancer through PPAR signaling pathway. PPAR Res. 2017;2017:1–15.

    Google Scholar 

  72. Gallagher EJ, Leroith D. Diabetes, cancer, and metformin: connections of metabolism and cell proliferation. Ann N Y Acad Sci. 2011;1243(1):54–68.

    Article  CAS  PubMed  Google Scholar 

  73. Zanger UM, Schwab M. Cytochrome P450 enzymes in drug metabolism: regulation of gene expression, enzyme activities, and impact of genetic variation. Pharmacol Ther. 2013;138(1):103–41.

    Article  CAS  PubMed  Google Scholar 

  74. Schultze SM, Hemmings BA, Niessen M, Tschopp O. PI3K/AKT, MAPK and AMPK signalling: protein kinases in glucose homeostasis. Expert Rev Mol Med. 2012;14:e1.

    Article  PubMed  CAS  Google Scholar 

  75. Holcombe M, Adra S, Bicak M, Chin S, Coakley S, Graham AI, Green J, Greenough C, Jackson D, Kiran M. Modelling complex biological systems using an agent-based approach. Integr Biol. 2012;4(1):53–64.

    Article  CAS  Google Scholar 

  76. Papaetis GS, Papakyriakou P, Panagiotou TN. Central obesity, type 2 diabetes and insulin: exploring a pathway full of thorns. Arch Med Sci. 2015;11(3):463–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Barrat A, Barthelemy M, Vespignani A. The architecture of complex weighted networks: measurements and models. In: Large scale structure and dynamics of complex networks: from information technology to finance and natural science. edn. New Jersey: World Scientific; 2007. p. 67–92.

    Chapter  Google Scholar 

  78. Wang L, Hou E-T, Wang L-J, Wang Y-J, Yang L-J, Zheng X-H, Xie G-Q, Sun Q, Liang M-Y, Tian Z-M. Reconstruction and analysis of correlation networks based on GC-MS metabolomics data for young hypertensive men. Anal Chim Acta. 2015;854:95–105.

    Article  CAS  PubMed  Google Scholar 

  79. Stumvoll M, Häring H. The peroxisome proliferator-activated receptor-γ2 Pro12Ala polymorphism. Diabetes. 2002;51(8):2341–7.

    Article  CAS  PubMed  Google Scholar 

  80. Black MH, Wu J, Takayanagi M, Wang N, Taylor KD, Haritunians T, Trigo E, Lawrence JM, Watanabe RM, Buchanan TA. Variation in PPARG is associated with longitudinal change in insulin resistance in mexican americans at risk for type 2 diabetes. J Clin Endocrinol Metab. 2015;100(3):1187–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Konheim YL, Wolford JK. Association of a promoter variant in the inducible cyclooxygenase-2 gene (PTGS2) with type 2 diabetes mellitus in Pima Indians. Hum Genet. 2003;113(5):377–81.

    Article  CAS  PubMed  Google Scholar 

  82. Nitz I, Fisher E, Grallert H, Li Y, Gieger C, Rubin D, Boeing H, Spranger J, Lindner I, Schreiber S. Association of prostaglandin E synthase 2 (PTGES2) Arg298His polymorphism with type 2 diabetes in two German study populations. J Clin Endocrinol Metab. 2007;92(8):3183–8.

    Article  CAS  PubMed  Google Scholar 

  83. Zain M, Awan FR, Najam SS, Islam M, Khan AR, Bilal A, Bellili N, Marre M, Roussel R, Fumeron F. Association of ACACB gene polymorphism (rs2268388, G> A) with type 2 diabetes and end stage renal disease in Pakistani Punjabi population. Meta Gene. 2017;12:109–12.

    Article  Google Scholar 

  84. Riancho JA, Vázquez L, García Pérez MA, Sainz J, Olmos JM, Hernández JL, Pérez López J, Amado JA, Zarrabeitia MT, Cano A. Association of ACACB polymorphisms with obesity and diabetes. Mol Genet Metab. 2011;104(4):670–6.

    Article  CAS  PubMed  Google Scholar 

  85. Coto E, Díaz Corte C, Tranche S, Gómez J, Alonso B, Iglesias S, Reguero JR, López Larrea C, Coto Segura P. Gene variants in the NF-KB pathway (NFKB1, NFKBIA, NFKBIZ) and their association with type 2 diabetes and impaired renal function. Hum Immunol. 2018;79(6):494–8.

    Article  CAS  PubMed  Google Scholar 

  86. Stentz FB, Kitabchi AE. Transcriptome and proteome expressions involved in insulin resistance in muscle and activated T-lymphocytes of patients with type 2 diabetes. Genomics Proteomics Bioinformatics. 2007;5(3–4):216–35.

    Article  CAS  PubMed  Google Scholar 

  87. Tanti JF, Jager J. Cellular mechanisms of insulin resistance: role of stress-regulated serine kinases and insulin receptor substrates (IRS) serine phosphorylation. Curr Opin Pharm. 2009;9(6):753–62.

    Article  CAS  Google Scholar 

  88. Banks AS, Mcallister FE, Camporez JPG, Zushin PJH, Jurczak MJ, Laznikbogoslavski D, Shulman GI, Gygi SP, Spiegelman BM. An Erk/Cdk5 axis controls the diabetogenic actions of PPARγ. Nature. 2015;517(7534):391–5.

    Article  CAS  PubMed  Google Scholar 

  89. Fanale D, Amodeo V, Caruso S. The interplay between metabolism, PPAR signaling pathway, and cancer. PPAR Res. 2017;2017:1–2.

    Article  CAS  Google Scholar 

  90. Cui X, Shen YM, Jiang S, Qian DW, Shang EX, Zhu ZH, Duan JA. Comparative analysis of the main active components and hypoglycemic effects after the compatibility of Scutellariae Radix and Coptidis Rhizoma. J Sep Sci. 2019;42(8):1520–7.

    Article  CAS  PubMed  Google Scholar 

  91. Tahtah Y, Kongstad KT, Wubshet SG, Nyberg NT, Jønsson LH, Jäger AK, Qinglei S, Staerk D. Triple aldose reductase/α-glucosidase/radical scavenging high-resolution profiling combined with high-performance liquid chromatography–high-resolution mass spectrometry–solid-phase extraction–nuclear magnetic resonance spectroscopy for identification of antidiabetic constituents in crude extract of Radix Scutellariae. J Chromatogr. 2015;1408:125–32.

    Article  CAS  Google Scholar 

  92. Zhang X, Liu S, Xing J, Pi Z, Liu Z, Song F. Systematic study on metabolism and activity evaluation of Radix Scutellaria extract in rat plasma using UHPLC with quadrupole time-of-flight mass spectrometry and microdialysis intensity-fading mass spectrometry. J Sep Sci. 2018;41(7):1704–10.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

We are very grateful to Shaanxi Qinxi agricultural and Forestry Development Co., Ltd. (Baoji 721600, China) for technical support.

Funding

This study is supported and sponsored by the Shaanxi Science and Technology Commission (No. 2018JQ2064), Shaanxi Education Commission (No.16JS008), the Key Project of Baoji University of Arts and Sciences (No. ZK16053), the Project of Shuangbai & Industry-University-Research (No. 209011024), and the Innovation and Entrepreneurial Programme for Undergraduate (No. 2019XJ090).

Author information

Authors and Affiliations

Authors

Contributions

LW and NT collected data for the study. NT and WD performed data analysis. LW and JH conceived and designed the study. LW and HW wrote the manuscript. XW revised the manuscript. All authors have read and approved the manuscript.

Corresponding authors

Correspondence to Le Wang or Xiaoling Wang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1.

Total ions chromatogram (TIC) of the n-butanol part of ethanol extract from Radix Scutellariae by UPLC-TripleTOF.

Additional file 2.

Primary fragmentation pathways of compound 14 (oroxylin A-7-O-β-D-Glucuronopyranoside), 16 (wogonoside), 17 (4′-hydroxywogonin), 23 (5,8,2′-trihydroxy-7-methoxyflavone), 26 (wogonin) and 28 (chrysin).

Additional file 3.

α-glucosidase inhibition curves of crude extract of Radix scutellariae and acarbose.

Additional file 4.

Information of nodes in the compound-target-pathway (CTP) network.

Additional file 5.

Degree values of CTP network, node strength and centrality indices of PPI network.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, L., Tan, N., Wang, H. et al. A systematic analysis of natural α-glucosidase inhibitors from flavonoids of Radix scutellariae using ultrafiltration UPLC-TripleTOF-MS/MS and network pharmacology. BMC Complement Med Ther 20, 72 (2020). https://doi.org/10.1186/s12906-020-2871-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12906-020-2871-3

Keywords