Network Pharmacology and Molecular Docking to Explore the Mechanism of Kangxian Decoction for Epilepsy

Abstract

Purpose. Kangxian decoction (KXD) has been used in clinical practice to treat epilepsy. The purpose of this study was to explore the active components of KXD and clarify its antiepileptic mechanism through network pharmacology and molecular docking. Methods. The components of KXD were collected from the Encyclopedia of Traditional Chinese Medicine (ETCM) database and the literature was searched. Then, active ingredients were screened by SwissADME and potential targets were predicted by the SwissTargetPrediction database. Epilepsy-related differentially expressed genes were downloaded from the Gene Expression Omnibus database. A component-target-pathway network was constructed with Cytoscape. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis and protein‒protein interaction network analysis revealed the potential mechanism and critical targets. Receiver operating characteristic (ROC) curves and box plots in microarray data validated the good diagnostic value and significant differential expression of these critical genes. Molecular docking verified the association between active ingredients and essential target proteins. Results. In our study, we screened the important compounds of KXD for epilepsy, including quercetin, baicalin, kaempferol, yohimbine, geissoschizine methyl ether, baicalein, etc. KXD may exert its therapeutic effect on epilepsy through the following targets: PTGS2, MMP9, CXCL8, ERBB2, and ARG1, acting on the following pathways: neuroactive ligand-receptor interactions, arachidonic acid metabolism, IL-17, TNF, NF-kappa B, and MAPK signaling pathways. The molecular docking results showed that the active ingredients in KXD exhibited good binding ability to the key targets. Conclusion. In this study, we explored the possibility that KXD for epilepsy may act on multiple targets through multiple active ingredients, involving neurotransmitters and neuroinflammatory pathways, providing a theoretical basis for subsequent clinical and experimental studies that will help develop effective new drugs to treat epilepsy.

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.

1. Introduction

Epilepsy is considered one of the most common serious brain disorders and affects more than 70 million people worldwide. It is not a disorder with a single presentation or cause, but a combination of symptoms with multiple risk factors and a strong genetic predisposition. Epilepsy has a bimodal distribution of prevalence, with infants and older populations at the highest risk [1]. Anxiety, depression, peptic ulcers, dementia, heart disease, and migraines are eight times more common in people with epilepsy than in the general population, leading to a high burden of comorbidities in people with epilepsy [2]. Approximately 80% of people with epilepsy live in low- and middle-income countries. The burden on people with epilepsy can be reduced through effective treatments [3]. However, although more than 20 antiepileptic drugs can be used for symptomatic treatment of seizures, approximately 1/3 of epileptic patients have seizures that are ineffective to drug therapy [4].

Traditional Chinese medicine (TCM) has been used for thousands of years as a unique treatment modality with the advantages of clinical efficacy, few side effects, and low drug resistance. In recent years, TCM research as a new treatment strategy has gained increased attention. Chinese medicine has a long history of treating epilepsy. Different herbs can exert antiepileptic effects through various mechanisms such as enhancing GABAergic effects, modulating sodium and NMDA channels, and exerting antioxidant, anti-inflammatory, and neuroprotective properties [5].

Kangxian decoction (KXD) is composed of six traditional Chinese medicines, including Acorus tatarinowii (Chinese name: Shichangpu), Pinellia ternata (Chinese name: Banxia), Panax notoginseng (Chinese name: Sanqi), Gastrodia elata (Chinese name: Tianma), Uncaria (Chinese name: Gouteng) and Pheretima (Chinese name: Dilong), which have the effect of eliminating phlegm and promoting blood circulation. The names of the plants involved are consistent with “The Plant List” (https://www.theplantlist.org/). Previous clinical and experimental studies have shown that KXD significantly prolongs the seizure incubation period, reduces seizure level, duration, and number of seizures in epileptogenic rats, and achieves antiapoptotic effects by reducing hippocampal Glu content, increasing GABA content, and promoting Bcl-2 protein expression [6–9]. In addition, traditional Chinese medicine relies on multicomponent, multitarget, and multipathway methods to treat complex diseases that are too complicated to analyze by conventional experimental methods. Therefore, we need innovative methods to study the pharmacological mechanism of KXD. As a new method, network pharmacology, combined with system network analysis and pharmacology, can clarify the synergistic effect and potential mechanism of the compound-target-disease network at the molecular level [10, 11]. This is consistent with the multicomponent, multitarget, and multipathway mechanism of action of TCM compounds and is in line with the principles of the holistic view and dialectical treatment in TCM. In particular, the application of network pharmacology in TCM provides researchers with a novel opportunity to understand TCM systematically, which may create new directions for studying potential pharmacological mechanisms and safety evaluation of TCM [12]. Molecular docking refers to the spatial docking of a small molecule into a macromolecule and obtaining complementary values for the binding site, which is used for structure-based drug design [13]. At present, although some herbal medicines are used to treat epilepsy, there is a lack of systematic and comprehensive observation of the pathways of drug intervention in the disease network. Therefore, this study investigated the molecular mechanism of action of KXD in the treatment of epilepsy using network pharmacology and molecular docking, which will provide a reference for future pharmacological and clinical studies and help develop new effective drugs for the treatment of epilepsy.

2. Materials and Methods

2.1. Screening of Active Ingredients and Potential Targets in KXD

The components in the KXD were collected from the Encyclopedia of Traditional Chinese Medicine (ETCM) database (https://www.tcmip.cn/ETCM/) [14]. To ensure the comprehensiveness and accuracy of the data, we also identified reported active ingredients from the PubMed, Web of Science, and China National Knowledge Infrastructure (CNKI) databases. The collected chemical components were individually entered into the PubChem database (https://pubchem.ncbi.nlm.nih.gov/) [15] to obtain the CAS number and canonical SMILES. Canonical SMILESs of the acquired compounds were put into the SwissADME platform (https://www.swissadme.ch/) [16]. The screening of the active ingredients was based on the following characteristics: (i) gastrointestinal absorption (GI absorption) was “High,” indicating good oral bioavailability and absorption of the ingredient; (ii) at least two of the five categories (Lipinski, Ghose, Veber, Egan, Muegge) were set to “Yes,” indicating that the compound has good drug-like properties [17]. Although some of the ingredients did not satisfy the above criteria, they were still included if their good pharmacological properties were confirmed through a literature review. Finally, the canonical SMILES of the collected active ingredients were entered into the SwissTargetPrediction database (https://www.swisstargetprediction.ch/) [18] to obtain the corresponding target genes. The criterion for target screening was a probability greater than 0.1.

2.2. Identification of Differentially Expressed Genes

The expression profile data of GSE143272 and the microarray platform GPL10558 were downloaded from the GEO database (https://www.ncbi.nlm.nih.gov/geo/). The expression profile data contained 34 untreated epilepsy patient samples (blood) and 51 normal samples (blood). Differentially expressed genes (DEGs) from the epilepsy samples and normal samples were screened using the limma package in R 4.0.2 software based on |logFC|>0.3, adj. $P$ value < 0.05. Visualization of DEGs was performed using the ggplot2 package as a volcano plot. From intersecting drug target genes with DEGs using the VennDiagram package, we eventually obtained KXD and epilepsy common associated genes.

2.3. Constructing Protein‒Protein Interaction (PPI) Networks and Screening Key Genes

KXD and epilepsy common-related genes were used to construct the PPI network using the STRING database (https://string-db.org/cgi/input.pl) [19]. The PPI network was obtained with default settings. The PPI network from the STRING database was then imported into Cytoscape 3.7.2 to research the critical genes. To ensure the accuracy of the results, we used two different plugins of Cytoscape to obtain two key subnetworks separately. First, using the MCODE plug-in in Cytoscape, the PPI network was further analyzed to obtain the potential module. Another way to retrieve the key subnetwork was to use the cytoHubba plug-in and apply the MCC algorithm to analyze the top 5 genes in the PPI network to obtain the second key subnetwork. Finally, we screened critical genes from the two subnetworks.

2.4. Receiver Operating Characteristic (ROC) Curves and Box Plots for Hub Genes

To assess the diagnostic ability of each hub gene for epilepsy, ROC curves were plotted, and the area under the curve (AUC) values were calculated using the pROC package in R software. The differential expression of hub genes between normal and epilepsy samples was assessed using box plots.

2.5. Pathway Enrichment Analysis

The cluster profiler package in R software was used to perform the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of genes commonly associated with KXD and epilepsy. $P<0.05$ was considered significant.

2.6. Constructing the Component-Target-Pathway Network

Cytoscape was used to build a network of relationships between active ingredients, common targets for KXD and epilepsy, and pathways. In the component-target-pathway network, “nodes” represent components, targets, and pathways, and “edges” represent the relationships between them. The “degree” parameter represents the number of connections between nodes in the network and is used to assess the important components and targets.

2.7. Molecular Docking

To validate the component-target association, the key target genes were selected for subsequent molecular docking analysis. The three-dimensional (3D) structure of the target protein was downloaded from the RCSB PDB database (https://www.rcsb.org/), and the two-dimensional structure of the molecular ligand was downloaded from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/). The proteins and molecules were preprocessed with PyMOL 2.4.0 software and saved in pdb format. The AutoDock Tools 1.5.6 hydrogenated and calculated the charges, selected the proteins as receptors and the molecular compounds as ligands, and then set up the docking parameters. The Grid Box was used to generate the docking area, which was defined as a 40 × 40 × 40 three-dimensional grid centered on the ligand binding site, with a grid space of 0.375 Å. Finally, the docking of the receptor protein and the molecular ligand was completed using AutoDock Vina 1.1.2, and the binding energy was obtained. The optimal conformation was selected and visualized using PyMOL. The entire network pharmacology and molecular docking process are illustrated in Figure 1.

[figure(s) omitted; refer to PDF]

3. Results

3.1. Screening of Active Ingredients and Potential Targets in KXD

After SwissADME screening and removing some compounds without targets, 49 active ingredients from six herbs were identified in the KXD formulation. Including 6 from Acorus tatarinowii (Chinese name: Shichangpu), 11 from Pinellia ternata (Chinese name: Banxia), 10 from Panax notoginseng (Chinese name: Sanqi), 8 from Gastrodia elata (Chinese name: Tianma), 18 from Uncaria (Chinese name: Gouteng), and 6 from Pheretima (Chinese name: Dilong). Notably, some of these ingredients were present in at least two herbs. Details of the active ingredients are shown in Table 1.

Table 1

Active ingredients of KXD.

Herb	ID	Molecule name	Molecular formula	CAS no.
BX: Pinellia ternata	BX1	Baicalin	C21H18O11	21967-41-9
	BX2	Baicalein	C22H20O11	491-67-8
	BX3	24-Ethylcholest-4-en-3-one	C29H48O	N/A
	BX4	Cavidine	C21H23NO4	32728-75-9
	BX5	Coniferin	C16H22O8	531-29-3
	BX6	Ephedrine	C10H15NO	299-42-3
	BX7	Coniine	C8H17N	458-88-8
DL: Pheretima	DL1	Succinic acid	C4H6O4	110-15-6
	DL2	Xanthine	C5H4N4O2	69-89-6
	DL3	Adenine	C5H5N5	73-24-5
	DL4	Valine	C5H11NO2	72-18-4
	DL5	Leucine	C6H13NO2	61-90-5
GT: Uncaria	GT1	Akuammicine	C20H22N2O2	639-43-0
	GT2	Angustidine	C19H15N3O	40217-50-3
	GT3	Corynantheine	C22H26N2O3	18904-54-6
	GT4	Corynoxeine	C22H26N2O4	630-94-4
	GT5	Dihydrocorynantheine	C22H28N2O3	50439-68-4
	GT6	Geissoschizine methyl ether	C22H26N2O3	60314-89-8
	GT7	Hirsuteine	C22H26N2O3	35467-43-7
	GT8	Isocorynoxeine	C22H26N2O4	51014-29-0
	GT9	Isopteropodine	C21H24N2O4	5171-37-9
	GT10	Isorhynchophylline	C21H26N2O4	144525-05-3
	GT11	Rhynchophylline	C22H28N2O4	76-66-4
	GT12	Vallesiachotamine	C21H22N2O3	5523-37-5
	GT13	Corynoxine	C22H28N2O4	6877-32-3
	GT14	Yohimbine	C21H26N2O3	146-48-5
SQ: Panax notoginseng	SQ1	Liquiritigenin	C15H12O4	578-86-9
	SQ2	Panaxynol	C17H24O	81203-57-8
	SQ3	Ginsenoside rh2	C36H62O8	78214-33-2
	SQ4	Dencichine	C5H8N2O5	5302-45-4
	SQ5	Methyl palmitate	C17H34O2	112-39-0
	SQ6	Panaxytriol	C17H26O3	N/A
SCP: Acorus tatarinowii	SCP1	8-Isopentenyl-kaempferol	C20H18O6	28610-31-3
	SCP2	Beta-asarone	C12H16O3	5273-86-9
	SCP3	Alpha-asarone	C12H16O3	2883-98-9
	SCP4	Methyl isoeugenol	C11H14O2	93-16-3
TM: Gastrodia elata	TM1	Vanillin	C8H8O3	121-33-5
	TM2	4-Ethoxymethylphenyl-4′-hydroxybenzylether	C16H18O3	77160-41-9
	TM3	Gastrodin	C13H18O7	62499-27-8
	TM4	4-(4′-hydroxybenzyloxy)benzyl methyl ether	C15H16O3	N/A
	TM5	Gastrodioside	C20H24O8	77162-64-2
	TM6	N6-(4-hydroxybenzyl)-adenosine	C17H19N5O5	110505-75-4
BX、GT、SQ、TM	D	Daucosterol	C35H60O6	474-58-8
BX、GT、SQ	B	Beta-sitosterol	C29H50O	83-46-5
BX、SQ	S	Stigmasterol	C29H48O	83-48-7
BX、SCP	C	Cycloartenol	C30H50O	469-38-5
DL、TM	P	Palmitic acid	C16H32O2	57-10-3
GT、SQ	Q	Quercetin	C15H10O7	117-39-5
GT、SCP	K	Kaempferol	C15H10O6	520-18-3

According to the target screening, there were 149 targets in Shichangpu, 212 targets in Banxia, 292 targets in Sanqi, 159 targets in Tianma, and 517 targets in Gouteng in KXD. After merging and removing duplicates, 705 target genes corresponding to the 49 active ingredients in KXD were finally identified.

3.2. Identification of Epilepsy-Related DEGs

We defined the filtering criteria as |logFC|>0.3, adj. $P$ value < 0.05. Finally, 312 differentially expressed genes associated with epilepsy were obtained, including 189 upregulated and 123 downregulated genes, and the results are presented in a volcano plot (Figure 2). For intersecting drug target genes with DEGs, we eventually obtained 19 genes commonly associated with KXD and epilepsy. The results are shown in a Venn diagram (Figure 3).

[figure(s) omitted; refer to PDF]

3.3. Construction of a Network of Component-Target-Pathway

To further understand the crossed genes, we performed a KEGG enrichment analysis. A total of 19 KEGG pathways were screened based on $p$ values < 0.05 (Figure 4). Details of the KEGG enrichment analysis can be found in Table 2. Major pathways included the sphingolipid signaling pathway, arachidonic acid metabolism, neuroactive ligand-receptor interaction, IL-17, TNF, NF-kappa B, and MAPK signaling pathways. We also identified 32 active ingredients targeting 19 intersecting targets from all active ingredients of KXD. To investigate the association between components, targets, and pathways, we constructed a component-target-pathway network (Figure 5). In general, a single ingredient can target multiple target genes, and a single target can be affected by multiple ingredients. The larger the area is, the larger the degree value, meaning the more important it is likely to be in the network. The network diagram showed that the mechanism of action of KXD has multicomponent, multitarget, and multipathway characteristics.

[figure(s) omitted; refer to PDF]

Table 2

Details of KEGG enrichment analysis.

ID	Description	$p$ value	Gene ID
hsa05219	Bladder cancer	9.18E-05	ERBB2 MMP9 CXCL8
hsa04657	IL-17 signaling pathway	0.001072	PTGS2 MMP9 CXCL8
hsa04071	Sphingolipid signaling pathway	0.002116	ADORA3 ASAH1 S1PR5
hsa05206	MicroRNAs in cancer	0.00413	CDC25B ERBB2 PTGS2 MMP9
hsa04913	Ovarian steroidogenesis	0.00552	PTGS2 AKR1C3
hsa04080	Neuroactive ligand-receptor interaction	0.007162	PTGER2 PTGDR ADORA3 S1PR5
hsa00590	Arachidonic acid metabolism	0.007821	PTGS2 AKR1C3
hsa05163	Human cytomegalovirus infection	0.012503	PTGER2 PTGS2 CXCL8
hsa05215	Prostate cancer	0.018985	ERBB2 MMP9
hsa01522	Endocrine resistance	0.019356	ERBB2 MMP9
hsa05146	Amoebiasis	0.020867	ARG1 CXCL8
hsa04064	NF-kappa B signaling pathway	0.021642	PTGS2 CXCL8
hsa04668	TNF signaling pathway	0.024856	PTGS2 MMP9
hsa04010	MAPK signaling pathway	0.025435	CDC25B ERBB2 CACNA2D2
hsa04670	Leukocyte transendothelial migration	0.025689	ITK MMP9
hsa04936	Alcoholic liver disease	0.038502	CPT1B CXCL8
hsa04921	Oxytocin signaling pathway	0.044611	CACNA2D2 PTGS2
hsa00220	Arginine biosynthesis	0.047584	ARG1
hsa05161	Hepatitis B	0.048873	MMP9 CXCL8

[figure(s) omitted; refer to PDF]

3.4. Identification of Hub Genes and Important Components

The 19 diseases and drug common target genes were imported into the STRING database, and tsv files were downloaded and imported into Cytoscape to construct PPI networks. Using the cytoHubba and MCODE plugins, we screened two identical subnetworks containing five hub genes, including MMP9, PTGS2, CXCL8, ERBB2, and ARG1 (Figures 6(a)–6(c)). We also identified 21 active ingredients targeting these hub genes through the component-target-pathway network. Based on literature reports of proven antiepileptic pharmacological effects, we screened six important components and validated them for subsequent molecular docking. These components were kaempferol, quercetin, baicalein, baicalin, yohimbine, and geissoschizine methyl ether.

[figure(s) omitted; refer to PDF]

3.5. Verification of Hub Genes Using ROC Curves and Box Plot Analysis

ROC curves were applied to evaluate the diagnostic value of the hub genes for epilepsy in the GSE143272 dataset. The AUC values of PTGS2, MMP9, CXCL8, ERBB2, and ARG1 were all greater than 0.65, indicating a good diagnostic value (Figure 7). Box plots indicated that all five hub genes were differentially expressed in the microarray dataset (Figure 8).

[figure(s) omitted; refer to PDF]

3.6. Molecular Docking Results

The five hub genes were selected for molecular docking analysis. Molecular docking was used to verify whether certain important compounds had a significant role in the regulation of the key target proteins. These important compounds were quercetin, baicalin, kaempferol, yohimbine, geissoschizine methyl ether, and baicalein. The interactions between the optimized conformations of certain important components and the key targets are shown visually in Figure 9. The docking scores are shown in Table 3. A binding energy of less than “−5” means that the compound and the target are well bound to each other [20]. The highest docking score was obtained for PTGS2 with baicalin at -10.1 kcal·mol⁻¹ and the lowest docking score was obtained for CXCL8 with quercetin at -5.3 kcal·mol⁻¹. Molecular docking showed that both these active compounds and the target proteins have good binding properties.

[figure(s) omitted; refer to PDF]

Table 3

Molecular docking score (kcal/mol).

Target name	PBD ID	Compounds
		Quercetin	Baicalin	Kaempferol	Yohimbine	Baicalein	Geissoschizine methyl ether
ARG1	3THJ	−7.9	−8.3	−7.7	−6.9	−7.9	−6.5
PTGS2	5IKT	−9.6	−10.1	−9.4	−9.2	−9.2	−7.3
MMP9	5TH6	−7.4	−8.1	−7.3	−8.1	−8.7	−7.4
CXCL8	1IKL	−5.3	−6.9	−5.6	−6.4	−5.8	−6.5
ERBB2	1MFL	−5.8	−6.6	−5.6	−5.7	−5.8	−5.8

4. Discussion

Pharmaco-resistance to antiepileptic drugs (AEDs) remains a primary unresolved therapeutic issue. Different antiepileptic drugs have been reported to cause serious side effects, with the development of pharmaco-resistance [21]. KXD is a Chinese herbal formula for the treatment of epilepsy, the effectiveness of which has been shown in our previous clinical research [6, 22, 23]. In our study, we screened the important compounds of KXD for epilepsy, including quercetin, baicalin, kaempferol, yohimbine, geissoschizine methyl ether, baicalein, etc.

Quercetin and kaempferol, both phytoflavonoids, are important dietary components of fruits and vegetables and have various biological effects, including antioxidant, anti-inflammatory, and neuroprotective activities. In vivo studies have shown that these phytoflavonoids can modulate the immune response, becoming an important treatment alternative for epilepsy [24]. GABA exists widely in the mammalian brain and is the main inhibitory neurotransmitter. Quercetin pretreatment exerted an inhibitory effect on GABAA receptor gene expression in a KA epilepsy model [25]. Silva et al. summarized the therapeutic effects of kaempferol in central nervous system disorders over the past ten years and found that kaempferol has significant anti-inflammatory pharmacological effects and exerts antioxidant effects through inhibition of MMP9 metalloproteinases, exhibiting multipotential neuroprotective effects [26]. The component-target-pathway network in this study also revealed that kaempferol may exert anti-inflammatory and neuroprotective effects by targeting MMP9. Yohimbine is a potent α2-adrenoceptor antagonist with anticonvulsant effects at low doses [27]. Geissoschizine methyl ether is a promising antiepileptic drug obtained from indole alkaloids isolated from Uncaria that can inhibit multiple neuronal channels [28]. Baicalein, as an important flavonoid, has various pharmacological effects, including anti-apoptotic, anti-inflammatory, antioxidant, and neuroprotective activities [29, 30]. Studies have demonstrated the neuroprotective effect of baicalin on PTZ-induced seizures [31].

The PPI network identified PTGS2, MMP9, CXCL8, ERBB2, and ARG1 as key genes. Inflammatory pathways play an essential role in the pathophysiology of epilepsy [32]. PTGS2, also known as COX2, is a major proinflammatory mediator associated with the severity and recurrence of epilepsy [33]. Epileptogenesis is associated with increased matrix metalloproteinases (MMPs) [34]. Tao et al. found that serum levels of MMP9 were higher in epileptic patients than in normal individuals and may serve as a biomarker to differentiate the cause of epilepsy [35]. This is consistent with the finding in this study that MMP9 was differentially expressed in the microarray dataset. Research has demonstrated that the CXCL8 (IL-8) signaling pathway can be activated during seizures [36]. ARG1 deficiency can lead to seizures, neurodevelopmental disorders, etc. [37]. Lamie et al. found that ERBB2 inhibitors could be a potential anticonvulsant drug [38]. These 5 key genes also demonstrated good diagnostic ability in the dataset using ROC curves and box plots.

Neuroinflammation-related pathways could be biomarkers and important therapeutic targets for epilepsy [39]. Glutamate and gamma-aminobutyric acid (GABA) are the main neurotransmitters that play key roles in the pathophysiology of epilepsy [40]. Numerous studies have shown that inflammatory responses and neurotransmitters in the brain play an important pathophysiological role in epilepsy [39, 40]. KEGG enrichment analysis enriched a total of 19 pathways involving neurotransmitters, including neuroactive ligand-receptor interactions and inflammatory responses, including arachidonic acid metabolism and the IL-17, TNF, MAPK, and NF-kappaB signaling pathways. Research has shown that abnormal interactions between neuroactive ligands and receptors can increase susceptibility to seizures [41]. Epileptogenesis is a complex pathological process involving mainly inflammatory responses and neuronal death [42]. The NF-κB pathway is associated with inflammation and oxidative stress, plays a key role in epileptogenesis, and can regulate neuronal survival and apoptosis [43]. The IL-17 signaling pathway and TNF signaling pathway play a critical role in inducing blood-brain barrier disruption and CNS autoinflammation, which are closely associated with seizures [44]. Cyclooxygenase (COX) converts arachidonic acid into prostaglandins, and activation of COX-2 is thought to be a trigger for neuroinflammation in the brain [45]. Finally, we used molecular docking to verify the association between active ingredients and key targets, and the results of molecular docking were consistent with the predicted interactions between targets and ingredients.

Systems pharmacology provided a robust approach to exploring the complex interactions between the ingredients and targets of KXD. In KXD, Acorus tatarinowii combined with Pinellia ternata dissolves turbidity and dispels phlegm, while Gastrodia elata and Uncaria pacify the liver and calm the wind, and Panax notoginseng combined with Pheretima has a meridian-clearing effect. This study revealed that different herbs have the same composition and that different components can act on the same target, allowing the individual herbs in KXD to better exert synergistic or antagonistic effects.

In this study, epilepsy-associated genes with differential expression were identified using microarray datasets, and the screened hub genes were validated with both ROC curves and box plots to enhance the accuracy of the prediction results. However, our study has several limitations. First, we need a more comprehensive database of TCM target genes to make the results of the network pharmacological analysis more reliable. Second, experimental validation and multidisciplinary collaboration are still needed to understand the exact mechanism of KXD in epilepsy.

5. Conclusion

In this study, we explored the multicomponent, multitarget, and multichannel characteristics of KXD-mediated epilepsy treatment through network pharmacology and molecular docking. The main components of KXD for epilepsy may be quercetin, baicalin, kaempferol, yohimbine, geissoschizine methyl ether, baicalein, etc. In conclusion, KXD may exert its therapeutic effect on epilepsy through the following targets: PTGS2, MMP9, CXCL8, ERBB2, and ARG1, acting on the following pathways: neuroactive ligand-receptor interactions, arachidonic acid metabolism, IL-17, TNF, NF-kappaB, and MAPK signaling pathways. This study preliminarily illustrated the mechanism of KXD in the treatment of epilepsy, providing a theoretical basis for further experimental validation and contributing to the development of effective drugs to treat epilepsy.

Disclosure

Weitao Wang, Yongquan Zhang, and Yibing Yang are the joint first authors of this paper.

Authors’ Contributions

Weitao Wang was responsible for writing the original manuscript. Yongquan Zhang and Yibing Yang collected and analyzed the data, and Lian Gu was responsible for supervising and revising the manuscript. Weitao Wang, Yongquan Zhang, and Yibing Yang represent equality of contribution to the work and should be considered as joint first authors.

Acknowledgments

The authors thank the authors for the development of the drug database, the sharing of microarray data, and the various analysis software. This study was funded by the National Natural Science Foundation of China (Grant no. 82160849), Guangxi Traditional Chinese Medicine Appropriate Technology Development and Promotion Project (Grant no. GZSY21-39), and Chongzuo City Science and Technology Plan Project (Grant no. Chongke FC2019009).

References

[1] R. D. Thijs, R. Surges, T. J. O'brien, J. W. Sander, "Epilepsy in adults," Lancet, vol. 393 no. 10172, pp. 689-701, DOI: 10.1016/S0140-6736(18)32596-0, 2019.

[2] M. R. Keezer, S. M. Sisodiya, J. W. Sander, "Comorbidities of epilepsy: current concepts and future perspectives," The Lancet Neurology, vol. 15 no. 1, pp. 106-115, DOI: 10.1016/s1474-4422(15)00225-2, 2016.

[3] E. Beghi, "The epidemiology of epilepsy," Neuroepidemiology, vol. 54 no. 2, pp. 185-191, DOI: 10.1159/000503831, 2020.

[4] W. Löscher, H. Potschka, S. M. Sisodiya, A. Vezzani, "Drug resistance in epilepsy: clinical impact, potential mechanisms, and new innovative treatment options," Pharmacological Reviews, vol. 72 no. 3, pp. 606-638, DOI: 10.1124/pr.120.019539, 2020.

[5] C. H. Lin, C. L. Hsieh, "Chinese herbal medicine for treating epilepsy," Frontiers in Neuroscience, vol. 15,DOI: 10.3389/fnins.2021.682821, 2021.

[6] Y. Zhang, W. Tan, H. Lu, "Clinical study on the treatment of idiopathic generalized tonic-clonic seizure type epilepsy with antiepileptic decoction," Liaoning Journal of Traditional Chinese Medicine, vol. 35 no. 8, pp. 1196-1197, 2008.

[7] Y. Zhang, G. Yang, "Effects of antiepileptic decoction on amino acid content in the brain of rats with experimental epilepsy," Lishizhen Medicine and Materia Medica Research, vol. 23 no. 3, pp. 672-673, 2012.

[8] W. Tan, Y. Zhang, H. Lu, "Effect of antiepileptic decoction on pentazocine induced epilepsy model in mice," Western journal of chinese medicine, vol. 24 no. 11, pp. 18-19, 2011.

[9] W. Tan, S. OuYang, Y. Zhang, "Effect of antiepileptic decoction on Bcl-2 protein expression in pentazocine epileptogenic rats," Journal of Emergency in Traditional Chinese Medicine, vol. 21 no. 8, pp. 1251-1253, 2012.

[10] X. Li, H. Yang, J. Xiao, J. Zhang, J. Zhang, M. Liu, Y. Zheng, L. Ma, "Network pharmacology based investigation into the bioactive compounds and molecular mechanisms of Schisandrae Chinensis Fructus against drug-induced liver injury," Bioorganic Chemistry, vol. 96,DOI: 10.1016/j.bioorg.2019.103553, 2020.

[11] T. Huang, Z. Ning, D. Hu, M. Zhang, L. Zhao, C. Lin, L. L. D. Zhong, Z. Yang, H. Xu, Z. Bian, "Uncovering the mechanisms of Chinese herbal medicine (MaZiRenWan) for functional constipation by focused network pharmacology approach," Frontiers in Pharmacology, vol. 9,DOI: 10.3389/fphar.2018.00270, 2018.

[12] T. Zhang, L. Pan, Y. Cao, N. Liu, W. Wei, H. Li, "Identifying the mechanisms and molecular targets of yizhiqingxin formula on alzheimer’s disease: coupling network pharmacology with GEO database," Pharmacogenomics and Personalized Medicine, vol. 13, pp. 487-502, DOI: 10.2147/pgpm.s269726, 2020.

[13] S. Saikia, M. Bordoloi, "Molecular docking: challenges, advances and its use in drug discovery perspective," Current Drug Targets, vol. 20 no. 5, pp. 501-521, DOI: 10.2174/1389450119666181022153016, 2019.

[14] H. Y. Xu, Y. Q. Zhang, Z. M. Liu, T. Chen, C. Y. Lv, S. H. Tang, X. B. Zhang, W. Zhang, Z. Y. Li, R. R. Zhou, H. J. Yang, X. J. Wang, L. Q. Huang, "ETCM: an encyclopaedia of traditional Chinese medicine," Nucleic Acids Research, vol. 47 no. D1, pp. D976-d982, DOI: 10.1093/nar/gky987, 2019.

[15] H. Dashti, J. R. Wedell, W. M. Westler, J. L. Markley, H. R. Eghbalnia, "Automated evaluation of consistency within the PubChem Compound database," Scientific Data, vol. 6 no. 1,DOI: 10.1038/sdata.2019.23, 2019.

[16] A. Daina, O. Michielin, V. Zoete, "SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules," Scientific Reports, vol. 7 no. 1,DOI: 10.1038/srep42717, 2017.

[17] H. Zhang, W. Dan, Q. He, J. Guo, S. Dai, X. Hui, P. Meng, Q. Cao, W. Yun, X. Guo, "Exploring the biological mechanism of huang yam in treating tumors and preventing antitumor drug-induced cardiotoxicity using network pharmacology and molecular docking technology," Evidence-based Complementary and Alternative Medicine, vol. 2021,DOI: 10.1155/2021/9988650, 2021.

[18] A. Daina, O. Michielin, V. Zoete, "SwissTargetPrediction: updated data and new features for efficient prediction of protein targets of small molecules," Nucleic Acids Research, vol. 47 no. W1, pp. W357-w364, DOI: 10.1093/nar/gkz382, 2019.

[19] D. Szklarczyk, A. L. Gable, D. Lyon, A. Junge, S. Wyder, J. Huerta-Cepas, M. Simonovic, N. T. Doncheva, J. H. Morris, P. Bork, L. J. Jensen, C. Mering, "STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets," Nucleic Acids Research, vol. 47 no. D1, pp. D607-d613, DOI: 10.1093/nar/gky1131, 2019.

[20] T. Gaillard, "Evaluation of AutoDock and AutoDock Vina on the CASF-2013 benchmark," Journal of Chemical Information and Modeling, vol. 58 no. 8, pp. 1697-1706, DOI: 10.1021/acs.jcim.8b00312, 2018.

[21] B. Chen, H. Choi, L. J. Hirsch, A. Katz, A. Legge, R. Buchsbaum, K. Detyniecki, "Psychiatric and behavioral side effects of antiepileptic drugs in adults with epilepsy," Epilepsy and Behavior, vol. 76, pp. 24-31, DOI: 10.1016/j.yebeh.2017.08.039, 2017.

[22] W. Tan, Y. Zhang, H. Lu, "Clinical study of antiepileptic decoction in the treatment of generalized tonic-clonic seizure type of symptomatic epilepsy," Western journal of chinese medicine, vol. 29 no. 03, pp. 76-78, 2016.

[23] W. Wang, Clinical Efficacy of Antiepileptic Decoction Combined with Citalopram in the Treatment of Epilepsy Combined with Anxiety and Depressive symptoms, 2016.

[24] H. Ahmed, M. A. Khan, S. A. Ali Zaidi, S. Muhammad, "In silico and in vivo: evaluating the therapeutic potential of kaempferol, quercetin, and catechin to treat chronic epilepsy in a rat model," Frontiers in Bioengineering and Biotechnology, vol. 9,DOI: 10.3389/fbioe.2021.754952, 2021.

[25] S. Moghbelinejad, Z. Rashvand, F. Khodabandehloo, G. Mohammadi, M. Nassiri-Asl, "Modulation of the expression of the GABAA receptor β 1 and β 3 subunits by pretreatment with quercetin in the KA model of epilepsy in mice: the effect of quercetin on GABAA receptor beta subunits," Journal of Pharmacopuncture, vol. 19 no. 2, pp. 163-166, DOI: 10.3831/KPI.2016.19.018, 2016.

[26] J. Silva Dos Santos, J. P. Gonçalves Cirino, P. De Oliveira Carvalho, M. M. Ortega, "The pharmacological action of kaempferol in central nervous system diseases: a review," Frontiers in Pharmacology, vol. 11,DOI: 10.3389/fphar.2020.565700, 2020.

[27] B. Payandemehr, A. Bahremand, A. Ebrahimi, S. E. Nasrabady, R. Rahimian, T. Bahremand, M. Sharifzadeh, A. R. Dehpour, "Protective effects of lithium chloride on seizure susceptibility: involvement of α 2-adrenoceptor," Pharmacology Biochemistry and Behavior, vol. 133, pp. 37-42, DOI: 10.1016/j.pbb.2015.03.016, 2015.

[28] Z. Q. Xie, X. T. Tian, Y. M. Zheng, L. Zhan, X. Q. Chen, X. M. Xin, C. G. Huang, Z. B. Gao, "Antiepileptic geissoschizine methyl ether is an inhibitor of multiple neuronal channels," Acta Pharmacologica Sinica, vol. 41 no. 5, pp. 629-637, DOI: 10.1038/s41401-019-0327-4, 2020.

[29] W. Z. Zhao, H. T. Wang, H. J. Huang, Y. L. Lo, A. M. Y. Lin, "Neuroprotective effects of baicalein on acrolein-induced neurotoxicity in the nigrostriatal dopaminergic system of rat brain," Molecular Neurobiology, vol. 55 no. 1, pp. 130-137, DOI: 10.1007/s12035-017-0725-x, 2018.

[30] S. Yang, H. Wang, Y. Yang, R. Wang, Y. Wang, C. Wu, G. Du, "Baicalein administered in the subacute phase ameliorates ischemia-reperfusion-induced brain injury by reducing neuroinflammation and neuronal damage," Biomedicine & Pharmacotherapy, vol. 117,DOI: 10.1016/j.biopha.2019.109102, 2019.

[31] J. Yang, Z. Jia, Z. Xiao, J. Zhao, Y. Lu, L. Chu, H. Shao, L. Pei, S. Zhang, Y. Chen, "Baicalin rescues cognitive dysfunction, mitigates neurodegeneration, and exerts anti-epileptic effects through activating TLR4/MYD88/caspase-3 pathway in rats," Drug Design, Development and Therapy, vol. 15, pp. 3163-3180, DOI: 10.2147/dddt.s314076, 2021.

[32] D. Boison, "New insights into the mechanisms of the ketogenic diet," Current Opinion in Neurology, vol. 30 no. 2, pp. 187-192, DOI: 10.1097/wco.0000000000000432, 2017.

[33] A. Dey, X. Kang, J. Qiu, Y. Du, J. Jiang, "Anti-inflammatory small molecules to treat seizures and epilepsy: from bench to bedside," Trends in Pharmacological Sciences, vol. 37 no. 6, pp. 463-484, DOI: 10.1016/j.tips.2016.03.001, 2016.

[34] D. W. Broekaart, A. Bertran, S. Jia, A. Korotkov, O. Senkov, A. Bongaarts, J. D. Mills, J. J. Anink, J. Seco, J. C. Baayen, S. Idema, E. Chabrol, A. J. Becker, W. J. Wadman, T. Tarrago, J. A. Gorter, E. Aronica, R. Prades, A. Dityatev, E. A. van Vliet, "The matrix metalloproteinase inhibitor IPR-179 has antiseizure and antiepileptogenic effects," The Journal of clinical investigation, vol. 131 no. 1,DOI: 10.1172/jci138332, 2021.

[35] H. Tao, Y. Gong, Q. Yu, H. Zhou, Y. Liu, "Elevated serum matrix metalloproteinase-9, interleukin-6, hypersensitive C-reactive protein, and homocysteine levels in patients with epilepsy," Journal of Interferon and Cytokine Research, vol. 40 no. 3, pp. 152-158, DOI: 10.1089/jir.2019.0137, 2020.

[36] R. Di Sapia, T. S. Zimmer, V. Kebede, S. Balosso, T. Ravizza, D. Sorrentino, M. A. M. Castillo, L. Porcu, F. Cattani, A. Ruocco, E. Aronica, M. Allegretti, L. Brandolini, A. Vezzani, "CXCL1-CXCR1/2 signaling is induced in human temporal lobe epilepsy and contributes to seizures in a murine model of acquired epilepsy," Neurobiology of Disease, vol. 158,DOI: 10.1016/j.nbd.2021.105468, 2021.

[37] J. Kido, S. Matsumoto, E. Takeshita, C. Hayakawa, K. Yamada, J. Kagawa, Y. Nakajima, T. Ito, H. Iijima, F. Endo, K. Nakamura, "Current status of surviving patients with arginase 1 deficiency in Japan," Molecular Genetics and Metabolism Reports, vol. 29,DOI: 10.1016/j.ymgmr.2021.100805, 2021.

[38] P. F. Lamie, A. M. El-Kalaawy, N. S. Abdel Latif, L. A. Rashed, J. N. Philoppes, "Pyrazolo[3,4-d]pyrimidine-based dual EGFR T790M/HER2 inhibitors: design, synthesis, structure-activity relationship and biological activity as potential antitumor and anticonvulsant agents," European Journal of Medicinal Chemistry, vol. 214,DOI: 10.1016/j.ejmech.2021.113222, 2021.

[39] A. Vezzani, S. Balosso, T. Ravizza, "Neuroinflammatory pathways as treatment targets and biomarkers in epilepsy," Nature Reviews Neurology, vol. 15 no. 8, pp. 459-472, DOI: 10.1038/s41582-019-0217-x, 2019.

[40] E. Akyuz, A. K. Polat, E. Eroglu, I. Kullu, E. Angelopoulou, Y. N. Paudel, "Revisiting the role of neurotransmitters in epilepsy: an updated review," Life Sciences, vol. 265,DOI: 10.1016/j.lfs.2020.118826, 2021.

[41] X. Liu, S. Ou, T. Xu, S. Liu, J. Yuan, H. Huang, L. Qin, H. Yang, L. Chen, X. Tan, Y. Chen, "New differentially expressed genes and differential DNA methylation underlying refractory epilepsy," Oncotarget, vol. 7 no. 52, pp. 87402-87416, DOI: 10.18632/oncotarget.13642, 2016.

[42] Z. Xiao, J. Peng, L. Yang, H. Kong, F. Yin, "Interleukin-1 β plays a role in the pathogenesis of mesial temporal lobe epilepsy through the PI3K/Akt/mTOR signaling pathway in hippocampal neurons," Journal of Neuroimmunology, vol. 282, pp. 110-117, DOI: 10.1016/j.jneuroim.2015.04.003, 2015.

[43] M. Cai, W. Lin, "The function of NF-kappa B during epilepsy, a potential therapeutic target," Frontiers in Neuroscience, vol. 16,DOI: 10.3389/fnins.2022.851394, 2022.

[44] G. Vitaliti, P. Pavone, S. Marino, M. A. N. Saporito, G. Corsello, R. Falsaperla, "Molecular mechanism involved in the pathogenesis of early-onset epileptic encephalopathy," Frontiers in Molecular Neuroscience, vol. 12,DOI: 10.3389/fnmol.2019.00118, 2019.

[45] X. Zhu, Y. Yao, J. Yang, J. Zhengxie, X. Li, S. Hu, A. Zhang, J. Dong, C. Zhang, G. Gan, "COX-2-PGE2 signaling pathway contributes to hippocampal neuronal injury and cognitive impairment in PTZ-kindled epilepsy mice," International Immunopharmacology, vol. 87,DOI: 10.1016/j.intimp.2020.106801, 2020.