[논문 리뷰] Numerical study of cancer cell invasion dynamics using adaptive mesh refinement: the urokinase model

Niklas Kolbe, Jana Kat외 3명|arXiv (Cornell University)|2014. 08. 04.

Mathematical Biology Tumor Growth참고 문헌 53인용 수 3

한 줄 요약

이 논문은 고정밀도와 낮은 계산 비용을 동시에 확보하기 위해 h-형태의 메쉬 세분화를 사용하는 적응형 메쉬 세분화(AMR) 유한체적 방법을 제안한다. 이 방법은 화학성향성 및 단백분해효소 환경에서 병합 및 형성되는 세포 집합체와 같은 복잡한 역학을 정확하게 해석하면서도 안정적인 이차수 수렴성을 달성한다. 이는 기울기 기반의 세분화 및 부드러운 세분화 규칙을 도입함으로써 달성된다.

ABSTRACT

In the present work we investigate the chemotactically and proteolytically driven tissue invasion by cancer cells. The model employed is a system of advection-reaction-diffusion equations that features the role of the serine protease urokinase-type plasminogen activator. The analytical and numerical study of this system constitutes a challenge due to the merging, emerging, and travelling concentrations that the solutions exhibit. Classical numerical methods applied to this system necessitate very fine discretization grids to resolve these dynamics in an accurate way. To reduce the computational cost without sacrificing the accuracy of the solution, we apply adaptive mesh refinement techniques, in particular h-refinement. Extended numerical experiments exhibit that this approach provides with a higher order, stable, and robust numerical method for this system. We elaborate on several mesh refinement criteria and compare the results with the ones in the literature. We prove, for a simpler version of this model, $L^\infty$ bounds for the solutions, we study the stability of its conditional steady states, and conclude that it can serve as a test case for further development of mesh refinement techniques for cancer invasion simulations.

연구 동기 및 목표

• 병합 및 형성되는 농도 집합체와 같은 복잡한 역학으로 인해 미세 격자에서 암세포 침투 시뮬레이션의 계산 비용이 높아지는 문제를 해결한다.
• 반응-확산-대류 시스템에 대한 정확도를 유지하면서 계산 부담을 줄이는 강력한 고차수 수치 해법을 개발한다.
• 향후 메쉬 세분화 기법을 평가하기 위한 시험 케이스로, 로지스틱 성장이 포함된 단순화된 화학성향-접착성 이동 모델을 정당화한다.
• 완전 모델과 단순 모델의 해의 안정성과 정성적 행동을 조사하여 방법론적 검증을 지원한다.
• 적응형 메쉬 세분화와 IMEX 시간 적분을 통합하여 2차원 암세포 침투 역학의 효율적이고 확장 가능한 시뮬레이션을 가능하게 한다.

제안 방법

• 암세포 침투를 모델링하는 변동-반응-확산 시스템을 해결하기 위해 고차수의 유한체적 방법과 IMEX3 시간 적분을 사용한다.
• 셀 이분할을 통한 h-세분화를 구현하여, 해의 기울기에 기반해 계산 자원을 필요로 하는 영역에 집중한다.
• 부드러운 세분화 제약 조건을 적용하여, 이웃 셀 간의 세분화 수준 차이가 최대 1단계 이내로 제한되도록 한다. 이는 안정성과 효율성을 보장한다.
• 메쉬 세분화 및 희귀화 결정의 주요 오차 추정자로 암세포 밀도의 기울기를 사용한다.
• 완전 모델과 단순 모델 모두에 대해 해의 $L^{inity}$ 유계성을 증명하여 수치적 안정성과 수렴성을 보장한다.
• 1차원 및 2차원 실험에서 적응형 격자와 균일 격자를 비교하여 정확도 및 효율성 향상을 검증한다.

실험 결과

연구 질문

• RQ1 복잡한 농도 역학을 포함한 암세포 침투 시뮬레이션에서, 적응형 메쉬 세분화가 정확도를 유지하면서 계산 비용을 크게 줄일 수 있는가?
• RQ2 병합 및 형성되는 농도 패턴을 해석할 때, 어떤 세분화 지표(예: 암세포 기울기)가 가장 정확하고 안정적인 결과를 낳는가?
• RQ3 로지스틱 성장이 포함된 단순화된 화학성향-접착성 이동 모델이 완전한 우로키나제 모델과 유사한 정성적 역학을 보이는가? 이는 향후 테스트 케이스로 사용할 수 있는가?
• RQ4 단순 모델이 $L^{inity}$ 유계 해를 유지하는 조건는 무엇인가? 이는 수치적 안정성을 보장한다.
• RQ5 세분화 전략의 선택(예: 부드러운 세분화)이 수치 해의 수렴성과 강건성에 어떤 영향을 미치는가?

주요 결과

• 기울기 기반 세분화 및 부드러운 세분화 규칙을 적용한 적응형 메쉬 세분화 방법은 수치 실험에서 안정적인 이차수 수렴성을 달성한다.
• 균일한 미세 격자 대비 계산 비용을 크게 감소시키면서도, 병합 및 형성되는 세포 집합체와 같은 복잡한 역학을 정확하게 해석한다.
• 로지스틱 성장이 포함된 단순화된 화학성향-접착성 이동 모델은 완전한 우로키나제 모델과 유사한 일시적 행동(병합 및 형성되는 농도 집합체)을 보이며, 이는 테스트 케이스로의 타당성을 검증한다.
• 완전 모델과 단순 모델 모두에 대해 해의 분석적 $L^{inity}$ 유계성이 증명되었으며, 이는 해의 안정성 확보와 신뢰할 수 있는 수치 시뮬레이션을 가능하게 한다.
• 수치 결과에 따르면, 암세포 전구는 $t=200$ 시점에 세포외 기질을 분해하면서 비균일한 클러스터를 형성하며 확장되며, 정적인 상태는 관찰되지 않는다.
• 가장 우수한 성능을 보인 세분화 전략은 암세포 기울기를 추정자로 사용하고, 이웃 셀 간의 최대 한 단계의 세분화 수준 차이를 제한하는 방식으로, 이는 강건성과 정확성을 보장한다.