《Environmental Science & Technology》：Integrated New Approach Methodologies Reveal the Potential Role of 2,7-Dibromocarbazole in Parkinson’s Disease via Monoamine Oxidase B Inhibition and Dopaminergic Dysfunction
编辑推荐：
作为编辑部，我们向读者推荐这篇研究。该文整合网络毒理学、分子对接、分子动力学模拟、转录组学及线虫在体实验等多种新方法体系（NAMs），系统阐述了环境污染物多卤代咔唑（PHCZs），特别是2,7-二溴咔唑（2,7-BCZ），通过抑制单胺氧化酶B（MAOB）、促进α-突触核蛋白聚集、扰乱多巴胺能突触功能等通路，在帕金森病（PD）发病机制中的潜在作用，为评估新兴污染物的神经毒性风险提供了创新性的多组学整合研究范例。
引言：环境中的潜在神经毒物
多卤代咔唑（Polyhalogenated carbazoles, PHCZs）是一类具有二噁英样效应的持久性有机污染物，近年来因其广泛的环境分布和潜在的毒性而受到关注。它们已在农田、淡水系统、室内灰尘、PM2.5颗粒物乃至北极土壤和人体尿液等多种环境介质中被检出。在中国农田土壤中，3,6-二氯咔唑（3,6-CCZ）是主要同系物，而在北京大气PM2.5中，3-氯咔唑（3-CCZ）、3-溴咔唑（3-BCZ）和3,6-CCZ占据了PHCZs总量的93%。特别值得注意的是，电子垃圾处理场土壤中的PHCZs浓度远高于其他污染物。尽管已有研究揭示了PHCZs的芳烃受体（AhR）激活、斑马鱼发育毒性以及对土壤微生物过程的影响，但它们与神经退行性疾病，特别是帕金森病（Parkinson‘s disease, PD）之间的具体联系和机制仍不清楚。本研究旨在利用整合新方法体系（New Approach Methodologies, NAMs），系统评估PHCZs的神经毒性及其潜在机制。
PHCZs的理化性质与环境毒性预测
研究首先评估了13种常见PHCZs的潜在环境毒性。通过ADMETlab平台进行的毒性预测显示，其脂溶性（以LogP表示）并未与分子量呈现一致的正相关。总体而言，PHCZs的毒性随分子量增加而增强，但3-BCZ和3,6-碘咔唑（3,6-ICZ）的生物浓缩因子（Bioconcentration Factor, BCF）相对较低。
识别PHCZs在PD发病中的潜在靶点
为了识别PHCZs在PD中的潜在分子靶点，研究通过GeneCards和OMIM数据库获取了PD相关基因，并利用SuperPRED、NetInfer和SEA等在线工具预测了PHCZs的潜在靶点基因。两者取交集后，获得了68个候选靶基因。基因富集分析表明，这些重叠靶点主要与膜电位调节、轴突远端功能、突触后神经递质受体活性、神经活性配体-受体相互作用、多巴胺能突触等通路相关。进一步构建这68个靶点的蛋白质-蛋白质相互作用（Protein-Protein Interaction, PPI）网络，识别出八个枢纽靶点：STAT3、APP、GSK3B、DRD2、SIRT1、MAOB、HTT和MAOA。
通过分子对接验证PHCZs与核心靶点的直接相互作用
为探究PHCZs与这八个枢纽靶点的结合亲和力，研究对所有PHCZs与这些蛋白质进行了分子对接。结果显示，单胺氧化酶B（Monoamine oxidase B, MAOB）对PHCZs的平均结合亲和力最强（平均Vina评分 = -7.8 kcal/mol）。虽然在所有PHCZs中，1-溴-3,6-二氯咔唑（1-B-3,6-CCZ）显示出最高的结合亲和力，但考虑到二卤代PHCZs在环境中更为常见和丰富，研究选择了2,7-二溴咔唑（2,7-BCZ）进行后续深入分析。分子对接表明，2,7-BCZ对MAOB具有最高的结合亲和力（-9.4 kcal/mol）。
分子动力学模拟与体外验证
比较分析发现，2,7-BCZ与已知的MAOB抑制剂A1IDI在酶的活性位点占据了几乎相同的结合构象。进一步的分子动力学（Molecular Dynamics, MD）模拟（持续50 ns）表明，MAOB–2,7-BCZ复合物在回转半径（Rg）、均方根波动（RMSF）和溶剂可及表面积（SASA）方面均保持结构稳定。基于MD轨迹的MM/PBSA结合自由能计算得出总结合自由能为-35.84 kJ/mol，证实了相互作用的稳定性和有利性。可视化分析显示，2,7-BCZ通过疏水和极性相互作用占据催化口袋。分子动力学模拟还显示出MAOB与2,7-BCZ之间持续存在氢键作用。定量MAOB检测进一步表明，暴露于2,7-BCZ会诱导MAOB蛋白表达增加，这可能是细胞对2,7-BCZ抑制MAOB酶活性的一种代偿性反应。
在体验证：2,7-BCZ高度促进α-突触核蛋白聚集
为了评估2,7-BCZ是否会导致PD样病理变化，研究使用了表达人源α-突触核蛋白与黄色荧光蛋白融合体（Punc-54::α-synuclein::YFP）的转基因秀丽隐杆线虫（C. elegans）品系NL5901。暴露于1 μM 2,7-BCZ显著促进了线虫体壁肌肉细胞中α-突触核蛋白的聚集。虽然观察到的聚集发生在肌肉细胞而非神经元中，但此反应反映了与PD更广泛发病机制相关的早期蛋白质毒性和细胞应激反应。
2,7-BCZ处理的SH-SY5Y细胞的转录组学特征
为探究其分子机制，研究对暴露于100 μM 2,7-BCZ的人神经母细胞瘤SH-SY5Y细胞进行了转录组测序。主成分分析（Principal Component Analysis, PCA）揭示了全基因组水平的显著转录组改变。差异表达分析共识别出9126个显著失调的基因。基因富集分析显示，这些差异表达基因主要涉及染色体区域、钙黏蛋白结合、肌动蛋白细胞骨架调控以及多种疾病的神经退行性变等通路。
PD与2,7-BCZ暴露的共享转录特征
为了探索2,7-BCZ与PD之间的关联，研究分析了来自GEO数据库的五个PD患者与健康对照者黑质微阵列数据集，共计90个样本。在去除批次效应后，差异表达分析识别出17个显著失调的基因。这些基因主要富集于多巴胺生物合成过程、神经元投射末端、多巴胺结合、PD、酪氨酸代谢和突触小泡循环等通路。值得注意的是，CLSTN2、CBLN1、AGTR1、DLK1和DDC这五个基因在PD患者数据集和2,7-BCZ处理的SH-SY5Y细胞数据集中均表现出一致的下调，提示了暴露与疾病之间共享的转录扰动。DDC编码芳香族氨基酸脱羧酶，该酶负责将L-DOPA转化为多巴胺，其表达降低意味着多巴胺合成能力下降。
2,7-BCZ在PD中毒性机制的综合分析
为阐明2,7-BCZ可能导致PD的途径，研究交叉分析了三个通路集合：S1（重叠靶基因的富集通路）、S2（PD转录组中的差异表达基因富集通路）和S3（2,7-BCZ暴露细胞转录组中的差异表达基因富集通路）。交叉得到的通路包括多巴胺能突触、多细胞生物生长、前脑发育、神经系统发育调节、氨基酸代谢过程、突触组织调节、突触结构或活性调节、细胞大小调节、小脑发育和中脑发育。这些发现强调了2,7-BCZ可能在PD背景下损害的关键神经发育和神经传递相关过程。
讨论与总结
本研究在NAMs框架内，综合运用毒性预测、网络毒理学、分子亲和力评估、在体验证和转录组学分析，系统评估了PHCZs与PD相关的神经毒性。环境分析与毒性预测表明，卤素类型和分子量影响PHCZs的脂溶性和生物累积性。网络毒理学识别出PHCZs与PD共享的八个枢纽靶点，其中MAOB与PHCZs的结合亲和力最高。鉴于其环境普遍性及对MAOB的高亲和力，2,7-BCZ被选作深入研究的对象。比较分析与分子动力学模拟支持2,7-BCZ作为MAOB抑制剂的角色。线虫在体实验证实2,7-BCZ显著加剧了PD病理标志α-突触核蛋白的聚集。SH-SY5Y细胞的转录组分析揭示了其对多巴胺能突触功能、突触组织和神经发育通路的广泛影响。与PD患者数据的整合分析发现了共享的下调基因集，特别是与多巴胺合成相关的DDC。综合通路分析进一步汇聚于多巴胺能突触调节、前脑与小脑发育以及突触结构组织等关键过程。这些结果表明，2,7-BCZ可能通过结合并抑制MAOB、破坏多巴胺代谢稳态、促进α-突触核蛋白聚集，并导致关键神经通路的转录失调，从而在PD发病机制中扮演潜在的危险角色。本研究采用节约动物的NAMs，强调了监测环境中新兴PHCZs污染物的必要性，并凸显了多组学整合方法在评估其神经毒性风险方面的价值。
