Center, to the coordinate origin, and orient, by principal axes, the coordinates of a given PDB structure or xyz vector.
orient.pdb(pdb, atom.subset = NULL, verbose = TRUE)
a pdb data structure obtained from read.pdb
or a vector of ‘xyz’ coordinates.
a subset of atom positions to base orientation on.
print dimension details.
Returns a numeric vector of re-oriented coordinates.
Grant, B.J. et al. (2006) Bioinformatics 22, 2695--2696.
# NOT RUN {
# PDB server connection required - testing excluded
pdb <- read.pdb( "1bg2" )
xyz <- orient.pdb(pdb)
#write.pdb(pdb, xyz = xyz, file = "mov1.pdb")
# Based on C-alphas
inds <- atom.select(pdb, "calpha")
xyz <- orient.pdb(pdb, atom.subset=inds$atom)
#write.pdb(pdb, xyz = xyz, file = "mov2.pdb")
# Based on a central Beta-strand
inds <- atom.select(pdb, resno=c(224:232), elety='CA')
xyz <- orient.pdb(pdb, atom.subset=inds$atom)
#write.pdb(pdb, xyz = xyz, file = "mov3.pdb")
# }
Run the code above in your browser using DataLab