《Toxicology and Applied Pharmacology》：Network toxicology, molecular docking, and multi-level bioinformatics integration reveal the targets and mechanisms linking benzo[
a]pyrene exposure to ischemic stroke
编辑推荐：
为了解决环境污染物苯并[a]芘（BaP）暴露与缺血性脑卒中（IS）风险之间的潜在机制不清问题，研究者整合网络毒理学、生物信息学与分子对接技术，开展主题研究。结果显示，MMP9和PTGS2为关键靶点，介导炎症反应、免疫浸润与脂质代谢，构建的列线图模型展现出稳健预测能力。该研究为理解环境毒物在IS中的作用提供了新视角，并提示了潜在的诊断生物标志物。
背景引入：隐藏的健康威胁与新探索
在我们生活的环境中，空气、土壤乃至食物中都可能潜藏着一种名为苯并[a]芘（Benzo[a]pyrene， BaP）的多环芳烃类环境污染物。人们主要通过呼吸和饮食途径接触到它。这种物质已被证实具有致癌、致突变和免疫毒性，而越来越多的流行病学证据提示，长期或短期暴露于BaP均与缺血性脑卒中（Ischemic Stroke， IS）风险的升高存在关联。然而，尽管关联性被观察到，环境污染物与脑卒中这类复杂神经系统疾病之间的“分子桥梁”却一直不甚清晰。传统风险因素如高血压、动脉粥样硬化已被广泛认知，但环境毒物如何从分子层面“攻击”血管与神经系统，进而触发或加剧脑卒中进程，是目前神经毒理学与脑血管病研究领域亟待解决的关键问题。针对这一空白，研究人员开展了一项创新性的整合研究，旨在系统揭示BaP诱发IS的潜在分子靶点与作用机制。这项研究最终发表于《Toxicology and Applied Pharmacology》杂志。
关键研究方法概览
本研究采用了一套整合多源数据与计算模拟的系统性方法。研究人员首先从STITCH、ChEMBL等数据库获取BaP的潜在靶点，同时从GeneCards、OMIM、TTD数据库收集IS相关靶点，并通过韦恩图取交集得到两者的共同靶点，构建了“BaP-靶点-IS”调控网络。从基因表达综合数据库（GEO）下载IS相关的转录组数据集（GSE58294、GSE22255），进行整合与批次效应校正后，通过差异表达分析筛选差异基因，并与共同靶点取交集得到关键基因。利用反转录定量聚合酶链式反应（RT-qPCR）在临床样本（来自宜兴市人民医院神经外科的12名IS患者和12名健康对照者外周血）中对关键基因表达进行验证。通过列线图构建风险预测模型，并利用独立数据集（GSE16561）进行外部验证。进一步利用单样本基因集富集分析（ssGSEA）评估免疫细胞浸润情况。最后，采用基于CB-Dock2平台的分子对接模拟，预测BaP与关键靶蛋白的空间结合模式与亲和力。
研究结果
1. 技术路线图
研究的技术路线图清晰地展示了从靶点获取、网络构建、数据分析到实验验证与模型预测的完整逻辑流程，体现了多技术整合的研究思路。
2. BaP与IS共同靶点的识别与网络构建
通过数据库整合分析，研究人员共获得113个BaP相关靶点和412个IS相关靶点，二者取交集得到21个共同靶点。围绕这21个靶点，成功构建了连接BaP暴露、分子靶点与IS的分子调控网络，初步勾勒出BaP可能影响IS的分子谱图。
3. 蛋白质相互作用网络与功能富集分析
对21个共同靶点进行蛋白质相互作用（PPI）网络分析，发现基质金属蛋白酶9（MMP9）和前列腺素内过氧化物合酶2（PTGS2）以及雌激素受体1（ESR1）是网络中连接度最高的核心节点。功能富集分析显示，这些共同靶点显著富集于炎症反应、脂质代谢相关通路（如脂质与动脉粥样硬化通路、IL-17信号通路）以及血管过程等生物学过程中。
4. 关键基因的筛选与验证
在整合的转录组数据中，共鉴定出193个差异表达基因（DEGs）。将其与上述21个共同靶点取交集，最终锁定MMP9和PTGS2为两个关键基因。分析显示，与对照组相比，这两个基因在IS样本中均显著上调。受试者工作特征曲线（ROC）分析表明，MMP9和PTGS2对IS具有一定的诊断区分能力，曲线下面积（AUC）分别为0.716和0.718。利用临床外周血样本进行的RT-qPCR实验进一步证实，IS患者中MMP9和PTGS2的表达水平显著高于健康对照。
5. 列线图风险预测模型的构建与验证
基于MMP9和PTGS2的表达水平，研究人员构建了一个用于预测IS风险的列线图模型。该模型在训练集（组合数据集）中表现出良好的预测性能（AUC=0.753），并在独立验证集GSE16561中取得了出色的外部验证效果（AUC=0.906）。校准曲线、决策曲线分析和临床影响曲线均表明该模型具有良好的校准度和潜在的临床应用价值。
6. 免疫浸润与关键基因的关联分析
免疫浸润分析发现，IS患者外周血中多种免疫细胞的浸润水平发生改变。相关性分析显示，MMP9的表达与中性粒细胞浸润呈显著正相关，而PTGS2的表达则与嗜酸性粒细胞浸润呈显著正相关。这表明关键基因的表达变化与特定的免疫微环境重塑密切相关。
7. 关键基因的基因集富集分析
单基因基因集富集分析揭示了MMP9高表达样本富集于白细胞介素-6（IL-6）信号通路、核因子-κB（NF-κB）信号通路和脂多糖诱导的炎症反应等通路。而PTGS2高表达样本则富集于转化生长因子-β1（TGF-β1）信号、肿瘤坏死因子（TNF）信号通路和炎症反应等。这进一步证实了这两个关键基因在驱动炎症反应中的核心作用。
8. 分子对接验证
分子对接模拟显示，BaP能够与MMP9和PTGS2的蛋白结构发生稳定的结合，其对接评分（结合自由能的近似值）分别为-9.9 kcal/mol和-10.9 kcal/mol。结合模式图显示，BaP主要通过疏水相互作用和π-π堆积（与MMP9）或传统氢键和疏水接触（与PTGS2）与靶蛋白结合，从结构层面支持了BaP可能直接作用于这两个靶点的可能性。
研究结论与重要意义
本研究的结论部分对以上发现进行了系统性总结与讨论，强调了其在理论、方法与实践层面的多重意义。
研究提出，环境毒物苯并[a]芘（BaP）可能通过调控MMP9和PTGS2这两个关键靶点，影响炎症反应、免疫细胞浸润（特别是中性粒细胞和嗜酸性粒细胞）以及脂质代谢等核心生物学过程，从而参与缺血性脑卒中（IS）的发病机制。MMP9作为一种能降解细胞外基质、破坏血脑屏障并促进神经炎症的酶，其上调可能介导了BaP引发的血管损伤。PTGS2（又称环氧合酶-2， COX-2）是炎症反应的关键介质，其激活会促进前列腺素等炎性介质的产生，加剧缺血后的炎症级联反应。二者可能通过协同作用，共同破坏了血管稳态并放大了炎症损伤。
这一发现的理论意义在于，它将环境暴露这一外部风险因素与脑卒中的内在分子病理网络直接连接起来，为“环境-基因”交互作用导致复杂疾病的假说提供了具体的分子证据和机制解释。它超越了传统的单一受体（如芳香烃受体， AhR）通路认知，提示了BaP可能通过多靶点、多途径的非基因组机制影响疾病易感性。
在方法论上，本研究成功示范了一套可复用的研究范式：整合毒理学与疾病数据库进行靶点预测与网络构建→利用公共组学数据筛选与验证关键分子→结合临床样本进行初步验证→通过计算模拟（分子对接）提供结构佐证→构建并验证临床应用模型（列线图）。这套“数据驱动-实验验证-模型预测”的整合策略，为研究其他环境神经毒物与神经系统疾病的关联提供了方法学参考。
研究的实践意义突出体现在两个方面：一是提供了潜在的生物标志物，MMP9和PTGS2在IS患者外周血中的上调及其诊断价值，提示它们可作为评估BaP暴露人群IS风险或辅助IS诊断的潜在指标。二是构建了一个具有良好预测性能的列线图模型，为临床医生提供了一种直观、量化的工具，用于对高危个体进行风险分层和早期预警。此外，研究结果也暗示，针对MMP9或COX-2（PTGS2编码的蛋白）的抑制剂，未来或可作为预防或治疗BaP相关IS的潜在干预策略。
当然，论文也客观指出了研究的局限性，例如使用的转录组数据来源于外周血而非脑组织，样本量相对有限，分子对接预测的结合模式需要后续功能实验验证，以及尚未建立BaP暴露剂量与关键基因表达或IS风险之间的定量关系等。这些都为未来的研究指明了方向，包括开展更大规模的前瞻性队列研究、利用细胞和动物模型深入探究BaP对靶点的功能调控、以及评估靶向干预策略的有效性等。
总而言之，这项研究通过创新的整合分析，系统揭示了环境污染物苯并[a]芘诱发缺血性脑卒中的潜在分子靶点与机制，将MMP9和PTGS2确立为连接环境暴露与脑血管事件的关键分子桥梁。它不仅增进了我们对环境因素在脑卒中病因学中作用的理解，也为该领域的风险预测、生物标志物发现和潜在治疗靶点探索提供了新的科学依据。
