Nathan T. Jacobs

Fiber angle and aspect ratio influence the shear mechanics of oriented electrospun nanofibrous scaffolds.

Fibrocartilages, including the knee meniscus and the annulus fibrosus (AF) of the intervertebral disc, play critical mechanical roles in load transmission across joints and their function is dependent upon well-defined structural hierarchies, organization, and composition. All, however, are compromised in the pathologic transformations associated with tissue degeneration. Tissue engineering strategies that address these key features, for example, aligned nanofibrous scaffolds seeded with mesenchymal stem cells (MSCs), represent a promising approach for the regeneration of these fibrous structures. While such engineered constructs can replicate native tissue structure and uniaxial tensile properties, the multidirectional loading encountered by these tissues in vivo necessitates that they function adequately in other loading modalities as well, including shear. As previous findings have shown that native tissue tensile and shear properties are dependent on fiber angle and sample aspect ratio, respectively, the objective of the present study was to evaluate the effects of a changing fiber angle and sample aspect ratio on the shear properties of aligned electrospun poly(ε-caprolactone) (PCL) scaffolds, and to determine how extracellular matrix deposition by resident MSCs modulates the measured shear response. Results show that fiber orientation and sample aspect ratio significantly influence the response of scaffolds in shear, and that measured shear strains can be predicted by finite element models. Furthermore, acellular PCL scaffolds possessed a relatively high shear modulus, 2-4 fold greater than native tissue, independent of fiber angle and aspect ratio. It was further noted that under testing conditions that engendered significant fiber stretch, the aggregate resistance to shear was higher, indicating a role for fiber stretch in the overall shear response. Finally, with time in culture, the shear modulus of MSC laden constructs increased, suggesting that deposited ECM contributes to the construct shear properties. Collectively, these findings show that aligned electrospun PCL scaffolds are a promising tool for engineering fibrocartilage tissues, and that the shear properties of both acellular and cell-seeded formulations can match or exceed native tissue benchmarks.

Elastic, permeability and swelling properties of human intervertebral disc tissues: A benchmark for tissue engineering.

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 discs mechanical behavior and transport.

Biaxial Tension of Fibrous Tissue: Using Finite Element Methods to Address Experimental Challenges Arising From Boundary Conditions and Anisotropy

Planar biaxial tension remains a critical loading modality for fibrous soft tissue and is widely used to characterize tissue mechanical response, evaluate treatments, develop constitutive formulas, and obtain material properties for use in finite element studies. Although the application of tension on all edges of the test specimen represents the in situ environment, there remains a need to address the interpretation of experimental results. Unlike uniaxial tension, in biaxial tension the applied forces at the loading clamps do not transmit fully to the region of interest (ROI), which may lead to improper material characterization if not accounted for. In this study, we reviewed the tensile biaxial literature over the last ten years, noting experimental and analysis challenges. In response to these challenges, we used finite element simulations to quantify load transmission from the clamps to the ROI in biaxial tension and to formulate a correction factor that can be used to determine ROI stresses. Additionally, the impact of sample geometry, material anisotropy, and tissue orientation on the correction factor were determined. Large stress concentrations were evident in both square and cruciform geometries and for all levels of anisotropy. In general, stress concentrations were greater for the square geometry than the cruciform geometry. For both square and cruciform geometries, materials with fibers aligned parallel to the loading axes reduced stress concentrations compared to the isotropic tissue, resulting in more of the applied load being transferred to the ROI. In contrast, fiber-reinforced specimens oriented such that the fibers aligned at an angle to the loading axes produced very large stress concentrations across the clamps and shielding in the ROI. A correction factor technique was introduced that can be used to calculate the stresses in the ROI from the measured experimental loads at the clamps. Application of a correction factor to experimental biaxial results may lead to more accurate representation of the mechanical response of fibrous soft tissue.

Human L3L4 intervertebral disc mean 3D shape, modes of variation, and their relationship to degeneration.

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.

Effect of orientation and targeted extracellular matrix degradation on the shear mechanical properties of the annulus fibrosus.

The intervertebral disc experiences combinations of compression, torsion, and bending that subject the disc substructures, particularly the annulus fibrosus (AF), to multidirectional loads and deformations. Combined tensile and shear loading is a particularly important loading paradigm, as compressive loads place the AF in circumferential hoop tension, and spine torsion or bending induces AF shear. Yet the anisotropy of AF mechanical properties in shear, as well as important structure-function mechanisms governing this response, are not well-understood. The objective of this study, therefore, was to investigate the effects of tissue orientation and enzymatic degradation of glycosaminoglycan (GAG) and elastin on AF shear mechanical properties. Significant anisotropy was found: the circumferential shear modulus, Gθz, was an order of magnitude greater than the radial shear modulus, Grθ. In the circumferential direction, prestrain significantly increased the shear modulus, suggesting an important role for collagen fiber stretch in shear properties for this orientation. While not significant and highly variable, ChABC treatment to remove GAG increased the circumferential shear modulus compared to PBS control (p=0.15). Together with the established literature for tensile loading of fiber-reinforced GAG-rich tissues, the trends for changes in shear modulus with ChABC treatment reflect complex, structure-function relationships between GAG and collagen that potentially occur over several hierarchical scales. Elastase digestion did not significantly affect shear modulus with respect to PBS control; further contributing to the notion that circumferential shear modulus is dominated by collagen fiber stretch. The results of this study highlight the complexity of the structure-function relationships that govern the mechanical response of the AF in radial and circumferential shear, and provide new and more accurate data for the validation of material models and tissue-engineered disc replacements.

Validation and application of an intervertebral disc finite element model utilizing independently constructed tissue-level constitutive formulations that are nonlinear, anisotropic, and time-dependent.

Finite element (FE) models are advantageous in the study of intervertebral disc mechanics as the stress-strain distributions can be determined throughout the tissue and the applied loading and material properties can be controlled and modified. However, the complicated nature of the disc presents a challenge in developing an accurate and predictive disc model, which has led to limitations in FE geometry, material constitutive models and properties, and model validation. The objective of this study was to develop a new FE model of the intervertebral disc, to validate the models nonlinear and time-dependent responses without tuning or calibration, and to evaluate the effect of changes in nucleus pulposus (NP), cartilaginous endplate (CEP), and annulus fibrosus (AF) material properties on the disc mechanical response. The new FE disc model utilized an analytically-based geometry. The model was created from the mean shape of human L4/L5 discs, measured from high-resolution 3D MR images and averaged using signed distance functions. Structural hyperelastic constitutive models were used in conjunction with biphasic-swelling theory to obtain material properties from recent tissue tests in confined compression and uniaxial tension. The FE disc model predictions fit within the experimental range (mean ± 95% confidence interval) of the discs nonlinear response for compressive slow loading ramp, creep, and stress-relaxation simulations. Changes in NP and CEP properties affected the neutral-zone displacement but had little effect on the final stiffness during slow-ramp compression loading. These results highlight the need to validate FE models using the discs full nonlinear response in multiple loading scenarios.

Toward tactilely transparent gloves: Collocated slip sensing and vibrotactile actuation

Tactile information plays a critical role in the human ability to manipulate objects with ones hands. Many environments require the use of protective gloves that diminish essential tactile feedback. Under these circumstances, seemingly simple tasks such as picking up an object can become very difficult. This paper introduces the SlipGlove, a novel device that uses an advanced sensing and actuation system to return this vital tactile information to the user. Our SlipGlove prototypes focus on providing tactile cues associated with slip between the glove and a contact surface. Relative motion is sensed using optical mouse sensors embedded in the gloves surface. This information is conveyed to the wearer via miniature vibration motors placed inside the glove against the wearers skin. The collocation of slip sensing and tactile feedback creates a system that is natural and intuitive to use. We report results from a human subject study demonstrating that the SlipGlove allows the wearer to approach the capabilities of bare skin in detecting and reacting to fingertip slip. Users of the SlipGlove also had significantly faster and more consistent reaction to fingertip slip when compared to a traditional glove design. The SlipGlove technology allows us to enhance human perception when interacting with real environments and move toward the goal of a tactilely transparent glove.

Novel human intervertebral disc strain template to quantify regional three-dimensional strains in a population and compare to internal strains predicted by a finite element model

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.

A model for intervertebral disc shape variation in a population including the effect of degeneration

The 3D shape of the intervertebral disc (IVD) is modeled in this work from a population of samples. Heretofore, 3D models of the disc have been created from individual samples only. L3/L4 discs were segmented from 14 magnetic resonance images. The mean shape was calculated using a signed distance function representation, and principal components analysis was used to re-express the residual variability with an orthonormal basis. The first basis vector, which represents lateral bulging of the disc, expresses the majority of the populations shape variation. This vector was also significantly correlated with Pfirrmann grade (p<0.01), a measure of disc degeneration. The IVD shape model can thus be used to generate an average IVD shape for the population, as well as shapes for specific degenerative grades.

Material Properties Over-Estimated by Boundary Conditions in Biaxial Tension Can Be Corrected Using Finite Element Analysis

Tissue modeling requires an appropriate stress-strain constitutive relationship and a corresponding set of material properties. It is often the goal of experimental studies to determine these material properties. Uniaxial tension experiments are simple in experimental approach and the interpretation of results is straightforward, prompting its use in several studies. However, the freely contracting lateral edges observed in this loading modality do not mimic the in situ conditions of many fiber-reinforced soft tissues and the unconstrained fibers may also create errors.Copyright