Towards advanced bioprocess optimization: A multiscale modelling approach.

Mammalian cells produce up to 80 % of the commercially available therapeutic proteins, with Chinese Hamster Ovary (CHO) cells being the primary production host. Manufacturing involves a train of reactors, the last of which is typically run in fed-batch mode, where cells grow and produce the required protein. The feeding strategy is decided a priori, from either past operations or the design of experiments and rarely considers the current state of the process. This work proposes a Model Predictive Control (MPC) formulation based on a hybrid kinetic-stoichiometric reactor model to provide optimal feeding policies in real-time, which is agnostic to the culture, hence transferable across CHO cell culture systems. The benefits of the proposed controller formulation are demonstrated through a comparison between an open-loop simulation and closed-loop optimization, using a digital twin as an emulator of the process.

Osmotic behaviour of human mesenchymal stem cells: Implications for cryopreservation.

Aimed at providing a contribution to the optimization of cryopreservation processes, the present work focuses on the osmotic behavior of human mesenchymal stem cells (hMSCs). Once isolated from the umbilical cord blood (UCB) of three different donors, hMSCs were characterized in terms of size distribution and their osmotic properties suitably evaluated through the exposure to hypertonic and isotonic aqueous solutions at three different temperatures. More specifically, inactive cell volume and cell permeability to water and di-methyl sulfoxide (DMSO) were measured, being cell size determined using impedance measurements under both equilibrium and dynamic conditions. Experimental findings indicate that positive cell volume excursions are limited by the apparent increase of inactive volume, which occurs during both the shrink-swell process following DMSO addition and the subsequent restoration of isotonic conditions in the presence of hypertonic solutions of impermeant or permeant solutes. Based on this evidence, hMSCs must be regarded as imperfect osmometers, and their osmotic behavior described within a scenario no longer compatible with the simple two-parameter model usually utilized in the literature. In this respect, the activation of mechano-sensitive ion-channels seemingly represents a reasonable hypothesis for rationalizing the observed osmotic behavior of hMSCs from UCB.

A novel population balance model to investigate the kinetics of in vitro cell proliferation: part I. Model development.

In biotechnology and biomedicine reliable models of cell proliferation kinetics need to capture the relevant phenomena taking place during the mitotic cycle. To this aim, a novel mathematical model helpful to investigate the intrinsic kinetics of in vitro culture of adherent cells up to confluence is proposed in this work. Specifically, the attention is focused on the simulation of proliferation (increase of cell number) and maturation (increase of cell size and DNA content) till contact inhibition eventually takes place inside a Petri dish. Accordingly, the proposed model is based on a population balance (PB) approach that allows one to quantitatively describe cell cycle progression through the different phases the cells of the entire population experienced during their own life. In particular, the proposed model has been developed as a 2D, multi-staged, and unstructured PB, by considering a different sub-population of cells for any single phase of the cell cycle. These sub-populations are discriminated through cellular volume and DNA content, that both increase during the mitotic cycle. The adopted mathematical expressions of the transition rates between two subsequent phases and the temporal increase of cell volume and DNA content are thoroughly analyzed and discussed with respect to those ones available in the literature. Specifically, the corresponding uncertainties and pitfalls are pointed out, by also taking into account the difficulties and the limitations involved in the quantitative measurements currently practicable for these biological systems. A novel mathematical expression for contact inhibition in line with the PB model developed is also formulated, along with a proper comparison between modeled and measurable DNA distributions. The strategy for a reliable, independent tuning of the adjustable parameters involved in the proposed model along with its numerical solution is outlined in Part II of this work, where it is also shown how it can be profitably used to gain a deeper insight into the phenomena involved during cell cultivation under microgravity conditions.

A novel population balance model to investigate the kinetics of in vitro cell proliferation: part II. Numerical solution, parameters' determination, and model outcomes.

Based on the general theoretical model developed in Part I of this work, a series of numerical simulations related to the in vitro proliferation kinetics of adherent cells is here presented. First the complex task of assigning a specific value to all the parameters of the proposed population balance (PB) model is addressed, by also highlighting the difficulties arising when performing proper comparisons with experimental data. Then, a parametric sensitivity analysis is performed, thus identifying the more relevant parameters from a kinetics perspective. The proposed PB model can be adapted to describe cell growth under various conditions, by properly changing the value of the adjustable parameters. For this reason, model parameters able to mimic cell culture behavior under microgravity conditions are identified by means of a suitable parametric sensitivity analysis. Specifically, it is found that, as the volume growth parameter is reduced, proliferation slows down while cells arrest in G0/G1 or G2/M depending on the initial distribution of cell population. On the basis of this result, model capabilities have been tested by means of a proper comparison with literature experimental data related to the behavior of synchronized and not-synchronized cells under micro- and standard gravity levels.

A simulation model for stem cells differentiation into specialized cells of non-connective tissues.

A novel mathematical model to simulate stem cells differentiation into specialized cells of non-connective tissues is proposed. The model is based upon material balances for growth factors coupled with a mass-structured population balance describing cell growth, proliferation and differentiation. The proposed model is written in a general form and it may be used to simulate a generic cell differentiation pathway during in vitro cultivation when specific growth factors are used. Literature experimental data concerning the differentiation of central nervous stem cells into astrocytes are successfully compared with model results, thus demonstrating the validity of the proposed model as well as its predictive capability. Finally, sensitivity analysis of model parameters is also performed in order to clarify what mechanisms most strongly influence differentiation and cell types distribution.