12.12 D3Docking

还是想交一批对接任务,把两千多个成分都包进去。
但是口袋不能太多,目前想的是库中现在有62个蛋白,把阳性化合物/老药打分最高的那个口袋拿出来。
定位到smina-dock_multi.sh里选择口袋的部分
找一下老药的对接结果

1212

  • # 先交一个理论上最快的:3CL和PLpro 打分最高构象的口袋

  • # /home/zjxu/zzy/hsbd/dock/1212/test3-input-sdf
  • # 输入是2049个sdf,完全按照D3Docking 的脚本执行
  • (base) [zjxu@R820 test3-input-sdf]$ nohup sh do.sh &
  • [2] 68921

  • # /home/zjxu/zzy/hsbd/dock/1212/smina-dock-directly
  • # 直接交了两个smina 任务,输入就直接是2049个一起的mol2文件
  • (base) [zjxu@R820 ORF1ab_819-2763_Papain-like_proteinase+Dimer-Pocket1]$ nohup /home/zjxu/prog/smina-code/build/linux/release/smina --seed 0 --custom_scoring /home/zjxu/xbzhang/2019-nCov-final/score_add_Mg-O.txt --config /home/zjxu/xbzhang/2019-nCov-final/D3Pocket-v1/ORF1ab_819-2763_Papain-like_proteinase+Dimer/pocket-1.txt -r /home/zjxu/xbzhang/2019-nCov-final/preparation/protein-v1/ORF1ab_819-2763_Papain-like_proteinase+Dimer.pdbqt -l ../hsbd-2049mol-2d-3d.mol2 -o inputmol2-dockP1.sdf --cpu 8 &
  • [2] 64489

  • (base) [zjxu@R820 ORF1ab_3264-3569_3C-like_proteinase+Dimer-Pocket1]$ nohup /home/zjxu/prog/smina-code/build/linux/release/smina --seed 0 --custom_scoring /home/zjxu/xbzhang/2019-nCov-final/score_add_Mg-O.txt --config /home/zjxu/xbzhang/2019-nCov-final/D3Pocket-v1/ORF1ab_3264-3569_3C-like_proteinase+Dimer/pocket-1.txt -r /home/zjxu/xbzhang/2019-nCov-final/preparation/protein-v1/ORF1ab_3264-3569_3C-like_proteinase+Dimer.pdbqt -l ../hsbd-2049mol-2d-3d.mol2 -o inputmol2-dockP1.sdf --cpu 8 &
  • [3] 75019

1213

  • # 85.23 /home/databank_70t/zzy/project/hsbd/D3Targets-dock/workdir/smina-dock-directly/
  • # 直接交了两个smina 任务,输入就直接是2049个一起的mol2文件
  • nohup smina --cpu 30 --seed 0 --custom_scoring /home/databank_70t/zzy/database/2019-nCov-final/score_add_Mg-O.txt --config /home/databank_70t/zzy/database/2019-nCov-final/D3Pocket-v1/ORF1ab_819-2763_Papain-like_proteinase+Dimer/pocket-1.txt -r /home/databank_70t/zzy/database/2019-nCov-final/preparation/protein-v1/ORF1ab_819-2763_Papain-like_proteinase+Dimer.pdbqt -l ../hsbd-2049mol-2d-3d.mol2 -o dock-PLpro-hsbd2049.sdf &

  • nohup smina --cpu 30 --seed 0 --custom_scoring /home/databank_70t/zzy/database/2019-nCov-final/score_add_Mg-O.txt --config /home/databank_70t/zzy/database/2019-nCov-final/D3Pocket-v1/ORF1ab_3264-3569_3C-like_proteinase+Dimer/pocket-1.txt -r /home/databank_70t/zzy/database/2019-nCov-final/preparation/protein-v1/ORF1ab_3264-3569_3C-like_proteinase+Dimer.pdbqt -l ../hsbd-2049mol-2d-3d.mol2 -o dock-3CLpro-hsbd2049.sdf &
  • # 85.23 D3Docking 已部署完成 /home/databank_70t/zzy/project/hsbd/D3Targets-dock/workdir/
  • # mgtools 安装位置 /home/databank_70t/zzy/software/mgtools
  • # rdkit pip 安装
  • # /home/databank_70t/zzy/project/hsbd/D3Targets-dock/workdir/3cl-pl-all-conf-best-pkt 任务正在进行

  • # 85.24 smina 环境下差一个 rdkit
  • # 需要尝试关掉sleep ,速度会起飞
  • # mgtools 通过conda 安装(~/.condarc 复制的85.12)
  • # smina 本身自带
  • /home/databank/zzy/project/hsbd/D3Targets-dock/workdir/