Гормональная регуляция динамики заболеваемости раком молочной железы: математический анализ, объясняющий крючок Клеммесена

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

• ID статьи: 2511.19964
• Название: Hormonal Regulation of Breast Cancer Incidence Dynamics: A Mathematical Analysis Explaining the Clemmesen's Hook
• Авторы: Navid Mohammad Mirzaei, Wan Yang (Колумбийский университет)
• Классификация: q-bio.PE (Количественная биология - популяции и эволюция)
• Дата публикации: 26 ноября 2025 г.
• Ссылка на статью: https://arxiv.org/abs/2511.19964v1

Аннотация

Данное исследование посвящено математическому моделированию феномена "крючка Клеммесена" в заболеваемости раком молочной железы - замедления и последующего восстановления темпов заболеваемости в возрастном диапазоне 45-55 лет, соответствующем менопаузе. На основе расширенной многостадийной модели клональной экспансии опухоли (MSCE-T) авторы посредством строгого теоретического анализа и численных экспериментов доказывают, что данный феномен может быть количественно объяснен временной специфичностью изменений баланса пролиферации и апоптоза в популяции ранних мутантных клеток. Эти изменения соответствуют снижению уровня прогестерона вокруг менопаузы и переходу к альтернативным путям роста после менопаузы (таким как сигналы роста, происходящие из жировой ткани). Исследование также демонстрирует, что изменения поздней клональной динамики и эффекты обнаружения неспособны эффективно воспроизвести этот немонотонный паттерн заболеваемости.

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

1. Основная проблема

Рак молочной железы является наиболее распространенным раком и основной причиной смертности у женщин во всем мире. Кривая возрастной заболеваемости демонстрирует уникальный феномен "крючка Клеммесена" в период менопаузального перехода (приблизительно 45-55 лет): темп прироста заболеваемости сначала замедляется, а затем восстанавливается. Этот феномен наблюдается в различных популяциях и когортах, что предполагает наличие дифференциальных факторов риска вокруг менопаузы.

2. Значимость проблемы

Понимание биологического механизма крючка Клеммесена имеет решающее значение, поскольку:

• Он раскрывает, как гормональные изменения при менопаузе влияют на риск рака молочной железы
• Может служить основой для стратегий профилактики и вмешательства
• Способствует пониманию причинной роли гормонов в развитии рака

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

• Биологические доказательства: Недавние исследования указывают на прогестерон (а не эстроген) как основной канцерогенный фактор, но отсутствует количественное объяснение
• Традиционные модели MSCE: Неспособны захватить немонотонные тренды, требуют введения временно-зависимых параметров или динамики обнаружения
• Неясные механизмы: Механистическая основа крючка Клеммесена остается неясной, отсутствует математическая база, связывающая гормональный переход с кривой заболеваемости

4. Исследовательская мотивация

На основе новых эпидемиологических доказательств (количество овуляторных циклов, а не кумулятивное воздействие эстрогена является основным определяющим фактором) авторы выдвигают гипотезу, что изменения уровня прогестерона являются основным движущим фактором крючка Клеммесена, и проверяют это посредством математического моделирования.

Основные вклады

1. Теоретическая аналитическая база: Строгий математический анализ модели MSCE-T с выведением достаточных условий для возникновения крючка Клеммесена
2. Идентификация механизма: Доказательство того, что временная специфичность изменений соотношения пролиферации/апоптоза (α₁/β₁) в ранних мутантных клетках (первая стадия) может объяснить феномен, тогда как поздняя динамика и эффекты обнаружения не могут
3. Теорема нечувствительности к параметрам: Посредством представления интегралов Вольтерра и анализа производных Гато доказывается, что изменения порога обнаружения опухоли (Mt) и скорости пролиферации злокачественных клеток (α₃) оказывают слабое влияние на заболеваемость в среднем возрасте
4. Биологическая верификация: Связь математических результатов с прогестерон-управляемыми механизмами пролиферации, обеспечивающая количественные доказательства в поддержку гипотезы прогестерона
5. Численная верификация: Использование данных трех когорт рождения из базы данных SEER для верификации модели, подтверждающее, что умеренные изменения α₁ и β₁ достаточны для воспроизведения крючка Клеммесена

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

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

Определить, при каких условиях параметров модель MSCE-T может производить немонотонную кривую возраст-заболеваемость, подобную крючку Клеммесена:

• Входные данные: Параметры модели (скорости пролиферации αᵢ, скорости апоптоза βᵢ, скорости мутаций μᵢ) и их временная зависимость
• Выходные данные: Функция риска рака h(t) и ее производная h'(t)
• Ограничения: Параметры должны находиться в биологически обоснованном диапазоне

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

1. Основы модели MSCE-T

Рассматривается процесс развития рака молочной железы с тремя лимитирующими скорость мутациями в генах-драйверах, описываемый следующей системой ОДУ:

ẋ₁ = μ₀N₀x₁(x₃ - 1)
ẋ₂ = -μ₀N₀x₄
ẋ₃ = β₁ - (α₁ + β₁ + μ₁)x₃ + μ₁x₃x₅ + α₁x₃²
ẋ₄ = -(α₁ + β₁ + μ₁)x₄ + μ₁x₄x₅ + μ₁x₃x₆ + 2α₁x₃x₄
ẋ₅ = β₂ - (α₂ + β₂ + μ₂)x₅ + μ₂f(t)x₅ + α₂x₅²
ẋ₆ = -(α₂ + β₂ + μ₂)x₆ + μ₂f'(t)x₅ + μ₂f(t)x₆ + 2α₂x₆x₅

где:

• x₁, x₃, x₅ - вероятности выживания на каждой стадии
• h(t) = x₂(t) - функция риска рака
• x₄ = ẋ₃, x₆ = ẋ₅
• f(t) - интегрирует эффекты обнаружения

2. Функция обнаружения

f(t) = 1 - (1 - e^(-α₃t))^(Mt-1)

где α₃ - чистая скорость пролиферации опухолевых клеток, Mt - количество злокачественных клеток, необходимое для обнаружения в момент времени t.

3. Ключевые параметры

Определяются ключевые составные параметры:

c(t) = α₁(2x₃(t) - 1) - β₁ + μ₁(x₅(t) - 1)
q(t) = -(α₂ + β₂ + μ₂) + μ₂f(t) + 2α₂x₅(t)

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

1. Представление интегралом Вольтерра (раздел 4)

Посредством метода интегрирующего множителя решаются x₆ и x₄, получая:

x₄(t) = ∫₀ᵗ K(t,r)f'(r)dr

где ядро функции:

K(t,r) = μ₁μ₂x₅(r)∫ᵣᵗ x₃(s)exp(∫ᵣˢ(c+q)du)ds

Это представляет h'(t) = -μ₀N₀x₄(t) как преобразование Вольтерра f'(r), обеспечивая основу для последующего анализа чувствительности.

2. Теорема нечувствительности эффектов обнаружения (раздел 5)

Теорема 5.3 (нечувствительность масштабирования Mt): Для окна среднего возраста ta, tb, если Mt масштабируется как cMt (c - константа), то:

sup_{t∈[ta,tb]} |hc(t) - h(t)| ≤ |c-1|C(ε₁ + ε₂ + 1/pmin)

где ε₁, ε₂ ≪ 1 в среднем возрасте, доказывая слабое влияние изменений порога обнаружения.

Теорема 5.4 (нечувствительность изменений α₃):

sup_{t∈[ta,tb]} |h̃(t) - h(t)| ≤ C(ε₁ + ε₂)(||α̃₃ - α₃||∞ + ||α̃'₃ - α'₃||∞)

Даже при значительных изменениях α₃ влияние остается слабым при малых ε₁, ε₂.

3. Условия спада и восстановления (раздел 6)

Предложение 6.1 (необходимое условие спада): Для c(t) < 0 (соответствующего замедлению заболеваемости) требуется:

μ₁(t) ≥ (β₁ - α₁(2x₃(t) - 1))/(1 - x₅(t))

однако это требует экстремального увеличения μ₁, биологически необоснованного.

Предложение 6.5 (спад через α₁): При снижении α₁ на Δα в ta, tb:

cnew(t) - cbase(t) = -Δα(2x₃(t) - 1)

При x₃(t) > 1/2 умеренное снижение α₁ достаточно для c(t) < 0.

Предложение 6.7 (восстановление через β₁): При снижении β₁ на Δβ после tb существует δ > 0 такое, что:

cnew(t) > 0, ∀t ∈ (tb, tb + δ]

Ключевое понимание: требуется снижение соотношения α₁/β₁ в среднем возрасте (соответствующее снижению прогестерона), а затем восстановление этого соотношения после менопаузы (соответствующее сигналам роста из жировой ткани).

4. Доказательство доминирования c(t) над q(t) (раздел 7)

Теорема 7.5: Когда ∫W+ |f'(r)|dr ≪ 1 (малые изменения f в среднем возрасте):

sup_{t∈W+} |δh'q(t)|/|δh'c(t)| = o(1)

Схема доказательства:

• Возмущение q ограничено интегралом f' (предложение 7.2)
• Возмущение c имеет "гомогенный путь x₄", не ограниченный f' (лемма 7.3)
• Следовательно, изменения c(t) доминируют в динамике среднего возраста

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

Наборы данных

1. База данных NORDCAN: Данные заболеваемости раком молочной железы из Дании, Швеции, Финляндии, Норвегии (возрастные группы 5 лет)
2. База данных SEER: Данные заболеваемости раком молочной железы в США (годовые возрастные интервалы)
3. Выбранные когорты: Когорты рождения 1935-1939, 1940-1944, 1945-1949 (уровень скрининга <29%, снижение смещения обнаружения)

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

• Качество соответствия функции риска модели h(t) наблюдаемой заболеваемости
• Биологическая обоснованность величины изменений параметров
• Способность воспроизвести морфологические характеристики крючка Клеммесена (спад и восстановление)

Методы сравнения

• Сравнение эффектов изменения различных параметров (μ₀, μ₁, μ₂, α₁, α₂, β₁, β₂)
• Верификация того, какие изменения параметров могут производить крючковидную форму в обоснованном диапазоне

Детали реализации

1. Трехэтапная подгонка:
   • Этап I (0 до ta): Оценка всех параметров
   • Этап II (ta до tb): Фиксация других параметров, изменение только одного параметра θᵢ
   • Этап III (tb до 60 лет): Фиксация других параметров, изменение только θⱼ
2. Диапазон поиска параметров: 0, 2θᵢ (допускается снижение параметра до 0 или увеличение в 2 раза)
3. Алгоритм оптимизации: Гибридный генетический алгоритм Matlab (глобальная оптимизация)
4. Выбор окна:
   • Когорта 1935-1939: 49, 54 лет
   • Когорта 1940-1944: 46, 51 лет
   • Когорта 1945-1949: 48, 54 лет

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

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

1. Наблюдения данных (рисунок 1)

Все пять стран (Дания, Швеция, Финляндия, Норвегия, США) во всех когортах рождения демонстрируют явный крючок Клеммесена:

• Наклон заболеваемости резко снижается в окне 45-55 лет, затем восстанавливается
• Феномен демонстрирует универсальность и согласованность

2. Результаты подгонки параметров (рисунок 2)

Успешное воспроизведение крючка Клеммесена для трех когорт посредством умеренных изменений α₁ и β₁:

Характеристики кривой:

• Оранжевая часть: период роста до крючка (базовые значения α₁ и β₁)
• Синяя часть: период спада (снижение соотношения α₁/β₁)
• Зеленая часть: период восстановления (восстановление соотношения α₁/β₁)

3. Анализ сравнения параметров

Численные эксперименты подтверждают:

• α₁ и β₁: Умеренные изменения (в обоснованном биологическом диапазоне) достаточны для производства крючковидной формы
• μ₁: Требуются экстремальные изменения (выходящие за пределы обоснованного диапазона) для производства эффекта
• α₂, β₂, μ₂ (связанные с q(t)): Неспособны производить крючковидную форму в обоснованном диапазоне
• μ₀: Эффект линейный, не соответствует нелинейным характеристикам крючка

Абляционные эксперименты

Верификация вклада каждого компонента посредством теоретического анализа (разделы 6-7) и численных экспериментов:

1. Динамика обнаружения (f(t)):
   • Теоремы 5.3 и 5.4 доказывают слабое влияние изменений Mt и α₃
   • В среднем возрасте ε₁, ε₂ ≈ O(e⁻⁴⁵) ≪ 1
   • Вывод: Эффекты обнаружения не являются основной причиной крючка
2. Поздняя клональная динамика (q(t)):
   • Теорема 7.5 доказывает, что ответ на возмущения q ограничен ∫|f'|
   • Численные эксперименты показывают, что изменения α₂, β₂, μ₂ не могут производить крючковидную форму
   • Вывод: Поздняя динамика недостаточна для объяснения крючка
3. Ранняя клональная динамика (c(t)):
   • Предложения 6.5 и 6.7 обеспечивают достаточные условия
   • Численные эксперименты подтверждают эффективность умеренных изменений α₁ и β₁
   • Вывод: Баланс ранней пролиферации/апоптоза является ключевым

Экспериментальные находки

1. Необходимое количество стадий мутаций: Требуется как минимум трехстадийная модель для воспроизведения наблюдаемой кривой возраст-заболеваемость (согласуется с находками Tomasetti и др.)
2. Время изменения параметров:
   • Средний возраст (45-55 лет): α₁/β₁ должно снижаться
   • После менопаузы: α₁/β₁ должно восстанавливаться или увеличиваться
   • Временное окно высоко согласуется с периодом менопаузального перехода
3. Биологическое объяснение:
   • Снижение α₁/β₁ соответствует снижению уровня прогестерона и уменьшению менструальных циклов
   • Восстановление α₁/β₁ соответствует доминирующей роли сигналов роста из жировой ткани
4. Устойчивость: Феномен согласован в различных странах и когортах, указывая на универсальность лежащего в основе механизма

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

1. Эпидемиологические исследования

• Clemmesen (1965): Первоначальное описание крючка
• Anderson и др. (2010): Наблюдение подобного паттерна в раке молочной железы у мужчин
• Gleason и др. (2012): Стратифицированное исследование по статусу ER/PR
• Collaborative Group (2012): Крупный метаанализ, подтверждающий связь количества менструальных циклов с риском

2. Механизмы гормонов и рака молочной железы

• Coelingh Bennink и др. (2023): Предложение о том, что прогестерон (а не эстроген) является основным канцерогенным фактором
• Kim & Munster (2025): Обзор разрешающей роли эстрогена
• An и др. (2022): Молекулярный механизм, посредством которого прогестерон способствует развитию рака молочной железы через GPR126
• Antoine и др. (2016): Связь гормональной заместительной терапии с заболеваемостью

3. Математическое моделирование

• Armitage & Doll (1954): Основополагающая работа по многостадийной теории канцерогенеза
• Moolgavkar и др. (1979, 1981): Основа модели MSCE
• Brouwer и др. (2018): Расширение MSCE с временно-зависимыми параметрами
• Meza и др. (2008): Применение к анализу рака легких
• Mirzaei & Yang (2025): Модель MSCE-T (основа данной работы)

Преимущества данной работы

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

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

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

1. Механизм крючка Клеммесена: Управляется временной специфичностью изменений баланса пролиферации/апоптоза в ранних мутантных клетках (первая стадия) вокруг менопаузы
2. Биологическое соответствие:
   • Период спада: Снижение уровня прогестерона → снижение α₁/β₁ → ослабление промотирующего действия
   • Период восстановления: Сигналы роста из жировой ткани → восстановление α₁/β₁ → новые промотирующие пути
3. Исключение других механизмов:
   • Эффекты обнаружения оказывают слабое влияние (теоремы 5.3, 5.4)
   • Поздняя клональная динамика неспособна производить крючковидную форму (теорема 7.5)
   • Изменения скорости мутаций требуют необоснованных величин (предложение 6.1)
4. Минимальное количество стадий: Требуется как минимум три стадии мутаций для воспроизведения наблюдаемой кривой

Ограничения

1. Упрощения модели:
   • Параметры рассматриваются как детерминированные и кусочно-непрерывные, игнорируя стохастичность гормональных или генетических факторов
   • Не учитываются влияния микроокружения (окислительный стресс, иммунная активность, жесткость матрикса)
   • Не включена тканевая гетерогенность
2. Ограничения данных:
   • Выбор когорт с низким уровнем скрининга (<29%) для снижения смещения обнаружения, но остаточное смешивание все еще возможно
   • Возрастная группировка (5 лет) может маскировать тонкую динамику
3. Оценка параметров:
   • Стратегия кусочной подгонки может вводить граничные эффекты
   • Биологическая интерпретация параметров требует дополнительной экспериментальной верификации
4. Универсальность:
   • Модель ориентирована на рак молочной железы у женщин, применимость к мужчинам или другим гормонально-зависимым раковым заболеваниям неизвестна
   • Различные подтипы (ER+/ER-) могут иметь различную динамику

Будущие направления

1. Расширение модели:
   • Интеграция индивидуальной вариабельности гормональных траекторий
   • Включение факторов микроокружения и иммунитета
   • Учет генетической предрасположенности
2. Экспериментальная верификация:
   • Верификация посредством in vitro экспериментов влияния прогестерона на пролиферацию/апоптоз клеток на разных стадиях
   • Измерение релевантных биомаркеров до и после менопаузы
3. Клиническое применение:
   • Оптимизация стратегий скрининга на основе модели
   • Оценка профилактической эффективности гормональных вмешательств
   • Индивидуализированное предсказание риска
4. Анализ между видами рака:
   • Применение к другим гормонально-зависимым раковым заболеваниям (рак предстательной железы, рак яичников)
   • Исследование общих механизмов гормональной регуляции

Углубленная оценка

Сильные стороны

1. Методологическая строгость (★★★★★):
   • Полные математические выведения (позитивность, инвариантность, корректность)
   • Строгий анализ чувствительности (производные Гато, представление Вольтерра)
   • Четкая характеризация необходимых и достаточных условий
2. Теоретический вклад (★★★★★):
   • Первое количественное объяснение крючка Клеммесена
   • Установление математического моста между гормонами и заболеваемостью раком
   • Доказательство нечувствительности эффектов обнаружения (контринтуитивно, но важно)
3. Биологические инсайты (★★★★☆):
   • Количественные доказательства в поддержку гипотезы прогестерона
   • Идентификация ключевой действующей стадии (ранняя vs поздняя)
   • Согласованность с последними эпидемиологическими доказательствами
4. Экспериментальная верификация (★★★★☆):
   • Данные из разных стран и когорт
   • Высокое согласие между теоретическими предсказаниями и численными экспериментами
   • Изменения параметров находятся в биологически обоснованном диапазоне
5. Ясность изложения (★★★★☆):
   • Четкая структура, логичное построение
   • Подробные математические выведения
   • Четкая биологическая мотивация

Недостатки

1. Предположения модели (★★★☆☆):
   • Предположение о кусочно-непрерывности параметров может быть чрезмерно идеализировано
   • Игнорирование индивидуальной гетерогенности и стохастичности
   • Факторы микроокружения не включены
2. Дизайн экспериментов (★★★☆☆):
   • Стратегия кусочной подгонки может вводить смещения
   • Отсутствие независимого набора верификации
   • Отсутствие сравнения с другими методами моделирования
3. Биологическая верификация (★★☆☆☆):
   • Отсутствие прямых экспериментальных доказательств (клеточные эксперименты, модели на животных)
   • Биологическая интерпретация параметров опирается на косвенные выводы
   • Отсутствие различия между подтипами ER+ и ER-
4. Клиническая применимость (★★★☆☆):
   • Путь от модели к клиническому применению неясен
   • Индивидуализированное предсказание требует дополнительной работы
   • Оценка стратегий вмешательства ограничена

Влияние

1. Академическое влияние (★★★★★):
   • Заполнение долгосуществующего пробела в знаниях
   • Предоставление новой парадигмы для исследований гормонов и рака
   • Вероятное стимулирование последующих моделирующих и экспериментальных исследований
2. Практическая ценность (★★★★☆):
   • Теоретическая основа для оптимизации стратегий скрининга
   • Поддержка профилактических стратегий гормонального вмешательства
   • Содействие стратификации риска
3. Воспроизводимость (★★★★☆):
   • Полные математические выведения, поддающиеся верификации
   • Общедоступные данные
   • Частично открытый код (см. ссылку 17)
   • Полная воспроизводимость требует значительных навыков в математике и программировании

Сценарии применения

1. Прямое применение:
   • Эпидемиологические исследования рака молочной железы
   • Оценка риска рака, связанного с менопаузой
   • Анализ риск-пользы гормональной заместительной терапии
2. Расширяемое применение:
   • Другие гормонально-зависимые раковые заболевания (рак предстательной железы, рак яичников, рак эндометрия)
   • Другие заболевания, демонстрирующие немонотонные возрастные паттерны
   • Исследования жизненного пути в эпидемиологии
3. Методологические заимствования:
   • Моделирование многостадийных процессов заболевания
   • Разделение эффектов обнаружения и истинного риска
   • Применение метода интегралов Вольтерра в биологии

Ключевые ссылки

1. Clemmesen J (1965): Основополагающая работа, первоначально описавшая крючок
2. Coelingh Bennink HJT и др. (2023): Последний обзор гипотезы прогестерона как канцерогена
3. Moolgavkar SH & Knudson AG (1981): Классическая работа по модели MSCE
4. Tomasetti C и др. (2015): Оценка количества мутаций генов-драйверов
5. Mirzaei NM & Yang W (2025): Предыдущая работа по модели MSCE-T

Общая оценка: Это высококачественная междисциплинарная исследовательская работа, объединяющая строгий математический анализ с важной биомедицинской проблемой. Теоретические выведения прочны, биологические инсайты глубоки, и работа предоставляет новую количественную базу для понимания возрастного паттерна заболеваемости раком молочной железы. Хотя в экспериментальной верификации и клиническом переводе остается место для улучшения, академическая ценность и потенциальное влияние значительны. Особенно примечательно, что посредством математических методов работа четко исключает эффекты обнаружения и позднюю динамику, обеспечивая сильные количественные доказательства в поддержку гипотезы прогестерона.