Bioinformatics ToolboxPrevious page   Next page

atomiccomp

Calculate atomic composition of a protein

Syntax 

Atoms = atomiccomp(SeqAA)

Arguments

SeqAA

Amino acid sequence. Enter a character string or vector of integers from the table . You can also enter a structure with the field Sequence.

Description

Atoms = atomiccomp(SeqAA) counts the type and number of atoms in an amino acid sequence (SeqAA) and returns the counts in a 1-by-1 structure (Atoms) with fields C, H, N, O, and S.

Examples

Get an amino acid sequence from the Protein Sequence Database (PIR-PSD) and count the atoms in the sequence.

pirdata = getpir('cchu','SequenceOnly',true);
mwcchu = atomiccomp(pirdata)

mwcchu = 
    C: 526
    H: 845
    N: 143
    O: 149
    S: 6

mwcchu.C

ans=
      526

See Also

Bioinformatics Toolbox functions aacount, msalign

Previous page aminolookup basecountNext page

© 1994-2005 The MathWorks, Inc.