MEK2 Glu207Lys突变在CFC4综合征中的机制性后果：融合调节性A-helix的结构建模、重复分子动力学模拟与MEK抑制剂选择性解析

《Journal of Molecular Graphics and Modelling》：In silico insights into the mechanistic consequences of the MEK2 Glu207Lys mutation in the CFC4 syndrome

　　

计算机模拟（in silico）研究揭示了MEK2 Glu207Lys突变在CFC4综合征（cardiofaciocutaneous syndrome type 4, CFC4）中的机制性后果。研究人员构建了包含调节性A-helix的MEK2结构模型；Glu207Lys破坏Arg53-Glu207盐桥（salt bridge），使MEK2的A-helix失稳；分子动力学（molecular dynamics, MD）模拟显示，在六种MEK抑制剂中，结合模式具有同种型（isoform）和突变依赖性；trametinib成为针对MEK2 Glu207Lys最有前景的抑制剂；trametinib缺乏选择性支持基于同种型/突变信息的抑制剂设计。

研究背景：CFC4综合征（cardiofaciocutaneous syndrome type 4, CFC4）是RASopathies中由生殖系突变导致Ras/MAPK信号通路失调的罕见发育障碍，常与BRAF、KRAS、MAP2K1和MAP2K2等基因突变相关；MAP2K2编码MEK2（mitogen-activated protein kinase kinase 2）。MEK1与MEK2激酶域序列相似性近90%，但MEK2缺乏包含调节环的高分辨率晶体结构，限制了对突变诱导构象变化及选择性抑制剂结合的分析。MEK1与MEK2功能并不冗余，CFC4由突变MEK2驱动而MEK1保持野生型（wild-type, WT）且具生理必要性，长期使用非选择性MEK1/2抑制剂存在不必要的MEK1抑制与毒性风险。Glu207Lys位于或邻近负调节域helix A，替换带负电谷氨酸为带正电赖氨酸，可能破坏锚定helix A的Arg53-Glu207盐桥；其MEK1类似物Glu203Lys有实验表征，为模拟预测提供外部基准。因此研究人员构建含helix A的MEK2结构，比较MEK1、MEK2 WT与MEK2 Glu207Lys的突变机制和六种MEK抑制剂结合，以指导同种型和突变选择性抑制剂设计。

研究开展与结论意义：研究人员用同源建模、诱导契合对接（induced fit docking, IFD）和重复分子动力学（molecular dynamics, MD）模拟，发现Glu207Lys可重复破坏Arg53-Glu207盐桥并使helix A失稳，激活环Ser222/Ser226局部相互作用网络重组但未表现一致失稳；在六种抑制剂中曲美替尼（trametinib）在突变体中保留最完整、可重复的相互作用指纹和最低MM-GBSA预测结合自由能，但缺乏MEK2相对MEK1选择性。该研究为CFC4相关MEK2变异的变构调节机制和抑制剂选择提供计算框架。

关键方法：研究人员以MEK2晶体结构（PDB ID: 1S9I）和MEK1结构（PDB ID: 3EQI）拼接构建含helix A的嵌合模板，用MODELLER生成WT与Glu207Lys模型，按最低DOPE（Discrete Optimized Protein Energy）评分选模，并用MolProbity Ramachandran分析验证；以LigPrep和Protein Preparation Wizard准备配体与蛋白，采用IFD对接六种MEK抑制剂；用Desmond、OPLS4力场和TIP3P水模型进行NPT系综MD，300 K、1.01325 bar、2 fs步长、1000 ns，每个系统三个独立重复；用VMD和MDTraj分析轨迹，包括Shrake-Rupley SASA（solvent-accessible surface area，溶剂可及表面积）、Baker-Hubbard氢键和残基距离，并用MM-GBSA（molecular mechanics generalized Born surface area）计算ΔG_bind。本研究无样本队列来源。

研究结果：2.1 Construction of WT and Glu207Lys mutant MEK2 models incorporating helix A：研究人员将MEK1的helix A移植到MEK2催化核心并校正序列差异，生成WT和Glu207Lys模型；Ramachandran分析显示WT模型94.6%残基位于 favored 区、4个离群残基均位于非目标环区，突变模型为94.5%且同一环区离群，提示突变未引入额外局部立体化学张力；MEK1与MEK2叠合总体RMSD（root mean square deviation，均方根偏差）为0.432 ?。2.2 MEK2 Glu207Lys mutation – in silico insight into the gain of function mechanism：重复MD显示突变体helix A骨架RMSD均值2.07 ?（WT）对3.73 ?（突变体），Arg53(N)-Glu207(O)平均距离7.75 ?对17.73 ?，突变体helix A远离激酶核心并部分解旋，说明电荷反转破坏盐桥锚定。Dynamics of the activation loop near the phosphorylation sites (Ser222/Ser226)：激活环骨架RMSD未显示一致突变效应，且差异可追溯至能量最小化后初始构象差异；Ser222在突变体三个重复中均更暴露于溶剂，Ser226无一致趋势；氢键网络发生局部重塑，Ser222获得Glu106和Ile107新接触并失去与Ser226直接氢键，Ser226失去Val228/Gly229接触并获得Glu209接触，作者将此类氢键观察表述为描述性和方向性。2.3 Docking Analysis of known MEK inhibitors：IFD显示常见接触包括Phe209/213的π-π堆积、Asp208/212氢键、Asn78/82氢键、Val131卤键、Lys97/101氢键/卤键/阳离子-π及Ser212/216氢键；binimetinib和selumetinib在突变体中采取更浅且倒置的姿态，trametinib和mirdametinib在突变体中结合模式几乎相同，cobimetinib和CI-1040姿态相似但卤代芳香环角度不同。2.4 MD Analysis：Protein-ligand interaction fingerprints显示Lys97/101氢键在MEK1中最可重复，在MEK2中呈配体和重复依赖；突变体中binimetinib、mirdametinib、selumetinib和trametinib保留该接触但稳定性低于MEK1和WT MEK2；Phe209/213的π-π接触在突变体中仅trametinib于三个重复中恢复至中等稳定性；Val127/131卤键在突变体中部分恢复，trametinib三个重复均恢复但一个不稳定，selumetinib在两种MEK2背景中均缺失；trametinib在突变体中恢复Lys101、Phe213、Val131三个核心接触，为最完整可重复指纹。MM-GBSA binding free energy：多数配体在MEK1中ΔG_bind最有利、突变体最不利，例外为CI-1040和trametinib；trametinib在三系统中ΔG_bind显著更负，突变体为?93.72±3.39 kcal/mol，MEK1为?88.72±7.15 kcal/mol，WT为?82.14±2.05 kcal/mol；cobimetinib的WT至突变体位移最大（?55.00±10.90至?45.26±9.05 kcal/mol），作者强调该值为半定量相对指标。Ligand RMSD and SASA：多数重复形成紧凑稳定簇，例外包括selumetinib和CI-1040的MEK2 WT重复及binimetinib、cobimetinib、mirdametinib、selumetinib的突变体重复；CI-1040和trametinib最一致，binimetinib与selumetinib的RMSD位移方向相反。Conformational Stability and Torsional Analysis：配体扭转角PCA（principal component analysis，主成分分析）显示MEK1中多数配体高度刚性，突变体使binimetinib、mirdametinib、selumetinib和cobimetinib构象异质性增加，CI-1040和trametinib仍保持紧凑突变体簇，trametinib簇最紧密且突变体重新占据两个类MEK1紧凑簇。

讨论与结论：讨论指出trametinib的优势排序与已报道的MEK1 Glu203Lys类似突变的体外效价变化大致一致，IC₅₀（half-maximal inhibitory concentration，半数抑制浓度）对trametinib约增加20倍、binimetinib约70倍、cobimetinib近140倍；仅保留相互作用网络而缺乏姿态稳定性并不保证效价保留。trametinib的显著局限是缺乏MEK2相对MEK1偏好，长期给药可能带来不必要的MEK1相关毒性；因此反对对CFC4相关MEK2突变谱采用“一刀切”抑制剂策略，支持按变异体进行突变和同种型选择性指导的抑制剂选择，并建议开展Glu207Lys功能获得机制的生化确认及trametinib、mirdametinib和cobimetinib对MEK2 Glu207Lys与WT MEK2和MEK1的IC₅₀检测。研究结论部分翻译：该研究提供了一个基于MD的分子框架，用于理解MEK2 Glu207Lys突变如何重塑激酶的调节架构及其变构抑制剂结合景观。MEK2 Glu207Lys突变体可重复地破坏锚定自抑制helix A的Arg53-Glu207盐桥，但未对相邻激活环产生同样一致的影响；在该区域，Ser222/Ser226动态发生局部重组，而非一致性失稳。为识别五种已批准MEK抑制剂和研究性化合物CI-1040中哪些最可能保留对Glu207Lys突变体的活性，研究人员比较了其在MEK1 WT、MEK2 WT和MEK2 Glu207Lys中的相互作用占据、MM-GBSA能量学和构象稳定性。trametinib在每一项单独分析中均成为最有前景的化合物，在MEK2 Glu207Lys中持续保留类WT特征。Mirdametinib和cobimetinib也在突变体中保留了核心相互作用网络，但结合姿态稳定性明显较低：mirdametinib的RMSD/PCA离散度更大，cobimetinib具有显著重复异质性及最大的不利MM-GBSA位移。该排序与已报道的MEK1 Glu203Lys类似突变的体外效价变化大致一致；其中IC₅₀对trametinib约增加20倍，对binimetinib约增加70倍，对cobimetinib增加近140倍；后者说明，仅保留相互作用网络而缺乏伴随的姿态稳定性，并不能保证效价得以保留。trametinib的一个显著局限是其对MEK2相对于MEK1缺乏偏好，长期给药可能转化为不必要的MEK1相关毒性。结合MEK1 Glu203Lys所报道的化合物特异性IC₅₀变化，这反对针对CFC4相关MEK2突变谱采用“一刀切”的抑制剂策略，而支持针对每个变异体采用突变和同种型选择性指导的抑制剂选择。直接的实验随访将在两方面加强这些结论：Glu207Lys功能获得机制本身的生化确认，以及trametinib、mirdametinib和cobimetinib对MEK2 Glu207Lys与WT MEK2和MEK1的IC₅₀检测，以确定突变体效价和MEK2选择性。