억제자 회로의 확률적 시뮬레이션

초록

억제자 회로는 세 개의 전사인자(억제제)가 서로를 순환적으로 음성 조절하는 구조를 가진다. 이 회로는 대장균(E. coli) 플라스미드에 인공적으로 구축되었으며, 세 억제제 농도가 주기적으로 변동하는 진동 현상을 보인다. 억제제와 결합 부위가 낮은 복제 수를 갖는 경우가 많아 진동이 잡음에 의해 불규칙해지므로, 결정론적 속도 방정식만으로는 충분히 설명할 수 없다. 본 연구에서는 마스터 방정식과 몬테카를로 시뮬레이션을 이용해 억제자 회로를 확률적으로 분석하였다. 결과는 변동성이 진동이 나타나는 조건 범위와 진폭·주기를 결정론적 예측과 다르게 변화시킨다는 것을 보여준다. 모든 자유 단백질, 플라스미드, 결합된 단백질이 고복제 수에 이를 때만 두 접근법이 일치한다. 또한, 협동 결합 및 결합 억제제의 분해와 같은 미세한 메커니즘이 진동의 존재와 특성에 큰 영향을 미침을 확인하였다.

상세 요약

억제자 회로는 2000년 초반에 합성생물학 분야에서 최초로 성공적인 인공 유전 회로로 보고된 사례이며, 세 개의 전사 억제제(LacI, TetR, λ cI 등)가 서로를 억제하는 고리 구조를 가지고 있다. 이와 같은 순환 억제 메커니즘은 이론적으로는 안정적인 한계 주기(limit‑cycle) 진동을 생성하지만, 실제 세포 내에서는 각 구성 요소가 수십 개 이하의 복제 수를 갖는 경우가 많다. 복제 수가 낮을 때는 분자 간 충돌, 전사·번역의 확률적 시작·정지, 그리고 억제제와 오페레톤 결합/해리 과정에서 발생하는 내재적 잡음이 시스템의 동역학에 큰 영향을 미친다.

본 논문은 두 가지 확률적 접근법을 사용한다. 첫 번째는 마스터 방정식을 통해 모든 가능한 상태(자유 단백질, 플라스미드 복제 수, 결합된 억제제 수 등)의 전이 확률을 기술하고, 이를 수치적으로 풀어 시간에 따른 확률분포를 얻는다. 두 번째는 Gillespie 알고리즘 기반의 직접적인 몬테카를로 시뮬레이션으로, 개별 반응 이벤트를 실제 발생 순서대로 재현한다. 두 방법 모두 결정론적 ODE 모델이 예측하는 일정하고 규칙적인 진동과는 달리, 진동의 존재 여부가 파라미터(전사율, 번역율, 결합 해리 상수, 단백질 분해율 등)에 따라 크게 변동한다는 점을 보여준다. 특히, 전사·번역 효율이 낮고 단백질 분해가 빠른 조건에서는 진동이 사라지고, 대신 큰 폭의 플럭투에이션이 관찰된다. 반대로, 전사·번역 효율이 충분히 높고 분해가 느린 경우에만 비교적 규칙적인 진동이 나타나며, 이때도 주기와 진폭은 확률적 변동에 의해 평균값 주변에서 흔들린다.

흥미로운 점은 ‘협동 결합(cooperative binding)’과 ‘결합 억제제의 분해(bound‑repressor degradation)’가 시스템의 동역학에 미치는 영향이다. 협동 결합을 도입하면 억제제-오페레톤 결합이 비선형적인 Hill 계수를 갖게 되어, 억제 효과가 급격히 전이한다. 이는 결정론적 모델에서는 진동 영역을 넓히는 것으로 보이지만, 확률적 시뮬레이션에서는 오히려 잡음에 의해 진동이 끊어지는 경우가 늘어난다. 반면, 결합된 억제제가 분해될 수 있다는 가정을 추가하면, 억제제의 ‘잠재적 저장소’ 역할이 감소하고, 시스템이 더 빨리 회복하여 진동이 지속될 가능성이 높아진다. 이는 실제 세포 내에서 DNA‑결합 단백질이 복합체 형태로도 분해될 수 있다는 실험적 증거와 일치한다.

결론적으로, 이 연구는 억제자 회로와 같은 저복제수 유전 회로를 설계·해석할 때 반드시 확률적 모델을 병행해야 함을 강조한다. 결정론적 접근은 고복제수(예: 플라스미드 수가 수백 이상) 상황에서는 유용하지만, 실제 실험 조건에서는 잡음이 진동의 존재 자체를 좌우한다. 따라서 합성생물학 설계자는 파라미터 튜닝 시 ‘복제 수·분해율·협동성’과 같은 미세 조정 요소를 고려해, 원하는 동적 특성을 얻기 위한 ‘노이즈 관리’ 전략을 수립해야 한다.