atomiccomp - Calculate atomic composition of protein

Syntax

Arguments

SeqAA

Amino acid sequence. Enter a character string or vector of integers from the table Mapping Amino Acid Letter Codes to Integers. You can also enter a structure with the field Sequence.

Description

NumberAtoms = 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 (NumberAtoms) with fields C, H, N, O, and S.

Examples

1. Retrieve an amino acid sequence from the NCBI GenPept database.

   rhodopsin = getgenpept('NP_000530');

2. Count the atoms in the sequence.

   rhodopsinAC = atomiccomp(rhodopsin)
   
   rhodopsinAC = 
   
       C: 1814
       H: 2725
       N: 423
       O: 477
       S: 25
   

3. Count the number of carbon atoms in the sequence.

   rhodopsinAC.C
   
   ans =
   
           1814

See Also

aacount | molweight | proteinplot

© 1984-2012- The MathWorks, Inc.    -   Site Help   -   Patents   -   Trademarks   -   Privacy Policy   -   Preventing Piracy   -   RSS