Термодинамика многоядерных клеток: системное моделирование самоорганизующегося саркомера с хаотической динамикой параметра порядка Текст научной статьи по специальности «Физика»

ИНФОРМАЦИОННЫЕ ТЕХНОЛОГИИ

Г. П. Быстрай, А. В. Богинич

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

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

Введение

Мышечное волокно представляет собой сложную многоядерную клетку, содержащую одну-две тысячи более тонких вытянутых волоконец (миофибрилл) диаметром 1-2 мкм, состоящих из элементарных сократительных единиц — саркомеров, длина которых в покоящейся мышце составляет ~2,2 мкм, т. е. это достаточно протяженные объекты, для моделирования деформации которых можно использовать феноменологический подход [1]. Толстые и тонкие нити саркомеров образованы из сократительных (миозин и актин) и Ca+2 — чувствительных регуляторных белков.

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

Общепринято считать, что в процессе цикла укорочения-удлинения саркомера нити сократительных белков актина и миозина не изменяют свою длину, а лишь скользят относительно друг друга [2, 3]. Имеются многочисленные экспериментальные доказательства того, что внутренние напряжения в саркомере формируются в результате химических реакций с участием Ca, Mg, ATP и др. [4].

Проблема моделирования сокращений саркомера в условиях скольжения усугубляется следующим. Экспериментально [5, 6] зафиксирована самовоз-буждающаяся хаотическая динамика сокращений, что, вероятно, и обусловливает дрожание фибриол и мышц в целом. Вопросы нелинейного поведения мышц (саркомеров) в нестационарных условиях остаются открытыми [1]. Нелинейные процессы, в отличие от линейных, имеют несколько стационарных состояний, которые должны обусловливать ступенчатый характер деформации саркомеров [7, 8]. При значительных деформациях таких ступенек должно

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

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

Динамика линейного сокращения саркомера

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

** = ,

где <зл — напряжения (поток импульса), — величина деформации (термо-

динамическая сила), ц — модуль упругости. Законом (1) описываются малые линейные трехразмерные деформации (по осям х, у, г). Постулируется, что развиваемые напряжения не зависят от поперечной координаты саркомера. Тогда динамическое уравнение одноразмерного продольного сокращения относительно по оси растяжения в линейной задаче в условиях действия внешнего напряжения а примет вид

Шг д0

— = -ф----, (2)

Ш дг

где 0 = -ае г +1 иг2,

е 2

следовательно:

1 Шг

= -иг + а; (3)

ф ш

здесь г — модуль относительной одноразмерной деформации вдоль протяженной цепи полимера, О — потенциальная функция, упругая составляющая которой пропорциональной квадрату деформации и работе внешних сил при растяжении, ф — некоторая размерная константа. Это означает, что рассогласование внутренних напряжений а = иг и внешних ае определяется скоростью деформации ёг / Ш с точностью до постоянной ф.

Будем предполагать, что в случае вязкоупругих напряжений саркомер с малой скоростью одноразмерных деформаций при а «ае подчиняется уравнению Максвелла [9]

Ша•

а +х—'- = Пог, (4)

т

где т — время релаксации внутренних напряжений саркомера, п0 — коэффициент динамической вязкости. При моделировании на основе (4) предполагается, что саркомер по своим деформационным свойствам является промежуточным между ньютоновской жидкостью (т/ Лt << 1) и твердым телом (т / At >> 1), подчиняющимся закону Гука (1) [9]. В последнее время процессы, описываемые уравнениями типа (4), получили название локально-неравновесных [13].

Связь с классической теорией. В литературе описывается несколько математических моделей мышечного сокращения [1-3]. В теории Хаксли и ряде других моделей постулируется линейная зависимость между силой и длиной (координатой мостика) саркомера. Такая зависимость следует из того, что данные модели опираются на классический закон Гука, который в свою очередь и дает линейную зависимость между силой и величиной упругой деформации.

В кинетической теории В. И. Дещеревского [3] рассмотрена модель работы мостика с трехстадийным кинетическим циклом, включающим одно замкнутое состояние а мостика и два разомкнутых — тянущее п при х > 0 и тормозящее т при х < 0 (рис. 1).

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

между стадиями цикла (к1, к2, /5) не зависят от координаты мостика.

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

Шп , . . п

—г- = к,(а - п - т) - —,

Ш 1V ' 5’ (5)

Шт п , ^

ш-=---------к2 т,

Ш 5 2

где а — полное число доступных для замыкания мостиков при длине 0,5 саркомера [1, 3], — скорость скольжения нитей, 5 — длина зоны, в которой

мостик развивает тянущую силу, к1 — среднее значение константы скорости замыкания свободных мостиков, к2 — константа скорости распада тормозящих мостиков.

Рис. 1. Трехстадийный кинетический цикл модели В. И. Дещеревского

Покажем, ввиду важности, что линейная модель малых деформаций (3) и

(4), являющаяся базовой для построения нелинейной модели саркомера, соответствует модели Дещеревского (5). Будем считать, что величина деформации пропорциональна числу тянущих мостиков s = L • n (L = const — де-

формационный параметр одного мостика). Тогда первое уравнение системы

(5) примет вид

^ = _^к +у^ + Lkl(a-т). (6)

Из сравнения (6) и (3) следует, что между модулем упругости До и константами рис. 1 существует взаимозависимость

Л ЬкЛа-т)

Д0- = ^ +-^ , а = —----------. (7Х (8)

Рассмотрим случай малой скорости деформации в , в котором внутренние и внешние напряжения равны ае ^аг. = а. Выражение (8) показывает, что с

увеличением числа тормозящих мостиков т напряжения уменьшаются. Таким образом, уравнения (3) и (6) идентичны с точностью до приведенных коэффициентов. Выразим т из (8): т = а-а-/Ьк1. Учитывая временные зависимости внутренних напряжений а = а(г) и числа тормозящих мостиков т = т(/) в саркомере, получаем

йт - йа йа

= 0. (9)

Подставим (9) во второе уравнение системы (5), а также учтем, что в соответствии с уравнением (3) справедливо уравнение

1 • о

8 =------8 +------. (10)

^0 ^0

В результате получаем уравнение типа уравнения Максвелла для внутренних напряжений, несущественно отличающееся от (4):

йо *

о + т— = /(По)8 + А . (11)

ш

Здесь т = ф / ф — время релаксации внутренних напряжений, /(п0) = Вк1 /Д05фф — функция коэффициента вязкости, А = к1к2аЬ/- — константа, - = д0к2 / к1 + $ / 5 . Если к1к2аЬ <<ф, то константой А можно

пренебречь, что и приводит к уравнению (4).

Такая модель соответствует следующим представлениям. Пусть сарко-

мер имел некоторую первоначальную длину, которой соответствовало неко-

торое число тормозящих мостиков. Если в момент времени t = 0 скачком изменить его длину /ь растянув саркомер или позволив ему сжаться до длины /2, то этому будет соответствовать другое число тормозящих мостиков. Из (5) и (10) следует, что количество тормозящих мостиков, а с ним и напряжения будут стремиться к своему новому значению по экспоненциальному закону:

—к2г —г / т

т = т0е 2 , о = о0е .

В итоге мы показали, что система линейных уравнений (3), (4), описывающая линейные деформации и релаксацию вязких напряжений в саркоме-

ре, эквивалентна системе уравнений Дещеревского (5), в которой релаксация напряжений напрямую не заложена.

Динамика нелинейного сокращения саркомера

Нелинейная модель. Известно, что свободная энергия гидролиза АТФ преобразуется в энергию упругой деформации молекул фермента, которая затем может быть использована для совершения механической работы. Описываемый процесс динамики мышечных белков при значительных (средних) деформациях не является линейным, поэтому при моделировании мы не можем использовать линейный закон Гука. Для нелинейных процессов одноразмерной деформации по оси растяжения представим, следуя [9], в виде полинома

д(в) = до(1 -кв + %в2), (12)

где д0 — коэффициент упругости для линейных систем, а к и X — некоторые коэффициенты [9], характеризующие зависимость модуля упругости в направлении оси растяжения от величины деформации.

Подставив выражение для д (в) в исходное динамическое уравнение (3),

получим для саркомера термодинамическое уравнение с кубической правой частью

Жв /3т 2 ч

— = -ф(ДоХв -Д0кв + Д0В-О , (13)

ж

или, введя в модель потенциальную функцию О — свободную энергию, связанную с нелинейной упругой деформацией, получим

Жв д0 _ д0хв4 д0 кв3 1 2

— = -ф---, 0 = —-------—-----+— д 0в2-ав. (14)

Ж дв 4 3 2 0

Линейная часть этого потенциала содержит квадратичную функцию. Для стационарных деформаций в результате из (15) следует нелинейный закон деформации саркомера

ае =Д0в-Д0 кв2 +Д0Хв3, (15)

в котором внешние напряжения ае и внутренние (правая часть (15)) равны

друг другу. При малых деформациях (в линейном случае) из (15) получаем закон Гука.

Для уравнения (14) можно выделить некоторую критическую точку, в которой происходит смена режима работы саркомера. Использовав значения переменной и параметров в этой точке, можно перейти к относительному виду уравнения (14) [14]. Для этого умножим левую и правую части дифференциального уравнения (14) на (2ддв3 )-1. Получаем

Жв* ( 3 1- 2 1 Л

*3 к *2 1 * *

в--------------в +----------2в -а

V Хвс Хвс у

(16)

где вс — значение модуля деформации в некоторой критической точке

* Л

в = вс (в = 1). Тем самым мы ввели масштаб деформации. Приведем следующие обозначения:

* £ * 2 t * О О

s = . t = Д0ХФМ = —^ t. О =------------3 =— ■

£с t0 №с Ос

Величина ос =д0в с характеризует значение обобщенного внешнего поля в

этой критической точке. Чтобы перейти к канонической форме, необходимо переобозначить переменную и параметры уравнения (8), что приведет к исчезновению части членов (в соответствии с теоремой Тома [10, 11]); в данном случае исчезнет квадратичный член [8, 12]. Следуя [8], перейдем к новой пе-

* *

ременной п = £ - £0 , характеризующей отклонение приведенной величины деформации саркомера s* = s / sc от некоторого среднего значения

*

е0 = е0/еc = k/3xsc ■ Тогда уравнение (16) при условии равенства в критической точке a = b = 0 и, следовательно, п = 0 примет вид

— = -(П3 + ац + Ъ*), (17)

dt

где канонические управляющие параметры a и b представлены выражениями:

а* = -3(s02 -1) , Ъ = -о* + 3s0 - 2s03 ■

В нашей модели параметр Ъ* является параметром рассогласования при**

- 2s 0 ^&0

* * * *3

ложенного внешнего поля о и собственного поля саркомера о. = 3s0 - 2s

1 * * *

Ь = -а + а.

Как одно из следствий нелинейной модели получаем, что при продвижении нитей навстречу друг другу в результате скачка длины п в последовательном упругом элементе замкнутого мостика возникает упругая деформация (сжатие):

* * * з

а = а + а п + п . О8)

Особенностью уравнения (18) является то, что в его трижды вырожденной

* * *-1

критической точке п =а =Ь =0 (в = в0 =а = 1) потенциал из одноямного

становится двухъямным, и рассматриваемая система переходит в режим бистабильного поведения [10]. Таким образом, имеет место потеря устойчивости, но и всегда есть возможность стабилизации в одном из минимумов потенциала. Именно такие уравнения (18) и изучаются в теории катастроф и нелинейной динамике [9, 10].

Для двухъямного потенцила (а* < 0) имеем для указанных знаков Ь* нелинейную модель деформации саркомера и ее асимптотические приближения — линейную модель и модель Дещеревского (табл.).

Следовательно, одно замкнутое состояние саркомера разбивается на два состояния в зависимости от знака развиваемой силы — тянущее при

П = в* -в0 > 0 и тормозящее при п = в* -в0 < 0 (рис. 3). Значение в0 определяется следующим образом: для стационарного состояния полагаем Ь* = 0 , тогда

п(п2 + а*) = 0 . (19)

*

Рис. 3. Напряжения а , развиваемые замкнутыми мостиками в саркомере, в зависимости от удлинения (укорочения). При а* = 0,02 имеем п1 = >/- а* =

= 0,141, п2 = -V- а* = 0, 141. Связь с моделью Дещеревского: п1 =5 / вс соответствует в = 5 — началу тянущей зоны; п1 = -5 / вс соответствует в = 5 —

началу тормозящей зоны

Решив (19), получаем два ненулевых корня

V* *

- а , П2 = -V- а .

Подставив решения в п = в - в0, получаем два выражения

* */** * / *

в, =в0 +>/-а , в2 =в0 -у- а , из сложения которых находим полезную

формулу для s0 как некоторого среднего значения деформации саркомера:

Таким образом, среднее значение деформации саркомера определяется как среднее арифметическое между деформациями в минимумах потенциала

при Ь* = 0. Следует отметить, что для локально-равновесных процессов, в отличие от локально-неравновесных [13], соответствие между параметром по-

•

рядка, потенциальной функцией и скоростью изменения параметра порядка п

для любых моментов времени всегда однозначно.

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

* *

динамических уравнений для параметра порядка п = е -ео и приведенной

скорости изменения энтропии в*^ + т г). Последняя является знакопеременной потенциальной функцией, и для нее справедливо градиентное уравнение в случае локально-неравновесных процессов

здесь тг =тг /10 — время релаксации скорости изменения энтропии, или время релаксации потока одного из самых длительных неравновесных процессов. В системе уравнений (21), (22) уже два параметра порядка, п(0 и

в*(^). Представление релаксации потенциальной функции в форме (22) для саркомера является в модели для локально-неравновесных процессов более общим, чем просто релаксация внутренних напряжений. Это связано с тем, что в реальных процессах вязкоупругости релаксируют не только как в уравнении Максвелла поток ст (внутренние напряжения), но и термодинамическая сила — величина деформации е. Поэтому такое представление, т. е. приближение в рамках локального неравновесия, является более полным, чем уравнение Максвелла (4).

Дифференцировав (22) по п и подставив полученное выражение в (21), получаем дифференциальное уравнение второго порядка для локальнонеравновесных процессов в саркомере

4 2

Сц = дв*

(21)

С дп

(22)

(такие процессы характеризуются временем релаксации т [13]). Здесь f — обобщенная сила саркомера, потенциальная энергия которого дается двухъ-

ямным потенциалом О*(Т) [8]. Член с тг можно не учитывать, когда время релаксации потенциальной функции О (изменения энтропии) существенно меньше времени действия внешних сил тг / Ат << 1. Следует обратить внимание на то, что в нелинейном уравнении (15) сила f, параметр порядка п,

дп д 2п

скорость его изменения — и член тг —— определены в один и тот же мо-

дt дt

мент времени t. Это означает, что эффект последействия в такой модели отсутствует.

*

Можно рассмотреть частный случай, когда внешняя сила а изменяется по гармоническому закону. Предположим, что сила, действующая на актин, равна по величине и противоположна по знаку силе, действующей на миозин. Вероятно, эту силу (напряжения) можно связать, например, с возникновением и исчезновением кулоновского заряда, который взаимодействует с зарядом, закрепленным на актине. Последнее означает, что управляющие параметры в (23) можно представить в условиях периодического действия кроветока в виде

у * * * *

Ь =-а + а, = -а0ео8ш1,

здесь ш — циклическая частота изменения темпов подачи кислорода с кровью

* * * * а ; при t = 0 а0 = а -аж.

Учет последействия. Если следовать [2], то для саркомеров — достаточно протяженных цепей полимеров — следует учитывать также последействие в связи с конечной скоростью распространения возмущений, связанной в том числе со скольжением нитей белков актина и миозина. При этом если внешнее воздействие представлено в момент времени t -т, то обобщенная

сила / = -дО* / дп задана при последействии в виде /(Т -т) . При этом

дп д 2п

— и инерционный тг—т

дт г дt2

мени t. В результате приходим к уравнению второго порядка для п с запаздыванием:

дп д 2п

+Тг "дт2

диссипативный — и инерционный т г —— члены определены в момент вре-

• + тг ттт = / (т - т) , (24)

здесь т — время последействия. Моделирование последействия в самом простом случае подразумевает разложение силы / (Т -т) в ряд Тейлора; ограничиваясь двумя слагаемыми (не доказывая строго в математическом отношении сходимость такого разложения), представим:

( д/(т) &п + д/ Л

/ (т -т) = / (т) -т

где

- + -

(25)

дО у 3 * ч *

/ (Т) =------= -(п + ап) + а0ео8 шТ. (26)

дп

В системе нелинейных уравнений (24)-(26) содержится уже три параметра порядка, п(0 , G\t) и f (t -т).

Обобщенная диссипативная модель

В результате совместного рассмотрения уравнений (16)-(18) для деформации саркомеров, находящихся в локально-неравновесном состоянии, с запаздыванием при средних деформациях каноническое однородное уравнение второго порядка для величины деформации п [3]:

•• •

Tr n + r(t) П + П + а П = °о cos ®t, (27)

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

r(t) = 1 - т(3п2 + а*) > 0 , а0* = о0(1 + Tratg(rat)).

Параметр т = т /10 — приведенное время последействия, тг = тr /10 — приведенное время релаксации потенциальной функции. Это время определяется наибольшим временем релаксации термодинамического потока или термодинамической силы.

В правой части уравнения (27) стоит периодическая сила, действующая на саркомер, присутствие которой мы связываем с кроветоком. Переменный

параметр y = T(t)/тr является в такой модели параметром диссипации (коэффициентом трения), а описываемый процесс — диссипативным. Справедливость уравнения (27), дающего при некоторых значениях управляющих параметров хаотические решения, означает, во-первых, что установление равновесия в саркомере происходит не монотонно, а в колебательном режиме с собственными частотами хаотических колебаний, определяемыми физическими параметрами саркомера (актин-миозиновой системы). Поэтому в реальном саркомере будет целый спектр собственных частот. Во-вторых, реально мостики нелинейно взаимодействуют с актином, в условиях последействия и релаксации результатом будут хаотические колебания такого маятника, следовательно, в реальном саркомере существуют малые высококоча-стотные пульсации, при которых энергия переходит в тепло. Это означает, что энергия диссипирует в мелкомасштабных пульсациях, как и при турбулентности. При этом малые частоты соответствуют масштабу энергии. Своеобразие заключается в том, что пульсации являются зависимыми от времени — каждое последующее неравновесное состояние зависит от предыдущего.

Результаты экспериментов. Результаты численных расчетов

В работах [5, 6] были проведены эксперименты по сокращению саркомеров. В ходе эксперимента было выяснено, что толстый филамент запирательной мышцы моллюска способен изменять свою длину в диапазоне физических нагрузок. Это свидетельствует о нелинейности данной системы.

Детальный анализ динамики саркомеров продемонстрировал отсутствие гладкости процессов удлинения и укорочения препаратов. Трассы изменения длины саркомеров как на фазе удлинения, так и на фазе укорочения имели хаотическую ступенчатообразную форму (рис. 4а), часто с высокой регулярностью чередования периодов движения и остановок. Размер ступеней не является случайной величиной, он оказался кратным определенной величине и не зависел от направленности изменения длины саркомеров. Конечно, теорети-

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

Рис. 4. Экспериментальные результаты по сокращению саркомера (а), которые сравниваются с теоретической кривой (б): a* = -1,5 , ю = 2,6 ,

т = 0,216 , п(0)=0,3, о0* = 1,8

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

Автономная система уравнений для нелинейных колебаний длины саркомера

Численные решения нелинейного уравнения для п (27), которое представлялось системой трех нелинейных дифференциальных уравнений:

1=Y і

< Y = -—(Г(X)У - H0* 008(7) + (X3 + а*X)), (28)

• т^

7 = ю,

где

X = п, Х = дп, Г(X) = 1 - т(3X2 + а*), Н 0* = Н0*(1 + тш^) , д

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

Рис. 5. Динамика самовозбуждения саркомеров с трансформацией периодических воздействий кроветока на хаотические параметры порядка (а) в модели одноразмерной нелинейной деформации саркомера при наличии релаксации и последействия; фазовый портрет (б): t = ^ = 0,34 с

Странный аттрактор

Как видно из рисунков, при некоторых значениях параметров в решениях уравнения (27) возникает гомо- и гетерофазный хаос, реализующийся по типу странного аттрактора. В такой нелинейной термодинамической системе параметр порядка «мечется» между тремя стационарными состояниями (фазами)

п+

= V- а* , п0 = 0 П- = -V-

а

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

*

(при а < 0, рис. 6). Ниши функции соответствуют состояниям 1 и 2 замкнуто-

1—1 * * го мостика. Параметр порядка п = 8 -£0 характеризует изменение механической координаты для саркомера.

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

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

*

Рис. 6. Хаотическая динамика изменения потенциальной функции О* (производства энтропии) при одноразмерной деформации саркомера с релаксацией и последействием (а, б) и параметра порядка п (в)- Положительные значения О* соответствуют асимптотически устойчивым состояниям по Ляпунову

О* > 0, ёО/& < 0 [15], отрицательные — структурной устойчивости Заключение

При разработке модели предполагалось, что деформация саркомеров как сложных конструктов миозиновых и актиновых нитей не подчиняется закону Гука. Закон деформации для саркомеров выбирался нелинейным, при этом модуль упругости являлся функцией величины деформации- Была выдвинута гипотеза, что саркомеры обладают эффектом вязкоупругости (подчиняются закону Максвелла), когда имеет место релаксация тензора внутренних напряжений со временем релаксации тг. Была установлена связь линейной модели малых деформаций саркомеров с классической теорией, в частности с линейной моделью В. И. Дещеревского. Далее считали эффект последействия для протяженных в одном направлении нитей, когда после начала действия обобщенной силы формирование поля скоростей деформации и поля «ускорений» задерживается во времени. Вводятся время последействия т и возмущающая сила, изменяющаяся по гармоническому закону, которую мы связываем с периодическим кроветоком. В общем случае временной вид возмущающей силы может задаваться, что должно соответствовать экспериментам.

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

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

Решения полученного уравнения реализуются по типу странного аттрактора, отличающегося от странного аттрактора Лоренца увеличенным числом аттракторов: в случае постоянно действующей гармонической силы хаотические решения наблюдаются между тремя ступенями (верхней, средней и нижней). Последнее соответствует экспериментам по сокращению саркомеров [9, 17, 18] при средних деформациях.

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

ЛИТЕРАТУРА

1. Рубин А. Б. Биофизика. — М.: Наука, 2000. — Т. 2. — 467 с.

2. Хилл А. Механика мышечного сокращения. — М.: Мир, 1972. — 182 с.

3. Дещеревский В. И. Математические модели мышечного сокращения. — М.: Наука, 1977. — 150 с.

4. Бендолл Дж. Мышцы, молекулы и движение. — М.: Мир, 1970. — 256 с.

5. Blyakhman F., Tourovskaya A, Pollack G. H. Quanta! sarcomere-length changes in relaxed single myofibrils // Biophysical Journ. — 2001. — August. — Vol. 81. — P. 1093.

6. Nagornyak Е., Blyakhman F., Pollack G. H. Effect of sarcomere length on step size in relaxed rabbit psoas muscle // Journ. of muscle research and cell motility. — 2004. — Vol. 25. — P. 37-43.

7. Быстрай Г. П., Ворох А. С., Андреев С. В. Детерминированный хаос в динамике тока одиночных ионных каналов биомембран // Биофизика. — 2005. — Т. 50, № 5. — С. 851-861.

8. Быстрай Г. П., Пивоваров Д. В. Неравновесные системы. — Свердловск: УрГУ, 1989. — 187 c.

9. Быстрай Г. П., Макаров Л. В., Шилин Г. Ф. Неравновесная термодинамика процессов горного производства. — М.: Недра, 1991. — 119 с.

10. Гилмор Р. Прикладная теория катастроф. — М.: Мир, 1984. — T. 1. — 350 c.

11. Арнольд В. И. Теория катастроф. — М.: Наука, 1990. — 115 c.

12. Быстрай Г. П., Студенок С. И., Иванова С. И. Детерминированная модель го-мофазных и гетерофазных флуктуаций в системе «жидкость-пар» // ТВТ. — 2002. — Т. 40, № 5. — С. 779.

13. Соболев С. Л. Локально-неравновесные модели процессов переноса // УФН. — Т. 167, № 10. — C. 1095.

14. Быстрай Г. П. Метод функций Ляпунова в анализе открытых термодинамических систем // Вестн. кибернетики. — Тюмень: Изд-во ИПОС СО РАН, 2005. — № 4. — С. 122-137.

15. Меркин Д. Р. Введение в теорию устойчивости движения. — М.: Наука, 1987. — 297 c.

THERMODYNAMICS OF POLYNUCLEAR CELLS: SYSTEMIC SIMULATION OF SELF-ORGANIZING SARCOMETER WITH RANDOM DYNAMICS OF ORDER PARAMETER

The article is devoted to analysis of general principles of self-organization under strain of sarcometers in non-linear model with random dynamics of order parameter’s behaviour. It was supposed that sarcometers possess not only the effect of viscoelasticity, but also that of secondary action. As a result, under minor strains for order parameter, a second-order homogeneous non-linear differential equation was obtained, which solutions give a strange attractor under self-excitation. Within the limits of the model, the authors cite evaluation results of its adequacy with the experiment, suggesting a number of criteria for further studying of sarcometers’ systems — myofibrils.