느린 이동이 변동 환경에서 협동형 항생제 내성을 소멸시킨다

초록

본 연구는 2차원 격자형 메타인구에서 저속 이동과 시간 변동 환경이 협동형 항생제 내성(R)과 민감성(S) 균주의 경쟁에 미치는 영향을 분석한다. 환경 변동과 인구 병목 현상이 동시에 발생할 때, 일정 범위의 비제로 이동이 내성 균주의 전멸 확률을 크게 높이며, 최적의 환경 주기와 이동률을 제시한다.

상세 분석

논문은 먼저 각 격자점(데메)을 한정된 자원(K)과 일정한 항생제 투입 하에 민감균(S)과 내성균(R) 두 종류가 경쟁하는 모델로 설정한다. R 균은 일정 수(N_th) 이상이 되면 β‑락타마제와 같은 효소를 분비해 약물을 무력화하고, 이때 S도 보호를 받아 성장률이 a만큼 감소한다. 반면 R은 대사 비용 s를 부담한다. 저자들은 Moran 과정을 기반으로 한 birth‑death 과정을 도입하고, 인구 규모에 따라 두 가지 이동 메커니즘(밀도‑의존·밀도‑비의존)을 구현한다. 환경 변동은 전역적으로 K(t)가 ‘가혹’과 ‘온화’ 상태를 주기적으로 전환함으로써 모델링되며, 이는 인구 병목을 유발한다. 핵심 분석은 (1) 이동률 m이 너무 작으면 각 데메가 독립적으로 붕괴돼 R이 국소적으로 사라지지만 전체 메타인구에서는 재이주가 부족해 회복이 어려워 전멸 가능성이 높아진다. (2) 반대로 m이 너무 크면 메타인구가 사실상 완전 혼합 상태가 되어, 기존 연구와 일치하게 R과 S가 공존한다. 저자들은 수학적 확률론적 접근과 대규모 Gillespie 시뮬레이션을 결합해, 변동 주기 τ와 이동률 m이 특정 영역(예: τ≈1030 세대, m≈10⁻³10⁻²)에서 R의 전멸 확률이 0.95 이상임을 확인한다. 특히 ‘느린‑비제로’ 이동이 병목 시점에 인구를 재분배시켜, R이 임계치 N_th 이하로 떨어지는 순간을 확대시켜 전멸을 촉진한다. 또한, 환경 변동이 R의 성장 이점을 일시적으로 억제하는 시점과 이동에 의한 재분포가 일치할 때, 전멸 시간이 최소화되는 최적 조건을 도출한다. 저자들은 1차원 격자와 다양한 K(t) 파형(사각파, 정현파)에서도 동일한 현상이 나타남을 보이며, 모델의 강건성을 검증한다.