Strain Concentrations Surrounding an Ellipsoid Model of Lacunae and Osteocytes

Direct cell sensing of tissue matrix strains is one possible signaling mechanism for mechanically mediated bone adaptation. We utilized homogenization theory lo estimate bone tissue matrix strains surrounding osteocytes using two sets of models. The first set of models estimated the strain levels surrounding the lacunae and canaliculi, taking into account variations in lamellar properties. The second set estimated strain levels in the osteocyte and the surrounding matrix for different cellular mechanical properties. The results showed that the strain levels found in and surrounding osteocytes, 1700 to 2700 microstrain (denoted as μe; 1 =.0001% strain), were significantly greater than the trabecular tissue level strains of [1325 μe, 287 μe, 87 μe] used for model input. Variation in lamellar properties did not affect strain levels, except at lamellar boundaries. Strain in and surrounding the osteocyte was not significantly affected by cellular stiffness ranging between 28 and 28,000 Pascals (Pa). Strain levels surrounding lacunae and canaliculi were approximately equivalent.     

Strain Concentrations Surrounding an Ellipsoid Model of Lacunae and Osteocytes

Direct cell sensing of tissue matrix strains is one possible signaling mechanism for mechanically mediated bone adaptation. We utilized homogenization theory to estimate bone tissue matrix strains surrounding osteocytes using two sets of models. The first set of models estimated the strain levels surrounding the lacunae and canaliculi, taking into account variations in lamellar properties. The second set estimated strain levels in the osteocyte and the surrounding matrix for different cellular mechanical properties. The results showed that the strain levels found in and surrounding osteocytes, 1700 to 2700 microstrain (denoted as μe; 1 μe =.0001% strain), were significantly greater than the trabecular tissue level strains of {1325 μe, 287 μe, 87 μe} used for model input. Variation in lamellar properties did not affect strain levels, except at lamellar boundaries. Strain in and surrounding the osteocyte was not significantly affected by cellular stiffness ranging between 28 and 28, 000 Pascals (Pa). Strain levels surrounding lacunae and canaliculi were approximately equivalent.     

A finite element model predicts the mechanotransduction response of tendon cells to cyclic tensile loading

The importance of fluid-flow-induced shear stress and matrix-induced cell deformation in transmitting the global tendon load into a cellular mechanotransduction response is yet to be determined. A multiscale computational tendon model composed of both matrix and fluid phases was created to examine how global tendon loading may affect fluid-flow-induced shear stresses and membrane strains at the cellular level. The model was then used to develop a quantitative experiment to help understand the roles of membrane strains and fluid-induced shear stresses on the biological response of individual cells. The model was able to predict the global response of tendon to applied strain (stress, fluid exudation), as well as the associated cellular response of increased fluid-flow-induced shear stress with strain rate and matrix-induced cell deformation with strain amplitude. The model analysis, combined with the experimental results, demonstrated that both strain rate and strain amplitude are able to independently alter rat interstitial collagenase gene expression through increases in fluid-flow-induced shear stress and matrix-induced cell deformation, respectively.     

Strain transfer through the aortic valve

The complex structural organization of the aortic valve (AV) extracellular matrix (ECM) enables large and highly nonlinear tissue level deformations. The collagen and elastin (elastic) fibers within the ECM form an interconnected fibrous network (FN) and are known to be the main load-bearing elements of the AV matrix. The role of the FN in enabling deformation has been investigated and documented. However, there is little data on the correlation between tissue level and FN-level strains. Investigating this correlation will help establish the mode of strain transfer (affine or nonaffine) through the AV tissue as a key feature in microstructural modeling and will also help characterize the local FN deformation across the AV sample in response to applied tissue level strains. In this study, the correlation between applied strains at tissue level, macrostrains across the tissue surface, and local FN strains were investigated. Results showed that the FN strain distribution across AV samples was inhomogeneous and nonuniform, as well as anisotropic. There was no direct transfer of the deformation applied at tissue level to the fibrous network. Loading modes induced in the FN are different than those applied at the tissue as a result of different local strains in the valve layers. This nonuniformity of local strains induced internal shearing within the FN of the AV, possibly exposing the aortic valve interstitial cells (AVICs) to shear strains and stresses.     

Transfer of macroscale tissue strain to microscale cell regions in the deformed meniscus

Cells within fibrocartilaginous tissues, including chondrocytes and fibroblasts of the meniscus, ligament, and tendon, regulate cell biosynthesis in response to local mechanical stimuli. The processes by which an applied mechanical load is transferred through the extracellular matrix to the environment of a cell are not fully understood. To better understand the role of mechanics in controlling cell phenotype and biosynthetic activity, this study was conducted to measure strain at different length scales in tissue of the fibrocartilaginous meniscus of the knee joint, and to define a quantitative parameter that describes the strain transferred from the far-field tissue to a microenvironment surrounding a cell. Experiments were performed to apply a controlled uniaxial tensile deformation to explants of porcine meniscus containing live cells. Using texture correlation analyses of confocal microscopy images, two-dimensional Lagrangian and principal strains were measured at length scales representative of the tissue (macroscale) and microenvironment in the region of a cell (microscale) to yield a strain transfer ratio as a measure of median microscale to macroscale strain. The data demonstrate that principal strains at the microscale are coupled to and amplified from macroscale principal strains for a majority of cell microenvironments located across diverse microstructural regions, with average strain transfer ratios of 1.6 and 2.9 for the maximum and minimum principal strains, respectively. Lagrangian strain components calculated along the experimental axes of applied deformations exhibited considerable spatial heterogeneity and intersample variability, and suggest the existence of both strain amplification and attenuation. This feature is consistent with an in-plane rotation of the principal strain axes relative to the experimental axes at the microscale that may result from fiber sliding, fiber twisting, and fiber-matrix interactions that are believed to be important for regulating deformation in other fibrocartilaginous tissues. The findings for consistent amplification of macroscale to microscale principal strains suggest a coordinated pattern of strain transfer from applied deformation to the microscale environment of a cell that is largely independent of these microstructural features in the fibrocartilaginous meniscus.     

Altered Virulence of Vaccine Strains of Measles Virus after Prolonged Replication in Human Tissue

To understand the molecular determinants of measles virus (MV) virulence, we have used the SCID-hu thymus/liver xenograft model (SCID-hu thy/liv) in which in vivo MV virulence phenotypes are faithfully duplicated. Stromal epithelial and monocytic cells are infected by MV in thymus implants, and virulent strains induce massive thymocyte apoptosis, although thymocytes are not infected. To determine whether passage of an avirulent vaccine strain in human tissue increases virulence, we studied a virus isolated from thymic tissue 90 days after infection with the vaccine strain Moraten (pMor-1) and a virus isolated from an immunodeficient child with progressive vaccine-induced disease (Hu2). These viruses were compared to a minimally passaged wild-type Edmonston strain (Ed-wt) and the vaccine strain Moraten. pMor-1, Hu2, and Ed-wt displayed virulent phenotypes in thymic implants, with high levels of virus being detected by 3 days after infection (10(5.2), 10(2.8), and 10(3. 4), respectively) and maximal levels being detected between 7 and 14 days after infection. In contrast, Moraten required over 14 days to grow to detectable levels. pMor-1 produced the highest levels of virus throughout infection, suggesting thymic adaptation of this strain. Similar to other virulent strains, Ed-wt, Hu2, and pMor-1 caused a decrease in the number of viable thymocytes as assessed by trypan blue exclusion and fluorescence-activated cell sorter analysis. Thymic architecture was also disrupted by these strains. Sequence analysis of the hemagglutinin (H) and matrix (M) genes showed no common changes in Hu2 and pMor-1. M sequences were identical in pMor-1 and Mor and varied in H at amino acid 469 (threonine to alanine), a position near the base of propeller 4 in the propeller blade/stem model of H structure. Further study will provide insights into the determinants of virulence.     

Subject-specific multiscale modeling of aortic valve biomechanics

A Finite Element workflow for the multiscale analysis of the aortic valve biomechanics was developed and applied to three physiological anatomies with the aim of describing the aortic valve interstitial cells biomechanical milieu in physiological conditions, capturing the effect of subject-specific and leaflet-specific anatomical features from the organ down to the cell scale. A mixed approach was used to transfer organ-scale information down to the cell-scale. Displacement data from the organ model were used to impose kinematic boundary conditions to the tissue model, while stress data from the latter were used to impose loading boundary conditions to the cell level. Peak of radial leaflet strains was correlated with leaflet extent variability at the organ scale, while circumferential leaflet strains varied over a narrow range of values regardless of leaflet extent. The dependency of leaflet biomechanics on the leaflet-specific anatomy observed at the organ length-scale is reflected, and to some extent emphasized, into the results obtained at the lower length-scales. At the tissue length-scale, the peak diastolic circumferential and radial stresses computed in the fibrosa correlated with the leaflet surface area. At the cell length-scale, the difference between the strains in two main directions, and between the respective relationships with the specific leaflet anatomy, was even more evident; cell strains in the radial direction varied over a relatively wide range (\(0.36-0.87\)) with a strong correlation with the organ length-scale radial strain (\(R^{2}= 0.95\)); conversely, circumferential cell strains spanned a very narrow range (\(0.75-0.88\)) showing no correlation with the circumferential strain at the organ level (\(R^{2}= 0.02\)). Within the proposed simulation framework, being able to account for the actual anatomical features of the aortic valve leaflets allowed to gain insight into their effect on the structural mechanics of the leaflets at all length-scales, down to the cell scale.

     

Zonal changes in the three-dimensional morphology of the chondron under compression: the relationship among cellular, pericellular, and extracellular deformation in articular cartilage

The pericellular matrix (PCM) is a narrow region of tissue that completely surrounds chondrocytes in articular cartilage. Previous theoretical models of the "chondron" (the PCM with enclosed cells) suggest that the structure and properties of the PCM may significantly influence the mechanical environment of the chondrocyte. The objective of this study was to quantify changes in the three-dimensional (3D) morphology of the chondron in situ at different magnitudes of compression applied to the cartilage extracellular matrix. Fluorescence immunolabeling for type-VI collagen was used to identify the boundaries of the cell and PCM, and confocal microscopy was used to form 3D images of chondrons from superficial, middle, and deep zone cartilage in explants compressed to 0%, 10%, 30%, and 50% surface-to-surface strain. Lagrangian tissue strain, determined locally using texture correlation, was highly inhomogeneous and revealed depth-dependent compressive stiffness and Poisson's ratio of the extracellular matrix. Compression significantly decreased cell and chondron height and volume, depending on the zone and magnitude of compression. In the superficial zone, cellular-level strains were always lower than tissue-level strains. In the middle and deep zones, however, tissue strains below 25% were amplified at the cellular level, while tissue strains above 25% were decreased at the cellular level. These findings are consistent with previous theoretical models of the chondron, suggesting that the PCM can serve as either a protective layer for the chondrocyte or a transducer that amplifies strain, such that cellular-level strains are more homogenous throughout the tissue depth despite large inhomogeneities in local ECM strains.     

Strain transfer in the annulus fibrosus under applied flexion

A detailed understanding of the anatomical and mechanical environment in the intervertebral disc at the scale of the cell is necessary for the design of tissue engineering repair strategies and to elucidate the role of mechanical factors in pathology. The objective of this study was to measure and compare the macroscale to microscale strains in the outer annulus fibrosus in various cellular regions of intact discs over a range of applied flexion. Macroscale strains were measured on the annulus fibrosus surface, and contrasted to in situ microscale strains using novel confocal microscopy techniques for dual labeling of the cell and the extracellular matrix. Fiber oriented surface strains were significantly higher than in situ fiber strains, which implies a mechanism of load redistribution that minimizes strain along the fibers. Non-uniformity of the strains and matrix distortion occurred immediately and most interestingly varied little with increase in flexion (3–16°), suggesting that inter-fiber shear is important in the initial stages of strain redistribution. Fiber oriented intercellular strains were significantly larger and compressive compared to in situ strains in other regions of the extracellular matrix indicating that the mechanical environment in this region may be unique. Further examination of the structural morphology in this pericellular region is needed to fully understand the pathway of strain transfer from the tissue to the cell. This study provides new knowledge on the complex in situ micro-mechanical environment of the annulus fibrosus that is essential to understanding the mechanobiological behavior of this tissue.     

Tissue and cellular morphological changes in growth plate explants under compression

The mechanisms by which mechanical loading may alter bone development within growth plates are still poorly understood. However, several growth plate cell or tissue morphological parameters are associated with both normal and mechanically modulated bone growth rates. The aim of this study was to quantify in situ the three-dimensional morphology of growth plate explants under compression at both cell and tissue levels. Growth plates were dissected from ulnae of immature swine and tested under 15% compressive strain. Confocal microscopy was used to image fluorescently labeled chondrocytes in the three growth plate zones before and after compression. Quantitative morphological analyses at both cell (volume, surface area, sphericity, minor/major radii) and tissue (cell/matrix volume ratio) levels were performed. Greater chondrocyte bulk strains (volume decrease normalized to the initial cell volume) were found in the proliferative (35.4%) and hypertrophic (41.7%) zones, with lower chondrocyte bulk strains (24.7%) in the reserve zone. Following compression, the cell/matrix volume ratio decreased in the reserve and hypertrophic zones by 24.3% and 22.6%, respectively, whereas it increased by 35.9% in the proliferative zone. The 15% strain applied on growth plate explants revealed zone-dependent deformational states at both tissue and cell levels. Variations in the mechanical response of the chondrocytes from different zones could be related to significant inhomogeneities in growth plate zonal mechanical properties. The ability to obtain in situ cell morphometry and monitor the changes under compression will contribute to a better understanding of mechanisms through which abnormal growth can be triggered.     

Homogenization theory and digital imaging: A basis for studying the mechanics and design principles of bone tissue

Bone tissue is a complex multilevel composite which has the ability to sense ad respond to its mechanical environment. It is believed that bone cells called osteocytes within the bone matrix sense the mechanical environment and determine whether structural alterations are needed. At present it is not known, however, how loads are transferred from the whole bone level to cells. A computational procedure combining representative volume element (RVE) based homogenization theory with digital imaging is proposed to estimate strains at various levels of bone structure. Bone tissue structural organization and RVE based analysis are briefly reviewed. The digital image based computational procedure was applied to estimate strains in individual trabeculae (first-level microstructure). Homogenization analysis of an idealized model was used to estimate strains at one level of bone structure around osteocyte lacunae (second-level trabecular microstructure). The results showed that strain at one level of bone structure is amplified to a broad range at the next microstructural level. In one case, a zeor-level tensile principal strain of 495 muE engendered strains ranging between -1000 and 7000 muE in individual trabeculae (first-level microstructure). Subsequently, a first-level tensile principal strains of 1325 muE within an inidividual trabecula engendered strains ranging between 782 and 2530 muE around osteocyte lacunae. Lacunar orientation was found to influence strains around osteocyte lacunae much more than lacunar ellipticity. In conclusion, the computational procedure combining homogenization theory with digital imaging can proveide estimates of cell level strains within whole bones. Such results may be used to bridge experimental studies of bone adaptation at the whole bone and cell culture level. (c) 1994 John Wiley & Sons, Inc.     

Effects of non-linearity on cell–ECM interactions

Filamentous biopolymers such as F-actin, vimentin, fibrin and collagen that form networks within the cytoskeleton or the extracellular matrix have unusual rheological properties not present in most synthetic soft materials that are used as cell substrates or scaffolds for tissue engineering. Gels formed by purified filamentous biopolymers are often strain stiffening, with an elastic modulus that can increase an order of magnitude at moderate strains that are relevant to cell and tissue deformation in vivo. This review summarizes some experimental studies of non-linear rheology in biopolymer gels, discusses possible molecular mechanisms that account for strain stiffening, and explores the possible relevance of non-linear rheology to the interactions between cell and extracellular matrices.     

Macro- to Microscale Strain Transfer in Fibrous Tissues is Heterogeneous and Tissue-Specific

Mechanical deformation applied at the joint or tissue level is transmitted through the macroscale extracellular matrix to the microscale local matrix, where it is transduced to cells within these tissues and modulates tissue growth, maintenance, and repair. The objective of this study was to investigate how applied tissue strain is transferred through the local matrix to the cell and nucleus in meniscus, tendon, and the annulus fibrosus, as well as in stem cell-seeded scaffolds engineered to reproduce the organized microstructure of these native tissues. To carry out this study, we developed a custom confocal microscope-mounted tensile testing device and simultaneously monitored strain across multiple length scales. Results showed that mean strain was heterogeneous and significantly attenuated, but coordinated, at the local matrix level in native tissues (35–70% strain attenuation). Conversely, freshly seeded scaffolds exhibited very direct and uniform strain transfer from the tissue to the local matrix level (15–25% strain attenuation). In addition, strain transfer from local matrix to cells and nuclei was dependent on fiber orientation and tissue type. Histological analysis suggested that different domains exist within these fibrous tissues, with most of the tissue being fibrous, characterized by an aligned collagen structure and elongated cells, and other regions being proteoglycan (PG)-rich, characterized by a dense accumulation of PGs and rounder cells. In meniscus, the observed heterogeneity in strain transfer correlated strongly with cellular morphology, where rounder cells located in PG-rich microdomains were shielded from deformation, while elongated cells in fibrous microdomains deformed readily. Collectively, these findings suggest that different tissues utilize distinct strain-attenuating mechanisms according to their unique structure and cellular phenotype, and these differences likely alter the local biologic response of such tissues and constructs in response to mechanical perturbation.     

Effect of Bacillus subtilis BsuM restriction-modification on plasmid transfer by polyethylene glycol-induced protoplast fusion

Polyethylene glycol (PEG)-induced cell fusion is a promising method to transfer larger DNA from one cell to another than conventional genetic DNA transfer systems. The laboratory strain Bacillus subtilis 168 contains a restriction (R) and modification (M) system, BsuM, which recognizes the sequence 5'-CTCGAG-3'. To study whether the BsuM system affects DNA transfer by the PEG-induced cell fusion between R(+)M(+) and R(-)M(-) strains, we examined transfer of plasmids pHV33 and pLS32neo carrying no and eight BsuM sites, respectively. It was shown that although the transfer of pLS32neo but not pHV33 from the R(-)M(-) to R(+)M(+) cells was severely restricted, significant levels of transfer of both plasmids from the R(+)M(+) to R(-)M(-) cells were observed. The latter result shows that the chromosomal DNA in the R(-)M(-) cell used as the recipient partially survived restriction from the donor R(+)M(+) cell, indicating that the BsuM R(-)M(-) strain is useful as a host for accepting DNA from cells carrying a restriction system(s). Two such examples were manifested for plasmid transfer from Bacillus circulans and Bacillus stearothermophilus strains to a BsuM-deficient mutant, B. subtilis RM125.     

Macro- to Microscale Strain Transfer in Fibrous Tissues is Heterogeneous and Tissue-Specific

Mechanical deformation applied at the joint or tissue level is transmitted through the macroscale extracellular matrix to the microscale local matrix, where it is transduced to cells within these tissues and modulates tissue growth, maintenance, and repair. The objective of this study was to investigate how applied tissue strain is transferred through the local matrix to the cell and nucleus in meniscus, tendon, and the annulus fibrosus, as well as in stem cell-seeded scaffolds engineered to reproduce the organized microstructure of these native tissues. To carry out this study, we developed a custom confocal microscope-mounted tensile testing device and simultaneously monitored strain across multiple length scales. Results showed that mean strain was heterogeneous and significantly attenuated, but coordinated, at the local matrix level in native tissues (35–70% strain attenuation). Conversely, freshly seeded scaffolds exhibited very direct and uniform strain transfer from the tissue to the local matrix level (15–25% strain attenuation). In addition, strain transfer from local matrix to cells and nuclei was dependent on fiber orientation and tissue type. Histological analysis suggested that different domains exist within these fibrous tissues, with most of the tissue being fibrous, characterized by an aligned collagen structure and elongated cells, and other regions being proteoglycan (PG)-rich, characterized by a dense accumulation of PGs and rounder cells. In meniscus, the observed heterogeneity in strain transfer correlated strongly with cellular morphology, where rounder cells located in PG-rich microdomains were shielded from deformation, while elongated cells in fibrous microdomains deformed readily. Collectively, these findings suggest that different tissues utilize distinct strain-attenuating mechanisms according to their unique structure and cellular phenotype, and these differences likely alter the local biologic response of such tissues and constructs in response to mechanical perturbation.     

Biochemical properties of Fusobacterium naviforme and phenotypically similar isolates

Three strains which resemble the type strain of Fusobacterium naviforme (ATCC 25832) by morphological and physiological criteria were isolated from human clinical specimens. All were non-fermentative, produced indole and, in common with other members of the genus Fusobacterium , butyrate was a major end-product of metabolism. Glutamate dehydrogenase and 2-oxoglutarate reductase were present in both taxa, but the enzymes of the test strains migrated to only about half the distance of that of strain ATCC 25832. The latter contained meso -diaminopimelic acid as its peptidoglycan dibasic amino acid whereas the test strains possessed meso -lanthionine. The wide divergence in DNA base composition between strain ATCC 25832 (49 mol% G + C) and the clinical isolates ( ca 30–31 mol% G + C) was reflected in their low DNA-DNA homology ( ca 5–15%). The present study therefore revealed major differences between F. naviforme (ATCC 25832) and the new isolates and indicate that the latter may belong to a hitherto undescribed taxon within the genus Fusobacterium.     

Multiscale mechanisms of tendon fatigue damage progression and severity are strain and cycle dependent

Tendinopathies are common chronic injuries that occur when damage accumulation caused by sub-rupture fatigue loading outpaces repair. Studies have linked fatigue loading with various mechanical, structural, and biological changes associated with pathology. However, the multiscale progression of damage accumulation with respect to area, severity and the distinct contributions of strain level and number of cycles has not been fully elucidated. The objective of this study was to investigate multiscale mechanisms underlying fatigue damage accumulation and their effect on the cellular environment. Using an in situ model in rat tail tendon (RTT), fatigue loading was applied at various strains and cycle numbers to induce fatigue damage. Pre- and post- fatigue diagnostic mechanical testing, second harmonic generation (SHG) imaging, and transmission electron microscope (TEM) imaging were used to investigate extracellular and cellular damage modes at multiple scales. Fatigue loading at strains at or below 1.0% resulted in no significant changes in SHG damage area or severity and no changes in collagen fibril or cell morphology compared with controls. Fatigue loading at strains above 1.5% resulted in greater mechanical changes correlated with increased damage area measured by SHG and collagenous damage observed by TEM. Increased cycles at high strain further altered mechanical properties, increased structural damage severity (but not area), and altered TEM collagen rupture patterns. Cell morphology was similarly progressively affected with increased strain and cycle number. These damage mechanisms that may trigger degenerative changes characteristic of tendinopathy could be targeted as a part of prevention or therapy.     

Finite element modeling of 3D human mesenchymal stem cell-seeded collagen matrices exposed to tensile strain

The use of human mesenchymal stem cells (hMSCs) in tissue engineering is attractive due to their ability to extensively self-replicate and differentiate into a multitude of cell lineages. It has been experimentally established that hMSCs are influenced by chemical and mechanical signals. However, the combined chemical and mechanical in vitro culture conditions that lead to functional tissue require greater understanding. In this study, finite element models were created to evaluate the local loading conditions on bone marrow-derived hMSCs seeded in three-dimensional collagen matrices exposed to cyclic tensile strain. Mechanical property and geometry data used in the models were obtained experimentally from a previous study in our laboratory and from mechanical testing. Eight finite element models were created to simulate three-dimensional hMSC-seeded collagen matrices exposed to different levels of cyclic tensile strain (10% and 12%), culture media (complete growth and osteogenic differentiating), and durations of culture (7 and 14 days). Through finite element analysis, it was determined that globally applied uniaxial tensile strains of 10% and 12% resulted in local strains up to 18.3% and 21.8%, respectively. Model results were also compared to experimental studies in an attempt to explain observed differences between hMSC response to 10% and 12% cyclic tensile strain.     

Deoxyribonucleoside triphosphate pool levels in three cell strains of human chromosome instability syndromes: ataxia telangiectasia (GM2052), Bloom's syndrome (GM1492), and Fanconi's anemia (GM368)

Deoxyribonucleoside triphosphate (dNTP) pool sizes were determined in cell strains derived from patients with the genetic diseases ataxia telangiectasia (GM2052), Bloom's syndrome (GM1492), and Fanconi's anemia (GM368), and were compared to the dNTP pools in a normal human fibroblast cell strain (253/79). In addition, the effect of deoxythymidine on both dNTP pool levels and cell growth was examined. The three mutant cell strains differed only slightly from the normal cell strain. The cellular characteristics of the cell strains, such as chromosome instability, are apparently not an effect of dNTP pool imbalance.     

Multilevel finite element modeling for the prediction of local cellular deformation in bone

The underlying mechanisms by which bone cells respond to mechanical stimuli or how mechanical loads act on osteocytes housed in lacunae in bone are not well understood. In this study, a multilevel finite element (FE) approach is applied to predict local cell deformations in bone tissue. The local structure of the matrix dictates the local mechanical environment of an osteocyte. Cell deformations are predicted from detailed linear FE analysis of the microstructure, consisting of an arrangement of cells embedded in bone matrix material. This work has related the loads applied to a whole femur during the stance phase of the gait cycle to the strain of a single lacuna and of canaliculi. The predicted bone matrix strains around osteocyte lacunae and canaliculi were nonuniform and differed significantly from the macroscopically measured strains. Peak stresses and strains in the walls of the lacuna were up to six times those in the bulk extracellular matrix. Significant strain concentrations were observed at sites where the process meets the cell body.