On-line version ISSN 1678-4448
Braz. J. Phys. vol.38 no.2 São Paulo June 2008
The effect of correlated noise in a Gompertz tumor growth model
Anita BeheraI; S. Francesca C. O'RourkeII
IQueen's University of Belfast, BT9 7AB, United Kingdom
IISchool of Mathematics & Physics Queen's University of Belfast, BT7 1NN, United Kingdom
We study the effect of noise in an avascular tumor growth model. The growth mechanism we consider is the Gompertz model. The steady state probability distributions and average population of tumor cells are analyzed within the Fokker-Planck formalism to investigate the importance of additive and multiplicative noise. We consider the effect of correlation on tumor growth for both the case of nonzero and zero correlation time. It is observed that the Gompertz model, driven by correlated noise exhibits a stochastic resonance and phase transition. This behaviour is attributed to multiplicative noise. In the case of nonzero correlation time, it is found that the correlation strength and correlation time have opposite effects on the steady state probability distribution. The Gompertz model simulations are also shown to be in qualitative agreement with another similiar non-bistable system, the logistic model.
Keywords: Fluctuation phenomena; Random processes; Noise; Brownian motion
Much of the attention in the recent years has been directed towards nonlinear physics and its application to uncover biological complexities. Studies have confirmed the role of noise in the nonlinear stochastic systems . In recent years the Fokker-Planck equation has become one of the important approaches in the studies of nonlinear dynamics based on stochastic form [2,3].
Tumor growth is a complex process. It has been a challenge for many years to search for a suitable growth law of tumors. Mathematical models based on mathematical equations such as the Gompertzian, logistic and exponential equations are used as a basic tool for describing avascular tumor growth [4-6]. It has been realized that tumor growth is governed by environmental fluctuations . It has also been shown that quite often there are discrepancies between theoretical predictions and clinical data due to more or less intense environmental fluctuations . Ferreira et al  analyzed the effect of distinct chemotheraputic strategies for the growth of avascular tumors. This study confirmed that an environment like chemotherapy affects tumor growth behaviour and may lead to morphological transitions under certain conditions.
Stochastic modelling using the Fokker-Planck equation has been recently applied to the logistic equation to study the steady state properties of avascular tumor growth driven by Gaussian correlated additive and multiplicative white noise by . This work has been extended to consider the nonzero correlation time between additive and multiplicative noise by . However among avascular tumor growth power laws the Gompertz model has been the most broadly and successfully applied to fit the experimental data [11-13]. Other successful applications in the literature of the Gompertz equation to cancer cell growth include [14-16]. These models are however deterministic models. In this article we wish to consider the Gompertz model from a stochastic viewpoint by examining the steady state properties of Gompertzian tumor cell growth driven by correlated noise. The Gompertz model has not been analyzed in the literature before in this context. In addition we extend our study to consider the nonzero correlation time for the Gompertz model.
The plan of this paper is as follows. Calculational details outlining the theoretical formulation of the Fokker-Planck equation applied to avascular tumor growth laws are presented in Sec.II. Sec. III presents the results for the steady state properties and average cell populations of the Gompertz model driven by cross-correlated additive and multiplicative noise for the case of zero correlation time between noises. In Sec. IV we consider nonzero correlation time and present the results for the steady state properties of the Gompertz model driven by cross-correlated noise for the case of nonzero correlation time. We use our model to simulate the behaviour of the logistic model of , , and  and compare with the results obtained from the Gompertz model.
Finally in Sec. V we discuss our conclusions.
II. THEORETICAL FORMULATION BASED ON THE FOKKER-PLANCK EQUATION
Avascular tumor growth laws may be described by a single deterministic differential equation,
where f(x) describes the tumor cell growth dynamics and x0 is the initial tumor size.
Typically the modelling assumptions that Eq.(1) are based on, are (i) the tumor only contains one cell type, (ii) is spatially independent, (iii) does not explicitly mention nutrients, growth factors or host vasculature, (iv) tumor volume is proportional to x(t), the number of tumor cells at time t. The Gompertz law may be modelled by taking
where x is the tumor cell number, b the cell decay rate, k is carrying capacity, where k = , and a is the cell growth rate. Gompertz growth is a result of two classes of competitive processes, the first process simulates growth and the second phase constrains growth at the saturation stage. Another alternative growth law which behaves similarly to the Gompertz model is the logistic model. The logistic model is governed by the equation where,
The solution of this equation belongs to same class of sigmoidal functions as the Gompertz model. It models growth exponentially in the early stages but eventually saturates due to the quadratic term in the above equation.
Eq. (1) can be generalized to consider stochastic effects due to external factors such as temperature, drugs, radiotherapy etc, by introducing, (i) Gaussian multiplicative noise to represent the effect of the treatment by altering the tumor cell dynamics (ii) a negative additive Gaussian noise which may represent fluctuations due to the treatment resulting in cell death. Here we have implicitly assumed both the multiplicative and additive noise are correlated since they have a common origin. Modelling this stochastic behaviour yields,
where e(t) and G(t) are Gaussian multiplicative and additive noises respectively with the following properties:
where the parameters D and a represent intensity of the multiplicative and additive noises respectively, and l denotes the strength of correlation between e(t) and G(t) with 0 < l < 1. Eq.(8) describes zero correlation time between additive and multiplicative noises. To study the effect of non-zero correlation time t, the statistical property áe(t)G(t')ñ is given by the more generalized form,
Eq.(9) reduces to Eq.(8) when t ® 0. According to the Langevin Eq.(4), we can derive the Fokker-Planck Equation for the positive values of x , which is given by,
where p(x,t) is the probability distribution function, A(x) and B(x) are the drift and diffusion coefficients respectively defined as
According to the reflecting boundary condition, the steady-state probability distribution function (SPDF) of Eq.(10) can be obtained from , . This is given by,
where N is the normalisation constant. Eq.(13) may be solved numerically to obtain the steady state probability distribution pst(x) and average cell population áxñ for the Gompertz and logistic models.
III. EFFECT OF ZERO CORRELATION TIME ON AVASCULAR TUMOR GROWTH LAWS
In this section we study the properties of tumor cell growth using the Gompertz law in the presence of cross correlated additive and multiplicative noises for the case of zero correlation time between noises. Although of course strictly speaking the correlation time is never actually equal to zero, we adopt this approach as a first approximation in considering the biological system driven by noise. This corresponds to the case where we deal only with correlated Gaussian white noise. This is valid biologically when the time scale for the correlation is much shorter than the time scale for the relaxation of the driven process.
In our tumor growth modelling, we assume that the stochastic characteristics arise from both external and internal factors. We further assume that external factors represent fluctuations due to some external factors such as treatments, i.e, chemotherapy or radiotherapy. These treatments affect the tumor's net growth rate generating multiplicative noise and at the same time they restrain the number of tumor cells which give rise to additive noise [9,10].
A. Gompertzian growth
The results of our calculations for the steady state probability pst(x) as a function of cell number x are presented in Figs. 1a and 1b. In particular, we have studied the effect of the intensity of correlation strength l between the additive and multiplicative noise intensities on pst(x). It is evident from the Fig. 1a that with the increase of l, pst(x) decreases at smaller x values ie, x < 5. However, the peak of pst(x) increases with higher values of l and at higher x values, pst(x) for different l, merge with each other and die out. This implies that higher values of l promote cell growth. In Fig. 1b we have studied the effect of negative correlation strength on pst(x). It is interesting to note that as the value of l decreases, pst(x) increases at small x, and decreases at large x. This implies that decreasing l causes the tumor cells to disappear. In other words the distribution of the cell population which was mainly peaked about zero (for smaller values of l) signifies high extinction rates, and moves towards zero with the decrease of the correlation parameter l. Similar behaviour has also been obtained in the case of logistic growth . Comparing Figs. 1a and 1b, we found that a change in sign concerning the correlation between additive and multiplicative noises has both positive and negative effects on tumor growth. Indeed, a negative value of l (which corresponds to a shift to lower values of the peak position of the pst(x)) has a meaning in this context. It corresponds to a change in sign of the correlation between additive and multiplicative noise, which corresponds to a better interpretation of the Eq. (4). In fact, in Eq. (4) the additive and multiplicative noise have opposite effects (when one is positive the other is negative). Therefore, it seems that a positive l indicates a negative feedback between the two noises. On the contrary a negative l, or alternatively a change in sign of either of the two noise terms in Eq.(4), may be correct, to simulate positive feedback of two effects in the tumor treatment. In Fig. 2, we present our results for average cell population áxñ as a function of the intensity of multiplicative noise D for different values of additive noise intensity a at a fixed l = 0.0. This is given by the equation,
We observe from Fig. 2 that áxñ increases first, and then decreases with the intensity of the multiplicative noise, showing a typical stochastic resonance characteristic. This means that the appropriate intensity of multiplicative noise is suitable for tumor cell growth and extra noise may restrain tumor cell growth. It is evident from Fig. 2 that as a increases the peaks die out and become flat. In Fig. 4, we also observe the same qualitative behaviour for the logistic model.
In Fig. 3a we present the results for pst(x) as a function of x at various values of a. This figure shows that the position of the peaks of the distribution remains unchanged as the intensity of the additive noise increases. However, the height of the peaks decrease as a increases. In Fig. 5a, we also observe the same qualitative behaviour for the logistic model.
We also show in Fig. 3b the average cell population áxñ as a function of a for different values of multiplicative noise intensity D. We observe that áxñ decreases and then increases showing a negative stochastic resonance-like characteristic at the values of D between 0.3 and 0.9. This is caused by the amplification of a weak periodic signal by means of noise. The resonances die out with increasing D. The critical phase transition occurs at the point of inflection, , from the deterministic Eq. (2) due to the sigmoidal nature of the Gompertz model. For comparison we show the results of the logistic model in Fig. 5(b). These results show a similiar trend to the Gompertz model. It may be noted that as the point of inflection, x = a/2b in the logistic model, (obtained from the deterministic Eq. (3)), is reached more slowly than in the Gompertz model, the resonances in Fig. 5b are more pronounced than in the Gompertz model in Fig. 3b. These resonances also die out with increasing D as in the Gompertz model. In both models, multiplicative noise induces a phase transition and resonance in tumor growth. This may require further study in order to understand the occurance of stochastic resonance-like phenomena and the underlying mechanism(s) in the present non-bistable systems.
IV. EFFECT OF FINITE CORRELATION TIME ON AVASCULAR TUMOR GROWTH MODELS
In this section we study the stationary properties and average populations of tumor cell growth using the Gompertz and logistic laws in the presence of cross correlated additive and multiplicative noises for the case of nonzero correlated time between noises. This reflects the inclusion of coloured noise into our system and is more appropriate in representing a non-sharp separation of time scales and a greater precision in analysis of the stochastic process. This generalization relaxes the restriction that the correlation time is zero which is the case when considering Gaussian white noise as in Sec. III above. In modelling the Gompertz and logistic growth laws driven by cross-correlated noises for the case of nonzero correlation time, we implement Eq. (9) instead of Eq. (8) in our probability distribution calculation (Eq. 13).
A. Gompertzian growth
First we study the effect of l on pst(x) as function of x for a fixed value of correlation time t = 0.3. This is presented in Fig.6(a). This figure shows that pst(x) as a function of x decreases at smaller values of x i.e, below x around 7.5 and then increases as we increase l. For the higher values of x the curves of different l merge with each other and die out. On comparing these results with the results presented in Fig. 1a for Gompertz growth model in the case of zero correlation time, we observe that, as expected there is no change in qualitative nature of the curves and these results follow the similar trend as the previous ones. However, here we are interested in studying how both the correlation strength l and the correlation time t affect the probability distribution. In Fig. 6b, we present the results of pst(x) as a function of x for a fixed value of correlation strength l = 0.2 and varying the correlation time t. As shown in the figure the probability as a function of x decreases as we increase t. On comparing Figs. 6a and 6b, we find that the correlation strength l and the correlation time t play the opposite role on the steady state probability i.e., pst(x) increases with the increase of l and decreases with the increase of t.
To see the effect of l and t on the cell population, we examine in Figs. 7a and 7b, the average cell population áxñ, as a function of t for various values of l and as a function of l for different values of t, respectively. From these figures, we observe that áxñ as a function of t increases with increasing l and, viewed as a function of l, decreases with increasing values of t. This also confirms that the correlation strength and the correlation time have the opposite effect on the average cell population in the case of Gompertz growth.
In Fig. 8 we examine the effect of t on the average cell population áxñ as a function of additive noise intensity a for different values of multiplicative noise intensity D at a fixed l = 0.1. áxñ as a function of a decreases for some fixed values of D and then increases showing a negative resonance-like characteristic and these resonance-like characteristics die out as D is increased. However, there is no change in qualitative nature of the curve with the change in t, i.e., the correlation time has no effect on the resonance-like characteristic (see left and right figures, these are for different t = 0.2 and 0.6). Similiar results have been obtained for the logistic models.
B. Logistic growth
We performed similar calculations for the nonzero correlation time for the logistic growth. As expected, since, logistic growth law predicts similiar behaviour to the Gompertz growth law, we obtained qualitatively similar behaviour of the curves in case of the logistic growth laws (see Figs. 7a and 7b). Similar calculations  have also been carried out in the case of logistic growth. However, they have shown that the correlation strength and the non-zero correlation time have the opposite effect on pst(x) contrary to our findings in the present calculations. The discrepancy between our results and those of  in this section are reflected in the fact that  have based their model on the original model of .
We also studied the effect of nonzero correlation time on average cell population (Figs. 10a and 10b). We observe that the correlation strength and the correlation time play opposite roles, which is similar to the case of the Gompertz growth law, in Sec. IV.A. We also find the average cell population increases with increasing l and decreases with increasing t. These results are also in qualitative disagreement with the calculations of  for reasons already outlined above.
As in the Gompertz model, we study the effect of t at a fixed l = 0.1 on the average cell population áxñ as a functon of additive noise intensity a for different values of multiplicative noise intensity D in the case of logistic growth model and the results are shown in Fig.11. As one can see these results are analogous to the results obtained in the case of zero correlation time in the logistic model and show similar resonance-like characteristics. However, as in the case of Gompertz growth law here we also observe that there is no effect of t on the resonance.
We have considered the steady state properties of tumor cell growth, average cell populations and the effect of correlated noise using the Gompertz model. Our results were compared with existing calculations of the logistic model [9,19] and . In each of the Gompertz and logistic models we found that the both negative and positive correlation between additive and multiplicative noise were important in predicting the likely outcome of treatment protocols. From our results it is also observed that an inappropriate noise intensity can lead to the malignant growth of tumor cells. Another interesting result we found, is that multiplicative noise induces a phase transition and resonance in tumor growth. This behaviour arises due to the sigmoidal nature of both the Gompertz and logistic models which both have points of inflection that are always a fixed proportion of their asymptotic values. It would be interesting to pursue this study further, in order to understand the underlying mechanism(s) for the negative stochastic resonance like characteristic in the Gompertz and logistic non-bistable models. We have studied the stationary properties of the Gompertz and logistic tumor cell growth models in the case of nonzero correlation time between additive and multiplicative noises. We found that the correlation time and the correlation strength play opposite roles on the steady state probability distribution in case of both the Gompertz and logistic model. As we increased the correlation strength, the steady state probability and the average cell population increased. This scenario is reversed if the correlation time is increased. Similarly to the case of nonzero correlation time, we observed a phase transition and resonance like characteristics in each of the Gompertz and logistic models. These phenomena are unaffected with the change of correlation time. It is worth mentioning here that in the logistic model we found our calculations significantly differed from those of both  and . These results can be reproduced if one considers a negative correlation strength (l) between additive and multiplicative noise [19, 20]. Finally there were no significant differences in the interpretation of the results between our simulations for the Gompertz and logistic models. This is probably due to the fact that these growth models belong to the same class of sigmoid function. These results thus provide strong evidence of a universal behaviour of sigmoid laws which are present in several problems of biological growth.
We thank Dr Sahoo for useful discussions. One of us A. Behera acknowledges financial assistance for an APG Studentship awarded by Queen's University, Belfast. SFC O'Rourke acknowledges current support from the Leverhulme Trust (Grant No. F/00203/K).
 L. Gammaitoni, P. Hanggi, P. Jung, and F. Marchesoni, Rev. Mod. Phys. 70, 223 (1998). [ Links ]
 Ya. Jia and Jia-rong Li, Phys. Rev. Lett. 78, 994 (1997). [ Links ]
 A. A. Zaikin, J. Kurths, and L. Schimansky-Geie, Phys. Rev. Lett. 85, 227 (2000). [ Links ]
 J. D. Murray, Mathematical Biology I. Springer: Berlin. (2002). [ Links ]
 M. Molski and J. Konarski, Phys. Rev. E. 68, 021916 (2003). [ Links ]
 S. C. Ferreira Jr., M. L. Matrins, and M. J. Vilela, Phy. Rev. E 67, 051914 (2003). [ Links ]
 A. Bru, S. Albertos, J. A. Lopez, I. Asenjo-Garcia, and Bru, Phys. Rev. Lett. 92, 238101 (2004). [ Links ]
 G. Albano and V. Giorno, J. Theo. Biol. 242, 329-336 (2006). [ Links ]
 B. Q. Ai, X. J. Wang, G. T. Liu, and L. G. Liu, Phys. Rev. E. 67, 022903 (2003). [ Links ]
 D. C. Mei, C. W. Xei, and L. Zhang, Eur. Phys. J. B. 41, 107 (2004). [ Links ]
 M. Marusic, Z. Bajzer, J. P. Freyer, and S. Vuk-Pavlovic, Cell. Prolif.27, 73 (1994). [ Links ]
 G. G. Steel, Basic Clinical Radiobiology. Arnold: third edition. (2002). [ Links ]
 M. A. A. Castro, F. Klamt, V. A.Grieneisen, I. Grivicih, and J. C. F. Moreira, Cell. Prolif. 36, 65 (2003). [ Links ]
 M. Gyllenberg and G. F. Webb, Growth, Dev. and Aging. 53, 25-33 (1989). [ Links ]
 A. R. Kansal, S. Torquato, E. A. Chiocca, and T. S. Deisboeck, J. Theor. Biol. 207, 431 (2000). [ Links ]
 F. Kozusko and Z. Bajzer, Math. Biosci. 185, 153 (2003). [ Links ]
 H. Risken, The Fokker-Planck Equation (Springer: Berlin), 96 -99 (1996). [ Links ]
 C. W. Gardiner, Handbook of Stochastic Methods. (Spinger: Berlin), 108 - 125 (2003). [ Links ]
 A. Behera and S. F. C. O'Rourke, Phys. Rev. E. 77, 013901 (2008). [ Links ]
 B. Q. Ai, X. J. Wang, and Li. Q. Liu, Phys. Rev. E. 77, 013902 (2008). [ Links ]
Received on 3 April, 2008