-
Notifications
You must be signed in to change notification settings - Fork 14
Description
Hi. I would like to write something about error I faced : "SOAP-potential file not found". This maybe helps someone.
The codes were:
../bin/frodock 3hfl_fv_ASA.pdb 3hfl_ly2_ASA.pdb -w 3hfl_fv_W.ccp4 -e 3hfl_fv_E.ccp4 --th 10 -d 3hfl_fv_DS.ccp4,3hfl_ly2_DS.ccp4 -o dock.dat
The outputs:
frodock> Name H convention: pdb
frodock> Force Field: CHARMM
frodock>
frodock> Reading receptor PDB
frodock> Receptor PDB center: -0.598341 -64.294624 -5.272957
frodock> Receptor accessibility read from receptor PDB
frodock> ...occupancy column seems to have ASA values
frodock> Creating Accessibility map from receptor PDB...
frodock> Accessibility map created from receptor PDB
frodock> Weights SOAP 0.700000 Tobi 15.000000
frodock> Reading SOAP-potential file: ../bin/soap.bin
frodock> �� SOAP-potential file not found
ERROR: ��
Then dock.dat file becomes empty, which breaks the code.
You need to add one line to eliminate the error.
- go to utils/frodock.py
- Add line just after frodock_gcc = '%s/bin/frodock_gcc' % FRODOCK (line 37)
write : soap='%s/bin/soap.bin' % FRODOCK - Go to line frodock execution line (it was 46, but now 47)
edit that code like that:
os.system(frodock_docking + ' receptor_ASA.pdb target_ASA.pdb -w receptor_W.ccp4 -e receptor_E.ccp4 '
'--th 10 -d receptor_DS.ccp4,target_DS.ccp4 -t E -o dock.dat --around %s -s %s' %(site, soap))
I hope it helps someone.