CIFSUP, is a program to superimpose one set of atomic coordinates onto another.
These coordinates may be low-resolution bead models and/or high resolution NMR
or x-ray crystal structures.
The program represents each input structure as an ensemble of points, then
minimised the distance between them, where the distance measure is defined by
the selected method. The following methods are
implemented:
The movable structure to align against the static template,
either in PDB or mmCIF format.
CIFSUP recognizes the following command-line options.
Short Option
Long Option
Description
--template-model=<ID>
Model ID in the template file; default: 1
--movable-model=<ID>
Model ID in the movable file; default: 1
--selection=<S>
One of ALL (all atoms), BACKBONE (only the backbone/CA atoms),
REGRID (convert both inputs to DAM) or SHELL (use only the outermost
atoms); default: ALL
Specify an output FILE where to write the moved and rotated model to.
By default this corresponds to the basename of the second file plus r.pdb. See also below.
CIFSUP does not have any runtime output. If the output file
is undefined, the transformed coordinates will be written to a (new) file path created
from the basename of the MOVABLE argument plus r.cif.
CIFSUP requires two atomic coordinate files as input, a template structure and
a target structure for alignment/superposition. Both may be either in PDB or
mmCIF format.
CIFSUP creates a single output file for the superimposed target structure. By
default the output filename is taken from the input MOVABLE
structure but is appended with 'r', eg.
$ CIFSUP file1.pdb file2.pdb
yields file2r.cif as output. A different filename, and file format, can
be specified with the --output option.
Use CIFSUP to align a dummy atom model on a crystal structure, using the
options to select the method that corresponds to
slow mode in SUPCOMB and writing the output in PDB format:
Use CIFSUP to align a dummy atom model on a crystal structure, using the
options to select the method that corresponds to
slow mode in SUPALM and writing the output in PDB format:
As above, but use --method=rmsd. Note that this requires a one-to-one
correspondence of ATOMS. This is generally not the case for dummy atom or
dummy residue models.
To superimpose 10 structures onto a single template, on linux, in bash syntax:
$ for i in `seq 1 10` ; do CIFSUP template.pdb file-$i.pdb; done
If one has unnumbered files to superimpose, on linux, in bash syntax:
$ for file in "a.pdb b.pdb c.pdb d.pdb" ; do CIFSUP template.pdb $file; done
where template.pdb is the reference/template structure.
Note that when one has a large collection of structures for superposition, the program DAMAVER is optimised for a fast superposition of multiple structures, automatically choosing the most representative structure.