mek与抑制剂u0126复合物的分子动力学模拟

《mek与抑制剂u0126复合物的分子动力学模拟》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、MEK与抑制剂U0126复合物的分子动力学模拟采用分子对接模拟U0126与MEK结合而形成的复合物，然后对复合物进行10.0ns的分子动力学模拟。结果表明，活性空腔氨基酸残基间的相互作用和氢键对抑制剂的识别与稳定性都有重要影响。关键词：U0126；MEK；分子动力学　　丝裂原活化蛋白激酶激酶(MAPKK或MEK)，是信号转导途径的重要激酶，在信号转导途径中起着承上启下的重要作用。MEK能够整合上游信号并通过磷酸化丝裂原活化蛋白激酶(MAPK)第七和第八亚结构域之间的丝/苏氨酸和酪氨酸残基，从而激活MAPK[1]。

2、　　MEK特异性抑制剂U0126是一个新的、可选择、有效的MEK抑制剂。U0126是一个有机物(1，4-二氨基-2，3-氰基-1，4-双[2-氨基苯基硫代]丁二烯)，是目前应用最广泛、最有效的一种抑制剂。U0126可以通透细胞，与MEK非竞争结合，抑制其催化活力，从而阻止MEK的磷酸化和激活，进而阻止下游MAPK的活化。而其它一些抑制剂只能与没有活化的MEK结合，从而限制它们的应用。　　研究U0126与MEK的作用模式，通过酶晶体结构解析和活性中心结构域及其与抑制剂结合位点的分析，确定与结合相关的关键氨基酸残基，

3、有助于设计和改造出具有更高活性和效果的经济适用的MEK抑制剂。有利于研究病菌的发生、发展和侵入机制，对研究海洋中的病原真菌和寄主植物之间的相互作用具有重要的理论意义[2]。　　利用分子对接方法来模拟U0126与MEK的复合物构象，并进行10ns的分子动力学模拟，通过对模拟结果的分析，得到U0126与MEK之间的作用关系，确定活性空腔内与配体结合的关键氨基酸残基以及形成的氢键对复合物稳定性的影响。　　1　材料与方法　　1.1　材料　MEK的晶体结构(PDBID:3EHQ)来自于蛋白质晶体数据库（ProteinDat

4、aBank），手动删除除蛋白质分子以外其他所有结构，对缺失的蛋白质残基用SProject)，并通过Pymol[5]保存为PDB格式，作为分子对接时配体的初始结构[3]。　　1.2　分子对接　首先采用AutoDock4.2[6]程序包对MEK和抑制剂U0126进行预处理。对接时受体大分子的格点盒子大小为124×124×124，格点间距为0.0375nm。运用拉马克（Lamarckian）遗传算法，将局部能量搜索与遗传算法相结合,以半经验势函数作为能量打分函数，对小分子构象和位置进行全局搜索，对MEK和抑制剂U012

5、6进行10次的独立对接实验。最后依据成簇分析情况和最低对接能来选取对接复合物作为分子动力学模拟的初始构象[4]。　　1.3　MD模拟　本实验运用gromacs程序包选用Opls-aa全原子力场，对复合物在水溶液中进行MD模拟。将复合物置于一个立方体水盒中，质心置于立方体盒子中心。然后在溶质外围加上1.0000nm的水分子层,采用的水分子为SPCE模型。消除分子间的高能碰撞采用3000步最陡下降法和2000步共轭梯度法进行能量最小化的方法，然后进行100ps的NVT和100ps的NPT模拟，体系平衡后进行10nsM

6、D模拟。处理长程静电作用采用周期性边界条件和ParticleMeshE,使用SHAKE选项用于约束和氢原子连接的键。系统的温度和压力则通过V-rescale和Parrinello-Rahman实现控制。模拟中采用的积分步长为2fs,非键连接表每5步更新一次。每2ps保存一次轨迹结构用于后期数据处理。　　2　结果与讨论　　2.1　MD轨迹稳定性分析　MD模拟轨迹的均方根偏差(rootmeansquaredeviation，RMSD)是衡量体系是否稳定的重要依据，为检测MD轨迹的稳定性,对U0126与MEK复合物的中

7、MEK的骨架原子均方根偏差进行分析。图1，为MEK的骨架原子RMSD随时间变化图。如图所示，体系在开始阶段，RMSD呈上升趋势；经过大约0.8nsMD模拟后，RMSD趋于稳定，MEK基本达到平衡状态；在大约2.0ns以后，RMSD值稳定在0.1700nm，MEK达到平衡状态。此外，在1.8ns的MD之后，体系势能也趋向最小值而达到稳定，从而为U0126与MEK复合物能够稳定存在提供有力的证据[6]。　　　　图1MEK酶骨架原子的RMSD变化图　　2.2　平均构象的活性空腔分析　通过Gromacs自带程序提取2-5

8、ns模拟轨迹的平均结构，作为U0126与MEK相互作用的代表性构像。然后，使用Ligplot程序画出U0126与MEK相互作用的平面结构图。MEK活性位点处残基与底物U0126的相互作用如图2所示，在MEK的活性空腔中，Phe129、Met143、Ile141、Leu115、Leu215、Ile216和Met219七个氨基酸残基的疏水侧链都能和U0126形成疏水作用，使得