Modeling sample/patient-specific structural and diffusional responses of cartilage using DT-MRI

Quantifying fat zonation in liver lobules: an integrated multiscale in silico model combining disturbed microperfusion and fat metabolism via a continuum biomechanical bi-scale, tri-phasic approach

Metabolic zonation refers to the spatial separation of metabolic functions along the sinusoidal axes of the liver. This phenomenon forms the foundation for adjusting hepatic metabolism to physiological requirements in health and disease (e.g., metabolic dysfunction-associated steatotic liver disease/MASLD). Zonated metabolic functions are influenced by zonal morphological abnormalities in the liver, such as periportal fibrosis and pericentral steatosis. We aim to analyze the interplay between microperfusion, oxygen gradient, fat metabolism and resulting zonated fat accumulation in a liver lobule. Therefore we developed a continuum biomechanical, tri-phasic, bi-scale, and multicomponent in silico model, which allows to numerically simulate coupled perfusion-function-growth interactions two-dimensionally in liver lobules. The developed homogenized model has the following specifications: (i) thermodynamically consistent, (ii) tri-phase model (tissue, fat, blood), (iii) penta-substances (glycogen, glucose, lactate, FFA, and oxygen), and (iv) bi-scale approach (lobule, cell). Our presented in silico model accounts for the mutual coupling between spatial and time-dependent liver perfusion, metabolic pathways and fat accumulation. The model thus allows the prediction of fat development in the liver lobule, depending on perfusion, oxygen and plasma concentration of free fatty acids (FFA), oxidative processes, the synthesis and the secretion of triglycerides (TGs). The use of a bi-scale approach allows in addition to focus on scale bridging processes. Thus, we will investigate how changes at the cellular scale affect perfusion at the lobular scale and vice versa. This allows to predict the zonation of fat distribution (periportal or pericentral) depending on initial conditions, as well as external and internal boundary value conditions.

Microstructural and mechanical insight into atherosclerotic plaques: an ex vivo DTI study to better assess plaque vulnerability

Non-invasive microstructural characterisation has the potential to determine the stability, or lack thereof, of atherosclerotic plaques and ultimately aid in better assessing plaques' risk to rupture. If linked with mechanical characterisation using a clinically relevant imaging technique, mechanically sensitive rupture risk indicators could be possible. This study aims to provide this link-between a clinically relevant imaging technique and mechanical characterisation within human atherosclerotic plaques. Ex vivo diffusion tensor imaging, mechanical testing, and histological analysis were carried out on human carotid atherosclerotic plaques. DTI-derived tractography was found to yield significant mechanical insight into the mechanical properties of more stable and more vulnerable microstructures. Coupled with insights from digital image correlation and histology, specific failure characteristics of different microstructural arrangements furthered this finding. More circumferentially uniform microstructures failed at higher stresses and strains when compared to samples which had multiple microstructures, like those seen in a plaque cap. The novel findings in this study motivate diagnostic measures which use non-invasive characterisation of the underlying microstructure of plaques to determine their vulnerability to rupture.

Magnetic Resonance Imaging-based biomechanical simulation of cartilage: A systematic review

MRI-based mathematical and computational modeling studies can contribute to a better understanding of the mechanisms governing cartilage's mechanical performance and cartilage disease. In addition, distinct modeling of cartilage is needed to optimize artificial cartilage production. These studies have opened up the prospect of further deepening our understanding of cartilage function. Furthermore, these studies reveal the initiation of an engineering-level approach to how cartilage disease affects material properties and cartilage function. Aimed at researchers in the field of MRI-based cartilage simulation, research articles pertinent to MRI-based cartilage modeling were identified, reviewed, and summarized systematically. Various MRI applications for cartilage modeling are highlighted, and the limitations of different constitutive models used are addressed. In addition, the clinical application of simulations and studied diseases are discussed. The paper's quality, based on the developed questionnaire, was assessed, and out of 79 reviewed papers, 34 papers were determined as high-quality. Due to the lack of the best constitutive models for various clinical conditions, researchers may consider the effect of constitutive material models on the cartilage disease simulation. In the future, research groups may incorporate various aspects of machine learning into constitutive models and MRI data extraction to further refine the study methodology. Moreover, researchers should strive for further reproducibility and rigorous model validation and verification, such as gait analysis.

Anisotropic Diffusivity Tensor in Articular Cartilage: Effective Medium Approach

Due to the avascular nature of articular cartilage, molecular transport occurs via interstitial fluid flow as well as via diffusion. Diffusion in cartilage has been studied experimentally, but no mathematical models have been developed to interpret the experimental results and the observed isotropy or anisotropy in the different cartilage zones. Here, we propose a model for the determination of the diffusivity tensor of uncharged macromolecules in articular cartilage, accounting for the inhomogeneity and anisotropy arising from fiber arrangement, volumetric fraction, and radius. We study a representative element of volume (REV) comprising a fiber surrounded by fluid-saturated proteoglycan matrix. The REV permeability tensor is evaluated using a previously developed model, while the REV diffusivity tensor is obtained by incorporating the hydrodynamic effect and the steric effect of the fiber-reinforced matrix. Both effects are represented by anisotropic second-order tensors. The overall diffusivity tensor is obtained as the averaging integral of the REV diffusivity, weighted by the probability distribution of fiber orientation. The model's predictions of the trend of the magnitude of the diffusivity of spheroidal macromolecules as a function of molecular radius agree with published experimental results. For large linear macromolecules, the model underestimates the diffusivity magnitude (i.e., the equivalent isotropic diffusivity). The model correctly predicts the anisotropic behavior for linear macromolecules, although it underestimates the numerical value of the diffusivity anisotropy ratio of large linear macromolecules in the superficial zone, and overestimates it in the deep zone. In summary, this model constitutes a first step toward understanding the relation between diffusivity and permeability in articular cartilage.

On incorporating osmotic prestretch/prestress in image-driven finite element simulations of cartilage

Medical imaging performed in vivo captures geometries under Donnan osmotic loading, even when the articulating joint is otherwise mechanically unloaded. Hence patient-specific finite element (FE) models constructed from such medical images of cartilage represent osmotically induced prestretched/prestressed states. When applying classical modeling approaches to patient-specific simulations of cartilage a theoretical inconsistency arises: the in-vivo imaged geometry (used to construct the model) is not an unloaded, stress-free reference configuration. Furthermore when fitting nonlinear constitutive models that include osmotic swelling (to obtain material parameters), if one assumes that experimental data-generated from osmotically loaded cartilage-begin from a stress-free reference configuration the fitted stress-stretch relationship (parameters) obtained will actually describe a different behavior. In this study we: (1) establish a practical computational method to include osmotically induced prestretch/prestress in image-driven simulations of cartilage; and (2) investigate the influence of considering the prestretched/prestressed state both when fitting fiber-reinforced, biphasic constitutive models of cartilage that include osmotic swelling and when simulating cartilage responses. Our results highlight the importance of determining the prestretched/prestressed state within cartilage induced by osmotic loading in the imaged configuration prior to solving boundary value problems of interest. With our new constitutive model and modeling methods, we aim to improve the fidelity of FE-based, patient-specific biomechanical simulations of joints and cartilage. Improved simulations can provide medical researchers with new information often unavailable in a clinical setting, information that may contribute to better insight into the pathophysiology of cartilage diseases.

Hydrogels for Cartilage Regeneration, from Polysaccharides to Hybrids

The aims of this paper are: (1) to review the current state of the art in the field of cartilage substitution and regeneration; (2) to examine the patented biomaterials being used in preclinical and clinical stages; (3) to explore the potential of polymeric hydrogels for these applications and the reasons that hinder their clinical success. The studies about hydrogels used as potential biomaterials selected for this review are divided into the two major trends in tissue engineering: (1) the use of cell-free biomaterials; and (2) the use of cell seeded biomaterials. Preparation techniques and resulting hydrogel properties are also reviewed. More recent proposals, based on the combination of different polymers and the hybridization process to improve the properties of these materials, are also reviewed. The combination of elements such as scaffolds (cellular solids), matrices (hydrogel-based), growth factors and mechanical stimuli is needed to optimize properties of the required materials in order to facilitate tissue formation, cartilage regeneration and final clinical application. Polymer combinations and hybrids are the most promising materials for this application. Hybrid scaffolds may maximize cell growth and local tissue integration by forming cartilage-like tissue with biomimetic features.

Reactive Constrained Mixtures for Modeling the Solid Matrix of Biological Tissues

This article illustrates our approach for modeling the solid matrix of biological tissues using reactive constrained mixtures. Several examples are presented to highlight the potential benefits of this approach, showing that seemingly disparate fields of mechanics and chemical kinetics are actually closely interrelated and may be elegantly expressed in a unified framework. Thus, constrained mixture models recover classical theories for fibrous materials with bundles oriented in different directions or having different reference configurations, that produce characteristic fiber recruitment patterns under loading. Reactions that exchange mass among various constituents of a mixture may be used to describe tissue growth and remodeling, which may also alter the material's anisotropy. Similarly, reactions that describe the breaking and reforming of bonds may be used to model free energy dissipation in a viscoelastic material. Therefore, this framework is particularly well suited for modeling biological tissues.

Patient-specific fracture risk assessment of vertebrae: A multiscale approach coupling X-ray physics and continuum micromechanics

While in clinical settings, bone mineral density measured by computed tomography (CT) remains the key indicator for bone fracture risk, there is an ongoing quest for more engineering mechanics-based approaches for safety analyses of the skeleton. This calls for determination of suitable material properties from respective CT data, where the traditional approach consists of regression analyses between attenuation-related grey values and mechanical properties. We here present a physics-oriented approach, considering that elasticity and strength of bone tissue originate from the material microstructure and the mechanical properties of its elementary components. Firstly, we reconstruct the linear relation between the clinically accessible grey values making up a CT, and the X-ray attenuation coefficients quantifying the intensity losses from which the image is actually reconstructed. Therefore, we combine X-ray attenuation averaging at different length scales and over different tissues, with recently identified 'universal' composition characteristics of the latter. This gives access to both the normally non-disclosed X-ray energy employed in the CT-device and to in vivo patient-specific and location-specific bone composition variables, such as voxel-specific mass density, as well as collagen and mineral contents. The latter feed an experimentally validated multiscale elastoplastic model based on the hierarchical organization of bone. Corresponding elasticity maps across the organ enter a finite element simulation of a typical load case, and the resulting stress states are increased in a proportional fashion, so as to check the safety against ultimate material failure. In the young patient investigated, even normal physiological loading is probable to already imply plastic events associated with the hydrated mineral crystals in the bone ultrastructure, while the safety factor against failure is still as high as five. Copyright © 2016 John Wiley & Sons, Ltd.

Investigation of the mechanical behavior of kangaroo humeral head cartilage tissue by a porohyperelastic model based on the strain-rate-dependent permeability

Solid-interstitial fluid interaction, which depends on tissue permeability, is significant to the strain-rate-dependent mechanical behavior of humeral head (shoulder) cartilage. Due to anatomical and biomechanical similarities to that of the human shoulder, kangaroos present a suitable animal model. Therefore, indentation experiments were conducted on kangaroo shoulder cartilage tissues from low (10(-4)/s) to moderately high (10(-2)/s) strain-rates. A porohyperelastic model was developed based on the experimental characterization; and a permeability function that takes into account the effect of strain-rate on permeability (strain-rate-dependent permeability) was introduced into the model to investigate the effect of rate-dependent fluid flow on tissue response. The prediction of the model with the strain-rate-dependent permeability was compared with those of the models using constant permeability and strain-dependent permeability. Compared to the model with constant permeability, the models with strain-dependent and strain-rate-dependent permeability were able to better capture the experimental variation at all strain-rates (p < 0.05). Significant differences were not identified between models with strain-dependent and strain-rate-dependent permeability at strain-rate of 5 × 10(-3)/s (p = 0.179). However, at strain-rate of 10(-2)/s, the model with strain-rate-dependent permeability was significantly better at capturing the experimental results (p < 0.005). The findings thus revealed the significance of rate-dependent fluid flow on tissue behavior at large strain-rates, which provides insights into the mechanical deformation mechanisms of cartilage tissues.

Techniques and applications of in vivo diffusion imaging of articular cartilage

Early in the process of osteoarthritis (OA) the composition (water, proteoglycan [PG], and collagen) and structure of articular cartilage is altered leading to changes in its mechanical properties. A technique that can assess the composition and structure of the cartilage in vivo can provide insight in the mechanical integrity of articular cartilage and become a powerful tool for the early diagnosis of OA. Diffusion tensor imaging (DTI) has been proposed as a biomarker for cartilage composition and structure. DTI is sensitive to the PG content through the mean diffusivity and to the collagen architecture through the fractional anisotropy. However, the acquisition of DTI of articular cartilage in vivo is challenging due to the short T2 of articular cartilage (∼40 ms at 3 Tesla) and the high resolution needed (0.5-0.7 mm in plane) to depict the cartilage anatomy. We describe the pulse sequences used for in vivo DTI of articular cartilage and discus general strategies for protocol optimization. We provide a comprehensive review of measurements of DTI of articular cartilage from ex vivo validation experiments to its recent clinical applications.

Stress distribution and consolidation in cartilage constituents is influenced by cyclic loading and osteoarthritic degeneration

The understanding of load support mechanisms in cartilage has evolved with computational models that better mimic the tissue ultrastructure. Fibril-reinforced poroelastic models can reproduce cartilage behaviour in a variety of test conditions and can be used to model tissue anisotropy as well as assess stress and pressure partitioning to the tissue constituents. The goal of this study was to examine the stress distribution in the fibrillar and non-fibrillar solid phase and pressure in the fluid phase of cartilage in axisymmetric models of a healthy and osteoarthritic hip joint. Material properties, based on values from the literature, were assigned to the fibrillar and poroelastic components of cartilage and cancellous and subchondral compact bone regions. A cyclic load representing walking was applied for 25 cycles. Contact stresses in the fibrillar and non-fibrillar solid phase supported less than 1% of the contact force and increased only minimally with load cycles. Simulated proteoglycan depletion increased stresses in the radial and tangential collagen fibrils, whereas fibrillation of the tangential fibrils resulted in increased compressive stress in the non-fibrillar component and tensile stress in the radial fibrils. However neither had an effect on fluid pressure. Subchondral sclerosis was found to have the largest effect, resulting in increased fluid pressure, non-fibrillar compressive stress, tangential fibril stress and greater cartilage consolidation. Subchondral bone stiffening may play an important role in the degenerative cascade and may adversely affect tissue repair and regeneration treatments.