
- •Содержание
- •Глава 1. Обзор литературы 11
- •Глава 2. Построение модели. 39
- •Глава 3. Результаты численного моделирования. Активность одиночного RyR-канала при стационарных условиях 79
- •Введение
- •Глава 1. Обзор литературы
- •1.1 Механизмы сокращения клеток сердечной мышцы
- •1.2 Рианодиновый рецептор – основной элемент управления кальциевой динамикой в клетке
- •1.3 Эксперименты по изучению изолированных RyR-каналов
- •1.4 Модели функционированияRyR-каналов
- •Стохастическая динамика и электронно-конформационные взаимодействия в белках
- •1.7 Модели «общего пула»
- •1.8. Теория локального контроля
- •1.9 Моделирование активности клеток водителей сердечного ритма
- •1.9.1 Современные представления об авторитмической активности пейсмейкеров
- •1.9.3 Модель Мальцева-Лакатты
- •Глава 2. Построение модели.
- •2.1 Электронно-конформационная модель RyR-канала
- •2.1.1 Гамильтониан канала
- •2.1.2. Конформационный потенциал
- •2.1.3 Влияние уровняtrans[Ca] на форму конформационного потенциала RyR-канала
- •2.1.4. Структурные изменения канала в электронно-конформационной модели
- •2.1.5 Динамика конформационной координаты
- •2.1.6 Динамика электронной степени свободы
- •2.1.7 Инактивационое состояние RyR-канала
- •2.1.9 Эффекты туннелирования
- •2.1.10 Проницаемость RyR-канала
- •2.2.1 Электронно-конформационная модель решетки RyR-каналов
- •2.2.1.1 Гамильтониан решетки RyR-каналов
- •2.2.2 Схема динамики RyR-каналов в решетке высвобождающей единицы
- •2.2.3 Сопряжение динамики RyR-каналов с динамикой кальция в отделах высвобождающей единицы
- •2.3 Методы численной реализации модели
- •2.3.1 Метод Эйлера-Марайамы
- •2.3.2 Реализация электронных и туннельных переходов. Метод Монте-Карло
- •2.3.3 Численная схема для эк-модели RyR-канала
- •2.4 Описание программного комплекса
- •2.5 Заключение
- •Глава 3. Результаты численного моделирования. Активность одиночного RyR-канала при стационарных условиях
- •3.1 Анализ временных зависимостей конформационной координатыQ
- •3.2 Медленная конформационная динамика RyR-канала
- •3.2.1 Параметр эффективного трения г. Конформационная динамика RyR-канала
- •3.2.2 Влияние коэффициента упругости каналаK на форму конформационного потенциала
- •3.2.3 Зависимость конформационного потенциала от параметра электронно-конформационного взаимодействияа
- •3.3 Стохастическая динамика RyR-канала. Быстрые переходы
- •3.3.1 Кинетические характеристики динамики RyR-канала
- •3.3.2 Зависимость вероятности электронных переходов отcis[Ca]
- •3.4 Активация одиночного канала
- •3.5 Исследование процесса закрытия RyR-канала
- •3.6 Процесс адаптации RyR-каналов к продолжительной стимуляции
- •3.7 Динамика одиночного RyR-канала при установившемся уровне cis[Ca]
- •3.7.1 Зависимость активности RyR-канала от времени
- •3.7.2 Зависимость активности RyR-канала от уровня cis[Ca]
- •3.8 Заключение
- •4.1 Анализ модели высвобождающей единицы
- •4.1.1 Процессы открытия и закрытия каналов в высвобождающих единицах.
- •4.1.2 Анализ кооперативной динамики RyR-каналов в кластере
- •4.2.1 Высвобождающая единица как самоподдерживающийся кальциевый осциллятор
- •4.2.3 Влияние взаимодействия междуRyR-каналами на стабильность осцилляций системы
- •4.2.3 Эффект случайной остановки автоколебаний
- •4.2.3.1 Форма и устойчивость кластеров открытых каналов
- •4.2.3.2 Характерное время перехода в стационарное состояние
- •4.3 Заключение
- •Заключение
- •Список литературы
- •Основные публикации по теме диссертации
1.2 Рианодиновый рецептор – основной элемент управления кальциевой динамикой в клетке
RyR-канал является гигантским макромолекулярным комплексом (рис. 1.4), тетрамером, объединяющим четыре субъединицы (полипептидные цепи) с молекулярной массой 567 кДа каждая [2, 13]. Активность RyR модулируется под действием растительного алкалоида рианодина из коры Ryania speciosa, что и определило его название «рианодиновый рецептор». Существуют три изоформы RyR-каналов: RyR1 (клетки скелетной мускулатуры и мозжечка), RyR2 (клетки сердечной мышцы и мозга), RyR3 (клетки гладкой мускулатуры). Они различаются по своей структуре и динамическим параметрам.
Как
уже было сказано выше,RyR-каналы
располагаются на мембране СР; при
переходе в открытое состояние через
них происходит высвобождение ионов
Са2+
из просвета ТЦ (люмена).
Согласно данным электронной микроскопии
[12], RyR-канал имеет форму четырехлистника
со стороной 27 нм (рис. 1.4).
В процессе открытия RyR-канала изменяется его конформационное (структурное) состояние: данный белок изменяется по форме и структуре, переходя в проводящее состояние (рис. 1.5).
На рисунке 1.6 приведены результаты экспериментальных исследований структуры мембраны СР [13]; видно, что на мембранах RyR-каналы объединены в кластеры (решетки) (от 100 до 300 каналов в каждой). Каналы в решетке упорядочены и находятся под определенным углом друг к другу, также экспериментальные данные многих исследовательских групп говорят о существовании взаимодействия между каналами в кластере благодаря наличию между ними специфических связывающих белков [12-15].
Из всего вышесказанного можно сделать следующие выводы:
RyR-каналы имеют сложную молекулярную структуру и, соответственно, огромное количество степеней свободы.
На мембране СР каналы образуют регулярную структуру решетки. Моделирование динамики RyR-каналов усложняется тем, что при разработке теории кластера RyR-каналов необходимо учесть взаимодействия между ними.
1.3 Эксперименты по изучению изолированных RyR-каналов
Изучая динамические особенности изолированных RyR-каналов независимо друг от друга несколько авторов исследовало поведение RyR-каналов в лабораторных условиях в липидных бислоях при постоянных концентрациях ионов Са2+ в растворе в условиях отсутствия кальциевой стимуляции [17-19]. Схема типичной установки представлена на рис. 1.7. В экспериментах раствор разделяется на две части липидным бислоем, имитирующим клеточную мембрану, в котором закрепляется RyR-канал. Та часть, которая соответствует люменальной стороне канала называется trans, цитозольной стороне – cis.
В работах [17-19] помощью подобных установок была проведена серия экспериментов по регистрации ионных токов, проходящих через одиночный канал при различных концентрациях Са2+ в растворе. Отсутствие тока указывало на нахождении канала в закрытом состоянии, а их присутствие – на его открытие.
Первым важнейшей
особенностью RyR-канала как нанообъекта
является стохастический
характер процессов открытия/закрытия
при постоянных значениях концентраций
и
[17-19].
Важной кинетической
характеристикой поведения RyR-канала
является вероятность пребывания канала
в открытом состоянии (),
равная отношению суммы времен пребывания
в открытом состоянии к длительности
эксперимента:
.
При разбиении сигналов на небольшие
промежутки времениТ<500
мс в работе [19] исследовалась временная
зависимость вероятности пребывания
канала в открытом состоянии. Исследования
показали неоднородность функции
(рис.1.8).
График зависимости
активности канала от времени можно
разбить на участки (так называемые пучки
(bursts))
с различной вероятностью нахождения
канала в открытом состоянии. Данный
эффект получил название модовой
проводимости
(Modal
Gating)
[21-23]. Так в работе [21] выявлены следующие
моды проводимости RyR-канала: H
(high
)
– мода высокой активности,
>0.1,L
(low
)
– мода низкой активности, 0<
<0.1,I
(inactivated
mode)
– мода нулевой активности,
=0
(рис. 1.9).
Исследование
зависимости вероятности пребывания
канала в открытом состоянии от значения
концентрации Са2+
в trans-части
показало [15, 24, 25], что в ответ на повышение
увеличивалась частота открытий канала
и вероятность
(рис. 1.10). Из данного экспериментального
факта можно сделать непосредственный
вывод, что вероятность открытия RyR-канала
увеличивается с ростом концентрации
Са2+
в люмене.
В ряде работ [15,
21, 24], исследовавших влияние концентрации
на активность RyR-канала, показано, что
при повышении концентрации Са2+
в цитозоле (< 10 мкM)
наблюдается резкое повышение частоты
открытий канала и, соответственно,
увеличение вероятности открытия
одиночного канала. Далее, при достижении
определенного значения концентрации
(> 100 мкM)
происходит спад
,
то есть наблюдается процесс инактивации
RyR-канала при повышенной концентрации
цитозольного кальция (рис.1.11).
Уменьшение значения
при высоких значениях cis[Ca] объясняется
многими исследователями [12, 15, 21, 26] тем
фактом, что, являясь Са2+-активируемым
рецептором, RyR-канал имеет на cis-стороне
два типа активных центров, способных
связаться с ионами Са2+:
активационный и инактивационный [26-28].
На рисунке 1.12 изображена классическая
схема активации/инактивации канала
[27].
Сложная регуляция
RyR-канала кальцием впервые была обнаружена
в работе [29] при проведении экспериментов
в плоском липидном бислое. Показано,
что резкое увеличение cis[Ca]
(~ 1 мс) сначала быстро активирует одиночный
RyR-канал до высокого уровня вероятности
,
а затем значение
медленно уменьшается со временем при
фиксированной концентрации Са2+
(рис. 1.13). Этот эффект получил название
адаптации
канала к продолжительной стимуляции.
Следует отметить,
что в экспериментальных работах повышение
уровня cis[Ca] проводилось двумя способами:
резкое повышение с помощью лазерного
флэш-фотолиза [29-33] и постепенное повышение
в растворе вблизи RyR-канала [34, 35]. Второй
способ подразумевает градуальное
повышение уровня
вблизи активных частей канала вследствие
диффузионных процессов.
Результаты
экспериментов [35] по наблюдению эффекта
адаптации при градуальном повышении
говорят о том, что только в 30% реализаций
обнаруживался эффект адаптации, в
остальных 70% случаях понижения вероятности
пребывания канала в открытом состоянии
со временем не происходило (рис. 1.14). На
основании этих фактов был сделан вывод:
понижение
со временем зависит от скорости увеличения
уровня
.
На основе приведенных данных можно сделать следующие выводы:
Результаты экспериментальных исследований указывают на стохастический характер динамики RyR-каналов при постоянных значениях
и
.
Зависимость вероятности
от уровня
имеет колоколообразный вид (при
(0, 100) мкМ
возрастает, при
>100 мкМ
уменьшается).
Вероятность открытия RyR-каналов прямо пропорциональна концентрации Са2+ в люмене.
При резком повышении значении
наблюдается эффект сначала резкого повышения
,а затем медленного спада значения
с течением времени до стационарного уровня – адаптация к продолжительной стимуляции.