두 유전자 모듈의 불린과 연속 동역학 비교

초록

두 개 유전자로 구성된 간단한 회로(2~3개의 조절 연결)를 대상으로, mRNA·단백질 농도의 연속 미분 방정식 모델과 불린 스위치 모델을 동시에 분석하였다. 일반화된 방법론을 이용해 다양한 연속 모델과 조절 함수 형태를 한 프레임워크 안에서 비교하고, 진동 발생 조건을 시간 스케일, 협동성 지수, 상호작용 부호와 연관시켰다. 불린 모델에서 진동이 나타나는 경우라도 연속 모델에서는 사라질 수 있고, 반대로 연속 모델만의 진동도 존재함을 확인하였다.

상세 분석

본 연구는 두 유전자로 이루어진 최소 회로를 대상으로, 전통적인 불린 모델과 연속적인 미분 방정식 모델을 동시에 적용함으로써 두 접근법 사이의 일치·불일치를 체계적으로 규명하였다. 먼저, 각 유전자는 전사·번역 과정을 통해 mRNA와 단백질 농도로 표현되며, 이들의 시간 변화는 일반적인 Hill‑type 조절 함수를 포함한 연속 미분 방정식으로 기술된다. 여기서 시간 상수(전사·번역 속도)와 Hill 계수(협동성)는 모델의 핵심 파라미터로, 이들의 비율이 진동 존재 여부를 결정한다는 점을 강조한다.

불린 모델에서는 각 유전자를 ‘켜짐(1)’ 혹은 ‘꺼짐(0)’의 이산 상태로 단순화하고, 조절 연결은 논리 게이트(AND, OR, NOT 등)로 매핑한다. 이산적 전이 규칙에 따라 시스템은 고정점 혹은 주기적 궤도(진동)를 형성한다. 저자들은 ‘일반화된 방법(generalized method)’이라 명명한 절차를 도입해, 연속 모델의 파라미터 공간을 불린 모델의 논리 구조와 직접 연결시켰다. 구체적으로, 연속 모델의 비선형성(협동성)과 시간 지연(전사·번역 지연)이 충분히 강하면, 연속 시스템은 불린 시스템이 예측한 진동과 동일한 위상 공간 구조를 보인다.

하지만 중요한 차이가 존재한다. 첫째, 연속 모델에서는 시간 스케일이 불균형할 경우(예: 한 유전자의 단백질 분해가 매우 빠른 경우) 불린 모델이 예측한 진동이 소멸하고 고정점으로 수렴한다. 둘째, 조절 함수의 부호(억제 vs. 활성)와 연결 형태에 따라, 동일한 논리 회로라도 연속 동역학에서는 비대칭적인 피드백 루프가 형성되어 새로운 진동 영역이 나타날 수 있다. 셋째, Hill 계수가 1에 가까운 경우(비협동성)에는 연속 모델이 선형에 가까워져 불린 모델이 요구하는 비선형 임계값을 충족하지 못하고, 따라서 진동이 억제된다.

결과적으로, 저자들은 “불린 진동 ⇒ 연속 진동”은 충분히 높은 협동성, 적절한 시간 상수 비율, 그리고 부호가 맞는 피드백 루프가 필요함을 수학적으로 증명했으며, 반대로 “연속 진동 ⇒ 불린 진동”은 반드시 성립하지 않으며, 연속 모델 고유의 연속적 피드백(예: 미세한 양의 억제)가 불린 모델에서는 무시되기 때문에 발생한다는 점을 강조한다. 이러한 분석은 유전자 회로 설계 시, 단순 불린 논리만으로는 충분히 예측할 수 없는 동적 특성을 이해하는 데 중요한 지침을 제공한다.