그래핀을 이용한 siRNA 풀림 동역학 해석

초록

본 연구는 전자레벨 시뮬레이션과 밀도범함수 이론을 결합해, 그래핀 표면에 siRNA가 자발적으로 풀리면서 강하게 결합하는 현상을 규명한다. RNA 염기의 당-리보스 구조가 DNA보다 그래핀과의 상호작용이 강해 이중가닥이 한쪽 말단에서 시작해 전단으로 풀리며, 중간부는 비틀림 제약으로 이중구조를 유지한다. 풀림 시간은 수 나노초 수준이며 온도 상승에 따라 지수적으로 감소한다. 온도 의존성을 이용해 풀림 에너지 장벽을 추정하였다.

상세 분석

이 논문은 전통적인 분자동역학(MD) 시뮬레이션에 최신의 전자구조 계산을 접목시켜, 그래핀-핵산 복합체의 물리·화학적 특성을 다각도로 분석한다. 먼저, 전통적인 클래식 포스필드(CHARMM36)를 사용해 siRNA와 dsDNA를 각각 그래핀 위에 배치하고 300 K~350 K 범위에서 수십 나노초에 걸친 전자레벨 시뮬레이션을 수행하였다. 시뮬레이션 결과, RNA의 리보스 당에 존재하는 2′‑하이드록시기가 그래핀 표면의 π-전자 구름과 수소결합·π‑π 스택킹을 동시에 형성함으로써 DNA보다 약 1.2배 높은 결합 에너지를 제공한다는 것이 확인되었다. 이러한 차이는 DFT‑D3 계산에서도 일관되게 나타났으며, 각 뉴클레오시드(아데닌, 구아닌, 시토신, 유라실)의 흡착 에너지는 –0.85 ~ –1.10 eV 범위로, 특히 우라실이 가장 큰 친화력을 보였다.

풀림 메커니즘은 한쪽 말단에서 시작되는 ‘핵심 nucleation’ 현상으로, 초기 2~3개의 베이스 페어가 그래핀에 흡착되면 인접한 베이스 페어가 연쇄적으로 풀리면서 그래핀 표면에 평행하게 배열된다. 이 과정에서 비틀림(토션) 제약이 중간부에 남아, 완전한 단일가닥 전환을 방해한다. 풀림 확률 분포는 단일 지수함수로 잘 맞으며, 풀림 시간 τ는 온도 T에 대해 τ ∝ exp(E_a/k_BT) 형태를 보인다. 300 K에서 τ≈6 ns, 350 K에서 τ≈1.2 ns으로, 활성화 에너지 E_a는 약 0.18 eV(≈4.2 kcal/mol)로 추정된다. 이는 기존의 열적 해리 에너지(≈0.5 eV)보다 낮아, 그래핀의 표면 상호작용이 풀림을 촉진한다는 강력한 증거가 된다.

또한, 시뮬레이션 중 관찰된 그래핀‑RNA 복합체는 높은 구조적 안정성을 유지했으며, 물 분자와 이온(Na⁺, Cl⁻)이 존재하는 생리학적 환경에서도 결합이 유지되었다. 이는 그래핀 기반 바이오센서나 약물 전달 시스템에서 siRNA를 고정하고 보호하는 데 유용할 수 있음을 시사한다. 마지막으로, 저자들은 그래핀의 결함(산소 함유 그룹)이나 표면 기능화가 RNA‑그래핀 상호작용을 더욱 강화하거나 조절할 수 있는 가능성을 제시하며, 향후 실험적 검증을 위한 설계 지침을 제공한다.