BioPhysConnectoR (version 1.6-10)

aa2num: Conversion of Amino Acids into Integer Values

Description

A sequence containing standard amino acids is converted into a sequence of integer values. An offset for the code can be specified.

Usage

aa2num(seq, offset = 0, code = 0:19, verbose = FALSE)

Arguments

seq

string vector containing a sequence of standard amino acids

offset

offset, added to the specified code

code

integer vector of the values to be assigned to the amino acids

verbose

logical, if TRUE the specific encoding is printed

Value

Returns a sequence of integer values according to the given numerical code.

Details

The default values for the code are the integer values in the range 0 to 19. It is possible to use another numerical code for the 20 amino acids. The default coding order of the amino acids is: CYS, MET, PHE, ILE, LEU, VAL, TRP, TYR, ALA, GLY, THR, SER, ASN, GLN, ASP, GLU, HIS, ARG, LYS, PRO. Each amino acid is converted to a single element according to the code-vector. The offset offset is added. Both sequences with three-letter and one-letter code can be handled.

References

Grant, Rodrigues, ElSawy, McCammon, Caves, (2006) Bioinformatics 22, 2695--2696.

See Also

show.code

Examples

seq<-c("MET", "GLY", "PRO", "LYS", "ASN")
aa2num(seq)