🩺 Медицина · Математична онкологія
📅 Липень 2026⏱ 12 хв🟡 Середній рівень · Останнє оновлення: 3 липня 2026 р.

Ріст пухлини: аваскулярна та васкулярна стадії

Крива росту твердої пухлини — це не проста експонента, а історія у двох актах. У першому акті обмежений дифузією аваскулярний вузол росте, зупиняється і поступово порожніє зсередини. У другому — він вербує власне кровопостачання через ангіогенез, і ріст відновлюється на зовсім іншій основі. Математика, що розділяє ці два режими — логістичний та Гомпертців закони, реакційно-дифузійні градієнти кисню та ангіогенне перемикання — становить основу сучасної математичної онкології та лежить в основі комп'ютерних симуляцій пухлин, що використовуються у розробці ліків.

1. Від однієї клітини до маси

Рак починається з однієї клітини, що накопичила мутації, які вимикають нормальний контроль над проліферацією — інактивовані пухлинні супресори, такі як p53 і Rb, активовані онкогени, такі як KRAS або MYC, та уникнення апоптозу. Без стримувальних факторів ця клітина ділиться приблизно кожні 24–48 годин. Приблизно через 30 подвоєнь клон досягає приблизно 10⁹ клітин — близько одного грама, або одного кубічного сантиметра, найменшої маси, яку надійно виявляють методи візуалізації. Ще два подвоєння — і він досягає 10¹⁰–10¹¹ клітин, розміру, який зазвичай пов'язаний із клінічними симптомами або метастатичним поширенням.

Ключове розуміння математичної онкології полягає в тому, що швидкість росту не є сталою протягом цього шляху. Спочатку поживні речовини та кисень вільно дифундують із навколишньої тканини, і клітини діляться майже з максимальною біологічною швидкістю. Але дифузія має жорстку фізичну межу — приблизно 100–200 мкм, ширина кількох сотень клітин — за межами якої концентрація кисню падає нижче того, що потрібно мітохондріям. Після цієї точки ріст без власного кровопостачання стає самообмеженим.

2. Закони росту: експоненційний, логістичний, Гомпертца

Для опису розміру пухлини в часі зазвичай застосовують три закони росту, кожен з яких дедалі точніше відображає емпіричне сповільнення.

Експоненційний ріст (без обмежень): dN/dt = r·N → N(t) = N₀·e^(rt) Логістичний ріст (граничність K): dN/dt = r·N·(1 − N/K) Швидкість росту падає лінійно з наближенням N до K. Ріст за Гомпертцем (найкраще емпіричне наближення для більшості твердих пухлин): dN/dt = r·N·ln(K/N) Еквівалентно: N(t) = K·exp[ −ln(K/N₀)·e^(−rt) ] Ключова відмінність: ріст за Гомпертцем сповільнюється набагато швидше за логістичний зі збільшенням N, оскільки ln(K/N) → 0 різкіше, ніж (1 − N/K) — це дає характерну сигмоїдну криву з дуже довгим, повільно наближуваним плато, що відповідає спостережуваним даним про час подвоєння об'єму пухлини.

Емпірично час подвоєння об'єму пухлини подовжується з ростом самої пухлини — вузол 0,5 см може подвоюватися за тижні, тоді як маса 5 см тієї ж гістології подвоюється набагато повільніше. Саме таке сповільнення передбачає рівняння Гомпертца, і саме тому модель Гомпертца (запропонована Бенджаміном Гомпертцем 1825 року для актуарної смертності, пізніше застосована до кінетики пухлин Анною Кейн Лерд у 1964 році) залишається стандартним описовим законом у клінічній онкології, плануванні фракціонування променевої терапії та моделях планування хіміотерапії.

3. Аваскулярний сфероїд

Перш ніж вербувати кровоносні судини, пухлина росте як аваскулярний вузол, який у лабораторії моделюють через багатоклітинний пухлинний сфероїд (MCTS) — вільно плаваючу кульку ракових клітин, вирощену в культурі, що самоорганізується у ту саму шарувату структуру, яка спостерігається в ранніх пухлинах in vivo. Сфероїди — це робочий інструмент дослідження аваскулярної біології пухлин саме тому, що їхній ріст керується майже виключно фізикою дифузії, а не кровотоком, що робить їх математично зручними для аналізу.

Радіус сфероїда R(t) спочатку росте експоненційно, потім переходить до приблизно лінійного росту, і врешті виходить на плато при стаціонарному радіусі R∞, зазвичай 200–500 мкм — визначеному балансом між проліферувальним краєм та некротичним ядром, що не додає чистого об'єму. Це плато є прямою експериментальною ознакою дифузійної межі, описаної нижче.

4. Дифузійні межі та три зони

Транспорт кисню крізь аваскулярну тканину описується реакційно-дифузійним рівнянням. У стаціонарному стані, за приблизно сталої швидкості споживання Γ (доки кисень не вичерпається):

Стаціонарна дифузійно-споживча рівновага (сферична симетрія): D·∇²c − Γ = 0 D·(1/r²)·d/dr(r²·dc/dr) = Γ Розв'язок за умови c(R) = c₀ (концентрація на поверхні): c(r) = c₀ − (Γ/6D)·(R² − r²) Критичний гіпоксичний радіус R_diff (відстань, на яку може проникнути кисень до c → 0), за оригінальним аналізом капілярів Крога 1919 року: R_diff = √(6·D·c₀ / Γ) ≈ 100–200 мкм для O₂ у тканині D (коефіцієнт дифузії кисню) ≈ 1–2 × 10⁻⁵ см²/с

Це формує класичну тришарову структуру аваскулярної пухлини, яку видно на гістологічних зрізах необроблених вузлів і сфероїдів однаково:

Чому це обмежує ріст: Об'єм масштабується як R³, але об'єм проліферувального шару масштабується лише як R² (площа поверхні × фіксована товщина шару), коли R перевищує R_diff. Тому частка активно ділених клітин падає як 1/R, що й зумовлює сповільнення за Гомпертцем і пояснює, чому необроблений аваскулярний вузол не може перевищити приблизно 1–2 мм у діаметрі без нової судинної мережі.

5. Ангіогенне перемикання

Перехід від сплячого, обмеженого за розміром аваскулярного вузла до агресивно зростаючої васкуляризованої пухлини називають ангіогенним перемиканням — концепцію, яку сформулював Джуда Фолкман у 1970-х роках. Це зсув локального балансу між про- та антиангіогенними сигнальними молекулами, зумовлений насамперед гіпоксією в зростаючому некротичному ядрі пухлини.

Каскад сигналізації, зумовлений гіпоксією: Низький O₂ → стабілізація HIF-1α (фактор, індукований гіпоксією) → HIF-1α зв'язується з елементами відповіді на гіпоксію (HRE) → транскрипція VEGF-A ↑ (до 30 разів) → також: FGF-2, PDGF, ангіопоетин-2, IL-8 ↑ Умова перемикання (чистий ангіогенний баланс B): B = Σ(проангіогенні) − Σ(антиангіогенні) = [VEGF, FGF-2, PDGF, ...] − [тромбоспондин-1, ангіостатин, ендостатин, ...] Перемикання спрацьовує, коли B > 0 стійко утримується на межі пухлини.

За нормоксії HIF-1α безперервно гідроксилюється прогідроксилазами та розщеплюється через убіквітиновий шлях фон Гіппеля–Ліндау (VHL). За гіпоксії гідроксилювання зупиняється, HIF-1α накопичується і запускає транскрипцію VEGF-A — найважливішого проангіогенного фактора в біології пухлин та мішені препаратів, таких як бевацизумаб (Авастин). Ангіогенне перемикання зазвичай відбувається, коли вузол наближається до свого дифузійно-обмеженого плато, приблизно 1–2 мм у діаметрі — що пояснює, чому аваскулярна фаза, хоч і необхідна, зазвичай коротка й клінічно непомітна.

6. Ангіогенез: проростання нових судин

Щойно концентрація VEGF перевищує локальний поріг, сусідні ендотеліальні клітини, що вистилають наявні капіляри, реагують:

  1. Дестабілізація судини: Ангіопоетин-2 послаблює контакти перицитів і руйнує базальну мембрану навколо материнської судини, готуючи її до проростання.
  2. Вибір верхівкової клітини: Ендотеліальна клітина з найвищою сигналізацією через рецептор VEGF (VEGFR-2) стає мігруючою "верхівковою клітиною", що випускає філоподії вгору за градієнтом VEGF, через латеральне гальмування Notch/Dll4, яке пригнічує сусідні клітини від того, щоб теж стати верхівковими.
  3. Проліферація стеблових клітин: Клітини, що слідують за верхівковою клітиною, проліферують і подовжуються, формуючи порожнистий паросток, що хемотактично тягнеться до пухлини.
  4. Анастомоз і формування просвіту: Паростки зливаються із сусідніми паростками або наявними судинами, утворюючи замкнені петлі; відкривається просвіт і починається кровотік, що доставляє кисень і поживні речовини безпосередньо в масу пухлини.

Саму міграцію верхівкової клітини добре описує зміщений випадковий рух вгору за хімічним градієнтом — процес хемотаксису, змодельований у межах підходу Келлера–Сегела, тісно пов'язаного з реакційно-дифузійною математикою, що використовується для морфогенезу та бактеріального хемотаксису в інших статтях сайту:

Хемотаксис за Келлером–Сегелом (щільність ендотеліальних клітин n, VEGF c): ∂n/∂t = D_n·∇²n − ∇·(χ·n·∇c) (дифузія − хемотактичний дрейф) ∂c/∂t = D_c·∇²c + продукція(пухлина) − розпад·c χ = коефіцієнт хемотактичної чутливості; дрейфовий доданок −∇·(χ·n·∇c) тягне щільність ендотеліальних клітин угору за градієнтом VEGF, породжуючи розгалужені, деревоподібні судинні мережі, а не прямі лінії — що добре відповідає спостережуваній судинній системі пухлин.

7. Васкулярний ріст та його хаос

Після перфузії пухлина вже не обмежена дифузією, і швидкість її росту різко зростає — але побудована нею судинна мережа зовсім не схожа на здорову тканину. Судини пухлини:

Із встановленим кровопостачанням ріст уже не обмежений геометрично дифузією і в принципі може знову наблизитися до експоненційних темпів — хоча закон Гомпертца, як і раніше, добре описує загальну криву об'єму в більшості клінічних серій, оскільки нові внутрішні гіпоксичні осередки постійно відтворюють зменшену версію аваскулярної межі по всій зростаючій масі. Ця вкладена, самоподібна обмежувальна поведінка — одна з причин, чому кінетика за Гомпертцем залишається настільки стійким емпіричним наближенням для майже всіх типів твердих пухлин і видів, від мишачих ксенотрансплантатів до даних радіологічного спостереження у людини.

8. Комп'ютерне моделювання росту пухлини

Оскільки реальні пухлини не можна безперервно спостерігати у тривимірному просторі всередині живого пацієнта, комп'ютерні моделі є ключовим інструментом математичної онкології. Домінують дві родини моделей:

Мінімальне правило оновлення гібридної моделі за один крок часу, для клітини i: 1. Виміряти локальний рівень кисню c_i з дифузійного поля 2. Якщо c_i > c_prolif: клітина може поділитися (ймовірність ∝ r·Δt) 3. Інакше якщо c_i > c_death: клітина переходить у спокій (без поділу) 4. Інакше: клітина гине → видаляється, додає некротичний детрит 5. Оновити поле кисню: розв'язати ∂c/∂t = D∇²c − Γ·(живі клітини) + джерело(судини) 6. Якщо VEGF(x) > порогу поблизу судини: породити/подовжити верхівкову клітину (крок хемотаксису Келлера–Сегела, див. вище)

Ці симуляції використовують для перевірки сценаріїв "що якщо", які були б неетичними або неможливими у пацієнтів — наприклад, щоб передбачити, як змінюється крива росту пухлини за різних графіків фракціонування променевої терапії, або як виникають резистентні субклони і просторово розділяються під тиском хіміотерапії, задовго до того, як пацієнту буде введено першу дозу препарату в клінічному дослідженні.

9. Терапевтичні наслідки

Кожен етап цієї історії росту є мішенню для ліків, які активно застосовують у клініці сьогодні:

Розуміння росту пухлини як двофазного процесу "дифузія-потім-перфузія" — а не єдиної експоненційної кривої — дозволяє онкологам та обчислювальним біологам кількісно оцінювати, коли пухлина, ймовірно, прискориться, і в якому місці її біології конкретна терапія насправді втручається.