PCK2作为结节性硬化症候选代谢节点的整合转录组学与结构研究

《Current Issues in Molecular Biology》：PCK2 Emerges as a Candidate Metabolic Node in Tuberous Sclerosis Complex: An Integrative Transcriptomic and Structural Study

　　

由TSC1/TSC2突变驱动的结节性硬化症（TSC）导致持续的mTORC1激活和频繁的癫痫发作。mTORC1过度活跃如何引发神经元过度兴奋仍知之甚少。本研究旨在揭示TSC相关癫痫的核心基因和候选治疗药物。研究人员将多组学生物信息学方法（包括基因集变异分析、加权基因共表达网络分析、随机森林和SHapley Additive exPlanations）应用于TSC转录组数据集。经过跨数据集验证后，预测了上游转录因子，并通过分子对接和分子动力学模拟评估了候选药物。内质网应激通路被富集；PCK2被识别为关键核心基因，ATF4为其上游调控因子。卡马西平被预测可与PCK2的Gly254、Gly255和Ser257结合，形成相对稳定的复合物，结合自由能为?37.61 kcal/mol。PCK2可能连接TSC癫痫中的mTORC1功能障碍、内质网应激、线粒体代谢和神经元损伤。除了钠通道作用外，卡马西平可能靶向PCK2，这有待实验验证。

研究背景与科学问题

结节性硬化症（TSC）是一种常染色体显性遗传的多系统疾病，约1/6000新生儿受累，由TSC1或TSC2基因致病性变异引起。其最突出的神经系统表现是癫痫，发生于80%–90%的患者中。TSC1和TSC2蛋白形成的复合物抑制mTOR复合物1（mTORC1），而TSC1/2复合物的缺失导致mTORC1过度激活，这是TSC大部分临床表现的基础。尽管mTORC1抑制剂（如雷帕霉素和依维莫司）对TSC相关肿瘤效果显著，但对神经症状的疗效有限——即使依维莫司获批用于TSC难治性癫痫，许多患者仍经历每日或近乎每日的癫痫发作，且对神经认知缺陷无显著改善。因此，亟需探索mTORC1的下游机制，特别是内质网（ER）应激介导的神经元炎症和细胞死亡在TSC相关癫痫中的作用。mTORC1过度激活导致蛋白质合成过量，进而增加未折叠或错误折叠蛋白的积累，诱发ER应激。ER应激通过未折叠蛋白反应（UPR）的三条分支（PERK/eIF2α通路、ATF6通路和IRE1/XBP1通路）感知并响应。持续或未解决的ER应激可转化为“适应不良”或“终末”UPR，过度激活炎症和细胞死亡通路，参与多种疾病发病。虽然ER应激升高已在TSC细胞模型中报道，但其在TSC患者致痫组织（如皮质结节）中的具体作用尚未系统阐明。基于此，研究人员假设ER应激可能通过特定基因调控网络介导TSC相关癫痫中的神经元功能障碍。

主要研究方法

研究人员整合了来自基因表达综合数据库（GEO）的转录组数据集，包括TSC伴癫痫患者的皮质结节样本以及诱导多能干细胞（iPSC）来源的神经前体细胞（NPC）数据集。采用基因集变异分析（GSVA）评估ER应激相关通路富集情况，利用加权基因共表达网络分析（WGCNA）和机器学习方法（随机森林RF和SHapley Additive exPlanations SHAP）筛选关键基因。通过四个在线数据库预测上游转录因子（TF），并使用HOCOMOCO数据库进行评分验证。最后，通过比较毒理基因组学数据库（CTD）筛选靶向PCK2的小分子，运用CB-DOCK2进行分子对接，并使用GROMACS进行全原子分子动力学（MD）模拟，通过MM/GBSA方法计算结合自由能。

研究结果

3.1 加权基因共表达网络分析（WGCNA）： 研究人员整合GSE62019和GSE16969两个数据集（共14个样本，7例TSC vs. 7例对照），选择软阈值β=12构建共表达网络，得到7个模块。其中黄色模块与TSC相关性最高（r=0.89，p=2×10^-5）。留一法敏感性分析显示PCK2所在模块与疾病的关联在所有迭代中均为正相关，PCK2保持高基因显著性（GS）和模块成员度（MM）值。全样本分析中，调整数据集后PCK2所在模块与疾病的关联仍显著（β=0.478，p=3.72×10^-5，FDR=2.98×10^-4）。

3.2 差异表达分析与功能富集分析： 鉴定出1795个差异表达基因（DEG），其中上调1007个、下调788个。GSVA分析显示，TSC组中“ER应激反应的调控”等生物过程显著上调。

3.3 核心基因筛选与验证： 将黄色模块核心基因与DEGs交集，再与GSE239412和GSE277844数据集的DEGs重叠，获得34个共同上调基因。随机森林分析排名前十的基因包括PCK2，SHAP分析也识别PCK2为重要特征。五个重叠基因（含PCK2）在独立数据集GSE93325中验证，仅PCK2表达趋势与前期分析一致。

3.4 转录因子预测与相关性分析： 通过四个数据库预测PCK2的上游TF，经HOCOMOCO评分筛选出ATF4和KLF6。相关性分析显示ATF4与PCK2呈最强正相关（R=0.77，p<0.05），表明ATF4可能是PCK2的上游调控因子，且该调控关系已有文献报道。

3.5 小分子优先化： 通过CTD数据库筛选，保留两种临床批准药物卡马西平和丙戊酸进行后续分析。

3.6 蛋白质结构优化： 由于缺乏适合分子对接的PCK2解析结构，研究人员对Chai-1预测的蛋白结构进行了100 ns MD模拟优化。RMSD曲线显示60 ns后趋于稳定（约0.3 nm），表明蛋白构象达到平衡。

3.7 分子对接分析： CB-DOCK2对接结果显示卡马西平与PCK2的对接分数（?9.4）远优于丙戊酸（?5.0），提示更有利的相互作用。

3.8 分子动力学模拟： PCK2-卡马西平复合物的MD模拟显示，RMSD在前40 ns波动后逐渐稳定于约0.15 nm；整体RMSF较低；回旋半径（Rg）维持在2.5 nm以内；溶剂可及表面积（SASA）稳定在250–275 nm²；氢键数量稳定在1–2个。100 ns结构的相互作用可视化显示卡马西平与PCK2的Ser257、Gly254和Gly255形成氢键。MM/GBSA计算预测结合自由能（ΔE_Bind）为?37.61 kcal/mol。

讨论与结论

TSC患者中约83.6%发生癫痫，mTORC1持续激活是其核心分子基础。本研究发现ER应激相关生物过程在TSC中显著激活，PCK2被识别为跨数据集一致失调的候选基因，且在癫痫发作前的样本中即已存在，提示PCK2失调可能与TSC1/2缺陷本身相关而非癫痫继发改变。ATF4被预测为PCK2上游TF，提示存在潜在的ER应激-ATF4-PCK2调控轴。PCK2编码线粒体磷酸烯醇式丙酮酸羧激酶，催化草酰乙酸转化为磷酸烯醇式丙酮酸，是连接线粒体碳代谢、生物合成和氧化还原稳态的重要代谢节点。在营养剥夺或ER应激下，ATF4诱导的PCK2表达最初可能维持细胞存活，但持续激活可能导致病理性代谢重编程，包括碳通量紊乱、氧化还原失衡、乳酸代谢异常和炎症信号增强，这些因素可能协同促进癫痫网络的形成。

药物筛选发现卡马西平可能是PCK2的结合分子。分子对接和MD模拟从计算层面支持两者形成相对稳定的复合物。然而，这些计算预测不能证明生理条件下的直接结合或确定卡马西平对PCK2的功能效应（抑制、激活或变构调节），仍需生物物理结合实验、酶学研究和细胞实验验证。此外，卡马西平诱导CYP3A4，而与TSC常用药依维莫司联用时可能加速后者代谢，这一潜在药代动力学相互作用值得关注。

研究局限性： 分析基于公共转录组数据集，存在组织来源、脑区、年龄、癫痫状态和用药史等方面的异质性；ATF4-PCK2关系尚未在TSC神经细胞中直接验证；PCK2失调是因果性、贡献性还是代偿性尚不清楚；分子对接仅为计算机预测，不能证明直接结合或体内功能效应。

研究结论： 本研究将PCK2确定为TSC相关癫痫的潜在关键基因，提出ER应激和ATF4介导的转录可能参与其失调。PCK2可能连接TSC/mTOR异常与线粒体代谢重编程、神经炎症和神经元功能障碍。卡马西平除已知的钠通道作用外，还可能通过与PCK2相互作用发挥作用，为机制性药物重定位提供了可能依据。这些结论仍属假说生成性质，需在细胞模型、动物模型和临床样本中验证。