Бесплатный автореферат и диссертация по биологии на тему

Методы моделирования наднуклеосомной структуры хроматина и расчета спектров малоуглового рассеяния нейтронов

ВАК РФ 03.01.02, Биофизика

Автореферат диссертации по теме "Методы моделирования наднуклеосомной структуры хроматина и расчета спектров малоуглового рассеяния нейтронов"

005043470

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

Илатовский Андрей Владимирович

МЕТОДЫ МОДЕЛИРОВАНИЯ НАДНУКЛЕОСОМНОЙ СТРУКТУРЫ ХРОМАТИНА И РАСЧЁТА СПЕКТРОВ МАЛОУГЛОВОГО РАССЕЯНИЯ НЕЙТРОНОВ

03.01.02 — биофизика

АВТОРЕФЕРАТ диссертации на соискание учёной степени кандидата физико-математических наук

1 7 мд? 2012

Санкт-Петербург — 2012

005043470

Работа выполнена в лаборатории биофизики макромолекул Отделения молекулярной и радиационной биофизики ФГБУ «Петербургский институт ядерной физики им. Б.П.Константинова» (г. Гатчина).

Научный руководитель: кандидат физико-математических наук

Исаев-Иванов Владимир Васильевич.

Официальные оппоненты: Сергей Борисович Вахрушев,

доктор физико-математических наук, профессор, ФГБУН Физико-технический институт им. А. Ф. Иоффе РАН, заведующий лабораторией;

Тулуб Александр Александрович, доктор физико-математических наук, Санкт-Петербургский государственный университет, ведущий научный сотрудник.

Ведущая организация: . Объединённый институт

ядерных исследований, , Дубна.

Защита состоится 23 мая 2012 г. в 16 часов на заседании диссертационного совета Д 212.229.25 при ФГБОУ ВПО «Санкт-Петербургский государственный политехнический университет» по адресу:

195251, Санкт-Петербург, ул. Хлопина, д. 5, факультет медицинской физики и биоинженерии, ауд. 305.

С диссертацией можно ознакомиться в фундаментальной библиотеке ФГБОУ ВПО «Санкт-Петербургский государственный политехнический университет».

Автореферат диссертации разослан 2^0. апреля 2012 г.

Учёный секретарь диссертационного совета

доктор физико-математических наук "Л Власова Ольга Леонардовна

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ Актуальность исследования

Ядерная ДНК в геномах эукариот представляет собой нуклеопроте-иновый комплекс — хроматин. Для компактного хранения генетической информации и обеспечения доступа к ней требуется его определённая организация. В середине 1970-х гг. было известно, что хроматин состоит из ДНК и гистонов Н1, Н2А, Н2В, НЗ и Н4. На основании имеющихся данных была выдвинута гипотеза о существовании структурного элемента хроматина, которая позднее была подтверждена экспериментально. Обнаруженные частицы назвали нуклеосомами. В 1997 г. была получена полноатомная структура нуклеосомы.

Достигаемая с помощью иуклеосомной организации степень компак-тизации хроматина является недостаточной для объяснения наблюдаемой плотности упаковки ДНК в ядре, что указывает на наличие наднук-леосомпых структур хроматина. Наблюдаемые при помощи электронной микроскопии нити хроматина диаметром 30 им считаются следующим уровнем компактизации ДНК. За прошедшие десятилетия нити хроматина изучались множеством экспериментальных методов, однако полученные результаты не позволяют сделать однозначных выводов об их структуре, а ряд исследователей ставит под сомнение само существование регулярной упаковки нуклеосом в 30-нанометровые нити.

Определённый прогресс достигнут в исследовании структуры хроматина в масштабе целого ядра. Применение биофизических и биохимических методов позволило установить двухуровневую фрактальную организацию хроматина и, в ряде случаев, построить трёхмерные модели исследуемых геномов. Тем не менее единая молекулярная модель структурной организации хроматина, охватывающая все уровни компактизации ДНК, от нуклеосомного до геномного, пока не создана.

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

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

возможно наблюдать рассеяние нейтронов на различных частях нук-леопротеидпых комплексов, например, при концентрации тяжёлой воды в растворе порядка 65 % спектр МУРН будет обусловлен преимущественно белковой компонентой. В-третьих, спектры МУРН охватывают практически полный интервал размеров наднуклеосомных структур хроматина, от нуклеосомных до геномных. Указанные преимущества использования МУРН для исследования структуры хроматина и других биомакромолекулярных комплексов обуславливают большое число нейтронных экспериментов для изучения биологических объектов.

Существует несколько методов расчёта спектров малоуглового рассеяния нейтронов и рентгеновского излучения для молекулярных моделей, разработаны соответствующие компьютерные программы. В рамках аналитического подхода может быть использовано полишаровое представление объекта, расчёт спектров МУРН в таком случае осуществляется при помощи формулы Дебая. Другим приближением является представление однородного объекта в виде радиальной оболочки, которая аппроксимируется сферическими функциями, что позволяет вычислить спектр рассеяния. Кроме того, в ряде работ разработаны методы непосредственного вычисления функции распределения парных расстояний (ФРПР) и получения спектров МУРН с помощью Фурье-преобразования. Общим недостатком этих методов является их неприменимость для систем большого размера, таких как сложные био-макромолекулярные комплексы, состоящие из миллионов атомов.

Таким образом, разработка методов молекулярного моделирования наднуклеосомной структуры хроматина геномного размера (порядка 106 нуклеосом, или Ю10 атомов) и соответствующих методов расчёта спектров МУРН, применимых для систем такого размера, является актуальной биофизической задачей. Её решение позволит рассчитывать спектры МУРН и другие свойства моделей наднуклеосомной структуры хроматина для исследования влияния различных параметров на экспериментально наблюдаемые характеристики хроматина.

Цель и задачи исследования

Целью настоящей работы являлась разработка методов молекулярного моделирования наднуклеосомной структуры хроматина геномного размера и расчёта спектров МУРН для получаемых моделей нуклеопро-теидных комплексов.

Для достижения этой цели были поставлены следующие задачи: 1. разработать метод построения физической модели нуклеосомы и со-

ответствующий метод расчёта спектров МУРН, применимый для нук-леопротеидных комплексов геномного размера;

2. определить геометрические параметры межнуклеосомного интерфейса и разработать методы построения геометрических моделей над-нуклеосомных структур геномного размера;

3. разработать методы построения фрактальных наднуклеосомных структур с заданным значением фрактальной размерности;

4. рассчитать спектры МУРН для построенных моделей наднуклеосом-ной структуры хроматина и сопоставить с полученными экспериментально.

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

1. Разработанный метод расчёта, спектров МУРН, основанный на методе Монте-Карло, применим для нуклеопротеидных комплексов, состоящих из ~ 1010 атомов.

2. Разработанный алгоритм фрактальной генерации наднуклеосомных структур позволяет создавать фрактальные структуры с заданной фрактальной размерностью в интервале 2 < ¿3 <3.

3. Спектры МУРН, рассчитанные для фрактальных моделей паднук-леосомной структуры хроматина высоких порядков упаковки, имеющих две фрактальные размерности, удовлетворительно описывают экспериментально наблюдаемые спектры МУРН, а именно имеют степенной характер зависимости интенсивности рассеяния от величины вектора рассеяния (I ос в~с1') и точку кроссовера, определяющую переход от одной фрактальной размерности с?/ к другой.

Научная новизна

1. Разработан метод построения фрактальных моделей наднуклеосом-пой структуры хроматина с заданным значением фрактальной размерности или двумя различными значениями фрактальной размерности.

2. Разработан метод расчёта спектров МУРН, применимый для нуклеопротеидных комплексов геномного размера.

3. Построена фрактальная молекулярная модель наднуклеосомной структуры хроматина геномного размера 106 мононуклеосомпых частиц) с двумя фрактальными размерностями, описывающая двухуровневую фрактальную организацию хроматина.

Теоретическая и практическая значимость

Разработанные методы построения молекулярных моделей наднуклеосомной структуры хроматина могут быть использованы при расчёте

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

Разработанные методы расчёта спектров МУРН могут применяться как для моделей наднуклеосомной структуры хроматина, полученных в данной работе, так и для полноатомпых моделей любых биомакромоле-кулярных или полимерных комплексов с большим количеством атомов, где невозможно применение других методов. Апробация работы

Материалы диссертации были представлены на VIII Национальной конференции «Рентгеновское, синхротронное излучения, нейтроны и электроны для исследования наносистем и материалов» (Курчатовский институт, Москва, 2011), на XV Международной Пущинской школе-конференции молодых учёных «Биология — наука XXI века» (Пущино, 2011), на XI Международной конференции по рассеянию рентгеновского излучения и нейтронов на поверхностях SXNS-11 (Северо-Западный университет, Эванстон, Иллинойс, США, 2010), на XX и XXI Совещаниях по использованию рассеяния нейтронов в исследованиях конденсированного состояния (ПИЯФ, Гатчина, 2008, Курчатовский институт, Москва, 2010), на Санкт-Петербургском семинаре по компьютерной биологии (СПбГПУ, Санкт-Петербург, 2009), на IX конференции молодых учёных Отделения молекулярной и радиационной биофизики ПИЯФ РАН (ПИЯФ, Гатчина, 2008), на XVI и XVII Политехнических симпозиумах «Молодые учёные — промышленности Северо-Западного региона» (СПбГПУ, Санкт-Петербург, 2007, 2008).

Начальный этап работы (магистерская диссертация) был отмечен медалью «За лучшую научную студенческую работу» по итогам Всероссийского открытого конкурса на лучшую научную работу студентов вузов по естественным, техническим и гуманитарным наукам в 2008 г. Публикации и личный вклад автора

По материалам диссертации опубликовано 13 работ, в том числе 4 в рецензируемых научных журналах из перечня изданий, рекомендованных Высшей аттестационной комиссией для публикации материалов диссертаций на соискание учёной степени кандидата наук (полный список приведён на странице 18).

Автор принимал непосредственное участие на всех этапах работы. Личный вклад автора в получение результатов, выносимых на защи-

ту, является определяющим. Вклад автора в совместных публикациях с научным руководителем и/или другими сотрудниками лаборатории является основным (анализ литературы, разработка, апробация и применение методов, расчёты и анализ результатов, написание статей), за исключением публикаций [3, 6, 9], в которых автором была выполнена лишь часть исследования.

Объём и структура диссертации

Диссертационная работа изложена на 128 страницах машинописного текста и включает введение, три главы (обзор литературы, методы, результаты и обсуждение) и выводы. Материал иллюстрирован 26 рисунками и 7 таблицами. Библиографический указатель содержит 120 источников.

СОДЕРЖАНИЕ РАБОТЫ

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

1. Обзор литературы

Обобщены современные представления о структурной организации хроматина на различных уровнях компактизации ДНК. На основе анализа экспериментальных данных был сделан вывод о том, что только иуклеосомпый уровень организации исследован с удовлетворительной полнотой, в то время как строение 30-нанометровой нити хроматина остаётся открытым вопросом. Рассмотренные данные о структуре хроматина на уровне целого ядра говорят в пользу двухуровневой фрактальной организации хроматина. Описаны основные принципы МУРН, способы и практические методы вычисления спектров МУРН.

2. Методы

2.1 Геометрическая модель нуклеосомы

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

2.2 Анализ GC-состава сайтов образования нуклеосом

Разработанный алгоритм анализа GC-состава сайтов образования нуклеосом (СОН) в последовательности генома человека hgl8 состоит в построении глобального (для целых хромосом) распределения GC-состава и последующем сравнении локальных (для небольших участков) распределений GC-состава с глобальным при помощи U-теста Манна — Уитни для определения участков последовательности ДНК с аномальным (существенно отличающимся от среднего) GC-составом.

2.3 Геометрическая модель наднуклеосомной структуры

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

В режиме одноцепочечной генерации построение наднуклеосомной структуры начинается с первой нуклеосомы, имеющей фиксированные координаты. В работе были разработаны два режима одноцепочечной генерации: свободная генерация и генерация в объёме. В режиме свободной генерации все параметры интерфейса выбираются случайным образом, при этом длина линкера L распределена в интервале [Irnin! ЬтахЬ все углы — в интервале [0; 27г]. Режим генерации в объёме дополнительно накладывает ограничение в виде кубической или сферической стенки, за пределами которой генерация запрещена. Кроме того, в подрежиме локально-регулярной генерации возможно получать регулярные участки максимальной протяжённостью Nr нуклеосом с фиксированными значениями параметров межнуклеосомного интерфейса.

2.3.1 Режим фрактальной генерации

Одним из методов построения фракталов является преобразование самоподобия. Фрактал 0-го порядка (инициатор) представляет собой отрезок длиной Л, который на первом шаге построения заменяется на генератор — ломаную, состоящую из п сегментов длиной А// (фрактал 1-го порядка). Далее каждый сегмент генератора заменяется на целый генератор. Для идеального фрактала (из бесконечного числа сегментов)

фрактальная размерность <1 определяется как размерность подобия

А 1ПП (Л\

= ь? (1)

Если на каждом шаге генератор произвольно изменяется (при постоянных п и /), то получается случайный фрактал, статистические характеристики которого близки к свойствам идеального фрактала.

С помощью данного метода осуществляется построение фрактальной модели наднуклеосомпой структуры хроматина. Генерация заканчивается фракталом порядка к = К: структура состоит из N = пК точек, соединённых сегментами длиной I — \/}к. N моделей нуклеосом с эффективным размером гтах помещаются в декартовы координаты N точек структуры. Фрактальная размерность структуры может быть определена по формуле (1) или получена из ожидаемой зависимости среднего радиуса гирации гд от количества точек Аг3:

(гэ> сх [Мд]1/а' . (2)

2.4 Полишаровая физическая модель нуклеосомы

Аналитически спектр МУРН может быть рассчитан для совокупности однородных шаров с помощью формулы Дебая:

= Е Е МШ^-Щр1, ГУ = ¡ту -г, I, (3)

г~1 ¿ = 1

где {/¿(я)}^! — амплитуды рассеяния отдельных шаров, зависящие от их плотности длин рассеяния (ПДР) и объёма. Для расчётов с помощью формулы Дебая (3) необходимо построить полишаровую модель цукле-осомы. Для белков и ДНК известны средние значения ПДР. МУРН является низкоразрешающим 10 А) методом, что позволяет строить упрощённые модели. Исходя из этих данных была выбрана следующая схема построения модели нуклеосомы: гистоновый кор моделируется одним шаром, ДНК — цепочкой шаров, образующих спираль.

2.5 Полноатомная физическая модель нуклеосомы Автором была разработана оригинальная схема вычислений спектра

МУРН с использованием полноатомной структуры (рис. 1). На первом этапе происходит анализ структуры из банка РОВ (п. 2.1) и строятся списки атомов и остатков. Последний используется при проверке на самопересечения во время генерации наднуклеосомной структуры (п. 2.3).

Анализ геометрии |

Генерация наднуклеосомной структуры

Список нуклеосом ■

Двумерная гистограмма

Определение статистических весов ПР,

Одномерная гистограмма

ФРПРОЧу

<Сиектр МУРН)

Рис. 1. Блок-схема расчётов спектров МУРН с использованием полноатомной модели нуклеосомы. Символами I 1, I I и <-> обозначены

входные данные/параметры, вычислительные процессы, промежуточные и выходные данные соответственно

Главным элементом вычислений является расчёт ФРПР, основанный на методе Монте-Карло (МК). ФРПР в пределе учитывает все парные расстояния (ПР) между элементами структуры (ср. с формулой Дебая (3)). Статистическая природа ФРПР позволяет ограничиться некоторым (достаточно большим) числом ПР при сохранении адекватности отражения строения структуры. На каждой МК-итерации случайно выбирается пара атомов, вычисляется ПР между ними, его значение и типы остатков записываются в двумерную гистограмму.

В режиме вычисления ФРПР однородной частицы (ФРПРОЧ) статистический вес каждого ПР равен 0 или 1 в зависимости от выбранных остатков, т. е. возможно вычислить ФРПРОЧ для белковой компоненты или ДНК. В режиме вычисления спектра МУРН статистический вес ПР определяется типом остатков с учётом концентрации тяжёлой воды в растворе 1У0 и вероятности обмена атомов водорода на дейтерий Рц/в ■

«ч, = п А* = ^ + - Ь?)^ = Рн+и/р(№ _,н)>

1,3

где ЬЬР, VI — длина рассеяния в Н20 и Б20 и объём остатка г-го типа, рп и рс — ПДР Н20 и Б20. Спектр МУРН получается Фурье-преобразованием одномерной гистограммы.

2.6 Модель однородного шара в методе Монте-Карло

Для сравнения двух физических моделей нуклеосомы, полишаровой (п. 2.4) и полноатомной (п. 2.5), а также для выяснения необходимого количества МК-итераций был разработан способ расчёта ФРПР для полишаровой модели. Пусть множество точек Q равномерно распределено в однородном шаре с радиусом Rq. Тогда сферические координаты точки Qi определяются выражениями: г = Rq {/тГ, cos 0 = — 1, <р = 27Г7з, где {7fc}I=i — случайные величины, равномерно распределённые в интервале [0; 1]. При вычислениях ФРПР на каждом шаге но выбору атома происходит генерация точек внутри шаров, и эти координаты используются для вычисления ПР.

3. Результаты и обсуждение

3.1 Анализ геометрии нуклеосомных структур из банка PDB

Для оценки структурного разнообразия мононуклеосом операции по определению геометрических параметров (п. 2.1) были проведены для 33 структур из банка PDB. Оказалось, что структуры имеют весьма схожее строение, типичная нуклеосома состоит из 12121 атома, 146 и. н. ДНК образуют 1,8 витка суперспирали диаметром 84,0 А и шагом 25,8 А. Нуклеосома 1КХ5 была выбрана для дальнейших вычислений как репрезентативный представитель данного класса структур, при этом она имеет максимальные разрешение и количество атомов.

3.2 Анализ GC-состава СОН

Для расчёта распределения длин линкеров с помощью специального алгоритма (п. 2.2) был проведён анализ последовательностей СОН. Из последовательности генома были извлечены фрагменты длиной 900 п. н., содержащие 438 599 СОН, и проанализировано относительное содержание участков с аномальным GC-составом. Оказалось, что в последовательностях СОН имеется чёткий сигнал — изменение профиля GC-состава в сторону больших значений. Однако за границами собственно пуклеосом такие сигналы отсутствуют, что не позволяет использовать данный анализ для расчёта распределения длин линкеров.

3.3 Расчёты спектров МУРН с использованием полишаровой модели нуклеосомы

На основе геометрических параметров нуклеосомы 1КХ5 (п. 3.1) была построена полишаровая модель (п. 2.4). С использованием формулы Дебая (3) были рассчитаны спектры МУРН и ФРПР для полишаровой модели нуклеосомы и построенного на её основе ряда структур длиной

до 15 нуклеосом. Получаемые в рамках этого подхода данные самосогласованы и достаточны для исследования моделей наднуклеосомных структур небольшого размера. Однако метод обладает рядом недостатков, главный из которых — время счёта, что делает его непригодным для расчётов спектров МУРН для структур геномного размера.

Для решения этой проблемы был разработан оригинальный метод расчёта спектров МУРН, основанный на методе Монте-Карло (п. 2.5). С применением модели однородного шара (п. 2.6) для ряда структур были рассчитаны спектры МУРН и сопоставлены с полученными при использовании формулы Дебая (3). Анализ результатов позволил сформулировать критерии выбора параметров МК-моделирования ФРПР.

3.4 Расчёты спектров МУРН с использованием полноатомной модели нуклеосомы

3.4.1 Мононуклеосома

При построении полноатомной модели нуклеосомы (п. 2.5) были использованы 33 структуры из банка РОВ, после чего был проанализирован ряд свойств, важных для МУРН (ПДР, точки компенсации). Как и в случае геометрических параметров (п. 3.1), проанализированные структуры оказались весьма похожи и для дальнейших расчётов была выбрана модель на основе структуры 1КХ5.

3.4.2 Наднуклеосомные структуры, полученные в режиме од-ноцепочечной генерации

На основе модели нуклеосомной частицы в соответствии с выбранным режимом генерации (п. 2.3) были построены модели наднуклеосомных структур хроматина. Во всех режимах генерации к случайно растущей цепи нуклеосом применялся ряд ограничений. При свободной генерации длина линкера должна была лежать в определённом диапазоне, также были запрещены самопересечения структуры. Полученные в этом режиме структуры из 106 нуклеосом имеют характерные линейные размеры порядка 105 А и, следовательно, обладают очень низкой плотностью, что послужило причиной введения ограничивающей стенки, т. е. режима генерации в объёме.

Так как эффективный размер используемой модели мононуклеосо-мы составляет около 100 А, куб с ребром 104 А может содержать до ~ 106 нуклеосом. Действительно, в случае применения этой кубической стенки, удельное (в расчёте на одну нуклеосому) время генерации существенно возрастает по мере приближения к предельному значению.

Рис. 2. Свойства моделей наднуклеосомных структур, полученных в режиме одноцепочечной генерации. А. Спектры МУРН для двух моделей наднуклеосомных структур для указанных на графике значений И^д. Обе серии спектров приведены в относительных единицах (о. е.) для помещения на один график. Дополнительная шкала показывает размеры в А в прямом пространстве. В. Зависимость среднеквадратичного расстояния между концами цепи от длины цепи для структур из 106 нуклеосом. Горизонтальная линия указывает размер стенки. Штриховой линией показана линейная аппроксимация (ЛА) с коэффициентом углового наклона (КУН) 0,5

Типичные спектры МУРН для наднуклеосомных структур имеют нуклеосомную (в > 0,02 А-1) и наднуклеосомпую (в < 0,02 А"1) области (рис. 2А). В последней спектры независимы от ДОц, а в нуклеосом-иой части имеют ряд И^Ь-зависимых особенностей, аналогичных тем, что наблюдались для мононуклеосомы. Кроме того, в спектре белковой компоненты (при Ш-о — 0,65) имеется пик 0,06 А-1, обусловленный узким интервалом длин линкеров.

В наднуклеосомной области спектр МУРН является линейным с КУН —2,0 (рис. 2А), т.е. структура проявляет свойства фрактала размерности 2,0, что близко к экспериментально наблюдаемому спектру, полученному в 40 % Б20 (КУН -2,2). Однако экспериментальные спектры в 64 % и >99 % Б20 имеют отличные значения КУН, равные -2,5 и —2,4 соответственно. Это не может быть объяснено в рамках простой наднуклеосомной модели. Относительная интенсивность пика 0,06 А-1 выше в расчётных спектрах по сравнению с экспериментальными — | должен существовать неучтённый фактор, способный изменить ФРПР белковой компоненты структуры, например, ненуклеосомный белок.

Наднуклеосомные структуры, полученные в режиме свободной генерации, имеют чрезвычайно низкую плотность, порядка 10® нуклеосом

в кубе с ребром 105 А, тогда как ~ 2 х 109 и. н. полного генома курицы (« 107 нуклеосом) упакованы в ядре диаметром около 3 мкм. Для увеличения плотности структуры был введён режим генерации в объёме. Единственным отличием между рассчитанными спектрами МУРН для двух типов наднуклеосомных структур оказалась серия пиков, обусловленная форм-фактором ограничивающей стенки (рис. 2А). Более детальный анализ показал (рис. 2Б), что внутри ограниченного объёма свойства структур совпадают, отличия наблюдаются только на размерах, близких к размеру стенки. Таким образом, несмотря на существенно большую (на три порядка) среднюю плотность структуры, ограничивающая стенка является недостаточным фактором для изменения фрактальных свойств системы.

Для выяснения влияния регулярной упаковки нуклеосом на спектры МУРН были созданы локально-регулярные модели наднуклеосомной структуры из 106 нуклеосом с использованием регулярных структур, соответствующих моделям «зигзаг» и «соленоид», длиной до ЛГГ = 20 нуклеосом. Их спектры в наднуклеосомной области являются линейными с КУН —2, и только в нуклеосомной области наблюдается ряд особенностей, обусловленных наличием регулярных структур.

3.4.3 Наднуклеосомные структуры, полученные в режиме фрактальной генерации

С помощью алгоритма фрактальной генерации (п. 2.3.1) для ряда значений межнуклеосомного расстояния I и фрактальной размерности ¿3 была получена серия структур из 220 нуклеосом. Рассчитанные зависимости радиуса гирации {гд(А^)) и спектры МУРН имели степенной характер, что позволило оценить фрактальные размерности йд и ¿[ по формуле (2) и из выражения I ос Было установлено, что значение I = 300 А обеспечивает хорошее согласие всех трёх оценок фрактальной размерности: КУН корреляций <1д/с18, ¿¡/с1в и 3*1/<1д составили 0,987 (Л2 = 0,987), 1,007 (Л2 = 0,984) и 1,021 (Я2 = 0,997) соответственно.

График зависимости {гд{Мд)) для типичной фрактальной структуры с размерностью <1В = 2,4 приведён на рис. ЗА. В начале наблюдается относительно короткий участок (Ид < 10) с «КУН» свыше 0,5, обусловленный эффективным отталкиванием нуклеосом. За исключением этого участка зависимость (гд(Мд)) является линейной с КУН 0,42, который соответствует значению размерности <1д = 2,38, что весьма близко к ожидаемому значению с18 = 2,4.

Н

Рис. 3. Свойства двух фрактальных моделей наднуклеосомных структур из 220 нуклеосом, фрактальные размерности: с1я = 2,4 (А, Б) и £ {2; 3} (В, Г). Пунктирные линии показывают ЛА. А, В. Зависимость (гэ(Лгэ)): символами • обозначены измеренные значения; сплошными линиями показаны значения, соответствующие штриховые линии, соответствующие ¿д 6 {1; 2; З},1 приведены для сравнения. А. ЛА у = а + @х; КУН 0,42 (Л2 = 0,999) соответствует значению ¿д = 2,38. В. ТК по построению (отмечена стрелкой — 256 нуклеосом (1697 А). ЛА у = а + х/в,а (Л2 = 0,999) указывают на другую ТК (отмечена стрелкой 4-) — 584 нуклеосом (2933 А). Б, Г. Спектры МУРН: показаны сплошными линиями для указанных значений И7^■ В. ЛА у = а + /Зж: КУН -2,45 (Л2 = 0,999) соответствует значению й/ = 2,45. Г. ЛА у = а - 2х (Л2 = 0,998) и у = а - Зж (Л2 = 0,991)

Как и для структур, полученных в режиме одноцепочечной генерации, спектры МУРН для фрактальных моделей имеют две области: наднуклеосомную (в < 0,02 А"1) и нуклеосомную (рис. ЗБ). В над-нуклеосомной области спектры являются независимыми от значения \Уп; в нуклеосомной области, напротив, наблюдаются особенности для спектров белковой части = 0,65) и ДНК (И^о = 0,40), обусловленные локально-регулярной упаковкой нуклеосом, которая возникла вследствие работы алгоритма генерации. Наднуклеосомная часть спектров является линейной при 5 < 0,01 А-1 с КУН —2,45, что соответствует размерности = 2,45, при этом ¿в = 2,4 и <1д = 2,38.

Из экспериментальных спектров МУРН для целых ядер были получены оценки фрактальной размерности хроматина ¿¡, находящиеся в диапазоне от до С помощью других экспериментальных методов были сделаны независимые оценки фрактальной размерности хроматина: от 2,2 (диффузионные исследования) до 3 (метод НьС). Это послужило отправной точкой для генерации наднуклеосомной структуры с двумя различными фрактальными размерностями

Наднуклеосомная структура из 220 нуклеосом с фрактальными размерностями = 3 и с^г = 2 весьма похожа на структуру с фрактальной размерностью <1Б = 2,4. Анализ зависимости (Гд{Ыд)) (рис. ЗВ) выявил двухуровневую организацию структуры с точкой кроссовера (ТК) 584 нуклеосом (2933 А), которая смещена от ТК по построению (256 нуклеосом, 1697 А) на лишь немногим более одного порядка генерации.

Спектр МУРН для структуры с кроссовером (рис. ЗГ) схож со спектром для наднуклеосомной структуры с фрактальной размерностью ¿5 = 2,4 и состоит из наднуклеосомной части (в < 0,01 А-1) с I ос и нуклеосомной части с характерными особенностями, обусловленными локально-регулярной структурой. В наднуклеосомной области видны два участка, соответствующие фрактальным размерностям й/ я», ¿3.

Наднуклеосомная часть спектров МУРН, рассчитанных для фрактальных моделей наднуклеосомной структуры хроматина, особенно спектр для структуры с двумя различными значениями фрактальной размерности — 3 (для больших размеров) и 2 (для меньших размеров), находится в качественном согласии с теми спектрами МУРН, которые были получены в нашей лаборатории для целых ядер.

На рис. 4 приведён спектр МУРН, полученный для ядер лимфоцитов крысы при концентрации тяжёлой воды в растворе И^о = 0,4. Спектр имеет два линейных участка в двойном логарифмическом масштабе с КУН -2,9 и —2,1 и ТК 1,9 х 103 А (рис. 4А). Предварительные расчёты показали, что спектры, рассчитанные для фрактальной наднуклеосомной модели из 220 нуклеосом для гтах = 120 А, £ — 300 А, п = 2, кс = 17 с двумя фрактальными размерностями = 2,9 и £¿¡,2. = 2,1, имеют в наднуклеосомной области (в < 0,01 А-1) два линейных участка с ¿1 ~ <18 и близкой ТК (рис. 4Б) и удовлетворительно описывают экспериментальные данные (Д1 = 0,01, рис. 4В).

Данные диффузионных исследований хроматина (характерные размеры <100 нм) говорят в пользу низких значений фрактальной размерности (2,2—2,6), в то время как данные метода НьС (разрешение

Рис. 4. Сравнение спектров МУРН. Интенсивность в о. е., на одном графике показаны: ЛА экспериментальных данных (А), Л А расчётного спектра (Б), наложение спектров (В). Символами ■ показаны экспериментальные данные, сплошными линиями — расчётный спектр, штриховыми линиями — ЛА, стрелками Т отмечены ТК. А. ЛА у = а + /?х: КУН -2,90 ± 0,04 {Я2 = 0,990) и -2,10 ± 0,02 (Я2 = 0,992). ТК 1,9 х 103 А. Б. ЛА у = а + с1лх (Я2 = 0,999) и у = а + <1лх (Я2 = 0,997). ТК 2,1 х 103 А. В. Нормировка выбрана так, чтобы минимизировать значение

~ 5 х 103 нуклеосом) описываются моделью фрактальной глобулы (размерность 3). Таким образом, единая фрактальная модель надиукле-осомной структуры хроматина с кроссовером способна описать данные, полученные с помощью трёх независимых экспериментальных методов.

I

Основные результаты и выводы

1. Разработан новый метод расчёта спектров МУРН, основанный на МК-моделировании распределения парных расстояний.

2. Установлено, что разработанный метод расчёта спектров МУРН применим для систем, состоящих из ~ Ю10 атомов.

3. Разработан алгоритм генерации фрактальных структур на основе преобразования самоподобия.

4. Показано, что разработанный алгоритм фрактальной генерации над-нуклеосомных структур позволяет создавать фрактальные структуры с заданным значением фрактальной размерности в интервале 2 < <1В < 3 или с двумя различными значениями фрактальной размерности (структуры с кроссовером).

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

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

Список опубликованных работ по теме диссертации

1. Ilatovskiy A. et al. SANS Spectra of the Fractal Supernucleosomal Chromatin Structure Models // J. Physi: Conf. Ser. 2012. V. 351. P. 012007.

2. Ilatovskiy A. et al. Modeling and SANS Spectra of Chromatin Supernucleosomal Structures at Genome Scale // J. Appl. Phys. 2011. V. 110. P. 102217.

3. Руденко E. u др. Анализ молекулярных механизмов растяжения коротких фрагментов двухнитевых ДНК // НТВ СПбГПУ. 2011. Т. 134. С. 147-154.

4. Ilatovskiy A., Petukhov М. Genome-Wide Search for Local DNA Segments with Anomalous GC-Content // J. Comput. Biol. 2009. V. 16. P. 555-564.

5. Платовский А. и др. Спектры MPH для локально-регулярных моделей супернуклеосомных структур хроматина хромосомного размера // Материалы VIII Национальной конференции РСНЭ-НБИК. М., 2011. С. 445.

6. Лебедев Д. и др. Использование методов нейтронного рассеяния в изучении структуры и динамики нуклеопротеидных комплексов // Материалы VIII Национальной конференции РСНЭ-НБИК. М., 2011. С. 31.

7. Платовский А. и др. Расчёт спектров МРН для моделей супернуклеосом-ной структуры хроматина // Материалы 15-й Международной Пущинской школы-конференции молодых учёных. Пущино, 2011. С. 118.

8. Ilatovskiy A. et al. SANS Spectra and Modeling of Chromatin Supernucleosomal Structures at Genome Scale // Proceedings of the 11th Intl. Conference on Surface X-ray and Neutron Scattering. USA, Evanston, 2010. PO-A-45.

9. Ilatovskiy A. et al. Molecular Dynamics and Large-Scale Conformational Flexibility of DNA-Binding Proteins // Proceedings of the 11th Intl. Conference on Surface X-ray and Neutron Scattering. USA, Evanston, 2010. PO-A-46.

10. Илатовский А. Моделирование супернуклеосомной структуры хроматина // Материалы 13-й Санкт-Петербургской ассамблеи молодых учёных и специалистов. СПб., 2008. С. 37.

11. Илатовский А., Петухов М. Расчёт спектров МРН для супернуклеосомной структуры хроматина // Материалы конференций политехнического симпозиума. СПб., 2008. С. 115-116.

12. Ilatovskiy A. et al. Modelling of the Supernucleosomal Chromatin Structure // Proceedings of XX Workshop on Neutron Scattering Application for Condensed Matter Investigations. Russia, Gatchina, 2008. P. 165.

13. Илатовский А., Петухов M. Моделирование супернуклеосомной структуры хроматина // Материалы конференций политехнического симпозиума. СПб., 2007. С. 119-120.

Работа выполнена при финансовой поддержке РАН (проект 4.7), Правительства Санкт-Петербурга (ПГ № 2.6/25-04/05), Министерства образования и иауки РФ (РНП 2.2.1.1.4663, ГК № 02.740.11.5223, № 11.519.11.2002) с использованием научного оборудования ЦКП «Аналитический цептр нано- и биотехнологий ГОУ СПбГПУ».

Отпечатано в типографии ФГБУ «ПИЯФ»

188300, Гатчина Ленинградской обл., Орлова роща Зак. 92, тир. 100, уч.-изд. л. 1; 13.04.2012 г.

Текст научной работыДиссертация по биологии, кандидата физико-математических наук, Илатовский, Андрей Владимирович, Санкт-Петербург

61 12-1/1057

ФГБУ «Петербургский институт ядерной физики им. Б.П.Константинова»

ИЛАТОВСКИЙ Андрей Владимирович

МЕТОДЫ МОДЕЛИРОВАНИЯ НАДНУКЛЕОСОМНОЙ СТРУКТУРЫ ХРОМАТИНА

И РАСЧЁТА СПЕКТРОВ МАЛОУГЛОВОГО РАССЕЯНИЯ НЕЙТРОНОВ

Диссертация на соискание учёной степени кандидата физико-математических наук

03.01.02 — биофизика

Научный руководитель: кандидат физико-математических наук Исаев-Иванов Владимир Васильевич

Санкт-Петербург — 2012

И си

Содержание

Список сокращений 4

Введение 5

1 Обзор литературы 12

1.1 Структура нуклеосомы .......................12

1.2 Компактизация хроматина.....................18

1.3 Структура нити хроматина.....................21

1.4 Структура хроматина на уровне ядра...............27

1.5 Малоугловое рассеяние нейтронов.................31

2 Методы 3-4

2.1 Геометрическая модель нуклеосомы................34

2.2 Алгоритм анализа СС-состава сайтов образования нуклеосом . 4-1

2.2.1 Определение последовательностей

сайтов образования нуклеосом...............41

2.2.2 Биномиальный тест .....................4^

2.2.3 и-тест Манна — Уитни ...................43

2.3 Геометрическая модель наднуклеосомной структуры.......44

2.3.1 Параметризация межнуклеосомного интерфейса.....44

2.3.2 Режим одноцепочечной генерации

наднуклеосомной структуры ................48

2.3.2.1 Подрежим локально-регулярной генерации ... 49

2.3.3 Режим фрактальной генерации . ..............50

2.4 Полишаровая физическая модель нуклеосомы..........52

2.5 Полноатомная физическая модель нуклеосомы..........57

2.5.1 Построение списка атомов нуклеосомы ..........59

2.5.1.1 Регуляризированная модель нуклеосомы .... 60

2.5.2 Построение списка остатков нуклеосомы .........60

2.5.3 Монте-Карло моделирование распределения парных расстояний и двумерная гистограмма...........61

2.5.4 Параметризация плотности длин рассеяния........62

2.5.5 Одномерная гистограмма и расчёт спектров малоуглового рассеяния нейтронов.............6-4

2.6 Модель однородного шара в методе Монте-Карло........65

3 Результаты и обсуждение 70

3.1 Анализ геометрии мононуклеосомных структур из банка РБВ . 70

3.2 Анализ СС-состава сайтов образования нуклеосом........75

3.3 Расчёты спектров малоуглового рассеяния нейтронов

с использованием полишаровой модели нуклеосомы.......78

3.3.1 Расчёты с помощью формулы Дебая............78

3.3.1.1 Мононуклеосома..................79

3.3.1.2 Наднуклеосомные структуры...........81

3.3.1.3 Анализ применения формулы Дебая.......85

3.3.2 Переход к методу Монте-Карло...............85

3.3.2.1 Мононуклеосома..................86

3.3.2.2 Наднуклеосомные структуры...........88

3.3.2.3 Анализ применения метода Монте-Карло .... 88

3.4 Расчёты спектров малоуглового рассеяния нейтронов

с использованием полноатомной модели нуклеосомы.......90

3.4.1 Мононуклеосома.......................90

3.4.2 Наднуклеосомные структуры, полученные в режиме од-

ноцепочечной генерации...................98

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

подрежиме локально-регулярной генерации . . .102

3.4.3 Наднуклеосомные структуры, полученные в режиме фрактальной генерации...................103

Основные результаты и выводы 113

Список опубликованных работ по теме диссертации 1Ц

Литература 116

Список сокращений

ACM атомно-силовая микроскопия

ДНК дезоксирибонуклеиновая кислота

днДНК двунитевая ДНК

КУН коэффициент углового наклона

МК Монте-Карло

МУРН малоугловое рассеяние нейтронов

о. е. относительные единицы

ПДР плотность длин рассеяния

п. н. пара нуклеотидов

ПР парное расстояние

ФРПР функция распределения ПР

ФРПРОЧ ФРПР однородной частицы

ЭМ электронная микроскопия

FRAP fluorescence recovery after photobleaching, восстановление флуоресценции после фотовыжигания RMSD root mean square deviation, среднеквадратичное отклонение; для двух векторов и {х2,i}i=i определяется как

Введение

Актуальность исследования

Современный арсенал экспериментальных биофизических методов исследования позволяет отвечать на многие вопросы о функционировании живой клетки. Тем не менее, наши знания о структурной организации биологических объектов по-прежнему недостаточны. Например, длина ДНК генома человека составляет около 2 м, в то время как характерный размер клеточных ядер — порядка 10 мкм, что означает весьма высокую степень упаковки ДНК.

Ядерная ДНК в геномах эукариот представляет собой нуклеопротеиновый комплекс — хроматин. Для компактного хранения генетической информации и обеспечения доступа к ней требуется его определённая организация. В середине 1970-х гг. было известно, что хроматин состоит из ДНК и гистонов HI, Н2А, Н2В, НЗ и Н4 (Kornberg and Thomas, 1974). На основании имеющихся биохимических данных Корнбергом (Kornberg, 1974) была выдвинута гипотеза о существовании структурного элемента хроматина, которая позднее была подтверждена экспериментально (Oudet et al., 1975). Обнаруженные частицы назвали нуклеосомами. В 1997 г. была получена полноатомная структура нуклеосомы (Luger et al., 1997).

Достигаемая с помощью нуклеосомной организации степень компак-тизации хроматина является недостаточной для объяснения наблюдаемой плотности упаковки ДНК в ядре, что указывает на наличие наднуклеосомных структур хроматина (Woodcock and Ghosh, 2010). Наблюдаемые при помощи электронной микроскопии 30-нанометровые нити хроматина считаются следующим уровнем компактизации ДНК (Oudet et al., 1975; Finch and Klug, 1976). За прошедшие десятилетия нити хроматина исследовались

множеством экспериментальных методов, однако полученные результаты не позволяют сделать однозначных выводов об их структуре, а ряд исследователей ставит под сомнение само существование регулярной упаковки нуклеосом в 30-нанометровые нити (van Holde and Zlatanova, 1995; Tremethick, 2007; Fussner et al., 2011).

Определённый прогресс достигнут в исследовании структуры хроматина в масштабе целого ядра. Применение биофизических и биохимических методов позволило установить двухуровневую фрактальную организацию хроматина (McNally and Mazza, 2010) и, в ряде случаев, построить трёхмерные модели исследуемых геномов (Duan et al., 2010). Тем не менее единая молекулярная модель структурной организации хроматина, охватывающая все уровни компактизации ДНК, от нуклеосомного до геномного, пока не создана.

В нашей лаборатории получены данные о малоугловом рассеянии нейтронов (МУРН) на различных клеточных ядрах (Lebedev et al., 2005; Исаев-Иванов и dp., 2010). Графики экспериментальных спектров МУРН имеют линейные участки в двойном логарифмическом масштабе для широкого интервала значений величины векторов рассеяния, что указывает на фрактальную организацию хроматина (Schmidt, 1989).

МУРН обладает рядом преимуществ перед другими экспериментальными методами исследования наднуклеосомной структуры хроматина. Во-первых, он позволяет получить данные о структуре хроматина в нативном ядре без применения биохимических техник выделения отдельных нитей хроматина. Во-вторых, при использовании метода вариации контраста возможно наблюдать рассеяние нейтронов на различных частях нуклеопротеидных комплексов, например, при концентрации тяжёлой воды в растворе порядка 65 % спектр МУРН будет обусловлен преимущественно белковой компонентой. В-третьих, спектры МУРН охватывают практически полный интервал размеров наднуклеосомных структур хроматина, от нуклеосомных до геномных.

Указанные преимущества использования МУРН для исследования структуры хроматина и других биомакромолекулярных комплексов обуславливают большое число нейтронных экспериментов для изучения биологических объектов (Harroun et al, 2006).

Существует несколько методов расчёта спектров малоуглового рассеяния нейтронов и рентгеновского излучения для молекулярных моделей, разработаны соответствующие компьютерные программы.

В рамках аналитического подхода может быть использовано полишаровое представление объекта, расчёт спектров МУРН в таком случае осуществляется при помощи формулы Дебая (Debye, 1915; Свергун и Фейгин, 1986). Такой метод расчёта спектров рассеяния применяется в ряде программ: DALAI_GA (Chacón et al., 1998, 2000), SAXS3D (Walther et al., 2000).

Другим приближением является представление однородного объекта в виде радиальной оболочки, которая аппроксимируется сферическими функциями, что позволяет непосредственно вычислить спектр рассеяния (Stuhrmann, 1970). Этот подход используются в программе SASHA (Svergun et al., 1996, 1997). Аналогично приближению радиальной оболочки, в программах DAMMIN, CRYSOL и CRYSON используется представление форм-факторов отдельных атомов в виде разложения по сферическим функциям, что позволяет вычислять спектры рассеяния для полноатомных моделей (Svergun et al., 1995, 1998; Svergun, 1999).

Кроме того, в ряде работ разработаны методы непосредственного вычисления функции распределения парных расстояний (ФРПР) и получения спектров МУРН с помощью Фурье-преобразования (Hammermann et al., 2000).

Общим недостатком этих методов является их неприменимость для систем большого размера, таких как сложные биомакромолекулярные комплексы, состоящие из миллионов атомов.

Таким образом, разработка методов молекулярного моделирования наднуклеосомной структуры хроматина геномного размера (порядка 106 нуклеосом или Ю10 атомов) и соответствующих методов расчёта спектров МУРН, применимых для систем такого размера, является актуальной биофизической задачей. Её решение позволит рассчитывать спектры МУРН и другие свойства моделей наднуклеосомной структуры хроматина для исследования влияния различных геометрических и физических параметров на экспериментально наблюдаемые характеристики хроматина.

Цель и задачи исследования

Целью настоящей работы являлась разработка методов молекулярного моделирования наднуклеосомной структуры хроматина геномного размера и расчёта спектров МУРН для получаемых моделей нуклеопротеидных комплексов.

Для достижения этой цели были поставлены следующие задачи:

1. разработать метод построения физической модели нуклеосомы и соответствующий метод расчёта спектров МУРН, применимый для нуклеопротеидных комплексов геномного размера;

2. определить геометрические параметры межнуклеосомного интерфейса и разработать методы построения геометрических моделей наднуклеосомных структур геномного размера;

3. разработать методы построения фрактальных наднуклеосомных структур с заданным значением фрактальной размерности;

4. рассчитать спектры МУРН для построенных моделей наднуклеосомной структуры хроматина и сопоставить с полученными экспериментально.

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

1. Разработанный метод расчёта спектров МУРН, основанный на методе Монте-Карло, применим для нуклеопротеидных комплексов, состоящих из ~ Ю10 атомов.

2. Разработанный алгоритм фрактальной генерации наднуклеосомных структур позволяет создавать фрактальные структуры с заданной фрактальной размерностью в интервале 2 < < 3.

3. Спектры МУРН, рассчитанные для фрактальных моделей наднуклеосомной структуры хроматина высоких порядков упаковки, имеющих две фрактальные размерности, удовлетворительно описывают экспериментально наблюдаемые спектры МУРН, а именно имеют степенной характер зависимости интенсивности рассеяния от величины вектора рассеяния (I ос

и точку кроссовера, определяющую переход от одной фрактальной размерности (¿/ к другой.

Научная новизна

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

2. Разработан метод расчёта спектров МУРН для молекулярных моделей, применимый для нуклеопротеидных комплексов геномного размера (порядка 106 нуклеосом « Ю10 атомов).

3. Построена фрактальная молекулярная модель наднуклеосомной структуры хроматина геномного размера 106 мононуклеосомных частиц) с двумя фрактальными размерностями, описывающая двухуровневую фрактальную организацию хроматина.

Теоретическая и практическая значимость

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

Разработанные методы расчёта спектров МУРН могут применяться как для моделей наднуклеосомной структуры хроматина, полученных в данной работе, так и для вычисления спектров МУРН для полноатомных моделей любых биомакромолекулярных или полимерных комплексов с большим количеством атомов, где невозможно применение других методов.

Апробация работы

Материалы диссертации были представлены на VIII Национальной конференции «Рентгеновское, синхротронное излучения, нейтроны и электроны для исследования наносистем и материалов» РСНЭ-НБИК

2011 (Курчатовский институт, Москва, 2011), на XV Международной Пущинской школ е-конференции молодых учёных «Биология — наука XXI века» (Пущино, 2011), на XI Международной конференции по рассеянию рентгеновского излучения и нейтронов на поверхностях БХИБ-П (СевероЗападный университет, Эванстон, Иллинойс, США, 2010), на XX и XXI Совещаниях по использованию рассеяния нейтронов в исследованиях конденсированного состояния РНИКС-2008 и РНИКС-2010 (ПИЯФ, Гатчина, 2008, Курчатовский институт, Москва, 2010), на Санкт-Петербургском семинаре по компьютерной биологии (СПбГПУ, Санкт-Петербург, 2009), на IX конференции молодых учёных Отделения молекулярной и радиационной биофизики ПИЯФ РАН (ПИЯФ, Гатчина, 2008), на XVI и XVII Политехнических симпозиумах «Молодые учёные — промышленности Северо-Западного региона» (СПбГПУ, Санкт-Петербург, 2007, 2008).

Начальный этап работы (диссертация на соискание степени магистра) был отмечен медалью «За лучшую научную студенческую работу» по итогам Всероссийского открытого конкурса на лучшую научную работу студентов вузов по естественным, техническим и гуманитарным наукам в 2008 г.

Публикации и личный вклад автора

По материалам диссертации опубликовано 13 работ, в том числе 4 в рецензируемых научных журналах из перечня изданий, рекомендованных Высшей аттестационной комиссией для публикации материалов диссертаций на соискание учёной степени кандидата наук (полный список приведён на страницах 1Ц—115).

Автор принимал непосредственное участие как в постановке задач, так и в их решении и обсуждении полученных результатов. Личный вклад автора в получение результатов, выносимых на защиту, является определяющим. Вклад автора в совместных публикациях с научным руководителем и/или другими сотрудниками лаборатории является основным (анализ литературы, разработка, апробация и применение методов, расчёты и анализ результатов, написание статей), за исключением публикаций [3, 6, 9], в которых автором была выполнена лишь часть исследования.

Объём и структура диссертации

Диссертационная работа изложена на 128 страницах машинописного текста и включает введение, три главы (обзор литературы, методы, результаты и обсуждение) и выводы. Материал иллюстрирован 26 рисунками и 7 таблицами. Библиографический указатель содержит 120 источников.

Финансовая поддержка

Работа выполнена при финансовой поддержке Российской академии наук (комплексная программа Президиума РАН «Физика элементарных частиц», подпрограмма «Нейтронная физика», направление «Исследование структуры, динамики и неординарных свойств вещества нейтронными методами», проект 4.7: «Изучение структуры и динамики молекулярных и надмолекулярных биологических комплексов методами рассеяния нейтронов»), Правительства Санкт-Петербурга (персональный грант № 2.6/25-04/05), Министерства образования и науки Российской Федерации (РНП 2.2.1.1.4663, Государственные контракты № 02.740.11.5223 и № 11.519.11.2002) с использованием научного оборудования ЦКП «Аналитический центр нано- и биотехнологий ГОУ СПбГПУ».

1. Обзор литературы

1.1 Структура нуклеосомы

Одни из первых доказательств существования структурного элемента хроматина были получены с помощью экспериментов по электрофорети-ческому разделению смесей гистоновых белков и исследования спектров рассеяния рентгеновских лучей на смесях днДНК и гистонов (Romberg and Thomas, 1974; Kornberg, 1974). На электрофореграммах смеси гистонов НЗ + Н4 наблюдалось 8 полос, сохраняющихся в широком интервале концентраций. После соотнесения молекулярных масс и подвижности частиц стало ясно, что в растворе присутствует гетеротетрамер (НЗ)2(Н4)2 и продукты его диссоциации вплоть до отдельных гистоновых белков. Аналогичным образом было установлено, что гистоны Н2А и Н2В образуют гетеродимер. Следующим этапом исследования стало получение спектров рассеяния рентгеновского излучения для растворов днДНК в присутствии тех или иных гистонов и их сравнение со спектрами для образцов хроматина, выделенных из клеток (Kornberg and Thomas, 1974). Оказалось, что для воспроизведения спектров (и, соответственно, структуры) хроматина необходимо присутствие четырёх гистонов Н2А, Н2В, НЗ и Н4, но не HI.

На основании имевшихся цитологических,