get.pdb(ids, path = ".", URLonly=FALSE, overwrite = FALSE, gzip = FALSE, split = FALSE, format = "pdb", verbose = TRUE, ncore = 1, ...)
pdbsplit
funciton will be called to split pdb files
into separated chains. ncore>1
requires package parallel installed. pdbsplit
function. Downloads PDB coordinate files from the RCSB Protein Data Bank.
This is a basic function to automate file download from the PDB.
Grant, B.J. et al. (2006) Bioinformatics 22, 2695--2696.
For a description of PDB format (version3.3) see: http://www.wwpdb.org/documentation/format33/v3.3.html.
# PDB server connection required - testing excluded ## PDB file paths get.pdb( c("1poo", "1moo"), URLonly=TRUE )[1] "http://www.rcsb.org/pdb/files/1poo.pdb" [2] "http://www.rcsb.org/pdb/files/1moo.pdb"## These URLs can be used by 'read.pdb' pdb <- read.pdb( get.pdb("5p21", URL=TRUE) ) summary(pdb)Call: read.pdb(file = get.pdb("5p21", URL = TRUE)) Total Models#: 1 Total Atoms#: 1567, XYZs#: 4701 Chains#: 1 (values: A) Protein Atoms#: 1323 (residues/Calpha atoms#: 166) Nucleic acid Atoms#: 0 (residues/phosphate atoms#: 0) Non-protein/nucleic Atoms#: 244 (residues: 213) Non-protein/nucleic resid values: [ GNP (1), HOH (211), MG (1) ] + attr: atom, xyz, seqres, helix, sheet, calpha, remark, call## Download PDB file ## get.pdb("5p21")