Advertisement
ProzacR

run_qvina

Oct 8th, 2020
1,168
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Bash 0.75 KB | None | 0 0
  1. #!/bin/bash
  2.  
  3. program="/home/vytautas/bin/qvina-master/bin/qvina-w"
  4. receptors="Cxs/"
  5. ligands="linear_terpenes/"
  6. #ligandai yra: *[0-9].pdbqt
  7.  
  8. for receptor_pdbqt in ${receptors}*.pdbqt;
  9. do
  10.         #echo $receptor_pdbqt
  11.     for file in ${ligands}*[0-9].pdbqt;
  12.     do
  13.         file $file;
  14.         #$program --receptor $receptor_pdbqt --ligand $file --center_x -45 --center_y 3 --center_z 20 --size_x 60 --size_y 40 --size_z 60 --out ${file%.pdbqt%/}_${receptor_pdbqt#Cxs/}_out.pdbqt;
  15.     done
  16.         for file in ${ligands}*_${receptor_pdbqt#Cxs/}_out.pdbqt;
  17.     do
  18.         #echo $file;
  19.         grep -m 1 "REMARK  Name =" $file | awk '{print $4}' | tr -d '\n' && echo -n ", " && grep -m 1 "REMARK VINA RESULT" $file | awk '{print $4}';
  20.         done > ${receptor_pdbqt#Cxs/}_rez.csv
  21. done
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement