A tissue adaptation model based on strain-dependent collagen degradation and contact-guided cell traction

Signatures of disease progression in knee osteoarthritis: insights from an integrated multi-scale modeling approach, a proof of concept

Introduction: Knee osteoarthritis (KOA) is characterized by articular cartilage degeneration. It has been widely accepted that the mechanical joint environment plays a significant role in the onset and progression of this disease. In silico models have been used to study the interplay between mechanical loading and cartilage degeneration, hereby relying mainly on two key mechanoregulatory factors indicative of collagen degradation and proteoglycans depletion. These factors are the strain in collagen fibril direction (SFD) and maximum shear strain (MSS) respectively. Methods: In this study, a multi-scale in silico modeling approach was used based on a synergy between musculoskeletal and finite element modeling to evaluate the SFD and MSS. These strains were evaluated during gait based on subject-specific gait analysis data collected at baseline (before a 2-year follow-up) for a healthy and progressive early-stage KOA subject with similar demographics. Results: The results show that both SFD and MSS factors allowed distinguishing between a healthy subject and a KOA subject, showing progression at 2 years follow-up, at the instance of peak contact force as well as during the stance phase of the gait cycle. At the peak of the stance phase, the SFD were found to be more elevated in the KOA patient with the median being 0.82% higher in the lateral and 0.4% higher in the medial compartment of the tibial cartilage compared to the healthy subject. Similarly, for the MSS, the median strains were found to be 3.6% higher in the lateral and 0.7% higher in the medial tibial compartment of the KOA patient compared to the healthy subject. Based on these intersubject SFD and MSS differences, we were additionally able to identify that the tibial compartment of the KOA subject at risk of progression. Conclusion/discussion: We confirmed the mechanoregulatory factors as potential biomarkers to discriminate patients at risk of disease progression. Future studies should evaluate the sensitivity of the mechanoregulatory factors calculated based on this multi-scale modeling workflow in larger patient and control cohorts.

Computational model of articular cartilage regeneration induced by scaffold implantation in vivo

Computational models allow to explain phenomena that cannot be observed through an animal model, such as the strain and stress states which can highly influence regeneration of the tissue. For this purpose, we have developed a simulation tool to determine the mechanical conditions provided by the polymeric scaffold. The computational model considered the articular cartilage, the subchondral bone, and the scaffold. All materials were modeled as poroelastic, and the cartilage had linear-elastic oriented collagen fibers. This model was able to explain the remodeling process that subchondral bone goes through, and how the scaffold allowed the conditions for cartilage regeneration. These results suggest that the use of scaffolds might lead the cartilaginous tissue growth in vivo by providing a better mechanical environment. Moreover, the developed computational model demonstrated to be useful as a tool prior experimental in vivo studies, by predicting the possible outcome of newly proposed treatments allowing to discard approaches that might not bring good results.

A linear viscoelastic model to understand skin mechanical behaviour and for cosmetic formulation design

OBJECTIVES

Skin is arguably one of the most important organs that plays an active role in our everyday biological functions after brain. Owing to the wide range of applications in medicine, cosmetics industry and more recently robotics, skin research has gained tremendous attention with respect to its mechanical behaviour. Various macro modelling approaches are available for modelling skin's mechanical behaviour. The objective of this paper is to study skin's mechanical property change with age and demonstrate anti-ageing effects of cosmetic formulations from skin mechanical property change perspective.

METHODS

In this study, skin's mechanical behaviour was modelled using a 1D linear viscoelastic phenomenological model and the model was validated using two sets of experimentally observed skin data (strain, stress relaxation and cyclical loading). The model was further modified to study the effect of the presence of a thin layer of cosmetic polymer and to demonstrate anti-ageing effects of the cosmetic polymer from the perspective of change in the mechanical behaviour of skin with cosmetic layer.

RESULTS

The estimated values of skin mechanical properties from the model agree with those in literature. The extracted model features show good correlation with skin age (viscosity and time constant). The results from our model indicate that the cosmetic polymers enhance the mechanical properties of skin significantly.

CONCLUSIONS

This work will find its applications in designing and testing anti-ageing formulations. This model can be used to filter various combinations of cosmetic formulations by knowing the mechanical response of polymer on skin, thereby accelerating the product development.

Predicting and understanding collagen remodeling in human native heart valves during early development

The hemodynamic functionality of heart valves strongly depends on the distribution of collagen fibers, which are their main load-bearing constituents. It is known that collagen networks remodel in response to mechanical stimuli. Yet, the complex interplay between external load and collagen remodeling is poorly understood. In this study, we adopted a computational approach to simulate collagen remodeling occurring in native fetal and pediatric heart valves. The computational model accounted for several biological phenomena: cellular (re)orientation in response to both mechanical stimuli and topographical cues provided by collagen fibers; collagen deposition and traction forces along the main cellular direction; collagen degradation decreasing with stretch; and cell-mediated collagen prestretch. Importantly, the computational results were well in agreement with previous experimental data for all simulated heart valves. Simulations performed by varying some of the computational parameters suggest that cellular forces and (re)orientation in response to mechanical stimuli may be fundamental mechanisms for the emergence of the circumferential collagen alignment usually observed in native heart valves. On the other hand, the tendency of cells to coalign with collagen fibers is essential to maintain and reinforce that circumferential alignment during development. STATEMENT OF SIGNIFICANCE: The hemodynamic functionality of heart valves is strongly influenced by the alignment of load-bearing collagen fibers. Currently, the mechanisms that are responsible for the development of the circumferential collagen alignment in native heart valves are not fully understood. In the present study, cell-mediated remodeling of native human heart valves during early development was computationally simulated to understand the impact of individual mechanisms on collagen alignment. Our simulations successfully predicted the degree of collagen alignment observed in native fetal and pediatric semilunar valves. The computational results suggest that the circumferential collagen alignment arises from cell traction and cellular (re)orientation in response to mechanical stimuli, and with increasing age is reinforced by the tendency of cells to co-align with pre-existing collagen fibers.

Strain mediated enzymatic degradation of arterial tissue: Insights into the role of the non-collagenous tissue matrix and collagen crimp

Collagen fibre remodelling is a strain dependent process which is stimulated by the degradation of existing collagen. To date, literature has focussed on strain dependent degradation of pure collagen or structurally simple collagenous tissues, often overlooking degradation within more complex, heterogenous soft tissues. The aim of this study is to identify, for the first time, the strain dependent degradation behaviour and mechanical factors influencing collagen degradation in arterial tissue using a combined experimental and numerical approach. To achieve this, structural analysis was carried out using small angle light scattering to determine the fibre level response due to strain induced degradation. Next, strain dependent degradation rates were determined from stress relaxation experiments in the presence of crude and purified collagenase to determine the tissue level degradation response. Finally, a 1D theoretical model was developed, incorporating matrix stiffness and a gradient of collagen fibre crimp to decouple the mechanism behind strain dependent arterial degradation. SALS structural analysis identified a strain mediated degradation response in arterial tissue at the fibre level not dissimilar to that found in literature for pure collagen. Interestingly, two distinctly different strain mediated degradation responses were identified experimentally at the tissue level, not seen in other collagenous tissues. Our model was able to accurately predict these experimental findings, but only once the load bearing matrix, its degradation response and the gradient of collagen fibre crimp across the arterial wall were incorporated. These findings highlight the critical role that the various tissue constituents play in the degradation response of arterial tissue.

STATEMENT OF SIGNIFICANCE

Collagen fibre architecture is the dominant load bearing component of arterial tissue. Remodelling of this architecture is a strain dependent process stimulated by the degradation of existing collagen. Despite this, degradation of arterial tissue and in particular, arterial collagen, is not fully understood or studied. In the current study, we identified for the first time, the strain dependent degradation response of arterial tissue, which has not been observed in other collagenous tissues in literature. We hypothesised that this unique degradation response was due to the complex structure observed in arterial tissue. Based on this hypothesis, we developed a novel numerical model capable of explaining this unique degradation response which may provide critical insights into disease development and aid in the design of interventional medical devices.

Modelling The Combined Effects Of Collagen and Cyclic Strain On Cellular Orientation In Collagenous Tissues

Adherent cells are generally able to reorient in response to cyclic strain. In three-dimensional tissues, however, extracellular collagen can affect this cellular response. In this study, a computational model able to predict the combined effects of mechanical stimuli and collagen on cellular (re)orientation was developed. In particular, a recently proposed computational model (which only accounts for mechanical stimuli) was extended by considering two hypotheses on how collagen influences cellular (re)orientation: collagen contributes to cell alignment by providing topographical cues (contact guidance); or collagen causes a spatial obstruction for cellular reorientation (steric hindrance). In addition, we developed an evolution law to predict cell-induced collagen realignment. The hypotheses were tested by simulating bi- or uniaxially constrained cell-populated collagen gels with different collagen densities, subjected to immediate or delayed uniaxial cyclic strain with varying strain amplitudes. The simulation outcomes are in agreement with previous experimental reports. Taken together, our computational approach is a promising tool to understand and predict the remodeling of collagenous tissues, such as native or tissue-engineered arteries and heart valves.

New algorithm for simulation of proteoglycan loss and collagen degeneration in the knee joint: Data from the osteoarthritis initiative

Osteoarthritis is a harmful joint disease but prediction of disease progression is problematic. Currently, there is only one modeling framework which can be applied to predict the progression of knee osteoarthritis but it only considers degenerative changes in the collagen fibril network. Here, we have developed the framework further by considering all of the major tissue changes (proteoglycan content, fluid flow, and collagen fibril network) occurring in osteoarthritis. While excessive levels of tissue stresses controlled degeneration of the collagen fibril network, excessive levels of tissue strains controlled the decrease in proteoglycan content and the increase in permeability. We created four knee joint models with increasing degrees of complexity based on the depth-wise composition. Models were tested for normal and abnormal, physiologically relevant, loading conditions in the knee. Finally, the predicted depth-wise compositional changes from each model were compared against experimentally observed compositional changes in vitro. The model incorporating the typical depth-wise composition of cartilage produced the best match with experimental observations. Consistent with earlier in vitro experiments, this model simulated the greatest proteoglycan depletion in the superficial and middle zones, while the collagen fibril degeneration was located mostly in the superficial zone. The presented algorithm can be used for predicting simultaneous collagen degeneration and proteoglycan loss during the development of knee osteoarthritis. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:1673-1683, 2018.

A homeostatic-driven turnover remodelling constitutive model for healing in soft tissues

Remodelling of soft biological tissue is characterized by interacting biochemical and biomechanical events, which change the tissue's microstructure, and, consequently, its macroscopic mechanical properties. Remodelling is a well-defined stage of the healing process, and aims at recovering or repairing the injured extracellular matrix. Like other physiological processes, remodelling is thought to be driven by homeostasis, i.e. it tends to re-establish the properties of the uninjured tissue. However, homeostasis may never be reached, such that remodelling may also appear as a continuous pathological transformation of diseased tissues during aneurysm expansion, for example. A simple constitutive model for soft biological tissues that regards remodelling as homeostatic-driven turnover is developed. Specifically, the recoverable effective tissue damage, whose rate is the sum of a mechanical damage rate and a healing rate, serves as a scalar internal thermodynamic variable. In order to integrate the biochemical and biomechanical aspects of remodelling, the healing rate is, on the one hand, driven by mechanical stimuli, but, on the other hand, subjected to simple metabolic constraints. The proposed model is formulated in accordance with continuum damage mechanics within an open-system thermodynamics framework. The numerical implementation in an in-house finite-element code is described, particularized for Ogden hyperelasticity. Numerical examples illustrate the basic constitutive characteristics of the model and demonstrate its potential in representing aspects of remodelling of soft tissues. Simulation results are verified for their plausibility, but also validated against reported experimental data.

In Vitro Bone Cell Models: Impact of Fluid Shear Stress on Bone Formation

This review describes the role of bone cells and their surrounding matrix in maintaining bone strength through the process of bone remodeling. Subsequently, this work focusses on how bone formation is guided by mechanical forces and fluid shear stress in particular. It has been demonstrated that mechanical stimulation is an important regulator of bone metabolism. Shear stress generated by interstitial fluid flow in the lacunar-canalicular network influences maintenance and healing of bone tissue. Fluid flow is primarily caused by compressive loading of bone as a result of physical activity. Changes in loading, e.g., due to extended periods of bed rest or microgravity in space are associated with altered bone remodeling and formation in vivo. In vitro, it has been reported that bone cells respond to fluid shear stress by releasing osteogenic signaling factors, such as nitric oxide, and prostaglandins. This work focusses on the application of in vitro models to study the effects of fluid flow on bone cell signaling, collagen deposition, and matrix mineralization. Particular attention is given to in vitro set-ups, which allow long-term cell culture and the application of low fluid shear stress. In addition, this review explores what mechanisms influence the orientation of collagen fibers, which determine the anisotropic properties of bone. A better understanding of these mechanisms could facilitate the design of improved tissue-engineered bone implants or more effective bone disease models.

Micromechanical Modeling Study of Mechanical Inhibition of Enzymatic Degradation of Collagen Tissues

This study investigates how the collagen fiber structure influences the enzymatic degradation of collagen tissues. We developed a micromechanical model of a fibrous collagen tissue undergoing enzymatic degradation based on two central hypotheses. The collagen fibers are crimped in the undeformed configuration. Enzymatic degradation is an energy activated process and the activation energy is increased by the axial strain energy density of the fiber. We determined the intrinsic degradation rate and characteristic energy for mechanical inhibition from fibril-level degradation experiments and applied the parameters to predict the effect of the crimped fiber structure and fiber properties on the degradation of bovine cornea and pericardium tissues under controlled tension. We then applied the model to examine the effect of the tissue stress state on the rate of tissue degradation and the anisotropic fiber structures that developed from enzymatic degradation.

A computational analysis of cell-mediated compaction and collagen remodeling in tissue-engineered heart valves

One of the most critical problems in heart valve tissue engineering is the progressive development of valvular insufficiency due to leaflet retraction. Understanding the underlying mechanisms of this process is crucial for developing tissue-engineered heart valves (TEHVs) that maintain their functionality in the long term. In the present study, we adopted a computational approach to predict the remodeling process in TEHVs subjected to dynamic pulmonary and aortic pressure conditions, and to assess the risk of valvular insufficiency. In addition, we investigated the importance of the intrinsic cell contractility on the final outcome of the remodeling process. For valves implanted in the aortic position, the model predictions suggest that valvular insufficiency is not likely to occur as the blood pressure is high enough to prevent the development of leaflet retraction. In addition, the collagen network was always predicted to remodel towards a circumferentially aligned network, which is corresponding to the native situation. In contrast, for valves implanted in the pulmonary position, our model predicted that there is a high risk for the development of valvular insufficiency, unless the cell contractility is very low. Conversely, the development of a circumferential collagen network was only predicted at these pressure conditions when cell contractility was high. Overall, these results, therefore, suggest that tissue remodeling at aortic pressure conditions is much more stable and favorable compared to tissue remodeling at pulmonary pressure conditions.