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

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

安装软件

AutoDock Vina下载地址

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

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

结果分析

参考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 人正在系统学习中