miRNA에 의한 유전 진동기의 확률적 동기화와 스트레스 조절 메커니즘

초록

본 연구는 miRNA가 확산 결합을 통해 두 개의 동일한 유전 진동기를 연결할 때, 동기화와 스트레스 유도라는 두 상반된 역할을 수행함을 확률적 마스터 방정식과 Gillespie 알고리즘으로 분석한다. coupling constant ε가 약 3일 때 동기화가 최강이며, 그 이상에서는 miRNA에 의한 스트레스가 동기화를 파괴한다. miRNA의 내부 합성량(k₁₉)과 분해율(k₂₀), 복합체 형성·분해 속도(k₂₁, k₂₂)도 ε에 영향을 주어 최적 조건을 제시한다.

상세 분석

이 논문은 Vilar 모델을 확장하여 miRNA(microRNA)를 포함한 유전 진동 회로를 구성하고, 두 회로를 miRNA의 확산을 매개로 한 diffusive coupling으로 연결한다. 시스템은 총 22개의 화학 종과 21개의 반응으로 정의되며, 각 반응의 확률적 발생은 Gillespie SSA를 통해 시뮬레이션된다. 핵심 변수는 coupling constant ε(=ε′)이며, 이는 miRNA가 한 세포에서 다른 세포로 이동하는 속도를 나타낸다. ε=0일 때 두 진동기는 완전히 독립적인 시간 시계열을 보이며, 2차원 재발 플롯에서 대각선과 무관한 점들이 무작위로 분포한다. ε를 증가시켜 3에 도달하면 A와 R 단백질의 농도 시계열이 서로 강하게 상관관계를 보이고, 재발 플롯의 점들이 대각선에 좁게 집중한다. 이는 stochastic noise에도 불구하고 동기화가 유지된다는 것을 의미한다. 그러나 ε를 15까지 높이면 다시 상관성이 사라지고, 점들이 대각선에서 흩어져 스트레스 효과가 동기화를 압도함을 보여준다.

무작위 coupling(ε= C·r, r은 난수) 실험에서는 C=8일 때 동기화가 발생한다. 이는 일정한 ε보다 더 큰 변동성을 필요로 함을 시사한다. miRNA 자체가 스트레스 신호로 작용할 수 있기 때문에, 확산 강도가 과도하면 세포 내 스트레스 반응이 과다 활성화되어 진동기의 자체 피드백 루프가 붕괴한다.

miRNA 합성 속도 k₁₉를 18~22 사이에서 변화시켰을 때, 최적 ε는 평균적으로 증가하지만 k₁₉가 지나치게 크면 진동 자체가 소멸한다. miRNA 분해 속도 k₂₀는 ε와 비선형 관계를 보이며, 중간값에서 ε가 최대가 된다. 복합체 C_RISC 형성 속도 k₂₁이 증가하면 ε도 상승하고, 복합체 분해 속도 k₂₂가 증가하면 ε는 감소한다. 이는 miRNA‑mRNA 상호작용의 동역학이 동기화 효율에 직접적인 조절 역할을 함을 의미한다.

결과적으로, miRNA는 두 가지 상반된 기능을 동시에 수행한다. 적절한 농도와 확산 강도에서는 세포 간 정보 전달 매개체로서 동기화를 촉진하지만, 과도한 농도나 강도에서는 스트레스 신호가 과다하게 전달되어 진동기의 정상적인 주기성을 방해한다. 이러한 경쟁 관계는 circadian clock과 같은 생체 리듬 조절에 중요한 의미를 가진다.