PDB File Coordinate Superposition

Usage

pdbfit(...)
"pdbfit"(pdb, inds = NULL, ...)
"pdbfit"(pdbs, inds = NULL, outpath = NULL, ...)

Arguments

pdb
a multi-model pdb object of class "pdb", as obtained from read.pdb.
pdbs
a list of class "pdbs" containing PDB file data, as obtained from read.fasta.pdb or pdbaln.
inds
a list object with a ‘xyz’ component with indices that selects the coordinate positions (in terms of x, y and z elements) upon which fitting should be based. This defaults to all equivalent non-gap positions for function pdbfit.pdbs, and to all calpha atoms for function pdbfit.pdb.
outpath
character string specifing the output directory for optional coordinate file output. Note that full files (i.e. all atom files) are written, seebelow.
...
extra arguments passed to fit.xyz function.

Description

Protein Databank Bank file coordinate superposition with the Kabsch algorithm.

Details

The function pdbfit is a wrapper for the function fit.xyz, wherein full details of the superposition procedure are documented.

Input to pdbfit.pdbs should be a list object obtained with the function read.fasta.pdb or pdbaln. See the examples below.

For function pdbfit.pdb the input should be a multi-model pdb object with multiple (>1) frames in the ‘xyz’ component.

The reference frame for supperposition (i.e. the fixed structure to which others are superposed) is the first entry in the input "pdbs" object. For finer control use fit.xyz.

Value

Returns moved coordinates.

References

Grant, B.J. et al. (2006) Bioinformatics 22, 2695--2696.

Kabsch Acta Cryst (1978) A34, 827--828.

Examples

#files <- get.pdb(c("4q21","5p21"), URLonly=TRUE) files <- get.pdb(c("4q21","5p21"), path=tempdir(), overwrite=TRUE) pdbs <- pdbaln(files)
Reading PDB files: /tmp/RtmpTDihxb/4q21.pdb /tmp/RtmpTDihxb/5p21.pdb .. Extracting sequences pdb/seq: 1 name: /tmp/RtmpTDihxb/4q21.pdb pdb/seq: 2 name: /tmp/RtmpTDihxb/5p21.pdb
xyz <- pdbfit(pdbs) # Superpose again this time outputing all-atom PDBs to disc #xyz <- pdbfit( pdbs, outpath="fitted" )

See also

pdbaln, read.fasta.pdb, fit.xyz, rmsd, read.pdb

Author

Barry Grant