Fluid pressure relaxation depends upon osteonal microstructure: modeling an oscillatory bending experiment.

When bone is mechanically loaded, bone fluid flow induces shear stresses on bone cells that have been proposed to be involved in bone's mechanosensory system. To investigate bone fluid flow and strain-generated potentials, several theoretical models have been proposed to mimic oscillatory four-point bending experiments performed on thin bone specimens. While these previous models assume that the bone fluid relaxes across the specimen thickness, we hypothesize that the bone fluid relaxes primarily through the vascular porosity (osteonal canals) instead and develop a new poroelastic model that integrates the microstructural details of the lacunar-canalicular porosity, osteonal canals, and the osteonal cement lines. Local fluid pressure profiles are obtained from the model, and we find two different fluid relaxation behaviors in the bone specimen, depending on its microstructure: one associated with the connected osteonal canal system, through which bone fluid relaxes to the nearby osteonal canals; and one associated with the thickness of a homogeneous porous bone specimen (approximately 1 mm in our model), through which bone fluid relaxes between the external surfaces of the bone specimen at relatively lower loading frequencies. Our results suggest that in osteonal bone specimens the fluid pressure response to cyclic loading is not sensitive to the permeability of the osteonal cement lines, while it is sensitive to the applied loading frequency. Our results also reveal that the fluid pressure gradients near the osteonal canals (and thus the fluid shear stresses acting on the nearby osteocytes) are significantly amplified at higher loading frequencies.     

Osteoblasts respond to pulsatile fluid flow with short-term increases in PGE(2) but no change in mineralization.

Although there is no consensus as to the precise nature of the mechanostimulatory signals imparted to the bone cells during remodeling, it has been postulated that deformation-induced fluid flow plays a role in the mechanotransduction pathway. In vitro, osteoblasts respond to fluid shear stress with an increase in PGE(2) production; however, the long-term effects of fluid shear stress on cell proliferation and differentiation have not been examined. The goal of this study was to apply continuous pulsatile fluid shear stresses to osteoblasts and determine whether the initial production of PGE(2) is associated with long-term biochemical changes. The acute response of bone cells to a pulsatile fluid shear stress (0.6 +/- 0.5 Pa, 3.0 Hz) was characterized by a transient fourfold increase in PGE(2) production. After 7 days of static culture (0 dyn/cm(2)) or low (0.06 +/- 0.05 Pa, 0.3 Hz) or high (0.6 +/- 0.5 Pa, 3.0 Hz) levels of pulsatile fluid shear stress, the bone cells responded with an 83% average increase in cell number, but no statistical difference (P > 0.53) between the groups was observed. Alkaline phosphatase activity per cell decreased in the static cultures but not in the low- or high-flow groups. Mineralization was also unaffected by the different levels of applied shear stress. Our results indicate that short-term changes in PGE(2) levels caused by pulsatile fluid flow are not associated with long-term changes in proliferation or mineralization of bone cells.     

Temporal gradients in shear stimulate osteoblastic proliferation via ERK1/2 and retinoblastoma protein

Bone cells are subject to interstitial fluid flow (IFF) driven by venous pressure and mechanical loading. Rapid dynamic changes in mechanical loading cause transient gradients in IFF. The effects of pulsatile flow (temporal gradients in fluid shear) on rat UMR106 cells and rat primary osteoblastic cells were studied. Pulsatile flow induced a 95% increase in S-phase UMR106 cells compared with static controls. In contrast, ramped steady flow stimulated only a 3% increase. Similar patterns of S-phase induction were also observed in rat primary osteoblastic cells. Pulsatile flow significantly increased relative UMR106 cell number by 37 and 62% at 1.5 and 24 h, respectively. Pulsatile flow also significantly increased extracellular signal-regulated kinase (ERK1/2) phosphorylation by 418%, whereas ramped steady flow reduced ERK1/2 activation to 17% of control. Correspondingly, retinoblastoma protein was significantly phosphorylated by pulsatile fluid flow. Inhibition of mitogen-activated protein (MAP)/ERK kinase (MEK)1/2 by U0126 (a specific MEK1/2 inhibitor) reduced shear-induced ERK1/2 phosphorylation and cell proliferation. These findings suggest that temporal gradients in fluid shear stress are potent stimuli of bone cell proliferation.     

Influence of interstitial bone microcracks on strain-induced fluid flow

It is well known that microcracks act as a stimulus for bone remodelling, initiating resorption by osteoclasts and new bone formation by osteoblasts. Moreover, microcracks are likely to alter the fluid flow and convective transport through the bone tissue. This paper proposes a quantitative evaluation of the strain-induced interstitial fluid velocities developing in osteons in presence of a microcrack in the interstitial bone tissue. Based on Biot theory in the low-frequency range, a poroelastic model is carried out to study the hydro-mechanical behaviour of cracked osteonal tissue. The finite element results show that the presence of a microcrack in the interstitial osteonal tissue may drastically reduce the fluid velocity inside the neighbouring osteons. This fluid inactive zone inside osteons can cover up to 10% of their surface. Consequently, the fluid environment of bone mechano-sensitive cells is locally modified.     

Fluid pressure gradients, arising from oscillations in intramedullary pressure, is correlated with the formation of bone and inhibition of intracortical porosity

Fluid flow that arises from the functional loading of bone tissue has been proposed to be a critical regulator of skeletal mass and morphology. To test this hypothesis, the bone adaptive response to a physiological fluid stimulus, driven by low magnitude, high frequency oscillations of intramedullary pressure (ImP), were examined, in which fluid pressures were achieved without deforming the bone tissue. The ulnae of adult turkeys were functionally isolated via transverse epiphyseal osteotomies, and the adaptive response to four weeks of disuse (n=5) was compared to disuse plus 10 min per day of a physiological sinusoidal fluid pressure signal (60 mmHg, 20Hz). Disuse alone resulted in significant bone loss (5.7+/-1.9%, p< or =0.05), achieved by thinning the cortex via endosteal resorption and an increase in intracortical porosity. By also subjecting bone to oscillatory fluid flow, a significant increase in bone mass at the mid-diaphysis (18.3+/-7.6%, p<0.05), was achieved by both periosteal and endosteal new bone formation. The spatial distribution of the transcortical fluid pressure gradients (inverted Delta P(r)), a parameter closely related to fluid velocity and fluid shear stress, was quantified in 12 equal sectors across a section at the mid-diaphyses. A strong correlation was found between the inverted Delta P(r) and total new bone formation (r=0.75, p=0.01); and an inverse correlation (r=-0.75, p=0.01) observed between inverted Delta P(r) and the area of increased intracortical porosity, indicating that fluid flow signals were necessary to maintain bone mass and/or inhibit bone loss against the challenge of disuse. By generating this fluid flow in the absence of matrix strain, these data suggest that anabolic fluid movement plays a regulatory role in the modeling and remodeling process. While ImP increases uniformly in the marrow cavity, the distinct parameters of fluid flow vary substantially due to the geometry and ultrastructure of bone, which ultimately defines the spatial non-uniformity of the adaptive process.     

Bone tissue engineering: the role of interstitial fluid flow

It is well established that vascularization is required for effective bone healing. This implies that blood flow and interstitial fluid (ISF) flow are required for healing and maintenance of bone. The fact that changes in bone blood flow and ISF flow are associated with changes in bone remodeling and formation support this theory. ISF flow in bone results from transcortical pressure gradients produced by vascular and hydrostatic pressure, and mechanical loading. Conditions observed to alter flow rates include increases in venous pressure in hypertension, fluid shifts occurring in bedrest and microgravity, increases in vascularization during the injury-healing response, and mechanical compression and bending of bone during exercise. These conditions also induce changes in bone remodeling. Previously, we hypothesized that interstitial fluid flow in bone, and in particular fluid shear stress, serves to mediate signal transduction in mechanical loading- and injury-induced remodeling. In addition, we proposed that a lack or decrease of ISF flow results in the bone loss observed in disuse and microgravity. The purpose of this article is to review ISF flow in bone and its role in osteogenesis.     

Influence of cortical canal architecture on lacunocanalicular pore pressure and fluid flow

Bone is a dynamic tissue that undergoes structural modification in response to its mechanical environment, but how bone cells sense and respond to loading conditions remains incompletely understood. Current theories focus on strain-induced fluid flow for the primary means of mechanotransduction. To examine the influence of age-related cortical rarefaction on lacunocanalicular fluid characteristics, coupled fluid flow and mechanical computational models of bone specimens representing young, mid-age and aged samples were derived artificially from the same original micro-computed tomography image data. Simulated mechanical loading was applied to the bone models to induce pressure-driven interstitial fluid flow. Results demonstrated a decrease in pore pressure and fluid velocity magnitudes with age as a result of increased cortical porosity. Mean canal separation, as opposed to canal size, was implicated as a primary factor affecting age-related fluid dynamics. Future investigations through refinement of the model may implicate fluid stasis or inadequate nutrient transport experienced by osteocytes as a key factor in the initiation of cortical remodelling events.     

Influence of cortical canal architecture on lacunocanalicular pore pressure and fluid flow

Bone is a dynamic tissue that undergoes structural modification in response to its mechanical environment, but how bone cells sense and respond to loading conditions remains incompletely understood. Current theories focus on strain-induced fluid flow for the primary means of mechanotransduction. To examine the influence of age-related cortical rarefaction on lacunocanalicular fluid characteristics, coupled fluid flow and mechanical computational models of bone specimens representing young, mid-age and aged samples were derived artificially from the same original micro-computed tomography image data. Simulated mechanical loading was applied to the bone models to induce pressure-driven interstitial fluid flow. Results demonstrated a decrease in pore pressure and fluid velocity magnitudes with age as a result of increased cortical porosity. Mean canal separation, as opposed to canal size, was implicated as a primary factor affecting age-related fluid dynamics. Future investigations through refinement of the model may implicate fluid stasis or inadequate nutrient transport experienced by osteocytes as a key factor in the initiation of cortical remodelling events.     

Effect of oscillating fluid flow stimulation on osteocyte mRNA expression

Structural adaptation of the bone tissue is mediated by loading-induced interstitial fluid flow within the bone microstructure. Within this framework, osteocytes fulfill the central mechanotransductive role in the bone remodeling process. While osteocytes have been demonstrated to be exquisitely sensitive to various forms of fluid flow stimulus in vitro, the effect of different oscillating fluid flow (OFF) parameters on osteocyte activity has yet to be systematically characterized. In this study, we investigate the effect of three OFF parameters on osteocyte activity in vitro and hypothesize that COX-2, RANKL, and OPG mRNA expression in osteocytes are sensitive to the OFF parameters: peak shear stress amplitude (0.5 Pa, 1 Pa, 2 Pa, and 5 Pa), oscillating frequency (0.5 Hz, 1 Hz, and 2 Hz), and total flow duration (1 h, 2 h, and 4 h). Our findings demonstrate that COX-2 mRNA levels are elevated in osteocytes subjected to higher peak shear stress amplitudes and longer flow durations, while RANKL/OPG mRNA levels decreased to a minimum threshold in response to higher peak shear stress amplitudes, faster oscillating frequencies, and longer flow durations. These findings suggest that dynamic fluid flow with higher peak shear stress amplitudes, faster oscillating frequencies, and longer loading durations provide the best conditions for promoting bone formation.     

Experimental measurement of dynamic fluid shear stress on the aortic surface of the aortic valve leaflet

Aortic valve (AV) calcification is a highly prevalent disease with serious impact on mortality and morbidity. Although exact causes and mechanisms of AV calcification are unclear, previous studies suggest that mechanical forces play a role. Since calcium deposits occur almost exclusively on the aortic surfaces of AV leaflets, it has been hypothesized that adverse patterns of fluid shear stress on the aortic surface of AV leaflets promote calcification. The current study characterizes AV leaflet aortic surface fluid shear stresses using Laser Doppler velocimetry and an in vitro pulsatile flow loop. The valve model used was a native porcine valve mounted on a suturing ring and preserved using 0.15% glutaraldehyde solution. This valve model was inserted in a mounting chamber with sinus geometries, which is made of clear acrylic to provide optical access for measurements. To understand the effects of hemodynamics on fluid shear stress, shear stress was measured across a range of conditions: varying stroke volumes at the same heart rate and varying heart rates at the same stroke volume. Systolic shear stress magnitude was found to be much higher than diastolic shear stress magnitude due to the stronger flow in the sinuses during systole, reaching up to 20 dyn/cm² at mid-systole. Upon increasing stroke volume, fluid shear stresses increased due to stronger sinus fluid motion. Upon increasing heart rate, fluid shear stresses decreased due to reduced systolic duration that restricted the formation of strong sinus flow. Significant changes in the shear stress waveform were observed at 90 beats/min, most likely due to altered leaflet dynamics at this higher heart rate. Overall, this study represents the most well-resolved shear stress measurements to date across a range of conditions on the aortic side of the AV. The data presented can be used for further investigation to understand AV biological response to shear stresses.     

A novel in vitro loading system to produce supraphysiologic oscillatory fluid shear stress

A multi-well fluid loading (MFL) system was developed to deliver oscillatory subphysiologic to supraphysiologic fluid shear stresses to cell monolayers in vitro using standard multi-well culture plates. Computational fluid dynamics modeling with fluid-structure interactions was used to quantify the squeeze film fluid flow between an axially displaced piston and the well plate surface. Adjusting the cone angle of the piston base modulated the fluid pressure, velocity, and shear stress magnitudes. Modeling results showed that there was near uniform fluid shear stress across the well with a linear drop in pressure across the radius of the well. Using the MFL system, RAW 264.7 osteoclastic cells were exposed to oscillatory fluid shear stresses of 0, 0.5, 1.5, 4, 6, and 17 Pa. Cells were loaded 1 h per day at 1 Hz for two days. Compared to sub-physiologic and physiologic levels, supraphysiologic oscillatory fluid shear induced upregulation of osteoclastic activity as measured by tartrate-resistant acid phosphatase activity and formation of mineral resorption pits. Cell number remained constant across all treatment groups.     

Nitric oxide production by bone cells is fluid shear stress rate dependent

Shear stress due to mechanical loading-induced flow of interstitial fluid through the lacuno-canalicular network is a likely signal for bone cell adaptive responses. Moreover, the rate (determined by frequency and magnitude) of mechanical loading determines the amount of bone formation. Whether the bone cells' response to fluid shear stress is rate dependent is unknown. Here we investigated whether bone cell activation by fluid shear stress is rate dependent. MC3T3-E1 osteoblastic cells were subjected for 15 min to fluid shear stress of varying frequencies and amplitudes, resulting in peak fluid shear stress rates ranging from 0 to 39.6 Pa-Hz. Nitric oxide production, a parameter for bone cell activation, was found to be linearly dependent on the fluid shear stress rate; the slope was steepest at 5 min (0.11 Pa-Hz(-1)) and decreased to 0.03 Pa-Hz(-1) at 15 min. We conclude that the fluid shear stress rate is an important parameter for bone cell activation.     

The effect of asymmetry in abdominal aortic aneurysms under physiologically realistic pulsatile flow conditions

In the abdominal segment of the human aorta under a patient's average resting conditions, pulsatile blood flow exhibits complex laminar patterns with secondary flows induced by adjacent branches and irregular vessel geometries. The flow dynamics becomes more complex when there is a pathological condition that causes changes in the normal structural composition of the vessel wall, for example, in the presence of an aneurysm. This work examines the hemodynamics of pulsatile blood flow in hypothetical three-dimensional models of abdominal aortic aneurysms (AAAs). Numerical predictions of blood flow patterns and hemodynamic stresses in AAAs are performed in single-aneurysm, asymmetric, rigid wall models using the finite element method. We characterize pulsatile flow dynamics in AAAs for average resting conditions by means of identifying regions of disturbed flow and quantifying the disturbance by evaluating flow-induced stresses at the aneurysm wall, specifically wall pressure and wall shear stress. Physiologically realistic abdominal aortic blood flow is simulated under pulsatile conditions for the range of time-average Reynolds numbers 50 < or = Rem < or = 300, corresponding to a range of peak Reynolds numbers 262.5 < or = Repeak < or = 1575. The vortex dynamics induced by pulsatile flow in AAAs is depicted by a sequence of four different flow phases in one period of the cardiac pulse. Peak wall shear stress and peak wall pressure are reported as a function of the time-average Reynolds number and aneurysm asymmetry. The effect of asymmetry in hypothetically shaped AAAs is to increase the maximum wall shear stress at peak flow and to induce the appearance of secondary flows in late diastole.     

The effects of external compression on venous blood flow and tissue deformation in the lower leg

External pneumatic compression of the lower legs is effective as prophylaxis against deep vein thrombosis. In a typical application, inflatable cuffs are wrapped around the patient's legs and periodically inflated to prevent stasis, accelerate venous blood flow, and enhance fibrinolysis. The purpose of this study was to examine the stress distribution within the tissues, and the corresponding venous blood flow and intravascular shear stress with different external compression modalities. A two-dimensional finite element analysis (FEA) was used to determine venous collapse as a function of internal (venous) pressure and the magnitude and spatial distribution of external (surface) pressure. Using the one-dimensional equations governing flow in a collapsible tube and the relations for venous collapse from the FEA, blood flow resulting from external compression was simulated. Tests were conducted to compare circumferentially symmetric (C) and asymmetric (A) compression and to examine distributions of pressure along the limb. Results show that A compression produces greater vessel collapse and generates larger blood flow velocities and shear stresses than C compression. The differences between axially uniform and graded-sequential compression are less marked than previously found, with uniform compression providing slightly greater peak flow velocities and shear stresses. The major advantage of graded-sequential compression is found at midcalf. Strains at the lumenal border are approximately 20 percent at an external pressure of 50 mmHg (6650 Pa) with all compression modalities.     

The mechanobiological effects of periosteal surface loads

We have developed an improved mechanobiological model of bone morphogenesis and functional adaptation that includes the influences of periosteum tension and pressure on bone formation and resorption. Previous models assumed that periosteal and endosteal bone deposition and resorption rates are governed only by the local intracortical daily stress or strain stimulus caused by cyclic loading. The new model incorporates experimental findings that pressures on periosteal surfaces can impede bone formation or induce bone resorption, whereas periosteal tensile strains perpendicular to bone surfaces can impede bone resorption or induce bone formation. We propose that these effects can produce flattened or concave bone surfaces in regions of periosteal pressure and bone ridges in regions of periosteal tension. The model was implemented with computer simulations to illustrate the role of adjacent muscles on the development of the triangular cross-sectional geometry of the rat tibia. The results suggest that intracortical stresses dictate bone size, whereas periosteal pressures may work in combination with intracortical stresses and other mechanobiological factors in the development of local bone cross-sectional shapes.     

Cells in fluidic environments are sensitive to flow frequency

Virtually all cells accommodate to their mechanical environment. In particular, cells subject to flow respond to rapid changes in fluid shear stress (SS), cyclic stretch (CS), and pressure. Recent studies have focused on the effect of pulsatility on cellular behavior. Since cells of many different tissue beds are constantly exposed to fluid flows over a narrow range of frequencies, we hypothesized that an intrinsic flow frequency that is optimal for determining cell phenotype exists. We report here that cells from various tissue beds (bovine aortic endothelial cells (BAEC), rat small intestine epithelial cells (RSIEC), and rat lung epithelial cells (RLEC)) proliferate maximally when cultured in a perfusion bioreactor under pulsatile conditions at a specific frequency, independent of the applied SS. Vascular endothelial and pulmonary epithelial cell proliferation peaked under 1 Hz pulsatile flow. In contrast, proliferation of gastrointestinal cells, which in their physiological context are subject to no flow or higher wavelength signal, was maximum at 0.125 Hz or under no flow. Moreover, exposure of BAEC to pulsatile flow of varying frequency influenced their nitric oxide synthase activity and prostacyclin production, which reached maximum values at 1 Hz. Notably, the "optimal" frequencies for the cell types examined correspond to the physiologic operating range of the organs from where they were initially derived. These findings suggest that frequency, independent of shear, is an essential determinant of cell response in pulsatile environments.     

A fluid--structure interaction finite element analysis of pulsatile blood flow through a compliant stenotic artery.

A new model is used to analyze the fully coupled problem of pulsatile blood flow through a compliant, axisymmetric stenotic artery using the finite element method. The model uses large displacement and large strain theory for the solid, and the full Navier-Stokes equations for the fluid. The effect of increasing area reduction on fluid dynamic and structural stresses is presented. Results show that pressure drop, peak wall shear stress, and maximum principal stress in the lesion all increase dramatically as the area reduction in the stenosis is increased from 51 to 89 percent. Further reductions in stenosis cross-sectional area, however, produce relatively little additional change in these parameters due to a concomitant reduction in flow rate caused by the losses in the constriction. Inner wall hoop stretch amplitude just distal to the stenosis also increases with increasing stenosis severity, as downstream pressures are reduced to a physiological minimum. The contraction of the artery distal to the stenosis generates a significant compressive stress on the downstream shoulder of the lesion. Dynamic narrowing of the stenosis is also seen, further augmenting area constriction at times of peak flow. Pressure drop results are found to compare well to an experimentally based theoretical curve, despite the assumption of laminar flow.     

Traction Forces of Endothelial Cells under Slow Shear Flow

Endothelial cells are constantly exposed to fluid shear stresses that regulate vascular morphogenesis, homeostasis, and disease. The mechanical responses of endothelial cells to relatively high shear flow such as that characteristic of arterial circulation has been extensively studied. Much less is known about the responses of endothelial cells to slow shear flow such as that characteristic of venous circulation, early angiogenesis, atherosclerosis, intracranial aneurysm, or interstitial flow. Here we used a novel, to our knowledge, microfluidic technique to measure traction forces exerted by confluent vascular endothelial cell monolayers under slow shear flow. We found that cells respond to flow with rapid and pronounced increases in traction forces and cell-cell stresses. These responses are reversible in time and do not involve reorientation of the cell body. Traction maps reveal that local cell responses to slow shear flow are highly heterogeneous in magnitude and sign. Our findings unveil a low-flow regime in which endothelial cell mechanics is acutely responsive to shear stress.     

Micromechanically based poroelastic modeling of fluid flow in Haversian bone

To explore the hypothesis that load-induced fluid flow in bone is a mechano-transduction mechanism in bone adaptation, unit cell micro-mechanical techniques are used to relate the microstructure of Haversian cortical bone to its effective poroelastic properties. Computational poroelastic models are then applied to compute in vitro Haversian fluid flows in a prismatic specimen of cortical bone during harmonic bending excitations over the frequency range of 10(0) to 10(6) Hz. At each frequency considered, the steady state harmonic response of the poroelastic bone specimen is computed using complex frequency-domain finite element analysis. At the higher frequencies considered, the breakdown of Poisueille flow in Haversian canals is modeled by introduction of a complex fluid viscosity. Peak bone fluid pressures are found to increase linearly with loading frequency in proportion to peak bone stress up to frequencies of approximately 10 kHz. Haversian fluid shear stresses are found to increase linearly with excitation frequency and loading magnitude up until the breakdown of Poisueille flow. Tan delta values associated with the energy dissipated by load-induced fluid flow are also compared with values measured experimentally in a concurrent broadband spectral analysis of bone. The computational models indicate that fluid shear stresses and fluid pressures in the Haversian system could, under physiologically realistic loading, easily reach the level of a few Pascals, which have been shown in other works to elicit cell responses in vitro.     

The congenital bicuspid aortic valve can experience high-frequency unsteady shear stresses on its leaflet surface

The bicuspid aortic valve (BAV) is a common congenital malformation of the aortic valve (AV) affecting 1% to 2% of the population. The BAV is predisposed to early degenerative calcification of valve leaflets, and BAV patients constitute 50% of AV stenosis patients. Although evidence shows that genetic defects can play a role in calcification of the BAV leaflets, we hypothesize that drastic changes in the mechanical environment of the BAV elicit pathological responses from the valve and might be concurrently responsible for early calcification. An in vitro model of the BAV was constructed by surgically manipulating a native trileaflet porcine AV. The BAV valve model and a trileaflet AV (TAV) model were tested in an in vitro pulsatile flow loop mimicking physiological hemodynamics. Laser Doppler velocimetry was used to make measurements of fluid shear stresses on the leaflet of the valve models using previously established methodologies. Furthermore, particle image velocimetry was used to visualize the flow fields downstream of the valves and in the sinuses. In the BAV model, flow near the leaflets and fluid shear stresses on the leaflets were much more unsteady than for the TAV model, most likely due to the moderate stenosis in the BAV and the skewed forward flow jet that collided with the aorta wall. This additional unsteadiness occurred during mid- to late-systole and was composed of cycle-to-cycle magnitude variability as well as high-frequency fluctuations about the mean shear stress. It has been demonstrated that the BAV geometry can lead to unsteady shear stresses under physiological flow and pressure conditions. Such altered shear stresses could play a role in accelerated calcification in BAVs.