Ion-Modulated Polyelectrolyte Complexation of DNA and Polyacrylic Acid from Molecular Dynamics Simulations

기술 요약: DNA와 폴리아크릴산의 이온 매개 폴리에lectrolyte 복합체 형성

문제 제기
동일한 전하를 가진 폴리에lectrolyte 간의 복합체 형성은 전형적인 정전기적 직관에 도전한다. 정전기적 직관은 유사한 전하를 가진 거대 분자 사이의 반발력을 예측하기 때문이다. 다가 양이온에 의한 DNA 응축과 같이 특정 시스템에서 counterion 상관관계(counterion correlations)나 전하 패치성(charge patchiness)이 인력을 매개하는 메커니즘은 이미 알려져 있으나, 화학적으로 상이한 동일 전하 고분자 간의 결합은 여전히 이해가 부족한 상태이다. 특히, 준강직한 생체 고분자(DNA)와 유연한 음이온 합성 고로머(폴리아크릴산, PAA)의 계면에서 이온의 원자 수준 조직화는 체계적으로 조사된 바 없다. 본 연구는 수용액 내에서 두 음전하 고분자 사이의 유효 인력과 가교(bridging)를 이온의 원자가(valency)와 이온 크기가 어떻게 조절하는지에 대한 이해의 공백을 다룬다.

방법론
저자들은 세 가지 서로 다른 이온 환경(단가 NaCl, 이가 MgCl₂, CaCl₂)에서 DNA–PAA 상호작용을 조사하기 위해 원자 수준의 분자 역학(MD) 시뮬레이션을 채택하였다.

• 시스템 구축: 시뮬레이션에는 24-base-pair B-DNA 이중나선(순전하 −46 e⁻)과 50% 탈양성자화된 40개 단량체의 단일 PAA 사슬(순전하 −20 e⁻)을 사용하였다. 시스템은 0.15 M의 생리적 이온 강도에서 TIP3P 물로 용매화되었다.
• 포스 필드(Force Fields): DNA는 AMBER parBSC1 포스 필드를 사용하였고, PAA는 GAFF(Generalized Amber Force Field)를 사용하였다. 이온 파라미터는 TIP3P 물에 최적화된 Joung 및 Cheatham 세트를 사용하였다.
• 시뮬레이션 프로토콜: 생산 실행(production runs)은 300 K 및 1 bar에서 NPT 앙상블로 수행되었으며, DNA 전용 시스템은 500 ns 동안, DNA–PAA 복합체는 1 μs 동안 시뮬레이션되었고, 통계적 신뢰성을 확보하기 위해 각 염 조건에 대해 세 개의 독립적인 레플리카(replica)를 수행하였다.
• 분석 프레임워크:
  • 상태 분류: 중원자 거리 컷오프(≤ 0.6 nm)와 최소 2개의 수소 결합 존재를 기준으로 "결합(bound)" 상태를 정의하는 시간 분해 기법을 사용하였으며, 일시적인 접촉과 지속적인 복합체 형성을 구분하기 위해 100 ns 이상의 지속성을 요구하였다.
  • 구조적 지표: 결합 및 비결합 상태에 따라 조건부로 RMSD, 회전 반경(Rg), 관성 주모멘트(principal moments of inertia), 비구형도(asphericity)를 계산하였다.
  • 상호작용 분석: 단거리 상호작용 에너지(Lennard-Jones 및 Coulombic)를 분해하였다. 이온 매개 접촉은 [DNA–이온–PAA] 삼중항(triplet)을 0.3 5 nm 컷오프로 식별하여 정량화하였다. 이온의 DNA 인산염 및 PAA 카르복실레이트 주변의 라디알 분포 함수(RDF)를 계산하였다.

주요 결과
시뮬레이션은 이온의 원자가와 수화 특성에 따른 결합 강도와 구조적 조직화의 명확한 계층 구조를 보여주었다:

1. 복합체 안정성 및 결합 분율:

   • CaCl₂: 가장 안정적인 복합체를 생성하였으며, 평균 결합 분율($f_{bound}$)은 0.89 ± 0.06이었다. DNA–PAA 복합체는 시뮬레이션 기간 거의 전체 동안 결합 상태를 유지하였다.
   • MgCl₂: 중간 정도의 안정성($f_{bound}$ = 0.51 ± 0.21)을 보였으며, 더 역동적이고 간헐적인 탈착을 나타냈다.
   • NaCl: 가장 높은 변동성($f_{bound}$ = 0.53 ± 0.28)을 보였으며, 결합과 비결합 상태 사이의 빈번한 전이가 발생하였고, 일부 레플리카에서는 장기간의 해리 현상이 관찰되었다.

2. 이온 배위 및 가교 메커니즘:

   • Ca²⁺: 지속적인 분자 교차 결합제(cross-linker) 역할을 하였다. RDF는 DNA 인산염 및 PAA 카르복실레이트 모두와 강한 내권(inner-sphere) 배위를 나타냈다. Ca²⁺는 결합된 프레임당 평균 2.17개의 가교를 형성하였으며, 결합 시간의 약 72% 동안 적어도 하나의 가교를 유지하였다.
   • Mg²⁺: 더 약하고 일시적인 가교(평균 1.08개 가교; 결합 시간의 약 44%)를 보였다. 단단하게 결합된 수화 껍질(hydration shell)은 외권(outer-sphere) 배위를 선호하게 하여, 안정적인 백본 고정보다는 더 역동적인 홈(groove) 관련 상호작용을 유도하였다.
   • Na⁺: 주로 정전기적 차폐를 제공하였으며, 가교 형성은 무시할 만한 수준이었다(~10%의 결합 시간). NaCl에서의 결합은 직접적인 이온 가교보다는 차폐와 약한 반데르발스 접촉에 의해 구동되었다.

3. 구조적 반응:

   • DNA: CaCl₂에서 결합된 DNA는 비결합 세그먼트(Rg ≈ 2.5 nm) 및 DNA 단독 상태(DNA alone)에 비해 더 큰 회전 반경(Rg ≈ 3.0 nm)을 보였는데, 이는 강한 교차 결합과 정렬로 인해 DNA 프레임워크가 확장되었음을 나타낸다. 반면, MgCl₂와 NaCl은 상태에 따른 DNA 구조 변화가 미미했다.
   • PAA: CaCl₂에서 PAA는 복합체 내에서 압축되어 더 응축된 형태를 취하였다.
   • 배향: Ca²⁺는 뚜렷한 방향성 정렬을 유도하였는데, DNA의 장축과 PAA의 단축이 안정적인 저각(low-angle) 구성을 채택하였다($\cos\theta$ 분포의 피크). Mg²⁺와 Na⁺는 더 유연하고 질서가 낮은 결합을 보였다.

4. 결합 모드:

   • Ca²⁺: 지속적인 백본(backbone) 결합에 의해 지배되며, 가교는 거의 전적으로 인산염 백본에 국한되었다.
   • Mg²⁺: 주요 홈(major groove) 및 부수 홈(minor groove) 내에서의 일시적 상호작용을 촉진하여, PAA가 Ca²⁺ 시스템보다 DNA 표면에 더 깊이 침투할 수 있게 하였다.
   • Na⁺: 접촉이 드물고 주로 백본에 국한되었으며 안정적인 가교는 없었다.

의의 및 주장
본 논문은 이온 특유의 성질(원자가, 수화 껍질의 가변성, 이동성)이 어떻게 동일 전하 폴리에lectrolyte의 자기 조립을 지배하는지에 대한 분자 수준의 프레임워크를 제공한다고 주장한다. 저자들은 다가 이온이 서로 대체 가능한 것이 아님을 강조한다. 즉, 가변적인 수화 껍질과 유연한 배위를 가진 Ca²⁺는 진정한 교차 결합제로서 복합체를 안정화하고 고분자 토폴로지를 재구성하는 반면, Mg²⁺는 일시적인 상호작용을 통해 중간 정도의 안정화를 제공하고, Na⁺는 차폐에 의존한다.

또한, 본 연구는 역동적인 시스템을 정확하게 특성화하기 위해서는 결합 구간과 비결합 구간을 분리하는 것이 필수적임을 강조한다. 단순한 시간 평균 지표는 강한 염 의존적 행동을 가릴 수 있기 때문이다. 나아가, 본 결과는 강한 국부적 교차 결합(Ca²⁺에서 관찰된 바와 같이)이 역설적으로 DNA 사슬의 전체 크기를 증가(확장)시킬 수 있으며, 이는 복합체를 단순히 응축으로만 해석하는 관점에 도전한다. 이러한 발견은 다가 이온 환경에서 유전자 전달, 바이오센싱 및 반응형 생체 재료 설계와 관련된 핵산-고분자 조립체를 이해하는 기초를 제공한다.