Анализ результатов моделирование перехода А-формы ДНК в В-форму в воде.

Начните просмотр результатов с изучения файла mdrun_mpi.out...

Занесите в отчет информацию о времени моделирования, продолжительности моделирования и количество использованых процессоров.

1. Любой анализ начинают с визуального анализа движений молекул. При вопросе о выводк групп выберите DNA.

      trjconv -f b_md.xtc -s b_md.tpr -o b_pbc_1.pdb -skip 20 -pbc mol
      

   Откройте ваш b_pbc_1.pdb в PyMol , не забудьте включить анимацию. Если Вас не устроил результат визуализации попробуйте:

      trjconv -f b_md.xtc -s b_md.tpr -o b_fit_1.pdb -skip 20 -fit rot+trans
      

   Занесите ваши наблюдения в журнал. Отметье в какой момент происходит образование B-формы. Для опредления соотвествия между номером модели и временем моделирования найдите в pdb файле выражение "MODEL 50" (цифра 50 это и есть номер модели) и двумя строчками выше будет упосинание о времени в моедлировании.
2. Определите средне-квадратичное отколнение в ходе моделирования. Так как у нас происходит конформационный переход сначала расчитаем отклонение в ходе все симуляции относительно стартовой структуры.

      g_rms -f dna_md.xtc -s dna_md.tpr -o rms_1
      

   И относительно каждой предидущей структуры на растоянии 400 кадров. Если ближе к концу закончился конформационный переход, то отколнение должно уменьшаться.

      g_rms -f dna_md.xtc -s dna_md.tpr -o rms_2 -prev 400
      

3. Определите изменение гидрофобной и гидрофильной поверхности в ходе конформационного перехода.

      g_sas -f dna_md.xtc -s dna_md.tpr -o sas_dna.xvg
      

   Постройте в Excell зависимость изменения гидрофобной гидрофильной поверхностей доступных растворителю от времени. Сделайте вывод о возможных причинах конформационного перехода из В в А.
4. Традиционным анализом для ДНК является расчет колчества образуемых водородных связей. Если мы будем исследовать связи между ДНК и ДНК, то это будут водородные связи между цепями ДНК. Для конца траектории:

      g_hbond  -f dna_md.xtc -s dna_md.tpr -num hbond_dna
      

   Не менее интересно будет изучить количество вдородных связей ДНК-Вода

      g_hbond  -f dna_md.xtc -s dna_md.tpr -num hbond_dna_sol 
      

   Постройте зависимости в Excel и добавтье к отчету. Сделайте вывод об количестве водородных связей в этом дуплексе, соотвествует ли это каноническим представлениям? Как влияет конформационный переход на количество водородных связей с водой?. делайте вывод о возможных причинах конформационного перехода из В в А.

-- AndreyGolovin - 2011-04-07