布洛芬结合并抑制微管亲和力调节激酶4（MARK4）的生物物理与计算证据

《Biomolecules》：Biophysical and Computational Evidence for Ibuprofen Binding to and Inhibition of MARK4

　　

MARK4是一种丝氨酸/苏氨酸激酶，与神经退行性疾病和癌症密切相关，使其成为治疗干预和药物重定位的有吸引力的靶点。在本研究中，研究人员采用稳态荧光和紫外-可见光谱、圆二色光谱、分子对接以及分子动力学模拟等综合实验与计算方法，研究了布洛芬对MARK4的抑制作用及其相互作用。ATPase实验表明，布洛芬对MARK4的ATPase活性产生剂量依赖性抑制。MARK4的ATPase活性被抑制超过93%，半数抑制浓度（IC50）为4.02 ± 0.33 μM。荧光滴定实验显示，加入布洛芬后MARK4的内源荧光发生猝灭。热力学分析表明结合是自发的，具有有利的熵变，并且疏水作用力参与了结合过程。观察到的椭圆率降低表明布洛芬结合后MARK4的二级结构发生了改变。分子对接得到的结合能为?7.0 kcal/mol。布洛芬被对接于ATP或抑制剂结合位点内，形成两个氢键。分子动力学模拟证实了MARK4-布洛芬复合物在100 ns内的动态稳定性。与游离蛋白相比，该复合物保持了整体紧凑性，均方根偏差（RMSD）和均方根波动（RMSF）仅有微小增加。这些发现证明了布洛芬抑制MARK4 ATPase活性的潜力，这可能进一步指导药物发现。

一、研究背景与研究目的

微管亲和力调节激酶（MARKs）属于丝氨酸/苏氨酸激酶亚家族，在细胞骨架动力学、神经元极性和细胞内信号传导中发挥重要作用。其中，MARK4是微管稳定性和tau蛋白磷酸化的关键调节因子，其功能失调与阿尔茨海默病等神经退行性疾病以及肿瘤进展、代谢紊乱和炎症反应密切相关，因此成为有前景的治疗靶点。布洛芬作为广泛使用的非甾体抗炎药（NSAID），传统上通过抑制环氧化酶（COX）发挥作用，但流行病学研究表明长期使用可能降低阿尔茨海默病风险。尽管已有多种MARK4抑制剂如黄酮类芹菜素和临床药物加兰他敏被报道，但它们普遍存在生物利用度差、脱靶效应或缺乏临床验证等问题。相比之下，布洛芬具有良好的药代动力学和安全性特征。然而，布洛芬是否能够直接与MARK4相互作用并抑制其活性尚不明确。为此，研究人员采用生物物理和计算方法，系统研究了布洛芬对MARK4的结合与抑制作用。

二、主要技术方法

研究人员使用了以下关键技术方法：ATPase抑制实验采用孔雀绿磷酸盐检测法测定布洛芬对MARK4 ATPase活性的抑制效果；紫外-可见吸收光谱监测配体结合引起的蛋白质发色团变化；稳态荧光光谱在三个温度下（298 K、303 K和310 K）定量分析结合常数和猝灭机制；远紫外圆二色谱评估布洛芬诱导的MARK4二级结构变化；分子对接利用AutoDock Vina将布洛芬对接至MARK4晶体结构（PDB ID: 5ES1）的活性位点；分子动力学模拟在GROMACS中使用Amber99sb-ILDN力场进行100 ns三重复模拟，分析复合物的稳定性与动态行为。

三、研究结果

3.1 布洛芬抑制MARK4 ATPase活性

通过孔雀绿ATPase抑制实验，研究人员发现布洛芬对MARK4 ATPase活性呈剂量依赖性抑制。在0.1 μM时抑制率为3.49 ± 1.24%，在100 μM时达到93.85 ± 1.17%。经非线性回归拟合，平均半数抑制浓度（IC₅₀）为4.02 ± 0.33 μM。考虑到酶浓度（5 μM）与IC₅₀值相近，采用Morrison方程分析得到抑制常数（K_i）为2.29 ± 0.31 μM，表明布洛芬在低微摩尔范围内结合MARK4，符合紧结合抑制特征。

3.2 紫外-可见光谱

MARK4与递增浓度布洛芬（2–10 μM）孵育后，吸收光谱显著增强，表明形成了复合物。双倒数作图法得到解离常数（K_d）为1.295 × 10^?5 M，结合常数（K_b）为7.71 × 10⁴ M^?1，吉布斯自由能变（ΔG）为?6.663 kcal/mol，表明结合自发进行。Hill系数接近1，提示近非协同性结合且存在单一主导结合位点。

3.3 稳态荧光光谱

布洛芬引起MARK4内源荧光的浓度依赖性猝灭。Stern-Volmer常数（K_SV）范围为2.90 ± 0.19 × 10⁴ M^?1至4.17 ± 0.25 × 10⁴ M^?1，双分子猝灭速率常数（K_q）远高于扩散控制限，表明为静态猝灭机制。结合常数随温度升高而增大（4.80 ± 0.78 × 10⁴ M^?1至5.76 ± 0.84 × 10⁴ M^?1），结合位点数（n）接近1。van't Hoff分析显示焓变（ΔH°）和熵变（ΔS°）均为正值，表明疏水作用是驱动结合的主要力量。所有温度下的吉布斯自由能变（ΔG°）均为负值（?6.75至?6.38 kcal/mol），证实结合自发进行。

3.4 圆二色光谱

MARK4的远紫外圆二色谱在208 nm附近呈现特征性负峰。加入布洛芬后，负带强度降低，表明MARK4的二级结构发生变化，α-螺旋含量减少，提示配体结合引起了构象调整。

3.5 分子对接

分子对接验证显示，重新对接的抑制剂与原始晶体结构的均方根偏差（RMSD）仅为0.40 ?，证明方法可靠。布洛芬以?7.0 kcal/mol的结合能对接于ATP/抑制剂结合口袋，与Asp196形成两个氢键（距离分别为3.22 ?和3.20 ?），并通过与Tyr134、Gly138、Glu133、Ala195等的范德华接触以及与Leu185、Ala135、Val116、Val70、Ile62等的疏水相互作用稳定结合。值得注意的是，催化关键残基Lys85也参与了疏水相互作用，提示布洛芬可能干扰ATP识别。

3.6 分子动力学模拟

100 ns分子动力学模拟结果显示：apo MARK4的平均RMSD为0.205 ± 0.025 nm，布洛芬复合物略增至0.256 ± 0.042 nm，抑制剂复合物为0.244 ± 0.036 nm，表明两种配体均未破坏蛋白稳定性。残基水平均方根波动（RMSF）显示大部分残基波动低于0.2 nm，柔性增加主要局限于loop区。回转半径（R_g）在apo蛋白（2.062 ± 0.009 nm）、布洛芬复合物（2.058 ± 0.016 nm）和抑制剂复合物（2.051 ± 0.007 nm）间几乎不变，表明整体紧凑性保持。溶剂可及表面积（SASA）从apo的165.748 ± 2.256 nm²轻微增加至布洛芬复合物的167.116 ± 2.143 nm²。二级结构分析显示各组分含量变化极小。氢键分析表明，布洛芬复合物平均形成0.73 ± 0.36个氢键，远低于抑制剂复合物的2.59 ± 0.48个，其中Glu133的氢键占有率最高（22.5%）。主成分分析和自由能景观图显示布洛芬复合物与apo蛋白的全局运动模式相似，但呈现两个能量最小值，表明适度的构象异质性。

四、结论与意义

本研究通过整合实验与计算方法，首次证明布洛芬能够结合MARK4并抑制其ATPase活性。酶学实验确立了剂量依赖性抑制，IC₅₀为4.02 ± 0.33 μM。紫外-可见和荧光光谱证实了自发复合物形成，且疏水作用为主要驱动力。圆二色谱揭示了布洛芬诱导的MARK4二级结构改变。分子对接支持布洛芬占据ATP/抑制剂结合口袋，分子动力学模拟证实复合物在100 ns内保持动态稳定。这些发现表明布洛芬可能作为一种MARK4抑制剂，为进一步的药物发现提供了基础。未来研究应采用正交激酶实验和阳性对照确认布洛芬对MARK4激酶活性的影响，并在细胞模型中验证这些发现，同时评估其对相关酶的选择性，探索基于结构的修饰以提高活性和结合亲和力。