- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
MEK及抑制剂U0126复合物分子动力学模拟.
MEK及抑制剂U0126复合物分子动力学模拟.
摘 要:采用分子对接模拟U0126与MEK结合而形成的复合物,然后对复合物进行10.0 ns的分子动力学模拟。结果表明,活性空腔氨基酸残基间的相互作用和氢键对抑制剂的识别与稳定性都有重要影响。
关键词:U0126;MEK;分子动力学
丝裂原活化蛋白激酶激酶(MAPKK或MEK),是信号转导途径的重要激酶,在信号转导途径中起着承上启下的重要作用。MEK能够整合上游信号并通过磷酸化丝裂原活化蛋白激酶(MAPK)第七和第八亚结构域之间的丝/苏氨酸和酪氨酸残基,从而激活MAPK[1]。
MEK特异性抑制剂U0126是一个新的、可选择、有效的MEK抑制剂。U0126是一个有机物(1,4-二氨基-2,3-氰基-1,4-双[2-氨基苯基硫代]丁二烯),是目前应用最广泛、最有效的一种抑制剂。U0126可以通透细胞,与MEK非竞争结合,抑制其催化活力,从而阻止MEK的磷酸化和激活,进而阻止下游MAPK的活化。而其它一些抑制剂只能与没有活化的MEK结合,从而限制它们的应用。
研究U0126与MEK的作用模式,通过酶晶体结构解析和活性中心结构域及其与抑制剂结合位点的分析,确定与结合相关的关键氨基酸残基,有助于设计和改造出具有更高活性和效果的经济适用的MEK抑制剂。有利于研究病菌的发生、发展和侵入机制,对研究海洋中的病原真菌和寄主植物之间的相互作用具有重要的理论意义[2]。
利用分子对接方法来模拟U0126与MEK的复合物构象,并进行10 ns的分子动力学模拟,通过对模拟结果的分析,得到U0126与MEK之间的作用关系,确定活性空腔内与配体结合的关键氨基酸残基以及形成的氢键对复合物稳定性的影响。
1 材料与方法
1.1 材料 MEK的晶体结构(PDB ID: 3EHQ)来自于蛋白质晶体数据库(Protein Data Bank),手动删除除蛋白质分子以外其他所有结构,对缺失的蛋白质残基用Swiss-PdbViewer程序包[4]补齐。配体小分子U0126的晶体结构来自于小分子数据库(The PubChem Project),并通过Pymol[5]保存为PDB格式,作为分子对接时配体的初始结构[3]。
1.2 分子对接 首先采用AutoDock4.2[6]程序包对MEK和抑制剂U0126进行预处理。对接时受体大分子的格点盒子大小为124×124 ×124,格点间距为0.0375 nm。运用拉马克(Lamarckian)遗传算法,将局部能量搜索与遗传算法相结合,以半经验势函数作为能量打分函数,对小分子构象和位置进行全局搜索,对MEK和抑制剂U0126进行10次的独立对接实验。最后依据成簇分析情况和最低对接能来选取对接复合物作为分子动力学模拟的初始构象[4]。
1.3 MD模拟 本实验运用gromacs程序包选用Opls-aa全原子力场,对复合物在水溶液中进行MD模拟。将复合物置于一个立方体水盒中,质心置于立方体盒子中心。然后在溶质外围加上1.0000 nm的水分子层,采用的水分子为SPCE模型。消除分子间的高能碰撞采用3 000步最陡下降法和2 000步共轭梯度法进行能量最小化的方法,然后进行100 ps的NVT和100 ps的NPT模拟,体系平衡后进行10 nsMD模拟。处理长程静电作用采用周期性边界条件和Particle Mesh Ewald(PME)方法。非键截距为0.9000 nm,使用SHAKE选项用于约束和氢原子连接的键。系统的温度和压力则通过V-rescale和Parrinello-Rahman实现控制。模拟中采用的积分步长为2fs,非键连接表每5步更新一次。每2 ps保存一次轨迹结构用于后期数据处理。 论文代写
2 结果与讨论
2.1 MD轨迹稳定性分析 MD模拟轨迹的均方根偏差(root mean square deviation,RMSD)是衡量体系是否稳定的重要依据,为检测MD轨迹的稳定性,对U0126与MEK复合物的中MEK的骨架原子均方根偏差进行分析。图1,为MEK的骨架原子RMSD随时间变化图。如图所示,体系在开始阶段,RMSD呈上升趋势;经过大约0.8 ns MD模拟后,RMSD趋于稳定,MEK基本达到平衡状态;在大约2.0 ns以后,RMSD值稳定在0.1700 nm,MEK达到平衡状态。此外,在1.8 ns的MD之后,体系势能也趋向最小值而达到稳定,从而为U0126与MEK复合物能够稳定存在提供有力的证据[6]。
图1 MEK酶骨架原子的RMSD变化图
2.2 平均构象的活性空腔分析 通过Gromacs自带程序提取2-5ns模拟轨迹的平均结构,作为U0126与MEK相互作用的代表性构像。然后,使用Ligplo
文档评论(0)