《AutoDock_ADT教程_AutoDockTools_分子对接_中文》由会员分享,可在线阅读,更多相关《AutoDock_ADT教程_AutoDockTools_分子对接_中文(38页珍藏版)》请在金锄头文库上搜索。
1、1分子对接使用 AutoDock 和 AutoDock Tools一、分子对接简介及软件介绍二、对接准备及对接操作三、结果分析一、分子对接简介及软件介绍1分子对接理论基础所谓分子对接就是两个或多个分子之间通过几何匹配和能量匹配而相互识别的过程。分子对接在酶学研究以及药物设计中具有十分重要的意义。在酶激活剂、酶抑制剂与酶相互作用以及药物分子产生药理反应的过程中,小分子(通常意义上的Ligand)与靶酶(通常意义上的Receptor )相互互结合,首先就需要两个分子充分接近,采取合适的取向,使两者在必要的部位相互契合,发生相互作用,继而通过适当的构象调整,得到一个稳定的复合物构象。通过分子对接确定
2、复合物中两个分子正确的相对位置和取向,研究两个分子的构象,特别是底物构象在形成复合物过程中的变化,是确定酶激活剂、抑制剂作用机制以及药物作用机制,设计新药的基础。分子对接计算是把配体分子放在受体活性位点的位置,然后按照几何互补、能量互补化学环境互补的原则来实时评价配体与受体相互作用的好坏,并找到两个分子之间最佳的结合模式。分子对接最初思想起源于Fisher E.的“锁和钥匙模型” (图 ) ,认为“锁”和“钥匙”的相识别的首要条件是他们在空间形状上要互相匹配。然而,配体和受体分子之间的识别要比“锁和钥匙”模型复杂的多。首先,配体和受体分子的构象是变化的,而不是刚性的,配体和受体在对接过程中互相
3、适应对方,从而达到更完美的匹配。其次,分子对接不但要满足空间形状的匹配,还要满足能量的匹配。配体和受体之间的通过底物分子与靶酶分子能否结合以及结合的强度最终是由形成此复合物过程的结合自由能变化G bind所决定的。互补性(complementarity)和预组坦织( pre-organization)是决定分子对接过程的两个重要原则,前者决定识别过程的选择性,而后者决定识别过理的结合能力。互补性包括空间结构的互补性和电学性质的互补性。1958年Koshland提出了分子识别过程中的诱导契合(induced fit)概念,指出配体与受体相互结合时,受体将采取一个能同底物达到最佳结合的构象(图1)
4、。而受体与配体分子在识别之前将受体中容纳配体的环境组织2的越好,其溶剂化能力越低,则它们的识别效果越佳,形成的复合物也就越稳定。图 1:“锁和钥模型(a) ”和“诱导契合模型(b)”示意图2分子对接常用方法分子对接方法根据不同的简化程度可以大致分为以下三类:(1)刚性对接;(2)半柔性对接;(3)柔性对接。刚性对接指在对接的过程中,研究体系的构象不发生变化;半柔性对接是指在对接过程中,研究体系尤其是配体的构象允许在一定的范围内变化;而柔性对接则是指在对接过程中,研究体系的构象基本上是可以自由变化的。当然,这只是一种简单的分类方法,而在很多分子对接程序中,实际上采取了多种处理方法。在这些分子对接
5、方法中,刚性对接适合考察比较大的体系,比如蛋白质和蛋白质以及蛋白质和核酸之间的相互作用,他计算较为简单,原理也相对简单,主要是考虑构象之间的契合程度。半柔性对接适合于处理小分子和大分子之间的对接。在对接过程中,小分子的构象一般是可以变化的,但大分子则是刚性的。由于小分子相对较小,因此在一定程度考察柔性的基础上,还可以保持较高的计算效率,在药物设计尤其在基于分子数据库的虚拟筛选过程中,一般采用半柔性的分子对接方法。柔性对接方法一般用于精确考察分子之间的识别情况,由于在计算过程中体系的构象是可以变化的,因此柔性对接在提高了对接准确性的同时却需要耗费较长的计算时间。上述各种几何优化方法可以获得分子对
6、应与初始态的优势构象,但实际上这样的构象可以有很多,一般认为自由能最小的构象存在的概率高,全局极小可能是比较重要的构象。同时,配体和受体活性部位结合时的构象不一定是全局极小构象,而可能是某一低能构象(药效构象) 。不管是寻找分子的全局极小构象还是药效构象,均要使用构象搜索方法。分子对接的目的是找到底物分子和受体分子的最佳结合位置。因此,分子对接面临的最重要的问题就是如何找到最佳的结合位置以及如何评价对接分子之间的结合强度。当然,这两个问题也是3相互关联的。如何找到最佳的结合位置就要牵涉到构象搜索方法。常用的构象搜索方法有系统搜索法和非系统搜索法。系统搜索法通过系统的改变每一个扭转角产生所有可能
7、的构象,从中挑选出能量较低的构象,但计算量非常大。所以通常使用非系统搜索法来寻找能量较低构象,常用方法有:(1)分子动力学方法(Molecular dynamics,MD ) ;(2)随机搜索(Random search) ;(3)遗传算法(Genetic algorithm,GA ) ;(4 )距离几何算法(Distance geometry,DG)等。随机搜索又包含:(1)完全随机算法;(2)蒙特卡罗法(Monte Carlo,MC) ;(3)模拟退火法(Simulated annealing,SA) 。分子对接的方法很多,下表(表 1)列出了针对不同对接体系的常用对接方法:表 1:常用的
8、分子对接方法对接类型 对接方法类型 对接方法Conformational 构象Fragmentation 片段生长Systematic系统方法Database 数据路Monte Carlo (MC) 蒙特卡罗Genetic algorithm (GA) 遗传算法Random/stochastic随机方法Tabu Search 禁忌搜索Molecular dynamics (MD) 分子动力学Flexible-Ligand Docking柔性配体对接Simulation methods模拟方法 Energy minimization 能量最小化Molecular dynamics (MD) 分子
9、模拟Monte Carlo (MC) 蒙特卡罗Rotamer libraries 旋转异构体库Protein-ensemble grids 蛋白集合栅格Flexible-Protein Docking柔性蛋白对接Soft-receptor modeling 软受体建模3分子对接软件第一个分子对接程序是 UCSF Kuntz 小组于 1982 年开发的 DOCK,早期的版本以刚性对接为主,从 4.0 版开始考虑配体的柔性。 像这样的半柔性(刚性受体柔性配体:rigid receptor-flexible ligand docking)对接程序还有 AutoDock、FlexX 等。同时考虑配体和
10、受体柔性的对接程序主要有 FlexiDock,它采用遗传算法来对配体和受体的结合构像进行优化,在初始条件较好的情况下,FlexiDock 可以比较精确的确定配体和受体的结合状态,但计算时间较长。下表为一些有代表性的分子对接软件(表 2)及软件的被引次数对比(图 2) 。表 2:有代表性的分子对接软件名称 优化方法 评价函数 速度 Protein-protein (peptide) docking Protein-ligand dockingDOCK 片断生长 分子力场、表面匹配得分、 化学环境匹配得分 快 4AutoDock 遗传算法 半经验自由能评价函数 一般 ICM-Docking 随机全
11、局优化 半经验自由能评价函数 快 GOLD 遗传算法 半经验自由能评价函数 快 FlexX 片断生长 半经验自由能评价函数 快 Affinity 蒙特卡罗/分子力 学/分子动力学 分子力场 慢 ZDock&RDock 几何匹配/分子动 力学 CAPRI* /分子力场 慢 FlexiDock 遗传算法 分子力场 慢 eHiTS 系统搜索 半经验自由能评价函数 快 Hex 几何匹配 CAPRI* 快 注:CAPRI ,Critical Assessment of PRediction of Interactions,相互作用预测临界评价图 2:分子对接软件 最常用的分子对接软件被引次数统计(ISI
12、 Web of Sciense,2005)4AutoDock 以及 ADT(AutoDock Tools )4.1.AutoDockAutoDock 是 The Scripps Research Institute 的 OIson 科研小组使用 C 语言开发的分子对接软件包,目前最新的版本为 4.01。AutoDock 其实是一个软件包,其中主要包含 AutoGrid 和AutoDock 两个程序。其中 AutoGrid 主要负责格点中相关能量的计算,而 AutoDock 则负责构象搜索及评价。5AutoDock 在早期版本中使用的是模拟退火算法(Simulated Annealing Alg
13、orithm)来寻找配体与受体最佳的结合位置状态,而从 3.0 版本开始使用一种改良的遗传算法,即拉马克遗传算法(Lamarckian Genetic Algorithm,LGA)。测试结果表明, LGA 比传统的遗传算法和模拟退火具有更高的效率。在 LGA 方法中,作者把遗传算法和局部搜索(Local search)结合在一起,遗传算法用于全局搜索,而局部搜索用于能量优化。LGA 算法引入了拉马克的遗传理论,这个操作过程可下图(图 3) 。图 3:LGA 算法操作过程图同时在AutoDock中配体和受体之间结合能力采用能量匹配来评价。在1.0和2.0 版本中,能量匹配得分采用简单的基于AMB
14、ER力场的非键相互作用能。非键相互作用来自于三部分的贡献:范德华相互作用,氢键相互作用以及静电相互作用。而在3.0 之后的版本中AutoDock提供了半经验的自由能计算方法来评价受体和配体之间的能量匹配。为了加快计算速度,AutoDock 采用格点对接的方法,但与 DOCK 中格点对接的处理方法有明显的区别。DOCK 中,格点上保存的不是能量,而是仅与受体有关的特征量。而在 AutoDock 中,格点上保存的是探针原子和受体之间的相互作用能。对于范德华相互作用的计算, 每个格点上保存的范德华能量的值的数目与要对接的配体上的原子类型(表 3)的数目一样。如果一个配件中含有 C、O 和 H 三种原
15、子类型,那么在每个格点上就需要用三个探针原子来计算探针原子与受体之间的范德华相互作用值。当配体和受体进行分子对接时,配体中某个原子和受体之间的相互作用能通过周围 8 个格点上的这种原子类型为探针的格点值用内插法得到。6表 3:AutoDock4 中的原子类型(*为默认在 gpf 中存在的原子类型)Atom TypeH Non H-bonding HydrogenHD* Donor 1 H-bond HydrogenHS Donor S Spherical HydrogenC* Non H-bonding Aliphatic CarbonA* Non H-bonding Aromatic Car
16、bonN* Non H-bonding NitrogenNA* Acceptor 1 H-bond NitrogenNS Acceptor S Spherical NitrogenOA* Acceptor 2 H-bonds OxygenOS Acceptor S Spherical OxygenF Non H-bonding FluorineMg Non H-bonding MagnesiumMG Non H-bonding MagnesiumP Non H-bonding PhosphorusSA* Acceptor 2 H-bonds SulphurS Non H-bonding SulphurCl Non H-bonding ChlorineCL Non H-bonding ChlorineCa Non H-bonding CalciumCA Non H-bonding CalciumMn Non H-bonding ManganeseMN Non H-bonding Mangane