Mechanical properties of a collagen fibril under simulated degradation

The surface charge of electroactive materials governs cell behaviour through its effect on protein deposition

The precise mechanisms underlying the cellular response to static electric cues remain unclear, limiting the design and development of biomaterials that utilize this parameter to enhance specific biological behaviours. To gather information on this matter we have explored the interaction of collagen type-I, the most abundant mammalian extracellular protein, with poly(vinylidene fluoride) (PVDF), an electroactive polymer with great potential for tissue engineering applications. Our results reveal significant differences in collagen affinity, conformation, and interaction strength depending on the electric charge of the PVDF surface, which subsequently affects the behaviour of mesenchymal stem cells seeded on them. These findings highlight the importance of surface charge in the establishment of the material-protein interface and ultimately in the biological response to the material. STATEMENT OF SIGNIFICANCE: The development of new tissue engineering strategies relies heavily on the understanding of how biomaterials interact with biological tissues. Although several factors drive this process and their driving principles have been identified, the relevance and mechanism by which the surface potential influences cell behaviour is still unknown. In our study, we investigate the interaction between collagen, the most abundant component of the extracellular matrix, and poly(vinylidene fluoride) with varying surface charges. Our findings reveal substantial variations in the binding forces, structure and adhesion of collagen on the different surfaces, which collectively explain the differential cellular responses. By exposing these differences, our research fills a critical knowledge gap and paves the way for innovations in material design for advanced tissue regeneration strategies.

Dissipation and recovery in collagen fibrils under cyclic loading: A molecular dynamics study

The hysteretic behavior exhibited by collagen fibrils, when subjected to cyclic loading, is known to result in both dissipation as well as accumulation of residual strain. On subsequent relaxation, partial recovery has also been reported. Cross-links have been considered to play a key role in overall mechanical properties. Here, we modify an existing coarse-grained molecular dynamics model for collagen fibril with initially cross-linked collagen molecules, which is known to reproduce the response to uniaxial strain, by incorporating reformation of cross-links to allow for possible recovery of the fibril. Using molecular dynamics simulations, we show that our model successfully replicates the key features observed in experimental data, including the movement of hysteresis loops, the time evolution of residual strains and energy dissipation, as well as the recovery observed during relaxation. We also show that the characteristic cycle number, describing the approach toward steady state, has a value similar to that in experiments. We also emphasize the vital role of the degree of cross-linking on the key features of the macroscopic response to cyclic loading.

Knee joint biomechanics and cartilage damage prediction during landing: A hybrid MD-FE-musculoskeletal modeling

Understanding the mechanics behind knee joint injuries and providing appropriate treatment is crucial for improving physical function, quality of life, and employability. In this study, we used a hybrid molecular dynamics-finite element-musculoskeletal model to determine the level of loads the knee can withstand when landing from different heights (20, 40, 60 cm), including the height at which cartilage damage occurs. The model was driven by kinematics-kinetics data of asymptomatic subjects at the peak loading instance of drop landing. Our analysis revealed that as landing height increased, the forces on the knee joint also increased, particularly in the vastus muscles and medial gastrocnemius. The patellar tendon experienced more stress than other ligaments, and the medial plateau supported most of the tibial cartilage contact forces and stresses. The load was mostly transmitted through cartilage-cartilage interaction and increased with landing height. The critical height of 126 cm, at which cartilage damage was initiated, was determined by extrapolating the collected data using an iterative approach. Damage initiation and propagation were mainly located in the superficial layers of the tibiofemoral and patellofemoral cartilage. Finally, this study provides valuable insights into the mechanisms of landing-associated cartilage damage and could help limit joint injuries and improve training programs.

Mineralized Collagen Fibrils: An Essential Component in Determining the Mechanical Behavior of Cortical Bone

Bone comprises mechanically different materials in a specific hierarchical structure. Mineralized collagen fibrils (MCFs), represented by tropocollagen molecules and hydroxyapatite nanocrystals, are the fundamental unit of bone. The mechanical characterization of MCFs provides the unique adaptive mechanical competence to bone to withstand mechanical load. The structural and mechanical role of MCFs is critical in the deformation mechanisms of bone and the marvelous strength and toughness possessed by bone. However, the role of MCFs in the mechanical behavior of bone across multiple length scales is not fully understood. In the present study, we shed light upon the latest progress regarding bone deformation at multiple hierarchical levels and emphasize the role of MCFs during bone deformation. We propose the concept of hierarchical deformation of bone to describe the interconnected deformation process across multiple length scales of bone under mechanical loading. Furthermore, how the deterioration of bone caused by aging and diseases impairs the hierarchical deformation process of the cortical bone is discussed. The present work expects to provide insights on the characterization of MCFs in the mechanical properties of bone and lays the framework for the understanding of the multiscale deformation mechanics of bone.

Surrogate modeling of articular cartilage degradation to understand the synergistic role of MMP-1 and MMP-9: a case study

A characteristic feature of arthritic diseases is cartilage extracellular matrix (ECM) degradation, often orchestrated by the overexpression of matrix metalloproteinases (MMPs) and other proteases. The interplay between fibril level degradation and the tissue-level aggregate response to biomechanical loading was explored in this work by a computational multiscale cartilaginous model. We considered the relative abundance of collagenases (MMP-1) and gelatinases (MMP-9) in surrogate models, where the diffusion (spatial distribution) of these enzymes and the subsequent, co-localized fibrillar damage were spatially randomized with Latin Hypercube Sampling. The computational model was constructed by incorporating the results from prior molecular dynamics simulations (tensile test) of microfibril degradation into a hyper-elastoplastic fibril-reinforced cartilage model. Including MMPs-mediated collagen fibril-level degradation in computational models may help understand the ECM pathomechanics at the tissue level. The mechanics of cartilage tissue and fibril show variations in mechanical integrity depending on the different combinations of MMPs-1 and 9 with a concentration ratio of 1:1, 3:1, and 1:3 in simulated indentation tests. The fibril yield (local failure) was initiated at 20.2 ± 3.0 (%) and at 23.0 ± 2.8 (%) of bulk strain for col 1:gel 3 and col 3: gel 1, respectively. The reduction in failure stress (global response) was 39.8% for col 1:gel 3, 37.5% for col 1:gel 1, and 36.7% for col 3:gel 1 compared with the failure stress of the degradation free tissue. These findings indicate that cartilage's global and local mechanisms of failure largely depend on the relative abundance of the two key enzymes-collagenase (MMP-1) and gelatinase (MMP-9) and the spatial characteristics of diffusion across the layers of the cartilage ECM.

Kinetic model description of dissipation and recovery in collagen fibrils under cyclic loading

Collagen fibrils, when subjected to cyclic loading, are known to exhibit hysteretic behavior with energy dissipation that is partially recovered on relaxation. In this paper, we develop a kinetic model for a collagen fibril incorporating presence of hidden loops and stochastic fragmentation as well as reformation of sacrificial bonds. We show that the model reproduces well the characteristic features of reported experimental data on cyclic response of collagen fibrils, such as moving hysteresis loops, time evolution of residual strains and energy dissipation, recovery on relaxation, etc. We show that the approach to the steady state is controlled by a characteristic cycle number for both residual strain as well as energy dissipation and is in good agreement with reported existing experimental data.

Multiscale Characterization of Type I Collagen Fibril Stress–Strain Behavior under Tensile Load: Analytical vs. MD Approaches

Type I collagen is one of the most important proteins in the human body because of its role in providing structural support to the extracellular matrix of the connective tissues. Understanding its mechanical properties was widely investigated using experimental testing as well as molecular and finite element simulations. In this work, we present a new approach for defining the properties of the type I collagen fibrils by analytically formulating its response when subjected to a tensile load and investigating the effects of enzymatic crosslinks on the behavioral response. We reveal some of the shortcomings of the molecular dynamics (MD) method and how they affect the obtained stress–strain behavior of the fibril, and we prove that not only does MD underestimate the Young’s modulus and the ultimate tensile strength of the collagen fibrils, but also fails to detect the mechanics of some stretching phases of the fibril. We prove that non-crosslinked fibrils have three tension phases: (i) an initial elastic deformation corresponding to the collagen molecule uncoiling, (ii) a linear regime related to the stretching of the backbone of the tropocollagen molecules, and (iii) a plastic regime dominated by molecular sliding. We also show that for crosslinked fibrils, the second regime can be subdivided into three sub-regimes, and we define the properties of each regime. We also prove, analytically, the alleged MD quadratic relation between the ultimate tensile strength of the fibril and the concentration of enzymatic crosslinks (β).

Exploring the Mechanical Properties and Performance of Type-I Collagen at Various Length Scales: A Progress Report

Collagen is the basic protein of animal tissues and has a complex hierarchical structure. It plays a crucial role in maintaining the mechanical and structural stability of biological tissues. Over the years, it has become a material of interest in the biomedical industries thanks to its excellent biocompatibility and biodegradability and low antigenicity. Despite its significance, the mechanical properties and performance of pure collagen have been never reviewed. In this work, the emphasis is on the mechanics of collagen at different hierarchical levels and its long-term mechanical performance. In addition, the effect of hydration, important for various applications, was considered throughout the study because of its dramatic influence on the mechanics of collagen. Furthermore, the discrepancies in reports of the mechanical properties of collagenous tissues (basically composed of 20-30% collagen fibres) and those of pure collagen are discussed.

The Effect of Enzymatic Crosslink Degradation on the Mechanics of the Anterior Cruciate Ligament: A Hybrid Multi-Domain Model

The anterior cruciate ligament’s (ACL) mechanics is an important factor governing the ligament’s integrity and, hence, the knee joint’s response. Despite many investigations in this area, the cause and effect of injuries remain unclear or unknown. This may be due to the complexity of the direct link between macro- and micro-scale damage mechanisms. In the first part of this investigation, a three-dimensional coarse-grained model of collagen fibril (type I) was developed using a bottom-up approach to investigate deformation mechanisms under tensile testing. The output of this molecular level was used later to calibrate the parameters of a hierarchical multi-scale fibril-reinforced hyper-elastoplastic model of the ACL. Our model enabled us to determine the mechanical behavior of the ACL as a function of the basic response of the collagen molecules. Modeled elastic response and damage distribution were in good agreement with the reported measurements and computational investigations. Our results suggest that degradation of crosslink content dictates the loss of the stiffness of the fibrils and, hence, damage to the ACL. Therefore, the proposed computational frame is a promising tool that will allow new insights into the biomechanics of the ACL.

Multiscale modeling of knee ligament biomechanics

Knee connective tissues are mainly responsible for joint stability and play a crucial role in restraining excessive motion during regular activities. The damage mechanism of these tissues is directly linked to the microscale collagen level. However, this mechanical connection is still unclear. During this investigation, a multiscale fibril-reinforced hyper-elastoplastic model was developed and statistically calibrated. The model is accounting for the structural architecture of the soft tissue, starting from the tropocollagen molecule that forms fibrils to the whole soft tissue. Model predictions are in agreement with the results of experimental and numerical studies. Further, damage initiation and propagation in the collagen fiber were computed at knee ligaments and located mainly in the superficial layers. Results indicated higher crosslink density required higher tensile stress to elicit fibril damage. This approach is aligned with a realistic simulation of a damaging process and repair attempt. To the best of our knowledge, this is the first model published in which the connective tissue stiffness is simultaneously predicted by encompassing the mesoscopic scales between the molecular and macroscopic levels.

Collagen denaturation is initiated upon tissue yield in both positional and energy-storing tendons

Tendons are collagenous soft tissues that transmit loads between muscles and bones. Depending on their anatomical function, tendons are classified as positional or energy-storing with differing biomechanical and biochemical properties. We recently demonstrated that during monotonic stretch of positional tendons, permanent denatured collagen begins accumulating upon departing the linear region of the stress-strain curve. However, it is unknown if this observation is true during mechanical overload of other types of tendons. Therefore, the purpose of this study was to investigate the onset of collagen denaturation relative to applied strain, and whether it differs between the two tendon types. Rat tail tendon (RTT) fascicles and rat flexor digitorum longus (FDL) tendons represented positional and energy-storing tendons, respectively. The samples were stretched to incremental levels of strain, then stained with fluorescently labeled collagen hybridizing peptides (CHPs); the CHP fluorescence was measured to quantify denatured collagen. Denatured collagen in both positional and energy-storing tendons began to increase at the yield strain, upon leaving the linear region of the stress-strain curve as the sample started to permanently deform. Despite significant differences between the two tendon types, it appears that collagen denaturation is initiated at tissue yield during monotonic stretch, and the fundamental mechanism of failure is the same for the two types of tendons. At tissue failure, positional tendons had double the percentage of denatured collagen compared to energy-storing tendons, with no difference between 0% control groups. These results help to elucidate the etiology of subfailure injury and rupture in functionally distinct tendons.

On collagen fiber morphoelasticity and homeostatic remodeling tone

A variety of biochemical and physical processes participate in the creation and maintenance of collagen in biological tissue. Under mechanical stimuli these collagen fibers undergo continuous processes of morphoelastic change. The model presented here is motivated by experimental reports of stretch-stabilization of the collagen fibers to enzymatic degradation. The fiber structure is modeled in terms of a fiber density evolution that is regulated by means of a fixed creation rate and a mechano-sensitive dissolution rate. The theory accounts for the possibly different natural configurations of the fiber unit constituents and the ground substance matrix. It also generalizes previous theoretical descriptions so as to account for finite survival times of the individual fiber units. Special consideration is given to steady state fiber-remodeling processes in which fiber creation and dissolution are in balance. Fiber assembly processes that involve prestretching the fiber constituents yield a homeostatic stress response with a characteristic fiber tone. Fiber density returns to homeostasis after mechanical disruption when sufficient time has passed.

Stress-Swelling Finite Element Modeling of Cervical Response With Homeostatic Collagen Fiber Distributions

During pregnancy, the cervix experiences significant mechanical property change due to tissue swelling, and to ongoing changes in the collagen content. In this paper, we model how these two effects contribute to cervical deformation as the pressure load on top of the cervix increases. The cervix and its surrounding supporting ligaments are taken into consideration in the resulting mechanical analysis. The cervix itself is treated as a multilayered tube-like structure, with layer-specific collagen orientation. The cervical tissue in each layer is treated in terms of a collagen constituent that remodels with time within a ground substance matrix that experiences swelling. The load and swelling are taken to change sufficiently slowly so that the collagen properties at any instant can be regarded as being in a state of homeostasis. Among other things, the simulations show how the luminal cross-sectional area varies along its length as a function of pressure and swelling. In general, an increase in pressure causes an overall shortening of the lumen while an increase in swelling has the opposite effect.

Comparison of the Structural Characteristics of Native Collagen Fibrils Derived from Bovine Tendons using Two Different Methods: Modified Acid-Solubilized and Pepsin-Aided Extraction

Native collagen fibrils (CF) were successfully extracted from bovine tendons using two different methods: modified acid-solubilized extraction for A-CF and pepsin-aided method for P-CF. The yields of A-CF and P-CF were up to 64.91% (±1.07% SD) and 56.78% (±1.22% SD) (dry weight basis), respectively. The analyses of both amino acid composition and sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) confirmed that A-CF and P-CF were type I collagen fibrils. Both A-CF and P-CF retained the intact crystallinity and integrity of type I collagen’s natural structure by FTIR spectra, circular dichroism spectroscopy (CD) and X-ray diffraction detection. The aggregation structures of A-CF and P-CF were displayed by UV–Vis. However, A-CF showed more intact aggregation structure than P-CF. Microstructure and D-periodicities of A-CF and P-CF were observed (SEM and TEM). The diameters of A-CF and P-CF are about 386 and 282 nm, respectively. Although both A-CF and P-CF were theoretically concordant with the Schmitt hypothesis, A-CF was of evener thickness and higher integrity in terms of aggregation structure than P-CF. Modified acid-solubilized method provides a potential non-enzyme alternative to extract native collagen fibrils with uniform thickness and integral aggregation structure.

Transcriptome Profile of Human Fibroblasts in an Ex Vivo Culture

Implantation of autologous fibroblasts is a method used to correct age-related changes in facial skin. The aim of this study was to establish the optimal population of cultured human fibroblasts according to the organization of the extracellular matrix in the dermis. Transcriptome profile analysis of cells derived from three consecutive passages indicated that fibroblasts after the second passage were the population with the greatest number of upregulated genes encoding the critical biological processes responsible for skin regeneration, such as extracellular matrix organization, collagen fibril organization, and cell adhesion. Furthermore, genes encoding proteinases responsible for the degradation of dermal extracellular matrix proteins were noticeably downregulated at this stage of culture. Autologous fibroblasts seem to be an optimal and safe biological filler for the renewal of all skin structures.

Cyclophilin B control of lysine post-translational modifications of skin type I collagen

Covalent intermolecular cross-linking of collagen is essential for tissue stability. Recent studies have demonstrated that cyclophilin B (CypB), an endoplasmic reticulum (ER)-resident peptidyl-prolyl cis-trans isomerase, modulates lysine (Lys) hydroxylation of type I collagen impacting cross-linking chemistry. However, the extent of modulation, the molecular mechanism and the functional outcome in tissues are not well understood. Here, we report that, in CypB null (KO) mouse skin, two unusual collagen cross-links lacking Lys hydroxylation are formed while neither was detected in wild type (WT) or heterozygous (Het) mice. Mass spectrometric analysis of type I collagen showed that none of the telopeptidyl Lys was hydroxylated in KO or WT/Het mice. Hydroxylation of the helical cross-linking Lys residues was almost complete in WT/Het but was markedly diminished in KO. Lys hydroxylation at other sites was also lower in KO but to a lesser extent. A key glycosylation site, α1(I) Lys-87, was underglycosylated while other sites were mostly overglycosylated in KO. Despite these findings, lysyl hydroxylases and glycosyltransferase 25 domain 1 levels were significantly higher in KO than WT/Het. However, the components of ER chaperone complex that positively or negatively regulates lysyl hydroxylase activities were severely reduced or slightly increased, respectively, in KO. The atomic force microscopy-based nanoindentation modulus were significantly lower in KO skin than WT. These data demonstrate that CypB deficiency profoundly affects Lys post-translational modifications of collagen likely by modulating LH chaperone complexes. Together, our study underscores the critical role of CypB in Lys modifications of collagen, cross-linking and mechanical properties of skin.

Deficiency of cyclophilin B (CypB), an endoplasmic reticulum-resident peptidyl-prolyl cis-trans isomerase, causes recessive osteogenesis imperfecta type IX, resulting in defective connective tissues. Recent studies using CypB null mice revealed that CypB modulates lysine hydroxylation of type I collagen impacting collagen cross-linking. However, the extent of modulation, the molecular mechanism and the effect on tissue properties are not well understood. In the present study, we show that CypB deficiency in mouse skin results in the formation of unusual collagen cross-links, aberrant tissue formation, altered levels of lysine modifying enzymes and their chaperones, and impaired mechanical property. These findings highlight an essential role of CypB in collagen post-translational modifications which are critical in maintaining the structure and function of connective tissues.

Effect of collagenase-gelatinase ratio on the mechanical properties of a collagen fibril: a combined Monte Carlo-molecular dynamics study

Loading in cartilage is supported primarily by fibrillar collagen, and damage will impair the function of the tissue, leading to pathologies such as osteoarthritis. Damage is initiated by two types of matrix metalloproteinases, collagenase and gelatinase, that cleave and denature the collagen fibrils in the tissue. Experimental and modeling studies have revealed insights into the individual contributions of these two types of MMPs, as well as the mechanical response of intact fibrils and fibrils that have experienced random surface degradation. However, no research has comprehensively examined the combined influences of collagenases and gelatinases on collagen degradation nor studied the mechanical consequences of biological degradation of collagen fibrils. Such preclinical examinations are required to gain insights into understanding, treating, and preventing degradation-related cartilage pathology. To develop these insights, we use sequential Monte Carlo and molecular dynamics simulations to probe the effect of enzymatic degradation on the structure and mechanics of a single collagen fibril. We find that the mechanical response depends on the ratio of collagenase to gelatinase—not just the amount of lost fibril mass—and we provide a possible mechanism underlying this phenomenon. Overall, by characterizing the combined influences of collagenases and gelatinases on fibril degradation and mechanics at the preclinical research stage, we gain insights that may facilitate the development of targeted interventions to prevent the damage and loss of mechanical integrity that can lead to cartilage pathology.

Biopolymer nanofibrils: structure, modeling, preparation, and applications

Biopolymer nanofibrils exhibit exceptional mechanical properties with a unique combination of strength and toughness, while also presenting biological functions that interact with the surrounding environment. These features of biopolymer nanofibrils profit from their hierarchical structures that spun angstrom to hundreds of nanometer scales. To maintain these unique structural features and to directly utilize these natural supramolecular assemblies, a variety of new methods have been developed to produce biopolymer nanofibrils. In particular, cellulose nanofibrils (CNFs), chitin nanofibrils (ChNFs), silk nanofibrils (SNFs) and collagen nanofibrils (CoNFs), as the four most abundant biopolymer nanofibrils on earth, have been the focus of research in recent years due to their renewable features, wide availability, low-cost, biocompatibility, and biodegradability. A series of top-down and bottom-up strategies have been accessed to exfoliate and regenerate these nanofibrils for versatile advanced applications. In this review, we first summarize the structures of biopolymer nanofibrils in nature and outline their related computational models with the aim of disclosing fundamental structure-property relationships in biological materials. Then, we discuss the underlying methods used for the preparation of CNFs, ChNFs, SNF and CoNFs, and discuss emerging applications for these biopolymer nanofibrils.