分子对接如何做啊?

分子对接如何做啊?

简单介绍使用 DSDP 程序进行分子对接的方法,详细内容见

对于要预测的每组复合物,需要建立一个独立的目录,其中存放配体和受体(需按规定命名)。 例如(如下是通过tree DSDP_dataset/ | head 查看的示例目录结构):

DSDP 的输入文件为 pdbqt 格式,可以通过 AutoDockTools 和 Openbabel 等生成。

在项目根目录运行如下命令:

参数解释:

本例使用 1a2b 蛋白。在工作目录(不一定是项目目录)准备受体和配体的 pdbqt 文件,命名为 protein.pdbqt 和 ligand.pdbqt。

运行

注:重对接的盒子信息(沿 x y z 轴的最小值和最大值)需要用户提供。 本例的盒子信息仅适用于本例的 1a2b 蛋白。

参数解释:

通过 --help 参数来查看参数的帮助信息。

以 PET 酶 IsPETase 为例,下载晶体结构,去水加氢、另存为 PET_protein.pdbqt,并与 PET_ligand.pdbqt 存放在 PET 目录下。通过上述命令进行盲对接,结果如图1所示。

可以看到,结合位点的搜索十分准确,且将底物定位到了SHD三联体,这与催化机理完全一致。但酯键的方向理应指向蛋白主链形成的氧负离子穴区域,这里所的结构朝向溶剂,可见盲对接的精度略有不足。

进一步在该位点的基础上进行重对接,如图2所示。(盒子中心为催化三联体,边长为 40 Å)

可以看到,重对接所得的构象中,羰基已经处于正确的朝向,指向Y58和M132的主链酰胺。以上测试证实了 DSDP 盲对接和重对接的可行性。

注意每次使用盲对接前都需要激活DSDP虚拟环境。