Towards unravelling the kinetics of an acute myeloid leukaemia model system under oxidative and starvation stress: a comparison between two- and three-dimensional cultures.

A great challenge when conducting ex vivo studies of leukaemia is the construction of an appropriate experimental platform that would recapitulate the bone marrow (BM) environment. Such a 3D scaffold system has been previously developed in our group [1]. Additionally to the BM architectural characteristics, parameters such as oxygen and glucose concentration are crucial as their value could differ between patients as well as within the same patient at different stages of treatment, consequently affecting the resistance of leukaemia to chemotherapy. The effect of oxidative and glucose stress-at levels close to human physiologic ones-on the proliferation and metabolic evolution of an AML model system (K-562 cell line) in conventional 2D cultures as well as in 3D scaffolds were studied. We observed that the K-562 cell line can proliferate and remain alive for 2 weeks in medium with glucose close to physiological levels both in 20 and 5% O2. We report interesting differences on the cellular response to the environmental, i.e., oxidative and/or nutritional stress stimuli in 2D and 3D. Higher adaptation to oxidative stress under non-starving conditions is observed in the 3D system. The glucose level in the medium has more impact on the cellular proliferation in the 3D compared to the 2D system. These differences can be of significant importance both when applying chemotherapy in vitro and also when constructing mathematical tools for optimisation of disease treatment.

Stem cell biomanufacturing under uncertainty: A case study in optimizing red blood cell production.

As breakthrough cellular therapy discoveries are translated into reliable, commercializable applications, effective stem cell biomanufacturing requires systematically developing and optimizing bioprocess design and operation. This article proposes a rigorous computational framework for stem cell biomanufacturing under uncertainty. Our mathematical tool kit incorporates: high-fidelity modeling, single variate and multivariate sensitivity analysis, global topological superstructure optimization, and robust optimization. The advantages of the proposed bioprocess optimization framework using, as a case study, a dual hollow fiber bioreactor producing red blood cells from progenitor cells were quantitatively demonstrated. The optimization phase reduces the cost by a factor of 4, and the price of insuring process performance against uncertainty is approximately 15% over the nominal optimal solution. Mathematical modeling and optimization can guide decision making; the possible commercial impact of this cellular therapy using the disruptive technology paradigm was quantitatively evaluated.

A mathematical model of subpopulation kinetics for the deconvolution of leukaemia heterogeneity.

Acute myeloid leukaemia is characterized by marked inter- and intra-patient heterogeneity, the identification of which is critical for the design of personalized treatments. Heterogeneity of leukaemic cells is determined by mutations which ultimately affect the cell cycle. We have developed and validated a biologically relevant, mathematical model of the cell cycle based on unique cell-cycle signatures, defined by duration of cell-cycle phases and cyclin profiles as determined by flow cytometry, for three leukaemia cell lines. The model was discretized for the different phases in their respective progress variables (cyclins and DNA), resulting in a set of time-dependent ordinary differential equations. Cell-cycle phase distribution and cyclin concentration profiles were validated against population chase experiments. Heterogeneity was simulated in culture by combining the three cell lines in a blinded experimental set-up. Based on individual kinetics, the model was capable of identifying and quantifying cellular heterogeneity. When supplying the initial conditions only, the model predicted future cell population dynamics and estimated the previous heterogeneous composition of cells. Identification of heterogeneous leukaemia clones at diagnosis and post-treatment using such a mathematical platform has the potential to predict multiple future outcomes in response to induction and consolidation chemotherapy as well as relapse kinetics.