《分子对接的原理,方法及应用.docx》由会员分享，可在线阅读，更多相关《分子对接的原理,方法及应用.docx（4页珍藏版）》请在三一文库上搜索。
1、分子对接的原理,方法及应用分子对接的原理 , 方法及应用(PPT 里弄一些分子对接的照片, 照片素材文件里有 )分子对接就是将已知三维结构数据库中的分子逐一放在靶标分子的活性位点处。通过不断优化受体化合物的位置、 构象、分子内部可旋转键的二面角与受体的氨基酸残基侧链与骨架 , 寻找受体小分子化合物与靶标大分子作用的最佳构象 , 并预测其结合模式、亲与力与通过打分函数挑选出接近天然构象的与受体亲与力最佳的配体的一种理论模拟分子间作用的方法。通过研究配体小分子与受体生物大分子的相互作用 , 预测其亲与力 , 实现基于结构的药物设计的一种重要方法。原理 :按照受体与配体的形状互补 , 性质互补原则
2、, 对于相关的受体按其三维结构在小分子数据库直接搜索可能的配体 , 并将它放置在受体的活性位点处 , 寻找其合理的放置取向与构象 , 使得配体与受体形状互补 , 性质互补为最佳匹配( 配体与受体结合时 , 彼此存在静电相互作用 , 氢键相互作用 , 范德华相互作用与疏水相互作用 , 配体与受体结合必须满足互相匹配原则 , 即配体与受体几何形状互补匹配 , 静电相互作用互补匹配 , 氢键相互作用互补匹配 , 疏水相互作用互补匹配 )目的 :找到底物分子与受体分子的最佳结合位置问题 :如何找到最佳的结合位置以及如何评价对接分子之间的结合强度方法 :1、 首先建立大量化合物的三维结构数据库2、 将库
3、中的分子逐一与靶分子进行“对接”3、 通过不断优化小分子化合物的位置以及分子内部柔性键的二面角子化合物与靶标大分子作用的最佳构象 , 计算其相互作用及结合能, 寻找小分4、 在库中所有分子均完成了对接计算之后 , 即可从中找出与靶标分子结合的最佳分子应用 :1) 直接揭示药物分子与靶点之间的相互作用方式2) 预测小分子与靶点蛋白结合时的构象3) 基于分子对接方法对化合物数据库进行虚拟筛选 , 用于先导化合物的发现分子对接的原理,方法及应用4) 预测化合物的亲与力及活性 , 用于先导化合物的优化分子对接思想来源于“锁与钥匙” , 但又比“锁与钥匙”复杂得多 , 表现在以下方面 :1) 药物分子与
4、靶酶分子就是柔性的 , 这样就要求在对接过程中要相互适应以达到最佳匹配2) 分子对接不仅要满足空间形状匹配 , 还要满足能量匹配 , 底物分子与靶酶分子能否结合以及结合的强度最终取决于形成此复合物进程的结合自由能。分子对接的种类 :1) 刚体对接 : 对接过程中 , 研究体系的构象不发生变化。 适合考察比较大的体系 , 如蛋白质与蛋白质间以及蛋白质与核酸间的对接2) 半柔性对接 : 对接过程中 , 研究体系尤其就是配体的构象允许在一定的范围内变化。适合处理大分子与小分子间对接 , 对接过程中 , 小分子的构象一般就是可以变化的 , 但大分就是刚性的3) 柔性对接 : 对接过程中 , 研究体系的
5、构象就是基本上可以自由变化的。 一般用于精确考虑分子间的识别情况 , 由于计算过程中体系的构象可以变化 , 所以计算破费最大分子对接中的问题 :如何找到最佳的结合位置。这牵涉到优化的问题 , 底物分子与受体分子都就是可以自由转动与平动的 , 同时两个分子自身的构象也存在变化 , 因此它们之间可能的结合方式就是非常复杂的 , 常用的有遗传算法、模拟退火以及禁忌搜索等如何确定对接分子间的结合强度。这涉及到底物分子与受体分子间结合能力的预测, 牵涉到结合自由能的计算。对接的基本类型 :1) 整体分子对接法 : 运用一种特定搜索算法考察配体分子在受体结合部位的能量 , 并找出对应给定评分函数的最优结合
6、方式2) 基于片断对接法 : 配体分子被视为若干片断结构的集合 , 先将其中一个或几个基本片断放入结合口袋 , 然后在活性部位构建分子的其余部分 , 最终得到理论上最优的结合方式几种有代表性的分子对接方法4、1 DOCKDOCK 就是 Kuntz 研究小组发展的分子对接程序 , 可能就是目前应用最为广泛的分子对接程序之一、 它能自动地模拟配体分子在受体活性位点的作用情况 , 并把理论预测最佳的方式记录下来。 而且该方法能够对配体的三维结构数据库进分子对接的原理,方法及应用行自动搜索 , 因此被广泛应用于基于受体结构的数据库搜索的药物设计中 , 并取得了巨大的成功。用 DOCK进行药物设计以及数
7、据库的搜索基本上可以分为下面几个步骤 : 配体与受体相互作用位点的确定 , 评分系统的生成 ,DOCK 计算及 DOCK 结果的处理与分析。活性位点的确定与表达就是 DOCK最重要的特点之一。活性位点特征的确定对于 DOCK研究就是非常重要的 , 因为配体分子与受体相互作用过程的模拟主要就就是参考几何位点的几何特征进行的。 在 DOCK中, 活性位点的确定通过 sphgen 程序来完成。 DOCK软件包中 sphgen 程序生成受体表面所有的凹陷的负像 , 并对这些负像进行聚类分析。在 DOCK程序中 , 表面点采用了 Richards 提出的模型。 在这些表面点的基础上 , 采用 sphge
8、n 程序生成了负像 , 它实际上由一些与分子表面点相切的圆球叠加而成。在生成负像的基础上 , 就可以进行配体分子与活性口袋之间的匹配。 在这里 , 配体也采用一组球集来表示 , 与负像不同的就是 , 配体所用的球集表示配体所占的空间区域。如果配体分子能与活性口袋形成比较好的匹配 , 那么配体的球集一定能与活性口袋中的负像形成好的叠合。 配体分子与负像之间的匹配原则就是基于配体与受体之间球集的内坐标的比较。按照匹配原则得到了配体与受体之间的匹配情况之后、 就要通过合理的得分函数来选择最优的结果。 DOCK提供了多种得分函数来评价配体与受体之间的结合情况 , 包括原子接触得分以及能量得分。 DOC
9、K 进行分子对接时 , 配体分子可以就是柔性的。对于柔性的分子 , 其键长与键角保持不变 , 但可旋转二面角就是可以发生变化的。在 DOCK中, 柔性分子的构象变化通过下面的操作实现 : 首先就是刚性片断的确定 , 然后就是构象搜索。构象搜索采用两种方法 : 一种就是锚优先搜索(anchor-first search),第二种方法就是同时搜索(simultaneous search)、4、2 AUTODOCKAUTODOCK 就是 Scripps 的 Olson 科研小组开发的分子对接软件包 , 最新的版本为 3、05,AUTODOCK采用模拟退火与遗传算法来寻找受体与配体最佳的结合位置,用半
10、经验的自由能计算方法来评价受体与配体之间的匹配情况。 在 AUTODOCK中,配体与受体之间结合能力采用能量匹配来评价。在 1、0 与 2、0 版本中 , 能量匹配得分采用简单的基于 AMBER力场的非键相互作用能。 非键相互作用来自于三部分的贡献 : 范得华相互作用 , 氢键相互作用 , 以及静电相互作用。在 3、0 版中 ,AUTODOCK提供了半经验的自由能计算方法来评价配体与受体之间的能量匹配 在最早的 AUTODOCK版本中 , 作者采用了模拟退火来优化配体与受体之间的结合。在 3、0 版本中 ,Morris 等发展了一种改良的遗传算法 , 即拉马克遗传算法 (LGA)。测试结果表明
11、 ,LGA 比传统的遗传算法比模拟退火具有更高的效率。在 LGA 方法中 , 作者把遗传算法与局部搜索 (local search) 结合在一起 , 遗传算法用于全局搜索 , 而局部搜索用于能量优化。在 AUTODOCK中, 局部搜索方法就是自适应的 , 它可以根据当前的能量调节步长大小。 LGA 算法引入了拉马克的遗传理论 ,LGA 最大的特点就就是通过进化映射 ( developmental mapping ) 把基因型转化为表现型而实现局部搜索与遗传算法的结合。 基因型空间通过遗传算子突变与交叉来定义 ; 而表现型则通过问题的解来定义 , 这里表示体系的能量得分。分子对接的原理,方法及应
12、用4、3 FIexXFIexX就是德国国家信息技术研究中心生物信息学算法与科学计算研究室的 Matthias Rarey 等发展的分子对接方法 , 现在己经作为 SYBYL分子模拟软件包中的一个模块实现了商业化。 FIexX 中结合了多种药物设计的方法进行配体与受体之间的对接。在 FIexX 中 , 配体与受体之间结合情况的评价采用了类似Bbhm提出的基于半经验方程的自由能评价方法。 在 FlexX 中 , 分子对接的流程主要分为下面的步骤 :(1) 核心结构确定 对接的第一步为核心片断 (base fragment) 的选择。核心片断就是指能对配体与受体之间相互作用起决定作用的基团 , 而且
13、核心片断的构象要尽量少一些。核心基团的正确选择对分子对接的计算结果有非常重要的影响。因为如果核心基团与受体之间不存在明显优势的相互作用时 , 则很难正确预测正确的结合模式。随着核心基团的增加 , 则核心基团与受体之间的相互作用也会相应增强 , 那么结合模式准确预测的机会就会大大增加。因此在选择核心基团时 , 核心基团包含的基团要尽量多一些。而且核心基团的构象数要尽量少一些。当核心基团选定以后 , 就可以把配体分子划分为多个片断。(2)(2)核心结构的放置当选择好了核心基团以后, 就要把核心基团放置在活性位点的正确部位。在放置核心基团的时候 ,FlexX 采用了一种形态聚类算法 (pose cl
14、ustering algorithm) 算法 ), 在这个算法中 , 一个核心基团可以瞧作为一个具有明确相互作用点的刚性物体 , 而受体分子的活性口袋也可以瞧作为一个具有明确相互作用点的刚性物体。 把核心基团放置在活性口袋中的过程就相当于把配体中的二个相互作用点叠合在活性位点的三个相互作用点上 ( 假设这三个点不共线 ) 。在匹配中 , 两个三角形三个顶点所具有的相互作用特征应该就是符合的 , 同时三角形对应边长的差别应该在一定的范围内。 放置核心结构的第一步就就是找出所有相匹配的这些三角形 , 而且对配体的位置进行坐标转化。当所有的转化完成以后 , 检查配体就是否与受体产生了碰撞 , 去掉一
15、些不合理的核心结构取向。对于得到这些核心结构的可能位置 , 通过核心结构空间的位置均方根位移 (r 、m、s、D ) 进行聚类分析。对应那些 r 、m、s、D 值小于一定阀值的空间位置进行归并 , 仅仅保留那些相差较大的空间位置。HIV 蛋白酶HIV 蛋白酶就是 HIV 产生的 , 这类蛋白酶专门负责病毒及其酶的前体清除工作 , 正就是在这种酶的作用下产生了成熟的、易感染的病毒。 Hiv 病毒复制周期正常运转与病毒毒粒成熟至关重要 , 就是病毒复制必需的酶 , 可以作为抗 HIV 的药物靶点。在寻找治疗艾滋病的方法中 ,HIV 蛋白酶抑制剂以其花费低廉、研制周期短、与目前生产水平接近以及较 HIV 疫苗安全等特点 , 在抗病毒的化学治疗中 , 越来越引起人们的注意 , 它就是近年来治疗艾滋病的方法中发展较快的方法之一。锌美洲铜斑蛇的毒液中发现的HIV 蛋白酶抑制剂含有Zn2+。已有的工作表明 ,Zn 2+ 直接与 2 个羧基络合 , 或与适当距离上带相反电荷的部分键合而发挥抑制作用。明明 明