厌倦了常规毒性机制研究？2025年环境健康领域期刊《Ecotoxicology and Environmental Safety》（IF = 6.1）最新研究给出完美示范：
王牌研究策略：计算生物学→ 机器学习 → 分子模拟 → 细胞实验
四位一体研究框架，全方位验证PANK2靶点
亮点一：多算法机器学习靶点筛选
采用LASSO回归、SVM-RFE和随机森林三重机器学习算法，从387个PFAS潜在作用靶点中精准筛选出关键基因PANK2，跨算法验证确保靶点可靠性！
亮点二：分子模拟揭示结合机制
分子对接显示PFOA/PFOS与PANK2结合能分别达-8.4/-9.6 kcal/mol！100ns分子动力学模拟表明复合物结构稳定，关键结合残基GLU338、TYR454、ARG407参与相互作用。
亮点三：体外实验验证功能影响
10μM PFOA/PFOS处理骨髓间充质干细胞：
➠ 显著抑制成骨分化（ALP、RUNX2、OCN基因下调）；
➠ 明显促进成脂分化（PPARγ、C/EBPα基因上调）；
➠ PANK2表达显著下降（mRNA和蛋白水平双重验证）。
亮点四：机制深度挖掘与临床意义
GSEA分析提示chemokine signaling通路参与调控；
ceRNA网络构建揭示miRNA-lncRNA上游调控轴；
为环境污染物所致骨质疏松的防治提供新靶点。
方法可迁移性强：
适用于环境毒理学、骨代谢疾病、计算生物学等多个研究方向；
完整的"计算预测-实验验证"研究范式，轻松复现6+SCI思路！
Highlights
1、网络毒理学和机器学习将PANK2确定为中枢骨毒性靶标
2、分子对接和动力学揭示了全氟辛酸/全氟辛烷磺酸与PANK2的强结合
3、全氟辛酸/全氟辛烷磺酸降低PANK2的表达，损害BMSCs中的骨脂肪生成平衡
4、GSEA和ceRNA分析揭示了PANK2与骨代谢之间的联系机制
背景介绍
近年来，全氟辛酸（PFOA）和全氟辛烷磺酸（PFOS）作为典型的全氟及多氟烷基物质（PFAS），因其在环境中的持久性、生物累积性以及对人体健康的潜在危害而受到广泛关注。研究表明，PFAS暴露与多种不良健康效应相关，包括代谢紊乱、内分泌干扰和免疫失调等。尤其值得注意的是，越来越多的流行病学证据提示PFOA和PFOS可能与骨密度降低、骨折风险增加以及骨质疏松等骨骼健康问题密切相关。然而，其具体的骨毒性机制尚不明确。
本研究通过整合网络毒理学、分子对接与动力学模拟以及体外细胞实验，系统探讨了PFOA和PFOS对骨代谢的毒性作用机制。研究发现，泛酸激酶2（PANK2）作为一个关键靶点，在PFAS暴露后表达显著下调，进而破坏骨髓间充质干细胞（BMSCs）的成骨-成脂分化平衡，最终导致骨质疏松的发生。该研究首次揭示了PANK2在PFAS诱导骨代谢紊乱中的核心作用，为环境污染物所致骨骼疾病的机制研究和靶向干预提供了新的理论依据。
主要研究结果
（一）PFOA/PFOS致骨质疏松作用靶点的鉴定
为了阐明PFAS暴露与骨质疏松症之间的潜在分子联系，研究首先确定了化合物相关靶标和疾病相关基因。从ChEMBL、STITCH和SwissTargetPrediction数据库中获得了387个PFOA和383个PFOS相关靶点，发现PFOA和PFOS在生物体系中的靶点几乎相同（图1A）。可能是由于它们的整体形状、分子大小和极性特征相似，使它们能够在生物系统中与相同的蛋白质发生相互作用。因此，合并PFOA和PFOS的靶标用于后续分析。
接下来，通过分析GEO数据库中的基因表达谱，获得了与骨质疏松症相关的DEGs，共鉴定出2785个DEGs，其中1475个上调，1310个下调（图1B）。通过将PFAS目标集与糖尿病相关的DEG交叉，鉴定了14个基因可能是连接PFAS暴露与骨质疏松发病机制的关键介质（图1C）。这14个候选基因的表达模式通过热图分析可视化，显示出在阿尔茨海默病与对照样品中不同的表达谱，这进一步支持了它们在疾病进展中的潜在相关性（图1D）。
图1. PFOA/PFOS诱导的骨质疏松症的靶点鉴定
（二）PFOA/PFOS在骨质疏松症中显著富集的代谢和脂质相关途径
通过KEGG和GO富集分析探索了14个PFAS-骨质疏松重叠基因的潜在生物学作用，通过WebGestalt进行的KEGG途径富集分析显示，这些基因主要参与脂肪酸的生物合成、降解和代谢，以及泛酸和辅酶A（CoA）的生物合成、铁凋亡和过氧化物酶体增殖物激活受体（PPAR）信号途径这些途径与脂质代谢和骨稳态密切相关，GO富集分析进一步揭示了PFAS基因在生物过程水平上主要参与代谢过程、生物调节和发育过程就细胞成分而言，它们主要定位于含蛋白质的复合物、染色体和囊泡，表明它们在细胞内能量和脂质处理中的潜在作用它们的分子功能主要与蛋白质结合、离子结合和核酸结合有关，这可能反映了它们参与信号转导和代谢调节（图2D）。
图2. PFOA/PFOS致骨质疏松症相关候选靶点的富集分析
（三）机器学习算法识别PANK 2为PFOA/PFOS诱导的骨质疏松症中的枢纽基因
为了在14个候选基因中优先考虑关键基因，采用了三种机器学习算法：LASSO回归，SVM-RFE和RF。LASSO回归中的系数轨迹图显示，随着λ的增加，基因系数逐渐缩小，交叉验证确定了使误差最小化的最佳λ（图3A-B）。随着所选特征数量的增加，SVMRFE模型的分类准确率和交叉验证误差都有所提高，在6个特征上两个指标都达到最佳值（图3C-D）。RF算法进一步使用准确度的平均下降和基尼指数来评估特征重要性。随着树的数量增加，模型误差稳定，表明性能稳健选择重要性得分> 1的基因作为最佳特征（图3E）。通过将三个模型的结果交叉，PANK 2是在所有算法中一致鉴定的唯一基因，强调其作为介导全氟辛酸和全氟辛烷磺酸骨毒性效应的强大枢纽基因的潜力，不同特征选择方法之间的这种收敛突出了PANK 2用于进一步功能研究的可靠性和相关性。
图3. 机器学习识别与PFOA/PFOS诱导的骨质疏松症相关的枢纽基因
（四）PFOA/PFOS与PANK形成稳定的相互作用
结果表明，PFOA和PFOS与PANK 2的结合能分别为-8.4 kcal/mol和-9.6 kcal/mol，表明它们与PANK 2的结合潜力很高，并通过多种非共价力相互作用，包括疏水作用、货车范德华力、氢键值得注意的是，全氟辛酸与PANK 2的ATP结合口袋内的残基形成稳定的相互作用，而全氟辛烷磺酸表现出更强的结合特征，与类似的功能位点结合。
为了评估小分子与目标蛋白的结合稳定性，对PANK 2-PFOA和PANK 2-PFOS复合物进行了分子动力学模拟，均方根偏差（RMSD）表明两种复合物在5 ns内达到平衡并保持稳定，波动在2.4 nm左右，表明有利的构象稳定性回转半径（Rg）和溶剂可及表面积（SASA）在模拟过程中仅显示出微小的波动，表明复合物经历了轻微的构象变化而没有显著的结构破坏氢键分析显示PANK 2-PFOA复合物保持0-3个氢键，通常约为1，而PANK 2-PFOS复合物形成0-6个氢键，平均约为3均方根波动（RMSF）分析显示，大多数氨基酸残基表现出低柔性，其值主要低于4 μm，表明相对刚性和稳定的复合物结构为了进一步评估构象稳定性，构建能量景观图。结果显示，PANK 2-PFOA和PANK 2-PFOS复合物在模拟期间占据低能稳定构象。总体而言，结果表明PFOA和PFOS均与PANK 2稳定结合，PFOS表现出略强的相互作用。
图4. 分子对接和分子动力学
（五）PFOA和PFOS可能通过下调PANK 2而破坏BMSCs的成骨-成脂平衡
为了评估PFOA和PFOS对BMSC的细胞毒性效应，基于先前发表的研究，用不同浓度处理细胞，在24 h时未观察到显著变化。48 h后，100μM全氟辛酸或全氟辛烷磺酸导致细胞存活率降低。到72 h时，10 μM和100 μM均显著降低了存活率。为了尽量减少细胞毒性干扰，同时保持可观察到的生物效应，在随后的实验中，将BMSC用10 μM PFOA或PFOS处理48 h。进一步分析骨髓间充质干细胞的多向分化潜能，发现全氟辛酸和全氟辛烷磺酸暴露会损害成骨分化，茜素红S染色显示钙沉积减少（图5A）；相反，成脂分化增强，油红O染色显示细胞内脂滴积累增加（图5 B）。
qRT-PCR分析证实了这些观察结果水平。成骨标志物的表达，包括与侏儒相关的标志物在PFOA或PFOS处理后，转录因子2（RUNX2）、碱性磷酸酶（ALP）、骨钙素（OCN）和骨基质（SP7）的表达显著下调（图5C）。RUNX 2是成骨的主要调节因子，ALP和OCN参与基质成熟和矿物质沉积，而SP 7是终末成骨细胞分化的关键。相反，关键的脂肪形成基因，如过氧化物酶体增殖物激活受体γ（PPARγ），CCAAT/增强子结合蛋白α（C/EBPα）和脂肪酸结合蛋白4（FABP 4）上调，表明PFOA和PFOS暴露破坏了BMSC中的骨-脂肪生成平衡（图5D）。PFOA和PFOS均能显著降低BMSCs中PANK 2的表达，提示PANK 2可能参与PFOS诱导的骨代谢紊乱与转录数据一致，蛋白质印迹显示在处理组中PANK 2蛋白表达降低，如免疫印迹结果和光密度定量所示（图5 F-G）。
图5. qRT-PCR验证实验
（六）GSEA和ceRNA网络揭示PANK 2的潜在调控机制
GSEA分析结果显示，PANK 2低表达时，在趋化因子信号通路、长时程抑制、嗅觉传导等信号通路中明显富集（图6A）。为了进一步研究PANK 2的上游调节机制，预测了可能靶向PANK 2的miRNAs，并鉴定了可能靶向这些miRNAs的相应lncRNAs所得的ceRNA网络表明PANK 2表达可以通过特异性miRNA-lncRNA相互作用来调节，强调了转录后调节在全氟辛酸/全氟辛烷磺酸暴露下介导PANK 2表达的潜在作用（图6B）。
图6. GSEA和ceRNA网络
总结
本研究通过整合网络毒理学、机器学习、分子模拟与体外实验，首次系统揭示了泛酸激酶2（PANK2）作为全氟辛酸（PFOA）和全氟辛烷磺酸（PFOS）诱导骨代谢紊乱的关键分子靶点。研究发现，PFOA和PFOS能够稳定结合PANK2蛋白并抑制其表达，进而破坏骨髓间充质干细胞（BMSCs）的成骨-成脂分化平衡，表现为成骨能力下降和成脂分化增强。机制上，PANK2可能通过调控辅酶A生物合成、能量代谢以及化学因子信号通路等途径影响骨骼稳态。该研究不仅为PFAS的骨毒性机制提供了novel的分子见解，也为预防和治疗环境污染物相关的骨骼疾病提供了潜在的治疗靶点。