需求描述
使用AutoDock CrankPep or ADCP进行蛋白质多肽对接
硬件及系统配置
自用电脑型号如下:
电脑:Precision Tower 7810 (Dell Inc.)
CPU : Intel® Xeon® CPU E5-2686 v4 @ 2.30GHz
GPU: NVIDIA GeForce GTX 1070
Linux版本:Linux version 6.4.0-150600.23.30-default
Opensuse版本:Opensuse 15.4
安装AutoDock CrankPep v1.0
尝试过安装AutoDock CrankPep v1.1,这个版本使用micromamba来安装,结果由于网络问题,没有安装成功。最后转而安装AutoDock CrankPep v1.0,从源文件安装
下载软件
可以从如下链接下载源文件:https://ccsb.scripps.edu/adcp/downloads/,由于我在linux电脑上面安装,所以使用ADFRsuite 1.0 Linux 64 tarball installer直接进行下载。
安装软件
#解压软件
tar zxvf ADFRsuite_x86_64Darwin_1.0.tar.gz
#建立安装目录
sudo mkdir /soft/ADFRsuite
#进入源文件目录并安装
cd ADFRsuite_x86_64Darwin_1.0
sudo ./install.sh -d /soft/ADFRsuite -c 0
#创建环境变量
export PATH=/soft/ADFRsuite/bin:$PATH
export LD_LIBRARY_PATH=/soft/ADFRsuite/lib:$LD_LIBRARY_PATH
多肽对接 – 教程(3Q47)
下载教程文件
浏览器打开 https://ccsb.scripps.edu/adcp/download/1063/
unzip ADCP_tutorial_data.zip
cd ADCP_tutorial_data/3Q47
预处理分子结构
# 分子结构加氢
/soft/ADFRsuite/bin/reduce 3Q47_rec.pdb > 3Q47_recH.pdb
/soft/ADFRsuite/bin/reduce 3Q47_pep.pdb > 3Q47_pepH.pdb
# 预处理分子结构,并且将结构转换为PDBQT格式
/soft/ADFRsuite/bin/prepare_receptor -r 3Q47_recH.pdb
/soft/ADFRsuite/bin/prepare_ligand -l 3Q47_pepH.pdb
生成目标文件(3Q47.trg)
/soft/ADFRsuite/bin/agfr -r 3Q47_recH.pdbqt -l 3Q47_pepH.pdbqt -asv 1.1 -o 3Q47
上面的命令会根据多肽的位置(3Q47_pepH.pdbqt)定义蛋白质的口袋,并且向外拓展0.4 纳米的区域。如果您的多肽并不在口袋内部的话,我们建议您首先构建一个虚拟的配体分子,并且将这个分子手动移动到蛋白质的口袋区域(借助pymol或者Maestro等工具)
多肽对接
/soft/ADFRsuite/bin/adcp -t 3Q47.trg -s npisdvd -N 20 -n 1000000 -o 3Q47_redocking -ref 3Q47_pepH.pdb
adcp软件的参数如下:
重要的参数有 -N 20 以及 -n 1000000,意思是:进行20次独立的搜索,每次搜索经历1000000步骤,这两个值越大,代表对接的精度越高。默认值是-N, —nbRuns 50 以及 -n, –numSteps 2.5 million
usage: usage: python runADCP.py -s GaRyMiChEL -t rec.trg -o outputAutoDock CrankPepoptional arguments:-h, --help show this help message and exit-v, --version show program's version number and exit-s SEQUENCE, --sequence SEQUENCEinitialize peptide from sequence, lower case for coiland UPPER case for helix-p PARTITION, --partition PARTITIONpartition for starting from a mixture of helix/coilconformation, percentage(helix)=partition/100 notethis option will overwrite the CaSe in sequence-i INPUT, --input INPUTuse conformation from pdb file as input-t TARGET, --target TARGETa zipped file prepared with AGFR describing thereceptor-n NUMSTEPS, --numSteps NUMSTEPSmax step for one replica-N NBRUNS, --nbRuns NBRUNSnumber of replicas-c MAXCORES, --maxCores MAXCORES-o JOBNAME, --jobName JOBNAME-y, --dryRun print the first adcp command line and exit-cyc, --cyclic option for cyclic peptide through backbone-cys, --cystein option for cyclic peptide through CYS-S-S-CYS-O, --overwriteFiles overwrite existing output files silently-S SEEDVALUE, --seed SEEDVALUEseed for random number generator-nc NC, --natContacts NCnative contacts cutoff used in the clustering-rmsd RMSD, --rmsd RMSDbackbone rmsd cutoff used in the clustering-ref REF, --ref REF reference peptide structure for calculating rmsd andfnc