蛋白-蛋白对接一直是分子模拟中非常重要同时非常难解决的问题,相较于小分子与蛋白之间的对接,蛋白-蛋白对接如今较不成熟,在蛋白-蛋白对接之前,最好能够搜集更多的文献进行支持,让模拟的结果不空洞,才能保证对接的准确性。
蛋白-蛋白对接如今比较著名的软件有Hex Protein Docking,ZDock,rDock,以及Rosetta等等。其中Hex Protein Docking算法复杂且每次耗时极长,ZDock一般作为前期的初对接,rDock为ZDock的升级算法,一般将ZDock对接后的前几个得分构象使用rDock进行进一步对接。今天想和大家介绍的是Rosetta的蛋白-蛋白对接模块,rosetta可谓一个蛋白质处理的全家桶,堪比Adobe在影像处理中的地位。
和其他对接一样,蛋白-蛋白对接也分为柔性对接和刚性对接两种。刚性对接中蛋白质的主链固定只进行移动或是转动,这种方法虽然精度较低,但是耗时短,计算量小,在实际操作中可以用蒙特卡洛算法产生的大量随机来达到期望的精度。柔性对接中两蛋白质的主链进行一定程度的摆动,更加符合实际情况,但是这种方法每次循环耗时长,不过可以用较小样本量达到预定的精度。
演示可以在安装文件的$ROSETTA3/demos/tutorials/Protein-Protein-Docking中查看,ROSETTA3为安装路径。
对接前结构准备
在进行对接前应对对接的位置有预先的了解并将要对接的蛋白放在预计的对接位置,手动生成预对接的结构文件。
Rosetta给的例子为Colicin-D与其抑制剂IMM之间的对接。需要提供两个蛋白结构的输入文件COL_D.pdb和IMM_D.pdb文件,以及1v74.pdb(预对接文件)文件在input_fi