A hyperelastic fibre-reinforced continuum model of healing tendons with distributed collagen fibre orientations

The alteration of the structure and macroscopic mechanical response of porcine patellar tendon by elastase digestion

Background: The treatment of patellar tendon injury has always been an unsolved problem, and mechanical characterization is very important for its repair and reconstruction. Elastin is a contributor to mechanics, but it is not clear how it affects the elasticity, viscoelastic properties, and structure of patellar tendon. Methods: The patellar tendons from six fresh adult experimental pigs were used in this study and they were made into 77 samples. The patellar tendon was specifically degraded by elastase, and the regional mechanical response and structural changes were investigated by: (1) Based on the previous study of elastase treatment conditions, the biochemical quantification of collagen, glycosaminoglycan and total protein was carried out; (2) The patellar tendon was divided into the proximal, central, and distal regions, and then the axial tensile test and stress relaxation test were performed before and after phosphate-buffered saline (PBS) or elastase treatment; (3) The dynamic constitutive model was established by the obtained mechanical data; (4) The structural relationship between elastin and collagen fibers was analyzed by two-photon microscopy and histology. Results: There was no statistical difference in mechanics between patellar tendon regions. Compared with those before elastase treatment, the low tensile modulus decreased by 75%-80%, the high tensile modulus decreased by 38%-47%, and the transition strain was prolonged after treatment. For viscoelastic behavior, the stress relaxation increased, the initial slope increased by 55%, the saturation slope increased by 44%, and the transition time increased by 25% after enzyme treatment. Elastin degradation made the collagen fibers of patellar tendon become disordered and looser, and the fiber wavelength increased significantly. Conclusion: The results of this study show that elastin plays an important role in the mechanical properties and fiber structure stability of patellar tendon, which supplements the structure-function relationship information of patellar tendon. The established constitutive model is of great significance to the prediction, repair and replacement of patellar tendon injury. In addition, human patellar tendon has a higher elastin content, so the results of this study can provide supporting information on the natural properties of tendon elastin degradation and guide the development of artificial patellar tendon biomaterials.

Adaptation of Fibril-Reinforced Poroviscoelastic Properties in Rabbit Collateral Ligaments 8 Weeks After Anterior Cruciate Ligament Transection

Ligaments of the knee provide stability and prevent excessive motions of the joint. Rupture of the anterior cruciate ligament (ACL), a common sports injury, results in an altered loading environment for other tissues in the joint, likely leading to their mechanical adaptation. In the collateral ligaments, the patterns and mechanisms of biomechanical adaptation following ACL transection (ACLT) remain unknown. We aimed to characterize the adaptation of elastic and viscoelastic properties of the lateral and medial collateral ligaments eight weeks after ACLT. Unilateral ACLT was performed in six rabbits, and collateral ligaments were harvested from transected and contralateral knee joints after eight weeks, and from an intact control group (eight knees from four animals). The cross-sectional areas were measured with micro-computed tomography. Stepwise tensile stress-relaxation testing was conducted up to 6% final strain, and the elastic and viscoelastic properties were characterized with a fibril-reinforced poroviscoelastic material model. We found that the cross-sectional area of the collateral ligaments in the ACL transected knees increased, the nonlinear elastic collagen network modulus of the LCL decreased, and the amount of fast relaxation in the MCL decreased. Our results indicate that rupture of the ACL leads to an early adaptation of the elastic and viscoelastic properties of the collagen fibrillar network in the collateral ligaments. These adaptations may be important to consider when evaluating whole knee joint mechanics after ACL rupture, and the results aid in understanding the consequences of ACL rupture on other tissues.

Mechano-biological and bio-mechanical pathways in cutaneous wound healing

Injuries to the skin heal through coordinated action of fibroblast-mediated extracellular matrix (ECM) deposition, ECM remodeling, and wound contraction. Defects involving the dermis result in fibrotic scars featuring increased stiffness and altered collagen content and organization. Although computational models are crucial to unravel the underlying biochemical and biophysical mechanisms, simulations of the evolving wound biomechanics are seldom benchmarked against measurements. Here, we leverage recent quantifications of local tissue stiffness in murine wounds to refine a previously-proposed systems-mechanobiological finite-element model. Fibroblasts are considered as the main cell type involved in ECM remodeling and wound contraction. Tissue rebuilding is coordinated by the release and diffusion of a cytokine wave, e.g. TGF-β, itself developed in response to an earlier inflammatory signal triggered by platelet aggregation. We calibrate a model of the evolving wound biomechanics through a custom-developed hierarchical Bayesian inverse analysis procedure. Further calibration is based on published biochemical and morphological murine wound healing data over a 21-day healing period. The calibrated model recapitulates the temporal evolution of: inflammatory signal, fibroblast infiltration, collagen buildup, and wound contraction. Moreover, it enables in silico hypothesis testing, which we explore by: (i) quantifying the alteration of wound contraction profiles corresponding to the measured variability in local wound stiffness; (ii) proposing alternative constitutive links connecting the dynamics of the biochemical fields to the evolving mechanical properties; (iii) discussing the plausibility of a stretch- vs. stiffness-mediated mechanobiological coupling. Ultimately, our model challenges the current understanding of wound biomechanics and mechanobiology, beside offering a versatile tool to explore and eventually control scar fibrosis after injury.

Insights into the mechanics of solid conical microneedle array insertion into skin using the finite element method

In order to develop optimum microneedle designs, researchers must first develop robust, repeatable and adaptable test methods which are representative of in vivo conditions. However, there is a lack of experimental tools which can accurately comparatively interrogate functional microneedle penetration of tissue. In this study, we seek to develop a state of the art finite element model of microneedle insertion into and penetration of human skin. The developed model employs a 3D hyperelastic, anisotropic pre-stressed multi-layered material which more accurately reflects in vivo skin conditions, while the microneedle is modeled as an array, which can capture the influence of adjacent microneedles on the overall response. Using the developed finite element model, we highlight the importance of accurate computational modeling which can decipher the mechanics of microneedle insertion, including the influence of its position within an array and how it correlates well with experimental observations. In particular, we have concluded that, for our model microneedle array, increasing skin pretension from 0 to 10% strain reduces the penetration force by 13%, ultimate local deformation about the microneedle by 22% and the ultimate penetration efficiency by 15%. We have also concluded that the presence of a base plate limits the penetration efficiency by up to 24%, while the penetration efficiency across a 5 × 1 microneedle array may vary by 27%. This model elucidates, for the first time, the combined effects of skin tension and needle geometry on accurately predicting microneedle penetration efficiency. STATEMENT OF SIGNIFICANCE: Microneedles arrays (MNAs) are medical devices with microscale protrusions, typically designed to penetrate the outermost layer of the skin, that upon optimisation, could lead to disruptive minimally-invasive disease management. However, the mechanics of MNA insertion are complex, due in part to a 'bed of nails' effect, and difficult to elucidate experimentally. Therefore, comparisons between designs, functional assessment of production batches and ultimately the likelihood of clinical translation are challenging to predict. Here, we have develop the most sophisticated in silico model of MNA insertion into pre-tensioned human skin to predict the extent of MNA penetration and therefore the likelihood of successful therapeutic delivery. Researchers can customise this model to predict the penetration efficiency of any MNA design.

Comparison of material models for anterior cruciate ligament in tension: from poroelastic to a novel fibril-reinforced nonlinear composite model

Computational models of the knee joint are useful for evaluating stresses and strains within the joint tissues. However, the outcome of those models is sensitive to the material model and material properties chosen for ligaments, the collagen reinforced tissues connecting bone to bone. The purpose of this study was to investigate different compositionally motivated material models and further to develop a model that can accurately reproduce experimentally measured stress-relaxation data of bovine anterior cruciate ligament (ACL). Tensile testing samples were extracted from ACLs of bovine knee joints (N = 10) and subjected to a three-step stress-relaxation test at the toe region. Data from the experiments was averaged and one average finite element model was generated to replicate the experiment. Poroelastic and different fibril-reinforced poro(visco)elastic material models were applied, and their material parameters were optimized to reproduce the experimental force-time response. Material models with only fluid flow mediated relaxation were not able to capture the stress-relaxation behavior (R² = 0.806, 0.803 and 0.938). The inclusion of the viscoelasticity of the fibrillar network improved the model prediction (R² = 0.978 and 0.976), but the complex stress-relaxation behavior was best captured by a poroelastic model with a nonlinear two-relaxation-time strain-recruited viscoelastic fibrillar network (R² = 0.997). The results suggest that in order to replicate the multi-step stress-relaxation behavior of ACL in tension, the fibrillar network formulation should include the complex nonlinear viscoelastic phenomena.

Understanding how reduced loading affects Achilles tendon mechanical properties using a fibre-reinforced poro-visco-hyper-elastic model

Understanding tendon mechanobiology is important for gaining insight into the development of tendon pathology and subsequent repair processes. The aim of this study was to investigate how experimentally observed mechanobiological adaptation of rat Achilles tendons translate to changes in constitutive mechanical properties and biomechanical behavior. In addition, we assessed the ability of the model to simulate tendon creep and stress-relaxation. A three dimensional finite element framework of rat Achilles tendon was implemented with a fibre-reinforced poro-visco-hyper-elastic constitutive model. Stress-relaxation and creep data from Achilles tendons of Sprague Dawley rats that had been subjected to both daily loading and a period of reduced loading were used to determine the constitutive properties of the tendons. Our results showed that the constitutive model captures creep and stress-relaxation data from rat Achilles tendons for both loaded and unloaded tendons with good accuracy (normalized root mean square error between model and experimental data were 0.010-0.027). Only when the model parameters were fitted to data from both mechanical tests simultaneously, were we able to also capture similar increase in elastic energy (increased stiffness) and decreased viscoelasticity in response to unloading, as was reported experimentally. Our study is the first to show that experimentally observed mechanobiological changes in tendon biomechanics, such as stiffness and viscoelasticity, can be designated to mechanical quantities in a constitutive model. Further investigation in this direction has potential to discriminate tissue components responsible for specific biomechanical response, and enable targeted treatment strategies for tendon health.

Evaluation of microstructurally motivated constitutive models to describe age-dependent tendon healing

Tendon injuries are common to all ages. Injured tendons typically do not recover full functionality. The amount and organization of tendon constituents dictate their mechanical properties. The impact of changes in these constituents during (patho)physiologic processes (e.g., aging and healing) are not fully understood. Toward this end, microstructurally motivated strain energy functions (SEFs) offer insight into underlying mechanisms of age-dependent healing. Several SEFs have been adapted for tendon; however, most are phenomenological. Therefore, the aims of this study are: (1) evaluate the descriptive capability of SEFs in age-dependent murine patellar tendon healing and (2) identify a SEF for implementation in a growth and remodeling (G&R) model. To accomplish these aims, models were fitted to patellar tendon tensile data from multiple age groups and post-injury timepoints. Model sensitivity to parameters and the determinability of the parameters were assessed. A two-way analysis of variance was used to identify changes in parameters and the feasibility of implementing each model into a G&R model is discussed. The evaluated SEFs exhibited adequate descriptive capability. Parameter determinability and sensitivity analysis, however, highlighted the need for additional data to inform and validate the models to increase physiologic relevance and enable G&R model formulation to determine underlying mechanisms of age-dependent healing. This work, as a first, evaluated changes in tendon mechanical properties both as functions of age and injury in an age-dependent manner using microstructurally motivated models, highlights inherent dependencies between parameters of widely used hyperelastic models, and identified unique post-injury behavior by the aging group compared to the mature and aged groups.

Modeling the effect of collagen fibril alignment on ligament mechanical behavior

Ligament mechanical behavior is primarily regulated by fibrous networks of type I collagen. Although these fibrous networks are typically highly aligned, healthy and injured ligament can also exhibit disorganized collagen architecture. The objective of this study was to determine whether variations in the collagen fibril network between neighboring ligaments can predict observed differences in mechanical behavior. Ligament specimens from two regions of bovine fetlock joints, which either exhibited highly aligned or disorganized collagen fibril networks, were mechanically tested in uniaxial tension. Confocal microscopy and FiberFit software were used to quantify the collagen fibril dispersion and mean fibril orientation in the mechanically tested specimens. These two structural parameters served as inputs into an established hyperelastic constitutive model that accounts for a continuous distribution of planar fibril orientations. The ability of the model to predict differences in the mechanical behavior between neighboring ligaments was tested by (1) curve fitting the model parameters to the stress response of the ligament with highly aligned fibrils and then (2) using this model to predict the stress response of the ligament with disorganized fibrils by only changing the parameter values for fibril dispersion and mean fibril orientation. This study found that when using parameter values for fibril dispersion and mean fibril orientation based on confocal imaging data, the model strongly predicted the average stress response of ligaments with disorganized fibrils ([Formula: see text]); however, the model only successfully predicted the individual stress response of ligaments with disorganized fibrils in half the specimens tested. Model predictions became worse when parameters for fibril dispersion and mean fibril orientation were not based on confocal imaging data. These findings emphasize the importance of collagen fibril alignment in ligament mechanics and help advance a mechanistic understanding of fibrillar networks in healthy and injured ligament.

Experimental evaluation of multiscale tendon mechanics

Tendon's primary function is a mechanical link between muscle and bone. The hierarchical structure of tendon and specific compositional constituents are believed to be critical for proper mechanical function. With increased appreciation for tendon importance and the development of various technological advances, this review paper summarizes recent experimental approaches that have been used to study multiscale tendon mechanics, includes an overview of studies that have evaluated the role of specific tissue constituents, and also proposes challenges/opportunities facing tendon study. Tendon has been demonstrated to have specific structural characteristics (e.g., multi-level hierarchy, crimp pattern, helix) and complex mechanical properties (e.g., non-linearity, anisotropy, viscoelasticity). Physical mechanisms including uncrimping, fiber sliding, and collagen reorganization have been shown to govern tendon mechanical responses under both static and dynamic loading. Several tendon constituents with relatively small quantities have been suggested to play a role in its mechanics, although some results are conflicting. Further research should be performed to understand the interplay and communication of tendon mechanical properties across levels of the hierarchical structure, and further show how each of these components contribute to tendon mechanics. The studies summarized and discussed in this review have helped elucidate important aspects of multiscale tendon mechanics, which is a prerequisite for analyzing stress/strain transfer between multiple scales and identifying key principles of mechanotransduction. This information could further facilitate interpreting the functional diversity of tendons from different species, different locations, and even different developmental stages, and then better understand and identify fundamental concepts related to tendon degeneration, disease, and healing. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:1353-1365, 2017.

Mechanobiological modelling of tendons: Review and future opportunities

Tendons are adapted to carry large, repeated loads and are clinically important for the maintenance of musculoskeletal health in an increasing, actively ageing population, as well as in elite athletes. Tendons are known to adapt to mechanical loading. Also, their healing and disease processes are highly sensitive to mechanical load. Computational modelling approaches developed to capture this mechanobiological adaptation in tendons and other tissues have successfully addressed many important scientific and clinical issues. The aim of this review is to identify techniques and approaches that could be further developed to address tendon-related problems. Biomechanical models are identified that capture the multi-level aspects of tendon mechanics. Continuum whole tendon models, both phenomenological and microstructurally motivated, are important to estimate forces during locomotion activities. Fibril-level microstructural models are documented that can use these estimated forces to detail local mechanical parameters relevant to cell mechanotransduction. Cell-level models able to predict the response to such parameters are also described. A selection of updatable mechanobiological models is presented. These use mechanical signals, often continuum tissue level, along with rules for tissue change and have been applied successfully in many tissues to predict in vivo and in vitro outcomes. Signals may include scalars derived from the stress or strain tensors, or in poroelasticity also fluid velocity, while adaptation may be represented by changes to elastic modulus, permeability, fibril density or orientation. So far, only simple analytical approaches have been applied to tendon mechanobiology. With the development of sophisticated computational mechanobiological models in parallel with reporting more quantitative data from in vivo or clinical mechanobiological studies, for example, appropriate imaging, biochemical and histological data, this field offers huge potential for future development towards clinical applications.

Comparison of structural anisotropic soft tissue models for simulating Achilles tendon tensile behaviour

The incidence of tendon injury (tendinopathy) has increased over the past decades due to greater participation in sports and recreational activities. But little is known about the aetiology of tendon injuries because of our limited knowledge in the complex structure-function relationship in tendons. Computer models can capture the biomechanical behaviour of tendons and its structural components, which is essential for understanding the underlying mechanisms of tendon injuries. This study compares three structural constitutive material models for the Achilles tendon and discusses their application on different biomechanical simulations. The models have been previously used to describe cardiovascular tissue and articular cartilage, and one model is novel to this study. All three constitutive models captured the tensile behaviour of rat Achilles tendon (root mean square errors between models and experimental data are 0.50-0.64). They further showed that collagen fibres are the main load-bearing component and that the non-collagenous matrix plays a minor role in tension. By introducing anisotropic behaviour also in the non-fibrillar matrix, the new biphasic structural model was also able to capture fluid exudation during tension and high values of Poisson׳s ratio that is reported in tendon experiments.