ДЕКОМПОЗИЦИЯ СИСТЕМ ЭНЗИМНОЙ КИНЕТИКИ С БЫСТРЫМИ И МЕДЛЕННЫМИ ПЕРЕМЕННЫМИ В ЗАДАЧЕ СУИЦИДНОГО СУБСТРАТА
- Авторы: Сметанников М.А.1
-
Учреждения:
- Самарский национальный исследовательский университет имени академика С.П. Королева
- Выпуск: Том 27, № 3 (2021)
- Страницы: 83-88
- Раздел: Статьи
- URL: https://journals.ssau.ru/est/article/view/10513
- DOI: https://doi.org/10.18287/2541-7525-2021-27-3-83-88
- ID: 10513
Цитировать
Полный текст
Аннотация
В данной статье рассматривается пример кооперативного явления, субстрат в котором известен как субстрат «смертника», или суицидный субстрат, потому что он связывается с активным ферментом как субстрат, но фермент превращает его в ингибитор, который необратимо инактивирует фермент. Таким образом, фермент «совершает самоубийство». Цель работы состоит в применении метода интегральных многообразий к редукции системы кинетики суицидного субстрата. В данной статье приводятся результаты применения методов декомпозиции и интегральных многообразий к системам кинетики суицидного субстрата и сравнения решений исходной и конечной систем. Сравнения решений для четырех уравнений приводятся графически, графики созданы посредством программы Microsoft Excel.
Полный текст
Введение
Среди главных инструментов прикладного анализа особая роль отводится методам возмущений, идея которых заключается в выделении основной и детализирующей структур сложной системы. При этом детализирующая структура рассматривается как возмущение основной. Анализ основной структуры сво- дится к рассмотрению существенно более простых моделей и часто сопровождается понижением раз- мерности модели. Поведение исходной модели изучается путём композиции результатов раздельного ис- следования основной структуры и детализирующих факторов [1–10]. В моделях химической кинетики наличие малого параметра связано с тем, что в химической системе одновременно происходят резко отличающиеся по скорости процессы. Имеется значительное число публикаций по теории и приложени- ям как методов упрощения моделей макроскопической кинетики [11–13], так и моделирования крити-
Сметанников М.А. Декомпозиция систем энзимной кинетики с быстрыми и медленными переменными ...
84 Smetannikov M.A. Decomposition of enzyme kinetics system with fast and slow variables in suicide substrate problem
ческих явлений [14–16]. При этом большое разнообразие задач сочетается со сравнительно небольшим арсеналом применяемых средств анализа и довольно распространенным мнением, что эти задачи не имеют ничего общего как по своей постановке, так и по методам решения. Понижение размерности моделей является основным приемом исследования сложных систем любой природы, а критические яв- ления исключительно важны и сами по себе, и как инструмент познания сложных процессов. Авторы книг [1] и [3], основываясь на геометрической теории сингулярных возмущений, предложили подход, позволяющий с единых позиций этой теории рассматривать и методы редукции кинетических систем, и методы математического моделирования критических явлений в таких системах. В энзимной кинети- ке важность субстратов смертников во многом определяется их возможностью нацелить специфичный фермент на инактивацию. Они особенно полезны при введении лекарственных средств, из соображения о том, что они не вредны в их общей форме и только обозначенный фермент может превратить их в их ингибирующую форму. В качестве примеров их важности субстраты смертников были исследованы для применения при лечении депрессии, эпилепсии и некоторых опухолей.
Постановка задачи
Ферментная система [2] включает в себя переменные S, P, E, X и Y , обозначающие субстрат, про- дукт, ингибитор и промежуточные продукты фермента и субстрата, после обезразмеривания системы сменяющиеся на s, ei, ξ, ζ.
Исходная система, рассматриваемая в данной статье, имеет вид:
ds(t) dt
= −s((ϵp + 1) − ϵpξ − (ϵp + 1)ζ −
(ϵp + 1)(η − ϕζ) ) +
ψ
ρ
ξ, (1.1)
1 + ρ
dη(t)
=
dt
pϕ
(1 + ϵp)(1 + ρ)
ξ, (1.2)
dξ(t)
ϵ
dt
= s((ϵp + 1) − ϵpξ − (ϵp + 1)ζ −
(ϵp + 1)(η − ϕζ)
ψ
) − ξ, (1.3)
с учётом замен:
И с начальными условиями:
dζ(t)
ϵ
dt
ϵp
=
(1 + ϵp)(1 + ρ)
σ = ϵp,
η = ϕζ + ψei.
ξ − ψζ (1.4)
s(0) = 1, ξ(0) = 0, ζ(0) = 0, ei(0) = 0.
Исходная и конечная системы с числовыми коэффициентами
При значениях коэффициентов ρ = 1 ,σ = 1 ,ψ = 3 ,p = 1 исходная система (1.1)–(1.4) выглядит
следующим образом:
ds(t)
3 16
4
(ϵ + 1)(η − 0, 125ζ)
dt = −s((ϵ + 1) − ϵξ − (ϵ + 1)ζ −
0, 75
) + 0, 25ξ,
dη(t)
=
dt
0, 09375
ξ,
(1 + ϵ)
dξ(t)
ϵ
dt
= s((ϵ + 1) − ϵξ − (ϵ + 1)ζ −
(ϵ + 1)(η − 0, 125ζ) 0, 75
) − ξ,
dζ(t)
ϵ =
dt
0, 75ϵ
(1 + ϵ)
ξ − ψζ.
Вид системы после применения метода интегральных многообразий:
3 5 13
1 2 7 2
21 2
w˙ 1 = w1w2 − 4 w1 + ϵ( 6 w1w2 − 16 w1 + 3 w1w2 − 8 w1 w2 + 32 w1 ),
1 3
1 2 3
9
17 2
17 2
w˙ 2 = − 8 w1w2 + 32 w1 + ϵ( 8 w1w2 − 16 w1w2 + 128 w1 + 192 w1 w2 − 256 w1 ),
ϵz˙ = B(w + ϵH(t, w, z, ϵ), t, ϵ)z.
Вестник Самарского университета. Естественнонаучная серия. 2021. Том 27, № 3. С. 83–88
Vestnik of Samara University. Natural Science Series. 2021, vol. 27, no. 3, pp. 83–88 85
Численное сравнение решений
На рис. 3.1–3.4 приведено численное сравнение решений исходной и конечной систем переменных при значении малого параметра ϵ = 0.1.
Рис. 3.1. Сравнение решений для первого уравнения задачи до и после построения интегрального многообразия при ϵ = 0, 1
Fig. 3.1. Difference between results of the first equation before using and after two methods ϵ = 0, 1
Рис. 3.2. Сравнение решений для второго уравнения задачи до и после построения интегрального многообразия при ϵ = 0, 1
Fig. 3.2. Difference between results of the second equation before and after using two methods ϵ = 0, 1
Рис. 3.3. Сравнение решений для третьего уравнения задачи до и после построения интегрального многообразия при ϵ = 0, 1
Fig. 3.3. Difference between results of the third equation before and after using two methods ϵ = 0, 1
Сметанников М.А. Декомпозиция систем энзимной кинетики с быстрыми и медленными переменными ...
86 Smetannikov M.A. Decomposition of enzyme kinetics system with fast and slow variables in suicide substrate problem
Рис. 3.4. Сравнение решений для четвертого уравнения задачи до и после построения интегрального многообразия при ϵ = 0, 1
Fig. 3.4. Difference between results of the fourth equation before and after using two methods ϵ = 0, 1
Заключение
В данной статье рассмотрены математические модели задачи динамики суицидного субстрата. К представленной в безразмерной форме математической модели применяется метод интегральных многообразий, существенно упрощающий сложность вычислительных операций. Сравнение численных решений исходной задачи и задачи пониженной размерности при значении малого параметра ϵ = 0, 1 приводится графически.
Об авторах
М. А. Сметанников
Самарский национальный исследовательский университет имени академика С.П. Королева
Автор, ответственный за переписку.
Email: ssmetannikoff@gmail.com
аспирант кафедры дифференциальных уравнений и теории управления
РоссияСписок литературы
- [1] Соболев В.А., Щепакина Е.А. Редукция моделей и критические явления в макрокинетике. Москва: ФИЗМАЛИТ, 2010. 320 с. URL: https://elibrary.ru/item.asp?id=21326259. EDN: https://elibrary.ru/ryrtfh.
- [2] Murray J.D. Mathematical Biology I. An Introduction. New York: Springer, 2001. 551 p. URL: https://booksee.org/book/1008393.
- [3] Воропаева Н.В., Соболев В.А. Геометрическая декомпозиция сингулярно возмущенных систем. Москва: ФИЗМАЛИТ, 2009. 256 с. URL: https://booksee.org/book/1500476.
- [4] Стрыгин В.В., Соболев В.А. Разделение движений методом интегральных многообразий. Москва: Наука, 1988. 256 с. URL: https://booksee.org/book/483890.
- [5] Гольдштейн В.М., Соболев В.А. Качественный анализ сингулярно возмущённых систем. Новосибирск: Ин-т математики АН СССР, Сиб. отд-ние, 1988. 154 с.
- [6] Щепакина Е.А. Интегральные многообразия, траектории-утки и тепловой взрыв // Вестник Самарского гос. университета. 1995. Спец. вып.
- [7] Shchpakina E., Sobolev V. Integral manifolds, canards and black swans // Nonlinear Analysis: Theory, Methods & Applications, 2001, vol. 44, issue 7, pp. 897–908. DOI: http://doi.org/10.1016/S0362-546X(99)00312-0.
- [8] Sobolev V.A. Integral manifolds and decomposition of singulary perturbed system // Systems & Control Letters. 1984. Vol. 5, Issue 3, pp. 169–179. DOI: http://doi.org/10.1016/S0167-6911(84)80099-7.
- [9] Митропольский Ю.А., Лыкова О.Б. Интегральные многообразия в нелинейной механике. Москва: Наука, 1973. URL: https://booksee.org/book/789024.
- [10] Knobloch H.-W., Aulbach B. Singular perturbations and integral manifolds // J. Math. Sci. 1984. V. 18. No. 5.
- [11] Seiler N., Jung M.J., Koch-Weser J. Enzyme-activated Irreversible Inhibitors. Elsevier/North-Holland, Oxford, 1978.
- [12] Walsh C.T. Suicide substrates, mechanism-based enzyme inactivators: recent developments // Annu. Rev. Biochem., 1984. Vol. 53. P. 493–535. DOI: http://doi.org/10.1146/annurev.bi.53.070184.002425.
- [13] Berding C., Keymer A.E., Murray J.D., Slater A.F.G. The population dynamics of acquired immunity to helminth infections // J. Theor. Biol., 1986. Vol. 122. Issue 4. P. 459–471. DOI: https://doi.org/10.1016/s0022-5193(86)80186-2.
- [14] Бобылев Н.А., Емельянов С.В., Коровин С.К. Геометрические методы в вариационных задачах. Москва: Магистр, 1998. 658 с.
- [15] Емельянов С.В., Коровин С.К., Мамедов И.Г. Структурные преобразования и пространственная декомпозиция дискретных регулируемых систем – метод квазирасщепления // Техн. кибернетика. 1986. № 6. С. 118–128.
- [16] Коровин С.К., Мамедов И.Г., Мамедова А.П. Равномерная по малому параметру устойчивость и стабилизация дискретных сингулярно возмущенных динамических систем // Техн. кибернетика. 1989. № 1. С. 21–29.