Математическая модель процесса биодеструкции парацетамола как открытой системы Текст научной статьи по специальности «Математика»

УДК 531/534:[57+61]

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОЦЕССА БИОДЕСТРУКЦИИ ПАРАЦЕТАМОЛА КАК ОТКРЫТОЙ СИСТЕМЫ

Е.В. Вихарева1, А.А. Селянинов2, Ю.Л. Данилов1, И.П. Рудакова1, Т.А. Нечеухина3,

И.Б. Ившина4, Ю.И. Няшин2

1 Кафедра фармацевтической технологии Пермской государственной фармацевтической академии, Россия, 614990, Пермь, ул. Ленина, 48, e-mail: perm@pfa.ru

2 Кафедра теоретической механики Пермского государственного технического университета, Россия, 614990, Пермь, Комсомольский проспект, 29, e-mail: rrn@theormech.pstu.ac.ru

3 Кафедра микробиологии и иммунологии Пермского государственного университета, Россия, 614990, Пермь, ул. Букирева, 15, e-mail: info@psu.ru

4 Институт экологии и генетики микроорганизмов УрО РАН, Россия, 614081, Пермь, ул. Голева, 13, e-mail: ivshina@iegm.ru

Аннотация. Разработана математическая модель процесса биодеструкции парацетамола с учетом особенностей открытой системы. Модель представляет собой четыре кинетических уравнения, два из них решены аналитически, остальные - численно. В виде линейной зависимости от времени определены коэффициенты к и к2 - «константы» скорости реакций биокаталитического окисления парацетамола и его первичного метаболита п-аминофенола соответственно. При этом расчетные данные согласуются с опытными, полученными в процессе биодеструкции парацетамола свободными и иммобилизованными клетками родококков. Математическая модель позволяет прогнозировать поведение исходного субстрата и продуктов его разложения при заданной схеме реакций и известных константах скоростей, а также окончание процесса биодеструкции парацетамола по убыли п-аминофенола.

Ключевые слова: парацетамол, п-аминофенол, биодеструкция, бактерии рода Rhodococcus, константа скорости реакции, кинетические показатели.

Введение

По данным авторов [3], процесс биодеструкции парацетамола с использованием клеток актинобактерий рода КЪойососст осуществляется по схеме, приведенной на рис. 1. При этом гидролиз парацетамола, последующее дезаминирование и окисление и-аминофенола протекают в открытой системе, в которой происходит обмен реагентов и продуктов реакций с окружающей средой (воздухом и культуральной жидкостью). В реальных процессах ферментации парацетамола значения констант скорости указанных реакций являются функциями времени. Данные показатели могут быть подобраны эмпирически или определены экспериментально [2].

Первичным продуктом разложения парацетамола является и-аминофенол, который обнаруживается в среде культивирования родококков в течение первых часов после начала процесса биодеструкции парацетамола и накапливается до конца его ферментации в культуральной жидкости (рис. 2).

Особо следует отметить, что и-аминофенол является наиболее токсичным (минимальная подавляющая концентрация составляет 12 мг/мл) для клеток родококков

© Вихарева Е.В., Селянинов А.А., Данилов Ю.Л., Рудакова И.П., Нечеухина Т.А., Ившина И.Б., Няшин Ю.И., 2008

HOH

т

CH3COOH

NHCOCH3

парацетамол

OH

O

Т

NH

COO H O O H

NH

«-аминофенол пирокатехин

NH

O

OH

О

ІІ

O

OH

гидрохинон

О

бензохинон

цис, цис-муконовая кислота

Рис. 1. Схема биодеструкции парацетамола с использованием клеток R. rubber ИЭГМ 77

/, сут

Рис. 2. Изменение содержания парацетамола (х1) и и-аминофенола (2) в процессе биодеструкции парацетамола клетками Я. вгуЖгороИ^' ИЭГМ 767

по сравнению с парацетамолом (минимальная подавляющая концентрация составляет 500 мг/мл) и другими продуктами его разложения. В связи с этим имеются различия в процессах биодеструкции и-аминофенола как единственного источника углерода и энергии для бактериальных клеток и как продукта разложения парацетамола. Данное обстоятельство необходимо учитывать при прогнозе продолжительности процесса биодеструкции парацетамола по убыли и-аминофенола. Возможность прогнозирования содержания продукта реакции во времени обеспечивает математическая модель, построенная на основании системы кинетических уравнений.

В связи с этим целью настоящих исследований является разработка математической модели процесса биодеструкции парацетамола, которая учитывает

особенности открытой системы и позволяет прогнозировать поведение объекта моделирования во времени.

Экспериментальное исследование процесса биодеструкции парацетамола

В лабораторных условиях эксперименты по биодеструкции парацетамола проводили в колбах Эрленмейера объемом 250 мл, закрепленных на платформе орбитального шейкера (160 об/мин) при 28 °С. В колбы помещали минеральную среду состава (г/л): K2HPO4 - 2,0; KH2PO4 - 2,0; KNO3 - 1,0; (NH4^SO4 - 2,0; NaCl - 1,0; М^04х7Н2О - 0,2; СаС12*2Н2О - 0,02; БеС1зх7Н2О - 0,001, парацетамол в виде фармацевтической субстанции в концентрации 0,2% (0,2 г на 100 мл среды) и бактериальную суспензию в концентрации 107 клеток/мл. В экспериментах по биодеструкции n-аминофенола в качестве исходного субстрата первоначальная масса его варьировалась от 0,075 до 0,14 г. Последнее значение соответствует максимальному количеству данного вещества, которое может образоваться при разложении 0,2 г парацетамола. В работе использовали предварительно отобранные путем скрининга наиболее активные в отношении парацетамола штаммы родококков R. erythropolis ИЭГМ 767 и R. ruber ИЭГМ 77 [4].

Для активизации оксигеназ бактериальные клетки до начала процесса биодеструкции выращивали в минеральной среде вышеуказанного состава в присутствии 0,1% фенола - структурного аналога парацетамола. В опытах по биодеструкции парацетамола использовали свободные и иммобилизованные (закрепленные на носителе) клетки родококков. В качестве носителя использовали криогель поливинилового спирта [7]. Продолжительность экспериментов составляла от 12 часов до 8 суток. Качественный состав продуктов разложения анализировали методом тонкослойной хроматографии. Количественные характеристики процесса биодеструкции парацетамола изучали по приросту первичного метаболита и-аминофенола и степени исчезновения парацетамола. Содержание и-аминофенола в культуральных жидкостях определяли фотометрическим методом с помощью спектрофотометра Lambda EZ 201, («Perkin-Elmer», США) при 450 нм по

модифицированной методике, основанной на реакции конденсации и-аминофенола с и-диметиламинобензальдегидом [3]. Количество остальных метаболитов, содержащихся в осадке (шламе), рассчитывали с учетом закона сохранения молекулярных масс веществ. С целью определения класса опасности шлама его токсичность (LD50) исследовали при пероральном введении белым беспородным мышам в дозах от 1000 до 4000 мг/кг массы животного [9].

Класс опасности культуральной жидкости и шлама определяли расчетным методом в соответствии с «Критериями отнесения опасных отходов к классу опасности для окружающей природной среды» [5] на основании известных токсикологических [6] и эколого-гигиенических [1] показателей и-аминофенола и значения LD50 осадка (шлама). Концентрацию парацетамола в культуральных жидкостях определяли методом высокоэффективной жидкостной хроматографии на хроматографе модели HP-1090 (Hewlett Packard, США) с использованием колонки ZORBAXExtend-C18 (4,6*250 мм). Эксперименты проводили в трехкратной повторности. Математическую обработку полученных результатов осуществляли с использованием компьютерных программ Excel2000 и Statistica (версия 3.4 для Windows).

Кинетические схемы процессов биодеструкции парацетамола

И М-АМИНОФЕНОЛА В КАЧЕСТВЕ ИСХОДНЫХ СУБСТРАТОВ

Процесс биодеструкции парацетамола ранее описан авторами [1] кинетическим уравнением, в котором коэффициент скорости реакции [2] биокаталитического окисления парацетамола представлял собой константу к1. Для построения более адекватной модели данного процесса предлагается считать к1 функцией времени.

-х, 7 . .

—і = -кіО) • хі, (1)

-ї

г

где х1 = 0 = х10 = 0,002 — - начальная концентрация парацетамола.

мл ]

Процесс биодеструкции и-аминофенола в качестве исходного субстрата может быть описан аналогичным кинетическим уравнением

—х

~Хг = - к2(ї) • Х2, (2)

где хА = х2 = 0,00075

2 її=0 2

мл

-ї

- начальная концентрация и-аминофенола,

соответствующая максимальному содержанию данного вещества как продукта разложения в процессе биодеструкции парацетамола. Необходимо отметить, что в концентрациях более 0,001 г/мл и-аминофенол вызывает гибель бактериальных клеток.

Уравнения (1) и (2) интегрируются методом разделения переменных.

Для уравнения (1)

х Сх 1 1 1

[ —1 = -[ к1 (г)С, отсюда 1п Ы 0 = — [ к1 (г)Сг ■ г.

у -у* « 1x1 / «I

Х0 10 1 0

1 Г

Для среднеинтегрального выражения -1 к1 (г)С используем линейную

г

0

к

аппроксимацию по времени в виде к[ = г + Ь .

Из физических соображений значения концентраций положительны, поэтому

модуль опускаем. Тогда 1п-Х- = -(Ь + кг)г и окончательно кинетическое уравнение

Х0 2

биодеструкции парацетамола принимает вид

-( Ь+1 г

Х = х”■е ' 2 ' . (3)

Аналогично для процесса биодеструкции и-аминофенола в качестве исходного субстрата

-І с+—ї Ії

г

t, сут

Рис. 3. Зависимость коэффициента к[ от времени в процессе биодеструкции парацетамола клетками R. ruber ИЭГМ 77 (медленный процесс).

к', сут-1

♦ Экспериментальные

точки

- Линейная

аппроксимация

• Среднее значение

k' = 0,22841 + 0,1459

г, сут

Рис. 4. Зависимость коэффициента к1' от времени в процессе биодеструкции парацетамола клетками Я. егугкгороИ^' ИЭГМ 767 (средний процесс).

К, ч-1

♦ Экспериментальные точки

---- Линейная

аппроксимация

• Среднее значение

k = 0,01871 + 0,0045

t, ч

Рис. 5. Зависимость коэффициента к[ от времени в процессе биодеструкции парацетамола иммобилизованными клетками Я. егугкгороИ^' ИЭГМ 767 (быстрый процесс)

В работе приведены определенный методом наименьших квадратов к

коэффициент к[ = — t + b линейной зависимости от времени и средние значения

«константы» к/ скорости реакции разложения парацетамола с использованием

свободных и иммобилизованных клеток родококков (рис. 3-5).

При этом следует отметить, что процесс биодеструкции парацетамола с использованием свободных клеток R. ruber ИЭГМ 77 (см. рис. 3), R. erythropolis ИЭГМ 767 (см. рис. 4) и иммобилизованных клеток R. erythropolis ИЭГМ 767 (см. рис. 5) протекает с разной скоростью и условно обозначен нами как медленный, средний и быстрый соответственно. Как видно из рис. 4, коэффициент к[ имеет выраженную линейную зависимость от времени только в условиях среднего процесса. Однако во всех остальных случаях авторами использованы также линейные зависимости к[(t),

которые позволяют получить соответствие теоретических и экспериментальных данных.

л

|н

,л

о

н

е

ф

о

н

и

к'2, сут-

t, сут а

к2 = 0,1329t + 0 ,1571

t, сут

б

Рис. 6. Зависимость концентрации и-аминофенола (а) и коэффициента к'2 (б) от времени в

процессе биодеструкции его как исходного субстрата клетками

Я. erythropolis ИЭГМ 767

Таблица

Параметры «констант» скорости реакций разложения парацетамола k1(t) и и-аминофенола k2(t)

Процесс биодеструкции парацетамола Значения параметров k1(t) и k2(t)

k b d c

Медленный (с использованием свободных клеток R. ruber ИЭГМ 77) -0,0104 [сутки-2] 0,5027 [сутки-1] 0,095 [сутки-2] -0,11 [сутки-1]

Средний (с использованием свободных клеток R. erythropolis ИЭГМ 767) 0,4568 [сутки-2] 0,1459 [сутки-1] 0,115 [сутки-2] -0,15 [сутки-1]

Быстрый (с использованием иммобилизованных клеток R. erythropolis ИЭГМ 767) 0,0374 [час-2] 0,0045 [час-1] 0,01 [час-2] -0,04 [час-1]

На рис. 6 представлено изменение содержания и-аминофенола и коэффициента к2 от времени в процессе биодеструкции и-аминофенола как исходного субстрата клетками R erythropolis ИЭГМ 767. Предложенное кинетическое уравнение (2) адекватно описывает процесс разложения и-аминофенола (см. рис. 6, а), а коэффициент к'2 имеет явно линейную зависимость от времени (см. рис. 6, б), как и коэффициент к[ в процессе биодеструкции парацетамола клетками R. erythropolis ИЭГМ 767 (см. рис. 4).

При переходе от среднеинтегрального выражения к «константе» скорости реакции к1 (t) получена ее зависимость от времени в виде к1 (t) = b + к • t для парацетамола и к2^) = с + d• t для и-аминофенола соответственно. Параметры ^(t) (коэффициенты b, к) и к2^) (коэффициенты c, d), полученные эмпирическим путем,

приведены в таблице.

Как видно из таблицы, для медленного (с использованием свободных клеток R ruber ИЭГМ 77) и среднего (с использованием свободных клеток R erythropolis ИЭГМ 767) процессов биодеструкции парацетамола параметры к2^) имеют близкие

значения. Параметры к2 (t) (с = 0,1571 и d/2 = 0,1329), полученные при разложении и-аминофенола в качестве исходного субстрата, достоверно близки к эмпирически найденным для к2 (t) в процессе биодеструкции парацетамола свободными клетками

R erythropolis ИЭГМ 767. В дальнейшем эмпирически найденные параметры к2^)

использовались авторами при разработке математической модели всего процесса биодеструкции парацетамола.

Математическая модель процесса биодеструкции парацетамола

Как было указано выше (см. рис. 1), процесс биодеструкции парацетамола представляет собой открытую систему. При этом скорость разложения парацетамола как исходного субстрата пропорциональна его концентрации [3]. Знак «-» в уравнении (1) означает убывание концентрации парацетамола со временем. Коэффициент к1

является линейной функцией времени. Кинетическое уравнение для парацетамола как для исходного субстрата сохраняет вид (1).

На первой стадии процесса биодеструкции парацетамола (концентрация х1) в результате реакции его гидролиза (+Н2О) образуются и-аминофенол (концентрация Х2) и уксусная кислота (концентрация Х3). Суммарная скорость образования данных

соединений увеличивается со скоростью разложения парацетамола (со знаком «+») и уменьшается за счет текущего разложения и-аминофенола. Кинетическое уравнение принимает вид

ё (х2 + х3)

= кО)• х1 -^0)• Х2 *

&

Поскольку в реакции гидролиза участвует вода (система открытая), при постановке задачи биодеструкции парацетамола в терминах концентраций следует учесть прирост массы воды (11,9 %) в виде поправочного коэффициента Кн2о. Расчёт

Хд0 может быть осуществлён на основании молекулярных масс, участвующих

в реакции соединений.

X

+ СНзСООН

МН2

Х1 Х2 Хз

(молекулярная масса) 151,2 + 18,0 ^ 109,1 + 64,0 . Отсюда коэффициент прироста массы за счет воды

18 0

Кн о = 1,0 +---------= 1,119 .

Н2<О 151,2

Для суммы концентраций х2 и х3 получаем кинетическое уравнение

= Кн2о • ^(0• Х1 -К(^• Х2 , (5)

ё (Х2 + Х3)

&

где Х2 I=0 = 0 Х3 1= 0 = 0.

Доля уксусной кислоты составляет 35,5 % от суммы молекулярных масс

парацетамола и воды. В связи с этим коэффициент убыли масс за счет уксусной

кислоты

К = 60,0 = 0 355 сн3соон 151,2 +18,0 , .

Кинетическое уравнение изменения содержания уксусной кислоты

(концентрация Х3 ) принимает вид

&Х

= Ксн3соон • Кн2о • к () • Х1 . (6)

Для прочих продуктов биодеструкции (в частности, пирокатехина и других в суммарной концентрации Х4 ) следует учесть прирост молекулярной массы кислорода и четырех протонов в количестве 33% от молекулярной массы и-аминофенола.

(молекулярная масса) 109,1 + 32,0 + 4,0 ^ 110,1 + 17,0 + 18,0.

Исходя из значений молекулярных масс соединений, участвующих в реакции окисления и-аминофенола, находим коэффициент прироста массы.

, 32,0 + 4,0 ,

Ко = 1 +-----------------= 1,330 .

°2 109,1

Кинетическое уравнение для остальных метаболитов принимает вид

йх4

йі

= Ко2 • К(і)• Х2

(7)

где Х41=0 = 0.

В результате система кинетических уравнений процесса биодеструкции парацетамола принимает вид

&Х1 ,, , Л

—1 = -(Ь + к •)-Х1, т

й (х2 + х3) йх?

йі

= Кн о • (Ь + к • і) • х1 - (с + й • і) • х

2^

(8)

йі

йх4

йі

= К

сн3соон ^'-н2о

= Ко2 • (с + й • і) • х2

Кн,о • (Ь + к •і) • х1

где х, = х0, х2 = х3 = х4 = 0.

^ 1 Іі=0 15 2 Іі=0 3 Іі= 0 4 Іі=0

Іі=0

і = 0

Іі=0

Система (8) включает четыре дифференциальных уравнения 1-го порядка с известными начальными условиями. Попытаемся решить данную систему аналитически. Первое уравнение в (8) проинтегрировано ранее (3),

Ь+кі |і

х = х10 • е ^ 2

Подставив х в третье уравнение системы (8), получаем

Шх3

Ш

— К

СН3С00Н -11-Н,0

КН20 • (Ъ + к • г) • X0 •е

Интегрируем это дифференциальное уравнение методом разделения переменных,

-"-3 А

^ Щх3 — КСН3С00Н • КН20 • Х0 • ^ (Ъ + к • *) • е

Ъ+-2-г | г

л.

отсюда

Х — К • К • Х

л3 -,'-СН3С00Н -''-Н20 л0

0 - е

(9)

Далее подставляем третье уравнение из системы (8) во второе уравнение системы (8),

— - К

• КН 0 • (Ъ+к • г) • х0 + КН 0 • (Ъ+к • г) • х0 - (с+ш • г) • х2 —

л ^НзС00Н

— (0 — КСН3С00Н ) • КН20 • (Ъ + к • г) • Х0 — (с + Щ г) • Х2'

Заменив х0 через зависимость (3), получаем неоднородное дифференциальное уравнение 0-го порядка для х2 с коэффициентом в виде линейной функции времени

Щх2

ш

+ (с + Щ • г) • х2 — (0 - КС

,) • КН20 • (Ъ+к • г) • Хо0 • е "СЪ+2 ‘к

(00)

Необходимо отметить, что после интегрирования уравнения (00) и подстановки результата в четвертое уравнение системы (8) для остальных метаболитов х4 можно получить дифференциальное уравнение 0-го порядка с разделяющимися переменными.

Для концентрации парацетамола х0 решение из системы кинетических уравнений (8) имеет естественный вид экспоненты (3). Для концентрации и-аминофенола х2 в процессе биодеструкции его как исходного субстрата решение

имеет такой же естественный вид (4).

Однако при разложении и-аминофенола в процессе биодеструкции парацетамола дифференциальное уравнение (00) не интегрируется в элементарных функциях. В связи с этим представляется необходимым решение системы (8) с применением численного метода. Для этого систему (8) преобразуем к стандартному виду для реализации на компьютере,

Лх0 ч

_! — -(Ъ + к •) •х,,

ш Шх

-Щ— — (1 — КСН3С00Н ) • КН20 • (Ъ + к • г) • Х0 - (с + Ш г) • Х2 ,

Шх3

(00)

— КС

• КН20 • (Ъ + к • г) • Xо,

Щх

—4 — К0 • (с+ш • г) • х шг 02

2

Первое и третье уравнения проинтегрированы ранее (3), (9). Поэтому, подставив полученную функцию х0(* ) во второе уравнение системы (00), получаем систему двух уравнений

где х2 г—0 — х4 г—0 — 0; /2, /4 - известные функции от времени г и концентрации и-аминофенола х2.

Одним из наиболее удобных и точных методов численного интегрирования дифференциальных уравнений является метод Рунге-Кутта, для которого используем стандартную программу интегрирования системы двух дифференциальных уравнений 0-го порядка.

Сравнение теоретических и экспериментальных данных

Результаты численного решения системы (02) для компонентов реакции х2 и х4 методом Рунге-Кутта и аналитические зависимости для компонента х0 в сравнении с экспериментальными данными представлены на рис. 7-9. Как видно из рис. 7-9, разработанная модель адекватно описывает процесс биодеструкции парацетамола как открытой системы. Содержание парацетамола х0 в процессе биодеструкции уменьшается сначала быстро, а затем медленно. При использовании клеток Я. гиЪег ИЭГМ 77 исчезновение парацетамола из среды культивирования родококков наблюдается через 8 суток (см. рис. 7), Я етугктороИя ИЭГМ 767 - через 4 суток (см. рис. 8), а при использовании иммобилизованных клеток Я. егугкгороНя ИЭГМ 767 -через 04 часов (см. рис. 9). Различие в скорости процесса биодеструкции парацетамола при использовании свободных и иммобилизованных клеток родококков зависит от каталитической активности клеточных ферментов, которая инициируется предварительной индукцией оксигеназного комплекса и в еще большей степени иммобилизацией бактериальных клеток.

Содержание и-аминофенола х2 как продукта разложения парацетамола сначала

достаточно быстро увеличивается, достигая максимальных значений в культуральной жидкости на 3-и сутки при использовании свободных клеток родококков (см. рис. 7, 8) или через 02 часов при использовании иммобилизованных клеток (рис. 8), а затем плавно уменьшается до нуля на 8-е сутки или через 25 часов соответственно. Разработанная модель дает возможность прогнозировать окончание процесса

(12)

Формализованная запись системы (02) принимает вид

(13)

t, сут

Рис. 7. Экспериментальные (точки) и теоретические (линии) зависимости от времени

концентрации парацетамола х1, и-аминофенола х2 и остальных продуктов разложения х4

в процессе биодеструкции парацетамола клетками R. ruber ИЭГМ 77 (медленный процесс). Пунктирными линиями обозначены время окончания процесса и соответствующее ему содержание и-аминофенола.

биодеструкции парацетамола по убыли и-аминофенола. Так, по нашим данным, парацетамол и и-аминофенол относятся к I классу опасности для окружающей природной среды, то есть являются чрезвычайно опасными отходами. Культуральная жидкость становится малоопасным отходом (IV класс опасности для окружающей природной среды) и может быть размещена на санитарном полигоне по достижении в ней концентрации и-аминофенола, равной 0,0005 г/мл. Следовательно, процесс разложения парацетамола при использовании свободных клеток родококков может быть завершен через 6 суток (см. рис. 7, 8), а при использовании иммобилизованных клеток - через 19 часов (см. рис. 9).

Содержание уксусной кислоты х3 как продукта реакции гидролиза парацетамола и практически не опасного в данной концентрации отхода в ходе экспериментов не регистрировалось, и на рис. 7-9 теоретические зависимости ее концентрации от времени не приводятся. Концентрация x4 остальных продуктов разложения

парацетамола, находящихся в осадке (шламе), заметно увеличивается, начиная с первых суток (см. рис. 7, 8) или нескольких часов (см. рис. 9) и достигает во всех случаях в конце процесса около 85% от начальной концентрации парацетамола. Образующийся осадок малотоксичен (LD50 > 4000 мг/кг), относится к V классу опасности (практически не опасен) для окружающей природной среды и наряду с культуральной жидкостью может быть размещен на санитарном полигоне.

Заключение

Таким образом, в результате проведенных исследований разработана математическая модель процесса биодеструкции парацетамола с учетом особенностей открытой системы. Модель представляет собой четыре кинетических уравнения, два из которых решены аналитически, остальные - численно. В виде линейной зависимости от времени определены коэффициенты k1 и k2 - «константы» скорости реакций

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

£

о „ Л и

к к

* о ч

0 со

^ 5

' Он

ей

&

с

х4

Х2

Х1

і, сут

Рис. 8. Экспериментальные (точки) и теоретические (линии) зависимости от времени концентрации парацетамола хь и-аминофенола х2 и остальных продуктов разложения х4 в процессе биодеструкции парацетамола клетками К. вгуЛгороШ ИЭГМ 767

(средний процесс).

л

Й

Ь ч о -$

к § ч м

о й

св О н ч

^ со

^ 5

сЗ Он

&

с

і, ч

Рис. 9. Экспериментальные (точки) и теоретические (линии) зависимости от времени концентрации парацетамола хь и-аминофенола х2 и остальных продуктов разложения х4 в процессе биодеструкции парацетамола иммобилизованными клетками К. вгуЛгороШ

ИЭГМ 767 (быстрый процесс).

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

Благодарности

Исследования поддержаны грантом РФФИ № 07-04-96038-р-Урал.

Список литературы

1. Вихарева, Е.В. Математическое моделирование процесса биодеструкции парацетамола актинобактериями рода Ккойососсш / Е.В. Вихарева, А.А. Селянинов, И.Б. Ившина, Ю.И. Няшин // Российский журнал биомеханики. - 2007. - Т. 11, № 2. - С. 93-100.

2. Денисов, Е.Т. Кинетика гомогенных химических реакций / Е.Т. Денисов. - М.: Высшая школа, 1988. - С. 390.

3. Ившина, И.Б. Деградация парацетамола с истекшим сроком годности свободными клетками актинобактерий / И.Б. Ившина, М.И. Рычкова, Е.В. Вихарева, Л. А. Чекрышкина, И.И. Мишенина // Катализ в промышленности. - 2006. - № 2. - С. 44- 49.

4. Каталог штаммов Региональной профилированной коллекции алканотрофных микроорганизмов / под ред. И.Б. Ившиной. - М.: Наука, 1994. - 163 с.

5. Колла, В.Э. Дозы лекарственных средств и химических соединений для лабораторных животных / В.Э. Колла, Б.Я. Сыропятов. - М.: Медицина, 1998. - 263 с.

6. Контроль химических и биологических параметров окружающей среды / под ред. Л.К. Исаева. -СПб.: Эколого-аналитический центр «Союз», 1998. - 896 с.

7. Лозинский, В. И. Криотропное гелеобразование растворов поливинилового спирта / В.И. Лозинский // Успехи химии. - 1998. - Т. 67, № 7. - С. 641-652.

8. Об утверждении критериев отнесения опасных отходов к классу опасности для окружающей природной среды: приказ Министерства природных ресурсов Российской Федерации № 511 от 15.06.01.

9. Прозоровский, В.В. Экспресс-метод определения средней эффективной дозы и ее ошибки / Прозоровский В.В., Прозоровская М.П., Демченко В.М. // Фармакология и токсикология. - 1978. - № 4. - С. 497-502.

A MATHEMATICAL MODEL OF PARACETAMOL BIODEGRADATION PROCESS AS AN OPEN SYSTEM

E.V. Vikhareva, A.A. Selyaninov, Yu.L. Danilov, I.P. Rudakova, T.A. Necheukhina,

I.B. Ivshina, Yu.I. Nyashin (Perm, Russia)

A mathematical model of the paracetamol biodegradation process was developed taking into account specific characteristics of an open system. The model consists of four kinetic equations; two of those were solved analytically. At the same time, the other equations were solved numerically. The coefficients k1 and k2 (“constants” of paracetamol and its

primary metabolite ^-aminophenol biocatalytic oxidation rate) were determined as linear dependence on time. At that, theoretical data are in a good agreement with experimental ones obtained during the paracetamol biodegradation process by free and immobilized Rhodococcus cells. The mathematical model enables us to predict the initial substrate and its metabolites behavior according to a given reaction outline and known reaction rate constants. Moreover, it allows prospecting the termination of paracetamol biodegradation process by the ^-aminophenol loss.

Key words: paracetamol, ^-aminophenol, biodegradation, Rhodococcus bacteria, reaction rate constant, kinetic factors.

Получено 01 апреля 2008