- Research article
- Open Access
Systems pharmacology-based approach to investigate the mechanisms of Danggui-Shaoyao-san prescription for treatment of Alzheimer’s disease
BMC Complementary Medicine and Therapies volume 20, Article number: 282 (2020)
Alzheimer’s disease (AD) is the most common cause of dementia in the elderly, characterized by a progressive and irreversible loss of memory and cognitive abilities. Currently, the prevention and treatment of AD still remains a huge challenge. As a traditional Chinese medicine (TCM) prescription, Danggui-Shaoyao-san decoction (DSS) has been demonstrated to be effective for alleviating AD symptoms in animal experiments and clinical applications. However, due to the complex components and biological actions, its underlying molecular mechanism and effective substances are not yet fully elucidated.
In this study, we firstly systematically reviewed and summarized the molecular effects of DSS against AD based on current literatures of in vivo studies. Furthermore, an integrated systems pharmacology framework was proposed to explore the novel anti-AD mechanisms of DSS and identify the main active components. We further developed a network-based predictive model for identifying the active anti-AD components of DSS by mapping the high-quality AD disease genes into the global drug-target network.
We constructed a global drug-target network of DSS consisting 937 unique compounds and 490 targets by incorporating experimental and computationally predicted drug–target interactions (DTIs). Multi-level systems pharmacology analyses revealed that DSS may regulate multiple biological pathways related to AD pathogenesis, such as the oxidative stress and inflammatory reaction processes. We further conducted a network-based statistical model, drug-likeness analysis, human intestinal absorption (HIA) and blood-brain barrier (BBB) penetration prediction to uncover the key ani-AD ingredients in DSS. Finally, we highlighted 9 key ingredients and validated their synergistic role against AD through a subnetwork.
Overall, this study proposed an integrative systems pharmacology approach to disclose the therapeutic mechanisms of DSS against AD, which also provides novel in silico paradigm for investigating the effective substances of complex TCM prescription.
Alzheimer’s disease (AD) is a progressive neurodegenerative disease, which starts with mild memory loss and gradually results in severe impairment of broad executive and cognitive functions [1,2,3]. According to the data from World Health Organization in 2019, approximately 50 million people are suffering from dementia, while AD contributes to 60–70% of these cases . More alarming is a continuous rise in AD cases year by year, leading to increasing disability and high socioeconomic burden. Unfortunately, currently there are only four available anti-AD drugs (Donepezil, Rivastigmine, Galantamine, Memantine) introduced into the market, and no approved disease-modifying treatments (DMTs) exist for AD.
It is worth noting that more and more scientists have turned their attention to the design of multi-targeted drugs instead of the traditional “one target-one drug” perspective, which could intervene the complex AD pathogenesis via multiple molecular targets . As a complementary and alternative medicine, traditional Chinese medicine (TCM), especially the TCM prescription, has been widely applied in Asian countries for the prevention and treatment of complex and degenerative diseases, including AD . Generally, a herbal formula usually comprises multiple components which could act on multiple targets and exert its synergistic effects on disease . In spite of lacking effective options for AD, lots of classical TCM prescriptions (e.g. SuHeXiang Wan, TongLuoJiuNao) have potential as therapeutic drugs for alleviating AD symptoms [8,9,10]. Danggui-Shaoyao-san (DSS), also known as Toki-shakuyaku-san in Japan and Dangguijakyakan in Korea, is a herb combination employed to improve cognitive function, which consists of six Chinese herbs including Paeoniae Radix Alba (BaiShao, BS), Angelica Sinensis Radix (DangGui, DG), Atractylodis Macrocephalae Rhizoma (BaiZhu, BZ), Chuanxiong Rhizoma (ChuanXiong, CX), Alismatis Rhizoma (ZeXie, ZX), and Poria (FuLing, FL). According to a meta-analysis conducted on randomized controlled trials (RCTs), DSS has a positive effect on the scores of Mini-Mental State Examination (MMSE) and activities of daily living (ADL) of AD patients . Clinical study found that the regional cerebral blood flow (rCBF) of AD patients in the posterior cingulate were significantly higher and orientation to place tended to improve after treatment with DSS . Besides, previous pharmacological studies have deciphered that DSS could exert potentially therapeutic effects for AD via multiple mechanisms such as disrupting the aggregation of Aβ, attenuating inflammatory reaction, and adjusting the mitochondrial membrane permeability (Table 1). For instance, a study in 2014 has reported that DSS could improve learning and memory capacity in female SAMP8 mice through modulating estrogen, nitric oxide, and glycine in plasma or hippocampal tissue . Furthermore, Kou et al. demonstrated that DSS ameliorated memory dysfunction and protected the ultrastructure of cortex changed by aging, which may be beneficial for the treatment of AD . Despite current reported therapeutic mechanisms, the underlying mechanism of actions (MOAs) and active components of DSS against AD still remain indistinct. Thus, it is necessary to comprehensively investigate the molecular mechanism of DSS for treatment of AD and deeply understand their synergistic effects on multiple pathways.
Systems pharmacology is an emerging discipline which combines in silico network-based tools and experimental assays, aiming to elucidate the changes in the functions and reactions of human body induced by medicines . Based on the holistic principle, the herbal formulas in TCM show the characteristics of multi-component and multi-target in therapy in contrast to synthetic drugs. Due to the complexity of ingredients and targets, conventional experimental approaches are time-consuming and expensive for TCM research. Systems pharmacology offers effective strategy to explore the multi-component network target research pattern of Chinese herbal medicine and formula [16,17,18].
In this study, we attempted to systematically explore the MOAs of DSS for treatment of AD through an integrated systems pharmacology framework (Fig. 1). Firstly, we collected the chemical ingredients of DSS with known protein targets from our previous integrated natural products database . We further computationally predicted the putative targets of DSS via a network-based inference method. Subsequently, a global drug-target network of DSS was constructed by incorporating the known and predicted DTIs. Furthermore, multi-level systems pharmacology analysis methods, including molecular-function analysis, biological process analysis, target-function modules analysis, and KEGG pathway enrichment were performed to elucidate the MOAs of DSS against AD. Additionally, we developed a network-based model for identifying the active anti-AD components of DSS by mapping the high-quality AD disease genes into the global drug-target network. After integrating high-performance liquid chromatography (HPLC) analysis data, drug-likeness analysis, human intestinal absorption (HIA) and blood-brain barrier (BBB) penetration assessment, and network statistical model prediction results, we highlighted the key anti-AD ingredients in DSS and illustrated the specific synergistic mechanisms through subnetwork analysis.
Manual curation of genes associated with AD
The genes associated with AD were collected from six authoritative databases: 1) the Malacards database (https://www.malacards.org); 2) the Comparative Toxicogenomics Database (CTD) ; 3) DisGeNET ; 4) the GWAS Catalogue ; 5) the Human Gene Mutation Database (HGMD) database ; and 6) AlzBase database (http://alz.big.ac.cn/alzBase/summary/Gene). Only genes labeled with “Alzheimer’s disease” were extracted from the databases mentioned above. Additionally, for the AlzBase database, the top 100 genes were preserved for further investigation. Finally, 299 AD-related genes (Supplementary Table S1) were integrated after removing the duplicates.
Collection of herbal ingredients and known protein targets
The ingredients and their protein targets in each herb of DSS were collected from our previous integrated database , which includes experimental validated DTIs of natural products extracting from over 2000 literatures and five authoritative compound-protein interaction databases: ChEMBL (v21) , BindingDB  (accessed in Sep. 2017), STITCH , the Herbal Ingredients’ Target Database (HIT) , and the Traditional Chinese Medicine Integrated Database (TCMID) . Consequently, a total of 1042 herbal ingredients in DSS were obtained. The corresponding number of the collected compounds in DG, CX, FL, ZX, BS, BZ is 549, 351, 84, 42, 125, 170, respectively. The detailed structural information of the 1042 herbal ingredients is provided in Supplementary Table S2.
Network-based target prediction for DSS
In this study, target fishing is carried out to predict targets based on known DTIs via a balanced substructure-drug-target network-based inference (bSDTNBI) method . We previously developed predictive network models to identify new targets of natural products with bSDTNBI method, which prioritizes potential targets utilizing resource-diffusion processes for both known and new natural products . During the process, four parameters (α = β = 0.1, γ = − 0.5, and k = 2) of bSDTNBI were adopted. The first parameter α was introduced to balance the initial resource allocation of different node types, while β was utilized to adjust weighted values of different edge types. The third parameter γ was applied to balance the effect of hub nodes in resource diffusion processes, and the last parameter κ represented the number of resource-diffusion processes. We subsequently calculated four types of molecular fingerprints for each compound based on PaDEL-Descriptor (version 2.18) , including MACCS, PubChem, Substructure (FP4) and Klekota-Roth (KR). Compared with the other three generated predictive models, bSDTNBI_KR performed best given its highest values of the area under the receiver operating characteristic curve (AUC = 0.959). Eventually, the bSDTNBI_KR predictive model was chosen to identify new targets of each natural product. In this study, the top 20 putative targets for each compound with known targets were selected (Supplementary Table S3).
Identification of key anti-AD ingredients in DSS
The key anti-AD ingredients in DSS were determined by four steps. Firstly, we obtained the main components in DSS according to the HPLC analysis results reported in previous references [33, 34]. Subsequently, we applied the Lipinski’s Rule of Five (molecular weight < 500, AlogP < 5, H-bond donor < 5, H-bond acceptor < 10)  to assess the drug-likeness of these components. Thirdly, a machine learning-based ADMET Simulator tool  was further utilized to evaluate their BBB penetration and HIA properties. The 166-bit fingerprint descriptor MACCS was used for chemical structure representation. Finally, we built a network-based statistical model to prioritize anti-AD compounds in DSS based on drug–target network and AD-related genes. We hypothesize that a compound would exert high potential for the treatment of AD if its targets tend to be encoded by AD-related genes . Fisher’s exact test was applied to calculate the statistical significance of the enrichment of AD-related proteins in target profiles of each ingredient in DSS. The P-values were corrected by Benjamin–Hochberg method . Predicted compound–AD pairs with adjust-P (q) value lower than 0.05 were regarded as significant.
Network construction and enrichment analysis
The drug-target (D-T) network and target-function (T-F) network was constructed to explore the molecular mechanism of DSS for treatment of AD. Networks were generated and analyzed by Cytoscape (v3.2.1, http://www.cytoscape.org/) and Gephi (v0.9.2, https://gephi.org/). In each graphical network, compounds or genes or functional modules were presented by nodes, while interactions were encoded by edges. The degree of each node was calculated, which represents the number of edges linked to it, characterizing the most important nodes in a network. Gene ontology (GO) term and KEGG pathway enrichment analysis were performed by DAVID (The Database for Annotation, Visualization and Integrated Discovery, http://david.abcc.ncifcrf.gov) and Omicshare webservers. Significantly enriched pathways in drug-target set comparing to the genome background were defined by hypergeometric test. The calculated p-value (q) was gone through FDR Correction and q < 0.05 was considered as statistically significant. The statistical analysis was performed by Python platform (v3.2, http://www.python.org/).
Overlap analysis of herbal ingredients and targets in DSS
“Jun-Chen-Zuo-Shi”, also known as the “sovereign-minister-assistant-courier” principle , is a key theory of TCM which guides physicians to formulate herbal medicine . This principle represents the different roles played by each herb in a formula as well as its change rules on effect after compatibility of medicines. In the DSS prescription, BS, DG and CX belong to Jun-Chen class (primary role), while FL, ZX and BZ serve as the Zuo-Shi class (subsidiary role). To figure out the specific formulate principle of DSS for treatment of AD at a molecular level, overlap analyses of herbal ingredients and their targets were performed. Figure 2a indicates the relationship of the six herbs based on the number of overlapped ingredients between each two herbs. Obviously, multiple herbs share large common compounds. Among the six herbs in DSS, DG and CX own the highest number of duplicate compounds (n = 48), followed by DG and BZ (n = 27), CX and BZ (n = 19). Meanwhile, Fig. 2b shows these six herbs in DSS covered 491 targets in total. Among them, DG has the greatest number of targets. Interestingly, we found that 79 common targets existed in all of these six herbs simultaneously, which implied DSS might magnify its therapeutic effect via modulating these common targets. These targets include estrogen receptor, peroxisome proliferator-activated receptor (PPAR), inflammatory factor, transcription factor, cytochrome P450 family and lipoxygenases. Taking PPARA as an example, a recent study has demonstrated that PPARA is an important factor regulating autophagy in the clearance of Aβ. Activation of PPARA-mediated autophagy could reduce AD-like pathology and cognitive decline in a murine model . The detailed information of the common targets can be found in Supplementary Table S4.
Construction of drug-target network for DSS
To develop a global drug-target network for DSS, all the chemical ingredients in DSS were mapped into a known DTI network of natural products to acquire experimental validated DTIs from our previous integrated database . Moreover, the predicted DTIs were obtained with the bSDTNBI method. Finally, the network of DSS contains 19,293 DTIs connecting 937 unique compounds with 490 targets (37 AD targets and 453 non-AD targets) (Supplementary Table S5).
As illustrated in Fig. 3, most compounds are connected to multiple targets. The average target degree (D) for a compound is 20.6, while the average drug degree (K) of a target is 39.4. Among the 937 compounds, the top 10 with largest target degree (D) are quercetin (CID5280343, D = 100), apigenin (CID5280443, D = 70), luteolin (CID5280445, D = 59), kaempferol (CID5280863, D = 47), caffeic acid (CID689043, D = 43), gallic acid (CID370, D = 40), DL-glutamic acid (CID611, D = 37), Gamma-aminobutyric acid (CID119, D = 36), emodin (CID3220, D = 35) and scopoletin (CID5280460, D = 32). Within these compounds, 4 out of them have AD-labeled target degree greater than 5, including quercetin, apigenin, luteolin and Gamma-aminobutyric acid. Indeed, evidence is emerging that these compounds are highly correlated with AD treatment. For example, quercetin has been reported to reverse histological hallmarks of AD and exert a protective effect on cognitive and emotional function in aged 3xTg-AD mice . Moreover, apigenin owns potent neuroprotective properties to protect iPSC-derived AD neurons via reducing the frequency of spontaneous Ca2+ signals and significantly reducing caspase-3/7 mediated apoptosis . Meanwhile, 16 proteins are targeted by more than 300 compounds (K > 300). Among them, MAPT (K = 905), ACHE (K = 448), MAPK1 (K = 435) and PPARG (K = 320) are AD-related targets, which have been confirmed to play a vital role in AD pathobiology. For instance, a previous research has indicated that MAPK1 (ERK2) is positioned to phosphorylate normal tau which could be involved in the generation of paired helical filaments in AD . Besides, arachidonate 15-lipoxygenase (ALOX15), a protein targeted by 276 compounds (K = 276) across all the 6 herbs in DSS, is an important enzyme to regulate the homeostasis of DHA metabolism and esterification that are associated with cognitive function . A recent study found that the mean densities of ALOX15 were improved by DSS treatment in both hippocampus and cortex of APP/PS1 mice . These findings suggest the potential therapeutic targets of DSS for exerting anti-AD effects.
Target-function network (T-F network)
Different network regions may underlie different biological pathways, processes or cellular localizations. Here we extracted the targets of DSS with degree ≥40 to construct a target-function (T-F) network. The T-F network demonstrates AD-related biological processes and associated targets based on the functional annotation bioinformatics microarray (DAVID 6.8) analysis . As depicted in Fig. 4, this network consists of 209 target-function pairs linking 66 targets with 11 AD-related functional modules. The average number of functional modules for each target has reached up to 6, and 34 targets are associated with more than 4 functional modules. These modules include oxidative stress, metabolic process, inflammatory response, transport process, transcription process, hormonal regulation, biosynthetic process, catabolic process, cell proliferation and differentiation, and cell death, indicating the potential AD-related pathological process mediating by DSS. For instance, the T-F network suggests that DSS may regulate oxidative stress via reducing oxidation and regulation of nitric oxide. AD is highly correlated with oxidative stress and increased production of reactive oxygen species (ROS) has direct effect on synaptic activity and neurotransmission in neurons, hence leading to cognitive dysfunction . Previous in vivo study had confirmed that DSS could significantly relieve oxidative stress, such as reducing the Aβ25–35-induced lipid peroxidation and possessing protective effect against neuronal damage .
GO term and KEGG pathway enrichment analysis
In this study, GO term enrichment analysis was performed on the 73 targets with degree equal or larger than 40. Figure 5a shows the significantly enriched terms (q < 0.05) in biological process (BP) and molecular function (MF) categories. Obviously, most of the biological processes are related with AD, such as the negative regulation of inflammatory response (BP, GO:0050728)  and oxidation reduction process (BP, GO:0055114) . Previous literature has revealed that anti-inflammatory drugs could significantly reduce the incidence of AD . Moreover, accumulating evidence demonstrates that anti-oxidations are potentially promising strategies for safe and efficient treatment of AD . Similarly, it is worthy to note that a great number of targets are associated with a variety of MF terms, which plays a role in the pathogenesis of AD. For example, carbonate dehydratase (MF, GO:0004089), has been indicated to serve as a novel therapy for AD . In addition, transcription factor (MF, GO:0008134) such as XBP1s, was also reported to restore hippocampal synaptic plasticity and memory by control of the Kalirin-7 pathway in Alzheimer model which may be therapeutic implications in AD pathology . The detailed information of GO terms analysis data were provided in Supplementary Table S6.
Generally speaking, drug action does not only have regulatory effects on their related targets, but also play a role for metabolic enzymes and downstream proteins, as well as some pathways that are related with the specific disease . To identify the significant pathways involved in DSS for treatment of AD, targets (K ≥ 40) in the drug-target network were mapped onto their related pathways based on KEGG annotation (Supplementary Table S7). Figure 5b showcases the top 10 significant KEGG pathways (q < 0.05), and some of them are intensively correlated with AD treatment. For instance, PPAR signaling pathway (ko03320) has been demonstrated to increase phosphorylation of fibroblast growth factor 14 in the Tg2576 mouse model of AD . Furthermore, a recent study shows that bile secretion (ko04976) is associated with “A/T/N”(amyloid, tau, and neurodegeneration) AD biomarkers, providing further support for a role of bile acid pathway in AD pathophysiology . Taken together, these results indicate the novel pathways and biological processes that may be regulated by DSS for exerting anti-AD effects, which is worthy of further research by in vivo and in vitro experiments.
Integrated pathway analysis of DSS related to AD pathogenesis
In this analysis, protein targets of DSS were mapped into the pathways that directly related to the pathological process of AD to construct an AD-integrated pathway (Supplementary Table S8). Several functional modules are involved in the integrated pathway, such as inflammation, oxidative stress, apoptosis, and neurotrophin modulation. Here, four representative modules (Fig. 6) were chosen to illustrate the potential therapeutic mechanisms of DSS on the metabolic pathways towards AD.
Anti-oxidative stress module
Oxidative stress has been considered as a contributing factor in the pathogenesis and progression of AD. The increased ROS directly affects synaptic activity and neurotransmission in neurons leading to cognitive dysfunction . There is increasing evidence that hypoxia-inducible factor 1 (HIF-1) possesses neuroprotection in the setting of neuronal insults. HIF-1 promotes glycolysis and glucose metabolism, thus helping to produce reductive equivalents of NADH/NADPH which counters oxidative stress . Figure 6 shows that DSS prescription could act on multiple targets on HIF-1 signal pathway (e.g. HIF-1α, IGF1R, ERBB2), indicating its potential pathway of oxidative stress regulation for AD treatment and prevention.
Neuroinflammation regulation module
Neuroinflammation is an important process in neurodegeneration in AD, involving in a vicious cycle of amyloid deposition, neuronal damage, tangle formation, and death . Evidences indicate that Aβ-induced neuroinflammation is associated with high levels of pro-inflammatory cytokines, including TNF-α, interleukin-6 (IL-6), and interleukin-1β (IL-1β) . P38 mitogen-activated protein kinase (MAPK) is a crucial target for chronic inflammatory diseases . Previous studies imply that MAPK is significantly involved in glial activation and subsequent neuroinflammation, resulting in chronic neurotoxicity . As exhibited in Fig. 6, compounds in DSS target to several key proteins in the MAPK signaling pathway, suggesting the important pathway that DSS may participate for the regulation of neuroinflammation and amyloid genesis.
Modulation functions related to NF-κB signaling
Nuclear factor-kappa B (NF-κB) plays pivotal role in gene regulation and implicates in oxidative stress, apoptosis, and inflammation . NF-κB can be activated in most cell types, including neurons, astrocytes, microglia, oligodendrocytes and endothelial cells of neurovascular and cerebrovascular units. Recent researches demonstrated that NF-ĸB signaling pathway plays key physiological role in central nervous system, which serves important functions in cellular responses to neuronal injury and synaptic plasticity . Figure 6 illustrated that lots of compounds in DSS act on the targets of NF-ĸB signaling pathway, indicating the potential AD therapeutic mechanism mediated by DSS. Recent in vitro study showed that DSS has a protective effect on neuroinflammation in lipopolysaccharide (LPS)-stimulated BV-2 microglia cells through the TLRs/NF-κB signaling pathway .
Neurotrophin modulation module
Neurotrophic factors, such as brain-derived neurotrophic factor (BDNF) and neurotrophin 3 (NT3) are critical to the maintenance of the nervous system in the adult brain . Growing evidences suggested that up-regulation of BDNF can relieve cognitive impairments and learning deficits in AD , while NT3 can promote the proliferation and differentiation of bone marrow-derived NSCs into cholinergic neurons and elevate the levels of acetylcholine (ACh) . Figure 6 shows that protein targets of DSS are enriched in the neurotrophin signal pathway, which suggests the critical role that DSS may play on neurotrophin modulation.
Prescription simplification: uncovering the key anti-AD ingredients in DSS
The complexity of ingredients as well as their intricate target network in TCM prescription is always a huge challenge that hinders the translation from pre-clinical testing results to clinical outcome. Although previous HPLC analysis had preliminary estimated the main active components in DSS [33, 34], the specific ingredients that exert therapeutic effect against AD are difficult to determine. In this section, we further narrow up the study scope of compounds in DSS, for identifying the key anti-AD ingredients and optimum drug combinations. After the four-step screening (see Methods), nine key anti-AD ingredients in DSS were finally highlighted, including 3-butylphthalide, ligustilide, senkyunolide I, Z-butylidenephthalide, senkyunolide A, atractylenolide I, tetramethylpyrazine, ferulic acid, and gallic acid (Table 2). According to the previous reported HPLC analysis results, all of the nine ingredients are the major constituents in DSS extract (Table 2) [33, 34]. To investigate the basis of compositions against AD, we searched the previous literature for these 9 ingredients and surprisingly found that 6 of them were confirmed to possess anti-AD effects (Table 2). For instance, tetramethylpyrazine (TMP), one of the major bioactive compounds purified from CX, shows great potential against AD, which could restore the function of cholinergic neurons and protect against memory loss in AD rat model . Besides, ligustilide was reported to have neuroprotective effect via inducting α-Secretase processing of both APP and Klotho in AD mouse . The high hit rate of our prediction suggests the AD therapy potential of the rest three unreported ingredients (atractylenolide I, senkyunolide I, senkyunolide A), which deserve to be validated by experimental assays in the future. Interestingly, we also found that these 9 components are mainly derived from three herbs (DG, CX, BS) in DSS, which serve as “JUN” (BS) and “CHEN” (DG, CX) role, indicating the rationality of DSS in formulation.
To explore the MOAs for the nine key ingredients toward AD, a subnetwork was constructed by extracting the DTIs of these compounds from the global drug-target network of DSS. As illustrated in Fig. 7, the specific D-T subnetwork is composed of 194 known DTIs and 14 predicted DTIs connecting 9 compounds with 62 target proteins. The 9 ingredients are connected to multiple targets with the average degree of 3.4. Among them, gallic acid (CID370) has the highest number of target connections, binding with 37 known targets and 3 predicted targets, which indicate its potential multi-target anti-AD mechanism. For example, MMP7, one of the three predicted targets, has been demonstrated to play an important role in the defensive mechanism against the aggregation of Aβ1–42, giving rise to the pathology of AD . Remarkably, the 9 ingredients in DSS share multiple overlapped protein targets. Five proteins, including TSHR, TOP1, CYP3A4, LMNA and MAPT, are simultaneously targeted by all the 9 ingredients, suggesting that DSS may mainly exert its anti-AD effect through the synergistic effect of these components. Taking MAPT as an example, published studies have revealed that MAPT gene polymorphisms could increase the AD risk, which is highly related with the deposition of Aβ proteins [68,69,70].
Overall, the prescription simplification strategy applied here provides a new perspective for uncovering the pivotal constituents in DSS and their specific molecular mechanism for treatment of AD.
To date, the prevention and treatment of AD remains a huge challenge. In last two decades, although more than 400 clinical trials have been conducted for AD treatment, over 99% failed in clinical trials owing to the adverse events and limited benefits . Therefore, it is urgent to find a new and effective therapy for AD from a broader perspective. TCM prescriptions emphasize the rebalancing of interactions among various disease-related factors within the abnormal body. As one of the representative prescriptions, DSS has been demonstrated to be effective for improving cognitive function of AD patients . However, since TCM prescriptions contain complex chemical components interacted with multiple targets, the specific pharmacological mechanisms of DSS against AD are still difficult to illustrate.
In this study, we developed an integrative systems pharmacology framework to decipher the MOAs of DSS against AD. We first performed an overlapped analysis among the herbal ingredients as well as potential targets and constructed a global drug-target network of DSS. We further explored the potential MOAs of DSS for AD treatment through multiple enrichment analyses, including target-function modules analysis, GO terms annotation, KEGG pathway enrichment, and AD-integrated pathway. Finally, we proposed a prescription simplification strategy, which integrates HPLC analysis results, drug-likeness analysis, pharmacological properties prediction, and network-based statistical model, to uncover the key anti-AD compounds in DSS. We highlighted 9 key ingredients in DSS that play synergistic role against AD in the subnetwork.
Systems pharmacology analysis suggests that DSS participates in the regulation of several important biological pathways related to AD pathogenesis, such as the oxidative stress, inflammatory reaction processes and neurotrophin signaling pathway (Fig. 4 and 5a). Recently, Huang et al.  found that DSS extract prominently lessened the abnormal activity of ROS and malondialdehyde (MDA) level, and increased the intracellular superoxide dismutase (SOD) and glutathione (GSH) level in APP/PS1 mice. Moreover, following treatment with DSS, the levels of inflammatory factors, including leukotriene B4 (LTB4), prostaglandin E2 (PGE2) and thromboxane B2 (TXB2), were significantly decreased in hippocampus and cortex of APP/PS1 mice . These findings suggested that DSS play a positive effective role in ameliorating oxidative stress and neuroinflammation and finally ameliorating cognition deficits in APP/PS1 mice , which are consistent with our predictions in this study.
However, several shortcomings should be recognized in the presented study. First, although we have integrated a large number of protein targets from published literatures, publicly available databases and network-based inference method, the incompleteness of current DTI networks may still exist. Second, as the intrinsic interactions and level of expression of multiple ingredients in the mixture of different plants are complicated, current study could not entirely reflect the actual functional actions of DSS in the human body. Finally, although recent experimental study had validated part of the predicted biological processes regulating by DSS, such as oxidative stress and inflammation , the remaining predicted targets and anti-AD mechanisms of DSS are necessary to be further validated by more comprehensive and in-depth wet-lab experiment in the future.
The present study provides potential strategies for exploring the therapeutic mechanism of DSS against AD, which also indicates that systems pharmacology could be an effective infrastructure to decipher the compatibility and MOAs of the complex components in TCM prescription.
Availability of data and materials
The data is available and will be provided upon request to the corresponding author.
Traditional Chinese medicine
Human intestinal absorption
Paeoniae Radix Alba, BaiShao
Angelica Sinensis Radix, DangGui
Atractylodis Macrocephalae Rhizoma, BaiZhu
Chuanxiong Rhizoma, ChuanXiong
Alismatis Rhizoma, ZeXie
Randomized controlled trials
Mini-Mental State Examination
Activities of daily living
Regional cerebral blood flow
Mechanism of actions
High-performance liquid chromatography
Comparative Toxicogenomics Database
Human Gene Mutation Database
Herbal Ingredients Target Database
Traditional Chinese Medicine Integrated Database
A balanced substructure-drug-target network-based inference
Peroxisome proliferator-activated receptor
Reactive oxygen species
Hypoxia-inducible factor 1
Mitogen-activated protein kinase
Nuclear factor-kappa B
Brain-derived neurotrophic factor
Masters CL, Bateman R, Blennow K, Rowe CC, Sperling RA, Cummings JL. Alzheimer's disease. Nat Rev Dis Primers. 2015;1:15056.
Hardy J, Selkoe DJ. The amyloid hypothesis of Alzheimer's disease: progress and problems on the road to therapeutics. Science. 2002;297(5580):353–6.
Braak H, Braak E. Neuropathological stageing of Alzheimer-related changes. Acta Neuropathol. 1991;82(4):239–59.
World Health Organization Dementia factsheet http://www.who.int/mediacentre/factsheets/fs362/en/ 2019.
Fang JS, Li YJ, Liu R, Pang XC, Li C, Yang R, et al. Discovery of multitarget-directed ligands against Alzheimer's disease through systematic prediction of chemical-protein interactions. J Chem Inf Model. 2015;55(1):149–64.
Cai HB, Luo YX, Yan X, Ding P, Huang YJ, Fang SH, et al. The mechanisms of Bushen-Yizhi formula as a therapeutic agent against Alzheimer's disease. Sci Rep. 2018;8(1):3104.
Jiang WY. Therapeutic wisdom in traditional Chinese medicine: a perspective from modern science. Discov Med. 2005;5(29):455–61.
Wang TDY, Wu ZR, Sun LX, Li WH, Liu GX, Tang Y. A computational systems pharmacology approach to investigate molecular mechanisms of herbal formula Tian-ma-Gou-Teng-Yin for treatment of Alzheimer's disease. Front Pharmacol. 2018;9:668.
Jeon S, Hur J, Jeong HJ, Koo BS, Pak SC. SuHeXiang wan essential oil alleviates amyloid beta induced memory impairment through inhibition of tau protein phosphorylation in mice. Am J Chin Med. 2011;39(5):917–32.
Gao AS, Du H, Gao Q, Wang YJ, Wang XY, Liu JW, et al. Tong Luo Jiu Nao, a Chinese medicine formula, reduces inflammatory stress in a mouse model of Alzheimer's disease. Am J Alzheimers Dis Other Dement. 2016;31(5):413–21.
Kim Y, Cho SH. Danggui-Shaoyao-san for dementia: a PRISMA-compliant systematic review and meta-analysis. Medicine. 2020;99(4):e18507.
Matsuoka T, Narumoto J, Shibata K, Okamura A, Taniguchi S, Kitabayashi Y, et al. Effect of Toki-shakuyaku-san on regional cerebral blood flow in patients with mild cognitive impairment and Alzheimer's disease. Evid Based Complement Alternat Med. 2012;2012:245091.
Huang Y, Hu ZY, Yuan H, Shu L, Liu G, Qiao SY, et al. Danggui-Shaoyao-san improves learning and memory in female SAMP8 via modulation of estradiol. Evid Based Complement Alternat Med. 2014;2014:327294.
Egashira N, Iwasaki K, Akiyoshi Y, Takagaki Y, Hatip-Al-Khatib I, Mishima K, et al. Protective effect of Toki-shakuyaku-san on amyloid beta25-35-induced neuronal damage in cultured rat cortical neurons. Phytother Res. 2005;19(5):450–3.
Guo PF, Cai CP, Wu XQ, Fan XD, Huang W, Zhou JW, et al. An insight into the molecular mechanism of Berberine towards multiple Cancer types through systems pharmacology. Front Pharmacol. 2019;10:857.
Cai FF, Zhou WJ, Wu R, Su SB. Systems biology approaches in the study of Chinese herbal formulae. Chin Med. 2018;13:65.
Li S, Zhang B, Zhang NB. Network target for screening synergistic drug combinations with application to traditional Chinese medicine. BMC Syst Biol. 2011;5(Suppl 1(Suppl 1)):S10.
Wu QH, Cai CP, Guo PF, Chen ML, Wu XQ, Zhou JW, et al. Identification and Mechanism Exploration of Hepatotoxic Ingredients in Traditional Chinese Medicine. Front Pharmacol. 2019;10:458.
Fang JS, Cai CP, Wang Q, Lin P, Zhao Z, Cheng FX. Systems pharmacology-based discovery of natural products for precision oncology through targeting Cancer mutated genes. CPT Pharmacometrics Syst Pharmacol. 2017a;6(3):177–87.
Davis AP, Grondin CJ, Johnson RJ, Sciaky D, McMorran R, Wiegers J, et al. The comparative Toxicogenomics database: update 2019. Nucleic Acids Res. 2019;47(D1):D948–54.
Pinero J, Bravo À, Queralt-Rosinach N, Gutiérrez-Sacristán A, Deu-Pons J, Centeno E, et al. DisGeNET: a comprehensive platform integrating information on human disease-associated genes and variants. Nucleic Acids Res. 2017;45(Database issue):D833–9.
Welter D, MacArthur J, Morales J, Burdett T, Hall P, Junkins H, et al. The NHGRI GWAS catalog, a curated resource of SNP-trait associations. Nucleic Acids Res. 2014;42(Database issue):1001–6.
Stenson PD, Mort M, Ball EV, Evans K, Hayden M, Heywood S, et al. The human gene mutation database: towards a comprehensive repository of inherited mutation data for medical research, genetic diagnosis and next-generation sequencing studies. Hum Genet. 2017;136(6):665–77.
Huang YJ, Fang JS, Lu WQ, Wang ZH, Wang Q, Hou Y, et al. A systems pharmacology approach uncovers Wogonoside as an angiogenesis inhibitor of triple-negative breast Cancer by targeting hedgehog signaling. Cell Chem Biol. 2019;26(8):1143–58.
Gaulton A, Hersey A, Nowotka M, Bento AP, Chambers J, Mendez D, et al. The ChEMBL database in 2017. Nucleic Acids Res. 2017;45(Database issue):D945–54.
Gilson MK, Liu T, Baitaluk M, Nicola G, Hwang L, Chong J. BindingDB in 2015: a public database for medicinal chemistry, computational chemistry and systems pharmacology. Nucleic Acids Res. 2016;44(Database issue):D1045–53.
Szklarczyk D, Santos A, Mering CV, Jensen LJ, Bork P, Kuhn M. STITCH 5: augmenting protein–chemical interaction networks with tissue and affinity data. Nucleic Acids Res. 2016;44(Database issue):D380–4.
Ye H, Ye L, Kang H, Zhang D, Tao L, Tang K, et al. HIT: linking herbal active ingredients to targets. Nucleic Acids Res. 2011;39(Database issue):D1055.
Huang L, Xie DL, Yu YR, Liu HL, Shi Y, Shi TL, et al. TCMID 2.0: a comprehensive resource for TCM. Nucleic Acids Res. 2018;46(D1):D1117–20.
Wu ZR, Lu WQ, Wu D, Luo AQ, Bian HP, Li J, et al. In silico prediction of chemical mechanism-of-action via an improved network-based inference method. Br J Pharmacol. 2016;173(23):3372–85.
Fang JS, Wu ZR, Cai CP, Wang Q, Tang Y, Cheng FX. Quantitative and systems pharmacology. 1. In Silico prediction of drug-target interaction of natural products to enable of new targeted Cancer therapy. J Chem Inf Model. 2017;57(11):2657–71.
Yap CW. PaDEL-descriptor: an open source software to calculate molecular descriptors and fingerprints. J Comput Chem. 2011;32(7):1466–74.
Xu F, Peng DY, Tao CL, Yin DK, Kou JP, Zhu DN, et al. Anti-depression effects of Danggui-Shaoyao-san, a fixed combination of traditional Chinese medicine, on depression model in mice and rats. Phytomedicine. 2011;18(13):1130–6.
Chen L, Qi J, Chang YX, Zhu D, Yu B. Identification and determination of the major constituents in traditional Chinese medicinal formula Danggui-Shaoyao-san by HPLC-DAD-ESI-MS/MS. J Pharm Biomed Anal. 2009;50(2):127–37.
Lipinski CA, Lombardo F, Dominy BW, Feeney PJ. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv Drug Deliv Rev. 2001;46(1–3):3–26.
Shen J, Cheng FX, Xu Y, Li WQ, Tang Y. Estimation of ADME properties with substructure pattern recognition. J Chem Inf Model. 2010;50(6):1034–41.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser A. 1995;57(1):289–300.
Li S, Zhang B. Traditional Chinese medicine network pharmacology: theory, methodology and application. Chin J Nat Med. 2013;11(2):110–20.
Qiu J. Traditional medicine: a culture in the balance. Nature. 2007;448(7150):126–8.
Luo R, Su LY, Li G, Yang J, Liu Q, Yang LX, et al. Activation of PPARA-mediated autophagy reduces Alzheimer disease-like pathology and cognitive decline in a murine model. Autophagy. 2020;16(1):52–69.
Sabogal-Guáqueta AM, Muñoz-Manco JI, Ramírez-Pineda JR, Lamprea-Rodriguez M, Osorio E, Cardona-Gómez GP. The flavonoid quercetin ameliorates Alzheimer's disease pathology and protects cognitive and emotional function in aged triple transgenic Alzheimer's disease model mice. Neuropharmacology. 2015;93:134–45.
Balez R, Steiner N, Engel M, Muñoz SS, Lum JS, Wu Y, et al. Neuroprotective effects of apigenin against inflammation, neuronal excitability and apoptosis in an induced pluripotent stem cell model of Alzheimer's disease. Sci Rep. 2016;6:31450.
Trojanowski JQ, Mawal-Dewan M, Schmidt ML, Martin J, Lee VM. Localization of the mitogen activated protein kinase ERK2 in Alzheimer's disease neurofibrillary tangles and senile plaque neurites. Brain Res. 1993;618(2):333–7.
Lukiw WJ, Cui JG, Marcheselli VL, Bodker M, Botkjaer A, Gotlinger K, et al. A role for docosahexaenoic acid-derived neuroprotectin D1 in neural cell survival and Alzheimer disease. J Clin Invest. 2005;115(10):2774–83.
Huang JW, Wang XY, Xie LY, Wu MG, Zhao W, Zhang YB, et al. Extract of Danggui-Shaoyao-san ameliorates cognition deficits by regulating DHA metabolism in APP/PS1 mice. J Ethnopharmacol. 2020;253:112673.
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, et al. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 2003;4(5):P3.
Tönnies E, Trushina E. Oxidative stress, synaptic dysfunction, and Alzheimer's disease. J Alzheimers Dis. 2017;57(4):1105–21.
Shadfar S, Hwang CJ, Lim MS, Choi DY, Hong JT. Involvement of inflammation in Alzheimer's disease pathogenesis and therapeutic potential of anti-inflammatory agents. Arch Pharm Res. 2015;38(12):2106–19.
Yu H, Yuan B, Chu QB, Wang CY, Bi H. Protective roles of isoastilbin against Alzheimer's disease via Nrf2-mediated antioxidation and anti-apoptosis. Int J Mol Med. 2019;43(3):1406–16.
Supuran CT. Carbonic anhydrases: novel therapeutic applications for inhibitors and activators. Nat Rev Drug Discov. 2008;7(2):168–81.
Cissé M, Duplan E, Lorivel T, Dunys J, Bauer C, Meckler X, et al. The transcription factor XBP1s restores hippocampal synaptic plasticity and memory by control of the Kalirin-7 pathway in Alzheimer model. Mol Psychiatry. 2017;22(11):1562–75.
Li XX, Xu X, Wang J, Yu H, Wang X, Yang HJ, et al. A system-level investigation into the mechanisms of Chinese traditional medicine: compound Danshen formula for cardiovascular disease treatment. PLoS One. 2012;7(9):e43918.
Hsu W-CJ, Wildburger NC, Haidacher SJ, Nenov MN, Folorunso O, Singh AK, et al. PPARgamma agonists rescue increased phosphorylation of FGF14 at S226 in the Tg2576 mouse model of Alzheimer's disease. Exp Neurol. 2017;295:1–17.
Nho K, Kueider-Paisley A, MahmoudianDehkordi S, Arnold M, Risacher SL, Louie G, et al. Altered bile acid profile in mild cognitive impairment and Alzheimer's disease: relationship to neuroimaging and CSF biomarkers. Alzheimers Dement. 2019;15(2):232–44.
Iyalomhe O, Swierczek S, Enwerem N, Chen Y, Adedeji MO, Allard J, et al. The role of hypoxia-inducible factor 1 in mild cognitive impairment. Cell Mol Neurobiol. 2017;37(6):969–77.
Calsolaro V, Edison P. Neuroinflammation in Alzheimer's disease: current evidence and future directions. Alzheimers Dement. 2016;12(6):719–32.
Luo YX, Li DL, Liao YF, Cai CP, Wu QH, Ke HZ, et al. Systems pharmacology approach to investigate the mechanism of Kai-Xin-san in Alzheimer's disease. Front Pharmacol. 2020;11:381.
Lee JK, Kim NJ. Recent advances in the inhibition of p38 MAPK as a potential strategy for the treatment of Alzheimer's disease. Molecules. 2017;22(8):1287.
Chen CH, Zhou W, Liu S, Deng Y, Cai F, Tone M, et al. Increased NF-κB signalling up-regulates BACE1 expression and its therapeutic potential in Alzheimer's disease. Int J Neuropsychopharmacol. 2012;15(1):77–90.
Sivandzade F, Prasad S, Bhalerao A, Cucullo L. NRF2 and NF-қB interplay in cerebrovascular and neurodegenerative disorders: molecular mechanisms and possible therapeutic approaches. Redox Biol. 2019;21:101059.
Ding RR, Chen W, Guo CY, Liao WT, Yang X, Liao FE, et al. Dangguishaoyao-san attenuates LPS-induced neuroinflammation via the TLRs/NF-κB signaling pathway. Biomed Pharmacother. 2018;105:187–94.
Chen XQ, Sawa M, Mobley WC. Dysregulation of neurotrophin signaling in the pathogenesis of Alzheimer disease and of Alzheimer disease in Down syndrome. Free Radic Biol Med. 2018;114:52–61.
Zhang L, Fang Y, Lian YJ, Chen Y, Wu TW, Zheng YK, et al. Brain-derived neurotrophic factor ameliorates learning deficits in a rat model of Alzheimer's disease induced by aβ1-42. PLoS One. 2015;10(4):e0122415.
Yan YH, Li SH, Gao Z, Zou SF, Li HY, Tao ZY, et al. Neurotrophin-3 promotes proliferation and cholinergic neuronal differentiation of bone marrow- derived neural stem cells via notch signaling pathway. Life Sci. 2016;166:131–8.
Lu F, Li X, Li W, Wei K, Yao Y, Zhang QL, et al. Tetramethylpyrazine reverses intracerebroventricular streptozotocin-induced memory deficits by inhibiting GSK-3β. Acta Biochim Biophys Sin. 2017;49(8):722–8.
Kuang X, Zhou HJ, Thorne AH, Chen XN, Li LJ, Du JR. Neuroprotective effect of Ligustilide through induction of α-Secretase processing of both APP and Klotho in a mouse model of Alzheimer's disease. Front Aging Neurosci. 2017;9:353.
Taniguchi M, Matsuura K, Nakamura R, Kojima A, Konishi M, Akizawa T. MMP-7 cleaves amyloid β fragment peptides and copper ion inhibits the degradation. Biometals. 2017;30(5):797–807.
Kakati T, Kashyap H, Bhattacharyya DK. THD-module extractor: an application for CEN module extraction and interesting gene identification for Alzheimer's disease. Sci Rep. 2016;6:38046.
Zhou F, Wang D. The associations between the MAPT polymorphisms and Alzheimer's disease risk: a meta-analysis. Oncotarget. 2017;8(26):43506–20.
Desikan RS, Schork AJ, Wang Y, Witoelar A, Sharma M, McEvoy LK, et al. Genetic overlap between Alzheimer's disease and Parkinson's disease at the MAPT locus. Mol Psychiatry. 2015;20(12):1588–95.
Dey A, Bhattacharya R, Mukherjee A, Pandey DK. Natural products against Alzheimer's disease: Pharmaco-therapeutics and biotechnological interventions. Biotechnol Adv. 2017;35(2):178–216.
This work was supported by the Provincial Key R&D Program of Hainan (No. ZDYF2019196), National Natural Science Foundation of China (No. 81903912, 81674040 and 81673627), the Youth Scientific Research Training Project of GZUCM (No. 2019QNPY05), China Postdoctoral Science Foundation funded project (No. 2019 M662878) and the Guangzhou Science Technology and Innovation Commission Technology Research Projects (No. 201805010005). The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Detailed information of the 299 genes related with AD.
Detailed structural information of the 1,042 herbal ingredients in DSS.
The top 20 putative targets for each compound in DSS with known targets.
The 79 common targets exist in all the six herbs of DSS.
Detail information for the drug-target network of DSS.
Detailed information of GO terms enrichment analysis data of targets with degree equal or larger than 40.
Detailed information of KEGG pathway enrichment analysis data of targets with degree equal or larger than 40.
Detailed information of the integrated pathway on DSS for treatment of AD.
About this article
Cite this article
Wu, Q., Chen, Y., Gu, Y. et al. Systems pharmacology-based approach to investigate the mechanisms of Danggui-Shaoyao-san prescription for treatment of Alzheimer’s disease. BMC Complement Med Ther 20, 282 (2020). https://doi.org/10.1186/s12906-020-03066-4
- Systems pharmacology
- Alzheimer’s disease
- Mechanism of action