Abstract Purpose This study was designed to explore the underlying mechanism of action for a Fuzheng Quxie Decoction (FQD) in Alzheimer’s disease (AD), to validate its neuroprotective effects, and to provide experimental support for its predicted mechanism of action. Methods An integrative approach to network pharmacology was performed to predict the mechanism of action for treatment of AD with FQD. The predicted mechanism was validated in SAMP8 mice. Results With predicted putative FQD targets and a collection of AD-related genes, 245 possible regulatory targets of FQD were identified for the treatment of AD. Pathway-enrichment analysis for the possible regulatory targets indicated that vascular endothelial growth factor (VEGF) and VEGF-receptor signaling were pivotal in the treatment of AD with FQD. In vivo experiments confirmed the neuroprotective effect and the predicted mechanism of action for treatment of AD with FQD. Conclusion This study contributes to an understanding of the neuroprotective effect of FQD and its potential mechanism of action for the treatment of AD. Keywords: Alzheimer’s disease, Fuzheng Quxie Decoction, network pharmacology, microcerebrovascular, vascular endothelial growth factor, hypoxia-induced factor 1α Introduction Alzheimer’s disease (AD) is a progressive neurodegenerative disease in aging populations that manifests as short- and long-term memory loss, learning and memory retention deficits, cognitive impairment, and behavioral and mental disability.[41]^1 By the year 2030, the number of AD patients is projected to reach 70 million,[42]^2 which represents a huge burden to society and the economy. Current commonly used drugs for treatment of AD, eg, donepezil and memantine, target the cholinergic system and N-methyl-d-aspartate receptors. They can only relieve symptoms temporarily, rather than halt or slow disease progression.[43]^3 The pathogenesis of AD is intricate and not clear. Although senile plaques and neurofibrillary tangles, mainly composed of β-amyloid (Aβ) and tau, respectively,[44]^4^,[45]^5 have been found to be the primary pathological features in the brains of AD patients, the development of different drugs targeting Aβ or tau has been unsuccessful.[46]^6^–[47]^11 Other proposed hypotheses for AD pathogenesis, such as oxidative stress and metalion accumulation, have been demonstrated to play a key role in blood–brain barrier (BBB) disruption, with BBB disruption possibly upstream of Aβ deposition and/or tau hyperphosphorylation in the pathogenic cascade that results in AD onset.[48]^12^,[49]^13 Based on traditional Chinese medicine (TCM) theory, AD is believed to be caused by qi deficiency, blood stasis, and toxins. Qi deficiency and blood stasis are very common in the elderly. This condition results in metabolic waste that cannot be timely removed and with stasis is transformed into toxins. Qi deficiency, blood stasis, and toxins can be respectively interpreted into modern concepts of aging symptoms, hyperlipemia, and Aβ and tau-protein deposition. Accumulating data have demonstrated that medicinal Chinese herbs or herbal formulas, such as Radix Ginseng,[50]^14^,[51]^15 Rhizoma Ligustici Chuanxiong, Rhizoma Coptidis,[52]^16^,[53]^17 Qifu Yin,[54]^18 and Huannao Yicong Decoction,[55]^19^–[56]^22 tonify qi, dissolve blood stasis, and remove toxins to delay onset and progression of AD, and these have been widely used in clinics to treat AD in China. The herbal formula Fuzheng Quxie Decoction (FQD) was prescribed to patients with AD or cognitive impairment in accordance with TCM theory by the leader of this study – Hao Li. FQD has been shown to alleviate AD symptoms and slow disease progression.[57]^16 The constitutive herbs of FQD are Radix Ginseng, Rhizoma Ligustici Chuanxiong, and Rhizoma Coptidis and are used in different TCM areas. FQD has been used widely in clinics to treat AD in China. Compared with other herbal formulas for AD treatment, FQD has the advantages of simple composition by virtue of optimized compatibility of herbs, economy, and convenience of preparation, with patient compliance comparatively high. Since a Chinese herbal formula consists of multiple compounds targeting different molecules, it is difficult to elucidate the underlying mechanistic basis for therapeutic disease treatment. However, it has become possible to interpret a Chinese herbal formula with the development of technologies and methodologies, such as network pharmacological approach which possesses the characteristics of holism similar to TCM. By means of network pharmacology, active compounds in a Chinese herbal formula can be identified and their potential targets defined. As such, the identification of the mechanistic basis for disease treatment using the Gene Ontology (GO) database and pathway-enrichment analysis is possible in the foreseeable future.[58]^23 Our previous study indicated that FQD was able to reduce tau hyperphosphorylation,[59]^24 but the underlying mechanistic basis was not determined. Therefore, in this study, active components of FQD were screened and possible targets and regulatory pathways predicted by network pharmacology. FQD-predicted mechanistic basis of action was validated in vivo, with results showing that FQD has a neuroprotective role in the treatment of AD by regulating predicted signaling pathways. Methods Prediction of FQD mechanism of action for AD treatment Active-compound screening The ingredients in three herbs contained in FQD were retrieved from the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP; [60]http://lsp.nwu.edu.cn/tcmsp.php),[61]^25 the TCM Database@ Taiwan ([62]http://tcm.cmu.edu.tw),[63]^26 the Chinese Academy of Sciences Chemistry Database ([64]http://www.organchem.csdb.cn/scdb/default.asp), and the Traditional Chinese Medicines Integrated Database (TCMID).[65]^27 All the ingredients were combined to construct an FQD-ingredient database after removal of duplications. Since glycosides can be deglycosylated after oral administration,[66]^28 their deglycosylated forms were included in the database. Active FQD compounds were screened by absorption, distribution, metabolism, and excretion (ADME) parameters, including oral bioavailability (OB), drug likeness (DL), and BBB permeability (BBBp). OB represents the percentage of an orally administered drug that distributes via systemic circulation.[67]^29 DL is an index of the structural similarities between an unknown compound and the compounds with the known structures in the Drug Bank database. Both OB and DL are key parameters for drug discovery. BBBp is the ability of a compound to penetrate the BBB. Compounds with BBBp <−0.3 do not penetrate the BBB, those with BBBp −0.3 to 0.3 have moderate penetration, and those with BBBp >0.3 have strong penetration.[68]^30 The BBBp of compounds in FQD was assessed because access to the brain is traditionally essential for treatment of AD. In this study, OB, DL, and BBBp were evaluated with the TCMSP calculation system. Active compounds were identified that had thresholds of OB ≥30%, DL ≥0.18, and BBBp ≥−0.3, as described by Liu et al.[69]^30 Neuroprotective compounds that had been previously reported to exert an effect on the central nervous system were also considered active compounds, even if they did not meet the aforementioned criteria. Target prediction After active compounds had been identified, their targets were subsequently predicted by the online target-prediction systems TCMSP, TCMID, SwissTargetPrediction ([70]http://www.swisstargetprediction.ch),[71]^31 and STITCH ([72]http://stitch.embl.de).[73]^32 Each of these is mature and reliable, with integrated and optimized algorithms. Possible targets were calculated and predicted by uploading the structure of the compounds. Targets with prediction scores >0.7 were considered putative targets. The names of the targets were converted into gene symbols via UniProt ([74]https://www.uniprot.org)[75]^33 to smooth subsequent network construction and analysis. Collection of AD-related genes There are several online databases of disease-related genes based on the published literature. In this study, AD-related genes were collected from Online Mendelian Inheritance in Man ([76]https://omim.org),[77]^34 MalaCards ([78]http://www.malacards.org),[79]^35 and the Kyoto Encyclopedia of Genes and Genomes (KEGG, [80]http://www.genome.jp/kegg).[81]^36 After redundancy removal, all genes retrieved were combined into a database of AD-related genes. Construction of protein–protein interaction network Putative targets of FQD and AD-related genes were intersected to identify shared genes, which were defined as direct regulatory targets (DRTs) of FQD for the treatment of AD. A protein–protein interaction (PPI) network of DRTs was constructed with the BisoGenet 3.0 plug-in for Cytoscape 3.6.0[82]^37 based on the online integrated PPI databases.[83]^38 In-house PPI databases, including the IntAct Molecular Interaction Database, Biological General Repository for Interaction Datasets, Database of Interacting Proteins, Human Protein Reference Database, Molecular Interaction Database, and Binding Database, were selected to predict interactions among proteins. The network was expressed with nodes and edges, which represented genes and their interactions, respectively. The PPI network was analyzed by one of its topological features – degree centrality (DC),[84]^39 which was calculated by the Centiscape 2.2.0 plug-in for Cytoscape. Genes represented by nodes with DC values greater than twice the average were considered to play a more important role in the network. These were chosen as possible regulatory targets (PRTs) of FQD for the treatment of AD. Enrichment analysis of PRTs In order to gain insight into the function of PRTs, their genes were annotated by GO enrichment analysis with FunRich 3.0,[85]^40 which describes the functions of genes based on biological process, cellular component, and molecular function.[86]^41 The most significantly enriched terms of each category were visualized with OmicShare ([87]http://www.omicshare.com). Further, to understand PRT-signaling pathways, pathway-enrichment analysis was performed and visualized with FunRich 3.0. Enriched signaling pathways were considered those that FQD primarily regulated in the treatment of AD. Predicted mechanisms were validated in vivo. Experimental validation Drugs and reagents for animal experimentation The three medicinal herbs of FQD were purchased from Hebei Shennong (Beijing, China) Pharmaceutical. Extract preparation was carried out by the Department of Pharmaceutical Preparation at Xiyuan Hospital. Briefly, all herbs were combined and boiled for 30 minutes in distilled water, after which the liquid was collected. This was repeated and the liquids from each boiling mixed. An FQD extract was obtained by reduced-pressure distillation of the mixed liquid. There was 3.71 g of crude drug per gram of the extract. Components were detected by HPLC as described by Yang et al.[88]^24 Memantine was obtained from the National Institute for Food and Drug Control (batch number H20120268; Beijing, China), produced by Lundbeck (Copenhagen, Denmark). When used, the extract and memantine were dissolved separately in distilled water at intended concentrations. Rabbit polyclonal anti-vascular endothelial growth factor (VEGF) antibody and rabbit polyclonal anti-hypoxia-induced factor 1α (HIF1α) antibody were provided by Abcam (Cambridge, UK). Mouse monoclonal anti-β-actin antibody was provided by ImmunoWay Biotechnology (Plano, TX, USA). Goat antirabbit/mouse IgG antibodies conjugated with HRP were provided by TDY Biotech (Beijing, China). The streptavidin–biotin complex kit, diaminobenzidine-substrate solution and enhanced chemiluminescence reagent were purchased from Thermo Fisher Scientific (Waltham, MA, USA). RIPA lysis buffer and a Micro BCA protein-assay kit were purchased from Sinoble Biotechnology Center (Beijing, China). Saline, paraformaldehyde, paraffin, H&E kit, cresyl violet, PBS, BSA, and Tris-buffered saline Tween (TBST) were purchased from Beyotime Biotechnology (Haimen, China). Animal grouping and drug administration Animal experiments were carried out on male SAMP8 mice. SAMP8 mice spontaneously develop accelerated senescence,[89]^42 showing symptomatic characteristics typical of AD, such as early memory and learning deficits,[90]^43^,[91]^44 disturbed BBB function, accumulation of Aβ,[92]^45 hyperphosphorylation of tau,[93]^46 and marked oxidant-stress pathology.[94]^47 These mice exhibit a similar pathological cascade to AD patients. Successful treatment of AD should improve memory, alleviate BBB dysfunction, and reduce oxidant stress in SAMP8 mice.[95]^48 In this study, FQD was predicted to regulate the VEGF pathway (see Results section), which is crucial for normal BBB function and prevention of the flux of Aβ into the brain. For these reasons, SAMP8 mice were chosen as an appropriate model of AD. A total of 40 male SAMP8 mice 3 months of age were housed in normal cages and kept in a controlled environment (22°C, 50% humidity, 12-hour light/dark cycle) with free access to water and food. After 1 week of adaptation, the mice were divided into four groups randomly: a memantine group, FQD high-dose group, FQD low-dose group, and a model group. Mice were treated orally once a day for 12 consecutive weeks with memantine at 2.6 mg/kg body weight, FQD at 2.6 g crude drug/kg body weight, FQD at 1.3 g crude drug/kg body weight, or distilled water of the same volume, respectively. Ten SAMR1 mice of the same age and genetic background as the SAMP8 mice served as a control group and were kept under the same conditions and given distilled water of the same volume. Two mice died when given the memantine gavage, and two mice died when given the high FQD dose due to operational mistakes. Animal experimentation was approved by the ethics committee of Xiyuan Hospital, China Academy of Chinese Medical Sciences (permit CACMS/20141220×21). All animal procedures were carried out in accordance with the Guide for the Care and Use of Laboratory Animals (8th edition) of the National Research Council. Morris water-maze test The Morris water maze (MWM) test was conducted after 12 weeks of treatment.[96]^49 The MWM apparatus, provided by Zhongshi Dichuang Technology (Beijing, China), consists of a circular water tank 120 cm in diameter and 50 cm deep, a platform 10 cm in diameter and 20 cm in height, and a computerized autoanalyzing camera system above the tank. The tank was divided into four quadrants, and the platform was placed at the very center of the fourth quadrant and 1 cm underwater. The water temperature was kept at 21°C±1°C. Four signs of different shapes (square, triangle, circle, and rhombus) were placed on the inner wall of each quadrant. Ink was used to stain the water black, ensuring that mice could not see the platform. For the first 5 days, each mouse was placed into water at different entry points facing the tank wall (as shown in [97]Figure 1) four times each day and allowed to search for the platform as training. The time required for each mouse to find the platform (latency to escape) was recorded automatically by the MWM software. If a mouse failed to reach the platform within 60 seconds, it was guided to the platform. All mice were allowed to stay on the platform for 10 seconds, after which they were dried and returned to their cages. Figure 1. Figure 1 [98]Open in a new tab The position of platform and entry points in the MWM test. Abbreviation: MWM, Morris water maze. On the sixth day, the platform was removed and spatial-probe trail analysis conducted. Each mouse was placed into water at a novel point facing the tank wall. During a 60-second interval, mice were permitted to swim freely. The time and distance spent in the tank and in the fourth quadrant were recorded as indices. Passive-avoidance test On the third day after the MWM test had been completed, a passive-avoidance test was performed. The experimental apparatus, provided by Zhongshi Dichuang Technology, consists of a well-illuminated square compartment and a dark square compartment (15×15×15 cm^3) with a removable roof and metal-grid floor for the administration of electric shocks. Between the two compartments there was a guillotine door. On the first day, each mouse was put into the illuminated compartment with free access to the dark compartment for 3 minutes to adapt to the environment. When the training began, the mice were given a weak electric shock to the foot (0.5 mA) as long as they were in the dark compartment. Then, the mice were returned to their cages. In order to test their memory, 24 hours later the same procedure was conducted without adaptation. The avoidance latency of mice entering the dark compartment for the first time was recorded and the number of times mice entered the dark compartment recorded as error times. If a mouse did not enter the dark compartment within the 5-minute test period, the test was terminated and avoidance latency recorded as 300 seconds. The training and test were conducted in an isolated experimental room with dim light. H&E and Nissl staining After the behavioral study, all the mice were perfused transcardially with normal saline under deep anesthesia. Brains of the mice were removed after perfusion. Cortices were dissected from five of the left hemispheres and stored at −80°C until used for Western blot (WB) analysis. The right hemispheres and another five of the left hemispheres were postfixed with 4% paraformaldehyde overnight. All dissections were completed on ice. The fixed hemispheres were paraffin-embedded and sliced into 20 μm-thick coronal sections, which were used for H&E staining, Nissl staining and immunohistochemistry. Five sections from each group were stained with H&E to observe the morphology of the cerebral microvascular endothelial cells microscopically. Another five sections were stained with 0.75% cresyl violet for neuron analysis. Neurons in the hippocampus CA1 area were counted microscopically. Immunohistochemistry A series of sections containing the hippocampus and cortex from five of the left hemispheres of each group were deparaffinized, blocked with sheep serum, and reacted with anti-VEGF antibody and the biotinylated secondary antibody (anti-IgG antibody), with PBS washes between each step. Then, the prepared streptavidin–biotin complex solution was added, followed by diaminobenzidine-substrate solution. After tissue samples had been stained and mounted, three fields for each slide were assessed and the number of positive cells determined microscopically. Western blot analysis The prepared cortices dissected from the left hemispheres of each group were lysed in RIPA buffer and then centrifuged at 13,000 rpm for 20 minutes at 4°C in a refrigerated centrifuge (Thermo Fisher Scientific). The supernatants were extracted and protein concentrations determined with the Micro BCA kit. The same amounts of protein were resolved and separated on 8% SDS-PAGE and transferred to polyvinylidene difluoride (PVDF) membranes (Thermo Fisher Scientific). After being blocked with 3% BSA-TBST for 30 minutes, PVDF membranes were incubated overnight at 4°C with polyclonal rabbit anti-VEGF antibody and anti-HIF1α antibody (1:1,000) in 3% BSA-TBST. The PVDF membranes were washed with TBST and incubated for 40 minutes at room temperature with anti-IgG antibody conjugated with HRP (1:10,000) dissolved in 1% BSA-TBST. Protein bands were detected with enhanced-chemiluminescence reagents and visualized with an imaging system (MicroChemi). Visualized protein bands were quantified with Image-Pro Plus 5.0 software. Statistical analysis All data are expressed as mean ± SEM and were analyzed using GraphPad Prism 5.0 (GraphPad software, San Diego, CA, USA). One-way ANOVA was used for intergroup comparisons. Least significant difference and Dunnett’s C test were used to assess the heterogeneity of variance. Differences were considered significant when P<0.05. Results Mechanism prediction by network pharmacology FQD-ingredient and active-compound database The ingredients of each herb contained in FQD were retrieved from online databases, including TCMSP, TCM Database@ Taiwan, Chemistry Database, and TCMID. There were 199, 189, and 33 ingredients for Radix Ginseng, Rhizoma Ligustici Chuanxiong, and Rhizoma Coptidis, respectively. After redundancies had been deleted, a total of 401 ingredients were used to construct the FQD-ingredient database ([99]Table S1). Active compounds were ingredients that were able to be absorbed into the bloodstream and exert some function by protein regulation during disease treatment. Active compounds were screened based on ADME parameters, with 21, 9, and 7 compounds identified for each herb, respectively. Further, compounds known to be effective in treating central nervous system disease were added. The final total of 36 active FQD compounds is shown in [100]Table 1. Table 1. Active compounds in FQD OB (%) DL BBBp Source Target number __________________________________________________________________ R-canadine 55.36687 0.77465 0.57189 Screened 56 Alexandrine_qt[101]^a 36.913906 0.75268 0.87952 Screened 28 Aposcopolamine 66.646907 0.21999 0.40163 Screened 8 Arachidonate 45.57325 0.20491 0.57974 Screened 56 Berberine 36.86125 0.77665 0.56718 Screened 52 Berberrubine 35.73551 0.7269 0.17408 Screened 40 Berlambine 36.6809 0.81596 0.17347 Screened 47 β-sitosterol 36.913906 0.75123 0.98588 Screened 72 Celabenzine 101.8826 0.48772 0.05084 Screened 11 Deoxyharringtonine 39.27444 0.8116 −0.24531 Screened 2 Diop 43.593325 0.39247 0.26014 Screened 8 Epiberberine 43.09233 0.7761 0.39823 Screened 35 Ferulic acid 39.56 0.06 −0.03 Literature[102]^77–[103]^79 56 Frutinone A 65.903731 0.34184 0.4604 Screened 27 Fumarine 59.262505 0.82694 −0.13192 Screened 46 Ginsenoside Rb1 6.24 0.04 −5.6 Literature[104]^80–[105]^84 7 Ginsenoside Re 4.27 0.12 −4.39 Literature[106]^85–[107]^87 4 Ginsenoside Rg1 10.04 0.28 −3.5 Literature[108]^88–[109]^90 5 Ginsenoside Rg5_qt 39.563071 0.78506 0.21243 Screened 0 Ginsenoside Rh4_qt 31.112148 0.77829 −0.18125 Screened 5 Girinimbine 61.215303 0.31484 1.22221 Screened 23 Gomisin B 31.990424 0.82858 0.17993 Screened 15 Inermin 65.830931 0.53754 0.36121 Screened 35 Ligustilide 23.5 0.07 1.2 Literature[110]^91–[111]^93 6 Ligustrazine 20.01 0.03 1.05 Literature[112]^94–[113]^96 6 Magnograndiolide 63.70888 0.18833 −0.24386 Screened 7 Malkangunin 57.713844 0.62642 −0.17138 Screened 5 Mandenol 41.9962 0.19321 1.13821 Screened 8 Myricanone 40.597566 0.51262 −0.07932 Screened 25 Palmatine 64.60111 0.64524 0.37072 Screened 45 Panaxadiol 33.087961 0.79404 0.22661 Screened 3 Perlolyrine 65.947753 0.2747 0.1464 Screened 12 Senkyunone 47.663947 0.24435 0.5013 Screened 6 Stigmasterol 43.829852 0.75665 1.00045 Screened 61 Suchilactone 57.518824 0.55573 0.27669 Screened 32 Wallichilide 42.310679 0.70639 0.73257 Screened 7 [114]Open in a new tab Note: ^a Compounds with the suffix “_qt” refer to the deglycosylated form of their prototypes. Abbreviations: BBBp, blood–brain barrier permeability; DL, drug likeness; FQD, Fuzheng Quxie Decoction; OB, oral bioavailability. Potential targets of FQD Target fishing is the crucial step in mechanism prediction. Each screened active compound was able to interact with specific proteins based on its structure, with the possibility of a protective effect for some pathological process. By uploading the structure of each active compound to the online prediction systems/websites, potential targets/proteins for all the active compounds were acquired with one exception. Predicted targets for each compound are shown in [115]Table 1. Ultimately, 236 predicted targets (putative targets) were obtained after deletion of duplicates (shown in [116]Table S2). Some of the putative targets are related to pathological processes of AD such as acetylcholinesterase, caspase 3, LDLR, brain-derived neurotrophic factor (BDNF), and iNOS. To interpret function, further analyses were performed. AD-related genes AD-related genes were retrieved from Online Mendelian Inheritance in Man, MalaCards, and KEGG to profile proteins that might be involved in the pathogenesis of AD. Of these, 127, 293, and 171 genes were obtained from the three databases, respectively. All these were merged to construct an AD-related gene database comprised of 456 genes in total ([117]Table S3). Determination of DRTs and PRTs by PPI-network construction and analysis After intersection of the putative FQD targets and AD-related genes, 31 DRTs were defined and considered the key proteins that FQD can regulate in AD ([118]Figure 2 and [119]Table 2). With 31 DRTs as seeds, a PPI network with 1,814 nodes and 34,063 edges was constructed. The BisoGenet plug-in was used to calculate the DC of each node in the network, with an average DC value for all nodes of 37. There were 245 nodes with greater DC value than twice the average, and the genes represented by these nodes can be regarded as PRTs of FQD for treatment of AD ([120]Figure 3 and [121]Table S4). Figure 2. Figure 2 [122]Open in a new tab Intersections of putative FQD targets and AD-related genes; 32 shared genes were obtained, which were defined as DRTs. Abbreviations: AD, Alzheimer’s disease; DRTs, direct regulatory targets; FQD, Fuzheng Quxie Decoction. Table 2. Shared genes of putative FQD targets and AD-related genes (DRTs) Target/protein name UniProt ID Gene __________________________________________________________________ ATP-binding cassette subfamily A member 1 ABCA1_HUMAN ABCA1 Acetylcholinesterase ACES_HUMAN ACHE Apolipoprotein E APOE_HUMAN APOE Cholinesterase CHLE_HUMAN BCHE Brain-derived neurotrophic factor BDNF_HUMAN BDNF Voltage-dependent L-type calcium channel subunit α[1S] CAC1S_HUMAN CACNA1S Calmodulin 1 CALM_HUMAN CALM1 Caspase 3 CASP3_HUMAN CASP3 C-C chemokine receptor type 2 CCR2_HUMAN CCR2 Choline O-acetyltransferase CLAT_HUMAN CHAT Neuronal acetylcholine receptor protein, α[7] chain ACHA7_HUMAN CHRNA7 Cytochrome C, somatic CYC_HUMAN CYCS Cytochrome P450 2D6 CP2D6_HUMAN CYP2D6 24-Dehydrocholesterol reductase DHC24_HUMAN DHCR24 Dual-specificity tyrosine-phosphorylation-regulated kinase 1A DYR1A_HUMAN DYRK1A Glycogen synthase kinase 3β GSK3B_HUMAN GSK3B Heme oxygenase 1 HMOX1_HUMAN HMOX1 5-Hydroxytryptamine 2A receptor 5HT2A_HUMAN HTR2A Low-density-lipoprotein receptor LDLR_HUMAN LDLR Low-density-lipoprotein receptor-related protein 8 LRP8_HUMAN LRP8 Amine oxidase (flavin-containing) A AOFB_HUMAN MAOB Microtubule-associated protein 2 MTAP2_HUMAN MAP2 Mitogen-activated protein kinase 1 MK01_HUMAN MAPK1 Mitogen-activated protein kinase 3 MK03_HUMAN MAPK3 Microtubule-associated protein tau TAU_HUMAN MAPT Nitric oxide synthase, inducible NOS2_HUMAN NOS2 Nitric oxide synthase, endothelial NOS3_HUMAN NOS3 Urokinase-type plasminogen activator UROK_HUMAN PLAU Prostaglandin G/H synthase 2 PGH2_HUMAN PTGS2 Sodium-dependent serotonin transporter SC6A4_HUMAN SLC6A4 Tumor necrosis factor TNFA_HUMAN TNF Very-low-density-lipoprotein receptor VLDLR_HUMAN VLDLR [123]Open in a new tab Abbreviations: FQD, Fuzheng Quxie Decoction; DRTs, direct regulatory targets. Figure 3. [124]Figure 3 [125]Open in a new tab PPI network of DRTs and PRTs. Notes: (A) PPI network with 32 DRTs as seeds. Pink nodes represent the 32 DRTs, blue nodes represent the proteins that have been reported to have interactions with the DRTs, and yellow nodes represent 245 PRTs with a DC value greater than twice the average. (B) PPI network of the 245 chosen PRTs, with nine DRTs (pink nodes). Abbreviations: DC, degree centrality; DRTs, direct regulatory targets; PPI, protein–protein interaction; PRTs, possible regulatory targets. GO and pathway-enrichment analysis for PRTs GO enrichment analysis annotated all the PRTs for three aspects: biological process, cellular component, and molecular function. The most significantly enriched terms for these three categories are shown in [126]Figure 4. Most of the 245 PRTs were involved in cellular processes, metabolic processes, and signaling transduction, through binding of other proteins. Figure 4. [127]Figure 4 [128]Open in a new tab GO enrichment of 245 PRTs. Most significantly enriched terms are displayed. Abbreviations: GO, gene ontology; PRTs, possible regulatory targets. Pathway-enrichment analysis can assist in identification of major signaling pathways in which genes participate. Through pathway enrichment, 245 PRTs were enriched in the VEGF signaling network. Pathways of activator protein 1 (AP1) transcription factor, small mothers against decapentaplegic homolog 2/3 (SMAD2/3), transforming growth factor β (TGFβ), and integrin-linked kinase were also significantly enriched, as shown in [129]Figure 5. VEGF promotes the proliferation and migration of the endothelium, playing a vital role in vascular permeability, especially the permeability of the BBB. Vascular impairment and BBB compromise have been widely acknowledged as important contributors to AD pathology. Regulating the expression of VEGF and involved pathways may ameliorate cerebrovascular function and thus protect neural function. Based on this network pharmacological analysis, it is hypothesized that FQD treatment of AD is by regulation of the VEGF-signaling pathway. As HIF1α is a transcription factor that regulates the expression of VEGF, a role for HIF1α was investigated. Figure 5. [130]Figure 5 [131]Open in a new tab The top-ten enriched pathways of the 245 PRTs analyzed by FunRich 3.0. Abbreviations: PRTs, possible regulatory targets; VEGF, vascular endothelial growth factor; VEGFR, VEGF-receptor. Animal experimentation Effect of FQD on spatial learning and memory ability In the MWM test, when compared to the control group, the ratio of time spent in the platform quadrant to total swimming (time ratio) decreased significantly in the model group, which was also observed for distance ratio (P<0.01). When compared to the model group, both the time ratio and distance ratio increased in drug-administered groups, which was especially significant in the FQD high-dose group (P<0.01 and P<0.05, respectively; [132]Figure 6). Figure 6. [133]Figure 6 [134]Open in a new tab Effect of Fuzheng Quxie Decoction (FQD) on spatial learning and memory ability of SAMP8 mice. Notes: (A) Ratio of time in platform quadrant to total swimming time; (B) ratio of swimming distance in platform quadrant to total swimming distance. (^▲▲P<0.01 vs control group; *P<0.05, **P<0.01 vs model group; n=10). Effect of FQD on memory-retention ability In the passive-avoidance test, when compared to the control group, avoidance latency was shortened (P<0.01), while error times decreased significantly (P<0.05) in the model group. Compared with the model group, avoidance latency was shortened significantly in the memantine and FQD low-dose groups (P<0.01 and P<0.05, respectively). For all drug-administration groups, error times decreased significantly (P<0.05; [135]Figure 7). Figure 7. [136]Figure 7 [137]Open in a new tab Effect of Fuzheng Quxie Decoction (FQD) on memory-retention ability of SAMP8 mice. Notes: (A) Avoidance latency that mice spent to come into dark compartment for the first time; (B) error times in the dark compartment. (^▲P<0.05, ^▲▲P<0.01 vs control group; *P<0.05, **P<0.01 vs model group; n=10). Effect of FQD on cerebrovascular cortex H&E staining showed that in the cortices of the control group, structures of the cerebrovascular endothelial cells (CECs) were clear, microvascular walls were smooth and intact, and surrounding spaces able to be observed, while in the cortices of the model group, the shape and arrangement of CECs and vascular walls were deformed, and surrounding spaces increased. Moreover, bleeding was observed due to broken vascular walls (or the BBB). The number of capillaries decreased. After administration of memantine and FQD, there were significant improvements in the shape and arrangement of CECs, with vascular walls more intact and smooth, surrounding spaces decreased, and capillary numbers higher. The results are shown in [138]Figure 8. Figure 8. [139]Figure 8 [140]Open in a new tab Effect of Fuzheng Quxie Decoction (FQD) on cerebrovascular morphology in the cortex of SAMP8 mice. Note: H&E staining, magnification 100× and 400×. Effect of FQD on neuronal morphology in hippocampus Nissl staining showed that neurons in the hippocampal CA1 region of the control group were neat, the cytoplasm dark blue, and nuclei light blue, with clear Nissl bodies within them. However, compared to the control group, neurons were disordered, the number of neurons significantly decreased, and the boundaries of the Nissl bodies were blurred, with some nuclei or Nissl bodies pyknotic or even dissolved. When compared to the model group, neurons were arranged better, numbers of neurons and Nissl bodies increased, and nuclear pyknosis partially alleviated after drug administration. The results are shown in [141]Figure 9. Figure 9. [142]Figure 9 [143]Open in a new tab Effect of Fuzheng Quxie Decoction (FQD) on neuronal morphology in hippocampus CA1 region of SAMP8 mice. Note: Nissl staining, magnification 200×. Effect of FQD on VEGF expression With immunohistochemistry and WB, VEGF expression was assessed. On immunohistochemistry, the number of VEGF-positive cells decreased significantly (P<0.01) in the model group when compared to the control group. Compared to the model group, the number of VEGF-positive cells increased significantly in the memantine and FQD high-dose group (P<0.05 and P<0.01, respectively). However, on WB, there was no significant difference in VEGF expression in these groups. Results are shown in [144]Figure 10. Figure 10. [145]Figure 10 [146]Open in a new tab Effect of Fuzheng Quxie Decoction (FQD) on VEGF expression in brains of SAMP8 mice. Notes: (A) VEGF-positive cells in cortex by immunohistochemistry (magnification 400×); (B) quantification of VEGF-positive cells in cortex; (C) VEGF bands detected by Western blot analysis; (D) quantification of VEGF expression, normalized as ratio to β-actin. ^▲▲P<0.01 vs control group; *P<0.05, **P<0.01 vs model group (B and D). Effect of FQD on HIF1α expression On WB, HIF1α expression in the model group increased significantly (P<0.05) when compared to the control group. When compared to the model group, HIF1α expression decreased in the memantine and FQD low-dose groups (P<0.05), and the decrease was especially significant in FQD high-dose group (P<0.01). The results are shown in [147]Figure 11. Figure 11. Figure 11 [148]Open in a new tab Effect of Fuzheng Quxie Decoction (FQD) on hypoxia-induced factor 1α (HIF1α) expression in brains of SAMP8 mice. Notes: (A) HIF1α bands detected by Western blot analysis; (B) shows the quantification of HIF1α expression normalized as the ratio to β-actin. (^▲P<0.05 vs control group; *P<0.05, **P<0.01 vs model group; n=5). Discussion Predictions for the mechanisms of action of Chinese herbal formulas in disease are often assessed through network pharmacological approaches.[149]^50 Herein, the composition of FQD was collected and its predicted active compounds were assessed by ADME parameters. The targets of FQD were identified as acetylcholinesterase, apolipoprotein, low-density lipoprotein receptor (LDLR), caspase 3, brain-derived neurotrophic factor (BDNF), and inducible nitric oxide synthases (iNOS), which are involved in neural transmission, lipid metabolism, cellular apoptosis, and inflammation. Each of these is a key biomarker[150]^51 of AD pathogenesis.[151]^51^–[152]^55 It is likely that AD-symptom improvement could be accomplished by regulation of these biomarkers, which would demonstrate the reliability of this predictive network pharmacology. By GO and KEGG enrichment analysis, FQD was found primarily to regulate the VEGF pathway and may be a means by which to treat AD. Regulating the VEGF pathway has been proposed as an approach to enhance cerebral vasculature function and consequently alleviate AD symptoms.[153]^56^,[154]^57 With increasing age, cerebral blood flow is gradually reduced, and this reduction has been proposed as the basic reason for AD onset.[155]^58^,[156]^59 Reduced cerebral blood flow is attributed to vascular aging and reduced angiogenesis, in which VEGF plays a key role.[157]^60 VEGF interacts with VEGF-receptor in the cerebral endothelium and in neurons, not only acting as the main regulator of angiogenesis but also as a contributor to neurogenesis and neuroprotection.[158]^61 In older people, angiogenesis is impaired by reduced numbers of endothelial progenitor cells and/or reduced expression of VEGF.[159]^62^,[160]^63 VEGF expression in AD patients is reduced,[161]^56 with treatments aimed at VEGF beneficial to neuronal survival.[162]^64^–[163]^66 In this study, FQD treatment upregulated VEGF expression, which may facilitate angiogenesis and increase cerebral blood flow, as well as neurogenesis, with each beneficial to cognition. During hypoperfusion, HIF1α is activated to promote VEGF expression, but if HIF1α is overactivated, large amounts of VEGF will increase the permeability of the BBB. Herein, FQD downregulated HIF1α expression, which would maintain HIF1α and VEGF at stable levels. Besides Aβ and tau hyperphosphorylation in AD, cerebrovascular pathology has been proposed as an important contributor to the onset of AD,[164]^67 in that there is vascular luminal stenosis, decreased numbers of capillaries, thickened basement membranes, and increased BBB permeability. Several vascular risk factors contribute to AD pathology, such as hypertension, atherosclerosis, diabetes mellitus, smoking, and obesity, which attack the cerebral vasculature in the early stages of AD.[165]^68^–[166]^71 This study suggests that prevention or reversal of vascular pathological changes and development of drugs that improve cerebral vasculature may reduce the development of AD. SAMP8 mice have been used widely to study AD.[167]^47^,[168]^72 These mice have features of vascular pathology,[169]^73^,[170]^74 which permitted an assessment of the effect of FQD on regulatory vascular pathways. The memory capacity of mice was assessed by the MWM and passive-avoidance tests, which assess long-term memory and short-term memory, respectively.[171]^75^,[172]^76 Based on these results, FQD had a significant effect on short- and long-term memory in mice. Consistently with these findings, H&E and Nissl staining showed vasculoprotective and neuroprotective effects for FQD. The potential mechanism may be by regulation of the VEGF pathway. Future studies confirming these findings should be carried out in future. Conclusion Numerous compounds are found within each herb. As such, targets for Chinese herbal formulas are difficult to confirm. This study demonstrated the reliability and possibility of network pharmacology for interpretation of the mechanism of action for disease treatment by herbal formulae. Despite their reliability, the predicted results require subsequent validation in vivo. In this study, we not only demonstrated the neuroprotective effects of FQD but also validated the predicted mechanism of action underlying the effect of the herbal formula. As a preliminary study, only the regulatory effect of FQD on the expression of VEGF and HIF1α was explored. The role of other proteins or mRNAs in the VEGF pathway and the regulatory effect of FQD on them should be studied further in AD. Acknowledgments