성장 피드백과 확산이 만든 다중 안정성 암에서의 혼합 표현형과 시공간 가소성

초록

본 연구는 성장에 의한 자원 경쟁과 확산이 결합된 반응‑확산 모델을 통해, 기존의 이중 안정성 회로가 높은 자원 경쟁 상황에서 네 가지 고정점(쿼드‑안정성)으로 전이하고, 이로 인해 암 세포의 혼합 표현형과 복잡한 시공간 가소성이 발생함을 보여준다. 고정점의 상대적 안정성, 확산계수의 비대칭성, 그리고 성장 피드백 강도가 패턴 형성과 전이 속도를 조절한다는 점을 규명함으로써 전이와 종양 이질성에 대한 새로운 메커니즘적 통찰을 제공한다.

상세 분석

이 논문은 유전 회로의 비선형성, 성장에 의한 자원 경쟁(feedback), 그리고 공간 확산이라는 세 가지 핵심 요소를 통합한 수학적 프레임워크를 제시한다. 기본 회로는 전형적인 이중 안정성 스위치(positive feedback loop)이며, 세포 성장에 따라 내부 자원(예: 전사인자, ATP, 영양소)의 가용성이 감소하는 성장 피드백을 도입한다. 이 피드백은 회로의 전이율을 자원 의존적인 함수로 변환시켜, 자원 고갈 시 전이율이 급격히 감소하도록 만든다. 이러한 비선형 자원 의존성은 시스템이 특정 파라미터 구간에서 두 개의 안정 고정점 외에 추가적인 불안정 고정점들을 생성하게 하며, 결국 네 개의 고정점(두 안정, 두 불안정)으로 확장되는 쿼드‑안정성을 초래한다.

공간 차원에서는 반응‑확산 방정식을 적용해 각 세포군의 농도와 자원 농도의 시간‑공간 변화를 기술한다. 확산계수 D의 비대칭성(예: 전사인자와 자원의 확산 속도 차이)은 패턴 형성에 중요한 역할을 한다. 저확산 영역에서는 지역적 고정점이 유지되어 ‘핵’ 형태의 고정점 클러스터가 형성되고, 고확산 영역에서는 인접 클러스터 간 상호작용이 강화돼 전이 파동이 전파된다. 특히, 고정점의 비교적 안정성(라플라스 고유값의 실수부)과 확산계수의 비율이 특정 임계값을 초과하면, 전이 전파 속도가 급격히 증가하고, 이는 전이 전이(transition front)의 가속화와 유사한 현상을 만든다.

수치 시뮬레이션은 파라미터 스위프를 통해 ‘bistable → quad‑stable’ 전이를 명확히 보여준다. 자원 경쟁 강도(γ)와 성장 피드백 강도(α)를 증가시킬수록 고정점 간 거리(상태 공간에서의 차이)가 확대되고, 이는 혼합 표현형(예: 부분적으로 EMT를 획득한 세포)의 존재를 가능하게 한다. 또한, 비대칭 확산(D₁≠D₂)과 경계 조건(노이즈가 섞인 초기 조건) 변화가 패턴의 비대칭성을 유도해, 종양 내부에서 전이성 ‘핵’이 여러 개 동시 발생하는 현상을 재현한다.

이러한 결과는 암에서 관찰되는 다중 표현형(epithelial, mesenchymal, hybrid 등)과 전이 과정의 비정형적 확산을 설명하는 데 기여한다. 특히, 성장에 의해 유발되는 자원 경쟁이 회로의 다중 안정성을 촉진하고, 확산 비대칭성이 공간적 이질성을 강화한다는 점은 기존의 단순한 이중 안정성 모델이 놓친 중요한 메커니즘을 제시한다.