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

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

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

Здесь использованы те же обозначения, что и в системах (6) – (9) (индексы опущены): С₁ – множество клеток разного возраста до митоза, С_m – митотические клетки; C_a – клетки в анабиозе; (C_k X_l ) – ингибированные клетки разных стадий; M₁ , M₂ –субстраты.

Следует отметить, что двухстадийный цикл (фазы S и M) наблюдается на ранних стадиях развития зародышей пойкилотермных животных [15, 16]. На этом основании в этот период в качестве единичного интервала времени можно использовать длительность t_c клеточного цикла («детлаф»).

В предположении постоянства концентраций субстратов М₁ , М₂ кинетика цепного роста популяции, состоящей из особей С₁ и С_m , описывается системой:

Здесь c₁ , c_m – количества растущих и митотические клеток; a, b, p –коэффициенты автоингибирования, рождения и роста популяции в отсутствии ингибиторов. В коэффициенты р и b включены постоянные количества субстратов М₁ и М₂ . f - коэффициент размножения. Коэффициенты b_x и p_x – функции количества ингибиторов x₁ и x₂ :

Система уравнений (12) представляет собой закон обобщенного движения двухстадийной популяции в пространстве состояний. Преобразованная в виде:

система (12) определяет соотношение между интервалами биологического dc_i и физического dt времени.

В приближении квазистационарности для митотических клеток С_m система (12.1-12.2) сводится к одному уравнению:

Динамику численности популяции с₁ (t) в общем случае нельзя выразить в виде явной функции от времени t. Поэтому используют обратную функцию t(c₁ ), получаемую интегрированием (14) по с₁ :

Уравнение (16) в явном виде отображает физическое время на множество состояний популяции.

6. Интерпретация модели (I-компонент теории, interpretation).

На рис.1 приведены экспериментальные точки и графики функции (16), описывающие дрожжевых клеток в присутствии солей хрома и никеля [13]. При расчете графиков брали значения a, b, р , f, определенные по экспериментальным данным [17]. В пределах точности измерений расчетные кривые согласуются с экспериментом при измененении численности примерно на шесть порядков.

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

Рис.1. Экспериментальные точки и графики функции (16), описывающие рост пивных дрожжей при разных концентрациях (ммоль/л) солей хрома и никеля [13]: (1) c(Ni) = c(Cr) = 0.0, (2) c(Ni) = 0.5, (3) c(Cr) = 0.5, and (4) c(Ni) + c(Cr) = 0.5 + 0.5. Коэффициенты: a=1.25.10^-7 мл /ч, b=0.8ч^-1 , р=0.32 ч^-1 , f =2.

Уравнения (12a) можно представить в виде:

где величины K_c1 =(-p_x c₁ +f b c_m +w₁ )и K_cm = (p c₁ -b_x c_m -a c₁ c_m )представляют собой калибровочные коэффициенты для перехода от интервала физического времени dt к интервалам биологического времени dc_j .

На основе (17) получают соотношение между конечными временными интервалами:

D c1 =c1 dt, D cm =ИНТЕГРАЛ( Kcm dt ).

( 17а)

Калибровочные соотношения (17) обладают следующими свойствами:

1. Коэффициенты K_c1 и K_cm зависят от кинетических констант, характеризующих внутри- и внесистемные взаимодействия. Это определяет специфику биологического времени данного объекта.

2. Коэффициенты K_c1 и K_cm зависят от наблюдаемого состояния объекта, то есть изменяются при движении по фазовой траектории.

3. Коэффициенты K_c1 и K_cm неодинаковы для однотипных элементов данного уровня иерархии. Это означает, что собственное время «течет» с разной скоростью не только на разных уровнях биологической системы, но и в различных элементах одного уровня.

Приращение суммарной массы dm_p или численности dN_p популяции определяют интервал биологического времени популяции в целом. Для двухстадийной популяции dN_p = V(с₁ +с_m ), где V – объем системы. Связь между популяционным и физическим временем согласно (17) определяется соотношением:

Через длительность клеточного цикла t_c в физической шкале (в «детлафах») эта величина выразится в виде:

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

Приращение численности популяции D с₁₂ =c₂ -c₁ в единице обьема наблюдается за время D t₁₂ =t₂ -t₁ , согласно (16) равное:

где c₁ ,c₂ – численности в моменты t₁ ,t₂ .

Среднее число делений n₁₂ каждой из c₁ клеток за это время равно:

Следователъно, t _c можно оценить по формуле: