스트레인된 DNA의 국부 탄성, 원자 수준 시뮬레이션으로 밝히다

초록

본 연구는 AT·GC 교대 서열을 가진 14‑bp DNA에 일정한 힘과 토크를 가해 전자 수준 분자동역학 시뮬레이션을 수행하였다. 결과는 인장력이 증가하면 신축 강성이, 비틀림이 증가하면 전단·비틀림 강성이 각각 달라지며, 특히 AT 서열은 평균 트위스트가 실험값에 근접할 때 비틀림 강성이 최소가 되는 현상을 보였다.

상세 분석

이 논문은 DNA의 국부 탄성 특성이 외부 기계적 스트레스(인장·비틀림) 하에서 어떻게 변하는지를 원자 수준에서 정량화한다. 저자들은 14 bp 길이의 AT‑alternating과 GC‑alternating 두 종류의 서열을 선택했으며, 이는 한 회전 주기를 충분히 포함하면서도 계산 비용을 최소화한다. 외부 힘과 토크는 ‘steady‑load’ 방법을 이용해 특정 원자 그룹에 균등하게 분배하고, 전체 외부 힘·토크는 0이 되도록 보정하였다. 이렇게 하면 DNA 전체에 균일한 인장·비틀림 변형이 유도되며, 실제 세포 내 초코일링 상황을 모사할 수 있다.

시뮬레이션은 AMBER98 파라미터와 TIP3P 물 모델, PME 전기장 처리 등을 사용했으며, 내부 좌표 MD(ICMD) 기법으로 타임스텝을 0.01 ps까지 확대해 164 ns(≈215 만 프레임)의 장시간 데이터를 확보했다. 탄성 파라미터는 전통적인 웜‑라이크 체인(WLC) 모델의 영구 길이(l_s), 비틀림 영구 길이(l_t), 굽힘 영구 길이(l_b)로 정의하고, 각각 길이·트위스트·굽힘 각도의 분산을 통해 추정하였다.

길이 측정 방법은 세 가지(L′, L′′, L′′′)를 도입했는데, L′은 말단 기준 프레임 원점 간 직선 거리, L′′는 말단 프레임 원점을 잇는 3차원 지그재그 경로, L′′′는 3DNA 출력의 로컬 rise 합계이다. 결과적으로 L′이 가장 실험적 관측에 부합한다는 결론에 도달했으며, 이는 짧은 DNA에서 길이 변동이 주로 각도 평탄화에 의한 것임을 시사한다.

탄성 결과는 다음과 같다. 인장 강성(l_s)은 인장력뿐 아니라 비틀림 토크가 증가할 때도 상승했으며, 이는 ‘twist‑stretch coupling’이 존재함을 의미한다. 그러나 측정 방법에 따라 커플링 부호가 달라졌는데, L′′와 L′′′를 사용할 경우 실험과 반대 부호가 관찰되었다. 비틀림 강성(l_t)은 인장에 거의 영향을 받지 않지만, 토크에 따라 크게 변동했으며 AT‑서열은 평균 트위스트가 약 34°(실험값)일 때 최소값을 보였다. 이는 서열 특이적인 비틀림 탄성 최소화 현상으로, 초코일링에 민감한 전사 조절에 기여할 가능성이 있다. 굽힘 강성(l_b)은 AT‑서열에서만 비틀림·인장에 따라 변했으며, GC‑서열은 거의 일정했다. 이러한 차이는 염기쌍 간 스택 및 수소결합 패턴 차이에 기인할 것으로 추정된다.

생물학적 함의로는, 작은 토크와 인장력만으로도 특정 서열의 탄성 파라미터가 크게 변할 수 있기에, 전사인자·DNA 결합 단백질이 ‘기계적 알러스테리’를 통해 활성화·억제될 수 있음을 뒷받침한다. 또한, 실험적 초코일링 효과를 해석할 때 국부 탄성 변화를 고려해야 함을 강조한다.