swissdock 和autodock的分子模拟

2020-10-19  本文已影响0人  历经沧海一梦曲不终

方式1:SwissDock 

www.swissdock.ch/docking/view/swissdockd_X0ivbA_MD94LBUQL0XHMQPD2ZX9

根据提示上传大分子和小分子即可(有时网站不稳定,上传不了或者搜索不了都为正常情况)

1.      分析结果全部下载

2.      Open(用 chimerae打开即可)

方式2:autodock(该软件在python语言中运行)

准备工作:将文件

粘贴至新的文件夹中,双击adt.bat即可打开改路径下的文件。

1.    蛋白质的处理:(也可以用pymol处理)

①加H

②计算电荷

③蛋白质类型(atoms assign AD4)

Save文件为.pdbqt

2.     小分子处理:

①Ligant(注意更改识别文件)

②ligant-input-open(Mol2文件,该文件可直接从ZINC database中下载Mol2格式,也可以自己画图再用openbabel进行转换)

③Torsion Tree:detect-choose Torsions

Ligant-save

3.      Grid-macmolecule-open

Grid-set map type-open ligant

4.      设置对接范围 grid-grid box

根据文献确定活性位点,若没有参考资料则可将大分子都选中(切记不能框住小分子)

File-saving current

Grid-output(.gpf file)切记文件名加后缀.gpf

5.      Run autogrid(自动识别.gpf)

设置:输出文件在同一路径下的文件夹,后缀改名为.glg(该过程运行时会产生多个map文件,运行时间约10min)

Edit-delete(删除软件中显示的现有文件)

6.     分子对接

Docking- macmolecule-rigid(根据自己的需要选择)

Docking-ligant

运行算法:Docking-search parameters-genetic(默认,根据需求进行选择,其中local是最快的计算方式,但是需要重复多次进行运算,平均一次5min左右);

Output(什么算法导出什么文件),加后缀.dpf

运行run autodock-accept(该步骤生成.dlg的文件)

7.      查看对接结果:

Analyze-docking-open(打开.dlg文件)

加protein:Analyze-macmolecule

分析:Analyze-conformations-play,rankand energy

Show information / build H-binds

导出:write complex(加后缀.pdbqt)

用openbabel将格式转换为.pdb文件。

8.       Pymol软件可视化

打开上述.pdb文件

显示氢键:[S]显示sequence,选中小分子和蛋白位点(上述预测的结果)

Show sticks

Action-find-excluding main chain(显示氢键)

显示距离:atoms(and joints)/工具栏中wizard可计算与显示距离

更改背景为白色:工具栏中background-white

显示残基:selection resides条件下操作-label

and joints条件下control+残基可改变标记位置

右键Ray1200更改分辨率

Export

最终显示如下:

图片参考来源: https://www.wecomput.com/wp-content/uploads/2015/10/Docking.png

上一篇下一篇

猜你喜欢

热点阅读