- •Содержание
- •Глава 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 Заключение
- •Заключение
- •Список литературы
- •Основные публикации по теме диссертации
2.1.4. Структурные изменения канала в электронно-конформационной модели
На основе представлений ЭК-модели можно рассмотреть возможные варианты процессов открытия и закрытия канала при постоянстве формы КП, связанные с взаимодействием ионов Са2+ с активационным центром канала и последующим изменением конформационного состояния. Простейшие схемы открытия и закрытия канала показаны на рисунке 2.3.

Схема 1. Открытие канала. Вначале канал находится в закрытом состоянии на левой ветви конформационного потенциала в минимуме (С), при этом имеют место следующие изменения степеней свободы RyR-канала:
1 – ионы Са2+ связываются с активационным центром канала, происходит переход на более высокий энергетический уровень, то есть изменение электронной степени свободы и переход на правую ветвь КП.
2 – происходит изменение конформационной координаты Q: релаксация канала со временем в минимум правой ветви КП (О). Соответственно, канал оказывается в электронно и конформационно открытом состоянии.
Схема 2. Закрытие
канала.
В начальный
момент времени канал находится в открытом
состоянии в минимуме (О)
на ветви КП
.
Наблюдаются следующие изменения степеней
свободы RyR-канала:
1 – ионы Са2+
отсоединяются от активационного центра
канала, происходит переход на более
высокий энергетический уровень, т.е.
изменение электронной
степени свободы
и переход на ветвь КП
.
2 – происходит изменение конформационной координаты Q, релаксация в минимум левой ветви КП (С). В итоге RyR-канал электронно и конформационно закрыт.
Переходы между состояниями на вышеприведенных схемах в рамках рассматриваемой модели будут подробно описаны далее.
2.1.5 Динамика конформационной координаты
Изменение конформационной координаты канала Q можно описать с помощью уравнения Ланжевена [105]:
(2.7)
где M – параметр эффективной массы RyR-канала (для упрощения принимается М=1);
–переменная,
описывающая электронное состояние
канала (
– электронно-закрытое, и
– электронно-открытое);
– параметр конформационного «трения»
(параметр диссипации);
– сила случайных температурных флуктуаций
в конформационной динамике RyR-канала.
Если пренебречь
случайными тепловыми флуктуациями (
),
то уравнение Ланжевена является вторым
законом Ньютона. В этом случае (2.5)
описывает процессы релаксации
конформационной координаты по ветвям
конформационного потенциала к их
минимумам:
– для электронно-закрытого,
– для электронно-открытого канала
(рис.2.2).
2.1.6 Динамика электронной степени свободы
Как уже отмечалось ранее, процесс активации RyR-каналов напрямую зависит от концентрации Са2+ в cis-части и является результатом связывания ионов Са2+ с активационными центрами RyR-канала. Результаты экспериментальных исследований [106, 107] говорят о том, что переходы в открытое состояние наблюдаются при достаточно малом уровне Са2+ в cis-части (~0.1 мкM).
Электронная степень свободы, в отличие от конформационной, является дискретной и в ЭК-модели описывается как марковский процесс с двумя состояниями (0 и 1).
Описание электронных переходов в ЭК-модели строится на предположении, что они обусловлены взаимодействием Са2+ с активными центрами канала, а их интенсивность зависит от энергии ионов Са2+ и энергии канала в текущем конформационном состоянии.
По аналогии с
теорией неупругого рассеяния
предполагается, что взаимодействие
Са2+
с активными центрами рассматривается
как резонансное рассеяние. В связи с
этим, такие переходы носят резонансный
характер, причем вероятность перехода
может быть аппроксимирована спектральной
функцией Лоренца (формулой Брейта-Вигнера)
[109]:
(2.8)
Где Е
энергия ионов,
энергия возбуждения (разность энергий
двух ветвей конформационного потенциала),
полуширина резонансного пика,
амплитуда вероятности, зависящая от
концентрации ионов
в диадном пространстве.
Таким образом, описанные выше электронные переходы являются следствием взаимодействия активных центров канала с ионами Са2+, и их интенсивность зависит не только от концентрации данного субстрата, но и от энергии ионов Са2+, способных связываться с активными центрами.
