DOCK 3.7 2014/09/25 FXa Tutorial: Difference between revisions

From DISI
Jump to navigation Jump to search
Line 56: Line 56:


The matching spheres are used during docking to orient the ligand into the binding pocket.   
The matching spheres are used during docking to orient the ligand into the binding pocket.   
The matching spheres are produced by converting the crystalagrafic ligand atoms into sphere and also using the spheres from sphgen.  
The matching spheres are produced by converting the crystallographic ligand atoms into sphere and also using the spheres from sphgen.  
    
    
[[File:match.png|thumb|center|375px|1FJS Receptor surface with matching spheres]]
[[File:match.png|thumb|center|375px|1FJS Receptor surface with matching spheres]]


This image was made with chimera.  to visulize the sphere there should be no sphere with a radius of zero.   
This image was made with chimera.  To visualize the sphere there should be no sphere with a radius of zero.   
the recommendation is to change "0.00" to "0.70".  Also you should remove the header from the file.  
the recommendation is to change "0.00" to "0.70".  Also you should remove the header from the file.  


Line 67: Line 67:
[[File:lowdiele.png|thumb|center|375px|1FJS Receptor surface with low dielectric spheres.]]
[[File:lowdiele.png|thumb|center|375px|1FJS Receptor surface with low dielectric spheres.]]


This image read shows the sph in the rec.crg.pdb.  This is the file on which qnifft runs the PB calculation.
This image read shows the spheres in the rec.crg.pdb.  This is the file on which qnifft runs the PB calculation.


=== pose visualization ==
=== pose visualization ==

Revision as of 11:51, 26 September 2014

Written by Trent E Balius and Crystal Nguyen on 2014/09/25

This part of the tutorial uses the release version of DOCK3.7-beta.1.0.1.

see the following for requirements: Install_DOCK_3.7

You must use a version of python (v2.7) that has the follow:

  • numpy
  • scipy
  • matplotlib
  • mysql-python
  • biopython

If you are using a virtual enviorment for python do something like this:

 source /home/tbalius/zzz.virtualenv/virtualenv-1.9.1/myVE/bin/activate.csh

to leave use the deactivate command:

 deactivate

making ligand databases.

Get smiles:

for factor Xa it is part of the DUDE database:

  wget http://dude.docking.org/targets/fa10/actives_final.ism
  wget http://dude.docking.org/targets/fa10/decoys_final.ism

We can download the isomorphic smiles from the dude webpage.

here is a webserver to gerenate decoys.

 http://dude.docking.org/generate

MORE TO COME HERE ON DATABASE PREP...

Here is what I think is what needs to be done.

mkdir ligands
bash $DOCKBASE/ligand/generate/build_smiles_ligand.sh ../actives_final.ism
python $DOCKBASE/ligand/finish/db2end-makedata.py ligands
mkdir decoys
cd decoys
bash $DOCKBASE/ligand/generate/build_smiles_ligand.sh ../decoys_final.ism
python $DOCKBASE/ligand/finish/db2end-makedata.py decoys

Receptor Preparation

  1. Splite your pdb into a rec.pdb and xtal-lig.pdb.
    1. you can do this by running this script:

The program sphgen will produce spheres that fill the nooks and crannies on the protein surface.

1FJS Receptor surface with all spheres from sphgen

The matching spheres are used during docking to orient the ligand into the binding pocket. The matching spheres are produced by converting the crystallographic ligand atoms into sphere and also using the spheres from sphgen.

1FJS Receptor surface with matching spheres

This image was made with chimera. To visualize the sphere there should be no sphere with a radius of zero. the recommendation is to change "0.00" to "0.70". Also you should remove the header from the file.

A final important set of spheres are the low dielectric spheres. These spheres are used during the PB calculation to force the binding site to have low dielectric.

1FJS Receptor surface with low dielectric spheres.

This image read shows the spheres in the rec.crg.pdb. This is the file on which qnifft runs the PB calculation.

= pose visualization

use chimera:

This part of the tutorial uses the GIST development version of DOCK

Put your rec.pdb and xtal-lig.pdb in the same frame as you simulation.

This chimera script might help.