全氟烷基酸与PPARδ相互作用分子机制的计算生物学研究：结合模式、关键残基与结构-活性关系分析

《Toxics》：Study on the Molecular Mechanism of Interaction Between Perfluoroalkyl Acids and PPAR by Molecular Docking

本研究通过分子对接和动力学模拟，系统揭示了20种全氟烷基酸（PFAAs）与过氧化物酶体增殖物激活受体δ（PPARδ）的相互作用机制。研究发现PFAAs结合强度受分子量、氢键供体数和熔点显著影响，结合位点集中于TRP-256等关键残基。支链结构会因空间位阻减弱结合能力，分子动力学验证了结合稳定性，为PFASs的生态风险评估提供了分子层面依据。

　　

研究背景

全氟和多氟烷基物质（PFASs）作为具有高环境持久性和生物累积性的"永久化学品"，其毒性机制研究备受关注。全氟烷基酸（PFAAs）是PFASs的重要组成部分，具有肝毒性、内分泌干扰等多种毒效应。过氧化物酶体增殖物激活受体（PPAR）是核受体超家族重要成员，在调节能量代谢和脂质稳态中起核心作用。由于PFAAs与内源性脂肪酸结构相似，PPAR可能是其发挥脂毒性的重要靶点。

研究方法

本研究采用分子对接技术，以PDB ID 3OZ0的PPARδ晶体结构为受体，系统分析了20种结构多样的PFAAs（包括直链/支链羧酸和磺酸）的结合模式。通过AutoDock 4.2进行半经验自由能力场计算，结合拉马克遗传算法进行构象搜索。选取5种新型PFAAs作为验证组，采用GROMACS 2024.5进行100 ns分子动力学（MD）模拟，验证结合稳定性。

分子相互作用特征

所有PFAAs均能进入PPARδ结合口袋，主要作用力为氢键和卤键。关键结合残基集中在TRP-256、ASN-269、GLY-270等高频区域。全氟烷基磺酸（PFSAs）结合亲和力（平均结合能-3.97 kcal·mol^-1）显著优于全氟烷基羧酸（PFCAs）。典型配体PFBA与PPARδ结合能为-2.99 kcal·mol^-1，主要与TRP-256等残基相互作用；而支链代表3m-PFOA结合能减弱至-2.92 kcal·mol^-1。

结构-活性关系

相关性分析表明，结合能绝对值与分子量、氢键供体数、熔点呈显著正相关，与拓扑极性表面积负相关。直链PFCAs结合能随分子量增加而降低（R²=0.7715），说明空间位阻效应削弱结合能力。支链异构体（如3,3m₂-PFOA结合能-2.45 kcal·mol^-1）因空间位阻效应表现出比直链对应物更弱的结合能力。

验证组实验

5种验证组PFAAs的结合模式与实验组高度一致，高频残基重复率达83.33%。PFPA、PFOA等配体的结合强度变化规律符合实验组结论，进一步证实分子量增加和支链结构会降低结合亲和力。

分子动力学验证

PFBA、PFOA等体系的RMSD值在20 ns后稳定在0.28-0.45 nm范围，证实对接构象在溶剂化和蛋白质柔性条件下的动态稳定性。比较显示分子量增大导致结合能降低，支链结构（3,3m₂-PFOA）的RMSD波动范围大于直链对应物，从动力学角度验证了空间位阻效应。

研究意义

本研究首次系统揭示PFAAs通过PPARδ通路干扰代谢稳态的分子机制，为PFASs的混合暴露风险评估提供重要依据。发现支链修饰是降低PPARδ结合亲和力的有效分子修饰策略，为开发更安全的PFASs替代品提供计算设计思路。未来结合自由能计算将为此相互作用提供更定量的热力学证据。