디지털 시계 논리 모델로 보는 서카디언 시스템

초록

본 논문은 서카디언 시계의 복잡한 유전자 네트워크를 불리언 논리 모델로 단순화하여, 파라미터 공간을 유한하고 탐색 가능하게 만든다. 합성 데이터와 루시페이스 이미징 실험 데이터를 이용해 모델을 최적화하고, 위상 관계와 빛에 대한 동기화 능력을 기존 미분방정식 모델 수준으로 재현함으로써, 규제 구조 추론에도 성공하였다.

상세 분석

이 연구는 전통적인 미분방정식(ODE) 기반 서카디언 모델이 가진 파라미터 폭발 문제를 근본적으로 해결하고자 한다. ODE 모델은 연속적인 시간과 무한한 파라미터 공간을 전제로 하여, 최적화 과정에서 고성능 컴퓨팅 자원과 방대한 실험 데이터가 요구된다. 저자들은 이를 대신해 각 유전자의 발현 상태를 ‘켜짐(1)’ 혹은 ‘꺼짐(0)’으로 이산화하고, 유전자 간 상호작용을 논리 게이트(AND, OR, NOT) 형태로 표현한 불리언 네트워크를 제안한다. 이렇게 하면 시스템의 전체 상태 공간이 2^N (N은 유전자 수)으로 제한되고, 각 논리 규칙에 대한 파라미터는 진리표 형태의 유한 집합으로 축소된다.

모델 구축 단계에서는 두 가지 업데이트 스킴을 비교한다. 동기식 업데이트는 모든 노드가 동시에 이전 상태를 기반으로 새 상태를 계산하는 반면, 비동기식 업데이트는 무작위 순서로 하나씩 갱신한다. 저자는 동기식이 실험 데이터와의 비교에 더 적합함을 확인하고, 이를 기본 프레임워크로 채택한다. 파라미터 최적화는 비용 함수로서 시계열 데이터와 모델 출력 사이의 Hamming 거리 평균을 사용한다. 탐색 알고리즘으로는 전역 최적화를 위한 유전 알고리즘과, 작은 네트워크에 대해 전수 탐색을 병행한다. 특히, 파라미터 공간이 유한하기 때문에 최적해에 도달하는 데 필요한 연산량이 ODE 기반 방법에 비해 수십 배에서 수백 배까지 감소한다.

검증 단계에서는 기존에 널리 사용되는 Goodwin‑type, Kim‑Forger, 그리고 Leloup‑Goldbeter 모델에서 생성한 합성 시계열을 대상으로 불리언 모델을 피팅한다. 결과는 위상 차이, 진폭, 주기 등 주요 동역학 특성이 ODE 모델과 거의 일치함을 보여준다. 이어서 실제 생쥐 조직에서 얻은 PER2‑luciferase 발현 데이터를 적용했을 때도, 논리 모델은 정상적인 24시간 주기를 재현하고, 빛 펄스에 대한 위상 이동을 정확히 예측한다.

특히 빛에 의한 동기화 실험에서는 일조시간(daylength)을 변화시켜 다양한 광주기(LD 8:16, 12:12, 16:8 등)를 적용했을 때, 불리언 모델이 기존 ODE 모델이 제시한 복합적인 위상 응답 곡선을 거의 동일하게 재현한다. 이는 논리 모델이 비선형 피드백과 시간 지연을 충분히 포착함을 의미한다.

마지막으로, 모델 구조 탐색을 통해 실험 데이터에 가장 부합하는 규제 연결망을 자동으로 도출한다. 예를 들어, PER와 CRY 사이의 억제 관계가 필수적이며, LIGHT → PER 전이의 활성화가 핵심적인 입력임을 밝혀냈다. 이러한 구조 추론은 기존에 가설 기반으로만 제시되던 회로도를 데이터‑주도적으로 검증·수정할 수 있는 새로운 길을 연다.

전체적으로 이 논문은 불리언 논리 모델이 서카디언 시스템의 핵심 동역학을 정량적으로 재현하면서도, 파라미터 최적화와 회로 구조 추론을 효율적으로 수행할 수 있음을 입증한다. 다만, 연속적인 발현 수준을 이산화함에 따라 미세한 진폭 변동이나 온도 보상 같은 복합 현상을 완전히 포착하기는 어려우며, 이러한 한계는 하이브리드 모델(불리언+ODE)이나 다중 레벨 이산화 기법으로 보완될 수 있다.