시스플라틴에 의한 DNA 구조 변형의 분자동역학 연구

초록

본 연구는 AMBER 힘장을 이용해 10개의 염기로 이루어진 DNA 이중나선에 시스플라틴 3′‑모노아답트를 결합한 뒤, Pt‑N7 결합이 점차 형성되는 과정을 시뮬레이션하였다. 500‑800 ps 사이에 직선형 DNA가 굽은 형태의 비스‑아답트로 전환되며, Pt와 두 번째 구아닌 N7 사이의 인력이 구조 변화를 주도함을 확인했다. 토션 상수 0.2 kcal/mol 이상이면 염기쌍 수소결합이 파괴되는 한계가 존재한다. 최종적으로 실험에서 보고된 굽힘 구조와 거의 일치하는 안정적인 모델을 얻었다.

상세 분석

이 논문은 시스플라틴이 DNA에 결합하면서 일어나는 구조적 변형을 원자 수준에서 실시간으로 추적한 최초의 분자동역학(MD) 연구이다. AMBER 5.0과 AMBER 4.1 파라미터를 기반으로, Pt(II)‑구아닌 N7 결합을 점진적으로 강화하는 ‘힘 상수 스케일링’ 기법을 도입했다. 초기 3′‑모노아답트는 Pt‑N7 거리가 약 3.6 Å로 비결합 상태였으며, 물과 Na⁺ 이온을 포함한 명시적 용매 환경에서 3500개의 TIP3P 물분자와 16개의 Na⁺가 전기 중성을 유지하도록 배치되었다.

시뮬레이션은 2 fs 타임스텝, NTP 앙상블, Berendsen 온도·압력 조절을 사용했으며, SHAKE 제약을 Pt‑N 결합 제외하고 적용했다. 초기 10‑50 ps 내에 플라티네이티드 구아닌 염기쌍은 지역적 평형에 도달했지만, 전체 DNA 사슬의 굽힘은 500 ps 이후에 급격히 진행되었다. Pt‑N7(2) 사이의 인력이 Roll과 Tilt을 증가시켜 인접 염기쌍을 기울이게 하고, 그 결과 중심부 G5‑G6 쌍의 Rise와 Roll이 크게 변동하였다.

특히, X‑N‑Pt‑N 토션 상수를 0.2 kcal·mol⁻¹ 이하로 제한해야 H‑bond가 유지되는 것이 확인되었다. 상수를 초과하면 Shear·Stretch·Stagger 값이 0.5‑1 Å까지 증가하고, H‑bond 파괴가 일어나 구조가 비정상적으로 변형된다. 이는 토션 힘이 염기쌍 사이의 수소결합을 약화시켜, Pt‑N 결합 형성 과정에서 DNA가 과도하게 왜곡되는 것을 방지한다는 중요한 메커니즘을 제시한다.

에너지 분석에서는 Pt‑N7 거리 변화에 따른 변형 에너지를 DFT 기반 Morse 포텐셜과 비교했으며, 시뮬레이션 전 구간에서 변형 에너지가 Pt‑N 인력 에너지보다 작아 실제 화학 결합 형성에 충분한 구동력을 제공함을 입증했다. 이는 조화 포텐셜을 사용하면서도 점진적 힘 상수 증가가 물리적으로 타당함을 뒷받침한다.

결과적으로 1000‑1300 ps 사이에 전체 DNA가 약 55°‑60° 정도 굽힌 안정적인 구조에 도달했으며, 이는 X‑ray와 NMR 실험에서 보고된 굽힘 각도(≈57°‑78°)와 근접한다. 전반적인 시뮬레이션은 DNA‑시스플라틴 복합체의 굽힘 메커니즘을 세 단계(구조적 자유도, Pt‑N 인력, 토션 조절)로 설명하고, 향후 다른 금속제제와 DNA 상호작용을 예측하는 모델링 기반을 제공한다.