Влияние кислорода на распределение повреждений ДНК в 3D геноме и его корреляция с коэффициентом усиления кислородом при облучении высокой ЛПЭ

Основная информация

• ID статьи: 2503.21837
• Название: Impact of Oxygen on DNA Damage Distribution in 3D Genome and Its Correlation to Oxygen Enhancement Ratio under High LET Irradiation
• Авторы: Ankang Hu, Wanyi Zhou, Xiyu Luo, Rui Qiu, Junli Li
• Учреждение: Департамент инженерной физики Университета Цинхуа, Институт физики жидкостей Китайской академии инженерной физики
• Классификация: physics.bio-ph physics.med-ph
• Ссылка на статью: https://arxiv.org/abs/2503.21837

Аннотация

Изменение коэффициента усиления кислородом (OER) при различных линейных передачах энергии (ЛПЭ) в настоящее время не имеет полного механистического объяснения и механистической модели. Предыдущие исследования авторов выявили значительную корреляцию между распределением двухцепочечных разрывов (ДЦР) в 3D геноме и радиационно-индуцированной гибелью клеток, что обеспечивает ценные сведения для понимания кислородного эффекта. В данном исследовании предложена модель, которая представляет реакцию кислорода как вероятность индукции разрывов цепей ДНК и интегрирует её в монте-карловское моделирование трековой структуры для изучения влияния кислорода на распределение ДЦР в 3D геноме. Результаты показывают, что соотношение частот возникновения скопленных ДЦР в одном топологически ассоциированном домене (ТАД) (случай 2) и ДЦР в часто взаимодействующих ТАДах (случай 3) при аэробных и анаэробных условиях согласуется с тенденцией OER выживаемости клеток при различных ЛПЭ.

Исследовательский контекст и мотивация

Основные проблемы

1. Отсутствие механизма OER: Изменение коэффициента усиления кислородом (OER) в зависимости от линейной передачи энергии (ЛПЭ) не имеет полного механистического объяснения, особенно механизм снижения OER при облучении высокой ЛПЭ неясен
2. Ограничения традиционных моделей: Существующие модели основаны главным образом на эмпирической подгонке параметров и не имеют механистического объяснения с точки зрения распределения повреждений ДНК
3. Важность структуры 3D генома: Традиционные исследования сосредоточены на общем количестве ДЦР, игнорируя важность пространственного распределения ДЦР в 3D геноме

Значимость исследования

• Оптимизация лучевой терапии: Понимание механизма OER имеет решающее значение для оптимизации стратегий лучевой терапии, особенно при лечении гипоксических опухолей
• Разработка радиосенсибилизаторов: Обеспечивает теоретическую основу для разработки радиосенсибилизаторов на основе принципов кислородного эффекта
• Механизм FLASH-эффекта: Обеспечивает поддержку для понимания механизма кислородного эффекта в современной FLASH-лучевой терапии

Ограничения существующих методов

1. Эмпирические модели: Большинство моделей OER зависят от эмпирических параметров прямой подгонки данных выживаемости клеток
2. Упрощение химических реакций: Хотя некоторые модели пытаются включить коррекции реакций ДНК-радикалов с кислородом, эти корректировки остаются эмпирическими
3. Несоответствие масштабов: Кривые выхода ДЦР не совпадают с кривыми относительной биологической эффективности (ОБЭ), что указывает на необходимость анализа структуры генома в большем масштабе

Основной вклад

1. Инновационная механистическая модель: Предложена модель механистического объяснения OER на основе распределения ДЦР в 3D геноме, моделирующая реакцию кислорода как вероятность индукции разрывов цепей ДНК
2. Интеграция монте-карловского моделирования трековой структуры: Интегрирует кислородный эффект в монте-карловское моделирование трековой структуры для изучения влияния кислорода на распределение ДЦР в 3D геноме
3. Классификация повреждений на уровне ТАД: Устанавливает систему классификации трёх типов моделей распределения ДЦР на основе топологически ассоциированных доменов (ТАД)
4. Экспериментальная верификация: Расчётные результаты OER демонстрируют хорошее соответствие экспериментальным данным, подтверждая эффективность модели
5. Объяснение механизма высокой ЛПЭ: Обеспечивает потенциальное механистическое объяснение снижения OER при облучении высокой ЛПЭ

Подробное описание методов

Определение задачи

Входные данные: Радиационные частицы различной ЛПЭ (электроны, протоны, ионы He-3, ионы C-12) и их энергия Выходные данные: Модели распределения ДЦР в 3D геноме, частота возникновения различных типов повреждений, значения OER Ограничения: На основе модели ядра клетки IMR90, учитывающей структуру ТАД и взаимодействия

Архитектура модели

1. Модель выживаемости клеток 3D генома

На основе структуры ТАД распределение ДЦР разделено на три случая:

• Случай 1: Изолированные ДЦР - только один ДЦР в одном ТАД, без ДЦР в часто взаимодействующих с ним ТАДах
• Случай 2: Скопленные ДЦР в одном ТАД - два или более ДЦР в одном ТАД
• Случай 3: ДЦР в часто взаимодействующих ТАДах - один или более ДЦР в каждом из нескольких часто взаимодействующих ТАДов

Формула расчёта доли выживших клеток:

-ln S = p₁n₁(D) + p₂n₂(D) + p₃n₃(D)

где p₃ > p₂ > p₁, указывая на то, что случаи 3 и 2 имеют более высокую вероятность гибели клеток.

2. Механистическое моделирование кислородного эффекта

Влияние кислорода на радиационное повреждение осуществляется главным образом через реакцию с ДНК-радикалами, включая три конкурирующих процесса:

1. Фиксация кислородом: Реакция с кислородом с образованием пероксирадикалов
2. Естественная фиксация: Преобразование в повреждения, требующие биологического восстановления
3. Химическое восстановление: Восстановление антиоксидантами, такими как глутатион

Вероятность преобразования ДНК-радикала в повреждение:

p_lesion([O₂]) = (k₁·[O₂] + k₂)/(k₁·[O₂] + k₂ + k₃·[antioxidants])

3. Монте-карловское моделирование трековой структуры

• Физическое моделирование: Использование Geant4-DNA для моделирования событий отложения энергии
• Модель ядра клетки: Применение модели ядра клетки IMR90 на основе данных Hi-C, содержащей 15 282 ТАД
• Анализ повреждений: Использование алгоритма DBSCAN для анализа распределения и скопления ДЦР

Технические инновации

1. Многомасштабная интеграция: Объединение молекулярного масштаба окислительных химических реакций со структурой 3D генома на уровне ядра клетки
2. Вероятностная обработка: Упрощение сложных окислительных химических реакций до вероятности индукции разрывов цепей при отложении энергии p_DNA
3. Анализ на основе ТАД: Классификация повреждений на основе биологически значимой структуры ТАД, а не простого пространственного расстояния
4. Статистические методы: Применение предположения распределения Пуассона для обработки статистической вариативности результатов моделирования

Экспериментальная установка

Параметры моделирования

• Модель ядра клетки: Клетки IMR90, радиус 5,0 мкм, объёмная доля 0,15
• Параметры ТАД: 15 282 ТАД, средний радиус 89,12 нм, среднее содержание ДНК 397 кбп
• Типы облучения: Электроны 0,9 МэВ (низкая ЛПЭ), протоны 0,4-66,46 МэВ, ионы He-3 0,8-23,0 МэВ/u, ионы C-12 2,98-209,46 МэВ/u
• Параметры кислородного эффекта: Диапазон p_DNA 0,04-0,22, аэробные условия p_DNA=0,14, анаэробные условия p_DNA=0,09

Показатели оценки

• OER_α: Коэффициент усиления кислородом на основе параметра α
• OER₁₀: Коэффициент усиления кислородом при 10% выживаемости
• Соотношение выхода повреждений: Соотношение однозепочечных разрывов (ОЗР), ДЦР и различных случаев в аэробных/анаэробных условиях

Процедура моделирования

1. Физическое моделирование с использованием Geant4 для генерации событий отложения энергии
2. Определение потенциальных повреждений на основе энергетического порога и вероятности p_DNA
3. Случайное распределение на две цепи ДНК двойной спирали
4. Определение ДЦР: два или более разрыва цепей на противоположных цепях в пределах 3,2 нм
5. Анализ распределения ДЦР с использованием алгоритма DBSCAN
6. Статистика частоты возникновения трёх случаев и расчёт доли выживших клеток

Результаты экспериментов

Основные результаты

1. Влияние p_DNA на распределение повреждений

• Низкая ЛПЭ: Выход ОЗР линейно зависит от p_DNA, ДЦР и сложные повреждения показывают нелинейную зависимость
• Высокая ЛПЭ: Все типы повреждений демонстрируют более сложные нелинейные отношения
• Классическая формула Альпера: Эффективна только для простых разрывов цепей, отклоняется от предсказаний для сценариев сложных повреждений

2. Влияние кислорода на распределение повреждений при различных ЛПЭ

С увеличением ЛПЭ:

• Соотношение выхода ОЗР и ДЦР снижается
• Соотношение частоты возникновения случаев 2 и 3 значительно снижается
• Кривая соотношения случаев 2 и 3 высоко согласуется с тенденцией OER выживаемости клеток

3. Результаты расчёта OER

• OER_α: Остаётся относительно постоянным в диапазоне ЛПЭ десятков кэВ/мкм
• OER₁₀: Показывает небольшое увеличение в том же диапазоне ЛПЭ, в конечном итоге стремясь к OER_α
• Сравнение с экспериментом: Расчётная кривая OER демонстрирует хорошее соответствие экспериментальным данным

Ключевые находки

1. Механистическое объяснение: Облучение высокой ЛПЭ создаёт плотную ионизацию, приводящую к избыточным повреждениям, и даже если часть восстанавливается химически, оставшихся повреждений достаточно для индукции случаев 2 и 3
2. Важность масштаба: Распределение ДЦР в 3D геноме лучше предсказывает гибель клеток, чем общее количество ДЦР
3. Корреляция с ТАД: Классификация повреждений на основе ТАД высоко коррелирует с изменениями OER, подтверждая важность структуры 3D генома

Связанные работы

Традиционные исследования кислородного эффекта

1. Гипотеза фиксации кислородом: Кислород реагирует с ДНК-радикалами, образуя неремонтируемые пероксирадикалы
2. Монте-карловское моделирование трековой структуры: Включение кислородного эффекта в монте-карловское моделирование, но главным образом сосредоточено на общем количестве ДЦР
3. Специальные механизмы высокой ЛПЭ: Гипотеза о множественной ионизации, производящей кислород, но с ограниченным влиянием

Исследования повреждений ДНК в 3D геноме

• Моделирование ядра клетки на основе данных Hi-C
• Влияние структуры ТАД на восстановление ДНК
• Связь трёхмерной структуры хроматина с радиационной чувствительностью

Методы моделирования OER

• Эмпирическая подгонка формул
• Модели химической кинетики реакций
• Микродозиметрические методы

Выводы и обсуждение

Основные выводы

1. Новое механистическое объяснение: Обеспечивает потенциальное механистическое объяснение изменения OER в зависимости от ЛПЭ на основе распределения ДЦР в 3D геноме
2. Важность масштаба ТАД: Подтверждает важность распределения повреждений ДНК на уровне ТАД для предсказания гибели клеток
3. Характеристики высокой ЛПЭ: Облучение высокой ЛПЭ снижает относительную важность кислородного эффекта путём создания избыточных повреждений

Ограничения

1. Упрощение модели: Упрощение сложных окислительных химических реакций до вероятностной модели может упустить некоторые детали
2. Единственный тип клеток: Использование только модели клеток IMR90 требует проверки применимости к другим типам клеток
3. Разброс экспериментальных данных: Экспериментальные данные OER при высокой ЛПЭ сильно разбросаны, ограничивая количественное сравнение
4. Детали химических реакций: Игнорирование вторичных факторов, таких как производство кислорода при множественной ионизации при высокой ЛПЭ

Направления будущих исследований

1. Проверка на нескольких типах клеток: Расширение на различные клеточные линии для проверки универсальности модели
2. Уточнение химических реакций: Включение более детальных механизмов химических реакций
3. Применение к FLASH-эффекту: Применение к пониманию кислородного эффекта в FLASH-лучевой терапии
4. Клиническая трансляция: Исследование потенциального применения в клинических планах лучевой терапии

Глубокая оценка

Преимущества

1. Высокая инновационность: Впервые объединяет структуру 3D генома с механизмом OER, обеспечивая новую исследовательскую перспективу
2. Методологическая строгость: Многодисциплинарный подход, объединяющий монте-карловское моделирование трековой структуры, статистический анализ и биологическое моделирование
3. Экспериментальная верификация: Хорошее соответствие расчётных результатов экспериментальным данным повышает достоверность модели
4. Механистическое объяснение: Обеспечивает разумное объяснение давно отсутствовавшего механистического объяснения изменения OER при высокой ЛПЭ

Недостатки

1. Сложность модели: Включение нескольких масштабов и дисциплин усложняет верификацию и применение модели
2. Зависимость от параметров: Производительность модели зависит от нескольких эмпирических параметров, требующих калибровки на основе большего количества экспериментальных данных
3. Вычислительные затраты: Монте-карловское моделирование и анализ 3D генома требуют значительных вычислительных ресурсов
4. Биологические предположения: Классификация повреждений на основе структуры ТАД требует дополнительной биологической экспериментальной верификации

Влияние

1. Научный вклад: Обеспечивает новую теоретическую базу и направление исследований для области радиационной биологии
2. Потенциал применения: Может использоваться для оптимизации планов лучевой терапии, особенно при тяжёлоионной терапии и лечении гипоксических опухолей
3. Междисциплинарная ценность: Способствует интеграции исследований физики, биологии и медицины
4. Методологическая ценность: Обеспечивает пример моделирования, интегрирующего многомасштабные биофизические процессы

Применимые сценарии

1. Планирование лучевой терапии: Оптимизация планов лучевой терапии с учётом состояния оксигенации опухоли
2. Радиационная защита: Оценка биологического эффекта различных типов облучения
3. Разработка лекарств: Механистические исследования и прогнозирование эффективности радиосенсибилизаторов
4. Фундаментальные исследования: Углубленное понимание механизмов восстановления повреждений ДНК

Библиография

Статья цитирует 38 важных работ, охватывающих:

• Классическую теорию кислородного эффекта (формула Альпера, гипотеза фиксации кислородом)
• Методы монте-карловского моделирования трековой структуры
• Технологии моделирования 3D генома
• Экспериментальные данные OER
• Исследования механизмов восстановления повреждений ДНК

Общая оценка: Это высокоинновационная междисциплинарная исследовательская статья, успешно интегрирующая структуру 3D генома в исследование механизма OER и обеспечивающая разумное механистическое объяснение давно озадачивавшего область изменения OER при облучении высокой ЛПЭ. Несмотря на некоторые упрощения модели и потребность в верификации, её инновационный подход и важные находки открывают новые направления исследований в области радиационной биологии и лучевой терапии.