Next Article in Journal
Sub-Diffusion Two-Temperature Model and Accurate Numerical Scheme for Heat Conduction Induced by Ultrashort-Pulsed Laser Heating
Previous Article in Journal
Numerical Investigation of the Three-Dimensional HCIR Partial Differential Equation Utilizing a New Localized RBF-FD Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Migration and Proliferation Dichotomy: A Persistent Random Walk of Cancer Cells

Department of Mathematics, College of Science, Sultan Qaboos University, P.O. Box 36, Al Khodh, Muscat 123, Oman
Fractal Fract. 2023, 7(4), 318; https://doi.org/10.3390/fractalfract7040318
Submission received: 27 February 2023 / Revised: 2 April 2023 / Accepted: 4 April 2023 / Published: 7 April 2023
(This article belongs to the Section Life Science, Biophysics)

Abstract

:
A non-Markovian model of tumor cell invasion with finite velocity is proposed to describe the proliferation and migration dichotomy of cancer cells. The model considers transitions with age-dependent switching rates between three states: moving tumor cells in the positive direction, moving tumor cells in the negative direction, and resting tumor cells. The first two states correspond to a migratory phenotype, while the third state represents a proliferative phenotype. Proliferation is modeled using a logistic growth equation. The transport of tumor cells is described by a persistent random walk with general residence time distributions. The nonlinear master equations describing the average densities of cancer cells for each of the three states are derived. The present work also includes the analysis of models involving power law distributed random time, highlighting the dominance of the Mittag–Leffler rest state, resulting in subdiffusive behavior.

1. Introduction

Tumor cell invasion is a complex process that involves both the proliferation and migration of cancer cells. Proliferation refers to the ability of cancer cells to divide and grow, while migration refers to their ability to move and invade surrounding tissues. Researchers have devoted large efforts in investigating the development of malignant tumors, particularly those that are characterized by dysregulated cell migration and uncontrolled proliferation. Gliomas, a type of brain tumor, are known for their high rates of cell proliferation and motility, as well as their ability to invade normal tissue away from the multicellular tumor. Clinical studies and investigations have indeed shown that the proliferation rate of migratory cells in gliomas is often lower in the invading tumor tissue compared to the tumor core [1,2]. This relationship between cell migration and proliferation, known as the migration–proliferation dichotomy, has been shown to play a critical role in the invasiveness of malignant gliomas and is an important factor that must be taken into consideration when developing treatments for these tumors. Despite advances in medical technology and research, the switching between two phenotypes, as well as the invasion process, remain not fully understood. This presents a challenge in developing effective treatments and accurate prognoses for patients with this type of tumor. The mathematical modeling of the migration and proliferation dichotomy of cancer cells has been an active area of research in recent years, as it provides a valuable tool for understanding the complex interplay between these two processes in cancer progression and invasion. The migration–proliferation dichotomy in gliomas has been studied using a variety of approaches. These include a stochastic two-state switching model, where the migration–proliferation of cancer cells was formulated in terms of a continuous time random walk (CTRW) [3,4,5,6,7]; multiparametric modeling of the phenotype switching [8]; the simulation of multiscale glioma growth and invasion using an agent-based approach [9,10]; and numerical and analytical approaches, where the phenotype switching depends on oxygen in a threshold manner, which were developed in the framework of reaction–diffusion equations [11]. It was proposed by Khain et al. [12,13] that the motility of cancer cells is a function of their density. The phenomenon of phenotypic switching caused by the density effect was also discussed in [14,15]. In recent years, there has been a significant amount of research focused on the “Go-or-Grow” mechanism and its effect on glioma invasion. The role of the “Go-or-Grow” mechanism in glioma invasion, which is thought to be linked to oxygen availability in the environment of a growing tumor, was proposed by Hatzikirou et al. [16]. Multiscale modeling for acid- and vasculature-mediated glioma invasion and the effects of go-or-grow dichotomy and tissue anisotropyin were discussed in [17]. Mansury and colleagues [18] introduced the use of evolutionary game theory in an agent-based brain tumor model as a means of studying the evolution and dynamics of brain tumors. An alternative approach to studying the growth and progression of tumors was proposed by Böttger et al. [19], in which a model that specifically incorporates tumor cell dynamics was developed to analyze the interplay of proliferation and migration processes in detail. Recently, Godlweski et al. [20] identified a single microRNA that plays a key role in the regulation of glioma cell proliferation, migration, and responsiveness to glucose deprivation. New approaches to modeling glioma invasion have emerged in recent years, with sources [21,22,23] providing examples.
Previous research by Fedotov et al. [24] demonstrated that the use of the standard diffusion approximation for transport processes, in conjunction with logistic growth, results in an overestimation of the overall propagation rate in glioma. This finding was verified in [5,6]. The presence of anomalous diffusion in the stochastic movement of cells is frequently observed, thereby highlighting the significance of fractional diffusion equations and fractional derivative in modeling some biological phenomena, such as the cancer cell motility [25] and tumor growth [26]. The persistent random walk model constitutes a well-established framework that can be easily augmented to incorporate extensions, such as cell reactions [27,28]. However, most existing methods and techniques for addressing the issue of persistent random walks with reactions focus on Markovian switching between two states [29,30,31]. According to a study described in [32], cell motility can be explained by a persistent random walk model, which demonstrates that cells have a memory of their past velocities.
In [33], the authors proposed a model that incorporates a three-state switching mechanism and generates anomalous subdiffusion. The resulting governing equations for this formulation are
p + t + ν p + x = λ p + ( x , t ) + r + i 0 ( x , t ) ,
p t ν p x = λ p ( x , t ) + r i 0 ( x , t ) ,
p 0 t = λ p + ( x , t ) + λ p ( x , t ) ( 1 r 0 ) i 0 ( x , t ) ,
where p + ( x , t ) and p ( x , t ) are the probability density functions (PDFs) of particles moving with constant speed ν in the positive or negative direction for exponentially distributed running times with rate λ , and p 0 ( x , t ) is the PDF of a particle that has zero velocity for Mittag–Leffler distributed resting times. Here, r + , r , and r 0 describe three possible transitions that the particle can make from the rest, and i 0 ( x , t ) is the integral escape rate from the rest state and is defined as
i 0 ( x , t ) = τ 0 β D t 1 β p 0 ( x , t ) ,
where D t 1 β denotes the Riemann–Liouville derivative and is given by
D t 1 β p 0 ( x , t ) = 1 Γ ( β ) t 0 t p 0 ( x , t ) ( t t ) 1 β d t , ( 0 < β < 1 ) .
It is important to highlight that for non-Markovian switching states in random walk theory, a general formula for the escape rate i 0 ( x , t ) can be represented through convolution of the memory kernel and the density, which will be illustrated later.
Although the model (1)–(3) involves a non-Markovian anomalous rest state, the switching running state is Markovian. The main objective of this paper is to expand upon the work of the anomalous stochastic transport model of particles outlined in [33] by incorporating non-Markovian running states. Specifically, we will assume that the running time is arbitrarily distributed [34], rather than solely following an exponential distribution as described in [33]. A second extension of [33] is the implementation of a reaction term into the non-Markovian model. Another main objective of this paper is to extend the analysis of the migration–proliferation dichotomy model presented in [6]. For a biologically motivated reaction–diffusion model, it may be worthwhile to enhance the model by incorporating additional states for the cells or the particles. Our approach addresses this by by including an additional state for migrating cells, in addition to the existing migratory and proliferative states. We are interested in investigating the impact of considering the movement of cells in both positive and negative directions with finite velocities ± ν , as opposed to the movement of cells in a single direction with a fixed velocity, as previously considered in [5,6]. Our goal is to create a more realistic theory for the phenomenon of phenotype cell switching. Our model provides a description of the complex cell transport in the framework of a persistent random walk with memory effects, taking into account anomalous transport and the resulting slow mobility of cancer cells.

2. Three-State Proliferation–Migration Model Involving a Residence Time Variable

In this paper, we provide an analysis of the migration and proliferation of tumor cells using a three-state persistent random walk model on a one-dimensional space with continuous sample paths (no jumps). The model considers transitions between moving states with constant velocity (migratory phenotype) via an intermediate resting state with zero velocity (proliferating phenotype). More specifically, we suggest a model in which the process of tumor cell invasion consists of three possible states: two moving states, in which the tumor cells randomly move (migrate) with a constant velocity ν in the positive (state +) or negative (state −) direction but are unable to proliferate and can switch phenotypes; and a resting state (state 0), in which the tumor cells are allowed to proliferate but do not migrate and can also switch phenotypes.
One way to characterize the migration–proliferation dichotomy is through the use of a commonly used phenomenological model that employs reaction–diffusion equations. The basic structure of our model is as follows: The migratory phenotype cell remains in state + during a running time T + (state +’s residence time) or in state − during a running time T (state −’s residence time). In these moving states, the cell starts to move with constant speed ν at time t = 0 in the positive (state +) or negative (state −) direction, and after a random time (running time) it switches to a proliferating phenotype cell. Following a waiting time T 0 (state 0’s residence time) spent in state 0, the cell makes the choice to switch to some next state. Specifically, it either returns to a migratory phenotype (state + or state −) or remains at its resting state (state 0). Clearly, the fundamental aspect of this process can be represented by a three-state random process, where the cell exhibits one of three distinct states: positive directional movement of tumor cells, negative directional movement of tumor cells, or a state of rest for tumor cells. We assume that the residence times are random variables.
To describe the transition process between these states, we introduce the switching rate, β j ( τ ) , as a function of residence time τ such that [35]
β j ( τ ) = ψ j ( τ ) Ψ j ( τ ) , j = + , , 0 ,
where ψ j ( τ ) and Ψ j ( τ ) , j = + , , 0 are PDFs of the residence times and the corresponding survival functions, respectively. Ψ j ( τ ) are defined as
Ψ j ( τ ) = τ ψ j ( u ) d u , j = + , , 0 .
In this case, β + ( τ ) h and β ( τ ) h represent the conditional probability of transition from migrating states to a proliferating state in the small interval ( τ , τ + h ) , given that there is no transition up to time τ . The product β 0 ( τ ) h has a similar meaning for the transition from a proliferating state to one of migrating states. Thus, β + ( τ ) is the switching rate from state + to state 0, β ( τ ) is the switching rate from state − to state 0, and β 0 ( τ ) determines the switching rate from state 0 to either state + or state −. Since the switching rate β j ( τ ) is dependent on the duration of time the cell has stayed in a specific state, our approach has the advantage of being adaptable to various forms of the residence time probability distribution. This is noteworthy due to recent evidence indicating that the PDF of the residence time may not be exponential. In this study, we introduce three PDFs: p + ( x , t ) , p ( x , t ) and p 0 ( x , t ) . The functions p + ( x , t ) and p ( x , t ) represent the probability densities that the position of cells with a migratory phenotype (moving state) is within the interval ( x , x + d x ) at time t and moving with positive and negative velocities, respectively. The function p 0 ( x , t ) represents the probability density that the position of cells with a proliferating phenotype (resting state) is within the interval ( x , x + d x ) at time t and has zero velocity.
In this section, we aim to develop a migration and proliferation model for cancer cells and derive non-Markovian master equations to describe the dynamics of the cells transport process. Transition probabilities β + ( τ ) , β ( τ ) and β 0 ( τ ) can be introduced that depend on the residence time τ spent in the migrating and proliferating states, respectively [35]. It is useful to define the structured densities of cancer cells as dependent on the residence time τ [36]. Let n + ( x , t , τ ) and n ( x , t , τ ) be the densities of migrating cells at point x at time t, whose residence time in a migratory phenotype state (state + or state −, respectively) lies in the interval ( τ , τ + d τ ) . The density of cells in a proliferating state that corresponds to the aforementioned migrating cells is n 0 ( x , t , τ ) .
The balance equations for structured densities n + ( x , t , τ ) , n ( x , t , τ ) and n 0 ( x , t , τ ) can be written as follows.
For migrating states,
n + t + n + τ + ν n + x = β + ( τ ) n + ,
n t + n τ ν n x = β ( τ ) n .
For a proliferating state,
n 0 t + n 0 τ = β 0 ( τ ) n 0 + f ( p ) n 0 .
The nonlinear function f ( p ) is the proliferation rate, where p = p + + p + p 0 . For example, we can use the logistic growth for cell proliferation involving the cell proliferation rate U and the carrying capacity K such that [6]
f ( p ) = U 1 p K .
We assume that the initial conditions for the system are given by
n j ( x , 0 , τ ) = p j 0 ( x ) δ ( τ ) , j = + , , 0 ,
where p j 0 ( x ) is the initial densities of cancer cells ( 0 τ t ) . It corresponds to the case when the residence time of all cells at t = 0 equals to zero.
The boundary conditions at zero running/waiting time ( τ = 0 ) are
n + ( x , t , 0 ) = α + 0 t β 0 ( τ ) n 0 ( x , t , τ ) d τ ,
n ( x , t , 0 ) = α 0 t β 0 ( τ ) n 0 ( x , t , τ ) d τ ,
n 0 ( x , t , 0 ) = 0 t β + ( τ ) n + ( x , t , τ ) d τ + 0 t β ( τ ) n ( x , t , τ ) d τ + α 0 0 t β 0 ( τ ) n 0 ( x , t , τ ) d τ ,
where α + , α and α 0 describe the conditional transition probabilities from proliferation state to a state of positive velocity, negative velocity, or remaining at its resting state, respectively. Specifically, α + represents the probability of transitioning to a state of positive velocity, ν , while α represents the probability of transitioning to a state of negative velocity, ν . On the other hand, α 0 represents the transition probability that the resting cell remains at rest again, where α + + α + α 0 = 1 . Generally, the probabilities α + , α and α 0 can facilitate the formulation of the effect of the external force.
The mean densities at point x at time t, p j ( x , t ) , can be obtained by integrating the structured densities n j ( x , t , τ ) over residence time variable τ
p j ( x , t ) = 0 t n j ( x , t , τ ) d τ , j = + , , 0 .

3. Non-Markovian Three-State Model

The aim of this section is to set up a non-Markovian model for the migration and proliferation of cells by eliminating the residence time variable τ and find equations for p + ( x , t ) , p ( x , t ) and p 0 ( x , t ) by solving the partial differential Equations (8)–(10) together with the boundary conditions (13)–(15) at τ = 0 and initial condition (12) at t = 0 .
The product β j ( τ ) n j ( x , t , τ ) gives the phenotype switching rate corresponding to a particular residence time τ . If we denote the fluxes between migrating and proliferating cells by i + ( x , t ) , i ( x , t ) and i 0 ( x , t ) , then the switching term i j ( x , t ) can be obtained by integrating β j n j over variable τ from 0 to t:
i j ( x , t ) = 0 t β j ( τ ) n j ( x , t , τ ) d τ , j = + , , 0 .
It follows from Equations (13)–(15) and (17) that
n + ( x , t , 0 ) = α + i 0 ( x , t ) , n ( x , t , 0 ) = α i 0 ( x , t ) , n 0 ( x , t , 0 ) = i + ( x , t ) + i ( x , t ) + α 0 i 0 ( x , t ) .
Applying the method of characteristics, we obtain the following solutions to Equations (8)–(10):
For migrating states ( 0 τ < t )
n + ( x , t , τ ) = n + ( x v τ , t τ , 0 ) e 0 τ β + ( u ) d u ,
n ( x , t , τ ) = n ( x + v τ , t τ , 0 ) e 0 τ β ( u ) d u .
For a proliferating state ( 0 τ < t )
n 0 ( x , t , τ ) = n 0 ( x , t τ , 0 ) e 0 τ β 0 ( u ) d u + t τ t f ( p ( x , u ) ) d u .
It can be observed that all of the derived solutions (19)–(21) contain an exponential factor e 0 τ β j ( u ) d u that can be understood to represent the survival function Ψ j ( τ ) :
Ψ j ( τ ) = e 0 τ β j ( u ) d u , j = + , , 0 .
Considering the boundary conditions specified in Equations (13)–(15), the solutions represented by Equations (19)–(21) can be expressed in terms of the survival function from (22) and the fluxes between the migrating and proliferating states (switching terms) from (17) as follows:
n + ( x , t , τ ) = α + i 0 ( x ν τ , t τ ) Ψ + ( τ ) ,
n ( x , t , τ ) = α i 0 ( x + ν τ , t τ ) Ψ ( τ ) ,
n 0 ( x , t , τ ) = i + ( x , t τ ) + i ( x , t τ ) + α 0 i 0 ( x , t τ ) Ψ 0 ( τ ) e t τ t f ( p ( x , u ) ) d u .
It should be noted that the residence time PDF, ψ j ( τ ) , can be expressed in terms of the switching rate, β j ( τ ) , as follows [35]:
ψ j ( τ ) = β j ( τ ) e 0 τ β j ( u ) d u , j = + , , 0 .
The nonlinear master equations for unstructured densities p + ( x , t ) , p ( x , t ) and p 0 ( x , t ) can be obtained by differentiating (16) with respect to time t together with (8)–(10):
p + t = n + ( x , t , t ) 0 t β + ( τ ) n + ( x , t , τ ) d τ ν 0 t n + x d τ 0 t n + τ d τ ,
p t = n ( x , t , t ) 0 t β ( τ ) n ( x , t , τ ) d τ + ν 0 t n x d τ 0 t n τ d τ ,
p 0 t = n 0 ( x , t , t ) 0 t β 0 ( τ ) n 0 ( x , t , τ ) d τ + 0 t f ( p ) n 0 ( x , t , τ ) d τ 0 t n 0 τ d τ .
By using (13)–(17), we obtain a system of integro-differential equations for the PDFs p + ( x , t ) , p ( x , t ) and p 0 ( x , t ) in terms of the switching terms i + ( x , t ) , i ( x , t ) and i 0 ( x , t ) :
p + t + ν p + x = i + ( x , t ) + α + i 0 ( x , t ) ,
p t ν p x = i ( x , t ) + α i 0 ( x , t ) ,
p 0 t = f ( p ) p 0 + i + ( x , t ) + i ( x , t ) ( 1 α 0 ) i 0 ( x , t ) .
Note that the above Equations (30)–(32) can be obtained by using the Fourier–Laplace transform technique. The switching terms i + ( x , t ) and i ( x , t ) describe the average flux of cells from migrating states to proliferating state and i 0 ( x , t ) describes the average flux of cells from the proliferating state to migrating states.
The terms accounting for changes in state i + ( x , t ) , i ( x , t ) and i 0 ( x , t ) are expressed in terms of p + ( x , t ) , p ( x , t ) and p 0 ( x , t ) , respectively, as follows (see Appendix A for the details of the derivation):
i + ( x , t ) = 0 t H + ( t τ ) p + ( x ν ( t τ ) , τ ) d τ ,
i ( x , t ) = 0 t H ( t τ ) p ( x + ν ( t τ ) , τ ) d τ ,
i 0 ( x , t ) = 0 t H 0 ( t τ ) p 0 ( x , τ ) e τ t f ( p ( x , u ) ) d u d τ ,
where H j ( t ) is the memory kernel defined by its Laplace transform [37]:
H ˜ j ( s ) = ψ ˜ j ( s ) Ψ ˜ j ( s ) , j = + , , 0 .
Here, ψ ˜ j ( s ) and Ψ ˜ j ( s ) are the Laplace transforms of the residence time PDF ψ j ( τ ) and the survival function Ψ j ( τ ) , respectively. The nonlinear master Equations (30)–(32) together with interaction terms (33)–(35) are a generalization of a linear system of equations obtained in [33,38].
The Fourier–Laplace transform of the total density p ( x , t ) = p + ( x , t ) + p ( x , t ) + p 0 ( x , t ) , where the proliferating rate f ( n ) = 0 , can be written as (see Appendix B)
p ˜ ^ ( k , s ) = p ^ + 0 α + ψ ˜ + Ψ ˜ 0 + α + Ψ ˜ + α α + ψ ˜ 0 ψ ˜ Ψ ˜ + + α α + ψ ˜ 0 ψ ˜ + Ψ ˜ α 0 α + ψ ˜ 0 Ψ ˜ + + p ^ 0 α ψ ˜ Ψ ˜ 0 + α Ψ ˜ α α + ψ ˜ 0 ψ ˜ + Ψ ˜ + α α + ψ ˜ 0 ψ ˜ Ψ ˜ + α 0 α ψ ˜ 0 Ψ ˜ + p ^ 0 0 Ψ ˜ 0 + α + ψ ˜ 0 Ψ ˜ + + α ψ ˜ 0 Ψ ˜ 1 α + ψ ˜ 0 ψ ˜ + α ψ ˜ 0 ψ ˜ α 0 ψ ˜ 0 ,
where p ^ j 0 = p j 0 ( x ) e i k x d x , j = + , , 0 .
In what follows, we will examine various distributions for residence time PDFs, ψ j ( τ ) , including power law distributions.

4. Examples of Three-State Switching Model

In this section, two examples are considered, one of which is Markovian, while the other is non-Markovian.

4.1. Markovian Case

Let us consider the case when the switching rates β j ( τ ) are constant. This case corresponds to the exponential time PDFs ψ j ( τ ) = β j e β j τ , and the exponential survival functions Ψ j ( τ ) = e β j τ , for which H ˜ j ( s ) = β j and H ( τ ) = β j δ ( τ ) , j = + , , 0 . It follows from Equations (16) and (17) that the switching terms of cells are
i + ( x , t ) = β + p + ( x , t ) , i ( x , t ) = β p ( , t ) , i 0 ( x , t ) = β 0 p 0 ( x , t ) .
For example, if we consider only the simple case of a switching process for which α + = 1 / 2 , α = 1 / 2 and α 0 = 0 and assume β + ( τ ) = β ( τ ) = β and β 0 ( τ ) = β 0 , then Equations (30)–(32) can be reduced to a classical three-state Markovian model for the mean densities of cells of migratory phenotype p + ( x , t ) and p ( x , t ) and the mean density of cells of proliferating phenotype p 0 ( x , t ) :
p + t + ν p + x = β p + + 1 2 β 0 p 0 ,
p t ν p x = β p + 1 2 β 0 p 0 ,
p 0 t = f ( p ) p 0 + β p + + β p β 0 p 0 .
If there is no proliferating in state 0 (in the case when f ( p ) = 0 ), then Equations (39)–(41) become the same as the governing equations in [38], where self-reinforcement is taken into consideration and the single governing hyperbolic partial differential equation for the probability density of random walk position is derived, from which the second moment is obtained in the long-time limit.

4.2. Non-Markovian Case

Now, we consider another distribution of the residence times ψ j ( t ) that corresponds to the family of gamma distributions with parameters γ and β j , that is
ψ j ( τ ) = β j γ τ γ 1 e β j τ Γ ( γ ) , j = + , , 0 .
Notice immediately that when γ = 1 , then (42) reduces to the form of the exponential distribution PDF, covered in the previous subsection. Hence, exponential density is a special case of gamma density. For example, if the parameter γ = 2 , then the residence time PDF and the survival probability turn out to be
ψ j ( τ ) = β j 2 τ e β j τ , Ψ j ( τ ) = ( β j τ + 1 ) τ e β j τ , j = + , , 0 ,
and their corresponding Laplace transforms are
ψ ˜ j ( s ) = β j 2 β j + s 2 , Ψ ˜ j ( s ) = 2 β j + s β j + s 2 , j = + , , 0 .
In this case, the Laplace transform of the memory kernel is
H ˜ j ( s ) = β j 2 2 β j + s , j = + , , 0 .
Taking the inverse Laplace transform, the memory kernel has an exponential form
H j ( τ ) = β j 2 e 2 β j τ , j = + , , 0 .
Therefore, the switching terms (33)–(35) are now written as
i + ( x , t ) = β + 2 0 t e 2 β + τ p + ( x ν ( t τ ) , τ ) d τ ,
i ( x , t ) = β 2 0 t e 2 β τ p ( x + ν ( t τ ) , τ ) d τ ,
i 0 ( x , t ) = β 0 2 0 t e 2 β 0 τ p 0 ( x , τ ) e τ t f ( p ( x , u ) ) d u d τ .
In this case, it is evident that it is no longer Markovian. The defining characteristic of the flux from the proliferating state to the migratory states i 0 ( x , t ) in this non-Markovian case is its dependence on proliferating rate f ( p ) , which is reflected in its exponential term e τ t f ( p ( x , u ) ) d u , while this phenomenon does not exist in the Markovian case. Similar outcomes were observed in the transport and nonlinear reaction system in [39]. In this paper, the main attention is directed toward the anomalous case, in which the velocities of cell migration alternate at random times distributed in accordance with an exponential distribution, and the behavior of proliferation is characterized by rest times distributed following a power law distribution.

5. Anomalous Three-State Model with Mittag–Leffler Distributed Rest Times

In this section, we consider an anomalous switching case and determine an average position of a cancer cell, x ( t ) , and the mean squared displacement (MSD), x 2 ( t ) , when there is no proliferating process in state 0: f ( p ) = 0 . Given that the proliferation of cancer cells is influenced by numerous conditions, it is assumed that a characteristic scale for the proliferating residence time is not present [6]. As a result, the PDF for the residence time in the proliferating state exhibits a power law behavior. Thus, we introduce a switching term, i 0 ( x , t ) , that incorporates the Riemann–Liouville fractional derivative, which is equivalent to the waiting time for a cell in the rest state, with a corresponding PDF derived from a Mittag–Leffler function. Its survival function is expressed as
Ψ 0 ( τ ) = E β 0 τ τ 0 β 0 , 0 < β 0 < 1 ,
then
ψ 0 ( τ ) = d d τ Ψ 0 ( τ ) = d d t E β 0 τ τ 0 β 0 ,
where τ 0 is a parameter with units of time and E β 0 ( . ) is the one-parameter Mittag–Leffler function. The parameter β 0 is the measure of the strength of the resting state [6]. A decrease in β 0 increases the probability of a longer residence time in the proliferating state. As β 0 approaches 1, the Markovian case with exponentially distributed resting times, τ 0 1 e t / τ 0 , is recovered.
The residence times in the proliferating state will be distributed approximately as
ψ 0 ( τ ) τ 0 τ 1 + β 0 , 0 < β 0 < 1 ,
for large values of τ / τ 0 (as τ ), leading to the formation of “heavy” or power law tail distributions. Power law distributions have been widely observed in various empirical studies of stochastic processes with intricate underlying mechanisms. Additionally, numerous examples have demonstrated the significance of power law waiting times in real-world phenomena.
The Laplace transforms of residence time PDF ψ 0 ( τ ) and its corresponding survival function Ψ 0 ( τ ) are
ψ ˜ 0 ( s ) = 1 1 + ( s τ 0 ) β 0 , Ψ ˜ 0 ( s ) = τ 0 ( s τ 0 ) β 0 1 1 + ( s τ 0 ) β 0 .
The Laplace transform ψ ˜ 0 ( s ) corresponding to Equation (52) can be approximated by
ψ ˜ 0 ( s ) 1 τ 0 s β 0 , 0 < β 0 < 1 ,
for small s[40]. The mean waiting time τ = 0 + τ ψ 0 ( τ ) d τ is infinite in this case. The memory kernel H 0 ( t ) in terms of its Laplace transform is
H ˜ 0 ( s ) s ( τ 0 s ) β 0 .
In this situation, it is more convenient to find the switching term i 0 ( x , t ) using the Laplace transform version of Equation (35):
i ˜ 0 ( x , s ) = H ˜ 0 ( s ) p ˜ 0 ( x , s ) .
Inserting H ˜ 0 ( s ) into the above equation and then taking the inverse Laplace transform gives
i 0 ( x , t ) = 1 τ 0 β 0 D t 1 β 0 p 0 ( x , t ) ,
where D t 1 β is the Riemann–Liouville fractional derivative defined in (5). Contrary to the proliferating process, an average transport time is finite. Therefore, we assume that the residence time PDFs for migrating states are exponentially distributed with rate β
ψ + ( τ ) = ψ ( τ ) = β e β τ ,
where β is constant. The Laplace transform of ψ + ( τ ) and ψ ( τ ) are
ψ ˜ + ( s ) = ψ ˜ ( s ) = β β + s ,
and the Laplace transform of their corresponding survival functions Ψ + ( τ ) and Ψ ( τ ) are
Ψ ˜ + ( s ) = Ψ ˜ ( s ) = 1 β + s .
In order to comprehend the implications of incorporating a heavy-tailed waiting time distribution for the rest state, two basic scenarios are analyzed: a single velocity model and a symmetric dual velocity model, both of which possess a non-Markovian rest state. The first and second moments at the long-term limit are performed for each of these models.

5.1. Single Migration State Model

The aim of this subsection is to show that if the migrating cancer cells move with a constant velocity ν , then the mean cell position x ( t ) increases as t β 0 for 0 < β 0 < 1 . Consider the particular case when all cells at t = 0 start to move to the positive direction with the velocity ν from the point x = 0 . To achieve this, we set α + = 1 , α = 0 and α 0 = 0 , for which the initial density of cells in migrating states are p + 0 ( x ) = δ ( τ ) and p 0 ( x ) = 0 , and the initial density of cells in proliferating state, p 0 0 ( x ) , is zero. Then, p ^ + 0 ( k ) = 1 , p ^ 0 ( k ) = 0 and p ^ 0 0 ( k ) = 0 . The system of integro-differential equations for the PDFs p + ( x , t ) and p 0 ( x , t ) in terms of the switching term i 0 ( x , t ) is
p + t + ν p + x = β p + ( x , t ) + 1 τ 0 β 0 D t 1 β 0 p 0 ( x , t ) ,
p 0 t = β p + ( x , t ) 1 τ 0 β 0 D t 1 β 0 p 0 ( x , t ) .
It follows from Equation (37) that
p ˜ ^ ( k , s ) = p ^ + 0 α + ψ ˜ + Ψ ˜ 0 + α + Ψ ˜ + 1 α + ψ ˜ 0 ψ ˜ + .
In the limit s 0 , we find from (54), (59) and (60) that
p ˜ ^ ( k , s ) β τ 0 β 0 s 1 β 0 ( i k v + β τ 0 β 0 s β 0 ) .
To find the average position of cancer cells for large time asymptotic t , we use the Laplace transform of x ( t ) defined by
x ( s ) = i k p ˜ ^ ( k , s ) | k = 0 ,
where p ˜ ^ ( k , s ) = p ˜ ^ + ( k , s ) + p ˜ ^ ( k , s ) + p ˜ ^ 0 ( k , s ) is the Fourier–Laplace transform of the total density of cells p = p + + p + p 0 .
By using Equation (65), we find
x ( s ) ν s 1 + β 0 β τ 0 β 0 .
Finally, taking the inverse Laplace transform, the average position of cells is
x ( t ) ν Γ ( 1 + β 0 ) β τ 0 β 0 t β 0 , 0 < β 0 < 1 ,
which is sublinear. The average position of cancer cells was estimated in [6] using a different idea, yet the same anomalous behavior was obtained. By employing a similar approach, the MSD, x 2 ( t ) , of this model can be determined and show that x 2 ( t ) t 2 β 0 . The anomalous advection reflects the memory effect associated with the slow motion of cells due to the power law of the residence time distribution for the proliferating state with an infinite mean residence time (52).

5.2. Symmetric Active States Model

The objective of this subsection is to demonstrate by calculation the MSD that, upon expanding the previously considered model by incorporating an additional active state with a velocity of ν , results in subdiffusion in the long-time asymptotic limit.
By introducing the flux
J ( x , t ) = ν p + ( x , t ) ν p ( x , t ) ,
Equations (30)–(32) can be combined and reduced to a system of governing equations for the total density function p ( x , t ) = p + ( x , t ) + p ( x , t ) + p 0 ( x , t ) and p 0 ( x , t ) .
By adding (30)–(32), we obtain
p t = J x .
Now multiplying (30) and (31) by ν and subtracting the resulting equations, and by setting α + = 1 / 2 , α = 1 / 2 and α 0 = 0 , we obtain
J t = ν 2 p x + ν 2 p 0 x ν i + ( x , t ) i ( x , t ) .
Then, by differentiating (69) with respect to t and (70) with respect to x and eliminating the flux 2 J / t x , and combining it with (32), we arrive at the governing equations
2 p t 2 = ν 2 2 p x 2 ν 2 2 p 0 x 2 ν x i + ( x , t ) i ( x , t ) ,
p 0 t = i + ( x , t ) + i ( x , t ) i 0 ( x , t ) .
We consider now the symmetrical initial conditions for which the cells start to move from the point x = 0 at t = 0 as follows: p + 0 ( x ) = δ ( τ ) / 2 and p 0 ( x ) = δ ( τ ) / 2 , and the initial density of cells in proliferating state, p 0 0 ( x ) , is zero. Then p ^ + 0 ( k ) = 1 / 2 , p ^ 0 ( k ) = 1 / 2 and p ^ 0 0 ( k ) = 0 . From (71) and (72), and by using β ( p + + p ) = β p β p 0 , we obtain a reduced system of integro-differential equations for the PDFs p ( x , t ) and p 0 ( x , t ) :
2 p t 2 = ν 2 2 p x 2 ν 2 2 p 0 x 2 β p t ,
p 0 t = β p ( x , t ) β p 0 ( x , t ) 1 τ 0 β 0 D t 1 β 0 p 0 ( x , t ) .
Performing the Fourier–Laplace transform to (73) and (74), we obtain
p ˜ ^ ( k , s ) = s + β s 2 + s β + ν 2 k 2 β ν 2 k 2 s + β + τ 0 β 0 s 1 β 0 .
In the limit s 0 ,
p ˜ ^ ( k , s ) 1 s + ν 2 β 2 τ 0 β 0 k 2 s 1 β 0 .
Mathematically, the MSD is defined by
x 2 ( t ) = x 2 ( t ) p ( x , t ) d x ,
which in the Fourier–Laplace space reads
x 2 ( s ) = 2 k 2 p ˜ ^ ( k , s ) | k = 0 .
Calculating the MSD at the long time limit using the above formula, we find
x 2 ( s ) 2 ν 2 s 1 + β 0 β 2 τ 0 β 0 ,
and subsequent Laplace inversion, we obtain
x 2 ( t ) 2 ν 2 ( 1 + β 0 ) β 2 τ 0 β 0 t β 0 , 0 < β 0 < 1 .
Once more, it is observed that the power law PDF for the residence time in the proliferating state results in subdiffusive behavior. The calculation of the second moment is performed utilizing the Laplace technique in [33]; however, the same anomalous behavior is observed.

6. Conclusions

We developed a three-state non-Markovian model to characterize the migration–proliferation dichotomy of tumor cells with general residence time distributions. Specifically, the model differentiates between a migratory state, where cells move without undergoing proliferation, and a proliferative state, where they do not migrate but instead undergo division. The objective of this study was to provide a mesoscopic description of the anomalous transport and reactions associated with the migration–proliferation dichotomy mechanism. This work further analyzed the two-state non-Markovian model for the migration–proliferation dichotomy of cancer cells introduced in [6] by incorporating an additional state for migrating cells. It also extended the stochastic transport of particles outlined in [33]. In comparison to the work of Han et al. (2021), we enhanced the model by introducing non-Markovian running states. In particular, we assumed that the running time follows an arbitrary distribution instead of being limited to an exponential distribution as stated in [33]. Additionally, we incorporated a reaction term into the non-Markovian model as another extension of [33]. The developed probabilistic approach for the transport of tumor cells with directional in migration speed and a rest state was described using a persistent random walk with an arbitrary residence time distribution, while the proliferation rate was modeled using a nonlinear function of the densities of all three cell types. A set of balance equations for cancer cells with two phenotypes, which can switch randomly between cell proliferation and migration, was derived with a non-local switching term involving the Riemann–Liouville fractional derivative. To achieve this, the Markovian model was employed as the initial approach, under the premise that the transition probabilities depend on the residence time variable.
To investigate the nature of our model in the anomalous case involving power law residence time distributions, we calculated the average position of cancer cells in a single migration state model that corresponds to sublinear growth in time: x ( t ) t β 0 , for 0 < β 0 < 1 . The MSD for symmetric bi-directional migration model was estimated and showed that x 2 ( t ) t β 0 , for 0 < β 0 < 1 . We demonstrated through analytical methods that if cells move in a migratory state with constant velocity ν , using power law residence time distributions for the proliferative state (the Mittag–Leffler distributed waiting times for rests) results in subdiffusive behavior in the long time asymptotic limit. This describes the impact of an anomalous rest state on persistent random walks of cells with finite velocity: the longer a cell survives in a proliferative state, the smaller the probability of switching to a migratory state.
Our proposed model offers a general framework for analyzing reaction–transport systems characterized by non-Markovian and anomalous transitions between active and inactive states. This model is suitable for bi-directional intracellular transport [41] that involves a resting state for power law distributed times. Recently, it was experimentally found that the residence time of active movement follows a power law distribution [42]. The presence of subdiffusive behavior in a stochastic transport system due to anomalous switching constitutes a beneficial aspect of the proposed model, enhancing its potential applications, particularly in describing biological or ecological movement, including the fields of population theory and cellular biology. Our intention is to apply the model to intracellular transport.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The author expresses his gratitude to Sergei Fedotov for the valuable discussions.

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CTRWContinuous Time Random Walk
PDFProbability Density Function
MSDMean Squared Displacement

Appendix A. Derivation of The Switching Terms

In this appendix, we demonstrate the process of deriving integral equations for i j ( x , t ) and the unstructured density p j ( x , t ) , and determining the structure of the switching terms i j ( x , t ) in terms of p j ( x , t ) . The derivation of (33)–(35) involves convolutions in Fourier–Laplace space via the transformation
F L { i j ( x , t ) } = i ˜ ^ j ( k , s ) = 0 i j ( x , t ) e s t + i k x d t d x j = + , , 0 ,
F L { p j ( x , t ) } = p ˜ ^ j ( k , s ) = 0 p j ( x , t ) e s t + i k x d t d x j = + , , 0 .
The substitution of (23)–(25) into (17) and (16), along with the initial condition (12), gives
i + ( x , t ) = α + 0 t ψ + ( τ ) i 0 ( x ν τ , t τ ) d τ + α + p + 0 ( x ν t ) ψ + ( t ) ,
i ( x , t ) = α 0 t ψ ( τ ) i 0 ( x + ν τ , t τ ) d τ + α p 0 ( x + ν t ) ψ ( t ) ,
i 0 ( x , t ) = 0 t ψ 0 ( τ ) i + ( x , t τ ) + i ( x , t τ ) + α 0 i 0 ( x , t τ ) e t τ t f ( p ( x , u ) ) d u d τ + p 0 0 ( x ) ψ 0 ( t ) e 0 t f ( p ( x , u ) ) d u ,
p + ( x , t ) = α + 0 t Ψ + ( τ ) i 0 ( x ν τ , t τ ) d τ + α + p + 0 ( x ν t ) Ψ + ( t ) ,
p ( x , t ) = α 0 t Ψ ( τ ) i 0 ( x + ν τ , t τ ) d τ + α p 0 ( x + ν t ) Ψ ( t ) ,
p 0 ( x , t ) = 0 t Ψ 0 ( τ ) i + ( x , t τ ) + i ( x , t τ ) + α 0 i 0 ( x , t τ ) e t τ t f ( p ( x , u ) ) d u d τ + p 0 0 ( x ) Ψ 0 ( t ) e 0 t f ( p ( x , u ) ) d u ,
where we use formula (26) to eliminate β j ( τ ) . Taking the Fourier–Laplace transform together with the shift theorem of Equations (A3), (A4), (A6) and (A7), we obtain
i ˜ ^ + ( k , s ) = i ˜ ^ 0 ( k , s ) + p ^ + 0 ( k ) α + ψ ˜ + ( s i k ν ) ,
i ˜ ^ ( k , s ) = i ˜ ^ 0 ( k , s ) + p ^ 0 ( k ) α ψ ˜ ( s + i k ν ) ,
p ˜ ^ + ( k , s ) = i ˜ ^ 0 ( k , s ) + p ^ + 0 ( k ) α + Ψ ˜ + ( s i k ν ) ,
p ˜ ^ ( k , s ) = i ˜ ^ 0 ( k , s ) + p ^ 0 ( k ) α Ψ ˜ ( s + i k ν ) .
Then, from (A9)–(A12), we obtain
i ˜ ^ + ( k , s ) = ψ ˜ + ( s i k ν ) Ψ ˜ + ( s i k ν ) p ˜ ^ + ( k , s ) ,
i ˜ ^ ( k , s ) = ψ ˜ ( s + i k ν ) Ψ ˜ ( s + i k ν ) p ˜ ^ ( k , s ) .
Inverse Fourier–Laplace transform of Equations (A13) and (A14) gives (33) and (34).
The expression for i 0 ( x , t ) in terms of p 0 ( x , t ) can be obtained by multiplying Equations (A5) and (A8) by e 0 t f ( p ( x , u ) ) d u , and taking the Laplace transform L { f } , we obtain
L { i 0 ( x , t ) e 0 t f ( p ( x , u ) ) d u } = p 0 0 ( x ) + L { i + ( x , t ) + i ( x , t ) + α 0 i 0 ( x , t ) e 0 t f ( p ( x , u ) ) d u } ψ ˜ 0 ( s ) ,
L { p 0 ( x , t ) e 0 t f ( p ( x , u ) ) d u } = p 0 0 ( x ) + L { i + ( x , t ) + i ( x , t ) + α 0 i 0 ( x , t ) e 0 t f ( p ( x , u ) ) d u } Ψ ˜ 0 ( s ) .
Then from (A15) and (A16), we obtain
L { i 0 ( x , t ) e 0 t f ( p ( x , u ) ) d u } = L { p 0 ( x , t ) e 0 t f ( p ( x , u ) ) d u } ψ ˜ 0 ( s ) Ψ ˜ 0 ( s ) .
Inverse Laplace transform gives
i 0 ( x , t ) e 0 t f ( p ( x , u ) ) d u = 0 t H 0 ( t τ ) p 0 ( x , τ ) e 0 τ f ( p ( x , u ) ) d u d τ .
From Equation (A18), we obtain (35).

Appendix B. Derivation of The Fourier–Laplace Transform of The Total Density

The purpose of this appendix is to determine the Fourier–Laplace transform of the total density p ( x , t ) for the case when the proliferating rate f ( p ) = 0 . The total density is given by the sum p ( x , t ) = p + ( x , t ) + p ( x , t ) + p 0 ( x , t ) .
If f ( p ) = 0 , then we obtain from (A15) and (A16)
i ˜ ^ 0 ( k , s ) = i ˜ ^ + ( k , s ) + i ˜ ^ ( k , s ) + α 0 i ˜ ^ 0 ( k , s ) + p ^ 0 0 ( k ) ψ ˜ 0 ( s ) ,
p ˜ ^ 0 ( k , s ) = p ˜ ^ + ( k , s ) + i ˜ ^ ( k , s ) + α 0 i ˜ ^ 0 ( k , s ) + p ^ 0 0 ( k ) Ψ ˜ 0 ( s ) .
Then from (A19) and (A20), we find
i ˜ ^ 0 ( k , s ) = ψ ˜ 0 ( s ) Ψ ˜ 0 ( s ) p ˜ ^ 0 ( k , s ) .
Substitution of (A13), (A14) and (A21) into (A11), (A12) and (A20) gives
p ˜ ^ + ( k , s ) = ψ ˜ 0 ( s ) Ψ ˜ 0 ( s ) p ˜ ^ 0 ( k , s ) + p ^ + 0 ( k ) α + Ψ ˜ + ( s i k ν ) ,
p ˜ ^ ( k , s ) = ψ ˜ 0 ( s ) Ψ ˜ 0 ( s ) p ˜ ^ 0 ( k , s ) + p ^ 0 ( k ) α Ψ ˜ ( s + i k ν ) ,
p ˜ ^ 0 ( k , s ) = ψ ˜ + ( s i k ν ) Ψ ˜ + ( s i k ν ) p ˜ ^ + ( k , s ) + ψ ˜ ( s + i k ν ) Ψ ˜ ( s + i k ν ) p ˜ ^ ( k , s ) + α 0 ψ ˜ 0 ( s ) Ψ ˜ 0 ( s ) p ˜ ^ 0 ( k , s ) + p ^ 0 0 ( k ) Ψ ˜ 0 ( s ) .
These three equations can be rewritten as
p ˜ ^ + ( k , s ) = α + ψ ˜ 0 ( s ) ψ ˜ + ( s i k ν ) p ˜ ^ + ( k , s ) + α + ψ ˜ 0 ( s ) Ψ ˜ + ( s i k ν ) ψ ˜ ( s + i k ν ) Ψ ˜ ( s + i k ν ) p ˜ ^ ( k , s ) + α + ψ ˜ 0 ( s ) Ψ ˜ + ( s i k ν ) p ^ 0 0 ( k ) + α + Ψ ˜ + ( s i k ν ) p ^ + 0 ( k ) + α + α 0 ψ ˜ 0 ( s ) Ψ ˜ + ( s i k ν ) ψ ˜ 0 ( s ) Ψ ˜ 0 ( s ) p ˜ ^ 0 ( k , s ) ,
p ˜ ^ ( k , s ) = α ψ ˜ 0 ( s ) ψ ˜ ( s + i k ν ) p ˜ ^ ( k , s ) + α ψ ˜ 0 ( s ) Ψ ˜ ( s + i k ν ) ψ ˜ + ( s i k ν ) Ψ ˜ + ( s i k ν ) p ˜ ^ + ( k , s ) + α ψ ˜ 0 ( s ) Ψ ˜ ( s + i k ν ) p ^ 0 0 ( k ) + α Ψ ˜ ( s + i k ν ) p ^ 0 ( k ) + α α 0 ψ ˜ 0 ( s ) Ψ ˜ ( s + i k ν ) ψ ˜ 0 ( s ) Ψ ˜ 0 ( s ) p ˜ ^ 0 ( k , s ) ,
p ˜ ^ 0 ( k , s ) = α + ψ ˜ 0 ( s ) ψ ˜ + ( s i k ν ) p ˜ ^ 0 ( k , s ) + α ψ ˜ 0 ( s ) ψ ˜ ( s + i k ν ) p ˜ ^ 0 ( k , s ) + α + ψ ˜ + ( s i k ν ) Ψ ˜ 0 ( s ) p ^ + 0 ( k ) + α ψ ˜ ( s + i k ν ) Ψ ˜ 0 ( s ) p ^ 0 ( k ) + Ψ ˜ 0 ( s ) p ^ 0 0 ( k ) + α 0 ψ ˜ 0 ( s ) p ^ 0 ( k ) .
From (A22) and (A23), we obtain
p ˜ ^ + ( k , s ) Ψ ˜ + ( s i k ν ) = α + p ˜ ^ ( k , s ) α Ψ ˜ ( s + i k ν ) α + p ^ 0 ( k ) + α + p ^ + 0 ( k ) ,
p ˜ ^ ( k , s ) Ψ ˜ ( s + i k ν ) = α p ˜ ^ + ( k , s ) α + Ψ ˜ + ( s i k ν ) α p ^ + 0 ( k ) + α p ^ 0 ( k ) ,
ψ ˜ 0 ( s ) p ˜ ^ 0 ( k , s ) Ψ ˜ 0 ( s ) = p ˜ ^ + ( k , s ) α + Ψ ˜ + ( s i k ν ) p ^ + 0 ( k ) ,
ψ ˜ 0 ( s ) p ˜ ^ 0 ( k , s ) Ψ ˜ 0 ( s ) = p ˜ ^ ( k , s ) α Ψ ˜ ( s + i k ν ) p ^ 0 ( k ) .
Substitution of (A29) and (A30) into (A25) gives
p ˜ ^ + ( k , s ) = α + ψ ˜ 0 ( s ) ψ ˜ + ( s i k ν ) p ˜ ^ + ( k , s ) + α ψ ˜ 0 ( s ) ψ ˜ ( s + i k ν ) p ˜ ^ + ( k , s ) α α + ψ ˜ 0 ( s ) Ψ ˜ + ( s i k ν ) ψ ˜ ( s + i k ν ) p ^ + 0 ( k ) α 0 α + ψ ˜ 0 ( s ) Ψ ˜ + ( s i k ν ) p ^ + 0 ( k ) + α + ψ ˜ 0 ( s ) Ψ ˜ + ( s i k ν ) p ^ 0 0 ( k ) + α + Ψ ˜ + ( s i k ν ) p ^ + 0 ( k ) + α 0 ψ ˜ 0 ( s ) p ^ + ( k ) + α α + ψ ˜ 0 ( s ) Ψ ˜ + ( s i k ν ) ψ ˜ ( s + i k ν ) p ^ 0 ( k ) ,
and substitution of (A28) and (A31) into (A26) gives
p ˜ ^ ( k , s ) = α ψ ˜ 0 ( s ) ψ ˜ ( s + i k ν ) p ˜ ^ ( k , s ) + α + ψ ˜ 0 ( s ) ψ ˜ + ( s i k ν ) p ˜ ^ ( k , s ) α α + ψ ˜ 0 ( s ) Ψ ˜ ( s + i k ν ) ψ ˜ + ( s i k ν ) p ^ 0 ( k ) α 0 α 0 ψ ˜ 0 ( s ) Ψ ˜ ( s + i k ν ) p ^ 0 ( k ) + α ψ ˜ 0 ( s ) Ψ ˜ ( s + i k ν ) p ^ 0 0 ( k ) + α Ψ ˜ ( s + i k ν ) p ^ 0 ( k ) + α 0 ψ ˜ 0 ( s ) p ^ ( k ) + α α + ψ ˜ 0 ( s ) Ψ ˜ ( s + i k ν ) ψ ˜ + ( s i k ν ) p ^ + 0 ( k ) .
It is convenient to write the functions without including their arguments. From (A27), (A32) and (A33), we find explicit expressions for p ˜ ^ + ( k , s ) , p ˜ ^ ( k , s ) and p ˜ ^ 0 ( k , s ) :
p ˜ ^ + = α + Ψ ˜ + α α + ψ ˜ 0 ψ ˜ Ψ ˜ + α 0 α + ψ ˜ 0 Ψ ˜ + p ^ + 0 + α α + ψ ˜ 0 ψ ˜ Ψ ˜ + p ^ 0 + α + ψ ˜ 0 Ψ ˜ + p ^ 0 0 1 α + ψ ˜ 0 ψ ˜ + α ψ ˜ 0 ψ ˜ α 0 ψ ˜ 0 ,
p ˜ ^ = α Ψ ˜ α α + ψ ˜ 0 ψ ˜ + Ψ ˜ α 0 α ψ ˜ 0 Ψ ˜ p ^ 0 + α α + ψ ˜ 0 ψ ˜ + Ψ ˜ p ^ + 0 + α ψ ˜ 0 Ψ ˜ p ^ 0 0 1 α + ψ ˜ 0 ψ ˜ + α ψ ˜ 0 ψ ˜ α 0 ψ ˜ 0 ,
p ˜ ^ 0 = α + ψ ˜ + Ψ ˜ 0 p ^ + 0 + α ψ ˜ Ψ ˜ 0 p ^ 0 + Ψ ˜ 0 p ^ 0 0 1 α + ψ ˜ 0 ψ ˜ + α ψ ˜ 0 ψ ˜ α 0 ψ ˜ 0 .
The Fourier–Laplace transform of the total density p ( x , t ) is p ˜ ^ ( k , s ) = p ˜ ^ + ( k , s ) + p ˜ ^ ( k , s ) + p ˜ ^ 0 ( k , s ) . Thus, using (A34)–(A36), we obtain (37).

References

  1. Giese, A.; Loo, M.A.; Tran, N.; Haskett, D.; Coons, S.W.; Berens, M.E. Dichotomy of astrocytoma migration and proliferation. Int. J. Cancer 1996, 67, 275–282. [Google Scholar] [CrossRef]
  2. Giese, A.; Bjerkvig, R.; Berens, M.E.; Westphal, M. Cost of migration: Invasion of malignant gliomas and implications for treatment. J. Clin. Oncol. 2003, 21, 1624–1636. [Google Scholar] [CrossRef] [PubMed]
  3. Iomin, A. Toy model of fractional transport of cancer cells due to self-entrapping. Phys. Rev. E 2006, 73, 061918. [Google Scholar] [CrossRef] [PubMed]
  4. Fedotov, S.; Iomin, A. Migration and Proliferation Dichotomy in Tumor-Cell Invasion. Phys. Rev. Lett. 2007, 98, 118101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Fedotov, S.; Iomin, A. Probabilistic approach to a proliferation and migration dichotomy in tumor cell invasion. Phys. Rev. E 2008, 77, 031911. [Google Scholar] [CrossRef] [Green Version]
  6. Fedotov, S.; Iomin, A.; Ryashko, L. Non-Markovian models for migration-proliferation dichotomy of cancer cells: Anomalous switching and spreading rate. Phys. Rev. E 2011, 84, 061131. [Google Scholar] [CrossRef] [Green Version]
  7. Iomin, A. Continuous time random walk and migration–proliferation dichotomy of brain cancer. Biophys. Rev. Lett. 2015, 10, 37–57. [Google Scholar] [CrossRef]
  8. Athale, C.; Mansury, Y.; Deisboeck, T.S. Simulating the impact of a molecular ‘decision-process’ on cellular phenotype and multicellular patterns in brain tumors. J. Theor. Biol. 2005, 233, 469–481. [Google Scholar] [CrossRef] [Green Version]
  9. Zhang, L.; Wang, Z.; Sagotsky, J.A.; Deisboeck, T.S. Multiscale agent-based cancer modeling. J. Math. Biol. 2009, 58, 545–559. [Google Scholar] [CrossRef]
  10. Zhang, L.; Chen, L.L.; Deisboeck, T.S. Multi-scale, multi-resolution brain cancer modeling. Math. Comp. Simul. 2009, 79, 2021–2035. [Google Scholar] [CrossRef] [Green Version]
  11. Kolobov, A.V.; Gubernov, V.V.; Polezhaev, A.A. Autowaves in the model of infiltrative tumour growth with migration-proliferation dichotomy. Math. Model. Nat. Phenom. 2011, 6, 27–38. [Google Scholar] [CrossRef]
  12. Khain, E.; Sander, L.M.; Stein, A.M. A model for glioma growth. Complexity 2005, 11, 53–57. [Google Scholar] [CrossRef] [Green Version]
  13. Khain, E.; Sander, L.M. Dynamics and pattern formation in invasive tumor growth. Phys. Rev. Lett. 2006, 96, 188103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Chauviere, A.; Preziosi, L.; Byrne, H. A model of cell migration within the extracellular matrix based on a phenotypic switching mechanism. Math. Med. Biol. 2010, 27, 255–281. [Google Scholar] [CrossRef] [PubMed]
  15. Tektonidis, M.; Hatzikirou, H.; Chauvière, A.; Simon, M.; Schaller, K.; Deutsch, A. Identification of intrinsic in vitro cellular mechanisms for glioma invasion. J. Theor. Biol. 2011, 287, 131–147. [Google Scholar] [CrossRef]
  16. Hatzikirou, H.; Basanta, D.; Simon, M. Schaller, K.; Deutsch, A. ‘Go or grow’: The key to the emergence of invasion in tumour progression? Math. Med. Biol. 2012, 29, 49–65. [Google Scholar] [CrossRef]
  17. Conte, M.; Surulescu, C. Mathematical modeling of glioma invasion: Acid-and vasculature mediated go-or-grow dichotomy and the influence of tissue anisotropy. Appl. Math. Comput. 2021, 407, 126305. [Google Scholar] [CrossRef]
  18. Mansury, Y.; Diggory, M.; Deisboeck, T.S. Evolutionary game theory in an agent-based brain tumor model: Exploring the ‘genotype–phenotype’ link. J. Theor. Biol. 2006, 238, 146–156. [Google Scholar] [CrossRef]
  19. Böttger, K.; Hatzikirou, H.; Chauviere, A.; Deutsch, A. Investigation of the migration/proliferation dichotomy and its impact on avascular glioma invasion. Math. Model. Nat. Phenom. 2012, 7, 105–135. [Google Scholar] [CrossRef] [Green Version]
  20. Godlewski, J.; Nowicki, M.O.; Bronisz, A.; Nuovo, G.; Palatini, J.; De Lay, M.; Van Brocklyn, J.; Ostrowski, M.C.; Chiocca, E.A.; Lawler, S.E. MicroRNA-451 regulates LKB1/AMPK signaling and allows adaptation to metabolic stress in glioma cells. Mol. Cell 2010, 37, 620–632. [Google Scholar] [CrossRef] [Green Version]
  21. Engwer, C.; Knappitsch, M.; Surulescu, C. A multiscale model for glioma spread including cell-tissue interactions and proliferation. Math. Biosci. Eng. 2016, 13, 443. [Google Scholar] [CrossRef] [PubMed]
  22. Hunt, A.; Surulescu, C. A multiscale modeling approach to glioma invasion with therapy. Vietnam J. Math. 2017, 45, 221–240. [Google Scholar] [CrossRef] [Green Version]
  23. Kumar, P.; Li, J.; Surulescu, C. Multiscale modeling of glioma pseudopalisades: Contributions from the tumor microenvironment. J. Math. Biol. 2021, 82, 1–45. [Google Scholar] [CrossRef] [PubMed]
  24. Fedotov, S.; Okuda, Y. Non-Markovian random processes and traveling fronts in a reaction-transport system with memory and long-range interactions. Phys. Rev. E 2002, 66, 021113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Huda, S.; Weigelin, B.; Wolf, K.; Tretiakov, K.V.; Polev, K.; Wilk, G.; Iwasa, M.; Emami, F.S.; Narojczyk, J.W.; Banaszak, M.; et al. Lévy-like movement patterns of metastatic cancer cells revealed in microfabricated systems and implicated in vivo. Nat. Commun. 2018, 9, 4539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Fritz, M.; Kuttler, C.; Rajendran, M.L.; Wohlmuth, B.; Scarabosio, L. On a subdiffusive tumour growth model with fractional time derivative. IMA J. Appl. Math. 2021, 86, 688–729. [Google Scholar] [CrossRef]
  27. Othmer, H.G.; Dunbar, S.R.; Alt, W. Models of dispersal in biological systems. J. Math. Biol. 1988, 26, 263–298. [Google Scholar] [CrossRef]
  28. Rossetto, V. The one-dimensional asymmetric persistent random walk. J. Stat. Mech. Theory Exp. 2018, 2018, 043204. [Google Scholar] [CrossRef] [Green Version]
  29. Méndez, V.; Fedotov, S.; Horsthemke, W. Reaction-Transport Systems: Mesoscopic Foundations, Fronts, and Spatial Instabilities; Springer Science & Business Media: Berlin, Germany, 2010. [Google Scholar]
  30. Bouin, E.; Calvez, V.; Nadin, G. Hyperbolic traveling waves driven by growth. Math. Models Methods Appl. Sci. 2014, 24, 1165–1195. [Google Scholar] [CrossRef] [Green Version]
  31. Méndez, V.; Campos, D.; Horsthemke, W. Growth and dispersal with inertia: Hyperbolic reaction-transport systems. Phys. Rev. E 2014, 90, 042114. [Google Scholar] [CrossRef] [Green Version]
  32. Selmeczi, D.; Mosler, S.; Hagedorn, P.H.; Larsen, N.B.; Flyvbjerg, H. Cell motility as persistent random motion: Theories from experiments. Biophys. J. 2005, 89, 912–931. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Han, D.; Alexandrov, D.V.; Gavrilova, A.; Fedotov, S. Anomalous Stochastic Transport of Particles with Self-Reinforcement and Mittag–Leffler Distributed Rest Times. Fractal Fract. 2021, 5, 221. [Google Scholar] [CrossRef]
  34. Fedotov, S. Single integrodifferential wave equation for a Lévy walk. Phys. Rev. E 2016, 93, 020101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Cox, D.R.; Miller, H.D. The Theory of Stochastic Processes; Routledge: Oxfordshire, UK, 2017. [Google Scholar]
  36. Vlad, M.O.; Ross, J. Systematic derivation of reaction-diffusion equations with distributed delays and relations to fractional reaction-diffusion equations and hyperbolic transport equations: Application to the theory of Neolithic transition. Phys. Rev. E 2002, 66, 061908. [Google Scholar] [CrossRef] [PubMed]
  37. Kenkre, V.M.; Montroll, E.W.; Shlesinger, M.F. Generalized master equations for continuous-time random walks. J. Stat. Phys. 1973, 9, 45–50. [Google Scholar] [CrossRef]
  38. Fedotov, S.; Han, D.; Ivanov, A.O.; da Silva, M.A.A. Superdiffusion in self-reinforcing run-and-tumble model with rests. Phys. Rev. E 2022, 105, 014126. [Google Scholar] [CrossRef]
  39. Fedotov, S.; Al-Shamsi, H.; Ivanov, A.; Zubarev, A. Anomalous transport and nonlinear reactions in spiny dendrites. Phys. Rev. E 2010, 82, 041103. [Google Scholar] [CrossRef] [Green Version]
  40. Metzler, R.; Klafter, J. The random walk’s guide to anomalous diffusion: A fractional dynamics approach. Phys. Rep. 2000, 339, 1–77. [Google Scholar] [CrossRef]
  41. Fedotov, S.; Falconer, S. Nonlinear degradation-enhanced transport of morphogens performing subdiffusion. Phys. Rev. E 2014, 89, 012107. [Google Scholar] [CrossRef] [Green Version]
  42. Fedotov, S.; Korabel, N.; Waigh, T.A.; Han, D.; Allan, V.J. Memory effects and Lévy walk dynamics in intracellular transport of cargoes. Phys. Rev. E 2018, 98, 042136. [Google Scholar] [CrossRef] [Green Version]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Al Shamsi, H. Migration and Proliferation Dichotomy: A Persistent Random Walk of Cancer Cells. Fractal Fract. 2023, 7, 318. https://doi.org/10.3390/fractalfract7040318

AMA Style

Al Shamsi H. Migration and Proliferation Dichotomy: A Persistent Random Walk of Cancer Cells. Fractal and Fractional. 2023; 7(4):318. https://doi.org/10.3390/fractalfract7040318

Chicago/Turabian Style

Al Shamsi, Hamed. 2023. "Migration and Proliferation Dichotomy: A Persistent Random Walk of Cancer Cells" Fractal and Fractional 7, no. 4: 318. https://doi.org/10.3390/fractalfract7040318

APA Style

Al Shamsi, H. (2023). Migration and Proliferation Dichotomy: A Persistent Random Walk of Cancer Cells. Fractal and Fractional, 7(4), 318. https://doi.org/10.3390/fractalfract7040318

Article Metrics

Back to TopTop