export PATH=${PATH}:/home/preps/golovin/progs/bin
for i in nag_h nag_oh nag_nh2 nag_ph ;do
obgen ${i}.smi > ${i}.mol
babel -imol ${i}.mol -opdb ${i}.pdb 
babel -ipdb ${i}.pdb -omop ${i}.mop -xk "PM6" 
MOPAC2009.exe ${i}.mop 
babel -imopout ${i}.out -opdb ${i}1.pdb
/home/preps/golovin/progs/bin/prepare_ligand4.py -l ${i}1.pdb -o ${i}.pdbqt
vina --config vina.cfg --receptor seq3.pdbqt --ligand ${i}.pdbqt --out ${i}_seq3.pdbqt --log ${i}_seq3.log
vina --config vina.cfg --receptor seq3_rigid.pdbqt --flex seq3_flex.pdbqt --ligand ${i}.pdbqt --out ${i}_seq3_fr.pdbqt --log ${i}_seq3_fr.log
done