Коэволюция фаг-бактерия: проточная модель 1 ресурс 1 пара фаг-бактерия
Популяционная, ресурсзависимая динамика одной пары взаимно специфичных штаммов фага и бактерии в проточном резервуаре может быть описана следующей сисемой дифференциальных уравнений:
* R'=Ro*(Rc-R)-V(R)(X+Z). * X'=V(R)*X-aXY-Ro*X. * Y'=wbZ-aXY-Ro*Y. * Z'=aXY-bZ-Ro*Z (12.4.1) где Rc- концентрация ресурса в резервуаре, а R- его концентрация в поступающем потоке; Ro -скорость протока; V(R)- функция, характеризующая ресурсзависимую скорость размножения бактерий, имеет ограничение сверху ln2 и положительную поизводную по ресурсу (предполагаемый качественный вид этой функции изображен на рис.12.4.1); X- концентрация свободных бактерий; Z- концентрация комплекса фаг-бактерия; Y- концентрация свободных вирулентных фагов; a- константа скорости адсорбции фага на бактерии; b -константа скорости лизиса зараженных бактерий; w- констаната скорости размножения фага (w=ln(K), где K- средний "выход" фага). Качественное исследование системы (12.4.1) показывает, что она имеет три особые точки:
a) (Rc,0,0,0)-тривиальная,
б) (V -1(Ro), R-V -1(Ro), 0, 0)-полутривиальная,
в) (R^, X^, Y^, Z^)-нетривиальная,
Ro*(Ro+b) где X^= ****************, a*[b*(w-1)-Ro]
V(R)-Ro Y^= *********, a
Ro*[V(R)-Ro] Z^= *****************, a*[b*(w-1)-Ro] а R^- корень уравнения
V(R)*(V(R)+b) Rc-R = **************** (12.4.3) a*(b*(w-1)-Ro)
Такой корень существует и единственен при b*(w-1)>Ro, как видно из геометрического представления уравнения (12.4.3) (рис.12.4.2).
Исследуем эти равновесные состояния на устойчивость.
1.Характеристический определитель точки (а) равен
(l+Ro) *(l+b+Ro)2*(l+(Ro-V(Rc)),
т.е. точка (а) устойчива при Ro E V(Rc) и неустойчива при Ro e V(Rc).
2. Точка (б) существует при Ro < V(R), соответствующий характеристический определитель равен
[l 2+ l*(Ro+X^*V'R (R)*(1+Ro)]*[l2 +l*(2*Ro+b+aX^) + (Ro2 +Ro*b-a*X^*(b*(w-1)-Ro))],
т.е. точка (б) устойчива при b(W-1)<Ro и неустойчива при
(Ro+b)*Ro < a*X^*[b*(w-1) - Ro].
Cледовательно, существует интервал
0 < a*X^*[b*(w-1)-Ro] < Ro*(Ro+b),
когда уже существует нетривиальная точка покоя (в), но полутривиальная еще устойчива. Oднако для существования устойчивого полутривиального равновесия должно выполняться условие
b*(w-1) < V(Rс), (12.4.4)
что маловероятно. Действительно, V(R) не может превышать величины ln2 (см. рис.12.4.1), тогда как b~0.2-0.3 и w=ln(K), где K - "выход фага", равный обычно десяткам и даже сотням дочерних фаговых частиц, образующихся в результате лизиса бактерии ( Адамс, 1961 ; Стент, 1965 ). Другими словами, точка (б) неустойчива по фагу и вряд ли вообще существует при биологически реалистичных значениях параметров размножения фага и бактерии.
3. Для существования нетривиальной точки равновесия (в) необходимо, чтобы V(R^) > Ro и b(w-1) > Ro. Характеристический определитель этой точки имеет вид
l 4+ A*l3 + B*l2 + C*l + D, где
A=3*Ro+b+a*X^+(X^+Z^)*V'R (R^),
B=V(R^)*V 'R(R^)*X^+(2*Ro+b+a*X^)*[Ro+(X^+Z^)*V'R (R^)]-a2*X*Y,
C=(2*Ro+b+a*X^)*V(R^)*V'R (R^)*X^ + (12.4.5)
+a*X*Y*[a*(b*(w-1)-2*Ro-(X^+Z^)*V 'R(R^)) + V(R^)*V'R (R^)],
D=a*X*Y*[a*(b*(w-1)-Ro)*(Ro+(X^+Z^)*V'R (Rv))+Ro*V(R^)*V'R(R^)].
Коэффициенты A и D положительны при любых значениях параметров, , функции V(R) и ее производной по ресурсу V'R (R), но B и C могут быть и отрицательны. B самом деле, функция V(R) задана так , (рис.12.4.1), что, например, при R-> бесконечности, V 'R(R)->0.
Рассмотрим такие большие R, при которых величиной V'R (R) можно пренебречь. Тогда
C~a2 *X^*Y^*[b*(w-1)-2*Ro] < 0
b*(w-1) при Ro > ********* 2и точка (а) будет в этом случае неустойчивой. Однако, если подставить в (12.4.5) значения X^, Y^, Z^ из (12.4.3), то нетрудно убедиться, что коэффициенты B и C будут положительны при выполнении неравенств
b*(w-1) > V(R) > 2*Ro. (12.4.6),
Поскольку явный вид функции V(R) и V'R(R) неизвестен, более сильные, нежели (12.4.6), достаточные условия устойчивости нетривиального равновесия получить не удается; по этой же причине не исследовался вопрос о существовании устойчивых циклов в данной модели. Тем не менее очевидно, что при биологически реалистичных (и оправданных самой постановкой задачи значениях параметров системы) наиболее вероятным устойчивым состоянием является нетривиальное равновесие с сохранением и фага, и бактерии, т.е. с невырождением системы как таковой.