【失败】Autodock Vina批量分子对接实验记录

• 安装软件
• 准备配体分子和受体蛋白pdb
• • 受体蛋白pbd结构获取：
  • 配体分子pbd结构获取：
• 预处理配体分子和受体蛋白
• 对接盒子(Grid Box)的设置
• • 结果分析

安装软件

AutoDock Vina下载地址

右键“此电脑” > 高级系统设置 > 环境变量 > 双击系统变量的path > 将vina和MGLTools的地址写入


设置工作空间：打开AutoDockTools

open babel下载地址

预处理配体分子和受体蛋白

用pymol打开蛋白受体，执行两个命令remove solvent(去除溶剂)和remove organic（去除小分子）即可。
安装pymol

操作蛋白分子加氢，参考link
操作小分子加氢 ，但是没有成功（如图）。而且有博主说obabel直接转换为pdbqt不能使用。所以就换了Meeko方式转换。先将mk_prepare_ligand.py放在工作文件夹中，再进行run。（此处自动为小分子加氢)

所以我参考利用pyMOL插件GetBox Plugin确定分子对接口袋 下载GetBox计算受体蛋白的对接位置。但是可能是因为我是直接在uniprot数据库搜索蛋白编号得到的蛋白，无法计算出对接位置。
然后我就把蛋白分子导入mgltools用grid box，手动将蛋白分子全包围了。参考分子对接教程 | (6) AutoDock对接操作与对接结果解读

结果分析

参考win10中使用AutoDock Vina进行批量对接的批处理程序进行批量处理分子对接，但是这样全包围分子对接，一个配体的运行时间几乎要2小时，代价太高。所以正在想其他方法（或者外包，毕竟我的专业和这个完全打不着）。
conf.txt

dock.bat

用这两个文件跑了很长时间得到一个结果（如下）

搜索后得到结果说明：

查看打分列表，了解各个构象的打分分布情况，比较各个化合物的打分差异。
RMSD：表示各构象相对于最佳构象（Pose 1）的重原子（除氢原子外的其他原子）空间位置差异程度。
RMSD_ub：RMSD upper bound，即构象间原子严格对应的RMSD（忽略对称性），常用于非对称分子；
RMSD_lb：RMSD lower bound，即构象间原子最佳RMSD（考虑对称性），常用于对称分子。
注意：在普通分析中，大可不必考虑该RMSD。

对接软件的打分函数给出来的一个参数。这个参数是希望用来表征配体与蛋白质口袋的结合强弱。可以认为一定程度上是为了类比delta G。
rmsd:root-mean-square deviation, 在vina里用于表明各个对接位置之间的差异。因为是以第一个位置作为参照，所以第一个位置的rmsd是0.
通常选用一个对接软件，及确定参数设置，需要做一步方法检验。常用的方法就是你对一些晶体结构中的配体用你的对接软件和参数方法做
对接，如果你对接出来的位置跟晶体结构中配体的位置的rmsd在2A范围内，基本认为你的对接软件及方法有一定的可靠性。希望有帮助
怎么确定这个affinity多大合适，完全与取决于你的体系和你想实现的目标。不同的体系都是不一样。从delta G的角度看，只要delta G小于零，
反应都是自发进行。但是作为一个药物设计，当然不能仅仅能结合就够。通常来说，绝对值越高越高。当然最后还是那句话，出发前，先要
知道自己的目的和参照，才能知道自己是否达到目标。

文章知识点与官方知识档案匹配，可进一步学习相关知识Python入门技能树预备知识常用开发工具210344 人正在系统学习中

来源：testuser100