Skip to main content

ORIGINAL RESEARCH article

Front. Pharmacol., 09 August 2022
Sec. Ethnopharmacology
This article is part of the Research Topic Functional lipidomics and metabolomics for investigating the pharmacological mechanism of Herbal Medicines and their functional compounds View all 7 articles

Therapeutic mechanism of Curcuma aromatica Salisb. rhizome against coronary heart disease based on integrated network pharmacology, pharmacological evaluation and lipidomics

Chenghao Fei&#x;Chenghao Fei1De Ji&#x;De Ji1Huangjin Tong&#x;Huangjin Tong2Yu LiYu Li1Lianlin SuLianlin Su1Yuwen QinYuwen Qin1Zhenhua Bian,Zhenhua Bian1,3Wei Zhang,Wei Zhang1,4Chunqin Mao
Chunqin Mao1*Lin Li
Lin Li1*Tulin Lu
Tulin Lu1*
  • 1School of Pharmacy, Nanjing University of Chinese Medicine, Nanjing, China
  • 2Department of Pharmacy, Affiliated Hospital of Integrated Traditional Chinese and Western Medicine, Nanjing University of Chinese Medicine, Nanjing, China
  • 3Department of Pharmacy, Wuxi Affiliated Hospital of Nanjing University of Chinese Medicine, Wuxi, China
  • 4School of Pharmacy, Anhui University of Chinese Medicine, Hefei, China

Curcuma aromatica Salisb. rhizome (CASR) has multifunctional characteristics worldwide and a long history of use as a botanical drug with. Currently, it is often used clinically to treat coronary heart disease (CHD) caused by blood stasis syndrome. However, the therapeutic mechanism of CASR in the treatment of CHD remains poorly understood. In study, the main chemical constituents of CASR were analyzed using UPLC-Q-TOF-MS/MS. Then, its potential therapeutic mechanism against CHD was predicted. Subsequently, pharmacological evaluation was performed using CHD rat model. Finally, a lipidomics approach was applied to explore the different lipid metabolites to verify the regulation of CASR on lipid metabolism disorders in CHD. A total of 35 compounds was identified from CASR. Seventeen active components and 51 potential targets related to CHD were screened by network pharmacology, involving 13 key pathways. In vivo experiments showed that CASR could significantly improve myocardial infarction, blood stasis, and blood lipid levels and regulate the PI3K/AKT/mTOR signaling pathway in CHD rats. Lipidomics further showed that CASR could regulate abnormal sphingolipid, glycerophospholipid, and glycerolipid metabolism in CHD rats. The therapeutic mechanism of CASR against CHD was initially elucidated and included the regulation of lipid metabolism. Its effects may be attributed to active ingredients, such as curzerene, isoprocurcumenol, and (+)-curcumenol. This study reveals the characteristics of multi-component and multi-pathway of CASR in the treatment of CHD, which provides a basis for the follow-up development and utilization of CASR.

1 Introduction

Coronary heart disease (CHD) is caused by various inducements based on abnormal lipid metabolism (An et al., 2020; Sun et al., 2021). Liquid chromatography-mass spectrometry (LC-MS) was used to detect endogenous metabolic molecular markers in patients with CHD, suggesting severe abnormal lipid metabolism (Cai et al., 2019). Currently, CHD is mainly treated through medication and surgery. The former includes nitrate drugs, antithrombotic drugs, calcium channel blockers, lipid-lowering drugs, and others, while the latter mainly includes percutaneous coronary intervention and coronary artery bypass grafting (Sipahi et al., 2014; Rosendorff et al., 2015).

Curcuma aromatica Salisb. rhizome (CASR) is obtained from plants of the genus Curcuma of the Zingiberaceae family. These plants are mainly found in the tropics and subtropics of South and Southeast Asia and largely cultivated in Bengal, China, and Sri Lanka (Sun et al., 2017). Although over 100 species reportedly exist worldwide, Curcuma aromatica Salisb. is primarily cultivated in China. It promotes blood circulation and reduces pain by regulating qi and relieving blood stasis. It is often used clinically in the treatment of CHD and other diseases (Ding et al., 2015; Jing et al., 2017). CASR mainly contains two categories of components: terpenoids (mainly sesquiterpenes) and curcuminoids (Sun et al., 2017). Modern pharmacological studies have shown that botanical drugs that promote blood circulation and remove blood stasis have many pharmacological effects, such as the regulation of lipid metabolism, dilation of coronary arteries, improvement of microcirculation, and regulation of hemodynamics (Pan et al., 2018). Curcuma oil has been shown to significantly regulate blood lipid indices such as total cholesterol (TC) and triglyceride (TG) and inflammatory factors such as interleukin-2 in atherosclerotic rats, which also provides a clinical basis for CASR in the treatment of CHD caused by atherosclerosis (Singh et al., 2015). However, the pharmacodynamics and therapeutic mechanisms of CASR in CHD treatment remain unclear. Therefore, to further understand the therapeutic mechanism of CASR in the treatment of CHD, we tested a hypothesis that CASR regulates lipid metabolism and improves myocardial infarction through its constituents in the plasma, to treat CHD.

Network pharmacology is a new drug design strategy based on the rapid advances in system biology and multidirectional pharmacology. Lipidomics is an important field of metabolomics. They are widely used to study lipids and their metabolic interactions (Li et al., 2014).

In this study, first the UPLC-Q-TOF-MS/MS technology was used to analyze the main chemical components of CASR, and then the potential therapeutic mechanism for CHD was predicted by network pharmacology. The therapeutic effect of CASR on CHD was preliminarily evaluated by pharmacological verification, and the expression of the components of PI3K/AKT/mTOR signaling pathway, which was previously predicted, was further verified using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). Finally, a lipidomics approach was applied to reveal the regulation of CASR on lipid metabolism disorders in CHD. The flow chart of this study is shown in Figure 1.

FIGURE 1
www.frontiersin.org

FIGURE 1. The flowchart of this study.

2 Materials and methods

2.1 Materials and chemicals

The CASR was collected from Zhejiang, China. The samples were identified as rhizomes of Curcuma aromatica Salisb. by Professor Jianwei Chen at the Nanjing University of Chinese Medicine. Vouchers for the sample are available at the Nanjing University of Chinese Medicine.

More details could be found in Supplementary File S1.

2.2 Chemical profiling of Curcuma aromatica Salisb. rhizome

2.2.1 Preparation of sample solution

The CASR slices and extracts were prepared as described in a study (Tong et al., 2021). Finally, CASR extracts containing 1.0 g/ml raw medicine was obtained.

The treatment method for CASR extracts was based on a previously published article (Zhu et al., 2016) with some modifications. After the elution sequence was completed, all eluents were collected and blended adequately. The mixture was centrifuged at 13,000 rpm for 5 min at 4°C, and the supernatant was diluted twice.

2.2.2 Preparation of standard solution

Known amounts of 10 standards (furanodiene, β-elemene, demethoxycurcumin, curzerene, germacrone, curcumenol, isocurcumenol, curdione, curcumin, and furanodienon) were weighed accurately and dissolved in methanol to form a mixed stock solution (approximately 5 μg/ml).

2.2.3 Liquid chromatographic and mass spectrometric conditions

Chromatographic and mass spectrometric conditions described in a previous study (Hao et al., 2019) were applied with some modifications. The optimized gradient elution was as follows: 0–1 min, 5–25% B; 1–3 min, 25–30% B; 3–13 min, 30–55% B; 13–15 min, 55–70% B; 15–25 min, 70–100% B; 25–28 min, 100–5% B; 28–30 min, 5% B. The flow rate was 0.3 ml/min and the sample injection volume was 1 µl. In addition, an electrospray ionization (ESI) source (AB SCIEX, CA) with positive mode was used to acquire mass spectra.

2.2.4 Identification of constituents in Curcuma aromatica Salisb. rhizome

The raw data were obtained using Analyst TF 1.6 software (AB SCIEX, United States). The methods used for the identification and statistical analysis of compounds were similar to the methods used in a previous article (Cai et al., 2019).

2.3 Network pharmacology analysis

2.3.1 Screening of active ingredients

The identified components were screened for absorption, distribution, metabolism, and extraction (ADME) parameters, including oral bioavailability (OB), drug-likeness (DL), and human Caucasian colon adenocarcinoma (Caco-2) permeability, provided by the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP, https://tcmspw.com/). The thresholds for OB, DL, and Caco-2 were set to 30%, 0.1, and −0.4, respectively (Wang et al., 2020). In addition, based on the literature, compounds with pharmacological activities were combined.

2.3.2 Acquisition of potential targets of Curcuma aromatica Salisb. rhizome against coronary heart disease

TCMSP and DrugBank (https://www.drugbank.com/) were used to search potential gene targets of the active compounds. The species was defined as homo sapiens. These target names were converted into official gene symbols using the UniProt database (http://www.uniprot.org/). Duplicates were removed after merging, and the target genes of the active ingredients in CASR were isolated.

The keywords “coronary heart disease” were put into the GeneCards database (http://www.genecards.org) and the Online Mendelian Inheritance in Man database (OMIM) to identify targets derived from the literature and tested using experiments. This process did not set any conditions. After removing duplicates, the related targets of CHD were isolated.

The overlapping targets of the active ingredients in CASR with the related targets of CHD were determined using Venny 2.1 online software (https://bioinfogp.cnb.csic.es/tools/venny/index.html), which were used as potential targets of CASR against CHD.

2.3.3 Gene ontology and kyoto encyclopedia of genes and genomes enrichment analysis

Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of the potential targets of CASR against CHD were performed using the Metascape database (http://metascape.org/). Official gene symbols were used and homo sapiens was selected. GO enrichment terms included biological processes (BP), molecular functions (MF), and cellular components (CC). A threshold of p < 0.01 was set, and the targets were sorted according to the gene number involved. OmicStudio (https://www.omicstudio.cn/tool) was used to screen and plot top-ranking biological processes and pathways.

2.3.4 Construction of the active ingredient-target-pathway network

The active ingredients, potential targets of CASR against CHD, and the KEGG pathways were imported to the Cytoscape (Version 3.8.2) software for the interaction network construction. Additionally, the CytoNCA plug-in was used to calculate topological parameters including “Degree,” “Closeness Centrality,” and “Betweenness Centrality”.

2.4 Pharmacological evaluation of Curcuma aromatica Salisb. rhizome against coronary heart disease in vivo

2.4.1 Animals and grouping

Male specific pathogen-free Sprague-Dawley rats (180–200 g) were obtained from the Experimental Animal Center, Hangzhou Medical College, Zhejiang, China (license number: SCXK (ZHE) 2019-0002). The animals were kept in a temperature (25 ± 2°C), humidity (60% ± 5%), and daylight (periodic 12 h light/dark cycle)-controlled room but with free access to standard diet and water. Following 1 week of acclimatization, the rats were divided randomly into a normal control group (NC) (n = 8), model group (M) (n = 15), positive drug (atorvastatin calcium tablets) control group (P) (n = 8), CASR high-dose treatment group (CASR-H) (n = 8), and CASR low-dose treatment group (CASR-L) (n = 8).

2.4.2 Experimental protocol and drug administration

The CHD model group of rats was fed an HFD combined with vitamin D3 and isoproterenol hydrochloride. The modeling method followed that of a previous study (Sun et al., 2019) with some modifications. The experimental period was 17 weeks. The NC group was fed a normal diet, and the other groups were fed an HFD. On the first day of modeling, except for the NC group, the other groups were orally administered with vitamin D3 (600,000 IU/kg). Additional vitamin D3 (100,000 IU/kg) was administered orally at the 2nd, 4th, 6th, and 8th week, and the same volume of normal saline was administered to the rats in the NC group. For the last 2 days of the 17th week, 30 min after oral administration, M, P, and CASR group of rats were subcutaneously injected (on the back) with isoproterenol hydrochloride (85 mg/kg), once each day. The NC group received the same volume of normal saline. The experimental protocol was in accordance with the Regulations of Experimental Animal Administration issued by the State Committee of Science and Technology of the People’s Republic of China.

From week 15 of modeling, the drug administration scheme was as follows: CASR-L group was orally administered with CASR extracts at a dose of 0.63 g/kg/d; CASR-H group was orally administered with CASR extracts at a dose of 0.95 g/kg/d; P group was orally administered with atorvastatin calcium tablets at a dose of 2.1 mg/kg/d; other groups were orally administered with the same amount of normal saline. The CASR dosage was converted equally according to the upper and lower limits stipulated in the Chinese Pharmacopoeia (2020 edition) and referred to previous article (Xie et al., 2020).

During the experimental period, six rats died. Among the 41 rats that survived, 8, 13, 7, 7, and 6 rats belonged to the NC, M, P, CASR-L, and CASR-H groups, respectively. The modelling rate exceeded 84%.

2.4.3 Sample collection and detection

The rats were euthanized 30 min after the final dose. Serum and plasma samples were immediately collected using the abdominal aortic method. Hemorheological parameters, such as whole blood viscosity (WBV) (200 s−1, 30 s−1, 5 s−1, 1 s−1) and erythrocyte aggregation index (EAI), were measured using the SA-9000 Automated Blood Rheology Analyzer (Beijing Succeeder, China). Serum TC, high-density lipoprotein cholesterol (HDL-C), and low-density lipoprotein cholesterol (LDL-C) levels were determined using Roche Cobas 8000 c702 (Roche Diagnostics, United States). Interleukin-6 (IL-6), oxidized low-density lipoprotein (OxLDL), endothelin (ET), and cardiac troponin I (cTn-I) levels in the serum were determined using the ELISA kits. Plasma was stored at −80°C for lipidomic analysis. Simultaneously, the heart was soaked in 4% PFA and embedded in paraffin. Hematoxylin and eosin (H&E) and Masson’s trichrome (MT) were applied to each of sections (5 µm thick). The histopathological samples were then examined under a light microscope (Leica, DFC-259, Germany).

2.4.4 Real-time quantitative PCR

The expression levels of PI3K, AKT, Bcl-2, TP53, NOS3, mTOR and SMS in the myocardial tissue were detected using RT-qPCR (Chen et al., 2022). β-Actin was the reference gene. Relative mRNA expression was calculated using 2−ΔΔCT method. Supplementary Table S1 lists the sequences of all primers.

2.5 Lipidomics analysis

2.5.1 Sample preparation

A plasma sample preparation strategy was employed based on a previous study (Shan et al., 2018). Internal standards (lysoPE (17:1), D5 TG (17:0–17:1–17:0), and PE (17:0), concentrations of approximately 2.5 μg/ml) were used. Finally, the supernatant was used for the analysis. Quality control (QC) samples were prepared simultaneously to ensure that the method and system were suitable.

2.5.2 Liquid chromatographic and mass spectrometric conditions

Lipidomics analysis conditions were based on a previous study (Shan et al., 2018) with some modifications. The optimized gradient elution was as follows: 0–2 min, 15–30% B; 2–2.5 min, 30–48% B; 2.5–11 min, 48–82% B; 11–11.5 min, 82–99% B; 11.5–12 min, 99% B; 12–13 min, 99–15% B; 13–15 min, 15% B. The column oven temperature was maintained at 65°C. The flow rate was 0.6 ml/min, and the sample injection volume was 3 µl (positive mode) and 2 µl (negative mode). An electrospray ionization (ESI) source was used for both the positive and negative modes. The scan range was 215–1800 m/z.

2.5.3 Lipid annotation and data processing

Lipid annotation and data processing methods were used, as described in previous studies (Shan et al., 2018; Yang et al., 2019). MetaboAnalyst 5.0 (https://www.metaboanalyst.ca/) was used for the multivariate analysis. Based on the logarithmic transformation (base 10) and Pareto scaling, statistical models were developed for the total ion chromatogram (TIC) normalized data (SERRF). Partial least-squares discriminant analysis (PLS-DA) was used to develop the classificatory/predictive models. Two factors were combined to select the significantly different lipids: (I) p value <0.05; (II) variable importance in projection (VIP) value >1.0.

2.6 Statistical analysis

Statistical analysis was performed using the SPSS software (version 23.0; IBM, United States). The results were expressed as mean ± standard deviation (SD). Differences between groups were analyzed using unpaired T tests. The significance level was set at p < 0.05.

3 Results

3.1 Constituent analysis in Curcuma aromatica Salisb. rhizome

A typical TIC of the compounds in CASR is shown in Supplementary Figure S1. Through a comparison with the standards and mass spectral fragmentation pathways in the literature, 35 compounds were identified in the CASR. Detailed information on the identified compounds can be found in Supplementary Table S2. These compounds mainly contained two monoterpenes, 31 sesquiterpenes (13 guaiane-type, 2 carabrane-type, 5 elemane-type, 10 germacrane-type, and 1 unique type), and two curcuminoids.

3.2 Prediction of the therapeutic mechanism of Curcuma aromatica Salisb. rhizome against coronary heart disease

3.2.1 Active compounds in Curcuma aromatica Salisb. rhizome

According to ADME screening, 19 candidate compounds were obtained. In addition, 10 compounds with pharmacological activities were combined. For example, curcumin (P23, OB = 4.37%, DL = 0.41, and Caco-2 = 0.35) could reduce lipid levels in the blood and inhibit CHD progression (Li et al., 2015), despite its relatively low OB value. A total of 29 candidate active compounds were screened (Table 1) for further analyses.

TABLE 1
www.frontiersin.org

TABLE 1. Potential active compounds of CASR in the treatment of CHD and their ADME parameters.

3.2.2 Potential targets of Curcuma aromatica Salisb. rhizome against coronary heart disease

Overall, 18 of the 29 candidate active compounds yielded 69 targets, according to TCMSP and DrugBank. Furthermore, 7,716 and 200 targets related to CHD were collected using GeneCards and OMIM, respectively. A total of 7,828 targets were collected after merging and removing the duplicates. The interaction of the targets of the active ingredients in CASR and CHD was considered the potential target of CASR against CHD. A total of 51 shared targets were identified.

3.2.3 GO and KEGG enrichment analysis

The 51 potential target genes of CASR against CHD were enriched and analyzed using the Metascape database. The results of GO and KEGG pathway analyses are shown in Figure 2. There were 33 main GO terms. Enriched biological processes included the circulatory system process, response to drug, response to steroid hormone, cellular response to organic cyclic compounds, and cellular response to nitrogen compounds; molecular functions included G protein-coupled amine receptor activity, nuclear receptor activity, neurotransmitter receptor activity, ubiquitin-like protein ligase binding, and amyloid-beta binding; cellular components included membrane rafts, postsynaptic membranes, dendrites, and endocytic vesicles (Figure 2A). KEGG pathway analysis indicated that CASR exerted regulatory effects through several pathways, including neuroactive ligand-receptor interaction, calcium signaling pathway, estrogen signaling pathway, serotonergic synapse, cholinergic synapse, PI3K-Akt signaling pathway, sphingolipid signaling pathway, platelet activation, bile secretion, retrograde endocannabinoid signaling, and non-alcoholic fatty liver disease (Figure 2B).

FIGURE 2
www.frontiersin.org

FIGURE 2. Analysis results of network pharmacology. (A) GO enrichment analysis. The horizontal axis represents the GO term, and the vertical axis represents the number of genes; (B) KEGG pathway enrichment analysis. The horizontal axis represents the rich factor, and the vertical axis represents the pathways. The bubble size represents the number of targets in the pathway. The bubble color indicates the magnitude of the p-value. The redder the color, the lower the p-value; (C) The active ingredient-target-pathway network. Blue represents the active components, yellow represents targets, and red represents the pathways.

3.2.4 Construction of the interaction network

Based on the results of the analysis, we predicted the specific mechanism of CASR against CHD. A global view of the active ingredient-target-pathway network, including 79 nodes (17 compounds, 51 potential targets, and 11 signaling pathways) and 197 edges, was constructed (Figure 2C). A complex relationship existed among the active ingredients in CASR, potential target genes, and the KEGG pathways related to them. Active ingredients, targets, pathways with degree, closeness centrality, and betweenness centrality ≥ their median values were screened as core elements of CASR against CHD. Finally, 8 core ingredients, 21 core targets, and 5 core pathways were screened. The results are listed in Table 2. The highest ranked compound was β-elemene. Certain target genes were found to be linked to various pathways. For example, CHRM1 participates in the neuroactive ligand-receptor interaction pathway, calcium signaling pathway, cholinergic synapse, and PI3K-Akt signaling pathway, whereas GABRA1 is involved in the neuroactive ligand-receptor interaction and retrograde endocannabinoid signaling pathway. These results suggest that CASR acts on CHD through a multitarget, multipathway, and integrative manner.

TABLE 2
www.frontiersin.org

TABLE 2. Core ingredients, targets and pathways of CASR against CHD.

3.3 Pharmacological verifications

3.3.1 Effect of Curcuma aromatica Salisb. rhizome on cardiac changes in coronary heart disease rats

Histopathological characteristics of CHD rats are shown in Figures 3A,B. In the NC group, the cardiomyocytes were dense and arranged neatly, with nuclei of the same size and shape. A normal space was observed in the myocardium. No inflammatory cell infiltration or myocardial fibrosis was observed. However, in the M group, the arrangement of cardiomyocytes was disordered, and the nuclei were deeply stained. The myocardial space had widened. There was a large amount of degeneration and necrosis of the myocardial fibers with inflammatory cell infiltration. Meanwhile, there was a significant increase in the amount of blue-stained collagen fibers, and extensive fibrosis were observed. In contrast, the degree of myocardial fibrosis, necrosis, and inflammatory cell infiltration in the P group was significantly lower than that of the model group. The therapeutic effect of the CASR groups was close to that of the P group and that of the CASR-H group was more significant.

FIGURE 3
www.frontiersin.org

FIGURE 3. Experimental validation results in vivo. (A) Heart histopathological changes of H&E staining (magnification = 20× and 500×). Nucleus was stained purple-blue and cytoplasm was stained red; (B) Heart histopathological changes of Masson staining (magnification = 20× and 500×). Collagen fibers appear in blue and muscle fiber appear in red; (C) Hemorheological parameters and serum biochemical indicators. All values were presented as mean ± SD (n = 6–8). All p-values were calculated by T test. *p < 0.05, **p < 0.01 compared with NC group, #p < 0.05, ##p < 0.01 compared with M group. NC: normal control group; M: model group; P: positive drug (atorvastatin calcium tablets) control group; CASR-H: Curcuma aromatica Salisb. rhizome high-dose treatment group; CASR-L: Curcuma aromatica Salisb. rhizome low-dose treatment group.

3.3.2 Effect of Curcuma aromatica Salisb. rhizome on hemorheological parameters in coronary heart disease rats

The hemorheological parameters of rats are shown in Figure 3C. The WBV and EAI were significantly higher in the M group (p < 0.01) than in the NC group, indicating that the rats in the M group were in a serious state of blood stasis. However, the hemorheological parameters of the rats in each treatment group decreased significantly (p < 0.01). The results of the P group were close to that of the NC group and that of the CASR-H group was more significant.

3.3.3 Effect of Curcuma aromatica Salisb. rhizome on serum biochemical indicators in coronary heart disease rats

The results for serum lipids and cardiac marker enzymes are summarized in Figure 3C. Serum concentrations of TC, HDL-C, LDL-C, IL-6, OxLDL, ET, and cTn-I were higher in the M group than in the NC group. Except for HDL-C, all other indices in each CASR-administered group showed a downward trend in comparison to those in the M group. Most of the changes were significant, and the curative effect of CASR was comparable to that of the positive drug.

3.3.4 Effect of Curcuma aromatica Salisb. rhizome on the PI3K/AKT/mTOR signaling pathway

To further investigate the mechanism of CASR in the treatment of CHD, the combined results of the KEGG enrichment analysis in network pharmacology prediction and the expression analysis of PI3K/AKT/mTOR signaling pathway components was verified. The expression of the key genes PI3K, AKT, Bcl-2, TP53, NOS3, and mTOR in this pathway was detected using RT-qPCR. The mTOR pathway was related to lipid metabolism. As shown in Supplementary Figure S2, CASR significantly downregulated the expression of PI3K, AKT, Bcl-2, TP53, NOS3, mTOR, and SMS in CHD rats.

3.4 Lipidomics analysis

3.4.1 Data quality assurance

The method was validated as described in a previous article (Yang et al., 2019). The relative standard deviations of the peak areas of the lipid internal standards were less than 15.83%, indicating that the experimental operation and instrument conditions were relatively stable.

3.4.2 Identification of lipid compounds

TICs of the samples from the positive and negative ion modes are shown in Supplementary Figure S3. The lipids were identified using MS-DIAL software and the LipidBlast database.

In the positive ion mode, 422 lipids were identified, mainly containing acylcarnitine (CAR), cholesteryl ester (CE), and diacylglycerol (DG). In the negative ion mode, 154 lipids were found, mainly containing bile acid (BA), ceramide (Cer), and cardiolipin (CL). To summarize, 576 lipids from plasma samples were identified.

3.4.3 Metabolic profile analysis of lipids

The datasets of each group were analyzed using PLS-DA, and each point in the figure represents a sample. As shown in Supplementary Figure S4A and Supplementary Figure S5B, the NC and M groups displayed significant separation, while CASR had a restorative effect, especially in the positive ion mode. The permutation test results showed that the PLS-DA model was not over-fitted (p < 0.05), validating the model (Supplementary Figure S4B; Supplementary Figure S5B). In addition, Supplementary Figure S4C; Supplementary Figure S5C showed obvious separation of samples between the NC and M groups. Supplementary Figure S4D; Supplementary Figure S5D showed same separation between the M and CASR-administered groups.

3.4.4 Comparative analysis of differential lipids

A total of 265 differential lipids with the criteria of p < 0.05 and VIP >1.0 were (NC group vs. M group) screened. Among them, 44 were regulated after CASR administration (28 in the positive ion mode and 16 in the negative ion mode). The 44 differentially expressed lipids are shown in Supplementary Table S3. Heat maps are also used to visually show the differences in lipids among the groups (Figures 4A,B).

FIGURE 4
www.frontiersin.org

FIGURE 4. Summary of differential lipids. (A) Heatmap of differential lipids in positive mode; (B) Heatmap of differential lipids in negative mode. Each column represents a sample and each row represents a differential lipid. Red changes to blue when the intensity becomes smaller, and vice versa. The darker the red, the larger the value; (C) Enrichment statistics plot of significantly regulated lipids (M vs. NC group); (D) Enrichment statistics plot of significantly regulated lipids (CASR vs. M group). Each bubble represented a significantly changed lipid group (p < 0.05), and the bubble size reflected the total number of lipids contained in each lipid group. Red represents increased lipids and blue represents decreased lipids. The purple bubbles have both increased and decreased lipids; (E) The metabolic network profile. The names of the possible metabolic pathways are denoted in the green, blue and orange dotted box. NC: normal control group; M: model group; P: positive drug (atorvastatin calcium tablets) control group; CASR: Curcuma aromatica Salisb. rhizome treatment group.

In addition, enrichment analysis was beneficial for obtaining mechanistic insights into the lipids. As shown in Figures 4C,D, the seven lipid classes in the M group were significantly different from those in the NC group. Except for the downregulation of LPC and FA, the trend for other lipid classes was upward as a whole. After CASR treatment, SM, LPC, FA, and OxPC levels were regulated (p < 0.05). SM 34:2; 2O was the most downregulated lipid in SM, LPC 20:1/0:0 was the most upregulated lipid in LPC, FA 20:5 was the most upregulated lipid in FA, and PC 32:0; 3O|PC 16:0_16:0; 3O was the most downregulated lipid in OxPC.

Further analysis performed on the metabolic network (Figure 4E) to understand the interaction between different metabolites showed that CASR in the CHD treatment mainly regulated the sphingolipid metabolism, glycerophospholipid metabolism, and glycerolipid metabolism.

3.5 Correlation analysis between core active ingredients and pharmacodynamic indicators

Among the eight core compounds listed in Table 2, only five were detected in the plasma of the CASR treatment groups. These active compounds were then correlated with pharmacodynamic indicators, such as hemorheological parameters and biochemical indicators. The in vivo detection and analysis methods for compounds in CASR were based on our previous study (Tong et al., 2021). The correlation network is shown in Figure 5. Curzerene, isoprocurcumenol, and (+)-curcumenol were negatively correlated with most pharmacodynamic indicators, including WBV, TC, and ET. Furanodiene was negatively correlated with the biochemical indicators. Curzerenone was negatively correlated with hemorheological parameters and positively correlated with biochemical indicators. β-elemene, demethoxycurcumin, and (+)-limonene were not detected.

FIGURE 5
www.frontiersin.org

FIGURE 5. Correlation analysis results between core active ingredients and pharmacodynamic indicators. Purple triangles represent biochemical factors; blue diamonds represent hemorheological parameters; yellow circles represent active compounds; the thickness of the line represents the magnitude of the correlation; solid and dashed line represents positive and negative correlations, respectively.

4 Discussion

In the network pharmacology analysis, GO enrichment analysis revealed that CASR against CHD is involved in a number of biological processes, such as circulatory processes, response to steroids, and phospholipase C activating G protein-coupled receptor signaling pathway, which are directly related to the pathogenesis of CHD (Jayalalitha et al., 2008; Morgan et al., 2020; Nagoor et al., 2020). The KEGG enrichment analysis also suggested that CASR regulates CHD through a variety of pathways. These pathways were partially confirmed by experiments, such as the neuroactive ligand-receptor interaction pathway (Wang et al., 2014), calcium signaling pathway, estrogen signaling pathway (Fortini et al., 2019), PI3K-Akt signaling pathway (Wang et al., 2019), sphingolipid signaling pathway and bile secretion. These results indicated that CASR in the treatment of CHD could protect nerves, repair vascular smooth muscle necrosis, alter lipid metabolism, and reduce inflammation.

In the “active ingredient-target-pathway” network, according to the order of the active compounds involved, there were eight compounds, each of which with a degree greater than the median. The core active ingredients β-elemene, curzerene, demethoxycurcumin, curzerenone, isoprocurcumenol, furanodiene, (+)-curcumenol, and (+)-limonene, interacted on an average with 10 targets. They have been shown to exert protective effects against cardiovascular diseases. For example, furanodiene, demethoxycurcumin, and (+)-curcumenol exhibit anti-inflammatory activities (Makabe et al., 2006; Li et al., 2019; Pintatum et al., 2020). Curzerenone and isoprocurcumenol inhibit platelet aggregation through the regulation of the MAKR and PI3K/AKT signaling pathways (Tong et al., 2021). The β-elemene could improve atherosclerosis by restoring nitric oxide (NO) levels and alleviating oxidative stress, which is associated with CHD (Liu et al., 2017). It is worth noting that although curcuminoids have anti-inflammatory, antioxidant, and hypolipidemic activities, they have poor solubility and low bioavailability, unless the drug delivery mode can be changed (Tagde et al., 2021). Previous studies have shown that furanodiene can induce apoptosis by caspase and p38 MAPK activation, and this process is also accompanied by the inhibition of Bcl-XL, AKT, Bad, and Bax (Aggarwal et al., 2013). Curcumenol significantly reduced LPS-induced production of NO and proinflammatory cytokines (IL-6 and TNF-α) and the expression of proinflammatory proteins iNOS and COX-2, which are associated with AKT-mediated inhibition of the NF-KB and p38 MAPK signaling pathways (Lo et al., 2015). The above mechanisms are the anti-inflammatory mechanisms of furanodiene and (+)-curcumenol respectively.

Each target interacts with an average of two compounds. For example, CHRM1 interacts with (+)-limonene, curzerene, furanodiene, germacrone, isoprocurcumenol, isocurcumenol, β-elemene, and curdione. Research has shown that CHRM1 plays a critical role in regulating atrial repolarization in IK, Ach. Furthermore, CHRM1 and CHRM2 regulate K+ channels through the action of G proteins. At the same time, CHRM1 increases the heart rate and contractility (Tai et al., 2019). Each signaling pathway was regulated by an average of seven targets. The PI3K/AKT pathway is an important intracellular signal transduction pathway that mediates a variety of biological effects, including regulation of vascular endothelial formation, inflammatory response, cardiomyocyte apoptosis, and autophagy, as well as the downstream mTOR signaling pathway to regulate lipid metabolism, which has a substantial impact on the development of cardiovascular diseases. Through the active ingredient-target-pathway network, it was found that the PI3K/AKT signaling pathway was enriched by eight targets such as NOS3, Bcl-2, and TP53 and targeted by 11 compounds including (+)-limonene, curzerene, furanodiene, germacrone, isoprocurcumenol, isocurcumenol, β-elemene, curdione, curcumol, procurcumenol, and demethoxycurcumin. Most of these compounds are the core ingredients listed in Table 2.

In the pharmacological evaluation, the results of the model group were consistent with those of many clinical studies (Ingelsson et al., 2007; Yao et al., 2017; Sun et al., 2021), indicating that the model was successfully created. CASR not only significantly improved the degree of myocardial fibrosis, inflammatory infiltration, and blood viscosity but also reduced the levels of IL-6, ET, TC, and LDL-C and increased the level of HDL-C. IL-6, and ET are related with vascular endothelial inflammation and functional damage (Hernandez et al., 2018). HDL-C can transport metabolic wastes, such as cholesterol in the blood, to the liver, thereby lowering cholesterol levels in the blood and reducing the occurrence of cardiovascular diseases. Arteriography showed that HDL-C levels were significantly negatively correlated with the degree of arterial stenosis (Yang et al., 2011). Therefore, HDL-C is an anti-atherosclerotic plasma lipoprotein cholesterol and a protective factor against CHD (Klancic et al., 2016). LDL-C is a risk factor for CHD, which transports lipids in the blood to vascular endothelial cells. Excessive lipid levels lead to excessive lipid deposition in vascular endothelial cells, aggravating the degree of atherosclerosis and increasing the risk of CHD (Holmes and Ala-Korpela, 2019). Studies have shown that statins can significantly reduce the level of LDL-C in patients and delay the progression of atherosclerosis, thereby inhibiting the progression of CHD. The results showed that the CASR high-dose administration group had similar therapeutic effect to the positive group. In summary, CASR can resist atherosclerosis and alleviate the progression of CHD. However, no differences in the levels of cTn-I, a biochemical marker of myocardial injury, were observed after administration, probably because it was less affected by CASR. The PI3K/AKT/mTOR is an important signal transduction pathway that regulates cardiac function (Wang et al., 2015). It can affect the occurrence and development of CHD by regulating inflammation, vascular endotheliogenesis, and cardiomyocyte apoptosis and by promoting glucose and lipid metabolism. Phosphatidylinositol 3-kinase (PI3K) is a lipid kinase that activates downstream targets. Protein kinase B (AKT) is a serine/threonine kinase that plays a key role in the regulation of cell metabolism, apoptosis, and vascular endotheliogenesis. B-cell lymphoma-2 (Bcl-2) is an oncogene and important regulator involved in the pathway of mitochondrial apoptosis. The tumor protein p53 (TP53) is a transcription factor that promotes apoptosis. The main function of nitric oxide synthase (NOS3) is to participate in the metabolism of arginine and proline and to catalyze the production of NO. Mammalian target of rapamycin (mTOR) is an atypical serine/threonine protein kinase that is mainly involved in many biological processes such as cell growth, proliferation, cell cycle, protein synthesis, and energy metabolism. The results of the RT-qPCR experiment preliminarily verified that CASR could regulate the PI3K/AKT/mTOR signaling pathway in CHD rats at the gene level.

CASR was proved to improve dyslipidemia of CHD model. The present results showed that after treatment with CASR, the lipids recalled included SM, LPC, FA, and OxPC. Among these lipid compounds, the most significant callbacks were observed in SM 34:2; 2O, LPC 20:1/0:0, FA 20:5, and PC 32:0; 3O|PC 16:0_16:0; 3O, which are presumed to be biomarkers for the regulation of lipid metabolism in CHD by CASR. The metabolic network revealed that CASR plays a major role in CHD treatment through three metabolic pathways: sphingolipid metabolism, glycerophospholipid metabolism, and glycerolipid metabolism. Noteworthily, the sphingolipid signaling pathway is a shared pathway between network pharmacology and lipidomics results. Sphingolipids are key components of organelles and cell membranes, which regulate cell growth, differentiation, and senescence (Mendelson et al., 2014). They include ceramide, sphingomyelin, and sphingosine 1-phosphate. Ceramide is the core of sphingolipid metabolism, and sphingolipid metabolism revolves around its synthesis and decomposition (Airola and Hannun, 2013). The level of SM, which is a downstream decomposition product, increased significantly in the CHD model group and was significantly reduced after CASR administration, similar to previous reports (Sun et al., 2021). Furthermore, it was proved that CASR could significantly downregulate the expression level of SMS gene in the myocardium of rats with CHD. Sphingomyelin synthase (SMS) is the rate-limiting enzyme in the sphingomyelin metabolic pathway, according to the KEGG database analysis. By reshaping the structure of lipid rafts, SM plays a significant role in maintaining the function of lipid rafts and causing inflammation (Xue et al., 2019). Glycerophospholipid metabolism is closely related to inflammation and oxidative stress. This metabolic pathway contains a variety of enzymatic lipid metabolites, such as PC, which can be hydrolyzed by phospholipase A2 to form LPC and polyunsaturated fatty acids (PUFAs). PUFAs can be further oxidized to prostaglandins, thromboxanes, prostacyclins and hydroxylated fatty acids by cyclooxygenase, lipoxygenase and cytochrome P450 enzymes. These enzymatic lipid metabolites play important roles in inflammatory responses and lipid oxidation. In this study, it was found that the content of LPC in plasma of rats in the model group was significantly lower than that in the normal group, but it was called back after being treated with CASR (Yin et al., 2013; Dennis and Norris, 2015; Lu et al., 2017). In addition, Triglyceride metabolism has been found to be associated with CHD in genome-wide studies of CHD, and its metabolic enzymes are widely expressed in ventricular and myocardial cells, which may be related to ischemia (Zhang et al., 2020). Lipolysis determines the inflammatory process in the development of atherosclerosis (Hasham and Pillarisetti, 2006). Under the condition of normalizing PPAR-α expression, glycerol metabolism protects the heart from lipid deposition (Rame et al., 2011). In this study, it was found that CASR could improve the abnormal metabolism of triglyceride induced by CHD, which was speculated to be related to the protection of myocardial cells from lipid accumulation. Combined with the results of the pharmacological verification, a comprehensive illustration of the therapeutic mechanism of CASR against CHD is shown in Figure 6.

FIGURE 6
www.frontiersin.org

FIGURE 6. Comprehensive analysis of the therapeutic mechanism of CASR against CHD. The yellow squares represent the key targets of CASR in the treatment of CHD. Pentagrams represent key pathways involved. The arrows represent activation effects, the T-shaped arrows represent inhibition effects, and the segments represent indirect effects.

In addition, through correlation analysis, it was found that curzerene, isoprocurcumenol and (+)-curcumenol reduced blood viscosity, improved blood lipid levels, and reduced vascular endothelial damage. Furanodiene was responsible only for regulating blood lipid levels and vascular endothelial function. Curzerenone was only responsible for reducing blood viscosity. The above results also show that the CASR has multicomponent and multieffect characteristics.

5 Conclusion

In view of the unclear therapeutic mechanism of CASR against CHD, this study clarified that CHD regulates the PI3K/AKT/mTOR signaling pathway. This in turn regulates downstream sphingolipid metabolism and improves abnormal lipid levels and hemorheology. A reduction in vascular endothelial injury, inflammatory reactions, and lipid deposition was observed using UPLC-Q-TOF-MS/MS, network pharmacology, lipidomics, and RT-qPCR, in the treatment of CHD with CASR. In addition, curzerene, isoprocurcumenol, and (+)-curcumenol were identified as the main active compounds and are expected to be potential novel drugs for treating CHD. The results also showed the multicomponent and multipathway characteristics of CASR in the treatment of CHD.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Ethics statement

The animal study was reviewed and approved by The Experimental Animal Ethics Committee of Nanjing University of Chinese Medicine.

Author contributions

CF performed the lipidomics and wrote the original draft; DJ created the rat model and edited the draft; HT carried out the mathematical analysis and validation; YL performed the bioinformatics analysis; LS applied the software and data visualization; YQ performed the histopathological and molecular biology experiments; ZB and WZ collected the samples and edited the draft. CM, TL, and LL conceived of the project and provided guidance. All data were generated in-house, and no paper mill was used. All authors agree to be accountable for all aspects of work ensuring integrity and accuracy.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 81973483, 82173975); the Postgraduate Research & Practice Innovation Program of Jiangsu Province (KYCX21-1755); and the Blue Project of Jiangsu Province.

Acknowledgments

The authors thank Dr. Zhou Peina and Dr. Luo Zichen for the guidance of molecular biology and lipidomics experiments. The authors also thank Prof. Shan Jinjun for providing the metabolomics detection platform and the following postgraduate students in our team: Zhang Jiuba, Xu Zhen, Xiao Xiaoyan, Zhang Qian, Wang Lu, Yu Mengting, Qu Lingyun, Jia Aomeng, Gong Jingwen, Chen Peng, and Chen Ming. The authors thank Wuhan Metware Biotechnology Co., Ltd. for the technical guidance provided.

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.

Publisher’s note

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.

Supplementary material

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

References

Aggarwal, B. B., Yuan, W., Li, S., and Gupta, S. C. (2013). Curcumin-free turmeric exhibits anti-inflammatory and anticancer activities: Identification of novel components of turmeric. Mol. Nutr. Food Res. 57, 1529–1542. doi:10.1002/mnfr.201200838

PubMed Abstract | CrossRef Full Text | Google Scholar

Airola, M. V., and Hannun, Y. A. (2013). “Sphingolipid metabolism and neutral sphingomyelinases,” in Sphingolipids: Basic science and drug development. Editors E. Gulbins, and I. Petrache (Vienna: Springer Vienna), 57–76.

CrossRef Full Text | Google Scholar

An, H., Kim, J., Gwon, M., Gu, H., Kim, H., Leem, J., et al. (2020). Beneficial effects of SREBP decoy oligodeoxynucleotide in an animal model of hyperlipidemia. Int. J. Mol. Sci. 21, 552. doi:10.3390/ijms21020552

CrossRef Full Text | Google Scholar

Cai, H., Xu, Y., Xie, L., Duan, Y., Zhou, J., Liu, J., et al. (2019). Investigation on spectrum-effect correlation between constituents absorbed into blood and bioactivities of baizhu shaoyao san before and after processing on ulcerative colitis rats by UHPLC/Q-TOF-MS/MS coupled with gray correlation analysis. Molecules 24, E940. doi:10.3390/molecules24050940

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Z., Zhou, L., Ge, Y., Chen, J., Du, W., Xiao, L., et al. (2022). Fuzi decoction ameliorates pain and cartilage degeneration of osteoarthritic rats through PI3K-Akt signaling pathway and its clinical retrospective evidence. Phytomedicine. 100, 154071. doi:10.1016/j.phymed.2022.154071

PubMed Abstract | CrossRef Full Text | Google Scholar

Dennis, E. A., and Norris, P. C. (2015). Eicosanoid storm in infection and inflammation. Nat. Rev. Immunol. 15, 511–523. doi:10.1038/nri3859

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, F., Zhang, Q., Ung, C. O. L., Wang, Y., Han, Y., Hu, Y., et al. (2015). An analysis of chemical ingredients network of Chinese herbal formulae for the treatment of coronary heart disease. PLoS One 10, e0116441. doi:10.1371/journal.pone.0116441

PubMed Abstract | CrossRef Full Text | Google Scholar

Fortini, F., Vieceli, D. S. F., Caliceti, C., Lambertini, E., Pannuti, A., Peiffer, D. S., et al. (2019). Estrogen-mediated protection against coronary heart disease: The role of the Notch pathway. J. Steroid Biochem. Mol. Biol. 189, 87–100. doi:10.1016/j.jsbmb.2019.02.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Hao, M., Ji, D., Li, L., Su, L., Zhang, J., Wang, Q., et al. (2019). Metabolic profiling analysis of three processed rhizomes of Curcuma wenyujin YH Chen et C. Ling by ultra-performance liquid chromatography/time-of-flight mass spectrometry. Pharmacogn. Mag. 15, 164. doi:10.4103/pm.pm_399_18

CrossRef Full Text | Google Scholar

Hasham, S. N., and Pillarisetti, S. (2006). Vascular lipases, inflammation and atherosclerosis. Clin. Chim. Acta. 372, 179–183. doi:10.1016/j.cca.2006.04.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Hernandez, M., Wicz, S., Santamaria, M. H., and Corral, R. S. (2018). Curcumin exerts anti-inflammatory and vasoprotective effects through amelioration of NFAT-dependent endothelin-1 production in mice with acute Chagas cardiomyopathy. Mem. Inst. Oswaldo Cruz 113, e180171. doi:10.1590/0074-02760180171

PubMed Abstract | CrossRef Full Text | Google Scholar

Holmes, M. V., and Ala-Korpela, M. (2019). What is ‘LDL cholesterol. Nat. Rev. Cardiol. 16, 197–198. doi:10.1038/s41569-019-0157-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Ingelsson, E., Schaefer, E. J., Contois, J. H., McNamara, J. R., Sullivan, L., Keyes, M. J., et al. (2007). Clinical utility of different lipid measures for prediction of coronary heart disease in men and women. J. Am. Med. Assoc. 298, 776–785. doi:10.1001/jama.298.7.776

CrossRef Full Text | Google Scholar

Jayalalitha, G., Shanthoshini, D. V., and Uthayakumar, R. (2008). Fractal model for blood flow in cardiovascular system. Comput. Biol. Med. 38, 684–693. doi:10.1016/j.compbiomed.2008.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Jing, Y., Jian, L., Li, L., Ning, Z., Xuyan, N., Xiaojuan, H., et al. (2017). Mechanism of herbal pairs with the properties of Qi-tonifying, blood activation, blood-stasis breaking in treating coronary heart disease. J. Traditional Chin. Med. 37, 269–278. doi:10.1016/S0254-6272(17)30054-7

CrossRef Full Text | Google Scholar

Klancic, T., Woodward, L., Hofmann, S. M., and Fisher, E. A. (2016). High density lipoprotein and metabolic disease: Potential benefits of restoring its functional properties. Mol. Metab. 5, 321–327. doi:10.1016/j.molmet.2016.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C., Miao, X., Li, F., Adhikari, B. K., Liu, Y., Sun, J., et al. (2019). Curcuminoids: Implication for inflammation and oxidative stress in cardiovascular diseases. Phytother. Res. 33, 1302–1317. doi:10.1002/ptr.6324

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M., Yang, L., Bai, Y., and Liu, H. (2014). Analytical methods in lipidomics and their applications. Anal. Chem. 86, 161–175. doi:10.1021/ac403554h

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Lu, Y., Sun, Y., and Zhang, Q. (2015). Effect of curcumin on permeability of coronary artery and expression of related proteins in rat coronary atherosclerosis heart disease model. Int. J. Clin. Exp. Pathol. 8, 7247–7253.

PubMed Abstract | Google Scholar

Liu, M., Chen, X., Ma, J., Hassan, W., Wu, H., Ling, J., et al. (2017). β-Elemene attenuates atherosclerosis in apolipoprotein E-deficient mice via restoring NO levels and alleviating oxidative stress. Biomed. Pharmacother. 95, 1789–1798. doi:10.1016/j.biopha.2017.08.092

PubMed Abstract | CrossRef Full Text | Google Scholar

Lo, J. Y., Kamarudin, M. N., Hamdi, O. A., Awang, K., and Kadir, H. A. (2015). Curcumenol isolated from Curcuma zedoaria suppresses Akt-mediated NF-κB activation and p38 MAPK signaling pathway in LPS-stimulated BV-2 microglial cells. Food Funct. 6, 3550–3559. doi:10.1039/c5fo00607d

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, J., Chen, B., Chen, T., Guo, S., Xue, X., Chen, Q., et al. (2017). Comprehensive metabolomics identified lipid peroxidation as a prominent feature in human plasma of patients with coronary heart diseases. Redox Biol. 12, 899–907. doi:10.1016/j.redox.2017.04.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Makabe, H., Maru, N., Kuwabara, A., Kamo, T., and Hirota, M. (2006). Anti-inflammatory sesquiterpenes from Curcuma zedoaria. Nat. Prod. Res. 20, 680–685. doi:10.1080/14786410500462900

PubMed Abstract | CrossRef Full Text | Google Scholar

Morgan, P. K., Fang, L., Lancaster, G. I., and Murphy, A. J. (2020). Hematopoiesis is regulated by cholesterol efflux pathways and lipid rafts: Connections with cardiovascular diseases. J. Lipid Res. 61, 667–675. doi:10.1194/jlr.TR119000267

PubMed Abstract | CrossRef Full Text | Google Scholar

Nagoor, M. M., Azimullah, S., Laham, F., Tariq, S., Goyal, S. N., Adeghate, E., et al. (2020). α-Bisabolol protects against β-adrenergic agonist-induced myocardial infarction in rats by attenuating inflammation, lysosomal dysfunction, NLRP3 inflammasome activation and modulating autophagic flux. Food Funct. 11, 965–976. doi:10.1039/c9fo00530g

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, X., Hao, E., Xie, J., Wei, W., Qin, J., Hou, X., et al. (2018). Molecular mechanism of blood-activating and stasis-removing Chinese medicine in regulating blood stasis syndrome. Chin. J. Exp. Traditional Med. Formulae 24, 227–234. doi:10.13422/j.cnki.syfjx.20182434

CrossRef Full Text | Google Scholar

Pintatum, A., Maneerat, W., Logie, E., Tuenter, E., Sakavitsi, M. E., Pieters, L., et al. (2020). In vitro anti-inflammatory, anti-oxidant, and cytotoxic activities of four Curcuma species and the isolation of compounds from Curcuma aromatica rhizome. Biomolecules 10, E799. doi:10.3390/biom10050799

PubMed Abstract | CrossRef Full Text | Google Scholar

Rame, J. E., Barouch, L. A., Sack, M. N., Lynn, E. G., Abu-Asab, M., Tsokos, M., et al. (2011). Caloric restriction in leptin deficiency does not correct myocardial steatosis: Failure to normalize PPAR{alpha}/PGC1{alpha} and thermogenic glycerolipid/fatty acid cycling. Physiol. Genomics 43, 726–738. doi:10.1152/physiolgenomics.00088.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosendorff, C., Lackland, D. T., Allison, M., Aronow, W. S., Black, H. R., Blumenthal, R. S., et al. (2015). Treatment of hypertension in patients with coronary artery disease: A scientific statement from the American heart association, American College of cardiology, and American society of hypertension. J. Am. Coll. Cardiol. 65, 1998–2038. doi:10.1016/j.jacc.2015.02.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Shan, J., Qian, W., Shen, C., Lin, L., Xie, T., Peng, L., et al. (2018). High-resolution lipidomics reveals dysregulation of lipid metabolism in respiratory syncytial virus pneumonia mice. RSC Adv. 8, 29368–29377. doi:10.1039/c8ra05640d

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh, V., Rana, M., Jain, M., Singh, N., Naqvi, A., Malasoni, R., et al. (2015). Curcuma oil attenuates accelerated atherosclerosis and macrophage foam-cell formation by modulating genes involved in plaque stability, lipid homeostasis and inflammation. Br. J. Nutr. 113, 100–113. doi:10.1017/S0007114514003195

PubMed Abstract | CrossRef Full Text | Google Scholar

Sipahi, I., Akay, M. H., Dagdelen, S., Blitz, A., and Alhan, C. (2014). Coronary artery bypass grafting vs percutaneous coronary intervention and long-term mortality and morbidity in multivessel disease: meta-analysis of randomized clinical trials of the arterial grafting and stenting era. JAMA Intern. Med. 174, 223–230. doi:10.1001/jamainternmed.2013.12844

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, G., Li, X., Wei, J., Zhang, T., Li, B., Chen, M., et al. (2021). Pharmacodynamic substances in Salvia miltiorrhiza for prevention and treatment of hyperlipidemia and coronary heart disease based on lipidomics technology and network pharmacology analysis. Biomed. Pharmacother. 141, 111846. doi:10.1016/j.biopha.2021.111846

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, H., Li, X., Zhang, A., Zhang, K., Yan, G., Han, Y., et al. (2019). Exploring potential biomarkers of coronary heart disease treated by Jing Zhi Guan Xin Pian using high-throughput metabolomics. RSC Adv. 9, 11420–11432. doi:10.1039/c8ra10557j

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, W., Wang, S., Zhao, W., Wu, C., Guo, S., Gao, H., et al. (2017). Chemical constituents and biological research on plants in the genus Curcuma. Crit. Rev. Food Sci. Nutr. 57, 1451–1523. doi:10.1080/10408398.2016.1176554

PubMed Abstract | CrossRef Full Text | Google Scholar

Tagde, P., Tagde, P., Islam, F., Tagde, S., Shah, M., Hussain, Z. D., et al. (2021). The multifaceted role of curcumin in advanced nanocurcumin form in the treatment and management of chronic disorders. Molecules 26, 7109. doi:10.3390/molecules26237109

PubMed Abstract | CrossRef Full Text | Google Scholar

Tai, J., Zou, J., Zhang, X., Wang, Y., Liang, Y., Guo, D., et al. (2019). Using network pharmacology to explore potential treatment mechanism for coronary heart disease using chuanxiong and jiangxiang essential oils in jingzhi guanxin prescriptions. Evid. Based. Complement. Altern. Med. 2019, 7631365. doi:10.1155/2019/7631365

CrossRef Full Text | Google Scholar

Tong, H., Yu, M., Fei, C., Ji, D., Dong, J., Su, L., et al. (2021). Bioactive constituents and the molecular mechanism of Curcumae Rhizoma in the treatment of primary dysmenorrhea based on network pharmacology and molecular docking. Phytomedicine. 86, 153558. doi:10.1016/j.phymed.2021.153558

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, D., Tian, L., Shi, C., Wei, Y., Wang, H., Liu, T., et al. (2020). Network pharmacology-based prediction of the active ingredients and mechanism of Shen Gui capsule for application to coronary heart disease. Comput. Biol. Med. 122, 103825. doi:10.1016/j.compbiomed.2020.103825

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Lu, L., Tan, Y., Jiang, L., Zhao, M., Gao, E., et al. (2019). GPR 30 reduces myocardial infarct area and fibrosis in female ovariectomized mice by activating the PI3K/AKT pathway. Life Sci. 226, 22–32. doi:10.1016/j.lfs.2019.03.049

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Li, C., Liu, Z., Shi, T., Wang, Q., Li, D., et al. (2014). DanQi Pill protects against heart failure through the arachidonic acid metabolism pathway by attenuating different cyclooxygenases and leukotrienes B4. BMC Complement. Altern. Med. 14, 67. doi:10.1186/1472-6882-14-67

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z. G., Wang, Y., Huang, Y., Lu, Q., Zheng, L., Hu, D., et al. (2015). bFGF regulates autophagy and ubiquitinated protein accumulation induced by myocardial ischemia/reperfusion via the activation of the PI3K/Akt/mTOR pathway. Sci. Rep. 5, 9287. doi:10.1038/srep09287

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, H., Su, D., Zhang, J., Ji, D., Mao, J., Hao, M., et al. (2020). Raw and vinegar processed Curcuma wenyujin regulates hepatic fibrosis via bloking TGF-β/Smad signaling pathways and up-regulation of MMP-2/TIMP-1 ratio. J. Ethnopharmacol. 246, 111768. doi:10.1016/j.jep.2019.01.045

PubMed Abstract | CrossRef Full Text | Google Scholar

Xue, J., Yu, Y., Zhang, X., Zhang, C., Zhao, Y., Liu, B., et al. (2019). Sphingomyelin synthase 2 inhibition ameliorates cerebral ischemic reperfusion injury through reducing the recruitment of toll‐like receptor 4 to lipid rafts. J. Am. Heart Assoc. 8, e012885. doi:10.1161/JAHA.119.012885

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, D., Liu, X., and Xiang, M. (2011). The correlation between lipids ratio and degree of coronary artery stenosis. High. Blood Press. Cardiovasc. Prev. 18, 53–56. doi:10.2165/11593480-000000000-00000

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, R., Zhang, Y., Qian, W., Peng, L., Lin, L., Xu, J., et al. (2019). Surfactant lipidomics of alveolar lavage fluid in mice based on ultra-high-performance liquid chromatography coupled to hybrid quadrupole-exactive orbitrap mass spectrometry. Metabolites 9, 80. doi:10.3390/metabo9040080

CrossRef Full Text | Google Scholar

Yao, W., Gao, Y., and Wan, Z. (2017). Serum metabolomics profiling to identify biomarkers for unstable Angina. Biomed. Res. Int. 2017, 7657306. doi:10.1155/2017/7657306

PubMed Abstract | CrossRef Full Text | Google Scholar

Yin, H., Zhou, Y., Zhu, M., Hou, S., Li, Z., Zhong, H., et al. (2013). Role of mitochondria in programmed cell death mediated by arachidonic acid-derived eicosanoids. Mitochondrion 13, 209–224. doi:10.1016/j.mito.2012.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y. Y., Zhao, Z. D., Kong, P. Y., Gao, L., Yu, Y. N., Liu, J., et al. (2020). A comparative pharmacogenomic analysis of three classic TCM prescriptions for coronary heart disease based on molecular network modeling. Acta Pharmacol. Sin. 41, 735–744. doi:10.1038/s41401-019-0352-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, T., Liu, X., Wang, X., Cao, G., Qin, K., Pei, K., et al. (2016). Profiling and analysis of multiple compounds in rhubarb decoction after processing by wine steaming using UHPLC-Q-TOF-MS coupled with multiple statistical strategies. J. Sep. Sci. 39, 3081–3090. doi:10.1002/jssc.201600256

PubMed Abstract | CrossRef Full Text | Google Scholar

Glossary

CHD Coronary heart disease

CASR Curcuma aromatica Salisb. rhizome

ESI electrospray ionization source

TIC total ion chromatogram

TC total cholesterol

TG triglyceride

LDL-C low-density lipoprotein cholesterol

HDL-C high-density lipoprotein cholesterol

IL-6 interleukin-6

OxLDL oxidized low-density lipoprotein

ET endothelin

HFD high-fat feed

ISO isoproterenol hydrochloride

cTn-I cardiac troponin I

MTBE methyl tertiary-butyl ether

PFA Paraformaldehyde

ADME absorption, distribution, metabolism, and extraction

OB oral bioavailability

DL drug-likeness

Caco-2 human caucasian colon adenocarcinoma

BP biological process biological process

MF molecular function molecular function

CC cellular component cellular component

NC normal control group

CASR-H Curcuma aromatica Salisb. rhizome high-dose treatment group

CASR-L Curcuma aromatica Salisb. rhizome low-dose treatment group

WBV whole blood viscosity

EAI erythrocyte aggregation index

H&E hematoxylin and eosin

MT Masson trichrome

GO gene ontology

KEGG Kyoto Encyclopedia of Genes and Genomes

BP biological process biological process

MF molecular function molecular function

CC cellular component cellular component

PI3K phosphoinositide 3-kinase

AKT protein kinase B

Bcl-2 B-cell CLL/lymphoma 2

TP53 tumor protein p53

NOS3 nitric oxide synthase 3

mTOR mammalian target of rapamycin

QC quality control

lysoPE lysophosphatidylethanolamine

PE phosphatidylethanolamine

PLS-DA partial least-squares discriminant analysis

PCA principal component analysis

SD standard deviation

CAR acylcarnitine

CE cholesteryl ester

Cer ceramide ceramide

DG diacylglycerol

HexCer hexosylceramide

LPC lysophosphatidylcholine

NAE N-acyl ethanolamine

PC phosphatidylcholine

PI phosphatidylinositol

CASE campesterol ester

SM sphingomyelin

BA bile acid

Cer ceramide ceramide

CL cardiolipin

FA free fatty acid; bileAcid, cholic acid

OxPC oxidized phosphatidylcholine

SMS sphingomyelin synthase

VIPd variable importance in the projection

Keywords: Curcuma aromatica Salisb. rhizome, coronary heart disease, therapeutic mechanism, network pharmacology, lipidomics, pharmacodynamic substances

Citation: Fei C, Ji D, Tong H, Li Y, Su L, Qin Y, Bian Z, Zhang W, Mao C, Li L and Lu T (2022) Therapeutic mechanism of Curcuma aromatica Salisb. rhizome against coronary heart disease based on integrated network pharmacology, pharmacological evaluation and lipidomics. Front. Pharmacol. 13:950749. doi: 10.3389/fphar.2022.950749

Received: 23 May 2022; Accepted: 14 July 2022;
Published: 09 August 2022.

Edited by:

Guangguo Tan, Fourth Military Medical University, China

Reviewed by:

Suling Zeng, University of California, United States
Qilong Li, Shenyang Agricultural University, China

Copyright © 2022 Fei, Ji, Tong, Li, Su, Qin, Bian, Zhang, Mao, Li and Lu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Chunqin Mao, 300555@njucm.edu.cn; Lin Li, lilin@njucm.edu.cn; Tulin Lu, ltl2021@njucm.edu.cn

These authors have contributed equally to this work and share first authorship

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.