압력이 작은 RNA의 돌출 염기 동역학에 미치는 영향

초록

분자 동역학 시뮬레이션을 통해 압력이 증가함에 따라 RNA의 비짝염기(돌출 염기) 루프 형성 자유에너지 장벽이 단조히 상승하고, 평균 첫 통과 시간(MFPT)은 2 kbar 이하에서는 서서히 늘어나지만 2 kbar를 초과하면 급격히 증가한다는 결과를 제시한다. 물의 수화층 변화가 이러한 현상에 기여함을 논의한다.

상세 요약

본 연구는 작은 RNA 구조에서 비짝염기(bulge base)가 루프 형태로 튀어나오는 전이 과정을 압력 변수와 연계시켜 정량적으로 분석하였다. 저자들은 AMBER 기반의 원자 수준 분자 동역학(MD) 시뮬레이션을 수행하고, 각 압력 조건(1 kbar, 2 kbar, 3 kbar, 4 kbar)에서 비짝염기의 회전 각도와 거리 좌표를 반응 좌표(reaction coordinate)로 정의하였다. 자유에너지 표면(FES)은 메타다이나믹스와 WHAM(Weighted Histogram Analysis Method)을 이용해 재구성했으며, 루프 형성에 해당하는 자유에너지 장벽(ΔG‡)이 압력 증가에 따라 1 kcal/mol 단위로 꾸준히 상승함을 확인하였다. 이는 고압 하에서 물 분자들이 RNA의 주 골격을 더욱 촘촘히 둘러싸면서 수소 결합 네트워크가 강화되고, 비짝염기의 탈출을 억제하는 메커니즘으로 해석된다.

MFPT 계산은 1차원 확산 방정식의 첫 통과 시간 해석식을 적용했으며, 확산 계수(D)는 압력에 따라 감소하는 경향을 보였다. 2 kbar 이하에서는 ΔG‡와 D의 변화가 완만해 MFPT는 약 10~30 ns 수준에서 서서히 증가했지만, 2 kbar를 초과하는 구간에서는 ΔG‡가 급격히 상승하고 D가 현저히 감소하면서 MFPT가 수백 나노초에서 마이크로초 수준으로 폭발적으로 늘어났다. 이러한 비선형 응답은 압력에 민감한 수화층 구조 전이가 임계점 역할을 함을 시사한다.

수화층 분석에서는 고압에서 첫 번째 수화 껍질이 더 밀집되고, 물 분자들의 재배열이 비짝염기의 접근 경로를 차단한다는 RDF(radial distribution function) 결과를 제시한다. 또한, 물의 전기적 특성(극성, 유전율) 변화가 RNA의 전하 분포에 미치는 영향을 고려해, 고압이 전반적인 전기적 스크리닝 효과를 강화시켜 비짝염기의 전자 구름을 안정화시키는 것으로 해석한다.

결과적으로, 압력은 RNA 구조 역학에 다중 단계 영향을 미치며, 특히 2 kbar라는 임계 압력 근처에서 수화층 재구성에 따른 자유에너지 장벽과 확산 동역학이 급변한다는 점이 핵심 통찰이다. 이러한 발견은 고압 환경(심해, 고압 실험실)에서 RNA‑리간드 상호작용을 설계하거나, 압력 감응형 바이오센서 개발에 중요한 기초 정보를 제공한다.