Multiscale mechanics of articular cartilage: potentials and challenges of coupling musculoskeletal, joint, and microscale computational models

Femtosecond Laser Maskless Optical Projection Lithography of Cartilage PCM Inspired 3D Protein Matrix to Chondrocyte Phenotype

Hydrogels containing chondrocytes have exhibited excellent potential in regenerating hyaline cartilage. However, chondrocytes are vulnerable to dedifferentiation during in vitro culture, leading to fibrosis and mechanical degradation of newly formed cartilage. It is proposed to modulate cartilage formation via the developed chondrocyte pericellular matrix (PCM) -like scaffolds for the first time, in which the S, M, and L-sized scaffolds are fabricated by femtosecond laser maskless optical projection lithography (FL-MOPL) of bovine serum albumin-glyceryl methacrylate hydrogel. Chondrocytes on the M PCM-like scaffold can maintain round morphology and synthesize extracellular matrix (ECM) to induce regeneration of hyaline cartilage microtissues by geometrical restriction. A series of M PCM-like scaffolds is fabricated with different stiffness and those with a high Young's modulus are more effective in maintaining the chondrocyte phenotype. The proposed PCM-like scaffolds are effective in modulating cartilage formation influenced by pore size, depth, and stiffness, which will pave the way for a better understanding of the geometric cues of mechanotransduction interactions in regulating cell fate and open up new avenues for tissue engineering.

Osteoarthritis year in review 2023: Biomechanics

Biomechanics plays a significant yet complex role in osteoarthritis (OA) onset and progression. Identifying alterations in biomechanical factors and their complex interactions is critical for gaining new insights into OA pathophysiology and identification of clearly defined and modifiable mechanical treatment targets. This review synthesized biomechanics studies from March 2022 to April 2023, from which three themes relating to human gait emerged: (1) new insights into the pathogenesis of OA using computational modeling and machine learning, (2) technology-enhanced biomechanical interventions for OA, and (3) out-of-lab biomechanical assessments of OA. We further highlighted future-focused areas which may continue to advance the field of biomechanics in OA, with a particular emphasis on exploiting technology to understand and treat biomechanical mechanisms of OA outside the laboratory. The breadth of studies included in this review highlights the complex role of biomechanics in OA and showcase numerous innovative and outstanding contributions to the field. Exciting cross-disciplinary efforts integrating computational modeling, mobile sensors, and machine learning methods show great promise for streamlining in vivo multi-scale biomechanics workflows and are expected to underpin future breakthroughs in the understanding and treatment of biomechanics in OA.

Contributions of cell behavior to geometric order in embryonic cartilage

During early development, cartilage provides shape and stability to the embryo while serving as a precursor for the skeleton. Correct formation of embryonic cartilage is hence essential for healthy development. In vertebrate cranial cartilage, it has been observed that a flat and laterally extended macroscopic geometry is linked to regular microscopic structure consisting of tightly packed, short, transversal clonar columns. However, it remains an ongoing challenge to identify how individual cells coordinate to successfully shape the tissue, and more precisely which mechanical interactions and cell behaviors contribute to the generation and maintenance of this columnar cartilage geometry during embryogenesis. Here, we apply a three-dimensional cell-based computational model to investigate mechanical principles contributing to column formation. The model accounts for clonal expansion, anisotropic proliferation and the geometrical arrangement of progenitor cells in space. We confirm that oriented cell divisions and repulsive mechanical interactions between cells are key drivers of column formation. In addition, the model suggests that column formation benefits from the spatial gaps created by the extracellular matrix in the initial configuration, and that column maintenance is facilitated by sequential proliferative phases. Our model thus correctly predicts the dependence of local order on division orientation and tissue thickness. The present study presents the first cell-based simulations of cell mechanics during cranial cartilage formation and we anticipate that it will be useful in future studies on the formation and growth of other cartilage geometries.

Computational assessment on the impact of collagen fiber orientation in cartilages on healthy and arthritic knee kinetics/kinematics

BACKGROUND

The inhomogeneous distribution of collagen fiber in cartilage can substantially influence the knee kinematics. This becomes vital for understanding the mechanical response of soft tissues, and cartilage deterioration including osteoarthritis (OA). Though the conventional computational models consider geometrical heterogeneity along with fiber reinforcements in the cartilage model as material heterogeneity, the influence of fiber orientation on knee kinetics and kinematics is not fully explored. This work examines how the collagen fiber orientation in the cartilage affects the healthy (intact knee) and arthritic knee response over multiple gait activities like running and walking.

METHODS

A 3D finite element knee joint model is used to compute the articular cartilage response during the gait cycle. A fiber-reinforced porous hyper elastic (FRPHE) material is used to model the soft tissue. A split-line pattern is used to implement the fiber orientation in femoral and tibial cartilage. Four distinct intact cartilage models and three OA models are simulated to assess the impact of the orientation of collagen fibers in a depth wise direction. The cartilage models with fibers oriented in parallel, perpendicular, and inclined to the articular surface are investigated for multiple knee kinematics and kinetics.

FINDINGS

The comparison of models with fiber orientation parallel to articulating surface for walking and running gait has the highest elastic stress and fluid pressure compared with inclined and perpendicular fiber-oriented models. Also, the maximum contact pressure is observed to be higher in the case of intact models during the walking cycle than for OA models. In contrast, the maximum contact pressure is higher during running in OA models than in intact models. Additionally, parallel-oriented models produce higher maximum stresses and fluid pressure for walking and running gait than proximal-distal-oriented models. Interestingly, during the walking cycle, the maximum contact pressure with intact models is approximately three times higher than on OA models. In contrast, the OA models exhibit higher contact pressure during the running cycle.

INTERPRETATION

Overall, the study indicates that collagen orientation is crucial for tissue responsiveness. This investigation provides insights into the development of tailored implants.

Osteoarthritis year in review 2021: mechanics

Osteoarthritis (OA) has a complex, heterogeneous and only partly understood etiology. There is a definite role of joint cartilage pathomechanics in originating and progressing of the disease. Although it is still not identified precisely enough to design or select targeted treatments, the progress of this year's research demonstrates that this goal became much closer. On multiple scales - tissue, joint and whole body - an increasing number of studies were done, with impressive results. (1) Technology based instrument innovations, especially when combined with machine learning models, have broadened the applicability of biomechanics. (2) Combinations with imaging make biomechanics much more precise & personalized. (3) The combination of Musculoskeletal & Finite Element Models yield valid personalized cartilage loads. (4) Mechanical outcomes are becoming increasingly meaningful to inform and evaluate treatments, including predictive power from biomechanical models. Since most recent advancements in the field of biomechanics in OA are at the level of a proof op principle, future research should not only continue on this successful path of innovation, but also aim to develop clinical workflows that would facilitate including precision biomechanics in large scale studies. Eventually this will yield clinical tools for decision making and a rationale for new therapies in OA.

Targeted mesenchymal stem cell therapy equipped with a cell-tissue nanomatchmaker attenuates osteoarthritis progression

Mesenchymal stem cells (MSCs) are at the forefront of research for a wide range of diseases, including osteoarthritis (OA). Despite having attracted the attention of orthopedists, current MSC therapy techniques are limited by poor MSC implantation in tissue defects and lack of lateral tissue integration, which has restricted the efficacy of cell therapy to alleviate OA symptoms only. Here, we developed targeted MSC therapy for OA cartilage using a cell-tissue matchmaking nanoconstruct (C-TMN). C-TMN, as an MSC vehicle, consists of a central iron oxide nanoparticle armed with two types of antibodies, one directed at the MSC surface and the other against articular cartilage. We treated rat OA articular cartilage with intra-articular injections of C-TMN with and without exogenous MSCs. We observed substantial improvements in both symptomatic and radiographic OA caused by C-TMN, which was independent of exogenous MSCs. This new approach could predict a promising future for OA management.

OA foundations - experimental models of osteoarthritis

Osteoarthritis (OA) is increasingly recognised as a disease of diverse phenotypes with variable clinical presentation, progression, and response to therapeutic intervention. This same diversity is readily apparent in the many animal models of OA. However, model selection, study design, and interpretation of resultant findings, are not routinely done in the context of the target human (or veterinary) patient OA sub-population or phenotype. This review discusses the selection and use of animal models of OA in discovery and therapeutic-development research. Beyond evaluation of the different animal models on offer, this review suggests focussing the approach to OA-animal model selection on study objective(s), alignment of available models with OA-patient sub-types, and the resources available to achieve valid and translatable results. How this approach impacts model selection is discussed and an experimental design checklist for selecting the optimal model(s) is proposed. This approach should act as a guide to new researchers and a reminder to those already in the field, as to issues that need to be considered before embarking on in vivo pre-clinical research. The ultimate purpose of using an OA animal model is to provide the best possible evidence if, how, when and where a molecule, pathway, cell or process is important in clinical disease. By definition this requires both model and study outcomes to align with and be predictive of outcomes in patients. Keeping this at the forefront of research using pre-clinical OA models, will go a long way to improving the quality of evidence and its translational value.

A Multidisciplinary Journey towards Bone Tissue Engineering

Millions of patients suffer yearly from bone fractures and disorders such as osteoporosis or cancer, which constitute the most common causes of severe long-term pain and physical disabilities. The intrinsic capacity of bone to repair the damaged bone allows normal healing of most small bone injuries. However, larger bone defects or more complex diseases require additional stimulation to fully heal. In this context, the traditional routes to address bone disorders present several associated drawbacks concerning their efficacy and cost-effectiveness. Thus, alternative therapies become necessary to overcome these limitations. In recent decades, bone tissue engineering has emerged as a promising interdisciplinary strategy to mimic environments specifically designed to facilitate bone tissue regeneration. Approaches developed to date aim at three essential factors: osteoconductive scaffolds, osteoinduction through growth factors, and cells with osteogenic capability. This review addresses the biological basis of bone and its remodeling process, providing an overview of the bone tissue engineering strategies developed to date and describing the mechanisms that underlie cell-biomaterial interactions.

A backward pre-stressing algorithm for efficient finite element implementation of in vivo material and geometrical parameters into fibril-reinforced mixture models of articular cartilage

Classical continuum mechanics has been widely used for implementation of the material models of articular cartilage (AC) mainly with the aid of the finite element (FE) method, which, in many cases, considers the stress-free configuration as the initial configuration. On the contrary, the AC experimental tests typically begin with the pre-stressed state of both material and geometrical properties. Indeed, imposing the initial pre-stress onto AC models with the in vivo values as the initial state would result in nonphysiologically expansion of the FE mesh due to the soft nature of AC. This change in the model configuration can also affect the material behavior kinematically in the mixture models of cartilage due to the intrinsic compressibility of the tissue. Although several different fixed-point backward algorithms, as the most straightforward pre-stressing methods, have already been developed to incorporate these initial conditions into FE models iteratively, such methods focused merely on the geometrical parameters, and they omitted the material variations of the anisotropic mixture models of AC. To address this issue, we propose an efficient algorithm generalizing the backward schemes to restore stress-free conditions by optimizing both the involving variables, and we hypothesize that it can affect the results considerably. To this end, a comparative simulation was implemented on an advanced and validated multiphasic model by the new and conventional algorithms. The results are in support of the hypothesis, as in our illustrative general AC model, the material parameters experienced a maximum error of 16% comparing to the initial in vivo data when the older algorithm was employed, and it led to a maximum variation of 44% in the recorded stresses comparing to the results of the new method. We conclude that our methodology enhanced the model fidelity, and it is applicable in most of the existing FE solvers for future mixture studies with accurate stress distributions.

Estimation of joint contact pressure in the index finger using a hybrid finite element musculoskeletal approach

The knowledge of local stress distribution in hand joints is crucial to understand injuries and osteoarthritis occurrence. However, determining cartilage contact stresses remains a challenge, requiring numerical models including both accurate anatomical components and realistic tendon force actuation. Contact forces in finger joints have frequently been calculated but little data is available on joint contact pressures. This study aimed to develop and assess a hybrid biomechanical model of the index finger to estimate in-vivo joint contact pressure during a static maximal strength pinch grip task. A finite element model including bones, cartilage, tendons, and ligaments was developed, with tendon force transmission based on a tendon-pulley system. This model was driven by realistic tendon forces estimated from a musculoskeletal model and motion capture data for six subjects. The hybrid model outputs agreed well with the experimental measurement of fingertip forces and literature data on the physiological distribution of tendon forces through the index finger. Mean contact pressures were 6.9 ± 2.7 MPa, 6.2 ± 1.0 MPa and 7.2 ± 1.3 MPa for distal, proximal interphalangeal and metacarpophalangeal joints, respectively. Two subjects had higher mean contact pressure in the distal joint than in the other two joints, suggesting a mechanical cause for the prevalence of osteoarthritis in the index distal joint. The inter-subject variability in joint contact pressure could be explained by different neuromuscular strategies employed for the task. This first application of an effective hybrid model to the index finger is promising for estimating hand joint stresses under daily grip tasks and simulating surgical procedures.

Use of Computational Modeling to Study Joint Degeneration: A Review

Osteoarthritis (OA), a degenerative joint disease, is the most common chronic condition of the joints, which cannot be prevented effectively. Computational modeling of joint degradation allows to estimate the patient-specific progression of OA, which can aid clinicians to estimate the most suitable time window for surgical intervention in osteoarthritic patients. This paper gives an overview of the different approaches used to model different aspects of joint degeneration, thereby focusing mostly on the knee joint. The paper starts by discussing how OA affects the different components of the joint and how these are accounted for in the models. Subsequently, it discusses the different modeling approaches that can be used to answer questions related to OA etiology, progression and treatment. These models are ordered based on their underlying assumptions and technologies: musculoskeletal models, Finite Element models, (gene) regulatory models, multiscale models and data-driven models (artificial intelligence/machine learning). Finally, it is concluded that in the future, efforts should be made to integrate the different modeling techniques into a more robust computational framework that should not only be efficient to predict OA progression but also easily allow a patient’s individualized risk assessment as screening tool for use in clinical practice.

Proteoglycan degradation mimics static compression by altering the natural gradients in fibrillar organisation in cartilage

Structural and associated biomechanical gradients within biological tissues are important for tissue functionality and preventing damaging interfacial stress concentrations. Articular cartilage possesses an inhomogeneous structure throughout its thickness, driving the associated variation in the biomechanical strain profile within the tissue under physiological compressive loading. However, little is known experimentally about the nanostructural mechanical role of the collagen fibrils and how this varies with depth. Utilising a high-brilliance synchrotron X-ray source, we have measured the depth-wise nanostructural parameters of the collagen network in terms of the periodic fibrillar banding (D-period) and associated parameters. We show that there is a depth dependent variation in D-period reflecting the pre-strain and concurrent with changes in the level of intrafibrillar order. Further, prolonged static compression leads to fibrillar changes mirroring those caused by removal of extrafibrillar proteoglycans (as may occur in aging or disease). We suggest that fibrillar D-period is a sensitive indicator of localised changes to the mechanical environment at the nanoscale in soft connective tissues.

Statement of Significance

Collagen plays a significant role in both the structural and mechanical integrity of articular cartilage, allowing the tissue to withstand highly repetitive loading. However, the fibrillar mechanics of the collagen network in cartilage are not clear. Here we find that cartilage has a spatial gradient in the nanostructural collagen fibril pre-strain, with an increase in the fibrillar pre-strain with depth. Further, the fibrillar gradient changes similarly under compression when compared to an enzymatically degraded tissue which mimics age-related changes. Given that the fibrils potentially have a finite capacity to mechanically respond and alter their configuration, these findings are significant in understanding how collagen may alter in structure and gradient in diseased cartilage, and in informing the design of cartilage replacements.

Numerical Study on Electromechanics in Cartilage Tissue with Respect to Its Electrical Properties

Hyaline cartilage undergoes many substantial age-related physiochemical and biomechanical changes that reduce its ability to overcome the effects of mechanical stress and injury. In quest of therapeutic options, magnetic stimulation and electrical stimulation (ES) have been proposed for improving tissue engineering approaches for the repair of articular cartilage. The aim of this study is to summarize in silico investigations involving induced electrical properties of cartilage tissue due to various biophysical stimuli along their respective mathematical descriptions. Based on these, a preliminary numerical study involving electromechanical transduction in bovine cartilage tissue has been carried out using an open source finite element computational software. The simulation results have been compared to experimental results from the literature. This study serves as a basis for further in silico studies to better understand the behavior of hyaline cartilage tissue due to ES and to find an optimal stimulation protocol for the cartilage regeneration. Moreover, it provides an overview of the basic models along with mathematical description and scope for future research regarding electrical behavior of the cartilage tissue using open source software.

Simulation of surface strain in tibiofemoral cartilage during walking for the prediction of collagen fiber orientation

The collagen fibers in the superficial layer of tibiofemoral articular cartilage exhibit distinct patterns in orientation revealed by split lines. In this study, we introduce a simulation framework to predict cartilage surface loading during walking to investigate if split line orientations correspond with principal strain directions in the cartilage surface. The two-step framework uses a multibody musculoskeletal model to predict tibiofemoral kinematics which are then imposed on a deformable surface model to predict surface strains. The deformable surface model uses absolute nodal coordinate formulation (ANCF) shell elements to represent the articular surface and a system of spring-dampers and internal pressure to represent the underlying cartilage. Simulations were performed to predict surface strains due to osmotic pressure, loading induced by walking, and the combination of both loading due to pressure and walking. Time-averaged magnitude-weighted first principal strain directions agreed well with split line maps from the literature for both the osmotic pressure and combined cases. This result suggests there is indeed a connection between collagen fiber orientation and mechanical loading, and indicates the importance of accounting for the pre-strain in the cartilage surface due to osmotic pressure.

Osteoarthritis as a disease of the cartilage pericellular matrix

Osteoarthritis is a painful joint disease characterized by progressive degeneration of the articular cartilage as well as associated changes to the subchondral bone, synovium, and surrounding joint tissues. While the effects of osteoarthritis on the cartilage extracellular matrix (ECM) have been well recognized, it is now becoming apparent that in many cases, the onset of the disease may be initially reflected in the matrix region immediately surrounding the chondrocytes, termed the pericellular matrix (PCM). Growing evidence suggests that the PCM - which along with the enclosed chondrocytes are termed the "chondron" - acts as a critical transducer or "filter" of biochemical and biomechanical signals for the chondrocyte, serving to help regulate the homeostatic balance of chondrocyte metabolic activity in response to environmental signals. Indeed, it appears that alterations in PCM properties and cell-matrix interactions, secondary to genetic, epigenetic, metabolic, or biomechanical stimuli, could in fact serve as initiating or progressive factors for osteoarthritis. Here, we discuss recent advances in the understanding of the role of the PCM, with an emphasis on the reciprocity of changes that occur in this matrix region with disease, as well as how alterations in PCM properties could serve as a driver of ECM-based diseases such as osteoarthritis. Further study of the structure, function, and composition of the PCM in normal and diseased conditions may provide new insights into the understanding of the pathogenesis of osteoarthritis, and presumably new therapeutic approaches for this disease.

How changes in interconnectivity affect the bulk properties of articular cartilage: a fibre network study

The remarkable compressive strength of articular cartilage arises from the mechanical interactions between the tension-resisting collagen fibrils and swelling proteoglycan proteins within the tissue. These interactions are facilitated by a significant level of interconnectivity between neighbouring collagen fibrils within the extracellular matrix. A reduction in interconnectivity is suspected to occur during the early stages of osteoarthritic degeneration. However, the relative contribution of these interconnections towards the bulk mechanical properties of articular cartilage has remained an open question. In this study, we present a simple 2D fibre network model which explicitly represents the microstructure of articular cartilage as collection of discrete nodes and linear springs. The transverse stiffness and swelling properties of this fibre network are studied, and a semi-analytic relationship which relates these two macroscopic properties via microscopic interconnectivity is derived. By comparing this derived expression to previously published experimental data, we show that although a reduction in network interconnectivity accounts for some of the observed changes in the mechanical properties of articular cartilage as degeneration occurs, a decrease in matrix interconnectivity alone do not provide a full account of this process.

Comparison between in vitro and in vivo cartilage overloading studies based on a systematic literature review

Methodological differences between in vitro and in vivo studies on cartilage overloading complicate the comparison of outcomes. The rationale of the current review was to (i) identify consistencies and inconsistencies between in vitro and in vivo studies on mechanically-induced structural damage in articular cartilage, such that variables worth interesting to further explore using either one of these approaches can be identified; and (ii) suggest how the methodologies of both approaches may be adjusted to facilitate easier comparison and therewith stimulate translation of results between in vivo and in vitro studies. This study is anticipated to enhance our understanding of the development of osteoarthritis, and to reduce the number of in vivo studies. Generally, results of in vitro and in vivo studies are not contradicting. Both show subchondral bone damage and intact cartilage above a threshold value of impact energy. At lower loading rates, excessive loads may cause cartilage fissuring, decreased cell viability, collagen network de-structuring, decreased GAG content, an overall damage increase over time, and low ability to recover. This encourages further improvement of in vitro systems, to replace, reduce, and/or refine in vivo studies. However, differences in experimental set up and analyses complicate comparison of results. Ways to bridge the gap include (i) bringing in vitro set-ups closer to in vivo, for example, by aligning loading protocols and overlapping experimental timeframes; (ii) synchronizing analytical methods; and (iii) using computational models to translate conclusions from in vitro results to the in vivo environment and vice versa. © 2018 The Authors. Journal of Orthopaedic Research® Published by Wiley Periodicals, Inc. J Orthop Res 9999:1-11, 2018.

The effect of constitutive representations and structural constituents of ligaments on knee joint mechanics

Ligaments provide stability to the human knee joint and play an essential role in restraining motion during daily activities. Compression-tension nonlinearity is a well-known characteristic of ligaments. Moreover, simpler material representations without this feature might give reasonable results because ligaments are primarily in tension during loading. However, the biomechanical role of different constitutive representations and their fibril-reinforced poroelastic properties is unknown. A numerical knee model which considers geometric and material nonlinearities of meniscus and cartilages was applied. Five different constitutive models for the ligaments (spring, elastic, hyperelastic, porohyperelastic, and fibril-reinforced porohyperelastic (FRPHE)) were implemented. Knee joint forces for the models with elastic, hyperelastic and porohyperelastic properties showed similar behavior throughout the stance, while the model with FRPHE properties exhibited lower joint forces during the last 50% of the stance phase. The model with ligaments as springs produced the lowest joint forces at this same stance phase. The results also showed that the fibril network contributed substantially to the knee joint forces, while the nonfibrillar matrix and fluid had small effects. Our results indicate that simpler material models of ligaments with similar properties in compression and tension can be used when the loading is directed primarily along the ligament axis in tension.

Multiphoton microscope measurement-based biphasic multiscale analyses of knee joint articular cartilage and chondrocyte by using visco-anisotropic hyperelastic finite element method and smoothed particle hydrodynamics method

The articular cartilage of a knee joint has a variety of functions including dispersing stress and absorbing shock in the tissue and lubricating the surface region of cartilage. The metabolic activity of chondrocytes under the cyclic mechanical stimulations regenerates the morphology and function of tissues. Hence, the stress evaluation of the chondrocyte is a vital subject to assess the regeneration cycle in the normal walking condition and predict the injury occurrence in the accidents. Further, the threshold determination of stress for the chondrocytes activation is valuable for development of regenerative bioreactor of articular cartilage. In this study, in both macroscale and microscale analyses, the dynamic explicit finite element (FE) method was used for the solid phase and the smoothed particle hydrodynamics (SPH) method was used for the fluid phase. In the homogenization procedure, the representative volume element for the microscale finite element model was derived by using the multiphoton microscope measured 3D structure comprising 3 different layers: surface, middle, and deep layers. The layers had different anisotropic structural and rigidity characteristics because of the collagen fiber orientation. In both macroscale and microscale FE analyses, the visco-anisotropic hyperelastic constitutive law was used. Material properties were identified by experimentally determined stress-strain relationships of 3 layers. With respect to the macroscale and microscale SPH models for non-Newtonian viscous fluid, the previous observation results of interstitial fluid and proteoglycan were used to perform parameter identifications. Biphasic multiscale FE and SPH analyses were conducted under normal walking conditions. Therefore, the hydrostatic and shear stresses occurring in the chondrocytes caused by the compressive load and shear viscous flow were evaluated. These stresses will be used to design an ex-vivo bioreactor to regenerate the damaged articular cartilage, where chondrocytes are seeded in the culture chamber. To know the stress occurred on and in the chondrocytes is vitally important not only to understand the normal metabolic activity of the chondrocyte but also to develop a bioreactor of articular cartilage regeneration as the knee joint disease treatment.

The Virtual Physiological Human: Ten Years After

Biomedical research and clinical practice are struggling to cope with the growing complexity that the progress of health care involves. The most challenging diseases, those with the largest socioeconomic impact (cardiovascular conditions; musculoskeletal conditions; cancer; metabolic, immunity, and neurodegenerative conditions), are all characterized by a complex genotype-phenotype interaction and by a "systemic" nature that poses a challenge to the traditional reductionist approach. In 2005 a small group of researchers discussed how the vision of computational physiology promoted by the Physiome Project could be translated into clinical practice and formally proposed the term Virtual Physiological Human. Our knowledge about these diseases is fragmentary, as it is associated with molecular and cellular processes on the one hand and with tissue and organ phenotype changes (related to clinical symptoms of disease conditions) on the other. The problem could be solved if we could capture all these fragments of knowledge into predictive models and then compose them into hypermodels that help us tame the complexity that such systemic behavior involves. In 2005 this was simply not possible-the necessary methods and technologies were not available. Now, 10 years later, it seems the right time to reflect on the original vision, the results achieved so far, and what remains to be done.

The potential for intercellular mechanical interaction: simulations of single chondrocyte versus anatomically based distribution

Computational studies of chondrocyte mechanics, and cell mechanics in general, have typically been performed using single cell models embedded in an extracellular matrix construct. The assumption of a single cell microstructural model may not capture intercellular interactions or accurately reflect the macroscale mechanics of cartilage when higher cell concentrations are considered, as may be the case in many instances. Hence, the goal of this study was to compare cell-level response of single and eleven cell biphasic finite element models, where the latter provided an anatomically based cellular distribution representative of the actual number of cells for a commonly used [Formula: see text] edge cubic representative volume in the middle zone of cartilage. Single cell representations incorporated a centered single cell model and eleven location-corrected single cell models, the latter to delineate the role of cell placement in the representative volume element. A stress relaxation test at 10% compressive strain was adopted for all simulations. During transient response, volume- averaged chondrocyte mechanics demonstrated marked differences (up to 60% and typically greater than 10%) for the centered single versus the eleven cell models, yet steady-state loading was similar. Cell location played a marked role, due to inhomogeneity of the displacement and fluid pressure fields at the macroscopic scale. When the single cell representation was corrected for cell location, the transient response was consistent, while steady-state differences on the order of 1-4% were realized, which may be attributed to intercellular mechanical interactions. Anatomical representations of the superficial and deep zones, where cells reside in close proximity, may exhibit greater intercellular interactions, but these have yet to be explored.

Systems approaches in osteoarthritis: Identifying routes to novel diagnostic and therapeutic strategies

Systems orientated research offers the possibility of identifying novel therapeutic targets and relevant diagnostic markers for complex diseases such as osteoarthritis. This review demonstrates that the osteoarthritis research community has been slow to incorporate systems orientated approaches into research studies, although a number of key studies reveal novel insights into the regulatory mechanisms that contribute both to joint tissue homeostasis and its dysfunction. The review introduces both top-down and bottom-up approaches employed in the study of osteoarthritis. A holistic and multiscale approach, where clinical measurements may predict dysregulation and progression of joint degeneration, should be a key objective in future research. The review concludes with suggestions for further research and emerging trends not least of which is the coupled development of diagnostic tests and therapeutics as part of a concerted effort by the osteoarthritis research community to meet clinical needs. © 2017 The Authors. Journal of Orthopaedic Research Published by Wiley Periodicals, Inc. on behalf of Orthopaedic Research Society. J Orthop Res 35:1573-1588, 2017.

In Vivo Multiscale and Spatially-Dependent Biomechanics Reveals Differential Strain Transfer Hierarchy in Skeletal Muscle

Biological tissues have a complex hierarchical architecture that spans organ to subcellular scales and comprises interconnected biophysical and biochemical machinery. Mechanotransduction, gene regulation, gene protection, and structure-function relationships in tissues depend on how force and strain are modulated from macro to micro scales, and vice versa. Traditionally, computational and experimental techniques have been used in common model systems (e.g., embryos) and simple strain measures were applied. But the hierarchical transfer of mechanical parameters like strain in mammalian systems is largely unexplored in vivo. Here, we experimentally probed complex strain transfer processes in mammalian skeletal muscle tissue over multiple biological scales using complementary in vivo ultrasound and optical imaging approaches. An iterative hyperelastic warping technique quantified the spatially-dependent strain distributions in tissue, matrix, and subcellular (nuclear) structures, and revealed a surprising increase in strain magnitude and heterogeneity in active muscle as the spatial scale also increased. The multiscale strain heterogeneity indicates tight regulation of mechanical signals to the nuclei of individual cells in active muscle, and an emergent behavior appearing at larger (e.g. tissue) scales characterized by dramatically increased strain complexity.

American Society of Biomechanics Clinical Biomechanics Award 2015: MRI assessments of cartilage mechanics, morphology and composition following reconstruction of the anterior cruciate ligament

BACKGROUND

The pathogenesis of osteoarthritis following anterior cruciate ligament (ACL) reconstruction is currently unknown. The study purpose was to leverage recent advances in quantitative and dynamic MRI to test the hypothesis that abnormal joint mechanics within four years of reconstruction is accompanied by evidence of early compositional changes in cartilage.

METHODS

Static MR imaging was performed bilaterally on eleven subjects with an ACL reconstruction (1-4years post-surgery) and on twelve healthy subjects to obtain tibial cartilage thickness maps. Quantitative imaging (mcDESPOT) was performed unilaterally on all subjects to assess the fraction of bound water in the tibial plateau cartilage. Finally, volumetric dynamic imaging was performed to assess cartilage contact patterns during an active knee flexion-extension task. A repeated-measures ANOVA was used to test for the effects of surgical reconstruction and location on cartilage thickness, bound water fractions, and contact across the medial and lateral tibia plateaus.

FINDINGS

No significant differences in cartilage thickness were found between groups. However, there was a significant reduction in the fraction of water bound by proteoglycan in the ACL reconstructed knees, most notably along the anterior portion of the medial plateau and the weight-bearing lateral plateau. During movement, reconstructed knees exhibited greater contact along the medial spine in the medial plateau and along the posterior aspect of the lateral plateau, when compared with their healthy contralateral knees and healthy controls.

INTERPRETATION

This study provides evidence that abnormal mechanics in anterior cruciate ligament reconstructed knees are present coincidently with early biomarkers of cartilage degeneration.

Modelling approaches for evaluating multiscale tendon mechanics

Tendon exhibits anisotropic, inhomogeneous and viscoelastic mechanical properties that are determined by its complicated hierarchical structure and varying amounts/organization of different tissue constituents. Although extensive research has been conducted to use modelling approaches to interpret tendon structure-function relationships in combination with experimental data, many issues remain unclear (i.e. the role of minor components such as decorin, aggrecan and elastin), and the integration of mechanical analysis across different length scales has not been well applied to explore stress or strain transfer from macro- to microscale. This review outlines mathematical and computational models that have been used to understand tendon mechanics at different scales of the hierarchical organization. Model representations at the molecular, fibril and tissue levels are discussed, including formulations that follow phenomenological and microstructural approaches (which include evaluations of crimp, helical structure and the interaction between collagen fibrils and proteoglycans). Multiscale modelling approaches incorporating tendon features are suggested to be an advantageous methodology to understand further the physiological mechanical response of tendon and corresponding adaptation of properties owing to unique in vivo loading environments.

Open Knee: Open Source Modeling and Simulation in Knee Biomechanics

Virtual representations of the knee joint can provide clinicians, scientists, and engineers the tools to explore mechanical functions of the knee and its tissue structures in health and disease. Modeling and simulation approaches such as finite element analysis also provide the possibility to understand the influence of surgical procedures and implants on joint stresses and tissue deformations. A large number of knee joint models are described in the biomechanics literature. However, freely accessible, customizable, and easy-to-use models are scarce. Availability of such models can accelerate clinical translation of simulations, where labor-intensive reproduction of model development steps can be avoided. Interested parties can immediately utilize readily available models for scientific discovery and clinical care. Motivated by this gap, this study aims to describe an open source and freely available finite element representation of the tibiofemoral joint, namely Open Knee, which includes the detailed anatomical representation of the joint's major tissue structures and their nonlinear mechanical properties and interactions. Three use cases illustrate customization potential of the model, its predictive capacity, and its scientific and clinical utility: prediction of joint movements during passive flexion, examining the role of meniscectomy on contact mechanics and joint movements, and understanding anterior cruciate ligament mechanics. A summary of scientific and clinically directed studies conducted by other investigators are also provided. The utilization of this open source model by groups other than its developers emphasizes the premise of model sharing as an accelerator of simulation-based medicine. Finally, the imminent need to develop next-generation knee models is noted. These are anticipated to incorporate individualized anatomy and tissue properties supported by specimen-specific joint mechanics data for evaluation, all acquired in vitro from varying age groups and pathological states.

Shape of chondrocytes within articular cartilage affects the solid but not the fluid microenvironment under unconfined compression

Metabolic activity of the chondrocytes in articular cartilage is strongly related to their zone-specific shape and the composition and mechanical properties of their surrounding extracellular matrix (ECM). However the mechanisms by which cell shape influences the response of the ECM microenvironment to mechanical loading is yet to be elucidated. This relationship was studied using a biphasic multiscale finite element model of different shaped chondrocytes in the superficial and deep zones of the ECM during unconfined stress relaxation. For chondrocytes in the superficial zone, increasing the cell's initial aspect ratio (length/height) increased the deformation and solid stresses of the chondrocyte and pericellular matrix (PCM) during the loading phase; for chondrocytes in the deep zone the effect of the cell shape on the solid microenvironment was time and variable dependent. However, for superficial and deep zone chondrocytes the cell shape did not affect the fluid pressure and fluid shear stress. These results suggest that mechanotransduction of chondrocytes in articular cartilage may be regulated through the solid phase rather than the fluid phase, and that high stresses and deformations in the solid microenvironment in the superficial zone may be essential for the zone-specific biosynthetic activity of the chondrocyte. The biphasic multiscale computational analysis suggests that maintaining the cell shape is critical for regulating the microenvironment and metabolic activity of the chondrocyte in tissue engineering constructs.

STATEMENT OF SIGNIFICANCE

We investigated the effect of chondrocyte shape on the cellular microenvironment using a biphasic multiscale finite element analysis. Our study showed that cell shapes affects the solid but not the fluid microenvironment of the chondrocyte, and that maintaining the cell shape is critical for regulating the microenvironment and metabolic activity of the chondrocyte in native cartilage and tissue engineering constructs. As far as we know, this is the first study on the mechanotransduction mechanisms by which cell shape influences the response of the microenvironment to mechanical loading. This study is important for understanding cell mechanobiology, not only for regulation of cell phenotype in tissue engineered constructs but, as important, for understanding changes in normal chondrocyte function after post-traumatic injury and in the initiation and progression of osteoarthritis.

Discrete Element Framework for Modelling Extracellular Matrix, Deformable Cells and Subcellular Components

This paper presents a framework for modelling biological tissues based on discrete particles. Cell components (e.g. cell membranes, cell cytoskeleton, cell nucleus) and extracellular matrix (e.g. collagen) are represented using collections of particles. Simple particle to particle interaction laws are used to simulate and control complex physical interaction types (e.g. cell-cell adhesion via cadherins, integrin basement membrane attachment, cytoskeletal mechanical properties). Particles may be given the capacity to change their properties and behaviours in response to changes in the cellular microenvironment (e.g., in response to cell-cell signalling or mechanical loadings). Each particle is in effect an ‘agent’, meaning that the agent can sense local environmental information and respond according to pre-determined or stochastic events. The behaviour of the proposed framework is exemplified through several biological problems of ongoing interest. These examples illustrate how the modelling framework allows enormous flexibility for representing the mechanical behaviour of different tissues, and we argue this is a more intuitive approach than perhaps offered by traditional continuum methods. Because of this flexibility, we believe the discrete modelling framework provides an avenue for biologists and bioengineers to explore the behaviour of tissue systems in a computational laboratory.

Modelling is an important tool in understanding the behaviour of biological tissues. In this paper we advocate a new modelling framework in which cells and tissues are represented by a collection of particles with associated properties. The particles interact with each other and can change their behaviour in response to changes in their environment. We demonstrate how the propose framework can be used to represent the mechanical behaviour of different tissues with much greater flexibility as compared to traditional continuum based methods.

Altered swelling and ion fluxes in articular cartilage as a biomarker in osteoarthritis and joint immobilization: a computational analysis

In healthy cartilage, mechano-electrochemical phenomena act together to maintain tissue homeostasis. Osteoarthritis (OA) and degenerative diseases disrupt this biological equilibrium by causing structural deterioration and subsequent dysfunction of the tissue. Swelling and ion flux alteration as well as abnormal ion distribution are proposed as primary indicators of tissue degradation. In this paper, we present an extension of a previous three-dimensional computational model of the cartilage behaviour developed by the authors to simulate the contribution of the main tissue components in its behaviour. The model considers the mechano-electrochemical events as concurrent phenomena in a three-dimensional environment. This model has been extended here to include the effect of repulsion of negative charges attached to proteoglycans. Moreover, we have studied the fluctuation of these charges owning to proteoglycan variations in healthy and pathological articular cartilage. In this sense, standard patterns of healthy and degraded tissue behaviour can be obtained which could be a helpful diagnostic tool. By introducing measured properties of unhealthy cartilage into the computational model, the severity of tissue degeneration can be predicted avoiding complex tissue extraction and subsequent in vitro analysis. In this work, the model has been applied to monitor and analyse cartilage behaviour at different stages of OA and in both short (four, six and eight weeks) and long-term (11 weeks) fully immobilized joints. Simulation results showed marked differences in the corresponding swelling phenomena, in outgoing cation fluxes and in cation distributions. Furthermore, long-term immobilized patients display similar swelling as well as fluxes and distribution of cations to patients in the early stages of OA, thus, preventive treatments are highly recommended to avoid tissue deterioration.

Multiscale cartilage biomechanics: technical challenges in realizing a high-throughput modelling and simulation workflow

Understanding the mechanical environment of articular cartilage and chondrocytes is of the utmost importance in evaluating tissue damage which is often related to failure of the fibre architecture and mechanical injury to the cells. This knowledge also has significant implications for understanding the mechanobiological response in healthy and diseased cartilage and can drive the development of intervention strategies, ranging from the design of tissue-engineered constructs to the establishment of rehabilitation protocols. Spanning multiple spatial scales, a wide range of biomechanical factors dictate this mechanical environment. Computational modelling and simulation provide descriptive and predictive tools to identify multiscale interactions, and can lead towards a greater comprehension of healthy and diseased cartilage function, possibly in an individualized manner. Cartilage and chondrocyte mechanics can be examined in silico, through post-processing or feed-forward approaches. First, joint-tissue level simulations, typically using the finite-element method, solve boundary value problems representing the joint articulation and underlying tissue, which can differentiate the role of compartmental joint loading in cartilage contact mechanics and macroscale cartilage field mechanics. Subsequently, tissue-cell scale simulations, driven by the macroscale cartilage mechanical field information, can predict chondrocyte deformation metrics along with the mechanics of the surrounding pericellular and extracellular matrices. A high-throughput modelling and simulation framework is necessary to develop models representative of regional and population-wide variations in cartilage and chondrocyte anatomy and mechanical properties, and to conduct large-scale analysis accommodating a multitude of loading scenarios. However, realization of such a framework is a daunting task, with technical difficulties hindering the processes of model development, scale coupling, simulation and interpretation of the results. This study aims to summarize various strategies to address the technical challenges of post-processing-based simulations of cartilage and chondrocyte mechanics with the ultimate goal of establishing the foundations of a high-throughput multiscale analysis framework. At the joint-tissue scale, rapid development of regional models of articular contact is possible by automating the process of generating parametric representations of cartilage boundaries and depth-dependent zonal delineation with associated constitutive relationships. At the tissue-cell scale, models descriptive of multicellular and fibrillar architecture of cartilage zones can also be generated in an automated fashion. Through post-processing, scripts can extract biphasic mechanical metrics at a desired point in the cartilage to assign loading and boundary conditions to models at the lower spatial scale of cells. Cell deformation metrics can be extracted from simulation results to provide a simplified description of individual chondrocyte responses. Simulations at the tissue-cell scale can be parallelized owing to the loosely coupled nature of the feed-forward approach. Verification studies illustrated the necessity of a second-order data passing scheme between scales and evaluated the role that the microscale representative volume size plays in appropriately predicting the mechanical response of the chondrocytes. The tools summarized in this study collectively provide a framework for high-throughput exploration of cartilage biomechanics, which includes minimally supervised model generation, and prediction of multiscale biomechanical metrics across a range of spatial scales, from joint regions and cartilage zones, down to that of the chondrocytes.

A multi-scale finite element model for investigation of chondrocyte mechanics in normal and medial meniscectomy human knee joint during walking

Mechanical signals experienced by chondrocytes (articular cartilage cells) modulate cell synthesis and cartilage health. Multi-scale modeling can be used to study how forces are transferred from joint surfaces through tissues to chondrocytes. Therefore, estimation of chondrocyte behavior during certain physical activities, such as walking, could provide information about how cells respond to normal and abnormal loading in joints. In this study, a 3D multi-scale model was developed for evaluating chondrocyte and surrounding peri- and extracellular matrix responses during gait loading within healthy and medial meniscectomy knee joints. The knee joint geometry was based on MRI, whereas the input used for gait loading was obtained from the literature. Femoral and tibial cartilages were modeled as fibril-reinforced poroviscoelastic materials, whereas menisci were considered as transversely isotropic. Fluid pressures in the chondrocyte and cartilage tissue increased up to 2MPa (an increase of 30%) in the meniscectomy joint compared to the normal, healthy joint. The elevated level of fluid pressure was observed during the entire stance phase of gait. A medial meniscectomy caused substantially larger (up to 60%) changes in maximum principal strains in the chondrocyte compared to those in the peri- or extracellular matrices. Chondrocyte volume or morphology did not change substantially due to a medial meniscectomy. Current findings suggest that during walking chondrocyte deformations are not substantially altered due to a medial meniscectomy, while abnormal joint loading exposes chondrocytes to elevated levels of fluid pressure and maximum principal strains (compared to strains in the peri- or extracellular matrices). These might contribute to cell viability and the onset of osteoarthritis.

In vivo cartilage strain increases following medial meniscal tear and correlates with synovial fluid matrix metalloproteinase activity

Meniscal tears are common injuries, and while partial meniscectomy is a frequent treatment option, general meniscus loss is a risk factor for the development of osteoarthritis. The goal of this study was to measure the in vivo tibiofemoral cartilage contact patterns in patients with meniscus tears in relation to biomarkers of cartilage catabolism in the synovial fluid of these joints. A combination of magnetic resonance imaging and biplanar fluoroscopy was used to determine the in vivo motion and cartilage contact mechanics of the knee. Subjects with isolated medial meniscus tears were analyzed while performing a quasi-static lunge, and the contralateral uninjured knee was used as a control. Synovial fluid was collected from the injured knee and matrix metalloproteinase (MMP) activity, sulfated glycosaminoglycan, cartilage oligomeric matrix protein, prostaglandin E2, and the collagen type II cleavage biomarker C2C were measured. Contact strain in the medial compartment increased significantly in the injured knees compared to contralateral control knees. In the lateral compartment, the contact strain in the injured knee was significantly increased only at the maximum flexion angle (105°). The average cartilage strain at maximum flexion positively correlated with total MMP activity in the synovial fluid. These findings show that meniscal injury leads to loss of normal joint function and increased strain of the articular cartilage, which correlated to elevated total MMP activity in the synovial fluid. The increased strain and total MMP activity may reflect, or potentially contribute to, the early development of osteoarthritis that is observed following meniscal injury.

On the relative relevance of subject-specific geometries and degeneration-specific mechanical properties for the study of cell death in human intervertebral disk models

Capturing patient- or condition-specific intervertebral disk (IVD) properties in finite element models is outmost important in order to explore how biomechanical and biophysical processes may interact in spine diseases. However, disk degenerative changes are often modeled through equations similar to those employed for healthy organs, which might not be valid. As for the simulated effects of degenerative changes, they likely depend on specific disk geometries. Accordingly, we explored the ability of continuum tissue models to simulate disk degenerative changes. We further used the results in order to assess the interplay between these simulated changes and particular IVD morphologies, in relation to disk cell nutrition, a potentially important factor in disk tissue regulation. A protocol to derive patient-specific computational models from clinical images was applied to different spine specimens. In vitro, IVD creep tests were used to optimize poro-hyperelastic input material parameters in these models, in function of the IVD degeneration grade. The use of condition-specific tissue model parameters in the specimen-specific geometrical models was validated against independent kinematic measurements in vitro. Then, models were coupled to a transport-cell viability model in order to assess the respective effects of tissue degeneration and disk geometry on cell viability. While classic disk poro-mechanical models failed in representing known degenerative changes, additional simulation of tissue damage allowed model validation and gave degeneration-dependent material properties related to osmotic pressure and water loss, and to increased fibrosis. Surprisingly, nutrition-induced cell death was independent of the grade-dependent material properties, but was favored by increased diffusion distances in large IVDs. Our results suggest that in situ geometrical screening of IVD morphology might help to anticipate particular mechanisms of disk degeneration.

The structure and function of the pericellular matrix of articular cartilage

Chondrocytes in articular cartilage are surrounded by a narrow pericellular matrix (PCM) that is both biochemically and biomechanically distinct from the extracellular matrix (ECM) of the tissue. While the PCM was first observed nearly a century ago, its role is still under investigation. In support of early hypotheses regarding its function, increasing evidence indicates that the PCM serves as a transducer of biochemical and biomechanical signals to the chondrocyte. Work over the past two decades has established that the PCM in adult tissue is defined biochemically by several molecular components, including type VI collagen and perlecan. On the other hand, the biomechanical properties of this structure have only recently been measured. Techniques such as micropipette aspiration, in situ imaging, computational modeling, and atomic force microscopy have determined that the PCM exhibits distinct mechanical properties as compared to the ECM, and that these properties are influenced by specific PCM components as well as disease state. Importantly, the unique relationships among the mechanical properties of the chondrocyte, PCM, and ECM in different zones of cartilage suggest that this region significantly influences the stress-strain environment of the chondrocyte. In this review, we discuss recent advances in the measurement of PCM mechanical properties and structure that further increase our understanding of PCM function. Taken together, these studies suggest that the PCM plays a critical role in controlling the mechanical environment and mechanobiology of cells in cartilage and other cartilaginous tissues, such as the meniscus or intervertebral disc.

A biphasic multiscale study of the mechanical microenvironment of chondrocytes within articular cartilage under unconfined compression

Computational analyses have been used to study the biomechanical microenvironment of the chondrocyte that cannot be assessed by in vitro experimental studies; yet all computational studies thus far have focused on the effect of zonal location (superficial, middle, and deep) on the mechanical microenvironment of chondrocytes. The aim of this paper was to study the effect of both zonal and radial locations on the biomechanical microenvironment of chondrocytes in inhomogeneous cartilage under unconfined stress relaxation. A biphasic multiscale approach was employed and nine chondrocytes in different locations were studied. Hyperelastic biphasic theory and depth-dependent aggregate modulus and permeability of articular cartilage were included in the models. It was found that both zonal and radial locations affected the biomechanical stresses and strains of the chondrocytes. Chondrocytes in the mid-radial location had increased volume during the early stage of the loading process. Maximum principal shear stress at the interface between the chondrocyte and the extracellular matrix (ECM) increased with depth, yet that at the ECM-pericellular matrix (PCM) interface had an inverse trend. Fluid pressure decreased with depth, while the fluid pressure difference between the top and bottom boundaries of the microscale model increased with depth. Regardless of location, fluid was exchanged between the chondrocyte, PCM, and ECM. These findings suggested that even under simple compressive loading conditions, the biomechanical microenvironment of the chondrocytes, PCM and ECM was spatially dependent. The current study provides new insight on chondrocyte biomechanics.

Automated generation of tissue-specific three-dimensional finite element meshes containing ellipsoidal cellular inclusions

Finite element analysis provides a means of describing cellular mechanics in tissue, which can be useful in understanding and predicting physiological and pathological changes. Many prior studies have been limited to simulations of models containing single cells, which may not accurately describe the influence of mechanical interactions between cells. It is desirable to generate models that more accurately reflect the cellular organisation in tissue in order to evaluate the mechanical function of cells. However, as the model geometry becomes more complicated, manual model generation can become laborious. This can be prohibitive if a large number of distinct cell-scale models are required, for example, in multiscale modelling or probabilistic analysis. Therefore, a method was developed to automatically generate tissue-specific cellular models of arbitrary complexity, with minimal user intervention. This was achieved through a set of scripts, which are capable of generating both sample-specific models, with explicitly defined geometry, and tissue-specific models, with geometry derived implicitly from normal statistical distributions. Models are meshed with tetrahedral (TET) elements of variable size to sufficiently discretise model geometries at different spatial scales while reducing model complexity. The ability of TET meshes to appropriately simulate the biphasic mechanical response of a single-cell model is established against that of a corresponding hexahedral mesh for an illustrative use case. To further demonstrate the flexibility of this tool, an explicit model was developed from three-dimensional confocal laser scanning image data, and a set of models were generated from a statistical cellular distribution of the articular femoral cartilage. The tools presented herein are free and openly accessible to the community at large.

High body mass index is associated with increased diurnal strains in the articular cartilage of the knee

OBJECTIVE

Obesity is an important risk factor for osteoarthritis (OA) and is associated with changes in both the biomechanical and inflammatory environments within the joint. However, the relationship between obesity and cartilage deformation is not fully understood. The goal of this study was to determine the effects of body mass index (BMI) on the magnitude of diurnal cartilage strain in the knee.

METHODS

Three-dimensional maps of knee cartilage thickness were developed from 3T magnetic resonance images of the knees of asymptomatic age- and sex-matched subjects with normal BMI (18.5-24.9 kg/m2) or high BMI (25-31 kg/m2). Site-specific magnitudes of diurnal cartilage strain were determined using aligned images recorded at 8:00 AM and 4:00 PM on the same day.

RESULTS

Subjects with high BMI had significantly thicker cartilage on both the patella and femoral groove, as compared to subjects with normal BMI. Diurnal cartilage strains were dependent on location in the knee joint, as well as BMI. Subjects with high BMI, compared to those with normal BMI, exhibited significantly higher compressive strains in the tibial cartilage. Cartilage thickness on both femoral condyles decreased significantly from the AM to the PM time point; however, there was no significant effect of BMI on diurnal cartilage strain in the femur.

CONCLUSION

Increased BMI is associated with increased diurnal strains in articular cartilage of both the medial and lateral compartments of the knee. The increased cartilage strains observed in individuals with high BMI may, in part, explain the elevated risk of OA associated with obesity or may reflect alterations in the cartilage mechanical properties in subjects with high BMI.

Micromechanical mapping of early osteoarthritic changes in the pericellular matrix of human articular cartilage

OBJECTIVE

Osteoarthritis (OA) is a degenerative joint disease characterized by the progressive loss of articular cartilage. While macroscale degradation of the cartilage extracellular matrix (ECM) has been extensively studied, microscale changes in the chondrocyte pericellular matrix (PCM) and immediate microenvironment with OA are not fully understood. The objective of this study was to quantify osteoarthritic changes in the micromechanical properties of the ECM and PCM of human articular cartilage in situ using atomic force microscopy (AFM).

METHOD

AFM elastic mapping was performed on cryosections of human cartilage harvested from both condyles of macroscopically normal and osteoarthritic knee joints. This method was used to test the hypotheses that both ECM and PCM regions exhibit a loss of mechanical properties with OA and that the size of the PCM is enlarged in OA cartilage as compared to normal tissue.

RESULTS

Significant decreases were observed in both ECM and PCM moduli of 45% and 30%, respectively, on the medial condyle of OA knee joints as compared to cartilage from macroscopically normal joints. Enlargement of the PCM, as measured biomechanically, was also observed in medial condyle OA cartilage, reflecting the underlying distribution of type VI collagen in the region. No significant differences were observed in elastic moduli or their spatial distribution on the lateral condyle between normal and OA joints.

CONCLUSION

Our findings provide new evidence of significant site-specific degenerative changes in the chondrocyte micromechanical environment with OA.

The effects of femoral graft placement on cartilage thickness after anterior cruciate ligament reconstruction

Altered joint motion has been thought to be a contributing factor in the long-term development of osteoarthritis after ACL reconstruction. While many studies have quantified knee kinematics after ACL injury and reconstruction, there is limited in vivo data characterizing the effects of altered knee motion on cartilage thickness distributions. Thus, the objective of this study was to compare cartilage thickness distributions in two groups of patients with ACL reconstruction: one group in which subjects received a non-anatomic reconstruction that resulted in abnormal joint motion and another group in which subjects received an anatomically placed graft that more closely restored normal knee motion. Ten patients with anatomic graft placement (mean follow-up: 20 months) and 12 patients with non-anatomic graft placement (mean follow-up: 18 months) were scanned using high-resolution MR imaging. These images were used to generate 3D mesh models of both knees of each patient. The operative and contralateral knee models were registered to each other and a grid sampling system was used to make site-specific comparisons of cartilage thickness. Patients in the non-anatomic graft placement group demonstrated a significant decrease in cartilage thickness along the medial intercondylar notch in the operative knee relative to the intact knee (8%). In the anatomic graft placement group, no significant changes were observed. These findings suggest that restoring normal knee motion after ACL injury may help to slow the progression of degeneration. Therefore, graft placement may have important implications on the development of osteoarthritis after ACL reconstruction.

Atomic force microscopy reveals regional variations in the micromechanical properties of the pericellular and extracellular matrices of the meniscus

Regional variations in the composition and architecture of the extracellular matrix (ECM) and pericellular matrix (PCM) of the knee meniscus play important roles in determining the local mechanical environment of meniscus cells. In this study, atomic force microscopy was used to spatially map the mechanical properties of matched ECM and perlecan-labeled PCM sites within the outer, middle, and inner porcine medial meniscus, and to evaluate the properties of the proximal surface of each region. The elastic modulus of the PCM was significantly higher in the outer region (151.4 ± 38.2 kPa) than the inner region (27.5 ± 8.8 kPa), and ECM moduli were consistently higher than region-matched PCM sites in both the outer (320.8 ± 92.5 kPa) and inner (66.1 ± 31.4 kPa) regions. These differences were associated with a higher proportion of aligned collagen fibers and lower glycosaminoglycan content in the outer region. Regional variations in the elastic moduli and some viscoelastic properties were observed on the proximal surface of the meniscus, with the inner region exhibiting the highest moduli overall. These results indicate that matrix architecture and composition play an important role in the regional micromechanical properties of the meniscus, suggesting that the local stress-strain environment of meniscal cells may vary significantly among the different regions.

Evaluation of a post-processing approach for multiscale analysis of biphasic mechanics of chondrocytes

Understanding the mechanical behaviour of chondrocytes as a result of cartilage tissue mechanics has significant implications for both evaluation of mechanobiological function and to elaborate on damage mechanisms. A common procedure for prediction of chondrocyte mechanics (and of cell mechanics in general) relies on a computational post-processing approach where tissue-level deformations drive cell-level models. Potential loss of information in this numerical coupling approach may cause erroneous cellular-scale results, particularly during multiphysics analysis of cartilage. The goal of this study was to evaluate the capacity of first- and second-order data passing to predict chondrocyte mechanics by analysing cartilage deformations obtained for varying complexity of loading scenarios. A tissue-scale model with a sub-region incorporating representation of chondron size and distribution served as control. The post-processing approach first required solution of a homogeneous tissue-level model, results of which were used to drive a separate cell-level model (same characteristics as the sub-region of control model). The first-order data passing appeared to be adequate for simplified loading of the cartilage and for a subset of cell deformation metrics, for example, change in aspect ratio. The second-order data passing scheme was more accurate, particularly when asymmetric permeability of the tissue boundaries was considered. Yet, the method exhibited limitations for predictions of instantaneous metrics related to the fluid phase, for example, mass exchange rate. Nonetheless, employing higher order data exchange schemes may be necessary to understand the biphasic mechanics of cells under lifelike tissue loading states for the whole time history of the simulation.

Recent advances in computational mechanics of the human knee joint

Computational mechanics has been advanced in every area of orthopedic biomechanics. The objective of this paper is to provide a general review of the computational models used in the analysis of the mechanical function of the knee joint in different loading and pathological conditions. Major review articles published in related areas are summarized first. The constitutive models for soft tissues of the knee are briefly discussed to facilitate understanding the joint modeling. A detailed review of the tibiofemoral joint models is presented thereafter. The geometry reconstruction procedures as well as some critical issues in finite element modeling are also discussed. Computational modeling can be a reliable and effective method for the study of mechanical behavior of the knee joint, if the model is constructed correctly. Single-phase material models have been used to predict the instantaneous load response for the healthy knees and repaired joints, such as total and partial meniscectomies, ACL and PCL reconstructions, and joint replacements. Recently, poromechanical models accounting for fluid pressurization in soft tissues have been proposed to study the viscoelastic response of the healthy and impaired knee joints. While the constitutive modeling has been considerably advanced at the tissue level, many challenges still exist in applying a good material model to three-dimensional joint simulations. A complete model validation at the joint level seems impossible presently, because only simple data can be obtained experimentally. Therefore, model validation may be concentrated on the constitutive laws using multiple mechanical tests of the tissues. Extensive model verifications at the joint level are still crucial for the accuracy of the modeling.

Subject-specific analysis of joint contact mechanics: application to the study of osteoarthritis and surgical planning

Advances in computational mechanics, constitutive modeling, and techniques for subject-specific modeling have opened the door to patient-specific simulation of the relationships between joint mechanics and osteoarthritis (OA), as well as patient-specific preoperative planning. This article reviews the application of computational biomechanics to the simulation of joint contact mechanics as relevant to the study of OA. This review begins with background regarding OA and the mechanical causes of OA in the context of simulations of joint mechanics. The broad range of technical considerations in creating validated subject-specific whole joint models is discussed. The types of computational models available for the study of joint mechanics are reviewed. The types of constitutive models that are available for articular cartilage are reviewed, with special attention to choosing an appropriate constitutive model for the application at hand. Issues related to model generation are discussed, including acquisition of model geometry from volumetric image data and specific considerations for acquisition of computed tomography and magnetic resonance imaging data. Approaches to model validation are reviewed. The areas of parametric analysis, factorial design, and probabilistic analysis are reviewed in the context of simulations of joint contact mechanics. Following the review of technical considerations, the article details insights that have been obtained from computational models of joint mechanics for normal joints; patient populations; the study of specific aspects of joint mechanics relevant to OA, such as congruency and instability; and preoperative planning. Finally, future directions for research and application are summarized.

Visualization of small lesions in rat cartilage by means of laboratory-based x-ray phase contrast imaging

Being able to quantitatively assess articular cartilage in three-dimensions (3D) in small rodent animal models, with a simple laboratory set-up, would prove extremely important for the development of pre-clinical research focusing on cartilage pathologies such as osteoarthritis (OA). These models are becoming essential tools for the development of new drugs for OA, a disease affecting up to 1/3 of the population older than 50 years for which there is no cure except prosthetic surgery. However, due to limitations in imaging technology, high-throughput 3D structural imaging has not been achievable in small rodent models, thereby limiting their translational potential and their efficiency as research tools. We show that a simple laboratory system based on coded-aperture x-ray phase contrast imaging (CAXPCi) can correctly visualize the cartilage layer in slices of an excised rat tibia imaged both in air and in saline solution. Moreover, we show that small, surgically induced lesions are also correctly detected by the CAXPCi system, and we support this finding with histopathology examination. Following these successful proof-of-concept results in rat cartilage, we expect that an upgrade of the system to higher resolutions (currently underway) will enable extending the method to the imaging of mouse cartilage as well. From a technological standpoint, by showing the capability of the system to detect cartilage also in water, we demonstrate phase sensitivity comparable to other lab-based phase methods (e.g. grating interferometry). In conclusion, CAXPCi holds a strong potential for being adopted as a routine laboratory tool for non-destructive, high throughput assessment of 3D structural changes in murine articular cartilage, with a possible impact in the field similar to the revolution that conventional microCT brought into bone research.

Diurnal variations in articular cartilage thickness and strain in the human knee

Due to the biphasic viscoelastic nature of cartilage, joint loading may result in deformations that require times on the order of hours to fully recover. Thus, cartilaginous tissues may exhibit cumulative strain over the course of each day. The goal of this study was to assess the magnitude and spatial distribution of strain in the articular cartilage of the knee with daily activity. Magnetic resonance (MR) images of 10 asymptomatic subjects (six males and four females) with mean age of 29 years were obtained at 8:00 AM and 4:00 PM on the same day using a 3T magnet. These images were used to create 3D models of the femur, tibia, and patella from which cartilage thickness distributions were quantified. Cartilage thickness generally decreased from AM to PM in all areas except the patellofemoral groove and was associated with significant compressive strains in the medial condyle and tibial plateau. From AM to PM, cartilage of the medial tibial plateau exhibited a compressive strain of -5.1±1.0% (mean±SEM) averaged over all locations, while strains in the lateral plateau were slightly lower (-3.1±0.6%). Femoral cartilage showed an average strain of -1.9±0.6%. The findings of this study show that human knee cartilage undergoes diurnal changes in strain that vary with site in the joint. Since abnormal joint loading can be detrimental to cartilage homeostasis, these data provide a baseline for future studies investigating the effects of altered biomechanics on diurnal cartilage strains and cartilage physiology.