PPARγ共价拮抗剂T0070907与GW9662在Cys285共价修饰前的非共价预反应稳定化机制——整合分子对接、50 ns分子动力学与MM-PBSA的比较研究

《Results in Chemistry》：Non-covalent pre-reactive stabilization of covalent antagonists prior to reactive cysteine modification in PPARγ

　　

过氧化物酶体增殖物激活受体γ（PPARγ）是参与代谢调控的关键核受体，也是代谢紊乱的重要治疗靶点。T0070907和GW9662是可共价且不可逆地修饰Cys285的PPARγ拮抗剂，但其非共价预反应稳定化仍未充分阐明。在此，研究人员整合分子对接、50 ns分子动力学（MD）模拟与MM-PBSA计算，以表征二者在PPARγ配体结合结构域（LBD）内的结合。对接显示两种配体均具有有利的预反应结合，GW9662得分略更有利。MD模拟揭示不同的动态稳定化模式：T0070907表现出较低的主链涨落和更频繁的瞬态氢键形成，而GW9662表现出更大的构象灵活性和以疏水作用为主的稳定化。MM-PBSA分析给出GW9662较T0070907更有利的预测结合自由能，这与其更强的范德华和静电相互作用相关，并补偿了其更高的极性溶剂化罚分。对Cys285–弹头距离的时间分辨分析进一步揭示通向预反应定位的不同时间路径。GW9662发生明显的晚期重排，随后在采样轨迹末端附近实现其Cys285邻近几何的局部稳定化；T0070907则先发生初始结构弛豫，之后维持相对持续的邻近。总体而言，结果表明两种拮抗剂在共价修饰前采用互补的稳定化策略：T0070907更强依赖方向性极性相互作用并伴随更大的动态稳定性，而GW9662表现出更有利的预测结合能量学和更紧密的预反应定位。这些发现为理解和优化共价PPARγ拮抗剂的预反应识别提供了可检验的计算假设。

研究背景方面，PPARγ是配体激活的核受体，参与脂肪生成、脂质代谢、葡萄糖稳态和炎症信号调控，其活性失调与2型糖尿病、肥胖和心血管疾病相关，因而成为代谢和炎症疾病的重要治疗靶点。早期噻唑烷二酮类全激动剂虽增强胰岛素敏感性，但受体重增加、液体潴留和心血管风险等不良反应限制；T0070907和GW9662作为共价拮抗剂，可不可逆修饰PPARγ LBD内保守Cys285。既有研究多集中于Cys285共价加合物及共价修饰后的功能后果，而对共价键形成前的早期识别、配体取向、稳定化和亲电体定位机制理解不足。成功共价抑制不仅取决于最终成键反应，也取决于配体能否通过非共价相互作用达到并维持有利的预反应取向，因此需要比较结构相关拮抗剂的预反应非共价稳定化网络，并区分动态姿态持续性与预测热力学结合有利性。研究人员采用整合计算工作流开展比较研究，得出两种配体共享关键锚定相互作用但采用互补稳定策略的结论，并提出可检验计算假设；其意义在于为理解共价PPARγ拮抗剂的预反应识别、选择性和优化提供结构-动态-能量框架。

关键技术方法概括（不超过250字）：该研究为计算研究，无实验样本队列；蛋白结构来自PDB ID 6ONI B链，配体三维结构来自PubChem SDF。研究人员先用CB-Dock2进行盲对接和结合腔评估，再用UCSF Chimera中的AutoDock Vina进行聚焦对接；用PLIP、ProteinsPlus、Schr?dinger Maestro和BIOVIA Discovery Studio Visualizer分析残基水平相互作用与可视化；用GROMACS 2025.2、CHARMM27力场、SwissParam配体拓扑、TIP3P水模型和0.15 M NaCl进行50 ns MD模拟，分析RMSD、RMSF、Rg和氢键；并进行Cys285 S^γ–弹头距离时间序列分析和MM-PBSA结合自由能计算（500帧）。

研究结果如下。3.1 T0070907在PPARγ配体结合结构域中的对接与残基水平相互作用分析：通过AutoDock Vina对接、RMSD聚类、三次独立对接及PLIP/ProteinsPlus/Maestro分析，Pose 1具有最低预测结合能（?4.80 kcal/mol）并与Pose 2形成低RMSD稳定簇；T0070907与His449形成氢键（2.79 ?）和π–π堆积（3.76 ?），与Tyr327形成卤键（2.43 ?），并与Phe282、Gln286、His449形成疏水/范德华接触，表明其在预反应非共价状态下具有中等稳定化并以Pose 1为代表性构象。3.2 GW9662在PPARγ配体结合结构域中的对接与残基水平相互作用分析：通过同样对接与相互作用分析，Pose 1预测结合能最低（?5.40 kcal/mol），Pose 2与之接近（?5.30 kcal/mol）；GW9662与His449形成氢键（2.82 ?）和π–π堆积（3.75 ?），与Tyr327形成卤键（2.40 ?），并与Phe282、Phe363、Gln286、His449形成更多疏水接触，表明其预反应非共价稳定化以疏水/芳香堆积为重要特征。3.3 T0070907与GW9662在PPARγ配体结合结构域中预反应非共价稳定模式的比较：通过整合对接能量、相互作用指纹、最短相互作用距离和三次独立对接均值分析，GW9662的Vina得分（?5.4 kcal/mol）较T0070907（?4.8 kcal/mol）略更有利；二者共享His449氢键/π–π和Tyr327卤键，但GW9662疏水接触更多，T0070907极性方向性相互作用更强；统计收敛支持Pose 1，且二者在未显式建模共价键时仍保持有利稳定，提示广泛非共价网络可能足以维持预反应识别。3.4 PPARγ–配体复合物稳定性的分子动力学模拟分析：通过50 ns MD的RMSD、RMSF、Rg和氢键分析，PPARγ–T0070907复合物平均RMSD为0.6 ± 0.08 nm，低于GW9662的0.8 ± 0.2 nm；二者RMSF均值相近（0.099 ± 0.049 nm与0.098 ± 0.053 nm），Rg稳定（1.956 ± 0.009 nm与1.950 ± 0.001 nm）；T0070907氢键数均值0.4 ± 0.6，高于GW9662的0.1 ± 0.3；重复模拟定性重现主要趋势，表明T0070907具有较低骨架涨落和较频繁瞬态氢键，GW9662更灵活并以疏水作用稳定。3.5 相对于Cys285的预反应距离分析：通过时间分辨的Cys285 S^γ–弹头距离轨迹分析，T0070907在初始约7 ns弛豫后维持相对持续邻近，总体平均距离7.2 ± 1.9 ?；GW9662在约35 ns后发生明显晚期重排，最终47–49 ns和49–50 ns平均距离分别为0.66 ± 0.06 nm和0.66 ± 0.07 nm，99.0%和95.1%帧位于0.50–0.80 nm操作预反应范围内；结果表明二者经不同时间路径采样Cys285邻近构象，支持短时间尺度局部预反应稳定而非全局平衡收敛。3.6 MM-PBSA结合自由能分析：通过从MD轨迹取500帧进行MM-PBSA及能量分量分解，GW9662总结合自由能（?7.23 kcal/mol）较T0070907（?4.07 kcal/mol）更有利；GW9662具有更强ΔE_vdW（?15.83 kcal/mol）和ΔE_elec（?5.21 kcal/mol），虽ΔG_polar罚分更高（+16.18 kcal/mol），但由有利气相相互作用补偿；T0070907为ΔE_vdW ?8.56、ΔE_elec ?3.87、ΔG_polar +9.62 kcal/mol；因未计?TΔS，结果为焓主导估计（ΔH_bind），用于相对比较；结论为GW9662热力学结合更有利，而动态稳定性与结合亲和力反映不同方面。3.7 分子对接与分子动力学模拟结果的整合：通过整合对接、MD和距离分析，二者共享His449与Tyr327锚定，但T0070907以方向性极性相互作用和较低RMSD/较高氢键频率为特征，GW9662以疏水堆积和更大动态灵活性为特征；二者均维持弹头靠近Cys285的预反应取向，表明非共价相互作用足以支持瞬态稳定配体结合并保持有利预反应取向，结构相似共价拮抗剂可采用不同稳定化策略。

讨论部分总结：研究强调预反应非共价稳定化可仅由非共价相互作用维持并控制化学修饰前初始配体定位；动态轨迹稳定性与MM-PBSA热力学结合有利性并不必然一致。局限包括：预反应状态以Cys285 S^γ与弹头空间邻近操作定义，未显式模拟化学过渡态、电荷转移或活化能垒；SwissParam拓扑基于CGenFF类比，缺乏卤素σ-hole等偏心电子参数；50 ns轨迹只反映短时间尺度相互作用持续性，不能外推全局平衡，尤其GW9662存在晚期重排；未计?TΔS，MM-PBSA值为焓近似。可检验假设为：改造T0070907吡啶环以增强疏水堆积并保留His449/Cys285附近方向性极性锚定，可能提高预反应结合及共价抑制效率k_inact/K_I，需IC₅₀和时间依赖k_inact/K_I生化实验验证。

研究结论（翻译）：本研究采用整合分子对接、MD模拟和MM-PBSA自由能计算的计算工作流，全面比较了两种已充分表征的PPARγ共价拮抗剂T0070907和GW9662的预反应识别与预反应非共价稳定化。评估亲核Cys285残基显式共价修饰前的非共价结合状态，两种配体均在PPARγ LBD内采用有效结合姿态，共享His449氢键和π–π堆积以及Tyr327卤键等关键锚定相互作用；但识别出不同的预反应非共价稳定化机制：对接和相互作用指纹分析显示GW9662主要由增强的疏水和芳香接触驱动，而T0070907主要依赖方向性极性相互作用。MD模拟支持这些结构预测，显示两种复合物在50 ns生产轨迹中保持动态稳定；PPARγ–T0070907复合物显示较低主链RMSD和较高瞬态氢键占据，反映持续极性相互作用网络；PPARγ–GW9662复合物显示更大动态灵活性和较少极性接触，符合疏水驱动的预反应非共价稳定化；这些发现进一步支持动态条件下PPARγ结合口袋内持续预反应配体定位。MM-PBSA自由能分析提供定量能量见解，显示GW9662具有更有利结合自由能，主要由更强范德华和静电相互作用贡献，尽管极性溶剂化罚分更高；相反，T0070907预测相互作用能较弱，但受益于较低溶剂化罚分和更频繁瞬态极性锚定；这些发现突出轨迹来源动态稳定性（RMSD持续性）与热力学结合有利性（MM-PBSA亲和力）之间的区别，表明结构稳定性不一定直接与预测热力学结合强度相关。总体而言，这些结果表明共价拮抗剂的稳健预反应非共价稳定化可仅通过非共价相互作用维持，并控制化学修饰前初始配体定位；整合对接-MD-MM-PBSA方法为表征预反应非共价稳定化机制提供稳健生物物理框架；关键的是，这些发现严格作为可检验计算假设提出，通过前瞻性实验动力学（k_inact/K_I）和生化测定（IC₅₀）验证预测相互作用网络，将为未来共价PPARγ调节剂的理性优化提供有价值见解。