Integrative Multi-cohort Transcriptomics and Network Pharmacology Analysis Reveals Key Network Nodes and Potential Drug Clues in PCOS Granulosa Cells

Background Polycystic ovary syndrome (PCOS) is a prevalent endocrine disorder with complex pathophysiology and limited therapeutic options. Identifying key molecular drivers and potential drug candidates is critical for improving clinical outcomes. Methods We integrated multi-cohort transcriptomics (GSE155489, GSE138518, GSE226146) with weighted gene co-expression network analysis (WGCNA), protein-protein interaction (PPI) network analysis, and drug repurposing. Differential expression analysis identified 1,039 DEGs, and WGCNA identified 10 PCOS-associated modules. Intersection of DEGs with module genes yielded 498 core candidate genes, which were subjected to functional enrichment, PPI network analysis, and connectivity map-based drug repurposing (CLUE/LINCS). Candidate drugs were further evaluated by molecular docking and ADMET prediction using a triple intersection strategy (hub genes, high differential expression, drug-target evidence). Results Functional enrichment revealed significant enrichment in cell adhesion and TGF-beta signaling. PPI network analysis identified CD44 as the top hub gene (degree=42). Drug repurposing identified 106 candidate drugs, including troglitazone and enzalutamide. Using the triple intersection strategy, five genes (ID2, NR4A1, GJA5, ID1, MYH11) were prioritized for molecular docking. GJA5 showed strong predicted binding affinity with flufenamic acid (−7.88 kcal/mol), and cytosporone B exhibited favorable druglikeness (0 Lipinski violations). Conclusion This study systematically characterizes PCOS-associated gene networks and provides a prioritized set of candidate targets and drugs through a purely computational framework. CD44 emerges as a key network node with potential relevance in PCOS pathophysiology. These findings offer testable hypotheses for future mechanistic studies and drug discovery efforts in PCOS.