Next Article in Journal
An Outlook on Hybrid Fractional Modeling of a Heat Controller with Multi-Valued Feedback Control
Previous Article in Journal
An H1-Galerkin Space-Time Mixed Finite Element Method for Semilinear Convection–Diffusion–Reaction Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Non-Markovian Persistent Random Walk Model for Intracellular Transport

1
Department of Mathematics, The University of Manchester, Oxford Road, Manchester M13 9PL, UK
2
Department of Mathematics, College of Science, Sultan Qaboos University, P.O. Box 36, Al Khodh, Muscat 123, Oman
3
Ural Mathematical Center, Ural Federal University, 620000 Ekaterinburg, Russia
*
Authors to whom correspondence should be addressed.
Fractal Fract. 2023, 7(10), 758; https://doi.org/10.3390/fractalfract7100758
Submission received: 8 September 2023 / Revised: 9 October 2023 / Accepted: 12 October 2023 / Published: 15 October 2023
(This article belongs to the Section Life Science, Biophysics)

Abstract

:
Transport of vesicles and organelles inside cells consists of constant-speed bidirectional movement along cytoskeletal filaments interspersed by periods of idling. This transport shows many features of anomalous diffusion. In this paper, we develop a non-Markovian persistent random walk model for intracellular transport that incorporates the removal rate of organelles. The model consists of two active states with different speeds and one resting state. The organelle transitions between states with switching rates that depend on the residence time the organelle spends in each state. The mesoscopic master equations that describe the average densities of intracellular transport in each of the three states are the main results of the paper. We also derive ordinary differential equations for the dynamics for the first and second moments of the organelles’ position along the cell. Furthermore, we analyse models with power-law distributed random times, which reveal the prevalence of the Mittag-Leffler resting state and its contribution to subdiffusive and superdiffusive behaviour. Finally, we demonstrate a non-Markovian non-additivity effect when the switching rates and transport characteristics depend on the rate of organelles removal. The analytical calculations are in good agreement with numerical Monte Carlo simulations. Our results shed light on the dynamics of intracellular transport and emphasise the effects of rest times on the persistence of random walks in complex biological systems.

1. Introduction

Transport processes within eukaryotic cells bear striking similarities to the intricate networks found in large cities. Cytoskeletal networks, analogous to roads and highways, crisscross cells, and facilitate the transport of various cargoes. Intracellular transport is a fundamental biological process involving the movement of molecules, proteins, and organelles from one location to another within cells. It relies on a combination of passive diffusion and ATP-driven movement along cytoskeletal filaments. This process is crucial to maintaining cellular structure, function, and survival, and is involved in essential cellular activities such as growth, division, and signaling. Dysregulation of intracellular transport has been implicated in numerous diseases, including cancer. Therefore, understanding the mechanisms and regulation of intracellular transport is of great importance for both basic cell biology research and clinical applications. In neurons, long-range intracellular transport is even more critical, as its malfunctioning has been linked to adult-onset neurodegenerative diseases, including Alzheimer’s, Parkinson’s, and Huntington’s disease [1].
Vesicles and organelles have been observed to undergo bidirectional motion along cytoskeletal filaments, periodically pausing for finite time intervals. These cargoes utilise molecular motors of opposite polarity in rapid succession [2,3], facilitating their transportation over long distances [4]. The two major families of motors involved in this process are kinesins and dyneins [5,6]. Kinesins transport vesicles along microtubules toward the plus ends, facilitating material transport from the cell interior toward the cortex (anterograde movement). Dyneins move material toward the microtubule minus ends, moving from the cell periphery to the cell interior (retrograde movement). Cytoplasmic dynein predominantly facilitates the minus-end directed long-range transport along microtubules (MTs) [7], and its activity is regulated by multifunctional adaptors such as dynactin [8].
Various modeling approaches have been reviewed to understand intracellular transport [2,9,10], including Markovian models such as reaction–diffusion models [11,12], Brownian ratchet models [13], random walk models [14,15], intermittent search processes [16,17], and exclusion models [18,19,20]. While most of these approaches treat microtubules implicitly, some models explicitly consider them, allowing the investigation of the influence of cytoskeleton topology on intracellular transport [21,22]. A Markovian model of particle diffusion with spatially varying diffusivity shows non-ergodic behaviour [23] frequently observed experimentally in intracellular transport, which was also modelled using non-Markovian fractional Brownian motion [24]. Several multi-state Markovian models exhibit anomalous diffusive regimes before transitioning into normal diffusion at longer times [22,25,26], or even subdiffusion due to topological trapping at longer times [21] or due to non-Markovian run-length-dependent detachment rate of cargo from a microtubule [27]. The effect of advection on diffusive particles was studied using the mobile–immobile model where the particle switch between mobile and immobile states with finite rates [28] and with power-law and mixed trapping time distributions [29]. The non-Markovian mobile-immobile model without advection was used to model the diffusion of excitons in layered perovskites and transition metal dichalcogenides [30].
Given the complexity of intracellular space, cytoskeleton topology, and various interactions within cells, it is not surprising that intracellular transport deviates from standard Brownian diffusion and displays anomalous diffusion [31,32,33,34,35,36,37]. In cellular and molecular biology, anomalous diffusion is also highly heterogeneous [38]. Experimental studies have observed different non-Brownian regimes of diffusion at various time scales [24,27,39,40,41,42,43,44,45,46,47,48,49,50,51,52]. The intracellular environment is highly crowded with macromolecules, subcellular compartments, and confinement domains, leading to subdiffusion at small time scales (less than tens of milliseconds) [53,54,55,56] which transitions to superdiffusion associated with the activity of molecular motors. Whether this superdiffusion persists at longer time scales or if it represents an intermediate regime before transitioning into further normal or subdiffusion regimes is still a subject of debate. Also, the long-time subdiffusion could originate from different mechanisms. At longer times, the influence of cytoskeleton topology [21] and the heterogeneity of the transport process come into play, leading to spurious subdiffusion [24]. In spite of this, the anomalous nature of intracellular transport suggests that non-Markovian models could be viable.
The non-Markovian nature of intracellular transport has been modelled using random walks [27], fractional Brownian motion [24,51,52] and Levy walks [57] and combinations of FBM and CTRW [46] and CTRW and Levy walk [49]. The later stochastic model consists of altering phases of active motion with constant velocity interchanged with periods of passive dynamics. Depending on the distribution of times in the two phases of motion, the model shows transitions from normal diffusive to superdiffusive and to ballistic behavior [58,59,60,61]. Recently, Han et al. [62] (see also [63]) proposed a three-state self-reinforcing random walk model involving Mittag–Leffler distributed rest:
p ± t ± ν p ± x = λ p ± ( x , t ) + r ± τ 0 β D t 1 β p 0 ( x , t ) ,
p 0 t = λ p + ( x , t ) + λ p ( x , t ) ( r + + r ) τ 0 β D t 1 β p 0 ( x , t ) .
Here, p + and p are PDFs of particles moving in the positive and negative direction. The running times are exponentially distributed with the rate λ . The PDF p 0 ( x , t ) describes the rest states drawn from the Mittag-Leffler distribution with parameter β and D t 1 β denotes the Riemann–Liouville fractional derivative given by:
D t 1 β p 0 ( x , t ) = 1 Γ ( β ) t 0 t p 0 ( x , t ) ( t t ) 1 β d t , ( 0 < β < 1 ) .
Three densities satisfy the following normalisation condition:
+ ( p + ( x , t ) + p ( x , t ) + p 0 ( x , t ) ) d x = 1 .
Parameters r + and r describe the probabilities of transitions from the rest. In the self-reinforcing random walk model [62], these probabilities depend on time and space. This three-state self-reinforcing random walk model is characterised by superdiffusive behaviour at intermediate times and by subdiffusion at a longer time scale.
Clearly, in reality, the speeds of cargoes in anterograde and retrograde directions are not equal, since they are powered by different motors. Therefore, in this paper, we analyse intracellular transport using a one-dimensional three-state persistent random walk model with different speeds of cargoes in anterograde and retrograde directions v + and v that generalises the system (1)–(3). We suggest a model in which cargo transitions between moving states and resting states with zero velocity are non-Markovian. We consider the case when the rate of transitions depends on the running time. In order to take into account the random duration of cargo trajectories observed in experiments [24,52], we introduce a removal rate that terminates trajectories at random times. This effect was not considered in a recent paper [64]. The constant removal rate corresponds to the exponentially distributed duration of trajectories. The structure of the paper is as follows. In Section 2, we set up the Three-State Model with a residence time variable and derive the equations for the structured densities in two active states and one passive state. In Section 3, the mesoscopic master equations for densities in all three states are introduced. Specific examples of the Markovian Three-State Model and the Three-State Model with Mittag-Leffler distributed resting times are considered, and equations are derived for the first and second moments. The analytical results of Section 3 are compared in Section 4, in which we perform Monte Carlo simulations. Finally, we discuss our results in Section 5 and give a summary and conclusions in Section 6.

2. Three-State Transport Model with a Residence Time Variable

We consider a model in which cargo randomly switches between three possible states: Two states in which the cargoes move with a constant speed ν + > 0 (state A) or ν > 0 (state R) and a resting state (state 0) in which the cargoes do not move. In this section, we aim to develop a stochastic intracellular transport model to describe the non-Markovian dynamics of intracellular transport.
We introduce the densities that are functions of running and residence time τ [65,66,67,68]. Consider n A ( x , t , τ ) and n R ( x , t , τ ) the densities of cargoes at point x at time t whose residence time in an active state (state A or state R, respectively) lies in the interval ( τ , τ + d τ ) . The density of cargoes in a resting state is n 0 ( x , t , τ ) . The balance equations for structured densities can be written as [66]:
n A t + n A τ + ν + n A x = λ A ( τ ) n A θ A n A ,
n R t + n R τ ν n R x = λ R ( τ ) n R θ R n R ,
n 0 t + n 0 τ = λ 0 ( τ ) n 0 θ 0 n 0 ,
where λ A ( τ ) , λ R ( τ ) and λ 0 ( τ ) are the switching rates that depend on τ and θ A , θ R and θ 0 are the constant removal rates. We assume that the residence time of all cargoes at t = 0 equals zero. In this case, the initial conditions are:
n j ( x , 0 , τ ) = p j 0 ( x ) δ ( τ ) , j = A , R , 0 ,
where p j 0 ( x ) represents the initial densities ( 0 τ t ) .
Three escape (switching) rates, λ A ( τ ) , λ R ( τ ) , λ 0 ( τ ) appearing in (5)–(7) can be defined as [69]:
λ j ( τ ) = ψ j ( τ ) Ψ j ( τ ) , j = A , R , 0 ,
where ψ j ( τ ) are PDFs of the residence times in corresponding states and Ψ j ( τ ) = τ ψ j ( u ) d u are the survival functions. Our approach is very flexible and can be used for various forms of the residence time probability distribution. This flexibility is important, since recently, the non-exponential PDF of the residence time was found in experiments [27].
The boundary conditions at zero running time ( τ = 0 ) can be written as follows:
n A ( x , t , 0 ) = α A 0 t λ A ( τ ) n A ( x , t , τ ) d τ + β A 0 t λ R ( τ ) n R ( x , t , τ ) d τ + γ A 0 t λ 0 ( τ ) n 0 ( x , t , τ ) d τ ,
n R ( x , t , 0 ) = α R 0 t λ A ( τ ) n A ( x , t , τ ) d τ + β R 0 t λ R ( τ ) n R ( x , t , τ ) d τ + γ R 0 t λ 0 ( τ ) n 0 ( x , t , τ ) d τ ,
n 0 ( x , t , 0 ) = α 0 0 t λ A ( τ ) n A ( x , t , τ ) d τ + β 0 0 t λ R ( τ ) n R ( x , t , τ ) d τ + γ 0 0 t λ 0 ( τ ) n 0 ( x , t , τ ) d τ .
These equations describe the switching process to the states A, R, and resting state 0. If the cargo moves in the positive direction, it can switch with rate λ A ( τ ) to the opposite direction with the probability α R , to the resting state with the probability α 0 , or continue in the same direction with the probability α A . The cargo moving in the negative direction can switch with rate λ B ( τ ) to the opposite direction with the probability β A , to the resting state with the probability β 0 , or continue in the same direction with the probability β R . Finally, for the resting cargo, it can switch with rate λ 0 ( τ ) to cargo moving in the positive direction with the probability γ A , to cargo moving in the negative direction with the probability γ R , or remain at rest again with the probability γ 0 . Obviously,
α A + α R + α 0 = 1 , β A + β R + β 0 = 1 , γ A + γ R + γ 0 = 1 .
The product λ j ( τ ) n j ( x , t , τ ) gives the escape rate corresponding to a particular residence time τ . If we denote the escape rates from the states A, R and 0 by i A ( x , t ) , i R ( x , t ) and i 0 ( x , t ) , correspondingly. They can be obtained by integrating λ j n j over variable τ from 0 to t [67]:
i j ( x , t ) = 0 t λ j ( τ ) n j ( x , t , τ ) d τ , j = A , R , 0 .
It follows from Equations (10)–(12) and (14) that
n A ( x , t , 0 ) = α A i A ( x , t ) + β A i R ( x , t ) + γ A i 0 ( x , t ) ,
n R ( x , t , 0 ) = α R i A ( x , t ) + β R i R ( x , t ) + γ R i 0 ( x , t ) ,
n 0 ( x , t , 0 ) = α 0 i A ( x , t ) + β 0 i R ( x , t ) + γ 0 i 0 ( x , t ) .
We solve Equations (5)–(7) using the method of characteristics: For active states ( 0 τ < t ) :
n A ( x , t , τ ) = n A ( x ν + τ , t τ , 0 ) e 0 τ λ A ( u ) d u e θ A τ .
n R ( x , t , τ ) = n R ( x + ν τ , t τ , 0 ) e 0 τ λ R ( u ) d u e θ R τ .
For resting state ( 0 τ < t )
n 0 ( x , t , τ ) = n 0 ( x , t τ , 0 ) e 0 τ λ 0 ( u ) d u e θ 0 τ .
One can see that all solutions (18)–(20) contain an exponential factor e 0 τ λ j ( s ) d s , which is recognised as the survival function Ψ j ( τ ) :
Ψ j ( τ ) = e 0 τ λ j ( u ) , j = A , R , 0 .
Using Equations (10)–(12), the solutions in (18) can be rewritten using the survival function from (21) and the escape rates (switching terms) as follows:
n A ( x , t , τ ) = α A i A ( x ν τ , t τ ) + β A i R ( x ν τ , t τ ) + γ A i 0 ( x ν τ , t τ ) × Ψ A ( τ ) e θ A τ ,
n R ( x , t , τ ) = α R i A ( x + ν τ , t τ ) + β R i R ( x + ν τ , t τ ) + γ R i 0 ( x + ν τ , t τ ) × Ψ R ( τ ) e θ R τ ,
n 0 ( x , t , τ ) = α 0 i A ( x , t τ ) + β 0 i R ( x , t τ ) + γ 0 i 0 ( x , t τ ) Ψ 0 ( τ ) e θ 0 τ .
Notice that the residence time PDF, ψ j ( τ ) , is related to the switching rate, λ j ( τ ) , as follows [69]:
ψ j ( τ ) = λ j ( τ ) e 0 τ λ j ( u ) d u , j = A , R , 0 .

3. Mesoscopic Master Equations

Now, we introduce the mesoscopic densities for all three states: p A ( x , t ) , p R ( x , t ) and p 0 ( x , t ) . The functions p A and p R are the densities for cargoes moving in a positive/negative direction and p 0 is the density of cargoes with a resting state. These densities 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 = A , R , 0 .
By differentiating (26) with respect to time t and using (5)–(7), we obtain the mesoscopic master equations for p A , p R and p 0 :
p A t = n A ( x , t , t ) 0 t λ A ( τ ) n A ( x , t , τ ) d τ ν + 0 t n A x d τ 0 t n A τ d τ θ A 0 t n A ( x , t , τ ) d τ ,
p R t = n R ( x , t , t ) 0 t λ R ( τ ) n R ( x , t , τ ) d τ + ν 0 t n R x d τ 0 t n R τ d τ θ R 0 t n R ( x , t , τ ) d τ ,
p 0 t = n 0 ( x , t , t ) 0 t λ 0 ( τ ) n 0 ( x , t , τ ) d τ 0 t n 0 τ d τ θ 0 0 t n 0 ( x , t , τ ) d τ .
Using (10)–(29), we rewrite the system of equations using the escape rates i A , i R , and i 0 :
p A t + ν + p A x = ( 1 α A ) i A ( x , t ) + β A i R ( x , t ) + γ A i 0 ( x , t ) θ A p A ,
p R t ν p R x = α R i A ( x , t ) ( 1 β R ) i R ( x , t ) + γ R i 0 ( x , t ) θ R p R ,
p 0 t = α 0 i A ( x , t ) + β 0 i R ( x , t ) ( 1 γ 0 ) i 0 ( x , t ) θ 0 p 0 .
The escape rates i A , i R , and i 0 can be expressed in terms of p A , p R and p 0 , respectively, as follows (see Appendix A for the details of the derivation):
i A ( x , t ) = 0 t H A ( t τ ) p A ( x ν + ( t τ ) , τ ) e θ A ( t τ ) d τ ,
i R ( x , t ) = 0 t H R ( t τ ) p R ( x + ν ( t τ ) , τ ) e θ R ( t τ ) d τ ,
i 0 ( x , t ) = 0 t H 0 ( t τ ) p 0 ( x , τ ) e θ 0 ( t τ ) d τ ,
where H i ( t ) is the memory kernel defined by its Laplace transform [70]:
H ˜ j ( s ) = ψ ˜ j ( s ) Ψ ˜ j ( s ) , j = A , R , 0 .
Here, ψ ˜ j ( s ) and Ψ ˜ j ( s ) are the Laplace transforms of the residence time PDF ψ j ( τ ) and the survival function Ψ j ( τ ) , respectively.
It is clear from (33)–(35) that the escape rates i A ( x , t ) , i R ( x , t ) and i 0 ( x , t ) depend on the rates of organelles removal. This is the non-Markovian non-additivity effect that leads to the dependence of the transport characteristics on the kinetics of organelles removal [70,71,72]. Note that the equations above can be applied for the analysis of the phenomenon of the migration–proliferation dichotomy in gliomas [73,74,75], anomalous transport in spiny dendrites [70]. It would be interesting to extend this model for 2 D [76,77,78,79].
Next, we will examine various distributions for residence time PDFs, ψ k ( τ ) , including power-law distributions. We will derive the equations for the moments of the random walk position
μ n ( t ) = x n p ( x , t ) d x ,
where p ( x , t ) = p A ( x , t ) + p R ( x , t ) + p 0 ( x , t ) .

3.1. Markovian Three-State Model

In this subsection, we consider the Markovian case when the switching rates for all three states, λ j ( τ ) , are constant. The residence time’s PDFs are exponential:
ψ j ( τ ) = λ j e λ j τ , j = A , R , 0 ,
with the Laplace transform ψ ˜ j ( s ) = λ j λ j + s . The escape rates i A ( x , t ) , i R ( x , t ) and i 0 ( x , t ) do not involve memory effects and can be written as follows:
i A ( x , t ) = λ A p A ( x , t ) , i R ( x , t ) = λ R p R ( x , t ) , i 0 ( x , t ) = λ 0 p 0 ( x , t ) .
The balance equations for the mean densities p A ( x , t ) , p R ( x , t ) and p 0 ( x , t ) are:
p A t + ν + p A x = ( 1 α A ) λ A p A + β A λ R p R + γ A λ 0 p 0 θ A p A ,
p R t ν p R x = α R λ A p A ( 1 β R ) λ R p R + γ R λ 0 p 0 θ R p R ,
p 0 t = α 0 λ A p A + β 0 λ R p R ( 1 γ 0 ) λ 0 p 0 θ 0 p 0 .
A similar Markovian model describing macroscopic intracellular transport of vesicles and organelles has been introduced in the classical paper [11]. Note that without the resting state, the system of (40), (41) can be reduced to the generalised telegraph equation [66].

3.2. Three-State Model with Mittag–Leffler Distributed Resting Times

Now, let us consider another example involving anomalous behavior. For simplicity, we assume that the removal rate θ 0 for the rest state is zero. We consider the case when the cargoes move along microtubules with running times that follow an exponential distribution with rate λ in the active states and pause for a Mittag–Leffler (ML) distributed residence time in the resting state. In this case, the PDF for the residence time exhibits a power law such that the mean residence time does not exist. The survival function can be written in terms of the one-parameter ML function, E μ ( . ) , as:
Ψ 0 ( τ ) = E μ τ τ 0 μ , 0 < μ < 1 ,
where τ 0 is a time scale. The Laplace transforms of the survival function Ψ 0 ( τ ) and residence time PDF ψ 0 ( τ ) are:
Ψ ˜ 0 ( s ) = τ 0 ( s τ 0 ) μ 1 1 + ( s τ 0 ) μ , ψ ˜ 0 ( s ) = 1 1 + ( s τ 0 ) μ .
Then, from (35), we obtain the rate i 0 as:
i 0 ( x , t ) = 1 τ 0 μ D t 1 μ p 0 ( x , t ) ,
where D t 1 μ is the Riemann–Liouville fractional derivative defined in (3). The system of integro-differential Equations (30)–(32) become:
p A t + ν + p A x = ( 1 α A ) λ p A + β A λ p R + γ A τ 0 μ D t 1 μ p 0 ( x , t ) θ A p A ,
p R t ν p R x = α R λ p A ( 1 β R ) λ p R + γ R τ 0 μ D t 1 μ p 0 ( x , t ) θ R p R ,
p 0 t = α 0 λ p A + β 0 λ p R ( 1 γ 0 ) 1 τ 0 μ D t 1 μ p 0 ( x , t ) .
The system of fractional PDEs (46)–(48) can be also rewritten in terms of Caputo fractional derivatives. In [62], the simplified version of (46)–(48) has been studied, in which the authors derived an equation with a telegraph operator; see Equation (8) in [62].

3.3. Moments Equations

In this subsection, we calculate the first two moments defined in (37). We put θ A = θ R = θ 0 = 0 and introduce the following functions:
m i j ( t ) = x i p J ( x , t ) d x i = 0 , 1 , 2 j = A , R , 0 .
The first moment, m 1 ( t ) , and second moment, m 2 ( t ) can be determined as:
m i ( t ) = m i A ( t ) + m i R ( t ) + m i 0 ( t ) , i = 1 , 2
where m i A ( t ) , m i R ( t ) , m i 0 ( t ) can be found from the system of equations:
d m i A d t = i ν + m i 1 A ( 1 α A ) λ m i A + β A λ m i R + γ A τ 0 μ D t 1 μ m i 0 ( t ) ,
d m i R d t = i ν m i 1 R + α R λ m i A ( 1 β R ) λ m i R + γ R τ 0 μ D t 1 μ m i 0 ( t ) ,
d m i 0 d t = α 0 λ m i A + β 0 λ m i R ( 1 γ 0 ) 1 τ 0 μ D t 1 μ m i 0 ( t ) .
From these equations, one can find the equations for the first and second moments
d m 1 d t = ν + m 0 A ν m 0 R , d m 2 d t = 2 ν + m 1 A 2 ν m 1 R .
Let us find the first moment m 1 . For simplicity, we consider the case when
α 0 = β 0 = 1 , α A = α R = β A = β R = 0 .
The equations for m 0 A , m 0 R and m 00 take the form:
d m 0 A d t = λ m 0 A + γ A τ 0 μ D t 1 μ m 00 ( t ) ,
d m 0 R d t = λ m 0 R + γ R τ 0 μ D t 1 μ m 00 ( t ) ,
d m 00 d t = λ m 0 A + λ m 0 R ( 1 γ 0 ) 1 τ 0 μ D t 1 μ m 00 ( t ) .
We take the Laplace transform of (56)–(58), set p A ( x , 0 ) = δ ( x ) , p R ( x , 0 ) = 0 and p 0 ( x , 0 ) = 0 and obtain, in the long-time limit ( s 0 ), expressions for m ^ 0 A ( s ) = 0 m 0 A ( t ) e s t d t and m ^ 0 R ( s ) = 0 m 0 R ( t ) e s t d t :
m ^ 0 A ( s ) = γ A λ τ 0 μ s 1 μ m ^ 00 ( s ) , m ^ 0 R ( s ) = γ R λ τ 0 μ s 1 μ m ^ 00 ( s ) .
From the normalisation condition (4), we find:
m ^ 0 A ( s ) + m ^ 0 R ( s ) + m ^ 00 ( s ) = s 1 ,
and in the limit s 0 , we obtain:
m ^ 00 ( s ) = 1 s γ A + γ R λ 1 ( τ 0 s ) μ .
Finally, from (59) and (61) we obtain the first moment in the long time limit
m 1 ( t ) = γ A v + γ R v λ Γ ( 1 + μ ) τ 0 μ t μ .
It follows from (62) that the anomalous rest state becomes dominant leading to subdiffusive motion for μ < 1 / 2 .
Now, let us find the second moment for the case when the first moment is zero; that is,
γ A = γ R = γ 2 , ν + = ν + = ν .
The main aim is to show again that the rest state is dominant in the long-time limit that leads to subdiffusion. The equation for the second moment is
d m 2 d t = 2 ν ( m 1 A m 1 R ) ,
where m 1 A , m 1 R obey the equations
d m 1 A d t = ν m 0 A λ m 1 A + γ 2 τ 0 μ D t 1 μ m 10 ( t ) ,
d m 1 R d t = ν m 0 R λ m 1 R + γ 2 τ 0 μ D t 1 μ m 10 ( t ) ,
d m 10 d t = λ m 1 A + λ m 1 R γ τ 0 μ D t 1 μ m 10 ( t ) .
We take the Laplace transform and obtain in the limit s 0
m ^ 2 ( s ) = 2 ν λ s ( m ^ 0 A ( s ) + m ^ 0 R ( s ) ) .
Performing straightforward calculations, we obtain
m ^ 2 ( s ) = 2 ν 2 γ s λ 2 1 ( τ 0 s ) μ .
Taking inverse Laplace transform, we obtain the second moment in the long time limit
m 2 ( t ) = 2 ν 2 γ λ 2 Γ ( 1 + μ ) τ 0 μ t μ .
Clearly, this formula gives us the subdiffusive motion.
In the next sections, we performed the Monte Carlo simulations.

4. Monte Carlo Simulations

Monte Carlo simulations of the random walk follow standard procedure:
(1)
Set initial conditions x 0 = 0 and t 0 = 0 . The initial state was randomly selected corresponding to v 0 = v + , v 0 = v or v 0 = 0 .
(2)
For the Markovian Three-State Model, generate an exponentially distributed random time τ 0 = 1 / λ ln ( 1 p ) where p is a uniformly distributed random number in [ 0 , 1 ) . In the Non-Markovian Three-State Model, resting times were drawn from the ML distribution, and running states were generated using exponentially distributed random times. For the resting state, generate the Mittag-Leffler random number τ using the Matlab mlrnd function (Guido Germano (2023). Mittag-Leffler random number generator (https://www.mathworks.com/matlabcentral/fileexchange/19392-mittag-leffler-random-number-generator (accessed on 1 September 2023)), MATLAB R2020b Central File Exchange).
(3)
Update position and time to x 1 = x 0 + v 0 τ 0 , t 1 = t 0 + τ 0 , respectively. Update state by randomly selecting new velocity v 1 = v + , v 1 = v or v 1 = 0 .
(4)
Repeat steps (2) and (3) until the predefined simulation time T m a x is reached. Thus, a single trajectory x ( t ) is generated.
(5)
Repeat steps (1)–(4) N = 10 4 times to generate an ensemble of N trajectories. Note that trajectories have random durations due to step (2).
(6)
The ensemble of trajectories is then analysed by calculating the distribution of positions at a given T m a x . We estimated these distributions using histograms. To quantify the anomalous diffusion in the Three-State Model, we calculated the moments of these distributions as a function of time. In particular, the first and the second moments were calculated as m 1 ( t ) = i = 1 N x ( i ) ( t ) , m 2 ( t ) = i = 1 N ( x ( i ) ( t ) ) 2 . Here, ( i ) denotes the index of a trajectory.
To achieve the simulation results for the Markovian Three-State Model, we compute the variance m 2 ( t ) m 1 2 ( t ) and distributions of walkers with parameters λ = 1 , v + = 1 and v = 1 . The results are shown in Figure 1. We chose equal probabilities of switching from resting to running state γ A = γ R = γ 0 = 1 / 3 . As expected, the PDFs of the walker’s positions follow the Gaussian distribution. The corresponding second moments show ballistic regimes at short times and asymptotic linear growth.
For the Non-Markovian Three-State model, we also find good agreement with the theoretical predictions. First, we calculated the first moment m 1 ( t ) of random walker position as a function of time with zero removal rates Q A = Q R = Q 0 = 0 , λ A = λ R = 1 / 3 , τ 0 = 1 , λ = 1 and v + = 1 , v = 0.1 . We considered two values of anomalous exponent μ = 0.3 and μ = 0.7 . As shown in Figure 2, the numerics agree with the analytical Equation (62) without fitting.
The ML-distributed resting states make the distributions of the walker’s position non-Gaussian. An example of the distributions of random walkers’ position with the resting states generated from the ML distribution with the scaling parameter μ = 0.7 is shown in Figure 3. At long times, the MSD is growing sub linearly with time, m 2 ( t ) t β with β = μ . At the same time, it was previously found that the two-state model with one running state and the ML resting state has the long-time behaviour of the variance m 2 ( t ) m 1 2 ( t ) t β with β = 2 μ . So, for μ < 0.5 the variance grows subdiffusively with time, and for μ > 0.5 the growth is superdiffusive.
We note that the two-state model is equivalent to the three-state model with the speed of one running state equal to zero. The following question arises: what would be the long-time behaviour of the three-state model with ML-distributed resting times with unequal retrograde and anterograde speeds? We calculated the variance of the position for the three-state model with ML-distributed resting times and showed that the model has subdiffusion or superdiffusion at long times depending on the parameters. At long times, the variance depends on cargo speeds and is well fitted by the power-law function:
m 2 ( t ) m 1 2 ( t ) t β , β = μ 1 + | v + v | v + + v .
For v + = v , we recover the subdiffusive behaviour β = μ < 1 predicted by the Equation (70). The dependence of β / μ 1 on the speed v + (with v kept constant) for the three-state model with ML-distributed resting state is shown in Figure 4a. The asymptotic behaviour of different models given by the dependence of exponents of MSDs or variances at long times, β , is shown in Figure 4a. In contrast to the three-state model with equal retrograde and anterograde cargo speeds v + = v , the three-state model with unequal retrograde and anterograde cargo speeds, v + v can have subdiffusive behaviour and superdiffusive behaviour.

5. Discussion

In cells, vesicles are transported by dyneins from the cell periphery to the cell interior (retrograde transport) and kinesins transport them from the cell interior toward the cortex (anterograde transport). Dyneins and kinesins move at different velocities. Thus, in contrast to previous works [62,64], here, we considered different cargo velocities in retrograde and anterograde directions. As a result, the first moment of the cargo position m 1 ( t ) is non-zero but grows linearly with time for the Markovian model or as a power law in the case of the non-Markovian model. In experiments, this could be used to distinguish between Markovian and non-Markovian dynamics.
The variance of cargo position m 2 ( t ) m 1 2 ( t ) bears additional information about the cargo dynamics. For the non-Markovian model with unequal velocities, the variance can grow subdiffusively and superdiffusively depending on the anomalous exponent μ of the distribution of the duration of resting states. Combining predictions for the first moment, variance, and distribution of cargo’s displacements, it would be possible to verify the non-Markovian three-state model using the ensemble of experimental trajectories.

6. Summary and Conclusions

We have introduced a non-Markovian persistent random walk model for intracellular transport that incorporates the removal rate and finite speed of propagation. This model includes random transitions between two active states (anterograde and retrograde movements) and one resting state, with the probability of switching depending on the amount of time the organelle spends in each state. We have derived new mesoscopic integro-differential equations for densities of organelles inside the cell by using the structural densities obeying Markov differential equations of the first order. We have found a non-Markovian non-additivity effect when the switching rates and transport characteristics depend on the rate of organelles removal. Ordinary differential equations describing the dynamics for the first and second moments of the organelles’ position along the cell have been derived. A model of stochastic transport of organelles which takes into account experimentally observed Mittag-Leffler resting time distribution has been analysed in detail. In particular, it determines the subdiffusive and superdiffusive behaviour of the first moment. Notice that the fractional Riemann–Liouville derivative naturally arises as the result of the Mittag-Leffler distribution of residence times in the resting state. For a general form of resting time distribution, one should not expect the escape rate to be expressed with any kind of fractional derivative. Rather, it will be given by an integral form defined by Equation (34) A numerical Monte Carlo algorithm has been developed to calculate moments and organelle mesoscopic densities in cells, taking into account non-Markovian effects caused by the dependence of transition rates on residence times.
In this paper, we propose the mathematical model of intracellular transport. This model is based on experimental observations (see, for example [2,6]) that intracellular transport demonstrates so-called stop-and-go behaviour and is characterised by long time sub-diffusive behaviour. In this work, we modelled the stop-and-go behaviour as switching from active to resting states and showed that the ML distributed time spent in resting states leads to sub-diffusion at long times. So, experimentally, it would be possible to test the predictions of our model, e.g., whether the resting times are described by the ML distribution. In fact, in our previous work [27,57], we found that the distribution of resting times indeed might be better fitted by a power law (notice that the long-time asymptotic of the ML distribution is t 1 μ ). However, more experimental data are needed to make this comparison more reliable.
In this paper, we used the assumption of constant cargo speeds. In experiments, intracellular organelles might be described based on the distribution of speed. Our constant speed assumption could be viewed as utilising the first moment of this distribution. In the future, it will be interesting to extend our work by taking into account the distribution of cargo speeds [51] and the heterogeneity of intracellular transport [24].

Author Contributions

Conceptualization, N.K., H.A.S. and S.F.; methodology, S.F., H.A.S. and A.O.I.; software, N.K.; validation, N.K., H.A.S. and S.F.; formal analysis, H.A.S. and S.F.; investigation, N.K.; writing—original draft preparation, N.K., H.A.S. and S.F.; writing—review and editing, A.O.I.; visualization, N.K.; supervision, N.K. and S.F.; project administration, N.K.; funding acquisition, N.K., S.F. and A.O.I. All authors have read and agreed to the published version of the manuscript.

Funding

N.K. and S.F. acknowledge financial support from EPSRC Grant No. EP/V008641/1. The research was partly supported by the Ural Mathematical Center, Project No. 075-02-2023-935 (AOI).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data supporting this study are included within the article and/or supporting materials.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Derivation of The Switching Terms

In this Appendix, we derive integral equations for i j ( x , t ) and the unstructured density p j ( x , t ) , and determine the structure of the switching functions 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 = A , R , 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 = A , R , 0 .
Substitution of (22)–(24) into (26) and (14), along with the initial condition (8), gives
i A = 0 t ψ A α A i A ( x ν + τ , t τ ) + β A i R ( x ν + τ , t τ ) + γ A i 0 ( x ν + τ , t τ ) e θ A τ d τ + p A 0 ( x ν + τ ) ψ A ( t ) e θ A τ ,
i R = 0 t ψ R α R i A ( x + ν τ , t τ ) + β R i R ( x + ν τ , t τ ) + γ R i 0 ( x + ν τ , t τ ) e θ R τ d τ + p R 0 ( x + ν τ ) ψ R ( t ) e θ R τ ,
i 0 ( x , t ) = 0 t ψ 0 ( τ ) α 0 i A ( x , t τ ) + β 0 i R ( x , t τ ) + γ 0 i 0 ( x , t τ ) e θ 0 τ d τ + p 0 0 ( x ) ψ 0 ( t ) e θ 0 τ ,
p A = 0 t Ψ A α A i A ( x ν + τ , t τ ) + β A i R ( x ν + τ , t τ ) + γ A i 0 ( x ν + τ , t τ ) e θ A τ d τ + p A 0 ( x ν + τ ) Ψ A ( t ) e θ A τ ,
p R = 0 t Ψ R α R i A ( x + ν τ , t τ ) + β R i R ( x + ν τ , t τ ) + γ R i 0 ( x + ν τ , t τ ) e θ R τ d τ + p R 0 ( x + ν τ ) Ψ R ( t ) e θ R τ ,
p 0 ( x , t ) = 0 t Ψ 0 ( τ ) α 0 i A ( x , t τ ) + β 0 i R ( x , t τ ) + γ 0 i 0 ( x , t τ ) e θ 0 τ d τ + p 0 0 ( x ) Ψ 0 ( t ) e θ 0 τ ,
where we use Formula (25) to eliminate β j ( τ ) . Taking the Fourier-Laplace transform together with shift theorem of Equations (A3)–(A8), we get
i ˜ ^ A ( k , s ) = α A i ˜ ^ A ( k , s ) + β A i ˜ ^ R ( k , s ) + γ A i ˜ ^ 0 ( k , s ) + p ^ A 0 ( k ) ψ ˜ A ( s i k ν + + θ A ) ,
i ˜ ^ R ( k , s ) = α R i ˜ ^ A ( k , s ) + β R i ˜ ^ R ( k , s ) + γ R i ˜ ^ 0 ( k , s ) + p ^ R 0 ( k ) ψ ˜ R ( s + i k ν + θ R ) ,
i ˜ ^ 0 ( k , s ) = α 0 i ˜ ^ A ( k , s ) + β 0 i ˜ ^ R ( k , s ) + γ 0 i ˜ ^ 0 ( k , s ) + p ^ 0 0 ( k ) ψ ˜ 0 ( s + θ 0 ) ,
p ˜ ^ A ( k , s ) = α A i ˜ ^ A ( k , s ) + β A i ˜ ^ R ( k , s ) + γ A i ˜ ^ 0 ( k , s ) + p ^ A 0 ( k ) Ψ ˜ A ( s i k ν + + θ A ) ,
p ˜ ^ R ( k , s ) = α R i ˜ ^ A ( k , s ) + β R i ˜ ^ R ( k , s ) + γ R i ˜ ^ 0 ( k , s ) + p ^ R 0 ( k ) Ψ ˜ R ( s + i k ν + θ R ) ,
p ˜ ^ 0 ( k , s ) = α 0 i ˜ ^ A ( k , s ) + β 0 i ˜ ^ R ( k , s ) + γ 0 i ˜ ^ 0 ( k , s ) + p ^ 0 0 ( k ) Ψ ˜ 0 ( s + θ 0 ) ,
Then from (A9)–(A13), we obtain
i ˜ ^ A ( k , s ) = ψ ˜ A ( s i k ν + + θ A ) Ψ ˜ A ( s i k ν + + θ A ) p ˜ ^ A ( k , s ) ,
i ˜ ^ R ( k , s ) = ψ ˜ R ( s + i k ν + θ R ) Ψ ˜ R ( s + i k ν + θ R ) p ˜ ^ R ( k , s ) ,
i ˜ ^ 0 ( k , s ) = ψ ˜ 0 ( s + θ 0 ) Ψ ˜ 0 ( s + θ 0 ) p ˜ ^ 0 ( k , s ) .
Inverse Fourier-Laplace transform of Equations (A15)–(A17) gives (33)–(35), respectively.

References

  1. Brady, S.T.; Morfini, G.A. Regulation of motor proteins, axonal transport deficits and adult-onset neurodegenerative diseases. Neurobiol. Dis. 2017, 105, 273–282. [Google Scholar] [CrossRef] [PubMed]
  2. Appert-Rolland, C.; Ebbinghaus, M.; Santen, L. Intracellular transport driven by cytoskeletal motors: General mechanisms and defects. Phys. Rep. 2015, 593, 1–59. [Google Scholar] [CrossRef]
  3. Hancock, W.O. Bidirectional cargo transport: Moving beyond tug of war. Nat. Rev. Mol. Cell Biol. 2014, 15, 615–628. [Google Scholar] [CrossRef] [PubMed]
  4. Schliwa, M.; Woehlke, G. Molecular motors. Nature 2003, 422, 759–765. [Google Scholar] [CrossRef]
  5. Vale, R.D. The molecular motor toolbox for intracellular transport. Cell 2003, 112, 467–480. [Google Scholar] [CrossRef] [PubMed]
  6. Hirokawa, N. Kinesin and dynein superfamily proteins and the mechanism of organelle transport. Science 1998, 279, 519–526. [Google Scholar] [CrossRef]
  7. Allan, V. One, two, three, cytoplasmic dynein is go! Science 2014, 345, 271–272. [Google Scholar] [CrossRef]
  8. Kardon, J.R.; Vale, R.D. Regulators of the cytoplasmic dynein motor. Nat. Rev. Mol. Cell Biol. 2009, 10, 854–865. [Google Scholar] [CrossRef]
  9. Bressloff, P.C.; Newby, J.M. Stochastic models of intracellular transport. Rev. Mod. Phys. 2013, 85, 135. [Google Scholar] [CrossRef]
  10. Briane, V.; Vimond, M.; Kervrann, C. An overview of diffusion models for intracellular dynamics analysis. Briefings Bioinform. 2020, 21, 1136–1150. [Google Scholar] [CrossRef]
  11. Smith, D.A.; Simmons, R.M. Models of motor-assisted transport of intracellular particles. Biophys. Journal 2001, 80, 45–68. [Google Scholar] [CrossRef]
  12. Bressloff, P.C. Stochastic Processes in Cell Biology; Springer: Berlin/Heidelberg, Germany, 2014; Volume 41. [Google Scholar]
  13. Jülicher, F.; Ajdari, A.; Prost, J. Modeling molecular motors. Rev. Mod. Phys. 1997, 69, 1269. [Google Scholar] [CrossRef]
  14. Klumpp, S.; Lipowsky, R. Cooperative cargo transport by several molecular motors. Proc. Natl. Acad. Sci. USA 2005, 102, 17284–17289. [Google Scholar] [CrossRef] [PubMed]
  15. Hafner, A.E.; Santen, L.; Rieger, H.; Shaebani, M.R. Run-and-pause dynamics of cytoskeletal motor proteins. Sci. Rep. 2016, 6, 37162. [Google Scholar] [CrossRef] [PubMed]
  16. Loverdo, C.; Bénichou, O.; Moreau, M.; Voituriez, R. Enhanced reaction kinetics in biological cells. Nat. Phys. 2008, 4, 134–137. [Google Scholar] [CrossRef]
  17. Newby, J.; Bressloff, P.C. Random intermittent search and the tug-of-war model of motor-driven transport. J. Stat. Mech. Theory Exp. 2010, 2010, P04014. [Google Scholar] [CrossRef]
  18. Chou, T.; Mallick, K.; Zia, R.K. Non-equilibrium statistical mechanics: From a paradigmatic model to biological transport. Rep. Prog. Phys. 2011, 74, 116601. [Google Scholar] [CrossRef]
  19. Schadschneider, A.; Chowdhury, D.; Nishinari, K. Stochastic Transport in Complex Systems: From Molecules to Vehicles; Elsevier: Amsterdam, The Netherlands, 2010. [Google Scholar]
  20. Pinkoviezky, I.; Gov, N. Transport dynamics of molecular motors that switch between an active and inactive state. Phys. Rev. E 2013, 88, 022714. [Google Scholar] [CrossRef]
  21. Ando, D.; Korabel, N.; Huang, K.C.; Gopinathan, A. Cytoskeletal network morphology regulates intracellular transport dynamics. Biophys. J. 2015, 109, 1574–1582. [Google Scholar] [CrossRef]
  22. Hafner, A.E.; Rieger, H. Spatial cytoskeleton organization supports targeted intracellular transport. Biophys. J. 2018, 114, 1420–1432. [Google Scholar] [CrossRef]
  23. Cherstvy, A.G.; Chechkin, A.V.; Metzler, R. Particle invasion, survival, and non-ergodicity in 2D diffusion processes with space-dependent diffusivity. Soft Matter. 2014, 10, 1591–1601. [Google Scholar] [CrossRef] [PubMed]
  24. Korabel, N.; Taloni, A.; Pagnini, G.; Allan, V.; Fedotov, S.; Waigh, T.A. Ensemble heterogeneity mimics ageing for endosomal dynamics within eukaryotic cells. Sci. Rep. 2023, 13, 8789. [Google Scholar] [CrossRef] [PubMed]
  25. Shaebani, M.R.; Sadjadi, Z.; Sokolov, I.M.; Rieger, H.; Santen, L. Anomalous diffusion of self-propelled particles in directed random environments. Phys. Rev. E 2014, 90, 030701. [Google Scholar] [CrossRef] [PubMed]
  26. Klein, S.; Appert-Rolland, C.; Santen, L. Fluctuation effects in bidirectional cargo transport. Eur. Phys. J. Spec. Top. 2014, 223, 3215–3225. [Google Scholar] [CrossRef]
  27. Korabel, N.; Waigh, T.A.; Fedotov, S.; Allan, V.J. Non-Markovian intracellular transport with sub-diffusion and run-length dependent detachment rate. PLoS ONE 2018, 13, e0207436. [Google Scholar] [CrossRef]
  28. Doerries, T.J.; Metzler, R.; Chechkin, A.V. Emergent anomalous transport and non-Gaussianity in a simple mobile–immobile model: The role of advection. New J. Phys. 2023, 25, 063009. [Google Scholar] [CrossRef]
  29. Doerries, T.J.; Chechkin, A.V.; Schumer, R.; Metzler, R. Rate equations, spatial moments, and concentration profiles for mobile-immobile models with power-law and mixed waiting time distributions. Phys. Rev. E 2022, 105, 014105. [Google Scholar] [CrossRef]
  30. Kurilovich, A.A.; Mantsevich, V.N.; Mardoukhi, Y.; Stevenson, K.J.; Chechkin, A.V.; Palyulin, V.V. Non-Markovian diffusion of excitons in layered perovskites and transition metal dichalcogenides. Phys. Chem. Chem. Phys. 2022, 24, 13941–13950. [Google Scholar] [CrossRef]
  31. Bouchaud, J.P.; Georges, A. Anomalous diffusion in disordered media: Statistical mechanisms, models and physical applications. Phys. Rep. 1990, 195, 127–293. [Google Scholar] [CrossRef]
  32. 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]
  33. Metzler, R.; Klafter, J. The restaurant at the end of the random walk: Recent developments in the description of anomalous transport by fractional dynamics. J. Phys. Math. Gen. 2004, 37, R161. [Google Scholar] [CrossRef]
  34. Barkai, E.; Garini, Y.; Metzler, R. Strange kinetics of single molecules in living cells. Phys. Today 2012, 65, 29–35. [Google Scholar] [CrossRef]
  35. Waigh, T.A. The Physics of Living Processes: A Mesoscopic Approach; John Wiley & Sons: Hoboken, NJ, USA, 2014. [Google Scholar]
  36. Zaburdaev, V.; Denisov, S.; Klafter, J. Lévy walks. Rev. Mod. Phys. 2015, 87, 483. [Google Scholar] [CrossRef]
  37. Metzler, R.; Jeon, J.H.; Cherstvy, A.G.; Barkai, E. Anomalous diffusion models and their properties: Non-stationarity, non-ergodicity, and ageing at the centenary of single particle tracking. Phys. Chem. Chem. Phys. 2014, 16, 24128–24164. [Google Scholar] [CrossRef] [PubMed]
  38. Waigh, T.A.; Korabel, N. Heterogeneous anomalous transport in cellular and molecular biology. Rep. Prog. Phys. 2023; in press. [Google Scholar]
  39. Salman, H.; Gil, Y.; Granek, R.; Elbaum, M. Microtubules, motor proteins, and anomalous mean squared displacements. Chem. Phys. 2002, 284, 389–397. [Google Scholar] [CrossRef]
  40. Caspi, A.; Granek, R.; Elbaum, M. Diffusion and directed motion in cellular transport. Phys. Rev. E 2002, 66, 011916. [Google Scholar] [CrossRef]
  41. Weiss, M.; Elsner, M.; Kartberg, F.; Nilsson, T. Anomalous subdiffusion is a measure for cytoplasmic crowding in living cells. Biophys. J. 2004, 87, 3518–3524. [Google Scholar] [CrossRef]
  42. Kulkarni, R.P.; Castelino, K.; Majumdar, A.; Fraser, S.E. Intracellular transport dynamics of endosomes containing DNA polyplexes along the microtubule network. Biophys. J. 2006, 90, L42–L44. [Google Scholar] [CrossRef]
  43. Kulić, I.M.; Brown, A.E.; Kim, H.; Kural, C.; Blehm, B.; Selvin, P.R.; Nelson, P.C.; Gelfand, V.I. The role of microtubule movement in bidirectional organelle transport. Proc. Natl. Acad. Sci. USA 2008, 105, 10011–10016. [Google Scholar] [CrossRef]
  44. Bruno, L.; Levi, V.; Brunstein, M.; Desposito, M.A. Transition to superdiffusive behavior in intracellular actin-based transport mediated by molecular motors. Phys. Rev. E 2009, 80, 011912. [Google Scholar] [CrossRef] [PubMed]
  45. Robert, D.; Nguyen, T.H.; Gallet, F.; Wilhelm, C. In vivo determination of fluctuating forces during endosome trafficking using a combination of active and passive microrheology. PLoS ONE 2010, 5, e10046. [Google Scholar] [CrossRef] [PubMed]
  46. Tabei, S.A.; Burov, S.; Kim, H.Y.; Kuznetsov, A.; Huynh, T.; Jureller, J.; Philipson, L.H.; Dinner, A.R.; Scherer, N.F. Intracellular transport of insulin granules is a subordinated random walk. Proc. Natl. Acad. Sci. USA 2013, 110, 4911–4916. [Google Scholar] [CrossRef] [PubMed]
  47. Chen, K.; Wang, B.; Granick, S. Memoryless self-reinforcing directionality in endosomal active transport within living cells. Nat. Mater. 2015, 14, 589–593. [Google Scholar] [CrossRef]
  48. Reverey, J.F.; Jeon, J.H.; Bao, H.; Leippe, M.; Metzler, R.; Selhuber-Unkel, C. Superdiffusion dominates intracellular particle motion in the supercrowded cytoplasm of pathogenic Acanthamoeba castellanii. Sci. Rep. 2015, 5, 11690. [Google Scholar] [CrossRef]
  49. Song, M.S.; Moon, H.C.; Jeon, J.H.; Park, H.Y. Neuronal messenger ribonucleoprotein transport follows an aging Lévy walk. Nat. Commun. 2018, 9, 344. [Google Scholar] [CrossRef]
  50. Flores-Rodriguez, N.; Rogers, S.S.; Kenwright, D.A.; Waigh, T.A.; Woodman, P.G.; Allan, V.J. Roles of dynein and dynactin in early endosome dynamics revealed using automated tracking and global analysis. PLoS ONE 2011, 6, e24479. [Google Scholar] [CrossRef]
  51. Han, D.; Korabel, N.; Chen, R.; Johnston, M.; Gavrilova, A.; Allan, V.J.; Fedotov, S.; Waigh, T.A. Deciphering anomalous heterogeneous intracellular transport with neural networks. Elife 2020, 9, e52224. [Google Scholar] [CrossRef]
  52. Korabel, N.; Han, D.; Taloni, A.; Pagnini, G.; Fedotov, S.; Allan, V.; Waigh, T.A. Local analysis of heterogeneous intracellular transport: Slow and fast moving endosomes. Entropy 2021, 23, 958. [Google Scholar] [CrossRef]
  53. Banks, D.S.; Fradin, C. Anomalous diffusion of proteins due to molecular crowding. Biophys. J. 2005, 89, 2960–2971. [Google Scholar] [CrossRef]
  54. Weiss, M. Crowding, diffusion, and biochemical reactions. Int. Rev. Cell Mol. Biol. 2014, 307, 383–417. [Google Scholar] [PubMed]
  55. Waigh, T.A. Microrheology of complex fluids. Rep. Prog. Phys. 2005, 68, 685. [Google Scholar] [CrossRef]
  56. Waigh, T.A. Advances in the microrheology of complex fluids. Rep. Prog. Phys. 2016, 79, 074601. [Google Scholar] [CrossRef] [PubMed]
  57. 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]
  58. Zumofen, G.; Klafter, J. Laminar–localized-phase coexistence in dynamical systems. Phys. Rev. E 1995, 51, 1818. [Google Scholar] [CrossRef]
  59. Thiel, F.; Schimansky-Geier, L.; Sokolov, I.M. Anomalous diffusion in run-and-tumble motion. Phys. Rev. E 2012, 86, 021117. [Google Scholar] [CrossRef]
  60. Portillo, I.G.; Campos, D.; Méndez, V. Intermittent random walks: Transport regimes and implications on search strategies. J. Stat. Mech. Theory Exp. 2011, 2011, P02033. [Google Scholar] [CrossRef]
  61. Klafter, J.; Sokolov, I.M. First Steps in Random Walks: From Tools to Applications; OUP Oxford: Oxford, UK, 2011. [Google Scholar]
  62. 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]
  63. Fedotov, S.; Han, D.; Ivanov, A.O.; Da Silva, M.A. Superdiffusion in self-reinforcing run-and-tumble model with rests. Phys. Rev. E 2022, 105, 014126. [Google Scholar] [CrossRef]
  64. Al Shamsi, H. Migration and Proliferation Dichotomy: A Persistent Random Walk of Cancer Cells. Fractal Fract. 2023, 7, 318. [Google Scholar] [CrossRef]
  65. 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]
  66. Mendez, V.; Fedotov, S.; Horsthemke, W. Reaction-Transport Systems: Mesoscopic Foundations, Fronts, and Spatial Instabilities; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  67. Fedotov, S. Single integrodifferential wave equation for a Lévy walk. Phys. Rev. E 2016, 93, 020101. [Google Scholar] [CrossRef] [PubMed]
  68. Giona, M.; Cairoli, A.; Klages, R. Extended Poisson-Kac theory: A unifying framework for stochastic processes with finite propagation velocity. Phys. Rev. X 2022, 12, 021004. [Google Scholar] [CrossRef]
  69. Cox, D.R.; Miller, H.D. The theory of Stochastic Processes; Routledge: Cambridge, MA, USA, 2017. [Google Scholar]
  70. 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] [PubMed]
  71. Henry, B.; Langlands, T.; Wearne, S. Anomalous diffusion with linear reaction dynamics: From continuous time random walks to fractional reaction-diffusion equations. Phys. Rev. E 2006, 74, 031116. [Google Scholar] [CrossRef] [PubMed]
  72. Angstmann, C.N.; Donnelly, I.C.; Henry, B.I. Continuous time random walks with reactions forcing and trapping. Math. Model. Nat. Phenom. 2013, 8, 17–27. [Google Scholar] [CrossRef]
  73. 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]
  74. 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]
  75. Iomin, A. Continuous time random walk and migration–proliferation dichotomy of brain cancer. Biophys. Rev. Lett. 2015, 10, 37–57. [Google Scholar] [CrossRef]
  76. Pogorui, A.A.; Rodríguez-Dagnino, R.M. Isotropic random motion at finite speed with K-Erlang distributed direction alternations. J. Stat. Phys. 2011, 145, 102–112. [Google Scholar] [CrossRef]
  77. Taylor-King, J.P.; Klages, R.; Fedotov, S.; Van Gorder, R.A. Fractional diffusion equation for an n-dimensional correlated Lévy walk. Phys. Rev. E 2016, 94, 012104. [Google Scholar] [CrossRef] [PubMed]
  78. Santra, I.; Basu, U.; Sabhapandit, S. Run-and-tumble particles in two dimensions: Marginal position distributions. Phys. Rev. E 2020, 101, 062120. [Google Scholar] [CrossRef] [PubMed]
  79. Di Crescenzo, A.; Iuliano, A.; Mustaro, V. On some finite-velocity random motions driven by the geometric counting process. J. Stat. Phys. 2023, 190, 44. [Google Scholar] [CrossRef]
Figure 1. Markovian three–state model. (a) The distribution of random walkers is calculated at T m a x = 10 with zero removal rates Q A = Q R = Q 0 = 0 (circles) and Q A = Q R = Q 0 = 0.1 (squares). Other parameters are given in the text. The dashed line represents the Gaussian distribution with σ 2 = 2 D T m a x . (b) The corresponding second moment m 2 ( t ) as a function of time (blue and red curves) grows ballistically for short times and switches to linear growth with the effective diffusion coefficient D.
Figure 1. Markovian three–state model. (a) The distribution of random walkers is calculated at T m a x = 10 with zero removal rates Q A = Q R = Q 0 = 0 (circles) and Q A = Q R = Q 0 = 0.1 (squares). Other parameters are given in the text. The dashed line represents the Gaussian distribution with σ 2 = 2 D T m a x . (b) The corresponding second moment m 2 ( t ) as a function of time (blue and red curves) grows ballistically for short times and switches to linear growth with the effective diffusion coefficient D.
Fractalfract 07 00758 g001
Figure 2. Non–Markovian three–state model. The first moment m 1 ( t ) of random walkers position as a function of time calculated with zero removal rates Q A = Q R = Q 0 = 0 , λ A = λ R = λ 0 = 1 / 3 , μ = 0.3 , μ = 0.7 , τ 0 = 1 , λ = 1 and v + = 1 , v = 0.1 . The dashed lines correspond to the analytical Equation (62) without fitting.
Figure 2. Non–Markovian three–state model. The first moment m 1 ( t ) of random walkers position as a function of time calculated with zero removal rates Q A = Q R = Q 0 = 0 , λ A = λ R = λ 0 = 1 / 3 , μ = 0.3 , μ = 0.7 , τ 0 = 1 , λ = 1 and v + = 1 , v = 0.1 . The dashed lines correspond to the analytical Equation (62) without fitting.
Fractalfract 07 00758 g002
Figure 3. Non–Markovian three–state model. (a) The distribution of random walkers is calculated at T m a x = 10 with zero removal rates Q A = Q R = Q 0 = 0 (circles) and Q A = Q R = Q 0 = 0.1 (squares), τ 0 = 1 , λ = 1 and v + = 1 , v = 0.1 . The dashed line shows the Gaussian distribution for comparison. (b) The corresponding mean squared displacements (blue and red curves) grow ballistically for short times and switch to sub-linear growth m 2 ( t ) t β with β = μ and the exponent μ = 0.7 .
Figure 3. Non–Markovian three–state model. (a) The distribution of random walkers is calculated at T m a x = 10 with zero removal rates Q A = Q R = Q 0 = 0 (circles) and Q A = Q R = Q 0 = 0.1 (squares), τ 0 = 1 , λ = 1 and v + = 1 , v = 0.1 . The dashed line shows the Gaussian distribution for comparison. (b) The corresponding mean squared displacements (blue and red curves) grow ballistically for short times and switch to sub-linear growth m 2 ( t ) t β with β = μ and the exponent μ = 0.7 .
Fractalfract 07 00758 g003
Figure 4. Three–state model with ML–distributed resting times can have subdiffusive or superdiffusive asymptotic behaviour. (a) The dependence β / μ 1 as a function of the speed v + (with v kept constant) for the three-state model with ML-distributed resting state. The dashed line corresponds to Equation (71) and symbols show the results of numerical simulations. (b) The diagram of the asymptotic behaviour of different models is given by the dependence of exponents of MSDs or variances at long times, β , as functions of the exponent of ML distribution of resting states, μ . The line with the biggest slope corresponds to the two-state model (that is, the three-state model with either v + = 0 or v = 0 ) with ML-distributed resting states. It has subdiffusive behaviour β < 1 for μ < 1 / 2 and superdiffusive behaviour β > 1 for μ > 1 / 2 . The lowest line corresponds to the three-state model with equal retrograde and anterograde cargo speeds v + = v . The line with an intermediate slope corresponds to the three-state model with unequal retrograde and anterograde cargo speeds, v + = 1 and v = 0.25 .
Figure 4. Three–state model with ML–distributed resting times can have subdiffusive or superdiffusive asymptotic behaviour. (a) The dependence β / μ 1 as a function of the speed v + (with v kept constant) for the three-state model with ML-distributed resting state. The dashed line corresponds to Equation (71) and symbols show the results of numerical simulations. (b) The diagram of the asymptotic behaviour of different models is given by the dependence of exponents of MSDs or variances at long times, β , as functions of the exponent of ML distribution of resting states, μ . The line with the biggest slope corresponds to the two-state model (that is, the three-state model with either v + = 0 or v = 0 ) with ML-distributed resting states. It has subdiffusive behaviour β < 1 for μ < 1 / 2 and superdiffusive behaviour β > 1 for μ > 1 / 2 . The lowest line corresponds to the three-state model with equal retrograde and anterograde cargo speeds v + = v . The line with an intermediate slope corresponds to the three-state model with unequal retrograde and anterograde cargo speeds, v + = 1 and v = 0.25 .
Fractalfract 07 00758 g004
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

Korabel, N.; Al Shamsi, H.; Ivanov, A.O.; Fedotov, S. Non-Markovian Persistent Random Walk Model for Intracellular Transport. Fractal Fract. 2023, 7, 758. https://doi.org/10.3390/fractalfract7100758

AMA Style

Korabel N, Al Shamsi H, Ivanov AO, Fedotov S. Non-Markovian Persistent Random Walk Model for Intracellular Transport. Fractal and Fractional. 2023; 7(10):758. https://doi.org/10.3390/fractalfract7100758

Chicago/Turabian Style

Korabel, Nickolay, Hamed Al Shamsi, Alexey O. Ivanov, and Sergei Fedotov. 2023. "Non-Markovian Persistent Random Walk Model for Intracellular Transport" Fractal and Fractional 7, no. 10: 758. https://doi.org/10.3390/fractalfract7100758

APA Style

Korabel, N., Al Shamsi, H., Ivanov, A. O., & Fedotov, S. (2023). Non-Markovian Persistent Random Walk Model for Intracellular Transport. Fractal and Fractional, 7(10), 758. https://doi.org/10.3390/fractalfract7100758

Article Metrics

Back to TopTop