突变诱导的pocket失活：Ser353/Pro245如何改变KCa2.2与KCa3.1的配体选择性

《Archiv der Pharmazie-Chemistry in Life Sciences》：Mutation-Induced Pocket Deactivation: How Ser353/Pro245 Alters KCa2.2 vs. KCa3.1 Ligand Selectivity

　　

KCa2.2和KCa3.1通道是膜电位和钙信号传导的基本调节因子，也是治疗脊髓小脑性共济失调和癌症等疾病的有前景靶点。为充分利用其治疗潜力并继续研究其病理生理作用，开发针对这两种通道的选择性调节剂至关重要。本研究展示了一项计算研究，旨在识别两种近期报道的KCa2.2调节剂（即N-(2,1,3-苯并噁二唑-4-基)-3-(4-甲氧基苯-1-磺酰胺基)苯甲酰胺和N-(2,1,3-苯并噁二唑-4-基)-4-(三氟甲基)苯甲酰胺）选择性背后的分子决定因素。研究人员利用结合了计算机模拟诱变、分子动力学模拟（MD）和蛋白-配体对接的方案，分析了这些配体靶向的口袋。研究人员发现，Ser353/Pro245替换是KCa2.2和KCa3.1通道中不同口袋形状的主要驱动因素，最终决定了调节剂的选择性。该方法提供了关于该结合位点在钾通道亚型间结构差异的新见解，提出了靶向该口袋的调节剂的潜在选择性决定因素。

**研究背景、问题与目的**
钾离子（K⁺）选择性通道通过调节细胞对K⁺离子的通透性来塑造膜电位，从而控制激素分泌、细胞兴奋性、上皮运输和细胞体积变化等基本生理功能。其中，小电导钙激活钾通道（KCa2.1、KCa2.2、KCa2.3，也称SK1-3、KCNN1-3）和中电导钙激活钾通道（KCa3.1，也称IK、SK4、KCNN4）在胞质Ca²⁺浓度变化时其K⁺电导发生改变，建立了局部离子浓度的紧密偶联。这些通道由同源四聚体构成，每个亚基包含6个跨膜螺旋（S1-S6）和3个额外的胞内螺旋（HA、HB、HC）。钙敏感性由辅助蛋白钙调蛋白（CaM）介导，其C-末端叶（C-Lobe）与每个通道亚基的HA和HB螺旋组成性结合。
KCa2.2和KCa3.1通道的调节具有广泛治疗潜力：KCa2.2抑制剂对心房颤动和阿尔茨海默病有效，分子AP30663已完成II期临床试验；KCa2.2激活剂可用于治疗脊髓小脑性共济失调和酒精依赖，或在卒中发作中作为神经保护剂；KCa3.1抑制剂在镰状细胞性贫血和癌症治疗中显示出前景，阻断剂senicapoc接近与perampanel联合用于胶质母细胞瘤的临床试验；KCa3.1激活剂在治疗高血压和囊性纤维化中具有潜力。然而，KCa2.2和KCa3.1的多样病理学意义要求开发选择性调节剂，以充分利用其治疗潜力并研究其功能。例如，KCa2.2激活剂对中枢神经系统疾病有益，但若对外周KCa3.1通道缺乏选择性，可能导致低血压等血管副作用；反之，抑制KCa3.1治疗胶质母细胞瘤时，若同时阻断KCa2.2通道，可能导致癫痫和神经退行性变。
近期，通过冷冻电镜（cryo-EM）解析了KCa2.2通道与新型小分子抑制剂（化合物1，N-(2,1,3-苯并噁二唑-4-基)-3-(4-甲氧基苯-1-磺酰胺基)苯甲酰胺）和激活剂（化合物4，N-(2,1,3-苯并噁二唑-4-基)-4-(三氟甲基)苯甲酰胺）的复合物结构，两者均结合于先前未描述的S5-Phelix-S6口袋。化合物1抑制KCa2.2的IC50为69 nM，且对KCa3.1具有10倍选择性（IC50 660 nM）。基于此，研究人员旨在确定该S5-Phelix-S6口袋是否可用于开发新型KCa3.1调节剂，并深入理解KCa2.2靶向化合物选择性的分子基础。通过结合计算模拟方法，研究人员揭示了KCa2.2与KCa3.1的S5-Phelix-S6口袋关键差异，并识别出可能决定调节剂选择性的特定残基。
本研究发表在《Archiv der Pharmazie-Chemistry in Life Sciences》。

**主要关键技术方法**
研究人员采用计算机模拟诱变（in silico mutagenesis）、分子动力学模拟（MD）和蛋白-配体对接（分子对接，使用GOLD软件）的组合方法。首先，利用KCa2.2和KCa3.1的cryo-EM结构（PDB ID: 9O53、9O5O、6CNM、9OA8）进行序列和结构比较。随后，对KCa3.1通道进行Pro245Ser计算机模拟诱变，并对野生型（WT）和突变型（P245S）系统进行三重复500 ns的MD模拟。最后，从MD轨迹中提取代表性构象，进行ensemble对接（综合对接）以评估配体结合能力。

**研究结果**
**2.1 S5-Phelix-S6口袋在KCa2.2与KCa3.1之间的比较**
通过对比KCa2.2（化合物1结合态）和KCa3.1（闭合态）中S5-Phelix-S6口袋的残基，发现23个残基中有9个序列替换，但最关键的是Trp322/Trp216的构象差异：在KCa2.2中，Trp322侧链朝向胞外，允许配体结合；在KCa3.1中，Trp216侧链朝向胞内，阻塞口袋，导致化合物1和4无法结合。此外，Ser353/Pro245替换位于Trp322/Trp216附近，计算机模拟将Ser353突变为Pro后，与Trp322产生严重冲突，提示Pro245可能限制Trp216的构象。该Ser/Pro替换在物种间保守，且与通道最大开放概率（POMAX）差异相关（KCa2.2为0.8，KCa3.1为0.1-0.2）。
结论：Trp322/Trp216构象差异和Ser353/Pro245替换是决定化合物1和4选择性的关键因素。

**2.2 计算机模拟诱变与分子动力学分析**
**2.2.1 蛋白稳定性与折叠**
RMSD和回转半径（RoG）分析表明，Pro245Ser突变未引起KCa3.1通道不稳定或去折叠。
**2.2.2 Trp216构象变化**
在KCa3.1_open_P245S模拟中，Trp216在一条链中发生自发构象转换，从初始的向下构象（构象a）转变为向上构象（构象b和c），与KCa2.2中Trp322的构象高度相似；而在WT模拟中，Trp216保持向下构象。该转换仅在三分之一的重复中出现，表明其低概率性。
结论：Pro245Ser突变使Trp216能够发生构象变化，从“失活”态转变为“接受”态。
**2.2.3 Trp216相互作用分析**
在WT系统中，Trp216通过疏水接触和π-π堆积与Leu213、Phe248、Leu249和Leu268稳定；在P245S突变系统中，这些相互作用丧失，转而与Ser245和Leu241形成氢键，稳定了向上构象。
结论：Pro245Ser突变通过改变Trp216与P-Helix的距离，使其能够采用新的稳定构象。

**2.3 分析KCa3.1_open_P245S系统中的S5-Phelix-S6口袋**
对KCa3.1_open_P245S轨迹进行口袋体积分析，发现口袋体积从396 ?3增加并稳定在约700 ?3，最大体积1240 ?3，与KCa2.2中口袋体积（1403 ?3）相当。对化合物4的ensemble对接结果表明，80%的对接姿态与cryo-EM结合模式一致，所有前十排名姿态均符合实验模式，且最佳结合发生在Trp216向上构象（构象c）的帧中。对化合物1的对接显示，只有4个姿态与核心骨架结合模式类似，但甲氧基苄基部分偏离，这与其抑制剂特性（需稳定闭合态）一致。
结论：Pro245Ser突变使KCa3.1的S5-Phelix-S6口袋能够容纳化合物4，且该口袋的构象变化与配体结合能力直接相关。

**讨论与结论**
本研究通过计算模拟揭示了KCa2.2和KCa3.1通道S5-Phelix-S6口袋的差异，提出Ser353/Pro245替换是化合物1和4选择性的主要决定因素。在KCa3.1中，Pro245限制了Trp216的构象，使其处于向下位置，阻碍配体结合；而Pro245Ser突变使Trp216能够采用向上构象，恢复口袋的“接受”状态。计算结果强调了该结合位点在设计KCa2.2选择性调节剂中的潜力，但仍需实验验证。
结论部分翻译：总之，本研究提供了关于KCa2.2和KCa3.1通道S5-Phelix-S6口袋配体选择性的机制见解，提出Ser353/Pro245替换是主要的选择性决定因素。虽然需要实验验证来确认这些计算预测，但本研究进一步强化了靶向该位点以开发对KCa3.1通道无活性的KCa2.2调节剂的潜力。