mRNA 번역의 반마코프 모델링과 단백질 생산 분포

초록

이 논문은 리보솜 이동 단계 사이의 비지수적 대기 시간을 반영한 배제 과정 모델을 제시한다. 실험 기반 파라미터로 시뮬레이션을 수행해 하나의 mRNA가 생산하는 단백질 수 E의 분포 P(E)를 구했으며, 작은 E 구간에서 기하분포와 차이가 있음을 발견했다. 또한 첫 번째 E개의 단백질 생성 시간을 이용해 P(E)를 연결시키는 해석 가능한 반마코프 모델을 제안한다.

상세 분석

본 연구는 기존의 전통적인 TASEP(전이적 비대칭 단순 배제 과정) 모델이 리보솜의 이동을 포아송 과정으로 가정하는 한계를 극복하고자 한다. 실제 실험에서 리보솜이 코돈을 읽고 이동하는 시간은 지수분포가 아니라 넓은 꼬리를 가진 분포를 보이며, 이는 번역 효율과 변동성에 중요한 영향을 미친다. 저자들은 이러한 비지수적 대기 시간을 반마코프 프로세스로 모델링하여, 각 리보솜이 다음 코돈으로 전진하기 전의 ‘대기 시간’이 임의의 확률밀도함수 ψ(t) 를 따른다고 가정한다. 이때 배제 효과는 인접한 리보솜 사이의 최소 거리(보통 10코돈)를 유지하도록 구현된다.

시뮬레이션에서는 실제 세포 내 mRNA 길이와 리보솜 입사율, 탈락율, 그리고 ψ(t) 로부터 추정된 평균 대기 시간 등을 파라미터로 사용하였다. 결과적으로 단일 mRNA가 생산하는 단백질 수 E에 대한 확률분포 P(E)를 수천 번의 시뮬레이션으로 추정했으며, 특히 E가 작을 때 P(E)가 기하분포(전통적인 마코프 모델이 예측하는 형태)와 현저히 다름을 확인했다. 이는 초기 번역 단계에서 대기 시간의 변동성이 크게 작용함을 의미한다.

이를 설명하기 위해 저자들은 ‘첫 번째 E개의 단백질이 생성될 때까지의 총 시간’ T_E 를 정의하고, T_E 의 분포를 구한 뒤 P(E) 와의 관계식을 도출하였다. 반마코프 모델에서는 각 번역 단계가 독립적인 대기 시간 ψ(t) 를 갖지만, 전체 과정은 비마코프적 특성을 유지한다. 따라서 P(E) 를 정확히 계산하려면 ψ(t) 의 Laplace 변환과 같은 수학적 도구가 필요하며, 저자들은 이를 통해 폐쇄형 해를 얻는 간소화된 모델을 제시한다. 이 모델은 실제 시뮬레이션 결과와 매우 높은 일치를 보이며, 복잡한 수치 시뮬레이션 없이도 번역 변동성을 예측할 수 있게 한다.

연구의 핵심 기여는 (1) 비지수적 대기 시간을 명시적으로 포함한 반마코프 배제 모델을 제시한 점, (2) 작은 단백질 수 구간에서 기존 마코프 모델이 과소평가하는 변동성을 정량화한 점, (3) 첫 번째 E개의 단백질 생성 시간과 전체 생산 분포를 연결하는 해석적 프레임워크를 제공한 점이다. 이러한 접근은 세포 수준에서 단백질 발현의 잡음과 조절 메커니즘을 이해하는 데 새로운 시각을 제공한다.