- 1Department of Pharmacy, Wenzhou TCM Hospital of Zhejiang Chinese Medical University, Wenzhou, China
- 2Department of Rehabilitation Medicine, Wenzhou TCM Hospital of Zhejiang Chinese Medical University, Wenzhou, China
- 3Wuyanling National Natural Reserve Administrative of Zhejiang, Wenzhou, China
- 4Department of Laboratory Medicine, Wenzhou TCM Hospital of Zhejiang Chinese Medical University, Wenzhou, China
- 5Department of Chinese Internal Medicine, Wenzhou Integrated Traditional Chinese and Western Medicine Hospital of Zhejiang Chinese Medical University, Wenzhou, China
- 6Department of Gynaecology, Wenzhou TCM Hospital of Zhejiang Chinese Medical University, Wenzhou, China
Objective: The incidence of Polycystic Ovary Syndrome (PCOS) is increasing annually. This study aims to investigate the therapeutic mechanisms of Yikang Decoction (YKD) in the treatment of PCOS through the integration of GEO datasets, network pharmacology, and dynamic simulation.
Methods: Active ingredients of YKD and their targets were collected from the Traditional Chinese Medicine Systems Pharmacology (TCMSP) platform. Disease-relevant targets for PCOS were retrieved from several databases, including GeneCards, OMIM, PharmGKB, DrugBank, and GEO. The underlying pathways associated with the overlapping targets between YKD and PCOS were identified using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. The mechanisms of interaction between the core targets and components were further explored through molecular docking and molecular dynamics simulations (MD).
Results: 139 potential active components and 315 targets of YKD were identified. A topological analysis of the PPI network revealed 10 core targets. These targets primarily participated in the regulation of biological processes, including cell metabolism, apoptosis, and cell proliferation. The pathways associated with treating PCOS encompassed PI3K-Akt signaling pathway, Lipid and atherosclerosis, MAPK signaling pathways, and Endocrine resistance signaling pathways. Moreover, molecular docking and MD have been shown to reveal a good binding capacity between active compounds and screening targets.
Conclusion: This study systematically investigates the multi-target mechanisms of YKD in the treatment of PCOS, with preliminary verification provided through molecular docking and MD. The findings offer compelling evidence supporting the efficacy of YKD in treating PCOS.
1 Introduction
Polycystic ovary syndrome (PCOS) is a common condition affecting women of reproductive age. It is characterized by a range of abnormalities that impact multiple systems, including the reproductive, metabolic, and psychological components, throughout the entire life cycle (1, 2). Epidemiologic surveys estimate that the global prevalence of PCOS is approximately 8–13% (3). In China, the prevalence is about 5.6% and shows a rising trend annually (4). The clinical manifestations of PCOS are highly variable, but are primarily marked by hyperandrogenism, polycystic ovarian changes, and irregular menstruation or amenorrhea. These manifestations are frequently accompanied by varying degrees of insulin resistance (IR), metabolic issues such as obesity, and mood disorders, including anxiety and depression (5, 6). The symptoms of PCOS not only cause physical and psychological distress for patients but also present significant challenges to their fertility. Moreover, they impose a considerable burden on families and society.
Currently, the management of PCOS primarily emphasizes symptomatic relief through exercise, dietary changes, ovulation-inducing medications, and assisted reproductive technologies (7). Despite these efforts, challenges remain, including extended treatment durations, multiple complications, and less-than-ideal therapeutic outcomes (8). Research has shown that oral contraceptives can have adverse effects such as hormone disruption, vertigo, migraines, and irregular blood pressure and blood sugar levels (9). Moreover, metformin, often prescribed for PCOS, may cause gastrointestinal issues including abdominal pain, diarrhea, nausea, and vomiting (10). It may also result in elevated aminotransferase levels, which can indicate potential liver function problems. Long-term use of these medications can lead to endocrine disruption, osteoporosis, and liver dysfunction, as well as emotional and mental health issues that may negatively impact patient adherence (11). Consequently, research into PCOS is crucial for its medical, economic, and social implications, as it plays a significant role in improving patient well-being and advancing societal progress.
Traditional Chinese Medicine (TCM) offers distinct advantages in the treatment of PCOS, including holistic management, evidence-based approaches, safe medications, and individualized treatment plans (12, 13). Yikang Decoction (YKD) is an empirical formula created by Sun Yun, a renowned Ouyue doctor from Wenzhou City, China. The formula comprises 13 herbs: Radix Lithospermi (Zicao, ZC), Radix Bupleuri (Chaihu, CH), Radix Rubiae (Qiancao, QC), Radix Salviae liguliobae (Danshen, DS), Cortex Moutan (Mudanpi, MDP), Folium Eriobotryae (Pipaye, PPY), Gentiana scabra var. buesgeri (Longdan, LD), Radix Cyathulae (Chuanniuxi, CNX), Rhizoma Cyperi Preparata (Cuxiangfu, CXF), Rhizoma Dioscoreae (Shanyao, SY), Fructus Gardeniae (Zhizi, ZZ), Radix Puerariae (Gegen, GG), and Semen Cuscutae (Tusizi, TSZ). This formulation is known for its effects of clearing liver fire, activating blood circulation, removing blood stasis, and dredging the veins. Our research group has found that YKD can effectively improve patients’ hormone levels and clinical symptoms (14, 15). However, the precise biomolecular pharmacological effects of this compound remain unknown due to its complex chemical nature.
In recent years, the rapid advancement of bioinformatics has introduced significant ideas for systematically studying the mechanisms of action in TCM compounding, garnering widespread attention and application. Network pharmacology can be employed in researching TCM compounding to reveal the interactions among various drug components and their associations with targets in the human body, thereby contributing to the comprehensive understanding of the mechanisms of action of Chinese medicines (16). This approach aligns perfectly with the principle of holistic treatment in TCM, highlighting the complex nature of multi-component, multi-target, and multi-pathway interactions inherent in compound prescriptions. Molecular docking technology facilitates the analysis of interactions and binding between active components and target proteins, aiding in the understanding of their effectiveness and potential mechanisms of action (17). Similarly, MD simulations can model binding interactions and study the movement behavior and structural changes of active ingredients and target proteins over different time scales, further elucidating the pharmacological actions and properties of these ingredients (18). This study utilizes a network pharmacology approach to elucidate the mechanism of action of YKD on PCOS. It is validated through molecular docking and MD simulations, offering new research directions and insights for the clinical treatment of PCOS.
2 Materials and methods
Figure 1 illustrates the process employed in this research study.
2.1 Data download and organization from GEO
The gene expression profile dataset GSE80432 was retrieved from the GEO database, which is affiliated with the National Center for Biotechnology Information. This dataset, obtained using ‘Polycystic Ovary Syndrome’ as the search keyword, includes gene chip data from a total of 16 samples.
2.2 Data analysis using Perl and R languages
Background correction and matrix data normalization were carried out based on the characteristics of the data samples. Gene differential analysis of the microarray data was performed using the limma R package. The filtering criteria for identifying significantly different genes were set with a p-value threshold of <0.05, |logFC| ≥ 0.5 (FC is fold change), and the differentially expressed genes (DEGs) were identified. The ggplot2 package was applied to create a volcano plot of the microarray data, while the pheatmap package was used to draw the heat map of the gene chip.
2.3 Active ingredients and targets prediction
Leveraging the TCMSP database,1 effective components and their target proteins were identified by applying the criteria of Oral Bioavailability (OB) ≥ 30% and Drug-Likeness (DL) ≥ 0.18.
2.4 Construction of the “YKD-active ingredients-targets of action” network diagram
The active ingredients and targets information were organized in Excel to create the network and type files. This data was subsequently imported into Cytoscape 3.10.1 to construct the ‘YKD-active ingredients - targets network and analyze its properties.
2.5 Collection of PCOS-related targets
Keyword searches for “polycystic ovary syndrome” were conducted across several databases, including GeneCards,2 OMIM,3 PharmGKB,4 and DrugBank.5 The resulting data were then consolidated and duplicates were removed. We subsequently integrated these results with differential expression genes from GEO, further eliminating any duplicate entries to create a comprehensive inventory of PCOS.
2.6 Mapping YKD-PCOS protein interactions and core targets screening
A comparison was conducted between the target genes associated with YKD and PCOS using the Venn platform, which produced a Venn diagram. The STRING database6 was utilized to construct the Protein–Protein Interaction (PPI) network for the intersecting genes, focusing on the “Homo sapiens” species with a confidence level set at 0.400. This network was then imported into Cytoscape version 3.10.1. Within Cytoscape, the CytoNCA module was employed to analyze six metrics: Betweenness Centrality (BC), Closeness Centrality (CC), Degree Centrality (DC), Eigenvector Centrality (EC), Local Average Connectivity (LAC), and Network Centrality (NC). Genes that surpassed the median value for these metrics were selected, and this subset was used to construct the core target network.
2.7 Enrichment analysis
The overlapping genes of YKD-PCOS were analyzed using GO and KEGG enrichment analyses. These analyses were performed using the “Bioconductor” package in R software and the bioinformatics platform.7 The results were visualized through bar charts and bubble charts.
2.8 Construction of the YKD active ingredients-intersecting gene-PCOS-KEGG pathways network map
The network map illustrating the interaction between YKD active ingredients and gene-PCOS-KEGG pathways was constructed using Cytoscape 3.10.1. Network properties were analyzed using the “Network Analyzer” module. Genes with scores exceeding the mean were selected as key targets based on degree sorting.
2.9 Molecular docking
The mol2 files for the primary active ingredients were sourced from PubChem. Core target proteins were screened using the UniProt8 and PDB9 databases. Molecular docking was conducted with AutoDock Vina, and the results were visualized using PyMoL. Heat mapping was performed using GraphPad Prism 10.1.2.
2.10 Molecular dynamics simulations
Molecular dynamics simulations of the complexes were conducted using the GROMACS 2020.3 software suite, with the Amber 99SB force field and the TIP3P water model. The simulation spanned 100,000 ps. Following docking, hydrogenation of the small molecule ligands was performed using Avogadro 1.2.0. The hydrogenated ligands were then uploaded to the CGenFF website to generate a topology file. This topology file was combined with the protein topology file to create the comprehensive topology files for the system. A dodecahedral box, utilizing the Simple Point Charge (SPC) water model, was used to enclose the complex. The distance between the complex and the boundary of the enclosure was set to 1 nm. Sodium or chloride ions were added to neutralize the system’s charge. After energy minimization, the system was equilibrated at a temperature of 310 K and a pressure of 1 bar. A subsequent simulation lasting 100 ns was performed. The trajectory files were analyzed to determine system parameters such as Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF), Radius of Gyration (Rg), Solvent Accessible Surface Area (SASA), and Principal Component Analysis (PCA). The results were visualized using Origin 2021.
3 Results
3.1 Screening of DEGs
Following data preprocessing, as illustrated in the volcano plot (Figure 2A), we identified 155 differentially expressed genes in PCOS tissues, with 58 genes showing down-regulation and 97 genes showing up-regulation. The top 20 genes with the most significant up-regulation and down-regulation were selected and analyzed (Figure 2B).
3.2 Active ingredients and action targets of YKD
Through a search of the TCMSP database and integration with the UniProt database to convert standard gene symbols, we identified a total of 139 active ingredients. These ingredients were associated with 315 target points of action, underscoring the molecular complexity of the formula and its potential for multiple mechanisms of action. The network diagram illustrating the “active ingredients - targets” relationships, constructed based on these data, is presented in Figure 3.
3.3 PCOS disease targets
Through systematic mining of GeneCards, OMIM, PharmGKB, DrugBank, and other databases, and after removing duplicates, a total of 1,929 disease targets were identified (Figure 4A). Comparing these targets with those associated with the sexual components of the group formula revealed 165 common genes (Figure 4B).
Figure 4. (A) Venn diagram of “PCOS” targets; (B) Venn diagram of “YKD-PCOS” targets; (C) PPI network diagram of “YKD-PCOS” targets; (D) Results of key targets screening.
3.4 Constructing a PPI network and identifying crucial targets through screening
The PPI network graph was constructed using the STRING platform, incorporating 165 intersecting genes from YKD and PCOS. A confidence level of 0.400 was specified, and free nodes were eliminated. Online analysis revealed that the network comprised 165 nodes and 1,611 edges, with an enrichment p-value of <1.0e-16, an average node degree of 51.9, and an average local clustering coefficient of 0.675. The PPI results were subsequently imported into Cytoscape 3.10.1 for topological analysis, and the network graph was plotted based on the degree values (Figure 4C). The analysis showed that the mean degree centrality of this network is 29.11.
The CytoNCA module employed six indicators BC, CC, DC, EC, LAC, and NC-to identify key genes. After three rounds of optimization, genes with values surpassing the median were selected as key genes. Ultimately, 10 essential genes were identified: TP53, VEGFA, MAPK14, CREB1, ESR1, STAT3, CCND1, AKT1, NOS3, and EGFR. These genes represent 6.06% of the intersecting genes, as shown in Figure 4D.
3.5 GO and KEGG enrichment analysis
The 165 intersecting genes were analyzed for GO enrichment with significance thresholds set at p < 0.05 and q < 0.05. This analysis identified 3,002 potential therapeutic targets for YKD, categorized into 2,708 entries for Biological Processes (BP), 85 entries for Cellular Components (CC), and 209 entries for Molecular Functions (MF). The BPs associated with these potential therapeutic targets included responses to xenobiotic stimuli, epithelial cell proliferation, responses to nutrient levels, and responses to peptide hormones. The GO enrichment analysis prioritized the top 10 entries based on their q-values, ranging from smallest to largest, with the results summarized in Figures 5A,B.
Figure 5. Results of GO and KEGG analysis. (A) Bar plot of GO enrichment analysis; (B) Bubble plot of GO enrichment analysis; (C) Bubble plot of KEGG enrichment analysis.
KEGG pathway analysis was conducted employing R and the bioinformatics platform available (see Footnote 7). A total of 186 pathways were identified as enriched with common targets. These pathways were then ranked by their q-values, from smallest to largest. From this ranking, 30 pathways were selected for the creation of a KEGG pathway bubble map, as shown in Table 1. In this map, the size of each bubble represents the number of associated genes, and the color intensity reflects the level of significance in gene enrichment. The primary signaling pathways identified include the PI3K-Akt signaling pathway, the lipid and atherosclerosis signaling pathway, the MAPK signaling pathway, and endocrine resistance, as depicted in Figure 5C.
3.6 Construction of the YKD-PCOS-KEGG pathways network map
The active ingredients of YKD, along with the intersecting targets, the 10 most central targets of YKD-PCOS, and 30 KEGG pathways, were used to construct a network diagram with Cytoscape 3.10.1 (Figure 6).
Figure 6. Network diagram of “YKD-PCOS-KEGG.” Pathway. Note: orange rectangular nodes represent drugs, rectangular nodes represent active ingredients, purple octagonal nodes represent KEGG pathways, blue circular nodes represent core targets and green rectangular nodes represent intersecting targets.
3.7 Molecular docking results
After identifying the core genes using Cytoscape, target proteins (6P8G, 2IOG, 2YIX, 6TLC, 8F0W) were retrieved from the PDB database based on the corresponding genes associated with the active ingredients. The core active ingredients—Chryseriol, Cryptotanshinone, Diosgenin, Formononetin, and Tanshinone IIA—were used as ligand molecules. The spatial coordinates for docking, as well as the details of the targets and compounds, are provided in Table 2. The ligands and receptors were docked separately, and the binding energy values are illustrated in Figure 7A. The results indicated that Cryptotanshinone_2IOG had the lowest binding energy of −10.2 kcal·mol−1. Based on the principle of lowest energy, the top candidate from each group was selected for visualization, analysis, and MD. Figures 7B–F present the 3D and 2D plots of the core components and their corresponding core targets visualized after molecular docking.
Figure 7. (A) Heat map of molecular docking score. Docking results of Binding energy(kcal·mol−1) of key targets and active compounds of herbs in 2D and 3D images: (B) Diosgenin_6P8G; (C) Cryptotanshinone_2IOG; (D) Formononetin_2YIX; (E) Diosgenin_6TLC; (F) Tanshinone IIA_8F0W.
3.8 MD results
Based on the molecular docking results, the following protein-ligand complex systems were analyzed using 100 ns MD simulations: Diosgenin_6P8G, Cryptotanshinone_2IOG, Formononetin_2YIX, Diosgenin_6TLC, and Tanshinone IIA_8F0W. After performing the 100 ns MD simulations, trajectory files were generated, and the stability of these complex systems was assessed through various analyses, including RMSD, RMSF, Rg, SASA, H-bond changes, and PCA. The detailed results of these analyses are presented in Figure 8.
Figure 8. (A) RMSD; (B) RMSF; (C) Rg; (D) SASA; (E) NHB; (F) PCA plots during of molecular dynamics simulation.
4 Discussion
PCOS is a prevalent gynecologic endocrine disorder that presents significant challenges for both families and society. With declining fertility rates and a trend toward delayed marriage, PCOS is becoming a major global public health concern (19). The disorder can lead to insulin resistance, hyperlipidemia, and metabolic syndrome (MS), which in turn increase the risk of glucose intolerance, type 2 diabetes mellitus, cardiovascular disease, and neoplasms in women (20, 21). Additionally, PCOS is associated with dyslipidemia, as well as emotional and psychosexual disorders, posing a serious threat to overall health and well-being. TCM, which has undergone extensive development over thousands of years, has had a profound and far-reaching impact both within China and globally. TCM offers distinct advantages in managing various illnesses, particularly chronic and challenging conditions (22, 23). Currently, our research group has not conducted in-depth basic research beyond evaluating clinical efficacy and animal experiments. Given the multi-component and multi-target nature of TCM, a comprehensive research approach is required to analyze the mechanisms of YKD in treating PCOS. Consequently, we employed network pharmacology combined with molecular docking and MD simulations to explore the mechanisms of YKD’s potential active components in the treatment of PCOS.
Following a network pharmacological analysis, this study identified 139 probable active components and proposed that YKD may exert its therapeutic effects on PCOS through the actions of Chryseriol, Cryptotanshinone, Diosgenin, Formononetin, and Tanshinone IIA. Chryseriol, sourced from Rhizoma Cyperi Preparata, possesses a range of biological activities, including antioxidant, anti-inflammatory, and hormone-regulating effects, making it a valuable candidate for studying and treating endocrine system diseases (24). Notably, Chryseriol has demonstrated phytoestrogen-like activity, which has garnered interest in the context of gynecological disorders such as menstrual irregularities and PCOS (25). Additionally, Cryptotanshinone and Tanshinone IIA are derived from Radix Salviae liguliobae.
Cryptotanshinone, a diterpene quinone compound, has been shown to modulate estrogen receptors and potentially influence estrogen activity (26). This has made it a focal point in studies related to breast cancer, endometriosis, and other estrogen-related diseases. Besides, research indicates that Cryptotanshinone may impact insulin signaling and glucose metabolism, generating interest in its potential applications for diabetes and metabolic disease research (27). Tanshinone IIA inhibits estrogen synthesis and promotes luteinizing hormone production, which may help restore the balance between estrogen and luteinizing hormone, thereby improving the PCOS (28). Furthermore, Tanshinone IIA enhances ovarian function and promotes oocyte development in mice, contributing to the restoration of normal ovulatory cycles (29). Diosgenin, derived from Rhizoma Dioscoreae, has been found to exhibit phytoestrogen-like effects in menopausal syndrome and estrogen-related disorders. It may help relieve symptoms of menopausal syndrome and ameliorate issues associated with estrogen deficiency (30). Moreover, Diosgenin has been investigated for its potential in treating PCOS, as it may modulate ovarian function and androgen levels, thereby improving symptoms and signs in PCOS patients (31). Formononetin, derived from Pueraria Mirifica, has demonstrated in some animal and in vitro studies the ability to mimic the biological activity of estrogen, modulate the insulin signaling pathway, ameliorate insulin resistance, and regulate the secretion of ovarian hormones. This suggests that Formononetin plays a modulating role in PCOS patients with abnormal estrogen levels (32, 33). Notably, Formononetin possesses antioxidant, anti-inflammatory, and anti-metabolic syndrome properties, which may positively affect symptoms and metabolic abnormalities in PCOS patients (34). The study of the PPI interaction network and CytoNCA findings revealed that the primary targets of YKD for the treatment of PCOS are CCND1, ESR1, MAPK14, STAT3, and TP53. Notably, the CCND1 protein plays a critical role in regulating biological processes such as cell growth, proliferation, and differentiation. It is also involved in regulating cell regeneration and repair processes in normal tissues, which is essential for maintaining the physiological functions of these tissues and organs (35). In patients with PCOS, ovarian follicle growth is often impaired, primarily due to an imbalance between the proliferation and apoptosis of ovarian granulosa cells. CCND1 influences the normal growth and maturation of follicles by controlling cell proliferation and differentiation. It promotes cell proliferation by interacting with CDK4/6, which facilitates the phosphorylation of the Rb protein and releases the E2F transcription factor, thereby enabling cells to transition from the G1 phase to the S phase (36). Estrogen receptor alpha (ERα), encoded by the ESR1 gene, plays a crucial regulatory role in ovarian tissues. It is involved in regulating growth, differentiation, and hormone sensitivity. Research on PCOS suggests that the ESR1 gene and the estrogen receptor it encodes may significantly contribute to the pathogenesis of PCOS (37). Polymorphisms in the ESR1 gene and abnormal expression of the estrogen receptor have been associated with the pathophysiological characteristics of PCOS patients. Certain polymorphic variants in the ESR1 gene have been identified in some PCOS patients, potentially affecting estrogen receptor function and estrogen signaling (38). Additionally, some studies have indicated that the expression levels of estrogen receptors in the ovarian tissues of PCOS patients may be altered, potentially contributing to impaired follicular development and abnormal ovulation. Beyond directly affecting the biological functions of ovarian tissues, the ESR1 gene and its encoded estrogen receptor may also regulate the balance of estrogen and other hormones by influencing endocrine system function, thereby impacting the endocrine metabolic status of PCOS patients (39). MAPK14, also known as p38 MAP kinase, is a member of the mitogen-activated protein kinase (MAPK) family, which regulates various biological processes, including cell proliferation, differentiation, apoptosis, and inflammatory responses. Studies related to PCOS suggest that MAPK14 may be involved in the disease’s pathogenesis (40). Research has shown that in PCOS, a common gynecological condition caused by hormonal imbalance and abnormal proliferation of ovarian tissues—there is often increased activity and expression of MAPK14 in the ovarian tissues of affected patients (41). STAT3, a crucial transcription factor, plays a significant role in regulating various cellular processes, including cell proliferation, differentiation, apoptosis, and inflammatory responses (42). The study revealed a notable increase in both the activity and expression levels of STAT3 in the ovarian tissues of rats with PCOS (43, 44). This increase may lead to excessive cell proliferation and impaired follicular development, resulting in disrupted ovulation and the formation of cystic ovaries. Overexpression of STAT3 contributes to ovarian dysfunction and insulin resistance in mice, further exacerbating pathological changes (45). The primary functions of TP53 include regulating the cell cycle, inducing apoptosis, and maintaining genome stability (46). Research has shown that TP53 expression is significantly reduced in the ovarian tissues of individuals with PCOS. This reduction may lead to increased cell growth, decreased apoptosis, and impaired formation of ovarian follicles (47). Lower TP53 activity can disrupt the cell cycle and diminish the repair of DNA damage, resulting in abnormal ovarian tissue growth and the development of polycystic ovaries. In addition, TP53 is involved in the pathogenesis of PCOS by affecting the inflammatory response and metabolic regulation of ovarian tissues, among other factors (48).
Furthermore, to provide a clearer understanding of the current research on these components and their key mechanisms of action, we have summarized the relationships between five active compounds and their respective targets, based on findings from previously published literature.
Previous research demonstrates that Chrysoeriol effectively alleviates skin inflammation by modulating STAT3 transcription, leading to reduced mRNA levels of the pro-inflammatory cytokines IL-6, IL-1β, and TNF-α (49). Chrysoeriol also inhibits the excessive proliferation of rheumatoid arthritis fibroblast-like synoviocytes (RA-FLS) by preventing the activation and phosphorylation of STAT3 (Tyr705), lowering STAT3 nuclear levels, and downregulating the expression of Bcl-2 and Mcl-1, both of which are regulated by STAT3 (50). Additionally, studies have indicated that the combined administration of Cryptotanshinone and Temozolomide has a significantly stronger inhibitory effect on cell proliferation than when each drug is used individually. Western blot analysis confirmed that the combination therapy markedly reduced STAT3 expression in glioblastoma multiforme cells compared to Temozolomide alone (51). Cryptotanshinone exhibits lethal effects on ER-positive breast cancer cells by modulating their proliferation and migration, primarily through the inhibition of ESR1 (52). Additionally, diosgenin directly inhibits the phosphorylation of STAT3 and its transcriptional regulation of NPC1L1 expression, leading to a significant reduction in NPC1L1 levels. This downregulation hinders intestinal cholesterol absorption, thereby preventing the formation of cholesterol gallstones (53). In food allergy responses, elevated levels of IgE and hyperresponsive mast cells and basophils are critical factors. Validation through qPCR demonstrates that Formononetin suppresses IgE-mediated mast cell activation by upregulating TP53 gene expression and downregulating the expression of STAT3 and CCND1, ultimately reducing the severity of food allergies (54). The proposed anti-allodynic and anti-neuroinflammatory effects of sodium tanshinone IIA sulfonate (STS), a derivative of tanshinone IIA, have been observed in neuropathic pain models, likely involving the miR-125b-5p/STAT3 axis (55). In conclusion, the clear correlation between these components and their specific targets underscores the significance of our study on YKD, warranting further investigation.
The results of GO analysis for relevant targets indicate that YKD demonstrates efficacy in PCOS by influencing various biological processes, cellular components, and molecular functions. Among these aspects, biological processes showed the highest level of enrichment, representing 90.21% of the GO analysis results. Key biological processes identified include the response to xenobiotic stimuli, epithelial cell proliferation, response to nutrient levels, and response to peptide hormones. GO analysis is instrumental in elucidating the molecular mechanisms underlying PCOS and enhancing our understanding of its onset and progression. By examining significantly enriched GO terms, researchers can pinpoint critical biological processes, molecular functions, and cellular components involved in PCOS. This information supports further functional studies, facilitates the identification of potential therapeutic targets, and aids in the development of novel treatment strategies.
The KEGG analysis indicates that YKD treatment for PCOS may influence several signaling pathways, including the PI3K-Akt signaling pathway, lipid and atherosclerosis signaling pathway, MAPK signaling pathway, and endocrine resistance pathway. In PCOS patients, the PI3K-Akt signaling pathway may be abnormally activated or dysfunctional, leading to metabolic disorders, ovarian hormonal imbalances, and the formation of ovarian polycysts (56). Specifically, PI3K (phosphatidylinositol 3-kinase) is the initiating enzyme of the PI3K-Akt signaling pathway, and its primary role is to transmit external growth factor signals into the cell. In patients with PCOS, the expression and activity of PI3K are increased, leading to over-activation of the PI3K-Akt signaling pathway. Akt (protein kinase B), a key component of this pathway, regulates several downstream signaling molecules, including mTOR, GSK3β, and FOXO, after phosphorylation. In PCOS patients, elevated phosphorylation levels of Akt exacerbate the activity of the PI3K-Akt signaling pathway. Excessive activation of this pathway can impact various biological processes within ovarian tissue, such as follicular development, ovarian hormone synthesis and secretion, and metabolic regulation (57). These effects contribute to features such as ovarian polycysticity, abnormal estrogen levels, and insulin resistance in PCOS patients. PCOS is a common endocrine disorder frequently associated with metabolic abnormalities, including insulin resistance, high cholesterol, and hypertriglyceridemia. These metabolic issues can increase the risk of atherosclerosis. There is a complex relationship between lipid metabolism, atherosclerotic signaling pathways, and PCOS. Abnormalities in lipid metabolism are common among PCOS patients (58). Due to insulin resistance, patients with PCOS often exhibit dyslipidemia, such as elevated cholesterol, high triacylglycerols, low-density lipoprotein (LDL), and very low-density lipoprotein (VLDL) (59). These lipid imbalances can disrupt metabolism, leading to lipid accumulation in blood vessels and potential atherosclerotic plaque formation. Atherosclerosis is a chronic inflammatory disease involving multiple inflammatory cells and factors (60). In PCOS patients, chronic low-grade inflammation, associated with metabolic abnormalities and abnormal androgen levels, may contribute to atherosclerosis development (61). Consequently, weight gain and abdominal obesity in PCOS patients can trigger an inflammatory response and insulin resistance, further increasing the risk of atherosclerosis (62). In PCOS patients, the imbalance between estrogen and androgen levels leads to abnormal activation of the MAPK signaling pathway in ovarian and other related cells. This pathway is involved in regulating physiological processes such as follicular development, androgen synthesis, and the proliferation and apoptosis of ovarian cells (63). Consequently, this can lead to abnormal ovarian function and cyst formation, contributing to PCOS development. Moreover, the chronic low-grade inflammatory state in PCOS patients is linked to the MAPK signaling pathway, which plays a significant role in the inflammatory response. Activated MAPK signaling can promote the release of inflammatory factors, exacerbating the inflammatory response and affecting the body’s metabolic and endocrine balance. These changes are closely associated with the pathogenesis of PCOS. Endocrine resistance, a common phenomenon where the body’s response to endocrine therapy weakens or fails, reduces the effectiveness of treatments and complicates disease management (64). PCOS patients often exhibit insulin resistance and elevated androgen levels, which can lead to resistance to certain endocrine therapies (65). In summary, multiple intricate signaling pathways play crucial roles in the pathological progression of PCOS, including those involved in cell survival, proliferation, differentiation, metabolic regulation, inflammatory response, and endocrine regulation. A comprehensive examination and understanding of these pathways will enhance the treatment strategies for PCOS.
Combining the perspectives of semi-flexible docking and energy stabilization, the visualization revealed that Chryseriol binds to CCND1 (PDB ID: 6P8G) via residue LYS-238, Cryptotanshinone binds to ESR1 (PDB ID: 2IOG) via residue TPR-392, Diosgenin binds to MAPK14 (PDB ID: 2YIX) via residue GLU-71, Formononetin binds to STAT3 (PDB ID: 6TLC) via residues GLN-644 and GLY-656, and Tanshinone IIA binds to STAT3 (PDB ID: 6TLC) via residues LEU-1568 and GLY-656. The binding energies for these interactions were − 8.9 kcal·mol−1, −10.1 kcal·mol−1, −9.4 kcal·mol−1, −7.9 kcal·mol−1, and − 8.0 kcal·mol−1, respectively. The binding energies were all <−7.5 kcal·mol−1, indicating strong affinity and stability. These findings suggest that these compounds may play significant roles in YKD’s pharmacological effects. However, the results of molecular docking might not always reflect the most accurate docking scenarios. To address this, we performed MD simulations on the complex systems. Using Gromacs, we conducted MD simulations on the most favorable conformations obtained from molecular docking. Through various analyses, including Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF), Radius of Gyration (Rg), Solvent Accessible Surface Area (SASA), Number of H-bonds (NHB), and Principal Component Analysis (PCA), we obtained further evidence supporting the stability and strong binding of Diosgenin with CCND1 (6P8G), Cryptotanshinone with ESR1 (2IOG), Formononetin with MAPK14 (2YIX), Diosgenin with STAT3 (6TLC), and Tanshinone IIA with TP53 (8F0W).
5 Summary
This study comprehensively investigated the potential of YKD for treating PCOS through network pharmacology, molecular docking, and MD techniques. The results suggest that YKD may exert its therapeutic effects on PCOS through key components such as Chryseriol, Cryptotanshinone, Formononetin, Tanshinone IIA, and Diosgenin. The MD data revealed significant associations between core genes and these components, reinforcing YKD’s therapeutic impact on PCOS. However, further clinical studies and experimental validation are necessary to fully elucidate YKD’s mechanisms in PCOS treatment. Overall, our findings validate the reliability of our molecular docking results and highlight YKD’s potential as a promising treatment for PCOS.
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 author.
Author contributions
JM: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. LG: Data curation, Writing – review & editing. XL: Data curation, Writing – review & editing. WC: Investigation, Writing – review & editing. LC: Formal analysis, Writing – review & editing. MC: Data curation, Validation, Writing – review & editing. YS: Funding acquisition, Supervision, Validation, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was funded by Wenzhou major science and technology projects (No. ZY2022014), and Basic Public Welfare Research Program of Wenzhou City (No. Y20240927).
Acknowledgments
The authors thank the members of their laboratory and their collaborators for their research work.
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.
Footnotes
1. ^http://lsp.nwu.edu.cn/tcmsp
2. ^https://www.genecards.org/
References
1. McCartney, CR, and Marshall, JCClinical Practice. Polycystic ovary syndrome. N Engl J Med. (2016) 375:54–64. doi: 10.1056/NEJMcp1514916
2. Wang, Z, Jukic, AMZ, Baird, DD, Wilcox, AJ, Li, H, Curry, CL, et al. Irregular cycles, ovulatory disorders, and Cardiometabolic conditions in a US-based digital cohort. JAMA Netw Open. (2024) 7:e249657. doi: 10.1001/jamanetworkopen.2024.9657
3. Liu, Q, Zhu, Z, Kraft, P, Deng, Q, Stener-Victorin, E, and Jiang, X. Genomic correlation, shared loci, and causal relationship between obesity and polycystic ovary syndrome: a large-scale genome-wide cross-trait analysis. BMC Med. (2022) 20:66. doi: 10.1186/s12916-022-02238-y
4. Dobbin, ZC, and Landen, CN. The importance of the PI3K/AKT/MTOR pathway in the progression of ovarian cancer. Int J Mol Sci. (2013) 14:8213–27. doi: 10.3390/ijms14048213
5. Joham, AE, Norman, RJ, Stener-Victorin, E, Legro, RS, Franks, S, Moran, LJ, et al. Polycystic ovary syndrome. Lancet Diabetes Endocrinol. (2022) 10:668–80. doi: 10.1016/s2213-8587(22)00163-2
6. Novakovic, S, Jakovljevic, V, Jovic, N, Andric, K, Milinkovic, M, Anicic, T, et al. Exploring the Antioxidative effects of ginger and cinnamon: a comprehensive review of evidence and molecular mechanisms involved in polycystic ovary syndrome (PCOS) and other oxidative stress-related disorders. Antioxidants. (2024) 13:392–11. doi: 10.3390/antiox13040392
7. Legro, RS, Arslanian, SA, Ehrmann, DA, Hoeger, KM, Murad, MH, Pasquali, R, et al. Diagnosis and treatment of polycystic ovary syndrome: an Endocrine Society clinical practice guideline. J Clin Endocrinol Metab. (2013) 98:4565–92. doi: 10.1210/jc.2013-2350
8. Stener-Victorin, E, and Deng, Q. Epigenetic inheritance of polycystic ovary syndrome - challenges and opportunities for treatment. Nat Rev Endocrinol. (2021) 17:521–33. doi: 10.1038/s41574-021-00517-x
9. Vrbíková, J, and Cibula, D. Combined oral contraceptives in the treatment of polycystic ovary syndrome. Hum Reprod Update. (2005) 11:277–91. doi: 10.1093/humupd/dmi005
10. Miralles-Linares, F, Puerta-Fernandez, S, Bernal-Lopez, MR, Tinahones, FJ, Andrade, RJ, and Gomez-Huelgas, R. Metformin-induced hepatotoxicity. Diabetes Care. (2012) 35:e21. doi: 10.2337/dc11-2306
11. Ee, C, Pirotta, S, Mousa, A, Moran, L, and Lim, S. Providing lifestyle advice to women with PCOS: an overview of practical issues affecting success. BMC Endocr Disord. (2021) 21:234. doi: 10.1186/s12902-021-00890-8
12. Ong, M, Peng, J, Jin, X, and Qu, X. Chinese herbal medicine for the optimal management of polycystic ovary syndrome. Am J Chin Med. (2017) 45:405–22. doi: 10.1142/s0192415x17500252
13. Zhang, Q, Ren, J, Wang, F, Li, M, Pan, M, Zhang, H, et al. Chinese herbal medicine alleviates the pathogenesis of polycystic ovary syndrome by improving oxidative stress and glucose metabolism via mitochondrial Sirtuin 3 signaling. Phytomedicine. (2023) 109:154556. doi: 10.1016/j.phymed.2022.154556
14. Wang, W, and Ma, D. Intervention role of YiKang decoction in the treatment of PCOS rats model. J Liaoning Univ Tradit Chin Med. (2010) 12:184–6. doi: 10.13194/j.jlunivtcm.2010.03.186.wangw.07
15. Sun, Y . Theoretical and clinical study on treatment of adolescent polycystic ovary syndrome by Yi Kang decoction. J Zhejiang Univ Tradit Chin Med. (2018) 42:24–8. doi: 10.16466/j.issn1005-5509.2018.01.004
17. Pinzi, L, and Rastelli, G. Molecular docking: shifting paradigms in drug discovery. Int J Mol Sci. (2019) 20:4331–354. doi: 10.3390/ijms20184331
18. Chen, Y, Tang, H, Li, H, Yin, Y, Song, W, Guo, H, et al. Molecular-level insight into the behavior of metal cations and organic matter during the aggregation of polystyrene nanoplastics. J Hazard Mater. (2024) 473:134665. doi: 10.1016/j.jhazmat.2024.134665
19. Stener-Victorin, E, Teede, H, Norman, RJ, Legro, R, Goodarzi, MO, Dokras, A, et al. Polycystic ovary syndrome. Nat Rev Dis Primers. (2024) 10:27. doi: 10.1038/s41572-024-00511-3
20. Zhu, T, and Goodarzi, MO. Causes and consequences of polycystic ovary syndrome: insights from Mendelian randomization. J Clin Endocrinol Metab. (2022) 107:e899–911. doi: 10.1210/clinem/dgab757
21. Ho, B, Ryback, D, Benson, B, Mason, CN, Torres, PJ, Quinn, RA, et al. Gut metabolites are more predictive of disease and cohoused states than gut bacterial features in a polycystic ovary syndrome-like mouse model. mSystems. (2021) 6:e0114920. doi: 10.1128/mSystems.01149-20
22. Huang, K, Zhang, P, Zhang, Z, Youn, JY, Wang, C, Zhang, H, et al. Traditional Chinese medicine (TCM) in the treatment of COVID-19 and other viral infections: efficacies and mechanisms. Pharmacol Ther. (2021) 225:107843. doi: 10.1016/j.pharmthera.2021.107843
23. Zhai, Z, Su, PW, Ma, LY, Yang, H, Wang, T, Fei, ZG, et al. Progress on traditional Chinese medicine in treatment of ischemic stroke via the gut-brain axis. Biomed Pharmacother. (2023) 157:114056. doi: 10.1016/j.biopha.2022.114056
24. Aboulaghras, S, Sahib, N, Bakrim, S, Benali, T, Charfi, S, Guaouguaou, FE, et al. Health benefits and pharmacological aspects of Chrysoeriol. Pharmaceuticals (Basel). (2022) 15:973–99. doi: 10.3390/ph15080973
25. Aćimović, M, Jeremić, K, Salaj, N, Gavarić, N, Kiprovski, B, Sikora, V, et al. Marrubium vulgare L.: a phytochemical and pharmacological overview. Molecules. (2020) 25:2898–922. doi: 10.3390/molecules25122898
26. Li, S, Wang, H, Hong, L, Liu, W, Huang, F, Wang, J, et al. Cryptotanshinone inhibits breast cancer cell growth by suppressing estrogen receptor signaling. Cancer Biol Ther. (2015) 16:176–84. doi: 10.4161/15384047.2014.962960
27. Tian, J, Zhao, Y, Liu, Y, Liu, Y, Chen, K, and Lyu, S. Roles and mechanisms of herbal medicine for diabetic cardiomyopathy: current status and perspective. Oxidative Med Cell Longev. (2017) 2017:8214541. doi: 10.1155/2017/8214541
28. Jin, J, Hu, QY, Xu, WW, Zhu, WJ, Liu, B, Liu, J, et al. Tanshinone IIA attenuates estradiol-induced polycystic ovarian syndrome in mice by ameliorating FSHR expression in the ovary. Exp Ther Med. (2019) 17:3501–8. doi: 10.3892/etm.2019.7352
29. Bai, L, He, GZ, Gao, CH, Yang, H, Li, MX, Huang, YL, et al. Tanshinone IIA enhances the ovarian reserve and attenuates ovarian oxidative stress in aged mice. Vet Med Sci. (2022) 8:1617–25. doi: 10.1002/vms3.811
30. Kronenberg, F, and Fugh-Berman, A. Complementary and alternative medicine for menopausal symptoms: a review of randomized, controlled trials. Ann Intern Med. (2002) 137:805–13. doi: 10.7326/0003-4819-137-10-200211190-00009
31. Sirotkin, AVJIJRMG . Action of diosgenin and diosgenin-containing plants on health and female reproduction. Int J Reprod Contracept Obstet Gynecol. (2021) 7:018–24.
32. Waris, G, Ahmed, A, Rukhsar, A, Zahra, FT, and Malik, A. Zahra F-t and Malik a. Formononetin as promising therapeutic intervention for restoring ovarian, uterine, and Hepato-renal functions in Letrozole-Induced,Polycystic ovarian syndrome Sprague Dawley rats. JPTCP. (2023) 30:727–35. doi: 10.53555/jptcp.v30i19.3749
33. Joshi, M, Shankar, R, Pathak, K, and Yadav, R. Polycystic ovarian syndrome: a review covering phytoconstituents for its outstrip management. Modern Chinese Med. (2021) 1:100011. doi: 10.1016/j.prmcm.2021.100011
34. Machado Dutra, J, Espitia, P, and Andrade Batista, R. Formononetin: biological effects and uses–a review. Food Chem. (2021) 359:129975. doi: 10.1016/j.foodchem.2021.129975
35. Asghar, U, Witkiewicz, AK, Turner, NC, and Knudsen, ES. The history and future of targeting cyclin-dependent kinases in cancer therapy. Nat Rev Drug Discov. (2015) 14:130–46. doi: 10.1038/nrd4504
36. Herrera-Abreu, MT, Palafox, M, Asghar, U, Rivas, MA, Cutts, RJ, Garcia-Murillas, I, et al. Early adaptation and acquired resistance to CDK4/6 inhibition in estrogen receptor-positive breast Cancer. Cancer Res. (2016) 76:2301–13. doi: 10.1158/0008-5472.Can-15-0728
37. Liu, F, Tian, L, Tan, J, Li, Z, Qin, H, Xu, D, et al. Identification of a novel ESR1 mutation in a Chinese PCOS woman with estrogen insensitivity in IVF treatment. Reprod Biol Endocrinol. (2022) 20:157. doi: 10.1186/s12958-022-01029-7
38. Jiao, X, Chen, W, Zhang, J, Wang, W, Song, J, Chen, D, et al. Variant alleles of the ESR1, PPARG, HMGA2, and MTHFR genes are associated with polycystic ovary syndrome risk in a Chinese population: a case-control study. Front Endocrinol (Lausanne). (2018) 9:504. doi: 10.3389/fendo.2018.00504
39. Tang, ZR, Zhang, R, Lian, ZX, Deng, SL, and Yu, K. Estrogen-receptor expression and function in female reproductive disease. Cells. (2019) 8:1–15. doi: 10.3390/cells8101123
40. Zhu, W, Han, B, Fan, M, Wang, N, Wang, H, Zhu, H, et al. Oxidative stress increases the 17,20-lyase-catalyzing activity of adrenal P450c17 through p38α in the development of hyperandrogenism. Mol Cell Endocrinol. (2019) 484:25–33. doi: 10.1016/j.mce.2019.01.020
41. Mao, Z, Li, T, Zhao, H, Qin, Y, Wang, X, and Kang, Y. Identification of epigenetic interactions between microRNA and DNA methylation associated with polycystic ovarian syndrome. J Hum Genet. (2021) 66:123–37. doi: 10.1038/s10038-020-0819-6
42. Zou, S, Tong, Q, Liu, B, Huang, W, Tian, Y, and Fu, X. Targeting STAT3 in Cancer immunotherapy. Mol Cancer. (2020) 19:145. doi: 10.1186/s12943-020-01258-7
43. Wang, H, Feng, X, Wang, T, Pan, J, Zheng, Z, Su, Y, et al. Role and mechanism of the p-JAK2/p-STAT3 signaling pathway in follicular development in PCOS rats. Gen Comp Endocrinol. (2023) 330:114138. doi: 10.1016/j.ygcen.2022.114138
44. Sudhakaran, G, Priya, PS, Jagan, K, Haridevamuthu, B, Meenatchi, R, and Arockiaraj, J. Osteoporosis in polycystic ovary syndrome (PCOS) and involved mechanisms. Life Sci. (2023) 335:122280. doi: 10.1016/j.lfs.2023.122280
45. Hong, G, Wu, H, Ma, ST, and Su, Z. Catechins from oolong tea improve uterine defects by inhibiting STAT3 signaling in polycystic ovary syndrome mice. Chin Med. (2020) 15:125. doi: 10.1186/s13020-020-00405-y
46. Olivier, M, Hollstein, M, and Hainaut, P. TP53 mutations in human cancers: origins, consequences, and clinical use. Cold Spring Harb Perspect Biol. (2010) 2:a001008. doi: 10.1101/cshperspect.a001008
47. Yumiceba, V, López-Cortés, A, Pérez-Villa, A, Yumiseba, I, Guerrero, S, García-Cárdenas, JM, et al. Oncology and pharmacogenomics insights in polycystic ovary syndrome: an integrative analysis. Front. Front Endocrinol. (2020) 11:585130. doi: 10.3389/fendo.2020.585130
48. Vousden, KH, and Prives, C. Blinded by the light: the growing complexity of p53. Cell. (2009) 137:413–31. doi: 10.1016/j.cell.2009.04.037
49. Law, SK, Wu, XX, Jiang, Z, Tong, CWS, Chow, WYL, and Au, DCT. Pharmacological activities of Lonicerae japonicae flos and its derivative-"Chrysoeriol" in skin diseases. Molecules. (2024) 29:3390–408. doi: 10.3390/molecules29091972
50. Wu, JY, Chen, YJ, Fu, XQ, Li, JK, Chou, JY, Yin, CL, et al. Chrysoeriol suppresses hyperproliferation of rheumatoid arthritis fibroblast-like synoviocytes and inhibits JAK2/STAT3 signaling. BMC Complement Med Ther. (2022) 22:73. doi: 10.1186/s12906-022-03553-w
51. Zhu, S, Guo, J, Yu, L, Liu, J, Chen, J, Xin, J, et al. Synergistic effect of cryptotanshinone and temozolomide treatment against human glioblastoma cells. Sci Rep. (2023) 13:21835. doi: 10.1038/s41598-023-48777-z
52. Li, H, Gao, C, Liang, Q, Liu, C, Liu, L, Zhuang, J, et al. Cryptotanshinone is a intervention for ER-positive breast Cancer: an integrated approach to the study of natural product intervention mechanisms. Front Pharmacol. (2021) 11:592109. doi: 10.3389/fphar.2020.592109
53. Shen, W, Shao, W, Wang, Q, Wang, B, Zhao, G, Gu, A, et al. Dietary diosgenin transcriptionally down-regulated intestinal NPC1L1 expression to prevent cholesterol gallstone formation in mice. J Biomed Sci. (2023) 30:44. doi: 10.1186/s12929-023-00933-3
54. Musa, I, Wang, ZZ, Yang, N, and Li, XM. Formononetin inhibits IgE by huPlasma/PBMCs and mast cells/basophil activation via JAK/STAT/PI3-Akt pathways. Front Immunol. (2024) 15:1427563. doi: 10.3389/fimmu.2024.1427563
55. Zeng, J, Gao, WW, Yang, H, Wang, YN, Mei, Y, Liu, TT, et al. Sodium tanshinone IIA sulfonate suppresses microglia polarization and neuroinflammation possibly via regulating miR-125b-5p/STAT3 axis to ameliorate neuropathic pain. Eur J Pharmacol. (2024) 972:176523. doi: 10.1016/j.ejphar.2024.176523
56. Shi, X, Wang, J, Lei, Y, Cong, C, Tan, D, and Zhou, X. Research progress on the PI3K/AKT signaling pathway in gynecological cancer (review). Mol Med Rep. (2019) 19:4529–35. doi: 10.3892/mmr.2019.10121
57. Li, T, Mo, H, Chen, W, Li, L, Xiao, Y, Zhang, J, et al. Role of the PI3K-Akt signaling pathway in the pathogenesis of polycystic ovary syndrome. Reprod Sci. (2017) 24:646–55. doi: 10.1177/1933719116667606
58. Ferrer, MJ, Silva, AF, Abruzzese, GA, Velázquez, ME, and Motta, AB. Lipid metabolism and relevant disorders to female reproductive health. Curr Med Chem. (2021) 28:5625–47. doi: 10.2174/0929867328666210106142912
59. Pazderska, A, Kyaw Tun, T, Phelan, N, McGowan, A, Sherlock, M, Behan, L, et al. In women with PCOS, waist circumference is a better surrogate of glucose and lipid metabolism than disease status per se. Clin Endocrinol. (2018) 88:565–74. doi: 10.1111/cen.13542
60. Lu, H, and Daugherty, A. Atherosclerosis. Arterioscler Thromb Vasc Biol. (2015) 35:485–91. doi: 10.1161/atvbaha.115.305380
61. Rudnicka, E, Suchta, K, Grymowicz, M, Calik-Ksepka, A, Smolarczyk, K, Duszewska, AM, et al. Chronic low grade inflammation in pathogenesis of PCOS. Int J Mol Sci. (2021) 22:3789–801. doi: 10.3390/ijms22073789
62. Ibáñez, L, and de Zegher, F. Adolescent PCOS: a postpubertal central obesity syndrome. Trends Mol Med. (2023) 29:354–63. doi: 10.1016/j.molmed.2023.02.006
63. Hu, M, Zhang, Y, Li, X, Cui, P, Li, J, Brännström, M, et al. Alterations of endometrial epithelial-mesenchymal transition and MAPK signalling components in women with PCOS are partially modulated by metformin in vitro. Mol Hum Reprod. (2020) 26:312–26. doi: 10.1093/molehr/gaaa023
64. Wang, J, Yin, T, and Liu, S. Dysregulation of immune response in PCOS organ system. Front Immunol. (2023) 14:1169232. doi: 10.3389/fimmu.2023.1169232
Keywords: Yikang decoction, polycystic ovary syndrome, network pharmacology, molecular docking, molecular dynamics simulation
Citation: Miao J, Gao L, Liu X, Cai W, Chen L, Chen M and Sun Y (2024) Exploring the therapeutic mechanisms of Yikang decoction in polycystic ovary syndrome: an integration of GEO datasets, network pharmacology, and molecular dynamics simulations. Front. Med. 11:1455964. doi: 10.3389/fmed.2024.1455964
Edited by:
Subhash C. Mandal, Directorate of Drugs Control, Government of West Bengal, IndiaReviewed by:
Prasoon Gupta, Indian Institute of Integrative Medicine (CSIR), IndiaAsis Bala, Ministry of Science and Technology, India
Copyright © 2024 Miao, Gao, Liu, Cai, Chen, Chen and Sun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yun Sun, sunyun7304@outlook.com