RNA 접힘 속도 예측을 위한 접촉 대기 시간 지표

초록

본 논문은 RNA 구조 형성에 기여하는 염기 접촉의 에너지 차이를 반영한 ‘접촉‑대기‑시간’(CWT) 지표를 제안한다. 실험적으로 측정된 여러 RNA의 접힘 속도와 CWT 값을 비교한 결과, 로그 변환된 접힘 속도와의 상관계수가 –0.95로 매우 높은 음의 상관관계를 보였다. 이를 바탕으로 박테리아와 진핵의 그룹 II 인트론 접힘 속도를 예측했으며, 향후 구조 RNA 예측에도 활용 가능함을 제시한다.

상세 분석

본 연구는 RNA 2차 구조 형성 과정에서 개별 염기쌍이 차지하는 시간적 역할을 정량화하려는 시도로, 기존의 단순한 접촉 수 기반 지표가 갖는 한계를 극복하고자 했다. 저자들은 각 염기쌍이 형성되는 데 필요한 평균 대기 시간을 ‘접촉‑대기‑시간(CWT)’이라는 새로운 메트릭으로 정의하였다. 이때 대기 시간은 해당 염기쌍의 자유 에너지 변화(ΔG)를 기반으로 Boltzmann 인자를 적용해 추정했으며, GC, AU, GU 등 서로 다른 결합 유형에 대해 서로 다른 에너지 파라미터를 할당하였다. 이렇게 함으로써, 강한 결합을 형성하는 접촉은 짧은 대기 시간을, 약한 결합은 긴 대기 시간을 갖게 되어 전체 구조 형성 속도에 미치는 기여도를 보다 정밀하게 반영한다.

실험 데이터는 기존 문헌에서 보고된 10여 종의 RNA(리보자임, 티머, 스위치 등)와 그들의 접힘 속도(k_f)를 사용하였다. 각 RNA에 대해 2차 구조를 예측하고, 모든 가능한 염기쌍에 대해 ΔG 값을 계산한 뒤 CWT를 산출하였다. 결과적으로 로그 변환된 접힘 속도와 CWT 사이의 피어슨 상관계수는 –0.95(p ≪ 0.01)로, CWT가 접힘 속도를 강하게 예측한다는 것을 보여준다. 특히, 기존의 평균 접촉 수나 단순한 길이 기반 지표와 비교했을 때, CWT는 예측 정확도가 현저히 높았다.

또한 저자들은 CWT를 이용해 아직 실험적으로 측정되지 않은 박테리아와 진핵의 그룹 II 인트론 5종의 접힘 속도를 예측하였다. 예측값은 기존의 이론적 모델이 제시한 범위와 일치했으며, 특히 대형 인트론에서 CWT가 예측 오차를 크게 감소시켰다. 이는 CWT가 길이가 길고 복잡한 RNA에서도 적용 가능함을 시사한다.

한계점으로는 ΔG 파라미터가 실험 조건(이온 강도, 온도 등)에 따라 변동될 수 있다는 점과, 3차원 구조 형성 과정에서 발생하는 비평형 현상을 완전히 포착하지 못한다는 점을 들 수 있다. 향후 연구에서는 조건별 ΔG 보정과 동역학 시뮬레이션을 결합해 CWT의 정밀도를 더욱 향상시킬 필요가 있다.

전반적으로 본 논문은 RNA 접힘 동역학을 이해하고 예측하는 데 있어 에너지 기반의 시간 메트릭을 도입함으로써, 구조 생물학 및 합성 RNA 설계 분야에 새로운 도구를 제공한다는 점에서 의미가 크다.