Difference between revisions of "DOCK 3.7 2015/04/15 abl1 Tutorial"

Jump to: navigation, search
Line 65: Line 65:
  wget http://autodude.docking.org/abl1/actives_final.ism
  wget http://autodude.docking.org/abl1/actives_final.ism
  mv actives_final.ism ligands.ism
  mv actives_final.ism ligands.ism
== download enrichment databases ==
== run be_blasti.py==
== run be_blasti.py==

Revision as of 14:11, 15 April 2015

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

This is for a Linux environment and the scripts assume that you are running on SGE queueing system.

set up directories and get databases

Create directory called "RotationProject"

create a python file called "autodude_db_download.py"

# this gets the database from the autodude webpage

import sys, os
import urllib

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

print "url = " + url


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


for line in splitpage:
   if system in line: 
      file = line.replace('"',' ').split()[2]
      print url+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

make directories for ligands and decoys and move the corresponding files into those directories

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

download the ligand and decoy isomeric 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

run be_blasti.py

run blastermaster.py

run enrichment calucaltions