5fd21b7250c455126.jpg_fo742.jpg

何对接计算都有两个相互矛盾的方面需要平衡:在尽可能精准的计算与合理(有限)的计算资源之间达到一个平衡。理想的步骤是通过搜索整个系统可能的自由度,在底物和目标蛋白的结合能中找到全局能量极小值。然而这样的工作只能在大型的工作站上实现,并且耗费和结构生物学家进行晶体结构修饰相当的时间。为了解决这一问题,很多对接软件简化了对接的步骤。Autodock 通过两种方法的结合使用解决了以上的问题:快速的基于格点能量的计算方法(rapid grid-based energy evaluation)和有效的扭转自由度搜索方法(efficient search of torisional freedom)。
AutoDock软件由 AutoGrid 和 AutoDock两个程序组成。其中 AutoGrid 主要负责格点中相关能量的计算,而 AutoDock 则负责构象搜索及评价。

AutoDock 在早期版本中使用的是模拟退火算法(Simulated Annealing Algorithm)来寻找配体与受体最佳的结合位置状态,而从 3.0 版本开始使用一种改良的遗传算法,即拉马克遗传算法(Lamarckian Genetic Algorithm,LGA)。测试结果表明,LGA 比传统的遗传算法和模拟退火具有更高的效率。在 LGA 方法中,作者把遗传算法和局部搜索(Local search)结合在一起,遗传算法用于全局搜索,而局部搜索用于能量优化。LGA 算法引入了拉马克的遗传理论,这个操作过程如下图
52_648_6176f990038a8a3.gif
同时在AutoDock中配体和受体之间结合能力采用能量匹配来评价。 在1.0和2.0版本中,能量匹配得分采用简单的基于AMBER力场的非键相互作用能。非键相互作用来自于三部分的贡献:范德华相互作用,氢键相互作用以 及静电相互作用。而在3.0之后的版本中AutoDock提供了半经验的自由能计算方法来评价受体和配体之间的能量匹配。
为了加快计算速度,AutoDock 采用格点对接的方法,但与 DOCK中格点对接的处理方法有明显的区别。DOCK 中,格点上保存的不是能量,而是仅与受体有关的特征量。而在 AutoDock 中,格点上保存的是探针原子和受体之间的相互作用能。
对于范德华相互作用的计算, 每个格点上保存的范德华能量的值的数目与要对接的配体上的原子类型(表 3)的数目一样。如果一个配件中含有 C、O 和 H 三种原子类型,那么在每个格点上就需要用三个探针原子来计算探针原子与受体之间的范德华相互作用值。当配体和受体进行分子对接时,配体中某个原子和受体之 间的相互作用能通过周围 8 个格点上的这种原子类型为探针的格点值用内插法得到。

发表回复

您的电子邮箱地址不会被公开。 必填项已用*标注

Captcha Code