암세포 집단에서 VEGF‑A 생산의 결정론적·확률론적 메커니즘과 협동 행동

초록

본 논문은 저산소증에 의해 유도되는 혈관신생 과정에서 HIF‑1α와 VEGF‑A의 상호작용을 수학적으로 모델링한다. 평균장(mean‑field) 방정식으로 결정론적 거동을 분석하고, Gillespie 알고리즘을 이용한 확률론적 시뮬레이션으로 유한 세포군에서 나타나는 협동 현상을 탐구한다. 결과는 정상, 저산소, 세포 사멸이라는 세 가지 거시적 상태를 설명하고, 잡음이 세포 간 산소 경쟁을 매개해 협동적 생존을 촉진함을 보여준다.

상세 분석

이 연구는 종양 세포 내부와 외부의 산소 흐름을 9개의 화학 종(O₂, W, Wₐ, HIF‑α, HIF‑β, VEGFᵢ, VEGF₀, E₀, Eᵢ)으로 추상화하고, 각각의 반응식을 (1)(9)로 제시한다. 평균장 접근법에서는 van Kampen 전개를 이용해 각 종의 농도를 연속적인 변수(O, W, Wₐ, Hα, Hβ, Vi, V0, Ei)로 변환하고, 8개의 상호 연결된 상미분 방정식(11)(18)을 도출한다. 보존법칙 O+W+V₀+Vᵢ+E₀+Wₐ+Eᵢ+Hα+Hβ=1을 통해 시스템 차원을 하나 줄였다. 고정점 해석에서는 세 가지 가능한 평형 상태를 찾는다. (i) 정상 상태: 산소가 충분히 공급되어 Wₐ와 Hα가 낮은 수준을 유지하고 VEGF 생산이 최소화된다. (ii) 저산소 상태: O가 감소하면서 Wₐ가 활성화되고 Hα가 축적, Hα‑Hβ 복합체 형성으로 VEGFᵢ가 증가, 외부 VEGF₀가 산소 흡수를 촉진한다. (iii) 사멸 상태: O가 거의 없고 모든 종이 빈 공간(E₀, Eᵢ)으로 전환된다. 선형 안정성 분석을 통해 각 고정점의 매력 영역을 파라미터 a, b, c, d, e, f, gα, gβ에 따라 구분한다.

확률론적 측면에서는 Gillespie 알고리즘을 적용해 유한 N(총 분자 수)에서 발생하는 인구통계적 잡음(demographic noise)을 고려한다. 시뮬레이션 결과는 평균장 예측과 달리, 작은 세포 집단에서 산소 공급이 제한될 때 서로 다른 세포가 VEGF₀를 분비함으로써 공동으로 산소를 끌어들이는 ‘협동 현상’이 나타난다. 특히, 두 번째 세포가 첫 번째 세포의 VEGF₀에 의해 산소를 획득하면 Hα 감소와 VEGF 생산 억제가 일어나, 전체 집단의 생존 확률이 비선형적으로 증가한다. 이러한 현상은 연속체(N→∞) 한계에서는 사라지며, 잡음이 시스템을 다중안정성(multistability) 영역으로 이동시키는 역할을 한다는 점에서 중요한 생물학적 의미를 가진다.

결과적으로, 모델은 저산소 환경에서 종양 세포가 어떻게 VEGF‑A를 조절해 혈관신생을 유도하고, 동시에 인구통계적 변동이 세포 간 협동을 촉진해 집단 수준의 적응을 가능하게 하는지를 정량적으로 설명한다. 이는 항‑VEGF 치료 전략을 설계할 때, 세포 간 변동성과 협동 메커니즘을 고려해야 함을 시사한다.