3.1. Identification of compounds from PZE
The base peak chromatograms (BPC) of PZE in positive and negative modes are exhibited in Fig. 2. A total of 37 compounds were characterized using a self-constructed database analysis strategy, and the properties of their candidate compounds are shown in Fig. 2. It contained information on retention time, molecular formulas, mass weight, mass errors and fragments of major components, as shown in Table 1.
Insert Table 1. and Fig. 2.
3.2. Identifcation of compounds from PZE
To date, 37 kinds of constituents have been identified from the PZE, among them, xanthones, glycolipids and acids are the main types. All compounds are summarized and recorded in Table 1. Herein, we provide the analytical process of two kinds of representative constituents as the internal standard substances recorded in the Chinese Pharmacopoeia. As we can see from the Fig. 3, polygalaxanthone III (M2) was a classic xanthone from P. tenuifolia with molecular formula of C25H28O15 as inferred from the [M-H]− (m/z: 567.1342). In the second stage mass spectrum detected, the daughter ions at m/z 59.0143 (C2H3O2−), 89.0248 (C3H5O3−), 271.0254 (C14H7O6−), 273.0415 (C14H9O6−), 315.0504 (C16H11O7−), 345.0604 (C17H13O8−), 417.0816 (C20H17O10−), 435.0924 (C20H19O11−), 447.0915 (C21H19O11−) in negative mass mode provided the evidence for the identification of the component.
Insert Fig. 3.
As for the constituents from Z. jujuba, typical compound spinosin (M28) recorded in the Chinese Pharmacopoeia as internal standard substance was detected with the chemical formula of C28H32O15 based on the [M + H]+ ion at m/z 609.1784 and [M–H]− ion at m/z 607.1675. In the second stage mass spectrum (Fig. 4), the daughter ions at m/z 297.0752 (C17H13O5+), 327.0854 (C12H23O10+), 351.0847 (C16H15O9+), 411.1058 (C22H19O8+), 429.1156 (C22H21O9+), 447.1262 (C22H23O10+) in positive mass mode and 427.1041 (C22H19O9−) in positive mass mode give the inspiration of the spectrum analysis.
Insert Fig. 4.
3.3. Acquisition of targets for PZE ingredients and anti-epilepsy
The results of the network construction showed: Screening of 24 core anti-epilepsy components from andidate ingredients by preliminary search. Generally, a total of 527 ingredient-related targets were identified from the SwissTargetPrediction databases. A total of 450 epilepsy-related targets were identified from DisGeNET (Score-gda ≥ 0.05) and GeneCards (Relevance score ≥ 9.21)
Insert Fig. 5.
3.4. Protein-protein interaction (PPI) network analysis results
To obtain the PPI function network, the top targets were entered into the STRING data platform. Of the 55 targets that were correlated, 48 targets had average interactions (with a confidence score ≥ 0.4) in the PPI map generated by the STRING database. The targets of CA5B、CA5A、CA6、PTPRF、PTPN1 and ACP1 targets that were not bound to the lattice were dislodged, and the Cytoscape software was used to own the topology analysis to screen for significant central targets in the network. As shown in Fig. 6, the node size and color settings mirror the size, while the edge thickness and color settings reflect the combined score. When the color of the node closer to red, the corresponding value is higher. The results of the topological analysis displayed that the three targets with a higher degree value among the anti-epilepsy targets of neuronal acetylcholine receptor subunit α-4 (SLC6A4, degree = 25), neuronal acetylcholine receptor subunit α-4 (CHRNA4, degree = 22) and monoamine oxidase A (MAOA, degree = 21), which played a major role in the network.
Insert Fig. 6.
3.5. Molecular mechanism of anti-epilepsy the effect of PZE
We carried out a Metascape enrichment analysis of the 55 duplicate genes for GO and KEGG and achieved 583 biological processes (BP), 48 cellular ingredients (CC) and 68 molecular functions (MF) under the GO entry. The p-value determines the result of the correlation test [53], sort the top 20 of each category by - LogP value from large to small (Fig. 7). The top rank belonged to the vasoconstriction, regulation of tube diameter and blood vessel diameter maintenance enrichment in BP’s, synaptic membrane, dendrite and integral component of synaptic membrane enrichment in CC’s and G protein-coupled amine receptor activity, carbonate dehydratase activity, neurotransmitter receptor activity in MF’s.
The KEGG enrichment analysis of potential anti-epilepsy targets producted 32 signal pathways, which were ranked on the basis of the -log10p value from small to large, of which 20 pathways were relative to epilepsy (Fig. 8). Of the pathways related to epilepsy, Nitrogen metabolism has a high degree of enrichment, and can be considered as the main anti-epilepsy pathway. The significance of the count and the -log10p value represent a highly correlated way in which active constitutuents display their antiepileptic activities.
Insert Fig. 7. and Fig. 8.
3.6. Construction of the ingredient-target-pathway network
We adopt a network of 99 nodes and 240 edges by Cytoscape software (Fig. 6). In particular, the nodes consisted of 24 kinds of puffball constituents, 55 kinds of potential PZE targets, and 20 signal paths. It can be seen that each active component corresponds to multiple targets in multiple pathways (Fig. 9), which fully mirrors the multi-ingredient, multi-target and multi-pathway mechanisms of the puffball anti-epilepsy activity. With respect to the component attribute nodes in the network graph, it has been reported that M17 (Tetrahydrocolumbamine, degree value 12) and M15 (O-hydroxybenzoic acid, degree value 11) had the highest degree value; each of these components can be thought of as the central nodes in the network with a major anti-epilepsy effect.
Insert Fig. 9.
3.7. Molecular docking analysis
Under the prediction of molecular simulation docking methods, it can be seen that the key components screened by the pharmacology of the network have a good interplay with the predicted targets. They have a high binding energy and can successfully penetrate into the binding pocket of the target molecule in order to interact with key handicaps in the form of intermolecular hydrogen bonds. The cinnamic acids showed a good interaction with monoamine oxidase A (MAOA) in terms of their Structure-Activity Relationship, xanthones showed favorable interaction with nicotinic-acetylcholine receptor α4 (nAChR4), and saponins displayed preferable interaction with 5-HT transporter (SERT). Molecular construction and docking analysis between 2Z5Y and M16, M18, M20, and M37 were performed; between 6AWN and M17, M22, M25, M35; and between 5KXI and M18, M20, M31, M32-all of which are listed below. The lowest energy predicted complexes are stable intermolecular hydrogen bonding and stacking interactions. Interactions in the complex are decided by steric hindrance and the functional group category. We can see from Table 2 that the docking results are basically in agreement with the structure-function relationship, the substituted cinnamyl group on the phenolic hydroxyl group of vanillyl moiety helped enormously in promoting the binding affinity of the enzyme-inhibitor complexes, the ester linkage played key roles in the interaction with the protein. In terms of the interaction energy, the result displayed that compound M16 and M20 outperformed M18 and M37 in 5Z5Y, M22 and M25 outperformed M17 and M35 in 6AWN, M20 outperformed M18, M31 and M32 in 5KXI. Additionally, M20 displayed the best interaction with the protein 5Z5Y, which was partially different with the result of the enzyme inhibition assay, in general, active 5Z5Y residues comprising TYR 444, ILE 335, ILE 180 can interact with M16, M18, and M37 to a different extent (Fig. 10) [54–56]; in 6AWN that M25 showed the best interaction with the protein than M17, M22 and M25, the active residues of including ASP 98, SER 336, ALA 96, TYR 176, SER 438 can interact with M17, M22 and M25 in different degree (Fig. 10); in 5KXI that M20 demonstrated the best interaction with the protein than M18, M31 and M32, the active residues of including TRP 156, LEU 121, TYR 204, CYS 199 can interact with M18, M31 and M32 in different degree (Fig. 10).
Insert Fig. 10. and Fig. 11.
3.8. Anti-epilepsy Pathway
This study divulges that PZE activation of nAChRs can repair synaptic related protein levels. nAChRs belong to ligand gated ion channel coupled receptors, and form a centrally located ion channel, which primarily modulates Ca2+ plasma flow in and out of the cell. Calcium ions, acting as cytoplasmic second messengers, flow into cells via α4 nAChRs and connect with CaM to form an active CaM complex. The CaM complex may conjugate CaMK II, which regulates the spatial conformation of CaMK II and phosphorylates the CaMK II protein. Phosphorylated CaMK II enters the nucleus to catalyse the phosphorylation of Extracellular signal-regulated kinase (ERK) and cAMP responsive element binding protein (CREB). ERK and CREB are downstream of the α4 nAChR signal transduction pathway and serve as transcription factors that regulate neuron growth and development as well as the formation of long term memory (Fig. 12).
Insert Figue 12.
3.9. The efficacy and latency of PTZ
Convulsions were caused by an intraperitoneal (i.p.) injection of PTZ (10 mg/kg) dissolved in saline. Compared with PTZ group, after different degrees of PZE treatment, the latency time of PZEH, PZEM and PZEL was significantly prolonged in Fig. 13.
Insert Figue 13.
3.10. The behavioral animal tests
To assess anxiety-like and depression-like behaviors in PTZ-induced epilepsy mice using the EPMT, OFT, FST and TST. After 15 days of PZE intervention, the EPMT results showed that compared with NS group, although the number and timing of entry into the open arms of PTZ-induced epilepsy model mice was significantly lower than that of NS group mice. In comparison to the PTZ group, the number and d timing of mice entering the open arm was increased in the PZEH, PZEM, and PZEL groups (Fig. 14C&D), the anxiety and depression-like behaviors were improved, it indicated that PTZ could induce anxiety and depression-like behavior in mice. The OFT result indicated that compared with the NS group, there showed striking differences in the total movement distance (Fig. 15E), the number of entering enter central entry (Fig. 15F), the residence time in the central area (Fig. 15G) and the number of feces (Fig. 15H) in mice with epilepsy induced by PTZ, the intervention of PZEH and PZEM group has significant improvement effect on the above indexes. FST result manifested that compared with NS group, the time of immobility of PTZ group was prolonged and the depression was obvious, when given PZE, the immobility time of PZEH, PZEM, and PZEL group decreased, the difference was statistically remarkable (Fig. 16I, P < 0.05). The difference in TST resting time of mice in each group was statistically significant, the resting time of PTZ group was higher than that of NS group. After PZE treatment, the immobility time of PZEH, PZEM, and PZEL group decreased (Fig. 16J). These findings suggested that PZE can ameliorate the depressive response in mice and has antidepressant effects.
Insert Figue 14., Figue 15. and Figue 16.
3.11. Histological observation of brain tissue
3.11.1. Hematoxylin and eosin (HE) staining
HE staining was used to determine the lesion effect of hippocampal neuros in mice exposed to BPA. The cells of the hippocampus cornuammon1 (CA1), and cornu ammon 3 (CA3), dentate gyrus (DG), and Temporal cortex (TeA) of mice in the NS group were arranged in a regular and compact manner, as shown in Fig. 17 by light microscopy. In the NS group, many pyramidal cells appeared smaller with small condensed nuclei compared to the mice. In contrast, the hippocampal structure of mice in the PZEH and PZEM groups was essentially normal, and nerve cell degeneration in the large hippocampal region was less than in the PTZ group.
Insert Figue 17.
3.11.2. Immunofluorescence staining
To investigate the role of CHRNA4/CaMK II in Epilepsy, we first examined the expression of CHRNA4/CaMK II in the hippocampus tissues through immunofluorescence (Fig. 18). Figure 18 is showing that the protein of CHRNA4 and CaMK II were mainly expressed within TeA, CA1 and CA3 regions. Compared with the NS group, the proportion of CHRNA4 positive cells decreased significantly in PTZ group, the level of CHRNA4 in hippocampus increased rapidly in PZEH and PZEM group after PTZ-induced epilepsy, which is approximate to control (Fig. 18K). Oppositely, compared with the NS group, the proportion of CaMK II positive cells increased significantly in PTZ group, The level of CaMK II in hippocampus decreased rapidly in PZEH and PZEM group after PTZ-induced epilepsy, which is approximate to control (Fig. 18L). These experimental results in vitro and in vivo showed that PZE could reduce the epilepsy response and reduce the loss of neuronal cells by inhibiting the CHRNA4/CaMK II signal pathway.
Insert Figue 18
3.11.3. Western blot analysis
To explore the effect of PZE on the CHRNA4/CaMKП signaling pathway, Western blot was used to determine the protein levels of CHRNA4, CaMKП, ERK1/2, CREB1 in hippocampus tissues (Fig. 19). Further analysis of relative protein expression was performed using the optical density of the band. We found that the protein levels of CaMK II and CREB1 were significantly decreased in the PTZ group and that these increases were blunted by treatment with PZE (Fig. 19N & P, p < 0.01). Upon further comparison, it was found that the protein levels of CaMK II and CREB1 in the PZEH group were significantly higher than those in the PZEL group (p < 0.01). CHRNA4 and ERK1/2 expression levels in the PTZ group were both higher than those in the NS group, upon further comparison, we found that the protein levels of CHRNA4 and ERK1/2 in the PZEH group and in the PZEM group were strikingly lower than those in the PTZ group (Fig. 19M & O, p < 0.01). Together, these data suggested that PZE may protect against neuronal injury. The results of western blot analysis reveal that PZE can regulate the expression of the CHRNA4/CaMKП signaling pathway on epilepsy induced by PTZ in mice.
Insert Fig. 19.