Description
Returns a table with the amino acid composition of the target protein
Usage
aa.comp(target, uniprot = TRUE)
Arguments
target
a character string specifying the UniProt ID of the protein of interest or, alternatively, the sequence of that protein.
uniprot
logical, if TRUE the argument ‘target’ should be an ID.
Value
Returns a dataframe with the absolute frequency of each type of residue found in the target peptide.
See Also
aa.at(), is.at(), renum.pdb(), renum.meto(), renum()
Details
The ptm package offers a set of ancillary functions aimed to carry out rutinary work, which may be needed when more elaborated analysis are required. Among these ancillary function are:
- aa.at
- is.at
- aa.comp (the current document)
- renum.pdb
- renum.meto
- renum
The function aa.comp() provides us the amino acid composition of a given sequence.
library(knitr)
mytable <- aa.comp('P01009')
mytable$percent <- round(100 * mytable$frequency/sum(mytable$frequency), 1)
kable(mytable)
aa | frequency | percent |
---|---|---|
A | 26 | 6.2 |
R | 7 | 1.7 |
N | 19 | 4.5 |
D | 24 | 5.7 |
C | 3 | 0.7 |
Q | 18 | 4.3 |
E | 32 | 7.7 |
G | 24 | 5.7 |
H | 13 | 3.1 |
I | 20 | 4.8 |
L | 51 | 12.2 |
K | 34 | 8.1 |
M | 10 | 2.4 |
F | 27 | 6.5 |
P | 19 | 4.5 |
S | 25 | 6.0 |
T | 30 | 7.2 |
W | 3 | 0.7 |
Y | 6 | 1.4 |
V | 27 | 6.5 |
Above we have passed the UniProt ID of the protein of interest, but we can also provide any customized sequence:
library(knitr)
mytable <- aa.comp('MAARESAPFDR', uniprot = FALSE)
mytable$percent <- round(100 * mytable$frequency/sum(mytable$frequency), 1)
kable(mytable)
aa | frequency | percent |
---|---|---|
A | 3 | 27.3 |
R | 2 | 18.2 |
N | 0 | 0.0 |
D | 1 | 9.1 |
C | 0 | 0.0 |
Q | 0 | 0.0 |
E | 1 | 9.1 |
G | 0 | 0.0 |
H | 0 | 0.0 |
I | 0 | 0.0 |
L | 0 | 0.0 |
K | 0 | 0.0 |
M | 1 | 9.1 |
F | 1 | 9.1 |
P | 1 | 9.1 |
S | 1 | 9.1 |
T | 0 | 0.0 |
W | 0 | 0.0 |
Y | 0 | 0.0 |
V | 0 | 0.0 |