Документ взят из кэша поисковой машины. Адрес оригинального документа : http://kodomo.cmm.msu.su/FBB/year_02/term_6/block_5/md_task1.doc Дата изменения: Fri Apr 22 10:35:32 2005 Дата индексирования: Tue Oct 2 00:28:46 2012 Кодировка: koi8-r

Задание 1
Задача первого занятия: провести симуляцию молекулярной динамики
предложенного белка, в свободном состоянии используя программный пакет
Gromacs. Всё подготовительные операции вы будете проводить на сервере
kodomo-count.cmm.msu.su

1) Загрузите предложенный вам файл из банка третичных структур PDB
(www.pdb.org) . Используя FTP помесите полуученый вами файл в рабочую
поддиректорию free.
2) Используя любой текстовый редактор создайте pdb файл содержащий записи
ATOM только Вашего белка.
3) Используя программу pdb2gmx получите файл топологии (*top) и файл
координат (*gro) вашего белка используя полученный в пункте 2. pdb
файл. Рекомендуемый синтаксис запуска программ предложен в подсказках.
4) Воспользуйтесь программой editconf для изменения параметров Вашей
ячейки. Установите отступ от края белка в 10 А.
5) Поместите файл параметров МД для минимизации энергии (em.mdp) в
рабочую директорию. Проведите минимизацию энергии структуры белка,
используя препроцессор grompp и саму программу процессор mdrun.
Установите и занесите в файл протокола изменение потенциальной энергии
в ходе минимизации энергии структуры.
6) Поместите в ячейку молекулы растворителя (программа genbox).
Используйте модель воды spc.
7) Если общий заряд системы не равен 0 (см. выдачу grompp п.5), то
нейтрализуйте общий заряд системы. Для этого, используя файл em.mdp
(пункт 5) с помощью программы grompp создайте файл бинарной топологии
(tpr). Используйте полученный файл как входной файл для программы
genion. Эта программа заменяет молекулы воды на ионы. Количество и тип
добавляемых ионов регулируется флагами -np или -nn.
8) Поправьте количество молекул воды в файле топологии (секция [
molecules ] ) и добавьте запись о ионах.
9) Проведите минимизацию энергии полученной системы (см. пункт 5).
10) Проведите «утряску воды». Для этого скопируйте в рабочую директорию
файл pr.mdp и отредактируйте его в секции контроля температуры. Вам
необходимо, что бы имя иона совпадало с тем, который вы использовали
при нейтрализации заряда. Получите файл бинарной топологии используя
этот файл и запустите mdrun.
11) Скопируйте в рабочую директорию файл md.mdp, и отредактируйте его в
секции контроля температуры. Получите файл бинарной топологии,
используя этот файл. Внимательно изучите файл mdout.mdp и если вы
согласны с параметрами, то запустите саму молекулярную динамику.



Подсказки.

Все программы пакета GROMACS имеют общую идеологию: на вход вы подаёте
входные файлы и на выходе получаете результат работы программы. Если вы
не указали, какую либо обязательную опцию, то программа выдаст сообщение
об ошибке. Если эта опция была не обязательная, то программа использует
значения по умолчанию. ВНИМАНИЕ! В большинстве случаев информации в
сообщениях об ошибках достаточно для того, что бы найти и исправить эту
ошибку.
Узнать опции каждой программы в GROMACS, можно запустив эту программу с
флагом -h (editconf -h).
Файлы параметров МД (*.mdp) вы можете найти в поддиректории md_files.

Рекомендуемый синтаксис запуска программ в пунктах задания.
Используя этот синтаксис Вы можете заменить ${i} на ID вашего белка или
присвоить значение переменной i используя команду : export i="mi_ID", где
my_ID это ID вашего белка.
2) Удобно убрать ненужную информацию из pdb файла с помощью grep.
grep "^ATOM" pdb${i}.ent > ${i}.pdb
3) Использование pdb2gmx.
pdb2gmx -f ${i}.pdb -o ${i} -p ${i} -ff gmx
4) Использование editconf.
editconf -f ${i} -o ${i} -d 1
5) Использование grompp и mdrun.
grompp -c ${i} -p ${i} -o ${i}_em1 -f ./em -v
mdrun -s ${i}_em1 -o ${i}_em1 -c ${i}_aem1 -v
6) Использование genbox.
genbox -cp ${i}_aem1 -o ${i}_s -p ${i} -cs spc216
7) Использование genion.
grompp -c ${i}_s -p ${i} -o ${i}_s -f ./em -v
genion -s ${i}_s -o ${i}_si -np X (or -nn X)
8) Отредактируйте файл ${i}.top (vi или mcedit)Поправьте количество
молекул воды в секции [ molecules ] и добавьте запись о ионах. Пример:
было стало
|[ molecules ] |[ molecules ] |
|; Compound #mols |; Compound #mols |
|Protein_A 1 |Protein_A 1 |
|SOL 15750 |SOL 15747 |
| |Na 3 |


9) Минимизация энергия соляного раствора.
grompp -c ${i}_si -p ${i} -o ${i}_em2 -f ./em -v
mdrun -s ${i}_em2 -o ${i}_em2 -c ${i}_aem2 -v
10) Исправление pr.mdp :
было стало
|; Berendsen temperature |; Berendsen temperature |
|coupling |coupling |
|Tcoupl = berendsen |Tcoupl = berendsen |
|tau_t = 0.1 0.1 0.1 |tau_t = 0.1 0.1 0.1 |
|tc-grps = SOL Protein |tc-grps = SOL Protein |
|Na |Cl |
|ref_t = 300 300 300 |ref_t = 300 300 300 |


«Утряска воды»:
grompp -c ${i}_si -p ${i} -o ${i}_pr -f ./pr -v
mdrun -s ${i}_pr -o ${i}_pr -c ${i}_apr -v
11) Молекулярная динамика на кластере SCI:
grompp -c ${i}_apr -p ${i} -o ${i}_md -f ./md -v
mpirun -np 1 -as single /home/golovin/progs/gmx-sci/i686-pc-linux-
gnu/bin/mdrun -s ${i}_md.tpr -o ${i}_md -c ${i}_amd -v