Yuming Wang1†
Xue Li1†Min Qi2Xiaokai Li1Fangfang Zhang1Yuyu Wang1Junke Wu1Lexin Shu1
Simiao Fan1Yunfei Li1*
Yubo Li1*- 1School of Chinese Materia, Tianjin University of Traditional Chinese Medicine, Tianjin, China
- 2TIPRHUYA Advancing Innovative Medicines Ltd., Tianjin, China
Introduction: YiYiFuZi powder (YYFZ) is a classical formula in Chinese medicine, which is commonly used clinically for the treatment of Chronic Heart Disease (CHD), but it’s pharmacological effects and mechanism of action are currently unclear.
Methods: An adriamycin-induced CHD model rat was established to evaluate the pharmacological effects of YYFZ on CHD by the results of inflammatory factor level, histopathology and echocardiography. Metabolomic studies were performed on rat plasma using UPLC-Q-TOF/MS to screen biomarkers and enrich metabolic pathways; network pharmacology analysis was also performed to obtain the potential targets and pathways of YYFZ for the treatment of CHD.
Results: The results showed that YYFZ significantly reduced the levels of TNF-α and BNP in the serum of rats, alleviated the disorder of cardiomyocyte arrangement and inflammatory cell infiltration, and improved the cardiac function of rats with CHD. The metabolomic analysis identified a total of 19 metabolites, related to amino acid metabolism, fatty acid metabolism, and other metabolic pathways. Network pharmacology showed that YYFZ acts through PI3K/Akt signaling pathway, MAPK signaling pathway and Ras signaling pathway.
Discussion: YYFZ treatment of CHD modulates blood metabolic pattern and several protein phosphorylation cascades but importance specific changes for therapeutic effect require further studies.
1 Introduction
Chronic heart disease (CHD) is a complex clinical syndrome in which structural or functional abnormalities of the heart cause a range of symptoms of diastolic and systolic dysfunction, such as dyspnea, ankle swelling, and fatigue. CHD is the end-stage of heart disease and a major cause of death from cardiovascular disease (Tsutsui, 2022). The main drugs used to treat heart failure are diuretics, angiotensin-converting enzyme (ACE) inhibitors, angiotensin receptor blockers (ARBs), β-blockers, and salt corticosteroid receptor antagonists (MRAs), but the morbidity and mortality rates of the disease remain high (Roger, 2021). Traditional Chinese medicine (TCM) has been used for thousands of years for the treatment of CHD with good results and multi-target advantages, and many studies are showing the use of TCM for the improvement of CHD, which makes the clinical application of TCM in the treatment of CHD possible (Du et al., 2021; Leung et al., 2021; Meng et al., 2021; Liao et al., 2022; Liu et al., 2022).
In Chinese medicine, CHD belongs to the category of “heart paralysis” and “chest paralysis,” and the main pathogenesis of CHD is a deficiency of positive energy and paralysis by cold and dampness. This formula is from the book “Jin Gui Yao Lve.” Coix seed has multiple pharmacological effects such as anti-tumor, improving body immunity, hypoglycemia, anti-inflammatory and analgesic, and regulating lipid metabolism, and its important pharmacological activities in the treatment of cervical cancer, lung cancer, and gastrointestinal tract tumors have been confirmed in clinical practice (Ni et al., 2021; Zhou et al., 2021; Sui and Xu, 2022; Yang et al., 2022). In a related experimental study, Coix seed was found to improve Th1/Th2 cytokines in mice to restore immune homeostasis after administration (Wang H. et al., 2022). Fuzi is commonly used in clinical practice for the treatment of heart diseases such as heart failure, and its main chemical constituents are alkaloids with cardiac analgesic, anti-myocardial ischemic, anti-arrhythmic, and improving the energy metabolism of cardiomyocytes (Zhang et al., 2017; Yan et al., 2020; Chen et al., 2022; Tai et al., 2022).
Metabolomics and network pharmacology are effective tools to elucidate the potential mechanisms of TCM compounding. In recent years, metabolomics has been widely used in biomedical research, and its target is endogenous metabolites with small molecular weight in the organism, and by detecting the changes of endogenous metabolites after the disturbance of the organism, the differential metabolites associated with the disturbance can be identified, and the metabolic pathways can be elucidated in combination with bioinformatics analysis, and then the biological metabolic mechanism can be analyzed (Li et al., 2020; Cui et al., 2021; Nicholson, 2021; Sun C. et al., 2021; Buergel et al., 2022; Santos-Gallego et al., 2022). Network pharmacology is a joint application of bioinformatics, systems biology and multidirectional pharmacology to study drugs. The “multi-target” and “synergistic mechanism” emphasized by cyber pharmacology are in line with the “holistic concept” and “diagnosis and treatment” emphasized by TCM theory. “Through the analysis of genes, proteins, diseases, drugs and other real data obtained from databases and experiments, the intrinsic relationship between multi-component and multi-target effects of TCM on the body can be explored from a systematic and holistic perspective, which is important for explaining the potential mechanisms of TCM treatment (Li et al., 2019; Yao et al., 2020; Li et al., 2021; Shi et al., 2022; Zhou et al., 2022). The combination of network pharmacology and metabolomics can complement each other to reveal the biological significance and mechanism of drug action more comprehensively.
In this study, we first established a rat model of CHD suitable for testing of pharmacological effects and metabolomic analysis on plasma samples to reveal potential biomarkers and metabolic pathways; secondly, we performed network pharmacological analysis on the in vivo components of YYFZ to predict the targets and pathways of action as well as allowing the in-depth study of the material basis and mechanism of action of YYFZ in the treatment of CHD.
2 Materials and methods
2.1 Instruments and reagents
Adriamycin hydrochloride was purchased from Solarbio (China). Purified water was purchased from Watson’s (China). Isoflurane was purchased from Rwd Life Science Co., Ltd. (China). Saline for injection was purchased from SSY Group Limited. Captopril was purchased from Shanghai Xudong Haipu Pharmaceutical Co., Ltd. (China). Coix Seed was purchased from Hunan Yaoshengtang Traditional Chinese Medicine Technology Co., Ltd. (China). Fuzi was purchased from Hebei Meiwei Pharmaceutical Co., Ltd. (China).
Vevo small animal ultrasound imager was purchased from Visualsonics Co., Ltd. (United States). ALLLEGRATM-64R high-speed frozen centrifuge was purchased from Beckman Co., Ltd. (United States). TecanInfiniteF50 enzyme labeler was purchased from Tecan, Co., Ltd. (Switzerland). UPLC/Q-TOF/MS was purchased from Waters Corporation Co., Ltd. (United States). The rat TNF-α kit and rat BNP kit were purchased from Nanjing Jiancheng Biotechnology Co., Ltd. (China). The pathological sections of heart tissue were processed by Hunan Fenghui Biotechnology Co., Ltd. (China).
2.2 Preparation of YYFZ decoction
Weighed coix seeds and Fuzi, mixed in the ratio of 5:3, placed in a round bottom flask, added 10 times the amount of water and soak for 45 min. The extract was then extracted by reflux extraction method for 1 h. The extract was filtered through 3 layers of gauze. The remaining drug residue was added with 8 times the amount of water and extracted for 45 min. The two filtrates were combined and concentrated to a viscous infusion (concentration of 1 g/mL in terms of raw drug). The samples were stored in a refrigerator at 4°C and set aside.
2.3 Animals
Sixty SPF-grade male Wistar rats (190–220 g, purchased from Beijing Spelford Biotechnology Co., Ltd., license number: SCXK (Beijing) 2019–0010) were selected, and the rats were acclimatized and fed at a temperature of 20°C–26°C, relative humidity of 40%–70%, ventilation of 10–15 times/h, and light of 12 h light/dark for 7 days, and all animals The feeding process was carried out by the operating procedures for clean grade laboratory animals. The study was approved by the Animal Ethics Committee of Tianjin University of Traditional Chinese Medicine (TCM-LAEC2021241), and ethical norms were followed in handling animals during the experiment to minimize animal numbers and suffering.
2.4 Establishment of CHD rat models and grouping
After 1 week of acclimatization feeding, 10 rats were randomly selected as the NS group (control group), and the remaining rats were injected intraperitoneally with 1.25 mg/mL aqueous solution of adriamycin hydrochloride at a dose of 1.25 mg/kg twice a week for 8 weeks, with a cumulative dose of 18 mg/kg; the NS group was injected with the same dose of saline. After the end of modeling, the rats were randomly divided into 5 groups: Model group, Captopril group, YYFZ-H group, YYFZ-M group, YYFZ-L group, 10 rats in each group and treated according to the dosing. YYFZ-H, YYFZ-M, and YYFZ-L groups were treated with 5.25 g/kg, 2.63 g/kg and 1.31 g/kg by gavage. The Captopril group was given 6.75 mg/kg by gavage. The NS and Model groups were given the same volume of pure water as a control and administered once a day for a treatment period of 30 days.
2.5 Observation and detection of indicators
2.5.1 General morphological observation
From the beginning of the modeling, the experimental animals were closely observed and recorded for weight, coat condition, mental status, respiratory system condition and survival.
2.5.2 Evaluation of rat CHD model
At the end of the modeling, all rats fasted for 12 h. Ten rats were randomly selected from the NS group and the model group, anesthetized with isoflurane. Cardiac function was measured using a vevo small animal ultrasound imager, with cardiac ejection fraction (EF) less than 55% as the modeling criterion. At the same time, blood was taken from the inner canthus of the rat’s eye, placed at room temperature for 30 min, and then centrifuged for 10 min at 4°C and 3,500 rpm/min in a freezing centrifuge for the detection of BNP and TNF-α.
2.5.3 Detection of cardiac function indexes
The rats were anesthetized with a VMR small animal respiratory anesthesia machine, skin was prepared on the chest using hair removal cream, coupling agent was applied, and an ultrasound probe was used to perform ultrasound examination of the rat’s heart from the long axis of the parasternal bone. In M ultrasound mode, and the rat’s EF, fractional shortening (FS), left ventricular posterior wall end-systolic thickness (LVPW; s), left ventricular posterior wall end-diastolic thickness (LVPW; d), left ventricular end-systolic internal diameter (LVID; s), left ventricular internal diameter end-diastolic internal diameter (LVID; d), left ventricular end-systolic septal thickness (IVS; s), Left ventricular end-diastolic septal thickness (IVS; d), and the average of 3 cardiac cycles were taken for each sample.
2.6 Observation and detection of indicators
At the end of day 30 dosing, all rats fasted without water for 12 h. After anesthetizing the rats, the cardiac function indexes were first tested using VMR, and then subsequent operations were performed under anesthesia.
2.6.1 Collection of blood samples
Blood was collected from rats through the abdominal aorta to the maximum extent possible, and half of the whole blood samples were placed in heparin anticoagulation tubes and the other half in EP tubes, which were left for 30 min and then pretreated. The plasma samples in the anticoagulation tubes were centrifuged at 4°C, 3,500 rpm for 10 min, and the supernatant was frozen in a −80°C refrigerator for subsequent metabolomic analysis; the blood samples in the EP tubes were centrifuged at 4°C, 3,500 rpm for 10 min, and the supernatant was frozen in a −80°C refrigerator for subsequent serum biochemical analysis.
2.6.2 Heart index
The hearts of rats were removed and weighed after blood sampling from the abdominal aorta, and the organ index was calculated: organ index = organ (g)/body weight (kg).
2.6.3 Detection of serum biochemical indexes
The levels of BNP and TNF-α in rat serum were measured by ELISA kits.
2.6.4 Cardiac histopathology
The rat heart was taken and placed in a 4% tissue fixating solution for 24 h, then dehydrated, paraffin-embedded wax blocks were made, which were successively sliced, dissected, baked, and stained by HE. Morphological changes in rat cardiomyocytes were observed under a light microscope.
2.7 Metabolomics analysis
2.7.1 Sample preparation
Frozen rat plasma samples were removed from the −80°C refrigerator and thawed first in a 4°C refrigerator. 100 μL of plasma was pipetted from the samples, 300 μL of chromatographically pure acetonitrile was added at a volume ratio of 1:3, vortexed and mixed for 5 min, then centrifuged in a refrigerated centrifuge at 13,000 rpm/min for 15 min, and the supernatant was aspirated for UPLC-Q-TOF/MS analysis.
2.7.2 Preparation of QC sample
10 μL of each plasma sample was taken in an EP tube, vortexed for 1 min, and then centrifuged at 4°C at 13,200*g for 15 min, and the supernatant was aspirated for methodological investigation.
2.7.3 Instrument conditions
In this experiment, UPLC-Q-TOF/MS was used to characterize metabolites. UPLC separation was performed on an ACQUITY UPLC BEH C18 column (2.1 mm × 100 mm, 1.7 μm, Waters Co., United States) at 45°C. The mobile phase consists of water (A) and acetonitrile (B) (both containing 0.1% formic acid). The gradient elution procedure is as follows: 0–0.5 min, 1%B; 0.5–2 min, 1%-50%B; 2–9 min, 50%-99%B; 9–10 min, 99%B; 10–10.5 min, 99%-1%B; 10.5–12 min, 1%B, the flow rate was 0.3 mL/min. The injection volume was 5 μL.
After separation, mass spectra were detected and analyzed using an electrospray ion source (ESI) in positive and negative ionization modes. Ion source parameters are set as follows: The capillary voltage is 3.0 kV, the drying gas temperature is 325°C, the atomized gas pressure is 310 kPa, the drying gas flow is 0.26 mL/min, the desorption gas flow is 600 L/h, the source temperature is 120°C, the desorption temperature is 350°C, and the cone gas flow is 50 L/h.
2.7.4 Data processing
Raw data were exported by Masslynx4.1 (Waters, United States) software, after which the data were imported into SIMCA 14.1 statistical software (Umetrics Corporation, Sweden) for multivariate statistical analysis. Subsequently, SPSS 26.0 was applied for statistical tests and the appropriate test was selected to determine whether metabolites changed significantly (p < 0.05) in the statistical analysis. The metabolites screened by the above analysis are input into the MetaboAnalyst platform(https://www.metaboanalyst.ca/) for cluster analysis and Mate-Pa analysis.
2.8 Network pharmacology
2.8.1 Blood-inlet-components in YYFZ
Based on the incoming components of YYFZ identified in our previous study (Supplementary Material), the information on these components was retrieved in TCMSP (http://lsp.nwsuaf.edu.cn/tcmsp.php) for relevant target information. On the other hand, the molecular structure formula was drawn in Chemdraw, saved in.sdf format and uploaded to PharmaMapper (http://lilab-ecust.cn/pharmmapper/index.html) for predicting potential targets, and SwissTargetPrediction (http://www.swisstargetprediction.ch/) for predicting the results. As a supplement, the targets of all components were combined and de-weighted after conversion to Genename by UniProt (https://www.UniProt.org/) (UniProt Consortium, 2013).
2.8.2 Establishment of CHD disease target library
Search by Genecards (https://www.genecards.org/), TTD data (http://bidd.nus.edu.sg/group/cjttd/) (Wang et al., 2020), OMIM (https://omim.org/) (Amberger and Hamosh, 2017) databases The targets related to the pathogenesis of “chronic heart disease” were combined by removing duplicate targets from several databases.
2.8.3 Visualization and analysis of active ingredient-target-disease network
PPI (protein-protein interaction) maps were created through protein interaction analysis via STRING database and imported into Cytoscape 3.8.2 for intra-network visualization of PPI data between potential therapeutic targets of diseases and mapping of active ingredient-target-disease networks.
KEGG pathway enrichment analysis was performed on the network nodes using DAVID database to obtain the functional pathways involved in disease modulation by each component-acting target of YYFZ bulk, and the enriched pathways were visualized using the advanced bubble map function of the Omicshare platform.
3 Results
3.1 Pharmacodynamic study of YYFZ on CHD rat model
3.1.1 General conditions of rats
During the modeling of CHD rats, the rats in the NS group grew naturally, fed and drank normally, and were in good condition, and there was no death. During the modeling process, the rats’ body weight growth slowed down significantly, and they lost their body weight and hair. After the fourth week, the rats gradually appeared ascites, abdominal distension, panting and shortness of breath, and individual rats lost their hair in flakes; until the seventh week, the rats had obvious abdominal distension, diarrhea, watery stool, filthy perianal area, shortness of breath, vertical hair, loose fur, lack of blood in the feet and paws and ears, mental inactivity, and increased secretion around the eyes and nose. After the start of administration, the condition of the rats in the Captopril group, YYFZ-H, YYFZ-M, and YYFZ-L improved, while no improvement was observed in the Model group. In addition, hepatomegaly was seen in the Model group at autopsy, and some of them had renal edema with mesenteric adhesions.
3.1.2 Weight gain of rats
After the beginning of modeling, it can be seen from (Figure 1A) that the weight gain of rats in the model group slowed down significantly, and there was a significant difference in body weight after the second injection (p < 0.01) until the end of the experiment. During the modeling period, the body weight of the NS group was stable, while that of the model group was slow, and occasionally fluctuated in the later period, but there was almost no obvious increase. The rats in the administration group showed no significant difference in body weight during the modeling period compared to NS (p < 0.05) and showed a significant increase after administration (p < 0.05). At the late stage of the experiment, the weight changes of rats were difficult to directly reflect the drug effect because of the severe abdominal ascites of some rats. Therefore, in this study, the therapeutic effect of the drug was reflected by calculating the body weight growth rate of rats, and the body weight growth rate % = (body weight without ascites—body weight before modeling)/body weight before modeling * 100%. It can be seen from Figure 1C that the body weight growth rate of the Model group and drug administration group was significantly different from that of the NS group (p < 0.05), while compared with the Model group, the body weight of Captopril group, YYFZ-H group and YYFZ-M group was significantly increased (p < 0.05).
FIGURE 1. (A) Weight change trend of rats; (B) Ascites of rats in each group; (C) Weight growth of rats in each group; (D) Heart index of rats in each group.
(Compared with NS group, n = 10, ‾x ± SD *: p<0.05, **: p < 0.01; compared with Model group, +: p < 0.05, ++: p < 0.01)
3.1.3 Ascites in rats
At the late stage of modeling, some rats showed significant signs of ascites; the ascites condition of the rats in the administered group did not continue to worsen as the administration treatment began. As seen in Figure 1B, the ascites of rats in Captopril, YYFZ-H and YYFZ-M groups appeared significantly different from the Model group (p< 0.05), and there was no significant difference in the low-dose group, suggesting that Captopril and high and medium doses of YYFZ have better efficacy.
3.1.4 Rat heart function
At the end of modeling, to determine the success of the model, we performed cardiac function tests on the rats. As seen in Table 1, the EF and FS in the Model group were significantly lower than those in the NS group (p < 0.05), indicating the success of the model. The rats were grouped after modeling, and cardiac function tests were performed after grouping. As seen in Table 2, compared with the NS group, LVPW; s and LVPW; d was significantly lower in the Model and drug administration groups (p < 0.05), with LVPW; d significantly higher in the Captopril group than in the Model group (p < 0.05). IVSs and IVSd in the Model group and YYFZ-H, YYFZ-M and YYFZ-L groups were significantly lower than those in the NS group (p < 0.05), with significant differences between IVSs in the YYFZ-H and YYFZ-M groups and the Model group (p < 0.05). The IVSs and IVSd in the Captopril group were not statistically different from those in the NS group and were significantly higher than those in the Model group (p < 0.05).
TABLE 1. Comparison of cardiac structural and functional indices between modeled rats and blank rats.
TABLE 2. Cardiac indexes (LVPW;s, LVPW;d, IVSs, IVSd) of rats in each group before drug administration.
Cardiac function tests were performed on the rats at the end of the dosing. As shown in Table 3 and Figure 2. The EF and FS in the Model group were significantly lower than those in the NS group (p < 0.05), indicating that the systolic-diastolic function of the heart was reduced in the Model group. The EF and FS in the Captopril and YYFZ-H groups were significantly higher than those in the Model group. Compared with the NS group, LVID;d was significantly lower (p < 0.05) and LVID;s was significantly higher (p < 0.05) in the Model group. LVID;d, LVID;s in the YYFZ-M and YYFZ-L groups were significantly lower than those in the NS and Model groups (p < 0.05). LVID;d, LVID;s in the Captopril and YYFZ-H groups were not significantly different from the NS group, but were significantly different from the Model group (p < 0.05).
TABLE 3. Comparison of cardiac function indices (FS, EF, LVID;s, LVID;d) in rats in each group after drug administration.
FIGURE 2. Echocardiography of rats after treatment ((A): NS group; (B) Model group; (C) Captopril group; (D)YYFZ-H group; (E) YYFZ-M group; (F) YYFZ-L group).
3.1.5 Rat cardiac index
After weighing and comparing, we found that there was no significant difference in the heart weight of the rats in each group, presumably related to the larger size of the rats in the NS group. After conversion into the cardiac index, it can be seen from Figure 1D that the cardiac index of the Model group and each dosing group was significantly higher than that of the NS group (p < 0.05), and the mean value of the cardiac index of each dosing group was lower than that of Model group, but only Captopril group and YYFZ-H group were significantly lower than that of Model group (p < 0.05), and the difference between the cardiac index of the remaining dosing groups and Model group was not The difference between the cardiac index of the remaining administration groups and the Model group was not statistically significant.
3.1.6 Detection of serum BNP and TNF-α in rats
To assess the rat modeling situation, this study used Elisa kits to detect BNP and TNF-α contents in rat serum, as shown in Figure 3C, there was a significant increase in BNP content in the serum of model rats about 2.16 times higher than that of NS group, and the difference was statistically significant (p < 0.01); TNF-α in model rats was about 1.6 times higher than that of NS rats, and the difference was statistically significant (p< 0.01).
FIGURE 3. (A) Reduction of BNP after treatment; (B) Reduction of TNF after treatment; (C) BNP and TNF-α before treament.
After the administration, the serum BNP and TNF-α levels in each group are shown in Figures 3A, B. The serum BNP levels in the Model, YYFZ-M and YYFZ-L groups were significantly higher than those in the NS group (p < 0.05), with the serum BNP in the YYFZ-M group being significantly lower than that in the Model group (p < 0.05) and the YYFZ-L group not significantly different from that in the Model group. The differences between the Captopril and YYFZ-H groups were statistically significantly lower than those in the Model group (p < 0.05) and the NS group (p < 0.05). The TNF-α levels in the Captopril, YYFZ-H and YYFZ-M groups were significantly lower than those in the Model group (p < 0.05), and there was no significant difference with the NS group. (Compared with the NS group, n = 10, ‾x ± SD *: p <0.05, **: p<0.01; compared with the Model group, +: p<0.05, ++: p<0.01).
3.1.7 Pathological results of myocardial tissue in rats
The results of histopathological HE (hematoxylin-eosin) staining of rat myocardium are shown in Figure 4: In the NS group, myocardial cells were clearly arranged and tightly packed, with a few myocardial transverse breaks, and vascular proliferation and inflammatory cell infiltration were not obvious. The myocardial cells in the Model group were unevenly stained, with relatively blurred transverse lines, reduced volume of some cells, and deepened staining of the envelope; the myocardial interstitium was obviously edematous, with loosely arranged cells and fat vacuoles within the cells (yellow arrows), an increased number of small blood vessels between the myocardium (red arrows), and a small number of lymphocytes between the myocardium. Captopril group and YYFZ-H group: myocardial cell transverse lines were clear, there was a small amount of edema in myocardial interstitium, and the area of lax cell arrangement was reduced compared with the model group. YYFZ-M group and YYFZ-L group: myocardial cell transverse lines were clear, the cell volume was reduced, the arrangement was more laxer, there was edema in the cell interstitium, and the number of blood vessels was increased.
FIGURE 4. Pathological sections of myocardial tissue of rats in each group ((A): NS group; (B) Model group; (C) Captopril group; (D) YYFZ-H group; (E) YYFZ-M group; (F) YYFZ-L group).
3.2 Metabolomics analysis results
3.2.1 Methodological investigation
After collecting the data by UPLC-Q-TOF/MS analysis, we randomly selected 20 peaks from the obtained QC sample profiles (as shown in Figure 5A) and calculated the RSD values of their peak areas and retention times. RSD<16.17%, retention time RSD<0.92%; sample stability study in the sample at 0 h, 6 h, 12 h, 18 h, 24 h time points in the peak area RSD<12.45%, retention time RSD<0.59%. The above indicates that the instrument precision, method precision and sample stability are good, indicating that the method is reliable and can be followed up.
FIGURE 5. (A) BPI plots of QC samples; (B) Marker metabolic pathway analysis (A: phenylalanine, tyrosine and tryptophan biosynthesis; B: linolenic acid metabolism; C: phenylalanine metabolism; D: arachidonic acid metabolism); (C) intergroup volcano plot analysis between NS and Model groups.
3.2.2 Metabolomics data pre-processing and multivariate statistical analysis
The raw data were exported from Masslynx4.1 software, and the data were subjected to a normality test and chi-square test after 80% modification, and suitable methods were selected to compare the NS and Model groups between groups, and p-values and fold change (FC) were calculated, and the substances and p-values and log2 (FC) values were entered into the Wukong data analysis cloud platform for volcano plot analysis to visually display the NS and Model groups group in the different substances. As shown in Figure 5C, the horizontal coordinate is log2 (FC), the greater the difference the more distant the metabolite is distributed, and the vertical coordinate is the p-value, the greater the difference the more distant the substance is from the horizontal axis, so the substance distributed in the upper left right corner is usually considered as a potential difference marker. Each of these points represents data, and the gray color is for substances with no significant changes, while the red and blue colors are for substances with significant elevated changes, which are potential differential substances that we need to pay attention to in the subsequent processing.
Multivariate statistical analysis was used to discriminate the changes of metabolites within the plasma of healthy and model rats. For this, we first set up an unsupervised PCA model as shown in Figure 6C. Each point in the figure represents a sample, respectively, and the color is used to distinguish the rat groups, and the farther the distance between the points of different groups indicates the greater difference. We can see that there is a certain convergence in the distribution of NS and Model groups, and Model and drug administration groups in the PCA model. Therefore, to further differentiate, we further established a supervised OPLS-DA model (Figure 6A), and we can see that the NS group and Model group are distributed on both sides in the OPLS-DA model, with obvious distinction.
FIGURE 6. (A) OPLS-DA; (B) CHD marker ROC analysis (C) PCA.
3.2.3 Identification of CHD biomarkers
The SPSS26.0 software was used to test for normality and the chi-square test for the differential substances. t-test and Mann-Whitney u-test were selected to test for significance according to whether the data distribution was normal and chi-square, and the differential substances with significant changes (p < 0.05) were considered as differential markers of heterogeneous disease and treatment. Finally, we obtained 732 NS-CHD differential substances, respectively. We entered m/z values in the HMDB database to search, selected [M + H]+, [M + K]+, [M + Na]+ in ion addition mode, and identified metabolites by MS/MS analysis, metabolite database information, and literature information. We finally identified 19 markers, which were leukotriene F4, arachidonic acid, lysoSM(d18:0), tryptophan, phenylpyruvic acid, S-lactoylglutathione, phenylalanine, ornithine, nicotinamide-N-oxide, fumaric acid, adipic acid, D-pantothenic acid, PC(14:1/18:4), tetradecanedioic acid, 3-deoxy glucosone, α-linolenic acid, L-octanoylcarnitine, lysoPE (20:1(11Z)/0:0), lysoPC (17:0/0:0), see Table 4 for specific information. 16 markers were upregulated and 3 markers were downregulated in the Model group.
TABLE 4. Differential metabolite information.
3.2.4 Hierarchical clustering analysis of CHD biomarkers
To observe the changes of markers in each group more intuitively, this study used hierarchical clustering analysis as shown in Figure 7, which can see the distribution of biomarkers in each group. The shades of color in the figure respond to the magnitude of the values, each row represents a metabolite, and each column represents the content of the whole group of samples; the left bifurcation is the cluster analysis of the substances, where the more clustering levels of the substances indicate the higher similarity of the substances, and there may be similar variations in the source and metabolic pathways. We can observe a significant change in plasma marker content in the Model group compared to the NS group, suggesting that the markers have some discriminatory ability. Metabolite levels in the captopril and YYFZ-H groups were closer to those in the NS group, and metabolite trends were more consistent, suggesting that high-dose YYFZ significantly modulates CHD disease-related markers.
FIGURE 7. Heat map of serum metabolite changes from minimum (dark blue) to maximum (dark red) in NS, Model, YYFZ-H, YYFZ-M, YYFZ-L and Captopril groups.
3.2.5 CHD biomarker ROC analysis
ROC curves are widely used for the evaluation of the sensitivity and specificity of markers and can screen for more diagnostic biomarkers. After the above study, we identified markers with diagnostic potential by ROC analysis of the above 19 biomarkers by SPSS 26.0, and the area under the curve was used to evaluate the diagnostic ability of the markers. Usually, an area under the curve>0.5 indicates good discriminatory ability. As shown in Figure 6B, we can see that each curve represents a substance and the area under the curve is distributed between 1 and 0.87 (95% confidence interval), indicating that the 19 markers have the diagnostic ability.
3.2.6 CHD biomarker pathway analysis
To further speculate on the mechanism of CHD metabolic disorders, this study performed metabolic pathway analysis (MetPA) on 19 markers. The MetPA database is a visual metabolic pathway analysis database (www.metaboanalyst.ca), which combines pathway enrichment analysis and topological analysis to assist in screening relevant metabolic pathways. We obtained 16 metabolic pathways (Figure 5B), including phenylalanine, tyrosine and tryptophan biosynthesis, alpha-linolenic acid metabolism, phenylalanine metabolism, and arachidonic acid metabolism are the main metabolic pathways involved, suggesting that lipid metabolism and amino acid metabolism are the pathways to focus on in CHD.
3.3 Network pharmacology
3.3.1 YYFZ blood entry components and CHD targets
The database was used to predict the YYFZ inlet component targets, and finally, 228 YYFZ component targets were obtained. 1488 CHD disease targets were obtained. Venny enrichment analysis was performed between CHD inlet component targets and disease targets, as shown in Figure 8A, and 44 YYFZ targets for CHD was obtained.
FIGURE 8. (A) Venn diagram of YYFZ component and CHD target; (B) PPI diagram; (C) Drug-active component-target network expressing the interaction between YYFZ inlet component, protein target and CHD pathway (red-inlet component, blue-target, yellow-pathway).
3.3.2 Visualization and analysis of active ingredient-target-disease network
The potential targets of CHD and YYFZ were imported into the String database for protein interaction analysis and the PPI network was constructed, and the results showed 42 nodes with 176 edges (Figure 8B). The results in.tsv format were imported into Cytoscape 3.8.2 software for network analysis, MAPK14, EGFR, HSP90AA1, MTOR, ESR1, and IGF1 had high Dgree values, which might be the key targets of YYFZ for CHD treatment.
The Cytoscape software was used to construct the drug-active-ingredient-target network (Figure 8C). The network contains 70 nodes and 90 edges, and the larger the node degree value, the more important it is in the network. The core nodes were selected based on the network topological features such as node degree values. Among them, senbusine B, sitosterol, hypaconitine, phenylalanine, and deoxyaconitine may be the pharmacodynamic components of YYFZ for CHD.
To further discover the mechanism of action of YYFZ in the treatment of CHD, GO enrichment analysis was performed using the DAVID database for the common targets obtained above (Figure 9), and a total of 162 enrichment results were obtained for biological process (BP) screening at p < 0.05, mainly involving positive regulation of protein kinase B signaling, positive regulation of smooth muscle cell proliferation, phosphatidylinositol 3-kinase signaling. Cellular component (CC) obtained 36 enrichment results, mainly involving the extracellular region, lysosomal lumen, perinuclear region of cytoplasm, etc.; molecular function (MF) obtained 37 enrichment results, mainly involving identical protein binding, Enzyme binding, Protein homodimerization activity, etc. The KEGG pathway analysis is shown in Figure 10, excluding the pathways unrelated to CHD. signaling pathway, PI3K-Akt signaling pathway, MAPK signaling pathway, etc.
FIGURE 9. GO function enrichment map.
FIGURE 10. KEGG pathway enrichment analysis.
4 Discussion
In this experiment, the CHD model of rats was established by intraperitoneal injection of adriamycin, and the rats were treated with captopril and different doses of YYFZ, and the efficacy of YYFZ was comprehensively evaluated by observing the biochemical indexes and cardiac function indexes of rats and myocardial histopathological sections. Adriamycin belongs to anthracycline antibiotics, which are widely used in clinical practice for the treatment of various malignant tumors (Sabatino et al., 2020; Younis et al., 2021; Xu et al., 2022). Cardiotoxicity is one of its main side effects, which can be manifested as irreversible dose-dependent cardiomyopathy and CHD, so a heart failure model was prepared using this side effect. TNF-α is an inflammatory factor that plays an important role in promoting myocardial remodeling and inhibiting myocardial contraction, as well as in increasing endothelial and myocardial apoptosis, and there is abundant evidence that inhibition or reduction has a protective effect in heart failure models (Lee et al., 2019; Reina-Couto et al., 2021; Szabo et al., 2021; Zhong et al., 2022). BNP is a marker of cardiac insufficiency, and many studies have shown that plasma BNP levels are significantly elevated in heart failure patients, which can be used for the diagnosis of heart failure (Ichiki et al., 2013; McDonald et al., 2018; Rørth et al., 2020). We found that the YYFZ-H group had a better therapeutic effect from the perspective of biochemical indexes, significantly regulating serum BNP and TNF-α and improving cardiac function (p < 0.05). The pathological results showed that the YYFZ dose groups could improve myocardial cell disorder and reduce inflammatory cell infiltration. Echocardiographic determination of cardiac function in rats has been widely accepted because it provides a convenient, reliable, and non-invasive method (Zeng et al., 2019; Quagliariello et al., 2021; Wang X. et al., 2022). Generally, EF is measured as an indicator of cardiac function, FS is also one of the commonly used indicators, and the magnitude of the two values shows a correlation. In the present study, we observed the long axis of the sternum in rats and used EF as an indicator to determine the modeling status of rats and the effect of YYFZ treatment. We found that EF and FS were significantly retraced in the administered group, and there were some differences in the mean values of other functional indicators LVPW;s, LVPW;d, LVID;s, LVID;d, IVS;s, IVS;d, etc. Moreover, the cardiac index of rats in the YYFZ-H and YYFZ-M groups was lower than that of the model group, and ascites were also significantly reduced, suggesting that YYFZ could improve cardiac function and alleviate the symptoms of CHD in rats with CHD.
In this study, we used non-targeted metabolomics techniques to study serum metabolites in CHD rats and identified 19 CHD-related metabolites, seven of which were regulated by YYFZ including S-lactoylglutathione, phenylalanine, fumaric acid, adipic acid, tetradecanedioic acid, α-linolenic acid, and lysoPE (20:1(11Z)/0:0), and the metabolic pathways involved mainly include four categories: amino acid metabolism (phenylalanine, tyrosine and tryptophan biosynthesis; phenylalanine metabolism; arginine biosynthesis; alanine, aspartate and glutamate metabolism; arginine and proline metabolism; tryptophan metabolism; tyrosine metabolism), lipid metabolism (α-linolenic acid metabolism, arachidonic acid metabolism, biosynthesis of unsaturated fatty acids, glycerophospholipid metabolism), energy metabolism (pyruvate metabolism, pantothenate and CoA biosynthesis, citrate cycle (TCA cycle)), glutathione metabolism.
It has been confirmed that myocardial infarction and heart failure can lead to abnormalities in amino acid metabolism such as branched-chain amino acids, taurine or glutamine, and abnormalities in the metabolism of certain amino acids, such as branched-chain amino acids, can activate the mTOR signaling pathway, which accelerates myocardial remodeling and leads to the development of heart failure after myocardial infarction. In addition, elevated leucine levels may activate mTOR, and inhibition of mTOR in a heart failure model improves cardiac function. When mitochondrial dysfunction and phenylalanine are increased in the heart failure state, tetrahydrobiopterin consumption is increased, resulting in decreased nitric oxide production and causing heart failure (Nishijima et al., 2011; Hong et al., 2019). Phenylalanine may also be associated with increased aromatic amino acids, and increased protein breakdown in muscle and impaired liver function can cause accumulation of aromatic group amino acids. Dysregulation of the hepatic urea cycle can lead to a further increase in ornithine, a precursor to the formation of polyamines. Polyamines are also involved in the development of cardiac disease processes, and in models of myocardial hypertrophy (Chen et al., 2019). This further reflects the regulatory role of amino acid metabolism in heart failure, but the exact mechanisms still need to be investigated in depth (Wang et al., 2016; Dimou et al., 2022).
Disorders of lipid metabolism play an extremely important role in the process of CHD. It has been shown that patients with CHD often have abnormal lipid levels, and that disorders of lipid metabolism further contribute to the development of CHD. In this experiment, linoleic acid was significantly elevated in the model group compared to the blank group, and the molecular structure of linoleic acid contains double bonds that are susceptible to oxidative stress, causing vasodilator dysfunction and endothelial damage (Asselin et al., 2013; Djuricic and Calder, 2021). Phospholipids are the main constituent molecules of biological membrane structure and assist in the function of the organism at the cellular level. Arachidonic acid is a widely distributed unsaturated fatty acid in living organisms, and in pathological states, it is broken down into biologically active derivatives that are involved in inflammatory responses, apoptosis, and other biological processes (Sonnweber et al., 2018). Prostaglandins (PGs) are produced by the metabolism of arachidonic acid through the COX pathway and are further metabolized to produce a variety of active substances that regulate the process of vascular remodeling and generation and can control vascular tone and influence the course of CHD (Mitchell et al., 2021; Wan et al., 2021).
The heart is the organ of the body with the greatest energy demand, and the diastolic and contractile movements of the heart muscle depend on the energy produced by the cells. Under normal conditions, 60%–80% of the heart’s energy is derived from fatty acid oxidation and 10%–20% from glucose metabolism. In contrast, when cardiac function is abnormal, the metabolic substrate is changed from glucose to fatty acids, resulting in increased levels of free fatty acids in serum and myocardial tissue. It has been shown that lipid deposition is observed in rat myocardium high in free fatty acids and that free fatty acids may contribute to impaired cardiac function via PI3K-Akt-GLUT4 and AMPK-eNOS pathways (Han et al., 2018). Acylcarnitine is mainly found in muscle tissues such as cardiac muscle and is involved in fatty acid metabolism and amino acid metabolism. Short-chain acylcarnitines are associated with branched-chain and aromatic amino acid metabolism, and the buildup of long-chain acylcarnitines results from impaired fatty acid oxidation (Roussel et al., 2015). One study found that acylcarnitine decreased in a guinea pig model of compensatory hypertrophy, while fatty acids gradually increased during the process. Long-chain acyl CoA dehydrogenase knockout mice were cardiac hypertrophy accompanied by elevated triglyceride levels, and after carnitine administration, triglyceride levels returned to normal in the myocardium of long-chain acyl CoA dehydrogenase knockout mice (Foster et al., 2016). The above suggests that amino acid metabolism, lipid metabolism and energy metabolism are closely related to the pathogenesis of CHD and that YYFZ acts to inhibit inflammation and reduce intravascular lipid accumulation by regressing lipid and amino acid metabolism.
Hypaconitine, deoxyaconitine and senbusine B are all diterpenoid diester alkaloids with analgesic, anti-inflammatory and antipyretic effects (Zhang et al., 2022). YYFZ is composed of aconitum and Coicis seed. Aconitoid alkaloids in aconitum aconitine are the main active components of YYFZ, while aconitine alkaloids can be hydrolyzed under heating conditions. It is speculated that hypaconitine, deoxyaconitine and senbusine B are the hydrolyzed alkaloids that play a role. Sitosterol is the most abundant sterol in Coix seed, which can interfere with a variety of cell signaling pathways, including cell cycle, cell apoptosis, survival, invasion, proliferation, metastasis, anti-inflammatory, angiogenesis, and cardiac protection (Khan et al., 2022).MAPK14 (p38α), a member of the mitogen-activated protein kinase (MAPK) family, is the most abundant and well-characterized of the four isoforms of p38MAPK and plays a central role in the initiation of stress-activated pro-inflammatory responses (Cuenda and Rousseau, 2007). p38MAPK14 acts downstream on transcription factors and protein kinases to regulate cellular stress and inflammatory responses. It has been demonstrated that p38MAPK plays an important role in myocarditis and myocardial remodeling (Niu et al., 2017; Turner and Blythe, 2019). EGFR is a receptor-type tyrosine kinase, which is important for embryonic heart development and maintenance of adult heart function. EGFR signaling pathway is a complex network system, including three main pathways: Ras/Raf/MAPK pathway, PI3K/AKT pathway, JAK pathway, and STAT pathway (Iwamoto and Mekada, 2006). mTOR is an atypical serine/threonine protein kinase, which can participate in gene transcription and protein expression by phosphorylating its downstream target proteins, and then affect biological activities such as apoptosis. mTOR signaling pathway is one of the effective ways to treat heart failure (Shi et al., 2020). IGF1 is a key hormone that regulates the growth of cardiomyocytes and physiological cardiac hypertrophy. IGF1 is closely related to cardiac hypertrophy and heart failure. In cardiomyocytes, after IGF1 binds to its receptor, the receptor itself is phosphorylated and further activates the MAPK pathway and PI3K pathway (Tao et al., 2017). The Ras family of small guanosine triphosphate (GTP)-binding proteins (G proteins) is one of the major components of intracellular signaling required for normal heart growth and also plays a key role in the development of heart failure. Ras regulates multiple downstream signaling pathways mainly involving MAPK as well as the PI3K/Akt/mTOR pathway (Ramos-Kuri et al., 2021). Overactivation of the JNK and p38 signaling pathways play a key role in the CHD Ras can activate the JNK and p38 signaling pathways through the PI3K/AKT/mTOR pathway. The PI3K-Akt signaling pathway has an important role in the pathogenesis of heart disease, not only regulating the survival and function of cardiomyocytes, but also influencing the proliferation, migration, and apoptosis of vascular smooth muscle cells through the regulation of the pathogenesis and development of CHD (Magaye et al., 2021; Qin et al., 2021; Sun G. et al., 2021). These studies suggest that Ras/MEK/ERK and PI3K/Akt/mTOR signaling control the interconnection of several levels of protein synthesis, cardiomyocyte growth, and remodeling. Ras and its downstream pathways (MAPK pathway, PI3K pathway) may be important pathways for CHD treatment by YYFZ.
5 Conclusion
In this study, the pharmacological effects and pharmacodynamic mechanisms of YYFZ in the treatment of adriamycin-induced CHD rats were investigated based on metabolomics combined with network pharmacology. The YYFZ group was able to improve the disorder of cardiomyocyte arrangement and reduce inflammatory cell infiltration, and the cardiac index of rats in the YYFZ-H and YYFZ-M groups was lower than that of the Model group, suggesting that YYFZ was able to improve the cardiac function and alleviate the symptoms of CHD in CHD rats. Metabolomics identified 19 biomarkers related to the pathogenesis of CHD, mainly involving the improvement of amino acid metabolic pathways, fatty acid metabolism, and energy metabolism. In addition, network pharmacology analysis showed that YYFZ may act through MAPK14, EGFR, HSP90AA1, MTOR, ESR1, IGF1, and other protein targets in the treatment of CHD. Combining the results of metabolomics and network pharmacology pathway enrichment analysis suggested that YYFZ may further regulate PI3K/Akt pathway and MAPK pathway by modulating amino acid metabolism and fatty acid metabolism. YYFZ, as a classical TCM formula, has been used to treat CHD from multiple links and targets.
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Ethics statement
The animal study was reviewed and approved by the Animal Ethics Committee of Tianjin University of Traditional Chinese Medicine.
Author contributions
YmW and XL conceived, designed the experiments. XL carried out the experiments and wrote the manuscript. MQ, XkL, FZ, YyW, JW, LS, and SF participated in the experiments. YfL and YbL reviewed and edited the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (No. 81903938) and Tianjin Development Program for Innovation.
Conflict of interest
Author MQ was employed by TIPRHUYA Advancing Innovative Medicines Ltd.
The remaining 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/fmolb.2023.1203208/full#supplementary-material
References
Amberger, J. S., and Hamosh, A. (2017). Searching online mendelian inheritance in man (OMIM): A knowledgebase of human genes and genetic phenotypes. Curr. Protoc. Bioinforma. 58, 1. doi:10.1002/cpbi.27-1.2.12
Asselin, C., Ducharme, A., Ntimbane, T., Ruiz, M., Fortier, A., Guertin, M. C., et al. (2013). Circulating levels of linoleic acid and HDL-cholesterol are major determinants of 4-hydroxynonenal protein adducts in patients with heart failure. Redox Biol. 2, 148–155. doi:10.1016/j.redox.2013.12.009
Buergel, T., Steinfeldt, J., Ruyoga, G., Pietzner, M., Bizzarri, D., Vojinovic, D., et al. (2022). Metabolomic profiles predict individual multidisease outcomes. Metabolomic profiles predict Individ. multidisease outcomes 28 (11), 2309–2320. doi:10.1038/s41591-022-01980-3
Chen, H. Y., Jia, X. L., Zhao, S. Q., Zheng, W. H., Mei, Z. G., Yang, H. W., et al. (2019). Dual role of polyamines in heart ischemia/reperfusion injury through regulation of mitochondrial permeability transition pore. Sheng Li Xue Bao 71 (5), 681–688. doi:10.13294/j.aps.2019.0057
Chen, X., Chen, Y., Xie, S., Wang, X., Wu, Y., Zhang, H., et al. (2022). The mechanism of Renshen-Fuzi herb pair for treating heart failure-Integrating a cardiovascular pharmacological assessment with serum metabolomics. Front. Pharmacol. 13, 995796. doi:10.3389/fphar.2022.995796
Cuenda, A., and Rousseau, S. (2007). p38 MAP-kinases pathway regulation, function and role in human diseases. Biochim. Biophys. Acta 1773 (8), 1358–1375. doi:10.1016/j.bbamcr.2007.03.010
Cui, H., Chen, Y., Li, K., Zhan, R., Zhao, M., Xu, Y., et al. (2021). Untargeted metabolomics identifies succinate as a biomarker and therapeutic target in aortic aneurysm and dissection. Eur. Heart J. 42 (42), 4373–4385. doi:10.1093/eurheartj/ehab605
Dimou, A., Tsimihodimos, V., and Bairaktari, E. (2022). The critical role of the branched chain amino acids (BCAAs) catabolism-regulating enzymes, branched-chain aminotransferase (BCAT) and branched-chain α-keto acid dehydrogenase (BCKD), in human pathophysiology. Int. J. Mol. Sci. 23 (7), 4022. doi:10.3390/ijms23074022
Djuricic, I., and Calder, P. C. (2021). Beneficial outcomes of omega-6 and omega-3 polyunsaturated fatty acids on human health: An update for 2021. Nutrients 13 (7), 2421. doi:10.3390/nu13072421
Du, Z. Y., Shu, Z. L., Li, C., Song, X. M., Ma, X. L., Liao, L. X., et al. (2021). Baoyuan decoction alleviates myocardial infarction through the regulation of metabolic dysfunction and the mitochondria-dependent caspase-9/3 pathway. Acupunct. Herb. Med. 1 (1), 49–58. doi:10.1097/hm9.0000000000000003
Foster, D. B., Liu, T., Kammers, K., O'Meally, R., Yang, N., Papanicolaou, K. N., et al. (2016). Integrated omic analysis of a Guinea pig model of heart failure and sudden cardiac death. J. Proteome Res. 15 (9), 3009–3028. doi:10.1021/acs.jproteome.6b00149
Han, L., Liu, J., Zhu, L., Tan, F., Qin, Y., Huang, H., et al. (2018). Free fatty acid can induce cardiac dysfunction and alter insulin signaling pathways in the heart. Lipids Health Dis. 17 (1), 185. doi:10.1186/s12944-018-0834-1
Hong, F. F., Liang, X. Y., Liu, W., Lv, S., He, S. J., Kuang, H. B., et al. (2019). Roles of eNOS in atherosclerosis treatment. Inflamm. Res. 68 (6), 429–441. doi:10.1007/s00011-019-01229-9
Ichiki, T., Huntley, B. K., and Burnett, J. C. (2013). BNP molecular forms and processing by the cardiac serine protease corin. Adv. Clin. Chem. 61, 1–31. doi:10.1016/b978-0-12-407680-8.00001-4
Iwamoto, R., and Mekada, E. (2006). ErbB and HB-EGF signaling in heart development and function. Cell. Struct. Funct. 31 (1), 1–14. doi:10.1247/csf.31.1
Khan, Z., Nath, N., Rauf, A., Emran, T. B., Mitra S, etc., Islam, F., et al. (2022). Multifunctional roles and pharmacological potential of β-sitosterol: Emerging evidence toward clinical applications. Chem. Biol. Interact. 365, 110117. doi:10.1016/j.cbi.2022.110117
Lee, H. C., Shiou, Y. L., Jhuo, S. J., Chang, C. Y., Liu, P. L., Jhuang, W. J., et al. (2019). The sodium-glucose co-transporter 2 inhibitor empagliflozin attenuates cardiac fibrosis and improves ventricular hemodynamics in hypertensive heart failure rats. Cardiovasc Diabetol. 18 (1), 45. doi:10.1186/s12933-019-0849-6
Leung, A. Y. L., Chen, H., Jia, Z., Li, X., and Shen, J. (2021). Study protocol: Traditional Chinese Medicine (TCM) syndrome differentiation for heart failure patients and its implication for long-term therapeutic outcomes of the Qiliqiangxin capsules. Chin. Med. 16 (1), 103. doi:10.1186/s13020-021-00515-1
Li, J., Guasch-Ferré, M., Chung, W., Ruiz-Canela, M., Toledo, E., Corella, D., et al. (2020). The Mediterranean diet, plasma metabolome, and cardiovascular disease risk. Eur. Heart J. 41 (28), 2645–2656. doi:10.1093/eurheartj/ehaa209
Li, L., Li, S. H., Jiang, J. P., Liu, C., and Ji, L. L. (2021). Investigating pharmacological mechanisms of andrographolide on non-alcoholic steatohepatitis (nash): A bioinformatics approach of network pharmacology. Chin. Herb. Med. 13 (3), 342–350. doi:10.1016/j.chmed.2021.05.001
Li, P., Chen, J., Zhang, W., Li, H., Wang, W., and Chen, J. (2019). Network pharmacology based investigation of the effects of herbal ingredients on the immune dysfunction in heart disease. Pharmacol. Res. 141, 104–113. doi:10.1016/j.phrs.2018.12.016
Liao, M., Xie, Q., Zhao, Y., Yang, C., Lin, C., Wang, G., et al. (2022). Main active components of Si-Miao-Yong-An decoction (SMYAD) attenuate autophagy and apoptosis via the PDE5A-AKT and TLR4-NOX4 pathways in isoproterenol (ISO)-induced heart failure models. Pharmacol. Res. 176, 106077. doi:10.1016/j.phrs.2022.106077
Liu, Y. J., Li, Z., Wang, X., Ni, T. Y., Ma, M., He, Y. Y., et al. (2022). Effects of adjuvant Chinese patent medicine therapy on major adverse cardiovascular events in patients with coronary heart disease angina pectoris: A population-based retrospective cohort study. Acupunct. Herb. Med. 2 (2), 109–117. doi:10.1097/hm9.0000000000000028
Magaye, R. R., Savira, F., Hua, Y., Xiong, X., Huang, L., Reid, C., et al. (2021). Attenuating PI3K/Akt-mTOR pathway reduces dihydrosphingosine 1 phosphate mediated collagen synthesis and hypertrophy in primary cardiac cells. Int. J. Biochem. Cell. Biol. 134, 105952. doi:10.1016/j.biocel.2021.105952
McDonald, K., Troughton, R., Dahlström, U., Dargie, H., Krum, H., van der Meer, P., et al. (2018). Daily home BNP monitoring in heart failure for prediction of impending clinical deterioration: Results from the HOME HF study. Eur. J. Heart Fail 20 (3), 474–480. doi:10.1002/ejhf.1053
Meng, H., Du, Z., Lu, W., Wang, Q., Sun, X., Jiang, Y., et al. (2021). Baoyuan decoction (BYD) attenuates cardiac hypertrophy through ANKRD1-ERK/GATA4 pathway in heart failure after acute myocardial infarction. Phytomedicine 89, 153617. doi:10.1016/j.phymed.2021.153617
Mitchell, J. A., Kirkby, N. S., Ahmetaj-Shala, B., Armstrong, P. C., Crescente, M., Ferreira, P., et al. (2021). Cyclooxygenases and the cardiovascular system. Pharmacol. Ther. 217, 107624. doi:10.1016/j.pharmthera.2020.107624
Ni, C., Li, B., Ding, Y., Wu, Y., Wang, Q., Wang, J., et al. (2021). Anti-cancer properties of coix seed oil against HT-29 colon cells through regulation of the PI3K/AKT signaling pathway. Foods 10 (11), 2833. doi:10.3390/foods10112833
Nicholson, J. K. (2021). Molecular phenomic approaches to deconvolving the systemic effects of SARS-CoV-2 infection and post-acute COVID-19 syndrome. Phenomics 1 (4), 143–150. doi:10.1007/s43657-021-00020-3
Nishijima, Y., Sridhar, A., Bonilla, I., Velayutham, M., Khan, M., Terentyeva, R., et al. (2011). Tetrahydrobiopterin depletion and NOS2 uncoupling contribute to heart failure-induced alterations in atrial electrophysiology. Cardiovasc Res. 91 (1), 71–79. doi:10.1093/cvr/cvr087
Niu, L., Li, C., Wang, Z., Xu, H., and An, X. (2017). Effects of the MAPK pathway and the expression of CAR in a murine model of viral myocarditis. Exp. Ther. Med. 13 (1), 230–234. doi:10.3892/etm.2016.3909
Qin, W., Cao, L., and Massey, I. Y. (2021). Role of PI3K/Akt signaling pathway in cardiac fibrosis. Mol. Cell. Biochem. 476 (11), 4045–4059. doi:10.1007/s11010-021-04219-w
Quagliariello, V., De Laurentiis, M., Rea, D., Barbieri, A., Monti, M. G., Carbone, A., et al. (2021). The SGLT-2 inhibitor empagliflozin improves myocardial strain, reduces cardiac fibrosis and pro-inflammatory cytokines in non-diabetic mice treated with doxorubicin. Cardiovasc Diabetol. 20 (1), 150. doi:10.1186/s12933-021-01346-y
Ramos-Kuri, M., Meka, S. H., Salamanca-Buentello, F., Hajjar, R. J., Lipskaia, L., and Chemaly, E. R. (2021). Molecules linked to Ras signaling as therapeutic targets in cardiac pathologies. Biol. Res. 4 (1), 23. doi:10.1186/s40659-021-00342-6
Reina-Couto, M., Pereira-Terra, P., Quelhas-Santos, J., Silva-Pereira, C., Albino-Teixeira, A., and Sousa, T. (2021). Inflammation in human heart failure: Major mediators and therapeutic targets. Front. Physiol. 12, 746494. doi:10.3389/fphys.2021.746494
Roger, V. L. (2021). Epidemiology of heart failure: A contemporary perspective. Circ. Res. 128 (10), 1421–1434. doi:10.1161/CIRCRESAHA.121.318172
Rørth, R., Jhund, P. S., Yilmaz, M. B., Kristensen, S. L., Welsh, P., Desai, A. S., et al. (2020). Comparison of BNP and NT-proBNP in patients with heart failure and reduced ejection fraction. Circ. Heart Fail 13 (2), e006541. doi:10.1161/CIRCHEARTFAILURE.119.006541
Roussel, J., Thireau, J., Brenner, C., Saint, N., Scheuermann, V., Lacampagne, A., et al. (2015). Palmitoyl-carnitine increases RyR2 oxidation and sarcoplasmic reticulum Ca2+ leak in cardiomyocytes: Role of adenine nucleotide translocase. Biochim. Biophys. Acta 1852 (5), 749–758. doi:10.1016/j.bbadis.2015.01.011
Sabatino, J., De Rosa, S., Tammè, L., Iaconetti, C., Sorrentino, S., Polimeni, A., et al. (2020). Empagliflozin prevents doxorubicin-induced myocardial dysfunction. Cardiovasc Diabetol. 19 (1), 66. doi:10.1186/s12933-020-01040-5
Santos-Gallego, C. G., Mayr, M., and Badimon, J. (2022). SGLT2 inhibitors in heart failure: Targeted metabolomics and energetic metabolism. Circulation 146 (11), 819–821. doi:10.1161/CIRCULATIONAHA.122.060805
Shi, B., Huang, Y., Ni, J., Chen, J., Wei, J., Gao, H., et al. (2020). Qi Dan Li Xin pill improves chronic heart failure by regulating mTOR/p70S6k-mediated autophagy and inhibiting apoptosis. Sci. Rep. 10 (1), 6105. doi:10.1038/s41598-020-63090-9
Shi, D., Liu, L., Li, H., Pan, D., Yao, X., Xiao, W., et al. (2022). Identifying the molecular basis of Jinhong tablets against chronic superficial gastritis via chemical profile identification and symptom-guided network pharmacology analysis. J. Pharm. Anal. 12 (1), 65–76. doi:10.1016/j.jpha.2021.01.005
Sonnweber, T., Pizzini, A., Nairz, M., Weiss, G., and Tancevski, I. (2018). Arachidonic acid metabolites in cardiovascular and metabolic diseases. Int. J. Mol. Sci. 19 (11), 3285. doi:10.3390/ijms19113285
Sui, Y., and Xu, D. (2022). Isolation and identification of anti-inflammatory and analgesic polysaccharides from Coix seed (Coix lacryma-jobi L.var. Ma-yuen (Roman) Stapf). Nat. Prod. Res. 30, 1–10. doi:10.1080/14786419.2022.2162896
Sun, C., Zeng, G., Wang, T., Ren, H., An, H., Lian, C., et al. (2021). Astragaloside IV ameliorates myocardial infarction induced apoptosis and restores cardiac function. Front. Cell. Dev. Biol. 9, 671255. doi:10.3389/fcell.2021.671255
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
Szabo, T. M., Frigy, A., and Nagy, E. E. (2021). Targeting mediators of inflammation in heart failure: A short synthesis of experimental and clinical results. Int. J. Mol. Sci. 22 (23), 13053. doi:10.3390/ijms222313053
Tai, C. J., El-Shazly, M., Yang, Y. H., Tsai, Y. H., Csupor, D., Hohmann, J., et al. (2022). The effectiveness of Fuzi in combination with routine heart failure treatment on chronic heart failure patients. J. Ethnopharmacol. 289, 115040. doi:10.1016/j.jep.2022.115040
Tao, L., Jia, L., Li, Y., Song, C., and Chen, Z. (2017). Recent advances of adapter proteins in the regulation of heart diseases. Heart Fail Rev. 22 (1), 99–107. doi:10.1007/s10741-016-9582-3
Tsutsui, H. (2022). Recent advances in the pharmacological therapy of chronic heart failure: Evidence and guidelines. Pharmacol. Ther. 238, 108185. doi:10.1016/j.pharmthera.2022.108185
Turner, N. A., and Blythe, N. M. (2019). Cardiac fibroblast p38 MAPK: A critical regulator of myocardial remodeling. J. Cardiovasc Dev. Dis. 6 (3), 27. doi:10.3390/jcdd6030027
UniProt Consortium (2013). Activities at the universal protein resource (UniProt). Nucleic Acids Res. 42, D191–D198. doi:10.1093/nar/gkt1140
Wan, Q., Kong, D., Liu, Q., Guo, S., Wang, C., Zhao, Y., et al. (2021). Congestive heart failure in COX2 deficient rats. Sci. China Life Sci. 64 (7), 1068–1076. doi:10.1007/s11427-020-1792-5
Wang, H., Yin, H., Zhong, Y., Hu, J., Xia, S., Wang, Z., et al. (2022). Polysaccharides from fermented coix seed modulates circulating nitrogen and immune function by altering gut microbiota. Curr. Res. Food Sci. 5, 1994–2003. doi:10.1016/j.crfs.2022.10.007
Wang, W., Zhang, F., Xia, Y., Zhao, S., Yan, W., Wang, H., et al. (2016). Defective branched chain amino acid catabolism contributes to cardiac dysfunction and remodeling following myocardial infarction. Am. J. Physiol. Heart Circ. Physiol. 311 (5), H1160–H1169. doi:10.1152/ajpheart.00114.2016
Wang, X., Simayi, A., Fu, J., Zhao, X., and Xu, G. (2022). Resveratrol mediates the miR-149/HMGB1 axis and regulates the ferroptosis pathway to protect myocardium in endotoxemia mice. Am. J. Physiol. Endocrinol. Metab. 323 (1), E21–E32. doi:10.1152/ajpendo.00227.2021
Wang, Y., Zhang, S., Li, F., Zhou, Y., Zhang, Y., Wang, Z., et al. (2020). Therapeutic target database 2020: Enriched resource for facilitating research and early development of targeted therapeutics. Nucleic Acids Res. 48 (1), D1031–D1041. doi:10.1093/nar/gkz981
Xu, L., He, D., Wu, Y., Shen, L., Wang, Y., and Xu, Y. (2022). Tanshinone IIA inhibits cardiomyocyte apoptosis and rescues cardiac function during doxorubicin-induced cardiotoxicity by activating the DAXX/MEK/ERK1/2 pathway. Phytomedicine 107, 154471. doi:10.1016/j.phymed.2022.154471
Yan, P., Mao, W., Jin, L., Fang, M., Liu, X., Lang, J., et al. (2020). Crude radix aconiti lateralis preparata (Fuzi) with Glycyrrhiza reduces inflammation and ventricular remodeling in mice through the TLR4/NF-κB pathway. Mediat. Inflamm. 2020, 5270508. doi:10.1155/2020/5270508
Yang, J., Liu, Y., Lu, S., Sun, X., Yin, Y., Wang, K., et al. (2022). Coix seed oil regulates mitochondrial functional damage to induce apoptosis of human pancreatic cancer cells via the PTEN/PI3K/AKT signaling pathway. Mol. Biol. Rep. 9 (7), 5897–5909. doi:10.1007/s11033-022-07371-8
Yao, Y. X., Xu, Y. L., Liu, B. W., Yang, H., Li, S. Y., Zhazo, L. L., et al. (2020). Identification and verification of effective components of Huanghuai for dysfunctional uterine bleeding based on network pharmacology and molecular docking. Chin. Herb. Med. 13 (2), 177–188. doi:10.1016/j.chmed.2020.08.003
Younis, N. N., Salama, A., Shaheen, M. A., and Eissa, R. G. (2021). Pachymic acid attenuated doxorubicin-induced heart failure by suppressing miR-24 and preserving cardiac junctophilin-2 in rats. Int. J. Mol. Sci. 22 (19), 10710. doi:10.3390/ijms221910710
Zeng, Z., Wang, Q., Yang, X., Ren, Y., Jiao, S., Zhu, Q., et al. (2019). Qishen granule attenuates cardiac fibrosis by regulating TGF-β/Smad3 and GSK-3β pathway. Phytomedicine 62, 152949. doi:10.1016/j.phymed.2019.152949
Zhang, J., Li, D., Zhong, D., Zhou, Q., Yin, Y., Gao, J., et al. (2022). Processed lateral root of aconitum carmichaelii debx.: A review of cardiotonic effects and cardiotoxicity on molecular mechanisms. Front. Pharmacol. 13, 1026219. doi:10.3389/fphar.2022.1026219
Zhang, L., Lu, X., Wang, J., Li, P., Li, H., Wei, S., et al. (2017). Zingiberis rhizoma mediated enhancement of the pharmacological effect of aconiti lateralis radix praeparata against acute heart failure and the underlying biological mechanisms. Biomed. Pharmacother. 96, 246–255. doi:10.1016/j.biopha.2017.09.145
Zhong, Y., Chen, L., Li, M., Chen, L., Qian, Y., Chen, C., et al. (2022). Dangshen erling decoction ameliorates myocardial hypertrophy via inhibiting myocardial inflammation. Front. Pharmacol. 12, 725186. doi:10.3389/fphar.2021.725186
Zhou, J., Qian, M., Li, C., Zhou, W., Yao, L., Chen, Y., et al. (2022). Metabolomics and integrated network pharmacology analysis reveal that ginkgolides act as potential active anticancer components by regulating one-carbon metabolism. J. Ethnopharmacol. 298, 115609. doi:10.1016/j.jep.2022.115609
Keywords: metabolomics, biomarker, chronic heart disease, YiYiFuZi powder, pathways
Citation: Wang Y, Li X, Qi M, Li X, Zhang F, Wang Y, Wu J, Shu L, Fan S, Li Y and Li Y (2023) Pharmacological effects and mechanisms of YiYiFuZi powder in chronic heart disease revealed by metabolomics and network pharmacology. Front. Mol. Biosci. 10:1203208. doi: 10.3389/fmolb.2023.1203208
Received: 10 April 2023; Accepted: 12 June 2023;
Published: 23 June 2023.
Edited by:
Ryszard Tomasz Smolenski, Medical University of Gdansk, PolandReviewed by:
Wanling Xuan, University of South Florida, United StatesXiao Li, Shandong Provincial Qianfoshan Hospital, China
Copyright © 2023 Wang, Li, Qi, Li, Zhang, Wang, Wu, Shu, Fan, Li and Li. 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: Yunfei Li, liyunfei@tjutcm.edu.cn; Yubo Li, yaowufenxi001@sina.com
†These authors have contributed equally to this work and share first authorship