Time Evolution of Deformation in a Human Cartilage Under Cyclic Loading

The future of in-field sports biomechanics: wearables plus modelling compute real-time in vivo tissue loading to prevent and repair musculoskeletal injuries

This paper explores the use of biomechanics in identifying the mechanistic causes of musculoskeletal tissue injury and degeneration. It appraises how biomechanics has been used to develop training programmes aiming to maintain or recover tissue health. Tissue health depends on the functional mechanical environment experienced by tissues during daily and rehabilitation activities. These environments are the result of the interactions between tissue motion, loading, biology, and morphology. Maintaining health of and/or repairing musculoskeletal tissues requires targeting the "ideal" in vivo tissue mechanics (i.e., loading and deformation), which may be enabled by appropriate real-time biofeedback. Recent research shows that biofeedback technologies may increase their quality and effectiveness by integrating a personalised neuromusculoskeletal modelling driven by real-time motion capture and medical imaging. Model personalisation is crucial in obtaining physically and physiologically valid predictions of tissue biomechanics. Model real-time execution is crucial and achieved by code optimisation and artificial intelligence methods. Furthermore, recent work has also shown that laboratory-based motion capture biomechanical measurements and modelling can be performed outside the laboratory with wearable sensors and artificial intelligence. The next stage is to combine these technologies into well-designed easy to use products to guide training to maintain or recover tissue health in the real-world.

Osteochondral junction leakage and cartilage joint lubrication

BACKGROUND AND OBJECTIVES

Previous studies have shown that there is potentially interstitial fluid exchange between cartilage tissue and the subarticular spongiosa region in the case of injury or disease (e.g., osteoarthritis and osteoporosis). Interstitial flow is also required for cartilage lubrication under joint load. A key question then is how cartilage lubrication is modified by increased interstitial fluid leakage across the osteochondral junction. Thus, the purpose of this study is to develop a numerical model to investigate changes in cartilage lubrication with changes in osteochondral junction leakage.

METHODS

The multi-phase coupled model includes domains corresponding to the contact gap, cartilage tissue and subchondral bone plate region (ScBP). Each of these domains are treated as poroelastic systems, with their coupling implemented through mass and pressure continuity. The effects of osteochondral junction leakage on lubrication were investigated with a parametric study on the relative permeability between the ScBP and cartilage tissue.

RESULTS

Significant effects of ScBP permeability were predicted, especially during the early stage of the junction leakage development (early stage of the disease). There is a significant reduction in mixed-mode lubrication duration under the effect of increased junction leakage (the cartilage tissue mixed-mode lubrication duration is about 33% decrease for a relative permeability ratio of 0.1 between ScBP and cartilage tissue, and about 52% decrease under the osteoarthritis condition). In addition, the time for cartilage to reach steady-state consolidation is significantly reduced when ScBP permeability increases (the consolidation time reduces from roughly 2 h to 1.2 h when the relative permeability ratio increases from 0.001 to 0.1, and it reduces to 0.8 h for an advanced osteoarthritis condition). It is predicted that the initial friction coefficient could increase by over 60% when the ScBP permeability is consistent with an advanced osteoarthritis (OA) condition.

CONCLUSION

Increased osteochondral junction leakage induced by joint injury and disease could result in increased cartilage surface wear rates due to more rapid interstitial fluid depressurization within articular cartilage.

A review on design of scaffold for osteoinduction: Toward the unification of independent design variables

Biophysical stimulus quantifies the osteoinductivity of the scaffold concerning the mechanoregulatory mathematical models of scaffold-assisted cellular differentiation. Consider a set of independent structural variables ($) that comprises bulk porosity levels ([Formula: see text]) and a set of morphological features of the micro-structure ([Formula: see text]) associated with scaffolds, i.e., [Formula: see text]. The literature suggests that biophysical stimulus ([Formula: see text]) is a function of independent structural variables ($). Limited understanding of the functional correlation between biophysical stimulus and structural features results in the lack of the desired osteoinductivity in a scaffold. Consequently, it limits their broad applicability to assist bone tissue regeneration for treating critical-sized bone fractures. The literature indicates the existence of multi-dimensional independent design variable space as a probable reason for the general lack of osteoinductivity in scaffolds. For instance, known morphological features are the size, shape, orientation, continuity, and connectivity of the porous regions in the scaffold. It implies that the number of independent variables ([Formula: see text]) is more than two, i.e., [Formula: see text], which interact and influence the magnitude of [Formula: see text] in a unified manner. The efficiency of standard engineering design procedures to analyze the correlation between dependent variable ([Formula: see text]) and independent variables ($) in 3D mutually orthogonal Cartesian coordinate system diminishes proportionally with the increase in the number of independent variables ([Formula: see text]) (Deb in Optimization for engineering design-algorithms and examples, PHI Learning Private Limited, New Delhi, 2012). Therefore, there is an immediate need to devise a framework that has the potential to quantify the micro-structural's morphological features in a unified manner to increase the prospects of scaffold-assisted bone tissue regeneration.

Injury-related cell death and proteoglycan loss in articular cartilage: Numerical model combining necrosis, reactive oxygen species, and inflammatory cytokines

Osteoarthritis (OA) is a common musculoskeletal disease that leads to deterioration of articular cartilage, joint pain, and decreased quality of life. When OA develops after a joint injury, it is designated as post-traumatic OA (PTOA). The etiology of PTOA remains poorly understood, but it is known that proteoglycan (PG) loss, cell dysfunction, and cell death in cartilage are among the first signs of the disease. These processes, influenced by biomechanical and inflammatory stimuli, disturb the normal cell-regulated balance between tissue synthesis and degeneration. Previous computational mechanobiological models have not explicitly incorporated the cell-mediated degradation mechanisms triggered by an injury that eventually can lead to tissue-level compositional changes. Here, we developed a 2-D mechanobiological finite element model to predict necrosis, apoptosis following excessive production of reactive oxygen species (ROS), and inflammatory cytokine (interleukin-1)-driven apoptosis in cartilage explant. The resulting PG loss over 30 days was simulated. Biomechanically triggered PG degeneration, associated with cell necrosis, excessive ROS production, and cell apoptosis, was predicted to be localized near a lesion, while interleukin-1 diffusion-driven PG degeneration was manifested more globally. Interestingly, the model also showed proteolytic activity and PG biosynthesis closer to the levels of healthy tissue when pro-inflammatory cytokines were rapidly inhibited or cleared from the culture medium, leading to partial recovery of PG content. The numerical predictions of cell death and PG loss were supported by previous experimental findings. Furthermore, the simulated ROS and inflammation mechanisms had longer-lasting effects (over 3 days) on the PG content than localized necrosis. The mechanobiological model presented here may serve as a numerical tool for assessing early cartilage degeneration mechanisms and the efficacy of interventions to mitigate PTOA progression.

Computational Modelling for Managing Pathways to Cartilage Failure

Over several decades the perception and therefore description of articular cartilage changed substantially. It has transitioned from being described as a relatively inert tissue with limited repair capacity, to a tissue undergoing continuous maintenance and even adaption, through a range of complex regulatory processes. Even from the narrower lens of biomechanics, the engagement with articular cartilage has changed from it being an interesting, slippery material found in the hostile mechanical environment between opposing long bones, to an intriguing example of mechanobiology in action. The progress revealing this complexity, where physics, chemistry, material science and biology are merging, has been described with increasingly sophisticated computational models. Here we describe how these computational models of cartilage as an integrated system can be combined with the approach of structural reliability analysis. That is, causal, deterministic models placed in the framework of the probabilistic approach of structural reliability analysis could be used to understand, predict, and mitigate the risk of cartilage failure or pathology. At the heart of this approach is seeing cartilage overuse and disease processes as a 'material failure', resulting in failure to perform its function, which is largely mechanical. One can then describe pathways to failure, for example, how homeostatic repair processes can be overwhelmed leading to a compromised tissue. To illustrate this 'pathways to failure' approach, we use the interplay between cartilage consolidation and lubrication to analyse the increase in expected wear rates associated with cartilage defects or meniscectomy.

Understanding the Influence of Local Physical Stimuli on Chondrocyte Behavior

Investigating the mechanobiology of chondrocytes is challenging due to the complex micromechanical environment of cartilage tissue. The innate zonal differences and poroelastic properties of the tissue combined with its heterogeneous composition create spatial- and temporal-dependent cell behavior, which further complicates the investigation. Despite the numerous challenges, understanding the mechanobiology of chondrocytes is crucial for developing strategies for treating cartilage related diseases as chondrocytes are the only cell type within the tissue. The effort to understand chondrocyte behavior under various mechanical stimuli has been ongoing over the last 50 years. Early studies examined global biosynthetic behavior under unidirectional mechanical stimulus. With the technological development in high-speed confocal imaging techniques, recent studies have focused on investigating real-time individual and collective cell responses to multiple / combined modes of mechanical stimuli. Such efforts have led to tremendous advances in understanding the influence of local physical stimuli on chondrocyte behavior. In addition, we highlight the wide variety of experimental techniques, spanning from static to impact loading, and analysis techniques, from biochemical assays to machine learning, that have been utilized to study chondrocyte behavior. Finally, we review the progression of hypotheses about chondrocyte mechanobiology and provide a perspective on the future outlook of chondrocyte mechanobiology.

Intra-Articular Mesenchymal Stem Cell Injection for Knee Osteoarthritis: Mechanisms and Clinical Evidence

Knee osteoarthritis presents higher incidences than other joints, with increased prevalence during aging. It is a progressive process and may eventually lead to disability. Mesenchymal stem cells (MSCs) are expected to repair damaged issues due to trilineage potential, trophic effects, and immunomodulatory properties of MSCs. Intra-articular MSC injection was reported to treat knee osteoarthritis in many studies. This review focuses on several issues of intra-articular MSC injection for knee osteoarthritis, including doses of MSCs applied for injection and the possibility of cartilage regeneration following MSC injection. Intra-articular MSC injection induced hyaline-like cartilage regeneration, which could be seen by arthroscopy in several studies. Additionally, anatomical, biomechanical, and biochemical changes during aging and other causes participate in the development of knee osteoarthritis. Conversely, appropriate intervention based on these anatomical, biomechanical, biochemical, and functional properties and their interactions may postpone the progress of knee OA and facilitate cartilage repair induced by MSC injection. Hence, post-injection rehabilitation programs and related mechanisms are discussed.

Effect of Articular Surface Compression on Cartilage Extracellular Matrix Deformation

Early stage osteoarthritis is characterized by disruption of the superficial zone (SZ) of articular cartilage, including collagen damage and proteoglycan loss, resulting in "mechanical softening" of the extracellular matrix (ECM). The role of the SZ in controlling fluid exudation and imbibition during loading and unloading, respectively, was studied using confined creep compression tests. Bovine osteochondral (OC) plugs were subjected to either a static (88 kPa) or cyclic (0-125 kPa at 1 Hz) compressive stress for five minutes, and the cartilage deformation and recovery were measured during tissue loading and unloading, respectively. During unloading, the articular surface of the cartilage was either loaded with a small 1% tare load (∼1 kPa) applied through a porous load platen (covered), or completely unloaded (uncovered). Then the SZ (∼10%) of the cartilage was removed and the creep tests were repeated. Randomized tests were performed on each OC specimen to assess variability within and between plugs. Static creep strain was always greater than cyclic creep strain except at the beginning of loading (10-20 cycles). Uncovering the articular surface after creep deformation resulted in faster thickness recovery compared to the covered recovery. Removal of the SZ resulted in increased static and cyclic creep strains, as well as an increase in the cyclic peak-to-peak strain envelope. Our results indicate that an intact SZ is essential for normal cartilage mechanical function during joint motion by controlling fluid exudation and imbibition, and concomitantly ECM deformation and recovery, when loaded and unloaded, respectively.

Vibrometry as a noncontact alternative to dynamic and viscoelastic mechanical testing in cartilage

Physiological loading of knee cartilage is highly dynamic and may contribute to the progression of osteoarthritis. Thus, an understanding of cartilage's dynamic mechanical properties is crucial in cartilage research. In this study, vibrometry was used as a fast (2 h), noncontact and novel alternative to the slower (30 h), traditional mechanical and biochemical assays for characterization of cartilage from the condyle, patella, trochlear groove and meniscus. Finite-element models predicted tissue resonant frequencies and bending modes, which strongly correlated with experiments (R2 = 0.93). Vibrometry-based viscoelastic properties significantly correlated with moduli from stress relaxation and creep tests, with correlation strengths reaching up to 0.78. Loss modulus also strongly correlated with glycosoaminoglycan (GAG) content. Dynamic properties measured by vibrometry significantly differed among various knee cartilages, ranging between 6.1 and 56.4 MPa. Interestingly, meniscus viscoelastic properties suggest that contrary to common belief, it may lack shock absorption abilities; instead, condylar hyaline cartilage may be a better shock absorber. These data demonstrate for the first time that vibrometry is a noncontact approach to dynamic mechanical characterization of hyaline and fibrocartilage cartilage with concrete relationships to standard quasi-static mechanical testing and biochemical composition. Thus, with a single tool, vibrometry greatly facilitates meeting multiple regulatory recommendations for mechanical characterization of cartilage replacements.

Reversible changes in the 3D collagen fibril architecture during cyclic loading of healthy and degraded cartilage

Biomechanical changes to the collagen fibrillar architecture in articular cartilage are believed to play a crucial role in enabling normal joint function. However, experimentally there is little quantitative knowledge about the structural response of the Type II collagen fibrils in cartilage to cyclic loading in situ, and the mechanisms that drive the ability of cartilage to withstand long-term repetitive loading. Here we utilize synchrotron small-angle X-ray scattering (SAXS) combined with in-situ cyclic loading of bovine articular cartilage explants to measure the fibrillar response in deep zone articular cartilage, in terms of orientation, fibrillar strain and inter-fibrillar variability in healthy cartilage and cartilage degraded by exposure to the pro-inflammatory cytokine IL-1β. We demonstrate that under repeated cyclic loading the fibrils reversibly change the width of the fibrillar orientation distribution whilst maintaining a largely consistent average direction of orientation. Specifically, the effect on the fibrillar network is a 3-dimensional conical orientation broadening around the normal to the joint surface, inferred by 3D reconstruction of X-ray scattering peak intensity distributions from the 2D pattern. Further, at the intrafibrillar level, this effect is coupled with reversible reduction in fibrillar pre-strain under compression, alongside increase in the variability of fibrillar pre-strain. In IL-1β degraded cartilage, the collagen rearrangement under cyclic loading is disrupted and associated with reduced tissue stiffness. These finding have implications as to how changes in local collagen nanomechanics might drive disease progression or vice versa in conditions such as osteoarthritis and provides a pathway to a mechanistic understanding of such diseases.

Statement of significance

Structural deterioration in biomechanically loaded musculoskeletal organs, e.g., joint osteoarthritis and back pain, are linked to breakdown and changes in their collagen-rich cartilaginous tissue matrix. A critical component enabling cartilage biomechanics is the ultrastructural collagen fibrillar network in cartilage. However, experimental probes of the dynamic structural response of cartilage collagen to biomechanical loads are limited. Here, we use X-ray scattering during cyclic loading (as during walking) on joint tissue to show that cartilage fibrils resist loading by a reversible, three-dimensional orientation broadening and disordering mechanism at the molecular level, and that inflammation reduces this functionality. Our results will help understand how changes to small-scale tissue mechanisms are linked to ageing and osteoarthritic progression, and development of biomaterials for joint replacements.

Hybrid fluorescence-AFM explores articular surface degeneration in early osteoarthritis across length scales

Atomic force microscopy (AFM) has become a powerful tool for the characterization of materials at the nanoscale. Nevertheless, its application to hierarchical biological tissue like cartilage is still limited. One reason is that such samples are usually millimeters in size, while the AFM delivers much more localized information. Here a combination of AFM and fluorescence microscopy is presented where features on a millimeter sized tissue sample are selected by fluorescence microscopy on the micrometer scale and then mapped down to nanometer precision by AFM under native conditions. This served us to show that local changes in the organization of fluorescent stained cells, a marker for early osteoarthritis, correlate with a significant local reduction of the elastic modulus, local thinning of the collagen fibers, and a roughening of the articular surface. This approach is not only relevant for cartilage, but in general for the characterization of native biological tissue from the macro- to the nanoscale. STATEMENT OF SIGNIFICANCE: Different length scales have to be studied to understand the function and dysfunction of hierarchically organized biomaterials or tissues. Here we combine a highly stable AFM with fluorescence microscopy and precisely motorized movement to correlate micro- and nanoscopic properties of articular cartilage on a millimeter sized sample under native conditions. This is necessary for unraveling the relationship between microscale organization of chondrocytes, micrometer scale changes in articular cartilage properties and nanoscale organization of collagen (including D-banding). We anticipate that such studies pave the way for a guided design of hierarchical biomaterials.

A Probabilistic Failure Risk Approach to The Problem of Articular Cartilage Lubrication

BACKGROUND AND OBJECTIVE

The geometrical and mechanical properties that characterise the cartilage contact gap are uncertain and spatially varied. To date the effects of such uncertainties on cartilage lubrication have not been explored. Using a probabilistic approach, the purpose of this study is to numerically investigate the influence of surficial cartilage glycoaminoglycan (GAG) content on joint lubrication behaviour. Gap asperity stiffness and polymer brush border (PBB) thickness are affected by the uncertainty of surficial GAG concentration, and so their correlated effects in maintaining hydrodynamic joint lubrication are investigated.

METHODS

Correlated sampling data are first generated by Monte Carlo simulation. These data are used as inputs for the cartilage contact model, which includes three distinctive features of cartilage tissue (tension-compression nonlinearity, aggrecan dependent permeability and compressive modulus) and fluid flow resistance effects of PBB on cartilage surface. The degree of hydrodynamic lubrication after thirty minutes of constant loading is used as an indicator for assessing the lubrication performance at the contact interface.

RESULTS

The increase of PBB thickness with GAG concentration enhances the hydrodynamic lubrication component in the cartilage contact gap, whereas increasing the asperity stiffness with GAG concentration impairs hydrodynamic lubrication. GAG loss rate increases with the rise of GAG concentration. More aggrecan shedding through the surface could result in a thicker and denser PBB, and therefore enhance the lubrication performance in mixed-mode regime. On the other hand, higher GAG content makes the asperities stiffer, which may impede contact gap closure, and thus encourage gap fluid loss and impair the lubrication performance of cartilage.

CONCLUSION

The lubrication performance of cartilage varies with the physiological conditions of the joint. Since a range of variables are internally related, the outcomes on joint lubrication are difficult to predict. A probabilistic approach accounting for the uncertainties can potentially result in more accurate evaluations of joint lubrication performance.

The investigation of bone fracture healing under intramembranous and endochondral ossification

After trauma, fractured bone starts healing directly through bone union or indirectly through callus formation process. Intramembranous and endochondral ossification are two commonly known mechanisms of indirect healing. The present study investigated the bone fracture healing under intramembranous and endochondral ossification by developing theoretical models in conjunction with performing a series of animal experiments. Using experimentally determined mean bone densities in sheep tibia stabilized by the Locking Compression Plate (LCP) fixation system, the research outcomes showed that intramembranous and endochondral ossification can be described by Hill Function with two unique sets of function parameters in mechanical stimuli mediated fracture healing. Two different thresholds exist within the range of mechanical simulation index which could trigger significant intramembranous and endochondral ossification, with a relatively higher bone formation rate of endochondral ossification than that of intramembranous ossification. Furthermore, the increase of flexibility of the LCP system and the use of titanium LCP could potentially promote uniform bone formation across the fracture gap, ultimately better healing outcomes.

Compressive stress relaxation behavior of articular cartilage and its effects on fluid pressure and solid displacement due to non-Newtonian flow

In this study, we investigate the effects of the power-law index and permeability parameter on the deformation of soft tissue (articular cartilage) which is bathed in the non-Newtonian fluid under stress-relaxation in compression. Ramp displacement is imposed on the surface of hydrated soft tissue. Deformation of the tissue and the fluid pressure is examined for the fast and slow rate of compression. We have employed a linear biphasic mixture theory to develop a mathematical model for compressive stress-relaxation behavior of articular cartilage for non-Newtonian flow. Numerical results indicate that shear-thinning fluids induce less solid deformation and exhibit more fluid pressure as compared to shear-thickening fluids for fast and slow rate of compression. The results also show that linear permeability induces more deformation as compared to strain-dependent nonlinear permeability due to viscoelastic nature of articular cartilage.

The status and challenges of replicating the mechanical properties of connective tissues using additive manufacturing

The ability to fabricate complex structures via precise and heterogeneous deposition of biomaterials makes additive manufacturing (AM) a leading technology in the creation of implants and tissue engineered scaffolds. Connective tissues (CTs) remain attractive targets for manufacturing due to their "simple" tissue compositions that, in theory, are replicable through choice of biomaterial(s) and implant microarchitecture. Nevertheless, characterisation of the mechanical and biological functions of 3D printed constructs with respect to their host tissues is often limited and remains a restriction towards their translation into clinical practice. This review aims to provide an update on the current status of AM to mimic the mechanical properties of CTs, with focus on arterial tissue, articular cartilage and bone, from the perspective of printing platforms, biomaterial properties, and topological design. Furthermore, the grand challenges associated with the AM of CT replacements and their subsequent regulatory requirements are discussed to aid further development of reliable and effective implants.

A probabilistic-based approach for computational simulation of bone fracture healing

BACKGROUND AND OBJECTIVE

It is widely known that bone fracture healing is affected by mechanical factors such as fracture geometry, fixation configuration and post-operative weight bearing loading. However, there are several uncertainties associated with the magnitude of the mechanical factors affecting bone healing as it is challenging to adjust and control them in clinical practice. The current bone fracture healing investigations mainly adopt a deterministic approach for identifying the optimal mechanical conditions for a favourable bone healing outcome. However, a probabilistic approach should be used in the analysis to incorporate such uncertainties for prediction of bone healing.

METHODS

In this study we developed a probabilistic-based computational model to predict the probability of delayed healing or non-union under different fracture treatment mechanical conditions for fractures stabilised by locking plates.

RESULTS

The results show that there is a strong positive linear correlation between the mechanical stimulations (S) in the fracture gap and the magnitude of weight bearing, the bone-plate distance (BPD) and the plate working length (WL), whereas the fracture gap size has a highly negative and nonlinear correlation with S. While the results show that fracture mechanical microenvironment is more sensitive to the uncertainties in WL compared to BPD, the uncertainty associated with the magnitude of WL is very low and can be resulted from implant manufacturing tolerance. However, there is a high uncertainty associated with the magnitude of BPD as it cannot be accurately adjusted during the surgery. The results show that the tissue differentiation at the far cortex of fracture gap is more sensitive to the variation of BPD compared with that at the near cortex. The probability of delayed healing (fibrous tissue formation) at far cortex is increased from 0% to 40% when coefficient of variation (COV) of BPD rises from 0.1 to 0.9 (for average BPD = 2 mm, WL = 65 mm, fracture gap size = 3 mm and Weight bearing = 150 N). Further, both near and far cortex of fracture site are sensitive to the variation in weight bearing loading.

CONCLUSIONS

The developed probabilistic model may lead to useful guidelines that could help orthopaedic surgeons identify how reliable a specific fracture treatment strategy is.

Propagation of microcracks in collagen networks of cartilage under mechanical loads

OBJECTIVE

We recently demonstrated that low-energy mechanical impact to articular cartilage, usually considered non-injurious, can in fact cause microscale cracks (widths <30μm) in the collagen network of visually pristine human cartilage. While research on macro-scale cracks in cartilage and microcracks in bone abounds, how microcracks within cartilage initiate and propagate remains unknown. We quantified the extent to which microcracks initiate and propagate in the collagen network during mechanical loading representative of normal activities.

DESIGN

We tested 76 full-thickness, cylindrical osteochondral plugs. We imaged untreated specimens (pristine phase) via second harmonic generation and assigned specimens to three low-energy impact groups (none, low, high), and thereafter to three cyclic compression groups (none, low, high) which simulate walking. We re-imaged specimens in the post-impact and post-cyclic compression phases to identify and track microcracks.

RESULTS

Microcracks in the network of collagen did not present in untreated controls but did initiate and propagate under mechanical treatments. We found that the length and width of microcracks increased from post-impact to post-cyclic compression in tracked microcracks, but neither depth nor angle presented statistically significant differences.

CONCLUSIONS

The microcracks we initiated under low-energy impact loading increased in length and width during subsequent cyclic compression that simulated walking. The extent of this propagation depended on the combination of impact and cyclic compression. More broadly, the initiation and propagation of microcracks may characterize pathogenesis of osteoarthritis, and may suggest therapeutic targets for future studies.

The Effects of Dynamic Loading on Bone Fracture Healing Under Ilizarov Circular Fixators

Early weight bearing appears to enhance bone fracture healing under Ilizarov circular fixators (ICFs). However, the role of early weight bearing in the healing process remains unclear. This study aims to provide insights into the effects of early weight bearing on healing of bone fractures stabilized with ICFs, with the aid of mathematical modeling. A computational model of fracture site was developed using poro-elastic formulation to simulate the transport of mesenchymal stem cells (MSCs), fibroblasts, chondrocytes, osteoblasts, osteogenic growth factor (OGF), and chondrogenic growth factor (CGF) and MSC differentiation during the early stage of healing, under various combinations of fracture gap sizes (GS), ICF wire pretension forces, and axial loads. 1 h of physiologically relevant cyclic axial loading followed by 23 h of rest in the post-inflammation phase (i.e., callus with granulation tissue) was simulated. The results show that physiologically relevant dynamic loading could significantly enhance cell and growth factor concentrations in the fracture site in a time and spatially dependent manner. 1 h cyclic loading (axial load with amplitude, PA, of 200 N at 1 Hz) increased the content of chondrocytes up to 37% (in all zones of callus), CGF up to 28% (in endosteal and periosteal callus) and OGF up to 50% (in endosteal and cortical callus) by the end of the 24 h period simulated. This suggests that the synergistic effect of dynamic loading-induced advective transport and mechanical stimuli due to early weight bearing is likely to enhance secondary healing. Furthermore, the study suggests that relatively higher PA values or lower ICF wire pretension forces or smaller GS could result in increased chondrocyte and GF content within the callus.

Effects of creep and creep-recovery on ratcheting strain of articular cartilage under cyclic compression

Since the accumulations of ratcheting strain combined with creep deformation, which are produced in normal activities, can accelerate the fatigue damage of cartilage in joint, the creep-ratcheting and creep-recovery-ratcheting behaviors of articular cartilage are experimentally investigated under creep-fatigue loads. The effect of pre-creep on ratcheting behavior of cartilage was probed firstly and it is found that the initial ratcheting strain of cartilage presents the larger value (30% and 35%) due to its pre-creep deformation in spite of the short pre-creep time applied. With the increasing pre-creep time the ratcheting strain of sample increases while the ratcheting strain rate decreases. The effects of pre-creep and recovery on ratcheting behavior of cartilage were also investigated and it is noted that the strain of cartilage increases fast at first, decreases partly and then changes periodically with cyclic stress. The ratcheting strain evolutions of different layers are not coincident for cartilage sample after the pre-creep or pre-creep-recovery and the ratcheting strains of different layers with pre-creep are larger than them with pre-creep-recovery. Finally the creep-ratcheting strain of cartilage with different peak-holding time during cyclic compression was studied and it is found that the creep-ratcheting strain with peak-holding time is significantly greater than that without peak-holding time. The creep-ratcheting strain increases with increase of peak-holding time or stress amplitude, while it reduces with rising of stress rate. The creep-ratcheting behavior of cartilage is dependent on its depth. These findings point out that the accumulated deformations, including creep deformation and ratcheting deformation, can accelerate the cartilage's damage.

Mechanical and biological performance of printed alginate/methylcellulose/halloysite nanotube/polyvinylidene fluoride bio-scaffolds

Use of artificial cartilage due to its poor regenerative characteristics is a challenging issue in the field of tissue engineering. In this regard, three-dimensional printing (3D) technique because of its perfect structural control is one of the best methods for producing biological scaffolds. Proper biomaterials for cartilage repairs with good mechanical and biological properties and the high ability for 3D printing are limited. In this paper, a novel biomaterial consisting of Alginate (AL), Methylcellulose (MC), Halloysite Nanotube (HNT), and Polyvinylidene Fluoride (PVDF) was printed and characterized for cartilage scaffold applications. Calcium chloride (CaCl₂) was used as a crosslinker for biomaterial after printing. Scanning Electron Microscopy (SEM), Energy-Dispersive X-Ray Spectroscopy (EDX), X-Ray Diffraction (XRD), Fourier-Transform Infrared Spectroscopy (FT-IR), Differential Scanning Calorimetry (DSC), tensile and compressive tests, chondrocytes seeding, cells staining, and MTT assay were carried out in the present work. The results show that in constant concentrations of AL, MC, and PVDF (40 mg/ml AL, 30 mg/ml MC, and 1% PVDF) when concentration of HNT increased from 20 mg/ml (S2) to 40 mg/ml (S14) tensile strength increased from 164 up to 381 kPa and compressive stress increased from 426 up to 648 kPa. According to spectroscopy and calorimetry results, Biomaterial shows an amorphous structure with good miscibility and a high percentage of water in its structure. PVDF reduces mechanical properties by 7% while increases cell viability by 8.75%. Histological studies and MTT assay results showed a high improvement in the percentage of living cells at the first 4 days of cell cultivation.

Functional MRI can detect changes in intratissue strains in a full thickness and critical sized ovine cartilage defect model

Functional imaging of tissue biomechanics can reveal subtle changes in local softening and stiffening associated with disease or repair, but noninvasive and nondestructive methods to acquire intratissue measures in well-defined animal models are largely lacking. We utilized displacement encoded MRI to measure changes in cartilage deformation following creation of a critical-sized defect in the medial femoral condyle of ovine (sheep) knees, a common in situ and large animal model of tissue damage and repair. We prioritized visualization of local, site-specific variation and changes in displacements and strains following defect placement by measuring spatial maps of intratissue deformation. Custom data smoothing algorithms were developed to minimize propagation of noise in the acquired MRI phase data toward calculated displacement or strain, and to improve strain measures in high aspect ratio tissue regions. Strain magnitudes in the femoral, but not tibial, cartilage dramatically increased in load-bearing and contact regions especially near the defect locations, with an average 6.7% ± 6.3%, 13.4% ± 10.0%, and 10.0% ± 4.9% increase in first and second principal strains, and shear strain, respectively. Strain heterogeneity reflected the complexity of the in situ mechanical environment within the joint, with multiple tissue contacts defining the deformation behavior. This study demonstrates the utility of displacement encoded MRI to detect increased deformation patterns and strain following disruption to the cartilage structure in a clinically-relevant, large animal defect model. It also defines imaging biomarkers based on biomechanical measures, in particular shear strain, that are potentially most sensitive to evaluate damage and repair, and that may additionally translate to humans in future studies.

A computational algorithm to simulate disorganization of collagen network in injured articular cartilage

Cartilage defects are a known risk factor for osteoarthritis. Estimation of structural changes in these defects could help us to identify high risk defects and thus to identify patients that are susceptible for the onset and progression of osteoarthritis. Here, we present an algorithm combined with computational modeling to simulate the disorganization of collagen fibril network in injured cartilage. Several potential triggers for collagen disorganization were tested in the algorithm following the assumption that disorganization is dependent on the mechanical stimulus of the tissue. We found that tensile tissue stimulus alone was unable to preserve collagen architecture in intact cartilage as collagen network reoriented throughout the cartilage thickness. However, when collagen reorientation was based on both tensile tissue stimulus and tensile collagen fibril strains or stresses, the collagen network architecture was preserved in intact cartilage. Using the same approach, substantial collagen reorientation was predicted locally near the cartilage defect and particularly at the cartilage-bone interface. The developed algorithm was able to predict similar structural findings reported in the literature that are associated with experimentally observed remodeling in articular cartilage. The proposed algorithm, if further validated, could help to predict structural changes in articular cartilage following post-traumatic injury potentially advancing to impaired cartilage function.

Osteoarthritis year in review 2016: mechanics

Inappropriate biomechanics, namely wear-and-tear, has been long believed to be a main cause of osteoarthritis (OA). However, this view is now being re-evaluated, especially when examined alongside mechanobiology and new biomechanical studies. These are multiscale experimental and computational studies focussing on cell- and tissue-level mechanobiology through to organ- and whole-body-level biomechanics, which focuses on the biomechanical and biochemical environment of the joint tissues. This review examined papers from April 2015 to April 2016, with a focus on multiscale experimental and computational biomechanical studies of OA. Assessing the onset or progression of OA at organ- and whole-body-levels, gait analysis, medical imaging and neuromusculoskeletal modelling revealed the extent to which tissue damage changes the view of inappropriate biomechanics. Traditional gait analyses studies reported that conservative treatments can alter joint biomechanics, thereby improving pain and function experienced by those with OA. Results of animal models of OA were consistent with these human studies, showing interactions among bone, cartilage and meniscus biomechanics and the onset and/or progression OA. Going down size scales, experimental and computational studies probed the nanosize biomechanics of molecules, cells and extracellular matrix, and demonstrated how the interactions between biomechanics and morphology affect cartilage dynamic poroelastic behaviour and pathways to OA. Finally, integration of multiscale experimental data and computational models were proposed to predict cartilage extracellular matrix remodelling and the development of OA. Summarising, experimental and computational methods provided a nuanced biomechanical understanding of the sub-cellular, cellular, tissue, organ and whole-body mechanisms involved in OA.

An analysis of changes in in vivo cartilage thickness of the healthy ankle following dynamic activity

Abnormal cartilage loading after injury is believed to be an important factor leading to post-traumatic ankle osteoarthritis. Due to the viscoelastic behavior of cartilage, it is possible to measure localized cartilage strains from changes in thickness following dynamic activities. However, there are limited data characterizing in vivo cartilage mechanics under physiological loading conditions in the healthy ankle. Therefore, the objective of this study was to directly measure in vivo cartilage strains in the healthy ankle joint in response to a dynamic hopping exercise. Ten healthy subjects with no history of ankle injury underwent magnetic resonance imaging before and after a single-leg hopping exercise. Bony and articular cartilage surfaces were created from these images using solid modeling software. Pre-exercise and post-exercise models were then registered to each other, and site-specific cartilage strains (defined as the normalized changes in cartilage thickness) were calculated at grid points spanning the articular surfaces. The effects of both location and exercise on strain were tested using a two-way repeated measures analysis of variance. We did not detect any significant interaction effect between location and exercise for either tibial or talar cartilage. However, hopping resulted in significant decreases in tibial (p<0.05) and talar (p<0.05) cartilage thicknesses, corresponding to strains of 3% and 2%, respectively. Additionally, pre-exercise cartilage thickness varied significantly by location in the talus (p<0.05), but not in the tibia. These strain data may provide important baseline information for future studies investigating altered biomechanics in those at high risk for the development of post-traumatic ankle osteoarthritis.

Osteoarthritis year in review 2015: mechanics

Motivated by the conceptual framework of multi-scale biomechanics, this narrative review highlights recent major advances with a focus on gait and joint kinematics, then tissue-level mechanics, cell mechanics and mechanotransduction, matrix mechanics, and finally the nanoscale mechanics of matrix macromolecules. A literature review was conducted from January 2014 to April 2015 using PubMed to identify major developments in mechanics related to osteoarthritis (OA). Studies of knee adduction, flexion, rotation, and contact mechanics have extended our understanding of medial compartment loading. In turn, advances in measurement methodologies have shown how injuries to both the meniscus and ligaments, together, can alter joint kinematics. At the tissue scale, novel findings have emerged regarding the mechanics of the meniscus as well as cartilage superficial zone. Moving to the cell level, poroelastic and poro-viscoelastic mechanisms underlying chondrocyte deformation have been reported, along with the response to osmotic stress. Further developments have emerged on the role of calcium signaling in chondrocyte mechanobiology, including exciting findings on the function of mechanically activated cation channels newly found to be expressed in chondrocytes. Finally, AFM-based nano-rheology systems have enabled studies of thin murine tissues and brush layers of matrix molecules over a wide range of loading rates including high rates corresponding to impact injury. With OA acknowledged to be a disease of the joint as an organ, understanding mechanical behavior at each length scale helps to elucidate the connections between cell biology, matrix biochemistry and tissue structure/function that may play a role in the pathomechanics of OA.

Predicting Knee Osteoarthritis

Treatment options for osteoarthritis (OA) beyond pain relief or total knee replacement are very limited. Because of this, attention has shifted to identifying which factors increase the risk of OA in vulnerable populations in order to be able to give recommendations to delay disease onset or to slow disease progression. The gold standard is then to use principles of risk management, first to provide subject-specific estimates of risk and then to find ways of reducing that risk. Population studies of OA risk based on statistical associations do not provide such individually tailored information. Here we argue that mechanistic models of cartilage tissue maintenance and damage coupled to statistical models incorporating model uncertainty, united within the framework of structural reliability analysis, provide an avenue for bridging the disciplines of epidemiology, cell biology, genetics and biomechanics. Such models promise subject-specific OA risk assessment and personalized strategies for mitigating or even avoiding OA. We illustrate the proposed approach with a simple model of cartilage extracellular matrix synthesis and loss regulated by daily physical activity.