RNA 머리카락 구조와 안정성 예측 모델

초록

본 연구는 기존 3‑비드 거친 입자 모델을 개선하여, 단일 및 이가 이온 환경에서 RNA 머리카락(77 nt 이하)의 3차원 구조, 유연성 및 열안정성을 시퀀스만으로 정확히 예측한다. 새로운 암시적 전해질 전위와 서열 의존 코액시얼 스태킹 포텐셜을 도입해 실험 데이터와 비교했을 때 구조 RMSD와 용융 온도 예측 정확도가 크게 향상되었다.

상세 분석

이 논문은 RNA 3차원 구조와 동역학을 저해상도 코스그레인 모델로 다루는 접근법을 한 단계 끌어올렸다. 기존 모델은 단일 이온(예: Na⁺, K⁺)에 대한 암시적 전해질 전위를 사용했으나, 실제 세포 내 환경은 Mg²⁺와 같은 이가 이온이 중요한 역할을 한다는 점을 간과했다. 저자들은 전기적 상호작용을 Debye‑Hückel 형태에서 Poisson‑Boltzmann 기반의 비선형 전위로 확장하고, 이가 이온 농도에 따라 변하는 스크리닝 길이를 동적으로 적용하였다. 또한, 코액시얼 스태킹을 서열 의존적으로 모델링함으로써, 두 줄기 사이의 비스테이블 스택이 실제 실험에서 관찰되는 것처럼 정확히 재현되었다.

모델 검증은 세 가지 축으로 진행되었다. 첫째, 77 nt 이하의 머리카락 구조(내부 루프, 버지 루프 포함)를 실험적 NMR·X‑ray 데이터와 비교했으며, RMSD 평균이 기존 3.2 Å에서 2.1 Å로 감소했다. 둘째, 다양한 길이와 위치의 버지 루프를 가진 머리카락의 유연성을 전자 현미경 기반 FRET 실험과 비교했을 때, B‑factor와 RMSF 패턴이 높은 상관관계를 보였다. 셋째, 열안정성(용융 온도, Tₘ)을 예측하기 위해 온도 의존 자유에너지 함수를 도입했으며, Mg²⁺ 농도 0.5–5 mM 범위에서 실험값과 평균 오차가 1.8 °C에 불과했다.

특히, 이가 이온 효과를 정량화한 부분이 눈에 띈다. Mg²⁺는 전하 중화와 동시에 특정 부위의 구조를 고정시키는 ‘특이적 결합’ 효과를 갖는데, 저자들은 이를 ‘가상 결합 포텐셜’로 구현해 실제 실험에서 관찰되는 ‘Mg²⁺‑촉진 스태킹’ 현상을 재현했다. 또한, 전해질 농도에 따른 스크리닝 길이 변화를 통해 저이온(Na⁺)과 고이온(Mg²⁺) 혼합 용액에서도 모델이 일관된 예측을 제공한다는 점을 입증했다.

이러한 개선은 계산 비용을 크게 증가시키지 않으며, 전체 시뮬레이션은 수백 나노초 수준의 MD 시간에 10⁴ CPU‑hour 이하로 수행된다. 따라서 대규모 RNA 설계와 변이 효과 예측에 실용적인 도구로 활용될 수 있다.