# A Network Pharmacology to Explore the Mechanism ofCalculus Bovis in the Treatment of Ischemic Stroke

Authors: Fangchen Liu; Ling Li; Jian Chen; Ying Wu; Yongbing Cao; Ping Zhong
Journal: BioMed Research International
 (2021)

---

## Abstract

Background. Calculus Bovis is a valuable Chinese medicine, which is widely used in the clinical treatment of ischemic stroke. The present study is aimed at investigating its target and the mechanism involved in ischemic stroke treatment by network pharmacology. Methods. Effective compounds of Calculus Bovis were collected using methods of network pharmacology and using the Bioinformatics Analysis Tool for Molecular Mechanism of Traditional Chinese Medicine (BATMAN-TCM) and the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP). Potential compound targets were searched in the TCMSP and SwissTargetPrediction databases. Ischemic stroke-related disease targets were searched in the Drugbank, DisGeNet, OMIM, and TTD databases. These two types of targets were uploaded to the STRING database, and a network of their interaction (PPI) was built with its characteristics calculated, aiming to reveal a number of key targets. Hub genes were selected using a plug-in of the Cytoscape software, and Gene Ontology (GO) biological processes and pathway enrichment analyses of Kyoto Encyclopedia of Genes and Genomes (KEGG) were conducted using the clusterProfiler package of R language. Results. Among 12 compounds, deoxycorticosterone, methyl cholate, and biliverdin were potentially effective components. A total of 344 Calculus Bovis compound targets and 590 ischemic stroke targets were found with 92 overlapping targets, including hub genes such as TP53, AKT, PIK2CA, MAPK3, MMP9, and MMP2. Biological functions of Calculus Bovis are associated with protein hydrolyzation, phosphorylation of serine/threonine residues of protein substrates, peptide bond hydrolyzation of peptides and proteins, hydrolyzation of intracellular second messengers, antioxidation and reduction, RNA transcription, and other biological processes. Conclusion. Calculus Bovis may play a role in ischemic stroke by activating PI3K-AKT and MAPK signaling pathways, which are involved in regulating inflammatory response, cell apoptosis, and proliferation.


---

## Body

## 1. Introduction

Stroke is an acute cerebrovascular disease with typical clinical manifestations of sudden weakness in one side of the face, arms, or legs; sudden faintness; and unconsciousness. Ischemic stroke, the most common form of stroke, accounts for 70-80% of the total number of cases among stroke patients [1]. In China, deaths due to cerebrovascular diseases accounted for more than 20% of the total deaths in 2018 [2]. The rehabilitation of patients with ischemic stroke is often ineffective which brings a heavy burden to society and families. Currently, tissue-type plasminogen activator (tPA) is the only approved treatment for acute ischemic stroke [3–5]. However, its clinical application is greatly limited due to the narrow treatment time window, high bleeding risk, and many contraindications [6]. In China, stroke has been managed with herbs or other Chinese methods for thousands of years. Chinese herbal medicine are now widely accepted as the main complementary treatment in East Asia, North America, and Europe because of their good therapeutic effect, low toxicity, and low cost [7–9].Calculus Bovis, one of the most commonly used Chinese herbs for stroke, has been used for over 2,000 years in China. It was first described in “Shen Nong Ben Cao Jing” as a medication with a bitter taste and cooling nature [10]. And it has been applied in conditions like loss of consciousness due to stroke, epilepsy, mania, and other mental disorders. It was shown that Calculus Bovis protects the brain through its anti-inflammatory, antiapoptosis [11], antilipid peroxidation [12], and antioxidative stress effects [13]. It is well known that herbs have multiple ingredients targeting multiple sites and multiple pathways [14, 15]. Currently, Calculus Bovis and its formulas are widely used to treat ischemic stroke, but the mechanisms underlying its therapeutic effect have not been studied intensively.Network pharmacology for Chinese herbs is developed to decipher interactions between herbs and diseases at a system level by analysing the network between herbs, compounds, targets, diseases, and genes [16–19]. In the present study, our aim is to reveal the underlying mechanisms of Calculus Bovis in managing ischemic stroke by network pharmacology methods, which will lay the foundation for future pharmacological and clinical studies on ischemic stroke. The protocol of our experimental procedures is shown in Figure 1.Figure 1
The schematic map of the present study.


## 2. Materials and Methods

### 2.1. Data Acquisition

#### 2.1.1. Prediction of Compounds ofCalculus Bovis and Their Targets

Compounds ofCalculus Bovis were collected from the herbal platform TCMSP and BATMAN-TCM. The TCMSP (https://tcmspw.com/tcmsp.php) is a systems pharmacology platform for herbs providing information about compounds and their targets [20]. The BATMAN-TCM (http://bionet.ncpsb.org/batman-tcm/) is an online bioinformatics analysis tool comprised of functions like target prediction for herbs and target analysis [21]. When “NIU HUANG (Calculus Bovis)” was typed in the “Cluster name,” “Score cutoff” value was set at 20, and “Adjusted p value” was set at 0.05, compounds of Calculus Bovis and their targets would be displayed. In addition, potential targets could be searched in the TCMSP and SwissTargetPrediction (http://www.swisstargetprediction.ch/) databases [22] to further confirm the targets of compounds derived from Calculus Bovis. Names of target proteins were translated into gene names in the UniProt (http://www.uniprot.org/) database. If there was overlap in their target genes, the duplicates were deleted. Similarly, when the gene names of the protein targets were not found in the Uniprot database, they were deleted. SMILES IDs of compounds contained in Calculus Bovis were searched in the PubChem (https://pubchem.ncbi.nlm.nih.gov/) database, and their targets were predicted using SwissTargetPrediction after setting “Homo sapiens.” After collecting targets from the TCMSP, BATMAN-TCM, and SwissTargetPrediction databases, the duplicates were deleted.


#### 2.1.2. Prediction of Pharmacodynamics

In pharmacological studies, absorption, distribution, metabolism, and excretion (ADME) are key indices for identifying specific drugs [23]. Herein, 2 key parameters related to ADME, namely, oral bioavailability (OB) and drug-like activities (DL), were analyzed to explore potential bioactive compounds in Calculus Bovis. Based on the content of known compounds, OB and DL were set at ≥15% and ≥0.1, respectively. It has been reported that ergosterol (MOL000298; OB: 14.29%: DL: 0.72) is an indispensable compound of Calculus Bovis [24, 25] and was included in the present study. All compounds included in the present study were supported by the literature.


#### 2.1.3. Collection of Disease Targets of Ischemic Stroke

Key words such as “ischemic stroke,” “cerebral ischemic stroke,” and “brain ischemia” were used, and “Homo sapiens” was selected for species. They were searched in the Drugbank [26], DisGeNet [27], OMIM [28], and TTD databases [29], and duplicate genes were deleted.


#### 2.1.4. Venny Plotting

Both compound targets and disease targets were uploaded to the website of Venny 2.1.0 (https://bioinfogp.cnb.csic.es/tools/venny/); overlapping genes were the potential targets of bioactive compounds, and they interact with the body in ischemic stroke.


#### 2.1.5. Protein-Protein Interaction (PPI)

The STRING (https://string-db.org) database has collected a large number of well-known or predicted protein-protein interaction results [30]. Overlapping genes from Venny plots were uploaded to the database and “Homo sapiens” was selected for species, high confidence (0.700) was set for the minimum required interaction score, and irrelevant targets were concealed. As a result, a network map showing interactions between individual targets was rendered.


### 2.2. Network Construction and Hub Gene Selection

Network analysis facilitates interpretation of relationships between herbs, compounds, diseases, and genes. In the present study, two networks were constructed using Cytoscape 3.7.0 (https://cytoscape.org/) [31]: (1) a network of Calculus Bovis compounds and compound targets and (2) a PPI network of Calculus Bovis compound targets and a PPI network of Calculus Bovis treating ischemic stroke targets after connecting Calculus Bovis compound targets and disease targets. Topographical analysis for networks was completed using the NetworkAnalyzer tool in Cytoscape. Overlapping genes, also named Hub genes, were selected from the PPI networks using the cytoHubba plug-in in Cytoscape and 3 algorithms were used in the calculations. The latter included degree of freedom, Maximum Neighborhood Component (MNC), and Maximal Clique Centrality (MCC).


### 2.3. Gene Functions and Pathway Enrichment Analysis

Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted using the clusterProfiler package of R (R 4.0.2 for Windows) to identify biological processes and pathways with systemic involvement. Biological processes and pathways with a significant difference were selected, and their numbers of enrichment as well their asp values were ranked. The top 20 results from the GO enrichment and KEGG pathway enrichment analyses were presented. Visualization of these pathways with a p value < 0.05 was completed using the R software package.


## 2.1. Data Acquisition

### 2.1.1. Prediction of Compounds ofCalculus Bovis and Their Targets

Compounds ofCalculus Bovis were collected from the herbal platform TCMSP and BATMAN-TCM. The TCMSP (https://tcmspw.com/tcmsp.php) is a systems pharmacology platform for herbs providing information about compounds and their targets [20]. The BATMAN-TCM (http://bionet.ncpsb.org/batman-tcm/) is an online bioinformatics analysis tool comprised of functions like target prediction for herbs and target analysis [21]. When “NIU HUANG (Calculus Bovis)” was typed in the “Cluster name,” “Score cutoff” value was set at 20, and “Adjusted p value” was set at 0.05, compounds of Calculus Bovis and their targets would be displayed. In addition, potential targets could be searched in the TCMSP and SwissTargetPrediction (http://www.swisstargetprediction.ch/) databases [22] to further confirm the targets of compounds derived from Calculus Bovis. Names of target proteins were translated into gene names in the UniProt (http://www.uniprot.org/) database. If there was overlap in their target genes, the duplicates were deleted. Similarly, when the gene names of the protein targets were not found in the Uniprot database, they were deleted. SMILES IDs of compounds contained in Calculus Bovis were searched in the PubChem (https://pubchem.ncbi.nlm.nih.gov/) database, and their targets were predicted using SwissTargetPrediction after setting “Homo sapiens.” After collecting targets from the TCMSP, BATMAN-TCM, and SwissTargetPrediction databases, the duplicates were deleted.


### 2.1.2. Prediction of Pharmacodynamics

In pharmacological studies, absorption, distribution, metabolism, and excretion (ADME) are key indices for identifying specific drugs [23]. Herein, 2 key parameters related to ADME, namely, oral bioavailability (OB) and drug-like activities (DL), were analyzed to explore potential bioactive compounds in Calculus Bovis. Based on the content of known compounds, OB and DL were set at ≥15% and ≥0.1, respectively. It has been reported that ergosterol (MOL000298; OB: 14.29%: DL: 0.72) is an indispensable compound of Calculus Bovis [24, 25] and was included in the present study. All compounds included in the present study were supported by the literature.


### 2.1.3. Collection of Disease Targets of Ischemic Stroke

Key words such as “ischemic stroke,” “cerebral ischemic stroke,” and “brain ischemia” were used, and “Homo sapiens” was selected for species. They were searched in the Drugbank [26], DisGeNet [27], OMIM [28], and TTD databases [29], and duplicate genes were deleted.


### 2.1.4. Venny Plotting

Both compound targets and disease targets were uploaded to the website of Venny 2.1.0 (https://bioinfogp.cnb.csic.es/tools/venny/); overlapping genes were the potential targets of bioactive compounds, and they interact with the body in ischemic stroke.


### 2.1.5. Protein-Protein Interaction (PPI)

The STRING (https://string-db.org) database has collected a large number of well-known or predicted protein-protein interaction results [30]. Overlapping genes from Venny plots were uploaded to the database and “Homo sapiens” was selected for species, high confidence (0.700) was set for the minimum required interaction score, and irrelevant targets were concealed. As a result, a network map showing interactions between individual targets was rendered.


## 2.1.1. Prediction of Compounds ofCalculus Bovis and Their Targets

Compounds ofCalculus Bovis were collected from the herbal platform TCMSP and BATMAN-TCM. The TCMSP (https://tcmspw.com/tcmsp.php) is a systems pharmacology platform for herbs providing information about compounds and their targets [20]. The BATMAN-TCM (http://bionet.ncpsb.org/batman-tcm/) is an online bioinformatics analysis tool comprised of functions like target prediction for herbs and target analysis [21]. When “NIU HUANG (Calculus Bovis)” was typed in the “Cluster name,” “Score cutoff” value was set at 20, and “Adjusted p value” was set at 0.05, compounds of Calculus Bovis and their targets would be displayed. In addition, potential targets could be searched in the TCMSP and SwissTargetPrediction (http://www.swisstargetprediction.ch/) databases [22] to further confirm the targets of compounds derived from Calculus Bovis. Names of target proteins were translated into gene names in the UniProt (http://www.uniprot.org/) database. If there was overlap in their target genes, the duplicates were deleted. Similarly, when the gene names of the protein targets were not found in the Uniprot database, they were deleted. SMILES IDs of compounds contained in Calculus Bovis were searched in the PubChem (https://pubchem.ncbi.nlm.nih.gov/) database, and their targets were predicted using SwissTargetPrediction after setting “Homo sapiens.” After collecting targets from the TCMSP, BATMAN-TCM, and SwissTargetPrediction databases, the duplicates were deleted.


## 2.1.2. Prediction of Pharmacodynamics

In pharmacological studies, absorption, distribution, metabolism, and excretion (ADME) are key indices for identifying specific drugs [23]. Herein, 2 key parameters related to ADME, namely, oral bioavailability (OB) and drug-like activities (DL), were analyzed to explore potential bioactive compounds in Calculus Bovis. Based on the content of known compounds, OB and DL were set at ≥15% and ≥0.1, respectively. It has been reported that ergosterol (MOL000298; OB: 14.29%: DL: 0.72) is an indispensable compound of Calculus Bovis [24, 25] and was included in the present study. All compounds included in the present study were supported by the literature.


## 2.1.3. Collection of Disease Targets of Ischemic Stroke

Key words such as “ischemic stroke,” “cerebral ischemic stroke,” and “brain ischemia” were used, and “Homo sapiens” was selected for species. They were searched in the Drugbank [26], DisGeNet [27], OMIM [28], and TTD databases [29], and duplicate genes were deleted.


## 2.1.4. Venny Plotting

Both compound targets and disease targets were uploaded to the website of Venny 2.1.0 (https://bioinfogp.cnb.csic.es/tools/venny/); overlapping genes were the potential targets of bioactive compounds, and they interact with the body in ischemic stroke.


## 2.1.5. Protein-Protein Interaction (PPI)

The STRING (https://string-db.org) database has collected a large number of well-known or predicted protein-protein interaction results [30]. Overlapping genes from Venny plots were uploaded to the database and “Homo sapiens” was selected for species, high confidence (0.700) was set for the minimum required interaction score, and irrelevant targets were concealed. As a result, a network map showing interactions between individual targets was rendered.


## 2.2. Network Construction and Hub Gene Selection

Network analysis facilitates interpretation of relationships between herbs, compounds, diseases, and genes. In the present study, two networks were constructed using Cytoscape 3.7.0 (https://cytoscape.org/) [31]: (1) a network of Calculus Bovis compounds and compound targets and (2) a PPI network of Calculus Bovis compound targets and a PPI network of Calculus Bovis treating ischemic stroke targets after connecting Calculus Bovis compound targets and disease targets. Topographical analysis for networks was completed using the NetworkAnalyzer tool in Cytoscape. Overlapping genes, also named Hub genes, were selected from the PPI networks using the cytoHubba plug-in in Cytoscape and 3 algorithms were used in the calculations. The latter included degree of freedom, Maximum Neighborhood Component (MNC), and Maximal Clique Centrality (MCC).


## 2.3. Gene Functions and Pathway Enrichment Analysis

Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted using the clusterProfiler package of R (R 4.0.2 for Windows) to identify biological processes and pathways with systemic involvement. Biological processes and pathways with a significant difference were selected, and their numbers of enrichment as well their asp values were ranked. The top 20 results from the GO enrichment and KEGG pathway enrichment analyses were presented. Visualization of these pathways with a p value < 0.05 was completed using the R software package.


## 3. Results and Analysis

A total of 12 bioactive compounds were found after ADME screening, and all of them have been verified in other studies.


### 3.1.Calculus Bovis Compound-Target Network

A total of 344 targets were found by 12 bioactive compounds ofCalculus Bovis. Details of these ingredients are listed in Table 1, and the map of the compound-target network is shown in Figure 2. Circles represent compounds of Calculus Bovis, and triangles represent the targets; their colors are darkened as their degrees of freedom increased. There were 362 nodes and 706 edges with a network density of 0.011 and a network diameter of 5. Details of these parameters are listed in Table 2.Table 1
Active ingredients ofCalculus Bovis.
Mol IDMol name2D structureOB (%)DLMOL000263Oleanolic acid29.020.76MOL000298Ergosterol14.290.72MOL000511Ursolic acid16.770.75MOL008834Bilicerdin23.790.75MOL0088353-[2-[[3-(2-Carboxyethyl)-4-methyl-5-[(E)-(4-methyl-5-oxo-3-vinyl-2-pyrrolylidene)methyl]-1H-pyrrol-2-yl]methyl]-4-methyl-5-[(Z)-(3-methyl-5-oxo-4-vinyl-2-pyrrolylidene)methyl]-1H-pyrrol-3-yl]propanoic acid16.530.75MOL008838Methyl (4R)-4-[(3R,5S,7S,8R,9S,10S,12S,13R,14S,17R)-3,7,12-trihydroxy-10,13-dimethyl-2,3,4,5,6,7,8,9,11,12,14,15,16,17-tetradecahydro-1H-cyclopenta[a]phenanthren-17-yl]pentanoate32.320.76MOL008839Methyl desoxycholate34.630.73MOL0088402-[(3alpha,12alpha-Dihydroxy-24-oxo-5beta-cholan-24-yl)amino]ethanesulfonic acid15.920.87MOL008843Cherianoine27.320.12MOL008846ZINC0128036546.380.49MOL000953CLR37.870.68MOL009807CHD22.170.72Abbreviations: OB = oral bioavailability; DL = drug likeness; Mol = molecular.Figure 2
Calculus Bovis compound-target network. Note: circles represent compounds, triangles represent targets, and their colors darken as their degrees of freedom increase.Table 2
Calculus Bovis compound-candidate target network parameters.
Network parametersValuesNumber of nodes362Network density0.011Network diameter5Network heterogeneity3.212Average number of neighbors3.901Characteristic path length3.236Shortest paths130682 (100%)Network centralization0.304


### 3.2. Disease Targets of Ischemic Stroke

Using key words listed in Section2.1.3, 74 disease targets were found in the Drugbank database, 313 were found in the DisGeNet database after screening for values above the average, 322 were found in the OMIM database, and 10 were found in the TTD database. A total of 590 targets were found after deleting gene duplicates.


### 3.3. Prediction ofCalculus Bovis Targets in Ischemic Stroke

After entering compound targets and disease targets to Venny 2.1.0 in the form of gene names, 92 overlapping genes were found (see Figure3). These genes were the shared targets of the bioactive compounds and disease targets in ischemic stroke.Figure 3
Venny diagram ofCalculus Bovis targets and ischemic stroke disease targets. CB: Calculus Bovis target; IS: ischemic stroke disease target.


### 3.4. Construction of PPI Networks

#### 3.4.1. The PPI Network ofCalculus Bovis Targets

PPI networks have been widely used in studying protein-protein interactions in different diseases. To construct the PPI network ofCalculus Bovis targets, 12 compounds were connected with 344 targets in the TCMSP and the SwissTargetPrediction databases, and they were imported into the STRING database (species: Homo sapiens; minimum required interaction score: high confidence (0.700)), and the PPI network was visualized after reconstructing it with Cytoscape (version 3.7.0). As shown in Figure 4, colors darkened as the degree of freedom increased. This PPI network contained 322 nodes and 2195 edges with a diameter of 8 and an average number of neighbors of 13.634. It showed that TP53 (degree=82), AKT1 (degree=67), MAPK1 (degree = 67), PIK3CA (degree=63), SRC (degree=58), MAPK3 (degree=56), VEGFA (degree=52), HSP90AA1 (degree=48), JUN (degree=46), EGFR (degree=46), MAPK8 (degree=45), TNF (degree=40), and CASP3 (degree=40) were the key nodes of this PPI network.Figure 4
Protein-protein interaction (PPI) network ofCalculus Bovis compound targets. Darkness of colors and sizes of circles represent the degree of freedom as the scale indicates.


#### 3.4.2. The PPI Network ofCalculus Bovis-Ischemic Stroke Targets and Hub Genes

To explore the potential therapeutic mechanisms ofCalculus Bovis in managing ischemic stroke, 92 shared targets by Calculus Bovis compounds and ischemic stroke were connected and imported into the STRING database as shown in Figure 5(a). The PPI network of Calculus Bovis-ischemic stroke targets were constructed by visualizing the results using Cytoscape (Figure 5(b)). This PPI network had 83 nodes and 403 edges with a network diameter of 6 and an average number of neighbors of 9.711. The top 10 targets with the highest degree of freedom were TP53 (degree=33), AKT1 (degree=30), MAPK1 (degree=29), VEGFA (degree=27), TNF (degree=25), PIK3CA (degree=25), MAPK3 (degree=24), MMP9 (degree=22), PTGS2 (degree=22), and IL-1B (degree=20). A total of 6 Hub genes were found using 3 algorithms: degree of freedom, Maximum Neighborhood Component (MNC), and Maximal Clique Centrality (MCC) (Figure 5(c)). These genes were TP53, AKT1, PIK3CA, MAPK3, MMP9, and MMP2. Their network of interactions is shown in Figure 5(d). In addition, TP53, AKT1, MAPK1, VEGFA, PIK3CA, and MAPK3 were among the top 10 candidates ranked by degree of freedom.Figure 5
The PPI network ofCalculus Bovis compound targets against ischemic stroke (IS): (a) the original PPI data generated from the STRING database showing detailed interactions of the targets; (b) the PPI network constructed using Cytoscape (version 3.7.0); (c) the top 15 genes were calculated from the PPI network by the degree of freedom, MNC, and MCC, and the overlapping genes were then screened by Venn diagrams; (d) the PPI network of hub genes.
(a)(b)(c)(d)


### 3.5. GO Enrichment Analysis

To further investigate the underlying mechanism ofCalculus Bovis in managing ischemic stroke, 344 targets and 92 shared targets were collected for GO enrichment analysis.


#### 3.5.1. GO Enrichment Analysis forCalculus Bovis Targets

Twenty results were selected based on theirp values and their numbers of enrichment. They were primarily involved in ion channel activity (GO:0005216), channel activity (GO:0015267), passive transmembrane transporter activity (GO:0022803), gated channel activity (GO:0022836), cation channel activity (GO:0005261), metal ion transmembrane transporter activity (GO:0046873), neurotransmitter receptor activity (GO:0030594), ligand-gated ion channel activity (GO:0015276), ligand-gated channel activity (GO:0022834), postsynaptic neurotransmitter receptor activity (GO:0098960), extracellular ligand-gated ion channel activity (GO:0005230), neurotransmitter receptor activity involved in regulation of postsynaptic membrane potential (GO:0099529), transmitter-gated ion channel activity (GO:0022824), transmitter-gated channel activity (GO:0022835), and steroid dehydrogenase activity (GO:0016229). Details are listed in Table 3. Results are presented using a bubble plot and a column chart using the R package.Table 3
GO enrichment analysis results ofCalculus Bovis targets.
No.IDDescriptionCountp value1GO:0005216Ion channel activity141.15918×10−102GO:0015267Channel activity143.57104×10−103GO:0022803Passive transmembrane transporter activity143.67058×10−104GO:0022836Gated channel activity137.88226×10−115GO:0005261Cation channel activity111.25701×10−86GO:0046873Metal ion transmembrane transporter activity111.85161×10−77GO:0030594Neurotransmitter receptor activity91.68357×10−108GO:0015276Ligand-gated ion channel activity81.76915×10−89GO:0022834Ligand-gated channel activity81.76915×10−810GO:0098960Postsynaptic neurotransmitter receptor activity73.22861×10−911GO:0005230Extracellular ligand-gated ion channel activity61.78299×10−712GO:0099529Neurotransmitter receptor activity involved in regulation of postsynaptic membrane potential56.54501×10−713GO:0022824Transmitter-gated ion channel activity51.78703×10−614GO:0022835Transmitter-gated channel activity51.78703×10−615GO:0016229Steroid dehydrogenase activity45.40×10−616GO:0016628Oxidoreductase activity, acting on the CH-CH group of donors, NAD or NADP as acceptor41.58×10−617GO:1904315Transmitter-gated ion channel activity involved in regulation of postsynaptic membrane potential41.78589×10−518GO:0022851GABA-gated chloride ion channel activity39.83×10−619GO:0042166Acetylcholine binding32.32×10−520GO:0015271Outward rectifier potassium channel activity31.56×10−5(1) Bubble Plot. In the bubble plot, letters on the left are GO names, numbers on the bottom are the proportions of genes, sizes of circles indicate the numbers of enriched genes, and colors reflect p values. The redder the colors are, the more enriched the genes are, and the smaller the p values are (Figure 6).Figure 6
Bubble plot of GO enrichment analysis ofCalculus Bovis targets. Bubble plot: letters on the left are GO names, numbers on the bottom are the proportions of genes, sizes of the circles indicate the numbers of enriched genes, and colors reflect p values. The redder the colors are, the more enriched the genes are, and the smaller the p values are.(2) Column Chart. In the column chart, letters on the left are GO names, numbers on the bottom are the numbers of genes enriched on GO, and p reflects significance of enrichment. The redder the colors are, the more enriched the genes are, and the smaller the p values are (Figure 7).Figure 7
Column chart of GO functional enrichment analysis ofCalculus Bovis targets. Column chart: letters on the left are GO names, numbers on the bottom are the numbers of genes enriched on GO, and p reflects the significance of enrichment. The redder the colors are, the more enriched the genes are, and the smaller the p values are.


#### 3.5.2. GO Enrichment Analysis of Shared Targets ofCalculus Bovis and Ischemic Stroke

GO enrichment analysis was performed against the 92 shared targets. The top 20 were selected based on theirp values and numbers of enrichment, including endopeptidase activity (GO:0004175); protein serine/threonine kinase activity (GO:0004674); metallopeptidase activity (GO:0008237); phosphoric diester hydrolase activity (GO:0008081); heme binding (GO:0020037); tetrapyrrole binding (GO:0046906); oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen (GO:0016705); serine-type peptidase activity (GO:0008236); serine hydrolase activity (GO:0017171); carboxylic acid binding (GO:0031406); organic acid binding (GO:0043177); metalloendopeptidase activity (GO:0004222); RNA polymerase II transcription factor binding (GO:0001085); monocarboxylic acid binding (GO:0033293); and fatty acid binding (GO:0005504) (see Table 4). Results were visualized as a bubble plot (Figure 8) and a column chart (Figure 9) using the R package.Table 4
GO enrichment analysis results of shared targets byCalculus Bovis and ischemic stroke.
No.IDDescriptionCountp value1GO:0004175Endopeptidase activity155.14505×10−92GO:0004674Protein serine/threonine kinase activity122.81896×10−63GO:0008237Metallopeptidase activity94.20451×10−74GO:0008081Phosphoric diester hydrolase activity83.02315×10−85GO:0020037Heme binding85.05803×10−76GO:0046906Tetrapyrrole binding88.71647×10−77GO:0016705Oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen81.74788×10−68GO:0008236Serine-type peptidase activity84.78249×10−69GO:0017171Serine hydrolase activity85.61466×10−610GO:0031406Carboxylic acid binding87.36664×10−611GO:0043177Organic acid binding81.14448×10−512GO:0004222Metalloendopeptidase activity71.07557×10−613GO:0001085RNA polymerase II transcription factor binding71.62174×10−514GO:0033293Monocarboxylic acid binding69.84665×10−715GO:0005504Fatty acid binding58.40369×10−716GO:0043560Insulin receptor substrate binding42.19487×10−717GO:0036041Long-chain fatty acid binding46.57868×10−718GO:00041143′,5′-Cyclic-nucleotide phosphodiesterase activity41.27273×10−519GO:0004112Cyclic-nucleotide phosphodiesterase activity41.47051×10−520GO:00163031-Phosphatidylinositol-3-kinase activity31.58908×10−5Figure 8
Bubble plot of GO enrichment analysis of shared targets byCalculus Bovis and ischemic stroke.Figure 9
Column chart of GO enrichment analysis of shared targets byCalculus Bovis and ischemic stroke.


### 3.6. KEGG Pathway Enrichment Analysis

KEGG pathway enrichment analysis was performed for 344Calculus Bovis targets and 92 shared targets.


#### 3.6.1. KEGG Pathway Enrichment Analysis forCalculus Bovis Targets

Eighteen results were selected based on theirp values and their numbers of enrichment. They were mainly involved in neuroactive ligand-receptor interaction (hsa04080), Alzheimer disease (hsa05010), calcium signaling pathway (hsa04020), nicotine addiction (hsa05033), steroid hormone biosynthesis (hsa00140), GABAergic synapse (hsa04727), morphine addiction (hsa05032), cholinergic synapse (hsa04725), retrograde endocannabinoid signaling (hsa04723), and cellular senescence (hsa04218) (see Table 5).Table 5
KEGG pathway enrichment analysis results ofCalculus Bovis targets.
No.IDDescriptionCountp value1hsa04080Neuroactive ligand-receptor interaction141.03×10−82hsa05010Alzheimer disease94.14×10−63hsa04020Calcium signaling pathway63.4×10−54hsa05033Nicotine addiction52.08104×10−45hsa00140Steroid hormone biosynthesis52.31001×10−46hsa04727GABAergic synapse53.04528×10−47hsa05032Morphine addiction53.73161×10−48hsa04725Cholinergic synapse56.29559×10−49hsa04723Retrograde endocannabinoid signaling51.033601×10−310hsa04218Cellular senescence51.40748×10−311hsa04330Notch signaling pathway41.594499×10−312hsa04115p53 signaling pathway42.068787×10−313hsa05204Chemical carcinogenesis42.107734×10−314hsa04540Gap junction42.532292×10−315hsa04970Salivary secretion42.652055×10−316hsa04914Progesterone-mediated oocyte maturation43.294802×10−317hsa05030Cocaine addiction33.367718×10−318hsa04340Hedgehog signaling pathway33.567181×10−3(1) Bubble Plot. In the bubble plot, letters on the left are KEGG names, numbers on the bottom are the proportions of genes, sizes of circles indicate the numbers of enriched genes, and colors reflect p values. The redder the colors are, the more enriched the genes are, and the smaller the p values are (Figure 10).Figure 10
Bubble plot of KEGG pathway enrichment analysis ofCalculus Bovis targets. Bubble plot: letters on the left are KEGG names, numbers on the bottom are the proportions of genes, sizes of the circles indicate the numbers of enriched genes, and colors reflect p values. The redder the colors are, the more enriched the genes are, and the smaller the p values are.(2) Column Chart. In the column chart, letters on the left are KEGG names, numbers on the bottom are the numbers of genes enriched on KEGG, columns represent genes enriched on KEGG, and p reflects significance of enrichment. The redder the colors are, the more enriched the genes are, and the smaller the p values are (Figure 11).Figure 11
Column chart of KEGG pathway enrichment analysis ofCalculus Bovis targets. Column chart: letters on the left are KEGG names, numbers on the bottom are the numbers of genes enriched on KEGG, columns represent genes enriched on KEGG, and p reflects the significance of enrichment. The redder the colors are, the more enriched the genes are, and the smaller the p values are.


#### 3.6.2. KEGG Pathway Enrichment Analysis for Shared Targets

KEGG pathway enrichment analysis was performed on 92 shared targets, and 145 pathways were found. The top 20 candidates were selected based on theirp values and numbers of enrichment, and they were involved in the PI3K-AKT signaling pathway (hsa04151), human papillomavirus infection (hsa05165), Kaposi sarcoma-associated herpesvirus infection (hsa05167), human cytomegalovirus infection (hsa05163), microRNAs in cancer (hsa05206), fluid shear stress and atherosclerosis (hsa05418), proteoglycans in cancer (hsa05205), AGE-RAGE signaling pathway in diabetic complications (hsa04933), cAMP signaling pathway (hsa04024), prostate cancer (hsa05215), TNF signaling pathway (hsa04668), platelet activation (hsa04611), EGFR tyrosine kinase inhibitor resistance (hsa01521), small cell lung cancer (hsa05222), and endocrine resistance (hsa01522) (Table 6). Results were visualized as a bubble plot (Figure 12) and a column chart (Figure 13) using the R package.Table 6
KEGG pathway enrichment analysis for shared targets.
No.IDDescriptionCountp value1hsa04151PI3K-AKT signaling pathway242.12708×10−132hsa05165Human papillomavirus infection202.45479×10−103hsa05167Kaposi sarcoma-associated herpesvirus infection191.3448×10−134hsa05163Human cytomegalovirus infection182.23341×10−115hsa05206MicroRNAs in cancer184.23388×10−96hsa05418Fluid shear stress and atherosclerosis178.25476×10−147hsa05205Proteoglycans in cancer174.88669×10−118hsa04933AGE-RAGE signaling pathway in diabetic complications166.21323×10−159hsa04024cAMP signaling pathway161.04013×10−910hsa05215Prostate cancer141.54571×10−1211hsa04668TNF signaling pathway141.16611×10−1112hsa04611Platelet activation136.34665×10−1013hsa01521EGFR tyrosine kinase inhibitor resistance123.67514×10−1114hsa05222Small cell lung cancer122.32309×10−1015hsa01522Endocrine resistance124.93414×10−1016hsa04919Thyroid hormone signaling pathway125.80768×10−917hsa05212Pancreatic cancer114.36655×10−1018hsa05210Colorectal cancer111.70773×10−919hsa04657IL-17 signaling pathway114.49095×10−920hsa05220Chronic myeloid leukemia107.44115×10−9Figure 12
Bubble plot of KEGG pathway enrichment analysis of shared targets.Figure 13
Column chart of KEGG pathway enrichment analysis of shared targets.


## 3.1.Calculus Bovis Compound-Target Network

A total of 344 targets were found by 12 bioactive compounds ofCalculus Bovis. Details of these ingredients are listed in Table 1, and the map of the compound-target network is shown in Figure 2. Circles represent compounds of Calculus Bovis, and triangles represent the targets; their colors are darkened as their degrees of freedom increased. There were 362 nodes and 706 edges with a network density of 0.011 and a network diameter of 5. Details of these parameters are listed in Table 2.Table 1
Active ingredients ofCalculus Bovis.
Mol IDMol name2D structureOB (%)DLMOL000263Oleanolic acid29.020.76MOL000298Ergosterol14.290.72MOL000511Ursolic acid16.770.75MOL008834Bilicerdin23.790.75MOL0088353-[2-[[3-(2-Carboxyethyl)-4-methyl-5-[(E)-(4-methyl-5-oxo-3-vinyl-2-pyrrolylidene)methyl]-1H-pyrrol-2-yl]methyl]-4-methyl-5-[(Z)-(3-methyl-5-oxo-4-vinyl-2-pyrrolylidene)methyl]-1H-pyrrol-3-yl]propanoic acid16.530.75MOL008838Methyl (4R)-4-[(3R,5S,7S,8R,9S,10S,12S,13R,14S,17R)-3,7,12-trihydroxy-10,13-dimethyl-2,3,4,5,6,7,8,9,11,12,14,15,16,17-tetradecahydro-1H-cyclopenta[a]phenanthren-17-yl]pentanoate32.320.76MOL008839Methyl desoxycholate34.630.73MOL0088402-[(3alpha,12alpha-Dihydroxy-24-oxo-5beta-cholan-24-yl)amino]ethanesulfonic acid15.920.87MOL008843Cherianoine27.320.12MOL008846ZINC0128036546.380.49MOL000953CLR37.870.68MOL009807CHD22.170.72Abbreviations: OB = oral bioavailability; DL = drug likeness; Mol = molecular.Figure 2
Calculus Bovis compound-target network. Note: circles represent compounds, triangles represent targets, and their colors darken as their degrees of freedom increase.Table 2
Calculus Bovis compound-candidate target network parameters.
Network parametersValuesNumber of nodes362Network density0.011Network diameter5Network heterogeneity3.212Average number of neighbors3.901Characteristic path length3.236Shortest paths130682 (100%)Network centralization0.304


## 3.2. Disease Targets of Ischemic Stroke

Using key words listed in Section2.1.3, 74 disease targets were found in the Drugbank database, 313 were found in the DisGeNet database after screening for values above the average, 322 were found in the OMIM database, and 10 were found in the TTD database. A total of 590 targets were found after deleting gene duplicates.


## 3.3. Prediction ofCalculus Bovis Targets in Ischemic Stroke

After entering compound targets and disease targets to Venny 2.1.0 in the form of gene names, 92 overlapping genes were found (see Figure3). These genes were the shared targets of the bioactive compounds and disease targets in ischemic stroke.Figure 3
Venny diagram ofCalculus Bovis targets and ischemic stroke disease targets. CB: Calculus Bovis target; IS: ischemic stroke disease target.


## 3.4. Construction of PPI Networks

### 3.4.1. The PPI Network ofCalculus Bovis Targets

PPI networks have been widely used in studying protein-protein interactions in different diseases. To construct the PPI network ofCalculus Bovis targets, 12 compounds were connected with 344 targets in the TCMSP and the SwissTargetPrediction databases, and they were imported into the STRING database (species: Homo sapiens; minimum required interaction score: high confidence (0.700)), and the PPI network was visualized after reconstructing it with Cytoscape (version 3.7.0). As shown in Figure 4, colors darkened as the degree of freedom increased. This PPI network contained 322 nodes and 2195 edges with a diameter of 8 and an average number of neighbors of 13.634. It showed that TP53 (degree=82), AKT1 (degree=67), MAPK1 (degree = 67), PIK3CA (degree=63), SRC (degree=58), MAPK3 (degree=56), VEGFA (degree=52), HSP90AA1 (degree=48), JUN (degree=46), EGFR (degree=46), MAPK8 (degree=45), TNF (degree=40), and CASP3 (degree=40) were the key nodes of this PPI network.Figure 4
Protein-protein interaction (PPI) network ofCalculus Bovis compound targets. Darkness of colors and sizes of circles represent the degree of freedom as the scale indicates.


### 3.4.2. The PPI Network ofCalculus Bovis-Ischemic Stroke Targets and Hub Genes

To explore the potential therapeutic mechanisms ofCalculus Bovis in managing ischemic stroke, 92 shared targets by Calculus Bovis compounds and ischemic stroke were connected and imported into the STRING database as shown in Figure 5(a). The PPI network of Calculus Bovis-ischemic stroke targets were constructed by visualizing the results using Cytoscape (Figure 5(b)). This PPI network had 83 nodes and 403 edges with a network diameter of 6 and an average number of neighbors of 9.711. The top 10 targets with the highest degree of freedom were TP53 (degree=33), AKT1 (degree=30), MAPK1 (degree=29), VEGFA (degree=27), TNF (degree=25), PIK3CA (degree=25), MAPK3 (degree=24), MMP9 (degree=22), PTGS2 (degree=22), and IL-1B (degree=20). A total of 6 Hub genes were found using 3 algorithms: degree of freedom, Maximum Neighborhood Component (MNC), and Maximal Clique Centrality (MCC) (Figure 5(c)). These genes were TP53, AKT1, PIK3CA, MAPK3, MMP9, and MMP2. Their network of interactions is shown in Figure 5(d). In addition, TP53, AKT1, MAPK1, VEGFA, PIK3CA, and MAPK3 were among the top 10 candidates ranked by degree of freedom.Figure 5
The PPI network ofCalculus Bovis compound targets against ischemic stroke (IS): (a) the original PPI data generated from the STRING database showing detailed interactions of the targets; (b) the PPI network constructed using Cytoscape (version 3.7.0); (c) the top 15 genes were calculated from the PPI network by the degree of freedom, MNC, and MCC, and the overlapping genes were then screened by Venn diagrams; (d) the PPI network of hub genes.
(a)(b)(c)(d)


## 3.4.1. The PPI Network ofCalculus Bovis Targets

PPI networks have been widely used in studying protein-protein interactions in different diseases. To construct the PPI network ofCalculus Bovis targets, 12 compounds were connected with 344 targets in the TCMSP and the SwissTargetPrediction databases, and they were imported into the STRING database (species: Homo sapiens; minimum required interaction score: high confidence (0.700)), and the PPI network was visualized after reconstructing it with Cytoscape (version 3.7.0). As shown in Figure 4, colors darkened as the degree of freedom increased. This PPI network contained 322 nodes and 2195 edges with a diameter of 8 and an average number of neighbors of 13.634. It showed that TP53 (degree=82), AKT1 (degree=67), MAPK1 (degree = 67), PIK3CA (degree=63), SRC (degree=58), MAPK3 (degree=56), VEGFA (degree=52), HSP90AA1 (degree=48), JUN (degree=46), EGFR (degree=46), MAPK8 (degree=45), TNF (degree=40), and CASP3 (degree=40) were the key nodes of this PPI network.Figure 4
Protein-protein interaction (PPI) network ofCalculus Bovis compound targets. Darkness of colors and sizes of circles represent the degree of freedom as the scale indicates.


## 3.4.2. The PPI Network ofCalculus Bovis-Ischemic Stroke Targets and Hub Genes

To explore the potential therapeutic mechanisms ofCalculus Bovis in managing ischemic stroke, 92 shared targets by Calculus Bovis compounds and ischemic stroke were connected and imported into the STRING database as shown in Figure 5(a). The PPI network of Calculus Bovis-ischemic stroke targets were constructed by visualizing the results using Cytoscape (Figure 5(b)). This PPI network had 83 nodes and 403 edges with a network diameter of 6 and an average number of neighbors of 9.711. The top 10 targets with the highest degree of freedom were TP53 (degree=33), AKT1 (degree=30), MAPK1 (degree=29), VEGFA (degree=27), TNF (degree=25), PIK3CA (degree=25), MAPK3 (degree=24), MMP9 (degree=22), PTGS2 (degree=22), and IL-1B (degree=20). A total of 6 Hub genes were found using 3 algorithms: degree of freedom, Maximum Neighborhood Component (MNC), and Maximal Clique Centrality (MCC) (Figure 5(c)). These genes were TP53, AKT1, PIK3CA, MAPK3, MMP9, and MMP2. Their network of interactions is shown in Figure 5(d). In addition, TP53, AKT1, MAPK1, VEGFA, PIK3CA, and MAPK3 were among the top 10 candidates ranked by degree of freedom.Figure 5
The PPI network ofCalculus Bovis compound targets against ischemic stroke (IS): (a) the original PPI data generated from the STRING database showing detailed interactions of the targets; (b) the PPI network constructed using Cytoscape (version 3.7.0); (c) the top 15 genes were calculated from the PPI network by the degree of freedom, MNC, and MCC, and the overlapping genes were then screened by Venn diagrams; (d) the PPI network of hub genes.
(a)(b)(c)(d)


## 3.5. GO Enrichment Analysis

To further investigate the underlying mechanism ofCalculus Bovis in managing ischemic stroke, 344 targets and 92 shared targets were collected for GO enrichment analysis.


### 3.5.1. GO Enrichment Analysis forCalculus Bovis Targets

Twenty results were selected based on theirp values and their numbers of enrichment. They were primarily involved in ion channel activity (GO:0005216), channel activity (GO:0015267), passive transmembrane transporter activity (GO:0022803), gated channel activity (GO:0022836), cation channel activity (GO:0005261), metal ion transmembrane transporter activity (GO:0046873), neurotransmitter receptor activity (GO:0030594), ligand-gated ion channel activity (GO:0015276), ligand-gated channel activity (GO:0022834), postsynaptic neurotransmitter receptor activity (GO:0098960), extracellular ligand-gated ion channel activity (GO:0005230), neurotransmitter receptor activity involved in regulation of postsynaptic membrane potential (GO:0099529), transmitter-gated ion channel activity (GO:0022824), transmitter-gated channel activity (GO:0022835), and steroid dehydrogenase activity (GO:0016229). Details are listed in Table 3. Results are presented using a bubble plot and a column chart using the R package.Table 3
GO enrichment analysis results ofCalculus Bovis targets.
No.IDDescriptionCountp value1GO:0005216Ion channel activity141.15918×10−102GO:0015267Channel activity143.57104×10−103GO:0022803Passive transmembrane transporter activity143.67058×10−104GO:0022836Gated channel activity137.88226×10−115GO:0005261Cation channel activity111.25701×10−86GO:0046873Metal ion transmembrane transporter activity111.85161×10−77GO:0030594Neurotransmitter receptor activity91.68357×10−108GO:0015276Ligand-gated ion channel activity81.76915×10−89GO:0022834Ligand-gated channel activity81.76915×10−810GO:0098960Postsynaptic neurotransmitter receptor activity73.22861×10−911GO:0005230Extracellular ligand-gated ion channel activity61.78299×10−712GO:0099529Neurotransmitter receptor activity involved in regulation of postsynaptic membrane potential56.54501×10−713GO:0022824Transmitter-gated ion channel activity51.78703×10−614GO:0022835Transmitter-gated channel activity51.78703×10−615GO:0016229Steroid dehydrogenase activity45.40×10−616GO:0016628Oxidoreductase activity, acting on the CH-CH group of donors, NAD or NADP as acceptor41.58×10−617GO:1904315Transmitter-gated ion channel activity involved in regulation of postsynaptic membrane potential41.78589×10−518GO:0022851GABA-gated chloride ion channel activity39.83×10−619GO:0042166Acetylcholine binding32.32×10−520GO:0015271Outward rectifier potassium channel activity31.56×10−5(1) Bubble Plot. In the bubble plot, letters on the left are GO names, numbers on the bottom are the proportions of genes, sizes of circles indicate the numbers of enriched genes, and colors reflect p values. The redder the colors are, the more enriched the genes are, and the smaller the p values are (Figure 6).Figure 6
Bubble plot of GO enrichment analysis ofCalculus Bovis targets. Bubble plot: letters on the left are GO names, numbers on the bottom are the proportions of genes, sizes of the circles indicate the numbers of enriched genes, and colors reflect p values. The redder the colors are, the more enriched the genes are, and the smaller the p values are.(2) Column Chart. In the column chart, letters on the left are GO names, numbers on the bottom are the numbers of genes enriched on GO, and p reflects significance of enrichment. The redder the colors are, the more enriched the genes are, and the smaller the p values are (Figure 7).Figure 7
Column chart of GO functional enrichment analysis ofCalculus Bovis targets. Column chart: letters on the left are GO names, numbers on the bottom are the numbers of genes enriched on GO, and p reflects the significance of enrichment. The redder the colors are, the more enriched the genes are, and the smaller the p values are.


### 3.5.2. GO Enrichment Analysis of Shared Targets ofCalculus Bovis and Ischemic Stroke

GO enrichment analysis was performed against the 92 shared targets. The top 20 were selected based on theirp values and numbers of enrichment, including endopeptidase activity (GO:0004175); protein serine/threonine kinase activity (GO:0004674); metallopeptidase activity (GO:0008237); phosphoric diester hydrolase activity (GO:0008081); heme binding (GO:0020037); tetrapyrrole binding (GO:0046906); oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen (GO:0016705); serine-type peptidase activity (GO:0008236); serine hydrolase activity (GO:0017171); carboxylic acid binding (GO:0031406); organic acid binding (GO:0043177); metalloendopeptidase activity (GO:0004222); RNA polymerase II transcription factor binding (GO:0001085); monocarboxylic acid binding (GO:0033293); and fatty acid binding (GO:0005504) (see Table 4). Results were visualized as a bubble plot (Figure 8) and a column chart (Figure 9) using the R package.Table 4
GO enrichment analysis results of shared targets byCalculus Bovis and ischemic stroke.
No.IDDescriptionCountp value1GO:0004175Endopeptidase activity155.14505×10−92GO:0004674Protein serine/threonine kinase activity122.81896×10−63GO:0008237Metallopeptidase activity94.20451×10−74GO:0008081Phosphoric diester hydrolase activity83.02315×10−85GO:0020037Heme binding85.05803×10−76GO:0046906Tetrapyrrole binding88.71647×10−77GO:0016705Oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen81.74788×10−68GO:0008236Serine-type peptidase activity84.78249×10−69GO:0017171Serine hydrolase activity85.61466×10−610GO:0031406Carboxylic acid binding87.36664×10−611GO:0043177Organic acid binding81.14448×10−512GO:0004222Metalloendopeptidase activity71.07557×10−613GO:0001085RNA polymerase II transcription factor binding71.62174×10−514GO:0033293Monocarboxylic acid binding69.84665×10−715GO:0005504Fatty acid binding58.40369×10−716GO:0043560Insulin receptor substrate binding42.19487×10−717GO:0036041Long-chain fatty acid binding46.57868×10−718GO:00041143′,5′-Cyclic-nucleotide phosphodiesterase activity41.27273×10−519GO:0004112Cyclic-nucleotide phosphodiesterase activity41.47051×10−520GO:00163031-Phosphatidylinositol-3-kinase activity31.58908×10−5Figure 8
Bubble plot of GO enrichment analysis of shared targets byCalculus Bovis and ischemic stroke.Figure 9
Column chart of GO enrichment analysis of shared targets byCalculus Bovis and ischemic stroke.


## 3.5.1. GO Enrichment Analysis forCalculus Bovis Targets

Twenty results were selected based on theirp values and their numbers of enrichment. They were primarily involved in ion channel activity (GO:0005216), channel activity (GO:0015267), passive transmembrane transporter activity (GO:0022803), gated channel activity (GO:0022836), cation channel activity (GO:0005261), metal ion transmembrane transporter activity (GO:0046873), neurotransmitter receptor activity (GO:0030594), ligand-gated ion channel activity (GO:0015276), ligand-gated channel activity (GO:0022834), postsynaptic neurotransmitter receptor activity (GO:0098960), extracellular ligand-gated ion channel activity (GO:0005230), neurotransmitter receptor activity involved in regulation of postsynaptic membrane potential (GO:0099529), transmitter-gated ion channel activity (GO:0022824), transmitter-gated channel activity (GO:0022835), and steroid dehydrogenase activity (GO:0016229). Details are listed in Table 3. Results are presented using a bubble plot and a column chart using the R package.Table 3
GO enrichment analysis results ofCalculus Bovis targets.
No.IDDescriptionCountp value1GO:0005216Ion channel activity141.15918×10−102GO:0015267Channel activity143.57104×10−103GO:0022803Passive transmembrane transporter activity143.67058×10−104GO:0022836Gated channel activity137.88226×10−115GO:0005261Cation channel activity111.25701×10−86GO:0046873Metal ion transmembrane transporter activity111.85161×10−77GO:0030594Neurotransmitter receptor activity91.68357×10−108GO:0015276Ligand-gated ion channel activity81.76915×10−89GO:0022834Ligand-gated channel activity81.76915×10−810GO:0098960Postsynaptic neurotransmitter receptor activity73.22861×10−911GO:0005230Extracellular ligand-gated ion channel activity61.78299×10−712GO:0099529Neurotransmitter receptor activity involved in regulation of postsynaptic membrane potential56.54501×10−713GO:0022824Transmitter-gated ion channel activity51.78703×10−614GO:0022835Transmitter-gated channel activity51.78703×10−615GO:0016229Steroid dehydrogenase activity45.40×10−616GO:0016628Oxidoreductase activity, acting on the CH-CH group of donors, NAD or NADP as acceptor41.58×10−617GO:1904315Transmitter-gated ion channel activity involved in regulation of postsynaptic membrane potential41.78589×10−518GO:0022851GABA-gated chloride ion channel activity39.83×10−619GO:0042166Acetylcholine binding32.32×10−520GO:0015271Outward rectifier potassium channel activity31.56×10−5(1) Bubble Plot. In the bubble plot, letters on the left are GO names, numbers on the bottom are the proportions of genes, sizes of circles indicate the numbers of enriched genes, and colors reflect p values. The redder the colors are, the more enriched the genes are, and the smaller the p values are (Figure 6).Figure 6
Bubble plot of GO enrichment analysis ofCalculus Bovis targets. Bubble plot: letters on the left are GO names, numbers on the bottom are the proportions of genes, sizes of the circles indicate the numbers of enriched genes, and colors reflect p values. The redder the colors are, the more enriched the genes are, and the smaller the p values are.(2) Column Chart. In the column chart, letters on the left are GO names, numbers on the bottom are the numbers of genes enriched on GO, and p reflects significance of enrichment. The redder the colors are, the more enriched the genes are, and the smaller the p values are (Figure 7).Figure 7
Column chart of GO functional enrichment analysis ofCalculus Bovis targets. Column chart: letters on the left are GO names, numbers on the bottom are the numbers of genes enriched on GO, and p reflects the significance of enrichment. The redder the colors are, the more enriched the genes are, and the smaller the p values are.


## 3.5.2. GO Enrichment Analysis of Shared Targets ofCalculus Bovis and Ischemic Stroke

GO enrichment analysis was performed against the 92 shared targets. The top 20 were selected based on theirp values and numbers of enrichment, including endopeptidase activity (GO:0004175); protein serine/threonine kinase activity (GO:0004674); metallopeptidase activity (GO:0008237); phosphoric diester hydrolase activity (GO:0008081); heme binding (GO:0020037); tetrapyrrole binding (GO:0046906); oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen (GO:0016705); serine-type peptidase activity (GO:0008236); serine hydrolase activity (GO:0017171); carboxylic acid binding (GO:0031406); organic acid binding (GO:0043177); metalloendopeptidase activity (GO:0004222); RNA polymerase II transcription factor binding (GO:0001085); monocarboxylic acid binding (GO:0033293); and fatty acid binding (GO:0005504) (see Table 4). Results were visualized as a bubble plot (Figure 8) and a column chart (Figure 9) using the R package.Table 4
GO enrichment analysis results of shared targets byCalculus Bovis and ischemic stroke.
No.IDDescriptionCountp value1GO:0004175Endopeptidase activity155.14505×10−92GO:0004674Protein serine/threonine kinase activity122.81896×10−63GO:0008237Metallopeptidase activity94.20451×10−74GO:0008081Phosphoric diester hydrolase activity83.02315×10−85GO:0020037Heme binding85.05803×10−76GO:0046906Tetrapyrrole binding88.71647×10−77GO:0016705Oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen81.74788×10−68GO:0008236Serine-type peptidase activity84.78249×10−69GO:0017171Serine hydrolase activity85.61466×10−610GO:0031406Carboxylic acid binding87.36664×10−611GO:0043177Organic acid binding81.14448×10−512GO:0004222Metalloendopeptidase activity71.07557×10−613GO:0001085RNA polymerase II transcription factor binding71.62174×10−514GO:0033293Monocarboxylic acid binding69.84665×10−715GO:0005504Fatty acid binding58.40369×10−716GO:0043560Insulin receptor substrate binding42.19487×10−717GO:0036041Long-chain fatty acid binding46.57868×10−718GO:00041143′,5′-Cyclic-nucleotide phosphodiesterase activity41.27273×10−519GO:0004112Cyclic-nucleotide phosphodiesterase activity41.47051×10−520GO:00163031-Phosphatidylinositol-3-kinase activity31.58908×10−5Figure 8
Bubble plot of GO enrichment analysis of shared targets byCalculus Bovis and ischemic stroke.Figure 9
Column chart of GO enrichment analysis of shared targets byCalculus Bovis and ischemic stroke.


## 3.6. KEGG Pathway Enrichment Analysis

KEGG pathway enrichment analysis was performed for 344Calculus Bovis targets and 92 shared targets.


### 3.6.1. KEGG Pathway Enrichment Analysis forCalculus Bovis Targets

Eighteen results were selected based on theirp values and their numbers of enrichment. They were mainly involved in neuroactive ligand-receptor interaction (hsa04080), Alzheimer disease (hsa05010), calcium signaling pathway (hsa04020), nicotine addiction (hsa05033), steroid hormone biosynthesis (hsa00140), GABAergic synapse (hsa04727), morphine addiction (hsa05032), cholinergic synapse (hsa04725), retrograde endocannabinoid signaling (hsa04723), and cellular senescence (hsa04218) (see Table 5).Table 5
KEGG pathway enrichment analysis results ofCalculus Bovis targets.
No.IDDescriptionCountp value1hsa04080Neuroactive ligand-receptor interaction141.03×10−82hsa05010Alzheimer disease94.14×10−63hsa04020Calcium signaling pathway63.4×10−54hsa05033Nicotine addiction52.08104×10−45hsa00140Steroid hormone biosynthesis52.31001×10−46hsa04727GABAergic synapse53.04528×10−47hsa05032Morphine addiction53.73161×10−48hsa04725Cholinergic synapse56.29559×10−49hsa04723Retrograde endocannabinoid signaling51.033601×10−310hsa04218Cellular senescence51.40748×10−311hsa04330Notch signaling pathway41.594499×10−312hsa04115p53 signaling pathway42.068787×10−313hsa05204Chemical carcinogenesis42.107734×10−314hsa04540Gap junction42.532292×10−315hsa04970Salivary secretion42.652055×10−316hsa04914Progesterone-mediated oocyte maturation43.294802×10−317hsa05030Cocaine addiction33.367718×10−318hsa04340Hedgehog signaling pathway33.567181×10−3(1) Bubble Plot. In the bubble plot, letters on the left are KEGG names, numbers on the bottom are the proportions of genes, sizes of circles indicate the numbers of enriched genes, and colors reflect p values. The redder the colors are, the more enriched the genes are, and the smaller the p values are (Figure 10).Figure 10
Bubble plot of KEGG pathway enrichment analysis ofCalculus Bovis targets. Bubble plot: letters on the left are KEGG names, numbers on the bottom are the proportions of genes, sizes of the circles indicate the numbers of enriched genes, and colors reflect p values. The redder the colors are, the more enriched the genes are, and the smaller the p values are.(2) Column Chart. In the column chart, letters on the left are KEGG names, numbers on the bottom are the numbers of genes enriched on KEGG, columns represent genes enriched on KEGG, and p reflects significance of enrichment. The redder the colors are, the more enriched the genes are, and the smaller the p values are (Figure 11).Figure 11
Column chart of KEGG pathway enrichment analysis ofCalculus Bovis targets. Column chart: letters on the left are KEGG names, numbers on the bottom are the numbers of genes enriched on KEGG, columns represent genes enriched on KEGG, and p reflects the significance of enrichment. The redder the colors are, the more enriched the genes are, and the smaller the p values are.


### 3.6.2. KEGG Pathway Enrichment Analysis for Shared Targets

KEGG pathway enrichment analysis was performed on 92 shared targets, and 145 pathways were found. The top 20 candidates were selected based on theirp values and numbers of enrichment, and they were involved in the PI3K-AKT signaling pathway (hsa04151), human papillomavirus infection (hsa05165), Kaposi sarcoma-associated herpesvirus infection (hsa05167), human cytomegalovirus infection (hsa05163), microRNAs in cancer (hsa05206), fluid shear stress and atherosclerosis (hsa05418), proteoglycans in cancer (hsa05205), AGE-RAGE signaling pathway in diabetic complications (hsa04933), cAMP signaling pathway (hsa04024), prostate cancer (hsa05215), TNF signaling pathway (hsa04668), platelet activation (hsa04611), EGFR tyrosine kinase inhibitor resistance (hsa01521), small cell lung cancer (hsa05222), and endocrine resistance (hsa01522) (Table 6). Results were visualized as a bubble plot (Figure 12) and a column chart (Figure 13) using the R package.Table 6
KEGG pathway enrichment analysis for shared targets.
No.IDDescriptionCountp value1hsa04151PI3K-AKT signaling pathway242.12708×10−132hsa05165Human papillomavirus infection202.45479×10−103hsa05167Kaposi sarcoma-associated herpesvirus infection191.3448×10−134hsa05163Human cytomegalovirus infection182.23341×10−115hsa05206MicroRNAs in cancer184.23388×10−96hsa05418Fluid shear stress and atherosclerosis178.25476×10−147hsa05205Proteoglycans in cancer174.88669×10−118hsa04933AGE-RAGE signaling pathway in diabetic complications166.21323×10−159hsa04024cAMP signaling pathway161.04013×10−910hsa05215Prostate cancer141.54571×10−1211hsa04668TNF signaling pathway141.16611×10−1112hsa04611Platelet activation136.34665×10−1013hsa01521EGFR tyrosine kinase inhibitor resistance123.67514×10−1114hsa05222Small cell lung cancer122.32309×10−1015hsa01522Endocrine resistance124.93414×10−1016hsa04919Thyroid hormone signaling pathway125.80768×10−917hsa05212Pancreatic cancer114.36655×10−1018hsa05210Colorectal cancer111.70773×10−919hsa04657IL-17 signaling pathway114.49095×10−920hsa05220Chronic myeloid leukemia107.44115×10−9Figure 12
Bubble plot of KEGG pathway enrichment analysis of shared targets.Figure 13
Column chart of KEGG pathway enrichment analysis of shared targets.


## 3.6.1. KEGG Pathway Enrichment Analysis forCalculus Bovis Targets

Eighteen results were selected based on theirp values and their numbers of enrichment. They were mainly involved in neuroactive ligand-receptor interaction (hsa04080), Alzheimer disease (hsa05010), calcium signaling pathway (hsa04020), nicotine addiction (hsa05033), steroid hormone biosynthesis (hsa00140), GABAergic synapse (hsa04727), morphine addiction (hsa05032), cholinergic synapse (hsa04725), retrograde endocannabinoid signaling (hsa04723), and cellular senescence (hsa04218) (see Table 5).Table 5
KEGG pathway enrichment analysis results ofCalculus Bovis targets.
No.IDDescriptionCountp value1hsa04080Neuroactive ligand-receptor interaction141.03×10−82hsa05010Alzheimer disease94.14×10−63hsa04020Calcium signaling pathway63.4×10−54hsa05033Nicotine addiction52.08104×10−45hsa00140Steroid hormone biosynthesis52.31001×10−46hsa04727GABAergic synapse53.04528×10−47hsa05032Morphine addiction53.73161×10−48hsa04725Cholinergic synapse56.29559×10−49hsa04723Retrograde endocannabinoid signaling51.033601×10−310hsa04218Cellular senescence51.40748×10−311hsa04330Notch signaling pathway41.594499×10−312hsa04115p53 signaling pathway42.068787×10−313hsa05204Chemical carcinogenesis42.107734×10−314hsa04540Gap junction42.532292×10−315hsa04970Salivary secretion42.652055×10−316hsa04914Progesterone-mediated oocyte maturation43.294802×10−317hsa05030Cocaine addiction33.367718×10−318hsa04340Hedgehog signaling pathway33.567181×10−3(1) Bubble Plot. In the bubble plot, letters on the left are KEGG names, numbers on the bottom are the proportions of genes, sizes of circles indicate the numbers of enriched genes, and colors reflect p values. The redder the colors are, the more enriched the genes are, and the smaller the p values are (Figure 10).Figure 10
Bubble plot of KEGG pathway enrichment analysis ofCalculus Bovis targets. Bubble plot: letters on the left are KEGG names, numbers on the bottom are the proportions of genes, sizes of the circles indicate the numbers of enriched genes, and colors reflect p values. The redder the colors are, the more enriched the genes are, and the smaller the p values are.(2) Column Chart. In the column chart, letters on the left are KEGG names, numbers on the bottom are the numbers of genes enriched on KEGG, columns represent genes enriched on KEGG, and p reflects significance of enrichment. The redder the colors are, the more enriched the genes are, and the smaller the p values are (Figure 11).Figure 11
Column chart of KEGG pathway enrichment analysis ofCalculus Bovis targets. Column chart: letters on the left are KEGG names, numbers on the bottom are the numbers of genes enriched on KEGG, columns represent genes enriched on KEGG, and p reflects the significance of enrichment. The redder the colors are, the more enriched the genes are, and the smaller the p values are.


## 3.6.2. KEGG Pathway Enrichment Analysis for Shared Targets

KEGG pathway enrichment analysis was performed on 92 shared targets, and 145 pathways were found. The top 20 candidates were selected based on theirp values and numbers of enrichment, and they were involved in the PI3K-AKT signaling pathway (hsa04151), human papillomavirus infection (hsa05165), Kaposi sarcoma-associated herpesvirus infection (hsa05167), human cytomegalovirus infection (hsa05163), microRNAs in cancer (hsa05206), fluid shear stress and atherosclerosis (hsa05418), proteoglycans in cancer (hsa05205), AGE-RAGE signaling pathway in diabetic complications (hsa04933), cAMP signaling pathway (hsa04024), prostate cancer (hsa05215), TNF signaling pathway (hsa04668), platelet activation (hsa04611), EGFR tyrosine kinase inhibitor resistance (hsa01521), small cell lung cancer (hsa05222), and endocrine resistance (hsa01522) (Table 6). Results were visualized as a bubble plot (Figure 12) and a column chart (Figure 13) using the R package.Table 6
KEGG pathway enrichment analysis for shared targets.
No.IDDescriptionCountp value1hsa04151PI3K-AKT signaling pathway242.12708×10−132hsa05165Human papillomavirus infection202.45479×10−103hsa05167Kaposi sarcoma-associated herpesvirus infection191.3448×10−134hsa05163Human cytomegalovirus infection182.23341×10−115hsa05206MicroRNAs in cancer184.23388×10−96hsa05418Fluid shear stress and atherosclerosis178.25476×10−147hsa05205Proteoglycans in cancer174.88669×10−118hsa04933AGE-RAGE signaling pathway in diabetic complications166.21323×10−159hsa04024cAMP signaling pathway161.04013×10−910hsa05215Prostate cancer141.54571×10−1211hsa04668TNF signaling pathway141.16611×10−1112hsa04611Platelet activation136.34665×10−1013hsa01521EGFR tyrosine kinase inhibitor resistance123.67514×10−1114hsa05222Small cell lung cancer122.32309×10−1015hsa01522Endocrine resistance124.93414×10−1016hsa04919Thyroid hormone signaling pathway125.80768×10−917hsa05212Pancreatic cancer114.36655×10−1018hsa05210Colorectal cancer111.70773×10−919hsa04657IL-17 signaling pathway114.49095×10−920hsa05220Chronic myeloid leukemia107.44115×10−9Figure 12
Bubble plot of KEGG pathway enrichment analysis of shared targets.Figure 13
Column chart of KEGG pathway enrichment analysis of shared targets.


## 4. Discussion

Network pharmacology was used in this study to investigateCalculus Bovis itself and its potential mechanism for the treatment of ischemic stroke through compound target network construction, PPI network analysis, GO enrichment analysis, and KEGG pathway analysis.Network analysis of compound targets showed that deoxycorticosterone (MOL008846), methyl cholate (MOL008838), and biliverdin (MOL008834) had the most connections with these targets, suggesting that these 3 compounds might be the key compounds ofCalculus Bovis.Deoxycorticosterone is a type of steroid hormone possessing activities of the mineralocorticoid and serves as the precursor of aldosterone. It is involved in water and salt metabolism, playing an important role in electrolyte balance and in the volume of body fluid [32]. It has been reported that deoxycorticosterone and its derivatives—neurosteroids transformed in the fetal brain—protect the central nervous system. Inhibiting the production of neurosteroids increases basal cell death [33]. Neuroactive steroid hormones are involved in the regulation of diverse physiological functions, such as cell differentiation, neuroprotection, memory reinforcement, and amelioration of anxiety and pressure [34]. Methyl cholate is the methyl ester form of cholic acid, inhibiting the synthesis of cholesterol [35]. It has been reported that methyl cholate suppresses the growth of certain Gram-positive and Gram-negative bacteria [36] and has a good anti-inflammatory effect [37]. Biliverdin is a type of bile pigment, an oxidized product of heme. Emerging evidence has shown that biliverdin is an endogenous antioxidant facilitating the restoration of the tissue oxidation-reduction environment [38]. In the middle cerebral artery occlusion (MCAO) model, it significantly decreased the infarct area and the production of peroxides in the cortex [39]. These indicate that biliverdin plays a pivotal role in mitigating ischemic brain injury through its antioxidative stress effect. In addition, a single target was regulated by multiple compounds as shown in our network. Protein-tyrosine phosphatase 1B (PTPN1) was subject to the regulation of deoxycorticosterone, oleanolic acid, ergosterol, ursolic acid, biliverdin, bilirubin, methyl cholate, methyl desoxycholate, taurodeoxychloic acid, and others. PTPs are involved in regulating differentiation and survival of neurons and have also been reported to be a new target of antiplatelet therapy [40]. PTPN1 is a negative regulator of the leptin and insulin signaling pathways, and PTP1B knockout mice are exempt from obesity and diabetes, both of which are risk factors of ischemic stroke [41].Similarly, other targets like 11-beta-hydroxysteroid dehydrogenase 1 (HSD11B1), dual-specificity phosphatase Cdc25A (CDC25A), cytochrome P450 19A1 (CYP19A1), progesterone receptor (PGR), and androgen receptor (AR) were also regulated by two or more compounds. The present study revealed not only relationships betweenCalculus Bovis compounds and their targets but also the potential pharmacological effects of Calculus Bovis, which reflects the multicompound and multitarget theory of modern drugs.Furthermore, the PPI network demonstrated information not only about protein homology and coexpression but also about protein-protein interactions. Our PPI analysis showed thatCalculus Bovis influences ischemic stroke through its impact on a complex biological network, including TP53, AKT1, MAPK1, VEGFA, TNF, PIK3CA, MAPK3, MMP9, PTGS2, and IL1B. Hub gene screening revealed that TP53, AKT1, PIK3CA, MAPK3, MMP9, and MMP2 were essential in this process. The above potential targets for the action of Calculus Bovis in the treatment of ischemic stroke are our first discoveries.TP53, cellular tumor antigen p53, is a tumor suppressing gene. It promotes cell apoptosis, increases gene stability, and suppresses tumorigenesis [42]. It has been reported that methylation of the TP53 promoter was increased in ischemic stroke, and this increase was associated with the thickness of the intima of the carotid artery, degree of atherosclerosis of the carotid artery, and levels of homocysteine in the peripheral blood [43]. More evidence showed that TP53 induced glycolysis, and apoptosis regulator (TIGAR) suppressed glycolysis, increased pentose-phosphate pathway flux, and maintained the function of mitochondria. As a result, it protected the brain from ischemic injury [44]. The Tp53 Arg/Arg genotype has been considered a genetic marker for predicting poor prognosis after ischemic stroke [45]. AKT1 (serine/threonine-protein kinase AKT) codes for the serine/threonine-protein kinase which regulates apoptosis proteins and transcription factors. It is the key player in regulating cell survival, growth, apoptosis, and proliferation in the presence of growth factors and external stimuli, especially brain ischemia and reperfusion injury. AKT1 gene variance is closely related to metabolic syndrome, a risk factor of stroke [46].AKT/PKB was involved in brain ischemia, and its activity was related to the extent of ischemic injury. Activation of AKT was the key factor in determining survival of neurons after ischemic injury [47]. PIK3CA codes for the p110 subunit of the phosphatidylinositol 3-kinases (PI3Ks), and its mutations decrease cell apoptosis and increase the activity of downstream PI3Ks. It is known that activation of PI3K/AKT promotes repair of neural injury due to ischemic stroke [48, 49] and angiogenesis in the hypoxic environment in vitro [50]. Consequently, it protects the rat brain from inflammation resulting from ischemia-reperfusion injury [51]. Activation of the TrkB/PI3K/AKT pathway also increases activation of Nrf2 and its translocation to the nucleus, which plays a pivotal role in protecting the central nervous system from oxidative stress [52].MAPK (mitogen-activated protein kinase) is involved in reaction to physiological and pathological stimuli, such as cytokines, neurotransmitters, hypoxia, and hypoglycemia [53]. MAPK3 (MAP kinase ERK1) and MAPK1 (MAP kinase ERK2) form ERK1/2, a subfamily of MAPK. They play important roles in cell proliferation, differentiation, migration, invasion, apoptosis, and other biological processes. In the H2O2 induced PC12 cell injury model, activation of the AKT and ERK1/2 pathways leads to an antioxidation effect [54]. Levels of ERK1 and ERK2 were increased after ischemic stroke onset. In vitro research further revealed that activation of ERK1/2 increased neuronal apoptosis, indicating that they are important targets for ischemic stroke treatment [55].GO and KEGG pathway enrichment analyses forCalculus Bovis targets showed that Calculus Bovis was closely related to ion channel activity, neurotransmitter receptor activity, and other physiological functions. Pathway enrichment analysis demonstrated its involvement in a number of pathways in the central nervous system, calcium-related signaling, and so on. These were consistent with the analgesic and antiepilepsy effects of Calculus Bovis and its clinical application in these fields, which lays the theoretical foundation for managing ischemic stroke using Calculus Bovis.GO enrichment and KEGG pathway analyses against the shared 92 targets showed thatCalculus Bovis was closely related to hydrolyzation of proteins, phosphorylation of serine/threonine residues of protein substrates, peptide bond hydrolyzation of peptides and proteins, hydrolyzation of second intracellular messengers, antioxidation and reduction, RNA transcription, and other biological processes. Among them, 15 (16.3%) were related to endopeptidase activities of the matrix metalloproteinase (MMPs) family (MMP1, MMP2, MMP3, MMP8, MMP9, MMP10, and MMP12) and the cysteine-containing aspartate-specific peptidase family (CASP1, CASP2, and CASP3). MMPs are a type of highly conserved proteinase in nature, belonging to the family of zinc-dependent endopeptidases. They are capable of degrading extracellular matrix, growth factors, cytokines, and cell adhesion molecules. They are indispensable in ECM and tissue remodeling, angiogenesis, immune reactions, inflammation, cell migration, proliferation, cell apoptosis, and other physiological and pathological processes [56]. Their mRNAs were dramatically increased in the cortex of the mouse ischemia model in which thrombus was induced [57]. Decreasing levels of MMPs significantly ameliorated transgression of neutrophils, which resulted in neuroprotection against ischemic stroke [58]. Targeted inhibition of the PI3K/AKT/MMP-9 signaling pathway suppressed tumor invasion and metastasis [59], but this has not been tested in neuronal cells. Caspases are closely related to apoptosis of eukaryocytes and regulation of cell proliferation as well as differentiation. Caspase 3 is the most important and indispensable one in the cascade of cell apoptosis. In brain ischemia, levels of CASP3 mRNA and protein were increased, and its activity was significantly reversed by a CASP3 inhibitor, preventing the hydrolyzation of poly(ADP-ribose) polymerase. Consequently, apoptosis was suppressed and neurological functions improved [60]. Neuroinflammation is a key pathological process, in which CASP1-activated inflammasomes play an essential role. It has been reported that CASP1 was increased in the mouse brain after cerebral ischemia, which was suppressed by a CASP1 inhibitor through decreasing the activation of microglial cells, protecting the brain from ischemic injury. These indicate that CASP1 is a potential drug target for ischemic stroke management [61]. Both in vitro and in vivo experiments have shown that activation of the PI3K/AKT pathway increases the phosphorylation of Bad and decreases the level of caspase-3, through which apoptosis is suppressed [62].Another 12 genes (13.0%) are related to activities of protein serine/threonine kinases, consistent with the findings of the PI3K-AKT signaling pathway (24 enriched genes, 27.2%) in the KEGG pathway analysis. AKT is also named serine/threonine kinase, whose activation is key for neuronal survival [63]. It is known that phosphorylated AKT was decreased in the infarct area after ischemia-reperfusion injury and increased in the penumbra after reperfusion. An inhibitor of PI3K decreased the level of phosphorylated AKT and increased the infarct area, indicating that PI3K/AKT is involved in the pathogenesis of brain ischemia and activation of AKT increases neuronal survival [64]. Among them, MAPK3 (ERK1) is not only a Hub gene involved in regulating the activity of protein serine/threonine kinases as shown by GO analysis but is also enriched in the PI3K-AKT signaling pathway. These suggest that MAPK3 is an important component of this network. Therefore, Calculus Bovis might protect the brain from ischemic stroke through its anti-inflammatory and antiapoptosis effects which are accomplished by interacting with the abovementioned key genes and pathways.However, there are some limitations in this study. Access to databases of Chinese herbal medicines and disease targets is relatively limited, and there may be some data selection bias. Meanwhile, due to the fact that this study was based on data analysis, the conclusion needs to be confirmed by further research.


## 5. Conclusion

Bioactive compounds, potential targets, and underlying mechanisms ofCalculus Bovis were examined using network pharmacology methods. KEGG pathway analysis showed that the PI3K/AKT and the MAPK signaling pathways were the key targets for ischemic stroke treatment. The effect of Calculus Bovis was achieved through directly or indirectly regulating the above targets and pathways. Our results confirmed that Calculus Bovis was a multicompound and multitarget drug with a multisystem character in the treatment of ischemic stroke, which laid the theoretical foundation for the development of drugs and therapeutic methods based on the results of Calculus Bovis in the future.



---

*Source: 6611018-2021-03-11.xml*