Phenomenological modelling and simulation of cell clusters in 3D cultures

Extracellular matrix density regulates the formation of tumour spheroids through cell migration

In this work, we show how the mechanical properties of the cellular microenvironment modulate the growth of tumour spheroids. Based on the composition of the extracellular matrix, its stiffness and architecture can significantly vary, subsequently influencing cell movement and tumour growth. However, it is still unclear exactly how both of these processes are regulated by the matrix composition. Here, we present a centre-based computational model that describes how collagen density, which modulates the steric hindrance properties of the matrix, governs individual cell migration and, consequently, leads to the formation of multicellular clusters of varying size. The model was calibrated using previously published experimental data, replicating a set of experiments in which cells were seeded in collagen matrices of different collagen densities, hence producing distinct mechanical properties. At an initial stage, we tracked individual cell trajectories and speeds. Subsequently, the formation of multicellular clusters was also analysed by quantifying their size. Overall, the results showed that our model could accurately replicate what was previously seen experimentally. Specifically, we showed that cells seeded in matrices with low collagen density tended to migrate more. Accordingly, cells strayed away from their original cluster and thus promoted the formation of small structures. In contrast, we also showed that high collagen densities hindered cell migration and produced multicellular clusters with increased volume. In conclusion, this model not only establishes a relation between matrix density and individual cell migration but also showcases how migration, or its inhibition, modulates tumour growth.

Integration of in vitro and in silico Models Using Bayesian Optimization With an Application to Stochastic Modeling of Mesenchymal 3D Cell Migration

Cellular migration plays a crucial role in many aspects of life and development. In this paper, we propose a computational model of 3D migration that is solved by means of the tau-leaping algorithm and whose parameters have been calibrated using Bayesian optimization. Our main focus is two-fold: to optimize the numerical performance of the mechano-chemical model as well as to automate the calibration process of in silico models using Bayesian optimization. The presented mechano-chemical model allows us to simulate the stochastic behavior of our chemically reacting system in combination with mechanical constraints due to the surrounding collagen-based matrix. This numerical model has been used to simulate fibroblast migration. Moreover, we have performed in vitro analysis of migrating fibroblasts embedded in 3D collagen-based fibrous matrices (2 mg/ml). These in vitro experiments have been performed with the main objective of calibrating our model. Nine model parameters have been calibrated testing 300 different parametrizations using a completely automatic approach. Two competing evaluation metrics based on the Bhattacharyya coefficient have been defined in order to fit the model parameters. These metrics evaluate how accurately the in silico model is replicating in vitro measurements regarding the two main variables quantified in the experimental data (number of protrusions and the length of the longest protrusion). The selection of an optimal parametrization is based on the balance between the defined evaluation metrics. Results show how the calibrated model is able to predict the main features observed in the in vitro experiments.

Optimization of cell seeding in a 2D bio-scaffold system using computational models

The cell expansion process is a crucial part of generating cells on a large-scale level in a bioreactor system. Hence, it is important to set operating conditions (e.g. initial cell seeding distribution, culture medium flow rate) to an optimal level. Often, the initial cell seeding distribution factor is neglected and/or overlooked in the design of a bioreactor using conventional seeding distribution methods. This paper proposes a novel seeding distribution method that aims to maximize cell growth and minimize production time/cost. The proposed method utilizes two computational models; the first model represents cell growth patterns whereas the second model determines optimal initial cell seeding positions for adherent cell expansions. Cell growth simulation from the first model demonstrates that the model can be a representation of various cell types with known probabilities. The second model involves a combination of combinatorial optimization, Monte Carlo and concepts of the first model, and is used to design a multi-layer 2D bio-scaffold system that increases cell production efficiency in bioreactor applications. Simulation results have shown that the recommended input configurations obtained from the proposed optimization method are the most optimal configurations. The results have also illustrated the effectiveness of the proposed optimization method. The potential of the proposed seeding distribution method as a useful tool to optimize the cell expansion process in modern bioreactor system applications is highlighted.