 |
Открытый исследовательский проект, к которому может присоединится каждый ... |
-
sad-RNAInSpace (Semi-automatic designing RNA in space) - программное обеспечение для полуавтоматического конструирования[1] РНК в пространстве. Обеспечивает 3D визуализацию структуры РНК, позволяет её изменять и с помощью связи с модулем RNAWorld позволяет автоматизировать некоторые этапы сворачивания РНК.
Данный проект развивается участниками Викиверситета для проведения исследований в трех направлениях: геномика, искусственный интеллект, гибкое программное обеспечение высокого качества.
Все авторские права самого программного обеспечения (sad-RNAInSpace и RNAWorld) принадлежат основателю проекта Сергею Яковлеву, но участникам проекта программное обеспечение может быть предоставлено для реализации научных и образовательных целей.
[править] Результаты моделирования
[править] Важные инструменты
[править] Наиболее важные выводы
Эти выводы получены как обобщение ряда экспериментов, но не обладают строгим доказательством. Тем не менее они призваны направлять последующие исследования.
- Для односпиральных цепей РНК форму при сворачивании предопределяют водородные связи. То есть достаточно добиться образования водородных связей и односпиральная цепь РНК будет соответствовать тому, что наблюдается в природе лишь с небольшими погрешностями в расположении боковых цепей нуклеотидов, которые не обладают водородными связями (как правило петля фрагмента).
- Для многоспиральных цепей РНК форма при сворачивании не предопределяется водородными связями. То есть все водородные связи принципиально можно образовать при таких формах цепи РНК, которых не существует в природе, или которые существуют лишь в процессе сворачивания. Тем не менее таких вариантов существует не более десятка, в основном как различные виды симметрии или недосвернутого состояния. И кроме того, окончательно сформировать все водородные связи у неправильных структур так пока и не удалось, хотя 80% были сформированы, а оставшиеся были в очень близком удалении. Это может указывать на то, что неправильные структуры хотя и находятся близко к правильным (в энергетическом смысле, но не геометрическом), тем не менее водорные связи у них окончательно образовать сильно затруднительно.
- При моделировании сворачивания у РНК часто раньше соединяются концы фрагмента, а затем по направлению к петли (середине фрагмента). Это создает принципиальную невозможность для прямых расчетов, так как невозможно зафиксировать концы и рассчитывать изменения внутри петли, так чтобы не происходило движения на концах или не допустить разрыва цепи.
- Поэтому важно удержать концы развернутыми до того как образуются водородные связи начиная с внутренней части фрагмента
- Возможно, концы цепей надо нарастить временными "опалубочными", эка-теломерными сегментами нуклеотидов. С помощью которых удастся управлять правильным формированием 3D конфигурации всей цепи. Возможно, в процессе построения 3D структуры, при репликации--размножении цепи "теломерные" звенья на концах цепи отщепляясь, позволяют контролировать 3D структуры? — Эта реплика добавлена участником 95.25.203.97 (о · в)
- Проблема сейчас заключается в том, что нельзя достаточно точно (грубо можно) соединить две спирали между которыми есть неканонические водородные связи, которые собственно и стабилизируют и определяют жесткое 3D положение многофрагментной РНК (например, рибозима). Проблема сугубо вычислительная (математическая, геометрическая), но достаточно сложная (нужен профессиональный математик!). Концы многоспиральной цепи в этом не принимают заметного участия, а для односпиральной цепи, конечно, концы цепи важны - но тут проблем не возникает. Поэтому "опалубочные" концы ничего нового не дадут при моделировании (возможно, в отличии от процесса "в пробирке"). --S.J. 06:36, 2 августа 2011 (UTC)
[править] Базы данных углов РНК
- рРНК Haloarcula marismortui (1FFK)
- рРНК Haloarcula marismortui (1JJ2)
[править] Программная разработка
-
-
Предметные причины создания описываются в статье «Геномика бросает вызов искусственному интеллекту».
SAD-RNAInSpace 1.0 - сделан реинжениринг программного обеспечения VMD, создана графическая библиотека визуализации молекул OpenGLPlus, и совместно с модулями RNAWorld включено в данное программное обеспечение. В статье Модуль RNAWorld описывается основной модуль, осуществляющий расчеты по сворачиванию РНК.
[править] Разработка теоретических основ
Не следует понимать данную часть как „разработку теории на бумаге“. Тут скорее наоборот, каждое теоретическое положение проходит экспериментальную проверку in silico со сверкой полученных результатов с биологически (in vitro) известными данными. Если такое положение себя зарекомендовало далее ложится в основу дальнейших гипотез. Иначе признается не удовлетворительным или не эффективным, выводы о чем не менее важны, так как позволяют эффективно направлять ход исследований.
[править] Подход № 1. Обучение с подкреплением
Это подход заморожен — признан недостаточным, детали см. в архиве Модуль RNAWorld/Архив#Версия RNAFoldingAI 0.1.
[править] Подход № 2. Поиск на дереве — понятие о параллельно идущих процессах
Это похоже на игру в шахматы, при том что правила игры одному из играющих неизвестны.
Эта задача может быть определена как игра с природой — игра, в которой имеется только один игрок, причем исход ее зависит не только от его решений, но и от состояния “природы”, то есть не от сознательно противодействующего противника, но от объективной, невраждебной действительности.
Разработка и проверка методов производится на основании данных о рибозиме NC_003540.
[править] Разработанные методы
[править] Подход № 3. Целенаправленное построение многоспиральных цепей РНК
Указанные здесь поднаправления заморожены, так как для многоспиральных цепей не ведут к решению.
[править] Подход № 4. Разрезание многоспиральной цепи РНК, и приведение проблематики к задаче докинга
[править] Биоинформационные ньюансы
[править] Отсутствие симметрии в моделировании
Следует иметь введу, что в третичной структуре отсутствует симметричность, в отличии от вторичной. То есть, последовательность ga≠ag или gagg≠ggag. Поэтому атомные цепи будут разные — координаты атомов будут разные, несмотря на то, что повороты одни и же, для соответствующих нуклеотидов. А так же их энергетические оценки будут существенно разные, так как они зависят от координат атомов.
[править] Последовательность сворачивания
 |
Во-первых, сворачивание начинается со стартовой конформации, поддерживаемой рибосомой. Во-вторых, поиск путей сворачивания идет последовательно с растущего N-концевого участка. В-третьих, в процессе сворачивания С-конец фиксирован.[2] |
 |
[править] Гибкость белковой цепи и РНК-цепи
Колебания валентных связей вообще не вносят вклада в гибкость цепи. Колебания валентных углов вносят малый вклад с амплитудой ≈5°. Гибкость цепи в основном обеспечивается вращением (не полностью свободным) вокруг валентных связей. Это вращение описывается как взаимное расположение атомных группировок, связанных валентной связью через двугранный угол.
Свобода вращение зависит от гибридизации атомов:
- для двух sp3- гибридизированных атомов (например, связи H3C-CH3, CH2C-CH2C) размер флуктуаций ≈15-20°. Такие связи характерны для боковых групп.
- для двух sp2- гибридизированных атомов (например, связи C-N) размер флуктуаций ≈10°
- для связи, соединяющей sp3- гибридизированный и sp2- гибридизированный атомы (например, N-C`, C`-C в главной цепи) существует почти свободное вращение
[править] Математический подход
У этого разрабатываемого подхода есть некоторая аналогия со следующим подходом:
Схематично такая процедура может быть описана следующим образом. Ландшафт энергии системы заливается водой, вода образует лужи вокруг локальных минимумов. При повышении уровня воды лужи начинают сливаться, пока не останется всего одна лужа. Результатом процедуры является направленное дерево бассейнов (луж), причем минимальным лужам ставится в соответствие их глубина (значения локальных минимумов энергии), более крупным лужам — уровень энергии, на которых эти лужи сливаются (величина барьера активации между бассейнами). В результате мы получаем дерево бассейнов и функцию на дереве бассейнов — распределение энергий локальных минимумов и барьеров активации. После этого строится модель иерархической динамики, основанная на дереве бассейнов и функции барьеров активации на нём. Это и есть модель межбассейновой кинетики. [3]
Но наш кибернетически-геометрический подход выгодно отличает, отсутствие сложной математики, благодаря чему его можно реализовать на практике.
Близкие математические работы:
- С. В. Козырев, Методы и приложения ультраметрического и p-адического анализа: от теории всплесков до биофизики
- p-Адическая математическая физика: основные конструкции и применения к сложным и наноскопическим системам, 2008
- Молекулярная динамика биополемеров
- Predicting protein structures with a multiplayer online game
[править] Научные дискуссии
- Дискуссия: Биологические аспекты сворачивания РНК
- Дискуссия: ИИ - аспекты при исследовании сворачивания РНК
Частично содержат такие же выводы как были получены в данном проекте.
-
- успешное предсказание РНК третичной структуры остается исключительными подвигами :)
|