Extrafibrillar Matrix Research Articles

Compressive behaviour of uniaxially aligned individual mineralised collagen fibres at the micro- and nanoscale

The increasing incidence of osteoporotic bone fractures makes fracture risk prediction an important clinical challenge. Computational models can be utilised to facilitate such analyses. However, they critically depend on bone’s underlying hierarchical material description. To understand bone’s irreversible behaviour at the micro- and nanoscale, we developed an in situ testing protocol that allows us to directly relate the experimental data to the mechanical behaviour of individual mineralised collagen fibres and its main constitutive phases, the mineralised collagen fibrils and the mineral nanocrystals, by combining micropillar compression of single fibres with small angle X-ray scattering (SAXS) and X-ray diffraction (XRD). Failure modes were assessed by SEM. Strain ratios in the elastic region at fibre, fibril and mineral levels were found to be approximately 22:5:2 with strain ratios at the point of compressive strength of 0.23 ± 0.11 for fibril-to-fibre and 0.07 ± 0.01 for mineral-to-fibre levels. Mineral-to-fibre levels showed highest strain ratios around the apparent yield point, fibril-to-fibre around apparent strength. The mineralised collagen fibrils showed a delayed mechanical response, contrary to the mineral phase, which points towards preceding deformations of mineral nanocrystals in the extrafibrillar matrix. No damage was measured at the level of the mineralised collagen fibre which indicates an incomplete separation of the mineral and collagen, and an extrafibrillar interface failure. The formation of kink bands and the gradual recruitment of fibrils upon compressive loading presumably led to localised strains. Our results from a well-controlled fibrillar architecture provide valuable input for micromechanical models and computational non-linear bone strength analyses that may provide further insights for personalised diagnosis and treatment as well as bio-inspired implants for patients with bone diseases. Statement of significanceMusculoskeletal diseases such as osteoporosis, osteoarthritis or bone cancer significantly challenge health care systems and make fracture risk prediction and treatment optimisation important clinical goals. Computational methods such as finite element models have the potential to optimise analyses but highly depend on underlying material descriptions. We developed an in situ testing set-up to directly relate experimental data to the mechanical behaviour of bone’s fundamental building block, the individual mineralised collagen fibre and its main constituents. Low multilevel strain ratios suggest high deformations in the extrafibrillar matrix and energy dissipation at the interfaces, the absence of damage indicates both an incomplete separation between mineral and collagen and an extrafibrillar interface failure. The formation of kink bands in the fibril-reinforced composite presumably led to localised strains. The deformation behaviour of a well-controlled fibrillar architecture provides valuable input for non-linear bone strength analyses.

Bioinspired design of hybrid composite materials

ABSTRACTMimicking the natural design motifs of structural biological materials is a promising approach to achieve a unique combination of strength and toughness for engineering materials. In this study, we proposed a 2D computational model, which is a two-hierarchy hybrid composite inspired by the ultrastructural features of bone. The model is composed of alternating parallel array of two subunits (A & B) mimicking ‘mineralized collagen fibril’ and ‘extrafibrillar matrix’ of bone at ultrastructural level. The subunit-A is formed by short stiff platelets embedded within a soft matrix. The subunit-B consists of randomly distributed stiff grains bonded by a thin layer of tough adhesive phase. To assess the performance of the bioinspired design, a conventional unidirectional long-fiber composite made with the same amount of hard and soft phases was studied. The finite element simulation results indicated that the toughness, strength and elastic modulus of the bioinspired composite was 312%, 83%, and 55% of that of the conventional composite, respectively. The toughness improvement was attributed to the prevalent energy-dissipating damage of adhesive phase in subunit-B and crack-bridging by subunit-A, the two major toughening mechanisms in the model. This study exemplifies some insights into natural design of materials to gain better material performance.

Computational investigation of ultrastructural behavior of bone using a cohesive finite element approach.

Bone ultrastructure at sub-lamellar length scale is a key structural unit in bone that bridges nano- and microscale hierarchies of the tissue. Despite its influence on bulk response of bone, the mechanical behavior of bone at ultrastructural level remains poorly understood. To fill this gap, in this study, a two-dimensional cohesive finite element model of bone at sub-lamellar level was proposed and analyzed under tensile and compressive loading conditions. In the model, ultrastructural bone was considered as a composite of mineralized collagen fibrils (MCFs) embedded in an extrafibrillar matrix (EFM) that is comprised of hydroxyapatite (HA) polycrystals bounded via thin organic interfaces of non-collagenous proteins (NCPs). The simulation results indicated that in compression, EFM dictated the pre-yield deformation of the model, then damage was initiated via relative sliding of HA polycrystals along the organic interfaces, and finally shear bands were formed followed by delamination between MCF and EFM and local buckling of MCF. In tension, EFM carried the most of load in pre-yield deformation, and then an array of opening-mode nano-cracks began to form within EFM after yielding, thus gradually transferring the load to MCF until failure, which acted as crack bridging filament. The failure modes, stress-strain curves, and in situ mineral strain of ultrastructural bone predicted by the model were in good agreement with the experimental observations reported in the literature, thus suggesting that this model can provide new insights into sub-microscale mechanical behavior of bone.

Second-harmonic patterned polarization-analyzed reflection confocal microscopy of stromal collagen in benign and malignant breast tissues

We present the results of polarimetric analysis of collagen on varying pathologies of breast tissues using second-harmonic patterned polarization-analyzed reflection confocal (SPPARC) microscopy. Experiments are conducted on a breast tissue microarray having benign tissues (BT), malignant invasive lobular carcinoma (ILC), and benign stroma adjacent to the malignant tissues (called the benign adjacent tissue, or BAT). Stroma in BAT and ILC exhibit the largest parameter differences. We observe that stromal collagen readings in ILC show lower depolarization, lower diattenuation and higher linear degree-of-polarization values than stromal collagen in BAT. This suggests that the optical properties of collagen change most in the vicinity of tumors. A similar trend is also exhibited in the non-collagenous extrafibrillar matrix plus cells (EFMC) region. The three highlighted parameters show greatest sensitivity to changes in the polarization response of collagen between pathologies.

Swelling of fiber-reinforced soft tissues is affected by fiber orientation, fiber stiffness, and lamella structure.

Native and engineered fiber-reinforced tissues are composites comprised of stiff collagen fibers embedded within an extrafibrillar matrix that is capable of swelling by absorbing water molecules. Tissue swelling is important for understanding stress distributions between collagen fibers and extrafibrillar matrix, as well as for understanding mechanisms of tissue failure. The swelling behavior of fiber-reinforced tissues in the musculoskeletal system has been largely attributed to the glycosaminoglycan content. Recent work demonstrated anisotropy in the swelling response of the annulus fibrosus in the intervertebral disc. It is well known that collagen fiber orientation affects elastic behavior, but the effect of collagen fiber network on tissue swelling behavior is not well understood. In this study, we developed three series of models to evaluate the effect of collagen fiber orientation, fiber network architecture (i.e., single or multi-fiber families within a layer), and fiber stiffness on bulk tissue swelling, which was simulated by describing the extrafibrillar matrix as a triphasic material, as proposed by Lai et al. Model results were within one standard deviation of reported mean values for changes in tissue volume, width, and thickness under free swelling conditions. The predicted swelling response of single-fiber family structures was highly dependent on fiber orientation and the number of lamellae in the bulk tissue. Moreover, matrix swelling resulted in tissue to twist, which reduced fiber deformations, demonstrating a balance between fiber deformation and matrix swelling. Large changes in fiber stiffness (20 × increase) had a relatively small effect on tissue swelling (~ 2% decrease in swelling). In conclusion, fiber angle, fiber architecture (defined as single- versus multiple fiber families in a layer), and the number of layers in a single fiber family structure directly affected tissue swelling behavior, including fiber stretch, fiber reorientation, and tissue deformation. These findings support the need to develop computational models that closely mimic the native architecture in order to understand mechanisms of stress distributions and tissue failure.

Effect of modifications in mineralized collagen fibril and extra-fibrillar matrix material properties on submicroscale mechanical behavior of cortical bone.

A key length scale of interest in assessing the fracture resistance of bone is the submicroscale which is composed of mineralized collagen fibrils (MCF) and extra-fibrillar matrix (EFM). Although the processes through which the submicroscale constituents of bone contribute to the fracture resistance in bone have been identified, the extent of the modifications in submicroscale mechanical response due to the changes in individual properties of MCFs and EFM has not been determined. As a result, this study aims to quantify the influence of individual MCF and EFM material property modifications on the mechanical behavior (elastic modulus, ultimate strength, and resistance to failure) of bone at the submicroscale using a novel finite element modeling approach that incorporate 3D networks of MCFs with three different orientations as well as explicit representation of EFM. The models were evaluated under tensile loading in transverse (representing MCF separation) and longitudinal (representing MCF rupture) directions. The results showed that the apparent elastic modulus at the submicroscale under both loading directions for all orientations was only affected by the change in the elastic modulus of MCFs. MCF separation and rupture strengths were mainly dependent on the ultimate strength of EFM and MCFs, respectively, with minimal influence of other material properties. The extent of damage during MCF separation increased with increasing ultimate strength of EFM and decreased with increasing fracture energy of EFM with minimal contribution from elastic modulus of MCFs. For MCF rupture, there was an almost one-to-one linear relationship between the percent change in fracture energy of MCFs and the percent change in the apparent submicroscale fracture energy. The ultimate strength and elastic modulus of MCFs had moderate to limited influence on the MCF rupture fracture energy. The results of this study quantified the extent of changes that may be seen in the energy dissipation processes during MCF rupture and separation relative to the changes in the individual constituents of the tissue. This new knowledge significantly contributes to improving the understanding of how the material property alterations at the submicroscale that can occur due to diseases, age-related changes, and treatments affect the fracture processes at larger length scales.

Mineralized collagen fibril network spatial arrangement influences cortical bone fracture behavior

Bone is a hierarchical material exhibiting different fracture mechanisms at each length scale. At the submicroscale, the bone is composed of mineralized collagen fibrils (MCF). At this scale, the fracture processes in cortical bone have not been extensively studied in the literature. In this study, the influence of MCF size and orientation on the fracture behavior of bone under both transverse and longitudinal loading was investigated using novel 3D models of MCF networks with explicit representation of extra-fibrillar matrix. The simulation results showed that separation between MCFs was the main cause of damage and failure under transverse loading whereas under longitudinal loading, the main damage and failure mechanism was MCF rupture. When the MCF network was loaded in the transverse direction the mechanical properties increased as the orientation of fibrils deviated farther from the main fibril orientation whereas the opposite trend was observed under longitudinal loading. The fracture energy was much larger in longitudinal than transverse loading. MCF diameter variation did not affect the mechanical properties under longitudinal loading but led to higher mechanical properties with increasing MCF diameter under transverse loading. The new modeling framework established in this study generate unique information on the effect of MCF network spatial arrangement on the fracture behavior of bone at the submicroscale which is not currently possible to measure via experiments. This unique information may improve the understanding of how structural alterations at the submicroscale due to disease, age-related changes, and treatments affect the fracture processes at larger length scales.

Second-harmonic patterned polarization-analyzed reflection confocal microscope.

We introduce the second-harmonic patterned polarization-analyzed reflection confocal (SPPARC) microscope-a multimodal imaging platform that integrates Mueller matrix polarimetry with reflection confocal and second-harmonic generation (SHG) microscopy. SPPARC microscopy provides label-free three-dimensional (3-D), SHG-patterned confocal images that lend themselves to spatially dependent, linear polarimetric analysis for extraction of rich polarization information based on the Mueller calculus. To demonstrate its capabilities, we use SPPARC microscopy to analyze both porcine tendon and ligament samples and find differences in both circular degree-of-polarization and depolarization parameters. Moreover, using the collagen-generated SHG signal as an endogenous counterstain, we show that the technique can be used to provide 3-D polarimetric information of the surrounding extrafibrillar matrix plus cells or EFMC region. The unique characteristics of SPPARC microscopy holds strong potential for it to more accurately and quantitatively describe microstructural changes in collagen-rich samples in three spatial dimensions.

Annulus fibrosus functional extrafibrillar and fibrous mechanical behaviour: experimental and computational characterisation.

The development of current surgical treatments for intervertebral disc damage could benefit from virtual environment accounting for population variations. For such models to be reliable, a relevant description of the mechanical properties of the different tissues and their role in the functional mechanics of the disc is of major importance. The aims of this work were first to assess the physiological hoop strain in the annulus fibrosus in fresh conditions (n = 5) in order to extract a functional behaviour of the extrafibrillar matrix; then to reverse-engineer the annulus fibrosus fibrillar behaviour (n = 6). This was achieved by performing both direct and global controlled calibration of material parameters, accounting for the whole process of experimental design and in silico model methodology. Direct-controlled models are specimen-specific models representing controlled experimental conditions that can be replicated and directly comparing measurements. Validation was performed on another six specimens and a sensitivity study was performed. Hoop strains were measured as 17 ± 3% after 10 min relaxation and 21 ± 4% after 20–25 min relaxation, with no significant difference between the two measurements. The extrafibrillar matrix functional moduli were measured as 1.5 ± 0.7 MPa. Fibre-related material parameters showed large variability, with a variance above 0.28. Direct-controlled calibration and validation provides confidence that the model development methodology can capture the measurable variation within the population of tested specimens.

Nanoscale deformation mechanisms and yield properties of hydrated bone extracellular matrix.

Bone's hierarchical architecture combines toughness and strength. To understand the governing deformation mechanisms, its postyield behaviour was analyzed at the microscale. Micropillars were compressed in physiological solution; an increased strength compared to macroscale and an influence of hydration was found. Transmission electron microscopy revealed cracks forming in regions with adverse fibril orientation in transverse pillars. In axial pillars kink bands were observed and shear cracks emerged at the interface of ordered and disordered regions. It was estimated that bone's fracture toughness for splitting between fibrils is significantly smaller than on the macroscale. This study underlines the importance of bone's hierarchical microstructure and the need to study structure-property relationships on all length scales.

An improved interfacial bonding model for material interface modeling

An improved interfacial bonding model was proposed from potential function point of view to investigate interfacial interactions in polycrystalline materials. It characterizes both attractive and repulsive interfacial interactions and can be applied to model different material interfaces. The path dependence of work-of-separation study indicates that the transformation of separation work is smooth in normal and tangential direction and the proposed model guarantees the consistency of the cohesive constitutive model. The improved interfacial bonding model was verified through a simple compression test in a standard hexagonal structure. The error between analytical solutions and numerical results from the proposed model is reasonable in linear elastic region. Ultimately, we investigated the mechanical behavior of extrafibrillar matrix in bone and the simulation results agreed well with experimental observations of bone fracture.

Contribution of extrafibrillar matrix to the mechanical behavior of bone using a novel cohesive finite element model

The mechanical behavior of bone is determined at all hierarchical levels, including lamellae (the basic building block of bone) that are comprised of mineralized collagen fibrils and extrafibrillar matrix. The mechanical behavior of mineralized collagen fibrils has been investigated intensively using both experimental and computational approaches. Yet, the contribution of the extrafibrillar matrix to bone mechanical properties is poorly documented. In this study, we intended to address this issue using a novel cohesive finite element (FE) model, in conjunction with the experimental observations reported in the literature. In the FE model, the extrafibrillar matrix was considered as a nanocomposite of hydroxyapatite (HA) crystals bounded through a thin organic interface modeled as a cohesive interfacial zone. The parameters required by the cohesive FE model were defined based on the experimental data reported in the literature. This hybrid nanocomposite model was tested in two loading modes (i.e. tension and compression) and under two hydration conditions (i.e. wet and dry). The simulation results indicated that (1) the failure modes of the extrafibrillar matrix predicted using the cohesive FE model were closely coincided with those experimentally observed in tension and compression tests; (2) the pre-yield deformation (i.e. internal strain) of HA crystals with respect to the applied strain was consistent with that obtained from the synchrotron X-ray scattering measurements irrespective of the loading modes and hydration status; and (3) the mechanical behavior of the extrafibrillar matrix was dictated by the properties of the organic interface between the HA crystals. Taken together, we postulate that the extrafibrillar matrix plays a major role in the pre-yield deformation and the failure mode of bone, thus, giving rise to important insights in the ultrastructural origins of bone fragility.

Microscopic and ultrastructural modifications of postmenopausal atrophic vaginal mucosa after fractional carbon dioxide laser treatment.

Vaginal atrophy occurring during menopause is closely related to the dramatic decrease in ovarian estrogens due to the loss of follicular activity. Particularly, significant changes occur in the structure of the vaginal mucosa, with consequent impairment of many physiological functions. In this study, carried out on bioptic vaginal mucosa samples from postmenopausal, nonestrogenized women, we present microscopic and ultrastructural modifications of vaginal mucosa following fractional carbon dioxide (CO2) laser treatment. We observed the restoration of the vaginal thick squamous stratified epithelium with a significant storage of glycogen in the epithelial cells and a high degree of glycogen-rich shedding cells at the epithelial surface. Moreover, in the connective tissue constituting the lamina propria, active fibroblasts synthesized new components of the extracellular matrix including collagen and ground substance (extrafibrillar matrix) molecules. Differently from atrophic mucosa, newly-formed papillae of connective tissue indented in the epithelium and typical blood capillaries penetrating inside the papillae, were also observed. Our morphological findings support the effectiveness of fractional CO2 laser application for the restoration of vaginal mucosa structure and related physiological trophism. These findings clearly coupled with striking clinical relief from symptoms suffered by the patients before treatment.

Molecular Dynamics Investigation on Shearing between Osteopontin and Hydroxyapatite in Biological Materials

Bone, a hard biological material, possesses a combination of high stiffness and toughness, even though the main basic building blocks of bone are simply mineral platelets and protein molecules. Bone has a very complex microstructure with at least seven hierachical levels. This unique material characteristic attracts great attention, but the deformation mechanisms in bone have not been well understood. Simulation at nanolength scale such as molecular dynamics (MD) is proven to be a powerful tool to investigate bone nanomechanics for developing new artificial biological materials. This study focuses on the ultra large and thin layer of extrafibrillar protein matrix (thickness = ~ 1 nm) located between mineralized collagen fibrils (MCF). Non-collagenous proteins such as osteopontin (OPN) can be found in this protein matrix, while MCF consists mainly of hydroxyapatite (HA) nanoplatelets (thickness = 1.5 4.5 nm). By using molecular dynamics method, an OPN peptide was pulled between two HA mineral platelets with water in presence. Periodic boundary condition (PBC) was applied. The results indicate that the mechanical response of OPN peptide greatly depends on the attractive electrostatics interaction between the acidic residues in OPN peptide and HA mineral surfaces. These bonds restrict the movement of OPN peptide, leading to a high energy dissipation under shear loading.

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−16m4/Ns) higher than the NP and CEP (~5.5 E−16m4/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.

Mechanical properties of the extra‐fibrillar matrix of human annulus fibrosus are location and age dependent

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.

Extra-fibrillar matrix mechanics of annulus fibrosus in tension and compression

The annulus fibrosus (AF) of the disk is a highly nonlinear and anisotropic material that undergoes a complex combination of loads in multiple orientations. The tensile mechanical behavior of AF in the lamellar plane is dominated by collagen fibers and has been accurately modeled using exponential functions. On the other hand, AF mechanics perpendicular to the lamella, in the radial direction, depend on the properties of the ground matrix with little to no fiber contribution. The ground matrix is mainly composed of proteoglycans (PG), which are negatively charged macromolecules that maintain the tissue hydration via osmotic pressure. The mechanical response of the ground matrix can be divided in the contribution of osmotic pressure and an elastic solid part known as extra-fibrillar matrix (EFM). Mechanical properties of the ground matrix have been measured using tensile and confined compression tests. However, EFM mechanics have not been measured directly. The objective of this study was to measure AF nonlinear mechanics of the EFM in tension and compression. To accomplish this, a combination of osmotic swelling and confined compression in disk radial direction, perpendicular to the lamella, was used. For this type of analysis, it was necessary to define a stress-free reference configuration. Thus, a brief analysis on residual stress in the disk and a procedure to estimate the reference configuration are presented. The proposed method was able to predict similar swelling deformations when using different loading protocols and models for the EFM, demonstrating its robustness. The stress-stretch curve of the EFM was linear in the range 0.9 < λ₃ < 1.3 with an aggregate modulus of 10.18±3.32 kPa; however, a significant nonlinearity was observed for compression below 0.8. The contribution of the EFM to the total aggregate modulus of the AF decreased from 70 to 30% for an applied compression of 50% of the initial thickness. The properties obtained in this study are essential for constitutive and finite element models of the AF and disk and can be applied to differentiate between functional degeneration effects such as PG loss and stiffening due to cross-linking.

Human annulus fibrosus material properties from biaxial testing and constitutive modeling are altered with degeneration

The annulus fibrosus (AF) of the intervertebral disk undergoes large and multidirectional stresses and strains. Uniaxial tensile tests are limited for measuring AF material properties, because freely contracting edges can prevent fiber stretch and are not representative of in situ boundary conditions. The objectives of this study were to measure human AF biaxial tensile mechanics and to apply and validate a constitutive model to determine material properties. Biaxial tensile tests were performed on samples oriented along the circumferential-axial and the radial-axial directions. Data were fit to a structurally motivated anisotropic hyperelastic model composed of isotropic extra-fibrillar matrix, nonlinear fibers, and fiber-matrix interactions (FMI) normal to the fibers. The validated model was used to simulate shear and uniaxial tensile behavior, to investigate AF structure-function, and to quantify the effect of degeneration. The biaxial stress-strain response was described well by the model (R (2)>0.9). The model showed that the parameters for fiber nonlinearity and the normal FMI correlated with degeneration, resulting in an elongated toe-region and lower stiffness with degeneration. The model simulations in shear and uniaxial tension successfully matched previously published circumferential direction Young's modulus, provided an explanation for the low values in previously published axial direction Young's modulus, and was able to simulate shear mechanics. The normal FMI were important contributors to stress and changed with degeneration, therefore, their microstructural and compositional source should be investigated. Finally, the biaxial mechanical data and constitutive model can be incorporated into a disk finite element model to provide improved quantification of disk mechanics.

Modeling interlamellar interactions in angle-ply biologic laminates for annulus fibrosus tissue engineering

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 10weeks 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.

Homologous structure–function relationships between native fibrocartilage and tissue engineered from MSC-seeded nanofibrous scaffolds

Understanding the interplay of composition, organization and mechanical function in load-bearing tissues is a prerequisite in the successful engineering of tissues to replace diseased ones. Mesenchymal stem cells (MSCs) seeded on electrospun scaffolds have been successfully used to generate organized tissues that mimic fibrocartilages such as the knee meniscus and the annulus fibrosus of the intervertebral disc. While matrix deposition has been observed in parallel with improved mechanical properties, how composition, organization, and mechanical function are related is not known. Moreover, how this relationship compares to that of native fibrocartilage is unclear. Therefore, in the present work, functional fibrocartilage constructs were formed from MSC-seeded nanofibrous scaffolds, and the roles of collagen and glycosaminoglycan (GAG) in compressive and tensile properties were determined. MSCs deposited abundant collagen and GAG over 120 days of culture, and these extracellular molecules were organized in such a way that they performed similar mechanical functions to their native roles: collagen dominated the tensile response while GAG was important for compressive properties. GAG removal resulted in significant stiffening in tension. A similar stiffening response was observed when GAG was removed from native inner annulus fibrosus, suggesting an interaction between collagen fibers and their surrounding extrafibrillar matrix that is shared by both engineered and native fibrocartilages. These findings strongly support the use of electrospun scaffolds and MSCs for fibrocartilage tissue engineering, and provide insight on the structure–function relations of both engineered and native biomaterials.