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.     

Heterogeneous perfusion is a consequence of uniform shear stress in optimized arterial tree models

Using optimized computer models of arterial trees we demonstrate that flow heterogeneity is a necessary consequence of a uniform shear stress distribution. Model trees are generated and optimized under different modes of boundary conditions. In one mode flow is delivered to the tissue as homogeneously as possible. Although this primary goal can be achieved, resulting shear stresses between blood and the vessel walls show very large spread. In a second mode, models are optimized under the condition of uniform shear stress in all segments which in turn renders flow distribution heterogeneous. Both homogeneous perfusion and uniform shear stress are desirable goals in real arterial trees but each of these goals can only be approached at the expense of the other. While the present paper refers only to optimized models, we assume that this dual relation between the heterogeneities in flow and shear stress may represent a more general principle of vascular systems.     

Computational models for wall shear stress estimation in scaffolds: a comparative study of two complete geometries

Fluid mechanical stimuli are known to upregulate cell differentiation and matrix formation. Since wall shear stress plays an important role various studies tried to estimate the scaffold fluid dynamic environment. However, because of the geometrical complexity, nearly all studies created their CFD model based on a submodel of the entire scaffold assuming that the model covers heterogeneity sufficiently. However to the authors' knowledge no study exist providing guidelines in this matter. In a previous study we demonstrated that submodels are influenced by the boundary conditions, inevitable when flow channels are chopped off. For the current study we therefore developed μCT based models of two complete scaffold geometries (one titanium and one hydroxyapatite). Imposing a 0.04 ml/min flow rate resulted in a surface area averaged wall shear stress of 1.41 mPa for titanium and 1.09 mPa for hydroxyapatite. In order to get insight in required model size we subdivided the domain in regions of different size. From our results we propose a model size between 6 and 10 times the average pore size. The wall shears stress should be calculated on a region at least one pore size away from the boundaries. These guidelines could be of use for computationally more costly simulations where it is not possible to simulate the complete scaffold domain.     

Three-dimensional macro-scale assessment of regional and temporal wall shear stress characteristics on aortic valve leaflets

The aortic valve (AV) achieves unidirectional blood flow between the left ventricle and the aorta. Although hemodynamic stresses have been shown to regulate valvular biology, the native wall shear stress (WSS) experienced by AV leaflets remains largely unknown. The objective of this study was to quantify computationally the macro-scale leaflet WSS environment using fluid–structure interaction modeling. An arbitrary Lagrangian–Eulerian approach was implemented to predict valvular flow and leaflet dynamics in a three-dimensional AV geometry subjected to physiologic transvalvular pressure. Local WSS characteristics were quantified in terms of temporal shear magnitude (TSM), oscillatory shear index (OSI) and temporal shear gradient (TSG). The dominant radial WSS predicted on the leaflets exhibited high amplitude and unidirectionality on the ventricularis (TSM>7.50 dyn/cm², OSI < 0.17, TSG>325.54 dyn/cm² s) but low amplitude and bidirectionality on the fibrosa (TSM < 2.73 dyn/cm², OSI>0.38, TSG < 191.17 dyn/cm² s). The radial WSS component computed in the leaflet base, belly and tip demonstrated strong regional variability (ventricularis TSM: 7.50–22.32 dyn/cm², fibrosa TSM: 1.26–2.73 dyn/cm²). While the circumferential WSS exhibited similar spatially dependent magnitude (ventricularis TSM: 1.41–3.40 dyn/cm², fibrosa TSM: 0.42–0.76 dyn/cm²) and side-specific amplitude (ventricularis TSG: 101.73–184.43 dyn/cm² s, fibrosa TSG: 41.92–54.10 dyn/cm² s), its temporal variations were consistently bidirectional (OSI>0.25). This study provides new insights into the role played by leaflet–blood flow interactions in valvular function and critical hemodynamic stress data for the assessment of the hemodynamic theory of AV disease.     

Quantifying effects of plaque structure and material properties on stress distributions in human atherosclerotic plaques using 3D FSI models

BACKGROUND: Atherosclerotic plaques may rupture without warning and cause acute cardiovascular syndromes such as heart attack and stroke. Methods to assess plaque vulnerability noninvasively and predict possible plaque rupture are urgently needed. METHOD: MRI-based three-dimensional unsteady models for human atherosclerotic plaques with multi-component plaque structure and fluid-structure interactions are introduced to perform mechanical analysis for human atherosclerotic plaques. RESULTS: Stress variations on critical sites such as a thin cap in the plaque can be 300% higher than that at other normal sites. Large calcification block considerably changes stress/strain distributions. Stiffness variations of plaque components (50% reduction or 100% increase) may affect maximal stress values by 20-50%. Plaque cap erosion causes almost no change on maximal stress level at the cap, but leads to 50% increase in maximal strain value. CONCLUSIONS: Effects caused by atherosclerotic plaque structure, cap thickness and erosion, material properties, and pulsating pressure conditions on stress/strain distributions in the plaque are quantified by extensive computational case studies and parameter evaluations. Computational mechanical analysis has good potential to improve accuracy of plaque vulnerability assessment.     

Microfluidic 3D cell culture: from tools to tissue models

1. Download : Download high-res image (132KB)
2. Download : Download full-size image

     

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.     

Computational models in systems biology: standards,dissemination, and best practices

1. Download : Download high-res image (94KB)
2. Download : Download full-size image

     

3D reconstruction techniques of human coronary bifurcations for shear stress computations

Background

Heterogeneity in plaque composition in human coronary artery bifurcations is associated with blood flow induced shear stress. Shear stress is generally determined by combing 3D lumen data and computational fluid dynamics (CFD). We investigated two new procedures to generate 3D lumen reconstructions of coronary artery bifurcations for shear stress computations.

Methods

We imaged 10 patients with multislice computer tomography (MSCT) and intravascular ultrasound (IVUS). The 3D reconstruction of the main branch was based on the fusion of MSCT and IVUS. The proximal part of side branch was reconstructed using IVUS data or MSCT data, resulting in two different reconstructions of the bifurcation region. The distal part of the side branch was based on MSCT data alone. The reconstructed lumen was combined with CFD to determine the shear stress. Low and high shear stress regions were defined and shear stress patterns in the bifurcation regions were investigated.

Results

The 3D coronary bifurcations were successfully generated with both reconstruction procedures. The geometrical features of the bifurcation region for the two reconstruction procedures did not reveal appreciable differences. The shear stress maps showed a qualitative agreement, and the low and high shear stress regions were similar in size and average shear stress values were identical. The low and high shear stress regions showed an overlap of approximately 75%.

Conclusion

Reconstruction of the side branch with MSCT data alone is an adequate technique to study shear stress and wall thickness in the bifurcation region. The reconstruction procedure can be applied to further investigate the effect of shear stress on atherosclerosis in coronary bifurcations.     

Axial stent strut angle influences wall shear stress after stent implantation: analysis using 3D computational fluid dynamics models of stent foreshortening

Introduction  

The success of vascular stents in the restoration of blood flow is limited by restenosis. Recent data generated from computational fluid dynamics (CFD) models suggest that the vascular geometry created by an implanted stent causes local alterations in wall shear stress (WSS) that are associated with neointimal hyperplasia (NH). Foreshortening is a potential limitation of stent design that may affect stent performance and the rate of restenosis. The angle created between axially aligned stent struts and the principal direction of blood flow varies with the degree to which the stent foreshortens after implantation.     

Computational fluid modeling and performance analysis of a bidirectional rotating perfusion culture system

A myriad of bioreactor configurations have been investigated as extracorporeal medical support systems for temporary replacement of vital organ functions. In recent years, studies have demonstrated that the rotating bioreactors have the potential to be utilized as bioartificial liver assist devices (BLADs) owing to their advantage of ease of scalability of cell‐culture volume. However, the fluid movement in the rotating chamber will expose the suspended cells to unwanted flow structures with abnormally high shear conditions that may result in poor cell stability and in turn lower the efficacy of the bioreactor system. In this study, we compared the hydrodynamic performance of our modified rotating bioreactor design with that of an existing rotating bioreactor design. Computational fluid dynamic analysis coupled with experimental results were employed in the optimization process for the development of the modified bioreactor design. Our simulation results showed that the modified bioreactor had lower fluid induced shear stresses and more uniform flow conditions within its rotating chamber than the conventional design. Experimental results revealed that the cells within the modified bioreactor also exhibited better cell‐carrier attachment, higher metabolic activity, and cell viability compared to those in the conventional design. In conclusion, this study was able to provide important insights into the flow physics within the rotating bioreactors, and help enhanced the hydrodynamic performance of an existing rotating bioreactor for BLAD applications. © 2013 American Institute of Chemical Engineers Biotechnol. Prog., 29:1002–1012, 2013     

Metabolic response of biofilm to shear stress in fixed-film culture

AIMS: In a biofilm reactor, detachment force resulting from hydraulic shear is a major factor that determines the formation and structure of steady state biofilm. The metabolic response of biofilm to change in shear stress was therefore investigated. METHODS AND RESULTS: A conventional annular reactor made of PVC was used, in which shearing over the rotating disc surface was strictly defined. Results from the steady state aerobic biofilm reactor showed that the biofilm structure (density and thickness) and metabolic behaviour (growth yield and dehydrogenase activity) were closely related to the shear stress exerted on the biofilm. Smooth, dense and stable biofilm formed at relatively high shear stress. Higher dehydrogenase activity and lower growth yield were obtained when the shear stress was raised. Growth yield was inversely correlated with the catabolic activity of biofilm. The reduced growth yield, together with the enhanced catabolic activity, suggests that a dissociation of catabolism from anabolism would occur at high shear stress. CONCLUSION: Biofilms may respond to shear stress by regulating metabolic pathways associated with the substrate flux flowing between catabolism and anabolism. A biological phenomenon, besides a simple physical effect, is underlying the observed relation between the shear stress and resulting biofilm structure. SIGNIFICANCE AND IMPACT OF THE STUDY: A hypothesis is proposed that the shear-induced energy spilling would be associated with a stimulated proton translocation across the cell membrane, which favours formation of a stronger biofilm. This research may provide a basis for experimental data on biofilm obtained at different shear stresses to be interpreted in relation to energy.     

Circumferential vascular deformation after stent implantation alters wall shear stress evaluated with time-dependent 3D computational fluid dynamics models.

The success of vascular stents in the restoration of blood flow is limited by restenosis. Recent data generated from computational fluid dynamics (CFD) models suggest that stent geometry may cause local alterations in wall shear stress (WSS) that have been associated with neointimal hyperplasia and subsequent restenosis. However, previous CFD studies have ignored histological evidence of vascular straightening between circumferential stent struts. We tested the hypothesis that consideration of stent-induced vascular deformation may more accurately predict alterations in indexes of WSS that may subsequently account for histological findings after stenting. We further tested the hypothesis that the severity of these alterations in WSS varies with the degree of vascular deformation after implantation. Steady-state and time-dependent simulations of three-dimensional CFD arteries based on canine coronary artery measurements of diameter and blood flow were conducted, and WSS and WSS gradients were calculated. Circumferential straightening introduced areas of high WSS between stent struts that were absent in stented vessels of circular cross section. The area of vessel exposed to low WSS was dependent on the degree of circumferential vascular deformation and axial location within the stent. Stents with four vs. eight struts increased the intrastrut area of low WSS in vessels, regardless of cross-sectional geometry. Elevated WSS gradients were also observed between struts in vessels with polygonal cross sections. The results obtained using three-dimensional CFD models suggest that changes in vascular geometry after stent implantation are important determinants of WSS distributions that may be associated with subsequent neointimal hyperplasia.     

Accuracy and reproducibility of CFD predicted wall shear stress using 3D ultrasound images

Computational fluid dynamics (CFD) flow simulation techniques have the potential to enhance our understanding of how haemodynamic factors are involved in atherosclerosis. Recently, 3D ultrasound has emerged as an alternative to other 3D imaging techniques, such as magnetic resonance angiography (MRA). The method can be used to generate realistic vascular geometry suitable for CFD simulations. In order to assess accuracy and reproducibility of the procedure from image acquisition to reconstruction to CFD simulation, a human carotid artery bifurcation phantom was scanned three times using 3D ultrasound. The geometry was reconstructed and flow simulations were carried out on the three sets as well as on a model generated using computer aided design (CAD) from the geometric information given by the manufacturer. It was found that the three reconstructed sets showed good reproducibility as well as satisfactory quantitative agreement with the CAD model. Analyzing two selected locations probably representing the 'worst cases,' accuracy comparing ultrasound and CAD reconstructed models was estimated to be between 7.2% and 7.7% of the maximum instantaneous WSS and reproducibility comparing the three scans to be between 8.2% and 10.7% of their average maximum.     

Development of 3D hydrogel culture systems with on-demand cell separation

Recently there has been an increased interest in the effects of paracrine signaling between groups of cells, particularly in the context of better understanding how stem cells contribute to tissue repair. Most current 3D co-culture methods lack the ability to effectively separate two cell populations after the culture period, which is important for simultaneously analyzing the reciprocal effects of each cell type on the other. Here, we detail the development of a 3D hydrogel co-culture system that allows us to culture different cell types for up to 7 days and subsequently separate and isolate the different cell populations using enzyme-sensitive glues. Separable 3D co-culture laminates were prepared by laminating PEG-based hydrogels with enzyme-degradable hydrogel adhesives. Encapsulated cell populations exhibited good segregation with well-defined interfaces. Furthermore, constructs can be separated on-demand upon addition of the appropriate enzyme, while cell viability remains high throughout the culture period, even after laminate separation. This platform offers great potential for a variety of basic cell signaling studies as the incorporation of an enzyme-sensitive adhesive interface allows the on-demand separation of individual cell populations for immediate analysis or further culture to examine persistence of co-culture effects and paracrine signaling on cell populations. See accompanying commentary by Danielle R. Bogdanowicz and Helen H. Lu DOI: 10.1002/biot.201300054     

Plantar shear stress distributions: comparing actual and predicted frictional forces at the foot-ground interface

Plantar shear stresses are believed to play a major role in diabetic ulceration. Due to the lack of commercial devices that can measure plantar shear distribution, a number of mathematical models have been developed to predict plantar frictional forces. This study assessed the accuracy of these models using a custom-built platform capable of measuring plantar stresses simultaneously. A total of 48 (38 healthy and 10 diabetic) human subjects (75+/-20 kg, 41+/-20 years, 32 males, 16 females) were recruited in the study. Plantar force data were collected for 2s at 50 Hz. Two models (M1 and M2) reported in the literature by different groups were used to predict local shear stresses. Root mean squared errors (RMSE) were calculated to compare model data with the actual data, focusing on three parameters: location, magnitude and timing of peak shear components. RMSE values of estimated peak AP and ML shear locations were 3.1 and 2.2 cm for M1 and 3.1 and 2.1cm for M2, respectively. Magnitude RMS error values for M1 were found to be 86.6 kPa in AP shear and 38.5 kPa in ML shear, whereas these values were determined to be 97.8 and 63.5 kPa, respectively by M2. Time to peak shear RMSE values averaged 17.2% in terms of the gait duration. In conclusion, distribution of plantar shear should be measured rather than predicted, particularly if one is interested in the magnitudes of shear components.