atomiccomp

Вычислите атомарный состав белка

Синтаксис

NumberAtoms = atomiccomp(SeqAA)

Аргументы

SeqAA

Аминокислотная последовательность. Введите вектор символов или вектор целых чисел из таблицы Сопоставление кодов Буквы аминокислот с целыми числами. Вы также можете ввести структуру с полем Sequence.

Описание

NumberAtoms = atomiccomp(SeqAA) подсчитывает тип и количество атомов в аминокислотной последовательности (SeqAA) и возвращает счетчики в 1-by- 1 структура (NumberAtoms) с полями C, H, N, O, и S.

Примеры

свернуть все

Извлеките аминокислотную последовательность из базы данных NCBI GenPept.

rhodopsin = getgenpept('NP_000530')
rhodopsin = 

  struct with fields:

                LocusName: 'NP_000530'
      LocusSequenceLength: '348'
     LocusNumberofStrands: ''
            LocusTopology: 'linear'
        LocusMoleculeType: ''
     LocusGenBankDivision: 'PRI'
    LocusModificationDate: '07-AUG-2015'
               Definition: 'rhodopsin [Homo sapiens].'
                Accession: 'NP_000530'
                  Version: 'NP_000530.1'
                       GI: '4506527'
                  Project: []
                   DBLink: ' DBSOURCE    REFSEQ: accession NM_000539.3'
                 Keywords: 'RefSeq.'
                  Segment: []
                   Source: 'Homo sapiens  human '
           SourceOrganism: [4×65 char]
                Reference: {1×18 cell}
                  Comment: [26×67 char]
                 Features: [219×74 char]
                      CDS: [1×1 struct]
                 Sequence: 'mngtegpnfyvpfsnatgvvrspfeypqyylaepwqfsmlaaymfl...'
                SearchURL: 'http://www.ncbi.nlm.nih.gov/entrez/viewer.fcgi...'
              RetrieveURL: 'http://eutils.ncbi.nlm.nih.gov/entrez/eutils/e...'

Подсчитайте атомы в последовательности.

rhodopsinAC = atomiccomp(rhodopsin)
rhodopsinAC = 

  struct with fields:

    C: 1814
    H: 2725
    N: 423
    O: 477
    S: 25

Найдите количество атомов углерода.

rhodopsinAC.C
ans =

        1814

См. также

| |

Представлено до R2006a