Elucidating the Therapeutic Mechanism of Puhuaiyin in Rosacea: An Integrated Study Based on Network Pharmacology, Bioinformatics, Molecular Docking, and Molecular Dynamics Simulations.
BACKGROUND Rosacea is a chronic inflammatory skin disorder with limited therapeutic options. Puhuaiyin (PHY), a traditional Chinese medicinal formula, shows clinical efficacy, but its multi-component mechanisms remain unclear. METHODS Chemical constituents of PHY were identified by UPLC-Q-TOF-MS. Network pharmacology was used to predict potential targets, which were intersected with rosacea-associated genes. Bioinformatics analyses (differential expression, WGCNA, and machine learning) were applied to the GEO dataset GSE65914 to refine core targets. Molecular docking and molecular dynamics simulations were conducted to validate the binding modes and stability between key active constituents and the core targets. RESULTS A total of 59 chemical constituents were identified in PHY, with five key active components subsequently screened: quercetin, emodin, kushenol N, physcion, and palmitic acid. Network pharmacology analysis revealed 44 intersecting targets, which were significantly enriched in inflammation-related pathways, such as MAPK, NF-κB, and JAK-STAT signaling. Integrated bioinformatics and machine learning approaches identified MMP9 and IL1B as core targets, both of which were markedly upregulated in rosacea lesions and demonstrated prominent diagnostic value (AUC = 0.999 for MMP9, 0.964 for IL1B). Molecular docking indicated strong binding affinity between the core components and MMP9/IL1B. Molecular dynamics simulations confirmed stable complex conformations over 200 ns, with MM/PBSA binding free energies of -15.54 Kcal/mol (quercetin-MMP9) and -15.57 Kcal/mol (quercetin- IL1B). DISCUSSION This study, through a multidisciplinary approach, systematically elucidates the "multi-component, multi-target, and multi-pathway" mode of action of PHY in the treatment of rosacea. However, the computational predictions remain to be further validated by in vivo and in vitro experiments. Future research should focus on verifying its therapeutic efficacy in animal or cellular models, as well as elucidating the regulatory effects of key active components on the MMP9 and IL1B targets. CONCLUSION These computational predictions suggest that PHY may exert therapeutic effects against rosacea via quercetin and other components targeting MMP9 and IL1B, thereby modulating MAPK, NF-κB, and JAK-STAT pathways. The proposed mechanisms include inhibition of inflammation, regulation of the immune microenvironment, attenuation of vascular dilation, and promotion of skin barrier recovery. These findings provide a theoretical basis for future experimental validation.