에피제네틱 풍경으로 본 세포 재프로그래밍 메커니즘과 핵심 전사인자 탐색

초록

본 연구는 스핀 글래스 물리학을 활용해 유전체 데이터를 기반으로 에피제네틱 풍경을 정량적으로 구축한다. 각 세포 운명은 동적 어트랙터(흡인점)로 모델링되며, 외부 신호에 따라 전환이 가능하다. 고차원 풍경에서 부분 재프로그래밍 세포는 여러 운명의 유전자 발현을 동시에 보이는 하이브리드 상태로 자연스럽게 발생한다는 예측을 제시하고, 기존 데이터 재분석을 통해 이를 실증한다. 또한 알려진 재프로그래밍 프로토콜을 재현하고, 새로운 세포 운명으로의 전환에 필요한 전사인자 후보를 제시한다.

상세 분석

이 논문은 전통적인 ‘에피제네틱 풍경’ 메타포를 수학적·물리적 모델로 구체화한다는 점에서 혁신적이다. 저자들은 먼저 각 전사인자와 그 조절 관계를 이진 스핀 변수로 매핑하고, 유전체 전사 데이터와 크로마틴 접근성 정보를 이용해 스핀-글래스 해밀토니안을 구성한다. 이 해밀토니안은 다중 최소 에너지 상태를 갖는 복잡한 에너지 지형을 형성하며, 각각의 최소점은 특정 세포 운명(예: 근육, 신경, 심근 등)의 전사 네트워크를 반영한다. 동역학은 마코프 체인 혹은 랜덤 워크 형태의 시뮬레이션으로 구현되어, 외부 신호(예: 특정 전사인자 과발현) 가해 시 시스템이 기존 어트랙터에서 새로운 어트랙터로 전이하는 과정을 재현한다.

특히 고차원 공간에서 에너지 장벽이 완전히 구분되지 않아, 두 개 이상의 어트랙터 사이에 ‘중간’ 상태가 존재한다는 점을 강조한다. 이러한 중간 상태는 부분 재프로그래밍 세포(partially reprogrammed cells)로 해석되며, 실제 실험에서 관찰되는 ‘중간형’ 혹은 ‘혼합형’ 세포와 일치한다. 모델은 이들 세포가 여러 운명의 마커 유전자를 동시에 발현하는 하이브리드 특성을 가질 것이라고 예측하고, 공개된 RNA‑seq 데이터셋을 재분석해 해당 현상을 통계적으로 검증한다.

또한, 알려진 Yamanaka 요인(Oct4, Sox2, Klf4, c‑Myc)에 의한 iPSC 유도 과정을 시뮬레이션했을 때, 초기 상태에서 목표 어트랙터(iPSC)로의 전이가 에너지 장벽을 넘어야 함을 확인한다. 이 과정에서 중간 어트랙터가 일시적으로 활성화되며, 이는 실험적으로 보고된 ‘전이 단계’ 세포와 일치한다.

마지막으로, 모델이 제공하는 ‘에너지 기울기’와 ‘전이 경로’ 정보를 이용해 아직 실험적으로 검증되지 않은 세포 운명(예: 췌장 베타세포, 심장 전도계 세포)으로의 재프로그래밍에 필요한 전사인자 조합을 제안한다. 후보 전사인자는 네트워크 중심성, 에너지 감소 기여도, 그리고 기존 데이터에서의 발현 패턴을 종합해 선정된다. 이러한 접근은 전통적인 실험적 스크리닝보다 효율적인 후보 탐색을 가능하게 한다.

전반적으로 이 연구는 복잡계 물리학과 시스템 생물학을 융합해 세포 운명의 다중 안정성, 전이 메커니즘, 그리고 재프로그래밍 전략을 정량적으로 설명한다는 점에서, 향후 맞춤형 세포 치료와 인공 조직 설계에 중요한 이론적 토대를 제공한다.