In regenerative medicine, computer models describing bioreactor processes can assist in designing optimal process conditions leading to robust and economically viable products. In this study, we started from a (3D) mechanistic model describing the growth of neotissue, comprised of cells, and extracellular matrix, in a perfusion bioreactor set-up influenced by the scaffold geometry, flow-induced shear stress, and a number of metabolic factors. Subsequently, we applied model reduction by reformulating the problem from a set of partial differential equations into a set of ordinary differential equations. Comparing the reduced model results to the mechanistic model results and to dedicated experimental results assesses the reduction step quality. The obtained homogenized model is 10 fold faster than the 3D version, allowing the application of rigorous optimization techniques. Bayesian optimization was applied to find the medium refreshment regime in terms of frequency and percentage of medium replaced that would maximize neotissue growth kinetics during 21 days of culture. The simulation results indicated that maximum neotissue growth will occur for a high frequency and medium replacement percentage, a finding that is corroborated by reports in the literature. This study demonstrates an in silico strategy for bioprocess optimization paying particular attention to the reduction of the associated computational cost.

Tissue engineering and regenerative medicine looks at improving or restoring biological tissue function in humans and animals. We consider optimising neotissue growth in a three-dimensional scaffold during dynamic perfusion bioreactor culture, in the context of bone tissue engineering. The goal is to choose design variables that optimise two conflicting objectives: (i) maximising neotissue growth and (ii) minimising operating cost. We make novel extensions to Bayesian multi-objective optimisation in the case of one analytical objective function and one black-box, i.e. simulation-based, objective function. The analytical objective represents operating cost while the black-box neotissue growth objective comes from simulating a system of partial differential equations. The resulting multi-objective optimisation method determines the trade-off in the variables between neotissue growth and operating cost. Our method outperforms the most common approach in literature, genetic algorithms, in terms of data efficiency, on both the tissue engineering example and standard test functions. The resulting method is highly applicable to real-world problems combining black-box models with easy-to-quantify objectives like cost.

Background: Human mesenchymal stromal cells (hMSCs) have become attractive candidates for advanced medical cell-based therapies. An in vitro expansion step is routinely used to reach the required clinical quantities. However, this is influenced by many variables including donor characteristics such as age and gender, and culture conditions, such as cell seeding density and available culture surface area. Computational modeling in general and machine learning in particular could play a significant role in deciphering the relationship between the individual donor characteristics and their growth dynamics. Methods: In this study, hMSCs obtained from 174 male and female donors, ages ranging from 3 to 64 with passage numbers ranging from 2 to 27 were studied. We applied a Random Forests (RF) technique to model the cell expansion procedure by predicting the population doubling time (PDT) for each passage, taking into account individual donor-related characteristics. Results: Using the RF model, the mean absolute error between model predictions and experimental results for the PDT in passage 1 to 4 is significantly lower compared to the errors obtained with theoretical estimates or historic data. Moreover, statistical analysis indicate that the PD and PDT in different age categories are significantly different, especially the youngest group (aged <10 yrs) compared to the other age groups. Discussion: In summary, we introduce a predictive computational model describing in vitro cell expansion dynamics based on individual donor characteristics, an approach that could greatly assist toward automation of a cell expansion culture process.