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.     

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.     

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.     

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.     

Quantification of fluid shear stress in bone tissue engineering scaffolds with spherical and cubical pore architectures

Recent studies have shown that mechanical stimulation, in the form of fluid perfusion and mechanical compression, can enhance osteogenic differentiation of mesenchymal stem cells and bone cells within tissue engineering scaffolds in vitro. The precise nature of mechanical stimulation within tissue engineering scaffolds is not only dictated by the exogenously applied loading regime, but also depends on the geometric features of the scaffold, in particular architecture, pore size and porosity. However, the precise contribution of each geometric feature towards the resulting mechanical stimulation within a scaffold is difficult to characterise due to the wide range of interacting parameters. In this study, we have applied a fluid–structure interaction model to investigate the role of scaffold geometry (architecture, pore size and porosity) on pore wall shear stress (WSS) under a range of different loading scenarios: fluid perfusion, mechanical compression and a combination of perfusion and compression. It is found that scaffold geometry (spherical and cubical pores), in particular the pore size, has a significant influence on the stimulation within scaffolds. Furthermore, we observed an amplified WSS within scaffolds under a combination of fluid perfusion and mechanical compression, which exceeded that caused by individual fluid perfusion or mechanical compression approximately threefold. By conducting this comprehensive parametric variation study, an expression was generated to allow the design and optimisation of 3D TE scaffolds and inform experimental loading regimes so that a desired level of mechanical stimulation, in terms of WSS is generated within the scaffold.     

Three-dimensional cell seeding and growth in radial-flow perfusion bioreactor for in vitro tissue reconstruction

Radial-flow perfusion bioreactor systems have been designed and evaluated to enable direct cell seeding into a three-dimensional (3-D) porous scaffold and subsequent cell culture for in vitro tissue reconstruction. However, one of the limitations of in vitro regeneration is the tissue necrosis that occurs at the central part of the 3-D scaffold. In the present study, tubular poly-L-lactic acid (PLLA) porous scaffolds with an optimized pore size and porosity were prepared by the lyophilization method, and the effect of different perfusion conditions on cell seeding and growth were compared with those of the conventional static culture. The medium flowed radially from the lumen toward the periphery of the tubular scaffolds. It was found that cell seeding under a radial-flow perfusion condition of 1.1 mL/cm2 x min was effective, and that the optimal flow rate for cell growth was 4.0 mL/cm2 x min. At this optimal rate, the increase in seeded cells in the perfusion culture over a period of 5 days was 7.3-fold greater than that by static culture over the same period. The perfusion cell seeding resulted in a uniform distribution of cells throughout the scaffold. Subsequently, the perfusion of medium and hence the provision of nutrients and oxygen permitted growth and maintenance of the tissue throughout the scaffold. The perfusion seeding/culture system was a much more effective strategy than the conventional system in which cells are seeded under a static condition and cultured in a bioreactor such as a spinner flask.     

Perfusion flow bioreactor for 3D in situ imaging: investigating cell/biomaterials interactions

The capability to image real time cell/material interactions in a three-dimensional (3D) culture environment will aid in the advancement of tissue engineering. This paper describes a perfusion flow bioreactor designed to hold tissue engineering scaffolds and allow for in situ imaging using an upright microscope. The bioreactor can hold a scaffold of desirable thickness for implantation (>2 mm). Coupling 3D culture and perfusion flow leads to the creation of a more biomimetic environment. We examined the ability of the bioreactor to maintain cell viability outside of an incubator environment (temperature and pH stability), investigated the flow features of the system (flow induced shear stress), and determined the image quality in order to perform time-lapsed imaging of two-dimensional (2D) and 3D cell culture. In situ imaging was performed on 2D and 3D, culture samples and cell viability was measured under perfusion flow (2.5 mL/min, 0.016 Pa). The visualization of cell response to their environment, in real time, will help to further elucidate the influences of biomaterial surface features, scaffold architectures, and the influence of flow induced shear on cell response and growth of new tissue.     

Epithelial Cell Repopulation and Preparation of Rodent Extracellular Matrix Scaffolds for Renal Tissue Development

This protocol details the generation of acellular, yet biofunctional, renal extracellular matrix (ECM) scaffolds that are useful as small-scale model substrates for organ-scale tissue development. Sprague Dawley rat kidneys are cannulated by inserting a catheter into the renal artery and perfused with a series of low-concentration detergents (Triton X-100 and sodium dodecyl sulfate (SDS)) over 26 hr to derive intact, whole-kidney scaffolds with intact perfusable vasculature, glomeruli, and renal tubules. Following decellularization, the renal scaffold is placed inside a custom-designed perfusion bioreactor vessel, and the catheterized renal artery is connected to a perfusion circuit consisting of: a peristaltic pump; tubing; and optional probes for pH, dissolved oxygen, and pressure. After sterilizing the scaffold with peracetic acid and ethanol, and balancing the pH (7.4), the kidney scaffold is prepared for seeding via perfusion of culture medium within a large-capacity incubator maintained at 37 °C and 5% CO₂. Forty million renal cortical tubular epithelial (RCTE) cells are injected through the renal artery, and rapidly perfused through the scaffold under high flow (25 ml/min) and pressure (~230 mmHg) for 15 min before reducing the flow to a physiological rate (4 ml/min). RCTE cells primarily populate the tubular ECM niche within the renal cortex, proliferate, and form tubular epithelial structures over seven days of perfusion culture. A 44 µM resazurin solution in culture medium is perfused through the kidney for 1 hr during medium exchanges to provide a fluorometric, redox-based metabolic assessment of cell viability and proliferation during tubulogenesis. The kidney perfusion bioreactor permits non-invasive sampling of medium for biochemical assessment, and multiple inlet ports allow alternative retrograde seeding through the renal vein or ureter. These protocols can be used to recellularize kidney scaffolds with a variety of cell types, including vascular endothelial, tubular epithelial, and stromal fibroblasts, for rapid evaluation within this system.     

Enhancement of cell growth in tissue-engineering constructs under direct perfusion: Modeling and simulation

Perfusion bioreactors improve mass transfer in cell-scaffold constructs. We developed a mathematical model to simulate nutrient flow through cellular constructs. Interactions among cell proliferation, nutrient consumption, and culture medium circulation were investigated. The model incorporated modified Contois cell-growth kinetics that includes effects of nutrient saturation and limited cell growth. Nutrient uptake was depicted through the Michaelis-Menton kinetics. To describe the culture medium convection, the fluid flow outside the cell-scaffold construct was described by the Navier-Stokes equations, while the fluid dynamics within the construct was modeled by Brinkman's equation for porous media flow. Effects of the media perfusion were examined by including time-dependant porosity and permeability changes due to cell growth. The overall cell volume was considered to consist of cells and extracellular matrices (ECM) as a whole without treating ECM separately. Numerical simulations show when cells were cultured subjected to direct perfusion, they penetrated to a greater extent into the scaffold and resulted in a more uniform spatial distribution. The cell amount was increased by perfusion and ultimately approached an asymptotic value as the perfusion rates increased in terms of the dimensionless Peclet number that accounts for the ratio of nutrient perfusion to diffusion. In addition to enhancing the nutrient delivery, perfusion simultaneously imposes flow-mediated shear stress to the engineered cells. Shear stresses were found to increase with cell growth as the scaffold void space was occupied by the cell and ECM volumes. The macro average stresses increased from 0.2 mPa to 1 mPa at a perfusion rate of 20 microm/s with the overall cell volume fraction growing from 0.4 to 0.7, which made the overall permeability value decrease from 1.35 x 10(-2)cm(2) to 5.51 x 10(-4)cm(2). Relating the simulation results with perfusion experiments in literature, the average shear stresses were below the critical value that would induce the chondrocyte necrosis.     

The effect of hydrodynamic shear on 3D engineered chondrocyte systems subject to direct perfusion

Bioreactors allowing direct-perfusion of culture medium through tissue-engineered constructs may overcome diffusion limitations associated with static culturing, and may provide flow-mediated mechanical stimuli. The hydrodynamic stress imposed on cells within scaffolds is directly dependent on scaffold microstructure and on bioreactor configuration. Aim of this study is to investigate optimal shear stress ranges and to quantitatively predict the levels of hydrodynamic shear imposed to cells during the experiments. Bovine articular chondrocytes were seeded on polyestherurethane foams and cultured for 2 weeks in a direct perfusion bioreactor designed to impose 4 different values of shear level at a single flow rate (0.5 ml/min). Computational fluid dynamics (CFD) simulations were carried out on reconstructions of the scaffold obtained from micro-computed tomography images. Biochemistry analyses for DNA and sGAG were performed, along with electron microscopy. The hydrodynamic shear induced on cells within constructs, as estimated by CFD simulations, ranged from 4.6 to 56 mPa. This 12-fold increase in the level of applied shear stress determined a 1.7-fold increase in the mean content in DNA and a 2.9-fold increase in the mean content in sGAG. In contrast, the mean sGAG/DNA ratio showed a tendency to decrease for increasing shear levels. Our results suggest that the optimal condition to favour sGAG synthesis in engineered constructs, at least at the beginning of culture, is direct perfusion at the lowest level of hydrodynamic shear. In conclusion, the presented results represent a first attempt to quantitatively correlate the imposed hydrodynamic shear level and the invoked biosynthetic response in 3D engineered chondrocyte systems.     

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.     

Influence of flow rate and scaffold pore size on cell behavior during mechanical stimulation in a flow perfusion bioreactor

Mechanically stimulating cell-seeded scaffolds by flow-perfusion is one approach utilized for developing clinically applicable bone graft substitutes. A key challenge is determining the magnitude of stimuli to apply that enhances cell differentiation but minimizes cell detachment from the scaffold. In this study, we employed a combined computational modeling and experimental approach to examine how the scaffold mean pore size influences cell attachment morphology and subsequently impacts upon cell deformation and detachment when subjected to fluid-flow. Cell detachment from osteoblast-seeded collagen-GAG scaffolds was evaluated experimentally across a range of scaffold pore sizes subjected to different flow rates and exposure times in a perfusion bioreactor. Cell detachment was found to be proportional to flow rate and inversely proportional to pore size. Using this data, a theoretical model was derived that accurately predicted cell detachment as a function of mean shear stress, mean pore size, and time. Computational modeling of cell deformation in response to fluid flow showed the percentage of cells exceeding a critical threshold of deformation correlated with cell detachment experimentally and the majority of these cells were of a bridging morphology (cells stretched across pores). These findings will help researchers optimize the mean pore size of scaffolds and perfusion bioreactor operating conditions to manage cell detachment when mechanically simulating cells via flow perfusion.     

Computational evaluation of oxygen and shear stress distributions in 3D perfusion culture systems: macro-scale and micro-structured models

We present a combined macro-scale/micro-scale computational approach to quantify oxygen transport and flow-mediated shear stress to human chondrocytes cultured in three-dimensional scaffolds in a perfusion bioreactor system. A macro-scale model was developed to assess the influence of the bioreactor design and to identify the proper boundary conditions for the micro-scale model. The micro-scale model based on a micro-computed tomography (microCT) reconstruction of a poly(ethylene glycol terephthalate)/poly(butylene terephthalate) (PEGT/PBT) foam scaffold, was developed to assess the influence of the scaffold micro-architecture on local shear stress and oxygen levels within the scaffold pores. Experiments were performed to derive specific oxygen consumption rates for constructs perfused under flow rates of 0.3 and 0.03 ml min(-1). While macro-scale and micro-scale models predicted similar average oxygen levels at different depths within the scaffold, microCT models revealed small local oxygen variations within the scaffold micro-architecture. The combined macro-scale/micro-scale approach indicated that 0.3 ml min(-1), which subjected 95% of the cells to less than 6.3 mPa shear, would maintain the oxygen supply throughout the scaffold above anoxic levels (>1%), with 99.5% of the scaffold supplied with 8-2% O(2). Alternatively, at 0.03 ml min(-1), the macro-scale model predicted 6% of the cells would be supplied with 0.5-1% O(2), although this region of cells was confined to the periphery of the scaffold. Together with local variations predicted by the micro-scale model, the simulations underline that in the current model system, reducing the flow below 0.03 ml min(-1) would likely have dire consequences on cell viability to pronounced regions within the engineered construct. The presented approach provides a sensitive tool to aid efficient bioreactor optimization and scaffold design.     

Shear stress modulates tumour cell adhesion to the endothelium

The adhesion of breast adenocarcinoma cells (MDA-MB-231) to human umbilical vein endothelial cells (HUVEC) was studied in whole blood and under varying flow conditions. This study was done on HUVEC either kept under static conditions or pre-conditioned in flow for 2 hours at a shear stress of 5 or 13 dyn/cm(2). Coverslips coated by HUVEC were placed in a parallel plate perfusion chamber and perfused at a shear rate of 300 or 1500 sec(-1) with heparin-anticoagulated blood containing 111In labelled MDA-MB-231 cells. We report here the optimal conditions for studying the adhesion of MDA-MB-231 to endothelial cells under shear constraints corresponding to those observed into small and medium sized arteries.     

Numerical simulation of sinusoidal flow in a straight elastic tube: effects of phase angles

Numerical simulations of flow in straight elastic (moving wall) tubes subjected to a sinusoidal pressure gradient were performed for conditions prevailing in large and medium sized arteries. The effects of varying the phase angle between the pressure gradient and the tube radius, the amplitude of wall motion, and the unsteadiness parameter (alpha) on flow rate and wall shear stress were investigated. Mean and peak flow rates and shear stresses were found to be strongly affected by the phase angle between the pressure gradient and the tube radius with greater sensitivity at higher diameter variation and higher alpha. In large artery simulations (alpha = 12), means flow rate was found to be 60% higher and peak flow rate to be 73% higher than corresponding rigid tube values for certain phase angles, while a threefold increase in mean wall shear stress and sevenfold increase in peak wall shear stress were observed in a sensitive phase angle range. Significant reversal in the wall shear stress direction occurred in the sensitive phase angle range even when there was negligible flow rate reversal. All effects were greatly diminished in simulations of medium sized vessels (alpha = 4). Some experimental evidence to support the predictions of a strong effect of phase angle on wall shear stress in large vessels is presented. Finally, physiological implications of the present work are discussed from a basis of aortic input impedance data, and a physical explanation for the extreme sensitivity of the flow field to small amplitude wall motion at high alpha is given.     

Theoretical analysis of insulin-dependent glucose uptake heterogeneity in 3D bioreactor cell culture

Three-dimensional (3D) cell cultures in bioreactors are becoming relevant as models for biological and physiological in vitro studies. In such systems, mathematical models can assist the experiment design that links the macroscopic properties to single-cell responses. We investigated the relationship between biochemical stimuli and cell response within a 3D cell culture in scaffold with heterogeneous porosity. Specifically, we studied the effect of insulin on the local glucose metabolism as a function of 3D pore size distribution. The multiscale mathematical model combines the mass transport within a 3D scaffold and a signaling pathways model. It considers the scaffold heterogeneity, and it describes spatiotemporal concentration of metabolites, biochemical stimuli, and cell density. The signaling model was integrated into this model, linking the local insulin concentration at cell membrane to the glucose uptake rate through glucose transporter type 4 (GLUT4) translocation from the cytosol to the cell membrane. The integrated model determines the cell response heterogeneities in a single channel, hence the biological response distribution in a 3D system. It also provides macroscopic outcomes to evaluate the feasibility of an experimental measurement of the system response. From our analysis, it became apparent that the flow rate is the most important operative variable, and that an optimum value ensures a fast and detectable cell response. This model on insulin-dependent glucose consumption rate offers insight into the cell metabolism physiology, which is a fundamental requirement for the study metabolic disorder such as Type 2 diabetes mellitus, in which the physiological insulin-dependent glucose metabolism is impaired.     

Effects of oxygen on engineered cardiac muscle

Concentration gradients associated with the in vitro cultivation of engineered tissues that are vascularized in vivo result in the formation of only a thin peripheral tissue-like region (e.g., approximately 100 microm for engineered cardiac muscle) around a relatively cell-free interior. We previously demonstrated that diffusional gradients within engineered cardiac constructs can be minimized by direct perfusion of culture medium through the construct. In the present study, we measured the effects of medium perfusion rate and local oxygen concentration (p(O2)) on the in vitro reconstruction of engineered cardiac muscle. Neonatal rat cardiomyocytes were seeded onto biodegradable polymer scaffolds (fibrous discs, 1.1 cm diameter x 2 mm thick, made of polyglycolic acid, 24 x 10(6) cells per scaffold). The resulting cell-polymer constructs were cultured for a total of 12 days in serially connected cartridges (n = 1-8), each containing one construct directly perfused with culture medium at a flow rate of 0.2-3.0 mL/min. In all groups, oxygen concentration decreased due to cell respiration, and depended on construct position in the series and medium flow rate. Higher perfusion rates and higher p(O2) correlated with more aerobic cell metabolism, and higher DNA and protein contents. Constructs cultured at p(O2) of 160 mm Hg had 50% higher DNA and protein contents, markedly higher expression of sarcomeric alpha-actin, better organized sarcomeres and cell junctions, and 4.5-fold higher rate of cell respiration as compared to constructs cultured at p(O2) of 60 mm Hg. Contraction rates of the corresponding cardiac cell monolayers were 40% higher at p(O2) of 160 than 60 mm Hg. The control of oxygen concentration in cell microenvironment can thus improve the structure and function of engineered cardiac muscle. Experiments of this kind can form a basis for controlled studies of the effects of oxygen on the in vitro development of engineered tissues.     

A mathematical model for fluid shear-sensitive 3D tissue construct development

This research studies dynamic culture for 3D tissue construct development with computational fluid dynamics. It proposes a mathematical model to evaluate the impact of flow rates and flow shear stress on cell growth in 3D constructs under perfusion. The modeling results show that dynamic flow, even at flow rate as low as 0.002 cm/s, can support much better mass exchange, higher cell number, and more even cell and nutrient distribution compared to static culture. Higher flow rate can further improve nutrient supply and mass exchange in the construct, promoting better nutritious environment and cell proliferation compared to lower flow rate. In addition, consideration of flow shear stress predicts much higher cell number in the construct compared to that without shear consideration. While the nutrient can dominate shear stress in influencing cell proliferation, the shear effect increases with flow rate. The proposed model helps tissue engineers better understand the cell-flow relationship at the molecular level during dynamic culture.     

Modeling fluid flow through irregular scaffolds for perfusion bioreactors

Direct perfusion of 3D tissue engineered constructs is known to enhance osteogenesis, which can be partly attributed to enhanced nutrient and waste transport. In addition flow mediated shear stresses are known to upregulate osteogenic differentiation and mineralization. A quantification of the hydrodynamic environment is therefore crucial to interpret and compare results of in vitro bioreactor experiments. This study aims to deal with the pitfalls of numerical model preparation of highly complex 3D bone scaffold structures and aims to provide more accurate wall shear stress (WSS) estimates. µCT imaging techniques were used to reconstruct the geometry of both a titanium (Ti) and a hydroxyapatite scaffold, starting from 430 images with a resolution of 8 µm. To tackle the tradeoff between model size and mesh resolution we selected two concentric regions of interest (cubes with a volume of 1 and 3.375 mm³, respectively) for both scaffolds. A flow guidance in front of the real inlet surface of the scaffold was designed to mimic realistic inlet conditions. With a flow rate of 0.04 mL/min perfused through a 5 mm diameter scaffold at an inlet velocity of 33.95 µm/s we obtained average WSSs of 1.10 and 1.46 mPa for the 1 mm³ and the 3.375 mm³ model of the hydroxyapatite scaffold compared to 1.40 and 1.95 mPa for the 1 mm³ model and the 3.375 mm³ model of the Ti scaffold, showing the important influence of the scaffold micro‐architecture heterogeneity and the proximity of boundaries. To assess that influence we selected cubic portions, of which the WSS data were analyzed, with the same size and the same location within both 1 and 3.375 mm³ cubic models. Varying the size of the inner portions simultaneously in both model selections gives a quantification of the sensitivity to boundary neighborhood. This methodology allows to get more insight in the complex concept of tissue engineering and will likely help to understand and eventually improve the fluid‐mechanical aspects. Biotechnol. Bioeng. 2009;103: 621–630. © 2009 Wiley Periodicals, Inc.