Integrated Metabolomics and Network Pharmacology to Establish the Action Mechanism of Qingrekasen Granule for Treating Nephrotic Syndrome

Nephrotic syndrome (NS) is a clinical syndrome resulting from abnormal glomerular permeability, mainly manifesting as edema and proteinuria. Qingrekasen granule (QRKSG), a Chinese Uyghur folk medicine, is a single-flavor preparation made from chicory (Cichorium intybus L.), widely used in treating dysuria and edema. Chicory, the main component in QRKSG, effectively treats edema and protects kidneys. However, the active components in QRKSG and its underlying mechanism for treating NS remain unclear. This study explored the specific mechanism and composition of QRKSG on an NS rat model using integrated metabolomics and network pharmacology. First, metabolomics explored the relevant metabolic pathways impacted by QRKSG in the treatment of NS. Secondly, network pharmacology further explored the possible metabolite targets. Afterward, a comprehensive network was constructed using the results from the network pharmacology and metabolomics analysis. Finally, the interactions between the active components and targets were predicted by molecular docking, and the differential expression levels of the target protein were verified by Western blotting. The metabolomics results showed “D-Glutamine and D-glutamate metabolism” and “Alanine, aspartate, and glutamate metabolism” as the main targeted metabolic pathways for treating NS in rats. AKT1, BCL2L1, CASP3, and MTOR were the core QRKSG targets in the treatment of NS. Molecular docking revealed that these core targets have a strong affinity for flavonoids, terpenoids, and phenolic acids. Moreover, the expression levels of p-PI3K, p-AKT1, p-mTOR, and CASP3 in the QRKSG group significantly decreased, while BCL2L1 increased compared to the model group. These findings established the underlying mechanism of QRKSG, such as promoting autophagy and anti-apoptosis through the expression of AKT1, CASP3, BCL2L1, and mTOR to protect podocytes and maintain renal tubular function.


INTRODUCTION
Nephrotic syndrome (NS) is a clinical condition caused by inflammation, oxidative stress, immune injury, and podocyte damage in the kidney (Imig and Ryan, 2013;Zhang et al., 2018). The disease is characterized by edema, proteinuria, hypoalbuminemia, and hyperlipidemia, due to abnormal glomerular permeability (Wang and Greenbaum, 2019;Tsuji et al., 2020). The current treatments for NS are extensive, repetitive, non-specific, and have side effects because of the understanding of the biological processes leading to NS (Saleem, 2019). Thus, NS is often accompanied by various complications, complicated treatment procedures, and unsatisfactory curative effects. Eventually, NS develops into renal failure, which negatively impacts livelihoods and the economy (Wang et al., 2000). Considering the far-reaching impacts of NS, it is essential to establish NS-related mechanisms and suitable drugs. Recently, many patients with kidney diseases have resorted to traditional Chinese medicines (TCM) (Zhao et al., 2015), whose biological activity and therapeutic effects are proven in vitro and animal experiments (Fang et al., 2013;Zhang et al., 2014). Chicory (Cichorium intybus L.), is a medicinal plant widely distributed in Europe, America, and Asia, and the main components of Qingrekasen Granule (QRKSG), listed in the Pharmacopoeia of the People's Republic of China (2020 edition). The plant functions as a diuretic, anti-inflammatory, digestive, cardiotonic, and liver tonic (Petrovic et al., 2004;Satmbekova et al., 2018). Moreover, the sesquiterpene lactones 8-Deoxylactucin, lactucin, and lactucopicrin in chicory extract inhibit the production of prostaglandin E2 (PGE2) and COX-2 protein expression (Cavin et al., 2005), hence, reducing incidences of edema and inflammatory pain (Wesolowska et al., 2006). For example, the chicory extract reduces carrageenan-induced paw edema (Rizvi et al., 2014). The extract significantly reduces α1microglobulin in urine, improves histological appearance, and repairs kidney damage (Pourfarjam et al., 2017). QRKSG is a granule made from chicory, and is approved by the China Food and Drug Administration, approval number (CFDA approval no. Z65020172). However, the functional mechanisms and value of QRKSG in NS are unknown.
Metabolomics and network pharmacology are two disciplines for understanding the mechanisms of drug action by linking the metabolites to disease targets. Metabolomics detect metabolic changes (Majumder et al., 2021) and the positive role of metabolites in physiology and disease development in organisms (Rinschen et al., 2019). In TCM, metabolomics revealed the mechanisms of various prescriptions, including the Qi-Dan Fang Granule (Wu et al., 2014) and the Danggui-Shaoyao-San (Wang et al., 2020). Meanwhile, network pharmacology predicts the possibility of drug compounddisease target combinations and pathway analysis (Xiong et al., 2019;Xiong et al., 2020). Thus, the approach is suitable for analyzing TCM preparations . For example, a new method of treating osteoarthritis has been established using network pharmacological analysis  thus, strengthening evidence-based research and promoting the development of high-quality Chinese medicine (Li and Zhang, 2013).
The findings in this research form a basis for clinical approval and QRKSG use for managing NS. The flow chart of the research process is shown in Figure 1.

In vivo Experiment
Thirty-two male SD (Sprague Dawley) rats (180-220 g) were obtained from Sibefu (Beijing) Biotechnology Co., Ltd. (YFY-2018-0,002, Beijing, China). The SD rats were housed in an environmentally conditioned room with a 12 h light/dark cycle, at 25 ± 2°C, and 50 ± 5% relative humidity. All rats accessed drinking water and feed ad libitum. The protocols in this study strictly followed the animal ethics standards of the Hubei University of Chinese Medicine. After 1 week of acclimation, the rats were randomly assigned into four groups, including the control (n 8) and ADR (7 mg/kg) treatments [three groups (n 24)].
The NS model was induced by injecting normal saline through the tail vein (Frenk et al., 1955;Korzets et al., 1997). A 24 h urine protein test was performed once a week for four experimental weeks. Mineral oil was added to the urine collection container to avoid evaporation. After successfully inducing the NS model, the ADR-treated rats were randomly assigned into three groups (n 8 per group): model, QRKSG (treated with a clinically equivalent dose of 1.62 g/kg/d), and benazepril (treated with a clinically equivalent dose of 0.9 mg/kg/d) groups. The control and model groups were treated with normal saline. All treatments were gavage administered once per day for 4 weeks, after which the animals were anesthetized, their kidneys were removed, and blood was collected from the abdominal aorta using a syringe.

Histopathological Analysis
One portion of the left kidney tissue from each rat was fixed in 4% paraformaldehyde impregnated wax, paraffin-embedded, and serially sectioned (4-7 μm). After staining, the kidney sections were imaged and photographed with a biomicroscope (Olympus BX53, Olympus Optical Co., Tokyo, Japan). Another portion of left kidney tissue (<1 mm 3 ) was fixed and post-fixed using 2.5% glutaraldehyde and 1% osmium tetroxide, respectively. The tissues were sequentially treated by ethanol acetone dehydration and acetone embedding infiltration overnight, embedded, sectioned (60-80 nm), uranium-lead double-stained, and finally observed by transmission electron microscopy (HT7700, Hitachi High-Technologies Corporation, Tokyo, Japan). Right kidney tissues were washed with saline, freeze-dried in liquid nitrogen on filter paper for 30 min, and stored in the refrigerator at −80°C.

Urine and Blood Chemistry
The collected blood samples were incubated at 4°C for 30 min, and the serum was separated by centrifugation (3,000 rpm for 10 min) for subsequent analysis. A biochemical analyzer (Chemray240, Rayto, Guangdong, China) detected albumin (ALB), total triglycerides (TG), creatinine (Cr), total cholesterol (TC), total protein (TP), and urea nitrogen (BUN) levels in the serum. In addition, the urine samples were centrifuged, and the collected supernatant was stored at −80°C for further analysis. The 24 h urine total protein concentration was determined using the BCA Protein Assay Kit and microplate reader (Epoch, Bio Tek, VT, United States) following the manufacturer's instructions.

Kidney Sample Preparation
The right kidney tissues were thawed at room temperature. Approximately 600 μL of the right kidney tissue homogenate was extracted with methanol, vortexed for 1 min, and ultrasonically centrifuged for 10 min in centrifuge tubes sitting in an ice bath to obtain the metabolites. Afterward, the supernatant was concentrated through centrifugation and evaporated to dryness. The internal standards: 10 μL D27myristic acid (0.75 mg/ml) and 40 μL methoxyamine pyridine solution (40 mg/ml) were then added and maintained for 90 min. Finally, 80 μL of MSTFA+1% TMCS (Sigma-Aldrich, MO, United States) and 20 μL n-hexane were added for derivatization and derivatization reaction termination. After centrifugation, 1 μL of the supernatant was transferred to a micro bottle for subsequent analysis. The GC-MS was performed as previously described (Cao et al., 2020), using Gas Chromatography-Mass Spectrometer (7,890B-5977B, Agilent Technologies, CA, United States) equipped with a DB-5MS quartz capillary column (0.25 mm × 30 mm, 0.25 μm). Analysis was performed using 1 µL injection volume with helium as the carrier gas (99.999%). The flow rate was 1.0 ml/min and a split ratio of 10:1. The column temperature was maintained at 60°C for 1 min, then increased at a rate of 10°C/ min up to 250°C and maintained at 250°C for 10 min. The injection temperature was 250°C, while that of the electron bombardment ion source was 230°C, with a solvent delay of 5.9 min and a scanning range of the mass spectrum at m/z 50.0-600.0.

Metabolites and Metabolic Pathway Analysis
The differential metabolites were identified by converting the GC-MS experimental data into compound retention time and relative peak area information using the Agilent Masshunter software. The online tool, MetaboAnalyst5.0 (Chong et al., 2018;Chong et al., 2019) (version 5.0), performed the principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA). The Pathway Analysis project in the online tool MetaboAnalyst5.0 analyzed metabolic pathways. Specified path analysis parameters included the visualization method using a scatter plot (testing significant features), the enrichment method using the hypergeometric test, the topology analysis approach using relative-betweenness centrality, and the reference metabolome technique using all compounds in the pathway library. Metabolites with p < 0.05 and VIP>1.0 were considered as potential biomarkers. The network analysis project in the MetaboAnalyst5.0 online tool performed the comprehensive metabolomic and network pharmacology analysis. Differential metabolites and differential metabolic genes were sequentially imported for analysis. Selected gene-metabolite interaction networks in the networks analysis options obtained the interaction between metabolites and targets.

UPLC-MS/MS Sample Preparation
A vacuum freeze dryer (Scientz-100F, Zhejiang, China) and grinder (MM 400, Retsch, Shanghai, China) were used to freeze-dry and crush the QRKSG. Afterward, 100 mg of the lyophilized powder was dissolved in methanol, vortexed six times, and refrigerated overnight at 4°C. The extracts were filtered for analysis on the UPLC-MS/MS system.

Network Pharmacology Analysis
First, the QRKSG biological targets were obtained from the SwissTargetPrediction database (http://www. swisstargetprediction.ch/). Secondly, the NS biological targets were obtained from GeneCards (https://www.genecards.org/), a searchable and integrated human gene database (Safran et al., 2010). Finally, the QRKSG target information for treating NS was obtained by interacting with the above two. The DAVID database (https://david.ncifcrf.gov/summary.jsp) was used as the functional annotation of gene ontologies (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways.

Comprehensive Analysis
The complex relationship of the metabolites-targets-pathways network was constructed using Cytoscape3.8.0 to determine the underlying mechanism of QRKSG on NS.

Molecular Docking
The AKT1, mTOR, CASP3, and BCL2L1 three-dimensional (3D) structures, in. pdb format were obtained from the RCSB PDB database (https://www. rcsb.org/). The AutodockTools1.5.6 software deleted water, cleaned the protein, changed the forcefield on the .pdb file, saved it in .pdbqt format, and was the receptor ready to complete. The ChemOffice and Openbabe 2. 3.1L preserved the 2D structure of the QRKSG active ingredient ligands and generated energy-minimized 3D conformations, saved in. pdb format. AutodockTools1.5.6 converted the .pdb files to .pdbqt format, and the ligands were ready for finishing.

Western Blot
The kidney tissues were thawed at 4°C in a refrigerator for protein extraction. Afterward, 20 mg of kidney tissues were homogenized with 200 μL lysis buffer. The total protein concentration was determined using the BCA kit (PICPI23225, Thermo Fisher Scientific, MA, United States) following the manufacturer's instructions. The proteins were separated using 10% SDS ammonium persulfate gel electrophoresis and transferred into Frontiers in Pharmacology | www.frontiersin.org December 2021 | Volume 12 | Article 765563 PVDF membranes. The membranes were blocked with 5% skimmed milk in TBST (Tris-buffered saline, 0.1% Tween 20) solution for 1 h at room temperature, followed by overnight incubation with the primary antibody PI3K (1:1,000 dilution), p-PI3K, AKT, p-AKT, mTOR, p-mTOR, BCL-2, Caspase-3 and reference Beta-actin (all 1:1,000 dilution) at 4°C. Next, the membranes were washed in TBST solution and incubated with Anti-Rabbit lgG conjugated to horseradish peroxidase (HRP) at 37°C for 1 h and washed in TSBT solution. The immune response area in the membranes was observed using the enhanced chemiluminescence (ECL) kit (WBKLS0100, Millipore, MA, United States) following the manufacturer's instruction. The density and immune response analysis was performed using the Tanon-5200 imaging system (Tanon Science and Technology Co. Ltd., Shanghai, China) and ImageJ Software (NIH, MD, United States).

Statistical Analysis
All quantitative group data were analyzed using a one-way analysis of variance (ANOVA) and presented as mean ± standard deviation. p < 0.05 was considered statistically significant.

Serum Biochemical Indexes and Urine Protein Quantitation
The TG, BUN, TC, and Cr levels in the model group serum increased, while ALB and TP decreased more than the control group. However, QRKSG treatment reduced the levels of TG, BUN, TC, and Cr in the serum, while the level of ALB and TP increased ( Figure 2). The 24 h proteinuria assessment in the rats revealed that the proteins in urine gradually increased from week 1 to week 4 after the model was established (Figure 3).
In week one, the urine protein content in the control group was significantly lower than the ADR-treated group, indicating that the NS was successfully induced. However, the proteinuria content was not significantly different between QRKSG and benazepril groups in the first week compared to the model group. Nevertheless, the proteinuria content was significantly higher in the model than in the QRKSG group in the 3rd and 4th weeks (p < 0.001). In addition, the proteinuria level was higher in the QRKSG than in the benazepril group in the 4th week. However, the urine protein levels were reduced in the QRKSG and benazepril groups compared to the control.

Histopathological Analysis
The pathological changes in the kidney tissues were observed under the light microscope ( Figure 4). The glomeruli and tubules structures were normal, and the control group had no visible damage. Nevertheless, the model group showed renal tubular abnormalities, dilation of the tubular lumen, and edema. However, the benazepril and QRKSG groups showed significantly reduced edema and morphological changes than the model group. Groups of rat glomerular podocytes were observed by transmission electron microscopy ( Figure 5). The results showed that foot processes in the normal group were neatly aligned, structurally intact, and the basement membranes were uniformly thick. The model group had numerous diffused fusions in the foot processes, and the structure of the podocytes changed, consistent with the structural changes of podocytes during NS development (Smoyer and Mundel, 1998). Moreover, the degree of fusion significantly reduced in the

Effects of QRKSG on Kidney Metabolism
The kidney samples obtained from the control, model, benazepril, and QRKSG groups were analyzed by GC-MS, and the total ion current chromatograms (TIC) are shown in Figure 6. The model and QRKSG groups showed a considerable degree of separation ( Figure 7C). The PLS-DA R 2 ( 0.963) and Q 2 ( 0.885) indicate an excellent PLS-DA model for discrimination and prediction ( Figure 7D). The kidney metabolism spectrum analysis by PCA and PLS-DA revealed that the metabolisms in QRKSG and benazepril groups ranged between the model and control groups. These results indicate that treatment caused normal metabolite levels in rat kidneys ( Figures 7A,B).

Metabolic Pathway Analysis
Twenty-four metabolites were identified as potential QRKSG biomarkers in the NS treatment at p < 0.05 and VIP > 1. Fifteen of the 24 metabolites were up-regulated in the QRKSG group, and nine were down-regulated than the model group ( Table 1). A heat map of the 24 metabolites showed that the differentially expressed metabolites in the QRKSG and benazepril groups changed significantly ( Figure 7E). Differential metabolic  pathway analysis was performed using the MetaboAnalyst5.0, with an impact value > 0.1 as the potential related metabolic pathway. The identified potential biomarkers mainly involved the metabolism of amino acids in the four main pathways: D-Glutamine and D-glutamate metabolism; alanine, aspartate, and glutamate metabolism; arginine and proline metabolism; and arginine biosynthesis ( Figure 7F; Table 2).

Component Analysis
Component analysis using UPLC-MS/MS identified 411 compounds, including 142 flavonoids, 149 phenolic acids, 16 terpenoids, 29 alkaloids, 27 lignans, and coumarins, 3 quinones, and 45 others were identified. Thirty components were relatively regulated in the QRKSG group based on their retention time and mass spectrum ( Table 3 and Supplementary Figure S1). Additionally, the size of the fan-shaped area for each compound in the fan diagram corresponded to the peak area (Supplementary Figure S2). The QRKSG group had the highest contents of Cynarin, Quercetin-4'-O-glucuronide, esculin, and Cichorium.

Metabolomics and Network Pharmacology
The network diagram revealed 34 node interactions (6 pathways, 5 metabolites, 15 enzymes, and 8 targets) (Figure 9). The greater the degree, the more the connected edges, hence the node was more important in the network. The average degree of the candidate components was 8.52, with AKT1 (degree 26), GAPDH (degree 21), CASP3 (degree 21), SRC (degree 18), BCL2L1 (degree 15), and mTOR (degree 13). AKT1, BCL2L1, and mTOR were involved in the PI3K-Akt signaling pathway, protecting against ADR-induced podocyte liver damage (Dai et al., 2010). Other targets, including the apoptosis-related proteins such as CASP3 and BCL2L1 (Pradeep and Srinivasan, 2018), were related to renal tubular cell apoptosis (Ahmed et al., 2019). Therefore, AKT1, BCL2L1, CASP3, and mTOR were the core QRKSG targets for treating NS. Moreover, L-Glutamic acid and ornithine (both degree 13) regulated multiple targets in the network diagram, indicating that QRKSG exerts therapeutic effects by regulating multiple metabolic pathways through multiple targets and metabolites.

DISCUSSION
Adriamycin-induced nephropathy (AN) has been extensively studied to unravel the underlying processes leading to kidney damage (Burke et al., 1977;Lu et al., 2017). In this study, the NS rat model was established by injecting ADR through the tail vein to induce NS, leading to an increased proteinuria content in the model than the control group. Renal histopathology results revealed significant morphological changes in the NS rat model manifesting as renal tubular abnormalities, dilation of kidney tubular lumen, edema, and diffuse fusion of foot processes, suggesting a successful NS modeling. The presence of proteins in the urine is a sign of podocyte damage, an indicator of glomerular filtration dysfunction, and damage of renal tubular reabsorption (Reiser et al., 2002). However, QRKSG treatment significantly reduced proteinuria, which manifested as reduced edema, morphological changes, and foot processes fusion in the renal tubes of rats in the QRKSG than the model group. Therefore, QRKSG protected rats against NS. The GC-MS metabolomics established the changes in renal metabolism, thus revealing the mechanism of QRKSG for treating NS. Twenty-four metabolic biomarkers were identified between the model and QRKSG groups. These metabolites were mainly involved in D-glutamine and D-glutamate metabolism; alanine, aspartate, glutamate metabolism; arginine and proline metabolism. The levels of L-glutamine and L-glutamate decreased after QRKSG treatment, indicating that the metabolism of D-glutamine and D-glutamate were interrupted. This observation is consistent with other findings where glutamine and phenylalanine levels were down-regulated in idiopathic NS (Sedic et al., 2014). Glutamine and glutamate are irreplaceable in FIGURE 9 | Metabolites-targets-pathways network.
Frontiers in Pharmacology | www.frontiersin.org December 2021 | Volume 12 | Article 765563 D-glutamine and D-glutamate metabolism. Glutamine is catalyzed to reversibly generate glutamate, which further generates α-ketoglutarate, participating in the tricarboxylic acid (TCA) cycle to produce adenosine triphosphate (ATP) (Bharadwaj et al., 2020). Thus, contributing to material and energy production (Hakimi et al., 2016). However, glutamine plays a different role as an essential amino acid in the cell cycle progression, activating the mammalian target of rapamycin (mTOR) and homeostasis activation of reactive oxygen species (ROS) (Still and Yuneva, 2017;Bernfeld and Foster, 2019). Glutamine also regulates the mTOR pathway and coordinates various environmental factors regulating cell growth and homeostasis (Still and Yuneva, 2017;Pradeep and Srinivasan, 2018). In addition, glutamine synthetase is important in renal ammonia metabolism, maintaining the acid-base balance in the whole body, essential for optimal health (Conjard et al., 2003;Verlander et al., 2013;Harris and Weiner, 2021). In mice, glutamate inhibits the fusion of podocytes with ADR, nephropathy, cell apoptosis, renin expression and reduces urine albumin (Gu et al., 2012). Besides, glutamate converts toxic ammonia into non-toxic glutamine, preventing ammonia poisoning (Tardito et al., 2015). Therefore, QRKSG probably metabolized glutamine and glutamate, reducing the urine protein in rats with NS, thereby protecting the kidneys of these rats. Network pharmacology analysis using UPLC-MS/MS predicted 146 targets and 153 pathways, with the PI3K-Akt signaling pathway, and the AGE-RAGE signaling pathway in diabetic complications, as the main pathways. AKT1, CASP3, BCL2L1, and mTOR were identified as potential QRKSG therapeutic targets for NS treatment. The PI3K/AKT signaling pathway is key in mammalian cell differentiation, proliferation, apoptosis, survival, autophagy and is a pivotal regulator of autophagy (Yan et al., 2018;Kim et al., 2020). PI3K phosphorylates AKT after activation, activates downstream signaling molecules such as mTOR, and exerts corresponding biological effects. mTOR is the core protein that determines the formation and maturation of autophagosomes . mTOR activation in the kidney promotes up-regulation of proinflammatory and pro-fibrotic factors, leading to tubular interstitial fibrosis and atrophy (Lieberthal and Levine, 2012). In contrast, inhibiting mTOR activation protects the kidneys against podocyte apoptosis . Inhibiting the PI3K/ Akt/mTOR pathway induces autophagy (Tsai et al., 2015), promoting cell health in response to various cell stresses and nutritional conditions and subsequently protecting cells from damage (Kitada et al., 2017). For example, mTOR activation negatively regulates autophagy in several types of cells (Lindqvist et al., 2015;Liu et al., 2018). Autophagy is key in preventing kidney diseases because it maintains homeostasis and the integrity of renal tubules (Peng et al., 2019). Thus, upregulating autophagy is an efficient treatment strategy for kidney diseases (Bhatia and Choi, 2020). Besides, autophagy substantially reduced renal inflammation and interstitial damage in a unilateral ureteral obstruction (UUO) model where mTOR inhibition significantly increased autophagy and suppressed UUO-induced inflammation (Du et al., 2017). In the model group, PI3K, AKT1, and mTOR expression levels increased in the QRKSG than the control group, indicating that the PI3K/Akt/mTOR signaling pathway was activated. However, QRKSG treatment inhibited the PI3K/Akt/mTOR signaling pathway, which indicates that QRKSG probably increased autophagy to prevent renal tubular interstitial fibrosis and atrophy. Bax, a member of the BCL-2 family, is up-regulated by p53 transactivation, while p53 transcriptionally suppresses BCL-2, thereby activating caspase-9 and subsequently caspase-3 (Wen et al., 2014;Wnek et al., 2017). Increasing the ratio of BCL-2 to Bax achieves the anti-apoptotic function of BCL-2 (Castro- Gonzalez et al., 2021), while the activation of caspase-3 promotes cell apoptosis (Sun et al., 2011). In addition, Srinivasan and Pradeep established that increasing the apoptosis mediator Bax decreased BCL-2 in the kidneys of diabetic rats (Pradeep and Srinivasan, 2018). Parallel to the present study, Osama et al. reported that paricalcitol and enalapril protective cells against streptozotocin (STZ)induced renal tubular cell apoptosis by reducing the caspase-3 expression and increasing BCL-2 expression (Pradeep and Srinivasan, 2018). Wu et al. also established that Fang Qi Dan down-regulates the expression of Bax protein while increasing the BCL-2 protein, which protects podocytes, and restores the selective filtration function of the glomerulus (Wu et al., 2014). The expression of CASP3 increased in the model group, while the expression level of BCL2L1 decreased (VS control group), indicating the increased ability of renal cell apoptosis. However, QRKSG treatment weakened the renal apoptotic capacity. These results indicate that QRKSG protects renal tubular cells against ADR-induced renal damage. The findings are consistent with previous research (Lieberthal and Levine, 2009;Kaushal et al., 2020;Zhou et al., 2020;Bhosale et al., 2021;El-Gowily et al., 2021), and the relevant target graphs are shown in Supplementary Figure S3. Therefore, QRKSG could promote autophagy and anti-apoptosis through the expression of AKT1, CASP3, BCL2L1, and mTOR, thereby conferring protection to the podocytes and maintaining renal function and tubular homeostasis.
Molecular docking revealed the mode of interaction between the active QRKSG components and target proteins in the kidneys, where most components are bound to core targets (AKT1, CASP3, BCL2L1, and mTOR). The comprehensive evaluation of the docking results initially revealed the QRKSG material basis FIGURE 11 | Western blot analysis of PI3K, AKT1, mTOR, BCL-2, and CASP3. (A) Representative images of PI3K, AKT1, mTOR, BCL-2, and CASP3 (B) PI3K, AKT1, and mTOR protein expressions (C) BCL-2 and CASP3 protein expressions. The data are presented as mean ± SD, **p < 0.05 shows comparison against the control group, while # p < 0.05 shows comparison against the model group.
Frontiers in Pharmacology | www.frontiersin.org December 2021 | Volume 12 | Article 765563 of treating NS. This observation has profound implications for exploring the molecular-level QRKSG mechanism of action.

CONCLUSION
A comprehensive analysis revealed 23 key targets related to five metabolites and six pathways. These potential targets were related to multiple QRKSG active components. Among the key targets, AKT1, MTOR, CASP3, and BCL2L1 were important for the QRKSG underlying mechanism in treating NS. Specifically, the potential QRKSG mechanisms include inhibiting cell apoptosis and increasing autophagy, D-Glutamine, and D-glutamate metabolism, reducing proteinuria and weakening kidney damage ( Figure 12). Therefore, these findings on the potential underlying QRKSG mechanism provide a theoretical basis for QRKSG clinical application.

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

ETHICS STATEMENT
The animal study was reviewed and approved by the Experimental Animal Ethics Committee of Hubei University of Chinese Medicine.

AUTHOR CONTRIBUTIONS
QY, JL, and LH designed the study; YD, DZ, SZ, PH, FZ, and YY performed the research; GM and FH collected and analyzed the data; YD drafted the paper.