RNA 변동과 조절 RNA가 부정 피드백 루프 동역학을 좌우한다

초록

이 연구는 mRNA의 무작위 변동이 번역 과정에서 증폭되어 단일 세포 수준에서 지속적인 단백질 펄스를 유발한다는 점을 밝힌다. 또한, 비암호화 조절 RNA가 mRNA와 결합해 번역을 억제함으로써 유전자 발현 변동성을 감소시키고, 인구 평균에서도 더 안정된 동작을 보인다는 결과를 제시한다.

상세 분석

본 논문은 부정 피드백 루프(negative feedback loop, NFL)를 구성하는 유전자‑mRNA‑단백질 삼각관계를 완전한 확률론적 모델로 구현하였다. 기존 연구가 전사 단계의 잡음만을 고려하거나 평균화된 ODE 접근에 머물렀던 것과 달리, 저자들은 Gillespie 알고리즘을 이용해 전사, mRNA 분해, 번역, 단백질 분해, 그리고 피드백 억제 반응을 모두 이산적 사건으로 시뮬레이션하였다. 핵심 발견은 mRNA의 수가 적은 경우(수십 개 이하) 전사 잡음이 큰 폭으로 발생하고, 이 변동이 번역 단계에서 단백질 수량으로 비선형적으로 확대된다는 점이다. 특히, 번역 효율이 높은 상황에서는 단일 mRNA 분자가 짧은 시간 내에 수백 개의 단백질을 생산할 수 있어, 순간적인 ‘펄스’가 연속적으로 발생한다. 이러한 펄스는 개별 세포에서는 지속적인 진동 형태로 관찰되지만, 인구 수준에서는 위상이 무작위로 섞여 평균적으로는 감쇠된 파형을 보인다.

또한, 최근 비암호화 조절 RNA(예: sRNA, miRNA)의 역할이 강조되는 흐름에 맞추어, 저자들은 조절 RNA(transcript, R)와 mRNA가 결합해 복합체를 형성하고 번역을 차단하는 메커니즘을 모델에 추가하였다. 결합 해리 상수와 복합체의 분해 속도를 파라미터화함으로써, 조절 RNA가 충분히 발현될 경우 mRNA의 가용량이 감소하고, 결과적으로 단백질 펄스의 진폭과 빈도가 크게 억제된다. 시뮬레이션 결과는 조절 RNA가 존재할 때 단일 세포 내 변동성(CV)이 30~40% 감소하고, 인구 평균에서도 진동이 거의 사라져 안정적인 고정점에 수렴함을 보여준다.

이러한 결과는 ‘RNA 수준의 잡음이 단백질 수준의 동역학을 결정한다’는 기존 가설을 강화하고, 조절 RNA가 잡음 억제 장치로 기능한다는 새로운 통찰을 제공한다. 모델 파라미터 스캔을 통해 전사율, 번역율, 피드백 강도, 조절 RNA 발현량이 각각 변동성에 미치는 상대적 기여도를 정량화했으며, 특히 피드백 강도가 약할 때 조절 RNA의 억제 효과가 가장 두드러진다는 점을 강조한다.