CC BY
37
ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ / PHYSICS AND MATHEMATICS
DOI: https://doi.org/10.23670/IRJ.2017.66.156 Абузярова Р.Т.1, Колпак Е.П.2
1ORCID: 0000-0002-3694-7832, студент, 2ORCID: 0000-0001-6956-4814, доктор физико-математических наук, Санкт-Петербургский государственный университет МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ЛЕЙКОЗА
Аннотация
Разработана математическая модель лейкоза на основании современных данных о кроветворении у млекопитающих. В модели рассматриваются три типа клеток, один из которых лейкозные. Взаимодействие клеток рассматривается как конкуренция за функциональное пространство делящихся клеток. За лейкозные клетки принимаются клетки, делящиеся с большей скоростью, чем все остальные. Модель замены лейкозных клеток клетками донора рассматривается как внедрение более активных клеток, чем лейкозные. Модель представлена задачей Коши для системы обыкновенных дифференциальных уравнений. Нарушение функций кроветворения сопоставляется с изменением ведущего параметра, переводящего систему из устойчивого стационарного состояния в неустойчивое.
Ключевые слова: лейкоз, математическое моделирование, устойчивость, дифференциальные уравнения, задача Коши, стационарная точка.
Abuzyarova R.T.1, Kolpak E.P.2
1ORCID: 0000-0002-3694-7832, Student, 2ORCID: 0000-0001-6956-4814, PhD in Physics and Mathematics, St. Petersburg State University MATHEMATICAL MODEL OF LEUKEMIA
Abstract
We developed a mathematical model of leukemia on the basis of modern data on hematopoiesis in mammals. Three types of cells are considered in the model, one of which is leukemia. Interaction of cells is studied as a competition for the functional space of dividing cells. Leukemia cells are cells dividing at a higher rate than all the others. The model of replacement of leukemic cells with donor cells is considered as the introduction of more active cells than leukemia ones. The model is represented by the Cauchy problem for a system of ordinary differential equations. Violation of hematopoiesis functions is correlated with a change in the leading parameter, which transfers the system from a stable steady state to an unstable one.
Keywords: leukemia, mathematical modeling, stability, differential equations, Cauchy problem, stationary point.
Введение
Одним из острых заболеваний является лейкоз. Ему подвержены все возрастные группы. Это заболевание трудноизлечимо, и может длиться, несмотря на лечение, годами. Очень высок процент летальных исходов. Болезнь протекает с осложнениями. Наиболее распространенным и дорогостоящим методом лечения является пересадка в больной организм донорского костного мозга. Однако чужеродные клетки начинает отторгать иммунная система и, в результате, выход в ремиссию часто сопровождается другими заболеваниями организма. Лечение осуществляется и введением в больной организм препаратов, помогающих справиться с заболеванием. Стоимость лечения препаратами составляет около 100 средних по России месячных зарплат. При благоприятном исходе лечения выход в ремиссию не означает, что не возникнет рецидив. Математическое моделирование лейкоза и вариантов его лечения могут помочь лечащим заболевание специалистам спланировать лечение и по его срокам и по его стоимости.
Кроветворение
Кроветворение - это процесс образования тромбоцитов, эритроцитов и лейкоцитов. Родоначальными клетками гемопоэза являются стволовые клетки (СК), благодаря которым обеспечивается стабильное кроветворение [1]. В результате деления СК могут образовывать два типа клеток: дочерние клетки, являющиеся идентичными родительской клетке, и клетки, которые уходят в созревающую фракцию и начинают путь дифференцировки. Острый лейкоз представляет собой клеточный клон, все элементы которого являются потомками одной клетки [2]. Первоначальная мутация происходит в стволовой клетке. Генетическая нестабильность приводит к дополнительным мутациям, в результате которых клетки образовавшегося клона перестают дифференцироваться, постепенно поражают костный мозг морфологически незрелыми кроветворными клетками с вытеснением ими нормальных клеток.
На поведение СК влияет их микроокружение или их ниша. Ниши представляют собой пространство для СК, в котором контролируется поведение СК и регулируется их численность через сигналы от клеток и неклеточных факторов. СК в зависимости от сигналов могут или пролиферировать, или находиться в состоянии покоя, или дифференцироваться [3]. Дифференцировка и пролиферация кроветворных клеток происходят параллельно. За счет ниш поддерживается равновесие между различными типами кроветворных клеток.
Лейкозные клетки располагаются в тех же нишах что и гемопоэтические СК и прогениторные клетки, в результате чего имеет место конкуренция за нишу. В частности, лейкозные клетки способны избыточно вырабатывать колониестимулирующий фактор, увеличивающий скорость их деления по сравнению со скоростью деления нормальных клеток (стимулятор миелопоэза), специфически сильнее действующий на лейкозные, чем на нормальные клетки-предшественники гемопоэза [1], [2]. В то же время злокачественные клетки, даже в малом количестве, обладают способностью избирательно угнетать пролиферацию и дифференцировку нормальных предшественников с помощью гуморальных ингибиторов
В результате в нише происходит конкуренция за функциональное пространство между различными клетками [3]. Конкурентное преимущество имеют клетки, которые быстрее пролиферируют, давая многочисленное потомство.
Скорость пролиферации СК регулируется сигналами, поступающими от клеток микроокружения. Нарушение функционального состояния ниши приводит к увеличению скорости генерации клеток или к ее уменьшению. Математическая модель лейкоза
Предлагаемые в литературных источниках математические модели [4], [5], представлены системами обыкновенных дифференциальных уравнений. В работах [6], [7] предлагается модель иммунного отклика на появление лейкозных клеток, а в [8] предлагаются варианты моделирования лечения лейкозов. Уровень заболеваний населения РФ заболеваниям крови отражен в [9]. За основу моделирования лейкоза взята концепция конкуренции всех типов клеток за функциональное пространство [10].
В математической модели предполагается, что лейкоз имеет иерархическую структуру. Рассматривается четыре типа клеток: здоровые, лейкозные, стволовые и клетки микроокружения, регулирующие пролиферацию стволовых
и
клеток. Лейкозные, здоровые и стволовые клетки объединены в один пул единичного объема. Если принять, что 5 -
К и,
количество стволовых клеток, находящихся в нише, " - здоровых, а - лейкозных, то их количество на момент
К + К + и^ < 1 „
начала кроветворения должно удовлетворять условию 5 " а . Распределение размножающихся
стволовых клеток между больными и здоровыми происходит под влиянием клеток микроокружения в пропорциях а 1 - & (0 < & < 1). Здоровые, больные и стволовые клетки могут самостоятельно размножаться с удельными
скоростями равными А, А2 и Аз соответственно. Скорости перехода здоровых и больных клеток в кровеносное
V, V- „ и
русло равны 1 и 2 соответственно. Скорость перехода клеток 5 в здоровые и больные регулируется клетками
и
микроокружения а . С учетом этих предположений математическая модель лейкоза представлена задачей Коши для системы четырех обыкновенных дифференциальных уравнений
аг
аи
Оии ,, ч
— = А (1 - ик - иа - К >Ъ + аГиаи5 - V1Uh
а =а2(1 -uh -иа -и5)иа + (1 -а)гиаи5 -
аг
аи„ 2
= А(1 -uh -иа -и5К -УиК<
аг
аи„
(1)
= -Аиа (Р- и5 )
аг
ин (г = 0) = к0 и,(г = 0) = и0 ^ = 0) = и0 иа = 0) = и°а
1_и _и _и
Множитель h а 5 в первых трех уравнениях представляет собой свободную от клеток часть пула, а
Р
параметр соответствует количеству стволовых клеток, находящихся в нише, в равновесном состоянии при естественной генерации всех типов клеток.
Все искомые величины в модели считаются безразмерными. Время г также считается безразмерным, однако, применительно к реальным процессам единица безразмерного времени может быть сопоставлена с сутками в реальном процессе.
и и уи и
Стволовые клетки 5 под влиянием клеток микроокружения а покидают нишу со скоростью ' 5 а и
И и
переходят в профилирующие здоровые ( 5) и лейкозные ( а) клетки. При увеличении количества стволовых клеток
Р
в нише выше порогового значения увеличивается и количество клеток микроокружения, вызывающих увеличение
уи.иа
скорости ухода стволовых клеток из ниши со скоростью ' 5 а, а при уменьшении количества стволовых клеток в нише скорость ухода будет уменьшаться.
-
аи
у = 0 и = 0 и = 0
В случае отсутствия ниши (в (1) '1 , 5 , а ) первые два уравнения в (1) принимают вид
аг
h _
= А(1 -иь -иа)иь ->
^ = А2(1 - К - иа )иа - V2иа> аг (2)
Эта система уравнений имеет две стационарные точки
1.
ик = 1
и,
0
л и. = 1 -2. ">, = 0, '
V-,
/2
В первой стационарной точке собственными значениями матрицы Якоби правой части уравнений (2) будут
■V,
Л VI
2
Л1 = и /1
И эта стационарная точка будет устойчивой, если выполняется неравенство /21 < /1>2 . То
есть, если удельная
и
1н
со временем
скорость генерации опухолевых клеток мала или скорость их ухода в систему велика, то клетки
вытеснят клетки и( . Если это неравенство не выполняется, то устойчивой будет вторая стационарная точка. Модель пересадки клеток
Одним из способов лечения является пересадка донорских клеток. Если их активность выше, чем активность клеток хозяина, то они постепенно вытеснят клетки хозяина. Как показывает практика клинических исследований, донорские клетки не всегда приживаются в новом организме. С учетом этого модель вытеснения клеток хозяина
донорскими клетками будет представлена системой дифференциальных уравнений
(и.
к
(И сИ
= /(1 - ик - иь )ик - Vик ,
= / (1 - Щ - иь К - Vиь - УЩиь,
(3)
в которой
/з
удельная скорость размножения донорских клеток,
V,
- скорость перехода донорских клеток в
кровеносное русло. Ингибирующее влияние клеток хозяина на клетки донора учитывается введением во второе
уравнение слагаемого с параметром У, определяющим интенсивность ингибирования.
Система уравнений (3) имеет три нетривиальные стационарные точки
1.
иь = 0, иь = 1 --1-
ик = 0, иь = 1
ин = —
3.
У
2.
/ VI" 1
/ Уз
, иь
/з ;
:1 " " VI / У
/ VI" 1
/ У3 ,
Третья стационарная точки существует, если выполняются неравенства
Г Л
У
1 </! < 1 1 -
V /1
В первой стационарной точке одно из собственных значений матрицы Якоби правой части уравнений (3) будет положительным, если выполняется неравенство
/
> 1 + у 1 -
V
Л
У .
Это неравенство для модели (3) является условием вытеснения донорскими клетками клеток хозяина.
Вторая стационарна точка будет устойчивой, если выполняется неравенство 1 < . В третьей стационарной точке характеристический полином матрицы Якоби правой части уравнений (3)
Р(Л) = Л + (/к + /зиь )Л - у/ии
имеет корни противоположных знаков. Поэтому эта стационарная точка будет неустойчивой. Таким образом, в модели (3) донорские клетки либо вытесняют клетки хозяина, либо их вытесняют клетки хозяина.
Модель (3) для случая наличия здоровых и лейкозных клеток принимает вид
/ин
Ж
/ис
ж
/иь
Ж
= М(1 - и, - ий - и
ь ин - Vин,
= М2(1 -ин -иа -иь)ис -Vис>
= Мз(1 - ин - ис - иь )иь-У(ин + ис )иь - уъиь •
Решение этих уравнений для случая
VI = V = V = 01 М = 1.0 М= 1.5 М = 2.0 У = 0.007
приведено на
рис. 1. Пунктирной линией отмечен момент времени (* 60) пересадки донорских клеток. В этом примере донорские клетки вытесняют и здоровые и лейкозные клетки.
1
0.8
^ 0.6-
, 0.4 0.2 0
Рис. 1 - Графики зависимости функций
0 50 100 150 200 250 300 350 400
1
иь )
ин )
(кривая 1),
ис (*)
(кривая 2),
(кривая 3) от * при условии,
и (*) = 0 * < 60 М = 10 М = 1.5
что ь при 1 < 60 для случая ^ / 2
Стационарные точки
Стационарные значения системы уравнений (1) и их устойчивость исследуются с применением численных
методов. Для случая констант
Мз = 20 а = 0.98 V = 0.5 = 0.5
на рис. 2 приведена зависимость
М= 10 ¡и2 = 1.0
, м = 1.0 = 15
кривая 2 - ^ , 72
областью неустойчивых (Л > 0) стационарных точек. В области устойчивых стационарных точек кривые 1 и 2 практически, совпадают.
стационарных значений и от значений параметра Р. Кривая 1 соответствует значениям
. Пунктирная линия отображает границу между областью устойчивых (Л < 0) и
0.7 0.6 0.5 0.4 0.3 0.2 0.1 0
КеЛ>0 1 КеЛ<0
- М1>М2 _
М1<М2 1,2 -
2
0 0.2 0.4 0.6 0.8
Р
и
Р
1.М= 10 М= 1.5
кривая 2
Рис. 2 - Графики зависимости стационарных значений , от параметра . Кривая 1
- М I.0, 1.0. Пунктирная линия отображает границу раздела между областями устойчивых (
Ке Л < 0) и неустойчивых (Л ^ 0) стационарных точек
Р
Стационарные точки системы уравнений (1) в зависимости от значений параметра могут быть как
устойчивыми, так и неустойчивыми. При малых значениях параметра Р в системе возникают колебания. Последнее интерпретируется как возникновение заболевания при недостатке стволовых клеток. На рис. 3 при значениях
9
Р = 0.2 ц = 1.0 ¡л2 = 1.5 !лъ = 20 а = 0.98 v = 0.5 v2 = 0.5
отражено изменение
uh (t)
параметров
и( (г) и (г) п „
( и времени. Стационарная точка в рассматриваемом случае не является устойчивой - в ее окрестности
возникают периодические колебания.
20
40
60
t
uh (t)
80
100
Рис. 3 - Графики зависимости функций
отклонениях от неустойчивой стационарной точки
1Ч U (t) . ЛЧ U, (t) t (кривая 1), s v 7 (кивая 2), ' от 1
от 1 при малых начальных
Модель лечения заболевания
L
Р
Возникновение заболевания в модели (1) рассматривается как изменение в момент времени 0 параметра переводящее устойчивое стационарное состояние в неустойчивое. Тогда под лечением понимается такое воздействие на этот параметр, при котором он возвращается в «нормальное» состояние. При этом введение препарата происходит
и t2 п АР
1 до момента времени 2 Пусть '
от момента времени 1 до момента времени 2. Пусть
Drug (t)
Р
отклонение параметра от исходного значения, а
закон введения препаратов, изменяющих значение возмущения . Тогда модель нарушения функционирования генерации клеток и лечения нарушений с учетом четвертого уравнения в (1) принимает вид
=-/иа (р + р и, )
dt
d АР dt
= -AfiDrug (t)
___ . Drug (t) = const.
Препараты могут подаваться непрерывно во времени ( °4 7 ) или периодически. Модель
периодического воздействия может быть представлена функцией
Drug (t) = D sin(®t) sin(®t) > 0 Drug(t) = 0 если sin>t) < 0
„ u, (t) u, (t)
Результаты моделирования для этой функции представлены на рис. 4 в виде зависимостей h 4 7, d 4 7 и
us (t) .
Для случая констант Р = 0.35, ! =10, ! =15, ! = 20, « = 0.98, v1 = 0.5, v2 = 0 5
стационарная
точка
u = 0.526 ^ = 0.031 и = 0.646 и = 0.359
[0, L 1. . A. _ t = t,
L ' UJ (на рис. 4 зона A ). В момент времени
является устойчивой на временном промежутке 0 параметр Р принимает значение °-20 . Стационарная точка
системы уравнений (1) при этом значении параметра Р будет неустойчивой (рис. 2). В системе возникают колебания (на рис. 4 зона В). То есть система переходит в неустойчивое состояние. С момента времени
t = t
до момента
времени
t = L
вводится препарат с
D = 0.12 и ^ 2^/20 (зона Drug на рис. 4). Постепенное уменьшение
Р
возмущения параметра переводит систему в область устойчивых стационарных точек.
U
0.6 0.5 0.4 0.3 0.2 0.1 0
0
50
100 t
150 t,
200
Рис. 4 - Графики зависимости функций и/ (() (кривая 1), ^ (() (кивая 2), (() (кривая 3) и (()
от
На промежутке
[0, к ] , лл (h, и ]
L 0 J система находится в равновесном состоянии (зона А), на промежутке 4 0 1J
в неустойчивом
состоянии (зона В), на промежутке (tl ' ] (зона ) происходит возврат системы в устойчивое состояние
Решение задачи Коши для систем дифференциальных уравнений осуществлялось с применением метода Рунге-Кутты в модификации Дорманда-Принса в среде программирования математического пакета MatLab. Использовалась
встроенная функция ode45. Параметры систем уравнений ^1, и , определяющие скорости реакций, как следует из результатов клинических исследований, могут отличаться в 2-3 раза. В этом диапазоне отличия параметров системы уравнений не являются жесткими. Тем не менее, поскольку решения строились в окрестности неустойчивых точек, результаты решения сопоставлялись с решениями, построенными с использованием функции ode23tb, предназначенной для решения «жестких» систем. Результаты, поостренными обеими функциями с точностью 10 -6 совпали в исследуемом диапазоне изменения параметров.
Заключение
Таким образом, разработанная модель объясняет механизм взаимодействия нескольких типов клеток как конкуренцию за функциональное пространство с вытеснением из него всех клеток более активными клетками, а нарушение функций кроветворения как переход системы из устойчивой области в неустойчивую.
Список литературы / References
1. Савостьянов Г. А. Возникновение стволовых клеток в развитии многоклеточности и их количественная характеристика / Г. А. Савостьянов // Цитология. - 2016. - Т. 58. - № 8. - С. 577-593.
2. Копнин Б.П. Современные представления о механизмах злокачественного роста: сходства и различия солидных опухолей и лейкозов / Б.П. Копнин // Клиническая онкогематология. - 2012. - № 5(3). - С. 165-185.
3. Васильев А.В. Ниши стволовых клеток и регенеративная медицина / А. В. Васильев, Е. А. Воротеляк, В. В. Терских // Российский физиологический журнал им. И.М. Сеченова. - 2016. - №103(3). - С. 241-261.
4. Fuentes-Gari M. A mathematical model of subpopulation kinetics for the deconvolution of leukaemia heterogeneity / M. Fuentes-Gari, R. Misener, D. Garcia-Munzer and others // Journal of the Royal Society Interface. - 2015. - V. 12. - № 108. 20150276. doi:10.1098/rsif.2015.0276.
5. Berezansky L. Stability and Controllability Issues in Mathematical Modeling of the Intensive Treatment of Leukemia / L. Berezansky, S. Bunimovich-Mendrazitsky, B. Shklyar // Journal of Optimization Theory and Applications. - 2015. - V. 167. - № 1. - P. 326-341. doi: 10.1007/s10957-015-0717-9.
6. Râdulescu I. R. A complex mathematical model with competition in Leukemia with immune response - An optimal control approach / I. R. Râdulescu, D. Cândea, A. Halanay // IFIP Advances in Information and Communication Technology. -2016. - V. 494. - P. 430-441.doi:10.1063/1.4765573
7. Bunimovich-Mendrazitsky C. Mathematical model of pulsed immunotherapy for bladder cancer / C. Bunimovich-Mendrazitsky, H. Byrne , L. Stone // Bulletin of mathematical biology. - 2008. -V.70. -№ 7. - P. 2055-2076. doi:10.1007/s11538-008-9344-z.
8. Clapp G. A review of mathematical models for leukemia and lymphoma / G. Clapp, D. Levy // Drug Discovery Today: Disease Models. - 2015, V. 16. - № 2. - P. 1-6.
9. Kolpak E. P. Neoplasm Morbidity among the Population of Russia / E. P. Kolpak, I. S. Frantsuzova, K. V. Kuvshinova // International journal of advanced biotechnology and research. - 2017. - V.8. - № 3 - P. 2315-2322.
10. Жукова И. В. Математические модели злокачественной опухоли / И. В. Жукова, E. П. Колпак // Вестник Санкт-Петербургского университета. Серия 10: Прикладная математика. Информатика. Процессы управления. - 2014. - № 3. - С. 5-18.
1
Список литературы на английском языке / References in English
1. Savostyanov G. A. Vozniknovenie stvolovyh kletok v razvitii mnogokletochnosti i ih kolichestvennaja harakteristika [The initiation of stem cells in the development of multicellularity and their quantitative characteristics] / G. A. Savostyanov // Citologiya [Cytology]. - 2016. - Vol. 58. - №8. -P. 577-593. [in Russian]
2. Kopnin B. P. Sovremennye predstavlenija o mehanizmah zlokachestvennogo rosta: shodstva i razlichija solidnyh opuholej i lejkozov [Modern views on the mechanisms of malignant growth: the similarities and differences of solid tumors and leukemias] / B. P. Kopnin // Klinicheskaya onkogematologiya [Clinical oncohematology]. - 2012. -№ 5(3). - P. 165-185. [in Russian]
3. Vasil'ev A. V. Nishi stvolovyh kletok i regenerativnaja medicina [Niches of stem cells and regenerative medicine] / A. V. Vasil'ev, E. A. Vorotelyak, V. V. Terskikh // Rossijskij fiziologicheskij zhurnal im. I. M. Sechenova [Russian physiological journal named after I. M. Sechenov]. -2016. -№103(3). - P. 241-261. [in Russian]
4. Fuentes-Gari M. A mathematical model of subpopulation kinetics for the deconvolution of leukaemia heterogeneity / M. Fuentes-Gari, R. Misener, D. Garcia-Munzer and others // Journal of the Royal Society Interface. - 2015. - V. 12. - № 108. 20150276. doi:10.1098/rsif.2015.0276.
5. Berezansky L. Stability and Controllability Issues in Mathematical Modeling of the Intensive Treatment of Leukemia / L. Berezansky, S. Bunimovich-Mendrazitsky, B. Shklyar // Journal of Optimization Theory and Applications. - 2015. - V. 167. - № 1. - P. 326-341. doi: 10.1007/s10957-015-0717-9.
6. Radulescu I. R. A complex mathematical model with competition in Leukemia with immune response - An optimal control approach / I. R. Radulescu, D. Candea, A. Halanay // IFIP Advances in Information and Communication Technology. -2016. - V. 494. - P. 430-441.doi:10.1063/1.47655733.
7. Bunimovich-Mendrazitsky C. Mathematical model of pulsed immunotherapy for bladder cancer / C. Bunimovich-Mendrazitsky, H. Byrne , L. Stone // Bulletin of mathematical biology. - 2008. -V.70. -№ 7. - P. 2055-2076. doi:10.1007/s11538-008-9344-z.
8. Clapp G. A review of mathematical models for leukemia and lymphoma / G. Clapp, D. Levy // Drug Discovery Today: Disease Models. - 2015, V. 16. - № 2. - P. 1-6.
9. Kolpak E. P. Neoplasm Morbidity among the Population of Russia / E. P. Kolpak, I. S. Frantsuzova, K. V. Kuvshinova // International journal of advanced biotechnology and research. - 2017. - V.8. - № 3 - P. 2315-2322.
10. Zhukova I. V. Matematicheskie modeli zlokachestvennoj opuholi [Mathematical models of malignant tumor] / I. V. Zhukova, E. P. Kolpak // Vestnik Sankt-Peterburgskogo universiteta. Serija 10: Prikladnaya matematika. Informatika. Processy upravlenija [Bulletin of Saint-Petersburg university. Series 10: Applied mathematics. Informatics. Control processes]. -2014. -№3. - P. 5-18. [in Russian]
