《Biomolecules》：Receptor-Specific Molecular Recognition of Maresin 1 and Docosahexaenoic Acid by the Nuclear Receptors PPARα, PPARγ, and RORα
编辑推荐：
摘要专业翻译
背景/目的：Maresin 1（MaR1）是一种来源于二十二碳六烯酸（DHA）的促消退介质，其与DHA在C7/C14羟基化模式和双键排布上存在差异。研究人员在代谢功能障碍相关脂肪性肝病（MASLD）背景下，比较了过氧化物酶体增殖物激活受体α（PP
摘要专业翻译
背景/目的：Maresin 1（MaR1）是一种来源于二十二碳六烯酸（DHA）的促消退介质，其与DHA在C7/C14羟基化模式和双键排布上存在差异。研究人员在代谢功能障碍相关脂肪性肝病（MASLD）背景下，比较了过氧化物酶体增殖物激活受体α（PPARα）、过氧化物酶体增殖物激活受体γ（PPARγ）和维甲酸相关孤儿受体α（RORα）对二者的识别。方法：研究人员通过重对接参考配体验证了Glide对接方法，并对每个复合物进行了四次独立的250 ns分子动力学模拟（24条轨迹；总计6.0 μs）。研究人员评估了结构稳定性、重复样本解析的接触模式、氧原子距离分布和构象采样；相对能量估算采用分子力学/广义Born表面积（MM-GBSA）计算获得。结果：重对接RMSD值≤1.5 ?，重复轨迹均未出现配体解离。在PPARα中，MaR1表现出反复出现的C7-OH–Lys358邻近接触以及可重复的羧酸根相关Tyr314和Ser280接触。在PPARγ中，MaR1羧酸根采样了螺旋12（helix 12）中的Tyr473和HIE323，且几何构型具有重复样本依赖性，而DHA则保留了以Arg288为中心的锚定模式。在RORα中，MaR1优先采样以Arg370为中心的排布，而DHA主要与Arg367结合；两种配体均表现出重复样本依赖性的替代极性接触。与DHA相比，MaR1在所有三对受体中均表现出更负的氢键、亲脂性和范德华贡献值，以及更低的数值化总MM-GBSA评分。结论：MaR1表现出与替代羧酸根定位相关的受体特异性极性接触再分布。尽管这些计算并不能确定结合亲和力或受体激活，但这些模式——而非保守的直接C7-OH/C14-OH识别基序——为匹配类似物研究和功能实验定义了可检验的相互作用位点。
论文主体内容解读
一、研究背景与科学问题
代谢功能障碍相关脂肪性肝病（MASLD）是系统性代谢紊乱在肝脏的表现，也是全球最常见的慢性肝病之一。该疾病以脂质代谢紊乱、胰岛素敏感性下降、线粒体功能障碍和炎症信号异常为特征，可进展为代谢功能障碍相关脂肪性肝炎（MASH）、纤维化、肝硬化和肝细胞癌。尽管resmetirom近期获批用于特定非肝硬化MASH患者，但MASLD的药物治疗选择仍然有限。
核受体在脂质代谢和肝脏代谢稳态中发挥关键作用。PPARα主要在肝细胞中表达，促进脂肪酸β氧化和脂质分解；PPARγ在生理条件下主要表达于脂肪组织，但在脂肪变性肝细胞中表达上调，参与脂质摄取、储存、葡萄糖稳态和胰岛素敏感性调控；RORα则通过单体识别ROR响应元件（RORE），调节肝脏脂质稳态，肝脏特异性敲除RORα可通过线粒体功能障碍加重实验性脂肪性肝炎。
DHA（C22:6 ω-3）可调节PPAR依赖性通路，并作为包括MaR1在内的特异性促消退介质（SPMs）的生物合成前体。MaR1是DHA衍生的脂质介质，具有抗炎、抗氧化和组织保护作用。Han等人已鉴定MaR1为RORα的内源性配体，并描述了促消退的RORα–12-脂氧合酶（12-LOX）–MaR1回路。尽管DHA和MaR1共享C22:6多不饱和骨架，但MaR1在C7和C14含有两个额外羟基，并具有重组的共轭三烯系统，这些差异可能影响配体识别。
然而，已有研究在不同生物学和方法学背景下进行，未能在统一、重复的分子动力学工作流程中直接比较MaR1和DHA在PPARα、PPARγ和RORα三种受体中的识别差异。因此，MaR1的额外极性官能团和改变的双键排布是否产生保守的识别基序或受体依赖性相互作用模式仍不清楚。研究人员假设，MaR1的结构差异将导致极性相互作用和能量贡献以受体依赖性方式重新分布。
二、研究设计与主要技术方法
研究人员采用了一套综合计算流程。首先，选取了三个高分辨率晶体结构作为受体模板：PPARα（PDB ID: 2P54；1.79 ?，共结晶GW590735和SRC1肽）、PPARγ（PDB ID: 6MS7；1.43 ?，共结晶(S)-VSP-77和PGC1共激活肽）和RORα（PDB ID: 1N83；1.63 ?，共结晶胆固醇）。使用Schr?dinger Suite 2020-3中的Protein Preparation Wizard制备蛋白结构，LigPrep制备配体（MaR1和DHA均以去质子化羧酸根形式，pH 7.4）。采用Glide SP–XP分层对接流程进行分子对接，并通过重对接共结晶参考配体验证协议（RMSD均≤1.5 ?）。随后使用Desmond进行四次独立的250 ns分子动力学模拟（每个复合物4条轨迹，共24条轨迹，约6.0 μs），OPLS3e力场，TIP3P水模型，0.15 M NaCl。轨迹分析包括Cα RMSD/RMSF稳定性评估、蛋白–配体接触占有率分析（氢键、离子、疏水和水桥）、残基–配体氧原子距离分布、自由能景观（FEL）和主成分分析（PCA），以及MM-GBSA能量分解计算。
三、主要研究结果
3.1 对接验证与初始结合模式：重对接RMSD值分别为GW590735（0.36 ?）、胆固醇（0.50 ?）和(S)-VSP-77（1.06 ?），均低于1.5 ?阈值。MaR1在所有三个受体中均获得比DHA更低的数值化XP GScore（PPARα：?10.174 vs ?5.422 kcal·mol?1；PPARγ：?10.502 vs ?9.434 kcal·mol?1；RORα：?10.165 vs ?9.058 kcal·mol?1）。对接结果显示MaR1具有更广泛的预测极性或离子相互作用网络。
3.2 重复分子动力学模拟支持稳定的受体–配体复合物：Cα RMSD在50–100 ns内达到稳定平台，均值范围为1.93 ± 0.25 ?（RORα/DHA）至2.49 ± 0.13 ?（PPARγ/DHA）。所有轨迹均未出现配体解离。平均α-螺旋含量为56.43 ± 0.45%至61.51 ± 0.37%，表明所有复合物在整个模拟时间内保持了核受体配体结合结构域（LBD）的α-螺旋折叠。
3.3 接触谱分析揭示受体特异性识别模式：在PPARα/MaR1中，Lys358、Tyr314、Ser280和Phe273在全部四个重复中均达到≥50%占有率标准，Lys358–C7-OH平均距离为3.12 ± 0.31 ?，Tyr314–羧酸根和Ser280–羧酸根距离分别为2.72 ± 0.13和2.87 ± 0.11 ?，均在4/4重复中可重复。在PPARα/DHA中，HIE440、Tyr314和Ser280在4/4重复中达到占有率标准，Tyr464/H12在3/4重复中达到标准。在PPARγ/MaR1中，Tyr473和HIE323在3/4重复中达到占有率标准，羧酸根邻近在Tyr473的2/4和HIE323的3/4重复中几何可重复。在PPARγ/DHA中，Arg288在4/4重复中达到97.9 ± 1.7%占有率和2.82 ± 0.12 ?羧酸根距离，而Tyr473仅显示重复依赖性的疏水接触。在RORα/MaR1中，Arg370在3/4重复中达到84.9 ± 30.2%占有率，羧酸根距离3.02 ± 0.53 ?；在RORα/DHA中，Arg367在3/4重复中达到72.7 ± 34.1%占有率，Arg370在2/4重复中达到标准。
3.4 FEL和PCA构象采样分析：在PPARγ中，MaR1产生相对紧凑的低表观相对自由能区域（全局最小值：RMSD = 0.193 nm，Rg = 1.899 nm），而DHA采样了更宽、更高RMSD的区域（RMSD = 0.258 nm，Rg = 1.916 nm），这是三对受体中最显著的配体依赖性差异。RORα和PPARα的配体相关差异较小。余弦含量分析显示PPARγ/DHA在所有四个重复中PC1和PC2均表现出升高的时间相关行为，而RORα/MaR1的高RMSD延伸主要由单个重复驱动。
3.5 MM-GBSA能量分析：MaR1在所有三对受体中均获得比DHA更低的数值化总MM-GBSA评分（RORα：?68.94 ± 2.45 vs ?61.91 ± 2.44 kcal·mol?1；PPARα：?76.49 ± 1.43 vs ?69.46 ± 2.85 kcal·mol?1；PPARγ：?76.95 ± 2.15 vs ?69.92 ± 2.84 kcal·mol?1）。MaR1表现出更负的氢键（?3.86至?3.50 vs ?3.34至?2.99 kcal·mol?1）、亲脂性（?27.44至?26.56 vs ?25.57至?24.69 kcal·mol?1）和范德华贡献（?55.51至?53.64 vs ?53.14至?51.27 kcal·mol?1），而DHA表现出更低的极性溶剂化贡献。
3.6 表观FEL最小值的代表性结构：RORα/DHA代表性结构显示Arg370邻近配体取向，但重复水平分析确定Arg367为更频繁再现的羧酸根锚点。PPARα/DHA快照将配体定位在涉及Tyr464、Tyr314和Ser280的经典极性区域。PPARγ/DHA保留Arg288中心排布。PPARγ/MaR1快照显示HIE323和Tyr473附近的羧酸根邻近。
四、讨论总结与结论翻译
讨论部分：本研究的核心发现是，MaR1骨架的结构差异——C7/C14羟基化和改变的双键位置与几何构型——与受体特异性的极性锚定再分布、构象采样和MM-GBSA能量组分变化相关。原子水平分析不支持跨三种受体的保守直接C7-OH/C14-OH识别基序。MaR1的适配以受体特异性羧酸根定位为特征，仅在PPARα中检测到C7-OH–Lys358的反复羟基特异性邻近。在PPARγ中，MaR1采样了涉及HIE323和AF-2/H12残基Tyr473的羧酸根排布，而DHA保持以Arg288为中心的锚定。在PPARα中，DHA形成涉及HIE440、Tyr314、Ser280和Tyr464/H12的可重复经典极性网络，而MaR1结合Tyr314和Ser280羧酸根相关接触与反复的C7-OH–Lys358邻近。在RORα中，MaR1优先采样Arg370，DHA主要结合Arg367，两种配体均接触Gln289。MM-GBSA分解显示MaR1在所有三对受体中均获得更低的数值化总评分和更负的氢键、亲脂性及范德华贡献。这些结果被解释为模型依赖的相对能量趋势，而非定量结合亲和力差异。
局限性：所有轨迹均源于单一选择的XP姿态，无法穷尽替代配体取向、缓慢H12转变或稀有受体状态。PPARγ/DHA和RORα/DHA仅返回一个XP姿态。仅使用分离的配体结合结构域，排除了DNA结合、共激活因子招募、翻译后修饰和PPAR–RXR异二聚化。模拟从具有已解析AF-2/H12构象的配体结合模板启动，可能限制了可访问的局部构象系综。未纳入apo受体和已确立的全激动剂或部分激动剂作为参考对照。
结论翻译：这项比较计算研究表明，MaR1和DHA在PPARα、PPARγ和RORα中通过受体特异性而非保守的识别模式被适配。在PPARα中，MaR1结合反复的C7-OH–Lys358邻近与Tyr314和Ser280相关羧酸根接触，而DHA结合涉及HIE440、Tyr314、Ser280和Tyr464/H12的经典极性网络。在PPARγ中，MaR1以重复样本依赖性几何构型采样HIE323和Tyr473，而DHA主要通过Arg288保持锚定。在RORα中，MaR1优先采样Arg370，DHA采样Arg367，尽管两种配体均可访问替代的精氨酸中心排布并接触Gln289。MaR1在所有三对受体中均获得更低的数值化总MM-GBSA评分和更负的氢键、亲脂性及范德华贡献。这些结果被解释为模型依赖的相对能量趋势，而非定量结合亲和力差异。原子水平分析未识别保守的直接C7-OH/C14-OH识别基序，反复的羟基特异性邻近仅限于PPARα中的C7-OH–Lys358。由于MaR1和DHA还在双键位置和几何构型上存在差异，需要匹配的结构类似物来解析每个分子特征的贡献。所识别的残基共同定义了在MASLD相关系统中进行匹配结合测量、定点突变和受体功能实验的可检验位点。
