基于极性杂环变构抑制致癌性JAK2 V617F假激酶域：热力学解卷积与跨结构域选择性谱

《Pharmaceuticals》：Allosteric Inhibition of the Oncogenic JAK2 V617F Pseudokinase Domain via Polar Heterocycles: Thermodynamic Deconvolution and Cross-Domain Selectivity Profiles

　　

专业摘要翻译 背景：由Janus激酶2（JAK2）假激酶（JH2）结构域的体细胞V617F突变驱动的组成型过度激活，破坏了调控相邻催化（JH1）结构域的关键自身抑制界面。由于标准疗法靶向JH1结构域上结构保守的正构ATP结合口袋，它们常诱发严重的全激酶组脱靶毒性。本研究建立了一个严格的生物物理学框架，以定义实现跨结构域选择性（避开催化JH1结构域而倾向于假激酶JH2结构域）所需的热力学和结构参数，并结合针对假激酶裂隙内V617F取代的突变特异性靶向。方法：研究人员利用卷积神经网络增强的分子对接技术，针对野生型JH1和突变型JH2晶体结构筛选了一个经计算过滤的药物样杂环化合物库。通过在独立的三次重复200 ns显式溶剂分子动力学轨迹中量化结构持久性和动态稳定性。结果：分子力学广义玻恩表面积（MM-GBSA）自由能解卷积确定候选化合物306为首选的计算预测突变选择性锚定物，其总结合自由能为-32.13 kcal/mol。该亲和力主要由强范德华相互作用（-41.11 kcal/mol）结合高度刚性、深层的自由能盆地驱动。动态互相关矩阵映射验证了候选化合物306通过机械同步P-loop、铰链区和DFG基序来抑制局部热波动。同时，卤代候选化合物253在突变结合口袋内预测了定向卤素介导的接触网络，这得益于深度学习评分函数对卤素-受体几何取向的参数化。结论：克服历史激酶数据库偏差需要靶向极性优化和形状互补堆积以实现稳定的蛋白质失活，为未来的实验验证提供了可靠的预测基础。

论文解读文章

研究背景与科学问题

Janus激酶2（JAK2）是非受体酪氨酸激酶家族的核心成员，负责转导红细胞生成素（EPO）、血小板生成素（TPO）等造血相关细胞因子信号。其C端串联排列着活性激酶结构域（JH1）和催化失活的假激酶结构域（JH2）。在生理状态下，JH2作为负调控模块，通过与JH1的物理相互作用抑制其基础磷酸化活性。然而，体细胞点突变V617F（第617位缬氨酸被苯丙氨酸取代）位于JH2的SH2-JH2连接区或结合口袋内，导致自身抑制界面破坏，引发JH1结构域不依赖配体的组成型过度活化，进而驱动下游STAT蛋白持续磷酸化和恶性基因转录，这是费城染色体阴性骨髓增殖性肿瘤（MPNs）的主要致病机制，在真性红细胞增多症中检出率超过95%。

当前临床一线治疗药物如ruxolitinib和fedratinib均为靶向JH1正构ATP结合口袋的I型或II型抑制剂。由于ATP结合位点在人类激酶组中具有极高的序列和结构同源性，这类药物常导致严重的脱靶毒性、剂量限制性骨髓抑制和贫血。此外，长期用药不可避免地诱导催化环上的继发耐药突变。因此，将治疗靶点从正构JH1转向变构JH2，特别是突变后的V617F JH2结构域，成为实现真正突变特异性治疗的关键突破。然而，历史上积累的生物活性数据（如ChEMBL数据库）几乎全部来自对JH1催化活性的测量，导致传统虚拟筛选算法偏向于发现正构激酶铰链结合剂，难以识别假激酶口袋所需的特殊拓扑特征。这一研究空白亟需一种能够实现跨结构域骨架跃迁（Cross-Domain Scaffold Hopping）的高级计算框架。

研究方法概述

研究人员首先构建了一个包含10万个新颖小分子的生成式化学库，并通过多层级计算机吸收、分布、代谢、排泄和毒性（ADMET）及合成可及性级联过滤器，将其缩减至376个精英结构支架。随后，利用XGBoost机器学习模型进行定量构效关系（QSAR）分类，并通过SHAP和平均杂质减少（MDI）特征重要性分析揭示决策逻辑。最终，采用卷积神经网络（CNN）增强的GNINA 2.0分子对接框架，并行针对野生型JH1（PDB: 6VN8）和突变型JH2（PDB: 8BAK）晶体结构进行筛选。对选定的四种化学原型（候选物227、253、299、306）执行了独立三次重复、每次200纳秒的显式溶剂分子动力学（MD）模拟，并利用MM-GBSA进行热力学解卷积、二维自由能景观（FEL）分析、动态互相关矩阵（DCCM）及拉马钱德兰图评估。

研究结果

2.1 机器学习性能指标与加性特征归因地形图：优化的XGBoost模型在独立测试集上实现了0.884的ROC-AUC，灵敏度0.7867，特异性0.8531，马修斯相关系数（MCC）为0.6405。MDI分析显示模型优先使用摩根指纹位点（MorganFP）进行早期结构路由，而SHAP分析表明连续理化描述符（如RDKit_Chi0）主导最终概率输出的精细调整，形成二分学习策略。

2.2 平行分子对接与跨结构域选择性谱：CNN增强对接结果显示，候选物306和253展现出对突变型JH2的绝对拓扑偏好。候选物306在突变JH2中获得Vina结合亲和力-8.45 kcal/mol、CNN姿势置信度0.9349，而在野生型JH1中遭遇+2.65 kcal/mol的不利对接惩罚。候选物253在突变JH2中获得CNN置信度0.9812，但在JH1中产生+21.52 kcal/mol的严重空间冲突。对照候选物227表现出无差别的高亲和力，而候选物299则完全倒转了选择性向量，偏好野生型JH1。

2.3 多重复轨迹解卷积与热力学能量景观分析：在突变JH2口袋内，选择性候选物306和253表现出快速的结构松弛和低均方根偏差（RMSD，约2.56±0.61 ?），并建立了紧密的疏水界面，RMSD与溶剂可及表面积（SASA）呈显著负相关。相反，非选择性对照候选物227在突变JH2中经历了完全的解离，RMSD峰值达57.05 ?。二维吉布斯自由能景观证实，选择性候选物形成了狭窄、深层的能量井，而非选择性候选物的能量景观则宽泛而碎片化。

2.4 终点MM-GBSA热力学解卷积与结合亲和力谱：在突变JH2中，候选物306的总结合自由能（ΔG_bind）最优，达到-32.13 kcal/mol，主要由极强的范德华力贡献（ΔE_vdW = -41.11 kcal/mol）驱动，有效平衡了极性溶剂化惩罚（ΔG_{polar_solv} = +14.71 kcal/mol）。候选物253的ΔG_bind为-26.88 kcal/mol。尽管对照候选物227的静态焓值看似有利，但其动态轨迹已证明其实际脱离口袋，揭示了连续介质溶剂化模型的局限性。

2.5 动态互相关矩阵（DCCM）解卷积与环锁定机制：DCCM分析显示，选择性候选物306在突变JH2中将反相关节点向量压制至极低水平（仅5对），同时维持284个相关残基向量，实现了结构的刚性锁定。候选物253同样限制了反相关运动。而非选择性候选物227在突变JH2中显示出碎片化的互相关拓扑，大量反相关节点（118对）反映了配体的滑动和逃逸。在野生型JH1中，所有复合物均表现出高灵活性的协同耦合。

2.6 全局热力学自由能景观（FEL）与拓扑盆地映射：选择性候选物306的全局能量最小值稳定在极窄的空间坐标（RMSD约1.52-1.89 ?，R_g约19.15-19.37 ?），表明其将突变假激酶结构域冻结在紧凑的非活性构象中。候选物253也表现出类似的深层热力学锁定。相比之下，非选择性候选物227在突变JH2中的能量景观缺乏封闭的陡峭能量阱，允许复合物轻易逃逸。在野生型JH1中，所有复合物均呈现宽阔但稳定的能量通道，表明活性口袋容纳依赖高构象适应性。

2.7 蛋白质主链二面角稳定性与拉马钱德兰相空间分析：拉马钱德兰图分析表明，所有复合物的主链φ和ψ二面角均保持在高度允许区域，未引入非生理性扭曲。选择性候选物306和253保持了正常的整体骨架柔韧性（约14.58°-14.74°），而非选择性候选物227在突变JH2中因与Phe617的空间冲突而表现出压缩的ψ方差（13.43°），揭示了导致其解离的结构应变。

结论与意义

本研究成功建立了一个用于跨结构域骨架跃迁的生物物理学框架，以规避历史上针对正构JH1结构域的激酶组结合偏差。通过将治疗焦点转向V617F突变JH2假激酶结构域的结构改变拓扑，研究人员表征了实现长期自身抑制所需的物理和热力学参数。多重复显式分子动力学模拟和二维吉布斯自由能景观分析证实，实现选择性变构动力学捕获依赖于紧密的形状互补堆积，而非松散的几何口袋进入。

通过MM-GBSA自由能解卷积和动态互相关矩阵分析，候选物306被验证为精英突变选择性锚定物，其通过优化的范德华界面和机械锁定侧翼P-loop、铰链区及DFG环架构，实现了-32.13 kcal/mol的结合自由能。同时，候选物253展示了另一种由高度定向的卤素介导非共价网络驱动的变构稳定机制，有效密封了突变裂隙以防止水渗漏。相比之下，非选择性内部对照在变构腔室内经历了完全的结构崩溃和逃逸，验证了结构过滤器的严格性。这些见解为设计下一代突变特异性JAK2候选结构提供了去风险的原子蓝图，旨在最小化与传统正构疗法相关的脱靶毒性和贫血的结构决定因素。