收藏切换
AI-Powered Discovery of Rosmarinic Acid as a Novel Ferroptosis Inhibitor for Ulcerative Colitis via Targeting the ALOX15–VDAC1 Axis
收藏切换
PDF
Guixuan Fang1, Jingli Ao1, Xiao Li1, Weiwei Wang1, Yamina Alioui1, Yang Huang1, Jiali Peng1, Bufu Tang2, *, Lianxiang Luo1, *
Research. Vol 9 Article ID 1352
Less
收藏切换
Research. Vol 9 Article ID 1352
• RESEARCH ARTICLE •
AI-Powered Discovery of Rosmarinic Acid as a Novel Ferroptosis Inhibitor for Ulcerative Colitis via Targeting the ALOX15–VDAC1 Axis
Full
Guixuan Fang1, Jingli Ao1, Xiao Li1, Weiwei Wang1, Yamina Alioui1, Yang Huang1, Jiali Peng1, Bufu Tang2, *, Lianxiang Luo1, *
Affiliations
  • 1School of Ocean and Tropical Medicine, Guangdong Medical University, Zhanjiang, Guangdong 524023, China.
  • 2Department of Interventional Radiology, Zhongshan Hospital, Fudan University, Shanghai 200032, China.
About Author:

† These authors contributed equally to this work.

Published: 2026-07-20 doi: 10.34133/research.1352
Outline
收藏切换

Ulcerative colitis (UC) is a chronic inflammatory bowel disease with limited mechanism-based therapies, and ferroptosis—an iron-dependent lipid peroxidation-driven cell death—has emerged as a critical pathogenic event. Using a machine learning model that combines molecular fingerprints and transfer learning representations to screen for ALOX15 inhibitors, we identified rosmarinic acid (RA) from natural product libraries. RA exhibited high-affinity binding to ALOX15 and maintained stable complex conformation in molecular dynamics simulations. In dextran sulfate sodium-induced colitis mice and lipopolysaccharide-treated intestinal epithelial cells (IEC-6 and mouse colonic epithelial cell), RA alleviated disease activity index, colon shortening, histological damage, and barrier protein loss (ZO-1, Occludin, and Claudin-1). Mechanistically, RA suppressed ferroptosis by restoring GPX4, SLC7A11, and FTH1; reducing ACSL4, COX-2, and lipid reactive oxygen species (ROS); and lowering intracellular Fe2+ and malondialdehyde. Crucially, RA disrupted the ALOX15–VDAC1 protein complex, preserved mitochondrial membrane potential, and reduced mitochondrial ROS, effects that were abrogated by ALOX15 knockdown or VDAC1 knockdown. Serum metabolomics further revealed that RA remodels lipid metabolism pathways (arachidonic acid and linoleic acid metabolism) linked to ferroptosis. Collectively, this artificial intelligence-driven discovery positions RA as a novel ferroptosis inhibitor that targets the ALOX15–VDAC1 axis, offering a promising therapeutic candidate for UC.

Guixuan Fang, Jingli Ao, Xiao Li, Weiwei Wang, Yamina Alioui, Yang Huang, Jiali Peng, Bufu Tang, Lianxiang Luo. AI-Powered Discovery of Rosmarinic Acid as a Novel Ferroptosis Inhibitor for Ulcerative Colitis via Targeting the ALOX15–VDAC1 Axis[J]. Research, 2026 , 9 (7) : 1352 - . DOI: 10.34133/research.1352
Ulcerative colitis (UC) is a chronic relapsing inflammatory bowel disease characterized by mucosal inflammation of the colorectal epithelium. Its global burden continues to rise, particularly in newly industrialized countries, where incidence rates have increased markedly in recent decades [1]. Clinically, UC presents with heterogeneous manifestations; rectal bleeding or hematochezia is a hallmark feature, reported in nearly 90% of cases, and is frequently accompanied by abdominal pain, fatigue, and systemic symptoms that markedly impair patients’ health-related quality of life. At the mechanistic level, dysfunction of the intestinal epithelial barrier is widely recognized as an early and pivotal event in UC pathogenesis. Despite considerable advances, the molecular underpinnings of UC disease progression remain incompletely understood, thereby limiting the rational development of mechanism-based therapeutics [2,3].
In parallel, machine learning-assisted virtual screening has emerged as a useful computational strategy for prioritizing bioactive compounds from large chemical libraries. Such approaches can improve the efficiency of candidate prioritization and support early-stage lead identification. Notably, machine learning-based screening approaches have shown utility in exploring the structural complexity and pharmacological diversity of natural products [4,5]. These advances collectively support the application of machine learning-assisted virtual screening as a candidate-prioritization strategy for identifying bioactive natural products of potential relevance to UC.
Ferroptosis is a regulated, iron-dependent form of cell death characterized by the accumulation of lipid peroxides and disrupted metabolism of polyunsaturated fatty acids [6–8]. Accumulating evidence indicates that ferroptosis is implicated in the pathogenesis of numerous diseases [9], including malignant tumors, metabolic disorders such as diabetes, and autoimmune conditions such as psoriasis [10–17]. Importantly, experimental inhibition of ferroptosis has been shown to attenuate disease progression in animal models. Of particular relevance to UC, emerging studies indicated that ferroptotic death of intestinal epithelial cells contributes critically to barrier disruption and mucosal inflammation, thereby accelerating disease onset and exacerbation.
Rosmarinic acid (RA) is a naturally occurring polyphenol [18] with a favorable safety profile, which has demonstrated therapeutic potential in a variety of cancers [19–27] and exhibits gastroprotective effects in animal models [20–28]. Recent studies also showed that RA can reverse chemotherapy resistance and mitigate treatment-induced toxicity [29–31]. Despite these promising findings, its therapeutic potential and underlying mechanisms in UC remain to be established.
In this study, we employed a comprehensive multidisciplinary strategy. First, we applied a machine learning model to perform virtual screening of a natural product library, which rapidly identified RA as a promising candidate with high-affinity binding potential to ALOX15. Thereafter, through molecular docking, molecular dynamics (MD) simulations, and a series of in vitro and in vivo experiments, we found that RA alleviates UC through inhibition of ALOX15 and suppression of ferroptosis in intestinal epithelial cells. This work not only identifies a novel drug candidate for UC therapy but also underscores the practical value of an artificial intelligence-driven, experimentally validated research paradigm in translational medicine.
Analysis of differentially expressed genes (DEGs) between UC and control samples in the GSE117993 dataset identified 1,262 DEGs meeting the criteria of |log2FC| > 1 and P < 0.05 (Fig. 1A and B) including the core gene ALOX15. Furthermore, relative to controls, ALOX15 expression was markedly elevated in the UC group (Fig. 1C). Receiver operating characteristic (ROC) curve analysis demonstrated that ALOX15 expression effectively discriminates UC samples from healthy controls, with an area under the curve (AUC) of 0.771 (Fig. 1D), indicating that it has good diagnostic potential. Gene set enrichment analysis (GSEA) revealed that high ALOX15 expression is associated with substantial up-regulation of inflammation-related pathways (especially cytokine signaling), while basal metabolism and protein synthesis-related pathways are markedly down-regulated (Fig. 1E). In the ALOX15 high-expression group, the pathway of positive enrichment includes cytokine–cytokine receptor interaction, which exhibited the highest enrichment, indicating activation of a cytokine- and receptor-mediated inflammatory signaling network. The interleukin-17 (IL-17) signaling pathway is also substantial ly enriched, suggesting that Th17 cell activity and IL-17 secretion in UC patients are associated with ALOX15 expression levels. The crosstalk between viral proteins and cytokine receptor pathways also implies that elevated ALOX15 expression may be linked to altered viral–cytokine interactions and disrupted intestinal barrier and mucosal immune function in UC.
Gene set variation analysis (GSVA) further revealed that high ALOX15 expression is associated with up-regulation of lipid biosynthesis and coenzyme A, while the down-regulated metabolic pathway involves bile acid, carbohydrate metabolism, and mitochondrial energy metabolism (Fig. 1F). Up-regulated pathways include sphingolipid biosynthesis, fatty acid biosynthesis, and unsaturated fatty acid biosynthesis. Notably, in both GSEA and GSVA, the oxidative phosphorylation pathway was negatively correlated with high ALOX15 expression, suggesting an association with mitochondrial dysfunction in the inflammatory microenvironment. This metabolic reprogramming pattern may shed new light on the mechanistic involvement of ALOX15 in the pathogenesis and progression of UC.
Prior to modeling, standardized cleaning of the ALOX15 inhibitor data sourced from ChEMBL was performed. Simplified Molecular-Input Line-Entry System (SMILES) validity was first checked using RDKit, and invalid structures were removed. Subsequently, duplicate molecular entries were aggregated by mean activity values, reducing the dataset from 904 original records to 698 molecules suitable for modeling. Half maximal inhibitory concentration (IC50) values were converted to pIC50 [−log10(IC50)], and active compounds were defined as those with pIC50 ≥ 5.2 (corresponding to IC50 ≤ 6.3 μM) to construct a binary classification task. The resulting dataset showed a near-balanced class distribution in the scaffold-based train-development subset (positive prevalence ≈ 50.7%), whereas the independent scaffold-disjoint test set contained 63.3% active compounds. Detailed statistics of scaffold-level dataset partitioning are provided in Table S3. The activity distribution of the 698 ChEMBL compounds is shown in Fig. 2A.
The molecular features comprised a combination of molecular fingerprints and transfer learning representations (Morgan + Molecular ACCess System [MACCS] + pretrained representation). Specifically, approximately 3 million drug-like small molecules were randomly collected from the ZINC database for unsupervised pretraining without using any ALOX15 activity annotations. Because the pretraining stage was conducted exclusively on unlabeled external molecules, no activity information from the downstream ALOX15 dataset was involved during representation learning.
For all molecules, Morgan fingerprints (radius 2, 2,048 bits) and MACCS fingerprints (167 bits) were computed, followed by median imputation, low-variance filtering (threshold 0.16), and standardization. Finally, incremental principal component analysis (IPCA) was applied to compress the features into a 128-dimensional latent space. For the downstream ALOX15 classification task, each molecule was processed with the same fingerprint extraction and preprocessing pipeline, and the resulting features were concatenated with the 128-dimensional pretrained representation to form the final model input. This hybrid feature strategy enabled the model to capture both target-specific local structural signals and global chemical space priors derived from large-scale pretraining. This strategy substantially improves generalization robustness when using limited labeled samples.
Among the benchmarked models, ExtraTrees demonstrated superior performance in terms of generalization capability and stability (Fig. 2B and C). After single-model nested optimization, the final ExtraTrees model achieved ROC-AUC = 0.8584 (95% confidence interval [CI]: 0.829 to 0.889) (Fig. 2D) and precision-recall (PR)-AUC = 0.861 (95% CI: 0.826 to 0.901) (Fig. 2E) on scaffold 5-fold out-of-fold (OOF) validation. Additional OOF metrics included F1 = 0.776, Matthews correlation coefficient (MCC) = 0.574, balanced accuracy = 0.786, precision = 0.822, and recall = 0.736 (threshold = 0.48, selected exclusively from OOF predictions). The 5-fold cross-validation ROC curves (Fig. 2F) and PR curves (Fig. 2G) showed tightly clustered profiles with minimal AUC variation across folds, indicating good predictive consistency and stability across different structural subsets. Class distributions across the outer validation folds remained relatively balanced (Table S5). On the independent scaffold test set, the model reached ROC-AUC = 0.9510 (95% CI: 0.898 to 0.994) (Fig. 2H) and PR-AUC = 0.9783 (95% CI: 0.952 to 0.995) (Fig. 2I). Corresponding test-set metrics included F1 = 0.938, MCC = 0.8459, balanced accuracy = 0.9328, precision = 0.9783, and recall = 0.9000, using the unchanged OOF-derived threshold (0.48). These results support the model’s discriminative capability under scaffold-disjoint evaluation settings. However, because the independent holdout subset contained a relatively limited number of scaffold series, the broader assessment of scaffold-level generalization primarily relied on scaffold-grouped OOF validation across the train-development subset. Nevertheless, because the independent scaffold test set was relatively small (n = 79; 3 Murcko scaffolds), the observed test performance should be interpreted cautiously. In addition, Y-randomization analysis (499 label permutations) yielded a null ROC-AUC distribution centered near random expectation (mean ROC-AUC = 0.494) (Table S8), substantially lower than the observed model performance, supporting that the predictive signal was unlikely to arise from chance correlations.
Following priority screening of the candidate antioxidant small molecule library using the final model, 41 small molecules predicted to be active were obtained (Table S1). Molecular docking was performed on these small molecules, and the docking results (Table S2) showed the docking scores of the candidate molecules, among which RA achieved the highest docking score of −7.22 kcal/mol.
To gain deeper insight into the metabolic alterations induced by RA treatment in C57 mice, differential metabolite analysis was performed.
The volcano plot revealed 29 markedly up-regulated metabolites and 86 markedly down-regulated metabolites (Fig. 3A). Between the RA and dextran sulfate sodium (DSS) groups, 48 major metabolites exhibiting significant differences were screened out (variable importance in the projection [VIP] > 1, P < 0.05) (Fig. 3B). Classification of the markedly altered metabolites showed that within the major categories of glycerophospholipids and fatty acids, the number of down-regulated metabolites exceeded that of up-regulated ones (Fig. 3C). Small Molecule Pathway Database (SMPDB) enrichment analysis indicated that RA influenced lipid metabolism, including alpha linolenic acid and linoleic acid metabolism, mitochondrial beta-oxidation of long-chain saturated fatty acids, and arachidonic acid metabolism (Fig. 3D). Additionally, RA was associated with caffeine metabolism, suggesting a potential link to drug metabolism. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis corroborated these findings and further revealed associations between RA and pathways such as Necroptosis and Autophagy (Fig. 3E). Furthermore, substantial changes in several metabolites associated with reactive oxygen species (ROS) suppression were also observed following RA administration (Fig. 3F). Therefore, we propose that RA is closely associated with lipid peroxidation and ferroptosis in UC.
To elucidate the molecular basis of ligand recognition, molecular docking simulations were conducted to characterize the binding interface between ligands and target proteins. This method provides insight for drug design and is a key tool for early screening of candidate drugs. Molecular docking predicted the high-affinity binding mode of RA in the active site of ALOX15, with a GlideScore of −7.219 kcal/mol. As shown in Fig. 4A, RA forms hydrogen bonds with TRP594 and PRO647, and additional interactions including hydrogen bonding and a critical salt bridge with ARG599 were identified, collectively indicating a stable binding interaction between RA and ALOX15. As shown in Fig. 4B, the root-mean-square deviation (RMSD) analysis of the 100-ns MD trajectory demonstrates that the average RMSD of the ALOX15–RA complex is 0.1846 nm, indicating that the overall conformational fluctuation of the protein after ligand binding is relatively small, which shows that the ALOX15–RA complex maintains high structural integrity throughout the simulation process. The root-mean-square fluctuation (RMSF) diagram of protein residues (Fig. 4C) further reveals the substantial limitation of atomic movement in the ALOX15 binding pocket, with an average RMSF of 0.1240 nm. Notably, the RMSF value at the key active-site residue PRO647 was markedly below the mean, indicating that the binding site has a high degree of structural rigidity. The analysis of the radius of gyration (Rg) (Fig. 4D) shows that the average Rg value of the complex is 2.4019 nm, and the fluctuation is less than 0.08 nm throughout the simulation, confirming that the protein maintains a tight folded configuration without obvious structural loosening. Free energy landscape (FEL) analysis revealed that through the 3-dimensional drawing of RMSD–Rg free energy value (Fig. 4E and F), the conformational space of the ALOX15–RA complex presents a unique deep and wide energy well distribution, which may be an important dynamic basis for supporting a favorable free energy of binding. To assess the interaction and stability between RA and ALOX15, we employed dimethyl sulfoxide (DMSO) or protein extracts from RA-treated mouse colonic epithelial cells (MCECs) for cellular thermal shift assay (CETSA)–Western blot analysis. The results showed that ALOX15 exhibits high thermal stability in the RA treatment group. Drug affinity responsive target stability (DARTS) experiments showed that RA inhibits the protein degradation of ALOX15 through affinity protease. In addition, a competitive pull-down experiment confirmed the interaction between RA and ALOX15 protein (Fig. 4G to I). Collectively, these computational and in vitro binding data support ALOX15 as a key functional protein regulated by RA.
To investigate the therapeutic potential of RA in UC, a mouse model was established using 3% DSS. Compared with the DSS group, both pretreatment and combined treatment with RA (100 mg/kg) substantially improved disease severity by attenuating weight loss, reducing colon shortening, and lowering disease activity index (DAI) score (Fig. 5A to C). Histological examination (hematoxylin and eosin [H&E] staining) showed that RA substantially improved colonic pathology and reduced inflammatory cell infiltration, mucosal damage, and goblet cell loss (Fig. 5D). Immunohistochemistry further demonstrated that RA treatment effectively down-regulated the expression of ALOX15 and up-regulated key intestinal barrier proteins, including ZO-1, Occludin, Claudin-1, and MUC-2 (Fig. 5D and E).
To further characterize the underlying mechanism, protein expression related to ferroptosis and epithelial barrier function was analyzed. Western blot results showed that DSS challenge promoted ferroptotic cell death, as indicated by increased ACSL4, ALOX15, and COX-2 alongside decreased GPX4, FTH1, and SLC7A11. RA administration effectively counteracted these protein expression changes (Fig. 5F). Additionally, tight junction proteins including ZO-1, Occludin, and Claudin-1, which were markedly reduced by DSS, were restored upon RA treatment, reflecting restoration of epithelial barrier function (Fig. 5G).
The capacity of RA to counteract ferroptosis was assessed through a battery of cell-based experiments. A Cell Counting Kit-8 (CCK-8) viability assay first established that RA concentrations up to 30 μM had no measurable effect on MCEC survival over 24 h (Fig. 6A), ruling out confounding cytotoxic effects. Protein expression profiling in both MCEC and IEC-6 cells revealed that RA dose-dependently reversed the lipopolysaccharide (LPS)-driven suppression of FTH1, SLC7A11, and GPX4, while simultaneously attenuating the LPS-induced rise in ACSL4, ALOX15, and COX-2 (Fig. 6D). On a biochemical level, LPS exposure substantially increased intracellular Fe2+, lipid ROS, and malondialdehyde (MDA), all of which are hallmark indicators of ferroptotic stress, while glutathione (GSH) levels dropped accordingly. Coadministration of RA brought all these parameters back toward control values (Fig. 6B to E). Taken together, the data point to RA as an effective suppressor of LPS-triggered ferroptosis across both cell lines tested.
To confirm that the protective effect of RA is mediated through ALOX15, experiments were conducted in ALOX15 knockdown (KD) cells. ALOX15 depletion abolished the ability of RA to rescue LPS-induced GPX4 down-regulation (Fig. 7A to E). Furthermore, in ALOX15-KD cells, the drug’s ability to inhibit intracellular Fe2+ accumulation was substantially diminished (Fig. 7B). Taken together, these data confirm that RA’s anti-ferroptotic protection is largely abrogated after ALOX15 depletion, establishing ALOX15 as a major effector of RA-mediated cytoprotection
To determine whether ferroptosis occurs in LPS-induced MCECs, we examined 3 key hallmarks: lipid peroxidation, Fe2+ accumulation, and antioxidant system dysfunction. Flow cytometry analysis revealed that lipid peroxidation levels were markedly elevated in cells treated with the ALOX15 inhibitor combined with LPS (Fig. 8A). Consistently, intracellular Fe2+ levels were also increased under the same conditions (Fig. 8B). Western blot analysis further showed that combined ALOX15 inhibitor and LPS treatment elevated ALOX15 protein levels while reducing GPX4 expression (Fig. 8C), collectively confirming the occurrence of ferroptosis in these cells.
Notably, ALOX15 protein levels were elevated in LPS-treated MCECs even in the presence of the ALOX15 inhibitor, suggesting compensatory up-regulation under inflammatory conditions. Importantly, cotreatment with RA markedly reversed lipid peroxidation, Fe2+ accumulation, and ferroptosis-related protein expression changes induced by PD146176 and LPS (Fig. 8A to C). Immunofluorescence staining and mitochondrial membrane potential (MMP) assays showed consistent trends (Fig. 8D and E), further supporting these findings. Taken together, these results identify ALOX15 as a critical mediator of RA’s protective action against LPS-induced ferroptosis in MCECs, with RA exerting its effect through inhibitory modulation of ALOX15 activity.
To characterize overall protein expression differences between sample groups, principal component analysis (PCA) revealed clear separation between the RA treatment group and the LPS group (Fig. 9A). To further investigate the mechanistic basis of these proteomic changes, we compared protein expression profiles in RA-treated and untreated MCECs. Using thresholds of P < 0.05 and |log₂(fold change)| > 2, a total of 65 differentially expressed proteins (DEPs) were identified (Fig. 9B). Hierarchical clustering analysis displayed in a heatmap revealed substantial differences in the top 30 DEPs between the LPS and RA treatment groups (Fig. 9C).
Following RA treatment, Gene Ontology (GO) analysis showed that the most notably enriched biological processes included cytoplasmic translation, cellular response to ROS, and positive regulation of mitophagy (Fig. 9D), suggesting a functional link between RA treatment and modulation of oxidative stress and protein homeostasis. Notably, KEGG pathway analysis showed that DEPs were markedly enriched in neutrophil extracellular trap formation, necroptosis, and oxidative phosphorylation pathways (Fig. 9E), indicating a potential regulatory role of RA in ferroptosis in UC. Furthermore, marked expression changes were observed in several proteins associated with ROS suppression following treatment (Fig. 9F).
To further confirm that RA protects against mitochondrial dysfunction and ferroptosis via the ALOX15–VDAC1 axis, VDAC1 KD and ALOX15 overexpression rescue assays were performed in MCECs.
Western blot confirmed that VDAC1 was efficiently knocked down (Fig. 10A). In VDAC1-KD cells, LPS-induced MMP collapse and adenosine triphosphate (ATP) reduction were markedly alleviated, and RA could not further enhance these protective effects (Fig. 10B and C), indicating that VDAC1 is an essential downstream mediator of RA-mediated mitochondrial protection.
In the ALOX15 overexpression rescue assay, ALOX15 overexpression markedly reversed RA-induced VDAC1 down-regulation (Fig. 10D). Meanwhile, the recovery of MMP and elevation of ATP levels by RA were also abolished by ALOX15 overexpression (Fig. 10E and F).
To further verify whether RA regulates the ALOX15–VDAC1 axis by interfering with the endogenous protein–protein interaction between ALOX15 and VDAC1, co-immunoprecipitation (Co-IP) and immunofluorescence colocalization assays were performed in MCECs. Co-IP results revealed that LPS stimulation markedly enhanced the physical binding between endogenous ALOX15 and VDAC1, whereas RA pretreatment substantially reduced their coprecipitation levels. Consistently, immunofluorescence double staining combined with confocal microscopy showed elevated colocalization of ALOX15 and VDAC1 upon LPS stimulation. RA administration effectively decreased overlapping fluorescent signals, confirming that RA impairs formation of the ALOX15–VDAC1 protein complex. These results demonstrate that RA inhibits ferroptosis and enhances mitochondrial function in intestinal epithelial cells by regulating the ALOX15–VDAC1 signaling axis.
Previous studies have shown that VDAC1 promotes mitochondrial dysfunction during ferroptosis. Therefore, we sought to determine whether VDAC1 is a key component of this inflammatory injury model. We first evaluated the effect of RA on ATP generation. As shown in Fig. 11C, ATP levels were markedly reduced in LPS-treated cells compared with the control group. Given that MMP decline occurs at an early stage of mitochondrial damage, MMP was assessed using 5,5′,6,6′-tetrachloro-1,1′,3,3′-tetraethylbenzimidazolylcarbocyanine iodide (JC-1) staining. Compared with the control group, LPS treatment led to a marked reduction in the red/green fluorescence ratio (Fig. 11E), indicating mitochondrial depolarization.
Mitochondria are widely recognized as the primary source of intracellular ROS. Using the Mito-ROS assay, mitochondrial ROS levels were found to increase sharply following LPS treatment, an effect markedly counteracted by RA (Fig. 11A). Mito-TEMPO is a mitochondria-targeted superoxide dismutase mimetic that accumulates within the inner mitochondrial membrane via triphenylphosphonium-mediated uptake, thereby reducing LPS-induced mitochondrial ROS accumulation. Consistent with these findings, mitochondrial GPX4 protein expression was progressively reduced in LPS-damaged cells. Western blot analysis revealed that VDAC1 protein levels were markedly up-regulated in damaged cells, and both RA and Mito-TEMPO treatment partially reversed these changes (Fig. 11D). MitoTracker Red CMXRos staining further supported these observations (Fig. 11B). Taken together, these results demonstrate that both RA and Mito-TEMPO promote VDAC1 and its oligomer degradation, thereby contributing to mitochondrial homeostasis.
This study delineates a coherent mechanistic pathway through which RA alleviates UC. The ALOX15–VDAC1 axis emerges as a promising therapeutic target through which RA exerts its effects via ferroptosis suppression. This multidisciplinary work, integrating bioinformatics discovery, artificial intelligence-driven screening, and rigorous experimental validation, provides a framework linking target identification to functional outcomes.
Our investigation was motivated by the unmet clinical need in UC and the emerging pathogenic role of ferroptosis. Bioinformatics analysis of patient-derived data identified ALOX15 as a markedly up-regulated ferroptosis-related gene in UC with important diagnostic potential [2,5]. This finding provides a robust clinical and molecular rationale for targeting ALOX15 as a central therapeutic node, supporting a causal role in disease progression.
To identify compounds targeting this pathway, we employed a robust computational strategy. A machine learning model prioritized RA from a natural compound library [32–43], consistent with and extending our previously established protocol for ALOX15 inhibitor discovery [44]. Molecular docking and MD simulations structurally validated this prediction, revealing a stable and high-affinity binding mode in the active site of ALOX15. Critically, this predicted interaction was experimentally confirmed under physiological conditions via CETSA, DARTS, and Co-IP assays, demonstrating that their interaction improves the stability of ALOX15 protein. This integrated pipeline, from in silico prediction to biophysical and cellular validation, represents a robust framework for natural product-based drug discovery and establishes ALOX15 as a direct and actionable molecular target of RA.
Following candidate selection by the machine learning model, we experimentally confirmed that RA effectively inhibits ferroptosis and protects mitochondrial function in vitro and in vivo [45].
Having confirmed target engagement, we demonstrated the functional consequences of this interaction. RA inhibited ferroptosis in a concentration-dependent manner in vitro, restoring core anti-ferroptotic regulators (GPX4 and SLC7A11) and suppressing pro-ferroptotic drivers (ACSL4 and COX-2) [7,8,17]. Crucially, ALOX15 KD abolished the protective effects of RA, providing clear genetic evidence that ALOX15 is essential for RA-mediated cytoprotection, establishing a mechanistic relationship rather than a simple correlation. The in vivo relevance of this axis was further confirmed. RA ameliorated DSS-induced colitis by up-regulating tight junction proteins (ZO-1, Occludin, and Claudin-1), reducing disease activity scores, and restoring the integrity of the intestinal barrier. By preventing ferroptosis in epithelial cells [19], RA preserves the epithelial barrier, which represents a central component of UC pathogenesis.
Our research further elucidates the key downstream pathways that connect ALOX15 inhibition and mitochondrial homeostasis. Consistent with the proteomic findings, serum metabolomic analysis demonstrated that RA extensively remodels lipid metabolism pathways (including arachidonic acid metabolism, fatty acid oxidation, and linolenic acid metabolism), thereby linking its therapeutic effects to reduced lipid peroxidation and suppressed ferroptosis in colitis. Proteomic analysis and subsequent verification demonstrated that RA-mediated ALOX15 inhibition leads to VDAC1 down-regulation. As a mitochondrial gatekeeper, VDAC1 dysregulation is a known factor in ferroptotic damage [9,42]. Our findings demonstrate that RA treatment effectively reversed the up-regulation of VDAC1, while preserving MMP, suppressing mitochondrial ROS production, and restoring ATP levels. Furthermore, we provide direct genetic evidence through VDAC1 KD and ALOX15 overexpression rescue assays: VDAC1 KD abolished the mitochondrial protective effect of RA, while ALOX15 overexpression reversed RA-mediated regulation of VDAC1 and mitochondrial function, thereby confirming that the ALOX15–VDAC1 axis is an indispensable pathway for the anti-UC effect of RA. These findings support a mechanistic model in which upstream inhibition of ALOX15-driven lipid peroxidation reduces the direct damage to the mitochondrial membrane, thus preventing VDAC1 dysregulation and interrupting ROS-induced mitochondrial damage and further ROS production and metabolic collapse. The comparable effects of RA and the mitochondria-targeted antioxidant Mito-TEMPO further underscore the central importance of mitochondrial recovery in the mechanism of action of RA.
Further Co-IP and confocal colocalization data with quantitative Pearson coefficient calculation demonstrated that LPS promotes endogenous ALOX15–VDAC1 complex formation, while RA disrupts their physical interaction to down-regulate VDAC1 expression and ameliorate mitochondrial dysfunction, further consolidating the regulatory role of RA within the ALOX15–VDAC1 axis.
In summary, this study identifies the previously uncharacterized RA–ALOX15–VDAC1 axis as a mechanistic basis for RA’s therapeutic effects in UC (Fig. 12). Our work extends the traditional view of RA as a broad-spectrum antioxidant, repositioning it as a selective inhibitor of specific pro-ferroptotic enzymes. This discovery adds a new dimension to the growing body of evidence that natural products can regulate ferroptosis in human disease. Despite these findings, several limitations merit acknowledgement. Notably, this study lacks an independent RA-alone control group in in vivo experiments. Given the well-documented biosafety of RA, no appreciable histological injury or inflammatory abnormality was observed in healthy mice administered moderate doses of RA in previous studies. Nevertheless, the absence of a single RA-treated group for blank comparison is an undeniable limitation, and further research will supplement this group to clarify the direct biological effect of RA on normal intestinal tissue. The use of gene KD models has inherent limitations; future research using the inducible cell-specific knockout model will provide more accurate insights. The translational potential of RA also needs to be evaluated in human-relevant systems, including patient-derived tissues and organoid models, such as intestinal organoid models. The promising preclinical efficacy of RA necessitates a critical evaluation of its translational pathway. Pharmacokinetically, RA exhibits moderate oral bioavailability, with its glycosidic derivatives and colonic microbiota metabolism potentially enhancing its local concentration in the gut—a property potentially advantageous for UC treatment [46]. Its established safety profile, supported by a long history of dietary consumption and a high LD₅₀ in rodents [47–49], positions it favorably for further development, possibly as a nutraceutical or dietary adjunct. However, key challenges must be addressed to advance its therapeutic application. These include optimization of formulation strategies to improve its stability and targeted colonic delivery, rigorous long-term toxicity evaluation, and determination of whether the effective concentrations achieved in vitro and in murine models are translatable to human pathophysiology. Future work should prioritize validation in human intestinal organoids or primary epithelial cells from UC patients, and explore potential synergies with current first-line therapies such as mesalazine. Finally, future research should explore the crosstalk between RA-regulated ferroptosis and other cell death pathways.
In addition, several limitations related to the machine learning-based virtual screening framework should be acknowledged. Although a strict scaffold-disjoint validation strategy was employed to minimize optimistic bias and prevent structural leakage, the independent holdout subset ultimately contained a relatively limited number of Murcko scaffold series due to the highly imbalanced scaffold-size distribution of the curated ChEMBL dataset. Broader scaffold-level generalization was therefore primarily assessed through repeated scaffold-grouped OOF validation across the train-development subset rather than by the independent holdout evaluation alone. Moreover, the machine learning framework presented in this study should be interpreted as a candidate-prioritization and virtual-screening tool rather than definitive mechanistic evidence. While the model demonstrated stable discriminative performance under scaffold-disjoint evaluation settings, computational predictions alone cannot establish direct biochemical causality or guarantee generalization to entirely unseen chemical spaces. Additional large-scale datasets, external validation cohorts, and prospective experimental evaluations will be necessary to further improve the robustness and translational applicability of the screening framework. Despite these future directions, our findings highlight that the ALOX15–VDAC1 axis represents a promising therapeutic target for future drug development in inflammatory bowel disease.
We demonstrated that RA, identified via machine learning as an ALOX15 inhibitor, alleviates UC by modulating ALOX15 expression, down-regulating VDAC1, and thereby inhibiting VDAC1-mediated ferroptosis and mitochondrial dysfunction. Metabolomics profiling further demonstrated that RA regulates lipid metabolism and oxidative stress to alleviate colitis. Additionally, RA was confirmed to suppress ALOX15–VDAC1 complex formation, thereby blocking downstream mitochondrial damage and ferroptosis in intestinal epithelial cells. These results highlight the therapeutic potential of RA and identify the ALOX15–VDAC1 pathway as a promising target for the treatment of inflammatory bowel disease.
ALOX15 inhibitor records were collected from ChEMBL (n = 904). Molecular structures were standardized and validated using RDKit, and invalid SMILES were excluded. Duplicate molecular entries were aggregated by mean activity values, yielding 698 unique compounds for downstream modeling. IC50 values were converted to pIC50, and binary labels were assigned using pIC50 ≥ 5.2 as the activity threshold. Features were generated from Morgan fingerprints, MACCS fingerprints, and pretrained molecular representations. The pretrained representation was constructed using approximately 3 million drug-like molecules collected from the ZINC database without using any ALOX15 activity annotations. The pipeline proceeded as follows: Morgan fingerprints (radius = 2, 2,048 bits) and MACCS keys (167 bits) were computed for each pretraining molecule and concatenated into an initial feature vector; missing values were imputed using the median strategy; columns with variance below 0.16 were discarded; features were standardized to zero mean and unit variance using StandardScaler; and finally, IPCA was performed with a batch size of 4,096 to reduce the feature dimensionality to 128. The fitted imputer, variance selector, scaler, and IPCA projector were serialized as the pretraining preprocessor artifact.
During downstream ALOX15 classification modeling, each molecule’s Morgan and MACCS fingerprints were first computed, then passed through the imputation, filtering, and standardization steps of the pretrained preprocessor, and finally projected via the frozen IPCA model to yield a 128-dimensional pretrained representation. The final model input was formed by concatenating the filtered raw fingerprints with this 128-dimensional representation. The rationale for this hybrid feature design is twofold: Morgan and MACCS fingerprints capture local substructural and functional group information directly relevant to target-specific activity, whereas the IPCA-reduced representation encodes the global manifold structure of the chemical space learned from large-scale data, providing informative priors for scaffold regions underrepresented in the training set. Their combination is particularly advantageous for improving generalization under strict scaffold extrapolation.
To ensure rigorous generalization assessment and minimize optimistic bias caused by structural similarity, a strict scaffold-disjoint validation framework was adopted. Molecules were grouped according to Bemis–Murcko scaffolds, followed by scaffold-level train-development/test partitioning (9:1), in which scaffold groups were greedily assigned to the independent test subset until the predefined test size threshold was reached, ensuring zero scaffold overlap between the training and independent test subsets. Within the train-development subset, 5-fold scaffold GroupKFold cross-validation was further applied, such that molecules sharing the same scaffold were always assigned to the same fold. The detailed scaffold distribution across outer validation folds is summarized in Table S4. During benchmark comparison, all candidate classifiers were evaluated under the identical scaffold partitioning and outer-fold cross-validation protocol. For the finalized ExtraTrees framework, hyperparameter optimization was performed exclusively within nested inner-fold cross-validation restricted to each outer training subset, whereas outer validation folds and the independent scaffold-disjoint test set were reserved strictly for unbiased performance estimation and never participated in model selection, threshold optimization, or feature engineering. The optimal hyperparameters identified across outer folds are summarized in Table S7. The final classifier deployed for holdout testing and virtual screening was a 5-fold scaffold ensemble, formed by averaging positive-class probabilities from the 5 outer-fold ExtraTrees models trained under nested cross-validation.
Model performance was evaluated using both scaffold-grouped OOF predictions and an independent scaffold-disjoint test set, including ROC-AUC, PR-AUC, F1 score, balanced accuracy, precision, recall, and MCC. Decision thresholds were determined exclusively from OOF predictions and subsequently applied unchanged to the independent test set, thereby avoiding threshold-selection bias during holdout evaluation. To quantify statistical uncertainty, stratified bootstrap resampling (2,000 iterations) was further applied to estimate 95% CIs for ROC-AUC and PR-AUC. Collectively, this validation framework ensured that structurally related molecules were not simultaneously exposed across training, validation, and independent test subsets, thereby minimizing potential information leakage during model development and evaluation. In addition, Y-randomization analysis (499 label permutations) was performed to evaluate the likelihood of chance correlations. After model comparison, ExtraTrees was selected as the final classifier and used for prospective virtual screening of an external candidate library, followed by docking-based prioritization and experimental follow-up.
The bulk RNA-seq dataset GSE117993, related to pediatric UC, was retrieved from the Gene Expression Omnibus (GEO) database, comprising 55 healthy control and 43 UC patient samples. To ensure analytical rigor, only UC and control samples were retained; Crohn’s disease samples were excluded. Duplicate entries were removed. For each sample, the GEO accession, GEO Platform (GPL), and disease status were recorded.
DEGs associated with UC were identified from the merged expression matrix using the R limma package, with significance cutoffs of false discovery rate (FDR) < 0.05 and |log₂FC| > 1. The diagnostic capacity of ALOX15 was evaluated through ROC curve analysis and AUC calculation employing the pROC package.
Using the median expression of ALOX15 as a cutoff, samples were categorized into high- and low-expression groups. Differential enrichment across 114 metabolic pathways was assessed by integrating the GSVA algorithm with the limma package (|t| > 0 defined as the significance criterion). Simultaneously, GSEA based on KEGG pathways was applied to DEGs (P < 0.05) using clusterProfiler (P < 0.05).
Proteomics data were processed using DIA-NN (v1.8.1) with a library-free search against UniProtKB, with quantification via the MaxLFQ algorithm. The limma package facilitated DEP identification (P < 0.05, |log₂FC| > 2). Functional enrichment (GO/KEGG) was performed with clusterProfiler, and protein feature diagrams were generated using drawProteins and ggplot2.
Serum metabolic profile data were obtained from C57 mice. The metabolic data were log-transformed and mean-centered using the SIMCA software (V16.0.2), followed by PCA and orthogonal partial least squares discriminant analysis. Metabolites with P < 0.05 and a VIP greater than 1 were selected as differential metabolites. Visualization was performed using R (v4.5.0). Subsequently, the differential metabolites were subjected to KEGG and SMPDB enrichment analysis through MetaboAnalyst.
C57BL/6 mice (n = 5 per group) were acclimated for 2 weeks prior to random allocation into 3 groups: a control group, a DSS group (DSS sourced from MP Biomedicals, UK), and a DSS+RA group (RA obtained from Shaanxi Tianduoli Biological Technology Co., Ltd., China). The administered dose of RA (100 mg/kg) was determined based on previously published pharmacological studies focusing on RA. This moderate dose within the commonly used effective range (50 to 150 mg/kg) was selected to ensure stable anti-inflammatory and antioxidant activity with negligible systemic toxicity, consistent with prior in vivo investigations of RA in murine inflammatory models [29,50,51]. The DSS+RA group received daily intragastric administration of RA (100 mg/kg) throughout the experimental period (days 0 to 10). Colitis was induced in both the DSS and DSS+RA groups by providing 3% (w/v) DSS in drinking water from day 3 to day 10. Body weight and DAI were recorded daily. On day 11, mice were euthanized, blood was collected via the retro-orbital plexus, and colons were harvested. Colon segments were either fixed in 4% paraformaldehyde (PFA) or snap-frozen. Serum, prepared by centrifuging blood at 12,000 rpm for 10 min (4 °C), was aliquoted and stored at −80 °C.
H&E staining was performed on colon tissues fixed in 4% PFA (Macklin, China). Subsequent processing included paraffin embedding, sectioning, deparaffinization, rehydration and staining according to standard protocols.
For immunohistochemistry (IHC) staining, tissue sections were deparaffinized, rehydrated, and subjected to heat-mediated antigen retrieval. To block endogenous peroxidase activity, sections were incubated with 3% H₂O₂ (Nanyue, China) for 15 min. Nonspecific binding was minimized by incubating sections in 3% goat serum (Sangon Biotech, China) for 30 min at room temperature. Primary antibodies targeting ZO-1, Occludin, Claudin-1, MUC-2, and ALOX15 (ABclonal, China), diluted in 5% bovine serum albumin (BSA), were then applied overnight at 4 °C. Following thorough washing, horseradish peroxidase (HRP)-conjugated secondary antibodies (1:5,000; CST, China) were applied for 30 min. Color development was carried out with DAB solution, nuclei were counterstained with hematoxylin, and coverslips were mounted using neutral resin (Solarbio, China).
The mouse colon epithelial cell line MCEC and rat intestinal epithelial cell line IEC-6 (Fenghui Biotechnology Co., Ltd., China) were cultured in Dulbecco's Modified Eagle Medium (Gibco, USA) supplemented with 10% fetal bovine serum at 37 °C under 5% CO₂. For experiments, cells were seeded in 6-well plates (3 × 106 cells/well) and treated for 24 h with RA (0, 10, 20, or 30 μM) or LPS (1 μg/ml; Sigma-Aldrich, USA). The LPS+RA group received 2 h of RA pretreatment prior to LPS exposure. For ALOX15 inhibition experiments, cells were pretreated with PD146176 (10 μM) for 1 h before LPS stimulation.
MCECs were seeded in 96-well plates (4 × 103 cells/well). Following the treatment period, 10 μl of CCK-8 reagent (NCM Biotech, China) was added per well and incubated for 3 h. Absorbance at 450 nm was measured to determine cell viability.
Cells were loaded with 1 μM FerroOrange (Dojindo, Japan) in serum-free medium for 20 min at 37 °C and washed twice with phosphate-buffered saline (PBS) (Sangon Biotech, China), and fluorescence intensity was quantified using a Cytation 5 microplate reader (BioTek, USA).
Cells were stained with 10 μM C11-BODIPY581/591 (ABclonal, China) for 30 min at 37 °C, washed with PBS, and analyzed by flow cytometry (DXP Athena, Cytek Biosciences, USA).
Intracellular GSH levels were quantified using a DTNB colorimetric assay kit (Beyotime, cat. no. S0053), and MDA content was measured via a thiobarbituric acid method kit (Beyotime, S0131S), strictly following the manufacturer’s protocols.
Cells were lysed in radioimmunoprecipitation assay buffer supplemented with protease inhibitors and phenylmethylsulfonyl fluoride (Solarbio, China). After protein quantification, equal amounts of lysate were separated by SDS-PAGE (Sangon Biotech, China) and transferred to nitrocellulose membranes. Membranes were blocked with 5% BSA (Solarbio, China) for 30 to 60 min, then probed overnight at 4 °C with primary antibodies against GPX4, FTH1, SLC7A11 (CST, China), COX-2, ACSL4, ALOX15, Occludin, Claudin-1, and ZO-1 (ABclonal, China). After incubation with HRP-conjugated secondary antibodies (1:5,000; CST, China), bands were visualized using BeyoECL Moon reagent (Beyotime Biotech, China). GAPDH and β-actin served as loading controls.
MCEC suspensions were subjected to a thermal gradient (40 to 80 °C, 10 min), chilled on ice (5 min), and centrifuged (12,000×g, 10 min, 4 °C). The resulting supernatants were subjected to Western blot analysis.
MCEC lysates were incubated with RA for 1 h at room temperature, followed by proteolytic digestion with Pronase (Roche, Switzerland) at an enzyme-to-protein ratio of 1:300 for 40 min at 37 °C. Reactions were terminated by boiling for 3 min prior to Western blot analysis.
Following fixation (4% formaldehyde) and permeabilization (0.5% Triton X-100), cells were incubated with primary antibodies at 4 °C overnight. After washing, fluorescent secondary antibodies and 4′,6-diamidino-2-phenylindole (DAPI) were applied for signal detection and nuclear counterstaining, respectively. Images were acquired using an FV3000 confocal microscope (Olympus, Japan).
Cell lysates were incubated with an anti-ALOX15 antibody. The immunoprecipitated complex was collected, washed, and analyzed by Western blot.
MMP was assessed using JC-1 staining according to the manufacturer’s instructions, with subsequent imaging by fluorescence microscopy.
Compound libraries were prepared using LigPrep (Schrödinger) at pH 7.0 ± 2.0. The ALOX15 protein structure (AlphaFold ID: P16050) was prepared using the Protein Preparation Wizard. Docking was performed in extra-precision mode with a defined receptor grid (center coordinates: 4.03, 6.16, 4.3 Å; grid dimensions: 36 Å) enclosing the active site.
The ALOX15–RA complex was subjected to a 100-ns MD simulation using the AMBER99SB-ILDN force field. Ligand parameters were generated with Sobtop/GAFF. The system was solvated, neutralized, energy-minimized, and equilibrated before the production run. Trajectories were analyzed for RMSD, RMSF, Rg, and FEL.
MCECs were transfected with VDAC1 siRNA (5′-UAAAUUCCAAUCCAUUCUC-3′) or negative control siRNA for 48 h. The KD efficiency was verified by Western blot. Transfected cells were allocated into LPS (1 μg/ml, 24 h) and LPS + RA (30 μM, 2 h pretreatment) groups, and MMP and ATP levels subsequently measured.
MCECs were transfected with ALOX15 overexpression plasmid or empty vector for 48 h. Overexpression efficiency was confirmed by Western blot. Cells were allocated into empty vector + LPS, empty vector + LPS + RA, and ALOX15 overexpression + LPS + RA groups. LPS (1 μg/ml) and RA (30 μM) were used for 24 h treatment with 2 h pretreatment. VDAC1 expression, MMP, and ATP levels were measured.
All quantitative data are expressed as mean ± standard deviation (SD). Statistical analyses were performed using GraphPad Prism (version 9.5). The normality of data distribution was assessed using the Shapiro–Wilk test. For comparisons between 2 groups, an unpaired 2-tailed Student t test was applied. For comparisons among 3 or more groups, one-way analysis of variance (ANOVA) followed by Tukey’s post hoc test for multiple comparisons was used. For data involving repeated measurements over time (e.g., body weight and disease activity index), 2-way repeated-measures ANOVA followed by Šidák’s multiple comparisons test was employed. A P value < 0.05 was considered statistically significant. The specific statistical tests used for each dataset are detailed in the corresponding figure legends.
All animal procedures were reviewed and approved by the Institutional Animal Care and Use Committee of Guangdong Medical University (approval no. GDY2502901). Experiments involving animals were carried out in strict accordance with national regulations and the guidelines established by the laboratory animal center at Guangdong Medical University.
1.
Che S, Qin B, Wu K, Zhu M, Hu H, Peng C, Wang Z, Yin Y, Xia Y, Wu M.EGCG drives gut microbial remodeling-induced epithelial GPR43 activation to lessen Th1 polarization in colitis. Redox Biol. 2024;75:Article 103291.
2.
Buie MJ, Quan J, Windsor JW, Coward S, Hansen TM, King JA, Kotze PG, Gearry RB, Ng SC, Mak JW, et al.Global hospitalization trends for Crohn’s disease and ulcerative colitis in the 21st century: A systematic review with temporal analyses. Clin Gastroenterol Hepatol. 2023;21(9):2211–2221.
3.
Zhou Y, Yang M, Yan X, Zhang L, Lu N, Ma Y, Zhang Y, Cui M, Zhang M, Zhang M.Oral nanotherapeutics of andrographolide/carbon monoxide donor for synergistically anti-inflammatory and pro-resolving treatment of ulcerative colitis. ACS Appl Mater Interfaces. 2023;15(30):36061–36075.
4.
McCoubrey LE, Seegobin N, Sangfuang N, Moens F, Duyvejonck H, Declerck E, Dierick A, Marzorati M, Basit AW.The colon targeting efficacies of mesalazine medications and their impacts on the gut microbiome. J Control Release. 2024;369:630–641.
5.
Cui Y, Chen Y, Li H, Zhang W, Wang X, Xia M, Gan N, Zhou Y, Li M, Zhang H, et al.PCSK9 promotes atherosclerotic plaque instability by inducing VSMC Ferroptosis through the YAP1–NUPR1 axis. Research. 2025;8:0922.
6.
Dixon SJ, Lemberg KM, Lamprecht MR, Skouta R, Zaitsev EM, Gleason CE, Patel DN, Bauer AJ, Cantley AM, Yang WS, et al.Ferroptosis: An iron-dependent form of nonapoptotic cell death. Cell. 2012;149(5):1060–1072.
7.
Dixon SJ, Stockwell BR.The hallmarks of ferroptosis. Annu Rev Cancer Biol. 2019;3:35–54.
8.
Li J, Cao F, Yin HL, Huang ZJ, Lin ZT, Mao N, Sun B, Wang G.Ferroptosis: Past, present and future. Cell Death Dis. 2020;11(2):88.
9.
Yao MY, Liu T, Zhang L, Wang MJ, Yang Y, Gao J.Role of ferroptosis in neurological diseases. Neurosci Lett. 2021;747:Article 135614.
10.
Hu M, Jiang W, Ye C, Hu T, Yu Q, Meng M, Sun L, Liang J, Chen Y.Honokiol attenuates high glucose-induced peripheral neuropathy via inhibiting ferroptosis and activating AMPK/SIRT1/PGC-1α pathway in Schwann cells. Phytother Res. 2023;37(12):5787–5802.
11.
Zhao QX, Yan SB, Wang F, Li XX, Shang GK, Zheng ZJ, Xiao J, Lin ZW, Li CB, Ji XP.STING deficiency alleviates ferroptosis through FPN1 stabilization in diabetic kidney disease. Biochem Pharmacol. 2024;222:Article 116102.
12.
Du WW, Zhou C, Yang H, Wen S, Chen Y, Chen EX, Yang XH, Li F, Du KY, Yuan H, et al.Aggravated ulcerative colitis via circNlgn-mediated suppression of nuclear actin polymerization. Research. 2024;7:0441.
13.
Fang X, Liu H, Du Y, Jiang L, Gao F, Wang Z, Chi Z, Shi B, Zhao X.Bacillus siamensis targeted screening from highly colitis-resistant pigs can alleviate ulcerative colitis in mice. Research. 2024;7:0415.
14.
Guo Y, Du J, Xiao C, Xiang P, Deng Y, Hei Z, Li X.Inhibition of ferroptosis-like cell death attenuates neuropathic pain reactions induced by peripheral nerve injury in rats. Eur J Pain. 2021;25(6):1227–1240.
15.
Wang L, Wang FS.The role of innate lymphoid cells in immune-mediated liver diseases. Signal Transduct Target Ther. 2022;7(1):237.
16.
Li C, Gao P, Zhuang F, Wang T, Wang Z, Wu G, Zhou Z, Xie H, Xie D, Zhao D, et al.Inhibition of ALOX12–12-HETE alleviates lung ischemia–reperfusion injury by reducing endothelial ferroptosis-mediated neutrophil extracellular trap formation. Research. 2024;7:0473.
17.
Doll S, Proneth B, Tyurina YY, Panzilius E, Kobayashi S, Ingold I, Irmler M, Beckers J, Aichler M, Walch A, et al.ACSL4 dictates ferroptosis sensitivity by shaping cellular lipid composition. Nat Chem Biol. 2017;13(1):91–98.
18.
Zhang K, Yang X, Wang Y, Yu Y, Huang N, Li G, Li X, Wu JC, Yang S.Artificial intelligence in drug development. Nat Med. 2025;31(1):45–59.
19.
Li M, Yu X, Liu Y, Ouyang S, Wu L, Chen X, Yu H, Chen H, Lian S, Li Z, et al.KRAS/ABHD17C/ALOX15B axis promotes pancreatic cancer progression via ferroptosis evasion. Adv Sci. 2025;12(35):Article e04470.
20.
Zhang W, Li X, Wang Y, Jiang Y, Lu Y, Tu P.Ethanol extract of Corydalis decumbens exerts neuroprotective effects against cerebral ischemia-reperfusion injury by targeting ALOX15 to suppress ferroptosis. Phytomedicine. 2025;149:Article 157546.
21.
Noor S, Mohammad T, Rub MA, Raza A, Azum N, Yadav DK, Hassan MI, Asiri AM.Biomedical features and therapeutic potential of rosmarinic acid. Arch Pharm Res. 2022;45(4):205–228.
22.
Deng S, Cao H, Li T, Wang X, Meng J, Zeng T, Zhang D, Zhang S, Wang G, Liu R, et al.Lachnospiraceae-bacterium alleviates ischemia-reperfusion injury in steatotic donor liver by inhibiting ferroptosis via the Foxo3-Alox15 signaling pathway. Gut Microbes. 2025;17(1):2460543.
23.
Han YH, Kee JY, Hong SH.Rosmarinic acid activates AMPK to inhibit metastasis of colorectal cancer. Front Pharmacol. 2018;5(9):68.
24.
Ye H, Wu L, Liu YM, Zhang JX, Hu HT, Dong ML, Ren J.Wogonin attenuates septic cardiomyopathy by suppressing ALOX15-mediated ferroptosis. Acta Pharmacol Sin. 2025;46(9):2407–2422.
25.
Esmeeta A, Adhikary S, Dharshnaa V, Swarnamughi P, Ummul Maqsummiya Z, Banerjee A, Pathak S, Duttaroy AK.Plant-derived bioactive compounds in colon cancer treatment: An updated review. Biomed Pharmacother. 2022;153:Article 113384.
26.
Konstantinou EK, Panagiotopoulos AA, Argyri K, Panoutsopoulos GI, Dimitriou M, Gioxari A.Molecular pathways of rosmarinic acid anticancer activity in triple-negative breast cancer cells: A literature review. Nutrients. 2023;16(1):2.
27.
Yang M, Wang T, Shao J, Ran X, Xiao R, Zhao R, Wu C, Ji M, Tian W, Sun H, et al.(+)-JQ-1 alleviates cardiac injury in myocardial infarction by inhibiting ferroptosis through the NAMPT/SIRT1 pathway. Cell Death Dis. 2025;16(1):548.
28.
Bouammali H, Zraibi L, Ziani I, Merzouki M, Bourassi L, Fraj E, Challioui A, Azzaoui K, Sabbahi R, Hammouti B, et al.Rosmarinic acid as a potential source of natural antioxidants and anticancer agents: A molecular docking study. Plants. 2023;13(1):89.
29.
Wang GY, Chen SY, Chen YY, Hong CJ, Hsu YH, Yen GC.Protective effect of rosmarinic acid-rich Trichodesma khasianum Clarke leaves against ethanol-induced gastric mucosal injury in vitro and in vivo. Phytomedicine. 2021;80:Article 153382.
30.
Liu J, Liu Y, Huang C, He C, Yang T, Ren R, Xin Z, Wang X.Quercetin-driven Akkermansia muciniphila alleviates obesity by modulating bile acid metabolism via an ILA/m6A/CYP8B1 signaling. Adv Sci. 2025;12(12):Article e2412865.
31.
Zhu L, Tan Q, Wang Y, Hong L, Chen C, Kong L, Luo J.Artemisitene triggers calcium-dependent ferroptosis by disrupting the LSH-EWSR1 interaction in colorectal cancer. Redox Biol. 2026;89:Article 103950.
32.
Jia J, Liu Y, Wang D, Pan Z, Zheng Q, Lu J, Liang C, Li D.Ulcerative colitis: Signaling pathways, therapeutic targets and interventional strategies. Signal Transduct Target Ther. 2026;11(1):51.
33.
Vamathevan J, Clark D, Czodrowski P, Dunham I, Ferran E, Lee G, Li B, Madabhushi A, Shah P, Spitzer M, et al.Applications of machine learning in drug discovery and development. Nat Rev Drug Discov. 2019;18(6):463–477.
34.
Rodríguez-Pérez R, Miljković F, Bajorath J.Machine learning in chemoinformatics and medicinal chemistry. Annu Rev Biomed Data Sci. 2022;10(5):43–65.
35.
Yuan M, Zhou X, Xu W, Liu J, Han T, Huang Y, Liao F, Chang J, Miao C.Wogonin ameliorates osteoarthritis by targeting the STAT3/PIM1 axis to attenuate chondrocyte senescence. Phytomedicine. 2026;153:Article 157986.
36.
Hou Y, Wang Z, Wang W, Tang Q, Cai Y, Yu S, Wang J, Yan X, Wang G, Lobie PE, et al.AI-identified CD133-targeting natural compounds demonstrate differential anti-tumor effects and mechanisms in pan-cancer models. EMBO Mol Med. 2025;17(11):2932–2965.
37.
Wang W, Wang X, Zhou ZD, Deng XB, Luo XW, Wang Q, Zhang GY, Zhang Z, Li ZS, He XJ, et al.7,8-Dihydroxyflavone alleviates ulcerative colitis by binding keap1 to activate an slc7a11-dependent inhibition of NF-κB p65 signaling in macrophages. J Adv Res. 2026;S2090-1232(26):00103-7.
38.
Zhang J, Lin X, Song H, Xie Y, Hou S, Huang S, Du X, Li H, Han Y, Liang J, et al.Dendrobium officinale polysaccharide inhibits M1 macrophage polarization via activating SENP1-SIRT3 signaling and alleviates ulcerative colitis. Food Res Int. 2025;221(Pt 4):Article 117582.
39.
Yang WS, Stockwell BR.Ferroptosis: Death by lipid peroxidation. Trends Cell Biol. 2016;26(3):165–176.
40.
Stockwell BR, Jiang X.The chemistry and biology of ferroptosis. Cell Chem Biol. 2020;27(4):365–375.
41.
Sun Q, Tu K, Xu Q, Yan L, Yang S, Wang J, Lv L, Liu H, Cai L.Berberine suppresses colorectal cancer progression by inducing ferroptosis-mediated energy metabolism disorders. J Adv Res. 2025;S2090-1232(25):00820-3.
42.
Liao Y, Cao P, Luo L.Development of novel ALOX15 inhibitors combining dual machine learning filtering and fragment substitution optimisation approaches, molecular docking and dynamic simulation methods. J Enzyme Inhib Med Chem. 2024;39(1):2301756.
43.
Deng L, Tian W, Luo L.Application of natural products in regulating ferroptosis in human diseases. Phytomedicine. 2024;128:Article 155384.
44.
Guo X, Zhao X, Lu X, Zhao L, Zeng Q, Chen F, Zhang Z, Xu M, Feng S, Fan T, et al.A deep learning-driven discovery of berberine derivatives as novel antibacterial against multidrug-resistant Helicobacter pylori. Signal Transduct Target Ther. 2024;9(1):183.
45.
Huang Y, Chen H, Wu J, Huang X, Zhang Y, Ren A, Wu X, Huang C, Xie W, Zu D, et al.Machine learning-driven discovery of aniline phenylacetate derivatives as LasR-targeted biofilm inhibitors against Pseudomonas aeruginosa infections. Acta Pharm Sin B. 2026;doi.org/10.1016/j.apsb.2026.03.048
46.
Rousseau PN, Bazin PL, Steele CJ.Pontine functional connectivity gradients. Cerebellum. 2025;25(1):1.
47.
He Y, Wu S, Li J, Chen S, Chen S, Zhang Z, He B, Hong Y, Sun C, Kai G.Artificial intelligence in traditional Chinese medicine: Unraveling herbal medicine’s mechanisms. Research. 2026;9:1224.
48.
Huang J, Wang J, Chai Z, Wang Z, Wang H, Su Y.The gut microbiota in neurodegenerative diseases: Current landscape and the potential of multiomics. Nutrients. 2023;15(15):3366.
49.
Allen HK, Levine UY, Looft T, Bandrick M, Casey TA.Treatment, promotion, commotion: Antibiotic alternatives in food-producing animals. Trends Microbiol. 2013;21(3):114–119.
50.
Ruiz-Malagón AJ, Rodríguez-Sojo MJ, Hidalgo-García L, Molina-Tijeras JA, García F, Pischel I, Romero M, Duarte J, Diez-Echave P, Rodríguez-Cabezas ME, et al.The antioxidant activity of Thymus serpyllum extract protects against the inflammatory state and modulates gut dysbiosis in diet-induced obesity in mice. Antioxidants. 2022;11(6):1073.
51.
Zhou B, Wang L, Yang S, Liang Y, Zhang Y, Pan X, Li J.Rosmarinic acid treatment protects against lethal H1N1 virus-mediated inflammation and lung injury by promoting activation of the h-PGDS-PGD2-HO-1 signal axis. Chin Med. 2023;18(1):139.
Year 2026 volume 9 Issue 7
PDF
101
56
Cite this Article
BibTeX
Article Info
doi: 10.34133/research.1352
  • Receive Date:2026-03-11
  • Online Date:2026-10-15
  • Published:2026-07-20
Article Data
Affiliations
History
  • Received:2026-03-11
  • Revised:2026-06-16
  • Accepted:2026-06-16
Affiliations
    1School of Ocean and Tropical Medicine, Guangdong Medical University, Zhanjiang, Guangdong 524023, China.
    2Department of Interventional Radiology, Zhongshan Hospital, Fudan University, Shanghai 200032, China.

Corresponding:

* Address correspondence to: (B.T.);
* Address correspondence to: (L.L.)
References
Share
https://castjournals.cast.org.cn/joweb/research/EN/10.34133/research.1352
Share to
QR

Scan QR to access full text

Cite this article
BibTeX
Citations
表12种不同金属材料的力学参数
科
Family
属数
Number of
genus
种数
Number of
species
占总种数比例
Percentage of
total species (%)
属
Genus
种数
Number of
species
占总种数比例
Percentage of total
species (%)
鹅膏菌科Amanitaceae 2 11 5.26 鹅膏菌属 Amanita 10 4.78
小菇科 Mycenaceae 2 12 5.74 丝盖伞属 Inocybe 5 2.39
多孔菌科 Polyporaceae 8 14 6.70 蜡蘑属 Laccaria 5 2.39
红菇科 Russulaceae 3 23 11.00 小皮伞属 Marasmius 6 2.87
小菇属 Mycena 11 5.26
光柄菇属 Pluteus 5 2.39
红菇属 Russula 17 8.13
栓菌属 Trametes 5 2.39
关闭全屏
  • BibTeX
  • EndNote
  • RefWorks
  • TxT