[논문 리뷰] Dynamics of a kinetic model describing protein exchanges in a cell population

Pierre Magal, Gaël Raoul|arXiv (Cornell University)|2015. 11. 09.

Mathematical Biology Tumor Growth참고 문헌 19인용 수 9

한 줄 요약

이 논문은 P-글리코프로테인(P-gp) 수준에 따라 구조화된 세포 집단에서 단백질 교환 동역학을 위한 운동론적 모델을 개발하며, 빈번한 교환이 발생할 경우 유체역학적 근사가 두 개의 결합된 상미분방정식으로 기술되며, 해가 유일한 평형 상태로 지수적으로 수렴함을 보여준다. 분석은 수렴성과 유체역학적 근사를 증명하기 위해 워샤르슈타인 거리 추정을 사용하며, 세포 간 단백질 이동을 통한 화학요법 내성 메커니즘을 엄밀한 수학적 프레임워크로 이해할 수 있게 한다.

ABSTRACT

We consider a cell population structured by a positive real number, describing the number of P-glycoproteins carried by the cell. We are interested in the effect of those proteins on the growth of the population: those proteins are indeed involve in the resistance of cancer cells to chemotherapy drugs. To describe this dynamics, we introduce a kinetic model. We then introduce a rigorous hydrodynamic limit, showing that if the exchanges are frequent, then the dynamics of the model can be described by a system of two coupled differential equations. Finally, we also show that the kinetic model converges to a unique limit in large times. The main idea of this analysis is to use Wasserstein distance estimates to describe the effect of the kinetic operator, combined to more classical estimates on the macroscopic quantities.

연구 동기 및 목표

• 세포 간 상호작용을 통해 암세포 간 P-글리코프로테인(P-gp) 교환의 동역학을 모델링하기 위해.
• 빈번한 교환 조건 하에서 운동론적 모델의 유체역학적 근사를 엄밀히 유도하기 위해.
• 운동론적 모델이 유일한 평형 상태로 장기적으로 수렴하는 것을 증명하기 위해.
• P-gp 수준이 세포 집단의 성장과 화학요법 내성에 미치는 영향을 분석하기 위해.
• 교환 메커니즘을 갖는 구조화된 인구 동역학을 연구하기 위해 워샤르슈타인 거리를 사용한 수학적 도구를 수립하기 위해.

제안 방법

• P-gp 농도에 따라 구조화된 세포 집단을 위한 운동론적 방정식을 수립하며, 쌍방향 교환을 모델링하는 충돌 연산자를 포함한다.
• 생물학적 일관성을 확보하기 위해 (비음수, 유계 이행) 보내는 확률의 커플러션을 통해 교환 연산자 K(x, x₁, x₂)를 유도한다.
• 워샤르슈타인 거리(W₂) 추정을 적용하여 교환 연산자의 수축 성질을 증명하고, 이로 인해 유체역학적 근사를 가능하게 한다.
• 해가 유일한 평형 상태로 지수 수렴함을 보이기 위해 W₁-워샤르슈타인 거리를 사용한다.
• 총 세포 수를 나타내는 매크로스코픽 변수 N(t)를 도입하고, 밀도 n(t,x)와 ODE 시스템을 통해 결합한다.
• 약한-* 수렴 추론과 L¹ 및 W₁ 공간 내 추정을 활용하여, 전체 해 n(t,x)가 평형 밀도와 총 세포 수의 곱으로 수렴함을 증명한다.

실험 결과

연구 질문

• RQ1 세포 간 빈번한 단백질 교환이 세포 집단 내 P-글리코프로테인 농도의 장기적 분포에 어떤 영향을 미치는가?
• RQ2 유체역학 스케일링 하에서 P-gp 교환의 운동론적 모델이 매크로스코픽 ODE 시스템으로 근사될 수 있는가?
• RQ3 운동론적 모델이 유일한 평형 상태로 수렴하는 데 필요한 조건은 무엇인가?
• RQ4 워샤르슈타인 거리 추정은 구조화된 인구 모델에서 수렴성과 유체역학적 근사를 입증하는 데 어떻게 기여하는가?
• RQ5 P-gp 수준은 세포 집단의 성장률과 내성에 어떤 역할을 하는가?

주요 결과

• 운동론적 모델은 W₁-워샤르슈타인 거리에서 유일한 평형 상태로 지수적으로 수렴하며, 수렴 속도는 성장률 함수 α(x)의 최솟값과 모델 파라미터에 의존한다.
• 빈번한 교환 조건 하에서 운동론적 모델은 총 세포 수 N(t)과 P-gp 밀도 n(t,x)에 대한 두 개의 결합된 상미분방정식으로 기술되는 유체역학적 근사를 갖는다.
• 평형 상태 해 n(t,x)는 R₊ 위의 확률측도 공간에서 유일한 극한으로 수렴하며, 수렴 속도는 r, min α(x), 그리고 ‖α′‖∞를 포함한 지수항으로 유계된다.
• 총 세포 수 N(t)는 유한한 극한 ¯N = (1/β)(r + ∫α(x)¯n(x)dx)로 수렴하며, 수렴 속도는 min(r + min α(x) − C₁‖α′‖∞, β¯N/4)이다.
• 전체 해 n(t,x)는 약한-* 수렴하여 ¯N¯n(x)로 수렴하며, 이는 W₁ 추정과 C⁰(R₊) 공간 내 근사화를 통해 증명된다.
• 분석은 워샤르슈타인 거리가 운동론적 연산자의 영향을 정량화하는 데서 중요한 역할을 함을 엄밀히 입증하며, 이는 유체역학적 근사와 장기 수렴 결과를 모두 가능하게 한다.