Kinetic heterogeneity of an experimental tumour revealed by brdurd incorporation and mathematical modelling

In the present paper we propose a method of analysis of the cell kinetic characteristics of in vivo experimental tumours, that uses DNA-BrdUrd flow cytometry data at various times after the BrdUrd injection and mathematical modelling. The model of the cell population takes into account the cell-to-cell heterogeneity of the progression rate across cell cycle phases within the tumour, and assumes a strict correlation between the durations of S and G2M phases. The model also allows for a nonconstant DNA synthesis rate across S phase. In addition, the measurement process is modeled, considering the possibility of nonimpulsive labelling and providing a representation of the time course of the bivariate DNA-BrdUrd fluorescence distribution. Sequential DNA-BrdUrd distributions were obtained in vivo from a human ovarian carcinoma transplanted in mice and, for comparison, in vitro from a cell line of the same origin. From these data, that included the fractional density and the mean BrdUrd-fluorescence of BrdUrd-positive cells as a function of the DNA-fluorescence, kinetic parameters such as the potential doubling time (T_pot) and the mean and variance of the transit times in S and G2M phases, were estimated. This study revealed the presence of a substantial heterogeneity in S and G2M phases within the in vivo cell population and of a lower heterogeneity in the in vitro population. Moreover, our analysis suggests a nonnegligible effect of the BrdUrd pharmacokinetics in the in vivo cell labelling.