Search

Determining tissue biomechanical material properties from mechanical test data is frequently required in a variety of applications. However, the validity of the resulting constitutive model parameters is the subject of debate in the field. Parameter optimization in tissue mechanics often comes down to the “identifiability” or “uniqueness” of constitutive model parameters; however, despite advances in formulating complex constitutive relations and many classic and creative curve-fitting approaches, there is currently no accessible framework to study the identifiability of tissue material parameters. Our objective was to assess the identifiability of material parameters for established constitutive models of fiber-reinforced soft tissues, biomaterials, and tissue-engineered constructs and establish a generalizable procedure for other applications. To do so, we generated synthetic experimental data by simulating uniaxial tension and compression tests, commonly used in biomechanics. We then fit this data using a multi-start optimization technique based on the nonlinear least-squares method with multiple initial parameter guesses. We considered tendon and sclera as example tissues, using constitutive models that describe these fiber-reinforced tissues. We demonstrated that not all the model parameters of these constitutive models were identifiable from uniaxial mechanical tests, despite achieving virtually identical fits to the stress-stretch response. We further show that when the lateral strain was considered as an additional fitting criterion, more parameters are identifiable, but some remain unidentified. This work provides a practical approach for addressing parameter identifiability in tissue mechanics.

The knee menisci are comprised of two orthogonal collagenous networks – circumferential and radial – that combine to enable efficient load bearing by the tissue in adults. Here, we assessed how the structural and functional characteristics of these networks developed over the course of skeletal maturation and determined the role of these fiber networks in defect tolerance with tissue injury. Imaging of the radial tie fiber (RTF) collagen structure in medial bovine menisci from fetal, juvenile, and adult specimens showed increasing heterogeneity, anisotropy, thickness, and density with skeletal development. Mechanical analysis showed that the tensile modulus in the radial direction did not change with skeletal development, though the resilience (in the radial direction) increased and the tolerance to defects in the circumferential direction decreased, in adult compared to fetal tissues. This loss of defect tolerance correlated with increased order in the RTF network in adult tissue. These data provide new insights into the role of the radial fiber network in meniscus function, will lead to improved clinical decision-making in the presence of a tear, and may improve engineering efforts to reproduce this critical load-bearing structure in the knee.

The meniscus plays a central load bearing role in the knee joint. Unfortunately, meniscus injury is common and can lead to joint degeneration and osteoarthritis. In small animal models, progressive degenerative changes occur with unloading of the meniscus via destabilization of the medial meniscus (DMM). However, few large animal models of DMM exist and the joint-wide initiation of disease has not yet been defined in these models. Thus, the goal of this study was to develop and validate a large-animal model of surgically-induced destabilization of the medial meniscus and to use multi-modal (mechanical, histological, and MRI) and multi-scale (joint to tissue level) quantitative measures to evaluate degeneration in both the meniscus and cartilage. DMM was achieved using an arthroscopic approach in thirteen Yucatan minipigs. One month after DMM, joint contact area decreased and peak pressure increased, indicating altered load transmission as a result of meniscus destabilization. By three months, the joint had adapted to the injury and load transmission patterns were restored to baseline, likely due to the formation and maturation of a fibrovascular scar at the anterior aspect of the meniscus. Despite this, we found a decrease in the indentation modulus of the tibial cartilage and an increase in cartilage histopathology scores at one month compared to Sham operated animals; these deleterious changes persisted through three months. Over this same time course, meniscus remodeling was evident through decreased proteoglycan staining in DMM compared to Sham menisci at both one and three months. These findings support that arthroscopic DMM results in joint degeneration in the Yucatan minipig and provides a new large animal test bed in which to evaluate therapeutics and interventions to treat post-traumatic osteoarthritis (PTOA) that originates from meniscal injury.

In fibrous tissues, prestressed boundary constraints at bone interfaces instil residual strain throughout the tissue, even when unloaded. For example, internal swelling pressures in the central nucleus pulposus of the intervertebral disc generate prestrain in the outer annulus fibrosus. With injury and depressurization, these residual strains are lost. Here we show that the loss of residual strains in the intervertebral disc alters the microenvironment and instigates aberrant tissue remodelling and the adoption of atypical cellular phenotypes. By using puncture surgery of the annulus fibrosus in rabbits, ex vivo puncture experiments and electrospun nanofibrous scaffolds recapitulating these evolving boundary constraints, we show that the loss of residual strain promotes short-term apoptosis and the emergence of a fibrotic phenotype. We also show that local fibre organization and cellular contractility mediate this process and that the aberrant cellular changes could be abrogated by targeting the cell-mechanosensing machinery with small molecules. Our findings indicate that injury to dense connective tissues under prestrain alters boundary constraints and residual strain; this leads to aberrant mechanosensing, which in turn promotes disease progression.

Tendinopathy, degeneration of the tendon that leads to pain and dysfunction, is common in both sports and occupational settings, but multi-scale mechanisms for tendinopathy are still unknown. We recently showed that micro-scale sliding (shear) is responsible for both load transfer and damage mechanisms in the rat tail tendon; however, the rat tail tendon is a specialized non-load-bearing tendon, and thus the load transfer and damage mechanisms are still unknown for load-bearing tendons. The objective of this study was to investigate the load transfer and damage mechanisms of load-bearing tendons using the rat plantaris tendon. We demonstrated that micro-scale sliding is a key component for both mechanisms in the plantaris tendon, similar to the tail tendon. Namely, the micro-scale sliding was correlated with applied strain, demonstrating that load was transferred via micro-scale sliding in the plantaris and tail tendons. In addition, while the micro-scale strain fully recovered, the micro-scale sliding was non-recoverable and strain-dependent, and correlated with tissue-scale mechanical parameters. When the applied strain was normalized, the % magnitudes of non-recoverable sliding was similar between the plantaris and tail tendons. Statement of clinical significance: Understanding the mechanisms responsible for the pathogenesis and progression of tendinopathy can improve prevention and rehabilitation strategies and guide therapies and the design of engineered constructs. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 37:1827-1837, 2019.

Degeneration of the intervertebral discs is a progressive cascade of cellular, compositional and structural changes that is frequently associated with low back pain. As the first signs of disc degeneration typically arise in the disc’s central nucleus pulposus (NP), augmentation of the NP via hydrogel injection represents a promising strategy to treat early to mid-stage degeneration. The purpose of this study was to establish the translational feasibility of a triple interpenetrating network hydrogel composed of dextran, chitosan, and teleostean (DCT) for augmentation of the degenerative NP in a preclinical goat model. Ex vivo injection of the DCT hydrogel into degenerated goat lumbar motion segments restored range of motion and neutral zone modulus towards physiologic values. To facilitate non-invasive assessment of hydrogel delivery and distribution, zirconia nanoparticles were added to make the hydrogel radiopaque. Importantly, the addition of zirconia did not negatively impact viability or matrix producing capacity of goat mesenchymal stem cells or NP cells seeded within the hydrogel in vitro. In vivo studies demonstrated that the radiopaque DCT hydrogel was successfully delivered to degenerated goat lumbar intervertebral discs, where it was distributed throughout both the NP and annulus fibrosus, and that the hydrogel remained contained within the disc space for two weeks without evidence of extrusion. These results demonstrate the translational potential of this hydrogel for functional regeneration of degenerate intervertebral discs. Statement of Significance The results of this work demonstrate that a radiopaque hydrogel is capable of normalizing the mechanical function of the degenerative disc, is supportive of disc cell and mesenchymal stem cell viability and matrix production, and can be maintained in the disc space without extrusion following intradiscal delivery in a preclinical large animal model. These results support evaluation of this hydrogel as a minimally invasive disc therapeutic in long-term preclinical studies as a precursor to future clinical application in patients with disc degeneration and low back pain.

A buffer solution is often used to maintain tissue hydration during mechanical testing. The most commonly used buffer solution is a physiological concentration of phosphate buffered saline (PBS); however, PBS increases the tissue’s water content and decreases its tensile stiffness. In addition, solutes from the buffer can diffuse into the tissue and interact with its structure and mechanics. These bathing solution effects can confound the outcome and interpretation of mechanical tests. Potential bathing solution artifacts, including solute diffusion and the effect on mechanical properties, are not well understood. The objective of this study was to measure the effects of long-term exposure of rat tail tendon fascicles to several concentrations (0.9% to 25%) of NaCl, sucrose, polyethylene glycol (PEG), and SPEG (NaCl + PEG) solutions on water content, solute diffusion, and mechanical properties. We found that with an increase in solute concentration the apparent water content decreased for all solution types. Solutes diffused into the tissue for NaCl and sucrose, however, no solute diffusion was observed for PEG or SPEG. The mechanical properties changed for both of NaCl solutions, in particular after long-term (8 hr) incubation the modulus and equilibrium stress decreased compared to short-term (15 min) for 25% NaCl, and the cross sectional area increased for 0.9% NaCl. However, the mechanical properties were unchanged for both PEG and SPEG except for minor alterations in stress relaxation parameters. This study shows that NaCl and sucrose buffer solutions are not suitable for long-term mechanical tests. We therefore propose using PEG or SPEG as alternative buffer solutions that after long-term incubation can maintain tissue hydration without solute diffusion and produce a consistent mechanical response.

Collagen fibrils in tendon are believed to be discontinuous and transfer tensile loads through shear forces generated during interfibrillar sliding. However, the structures that transmit these interfibrillar forces are unknown. Various extrafibrillar tissue components (e.g., glycosaminoglycans, collagens XII and XIV) have been suggested to transmit interfibrillar loads by bridging collagen fibrils. Alternatively, collagen fibrils may interact directly through physical fusions and interfibrillar branching. The objective of this study was to test whether extrafibrillar proteins are necessary to transmit load between collagen fibrils or if interfibrillar load transfer is accomplished directly by the fibrils themselves. Trypsin digestions were used to remove a broad spectrum of extrafibrillar proteins and measure their contribution to the multiscale mechanics of rat tail tendon fascicles. Additionally, images obtained from serial block-face scanning electron microscopy were used to determine the three-dimensional fibrillar organization in tendon fascicles and identify any potential interfibrillar interactions. While trypsin successfully removed several extrafibrillar tissue components, there was no change in the macroscale fascicle mechanics or fibril:tissue strain ratio. Furthermore, the imaging data suggested that a network of smaller diameter fibrils (

Geometry is an important indicator of disc mechanical function and degeneration. While the geometry and associated degenerative changes in the nucleus pulposus and the annulus fibrosus are well-defined, the geometry of the cartilage endplate (CEP) and its relationship to disc degeneration are unknown. The objectives of this study were to quantify CEP geometry in three dimensions using an MRI FLASH imaging sequence and evaluate relationships between CEP geometry and age, degeneration, spinal level, and overall disc geometry. To do so, we assessed the MRI-based measurements for accuracy and repeatability. Next, we measured CEP geometry across a larger sample set and correlated CEP geometric parameters to age, disc degeneration, level, and disc geometry. The MRI-based measures resulted in thicknesses (0.3-1 mm) that are comparable to prior measurements of CEP thickness. CEP thickness was greatest at the anterior/posterior (A/P) margins and smallest in the center. The CEP A/P thickness, axial area, and lateral width decreased with age but were not related to disc degeneration. Age-related, but not degeneration-related, changes in geometry suggest that the CEP may not follow the progression of disc degeneration. Ultimately, if the CEP undergoes significant geometric changes with aging and if these can be related to low back pain, a clinically feasible translation of the FLASH MRI-based measurement of CEP geometry presented in this study may prove a useful diagnostic tool. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1410-1417, 2016.

The inter-lamellar matrix (ILM) has an average thickness of less than 30 µm and lies between adjacent lamellae in the annulus fibrosus (AF). The microstructure and composition of the ILM have been studied in various anatomic regions of the disc; however, their contribution to AF mechanical properties and structural integrity is unknown. It was suggested that the ILM components, mainly elastic fibers and cross-bridges, play a role in providing mechanical integrity of the AF. Therefore, the manner in which they respond to different loadings and stabilize adjacent lamellae structure will influence AF tear formation and subsequent herniation. This review paper summarizes the composition, microstructure, and potential role of the ILM in the progression of disc herniation, clarifies the micromechanical properties of the ILM, and proposes critical areas for future studies. There are a number of unknown characteristics of the ILM, such as its mechanical role, impact on AF integrity, and ultrastructure of elastic fibers at the ILM-lamella boundary. Determining these characteristics will provide important information for tissue engineering, repair strategies, and the development of more-physiological computational models to study the initiation and propagation of AF tears that lead to herniation and degeneration. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1307-1315, 2016.

Tendon’s viscoelastic behaviors are important to the mechanical function and mechanobiology. When loaded in longitudinal tension, tendons often have a very large Poisson’s ratio (ν > 2) that exceeds the limit of incompressibility for isotropic material (ν = 0.5), indicating that tendon experiences volume loss, inducing poroelastic fluid exudation in the transverse direction. Therefore, transverse poroelasticity is an important contributor to tendon material behavior. Tendon hydraulic permeability which is required to evaluate the fluid flow contribution to viscoelasticity, is mostly unavailable, and where available, varies by several orders of magnitude. In this manuscript, we quantified the transverse poroelastic material parameters of rat tail tendon fascicles by conducting transverse osmotic loading experiments, in both tension and compression. We used a multi-start optimization method to evaluate the parameters based on biphasic finite element modeling. Our tendon samples had a transverse hydraulic permeability of 10(−4) to 10(−5) mm(4) (Ns)(−1) and showed a significant tension-compression nonlinearity in the transverse direction. Further, using these results, we predict hydraulic permeability during longitudinal (fiber-aligned) tensile loading, and the spatial distribution of fluid flow during osmotic loading. These results reveal novel aspects of tendon mechanics and can be used to study the physiomechanical response of tendon in response to mechanical loading.

Here we develop a three-dimensional analytic model for MR image contrast of collagen lamellae in the annulus fibrosus of the intervertebral disc of the spine, based on the dependence of the MRI signal on collagen fiber orientation. High-resolution MRI scans were performed at 1.5 and 7 T on intact whole disc specimens from ovine, bovine, and human spines. An analytic model that approximates the three-dimensional curvature of the disc lamellae was developed to explain inter-lamellar contrast and intensity variations in the annulus. The model is based on the known anisotropic dipolar relaxation of water in tissues with ordered collagen. Simulated MRI data were generated that reproduced many features of the actual MRI data. The calculated inter-lamellar image contrast demonstrated a strong dependence on the collagen fiber angle and on the circumferential location within the annulus. This analytic model may be useful for interpreting MR images of the disc and for predicting experimental conditions that will optimize MR image contrast in the annulus fibrosus.

Tissue strain is an important indicator of mechanical function, but is difficult to noninvasively measure in the intervertebral disc. The objective of this study was to generate a disc strain template, a 3D average of disc strain, of a group of human L4-L5 discs loaded in axial compression. To do so, magnetic resonance images of uncompressed discs were used to create an average disc shape. Next, the strain tensors were calculated pixel-wise by using a previously developed registration algorithm. Individual disc strain tensor components were then transformed to the template space and averaged to create the disc strain template. The strain template reduced individual variability while highlighting group trends. For example, higher axial and circumferential strains were present in the lateral and posterolateral regions of the disc, which may lead to annular tears. This quantification of group-level trends in local 3D strain is a significant step forward in the study of disc biomechanics. These trends were compared to a finite element model that had been previously validated against the disc-level mechanical response. Depending on the strain component, 81-99% of the regions within the finite element model had calculated strains within one standard deviation of the template strain results. The template creation technique provides a new measurement technique useful for a wide range of studies, including more complex loading conditions, the effect of disc pathologies and degeneration, damage mechanisms, and design and evaluation of treatments. © 2015 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1264-1273, 2016.

Viscoelastic mechanical properties are frequently altered after tendon injuries and during recovery. Therefore, non-invasive measurements of shear viscoelastic properties may help evaluate tendon recovery and compare the effectiveness of different therapies. The objectives of this study were to describe an elastography method for measuring localized viscoelastic properties of tendons and to discuss the initial results in healthy and injured human Achilles and semitendinosus tendons. The technique used an external actuator to generate the shear waves in the tendon at different frequencies and plane wave imaging to measure shear wave displacements. For each of the excitation frequencies, maps of direction-specific wave speeds were calculated using local frequency estimation. Maps of viscoelastic properties were obtained using a pixel-wise curve fit of wave speed and frequency. The method was validated by comparing measurements of wave speed in agarose gels with those obtained using magnetic resonance elastography. Measurements in human healthy Achilles tendons revealed a pronounced increase in wave speed as a function of frequency, which highlights the importance of tendon viscoelasticity. Additionally, the viscoelastic properties of the Achilles tendon were larger than those reported for other tissues. Measurements in a tendinopathic Achilles tendon indicated that it is feasible to quantify local viscoelastic properties. Similarly, measurement in the semitendinosus tendon revealed substantial differences in viscoelastic properties between the healthy and contralateral tendons. Consequently, this technique has the potential to evaluate localized changes in tendon viscoelastic properties caused by injury and during recovery in a clinical setting.

Changes in tendon viscoelastic properties are observed after injuries and during healing as a product of altered composition and structure. Continuous Shear Wave Elastography is a new technique measuring viscoelastic properties of soft tissues using external shear waves. Tendon has not been studied with this technique, therefore, the aims of this study were to establish the range of shear and viscosity moduli in healthy Achilles tendons, determine bilateral differences of these parameters and explore correlations of viscoelasticity to plantar flexion strength and tendon area. Continuous Shear Wave Elastography was performed over the free portion of both Achilles tendons from 29 subjects. Isometric plantar flexion strength and cross sectional area were measured. The average shear and viscous moduli was 83.2kPa and 141.0Pa-s, respectively. No correlations existed between the shear or viscous modulus and area or strength. This indicates that viscoelastic properties can be considered novel, independent biomarkers. The shear and viscosity moduli were bilaterally equivalent (p=0.013,0.017) which allows determining pathologies through side-to-side deviations. The average bilateral coefficient of variation was 7.2% and 9.4% for shear and viscosity modulus, respectively. The viscoelastic properties of the Achilles tendon may provide an unbiased, non-subjective rating system of tendon recovery and optimizing treatment strategies.

Degeneration of the intervertebral discs is strongly implicated as a cause of low back pain. Since current treatments for discogenic low back pain show poor long-term efficacy, a number of new biological strategies are being pursued. For such therapies to succeed, it is critical that they be validated in conditions that mimic the unique biochemical microenvironment of the nucleus pulposus (NP), which include low oxygen tension. Therefore, the objective of this study was to investigate the effects of oxygen tension on NP cell functional extracellular matrix elaboration in 3D culture. Bovine NP cells were encapsulated in agarose constructs and cultured for 14 or 42 days in either 20% or 2% oxygen in defined media containing transforming growth factor beta-3. At each time point, extracellular matrix composition, biomechanics, and mRNA expression of key phenotypic markers were evaluated. Results showed that while bulk mechanics and composition were largely independent of oxygen level, low oxygen promoted improved restoration of the NP phenotype, higher mRNA expression of extracellular matrix and NP specific markers, and more uniform matrix elaboration. These findings indicate that culture under physiological oxygen levels is an important consideration for successful development of cell and growth factor-based regenerative strategies for the disc.

Degradation of the nucleus pulposus (NP) is an early hallmark of intervertebral disc degeneration. The capacity for endogenous regeneration in the NP is limited due to the low cellularity and poor nutrient supply of this avascular tissue. Towards restoring the NP, a number of biomaterials have been explored for cell delivery. These materials must support the NP cell phenotype while promoting the elaboration of an NP-like extracellular matrix in the shortest possible time. Our previous work with chondrocytes and mesenchymal stem cells demonstrated that hydrogels based on hyaluronic acid (HA) are effective at promoting matrix production and the development of functional material properties. However, this material has not been evaluated in the context of NP cells. Therefore, to test this material for NP regeneration, bovine NP cells were encapsulated in 1% w/vol HA hydrogels at either a low seeding density (20 × 106 cells/ml) or a high seeding density (60 × 106 cells/ml), and constructs were cultured over an 8 week period. These engineered NP cell-laden HA hydrogels showed functional matrix accumulation, with increasing matrix content and mechanical properties with time in culture at both seeding densities. Furthermore, encapsulated cells showed NP-specific gene expression profiles that were significantly higher than expanded NP cells prior to encapsulation, suggesting a restoration of phenotype. Interestingly, these levels were higher at the lower seeding density compared to the higher seeding density. These findings support the use of HA-based hydrogels for NP tissue engineering and cellular therapies directed at restoration or replacement of the endogenous NP.

It is common in biomechanics to use previously frozen tissues, where it is assumed that the freeze-thaw process does not cause consequential mechanical or structural changes. We have recently quantified multiscale tendon mechanics and damage mechanisms using previously frozen tissue, where damage was defined as an irreversible change in the microstructure that alters the macroscopic mechanical parameters. Because freezing has been shown to alter tendon microstructures, the objective of this study was to determine if freezing alters tendon multiscale mechanics and damage mechanisms. Multiscale testing using a protocol that was designed to evaluate tendon damage (tensile stress-relaxation followed by unloaded recovery) was performed on fresh and previously frozen rat tail tendon fascicles. At both the fascicle and fibril levels, there was no difference between the fresh and frozen groups for any of the parameters, suggesting that there is no effect of freezing on tendon mechanics. After unloading, the microscale fibril strain fully recovered, and interfibrillar sliding only partially recovered, suggesting that the tendon damage is localized to the interfibrillar structures and that mechanisms of damage are the same in both fresh and previously frozen tendons.

Despite current knowledge of tendon structure, the fundamental deformation mechanisms underlying tendon mechanics and failure are unknown. We recently showed that a shear lag model, which explicitly assumed plastic interfibrillar load transfer between discontinuous fibrils, could explain the multiscale fascicle mechanics, suggesting that fascicle yielding is due to plastic deformation of the interfibrillar matrix. However, it is unclear whether alternative physical mechanisms, such as elastic interfibrillar deformation or fibril yielding, also contribute to fascicle mechanical behavior. The objective of the current work was to determine if plasticity of the interfibrillar matrix is uniquely capable of explaining the multiscale mechanics of tendon fascicles including the tissue post-yield behavior. This was examined by comparing the predictions of a continuous fibril model and three separate shear lag models incorporating an elastic, plastic, or elastoplastic interfibrillar matrix with multiscale experimental data. The predicted effects of fibril yielding on each of these models were also considered. The results demonstrated that neither the continuous fibril model nor the elastic shear lag model can successfully predict the experimental data, even if fibril yielding is included. Only the plastic or elastoplastic shear lag models were capable of reproducing the multiscale tendon fascicle mechanics. Differences between these two models were small, although the elastoplastic model did improve the fit of the experimental data at low applied tissue strains. These findings suggest that while interfibrillar elasticity contributes to the initial stress response, plastic deformation of the interfibrillar matrix is responsible for tendon fascicle post-yield behavior. This information sheds light on the physical processes underlying tendon failure, which is essential to improve our understanding of tissue pathology and guide the development of successful repair.

Intervertebral disc mechanics are affected by both disc shape and disc degeneration, which in turn each affect the other; disc mechanics additionally have a role in the etiology of disc degeneration. Finite element analysis (FEA) is a favored tool to investigate these relationships, but limited data for intervertebral disc 3D shape has forced the use of simplified or single-subject geometries, with the effect of inter-individual shape variation investigated only in specialized studies. Similarly, most data on disc shape variation with degeneration is based on 2D mid-sagittal images, which incompletely define 3D shape changes. Therefore, the objective of this study was to quantify inter-individual disc shape variation in 3D, classify this variation into independently-occurring modes using a statistical shape model, and identify correlations between disc shape and degeneration. Three-dimensional disc shapes were obtained from MRI of 13 human male cadaver L3L4 discs. An average disc shape and four major modes of shape variation (representing 90% of the variance) were identified. The first mode represented disc axial area and was significantly correlated to degeneration (R(2)=0.44), indicating larger axial area in degenerate discs. Disc height variation occurred in three distinct modes, each also involving non-height variation. The statistical shape model provides an average L3L4 disc shape for FEA that is fully defined in 3D, and makes it convenient to generate a set of shapes with which to represent aggregate inter-individual variation. Degeneration grade-specific shapes can also be generated. To facilitate application, the model is included in this paper׳s supplemental content.

Intervertebral disc degeneration has been implicated in the etiology of low back pain; however, the current surgical strategies for treating symptomatic disc disease are limited. A variety of materials have been developed to replace disc components, including the nucleus pulposus (NP), the annulus fibrosus (AF) and their combination into disc-like engineered constructs. We have previously shown that layers of electrospun poly(ε-caprolactone) scaffold, mimicking the hierarchical organization of the native AF, can achieve functional parity with native tissue. Likewise, we have combined these structures with cell-seeded hydrogels (as an NP replacement) to form disc-like angle-ply structures (DAPS). The objective of this study was to develop a model for the evaluation of DAPS in vivo. Through a series of studies, we developed a surgical approach to replace the rat caudal disc with an acellular DAPS and then stabilized the motion segment via external fixation. We then optimized cell infiltration into DAPS by including sacrificial poly(ethylene oxide) layers interspersed throughout the angle-ply structure. Our findings illustrate that DAPS are stable in the caudal spine, are infiltrated by cells from the peri-implant space and that infiltration is expedited by providing additional routes for cell migration. These findings establish a new in vivo platform in which to evaluate and optimize the design of functional disc replacements.

The aim of functional tissue engineering is to repair and replace tissues that have a biomechanical function, i.e., connective orthopaedic tissues. To do this, it is necessary to have accurate benchmarks for the elastic, permeability, and swelling (i.e., biphasic-swelling) properties of native tissues. However, in the case of the intervertebral disc, the biphasic-swelling properties of individual tissues reported in the literature exhibit great variation and even span several orders of magnitude. This variation is probably caused by differences in the testing protocols and the constitutive models used to analyze the data. Therefore, the objective of this study was to measure the human lumbar disc annulus fibrosus (AF), nucleus pulposus (NP), and cartilaginous endplates (CEP) biphasic-swelling properties using a consistent experimental protocol and analyses. The testing protocol was composed of a swelling period followed by multiple confined compression ramps. To analyze the confined compression data, the tissues were modeled using a biphasic-swelling model, which augments the standard biphasic model through the addition of a deformation-dependent osmotic pressure term. This model allows considering the swelling deformations and the contribution of osmotic pressure in the analysis of the experimental data. The swelling stretch was not different between the disc regions (AF: 1.28±0.16; NP: 1.73±0.74; CEP: 1.29±0.26), with a total average of 1.42. The aggregate modulus (Ha) of the extra-fibrillar matrix was higher in the CEP (390kPa) compared to the NP (100kPa) or AF (30kPa). The permeability was very different across tissue regions, with the AF permeability (64 E(-16)m(4)/Ns) higher than the NP and CEP (~5.5 E(-16)m(4)/Ns). Additionally, a normalized time-constant (3000s) for the stress relaxation was similar for all the disc tissues. The properties measured in this study are important as benchmarks for tissue engineering and for modeling the disc's mechanical behavior and transport.

The mechanical properties of the human supraspinatus tendon (SST) are highly heterogeneous and may reflect an important adaptive response to its complex, multiaxial loading environment. However, these functional properties are associated with a location-dependent structure and composition that have not been fully elucidated. Therefore, the objective of this study was to determine the concentrations of types I, II and III collagen in six distinct regions of the SST and compare changes in collagen concentration across regions with local changes in mechanical properties. We hypothesized that type I collagen content would be high throughout the tendon, type II collagen would be restricted to regions of compressive loading and type III collagen content would be high in regions associated with damage. We further hypothesized that regions of high type III collagen content would correspond to regions with low tensile modulus and a low degree of collagen alignment. Although type III collagen content was not significantly higher in regions that are frequently damaged, all other hypotheses were supported by our results. In particular, type II collagen content was highest near the insertion while type III collagen was inversely correlated with tendon modulus and collagen alignment. The measured increase in type II collagen under the coracoacromial arch provides evidence of adaptation to compressive loading in the SST. Moreover, the structure-function relationship between type III collagen content and tendon mechanics established in this study demonstrates a mechanism for altered mechanical properties in pathological tendons and provides a guideline for identifying therapeutic targets and pathology-specific biomarkers.

The mechanical behavior of the annulus fibrosus (AF) of the intervertebral disc can be modeled as a mixture of fibers, extra-fibrillar matrix (EFM), ions, and fluid. However, the properties of the EFM have not been measured directly. We measured mechanical properties of the human EFM at several locations, determined the effect of age and degeneration, and evaluated whether changes in EFM properties correspond to AF compositional changes. EFM mechanical properties were measured using a method that combines osmotic loading and confined compression. AF samples were dissected from several locations, and mechanical properties were correlated with age, degeneration, and composition. EFM modulus was found to range between 10 and 50 kPa, increasing nonlinearly with compression magnitude and being highest in the AF outer-anterior region. EFM properties were not correlated with composition or degeneration. However, the EFM modulus, its relative contribution to tissue modulus, and model parameters were correlated with age. These measurements will result in more accurate predictions of deformations in the intervertebral disc. Additionally, parameters such as permeability and diffusivity used for biotransport analysis of glucose and other solutes depend on EFM deformation. Consequently, the accuracy of biotransport simulations will be greatly improved.

Background context Reduced spinal loading degrades the intervertebral disc and alters the muscle size. Purpose To determine the ability of high-frequency and low-intensity vibrations to maintain disc biomechanics and prevent muscle changes during hindlimb unloading. Study design Three groups of Sprague-Dawley rats were hindlimb unloaded for 4 weeks. In two hindlimb unloaded groups, unloading was interrupted for 15 min/d and the rats were positioned upright on a 90 Hz vertically oscillating plate or a sham control inactive plate. One author owns (provisional) patents regarding the application of vibrations to the musculoskeleton. Methods The motion segments L4–L5 were mechanically evaluated in compression-tension, axial creep, and torsion loading. In vivo microcomputed tomography was used to determine longitudinal psoas and paraspinal muscle area. This work was supported by National Institutes of Health, National Aeronautics and Space Administration (NASA), Alliance for Graduate Education and the Professoriate, and NASA-Harriett G. Jenkins Predoctoral and W. Burghardt Turner Fellowships. The author (SJ) holds (provisional) patents regarding the application of vibrations. Results There were no differences between the discs of uninterrupted unloading and sham animals and these groups were pooled. Compared with normally ambulating age-matched controls, hindlimb unloaded discs had altered properties in every loading modality. Psoas area of the unloaded rats increased at L4 and L5 and the paraspinal area decreased at L4. Vibrations (90 Hz) maintained compression-tension properties, partially maintained creep properties, but did not mitigate torsional weakening because of unloading. Low-intensity vibrations prevented the increase in psoas area but did not abate paraspinal muscle loss. Conclusions In support of clinical studies, unloading deconditioned the rodent disc and altered the muscle area. Although brief exposures to upright posture provided only limited benefits, low-intensity vibrations superimposed on upright posture served to preserve disc mechanics during unloading.

Low back pain is a significant socioeconomic burden and intervertebral disc degeneration has been implicated as a cause. A reliable animal model of disc degeneration is necessary to evaluate therapeutics, and functional metrics are essential to quantify their benefit. To this end, needle puncture injuries were created in the caudal intervertebral discs of mice to induce disc degeneration. Compression, torsion, and creep mechanics were assessed both immediately and after eight weeks to distinguish between the effects of injury and the subsequent reparative or degenerative response. Two needle sizes (29 and 26 gauge) were used to determine injury size-dependence. Compressive stiffness (62%), torsional stiffness (60%), and early damping stiffness (84%) decreased immediately after injury with the large needle (26G). These mechanical properties did not change over time despite structural and compositional changes. At 8 weeks following large needle injury, disc height decreased (37%), nucleus pulposus (NP) glycosaminoglycan content decreased (41%), and NP collagen content increased (45%). The small needle size had no significant effect on mechanics and did not initiate degenerative changes in structure and composition. Thus, the injection of therapeutics into the NP with a minimal needle size may limit damage due to the needle insertion. These findings, along with the wide commercial availability of mouse-specific biological probes, indicate that the mouse caudal disc model can be a powerful tool for investigating disc degeneration and therapy.

The annulus fibrosus (AF) of the intervertebral disk plays a critical role in vertebral load transmission that is heavily dependent on the microscale structure and composition of the tissue. With degeneration, both structure and composition are compromised, resulting in a loss of AF mechanical function. Numerous tissue engineering strategies have addressed the issue of AF degeneration, but few have focused on recapitulation of AF microstructure and function. One approach that allows for generation of engineered AF with appropriate (+/-)30° lamellar microstructure is the use of aligned electrospun scaffolds seeded with mesenchymal stem cells (MSCs) and assembled into angle-ply laminates (APL). Previous work indicates that opposing lamellar orientation is necessary for development of near native uniaxial tensile properties. However, most native AF tensile loads are applied biaxially, as the disk is subjected to multi-axial loads and is constrained by its attachments to the vertebral bodies. Thus, the objective of this study was to evaluate the biaxial mechanical response of engineered AF bilayers, and to determine the importance of opposing lamellar structure under this loading regime. Opposing bilayers, which replicate native AF structure, showed a significantly higher modulus in both testing directions compared to parallel bilayers, and reached ∼60% of native AF biaxial properties. Associated with this increase in biaxial properties, significantly less shear, and significantly higher stretch in the fiber direction, was observed. These results provide additional insight into native tissue structure-function relationships, as well as new benchmarks for engineering functional AF tissue constructs.

The annulus fibrosus (AF) of the intervertebral disc experiences cyclic tensile loading in vivo at various states of mechanical equilibrium. Disc degeneration is associated with alterations in the biochemical composition of the AF including decreased water content, decreased proteoglycan concentration, and increased collagen deposition that affect mechanical function of the AF in compression and shear. Such changes may also affect the dynamic viscoelastic properties of the AF and thus alter the disc's ability to dissipate energy under physiologic loading. The objectives of this study were to quantify the dynamic viscoelastic properties of human AF in circumferential tension and to determine the effect of degeneration on these properties. Nondegenerated and degenerated human AF tensile samples were tested in uniaxial tension over a spectrum of loading frequencies spanning 0.01Hz to 2Hz at several states of equilibrium strain to determine the dynamic viscoelastic properties (dynamic modulus, phase angle) using a linear viscoelastic model. The AF dynamic modulus increased at higher equilibrium strain levels. The AF behaved more elastically at higher frequencies with a decreased phase angle. Degeneration resulted in a higher dynamic modulus at all strain levels but had no effect on phase angle. The findings from this study elucidate the effect of degeneration on the dynamic viscoelastic properties of human AF and lend insight into the mechanical role of the AF in cyclic loading conditions.

The intervertebral disc maintains a balance between externally applied loads and internal osmotic pressure. Fluid flow plays a key role in this process, causing fluctuations in disc hydration and height. The objectives of this study were to quantify and model the axial creep and recovery responses of nondegenerate and degenerate human lumbar discs. Two experiments were performed. First, a slow compressive ramp was applied to 2000 N, unloaded to allow recovery for up to 24 hours, and re-applied. The linear-region stiffness and disc height were within 5% of the initial condition for recovery times greater than 8 hours. In the second experiment, a 1000 N creep load was applied for four hours, unloaded recovery monitored for 24 hours, and the creep load repeated. A viscoelastic model comprised of a “fast” and “slow” exponential response was used to describe the creep and recovery, where the fast response is associated with flow in the nucleus pulposus (NP) and endplate, while the slow response is associated with the annulus fibrosus (AF). The study demonstrated that recovery is 3-4X slower than loading. The fast response was correlated with degeneration, suggesting larger changes in the NP with degeneration compared to the AF. However, the fast response comprised only 10-15% of the total equilibrium displacement, with the AF-dominated slow response comprising 40-70%. Finally, the physiological loads and deformations and their associated long equilibrium times confirm that diurnal loading does not represent “equilibrium” in the disc, but that over time the disc is in steady-state.

Cross bridges are radial structures within the highly organized lamellar structure of the annulus fibrosus of the intervertebral disc that connect two or more non-consecutive lamellae. Their origin and function are unknown. During fetal development, blood vessels penetrate deep within the AF and recede during postnatal growth. We hypothesized that cross bridges are the pathways left by these receding blood vessels. Initially, the presence of cross bridges was confirmed in cadaveric human discs aged 25 and 53 years. Next, L1-L2 intervertebral discs (n=4) from sheep ranging in age from 75 days fetal gestation to adult were processed for paraffin histology. Mid-sagittal sections were immunostained for endothelial cell marker PECAM-1. The anterior and posterior AF were imaged using differential interference contrast microscopy, and the following parameters were quantified: total number of distinct lamellae, total number of cross bridges, percentage of cross bridges staining positive for PECAM-1, cross bridge penetration depth (% total lamellae), and PECAM-1 positive cross bridge penetration depth. Cross bridges were first observed at 100 days fetal gestation. The overall number peaked in neonates then remained relatively unchanged. The percentage of PECAM-1 positive cross bridges declined progressively from almost 100% at 100 days gestation to less than 10% in adults. Cross bridge penetration depth peaked in neonates then remained unchanged at subsequent ages. Depth of PECAM-1 positive cross bridges decreased progressively after birth. Findings were similar for both the anterior and posterior. The AF lamellar architecture is established early in development. It later becomes disrupted as a consequence of vascularization. Blood vessels then recede, perhaps due to increasing mechanical stresses in the surrounding matrix. In this study we present evidence that the pathways left by receding blood vessels remain as lamellar cross bridges. It is unclear whether the presence of cross bridges in the aging and degenerating intervertebral disc would be advantageous or detrimental, and this question should be addressed by future studies.

Intervertebral disc biodevices that employ motion-preservation strategies (e.g., nucleus replacement, total disc replacement and posterior stabilization devices) are currently in use or in development. However, their long-term performance is unknown and only a small number of randomized controlled trials have been conducted. In this article, we discuss the following biodevices: interbody cages, nuclear pulposus replacements, total disc replacements and posterior dynamic stabilization devices, as well as future biological treatments. These biodevices restore some function to the motion segment; however, contrary to expectations, the risk of adjacent-level degeneration does not appear to have been reduced. The short-term challenge is to replicate the complex biomechanical function of the motion segment (e.g., biphasic, viscoelastic behavior and nonlinearity) to improve the quality of motion and minimize adjacent level problems, while ensuring biodevice longevity for the younger, more active patient. Biological strategies for regeneration and repair of disc tissue are being developed and these offer exciting opportunities (and challenges) for the longer term. Responsible introduction and rigorous assessment of these new technologies are required. In this article, we will describe the properties of the disc, explore biodevices currently in use for the surgical treatment of low back pain (with an emphasis on lumbar total disc replacement) and discuss future directions for biological treatments. Finally, we will assess the challenges ahead for the next generation of biodevices designed to replace the disc.

Mechanical function of the annulus fibrosus of the intervertebral disc is dictated by the composition and microstructure of its highly ordered extracellular matrix. Recent work on engineered angle-ply laminates formed from mesenchymal stem cell (MSC)-seeded nanofibrous scaffolds indicates that the organization of collagen fibers into planes of alternating alignment may play an important role in annulus fibrosus tissue function. Specifically, these engineered tissues can resist tensile deformation through shearing of the interlamellar matrix as layers of collagen differentially reorient under load. In the present work, a hyperelastic constitutive model was developed to describe the role of interlamellar shearing in reinforcing the tensile response of biologic laminates, and was applied to experimental results from engineered annulus constructs formed from MSC-seeded nanofibrous scaffolds. By applying the constitutive model to uniaxial tensile stress–strain data for bilayers with three different fiber orientations, material parameters were generated that characterize the contributions of extrafibrillar matrix, fibers, and interlamellar shearing interactions. By 10 weeks of in vitro culture, interlamellar shearing accounted for nearly 50% of the total stress associated with uniaxial extension in the anatomic range of ply angle. The model successfully captured changes in function with extracellular matrix deposition through variations in the magnitude of model parameters with culture duration. This work illustrates the value of engineered tissues as tools to further our understanding of structure–function relations in native tissues and as a test-bed for the development of constitutive models to describe them.

Due to the inability of current clinical practices to restore function to degenerated intervertebral discs, the arena of disc tissue engineering has received substantial attention in recent years. Despite tremendous growth and progress in this field, translation to clinical implementation has been hindered by a lack of well-defined functional benchmarks. Because successful replacement of the disc is contingent upon replication of some or all of its complex mechanical behaviour, it is critically important that disc mechanics be well characterized in order to establish discrete functional goals for tissue engineering. In this review, the key functional signatures of the intervertebral disc are discussed and used to propose a series of native tissue benchmarks to guide the development of engineered replacement tissues. These benchmarks include measures of mechanical function under tensile, compressive and shear deformations for the disc and its substructures. In some cases, important functional measures are identified that have yet to be measured in the native tissue. Ultimately, native tissue benchmark values are compared to measurements that have been made on engineered disc tissues, identifying measures where functional equivalence was achieved, and others where there remain opportunities for advancement. Several excellent reviews exist regarding disc composition and structure, as well as recent tissue engineering strategies; therefore this review will remain focused on the functional aspects of disc tissue engineering.

BACKGROUND CONTEXT: Vertebral end plate support is necessary for successful lumbar total disc replacement (TDR) surgery. Failure to achieve anterior column support as a result of lumbar TDR device undersizing could lead to implant subsidence and fracture. PURPOSE: The purpose of the study was to examine the compressive biomechanical behavior of the vertebral end plate with varying sizes of disc replacement implants. STUDY DESIGN: The study design comprises a biomechanical investigation using a human cadaveric lumbar spine model. METHODS: Fifty-six vertebrae with intact posterior elements were prepared from 13 fresh frozen lumbar spines. Peripheral quantitative computed tomography was performed to assess regional bone density. Vertebrae were potted and subjected to nondestructive compression testing with a small, medium, and large custom-made implants with the footplate geometry of the ProDisc-LTDR (Synthes Spine, West Chester, PA, USA) system and having no keel. Failure testing was performed using the ProDisc-L implant with an intact keel. Pressure sensor film was used to assess contact pressure and distribution. RESULTS: Therewasalinearcorrelationbetweenpercentcoverageoftheendplateandimplant-end plate stiffness (p5.0001) and an inverse correlation with displacement (p5.01). The difference in implant-end plate stiffness between small-medium, medium-large, and small-large implants was 10.5% (p5.03), 10.2% (p5.02), and 19.6% (p!.0001), respectively. Failure analysis revealed similartrendsforimplantsizing,butonlybonedensitywasfoundtosignificantlycorrelatewithfailure properties (r50.76, p!.0001).Therewas a significant reduction in implant-end plate stiffness of 18% whenthekeelwasintactcomparedtowithoutthekeel(range6‐27%,p5.0008).Pressurefilmanalysis revealed that the implant was loaded peripherally and did not have central contact during nondestructiveloading.Therewasatrendtowardgreatercontactpressurewiththesmallimplantwhencompared with the medium implant (p5.06) and the large implant (p5.06). CONCLUSIONS: Although larger implants reduce end plate displacement, increase apparent implant-end plate stiffness, increase the implant-end plate contact area, and decrease the peak contact pressures, low bone density reduces failure properties. The keel introduces a reduction in stiffness to the implant-end plate interface, the clinical significance of which is currently

Modeling of connective tissues often includes collagen fibers explicitly as one of the components. These fibers can be oriented in many directions; therefore, several studies have considered statistical distributions to describe the fiber arrangement. One approach to formulate a constitutive framework for distributed fibers is to express the mechanical parameters, such as strain energy and stresses, in terms of angular integrals. These integrals represent the addition of the contribution of infinitesimal fractions of fibers oriented in a given direction. This approach leads to accurate results; however, it requires lengthy calculations. Recently, the use of generalized structure tensors has been proposed to represent the angular distribution in the constitutive equations of the fibers. Although this formulation is much simpler and fewer calculations are required, such structure tensors can only be used when all the fibers are in tension and the angular distribution is small. However, the amount of error introduced in these cases of non-tensile fiber loading and large angular distributions have not been quantified. Therefore, the objective of this study is to determine the range of values of angular distribution for which acceptable differences (less than 10%) between these two formulations are obtained. It was found, analytically and numerically, that both formulations are equivalent for planar distributions under equal-biaxial stretch. The comparison also showed, for other loading conditions, that the differences decrease when the fiber distribution is very small. Differences of less than 10% were usually obtained when the fiber distribution was very low (κ ≈ 0.03; κ ranges between 0 and 1/3, for aligned and isotropic distributed fibers, respectively). This range of angular distribution greatly limits the types of tissue that can be accurately analyzed using generalized structure tensors. It is expected that the results from this study guide the selection of a proper approach to analyze a particular tissue under a particular loading condition.

Successful engineering of load-bearing tissues requires recapitulation of their complex mechanical functions. Given the intimate relationship between function and form, biomimetic materials that replicate anatomic form are of great interest for tissue engineering applications. However, for complex tissues such as the annulus fibrosus, scaffolds have failed to capture their multi-scale structural hierarchy. Consequently, engineered tissues have yet to reach functional equivalence with their native counterparts. Here, we present a novel strategy for annulus fibrosus tissue engineering that replicates this hierarchy with anisotropic nanofibrous laminates seeded with mesenchymal stem cells. These scaffolds directed the deposition of an organized, collagen-rich extracellular matrix that mimicked the angle-ply, multi-lamellar architecture and achieved mechanical parity with native tissue after 10 weeks of in vitro culture. Furthermore, we identified a novel role for inter-lamellar shearing in reinforcing the tensile response of biologic laminates, a mechanism that has not previously been considered for these tissues.

Tendon exhibits nonlinear stress-strain behavior that may be partly due to movement of collagen fibers through the extracellular matrix. While a few techniques have been developed to evaluate the fiber architecture of other soft tissues, the organizational behavior of tendon under load has not been determined. The supraspinatus tendon (SST) of the rotator cuff is of particular interest for investigation due to its complex mechanical environment and corresponding inhomogeneity. In addition, SST injury occurs frequently with limited success in treatment strategies, illustrating the need for a better understanding of SST properties. Therefore, the objective of this study was to quantitatively evaluate the inhomogeneous tensile mechanical properties, fiber organization, and fiber realignment under load of human SST utilizing a novel polarized light technique. Fiber distributions were found to become more aligned under load, particularly during the low stiffness toe-region, suggesting that fiber realignment may be partly responsible for observed nonlinear behavior. Fiber alignment was found to correlate significantly with mechanical parameters, providing evidence for strong structure-function relationships in tendon. Human SST exhibits complex, inhomogeneous mechanical properties and fiber distributions, perhaps due to its complex loading environment. Surprisingly, histological grade of degeneration did not correlate with mechanical properties.

Background context Rat models with altered loading are used to study disc degeneration and mechano-transduction. Given the prominent role of mechanics in disc function and degeneration, it is critical to measure mechanical behavior to evaluate changes after model interventions. Axial compression mechanics of the rat disc are representative of the human disc when normalized by geometry, and differences between the lumbar and caudal disc have been quantified in axial compression. No study has quantified rat disc torsional mechanics. Purpose Compare the torsional mechanical behavior of rat lumbar and caudal discs, determine the contribution of combined axial load on torsional mechanics, and compare the torsional properties of rat discs to human lumbar discs. Study design Cadaveric biomechanical study. Methods Cyclic torsion without compressive load followed by cyclic torsion with a fixed compressive load was applied to rat lumbar and caudal disc levels. Results The apparent torsional modulus was higher in the lumbar region than in the caudal region: 0.081±0.026 (MPa/°, mean±SD) for lumbar axially loaded; 0.066±0.028 for caudal axially loaded; 0.091±0.033 for lumbar in pure torsion; and 0.056±0.035 for caudal in pure torsion. These values were similar to human disc properties reported in the literature ranging from 0.024 to 0.21MPa/°. Conclusions Use of the caudal disc as a model may be appropriate if the mechanical focus is within the linear region of the loading regime. These results provide support for use of this animal model in basic science studies with respect to torsional mechanics.

STUDY DESIGN Integrating theoretical and experimental approaches for annulus fibrosus (AF) functional tissue engineering. OBJECTIVE Apply a hyperelastic constitutive model to characterize the evolution of engineered AF via scalar model parameters. Validate the model and predict the response of engineered constructs to physiologic loading scenarios. SUMMARY OF BACKGROUND DATA There is need for a tissue engineered replacement for degenerate AF. When evaluating engineered replacements for load-bearing tissues, it is necessary to evaluate mechanical function with respect to the native tissue, including nonlinearity and anisotropy. METHODS Aligned nanofibrous poly-epsilon-caprolactone scaffolds with prescribed fiber angles were seeded with bovine AF cells and analyzed over 8 weeks, using experimental (mechanical testing, biochemistry, histology) and theoretical methods (a hyperelastic fiber-reinforced constitutive model). RESULTS The linear region modulus for phi = 0 degrees constructs increased by approximately 25 MPa, and for phi = 90 degrees by approximately 2 MPa from 1 day to 8 weeks in culture. Infiltration and proliferation of AF cells into the scaffold and abundant deposition of s-GAG and aligned collagen was observed. The constitutive model had excellent fits to experimental data to yield matrix and fiber parameters that increased with time in culture. Correlations were observed between biochemical measures and model parameters. The model was successfully validated and used to simulate time-varying responses of engineered AF under shear and biaxial loading. CONCLUSION AF cells seeded on nanofibrous scaffolds elaborated an organized, anisotropic AF-like extracellular matrix, resulting in improved mechanical properties. A hyperelastic fiber-reinforced constitutive model characterized the functional evolution of engineered AF constructs, and was used to simulate physiologically relevant loading configurations. Model predictions demonstrated that fibers resist shear even when the shearing direction does not coincide with the fiber direction. Further, the model suggested that the native AF fiber architecture is uniquely designed to support shear stresses encountered under multiple loading configurations.