整合网络毒理学与分子对接揭示PFOA暴露致男性生殖功能障碍的共享核心靶点及激素–NO–凋亡/炎症机制

《Basic and Clinical Andrology》：Integrative network toxicology analysis reveals mechanistic links between PFOA exposure and male reproductive dysfunction

目的：在全氟辛酸（perfluorooctanoic acid，PFOA）网络毒理学框架内，系统探究PFOA与男性生殖健康结局之间的分子机制，并比较勃起功能障碍（erectile dysfunction，ED）、男性不育（male infertility）和睾丸功能障碍（testicular dysfunction）之间的共享通路。方法：使用比较毒理基因组学数据库（Comparative Toxicogenomics Database，CTD）、SwissTargetPrediction和TargetNet识别PFOA相关靶点；从GeneCards收集疾病相关基因。使用STRING构建蛋白–蛋白相互作用（protein–protein interaction，PPI）网络，并在Cytoscape中分析以识别枢纽基因。进行基因本体（Gene Ontology，GO）和京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes，KEGG）富集分析，以表征受PFOA影响的生物学过程和通路。采用分子对接（molecular docking）评估PFOA与关键靶蛋白之间的结合亲和力。此外，研究人员进行了初步体外实验，以检测PFOA暴露后所选候选基因的表达模式。结果：共识别105个PFOA相关基因。交集分析分别得到15、25和42个与ED、男性不育和睾丸功能障碍相关的基因。10个核心基因（AR、ESR1、BCL2、NOS1、NOS2、NOS3、ACE、SIRT1、TLR9、SHBG）在所有表型中共享，并在PPI网络中显示高度连接性。富集分析提示其参与生殖系统发育、激素和核受体调控、钙和一氧化氮（nitric oxide，NO）信号、线粒体和氧化还原过程，以及凋亡和炎症相关通路。初步体外实验进一步验证了所选关键基因的表达变化。结论：该研究识别了可能介导PFOA相关男性生殖毒性的关键分子靶点和通路。结果支持一个工作假说，即PFOA相关男性生殖毒性可能涉及雄激素–雌激素受体信号改变、NO相关通路以及凋亡/炎症相关调控。这些发现为未来机制验证和风险评估提供了理论基础。

研究背景与问题：男性生殖健康已成为全球公共卫生问题，男性不育、勃起功能障碍（erectile dysfunction，ED）和睾丸功能障碍（testicular dysfunction）患病率近年上升，给个体、家庭和社会带来负担。环境因素中，化学污染物因普遍存在、持久、抗降解和生物累积而受到关注，尤其是内分泌干扰化学物（endocrine-disrupting chemicals，EDCs），其可模拟或阻断激素信号，干扰激素合成与代谢，影响精子发生、性功能以及睾丸结构和功能。全氟辛酸（perfluorooctanoic acid，PFOA）是典型全氟和多氟烷基物质（per- and polyfluoroalkyl substances，PFAS），化学稳定性强、环境持久并可在生物体累积，被称为“永久化学物”。动物和流行病学研究提示，PFOA暴露与血清睾酮下降、精子参数异常、睾丸和附睾绝对/相对重量改变及睾丸组织病理变化相关。然而，既往研究多描述PFOA暴露与生殖结局的关联，分子机制仍不充分；且多聚焦单一生殖终点，较少在统一框架内评估多个男性生殖表型。因此，研究人员采用网络毒理学（network toxicology）整合毒物、靶点、通路与疾病，探讨PFOA与ED、男性不育和睾丸功能障碍的共享机制，为风险评估、机制研究和公共卫生干预提供理论依据。

研究人员开展的研究、结论与意义：研究人员整合CTD、SwissTargetPrediction和TargetNet获取PFOA相关靶点，从GeneCards收集疾病相关基因，构建蛋白–蛋白相互作用（protein–protein interaction，PPI）网络，进行基因本体（Gene Ontology，GO）和京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes，KEGG）富集分析，并采用分子对接（molecular docking）评估结合亲和力，同时在Leydig细胞和大鼠海绵体内皮细胞（rat corpus cavernosum endothelial cells，CCECs）中进行初步体外验证。研究共识别105个PFOA相关基因，交集分析分别得到15、25和42个与ED、男性不育和睾丸功能障碍相关基因，并确定10个共享核心基因：AR、ESR1、BCL2、NOS1、NOS2、NOS3、ACE、SIRT1、TLR9、SHBG。富集分析提示生殖系统发育、激素和核受体调控、钙和一氧化氮（nitric oxide，NO）信号、线粒体和氧化还原过程、凋亡和炎症相关通路参与其中。分子对接显示除NOS1外，其余核心蛋白与PFOA具有较强结合倾向。体外实验验证了部分关键基因mRNA表达变化。研究结论支持工作假说：PFOA相关男性生殖毒性可能涉及雄激素–雌激素受体信号改变、NO相关通路以及凋亡/炎症相关调控，为未来机制验证和风险评估提供了候选靶点和理论框架。

关键技术方法（不超过250字）：研究采用网络毒理学、分子对接和体外验证。PFOA靶点经CTD、SwissTargetPrediction和TargetNet获取，疾病基因来自GeneCards；STRING构建PPI网络，Cytoscape进行拓扑分析和枢纽基因筛选；GO/KEGG富集分析；CB-Dock2进行分子对接；体外使用Leydig细胞和CCECs（均Procell，中国），设正常对照及5、10、20 μmol/L PFOA处理，RT-qPCR检测mRNA表达，相对表达采用2^?ΔΔCt法，统计使用GraphPad Prism 10和R 4.2.0。

研究结果：
PFOA的毒性预测：通过PubChem获取PFOA化学结构，并用admetSAR 3.0预测其生殖毒性，结果提示生殖毒性概率为30.1%，表明PFOA可能具有潜在生殖毒性。
PFOA对ED的潜在影响：获得105个PFOA相关基因和1,687个ED相关基因，交集识别15个PFOA–ED相关基因。PPI网络提示BCL2、激素受体相关节点和一氧化氮合酶（nitric oxide synthase，NOS）家族可能发挥重要作用。GO富集涉及性腺发育、性别分化、生殖结构形成、生殖系统发育、类固醇激素反应、钙离子转运调控和血压调节等；KEGG富集涉及精氨酸生物合成、精氨酸/脯氨酸代谢、缺氧诱导因子1（hypoxia-inducible factor 1，HIF-1）信号、松弛素信号、雌激素信号、apelin信号、流体剪切应力和动脉粥样硬化、晚期糖基化终产物–其受体（advanced glycation end-product–receptor for advanced glycation end-product，AGE–RAGE）信号及心脏肥厚相关通路，提示PFOA可能通过激素信号、钙稳态、氧化还原、NO合成、内皮依赖性血管舒张和血流动力学调控影响ED。
PFOA对男性不育的潜在影响：获得105个PFOA相关基因和2,315个男性不育相关基因，交集识别25个相关基因。PPI网络提示凋亡相关分子、激素受体相关节点和NOS家族为核心成分。GO富集涉及细菌来源分子和脂多糖（lipopolysaccharide，LPS）反应、炎症或免疫失调、性腺发育、第一性征发育、性别分化、节律过程、钙离子转运负调控、核受体活性、氧化还原酶活性、类固醇结合和血红素结合等；KEGG富集涉及精氨酸生物合成、精氨酸/脯氨酸代谢、氮代谢、肾素–血管紧张素系统、apelin信号、前列腺癌、小细胞肺癌、化学致癌、结核和恰加斯病等，提示PFOA可能通过NO合成与代谢稳态、睾丸微环境、炎症、免疫失调、DNA损伤反应和凋亡失衡影响男性不育。
PFOA对睾丸功能障碍的潜在影响：获得105个PFOA相关基因和3,064个睾丸功能障碍相关基因，交集识别42个相关基因。PPI分析提示BCL2、ESR1和SIRT1为关键网络节点。GO富集涉及类固醇激素反应、激素介导信号、细胞内受体信号、LPS和细菌来源分子免疫反应、酮体反应、小分子代谢调控、膜筏、核膜和外核膜、细胞器外膜、线粒体外膜、小窝富集膜域和真染色质区域等；分子功能涉及核受体活性、配体激活转录因子活性、DNA结合转录因子结合、RNA聚合酶II特异性转录因子结合、氧化还原酶活性和转录共激活因子结合。KEGG富集涉及精氨酸/脯氨酸代谢、精氨酸生物合成、apelin信号、凋亡、胰岛素抵抗、癌症相关通路和化学致癌–受体激活通路，提示PFOA可能通过代谢紊乱、凋亡激活和信号通路扰动损害睾丸组织。
PFOA对男性生殖健康的潜在影响：PFOA相关基因与所有表型相关基因集交集识别10个PFOA–男性生殖健康相关基因。PPI结果提示BCL2、AR、NOS家族（NOS1、NOS2、NOS3）和ESR1为中心网络节点。GO富集涉及性腺发育、男性第一性征发育、男性性别分化、生殖系统发育和钙离子转运负调控；细胞组分显著富集于线粒体；分子功能涉及钙离子结合、核受体活性、氧化还原酶活性、类固醇结合、雌激素受体结合和烟酰胺腺嘌呤二核苷酸磷酸（nicotinamide adenine dinucleotide phosphate，NADP）结合。KEGG富集涉及钙信号、apelin信号、雌激素信号、HIF-1信号、松弛素信号、精氨酸和脯氨酸代谢、精氨酸生物合成及化学致癌–受体激活通路。基因–通路关联分析中，NOS1、NOS2、NOS3、ACE、AR、ESR1和BCL2连接多个通路，提示NOS家族与钙信号、HIF-1信号和氨基酸代谢相关，AR和ESR1强调雄激素与雌激素受体信号，BCL2提示凋亡参与。
分子对接：对10个核心蛋白进行分子对接，结果显示除NOS1外，其余蛋白表现出较强结合倾向（Vina score < ?7 kcal/mol）。SHBG（?9.3 kcal/mol）、AR（?9.2 kcal/mol）、SIRT1（?9.1 kcal/mol）和ESR1（?8.9 kcal/mol）结合能最有利；NOS2（?8.4 kcal/mol）、NOS3（?8.3 kcal/mol）、NOS1（?4.8 kcal/mol）表现出不同程度结合；BCL2（?7.8 kcal/mol）、ACE（?8.1 kcal/mol）和TLR9（?8.0 kcal/mol）也显示有利结合。结果提示PFOA可能与性激素结合球蛋白、雄激素受体、雌激素受体、SIRT1、NOS、BCL2、ACE和TLR9等关键蛋白稳定结合，进而影响内分泌、NO、凋亡、肾素–血管紧张素系统和先天免疫。
功能分析：研究人员基于网络毒理学、分子对接和功能特征构建机制示意图，将10个共享关键基因分为五个功能模块：内分泌干扰、NO信号和血管调节、线粒体功能和代谢调节、凋亡调节、炎症/先天免疫反应。AR、ESR1和SHBG主要参与雄激素–雌激素平衡和激素生物利用度调控；NOS1、NOS2、NOS3和ACE主要富集于NO产生和血管稳态；SIRT1主要与线粒体稳态、氧化还原平衡和代谢调控相关；BCL2主要参与凋亡调控；TLR9与炎症反应和先天免疫激活密切相关。
RT-qPCR：在Leydig细胞和CCECs中检测AR、ESR1、NOS3、SIRT1、BCL2和TLR9 mRNA表达，PFOA浓度设为0、5、10和20 μmol/L。在Leydig细胞中，PFOA处理使AR、ESR1和TLR9 mRNA表达升高，NOS3和SIRT1无显著变化，BCL2在高浓度后呈下降趋势。在CCECs中，PFOA处理上调AR、ESR1和TLR9，下调NOS3和BCL2，SIRT1无显著变化。结果提示PFOA对关键基因表达的影响具有细胞类型特异性。

讨论部分总结：研究人员认为，PFOA相关男性生殖损伤不太可能依赖单一分子事件，而更可能由多条通路协同介导，包括内分泌失衡、NO信号异常以及凋亡–炎症失衡。第一，内分泌干扰可能是中心启动事件。AR和ESR1在三种表型中均处于共享核心位置，分子对接显示强结合倾向；SHBG结合能最有利之一，提示PFOA可能影响性激素的循环生物利用度和受体暴露水平。富集分析显示核心基因与核受体活性、配体激活转录因子活性、类固醇结合、激素介导信号和雌激素信号相关；体外实验中AR和ESR1在Leydig细胞和CCECs中mRNA表达升高，支持PFOA可能诱导内分泌相关异常激活，但这种上调不一定代表生理增强，可能是应激相关重塑或受体网络失衡激活。第二，NO信号异常和血管/微循环稳态受损可能是连接ED和睾丸功能障碍的重要桥梁。NOS家族在共享核心网络中突出，并与精氨酸生物合成、精氨酸/脯氨酸代谢、钙信号、HIF-1信号、松弛素信号和apelin信号相关。分子对接显示NOS2和NOS3具有较稳定结合倾向；体外CCECs中NOS3下降，而Leydig细胞中无明显变化，提示PFOA对NO相关损伤可能具有细胞类型特异性。第三，凋亡和炎症失调可能是共同终末效应通路。BCL2作为经典抗凋亡分子，在网络中高度连接并在对接中显示稳定结合；TLR9作为先天免疫识别分子，处于共享网络关键位置，其异常激活可能促进慢性炎症或免疫失衡。第四，SIRT1介导的线粒体和代谢失调可能提供背景机制。尽管体外实验中两种细胞SIRT1表达无显著变化，但研究人员认为这可能与功能活性、亚细胞定位或下游调控有关，当前处理时长和浓度可能不足以捕捉其动态变化，SIRT1可能更多作为背景调节因子。研究人员提出多通路协同模型：PFOA可能首先作用于AR/ESR1/SHBG相关内分泌轴，扰乱雌激素和雄激素信号及受体介导的转录调控；随后通过NOS家族、钙信号和精氨酸代谢影响NO产生及血管/微循环功能；在代谢应激和线粒体功能障碍背景下，进一步诱导炎症反应和凋亡调控紊乱，最终导致ED、男性不育和睾丸功能损害等多种表型。局限性包括：PFOA相关靶点主要来自数据库检索和算法预测，可能受数据覆盖、证据水平和模型偏倚影响，且缺乏剂量、暴露时长和组织特异性信息；疾病相关基因来自GeneCards并以相关性评分阈值过滤，可能排除低评分但潜在重要基因；分子对接仅提供结构层面结合倾向推断，需分子动力学模拟和进一步实验验证；未来需动物实验和大规模流行病学研究验证结论。

研究结论翻译：该研究在系统生物学层面识别了PFOA相关男性生殖健康损害的潜在关键靶点和通路，强调内分泌干扰信号在多个生殖终点中的核心作用。结果进一步提示，一氧化氮（NO）轴失调以及炎症/凋亡通路激活可能与激素紊乱协同，促进男性生殖损伤进展。尽管仍需进一步实验验证，这些发现提名了特定候选靶点，并为男性生殖健康的预防策略和公共卫生风险评估提供了框架。