Dynamic permeability of the lacunar-canalicular system in human cortical bone.

A new method for the experimental determination of the permeability of a small sample of a fluid-saturated hierarchically structured porous material is described and applied to the determination of the lacunar-canalicular permeability [Formula: see text] in bone. The interest in the permeability of the lacunar-canalicular pore system (LCS) is due to the fact that the LCS is considered to be the site of bone mechanotransduction due to the loading-driven fluid flow over cellular structures. The permeability of this space has been estimated to be anywhere from [Formula: see text] to [Formula: see text]. However, the vascular pore system and LCS are intertwined, rendering the permeability of the much smaller-dimensioned LCS challenging to measure. In this study, we report a combined experimental and analytical approach that allowed the accurate determination of the [Formula: see text] to be on the order of [Formula: see text] for human osteonal bone. It was found that the [Formula: see text] has a linear dependence on loading frequency, decreasing at a rate of [Formula: see text]/Hz from 1 to 100 Hz, and using the proposed model, the porosity alone was able to explain 86 % of the [Formula: see text] variability.

Tortuosity and the Averaging of Microvelocity Fields in Poroelasticity.

The relationship between the macro- and microvelocity fields in a poroelastic representative volume element (RVE) has not being fully investigated. This relationship is considered to be a function of the tortuosity: a quantitative measure of the effect of the deviation of the pore fluid streamlines from straight (not tortuous) paths in fluid-saturated porous media. There are different expressions for tortuosity based on the deviation from straight pores, harmonic wave excitation, or from a kinetic energy loss analysis. The objective of the work presented is to determine the best expression for tortuosity of a multiply interconnected open pore architecture in an anisotropic porous media. The procedures for averaging the pore microvelocity over the RVE of poroelastic media by Coussy and by Biot were reviewed as part of this study, and the significant connection between these two procedures was established. Success was achieved in identifying the Coussy kinetic energy loss in the pore fluid approach as the most attractive expression for the tortuosity of porous media based on pore fluid viscosity, porosity, and the pore architecture. The fabric tensor, a 3D measure of the architecture of pore structure, was introduced in the expression of the tortuosity tensor for anisotropic porous media. Practical considerations for the measurement of the key parameters in the models of Coussy and Biot are discussed. In this study, we used cancellous bone as an example of interconnected pores and as a motivator for this study, but the results achieved are much more general and have a far broader application than just to cancellous bone.

Analytical basis for the determination of the lacunar-canalicular permeability of bone using cyclic loading.

An analytical model for the determination of the permeability in the lacunar-canalicular porosity of bone using cyclic loading is described in this contribution. The objective of the analysis presented is to relate the lacunar-canalicular permeability to a particular phase angle that is measurable when the bone is subjected to infinitesimal cyclic strain. The phase angle of interest is the lag angle between the applied strain and the resultant stress. Cyclic strain causes the interstitial fluid to move. This movement is essential for the viability of osteocytes and is believed to play a major role in the bone mechanotransduction mechanism. However, certain bone fluid flow properties, notably the permeability of the lacunar-canalicular porosity, are still not accurately determined. In this paper, formulas for the phase angle as a function of permeability for infinitesimal cyclic strain are presented and mathematical expressions for the storage modulus, loss modulus, and loss tangent are obtained. An accurate determination of the PLC permeability will improve our ability to understand mechanotransduction and mechanosensory mechanisms, which are fundamental to the understanding of how to treat osteoporosis, how to cope with microgravity in long-term manned space flights, and how to increase the longevity of prostheses that are implanted in bone tissue.

Perspectives on biological growth and remodeling.

The continuum mechanical treatment of biological growth and remodeling has attracted considerable attention over the past fifteen years. Many aspects of these problems are now well-understood, yet there remain areas in need of significant development from the standpoint of experiments, theory, and computation. In this perspective paper we review the state of the field and highlight open questions, challenges, and avenues for further development.

Topological optimization in hip prosthesis design.

With particular interest on total hip arthroplasty (THA), optimization of orthopedic prostheses is employed in this work to minimize the probability of implant failure or maximize prosthesis reliability. This goal is often identified with the reduction of stress concentrations at the interface between bone and these devices. However, aseptic loosening of the implant is mainly influenced by bone resorption phenomena revealed in some regions of the femur when a prosthesis is introduced. As a consequence, bone resorption appears due to stress shielding, that is to say the decrease of the stress level in the implanted femur caused by the significant load carrying of the prosthesis due to its higher stiffness. A maximum stiffness topological optimization-based (TO) strategy is utilized for non-linear static finite element (FE) analyses of the femur-implant assembly, with the goal of reducing stress shielding in the femur and to furnish guidelines for re-designing hip prostheses. This is accomplished by employing an extreme accuracy for both the three- dimensional reconstruction of the femur geometry and the material properties maps assigned as explicit functions of the local densities.

Mechanosensation and fluid transport in living bone.

The mechanosensory mechanisms in bone include (i) the cell system that is stimulated by external mechanical loading applied to the bone; (ii) the system that transduces that mechanical loading to a communicable signal; and (iii) the systems that transmit that signal to the effector cells for the maintenance of bone homeostasis and for strain adaptation of the bone structure. The effector cells are the osteoblasts and the osteoclasts. These systems and the mechanisms that they employ have not yet been unambiguously identified. The candidate systems will be reviewed. It will be argued that the current theoretical and experimental evidence suggests that osteocytes are the principal mechanosensory cells of bone, that they are activated by shear stress from fluid flowing through the osteocyte canaliculi, and that the electrically coupled three-dimensional network of osteocytes and lining cells is a communications system for the control of bone homeostasis and structural strain adaptation. The movement of bone fluid from the region of the bone vasculature through the canaliculi and the lacunae of the surrounding mineralized tissue accomplishes three important tasks. First, it transports nutrients to the osteocytes in the lacunae buried in the mineralized matrix. Second, it carries away the cell waste. Third, the bone fluid exerts a force on the cell process, a force that is large enough for the cell to sense. This is probably the basic mechanotransduction mechanism in bone, the way in which bone senses the mechanical load to which it is subjected. The mechanisms of bone fluid flow are described with particular emphasis on mechanotransduction. Also described is the cell to cell communication by which higher frequency signals might be transferred, a potential mechanism in bone by which the small whole tissue strain is amplified so the bone cells can respond to it. One of the conclusions is that higher frequency low amplitude strains can maintain bone as effectively as low frequency high amplitude strains. This conclusion leads to a paradigm shift in how to treat osteoporosis and how to cope with microgravity.

Estimation of the effective transversely isotropic elastic constants of a material from known values of the material's orthotropic elastic constants.

A method is illustrated for determining the effective transversely isotropic (or isotropic) elastic constants from measured orthotropic elastic constants. This method consists of constructing upper and lower bounds on the effective transversely isotropic (or isotropic) elastic constants using the known orthotropic values. This method is illustrated using three sets of elastic constants for bone. Fortunately, the upper and lower bounds are very close. Thus very good approximations for the effective transversely isotropic (or isotropic) elastic constants for cortical and cancellous bone are obtained from previously published data on the orthotropic elastic constants for those tissue types. This work is undertaken to build a greater database for the transversely isotropic elastic constants of bone with the intention of employing them in a transversely isotropic model of bone poroelasticity. An interesting aspect of the present result is that the Voigt and Reuss bounds are very tight for these anisotropic materials. This is not always the case for these bounds.

A model for strain amplification in the actin cytoskeleton of osteocytes due to fluid drag on pericellular matrix.

A model is presented that provides a resolution to a fundamental paradox in bone physiology, namely, that the strains applied to whole bone (i.e., tissue level strains) are much smaller (0.04-0.3 percent) than the strains (1-10 percent) that are necessary to cause bone signaling in deformed cell cultures (Rubin and Lanyon, J. Bone Joint Surg. 66A (1984) 397-410; Fritton et al., J. Biomech. 33 (2000) 317-325). The effect of fluid drag forces on the pericellular matrix (PM), its coupling to the intracellular actin cytoskeleton (IAC) and the strain amplification that results from this coupling are examined for the first time. The model leads to two predictions, which could fundamentally change existing views. First, for the loading range 1-20MPa and frequency range 1-20Hz, it is, indeed, possible to produce cellular level strains in bone that are up to 100 fold greater than normal tissue level strains (0.04-0.3 percent). Thus, the strain in the cell process membrane due to the loading can be of the same order as the in vitro strains measured in cell culture studies where intracellular biochemical responses are observed for cells on stretched elastic substrates. Second, it demonstrates that in any cellular system, where cells are subject to fluid flow and tethered to more rigid supporting structures, the tensile forces on the cell due to the drag forces on the tethering fibers may be many times greater than the fluid shear force on the cell membrane.

Modeling tracer transport in an osteon under cyclic loading.

A mathematical model is developed to explain the fundamental conundrum as to how during cyclic mechanical loading there can be net solute (e.g., nutrient, tracer) transport in bone via the lacunar-canalicular porosity when there is no net fluid movement in the canaliculi over a loading cycle. Our hypothesis is that the fluid space in an osteocytic lacuna facilitates a nearly instantaneous mixing process of bone fluid that creates a difference in tracer concentration between the inward and outward canalicular flow and thus ensures net tracer transport to the osteocytes during cyclic loading, as has been shown experimentally. The sequential spread of the tracer from the osteonal canal to the lacunae is investigated for an osteon experiencing sinusoidal loading. The fluid pressure in the canaliculi is calculated using poroelasticity theory and the mixing process in the lacunae is then simulated computationally. The tracer concentration in lacunae extending radially from the osteonal canal to the cement line is calculated as a function of the loading frequency, loading magnitude, and number of loading cycles as well as the permeability of the lacunar-canalicular porosity. Our results show that net tracer transport to the lacunae does occur for cyclic loading. Tracer transport is found to increase with higher loading magnitude and higher permeability and to decrease with increasing loading frequency. This work will be helpful in designing experimental studies of tracer movement and bone fluid flow, which will enhance our understanding of bone metabolism as well as bone adaptation.

How is a tissue built?

Tissues change in many ways in the period that they are part of a living organism. They are created in fairly repeatable structural patterns, and we know that the patterns are due to both the genes and the (mechanical) environment, but we do not know exactly what part or percentage of a particular pattern to consider the genes, or the environment, responsible for. We do not know much about the beginning of tissue construction (morphogenesis) and we do not know the methods of tissue construction. When the tissue structure is altered to accommodate a new loading, we do not know how the decision is made for the structural reconstruction. We do know that tissues grow or reconstruct themselves without ceasing to continue with their structural function, but we do not understand the processes that permit them to accomplish this. Tissues change their structures to altered mechanical environments, but we are not sure how. Tissues heal themselves and we understand little of the structural mechanics of the process. With the objective of describing the interesting unsolved mechanics problems associated with these biological processes, some aspects of the formation, growth, and adaptation of living tissues are reviewed. The emphasis is on ideas and models. Beyond the objective is the hope that the work will stimulate new ideas and new observations in the mechanical and chemical aspects of developmental biology.

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.

Bone poroelasticity.

Poroelasticity is a well-developed theory for the interaction of fluid and solid phases of a fluid-saturated porous medium. It is widely used in geomechanics and has been applied to bone by many authors in the last 30 years. The purpose of this work is, first, to review the literature related to the application of poroelasticity to the interstitial bone fluid and, second, to describe the specific physical and modeling considerations that establish poroelasticity as an effective and useful model for deformation-driven bone fluid movement in bone tissue. The application of poroelasticity to bone differs from its application to soft tissues in two important ways. First, the deformations of bone are small while those of soft tissues are generally large. Second, the bulk modulus of the mineralized bone matrix is about six times stiffer than that of the fluid in the pores while the bulk moduli of the soft tissue matrix and the pore water are almost the same. Poroelasticity and electrokinetics can be used to explain strain-generated potentials in wet bone. It is noted that strain-generated potentials can be used as an effective tool in the experimental study of local bone fluid flow, and that the knowledge of this technique will contribute to the answers of a number of questions concerning bone mineralization, osteocyte nutrition and the bone mechanosensory system.

The effect of surface roughness on the stress adaptation of trabecular architecture around a cylindrical implant.

The effect of implant-bone bonding and the effect of implant surface roughness on bone remodeling near the bone-implant interface were studied by using a surface remodeling theory and the boundary element method. The study has shown that implant attachment plays an important role in bone remodeling near the implant. It has been observed in animal experiments and in clinical situations that the remodeled trabecular bone architecture around a cylindrical implant could vary, on one hand, from a hub surrounding the implant with a set of external spokes to, on the other hand, a hubless situation in which a set of spokes attach directly to the implant. It is shown here that the difference in these structures may be attributed to differences in implant attachment. The results show that the bone with perfect bonding or roller boundary condition without a gap remodeled to a hubless spoke trabecular bone architecture. On the other hand, the roller boundary condition with a specified gap yielded a spoke trabecular architecture with a hub or ring surrounding the implant. These quantitative results mirror the experimental and clinical observations. It is concluded that the hub is a consequence of the gap and not a consequence of the lack of friction between the implant and the bone.

Structural changes in living tissues.

Tissues change in many ways in the period that they are part of a living organism. Tissues are created in fairly repeatable structural patterns, and the patterns are due to both the genes and the (mechanical) environment, but we do not know exactly what part or percentage of a particular pattern to consider the genes, or the environment, responsible for. We do not know much about the beginning of tissue construction (morphogenesis) and also the methods of tissue construction. When the tissue structure is altered to accommodate a new loading, it is not known how the decision is made for the structural reconstruction. We know that tissues grow or reconstruct themselves without ceasing to continue with their structural function, but we do not understand the processes that permit them to accomplish this. Tissues change their structures to altered mechanical environments, but we are not sure how. Tissues heal themselves and we understand little of the structural mechanics of the process. With the objective of describing the interesting unsolved mechanics problems associated with these biological processes, some aspects of the formation, growth and adaptation of living tissues are reviewed. Beyond the objective is the hope that the work will stimulate new ideas and new observations in developmental biology.

Imposing thermodynamic restrictions on the elastic constant-fabric tensor relationship.

Zysset and Curnier (1995) rederive the general representation for the fourth rank elasticity tensor in terms of the fabric tensor given originally by Cowin (1985). Zysset and Curnier (1995) then obtain a form of the representation that is sufficiently restricted to satisfy the thermodynamic constraint that the compliance tensor be positive definite. Zysset and Curnier (1995, 1996) argue that their representation has an advantage over that of Cowin (1985) because it has, a priori, been sufficiently restricted to satisfy the thermodynamic constraint that the compliance tensor be positive definite. The purpose of this note is to point out that the limitations of the Zysset and Curnier (1995, 1996) sufficiency representation are quite severe, restricting the flexibility of the model by allowing only five of the nine normally independent orthotropic elastic coefficients to be independent. Potential users should consider employing the Cowin (1985) representation and subsequently imposing the condition that the elasticity tensor be positive definite.

Strain amplification in the bone mechanosensory system.

This article discusses the potential mechanisms by which the strain induced at the membrane of an osteocyte may be amplified from the strain experienced by the whole bone due to mechanical loading. These mechanisms address the question of how these mechanical load-induced small strains of (typically) about 0.1% (but up to 0.5%) applied to a whole bone are amplified to strains of 1% or larger at the membrane of the osteocyte buried in its lacuna in the bone matrix. The answer to this question is an important link in the mechanosensory system in bone and in relating in vitro cell studies to in vivo cellular response.

Electrical signal transmission in a bone cell network: the influence of a discrete gap junction.

A refined electrical cable model is formulated to investigate the role of a discrete gap junction in the intracellular transmission of electrical signals in an electrically coupled system of osteocytes and osteoblasts in an osteon. The model also examines the influence of the ratio q between the membrane's electrical time constant and the characteristic time of pore fluid pressure, the circular, cylindrical geometry of the osteon, and key simplifying assumptions in our earlier continuous cable model (see Zhang, D., S. C. Cowin, and S. Weinbaum. Electrical signal transmission and gap junction regulation in a bone cell network: A cable model for an osteon. Ann. Biomed. Eng. 25:379-396, 1997). Using this refined model, it is shown that (1) the intracellular potential amplitude at the osteoblastic end of the osteonal cable retains the character of a combination of a low-pass and a high-pass filter as the corner frequency varies in the physiological range; (2) the presence of a discrete gap junction near a resting osteoblast can lead to significant modulation of the intracellular potential and current in the osteoblast for measured values of the gap junction coupling strength; and (3) the circular, cylindrical geometry of the osteon is well simulated by the beam analogy used in Zhang et al.

On mechanosensation in bone under microgravity.

Recent research has illuminated the biological response of bone to mechanical loading at the cellular level, but the precise mechanosensory system that signals bone cells to deposit or resorb tissue has not been identified. The purpose of this paper is to describe the current status of this research and to suggest some possible mechanosensory systems by which bone cells might sense environmental loads. The question of whether the mechanosensory system in bone tissue is at the level of the cell or whether it is at the tissue level and involving the cells is considered here. More precisely, the following question is addressed: can an osteocyte or an osteoblast read the gravitational field changes directly (and independent of changes in its environment), or does it detect those changes indirectly from its environment by contact stresses as it must detect other changes in mechanical loading on the surface of the earth? Our strategies for coping with the decay of the musculoskeletal system in long term space flight are somewhat dependent upon the answer to this question.

The anisotropic Hooke's law for cancellous bone and wood.

A method of data analysis for a set of elastic constant measurements is applied to data bases for wood and cancellous bone. For these materials the identification of the type of elastic symmetry is complicated by the variable composition of the material. The data analysis method permits the identification of the type of elastic symmetry to be accomplished independent of the examination of the variable composition. This method of analysis may be applied to any set of elastic constant measurements, but is illustrated here by application to hardwoods and softwoods, and to an extraordinary data base of cancellous bone elastic constants. The solid volume fraction or bulk density is the compositional variable for the elastic constants of these natural materials. The final results are the solid volume fraction dependent orthotropic Hooke's law for cancellous bone and a bulk density dependent one for hardwoods and softwoods.