Empirical models of the proliferative response of cytokine-dependent hematopoietic cell lines

There is an expanding need for predictive mathematical models to accelerate the optimization of cell therapy culture processes. Here we demonstrate the ability of simple mathematical models to describe quantitatively the cytokine growth-rate dependence of two human hematopoietic cell lines, TF-1 and MO7e. These cells are immortal but depend on either interleukin-3 (IL-3) or granulocyte-macrophage colony stimulating factor (GM-CSF) for their continued survival and maximal proliferation. They are also responsive to interleukin-6 (IL-6) and exhibit saturation kinetics when these cytokines are limiting. A Monod-type relationship consistently failed to fit measured cytokine dose-proliferation response curves while a Hill-type relationship showed a good fit. Cytokine interactions were first modeled by modifying the Hill-function to include an interaction parameter, gamma. This model did not indicate either synergistic or even additive effects between IL-3 and GM-CSF. Based on the reported competition between IL-3 and GM-CSF for their common receptor (beta(c)) subunit, a competitive model was also developed. This model had no new parameters beyond those obtained from single cytokine cultures and provided improved prediction of the growth rates for both cell lines exposed to combinations of IL-3 and GM-CSF over a wide range of concentrations. As expected, the competitive model failed to fit the data for IL-6 in combination with either IL-3 or GM-CSF, since IL-6 signaling does not involve the beta(c) chain of the IL-3/GM-CSF receptors. Interestingly, the cell-specific rates of GM-CSF uptake and cell proliferation were found to be uncoupled processes. Taken together, these results illustrate the utility of appropriately designed empirical models to describe the proliferative responses of hematopoietic cells to cytokine stimulation.