MMTSB
Tool Set Documentation

Difference between revisions of "vicinity.pl"

From MMTSB
Jump to: navigation, search
Line 27: Line 27:
 
relief but no significant structural rearrangements beyond the immediate
 
relief but no significant structural rearrangements beyond the immediate
 
loop neighborhood. Force constants used in the second line of output
 
loop neighborhood. Force constants used in the second line of output
are by default 0.5 for residues within the soft cutoff and 10.0 kcal/mol
+
are by default 0.5 kcal/mol for residues within the soft cutoff and 10.0 kcal/mol
 
beyond. They can be changed by appending the new force constant to the cutoff
 
beyond. They can be changed by appending the new force constant to the cutoff
 
in the <B>-hard</B> and <B>-soft</B> options.<BR>
 
in the <B>-hard</B> and <B>-soft</B> options.<BR>

Revision as of 02:49, 28 July 2009

Usage

usage:   vicinity.pl [options] [file(s)]
options: [-l min:max[=min:max=...]]
         [-hard value[:force]]
         [-soft value[:force]]
         [-f filelist]
         [-[no]margin]

Show source


Description

This script is used for loop or fragment modeling to find residues that are located within a cutoff of a given residue subset in a larger protein. The output consists of two lines. The first line contains a residue specification for all residues that are within the cutoff (including the loop residues) that can be used with -l as command line argument in other MMTSB tools. The second line contains a residue specification including restraint force constants for use in minCHARMM.pl, for example, to run minimizations where only the neighboring residues of the modeled loop or fragment are kept flexible. The loop residues are not included in the second residue list.

One or more PDB files (with different conformations for the same protein) are read as input. The loop or fragment residues are specified with -l as a list of residue ranges. The cutoff is changed with -hard (default value: 16.0 A). A second, soft cutoff is set to 12.0 A by default and can be changed with -soft. Residues within this smaller cutoff are expected to be restrained with a smaller force constant while residues beyond the soft cutoff but within the main cutoff are restrained to a much larger extent to allow for stress relief but no significant structural rearrangements beyond the immediate loop neighborhood. Force constants used in the second line of output are by default 0.5 kcal/mol for residues within the soft cutoff and 10.0 kcal/mol beyond. They can be changed by appending the new force constant to the cutoff in the -hard and -soft options.

The -f option is available to read a list of files with different protein conformations from an external file (or standard input if - is given as the file name).

Options

-help 
usage information


Examples

vicinity.pl -l 71:81 1poa.exp.pdb
determines the residues within the default cutoffs of the loop at 71 through 81

1:19=21=28=39:40=42:44=46:55=57:103
16:19_10=21=28=39:40=42:44=46=48:49=52:53=55=57:61=63:64=103=1:15_0.5=47=50:51=54=62=65:70=82:102


vicinity.pl -l 71:81 -hard 12:1.0 -soft 10:0.2 1poa.exp.pdb
determine the residues within a 12 A from the loop at 71 through 81. The soft cutoff is also reduced to 10 A and modified force constants of 0.2 and 1.0 kcal/mol are used in the output.

1:15=47=50:51=54=62=64:102
1:2_1.0=15=47=50:51=54=62=64=102=3:14_0.2=65:70=82:101