Evidence-Based Complementary and Alternative Medicine
Volume 2019 (2019), Article ID 3104741, 7 pages
https://doi.org/10.1155/2019/3104741

Research Article

Uncovering the Pharmacological Mechanism of Chaibei Zhixian Decoction on Epilepsy by Network Pharmacology Analysis

Jian Zhang,¹ Chenglong Zheng,² Siyuan Yuan,¹ Xiaoke Dong,³ Le Wang,³ Yong Wang,⁴ Wei Wang,¹ Kuo Gao,³ and Jinmin Liu³

¹School of Traditional Chinese Medicine, Beijing University of Chinese Medicine, Beijing 100029, China
²Beijing Gulou Hospital of Traditional Chinese Medicine, Beijing, 10009, China
³Dongfang Hospital, Beijing University of Chinese Medicine, Beijing, 100078, China
⁴School of Life Science, Beijing University of Chinese Medicine, Beijing 100029, China

Correspondence should be addressed to Wei Wang; wangwei@bucm.edu.cn, Kuo Gao; linfengtingchan@foxmail.com, and Jinmin Liu; jmvip@vip.163.com

Received 20 February 2019; Revised 29 March 2019; Accepted 21 April 2019; Published 12 May 2019

Guest Editor: José C. T. Carvalho

Copyright © 2019 Jian Zhang 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.

Abstract.

Objective. Epilepsy is a neuronal disorder that is characterized by epileptic seizures and linked with abnormal neural functioning in the brain. Traditional Chinese medicine (TCM) formula Chaibei Zhixian decoction (CZD) has been widely used for epilepsy in China while the pharmacological mechanisms are still unclear. In the present study, systematic and comprehensive network pharmacology was utilized for the first time to reveal the potential pharmacological mechanisms of CZD on epilepsy. Methods. Traditional Chinese Medicine Systems Pharmacology (TCMSP) database and analysis platform was utilized for the development of an ingredients-targets database. After identifying epileptic targets of CZD, their interaction with other proteins was estimated based on protein-protein interaction network created from STITCH and gene ontology (GO) enrichment analysis utilizing Cytoscape-ClueGO plugin. Results. CZD formula was found to have 643 chemical ingredients, and the potential protein targets of these ingredients were 5230, as retrieved from TCMSP database. Twenty-six protein targets were found to be associated with epilepsy. Thirteen hub genes were regulated by CZD in epilepsy, including estradiol, ESR1, ESR2, SRC, CTNNB1, EP300, MAPK1, MAPK3, SP1, BRCA1, NCOA3, CHRM1, and GSK3B. The results of GO terms analysis showed that 8 GO terms were recovered in the form of 3 clusters, including negative regulation of protein kinase B signaling, positive regulation of interleukin-1 production, and microvillus assembly. Conclusions. Network pharmacology approach provides better understanding of the underlying pharmacological mechanisms of CZD on epilepsy. Estradiol, ESR1, ESR2, CTNNB1, EP300, MAPK1, MAPK3, BRCA1, and GSK3B are likely to be important molecules regulated by CZD in treatment of epilepsy. Negative regulation of protein kinase B signaling may play vital roles in the treatment of epilepsy by CZD.

1. Introduction

Epilepsy is a complex disorder involving neurological alterations that lead to the pathological development of recurrent seizures [1, 2]. Epilepsy affects millions of people worldwide and approximately one-third of patients suffer from cognitive impairment, particularly memory disruption [1, 3, 4]. First-line antiepileptic drugs have been given priority in the clinical treatment of epileptic seizures [1]. However, the risk of adverse effects from antiepileptic drugs is considerable and includes potential cognitive and behavioral effects [5]. Therefore, strategies that reduce the side effects of antiepileptic drugs or develop new drugs are urgently needed for epilepsy therapies.

Traditional Chinese medicine (TCM) has a long history in prevention and treatment of epilepsy in China [6, 7]. Chaibei Zhixian decoction (CZD), composed of Radix Bupleuri, Bulbus Fritillariae Thunbergii, Rhizoma Gastrodiae, Rhizoma Pinelliae, Rhizoma Acori Tatarinowii, Concha Ostreae, and Pheretima in a 4:3:5:3:3:10:2 ratio (Table 1), has been widely used in clinical treatment of epilepsy in China. Clinical study has shown that CZD is safe and effective for intractable epilepsy [8]. In addition, the combination of CZD with first-line antiepileptic drugs could reduce side effects and increase curative effects [9]. Some experimental studies have found CZD to have therapeutic effects on epilepsy by regulating multidrug resistance-associated protein 1, nuclear factor-kappa B, breast cancer resistance protein, and p-glycoprotein [10–13]. These studies all use traditional research method of single-drug, single-target, and single-pathway, but the TCM formula CZD has the characteristics of being multicomponent, multitarget, and multipathway. Thus, a new comprehensive and systematic evaluation of the pharmacological mechanism of CZD on epilepsy is critically needed.

Table 1: Pharmaceutical ingredients of Chaibei Zhixian decoction.

Latin name
Species
Family
Part used

Radix Bupleuri
Bupleurum Chinese DC.
Bupleurum scorzonerifolium Willd.
Umbelliferae
Roots

Bulbus Fritillariae Thunbergii
Fritillaria thunbergii Miq.
Liliaceae
Bulbs

Rhizoma Gastrodiae
Gastrodia elata Bl.
Orchidaceae
Rhizomes

Rhizoma Pinelliae
Pinellia ternata (Thunb.) Breit.
Araceae
Rhizomes

Rhizoma Acori Tatarinowii
Acorus tatarinowii Schott.
Araceae
Rhizomes

Concha Ostreae
Ostrea gigas Thunb.
Ostrea talienwhanensis Crosse
Ostrea rivularis Gould
Ostreidae
Concha

Pheretima
Pheretima aspergillum (E. Perrier)
Megascolecidae
Bodies

The ratio of these herbs was 4:3:5:3:3:10:2

Network pharmacology, including chemoinformatics, bioinformatics, network biology, and pharmacology, is a comprehensive method to uncover the bioactive components and potential mechanisms of TCM formulas from a systemic perspective [14]. In this study, the potential pharmacological mechanisms of CZD on epilepsy have been probed using network pharmacology, drug-target interaction databases, and a biological process analysis.

2. Methods

The first step of this study involved the retrieval of CZD constituents and their target proteins in Homo sapiens. Then the construction of CZD-target interaction network and its analysis was accomplished by using various GO terms. Finally, to assess the molecular mechanisms of CZD effects in epilepsy, Cytoscape along with its plugin ClueGO was utilized for GO enrichment analysis, followed by the analysis of biological processes.

2.1. Chemical Search and Their Target Retrieval

CZD contains Radix Bupleuri, Bulbus Fritillariae Thunbergii, Rhizoma Gastrodiae, Rhizoma Pinelliae, Rhizoma Acori Tatarinowii, Concha Ostreae, and Pheretima (Table 1). The chemical constituents present in these seven sources as well as the protein targets of these chemicals were searched via Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP, http://5th.tcmspw.com/tcmsp.php) [15]. The duplications in these chemical constituents from various sources and their targets were removed. From the retrieved targets, epilepsy-related targets were screened via Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.kegg.jp/) and employed for further analysis.

2.2. Conversion of Target Proteins into Network and Its Analysis

STITCH 5.0 database (http://stitch.embl.de/) [16] was utilized for analysis of the interaction among the identified protein targets to systematically investigate the mode of action of CZD. STITCH database has been furnished with extensive information regarding protein interactions. The main sources of information about these interactions are genomic model evaluations and high-throughput experimental outcomes. The information about > 9 million proteins from more than 2000 organisms has been added to this database. The mechanism of action of CZD and its essential pharmacodynamic constituents was assessed by developing a protein interaction network. The optional setting for network construction was set as follows: number of interactors = not more than 10; minimum required interaction score = 0.700 [16].

2.3. GO Terms Analysis through ClueGO Plugin

After identifying typical biological features of the protein targets, the ingredients of protein interaction network obtained from STITCH were used in ClueGO-based analysis, and GO enrichment analysis was introduced for the segmentation of target genes in a hierarchically arranged manner. ClueGO was utilized as a plugin of Cytoscape 3.4.0 software to construct [17], visualize, and evaluate protein target network and to study the biological pathways [18]. ClueGO analysis was conducted at a level of significance of 0.05. This study assumed that the network was of medium type. In contradiction to the detailed and global networks, medium network establishes GO terms belonging to GO levels 4-8, having a medium number of associated genes, and a medium percentage of uploaded genes found. Furthermore, two-sided hypergeometric test along with Bonferroni correction was used in network analysis. Lastly, organic layout algorithm was employed to visualize the functional network.

3. Results

3.1. Chemical Search and Their Target Retrieval

TCMSP search resulted in the retrieval of 643 chemical ingredients in the five herbs, Radix Bupleuri, Bulbus Fritillariae Thunbergii, Rhizoma Gastrodiae, Rhizoma Pinelliae, and Rhizoma Acori Tatarinowii, and two animals, Concha Ostreae and Pheretima. The protein targets of these 643 chemical ingredients recovered from TCMSP database were 5230 in number (Supplemental Table 1). Some studies revealed that CZD could be used for treating epilepsy [8–13]. After removing the repeated protein targets, 941 protein targets still remain (Supplemental Table 2). Among them, 26 protein targets of CZD were found to be associated with epilepsy through scanning TCMSP database (Supplemental Table 3), followed by their standardization through UniProt database mapping (http://www.uniprot.org/).

3.2. Conversion of Target Proteins into Network and Its Analysis

The systematically selected protein targets with a probabilistic confidence score of 0.700 were plotted as an interaction network (Figure 1) by using STITCH database (accessed in Oct 2018). The number of nodes and edges of this network were 21 and 48, respectively. Out of 21, there were 11 physical and 10 functional interactions. STITCH is a database of known and predicted interactions between chemicals and proteins. The interactions include direct (physical) and indirect (functional) associations; they stem from computational prediction, from knowledge transfer between organisms, and from interactions aggregated from other (primary) databases. The nodes and the edges represent protein/gene targets and their interactions, respectively. In case of random selection of nodes, the expected number of edges of the acquired PPIN was 23. The statistics of PPIN enrichment had a very small p-value (3.56E-06), indicating arbitrary nature of nodes and a significant number of edges.

Figure 1: PPIN (action view) showing CZD targets. The colored edges indicate the nature of action, as interpreted here: activation (), inhibition (), binding (), catalysis (), phenotype (), posttranslational modification (), reaction (), and transcriptional regulation (). The effects of action are represented by the following symbols: positive (), negative (), and unspecified (). AXIN1: axin 1; NCOA3: nuclear receptor coactivator 3; BRCA1: breast cancer 1, early onset; CHRM1: cholinergic receptor, muscarinic 1; CHRM5: cholinergic receptor, muscarinic 5; ESR1: estrogen receptor 1; CTNNB1: catenin (cadherin-associated protein), beta 1; CYP1A2: cytochrome P450, family 1, subfamily A, polypeptide 2; ESR2: estrogen receptor 2; HTR2A: 5-hydroxytryptamine (serotonin) receptor 2A; GSK3B: glycogen synthase kinase 3 beta; MAPK1: mitogen-activated protein kinase 1; MAPK3: mitogen-activated protein kinase 3; HTR2C: 5-hydroxytryptamine (serotonin); MAPT: microtubule-associated protein tau; SP1: Sp1 transcription factor; receptor 2C; NTF4: neurotrophin 4; EP300: E1A binding protein p300; NTRK2: neurotrophic tyrosine kinase, receptor, type 2; SLC6A4: solute carrier family 6 (neurotransmitter transporter, serotonin), member 4; SRC: v-src sarcoma (Schmidt-Ruppin A-2) viral oncogene homolog (avian).

The average node degree and clustering coefficient are the other important features of PPIN. Degree is a topological parameter that refers to the number of connections between a node and other nodes and can be used to describe the characteristics (particularly centrality) of nodes in the network. Node degree is a quantitative feature of a node that represents the number of linkages of a node in a network. A higher degree means a stronger correlation. The average node degree refers to the mean number of associations of a protein in a PPIN at threshold score, whereas the connectivity degree of PPIN nodes is indicated by the clustering coefficient. An increase in the clustering coefficient results in the increase in network connectivity. The average node degree and clustering coefficient values were 4.57 and 0.57, respectively. A node is termed as hub if it has number of linkage higher than its average node degree. Thirteen hubs have higher node degree than the average node degree, including estrogen receptor 1 (ESR1), estradiol, v-src sarcoma (Schmidt-Ruppin A-2) viral oncogene homolog (avian) (SRC), catenin beta 1 (CTNNB1), E1A binding protein p300 (EP300), estrogen receptor 1 (ESR2), mitogen-activated protein kinase 1 (MAPK1), Sp1 transcription factor (SP1), breast cancer 1, early onset (BRCA1), nuclear receptor coactivator 3 (NCOA3), cholinergic receptor, muscarinic 1 (CHRM1), glycogen synthase kinase 3 beta (GSK3B), and mitogen-activated protein kinase 3 (MAPK3) (Table 2). In addition, the functional proteins except NCOA3 and BRCA1 could be activated by CZD, which, on the other hand, inhibit CTNNB1, AXIN1, SRC, MAPT, SP1, and estradiol (Table 3). All functional proteins, as listed in Table 3, can bind with CZD. Moreover, the catalysis, posttranslational modifications, reactions, and expression are also found to be affected by CZD (Table 3).

Table 2: Node degree of the targets of CZD acquired via STITCH database.

Targets
Node Degree
Targets
Node Degree

ESR1
11
GSK3B
5

Estradiol
12
MAPK3
5

SRC
10
HTR2A
4

CTNNB1
8
MAPT
3

EP300
8
AXIN1
2

ESR2
8
CHRM5
2

MAPK1
8
HTR2C
2

SP1
8
NTRK2
2

BRCA1
7
SLC6A4
2

NCOA3
6
CYP1A2
1

CHRM1
5
NTF4
1

Table 3: Nature of action of functional targets of CZD acquired via STITCH.

Functional targets
Activation
Inhibition
Binding
Phenotype
Catalysis
Post-Trans. Mod.
Reaction
Expression
Score

CTNNB1
•
•
•

•
•
•

0.999

AXIN1
•
•
•

•
•
•

0.999

NCOA3

•

0.999

SRC
•
•
•

•
•
•
•
0.999

MAPT
•
•
•

•

•
0.999

SP1
•
•
•

•

•
0.999

BRCA1

•

0.999

NTF4
•

•

0.999

Estradiol
•
•
•

•

•
0.999

3.3. GO Terms through ClueGO Plugin

The annotation of biological functions was carried out by using GO terms and ClueGO plugin. This enrichment analysis of CZD targets resulted in the evolution of 8 GO terms which were ordered into 3 subgroups, including negative regulation of protein kinase B signaling, positive regulation of interleukin-1 production, and microvillus assembly (Table 4, Figure 2).

Table 4: Recovery of GO terms and the associated genes.

GO ID
GO Term
Term p Value (¤)
Group p Value (¤)
Associated Genes Found

51898
Negative regulation of protein kinase B signaling
750.0E-6 (4.5E-3)
750.0E-6 (750.0E-6)
PHLPP1, SLC9A3R1

32732
Positive regulation of interleukin-1 production
6.7E-6 (140.0E-6)
64.0E-6 (190.0E-6)
AZU1, HMGB1, TLR4

30033
Microvillus assembly
1.7E-6 (39.0E-6)
420.0E-6 (840.0E-6)
RAP1A, RAPGEF2, SLC9A3R1

with Bonferroni step down.

Figure 2: Targets involved in the biological effects. The most significant term in each stack is used to label the respective group. Node size is directly related to the term enrichment significance. The groups of GO terms having similar function are partially overlapped. [1-AIM2, 2-AZU1, 3-CALCA, 4-CARD8, 5-CASP1, 6-CASP5, 7-CCL19, 8-EGR1, 9-GSDMD, 10-HAVCR2, 11-HDAC2, 12-HMGB1, 13-HSPB1, 14-NOD1, 15-SMAD3, 16-TLR4, 17-ATP8B1, 18-EZR, 19-FSCN1, 20-FXYD5, 21-PLD1, 22-PRKCSH, 23-RAP1A, 24-RAPGEF2, 25-RAPGEF6, 26-SLC9A3R1, 27-PHLPP1].

4. Discussion

Epilepsy affects millions of people worldwide and approximately one-third of patients suffer from cognitive deficits. Due to the side effects of first-line antiepileptic drugs, more effective treatments are still needed. The TCM formula CZD not only is safe and effective for intractable epilepsy but also reduces side effects and increase curative effects when in combination with first-line antiepileptic drugs. However, the underlying mechanism of CZD on epilepsy is still unclear and remains unrevealed from a systemic point of view. Therefore, we adopted network pharmacology to further explore the mechanisms of CZD on epilepsy in this study. This systematic network pharmacology approach is a combination of various procedures, including retrieval of chemical ingredients of CZD, target search of these chemicals, development of network using these targets, and GO terms analysis. CZD formula was found to have 643 chemical ingredients, and the potential protein targets of these ingredients were 5230. Two aspects aroused our attention: first, 26 protein targets were found to be associated with epilepsy. Some of them are likely to be key molecules in the treatment of epilepsy with CZD. Second, GO terms analysis indicated that negative regulation of protein kinase B signaling, positive regulation of interleukin-1 production, and microvillus assembly have linkage with CZD treatment for epilepsy.

Network pharmacology analysis has shown that 13 hub genes were regulated by CZD in epilepsy, including estradiol, ESR1, ESR2, SRC, CTNNB1, EP300, MAPK1, MAPK3, SP1, BRCA1, NCOA3, CHRM1, and GSK3B. Among them, estradiol, ESR1, ESR2, CTNNB1, EP300, MAPK1, MAPK3, BRCA1, and GSK3B are closely related to epilepsy based on current studies. Thus, they are likely to be the main regulators of CZD in treatment of epilepsy. Estrogens affect neuronal excitability and have neuroprotective effects on seizure-induced hippocampal damage [19, 20]. Several studies have confirmed ESR was associated with epilepsy [21–24]. CTNNB1 has been implicated in epilepsy because of its altered postseizure expression [25, 26]. The dysfunction of CTNNB1-mediated signaling pathways leads to cortical malformation and increased seizure susceptibility [25]. EP300 may serve as potential targets for the treatment of epilepsy based on gene expression profile analysis of brain tissue of patients with epilepsy [27]. MAPK, as an important regulator of synaptic excitability, exerts an influence on epilepsy in animal models as well as human disease [27–30]. Interestingly, the component of CZD, gastrodin, has been reported to attenuate seizures by modulating the MAPK-associated inflammatory responses [31]. Variants in BRCA1-associated protein required for ATM activation-1 cause multifocal seizure syndrome [32, 33]. GSK3B activity protects neuronal networks from hyperactivation in response to epileptogenic stimuli [34]. Given that some main hubs of PPIN are closely related to epilepsy, they are likely to be important molecules regulated by CZD in treatment of epilepsy.

GO terms analysis revealed that negative regulation of protein kinase B signaling, positive regulation of interleukin-1 production, and microvillus assembly have linkage with CZD treatment for epilepsy. Protein kinase B, a serine/threonine-specific protein kinase, is involved in the regulation of binding phospholipids, phosphorylation, and ubiquitination but also modulates a wide array of cellular processes including cell apoptosis, metabolism, and proliferation [35]. Activation of the Akt signaling could alleviate neuronal apoptosis and oxidative stress [36]. The other two GO terms are positive regulation of interleukin-1 production and microvillus assembly. The former term represents a process in which the production of interleukin-1 is positively regulated, while the latter refers to the formation of microvillus. The positive regulation of interleukin-1 production has been found to be linked with three genes (i.e., AZU1, HMGB1, and TLR4). Microvillus assembly is regulated under the effect of three genes such as RAP1A, RAPGEF2, and SLC9A3R1. It has been reported that microvilli-like entities are linked with the inward movement of lethal capsulated neisseria meningitidis into vascular endothelial cells that may affect BBB resulting in seizures [37, 38].

In short, the present study has suggested various modes of CZD action against epilepsy, revealing that CZD profoundly enhances the performance of target genes involved in inhibiting epilepsy. The limitation of this study is that the bioactive components and targets found by network pharmacology analysis are the result of theoretical predictions and they should be verified by experiments. Further study will focus on using animal experiments and clinical trials to verify the hypothesis.

5. Conclusion

Network pharmacology analysis provides better understanding of the underlying pharmacological mechanisms of CZD on epilepsy. Our results revealed that estradiol, ESR1, ESR2, CTNNB1, EP300, MAPK1, MAPK3, BRCA1, and GSK3B are likely to be important molecules regulated by CZD in treatment of epilepsy. In addition, negative regulation of protein kinase B signaling may play vital roles in the treatment of epilepsy by CZD.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Kuo Gao, Jinmin Liu, and Wei Wang conceived and designed the project. Jian Zhang, Chenglong Zheng, Siyuan Yuan, Le Wang, and Xiaoke Dong implemented the methods and conducted the analysis. Kuo Gao, Jian Zhang, and Yong Wang drafted the manuscript. Jinmin Liu and Wei Wang revised the manuscript. All authors read and approved the final manuscript. Jian Zhang and Chenglong Zheng contributed equally to this work.

Acknowledgments

This study was financially supported by National Natural Science Foundation of China (nos. 81373586 and 81774277), Young Teacher Project of Beijing University of Chinese Medicine (2019-JYB-JS-095), and China Postdoctoral Science Foundation (2018M641286).

Supplementary Materials

Supplementary 1. Supplemental Table 1: 5230 protein targets from the TCMSP database.

Supplementary 2. Supplemental Table 2: nonrepetitive 941 protein targets from the TCMSP database.

Supplementary 3. Supplemental Table 3: 26 protein targets of CZD associated with epilepsy through scanning TCMSP database.

References

H. Yang, R. Zhang, C. Jia et al., “Neuronal protective effect of Songling Xuemaikang capsules alone and in combination with carbamazepine on epilepsy in kainic acid-kindled rats,” Pharmaceutical Biology, vol. 57, no. 1, pp. 22–28, 2019.

A. Bui, H. K. Kim, M. Maroso, and I. Soltesz, “Microcircuits in epilepsy: heterogeneity and hub cells in network synchronization,” Cold Spring Harbor Perspectives in Medicine, vol. 5, no. 11, 2015.

B. Bell, J. J. Lin, M. Seidenberg, and B. Hermann, “The neurobiology of cognitive disorders in temporal lobe epilepsy,” Nature Reviews Neurology, vol. 7, no. 3, pp. 154–164, 2011.

M. Maschio, L. Dinapoli, A. Fabi, D. Giannarelli, and T. Cantelmi, “Cognitive rehabilitation training in patients with brain tumor-related epilepsy and cognitive deficits: a pilot study,” Journal of Neuro-Oncology, vol. 125, no. 2, pp. 419–426, 2015.

G. Liu, N. Slater, and A. Perkins, “Epilepsy: treatment options,” American Family Physician, vol. 96, no. 2, pp. 87–96, 2017.

Z. Zhao, X. He, C. Ma et al., “Excavating anticonvulsant compounds from prescriptions of traditional chinese medicine in the treatment of epilepsy,” American Journal of Chinese Medicine, vol. 46, no. 4, pp. 707–737, 2018.

Y. Hijikata, A. Yasuhara, Y. Yoshida, and S. Sento, “Traditional Chinese medicine treatment of epilepsy,” The Journal of Alternative and Complementary Medicine, vol. 12, no. 7, pp. 673–677, 2006.

X. C. Zheng, S. F. Li, and J. M. Liu, “Clinical observation of single using Chaibei Zhixian decoction and combined medication on intractable epilepsy,” Tianjin Journal of Traditional Chinese Medicine, vol. 29, no. 03, pp. 224–227, 2012 (Chinese).

C. H. Wu and X. L. Liao, “Randomized parallel controlled study on complicated partial seizure of refractory epilepsy treated with chaibei zhixian decoction combined with western medicine,” Journal of Practical Traditional Chinese & Internal Medicine, vol. 32, no. 05, pp. 21–23, 2018 (Chinese).

X. H. Wang, Z. R. Yan, Q. Zhang et al., “Effects of different doses of Kangxianling Decoction on type I collagen gene expression in rats with doxorubicin-induced nephropathy,” China Journal of Traditional Chinese Medicine and Pharmacy, vol. 30, no. 06, pp. 2062–2065, 2015 (Chinese).

X. H. Wang, Z. R. Yan, Q. Zhang et al., “Effects of Chaibei Zhixian decoction on regulating the expression of P-gp/ Mdr1 in rat brain microvascular endothelial cells,” China Journal of Traditional Chinese Medicine and Pharmacy, vol. 31, no. 12, pp. 4961–4965, 2016 (Chinese).

Z. R. Yan, Q. Zhang, X. H. Wang et al., “Role of Chaibei Zhixian decoction on regulating the expression of BCRP and NF-κB p65 in rat brain microvascular endothelial cells,” Global Traditional Chinese Medicine, vol. 7, no. 04, pp. 241–246, 2014 (Chinese).

X. C. Zheng and J. M. Liu, “Effects of Chaibei Zhixian decoction on the expressions of P-glycoprotein and MDR1 in intractable epileptic rats,” Journal of Traditional Chinese Medicine, vol. 33, no. 11, pp. 4937–4940, 2018 (Chinese).

W. Song, S. Ni, Y. Fu, and Y. Wang, “Uncovering the mechanism of Maxing Ganshi Decoction on asthma from a systematic perspective: a network pharmacology study,” Scientific Reports, vol. 8, no. 1, Article ID 17362, 2018.

J. Ru, P. Li, J. Wang et al., “TCMSP: a database of systems pharmacology for drug discovery from herbal medicines,” Journal of Cheminformatics, vol. 6, no. 1, article 13, 2014.

M. Kuhn, D. Szklarczyk, S. Pletscher-Frankild et al., “STITCH 4: integration of protein-chemical interactions with user data,” Nucleic Acids Research, vol. 42, no. 1, pp. D401–D407, 2014.

G. Bindea, B. Mlecnik, H. Hackl et al., “ClueGO: a cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks,” Bioinformatics, vol. 25, no. 8, pp. 1091–1093, 2009.

X. Yao, H. Hao, Y. Li, and S. Li, “Modularity-based credible prediction of disease genes and detection of disease subtypes on the phenotype-gene heterogeneous network,” BMC Systems Biology, vol. 5, no. 1, article 79, 2011.

J. Velíšková, “Estrogens and epilepsy: why are we so excited?” Neuroscientist, vol. 13, no. 1, pp. 77–88, 2007.

D. A. Iacobas, S. Iacobas, N. Nebieridze, L. Velíšek, and J. Velíšková, “Estrogen protects neurotransmission transcriptome during status epilepticus,” Frontiers in Neuroscience, vol. 12, article 332, 2018.

C. F. Zorumski, S. M. Paul, Y. Izumi, D. F. Covey, and S. Mennerick, “Neurosteroids, stress and depression: potential therapeutic opportunities,” Neuroscience & Biobehavioral Reviews, vol. 37, no. 1, pp. 109–122, 2013.

L. Velíšek, N. Nebieridze, T. Chachua, and J. Velíšková, “Anti-seizure medications and Estradiol for Neuroprotection in epilepsy: The 2013 update,” Recent Patents on CNS Drug Discovery, vol. 8, no. 1, pp. 24–41, 2013.

H. E. Scharfman and N. J. MacLusky, “Sex differences in the neurobiology of epilepsy: a preclinical perspective,” Neurobiology of Disease, vol. 72, pp. 180–192, 2014.

D. H. Pham, C. C. Tan, C. C. Homan et al., “Protocadherin 19 (PCDH19) interacts with paraspeckle protein NONO to co-regulate gene expression with estrogen receptor alpha (ERα),” Human Molecular Genetics, vol. 26, no. 11, pp. 2042–2052, 2017.

V. E. Campos, M. Du, and Y. Li, “Increased seizure susceptibility and cortical malformation in β-catenin mutant mice,” Biochemical and Biophysical Research Communications, vol. 320, no. 2, pp. 606–614, 2004.

A. Rosiles, C. Rubio, C. Trejo, J. Gutierrez, L. Hernández, and C. Paz, “Participation of sox-1 expression and signaling of β-catenin in the pathophysiology of generalized seizures in cerebellum of rat,” CNS and Neurological Disorders - Drug Targets, vol. 15, no. 1, pp. 3–6, 2016.

D. Wang, X. Song, Y. Wang, X. Li, S. Jia, and Z. Wang, “Gene expression profile analysis in epilepsy by using the partial least squares method,” The Scientific World Journal, vol. 2014, Article ID 731091, 5 pages, 2014.

D.-M. Wu, Y.-T. Zhang, J. Lu, and Y.-L. Zheng, “Effects of microRNA-129 and its target gene c-Fos on proliferation and apoptosis of hippocampal neurons in rats with epilepsy via the MAPK signaling pathway,” Journal of Cellular Physiology, vol. 233, no. 9, pp. 6632–6643, 2018.

X. Zhu, Y. Chen, Y. Du, Q. Wan, Y. Xu, and J. Wu, “Astragaloside IV attenuates penicillin-induced epilepsy via inhibiting activation of the MAPK signaling pathway,” Molecular Medicine Reports, vol. 17, no. 1, pp. 643–647, 2017.

Z. Yang, J. Wang, C. Yu et al., “Inhibition of p38 MAPK signaling regulates the expression of EAAT2 in the brains of epileptic rats,” Frontiers in Neurology, vol. 9, article 925, 2018.

L. Chen, X. Liu, H. Wang, and M. Qu, “Gastrodin attenuates pentylenetetrazole-induced seizures by modulating the mitogen-activated protein kinase-associated inflammatory responses in mice,” Neuroscience Bulletin, vol. 33, no. 3, pp. 264–272, 2017.

S. Srivastava, H. E. Olson, J. S. Cohen et al., “BRAT1 mutations present with a spectrum of clinical severity,” American Journal of Medical Genetics Part A, vol. 170, no. 9, pp. 2265–2273, 2016.

S. A. Mundy, B. L. Krock, R. Mao, and J. J. Shen, “BRAT1-related disease-identification of a patient without early lethality,” American Journal of Medical Genetics Part A, vol. 170, no. 3, pp. 699–702, 2016.

M. Urbanska, P. Kazmierska-Grebowska, and T. Kowalczyk, “GSK3beta activity alleviates epileptogenesis and limits GluA1 phosphorylation,” EBioMedicine, vol. 39, pp. 377–387, 2019.

I. Hers, E. E. Vincent, and J. M. Tavaré, “Akt signalling in health and disease,” Cellular Signalling, vol. 23, no. 10, pp. 1515–1527, 2011.

A. Liu, M. Chu, and Y. Wang, “Up-regulation of trem 2 inhibits hippocampal neuronal apoptosis and alleviates oxidative stress in epilepsy via the PI3K/Akt pathway in mice,” Neuroscience Bulletin, 2019.

E. Eugène, I. Hoffmann, C. Pujol, P.-O. Couraud, S. Bourdoulous, and X. Nassif, “Microvilli-like structures are associated with the internalization of virulent capsulated Neisseria meningitidis into vascular endothelial cells,” Journal of Cell Science, vol. 115, no. 6, pp. 1231–1241, 2002.

M. R. Bower, M. Stead, R. S. Bower et al., “Evidence for consolidation of neuronal assemblies after seizures in humans,” The Journal of Neuroscience, vol. 35, no. 3, pp. 999–1010, 2015.