《Ecotoxicology and Environmental Safety》：Computational screening reveals class-wide CYP19A1 binding potential across per- and polyfluoroalkyl substances
编辑推荐：
本研究针对PFAS（全氟及多氟烷基物质）与关键雌激素合成酶CYP19A1（芳香酶）相互作用机制不清的问题，开发了基于药效团的虚拟筛选框架，对10,946种EPA名录PFAS进行了系统性评估。研究发现90.1%的PFAS具有CYP19A1结合潜力，其中7078种为高置信度结合物。分子动力学模拟表明PFHxS比PFOA结合更稳定，CETSA实验证实了PFHxS的直接结合。该研究挑战了现有逐项监管模式，为基于类别的PFAS生殖毒理学风险评估提供了机制支持。
在现代化工产品无处不在的今天，有一类被称为“永久化学品”的物质——全氟及多氟烷基物质（PFAS），正悄然成为全球性的环境健康挑战。从防水服装到不粘锅涂层，从消防泡沫到食品包装，PFAS的应用极其广泛，但它们的持久性和生物累积性也带来了严重的环境与健康风险。尤其令人担忧的是，越来越多的证据表明PFAS暴露与生殖健康问题密切相关，包括生育能力下降、性激素水平紊乱等。然而，面对超过12,000种已知的PFAS化合物，传统的毒理学测试方法显得力不从心，我们迫切需要更高效的方法来评估这类化学物质的潜在健康风险。
在这一背景下，CYP19A1（芳香酶）作为体内雄激素向雌激素转化的关键限速酶，成为了理解PFAS生殖毒性的重要突破口。雌激素平衡对维持生殖健康、胎儿发育和性分化至关重要，任何干扰都可能带来严重后果。尽管已有研究表明某些PFAS（如PFOS和PFOA）能够影响CYP19A1的功能，但大多数PFAS化合物是否以及如何与这一关键酶直接相互作用，仍然是一个未解之谜。
为了解决这一知识空白，发表在《Ecotoxicology and Environmental Safety》上的最新研究开发了一套创新的计算毒理学框架，首次系统评估了PFAS作为一类化学物质与CYP19A1的结合潜力。研究人员从美国环境保护署（EPA）化学品目录中收集了10,946种PFAS化合物，利用基于结构的药效团建模、分子对接、分子动力学模拟和细胞实验验证相结合的多层次策略，深入探索了PFAS-CYP19A1相互作用的分子基础。
研究主要采用了几个关键技术方法：首先，基于PFOS与CYP19A1的对接复合物生成并验证了10个药效团模型，从中筛选出3个高性能模型（AUC > 0.848，MCC 0.257-0.294）用于大规模虚拟筛选。其次，对代表性PFAS化合物进行了分子对接和100纳秒的分子动力学模拟，评估结合稳定性和相互作用模式。最后，通过细胞热位移分析（CETSA）实验验证了PFHxS与CYP19A1的直接结合。所有统计分析使用R 4.4.1和GraphPad Prism 10完成。
3.1. 分子对接揭示PFOS-CYP19A1结合相互作用
分子对接分析显示，PFOS与CYP19A1活性位点形成了六种有利相互作用，包括四个氢键和两个疏水接触，以及五个卤素键。关键相互作用残基包括ARG115、ILE133、ALA306等，表明PFAS化合物能够通过多种分子相互作用与CYP19A1结合腔稳定结合。
3.2. 从PFOS-CYP19A1相互作用生成药效团模型
从PFOS-CYP19A1结合复合物中生成了10个药效团模型，这些模型在特征数量和选择性值上各不相同，包含氢键受体（A）、疏水特征（F）和卤素键合位点（H）等不同组合，为后续虚拟筛选提供了多样化的分子识别模式。
3.3. 药效团模型验证识别三个高性能模型
通过包含12个阳性对照PFAS和600个诱饵化合物的验证数据集评估，Pharm02、Pharm03和Pharm04三个模型表现出最优性能（AUC > 0.8，MCC > 0.25），其中Pharm03的AUC达到0.936，MCC为0.294，表现出完美的灵敏度（1.000）和0.873的特异性。
3.4. 大规模虚拟筛选揭示PFAS化学品中广泛的CYP19A1结合潜力
大规模虚拟筛选结果显示，10,946种PFAS中有9861种化合物（90.1%）表现出可测量的CYP19A1结合潜力。三个药效团模型的交叉分析发现7078种化合物匹配所有三个模型，被视为高置信度结合物，表明CYP19A1结合是PFAS的一类广泛特性。
3.5. 代表性PFAS与CYP19A1的分子对接
对PFOA和PFHxS的分子对接比较显示，PFHxS形成了更广泛的相互作用网络（七个有利相互作用），而PFOA仅显示五个氢键和四个卤素键，没有疏水相互作用，提示磺酸基团可能比羧酸基团形成更强的结合相互作用。
3.6. 通过分子动力学模拟进行计算稳定性评估
100纳秒的分子动力学模拟表明，PFHxS-CYP19A1复合物比PFOA-CYP19A1复合物更稳定，PFHxS的配体RMSD维持在0.2-0.3纳米，而PFOA稳定在0.8纳米左右。吉布斯自由能景观分析进一步证实了PFHxS结合复合物的热力学稳定性更优。
3.7. 使用CETSA进行实验验证确认PFHxS-CYP19A1体外结合
细胞热位移分析（CETSA）实验证实，PFHxS处理显著提高了CYP19A1的热稳定性，在55-65°C的温度范围内，PFHxS处理组的可溶性CYP19A1比例显著高于DMSO对照组（p = 0.024），直接证明了PFHxS与CYP19A1的物理结合。
这项研究的讨论部分强调了几个重要发现。首先，90.1%的PFAS具有CYP19A1结合潜力的发现挑战了当前对PFAS内分泌干扰机制的理解，表明直接酶结合可能是之前未被充分认识的毒性通路。其次，分子动力学模拟揭示的磺酸基PFAS（如PFHxS）比羧酸基PFAS（如PFOA）结合更稳定的现象，为理解不同PFAS化合物的相对毒性提供了分子基础。此外，研究建立的计算-实验验证框架为高效评估其他PFAS的健康风险提供了可行方案。
该研究的结论部分明确指出，CYP19A1结合是PFAS的一类广泛特性，而非仅限于少数传统化合物。这一发现对PFAS的风险评估和监管政策具有深远意义，支持了将PFAS作为一类化学物质进行分组监管的科学基础，而非当前的逐项评估模式。鉴于CYP19A1在雌激素生物合成和生殖稳态中的核心作用，直接酶结合代表了PFAS生殖毒性的合理分子机制，为未来开发更安全的化学替代品和减轻PFAS引起的生殖健康影响提供了重要线索。
总的来说，这项研究不仅填补了PFAS内分泌干扰机制的知识空白，也为应对这一全球性环境健康挑战提供了新的科学工具和视角。随着PFAS污染的持续蔓延，此类研究对于保护公众生殖健康、指导化学品安全设计以及制定科学合理的环境政策都具有不可估量的价值。
