DOCK 3.7 2015/04/15 abl1 Tutorial

From DISI
Revision as of 21:05, 15 April 2015 by TBalius (talk | contribs)
Jump to navigation Jump to search

This tutoral use the 3.7.2 beta version of dock release on XXX.

set up directories and get databases

Create directory called "RotationProject"

create a file called "autodude_db_download.py"

# this get the database from the autodude webpage

import sys, os
import urllib

system = 'abl1'
url = 'http://autodude.docking.org/dude_e_db2/'

print "url = " + url

#page=requests.get(url)

webfile = urllib.urlopen(url)
page    = webfile.read()
webfile.close()

splitpage=page.split('\n')

for line in splitpage:
   if system in line: 
      file = line.replace('"',' ').split()[2]
      print url+file
      urllib.urlretrieve(url+file,file)

     # exit()

This python script will download the dockable db2 databases from the autodude webpage.python /mnt/nfs/home/rstein/RotationProject/autodude_db_download.py

make a subdirectory called databases:

mkdir databases

go inside.

cd databases


mkdir decoys 
mv decoys*db2.gz decoys
mkdir ligands 
mv ligands*db2.gz ligands

download the ligand isomorfic smiles file:

wget http://autodude.docking.org/abl1/decoys_final.ism
mv decoys_final.ism decoys.ism

note that the scripts expect the name to be decoys.ism, so we changed the name.

wget http://autodude.docking.org/abl1/actives_final.ism
mv actives_final.ism ligands.ism

download enrichment databases

run be_blasti.py

run blastermaster.py

run enrichment calucaltions