Computational Modeling to Predict the Temporal Regulation of Chondrocyte Metabolism in Response to Various Dynamic Compression Regimens

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

Articular cartilage experiences significant mechanical loads during daily activities. Healthy cartilage provides the capacity for load bearing and regulates the mechanobiological processes for tissue development, maintenance, and repair. Experimental studies at multiple scales have provided a fundamental understanding of macroscopic mechanical function, evaluation of the micromechanical environment of chondrocytes, and the foundations for mechanobiological response. In addition, computational models of cartilage have offered a concise description of experimental data at many spatial levels under healthy and diseased conditions, and have served to generate hypotheses for the mechanical and biological function. Further, modeling and simulation provides a platform for predictive risk assessment, management of dysfunction, as well as a means to relate multiple spatial scales. Simulation-based investigation of cartilage comes with many challenges including both the computational burden and often insufficient availability of data for model development and validation. This review outlines recent modeling and simulation approaches to understand cartilage function from a mechanical systems perspective, and illustrates pathways to associate mechanics with biological function. Computational representations at single scales are provided from the body down to the microstructure, along with attempts to explore multiscale mechanisms of load sharing that dictate the mechanical environment of the cartilage and chondrocytes.

Human amniotic membrane as a chondrocyte carrier vehicle/substrate: in vitro study

Human amniotic membrane (HAM) is an established biomaterial used in many clinical applications. However, its use for tissue engineering purposes has not been fully realized. A study was therefore conducted to evaluate the feasibility of using HAM as a chondrocyte substrate/carrier. HAMs were obtained from fresh human placenta and were process to produced air dried HAM (AdHAM) and freeze dried HAM (FdHAM). Rabbit chondrocytes were isolated and expanded in vitro and seeded onto these preparations. Cell proliferation, GAG expression and GAG/cell expression were measured at days 3, 6, 9, 12, 15, 21, and 28. These were compared to chondrocytes seeded onto plastic surfaces. Histological analysis and scanning electron microscopy was performed to observe cell attachment. There was significantly higher cell proliferation rates observed between AdHAM (13-51%, P=0.001) or FdHAM (18-48%, p = 0.001) to chondrocytes in monolayer. Similarly, GAG and GAG/cell expressed in AdHAM (33-82%, p = 0.001; 22-60%, p = 0.001) or FdHAM (41-81%, p = 0.001: 28-60%, p = 0.001) were significantly higher than monolayer cultures. However, no significant differences were observed in the proliferation rates (p = 0.576), GAG expression (p = 0.476) and GAG/cell expression (p = 0.135) between AdHAM and FdHAM. The histology and scanning electron microscopy assessments demonstrates good chondrocyte attachments on both HAMs. In conclusion, both AdHAM and FdHAM provide superior chondrocyte proliferation, GAG expression, and attachment than monolayer cultures making it a potential substrate/carrier for cell based cartilage therapy and transplantation.

Numerical assessment on the effective mechanical stimuli for matrix-associated metabolism in chondrocyte-seeded constructs

The self-regeneration capacity of articular cartilage is limited, due to its avascular and aneural nature. Loaded explants and cell cultures demonstrated that chondrocyte metabolism can be regulated via physiologic loading. However, the explicit ranges of mechanical stimuli that correspond to favourable metabolic response associated with extracellular matrix (ECM) synthesis are elusive.Unsystematic protocols lacking this knowledge produce inconsistent results. This study aims to determine the intrinsic ranges of physical stimuli that increase ECM synthesis and simultaneously inhibit nitric oxide (NO) production in chondrocyte–agarose constructs, by numerically reevaluating the experiments performed by Tsuang et al. (2008). Twelve loading patterns were simulated with poro-elastic finite element models in ABAQUS. Pressure on solid matrix, von Mises stress, maximum principle stress and pore pressure were selected as intrinsic mechanical stimuli.Their development rates and magnitudes at the steady state of cyclic loading were calculated with MATLAB at the construct level. Concurrent increase in glycosaminoglycan and collagen was observed at 2300 Pa pressure and 40 Pa/s pressure rate. Between 0–1500 Pa and 0–40 Pa/s, NO production was consistently positive with respect to controls, whereas ECM synthesis was negative in the same range. A linear correlation was found between pressure rate and NO production (R =0.77). Stress states identified in this study are generic and could be used to develop predictive algorithms for matrix production in agarose–chondrocyte constructs of arbitrary shape, size and agarose concentration. They could also be helpful to increase the efficacy of loading protocols for avascular tissue engineering.

How to optimize maturation in a bioreactor for vascular tissue engineering: focus on a decision algorithm for experimental planning

Bioreactors may become essential tools for developing tissue-engineered organs from cells, with or without scaffolds. Cells seeded in these devices are expected to grow and form a tissue. Up to now, one of the main challenges to successfully develop functional organs is the structural organization of the cells into the scaffold during maturation in the bioreactor. The maturation step is affected by a set of highly interlinked dynamical variables (flow, stress, pH, temperature, and growth factors) in such a way that fixing the optimal environmental conditions becomes very complex. This work focuses on how the experimental parameters in the bioreactor can be optimized through numerical modeling to maximize tissue growth. Genetic programming (GP) and Markov decision processes (MDPs) were used in synergy to generate and take full advantage of a model of the vascular construct growth. The approach consists in formulating a model through GP to explain the growth of the construct and using MDPs to come up with a strategy to yield the best results in the experimental runs. Construct growth was improved, and the regeneration process was better understood in numerical simulations that relied on this control system. Therefore, an advanced numerical controller of this type could become an effective and inexpensive tool for planning experimental work in tissue engineering.

Simulating the growth of articular cartilage explants in a permeation bioreactor to aid in experimental protocol design

Recently a cartilage growth finite element model (CGFEM) was developed to solve nonhomogeneous and time-dependent growth boundary-value problems (Davol et al., 2008, "A Nonlinear Finite Element Model of Cartilage Growth," Biomech. Model. Mechanobiol., 7, pp. 295-307). The CGFEM allows distinct stress constitutive equations and growth laws for the major components of the solid matrix, collagens and proteoglycans. The objective of the current work was to simulate in vitro growth of articular cartilage explants in a steady-state permeation bioreactor in order to obtain results that aid experimental design. The steady-state permeation protocol induces different types of mechanical stimuli. When the specimen is initially homogeneous, it directly induces homogeneous permeation velocities and indirectly induces nonhomogeneous solid matrix shear stresses; consequently, the steady-state permeation protocol is a good candidate for exploring two competing hypotheses for the growth laws. The analysis protocols were implemented through the alternating interaction of the two CGFEM components: poroelastic finite element analysis (FEA) using ABAQUS and a finite element growth routine using MATLAB. The CGFEM simulated 12 days of growth for immature bovine articular cartilage explants subjected to two competing hypotheses for the growth laws: one that is triggered by permeation velocity and the other by maximum shear stress. The results provide predictions for geometric, biomechanical, and biochemical parameters of grown tissue specimens that may be experimentally measured and, consequently, suggest key biomechanical measures to analyze as pilot experiments are performed. The combined approach of CGFEM analysis and pilot experiments may lead to the refinement of actual experimental protocols and a better understanding of in vitro growth of articular cartilage.

Theoretical analysis of engineered cartilage oxygenation: influence of construct thickness and media flow rate

A novel parallel-plate bioreactor has been shown to modulate the mechanical and biochemical properties of engineered cartilage by the application of fluid-induced shear stress. Flow or perfusion bioreactors may improve tissue development via enhanced transport of nutrients or gases as well as the application of mechanical stimuli, or a combination of these factors. The goal of this study was to complement observed experimental responses to flow by simulating oxygen transport within cartilage constructs of different thicknesses (250 microm or 1 mm). Using numerical computation of convection-diffusion equations, the evaluation of the tissue oxygenation is performed. Four culture conditions are defined based on tissue thickness and flow rates ranging from 0 to approximately 25 mL min(-1). Under these experimental conditions results show a mean oxygen concentration within the tissue varying from 0.01 to 0.19 mol m(-3) as a function of the tissue thickness and the magnitude of the applied shear stress. More generally, the influence of shear stress varying (via flow rate modification) from 10(-3) to 10 dynes cm(-2) on the tissue oxygenation is studied. The influence on the results of important physical parameters such as the maximal oxygen consumption rate of cells is discussed. Lastly, the importance of oxygen concentration in the lower chamber and its relevance to tissue oxygenation are highlighted by the model results.