Математическое моделирование и оптимальное управление Вестник Нижегородского университета им. Н.И. /Лобачевского, 2007, № 2, с. 186-191
УДК 517.925.54
КОНКУРЕНЦИЯ В ПРОЦЕССАХ САМОСБОРКИ ЛИНЕЙНЫХ ЦЕПЕЙ
© 2007 г. А.С. Мартынова
Нижегородский госуниверситет им. Н.И. Лобачевского martynova@vist.nnov.ru
Поступило вредокцию 19.02.2007
Обнаружено, что конкуренция в самосборке линейных цепей с элементами двух типов может порождать неустойчивость состояния равновесия. Рассмотрены два способа конкуренции: через интенсивности образования и через интенсивности разрыва связей. Гиперцикл конкуренции, напротив, сохраняет устойчивость состояния равновесия. Полученные результаты отражают некоторые особенности развивающихся систем различной природы.
Постановка задачи
Пусть имеются два типа элементов, обозначим их как тип 1 и тип 2, и достаточно много неразличимых копий каждого из них равномерно заполняет некоторый замкнутый объём. Независимо перемещаясь и сталкиваясь между собой, копии элементов могут взаимодействовать друг с другом, образуя направленные связи. После образования связи между элементами появляются компоненты или линейные цепи из двух звеньев, составленные из этих элементов. Эти цепи, в дальнейшем при взаимодействии с другими элементами, могут расти до сколь угодно длинных цепей. Будем считать число копий всех элементов настолько большим, что при взаимодействиях их друг с другом они могут образовывать бесконечные цепи. Цепи могут быть сколь угодно большими, но никогда не могут превратиться в замкнутые кольца, то есть зациклиться (рис. 1).
-0
0
i +1 i + 2
] -1 ]
Рис. 1
Будем предполагать, что образовавшиеся связи между элементами цепи не являются «вечными» и при самосборке возможны разрывы в цепях. Это предположение недопустимо для биополимеров, таких как ДНК и РНК, но является справедливым для химических полимеров. Однако если в полученной для химических полимеров модели принять параметры, описывающие разрывы в цепях, равными нулю, то получившаяся модель будет актуальна для описания и биополимеров.
Образования и разрывы связей между элементами в цепях будем описывать с помощью интенсивностей образования и разрыва связей -р и д. Отметим, что интенсивности образования и разрыва связей являются аналогами скоростей реакций в химических моделях.
Под конкуренцией в линейных цепях будем понимать такой тип взаимодействий между элементами самосборки, при котором:
- рост концентрации одного типа связи препятствует появлению или сохранению других типов связей;
и/или
- скорость появления одного типа связи увеличивается при увеличении концентрации этой связи в самосборке.
Последнее условие сродни условию автокатализа, которое заключается в том, что константа скорости появления некоторого вещества увеличивается при увеличении концентрации этого вещества в химической реакции [6].
Термин «конкуренция» введен из-за схожести процесса самосборки с рассмотренным выше типом взаимодействий - с математическими моделями естественного отбора [1]. Два различных условия существования конкуренции в самосборке дают два различных способа ее моделирования - через интенсивности разрыва д и через интенсивности образования связей р .
Отметим, что конкуренция является одним из видов объемных взаимодействий процессов самосборки, но отдельное ее исследование позволяет наглядно интерпретировать получаемые результаты [4].
Модель марковского процесса самосборки линейных цепей была предложена и исследована М.Л Таем [3]. Им же было найдено интегральное многообразие. Для моделирования конкуренции в линейных цепях, в рамках уже
(1)
существующей математической модели, достаточно лишь изменить вид интенсивностей образования и разрывов связей. Поэтому дальнейшее исследование динамики этого процесса будем проводить на найденном интегральном многообразии:
у11 = (—1 -у11 -уп) х х Ри( А - У11- У21)- ч- У11; у12 = (D1 - У11 - У12) х х р12(—2 - у12 - у22) - д12у12 ;
У21 = (D2 - У21 - У22) х х р21(D1 - у11 - у21) - д21 у21 ; у22 = (D2 - у21 - у22) х х Р22(D2 - У12 - У22) - 422У22 •
Отметим, что конкуренция (можно назвать ее материальной) за ресурсы (элементы самосборки) уже присутствовала в процессе самосборки линейных цепей. Действительно, объем, в котором происходит процесс самосборки, является замкнутым, а число копий всех элементов пусть и достаточно большое, но всё же ограничено. Однако материальная конкуренция не приводит к смене устойчивости состояния равновесия на интегральном многообразии [3] .
Конкуренция между гомо- и гетеросвязями
Гомосвязью будем называть связь, образованную между двумя элементами одного и того же типа, а связь между элементами разных типов назовем гетеросвязью. Таким образом, связь, образованная между элементами в компоненте 11 или 22, будет гомосвязью, а связь в компоненте 12 или 21 - гетеросвязью. Подобная же терминология используется и в химии: например, связи между одинаково и разноименно заряженными ионами называют гомо- и гетеро-полярными [6].
В случае конкуренции между гомо- и гетеросвязями процесс самосборки линейных цепей из двух элементов в фазовом пространстве
ф = {(У1ЬУ12,У2ЬУ22): 0 < У11 +
+ у12; у11 + у21 < А’0 < у12 + (2)
+ У22; У21 + У22 < —2 }
будет описываться с помощью следующей динамической системы:
У11 = (—1 - у11 - у12)(А - у11 - у21) х х р11(1 + а11у11) - д11(1 + в12у12 + в21у21)у11 ;
у12 = (А1 — у11 — у12)(А — у12 — у22) Х Х р12(1 + а12у12) — Ч12(1 + в11у11 + в22у22)у12 ; у21 = (В2 — у21 — у22)(А — у11 — у21) Х Х р21(1 + а21у21) — Ч21(1 + в11у11 + в22у22)у у22 = (В2 — у21 — у22)(В2 — у12 — у22) Х Х р22(1 + а22у22) — Ч22(1 + РиУ 12 + в21у21)у
(3)
21
22
где константы рп, Р12 , Р21 , Р22, ®1Ь «12,
а21, а22 > 0 , Чи, 412, 421, 422 , Аь в12, Р21, Р22 > 0 . Подчеркиванием сверху здесь и всюду ниже будем различать переменные, через которые моделируем конкуренцию.
Справедлива следующая теорема.
Теорема 1. Если выполняются равенства
А = А = А, ап = «12 = «21 = «22 = а, в11 = в12 = в21 = в22 = в, р11 = р12 = р22 = р21 = р, Ч11 = Ч12 = Ч22 = Ч21 = 4 ,
(4)
тогдо у системы (3) в фозовом простронстве
(2) существует симметричное состояние ров-
новесия (у, у, у, у), где 0 < у < —, которое
2
будет устойчивым, если
ав < —-у ,
2 у2
и неустойчивым при норушении условия (4).
(5)
Доказательство. 1) Для того чтобы (у, у, у, у) было состоянием равновесия, нужно, чтобы оно удовлетворяло уравнению
(— - 2у )2 р(1 + а у) = 4 у (1 + 2в у). (6)
При выполнении равенств (3) и при у^ = у все уравнения системы (2) примут вид:
уу = (— - 2у )2 р(1 + а у ) - 4 у (1 + 2в у ) , (7)
что и доказывает существование симметричного состояния равновесия у системы (3).
2) Докажем теперь устойчивость состояния равновесия (у, у, у, у).
Исследуем поведение интегральных кривых системы (3) в окрестности состояния равновесия: линеаризуем ее в окрестности (у, у, у, у) и составим характеристическое уравнение. Оно будет иметь вид:
Л — а Ь — —Ь 0
— Ь Л—а 0 —Ь
Ь — 0 а — Л Ь —
0 —Ь Ь — а — Л
= 0, (8)
где
— (Р — 2 у)(Р + 2ау 2)
а = — р----------------^< 0,
У
ь = — р(Р — 2 у )(1 + а)(1 + /р) < о
(1 + 2в У) ■
Значения а и Ь получены из дифференцирования правых частей системы (3) при подстановке значений производных в состоянии равновесия (у, у, у, у), с учётом равенств (4).
Упростим определитель в (6): вычтем последний столбец из первого, а третий столбец из второго. Затем к последней строчке прибавим первую строчку, а к третьей строчке вторую. Получим:
= 0.
Отсюда получаем
А(Л) = (Л — а )2((Л — а )2 — 4Ь2) = = (Л — а)2 (Л — а — 2Ь)(Л — а + 2Ь) = 0 .
Видим, что три первых характеристических корня будут строго отрицательными. Осталось только проверить знак последнего корня.
Л — а 0 Ь — 0
0 Л—а 0 — Ь
0 0 Л—а — 2Ь
0 0 — 2Ь Л—а
Лд — а — 2Ь — —
р( Р — 2 у)
зт-[(Р + 2а у2) х
у(1 + 2в у) х (1 + 2в у) — 2у (1 + ау )(1 + вР)] = р(Р — 2у) х [(р — 2у) — 2а/ у2 х
Замечание 1: Условия (5) автоматически выполнятся при а = 0 или в = 0, и, следовательно, в случае, когда выполнено лишь одно из условий конкуренции.
Замечание 2: Следует отметить, что выполнение условия (5) зависит от величины у в состоянии равновесия у. Это связано с большим числом параметров в системе (3), и этого можно избежать, если дополнить (3) необходимыми значениями параметров. С другой стороны, на у в условии (5) можно смотреть как на новый параметр системы (3), и он, к тому же, может быть найден численно на ЭВМ.
Кроме этого, справедливо
Следствие. При выполнении равенств (3) в
2
области параметров ав >> —- (5) у системы
Р 2
(3) существует неустойчивое симметричное состояние равновесия (у, у, у, у), где 0 < у < Рр.
Итак, аналитически обнаружили, что смена устойчивости симметричного состояния равновесия в процессе соперничества между гомо- и гетерополимерными структурами при «равных условиях борьбы» (4) будет происходить только при достаточно «сильной» конкуренции, выраженной и через интенсивности образования р, и через интенсивности разрыва д. В ходе численного эксперимента на ЭВМ (рис. 2) обнаружено, что нарушение устойчивости повлечет за собой появление в фазовом пространстве (2) еще двух состояний равновесия. Они будут устойчивы, и в одном из них преобладает концентрация гомосвязей, а в другом - гетеросвязей. Процесс самосборки линейных цепей в этом
у(1 + 2/ у)
х (Р — 2у)] = — Р(Р —/1(1 - 2 а/ у2) .
у(1 + 2/ у)
Видим, что в случае, когда (1 — 2 а/ у2) > 0 или а/ < —-^т, Л4 < 0, состояние равновесия
2 у 2
будет устойчивым. Однако когда а и / будут
а 1 2
строго положительны и а/ > т > —т , тогда
2 у Р2
Л4 > 0 и состояние равновесия будет неустойчивым.
Теорема 1 доказана.
Рис. 2. Фазовый портрет системы (3), полученный при следующих значениях параметров: Р = 1,а = 6,
в = 4, р = 1, д = 1. Крестиками отмечены два устойчивых состояния равновесия, биссектрисой -интегральное многообразие, на котором находится неустойчивое симметричное состояние равновесия
случае будет приходить к одному или другому «лидерству» в зависимости от выбора начальных условий. Нечто подобное можно увидеть в утверждении одной из гипотез о превосходстве существующего, единого для всей биосферы, генетического кода [5]. В этой гипотезе принимается, что все варианты кода были равноправны (условия (4)), но в результате взаимодействия между разными популяциями «случайно» (за счет выбора начальных условий) был выбран существующий сейчас код. Также существует другая точка зрения - гипотеза, утверждающая, что запоминание единого кода произошло за счет его превосходства над другими вариантами кода [7]. Другими словами, выбор кода стал результатом естественного отбора в эволюции. К сожалению, ни одну из этих гипотез невозможно проверить, т.к. для этого требуется провести очень длительный по времени эксперимент (биологическая эволюция предположительно длилась около 109 лет).
Г иперцикл конкуренции
Теперь будем считать, что все концентрации связей конкурируют между собой. Причем рост концентрации связи 11 приводит к увеличению интенсивности разрыва связи 12, а рост концентрации связи 12 - к росту интенсивности разрыва для связи 21. В свою очередь, рост концентрации связи 21 приводит к росту разрывов связей 22, рост концентрации которых увеличивает интенсивность разрыва связей 11. Схематично эти взаимодействия между концентрациями связей можно изобразить так:
Рис. 3
Из рисунка 3 видно, что теперь в самосборке устанавливаются конкурирующие взаимодействия, имеющие замкнутую (циклическую) структуру. Биофизик М. Эйген назвал подобную систему взаимодействий гиперциклом [7, 8]. Поэтому, по аналогии, назовем рассматриваемые в самосборке линейных цепей взаимодействия гиперциклом конкуренции.
Следует заметить, что гиперцикл конкуренции имеет существенное отличие от «классического» гиперцикла Эйгена - Шустера. «Класси-
ческий» гиперцикл возможен только для само-генерирующихся динамических систем, когда компоненты гиперцикла участвуют в синтезе друг друга. А гиперцикл конкуренции заключается в обратном - компоненты взаимодействия участвуют в разложении друг друга.
Перечислим основные преимущества гиперцикла конкуренции перед описанной выше конкуренцией:
— Г иперцикл поддерживает «оптимальную» конкуренцию между компонентами самосборки, когда «преимущество одних не идет в ущерб другим», и позволяет всем структурам развиваться одновременно.
— Гиперцикл объединяет все связи при самосборке в единую кооперативную систему, способную к слаженному развитию, когда преимущества одной связи могут быть использованы для развития остальных связей.
— Гиперцикл - это одна из разновидностей принципа обратной связи, позволяющей системе динамически «приспосабливаться» к условиям окружающей среды.
В случае гиперцикла конкуренции процесс самосборки линейных цепей в фазовом пространстве (1) будет описываться следующей динамической системой:
у11 = (Р1 — у11 — у12) р11(Р1 — у11 — у21) —
— д11у11(1 + в22 у22);
у12 = (Р1 — у11 — у12) р12( Р2 — у12 — у22) —
— д12 у12(1 + в11уиУ; у21 = (Р2 — у21 — у22) р21( Р1 — у11 — у21) —
— д21у21(1 + в12 у12);
у22 = (Р2 — у21 — у22) р22( Р2 — у12 — у22) — д22 у22(1 + в21у21)
где константы р11 , р12 , р21 , р22 , Я11 , д12 ,
Я21 , ^22, в11 , в 12, в22 строго положительны.
Справедлива следующая теорема.
Теорема 2. Если выполняются равенства
Р = Р2 = Р, в11 = в12 = в21 = в22 = в; р11 = р12 = р21 = р22 = р , ----- (8) д11 = д 12 = д21 = д22 = тогда у системы (7) в фазовом пространстве (2) существует устойчивое симметричное состоя,_________________________, п - Р
ние равновесия (у, у, у, у), где 0 < у < —.
(7)
Доказательство. 1) Для того чтобы (у, у, у, у ) было состоянием равновесия, нужно, чтобы оно удовлетворяло уравнению
(Р — 2у)2 р(1 + а у ) = ду (1 + в у ). (9)
При выполнении равенств (8) и при у^ = у все уравнения системы (7) примут вид
уи = (Р — 2у )2 р(1 + а у) — ду (1 + в у), (10)
что и доказывает существование симметричного состояния равновесия у системы (7).
2) Докажем теперь устойчивость состояния равновесия (у, у, у, у). Для этого исследуем поведение интегральных кривых системы (7) в окрестности состояния равновесия. Линеаризуем ее в окрестности (у, у, у, у) и составим характеристическое уравнение. Оно будет иметь вид:
А(Л) =
Л-а -Ь -Ь
- с
- Ь-с Л-а 0 -Ь
-Ь -с Л-а -Ь
0 -Ь - Ь-с Л-а
= 0,(11)
где
_D(Р - 2у) _
а = - р—-------— < 0,
с = - р
У
Г Р - 2 У )2 в
Ь = -р (D - 2у) < 0,
(1+в У)
< 0.
Значения а , Ь и с получены из дифференцирования правых частей системы (4) при подстановке значений производных в состоянии равновесия (у, у, у, у ), с учётом равенств (8).
Упростим определитель в (11), прибавив все его столбцы к последнему:
А(Л) = (Л - а - 2Ь - с)А1 (Л) = 0,
где
А1(Л)=
Л-а -Ь -Ь 1
- Ь -с Л-а 0 1
-Ь -с Л-а 1
0 - Ь - Ь - с 1
Определим знак первого характеристического корня:
Л = а + 2Ь + с = - -РР—2у) [(Р + 2 у) х У (1 + в У )
х (1 + в У) + У(D - 2У)в] = - IУ) х
У (1 + в У )
х [(D + 2У) + в У (D + 2У + D - 2У)] =
[(D + 2 У) + 2 У в D ] < 0.
У (1 + в У)
Получили, что Л1 < 0.
Преобразуем А1 (Л): в^ічтем последнюю
строчку из всех остальных строк:
А1(Л) =
Л - а 0 с
— Ь — с Л — а + Ь Ь + с
— Ь Ь — с Л — а + Ь + с
. (12)
Преобразуем определитель в (7) следующим образом: вычтем из третьей строчки первую, а затем прибавим к первому столбцу последний. В результате этих действий получим:
А1(Л) =
Л — а + с 0 0
0 с
Л-а+Ь Ь+с
Ь — с Л — а + Ь
= (Л — а + с) • А2 (Л),
где
А 2(Л) = (Л2 + R1Л + R2).
Rl = Ь - а,
^2 — а — 2аЬ + с — а (а — 2Ь) + с . Теперь можем определить знак второго характеристического корня:
Л = а — с=— р (Р - 2 у) х
у (1 + в у)
х [Р(1 + в у) — ув(Р — 2у)] =
р(Р — 2у) [р + 2ву2] < 0.
2
2
у (1 + в у)
Получили, что и Л2 < 0 .
Определим условия устойчивости полинома А 2 (Л) - условие отрицательности его характеристических корней, применив к нему критерий Раусса - Гурвица. Для выполнения критерия необходимо и достаточно, чтобы R1 >0 и R2 >0.
Определим их знаки:
& = Ь - а =
Р(Р - 2У )(Р - У )
У
> 0,
&2 — а(а — 2Ь) + с .
Определим знак выражения в круглых скоб-
ках в &2 : а — 2Ь = —
Р(Р - 2У)2 У
< 0 .
Получили, что все слагаемые в R2 будут строго положительны.
Итак, при всех в > 0 : R1 >0, R2 >0. Это и завершает доказательство устойчивости симметричного состояния равновесия (у, у, у, у).
Теорема 2 доказана.
Замечание: Теорема 2 будет справедлива и
для системы
(1З)
у11 = (Р1 — у11 — у12)р11( Р1 — у11 — у21) —
— д11у11(1 + в22 у22 ); у12 = (Р1 — у11 — у12 ) р12 (Р2 — у12 — у22 ) —
— д12 у12(1 + в11упХ у21 = (Р2 — у21 — у22) р21(Р1 — у11 — у21) —
— д21у21(1 + в12 у12);
у22 = (Р2 — у21 — у22 )р22 (Р2 — у12 — у22) —
— д22 у22(1 + в21у21),
когда гиперцикл конкуренции будет направлен в другую сторону (рис. 4).
Рис. 4
Действительно, с помощью равенств (8) система (13) примет вид (7), и, следовательно, для нее справедливы все выводы, полученные в ходе доказательства теоремы 2.
Численный эксперимент на ЭВМ при несимметричных параметрах взаимодействий показал, что при Р\ = Р2 траектории системы (7) приходят в устойчивое состояние равнове-
сия (уп, у12, у21, у22 ). И в этом состоянии будет преобладать та концентрация связи, у кото-
„ „ Я
рой соотношение интенсивностей — является
р
наименьшим.
Таким образом, получили, что гиперцикл не нарушает устойчивости состояния равновесия соответствующей динамической системы и поддерживает «оптимальную» конкуренцию между концентрациями связей, когда преимущество одних концентраций связей не приводит к уменьшению других концентраций.
Список литературы
1. Вольтерра В. Математическая теория борьбы за существование. - М.: Наука, 1976.
2. Леонтович А.М. Одна задача о самосборке отрезков // Проблема передачи информации. 1975. № 15. С. 97-106.
3. Тай М.Л. Исследование стохастической модели самосборки линейных цепей // Проблема передачи информации. 1979. Т. 15. № 4. С. 40-52.
4. Тай М.Л. Процессы самосборки с объемными взаимодействиями // Докл. АН СССР. 1990. Т. 314. № 1. С. 189-193.
5. Романовский Ю.М., Степанова Н.В., Чернав-ский Д.С. Математическая биофизика. - М.: Наука, 1984.
6. Эмануэль Н.М., Кнорре Д.Г. Курс химической кинетики. - М.: Высшая школа, 1984.
7. Эйген М. Самоорганизация материи и эволюция биологических макромолекул. - М.: Мир, 1982.
8. Эйген М., Шустер П. Гиперцикл. Принципы самоорганизации макромолекул. - М.: Мир, 1982.
COMPETITION IN THE PROCESSES OF SELF-ASSEMBLY OF LINEAR CHAINS
A. S. Martynova
It is found that competition in self-assembly of linear chains with two types of elements can generate instability of equilibrium. Two ways of competition are considered: via intensities of formation and via intensities of breakup of relationships. On the contrary, the hypercycle of competition preserves the stability of equilibrium. The results of this research reflect some features of developing systems of different nature.