1 Fracional Transpor o Tumor Cells ALEXANDER IOMIN Deparmen o Physics Technion Israel Insiue o Technology Technion Ciy Haia 3 ISRAEL Absrac: - The inl...

Author:
Claud Rich

0 downloads 98 Views 197KB Size

1 Introduction Mathematical modeling of tumor development is an important and new application of mathematical physics in biology and medicine. Although it is mainly aimed at diagnostics and treatments of cancers, the importance of tumor modeling for the understanding of cancer cell transport cannot be overestimated. Recent surveys describe different aspects of the modeling of tumors including solid tumors [1,2,3] interacting with the immune system [1], diffusive models related e.g. to brain tumors [4], process of tumor induced vascularization [5,6], and fractal geometry of pathological architecture of tumor [7], as well as chemotherapy strategies [1,2,3,4,6]. Tumor development consists of complicated processes with different stages (see e.g. [1]) where the tumor’s cell transport and their proliferation are the main contributors to the malignant neoplasm dissemination. Interplay between these two main processes of cell proliferation and transport leads to the essential complication of the mathematical modeling of the tumor growth [1,5]. In the present study, we focus primarily on the influence of the cell fission on transport properties through vessel network. It could be either vascular or lymphatic net. Since we do not specify a kind of tumor, we do not specify a kind of vessels. Our primary interest is concerned with the main stages of tumor development, which are cell fission and transport. A simple mathematical model of a

continuous time random walk (CTRW) [8] is proposed, using heuristic arguments on tumor development due to these two main stages. The constructed model is a modification of a so–called comb structure [9,10,11,12]. By virtue of this model we are able to show that the tumor development corresponds to fractional transport whose mathematical apparatus is well established (see e.g. [13,14,15,16,17]). Using this simplified approach of fractional transport, a possible answer to the question how the neoplasm cells appear arbitrarily far from the main (primary) tumor in the case of solid tumor [3] is proposed. The model can be considered as a possible mechanism for diffusive cancers [4] as well.

2 Main Assumptions for the CTRW First, we consider a simplified scheme of cell dissemination through the vessel network. We consider this process by means of the following two steps. The first step is the biological process of cell fission. The duration of this stage is T f . The second process is cell transport itself having a duration Tt . Therefore the cell dissemination is approximately characterized by the fission time T f and the transport time Tt . During the time scale T f the cells interact strongly [18], motility of the cells is weak, and there is no transport (approximately). The duration of T f could be arbitrarily large, and it reaches 107 ÷ 108 sec [19]. During the second time

Tt interaction between the cells is weak and

Laplace

motility of the cells is determined by the velocity V of either vascular or lymphatic flow through the vessel network. It is convenient to introduce a “jump” length X t as the distance which a cell

such that the Fourier-Laplace transforms could be applied

travels during the time Tt

X t = VTt .

(1) Hence, the cells form an initial packet of free spreading particles and the tumor development process consists of the following time consequences (2) T f (1)Tt (2)T f (3)…. There are different realizations of this chain of times, due to different duration of T f (i ) and Tt (i ) , where i = 1, 2,…. Therefore one comes to the conclusion that transport is characterized by random values T (i ) which are waiting times between any two successive jumps of random length X (i ) . This phenomenon is known as a continuous time random walk (CTRW) [8,13,14,15]. It arises as a result of a sequence of independent identically distributed random waiting times T (i ) , each having the same probability density function (PDF) w(t ), t > 0 with a mean waiting time ∞

T ≡ 〈t 〉 = ∫ tw(t )dt , 0

and a sequence of independent identically distributed random jumps, x = X (i ) , each having the same PDF λ (x) with the jump length variance ∞

σ 2 ≡ 〈 x 2 〉 = ∫ x 2 λ ( x)dx . −∞

Now we introduce the PDF P ( x, t ) of the particle to be in point x at the time t . Due to the probabilistic description that defines an appropriate relation between these three PDFs, P( x, t ), w(t ), λ ( x) (see e.g. [13,15]), one obtains the following integral equation for P ( x, t ) : ∞

P ( x, t ) = δ ( x) ∫ w(t ′)dt ′ + t

∫

∞

0

∞

w(t − t ′) ∫ λ (x − x′) P ( x′, t ′)dt ′dx′ −∞

(3)

with the initial condition P ( x, 0) = δ ( x) . It is worth stressing that the Fourier and Laplace ( Fˆ − Lˆ ) transforms play an important role in the CTRW, since a simple form relation between P( x, t ) and w(t ), λ ( x) takes place in the Fourier-

space

[8,13,14,15].

Suppose

that

P( x, t ), w(t ), λ ( x) are well behaved functions,

w( p ) = Lˆ[ w(t )]; λ (k ) = Fˆ [λ ( x)]; ˆ ˆ[ P( x, t )]. P(k , p) = FL

(4) Then we deduce the integral equation (3) to the Montroll-Weiss equation [13,15]:

P(k , p) =

1 − w( p) 1 ⋅ . p 1 − w( p)λ (k )

(5)

This is the main result in which we were able to establish a link between the tumor development and the CTRW process which is described by equation (5). In sequel we consider some examples of the CTRW dynamics that could be applied for different realization of tumor cell transport.

3 CTRW Equation First, we present familiar examples of fractional transport used for a variety of realizations in physics, chemistry, biology and so on (see e.g. recent surveys [13,17]). These examples are also relevant to the tumor development. We consider a situation when σ 2 is finite and corresponds to the following distribution

λ ( x) =

2 2 1 e − x / 4σ , 2 1/ 2 [4πσ ]

(6)

while T diverges and is described by a long–tailed waiting time PDF with an asymptotic behavior The Laplace w(t ) ∼ Aα (T /t )α , 0 < α < 1 . transform is w( p ) ∼ 1 − ( pT )α . Therefore, Eq. (5) reads P(k , p) = 1/ ( p + Dα p1−α k 2 ) , (7) where the generalized diffusion constant is now Dα = σ 2 /T α . The MSD is calculated from (7) via the following relation

〈 x 2 (t )〉 = Lˆ −1 lim[−(d 2 /dk 2 ) P(k , p)] , k →0

where Lˆ in

−1

means the Laplace inversion. It results

〈 x 2 (t )〉 = 2 Dα t α /Γ(1 + α ) , (8) where Γ( z ) is a gamma function [20]. Since α < 1 , this is subdiffusion. This is an example for tumor development with different rates of cell dissemination through the fractional net of vessels embedded in the three– dimensional (3d) space. This is relevant, e.g., for

description of both a diffusive cancer and a primary solid tumor.

3.1 Fractional Equation To obtain the fractional or CTRW equation which produces the solution (7), we introduced here the Riemann–Liouville fractional derivatives (see, for example, [13,16])

∂α ∂tα

f (t ) by means of the

Laplace transform ( 0 < α < 1 ): α

∂ (9) Lˆ[ α f (t )] = pα f ( p) − p1−α f (0+ ), ∂t where Lˆ[ f (t )] = f ( p ) , and it also implies ∂α [1]/∂t α = 0 [16]. Using this definition, we write the CTRW equation in the following form

∂P ∂1−α + α 1−α Lˆ FP ( x) P = 0 , ∂t ∂t (10) where the Fokker–Planck operator is diffusion

Lˆ FP ( x) = − D

∂2 . ∂x 2

(11)

For α = 1/ 2 the situation is the most simple, and the traps can be modeled by normal diffusion in the additional y direction. Therefore the fraction equation (10) corresponds to a so–called comb model.

Figure 1. A comb structure In this section we consider homogeneous convection in the x direction with a velocity v( y ) = V δ ( y ) instead of diffusion. A random walk on the comb structure is described by the distribution function P1 = P1 (t , x, y ) and the current

j = (−v( y ) P1 , − D

∂P1 ), ∂y

where ∞

P (t , x) = ∫ P1 (t , x, y ) dy , −∞

(12)

The Liouville equation

∂P1 + div j = 0 ∂t

(13)

yields the following Fokker–Planck equation

3.2 Comb Model Fractional transport of cells, namely subdiffusion, could be described in the framework of the comb model (or CTRW structure) [9]. The comb model shown in Fig. 1 is an example of subdiffusive 1d media where CTRW takes place along the x structure axis. Diffusion in the y direction plays the role of traps with the PDF of delay times of the form w(t ) ∼ 1/ (1 + t /T )3 / 2 . A special behavior of diffusion on the comb structure is that the displacement in the x –direction is possible only along the structure axis ( x -axis at y = 0 ). Therefore, cell motility is highly inhomogeneous in the y -direction, while diffusion coefficient in the

y –direction along the teeth is a constant Dyy = D .

∂P1 ∂P ∂2P − V δ ( y ) 1 − D 21 = 0 (14) ∂t ∂x ∂y with the initial conditions P1 (0, x, y ) = δ ( x)δ ( y ) and the boundary conditions on the infinities P1 (t , ±∞, ±∞) = P1′ (t , ±∞, ±∞) = 0 . Here and in sequel the primes mean the spatial derivatives. Applying the Laplace and Fourier transforms, this equation is solved exactly with the solution

P1 (t , x, 0) =

D1/ 2 x 2 0

v

πt

3

exp( − Dx 2 /V 2t ) ,

(15)

and P1 = 0 for x < 0 , since the distribution function must be positive. This solution describes diffusion of cells in the convection media with traps. It corresponds to the normal diffusion with the second moment

〈 x 2 (t )〉 =

V2 t, D

(16)

but the effective diffusion coefficient V 2 /D is determined by the external convective forcing V which is defined e.g. in (1). It is worth stressing that it is a nontrivial result, and therefore one should anticipate superdiffusion due to the inhomogeneous convection of the form V → vx s with s > 0 , and, correspondingly, subdiffusion when s < 0 [11,12].

4 Proliferation of cells In this section we consider a possible mechanism of tumor cell proliferation in the framework of the comb model. The total number of the transporting cells, ∞

N (t ) = ∫ dxP (t , x) ,

(17)

−∞

described by the Fokker–Planck equation (14) is conserved. Nevertheless, the process of proliferation should be taken into account not only by counting the waiting time PDF but also due to the fact that number of cells is not conserved dN /dt > 0 . Since, according the CTRW model, the transporting cells along the x axis do not proliferate, we introduce the proliferation rate as a change of the total number of cells with time dN /dt = g (t ) , where g (t ) is taken from empirical (clinical) data. It is convenient to present g (t ) as an integration over entire configuration space ( x, y ) . Then we have ∞ ∞ dN = g (t ) = ∫ ∫ dxdyG ( x, y, t ) , (18) −∞ −∞ dt where G ( x, y, t ) is an arbitrary function which

satisfies (18). For example, one could present it in the following convenient form

G ( x, y, t ) =

1 exp −( x 2 + y 2 ) / 4 Dg (t ) . (19) 4π D

Substituting Eqs. (18) and (19) in Eq. (14) one obtains

∂P1 ∂2P + δ ( y ) Lˆ FPP1 − D 21 = G ( x, y, t ) , ∂t ∂y

(20)

where the boundary conditions remain the same, but the initial condition is P (t = 0, x, y ) = 0 . This condition means that the population of cells in the system is just due to the proliferation. For the sake of clarity and simplicity we take here that

g (t ) = bt ,

while

∂ Lˆ FP ( x) = − D ∂x2 , 2

which

corresponds to the standard comb model [9,10,11]. Performing the Laplace transform, we obtain the solution in the following form:

P1 ( p, x, y ) = f ( p, x) exp[− | y |

p /D ] .

(21) Then integrating both the sides of the equation over

y;

∫

∞

−∞

dy , we arrive at the equation for

f ≡ f ( p, x) Df ′′ + 2 pD f = −

b d e −| x| p/D . 2 D1/ 2 dp p

(22)

Performing the Fourier transform of Eq. (22), we obtain, eventually, the solution in the Fourier– Laplace domain

f ( p, k ) =

2b 1 ⋅ . 2 2 2 π ( p + Dk ) Dk + 2 Dp

(23)

Since solution (23) is the product of two functions where f ( p, k ) = f 1( p, k ) f 2( p, k ) , ˆ ˆ f (t , x ) , one can use the properties f 1(2) = LF 1(2) of the Fourier and Laplace transforms for the convolutions. Therefore we obtain for the Laplace inversion −1 P1(t , k , 0) = Lˆ f ( p, k ) =

∫

t

0

f1 (t − τ , k ) f 2 (τ , k )dτ .

(24)

Here the functions f1(2) are

f1 (t , k ) = t exp[ Dtk 2 ] ,

(25)

1 f 2 (t , k ) = − 4π Dt Dk 2t1/ 2 2k 4t Dk 2 exp D erfc , (26) 1/ 2 4D 4D 2D where erfc ( z ) is the error function [20]. To perform the Fourier inverse transform, we consider two limits of (26). The first limit is the large scale asymptotic, when k → 0 and k 2t 1, ∀ t . In this case, only the first term in (26) should be taken into consideration. The Fourier inversion is carried out exactly. One arrives at the following expression:

P1 (t , x, 0) =

∫

t

0

dτ

x2 t2 exp − . x2 4 Dt

t −τ x2 exp − ∝ DDt 4 D (t − τ ) (27)

This result is valid for x D /D . As follows from Eq. (27), the rate of the cells dissemination on this

large scale asymptotic is of the order of 〈 x 2 (t )〉 ∼ t 5 / 2 . The second limit is k 4t 1 . In this case Eq. (26) reads f 2 (t , k ) ≈

4D 4 D 2k 4 t 3 / 2

. The Fourier

inversion is carried out by the stationary phase approximation that gives

∫

∞

−∞

1 − ikx − Dtk 2 16 D 7 / 2t 7 / 2 x2 − ≈ e e dk e xp[ ]. k4 4 Dt 2 x4

numerous and helpful discussions. This research was supported in part by the Israel Science Foundation and the Minerva Center for Nonlinear Physics of Complex Systems.

References: (28)

Finally, we obtain for x < D/D

32 D 3D (t − τ )9 / 2 0 t 3/ 2 x 4 x2 x2 t5 exp × exp − ∝ − 4 Dt . x6 4 D (t − τ )

Acknowledgments: - I thank E. Baskin for

t

P(t , x, 0) = ∫ dτ

(29)

Since 〈 x 2 (t )〉 ∼ ∞ , hence in this short range,

x < D/D , the cells spread due to proliferation and tumor size is determined by transport on the size scale larger than D/D .

5 Conclusion The main focus of the present study has been the influence of cell proliferation on transport properties through vessel’s network (either vascular or lymphatic). Two main stages have been highlighted: cell fission and transport with durations T f and Tt , correspondingly. Using these time scales, we were able to deduce a description of the tumor development to a CTRW process. A simple mathematical model is constructed, using heuristic arguments on the relation of tumor development and the CTRW whose mathematical apparatus is well established (see e.g. [13,14,15,17]. The constructed model is a modification of a so–called comb structure [9,10,11,12]. Using this simplified approach to the fractional transport of tumor cells, we can answer the question how the malignant neoplasm cells spread for both solid tumors and for diffusive cancers. We presented analytical solutions of the problem. To this end we consider one dimensional (1d) transport. A generalization of the analytical consideration on the 3d case is straightforward when an interaction between the degrees of freedom is absent. An important feature considered in the framework of the CTRW model is an essential enhancing of anomalous transport due to proliferation. Moreover in some cases it could be dominant.

[1] N. Bellomo and L. Preziosi, Math. Comp. Model., 2000, 413. [2] D. Ambrosi and L. Preziosi, Math. Models Methods Appl. Sci., 2002, 737. [3] S.E. Clare, F. Nakhils and J.C. Panetta, Breast Cancer Res. , 2000, 430. [4] K.R. Swanson, S. Bridge, J.D. Murray and E.S. Alvord Jr, J. Neurol. Scien., 2003, 1. [5] S.R. McDougall, A.R.A. Anderson, M.A.J. Chaplain and J.A. Sherratt, Bull. Math. Biol., 2002, 673. [6] L. Arakelyan, Z. Agur and V. Vainstein, Angiogenesis, 2002, 203. [7] J.W. Bahish, and R.K. Jain, Cancer Res., 2000, 3683. [8] E.W. Montroll and G.H. Weiss, J. Math. Phys., 1965, 167. [9] G.H. Weiss and S. Havlin, Physica A, , 1984, 474. [10] V.E. Arkhincheev and E.M. Baskin, Sov. Phys. JETP, 1991, 161. [11] E.M. Baskin and A. Iomin, Phys. Rev. Lett., 2004, 120603 [12] A. Iomin and E.M. Baskin, submitted to Phys. Rev. E, cond–mat/0405089. [13] R. Metzler and J. Klafter, Phys. Rep., 2000, 1. [14] E.W. Montroll and M.F. Shlesinger, in Studies in Statistical Mechanics, v. 11, Noth–Holland, Amsterdam 1984. [15] R. Gorenflo and F. Mainardi, in Processes with long range Correlations, Springer Verlag, Berlin, 2003. [16] F. Mainardi, Chaos, Sol. and Fract., 1996, 1461. [17] R. Hilfer, (editor) Fractional Calculus in Physics, World Scientific, Singapore, 2000. [18] A.A. Tempia-Caleira, et al, J. Surgical Oncology, , 2002, 93. [19] V. Kumar, R.S. Cotran and S.L. Robbins, Basic Pathology, Elsevier, Philadelphia, 2003. [20] E. Janke, F. Emde and F. Lösh, Tables of Higher Functions, McGraw-Hill, New York, 1960.

Our partners will collect data and use cookies for ad personalization and measurement. Learn how we and our ad partner Google, collect and use data. Agree & Close