ORIGINAL RESEARCH article

Front. Cardiovasc. Med., 03 December 2020

Sec. Atherosclerosis and Vascular Medicine

Volume 7 - 2020 | https://doi.org/10.3389/fcvm.2020.592201

Systematic Pharmacology and GEO Database Mining Revealed the Therapeutic Mechanism of Xuefu Zhuyu Decoration for Atherosclerosis Cardiovascular Disease

  • Center for Gene Diagnosis, Zhongnan Hospital of Wuhan University, Wuhan, China

Abstract

Background: Xuefu Zhuyu decoration (XFZYD), as a traditional Chinese compound recipe, has been used to treat atherosclerosis cardiovascular disease (ASCVD) for thousands of years in China, but its effective compounds and underlying treatment molecular mechanism remains promiscuous, which severely limits its clinical application.

Methods: The effective components and their targets of XFZYD were predicted and screened based on the Traditional Chinese Medicine System Pharmacology (TCMSP) database. The candidate therapeutic targets of ASCVD were screened by Pharmacogenomics Knowledgebase (PharmGKB) and Comparative Toxicogenomics Database (CTD). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses for target proteins were performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID) database. Differentially expressed genes were identified using the GEO2R online tool. Molecular docking was performed by Schrodinger software. To assess the efficacy of the prediction, human umbilical vein endothelial cells (HUVECs) treated with the effective compound of XFZYD were used as the in vitro model.

Results: A total of 108 effective compounds (including quercetin) and 137 candidate therapeutic targets were identified. Analyzing the relationships among effective compounds, candidate therapeutic targets, and signaling pathways, the therapy mechanisms of XFZYD were mainly reflected in the protection of vascular endothelium, anti-inflammatory, antioxidant stress, etc. Accordingly, we found the effective compound of XFZYD (quercetin) decreased intracellular adhesion molecule 1 (ICAM-1) and vascular cell adhesion molecule 1 (VCAM-1) expressions and pro-inflammatory cytokines in HUVECs treated with lipopolysaccharide (LPS), and reduced the adhesion function of HUVECs with monocytes. The inhibitor of the predicted target protein (PTGS2) could further reduce the expressions of VCAM-1, ICAM-1, and TNF-α induced by LPS, and inhibit the adhesion function of HUVECs with monocytes, while PTGS2 agonists partially counteracted the protective effect of quercetin.

Conclusions: In this study, the effective components and potential therapeutic targets of XFZYD for ASCVD treatment were explored from the perspective of systemic pharmacology. The effective component quercetin was verified to protect endothelial cells by reducing endothelial inflammatory response and impeding the attachment of monocytes against the predicted therapeutic target PTGS2.

Introduction

Atherosclerosis cardiovascular disease (ASCVD) was a systemic disease based on atherosclerosis, becoming a leading killer worldwide due to the high morbidity and mortality (1). Atherosclerosis is a pathological status characterized by fibrogenesis, chronic inflammation, lipid accumulation, and vascular wall immunity disorders (2, 3). As atherosclerotic plaques develop into advanced stages, brittle plaques tend to rupture (4), leading to acute cardiovascular events such as ischemic stroke and myocardial infarction.

At present, clinical medications ideas for ASCVD mainly focus on correcting of atherogenic dyslipidemia and platelet aggregation. The clinical efficacy of statins and aspirin has been well-established (5, 6). Although current atherosclerosis medications can partially relieve the symptoms of ASCVD patients, gastrointestinal reactions, liver and kidney damage, etc. occur constantly, due to individual differences, adverse effects of medicine, or usage without doctor's prescription, which greatly disturb the therapeutic efficacy.

Traditional Chinese Medicine (TCM) plays an important role in the adjuvant treatment of ASCVD (7). Compared with single drug prescriptions, TCM compound recipe has the characteristics of multicomponent, multitarget, and multipathway interaction. Xuefu zhuyu decoration (XFZYD), a famous herbal remedy, has been used to relieve symptoms in patients with ASCVD for thousand years in China with few adverse events (8). XFZYD is composed of 11 herbs: Semen Persicae, Flos Carthami, Rhizoma Chuanxiong, Radix Angelicae Sinensis, Radix Paeoniae Rubra, Radix Rehmanniae, Fructus Aurantii, Radix Bupleuri, Radix Platycodonis, Radix Achyranthis Bidentatae, and Radix Glycyrrhizae, and quercetin was shown to be one of the major components of XFZYD evaluated using ultrahigh-performance liquid chromatography with hybrid ion trap time-of-flight mass spectrometry (UHPLC-ESI-IT-TOF-MS) (9). These herbs have been shown to have antiatherosclerosis properties; for example, Rhizoma Chuanxiong was demonstrated to prevent atherosclerosis as well as ischemia–reperfusion injury in clinical trials (10). Radix Glycyrrhizae possessed important antioxidant activity and protective effect against the human lipoprotein oxidative system (11). Moreover, this formula has been proven reliable and effective for curing ASCVD (8, 12) and its risk factors such as hyperlipidemia (13) and hypertension (14). It has been documented to function as an anti-inflammatory agent by inhibiting the PI3K-AKT-mTOR pathway (15, 16) to increase coronary blood flow, to improve the cardiac microcirculation, to accommodate blood lipids (13), and to prevent platelet aggregation, maintaining vessel growth in physiological or repair range to avoid angiogenesis in the atheromatous plaque (17).

Despite increasing number of researches into the cardioprotective effects of XFZYD, the characterization of its effective compounds and the exact mechanisms underlying its therapeutic action are not fully understood. In the present study, the effective components, potential therapeutic targets, and therapeutic pathways of XFZYD for ASCVD were described, and the therapeutic mechanism of quercetin and therapeutic targets PTGS2 on human umbilical vein endothelial cells (HUVECs) were verified in vitro.

Materials and Methods

Identification of Effective Compounds of XFZYD

The effective components in XFZYD were identified from the Traditional Chinese Medicine System Pharmacology (TCMSP, http://lsp.nwu.edu.cn/browse.php) database. The database provides comprehensive information about ingredients in herbs including chemical structure, oral bioavailability (OB), half-life (HL), drug likeness (DL), drug targets, etc. The pharmacokinetic properties including absorption, distribution, metabolism, and excretion (ADME) are important contributors for bioactivities of drugs. In this study, three ADME-related parameters including OB ≥30%, HL ≥4, and DL ≥0.18 were employed to identify the potential effective compounds in XFZYD. As recommended by TCMSP, the compounds with OB ≥30% and HL ≥4 have good absorption and slow metabolism after oral administration, while the compounds with DL ≥0.18 were chemically suitable for drug development.

Prediction of Compound-Related Targets

The compound-related targets were predicted depending on chemical similarities and pharmacophore models via the TCMSP databases. TCMSP compound data were obtained from databases such as DrugBank, HIT, TTD, Pharmacogenomics Knowledgebase (PharmGKB), etc. All the targets obtained above were standardized as gene names and UniProt IDs by searching from UniprotKB database with “Homo sapiens” species.

Identification of ASCVD-Related Therapeutic Targets

These ASCVD-related therapeutic targets were mined from two databases including PharmGKB and Comparative Toxicogenomics Database (CTD). The key words were “atherosclerosis,” “coronary heart disease,” “angina,” “acute coronary syndrome,” “stroke,” “transient ischemic attacks,” and “peripheral arterial disease.” All the targets in PharmGKB and the top 200 targets in CTD based on inference score were selected, and the obtained targets were standardized as gene names and UniProt IDs by searching from UniprotKB database with “Homo sapiens” species.

Network Construction and Topological Analysis

The compound–target network of XFZYD were constructed by Cytoscape v3.7.1 software, which is a tool for analysis and visualization of the biological network (18). The topological analysis was performed by the Network Analyzer module of Cytoscape software. According to the topology of network, degree centrality (DC), betweenness centrality (BC), and closeness centrality (CC) are the most important parameters for measuring the criticality of a node in the network, as well as the important index for new drug discovery and target prediction.

KEGG Pathway Enrichment Analysis

The KEGG pathway enrichment analysis was carried out using the Database for Annotation, Visualization and Integrated Discovery (DAVID) database, which is an online biological knowledgebase and an analytic tool to extract biological information about gene functional classification, functional annotation, and enriched pathways (19). KEGG pathways with P < 0.05 were considered statistically significant.

GEO Database Validation

Peripheral blood RNA expression profiles of three atherosclerosis patients and three controls were obtained from GSE71226, and the differentially expressed genes (DEGs) in GSE71226 were identified using GEO2R (https://www.ncbi.nlm.nih.gov/geo/geo2r/). P < 0.05 and |logFC| ≥ 1.5 were the screening criteria.

Binding Capacity Between Effective Compounds and Key Targets by Molecular Docking

The X-ray crystal structures of the candidate therapeutic targets were taken from the RCSB PDB database, and all 3D structures of these components were obtained from the PubChem database. Molecular docking was performed by Schrodinger software. The compounds and target proteins were input and pretreated to perform the molecular docking command, and finally, the docking score was obtained. The magnitude of the absolute value of docking score is proportional to the strength of bonding.

Cell Culture and Treatment

HUVECs were supplemented with Dulbecco's modified Eagle medium (DMEM, Gibco) containing 10% fetal bovine serum (Gibco) and 1% penicillin/streptomycin (Gibco) in a CO2 incubator (5% CO2 at 37°C). HUVECs were stimulated with lipopolysaccharide (LPS) (10 μg/ml) (Sigma-Aldrich) in the presence or absence of quercetin (25, 50 μM) (Solarbio), celecoxib (PTGS2-selective inhibitor) (10 μM) (MCE, HY-14398), and rebamipide (PTGS2 agonists) (100 nM) (MCE, HY-B0360) for 24 h.

Reverse Transcriptase Polymerase Chain Reaction Analysis

Total RNA of HUVECs was isolated using E.Z.N.A.® HP Total RNA Kit (Omega, USA) according to the manufacturer's instructions. Nanodrop 2000 (Thermo Scientific, USA) was used to detect the concentration and purity of RNA, and complementary DNA (cDNA) was synthesized according to the manufacturer's instructions of PrimeScriptTM RT reagent Kit (Toyobo, Japan). The messenger RNA (mRNA) expression levels of vascular cell adhesion molecule 1 (VCAM-1) and intracellular adhesion molecule 1 (ICAM-1) were detected using quantitative real-time polymerase chain reaction (RT-qPCR) on the BioRad CFX96 (Bio-Rad, USA) with the Fast SYBR® Green PCR Master Mix-PE (Applied Biosystem, USA). Relative expressions were calculated as 2−ΔCt using glyceraldehyde 3-phosphate dehydrogenase (GAPDH) as a reference gene. Primers are listed in Supplementary Table 4.

Adhesion Assay

About 2.5 × 105/well HUVECs were seeded in 12-well plates overnight and treated with LPS (10 μg/ml) in the presence or absence of quercetin, celecoxib, and rebamipide for 24 h. THP-1 cells (5 × 105 cells/mL) were stained with 5 μM BCECF-AM (Beyotime) in the dark for 30 min and cocultured with HUVECs at 37°C for 60 min. Non-adherent cells were removed and washed with phosphate-buffered saline (PBS) for twice. The attached cells were observed with an inverted fluorescence microscope, and the pictures were captured in three random fields per well at 400 × magnification. The number of adherent cells was quantified by ImageJ software.

Statistical Analysis

All statistical analyses were performed on SPSS version 23.0 (SPSS, USA) and GraphPad Prism 8.0 (GraphPad Software, USA). Statistical significance of the differences was analyzed with in independent t-test or Mann–Whitney U-test. P < 0.05 (two-tailed) was considered statistically significant.

Results

Effective Components and Targets of XFZYD

A total of 129 compounds and 256 target proteins in 11 kinds of herbs were identified in XFZYD through the TCMSP database with the criteria of OB ≥30%, HL ≥4, and DL ≥0.18 (Supplementary Table 1).

ASCVD-Related Targets

ASCVD-related therapeutic targets were retrieved in two databases, including 73 in PharmGKB and 582 in CTD. After removing the repeated targets, a total of 620 ASCVD-related therapeutic targets were identified (Supplementary Table 2), 137 of which were overlapped with targets of XFZYD (Figure 1).

Figure 1

Compound–Target Network of XFZYD for the Treatment of ASCVD

The compound–target network was constructed to elaborate the multiplex interplay between compounds and their related targets of XFZYD at a systematic perspective. The compound–target network consists of 245 nodes (108 active components and 137 candidate therapeutic targets) and 945 edges (Figure 2). The compound quercetin has 90 targets, suggesting that it may be critical in the treatment of atherosclerosis (Table 1). Topological analysis was adopted to determine the core targets and compounds of XFZYD in the treatment of ASCVD with the screening criteria “DC ≥ 3, BC ≥ 0.000348, and CC ≥ 0.3754.” The topological network consists of 86 nodes (65 active components and 21 component targets) and 441 edges (Figure 3). From this result, it can be inferred that these high-degree compounds are likely to be the core pharmacodynamic substances in the XFZYD. Besides, the compound–target network illustrated the multicomponent, multitarget characteristics of XFZYD.

Figure 2

Table 1

Mol IDCompoundsOB (%)DLHLHerbTarget number
MOL000006Luteolin36.160.2515.94Flos Carthami
Radix Platycodonis
35
MOL000098Quercetin46.430.2814.4Radix Bupleuri
Flos Carthami
Radix Achyranthis
Bidentatae
Radix Glycyrrhizae
90
MOL000173Wogonin30.680.2317.75Radix Achyranthis
Bidentatae
24
MOL000239Jaranol50.830.2915.5Radix Glycyrrhizae5
MOL000296Hederagenin36.910.755.35Semen Persicae4
MOL000354Isorhamnetin49.60.3114.34Radix Bupleuri
Radix Glycyrrhizae
14
MOL000358Beta-sitosterol36.910.755.36Radix Paeoniae Rubra
Radix Angelicae Sinensis
Flos Carthami
Radix Achyranthis
Bidentatae
Semen Persicae
Fructus Aurantii
14
MOL000392Formononetin69.670.2117.04Radix Glycyrrhizae16
MOL000417Calycosin47.750.2417.1Radix Glycyrrhizae11
MOL000422Kaempferol41.880.2414.74Radix Bupleuri
Radix Glycyrrhizae
Flos Carthami
Radix Achyranthis
Bidentatae
33
MOL000433FA68.960.7124.81Rhizoma Chuanxiong2
MOL000449Stigmasterol43.830.765.57Radix Bupleuri
Radix Paeoniae Rubra
Radix Angelicae
Sinensis
Radix Rehmanniae
Flos Carthami
Radix Achyranthis
Bidentatae
8
MOL000493Campesterol37.580.714.71Semen Persicae3
MOL000497Licochalcone a40.790.2916.2Radix Glycyrrhizae19
MOL001323Sitosterol alpha143.280.785.64Semen Persicae1
MOL0013282,3-Didehydro GA7063.290.57.62Semen Persicae2
MOL0013292,3-Didehydro GA7788.080.537.6Semen Persicae1
MOL001340GA12084.850.458.4Semen Persicae1
MOL001352GA5464.210.5310.19Semen Persicae2
MOL001355GA6365.540.549.85Semen Persicae1
MOL001358Gibberellin 773.80.59.79Semen Persicae2
MOL001361GA8768.850.578.76Semen Persicae1
MOL0013683-O-p-Coumaroylquinic acid37.630.295.15Semen Persicae3
MOL001454Berberine36.860.786.57Radix Achyranthis
Bidentatae
7
MOL001458Coptisine30.670.869.33Radix Achyranthis
Bidentatae
5
MOL001484Inermine75.180.5411.72Radix Glycyrrhizae5
MOL001494Mandenol420.195.39Rhizoma Chuanxiong2
MOL001645Linoleyl acetate42.10.27.48Radix Bupleuri2
MOL001689Acacetin34.970.2417.25Radix Platycodonis14
MOL001792DFV32.760.1817.89Radix Glycyrrhizae6
MOL001924Paeoniflorin53.870.7913.88Radix Paeoniae Rubra2
MOL0021353,5,6,7-Tetramethoxy−2-(3,4,5-Trimethoxyphenyl) chromone40.60.514.39Rhizoma Chuanxiong13
MOL002140Perlolyrine65.950.2712.62Rhizoma Chuanxiong1
MOL002157Wallichilide42.310.716.85Rhizoma Chuanxiong2
MOL002311Glycyrol90.780.679.85Radix Glycyrrhizae8
MOL002341Hesperetin70.310.2715.78Fructus Aurantii3
MOL002565Medicarpin49.220.348.46Radix Glycyrrhizae9
MOL002695Lignan43.320.6514.88Flos Carthami2
MOL0027126-Hydroxykaempferol62.130.2714.29Flos Carthami5
MOL002714Baicalein33.520.2116.25Radix Achyranthis
Bidentatae
Flos Carthami
Radix Paeoniae Rubra
18
MOL002721Quercetagetin45.010.3113.82Flos Carthami3
MOL002773Beta-carotene37.180.584.36Flos Carthami21
MOL002897Epiberberine43.090.786.1Radix Achyranthis
Bidentatae
4
MOL003656Lupiwighteone51.640.3715.63Radix Glycyrrhizae9
MOL003847Inophyllum E38.810.8515.51Radix Achyranthis
Bidentatae
4
MOL0038967-Methoxy-2- methyl isoflavone42.560.216.89Radix Glycyrrhizae15
MOL004328Naringenin59.290.2116.98Radix Glycyrrhizae
Fructus Aurantii
26
MOL004580cis-Dihydroquercetin66.440.2714.51Radix Platycodonis3
MOL0045983,5,6,7-Tetramethoxy-2- (3,4,5-Trimethoxyphenyl) Chromone31.970.5915.54Radix Bupleuri3
MOL004609Areapillin48.960.4116.52Radix Bupleuri3
MOL004805(2S)-2-[4-Hydroxy-3- (3-methylbut-2-enyl)phenyl]−8,8-Dimethyl-2,3- dihydropyrano[2,3-f]chromen−4-one31.790.7214.82Radix Glycyrrhizae7
MOL004806Euchrenone30.290.5715.89Radix Glycyrrhizae4
MOL004808Glyasperin B65.220.4416.1Radix Glycyrrhizae10
MOL004810Glyasperin F75.840.5415.64Radix Glycyrrhizae10
MOL004814Isotrifoliol31.940.427.91Radix Glycyrrhizae8
MOL004815(E)-1-(2,4-dihydroxyphenyl)−3-(2,2-Dimethylchromen-6-yl) prop-2-en-1-one39.620.3516.16Radix Glycyrrhizae10
MOL004824(2S)-6-(2,4-dihydroxyphenyl)−2-(2-hydroxypropan-2-yl)-4 -methoxy-2,3-dihydrofuro [3,2-g]chromen-7-one60.250.634.31Radix Glycyrrhizae10
MOL004827Semilicoisoflavone B48.780.5517.02Radix Glycyrrhizae8
MOL004828Glepidotin A44.720.3516.09Radix Glycyrrhizae11
MOL004835Glypallichalcone61.60.1917.01Radix Glycyrrhizae13
MOL0048388-(6-Hydroxy-2-benzofuranyl)−2,2-dimethyl-5-chromenol58.440.388.71Radix Glycyrrhizae4
MOL004841Licochalcone B76.760.1917.02Radix Glycyrrhizae11
MOL004848Licochalcone G49.250.3215.75Radix Glycyrrhizae10
MOL004855Licoricone63.580.4716.37Radix Glycyrrhizae6
MOL004856Gancaonin A51.080.416.82Radix Glycyrrhizae8
MOL004857Gancaonin B48.790.4516.49Radix Glycyrrhizae10
MOL0048633-(3,4-Dihydroxyphenyl)-5,7- dihydroxy-8-(3-methylbut-2-enyl) chromone66.370.4115.81Radix Glycyrrhizae9
MOL0048645,7-Dihydroxy-3-(4-methoxyphenyl)−8-(3-methylbut-2-enyl)chromone30.490.4114.99Radix Glycyrrhizae10
MOL0048662-(3,4-Dihydroxyphenyl)-5,7- dihydroxy-6-(3-methylbut-2-enyl) chromone44.150.4116.77Radix Glycyrrhizae6
MOL004882Licocoumarone33.210.369.66Radix Glycyrrhizae5
MOL004883Licoisoflavone41.610.4216.09Radix Glycyrrhizae9
MOL004884Licoisoflavone B38.930.5515.73Radix Glycyrrhizae8
MOL004885Licoisoflavanone52.470.5415.67Radix Glycyrrhizae10
MOL004891Shinpterocarpin80.30.736.5Radix Glycyrrhizae12
MOL004898(E)-3-[3,4-Dihydroxy-5- (3-methylbut-2-enyl)phenyl]−1-(2,4-dihydroxyphenyl)prop−2-en-1-one46.270.3115.24Radix Glycyrrhizae8
MOL004903Liquiritin65.690.7417.96Radix Glycyrrhizae3
MOL004907Glyzaglabrin61.070.3521.2Radix Glycyrrhizae10
MOL004910Glabranin52.90.3116.24Radix Glycyrrhizae5
MOL004912Glabrone52.510.516.09Radix Glycyrrhizae10
MOL0049131,3-Dihydroxy-9-methoxy-6- benzofurano[3,2-c]chromenone48.140.438.87Radix Glycyrrhizae7
MOL0049141,3-Dihydroxy-8,9-dimethoxy-6-benzofurano[3,2-c]chromenone62.90.539.32Radix Glycyrrhizae6
MOL004915Eurycarpin A43.280.3717.1Radix Glycyrrhizae9
MOL004924(–)-Medicocarpin40.990.9513.2Radix Glycyrrhizae2
MOL004935Sigmoidin-B34.880.4114.49Radix Glycyrrhizae4
MOL004941(2R)-7-hydroxy-2-(4-hydroxyphenyl) chroman-4-one71.120.1818.09Radix Glycyrrhizae6
MOL004945(2S)-7-Hydroxy-2-(4-hydroxyphenyl)−8-(3-methylbut-2-enyl)chroman-4-one36.570.3217.95Radix Glycyrrhizae7
MOL004948Isoglycyrol44.70.846.69Radix Glycyrrhizae4
MOL004949Isolicoflavonol45.170.4215.55Radix Glycyrrhizae8
MOL004957HMO38.370.2116.56Radix Glycyrrhizae11
MOL0049591-Methoxyphaseollidin69.980.649.53Radix Glycyrrhizae14
MOL004961Quercetin der.46.450.3316.61Radix Glycyrrhizae9
MOL004988Kanzonol F32.470.899.98Radix Glycyrrhizae2
MOL0049896-Prenylated eriodictyol39.220.4116.52Radix Glycyrrhizae4
MOL0049917-Acetoxy-2-methylisoflavone38.920.2617.49Radix Glycyrrhizae12
MOL0049938-Prenylated eriodictyol53.790.415.7Radix Glycyrrhizae3
MOL005000Gancaonin G60.440.3916.13Radix Glycyrrhizae8
MOL005001Gancaonin H50.10.7816.64Radix Glycyrrhizae5
MOL005003Licoagrocarpin58.810.589.45Radix Glycyrrhizae14
MOL005007Glyasperins M72.670.5915.57Radix Glycyrrhizae12
MOL005008Glycyrrhiza flavonol A41.280.613.71Radix Glycyrrhizae8
MOL005012Licoagroisoflavone57.280.4919.64Radix Glycyrrhizae8
MOL005016Odoratin49.950.316.35Radix Glycyrrhizae10
MOL005017Phaseol78.770.589.64Radix Glycyrrhizae9
MOL005018Xambioona54.850.8714.5Radix Glycyrrhizae3
MOL005828Nobiletin61.670.5216.2Fructus Aurantii17
MOL006992(2R,3R)-4-Methoxyl-distylin59.980.315.08Radix Paeoniae Rubra5
MOL013187Cubebin57.130.6412.4Radix Bupleuri4
MOL013381Marmin38.230.314.68Fructus Aurantii2

The characteristics of effective compounds in Xuefu Zhuyu decoration (XFZYD).

Figure 3

KEGG Enrichment Analysis

KEGG pathway enrichment analysis was performed to elucidate related pathways of the 137 candidate therapeutic targets. The representative top 20 pathways based on the number of enriched genes as well as P-value are shown in Figure 4. These key targets were closely related to the tumor necrosis factor (TNF) signaling pathway, PI3K-Akt signaling pathway, vascular endothelial growth factor (VEGF) signaling pathway, Toll-like receptor signaling pathway, etc., participating in the process of atherosclerotic plaque formation, such as oxidative stress, inflammatory response, angiogenesis, etc.

Figure 4

Targets in the Intersection With GEO Database

With the analysis of GSE71226 microarray data of atherosclerosis, 673 differentially expressed genes were identified, among which 133 were upregulated and 540 were downregulated in the atherosclerosis group (Supplementary Table 3). There were six targets in the intersection with 137 candidate therapeutic targets of XFZYD, which included upregulated PTGS2, MMP9, and BCL2L1 and downregulated JUN, VEGFA, and CXCL2 in the atherosclerosis group (Table 2).

Table 2

GroupGene symbolGene titleP-valuelogFC
Upregulated genesBCL2L1BCL2-like 10.01633962.38
MMP9Matrix metallopeptidase 90.01871662.27
PTGS2Prostaglandin-endoperoxide synthase 20.03749681.75
Downregulated genesCXCL2C-X-C motif chemokine ligand 20.00806−1.59
JUNJun proto-oncogene, AP-1 transcription factor subunit0.0008752−1.99
VEGFAVascular endothelial growth factor A0.0011307−1.69

Differentially expressed genes (DEGs) that overlapped with potential therapeutic targets of Xuefu Zhuyu decoration (XFZYD).

Molecular Docking

To confirm the binding capacity between active compounds and key therapeutic targets, molecular docking using Schrodinger was performed. The target proteins' degrees ranked top 5 in compound–target network topology analysis, and six above DEGs overlapped with candidate therapeutic targets were selected as docking objects. The docking results are listed in Supplementary Table 5. There were eight targets binding to the quercetin in these 10 docking targets, indicating that quercetin might be the main component of XFZYD to exert antiatherosclerosis efficiency.

Quercetin Suppressed LPS-Induced Attachment of HUVECs

The mRNA levels of VCAM-1 and ICAM-1 in HUVECs were significantly increased by exposure to 10 μg/ml LPS (both P < 0.0001), which were decreased by quercetin (25 and 50 μM) (Figures 5A,B). Moreover, quercetin could significantly inhibit LPS-induced attachment of monocytes to HUVECs (Figures 5C,D), suggesting that quercetin may alleviate endothelial dysfunction by reducing adhesion of endothelial cells.

Figure 5

Quercetin Inhibited LPS-Induced Elevation of Pro-inflammatory Cytokines

In order to study the effect of quercetin on pro-inflammatory cytokines expression, HUVECs were stimulated with LPS in the presence or absence of quercetin (25, 50 μM) for 24 h. As shown in Figure 6, LPS treatment significantly increased the mRNA levels of IL-1β (P < 0.001), IL-6 (P < 0.0001), and TNF-α (P < 0.0001) in HUVECs, which were reversed by quercetin (25, 50 μM). These results suggested a strong inhibitory effect of quercetin on pro-inflammatory cytokines.

Figure 6

PTGS2 Might Be a Therapeutic Target of XFZYD

The inhibitor of PTGS2 celecoxib could decrease the expressions of VCAM-1, ICAM-1, and TNF-α in HUVECs induced by LPS, and significantly inhibit LPS-induced attachment of monocytes to HUVECs, while the agonists rebamipide attenuated the inhibition effect of quercetin on adhesion molecules VCAM-1, ICAM-1, and pro-inflammatory cytokine TNF-α as well as the number of monocytes attached to endothelial cells (Figure 7), suggesting that high-expression PTGS2 might be the target of quercetin and inhibition of PTGS2 may improve endothelial dysfunction by reducing adhesion of endothelial cells.

Figure 7

Discussion

ASCVD is a kind of complex and multifactorial disease and remains the leading cause of death worldwide (20). TCM has characteristics of multicomponent and multitarget, which can affect different biological processes to control symptoms and solve the fundamental problems.

In the present study, the effective compounds and candidate therapeutic targets in XFZYD for the treatment of ASCVD were 108 and 137, respectively. Moreover, 56.2% of 137 candidate therapeutic targets of XFZYD could be overlapped by at least two effective compounds, which demonstrated the effective compounds in XFZYD worked against ASCVD through a multitarget synergistic way. In addition, 94.44% of 108 effective compounds acted on at least two candidate therapeutic targets. Besides, quercetin was contained in four herbs (Radix Bupleuri, Flos Carthami, Radix Achyranthis Bidentatae, and Radix Glycyrrhizae), luteolin was contained in two herbs (Flos Carthami and Radix Platycodonis), and kaempferol was contained in four herbs (Radix Bupleuri, Radix Glycyrrhizae, Flos Carthami, and Radix Achyranthis Bidentatae). They acted on 90, 35, and 33 candidate therapeutic targets against ASCVD, respectively, suggesting that XFZYD is a combination of multiple herbs, multiple compounds, and multiple targets in the treatment of ASCVD. Since quercetin has the most predicted targets and was identified as the main component of XFZYD by UHPLC-ESI-IT-TOF-MS (9), it was selected as the tested effective compound.

LPS was considered as an important risk factor contributing to endothelial dysfunction, which was usually known as one of the early hallmarks of atherosclerosis (21, 22). Sequentially, monocyte adhesion to endothelial cells was a crucial step in the early stages of atherosclerosis development, in which ICAM-1 and VCAM-1 were crucially involved. Therefore, preventing endothelial injury has aroused considerable attention as a potential therapeutic target for cardiovascular disease treatment. In this study, quercetin could inhibit LPS-induced attachment of monocytes to HUVECs by downregulating expressions of VCAM-1 and ICAM-1 in vitro as well as suppress the expression of pro-inflammatory cytokines, which demonstrated the ability of quercetin to ameliorate endothelial dysfunction. Coincident with our results, quercetin was reported to remarkably suppress endothelial dysfunction and significantly attenuate atherosclerotic lesion progression in high-fat diet (HFD)-fed ApoE−/− mice (23). Likewise, it was demonstrated that quercetin had an inhibitory action on the expression of pro-inflammatory cytokines in hypercholesterolemic diet (HCD)-fed rats (24). Quercetin also exhibited the therapeutic ability through restraining endothelial dysfunction and vascular inflammation produced by LPS in vivo (25). Furthermore, previous studies have shown the therapeutic application of quercetin in the treatment of other diseases. Nanocapsulated quercetin significantly reduced the incidence of hepatocellular carcinoma in rats and decreased the production of TNF-α and IL-6 in liver induced by diethyl nitrosamine (26). Supplementation with quercetin were proved to improve hyperglycemia, dyslipidemia, and antioxidant status in type 2 diabetes (27). Quercetin also has protective effects on neuronal cell death, Aβ-induced oxidative stress, and memory degradation (28). Taken together, our prediction results were in line with previous reports and confirmed the reliability of the prediction on effective compounds of XFZYD against atherosclerosis.

Besides, KEGG results also suggested that XFZYD might reduce the inflammatory response in atherosclerosis principally by affecting TNF signaling pathway and Toll-like receptor (TLR) signaling pathway. As for TNF signaling pathway, its downstream gene IL-1β activates the nuclear factor-κB (NF-κB) signaling pathway (29), induces the production of various pro-inflammatory cytokines such as TNF-α and IL-6, and positively regulates the further activation of NF-κB, resulting in an inflammatory cascade amplification effect (30). TLRs were the well-defined pattern recognition receptors of immune system, participating in the chronic inflammation and immune response in atherosclerosis (31). TLRs engagement with their ligands stimulated pro-inflammatory cytokine production and foam cell generation, mediating the occurrence and development of coronary atherosclerotic plaque by regulating inflammation and immune response (32). Activation of the TLR signaling pathway lead to the production of multiple pro-inflammatory cytokines (IL-6, TNF-α), accelerating the pathological process of atherosclerosis (33, 34). In this study, the effective compound quercetin possessed a restrain function on pro-inflammatory cytokines produced by LPS-stimulated endothelial cells, which partly explained the therapeutic effect of XFZYD on atherosclerosis by inhibiting inflammatory pathways.

As for the predicted therapeutic targets, PTGS2 was selected for the verification, which was reported to be involved in the early atherosclerotic process, and PTGS2 was highly expressed in atherosclerotic lesions from both human (35, 36) and animals (37, 38). PTGS2 promoted early atherosclerotic lesion formation in low-density lipoprotein receptor deficient (LDLR–/–) mice in vivo (37). Selective inhibitor of PTGS2 celecoxib prevented the development of atherosclerotic lesions in the proximal aortas from ApoE–/– mice by reducing the expressions of ICAM-1 and VCAM-1 (39). Interestingly, Metzner et al. believed that PTGS2 inhibitor had a dual effect on atherosclerosis: promoting the occurrence and development of atherosclerosis in the early stage of atherosclerosis, but playing a protective role against atherosclerosis in the late stage of atherosclerosis (40). In present study, we found that the selective PTGS2 inhibitor celecoxib could reduce the expression of adhesion molecules VCAM-1 and ICAM-1 as well as pro-inflammatory cytokine TNF-α on HUVECs induced by LPS and decrease the amount of monocytes attached to HUVECs, while rebamipide (PTGS2 agonists) weakened the inhibitory effect of quercetin on cell adhesion of HUVECs and pro-inflammatory cytokine, implying that quercetin might play an antiatherosclerosis role by restraining PTGS2 to improve the dysfunction of vascular endothelial cells.

However, there were still some shortcomings in this research. Since the pathological development of atherosclerosis involves complex pathological processes, the mechanism predicted above of XFZYD in treating atherosclerosis still needs to be supplemented by in vivo and in vitro experiments. Here, we have verified only one effective component in vitro, and the other effective components still needed to be confirmed experimentally. In the future, we will continue to explore the mechanisms of XFZYD's other effective components in the treatment of atherosclerotic diseases, in the hope of providing new insights for clinical drug development.

Conclusions

In this study, the effective components and potential therapeutic targets of XFZYD for ASCVD treatment were explored from the perspective of systemic pharmacology, and the effective component quercetin was proved to protect injured endothelial cells and reduce endothelial inflammatory response in vitro. In addition, PTGS2 might be a therapeutic target of quercetin in XFZYD.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author contributions

FZ designed the experiments. FZ, BL, YX, and XZ wrote the manuscript and conducted the pharmacology information analyses. BL, BJ, and YZ visualized the data. BL and CW completed the in vitro experiments. All authors read and approved the final manuscript.

Funding

This work was supported by the Grant of National Natural Science Foundation of China (Grant Number 81871722).

Acknowledgments

This manuscript has been released as a preprint at Research Square (https://www.researchsquare.com/article/rs-21338/v1) (41).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2020.592201/full#supplementary-material

References

  • 1.

    JeeYJungKJLeeSBackJHJeeSHChoSI. Smoking and atherosclerotic cardiovascular disease risk in young men: the Korean Life Course Health Study. BMJ Open. (2019) 9:e024453. 10.1136/bmjopen-2018-024453

  • 2.

    LibbyPBuringJEBadimonLHanssonGKDeanfieldJBittencourtMSet al. Atherosclerosis. Nat Rev Dis Primers. (2019) 5:56. 10.1038/s41572-019-0106-z

  • 3.

    MarchioPGuerra-OjedaSVilaJMAldasoroMVictorVMMauricioMD. Targeting early atherosclerosis: a focus on oxidative stress and inflammation. Oxid Med Cell Longev. (2019) 2019:8563845. 10.1155/2019/8563845

  • 4.

    BadimonLVilahurG. Thrombosis formation on atherosclerotic lesions and plaque rupture. J Intern Med. (2014) 276:61832. 10.1111/joim.12296

  • 5.

    ZhangLHeSLiZGanXLiSChengXet al. Apolipoprotein E polymorphisms contribute to statin response in Chinese ASCVD patients with dyslipidemia. Lipids Health Dis. (2019) 18:129. 10.1186/s12944-019-1069-5

  • 6.

    MannuGSMacartneyALambertJRBettencourt-SilvaJHLawnMLyallHet al. The clinical utility of multiplate analyser measurement in platelet function testing following stroke and transient ischaemic attack. Eur J Haematol. (2015) 94:13844. 10.1111/ejh.12406

  • 7.

    RenYQiaoWFuDHanZLiuWYeWet al. Traditional Chinese medicine protects against cytokine production as the potential immunosuppressive agents in atherosclerosis. J Immunol Res. (2017) 2017:7424307. 10.1155/2017/7424307

  • 8.

    YangXXiongXYangGWangJ. Chinese patent medicine Xuefu Zhuyu capsule for the treatment of unstable angina pectoris: a systematic review of randomized controlled trials. Complement Ther Med. (2014) 22:3919. 10.1016/j.ctim.2014.01.003

  • 9.

    FuCXiaZLiuYLuHZhangZWangYet al. Qualitative analysis of major constituents from Xue Fu Zhu Yu Decoction using ultra high performance liquid chromatography with hybrid ion trap time-of-flight mass spectrometry. J Sep Sci. (2016) 39:345768. 10.1002/jssc.201600083

  • 10.

    GuoMLiuYShiD. Cardiovascular actions and therapeutic potential of tetramethylpyrazine. (Active component isolated from Rhizoma Chuanxiong): roles and mechanisms. Biomed Res Int. (2016) 2016:2430329. 10.1155/2016/2430329

  • 11.

    VisavadiyaNPSoniBDalwadiN. Evaluation of antioxidant and anti-atherogenic properties of Glycyrrhiza glabra root using in vitro models. Int J Food Sci Nutr. (2009) 60(Suppl. 2):13549. 10.1080/09637480902877998

  • 12.

    YangTLiXLuZHanXZhaoM. Effectiveness and safety of Xuefu Zhuyu decoction for treating coronary heart disease angina: a systematic review and meta-analysis. Medicine. (2019) 98:e14708. 10.1097/MD.0000000000014708

  • 13.

    WangSQiuXJ. The efficacy of Xue Fu Zhu Yu prescription for hyperlipidemia: a meta-analysis of randomized controlled trials. Complement Ther Med. (2019) 43:21826. 10.1016/j.ctim.2019.02.008

  • 14.

    WangPXiongXLiS. Efficacy and safety of a traditional Chinese herbal formula Xuefu Zhuyu decoction for hypertension: a systematic review and meta-analysis. Medicine. (2015) 94:e1850. 10.1097/MD.0000000000001850

  • 15.

    XingZXiaZPengWLiJZhangCFuCet al. Xuefu Zhuyu decoction, a traditional Chinese medicine, provides neuroprotection in a rat model of traumatic brain injury via an anti-inflammatory pathway. Sci Rep. (2016) 6:20040. 10.1038/srep20040

  • 16.

    WangJYangXChuFChenJHeQYaoKet al. The effects of Xuefu Zhuyu and Shengmai on the evolution of syndromes and inflammatory markers in patients with unstable angina pectoris after percutaneous coronary intervention: a randomised controlled clinical trial. Evid Based Complement Alternat Med. (2013) 2013:896467. 10.1155/2013/896467

  • 17.

    SongJChenWYWuLYZhengLPLinWGaoDet al. A microarray analysis of angiogenesis modulation effect of Xuefu Zhuyu Decoction on endothelial cells. Chin J Integr Med. (2012) 18:5026. 10.1007/s11655-012-1143-6

  • 18.

    ShannonPMarkielAOzierOBaligaNSWangJTRamageDet al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. (2003) 13:2498504. 10.1101/gr.1239303

  • 19.

    JiaoXShermanBTHuang daWStephensRBaselerMWLaneHCet al. DAVID-WS: a stateful web service to facilitate gene/protein list analysis. Bioinformatics. (2012) 28:18056. 10.1093/bioinformatics/bts251

  • 20.

    HerringtonWLaceyBSherlikerPArmitageJLewingtonS. Epidemiology of atherosclerosis and the potential to reduce the global burden of atherothrombotic disease. Circ Res. (2016) 118:53546. 10.1161/CIRCRESAHA.115.307611

  • 21.

    GimbroneMAJrGarcia-CardenaG. Endothelial cell dysfunction and the pathobiology of atherosclerosis. Circ Res. (2016) 118:62036. 10.1161/CIRCRESAHA.115.306301

  • 22.

    ShenYWardNCHodgsonJMPuddeyIBWangYZhangDet al. Dietary quercetin attenuates oxidant-induced endothelial dysfunction and atherosclerosis in apolipoprotein E knockout mice fed a high-fat diet: a critical role for heme oxygenase-1. Free Radic Biol Med. (2013) 65:90815. 10.1016/j.freeradbiomed.2013.08.185

  • 23.

    LuoMTianRLuN. Quercetin inhibited endothelial dysfunction and atherosclerosis in apolipoprotein E-deficient mice: critical roles for NADPH oxidase and heme oxygenase-1. J Agric Food Chem. (2020) 68:1087583. 10.1021/acs.jafc.0c03907

  • 24.

    BhaskarSSudhakaranPRHelenA. Quercetin attenuates atherosclerotic inflammation and adhesion molecule expression by modulating TLR-NF-kappaB signaling pathway. Cell Immunol. (2016) 310:13140. 10.1016/j.cellimm.2016.08.011

  • 25.

    LuNSuiYTianRPengYY. Inhibitive effects of quercetin on myeloperoxidase-dependent hypochlorous acid formation and vascular endothelial injury. J Agric Food Chem. (2018) 66:493340. 10.1021/acs.jafc.8b01537

  • 26.

    MandalAKGhoshDSarkarSGhoshASwarnakarSDasN. Nanocapsulated quercetin downregulates rat hepatic MMP-13 and controls diethylnitrosamine-induced carcinoma. Nanomedicine. (2014) 9:232337. 10.2217/nnm.14.11

  • 27.

    JeongSMKangMJChoiHNKimJHKimJI. Quercetin ameliorates hyperglycemia and dyslipidemia and improves antioxidant status in type 2 diabetic db/db mice. Nutr Res Pract. (2012) 6:2017. 10.4162/nrp.2012.6.3.201

  • 28.

    LiYLGuoHZhaoYQLiAFRenYQZhangJW. Quercetin protects neuronal cells from oxidative stress and cognitive degradation induced by amyloid beta-peptide treatment. Mol Med Rep. (2017) 16:15737. 10.3892/mmr.2017.6704

  • 29.

    RenHYHuangGLLiuWMZhangWLiuYSuGQet al. IL-1beta induced RXRalpha overexpression through activation of NF-kappaB signaling in gastric carcinoma. Biomed Pharmacother. (2016) 78:32934. 10.1016/j.biopha.2016.01.033

  • 30.

    LuoYZhengSG. Hall of Fame among pro-inflammatory cytokines: interleukin-6 gene and its transcriptional regulation mechanisms. Front Immunol. (2016) 7:604. 10.3389/fimmu.2016.00604

  • 31.

    SeneviratneANMonacoC. Role of inflammatory cells and toll-like receptors in atherosclerosis. Curr Vasc Pharmacol. (2015) 13:14660. 10.2174/15701611113116660160

  • 32.

    LinJKakkarVLuX. Essential roles of toll-like receptors in atherosclerosis. Curr Med Chem. (2016) 23:43154. 10.2174/0929867323666151207111408

  • 33.

    OspeltCGayS. TLRs and chronic inflammation. Int J Biochem Cell Biol. (2010) 42:495505. 10.1016/j.biocel.2009.10.010

  • 34.

    JoostenLAAbdollahi-RoodsazSDinarelloCAO'NeillLNeteaMG. Toll-like receptors and chronic inflammation in rheumatic diseases: new developments. Nat Rev Rheumatol. (2016) 12:34457. 10.1038/nrrheum.2016.61

  • 35.

    SongKLiLSunGWeiY. MicroRNA-381 regulates the occurrence and immune responses of coronary atherosclerosis via cyclooxygenase-2. Exp Ther Med. (2018) 15:455763. 10.3892/etm.2018.5947

  • 36.

    BeltonOByrneDKearneyDLeahyAFitzgeraldDJ. Cyclooxygenase-1 and −2-dependent prostacyclin formation in patients with atherosclerosis. Circulation. (2000) 102:8405. 10.1161/01.CIR.102.8.840

  • 37.

    BurleighMEBabaevVROatesJAHarrisRCGautamSRiendeauDet al. Cyclooxygenase-2 promotes early atherosclerotic lesion formation in LDL receptor-deficient mice. Circulation. (2002) 105:181623. 10.1161/01.CIR.0000014927.74465.7F

  • 38.

    LiLLiJYiJLiuHLeiH. Dose-Effect of irbesartan on cyclooxygenase-2 and matrix metalloproteinase-9 expression in rabbit atherosclerosis. J Cardiovasc Pharmacol. (2018) 71:8294. 10.1097/FJC.0000000000000544

  • 39.

    JacobSLaury-KleintopLLanza-JacobyS. The select cyclooxygenase-2 inhibitor celecoxib reduced the extent of atherosclerosis in apo E-/- mice. J Surg Res. (2008) 146:13542. 10.1016/j.jss.2007.04.040

  • 40.

    MetznerJPoppLMarianCSchmidtRManderscheidCRenneCet al. The effects of COX-2 selective and non-selective NSAIDs on the initiation and progression of atherosclerosis in ApoE-/- mice. J Mol Med. (2007) 85:62333. 10.1007/s00109-007-0162-9

  • 41.

    LiangBXiangYZhangX. Systematic pharmacology and GEO database mining revealed the mechanism of XFZYD therapy for ASCVD. Res Square. (2020). 10.21203/rs.3.rs-21338/v1. [Epub ahead of print].

Summary

Keywords

systematic pharmacology, Xuefu Zhuyu decoration, atherosclerosis cardiovascular disease, therapeutic mechanism, GEO database

Citation

Liang B, Xiang Y, Zhang X, Wang C, Jin B, Zhao Y and Zheng F (2020) Systematic Pharmacology and GEO Database Mining Revealed the Therapeutic Mechanism of Xuefu Zhuyu Decoration for Atherosclerosis Cardiovascular Disease. Front. Cardiovasc. Med. 7:592201. doi: 10.3389/fcvm.2020.592201

Received

06 August 2020

Accepted

05 November 2020

Published

03 December 2020

Volume

7 - 2020

Edited by

Hong Chen, Boston Children's Hospital and Harvard Medical School, United States

Reviewed by

Bo Zhu, Boston Children's Hospital and Harvard Medical School, United States; Yunzhou Dong, Boston Children's Hospital and Harvard Medical School, United States

Updates

Copyright

*Correspondence: Fang Zheng

This article was submitted to Atherosclerosis and Vascular Medicine, a section of the journal Frontiers in Cardiovascular Medicine

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics