aa2mass: Amino Acid Residues to Mass Converter

Description

Convert a sequence of amino acid residue names to mass.

Usage

aa2mass(pdb, inds=NULL, mass.custom=NULL, addter=TRUE, mmtk=FALSE)

Arguments

pdb

a character vector containing the atom names to convert to atomic masses. Alternatively, a object of type pdb can be provided.

inds

atom and xyz coordinate indices obtained from atom.select that selects the elements of pdb upon which the calculation should be based.

mass.custom

a list of amino acid residue names and their corresponding masses.

addter

logical, if TRUE terminal atoms are added to final masses.

mmtk

logical, if TRUE use the exact aminoacid residue masses as provided with the MMTK database (for testing purposes).

Value

Returns a numeric vector of masses.

Details

This function converts amino acid residue names to their corresponding masses. In the case of a non-standard amino acid residue name mass.custom can be used to map the residue to the correct mass. User-defined amino acid masses (with argument mass.custom) will override mass entries obtained from the database.

See examples for more details.

References

Grant, B.J. et al. (2006) Bioinformatics 22, 2695--2696.

Examples

Run this code

resi.names <- c("LYS", "ALA", "CYS", "HIS")
masses <- aa2mass(resi.names, addter=FALSE)

## Fetch atomic masses in a PDB object
pdb <- read.pdb("3dnd")
masses <- aa2mass(pdb)

## or
masses <- aa2mass(pdb$atom[1:10,"resid"])

## Dealing with unconventional residues
#pdb <- read.pdb("1xj0")

#mass.cust <- list("CSX"=122.166)
#masses <- aa2mass(pdb, mass.custom=mass.cust)

Run the code above in your browser using DataLab

Data Engineering and BI courses are free this week!