WWW.DISS.SELUK.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА
(Авторефераты, диссертации, методички, учебные программы, монографии)

 

На правах рукописи

Жмуров Артём Андреевич

Моделирование больших биомолекул и

биомолекулярных систем

с использованием графического

процессора

Специальность 05.13.18 – математическое

моделирование, численные методы и комплексы

программ

Автореферат диссертации на соискание учёной степени

кандидата физико-математических наук

Москва - 2011

Работа выполнена на кафедре вычислительной математики Московского физико-технического института (государственного университета)

Научный руководитель:

кандидат физико-математических наук, доцент Холодов Ярослав Александрович

Официальные оппоненты:

доктор физико-математически наук, профессор Ефремов Роман Гербертович доктор физико-математически наук, профессор, кандидат биологических наук Пантелеев Михаил Александрович

Ведущая организация:

Московкий государственный университет им. М.В. Ломоносова, Кафедра биофизики

Защита диссертации состоится “27” октября 2011 г. в 1200 час. на заседании диссертационного совета Д 212.156.05 при Московском физико-техническом институте (государственном университете) по адресу: 141700, г. Долгопрудный Московской обл., Институтский пер. д.9, ауд. 903 КПМ.

С диссертацией можно ознакомится в библиотеке МФТИ (ГУ).

Автореферат разослан “26” сентября 2011 г.

Учёный секретарь диссертационного совета Д 212.156.05 Федько О.С.

Общая характеристика работы

Актуальность темы. Большие белковые молекулы, образования и волокна играют важную роль в жизнедеятельности организма. Нити фибрина, образуя сложную ветвистую структуру, формируют сгусток и останавливают кровотечение. Различные патологии, влияющие на механические свойства фибрина, способны нарушить его функцию и привести как к образованию тромба и возможному инсульту или инфаркту, так и к повышеному кровотечению. Капсиды животных и растительных вирусов способны выдерживать колоссальное давление со стороны упакованного генетического материала. Их физические свойства, а также переходы между стабильным и нестабильным состояниями определяют жизненный цикл вирусов, в том числе созревание вируса и заражение клеток.

Современные эксперименты по изучению одиночных молекул такие как атомно-силовая микроскопия и оптические пинцеты широко используются для изучения механических свойств белковых волокон и капсул вирусов. Однако, ввиду сложности строения этих систем ( 103 105 частиц) и их больших размеров ( 50 200 нм), результаты подобных экспериментов сложно интерпретировать без понимания микромеханики исследуемых биомолекулярных систем.

Стандартные вычислительные методы молекулярной динамики (МД) в полноатомном разрешении широко используются для изучения поведения биомолекул. Но полноатомное моделирование в настоящее время ограничено размером молекулы в 10-50нм и длительностью процесса моделирования в 0, 1 10мс. Поэтому данный метод хорошо подходит только для моделирования равновесных процессов, а достижение биологически важного временного интервала от микросекунды до секунды практически невозможно даже для малых систем. В основе упрощённого метода молекулярного моделирования лежит предположение, что не все степени свободы биомолекулы одинаково важны, и за структурные единицы системы берутся не атомы системы, а их группа (например, одна аминокислота). Такой метод не только позволяет существенно сократить количество производимых вычислений, но и перейти к большему шагу по времени без потери численной стабильности системы. Графические процессоры (ГП), которые изначально были спроектированы для ускорения работы с трёхмерной графикой, способны выполнять многие вычислительные задачи, в том числе и те, которые не связаны с обработкой изображения. Программные платформы для современных ГП включают NVIDIA Compute Unied Device Architecture (CUDA) и Open Computing Language (OpenCL). CUDA, программная среда для параллельных вычислений является высокоуровневой программной платформой, расширяющей стандартные языки С и C++. Это позволяет разработчику реализовывать процедуры (ядра), которые могут выполняться одновременно во множестве независимых потоков на ГП. Комбинация упрощённого метода моделирования и высокой производительности современных ГП обладает огромным потенциалом в области молекулярного моделирования.

Ещё одним преимуществом ГП является их доступность: цена даже самых современных устройств не превышает 500 2000 долларов США, а обычный персональный компьютер, оборудованный двумячетырьмя такими устройствами, способен достигать производительности компьютерного кластера с 200 500 вычислительными ядрами.

Цели работы, объекты исследований. Объекты исследований диссертации - мономер и димер фибриногена, его фрагменты, а также капсид вируса HK97. Известные экспериментальные данные, полученные при помощи атомной силовой микроскопии, показали, что молекулы фибриногена обладают сложными механическими свойствами. Однако, эти данные не позволили точно определить источник динамических характеристик молекулы. Именно поэтому необходимо провести компьютерное молекулярное моделирование силовой денатурации фибриногена. К сожалению, ни один из существующих вычислительных методов не позволяет проводить молекулярное моделирование в условиях воздействия на молекулу, идентичных экспериментальным. Поэтому было решено адаптировать вычислительные методы для эффективной работы на ГП. Полученная реализация также была применена для объяснения динамических характеристик капсида вируса HK97.

Основной целью данной работы являлось моделирование экспериментов атомной силовой микроскопии (АСМ) на единичной молекуле. Для реализации была выбрана модель Самоорганизующегося Полимера (англ. Self Orginized Polymer, SOP). Использование модели SOP, полностью реализованной на ГП, дало возможность наблюдать микромолекулярную динамику белковых систем размером до 106 аминокислот.

Целью данной работы также являлась разработка программной реализации численных методов упрощённого моделирования белковых молекул с использованием графических процессоров. Такая комбинация позволяет моделировать большие белковые системы на экспериментальных временных интервалах без использования дорогих и сложных в обслуживании вычислительных кластеров.

Научная новизна. Данная тема является актуальной с точки зрения развития методологии комплексного исследования больших биомолекул на основе математического моделирования и вычислительных экспериментов. Представленная программная реализация для моделирования больших белковых систем на экспериментальных временных интервалах - единственная из известных на сегодняшний день. Результаты численных исследований, полученные для мономера, димера фибриногена и его фрагментов, позволили поновому интерпретировать экспериментальные данные и лучше понять микромеханику данного белка. Моделирование капсида HK показало, что динамические характеристики данной белковой системы зависят от скорости и геометрии воздействия.

Практическая ценность. Математическая модель и её численная и программная реализации были организованы как программный пакет SOP-GPU, который распространяется в виде открытого кода, сопровождается детальным руководством пользователя и не требует навыков программирования для своего использования. Таким образом, он может быть полезен другим научным группам для исследования интересующих их биомолекулярных систем.

Основные положения, выносимые на защиту:

1. Модификация математической модели самоорганизующегося полимера (SOP) и разработка численных методов для работы на графическом процессоре (ГП).

2. Реализация генераторов псевдослучайных чисел для использования в динамике Ланжевена на ГП.

3. Программный пакет SOP-GPU, размещённый в свободном доступе в сети интернет.

4. Сравнительный анализ экспериментальных данных АСМ и данных молекулярного моделирования силовой денатурации, полученных для мономера, димера фибриногена и различных фрагментов этого белка.

5. Различные динамические режимы, полученные при моделировании силовой индентации капсида вируса HK97.

Апробация работы. Полученные результаты были представлены на 238-ой и 240-ой национальных встречах американского химического общества (Вашингтон, 2009 и Бостон, 2010), на 12-ом и 13-ом конгрессах международного общества тромбоза и гемостаза (Бостон, 2009 и Киото, 2011), на научно-практической конференции “Вычисления с использованием графических процессоров в молекулярной биологии и биоинформатике” (Москва, 2010), на 53-ей научной конференции МФТИ (Долгопрудный, 2010), на конференция “Опыт и результаты исследований, проводимых под руководством приглашенных ученых-соотечественников” (Москва, 2011), на 18ой конференции “Математика. Компьютер. Образование.” (Дубна, 2011) и на 17-ой конференции по структуре и динамике биомолекул (Олбани, 2011).

Публикации. Материалы исследований опубликованы в 16-ти работах [1–16], в том числе семи [1–7] - в изданиях, рекомендованых ВАК РФ. Работы [3, 4] - в изданиях по специальности 05.13.18.

Структура и объём работы. Работа состоит из введения, шести глав, заключения, списка использованных источников, включающего 225 наименований, и трёх приложений. Диссертация изложена на 124-ёх страницах.

Основное содержание работы

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

В первой главе даётся краткий обзор методов молекулярного моделирования. Дан обзор основ молекулярного моделирования в явном и неявном растворителе, объяснены предпосылки использования упрощённых моделей. Приведена потенциальная функция модели самоорганизующегося полимера (Self Orginized Polymer, SOP), используемой в данной работе. Данная модель была разработана для описания механических свойств белков и моделирует процесс экспериментов на одиночных молекулах белка (таких как атомная силовая микроскопия и экспериментов при помощи оптических пинцетов). За центры взаимодействия в модели SOP берутся атомы C, а потенциальная энергия состояния белка зависит от координат {r} = r1, r2,..., rN этих атомов и задаётся выражением:

В уравнении 1, первое слагаемое описывает ковалентные связи в цепи при помощи конечно-растяжимого нелинейного эластичного потенциала (англ. Finitely Extensible Nonlinear Elastic, FENE) VF EN E.

Расстояние между соседними частицами i и j задаётся как rij, а rij - его равновесное значение (значение в нативном состоянии белка), R0 = 2 - индекс чувствительности к изменениям ковалентA ной связи. Второе слагаемое, описываемое потенциалом ЛеннардаДжонса (VN B T ), используется для описания нековалентных связей, стабилизирующих нативное состояние (нативных контактов). Значение параметра h = 0, 7 1, 6 ккал/моль характеризует силу нековалентных связей. Все прочие пары атомов (не связанные ковалентно и не формирующие нативный контакт) взаимодействуют согласно потенциалу VN B, благодаря которому становятся невозможными самопересечения цепи. Параметры r = 1.0 ккал/моль и = 3.8 постоянны и характеризуют силу и радиус отталкивания.

Основным преимуществом модели SOP является то, что несмотря на простоту силового поля, она достаточно точно описывает механическую денатурацию белков. Кроме того, благодаря использованию динамики Ланжевена в задемпфированом пределе, вычислительные процедуры остаются стабильными даже при использовании большого шага по времени в 40пс. Данная модель была успешно использована для достаточно больших систем без применения аппаратного ускорения. Однако моделирование больших систем на экспериментальных временных интервалах с использованием центрального процессора невозможно даже с применением модели SOP.

Поэтому, в данной работе модель была реализована на графических процессорах, что позволило обойти данный барьер.

Во второй главе приведен обзор публикаций по развитию графических процессоров (ГП) как альтернативных вычислительных устройств высокой производительности. Приведены примеры использования как современных ГП, так и устройств предыдущих поколений. Также во второй главе даётся обзор языка программирования CUDA, основ устройства аппаратной части и примеров его использования для научных вычислений. В последнем разделе второй главы приведены методы для оптимизации работы программы, реализованной на ГП. В биомолекулярном моделировании с использованием молекулярной динамики (МД) или динамики Ланжевена (ДЛ) взаимодействия между частицами описываются при помощи потенциальной функции (силового поля). Функциональная часть силового поля одинакова для всех частиц в системе, а уравнения движения всех частиц могут быть численно проинтегрированы независимо друг от друга. Благодаря этому, численные алгоритмы молекулярной динамики и динамики Ланжевена могут быть успешно реализованы на ГП.

В третьей главе предлагаются численные методы и алгоритмы, использованные при реализации модели SOP на ГП. Показаны два раздельных подхода к параллелизации на ГП: “Распараллеливание по частицам” и “Распараллеливание по парам взаимодействующих частиц”. Приведены процедуры создания списка Верле, численного решения уравнений движения. Далее полученная реализация была проверена на предмет аппаратных и программных ошибок. Для этого данные, полученные на графическом процессоре (ГП) были сравнены с данными, полученными на центральном процессоре (ЦП).

Результаты вычислений сравнивались также с данными, полученными аналитически. Было показано, что даже при использовании одинарной точности результаты численного эксперимента хорошо соответствуют аналитическому решению. Была проверена целесообразность использования схемы интегрирования более высокого порядка и показано, что в данном случае первого порядка точности достаточно даже при использовании шага по времени в 20-40пс. В работе впервые реализован подход “много запусков на ГП”. Данный подход позволяет полностью использовать вычислительные возможности ГП даже при моделировании сравнительно небольших систем.

В дальнейшем этот подход может быть легко расширен для использования в различных методах рассчёта термодинамики системы.

Рис. 1: Относительная производительность, или отношение скорости вычислений на ГП к скорости вычислений на ЦП (ускорение), полученное для приложения SOP-GPU как функция размера системы Ntot или количества траекторий s. Также показаны графики ускорения при использовании кеширования текстур и аппаратного ускорения математических функций. Для наглядности, график показан с логарифмическим масштабом по оси x.

Итоговая производительность разработанного приложения для ГП (программа SOP-GPU) была сравнена с оптимизированной версией реализации того же приложения для ЦП. При этом оба приложения были использованы для описания динамики Ланжевена одной и той же молекулярной системы - домена W W в состоянии равновесия (Рис. 1). Для полного задействования ресурсов ГП представленна вычислительная производительность программы SOP-GPU как функция от числа независимых траекторий, обсчитываемых одновременно на одном ГП (подход “много-траекторий-на-одном-ГП”).

Альтернативный этому вариант заключался в расчете одной траектории на одном ГП, но для систем разного размера (подход “однатраектория-на-одном-ГП”). Результаты показывают, что для небольшой системы из 34 аминокислотных остатков (домен W W ), использование одного ГП позволяет ускорить вычисления при получении трёх и более независимых траекторий, так как небольшие системы не могут полностью загрузить все параллельные вычислительные ресурсы ГП. Эквивалентную загрузку одного ГП обеспечивает запуск обсчёта одной траектории для системы из 102 частиц. В то время, как время моделирования на ЦП линейно возрастает с увеличением количества запускаемых потоков, нагрузка на ГП в таком режиме остаётся сублинейной (примерно постоянной) вплоть до 500 траекторий. Далее, ГП демонстрирует значительное превосходство в производительности над ЦП, достигая максимального ускорения в 8090 раз (Рис. 1). Для сравнения производительности были использованы ГП и ЦП одинакового технологического уровня (ГП NVIDIA GeForce GTX 295 и ЦП Quad Core Xeon 2, 66 ГГц).

В четвёртой главе приведено описание реализации генераторов псевдослучайных чисел (ГСЧ) на ГП. В силу того, что ГП не так давно стали широко использоваться в научных вычислениях, многие базовые вычислительные методы пока не реализованы и не оптимизированы для работы на этих устройствах. В этой главе рассматривается четыре различных генератора случайных чисел: линейный конгруэнтный генератор LCG, гибридный генератор “гибридный Таус”, генератор Ran2 и Фибоначчи с запаздыванием. В первой части главы даётся описание математических основ трёх методов. Далее показано, как эти методы могут быть эффективно адаптированы для работы на ГП.

Для разработки параллельных реализаций ГСЧ на ГП был использован метод разделения цикла генератора. Основная идея метода заключается в распределении одной последовательности ГСЧ, которая может быть рассмотрена как периодический цикл случайных чисел, между множеством потоков, каждый из которых производит подпоследовательность случайных чисел. Большинство ГСЧ, включая LCG, Ran2 и гибридный Таус, основаны на последовательном преобразовании текущего состояния, поэтому наиболее естественным способом получения разных последовательностей случайных чисел в разных потоках является сообщение потокам различных начальных значений состояния. При этом начальные состояния распределяются между потоками таким образом, чтобы исключить межпоточные корреляции случайных чисел. На этом основан подход “один-ГСЧ-на-поток”. Существуют ГСЧ, такие как вихрь Менсенна и алгоритм Фибоначчи, в которых можно перескочить вперёд по последовательности и вычислить случайное число (n+1), не вычисляя перед этим n-ное. Длина запаздывания, которая зависит от параметров ГСЧ, может быть выбрана больше числа потоков (которое в случае молекулярной динамики или динамики Ланжевена равно числу частиц). Тогда все N чисел могут быть получены одновременно, то есть. j-тый поток вычисляет j-тоe, (j + N )-тоe, (j + 2N )-тоe и так далее числа. В этом случае, только одно состояние ГСЧ используется для всех потоков, каждый из которых обновляет только один элемент этого состояния. Такой подход можно назвать “один-ГСЧна-все-потоки”. Реализации ГСЧ затем были проверены на предмет возникновения корреляций, также была проверена их производительность и требования к памяти.

В пятой главе описаны результаты численного моделирования силовой денатурации мономера, димера фибриногена и его фрагментов. Полученные результаты позволили по-новому интерпретировать доступные экспериментальные данные. Для исследования механизмов денатурации и для получения энергетических характеристик были использованы следующие системы: изолированный модуль, -модуль с 25-ю примыкающими аминокислотами -модуля, комплекс и -модулей (глобулярная часть домена D), фрагмент “двойной D-домен”, мономер и димер фибриногена. Для каждой системы были получены от 6 до 10 траекторий.

Экспериментальные и теоретические графики зависимости силы от растяжения для мономера, димера и олигомеров фибриногена обладали очень похожими характеристиками (Рис. 2 и 3). Статистический анализ, усреднённый по пяти траекториям моделирования, длиной 0, 2 секунды каждая, показал, что среднее значение расстояния между максимумами x 29нм для мономера фибриногена и x 27нм для димера, что хорошо соответствует значениям, полученным экспериментально. Среднее значение сил денатурации оказались равными, соответственно, f 134 и f 112пН для мономера и димера фибриногена. Эти значения несколько превосходят экспериментальные данные, что может быть следствием слишком Рис. 2: Экспериментальные результаты денатурация фибриногена при помощи АСМ. Графики зависимости силы от растяжения были получены при помощи скорости движения зонда в f = 1, 0мкм/сек для мономера фибриногена (панель (а)) и его олигомеров (панель (б)). Гистограммы распределения расстояний между пиками и сил денатурации показаны на панелях (в) и (г) соответственно. На панели (в) также показана суперпозиция трёх нормальных распределений (сплошная линия) расстояний между максимумами сил, полученных из результатов моделирования для переходов типа 1-3 (пунктирные линии), которые изображены на Рис. 3).

сильной стабилизации структуры белка в модели SOP. Главной целью моделирования было получение значений для расстояний между максимумами, значения которых хорошо соответствуют экспериментальным данным.

Основным результатом является то, что механизм денатурации фибрин(оген)а зависит от нескольких структурных элементов этой Рис. 3: Кривые зависимости силы от растяжения для мономеров (панель (а)) и димеров (панель (б)) фибриногена, полученные при помощи молекулярного моделирования с применением скорости зонда в f = 1, 0мкм/с. Точки на панелях (в) и (г) показывают расстояния между максимумами и силы денатурации, соответственно, сгруппированные для переходов типа 1-3. Максимы силы, соответствующие этим переходам, пронумерованы на графиках (панели (а) и (б).

комплексной молекулы. Главным образом, процесс ответа фибриногена на механическое воздействие состоит из обратимого процесса растяжения/сжатия свёрнутых -спиралей и последовательной денатурации глобулярных элементов -модуля. В каждом -модуле, первым событием денатурации является отделение C-терминального -луча (аминокислоты 380 392). Сразу вслед за этим, денатурирует центральная часть -модуля (аминокислоты 234 311), который при этом разделяется на два глобулярных фрагмента - Cтерминальный (аминокислоты 311 380) и N -терминальный (аминокислоты 139 234). В результате этого денатурация одного модуля происходит в три этапа, то есть в три перехода денатурации (типа 1, 2 и 3), что соответствует трём максимумам на графике зависимости силы от молекулярного растяжения.

В шестой главе приведены результаты моделирования силовой индентации капсида вируса HK97. Механический ответ капсида показывает стохастическое поведение при малой скорости движения зонда f = 2, 5мкм/с (Рис. 4а), но когда скорость увеличивается, зависимость становится более определённой (Рис. 4б). Значение силы при этом постоянно флуктуирует, что говорит о том, что упругость капсида - стохастическая характеристика. Это также подтверждается графиком зависимости коэффициента упругости K от смещения. Интересным результатом также является то, что с увеличением скорости увеличивается и K. Так, при f = 2, 5мкм/с K = 0, 01 0, 025, при f = 25мкм/с K = 0, 05 0, 075, а при f = 250мкм/с K = 0, 2 0, 35. Это указывает на то, что стенки капсида становятся более упругими при более сильном воздействии. Также, K зависит и от размеров зонда, то есть от геометри воздействия. “Быстрое проминание” (англ. buckling) возникает только при самой маленькой скорости воздействия в f = 2, 5мкм/с и при максимальном размере зонда, использованого в данной работе, R = 25нм. При быстром проминании капсида, K сначала увеличивается, а затем - быстро падает до нуля. Проминание при этом происходит при смещении Z 25нм, при котором K падает до нуля, что говорит о потере механического сопротивления капсидом Рис. 4: Силовая индентация капсида вируса HK97. Показана зависимость силы реакции капсида от смещения зонда. Для того, чтобы промоделировать условия эксперимента, на капсиду давил сферический зонд (серый шар) радиуса R = 5, 10 и 25нм. Зонд двигался вниз (по направлению серой стрелки) со скоростью в f = 2, 5мкм/с (а) и 25мкм/с (б). Также показаны промежуточные структуры, возникающие в процессе индентации. На вложенных графиках показана зависимость количества нативных контактов Q и коэффициент упругости K от смещения зонда Z.

(Рис. 4a). При f = 2, 5мкм/с и R = 10нм индентация капсида монотонна (быстрого проминания не происходит). А при f = 2, 5мкм/с и R = 5нм капсид разрывается, когда Z 55нм, а внешняя сила F 1нН (Рис. 4а). Разрыв происходит в непосредственной близости к зонду, посредством частичной денатурации формирующего капсид белка. Это также заметно на графике зависимости количества нативных контактов Q от смещения Z, на котором видно резкое падение (Рис. 4).

В заключении приведены основные результаты работы.

В приложениях 1 и 2 приведены псевдокоды реализации ГСЧ “гибридный Таус”, Ran2 и Фибоначчи с запаздыванием на ГП.

В приложении 3 показаны графики силовой денатурации, полученные группой Вайзела при проведении экспериментов АСМ по силовой денатурации мономеров и олигомеров фибриногена, которые сравнивались с полученными в диссертации результатами.

Основные результаты диссертации 1. Модификация математической модели самоорганизующегося полимера (SOP) для выполнения на графических процессорах (ГП).

Разработка вычислительных процедур для расчёта межчастичных сил, создания списков Верле, интегрирования уравнений движения. Разработка метода компьютерного моделирования одновременно на одном ГП большого числа траекторий движения молекул, который позволил эффективно использовать вычислительные возможности ГП и получить статистически важную выборку результатов за короткое время.

2. Разработка численных процедур для реализации генераторов псевдослучайных чисел на ГП. Полученные реализации применимы для динамики Ланжевена и других научных приложений с жёсткими требованиями к периоду последовательности случайных чисел, их качеству и к производительности программы.

3. Создание программной реализации модели самоорганизующегося полимера на ГП (пакет SOP-GPU) и её размещение в открытом доступе в сети интернет. Использования SOP-GPU даёт прирост производительности в 90 раз на ГП GeForce GTX 295.

4. Проведение вычислительных экспериментов для комплексного исследования силовой денатурации мономера и димера фибриногена и отдельных фрагментов этой молекулы. Результаты моделирования хорошо согласуются с результатам экспериментов атомной силовой микроскопии (АСМ). Сравнительный анализ экспериментальных данных и данных компьютерного молекулярного моделирования позволил объяснить микромолекулярный механизм денатурации фибриногена.

5. Созданная система компьютерного и иммитационного моделирования была применена для проведения численных экспериментов по силовой индентации капсида вируса HK97. Моделирование показало наличие различных динамических режимов механического ответа оболочки вируса на внешнее воздействие - режима быстрого проминания, упругой индентации и механического разрыва.

Список публикаций по теме диссертации 1. SOP-GPU: Accelerating biomolecular simulations in the centisecond timescale using graphics processors / A. Zhmurov, R. I. Dima, Y. Kholodov, V. Barsegov // Proteins. 2010. Vol. 78, no. 14.

P. 2984–2999.

2. Generation of random numbers on graphics processors: Forced indentation in silico of the bacteriophage HK97 / A. Zhmurov, K. Rybnikov, Y. Kholodov, V. Barsegov // J. Phys. Chem. B.

2011. Vol. 115, no. 18. P. 5278–5288.

3. Моделирование микромеханики биомолекул на графических процессорах с использованием динамики Ланжевена / А. А. Жмуров, В. А. Барсегов, С. В. Трифонов и др. // Мат.

Модел. 2011. Т. 23, № 10. С. 133–156.

4. Эффективные генераторы псевдослучайных чисел для молекулярного моделирования на графических процессорах / А. А. Жмуров, И. Морозов, Я. А. Холодов и др. // Компьют.

Исслед. и Модел. 2011. Т. 3, № 3. С. 296–311.

5. Exploring the mechanical stability of the C2 domains in human synaptotagmin 1 / L. Duan, A. Zhmurov, V. Barsegov, R. Dima // J. Phys. Chem. B. 2011. Vol. 115, no. 33. P. 10133–10146.

6. Bura E., Zhmurov A., Barsegov V. Nonparametric density estimation and optimal bandwidth selection for protein unfolding and unbinding data // J. Chem. Phys. 2009. Vol. 130, no. 1.

P. 015102.

7. Zhmurov A., Dima R. I., Barsegov V. Order statistics theory of unfolding of multimeric proteins // Biophys J. 2010. Vol. 99, no. 6. P. 1959–1968.

8. Molecular simulations of forced unfolding of brin monomers / V. Barsegov, A. Zhmurov, J. W. Weisel et al. // JTH. 2009.

Vol. 7, no. Suppl. 2. P. 1005–1006.

9. Barsegov V., Dima R., Zhmurov A. Molecular simulations of large-size protein assemblies on graphics processors // Abstracts of Papers of the American Chemical Society, http://oasys2.confex.com/acs/238nm/techprogram/P1294978.HTM.

10. Resolving brinogen nanomechanics using dynamic force measurements in vitro and in silico / A. Zhmurov, A. Brown, R. I. Litvinov et al. // J. Biomol. Struct. Dyn. 2011. Vol. 28, no. 6. P. 975–976.

11. Phase transition from -helices to -sheets in brinogen coiled coils / A. Zhmurov, A. Brown, R. I. Litvinov et al. // J. Biomol. Struct.

Dyn. 2011. Vol. 28, no. 6. P. 976–977.

12. Origins of the mechanical stability of the C2 domains in human synaptotagmin 1 / L. Duan, A. Zhmurov, V. Barsegov, R. I. Dima // J. Biomol. Struct. Dyn. 2011. Vol. 28, no. 6. P. 997–998.

13. Fibrin nanomechanics: structural underpinnings and unfolding mechanism underlying forced elongation of brin(ogen) monomers and polymers / A. Zhmurov, A. E. Brown, R. I. Litvinov et al. // JTH. 2011. Vol. 9, no. Suppl. 2. P. 831.

14. Модульная программная платформа для многоуровневого моделирования биомолекулярных систем / А. Жмуров, А. Алексеенко, В. Барсегов и др. // Математика.

http://www.mce.su/rus/archive/mce18/sect101359/doc100272/.

15. Разработка гибридной высокопроизводительной вычислительной платформы для прикладного компьютерного моделирования задач взаимодействия макромолекул / А. Жмуров, И. Коваленко, В. Барсегов и др. // Вычисления с использованием графических процессоров в молекулярной биологии и биоинформатике. Москва: 2010. С. 47–48.

16. Программно-аппаратная платформа для высокопроизводительных вычислений разномасштабных задач молекулярной динамики / С. Трифонов, А. Жмуров, В. Барсегов и др. // Труды 53-й научной конференции МФТИ. Долгопрудный: 2010. С. 31– Личный вклад соискателя в работах с соавторами заключается в разработке модификации математической молекулярной модели, разработке и реализации численных алгоритмов, проведении численных экспериментов и анализе результатов.

Жмуров Артём Андреевич Моделирование больших биомолекул и биомолекулярных систем с использованием графического Подписано в печать: 15.09. Усл.п.л. – 1,5 Заказ № 123 Тираж: 100 экз.

(499) 788-78-56 www.autoreferat.ru

Похожие работы:

«УДК 533.9 МАЛЮТИН Александр Евгеньевич ИОННЫЕ ИСТОЧНИКИ ПОЛЕВОГО ТИПА ИЗ УГЛЕРОДНЫХ МАТЕРИАЛОВ ДЛЯ МАСС-СПЕКТРОМЕТРИИ Специальность 01.04.08 – Физика плазмы АВТОРЕФЕРАТ диссертации на соискание учёной степени кандидата физико-математических наук Долгопрудный – Работа...»

«Вокин Алексей Иннокентьевич ЭКОЛОГИЯ ХАРИУСОВЫХ РЫБ (THYMALLIDAE) ГОРНЫХ ВОДОЕМОВ БАЙКАЛЬСКОЙ РИФТОВОЙ ЗОНЫ 03.00.16 – экология Автореферат диссертации на соискание ученой степени кандидата биологических наук Улан-Удэ – 2008 3 Работа выполнена на кафедре зоологии позвоночных и экологии и кафедре водных ресурсов ЮНЕСКО Иркутского государственного университета Научный руководитель : кандидат биологических наук, доцент Самусёнок Виталий Петрович Официальные оппоненты : доктор...»

«ЧЕСТНОВ ОЛЕГ ПЕТРОВИЧ НАУЧНОЕ ОБОСНОВАНИЕ СТРАТЕГИЧЕСКИХ НАПРАВЛЕНИЙ МЕЖДУНАРОДНОГО СОТРУДНИЧЕСТВА РОССИЙСКОЙ ФЕДЕРАЦИИ С ВСЕМИРНОЙ ОРГАНИЗАЦИЕЙ ЗДРАВООХРАНЕНИЯ 14.00.33 – Общественное здоровье и здравоохранение АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата медицинских наук Москва - 2008 г. 2 Работа выполнена в Федеральном государственном учреждении Центральный научно-исследовательский институт организации и информатизации здравоохранения Федерального...»

«ЦЫРО Светлана Геннадьевна РЕГИОНАЛЬНАЯ МОДЕЛЬ ДЛЯ ИССЛЕДОВАНИЯ ФИЗИЧЕСКИХ И ХИМИЧЕСКИХ СВОЙСТВ ВЗВЕШЕННЫХ ЧАСТИЦ В ЕВРОПЕ Специальность 25.00.30 – метеорология, климатология, агрометеорология Автореферат диссертации на соискание ученой степени кандидата физико-математических наук Санкт-Петербург 2008 Работа выполнена в Главной геофизической обсерватории им. А. И. Воейкова Научный руководитель : кандидат физико-математических...»

«ФЕДАШ Анатолий Владимирович Развитие методологии проектирования и обоснования функциональной структуры горнотехнических систем освоения георесурсного потенциала угольных месторождений Специальность: 25.00.21Теоретические основы проектирования горнотехнических систем Автореферат диссертации на соискание ученой степени доктора технических наук Новочеркасск – 2013 Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального...»

«ПОЖИДАЕВ Александр Петрович АЛГЕБРАИЧЕСКИЕ СИСТЕМЫ ЛИЕВА ТИПА 01.01.06 математическая логика, алгебра и теория чисел АВТОРЕФЕРАТ диссертации на соискание учёной степени доктора физико-математических наук Новосибирск 2010 Работа выполнена в Институте математики им. С.Л.Соболева Сибирского отделения Российской академии наук Научный консультант : доктор физико-математических наук, профессор Шестаков Иван Павлович Официальные оппоненты : доктор физико-математических наук,...»

«Даниленко Ольга Константиновна ПОВЫШЕНИЕ ЭФФЕКТИВНОСТИ ЛЕСОПОЛЬЗОВАНИЯ ПРИ ПОДГОТОВКЕ ЛОЖ ВОДОХРАНИЛИЩ (НА ПРИМЕРЕ БОГУЧАНСКОЙ ГЭС) 05.21.01 – Технология и машины лесозаготовок и лесного хозяйства АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Братск – 2008 2 Работа выполнена в Братском государственном университете. доктор технических наук, профессор Научный руководитель : Угрюмов Борис Иванович доктор технических наук, профессор Официальные...»

«УДК 577.214.4:575.22 ЧЕТВЕРИНА ДАРЬЯ АЛЕКСАНДРОВНА ОБНАРУЖЕНИЕ И ИЗУЧЕНИЕ НОВОГО ИНСУЛЯТОРА У DROSOPHILA MELANOGASTER Специальность 03.00.26 - молекулярная генетика АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата биологических наук Москва 2008 Работа выполнена в лаборатории Регуляции генетических процессов Института биологии гена РАН Научный руководитель : академик РАН, доктор биологических наук, профессор П.Г. Георгиев Официальные оппоненты : доктор...»

«Тараева Галина Рубеновна Семантика музыкального языка: конвенции, традиции, интерпретации Специальность 17.00.02 – Музыкальное искусство Автореферат диссертации на соискание ученой степени доктора искусствоведения Ростов-на-Дону – 2013 Работа выполнена в Ростовской государственной консерватории (академии) им. С. В. Рахманинова Официальные оппоненты : Казанцева Людмила Павловна, доктор искусствоведения, профессор кафедры истории и теории музыки Астраханской государственной...»

«Томин Павел Юрьевич МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ФИЛЬТРАЦИИ В ТРЕЩИНОВАТЫХ СРЕДАХ 05.13.18 — Математическое моделирование, численные методы и комплексы программ Автореферат диссертации на соискание ученой степени кандидата физико-математических наук Москва 2011 2 Работа выполнена на кафедре плазменной энергетики факультета проблем физики и энергетики Московского физико-технического института. Научный руководитель : кандидат физико-математических наук, доцент...»

«ТКАЧУК АРТЕМ ПЕТРОВИЧ РАЗРАБОТКА МЕТОДА СТАБИЛИЗАЦИИ ТРАНСГЕНОВ ПОСЛЕ ИХ ИНТЕГРАЦИИ В ГЕНОМ Специальность 03.01.07 – молекулярная генетика АВТОРЕФЕРАТ Диссертации на соискание ученой степени кандидата биологических наук Москва – 2010   Работа выполнена в группе биологии теломер Учреждения Российской академии наук Института биологии гена РАН Научный руководитель : кандидат биологических наук Савицкий Михаил Юрьевич...»

«Димитриев Юрий Олегович СОВРЕМЕННОЕ ЭКОЛОГИЧЕСКОЕ СОСТОЯНИЕ ФЛОРЫ УРБАНИЗИРОВАННЫХ ТЕРРИТОРИЙ ЛЕСОСТЕПНОЙ ЗОНЫ (НА ПРИМЕРЕ ГОРОДА УЛЬЯНОВСКА) 03.02.01 – Ботаника 03.02.08 – Экология АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата биологических наук Сыктывкар – 2011 Работа выполнена на кафедре ботаники ГОУ ВПО Ульяновский государственный педагогический университет имени И.Н. Ульянова кандидат биологических наук, доцент Научный руководитель : Масленников Андрей...»

«УДК 39 (575.1) (09) 641.55 (575.1) (09) ФАЙЗУЛЛАЕВА МАВЛЮДА ХАМЗАЕВНА НАЦИОНАЛЬНАЯ ПИЩА В ТРАДИЦИЯХ И ОБРЯДАХ НАСЕЛЕНИЯ СУРХАНСКОГО ОАЗИСА (ПЕРВАЯ ПОЛОВИНА XX ВЕКА) 07.00.07 – Этнография, этнология и антропология АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата исторических наук Ташкент – 2010 Работа выполнена на кафедре Всемирная история Термезского государственного университета доктор исторических наук,...»

«Полуэктова Мария Михайловна МЕТОД ОЦЕНКИ ЗАГРЯЗНЕНИЯ АТМОСФЕРНОГО ВОЗДУХА АВТОМОБИЛЬНЫМ ТРАНСПОРТОМ С ИСПОЛЬЗОВАНИЕМ ГЕОИНФОРМАЦИОННЫХ СИСТЕМ Специальность: 25.00.30 - метеорология, климатология, агрометеорология АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Санкт-Петербург 2009 Работа выполнена в государственном учреждении Главная геофизическая обсерватория им. А. И. Воейкова Научный руководитель : Заслуженный деятель науки РФ, доктор...»

«ФИЛИМОНОВА Наталья Владимировна ФРАЗЕОЛОГИЗМЫ, НОМИНИРУЮЩИЕ ЧЕЛОВЕКА ПО ЧЕРТАМ ХАРАКТЕРА, В РУССКОМ И НЕМЕЦКОМ ЯЗЫКАХ: СТРУКТУРНЫЙ И СЕМАНТИЧЕСКИЙ АСПЕКТЫ Специальность 10.02.20 – Сравнительно-историческое, типологическое и сопоставительное языкознание АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата филологических наук Челябинск, 2011 Работа выполнена в Государственном образовательном учреждении высшего профессионального образования Челябинский государственный...»

«Рак Наталья Семеновна ЭКОЛОГО-БИОЛОГИЧЕСКИЕ АСПЕКТЫ ТРОФИЧЕСКИХ СВЯЗЕЙ В СИСТЕМЕ РАСТЕНИЕ – ФИТОФАГ – ЭНТОМОФАГ В ЗАЩИЩЕННОМ ГРУНТЕ ЗАПОЛЯРЬЯ 03.02.08. – Экология Автореферат диссертации на соискание ученой степени доктора биологических наук Петрозаводск – 2012 Работа выполнена в Полярно-альпийском ботаническом саду им. Н. А. Аврорина Кольского научного центра РАН Консультанты: член-корреспондент РАН, профессор В. К. Жиров кандидат биологических наук Л. П. Красавина...»

«Окунев Борис Николаевич Энергосберегающие термодинамические циклы в химико-технологических системах Специальность 02.00.04 — физическая химия АВТОРЕФЕРАТ диссертации на соискание ученой степени доктора физико-математических наук Москва — 2013 год 2 Работа выполнена на кафедре химической технологии и новых материалов Химического факультета Московского государственного университета имени М.В. Ломоносова Официальные оппоненты : Агеев Евгений Петрович, доктор химических наук,...»

«БУКШУК НАТАЛЬЯ АЛЕКСАНДРОВНА ЭКОЛОГИЧЕСКИЕ ОСОБЕННОСТИ ЭНДЕМИЧНЫХ ГУБОК ОЗЕРА БАЙКАЛ: РАСПРЕДЕЛЕНИЕ И ЖИЗНЕННЫЕ ЦИКЛЫ 03.02.08 - экология АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата биологических наук Иркутск – 2014 Работа выполнена в Лаборатории биологии водных беспозвоночных Федерального государственного бюджетного учреждения науки Лимнологический институт Сибирского отделения Российской академии наук (ЛИН СО РАН), г. Иркутск. Научный доктор биологических...»

«Филимонов Иван Сергеевич КИНЕТИЧЕСКИЕ ЗАКОНОМЕРНОСТИ КАТАЛИЗА И ИНАКТИВАЦИИ ФЕРМЕНТА ПРОСТАГЛАНДИН Н СИНТАЗЫ Специальность 03.01.02. – Биофизика АВТОРЕФЕРАТ диссертации на соискание учёной степени кандидата физико-математических наук Москва – 2010 Работа выполнена в Международный учебно-научном биотехнологическом центре МГУ имени М.В.Ломоносова и Институте биохимической физики имени Н.М.Эмануэля РАН. Научный руководитель : доктор химических наук, профессор Вржещ Петр...»

«ДУХАНИН МИХАИЛ ЮРЬЕВИЧ НАУЧНО-ТЕХНИЧЕСКИЙ ПРОГРЕСС КАК ФАКТОР РОСТА ЭФФЕКТИВНОСТИ ПРОИЗВОДСТВА В МОЛОЧНОМ СКОТОВОДСТВЕ Специальность – 08.00.05 – экономика и управление народным хозяйством (экономика, организация и управление предприятиями, отраслями, комплексами – АПК и сельское хозяйство) Автореферат диссертации на соискание ученой степени кандидата экономических наук Москва – 2014 г. 1 Работа выполнена в Федеральном государственном бюджетном образовательном учреждении...»




























 
2014 www.av.disus.ru - «Бесплатная электронная библиотека - Авторефераты, Диссертации, Монографии, Программы»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.