Computational modeling of the expansion of human cord blood CD133+ hematopoietic stem/progenitor cells with different cytokine combinations

Motivation: Many important problems in cell biology require dense non-linear interactions between functional modules to be considered. The importance of computer simulation in understanding cellular processes is now widely accepted, and a variety of simulation algorithms useful for studying certain subsystems have been designed. Expansion of hematopoietic stem and progenitor cells (HSC/HPC) in ex vivo culture with cytokines and small molecules is a method to increase the restricted numbers of stem cells found in umbilical cord blood (CB), while also enhancing the content of early engrafting neutrophil and platelet precursors. The efficacy of the expanded product depends on the composition of the cocktail of cytokines and small molecules used for culture. Testing the influence of a cytokine or small molecule on the expansion of HSC/HPC is a laborious and expensive process. We therefore developed a computational model based on cellular signaling interactions that predict the influence of a cytokine on the survival, duplication and differentiation of the CD133+ HSC/HPC subset from human umbilical CB. Results: We have used results from in vitro expansion cultures with different combinations of one or more cytokines to develop an ordinary differential equation model that includes the effect of cytokines on survival, duplication and differentiation of the CD133+ HSC/HPC. Comparing the results of in vitro and in silico experiments, we show that the model can predict the effect of a cytokine on the fold expansion and differentiation of CB CD133+ HSC/HPC after 8-day culture on a 3D scaffol