OBJECTIVE To construct an intelligent pharmaceutical Q&A platform for precision medication with low “artificial intelligence （AI） hallucination”， aiming to enhance the accuracy， consistency， and traceability of medication consultations. METHODS Medication package inserts were batch-processed and converted into structured data through Python programming to build a local pharmaceutical knowledge base. The retrieval and question-answering processes were designed based on large language models， and system integration and localized deployment were completed on Dify platform. By designing typical clinical medication questions and comparing the output of the intelligent pharmaceutical Q&A platform with the online version of DeepSeek across dimensions such as peak time retrieval， half-life， and dosage adjustment reasoning for patients with renal impairment， the accuracy and reliability of its retrieval and reasoning results were evaluated. RESULTS The intelligent pharmaceutical Q&A platform， constructed based on local drug package inserts， achieved 100% accuracy in retrieval and reasoning for peak time， half-life， and dosage adjustment schemes. In comparison， the online version of DeepSeek demonstrated accuracies of 30%（6/20）， 50%（10/20）， and 38%（23/60） across these three dimensions， respectively. CONCLUSIONS The constructed intelligent pharmaceutical Q&A platform is capable of accurately retrieving and extracting information from the local knowledge base based on clinical inquiries， thereby avoiding the occurrence of AI hallucinations and providing reliable medication decision support for healthcare professionals.

Metabolic dysfunction-associated steatotic liver disease (MASLD) has become the most prevalent chronic liver disease worldwide, affecting approximately 32%-38% of the adult population and posing a growing public health burden. MASLD represents a continuous disease spectrum ranging from simple steatosis to metabolic dysfunction-associated steatohepatitis (MASH), progressive hepatic fibrosis, cirrhosis, and ultimately hepatocellular carcinoma (HCC). The pathological core of MASLD lies in disruption of hepatic lipid metabolic homeostasis, characterized by an imbalance among de novo lipogenesis, fatty acid β-oxidation, and very-low-density lipoprotein (VLDL)-mediated lipid export. This metabolic disequilibrium subsequently drives inflammatory injury and fibrotic progression. Among the multiple regulatory pathways involved, thyroid hormone (TH) signaling has emerged as a central regulator of hepatic metabolic homeostasis. The liver is a major peripheral target organ of TH action, where TH predominantly exerts its metabolic effects through thyroid hormone receptor β (TRβ). Large-scale epidemiological studies and meta-analyses have demonstrated that hypothyroidism is significantly associated with increased MASLD prevalence, more severe histological injury, and advanced hepatic fibrosis, suggesting that dysregulation of TH signaling may participate throughout the entire MASLD disease spectrum. At the molecular level, TH regulates hepatic lipid metabolism by coordinating suppression of lipogenesis, enhancement of mitochondrial fatty acid oxidation, and promotion of VLDL assembly and secretion through integrated genomic actions of the T3-TRβ axis and non-genomic signaling pathways. Across different stages of MASLD, TH signaling exerts stage-dependent protective effects. In the steatosis stage, TH improves metabolic flexibility by modulating insulin sensitivity, glucose metabolism, and lipid droplet clearance, thereby alleviating early lipotoxic stress. During progression to MASH, TH attenuates inflammatory amplification by improving mitochondrial homeostasis, suppressing activation of the NOD-like receptor family pyrin domain containing 3 (NLRP3) inflammasome, and modulating the gut-liver axis microenvironment. In advanced stages, TH signaling influences hepatic stellate cell activation and extracellular matrix deposition, partly through interaction with the transforming growth factor-β (TGF-β)/SMAD pathway, while alterations in intrahepatic TH availability, mediated by dynamic changes in iodothyronine deiodinase 1 (DIO1), contribute to fibrosis progression and hepatocellular dedifferentiation. In hepatocellular carcinoma, coordinated downregulation of TRβ and DIO1 establishes a tumor-associated hypothyroid state that promotes metabolic reprogramming and tumor progression. The clinical relevance of TH signaling in MASLD has been underscored by the recent approval of Resmetirom, a liver-targeted TRβ‑selective agonist, for the treatment of non-cirrhotic MASH with moderate-to-severe fibrosis (F2-F3). This approval represents a landmark transition from mechanistic understanding to metabolism-centered precision therapy in MASLD. Clinical trials have demonstrated that Resmetirom not only improves key histological endpoints, including MASH resolution and fibrosis regression, but also favorably modulates atherogenic lipid profiles, highlighting the therapeutic potential of selectively targeting hepatic TH pathways. This review systematically summarizes the multidimensional regulatory roles of TH across the MASLD disease spectrum and discusses emerging diagnostic and therapeutic implications of TH-based interventions, aiming to inform future mechanistic research and optimize clinical management strategies.

IgA nephropathy (IgAN) is the most common primary glomerulonephritis worldwide and a major cause of chronic kidney disease and kidney failure. IgAN exhibits marked heterogeneity in clinical presentation, histopathology, and pathogenic mechanisms, contributing to variable treatment responses and prognosisamong patients. Precise risk assessment and individualized intervention are therefore of critical importance. This review systematically traces the evolution of IgAN management from traditional risk stratification toward biomarker-driven precision medicine. We first review the clinical utility and limitations of established risk stratification tools, including the KDIGO guidelines, the Oxford MEST-C classification, and the International IgAN Prediction Tool. We then discuss emerging biomarkers closely linked to disease pathogenesis, including galactose-deficient IgA1 (Gd-IgA1), anti-Gd-IgA1 autoantibodies, B cell activating factor (BAFF), a proliferation-inducing ligand (APRIL), and complement components, as well as the targeted therapies they have informed. In addition, urinary biomarkers and multi-omics approaches show promise for dynamic disease monitoring and individualized risk stratification.

OBJECTIVE To establish a high-risk medication list and preventive and therapeutic measures for drug-induced hypofibrinogenemia， and to provide a reference for the prevention and treatment of this condition. METHODS By integrating domestic and international case reports， retrospective case-control studies， and spontaneous adverse drug reaction reporting databases， 19 domestically marketed high-risk drugs for drug-induced hypofibrinogenemia were identified. Based on the clinical characteristics and mechanisms of these drugs， relevant risk factors were systematically reviewed， and existing treatment options were summarized， leading to the preliminary development of recommended preventive and therapeutic measures. A two-round Delphi consultation was conducted to evaluate， revise， and ultimately reach consensus on the preliminary findings， using a mean importance score of ≥3.5 points for indicators and a coefficient of variation ＜0.3 as screening criteria. RESULTS The coefficient of expert authority for both rounds of expert consultation was 0.904. In the first round， the Kendall coordination coefficients （Kendall’s W ） for the high-risk medication list and the proposed preventive and therapeutic measures were 0.390 and 0.223 （ P ＜0.05）， respectively. In the second round， the Kendall’s W were 0.227 and 0.200 （ P ＜0.05）， respectively. After two rounds of expert consultation and discussion， 11 high-risk drugs for drug-induced hypofibrinogenemia， represented by hemocoagulase and certain anti-infective agents， were ultimately identified， along with 5 preventive and therapeutic measures spanning the entire process of “pre-medication assessment， intra-medication monitoring， and bleeding event management”. CONCLUSIONS This study has established a scientific and reliable high-risk medication list， and corresponding preventive and therapeutic measures for drug-induced hypofibrinogenemia， providing a theoretical basis and practical support for the early identification， stratified management， and precise intervention of this condition.

The study focuses on the concept of multifunctional traditional Chinese medicine (TCM) formulas and aims to evaluate the efficacy of the classical formula Xiaoyao San (逍遥散). Study employs the integrated evidence chain (Eff-iEC) method to organize, integrate, and evaluate its therapeutic efficacy in treating different diseases with the same therapy, and to investigate the feasibility of using Eff-iEC to evaluate the multifunctionality of TCM formulas. The evaluation covered Xiaoyao San's therapeutic effects on depression, premenstrual syndrome, chronic hepatitis, irritable bowel syndrome, dyspepsia, and menopausal syndrome. Concurrently, the Grading of Recommendations, Assessment, Development, and Evaluation (GRADE) system was used for evaluation, and authoritative medical documents were incorporated to corroborate the recognition of Xiaoyao San within the medical community. Depression and menopausal syndrome received higher ratings than other conditions in the Eff-iEC, GRADE, and Medical Community Recognition assessments. The Eff-iEC evidence grade for Xiaoyao San was rated as "High" or above for chronic hepatitis, irritable bowel syndrome, dyspepsia, and menopausal syndrome. Premenstrual syndrome received a "Moderate ⁺" rating. The GRADE evidence level was "Low-〇〇⨁⨁" for depression, premenstrual syndrome, and chronic hepatitis; "Moderate-〇⨁⨁⨁" for dyspepsia and menopausal syndrome; and "Very Low-〇〇〇⨁" for irritable bowel syndrome. Depression and menopausal syndrome had the highest inclusion frequency, appearing in all 4 categories. Premenstrual syndrome, chronic hepatitis, and dyspepsia are not recommended in Western medical guidelines, but they are included in TCM guidelines, the China National Basic Medical Insurance Drug List, and the China National Essential Drug List. Irritable bowel syndrome appears only in the China National Basic Medical Insurance Drug List and China National Essential Drug List. The evaluation results obtained using the Eff-iEC method align with Medical Community Recognition, providing an objective and comprehensive assessment of Xiaoyao San's efficacy. The findings suggest that Xiaoyao San has strong evidence for treating depression and menopausal syndrome. However, further experimental and clinical trials are needed to assess its efficacy in treating premenstrual syndrome, chronic hepatitis, irritable bowel syndrome, and dyspepsia. These results support the clinical efficacy and rational use of Xiaoyao San, expand the application scope of the Eff-iEC method, and offer valuable insights and methodological references for the comparative evaluation of multifunctional TCM formulas.

ObjectiveSpinal cord injury (SCI) directly impairs the regulatory function of the autonomic nervous system, induces intestinal dysfunction, and significantly reduces patients’ quality of life. Preclinical studies have shown that electroacupuncture (EA) therapy can regulate the brain-gut axis and is used to treat central nervous system diseases such as major depressive disorder, Alzheimer’s disease and Parkinson’s disease. Recent research has established that fecal microbiota transplantation (FMT) from EA-treated SCI rats restored intestinal motility and colonic morphology. However, it remains unclear whether the regulation of gut microbiota by EA therapy directly contributes to neural repair after SCI. This study aims to explore whether gut microbiota mediates the neuroprotective effect of EA in the treatment of SCI and its possible mechanism. MethodsThe study employed RNA transcriptome analysis of spinal cord tissue to characterize gene expression profiles and to identify key signaling pathways following EA treatment for SCI. Hematoxylin-Eosin (HE) staining and Nissl staining were used to observe the morphological changes in spinal cord tissue. Western blot (WB) and enzyme-linked immunosorbent assay (ELISA) were applied to detect the effects of EA on the expression of proteins related to nucleotide-binding domain leucine-rich repeat and pyrin domain-containing receptor 3 (NLRP3) -dependent pyroptosis. Using 16S rDNA sequencing, the study observed alterations in gut microbiota diversity and community composition in SCI rats. Prior to establishing SCI models, rats were pretreated with an antibiotic cocktail to induce gut dysbiosis, and the effects on intestinal function and spinal cord neural repair were evaluated. FMT was performed to investigate the regulatory effects of post-EA FMT on motor function, general status, liver and spleen indices, and NLRP3-mediated pyroptosis in SCI rats. ResultsEA improved motor function and reduced regulated neuronal cell death in SCI rats. Transcriptomic analysis demonstrated the activation of immune- and inflammation-related pathways post-SCI, including NOD-like receptors, nuclear factor-kappa B(NF-κB), and Toll-like receptor (TLR) pathways. EA primarily influenced intestinal inflammation and autoimmune functions. 16S rDNA sequencing illustrated that EA did not alter the diversity of gut microbiota. However, EA altered the gut microbiota composition in SCI rats, increasing Lactobacillus and Akkermansia genera while rebalancing the Firmicutes/Bacteroidetes ratio. Furthermore, depletion of gut microbiota by antibiotics disrupted the intestinal barrier, reduced the expression of intestinal barrier proteins Zonula Occludens-1 (ZO-1) and Occludin, elevated serum lipopolysaccharide-binding protein (LBP) levels, exacerbated spinal cord tissue damage, and hindered motor function recovery in SCI rats. FMT from donors treated with EA reduced LBP levels in the intestine, blood, and spinal cord of rats, inhibited the TLR4 myeloid differentiation primary response protein 88 (MyD88)-NF‑κB pathway and NLRP3-dependent pyroptosis, and improved motor function. On the other hand, FMT treatment resulted in decreased body weight and food intake, whereas FMT using EA-treated donors effectively alleviated these alterations. ConclusionEA effectively alleviated neuroinflammatory responses in rats with SCI, primarily through regulating the gut microbiota and suppressing the NLRP3-dependent pyroptosis signaling pathway.

The initiation of adaptive immune responses relies on the precise recognition and interpretation of antigenic information. In this process, the specific binding of T cell receptors (TCRs) to peptide-major histocompatibility complex (pMHC) molecules represents one of the key molecular events in the initiation of adaptive immune responses. Accordingly, the structural features of TCR-pMHC complexes provide a fundamental basis for dissecting antigen recognition mechanisms and support rational vaccine design, therapeutic target discovery in TCR-based immunotherapy, and TCR identification and optimization. However, experimental determination of TCR-pMHC structures remains costly, time-consuming, and limited in coverage, making computational approaches essential for rapidly obtaining reliable structural information. Computational methods for predicting the structures of TCR-pMHC complexes have advanced rapidly in recent years, driven by progress in deep learning-based modeling frameworks and the increasing availability of structural and sequence resources. Despite these developments, most existing tools do not adequately distinguish the key structural and biophysical differences between MHC class I (MHC-I) and MHC class II (MHC-II) complexes during model construction. As a consequence, their predictive performance differs substantially between class I and class II complexes. In general, structural predictions for class I complexes outperform those for class II complexes. This discrepancy may be related to several fundamental differences between the two systems, including the architecture of the peptide-binding groove, the distribution of peptide lengths, and the properties of peptide flanking residues (PFRs). Compared with MHC-I molecules, MHC-II molecules usually bind longer antigenic peptides, which typically range from 13 to 25 amino acids in length. PFRs at both termini of these peptides participate in regulating the overall conformation of TCR-pMHC class II complexes and exert a pronounced effect on the geometric and physicochemical characteristics of the TCR-pMHC binding interface. Furthermore, within the TCR recognition interface, the complementarity-determining regions (CDRs) consist of segments that differ markedly in conformational behavior. They commonly include regions that are relatively rigid and structurally stable, together with highly flexible segments exhibiting substantial conformational plasticity. These rigidity-flexibility features constitute an essential structural basis enabling TCRs to recognize diverse peptide-MHC ligands and to accommodate conformational heterogeneity at the interface. However, many current modeling tools, in an effort to enforce global conformational stability or reduce structural noise, tend to over-constrain intrinsically flexible regions. Such oversimplification may lead to inappropriate rigidification of flexible CDR loops, resulting in local structural distortions, compromised interface geometry, or even complete modeling failure for specific complexes. Against this background, the review approaches the field from the perspective of computational differences between MHC-I and MHC-II complexes. We first systematically organize and summarize available resources related to TCRs and pMHCs, including structural datasets, sequence databases, prediction tools, and benchmarking studies. We then focus on five representative tools capable of predicting both class I and class II complexes—AlphaFold2, AlphaFold3, TCRmodel2, tFold-TCR, and TCR-pHLA_ModellerS. After excluding structures present in the training sets of these tools, we constructed a benchmark dataset comprising 25 class I and 10 class II TCR-pMHC complexes in the bound state and conducted a systematic evaluation using this dataset. We first employ widely used general evaluation metrics, including All-Atom Root Mean Square Deviation (All-Atom RMSD), Backbone RMSD, Template Modeling score (TM-score), and DockQ, to assess the global conformational accuracy and interface modeling quality of class I and class II complexes. For class II complexes, we propose for the first time a peptide flanking residue deviation index, including the PFRs-Deviation Index (PFRs-DI), N-PFR-Deviation Index (N-PFR-DI), and C-PFR-Deviation Index (C-PFR-DI), to quantitatively characterize conformational deviations in PFRs. In addition, we propose the CDR conformational consistency index (CCC) designed to qualitatively evaluate the ability of prediction tools to capture TCR CDR conformational flexibility. These metrics collectively assess a tool’s ability to model both overall conformation and critical functional regions, thereby addressing the limitations of existing evaluation criteria that overemphasize global structure while inadequately capturing modeling quality in key functional areas. This establishes a unified analytical framework for MHC-I and MHC-II complexes to guide data resource selection, modeling strategy formulation, and evaluation system development. The framework further advances computational modeling and provides crucial support for multi-scale analysis of TCR-pMHC recognition mechanisms and their biological functions.

ObjectiveSpinal cord injury (SCI) directly impairs the regulatory function of the autonomic nervous system, induces intestinal dysfunction, and significantly reduces patients’ quality of life. Preclinical studies have shown that electroacupuncture (EA) therapy can regulate the brain-gut axis and is used to treat central nervous system diseases such as major depressive disorder, Alzheimer’s disease and Parkinson’s disease. Recent research has established that fecal microbiota transplantation (FMT) from EA-treated SCI rats restored intestinal motility and colonic morphology. However, it remains unclear whether the regulation of gut microbiota by EA therapy directly contributes to neural repair after SCI. This study aims to explore whether gut microbiota mediates the neuroprotective effect of EA in the treatment of SCI and its possible mechanism. MethodsThe study employed RNA transcriptome analysis of spinal cord tissue to characterize gene expression profiles and to identify key signaling pathways following EA treatment for SCI. Hematoxylin-Eosin (HE) staining and Nissl staining were used to observe the morphological changes in spinal cord tissue. Western blot (WB) and enzyme-linked immunosorbent assay (ELISA) were applied to detect the effects of EA on the expression of proteins related to nucleotide-binding domain leucine-rich repeat and pyrin domain-containing receptor 3 (NLRP3) -dependent pyroptosis. Using 16S rDNA sequencing, the study observed alterations in gut microbiota diversity and community composition in SCI rats. Prior to establishing SCI models, rats were pretreated with an antibiotic cocktail to induce gut dysbiosis, and the effects on intestinal function and spinal cord neural repair were evaluated. FMT was performed to investigate the regulatory effects of post-EA FMT on motor function, general status, liver and spleen indices, and NLRP3-mediated pyroptosis in SCI rats. ResultsEA improved motor function and reduced regulated neuronal cell death in SCI rats. Transcriptomic analysis demonstrated the activation of immune- and inflammation-related pathways post-SCI, including NOD-like receptors, nuclear factor-kappa B(NF-κB), and Toll-like receptor (TLR) pathways. EA primarily influenced intestinal inflammation and autoimmune functions. 16S rDNA sequencing illustrated that EA did not alter the diversity of gut microbiota. However, EA altered the gut microbiota composition in SCI rats, increasing Lactobacillus and Akkermansia genera while rebalancing the Firmicutes/Bacteroidetes ratio. Furthermore, depletion of gut microbiota by antibiotics disrupted the intestinal barrier, reduced the expression of intestinal barrier proteins Zonula Occludens-1 (ZO-1) and Occludin, elevated serum lipopolysaccharide-binding protein (LBP) levels, exacerbated spinal cord tissue damage, and hindered motor function recovery in SCI rats. FMT from donors treated with EA reduced LBP levels in the intestine, blood, and spinal cord of rats, inhibited the TLR4 myeloid differentiation primary response protein 88 (MyD88)-NF‑κB pathway and NLRP3-dependent pyroptosis, and improved motor function. On the other hand, FMT treatment resulted in decreased body weight and food intake, whereas FMT using EA-treated donors effectively alleviated these alterations. ConclusionEA effectively alleviated neuroinflammatory responses in rats with SCI, primarily through regulating the gut microbiota and suppressing the NLRP3-dependent pyroptosis signaling pathway.

The initiation of adaptive immune responses relies on the precise recognition and interpretation of antigenic information. In this process, the specific binding of T cell receptors (TCRs) to peptide-major histocompatibility complex (pMHC) molecules represents one of the key molecular events in the initiation of adaptive immune responses. Accordingly, the structural features of TCR-pMHC complexes provide a fundamental basis for dissecting antigen recognition mechanisms and support rational vaccine design, therapeutic target discovery in TCR-based immunotherapy, and TCR identification and optimization. However, experimental determination of TCR-pMHC structures remains costly, time-consuming, and limited in coverage, making computational approaches essential for rapidly obtaining reliable structural information. Computational methods for predicting the structures of TCR-pMHC complexes have advanced rapidly in recent years, driven by progress in deep learning-based modeling frameworks and the increasing availability of structural and sequence resources. Despite these developments, most existing tools do not adequately distinguish the key structural and biophysical differences between MHC class I (MHC-I) and MHC class II (MHC-II) complexes during model construction. As a consequence, their predictive performance differs substantially between class I and class II complexes. In general, structural predictions for class I complexes outperform those for class II complexes. This discrepancy may be related to several fundamental differences between the two systems, including the architecture of the peptide-binding groove, the distribution of peptide lengths, and the properties of peptide flanking residues (PFRs). Compared with MHC-I molecules, MHC-II molecules usually bind longer antigenic peptides, which typically range from 13 to 25 amino acids in length. PFRs at both termini of these peptides participate in regulating the overall conformation of TCR-pMHC class II complexes and exert a pronounced effect on the geometric and physicochemical characteristics of the TCR-pMHC binding interface. Furthermore, within the TCR recognition interface, the complementarity-determining regions (CDRs) consist of segments that differ markedly in conformational behavior. They commonly include regions that are relatively rigid and structurally stable, together with highly flexible segments exhibiting substantial conformational plasticity. These rigidity-flexibility features constitute an essential structural basis enabling TCRs to recognize diverse peptide-MHC ligands and to accommodate conformational heterogeneity at the interface. However, many current modeling tools, in an effort to enforce global conformational stability or reduce structural noise, tend to over-constrain intrinsically flexible regions. Such oversimplification may lead to inappropriate rigidification of flexible CDR loops, resulting in local structural distortions, compromised interface geometry, or even complete modeling failure for specific complexes. Against this background, the review approaches the field from the perspective of computational differences between MHC-I and MHC-II complexes. We first systematically organize and summarize available resources related to TCRs and pMHCs, including structural datasets, sequence databases, prediction tools, and benchmarking studies. We then focus on five representative tools capable of predicting both class I and class II complexes—AlphaFold2, AlphaFold3, TCRmodel2, tFold-TCR, and TCR-pHLA_ModellerS. After excluding structures present in the training sets of these tools, we constructed a benchmark dataset comprising 25 class I and 10 class II TCR-pMHC complexes in the bound state and conducted a systematic evaluation using this dataset. We first employ widely used general evaluation metrics, including All-Atom Root Mean Square Deviation (All-Atom RMSD), Backbone RMSD, Template Modeling score (TM-score), and DockQ, to assess the global conformational accuracy and interface modeling quality of class I and class II complexes. For class II complexes, we propose for the first time a peptide flanking residue deviation index, including the PFRs-Deviation Index (PFRs-DI), N-PFR-Deviation Index (N-PFR-DI), and C-PFR-Deviation Index (C-PFR-DI), to quantitatively characterize conformational deviations in PFRs. In addition, we propose the CDR conformational consistency index (CCC) designed to qualitatively evaluate the ability of prediction tools to capture TCR CDR conformational flexibility. These metrics collectively assess a tool’s ability to model both overall conformation and critical functional regions, thereby addressing the limitations of existing evaluation criteria that overemphasize global structure while inadequately capturing modeling quality in key functional areas. This establishes a unified analytical framework for MHC-I and MHC-II complexes to guide data resource selection, modeling strategy formulation, and evaluation system development. The framework further advances computational modeling and provides crucial support for multi-scale analysis of TCR-pMHC recognition mechanisms and their biological functions.

ObjectiveTraditional Chinese medicine (TCM) constitutes a valuable cultural heritage and an important source of antitumor compounds. Poria (Poria cocos (Schw.) Wolf), the dried sclerotium of a polyporaceae fungus, was first documented in Shennong’s Classic of Materia Medica and has been used therapeutically and dietarily in China for millennia. Traditionally recognized for its diuretic, spleen-tonifying, and sedative properties, modern pharmacological studies confirm that Poria exhibits antioxidant, anti-inflammatory, antibacterial, and antitumor activities. Pachymic acid (PA; a triterpenoid with the chemical structure 3β-acetyloxy-16α-hydroxy-lanosta-8,24(31)-dien-21-oic acid), isolated from Poria, is a principal bioactive constituent. Emerging evidence indicates PA exerts antitumor effects through multiple mechanisms, though these remain incompletely characterized. Neuroblastoma (NB), a highly malignant pediatric extracranial solid tumor accounting for 15% of childhood cancer deaths, urgently requires safer therapeutics due to the limitations of current treatments. Although PA shows multi-mechanistic antitumor potential, its efficacy against NB remains uncharacterized. This study systematically investigated the potential molecular targets and mechanisms underlying the anti-NB effects of PA by integrating network pharmacology-based target prediction with experimental validation of multi-target interactions through molecular docking, dynamic simulations, and in vitro assays, aimed to establish a novel perspective on PA’s antitumor activity and explore its potential clinical implications for NB treatment by integrating computational predictions with biological assays. MethodsThis study employed network pharmacology to identify potential targets of PA in NB, followed by validation using molecular docking, molecular dynamics (MD) simulations, MM/PBSA free energy analysis, RT-qPCR and Western blot experiments. Network pharmacology analysis included target screening via TCMSP, GeneCards, DisGeNET, SwissTargetPrediction, SuperPred, and PharmMapper. Subsequently, potential targets were predicted by intersecting the results from these databases via Venn analysis. Following target prediction, topological analysis was performed to identify key targets using Cytoscape software. Molecular docking was conducted using AutoDock Vina, with the binding pocket defined based on crystal structures. MD simulations were performed for 100 ns using GROMACS, and RMSD, RMSF, SASA, and hydrogen bonding dynamics were analyzed. MM/PBSA calculations were carried out to estimate the binding free energy of each protein-ligand complex. In vitro validation included RT-qPCR and Western blot, with GAPDH used as an internal control. ResultsThe CCK-8 assay demonstrated a concentration-dependent inhibitory effect of PA on NB cell viability. GO analysis suggested that the anti-NB activity of PA might involve cellular response to chemical stress, vesicle lumen, and protein tyrosine kinase activity. KEGG pathway enrichment analysis suggested that the anti-NB activity of PA might involve the PI3K/AKT, MAPK, and Ras signaling pathways. Molecular docking and MD simulations revealed stable binding interactions between PA and the core target proteins AKT1, EGFR, SRC, and HSP90AA1. RT-qPCR and Western blot analyses further confirmed that PA treatment significantly decreased the mRNA and protein expression of AKT1, EGFR, and SRC while increasing the HSP90AA1 mRNA and protein levels. ConclusionIt was suggested that PA may exert its anti-NB effects by inhibiting AKT1, EGFR, and SRC expression, potentially modulating the PI3K/AKT signaling pathway. These findings provide crucial evidence supporting PA’s development as a therapeutic candidate for NB.