RNA 분자의 평형 특성 및 힘에 의한 전개 경로

초록

이 논문은 Ising‑유사 모델을 이용해 단순 RNA 헤어핀과 236‑염기 길이의 Tetrahymena thermophila 리보자이머를 기계적 힘에 의해 전개시키는 과정을 분석한다. 정확한 자유에너지 지형과 상도표를 통해 헤어핀은 두 상태(접힌/펴진)로, 리보자이머는 중간 상태를 포함한 다중 단계 전개 경로를 보인다. 비평형 시뮬레이션 결과는 실험에서 관찰된 주요 전개 경로와 일치한다.

상세 분석

본 연구는 RNA의 구조역학을 정량적으로 이해하기 위해 기존의 전통적 분자동역학(MD) 접근법 대신, 이산적인 스핀 변수로 전환한 Ising‑유사 모델을 채택하였다. 각 뉴클레오티드 쌍을 ‘접힌’(스핀 = +1) 혹은 ‘펴진’(스핀 = ‑1) 상태로 표현하고, 인접한 염기쌍 사이의 상호작용을 조절함으로써 전체 사슬의 자유에너지 함수를 구성한다. 외부 힘은 각 뉴클레오티드에 작용하는 선형 포텐셜로 모델링되어, 힘‑길이 관계를 직접적으로 자유에너지에 포함시킨다. 이 모델은 해밀턴 행렬을 정확히 대각화할 수 있는 특성을 갖고 있어, 온도와 힘의 함수로서 정확한 상도표와 자유에너지 지형을 계산할 수 있다. 헤어핀에 대해서는 온도‑힘 평면에서 명확한 1차 상전이가 나타나며, 전이선은 두 개의 안정적인 최소값(접힌 상태와 완전 전개 상태) 사이를 오가게 된다. 이는 전통적인 두‑상 모델과 일치한다. 반면 236‑염기 리보자이머는 복잡한 접힘 네트워크를 가지고 있어, 자유에너지 지형에 다수의 얕은 최소값이 존재한다. 특히, 특정 구간(약 12–18 pN)에서 중간 단계가 안정화되며, 이는 실험적으로 보고된 부분 전개 구조와 부합한다. 비평형 Monte‑Carlo 시뮬레이션을 통해 전개 속도와 경로를 조사한 결과, 가장 빈번히 관찰되는 경로는 먼저 P2·P3·P4 영역이 풀리고, 이어서 P5·P6·P7이 차례로 전개되는 순서를 보인다. 이는 실험적 단일분자 힘신장 실험에서 보고된 ‘주요 전개 경로’와 일치한다. 모델의 장점은 파라미터가 적고 계산 비용이 낮아 대규모 RNA 시스템에도 적용 가능하다는 점이며, 한계는 실제 염기 간 상호작용의 비등방성 및 전해질 효과를 완전히 반영하지 못한다는 점이다. 전반적으로, Ising‑유사 모델은 RNA의 힘‑유도 전개 메커니즘을 정량적으로 예측하고, 중간 상태의 존재와 전개 경로 선택성을 설명하는 유용한 이론적 프레임워크를 제공한다.