背景：
之前有讲过可以通过pymol统计氢键的信息，
由于pymol对pdbqt的支持不好
ADT中的脚本pdbqt_to_pdb.py
只能实现pdbqts_to_pdb的功能
所以我写了个shell脚本，实现pdbqts_to_pdbs功能，用这个脚本前必须配置一些环境变量。
======================================pdbqts_to_pdbs.sh
#!/bin/sh

#Function : convert the pdbqt format of  ligand database to pdb format
#author: Chen Zhaoqiang
#email: 744891290@qq.com
#date:2013.10.09
#required:vina_split in vina  and pdbqt_to_pdb.py in adt
#usage: pdbqts_to_pdbs.sh --input xxx.pdbqt


if [ $# -lt 2 ];then
   echo "usage: pdbqts_to_pdbs.sh --input xxx.pdbqt\n";
fi

if [ -d temppp ];then
   cd ./temppp
   rm ./* -f
else
   mkdir temppp
   cd ./temppp
fi
cp ../$2 ./
vina_split --input $2
rm ./$2
for file in $( ls );
do
   echo $file
   pdbqt_to_pdb.py -f $file -o $file'.pdb'
   echo  $file
       rm $file


done

#modify the name to short

rename .pdbqt.out_ligand "" ./*
rename .pdbqt "" ./*
let count=1
for file in $( ls );
do
   echo -e "MODEL $count\n" >>$2
   cat $file >>$2
   echo -e "ENDMDL\n" >>$2
   rm $file
   let count+=1
done
let count-=1
echo $count
rename .pdbqt.out.pdbqt _out_$count.pdb ./$2
cp ./* ../
cd ..
rm -f temppp -r
===========================================================
把这个脚本放到path中，实现和pdbqt_to_pdb.py同样的地位

####

pdb的库中必须要有MODEL ENDMDL的标签，

chimera这个软件不识别pdb文件中*，*不能代替一般字母

如果这个脚本不能用，先测试vina_split这个命令能不能使用。