as.pdb. bio3d 2.3-0

Usage

as.pdb(...)
"as.pdb"(mol, ...)
"as.pdb"(prmtop, crd=NULL, inds=NULL, inds.crd=inds, ncore=NULL, ...)
"as.pdb"(pdb=NULL, xyz=NULL, type=NULL, resno=NULL, resid=NULL, eleno=NULL, elety=NULL, chain=NULL,  insert=NULL, alt=NULL, o=NULL, b=NULL, segid=NULL,  elesy=NULL, charge=NULL, verbose=TRUE, ...)

Arguments

...: arguments passed to and from functions.
mol: a list object of type "mol2" (obtained with read.mol2).
prmtop: a list object of type "prmtop" (obtained with read.prmtop).
crd: a list object of type "crd" (obtained with read.crd.amber).
inds: a list object of type "select" as obtained from atom.select. The indices points to which atoms in the PRMTOP object to convert.
inds.crd: same as the ‘inds’ argument, but pointing to the atoms in CRD object to convert. By default, this argument equals to ‘inds’, assuming the same number and sequence of atoms in the PRMTOP and CRD objects.
ncore: number of CPU cores used to do the calculation. ncore>1 requires package ‘parallel’ installed.
pdb: an object of class ‘pdb’ as obtained from read.pdb.
xyz: a numeric vector/matrix of Cartesian coordinates. If provided, the number of atoms in the new PDB object will be set to ncol(as.xyz(xyz))/3 (see as.xyz). If xyz is not provided the number of atoms will be based on the length of eleno, resno, or resid (in that order).
type: a character vector of record types, i.e. "ATOM" or "HETATM", with length equal to ncol(as.xyz(xyz))/3. Alternatively, a single element character vector can be provided which will be repeated to match the number of atoms.
resno: a numeric vector of residue numbers of length equal to ncol(as.xyz(xyz))/3.
resid: a character vector of residue types/ids of length equal to ncol(as.xyz(xyz))/3. Alternatively, a single element character vector can be provided which will be repeated to match the number of atoms.
eleno: a numeric vector of element/atom numbers of length equal to ncol(as.xyz(xyz))/3.
elety: a character vector of element/atom types of length equal to ncol(as.xyz(xyz))/3. Alternatively, a single element character vector can be provided which will be repeated to match the number of atoms.
chain: a character vector of chain identifiers with length equal to ncol(as.xyz(xyz))/3. Alternatively, a single element character vector can be provided which will be repeated to match the number of atoms.
insert: a character vector of insertion code with length equal to ncol(as.xyz(xyz))/3.
alt: a character vector of alternate record with length equal to ncol(as.xyz(xyz))/3.
o: a numeric vector of occupancy values of length equal to ncol(as.xyz(xyz))/3. Alternatively, a single element numeric vector can be provided which will be repeated for to match the number of atoms.
b: a numeric vector of B-factors of length equal to ncol(as.xyz(xyz))/3. Alternatively, a single element numeric vector can be provided which will be repeated to match the number of atoms.
segid: a character vector of segment id of length equal to ncol(as.xyz(xyz))/3. Alternatively, a single element character vector can be provided which will be repeated to match the number of atoms.
elesy: a character vector of element symbol of length equal to ncol(as.xyz(xyz))/3. Alternatively, a single element character vector can be provided which will be repeated to match the number of atoms.
charge: a numeric vector of atomic charge of length equal to ncol(as.xyz(xyz))/3.
verbose: logical, if TRUE details of the PDB generation process is printed to screen.

Description

Convert Tripos Mol2 format, or Amber parameter/topology and coordinate data to PDB format.

Details

This function converts Tripos Mol2 format, Amber formatted parameter/topology (PRMTOP) and coordinate objects, and vector data to a PDB object.

While as.pdb.mol2 and as.pdb.prmtop converts specific objects to a PDB object, as.pdb.default provides basic functionality to convert raw data such as vectors of e.g. residue numbers, residue identifiers, Cartesian coordinates, etc to a PDB object. When pdb is provided the returned PDB object is built from the input object with fields replaced by any input vector arguments. e.g. as.pdb(pdb, xyz=crd) will return the same PDB object, with only the Cartesian coordinates changed to crd.

Value

atom: a data.frame containing all atomic coordinate ATOM data, with a row per ATOM and a column per record type. See below for details of the record type naming convention (useful for accessing columns).
xyz: a numeric matrix of ATOM coordinate data of class xyz.
calpha: logical vector with length equal to nrow(atom) with TRUE values indicating a C-alpha “elety”.
call: the matched call.

References

Grant, B.J. et al. (2006) Bioinformatics 22, 2695--2696. http://ambermd.org/formats.html

Examples

## Vector(s) to PDB object
pdb <- as.pdb(resno=1:6, elety="CA", resid="ALA", chain="A")


 Summary of PDB generation:
 .. number of atoms in PDB determined by 'resno'

 ..        6 atom(s) from 'string' selection 
 ..        6 atom(s) in final combined selection 

 .. number of atoms in PDB:  6 
 .. number of calphas in PDB: 6 
 .. number of residues in PDB: 6 


pdb


 Call:  as.pdb.default(resno = 1:6, resid = "ALA", elety = "CA", chain = "A")

   Total Models#: 1
     Total Atoms#: 6,  XYZs#: 18  Chains#: 1  (values: A)

     Protein Atoms#: 6  (residues/Calpha atoms#: 6)
     Nucleic acid Atoms#: 0  (residues/phosphate atoms#: 0)

     Non-protein/nucleic Atoms#: 0  (residues: 0)
     Non-protein/nucleic resid values: [ none ]

   Protein sequence:
      AAAAAA

+ attr: atom, xyz, calpha, call


## Read a PRMTOP file
prmtop <- read.prmtop(system.file("examples/crambin.prmtop", package="bio3d"))

## Read Amber coordinates
crds <- read.crd(system.file("examples/crambin.inpcrd", package="bio3d"))

## Atom selection
ca.inds <- atom.select(prmtop, "calpha")


 Summary of PDB generation:
 .. number of atoms in PDB determined by 'xyz'

 ..       46 atom(s) from 'string' selection 
 ..       46 atom(s) in final combined selection 

 .. number of atoms in PDB:  639 
 .. number of calphas in PDB: 46 
 .. number of residues in PDB: 46 



## Convert to PDB format
pdb <- as.pdb(prmtop, crds, inds=ca.inds)


 Summary of PDB generation:
 .. number of atoms in PDB determined by 'xyz'

 ..       46 atom(s) from 'string' selection 
 ..       46 atom(s) in final combined selection 

 .. number of atoms in PDB:  46 
 .. number of calphas in PDB: 46 
 .. number of residues in PDB: 46 




## Read a single entry MOL2 file
## (returns a single object)
mol <- read.mol2( system.file("examples/aspirin.mol2", package="bio3d") )

## Convert to PDB
pdb <- as.pdb(mol)


 Summary of PDB generation:
 .. number of atoms in PDB determined by 'xyz'

 ..        0 atom(s) from 'string' selection 
 ..        0 atom(s) in final combined selection 

 .. number of atoms in PDB:  20 
 .. number of calphas in PDB: 0 
 .. number of residues in PDB: 1

Author

Lars Skjaerven

Convert to PDB format