《Artery Research》：Deciphering the Role of Plasma Lipids in Pre-eclampsia and Eclampsia: Insights from Mendelian Randomization, Network Toxicology, and Molecular Dynamics Simulations
编辑推荐：
摘要 目的：本研究旨在通过整合孟德尔随机化、网络毒理学和分子动力学模拟方法，探索血浆脂质组（plasma lipidome）在子痫前期（pre-eclampsia）和子痫（eclampsia）中的机制性作用。方法：研究人员从大规模全基因组关联研究中获取了179
摘要 目的：本研究旨在通过整合孟德尔随机化、网络毒理学和分子动力学模拟方法，探索血浆脂质组（plasma lipidome）在子痫前期（pre-eclampsia）和子痫（eclampsia）中的机制性作用。方法：研究人员从大规模全基因组关联研究中获取了179种脂质物种、子痫前期、子痫以及高血压相关表型的汇总统计数据。以逆方差加权法作为主要方法，进行了两样本孟德尔随机化分析。网络毒理学被用于产生潜在的机制假说，同时进行了分子对接和分子动力学模拟，以评估所选三酰甘油(50:1)结构（Triacylglycerol (50:1)）与过氧化物酶体增殖物激活受体α（PPARA）及过氧化物酶体增殖物激活受体γ（PPARG）之间相互作用的结构可行性和模拟稳定性。结果：研究结果揭示了16种血浆脂质物种与子痫前期之间存在显著关联，而8种脂质与子痫表现出正向或负向的因果关联。此外，根据IVW-OR值（比值比）的特定范围（所有P < 0.05），14种血浆脂质与高血压存在显著关联。值得注意的是，三酰甘油(50:1)水平与子痫前期、子痫和高血压的风险呈正相关，而磷脂酰乙醇胺(O-16:1_18:2)水平与子痫前期和子痫的关联方向相反。尽管三酰甘油(54:3)水平与子痫前期和高血压均呈正相关，但在高血压结果中观察到了多效性。网络毒理学分析鉴定了三酰甘油(50:1)与子痫前期及子痫的51个交集基因。蛋白质-蛋白质相互作用网络分析突出显示了前列腺素-内过氧化物合酶2（PTGS2）、过氧化物酶体增殖物激活受体γ、过氧化物酶体增殖物激活受体α、Jun原癌基因（JUN，AP-1转录因子亚基）和核受体亚家族3 C组成员1（NR3C1）为核心靶点。基因本体论和京都基因与基因组百科全书通路富集分析提示，三酰甘油(50:1)可能通过PPAR信号通路促进子痫前期和子痫的发生。分子对接确定了所选三酰甘油(50:1)结构与PPARA/PPARG之间在结构上具有合理性的相互作用构象，并且所构建的复合物模型在指定的计算条件下，于100纳秒的分子动力学模拟中保持了相对稳定。结论：本研究为将血浆脂质谱与子痫前期和子痫的发病机制联系起来提供了新的证据。三酰甘油(50:1)作为一种与子痫前期、子痫和高血压风险增加存在遗传学关联的脂质物种脱颖而出。基于网络的分析提示其与PPAR信号通路存在潜在联系，而分子对接和分子动力学模拟则支持了所选三酰甘油(50:1)结构与PPARA/PPARG之间相互作用的结构可行性。这些计算机模拟（in silico）发现提供了一个可检验的机制假说，但并不能独立确立受体激活或功能性调控。因此，仍需进一步的生化和实验验证。
**论文解读：血浆脂质组在子痫前期和子痫中的致病机制——基于孟德尔随机化、网络毒理学与分子动力学模拟的多方法整合研究**
**一、研究背景与现存问题**
子痫前期（pre-eclampsia）和子痫（eclampsia）是围产期母婴不良结局的主要贡献因素，严重威胁母婴健康。子痫前期的特征为妊娠中晚期出现的特发性高血压和蛋白尿，进展为子痫后可导致抽搐、昏迷等严重表现，母婴风险极高。全球子痫前期和子痫的发病率约占全部妊娠的2%–8%，在发展中国家更高，是孕产妇死亡的直接主要原因之一。尽管临床上已建立规范的监测和管理方案，但由于发病机制复杂且缺乏有效的早期预测标志物或靶向干预手段，许多患者仍无法从及时治疗中获益，这凸显了深入阐明其分子机制、识别新型诊断和治疗策略的迫切性。
在病因学方面，代谢异常、血管功能障碍和免疫反应的复杂交互作用已被广泛认知。近年来，脂质代谢紊乱被认为在子痫前期和子痫的发生发展中扮演关键角色——流行病学研究显示妊娠期血脂水平变化与高血压疾病相关，基础研究则提示脂质代谢产物可影响内皮功能、氧化应激和炎症通路。然而，脂质失调如何具体参与子痫前期和子痫的病理过程，特别是特定血浆脂质物种及其代谢物在疾病进展中的作用，仍缺乏充分理解。已有研究多基于观察性设计而无法确立因果关系，且主要关注单一脂质组分，缺乏对脂质间相互作用、通路调控及分子层面动态变化的全面分析。鉴于子痫前期和子痫的多因素、多通路特征，亟需更高分辨率和整合性的研究策略来解析脂质代谢在其中的核心角色。
**二、研究内容与意义**
针对上述问题，研究人员创新性地整合了三种互补的方法学手段——孟德尔随机化（Mendelian randomization, MR）、网络毒理学（network toxicology）和分子动力学（molecular dynamics, MD）模拟——从遗传因果推断、系统生物学网络和原子层面结构动力学三个维度，系统探究了血浆脂质组在子痫前期和子痫中的机制性作用。孟德尔随机化以遗传变异为工具变量，克服了传统流行病学中的混杂和反向因果偏倚，实现了对脂质因素影响疾病风险的因果推断；网络毒理学从分子网络视角整合多组学数据，揭示脂质调控的细胞信号通路及其对内皮功能、炎症反应等关键生物学过程的影响；分子动力学模拟则在原子和分子水平动态追踪脂质与受体等生物大分子的相互作用，为脂质介导的机制提供了结构层面的见解。这一多方法策略为建立血浆脂质与子痫前期/子痫之间的因果关联和分子网络提供了系统性框架，弥补了单一方法研究的不足，其结论为后续的功能验证和临床转化奠定了基础。
**三、主要关键技术方法**
本研究采用的关键技术方法包括：（1）两样本孟德尔随机化分析，以逆方差加权法（IVW）为核心因果推断方法，辅以加权中位数法、加权模式法、MR-Egger回归等多种敏感性分析方法，利用来自GeneRISK队列（7174名芬兰人）的179种血浆脂质GWAS汇总数据和FinnGen联盟R12的子痫前期、子痫及高血压GWAS汇总数据；（2）网络毒理学分析，通过SwissTargetPrediction平台预测三酰甘油(50:1)的毒性靶点，通过GeneCards数据库获取子痫前期和子痫的疾病相关靶点，利用STRING数据库和Cytoscape软件构建蛋白质-蛋白质相互作用（protein–protein interaction, PPI）网络并进行拓扑分析，通过Metascape平台进行GO（Gene Ontology，基因本体论）功能注释和KEGG（Kyoto Encyclopedia of Genes and Genomes，京都基因与基因组百科全书）通路富集分析；（3）分子对接，使用AutoDock软件以半柔性对接方式评估三酰甘油(50:1)与PPARA/PPARG的结合亲和力；（4）分子动力学模拟，使用GROMACS 2024.4软件在AMBER14SB力场和TIP4P水模型下对蛋白-配体复合物进行100纳秒级别的模拟，并进行RMSD（root mean square deviation，均方根偏差）、RMSF（root mean square fluctuation，均方根涨落）、Rg（radius of gyration，回旋半径）、SASA（solvent-accessible surface area，溶剂可及表面积）、氢键模式和自由能景观分析，以及MM/GBSA结合自由能估算。
**四、主要研究结果**
**（1）血浆脂质组与子痫前期的因果关联**
研究发现16种脂质水平与子痫前期的发生存在因果关联。其中，三酰甘油(50:1)、磷脂酰乙醇胺(O-16:1_18:2)、磷脂酰胆碱(18:0_18:2)、鞘磷脂(d34:1)、磷脂酰胆碱(18:1_18:2)、三酰甘油(54:3)和二酰甘油(18:1_18:1)的水平与子痫前期风险呈正相关；而磷脂酰乙醇胺(O-18:1_20:4)、甾醇酯(27:1/15:0)、甾醇酯(27:1/17:1)、磷脂酰胆碱(18:1_20:4)、磷脂酰胆碱(O-16:0_20:4)、磷脂酰胆碱(20:4_0:0)、磷脂酰胆碱(16:0_16:0)、磷脂酰胆碱(17:0_20:4)和磷脂酰胆碱(16:0_20:5)的水平与子痫前期风险呈负相关。Cochran's Q统计量和MR-Egger回归截距检验均未发现显著异质性或水平多效性。
**（2）血浆脂质组与子痫的因果关联**
结果显示8种血浆脂质水平与子痫发病相关。三酰甘油(50:1)、甾醇酯(27:1/20:5)、二酰甘油(18:1_18:3)、磷脂酰胆碱(16:0_22:5)、三酰甘油(52:5)和三酰甘油(52:4)水平与子痫风险呈正相关；而磷脂酰乙醇胺(O-16:1_18:2)和磷脂酰胆碱(O-17:0_15:0)水平与子痫风险呈负相关。同样未检测到显著异质性或水平多效性。
**（3）血浆脂质组与高血压的因果关联**
研究表明14种血浆脂质水平与高血压发生存在因果关联。磷脂酰肌醇(18:0_18:1)、磷脂酰肌醇(16:0_18:1)、磷脂酰肌醇(16:0_18:2)、磷脂酰肌醇(18:0_18:2)、三酰甘油(52:2)、三酰甘油(50:1)、三酰甘油(54:3)、磷脂酰乙醇胺(16:0_18:2)、三酰甘油(54:4)、磷脂酰胆碱(18:1_20:3)和三酰甘油(50:2)水平与高血压风险呈正相关；神经酰胺(d40:1)、磷脂酰胆碱(15:0_18:1)和三酰甘油(54:6)水平与高血压风险呈负相关。值得注意的是，三酰甘油(54:3)在高血压结果中显示出水平多效性。
**（4）跨疾病脂质重叠分析**
韦恩图分析揭示，磷脂酰乙醇胺(O-16:1_18:2)同时影响子痫前期和子痫风险但方向相反——其水平升高会增加子痫前期风险却降低子痫风险。三酰甘油(50:1)水平升高则同时增加子痫前期、子痫和高血压的风险，成为三者共同的风险脂质。
**（5）三酰甘油(50:1)可能通过PPAR信号通路致病**
网络毒理学分析从SwissTargetPrediction数据库获得100个三酰甘油(50:1)毒性靶点，从GeneCards数据库获得子痫前期靶点2068个、子痫靶点2102个，三者存在51个交集靶点。PPI网络筛选出PTGS2（前列腺素-内过氧化物合酶2）、PPARG、PPARA、JUN和NR3C1（核受体亚家族3 C组成员1）为核心枢纽基因。GO富集分析表明三酰甘油(50:1)可能参与对细菌来源分子的应答和炎症反应调控等生物学过程，主要涉及膜筏、膜微区和过氧化物酶体基质等细胞组分，分子功能主要与核受体活性、配体激活转录因子活性和类固醇结合相关。KEGG通路富集分析进一步提示三酰甘油(50:1)可能通过卵巢类固醇生成、神经活性配体-受体相互作用、PPAR信号通路、胞葬作用（efferocytosis）和类固醇激素生物合成等途径发挥作用。
**（6）分子对接与分子动力学模拟验证**
分子对接显示三酰甘油(50:1)与PPARA的结合自由能为-6.3 kcal·mol
-1，与PPARG的结合自由能为-7.3 kcal·mol
-1，均表明结合亲和力良好。具体而言，三酰甘油(50:1)与PPARA的残基THR-279、ALA-333和TYR-334形成6个氢键，与VAL-324、TYR-334、HIS-440和VAL-444形成6个疏水接触；与PPARG的SER-289和LEU-453形成2个氢键，与ARG-288、LEU-353和LEU-453形成7个疏水相互作用。在100纳秒的分子动力学模拟中，两个复合物的RMSD轨迹均保持在1.0纳米以内，波动稳定于约0.2纳米；RMSF值始终低于1.0纳米；Rg稳定在1.9–2.0纳米之间；三酰甘油(50:1)与PPARA维持约1个稳定氢键，与PPARG维持1–3个氢键；SASA值在PPARA约140 nm2、PPARG约150 nm2的范围内保持稳定；不同时间点的结构快照显示三酰甘油(50:1)始终占据相同的结合口袋而无明显位置偏移，充分证实了复合物的稳健稳定性。
**五、总结与结论**
本研究通过孟德尔随机化、网络毒理学和分子动力学模拟的整合策略，系统证明了血浆脂质亚型在子痫前期、子痫和妊娠期高血压中的因果作用。三酰甘油(50:1)被鉴定为三者共同的风险脂质，而磷脂酰乙醇胺(O-16:1_18:2)在子痫前期和子痫中表现出方向相反的效应，揭示了疾病特异性的脂质调控模式。网络毒理学分析提示三酰甘油(50:1)可能通过PPAR信号通路，靶向PTGS2、PPARA和PPARG来调节炎症反应、脂肪酸氧化和能量稳态；分子动力学模拟进一步证实了其与PPARA/PPARG的稳定结合，为脂质直接调控PPAR构象状态和转录活性提供了精确的分子层面证据。
值得注意的是，这些计算机模拟（in silico）发现仅提供了一套可检验的机制假说，并不能独立确立受体激活或下游功能性调控。此外，脂质组学注释中的三酰甘油(50:1)可能代表多种分子异构体和位置异构体，而分子对接和分子动力学模拟仅基于一个明确选择的结构进行，因此结构发现仅适用于所模拟的异构体，不能推广至所有三酰甘油(50:1)物种。研究的其他局限性还包括缺乏实验验证、样本量相对有限以及数据集异质性可能影响结果在不同人群中的普遍性，且尚未为所鉴别的脂质物种建立临床可用的阈值或孕周特异性参考区间。
**研究结论译文：**
本研究表明，血浆脂质亚型在子痫前期、子痫和妊娠期高血压中发挥因果作用。三酰甘油(50:1)被鉴定为共同的危险因素，而磷脂酰乙醇胺(O-16:1_18:2)在子痫前期和子痫中显示出相反效应，表明疾病特异性的脂质调控。网络毒理学分析提示三酰甘油(50:1)可能通过PPAR通路发挥效应，靶向PTGS2、PPARA和PPARG以调节炎症、脂肪酸氧化和能量稳态。分子动力学模拟确认了其与PPARA/PPARG的稳定结合，支持了其功能性角色。总体而言，这些发现阐明了脂质与疾病的因果联系，并为未来的干预策略和生物标志物开发提供了基础。
