ShadTask12

Version 4 (Andrey Golovin, 04.12.2013 18:12)

1 1 Andrey Golovin
2 1 Andrey Golovin
h1. Макромолекулярный докинг
3 1 Andrey Golovin
4 1 Andrey Golovin
Суть задания ознакомиться с программой ZDOCK и сделать предсказание о свзывании фрагмента антитела (VHH domain) c панкреатической альфа-амилазой.
5 1 Andrey Golovin
6 1 Andrey Golovin
7 1 Andrey Golovin
# Вся работа будет происходить на shadbox (cм. предыдущие практикумы).
8 1 Andrey Golovin
# Скопируем в рабочую директорию файлы amilase.pdb,camelid.pdb и uniCHARMM.
9 1 Andrey Golovin
<pre>
10 2 Andrey Golovin
cp ~/materials/zdock/* .
11 1 Andrey Golovin
</pre>
12 1 Andrey Golovin
# Обратите внимание, что в файле amilase.pdb есть некоторые водороды. Давайте добавим водороды к обеим pdb исполльзуя pdb2gmx. 
13 1 Andrey Golovin
Пример: 
14 1 Andrey Golovin
<pre>pdb2gmx -f мой.pdb -o мой_h.pdb -p</pre>
15 1 Andrey Golovin
Используйте силовое поле Gromos. Опция *-ignh* позволит избежать ошибки об *unknown atoms*.
16 1 Andrey Golovin
# С помощью утилиты *mark_sur* проведём препроцессинг файлов pdb. 
17 1 Andrey Golovin
Пример: 
18 1 Andrey Golovin
<pre>mark_sur my.pdb my_m.pdb</pre>
19 1 Andrey Golovin
# Постарайтесь разобраться в том, что делает утилита *mark_sur*, внесите ваше предположение в отчёт.
20 1 Andrey Golovin
# Прочитайте информацию об использовании *zdock* просто запустив его.
21 1 Andrey Golovin
 + Для успешной работы zrank и zdock из pdb файлов надо удалить строчки MODEL и TER +
22 1 Andrey Golovin
# Запустите сам процесс докинга.
23 1 Andrey Golovin
# Проведите предварительный анализ результатов докинга с помощью утилиты zrank: 
24 3 Andrey Golovin
25 4 Andrey Golovin
>Забавно, но для работы zrank  НАДО удалить третью строчку в zdock.out где три значения и все равны нулю
26 3 Andrey Golovin
27 1 Andrey Golovin
<pre>zrank zdock.out 1 100</pre>
28 1 Andrey Golovin
     или
29 1 Andrey Golovin
<pre>create.pl zdock.out
30 1 Andrey Golovin
ls complex*pdb >list
31 1 Andrey Golovin
zrank list
32 1 Andrey Golovin
</pre>
33 1 Andrey Golovin
# Визуально сравните результаты с известной структурой 1kxt.
34 1 Andrey Golovin
# * Если Вас результат не устроил, попробуйте опцию **-D** программы **zdock**
35 1 Andrey Golovin
# * Попробуйте найти в ваших  результатах структуру наиболее близкую к результату РСА. Используйте скрипт на bash  с использованием g_rms.