计算机药物重定位识别脾酪氨酸激酶抑制剂作为共济失调-毛细血管扩张症的潜在神经保护剂

《Chemical Biology & Drug Design》：In Silico Drug Repositioning Identifies SYK Kinase Inhibitors as Potential Neuroprotective Agents for Ataxia–Telangiectasia

　　

本研究采用药物重定位策略，通过整合基于卷积神经网络（CNN）的分子对接与SAveRUNNER网络医学分析平台，筛选2342种美国食品药品监督管理局（FDA）批准的药物库，以识别可抑制脾酪氨酸激酶（SYK）活性的临床已批准化合物。为评估候选药物在显式溶剂条件下于SYK催化口袋中的结合稳定性与构象持续性，研究人员进一步开展了长时程分子动力学（MD）模拟及模拟后重对接分析。该多级计算分析将维罗非尼和帕博西尼鉴定为最具前景的SYK抑制剂，二者在蛋白质松弛后均展现出持续结合亲和力与稳定的分子间相互作用。研究结果支持将FDA已批准药物重新用于调控SYK驱动的神经炎症通路的可行性，并为开发共济失调-毛细血管扩张症（AT） ladder的新型神经保护干预措施提供了机制框架。

本研究聚焦于共济失调-毛细血管扩张症（Ataxia-Telangiectasia, AT）这一罕见神经退行性疾病的神经炎症机制，通过计算生物学方法探索药物重定位策略的可行性。

研究背景方面，AT是一种由共济失调毛细血管扩张突变基因（ATM）突变导致的常染色体隐性遗传病，其特征为进行性神经元丢失和慢性神经炎症。ATM蛋白在DNA损伤应答（DNA damage response, DDR）中发挥关键作用，尤其参与双链断裂（double-strand breaks, DSBs）的修复。除中枢神经系统（central nervous system, CNS）功能障碍外，近年研究表明神经炎症积极参与AT的疾病进展：缺乏功能性ATM的小胶质细胞表现出内在激活状态及增强的促炎信号传导，先于神经元凋亡发生；AT患者来源的诱导性多能干细胞（iPSCs）研究证实，激活的小胶质细胞与小脑神经元共培养时可诱导细胞毒性效应。在介导该炎症反应的分子中，TYRO蛋白酪氨酸激酶结合蛋白（TYROBP，又称DAP12）在AT患者中显著过表达，其酪氨酸磷酸化形式可招募脾酪氨酸激酶（spleen tyrosine kinase, SYK），激活释放肿瘤坏死因子-α（TNF-α）和白介素-1β（IL-1β）等促炎细胞因子的效应分子。鉴于TYROBP-SYK信号轴在AT小胶质细胞驱动神经炎症中的核心作用，靶向抑制SYK成为减轻AT神经退行性成分的潜在治疗策略。

目前存在的问题在于，AT迄今尚无确定性治愈方法，治疗主要聚焦于多学科临床管理；传统从头药物开发因时间和成本限制对于罕见病往往不切实际。因此，研究人员开展了本项计算生物学研究，旨在利用整合计算策略识别能够抑制SYK的小分子化合物，从而调控TYROBP-SYK炎症轴，为AT的药物重定位提供理论依据。

研究人员采用的关键技术方法主要包括：基于卷积神经网络（CNN）的分子对接算法GNINA进行高通量虚拟筛选（HTVS）；SAveRUNNER网络医学平台进行基因组学分析以量化药物靶点与AT相关蛋白在蛋白质相互作用网络中的拓扑邻近性；长时程（200 ns）全原子分子动力学（MD）模拟评估配体-受体复合物的构象稳定性；模拟后聚类重对接（post-MD redocking）结合MM/PBSA（分子力学/泊松-玻尔兹曼表面积）方法计算结合自由能；主成分分析（PCA）和自由能景观（FEL）分析刻画蛋白质大尺度运动特征。

研究结果部分按以下小节展开：

结构依据与靶点表征方面，SYK是一种胞质非受体酪氨酸激酶，由N端两个Src同源2结构域（SH2 domains）——N-SH2和C-SH2，以及C端激酶结构域组成，其间穿插连接区A和B。在静息免疫细胞中，串联SH2结构域呈自抑制性的"连接肽-激酶三明治"构象；免疫受体激活后，SH2结构域与免疫受体酪氨酸激活基序（ITAMs）中磷酸化酪氨酸结合，触发构象重排使激酶结构域释放而激活SYK。在AT中，TYROBP下游SYK的异常激活维持了促炎信号状态，因此靶向SYK的ATP结合口袋可阻断病理性信号传导。

蛋白质模型验证方面，研究人员从PDB-Redo库获取SYK结构（PDB ID: 3EMG），该精修模型分辨率为2.60 ?，共晶配体实验测定抑制常数（Ki）为9 nM。使用GNINA算法进行共晶配体重对接后，根均方偏差（RMSD）为1.00 ?，CNN分值0.99，CNN_VS值7.77，对应预测解离常数（Kd）16.98 nM，与实验Ki值高度吻合。统计验证表明该对接方案的受试者工作特征曲线下面积（ROC-AUC）达0.90，富集因子EF1、EF5和EF10分别为10、12和7，显著优于随机选择。

虚拟筛选方面，对2342种FDA批准药物进行筛选，以CNN_VS ≥ 6.30（对应Kd ≤ 500 nM）为阈值，获得43种化合物。经关键残基相互作用分析——包括铰链区Ala451氢键、Met448/Met450/Pro455疏水或范德华相互作用，以及DFG基序Asp512氢键或静电相互作用——鉴定出奥氮平、莫洛替尼、维罗非尼、帕博西尼、伊马替尼、尼洛替尼、氟马替尼和核黄素八种满足全部结构标准的化合物。

SAveRUNNER基因组学分析层面，该平台基于网络医学原理，通过量化药物靶点与疾病相关蛋白在相互作用组中的拓扑邻近性预测药物-疾病关联。应用严格校正相似性阈值>0.96进行交叉验证后，仅维罗非尼、帕博西尼和伊马替尼同时满足GNINA结构驱动筛选和SAveRUNNER网络驱动筛选的双重标准。其余五种化合物因未达网络拓扑关联阈值而被排除，凸显了整合网络水平信息的重要性。

同源模建方面，由于可用共晶SYK结构不完整，研究人员采用同源模建构建全长模型。以UniProt ID P43405的人SYK氨基酸序列为输入，以实验解析结构为主要模板，经YASARA自动模建模块处理缺失环区，获得整体Z-分数为-0.547的全长模型用于后续MD模拟。

分子动力学模拟方面，三个候选药物-蛋白复合物均进行200 ns模拟。势能曲线显示所有系统在产博期围绕稳定均值波动，无系统性漂移。RMSD轨迹表明维罗非尼-SYK和伊马替尼-SYK复合物约10 ns后稳定，帕博西尼-SYK系统约110 ns后达平衡。RMSF分析显示结合位点残基（450-520区域，涵盖铰链区和DFG基序）波动显著降低，表明配体作为构象锁维持刚性结构锚定。溶剂可及表面积（SASA）计算证实蛋白质整体球状折叠被严格保持。氢键和疏水接触占据率监测显示，维罗非尼和帕博西尼对关键残基保持>50%的持续相互作用。

主成分分析与自由能景观方面，PCA特征值曲线证实前两个主成分（PC1和PC2）涵盖大部分协同运动。帕博西尼复合物表现卓越的构象锁定，PC1在~100 ns后进入平坦平台期；维罗非尼复合物在受限构象空间中探索离散亚态；伊马替尼则PC1和PC2均未能稳定，持续混沌漂移。FEL显示维罗非尼呈现受限明确的低能谷，帕博西尼具有深能阱且被陡峭能壁包围，而伊马替尼对应高度扰动的构象行为，与高RMSD和结合亲和力显著丧失一致。

重对接分析与MM/PBSA计算方面，从MD平衡段进行轨迹聚类提取代表性构象：维罗非尼10个、帕博西尼3个、伊马替尼11个聚类。维罗非尼和帕博西尼的平均CNN_VS分别为7.60±0.21和7.69±0.22，RMSD分别为1.34±0.44 ?和1.50±0.51 ?，MM/PBSA平均结合能分别为53.583 kcal/mol和52.448 kcal/mol。伊马替尼则CNN_VS降至2.16±1.47，RMSD高达10.53±1.66 ?，结合能仅15.218 kcal/mol，较前两者的降幅超过3.4倍。与已知SYK抑制剂基准数据集比较显示，维罗非尼和帕博西尼的后MD系综重对接分值落入实验验证的SYK抑制剂活性范围内。

血脑屏障穿透性方面，维罗非尼和帕博西尼虽为P-糖蛋白（P-gp）和乳腺癌耐药蛋白（BCRP）等外排转运体底物，但临床和临床前证据表明二者可在脑内达到药理学相关浓度：维罗非尼在BRAF突变黑色素瘤脑转移患者中显示显著颅内疗效；帕博西尼在选定脑肿瘤模型中显示出可测量的颅内抗肿瘤活性。加之AT病理背景下慢性神经炎症常伴随局部血脑屏障完整性破坏形成的"渗漏"屏障，可能进一步促进药物向受影响小脑区域的渗透。

结合模式与炎症信号影响方面，维罗非尼和帕博西尼均有效占据SYK ATP结合口袋，与铰链区和DFG基序保持稳定相互作用。这种催化裂隙的稳定占据作为ATP结合的竞争性阻断，有效阻止SYK自磷酸化及其后续激活。在TYROBP异常过表达的AT微环境中，该药理学阻断可减少活性氧（ROS）和TNF-α、IL-1β等促炎细胞因子的释放，系统性减轻小胶质细胞驱动的神经退行性变。

研究结论部分翻译如下：通过整合基于结构的分子建模与网络驱动的基因组学分析，本研究将SYK激酶抑制鉴定为减轻AT神经炎症成分的有力治疗策略。多级计算工作流证明，维罗非尼和帕博西尼展现出有效结合SYK催化位点所需的结构兼容性、动态稳定性和网络水平相关性。重要的是，长时程MD模拟和模拟后重对接的引入，使研究人员能够在构象松弛且具有生物学相关性的受体状态中验证配体结合，从而克服了静态对接模型的固有局限性。从转化医学角度，该药物重定位策略为AT等罕见病提供了实质性优势，因为传统从头药物开发往往因时间和成本限制而不切实际。尽管两种先导化合物均为已知的外排转运体底物，神经肿瘤学现有临床和临床前证据表明其可在中枢神经系统达到药理学有意义的浓度，这一能力可能因AT慢性神经炎症相关的局部血脑屏障破坏而进一步增强。当然，将维罗非尼和帕博西尼等强效抗癌药物重新用于神经退行性疾病的转化可行性，需要严格的实验验证和体内毒性管理。未来研究应采用体外酶学激酶测定实验验证靶点结合并确定SYK抑制的精确IC₅₀或Ki值；利用AT患者来源iPSC分化为小胶质细胞进行细胞功能测定；通过小胶质细胞-神经元共培养系统评估神经保护潜力；最终在ATM缺陷小鼠模型中进行体内验证，以确认靶向TYROBP-SYK轴的治疗可行性。尽管计算筛选存在固有限制，本研究为SYK靶向药物在AT中的临床调查提供了有力的机制和转化理论基础。