网络药理学：12、分子对接之快速实践流程总结：autodock/autodock vina+mgltools实现大分子蛋白（PDB ID已知）和小分子配体的分子对接

一、前言

本篇文章的步骤细节位于：http://t.csdnimg.cn/E8sHZ，即本专栏p11。
（想知道autodock和autodock vina的不同也可见该文。）
这里是快速实践的流程总结版（正文步骤都需要在autodock tools中进行）。
同时本文的第二、三步是autodock和autodock vina通用的，从第四步开始不同。

前置准备：

• 大分子蛋白的PDB格式文件：https://g6altair.sci.hokudai.ac.jp/g6/service/pocasa/
• pymol前置处理大分子蛋白，保存为PUB ID_pymoled.pdb格式文件

二、大分子蛋白准备

• 加载蛋白：选择File/Read Molecule后选择文件加载
• 加氢：选择Edit/Hydrogens/Add进行加氢，弹出两个页面都直接选YES。
• 计算总电荷：选择Edit/Chargs/Compute Gasteiger
• 指定原子类型：选择Edit–Atoms–Assign AD4 type
• 转为pdbqt文件：选择File/save/write PDBQT后直接点击OK

三、小分子配体准备（如果事先准备了小分子的pdbqt文件，可以跳过）

• 重新打开autodock tools软件，或者点击Edit/Delete/Delete All Molecules清空页面。
• 加载配体：选择Ligand/Input/Open（可能跳出弹窗，直接选择yes即可）
• 选择并判断配体的Root：选择Ligand> Torsion Tree > Choose Root”和“Ligand > Torsion Tree > Detect Root。
• 加载更多Root：选择Ligand > Torsion Tree > Show Root Expansion
• 查看可旋转的键：选择Ligand > Torsion Tree > Choose Torsion
• 保存成PDBQT：选择Ligand—Output > Save as PDBQT
• 随后点击Edit/Delete/Delete All Molecules清空软件页面

四、设置对接口袋（GPF文件和config.txt）

如下步骤autodock和autodock vina通用。

• 重新打开autodock tools软件，或者点击Edit/Delete/Delete All Molecules清空页面。
• 加载大分子蛋白：选择Grid > Macromolecule > Open（冒出三个弹窗，全都选Yes）
• 加载小分子配体：选择Grid > Set Map Types > Choose Ligand…或Grid > Set Map Types > Open Ligand…
• 移动小分子配体：……
• 设置并保存对接口袋：选择Grid > Grid Box…设置好口袋大小后点击Grid Options弹窗的File > Close saving current退出。

Autodock

• 保存GPF文件选择Grid > Output > Save GPF...，命名类似7xvy_0472.gpf
• 运行Grid：选择Run > Run AutoGrid
• 点击Edit/Delete/Delete All Molecules清空软件页面。

Autodock vina

• 保存config.txt文件：选择Docking > Output > Vina，得到config.txt，记得查看文件夹中是否真的存在该文件。因为报错只会在bat的终端显示，在mgltools软件页面可能看不到。

五、正式分子对接

Autodock

• 加载大分子蛋白：选择Docking > Macromolecule > Set Rigid Filename…
• 加载小分子配体：选择Docking > Ligand > Choose...或选择Docking > Ligand > Open...
• 设置算法：选择Docking > Search Parameters > Genetic Algorithm，默认设置，点击Accept。
• 导出DPF文件：选择Docking > Docking Parameters…，默认设置，点击Accept。
• 保存算法：选择Docking > Output > Lamarckian GA (4.2)…。文件名类似7xvy_0472.dpf
• 运行Dock：选择Run > Run AutoDock
• 点击Edit/Delete/Delete All Molecules清空软件页面

Autodock vina

可视化操作

• 运行config.txt文件：选择Run > Run Autodock vina，加载刚刚的config.txt，得到小分子配体命名_out.pdbqt文件。
• 点击Edit/Delete/Delete All Molecules清空软件页面。

命令行操作

vina --config config.txt --log 7xvy_0472.txt

执行分子对接并输出日志文件

六、结果分析

Autodock

• 打开.dlg文件：Analyze > Dockings > Open...。（可能有warning弹窗，直接选择ok就好）
• 显示 Receptor大分子蛋白：Analyze > Macromolecule > Open...
• 查看分子对接结果：Analyze > Conformations > Play, ranked by energy...。
• 保存成pdbqt格式文件：点击Write Complex图标，命名类似7xvy_0472.pdbqt

Autodock vina

可视化操作

• 清空页面：Edit > Delete > Delete All Molecules
• 分析结果：Analyze > Dockings > Open AutoDock vina result
• 转为pdbqt文件：选择File > save > write PDBQT后直接点击OK

命令行操作

直接打开日志文件即可。