УДК 535.341:535.36:535.361
DOI: 10.21779/2542-0321-2018-33-4-41-50
12 12 12 1 А.А. Муртазаева ' , К.М. Гираев ' , М.А. Магомедов ' , Н.А. Ашурбеков
Влияние злокачественных новообразований в тканях поджелудочной железы на оптические спектры пищеварительных ферментов
1 Дагестанский государственный университет; Россия, 367001, г. Махачкала, ул. М. Гаджиева, 43а;
2 Институт физики им. Х.И. Амирханова ДНЦ РАН; Россия, 367003, г. Махачкала, ул. М. Ярагского, 94; nashurb@mail.ru
Выполнены спектральные трансмиссионные и рефлектометрические измерения в интервале длин волн 250-2500 нм для проб желчи и панкреатического сока при различных заболеваниях поджелудочной железы. Определены спектры коэффициентов оптического поглощения, рассеяния и фактора анизотропии рассеяния с использованием инверсного метода Монте-Карло. Установлено, что исследуемые объекты представляют собой коллоиды среды, образованные как малыми рассеивателями, так и их крупными фракциями, сформированными в процессе окислительной агрегации биосред. Установлено, что по мере малигнизации для проб пищеварительных ферментов наблюдаются изменения в размерах и концентрации рассеивателей, что следует из динамики спектральной зависимости оптических коэффициентов.
Ключевые слова: желчь, панкреатический сок, процессы малигнизации, спектрофото-метрия, поглощение, рассеяние, фактор анизотропии, морфология, диагностика.
Введение
Известно, что аденокарцинома поджелудочной железы является одним из наиболее распространенных и агрессивных типов онкологических заболеваний [1, 2]. Данная ситуация вызвана множеством причин. Авторами работы [1] отмечено: «Из общего числа негативных факторов можно выделить длительную бессимптомность течения заболевания, сложности в проведении дифференциальной диагностики, высокую чувствительность к употреблению алкоголя, некачественному, нерегулярному и несбалансированному питанию, аутоимунные процессы. При этом, как правило, наблюдается вовлечение в процесс воспаления прилежащих органов с развитием тяжелых форм ассоциированной патологии». В этой связи развитие и совершенствование методов диагностики и контроля заболеваний поджелудочной железы, включая злокачественные новообразования, является важной задачей биофотоники.
На сегодняшний день методы спектрофотометрии интегрирующих сфер с последующим определением коэффициентов поглощения (да), рассеяния (p.s) и фактора анизотропии рассеяния (g) получили широкое распространение в диагностических исследованиях биологических объектов, поскольку позволяют оценить комплекс их морфо-функциональных, физиологических и биохимических свойств по мере развития различных заболеваний [3; 4]. Изучению оптико-спектральных характеристик биообъектов посвящено много работ, в которых приводятся результаты исследований тканей головного мозга, дыхательной, сердечно-сосудистой и пищеварительной систем. Системати-
зированные и обобщённые результаты этих исследований приведены в многочисленных статьях и обзорах [5; 6].
Целью данной работы является исследование воздействия злокачественных процессов на структурно-морфологические и биохимические свойства желчи и панкреатического сока. Для этого были выполнены спектрофотометрические измерения отражения и светопропускания в интервале длин волн 250-2500 нм и на основе метода Монте-Карло рассчитаны спектры коэффициентов оптического поглощения и анизотропии светорассеяния для пищеварительных ферментов при различных формах заболевания поджелудочной железы, включая процессы малигнизации.
Материалы и методы
Для проведения спектроскопических исследований в работе использовались образцы желчи и панкреатического сока, полученные в ГБУ Республики Дагестан «Республиканская клиническая больница» в ходе оперативного лечения и секционных исследований больных с за болеваниями поджелудочной железы различной степени тяжести. При этом патогистологические исследования биопсийного материала позволили отобрать по 7 образцов желчи и панкреатического сока для хронического панкреатита (биосреды при средней стадии поражения) и по 9 образцов желчи и панкреатического сока для случая аденокарциномы поджелудочной железы (образцы биосред при крайней стадии поражения).
Методика пробоподготовки заключалась в следующем. В целях минимизации спектральных вкладов поглощения воды образцы пищеварительных ферментов равномерно наносились на кварцевые стекла вдоль всей их поверхности и подвергались термическому (температура 40,0 + 2,0 °С) и вакуумному (давление 5х 10_6 мм рт. ст.) осушению. Процесс дегидратации последовательно повторялся до получения пленок биосред обоих видов толщиной примерно й = 100,0 + 10,0 мкм и площадью 2,5 X 4,0 см. При этом показатель преломления биосред принимался равным
1,4, а кварцевых пластин (Пд1а55) - 1,55.
Экспериментальные измерения спектров полного и диффузного светопропускания - и Та, а также диффузного отражения - для пищеварительных ферментов проводились в интервале длин волн 250-2500 нм на автоматизированном спектрофотометре иУ-3600 (8Ышаё2и, Япония) с интегрирующей сферой Ь18Я-3100 со спектральным разрешением ~0.1 нм.
Конечные данные спектрофотометрических коэффициентов , Та и определялись как р = Кд;(Л)_Ко(Л) и V = ^^ где адл) и ВД) - спектры пропускания и отражения проб биосред; и - спектры референсного сигнала, измеренные с кварцевой пластиной; Г0(А) - сигнал интегрирующей сферы с прикрытым входным и открытым выходным портами; Й0(А) - сигнал для сферы при открытых оптических портах.
Вычисление спектральной зависимости коэффициентов оптического поглощения (да), рассеяния и фактора анизотропии рассеяния (д) проводилось на основе инверсного метода численного моделирования Монте-Карло, разработанного и подробно описанного в работах [7-10]. Алгоритм используемого метода вкратце представляет собой последовательную реализацию следующих шагов:
1) вычисление начальных приближений для коэффициентов , и д на основе
г\6Х^ ГТ1&ХО гг1 СлР п п
данных кй , ^ и 1(1 с использованием соотношения [11].
gехр d <0.1
M'
ld
lnT?xpln0.05 „m
L ttwtx П^^Г
, рхп -,при RH < 0.1
(Ha+H's)d = { lnRedxp P d
21+s(RdXP~TtXP), при Redxp> 0.1
= -1ПТГ,
где = (1 — - коэффициент транспортного рассеяния, d - толщина образцов,
Tl бЛО rrt СлР rrt СлР rrt СлР
с - коллимированное пропускание, lc = lt ld ;
2) вычисление параметров Rdalc, T£alc и Tdalc на основе полученных значений коэффициентов , и g с использованием метода Монте-Карло;
3) построение целевой функции и ее минимизация: F = (—jcaïc^
-rdcaic)2 + (Redxp — Rdalc)2;
4) процедуру минимизации на основе симплексного метода Нелдера-Мида [12] до выполнения условия
|„ехр „calc\ \Texp_Tcalc\ , ехр „cale I
HL-ZIÉ_I + ¡1*__I + lRt ~Rt_I <0 01
„ехр ~ „ехр ~ Dexp 't 'd Kt
Для каждой исследуемой пробы пищеварительных ферментов вычисления выполнялись по несколько раз. Окончательный результат по спектрам оптического поглощения и анизотропии светорассеяния определялся путем усреднения серийных измерений: о = I—^^^ —, где п - количество серийных измерений, aj - значения коэффициентов , и g для i-й пробы, ô - среднее значение оптических коэффициентов, определяемое как - YH= i
Результаты и обсуждение
На рисунках 1 и 2 приводятся спектры коэффициентов диффузного отражения ('Rd), полного (Tt) и коллимированного (Гс) светопропускания пищеварительных ферментов печени и поджелудочной железы при средней и крайней стадиях поражения. Для формы спектрального контура параметров (Rd, Tt и Тс) обнаружены как общие закономерности, так и отличительные особенности. Так, например, в видимой области спектра трансмиссионные и отражательные коэффициентов характеризуются 3 минимумами, образованными, по-видимому, поглощением биллирубина и комплексом аминокислот, входящих в состав желчи и панкреатического сока [5, 6, 13, 14]. В ближней инфракрасной области спектра значения этих коэффициентов растут с образованием несущественных экстремумов в области полос поглощения воды [13, 14]. При этом в абсолютных пределах значения спектрофотометрических параметров для проб пищеварительных ферментов при средней стадии поражения составили -rt~0.15 ± 0.02 - 0.6 ± +0.1, Rd~0.05 ± 0.01 - 0.5 ± 0.1 и Гс~0.01 - 0.09. В то же время по мере развития процессов малигнизации в тканях поджелудочной железы значения Rd, Tt и Тс уменьшаются в среднем на 20 % для проб желчи и растут на 15 % для
панкреатического сока, что, как следствие, отражается на динамике спектральной зависимости оптических коэффициентов , и д.
0,6
0,5 -
0,4 -
0,3 -
0,2 -
0,1 -
0,0 ■ 0,6
0,5 0,4 Ь 0,3 • 0,2 0,1 0,0
Т.
0,09
0,07
-0,05 Т
0,03
1—1—I—1—I—1—I—1—I—1—I—1—I—1—I—1—Г
1 1 1 1 1 1 1 1 1 1 1 1 1 I 1 п
1 1 1 0 • ♦ 1 1 1
0,10
0,07 Т
0,04
0,01
250 500 750 1000 1250 1500 1750 2000 2250 2500
X, 11111
к
Т
Рис. 1. Спектры коэффициентов полного пропускания, полного отражения и коллимированного пропускания для проб желчи (1) и панкреатического сока (2)
при средней стадии поражения
0,6
0,4-
0,2 -
0,0 0,6
0,4-
0,2 -
0,0-
0,09
0,07
-0,04 Т
0,03
"1—1—I—1—I—1—I—1—I—1—I—1—I—1—I—1—I—1—Г
1 1 1 . | . | . | . | . | . г 1 •
• V........ - - к.|............. -
9 •
0,01 0,16
-0,13
-0,10 Т
-0,07
■0,04
250 500 750 1000 1250 1500 1750 2000 2250 2500
X, пш
Рис. 2. Спектры коэффициентов полного пропускания, полного отражения и коллимированного пропускания для проб желчи (1) и панкреатического сока (2)
при крайней стадии поражения
к
к
На рис. 3 и 4 показана спектральная зависимость коэффициентов поглощения (да), рассеяния (р.3) и фактора анизотропии рассеяния ($), рассчитанная на основе данных , и Тс при помощи инверсного метода МС для проб пищеварительных ферментов и соответствующих им форм заболеваний. Анализ полученных результатов показывает, что в спектрах оптических коэффициентов также имеется множество общих закономерностей. В частности, спектр поглощения для проб желчи характеризуется наличием явно выраженных экстремумов вблизи длин волн 280±5, 410±5 и 630±15 нм, соответствующих, скорее всего, полосам поглощения смеси желчных пигментов - билирубина и биливер-дина. Для проб панкреатического сока контур поглощения имеет максимум на длине волны 280±5 нм, образованный, по-видимому, поглощением комплекса аминокислот, входящих в состав панкреатического фермента [13-16]. С ростом длины волны коэффициент поглощения для обоих ферментов монотонно убывает, где в интервале длин волн 900-2500 нм достигает минимальных значений - 0.07±0.02 мм-1.
Следует заметить, что обезвоживание проб ферментов, проведенное для оценки спектральных вкладов эндогенных хромофоров, не выявило компонентов, соответствующих поглощению липидов (1220, 1730 и 1750 нм). При этом наблюдались спектральные минимумы на длинах волн 1450±20 и 1940±20 нм (полосы поглощения свободной воды) на фоне неспецифического поглощения продуктов окисления биохимических компонентов (холестерина, фосфолипидов, липидов и др.). Как показано в [14, 15, 17], «...продукты окисления могут быть образованы за счет вакуумно-термической дегидратации биосред, а также наличием в них элементов - 2п2+ и Си2+, способствующих процессу окисления». Это, по-видимому, легло в основу агрегации пищеварительных ферментов, из-за чего пленки исследуемых веществ представляют собой мутные сильнорассеивающие среды с выраженной коллоидной структурой.
Аналогичная спектральная зависимость наблюдалась для параметров анизотропии светорассеяния - и д, значения которых достигают максимума в видимой области спектра на участке длин волн 250-500 нм и с ростом длины волны монотонно спадают с образованием незначительных экстремумов вблизи длин волн 1450±20 и 1940±20 нм. В то же время в сравнении со спектром коэффициента поглощения для спектрального контура коэффициента рассеяния пищеварительных ферментов наблюдается обратная симметрия полос поглощения билирубина и аминокислот на длинах волн 280±5, 410±5 и 640±15 нм. Данный факт может быть обусловлен несколькими причинами, одной из которых является резкий рост мнимой составляющей комплексного показателя преломления среды (. п") в спектральной области сильного поглощения. Другим механизмом обнаруженной зависимости является эффект усиления оптического поглощения при больших значениях рассеяния света. Как показано в работах [14, 17], «. высокий уровень светорассеяния приводит к тому, что на длине свободного пробега фотоны взаимодействуют с большим числом эндогенных хромофоров, вследствие чего количество рассеянных фотонов резко падает». Данный эффект хорошо известен при оптических исследованиях других сильнорассеивающих объектов [9-10].
Применение методики среднеквадратического отклонения для выявления погрешности в определении оптических коэффициентов позволило установить, что максимальный статистический разброс значений наблюдался для коэффициента поглощения, который достигал о^а~40 % в видимой области спектра длин волн 280±5, 410±5 и 640±15 нм. В сравнении с этим статистический разброс значений для коэффициента рассеяния и фактора анизотропии был значительно ниже и не превышал порога в о^д < 10 % от величины полезного сигнала по всему спектральному диапазону проведенных исследований.
mm
mm
750 1 000 1 250 1 500 1 750 2000 2250 2500 X, nm
0,9
M
g
0,0
0,8
0,6
M
g
0,2
0,0
250
Рис. 3. Спектры оптических коэффициентов поглощения, рассеяния и фактора анизотропии рассеяния для проб желчи (1) и панкреатического сока (2) при средней стадии патологического
поражения
50
м„
Ms
mm'
5-=
0,05 50
m.
Ms
mm-
0,5
0,05 -L
0,9
5-, г.
-0,6
g
-0,3
0,6
0,4 g
0,2
250 500 750 1000 1250 1500 1750 2000 2250 2500
X, nm
0
8
0
Рис. 4. Спектры оптических коэффициентов поглощения, рассеяния и фактора анизотропии рассеяния для проб желчи (1) и панкреатического сока (2) при крайней стадии
патологического поражения
Вместе с этим оптическим спектрам пищеварительных ферментов и соответствующих им патологий свойственны явные отличительные особенности, что в принципе указывает на возможность использования методов спектрофотометрии в дифференциации
заболеваний поджелудочной железы. Исследуемые вещества обнаруживают существенные различия в значениях коэффициентов и которые при средней стадии поражения для проб желчи превышают значения таковых для проб панкреатического сока в 1,5 раза. При крайней стадии поражения эти различия увеличиваются в 2-3 раза. Вместе с этим по мере развития онкологической патологии для проб желчи наблюдается рост значений коэффициентов и до 1.5 раз в видимой области спектра. Для проб панкреатического сока при крайней стадии патологического поражения эта зависимость носит обратный характер - наблюдается уменьшение коэффициентов и до 2 раз. При этом спектры фактора анизотропии д для ферментов обоих видов менялись несущественно и с развитием патологии смещались в длинноволновую область спектра.
Особый интерес вызывает анализ параметров анизотропии рассеяния и д исследуемых ферментов с патологией. В частности, анализ коэффициента рассеяния посредством аппроксимации двухстепенной функцией типа = (а и № - безразмерные параметры: а является функцией концентрации рассеивателей и характеризует общий уровень рассеяния; № определяет средний размер рассеивающих частиц и формирует спектральный контур Д5(А)) позволил установить, что для исследуемых объектов с ростом длины волны возрастает роль рэлеевского рассеяния и снижается вклад рассеяния Ми. Результат такой аппроксимации показан на рис. 5, для проб желчи
Ь п / Ь сг\
в норме - и при крайней степени поражения (р.5~ ) принимает соответственно вид: М5Ь-П(Л) = 2.669 х 108А-3 293 + 68.517А"0124, М5Ь-СГ(Л) = 2.621 х 108А"3148 + 73.194А"0122.
'П.
Для панкреатического сока в норме - и при крайней степени поражения -
р сг
МГ :
дР-п(А) = 3.478 х 108Г 3-544 + 72.325Г0147, дР-сг(А) = 3.411 х 108Г 3-702 + 66.9 7 6А"0149.
Тот факт, что аппроксимацию коэффициента рассеяния провели при помощи функции двухстепенного вида, свидетельствует в пользу формирования спектров коэффициента как минимум двумя типами рассеивающих частиц. Причем первое слагаемое функции отвечает за светорассеяние, вызванное достаточно малыми частицами (близкими к рэлеевским рассеивателям), для которых волновой экспонент 4. Тогда как второй член функции соответствует крупным (рассеиватели Ми) частицам, для которых ж<1.
Ранее проведенные морфологические исследования методом РЭМ [17] показали, что пробы пищеварительных ферментов представляют собой коллоиды, для которых размер дисперсной фазы и их плотность зависят от вида биосред и формы заболевания поджелудочной железы. В частности, для исследуемых биосред в качестве малых рас-сеивателей выступают гранулы размером от 90-250 нм, а в качестве крупных частиц -их агломераты диаметром >2 мкм, образовавшиеся в результате окислительной агрегации ферментов. Вместе с этим анализ аппроксимирующих функций показывает, что развитие злокачественных процессов в тканях поджелудочной железы для проб желчи приводит к увеличению, а для проб панкреатического сока - к уменьшению размеров рассеивающих частиц при практически неизменной их концентрации. Обнаруженная морфологическая зависимость находит свое подтверждение как в динамике нисходящей с ростом длины волны спектральной зависимости фактора анизотропии рассеяния д, так и в результатах электронной микроскопии проб пищеварительных ферментов с патологией [17].
X, нм
Рис. 5. Аппроксимация двухстепенной функцией спектров рассеяния проб желчи (А) и панкреатического сока (В) в норме (кривая 1) и при крайней стадии патологического поражения (кривая 2). Пунктиром показана кривая аппроксимации
Выводы
Выполненные спектрально-оптические исследования для проб желчи и панкреатического сока при различных заболеваниях поджелудочной железы позволяют выделить следующие их структурно-морфологические и биохимические свойства.
1. Спектры коэффициента поглощения для проб желчи сформированы поглощением билирубина и биливердина на длинах волн 280±5, 410±5 и 630±15 нм, а для проб панкреатического сока - поглощением комплекса аминокислот на длине волны 280±5 нм. С ростом длины волны коэффициенты поглощения обоих ферментов убывают, где в интервале длин волн 900-2500 нм уменьшаются более чем в 100 раз.
2. Пробы пищеварительных ферментов представляют собой коллоиды среды, образованные как малыми рассеивателями, так и их крупными фракциями, сформированными в процессе окислительной агрегации биосред. Это подтверждается спектральной зависимостью фактора анизотропии д и коэффициента рассеяния для которых с ростом длины волны наблюдается увеличение вклада рэлеевского рассеяния и снижение рассеяния Ми.
3. Развитие аденокарциномы в тканях поджелудочной железы для проб желчи приводит к росту поглощения и рассеяния в 1.5 раза, в то время как для проб панкреа-
тического сока происходит уменьшение коэффициентов поглощения и рассеяния до 2-х раз.
Работа выполнена с использованием приборного парка ЦКП «Аналитическая спектроскопия» Дагестанского государственного университета.
Литература
1. Блохин Н.Н., Итин А.Б., Клименков А.А. Рак поджелудочной железы и внепече-ночных желчных путей. - М.: Медицина, 1982. - 271 с.
2. Маев И.В., Кучерявый Ю.А. Болезни поджелудочной железы: руководство. -М.: ГЭОТАР-МЕДИА, 2009. - 736 с.
3. Welch A.J., Van GemertM.J.C. Optical-Thermal Response of Laser Irradiation Tissues. Ed. by 2nd edition: © Springer Science+Business Media B.V. - 2011. - 958 p.
4. Тучин В.В. Оптика биологических тканей. Методы рассеяния света в медицинской диагностике. - М.: Физматлит, 2013. - 812 с.
5. Jacques S.L. Optical properties of biological tissues: a review // Phys. Med. Biol. 2013. - V. 58. - R. 37-61.
6. Bashkatov A.N., Genina E.A., Tuchin V.V. Tissue. Optical Properties // Handbook of Biomedical Optics. Ed. by D.A. Boas, C. Pitris and N. Ramanujam. Boca Raton, London, New York: CRC Teylor & Francisc Group, 2011. - P. 67-100.
7. Wang L.-H., Jacques S.L., Zheng L.-Q. MCML - Monte Carlo modeling of photon transport in multi-layered tissues // Computer Methods and Programs in Biomedicine. - 1995. - V. 47. - P. 131-146.
8. Wang L.-H., Jacques S.L., Zheng L.-Q. CONV - Convolution for responses to a finite diameter photon beam incident on multi-layered tissues // Computer Methods and Programs in Biomedicine. - 1997. - V. 54. - P. 141-150.
9. Башкатов А.Н., Генина Э.А., Козинцева М.Д. и др. Оптические свойства биологических тканей брюшины в спектральном диапазоне 350-2500 нм // Оптика и спектроскопия. - 2016. - Т. 120. - С. 6-14.
10. Башкатов А.Н., Генина Е.А., Кочубей В.И. и др. Оптические свойства тканей толстой кишки человека в спектральном диапазоне 350-2500 нм // Квантовая электроника. - 2014. - Т. 44. - С. 779-784.
11. Prahl S.A., Van Gemert M.J.C., Welch A.J. Determining of the optical properties of turbid media by using the adding-doubling method // Appl. Optics. - 1993. - V. 32. - P. 359568.
12. Банди Б. Методы оптимизации. - М.: Радио и связь, 1988. - 225 с.
13. Орегонский медицинский - биологический центр. Доступно: http://omlc.org/spectra/PhotochemCAD/html/136.html.
14. Гираев К.М., Ашурбеков Н.А., Магомедов М.А. и др. Влияние процессов ма-лигнизации в тканях поджелудочной железы на формирование оптических и инфракрасных спектров проб желчи и панкреатического сока // Биофизика. - 2017. - Т. 62. -С. 784-793.
15. Ленинджер А. Биохимия. - М.: Мир, 1976. - 958 с.
16. Nachabe R., Evers D.J., Hendriks B.H.W. et al. Effect of bile absorption coefficients on the estimation of liver tissue optical properties and related implications in discriminating healthy and tumorous samples // Biomedical Optics Express. - 2011. - V. 2. - P. 600614.
17. Гираев К.М., Ашурбеков Н.А., Магомедов М.А. и др. Влияние патологических процессов на оптические спектры поглощения и рассеяния проб желчи и панкреатического сока // Оптика и спектроскопия. - 2015. - Т. 119. - С. 162-171.
Поступила в редакцию 3 декабря 2018 г.
UDC 535.33/.34; 535.3
DOI: 10.21779/2542-0321-2018-33-4-41-50
The effect of malignant neoplasms in pancreatic tissues on the optical spectra
of digestive secretions
A.A. Murtazaeva1,2, K.M. Giraev1,2, М.А. Magomedov1,2, N.A. Ashurbekov1
1 Dagestan State University; Russia, 367001, Makhachkala, M. Gadzhiev st., 43a;
2 Amirkhanov Institute of Physics Dagestan Scientific Center of the Russian Academy of Sciences; Russia, 367003, Makhachkala, М.Yaragsky st., 94; nashurb@mail.ru
The spectra of transmittance and reflection coefficients are measured for samples of bile and pancreatic juice in various forms of pancreatic pathology in the wavelength range of 250-2500 nm. The spectra of optical coefficient of absorption, scattering and scattering anisotropy factors are determined using the inverse Monte Carlo method. It is found, that the researched objects are colloids consisting mainly of small particles and their large agglomerates, formed during dehydration and oxida-tive aggregation of biological media. It is shown, that the development of malignance for bile samples leads to an increase, and for pancreatic juice - to a decrease of the density and size of particles, which follows from the increase and decrease in the absorption and light scattering coefficients.
Keywords: bile, pancreatic juice, malignancy processes, spectrophotometry, absorption, scattering, anisotropy factor, morphology, diagnostics.
Received 3 December, 2018