МЕТОД
Количественное определение йода при помощи рентгеновской компьютерной томографии для дозиметрического обеспечения фотон-захватной терапии
1 Российский онкологический научный центр имени Н. Н. Блохина, Москва
2 Федеральный медицинский биофизический центр имени А. И. Бурназяна, Москва
3 Национальный исследовательский ядерный университет «МИФИ», Москва
4 Австрийский технологический институт, Вена, Австрия
Для корреспонденции: Алексей Андреевич Липенгольц
Каширское шоссе, д. 24, г. Москва, 115478; ur.liam@stlognepil
Эффективность и перспективность фотон-захватной терапии (ФЗТ) в лечении злокачественных опухолей, особенно опухолей мозга, показана во многих российских и зарубежных исследованиях [1, 2, 3, 4, 5]. Данный метод основан на преимущественном поглощении внешнего рентгеновского излучения ортовольтного диапазона (от 30 до 300 кэВ) дозоповышающими агентами (ДПА), локализованными в опухоли: это приводит к увеличению поглощаемой опухолью дозы излучения. ДПА в ФЗТ являются химические элементы с порядковым номером Z >52, которые обладают большей поглощающей способностью по сравнению с химическими элементами, составляющими биологические мягкие ткани (H, C, O, N и др.) [6].
Первые клинические испытания фазы I ФЗТ опухолей мозга были проведены в США на модифицированном компьютерном томографе в 1990-х гг. [2]. В настоящее время такие испытания проводят во Франции в Европейском центре синхротронного излучения [7]. Несмотря на обнадеживающие результаты, полученные в этих исследованиях, ФЗТ может так и остаться на этапе клинических испытаний фазы I и не перейти к следующим фазам, если не будет решена проблема получения количественных данных о распределении ДПА в опухоли и окружающих ее тканях до и во время сеанса облучения. Дело в том, что величина поглощенной дозы при ФЗТ определяется величиной концентрации ДПА в мишени и в зависимости от нее может увеличиваться в 1,5–5 раз по сравнению с дозой, поглощенной такой же, но не содержащей ДПА мишенью [8, 9, 10, 11, 12]. В связи с этим планирование облучения и контроль величины доставленной поглощенной дозы в ткани пациента, необходимые для продолжения клинических испытаний, невозможны без разработки метода получения количественных данных о содержании ДПА в опухолевых и здоровых тканях.
Планирование облучения в обоих упомянутых исследованиях осуществлялось без учета нахождения в мишени ДПА; режим и длительность облучения определяли исходя лишь из взаимодействия рентгеновского излучения с мягкими тканями, без учета дополнительного энерговыделения на атомах ДПА. Поэтому режим облучения в данных исследованиях был фракционированным, аналогичным для традиционной дистанционной лучевой терапии. Оценка увеличения поглощенной дозы за счет взаимодействия с ДПА производилась лишь позднее: на этапе анализа полученных результатов, с использованием усредненных по всем пациентам данных рентгеновской компьютерной томографии (КТ), полученных на этапе подготовки исследования. При расчете величины поглощенной дозы распределение ДПА в мишени вынужденно считалось однородным. Определение параметров оценки содержания ДПА, а также их влияния на итоговую величину поглощенной опухолью дозы в данных исследованиях не проводилось. Очевидно, что достоверное количественное определение ДПА в тканях пациента in vivo является основополагающим фактором эффективного и безопасного применения ФЗТ в лечении злокачественных опухолей.
КТ является наиболее вероятным методом количественного определения ДПА при ФЗТ. КТ и ФЗТ основаны на одних и тех же физических принципах, а именно: на поглощении рентгеновского излучения в веществе. КТ широко применяется в клинической практике и является одним из основных методов медицинской визуализации. Принципиальная возможность количественного определения ДПА с помощью КТ очевидна, исходя из физических основ КТ, и была экспериментально подтверждена рядом исследователей [13, 14, 15]. Однако метрологические характеристики метода количественного определения ДПА при помощи КТ до сих пор остаются неисследованными.
Целью данной работы являлись демонстрация принципиальной возможности определения йода в тканеэквивалентном фантоме при помощи КТ, оценка различных характеристик метода, а также оценка влияния на них значения напряжения рентгеновской трубки КТ-сканера и изменения рентгеновского полихроматического спектра в толще тканеэквивалентного фантома, вызванного неодинаковым ослаблением его компонент.
МАТЕРИАЛЫ И МЕТОДЫ
В исследовании был использован специально изготовленный прямоугольный полиэтиленовый фантом размером 134 × 134 × 63 мм. В нем были проделаны два перпендикулярных ряда отверстий от поверхности к центру (вглубь) для размещения микропробирок объемом 250 мкл с водными растворами йода, имевшими концентрацию йода от 0,5 до 50 мг/мл (рис. 1). В качестве йодсодержащего вещества использовался йопромид в форме лекарственного средства «Ультравист-370» (Bayer Schering Pharma AG, Германия). Приготовление йодсодержащих растворов осуществлялось путем разбавления исходного раствора «Ультрависта-370», содержавшего 370 мг/мл йода, при помощи автоматических дозаторов. Визуализация фантома с микропробирками осуществлялась на томографе Siemens Biograph 40 (Siemens, Германия). Сканирование проводилось при напряжениях рентгеновской трубки КТ-сканера 80, 100, 120 и 140 кВ и токе трубки 200 мА. Реконструкция изображений осуществлялась c помощью программного ядра B30f стандартного программного обеспечения томографа. Количественный анализ томограмм в стандартном формате DICOM проводился при помощи программы ImageJ (National Institutes of Health, США). Измерение средних значений рентгеновской плотности (в единицах Хаунсфилда, HU) йодсодержащих растворов, а также величин стандартного отклонения осуществлялось встроенными средствами программы ImageJ. Линейная аппроксимация полученных средних значений осуществлялась в программной среде R (R Foundation).
РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЯ
Зависимость значений рентгеноплотности водных растворов йода от его концентрации в них приведена на рис. 2. Зависимость носит линейный характер (R2 = 0,998) во всем диапазоне исследованных концентраций — от 0,5 до 50 мг/мл. По полученным экспериментальным точкам были построены градуировочные кривые. Отклонение рассчитанных по градуировочным кривым значений концентрации йода для исследуемых растворов йода приведены в таблице. В диапазоне концентраций от 10 до 50 мг/мл отклонение не превышало 5 %. Для концентрации 5 мг/мл отклонение составило 5–10 %. Наибольшее отклонение (до 80 %) наблюдалось для растворов с концентрацией йода менее 1 мг/мл. Наблюдалось незначительное изменение величины рентгеноплотности в HU в зависимости от местоположения микропробирки с раствором в фантоме. Градуировочные кривые, построенные для растворов, распологавшихся около поверхности, систематически оказались ниже градуировочных кривых, построенных для растворов, располагавшихся в центре (в глубине) фантома. Однако различие между расчетными концентрациями одного и того же раствора по этим градуировочным кривым не превысило 2 % для всех использованных напряжений КТ-сканера.
ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ
Впервые проведено исследование количественного определения ДПА на клиническом компьютерном томографе с оценкой характеристик предлагаемого метода, необходимого для фотон-захватной терапии. Была исследована зависимость значений рентгеноплотности водных растворов йода от напряжения на рентгеновской трубке КТ-сканера во время измерения и от местоположения пробирок с растворами (расстояния от поверхности) в тканеэквивалентном фантоме. Зависимость значений рентгеноплотности водных растворов йода от его концентрации носит линейный характер в диапазоне концентраций 0,5–50 мг/мл для всех использованных значений напряжения трубки КТ-сканера. Как и ожидалось, для напряжения 80 кВ наблюдается большая чувствительность количественного определения йода, чем для напряжений 100, 120 и 140 кВ. Однако наблюдаемая разница в чувствительности непринципиальна для ФЗТ и не влияет на точность количественного определения ДПА с помощью КТ (таблица), на которую в большей степени влияет разброс HU в области проведения анализа. По этой причине выбор оптимального напряжения для количественного определения ДПА может быть сделан исходя из других, более важных требований ФЗТ [6, 16].
Несмотря на существенное отклонение измеренных значений концентраций йода от истинных значений в диапазоне концентраций от 0,5 до 5 мг/мл (до 80 %), следует отметить, что при планировании и оценке поглощенных доз в ФЗТ для этих концентраций такое отклонение приемлемо. Как нами было показано ранее [16], увеличение поглощенной дозы при ФЗТ зависит от концентрации ДПА, и данная зависимость может быть описана линейным законом (R2 = 0,99764). По этой причине для вычисления поглощенной дозы при ФЗТ имеет значение абсолютная, а не относительная величина погрешности количественного определения ДПА. Наблюдаемая неопределенность в 80 % для концентраций 1 мг/мл соответствует абсолютной погрешности в 0,8 мг/мл, что, в свою очередь, приводит к относительному изменению поглощенной дозы при ФЗТ не более чем на 2,8 %. Такая ошибка при определении величины поглощенной дозы является приемлемой в лучевой терапии.
Таким образом, предлагаемый метод количественного определения йода при помощи КТ позволяет использовать получаемые количественные данные о распределении ДПА для планирования облучения и контроля поглощенной дозы в ФЗТ. Точность количественного определения ДПА при помощи КТ, вероятно, может быть улучшена путем разработки специализированных алгоритмов реконструкции изображений, нацеленных на получение изображений с меньшим контрастом и резкостью, но при этом и меньшим разбросом HU для однородного рентгеноплотного объекта. Степень ужесточения спектра рентгеновского излучения в исследованном объекте существенно не влияет на количественное определение йода. При количественном относительном определении ДПА при помощи КТ in vivo возможно расположение образцов сравнения в непосредственной близости с пациентом и не требует создания и размещения дополнительного антропоморфного фантома для градуировки.
ВЫВОДЫ
Экспериментально подтверждена принципиальная возможность количественного определения йода при помощи КТ в тканеэквивалентном фантоме. Метод показал себя достаточно надежным, пригодным для планирования и дозиметрического сопровождения ФЗТ. Метод обладает стабильными характеристиками в широком для ФЗТ диапазоне напряжений трубки КТ-сканера и концентраций ДПА, а также для исследования различных по протяженности и форме объектов. КТ является наиболее универсальным для ФЗТ методом количественного определения широкого спектра ДПА.