# NOT RUN {
## Read Amber PRMTOP and CRD files
prm <- read.prmtop(system.file("examples/crambin.prmtop", package="bio3d"))
crd <- read.crd(system.file("examples/crambin.inpcrd", package="bio3d"))
## Convert to PDB format
pdb <- as.pdb(prm, crd)
## Atom selection
ca.inds <- atom.select(prm, "calpha")
# }
Run the code above in your browser using DataLab