Feasibility of agent-based modelling of articular cartilage including a conceptual representation of its structure

Articular cartilage is a complex structure with an architecture in which fluid-swollen proteoglycans constrained within a 3D network of collagen fibrils. Because of the complexity of the cartilage structure, the relationship between its mechanical behaviours at the macroscale level and its components at the micro-scale level are not completely understood. The research objective in this thesis is to create a new model of articular cartilage that can be used to simulate and obtain insight into the micro-macro-interaction and mechanisms underlying its mechanical responses during physiological function. The new model of articular cartilage has two characteristics, namely: i) not use fibre-reinforced composite material idealization ii) Provide a framework for that it does probing the micro mechanism of the fluid-solid interaction underlying the deformation of articular cartilage using simple rules of repartition instead of constitutive / physical laws and intuitive curve-fitting. Even though there are various microstructural and mechanical behaviours that can be studied, the scope of this thesis is limited to osmotic pressure formation and distribution and their influence on cartilage fluid diffusion and percolation, which in turn governs the deformation of the compression-loaded tissue. The study can be divided into two stages. In the first stage, the distributions and concentrations of proteoglycans, collagen and water were investigated using histological protocols. Based on this, the structure of cartilage was conceptualised as microscopic osmotic units that consist of these constituents that were distributed according to histological results. These units were repeated three-dimensionally to form the structural model of articular cartilage. In the second stage, cellular automata were incorporated into the resulting matrix (lattice) to simulate the osmotic pressure of the fluid and the movement of water within and out of the matrix; following the osmotic pressure gradient in accordance with the chosen rule of repartition of the pressure. The outcome of this study is the new model of articular cartilage that can be used to simulate and study the micromechanical behaviours of cartilage under different conditions of health and loading. These behaviours are illuminated at the microscale level using the socalled neighbourhood rules developed in the thesis in accordance with the typical requirements of cellular automata modelling. Using these rules and relevant Boundary Conditions to simulate pressure distribution and related fluid motion produced significant results that provided the following insight into the relationships between osmotic pressure gradient and associated fluid micromovement, and the deformation of the matrix. For example, it could be concluded that: 1. It is possible to model articular cartilage with the agent-based model of cellular automata and the Margolus neighbourhood rule. 2. The concept of 3D inter connected osmotic units is a viable structural model for the extracellular matrix of articular cartilage. 3. Different rules of osmotic pressure advection lead to different patterns of deformation in the cartilage matrix, enabling an insight into how this micromechanism influences macromechanical deformation. 4. When features such as transition coefficient were changed, permeability (representing change) is altered due to the change in concentrations of collagen, proteoglycans (i.e. degenerative conditions), the deformation process is impacted. 5. The boundary conditions also influence the relationship between osmotic pressure gradient and fluid movement at the micro-scale level. The outcomes are important to cartilage research since we can use these to study the microscale damage in the cartilage matrix. From this, we are able to monitor related diseases and their progression leading to potential insight into drug-cartilage interaction for treatment. This innovative model is an incremental progress on attempts at creating further computational modelling approaches to cartilage research and other fluid-saturated tissues and material systems.

Near Infrared Spectroscopy Evaluation of Articular Cartilage : Decision-Making in Arthroplasty : Real-Time Assessment of Articular Cartilage

Current diagnostic methods for assessing the severity of articular cartilage degenerative conditions, such as osteoarthritis, are inadequate. There is also a lack of techniques that can be used for real-time evaluation of the tissue during surgery to inform treatment decision and eliminate subjectivity. This book, derived from Dr Afara’s doctoral research, presents a scientific framework that is based on near infrared (NIR) spectroscopy for facilitating the non-destructive evaluation of articular cartilage health relative to its structural, functional, and mechanical properties. This development is a component of the ongoing research on advanced endoscopic diagnostic techniques in the Articular Cartilage Biomechanics Research Laboratory of Professor Adekunle Oloyede at Queensland University of Technology (QUT), Brisbane Australia.

Load-unloading response of intact and artificially degraded articular cartilage correlated with near infrared (NIR) absorption spectra

The conventional mechanical properties of articular cartilage, such as compressive stiffness, have been demonstrated to be limited in their capacity to distinguish intact (visually normal) from degraded cartilage samples. In this paper, we explore the correlation between a new mechanical parameter, namely the reswelling of articular cartilage following unloading from a given compressive load, and the near infrared (NIR) spectrum. The capacity to distinguish mechanically intact from proteoglycan-depleted tissue relative to the "reswelling" characteristic was first established, and the result was subsequently correlated with the NIR spectral data of the respective tissue samples. To achieve this, normal intact and enzymatically degraded samples were subjected to both NIR probing and mechanical compression based on a load-unload-reswelling protocol. The parameter δ(r), characteristic of the osmotic "reswelling" of the matrix after unloading to a constant small load in the order of the osmotic pressure of cartilage, was obtained for the different sample types. Multivariate statistics was employed to determine the degree of correlation between δ(r) and the NIR absorption spectrum of relevant specimens using Partial Least Squared (PLS) regression. The results show a strong relationship (R(2)=95.89%, p<0.0001) between the spectral data and δ(r). This correlation of δ(r) with NIR spectral data suggests the potential for determining the reswelling characteristics non-destructively. It was also observed that δ(r) values bear a significant relationship with the cartilage matrix integrity, indicated by its proteoglycan content, and can therefore differentiate between normal and artificially degraded proteoglycan-depleted cartilage samples. It is therefore argued that the reswelling of cartilage, which is both biochemical (osmotic) and mechanical (hydrostatic pressure) in origin, could be a strong candidate for characterizing the tissue, especially in regions surrounding focal cartilage defects in joints.

Natural articular joints : model of lamellar-rollerbearing lubrication and the nature of the cartilage surface

The influence of pH on interfacial energy and wettability distributed over the phospholipid bilayer surface were studied, and the importance of cartilage hydrophobicity (wettability) on the coefficient of friction (f) was established. It is argued that the wettability of cartilage signifi antly depends on the number of phospholipid bilayers acting as solid lubricant; the hypothesis was proven by conducting friction tests with normal and lipid- depleted cartilage samples. A lamellar-roller-bearing lubrication model was devised involving two mechanisms: (i) lamellar frictionless movement of bilayers, and (ii) roller-bearing lubrication mode through structured synovial fluid, which operates when lamellar spheres, liposomes and macromolecules act like a roller-bearing situated between two cartilage surfaces in effective biological lubrication.

A study of the diffusion characteristics of normal, delipidized and relipidized articular cartilage using magnetic resonance imaging

This paper assesses the capacity to provide semipermeability of the synthetic layer of surface-active phospholipids created to replace the depleted surface amorphous layer of articular cartilage. The surfaces of articular cartilage specimens in normal, delipidized, and relipidized conditions following incubation in dipalmitoyl-phosphatidylcholine and palmitoyl-oleoyl-phosphatidylcholine components of the joint lipid mixture were characterized nanoscopically with the atomic force microscope and also imaged as deuterium oxide (D2O) diffused transiently through these surfaces in a magnetic resonance imaging enclosure. The MR images were then used to determine the apparent diffusion coefficients in a purpose-built MATLAB®-based algorithm. Our results revealed that all surfaces were permeable to D2O, but that there was a significant difference in the semipermeability of the surfaces under the different conditions, relative to the apparent diffusion coefficients. Based on the results and observations, it can be concluded that the synthetic lipid that is deposited to replace the depleted SAL of articular cartilage is capable of inducing some level of semipermeability.

Feasibility of ranking articular cartilage conditions with non-destructive near-infrared spectroscopy with extension to the Mankin grading system

Grading osteoarthritic tissue has, until now, been a laboratory process confined to research activities. This thesis establishes a scientific protocol that extends osteoarthritic tissue ranking to surgical practice. The innovative protocol, which now incorporates the structural degeneration of collagen, enhances the traditional Modified Mankin ranking system, enabling its application to real time decision during surgery. Because it is fast and without time consuming laboratory process, it would potentially enable the cataloguing of tissues in osteoarthritic joints in all compartments of diseased joints during surgery for epistemological study and insight into the manifestation of osteoarthritis across age, gender, occupation, physical activities and race.

An exploratory study of the potential of resurfacing articular cartilage with synthetic phospholipids

This thesis is aimed at further understanding the uppermost lipid-filled membranous layer (i.e. surface amorphous layer (SAL)) of articular cartilage and to develop a scientific framework for re-introducing lipids onto the surface of lipid-depleted articular cartilage (i.e. "resurfacing"). The outcome will potentially contribute to knowledge that will facilitate the repair of the articular surface of cartilage where degradation is limited to the loss of the lipids of the SAL only. The surface amorphous layer is of utmost importance to the effective load-spreading, lubrication, and semipermeability (which controls its fluid management, nutrient transport and waste removal) of articular cartilage in the mammalian joints. However, because this uppermost layer of cartilage is often in contact during physiological function, it is prone to wear and tear, and thus, is the site for damage initiation that can lead to the early stages of joint condition like osteoarthritis, and related conditions that cause pain and discomfort leading to low quality of life in patients. It is therefore imperative to conduct a study which offers insight into remedying this problem. It is hypothesized that restoration (resurfacing) of the surface amorphous layer can be achieved by re-introducing synthetic surface-active phospholipids (SAPL) into the joint space. This hypothesis was tested in this thesis by exposing cartilage samples whose surface lipids had been depleted to individual and mixtures of synthetic saturated and unsaturated phospholipids. The surfaces of normal, delipidized, and relipidized samples of cartilage were characterized for their structural integrity and functionality using atomic force microscope (AFM), confocal microscope (COFM), Raman spectroscopy, magnetic resonance imaging (MRI) with image processing in the MATLAB® environment and mechanical loading experiments. The results from AFM imaging, confocal microscopy, and Raman spectroscopy revealed a successful deposition of new surface layer on delipidized cartilage when incubated in synthetic phospholipids. The relipidization resulted in a significant improvement in the surface nanostructure of the artificially degraded cartilage, with the complete SAPL mixture providing better outcomes in comparison to those created with the single SAPL components (palmitoyl-oleoyl-phosphatidylcholine, POPC and dipalmitoyl-phosphatidylcholine, DPPC). MRI analysis revealed that the surface created with the complete mixture of synthetic lipids was capable of providing semipermeability to the surface layer of the treated cartilage samples relative to the normal intact surface. Furthermore, deformation energy analysis revealed that the treated samples were capable of delivering the elastic properties required for load bearing and recovery of the tissue relative to the normal intact samples, with this capability closer between the normal and the samples incubated in the complete lipid mixture. In conclusion, this thesis has established that it is possible to deposit/create a potentially viable layer on the surface of cartilage following degradation/lipid loss through incubation in synthetic lipid solutions. However, further studies will be required to advance the ideas developed in this thesis, for the development of synthetic lipid-based injections/drugs for treatment of osteoarthritis and other related joint conditions.

Application of near infrared (NIR) spectroscopy for determining the thickness of articular cartilage

The determination of the characteristics of articular cartilage such as thickness, stiffness and swelling, especially in the form that can facilitate real-time decisions and diagnostics is still a matter for research and development. This paper correlates near infrared spectroscopy with mechanically measured cartilage thickness to establish a fast, non-destructive, repeatable and precise protocol for determining this tissue property. Statistical correlation was conducted between the thickness of bovine cartilage specimens (n = 97) and regions of their near infrared spectra. Nine regions were established along the full absorption spectrum of each sample and were correlated with the thickness using partial least squares (PLS) regression multivariate analysis. The coefficient of determination (R2) varied between 53 and 93%, with the most predictive region (R2 = 93.1%, p < 0.0001) for cartilage thickness lying in the region (wavenumber) 5350–8850 cm−1. Our results demonstrate that the thickness of articular cartilage can be measured spectroscopically using NIR light. This protocol is potentially beneficial to clinical practice and surgical procedures in the treatment of joint disease such as osteoarthritis.

A comparison of the histochemical and image-derived proteoglycan content of articular cartilage

There are several methods for determining the proteoglycan content of cartilage in biomechanics experiments. Many of these include assay-based methods and the histochemistry or spectrophotometry protocol where quantification is biochemically determined. More recently a method based on extracting data to quantify proteoglycan content has emerged using the image processing algorithms, e.g., in ImageJ, to process histological micrographs, with advantages including time saving and low cost. However, it is unknown whether or not this image analysis method produces results that are comparable to those obtained from the biochemical methodology. This paper compares the results of a well-established chemical method to those obtained using image analysis to determine the proteoglycan content of visually normal (n=33) and their progressively degraded counterparts with the protocols. The results reveal a strong linear relationship with a regression coefficient (R2) = 0.9928, leading to the conclusion that the image analysis methodology is a viable alternative to the spectrophotometry.

Simultaneous magnetic resonance imaging and consolidation measurement of articular cartilage

Magnetic resonance imaging (MRI) offers the opportunity to study biological tissues and processes in a non-disruptive manner. The technique shows promise for the study of the load-bearing performance (consolidation) of articular cartilage and changes in articular cartilage accompanying osteoarthritis. Consolidation of articular cartilage involves the recording of two transient characteristics: the change over time of strain and the hydrostatic excess pore pressure (HEPP). MRI study of cartilage consolidation under mechanical load is limited by difficulties in measuring the HEPP in the presence of the strong magnetic fields associated with the MRI technique. Here we describe the use of MRI to image and characterize bovine articular cartilage deforming under load in an MRI compatible consolidometer while monitoring pressure with a Fabry-Perot interferometer-based fiber-optic pressure transducer.

MAS NMR measurements of intact articular bovine cartilage

Articular cartilage (AC), an avascular connective tissue lining articulating surfaces of the long bones, comprises extracellular biopolymers. In functionally compromised states such as osteoarthritis, thinned or lost AC causes reduced mobility and increased health-care costs. Understanding of the characteristics responsible for the load bearing efficiency of AC and the factors leading to its degradation are incomplete. DTI shows the structural alignment of collagen in AC [1] and T2 relaxation measurements suggest that the average director of reorientational motion of water molecules depends on the degree of alignment of collagen in AC [2]. Information on the nature of the chemical interactions involved in functional AC is lacking. The need for AC structural integrity makes solid state NMR an ideal tool to study this tissue. We examined the contribution of water in different functional ‘compartments’ using 1H-MAS, 13C-MAS and 13C-CPMAS NMR of bovine patellar cartilage incubated in D2O. 1H-MAS spectra signal intensity was reduced due to H/D exchange without a measureable redistribution of relative signal intensity. Chemical shift anisotropy was estimated by lineshape analysis of multiple peaks in the 1H-MAS spinning sidebands. These asymmetrical sidebands suggested the presence of multiple water species in AC. Therefore, water was added in small aliquots to D2O saturated AC and the influence of H2O and D2O on organic components was studied with 13C-MAS-NMR and 13C-CPMAS-NMR. Signal intensity in 13C-MAS spectra showed no change in relative signal intensity throughout the spectrum. In 13C-CPMAS spectra, displacement of water by D2O resulted in a loss of signal in the aliphatic region due to a reduction in proton availability for cross-polarization. These results complement dehydration studies of cartilage using osmotic manipulation [3] and demonstrate components of cartilage that are in contact with mobile water.

Solid-state MAS NMR measurements of intact articular bovine cartilage

Articular cartilage (AC), an avascular connective tissue lining articulating surfaces of the long bones, comprises extracellular biopolymers. In functionally compromised states such as osteoarthritis, thinned or lost AC causes reduced mobility and increased health-care costs. Understanding of the characteristics responsible for the load bearing efficiency of AC and the factors leading to its degradation are incomplete. DTI shows the structural alignment of collagen in AC [1] and T2 relaxation measurements suggest that the average director of reorientational motion of water molecules depends on the degree of alignment of collagen in AC [2]. Information on the nature of the chemical interactions involved in functional AC is lacking. The need for AC structural integrity makes solid state NMR an ideal tool to study this tissue. We examined the contribution of water in different functional ‘compartments’ using 1H-MAS, 13C-MAS and 13C-CPMAS NMR of bovine patellar cartilage incubated in D2O. 1H-MAS spectra signal intensity was reduced due to H/D exchange without a measureable redistribution of relative signal intensity. Chemical shift anisotropy was estimated by lineshape analysis of multiple peaks in the 1H-MAS spinning sidebands. These asymmetrical sidebands suggested the presence of multiple water species in AC. Therefore, water was added in small aliquots to D2O saturated AC and the influence of H2O and D2O on organic components was studied with 13C-MAS-NMR and 13C-CPMAS-NMR. Signal intensity in 13C-MAS spectra showed no change in relative signal intensity throughout the spectrum. In 13C-CPMAS spectra, displacement of water by D2O resulted in a loss of signal in the aliphatic region due to a reduction in proton availability for cross-polarization. These results complement dehydration studies of cartilage using osmotic manipulation [3] and demonstrate components of cartilage that are in contact with mobile water.

Near infrared spectroscopy for rapid determination of Mankin score components : a potential tool for quantitative characterization of articular cartilage at surgery

PURPOSE The purpose of this study was to demonstrate the potential of near infrared (NIR) spectroscopy for characterizing the health and degenerative state of articular cartilage based on the components of the Mankin score. METHODS Three models of osteoarthritic degeneration induced in laboratory rats by anterior cruciate ligament (ACL) transection, meniscectomy (MSX), and intra-articular injection of monoiodoacetate (1 mg) (MIA) were used in this study. Degeneration was induced in the right knee joint; each model group consisted of 12 rats (N = 36). After 8 weeks, the animals were euthanized and knee joints were collected. A custom-made diffuse reflectance NIR probe of 5-mm diameter was placed on the tibial and femoral surfaces, and spectral data were acquired from each specimen in the wave number range of 4,000 to 12,500 cm(-1). After spectral data acquisition, the specimens were fixed and safranin O staining (SOS) was performed to assess disease severity based on the Mankin scoring system. Using multivariate statistical analysis, with spectral preprocessing and wavelength selection technique, the spectral data were then correlated to the structural integrity (SI), cellularity (CEL), and matrix staining (SOS) components of the Mankin score for all the samples tested. RESULTS ACL models showed mild cartilage degeneration, MSX models had moderate degeneration, and MIA models showed severe cartilage degenerative changes both morphologically and histologically. Our results reveal significant linear correlations between the NIR absorption spectra and SI (R(2) = 94.78%), CEL (R(2) = 88.03%), and SOS (R(2) = 96.39%) parameters of all samples in the models. In addition, clustering of the samples according to their level of degeneration, with respect to the Mankin components, was also observed. CONCLUSIONS NIR spectroscopic probing of articular cartilage can potentially provide critical information about the health of articular cartilage matrix in early and advanced stages of osteoarthritis (OA). CLINICAL RELEVANCE This rapid nondestructive method can facilitate clinical appraisal of articular cartilage integrity during arthroscopic surgery.

Effect of partial H2O-D2O replacement on the anisotropy of transverse proton spin relaxation in bovine articular cartilage

Anisotropy of transverse proton spin relaxation in collagen-rich tissues like cartilage and tendon is a well-known phenomenon that manifests itself as the "magic-angle" effect in magnetic resonance images of these tissues. It is usually attributed to the non-zero averaging of intra-molecular dipolar interactions in water molecules bound to oriented collagen fibers. One way to manipulate the contributions of these interactions to spin relaxation is by partially replacing the water in the cartilage sample with deuterium oxide. It is known that dipolar interactions in deuterated solutions are weaker, resulting in a decrease in proton relaxation rates. In this work, we investigate the effects of deuteration on the longitudinal and the isotropic and anisotropic contributions to transverse relaxation of water protons in bovine articular cartilage. We demonstrate that the anisotropy of transverse proton spin relaxation in articular cartilage is independent of the degree of deuteration, bringing into question some of the assumptions currently held over the origins of relaxation anisotropy in oriented tissues.