Описание для программы g_hbond, которая производит компьютерный анализ водородных связей.
g_hbond - компьютерный анализ водородных связей /* для System, Water, SOL - вообщем для индексных групп, отвечающих за целы молекулы в-ва (готовые соединения).*/ OH и NH группы относятся к донорам, О всегда акцптор, N - акцептор по умолчанию, но это можно изменить исользуя -nitacc. Фиктивные водородные атомы предпологаются присоединееными к первому предшествующему неводородному атому.
Если вы установили -shell, вам зададут вопрос о добавлении индексных групп, которые должны вмещать ровно один атом. В таком случае рассматриваются только водородные связи между атомами в пределах shell дистанции из одного атома.
Также возможен анализ особых водородных связей с -sel. Это требует индексный файл, который болжен включать группы атомов триплетов Донор Водород Акцептор в следующем порядке:
[ selected ]
20 21 24
25 26 29
1 3 6
Замечу, что нет необходимости в разделительных линиях между триплетами. Каждый атом триплета определенной водородной связи будет проанализирован, замечу также что не проверяются типы атомов.
/* -ins меняет направление расчетного растворителя вставленного в водородные связи. */ В этом случае в добавленные группы должны быть отобранны определенные молекулы растворителя.
Водородная связь существует, тогда когда выполняется геометрический критерий водородной связи r <= rHB = 0,35 nm (значение первого минимума на ФРР) и угол a <= aHB = 30 градусов. Программа анализирует все водородные связи, существующие между двумя группами, которые должны быть или идентичными, или не пересекающимися /*неперекрывающимися*/ или особенными Донор-Водород-Аксцептор триплетами по следующему алгоритму:
-num: число водородных связей как функция от времени
-ac: среднее по всем автокорреляциям существующих функций (или 0 или 1) всех водородных связей.
-dist: распределение по расстоянию (r) Донор-Акцептор для всех водородных связей
-ang: распределение по углу Донор-Водород-Акцептор а для всех водородных связей
-hx: число водородных связей как функция времени, разделенных на группы n-n+i, где n и n+i означают номера молекул и i принадлежит от 0 до 6. Группа с i=6 влключает все водородные связи для i>6. Эти группы включают n-n+3, n-n+4 и n-n+5 водородных связей связанных со спиралями в протеинов.
-life: время жизни водороных связей вычисляется для усредненных по все автокорреляционым ф-циям существующих ф-ций (или 0 или 1) всех водородных связий: С(тау)=<Si(t)*Si(t+тау)>, где Si(t) = {0,1} для i водородной связи в момент времени t. Интеграл тауHB от С(тау) берется с грубой оценкой для среднего времени жизни водородной связи от 0 до бесконечности по d(тау). Оба интеграл и полная автокорреляционная ф-ция С(тау) будет выходными, так что можно будет проводить более сложный анализ (например, с использованием мульти-экспоненты) для получения более точных оценок для времени жищзни тауHB. Более полный анализ в работе. [135], один из более хороших обзоров у Luzar и Чендлер - анализ кинетики водородных связей [136, 137].
-hbm: карта жизни водородных связей будет сгенерирована на измеренях #H-bonds x #frames. Упорядочение такой же как в индексом файле (см. ниже), но вспять, это означает, что последние триплет в индексном файле соответствует первой строке существования на карте.
Порядок идентичный тому, что и в -hbn индексном файле.
-hbn: индексные группы выходя содержат проанализированные группы, все Донор-Водородные атомные пары и Акцепторные пары в этих группах, Донор-Водород-Акцептор триплеты в водородных связях между анализируемой группой и всеми растворенными атомами участвующих в вставке.
-dan: запись выходных номеров доноров и акцепторов проанализированных для каждого временного отрезка. Это особенно полезно, когда используется -shell
-nhbdist: вычисление числа HBonds от водородов в порядке сравнения результатов с спекстроскопией Рамана
Замечание: опции -ac, -life, -hbn и -hbm требуют наличия памяти пропорциональной общему числу донерного времени, общему числу акцепторного времени выбранной(ых) группы(групп).
Файлы:
-f traj.xtc Input Trajectory: xtc trr trj gro g96 pdb cpt
-s topol.tpr Input Run input file: tpr tpb tpa
-n index.ndx Input, Opt. Index file
-num hbnum.xvg Output xvgr/xmgr file // /* выходной 4 */файл с графиком колличества водородных связей за отрезок прогона
-g hbond.log Output, Opt. Log file // /* выходной 3 */ лог файл содержаший в себе три столбца номеров атомов: донер, водород, акцептор.
-ac hbac.xvg Output, Opt. xvgr/xmgr file // /* выходной 8 */ файл с графиком автокореляционной ф-ции водородных связей /* строится ACF число/число */
-dist hbdist.xvg Output, Opt. xvgr/xmgr file // /* выходной 5 */ файл с графиком распределения водородных связей (Водород-Акцептор (расстояние))
-ang hbang.xvg Output, Opt. xvgr/xmgr file // /* выходной 6 */ файл с графиком распределения водородных связей (Донор-Водород-Акцептор (угол))
-hx hbhelix.xvg Output, Opt. xvgr/xmgr file // /* выходной 7 */ 0 файл с графиком водородных связей за отрезок прогона
-hbn hbond.ndx Output, Opt. Index file // /* выходной 2 */ индексный файл содержаший индексные группы: выбранная группа при обработке файла траекторий, водородные доноры этой группы, акцепторы этой группы, список водородных связей
-hbm hbmap.xpm Output, Opt. X PixMap compatible matrix file // /* выходной 10 */ выходной файл с картой жизни водородных связей
-don donor.xvg Output, Opt. xvgr/xmgr file // /* выходной 1 */ файл с графиком числа доноров за отрезок прогона (свойства доноров)
-dan danum.xvg Output, Opt. xvgr/xmgr file // /* выходной 11 */ файл с графиком соотношения доноров и акцепторов
-life hblife.xvg Output, Opt. xvgr/xmgr file // /* выходной 9 */ файл с графиком времени жизни неразрущенных водородных связей /* строится процентами % */
-nhbdist nhbdist.xvg Output, Opt. xvgr/xmgr file
Дополнительные опции:
-[no]h bool no
-nice int 19
-b time 0 //первый фрагмент (отрезок) пс для чтения директории
-e time 0 //последний фрагмент (отрезок) пс для чтения директории
-dt time 0
-[no]ins bool no // анализ растворителя вставки
-a real 30 // размер угла (градусы, Акцептор - Донор - Водород), угол OOH
-r real 0.35 // радиус обрезания (нм, X - Акцептор, см.след.опцию), например растояние Roh
-[no]da bool yes // использовать растояние Донор-Акцептор (если истина) или Водород-Акцептор (если ложь)
-r2 real 0 // второй радиус обрезания. Главным образом полезет с -contact и -ac
-abin real 1 // распределение по углам Бинвинда (градусы)
-rbin real 0.005 // распределение по углам Бинвинда (нм)
-[no]nitacc bool yes // расценивать азотные атомы как акцепторы
-[no]contact bool no // не смотреть водородные связи, но просто связывать в пределах радиуса обрезания
-shell real -1 // когда > 0, вычисляются воодородные связи в пределах # нм shell около (вокруг) одной частицы
-fitstart real 1 // Время (ps) с когорого начинает работать подгоночная ф-ция для ACF в соответствии с передовой и обратной (остальной) оценочной постоянной для разрыва и формирования HB.
-temp real 298.15 // температура (K) для вычисления энергии Гиббса соответствующая разрыву и преобразованию HB
-smooth real -1 // Если >= 0, конец ACF будет сглажен выбранной подгоночной экспоненциальной функцией: y = A exp(-x/tau)
-dump int 0 // выбросить первые N водородных связей ACFs в одиночном /*единственном*/ xvg file для отладки
-max_hb real 0 // Теоретический максимум числа водородных связей использумемых для нормализации HB автокорреляционной ф-ции. Может быть использован в случае, если программа оценивает не правильно
-[no]merge bool yes // H-bonds между теми же донором и акцептором, но с другим водородом рассматриваемым как одиночая водородная связь. Главным образом важно для ACF.
-acflen int -1 // длинна для ACF, по умолчанию половина числа фрагментов.
-[no]normalize bool yes // норамализовать ACF
-P enum 0 // порядок полинома Лежандра для автокорреляционной ф-ции (0 указывает на none): 0,1,2 или 3
-fitfn enum none // подгоночная функция: none, exp, aexp, exp_exp, vac, exp5, exp7 or exp9
-ncskip int 0 // пропустить N точек в выходном файле автокорреляционной ф-ции
-beginfit real 0 // время, для начала подгоночной экспоненциальной ф-ции для CF
-endfit real -1 // время, для конца подгоночной экспоненциальной ф-ции для CF, -1 пока не достигнет конца.
========
Также рекомендую внимательно смотреть на консольный вывод программы, из него можно тоже под черпнуть важные вещи, как число водородных связей, общее число пар атомов, данные по ACF, а также ссылки на статьи по различным аспектам связанным с этой программой.
P.S. Если статья была вам полезна проголосуйте в опросе блога, этим вы скажите СПАСИБО Автору, также можно обсудить её на Форуме.
Назад к оглавлению : руководство по Gromacs
g_hbond - компьютерный анализ водородных связей /* для System, Water, SOL - вообщем для индексных групп, отвечающих за целы молекулы в-ва (готовые соединения).*/ OH и NH группы относятся к донорам, О всегда акцптор, N - акцептор по умолчанию, но это можно изменить исользуя -nitacc. Фиктивные водородные атомы предпологаются присоединееными к первому предшествующему неводородному атому.
Если вы установили -shell, вам зададут вопрос о добавлении индексных групп, которые должны вмещать ровно один атом. В таком случае рассматриваются только водородные связи между атомами в пределах shell дистанции из одного атома.
Также возможен анализ особых водородных связей с -sel. Это требует индексный файл, который болжен включать группы атомов триплетов Донор Водород Акцептор в следующем порядке:
[ selected ]
20 21 24
25 26 29
1 3 6
Замечу, что нет необходимости в разделительных линиях между триплетами. Каждый атом триплета определенной водородной связи будет проанализирован, замечу также что не проверяются типы атомов.
/* -ins меняет направление расчетного растворителя вставленного в водородные связи. */ В этом случае в добавленные группы должны быть отобранны определенные молекулы растворителя.
Водородная связь существует, тогда когда выполняется геометрический критерий водородной связи r <= rHB = 0,35 nm (значение первого минимума на ФРР) и угол a <= aHB = 30 градусов. Программа анализирует все водородные связи, существующие между двумя группами, которые должны быть или идентичными, или не пересекающимися /*неперекрывающимися*/ или особенными Донор-Водород-Аксцептор триплетами по следующему алгоритму:
-num: число водородных связей как функция от времени
-ac: среднее по всем автокорреляциям существующих функций (или 0 или 1) всех водородных связей.
-dist: распределение по расстоянию (r) Донор-Акцептор для всех водородных связей
-ang: распределение по углу Донор-Водород-Акцептор а для всех водородных связей
-hx: число водородных связей как функция времени, разделенных на группы n-n+i, где n и n+i означают номера молекул и i принадлежит от 0 до 6. Группа с i=6 влключает все водородные связи для i>6. Эти группы включают n-n+3, n-n+4 и n-n+5 водородных связей связанных со спиралями в протеинов.
-life: время жизни водороных связей вычисляется для усредненных по все автокорреляционым ф-циям существующих ф-ций (или 0 или 1) всех водородных связий: С(тау)=<Si(t)*Si(t+тау)>, где Si(t) = {0,1} для i водородной связи в момент времени t. Интеграл тауHB от С(тау) берется с грубой оценкой для среднего времени жизни водородной связи от 0 до бесконечности по d(тау). Оба интеграл и полная автокорреляционная ф-ция С(тау) будет выходными, так что можно будет проводить более сложный анализ (например, с использованием мульти-экспоненты) для получения более точных оценок для времени жищзни тауHB. Более полный анализ в работе. [135], один из более хороших обзоров у Luzar и Чендлер - анализ кинетики водородных связей [136, 137].
-hbm: карта жизни водородных связей будет сгенерирована на измеренях #H-bonds x #frames. Упорядочение такой же как в индексом файле (см. ниже), но вспять, это означает, что последние триплет в индексном файле соответствует первой строке существования на карте.
Порядок идентичный тому, что и в -hbn индексном файле.
-hbn: индексные группы выходя содержат проанализированные группы, все Донор-Водородные атомные пары и Акцепторные пары в этих группах, Донор-Водород-Акцептор триплеты в водородных связях между анализируемой группой и всеми растворенными атомами участвующих в вставке.
-dan: запись выходных номеров доноров и акцепторов проанализированных для каждого временного отрезка. Это особенно полезно, когда используется -shell
-nhbdist: вычисление числа HBonds от водородов в порядке сравнения результатов с спекстроскопией Рамана
Замечание: опции -ac, -life, -hbn и -hbm требуют наличия памяти пропорциональной общему числу донерного времени, общему числу акцепторного времени выбранной(ых) группы(групп).
Файлы:
-f traj.xtc Input Trajectory: xtc trr trj gro g96 pdb cpt
-s topol.tpr Input Run input file: tpr tpb tpa
-n index.ndx Input, Opt. Index file
-num hbnum.xvg Output xvgr/xmgr file // /* выходной 4 */файл с графиком колличества водородных связей за отрезок прогона
-g hbond.log Output, Opt. Log file // /* выходной 3 */ лог файл содержаший в себе три столбца номеров атомов: донер, водород, акцептор.
-ac hbac.xvg Output, Opt. xvgr/xmgr file // /* выходной 8 */ файл с графиком автокореляционной ф-ции водородных связей /* строится ACF число/число */
-dist hbdist.xvg Output, Opt. xvgr/xmgr file // /* выходной 5 */ файл с графиком распределения водородных связей (Водород-Акцептор (расстояние))
-ang hbang.xvg Output, Opt. xvgr/xmgr file // /* выходной 6 */ файл с графиком распределения водородных связей (Донор-Водород-Акцептор (угол))
-hx hbhelix.xvg Output, Opt. xvgr/xmgr file // /* выходной 7 */ 0 файл с графиком водородных связей за отрезок прогона
-hbn hbond.ndx Output, Opt. Index file // /* выходной 2 */ индексный файл содержаший индексные группы: выбранная группа при обработке файла траекторий, водородные доноры этой группы, акцепторы этой группы, список водородных связей
-hbm hbmap.xpm Output, Opt. X PixMap compatible matrix file // /* выходной 10 */ выходной файл с картой жизни водородных связей
-don donor.xvg Output, Opt. xvgr/xmgr file // /* выходной 1 */ файл с графиком числа доноров за отрезок прогона (свойства доноров)
-dan danum.xvg Output, Opt. xvgr/xmgr file // /* выходной 11 */ файл с графиком соотношения доноров и акцепторов
-life hblife.xvg Output, Opt. xvgr/xmgr file // /* выходной 9 */ файл с графиком времени жизни неразрущенных водородных связей /* строится процентами % */
-nhbdist nhbdist.xvg Output, Opt. xvgr/xmgr file
Дополнительные опции:
-[no]h bool no
-nice int 19
-b time 0 //первый фрагмент (отрезок) пс для чтения директории
-e time 0 //последний фрагмент (отрезок) пс для чтения директории
-dt time 0
-[no]ins bool no // анализ растворителя вставки
-a real 30 // размер угла (градусы, Акцептор - Донор - Водород), угол OOH
-r real 0.35 // радиус обрезания (нм, X - Акцептор, см.след.опцию), например растояние Roh
-[no]da bool yes // использовать растояние Донор-Акцептор (если истина) или Водород-Акцептор (если ложь)
-r2 real 0 // второй радиус обрезания. Главным образом полезет с -contact и -ac
-abin real 1 // распределение по углам Бинвинда (градусы)
-rbin real 0.005 // распределение по углам Бинвинда (нм)
-[no]nitacc bool yes // расценивать азотные атомы как акцепторы
-[no]contact bool no // не смотреть водородные связи, но просто связывать в пределах радиуса обрезания
-shell real -1 // когда > 0, вычисляются воодородные связи в пределах # нм shell около (вокруг) одной частицы
-fitstart real 1 // Время (ps) с когорого начинает работать подгоночная ф-ция для ACF в соответствии с передовой и обратной (остальной) оценочной постоянной для разрыва и формирования HB.
-temp real 298.15 // температура (K) для вычисления энергии Гиббса соответствующая разрыву и преобразованию HB
-smooth real -1 // Если >= 0, конец ACF будет сглажен выбранной подгоночной экспоненциальной функцией: y = A exp(-x/tau)
-dump int 0 // выбросить первые N водородных связей ACFs в одиночном /*единственном*/ xvg file для отладки
-max_hb real 0 // Теоретический максимум числа водородных связей использумемых для нормализации HB автокорреляционной ф-ции. Может быть использован в случае, если программа оценивает не правильно
-[no]merge bool yes // H-bonds между теми же донором и акцептором, но с другим водородом рассматриваемым как одиночая водородная связь. Главным образом важно для ACF.
-acflen int -1 // длинна для ACF, по умолчанию половина числа фрагментов.
-[no]normalize bool yes // норамализовать ACF
-P enum 0 // порядок полинома Лежандра для автокорреляционной ф-ции (0 указывает на none): 0,1,2 или 3
-fitfn enum none // подгоночная функция: none, exp, aexp, exp_exp, vac, exp5, exp7 or exp9
-ncskip int 0 // пропустить N точек в выходном файле автокорреляционной ф-ции
-beginfit real 0 // время, для начала подгоночной экспоненциальной ф-ции для CF
-endfit real -1 // время, для конца подгоночной экспоненциальной ф-ции для CF, -1 пока не достигнет конца.
========
Также рекомендую внимательно смотреть на консольный вывод программы, из него можно тоже под черпнуть важные вещи, как число водородных связей, общее число пар атомов, данные по ACF, а также ссылки на статьи по различным аспектам связанным с этой программой.
P.S. Если статья была вам полезна проголосуйте в опросе блога, этим вы скажите СПАСИБО Автору, также можно обсудить её на Форуме.
Назад к оглавлению : руководство по Gromacs
0 коммент.:
Отправить комментарий
Ваш комментарий будет "принят к сведению" и пойдет на рассмотрение.