Биологическое время и его моделирование в квазихимическом пространстве

Методология построения теории времени естественных объектов, детально изложена [1, 2]. В данной работе рассмотрены компоненты этой теории на примере клеточной популяции.

1. Структура клеточной популяции. (O -компонент теории, object).

Первая проблема состоит в выборе составных частей (элементов, таксонов, униструктур) объекта и определении их пространственных границ [3, 4].

Принципиальная трудность состоит в том, что строгий критерий членения объекта на части отсутствует. Выделение объектов природы – идеализация, связанная, главным образом, с особенностями информационных взаимодействий живых организмов со средой [5].

Объективным критерием членения объекта может служить принцип энергетической дифференцировки [6], который можно записать в виде:

Еj+1j,

где Еj , Еj+1 — прочность связи между элементами уровня j и между элементами более высокого уровня j+1 (схема 1).

Для клеточной популяции порядок элементов с усложнением уровня (цифры - номера j уровней) можно представить в виде:

1 молекулы(t m)I 2 органеллы(t o) I 3 клетки (t c) I 4 популяция(t p) I 5 экосистема(t e) (1)

Каждый уровень j разбивается на подмножества элементов (схема 1).

Схема 1. Иерархия подсистем естественных объектов

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

Ms = (Ms1, Ms2, ... Mse) Xs = (Xs1, X2s, ... Xse) (2) (3)

где Ms, Xs - векторы (наборы) субстратов (Ms1, Mse) и химических агентов (Xs1,Xse) для биологического вида S.

Многомерность структуры (1-3) создает слишком большие сложности не только для количественной, но и для качественной идентификации модели. Поэтому целесообразно исследовать агрегированные модели. В случае клеточной монокультуры весьма информативной оказывается структура (1) из трех уровней:

(3) клетки I (4) популяция I (5) экосистема (4)

Описательная модель клеточного цикла [7], начиная с митоза Cm, может быть представлена в виде разветвленной цепи:

(5a)

или в более сжатом варианте:

Cm ® f C1® C2® C3® Cm

(5b)

Здесь C1 - молодая клетка непосредственно после деления, f - коэффициент размножения, который в общем случае может отличаться от 2. С2, С3, Сm - фазы развития каждой из f клеток C1.

Фазы C1, C2, C3, Cm соответствуют стадиям клеточного цикла G1, S, G2, M (клетки четырех возрастов) и являются подмножествами 4-го уровня иерархии (4) - популяции. Популяция с наборами химических агентов (2) и (3) формирует 5-й уровень модели (4).

2. Пространство допустимых состояний популяции (S-компонент теории, states).

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

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

В модели (4) состояние задается числами всех типов клеток, химических и физических агентов среды. Уменьшение размерности модели достигается путем аггрегирования разных элементов одного уровня по тому или иному признаку.

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

(5c)

3. Изменение состояний популяции, «предвремя» (C-компонент теории, change).

Разбиение множества особей на возрастные подмножества зависит от критериев, заложенных в основу членения объекта по данному признаку, от разрешающей способности экспериментальной методики и других факторов. Это обстоятельство хорошо иллюстрируется общепринятым членением клеточной популяции на четыре стадии (фазы) развития клеток C1, С2, С3, Сm (G1, S, G2, M). Мало вероятно, что в природе существуют хотя бы две клетки, идентичных по наборам элементов всех уровней. Поэтому отнесение данной клетки к той или иной стадии проводится по степени соответствия индивидуальных характеристик этой клетки и усредненных характеристик стадий.

В результате метаболизма данная клетка практически непрерывно меняет набор элементов (1) и движется по своей индивидуальной траектории в определенном выше пространстве состояний (S-компонент теории). Эта траектория описывает все стадии развития данной клетки (движение в фазовом пространстве) от деления до деления. Последовательность точек этой фазовой траектории, отображенная на числовую ось, определяет предвремя индивидуальной клетки.

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

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

4. Часы и параметрическое время популяции (T-компонент теории, time).

Проблема выбора временных границ существования объекта тесно связана с C-компонентом теории, с определением масштабов предвремени, границ в пространстве.

Событие (момент времени) – переход объекта из одного состояния в другое. Тип события определяется уровнем, на котором оно происходит [2].

Предвремя природного объекта – отношение порядка (refl+trans) на множестве событий в этом объекте. Отображением этого порядка на числовой ряд получают числовое выражение предвремени.

Время – интервал между двумя наблюдаемыми событиями одного типа. Измеряется время числом событий между этими двумя событиями.

С биологической точки зрения возраст особи соответствует некоторой стадии индивидуального развития, начиная от зарождения. Стадия развития биологического объекта – синоним его возраста или биологического времени В общем случае возраст однозначно не связан с физическим временем. Исходя из этого и вводят понятие биологического времени.

Экспликацию понятия “биологическое время” целесообразно дать на примере клеточных популяций. Возраст популяции tp определяется средним числом делений Nd первичных клеток посева или, что то же самое, числом циклов в цепи (5): tp порядка tcNd, где t c - среднее время клеточного цикла.

Каждую из биохимических реакций можно трактовать как замену одних молекул другими на молекулярном уровне структурной иерархии (1). Такая замена осуществляется в результате атомно-молекулярных взаимодействий. Поэтому естественным эталоном времени на молекулярном уровне являются атомные часы. Соответственно временные интервалы между «химическими» событиями измеряются в шкале физического времени. Характерное время tm молекулярного уровня определяется скоростями ферментативных реакций. Интервалы между событиями на молекулярном уровне имеют значения порядка секунд.

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

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

Клетки формируются в результате образования наборов органелл и молекул и их трансформации. Временные интервалы клеток tc имеют значения порядка 104 с.

Популяция развивается в результате изменения наборов клеток в разных фазах развития. Взаимодействия между клетками носят диффузионно-обменный и контактный характер. Временные интервалы популяции tp порядка 105 с.

Динамика экосистемы зависит от запаса ресурсов, емкости, открытости а также от характера взаимодействия популяции со средой. Временные интервалы экосистемы te имеют значения порядка tp.

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

5. Закон обобщенного движения популяции (L-компонент теории, law).

Для вывода закона обобщенного движения популяции в пространстве состояний предлагается использовать единый «квазихимический» язык для описания событий на всех уровнях популяции. При этом все взаимодействия отображаются квазихимическими уравнениями [5, 8, 9]. Пространство состояний при таком описании можно назвать квазихимическим.

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

Популяция развивается в результате усвоения субстратов. Это развитие описывается системой квазихимических реакций (индекс s опускается):

C1+(M1,Me) ® C2 (p11),
C2+(M1,Me) ® C3 (p21)
...........................................

(6)

Cn+(M1,Me)® fC1 (pne(b)) ,

Набор кинетических констант (p11, pne) определяет кинетический вектор роста.

Действие токсикантов Хi, может проявляться на любой стадии роста и описывается сходным образом:

C1+X1® (C1X1) (d11)
C2+X2® (C2X1) (d12)
.........................................

(7)

Cn+Xt® (CnXt) (dnt)
X1,Xt® EE (rx1,rxt)

Здесь (CkXl) – дезактивированные клетки разных стадий, (X1,Xt) – вектор токсикантов; {dij} – матрица коэффициентов ингибирования, (rx1, rxt) - вектор скоростей притока токсикантов из среды.

Кроме того, следует учесть клеточные взаимодействия (автоингибирование). Например, зрелые особи Ci переводят в неактивное состояние молодых Cak:

Ci+C1® Ca1+Ci (ai1)
Ci+C2® Ca2+Ci (ai2)
.......................................

(8)

Ci+Cm® Cam+Ci (aim)

Взаимодействия клеток описываются кинетическим вектором автоингибировния (ai1,…, aim).

В открытых системах учитывается взаимодействие экосистемы со средой ЕЕ:

EЕ C1 w1
EЕ C2 w2

--------

(9)

EЕ Cn wn

где (w1,…, wn) - вектор скоростей притока особей С1,…, Сn из среды.

Система псевдохимических реакций (6)-(9) описывается стандартным образом [10-12] системой кинетических уравнений, представляющих собой закон обобщенного движения популяции в пространстве состояний:

dC/dt = S KC - S ACC - S DCX + S R + S W

(10)

Здесь C = (c1, c2, c3, cm) вектор количества клеток в разных фазах, K, A, D, R, W – матрицы кинетических параметров.

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

Рассмотрим редуцированную модель (1) из двух стадий – роста и деления, дополненных стадией самоингибирования:

C1+M1® Cm (p)
Cm+M2® fC1 (b)
C1® Cd (g)
C1+Cm® Ca+C1 (a) (11)
C1« EE (w1)
C1+X1« (C1X11) (d11)
C1+X2« (C1X12) (d12)
Cm+X1« (CmX21) (d21)
Cm+X2« (CmX22) (d22)

Здесь использованы те же обозначения, что и в системах (6) – (9) (индексы опущены): С1 – множество клеток