《Journal of Genetic Engineering and Biotechnology》：Phosphosite-specific co-regulation networks of MELK kinase: insights from integrative global phosphoproteomes
编辑推荐：
母体胚胎亮氨酸拉链激酶（Maternal embryonic leucine zipper kinase, MELK）是一种丝氨酸/苏氨酸激酶，在侵袭性癌症中频繁过表达，但其激活和信号传导特异性的精确机制仍不清楚。本文通过对3,825个全局人类磷酸化蛋白质组学
母体胚胎亮氨酸拉链激酶（Maternal embryonic leucine zipper kinase, MELK）是一种丝氨酸/苏氨酸激酶，在侵袭性癌症中频繁过表达，但其激活和信号传导特异性的精确机制仍不清楚。本文通过对3,825个全局人类磷酸化蛋白质组学数据集进行整合性荟萃分析，首次构建了MELK的磷酸化位点分辨共调控图谱。三个磷酸化位点——S356、S505和S529——作为主要调控节点浮现，表现出高检测频率和独特的共调控模式。S356和S505形成一个由MAPK（Mitogen-Activated Protein Kinase Kinase）、RSK（Ribosomal S6 Kinase）、CaMK（Calcium/Calmodulin-dependent protein Kinase）、Hippo相关激酶以及纺锤体检查点激酶（NEK4、TTK/MPS1）汇聚控制的紧密偶联的增殖-有丝分裂轴，而S529则作为部分拮抗性的应激和极性响应模块发挥作用。增殖标志物Ki-67（MKI67）的磷酸化位点在几乎所有增殖背景下均与三个MELK位点共变，从而在MELK活性与临床增殖标志物之间建立了直接的机制联系。共调控的上游激酶、磷酸酶、二元互作蛋白和下游底物的广泛网络进一步揭示了功能分离：S356/S505主要驱动细胞周期进程和染色质组织，而S529整合钙信号、代谢信号和细胞骨架极性信号。对TCGA（The Cancer Genome Atlas）队列的Kaplan–Meier生存分析进一步揭示，MELK、MKI67和有丝分裂检查点激酶TTK的高表达一致预测肺腺癌和肝细胞癌患者的总生存期和无病生存期较差，这强化了MELK活性与临床增殖标志物之间强烈的磷酸化动态偶联。通过证明MELK信号传导是通过模块化、位点特异性的磷酸化逻辑而非总蛋白丰度来协调的，这项工作为理解和治疗靶向这一神秘的致癌激酶建立了新范式。
研究背景与问题
母体胚胎亮氨酸拉链激酶（Maternal embryonic leucine zipper kinase, MELK）是一种保守的丝氨酸/苏氨酸激酶，属于AMPK（AMP-activated protein kinase）相关激酶亚家族，在细胞周期进程、增殖和应激响应中发挥关键作用。该激酶在乳腺癌、结直肠癌、胶质母细胞瘤和前列腺癌等多种侵袭性癌症中频繁过表达，与肿瘤增殖、治疗抵抗和不良预后密切相关。尽管MELK已成为一个有前景的治疗靶点，并且其抑制剂OTSSP167已进入早期临床试验，但关于其激活和信号传导特异性的精确分子机制，尤其是在磷酸化位点水平的调控，仍不清楚。现有数据库如PhosphoSitePlus虽然记录了MELK的磷酸化位点和相互作用，但缺乏来自差异性磷酸化蛋白质组学数据集的共表达分析。鉴于磷酸化可能因条件不同而变化并塑造信号传导结果，因此，深入解析MELK关键磷酸化位点的特异性调控网络，对于理解其在肿瘤发生发展中的作用、发现新的生物标志物和精准治疗机会至关重要。
研究目的与意义
本研究旨在通过整合分析大量公开的全局人类磷酸化蛋白质组学数据集，构建首个MELK磷酸化位点分辨的共调控图谱，以阐明其功能调控的分子逻辑。研究人员假设MELK的关键主要位点存在不同的磷酸化模式，这些模式可能与肿瘤表型、进展和治疗抵抗相关，而磷酸化位点分辨的相互作用图谱将揭示新的调控因子、生物标志物和精准肿瘤学机会。该研究发表在《Journal of Genetic Engineering and Biotechnology》期刊上，其意义在于超越对MELK总蛋白丰度的关注，从磷酸化动态角度揭示其作为模块化信号枢纽的复杂性，为开发位点特异性靶向疗法和动态磷酸化生物标志物提供理论基础。
主要技术方法
研究人员开展了一项大规模的整合性荟萃分析。首先，系统检索并筛选了3,825个人类细胞质谱磷酸化蛋白质组学公共数据集，从中识别出与MELK相关的数据。为确保可靠性，仅纳入定位概率≥75%、A-score≥13的I类磷酸化位点。通过构建定制化分析流程，将基因符号和磷酸化位点统一映射至标准数据库（如HGNC、UniProt）。利用检测频率分析确定了MELK的主要磷酸化位点。通过共现分析和共调控分析，探究了MELK不同磷酸化位点之间及其与其他蛋白磷酸化位点（PsOPs）的关联。进一步，整合了来自多个数据库的激酶-底物关系、蛋白-蛋白相互作用数据，并结合计算预测工具，构建了MELK磷酸化位点的上游激酶、磷酸酶、二元互作蛋白和下游底物的调控网络。最后，利用GEPIA2平台，基于TCGA数据库的RNA-seq和临床生存数据，对MELK及其相关基因（MKI67, TTK）在乳腺癌、肺腺癌和肝细胞癌中的预后价值进行了Kaplan-Meier生存分析。
研究结果
3.1. MELK的全局磷酸化蛋白质组图谱和主要磷酸化位点
通过对数据集的分析，共编译出68个独特的MELK磷酸化位点，其中21个在不同条件下呈现差异调控。磷酸化位点S356、S505和S529是检测频率最高的位点，被确定为主要调控位点。可视化分析显示，这些主要位点聚集在激酶结构域之外，提示它们在非催化调控功能（如相互作用、定位或底物特异性）中的重要性。
3.2. MELK磷酸化位点的共现模式
共现分析揭示了MELK磷酸化位点间存在显著的正向和负向相关性模式。S356和S505之间存在强烈的正向共现，表明这两个位点几乎同时被磷酸化。S529与S505、T415与T409之间也存在中等程度的正向关联。相反，T460与S356表现出强烈的负相关，S529与T466、T518呈中等负相关。这些模式支持了MELK调控中存在协同（S356-S505）和部分拮抗模块的假说。
3.3. 与MELK主要位点共调控的磷酸化蛋白（PsOPs）
研究人员鉴定出大量与MELK主要磷酸化位点高度共调控的PsOPs。通路富集分析显示，与S356共调控的PsOPs显著富集于染色质重塑、细胞周期调控、Notch/ErbB信号通路等。S505相关网络则广泛富集于DNA复制、细胞周期进程、AMPK/MTOR信号通路和多种癌症相关通路。S529相关的PsOPs则选择性富集于同源重组、范可尼贫血通路、黏着斑、肌动蛋白细胞骨架调控等通路。这表明S356、S505和S529形成了三个功能不同的MELK信号模块，分别侧重于染色质控制、增殖信号传导以及DNA修复-细胞骨架编程。
3.4. MELK磷酸化位点的预测上游激酶
预测分析发现了14个可能调控MELK主要磷酸化位点的上游激酶。S505似乎是多种激酶家族的汇聚点，包括MAPK、RSK、CaMK和Hippo相关激酶家族的成员。相比之下，S356的调控似乎更特异性地依赖于预测的有丝分裂检查点激酶NEK4和TTK（MPS1）。S529则由钙/钙调蛋白依赖性激酶和应激相关激酶调控。这揭示了MELK磷酸化受到一个分层级激酶网络的控制。
3.5. 与MELK主要磷酸化位点共调控的激酶
共调控分析识别出大量与MELK各主要位点共变的激酶。其中，CIT、EEF2K、TYK2、PI4KA和PRKAR1B五个激酶的磷酸化位点在所有三个MELK位点中均表现出一致的共调控，表明存在核心信号模块。此外，S356和S505共享29个共调控激酶，包括MAP3K1、MARK2、LIMK1等，定义了S356–S505共磷酸化轴。S529则与更多应激和极性响应激酶相关联。
3.6. 与主要磷酸化位点正向共调控的磷酸酶
研究鉴定了与MELK各主要位点正向共调控的磷酸酶。四个磷酸酶（PTPN2、MTMR2、PGAM1、CDC25C）在所有三个MELK位点中均出现，构成MELK去磷酸化的核心调控因子。不同位点还关联了特异的磷酸酶，包括酪氨酸磷酸酶、双特异性MAPK磷酸酶、细胞周期磷酸酶、脂质定向磷酸酶等，表明MELK的磷酸化状态受到细胞周期、MAPK、代谢和脂质信号网络中的位点特异性磷酸酶的精细调节。
3.7. 二元互作蛋白磷酸化位点共调控网络识别核心和位点特异性MELK互作蛋白
通过映射与MELK主要位点显著共变的互作蛋白磷酸化位点，构建了MELK的二元互作网络。一个核心发现是，七个MKI67（Ki-67）的磷酸化位点作为全位点共调控因子出现，将这一关键增殖标志物确立为MELK在所有激活状态下的核心伙伴。此外，不同位点有其独特的互作蛋白：S356与有丝分裂检查点激酶TTK、应激激酶MAP3K5等关联；S505与染色质调节因子EZH2、剪接因子CDC5L等关联；S529则与应激适应和极性相关蛋白如CAMK2G、GIGYF1等关联。这些互作将MELK定位为连接有丝分裂进程、RNA/染色质生物学、细胞骨架动力学和适应性信号传导的磷酸化调控枢纽。
3.8. 共调控蛋白中的下游底物
研究人员鉴定了166个与MELK主要位点共调控的潜在下游底物。一个核心底物三联体（CDC25C、ZNF395、RAPGEF6）在所有三个MELK位点中均被共调控。位点特异性底物分析显示，S356和S505主要正向调控细胞周期调节因子、染色质修饰因子等，驱动促癌程序；而S529则显示出混合的诱导/抑制效应，富含应激适应、凋亡抑制和细胞骨架极性相关底物。功能富集分析进一步证实了这种功能分离。
3.9. MELK、MKI67和TTK的生存分析
基于TCGA数据库的Kaplan-Meier生存分析显示，在肺腺癌（LUAD）和肝细胞癌（LIHC）中，MELK、MKI67和TTK的高表达均与较差的总生存期（OS）和无病生存期（DFS）显著相关。其中，在LIHC中，三者的预后关联最为强烈。然而，在乳腺癌（BRCA）队列中，这三个基因的表达与生存期无显著统计学关联。这表明MELK及其共调控网络在特定癌症类型中具有显著的临床预后价值。
讨论与结论总结
在讨论部分，研究人员指出本研究通过最大的整合性磷酸化蛋白质组学荟萃分析，揭示了MELK的功能是由三个主要调控磷酸化位点（S356、S505、S529）以模块化方式协调的。S356和S505构成一个紧密偶联的有丝分裂-增殖轴，而S529作为一个应激和极性响应传感器发挥部分拮抗作用。MKI67（Ki-67）被确定为跨所有位点的通用磷酸化动态伙伴，从而在机制上将MELK活性与细胞增殖联系起来。共调控的磷酸酶、底物和互作蛋白网络解释了MELK如何同时加速细胞周期进程并维持适应性检查点。生存分析将鉴定的磷酸化位点分辨的增殖模块转化为临床上相关的预后特征。这些发现阐明了MELK在癌症中背景依赖性必需性的长期悖论，为其致癌多能性提供了分子框架。
研究结论翻译
本研究通过对超过3,800个全局人类磷酸化蛋白质组学数据集的整合分析，首次建立了母体胚胎亮氨酸拉链激酶（MELK）的磷酸化位点特异性调控图谱。我们证明MELK信号传导并非由均匀的磷酸化控制，而是由三个主要调控位点S356、S505和S529运作，它们作为功能不同但相互关联的模块。S356和S505在MAPK、RSK、CaMK、Hippo和纺锤体检查点激酶（NEK4、TTK/MPS1）的汇聚控制下形成一个紧密偶联的有丝分裂-增殖轴，而S529则作为具有部分拮抗输出的应激和极性响应传感器。MKI67（Ki-67）作为跨所有三个位点的通用磷酸化动态伙伴出现，从而在机制上将MELK活性与细胞增殖联系起来，而CDC25C-S216以及广泛的磷酸酶、底物和二元互作蛋白网络揭示了MELK如何同时加速细胞周期进程并维持适应性检查点。重要的是，对TCGA队列的生存分析表明，高MELK和MKI67表达，连同共调控的有丝分裂激酶TTK，与肺腺癌和肝细胞癌的不良预后相关，从而将所鉴定的磷酸化位点分辨的增殖模块转化为临床相关的预后特征。这些发现阐明了MELK背景依赖性必需性在癌症中的长期悖论，解释了其在某些模型中尽管存在争议性的“非必需性”但仍与不良预后频繁相关的原因，并为其致癌多能性提供了潜在的分子框架。临床上，磷酸化位点分辨图提名了用于肿瘤增殖的动态磷酸化生物标志物，突出了与有丝分裂检查点抑制剂的合成致死机会，并强烈推动了开发位点选择性或模块特异性MELK疗法的动机。
