Mathematical Analysis of Two Competing Cancer Cell Migration Mechanisms Driven by Interstitial Fluid Flow

Abstract

Recent experimental work has revealed that interstitial fluid flow can mobilize two types of tumor cell migration mechanisms. One is a chemotactic-driven mechanism where chemokine (chemical component) bounded to the extracellular matrix (ECM) is released and skewed in the flow direction. This leads to higher chemical concentrations downstream which the tumor cells can sense and migrate toward. The other is a mechanism where the flowing fluid imposes a stress on the tumor cells which triggers them to go in the upstream direction. Researchers have suggested that these two migration modes possibly can play a role in metastatic behavior, i.e., the process where tumor cells are able to break loose from the primary tumor and move to nearby lymphatic vessels. In Waldeland and Evje (J Biomech 81:22–35, 2018), a mathematical cell–fluid model was put forward based on a mixture theory formulation. It was demonstrated that the model was able to capture the main characteristics of the two competing migration mechanisms. The objective of the current work is to seek deeper insight into certain qualitative aspects of these competing mechanisms by means of mathematical methods. For that purpose, we propose a simpler version of the cell–fluid model mentioned above but such that the two competing migration mechanisms are retained. An initial cell distribution in a one-dimensional slab is exposed to a constant fluid flow from one end to the other, consistent with the experimental setup. Then, we explore by means of analytical estimates the long-time behavior of the two competing migration mechanisms for two different scenarios: (i) when the initial cell volume fraction is low and (ii) when the initial cell volume fraction is high. In particular, it is demonstrated in a strict mathematical sense that for a sufficiently low initial cell volume fraction, the downstream migration dominates in the sense that the solution converges to a downstream-dominated steady state as time elapses. On the other hand, with a sufficiently high initial cell volume fraction, the upstream migration mechanism is the stronger in the sense that the solution converges to an upstream-dominated steady state.