Coupling curvature-dependent and shear stress-stimulated neotissue growth in dynamic bioreactor cultures: a 3D computational model of a complete scaffold

The main challenge in tissue engineering consists in understanding and controlling the growth process of in vitro cultured neotissues toward obtaining functional tissues. Computational models can provide crucial information on appropriate bioreactor and scaffold design but also on the bioprocess environment and culture conditions. In this study, the development of a 3D model using the level set method to capture the growth of a microporous neotissue domain in a dynamic culture environment (perfusion bioreactor) was pursued. In our model, neotissue growth velocity was influenced by scaffold geometry as well as by flow- induced shear stresses. The neotissue was modeled as a homogenous porous medium with a given permeability, and the Brinkman equation was used to calculate the flow profile in both neotissue and void space. Neotissue growth was modeled until the scaffold void volume was filled, thus capturing already established experimental observations, in particular the differences between scaffold filling under different flow regimes. This tool is envisaged as a scaffold shape and bioprocess optimization tool with predictive capacities. It will allow controlling fluid flow during long-term culture, whereby neotissue growth alters flow patterns, in order to provide shear stress profiles and magnitudes across the whole scaffold volume influencing, in turn, the neotissue growth.     

Flow rates in perfusion bioreactors to maximise mineralisation in bone tissue engineering in vitro

In bone tissue engineering experiments, fluid-induced shear stress is able to stimulate cells to produce mineralised extracellular matrix (ECM). The application of shear stress on seeded cells can for example be achieved through bioreactors that perfuse medium through porous scaffolds. The generated mechanical environment (i.e. wall shear stress: WSS) within the scaffolds is complex due to the complexity of scaffold geometry. This complexity has so far prevented setting an optimal loading (i.e. flow rate) of the bioreactor to achieve an optimal distribution of WSS for stimulating cells to produce mineralised ECM. In this study, we demonstrate an approach combining computational fluid dynamics (CFD) and mechano-regulation theory to optimise flow rates of a perfusion bioreactor and various scaffold geometries (i.e. pore shape, porosity and pore diameter) in order to maximise shear stress induced mineralisation. The optimal flow rates, under which the highest fraction of scaffold surface area is subjected to a wall shear stress that induces mineralisation, are mainly dependent on the scaffold geometries. Nevertheless, the variation range of such optimal flow rates are within 0.5–5 mL/min (or in terms of fluid velocity: 0.166–1.66 mm/s), among different scaffolds. This approach can facilitate the determination of scaffold-dependent flow rates for bone tissue engineering experiments in vitro, avoiding performing a series of trial and error experiments.     

3-D computational modeling of media flow through scaffolds in a perfusion bioreactor

Media perfusion bioreactor systems have been developed to improve mass transport throughout three-dimensional (3-D) tissue-engineered constructs cultured in vitro. In addition to enhancing the exchange of nutrients and wastes, these systems simultaneously deliver flow-mediated shear stresses to cells seeded within the constructs. Local shear stresses are a function of media flow rate and dynamic viscosity, bioreactor configuration, and porous scaffold microarchitecture. We have used the Lattice-Boltzmann method to simulate the flow conditions within perfused cell-seeded cylindrical scaffolds. Microcomputed tomography imaging was used to define the scaffold microarchitecture for the simulations, which produce a 3-D fluid velocity field throughout the scaffold porosity. Shear stresses were estimated at various media flow rates by multiplying the symmetric part of the gradient of the velocity field by the dynamic viscosity of the cell culture media. The shear stress algorithm was validated by modeling flow between infinite parallel plates and comparing the calculated shear stress distribution to the analytical solution. Relating the simulation results to perfusion experiments, an average surface shear stress of 5x10(-5)Pa was found to correspond to increased cell proliferation, while higher shear stresses were associated with upregulation of bone marker genes. This modeling approach can be used to compare results obtained for different perfusion bioreactor systems or different scaffold microarchitectures and may allow specific shear stresses to be determined that optimize the amount, type, or distribution of in vitro tissue growth.     

Cells in 3D matrices under interstitial flow: Effects of extracellular matrix alignment on cell shear stress and drag forces

Interstitial flow is an important regulator of various cell behaviors both in vitro and in vivo, yet the forces that fluid flow imposes on cells embedded in a 3D extracellular matrix (ECM), and the effects of matrix architecture on those forces, are not well understood. Here, we demonstrate how fiber alignment can affect the shear and pressure forces on the cell and ECM. Using computational fluid dynamics simulations, we show that while the solutions of the Brinkman equation accurately estimate the average fluid shear stress and the drag forces on a cell within a 3D fibrous medium, the distribution of shear stress on the cellular surface as well as the peak shear stresses remain intimately related to the pericellular fiber architecture and cannot be estimated using bulk-averaged properties. We demonstrate that perpendicular fiber alignment of the ECM yields lower shear stress and pressure forces on the cells and higher stresses on the ECM, leading to decreased permeability, while parallel fiber alignment leads to higher stresses on cells and increased permeability, as compared to a cubic lattice arrangement. The Spielman–Goren permeability relationships for fibrous media agreed well with CFD simulations of flow with explicitly considered fibers. These results suggest that the experimentally observed active remodeling of ECM fibers by fibroblasts under interstitial flow to a perpendicular alignment could serve to decrease the shear and drag forces on the cell.     

A dynamical study of the mechanical stimuli and tissue differentiation within a CaP scaffold based on micro-CT finite element models

The control of the mechanical stimuli transmitted to the cells is critical for the design of functional scaffolds for tissue engineering. The objective of this study was to investigate the dynamics of the mechanical stimuli transmitted to the cells during tissue differentiation in an irregular morphology scaffold under compressive load and perfusion flow. A calcium phosphate-based glass porous scaffold was used. The solid phase and the fluid flow within the pores were modeled as linear elastic solid material and Newtonian fluid, respectively. In the fluid model, different levels of viscosity were used to simulate tissue differentiation. Compressive strain of 0.5% and fluid flow with constant inlet velocity of 10 μm/s or constant inlet pressure of 3 Pa were applied. Octahedral shear strain and fluid shear stress were used as mechano-regulatory stimuli. For constant inlet velocity, stimuli equivalent to bone were predicted in 80% of pore volume for the case of low tissue viscosity. For the cases of high viscosity, fluctuations between stimuli equivalent to tissue formation and cell death were predicted due to the increase in the fluid shear stress when tissue started to fill pores. When constant pressure was applied, stimuli equivalent to bone were predicted in 62% of pore volume when low tissue viscosity was used and 42% when high tissue viscosity was used. This study predicted critical variations of fluid shear stress when cells differentiated. If these variations are not controlled in vitro, they can impede the formation of new matured tissue.     

Modeling and design of optimal flow perfusion bioreactors for tissue engineering applications

Perfusion bioreactors have been used in different tissue engineering applications because of their consistent distribution of nutrients and flow-induced shear stress within the tissue-engineering scaffold. A widely used configuration uses a scaffold with a circular cross-section enclosed within a cylindrical chamber and inlet and outlet pipes which are connected to the chamber on either side through which media is continuously circulated. However, fluid-flow experiments and simulations have shown that the majority of the flow perfuses through the center. This pattern creates stagnant zones in the peripheral regions as well as in those of high flow rate near the inlet and outlet. This non-uniformity of flow and shear stress, owing to a circular design, results in limited cell proliferation and differentiation in these areas. The focus of this communication is to design an optimized perfusion system using computational fluid dynamics as a mathematical tool to overcome the time-consuming trial and error experimental method. We compared the flow within a circular and a rectangular bioreactor system. Flow simulations within the rectangular bioreactor are shown to overcome the limitations in the circular design. This communication challenges the circular cross-section bioreactor configuration paradigm and provides proof of the advantages of the new design over the existing one.     

Computational modeling of flow-induced shear stresses within 3D salt-leached porous scaffolds imaged via micro-CT

Flow-induced shear stresses have been found to be a stimulatory factor in pre-osteoblastic cells seeded in 3D porous scaffolds and cultured under continuous flow perfusion. However, due to the complex internal structure of porous scaffolds, analytical estimation of the local shear forces is impractical. The primary goal of this work is to investigate the shear stress distributions within Poly(l-lactic acid) scaffolds via computation. Scaffolds used in this study are prepared via salt leeching with various geometric characteristics (80–95% porosity and 215–402.5 μm average pore size). High resolution micro-computed tomography is used to obtain their 3D structure. Flow of osteogenic media through the scaffolds is modeled via lattice Boltzmann method. It is found that the surface stress distributions within the scaffolds are characterized by long tails to the right (a positive skewness). Their shape is not strongly dependent on the scaffold manufacturing parameters, but the magnitudes of the stresses are. Correlations are prepared for the estimation of the average surface shear stress experienced by the cells within the scaffolds and of the probability density function of the surface stresses. Though the manufacturing technique does not appear to affect the shape of the shear stress distributions, presence of manufacturing defects is found to be significant: defects create areas of high flow and high stress along their periphery. The results of this study are applicable to other polymer systems provided that they are manufactured by a similar salt leeching technique, while the imaging/modeling approach is applicable to all scaffolds relevant to tissue engineering.     

Modeling evaluation of the fluid-dynamic microenvironment in tissue-engineered constructs: a micro-CT based model

Natural cartilage remodels both in vivo and in vitro in response to mechanical stresses, hence mechanical stimulation is believed to be a potential tool to modulate extra-cellular matrix synthesis in tissue-engineered cartilage. Fluid-induced shear is known to enhance chondrogenesis in engineered cartilage constructs. The quantification of the hydrodynamic environment is a condition required to study the biochemical response to shear of 3D engineered cell systems. We developed a computational model of culture medium flow through the microstructure of a porous scaffold, during direct- perfused culture. The 3D solid model of the scaffold micro-geometry was reconstructed from 250 micro-computed tomography (micro-CT) images. The results of the fluid dynamic simulations were analyzed at the central portions of the fluid domain, to avoid boundary effects. The average, median and mode shear stress values calculated at the scaffold walls were 3.48, 2.90, and 2.45 mPa respectively, at a flow rate of 0.5 cm(3)/min, perfused through a 15 mm diameter scaffold, at an inlet fluid velocity of 53 microm/s. These results were compared to results estimated using a simplified micro-scale model and to results estimated using an analytical macro-scale porous model. The predictions given by the CT-based model are being used in conjunction with an experimental bioreactor model, in order to quantify the effects of fluid-dynamic shear on the growth modulation of tissue-engineered cartilage constructs, to potentially enhance tissue growth in vitro.     

Prediction of the micro-fluid dynamic environment imposed to three-dimensional engineered cell systems in bioreactors

Bioreactors allowing culture medium perfusion overcome diffusion limitations associated with static culturing and provide flow-mediated mechanical stimuli. The hydrodynamic stress imposed to cells will depend not only on the culture medium flow rate, but also on the scaffold three-dimensional (3D) micro-architecture. We developed a CFD model of the flow of culture medium through a 3D scaffold of homogeneous geometry, with the aim of predicting the shear stress acting on cells as a function of parameters that can be controlled during the scaffold fabrication process, such as the scaffold porosity and the pore size, and during the cell culture, such as the medium flow rate and the diameter of the perfused scaffold section. We built three groups of models corresponding to three pore sizes: 50, 100 and 150 microm. Each group was made of four models corresponding to 59%, 65%, 77%, and 89% porosity. A commercial finite-element code was used to set up and solve the problem and to analyze the results. The mode value of shear stress varied between 2 and 5 mPa, and was obtained for a circular scaffold of 15.5 mm diameter, perfused by a flow rate of 0.5 ml/min. The simulations showed that the pore size is a variable strongly influencing the predicted shear stress level, whereas the porosity is a variable strongly affecting the statistical distribution of the shear stresses, but not their magnitude. Our results provide a basis for the completion of more exhaustive quantitative studies to further assess the relationship between perfusion, at known micro-fluid dynamic conditions, and tissue growth in vitro.     

A method for the design of 3D scaffolds for high-density cell attachment and determination of optimum perfusion culture conditions

The application of in vitro cultured cells in tissue engineering or drug screening, aimed at complex soft tissues such as liver, requires in vivo physiological function of the cultured cells. For this purpose, the scaffold in which cells are cultured should provide a microenvironment similar to an in vivo one with a three-dimensional extracellular matrix, a high supply capacity of O₂ and nutrients, and high cell density. In this paper, we propose a method to design (1) the geometry of the scaffold, with a surface/volume ratio optimized to allow high-density (5×10⁷ cells/mL) cell culture and (2) culture conditions that will supply optimal quantities of oxygen and nutrients. CFD modeling of mass transport was used to determine the shear stress as well as O₂ and glucose metabolism in the scaffold (20 mm width–35 mm length) for various flow rates. Validation of the model was done through comparison with flow resistance and micro-PIV experiments. CFD analysis showed the maximum metabolic rate densities for this scaffold are 6.04×10⁻³ mol/s/m³ for O₂ at 0.71 mL/min and 1.91×10⁻² mol/s/m³ for glucose at 0.35 mL/min.     

Modeling of flow‐induced shear stress applied on 3D cellular scaffolds: Implications for vascular tissue engineering

Novel tissue‐culture bioreactors employ flow‐induced shear stress as a means of mechanical stimulation of cells. We developed a computational fluid dynamics model of the complex three‐dimensional (3D) microstructure of a porous scaffold incubated in a direct perfusion bioreactor. Our model was designed to predict high shear‐stress values within the physiological range of those naturally sensed by vascular cells (1–10 dyne/cm²), and will thereby provide suitable conditions for vascular tissue‐engineering experiments. The model also accounts for cellular growth, which was designed as an added cell layer grown on all scaffold walls. Five model variants were designed, with geometric differences corresponding to cell‐layer thicknesses of 0, 50, 75, 100, and 125 µm. Four inlet velocities (0.5, 1, 1.5, and 2 cm/s) were applied to each model. Wall shear‐stress distribution and overall pressure drop calculations were then used to characterize the relation between flow rate, shear stress, cell‐layer thickness, and pressure drop. The simulations showed that cellular growth within 3D scaffolds exposes cells to elevated shear stress, with considerably increasing average values in correlation to cell growth and inflow velocity. Our results provide in‐depth analysis of the microdynamic environment of cells cultured within 3D environments, and thus provide advanced control over tissue development in vitro. Biotechnol. Bioeng. 2010; 105: 645–654. © 2009 Wiley Periodicals, Inc.     

Curvature- and fluid-stress-driven tissue growth in a tissue-engineering scaffold pore

Cell proliferation within a fluid-filled porous tissue-engineering scaffold depends on a sensitive choice of pore geometry and flow rates: regions of high curvature encourage cell proliferation, while a critical flow rate is required to promote growth for certain cell types. When the flow rate is too slow, the nutrient supply is limited; when it is too fast, cells may be damaged by the high fluid shear stress. As a result, determining appropriate tissue-engineering-construct geometries and operating regimes poses a significant challenge that cannot be addressed by experimentation alone. In this paper, we present a mathematical theory for the fluid flow within a pore of a tissue-engineering scaffold, which is coupled to the growth of cells on the pore walls. We exploit the slenderness of a pore that is typical in such a scenario, to derive a reduced model that enables a comprehensive analysis of the system to be performed. We derive analytical solutions in a particular case of a nearly piecewise constant growth law and compare these with numerical solutions of the reduced model. Qualitative comparisons of tissue morphologies predicted by our model, with those observed experimentally, are also made. We demonstrate how the simplified system may be used to make predictions on the design of a tissue-engineering scaffold and the appropriate operating regime that ensures a desired level of tissue growth.

     

Computational fluid dynamics modeling of steady-state momentum and mass transport in a bioreactor for cartilage tissue engineering

Computational fluid dynamics (CFD) models to quantify momentum and mass transport under conditions of tissue growth will aid bioreactor design for development of tissue-engineered cartilage constructs. Fluent CFD models are used to calculate flow fields, shear stresses, and oxygen profiles around nonporous constructs simulating cartilage development in our concentric cylinder bioreactor. The shear stress distribution ranges from 1.5 to 12 dyn/cm(2) across the construct surfaces exposed to fluid flow and varies little with the relative number or placement of constructs in the bioreactor. Approximately 80% of the construct surface exposed to flow experiences shear stresses between 1.5 and 4 dyn/cm(2), validating the assumption that the concentric cylinder bioreactor provides a relatively homogeneous hydrodynamic environment for construct growth. Species mass transport modeling for oxygen demonstrates that fluid-phase oxygen transport to constructs is uniform. Some O(2) depletion near the down stream edge of constructs is noted with minimum pO(2) values near the constructs of 35 mmHg (23% O(2) saturation). These values are above oxygen concentrations in cartilage in vivo, suggesting that bioreactor oxygen concentrations likely do not affect chondrocyte growth. Scale-up studies demonstrate the utility and flexibility of CFD models to design and characterize bioreactors for growth of tissue-engineered cartilage.     

Flow in the well: computational fluid dynamics is essential in flow chamber construction

A perfusion system was developed to generate well defined flow conditions within a well of a standard multidish. Human vein endothelial cells were cultured under flow conditions and cell response was analyzed by microscopy. Endothelial cells became elongated and spindle shaped. As demonstrated by computational fluid dynamics (CFD), cells were cultured under well defined but time varying shear stress conditions. A damper system was introduced which reduced pulsatile flow when using volumetric pumps. The flow and the wall shear stress distribution were analyzed by CFD for the steady and unsteady flow field. Usage of the volumetric pump caused variations of the wall shear stresses despite the controlled fluid environment and introduction of a damper system. Therefore the use of CFD analysis and experimental validation is critical in developing flow chambers and studying cell response to shear stress. The system presented gives an effortless flow chamber setup within a 6-well standard multidish.     

Permeability analysis of scaffolds for bone tissue engineering

Porous artificial bone substitutes, especially bone scaffolds coupled with osteobiologics, have been developed as an alternative to the traditional bone grafts. The bone scaffold should have a set of properties to provide mechanical support and simultaneously promote tissue regeneration. Among these properties, scaffold permeability is a determinant factor as it plays a major role in the ability for cells to penetrate the porous media and for nutrients to diffuse. Thus, the aim of this work is to characterize the permeability of the scaffold microstructure, using both computational and experimental methods. Computationally, permeability was estimated by homogenization methods applied to the problem of a fluid flow through a porous media. These homogenized permeability properties are compared with those obtained experimentally. For this purpose a simple experimental setup was used to test scaffolds built using Solid Free Form techniques. The obtained results show a linear correlation between the computational and the experimental permeability. Also, this study showed that permeability encompasses the influence of both porosity and pore size on mass transport, thus indicating its importance as a design parameter. This work indicates that the mathematical approach used to determine permeability may be useful as a scaffold design tool.     

A multiscale computational fluid dynamics approach to simulate the micro-fluidic environment within a tissue engineering scaffold with highly irregular pore geometry

Mechanical stimulation can regulate cellular behavior, e.g., differentiation, proliferation, matrix production and mineralization. To apply fluid-induced wall shear stress (WSS) on cells, perfusion bioreactors have been commonly used in tissue engineering experiments. The WSS on cells depends on the nature of the micro-fluidic environment within scaffolds under medium perfusion. Simulating the fluidic environment within scaffolds will be important for gaining a better insight into the actual mechanical stimulation on cells in a tissue engineering experiment. However, biomaterial scaffolds used in tissue engineering experiments typically have highly irregular pore geometries. This complexity in scaffold geometry implies high computational costs for simulating the precise fluidic environment within the scaffolds. In this study, we propose a low-computational cost and feasible technique for quantifying the micro-fluidic environment within the scaffolds, which have highly irregular pore geometries. This technique is based on a multiscale computational fluid dynamics approach. It is demonstrated that this approach can capture the WSS distribution in most regions within the scaffold. Importantly, the central process unit time needed to run the model is considerably low.

     

A “sweet-spot” for fluid-induced oscillations in the conditioning of stem cell-based engineered heart valve tissues

Fluid-induced shear stresses are involved in the development of cardiovascular tissues. In a tissue engineering framework, this stimulus has also been considered as a mechanical regulator of stem cell differentiation. We recently demonstrated that the fluid-oscillating effect in combination with a physiologically-relevant shear stress magnitude contributes to the formation of stem cell-derived de novo heart valve tissues. However, the range of oscillations necessary to induce favorable gene expression and engineered tissue formation is unknown. In this study, we took a computational approach to establish a range of oscillatory shear stresses that may optimize in vitro valvular tissue growth. Taking a biomimetic approach, three physiologically-relevant flow waveforms from the human: (i) aorta, (ii) pulmonary artery and (iii) superior vena cava were utilized to simulate pulsatile flow conditions within a bioreactor that housed 3 tissue specimens. Results were compared to non-physiological pulsatile flow (NPPF) and cyclic flexure-steady flow (Flex-Flow) conditions. The oscillatory shear index (OSI) was used to quantify the fluid-induced oscillations occurring on the specimen surfaces. The range of mean OSI under the physiological conditions investigated was found to be 0.18 ≤ OSI ≤ 0.23. On the other hand, NPPF and Flex-Flow environments yielded a mean OSI of 0.37 and 0.11 respectively, which were 46% higher and 45% lower than physiological conditions. Moreover, we subsequently conducted OSI-based human bone marrow stem cell (HBMSC) culture experiments which resulted in preferential valvular gene expression and phenotype (significant upregulation of BMP, KLF2A, CD31 and α-SMA using an OSI of 0.23 in comparison to a lower OSI of 0.10 or a higher OSI of 0.38; p < .05). These findings suggest that a distinct range or a “sweet-spot” for physiological OSI exists in the mechanical conditioning of tissue engineered heart valves grown from stem cell sources. We conclude that in vitro heart valve matrix development could be further enhanced by simultaneous exposure of the engineered tissues to physiologically-relevant magnitudes of both fluid-induced oscillations and shear stresses.     

Breaking symmetry in non-planar bifurcations: distribution of flow and wall shear stress

Non-planarity in blood vessels is known to influence arterial flows and wall shear stress. To gain insight, computational fluid dynamics (CFD) has been used to investigate effects of curvature and out-of-plane geometry on the distribution of fluid flows and wall shear stresses in a hypothetical non-planar bifurcation. Three-dimensional Navier-Stokes equations for a steady state Newtonian fluid were solved numerically using a finite element method. Non-planarity in one of the two daughter vessels is found to deflect flow from the inner wall of the vessel to the outer wall and to cause changes in the distribution of wall shear stresses. Results from this study agree to experimental observations and CFD simulations in the literature, and support the view that non-planarity in blood vessels is a factor with important haemodynamic significance and may play a key role in vascular biology and pathophysiology.     

A multiphysics/multiscale 2D numerical simulation of scaffold-based cartilage regeneration under interstitial perfusion in a bioreactor

In vitro tissue engineering is investigated as a potential source of functional tissue constructs for cartilage repair, as well as a model system for controlled studies of cartilage development and function. Among the different kinds of devices for the cultivation of 3D cartilage cell colonies, we consider here polymeric scaffold-based perfusion bioreactors, where an interstitial fluid supplies nutrients and oxygen to the growing biomass. At the same time, the fluid-induced shear acts as a physiologically relevant stimulus for the metabolic activity of cells, provided that the shear stress level is appropriately tuned. In this complex environment, mathematical and computational modeling can help in the optimal design of the bioreactor configuration. In this perspective, we propose a computational model for the simulation of the biomass growth, under given inlet and geometrical conditions, where nutrient concentration, fluid dynamic field and cell growth are consistently coupled. The biomass growth model is calibrated with respect to the shear stress dependence on experimental data using a simplified short-time analysis in which the nutrient concentration and the fluid-induced shear stress are assumed constant in time and uniform in space. Volume averaging techniques are used to derive effective parameters that allow to upscale the microscopic structural properties to the macroscopic level. The biomass growth predictions obtained in this way are significant for long times of culture.