Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.library.biophys.msu.ru/mce/20002402.htm
Дата изменения: Fri Jan 17 21:43:17 2003
Дата индексирования: Mon Oct 1 19:43:00 2012
Кодировка: Windows-1251
ПРОБЛЕМЫ МОДЕЛИРОВАНИЯ ФЕРМЕНТАТИВНОЙ РЕАКЦИИ ЭЛЕКТРИЧЕСКИМИ ЦЕПЯМИ

 

ПРОБЛЕМЫ МОДЕЛИРОВАНИЯ ФЕРМЕНТАТИВНОЙ РЕАКЦИИ

ЭЛЕКТРИЧЕСКИМИ ЦЕПЯМИ

 

Бирюк Н.Д., Ковалева Т.А., Юргелас В.В.

 

(Воронеж)

 

 

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

В современной биофизике [1, 2] активно используются достижения радиоэлектроники. Широко используется, например, понятие динамической модели биологического процесса. Эти динамические модели весьма похожи на нелинейные радиоцепи. Применяются также понятия биологического триггера, биологического автогенератора. Биохимическая реакция в длинном и узком сосуде распространяется вдоль длины сосуда волнами и описывается дифференциальными уравнениями в частных производных, аналогичными телеграфным уравнениям длинной линии. Математическое моделирование в биофизике широко используется. Однако математическим моделям в силу их абстрактности присущи известные ограничения, от которых свободны радиофизические модели. Последние позволяют с помощью приборов (осциллографов, самописцев) визуализировать некоторые скрытые особенности реакций.

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

 

(1)

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

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

 

(2)

Это типичная система уравнений для свободных процессов радиоцепей с сосредоточенными, постоянными параметрами, которую можно представить, как одно векторное дифференциальное уравнение первого порядка в компактном виде:

(3)
где - неизвестный вектор-столбец -го порядка , - квадратная матрица -го порядка. Элементы вектора являются аналогами компонент реакции, выражаемых количественно как превышения их текущих компонент над стационарными.

Уравнение (3) описывает свободные колебания системы, под которой в данном случае понимается радиоцепь. Для его решения принципиальное значение имеет отыскание корней характеристического уравнения:

(4)
где - единичная матрица -го порядка. Здесь целесообразно выделить следующие три случая: а) все корни уравнения (4) действительны и различны, т.е. ; б) все корни различны, но среди них есть одна или несколько пар комплексно сопряженных, например, , где ; в) среди корней есть один или несколько кратных, при этом сумма числа различных корней с учетом кратности (корень учитывается столько раз, какова его кратность) равна .

В случае а) каждому корню соответствует решение: