NAME
king-probe - Evaluate and visualize protein interatomic packing
DESCRIPTION
Syntax: probe input.pdb >> out.kin
-
or: probe [flags] "src pattern" ["target pattern"] pdbfiles... >> out.kin
Flags:
- -SElf
-
self intersection: src -> src (default)
- -Both
-
intersect both ways: src <=> targ
- -ONce
-
single intersection: src -> targ
- -OUt
-
external van der Waals surface of src (solvent contact surface)
- -AUTObondrot filename
-
read and process an autobondrot file
-
shortcuts:
-
<<NO FLAGS>>same as: -4H -mc -het -self "altA ogt33"
- -DEFAULTs
-
same as: <<NO FLAGS>>, but allows some other flags
-SCSurface same as: -drop -rad1.4 -out "not water"
- -EXPOsed
-
same as: -drop -rad1.4 -out (note: user supplies pattern)
- -ASurface
-
same as: -drop -rad0.0 -add1.4 -out "not water"
- -ACCESS
-
same as: -drop -rad0.0 -add1.4 -out (note: user supplies pattern)
-SCAN0 same as: -4H -mc -self "alta blt40 ogt33"
-SCAN1 same as: -4H -once "sc alta blt40 ogt33" "alta blt40 ogt65,(not water ogt33)"
- -DUMPAtominfo
-
count the atoms in the selection: src
-
(note that BOTH and ONCE require two patterns while
-
OUT, SELF and DUMPATOMINFO require just one pattern)
- -Implicit
-
implicit hydrogens
- -Explicit
-
explicit hydrogens (default)
- -DEnsity#
-
set dot density (default 16 dots/sq A)
- -Radius#.#
-
set probe radius (default 0.25 A)
- -ADDvdw#.#
-
offset added to Van der Waals radii (default 0.0)
-SCALEvdw#.# scale factor for Van der Waals radii (default 1.0)
- -COSCale#.#
-
scale C=O carbon Van der Waals radii (default 0.94)
- -SPike
-
draw spike instead of dots (default)
- -SPike#.#
-
set spike scale (default=0.5)
- -NOSpike
-
draw only dots
-HBRegular#.# max overlap for regular Hbonds(default=0.6)
-HBCharged#.# max overlap for charged Hbonds(default=0.8)
- -Keep
-
keep nonselected atoms (default)
- -DRop
-
drop nonselected atoms
- -LIMit
-
limit bump dots to max dist when kissing (default)
- -NOLIMit
-
do not limit bump dots
- -LENs
-
add lens keyword to kin file
- -NOLENs
-
do not add lens keyword to kin file (default)
- -MC
-
include mainchain->mainchain interactions
- -HETs
-
include dots to non-water HET groups (default)
- -NOHETs
-
exclude dots to non-water HET groups
- -WATers
-
include dots to water (default)
- -NOWATers
-
exclude dots to water
- -WAT2wat
-
show dots between waters
- -DUMPH2O
-
include water H? vectorlist in output
- -4H
-
extend bond chain dot removal to 4 for H (default)
- -3
-
limit bond chain dot removal to 3
- -2
-
limit bond chain dot removal to 2
- -1
-
limit bond chain dot removal to 1
-IGNORE "pattern" explicit drop: ignore atoms selected by pattern
- -DOCHO
-
recognize CH..O Hbonds
- -CHO#.#
-
scale factor for CH..O Hbond score (default=0.5)
- -PolarH
-
use short radii of polar hydrogens (default)
- -NOPolarH
-
do not shorten radii of polar hydrogens
-NOFACEhbond do not identify HBonds to aromatic faces
-Name "name" specify the group name (default "dots")
- -DOTMASTER
-
group name used as extra master={name} on lists
- -NOGroup
-
do not generate @group statement in .kin format output
- -KINemage
-
add @kinemage 1 statement to top of .kin format output
- -Countdots
-
produce a count of dots-not a dotlist
- -Unformated
-
output raw dot info
-
name:pat:type:srcAtom:targAtom:mingap:gap:spX: spY:spZ:spikeLen:score:stype:ttype:x:y:z:sBval:tBval:
- -OFORMAT
-
output dot info formatted for display in O
- -XVFORMAT
-
output dot info formatted for display in XtalView
- -ONELINE
-
output one line :contacts:by:severity:type:
- -GAPcolor
-
color dots by gap amount (default)
- -ATOMcolor
-
color dots by atom type
- -BASEcolor
-
color dots by nucleic acid base type
- -COLORBase
-
color dots by gap and nucleic acid base type
-OUTCOLor "name" specify the point color for -OUT (default "gray")
- -GAPWeight#
-
set weight for scoring gaps (default 0.25)
-BUMPWeight# set relative scale for scoring bumps (default 10.0)
- -HBWeight#
-
set relative scale for scoring Hbonds (default 4.0)
- -DIVLow#.#
-
Division for Bump categories (default -0.4)
- -DIVHigh#.#
-
Division for Contact categories (default 0.25)
- -MINOCCupancy#.#
-
Occupancy below this is same as zero (default 0.02)
- -ELEMent
-
add master buttons for different elements in kin output
- -NOHBOUT
-
do not output contacts for HBonds
- -NOCLASHOUT
-
do not output contacts for clashes
- -NOVDWOUT
-
do not output contacts for van der Waals interactions
- -ONLYBADOUT
-
onlybadout output bad clashes (severe overlap contacts)
- -SUMMARY
-
output summary list of contacts and clashes
- -ONELINE
-
output summary list on oneline
- -NOTICKs
-
do not display the residue name ticker during processing
- -STDBONDs
-
assume only standard bonding patterns in standard residues
- -NOPARENT
-
do not bond hydrogens based on table of parent heavy atoms
- -SEGID
-
use the PDB SegID field to descriminate between residues
- -OLDU
-
generate old style -u output: kissEdge2BullsEye, etc
- -VErbose
-
verbose mode (default)
- -REFerence
-
display reference string
- -CHANGEs
-
display a list of program changes
- -Quiet
-
quiet mode
- -Help
-
show expanded help notice (includes other flags)
- -VERSION
-
one line version to stdout
Pattern elements: (should be put in quotes on the command line)
- FILE#
-
within file #
- MODEL#
-
within model #
- CHAINaa
-
within chain a
- SEGaaaa
-
segment identifier aaaa (where _ represents blank)
- ALTa
-
alternate conformation a
- ATOMaaaa
-
atom name aaaa (where _ represents blank) (all 4 characters are used so H would be ATOM_H__)
- RESaaa
-
residue aaa
- #
-
residue #
- #a
-
residue #, insert a
- #-#
-
residue range # (insert codes ignored)
- a
-
residue type by one letter codes (eg. y)
- aaa
-
residue type by three letter codes (eg. tyr)
-
ALL,PROTEIN,MC,SC,BASE,ALPHA,BETA,NITROGEN,CARBON, OXYGEN,SULFUR,PHOSPHORUS,HYDROGEN,METAL,POLAR, NONPOLAR,CHARGED,DONOR,ACCEPTOR,AROMATIC,METHYL, HET,WATER,DNA,RNA
-
all or a subset of the atoms
- OLT#
-
Occupancy less than # (integer percent)
- OGT#
-
Occupancy greater than # (integer percent)
- BLT#
-
B-value less than # (integer)
- BGT#
-
B-value greater than # (integer)
- INSa
-
Insert code a (where _ represents blank)
- WITHIN #.# OF #.#, #.#, #.#
-
atoms within distance from point
-
Patterns can be combined into comma separated lists such as "trp,phe,tyr" meaning TRP or PHE or TYR.
-
Patterns that are sepatated by blanks must all be true such as "chainb 1-5" meaning residues 1 to 5 in chain B.
-
You can also group patterns with parenthesis, separate multiple patterns with | meaning 'or' and choose the complement with NOT as in "not file1" meaning not in file 1.
-
An autobondrot file is similar to other PDB input files but it includes information identifying atoms subject to rotations and other transformations.
-
Example autobondrot file fragment showing Calpha-Cbeta bond rotation and a periodic torsion penalty function for this rotation
- ATOM
-
1 CB TYR 61 34.219 17.937 4.659 1.00 0.00
- bondrot:chi1:78.7:
-
0:359:5:33.138:18.517: 5.531:34.219:17.937: 4.659
-
cos:-3:60:3: ATOM 1 1HB TYR 61 34.766 18.777 4.206 1.00 0.00 ATOM 1 2HB TYR 61 34.927 17.409 5.315 1.00 0.00 ATOM 1 CG TYR 61 33.836 16.989 3.546 1.00 0.00 ...
-
Autobondrot commands use colons to separate values Transformations: BONDROT:id:currAng:start:end:stepSz:x1:y1:z1:x2:y2:z2
- TRANS:
-
id:currpos:start:end:stepSz:x1:y1:z1:x2:y2:z2
- NULL
-
# dummy
- Bias functions:
-
COS:scale:phaseOffset:frequency POLY:scale:offset:polynomialDegree CONST:value
- Branching:
-
SAVE and RESTORE or "(" and ")"
-
(e.g. to rotate each Chi and the methyls for isoleucine the
-
sequence is: rotChi1/SAVE/rotChi2/rotCD1/RESTORE/rotCG2)
-
Set orientation: GO:angle1:angle2:... Include files: @filename Comments: # comment text
probe: version 2.13.110909, Copyright 1996-2011, J. Michael Word