흐름에 따라 달라지는 RNA 전개·재접힘, 비평형 우산 샘플링으로 밝히다

초록

본 연구는 비평형 우산 샘플링(NEUS) 기법을 이용해 명시적 흐름장 속에서 262‑뉴클레오타이드 RNA의 전개·재접힘 과정을 고해상도로 조사한다. NEUS를 병렬화하여 3.4 × 10¹³ 스텝(≈1.4 초)의 평균 첫 통과 시간을 갖는 극히 드문 전개 이벤트를 포착했으며, 흐름 속도에 따라 두 개의 경쟁적 전개 경로가 존재함을 확인하였다.

상세 분석

이 논문은 비평형 시스템에서의 장벽 극복을 목표로 하는 강화 샘플링 기법인 NEUS(Non‑Equilibrium Umbrella Sampling)를 구현하고, 이를 대규모 병렬 컴퓨팅 환경에 최적화한 점이 가장 큰 혁신이다. 기존의 평형 우산 샘플링은 자유에너지 표면을 재구성하는 데 초점을 맞추지만, 흐름과 같은 외부 구동이 존재하면 시스템은 비평형 정류 상태에 머무른다. 저자들은 이러한 비평형 상태에서도 “윈도우”를 정의하고, 각 윈도우에서 독립적인 마이크로트라젝터리를 실행한 뒤, 윈도우 간 전이 확률을 재귀적으로 업데이트함으로써 전체 전이 경로의 통계적 가중치를 정확히 추정한다.

구현상의 핵심은 (1) 각 윈도우의 초기 조건을 다른 윈도우에서 얻은 “재시작 포인트”로 설정해 샘플링 효율을 극대화하고, (2) 전이 이벤트가 발생할 때마다 전이 행렬을 즉시 갱신해 전역적인 흐름을 반영한다는 점이다. 이러한 설계는 통신 비용을 최소화하면서도 수천 개의 프로세스가 동시에 작업할 수 있게 해, 10⁶ ~ 10⁸ 단계 수준의 시뮬레이션을 실시간으로 수행할 수 있게 만든다.

모델 시스템은 262‑nt RNA를 코스코스(코스코스) 모델로 단순화하고, 주변 유체는 SRD(Stochastic Rotation Dynamics) 방식으로 명시적으로 시뮬레이션한다. SRD는 입자 기반 유체 역학을 저렴하게 구현하면서도 점성, 난류, 그리고 흐름에 의한 전단력을 자연스럽게 제공한다. 흐름 속도는 유체 입자에 일정한 평균 속도를 부여함으로써 조절했으며, 이는 RNA에 가해지는 전단 응력과 직접 연결된다.

시뮬레이션 결과, 흐름이 약할 때는 전개가 거의 일어나지 않으며, 전개가 일어나더라도 주로 내부 루프와 말단 헬릭스가 차례로 풀리는 “단일 경로”가 지배적이었다. 반면 흐름 속도가 증가하면 전단 응력이 특정 부위(예: 5′ 말단의 G‑U 쌍)를 먼저 파괴하고, 그 결과 두 개의 경쟁적 전개 경로가 나타났다. 경로 A는 말단 헬릭스가 먼저 풀리면서 전체 구조가 점진적으로 늘어나는 반면, 경로 B는 내부 스템(stem) 영역이 먼저 파괴되어 급격한 구조 변화를 일으킨다. 두 경로 모두 재접힘 단계에서 역전이 가능하지만, 경로 B는 재접힘 시 에너지 장벽이 더 높아 회복 확률이 낮다.

특히 NEUS를 이용해 평균 첫 통과 시간(MFPT)을 직접 측정했는데, 가장 느린 전개 사건의 MFPT는 약 1.4 초(3.4 × 10¹³ MD 스텝)로, 전통적인 분자 동역학으로는 실현 불가능한 시간 규모다. 이는 NEUS가 희귀 이벤트를 효율적으로 샘플링할 수 있음을 실증한다. 또한 전이 행렬을 통해 각 윈도우 간 전이 확률을 정량화함으로써, 흐름 강도에 따른 전개·재접힘 경로의 상대적 가중치를 정확히 추정했다.

결론적으로, 이 연구는 (i) 비평형 환경에서의 장벽 극복을 위한 샘플링 프레임워크를 제시하고, (ii) 흐름에 의한 전단 응력이 RNA 구조 역학에 미치는 구체적인 메커니즘을 밝히며, (iii) 대규모 병렬 시뮬레이션을 통해 실험적으로 접근하기 어려운 초희귀 전이 현상을 정량화했다는 점에서 의미가 크다. 향후 단백질·핵산 복합체, 나노채널을 통한 물질 이동 등 다양한 비평형 생물물리 시스템에 적용될 가능성이 높다.