PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of tbiomedBioMed CentralBiomed Central Web Sitesearchsubmit a manuscriptregisterthis articleTheoretical Biology & Medical ModellingJournal Front Page
 
Theor Biol Med Model. 2010; 7: 30.
Published online 2010 July 30. doi:  10.1186/1742-4682-7-30
PMCID: PMC2920266

Tracking tumor evolution via prostate-specific antigen: an individual post-operative study

Abstract

Background

The progress of the prostate-specific antigen (PSA) level after radical prostatectomy is observed for a patient in order to extract information about the mode of tumor cell growth. Although PSA values are determined routinely to find the doubling time of the prostate marker, to our knowledge, this analysis is the first in the literature.

Results

The prostate tumor marker values were determined regularly after the surgery and plotted on a logarithmic scale against time. An initial rapid-growth mode changed to a slower power-law regime within two years of surgery. Our analysis associates this observation with a transition in the growth mode from unrestricted growth of dispersed cells to their clumping into macroscopic structures.

Conclusions

Such studies may help determine the appropriate time window for postoperative therapies in order to increase the life expectancy of the patient.

Background

Cancer of the prostate gland is one of the most frequently diagnosed male illnesses and may lead to death of the patient. The carcinoma is routinely detected by a straight-forward blood test that measures a glycoprotein called prostate-specific antigen (PSA). At an early stage of cancer growth with a localized tumor, radical removal of the prostate gland has proved to be the optimum treatment. If the PSA value rises after radical prostatectomy, different alternatives for treatment are currently under debate. The doubling time (DT) of the PSA value is accepted as a strong prognostic factor for the risk of cancer death. In a group of 379 patients, almost no prostate cancer deaths were recorded within approximately 4 years of prostate removal for 3 < DT < 8.9 months, while some patients with DT < 3 months died within 1.5 years [1]. It is therefore reasonable to infer that findings during the last few years based on long-term statistics suggest a longer life expectancy for patients with postoperative radiotherapy that follows (within 6 months of) radical prostate surgery [2]. Transdermal radiotherapy is commonly applied after the PSA level reaches a threshold value. However, a wait-and-watch method may cost valuable time and the relevant moment for action may be missed regardless of how low the threshold value is set. Many authors have already suggested [3-5] that the entire course of tumor growth offers important information regarding the clinical strategies to be followed.

The purpose of this work is threefold. (a) To indicate the possibility of detecting fast (exponential) growth of the PSA score well before an arbitrary threshold value is reached, thus gaining time for deciding the therapies to be followed. In principle, this strategy is analogous to determining the DT. (b) To analyze PSA data in a way that reveals a sharp crossover from exponential to power-law growth. (c) To propose a simple model to explain the crossover to slower (power-law) growth.

The second, slower growth regime, we argue, arises because of the coalescence or "condensation" [6] of freely-dividing cancerous cells to form one or more compact tumors, with growth essentially confined to the edges or the surface [4,5,7]. It has been pointed out that at this stage (i) "sensitivity to anti-metabolic drugs decreases, ... (since the fraction of) tumor cells that are in the cell division cycle decreases" [3] and (ii) angiogenesis is expected to start [5].

Results

Case presentation

After a radical prostatovesiculectomy (pT2c N0 M0 G2, Gleason 3+4 = 7) applied to one of us (ME), PSA values were determined with state-of-the-art precision using constant laboratory conditions (Viollier, Brunngasse 6, CH-8401 Winterthur) at time intervals of initially 6 months (see Table Table1).1). The error in time measure was ± 0.5, while each PSA value was determined with an uncertainty of ± 0.002 μg/l. In Figure Figure11 we plot the values listed in the table as a function of time t in months after the surgery. The graph has a characteristic "U" shape, i.e., a shallow increase during the first two years and a steep rise in the last two.

Table 1
Measured PSA scores
Figure 1
Linear plot of PSA score vs. time. PSA values in μg/l are displayed as a function of measurement time in months after the surgery.

In order to determine the kinetics of cell growth underlying the PSA progression, we present the same PSA values as a function of time in Figure Figure2a,2a, but plotted on a semi-logarithmic scale. From the time of surgery until about 30 months thereafter we observe a clear linear increase. Analytically, this corresponds to "exponential growth", with the functional form ~exp(pt), where p is a constant rate of cell division. This can be determined at an early stage, before an arbitrary threshold value is attained.

Figure 2
Logarithmic plot of PSA score vs. time. PSA values plotted on a logarithmic scale against (a) linear time and (b) time on a logarithmic scale. Note that the straight line fit in panel (a) signals exponential growth, while that in panel (b) indicates "power-law" ...

Later values are noisier and remain below the straight line, implying a different, slower growth law. In Figure Figure2b2b we examine the departure from exponential behavior by displaying the PSA values vs. time on a log-log plot. Here a straight line signifies a power law (with the functional form ~tu). There is sharp crossover from exponential to power-law behavior at about two years after the surgery, rather than a gradual slowing down. From this point on, up to the last measured value, the PSA values grow as a power of the time elapsed after the crossover point.

Growth modes

We assume that the PSA value is linearly proportional to the number N of carcinoma cells and that, initially, each cell freely divides at a constant rate p (probability per unit time). For such unrestricted growth, the number of cells, N, increases by pN per unit time on average, i.e.,

equation image
(1)

At any time t from the start of the growth process, N is found to be

equation image
(2)

where N0 is the number at t = 0. The growth rate p is the only relevant parameter that has to be experimentally determined. In this type of growth, the birth of a new cell has no effect on that of subsequent cells.

As the number of malignant cells grows within the tissue, there must be a spontaneous formation of macroscopic clusters of cells (i.e., tumors) that mop up almost all the microscopic clusters, at a rather sharp transition point, the so-called percolation threshold [8]. We identify this threshold with the crossover observed in Figure Figure2b.2b. Once clusters are formed, growth is confined to the surface of the tumor [5] and the kinetic equations should only involve the number Ns of actively dividing cells in the surface layer.

The number of cells in the tumor is roughly N ~ RD where R is the average radius and D is the dimension of the cluster. The number in the surface layer will grow as Ns ~ RS, with S being the surface dimension [9]. In the most general case where D ≤3, Equation 1 is replaced by

equation image
(3)

where k1 and k2 are constant geometrical factors, leading to

equation image
(4)

Here An external file that holds a picture, illustration, etc.
Object name is 1742-4682-7-30-i5.gif, with N0 having the same meaning as in Equation 1, and c2 = pk2. Irrespective of the exact value of the exponent, for a sufficiently small c1, Equation 3 predicts power-law growth. For spheroidal tumors [7] with growth confined to the surface region, u = 3. The parameter u is related to the compactness or looseness of the clumps of malignant cells and has to be determined experimentally.

Data analysis

Figure Figure2a2a confirms exponential growth of the initial PSA scores with a rate of p = 0.090 ± 0.004 per month. Thus, the DT is about 8 months [DT in months is given by (ln 2)/p], or the PSA score increases by more than a factor of three within a year. In current practice, each datum point at a particular time is used by the physician to assess the patient's health condition and to decide upon further action. Although the absolute PSA values are much lower than the widely accepted threshold values for the recurrence of prostate cancer, the growth rate is alarmingly fast. Yet after about three years, the growth rate has slowed down and is seen to deviate from exponential behavior. Nevertheless, any arbitrarily set threshold will eventually be attained.

The log-log plot of the PSA values vs. time in Figure Figure2b2b shows a crossover from an exponential to a power law. The value of the exponent is u = 2.57 ± 0.07, very close to what would be predicted for percolation clusters [8,9] with growth confined to an outer shell.

The values for the coefficient p in the exponential form exp(pt) and the exponent u in the form tu are found from least squares fits to the linear parts of the plot in the respective cases. It should be emphesized that the main import of the paper is not the precise numerical values of p or u, although with a linear fit to the straight lines in the two plots these can be determined to an accuracy of two (one) significant digits (digit), respectively, with Pearson's correlation coefficients of r2 = 0.98 and 0.96. The point is that there is an unambiguous crossover, from a characteristically exponential to a characteristically power-law behavior.

Discussion

The "Gompertzian" [10] growth curve widely accepted as a "universal law" [3-7,11-14] for the growth of diverse populations including tumors and cell cultures [13] exhibits initial exponential growth, gradually slowing down and finally saturating (in vitro) to a constant value. We find that the data reported in Table Table11 are also fitted reasonably well by a Gompertzian (see Figure three of Ref [14]). It would be worthwhile to re-plot the data traditionally analyzed within this gradualist picture on a log-log scale, and see whether the same sharp crossover behavior is actually hiding there as well.

Once the switch to power-law growth occurs, signalling discrete, compact tumors, we can estimate their size assuming that the PSA value is linearly proportional to the total number of malignant cells and knowing the PSA value and tumor size obtained from magnetic-resonance imaging (MRI) at the time of the operation.

In the particular case under study, the gland, prior to its removal, had a diameter of less than 40 mm (Huber D: MRI Report. Zurich: Klinik Hirslanden; 2002) and the PSA value was 8.5 μg/l. The prostate was about 50% cancerous according to the post-operative biopsy, so we deduce that a PSA level of 0.485 μg/l, measured at t = 65, would correspond to a tumor (consisting only of cancerous cells) approximately 5 mm in diameter. If not one but two equally-sized compact clusters condense out of the scatter of individual cells, our naive calculation gives a diameter of approximately 4 mm for each of those tumors. Digital rectal examinations by three independent experts (Brodmann S; Riesterer O; Vollenweider P; 2008) as well as the MRI analysis (Hilfiker P: Medical Report. Zurich: MRI Bethanien; 2008) performed at t = 67 revealed two masses of about 4 mm in agreement with our prediction. Subsequently, the subject received 70 Gy of radiation therapy in the anastomosing region during weekdays of t = 68 and 69 in equal doses, with full bladder and an inflated (50 ml) rectal balloon. The volume was reduced after 46 Gy in order to spare the vicinity. Since the radiation therapy, the PSA values have remained around 0.16 μg/l.

Conclusions

We find that the PSA values of this patient after surgery follow an initial exponential growth curve, with a sharp crossover to a slower power-law regime at around 20 months. We argue that this may be due to the clumping of individual cancerous cells into macroscopic structures with distinct surfaces, to which subsequent growth is confined.

In the present case, after the onset of the power-law growth regime, it was possible to detect the clusters of cells, i.e., the tumors, via standard imaging techniques, and verify that their sizes coincided with predictions from the model.

Radiation therapy is not routinely applied after the surgery. The predictive power of our simple analysis, however, makes it highly worthwhile to monitor the PSA scores closely during the so-called wait-and-watch period no matter how low their absolute value is, as has been done in the present case, in order not to miss the optimum time window for post-operative therapy to increase the life expectancy of the cancer patient.

The doubling time is a widely-used characteristic of cell growth and a constant value testifies to exponential growth. The novel aspect of our analysis deals with the deviation from exponential behavior, which transforms into a subsequent power-law regime. Though mathematically convincing, this non-exponential late-stage growth behavior observed for one particular patient may not represent a universal phenomenon. However, monitoring the growth of the PSA value as a function of time may provide the opportunity for observing such a crossover, as for this patient. If such a crossover indeed occurs, it may indicate, as explained above, the formation of macroscopic masses within the tissue, and we suggest that this may be taken as a clue for deciding upon post-operative treatment.

Consent

Written informed consent for publication was obtained from the patient, who is one of the authors (ME), for publication of this case report and accompanying images. A copy of the written consent is available for review by the Editor-in-Chief of this journal.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

The authors contributed equally to this work, and read and approved the final manuscript.

Acknowledgements

ME thanks Drs. M Dubs and P Vollenweider for fruitful discussions during the active surveillance throughout the years before and after prostate surgery. AE would like to thank M Kardar for a helpful comment and also acknowledge partial support by the Turkish Academy of Sciences.

References

  • Freedland SJ, Humphreys EB, Mangold LA, Eisenberger M, Dorey FJ, Walsh PC, Partin AW. Death in patients with recurrent prostate cancer after radical prostatectomy: Prostate-specific antigen doubling time subgroups and their associated contributions to all-cause mortality. J Clin Oncol. 2007;25:1765–71. doi: 10.1200/JCO.2006.08.0572. [PubMed] [Cross Ref]
  • Bolla M, van Poppel H, Collette L, van Cangh P, Vekemans K, Da Pozzo L, de Reijke TM, Verbaeys A, Bosset J-F, van Velthoven R, Maréchal J-M, Scalliet P, Haustermans K, Piérart M. Postoperative radiotherapy after radical prostatectomy: a randomised controlled trial. Lancet. 2005;366:572–578. doi: 10.1016/S0140-6736(05)67101-2. [PubMed] [Cross Ref]
  • Schabel FM Jr. The use of tumor growth kinetics in planning "curative" chemotherapy of advanced solid tumors. Cancer Res. 1969;29:2384–2389. [PubMed]
  • Brú A, Albertos S, Subiza JL, Garcia-Asenjo JL, Brú I. The universal dynamics of tumor growth. Biophys J. 2003;85:2948–2961. doi: 10.1016/S0006-3495(03)74715-8. [PubMed] [Cross Ref]
  • Kohandel M, Kardar M, Milosevic M, Sivaloganathan S. Dynamics of tumor growth and combination of anti-angiogenic and cytotoxic therapies. Phys Med Biol. 2007;52:3665–3677. doi: 10.1088/0031-9155/52/13/001. [PubMed] [Cross Ref]
  • Torkington P. Kinetics of deterrence of Gompertzian growth. Bull Math Biol. 1983;45:21–31. [PubMed]
  • Delsanto PP, Guiot C, Degiorgis PG, Condat CA, Mansury Y, Deisboek TS. Growth model for multicellular tumor spheroids. App Phys Lett. 2004;85:4225–4227. doi: 10.1063/1.1812842. [Cross Ref]
  • Stauffer D, Aharony A. Introduction to Percolation Theory. London: Taylor and Francis; 1992.
  • Mandelbrot BB. The Fractal Geometry of Nature. New York: Macmillan; 1983.
  • Gompertz B. On the nature of the function expressive of the law of human mortality, and on a new mode of determining the value of life contingencies. Phil Trans R Soc. 1825;115:513–583. doi: 10.1098/rstl.1825.0026. [PMC free article] [PubMed] [Cross Ref]
  • DeWys WD. Studies correlating the growth rate of a tumor and its metastases and providing evidence for tumor-related systemic growth-related factors. Cancer Res. 1972;32:374–379. [PubMed]
  • Iwata K, Kawasaki K, Shigesada N. A dynamical model for the growth and size distribution of multiple metastatic tumors. J Theor Biol. 2000;203:177–186. doi: 10.1006/jtbi.2000.1075. [PubMed] [Cross Ref]
  • Norton L, Simon R, Brereton HD, Bogden AE. Predicting the course of Gompertzian growth. Nature. 1976;264:542–545. doi: 10.1038/264542a0. [PubMed] [Cross Ref]
  • Guiot C, Degiorgis PG, Delsanto PP, Gabriele P, Deisboek TS. Does tumor growth follow a "universal law"? J Theor Biol. 2005;225:147–151. doi: 10.1016/S0022-5193(03)00221-2. [PubMed] [Cross Ref]

Articles from Theoretical Biology & Medical Modelling are provided here courtesy of BioMed Central