First Published Online October 6, 2008 The Oncologist, Vol. 13, No. 10, 1046-1054, October 2008; doi:10.1634/theoncologist.2008-0075 © 2008 AlphaMed Press
Tumor Growth Rates Derived from Data for Patients in a Clinical Trial Correlate Strongly with Patient Survival: A Novel Strategy for Evaluation of Clinical Trial DataaMedical Oncology Branch, Center for Cancer Research, National Cancer Institute, National Institutes of Health, Bethesda, Maryland, USA; bDepartment of Biological Chemistry, Silberman Institute of Life Sciences, Hebrew University, Jerusalem, Israel; cHubert Department of Global Health, Rollins School of Public Health, Emory University, Atlanta, Georgia, USA; dHebrew University-Hadassah School of Public Health, Hebrew University, Ein Kerem Medical Centre, Jerusalem, Israel Correspondence: Tito Fojo, M.D., Ph.D., Medical Oncology Branch, National Cancer Institute, Center for Cancer Research, National Institutes of Health, Building 10, Room 12N226, 9000 Rockville Pike, Bethesda, Maryland 20892, USA. Telephone: 301-402-1357; Fax: 301-402-1608; e-mail: tfojo{at}helix.nih.gov Received March 27, 2008; accepted for publication August 13, 2008; first published online in THE ONCOLOGIST Express on October 6, 2008. Disclosure: The content of this article has been reviewed by independent peer reviewers to ensure that it is balanced, objective, and free from commercial bias. No financial relationships relevant to the content of this article have been disclosed by the authors, planners, or staff managers.
Purpose. The slow progress in developing new cancer therapies can be attributed in part to the long time spent in clinical development. To hasten development, new paradigms especially applicable to patients with metastatic disease are needed. Patients and Methods. We present a new method to predict survival using tumor measurement data gathered while a patient with cancer is receiving therapy in a clinical trial. We developed a two-phase equation to estimate the concomitant rates of tumor regression (regression rate constant d) and tumor growth (growth rate constant g). Results. We evaluated the model against serial levels of prostate-specific antigen (PSA) in 112 patients undergoing treatment for prostate cancer. Survival was strongly correlated with the log of the growth rate constant, log(g) (Pearson r = –0.72) but not with the log of the regression rate constants, log(d) (r = –0.218). Values of log(g) exhibited a bimodal distribution. Patients with log(g) values above the median had a mortality hazard of 5.14 (95% confidence interval, 3.10–8.52) when compared with those with log(g) values below the median. Mathematically, the minimum PSA value (nadir) and the time to this minimum are determined by the kinetic parameters d and g, and can be viewed as surrogates. Conclusions. This mathematical model has applications to many tumor types and may aid in evaluating patient outcomes. Modeling tumor progression using data gathered while patients are on study, may help evaluate the ability of therapies to prolong survival and assist in drug development.
The therapy of cancer continues to challenge oncologists. The pace of progress has often been slow, in part because of the time required to evaluate new therapies. With a few notable exceptions, such as the proteasome inhibitor bortezomib, the time from conception to approval for most new cancer therapeutics is at least 10 years, with a large part of this time spent in clinical development [1]. To reduce the time to approval, new paradigms to assess therapeutic efficacy are needed. Furthermore, because most new therapies are initially tried and approved in patients with metastatic disease, such novel paradigms must be especially applicable to this patient population. For patients with advanced prostate cancer, androgen deprivation remains the first therapeutic option, and usually results in a prompt and often sustained clinical response [2]. However, the disease can eventually evolve to androgen-independent prostate cancer (AIPC), and it is at this stage that current therapeutic efforts are increasingly focused [3–5]. Metastatic prostate cancer, and especially AIPC, behaves like other drug-resistant cancers and together with the tumor marker PSA is an excellent model for metastatic cancer, and for evaluating new strategies for disease assessment. The present study was undertaken to develop a rational basis for evaluating drug efficacy in patients participating in clinical trials using data gathered while the patient is enrolled in the study. We describe a novel paradigm for predicting drug efficacy using prostate cancer and PSA values as a model that may also be applicable to most cancers where tumor load can be assessed by serum or radiographic measurements.
Patient Characteristics The data for this analysis came from two clinical trials approved by the institutional review board of the National Cancer Institute. The primary objective of both trials was to determine whether novel combinations of chemotherapy produced sufficiently high clinical responses to warrant further investigation in patients with AIPC [6, 7]. All patients had metastatic AIPC and had failed to benefit from combined androgen blockade, as well as antiandrogen withdrawal. Therapies consisted of either thalidomide plus docetaxel (docetaxel as the control arm) or ketoconazole plus hydrocortisone and alendronate (ketoconazole plus hydrocortisone as the control arm). The overall survival time was calculated from the on-study date until the date of death. In total, 112 patients were included in the mathematical analysis, comprising all patients enrolled in the two studies. Thirty-two patients who enrolled in both studies are only included in the analysis once with the start date as the date of enrollment in their first protocol. All datapoints were obtained from patients during enrollment in clinical trials. None of the datapoints were obtained after patients were removed from clinical trials.
Mathematical Analysis The Regression–Growth Equation
We note that, for individuals in whom the data show a continuous decrease from the time of treatment, equation (1) can be replaced by the reduced form
is the time until regrowth reached the pretreatment value.
Relation Between g, d, and the Nadir and Time to Nadir
Differentiating equation (1) with respect to t, and setting the value of df/dt equal to zero, gives the time taken to reach the minimum as:
Substituting tmin for t in equation (1), one obtains
Data Analysis
Statistical Analyses
The median number of datapoints was 10 per patient and the median time over which data were collected was 227 days. Following the curve-fitting analysis of each patient dataset, we identified five groups among the 109 patients for which we had more than two time points as shown in Figure 2 and online supplementary Figure S1. The growth rate constants varied over a nearly 1,500-fold range, while the regression rate constants varied over a 50-fold range (Fig. 3A). Furthermore, the regression rate constants were consistently larger than the growth rate constants, with median values of 10–1.7 day–1 versus 10–2.5 day–1, respectively.
Correlations with Survival Time Figure 3B depicts the distributions of patient survival times and the two rate constants for the 99 patients in whom the PSA time courses had a g or d (or both) with an associated p < .05. Survival was more strongly correlated with the logarithm of the growth rate constant than with the logarithm of the regression rate constant. This result suggests that, while a given therapy may result in tumor reduction, the critical determinant in survival is whether or not the therapy alters the inherent growth rate of the tumor. Figure 3C depicts the correlations between patient survival times and the initial PSA level, the minimum PSA level, and the time to the PSA minimum. Figure 4 depicts Kaplan–Meier plots of fractional survival against time of survival for the upper and lower 50% of cases, in each case stratified by log(g), min, tmin, and the initial PSA value. Patients whose tumor growth rate constants were in the upper 50% had shorter survival times than patients with tumor growth rate constants in the lower 50%. The time to the minimum PSA level also had a strong impact on survival, comparing the upper 50% with the lower 50%, as did the nadir or minimum. An initial PSA signal in the upper 50% had a small, but statistically significant, impact on survival. These results are best understood when one recognizes the dependence of time to minimum and minimum on the growth rate constant (equation (5) and equation (6A), respectively), so that these values are surrogates of the growth rate constant. This is underscored by the results in online supplementary Figure S2.
Patient-to-Patient Variability in the Kinetic Parameters Figure 5A depicts the distribution of the growth rate constants (g) in the form of a histogram. Interestingly, when the clinical outcome was evaluated, we found that, among the 46 patients with the highest growth rate constants, 28 experienced progressive disease (PD), only four had a partial regression (PR), one had a PSA response, ten were scored as stable disease (SD), and three were not evaluable. In contrast, among the 46 patients with the lowest growth rate constants, there were only two patients with PD, 26 patients who experienced a PR, one who had a minimal response, an additional 11 with a PSA response, five demonstrating SD, and one that was not evaluable. The shaded columns depict the validation dataset (see below).
Validation Study The approach used in the present study was validated using data from an independent trial conducted by the Southwest Oncology Group, as depicted in Figure 5. The shorter lower line depicts the regression through the data obtained from the validation set. Neither the slope of this line nor its intercept is statistically different from the line above it representing the study group, supporting the validity of our analysis.
Clinicians and regulators have an interest in early prediction of treatment outcomes. Regulatory agencies often seek improved survival in the approval of a new cancer therapeutic. However, the effect of a therapeutic agent on survival has become much more difficult to define, especially as patients increasingly enroll in successive treatment regimens. The current study underscores the obvious: the growth of treatment-refractory cancer cells is responsible for the death of a patient. We show that survival predictions can be made in patients with metastatic disease using data gathered while a patient is enrolled in a clinical trial and undergoing treatment. While in the present study we have used prostate cancer as a model, the biology described is likely applicable to many cancers wherein measures of tumor load, including values such as radiographic measurements, are available. It has long been recognized that, at least during part of their growth, tumors follow exponential kinetics. Gompertzian tumor growth has been debated, and some have argued that it cannot apply to tumors over their entire life span [10, 11]. In our dataset (see online supplementary Fig. S1), only three cases of 112 (11, 47, and 77) show even a suggestion of the limited terminal plateau that the Gompertzian equation describes. While it may still be true that a Gompertzian equation might describe the entire growth history of a tumor, our results clearly demonstrate that, during the window of observation that our datasets comprise, the tumors that grow appear to be growing exponentially. One might consider the total course of growth of a tumor as consisting of a series of intervals each with a different growth rate constant. After all, why should the growth rate constant be fixed over the years during which a tumor grows? Patients with advanced disease may be in the final interval, such that during the brief period of a clinical trial we are able to measure a growth rate constant that in turn predicts survival. Furthermore, most models consider tumor growth in the absence of treatment. The mathematical equation used in this work recognizes that, especially during therapy, both tumor growth and regression occur simultaneously, and discerns their independent contributions to the measured growth. To our surprise, however, the regression portion of the curve, while needed to accurately describe the data, does not predict survival in these patients with prostate cancer. It is the growing (surviving) fraction that determines survival. The use of mathematics to describe tumor kinetics has been widely explored in prostate cancer because of the sensitivity and specificity of the tumor marker PSA. Two derivations, PSA doubling time (PSA-DT) and PSA velocity (PSAV), have received special attention [12–22]. The PSA-DT rests on the assumption that increases in PSA follow first-order kinetics and hence an exponential growth curve, so that a plot of the log of the PSA versus time produces a slope that should remain constant provided the patient is not receiving an effective therapy. PSA-DT has been advanced as a method to discern disease aggressiveness [13, 16, 17, 19]. Others have described its usefulness in predicting cancer-specific mortality [21]. Similar results have been reported for PSAV, the change in PSA over time [12, 20]. Like the growth rate constant described in this study, PSA-DT is a mathematical estimate of the rate of tumor growth. With three or more PSA measurements, the PSA slope, PSA-DT, is calculated using a least-squares regression formula and the natural log of PSA values in order to make the kinetic pattern more linear. But in most cases, even this natural log transformation does not conform to a purely linear kinetic pattern. We would argue that one explanation for the poor fit is the failure to model the two concurrent processes of regression and growth. The mathematical calculations used here discern these two independent variables. Our observations, correlating overall survival with a calculated growth rate constant, have precedence in the literature both in patients who have received definitive local therapy and are followed without treatment and in patients with metastatic disease receiving generally ineffective therapies. In both of these scenarios, tumor growth unopposed by any regression results in a PSA-DT that is comparable with the growth rate constant g described in the current study [21–24]. Precedence can also be found in patients with metastatic hormone-refractory prostate cancer, in whom PSAV has been shown to be associated with the time to death after treatment [25]. However, with the advent of more effective therapies, models that account for simultaneous tumor regression and growth are needed, because net growth is a complex interplay of cellular proliferation and necrosis/apoptosis/senescence. The small sample size and patient enrolment in clinical trials at a single institution are limitations of the present study. Positive attributes include the fact that all but six of the 1,369 available PSA values were used in the calculations, avoiding arbitrary selection of PSA measurements within an interval of time. Another positive attribute is that the growth rate constant is a continuum and survival correlates very highly across this continuum. We would also note that because patients with metastatic disease usually die from their cancer, the analysis in this study was likely not biased by death from causes other than cancer. We also recognize that the value of this approach will need larger trials for validation. For a therapy to alter survival it must fully eradicate the tumor, reset the pace by selecting for a clone with a slower growth rate constant, or alter the pretherapy growth rate constant by affecting the biology of all cells without killing them. Because the Response Evaluation Criteria in Solid Tumors define disease progression as an increase in a single dimension of 20% (a volumetric increase of 72.8% if one assumes a sphere), the tumor volume can nearly double before disease is scored as progressive. As depicted in online supplementary Figure S3, if our dataset was limited to PSA values that had risen no more than 72.8% above the starting value or the nadir, this limited dataset could also reliably predict survival. We note, however, that we cannot exclude in every case that therapy might alter the growth rate without "resetting it," and that the measured change in the growth rate might require continued therapy. We would also caution that the conclusions and their potential relevance to other cancers apply only to patients with metastatic disease who are treatable and for whom therapeutic interventions have limited efficacy. Finally we recognize that the therapies used in these patients included conventional cytotoxic agents. Similar studies will have to be conducted with novel targeted agents to determine whether the same principles apply. Examining the derived growth and regression rate constants leads one to several observations that match clinical intuition. First, the self-evident observation that patient survival relates inversely to the tumor growth rate. Second, the relationship between the minimal value of PSA and the growth rate constant implies that the rate of tumor growth determines the quality of the response to therapy. The minimum is lower if the growth rate constant is lower, in the face of at least modestly effective therapy. Nearly all of the partial and complete responses—as judged by the depth of tumor reduction—occurred in patients with slower growth rates. Another observation relative to the clinical heuristic is that the regression rate is faster than the growth rate. We can agree that, generally, tumors respond faster than they grow over time. Further, the regression rate has a much narrower range across the patient population than the growth rate. This implies that the biology of tumor regression may be more similar across patients than that of tumor growth—there are very few pathways for cell death, but a myriad of growth factors and signal transduction pathways to facilitate/encourage growth. A final implication of these studies is that a higher risk patient population with a more rapid growth rate and a higher hazard of earlier death (hazard ratio, 5.14) can be identified early on in the course of treatment. The analysis presented here underscores the obvious: a successful therapy must either eradicate the tumor in its entirety or significantly alter the growth rate constant. But the analysis here also demonstrates clearly that, given a reliable measure of tumor load, the growth rate constant can be estimated from data gathered while a patient is on study. We believe that this approach will not be limited to analyses where serum markers are measured. In an accompanying manuscript [26], we have applied the same methodology to analyze data obtained by radiographic measurement of renal cell carcinoma and have obtained similar results. For a homogenous patient population, it ought to be possible to generate a "mean" growth rate constant for a given disease after, for example, failure of second-line therapy. This in turn should allow more rapid discernment of those therapies that can prolong survival. Because the growth rate constant is a validated surrogate for survival, investigators should be able to predict the effect an experimental intervention will have on survival by examining the data harvested during the study. This information might streamline drug development and drug approval.
Conception/Design: Wilfred D. Stein, William Doug Figg, William Dahut, Susan E. Bates, Tito Fojo Provision of study materials: William Doug Figg, William Dahut, Doug Price Collection/assembly of data: William Doug Figg, William Dahut, Doug Price Data analysis: Wilfred D. Stein, William Doug Figg, William Dahut, Aryeh D. Stein, Moshe B. Hoshen, Susan E. Bates, Tito Fojo Manuscript writing: Wilfred D. Stein, William Doug Figg, William Dahut, Aryeh D. Stein, Moshe B. Hoshen, Doug Price, Susan E. Bates, Tito Fojo Final approval of manuscript: Wilfred D. Stein, William Doug Figg, William Dahut, Aryeh D. Stein, Moshe B. Hoshen, Doug Price, Susan E. Bates, Tito Fojo
The authors thank Maha Hussain (University of Michigan Comprehensive Cancer Center, Ann Arbor, MI) and Cathy Tangen (Fred Hutchinson Cancer Center, Seattle, WA) for their kind help and cooperation in providing the data from the Southwest Oncology Group study used to validate the methodology and for reading and commenting on the manuscript. This work was supported by the intramural program of the National Cancer Institute. For the convenience of future users of this approach, a file in Excel into which data on tumor sizes at given times can be uploaded and the parameters g and d extracted is available in the Supplementary Material as an Excel File.
This article has been cited by other articles:
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||