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.     

Computational modeling of combined cell population dynamics and oxygen transport in engineered tissue subject to interstitial perfusion

This work presents a computational model of tissue growth under interstitial perfusion inside a tissue engineering bioreactor. The model accounts both for the cell population dynamics, using a model based on cellular automata, and for the hydrodynamic microenvironment imposed by the bioreactor, using a model based on the Lattice-Boltzmann equation and the convection-diffusion equation. The conditions of static culture versus perfused culture were compared, by including the population dynamics along with oxygen diffusion, convective transport and consumption. The model is able to deal with arbitrary complex geometries of the spatial domain; in the present work, the domain modeled was the void space of a porous scaffold for tissue-engineered cartilage. The cell population dynamics algorithm provided results which qualitatively resembled population dynamics patterns observed in experimental studies, and these results were in good quantitative agreement with previous computational studies. Simulation of oxygen transport and consumption showed the fundamental contribution of convective transport in maintaining a high level of oxygen concentration in the whole spatial domain of the scaffold. The model was designed with the aim to be computationally efficient and easily expandable, i.e. to allow straightforward implementability of further models of complex biological phenomena of increasing scientific interest in tissue engineering, such as chemotaxis, extracellular matrix deposition and effect of mechanical stimulation.     

Computational modeling of combined cell population dynamics and oxygen transport in engineered tissue subject to interstitial perfusion

This work presents a computational model of tissue growth under interstitial perfusion inside a tissue engineering bioreactor. The model accounts both for the cell population dynamics, using a model based on cellular automata, and for the hydrodynamic microenvironment imposed by the bioreactor, using a model based on the Lattice–Boltzmann equation and the convection-diffusion equation. The conditions of static culture versus perfused culture were compared, by including the population dynamics along with oxygen diffusion, convective transport and consumption. The model is able to deal with arbitrary complex geometries of the spatial domain; in the present work, the domain modeled was the void space of a porous scaffold for tissue-engineered cartilage. The cell population dynamics algorithm provided results which qualitatively resembled population dynamics patterns observed in experimental studies, and these results were in good quantitative agreement with previous computational studies. Simulation of oxygen transport and consumption showed the fundamental contribution of convective transport in maintaining a high level of oxygen concentration in the whole spatial domain of the scaffold. The model was designed with the aim to be computationally efficient and easily expandable, i.e. to allow straightforward implementability of further models of complex biological phenomena of increasing scientific interest in tissue engineering, such as chemotaxis, extracellular matrix deposition and effect of mechanical stimulation.     

Effects of shear stress on articular chondrocyte metabolism

The articular cartilage of diarthrodial joints experiences a variety of stresses, strains and pressures that result from normal activities of daily living. In normal cartilage, the extracellular matrix exists as a highly organized composite of specialized macromolecules that distributes loads at the bony ends. The chondrocyte response to mechanical loading is recognized as an integral component in the maintenance of articular cartilage matrix homeostasis. With inappropriate mechanical loading of the joint, as occurs with traumatic injury, ligament instability, bony malalignment or excessive weight bearing, the cartilage exhibits manifestations characteristic of osteoarthritis. Breakdown of cartilage in osteoarthritis involves degradation of the extracellular matrix macromolecules and decreased expression of chondrocyte proteins necessary for normal joint function. Osteoarthritic cartilage often exhibits increased amounts of type I collagen and synthesis of proteoglycans characteristic of immature cartilage. The shift in cartilage phenotype in response to altered load yields a matrix that fails to support normal joint function. Mathematical modeling and experimental studies in animal models confirm an association between altered loading of diarthrotic joints and arthritic changes. Both types of studies implicate shear forces as a critical component in the destructive profile. The severity of cartilage destruction in response to altered loads appears linked to expression of biological factors influencing matrix integrity and cellular metabolism. Determining how shear stress alters chondrocyte metabolism is fundamental to understanding how to limit matrix destruction and stimulate cartilage repair and regeneration. At present, the precise biochemical and molecular mechanisms by which shear forces alter chondrocyte metabolism from a normal to a degenerative phenotype remain unclear. The results presented here address the hypothesis that articular chondrocyte metabolism is modulated by direct effects of shear forces that act on the cell through mechanotransduction processes. The purpose of this work is to develop critical knowledge regarding the basic mechanisms by which mechanical loading modulates cartilage metabolism in health and disease. This presentation will describe the effects of using fluid induced shear stress as a model system for stimulation of articular chondrocytes in vitro. The fluid induced shear stress was applied using a cone viscometer system to stimulate all the cells uniformly under conditions of minimal turbulence. The experiments were carried using high-density primary monolayer cultures of normal and osteoarthritic human and normal bovine articular chondrocytes. The analysis of the cellular response included quantification of cytokine release, matrix metalloproteinase expression and activation of intracellular signaling pathways. The data presented here show that articular chondrocytes exhibit a dose- and time-dependent response to shear stress that results in the release of soluble mediators and extracellular matrix macromolecules. The data suggest that the chondrocyte response to mechanical stimulation contributes to the maintenance of articular cartilage homeostasis in vivo.     

The effect of solid surface tension and exposure to elevated hydrodynamic shear on Pseudomonas fluorescens biofilms grown on modified titanium surfaces

The solid surface tension of titanium was varied by using organosilane monolayers of various terminations, minimising differences in other material properties. Both the quantity of Pseudomonas fluorescens biofilms grown on the modified surfaces, and the percentage of biofilm remaining after exposure to hydrodynamic shear stress, varied significantly as a function of solid surface tension. The quantity of biofilm was less on chloropropyl-terminated surfaces than on an alkyl-terminated surfaces. However, the percentage of biofilm remaining after exposure to hydrodynamic shear stress (which depends on the adhesion and cohesion strengths of the biofilm) was less for the alkyl-terminated surface than for the chloropropyl-terminated surface, for one of the two sample sets analysed. These results demonstrate the importance of differentiating between the quantity of biofilm on a surface and the adhesion and cohesion strength of the biofilm, and may help explain discrepancies in the existing literature regarding the effect of solid surface tension on the propensity of a surface for microfouling.     

The effect of solid surface tension and exposure to elevated hydrodynamic shear on Pseudomonas fluorescens biofilms grown on modified titanium surfaces

Abstract

The solid surface tension of titanium was varied by using organosilane monolayers of various terminations, minimising differences in other material properties. Both the quantity of Pseudomonas fluorescens biofilms grown on the modified surfaces, and the percentage of biofilm remaining after exposure to hydrodynamic shear stress, varied significantly as a function of solid surface tension. The quantity of biofilm was less on chloropropyl-terminated surfaces than on an alkyl-terminated surfaces. However, the percentage of biofilm remaining after exposure to hydrodynamic shear stress (which depends on the adhesion and cohesion strengths of the biofilm) was less for the alkyl-terminated surface than for the chloropropyl-terminated surface, for one of the two sample sets analysed. These results demonstrate the importance of differentiating between the quantity of biofilm on a surface and the adhesion and cohesion strength of the biofilm, and may help explain discrepancies in the existing literature regarding the effect of solid surface tension on the propensity of a surface for microfouling.     

Effects of shear stress on 3-D human mesenchymal stem cell construct development in a perfusion bioreactor system: Experiments and hydrodynamic modeling

Shear stress is an important biomechanical parameter in regulating human mesenchymal stem cell (hMSC) construct development. In this study, the biomechanical characteristics of hMSCs within highly porous 3-D poly (ethylene terephthalate) (PET) matrices in a perfusion bioreactor system were analyzed for two flow rates of 0.1 and 1.5 mL/min, respectively over a 20-day culture period. A 1.4 times higher proliferation rate, higher CFU-F formation, and more fibronectin and HSP-47 secretion at day 20 were observed at the flow rate of 0.1 mL/min compared to those at the flow rate of 1.5 mL/min. The higher flow rate of 1.5 mL/min upregulated osteogenic differentiation potential at day 20 as measured by the expression of alkaline phosphatase activity and calcium deposition in the matrix after 14 days osteogenic induction, consistent with those reported in literatures. Mathematical modeling indicated that shear stress existed in the range of 1 x 10(-5) to 1 x 10(-4) Pa in the constructs up to a depth of 70 microm due to flow penetration in the porous constructs. Analysis of oxygen transport in the constructs for the two flow rates yielded oxygen levels significantly higher than those at which cell growth and metabolism are affected (Jiang et al., 1996). This indicates that differences in convective transport have no significant influence on cell growth and metabolism for the range of flow rates studied. These results demonstrate that shear stress is an important microenvironment parameter that regulates hMSC construct development at a range significantly lower than those reported previously in the perfusion system.     

Uniform tissues engineered by seeding and culturing cells in 3D scaffolds under perfusion at defined oxygen tensions

In this work, we assessed whether culture of uniformly seeded chondrocytes under direct perfusion, which supplies the cells with normoxic oxygen levels, can maintain a uniform distribution of viable cells throughout porous scaffolds several milimeters in thickness, and support the development of uniform tissue grafts. An integrated bioreactor system was first developed to streamline the steps of perfusion cell seeding of porous scaffolds and perfusion culture of the cell-seeded scaffolds. Oxygen tensions in perfused constructs were monitored by in-line oxygen sensors incorporated at the construct inlet and outlet. Adult human articular chondrocytes were perfusion-seeded into 4.5 mm thick foam scaffolds at a rate of 1 mm/s. Cell-seeded foams were then either cultured statically in dishes or further cultured under perfusion at a rate of 100 microm/s for 2 weeks. Following perfusion seeding, viable cells were uniformly distributed throughout the foams. Constructs subsequently cultured statically were highly heterogeneous, with cells and matrix concentrated at the construct periphery. In contrast, constructs cultured under perfusion were highly homogeneous, with uniform distributions of cells and matrix. Oxygen tensions of the perfused medium were maintained near normoxic levels (inlet congruent with 20%, outlet > 15%) at all times of culture. We have demonstrated that perfusion culture of cells seeded uniformly within porous scaffolds, at a flow rate maintaining a homogeneous oxygen supply, supports the development of uniform engineering tissue grafts of clinically relevant thicknesses.     

The stability of predator-prey systems subject to the Allee effects

In recent years, many theoreticians and experimentalists have concentrated on the processes that affect the stability of predator-prey systems. But few papers have addressed the Allee effect with focus on the their stability. In this paper, we select two classical models describing predator-prey systems and introduce the Allee effects into the dynamics of both the predator and prey populations in these models, respectively. By combining mathematical analysis with numerical simulation, we have shown that the Allee effect may be a destabilizing force in predator-prey systems: the equilibrium point of the system could be changed from stable to unstable or otherwise, the system, even when it is stable, will take much longer time to reach the stable state. We also conclude that the equilibrium of the prey population will be enlarged due to the Allee effect of the predator, but the Allee effects of the prey may decrease the equilibrium value of the predator, or that of both the predator and prey. It should also be pointed out that the impact of the Allee effects of predator and prey due to different mechanisms on different predator-prey systems could also vary.     

IL-1beta synthesis by chondrocyte analyzed by 3D microscopy and flow cytometry: effect of Rhein

Several factors are known to be involved in the destruction of the articular cartilage. Interleukin-1 (IL-1) plays an important role in the pathogenesis of osteoarthritis (OA) either directly or through the stimulation of catabolic factors. The action of IL-1 on articular cartilage is multifaceted and it most likely plays an important role in the mechanism of cartilage destruction. IL-1 suppresses the synthesis of the cartilage matrix components and promotes the degradation of cartilage matrix macromolecules. Diacerein is an anthraquinone molecule that has been shown to reduce the severity of OA, both in man and in animal models. The present study was designed to evaluate in vitro effects of diacerein on IL-1beta expression in LPS or IL-1alpha stimulated chondrocytes. Intracellular IL-1beta production was analysed in articular chondrocytes cultured in monolayer or in alginate 3D-biosystems in the presence of lipopolysaccharide (LPS) or IL-1alpha, with or without diacerein. The results show that LPS and IL-1alpha increase intracellular IL-1beta and Diacerein inhibited LPS-induced and IL-1alpha induced IL-1beta production by articular chondrocytes. Moreover, the effect of mechanical stimulation was analysed. An inhibitory effect of DAR at therapeutic concentrations on IL-1beta production in articular chondrocytes is suggested.     

Engineered cartilage constructs subject to very low regimens of interstitial perfusion

We have studied an in vitro engineered cartilage model, consisting of bovine articular chondrocytes seeded on micro-porous scaffolds and perfused with very low regimens of interstitial flow. Our previous findings suggested that synthesis of sulphated glycosaminoglycans (sGAG) was promoted in this model, if the level of shear generated on cells was maintained below 10 mPa (0.1 dyn/cm2). Constructs were stimulated with a median shear stress of 1.2 and 6.7 mPa using two independent culture chambers. Quantification of the applied stresses and of oxygen consumption rates was obtained from computational modelling. Experimentally, we set a time zero reference at 24 hours after cell seeding and total culture time at two weeks. The cell metabolic activity, measured by MTT, was significantly lower in all constructs at two weeks (-73% in static controls, -66% in the 1.2 mPa group and -60% in the 6.7 mPa group) vs. the time zero group, and significantly higher (+33%) in the 7 mPa group vs. static controls. The ratio between synthesis of collagen type II/type I, measured by Western Blot, was significantly higher in the 1.2 mPa constructs (+109% vs. the 6.7 mPa group, +120% vs. the time zero group and +286% vs. static controls). A trend of decreased alpha-actin expression was observed with increased ratio of type II to type I collagen, in all groups. These results reinforce the notion that, at early time points in culture, hydrodynamic shear below 10 mPa may promote formation of extra-cellular matrix specific to hyaline cartilage in chondrocyte-seeded constructs.     

The competitive dynamics of metapopulations subject to the Allee-like effect

It is well recognized that individuals of many species can benefit from the presence of conspecifics, a concept broadly referred to as the Allee effect. At the metapopulation level, there is an analogous but essentially different phenomenon called the Allee-like effect that leads to metapopulation extinction thresholds at low habitat occupancy. But so far not adequate attention has been paid to this phenomenon. In this paper, the Allee-like effect is introduced into a metapopulation model of one species and also that of a three-state two-species competitive system. Phase plane analysis is used to investigate the dynamics of these models. We demonstrate that the Allee-like effect alone could lead to multiple stable states in three-state two-species competitive systems at the metapopulation level, and the number of stable states decrease as the Allee-like effect becomes more severe. Severe Allee-like effects may make coexistence impossible and may even lead to the extinction of both species even if their initial habitat occupancies are high and suitable habitats are enough. It is especially noticeable that depending on their initial conditions one species may exclude the other one that subjects to a weaker Allee-like effect than the former, while the second species always excludes the first one when both species are assumed to be in the absence of the Allee-like effect. We also investigate the habitat destructive effect on the Allee-like system mentioned above. Research indicates that the existence of the Allee-like effect makes a metapopulation more susceptible to habitat destruction. All in all, the Allee-like effect is probably a destabilizing factor that, together with habitat destruction, would affect the continuous existence of species. These conclusions may have important implications for conservation and metacommunity organization.     

Concentric cylinder bioreactor for production of tissue engineered cartilage: effect of seeding density and hydrodynamic loading on construct development

A concentric cylinder bioreactor has been developed to culture tissue engineered cartilage constructs under hydrodynamic loading. This bioreactor operates in a low shear stress environment, has a large growth area for construct production, allows for dynamic seeding of constructs, and provides for a uniform loading environment. Porous poly-lactic acid constructs, seeded dynamically in the bioreactor using isolated bovine chondrocytes, were cultured for 4 weeks at three seeding densities (60, 80, 100 x 10(6) cells per bioreactor) and three different shear stresses (imposed at 19, 38, and 76 rpm) to characterize the effect of chondrocyte density and hydrodynamic loading on construct growth. Construct seeding efficiency with chondrocytes is greater than 95% within 24 h. Extensive chondrocyte proliferation and matrix deposition are achieved so that after 28 days in culture, constructs from bioreactors seeded at the highest cell densities contain up to 15 x 10(6) cells, 2 mg GAG, and 3.5 mg collagen per construct and exhibit morphology similar to that of native cartilage. Bioreactors seeded with 60 million chondrocytes do not exhibit robust proliferation or matrix deposition and do not achieve morphology similar to that of native cartilage. In cultures under different steady hydrodynamic loading, the data demonstrate that higher shear stress suppresses matrix GAG deposition and encourages collagen incorporation. In contrast, under dynamic hydrodynamic loading conditions, cartilage constructs exhibit robust matrix collagen and GAG deposition. The data demonstrate that the concentric cylinder bioreactor provides a favorable hydrodynamic environment for cartilage construct growth and differentiation. Notably, construct matrix accumulation can be manipulated by hydrodynamic loading. This bioreactor is useful for fundamental studies of construct growth and to assess the significance of cell density, nutrients, and hydrodynamic loading on cartilage development. In addition, studies of cartilage tissue engineering in the well-characterized, uniform environment of the concentric cylinder bioreactor will develop important knowledge of bioprocessing parameters critical for large-scale production of engineered tissues.     

Anti-biofouling 3D porous systems: the blend effect of oxazoline-based oligomers on chitosan scaffolds

The production, characterization and anti-biofouling activity of 3D porous scaffolds combining different blends of chitosan and oxazoline-based antimicrobial oligomers is reported. The incorporation of ammonium quaternized oligo(2-oxazoline)s into the composition of the scaffold enhances the stability of the chitosan scaffold under physiological conditions as well as its ability to repel protein adsorption. The blended scaffolds showed mean pore sizes in the range of 18–32?μm, a good pore interconnectivity and high porosity, as well as a large surface area, ultimate key features for anti-biofouling applications. Bovine serum albumin (BSA) adhesion profiles showed that the composition of the scaffolds plays a critical role in the chitosan–oligooxazoline system. Oligobisoxazoline-enriched scaffolds (20%?w/w, CB8020) decreased protein adsorption (BSA) by up to 70%. Moreover, 1?mg of CB8020 was able to kill 99.9% of Escherichia coli cells upon contact, demonstrating its potential as promising material for production of tailored non-fouling 3D structures to be used in the construction of novel devices with applications in the biomedical field and water treatment processes.     

The effect of implantation on scaffoldless three-dimensional engineered bone constructs

Our laboratory has previously developed scaffoldless engineered bone constructs (EBC). Bone marrow stromal cells (BMSC) were harvested from rat femur and cultured in medium that induced osteogenic differentiation. After reaching confluence, the monolayer of cells contracted around two constraint points forming a cylinder. EBCs were placed in small diameter (0.5905 × 0.0625 in.) or large diameter (0.5905 × 0.125 in.) silicone tubing and implanted intramuscularly in the hind limb of a rat. Bone mineral content (BMC) of the EBC was analyzed before implantation and at 1 and 2 mo following implantation and compared to that of native femur bone at different stages of development. Negligible BMC was observed in E-20 femur or EBCs prior to implantation. One-month implantation in both small and large tubing increased BMC in the EBC. BMC of EBC from large tubing was greater than in 14 d rat neonatal femurs, but was 2% and 3% of BMC content in adult bone after 1 and 2 mo of implantation, respectively. Alizarine Red and osteopontin staining of the EBCs before and after implantation confirmed increased bone mineralization in the implanted EBCs. Implanted EBCs also had extensive vascularization. Our data suggest that BMSC can be successfully used for the generation of scaffoldless EBC, and this model can be potentially used for the generation of autologous bone transplants in humans.