![]()
|
LSQKAB (CCP4: Supported Program)NAMElsqkab - Work out and apply various transformations to coordinate files.SYNOPSISlsqkab xyzinm tomatch.brk xyzinf reference.brk
xyzout matched.brk Note: XYZINM has the aliases WORKCD and XYZINW DESCRIPTIONPrimary useOptimises the fit of a subset of atomic coordinates from one file (assigned to XYZINM) to the same subset of another file (assigned to XYZINF). The "working" coordinates in XYZINM are to be moved; th "reference" coordinates in XYZINF are fixed.The program assumes both sets of coordinates are in the PDB format. If there are two or more conformations, the first (labelled A) is chosen for comparison. The subset can be defined either by residue numbers with options for main chain only, side chain only or CAs, or by the position of the atom in the file. The centroid of the moving subset is first shifted to that of the fixed set and then a rotation matrix is derived using a method suggested by Kasbch (reference [1]). This rotation matrix is defined relative to the orthogonal coordinate system. The whole of the moving coordinate set can be transformed and output in an orientation matching that of the fixed reference set. The output file assigned to XYZOUT. You can give CELL dimensions to set up the CRYST1 and SCALEi lines, otherwise they will be copied from the XYZINM coordinate file. There are various output options described on the input cards. Secondary purpose - much better done in PDBSETThe rotating coordinate set contained in file XYZINM can be transformed by a given rotation matrix and translation (see ROTATE and TRANSLATE keywords), and then output in a new orientation to the file XYZOUT. KEYWORDED INPUTThe keywords are: FIT, MATCH, OUTPUT, RADIUS, ROTATE, TITLE, TRANSLATE, CELL, END. FIT [ ATOM | RESIDUE ] <arguments>Describes the rotating coordinate subset taken from XYZINM. Keyword must be followed by ATOM (alias WATOM) or RESIDUE (alias WRESIDUE).
For fitting the next card must be MATCH. MATCH <n1> TO <n2> [ CHAIN <chnam> ]Describes the reference coordinate subset taken from XYZINF.
CELL a b c [alpha beta gamma ncode]Default alpha beta gamma 90.0 90.0 90.0. Default ncode 1. Read cell dimensions and prepare CRYST1 & SCALE header records for XYZOUT. These will replace any cell dimensions read from th input working coordinates. The fractional translation vector will be given relative to these cell dimensions. ncode selects the orthogonalization matrix to be used. Code :- = 1 axes along a, c* x a, c* (Brookhaven standard, default) = 2 axes along b, a* x b, a* = 3 axes along c, b* x c, b* = 4 axes along a+b, c* x (a+b), c* = 5 axes along a*, c x a*, c ( Rollett ) = 6 axes along a, b*, a x b* = 7 axes along a*, b, a* x b (TNT convention,
RADIUS <radchk> [ <Xc> <Yc> <Zc> ]Use atoms in a sphere for fitting. radchk: All atoms more than radchk Angstroms from the reference centroid are excluded. radchk Xc Yc Zc: All atoms more than radchk Angstroms from the reference point Xc Yc Zc (in As) are excluded ROTATE [ EULERIAN / POLAR / DCS / MATRIX ]must be followed by EULERIAN or POLAR or DCS (direction cosines) or MATRix. All generate a matrix for rotation of XYZINM. EULERIAN alpha beta gamma POLAR omega phi kappa DCS dc1 dc2 dc3 kappa MATRIX nine elements: m11 m12 m13 m21 m22 m23 m31 m32 m33 TRANSLATE <t1> <t2> <t3>Must be followed by t1 t2 t3, the translation vector in Ås along PDB axes. Alternative: Input coordinates are rotated by matrix [R] then translated by vector
tr. TITLE <title>ENDTerminate input and begin program. REQUIRED FILES
KNOWN PROBLEMSThere has been a report that in some cases using the FIT RESIDUE CA options, the fitting has been ``reversed'' e.g. the commands FIT RESIDUE CA 7 TO 12 MATCH 7 TO 12result in the CA atom of residue 12 in the rotating set is fitted to the CA of residue 7 in the reference set and so on. The reported workaround is to split the FIT/MATCH commands, for example: FIT RESIDUE CA 7 TO 9 MATCH 7 TO 9 FIT RESIDUE CA 10 TO 12 MATCH 10 TO 12which apparently works in most cases. (Thanks to Toby Galbraith) REFERENCEEXAMPLESUNIX example script found in $CEXAM/unix/runnable/VMS example script found in $CEXAM/vms/Very complicated matchlsqkab xyzinf ../data/takaxp_model8_9.pdb xyzinm ../data/takap21.pdb xyzout {$SCRATCH}lsq_cat_residues.pdb DELTAS {$SCRATCH}lsq_cat_residues.delta << END-lsqkab title matching catalytic residues of 3A to 2.1A taka output XYZ output deltas fit RESIDU side 155 TO 156 CHAIN A MATCH RESIDU 155 TO 156 CHAIN A fit RESIDU side 210 TO 210 CHAIN A MATCH RESIDU 210 TO 210 CHAIN A fit RESIDU side 206 TO 206 CHAIN A MATCH RESIDU 206 TO 206 CHAIN A fit RESIDU side 122 TO 122 CHAIN A MATCH RESIDU 122 TO 122 CHAIN A fit RESIDU side 83 TO 83 CHAIN A MATCH RESIDU 83 TO 83 CHAIN A fit RESIDU side 168 TO 168 CHAIN A MATCH RESIDU 168 TO 168 CHAIN A fit RESIDU side 74 TO 74 CHAIN A MATCH RESIDU 74 TO 74 CHAIN A fit RESIDU side 35 TO 35 CHAIN A MATCH RESIDU 35 TO 35 CHAIN A fit RESIDU side 79 TO 79 CHAIN A MATCH RESIDU 79 TO 79 CHAIN A fit RESIDU side 340 TO 340 CHAIN A MATCH RESIDU 340 TO 340 CHAIN A fit RESIDU side 344 TO 344 CHAIN A MATCH RESIDU 344 TO 344 CHAIN A fit RESIDU side 296 TO 297 CHAIN A MATCH RESIDU 296 TO 297 CHAIN A fit RESIDU side 230 TO 233 CHAIN A MATCH RESIDU 230 TO 233 CHAIN A fit RESIDU side 209 TO 209 CHAIN A MATCH RESIDU 209 TO 209 CHAIN A end END-lsqkab |