bio3d (version 2.3-0)

read.prmtop: Read AMBER Parameter/Topology files

Description

Read parameter and topology data from an AMBER PrmTop file.

Usage

read.prmtop(file)
"print"(x, printseq=TRUE, ...)

Arguments

file

a single element character vector containing the name of the PRMTOP file to be read.

x

a PRMTOP structure object obtained from read.prmtop.

printseq

logical, if TRUE the residue sequence will be printed to the screen. See also pdbseq.

...

additional arguments to ‘print’.

Value

Returns a list of class ‘prmtop’ (inherits class ‘amber’) with components according to the flags present in the PrmTop file. See the AMBER documentation for a complete list of flags/components: http://ambermd.org/formats.html.Selected components:

Details

This function provides basic functionality to read and parse a AMBER PrmTop file. The resulting ‘prmtop’ object contains a complete list object of the information stored in the PrmTop file.

See examples for further details.

References

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

See Also

read.crd, read.ncdf, as.pdb, atom.select, read.pdb

Examples

## Not run: 
# ## Read a PRMTOP file
# prmtop <- read.prmtop(system.file("examples/crambin.prmtop", package="bio3d"))
# print(prmtop)
# 
# ## Explore prmtop file
# head(prmtop$MASS)
# head(prmtop$ATOM_NAME)
# 
# ## Read Amber coordinates
# crds <- read.crd(system.file("examples/crambin.inpcrd", package="bio3d"))
# 
# ## Atom selection
# ca.inds <- atom.select(prmtop, "calpha")
# 
# ## Convert to PDB format
# pdb <- as.pdb(prmtop, crds)
# pdb.ca <- as.pdb(prmtop, crds, inds=ca.inds)
# 
# ## Trajectory processing
# #trj <- read.ncdf("traj.nc", at.sel=ca.inds)
# 
# ## Convert to multimodel PDB format
# #pdb <- as.pdb(prmtop, trj[1:20,], inds=ca.inds, inds.crd=NULL)
# 
# ## RMSD of trajectory
# #rd <- rmsd(crds$xyz[ca.inds$xyz], traj, fit=TRUE)
# ## End(Not run)