Elastic deformation of mineralized collagen fibrils: an equivalent inclusion based composite model

Mineralized collagen fibrils are the basic building blocks of bone tissue at the supramolecular level. Several disease states, manipulation of the expression of specific proteins involved in biomineralization, and treatment with different agents alter the extent of mineralization as well as the morphology of mineral crystals which in turn affect the mechanical function of bone tissue. An experimental assessment of mineralized fibers' mechanical properties is challenged by their small size, leaving analytical and computational models as a viable alternative for investigation of the fibril-level mechanical properties. In the current study the variation of the elastic stiffness tensor of mineralized collagen fibrils with changing mineral volume fraction and mineral aspect ratios was predicted via a micromechanical model. The partitioning of applied stresses between mineral and collagen phases is also predicted for normal and shear loading of fibrils. Model predictions resulted in transversely isotropic collagen fibrils in which the modulus along the longer axis of the fibril was the greatest. All the elastic moduli increased with increasing mineral volume fraction whereas Poisson's ratios decreased with the exception of v12 (=v21). The partitioning of applied stresses were such that the stresses acting on mineral crystals were about 1.5, 15, and 3 times greater than collagen stresses when fibrils were loaded transversely, longitudinally, and in shear, respectively. In the overall the predictions were such that: (a) greatest modulus along longer axis; (b) the greatest mineral/collagen stress ratio along the longer axis of collagen fibers (i.e., greatest relief of stresses acting on collagen); and (c) minimal lateral contraction when fibers are loaded along the longer axis. Overall, the pattern of mineralization as put forth in this model predicts a superior mechanical function along the longer axis of collagen fibers, the direction which is more likely to experience greater stresses.     

Tensile behavior of cortical bone: dependence of organic matrix material properties on bone mineral content

A porous composite model is developed to analyze the tensile mechanical properties of cortical bone. The effects of microporosity (volksman's canals, osteocyte lacunae) on the mechanical properties of bone tissue are taken into account. A simple shear lag theory, wherein tensile loads are transferred between overlapped mineral platelets by shearing of the organic matrix, is used to model the reinforcement provided by mineral platelets. It is assumed that the organic matrix is elastic in tension and elastic-perfectly plastic in shear until it fails. When organic matrix shear stresses at the ends of mineral platelets reach their yield values, the stress-strain curve of bone tissue starts to deviate from linear behavior. This is referred as the microscopic yield point. At the point where the stress-strain behavior of bone shows a sharp curvature, the organic phase reaches its shear yield stress value over the entire platelet. This is referred as the macroscopic yield point. It is assumed that after macroscopic yield, mineral platelets cannot contribute to the load bearing capacity of bone and that the mechanical behavior of cortical bone tissue is determined by the organic phase only. Bone fails when the principal stress of the organic matrix is reached. By assuming that mechanical properties of the organic matrix are dependent on bone mineral content below the macroscopic yield point, the model is used to predict the entire tensile mechanical behavior of cortical bone for different mineral contents. It is found that decreased shear yield stresses and organic matrix elastic moduli are required to explain the mechanical behavior of bones with lowered mineral contents. Under these conditions, the predicted values (elastic modulus, 0.002 yield stress and strain, and ultimate stress and strain) are within 15% of experimental data.     

Elastic properties of woven bone: effect of mineral content and collagen fibrils orientation

Woven bone is a type of tissue that forms mainly during fracture healing or fetal bone development. Its microstructure can be modeled as a composite with a matrix of mineral (hydroxyapatite) and inclusions of collagen fibrils with a more or less random orientation. In the present study, its elastic properties were estimated as a function of composition (degree of mineralization) and fibril orientation. A self-consistent homogenization scheme considering randomness of inclusions’ orientation was used for this purpose. Lacuno-canalicular porosity in the form of periodically distributed void inclusions was also considered. Assuming collagen fibrils to be uniformly oriented in all directions led to an isotropic tissue with a Young’s modulus \(E = 1.90\) GPa, which is of the same order of magnitude as that of woven bone in fracture calluses. By contrast, assuming fibrils to have a preferential orientation resulted in a Young’s modulus in the preferential direction of 9–16 GPa depending on the mineral content of the tissue. These results are consistent with experimental evidence for woven bone in foetuses, where collagen fibrils are aligned to a certain extent.     

Fracture Mechanics of Collagen Fibrils: Influence of Natural Cross-Links

Tendons are important load-bearing structures, which are frequently injured in both sports and work. Type I collagen fibrils are the primary components of tendons and carry most of the mechanical loads experienced by the tissue, however, knowledge of how load is transmitted between and within fibrils is limited. The presence of covalent enzymatic cross-links between collagen molecules is an important factor that has been shown to influence mechanical behavior of the tendons. To improve our understanding of how molecular bonds translate into tendon mechanics, we used an atomic force microscopy technique to measure the mechanical behavior of individual collagen fibrils loaded to failure. Fibrils from human patellar tendons, rat-tail tendons (RTTs), NaBH₄ reduced RTTs, and tail tendons of Zucker diabetic fat rats were tested. We found a characteristic three-phase stress-strain behavior in the human collagen fibrils. There was an initial rise in modulus followed by a plateau with reduced modulus, which was finally followed by an even greater increase in stress and modulus before failure. The RTTs also displayed the initial increase and plateau phase, but the third region was virtually absent and the plateau continued until failure. The importance of cross-link lability was investigated by NaBH₄ reduction of the rat-tail fibrils, which did not alter their behavior. These findings shed light on the function of cross-links at the fibril level, but further studies will be required to establish the underlying mechanisms.     

Lamellar orientation in human cornea in relation to mechanical properties

We have applied wide-angle X-ray scattering to the human cornea in order to quantify the relative number of stromal collagen fibrils directed along the two preferred corneal lamellar directions: superior-inferior and nasal-temporal. The data suggest that, on average, the two directions are populated in equal proportion at the corneal centre. Here, approximately one-third of the fibrils throughout the stromal depth tend to lie within a 45 degrees sector of the superior-inferior meridian, and similarly for the nasal-temporal direction. However, in some eyes we have measured significant differences between the two preferential fibril populations, with some corneas exhibiting as much as 25% more collagen in one direction than the other. Based on these findings, a mechanical model of the normal cornea may be envisaged, whereby the fibril tension in the underlying "background" of isotropically arranged collagen helps to balance the intraocular pressure; while the extra preferentially aligned fibrils take up the additional tensile stress along the superior-inferior and nasal-temporal meridians exerted by the rectus muscles and the orbicularis. It is possible that, through a direct impact on the elastic modulus of the tissue, an imbalance of superior-inferior and nasal-temporal fibrils in some eyes might affect corneal shape.     

Modeling the tensile mechanical behavior of bone along the longitudinal direction

The tensile stress-strain behavior of bone along its longitudinal axis is modeled by using a simple shear-lag theory, wherein, stresses and strains in a unit cell consisting of an organic matrix reinforced by overlapped mineral platelets are computed. It is assumed that loads are transferred between overlapped mineral-platelets by shear in the organic matrix. The mechanical behavior of bone in which the matrix partially or completely debonds from the sides of the overlapped mineral platelets (after an ultimate interfacial shear stress value is exceeded) is modeled. It is shown that the tensile mechanical behavior of bone can be modeled only by assuming little or no debonding of the organic from the mineral. A physical phenomenon that explains the tensile behavior of bone is, after the interfacial shear stress has reached a constant value over the length of the mineral platelets, the collagen molecules/microfibrils (with the associated mineral platelets) move relative to one another. The tensile stress-strain curve of bovine bone is modeled using this model. The theory predicts the mechanical behavior of the tissue in the elastic, yield and post-yield region. The ultimate strain and strengths are not predicted in the present model.     

The relationship of mechanical properties to morphology in patellar tendon autografts after posterior cruciate ligament replacement in sheep.

In a sheep model the posterior cruciate ligament (PCL) was replaced by a patellar tendon autograft (PTAG) using the central one-third of the ipsilateral patellar tendon (PT). The sheep were sacrificed at 16, 26, 52 and 104 weeks postoperation. The PTAG, and, as controls, the contralateral PCL and PT were harvested. These were examined using biomechanical testing as well as light and transmission electron microscopy, including immunohistological techniques. The material properties (maximum stress, elastic modulus) were compared to the morphological features. The cellular distribution, the distribution of glycosaminoglycans (GAGs), the collagen fibril diameter and the occurrence of Type III collagen were studied. Prior to transplantation, the PTAG was shown to be superior in maximum stress (57.2 +/- 5.5 MPa vs 41.3 +/- 1.9 MPa) and elastic modulus (368.8 +/- 49.3 MPa vs 172.3 +/- 14.6 MPa) to the PCL. The early decline in material properties of the PTAG (maximum stress 22% and elastic modulus 42% of the control) after free grafting paralleled a cell- and capillary-rich PTAG tissue with remnants of necrosis and a poorly organized extracellular matrix. Two years after implantation, with progressive alignment of the tissue matrix, maximum stress and elastic modulus acquired approximately 60 and 70% of the control, respectively. However, there was also an evidence of degenerative changes characterized by acellular areas, loss of the normal bundling pattern of collagen fibers and abnormal accumulation of GAGs. Ultrastructurally, there was a predominant shift to thin collagen fibrils in the PTAG compared to PCL and PT, both consisting of thick and thin collagen fibrils. Thin fibrils were demonstrated to be, in part, split thick fibrils as well as newly formed fibrils. Most of these thin fibrils revealed a positive reaction with antibodies to Type III collagen.     

Viscoelastic properties of isolated collagen fibrils

Understanding the viscoelastic behavior of collagenous tissues with complex hierarchical structures requires knowledge of the properties at each structural level. Whole tissues have been studied extensively, but less is known about the mechanical behavior at the submicron, fibrillar level. Using a microelectromechanical systems platform, in vitro coupled creep and stress relaxation tests were performed on collagen fibrils isolated from the sea cucumber dermis. Stress-strain-time data indicate that isolated fibrils exhibit viscoelastic behavior that could be fitted using the Maxwell-Weichert model. The fibrils showed an elastic modulus of 123 ± 46 MPa. The time-dependent behavior was well fit using the two-time-constant Maxwell-Weichert model with a fast time response of 7 ± 2 s and a slow time response of 102 ± 5 s. The fibrillar relaxation time was smaller than literature values for tissue-level relaxation time, suggesting that tissue relaxation is dominated by noncollagenous components (e.g., proteoglycans). Each specimen was tested three times, and the only statistically significant difference found was that the elastic modulus is larger in the first test than in the subsequent two tests, indicating that viscous properties of collagen fibrils are not sensitive to the history of previous tests.     

Stress transfer in collagen fibrils reinforcing connective tissues: effects of collagen fibril slenderness and relative stiffness

Unlike engineering fibre composite materials which comprise of fibres that are uniform cylindrical in shape, collagen fibrils reinforcing the proteoglycan-rich (PG) gel in the extra-cellular matrices (ECMs) of connective tissues are taper-ended (paraboloidal in shape). In an earlier paper we have discussed how taper of a fibril leads to an axial stress up-take which differs from that of a uniform cylindrical fibre and implications for fibril fracture. The present paper focuses on the influence of fibre aspect ratio, q (slenderness), and Young's modulus (stiffness), relative to that of the gel phase, E(R), on the magnitude of the axial tensile stresses generated within a fibril and wider implications on failure at tissue level. Fibre composite models were evaluated using finite element (FE) and mathematical analyses. When the applied force is low, there is elastic stress transfer between the PG gel and a fibril. FE modelling shows that the stress in a fibril increases with E(R) and q. At higher applied forces, there is plastic stress transfer. Mathematical modelling predicts that the stress in a fibril increases linearly with q. For small q values, fibrils may be regarded as fillers with little ability to provide tensile reinforcement. Large q values lead to high stress in a fibril. Such high stresses are beneficial provided they do not exceed the fracture stress of collagen. Modulus difference regulates the strain energy release density, u, for interfacial rupture; large E(R) not only leads to high stress in a fibril but also insures against interfacial rupture by raising the value of u.     

Mechanical properties of physiological and pathological models of collagen peptides investigated via steered molecular dynamics simulations

In this work we used molecular simulations to investigate the elastic properties of collagen single chain and triple helix with the aim of understanding its features starting from first principles. We analysed ideal collagen peptides, homotrimeric and heterotrimeric collagen type I and pathological models of collagen. Triple helices were found much more rigid than single chains, thus enlightening the important role of interchain stabilizing forces, like hydrophobic interaction and hydrogen bonds. We obtained Young's moduli close to 4.5GPa for the ideal model of collagen and for the physiological heterotrimer, while the physiological homotrimer presented a Young's modulus of 2.51GPa, that can be related to a mild form of Osteogenesis Imperfecta in which only the homotrimeric form of collagen type I is produced. Otherwise, the pathological model (presenting a glycine to alanine substitution) showed an elastic modulus of 4.32GPa, thus only slightly lower than the ideal model. This suggests that this mutation only slightly affects the mechanical properties of the collagen molecule, but possibly acts on an higher scale, such as the packing of collagen fibrils.     

Effect of bone mineral content on the tensile properties of cortical bone: experiments and theory

The effect of mineral volume fraction on the tensile mechanical properties of cortical bone tissue is investigated by theoretical and experimental means. The mineral content of plexiform, bovine bone was lowered by 18% and 29% by immersion in fluoride solutions for 3 days and 12 days, respectively. The elastic modulus, yield strength and ultimate strength of bone tissue decreased, while the ultimate strain increased with a decrease in mineral content. The mechanical behavior of bone tissue was modeled by using a micromechanical shear lag theory consisting of overlapped mineral platelets reinforcing the organic matrix. The decrease in yield stress, by the 0.002 offset method, of the fluoride treated bones were matched in the theoretical curves by lowering the shear yield stress of the organic matrix. The failure criterion used was based on failure stresses determined from a failure envelope (Mohr's circle), which was constructed using experimental data. It was found that the model predictions of elastic modulus got worse with a decrease in mineral content (being 7.9%, 17.2% and 33.0% higher for the control, 3-day and 12-day fluoride-treated bones). As a result, the developed theory could not fully predict the yield strain of bones with lowered mineral content, being 12.9% and 21.7% lower than the experimental values. The predicted ultimate stresses of the bone tissues with lower mineral contents were within +/- 10% of the experimental values while the ultimate strains were 12.7% and 26.3% lower than the experimental values. Although the model developed in this study did not take into account the presence of hierarchical structures, voids, orientation of collagen molecules and micro cracks, it still indicated that the mechanical properties of the organic matrix depend on bone mineral content.     

Mechanical properties of mineralized collagen fibrils as influenced by demineralization

Dentin and bone derive their mechanical properties from a complex arrangement of collagen type-I fibrils reinforced with nanocrystalline apatite mineral in extra- and intrafibrillar compartments. While mechanical properties have been determined for the bulk of the mineralized tissue, information on the mechanics of the individual fibril is limited. Here, atomic force microscopy was used on individual collagen fibrils to study structural and mechanical changes during acid etching. The characteristic 67 nm periodicity of gap zones was not observed on the mineralized fibril, but became apparent and increasingly pronounced with continuous demineralization. AFM-nanoindentation showed a decrease in modulus from 1.5 GPa to 50 MPa during acid etching of individual collagen fibrils and revealed that the modulus profile followed the axial periodicity. The nanomechanical data, Raman spectroscopy and SAXS support the hypothesis that intrafibrillar mineral etches at a substantially slower rate than the extrafibrillar mineral. These findings are relevant for understanding the biomechanics and design principles of calcified tissues derived from collagen matrices.     

Advanced Glycation End-Products Reduce Collagen Molecular Sliding to Affect Collagen Fibril Damage Mechanisms but Not Stiffness

Advanced glycation end-products (AGE) contribute to age-related connective tissue damage and functional deficit. The documented association between AGE formation on collagens and the correlated progressive stiffening of tissues has widely been presumed causative, despite the lack of mechanistic understanding. The present study investigates precisely how AGEs affect mechanical function of the collagen fibril – the supramolecular functional load-bearing unit within most tissues. We employed synchrotron small-angle X-ray scattering (SAXS) and carefully controlled mechanical testing after introducing AGEs in explants of rat-tail tendon using the metabolite methylglyoxal (MGO). Mass spectrometry and collagen fluorescence verified substantial formation of AGEs by the treatment. Associated mechanical changes of the tissue (increased stiffness and failure strength, decreased stress relaxation) were consistent with reports from the literature. SAXS analysis revealed clear changes in molecular deformation within MGO treated fibrils. Underlying the associated increase in tissue strength, we infer from the data that MGO modified collagen fibrils supported higher loads to failure by maintaining an intact quarter-staggered conformation to nearly twice the level of fibril strain in controls. This apparent increase in fibril failure resistance was characterized by reduced side-by-side sliding of collagen molecules within fibrils, reflecting lateral molecular interconnectivity by AGEs. Surprisingly, no change in maximum fibril modulus (2.5 GPa) accompanied the changes in fibril failure behavior, strongly contradicting the widespread assumption that tissue stiffening in ageing and diabetes is directly related to AGE increased fibril stiffness. We conclude that AGEs can alter physiologically relevant failure behavior of collagen fibrils, but that tissue level changes in stiffness likely occur at higher levels of tissue architecture.     

Influence of fibril taper on the function of collagen to reinforce extracellular matrix

Collagen fibrils provide tensile reinforcement for extracellular matrix. In at least some tissues, the fibrils have a paraboloidal taper at their ends. The purpose of this paper is to determine the implications of this taper for the function of collagen fibrils. When a tissue is subjected to low mechanical forces, stress will be transferred to the fibrils elastically. This process was modelled using finite element analysis because there is no analytical theory for elastic stress transfer to a non-cylindrical fibril. When the tissue is subjected to higher mechanical forces, stress will be transferred plastically. This process was modelled analytically. For both elastic and plastic stress transfer, a paraboloidal taper leads to a more uniform distribution of axial tensile stress along the fibril than would be generated if it were cylindrical. The tapered fibril requires half the volume of collagen than a cylindrical fibril of the same length and the stress is shared more evenly along its length. It is also less likely to fracture than a cylindrical fibril of the same length in a tissue subjected to the same mechanical force.     

Sensitivity analysis and parametric study of elastic properties of an unidirectional mineralized bone fibril-array using mean field methods

The key parameters determining the elastic properties of an unidirectional mineralized bone fibril-array decomposed in two further hierarchical levels are investigated using mean field methods. Modeling of the elastic properties of mineralized micro- and nanostructures requires accurate information about the underlying topology and the constituents’ material properties. These input data are still afflicted by great uncertainties and their influence on computed elastic constants of a bone fibril-array remains unclear. In this work, mean field methods are applied to model mineralized fibrils, the extra-fibrillar matrix and the resulting fibril-array. The isotropic or transverse isotropic elastic constants of these constituents are computed as a function of degree of mineralization, mineral distribution between fibrils and extra-fibrillar matrix, collagen stiffness and fibril volume fraction. The linear sensitivity of the elastic constants was assessed at a default set of the above parameters. The strain ratios between the constituents as well as the axial and transverse indentation moduli of the fibril-array were calculated for comparison with experiments. Results indicate that the degree of mineralization and the collagen stiffness dominate fibril-array elasticity. Interestingly, the stiffness of the extra-fibrillar matrix has a strong influence on transverse and shear moduli of the fibril-array. The axial strain of the intra-fibrillar mineral platelets is 30–90% of the applied fibril strain, depending on mineralization and collagen stiffness. The fibril-to-fibril-array strain ratio is essentially ~1. This study provides an improved insight in the parameters, which govern the fibril-array stiffness of mineralized tissues such as bone.     

Hierarchical modeling of the elastic properties of bone at submicron scales: the role of extrafibrillar mineralization

We model the elastic properties of bone at the level of mineralized collagen fibrils via step-by-step homogenization from the staggered arrangement of collagen molecules up to an array of parallel mineralized fibrils. A new model for extrafibrillar mineralization is proposed, assuming that the extrafibrillar minerals are mechanically equivalent to reinforcing rings coating each individual fibril. Our modeling suggests that no more than 30% of the total mineral content is extrafibrillar and the fraction of extrafibrillar minerals grows linearly with the overall degree of mineralization. It is shown that the extrafibrillar mineralization considerably reinforces the fibrils’ mechanical properties in the transverse directions and the fibrils’ shear moduli. The model predictions for the elastic moduli and constants are found to be in a good agreement with the experimental data reported in the literature.     

Possible role of decorin glycosaminoglycans in fibril to fibril force transfer in relative mature tendons--a computational study from molecular to microstructural level

Experimental studies on immature tendons have shown that the collagen fibril net is discontinuous. Manifold evidences, despite not being conclusive, indicate that mature tissue is discontinuous as well. According to composite theory, there is no requirement that the fibrils should extend from one end of the tissue to the other; indeed, an interfibrillar matrix with a low elastic modulus would be sufficient to guarantee the mechanical properties of the tendon. Possible mechanisms for the stress-transfer involve the interfibrillar proteoglycans and can be related to the matrix shear stress and to electrostatic non-covalent forces. Recent studies have shown that the glycosaminoglycans (GAGs) bound to decorin act like bridges between contiguous fibrils connecting adjacent fibril every 64-68 nm; this architecture would suggest their possible role in providing the mechanical integrity of the tendon structure. The present paper investigates the ability of decorin GAGs to transfer forces between adjacent fibrils. In order to test this hypothesis the stiffness of chondroitin-6-sulphate, a typical GAG associated to decorin, has been evaluated through the molecular mechanics approach. The obtained GAG stiffness is piecewise linear with an initial plateau at low strains (<800%) and a high stiffness region (3.1 x 10(-11)N/nm) afterwards. By introducing the calculated GAG stiffness in a multi-fibril model, miming the relative mature tendon architecture, the stress-strain behaviour of the collagen fibre was determined. The fibre incremental elastic modulus obtained ranges between 100 and 475 MPa for strains between 2% and 6%. The elastic modulus value depends directly on the fibril length, diameter and inversely on the interfibrillar distance. In particular, according to the obtained results, the length of the fibril is likely to play the major role in determining stiffness in mature tendons.     

Viscoelastic properties of collagen: synchrotron radiation investigations and structural model

Collagen type I is the most abundant structural protein in tendon, skin and bone, and largely determines the mechanical behaviour of these connective tissues. To obtain a better understanding of the relationship between structure and mechanical properties, tensile tests and synchrotron X-ray scattering have been carried out simultaneously, correlating the mechanical behaviour with changes in the microstructure. Because intermolecular cross-links are thought to have a great influence on the mechanical behaviour of collagen, we also carried out experiments using cross-link-deficient tail-tendon collagen from rats fed with beta-APN, in addition to normal controls. The load-elongation curve of tendon collagen has a characteristic shape with, initially, an increasing slope, corresponding to an increasing stiffness, followed by yielding and then fracture. Cross-link-deficient collagen produces a quite different curve with a marked plateau appearing in some cases, where the length of the tendon increases at constant stress. With the use of in situ X-ray diffraction, it was possible to measure simultaneously the elongation of the collagen fibrils inside the tendon and of the tendon as a whole. The overall strain of the tendon was always larger than the strain in the individual fibrils, which demonstrates that some deformation is taking place in the matrix between fibrils. Moreover, the ratio of fibril strain to tendon strain was dependent on the applied strain rate. When the speed of deformation was increased, this ratio increased in normal collagen but generally decreased in cross-link-deficient collagen, correlating to the appearance of a plateau in the force-elongation curve indicating creep. We proposed a simple structural model, which describes the tendon at a hierarchical level, where fibrils and interfibrillar matrix act as coupled viscoelastic systems. All qualitative features of the strain-rate dependence of both normal and cross-link-deficient collagen can be reproduced within this model. This complements earlier models that considered the next smallest level of hierarchy, describing the deformation of collagen fibrils in terms of changes in their molecular packing.     

Elastic modulus of hard tissues

This work aims at evaluating the elastic modulus of hard biological tissues by considering their staggered platelet micro-structure. An analytical expression for the effective modulus along the stagger direction is formulated using three non-dimensional structural variables. Structures with a single staggered hierarchy (e.g. collagen fibril) are first studied and predictions are compared with the experimental results and finite element simulations from the literature. A more complicated configuration, such as an array of fibrils, is analyzed next. Finally, a mechanical model is proposed for tooth dentin, in which variations in the multi-scale structural hierarchy are shown to significantly affect the macroscopic mechanical properties.