Secangkir air putih dan puisi hari ini …

Sebenarnya mau menulis tentang hari AIDS dikaitkan dengan berita di kompas.com yang berubah judul di malam hari ini, namun tidak jadi karena harus merampungkan puisi hari ini supaya besok bisa ambil day off. Day off? Iyah … menghabiskan jatah cuti tahunan

Dan ufff … akhirnya kelar juga .. 3 puisi di hari ini … dan masih ada satu puisi yang belum selesai karena lisensi expired … (lol)

mae-pdb_output.sh
#!/bin/sh
export MYAPP=~/shared/apps/schrodinger/utilities/pdbconvert
if [ $1 = “enade” ]
then
sed ‘s/,/ /g’ report.csv | awk ‘{print $5}’ | grep -Ev ‘File’ > files.lst
mkdir result-temp
cd result-temp
cp ../*-out_pv-*.mae .
mv ../files.lst .
pdblist=files.lst
DIR=`cat $pdblist`
for enade in $DIR; do
$MYAPP -imae $enade -opdb $enade.pdb
cp $enade.pdb ../
done
cd ..
rm -R result-temp
else
echo “Convert mae to pdb from IFD output”
echo “Usage: ./mae-pdb_output.sh enade”
fi

pdbmol2_ligand.sh
#!/bin/sh
export MYAPP=/usr/local/moe-2008.10/bin/moebatch
if [ $1 = “enade” ]
then
mkdir result-temp
cd result-temp
cp ../*.mae.pdb .
ls *.mae.pdb > file.lst
sed ‘s/.mae.pdb//g’ file.lst > files.lst
pdblist=files.lst
DIR=`cat $pdblist`
for enade in $DIR; do
grep “UNK” $enade.mae.pdb > ligand_$enade
cp ligand_$enade temp.pdb
$MYAPP -script ~/scripts/pdbmol2.svl
rm temp.pdb
mv temp.mol2 ligand_$enade.mol2
cp ligand_$enade.mol2 ../
done
cd ..
rm -R result-temp
else
echo “Convert the ligand from pdb to mol2”
echo “Usage: ./pdbmol2_ligand.sh enade”
fi

pdbmol2_pocket_2.sh
#!/bin/sh
export MYAPP=/usr/local/moe-2008.10/bin/moebatch
if [ $1 = “enade” ]
then
mkdir result-temp
cd result-temp
cp ../*.mae.pdb .
ls *.mae.pdb > file.lst
sed ‘s/.mae.pdb//g’ file.lst > files.lst
pdblist=files.lst
DIR=`cat $pdblist`
for enade in $DIR; do
grep “A 328 \|A 329 \|A 332 \|A 333 \|A 336 \|A 337 \|A 340 \|A 456 \|A 460 ” $enade.mae.pdb > pocket_$enade
cp pocket_$enade temp.pdb
$MYAPP -script ~/scripts/pdbmol2.svl
rm temp.pdb
mv temp.mol2 pocket_$enade.mol2
cp pocket_$enade.mol2 ../
done
cd ..
rm -R result-temp
else
echo “Convert the pocket from pdb to mol2”
echo “Usage: ./pdbmol2_pocket.sh enade”
fi

~DosGil lapar betul sekarang…

Iklan