Временная иерархия в биологических системах 1 Масштабы времен

Временная иерархия в биологических системах 1. Масштабы времен при моделировании 2. Особенности биологических систем 3. Быстрая и медленная динамика 4. Метод квазистационарных концентраций 5. Теорема Тихонова 6. Обезразмеривание 7. Редукция Михаил Пантелеев 22 октября 2011

Задача 1 А l=100 km t=1 h B

Задача 2 v=100 km/h a=1 m/s 2 B

Задача 3 P=5%

Усложнение задачи: Биологические, социальные, экономические, экологические системы. Какие характерные времена тут?

А тут?

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

Быстрая и медленная динамика A B

Метод квазистационарных концентраций (Он же метод Боденштейна-Семенова, он же метод квазистационарных приближений) 1. Если в сложной реакции времена жизни промежуточных продуктов малы, то можно считать процесс стационарным и заменить дифференциальные уравнения алгебраическими. 2. Система дифференциальных уравнений, получаемых в результате применения закона действующих масс к элементарным стадиям, сводится к системе уравнений алгебраических, поскольку все производные концентраций промежуточных веществ по времени полагают равными нулю. 3. Обратите внимание на некоторую парадоксальность приема: формально мы приравниваем нулю скорости изменения САМЫХ БЫСТРО МЕНЯЮЩИХСЯ веществ.

Проблема размерности A С V=1 n. M/min V=1000 n. M/min B D

Проблемы с квазистационарностью 1. МКК удобен в простом химическом реакторе с цепочкой превращений и нормальной динамикой. В биохимических системах: а) много входов и выходов, б) каждое вещество участвует в куче сложных реакций, в) динамика может быть и колебательной, и хаотической, и любой непонятно заранее какой. 2. Тем более, что там есть масса переменных, не являющихся промежуточными продуктами. 3. И вообще полно процессов, не обязательно являющихся химией. 4. И не факт, что есть только два масштаба - быстрый и медленный. 5. Даже в простом реакторе применимость МКК в каждом конкретном случае нереально надежно оценить. 6. Это не мешает им пользоваться в биохимии, иногда даже успешно.

Регулярное и сингулярное возмущение

Теорема Тихонова (плюс расширение Васильевой) Полная система: Вырожденная система: 1. Пусть полная система (где u и v могут быть векторами; плюс начальные условия) имеет достаточно гладкие правые части и вообще имеет решение 2. Пусть вырожденная система тоже имеет решение 3. Если присоединенная система имеет непрерывное решение, которое является изолированным 4. Если начальные условия полной системы находятся в области притяжения этого решения вырожденной системы 5. То решение полной системы стремится к решению вырожденной системы при ε→ 0 (однако, для v это справедливо только при t > d > 0) 6. Причем толщина (масштаб) пограничного слоя имеет порядок ε

Обезразмеривание Процедура нормировки концентраций, времени, координат на некоторые "естественные" масштабы, в результате которой все переменные будут иметь близкий порядок. A+B C Concentrations at t=0: A=1 p. M B=100 p. M C=0 Velocities: k+1=10 n. M-1 min-1 k-1=1 min-1

Анализ временной иерархии и редукция 1. Написать систему уравнений. 2. Произвести обезразмеривание концентраций и времен на интересном для нас временном масштабе (обычно выбирают самый медленный, а в качестве естественной нормировки для концентраций выбирают интегралы системы). 3. Привести все к одинаковому порядку в правой части и выделить малые параметры в левой. Если нужно, для этого можно складывать и вычитать уравнения. Более того, это НУЖНО делать для выделения истинных переменных. В итоге на основании параметров оценить временную иерархию в системе. 4. Изучить динамику системы быстрых движений на предмет устойчивости. 5. Затем заняться системой медленных движений. 6. Проходя от самых быстрых масштабов к самым медленным, последовательно применить теорему, превращая по мере необходимости дифференциальные уравнения в алгебраические.

Преимущества строгих методов редукции 1. Они применимы для сложных систем любого устройства. Главное, чтобы их можно было описать ОДУ. 2. Мы не просто избавляемся от лишних ДУ, но еще и получаем четкую картину временных масштабов. 3. При этом обнаруживаются новые, априори неочевидные переменные, которые контролируют динамику системы. 4. В результате важные биологические результаты получаются с помощью чисто математических процедур, которые можно в значительной степени формализовать.

Проблемы: ДУ в ЧП

Задачи 1. 2. 3. 4. 5. Ингибирование первого порядка Превращение с необычной стехиометрией Трехступенчатое превращение Связывание-развязывание

Приложения

Теорема Тихонова

Теорема Тихонова: простая формулировка