Front. Pharmacol.Frontiers in PharmacologyFront. Pharmacol.1663-9812Frontiers Media S.A.111321310.3389/fphar.2023.1113213PharmacologyOriginal ResearchIntegrating strategies of metabolomics, network pharmacology, and experiment validation to investigate the processing mechanism of Epimedium fried with suet oil to warm kidney and enhance yangSun et al.10.3389/fphar.2023.1113213SunE.12†*HuangRan13†DingKe1WangLing1HouJian2TanXiaobin12WeiYingjie12FengLiang4JiaXiaobin14*1The Third Clinical Medical College, Nanjing University of Chinese Medicine, Nanjing, China2Key Laboratory of New Drug Delivery System of Chinese Materia Medica, Jiangsu Academy of Traditional Chinese Medicine, Nanjing, China3Affiliated Hospital of Nanjing University of Chinese Medicine, Nanjing, China4School of Traditional Chinese Pharmacy, China Pharmaceutical University, Nanjing, China
Edited by:Junmin Zhang, Lanzhou University, China
Reviewed by:Thomas Hsueh, Taipei City Hospital, Taiwan
Bing Chen, Fujian Medical University, China
*Correspondence: E. Sun, sune0825@163.com; Xiaobin Jia, jiaxiaobin2015@163.com
These authors have contributed equally to this work and share first authorship
This article was submitted to Drug Metabolism and Transport, a section of the journal Frontiers in Pharmacology
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.
Introduction: Epimedium, a traditional Chinese medicine (TCM) commonly used in ancient and modern China, is one of the traditional Chinese medicines clinically used to treat kidney yang deficiency syndrome (KYDS). There are differences in the efficacy of Epimedium before and after processing, and the effect of warming the kidney and enhancing yang is significantly enhanced after heating with suet oil. However, the active compounds, corresponding targets, metabolic pathways, and synergistic mechanism of frying Epimedium in suet oil to promote yang, remain unclear.
Methods: Herein, a strategy based on comprehensive GC-TOF/MS metabolomics and network pharmacology analysis was used to construct an “active compounds-targets-metabolic pathways” network to identify the active compounds, targets and metabolic pathways involved. Subsequently, the targets in kidney tissue were further validated by real-time quantitative polymerase chain reaction (RT-qPCR). Histopathological analysis with physical and biochemical parameters were performed.
Results: Fifteen biomarkers from urine and plasma, involving five known metabolic pathways related to kidney yang deficiency were screened. The network pharmacology results showed 37 active compounds (13 from Epimedium and 24 from suet oil), 159 targets, and 267 pathways with significant correlation. Importantly, integrated metabolomics and network pharmacologic analysis revealed 13 active compounds (nine from Epimedium and four from suet oil), 7 corresponding targets (ALDH2, ARG2, GSTA3, GSTM1, GSTM2, HPGDS, and NOS2), two metabolic pathways (glutathione metabolism, arginine and proline metabolism), and two biomarkers (Ornithine and 5-Oxoproline) associated with improved kidney yang deficiency by Epimedium fried with suet oil.
Discussion: These finds may elucidate the underlying mechanism of yang enhancement via kidney warming effects. Our study indicated that the mechanism of action mainly involved oxidative stress and amino acid metabolism. Here, we demonstrated the novel strategies of integrating metabolomics and network pharmacology in exploring of the mechanisms of traditional Chinese medicines.
Epimediummetabolomicsnetwork pharmacologyprocessing mechanismkidney yang deficiency syndromeNational Natural Science Foundation of China10.13039/501100001809Jiangsu Provincial Medical Youth Talent10.13039/5011000130591 Introduction
Epimedium, also known as Xianlingpi, is the dried leaves of a Berberis plant Epimedium brevicornum Maxim., Epimedium sagittatum (Sieb.et Zucc.) Maxim., Epimedium pubescens Maxim., or the dried leaves of Epimedium koreanum Nakai, which was first recorded in Shennong Materia Medica Classic and listed as a medium product, with functions of tonifying the liver and kidney, strengthening muscles and bones, and removing wind and dampness (Chinese Pharmacopoeia Commission [CPC], 2020). There are more than 20 processing methods for Epimedium, including stir-frying, wine roasting, salt roasting, ghee roasting, and suet oil roasting, among which suet oil roasting is the most commonly used and has been included in The Chinese Pharmacopoeia 2020 edition. According to traditional processing theory, suet oil is sweet and warm, possessing tonifying effects to treat deficiency syndromes. Epimedium fried with suet oil can play a synergistic role in strengthening its effect of warming the kidney and enriching the yang, mainly for impotence and infertility treatments. Clinical studies showed that Epimedium raw product focuses on dispelling rheumatism, and strengthening muscles and bones.
Kidney yang deficiency syndrome (KYDS) is a common syndrome type in Chinese medicine consultation, and it is also the basic syndrome type of various diseases. Insufficient kidney Yang may cause somnolence, cold, cold limbs, impotence, infertility in women, diarrhea, along with other symptoms of Yang deficiency (Wu et al., 2022). In the pharmacological research on kidney-yang deficiency, glucocorticoids, sex hormones, adenine, and kidney/bilateral adrenalectomy are often used to establish animal models of kidney-yang deficiency. Conventionally, they are supplemented with large doses of exogenous glucocorticoids (hydrocortisone) to further establish the model (Sun et al., 2017).
In previous studies, our research mainly focused on the mechanism of “heating” and “suet oil” processing factors from the perspective of the absorption and metabolism of flavonoids in Epimedium, indicating that there were differences in the absorption and metabolism of flavonoid glycosides in Epimedium (Chen et al., 2009). In the processing process, “heating” could change the content of main active flavonoids of Epimedium, hence produced more easily absorbed bioactive flavonoids (Chen et al., 2007; Sun et al., 2014). The excipient “suet oil” could further promote the formation of self-assembled micelles of Epimedium flavonoids in vivo, increased the solubility of active flavonoids, and improved the absorption of active flavonoids, to achieve the purpose of synergism (Gu et al., 2019). However, the effect of traditional Chinese medicine on human body is a biological process of interaction and integration between “intervention system (TCM) and response system (biological organism)" (Xiang et al., 2012). It is difficult to match corresponding disease targets through the disease database for TCM diseases, especially for deficiency diseases. Therefore, it is necessary to adopt new methods and strategies to systematically and comprehensively explore the processing mechanism of TCM.
Metabolomics can analyze all metabolites in the whole sample and comprehensively monitor the changes of metabolites in the body. Applying the methods and ideas of metabolomics to the processing of TCM can reveal the changes of metabolism pathway and action network of TCM before and after processing as a whole (Kui et al., 2022). Network pharmacology is a new method to systematically and integrally study drugs from the perspective of biological network stability. It can systematically study the molecular relationship between TCM components and complex diseases (Xie et al., 2019). Metabolomics can be used to construct the target database of disease metabolites, hence the combined application of metabolomics and network pharmacology is gradually increasing (Chen et al., 2019). The multi-omics technology in systems biology may elucidate the processing mechanism of TCM, from the active components to the target genes, proteins, metabolites and metabolic pathways of the body. Therefore, the integration of multi-omics technology is vital to clarify the mechanism of action in TCM and provide insights into the processing mechanism of Epimedium.
Epimedium treatment may cause the linkage effect of “active component group, target in vivo and biological metabolism”, so the changes of downstream biological metabolites directly reflect the status of upstream targets. Therefore, based on previous research, this study used metabonomics—Network pharmacology—Real-time quantitative polymerase chain reaction (RT qPCR) integration technology to obtain reasonable and objective data from biomarkers, metabolic pathways, active ingredient groups and potential targets, and build a “component—target—metabolism” network, thereby revealing the processing and synergistic mechanism of Epimedium fried with suet oil for warming the kidney and promoting yang on the whole. This study provides new insights for studying processing mechanism of traditional Chinese medicine.
2 Materials and methods2.1 Reagents and materials
Epimedium was purchased from Shanxi, identified as Epimedium brevicornum Maxim., courtesy of Professor Baolin Guo, Institute of Pharmaceutical Botany, Chinese Academy of Medical Sciences (batch number: 161231). Hydrocortisone injections were bought from Tianjin King York Group Co., Ltd. (batch number: 1407271; Tianjin, China). Suet oil was purchased from Inner Mongolia Xilin Gol Grassland Co., Ltd. (Inner Mongolia, China). Formalin solution was purchased from Shanghai Aladdin Bio-Chem Technology Co., LTD. (batch number: F140859, Shanghai, China). Hydrochloric Acid was purchased from Nanjing Aojia Chemical Co., Ltd. (batch number: 2013050203, Nanjing, China). 2-Chloro-L-phenylalanine (content ≥98%) was purchased from Shanghai Hengbo Biological Technology Co., Ltd. (Shanghai, China). Bis(trimethylsilyl) trifluoroacetamide (99% BSTFA +1% TMCS) was provided by Regis Technologies, Inc. (Morton Grove, United States). Sodium azide was bought from Shandong Xiya Chemical Industry Co., Ltd. (batch number: S7470, Jinan, China). Urease was bought from Sigma-Aldrich (batch number: SLBB0100V, St. Louis, MO, United States). Pyridine was purchased from Adamas Reagent, Ltd. (batch number: 030928, Shanghai, China). Chromatography grade methanol was purchased from TEDIA (batch number: 1603253, State of Ohio, United States). Ethanol (analytically pure) was purchased from Nanjing Chemical Reagent Co., Ltd. (batch number: 1603233, Nanjing, China). Trichloromethane (batch number: 10006818) and propan-2-ol (batch number: 80109218) were purchased from Sinopharm Chemical Reagent Co., Ltd. (Shanghai, China). HyPure TMMolecular Biology Grade Water was provided by HyClone (batch number: SH30538.02, Logan, United States).
Agilent 7890B gas chromatography system (Agilent, United States) coupled with a LECO Chroma TOF Pegasus HT mass spectrometer (LECO, United States) were employed for metabonomics analysis. Ultrapure water was prepared by a Milli-Q water purifier (Merck Millipore, Germany). Heraeus Fresco 17 Centrifuge (Thermo Scientific, United States). BPG-9050AH high temperature blast drying oven (HASUC, China). Vortex Genie® 2 Vortex instrument (Scientific Industries, United States). CX23 electron microscope (Olympus, Japan). TNG-T98 Frozen Concentration Centrifugal Dryer (Taicang Huamei, China). ELISA (Spectra max 190, Shanghai Meigu Molecular Instrument Co., Ltd., shanghai, china). Quantitative Real-time PCR LightCycler instrument was provided by Bio-Rad Laboratories, Inc. (California, United States). Servicebio® RT First Strand cDNA Synthesis Kit and 2×SYBR Green qPCR Master Mix (None ROX) were obtained from Wuhan servicebio technology Co., LTD. (Wuhan, China).
2.2 Preparation of Epimedium processing and extraction
The prepared Epimedium decoction pieces (1 kg) were selected and made three copies in parallel. One sample was evenly added with hot melted suet oil (200 g), and mixed well. Then, it was oven-baked at 170°C for 7 min to get the Epimedium fried with suet oil. The second sample was placed into the oven without oil to obtain a similarly heated Epimedium. The third decoction piece was reserved as raw control.
All three samples were immersed in boiling water (1:30) for 1 h before extraction. Then it was extracted three times, 30 min for each time, filtered through gauze, respectively. The next three filtrates were merged and concentrated under vacuum, thus extracting solution was obtained respectively. The extraction flow chart is shown in Supplementary Figure S1.
2.3 Animals and treatments
Thirty-six adult male Sprague-Dawley rats (body weight 230 ± 20 g, animal licence No. SCXK(HU)2014-0001) were provided by Animal Experiment Center of Soochow University (Suzhou, China), and raised in the Experimental Animal Center of Jiangsu Academy of Traditional Chinese Medicine (Nanjing, China). Rats were housed in a climate-controlled room (relative humidity of 45% and ambient temperature 22°C–24°C), kept on a 12 h/12 h light-dark cycle, with free access to food and drinking water.
After a 2-week habituation, all rats were randomly divided into 6 groups with 6 rats in each; Control group (N, No. N1∼N6), Model group (M, No. M1∼M6), Suet oil group (A, No. A1∼A6), Epimedium raw product group (B, No. B1∼B6), Epimedium heating product group (C, No. C1∼C6), Epimedium fried with suet oil group (D, No. D1∼D6). M, A, B, C, and D group rats were injected intraperitoneally with hydrocortisone once daily for the following 10 days. The injection dose was 15 mg kg−1 on the first day, the dose was halved to 7.5 mg kg−1 on the second day and subsequent 9 days. N group was injected intraperitoneally with an equal volume of 0.9% sodium chloride injection. A, B, C, and D group rats were administered via intragastric with suet oil emulsion (0.144 g kg−1, 20% of Epimedium content, dispersed in water with CMC-Na), the Epimedium raw product extract (0.72 g kg−1), the Epimedium heating product extract (0.72 g kg−1), and the Epimedium fried with suet oil extract (0.72 g kg−1), respectively. N and M groups were orally administered with the same volume of pure water. The intragastric administration was performed once a day for 10 days. Meanwhile, body weights of the rats in each group were recorded daily. All animal treatments and experiments were approved by the Animal Ethics Committee of Jiangsu Provincial Academy of Chinese Medicine, and were strictly performed in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals (AEWC-20200702-119).
2.4 Collection of samples
After the 10th day of administration, all rats were fasted but had free access to water for 12 h while urine samples were collected. The 12 h urines of all rats were collected into a 50 mL centrifuge tube containing 20 µL of 1% sodium azide (NaN3). After centrifugation at 4,500 rpm min−1 for 10 min in a 4°C low-temperature centrifuge, the supernatant was taken and stored at −80°C until metabonomics analysis. Afterwards, blood samples were drawn from the orbit and collected in two 2 mL centrifuge tubes, one of which contained 20 µL heparin sodium solution for preparing plasma samples, the other for preparing serum samples. Then, centrifugation was performed at 4,000 rpm min−1 for 10 min in a 4°C low-temperature centrifuge, and the supernatant was taken to obtain plasma samples and serum samples. The plasma and serum samples were stored at −80°C until plasma metabonomics analysis and serum biochemical analysis.
Biochemical indicators included the hypothalamus-pituitary-adrenal axis in kidney-yang deficiency rats (Chen and Wang, 2015; Liu et al., 2016; Zhang et al., 2020): CORT, ACTH and 17-OHCS; Hypothalamus-pituitary-thyroid axis related indicators: T3, T4, and TSH; Hypothalamus-pituitary-gonad axis related indicators: LH, T, and FSH.
Thereafter, all rats were euthanized in parallel. The kidney tissues were quickly excised, directly frozen in liquid nitrogen, and stored at −80°C for subsequent quantitative real-time PCR assays. Simultaneously, the hypothalamus, pituitary, thyroid, adrenal gland, testis and kidney of all rats were collected and fixed in formalin solution for histopathological analysis.
2.5 Preparation of metabonomic samples
The urine samples were thawed before analysis and 100 µL aliquots of urine sample were added to 20 µL urease (80 mg mL−1). Then the mixture was shaken and mixed for 30 s, and incubated in an oven at 37°C for 1 h. Subsequently, urine sample mixture and 100 µL aliquots of plasma sample were mixed with 350 µL methanol and 20 µL L-2-chlorophenylalanine. Post vortex-mixing for 30 s and centrifugation at 16000 rpm for 20 min at 4°C, the supernatant of urine sample and plasma sample were transferred into a 2 mL gas sample bottle (siliconized methane). Moreover, 11 µL of each urine sample and plasma sample were accurately pipetted and mixed as the urine and plasma quality control (QC) sample, respectively. After the supernatants of the urine sample and plasma samples were completely dry in a vacuum concentrator, 60 µL of methoxyamine salt reagent (methoxamine hydrochloride, dissolved in 20 mg mL−1 pyridine) was added and mixed well, and incubated in an oven at 80°C for 30 min. Thereafter, 80 µL of BSTFA (containing 1% TCMS, v/v) were quickly added to each sample. The mixture was incubated in an oven at 70°C for 2 h and cooled to 25°C. Finally, 10 µL saturated fatty acid methyl ester standard mixture was added and mixed well for GC-MS analysis.
2.6 Metabonomics analysis
The analysis of plasma and urine samples was performed on an Agilent 7890B gas chromatography system (Agilent, United States) coupled with a LECO Chroma TOF Pegasus HT mass spectrometer (LECO, United States), which was equipped with Agilent DB-5MS capillary column (30 m × 0.25 mm × 0.25 μm, J&W Scientific, Folsom, CA, United States). GC-TOF-MS analysis conditions were set as follows: the injection volume was 1 µL with non-split mode. Carrier gas was helium used with a column flow rate of 1.0 mL min−1, the column temperature maintained at 50°C for 1 min and then heated at a rate of 20°C·min−1 up to 310°C maintained for 6 min. The forward inlet purge flow rate was 3 mL min−1 and the forward inlet temperature was 280°C. Transmission line temperature and ion source temperature were 270°C and 220°C respectively. The ionization voltage was −70 eV. The scan range and speed were 50–500 m/z, 20 spectra/sec. The solvent delays of urine and plasma samples were 455 S and 366 S, respectively.
Random urine samples (N2), plasma samples (N5) and the internal standard sample (L-2-chlorophenylalanine) were tested according to the above conditions. Thereafter, a comparative analysis of the GC-TOF-MS total ion characteristics (TIC) between the random urine sample (N2), plasma samples (N5) and the internal standard sample were analyzed in order to corroborate and validate our experimental methodologies as well as the instrument platforms. The 45 urine samples (including 9 urine QC samples) and plasma samples (including 9 plasma QC samples) were also tested according to the above conditions, and TIC of urine and plasma GC-TOF-MS were obtained for all samples. The retention time of the internal standard (L-2-chlorophenylalanine) in each TIC was extracted to investigate the stability of the system.
First, the raw data collected by the GC-TOFMS system were preprocessed. Single data and single peak in the original data were filtered, the missing value in the original data were recoded, and the internal standard (IS) to were utilized to normalize each group of data. Meanwhile, both TOF 4.3X software (LECO Corporation, United States) and LECO Fiehn Rtx5 database were used for raw peak exact matching, baseline filtering and calibration, peak alignment, deconvolution analysis, peak identification and integration of peak area. Finally, the standardized data were imported into SIMCA software (V14.1, MKS Data Analytics Solutions, Umea, Sweden) for multivariate statistical analysis, such as, principal component analysis (PCA) and orthogonal projection to latent structures-discriminant analysis (OPLS-DA) analysis. Prior to PCA, all variables obtained from GC-MS data sets went through logarithmic (LOG) transformation and centralization (CTR) scaling. Meanwhile, prior to OPLS-DA, all variables obtained from GC-MS data sets underwent logarithmic (LOG) transformation and unit variance (UV) scaling. Afterwards, the validity of the OPLS-DA model was tested by 7-fold cross validation and permutation test (random 200 times). Multivariate statistical analysis was conducted, combined with traditional univariate statistical analysis, with variable importance in the projection of the first principal component of the OPLS-DA model greater than 1 (VIP >1), while the p-value of Student’s t-test is less than 0.05 (p < 0.05), was used to screen significantly different metabolites, and the result was visualized in the form of a volcano graph. Thereafter, significantly different metabolites were subjected to pathway analysis with KEGG (https://www.kegg.jp/) and MetaboAnalyst 5.0 (https://www.metaboanalyst.ca/), which is a web-based tool for visualization of metabolomics, to identify related metabolic pathways. Metabolic pathways with p < 0.05 or Impact >0.10 were screened as the most relevant metabolic pathways and may be used to distinguish the kidney-yang deficiency model group from the normal control group.
2.7 Network pharmacology
The active ingredients of Epimedium fried with suet oil mainly include two parts: the active ingredients of Epimedium flavonoids and the chemical ingredients of “suet oil” as the processing auxiliary material. Based on the previous research of the research group (Jiang et al., 2014; Li et al., 2020), 32 flavonoids were identified in raw and processed Epimedium products and 25 fatty acids were identified in suet oil, and made available via PubChem database (https://pubchem.ncbi.nlm.nih.gov/) to be downloaded in SDF format for the 2D or 3D structure of each component. The components were then screened with a definite 2D or 3D structure as the component library of Epimedium fried with suet oil. Thereafter, based on the PharmMapper database (http://www.lilab-ecust.cn/pharmmapper/), the species was set to “Homo Sapiens”, and the potential targets of Epimedium pilaris were predicted. According to the z score value, the top 15 targets of each component corresponding to the score were screened as the action targets of chemical components in Epimedium. Targets were canonicalized to standard gene names via the UniProt database (https://www.uniprot.org). The targets of chemical constituents in Epimedium fried with suet oil are imported into the STRING database (https://www.string-db.org/) to construct a PPI network, and then the targets are imported into Cytoscape 3.6.1 to construct an active ingredient-potential target network. Finally, these targets were uploaded to the Metascape database (https://metascape.org/gp/index) for KEGG pathway analysis.
2.8 Integrated metabonomics and network pharmacology analysis
The metabolic pathways of p < 0.05 or Impact >0.1 obtained from the metabonomics analysis were mapped to pathways acquired by the network pharmacology to uncover pathways in which the two correlates. Afterwards, target proteins and active compounds corresponding to Epimedium fried with suet oil from network pharmacology were reversely searched from the overlapping pathways, and finally the potential active compounds, target proteins, and metabolic pathways of Epimedium fried with suet oil affecting KYDS were obtained. On this basis, these potential targets were further verified via RT-qPCR, and differential target proteins were shown.
2.9 Quantitative real-time PCR analysis
Firstly, total RNA was isolated and extracted from flash-frozen kidney tissues. 100 mg of tissue was taken and added into a pre-cooled homogenized tube (containing 1 mL RNA extract). Grind thoroughly with the grinder until no tissue is visible, and then take the supernatant after centrifuging 10 min at 12,000 rpm. Then 250 µL of trichloromethane was added, the centrifuge tube was inverted for 15 s, mixed well, and left for 3 min. After separating the cores at 12,000 rpm at 4°C for 10 min, 400 µL of supernatant was transfered to a new centrifugal tube, adding 0.8 times the volume of propan-2-ol and mixing it upside-down. After being placed at −20°C for 15 min and centrifuged at 12,000°C for another 10 min, the white precipitate at the bottom of the tube was RNA. The liquid was removed by suction, and 1.5 mL of 75% ethanol was added to wash the precipitate, and then centrifuged at 4°C at 12,000 rpm for 5 min. After removing the liquid by suction, the centrifuge tube was placed on the ultra-clean table and blown for 3 min. The RNA was dissolved by adding 15 µL of RNase-free water and incubated at 55°C for 5 min. The concentration and purity of total RNA samples were determined by Nanodrop 2000. Thereafter, RNA was reversely transcribed to cDNA by using Servicebio®RT Enzyme Mix. Finally, RT-qPCR and data collection were performed with PCR Master Mix (Bio-rad). The samples were exposed to pre-denaturation at 95°C for 10 min, followed by 40 cycles of denaturation at 95°C for 15 s and annealing at 60°C for 30 s. The dissolution curve conditions were from 65°C to 95°C, and the fluorescence signal was collected for every 0.5°C increase in temperature. Target mRNA expression in each sample was normalized to the housekeeping gene (GAPDH) to normalize the starting cDNA levels. The 2−ΔΔCT method was used to calculate relative mRNA expression levels. The PCR primers were listed in Supplementary Table S1. The experiment was repeated thrice.
2.10 Statistical analysis
The data from body weight and biochemical indicators were screened by SPSS 26.0 statistical software package (IBM, Armonk, NY, United States) and GraphPad Prism 6.01 (GraphPad Software, Inc., San Diego, CA, United States) software with Student’s t-test or one-way analysis of variance (ANOVA). Data were presented as the mean ± standard deviation. In all experiments, confidence level was set at 95% and 99% to determine the significance of difference (p < 0.05, p < 0.01).
3 Results3.1 Weight changes
The weight changes of rats in each group were shown in Figure 1, indicating gradual increase of the body weight of rats in the N group from day 1 until day 10 of the experiment. After modeling, the body weight of rats in the M group decreased to the lowest at day 3, and gradually recovered to the original body weight at day 10. The weight of rats in all treatment groups showed initial decrease followed by a gradual recovery, and the weight changes of rats in the D group was approximate to that of N group.
Line chart of changes in body weight of rats during the experiment N, normal group; M, model group; A, suet oil group; B, Epimedium raw product group; C, Epimedium heating product group; D, Epimedium fried with suet oil group.
3.2 Biochemical and histopathological results
Figure 2 showed that compared with the N group, the levels of CORT, ACTH, 17-OHCS, T3, T4, TSH, LH, T, and FSH in the M group were significantly decreased (p < 0.05, p < 0.01). Apart from ACTH, the A group had a callback effect on other indicators, although there was no significant difference. The B group has significant callback to T3, T and FSH (p < 0.05). The C group had significant callback to 17-OHCS, T3, TSH (p < 0.05, p < 0.01). All biochemical indexes were significantly increased in the D group.
Biochemical indexes in kidney-yang deficiency rats N, normal group; M, model group; A, suet oil group; B, Epimedium raw product group; C, Epimedium heating product group; D, Epimedium fried with suet oil group. Significance compared to N group, *p < 0.05, **p < 0.01; Significance compared to M group, #p < 0.05, ##p < 0.01.
As shown in Figure 3, in the kidney-yang deficiency group M, multiple mild hemorrhages were observed in the hypothalamus tissue. Pituitary tissue was slightly bleeding, the boundary between nucleus and cytoplasm was unclear. In the thyroid tissue, a large number of follicular matrix was reduced, many epithelial cells were necrotic and exfoliated, and nuclear fragmentation was observed. Hemorrhage was observed locally in adrenal tissue, a large number of parenchymal cells had granular degeneration, cells were swollen, and the cytoplasm was loose and lightly stained, showing fine granularity. A small amount of eosinophilic serous material exuded from multiple interstitial areas of testis tissue. In the renal medulla, the renal tubules were significantly dilated, accompanied by mild connective tissue hyperplasia, and punctate infiltration of lymphocytes and macrophages. These were consistent with the results of Tang et al. (Xiang et al., 2012). Except for group A, the hypothalamus-pituitary-target gland axis (hypothalamus, pituitary, thyroid, adrenal gland, testis) in group B/C/D was improved to a certain extent after administration, and the improvement effect of group D was better than that of group B/C. Taken together, these results suggested that the hypothalamic, pituitary, thyroid, adrenal, testis and kidney tissue and cell structures could be altered in the kidney-yang deficiency induced by intraperitoneal injection of hydrocortisone, again demonstrating that the hypothalamic-pituitary gland—the target gland axis was inhibited. After the intervention treatment of Epimedium processed by heating and suet oil, the pathological changes of tissue cell structure caused by kidney yang deficiency was shown to have improved, which further revealed the mechanism of Epimedium fried with suet oil to warm the kidney and promote yang.
H&E staining of hypothalamus, pituitary, thyroid, adrenal, testis and kidney tissue sections (magnification ×200). N, normal group; M, model group; A, suet oil group; B, Epimedium raw product group; C, Epimedium heating product group; D, Epimedium fried with suet oil group. Hypothalamus: Black arrow represents the bleeding. Pituitary: Yellow arrow represents the bleeding. Black arrow represents the cell nucleus and the cytosolic boundary. Thyroid: Black arrow represents the follicular matrix. Yellow arrow represents epithelial cell necrosis. Red arrow represents nuclear fragmentation. Adrenal: Black arrow represents parenchymal cell degeneration. Testis: Black arrow represents eosinophilic serous material. Kidney: Black arrow represents the renal tubules. Yellow arrow represents the connective tissue. Red arrow represents the cellular infiltration.
3.3 Metabonomics results3.3.1 Metabolic profile analysis of urine and plasma in rats
The internal standard (L-2-chlorophenylalanine) retention time, the GC-TOFMS TIC of the random sample (N2, N5) and internal standard (L-2-chlorophenylalanine) were shown in Supplementary Tables S2, S3, Supplementary Figures S2–S5. The results showed that the experimental method was viable; the instrument platform was justified, and corroborated with the conditions in the experimental methodologies. Six groups of total ion current diagrams of the urine samples and serum samples analyzed by GC-MS were shown in Supplementary Figures S6, S7. A total of 759 active peaks were detected in urine and 261 active peaks in plasma.
3.3.2 Multivariate statistical analysis and potential biomarkers exploring
An unsupervised PCA was performed first to characterize between-group differences in rat urine and plasma samples, respectively. As displayed in Figure 4A, B, the urine and plasma samples of the rats in the control group and the model group were clustered into two categories on the PCA score plots, indicating that after modeling, the urine and plasma endogenous metabolites of the rats in the two groups existed significant differences. As shown in Figures 4C, D, the rat urine and plasma QC samples were well-aggregated, indicating that the detection system was stable and reliable, and the obtained differences could reflect the biological differences between samples. The urine and plasma samples of the A group, the B group, the C group, and the D group were all separated from the M group in the PCA score plots and approached the N group, indicating that each administration group could regulate the metabolic disorder of hydrocortisone-induced kidney-yang deficiency rats to normal to different degrees, and the regulation trend was as follows: D group > C group > B group > A group.
PCA scores plots of urine and plasma samples in rats N: Control group; M: Model group; A: Suet oil group; B: Epimedium raw product group; C: Epimedium heating product group; D: Epimedium fried with suet oil group. (N VS. M, (A) Urine, (B) Plasma; the PCA scores plots of six group, (C) Urine, (D) Plasma.)
To further determine the differences of endogenous metabolites in the urine and plasma of rats after administration, OPLS-DA analysis was performed on the M group and the N group, respectively. As shown in Figure 5A, B, both in the urine samples and in the plasma samples, the N group and the M group were respectively clustered into one category in the OPLS-DA scores plots, indicating a significant difference between the N group and the M group. Then, a random permutation test (n = 200) was performed under the established OPLS-DA model to evaluate the model’s interpretation rate (R2Y), predictive ability (Q2) and other parameters. The results were arranged in the experiment. R2 = 0.900, Q2 = −0.330 (urine, Figure 5C); R2 = 0.850, Q2 = −0.430 (plasma, Figure 5D). The results show that the model has good stability and prediction ability, and no overfitting occurs. Combined with VIP >1 and p < 0.05 in the Volcano plot (Figures 5E, F), the differential variables with the largest contribution were screened as biomarkers related to kidney-yang deficiency syndrome, and were identified by LECO Fiehn Rtx5 database and combined with HMDB databases. Twenty-eight relevant biomarkers (13 in urine and 15 in plasma) were finally screened, of which 15 were identified (5 in urine and 10 in plasma). The results are listed in Table 1. Compared with the N group, nine biomarkers were upregulated in the M group, four of which were Tyramine, Galacturonic acid, Gulonate, 3-Hydroxyproline in urine. And 5 in plasma were D-Mannose, Ornithine, Sucrose, 2,4,6-trimethylPyridine, Cholesterol. Six biomarkers in the M group were downregulated, 1 in urine was downregulated as Ribitol, and five in plasma were downregulated as 5-Oxoproline, Citric acid, L-Tryptophan, Glycerol, Threonic acid.
Screening of biomarkers in rats with kidney-yang deficiency. N VS. M, (A) Urine OPLS-DA scores plots, (B) Plasma OPLS-DA scores plots, (C) Urine OPLS-DA permutation test, (D) Plasma OPLS-DA permutation test, (E) Urine Volcano plot, (F) Plasma Volcano plot. Red represents upregulation of differential metabolites. Blue represents downregulation of differential metabolites. Gray represents metabolites with no statistical difference.
Biomarkers of urine and plasma in kidney-yang deficiency rats (n = 6, Mean ± SD).
Sample
Metabolite
tR/min
Formulas
Similarity
Datebase ID
VIP
P
Variation between groups
M vs. N
A vs. M
B vs. M
C vs. M
D vs. M
Urine
Ribitol
17.776
C5H12O5
864
HMDB0000508
2.103
0.022
↓*
↑
↑
↑
↑
Urine
Tyramine
22.732
C8H11NO
920
HMDB0000306
1.990
0.032
↑*
↑
↓
↓
↓
Urine
Galacturonic acid
23.108
C6H10O7
791
HMDB0002545
1.250
0.009
↑**
↓
↓
↓
↓
Urine
Gulonate
29.778
C6H12O7
909
HMDB0003290
2.468
0.022
↑*
↓
↓
↓
↓
Urine
3-Hydroxyproline
26.867
C5H9NO3
986
HMDB0245903
1.478
0.037
↑*
↓
↓
↓
↓
Plasma
5-Oxoproline
10.037
C5H7NO3
954
HMDB0000267
2.122
0.042
↓*
↑
↑
↑
↑
Plasma
D-Mannose
11.970
C6H12O6
942
HMDB0000169
2.097
0.046
↑*
↑
↓
↓
↓
Plasma
Ornithine
11.657
C5H12N2O2
935
HMDB0000214
2.202
0.027
↑*
↓
↓
↓
↓
Plasma
Citric acid
11.616
C6H8O7
916
HMDB0000094
1.918
0.021
↓*
↓
↑
↑
↑
Plasma
L-Tryptophan
13.690
C11H12N2O2
762
HMDB0000929
2.185
0.008
↓**
↑
↑
↑
↑
Plasma
Sucrose
15.340
C12H22O11
733
HMDB0000258
2.768
0.000
↑**
↓
↓
↓
↓
Plasma
2,4,6-trimethylPyridine
6.477
C8H11N
895
HMDB0245482
1.642
0.012
↑*
↓
↓
↓
↓
Plasma
Glycerol
8.229
C3H8O3
779
HMDB0000131
2.018
0.020
↓*
↑
↑
↑
↑
Plasma
Threonic acid
10.115
C4H8O5
745
HMDB0000943
2.184
0.010
↓*
↓
↑
↑
↑
Plasma
Cholesterol
19.184
C27H46O
872
HMDB0000067
2.427
0.002
↑**
↓
↓
↓
↓
‘↑’ and ‘↓’ represent compounds which are up- and downregulated in the M group compared with the N group or in the administration group (A/B/C/D) compared with the M group. *p < 0.05, **p < 0.01 (one-way ANOVA, with a Bonferroni correction). Similarity: Score the match between the substance and the peak detected by mass spectrometry.
3.3.3 Regulation of differential metabolites
Column charts were drawn according to the relative contents and change trends of the biomarkers in urine and plasma of rats in each group (Figure 6). After the treatment intervention, the B, C, and D group biomarker contents demonstrated certain normalizing trends, with that of the D group being the most significant. The Sucrose and Galacturonic acid could be significantly recalled in the A group (p < 0.05, p < 0.01), suggesting that the suet oil had some intervention effect on kidney-yang deficiency syndrome. The B group could significantly recall six differential metabolites, especially Ornithine and Sucrose (p < 0.01). The C group could significantly recall seven differential metabolites, especially Sucrose and 2,4,6-trimethyl-Pyridine (p < 0.01). The D group could significantly recall 13 differential metabolites, especially Tyramine, Ornithine, Sucrose, 2,4,6-trimethyl-Pyridine and Glycerol (p < 0.01). Compared with B and C groups, D group had the most differentmetabolites in callbacks, indicating that the effect of Epimedium fried with suet oil group was the best. It is possible that Epimedium fried with suet oil regulated kidney-yang deficiency in rats through the above-described metabolites, thereby enhancing the warming effect on the kidneys and strengthen yang.
The levels of potential biomarkers of kidney-yang deficiency in each group Compared with the N group, *p < 0.05, **p < 0.01; compared with the M group, #p < 0.05, ##p < 0.01.
3.3.4 Metabolic pathway analysis
The biomarkers in Table 1 were imported into MetaboAnalyst 5.0 (http://www.metaboanalyst.ca) for pathway analysis, and the metabolic pathways with p < 0.05 or Impact >0.1 were screened out as potential kidney-yang deficiency metabolic pathways. As shown in Figure 7, five related metabolic pathways were finally screened out, namely Galactose metabolism, Glutathione metabolism, Glycerolipid metabolism, Arginine and proline metabolism and Tryptophan metabolism.
Pathways analysis of potential biomarkers.
Analysis of potential biomarkers of kidney-yang deficiency and related metabolic pathways suggested that the processing excipient of suet oil might play a role in warming kidney and promoting Yang by regulating Galactose metabolism. Epimedium raw product and Epimedium heating product might improve kidney-yang deficiency through Galactose metabolism, Glutathione metabolism and Arginine and proline metabolism. Meanwhile, the metabolic pathways above were thel metabolic pathways of Epimedium fried with suet oil to improve the kidney yang deficiency syndrome in rats, which further explained that the processing mechanism of epimedium could enhance its effect of warming kidney and promoting Yang after heating with the auxiliary material suet oil.
3.4 Network pharmacological analysis
Thirty-seven compounds were screened out as the active compounds library of Epimedium fried with suet oil, including 13 from Epimedium and 24 from suet oil. The specific results are listed in Supplementary Table S4. Homo sapiens species were selected to draw the protein-protein interaction (PPI) network of 159 common targets (Supplementary Figure S8). Then in order to reflect the relationship between the active compounds of Epimedium and suet oil and target genes in Epimedium fried with suet oil, a drug-active compounds-target network was constructed, as shown in the Figure 8. KEGG enrichment analysis was performed on 159 potential targets of active components of Epimedium fried with suet oil using Metascape database, and 267 pathways were obtained (Supplementary Figure S9).
Pharmacology network of the “drug-active compounds-target” regulated by Epimedium fried with suet oil pink triangles indicate drugs, blue hexagons indicate active compounds, yellow diamonds indicate target genes.
3.5 Integrative analysis of the network pharmacology and metabolomics
The five key metabolic pathways in metabolomics were mapped to the pathways of network pharmacology for analysis as shown in the Figure 9. It was found that two metabolic pathways (glutathione metabolism, arginine and proline metabolism) overlapped with it, and the corresponding 13 active compounds (9 active compounds in Epimedium and 4 in suet oil) and 7 target genes were obtained by reverse mapping analysis. Thirteen active compounds were Quercetin, Epimedin A, Caohuoside B, Korepimedoside C, Hexandraside E, Epimedin B, Epimedoside E, Baohuoside I, 2″-O-Rhamnosylicariside II, Methyl tetradecanoate, Methyl (Z)-heptadec-10-enoate, Methyl (9E)-9-octadecenoate, Methyl (E)-nonadec-10-enoate. The seven target genes were acetaldehyde dehydrogenase 2 (ALDH2), arginase type II (ARG2), glutathione S-transferase A3 (GSTA3), glutathione S-transferase M1 (GSTM1), glutathione S-transferase M2 (GSTM2), hematopoietic prostaglandin D synthase (HPGDS), Nitric Oxide Synthase 2 (NOS2). The drug-active compounds-target genes-pathways network was constructed as shown in Figure 10.
Overlapping pathways between metabolomics and network pharmacology analysis.
The network of drug-active compounds-target genes-pathways.
3.6 Experiment validation
The RT-qPCR verification results were shown in Figure 11. Seven target genes were expressed in kidney tissue. Compared with the N group, the expression levels of two target genes (ARG2 and NOS2) in the M group were significantly increased (p < 0.01). The expression of five target genes was decreased, and which of three target genes (ALDH2, GSTA3, HPGDS) were significantly decreased (p < 0.01). Compared with the M group, the A, B, C and D groups had different degrees of callback effect on seven target genes. ALDH2, ARG2, GSTA3, HPGDS, NOS2 had a callback trend in the A group, but there was no significant difference. The expression of HPGDS target genes was significantly increased in the B group (p < 0.05). The expression of ARG2 and NOS2 were significantly decreased in the B and C groups (p < 0.05). The expression levels of GSTM1, GSTM2 and HPGDS were significantly increased (p < 0.05, p < 0.01), ARG2 and NOS2 expressions were significantly decreased (p < 0.01) in the D group.
Epimedium fried with suet oil regulates differential gene expression associated with kidney yang deficiency. All genes expression were examined by RT-qPCR and normalized to GAPDH expression. Data represent mean ± SD for at least three independent experiments. Compared with the N group, *p < 0.05, **p < 0.01; compared with the M group #p < 0.05, ##p < 0.01.
The callback results of the target genes in each administration group showed that the improvement effect of Epimedium fried with suet oil on kidney-yang deficiency was better than that of the other three administration groups. Notably, the five target genes with significant callback effects in Epimedium fried with suet oil were derived from the target genes corresponding to the active compounds of Epimedium (GSTM1, GSTM2 ARG2 and NOS2) and the target gene corresponding to the active compounds of suet oil (HPGDS). Furthermore, the synergistic mechanism of Epimedium fried with suet oil to warm the kidney and enhance yang was clarified.
3.7 Overall interactive network diagram of drugs, active compounds, targets, pathways, and metabolites
Figure 12 showed the overall interaction network diagram of drugs, active compounds, targets, pathways and metabolites. Integrating network pharmacology, metabolomics and RT-qPCR analysis, it was found that Epimedium fried with suet oil mainly affects 13 active compounds, seven targets, and two biomarkers in two metabolic pathways, thereby improving kidney-yang deficiency. Its mechanism of action may mainly involve amino acid metabolism and oxidative stress.
The overall interactive network diagram pink triangles indicate drugs, blue hexagons indicate active compounds, yellow diamonds indicate target genes, purple ellipses indicate pathways, green “V” indicates biomarks.
4 Discussion
Epimedium is a TCM for kidney-yang deficiency. The most common processing method is suet oil processing, which aims to enhance the kidney-warming and yang-boosting effects of Epimedium. Previous studies have shown that the model of kidney yang deficiency was established through glucocorticoid, and its mechanism was to simulate the pathological state when the hypothalamic-pituitary-target gland axis (adrenal gland, thyroid gland, gonad) was inhibited, thereby producing deficiency symptoms similar to “kidney yang deficiency” (Zhang et al., 2014). Kidney-yang deficiency is manifested in macroscopic characterization (e.g., weight loss), biochemical standards (e.g., decreased hormone levels of CORT, 17-OHCS, T3, T4, and T), and organ histopathology (e.g., pathological changes in the hypothalamus, pituitary gland, thyroid gland, adrenal gland, and testis). This study found that compared with the N group, the body weight of the M group rats was reduced, and the levels of hypothalamic-pituitary-adrenal axis hormones (CORT, ACTH, 17-OHCS), hypothalamic-pituitary-thyroid axis hormones (T3, T4, TSH) and hypothalamic-pituitary-gonadal axis hormones (LH, T, FSH) were inhibited. Pathological changes occurred in the tissues and cellular structures of hypothalamus, pituitary gland, thyroid gland, adrenal gland, testis and kidney. Compared with M group, in addition to the A group, B, C, and D groups could reverse this inhibition of different level and improve the pathological changes of the organization, and the reverse effect of Epimedium fried with suet oil was better than that of raw epimedium and heated epimedium, suggesting that the effect of warming the kidney and promoting yang of epimedium is enhanced after being heated and processed with suet oil.
Subsequently, in this study, we analyzed the processing synergistic mechanism of Epimedium fried with suet oil through metabonomics combined with network pharmacology. In the metabolomic results, 15 biomarkers related to kidney-yang deficiency were screened in urine and plasma, involving five metabolic pathways, namely Galactose metabolism, Glutathione metabolism, Glycerolipid metabolism, Arginine and proline metabolism, Tryptophan metabolism. Among them, Galactose metabolism was related to energy metabolism, while Glutathione metabolism, Glycerolipid metabolism, Arginine and proline metabolism, and Tryptophan metabolism were closely related to amino acid metabolism and oxidative stress. Galactose metabolism regulated by suet oil and Galactose metabolism, Glutathione metabolism, Arginine and proline metabolism regulated by Epimedium heating product were the metabolic pathways of Epimedium fried with suet oil group to improve kidney yang deficiency in rats. In the network pharmacology results, seven target genes (ALDH2, ARG2, GSTA3, HPGDS, GSTM1, GSTM2, and NOS2) and two linked metabolic pathways (glutathione metabolism, arginine and proline metabolism) were identified by the drug-active compounds-target genes-pathways network.
Glutathione comprises cysteine, glutamic acid and glycine. The main metabolic pathways of glutathione include glutathione producing glutathione disulfide (GSSG) under the action of glutathione peroxidase (GSHPx) and cross-linking complex under the action of glutathione transferase (GST). It plays a role in oxidative stress, participating in cell apoptosis and regulating signal transduction (Xu et al., 2014). Studies have found that genes related to kidney-yang deficiency in the aged are closely related to aging (Dong et al., 2013), and glutathione deficiency can induce oxidative stress in the aging process. Increasing the uptake of the glutathione, precursors glycine and cysteine, can restore the synthesis and concentration of glutathione and reduce oxidative stress and oxidative damage in the aging process (Sekhar et al., 2011). (Du et al., 2020) used LC-MS metabolomics method to identify 31 biomarkers related to kidney Yang deficiency from testicular tissue, and 26 of them could be significantly recovered after Gulingji intervention. It was speculated that they might improve kidney Yang deficiency syndrome by regulating glutathione metabolism, arginine and proline metabolism, glyoxylate and dicarboxylic acid metabolism and other pathways. In this study, it was found that Ornithine and 5-Oxoproline are both products of glutathione metabolism and biomarkers of kidney Yang deficiency. Epimedium fried with suet oil could play its role in warming kidney and helping Yang by regulating these two biomarkers of kidney Yang deficiency. GSTA3, GSTM1 and GSTM2 are transferases of glutathione metabolism, which can catalyze the cross-linking complex between glutathione and electron-philic complex to protect cells from free radical damage. HPGDS, also known as glutathione dependent prostaglandin D synthetase, is involved in the biosynthesis of prostaglandin in adult testis, thereby affecting the production of spermatogenesis and the normal function of testis (Körber and Goericke-Pesh, 2019). The active components of Epimedium fried with suet oil, Quercetin, Epimedin A, Caohuoside B and Korepimedoside C, mainly acted on glutathione sulfur transferases (GSTA3, GSTM1 and GSTM2). Methyl tetradecanoate, Methyl (Z)-heptadec-10-enoate, Methyl (9E)-9-octadecenoate and Methyl (E)-nonadec-10-enoate mainly acted on HPGDS, and then affected glutathione metabolism.
Arginine and proline metabolism is the core process of the biosynthesis and metabolism of arginine, ornithine, proline, citrulline and glutamate (Feng et al., 2022). Arginine is a semi-essential amino acid synthesized from glutamine, glutamate, and proline via the gut-kidney axis in humans and most mammals, and plays an important role in the human body as a precursor of various physiologically important substances (Ma et al., 2019). Several studies have shown that arginine and its derivatives may be involved in the pathogenesis of kidney diseases and are markers of pathological progression (Feng et al., 2020). Arginine is the precursor of the synthesis of ornithine, which participates in the ornithine cycle under the action of ornithine carbamoyl transferase, promotes the synthesis of urea and accelerates the degradation of blood ammonia. Arginine can form ornithine and urea under the action of arginase, and produce nitric oxide under the action of nitric oxide synthase (Castilho-Martins et al., 2015). These pathways are involved in biochemical processes such as ammonia detoxification, hormone secretion and immune system regulation (Chang et al., 2020). The level of ornithine increased in kidney-yang deficiency rats. It was speculated that the regulation of ornithine in Epimedium fried with suet oil group might be related to the adjustment of gut microbiota balance, the reduction of renal tubulointerstitial fibrosis and the improvement of renal function. In this study, Hexandraside E, the active component of Epimedium fried with suet oil, mainly acted on ALDH2. Epimedin B and Epimedoside E mainly acted on ARG2. Baohuoside I and 2″-O-Rhamnosylicariside II mainly acted on NOS2. Thus, they jointly regulated arginine and proline metabolism.
Galactose metabolism is closely related to glycolytic metabolic pathways, which can reflect the level of energy metabolism of the body and is important for stabilizing the structure of the whole branched carbohydrate chain (Guo et al., 2016; Jin et al., 2019). Galactose metabolism is not only important for energy production in cells, but also contributes to the modification of glycolipids and glycoproteins (Petry and Reichardt, 1998). Galactose metabolism may be associated with chronic inflammation. When galactose metabolism is disrupted, it also promotes the release some inflammatory factors to participate in the body’s inflammatory response (Nikolac Perkovic et al., 2014). It is found that the liver and kidney of KYDS rats may be damaged due to the effect of hydrocortisone, thus affecting galactose metabolism (Chen et al., 2018). In this paper, the rats with kidney yang deficiency have cold phobia, kidney damage and other deficiency symptoms, indicating that their energy metabolism is abnormal. After administration, the abnormal galactose metabolism in the body gradually returned to normal, the body’s energy metabolism also tended to return to normal, and the inflammation was improved.
Glycerolipid metabolism is one of the important ways of lipid metabolism in the body, and the core reaction is the acetylation of glycerol and the degradation of triacylglycerol (Li et al., 2022). Lipids are divided into eight groups, including fatty acids, triglycerides, glycerolipids, sphingolipids, glycolipids, polyketones, sterolipids, and isopentenol lipids (Fahy et al., 2005). Changes in lipid metabolites can lead to a series of pathophysiological phenomena, including obesity (Nam et al., 2018), inflammation etc, (Ralston et al., 2017). Disturbance of Glycerolipid metabolism was observed in the early, middle and late stages of Aristolochic acid nephropathy (Zhao et al., 2015). Glycerol, as an important small molecule, participates in the Glycerolipid metabolism (Liu et al., 2019). The results of this experiment showed that glycerol, a differential marker screened, had an effect on Glycerolipid metabolism. Compared with the Control group, glycerol content in the kidney-yang deficiency model group was significantly decreased. Different administration groups could adjust the contents to different degrees, especially the callback effect of the D group had the most significant.
Tryptophan, as an essential amino acid, is involved in tryptophan metabolism (Liu et al., 2022). Tryptophan metabolism is closely related to the central nervous system because tryptophan is involved in serotonin synthesis (Kałużna-Czaplińska et al., 2019). Tryptophan metabolism is related to the activity of inflammatory bowel disease and has significant effects on gastrointestinal physiology (Hua et al., 2019). Zhou et al. found that tryptophan metabolism was seriously disturbed in KYDS state in the study of the therapeutic effect of ShenQiWan on kidney yang deficiency (Zhou et al., 2016). Tong et al. found that tryptophan is a precursor of serotonin, which is closely related to KYDS metabolic disorders requiring the participation of intestinal flora (Tong et al., 2022). In this study, tryptophan metabolism was disturbed in the model group of kidney-yang deficiency rats. This is consistent with the results of the literature. Compared with the control group, the tryptophan level in the model group decreased significantly. Epimedium fried with suet oil could improve the disorder of tryptophan metabolism by significantly regulating the tryptophan metabolite, showing obvious anti KYDS effect.
To further verify the results, RT-qPCR was used to analyze the expression levels of the seven selected targets in the kidney tissues of the six groups of rats. GSTA3, GSTM1, GSTM2, and HPGDS belong to the glutathione sulfur transferase series, which are related to prostaglandin synthesis and participates in glutathione metabolism. GSTA3, GSTM1, and GSTM2 are widely present in various tissues and cells of the body, and are related to the occurrence of cell damage, oxidative stress, poisoning, aging and other disease processes (Chen et al., 2018). Glutathione thitransferase A3 (GSTA3), as one of the most important members of the glutathione transferase family, is involved in detoxication and cell protection (Ilic et al., 2010). GSTA3 has a protective effect on Tubular epithelial-mesenchymal transition in renal fibrosis (Xiao et al., 2016). Glutathione S transferase M1 (GSTM1) is an important member of phase II toxic metabolizing enzyme GST family, which plays a detoxification role by binding with glutathione (Mehdi et al., 2019). It is reported that GSTM1 can repair acute kidney injury by reducing endoplasmic reticulum and oxidative stress (Dai et al., 2021). Glutathione S-transferase M2 (GSTM2) is a protein involved in the detoxification of reactive oxygen species (Zhou et al., 2008). GSTM2 has the activity of prostaglandin E synthetase (PGES) and can participate in the synthesis of prostaglandin E2 (PGE2) (Beuckmann et al., 2000).
HPGDS is a member of the glutathione S-transferases (GSTs) family that is closely associated with inflammatory diseases (Huang et al., 2021). HPGDS is mainly expressed in hematopoietic cell lines, catalyzes the synthesis of prostaglandin D2, and reduces inflammation by regulating the apoptosis of T cells and B cells (Rittchen and Heinemann, 2019). It was found that icariin can increase the levels of glutathione sulfotransferase and glutathione, and decrease the levels of malondialdehyde and nitric oxide (Amanat et al., 2022). In this study, the expression levels of GSTA3, GSTM1, GSTM2 and HPDGS were decreased in kidney tissues after modeling, suggesting that kidney Yang deficiency was accompanied by apoptosis, oxidative stress injury and inflammation of kidney tissue cells, and then caused kidney injury. This is consistent with the kidney histopathological results described above. However, D group could reverse the above level expression. Therefore, it could be speculated that Epimedium fried with suet oil might alleviate kidney injury caused by apoptosis of kidney tissue by increasing the level of glutathione transferase, and then played the role of warming kidney and helping Yang.
ALDH2 is related to the body’s energy metabolism and carbohydrate metabolism (Guo et al., 2020). ALDH2 is an important endogenous cardioprotective factor in mitochondria (Xia et al., 2023). Studies have shown that acetaldehyde dehydrogenase ALDH2, as a key enzyme in alcohol metabolism in human body, may affect the occurrence and development of coronary heart disease by affecting the level of blood glucose metabolism (You et al., 2018). ALDH2 plays an important regulatory role in the occurrence and development of tumors, and has a key role in maintaining tumor cell stemness (Toledo-Guzmán et al., 2019). ALDH2 is related to glycolysis and affects the level of blood glucose metabolism, which is consistent with the results of kidney yang deficiency involving Galactose metabolism pathway in this study. In this study, it was found that the expression level of ALDH2 in the kidney tissue of the M group was significantly decreased. It was speculated that the expression level of these related target genes was decreased during the occurrence of kidney Yang deficiency, resulting in the inhibition of the body’s energy metabolism, which was consistent with the literature that kidney Yang deficiency caused the body’s low energy metabolism (Chen et al., 2018). After administration, the callback trend of D group was closest to that of the N group. This suggested that Epimedium fried with suet oil could promote the expression of ALDH2, improve the body’s energy metabolism, and then improve the state of kidney Yang deficiency.
ARG2 is mainly expressed in kidney, brain and prostate, and located in mitochondria. In clinical studies, arginase activation is closely related to the occurrence and development of heart, lung and kidney ischemia-reperfusion injury, hypertension, erectile dysfunction, atherosclerosis, aging and other diseases (Wu et al., 2015). Studies have shown that cell aging caused by aging is associated with increased intracellular ARG2 activity/expression (Meng and Xiong, 2016). Inhibition of arginase I leads to increased expression of NOS2, thus promoting the production of NO (Raber et al., 2012). In this study, it was found that ARG2 expression level was significantly increased in the M group, suggesting that kidney Yang deficiency might be related to the increased expression level of ARG2 in the kidney tissue, the decreased content of NO and the increased level of superoxide production, which may lead to renal endothelial dysfunction. However, Epimedium fried with suet oil could significantly reduce the ARG2 expression levels, indicating that it might play its role by reducing ARG2 expression levels and alleviating the endothelial dysfunction caused by its overexpression.
NOS2 is the rate-limiting enzyme of NO synthesis, which plays its biological role by producing NO and has important clinical significance in the occurrence and development of diseases (Mount and Power, 2006). Excessive production of NO by NOS2 will lead to cell damage and tissue necrosis, and further promote the occurrence and development of inflammatory diseases. NOS2 is an important intracellular messenger and molecular marker in the mechanism of inflammation (Jiang et al., 2016). In the kidney, NO is involved in the regulation of renal vascular resistance, glomerular filtration rate and maintenance of renal structural integrity (Da Silva et al., 2021). Many studies have shown that NOS2 plays a key role in the pathogenesis of metabolic nephropathy, and it has been found in the pathogenesis of chronic nephropathy that elevated NOS2 levels produce highly reactive NO, which is associated with chronic low-grade inflammatory states (Zhang et al., 2021). In this study, the expression level of NOS2 in the D group was significantly decreased. It suggested that Epimedium fried with suet oil could reduce the production of NO, reduce the inflammatory state of kidney tissue and maintain the integrity of kidney structure by reducing the overexpression of NOS2 in the kidney tissue.
5 Conclusion
In this study, 15 differential metabolites in urine and plasma were screened by GC-MS technique. The regulatory effect of Epimedium on abnormal metabolism of kidney Yang deficiency was related with Galactose metabolism, Glutathione metabolism, Glycerolipid metabolism, Arginine and proline metabolism, Tryptophan metabolism. The processing excipient “suet oil” regulated Galactose metabolism, and “heating” regulated Galactose metabolism, Glutathione metabolism, Arginine and proline metabolism. Together, the two factors enhanced the role of Epimedium fried with suet oil in warming kidney and enriching yang. This further explained the processing mechanism of synergistic effect of Epimedium processed by heating and suet oil. Based on the results of metabolomics and network pharmacology, combined with RT-PCR verification, the results showed that the 13 active components of Epimedium fried with suet oil played the role of warming kidney and promoting Yang by acting on the seven target genes, two pathways and two biomarkers of kidney Yang deficiency, which further elucidated the processing and synergistic mechanism of Epimedium fried with suet oil. In view of the characteristics of multi-component, multi-target and overall synergy of TCM, this study comprehensively characterized the mechanism of processing effect of Epimedium, which would help to systematically and profoundly reveal the processing mechanism of Epimedium fried with suet oil. At the same time, it provides a research strategy for studying the processing mechanism of TCM and the action mechanism of TCM in treating deficiency syndrome.
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 Animal Ethics Committee of Jiangsu Provincial Academy of Chinese Medicine.
Author contributions
ES: Conceptualization, funding acquisition, design, supervision. RH: Data processing analysis and writing. KD: Formal analysis. LW: Experiments and database retrieval and analysis. JH: Experiments and writing. XT: Conceptualization, investigation. YW: Project administration. LF: Resources. XJ: Conceptualization, resources, funding acquisition.
Funding
This work was supported by the National Natural Science Foundation of China (No. 81503253, 81573833), the Jiangsu Provincial Medical Youth Talent (No. QNRC2016635), Natural Science Foundation of Nanjing University of Chinese Medicine (No. XZR2020028), Medical Innovation Team Project of Jiangsu Province (No. CXTDB2017003), National Administration of traditional Chinese medicine key research laboratory capacity improvement project (No. 2022GJJZDYJS-01), “Major New Drug Innovation” Science and Technology major project (No. 2017ZX09301056), and 2021 Postgraduate Research and Practice Innovation Project of Jiangsu Province (No. SJCX21-0715, SJCX21-0717).
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.2023.1113213/full#supplementary-material
ReferencesAmanatS.ShalB.Kyoung-SeoE.AliH.KhanS. (2022). Icariin attenuates cyclophosphamide-induced cystitis via down-regulation of NF-кB and up-regulation of Nrf-2/HO-1 signaling pathways in mice model. Int. Immunopharmacol.106, 108604. 10.1016/j.intimp.2022.108604BeuckmannC. T.FujimoriK.UradeY.HayaishiO. (2000). Identification of mu-class glutathione transferases M2-2 and M3-3 as cytosolic prostaglandin E synthases in the human brain. Neurochem. Res.25 (5), 733–738. 10.1023/a:1007579507804Castilho-MartinsE. A.CanutoG.MuxelS. M.daSilvaM.Floeter-WinterL. M.Del-AguilaC. (2015). Capillary electrophoresis reveals polyamine metabolism modulation in Leishmania (Leishmania) amazonensis wild-type and arginase-knockout mutants under arginine starvation. Electrophoresis36 (18), 2314–2323. 10.1002/elps.201500114ChangH.LiuQ.BaiW. F.BaiY. C.JiaX. Y.GaoC. (2020). Protective effects of Amygdalus mongolica on rats with renal fibrosis based on serum metabolomics. J. Ethnopharmacol.257, 112858. 10.1016/j.jep.2020.112858ChenR. Q.WangJ.LiaoC. B.ZhangL.GuoQ.WangX. F. (2018). Exploring the biomarkers and therapeutic mechanism of kidney-yang deficiency syndrome treated by You-gui pill using systems pharmacology and serum metabonomics. RSC Adv.8 (2), 1098–1115. 10.1039/c7ra12451aChenR. Q.WangJ.ZhanR. H.ZhangL.WangX. F. (2019). Integrated systems pharmacology, urinary metabonomics, and quantitative real-time PCR analysis to uncover targets and metabolic pathways of the you-gui pill in treating kidney-yang deficiency syndrome. Int. J. Mol. Sci.20 (15), 3655. 10.3390/ijms20153655ChenS.WangX. F. (2015). Progress on metabolomics studies of the Kidney Yang deficiency syndrome. Acad. J. Guangdong Coll. Pharm.31 (06), 829–832. 10.3969/j.issn.1006-8783.2015.06.030ChenY.JiaX. B.CaiY.DingA. W. (2007). Temperature and heating time - two key factors affecting icariin content. Chin. Tradit. Pat. Med.2007 (07), 1038–1041. 10.3969/j.issn.1001-1528.2007.07.030ChenY.JiaX. B.TanX. B.FanC. Y.HuM. (2009). Absorption and metabolism of flavonoids in Herba Epimedii via rat intestinal perfusion model. Chin. J. Chin. Mater Med.34 (22), 2928–2931.Chinese Pharmacopoeia Commission [CPC] (2020). Pharmacopoeia of people’s Republic of China. Beijing: Science and Technology Press.Da-SilvaC. B. P.CeronC. S.MendesA. S.de-MartinisB. S.CastroM. M.TirapelliC. R. (2021). Inducible nitric oxide synthase (iNOS) mediates ethanol-induced redox imbalance and upregulation of inflammatory cytokines in the kidney. Can. J. Physiol. Pharmacol.99 (10), 1016–1025. 10.1139/cjpp-2021-0108DaiT.YangP.LiuB.ZhangS. C. (2021). Effect of paeonol on the survival of ischemia and reperfusion induced cardiomyocytes by regulating miR-122-5p/GSTM1 molecular Axis. China Pharm.24 (02), 232–236. 10.3969/j.issn.1008-049X.2021.02.005DongY.DingY.LiuP. Z.SongH. Y.ZhaoY. P.LiM. (2013). Investigation of the material basis underlying the correlation between presbycusis and kidney deficiency in traditional Chinese medicine via GC/MS metabolomics. Evid. Based Complement. Altern. Med.2013, 762092. 10.1155/2013/762092DuK.GaoX. X.FengY.LiJ.WangH.LvS. L. (2020). Integrated adrenal and testicular metabolomics revealed the protective effects of Guilingji on the Kidney-Yang deficiency syndrome rats. J. Ethnopharmacol.255, 112734. pii: S0378-8741(19)34453-8. 10.1016/j.jep.2020.112734FahyE.SubramaniamS.BrownH. A.GlassC. K.MerrillA. H.MurphyR. C. (2005). A comprehensive classification system for lipids. J. Lipid Res.46 (5), 839–861. 10.1194/jlr.E400004-JLR200FengQ.LiY. M.YangY. W.FengJ. F. (2020). Urine metabolomics analysis in patients with normoalbuminuric diabetic kidney disease. Front. Physiol.11, 578799. 10.3389/fphys.2020.578799FengQ. S.TongL.LuQ.LiuS.ZhaoL. S.XiongZ. L. (2022). 1H NMR serum metabolomics and its endogenous network pharmacological analysis of Gushudan on kidney-yang-deficiency-syndrome rats. Anal. Biochem.643, 114580. 10.1016/j.ab.2022.114580GuH. M.SunE.LiJ.HouJ.JiaX. B. (2019). Effect of processing excipient suet oil on formation and absorption of baohuoside -bile salt self-assembled micelles. Chin. J. Chin. Mater Med.44 (23), 5143–5150. 10.19540/j.cnki.cjcmm.20190916.311GuoG. L.WangS.SunP.LiX. L.WangX.CaoG. S. (2020). Research progress on processing mechanism of epimedium fried with sheep fat oil based on warming kidney and promoting yang. J. Liaoning Univ. Tradit. Chin. Med.22 (07), 1–5. 10.13194/j.issn.1673-842x.2020.07.001GuoH.NiuX.GuY.LuC.XiaoC.YueK. (2016). Differential amino acid, carbohydrate and lipid metabolism perpetuations involved in a subtype of rheumatoid arthritis with Chinese medicine cold pattern. Int. J. Mol. Sci.17 (10), 1757. 10.3390/ijms17101757HuaY. L.MaQ.LiW.ZhangX. S.ChengX. H.JiaY. Q. (2019). Metabolomics analysis of Pulsatilla decoction on treatment of wetness-heat-induced diarrhea in rats based on UPLC-Q/TOF-MS/MS. Biomed. Chromatogr.33 (11), e4629. 10.1002/bmc.4629HuangC.GeF.RenW.ZhangY.WuX.ZhangQ. (2021). Copy number variation of the HPGDS gene in the Ashidan yak and its associations with growth traits. Gene772, 145382. 10.1016/j.gene.2020.145382IlicZ.CrawfordD.VakhariaD.EgnerP. A.SellS. (2010). Glutathione-S-transferase A3 knockout mice are sensitive to acute cytotoxic and genotoxic effects of aflatoxin B1. Toxicol. Appl. Pharmacol.242 (3), 241–246. 10.1016/j.taap.2009.10.008JiangJ.FengL.LiJ.SunE.DingS. M.JiaX. B. (2014). Multielemental composition of suet oil based on quantification by ultrawave/ICP-MS coupled with chemometric analysis. Molecules19 (4), 4452–4465. 10.3390/molecules19044452JiangQ. Y.ZhengM. S.YangX. J.SunX. S. (2016). Analysis of molecular networks and targets mining of Chinese herbal medicines on anti-aging. BMC Complement. Altern. Med.16 (1), 520. 10.1186/s12906-016-1513-2JinH.MaN.LiX.KangM.GuoM.SongL. (2019). Application of GC/MS-Based metabonomic profiling in studying the therapeutic effects of aconitum carmichaeli with ampelopsis japonica extract on collagen-induced arthritis in rats. Molecules24 (10), 1934. 10.3390/molecules24101934Kałużna-CzaplińskaJ.GątarekP.ChirumboloS.ChartrandM. S.BjørklundG. (2019). How important is tryptophan in human health?Crit. Rev. Food Sci. Nutr.59 (1), 72–88. 10.1080/10408398.2017.1357534KörberH.Goericke-PeschS. (2019). Expression of prostaglandin (PG) D synthase lipocalin and hematopoietic type and PG D receptor during restart of spermatogenesis following downregulation using a slow release GnRH agonist implant in the dog. Cell Tissue Res.378 (2), 359–370. 10.1007/s00441-019-03059-xKuiX. T.WangJ.ChenW.WeiY. P.LvY. H.ZuX. P. (2022). Metabolomics and its application of in Chinese medicine processing. Chin. J. Chin. Mater Med.47 (03), 593–602. 10.19540/j.cnki.cjcmm.20210927.601LiB.MaJ.ZhuL.YanY.WuY.ZhengC. (2022). Metabolomic profiling and alleviation of oxidative stress of Huangjing Zanyu capsule treating oligoasthenospermia. Transl. Androl. Urol.11 (5), 659–674. 10.21037/tau-22-293LiM. Y.SunE.XuF. J.XuJ. D.JiaoX. B. (2020). Analysis changes of Epimedii Folium's flavonoids before and after processing based on UPLC-Q/TOF-MS. Chin. Tradit. Herb. Drugs51 (11), 2900–2907. 10.7501/j.issn.0253-2670.2020.11.007LiuH.YangL.WanC.LiZ.YanG.HanY. (2022). Evaluation of the pharmacological effects and exploration of the mechanism of traditional Chinese medicine preparation Ciwujia tablets in treating insomnia based on ethology, energy metabolism, and urine metabolomic approaches. Front. Pharmacol.13, 1009668. 10.3389/fphar.2022.1009668LiuQ.ZhaoH. W.ZhangA. H.SunH.ZhouX. H.NanY. (2016). Chinmedomics strategy to discover effective constituents and elucidate action mechanism of Nanshi capsule against kidney-yang deficiency syndrome. China J. Chin. Mater Med.41 (15), 2901–2914. 10.4268/cjcmm20161526LiuZ.ZhangC.LeeS.KimW.KlevstigM.HarzandiA. M. (2019). Pyruvate kinase L/R is a regulator of lipid metabolism and mitochondrial function. Metab. Eng.52, 263–272. 10.1016/j.ymben.2019.01.001MaP.ZhangZ. M.ZhouX. Y.LuoJ. K.LuH. M.WangY. (2019). Characterizing semen abnormality male infertility using non-targeted blood plasma metabolomics. PLoS One14 (7), e0219179. 10.1371/journal.pone.0219179MehdiW. A.MehdeA. A.YusofF.RausR. A.ResenA. K.GhazaliH. (2019). The endothelial nitric oxide synthase gene G894T, glutathione S-transferase (GSTM1 and GSTT1) polymorphisms as a risk factor in the patient with nephrolithiasis. Int. J. Biol. Macromol.140, 719–726. 10.1016/j.ijbiomac.2019.08.184MengM.XiongX. D. (2016). Vascular endothelial cell senescence and vascular disease. Chin. J. Biochem. Mol. Biol.32 (03), 253–259. 10.13865/j.cnki.cjbmb.2016.03.03MountP. F.PowerD. A. (2006). Nitric oxide in the kidney: Functions and regulation of synthesis. Acta Physiol. (Oxf)187 (4), 433–446. 10.1111/j.1748-1716.2006.01582.xNamM.ChoiM. S.ChoiJ. Y.KimN.KimM. S.JungS. (2018). Effect of green tea on hepatic lipid metabolism in mice fed a high-fat diet. J. Nutr. Biochem.51, 1–7. 10.1016/j.jnutbio.2017.09.002Nikolac-PerkovicM.Pucic-BakovicM.KristicJ.NovokmetM.HuffmanJ. E.VitartV. (2014). The association between galactosylation of immunoglobulin g and body mass index. Prog. Neuropsychopharmacol. Biol. Psychiatry.48, 20–25. 10.1016/j.pnpbp.2013.08.014PetryK. G.ReichardtJ. K. (1998). The fundamental importance of human galactose metabolism: Lessons from genetics and biochemistry. Trends Genet.14 (3), 98–102. 10.1016/s0168-9525(97)01379-6RaberP.OchoaA. C.RodríguezP. C. (2012). Metabolism of L-arginine by myeloid-derived suppressor cells in cancer: mechanisms of T cell suppression and therapeutic perspectives. Immunol. Invest.41 (6-7), 614–634. 10.3109/08820139.2012.680634RalstonJ. C.LyonsC. L.KennedyE. B.KirwanA. M.RocheH. M. (2017). Fatty acids and NLRP3 inflammasome-mediated inflammation in metabolic tissues. Annu. Rev. Nutr.37, 77–102. 10.1146/annurev-nutr-071816-064836RittchenS.HeinemannA. (2019). Therapeutic potential of hematopoietic prostaglandin D2 synthase in allergic inflammation. Cells8 (6), 619. 10.3390/cells8060619SekharR. V.PatelS. G.GuthikondaA. P.ReidM.BalasubramanyamA.TaffetG. E. (2011). Deficient synthesis of glutathione underlies oxidative stress in aging and can be corrected by dietary cysteine and glycine supplementation. Am. J. Clin. Nutr.94 (3), 847–853. 10.3945/ajcn.110.003483SunE.WeiY. J.ZhangZ. H.CuiL.XuF. J.JiaX. B. (2014). Processing Mechanism of Epimedium fried with suet oil based on absorption and metabolism of flavonoids. Chin. J. Chin. Mater Med.39 (03), 383–390. 10.4268/cjcmm20140306SunX. T.LiuX.LiQ. (2017). Research overview of animal model with kidney deficiency syndrome. Chin. J. Tradit. Chin. Med. Pharm.32 (02), 659–662.Toledo-GuzmánM. E.HernándezM. I.Gómez-GallegosÁ. A.Ortiz-SánchezE. (2019). ALDH as a stem cell marker in solid tumors. Curr. Stem Cell Res. Ther.14 (5), 375–388. 10.2174/1574888X13666180810120012TongL.FengQ.LuQ.ZhangJ.XiongZ. (2022). Combined 1H NMR fecal metabolomics and 16S rRNA gene sequencing to reveal the protective effects of Gushudan on kidney-yang-deficiency-syndrome rats via gut-kidney axis. J. Pharm. Biomed. Anal.217, 114843. 10.1016/j.jpba.2022.114843WuQ.ZhangX. X.JiaoT. T.ZhangX. J.WangS. Y.LiuY. H. (2022). Correlations of characteristics of kidney yang deficiency syndrome with clinical indicators and risk factors of chronic kidney disease. Chin. J. Exp. Tradit. Med. Form.28 (19), 204–213. 10.13422/j.cnki.syfjx.20221726WuZ. S.YuY.LiuC.XiongY. Y.MontaniJ. P.YangZ. H. (2015). Role of p38 mitogen-activated protein kinase in vascular endothelial aging: Interaction with arginase-II and S6K1 signaling pathway. Aging (Albany NY)7 (1), 70–81. 10.18632/aging.100722XiaG.XuJ.ChenM.JinJ.WangX.YeY. (2023). Loss of ALDH2 aggravates mitochondrial biogenesis disorder in cardiac myocytes induced by TAC. Biochem. Biophys. Res. Commun.639, 189–196. 10.1016/j.bbrc.2022.12.001XiangZ.WangX. Q.LiuT.LvD. (2012). Thoughts and exploration on studying pharmacological mechanism of traditional Chinese medicine using network biology approach. Chin. J. Chin. Mater Med.37 (02), 146–151. 10.4268/cjcmm20120206XiaoY.LiuJ.PengY.XiongX.HuangL.YangH. (2016). GSTA3 attenuates renal interstitial fibrosis by inhibiting TGF-beta-induced tubular epithelial-mesenchymal transition and fibronectin expression. PLoS One11 (9), e0160855. 10.1371/journal.pone.0160855XieJ.GaoS.LiL.XuY. L.GaoS. M.YuC. Q. (2019). Research progress and application strategy on network pharmacology in Chinese Materia medica. Chin. Tradit. Herb. Drugs50 (10), 2257–2265. 10.7501/j.issn.0253-2670.2019.10.001XuS. L.LiW. H.JiaB.CuiJ. M.ShiN. G. (2014). Study on the differentially expressed genes of aged kidney-yang deficiency syndrome constitution and verify genes with RT-PCR. Chin. J. Tradit. Chin. Med. Pharm.29 (05), 1568–1572.YouL.LiC.ZhaoJ.WangD. W.CuiW. (2018). Associations of common variants at ALDH2 gene and the risk of stroke in patients with coronary artery diseases undergoing percutaneous coronary intervention. Med. Baltim.97 (19), e0711. 10.1097/MD.0000000000010711ZhangL. Z.LangJ. L.JinL.JinL. H.CaoB. B.ShaoX. (2020). Effect of hirsutella sinensis fungus on the hypothalamic-pituitary-adrenal Axis in lewis rats with kidney-yang deficiency syndrome. Evid. Based Complement. Altern. Med.2020, 5952612. 10.1155/2020/5952612ZhangZ. M.YangL.WanY.JiangS.ShangE. X.QianD. W. (2021). The synergic renoprotective effect of Rehmanniae Radix Preparata and Corni Fructus on adenine-induced chronic kidney disease rats based on integrated plasma metabolomics and network pharmacology approach. Life Sci.278, 119545. 10.1016/j.lfs.2021.119545ZhangZ. Y.ChenB. J.ZhangQ.WangY. R.SuY. X. (2014). Establishment method and evaluation index of animal model of kidney yang deficiency syndrome. Fujian J. Tradit. Chin. Med.45 (02), 61–63. 10.13260/j.cnki.jfjtcm.010691ZhaoY. Y.WangH. L.ChengX. L.WeiF.BaiX.LinR. C. (2015). Metabolomics analysis reveals the association between lipid abnormalities and oxidative stress, inflammation, fibrosis, and Nrf2 dysfunction in aristolochic acid-induced nephropathy. Sci. Rep.5, 12936. 10.1038/srep12936ZhouS. G.WangP.PiR. B.GaoJ.FuJ. J.FangJ. (2008). Reduced expression of GSTM2 and increased oxidative stress in spontaneously hypertensive rat. Mol. Cell Biochem.309 (1-2), 99–107. 10.1007/s11010-007-9647-7ZhouX. H.ZhangA. H.WangL.TanY. L.GuanY.HanY. (2016). Novel chinmedomics strategy for discovering effective constituents from ShenQiWan acting on ShenYangXu syndrome. Chin. J. Nat. Med.14 (8), 561–581. 10.1016/S1875-5364(16)30067-X