Имея pdb-файл записи 1jw9 (1jw9.pdb), запустив на kodomo команду:
dsspcmbi 1jw9.pdb 1jw9.dssp
получим на выходе dssp (1jw9.dssp).
Выходной файл импортирую в excel (dssp.xls). Необходимая нам информация заключена всего в двух колонках: RESIDUE (где приведены номер остатка и
идентификатор цепи) и STRUCTURE, из содержимого каждой из ячеек которой нам нужен только первый символ. Используя данные выдачи dssp, получу информацию
об элементах вторичной структуры для каждой цепи.
Цепь B:
1 альфа-спираль: 6-11
2 альфа спираль: 13-16
3 альфа спираль: 22-31
4 альфа спираль: 41-53
5 альфа спираль: 86-97
6 альфа спираль: 112-120
7 альфа спираль: 132-145
8 альфа спираль: 174-177
9 альфа спираль: 194-213
1 бета-тяж: 33-37
2 бета-тяж: 57-61
3 бета-тяж: 102-106
4 бета-тяж: 124-127
5 бета-тяж: 149-155
6 бета-тяж: 158-164
7 бета-тяж: 222-227
8 бета-тяж: 232-237
Цепь D:
1 альфа-спираль: 8-14
2 альфа спираль: 28-36
3 альфа спираль: 40-45
1 бета-тяж: 2-6
2 бета-тяж: 18-21
3 бета-тяж: 52-55
4 бета-тяж: 58-59
5 бета-тяж: 71-75
Теперь посмотрим данные о вторичной структуре, приведенные авторами в pdb-файле.
Цепь B:
1 альфа-спираль: 5-12
2 альфа спираль: 12-17
3 альфа спираль: 21-32
4 альфа спираль: 40-54
5 альфа спираль: 67-73
6 альфа спираль: 78-82
7 альфа спираль: 85-98
8 альфа спираль: 111-121
9 альфа спираль: 131-146
10 альфа спираль: 172-178
11 альфа спираль: 193-214
1 бета-тяж: 33-37
2 бета-тяж: 57-61
3 бета-тяж: 102-106
4 бета-тяж: 124-127
5 бета-тяж: 149-155
6 бета-тяж: 158-164
7 бета-тяж: 222-227
8 бета-тяж: 232-237
Цепь D:
1 альфа-спираль: 7-15
2 альфа спираль: 27-37
3 альфа спираль: 39-46
1 бета-тяж: 2-6
2 бета-тяж: 18-21
3 бета-тяж: 52-55
4 бета-тяж: 58-59
5 бета-тяж: 71-75
Теперь сравним результаты. DSSP не нашел две альфа-спирали цепи B, которые числятся в pdb (выделены жирным выше - 67-73 и 78-82). В выдаче dssp эти две последовательности остатков считаются G = 3-helix (3/10 helix) . Все остальные элементы вторичной структуры полностью совпали, если учесть, что в dssp в альфа-спираль не включены первый и последний ее остатки (таков алгоритм программы).