УДК 538.9; 616-006
https://doi.org/! 0.18500/0869-6632-2020-28-1 -29-45
Нелинейная двумерная динамика взаимодействующих раковых клеток в условиях экстраклеточного поля
О. И. Артемова, В. Д. Кревчик, М. Б. Семенов
Пензенский государственный университет Россия, 440026 Пенза, ул. Красная, 40 E-mail: artyomovaolg@gmail.com, physics@pnzgu.ru, Misha29.02.l@gmail.com Автор для переписки Михаил Борисович Семенов, Misha29.02.l@gmail.com
Поступила в редакцию 21.10.2019, принята к публикации 6.12.2019, опубликована 26.02.2020
Цель настоящей работы - теоретически исследовать влияние экстраклеточного поля на антипараллельную и параллельную подвижность взаимодействующих раковых клеток в зависимости от внешних параметров экстраклеточной матрицы и параметров профиля метаболического потенциала в динамике рисков формирования и развития раковой опухоли. Методы. В данной работе используется метод квазиклассического инстантонного приближения (разреженного «газа» пар «инстантон - антиинстантон») для рассмотрения подвижности взаимодействующих раковых клеток в модельном 2D осцилляторном потенциале в условиях внешнего экстраклеточного поля. Результаты. В рамках построенной аналитической 20-модели показано, что микровзаимодействие через экстраклеточную матрицу возникающих раковых клеток посредством динамического метаболического профиля существенно влияет на динамику рисков формирования и развития раковой опухоли. Показано, что в зависимости от структуры 2D метаболического профиля возникает ряд характерных нелинейных особенностей типа 2D бифуркаций, биений, хаотизации, накладываемых на интегральные динамические кривые, напоминающие по виду функцию Гомпертца, описывающие вероятные риски формирования и развития раковой опухоли. Заключение. Теоретически исследован эффект неустойчивого излома в температурной зависимости антипараллельной динамики двух раковых клеток в 2D метаболическом потенциале. Показано, что эффект нарушения симметрии оказывается устойчивым для параллельного и неустойчивым для антипараллельного переноса. Выявлена сложная тонкая структура в бифуркационной области, обусловленная флуктуациями (биениями) для параллельного двумерного переноса раковых клеток. Для случая антипараллельного переноса в бифуркационной области становится существенным вклад 4, 6, 12, и т.д. пар траекторий, что напоминает один из режимов перехода к хаосу. Дополнительно исследовано взаимодействие пары клеток с колебательными модами экстраклеточной матрицы. Такое взаимодействие значительно и различным образом модифицирует процессы антипараллельного и параллельного двумерного переноса раковых клеток.
Ключевые слова: экстраклеточное поле, метаболический потенциал, бифуркационная и хаотическая динамика, модели формирования и развития раковой опухоли.
Образец цитирования: Артемова О.П., Кревчик В.Д., Семенов М.Б. Нелинейная двумерная динамика взаимодействующих раковых клеток в условиях экстраклеточного поля//Известия вузов. ПНД. 2020. Т. 28, № 1. С. 29-45. https://doi.org/10.18500/0869-6632-2020-28-l-29-45
Статья опубликована на условиях лицензии Creative Commons Attribution License (CC-BY 4.0).
Финансовая поддержка. Работа выполнена в рамках госзадания Минобрнауки РФ, грант 3.6321.2017/8.9 и 2019-0940 (0748-2020-0012).
©Артемова О.И., Кревчик В.Д., Семенов М.Б., 2020
29
https://doi.org/10.18500/0869-6632-2020-28-1-29-45
Nonlinear two-dimensional dynamics of interacting cancer cells under extracellular field conditions
O. I. Artemova, V.D. Krevchick, M. B. Semenov
Penza State University 40, Krasnaya str., Penza 440026, Russia E-mail: artyomovaolg@gmail.com, physics@pnzgu.ru, Misha29.02.1@gmail.com Correspondence should be addressed to Mikhail B. Semenov, Misha29.02.1@gmail.com Received 21.10.2019, accepted 6.12.2019, published 26.02.2020
Purpose of this work is to theoretically investigate the effect of extracellular field on antiparallel and parallel transfer of interacting cancer cells depending on external parameters of extracellular matrix and parameters of metabolic potential profile in the dynamics of cancer tumor formation and development risks. Methods. In this paper, we use the quasi-classical instanton approximation method (or method of the rarefied «gas» of the instanton - antiinstanton pairs) to consider the mobility of interacting cancer cells in the 2D oscillator potential model under conditions of an external extracellular field. Results. Within framework of the proposed analytical 2D-model, it has been shown that microinteraction in the extracellular matrix of emerging cancer cells through the dynamic metabolic profile significantly affects the dynamics of the formation risks and development of a cancer tumor. It is shown that, depending on the structure of the 2D metabolic profile, a number of characteristic nonlinear features arise, such as 2D bifurcations, beats, randomization superimposed on integral dynamic curves resembling in appearance of the Gompertz function that describe probable risks of the formation and development of a cancerous tumor. Conclusion. Effect of an unstable break in the temperature dependence of the antiparallel dynamics of two cancer cells in a 2D metabolic potential has been theoretically investigated. It is shown that effect of the symmetry breaking is stable for parallel and unstable for antiparallel transfer. A complex fine structure was revealed in the bifurcation region due to fluctuations (beats) for parallel two-dimensional transfer of cancer cells. For the case of antiparallel transfer in the bifurcation region, the contribution of 4, 6, 12, etc. pairs of trajectories becomes significant, which resembles one of the possible regimes of transition to chaos. Additionally, the interaction of a pair of cells with the vibrational modes of the extracellular matrix has been investigated. Such interaction modifies significantly and differently the processes of antiparallel and parallel two-dimensional transfers of cancer cells.
Key words: extracellular field, metabolic potential, bifurcation and chaotic dynamics, models of the formation and development of a cancer tumor.
Reference: Artemova O.I., Krevchik V.D., Semenov M.B. Nonlinear two-dimensional dynamics of interacting cancer cells under extracellular field conditions. Izvestiya VUZ. Applied Nonlinear Dynamics, 2020, vol. 28, no. 1, pp. 29-45. https://doi.org/10.18500/0869-6632-2020-28-1-29-45
This is an open access article distributed under the terms of Creative Commons Attribution License (CC-BY 4.0).
Acknowledgements. This work has been performed as part of the state task of the Ministry of Education and Science of the Russian Federation, grant 3.6321.2017/8.9 and 2019-0940 (0748-2020-0012).
Введение
В настоящее время публикуется достаточно много работ, утверждающих, что рак является «метаболическим» заболеванием («metabolic disease») [1-15]. В одной из недавних, активно цитируемых статей исследователей школы медицины из Йельского университета [15] говорится, что «пусковые механизмы» раковых метаболических изменений остаются на сегодняшний день неясными, и предлагаются различные модели для их объяснения. В частности, эти механизмы могут носить стохастический характер. При этом «стохастичность» метаболизма может проявляться как внутри раковых клеток, так и при их взаимодействии с «экстраклеточной матрицей» («extracellular matrix») [1,15] при реализации различных «метаболических сигнальных путей» («metabolic signaling pathways») мобильности раковых клеток. Сотрудники института рака Бостона совместно с исследователями-физиками Массачусетского технологического института недавно выполнили интересные исследования [1], показывающие, что метаболизм раковых
клеток обусловлен не только внутриклеточными онкогенными мутациями, но что на него существенное влияние оказывают различные физические факторы микроокружения раковых клеток или различные физические свойства «экстраклеточной матрицы». В частности, в выращиваемых искусственно 2Э- и ЗЭ-культурах раковых клеток исследовалась с помощью изотопов-маркеров клеточная подвижность как реакция на наличие питательных веществ, локальные изменения в экстраклеточной матрице, изменения температуры и др. Проводилось качественное сравнение с клеточной подвижностью при развитии опухолей в организме мышей. Сделан вывод о существенном влиянии упомянутых физических факторов на метаболизм раковых клеток и клеточную подвижность [1]. Эти исследования вдохновили нас на попытку использования нелинейной микромодели, описывающей подвижность взаимодействующих частиц в матрице-термостате, показавшей свою работоспособность при теоретическом исследовании температурной химической кинетики в 2Э-соединениях типа порфиринов [16]. Хорошо известно, что гем, один из представителей порфиринов, является простетической (небелковой) группой гемоглобина, а некоторые производные порфирина вызывают большой интерес с точки зрения терапии онкологических заболеваний. (Так, димегин, обладая высоким сродством к раковым клеткам, избирательно фотосенсибилизирует их. Благодаря этому последующее воздействие лазерного излучения приводит к редукции опухолевой ткани.) Взаимодействие заряженных частиц (протонов) в порфиринах описывается поверхностью потенциальной энергии (ППЭ), напоминающей 2Э осцилляторный потенциал, используемый и в нашей нелинейной динамической модели.
Таким образом, объектом исследования нашей нелинейной 2Э-модели является подвижность взаимодействующих раковых клеток без учета их внутренней структуры. А предметом исследования является учет влияния внешних физических факторов типа температуры и локальных параметров экстраклеточной матрицы на эту подвижность. Аналогом ППЭ в нашей модели выступает «метаболический потенциал» (МП).
В разделе 1 статьи приводится описание нелинейной динамической 2Э-модели с использованием осцилляторного «метаболического потенциала» с учетом параметров экстраклеточной матрицы. В разделе 2 рассматривается хаотизация подвижности раковых клеток в 2Э-модели антипараллельного переноса в «метаболическом потенциале». В разделе З теоретически исследованы режимы бифуркаций и биений подвижности раковых клеток в 2Э-модели параллельного переноса в МП с учетом двух типов экстраклеточной матрицы. Также сопоставляются факторы риска формирования раковых опухолей с теоретическими параметрами модели. В разделе 4 проводится сравнение полученных результатов аналитических 2Э-моделей антипараллельной и параллельной динамики с результатами моделирования других авторов по рассмотрению хаотической и бифуркационной динамики в системе «опухоль - иммунный кластер - вирус».
1. Описание модели
При построении динамической нелинейной 2Э-модели влияния экстраклеточного поля в динамике рисков формирования и развития раковой опухоли вводится понятие параллельных и антипараллельных координат метаболических путей роста [1,14]. Аналитически вводится 2Э-потенциал экстраклеточного поля для случаев параллельного и антипараллельного 2Э-переноса пары раковых клеток с учетом их взаимодействия через микроокружение. Направление этих координат коррелирует с наличием (или отсутствием в случае антипараллельного переноса) характерных градиентов метаболического потенциала [1], связанных с наличием клеточных пищевых ресурсов и т.д. В качестве одного из существенных внешних параметров, влияющих на динамику топологии 2Э-потенциала экстраклеточного поля, рассматривается фактор взаимодействия со стромальными клетками из-за конкуренции за ограниченные пищевые ресурсы и т.д.
В результате подобного взаимодействия в микроокружение вбрасываются секреты [1], влияющие на поведение стромальных клеток с формированием условий, подавляющих или ускоряющих возможное развитие опухоли. Среди модельных параметров экстраклеточной матрицы, обеспечивающих условия взаимодействия раковых клеток, вводятся колебательные моды экстраклеточной матрицы, обуславливающие ее физические свойства, коэффициенты взаимодействия с этими колебательными модами в линейном приближении, а также параметр квазиравновесной температуры экстраклеточной матрицы, влияющий на динамику рисков формирования и развития раковой опухоли. При моделировании параллельного 2Э-переноса в экстраклеточном поле рассматриваются два типа экстраклеточной матрицы, различающиеся знаком коэффициента взаимодействия раковых клеток, что существенно влияет на динамику рисков возникновения и формирования раковой опухоли.
Рассматриваются особенности 2Э антипараллельной динамики подвижности раковых клеток в экстраклеточном поле (МП) с учетом режимов неустойчивых бифуркаций и хаотизации, влияющих на динамику рисков возникновения и формирования раковой опухоли.
Двухклеточную подвижность в синхронном и асинхронном режимах предлагается рассмотреть в рамках теории 2Э подбарьерного переноса в одноинстантонном приближении с учетом достаточно успешной адаптации к системам типа порфиринов. Использование такой модели оправдано для случая 2Э антипараллельного протонного переноса с учетом неустойчивой точки бифуркации на температурной зависимости для этих соединений [16]. Далее определяется критическая температура Тс, соответствующая бифуркации подбарьерной траектории. Также исследуется эффект влияния локальной моды экстраклеточной матрицы-термостата на вероятность двумерного переноса в МП. При определенных значениях параметров оказывается важным вырождение траекторий антипараллельного переноса. Таким образом, следует учитывать эффект каскада бифуркаций, то есть появления четырех, шести, двенадцати и т.д. пар траекторий. Для случая параллельного переноса клеток бифуркация напоминает фазовый переход первого рода, тогда как для случая антипараллельного клеточного переноса бифуркационное поведение напоминает фазовый переход второго рода. Следует отметить, что предложенная модель позволяет объяснить экспериментальные данные по квантовым флуктуациям в двухпротонном туннелировании в порфиринах вблизи критической температуры, а также предсказать вероятные особенности рисков формирования и развития раковой опухоли в определенных условиях и при определенном профиле МП.
Рассмотрим две клетки, которые осуществляют перенос по метаболическим путям роста в независимых двухъямных потенциалах и(91) и и(92), представляемых как
где сумма а + Ь определяет длину связи в соответствующем фрагменте экстраклеточеной матрицы, Д/ = (1/2)ю2(Ь2 — а2) является смещением (параметром асимметрии метаболического потенциала), в(^) - ступенчатая функция, ю - частота. Масса клетки входит в определение (формально мы полагаем массу равной 1). Взаимодействие между двумя клетками рассматривается в диполь-дипольном приближении
2. Хаотизация подвижности раковых клеток в 2D-модели антипараллельного переноса в экстраклеточном поле
(1)
Vint (91,92)
а
(91 - 92)2,
(2)
где а - положительная константа, а* = (2а)/со2 является безразмерным параметром, а* < 1; и #2 ~ координаты метаболических путей роста клеток.
Для антипараллельного переноса двумерная поверхность метаболического потенциала с учетом взаимодействия клеток может быть определена как
Ua(qi,Q2) = 2Uaiq^q2) = (qi + a)2Q(—qi) + [-(б2 - a2) + (9l - 6)2] 0(9l)+
ОГ
(02 - а)2в(д2) + [-(62 - a2) + (g2 + b)2] 9(~q2) ~\(0i~ 02?. (3)
Потенциал (3) изображен на рис. .
При антипараллельном переносе двух клеток (формула (3) для поверхности метаболического потенциала) инстантонное действие, которое зависит от двух параметров гит (определяющих моменты времени прохождения клетками вершин МП по параллельным координатам метаболических путей роста), для = 0 (пренебрегая взаимодействием с экстраклеточной матрицей) вычисляется с помощью выражения
S=-
a)t(&2 - а2) со (a + Ъ)2 Г / 1
1 - а*
1 - а*
sinh (|е|л/1 ct*)
(1 - а*)3/2
sinh |е| +
cosh (ел/1 ct* + 1)
(1 - а*)3/2
sinh (р*л/Г
-1
cosh
((P*-t)VT
— cosh у 1 — а
cosh 8—1 sinh |3*
[cosh(|3*-т)+ cosh|3*] }. (4)
Рис. 1. Асимметричная поверхность метаболического потенциала (3) для случая антипараллельного клеточного переноса: a = 2, Ъ = 2.3, а* = 0.1; А и В обозначают исходное и конечное состояния взаимодействующих клеток
Fig. 1. Asymmetric surface of metabolic potential (3) for the case of antiparallel transfer: a = 2, b = 2.3, a* = 0.1; A and В denote the initial and final state of interacting cells
Параметры е и т найдем из следующей системы уравнений:
— sinh е [coth Р* + cosh т coth р* — sinh т] + --— sinh ^е/1 — а* j coth ^Р* /1 — а*^ —
— cosh (т\/1 — а*) coth ( р*\/ 1 — аЧ + sinh (т\/1 —
1
(1 + b*)(1 — а*) 1 — а
а*
+ -—+ (cosh е — 1)(sinh т coth р* — cosh т) + cosh е + 1
1 а*
cosh ( е\/1 — а*) + 1 sinh (ту/1 — а*) coth ( р*\/ 1 — а*) — cosh ( ту/ 1 —
— cosh (W 1 — а^ | =0. (5)
Для юр ^ 1 (предел сравнительно низких температур экстраклеточной матрицы), решение системы уравнений (5) может быть найдено с использованием теории возмущений (для малых параметров е) при заданных величинах параметров асимметрии МП (Ь — а)/(Ь + а) и а*. Для е = 0 действие ( ) примет вид
rn(b2 — а2) . S = —-Т^тг Arccosh
(1 — а*)3/2
coth
b — а , шРл/ 1 — а*
--sinh-
b + а 2
а
2
ю2р(а + b)2 m(b + а)2
2(1 — а2) + (1 — а*)3/2 Х
\ 2\ I/2'
— аа ( b — а
2 +
/ -2 юРУТ^СХ* + (b-ol) V 2 + (b + а)2 J
(6)
Так же, как и для случая параллельного переноса, при р > рс (при некоторой критической температуре) происходит смена режимов клеточной подвижности: с синхронного переноса (с одной экстремальной траекторией) на асинхронный (которому отвечают две отщепленные траектории). Однако, в отличие от параллельного переноса, данное отщепление происходит для любых значений параметров метаболического потенциала. При условии, что Р > Р с^ имеем 5*е=0 > 5*е=0. Тогда величина £е=о будет определять скорость переноса (или вероятного риска перспективного развития опухоли). А для р < рс будем иметь преобразование двух вырожденных отщепленных траекторий в единичную, ql = —52, которая соответствует синхронному антипараллельному переносу клеток.
Для случая одноклеточного переноса существует только одна траектория, которая минимизирует действие. Для переноса двух клеток, взаимодействующих друг с другом, существует два типа подбарьерных траекторий. А точнее, в зависимости от величины параметра Р, основной вклад в инстантонное действие определяется либо двукратно вырожденной, либо единичной траекторией. Также отметим, что при параллельном переносе, для р > рс, перенос оказывается асинхронным, Т1 = Т2. То есть клетки проходят верхушки барьера не одновременно. Расщепление единичной клеточной траектории (51 = 52) не происходит для малых значений параметров взаимодействия а* и для температур экстраклеточной матрицы, которые удовлетворяют условию Р < рс .То есть клетки проходят верхушки барьера по своим координатам метаболических путей роста МП в один момент времени (т1 = Т2), и перенос оказывается синхронным.
Тип взаимодействия, задаваемый соотношениями (2), (3), оказывается таким, что не влияет на движение вдоль координаты центра масс клеток, 51 = 52. По этой причине евклидово действие оказывается не зависимым от параметра взаимодействия, как для случая параллельного переноса. Численный анализ системы трансцендентных уравнений (5) выявил интересные особенности для переходной области между режимами переноса клеток, то есть тонкую структуру
0
4
х
*
X
X
(Tj+T^OÖ ß*= 18.83
- ::::::::::::::::::т Т.:::::::::::::::::
- ::::::т:
: 1 1 1 1 1 1 1 (Xj -т2)ю 1 1 1 1 1
-5 -3 0 3 5
с
Рис. 2. Численное решение трансцендентных уравнений (5), отвечающее режиму каскада бифуркаций или хаотизации
Fig. 2. Numerical solution of transcendental equations (5) corresponding to the regime of a cascade of bifurcations or randomization
вблизи первой бифуркационной точки для антипараллельного переноса. Результаты численного анализа представлены на рис. 2. Установлено, что в добавление к первой бифуркационной точке, характеризуемой двумя решениями (рис. 2, а), существуют дополнительные бифуркационные точки при более «низких» температурах экстраклеточной матрицы, то есть, например, 4 пары (рис. 2, b), 6 пар (рис. 2, с) и даже 12 пар дополнительных решений при ß* = 19.2009 (а* = 0.05) и т.д. Это явление похоже на множественные бифуркации или каскад бифуркаций. Такой эффект напоминает один из сценариев перехода к хаосу.
Хотя синхронный режим оказывается предпочтительным благодаря тому, что в этом случае действие принимает минимальное значение, в определенном диапазоне температур экстраклеточной матрицы эта величина действия оказывается сравнимой с теми значениями, что соответствуют каскадным решениям. В результате происходят флуктуации или биения нерегулярного характера в противоположность тому, что наблюдается для параллельного переноса. Антипараллельный 2D клеточный перенос, таким образом, характеризуется неустойчивостью перехода от синхронного к асинхронному поведению. Такие неустойчивости напоминают непрерывный фазовый переход второго рода, тогда как при параллельном переносе раковых клеток бифуркационный режим напоминает ступенчатый процесс типа фазового перехода первого рода. Найдено, что зависимости для ßc(a) и ac(ß) в случае антипараллельного переноса имеют такой же характер, как и в случае параллельного переноса.
Таким образом, выявлена достаточно сложная тонкая структура перехода для параллельного и антипараллельного переноса двух клеток с различными вырожденными подбарьерными траекториями, приводящими к биениям, каскадам бифуркаций.
Эффекты нарушения симметрии могут иметь место при сравнительно высоких температурах экстраклеточной матрицы, зависящих от частоты барьера МП. Например, для порфиринов критическая температура Тс составляет 200 К. Для раковых клеток она окажется заметно выше (до 310 К).
3. Режимы бифуркаций и биений подвижности раковых клеток в 2В-модели параллельного переноса в экстраклеточном поле с учетом двух типов экстраклеточной матрицы
Рассматриваются особенности 2Б параллельной динамики в экстраклеточном поле (МП) с учетом режимов устойчивых бифуркаций и биений двух типов, влияющих на динамику рисков возникновения и формирования раковой опухоли. Как и в предыдущем разделе рассмотрим подвижность двух раковых клеток, которые осуществляют перенос по параллельным координатам метаболических путей роста МП и клеток.
Двумерная поверхность МП для случая параллельного клеточного переноса, нормированная на со2, задается соотношением (рис. 3)
иМъ 02) = 2-Щ,92) = (51 + г) + [-(б2 - а2) + (91 - б)2] 9Ы+
О)
+ + а2)9(-д2) + [-(б2 - а2) + (д2 - Ъ)2] в(®) - у (дг - д2)2. (7)
Рис. 3. Асимметричная поверхность МП (7) для случая параллельного переноса (подвижности) раковых клеток: А и В обозначают исходное и конечное состояния клеток, соответственно
Fig. 3. Asymmetric surface of the metabolic potential (7) for the case of parallel transfer (motility) of cancer cells: A and В denote the initial and final state of the cells, respectively
Параметры «а» и «Ь» потенциала перенормируются во внешнем МП поле: а = а* — А,/®о> Ь = Ь* + (где X - параметр, определяющий напряженность МП, связанную с особенностями взаимодействия раковых и стромальных клеток); либо перенормируется безразмерный параметр Ь* = Ь/а, который слабо нелинейно зависит от напряженности МП.
Квазиклассическое (инстантонное) действие, которое с экспоненциальной точностью определяет вероятность 2Э параллельного клеточного переноса, рассчитывается по формуле
5 = 2а(а + Ь)("1 + хо)ш2 — 1 ш2(а + Ь)2(т1 + то)2 — ^ + Ь)°(Т1 —
в v у v 2 (ю2 - 2а)в
2ю4(а + b)2 ^л |"sin2(vraXl) + sin2(v„X2) + (sin(v„Xl) - sin(v„X2))2
E
в n=1
(vn + ®2 + U) (v£ + ю2 - 2a)
, (8)
где были введены следующие обозначения: е = е*ш = (т1 — т2)ю; т = 2т*ю = (т1 + т2)ю; в* = (Рю)/2; а* = (2а)/ш2; Ь* = Ь/а, т1 и т2 - центры 2Б-инстантона, то есть моменты мнимого времени проскока параллельно перемещающимися клетками в поле МП верхушки потенциального барьера МП (точки с нулевой координатой) вдоль соответствующей координаты метаболических путей роста; уга = (2пп)/в - мацубаровские частоты; в - обратная температура экстраклеточной матрицы.
В пределе слабой диссипации, когда взаимодействие с локальными колебательными модами экстраклеточной матрицы предполагается пренебрежимо малым, получим следующие результаты для вероятности 2Э параллельного клеточного переноса. В случае пренебрежения взаимодействием с экстраклеточной матрицей, то есть при = 0, действие ( ) как функция параметров е и т принимает вид
S =
(а + Ь)2ю
4аТ Т Л , 1 А , (Т -|8|)а* , „^u-lo*
1 + -- + v ' 17 + coth в* - sinh-1 в* х
а + b а + b\ 1 — а*/ 1 — а* х [cosh(e* - т) cosh е + cosh(e* - т) - cosh(e* - |е|)| - (1 - а)-3/2 (- coth (вVT—а*) +
+ sinh-1 (в^1-а*) {cosh ((в*-т) /1-а*) [cosh (е/1-а*-1)] +cosh ((в*-|е|)/1-а*)})
(9)
Как только траектория найдена, уравнения q1(T1) = 0, q2(T2) = 0 могут быть представлены в следующей форме:
sinh е [coshт coth в* - sinhт - coth в*] + --* sinh (е/1 - а*) х
х [cosh (ту/1 - а*) coth (в^л/Т-а*) - sinh (т\/1 - а*) + coth (в* уТ-а*)] = 0, 41
3--;----+ cosh е [sinh т coth в* — cosh т — 1] + sinh т coth в* — cosh т+
1 + b* 1 - а* L 1 J 1
+ -—^^ cosh (е/1 - а*) [sinh (т/1 - а*) coth (в*л/1 - а*) - cosh (т/1 - а*) +1] -1 — а
[sinh (т\/1 - а*) coth (в*/1 - а*) - cosh (ту/1 - а*)] = 0. (10)
1 — а
Численное решение системы (10) позволяет выявить бифуркацию 2Э клеточных траекторий, то есть при определенном значении обратной температуры в* экстраклеточной матрицы, либо параметра асимметрии потенциала МП, связанного с величиной напряженности приложенного поля МП Ь* = Ь/а, либо коэффициента взаимодействия а* = (2а)/ю2. (Знак коэффициента
взаимодеиствия раковых клеток зависит, в частности, от типа и параметров экстраклеточнои матрицы.) Численный анализ системы (10) позволяет также выявить тонкую структуру перехода в окрестности точки бифуркации, а именно режим биений для параллельного переноса (подвижности) раковых клеток. В итоге вероятность 2D параллельного клеточного переноса с экспоненциальной точностью определяется как Г = ехр(—5), где 5 задается выражением ( ), с учетом решения системы (10).
Простые аналитические решения уравнений (10) получаются в частном случае, когда
ю
е = (ti — т2)ю = 0, Vß, а < —,
ti = Т2 = = -1 Arccosh 2ю 2ю
1 — Ь*(Х) , ßrn
-——sinh—
1 + Ь*(Х) 2
+ 4
(11)
Однако полный анализ требует громоздкого численного расчета. При достаточно низких температурах, юр » 1, при 1 < Ь/а < 3, и
Ь — a 2а * 2(Ь — а)
< — < а* = v 7
2(Ь + а) - ю2 с 3Ь - а мы, наконец, получим решение с экспоненциальной точностью
exp (—т\/ 1 — а*) ~
3 -
4
1
1 + Ь*
1 а*
(1 — ^lAW1-^)
х^ 1 + (1 — а*)
*) 1/(1-^/1-0^)
1 а*
+ з —
1 + Ь*
1а
.у (1 — VT-«'
1
exp(—е) ~
3
1 + Ь*
1 а*
exp (т\/1 — а*) +
1 а*
(12)
Решение (12) справедливо при
ß > ßc =
С^Г—
а*
ю
(13)
Согласно предлагаемой модели оценивается пороговая вероятность формирования раковой опухоли выше 0.7 (70%), см. Таблицу. При реализации критических значений параметров
Таблица/Table
х
1
4
1
4
1
1
Факторы рисков формирования раковых опухолей Теоретические параметры Теоретическая зависимость
Повышенная субфебрильная температура (выше 37° C) Р* ~ 1/Т, р* > 4 Критическая температура, определяющая точку бифуркации, с последующим лавинообразным ростом подвижности раковых клеток
Уровень каспаз (3, 6, 8), увеличение активности на 25 % а* > 0.3 (1 > а* > 0) Параметр, определяющий интенсивность взаимодействия раковых клеток
Кольпоскопический индекс Рейда (выше 3) Хс («напряженность» метаболического потенциала), Ь*(Х)>1.5 Параметр асимметрии метаболического потенциала Ь* (X)
|3С, а*, \с осуществляется переход через точку бифуркации с последующим лавинообразным ростом подвижности раковых клеток с вероятностью 0.7 (70%).
Полученная зависимость вероятности 2Б параллельного клеточного переноса от напряженности МП-поля с учетом влияния двух локальных колебательных мод экстраклеточной матрицы позволяет проанализировать режим 2Б клеточных бифуркаций (смена режима переноса - клеточной подвижности - с синхронного на асинхронный), а также биений в окрестности точки бифуркации. Так, на рис. 4 после режима синхронного параллельного клеточного переноса с двумя характерными пиками точка излома отвечает точке бифуркации, а последующие осцилляции -биениям.
Наряду с режимом биений с «провалами» на зависимости вероятности 2Б параллельного клеточного переноса от напряженности МП-поля (см. рис. ) при увеличении температуры экстраклеточной матрицы и при частотах колебательных локальных мод, значительно меньших характерных частот двухъямного метаболического потенциала, может иметь место режим биений с резонансной структурой (рис. 5), связанной с интерференцией различных каналов клеточного переноса.
В следующем разделе представлены результаты модели параллельного 2D клеточного переноса с учетом смены знака взаимодействия клеток в экстраклеточной матрице при изменении условий клеточной конкуренции за источники питания, а также проведено сравнение результатов представленного модельного рассмотрения с результатами моделирования других авторов.
3.5 Ь(к)
Рис. 4. Зависимость вероятности 2D параллельного клеточного переноса от напряженности МП-поля с учетом точки бифуркации и режима биений
Fig. 4. Dependence of the 2D parallel cellular transfer probability on the strength of the MP-field, taking into account the bifurcation point and the beat mode
Рис. 5. Режим биений на зависимости вероятности 2D параллельного клеточного переноса от напряженности МП-поля: резонансы до точки бифуркации; резонансы и провалы выше точки бифуркации
Fig. 5. Beat mode on the 2D parallel cellular transfer probability dependence on the MP-field strength: resonances to the bifurcation point; resonances and dips above the bifurcation point
4. Сравнение результатов аналитических моделей 2В- антипараллельной и параллельной динамики (разделы 2 и 3) с результатами моделирования других авторов по рассмотрению хаотической и бифуркационной динамики в системе «опухоль - иммунный кластер - вирус»
В рамках модели, рассмотренной в разделе 3, в режиме слабого взаимодействия раковых клеток с экстраклеточной матрицей в левой части кривой зависимости вероятности 2Б параллельного переноса (клеточной подвижности) от напряженности МП наблюдается единичный пик (рис. 6), отвечающий случаю, когда МП при некотором значении параметра напряженности МП (X) становится симметричным. Излом в правой части модельной кривой отвечает устойчивой точке бифуркации при 2Т> параллельном переносе в режиме слабой диссипации: или пренебрежении взаимодействием раковых клеток с экстраклеточной матрицей.
На рис. представлены результаты моделирования параллельного 2Т> клеточного переноса с учетом смены знака взаимодействия клеток в экстраклеточной матрице при изменении условий клеточной конкуренции за источники питания. На рисунке рассмотрен режим двойной бифуркации (на модельной кривой точки бифуркации обозначены стрелками 1).
На рис. 8 представлены эффекты биений в окрестностях точек двойной бифуркации в случае смены знака взаимодействия клеток в экстраклеточной матрице (обозначены стрелками 2).
В пределе сильного взаимодействия раковых клеток с экстраклеточной матрицей (раздел 3) с учетом влияния на их подвижность двух локальных колебательных мод экстраклеточной матрицы получим модельную зависимость вероятности 2Б параллельного клеточного переноса с учетом режима бифуркации и биений, представленную на рис. 4 и 5. Правый фрагмент этой модельной зависимости (правее точки бифуркации) частично напоминает результат моделирования временной динамики плотности раковой опухоли с учетом режима бифуркаций (см. рис. 9 и 10 в сравнении с результатами работы [2]).
Г, а.и. 1.0
0.1
о
10 Ь(Х), а.и.
1
2 Ь(к), а.и.
Рис. 6. Зависимость вероятности 2D параллельного переноса (клеточной подвижности) от напряженности МП в пределе слабого взаимодействия раковых клеток с экстраклеточной матрицей
Fig. 6. Dependence of the 2D parallel transfer (cell motility) probability on the strength of the MP-field in the limit of weak interaction of cancer cells with extracellular matrix
Рис. 7. Режим двойной бифуркации (на модельной кривой точки бифуркаций обозначены стрелками 1) на зависимости вероятности параллельного 2D клеточного переноса с учетом смены знака взаимодействия клеток в экстраклеточной матрице
Fig. 7. Double bifurcation mode (bifurcation points on the model curve are indicated by arrows 1) on the parallel 2D cell transfer probability dependence taking into account the change in the sign of cell interaction in the extracellular matrix
Рис. 8. Режим биений (обозначены стрелками 2) на фоне двойной бифуркации (на модельной кривой точки бифуркаций обозначены стрелками 1) на зависимости вероятности параллельного 2D клеточного переноса с учетом смены знака взаимодействия клеток в экстраклеточной матрице
Fig. 8. Beat modes (indicated by arrows 2) against the background of double bifurcation (on the model curve, bifurcation points are indicated by arrows 1) on the parallel 2D cell transfer probability dependence taking into account the change in the sign of cell interaction in the extracellular matrix
Рис. 9. Сравнение полевой зависимости вероятности 2D параллельного клеточного переноса с учетом режима бифуркации и биений в случае более сильного взаимодействия раковых клеток с экстраклеточной матрицей с учетом влияния на их подвижность двух локальных колебательных мод экстраклеточной матрицы (кривая 1, см. также рис. 4) с результатом моделирования временной динамики плотности раковой опухоли с учетом режима бифуркаций (кривая 2) [2]
Fig. 9. Comparison of the field dependence of the 2D parallel cell transfer probability taking into account the bifurcation and beating regimes in the case of stronger interaction of cancer cells with the extracellular matrix, taking into account the effect of two local vibrational modes of the extracellular matrix on their mobility (curve 1, see also Fig. 4) with the result modeling the temporal dynamics of the cancerous tumor density, taking into account the bifurcation mode curve 2 [2]
1.0-
0.1-
о
4 5
10
b(X), a.u.
Рис. 10. Сравнение режима биений на зависимости вероятности 2D параллельного клеточного переноса от напряженности МП-поля (рис. 5; кривая 1) с модельной кривой динамики роста раковой опухоли по результатам работы [2], кривая 2
Fig. 10. Comparison of the beat mode on the 2D parallel cell transfer probability dependence on the strength of MP-field (Fig. 5; curve 1) with model curve of the cancer growth dynamics according to results of [2], curve 2
Хаотизация клеточной подвижности, представленная в работе [2], напоминает режим антипараллельной подвижности раковых клеток в пределе слабой диссипации, представленной в разделе 2.
Теоретически исследован эффект неустойчивого излома в температурной зависимости антипараллельной динамики двух раковых клеток в 2Т> метаболическом потенциале. Показано, что эффект нарушения симметрии оказывается устойчивым для параллельного и неустойчивым для антипараллельного переноса. Выявлена сложная тонкая структура в бифуркационной области, обусловленная флуктуациями (биениями) для параллельного двумерного переноса раковых клеток. Для случая антипараллельного переноса в бифуркационной области становится существенным вклад четырёх, шести, двенадцати и т.д. пар траекторий, что напоминает один из режимов перехода к хаосу. Дополнительно исследовано взаимодействие пары клеток с колебательными модами экстраклеточной матрицы. Такое взаимодействие значительно и различным образом модифицирует процессы антипараллельного и параллельного двумерного переносов раковых клеток. Взаимодействие с экстраклеточной матрицей не меняет динамику центра масс пары раковых клеток для случая антипараллельного переноса, тогда как в динамику для случая параллельного переноса вносит существенный вклад. Определена граница диапазона температур экстраклеточной матрицы, начиная с которой реализуются устойчивые двумерные синхронные корреляции всех рассмотренных типов. Найдено, что переход от синхронного к асинхронному режиму 2Б антипараллельного переноса раковых клеток в экстраклеточной матрице в условиях 2Т> метаболического потенциала напоминает фазовый переход второго рода.
Показано, что на риски возникновения и формирования раковой опухоли существенное влияние оказывает взаимодействие с экстраклеточной матрицей, то есть реализуется предел
Заключение
сильной диссипации, что совпадает с выводом работ [1, 14]. В рамках проведенного моделирования установлено, что взаимодействие раковых и стромальных клеток, существенно перенормирующее МП и способное менять знак взаимодействия раковых клеток между собой, при определенном критическом значении величины этого взаимодействия приводит к режиму бифуркаций подвижности раковых клеток и, как следствие, к лавинообразному росту раковой опухоли. Показано, что параметры МП и характеристики экстраклеточной матрицы (в том числе физические, что соответствует выводам работы [1]) существенно определяют как риски возникновения, так и характерную динамику роста раковой опухоли.
В результате проведенных количественных оценок в рамках предложенной теоретической модели сформулирован способ оценки рисков развития злокачественных новообразований, в котором в качестве факторов риска предложено выделять субфебрильную температуру, уровень каспаз, кольпоскопический индекс Рейда, порог вероятности формирования злокачественных новообразований. Рассчитана вероятность подвижности взаимодействующих раковых клеток в экстраклеточной матрице в метаболическом потенциале, зависящем от параметров, связанных с определяемыми факторами риска, при значениях которых: кольпоскопический индекс Рейда от 3 и более, уровень активности каспаз от 25 % и более, субфебрильная температура выше 37°C, порог вероятности формирования злокачественных новообразований более 70%, параметр взаимодействия раковых клеток между собой от 0.3 и более, когда может быть диагностирован ускоренный рост злокачественных новообразований.
Библиографический список
1. Muir A., Danai L.V. and Matthew G. Vander Heiden. Microenvironmental regulation of cancer cell metabolism: implications for experimental design and translational studies // Disease Models & Mechanisms. 2018. 11. dmm035758. doi:10.1242/dmm.035758 (pp. 1-12), Special collection: Cancer metabolism.
2. Eftimie R., Macnamara C.K., Dushoff J., Bramson J.L., David J.D. Earn bifurcations and chaotic dynamics in a tumour-immune-virus system // Math. Model. Nat. Phenom. 2016. Vol. 10, no. 14. P. 1-21.
3. Al-Mahdi A.M., Mustafa Q. Khirallah bifurcation analysis of a model of cancer // European Scientific Journal. January 2016 edition. Vol. 12, no. 3. P. 67-83.
4. Fadaka A., Ajiboye B., Ojo O., Adewale O., Olayide I., Emuowhochere R. Biology of glucose metabolization in cancer cells // Journal of Oncological Sciences. 2017. 3. P. 45-51.
5. Колобов А.В., Анашкина А.А., Губернов В.В., Полежаев А.А. Математическая модель роста опухоли с учетом дихотомии миграции и пролиферации // Компьютерные исследования и моделирование. 2009. Т. 1, № 4. P. 415-422.
6. Жукова И.В., Колпак ^.П.Математические модели злокачественной опухоли // Вестник Санкт-Петербургского университета, сер. 10. 2014. Вып. 3. С. 5-17.
7. Kolobov A.V., Gubernov V.V. and Polezhaev A.A. Autowaves in the model of infiltrative tumour growth with migration-proliferation dichotomy // Math. Model. Nat. Phenom. 2011. Vol. 6, no. 7. P. 27-38.
8. Altrock P.M., Liu L.L. and Michor F. The mathematics of cancer: Integrating quantitative models //Nature Reviews: Cancer. 2015. Vol. 15. P. 730-745.
9. Уханова Ю.Ю. Клинико-диагностическое значение маркеров биологических жидкостей в оценке риска пролиферации миомы матки // Диссертация на соискание ученой степени кандидата медицинских наук, Казань, 2016. 146 с.
10. Кучумов А.Г. Математическое моделирование и биомеханический подход к описанию развития, диагностики и лечения онкологических заболеваний // Российский журнал биомеханики. 2010. Т. 14, № 4 (50). C. 42-69.
11. Anastasiou D. Tumour microenvironment factors shaping the cancer metabolism landscape // British Journal of Cancer. 2017. Vol. 116. P. 277-286.
12. Колобов А.В. Моделирование роста и прогрессии опухоли с учетом ее пролифератив-ной и пространственной гетерогенности // Автореферат диссертации на соискание ученой степени кандидата физико-математических наук по специальности 03.00.02 «биофизика», Москва, ФиРАН, 2003.
13. Гомцян Д. (04.07.2017), Математическое исследование рака http://iteach.vspu.ru/07-2017/14329/
14. Boroughs L.K., DeBerardinis R.J.Metabolic pathways promoting cancer cell survival and growth // Nat. Cell Biol. 2015. Vol. 17, no. 4. P. 351-359.
15. Nagarajan A., Malvi P., Wajapeyee N. Oncogene-directed alterations in cancer cell metabolism // Trends Cancer. 2016. Vol. 2, no. 7. P. 365-377.
16. Леггет Э.Дж., Овчинников Ю.Н., Кревчик В.Д., Семенов М.Б. и др. Управляемое дисси-пативное туннелирование. Туннельный транспорт в низкоразмерных системах // Коллективная монография, посвященная памяти академика РАН А.И. Ларкина / Под редакцией Нобелевского лауреата по физике 2003 Э. Дж. Леггета. М.: Физматлит, 2011-2012. 600 с.
References
1. Muir A., Danai L.V. and Matthew G. Vander Heiden. Microenvironmental regulation of cancer cell metabolism: Implications for experimental design and translational studies. Disease Models & Mechanisms, 2018, 11, dmm035758. doi:10.1242/dmm.035758 (pp. 1-12), Special collection: Cancer metabolism.
2. Eftimie R., Macnamara C.K., Dushoff J., Bramson J.L., Earn D.J.D. Bifurcations and chaotic dynamics in a tumour-immune-virus system. Math. Model. Nat. Phenom., 2016, vol. 10, no. 14, pp. 1-21.
3. Al-Mahdi A.M., Khirallah M.Q. Bifurcation analysis of a model of cancer. European Scientific Journal, 2016, vol. 12, no. 3, pp. 67-83.
4. Fadaka A., Ajiboye B., Ojo O., Adewale O., Olayide I., Emuowhochere R. Biology of glucose metabolization in cancer cells. Journal of Oncological Sciences, 2017, 3, pp. 45-51.
5. Kolobov A.V., Anashkina A.A., Gubernov V.V., Polezhaev A.A. Mathematical model of tumor growth with migration and proliferation dichotomy. Computer Research and Modeling, 2009, vol. 1, no. 4, pp. 415-422 (in Russian).
6. Zhukova I.V., Kolpak E.P. Mathematical models of a malignant tumor. Bulletin of St. Petersburg University, ser. 10, 2014, no. 3, pp. 5-17.
7. Kolobov A.V., Gubernov V.V. and Polezhaev A.A. Autowaves in the model of infiltrative tumour growth with migration-proliferation dichotomy. Math. Model. Nat. Phenom., 2011, vol. 6, no. 7, pp. 27-38.
8. Altrock Ph.M., Liu L.L. and Michor F. The mathematics of cancer: Integrating quantitative models. Nature Reviews: Cancer, 2015, vol. 15, pp. 730-745.
9. Ukhanova Yu.Yu. Kliniko-Diagnostic Value of Markers for Biological Fluids in the Risk Assessment of Proliferation of Uterine Fibroids. Thesis for the degree of candidate of medical sciences, Kazan, 2016, 146 p. (in Russian).
10. Kuchumov A.G. Mathematical modeling and biomechanical approach to description of the development, diagnosis and treatment of cancer. Russian Journal of Biomechanics, 2010, vol. 14, no. 4 (50), pp. 42-69 (in Russian).
11. Anastasiou D. Tumour microenvironment factors shaping the cancer metabolism landscape. British Journal of Cancer, 2017, vol. 116, pp. 277-286.
12. Kolobov A.V. Modeling of the Growth and Progression of a Tumor Taking into Account its Proliferative and Spatial Heterogeneity. Abstract of dissertation for the degree of candidate of physical and mathematical sciences, specialty 03.00.02 «Biophysics», Moscow, FiRAN, 2003 (in Russian).
13. Gomtsyan D. (04.07.2017), The Mathematical Study of Cancer, http://iteach.vspu.ru/07-2017/14329/
14. Boroughs L.K., Ralph J. DeBerardinis Metabolic pathways promoting cancer cell survival and growth. Nat. Cell Biol., 2015, vol. 17, no. 4, pp. 351-359.
15. Nagarajan A., Malvi P., Wajapeyee N. Oncogene-directed alterations in cancer cell metabolism. Trends Cancer, 2016, vol. 2, no. 7, pp. 365-377.
16. Leggett A.J., Ovchinnikov Yu.N., Krevchik V.D., Semenov M.B. et al. Controllable dissipative tunneling. Tunnel transport in low-dimensional systems. collective monograph dedicated to the memory of Academician A.I. Larkin, edited by Nobel Laureate in Physics (2003) A.J. Leggett, M. Fizmatlit, 2011-2012. 600 p. (in Russian).
Артемова Ольга Игоревна - родилась в Пензе (1990). Окончила с красным дипломом Медицинский институт Пензенского государственного университета по специальности «Лечебное дело» (2013). Окончила ординатуру по специальности «Акушерство и гинекология» (2015), получила сертификат по направлению «Ультразвуковая диагностика». С 2015 года работает в должности врача акушера-гинеколога в ГБУЗ «Пензенский городской родильный дом». С 2019 года является соискателем на степень кандидата медицинских наук в ФГБОУ ВО «Самарский государственный медицинский университет» на кафедре «Акушерство и гинекология». Автор более 20 научных публикаций и участник международных конференций. Сфера научных интересов - неивазивная диагностика в акушерстве и гинекологии, патологии шейки матки.
Пензенский государственный университет Россия, 440026 Пенза, ул. Красная, 40 E-mail: artyomovaolg@gmail.com
Владимир Дмитриевич Кревчик - родился в Борисове (Беларусь, 1955). Окончил физический факультет Ташкентского государственного университета им. В.И. Ленина (1979). Защитил диссертацию на соискание ученой степени кандидата физико-математических наук на тему «К теории оптических явлений в массивных полупроводниках и размерно-квантованных системах с глубокими примесными центрами» (1983) и диссертацию на соискание ученой степени доктора физико-математических наук на тему «Акустостиму-лированные явления при релаксации фононных возбуждений в полупроводниках» (1994, Физико-технический институт им. С.В. Стародубцева, Ташкент). В 1999 году присвоена ученая степень доктора физико-математических наук на основании переаттестации (ВАК РФ, Москва). В 2000 году присвоено ученое звание профессора, а в 2009 году - звание «Заслуженный деятель науки РФ». Автор более 500 научных работ, из них 6 монографий, в области полупроводниковой наноэлектроники. Основатель и руководитель научно-педагогической школы «Физика электронных процессов в низкоразмерных системах». Подготовил 26 кандидатов и 2 докторов физико-математических наук.
Пензенский государственный университет Россия, 440026 Пенза, ул. Красная, 40 E-mail: physics@pnzgu.ru
Михаил Борисович Семенов - родился в 1960 году в Пензенской области. Окончил Московский государственный университет им. М.В. Ломоносова по специальности «Физика» (1985). До 1995 года работал на кафедре теоретической физики физического факультета МГУ в должности младшего научного сотрудника. Защитил диссертацию на соискание ученой степени кандидата физико-математических наук на тему «Квантовое тун-нелирование с диссипацией» (1991) и диссертацию на соискание ученой степени доктора физико-математических наукна тему «Квантовые механизмы управления параметрами ме-зоскопических систем» (2004). C 1995 года работает на кафедре «Физика» Пензенского государственного университета, с 2004 года - в должности профессора, c 2013 года - в должности заведующего кафедрой. Автор более 300 научных работ.
Пензенский государственный университет Россия, 440026 Пенза, ул. Красная, 40 E-mail: Misha29.02.1@gmail.com