《International Journal of Molecular Sciences》：Deciphering Site-Specific Regulatory Networks of the Kinesin Protein KIF21A Through Integrative Phosphoproteomic Analysis
编辑推荐：
KIF21A是驱动蛋白(Kinesin)-4家族马达蛋白成员，参与微管(Microtubule)动力学调节及细胞内运输，现有证据提示其在癌症进展中具有潜在作用。本研究通过对超过3825个人类磷酸化蛋白质组(Phosphoproteomics)数据集的整合分析，
KIF21A是驱动蛋白(Kinesin)-4家族马达蛋白成员，参与微管(Microtubule)动力学调节及细胞内运输，现有证据提示其在癌症进展中具有潜在作用。本研究通过对超过3825个人类磷酸化蛋白质组(Phosphoproteomics)数据集的整合分析，表征KIF21A的位点特异性磷酸化。在磷酸化蛋白质组研究中鉴定出KIF21A的三个主要磷酸化位点（S853、S1212和S1239）检测频率最高，并分析其共调控模式以识别潜在的激酶关联与功能网络。磷酸化位点S853与细胞骨架组织及皮质微管稳定复合物(Cortical Microtubule Stabilization Complexes, CMSCs)组分（包括KANK1 (S186)、PHLDB2 (S513, S42)和CLASP1 (S600, S572, S646)）表现出强关联，表明其在细胞骨架组织中的作用。上游激酶分析鉴定出潜在调节因子，如PAK2、RPS6KA1/A3、RPS6KB1、CHEK1/2及CDK18/16，其与KIF21A主要位点的关联具有位点特异性变异。有趣的是，KIF21A与候选激酶间的位点特异性相关性分析显示，KIF21A S1239磷酸化位点在多种癌症类型中与CDK18呈肿瘤特异性相关。功能富集揭示共调控磷酸化蛋白参与细胞骨架调节、细胞周期调节及致癌过程。泛癌(Pan-cancer)分析证明KIF21A在多种肿瘤类型中表达失调，在选定癌症中随分期上调。基因水平验证进一步支持这些发现，显示KIF21A与关键调节因子（如CTNND1和PTK2）及其他细胞骨架和癌症相关基因呈一致正相关。总之，本研究强调位点特异性磷酸化是KIF21A的关键调节机制，并提示其参与癌症中细胞骨架相关信号网络。
研究背景方面，驱动蛋白超家族蛋白(Kinesin Superfamily Proteins, KIFs)包含多种依赖微管并利用ATP水解驱动细胞器、蛋白复合物及mRNA转运的马达蛋白，目前已鉴定超45个成员并分为14类，其中Kinesin-4家族在细胞分裂、微管组织及信号传导中起关键作用。KIF21A作为Kinesin-4亚家族成员，具有N端马达域、卷曲螺旋 stalk 区及参与货物识别与自调节相互作用的C端WD40重复域，广泛表达于人体组织尤其是中枢与外周神经系统神经元，主要参与神经元顺行胞内运输，其自抑制调节对正常功能至关重要，该界面破坏会导致马达异常激活及微管动力学失调。既往研究明确KIF21A突变导致先天性眼外肌纤维化1型(Congenital Fibrosis of the Extraocular Muscles Type 1, CFEOM1)，且近期发现该域新型突变与周围神经病变相关，但其翻译后调节及磷酸化介导的控制仍不清楚，磷酸化可影响蛋白功能、互作及定位，而KIF21A的位点特异性磷酸化调节尚不明晰，因此研究人员开展本研究以探究其磷酸化调控网络。
为开展研究，研究人员用到的主要关键技术方法包括：通过文献调研组装并筛选3825个公开人类细胞磷酸化蛋白质组数据集（分为定性谱图数据集与定量差异数据集），采用内部磷酸化位点映射流程映射位点并标准化基因符号；通过检测频率排序鉴定主要磷酸化位点；采用“共罪关联”策略结合PhosphoSitePlus数据库注释共调控其他蛋白磷酸化位点(Phosphosites in Other Proteins, PsOPs)的生物功能；利用Enrichr进行功能富集分析；通过NetworKIN与AKID等工具预测上游激酶，整合实验验证数据库及Johnson等人(2023)报道的激酶关联；利用cProSite进行癌症特异性磷酸化蛋白质组相关性分析；通过HPRD、BioGRID等多数据库构建蛋白互作网络；利用GEPIA2与UALCAN进行泛癌基因表达及分期表达分析；基于UALCAN数据进行基因水平共表达相关性分析。
研究结果部分，各小节结论如下：
2.1. Compilation of Phosphoproteomics Datasets and Identification of Major Phosphosites of KIF21A：研究人员检查3825个公开人类全局细胞磷酸化蛋白质组研究，鉴定出定性谱图数据集51个位点、差异表达数据集26个位点，其中S1239、S1212、S853在两类数据集中检测频率均最高，被确定为主要磷酸化位点，且三者均位于蛋白中央区域（N端马达域与C端WD40重复域之间），定性数据集中10个位点未在PhosphoSitePlus注释。
2.2. Evolutionary Conservation of Predominant KIF21A Phosphosites Within the Kinesin-4 Family：S1239在KIF21B保守为S1215，S1212在KIF21B保守为S1193、在KIF7保守为S1252、在KIF27保守为S1265，KIF4A/B无保守；跨物种直系同源分析显示S853与S1239在物种间保守，支持其功能与进化重要性。
2.3. Co-Regulation Analysis of Predominant KIF21A Phosphosites：经Fisher精确检验(p<0.05)、频率≥10%、至少3篇独立发表、至少3种实验条件筛选，S853正共调控1075个PsOPs、负共调控156个；S1212正428个、负20个；S1239正1016个、负131个。正共调控中S853与S1239共享91个PsOPs，S1212与S1239共享50个，S1212与S853共享16个；负共调控中S853与S1239共享4个，S1212与S853共享1个。各位点顶分享共调控PsOPs具特异性，如S853正共调控最频为NCBP1(S22)、SLC9A1(S703)、REPS1(S709)，负共调控最频为ANK3(S4350)、ACIN1(S328)、ITGA4(S1021)。
2.4. Phosphosite-Specific Biological Functions of Co-Regulated PsOPs：共调控PsOPs主要涉及致癌、细胞骨架组织、细胞周期调节及细胞运动性，推断KIF21A位点特异性功能。
2.5. Functional Enrichment of Positively and Negatively Co-Regulated PsOPs：正共调控S853、S1239富集于染色质重塑、染色质组织、DNA损伤应答；S1212富集于微管聚合/解聚、微管锚定、细胞骨架组织。负共调控S853富集于DNA损伤应答、胚胎造血、蛋白定位建立调节；S1212富集于mRNA剪接、加工及剪接体功能；S1239富集于剪接体介导mRNA剪接、转酯反应及mRNA加工。
2.6. Site-Specific Upstream Kinase Co-Regulation Associated with Predominant KIF21A Phosphosites：S853正共调控上游调节因子含PAK2(S152,S141)、RPS6KA1(S380,S732,S363,S369)、RPS6KB1(S427)、CHEK1(S280)、RPS6KA3(S227,S375)、SGK1(S401)、WNK1(S167)、SNRK(S569)、CHEK2(S260)，负共调控含PRKD3(S216)；S1239正共调控含CDK18(S132,S74)、CDK16(S95,S138)；S1212无共调控上游激酶或预测激酶。
2.7. Cancer-Specific Correlation Between the KIF21A S1239 and Phosphosites in Predicted Upstream Regulators and Carcinogenesis-Associated PsOPs：cProSite分析显示多癌种（子宫、肺、肾、脑癌）中KIF21A S1239与CDK18(S74,S132)呈中强正相关；S1239与PAK2(S141)正共调控，S1212与WDR77(T5)在多癌种（脑、肾、肺鳞癌、子宫癌）显著表达相关，S853无cProSite数据未分析。
2.8. Binary and Complex Interactors Associated with KIF21A Predominant Phosphosites：S853正共调控PsOPs含14个二元互作子、114个复合互作子，负共调控含4个二元、15个复合；S1212正共调控含4个二元、45个复合，负共调控含4个复合；S1239正共调控含10个二元、99个复合，负共调控含1个二元、18个复合。
2.9. Phosphosite-Specific Interaction and Kinase Network Analysis of KIF21A Within Cortical Microtubule Stabilization Complexes：S853与CMSCs经典组分KANK1(S186)、ERC1(S37)、PHLDB2(S513,S42)、CLASP1(S600,S572,S646)强正共调控；S1212与PPFIA1(S708,S1172)、PPFIBP1(S794)、CLASP2(S952)、CLASP1(S1091)共调控；S1239与ERC1(S21,S17)、PPFIA1(S242)、CLASP1(S1091)、CLASP2(S455,S8)共调控。上游激酶CHEK2、RPS6KA1、RPS6KA3、RPS6KB1、SGK1、WNK1、PAK2与CMSCs蛋白KANK1(S186)磷酸化相关，提示共享激酶调节轴。
2.10. Pan-Cancer Gene Expression Profiling of KIF21A：GEPIA2分析显示14种癌症KIF21A表达显著改变，11种上调（乳腺浸润癌BRCA、胆管癌CHOL、弥漫性大B细胞淋巴瘤DLBC、肾嗜酸细胞瘤KICH、肾透明细胞癌KIRP、肺腺癌LUAD、肺鳞癌LUSC、嗜铬细胞瘤PCPG、前列腺腺癌PRAD、胃腺癌STAD、胸腺瘤THYM），3种下调（急性髓系白血病LAML、子宫癌肉瘤UCS、皮肤黑色素瘤SKCM）。UALCAN分期分析显示BRCA、STAD、LUSC中KIF21A表达随分期逐步升高。
2.11. Gene Expression Correlation Analysis of Proteins Co-Regulated with KIF21A：正共调控蛋白与KIF21A在基因水平一致性相关，负共调控极少相关；S853对应正共调控蛋白基因水平一致性最高，次为S1239、S1212。最高皮尔逊相关系数见于胸腺瘤中PTK2(0.77)与CTNND1(0.74)。RPS6KA3、CHEK1、PAK2、CDK16等上游激酶及CMSCs组分KANK1、ERC1在多癌种与KIF21A正基因相关。
讨论部分总结：研究人员指出KIF21A为Kinesin-4马达蛋白，参与顺行运输与皮质微管动力学调节，与BIG1互作维持高尔基体，耗竭损害溶酶体运输，突变致CFEOM1，在乳腺癌、肺腺癌、胰腺导管腺癌中起作用。本分析中S853富集于微管稳定通路，与KANK1-PHLDB2等CMSCs组分共调控，关联皮质微管稳定；S1239与CDK18呈肿瘤相关相关，CDK18失调涉代谢病、脑缺血等，需实验验证。S1239与PAK2(S141)在肺鳞癌正关，S1212与WDR77(T5)多癌正关，PAK2促癌增殖侵袭，WDR77过表达促增殖抑免疫，为预后标志物，提示S1212/S1239关联肿瘤信号。泛癌分析KIF21A在多瘤种失调上调且与分期相关，高表达或支持细胞骨架重塑促侵袭转移，早期胰腺导管腺癌高表达预后差可作标志物。基因水平正负共调控不对称提示磷酸化信号选择性反映于转录水平，PTK2、CTNND1强相关支持KIF21A参与癌症细胞骨架信号，上游激酶及CMSCs组分KANK1、ERC1正关提示其属细胞骨架重塑网络。本研究强调KIF21A位点特异性磷酸化的重要性，位点关联与激酶互作具细胞与疾病条件依赖性，保守与癌相关位点为信号网络关键节点，为探索磷酸化微调KIF21A活性及整合入细胞骨架与信号通路提供框架，弥补KIF21A磷酸化表征不足，生成假设引导未来实验验证调控机制与疾病干预。
研究结论部分翻译：本研究中，研究人员开展全面整合磷酸化蛋白质组分析以探究KIF21A（参与微管动力学与神经元功能的Kinesin-4马达蛋白）的位点特异性调控机制。分析提示KIF21A不同磷酸化位点可能关联不同细胞过程，如细胞骨架组织、细胞内运输及通过自抑制的调节，表明其调节可能具情境依赖性而非均匀。这些发现更广泛地阐释KIF21A如何通过位点特异性调节参与多细胞通路，需进一步研究实验验证观察结果并理解单个磷酸化位点如何影响KIF21A功能，采用诱变分析与功能测定等方法将助于阐明其在细胞过程与疾病状态的作用。
