Travel Depth - Revision history

Therese: 18 revisions

2012-10-08T20:33:13Z

18 revisions

← Older revision

Revision as of 20:33, 8 October 2012

Rgc: /* Pocket Picking */

2011-09-05T21:36:19Z

Pocket Picking

New page

I've been forwarding this email to various people, so I figured I should just put it here. This is documentation for the Travel Depth/CLIPPERS code developed in the Sharp lab.

{{TOCright}}

http://crystal.med.upenn.edu/travel_distance.html

Publications:

*http://dx.doi.org/10.1021/ci900397t
*http://dx.doi.org/10.1002/prot.22558
*http://dx.doi.org/10.1529/biophysj.108.135970
*http://dx.doi.org/10.1016/j.jmb.2006.07.022

= Pocket Examination =

Here's the brief rundown

setenv TDHOME /raid1/people/rgc/Source/travel_distance_code_3

goes in your .cshrc/.login file (and you have to have python+numpy in your path)

Then with a pdb file you run:

$TDHOME/src/tstTravelDepth.py $prefix.pdb

(if you want a finer grid size you can specify it after the pdb file,
the default is 1 angstrom/grid)

After this you'll have several new files, with which you can run the
pocket mapping. There are 2 sets of files, one with cavities
(completely separated from bulk solvent) and one where those have been
removed.

$TDHOME/src/tstConvexHull.py $TDHOME/bin/qhull $prefix.cav.tst
$TDHOME/src/tstPocketMap.py $prefix.cav.tst $prefix.cav.phi

for cavities or

$TDHOME/src/tstPocketMap.py $prefix.nocav.tst $prefix.nocav.phi

with no cavities.

The primary output of these files are $prefix.cav.tst.tree.tm3 files
which have data about all the pockets. These are tab-delimited files
you can load in excel/spreadsheet program.

To see the pockets/pick one manually, you can run the following script in pymol:

run ~rgc/Source/travel_distance_code_3/pymol/tstdisp.py

This gives you several new pymol commands to run. First

tstOpen $prefix.nocav.tst

and then:

tstDraw

which opens and draws the complete surface.

To change colors you can try for example:

tstDraw grad=mono

Now you can load a
specific pocket by number (the first column in the
$prefix.nocav.tst.tree.tm3 file), where X is your number:

tstDraw group=X

alternatively, you can run

tstDrawDeepest

which displays the deepest pocket (the pocket number is encoded as
part of the pymol name, right after the name of the file).
tstDrawParent goes up from the last drawn pocket, tstDrawParent N goes
up the tree N times. When you're happy with the pocket, you can look
at the $prefix.nocav.tst.tree.tm3 file and find that pocket id number
and look up the various properties it reports.

If you have a list of residues in the binding site, you want to format
it like this "ALA147A+GLY148A+HIS190A" or just "147A+148A+190A" or
just "147+148+190". Chain numbers should be used if the original file
had them, or left out if it didn't. Then you run

$TDHOME/src/tm3residues name ALA147+GLY148+HIS190 [any number of tm3 files]

or

$TDHOME/src/tm3residues num 147+148+190 [any number of tm3 files]

There will be a summary file name temp.output.txt which is basically
the rows of the various tm3 files for the pockets that matched the
best, with additional data about the overlap between them and the
binding site (I think, it's been awhile since I ran this). Anyway, if
you have problems, feel free to ask.

= Ligand Volume/Surface Area =

For those interested in finding the volume and surface area of ligands (or whole proteins, other macromolecules), this is easy by starting in a similar fashion (follow the first several steps above to get your environment setup correctly).

$TDHOME/src/tstTravelDepth.py $prefix.pdb

If you really have a ligand and want "more accurate" numbers, feel free to specify the grid spacing by:

$TDHOME/src/tstTravelDepth.py $prefix.pdb 0.5

where 0.5 is the number of Å per grid (lower is finer resolution). the next step is:

$TDHOME/src/tstSurfAreaVolume.py $prefix.nocav.tst [more *.nocav.tst files]

The output is a tab-delimited stream to standard out containing the volume and surface area of the entire molecule, be it ligand, protein, etc. The volume and surface area are based on the same radii definitions used in the surface generation, essentially Bondi radii (Bondi, A. (1964). van der Waals Volumes and Radii. J. Phys. Chem. 68, 441-451.), and due to this exotic radii may be missing and need to be added.

= Pocket Picking =

run

$TDHOME/src/tm3picker.py blah.tm3

to get a list of potential binding sites as one list of residues lining the active site per line.

$TDHOME/src/tm3pickerPdb.py blah.pdb blah.tm3

outputs one or more pdb files containing the residues in each pocket.

[[Category:Software]][[Category:Tutorials]]

Rgc at 00:13, 30 October 2010

2010-10-30T00:13:57Z

← Older revision		Revision as of 00:13, 30 October 2010
Line 108:		Line 108:

	The output is a tab-delimited stream to standard out containing the volume and surface area of the entire molecule, be it ligand, protein, etc. The volume and surface area are based on the same radii definitions used in the surface generation, essentially Bondi radii (Bondi, A. (1964). van der Waals Volumes and Radii. J. Phys. Chem. 68, 441-451.), and due to this exotic radii may be missing and need to be added.		The output is a tab-delimited stream to standard out containing the volume and surface area of the entire molecule, be it ligand, protein, etc. The volume and surface area are based on the same radii definitions used in the surface generation, essentially Bondi radii (Bondi, A. (1964). van der Waals Volumes and Radii. J. Phys. Chem. 68, 441-451.), and due to this exotic radii may be missing and need to be added.

			= Pocket Picking =

			run

			$TDHOME/src/tm3picker.py blah.tm3

			to get a list of potential binding sites as one list of residues lining the active site per line.

	[[Category:Software]][[Category:Tutorials]]		[[Category:Software]][[Category:Tutorials]]

Rgc: /* Ligand Volume/Surface Area */ explanation of environment setup

2010-03-29T20:51:16Z

Ligand Volume/Surface Area: explanation of environment setup

← Older revision		Revision as of 20:51, 29 March 2010
Line 95:		Line 95:
	= Ligand Volume/Surface Area =		= Ligand Volume/Surface Area =

	For those interested in finding the volume and surface area of ligands (or whole proteins, other macromolecules), this is easy by starting in a similar fashion		For those interested in finding the volume and surface area of ligands (or whole proteins, other macromolecules), this is easy by starting in a similar fashion (follow the first several steps above to get your environment setup correctly).

	$TDHOME/src/tstTravelDepth.py $prefix.pdb		$TDHOME/src/tstTravelDepth.py $prefix.pdb

Rgc: moving toc

2010-03-29T17:58:47Z

moving toc

← Older revision		Revision as of 17:58, 29 March 2010
Line 1:		Line 1:
	I've been forwarding this email to various people, so I figured I should just put it here. This is documentation for the Travel Depth/CLIPPERS code developed in the Sharp lab.		I've been forwarding this email to various people, so I figured I should just put it here. This is documentation for the Travel Depth/CLIPPERS code developed in the Sharp lab.

			{{TOCright}}

	http://crystal.med.upenn.edu/travel_distance.html		http://crystal.med.upenn.edu/travel_distance.html
Line 10:		Line 12:
	*http://dx.doi.org/10.1016/j.jmb.2006.07.022		*http://dx.doi.org/10.1016/j.jmb.2006.07.022

	~~{{TOCright}}~~

	= Pocket Examination =		= Pocket Examination =

Rgc: /* Ligand Volume/Surface Area */

2010-03-29T17:58:23Z

Ligand Volume/Surface Area

← Older revision		Revision as of 17:58, 29 March 2010
Line 101:		Line 101:
	$TDHOME/src/tstTravelDepth.py $prefix.pdb 0.5		$TDHOME/src/tstTravelDepth.py $prefix.pdb 0.5

	where 0.5 is the number of ~~grids~~ per Å. the next step is:		where 0.5 is the number of Å per grid (lower is finer resolution). the next step is:

	$TDHOME/src/tstSurfAreaVolume.py $prefix.nocav.tst [more *.nocav.tst files]		$TDHOME/src/tstSurfAreaVolume.py $prefix.nocav.tst [more *.nocav.tst files]

Rgc: /* Ligand Volume/Surface Area */ bondi radii description

2010-03-29T16:42:08Z

Ligand Volume/Surface Area: bondi radii description

← Older revision		Revision as of 16:42, 29 March 2010
Line 105:		Line 105:
	$TDHOME/src/tstSurfAreaVolume.py $prefix.nocav.tst [more *.nocav.tst files]		$TDHOME/src/tstSurfAreaVolume.py $prefix.nocav.tst [more *.nocav.tst files]

	The output is a tab-delimited stream to standard out containing the volume and surface area of the entire molecule, be it ligand, protein, etc.		The output is a tab-delimited stream to standard out containing the volume and surface area of the entire molecule, be it ligand, protein, etc. The volume and surface area are based on the same radii definitions used in the surface generation, essentially Bondi radii (Bondi, A. (1964). van der Waals Volumes and Radii. J. Phys. Chem. 68, 441-451.), and due to this exotic radii may be missing and need to be added.


	[[Category:Software]][[Category:Tutorials]]		[[Category:Software]][[Category:Tutorials]]

Rgc: adding toc

2010-03-29T16:40:24Z

adding toc

← Older revision		Revision as of 16:40, 29 March 2010
Line 9:		Line 9:
	*http://dx.doi.org/10.1529/biophysj.108.135970		*http://dx.doi.org/10.1529/biophysj.108.135970
	*http://dx.doi.org/10.1016/j.jmb.2006.07.022		*http://dx.doi.org/10.1016/j.jmb.2006.07.022

			{{TOCright}}

	= Pocket Examination =		= Pocket Examination =

Rgc: adding ligand volume description

2010-03-29T16:38:55Z

adding ligand volume description

@@ Line 88: / Line 88: @@
 binding site (I think, it's been awhile since I ran this). Anyway, if
 you have problems, feel free to ask.
 [[Category:Software]][[Category:Tutorials]]

Rgc: adding headings

2010-03-29T16:35:26Z

adding headings

← Older revision		Revision as of 16:35, 29 March 2010
Line 10:		Line 10:
	*http://dx.doi.org/10.1016/j.jmb.2006.07.022		*http://dx.doi.org/10.1016/j.jmb.2006.07.022

			= Pocket Examination =

	Here's the brief rundown		Here's the brief rundown