Задачи молекулярной механики мышечного сокращения тема автореферата и диссертации по механике, 01.02.08 ВАК РФ

Метальникова, Надежда Алексеевна АВТОР
кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Москва МЕСТО ЗАЩИТЫ
2013 ГОД ЗАЩИТЫ
   
01.02.08 КОД ВАК РФ
Диссертация по механике на тему «Задачи молекулярной механики мышечного сокращения»
 
Автореферат диссертации на тему "Задачи молекулярной механики мышечного сокращения"

МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ имени М.В. Ломоносова

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

МЕТАЛЬНИКОВА Надежда Алексеевна ЗАДАЧИ МОЛЕКУЛЯРНОЙ МЕХАНИКИ МЫШЕЧНОГО СОКРАЩЕНИЯ

Специальность 01.02.08- биомеханика

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

рэ—-

2 6 СЕН 2013

Ои:

Москва-2013

005533821

Работа выполнена на кафедре гидромеханики механико-математического факультета и в лаборатории биомеханики Института механики Московского государственного университета имени М.В. Ломоносова.

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

Цатурян А.К., доктор физико-математических наук, старший научный сотрудник, ведущий научный сотрудник НИИ механики МГУ имени М.В. Ломоносова

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

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

Штейн A.A., кандидат физико-математических

наук, ведущий научный сотрудник

НИИ механики МГУ имени М.В. Ломоносова

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

ФГБУН Институт иммунологии и физиологии Уральского отделения РАН (ИИФ УрО РАН)

Защита состоится 18 октября 2013 г. в 15 часов на заседании диссертационного совета Д.501.001.89 при Московском государственном университете имени М.В.Ломоносова по адресу: 119991, г. Москва, Ленинские горы, Главное здание МГУ, аудитория 16-10.

С диссертацией можно ознакомиться в Научной библиотеке механико-математического факультета МГУ имени М.В. Ломоносова.

Автореферат разослан "■ff)" сентября 2013 г.

Учёный секретарь

диссертационного совета Д.501.001.89, доктор физико-математических наук

Я

В.В. Измоденов

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы. Исследование механизма мышечного сокращения является одной из актуальных и важных проблем современной биомеханики. Мышцы представляют собой специализированные органы, способные генерировать активные механические усилия и совершать значительную работу. Большую часть объема мышечной клетки занимают миофибриллы -цилиндрические структуры, образованные одинаковыми повторяющимися элементами - саркомерами. Саркомер ограничен с двух сторон плотными Z-дисками. К этим дискам с обеих сторон крепятся тонкие актиновые нити. В середине саркомера располагаются толстые нити, построенные преимущественно из другого сократительного белка, миозина. Укорочение происходит в результате перемещения (скольжения) нитей друг относительно друга (H.E.Huxley, J.Hanson, 1954; A.F.Huxley, R. Niedergerke, 1954; A.F.Huxley, 1957; H.E.Huxley 1969; A.F.Huxley, R.M.Simmons, 1971). Вследствие такого скольжения уменьшается длина каждого саркомера и всей мышцы в целом. Перемещение нитей друг относительно друга происходит за счет циклического замыкания и размыкания контактов между нитями актина и миозина. Эти контакты формируются головками миозина, которые расщепляют АТФ и преобразуют освободившуюся химическую энергию в механическую, генерируя тянущее усилие и совершая работу. Миозиновые головки тянут нити актина к центру саркомера, что приводит к сокращению мышцы.

Более поздние эксперименты показали, что нити актина и миозина растяжимы (H.E.Huxley и др., 1994; К. Wakabayashi и др., 1994). Более того, оказалось, что, кроме растяжения актиновых нитей, происходит и изменение угла закрутки актиновой спирали (J. Bordas и др., 1999; А.К. Цатурян и др., 2005). Такие изменения угла составляют доли градуса на каждый мономер актина. Однако, поскольку их число в нити велико - 360-370 штук, то даже малые изменения угла спирали могут привести к большим поворотам конца нити. Было обнаружено, что прочное связывание миозиновых головок с актином изменяет спиральный угол актина и растягивает актиновые нити (А.К. Цатурян и др., 2005). Активация тонких нитей ионами Са2+ и пассивное растяжение расслабленных волокон в отсутствие Са2+ также приводят к изменению угла закрутки актиновой спирали (J. Bordas и др., 1999; Y. Takezawa и др., 1998). Хотя некоторые из имеющихся математических моделей сокращения учитывают растяжимость актиновых и миозиновых нитей (Д.А. Шестаков, А.К. Цатурян, 1998), ни одна из них не принимает во внимание закручивание актина. Выяснение вопроса о том, какое влияние это может оказать на работу мышцы, остается открытым.

Регуляция сокращения поперечно-полосатых мышц осуществляется путем изменения концентрации ионов Са2+ в мышечной клетке. Взаимодействуя с регуляторным белком тропонином (Тп) и связанным с ним тропомиозином (Тт), которые расположены на поверхности тонких актиновых нитей, Са2+ вызывает сокращение мышечной клетки. Удаление Са2+ приводит к расслаблению. Имеются экспериментальные данные, показывающие, что

эффективность кальциевой регуляции мышечного сокращения зависит от механических свойств Тт. В 2001 году D.A. Smith предложил рассматривать Тт как упругий стержень, который находится в поле электростатических сил и подвергается изгибу при присоединении тропонина или миозиновой головки к актину. Однако эта и последующие модели не учитывали спиральную форму тропомиозина и не позволяли объяснить характер зависимости активации мышечного сокращения от концентрации Са2+. Построение модели активации, учитывающей спиральную форму тропомиозина в недеформированном состоянии и влияние миозиновых головок на активацию сокращения, также актуально. Особенно важно это для сердечной мышцы, где Тп связывает только один ион кальция, но механическое напряжение резко возрастает в узком диапазоне изменения концентрации Са2+. Такую зависимость характеристик системы от управляющего параметра, в данном случае от концентрации Са2+, называют "кооперативной". Она обусловлена взаимодействием между соседними кальций-связывающими молекулами Тп, соединенными упругими молекулами Тт, и играет важную функциональную роль ("закон сердца" Франка-Старлинга).

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

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

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

• Построение механической модели активации сокращения мышц ионами кальция с учетом спиральной формы тропомиозина в недеформированном состоянии.

• Исследование влияния параметров модели на характер активации мышечного сокращения ионами кальция.

Научная новизна. В работе получены следующие новые результаты, выносимые на защиту:

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

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

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

• Предложена механическая модель, которая рассматривает тропомиозиновую цепь как упругий стержень, имеющий в недеформированном состоянии форму спирали и способный скользить без удлинения по поверхности нити актина. Модель описывает существенные черты Са2+ активации тонких нитей в сердечной мышце и объясняет высокую кооперативность связывания ТпС ионами Са2+ даже в отсутствие миозиновых головок.

• Модель также воспроизводит увеличение Са-чувствительности тонких нитей при присоединении миозиновых головок без существенного увеличения кооперативности. Кроме того, модель объясняет результаты экспериментов, в которых некоторые молекулы ТпС были заменены мутантными, нечувствительными к Са2+.

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

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

Апробация работы. Результаты, полученные в работе, были доложены на следующих научных конференциях: Конференция-конкурс молодых ученых НИИ механики МГУ (2009, 2011, 2012); 10-ая Всероссийская конференция по биомеханике "Биомеханика 2010" (Саратов, 2010), Международная конференция "Биологическая подвижность: от фундаментальных исследований до нанотехнологий" (Пущино, 2010), Конференция МГУ "Ломоносовские чтения" (2012), Международная конференция "Биологическая подвижность:

фундаментальные и прикладные исследования" (Пущино, 2012), 41-ая Европейская мышечная конференция (Родос, Греция, 2012).

Результаты работы были также доложены и обсуждены на специализированных научных семинарах: семинар по биофизике под руководством проф. Роджера Кука (Roger Cooke, UCSF, Сан-Франциско, Калифорния, США, 2012), семинар по механике сплошных сред под руководством акад. РАН А.Г. Куликовского, проф. В.П. Карликова, члена-корр. РАН О.Э. Мельника (НИИ механики МГУ, Москва, 2012,2013).

Публикации по теме диссертации. Основные результаты работы изложены в 10 научных публикациях, из которых 2 ([4], [10]) - статьи в журналах, входящих в Перечень ВАК на момент публикации, и 8 — тезисы докладов. В работах [4], [6], [8], [10], опубликованых в соавторстве с научным руководителем, автору принадлежит участие в постановке задачи, разработка алгоритмов решения, численное моделирование и анализ результатов. Все положения, выносимые на защиту, получены лично соискателем.

Структура работы. Диссертация содержит введение, четыре главы, заключение и список литературы. В работе имеется 49 рисунков, 2 таблицы и 133 библиографические ссылки. Общий объем диссертации составляет 115 страниц.

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

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

До недавнего времени полагали, что при таком механизме генерации движения, получившем название "скользящих нитей", длина самих нитей не меняется. Позднее выяснилось, что и актиновые, и миозиновые нити растяжимы, а актиновые нити, кроме того, могут закручиваться под влиянием различных факторов, таких как растягивающие напряжения, крутящие моменты, образование белковых комплексов. Хотя угол закрутки составляет доли градуса на мономер актина, закрутка всей нити, содержащей >300 мономеров, может быть существенной. Выделены факторы, влияющие на изменение закрутки нитей актина, и оценены величины некоторых параметров, влияющих на этот процесс.

Z-my\\л9\

иебулин

/ \

в-актины

Я-акт^

а) г-линия ¡^тонкие нити

Ч '—

саркомер

тропонин тропомиозин

б)тонкая нить

М-линия

миофибрилла

Б1, глобулярная головка

в) толстые нити

д) молекула миозина Э1-82 соединение

Рис. 1. Схема строения саркомера. Показана упаковка основных мышечных белков в тонкие (а, б) и толстые (в, д) нити. Число миозиновых молекул в ярусе толстой нити может меняться в зависимости от вида животного и типа мышц, в мышцах высших позвоночных оно равно трём.

Регуляция сокращения поперечно-полосатых мышц осуществляется путем изменения концентрации ионов Са2+ в цитоплазме. Взаимодействуя с регуляторными белками, тропонином (Тп) и тропомиозином (Тт), входящими в тонкие нити, Са2+ вызывает сокращение мышечной клетки. Удаление Са2+ приводит к расслаблению. Молекула Тш образована двумя перекрученными а-спиралями и имеет форму длинной спирали. Две таких молекулы располагаются напротив друг друга в желобах длинной актиновой спирали (рис. 2). Соседние молекулы Тт соединяются друг с другом "голова к хвосту" и образуют непрерывный изогнутый стержень длиной около 1 мкм, который тянется вдоль всей актиновой нити.

Актин Тропомиозин (Тт)

Рис. 2. Белки тонкой нити: актин, тропомиозин (Тт), тропонин (Тп).

С каждой молекулой Тт связана молекула Тп, состоящая из трех субъединиц - ТпС, Тп1 и ТпТ. В отсутствии Са2+ подвижный участок Тп1 присоединяется к актину и сдвигает весь Тп-Тгп комплекс так, что он блокирует участки связывания миозиновых головок на актине. При связывании Са2+ с ТпС мобильный участок Тп1 отсоединяется от актина и присоединяется к ТпС, а весь Тгп-Тп комплекс поворачивается относительно оси актиновой нити и освобождает её участок для присоединения миозиновых головок, т.е. переводит систему из блокированного состояния (В) в закрытое (С). Наконец, присоединение миозиновой головки к актину вызывает дальнейший сдвиг Тт спирали и переводит систему в открытое состояние (О).

Зависимость вероятности образования комплексов СаТпС, V, от концентрации Са2+, С, хорошо описывается формулой Хилла:

V = —— СП

М?о+[с]п '

с показателем кооперативности п (его также называют коэффициентом Хилла). Здесь С50 - концентрация, при которой половина молекул ТпС связана с Са2+. В ТпС сердечной мышцы имеется лишь один участок связывания с ионом Са2+. Несмотря на это, и а 2.5. Причины такой кооперативности (независимое связывание одного иона Са2+ с каждой молекулой ТпС соответствует п= 1), а также зависимости параметров, входящих в формулу (1), от механических величин представляют предмет интенсивной дискуссии.

В обзоре также рассмотрены современные данные о механических свойствах всех белковых компонент саркомера: актина, миозина, тропомиозина и др. Также описаны современные представления о молекулярных механизмах сокращения мышц и регуляции этого процесса. Особенное внимание уделено регуляции сокращения сердечной мышцы и процессам, лежащим в основе так называемого "закона сердца Франка-Старлинга". В последнем разделе главы 1 описаны существующие математические модели сокращения мышц и его регуляции ионами Са2+.

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

Глава 2 посвящена численному исследованию влияния закручивания и растяжения белковых нитей саркомера на характеристики сокращения мышц. Выведены основные уравнения, описывающие изменение угла закрутки актиновой нити ф^Л), перемещение нитей актина иа(ъ, I) и миозина удлинения саркомера и концентрации силогенерирующих актин-

миозиновых мостиков п(б, $ в мышечном волокне в зависимости от времени t и осевой координаты 5 (рис. 3). В разделе 2.1 дана постановка задачи о релаксации напряжения и других механических характеристик в ответ на

ступенчатую деформацию (рис. ЗА), а в разделах 2.2, 2.3 - приведены результаты ее численного решения.

Рис. 3. А: изменение удлинения (у) и ответ развиваемого напряжения (Р) на ступенчатую деформации мышечного волокна. Б: схема строения полусаркомера мышцы и обозначения используемых в модели величин.

Система уравнений, описывающих изменения угла ф, относительного перемещения нитей £/=иа-мт и концентрации мостиков п, имеет следующий вид:

ôn(s, t)

а2иМ~Рщи {..«Wo^feO-'*

ds2

с„ + с„

3n(s, t) ds

dt

= fc+( 1 - n(s, t)) - fe_n(s, t)

pn0U0(_s) (2)

d2<j>{s,t)

ds2

-—<P(s. t)

(-ZTX--^)

3n(s, t) as

d2U(s, t) as2

— ФоОО--Z-.

Здесь ca, cm - коэффициенты продольной податливости (растяжимости) актиновых и миозиновых нитей, кинетические константы /с+(£/ — С/0) =

-k(U-U0)h Zk(U-U0)h

k_(U — U0)rje квт и /c_(f/ — Uq) — е квт , описывают структурные изменения в актин-миозиновых мостиках в зависимости от относительного перемещения нитей (A.F.Huxley, R.M.Simmons, 1971); к- жесткость одного миозинового мостика; h — "шаг", т.е. осевое перемещение при переходе миозинового мостика из состояния, слабо присоединенного к актину, в прочно связанное силогенерирующее состояние; Я - константа, описывающая зависимость скорости отсоединения от относительного перемещения; rj -константа скорости перехода из "слабого", "нетянущего" в прочное, тянущее состояние в изометрическом сокращении; кв, Т - постоянная Больцмана и абсолютная температура; р - погонная жесткость присоединенных к актину миозиновых мостиков на единицу длины при максимально возможном их числе; щ - погонная концентрация присоединенных мостиков. Параметры у, L,

ô, ß, x> характеризуют раскрутку нити актина под действием растягивающей силы, жесткость актиновой нити на кручение, ее растяжимость и закрутку под действием присоединенных миозиновых мостиков, жесткость мостиков на кручение, а также крутящий момент, создаваемый этими мостиками в изометрических условиях, соответственно. Выражения для полного механического напряжения, развиваемого мышечным волокном в расчете на одну актиновую нить, F, и удлинения волокна, у, имеют вид:

Fit) = J f(t)ds,f = phn(s, t)n0 + pn0(i/m(s, t) - t/° (s) - Ua(s, t) + £/°(s)),

y = Fcas1 + Fcm(/ - s2) + ca J/2 aa(s, t)ds + Sn0 fSz n(s, t)ds + u\s . (3)

Граничные условия для угла: ф(0) = 0, M(s2) = 0 , 0|Sl_o = 0|Sl+o, ^Ist-o = M|s1+0> где М- момент внутренних сил, = —fl, а £2 — момент сил, создаваемый миозиновыми головками. Граничные условия для перемещения имеют вид = -caF\Sl - ôn0n\Si , = cmF\Sl - dn0n\S2,U(0) = 0 и

содержат функционал от искомого решения, F. В граничное условия также входят распределения перемещений нитей иао, ито и развиваемое напряжение F0 в ходе стационарного сокращения при постоянной длине, которые находятся из аналитического решения соответствующей задачи. Распределения искомых величин в ходе мгновенно-упругой реакции в предположении о том, что длительность деформирования много меньше характерного времени структурных изменений в актин-миозиновых мостиках fcj1 были найдены аналитически. В этом случае система уравнений (2) сводится к системе обыкновенных дифференциальных уравнений. Распределения, полученные в процессе мгновенно-упругой реакции, были использованы как начальные условия для моделирования процесса релаксации. После быстрой ступенчатой деформации сразу после по окончания мгновенно упругой реакции, напряжение в мышечном волокне релаксирует на миллисекундной временной шкале (A.F. Huxley, R.M. Simmons, 1971). При этом изменения перемещения и угла зависят не только от пространственной координаты, но и от времени. Кроме того, в зависимости от времени изменяется распределение величины п — концентрации "тянущих" миозиновых мостиков, от которой зависят и растяжение, и кручение нитей. Задача была решена численно методом параметрической прогонки с итерациями. Из сравнения с опубликованными экспериментальными данными полученных решений задач об установившемся сокращении и мгновенно-упругой реакции была проведена оценка параметров, которые не были измерены экспериментально: у> 0 - коэффициент, описывающий раскрутку нити актина под действием растягивающей силы; Q0 -"погонного" крутящего момента, создаваемого миозиновыми головками в изометрическом сокращении. Затем, используя полученные оценки этих параметров, были проведены численные расчеты для процесса релаксации мышечного волокна. Примеры моделирования ответов на растяжение мышцы

из состояния изометрического сокращения с _у0=1.62нм до у = 4 нм и её укорочения до у = -3 нм показаны на рис. 4, где представлен временной ход изменения напряжения, развиваемого волокном, и среднего значения концентрации тянущих мостиков п.

Рис. 4. Изменение силы Р, развиваемой волокном в пересчете на одну актиновую нить, (А) и средней концентрации мостиков п (Б) в ответ на ступенчатую деформацию. Пунктирные кривые соответствуют изменению силы при удлинении волокна до у = 4 нм; сплошные - укорочению до _у = -3нм. В: изменение угла поворота конца актиновой нити

Из графиков видно, что модель описывает экспериментально обнаруженную нелинейность релаксации напряжения в мышце - релаксация ускоряется при укорочении и замедляется при растяжении (раздел 2.4). Также оказалось, что при релаксации напряжения после растяжении или укорочения мышечного волокна угол закрутки актиновой спирали не возвращается к своему исходному значению, а продолжает изменяться в том же направлении, что и при мгновенно упругой реакции. Изменение угла поворота конца актиновой нити превышает 180° (рис. 4В).

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

Перемещение точки тропомиозиновой спирали имеет осевую, и(г), и угловую, (р(г), компоненты, зависящие от координаты г, измеряемой вдоль оси тонкой нити (рис. 5). Предполагая, что тропомиозин нерастяжим, а его перемещения малы, получено соотношение, связывающее между собой радиальную и осевую компоненты вектора перемещения и'(^) + а21р(р'{г) = 0. Чтобы упростить дальнейшие выкладки, вводится параметр д0 = 1 + а2-ф2, где

а - радиус цилиндрической поверхности по которой перемещается осевая линия тропомиозинового стержня, а гр - параметр спирали, определяемый так, что её полному обороту соответствует осевое смещение 2л:/1р (рис. 5).

Рис. 5. Модель Тт спирали, показаны базисные вектора лабораторной (ё;) и сопутствующей систем координат. А: актиновый цилиндр с Тт спиралью на его поверхности. Б: развертка поверхности цилиндра.

В модели предполагалось, что изгибная жесткость Тш не зависит от направления изгиба, т.е. модули изгибной жесткости одинаковы К1=К2=К. Для коротковолновых изгибов такое предположение неверно, однако в данной задаче оно, по-видимому, оправдано, так как полупериод закрутки двух альфа-спиралей в суперспираль Тт составляет 5.5 нм, что много меньше характерного радиуса изгиба всего тропомиозинового стержня (=20 нм). Поэтому можно считать, что энергия изгиба-кручения Тт принимает вид: е = 0.5(^(Дш2 + Ао>|) + СДш|), где Лсои Лсо2, Ла>3 - изменения компонент вектора кривизны-кручения Тт в сопутствующей системе координат. Изгибная жесткость Тт К много больше его жесткости при кручении б в силу особенностей молекулярной структуры Тт. В положении равновесия (С-состояние) Тт удерживается электростатической силой -аср, которая определяется ионными и электростатическими взаимодействиями, и пропорциональна смещению. Экспериментально было показано, что при присоединении миозиновой головки к актину Тт, смещаясь, скользит по поверхности актина, не прокатываясь. Таким образом, упругая энергия Е деформированного Тт определяется энергией изгиба и работой сил актин-Тт взаимодействия и для участка стержня между точками Ь и с равна:

(«Ш)2 + К + (2агРср'Ю)2^ (4)

Была введена безразмерная осевая координата: х = , где Е, = (4Ка2/а^дй)1^ ~ характерная длина, и интеграл энергии (4) преобразован в безразмерный вид, а из вариационного принципа получено уравнение равновесия Тт стержня как условие минимума энергии:

<р'у(х) - 4Р(р"(х) + 4<р{х) = 0, (5)

где Р = 2аф2^Кд0/а - безразмерное натяжение нити.

Присоединение Tnl или миозиновой головки к актину фиксирует точки тропомиозинового стержня на актине, задавая определенные значения угла ср, равные <рх или (рм = —Хф-т, соответственно, где х > 0, так как тропонин и миозин смещают тропомиозин в противоположных направлениях. Предполагается, что фиксация Тш не создает изгибающий момент, поэтому в точках фиксации и само угловое перемещение, и его первые две производные непрерывны. Тропониновые комплексы расположены на расстоянии 38.5 нм друг от друга вдоль оси тонкой нити и, соответственно, могут присоединять тропомиозиновый стержень к каждому седьмому мономеру актина, смещая тропомиозин в противоположную сторону. Миозиновая головка может присоединиться к любому актиновому мономеру, кроме тех, к которым могут присоединиться подвижные участки Tnl. На тонкой нити длиной 1.05 мкм имеются 27 Тп комплексов, в промежутках между которыми находятся 156 актиновых мономеров, доступных для миозиновых головок, итого 183=27+6x26 потенциальных участков закрепления стержня. Присоединение ТпС или миозиновой головки к актину или отсоединение от него представляют собой стохастический процесс, в котором вероятности присоединения и отсоединения зависят от концентрации ионов Са2+, а также от изменения упругой энергии Тш стержня, которое сопровождает такое присоединение или отсоединение. Решение уравнения (5) для произвольного числа точек закрепления требует достаточно большого объема вычислений. Поэтому мы использовали трехточечное приближение (раздел 3.2), считая, что изменение энергии деформированной спирали, связанное с ее присоединением или отсоединением от актина, такое же, как для бесконечно длинной спирали, закрепленной в точках, совпадающих с левой и правой ближайшими точками закрепления. Были получены решения уравнения (5) для бесконечно длинной тропомиозиновой спирали, закрепленной в одной, двух, трех или четырёх точках. Предварительные расчёты показали, что трехточечная аппроксимация приближает решение намного лучше, чем двухточечная, использовавшая в работах (D.A. Smith и др., 2001, 2003, 2011), а учёт еще одной, четвертой внешней точки крепления, наоборот, практически не влияет на результат, существенно увеличивая время вычислений. Энергия бесконечно длинного стержня, закрепленного Tnl в единственной точке при ср= <рг, равна Ет = 9о(1 + Ю- Поскольку присоединение и отсоединение Тш от актина происходят в результате химических реакций, в основе которых лежит броуновское движение молекул, мы ввели еще один безразмерный параметр у = a^cpr-Jg0(1 + Р)/квТ - отношение упругой энергии к тепловой. Раздел 3.3 посвящен определению условий на концах стержня, то есть определению состояний 1-го и 183-го участка закрепления Тт, а также введению корректировочной функции, которая позволила более точно вычислять изменения энергии Тш стерня конечной длины при отсоединении и присоединении крайних элементов.

В разделе 3.4 в соответствии с описанным ранее механизмом взаимодействия субъединиц Тп, Са2+ и актина рассмотрена следующая кинетическая схема (рис. 6).

Рис. 6. Кинетическая схема взаимодействия субъединиц тропонина, Са2+ и актина, Кса, К\, КА - константы равновесия соответствующих реакций.

Здесь Кса, К\ и КА - константы равновесия присоединения Са2+ к ТпС, переключающего сегмента Тп! к ТпС и присоединения ингибиторного домена Тп1 к актину, соответственно, а 0 < е < 1 параметр, характеризующий связывание Тп1 с ТпС в отсутствии Са2+. Схема предполагает, что ТпС может связать только один ион Са2+, как в сердечной мышце, и такое связывание способствует отсоединению Тп1 от актина и присоединению к ТпС. Обозначим С = Кса[Са2+], у, 2 — безразмерную концентрацию Са2+ и вероятности пребывания ТпС и ингибирующего участка Тп1 в свободном состоянии, соответственно. Тогда для стационарного процесса вероятности пребывания компонент тропонинового комплекса в состояниях СаТпС, СаТпСТпГ, ТпСТп1 равны Су, СугК1 и уегК[, соответственно, при условии, что ингибирующий участок Тп1 не присоединен к актину. Вероятность присоединения Тп1 к актину равна при этом гКА, а у и г удовлетворяют условиям:

у + уС + уСгК1 + уггК1 = 1,

г + гКА+ уегК, + уСгК1 = 1. (6)

Система (6) представляет собой уравнения химической кинетики (т.н. закон действующих масс), которые записаны не для концентраций веществ, а для вероятностей пребывания молекул в том или ином состоянии.

В разделе 3.5 рассмотрен метод Монте-Карло для численного решения задачи определения среднестатистических характеристик системы в зависимости от концентрации ионов кальция. Активация сокращения характеризуется вероятностями пребывания ТпС и Тп1 в свободном и связанном состоянии, а также вероятностью связывания миозиновых головок с актином. Изменения состояния системы имеют стохастический характер и при заданной величине безразмерной концентрации кальция С, определяются ее текущим состоянием, т.е. представляют собой стационарный марковский процесс. Поскольку даже в одном саркомере имеются тысячи тонких нитей, все экспериментально наблюдаемые характеристики являются средними по

ансамблю. Скорости присоединения и отсоединения от актина Тп1 или миозиновых головок определяются изменением упругой энергии Тш стержня вследствие такого перехода. Было сделано предположение, что связывание ингибирующего участка Тп1 с актином и присоединение миозиновой головки -обратимые стохастические процессы, такие, что отношение скоростей прямого и обратного переходов для Тп1 и миозина определяются соотношениями:

Р] = кт+/кт_ = гКАехр (-у(Я3 - Я2)) , р, = км+/км_ = рехр (-у(Е3 - Е2)) ,

где индексы "+" и "-" соответствуют присоединению и отсоединению, аГиМ-Тп1 и миозину (рис. 7), соответственно, р - доступность миозиновых головок для связывания с актином, Е3, Е2 - безразмерные (нормированные на Ет) упругие энергии тропомиозинового стержня предварительно рассчитанные по формуле (4) для всех возможных конфигураций, т.е. расстояний между двумя или тремя точками фиксации Тт стержня на актине молекулами Тп1 или миозиновыми головками.

Рис. 7. Схематическое изображение фрагмента Тт стержня (черные линии). Тропониновые комплексы показаны кругами (черные соответствуют присоединенным Тп1, серые неприсоединенные Тп1). Аналогично, присоединенные миозиновые головки показаны черными эллипсами, неприсоединенные - серыми. Левый пунктирный прямоугольник показывает сегмент тонкой нити до и после присоединения миозиновой головки к актину (константа скорости &м+)- Правый пунктирный прямоугольник показывает другой сегмент до и после отсоединения молекулы Тп1 от актина (константа скорости кт„1_).

Для определения средних стационарных характеристик системы, т.е. средней доли актин-связанных молекул Тп1, доли комплексов СаТпС и вероятности присоединения к актину миозиновых головок в зависимости от концентрации кальция С использован метод Монте-Карло. Основная проблема заключалась в большой размерности фазового пространства 2183 (183 элемента, каждый из которых может находиться в любом из двух состояний -присоединенном или отсоединенном). Для определения характеристик системы был применен видоизмененный метод Метрополиса-Гастингса, который позволил добиться эффективной сходимости средних характеристик при относительно небольшом числе псевдослучайных испытаний. Основная идея

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

Каждому элементу была поставлена в соответствие вероятность присоединения Ifpi, если он был отсоединен, и вероятность отсоединения (Pi)-1/3, если он был присоединен к актину. С учётом этих вероятностей производился псевдослучайный выбор элемента. После того, как элемент был выбран, решался вопрос о его переходе в новое состояние: вероятность присоединения вычислялась по формуле min[l, (р;)1^3], а отсоединения-min[l, (pi)-1/3] . в результате такого правила выбора вероятность присоединения элемента равнялась ^/p]min[l; lfpi\ , а отсоединения -l/^min[l; l/\fpi], соответственно. Таким образом, отношение вероятностей прямого и обратного перехода сохранялось равным ph как и должно быть.

В результате удовлетворительная сходимость средних характеристик системы была достигнута после 106-107 испытаний, что на много порядков меньше размерности системы. В разделе 3.6 по данным опубликованных экспериментов произведена оценка параметров, необходимых для моделирования кальциевой регуляции.

В Главе 4 представлены результаты моделирования кальциевой активации тонких нитей в отсутствии (раздел 4.1) и при наличии миозиновых головок (раздел 4.2), исследовано влияние параметров задачи на кальциевые кривые и на кажущийся коэффициент Хилла - раздел 4.3. Базовый набор параметров был следующим: К] = 1000, Кл = 200, у = 4, Я = 2 (§ = 19.25 нм), ß= 3.5, ^=0.4, е - 0. В отсутствии миозиновых головок при значениях / меньше 3.5, получить достаточно высокие значения коэффициента Хилла не удалось. Это означает, что высокая кооперативность активации тонких нитей требует, чтобы механическая энергия, связанная с присоединением Tnl к актину, была существенно больше тепловой. При больших значениях параметра у кальциевые кривые могут быть круче, чем представленные в работе. Расчётные кривые при А=2 (§= 19.25 нм), /3=3.5 показывают достаточно резкие и высококоперативные зависимости активациии от концентрации кальция (и » 2.5) даже при отсутствии миозиновых головок, что согласуется с данными (Y.B Sun и М. Irving, 2009, 2010), которые показали, что тонкие нити активируются кооперативно, даже если присоединение миозиновых головок заблокировано. Увеличение параметра ß приводит к снижению Са-чувствительности и к увеличению кооперативности: кажущийся коэффициент Хилла, п, увеличивается (рис. 8А). При постоянной характерной длине Тш стержня цепи (постоянная изгибная жесткость и а = const) ß зависит от параметра спирали актина у/. В модели, описанной ранее D.A. Smith и коллегами (2001, 2003, 2011), была использована теория прямолинейного упруго стержня, т.е. ß=0. Увеличение параметра ß до значения, соответствующего реальной спиральной структуре актина, приводит к более высокой кооперативности Са-активации тонких нитей, соответствующей

экспериментальным данным (п » 2.5). Изменение X характеризует изменения в изгибной жесткости Тт, К, или параметра а, который удерживает Тгп в С-состоянии. Увеличение X приводит к сдвигу кривой активации влево, т.е. увеличивается Са2+-чувствительность и незначительно уменьшается кооперативность (рис. 8Б). Так как X является отношением осевого расстояние между соседними молекулами тропонина (38,5 нм) к характерной длине Тт цепи, оно увеличивается, если изгибная жесткость Тт уменьшается. С другой стороны, увеличение коэффициента а, который характеризует актин-миозиновые электростатические силы взаимодействия, удерживающие Тт в С-состоянии, приводит к увеличению Я.

отсутствие миозиновых головок для различных значений параметров /? и X. А: Расчет для X = 2, значения /?, приведены на графике. Б: расчеты для значений /? = 3.5, величины X показаны на графике. Средние значения доли несвязанных с актином молекул Tnl показаны пунктирными линиями, а связанных с ионами Са2+ молекул ТпС — сплошными.

При присоединении миозиновых головок к актину расчетные кривые кальциевой активации смещаются влево, показывая повышенную чувствительность без существенного увеличения наклона кривой (рис. 9). Причина этого состоит в том, что при прочном связывании миозиновой головки с актином Тт стержень сдвигается и облегчает отсоединение соседних молекул Tnl от актина. В работах Y.B. Sun и M.Irving (2009, 2010) была измерена относительная доля молекул ТпС, связанных с Са2+ в состоянии ригора, в котором все миозиновые головки прочно соединены с актином, а Тт -неподвижен, и обнаружено, что [Са2+]50 в этом случае гораздо ниже, чем в присутствии АТФ. Кооперативность при этом также отсутствовала, п= 1. Для моделирования состояния ригора в настоящей работе задавали Кл = 0, чтобы отключить возможность связывания Tnl с актином. Описанные выше свойства кальциевой активации в состоянии ригора хорошо описываются моделью (рис. 9А).

о -■-■---.

0.0001 0.001 0.01 0.1 1 Нормализованная концентрация С=КГя[Са2Ч

Рис. 9. Моделирование влияния головок миозина на кривые Са-активации. А: средняя доля ТпС, связанных с ионами Са2+ при различных значениях доступности миозиновых головок, р. Расчет КА = 0 соответствует состоянию ригора, показан пунктирной линией. Б: средняя доля Тп1, не связанных с актином. /? = 3.5, X = 2 для всех приведенных расчетов.

Графики изменения вероятности пребывания миозиновых головок в присоединенном состоянии в зависимости от концентрации кальция и при различных значениях доступности миозиновых головок представлены на (рис. 10). Эти результаты хорошо согласуются с результатами экспериментов, в которых связывание миозиновых головок с актином было блокировано блебистатином или при помощи других биохимических воздействий (О.Р. Рагтап и др. 2010).

Рис. 10. Нормализованная средняя занятость мономеров актина миозиновыми головками в зависимости от безразмерной концентрации Са2+ при различных значениях доступности миозиновых головок, р. Вставка: ненормированная занятость актина миозиновыми головками при тех же условиях (/? = 3,5, Я = 2, значения р показаны на графике).

Результаты расчётов показывают, что предложенная нами модель впервые позволила объяснить недавно обнаруженную экспериментально

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

• Построена кинетическая модель сокращения саркомеров поперечнополосатых мышц, учитывающая кручение актиновых нитей под действием внешних сил и моментов, а также структурных перестроек. С использованием экспериментальных данных проведена оценка параметров модели, описывающих кручение актиновой нити под действием растягивающей силы и момента, создаваемого присоединенными миозиновыми головками.

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

• Построена математическая модель регуляции актин-миозинового взаимодействия, в которой тропомиозин рассматривается как упругий стержень, имеющий в недеформированном состоянии спиральную форму. Модель учитывает взаимодействие тропомиозина с актином и тропонином, а также связывание последнего с ионами Са2+ и актином.

• Разработана модификация алгоритма Метрополиса-Гастингса, позволяющая эффективно вычислять среднестатистические стационарные характеристики стохастической марковской системы высокой размерности.

• Проведены численные расчеты зависимостей величин, характеризующих активацию мышечного сокращения, от концентрации ионов Са2+. Исследовано влияние параметров модели на эти зависимости. Показано, что предложенная модель позволяет объяснить высокую кооперативность системы активации сокращения мышцы в отсутствии и в присутствии миозиновых головок.

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Метальникова H.A. Математическая модель механики саркомера мышцы с учетом растяжения и кручения актиновых нитей. // Труды конф.-конкурса молодых ученых НИИ механики МГУ 14-16 октября 2009 г, 2010 М., МГУ. С. 225-232.

2. Метальникова H.A. Учет растяжения и кручения актиновых нитей в модели механики саркомера поперечно-полосатых мышц. // Сб. "Материалы X Всероссийской конференции: Биомеханика 2010". Саратов. 2010. С. 116.

3. Metalnikova N.A. Twist and stretch of actin filaments in sarcomere: a model study. // Biological motility: from fundamental achievements to nanotechnologies. Pushchino. 2010. P. 173-174.

4. Метальникова H.A., Цатурян A.K. Математическая модель механики саркомера поперечно-полосатых мышц с учетом растяжения и кручения актиновых нитей. // Биофизика. 2010. Т. 55(5). С. 892-898.

5. Метальникова Н.А. Математическая модель влияния механических факторов на кальциевую регуляцию мышечного сокращения. // Труды конференции-конкурса молодых ученых НИИ механики МГУ 12-14 октября 2011 г., 2012 М.,МГУ. С. 179-185.

6. Метальникова Н.А., Цатурян А.К. Механическая модель кальциевой активации тонких нитей в саркомерах мышц. // Ломоносовские чтения. Тезисы докладов научной конференции. Секция механики. 16-25 апреля 2012г., Москва. Изд. МГУ. 2012. С. 118.

7. Metalnikova N.A. A mechanistic model of Ca-regulation of thin filaments in striated muscle. // Biological motility: Fundamental and Applied Science. Pushchino. 2012. P. 131-135.

8. Metalnikova N.A., Tsaturyan A.K. A mechanistic model of Ca-activation of the thin filaments in cardiac muscle. // J. Muscle Res. and Cell Motil. 2012. v 33., P. 241.

9. Метальникова H.A. Влияние механических факторов на регуляцю мышечного сокращения. // Труды конф.-конкурса молодых ученых НИИ механики МГУ 8-10 октября 2012 г. (в печати).

10. Metalnikova N.A., Tsaturyan A.K. A Mechanistic Model of Ca-Regulation of Thin Filaments in Cardiac Muscle. // Biophysical Journal. 2013., v. 105(4), P. 941-950.

Заказ № 45-Р/09/2013 Подписано в печать 11.09.13 Тираж 70 экз. Усл. пл. 1,0

ООО "Цифровичок", тел. (495) 797-75-76 www.cfr.ru; e-mailAnfo@cfr.ru

 
Текст научной работы диссертации и автореферата по механике, кандидата физико-математических наук, Метальникова, Надежда Алексеевна, Москва

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ «МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

имени М.В. ЛОМОНОСОВА»

МЕХАНИКО-МАТЕМАТИЧЕСКИИ ФАКУЛЬТЕТ Кафедра гидромеханики

04201361492

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

МЕТАЛЬНИКОВА Надежда Алексеевна

ЗАДАЧИ МОЛЕКУЛЯРНОЙ МЕХАНИКИ МЫШЕЧНОГО

СОКРАЩЕНИЯ

Специальность 01.02.08 - биомеханика

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

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

Цатурян Андрей Кимович

Москва 2013

ОГЛАВЛЕНИЕ

ОГЛАВЛЕНИЕ.................................................................................................................................1

ВВЕДЕНИЕ. ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ......................................................3

ГЛАВА 1. ОБЗОР ЛИТЕРАТУРЫ...........................................................................................9

1.1. Строение поперечно-полосатых мышц.....................................................9

1.2. Современные представления о молекулярных механизмах сокращения поперечно-полосатых мышц и его регуляции.......15

1.3. Механические свойства белковых компонент саркомера..........24

1.3.1. Механические свойства нитей актина и миозина....................24

1.3.2. Механичнские свойства титина..........................................................28

1.3.3. Механические свойства тропомиозина..........................................28

1.3.4. Факторы, влияющие на расстояние между актиновыми и миозиновыми нитями.............................................................................32

1.4. Особенность кальциевой регуляции сердечной и скелетной мышцы.......................................................................................................................33

1.5. Математические модели мышечного сокращения и кальциевой активации сокращения поперечно-полосатых мышц...................41

1.5.1. Учет деформируемости белковых нитей в кинетических моделях мышечного сокращения......................................................41

1.5.2. Кинетические модели кальциевой активации сокращения поперечно-полосатых мышц................................................................42

1.5.3. Модели регуляции, рассматривающие тропомиозин как упругий стержень.......................................................................................47

ГЛАВА 2. МОДЕЛИРОВАНИЕ КРУЧЕНИЯ АКТИНА.................................................54

Глава 2. Таблица обозначений..........................................................................................54

2.1. Постановка задачи..............................................................................................56

2.2. Установившееся сокращение и мгновенно-упругая реакция...59

2.3. Процесс релаксации...........................................................................................63

2.3.1. Оценка параметров...................................................................................64

2.3.2. Моделирование процесса релаксации............................................67

2.4. Заключение.............................................................................................................73

ГЛАВА 3. МОДЕЛИРОВАНИЕ КАЛЬЦИЕВОЙ РЕГУЛЯЦИИ..................................74

Глава 3. Таблица обозначений..........................................................................................74

3.1. Модель механики тропомиозинового тяжа.........................................76

3.2. Вычисление энергий тропомиозиновой цепи....................................81

3.3. Граничные условия, корректировочная функция...........................83

3.4. Математическая модель Са-регуляции тропонинового комплекса................................................................................................................85

3.5. Метод Монте-Карло..........................................................................................86

3.6. Оценка параметров............................................................................................88

ГЛАВА 4. РЕЗУЛЬТАТЫ МОДЕЛИРОВАНИЯ АКТИВАЦИИ ТОНКИХ НИТЕЙ.............................................................................................................................................91

4.1. Кальциевая активация при отсутствии миозиновых головок.91

4.2. Влияние миозиновых головок.....................................................................93

4.3. Влияния параметров модели характер кальциевой активации................................................................................................................99

4.4. Заключение..........................................................................................................104

ЗАКЛЮЧЕНИЕ. ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ.....................................105

СПИСОК ЛИТЕРАТУРЫ........................................................................................................106

ВВЕДЕНИЕ. ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы. Исследование механизма мышечного сокращения является одной из актуальных и важных проблем современной биомеханики. Мышцы представляют собой специализированные органы, способные генерировать активные механические усилия и совершать значительную работу. Большую часть объема мышечной клетки занимают миофибриллы -цилиндрические структуры, образованные одинаковыми повторяющимися элементами - саркомерами. Саркомер ограничен с двух сторон плотными Z-дисками. К этим дискам с обеих сторон крепятся тонкие актиновые нити. В середине саркомера располагаются толстые нити, построенные преимущественно из другого сократительного белка, миозина. Укорочение происходит в результате перемещения (скольжения) нитей друг относительно друга (Н.Е. Huxley, J. Hanson, 1954; A.F. Huxley, R. Niedergerke, 1954; A.F. Huxley, 1957; H.E. Huxley 1969; A.F. Huxley, R.M. Simmons, 1971). Вследствие такого скольжения уменьшается длина каждого саркомера и всей мышцы в целом. Перемещение нитей друг относительно друга происходит за счет циклического замыкания и размыкания контактов между нитями актина и миозина. Эти контакты формируются головками миозина, которые расщепляют АТФ и преобразуют освободившуюся химическую энергию в механическую, генерируя тянущее усилие и совершая работу. Миозиновые головки тянут нити актина к центру саркомера, что приводит к сокращению мышцы.

Более поздние эксперименты показали, что нити актина и миозина растяжимы (Н.Е. Huxley и др., 1994; К. Wakabayashi и др., 1994). Более того, оказалось, что, кроме растяжения актиновых нитей, происходит и изменение угла закрутки актиновой спирали (J. Bordas и др., 1999; А.К. Цатурян и др., 2005). Такие изменения угла составляют доли градуса на каждый мономер актина. Однако, поскольку их число в нити велико - 360-370 штук, то даже малые изменения угла спирали могут привести к большим поворотам конца

нити. Было обнаружено, что прочное связывание миозиновых головок с актином изменяет спиральный угол актина и растягивает актиновые нити (А.К. Цатурян и др., 2005). Активация тонких нитей ионами Са и пассивное растяжение расслабленных волокон в отсутствие Са2+ также приводят к изменению угла закрутки актиновой спирали (J. Bordas и др., 1999; Y. Takezawa и др., 1998). Хотя некоторые из имеющихся математических моделей сокращения учитывают растяжимость актиновых и миозиновых нитей (Д.А. Шестаков, А.К. Цатурян, 1998), ни одна из них не принимает во внимание закручивание актина. Выяснение вопроса о том, какое влияние это может оказать на работу мышцы, остается открытым.

Регуляция сокращения поперечно-полосатых мышц осуществляется путем изменения концентрации ионов Са2+ в мышечной клетке. Взаимодействуя с регуляторным белком тропонином (Тп) и связанным с ним тропомиозином (Тш), которые расположены на поверхности тонких актиновых нитей, Са вызывает сокращение мышечной клетки. Удаление Са приводит к расслаблению. Имеются экспериментальные данные, показывающие, что эффективность кальциевой регуляции мышечного сокращения зависит от механических свойств Тт. В 2001 году D.A. Smith предложил рассматривать Тт как упругий стержень, который находится в поле электростатических сил и подвергается изгибу при присоединении тропонина или миозиновой головки к актину. Однако эта и последующие модели не учитывали спиральную форму тропомиозина и не позволяли объяснить характер зависимости активации мышечного сокращения от концентрации Са . Построение модели активации, учитывающей спиральную форму тропомиозина в недеформированном состоянии и влияние миозиновых головок на активацию сокращения, также актуально. Особенно важно это для сердечной мышцы, где Тп связывает только один ион, но механическое напряжение быстро возрастает в узком диапазоне изменения концентрации Са2+. Такую зависимость характеристик системы от

управляющего параметра, в данном случае от концентрации Са2+, называют "кооперативной". Она обусловлена взаимодействием между соседними кальций-связывающими молекулами Тп, соединенными упругими молекулами Тш, и играет важную функциональную роль ("закон сердца" Франка-Старл инга).

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

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

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

• Построение механической модели активации сокращения мышц ионами кальция с учетом спиральной формы тропомиозина в недеформированном состоянии.

• Исследование влияния параметров модели на характер активации мышечного сокращения ионами кальция.

Научная новизна. В работе получены следующие новые результаты, выносимые на защиту:

• Предложена модель механики саркомера, учитывающая кручение нитей актина. Решена задача о релаксации механических характеристик после

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

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

• Предложена механическая модель, которая рассматривает тропомиозиновую цепь как упругий стержень, имеющий в недеформированном состоянии форму спирали и способный скользить без удлинения по поверхности нити актина. Модель описывает существенные черты Са активации тонких нитей в сердечной мышце и объясняет высокую кооперативность связывания ТпС ионами Са2+ даже в отсутствие миозиновых головок.

• Модель также воспроизводит увеличение Са-чувствительности тонких нитей при присоединении миозиновых головок без существенного увеличения кооперативности. Кроме того, модель объясняет результаты экспериментов, в которых некоторые молекулы ТпС были заменены мутантами, нечувствительными к Са2+.

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

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

Апробация работы. Результаты, полученные в работе, были доложены на следующих научных конференциях: Конференция-конкурс молодых ученых НИИ механики МГУ (2009, 2011, 2012); 10-ая Всероссийская конференция по биомеханике "Биомеханика 2010" (Саратов, 2010), Международная конференция "Биологическая подвижность: от фундаментальных исследований до нанотехнологий" (Пущино, 2010), Конференция МГУ "Ломоносовские чтения" (2012), Международная конференция "Биологическая подвижность: фундаментальные и прикладные исследования" (Пущино, 2012), 41-ая Европейская мышечная конференция (Родос, Греция, 2012).

Результаты работы были также доложены и обсуждены на специализированных научных семинарах: семинар по биофизике под руководством проф. Роджера Кука (Roger Cooke, UCSF, Сан-Франциско, Калифорния, США, 2012), семинар по механике сплошных сред под

руководством акад. РАН А.Г. Куликовского, проф. В.П. Карликова, члена-корр. РАН О.Э. Мельника (НИИ механики МГУ, Москва, 2012, 2013).

Публикации по теме диссертации. Основные результаты работы изложены в 10 научных публикациях, из которых 2 - статьи в журналах, входящих в Перечень ВАК на момент публикации, и 8 - тезисы докладов.

ГЛАВА 1. ОБЗОР ЛИТЕРАТУРЫ

1.1 Строение поперечно-полосатых мышц

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

Рассмотрим строение скелетной мышцы (см., например, (Кубасова и Цатурян, 2011)). Она состоит из пучков параллельно упакованных мышечных клеток, или волокон (Рис. 1.1). Характерный диаметр клетки 50-100 мкм, диапазон возможных длин очень широк. Каждое волокно представляет собой гигантскую многоядерную клетку, в которой ядра отнесены к поверхности, и до 80% объема занимают цилиндрические образования - миофибриллы. Диаметр отдельной миофибриллы порядка 1 мкм. Миофибрилла состоит из повторяющихся структурных образований - саркомеров, длина которых составляет 2-2.5 мкм. Благодаря этому, в миофибриллах скелетных мышц при большом увеличении микроскопа можно наблюдать правильное чередование более светлых и более темных участков. Так как отдельные миофибриллы расположены в волокне параллельно друг другу, а границы саркомеров соседних миофибрилл расположены в плоскостях, параллельных оси волокна, поперечная исчерченность с периодом, равным длине саркомера видна и на уровне целого волокна. Поэтому скелетные и сердечную мышцы называют поперечно-полосатыми. Саркомер ограничен с двух сторон

плотными образованиями, называемыми Z-дисками.

мышечная клетка

,! / : *; . • tч Л Mwlit^''

мышца

МЛ

I •

•> < ■ I

capucwep

JL

.Л I ......р»

■"•'"■■ ..У

Ч:

2 МКМ

15 нм

«юпгиент сапкгмера: ииозлновые гсловкр

тонкая а*ги новая >«ть

толстая рлиоэигаооя нить

Рис. 1.1. Структура скелетной мышцы. Указаны характерные размеры саркомера и головки миозиновои молекулы. Синим, красным и зелёным показаны миозиновые головки, принадлежащие трём толстым нитям, окружающим одну тонкую актиновую нить. Серым цветом показаны мономеры актина, доступные для миозиновых головок, голубым - актиновые мономеры, с которыми миозиные головки не могут связаться из-за присутствия молекул регуляторного белка тропонина (регуляторные белки на рисунке не показаны) (Кубасова и Цатурян, 2011).

В еаркомере присутствуют два различных семейства нитей, толстые и тонкие. Это было обнаружено с помощью электронной микроскопии в 1953 г. (Hanson и Huxley, 1953). Впоследствии было показано, что тонкие нити состоят, преимущественно, из белка актина. Мономеры актина молекулярной массой 42 кД (глобулярный, или G-актин) полимеризуются, образуя фибриллярный F-актин. Тонкая нить актина имеет спиральную структуру (Oda и др., 2009), (Fujii и др., 2010) и (Dominguez и Holmes, 2011), (Рис. 1.1),

которую часто представляют как левую спираль с полным периодом равным 13 единиц расположенных на 6 оборотах. Угол поворота соседних мономеров друг относительно друга примерно 167°, шаг между мономерами, измеренный по оси нити, равен примерно 2,75 нм, таким образом, полный период спирали равен примерно 36 нм (Holmes, Popp и др., 1990).

Тонкие нити поперечно-полосатых мышц имеют в своем составе и другие белки, в том числе гигантский белок небулин, который играет важную роль