바이러스 내부 전기역학: 비특이적 전기적 상호작용이 에너지와 압력에 미치는 영향

초록

이 논문은 단일가닥 RNA 바이러스와 이중가닥 DNA 박테리오파지를 대상으로, 비특이적 전기적 상호작용을 단순화된 물리 모델로 분석한다. 캡시드 크기, 전체 전하, 그리고 유전체 전하 분포 차이에 따라 전기적 자유에너지와 내부-외부 삼투압 차이가 어떻게 달라지는지를 정량적으로 제시한다.

상세 분석

본 연구는 바이러스 내부의 전기적 상호작용을 이해하기 위해, 먼저 캡시드가 거의 구형이며 icosahedral 대칭을 따른다는 전제를 둔다. 캡시드 표면에 존재하는 양성 아미노산 잔기와 음성 인산 그룹을 각각 전하 밀도 ρ⁺, ρ⁻ 로 모델링하고, 이를 연속적인 전하 분포로 치환한다. 이후 선형화된 Poisson‑Boltzmann 방정식(Debye‑Hückel 근사)을 적용해 전위 ψ(r) 를 구하고, 전기적 자유에너지 F_el = ½∫ρ(r)ψ(r)d³r 로 표현한다.

두 바이러스군의 차이는 유전체 전하의 공간 배치에 있다. ssRNA 바이러스는 유전체가 캡시드 내부 전체에 골고루 퍼져 있어, 전하 밀도가 거의 균일한 구형 전하 구름으로 근사할 수 있다. 이 경우 전기적 에너지는 F_el ∝ Q²/(ε R) · f(κR) 로, 여기서 Q 는 전체 유전체 전하, R 은 캡시드 반경, κ 는 이온 강도에 의한 역전파 길이, f 는 차폐 효과를 나타내는 차원 없는 함수이다. 반면 dsDNA 박테리오파지는 고밀도 DNA가 캡시드 내부에 나선형으로 감겨 있어, 전하가 거의 중심축에 집중된다. 이 경우 전기적 에너지는 F_el ∝ (λ² L)/(ε) · g(κa) 로, λ는 선전하 밀도, L 은 DNA 길이, a 는 DNA의 유효 반경이며, g는 원통형 전하에 대한 차폐 함수를 나타낸다.

이러한 차이는 내부 삼투압 Π 에도 직접적인 영향을 미친다. 전기적 자유에너지의 부피 미분 Π = −∂F_el/∂V 를 적용하면, ssRNA 바이러스는 캡시드 부피가 커질수록 전기적 압력이 급격히 감소해 약 10 kPa 수준에 머문다. 반면 dsDNA 박테리오파지는 DNA가 고밀도로 압축돼 있기 때문에, 압력은 수백 kPa까지 상승하고, 이는 캡시드가 내부 압력에 의해 파열될 위험을 의미한다.

모델은 몇 가지 제한점을 가진다. 첫째, 선형화된 PB 방정식은 높은 전하 밀도나 다가이온(멀티밸런스) 존재 시 비선형 효과를 무시한다. 둘째, DNA/RNA의 구조적 강성, 겹침, 그리고 이온 특이적 결합(예: Mg²⁺)은 전기적 차폐를 강화하거나 약화시킬 수 있지만, 현재 모델에서는 평균적인 이온 강도만 고려한다. 셋째, 캡시드 내부에 존재하는 단백질 ‘tails’와 같은 국소 전하는 구형 혹은 원통형 연속 모델에 포함되지 않아, 실제 전하 분포와는 차이가 있다. 그럼에도 불구하고, 스케일링 기반 접근법은 실험적으로 측정된 캡시드 크기와 전하량을 이용해 정량적 예측을 가능하게 하며, 바이러스 설계 및 항바이러스 전략 수립에 유용한 물리적 직관을 제공한다.