Unique osteogenic profile of bone marrow stem cells stimulated in perfusion bioreactor is Rho-ROCK-mediated contractility dependent

The fate determination of bone marrow mesenchymal stem/stromal cells (BMSC) is tightly regulated by mechanical cues, including fluid shear stress. Knowledge of mechanobiology in 2D culture has allowed researchers in bone tissue engineering to develop 3D dynamic culture systems with the potential for clinical translation in which the fate and growth of BMSC are mechanically controlled. However, due to the complexity of 3D dynamic cell culture compared to the 2D counterpart, the mechanisms of cell regulation in the dynamic environment remain relatively undescribed. In the present study, we analyzed the cytoskeletal modulation and osteogenic profiles of BMSC under fluid stimuli in a 3D culture condition using a perfusion bioreactor. BMSC subjected to fluid shear stress (mean 1.56 mPa) showed increased actomyosin contractility, accompanied by the upregulation of mechanoreceptors, focal adhesions, and Rho GTPase-mediated signaling molecules. Osteogenic gene expression profiling revealed that fluid shear stress promoted the expression of osteogenic markers differently from chemically induced osteogenesis. Osteogenic marker mRNA expression, type 1 collagen formation, ALP activity, and mineralization were promoted in the dynamic condition, even in the absence of chemical supplementation. The inhibition of cell contractility under flow by Rhosin chloride, Y27632, MLCK inhibitor peptide-18, or Blebbistatin revealed that actomyosin contractility was required for maintaining the proliferative status and mechanically induced osteogenic differentiation in the dynamic culture. The study highlights the cytoskeletal response and unique osteogenic profile of BMSC in this type of dynamic cell culture, stepping toward the clinical translation of mechanically stimulated BMCS for bone regeneration.

Measuring the density and viscosity of culture media for optimized computational fluid dynamics analysis of in vitro devices

Culture medium is frequently modelled as water in computational fluid dynamics (CFD) analysis of in vitro culture systems involving flow, such as bioreactors and organ-on-chips. However, culture medium can be expected to have different properties to water due to its higher solute content. Furthermore, cellular activities such as metabolism and secretion of ECM proteins alter the composition of culture medium and therefore its properties during culture. As these properties directly determine the hydromechanical stimuli exerted on cells in vitro, these, along with any changes during culture must be known for CFD modelling accuracy and meaningful interpretation of cellular responses. In this study, the density and dynamic viscosity of DMEM and RPMI-1640 media supplemented with typical concentrations of foetal bovine serum (0, 5, 10 and 20% v/v) were measured to serve as a reference for computational design analysis. Any changes in the properties of medium during culture were also investigated with NCI-H460 and HN6 cell lines. The density and dynamic viscosity of the media increased proportional to the % volume of added foetal bovine serum (FBS). Importantly, the viscosity of 5% FBS-supplemented RPMI-1640 was found to increase significantly after 3 days of culture of NCI-H460 and HN6 cell lines, with distinct differences between magnitude of change for each cell line. Finally, these experimentally-derived values were applied in CFD analysis of a simple microfluidic device, which demonstrated clear differences in maximum wall shear stress and pressure between fluid models. Overall, these results highlight the importance of characterizing model-specific properties for CFD design analysis of cell culture systems.

Predicting Bone Formation in Mesenchymal Stromal Cell-Seeded Hydrogels Using Experiment-Based Mathematical Modeling

In vitro bone formation by mesenchymal stromal cells encapsulated in type-1 collagen hydrogels is demonstrated after a 28-day in vitro culture period. Analysis of the hydrogels is carried out by X-ray microcomputed tomography, histology, and immunohistochemistry, which collectively demonstrates that bone formation in the hydrogels was quantifiably proportional to the initial collagen concentration, and subsequently the population density of seeded cells. This was established by varying the initial collagen concentration at a constant cell seeding density (3 × 10⁵ cells/0.3 mL hydrogel), and separately varying cell seeding density at a constant collagen concentration (1 mg/mL). Using these data, a mathematical model is presented for the total hydrogel volume and mineralization volume based on the observed linear contraction dynamics of cell-seeded collagen gels. The model parameters are fitted by comparing the predictions of the mathematical model for the hydrogel and mineralized volumes on day 28 with the experimental data. The model is then used to predict the hydrogel and mineralization volumes for a range of hydrogel collagen concentrations and cell seeding densities, providing comprehensive input/output descriptors for generating mineralized hydrogels for bone tissue engineering. It is proposed that this quantitative approach will be a useful tool for generating in vitro manufactured bone tissue, defining input parameters that yield predictable output measures of tissue maturation. Impact statement This article describes a simple yet powerful quantitative description of in vitro tissue-engineered bone by combining experimental data with mathematical modeling. The overall aim of the article is to examine what is currently known about cell-mediated collagen contraction, and demonstrate that this phenomenon can be exploited to tailor bone formation by choosing a specific set of input parameters in the form of cell seeding density and collagen hydrogel concentration. Our study utilizes a clinically relevant cell source (human mesenchymal stem cells) with a biomaterial that has received regulatory approval for use in humans (collagen type 1), and hence could be useful for clinical applications, as well as furthering our understanding of cell/extracellular matrix interactions in determining in vitro bone tissue formation.

Lattice and continuum modelling of a bioactive porous tissue scaffold

A contemporary procedure to grow artificial tissue is to seed cells onto a porous biomaterial scaffold and culture it within a perfusion bioreactor to facilitate the transport of nutrients to growing cells. Typical models of cell growth for tissue engineering applications make use of spatially homogeneous or spatially continuous equations to model cell growth, flow of culture medium, nutrient transport and their interactions. The network structure of the physical porous scaffold is often incorporated through parameters in these models, either phenomenologically or through techniques like mathematical homogenization. We derive a model on a square grid lattice to demonstrate the importance of explicitly modelling the network structure of the porous scaffold and compare results from this model with those from a modified continuum model from the literature. We capture two-way coupling between cell growth and fluid flow by allowing cells to block pores, and by allowing the shear stress of the fluid to affect cell growth and death. We explore a range of parameters for both models and demonstrate quantitative and qualitative differences between predictions from each of these approaches, including spatial pattern formation and local oscillations in cell density present only in the lattice model. These differences suggest that for some parameter regimes, corresponding to specific cell types and scaffold geometries, the lattice model gives qualitatively different model predictions than typical continuum models. Our results inform model selection for bioactive porous tissue scaffolds, aiding in the development of successful tissue engineering experiments and eventually clinically successful technologies.

Permeability versus Design in TPMS Scaffolds

Scaffolds for bone tissue engineering are porous structures that serve as support for cellular growth and, therefore, new tissue formation. The present work assessed the influence of the porous architecture of triply periodic minimal surface (TPMS) scaffolds on their macroscopic permeability behavior, combining numerical and experimental methods. The TPMS scaffolds considered were Schwartz D, Schwartz P, and Gyroid, which have been previously studied for bone tissue engineering, with 70% porosity. On the experimental side, these scaffolds were produced by MultiJet 3D printing and tested for fluid passage to calculate their permeability through Darcy’s Law. On the numerical side, finite element (FE) models of the scaffolds were simulated on ABAQUS^® for fluid passage under compression to assess potential fluid concentration spots. The outcomes revealed that the design of the unit cell had a noticeable effect on both calculated permeability and FE computed fluid flow velocity, regardless of the identical porosity, with the Gyroid scaffold having higher permeability and the Schwartz P a lower probability of fluid trapping. Schwartz D had the worst outcomes in both testing modalities, so these scaffolds would most likely be the last choice for promoting cell differentiation onto bone cells. Gyroid and Schwartz P would be up for selection depending on the application and targeted bone tissue.

The influence of hydrostatic pressure on tissue engineered bone development

The hydrostatic pressure stimulation of an appropriately cell-seeded porous scaffold within a bioreactor is a promising method for engineering bone tissue external to the body. We propose a mathematical model, and employ a suite of candidate constitutive laws, to qualitatively describe the effect of applied hydrostatic pressure on the quantity of minerals deposited in such an experimental setup. By comparing data from numerical simulations with experimental observations under a number of stimulation protocols, we suggest that the response of bone cells to an applied pressure requires consideration of two components; (i) a component describing the cell memory of the applied stimulation, and (ii) a recovery component, capturing the time cells require to recover from high rates of mineralisation.

Evaluation of the growth environment of a hydrostatic force bioreactor for preconditioning of tissue-engineered constructs

Bioreactors have been widely acknowledged as valuable tools to provide a growth environment for engineering tissues and to investigate the effect of physical forces on cells and cell-scaffold constructs. However, evaluation of the bioreactor environment during culture is critical to defining outcomes. In this study, the performance of a hydrostatic force bioreactor was examined by experimental measurements of changes in dissolved oxygen (O2), carbon dioxide (CO2), and pH after mechanical stimulation and the determination of physical forces (pressure and stress) in the bioreactor through mathematical modeling and numerical simulation. To determine the effect of hydrostatic pressure on bone formation, chick femur skeletal cell-seeded hydrogels were subjected to cyclic hydrostatic pressure at 0-270 kPa and 1 Hz for 1 h daily (5 days per week) over a period of 14 days. At the start of mechanical stimulation, dissolved O2 and CO2 in the medium increased and the pH of the medium decreased, but remained within human physiological ranges. Changes in physiological parameters (O2, CO2, and pH) were reversible when medium samples were placed in a standard cell culture incubator. In addition, computational modeling showed that the distribution and magnitude of physical forces depends on the shape and position of the cell-hydrogel constructs in the tissue culture format. Finally, hydrostatic pressure was seen to enhance mineralization of chick femur skeletal cell-seeded hydrogels.

Multiphase modelling of the effect of fluid shear stress on cell yield and distribution in a hollow fibre membrane bioreactor

We present a simplified two-dimensional model of fluid flow, nutrient transport and cell distribution in a hollow fibre membrane bioreactor, with the aim of exploring how fluid flow can be used to control the distribution and yield of a cell population which is sensitive to both fluid shear stress and nutrient concentration. The cells are seeded in a scaffold in a layer on top of the hollow fibre, only partially occupying the extracapillary space. Above this layer is a region of free-flowing fluid which we refer to as the upper fluid layer. The flow in the lumen and upper fluid layer is described by the Stokes equations, whilst the flow in the porous fibre membrane is assumed to follow Darcy’s law. Porous mixture theory is used to model the dynamics of and interactions between the cells, scaffold and fluid in the cell–scaffold construct. The concentration of a limiting nutrient (e.g. oxygen) is governed by an advection–reaction–diffusion equation in each region. Through exploitation of the small aspect ratio of each region and asymptotic analysis, we derive a coupled system of partial differential equations for the cell volume fraction and nutrient concentration. We use this model to investigate the effect of mechanotransduction on the distribution and yield of the cell population, by considering cases in which cell proliferation is either enhanced or limited by fluid shear stress and by varying experimentally controllable parameters such as flow rate and cell–scaffold construct thickness.

Mathematical model of growth factor driven haptotaxis and proliferation in a tissue engineering scaffold

Motivated by experimental work (Miller et al. in Biomaterials 27(10):2213-2221, 2006, 32(11):2775-2785, 2011) we investigate the effect of growth factor driven haptotaxis and proliferation in a perfusion tissue engineering bioreactor, in which nutrient-rich culture medium is perfused through a 2D porous scaffold impregnated with growth factor and seeded with cells. We model these processes on the timescale of cell proliferation, which typically is of the order of days. While a quantitative representation of these phenomena requires more experimental data than is yet available, qualitative agreement with preliminary experimental studies (Miller et al. in Biomaterials 27(10):2213-2221, 2006) is obtained, and appears promising. The ultimate goal of such modeling is to ascertain initial conditions (growth factor distribution, initial cell seeding, etc.) that will lead to a final desired outcome.

Modeling of cell cultures in perfusion bioreactors

Cultivating cells and tissues in bioreactors is a critical step in forming artificial tissues or organs prior to transplantation. Among various bioreactors, the perfusion bioreactor is known for its enhanced convection through the cell-scaffold constructs. Knowledge of mass transfer is essential for controlling the cell culture process; however, obtaining this information remains a challenging task. In this research, a novel mathematical model is developed to represent the nutrient transport and cell growth in a 3-D scaffold cultivated in a perfusion bioreactor. Numerical methods are employed to solve the equations involved, with a focus on identifying the effect of factors such as porosity, culturing time, and flow rate, which are controllable in the scaffold fabrication and culturing process, on cell cultures. To validate the new model, the results from the model simulations were compared to the experimental results extracted from the literature. With the validated model, further simulations were carried out to investigate the glucose and oxygen distribution and the cell growth within the cell-scaffold construct in a perfusion bioreactor, thus providing insight into the cell culture process.

Design Optimization of Scaffold Microstructures Using Wall Shear Stress Criterion Towards Regulated Flow-Induced Erosion

Tissue scaffolds aim to provide a cell-friendly biomechanical environment for facilitating cell growth. Existing studies have shown significant demands for generating a certain level of wall shear stress (WSS) on scaffold microstructural surfaces for promoting cellular response and attachment efficacy. Recently, its role in shear-induced erosion of polymer scaffold has also drawn increasing attention. This paper proposes a bi-directional evolutionary structural optimization (BESO) approach for design of scaffold microstructure in terms of the WSS uniformity criterion, by downgrading highly-stressed solid elements into fluidic elements and/or upgrading lowly-stressed fluidic elements into solid elements. In addition to this, a computational model is presented to simulate shear-induced erosion process. The effective stiffness and permeability of initial and optimized scaffold microstructures are characterized by the finite element based homogenization technique to quantify the variations of mechanical properties of scaffold during erosion. The illustrative examples show that a uniform WSS is achieved within the optimized scaffold microstructures, and their architectural and biomechanical features are maintained for a longer lifetime during shear-induced erosion process. This study provides a mathematical means to the design optimization of cellular biomaterials in terms of the WSS criterion towards controllable shear-induced erosion.