다중 균주 바이러스 변이 동역학의 전역 분석

초록

본 논문은 숙주 내에서 여러 바이러스 균주가 동시에 존재하고 서로 변이할 수 있는 모델을 제시한다. 변이가 없을 때는 가장 적합한 균주만이 지속되는 전역 안정성을 Lyapunov 함수로 증명한다. 작은 변이율을 도입하면 그래프의 연결성에 따라 다른 균주의 미세한 농도가 나타나지만, 기존의 최적 균주 평형은 여전히 전역적으로 안정한다는 결과를 Smith‑Waltman의 섭동 정리를 이용해 보인다.

상세 분석

논문은 먼저 기존의 단일 균주 모델을 다중 균주 형태로 확장한다. 각 균주는 감염 세포 내 복제율, 사멸률, 면역 반응에 대한 민감도가 서로 다르며, 이를 파라미터 벡터 β_i, δ_i 등으로 표현한다. 변이 과정은 유향 그래프 G 로 모델링되며, 정점은 균주, 간선은 변이 가능성을 나타낸다. 변이율은 작은 양수 μ_{ij} 로 가정하고, 전체 시스템은 비선형 미분방정식 집합으로 기술된다. 변이가 전혀 없을 경우, 각 균주가 독립적인 동역학을 보이며, 가장 큰 기본 재생산수 R_0 (또는 가장 큰 고유값)을 가진 균주가 전역적으로 지배한다. 저자들은 이 균주의 평형점 E* 에 대해 Lyapunov 함수를 구성하여 전역 안정성을 엄밀히 증명한다.

그 다음, 변이율을 작은 파라미터 ε 로 두고 시스템을 ε‑섭동 형태로 재구성한다. 여기서 핵심은 변이 그래프 G 의 강연결성 여부에 따라 평형점의 구조가 달라진다는 점이다. 만약 G 가 강연결이면 모든 균주가 어느 정도 존재하게 되며, 평형점은 원래 최적 균주의 농도에 작은 비율로 다른 균주의 농도가 추가된 형태가 된다. 반대로 그래프가 약연결이거나 일부 정점이 고립돼 있으면 해당 균주는 완전히 사라진다.

Smith‑Waltman의 섭동 정리를 적용하기 위해 저자들은 평형점 E* 의 비퇴화성(즉, Jacobian의 고유값이 0이 아닌)과 시스템의 지속성(positively invariant set) 조건을 확인한다. 이 조건이 만족되면, 충분히 작은 ε 에 대해 원래의 전역 안정 평형이 연속적으로 변형된 새로운 평형 E_ε 으로 존재하고, 그 안정성도 보존된다. 논문은 이를 수학적으로 엄밀히 증명하고, 변이율이 충분히 작을 때 전역 안정성이 유지된다는 정리를 정리한다.

마지막으로, 저자들은 몇 가지 수치 실험을 통해 이론적 결과를 검증한다. 변이 그래프의 구조를 바꾸고 ε 값을 조정하면서 시스템이 어떻게 새로운 평형으로 수렴하는지를 시뮬레이션하고, Lyapunov 함수 값이 시간에 따라 감소하는 모습을 확인한다. 이러한 실험은 이론적 분석이 실제 동역학에 적용 가능함을 보여준다.