Exploring the Potential Mechanism of Qi-Shen-Di-Huang Drug Formulary for Myasthenia Gravis (MG) based on UHPLC-QE-MS Network Pharmacology and Molecular Docking Techniques

Myasthenia gravis (MG) is a rare and refractory autoimmune disease, and Qi Shen Di Huang (QSDH) drug formulary is an in-hospital herbal decoction with proven clinical efficacy in treating MG. Currently, most of the research on the QSDH drug formulary has concentrated on its clinical efficacy, and there is a lack of systematic study on the material basis. The active compounds and their mechanism of action have not been entirely determined. Therefore, this study sought to identify the active compounds in the QSDH drug formulary and analyze the key targets and potential mechanisms. We used ultra-performance liquid chromatography Q Exactive-mass spectrometry (UHPLC-QE-MS) and Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP) database to identify and screen 85 active ingredients corresponding to 59 potential targets (17 herbs) associated with myasthenia gravis, and further identified AKT1 as the primary core target and the PI3K/AKT signaling pathway as the most substantial enriched pathway. Molecular docking and UPLC-MS analysis identified quercetin, luteolin, wogonin, kaempferol, laccasein, and epigallocatechin gallate are the core compounds of the QSDH drug formulary. In vivo rat studies showed that the QSDH drug formulary reduced Lennon's clinical score and decreased acetylcholine receptor antibody levels in peripheral blood rats with experimental autoimmune myasthenia gravis. In addition, the QSDH drug formulary downregulated P-PI3K/PI3K and P-Akt/Akt protein expression. Collectively, these findings describe the role and potential mechanism of the QSDH drug formulary in the treatment of MG, which exerts potential value by acting on AKT targets and regulating the PI3K/AKT signaling pathway and providing a theoretical reference for QSDH drug formulary application in the clinical treatment of MG.


Introduction
Myasthenia gravis (MG) is an autoimmune disease with lesions in the postsynaptic membrane of the neuromuscular junction (NMJ). Pathogenic antibodies include the common acetylcholine receptor antibody (AChR-ab) and the less common muscle-specific receptor tyrosine kinase antibody (Musk-ab), and the low-density lipoprotein receptor-related protein4 antibody (LRP4-ab) blocks the aggregation of acetylcholine receptors (AChR) and disrupts the function of AChR and the signaling between NMJ [1]. Fatigue and fluctuating muscle contraction are the core clinical symptoms, which may involve skeletal muscles throughout the body, including the eyes, medulla oblongata, respiratory, and extremity muscles. e clinical presentation also varies depending on the autoantibody type and a thymoma's presence [2]. e global prevalence of MG is (150∼250)/ million, and the prevalence of MG in China is about 0.68/per 100,000, slightly higher in women than in men. e whole body's skeletal muscle can be involved, and the symptoms are "light in the morning and heavy in the evening," aggravated by activity and relieved by rest [3]. e treatment of MG is still mainly composed of cholinesterase inhibitors, glucocorticoids, immunosuppressants, intravenous immunoglobulins, and plasma exchange but traditional Chinese medicine has unique advantages in the treatment of MG, improving clinical efficacy while reducing the dosage and side effects of western drugs and improving the immune balance of the body [4]. e QSDH drug formulary is based on the theory of "Fu Xie" and "Nao Sui" which proposes that the fundamental mechanism of this disease is "a deficiency of spleen and kidney and deficiency of brain marrow," and the treatment should be essential in the basic principle of "JianPi YiQi BuSui" [5]. e metabolic profile of TCM is the key to pharmacological research and clinical application. QSDH drug formulary consists of 17 herbs, including Astragalus, Radix Codonopsis, Atractylodes, Radix et Rhizoma, and others. e current research on the QSDH drug formulary is mainly focused on clinical efficacy, lacking systematic material basis research, in-depth research on its overall chemical composition, the material basis of pharmacological efficacy, activity, and its active ingredients and specific mechanisms not fully understood. TCM research needs some new approaches and various network technologies have been evolving with the development of the information age. Cyberpharmacology has recently been an emerging discipline based on systems biology, multidirectional pharmacology, and proteomics. It combines pharmacology and information networks to construct molecular biological networks of diseases and drugs by analyzing the whole network to elucidate the pharmacological effects of drugs and pathogenesis of diseases through high-throughput histological data analysis and search of network databases [6]. We can directly identify drug and disease targets through network pharmacology from a large amount of data and understand the mechanisms and pathways [7]. Ultrahigh-pressure liquid chromatography coupled with high-resolution mass spectrometry has the advantages of solid separation and structure identification, which can reasonably predict the structure of unknown compounds. e combination of UHPLC-MS and network pharmacology has dramatically promoted the novel pharmacological research and drug development of TCM. erefore, in our study (Figure 1), we used UHPLC-QE-MS to identify the active ingredients of the QSDH drug formulary and then screened the potential targets of action of the QSDH drug formulary for treating MG by network pharmacology. Finally, molecular docking techniques and animal experiments to validate the key targets and pathways of the QSDH drug formulary and its essential substances acting on MG.

Metabolites Extraction.
e samples were crushed with a mixer mill for 30 s at 65 Hz. 100 mg of sample was added to 500 μL of an extracted solution containing 1 μg/mL of internal standard and dissolved in 80% methanol. After the 30s vortex, the samples were homogenized at 45 Hz for 4 min and sonicated for 1 hour in an ice-water bath. After 1 hour in −40°C, the sample was centrifuged at 12000 rpm for 15 min at 4°C. Finally, the supernatant was obtained and put in a fresh 2 mL tube for LC-MS/MS analysis, taking 40 ul from each sample and pooling it as QSDH samples.

LC-MS/MS
Conditions. LC-MS/MS analysis was performed on an Agilent ultra-high performance liquid chromatography 1290 UPLC system with a Waters UPLC BEH C18 column (1.7 μm 2.1 * 100 mm). e flow rate was set at 0.4 mL/min, and the sample injection volume was set at 5 μL.

QSDH Drug Formulary Screening of Active Ingredients and Targets of Action.
According to the principles of toxic pharmacokinetics (ADME), the TCMSP [8] (https:// tcmspw.com/tcmspsearch.php) database was used according to the conditions of oral bioavailability (OB) ≥ 30 and drug-like likeness (DL) ≥ 0.18. e active ingredients of QSDH drug formulary detected by UHPLC-QE-MS were screened. e targets corresponding to the chemical composition of the active ingredients were retrieved through the TCMSP database, and the target names were transformed to gene abbreviations through the UniProt (https://www. uniprot.org/) database [9].

Construction Network of the "TCM Active Ingredients -
Intersection Target -Disease.". In order to determine the MG target of the QSDH drug formulary, the intersection target of the MG target and the corresponding target of the QSDH drug formulary active compound were obtained by Perl language analysis, and a Wayne diagram was drawn. e data of the intersection targets of the active compounds and their interactions were imported into Cytoscape3.8.2 software to construct a network of "TCM active componentsintersectiontarget-disease," and topological analysis was conducted to explore the pharmacological mechanism of the active components of QSDH prescription.

QSDH Drug Formulary Predictive Target Protein Interaction (PPI) Network Construction and Analysis.
After importing the intersection targets obtained above into the STRING database (https://www.string-db.org) [15], multiple proteins were selected to constrain Homo sapiens, hide-free nodes, and retain default values for other parameters to obtain the intersection target interaction network. At the same time, the data were set with medium confidence >0.4 to ensure the reliability of the analysis.

GO and KEGG Enrichment Analysis.
e intersection target was imported into the Hiplot database (https://hiplot. com.cn/, updated on May 4, 2021), and the species was Homo sapiens. e threshold P < 0.01 and the Gene Ontology, GO enrichment analysis, and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis were set. R 4.0.3 software was used to draw bubble charts and bar charts from obtaining the biological process (BP), molecular function (MF), cellular components (CC), and main signal pathways involved in the treatment of MG by QSDH drug formulary.

Molecular Docking Verification of Key Active Ingredients
and Core Targets. Based on the above research results, the network analyzer plug-in is used to analyze the network Evidence-Based Complementary and Alternative Medicine topology structure. Active components with the most significant degree value and the largest number of targets are selected as ligands, and core key genes are screened as receptors by the Cytohubba plug-in. e PubChem database (https://pubchem.ncbi.nlm.nih.gov/) downloads the 3D structure of active ingredients in the PDB database (http:// www1.rcsb.org/) to obtain the target protein 3D structure, AutoDockTools1.5.6, AutoDock Vina and PyMOL software were used for molecular structure processing and molecular docking.

Experimental Animal Selection.
Sixty female SPF Lewis rats, weighing 160∼180 g, aged 6∼8 weeks, were provided by Beijing Viton Lever Experimental Animal Technology Co., LTD., certificate No. SCXK (Beijing) 2016-0006. ey were fed in the SPF barrier laboratory of the Animal Experiment Center of Changchun University of Chinese Medicine with free drinking water, indoor keep (21 ± 2)°C, relative humidity 50%-60%, and day and night alternating light and dark for 12 h. e experimental procedure was reviewed and approved by the Experimental Animal Ethics Committee of the Changchun University of Chinese Medicine (NO.2021228).

Rat Model Preparation and Grouping.
Eight of 60 Lewis female rats were chosen randomly as the adjuvant control group, and the remaining 52 were referred to Baggi et al [16]. e EAMG rat model was established by the active immunization method. Each rat was injected with 200 μL of mixed emulsion, and Rα97-116 peptide dry powder and Mycobacterium tuberculosis H37RA powder were added to PBS. Each 200μL of PBS contained 50 μg of Rα97-116 dry powder and 1 mg of H37RA dry powder, and the peptide and dry powder were dissolved by ultrasonication on an ice box to make the peptide and dry powder fully fused in PBS until the solution was cloudy. e dissolved liquid is extracted into a 1 mL syringe. Another syringe extracts 200 μL of CFA, and a medical hose connects the two syringes to mix the PBS and CFA of the meaning peptide segment and bacteriophage powder to prevent protein denaturation, the whole process is carried out on the ice, repeatedly pumping two syringes back and forth until milky white liquid appears and the emulsion does not dissolve and does not spread, keeping the spherical shape floating on the water surface for a long time, then the mold-making drug is drug formularyted successfully. e equipped mold-making antigen emulsion was injected subcutaneously at five sites: the tail's base, both hind limbs' foot pads, and both sides of the back. e adjuvant control group was injected with a mixture of CFA and PBS in equal amounts, and the day of the first immunization was recorded as day 0. One booster immunization was given on day 30 and day 45, with the same preparation, injection volume, injection sites, and injection method as the first immunization, except that CFA was replaced with IFA. e adjuvant control group was injected with a mixture of IFA and PBS, with the same injection volume, injection sites, and injection method as the first immunization. e dose, injection site, and injection method were the same as the first immunization. e model was evaluated on day 60, and the change in body mass assessed the success, Lennon's clinical symptom grading method was used [17], and the difference in the level of AChR-Ab in the peripheral serum of the tail vein of rats. Furthermore, Lewis rats meeting the EAMG criteria were divided randomly into eight rats in the model group, the QSDH drug formulary low-dose group (4.8 g/kg/d), the QSDH drug formulary medium-dose group (9.6 g/kg/d), the QSDH drug formulary high-dose group (19.2 g/kg/d), and the prednisone acetate tablet group (5.4 mg/kg/d).

Western Blot.
e removed rat spleen tissue was cut and put into RIPA lysis solution, lysed on ice for 30 min at 12000 r/min, centrifuged at 4°C for 10 min, the supernatant was taken, the protein concentration was measured by BCA method, the sample protein concentration was adjusted to the same, the protein samples were separated by 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), and then transferred to polyvinylidene fluoride (PVDF) membranes. After the end of the transfer, PI3K, p-PI3K, Akt, p-AKT primary antibodies (1 : 500), and secondary antibodies were incubated with 5% skimmed milk for 1 hour at room temperature. e transferred protein strips were developed by exposure using a chemiluminescence imaging system, and the images were acquired by quantitative optical density analysis.

Statistical
Analysis. GraphPad Prism 5.0 software was used for statistical analysis, and measures were expressed as mean plus or minus standard deviation (x ± s) using oneway ANOVA. P < 0.05 indicates that the difference is statistically significant.

UPLC-QE-MS Detection for QSDH Drug Formulary.
A total of 930 compounds in 60 categories, including quinones, terpenoids, flavonoids, and phenols, were detected by UPLC-QE-MS. e corresponding pharmacokinetic parameters (ADME) were searched and downloaded from the TCMSP website. According to OB ≥ 30% and DL ≥ 0.18 as screening conditions, 85 active compounds were obtained, as shown in Table 1.

Screening for Targets Related to MG.
A total of 822 targets related to myasthenia gravis were obtained from the Gen-eCards, Drugbank, DisGeNET, TTD, and PharmGKB databases using the keyword "myasthenia gravis," as shown in

Construction of the "TCM Active Ingredients -Intersection
Target -Disease Network". Ased on Perl language analysis of 85 targets of active components of QSDH prescription and corresponding targets of MG disease, 59 intersection target proteins were obtained, considered as targets of QSDH drug  Figure 3. rough Cytoscape3.8.2, the "TCM active componentintersectiontarget-disease" network is constructed (Figure 4). e network has 112 nodes and 320 edges. e degree value represents the number of edges connected to this point. e larger the degree value is, the more extensive the effect of this node on other nodes is. e red triangle represents MG disease, the blue rectangle represents active chemical components of the QSDH drug formulary, and the orange ellipse represents the intersection target of active components of herbs and MG. Further topological analysis (based on degree and betweenness centrality of network nodes) shows that the key chemical components in this network are quercetin (MOL000098) and epigallocatechin-3-gallate (MOL006821), luteolin (MOL000006), wogonin (MOL000173), kaempferol (MOL000422), and fisetin (MOL013179).
is disease-TCM active ingredient-target network is characterized by complex components, multiple targets, and close interaction between components and targets, forming a complex network for the QSDH drug formulary treatment of MG.

Construction and Analysis of a PPI Network of Predictive Targets for MG Treatment with the QSDH Drug Formulary.
To better explain the mechanisms by which the QSDH drug formulary acts concerning MG, we further evaluated its 59 intersecting targets to analyze their relationship. e 59 intersecting targets were imported into the STRING database to obtain protein interaction information, which was visualized using Cytoscape 3.8.2 software ( Figure 5). ere are 58 nodes and 1172 edges in the intersection target network. e nodes represent the targets of the active ingredient of the QSDH drug formulary for the treatment of MG, and the larger the mutuality value between the targets, the darker the color of the nodes, indicating the more critical the nodes are. e CytoNCA plug-in was used further to perform core screening by degree centrality (DC), closeness centrality (CC), and betweenness centrality (BC) and to construct a protein interaction network by degree value, see Figure 6, to obtain the final core targets: AKT1, STAT3, IL1B, and TP53. ey suggested that the QSDH drug formulary may exert the therapeutic effect of MG by regulating these proteins.

GO Enrichment Analysis.
e 59 intersecting targets were entered into the Hiplot online database, set p-value<0.01 and q-value<0.05, e top 10 entries obtained from GO enrichment analysis were performed, and bubble plotted using R4.0.3 software. A total of 1959 items were significantly enriched (P < 0.05), of which 1834 were enriched in biological process (BP), 35 in cellular component (CC), and 90 in molecular function (MF), Figure 7. It is mainly related to T cell activation, nutrition level reaction, response to lipopolysaccharide, chemical stress reaction, physical stimulation, oxidative stress, cytokines receptors, cytokine activity, G proteincoupled receptor of neurotransmitter receptor activity, passing the gated ion channel activity involved in the regulation of the postsynaptic membrane potential, and neurotransmitter receptors postsynaptic activity.

KEGG Signaling Pathway Enrichment Analysis.
KEGG signaling pathway enrichment analysis of 59 intersecting targets was performed by Hiplot online database, and the results showed that they were mainly enriched in 157 signaling pathways (P < 0.05), and the top 30 pathways were bar charted by R 4.0.3 software (Figure 8). e signaling pathways closely related to MG mainly included the IL-17 signaling pathway, the HIF-1 signaling pathway, the PI3K/AKT signaling pathway, and  Evidence-Based Complementary and Alternative Medicine the mTOR signaling pathway. e above results suggest that QSDH can be used to treat myasthenia gravis by regulating multiple biological processes and the coordinated action of multiple signaling pathways.

Molecular Docking Verification.
Six central herbal core ingredients in the QSDH drug formulary were quercetin, epigallocatechin-3-gallate, luteolin, wogonin, kaempferol, fisetin, and epigallocatechin-3-gallate, respectively, according to the "active ingredient-intersectiontarget-disease" network analysis. eir binding ability was predicted with AKT1, the top core target screened by the CytoNCA plug-in. See Figure 9 for the chemical structure of the core active components of herbal medicine. Binding energy < −4.25 kcal/mol indicates some binding activity between the small ligand molecule and the receptor protein, binding energy < −5.0 kcal/mol indicates good binding activity between the two, and binding energy < −7.0 kcal/mol indicates that the ligand has binding activity with the receptor binding solid activity [18]. It was suggested that the more stable the binding conformation of the small molecule receptor to the ligand, the lower its energy and the greater the interaction produced [19]. According to the molecular docking binding energy scores (Table 2), it is clear that the binding energy of quercetin, luteolin, wogonin, kaempferol, fisetin, and the above five compounds with the critical core action target AKT1 are all less than −5.0 kcal/mol, indicating an excellent bonding activity between them. e binding energy of AKT1 to epigallocatechin-3-gallate was less than /−7.0 kcal/mol, indicating a solid binding activity between the two. Molecular docking visualization was performed separately using Pymol 2.4.2 software (Figure 10). Epigallocatechin-3-gallate formed hydrogen bonds with five amino acids, ASP-46, GLU-40, LYS-39, ALA-50, and GLN-47, near the active site to bind to AKT1. Figure 11, during the first 15d after the first immunization, there was still no significant change in Lennon's clinical symptom score. However, the EAMG-modeled rats began to show a gradual decrease in grip strength and activity. From the 15th day after the first immunization, the rats in the model group began to show significantly weaker bilateral forelimb grasping power, reduced antagonism, less flexible body than the adjuvant control group, and preferred to curl up in the corner of the cage. After the second booster immunization (the 30th day after the first immunization), the EAMG symptoms gradually increased (P < 0.001).

Changes in Serum AChR-Ab of EAMG Model Rats.
On day 60 after the first immunization, the serum AChR-Ab levels of rats were analyzed by ELISA. As shown in Figure 12, the serum AChR-Ab levels of rats in the EAMG modeling group were significantly higher than those in the adjuvant control group (P < 0.01).

Effect of QSDH Drug Formulary on Lennon's Clinical Symptom Score in EAMG Rats.
e EAMG model was successfully established and then randomized to group administration of the intervention, with Lennon clinical score tests performed every seven days during treatment. During the administration treatment, the clinical symptoms of rats in the EAMG model group did not improve as the experiment progressed. Even individual rats showed signs of paralysis, convulsions, loose skin, jaws clinging to the bottom of the cage, and inability to eat, and the Lennon clinical score was significantly different from that of the adjuvant control rats (P < 0.001), as shown in Figure 13. Indicates that the EAMG model was established stably and did not heal during the administration period. As the experiment continued, the clinical symptoms of the rats in the high, medium, and low doses of QSDH drug formulary and prednisone acetate administration groups were gradually relieved, and the Lennon clinical score tended to decrease and was significantly different from that of the EAMG model group (P < 0.001).

Effect of QSDH Drug Formulary on Serum AChR-Ab Levels in EAMG Rats.
At the end of the dosing cycle, the rats were anesthetized with gas anesthesia apparatus, taken through the abdominal aorta, and the serum was separated by centrifugation. e ELIAS method detected the peripheral serum AChR-Ab content of the rats in each experiment group. e serum AChR-Ab levels of rats in the EAMG model group were significantly higher than those in the adjuvant control group (P < 0.01). After the intervention of each dosing group, the serum AChR-Ab levels of rats in the QSDH drug formulary low, medium, and high-dose groups and the prednisone acetate tablet group were significantly reduced compared with those in the EAMG model group (P < 0.05, P < 0.01, Figure 14).

QSDH can Inhibit the PI3K/AKT Signaling Pathway.
Based on the core targets screened by network pharmacology and the results of the KEGG enrichment analysis, we further verified whether the QSDH drug formulary for MG treatment interferes with the PI3K/AKT signaling pathway. Western blot was used to detect the expression of the PI3K/ AKT pathway-related proteins in the spleen tissues of each group of rats. Detection of protein expression of PI3K, phosphorylated PI3K (p-PI3K), AKT, and phosphorylated AKT (p-AKT). As Figures 15(a)-15(b) show, the expression of pathway-related proteins p-PI3K and p-AKT in the spleen tissue of rats in the EAMG model group was significantly higher (P < 0.01) than in the adjuvant control group. e pathway-related proteins p-PI3K and p-AKT decreased to  Evidence-Based Complementary and Alternative Medicine different degrees after each dose group of QSDH drug formulary low, medium, and high and positive drug treatment (P < 0.05, P < 0.01).

Discussion
Network pharmacology is based on systems biology theory, integrates pharmacology, high-throughput sequencing, genomics, and other technologies, emphasizing the multichannel regulation of signaling pathways [20]. Our study explores the mechanism of action and potential therapeutic targets of the QSDH drug formulary for MG from the perspective of network pharmacology. Using ultra-performance liquid chromatography Q Exactive-mass spectrometry, network pharmacology, and molecular docking technology, 930 chemical components were analyzed from the QSDH drug formulary. e "active component-intersectiontarget-disease" network of QSDH drug formulary treatment for MG was constructed. Eightyfive active components and 322 potential targets of quercetin, luteolin, wogonin, kaempferol, fisetin, and epigallocatechin-3-gallate were screened. AKT1, STAT3, IL1B, and TP53 were selected as core targets through the PPI network, and AKT1 was the first and most core target. GO, and KEGG enrichment analysis further explored the biological pathways involved in the expression of potential target proteins, which were enriched in the IL-17 signaling pathway, the HIF-1 signaling pathway, the PI3K-Akt signaling pathway, and the mTOR signaling pathway. PI3K-AKT was the most abundant protein enrichment signaling pathway. e binding activity of the core active ingredient and core target protein in the QSDH drug formulary were verified. Finally, the rat experiment verified that AKT1 might be the main target of QSDH prescription in the treatment of MG, and QSDH prescription can play a role in the treatment of MG by inhibiting the PI3K/AKT signaling pathway. e active ingredient-intersection target network showed that quercetin, epigallocatechin-3-gallate, luteolin, wogonin, kaempferol, and fisetin matched more targets, which may be the core compounds and essential material basis for the QSDH drug formulary to exert its therapeutic effect on MG. Studies have shown that epigallocatechin-3gallate, quercetin, luteolin, wogonin, kaempferol, and fisetin have anti-inflammatory, antioxidant, immunomodulatory, and neuroprotective effects [21][22][23][24][25][26]. Epigallocatechin-3gallate promotes dendritic growth and synaptogenesis through selective inhibition of the AKT/mTOR signaling pathway, inhibits acetylcholinesterase activity, and effectively ameliorates neurological damage in neurodegenerative diseases. In addition, it can regulate immune and inflammatory responses by influencing various aspects of immune cell function in innate and adaptive immunity. It   can affect the differentiation of initial CD4+T cells to different effector subpopulations and thus treats autoimmune inflammatory diseases [27]. Quercetin has a wide range of pharmacological effects, including anti-inflammatory, antiviral, antibacterial, and immunomodulatory effects, and can significantly upregulate IFN-c in 1 and downregulate IL-4 expression in 2, thereby regulating immune function [28]. Luteolin has antitumor, anti-inflammatory, and cerebrovascular protective effects, and its related mechanism involves the regulation of PI3K/AKT, MAPK, and other signaling pathways as well as the expression of related cytokines and kinases [29]. An experimental study found that kaempferol acts the same way as the immunosuppressant cyclosporine A, mediating calcium-regulated neural phosphatase activity from inhibiting immune cell differentiation [30]. Fisetin reduces IL-17 production by CD4+T lymphocytes through inhibition of the PI3K/AKT/mTOR pathway and decreases 1/ 17 pro-inflammatory cytokine secretion in human peripheral blood mononuclear cells [31]. e above can indicate that the active ingredients in the QSDH drug formulary can regulate the subpopulation differentiation of CD4+T cells, alter the secretion of related cytokines between 1/ 2/ 17, and exert antiinflammatory, immunomodulatory, and neurological function protection effects. e core targets obtained from screening intersecting targets in the PPI network by Cytohubba and CytoNCA plug-ins, of which the top three occupied by AKT1, STAT3, and TP53, were predicted to be the possible main targets of action of the QSDH drug formulary for the treatment of MG. Studies have reported that AKT1 enhances the differentiation of CD8+T cells, improves the proliferation of CD8+T cells, and inhibits the over-activation of DCs, thus regulating the immune response [32].
e study showed increased expression of 1 and 17 cells and their associated cytokines IL-1, IL-6, IL-17, IFN-c, and TNF-α in the peripheral serum of MG patients [33]. e expression of IL-6 can be regulated by inhibiting the AKT/mTOR signaling pathway to improve the symptoms of muscle fatigue weakness in MG patients [34]. STAT3, which has immunomodulatory effects, modulates the production of 1 cellspecific cytokines by altering the balance between 1 and 2 cells, thereby altering immune function and inflammatory responses [35]. TP53 is currently the most widely studied oncogene. TP53 can be both an activator and an inhibitor of autophagy. Under stress conditions such as nutrient depletion or hypoxia, TP53 promotes autophagy activation by inhibiting the mTOR signaling pathway, thereby suppressing the inflammatory response. Inhibition of TP53 increases IL-6 gene expression while promoting cell proliferation and differentiation [36].
Based on the KEGG pathway enrichment analysis results, we can see that the QSDH drug formulary exerts its effect on the treatment of MG mainly through signaling pathways such as the IL-17 signaling pathway, the HIF-1 signaling pathway, the PI3K-AKT signaling pathway, and the mTOR signaling pathway. e PI3K-AKT signaling pathway regulates cell proliferation, differentiation, apoptosis, and glucose transport and can regulate T cells' development, stability, and function [37]. Studies suggest that PI3K-AKT-specific inhibitors significantly upregulate 1 and 17 cells in PBMCs of MG patients, and mTOR/HIF-1α activates the expression of glycolytic genes to provide energy for rapid activation of immune cells, suggesting that the PI3K/AKTsignaling pathway regulates immune function in MG [38]. Our in vivo experimental study in rats confirmed that the QSDH drug formulary could alleviate the symptoms of EAMG and reduce the level of acetylcholine receptor antibodies in the peripheral blood of EAMG rats by inhibiting the PI3K/AKT signaling pathway, thus demonstrating that the QSDH drug formulary for MG can work by inhibiting the PI3K/AKT signaling pathway. We selected the core active ingredients quercetin, epigallocatechin-3-gallate, luteolin, wogonin, kaempferol, and fisetin as ligands and the first core target AKT1 as the receptor for molecular docking. Furthermore, we verified the pharmacodynamic substance basis of the QSDH drug formulary for MG treatment. e results showed that the main active compounds could stably bind to the core target AKT1 through hydrogen bonding. Among them, epigallocatechin-3-gallate formed the strongest hydrogen bond to AKT1 with asp-46, Glu-40, LYS-39, AlA-50, and GLN-47 amino acids near the active site. It further indicates that the QSDH drug formulary can treat myasthenia gravis by acting on the AKT1 target protein.

Conclusions
In conclusion, 85 active pharmaceutical ingredients of the QSDH drug formulary were identified by UHPLC-QE-MS combined with the network pharmacology method. Preliminary prediction QSDH drug formulary can regulate the PI3K/AKT signaling pathway through quercetin, epigallocatechin-3-gallate, luteolin, wogonin, kaempferol, fisetin, and other active ingredients, acting on AKT1 target protein to inhibit inflammation, immune function, and cell metabolism, to achieve the purpose of treating MG. Provides new insights into the QSDH drug formulary's mechanism to alleviate MG disease's progression.  Data Availability e datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Disclosure
Tianying Chang and Qi Lu and yibin zhang are co-first authors.

Conflicts of Interest
e authors declare that they have no conflicts of interest regarding the publication of this paper.