Abstract

In traditional Chinese medicine (TCM), Acori Tatarinowii Rhizoma (ATR) is widely used to treat memory and cognition dysfunction. This study aimed to confirm evidence regarding the potential therapeutic effect of ATR on Alzheimer’s disease (AD) using a system network level based in silico approach. Study results showed that the compounds in ATR are highly connected to AD-related signaling pathways, biological processes, and organs. These findings were confirmed by compound-target network, target-organ location network, gene ontology analysis, and KEGG pathway enrichment analysis. Most compounds in ATR have been reported to have antifibrillar amyloid plaques, anti-tau phosphorylation, and anti-inflammatory effects. Our results indicated that compounds in ATR interact with multiple targets in a synergetic way. Furthermore, the mRNA expressions of genes targeted by ATR are elevated significantly in heart, brain, and liver. Our results suggest that the anti-inflammatory and immune system enhancing effects of ATR might contribute to its major therapeutic effects on Alzheimer’s disease.

1. Introduction

Alzheimer's disease (AD) is one of the most common age-related serious fatal neurodegenerative diseases, and its prominent feature is progressive cognitive, learning, and memory dysfunction. Gradually, bodily functions are lost, eventually leading to death [1]. According to a statistical report issued by Alzheimer's disease association in 2016, one new case of AD is expected to develop every 33 seconds, resulting in nearly 1 million new cases per year, and the estimated prevalence is expected to range from 11 million to 16 million by 2050 [2]. Therefore, prevention and treatment of AD are particularly important. Although AD is a global health problem, its pathology remains poorly understood [3]. About 70% of the risk is thought to be genetic, and many genes are usually involved. The disease process is associated with high density of senile plaques and neurofibrillary tangle in brain [4]. Currently, no drugs or supplements can reduce the processes [2].

Recent study shows that AD might result from a series of complex reasons, such as aging, immunity dysfunction, and metabolic disturbances [5, 6]. In order to understand the pathogenesis of complex diseases a new “multitarget, multidrug” model strategy of drug discovery was presented [7]. Interestingly, as a multidisciplinary system of multicomponent therapeutics, Traditional Chinese Medicine (TCM) Systems revealed a probability to explain the relationship between multicomponents and drug synergistic effects [8]. Acori Tatarinowii Rhizoma (ShiChangPu, ATR) is a well-known herb from Chinese traditional medicine (TCM) being used for neurodegenerative diseases in China at least 2,000 years. ATR was recorded in the Chinese pharmacopoeia (2015) as the official botanical source. It contains several phytochemicals, and the main ingredients include β-asarone, α-asarone, estragole, cis-methylisoeugenol, isoshyobunone, δ-cadinene, and methyleugenol. Some experimental researches by cellular and animal models support that ATR had anti-AD effects. Volatile oils of ATR were being used in AD for promoting neural progenitor proliferation, protecting β-amyloid induced neurotoxicity, and improving memory and cognitive function [911]. β-Asarone, a major essential oil component from ATR, has anti-AD effects by nerve growth factor signaling pathways, ameliorating oxidative stress, autophagy, and neuronal apoptosis in vivo [1214]. In addition, many studies had confirmed that ATR has therapeutic effects on nervous disorders, depression, epilepsy, and dementia [1519].

In recent years, systems pharmacology has emerged as a new field that integrates pharmacology, biochemistry, genomics, and bioinformatics techniques, to establishe a model based on computer network analysis and therapeutic targets prediction, and it highlights the holistic thinking of TCM [20]. It provides a new powerful tool to analyze and visualize the complex interaction data in herb, compound, target, and disease. To introduce the TCM application on AD by using systems pharmacology approach, Cai’s group used a systems pharmacology approach to reveal the underlying molecular mechanisms of BSYZ in treating AD [21]. Luo’s group established multiple herbs-compounds-targets-pathway-cooperation networks model to analyze the effects of Danggui-Shaoyao-San in AD [22]. Fang et al. performed literature mining of PubMed for top 10 anti-AD herbs (include ATR), based on systematic pharmacological analysis, and they interpreted the multiscale mechanisms of action of herbs in AD management [23]. Although these researches provided valuable support for the mechanisms of TCM in the treatment of AD, there is no analysis on the mechanism of the single herb, such as ATR, because they were based on the combined action of multiple herbs.

In this study, we discussed how system pharmacology demonstrated the relationship of ATR to treat AD from the system level. First, ingredients of ATR were collected from traditional Chinese medicine systems pharmacology (TCMSP, http://lsp.nwu.edu.cn/) database [24]. Oral bioavailability (OB), drug-likeness (DL), and blood-brain barrier (BBB) analysis were carried out to filter the compounds for subsequent analysis. Next, system pharmacology approach was used to explore the compounds-targets interactions and to establish the compounds-targets network and the compounds-targets-AD network. Finally, based on bioinformatics analysis, we illuminated the multiscale mechanisms of action of ATR on AD. A workflow of the systems pharmacological study was summarized as shown in Figure 1.

2. Material and Methods

2.1. Establishment of Database

All the compounds of ATR are collected by TCMSP database (a free phytochemical database of herbal medicine). Evaluating the absorption, distribution, metabolism, and excretion (ADME) parameters: oral bioavailability (OB), drug-likeness (DL), and blood-brain barrier (BBB) for all the natural compounds, the active compounds were collected to illustrate the proposed model.

2.2. OB Evaluation

OB shows the percentage of an active compound that reaches the systemic circulation, which defines one of the most important pharmacokinetic parameters to the convergence of the ADME process [25]. OB is essential to determine whether the chemical components of TCM have pharmacodynamic activity. The OB values are predicted by using the OBioavail 1.1 [26], which is a tool that contains 805 different structures of drugs and drug-like molecules and metabolism information. In this study, the OB threshold was set as ≥15%.

2.3. DL Prediction

DL helps filter out “drug-like" compounds in the traditional Chinese herbs, as DL is a qualitative concept used in drug design for an estimate on how “drug-like” a prospective compound is [27]. Compounds with high DL values are more likely to have certain biological properties that can increase the possibility for drug treatment. The calculation of DL value is based on Tanimoto coefficient formula in TCMSP database [28].

2.4. BBB Screening

The blood-brain barrier (BBB) is anatomically characterized by the presence of intercellular tight junctions between continuous nonfenestrated endothelial cells, whose normal function is to limit the passage of protein, potentially diagnostic and therapeutic agents into the brain parenchyma [29]. Understanding and evaluating the capacity of compounds of entering into the central nervous system, the compounds with BBB <-0.3 were considered as nonpenetrating (BBB-), those from -0.3 to +0.3 moderate penetrating (BBB±), and those >0.3 strong penetrating (BBB+) [8].

2.5. Target Fishing

To obtain active compounds of ATR target information, a comprehensive drug targeting protocol includes text mining database search and chemometric analysis was applied. First, the information on molecular targets was found from the TCMSP. Second, based on the set-wise chemical similarity among their ligands of active compounds, the targets were adopted by Similarity Ensemble Approach (SEA, http://sea.bkslab.org) [30] and the Binding Database (http://www.bindingdb.org) [31]. Third, the targets obtained from different databases were input to UniProt (http://www.uniprot.org/) [32] to make the targets names uniformly standardized. Finally, noise and error information was eliminated by PharGKB (https://www.pharmgkb.org/) [33], Therapeutic Target Database (DTT, http://bidd.nus.edu.sg/group/cjttd/) [34], and the Comparative Toxicogenomics Database (CTD, http://ctdbase.org/) [35].

2.6. Drug-Target Network

In order to understand the complex relationship between ATR active compounds and potential targets, a visual network is established by CytoScape_v3.4.0 [36]. This network is composed of node and edge. Nodes represent to molecules (compounds and targets), and edges indicate intermolecular interactions (compounds and targets interactions), namely, the connections between nodes.

2.7. Drug-Target-AD Network

DisGeNET database (http://www.disgenet.org/) [37] is a discovery platform that serves as one of the largest publicly available databases, containing a mass of collections of genes and variants associated with human diseases. It contains 561,119 gene-disease associations and 135,588 variant-disease associations. Based on expert curated repositories, GWAS catalogues, animal models, and the scientific literature, it can be used for different research purposes, including the generation of hypothesis on drug therapeutic action and drug adverse effects. In order to elucidate the action relations of ATR in Alzheimer's disease, the integration is performed by means of gene and AD vocabulary mapping and using the DisGeNET association type ontology to select the targets related to AD. The active compounds-targets-AD network is established by CytoScape_v3.4.0.

2.8. Target-Organ Location Network

BioGPS is a database for querying and organizing genetic annotation resources [38]. It provides gene expression abundance data in cells or tissues by microarrays analysis [39]. The target-organ location network is constructed with the use of the dataset: GeneAtlas U133A, gcrma (http://biogps.org). First, the mRNA expression patterns of each target gene are obtained in 84 parts of organ tissues. Second, average values for each gene are applied to computer. Third, genes will locate in the relevant organ tissues where the mRNA expression level is higher than the mean. Finally, a target-organ location network is constructed using Cytoscape 3.4.0 and organ-specific, Alzheimer's disease related gene overexpression data.

2.9. Gene Ontology (GO) Analysis

Gene Ontology (GO) is a framework for the model of biology. The GO defines gene classes that describe gene functions and relationships between these concepts. It covers three domains: Cellular Component (CC) explains where gene products are active; Molecular Function (MF) means molecular activities of gene products; Biological Process (BP) interprets the pathways and larger processes made up of the activities of multiple gene products [40]. In this study, GO terms with P values < 0.01 and Benjamini < 0.05 were employed and the data were collected by the DAVID 6.8 Gene Functional Classification Tool [41] (http://david.abcc.ncifcrf.gov/).

2.10. Network Pathway Analysis

KEGG pathway enrichment analysis provides not only pathway functional annotations of given gene set but also pathway enrichment analysis [42]. First, gene numbers of each pathway in the given gene set are calculated. Then, significantly enriched pathways in the given gene set compared to the genome background are defined by hypergeometric test. The calculating formula of P value isHere N is the number of all genes with KEGG pathway annotation (http://www.genome.jp/kegg/) [43] which represents the gene number of given genes set in N; M is the number of all genes that are annotated to the certain pathway; m is the gene number in the given gene set that is annotated to the certain pathways. Taking P values < 0.01 and FDR ≤ 0.05 as threshold, pathways that meet this condition were defined as significantly enriched pathways in the given genes set. In this research, pathway enrichment analysis was performed by using the OmicShare tools, a free online platform for data analysis (http://www.omicshare.com/tools/). In order to illuminate the action mechanisms of ATR in Alzheimer's disease, filtered target proteins were entered into the pathway map of Alzheimer's disease acquired from the KEGG.

3. Results

3.1. Identification of Active Compounds

As 105 compounds were identified belonging to ATR in TCMSP (as shown in Table S1), 30 compounds were selected with the conditions of OB≥15%, BBB≥-0.3, and DL≥0.1. Removing 2 compounds (4-aromadendrene, isocembrol) with no target information, 28 compounds with known target information were chosen for the following analysis. Another 4 compounds (β-asarone, γ-asarone, cis-methylisoeugenol, and methyleugenol) reported to have anti-inflammation, antioxidation, and antimicrobial effects were added additionally [13, 4446], and finally 32 compounds were analyzed as shown in Table 1.

Most of these 32 compounds showed various biological activities experimentally. For example, eudesmin (M23; DL=0.62, BBB=0.05, and OB=52.35) was experimentally identified which has anti-inflammation [47] and anticonvulsion effects [48]. A angiogenesis inhibitor of marmesin (M28; DL=0.18, BBB=0.07, and OB=50.28) was reported to have regulatory effect on endothelial cell fate and angiogenesis [49]. Antioxidation and anti-inflammation effect of patchoulene (M31; DL=0.11, BBB=2.17, and OB=49.06) was proved through inhibition of nuclear factor κB (NF-κB) and downregulation of COX-2/iNOS [50]. Bergapten (M12; DL=0.13, BBB=0.69, and OB=42.21) was reported to have anti-inflammatory properties [51] and the ability of ameliorating the insulin resistance [52]. Aminacrine (M9; DL=0.12, BBB=0.33, and OB=35.00) and its derivatives were experimentally identified to be reversible inhibitors of acetyl cholinesterase (AChE) [53]. α-Cedrene (M2; DL=0.10, BBB=2.16, and OB=55.56) was reported to exhibit antimicrobial activity [54]. Veraguensin (M7; DL=0.39, BBB=0.72, and OB=25.53) was reported to have analgesic and anti-inflammatory activities [55] and neuroprotective effects[56]. β-Asarone (M16; DL=0.06, BBB=1.24, and OB=35.61) increased striatal level of dopamine by enhancing dopa decarboxylase activity which improved behavioral competence in Parkinson's rat model [57]. γ-Asarone (M24; DL=0.06, BBB=1.33, and OB=22.76) prevented oxidative stress-induced cell injury in cultured astrocytes through Akt signaling activation [46]. Aristolone (M11; DL=0.13, BBB=1.54, and OB=45.31) was reported to have regulatory effects on lipid peroxidation and proliferation in human cancer cells [58]. α-Cubebene (M10; DL=0.11, BBB=2.1, and OB=16.73) was experimentally identified with significant antifungal activity inhibition of A. flavus growth [59]. β-Gurjunene (M15; DL=0.10, BBB=2.07, and OB=51.36) was experimentally identified with antioxidation activity and antimicrobial activity [60]. Isopimpinellin (M26; DL=0.17, BBB=0.5, and OB=25.93) was reported to exhibit antiallergic inflammation [61]. α-Longipinene (M3; DL=0.12, BBB=2.05, and OB=57.47) was reported to inhibit biofilm formation and hyphal growth in Candida albicans [62]. Anti-inflammation effect of aristolene (M11; DL=0.11, BBB=2.08, and OB=52.2) was proved through measuring secretion of the proinflammatory cytokines IL-6 and TNF-α by using a human-derived macrophage cell line [63]. Calarene (M20; DL=0.11, BBB=2.04, and OB=52.16) and caryophyllene oxide (M4; DL=0.13, BBB=1.76, and OB=32.67) from the essential oil of Patrinia scabiosaefolia were both found to have antineuroinflammation, anticancer, and antioxidation effects[64]. As mentioned above, many compounds of ATR, which are universal in herbs, plants, and fruits, contribute to neuroprotective, anti-inflammation, antimicrobial, and antioxidation effects.

3.2. Target Fishing

These 32 identified active compounds interact with 181 target proteins (as shown in Table S2) based on a target fishing technique [65]; that is, on average, they interact with 5.7 target genes, which did fully explain the multiple-target effects of pharmacology by ATR [66]. Among the 181 obtained targets entered to DisGeNET database for screening AD-related targets, only 97 potential targets (as shown in Table S3) were associated with AD, and they were reserved for further analysis.

3.3. Network Construction

To facilitate the visualization of multiple-target effects of ATR and AD interrelation, network analysis was performed based on the context of the whole human genome [67]. Two networks (A) C-T and (B) C-T-AD were constructed that might show multicompound and multitarget effects and they might explain why ATR can be used for the treatment of AD.

3.3.1. Compound-Target Network (C-T Network)

C-T interactions data between 181 targets and 32 compounds of ATR were collected by the method in Section 2.5. The network comprised 213 nodes and 509 edges. Circular nodes represented targets and triangle nodes represented compounds. Node size was regulated by degree centrality, resulting in an average degree of 4.70 nodes per target and 15.91 edges per compound, respectively (as shown in Figure 2). Among the 32 active compounds, 10 compounds with high-degree distributions were identified, and each of these compounds hit more than 21 potential targets. Usually, higher degree meant that the compounds had more important pharmacological effects on the potential treatment of AD [68]. For example, veraguensin (M7, degree =82) had the highest number of interactions with targets, followed by methyleugenol (M29, degree =47) and cis-methylisoeugenol (M21, degree = 43), indicating that those compounds might have more pharmacological effects than other compounds in this C-T network. Meanwhile, muscarinic acetylcholine receptor M3 (CHRM3, degree =20) showed the most intimate connections with compounds, followed by muscarinic acetylcholine receptor M1 (CHRM1, degree = 19) and prostaglandin G/H synthase 2 (PTGS2, degree=18). The C-T network results demonstrated the multitarget effects of active compounds by ATR.

3.3.2. Compound-Target-AD Network (C-T-AD Network)

Network pharmacological methods provided a visual approach to understanding the complex relationship between disease and therapeutic spots [69, 70]. In the present study, 97 potential targets of AD were selected form DisGeNET database and 30 involved compounds (M27 and M30 have no target) were used to construct the C-T-AD network for further cluster analysis. The C-T-AD network was built in the same way as C-T network, as shown in Figure 3. Veraguensin (M7, degree=42), cis-methylisoeugenol (M21, degree=28), methyleugenol (M29, degree =24), and β-asarone (M13, degree=20) might play an important role in the treatment of AD. CHRM1 (degree = 20), PTGS2 (degree = 19), CHRM2 (degree = 19), and ADRB2 (degree = 9) might be the important drug-targets. 26 active compounds connected with more than two targets and all of the 97 targets interact with more than one compound, indicating that many proteins associated with AD might display similar binding patterns with ligands.

3.4. Multiple-Target/Organ Cooperation for AD Treatment

Traditional Chinese medical workers consider the human body as an organic whole, and none of the organs in this whole organization are independent. Theses organs are structurally inseparable, functionally coordinated, and complementary to each other. They interact with each other in pathology. Therefore, in this study, the complex process of exploring ATR for AD was used by multiorgan cooperation and multitarget action.

3.4.1. GO Analysis

GOBP describes a series of events accomplished by one or more organized assemblies of molecular functions [40]. GOBP showed that these targets were enriched to 20 biological process terms, and the 181 genes were highly related to response to drug, signal transduction, inflammatory response, response to hypoxia, aging, memory, and so on (as shown in Figure 4). In detail, phosphatidylinositol-4, 5-bisphosphate 3-kinase catalytic subunit, gamma isoform (PIK3CG), toll-like receptor 4 (TLR4), tumor necrosis factor receptor superfamily member 1A (TNFRSF1A), protooncogene c-fos (FOS), nuclear factor erythroid 2-related factor 2 (NFE2L2), macrophage migration inhibitory factor (MIF), and nuclear factor NF-kappa-B p105 subunit (NFKB1) were involved. 14 genes were related to “inflammatory response,” including sodium-dependent serotonin transporter (SLC6A4), matrix metalloproteinase-2 (MMP-2), hypoxia-inducible factor 1-alpha (HIF1A), dopamine D2 receptor (DRD2), inducible nitric oxide synthase (NOS2), and cyclic AMP-responsive element-binding protein 1 (CREB1). 12 genes were related to “response to hypoxia,” such as NFE2L2, Alpha-1A adrenergic receptor (ADRA1A), neuromedin-K receptor (TACR3), 5-hydroxytryptamine 2A receptor (HTR2A), transcription factor p65 (RELA), 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGCR), sodium-dependent dopamine transporter (SLC6A3), FOS, and CREB1. 11 genes were related to “aging” including SLC6A4, HTR2A, PTGS2, DRD1, CREB1, neuronal acetylcholine receptor subunit alpha-7 (CHRNA7), glutamate receptor ionotropic, and NMDA 2A (GRIN2A). 11 genes were related to “memory.”

In conclusion, ATR is entirely possible as a therapeutical herb on AD, which is attributed to affect certain biological processes, such as aging, memory, anti-inflammation, and response to hypoxia.

3.4.2. Target Tissue Location

The mRNA specific expression patterns in tissues are important clues to the gene function [71]. Thus, they are of great significance to determine the tissue mRNA expression profiles of various genes at the organ level which is to explore the beneficial effects of ATR on AD [72]. In this paper, the results of mRNA expression were demonstrated by using BioGPS microarray analysis. 175 of the 181 targets had their expression profiles in 84 normal tissues (as shown in Table S3).

The 175 targets were generally found in human tissues but showed different levels of mRNA expression. We focused on some of these organizations, including brain, heart, kidney, liver, lung, and whole blood. A target was positioned to the organization where it had the highest mRNA expression level. For example, we compared the expression patterns of different tissues and observed that the 21 targets contained had higher mRNA expression levels in the brain than the remaining 84 normal tissues, 27 targets were positioned in heart, and, meanwhile, 32 targets located in whole blood were linked with tissues in any forms. The organization distribution network of 175 targets (as shown in Figure 5) is based on their expression patterns.

Results from our research suggested that the 21 high-abundant targets in brain were regarded as therapeutic targets for AD, and the other targets were also thought to be associated with AD. Our results implied that drug targets worked on these tissues rather than focusing on a specific one, and the whole blood serves as bridge between these tissues, so as to promote the coordination of organs to achieve the positive effect of AD. The effects of AD are not just in the brain but also in other organs. For example, AD can cause a series of bodily functions damage incidents such as cardiac autonomic dysfunction, urinary incontinence, and metabolic disorders [2, 73, 74]. Meanwhile, other organ lesions may also affect the progression of AD. For example, insulin resistance, hyperglycemia, hyperinsulinemia, and type 2 diabetes mellitus could provoke AD [75]. Coronary heart disease and heart surgery correlate with the development of AD cognitive dysfunction [76, 77]. Therefore, the prevention of AD requires not only attention to the brain but also the coordination of other organs. An interpretation of these findings is that the holistic view of TCMs diagnosis and treatment theory is correct.

3.4.3. Alzheimer's Disease-Relating Pathway

Through the C-T-AD network, a macroscopic visual of the relationship between ATR, targets, and AD was obtained, but the detailed mechanism of these remains unclear. Therefore, for determining the interactions among ATR and AD multiple targets, we did the following things. First, 97 potential targets of AD underwent pathway enrichment analysis, and a holistic picture of KEGG pathway annotation constructed by the use of the OmicShare tools was provided (as shown in Figure 6). Second, an “AD-pathway” was established based on Alzheimer's disease pathway which was provided by the KEGG pathway database (as shown in Figure 9).

(1) KEGG Pathway Annotation. KEGG pathway annotation can be divided into several modules, such as “Environmental Information Processing”, “organismal systems”, and “metabolism.” Two representative modules, “environmental information processing” and “organismal systems,” were described in detail to clarify the underlying mechanism.

Signal transduction in “environmental information processing” module was studied, and a target-pathway network (T-P) was displayed (as shown in Figure 7) whose node size was regulated by degree centrality. 12 signal transduction pathways were concentrated.

There were 4 inflammation-related pathways shown in signal transduction module, including TNF signaling pathway (degree=9), NF-kappa B signaling pathway (degree=6), HIF-1 signaling pathway (degree=9), and VEGF signaling pathway (degree=4), and these pathways indicated that the anti-inflammatory action is important for the treatment of AD. The synthesis of inflammatory critical cytokines and receptors could induce a wide range of intracellular biological processes including apoptosis and cell survival as well as inflammation and immunity [78, 79]. We found that γ-asarone (M24) targets TNF receptor superfamily member 1A (TNFRSF1A), cis-methylisoeugenol (M21) and veraguensin (M7) target NF-kappa B p65 subunit (RELA), bergapten (M12) targets NF-kappa B p105 subunit (NFKB1), and eudesmin (M23) targets toll-like receptor 4 (TLR4), and so on.

The G-protein-coupled receptor-triggered signaling cascade related pathways about cAMP signaling pathway (degree=20) had the maximum interaction targets, suggesting that it was closely related to AD. Research shows that activation of cAMP-dependent pathway may lead to phosphorylate alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) [80], tau hyperphosphorylation, and memory impairment [81].

Calcium signaling pathway also plays an important role (degree=17). Calcium can act in signal transduction resulting from activation of ion channels or as a second messenger caused by indirect signal transduction pathways such as G-protein-coupled receptors [82]. In neurons, calcium is important for the synchronization of neuronal electrical activity with mitochondrial energy metabolism when it increases in cytosolic and mitochondrial system [83], and the increasing of Ca2+ signals due to Ca2+ release has been implicated to play roles in synaptic plasticity and memory, neurotransmitter release, and neuronal excitability[84].

The signaling pathway mediated by small GTPase affects the structure and function of synapses, such as Ras signaling pathway (degree=8) and Rap1 signaling pathway (degree=7). For example, Ras protein family (Ras, Rap1, and Rap2), which belong to small GTPase, regulate the secretion of presynaptic neurotransmitters and the transshipment of the postsynaptic glutamate receptors in the cell membrane [85].

PI3K-Akt signaling pathway (degree=12) is important to regulate the cell cycle and promote growth and proliferation neural stem cells (NSCs) differentiate into motor neurons in adult [86] and promote neuroprotection from NSCs injury [87]. 2′-O-Methylisoliquiritigenin (M8) targets glycogen synthase kinase-3 beta (GSK3β). GSK3β plays an important role in the ultraphosphorylation of tau, which is one of the pathological features of AD. Increased GSK3β activity causes PI3K/Akt dysfunction, which regulates glucose metabolism in the brain and leads to tau hyperphosphorylation in the brain of AD patients [88]. Moreover, in long-term potentiation (LTP), the PI3K binds to an AMPA receptor in a conservative region to determine the receptor in the membrane, especially in the GluR subunit, which causes AMPA receptors to be inserted postsynaptically [89].

In addition, cGMP-PKG signaling pathway (degree=12), sphingolipid signaling pathway (degree=6), and AMPK signaling pathway (degree=5) are also thought to be related to AD.

Nervous system in organismal systems module was also studied, and a target-pathway network (T-P) was established (as shown in Figure 8). Node size was regulated by degree centrality. They display extremely significantly close functional correlation to the AD, including dopaminergic synapse, cholinergic synapse, neurotrophin signaling pathway, serotonergic synapse, and long-term potentiation (LTP). Some researches showed that AD was also recognized as a disease of synaptic failure. Beta-secretase 1 (BACE1) cleaved amyloid precursor protein (APP) products impact learning and memory through proteins localized on glutamatergic, GABAergic, and dopaminergic synapses [90]. Acetylcholine is combined with both muscarinic acetylcholine receptor M1 (CHRM1) and muscarinic acetylcholine receptor M2 (CHRM2). To protect the memory and attention deficit caused by the loss of cholinergic neurons [91], acetylcholinesterase (ACHE) inhibitor works by partially blocking the degradation of acetylcholine in the synapse and enabling more of the neurotransmitter to reach and activate cholinergic receptors [92]. Soluble oligomers of the amyloid-β peptide (AβO) decreased brain serotonin (5-HT) levels in mice, whereas treatment with 5-HT prevented AβO-induced microglial activation and increased TNF-α levels [93]. We found that bisasarcin (M17) targets BACE1 and ACHE. (-)-Caryophyllene oxide (M4), veraguensin (M7), and aminacrine (M9) target ACHE. Among the 32 identified active compounds, 19 compounds target CHRM1 and 18 compounds target CHRM2, respectively. (-)-Alloaromadendrene (M1), β-asarone (M13), cis-methylisoeugenol (M21), and methyleugenol (M29) target 5-hydroxytryptamine 2A receptor (HTR2A).

(2) AD-Pathway Network. Nine protein targets were mapped to the AD-pathway network, indicating that these proteins play an important role in AD (as shown in Figure 9). AD is associated with senile plaques and neurofibrillary tangles (NFTs). Amyloid-beta (Aβ), a major component of senile plaques, has various pathological effects on cell and organelle function. Intracellular Aβ may contribute to pathology by facilitating tau hyperphosphorylation, disrupting mitochondria function, and triggering calcium dysfunction. As displayed in Figure 9. Some target receptors involved in Aβ aggregation and tau hyperphosphorylation were affected by ATR components to produce neuroprotective functions.

Recent studies have shown that the soluble oligomers of Aβ polypeptides may be a pathogenic factor for the development of Alzheimer's disease [94]. The amyloidogenic pathway triggered from two sequential cleavages of APP by BACE1 resulted in the formation of Aβ peptides [95]. Based on the above findings, blocking the enzyme (BACE inhibitor) in theory would prevent the accumulation of Aβ and (per the amyloid hypothesis) may help slow or stop Alzheimer's disease. For example, Albany Molecular Research Inc. reported on their development of an experimental drug for the treatment of Alzheimer's disease called verubecestat (MK-8931), which is an inhibitor of BACE1[96]. AstraZeneca and Eli Lilly and Company announced an agreement to codevelop lanabecestat (AZD3293), which is an oral beta-secretase 1 cleaving enzyme (BACE) inhibitor [97], and the drug has been developed with the FDA’s fast-track designation. Our research data indicated that veraguensin (M7) and cis-methylisoeugenol (M21) target APP and bisasarcin (M17) can target BACE1 being used to treat AD through decreasing production of amyloid.

Excessive or abnormal tau phosphorylation destroys the cytoplasm function and interferes with axonal transport, resulting in cell death [98]. GSK3β has a great influence on activated overphosphorylation factor. 2′-O-Methylisoliquiritigenin (M8) targets GSK3β with AD treatment by preventing tau overphosphorylation.

As far as we know, glutamate-mediated excitotoxicity is one of the major reasons for chronic neurodegeneration. Aminacrine (M9) targets GRIN2A (NMDAR), and calamendiol (M19) and γ-asarone (M24) have shown connections with GRIA2 (GLuR). GRIA2 is a mediated glutamate receptor, which can mediate excitatory synaptic transmission in synaptic hyperactivity [99]. NMDA subtypes of glutamate receptors are agents that rely on the death of cells dependent on Ca2+ [100], which can lead to degenerative disorders due to excessive or inappropriate activity [101]. Therefore, GRIN2A and GRIA2 are the intrinsic neural protection targets for regulating Ca2+ concentration. Meanwhile, aminacrine (M9) targets voltage-dependent calcium channel subunit alpha-2/delta-1 (CACNA2D1), a protein in a voltage-dependent calcium channel mediates calcium ions into the membrane polarization. GRIN2B and GRM5 are associated with calcified neuropathy. Therefore, neuronal Ca2+ signaling is partially targeted through ATR for the treatment of AD.

In addition, TNFRSF1A (TNFR) participated in the MAPK signaling pathway and has been associated with inflammation, which is a risk factor for AD.

4. Discussion and Conclusion

TCM is a complex system composed of multicomponent and multitarget, and the function of each component is synergistic. The composition of TCM is so complex that it is difficult to study its effects from the mixture and to elaborate the mechanism of action from the perspective of modern medicine. The systematic pharmacological method provides new ideas and perspectives for the research of complicated Chinese medicine system. The pharmacology of TCM system has established a model of interaction of elements, such as drug molecules, protein targets, tissues, and organs, to clarify and predict the efficacy and toxicity of TCM [102]. In this study, a systematic pharmacological method was used to study the active ingredients of ATR. The potential targets of pharmacodynamic components, gene ontology, network analysis, organ location network, and pathway enrichment analysis were employed to explore the relationship between active ingredients and targets of Alzheimer's disease. Our result showed that ATR and AD are highly correlated in biological processes, organs, and signaling pathways. Our main findings are as follows.

(1) Through defining components in this system, we have studied systematic pharmacological methods to integrate multiple technologies, including ADME-system assessment, drug targeting, and target tissue distribution. 32 active compounds were detected. These compounds interact with 181 different targets by drug targeting. The resulting C-T network indicated that compounds veraguensin (M7), methyleugenol (M29), and cis-methylisoeugenol (M21) and targets CHRM3, CHRM1, and PTGS2 are the key factors to play an important role in the drug-target interaction network.

(2) The identified 97 targets associated with Alzheimer's disease are essential to understanding the pharmacological mechanisms of ATR. The C-T-AD network has clearly elucidated the pharmacological action and mechanism of ATR on Alzheimer's disease. For example, aristolone (M11) mainly targets at acetylcholine receptors CHRM1 and CHRM2 to treat AD through the cholinergic synapse pathway. The resulting network builds a more comprehensive visualization of the C-T-AD interaction pattern.

(3) TCM is the embodiment of holistic medicine, which is characterized as multicomponent treatment, multiobjective/approach regulation, and multiorgan cooperation. The response to inflammatory, hypoxia, aging, and memory mechanism of ATR is also illustrated by GO analysis. For instance, the 14 potential targets from the 7 compounds were enriched to the biological process of inflammatory response, which indicated that anti-inflammation may be the key to ATR treatment of AD. The target tissue location results had shown that 27 and 48 of 185 targets of ATR were highly expressed in heart and liver, respectively. These organs were closely related to AD, and ATR was closely related to these organs in the meridian theory.

(4) KEGG pathway annotation analysis showed that the potential targets of the compounds were clustered in 12 signaling pathways, mainly related to inflammation, calcium signaling pathway, G-protein mediated signaling pathway, and energy metabolism, and effects on the nervous system are mainly manifested in the dopaminergic synapse, cholinergic synapse, neurotrophin signaling pathway, serotonergic synapse, and LTP. Finally, pathway-mapping result showed that 9 protein targets were closely related to the ATR treatment AD.

In summary, this study not only conducts a comprehensive analysis of relevant AD through systematic pharmacological methods and discovers potential active compounds but also explains the treating mechanism of ATR in Alzheimer’s disease, widely applied to clarify the effectiveness and mystery of the medicinal plant ATR; what is more, it provides an example for the treatment of complex diseases in the near future. Although the results are interesting, further research requires that the support of experimental data includes the confirmation of drug dose relationships. In addition, the present study needs to test the molecular mechanisms of the active compounds in vivo to support further assessment of potential clinical applications.

Data Availability

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

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The research is supported by the open fund project of Provincial Education Department of Hunan (15K088), the National Natural Science Foundation of China (Grant no. 81774129), and Hunan Administration of Traditional Chinese Medicine Science Foundation (Grant no. 201825).

Supplementary Materials

Table S1 shows the structure of 105 compounds in Acori Tatarinowii Rhizoma (ShiChangPu, ATR). Compounds information was found in the TCMSP (http://lsp.nwu.edu.cn/) and was input to Pubchem (https://www.ncbi.nlm.nih.gov/pccompound/) for compounds names to be uniformly standardized. Table S2 shows the targets prediction of 32 identified active compounds of ATR. Molecular targets information was found in the TCMSP (http://lsp.nwu.edu.cn/), SEA (http://sea.bkslab.org) and the Binding Database (http://www.bindingdb.org). The targets were obtained from different databases and input to UniProt (http://www.uniprot.org/) to make the targets names uniformly standardized. Noise and error information was eliminated by PharGKB (https://www.pharmgkb.org/), Therapeutic Target Database (DTT, http://bidd.nus.edu.sg/group/cjttd/), and the Comparative Toxicogenomics Database (CTD, http://ctdbase.org/). Table S3 shows the potential targets of 32 identified active compounds of ATR were associated with AD. In order to elucidate the action relations of ATR in Alzheimer's disease, the integration was performed by means of gene and AD vocabulary mapping by using the DisGeNET association type ontology to select the targets related to AD. Table S4 shows the 175 targets mRNA expression profiles in 84 normal tissues. The targets organ location network was used with the dataset: GeneAtlas U133A, gcrma (http://biogps.org). (Supplementary Materials)