CC BY
30
Прикладные задачи
^^^^^^^^^^»нелинейной теории колебаний и вслн
Изв. вузов «ПНД», т. 18, № 3, 2010 УДК 537.86, 530.182
КОНКУРЕНЦИЯ В ДВУХКОМПОНЕНТНОЙ МОДЕЛИ АНСАМБЛЯ ИММУННЫХ Т-КЛЕТОК
М.В. Иванченко
Исследован процесс конкурентного взаимодействия в двухкомпонетной модели ансамбля иммунных Т-клеток, лежащий в основе селекции наиболее эффективных семейств (клонотипов) Т-клеток. Показано отсутствие периодических колебаний, определены границы между режимами сосуществования, частичного и полного вымирания клонотипов в пространстве параметров модели. Указаны границы применимости приближения среднего поля. Обсуждается биологический смысл полученных результатов.
Ключевые слова: Конкуренция, бифуркации, иммунная система.
Введение
Изучение конкуренции и самоорганизации в больших ансамблях нелинейных элементов со сложными связями является одной из фундаментальных проблем современной теоретической физики. Полученные результаты находят широкое применение в различных прикладных областях: в оптике [1], нелинейных сетях синхронизации [2], химии [3], экологии, экономике и социологии [4,5]. В последнее время методы исследования сложных радиофизических систем успешно распространяются на высоко актуальные биологические задачи, где явления конкуренции и самоорганизации присутствуют повсеместно. При этом преобладающая доля работ приходится на нейродинамику: одним из примеров является теория организации последовательной активности, механизмов памяти и распознавания на основе конкуренции в нейронных ансамблях [6,7].
Проблема регуляции иммунной системы является одним из актуальных вопросов современной иммунологии [8], и в настоящее время интерес к применению методов нелинейной динамики и статистической радиофизики в этой области переживает бурный рост [9]. На фоне различных задач ярко выделяется проблема механизмов саморегуляции численного и видового состава популяции Т-клеток, лимфоцитов, отвечающих за функционирование приобретенного иммунитета. Экспериментальные данные свидетельствуют, что у взрослых людей число Т-клеток находится
вблизи некоторого равновесного значения (порядка 1011-1012 Т-клеток, формирующих порядка 107-108 различных клонотипов). Существуют теоретические и экспериментальное указания о характере такой регуляции: в отсутствие инфекции Т-клетки конкурируют за «стимулы выживания», получаемые от специализированных клеток; полное их число сохраняется вблизи равновесного состояния [10,11]. Однако для надежной работы иммунной системы одного большого числа лимфоцитов недостаточно: необходимо, чтобы в ансамбле сосуществовали различные клонотипы Т-клеток, вызывающие специфическую иммунную реакцию на различные возможные патогены. Вопрос о том, какой механизм обеспечивает равномерное распределение защиты в пространстве патогенов, остается практически неисследованным.
Существенный прогресс был достигнут в [12], где была построена математическая модель, представляющая собой ансамбль взаимосвязанных стохастических процессов рождения и уничтожения Т-клеток различных клонотипов. Было показано, что равномерность покрытия пространства патогенов достигается за счет резкого ускорения вымирания тех клонотипов Т-клеток, которые интенсивно конкурируют за одни и те же области этого пространства.
Исследование многомерного марковского стохастического процесса рождения-уничтожения - исключительно сложная задача, по всей видимости, аналитически неразрешимая. Среди открытых принципиальных вопросов - существование (ква-зи)стационарного распределения вероятности, его характеристики, время вымирания отдельных клонотипов. Результаты работы [12] получены в приближении среднего поля, с помощью которого задача сводится к одномерной модели, что, в свою очередь, существенно ограничивает область применения. В связи с этим особый интерес представляет исследование многомерных динамических и стохастических моделей конкуренции Т-клеток, анализ реализуемых динамических режимов, установление закономерностей различных сценариев конкуренции. Шагом к достижению этой цели является настоящая работа, где рассматривается двумерная динамическая модель конкуренции.
В разделе 1 обсуждается двухкомпонентная математическая модель ансамбля Т-клеток. В разделе 2 изучаются глобальные свойства фазового портрета. В разделе 3 исследуются существование и устойчивость состояний равновевсия. В Заключении формулируются основные результаты исследования.
1. Математическая модель
Т-клетки, присутствующие в иммунной системе человека, представлены ансамблем различных семейств (клонотипов), каждый из которых характеризуется уникальным рецептором, активирующимся при контакте с определенным классом дендритных клеток. Клонотипы, имеющие общее подмножество активаторов, конкурируют таким образом за «стимулы выживания». Получившая «стимул выживания», Т-клетка делится; независимо от этого, клетка может погибнуть в некоторый момент времени. Простейшей стохастической моделью такого поведения ансамбля Т-клеток является марковский процесс рождения-уничтожения [12]. Конкуренция за доступ к
«стимулам выживания» описывается коэффициентами рождения (для одной клетки) следующим образом:
у = £ - = £ , (1)
4 4
где г - индекс клонотипа; Qi - множество «стимулов выживания», активирующих клетки данного клонотипа; щ - число клеток данного клонотипа; у - коэффициент активации; Нщ - суммарное число клеток конкурирующих клонотипов; пЯ - полное число клеток, конкурирующих за доступ к «стимулам выживания» с индексом д. Таким образом, коэффициент рождения тем больше, чем слабее конкуренция за «стимулы выживания». В данной модели предполагается, что вероятность гибели любой Т-клетки в единицу времени не зависит от числа клеток конкретного клонотипа, поэтому коэффициент уничтожения (для одной клетки) ц постоянен. Для г-го клонотипа как популяции коэффициенты рождения и уничтожения равны пУ и п^, соответственно.
В приближении среднего поля коэффициент рождения (1) принимает вид [12]
^ ~-:-> (2)
Н + Vi (п)
где (п) - средний размер клонотипа по ансамблю, а vi - эффективный коэффициент конкурнции данного клонотипа с остальными. Динамика математического ожидания числа Т-клеток в этом клонотипе описывается нелинейным обыкновенным дифференциальным уравнением первого порядка [13]
п=п4—YQ^ \ - м • (3)
\т + Vi (п) )
В зависимости от параметров уравнения, устойчивым может оказаться как нулевое состояние равновесия (отвечает вымиранию клонотипа), так и нетривиальное (выживание).
Недостатки приближения среднего поля очевидны уже из уравнения (3): величина среднего поля нечувствительна к вариациям щ. Более того, условием справедливости данной аппроксимации является малость дисперсии размеров клонотипов и их эффективных коэффициентов конкуренции по отношению к соответствующим средним величинам [12]. Легко видеть, что это условие нарушается, когда, например, вымирает существенная часть клонотипов, либо когда коэффициенты конкуренции сильно варьируют.
Вопрос о грубости (структурной устойчивости) среднеполевого решения по отношению к подобным отклонениям, если они малы, и о характере режимов, если они велики, остается по-прежнему чрезвычайно сложным в рамках исходной многомерной модели. Как указано выше, в настоящей работе ставится более простая задача и ищется ответ на данные вопросы в рамках модели конкуренции двух неидентичных (в том числе, сильно неидентичных) клонотипов на фоне среднеполевого описания влияния остальных клонотипов.
Используя методику, аналогичную предложенной в [12], можно получить следующие выражения для коэффициентов рождения в двумерном марковском процессе [14]:
, У$12 . У$1
м =-:-:—; г +
и + П2 + V (и) П1 + VI (и) (4)
, У$12 + У$2 ()
Л2 — -;—г—|--;—т,
Щ + П2 + V (и) П2 + V2 (и) где $12 - размер подмножества «стимулов выживания», общих для двух клонотипов; $1,2 - размеры подмножеств «стимулов выживания», общих для каждого из двух клонотипов и остального ансамбля; V, V!, V2 - коэффициенты конкуренции между первым и вторым, первым и остальными, вторым и остальными клонотипами, соответственно. Как указано выше, в отличие от среднеполевого описания, рассматривающего динамику только одного клонотипа в явном виде, данное приближение допускает сильную неидентичность между двумя клонотипами как в параметрах, так и в результирующих размерах.
В данной работе исследуется динамика математических ожиданий числа Т-клеток клонотипов. Соответсвующие уравнения имеют вид [13]
• , . У$12 . у$1
и 1 — П1 -Ц1 +------1-1- +
П1 + и2 + V (и) Щ + Vl (и)
• , . У$12 . У$2 и 2 — и2 -Ц2 +--:-:-ГГ +
и1 + и2 + V (и) и2 + V2 (и)
(5)
или после замены переменных £ — ц^, ц — ц2/ць Р1,2 — $12/($12 + $1,2), «1,2 — Ц1,2/(У($12 + $1,2))и1,2, а — Р1 /(ЦР2), V — ^1/(у($12 + $1)) (и), £1,2 — v1)2ц1/(Y($12 + $1)) (и) (знак «тильда» в дальнейшем для краткости опускаем) получаем
• ( 1 , Р1 , 1 - Р1
и 1 — ин -1+-------+
«1 + «2 + V «1 + Vl 7 (6)
• . 1 . ар2 а(1 - Р2)4 « 2 — Ц«2 -1 +-------+
«1 + «2 + V «2 + V2
Здесь все переменные и параметры являются безразмерными и неотрицательными: «1,2 > 0 - число Т-клеток каждого из клонотипов; ц - отношение коэффициента уничтожения для второго клонотипа к первому; ца - имеет смысл аналогичного отношения префакторов для коэффициентов рождения; 0 < р1,2 < 1 - доля одновременно распознаваемых патогенов (р1 — ар2); V, v1, v2 - среднеполевой вклад остальной части популяции Т-клеток в конкуренцию за общие и собственные патогены, соответственно. Обозначим левую часть первого уравнения как Р(«1,^), а второго - как $(п1,п2).
Поскольку система двумерная, квазипериодические и хаотические колебания отсутствуют. Как будет показано в дальнейшем, предельные циклы также не могут реализоваться. Легко видеть, что возможны следующие виды состояний равновесия: (а) тривиальное «1,2 — 0 (существует всегда); (б) частично тривиальные
и1 = 0, и2 = 0 и и1 = 0, и2 = 0; (в) нетривиальные и1,2 = 0. Таким образом, ключевым вопросом являются существование и устойчивость состояний равновесия, определяющие результат процесса конкуренции.
2. Глобальная динамика
Отсутствие предельных циклов вытекает из критерия Бендиксона-Дюлака. В самом деле, при выборе В(и1,и2) = 1/(и1и2) сумма (ВР)'и1 + (BQ)U2 является строго отрицательной для всех и1>2
(ВРи = - - ((-^^ + Т1-^ ) < 0,
1 и2 \(и1 + и2 + V)2 (и1 + Vl)2 )
(RQ)' = _Р2 1 - Р2 \ < 0
1 Ч)и2 и1 \,(и1 + и2 + + (и2 + V2)V
Кроме того, удается сконструировать поглощающую область, которая гарантирует отсутствие траекторий, уходящих на бесконечность и наличие по крайней мере одного устойчивого состояния равновесия внутри.
Правые части уравнения (1) ограничены снизу
1 , Р1 , 1 - Р1 ^ 1,1 - Р1 ^ 0
— 1 +---1--> -1 +--> 0 при и1 < 1 - р1 - v1,
и1 + и2 + V и1 + Vl и1 + Vl
ар2 а(1 — р2) а(1 — р2)
- 1 +-—-+ —-— > -1 + —-— > 0 при и2 < а(1 - р2) - V2•
и1 + и2 + V и2 + V2 и2 + V2
(8)
Таким образом,
{и 1 > 0 при и1 < 1 - р1 - v1, и2 > 0 при и2 < а(1 - р2) - v2• Следовательно, нижние границы поглощающей области задаются уравнениями
( С++ : и1 = шах{1 - р1 - Vl, 0}, (Ю)
[ С+ : и2 = шах{а(1 - Р2) - V2, 0}
Аналогично
Р1 1 _ Р1 Р1 1 _ Р1 1
- 1 + -—- + -— < -1 + —— + -— < -1 + -—7-Г < 0,
и\ + и2 + V П\ + V! П\ + V П\ + V! П\ + шш^, V!}
при П\ > 1 - шш^, V!},
аР2 а(1 — Р2) аР2 а(1 — Р2) а
- 1 +-—-+ —-— < -1 + + —-— < -1 +-р-7 < 0,
и! + П2 + V П2 + V2 П2 + V П2 + V2 «2 + шin{v, V2}
при и2 > а - шin{v, v2 }•
(11)
Отсюда
(и 1 < 0 при и1 > 1 - шт^, v1}, (12)
и2 < 0 при и2 > а - шin{v, v2},
и верхние границы поглощающей области задаются формулами
: u1 = 1 — min{v, vi}, L— : u2 = а — min{v, v2}.
(13)
Рис. 1. Поглощающая область при = 0.5, Р2 = 0.6, V = 0.3, VI = 0.4, VI = 0.2, ц =1, а = 5/6. Стрелки указывают направление потока
Таким образом, поглощающая область ограничена линиями £+'- (рис. 1), и все нетривиальные состояния равновесия находятся внутри нее. Если границы £+2 лежат в первом квадрате и не совпадают с осями координат, в системе существует, по крайней мере, одно устойчивое нетривиальное состояние равновесия.
3. Локальная динамика
3.1. Тривиальное состояние равновесия. Тривиальное состояние равновесия О0,0 : {«1,2 — 0} существует при любых параметрах системы. Его устойчивость определяется характеристическими показателями
, 1,^1 , 1 — Pi ^ п ^ *
А1 = — 1 +---1--< 0 при v1 > v*,
v v1
(P2 1 — P2 \
—1 + а--+ а- < 0 при v2 > v2,
v v2 2
(14)
где v* = v(1 — p1)/(v — p1), v2 = va(1 — p2)/(v — ар2). Устойчивость O°'° (A1)2 < 0) соответствует вымиранию обоих клонотипов. При v < P1 первый показатель заведомо положителен: A4 > 0 V v1, а при v < ар2 заведомо положительным становится второй: A2 > 0 V v2. Указанные ограничения сверху на v означают, что при
не слишком большой конкуренции с остальной популяцией один или оба клонотипа получают достаточно стимулов выживания от специфичных распознаваемых антигенов, и вымирания не происходит. Бифуркационная диаграмма для тривиального состояния равновесия в плоскости параметров {v1, v2} представлена на рис. 2.
Система (1) обладает инвариантными многообразиями M1 : {u1 = 0}
0 0 и M2 : {u2 = 0}. Они глобально устой-Рис. 2. Устойчивость O ' при p1 = 0.5,
Р2 =0.6, v = 0.3, И = 1, а = 5/6. S0'0-устойчи- чивы при P(u1,u2) < 0 V u1,2 > 0 и
вый узел, Sd0'0 - седло, U0,0 - неустойчивый узел Q(u1, u2) < 0 V u1 2 > 0,
V2 5.0
4.0
3.0
2.0
1.0
0
SdQ>°
£/0,0
0 1.0 2.0 3.0 4.0 5.0 Vi
соответственно. Это условие имеет вид
- 1 + * +1-Р1 < 0,
- 1 + 0Р2 + а(1 - Р2) < о,
V У2
что совпадает с условием отрицательности характеристических показателей для О0'0 (14). Когда одно из инвариантных многообразий становится устойчивым, траектории системы асимптотически приближаются к нему по соответствующей переменной. Когда они устойчивы одновременно, значения обеих переменных стремятся к нулю, и других состояний равновесия не существует. Биологической интерпретацией этих режимов является вымирание одного или обоих клонотипов, соответственно.
3.2. Частично тривиальные состояния равновесия. Исследуем существование и устойчивость частично тривиальных состояний равновесия О1'0 : {«1 = и\ =0,42 = 0} и О0'1 : [т =0,П2 = и2 = 0}.
Для О0'1 стационарное значение «2 удовлетворяет уравнению
-1 + -0+- + ^^ = 0, (16)
+ V «2 + V2
откуда следует
(«2)± = 1 (а - V - V2 ± д/(а - V - V2)2 - 4^(1 - а -—— - а—) ) . (17) 2 \ V V2 V у
а(1 - Р2)
Характеристические показатели имеют противоположные знаки, когда 1 -
V2
Р2
-а— < 0 и одинаковые - в обратном случае. Предположим, что а - V - V2 > 0,
V
1 - Р2 Р2 ^ 1 ,л Л + V2 V + V2 тогда а > V + v2 и 1 - а- - а— < 1 - (1 - р2)- - р2- =
v2 V v2 V
V V2
= -(1 - р2)--р2— < 0. Таким образом, (и2)- < 0, и только одно частично
V2 V
тривиальное состояние равновесия О0'1 : {«1 = 0,42 = (и£)+} может существовать.
Возьмем теперь (и2)+. Как показано выше, а - V - V2 становится отрицательным
, Ч 1 1 - Р2 Р2
(например, при увеличении V2) прежде чем 1 - а--а— становится положи-
V2 V
тельным. Это означает, что (и^)+ > 0 пока 1 - а-— - а— < 0. Следовательно,
V2 V
^ - а - а^ = ^ (18)
задает бифуркационное многообразие в пространстве параметров, соответствующее влипанию О0'1 в О0'0 (транскритическая бифуркация). По одну сторону существует единственное частично тривиальное состояние равновесия О0'1, по другую - ни
одного. Обратим внимание, что эта граница точно совпадает с изменением устойча-вости многообразия и знака А2 у О0'0 (14), (15). Аналогично, для О1'0 имеем
(«п+ — 2(1 - V - Vl + у (а - V - Vl)2 - 4^1(1 - - р1) ) (19)
и бифуркационное многообразие
РеХ - 1--р - р1 — ^ • (20)
Заметим, что О0,0 становится устойчивым после последовательного слияния с О1,0 и О0,1.
Исследуем теперь устойчивость частично тривиальных состояний равновесия. Характеристические показатели О0,1 задаются диагональными коэффициентами матрицы линеаризации
— Р'и1 — -1 + -Р- + ,
и1 «2 + V v1
\ ПГ (-11 P2V , (1 - Р2>2\ . / 1 . Р2 . 1 - Р2 \ 0
^2 — $и2 — ц -1 + а Л2 + а; , \2 <ц -1+а^--^а^- — 0.
и2 V («2+v)2 («2 + V2)2У V «2 + V «2 + V2 )
(21)
Легко видеть, что А, 2 — 0 в момент слияния с О0,0: «2 — 0. Поскольку А,1 —
1 . Р1 ,1 - Р1 ^ 1 . Р1 ,1 - Р1 — -1 +-----1--< -1 +---1--, его значение отрицательно в мо-
«2 + V V1 V v1
мент исчезновения О1,0 (20). Таким образом, О0,1 устойчиво при условии, что О1,0 не существует. Однако это условие лишь достаточное: бифуркация, приводящая к приобретению устойчивости, происходит до исчезновения О1,0. Решая систему уравнений
$(0,«2) — -1 + -От- + ^^ — 0,
«2 + V 4*2 + V2
Х1 — -1 + -Р+- + — 0.
«2 + V v1
можно получить точное выражение, отвечающее бифуркационному многообразию
Р^1 ( Р1 Vl
а —
-0,1 I . Vl - 1+ Р1 \ Vl - 1+ Р17 .
Рз1 : < ^ — v + -Р^-- } . (23)
аР2
Vl - 1 + Р1
Устойчивость О1,0 определяется следующими характеристическими показате-
лями:
I — Р/ — 1 + Plv + < _1 + Р1 + 1 - Р1 — 0
1 и1 («1 + V)2 («1 + Vl)2 - ^ + V + Vl
1 /О' ( 1 I Р2 ,1 - Р2 \ ^ ( 1 , Р2 ,1 - Р2 \ А2 — $и2 — ц -1 + а^—■--+ а- < ц -1 + а--+ а- .
и2 \ «I + V V2 У V V V2 /
Аналогично предыдущему случаю, О1'0 устойчиво, если О0'1 не существует. Бифуркационное многообразие приобретения устойчивости задается выражением
а^2
V
1,0
VI = V +
У2 - а(1 - Р2)
1 -
V2 - а(1 - Р2),
аp2V2
V2 - а(1 - Р2)
Р1
(25)
Итак, О1''0 и О0'1 меняют устойчивость в области своего существования (см. пример бифуркационной диаграммы на рис. 3). В частности, они становятся устойчивыми при увеличении V2 и V!, соответственно (при дальнейшем увеличении происходит влипание в тривиальное состояние равновесия). Особо следует обратить внимане на тот факт, что отсутствие одного из состояний гарантирует устойчивость другого (в области существования). Биологическая интерпретация данного эффекта такова: когда один из клонотипов испытывает слишком большую конкуренцию, он вымирает, и выживает только другой.
Обратимся наконец к вопросу о существовании и устойчивости нетривиальных состояний равновесия О1'1 : {и^2 = 0} системы (1). Они удовлетворяют системе уравнений
£0,0
■¿од / гех Ул) 1,0,* и 0,0 [ '-¿я- 50-1; -
2
5.0 4.0 3.0 2.0 1.0 о
0 1.0 2.0 3.0 4.0 5.0 VI
Рис. 3. Устойчивость частично тривиальных состояний равновесия в плоскости параметров {VI,V2} при р1 = 0.5, р2 = 0.6, V = 0.3, = 1, а = 5/6. 5и 50,1 - устойчивые узлы, в^'0 и в^'0 - седла, и1,0 и и0,1 - неустойчивые узлы
- 1 +
Р1
- 1 +
П1 + П2 + V аР2
+
+
1 - Р1 =
П1 + Vl
а(1 - Р2)
0,
(26)
0.
П1 + П2 + V П2 + V2
Решения отвечают точкам пересечения кривых и2 = (и1) и и1 = Т2(и2), задаваемых соответственно первым и вторым уравнениями,
и2 = ^Ы = -и1 - V +
и1 = ?2(и2) = -и2 + Р2а - V +
Р1(и1 + Vl) -1 + Р1 + и1 + Vl '
(1 - Р2)Р2а2
(27)
(-1 + Р2)а + и2 + V2
Функции ^1'2(и1'2) являются непрерывными внутри поглощающей области, то есть над прямыми £+2 (6). Поскольку все нетривиальные состояния равновесия принадлежат этой области, ограничимся исследованием этих функций в ней. Пусть и\ 2 являются решением системы
\и*2 = л (и )
\и*1 = ?2(и*2).
Тогда Щ является решением уравнения
^(щ) = ^(^(щ)) - щ =0, (29)
и [^"21^1)]^ = [^2(п2)]^3 [^(щ)]^ - 1. Подстановка (27) дает
гт ( )]/ 1 (1 - Р1)Р1 / 1
1^1(п1)1и1 = -1 — —--Т77 < -1,
1 и 1Ли1 (-1+ Р1 + П1 + У1)2
( )]/ 1 (1 - Р2)Р2а2 1
[•^М] = -1 - ---г-■-■-г^ < -1.
"2 (-(1 - Р2)а + П2 + У2)2
(30)
Таким образом, [^21^1)]^ > 0 и (29) имеет не более одного решения. Аналогичные рассуждения могут быть проведены по отношению к Г12(п2) = Т1(Т2(и2)) -п2 = 0.
Итак, система (1) не может иметь более одного нетривиального состояния равновесия. Если £+2 лежат в первом квадранте и не совпадают с осями координат (то есть 1 - р1 - у1 > 0, а(1 - р2) - у2 > 0), то О1'1 существует, единственно и устойчиво. В противном случае, ни существование, ни стабильность не гарантированы (устойчивыми могут быть частично тривиальное или тривиальное состояния равновесия). Наконец, когда инвариантные многообразия М122 становятся глобально устойчивыми (28) (одновременно с приобретением устойчивости О0'0), никаких иных состояний равновесия существовать не может.
Какая бифуркация приводит к устойчивости О0'1? (Пока оно неустойчиво, О1'1 обязательно существует и устойчиво. Здесь мы предполагаем, что О1'0 либо неустойчиво, либо не существует, поскольку устойчивость О1'0 и О0'1 взаимоисключающи. Аналогичный анализ можно провести в предположении устойчивости О1'0.)
Покажем, что О1'1 сливается с О0'1 как раз в тот момент, когда О0'1 приобретает устойчивость (23). В самом деле, предполагая п1 = 0 в (26)
Р1 1 - Р1
- 1 +-—-+-— = 0,
П1 + П2 + V Щ + V!
_1+ аР2 + а(1 - Р2) =0 (31)
П1 + П2 + V П2 + V2 '
получаем
п1 = 0,
Р1 1 - Р1
- 1 + + -— = 0,
П2 + V Vl
_ 1 + ар2 + а(1 - Р2) _
(32)
П2 + V П2 + V2
что идентично (22). Таким образом, О0'1 становится устойчивым в результате транскритической бифуркации с устойчивым узлом О1,1.
Убедимся, что после бифуркации формальное решение О1'1 имеет п1 < 0 и не существует в области неотрицательных размеров клонотипов. Вследствие большого числа параметров, переход между (1 - p1)/v1 = 1 (влипание в ось у) и (1 - Р1)^1 = 1 - Рl/v (М1 становится устойчивым) можно осуществить бесконечным числом путей (как и в ситуации с £+ и М2). Тем не менее, каждый из них содержит транскритическую бифуркацию О1'1 и О0'1. Будем, для определенности,
0
Рис. 4. Устойчивые состояния равновесия в плоскости параметров {VI, у2} при р1 = 0.5, р2 = 0.6, V = 0.3, = 1, а = 5/6. Вставки: примеры локальных фазовых портретов, по горизонтали - и1, по вертикали - и2
увеличивать VI при постоянных значениях остальных параметров. Тогда и* является функцией переменной VI.
Легко показать, что ду1 и* < 0 и, следовательно, и\ < 0 и О1'1 не существует после бифуркации. Координаты состояния равновесия удовлетворяют уравнению (29)
Г21(и*1, VI) = ЛОЛК, VI)) - и*, =0. (33)
Частные производные по V! равны
дv1 [^2(Л(и1, Vl)) - и1] = ди*?2 [ди*1 Flдv1 и1 + дv1 Л] - дv1 и1 =
= дvl и1 [ди1 ?2ди1 - 1 + ди1 ?2 дvl = 0.
(34)
Отсюда
„ * ди* ?2 дv*
дv1 и1 = 2 1
1 - ди*2 ди*1
(35)
Из (30) и ?1 = -^1(1 - Р1)/(и1 + Vl + Р1 - 1)2 < 0 следует 1 - ди* К ди* К < 0, ди* Т2 д^* > 0 и, окончательно, дV1 и1 < 0, что завершает доказательство.
В результате, типичное разбиение плоскости параметров Vl'2 на области устойчивости различных состояний равновесия имеет вид, представленный на рис. 4.
Заключение
В работе исследована динамика двухкомпонентной модели популяции Т-клеток. Показано, что периодические колебания в системе невозможны, и траектории асимптотически приближаются к состоянию равновесия. Определены области существования и устойчивости различных состояний равновесия. Установлено, что для всех
значений параметров лишь одно из них является устойчивым. В зависимости от конкретных значений, определяющих конкурентное преимущество одного из кло-нотипов над другим, или остальной части популяции (описываемой в среднеполе-вом приближении) над обоими, возможны: (1) сосуществование двух клонотипов, (2) вымирание одного из них, (3) вымирание обоих. Таким образом, среднеполевое приближение неприменимо при больших вариациях характеристик клонотипов (режимы вымирания одного или обоих), однако обладает свойством грубости по крайней мере относительно небольших вариаций характеристик второго кло-нотипа, выделенного из среднеполевого описания (режим сосуществования).
Полученные результаты дают представление о динамике многомерной стохастической модели в частности и регуляции численности и состава в популяции иммунных Т-клеток в целом: отсутствие периодических колебаний моментов случайных величин - размеров клонотипов, единственность и устойчивость (квази)ста-ционарного распределения вероятности размеров, границы областей вымирания отдельных клонотипов. Они поддерживают гипотезу о природе механизма регуляции численности и видового разнообразия Т-клеток, воспроизводят экспериментально наблюдаемые отсутствие периодических колебаний в составе ансамбля Т-клеток и быстрое возвращение к прежнему равновесному состоянию в постинфекционный период.
Представленный анализ, будучи ограниченным случаем двух клонотипов, закладывает, вместе с тем, основу для развития теории конкуренции в многокомпонентных моделях иммунной системы. Предварительные аналитические и численные результаты для трех- и N-компонентных моделей (последние уже не используют приближение среднего поля) указывают на то, что автоколебания не возникают и в таких системах, а состояние равновесия (будь то нетривиальное, тривиальное или одно из 2n — 2 частично тривиальных состояний) всегда существует и, по-видимому, единственно.
Автор благодарит К. Молина-Парис (C. Molina-Paris), Г. Лайт (G. Lythe) и Э. Стирк (E. Stirk) за плодотворные дискуссии, а также рецензента за конструктивные замечания и предложения, которые позволили существенно улучшить статью.
Библиографический список
1. Arecchi F.T. et al. Pattern formation and competition in nonlinear optics II Physics Reports. 1999. Vol. 318. P. 1.
2. Afraimovich VS., Nekorkin V.I., Osipov G.V., Shalfeev V.D. Stability, Structures and Chaos in Nonlinear Synchronization Networks. Singapore: World Scientific, 1994.
3. Kuramoto Y Chemical oscillations and turbulence. Tokio: Springer, 1984.
4. Короновский А.А., Трубецков Д.И. Нелинейная динамика в действии. Саратов: Изд-во ГосУНЦ «Колледж», 2002.
5. Murray J.D. Mathematical Biology. Berlin Heidelberg: Springer-Verlag, 2002.
6. Rabinovich M.I., Varona P., Selverston A.I. et al. Dynamical principles in neuroscience II Rev. Mod. Phys. 2006. Vol. 78. P. 1213.
7. Nekorkin V.I., Shapin D.S., Dmitrichev A.S. et al. Heteroclinic contours and self-replicated solitary waves in a reaction-diffUsion lattice with complex threshold excitation // Physica D. 2008. Vol. 237. P. 2463.
8. Jameson S.C. Maintainting the norm: T-cell homeostasis // Nature Reviews Immunology. 2002. Vol. 2. P. 547.
9. Perelson A. and Weisbuch G. Immunology for physicists // Reviews of Modern Physics. 1997. Vol. 69. P. 1219.
10. Freitas A.A. and Rocha B.B. Lymphocyte lifespans: homeostasis, selection and competition // Immun. Today. 1993. Vol. 14. P. 25.
11. de Boer R. and Perelson A. Competitive control of the self-renewing T cell repertoire // International Immunology. 1997. Vol. 9. P. 779.
12. Stirk E.R., Molina-Paris C., and van den Berg H.Stochastic niche structure and diversity maintenance in the T cell repertoire // J. Theor. Biol. 2008. Vol. 255. P. 237.
13. Gillespie D.T. Stochastic simulation of chemical kinetics // Annu. Rev. Phys. Chem. 2007. Vol. 58. P. 35.
14. Stirk E., Lythe G., van den Berg H., Hurst G. and Molina-Paris C. The limiting conditional distribution in a stochastic model of T cell repertoire maintenance // Math. Biosciences, 2010 (in print).
University of Leeds Поступила в редакцию 3.12.2009
United Kingdom После доработки 24.02.2010
COMPETITION IN THE TWO-COMPONENT MODEL OF THE IMMUNE T-CELL ENSEMBLE
M.V. Ivanchenko
We study the process of competition in the two-component model of the immune T-cells ensemble that underpins the selection mechanism of the most efficient T-cell species (clonotypes). We demonstrate the absence of periodic oscillations, determine the regions of coexistence, partial and mutual extinction of clonotypes. Applicability of the mean field approximation is analyzed. The biological implications of the results are discussed.
Keywords: Competition, bifurcations, immune system.
Иванченко Михаил Васильевич - родился в 1981 году. Окончил радиофизический факультет Нижегородского государственного университета им. Н.И. Лобачевского (2004), там же защитил кандидатскую диссертацию (2007). В настоящее время является сотрудником кафедры прикладной математики универ-систета Лидса (Великобритания). Область научных интересов - нелинейная динамика, математическая иммунология, локализация энергии в нелинейных решеточных системах, синхронизация, нейродинамика. Имеет около 35 печатных работ по данной тематике. Был удостоен медали Российской академии наук в конкурсе 2004 года на лучшие научные работы по направлению «Общая физика и астрономия».
Leeds, United Kingdom, LS2 9JT
University of Leeds, Department of Applied Mathematics E-mail: [email protected]
