Identification of inhibitors from a functional food-based plant Perillae Folium against hyperuricemia via metabolomics profiling, network pharmacology and all-atom molecular dynamics simulations

IntroductionHyperuricemia (HUA) is a metabolic disorder caused by purine metabolism dysfunction in which the increasing purine levels can be partially attributed to seafood consumption. Perillae Folium (PF), a widely used plant in functional food, has been historically used to mitigate seafood-induced diseases. However, its efficacy against HUA and the underlying mechanism remain unclear. MethodsA network pharmacology analysis was performed to identify candidate targets and potential mechanisms involved in PF treating HUA. The candidate targets were determined based on TCMSP, SwissTargetPrediction, Open Targets Platform, GeneCards, Comparative Toxicogenomics Database, and DrugBank. The potential mechanisms were predicted via Gene Ontology (GO) and Kyoto Gene and Genome Encyclopedia (KEGG) analyses. Molecular docking in AutoDock Vina and PyRx were performed to predict the binding affinity and pose between herbal compounds and HUA-related targets. A chemical structure analysis of PF compounds was performed using OSIRIS DataWarrior and ClassyFire. We then conducted virtual pharmacokinetic and toxicity screening to filter potential inhibitors. We further performed verifications of these inhibitors’ roles in HUA through molecular dynamics (MD) simulations, text-mining, and untargeted metabolomics analysis. ResultsWe obtained 8200 predicted binding results between 328 herbal compounds and 25 potential targets, and xanthine dehydrogenase (XDH) exhibited the highest average binding affinity. We screened out five promising ligands (scutellarein, benzyl alpha-D-mannopyranoside, elemol, diisobutyl phthalate, and (3R)-hydroxy-beta-ionone) and performed MD simulations up to 50 ns for XDH complexed to them. The scutellarein-XDH complex exhibited the most satisfactory stability. Furthermore, the text-mining study provided laboratory evidence of scutellarein’s function. The metabolomics approach identified 543 compounds and confirmed the presence of scutellarein. Extending MD simulations to 200 ns further indicated the sustained impact of scutellarein on XDH structure. ConclusionOur study provides a computational and biomedical basis for PF treating HUA and fully elucidates scutellarein's great potential as an XDH inhibitor at the molecular level, holding promise for future drug design and development.