3D finite element analysis of nutrient distributions and cell viability in the intervertebral disc: effects of deformation and degeneration

Effect of compressive and tensile forces on glucose concentration and cell viability within the intervertebral disc: A finite element study

Understanding the role of mechanical force on tissue nutrient transport is essential, as sustained force may affect nutrient levels within the disc and initiate disc degeneration. This study aims to evaluate the time-dependent effects of different compressive force amplitudes as well as tensile force on glucose concentration and cell viability within the disc. Based on the mechano-electrochemical mixture theory, a multiphasic finite element model of the lumbar intervertebral disc was developed. The minimum glucose concentration and minimum cell density in both normal and degenerated discs were predicted for different compressive force amplitudes, tensile force, and corresponding creep time. Under high compressive force, the minimum glucose concentration exhibited an increasing and then decreasing trend with creep time in the normal disc, whereas that of the degenerated disc increased, then decreased, and finally increased again. At steady state, a higher compressive force was accompanied by a lower glucose concentration distribution. In the degenerated disc, the minimum cell density was negatively correlated with creep time, with a greater range of affected tissue under a higher compressive force. For tensile force, the minimum glucose concentration of the degenerated disc raised over time. This study highlighted the importance of creep time, force magnitude, and force type in affecting nutrient concentration and cell viability. Sustained weight-bearing activities could deteriorate the nutrient environment of the degenerated disc, while tensile force might have a nonnegligible role in effectively improving nutrient levels within the degenerated disc.

New evidence on the controversy over the correlation between vertebral osteoporosis and intervertebral disc degeneration: a systematic review of relevant animal studies

OBJECTIVE

The effect of vertebral osteoporosis on disc degeneration remains controversial. The aim of this study was to conduct a systematic review and meta-analysis of relevant animal studies to shed more light on the effects and mechanisms of vertebral osteoporosis on disc degeneration and to promote the resolution of the controversy.

METHODS

The PubMed, Cochrane Library, and Embase databases were searched for studies that met the inclusion criteria. Basic information and data were extracted from the included studies and data were analyzed using STATA 15.1 software. This study was registered on INPLASY with the registration number INPLASY202370099 and https://doi.org/10.37766/inplasy2023.7.0099 .

RESULTS

A total of 13 studies were included in our study. Both animals, rats and mice, were covered. Meta-analysis results showed in disc height index (DHI) (P < 0.001), histological score (P < 0.001), number of osteoblasts in the endplate (P = 0.043), number of osteoclasts in the endplate (P < 0.001), type I collagen (P < 0.001), type II collagen (P < 0.001), aggrecan (P < 0.001), recombinant a disintegrin and metalloproteinase with thrombospondin-4 (ADAMTS-4) (P < 0.001), matrix metalloproteinase-1 (MMP-1) (P < 0.001), MMP-3 (P < 0.001), MMP-13 (P < 0.001), the difference between the osteoporosis group and the control group was statistically significant. In terms of disc volume, the difference between the osteoporosis group and the control group was not statistically significant (P = 0.459).

CONCLUSION

Our study shows that vertebral osteoporosis may exacerbate disc degeneration. Abnormal bone remodeling caused by vertebral osteoporosis disrupts the structural integrity of the endplate, leading to impaired nutrient supply to the disc, increased expression of catabolic factors, and decreased levels of type II collagen and aggrecan may be one of the potential mechanisms.

Glucose depletion decreases cell viability without triggering degenerative changes in a physiological nucleus pulposus explant model

Although the etiology of intervertebral disc degeneration is still unresolved, the nutrient paucity resulting from its avascular nature is suspected of triggering degenerative processes in its core: the nucleus pulposus (NP). While severe hypoxia has no significant effects on NP cells, the impact of glucose depletion, such as found in degenerated discs (0.2-1 mM), is still uncertain. Using a pertinent ex-vivo model representative of the unique disc microenvironment, the present study aimed, therefore, at determining the effects of "degenerated" (0.3 mM) glucose levels on bovine NP explant homeostasis. The effects of glucose depletion were evaluated on NP cell viability, apoptosis, phenotype, metabolism, senescence, extracellular matrix anabolism and catabolism, and inflammatory mediator production using fluorescent staining, RT-qPCR, (immuno)histology, ELISA, biochemical, and enzymatic assays. Compared to the "healthy" (2 mM) glucose condition, exposure to the degenerated glucose condition led to a rapid and extensive decrease in NP cell viability associated with increased apoptosis. Although the aggrecan and collagen-II gene expression was also downregulated, NP cell phenotype, and senescence, matrix catabolism, and inflammatory mediator production were not, or only slightly, affected by glucose depletion. The present study provided evidence for glucose depletion as an essential player in NP cell viability but also suggested that other microenvironment factor(s) may be involved in triggering the typical shift of NP cell phenotype observed during disc degeneration. The present study contributes new information for better understanding disc degeneration at the cellular-molecular levels and thus helps to develop relevant therapeutical strategies to counteract it.

ISSLS prize in clinical/bioengineering science 2024: How standing and supine positions influence nutrient transport in human lumbar discs?-A serial post-contrast MRI study evaluating interplay between convection and diffusion

PURPOSE

The intervertebral disc being avascular depends on diffusion and load-based convection for essential nutrient supply and waste removal. There are no reliable methods to simultaneously investigate them in humans under natural loads. For the first time, present study aims to investigate this by strategically employing positional MRI and post-contrast studies in three physiological positions: supine, standing and post-standing recovery.

METHODS

A total of 100 healthy intervertebral discs from 20 volunteers were subjected to a serial post-contrast MR study after injecting 0.3 mmol/kg gadodiamide and T1-weighted MR images were obtained at 0, 2, 6, 12 and 24 h. At each time interval, images were obtained in three positions, i.e. supine, standing and post-standing recovery supine. The signal intensity values at endplate zone and nucleus pulposus were measured. Enhancement percentages were calculated and analysed comparing three positions.

RESULTS

During unloaded supine position, there was slow gradual increase in enhancement reaching peak at 6 h. When the subjects assumed standing position, there was immediate loss of enhancement at nucleus pulposus which resulted in reciprocal increase in enhancement at endplate zone (washout phenomenon). Interestingly, when subjects assumed the post-standing recovery position, the nucleus pulposus regained the enhancement and endplate zone showed reciprocal loss (pumping-in phenomenon).

CONCLUSIONS

For the first time, present study documented acute effects of physiological loading and unloading on nutrition of human discs in vivo. While during rest, solutes diffused gradually into disc, the diurnal short loading and unloading redistribute small solutes by convection. Standing caused rapid solute depletion but promptly regained by assuming resting supine position.

Analysis of Extracellular ATP Distribution in the Intervertebral Disc

Progressive loss of proteoglycans (PGs) is the major biochemical change during intervertebral disc (IVD) degeneration. Adenosine triphosphate (ATP) as the primary energy source is not only critical for cell survival but also serves as a building block in PG synthesis. Extracellular ATP can mediate a variety of physiological functions and was shown to promote extracellular matrix (ECM) production in the IVD. Therefore, the objective of this study was to develop a 3D finite element model to predict extracellular ATP distribution in the IVD and evaluate the impact of degeneration on extracellular ATP distribution. A novel 3D finite element model of the IVD was developed by incorporating experimental measurements of ATP metabolism and ATP-PG binding kinetics into the mechano-electrochemical mixture theory. The new model was validated by experimental data of porcine IVD, and then used to analyze the extracellular distribution of ATP in human IVDs. Extracellular ATP was shown to bind specifically with PGs in IVD ECM. It was found that annulus fibrosus cells hydrolyze ATP faster than that of nucleus pulposus (NP) cells whereas NP cells exhibited a higher ATP release. The distribution of extracellular ATP in a porcine model was consistent with experimental data in our previous study. The predictions from a human IVD model showed a high accumulation of extracellular ATP in the NP region, whereas the extracellular ATP level was reduced with tissue degeneration. This study provides an understanding of extracellular ATP metabolism and its potential biological influences on the IVD via purinergic signaling.

Lactic acid promotes nucleus pulposus cell senescence and corresponding intervertebral disc degeneration via interacting with Akt

The accumulation of metabolites in the intervertebral disc is considered an important cause of intervertebral disc degeneration (IVDD). Lactic acid, which is a metabolite that is produced by cellular anaerobic glycolysis, has been proven to be closely associated with IVDD. However, little is known about the role of lactic acid in nucleus pulposus cells (NPCs) senescence and oxidative stress. The aim of this study was to investigate the effect of lactic acid on NPCs senescence and oxidative stress as well as the underlying mechanism. A puncture-induced disc degeneration (PIDD) model was established in rats. Metabolomics analysis revealed that lactic acid levels were significantly increased in degenerated intervertebral discs. Elimination of excessive lactic acid using a lactate oxidase (LOx)-overexpressing lentivirus alleviated the progression of IVDD. In vitro experiments showed that high concentrations of lactic acid could induce senescence and oxidative stress in NPCs. High-throughput RNA sequencing results and bioinformatic analysis demonstrated that the induction of NPCs senescence and oxidative stress by lactic acid may be related to the PI3K/Akt signaling pathway. Further study verified that high concentrations of lactic acid could induce NPCs senescence and oxidative stress by interacting with Akt and regulating its downstream Akt/p21/p27/cyclin D1 and Akt/Nrf2/HO-1 pathways. Utilizing molecular docking, site-directed mutation and microscale thermophoresis assays, we found that lactic acid could regulate Akt kinase activity by binding to the Lys39 and Leu52 residues in the PH domain of Akt. These results highlight the involvement of lactic acid in NPCs senescence and oxidative stress, and lactic acid may become a novel potential therapeutic target for the treatment of IVDD.

Does vertebral osteoporosis delay or accelerate lumbar disc degeneration? A systematic review

The effect of vertebral osteoporosis on disc degeneration is still debated. The purpose of this study was to provide a systematic review of studies in this area to further reveal the relationship between the two. Relevant studies were searched in electronic databases, and studies were screened according to inclusion and exclusion criteria, and finally, basic information of the included studies was extracted and summarized. This study was conducted according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) statement. A total of 34 publications spanning 24 years were included in our study. There were 19 clinical studies, including 12 prospective studies and 7 retrospective studies. Of these, 7 considered vertebral osteoporosis to be positively correlated with disc degeneration, 8 considered them to be negatively correlated, and 4 considered them to be uncorrelated. Two cadaveric studies were included, one considered the two to be negatively correlated and one considered them not to be correlated. Seven animal studies were included, of which five considered a positive correlation between vertebral osteoporosis and disc degeneration and two considered a negative correlation between the two. There were also 6 studies that used anti-osteoporosis drugs for intervention, all of them were animal studies. Five of them concluded that vertebral osteoporosis was positively associated with disc degeneration, and the remaining one concluded that there was no correlation between the two. Our systematic review shows that the majority of studies currently consider an association between vertebral osteoporosis and disc degeneration, but there is still a huge disagreement whether this association is positive or negative. Differences in observation time and follow-up time may be one of the reasons for the disagreement. A large number of clinical and basic studies are still needed in the future to further explore the relationship between the two.

Diffusivity of Human Cartilage Endplates in Healthy and Degenerated Intervertebral Disks

The cartilage endplates (CEPs) on the superior and inferior surfaces of the intervertebral disk (IVD), are the primary nutrient transport pathways between the disk and the vertebral body. Passive diffusion is responsible for transporting small nutrient and metabolite molecules through the avascular CEPs. The baseline solute diffusivities in healthy CEPs have been previously studied, however alterations in CEP diffusion associated with IVD degeneration remain unclear. This study aimed to quantitatively compare the solute diffusion in healthy and degenerated human CEPs using a fluorescence recovery after photobleaching (FRAP) approach. Seven healthy CEPs and 22 degenerated CEPs were collected from five fresh-frozen human cadaveric spines and 17 patients undergoing spine fusion surgery, respectively. The sodium fluorescein diffusivities in CEP radial and vertical directions were measured using the FRAP method. The CEP calcification level was evaluated by measuring the average X-ray attenuation. No difference was found in solute diffusivities between radial and axial directions in healthy and degenerated CEPs. Compared to healthy CEPs, the average solute diffusivity was 44% lower in degenerated CEPs (Healthy: 29.07 μm2/s (CI: 23.96-33.62 μm2/s); degenerated: 16.32 μm2/s (CI: 13.84-18.84 μm2/s), p < 0.001). The average solute diffusivity had an inverse relationship with the degree of CEP calcification as determined by the normalized X-ray attenuation values (ß = -22.19, R2 = 0.633; p < 0.001). This study suggests that solute diffusion through the disk and vertebral body interface is significantly hindered by CEP calcification, providing clues to help further understand the mechanism of IVD degeneration.

Dynamic loading leads to increased metabolic activity and spatial redistribution of viable cell density in nucleus pulposus tissue

Background

Nucleus pulposus (NP) cell density is orchestrated by an interplay between nutrient supply and metabolite accumulation. Physiological loading is essential for tissue homeostasis. However, dynamic loading is also believed to increase metabolic activity and could thereby interfere with cell density regulation and regenerative strategies. The aim of this study was to determine whether dynamic loading could reduce the NP cell density by interacting with its energy metabolism.

Methods

Bovine NP explants were cultured in a novel NP bioreactor with and without dynamic loading in milieus mimicking the pathophysiological or physiological NP environment. The extracellular content was evaluated biochemically and by Alcian Blue staining. Metabolic activity was determined by measuring glucose and lactate in tissue and medium supernatants. A lactate-dehydrogenase staining was performed to determine the viable cell density (VCD) in the peripheral and core regions of the NP.

Results

The histological appearance and tissue composition of NP explants did not change in any of the groups. Glucose levels in the tissue reached critical values for cell survival (≤0.5 mM) in all groups. Lactate released into the medium was increased in the dynamically loaded compared to the unloaded groups. While the VCD was unchanged on Day 2 in all regions, it was significantly reduced in the dynamically loaded groups on Day 7 (p ≤ 0.01) in the NP core, which led to a gradient formation of VCD in the group with degenerated NP milieu and dynamic loading (p ≤ 0.05).

Conclusion

It was demonstrated that dynamic loading in a nutrient deprived environment similar to that during IVD degeneration can increase cell metabolism to the extent that it was associated with changes in cell viability leading to a new equilibrium in the NP core. This should be considered for cell injections and therapies that lead to cell proliferation for treatment of IVD degeneration.

Controversies in spine research: Organ culture versus in vivo models for studies of the intervertebral disc

Intervertebral disc degeneration is a common cause of low back pain, the leading cause of disability worldwide. Appropriate preclinical models for intervertebral disc research are essential to achieving a better understanding of underlying pathophysiology and for the development, evaluation, and translation of more effective treatments. To this end, in vivo animal and ex vivo organ culture models are both widely used by spine researchers; however, the relative strengths and weaknesses of these two approaches are a source of ongoing controversy. In this article, members from the Spine and Preclinical Models Sections of the Orthopedic Research Society, including experts in both basic and translational spine research, present contrasting arguments in support of in vivo animal models versus ex vivo organ culture models for studies of the disc, supported by a comprehensive review of the relevant literature. The objective is to provide a deeper understanding of the respective advantages and limitations of these approaches, and advance the field toward a consensus with respect to appropriate model selection and implementation. We conclude that complementary use of several model types and leveraging the unique advantages of each is likely to result in the highest impact research in most instances.

Consolidating and re-evaluating the human disc nutrient microenvironment

Background

Despite exciting advances in regenerative medicine, cell‐based strategies for treating degenerative disc disease remain in their infancy. To maximize the potential for successful clinical translation, a more thorough understanding of the in vivo microenvironment is needed to better determine and predict how cell therapies will respond when administered in vivo.

Aims

This work aims to reflect on the in vivo nutrient microenvironment of the degenerating IVD through consolidating what has already been measured together with investigative in silico models.

Materials and Methods

This work uses in silico modeling, underpinned by more recent experimentally determined parameters of degeneration and nutrient transport from the literature, to re‐evaluate the current knowledge in terms of grade‐specific stages of degeneration.

Results

Through modeling only the metabolically active cell population, this work predicts slightly higher glucose concentrations compared to previous in silico models, while the predicted results show good agreement with previous intradiscal pH and oxygen measurements. Increasing calcification with degeneration limits nutrient transport into the IVD and initiates a build‐up of acidity; however, its effect is compensated somewhat by a reduction in diffusional distance due to decreasing disc height.

Discussion

This work advances in silico modeling through a strong foundation of experimentally determined grade‐specific input parameters. Taken together, pre‐existing measurements and predicted results suggest that metabolite concentrations may not be as critically low as commonly believed, with calcification not appearing to have a detrimental effect at stages of degeneration when cell therapies are an appropriate intervention.

Conclusion

Overall, our initiative is to provoke greater deliberation and consideration of the nutrient microenvironment when performing in vitro cell culture and cell therapy development. This work highlights urgency for robust experimental glucose measurements in healthy and degenerating IVDs, not only to validate in silico models but to significantly advance the field in fully elucidating the nutrient microenvironment and refining in vitro techniques to accelerate clinical translation.

Computational Modeling Intervertebral Disc Pathophysiology: A Review

Lower back pain is a medical condition of epidemic proportion, and the degeneration of the intervertebral disc has been identified as a major contributor. The etiology of intervertebral disc (IVD) degeneration is multifactorial, depending on age, cell-mediated molecular degradation processes and genetics, which is accelerated by traumatic or gradual mechanical factors. The complexity of such intertwined biochemical and mechanical processes leading to degeneration makes it difficult to quantitatively identify cause–effect relationships through experiments. Computational modeling of the IVD is a powerful investigative tool since it offers the opportunity to vary, observe and isolate the effects of a wide range of phenomena involved in the degenerative process of discs. This review aims at discussing the main findings of finite element models of IVD pathophysiology with a special focus on the different factors contributing to physical changes typical of degenerative phenomena. Models presented are subdivided into those addressing role of nutritional supply, progressive biochemical alterations stemming from an imbalance between anabolic and catabolic processes, aging and those considering mechanical factors as the primary source that induces morphological change within the disc. Limitations of the current models, as well as opportunities for future computational modeling work are also discussed.

A Proof-of-Concept Study Using Numerical Simulations of an Acoustic Spheroid-on-a-Chip Platform for Improving 3D Cell Culture

Microfluidic lab-on-chip devices are widely being developed for chemical and biological studies. One of the most commonly used types of these chips is perfusion microwells for culturing multicellular spheroids. The main challenge in such systems is the formation of substantial necrotic and quiescent zones within the cultured spheroids. Herein, we propose a novel acoustofluidic integrated platform to tackle this bottleneck problem. It will be shown numerically that such an approach is a potential candidate to be implemented to enhance cell viability and shrinks necrotic and quiescent zones without the need to increase the flow rate, leading to a significant reduction in costly reagents' consumption in conventional spheroid-on-a-chip platforms. Proof-of-concept, designing procedures and numerical simulation are discussed in detail. Additionally, the effects of acoustic and hydrodynamic parameters on the cultured cells are investigated. The results show that by increasing acoustic boundary displacement amplitude (d0), the spheroid's proliferating zone enlarges greatly. Moreover, it is shown that by implementing d0  = 0.5 nm, the required flow rate to maintain the necrotic zone below 13% will be decreased 12 times compared to non-acoustic chips.

Investigating the physiological relevance of ex vivo disc organ culture nutrient microenvironments using in silico modeling and experimental validation

BACKGROUND

Ex vivo disc organ culture systems have become a valuable tool for the development and pre-clinical testing of potential intervertebral disc (IVD) regeneration strategies. Bovine caudal discs have been widely selected due to their large availability and comparability to human IVDs in terms of size and biochemical composition. However, despite their extensive use, it remains to be elucidated whether their nutrient microenvironment is comparable to human degeneration.

AIMS

This work aims to create the first experimentally validated in silico model which can be used to predict and characterize the metabolite concentrations within ex vivo culture systems.

MATERIALS & METHODS

Finite element models of cultured discs governed by previously established coupled reaction-diffusion equations were created using COMSOL Multiphysics. Experimental validation was performed by measuring oxygen, glucose and pH levels within discs cultured for 7 days, in a static compression bioreactor.

RESULTS

The in silico model was successfully validated through good agreement between the predicted and experimentally measured concentrations. For an ex vivo organ cultured in high glucose medium (4.5 g/L or 25 mM) and normoxia, a larger bovine caudal disc (Cd1-2 to Cd3-4) had a central concentration of ~2.6 %O2, ~8 mM of glucose and a pH value of 6.7, while the smallest caudal discs investigated (Cd6-7 and Cd7-8), had a central concentration of ~6.5 %O2, ~12 mM of glucose and a pH value of 6.9.

DISCUSSION

This work advances the knowledge of ex vivo disc culture microenvironments and highlights a critical need for optimization and standardization of culturing conditions.

CONCLUSION

Ultimately, for assessment of cell-based therapies and successful clinical translation based on nutritional demands, it is imperative that the critical metabolite values within organ cultures (minimum glucose, oxygen and pH values) are physiologically relevant and comparable to the stages of human degeneration.

Growing a backbone - functional biomaterials and structures for intervertebral disc (IVD) repair and regeneration: challenges, innovations, and future directions

Back pain and associated maladies can account for an immense amount of healthcare cost and loss of productivity in the workplace. In particular, spine related injuries in the US affect upwards of 5.7 million people each year. The degenerative disc disease treatment almost always arises due to a clinical presentation of pain and/or discomfort. Preferred conservative treatment modalities include the use of non-steroidal anti-inflammatory medications, physical therapy, massage, acupuncture, chiropractic work, and dietary supplements like glucosamine and chondroitin. Artificial disc replacement, also known as total disc replacement, is a treatment alternative to spinal fusion. The goal of artificial disc prostheses is to replicate the normal biomechanics of the spine segment, thereby preventing further damage to neighboring sections. Artificial functional disc replacement through permanent metal and polymer-based components continues to evolve, but is far from recapitulating native disc structure and function, and suffers from the risk of unsuccessful tissue integration and device failure. Tissue engineering and regenerative medicine strategies combine novel material structures, bioactive factors and stem cells alone or in combination to repair and regenerate the IVD. These efforts are at very early stages and a more in-depth understanding of IVD metabolism and cellular environment will also lead to a clearer understanding of the native environment which the tissue engineering scaffold should mimic. The current review focusses on the strategies for a successful regenerative scaffold for IVD regeneration and the need for defining new materials, environments, and factors that are so finely tuned in the healthy human intervertebral disc in hopes of treating such a prevalent degenerative process.

Effect of intervertebral disc degeneration on mechanical and electric signals at the interface between disc and vertebra

Intervertebral disc (IVD) degeneration is significantly correlated with the changes in structure and material properties of adjacent vertebral bone, possibly through mechanical and electrical interactions. However, the mechanisms underlying the alteration of the mechanical and electrical environment at the disc-vertebra interface related with disc degeneration have not been well studied. The objective of this study was to numerically investigate the long-term distributions of mechanical and electrical signals on the disc-vertebra interface with disc degeneration. A three-dimensional finite element model of a human lumbar IVD was used to study the mechanical and electric signals at the interface between disc and vertebral body. The disc degeneration was simulated by reducing the nutrition levels on the nucleus pulposus (NP)-vertebra interface and on the annulus fibrosus (AF) periphery to 30% and 60% of its normal values, respectively. In the simulation, the total external mechanical load applied to the disc-vertebra segment was assumed unchanged during disc degeneration. The simulation results showed that the compressive stress of solid matrix changed by up to ~37 kPa on the NP-vertebra interface, while it increased by up to ~32 kPa on the AF-vertebra interface. The shear stress increased by up to ~37 kPa with disc degeneration. The absolute value of the electric potential on the disc-vertebra interface of the disc slightly decreased with the disc degeneration (~0.5 mV). The knowledge of these spatial and temporal variations of the mechanical stresses and electric potential on the disc-vertebra interface is important for understanding the vertebrae adaptation and remodeling during disc degeneration.

Combined Experimental Approach and Finite Element Modeling of Small Molecule Transport Through Joint Synovium to Measure Effective Diffusivity

Trans-synovial solute transport plays a critical role in the clearance of intra-articularly (IA) delivered drugs. In this study, we present a computational finite element model (FEM) of solute transport through the synovium validated by experiments on synovial explants. Unsteady diffusion of urea, a small uncharged molecule, was measured through devitalized porcine and human synovium using custom-built diffusion chambers. A multiphasic computational model was constructed and optimized with the experimental data to extract effective diffusivity for urea within the synovium. A monotonic decrease in urea concentration was observed in the donor bath over time, with an effective diffusivity found to be an order of magnitude lower in synovium versus that measured in free solution. Parametric studies incorporating an intimal cell layer with varying thickness and varying effective diffusivities were performed, revealing a dependence of drug clearance kinetics on both parameters. The findings of this study indicate that the synovial matrix impedes urea solute transport out of the joint with little retention of the solute in the matrix.

Interaction between Mesenchymal Stem Cells and Intervertebral Disc Microenvironment: From Cell Therapy to Tissue Engineering

Low back pain (LBP) in one of the most disabling symptoms affecting nearly 80% of the population worldwide. Its primary cause seems to be intervertebral disc degeneration (IDD): a chronic and progressive process characterized by loss of viable cells and extracellular matrix (ECM) breakdown within the intervertebral disc (IVD) especially in its inner region, the nucleus pulposus (NP). Over the last decades, innovative biological treatments have been investigated in order to restore the original healthy IVD environment and achieve disc regeneration. Mesenchymal stem cells (MSCs) have been widely exploited in regenerative medicine for their capacity to be easily harvested and be able to differentiate along the osteogenic, chondrogenic, and adipogenic lineages and to secrete a wide range of trophic factors that promote tissue homeostasis along with immunomodulation and anti-inflammation. Several in vitro and preclinical studies have demonstrated that MSCs are able to acquire a NP cell-like phenotype and to synthesize structural components of the ECM as well as trophic and anti-inflammatory mediators that may support resident cell activity. However, due to its unique anatomical location and function, the IVD presents distinctive features: avascularity, hypoxia, low glucose concentration, low pH, hyperosmolarity, and mechanical loading. Such conditions establish a hostile microenvironment for both resident and exogenously administered cells, which limited the efficacy of intradiscal cell therapy in diverse investigations. This review is aimed at describing the characteristics of the healthy and degenerated IVD microenvironment and how such features influence both resident cells and MSC viability and biological activity. Furthermore, we focused on how recent research has tried to overcome the obstacles coming from the IVD microenvironment by developing innovative cell therapies and functionalized bioscaffolds.

Transport of Vancomycin and Cefepime Into Human Intervertebral Discs: Quantitative Analyses

STUDY DESIGN

Simulation of antibiotics transport into human intervertebral disc with intravenous infusion.

OBJECTIVE

The objective of this study was to quantitatively investigate antibiotic concentrations in the disc.

SUMMARY OF BACKGROUND DATA

Intravenous infusion of antibiotics is typically used to treat intervertebral disc infection in clinics. However, it is difficult to evaluate the drug concentrations within discs in vivo.

METHODS

A computational model was used in this study. The variation of drug charge with pH was considered in the model. Thirty-minute infusions of two commonly used antibiotics in clinic-vancomycin and cefepime-were numerically investigated. Spatial and temporal concentration distributions of these drugs in both nondegenerated and moderately degenerated discs were calculated.

RESULTS

For intravenous infusion of 1 g vancomycin and 2 g cefepime in 30 minutes repeated every 12 hours, it was predicted that vancomycin concentration in the disc fluctuated between 17.0 and 31.0 times of its minimum inhibitory concentration (1 ug/mL) and cefepime concentration fluctuated between 1.1 and 4.2 times of its minimum inhibitory concentration (i.e., 8 ug/mL) in about 2 days. It was also found that vancomycin concentration in moderately degenerated disc was lower than that in the nondegenerated disc.

CONCLUSION

This study provides quantitative guidance on selecting proper dosage for treating disc infection. The method used in this study could be used to provide quantitative information on transport of other antibiotics and drugs in discs as well.

LEVEL OF EVIDENCE

N/A.

Effects of Oxygen Concentration and Culture Time on Porcine Nucleus Pulposus Cell Metabolism: An in vitro Study

Low back pain is a common ailment that affects millions of individuals each year and is linked to degeneration of the intervertebral discs in the spine. Intervertebral disc degeneration is known to result from an imbalance in anabolic and catabolic activity by disc cells. Due to the avascular nature of the intervertebral disc, oxygen deficiency may occur in the central nucleus pulposus (NP). The resulting hypoxia affects matrix regulation and energy metabolism of disc cells, although the mechanisms are not fully understood. This study investigates in vitro glucose consumption and gene expression by NP cells over time under varying oxygen tensions. Notochordal porcine NP cells were cultured in agarose discs at 21, 5, or 1% oxygen tension for 1, 5, or 10 days. The expression of 10 key matrix genes, as well as Brachyury (T), by NP cells was analyzed using RT-PCR. Glucose consumption was measured using a two-point method. Results show that culture time and oxygen tension significantly affect glucose consumption rates by porcine NP cells. There were also significant changes in T expression based on oxygen level and culture time. The 1% oxygen tension had a significantly higher T expression on day 10 than the other two groups, which may indicate a better maintenance of the notochordal phenotype. MMP 1 and 13 expression increased over time for all groups, while only the 5% group showed an increase over time for MMP 3. TIMP expression followed the direction of MMPs but to a lesser magnitude. Five percent and twenty-one percent oxygen tensions led to decreases in anabolic gene expression while 1% led to increases. Oxygen concentration and culture time significantly impacted glucose consumption rate and the gene expression of matrix regulatory genes with hypoxic conditions most accurately maintaining the proper NP phenotype. This information is valuable not only for understanding disc pathophysiology, but also for harnessing the potential of notochordal NP cells in therapeutic applications.

Effects of diurnal loading on the transport of charged antibiotics into intervertebral discs

The objective of this study was to quantitatively analyze the effect of diurnal loading on the transport of various charged antibiotics into negatively charged human intervertebral disc (IVD). Transport of charged antibiotics into a human lumbar disc was analyzed using a 3D finite element model. The valence (z) of the electrical charge of antibiotics varied from z = +2 (positively charged) to z = -2 (negatively charged). An uncharged antibiotic (z = 0) was used as a control. Cases with transient antibiotic concentration at disc boundaries [to mimic intravenous (IV) infusion] were simulated. Our results showed that diurnal compression increased the concentrations in the nucleus pulposus (NP) region, but degreased the concentrations in the annulus fibrosus (AF) region for all charged or non-charged drugs. The overall concentration (averaged over disc) increased with diurnal compression. The diurnal compression had more effects on negatively charged antibiotics than positively charged ones. For example, at day 5 with diurnal compression, the diurnal compression increased the concentration of negatively charged drug (z = -1) in NP by 18.3%, but only by 6.6% for positively charged one (z = +1). In AF, diurnal compression decreased the concentration by 13.2% for negatively charged drug (z = -1) versus 1.2% for positively charged one (z = +1). Note these percentages are the averaged values over day 5. This study provides quantitative information on understanding the mechanisms of charged drug transport in human IVDs.

Theoretical Evaluation and Experimental Validation of Localized Therapeutic Hypothermia Application to Preserve Residual Hearing After Cochlear Implantation

OBJECTIVES

Cochlear implantation surgery has been shown to result in trauma to inner ear sensory structures, resulting in loss of residual hearing. Localized therapeutic hypothermia has been shown in clinical care to be a neuroprotective intervention. Previously, we have shown in an experimental model that localized hypothermia protects cochlear hair cells and residual hearing function against surgical and cochlear implantation trauma. Using experimental temperature measurements carried out in human cadaver temporal bones and a finite element model of the inner ear, the present study examined the temperature distribution of a custom-designed hypothermia delivery system in the human inner ear organs.

DESIGN

The efficacy of the hypothermia probe and resulting heat distribution across human cochlea and surrounding tissues were modeled in three-dimensional in COMSOL. The geometry and dimensions of inner ear and temporal bones were derived from computed tomographic and magnetic resonance imaging images. Model predictions were compared with experimental observations from five human temporal bones.

RESULTS

In both the modeling and experimental studies, the cochlear temperature was lowered by 4 to 6 °C on the round window from a baseline of 37 °C within 16 to 18 minutes. The model simulations showed uniformly distributed cooling across the cochlea. This study provides insight for design, operation, and protocols for efficacious delivery of mild therapeutic hypothermia to the human cochlea that may significantly benefit patients undergoing surgical cochlear implantation by preserving residual hearing.

CONCLUSION

There was a close correlation between the results of the numerical simulations and experimental observations in this study. Our custom-designed system is capable of effectively providing mild therapeutic hypothermia locally to the human cochlea. When combined with results from in vivo animal experiments, the present study suggests that the application of localized therapeutic hypothermia may hold potential for patients with an aim to preserve residual hearing after cochlear implantation.

Application of Finite Element Analysis for Investigation of Intervertebral Disc Degeneration: from Laboratory to Clinic

Due to the ethical concern and inability to detect inner stress distributions of intervertebral disc (IVD), traditional methods for investigation of intervertebral disc degeneration (IVDD) have significant limitations. Many researchers have demonstrated that finite element analysis (FEA) is an effective tool for the research of IVDD. However, the specific application of FEA for investigation of IVDD has not been systematically elucidated before. In the present review, we summarize the current finite element models (FEM) used for the investigation of IVDD, including the poroelastic nonlinear FEM, diffusive-reactive theory model and cell-activity coupled mechano-electrochemical theory model. We further elaborate the use of FEA for the research of IVDD pathogenesis especially for nutrition and biomechanics associated etiology, and the biological, biomechanical and clinical influences of IVDD. In addition, the application of FEA for evaluation and exploration of various treatments for IVDD is also elucidated. We conclude that FEA is an excellent technique for research of IVDD, which could be used to explore the etiology, biology and biomechanics of IVDD. In the future, FEA may help us to achieve the goal of individualized precision therapy.

A 3-Dimensional Finite Element Analysis of Adjacent Segment Disk Degeneration Induced by Transforaminal Lumbar Interbody Fusion After Pedicle Screw Fixation

BACKGROUND

Transforaminal lumbar interbody fusion (TLIF) is an effective treatment of upper lumbar intervertebral disk herniation. However, its clinical efficacy for adjacent segment disk degeneration (ASDD) remains undefined. Therefore, the biomechanical evaluation of ASDD caused by TLIF after pedicle screw fixation (PSF) was explored via a 3-dimensional (3D) finite element analysis.

METHODS

Computed tomography images of a healthy male adult volunteer were used in this study. A L3-4 3D finite element model (model) was successfully constructed using Pro/E software, which was also used to establish the L4-5 of the TLIF, PSF, and PSF + TLIF models. Under the same loading conditions, the protrusion and retraction of the adjacent intervertebral disk and the stress distribution of the annulus fibrosis, facet joint, and articular process in the TLIF, PSF, and PSF + TLIF models were all compared.

RESULTS

Protrusion and retraction of the adjacent intervertebral disk were more notable in the PSF + TLIF model than in the PSF model under the same loading conditions. The stress of the annulus fibrosis of the PSF + TLIF model was stronger relative to that of the PSF model under flexion, extension, or lateral bending. The stress of the articular process of the PSF + TLIF model was also stronger than that of the PSF model under extension or lateral bending.

CONCLUSIONS

This study provides evidence that TLIF may aggravate ASDD after PSF. Furthermore, the findings provided in this report represent the theoretic basis for the clinical analysis of ASDD caused by TLIF after PSF.

Theoretical Explorations Generate New Hypotheses About the Role of the Cartilage Endplate in Early Intervertebral Disk Degeneration

Altered cell nutrition in the intervertebral disk (IVD) is considered a main cause for disk degeneration (DD). The cartilage endplate (CEP) provides a major path for the diffusion of nutrients from the peripheral vasculature to the IVD nucleus pulposus (NP). In DD, sclerosis of the adjacent bony endplate is suggested to be responsible for decreased diffusion and disk cell nutrition. Yet, experimental evidence does not support this hypothesis. Hence, we evaluated how moderate CEP composition changes related to tissue degeneration can affect disk nutrition and cell viability. A novel composition-based permeability formulation was developed for the CEP, calibrated, validated, and used in a mechano-transport finite element IVD model. Fixed solute concentrations were applied at the outer surface of the annulus and the CEP, and three cycles of daily mechanical load were simulated. The CEP model indicated that CEP permeability increases with the degeneration/aging of the tissue, in accordance with recent measurements reported in the literature. Additionally, our results showed that CEP degeneration might be responsible for mechanical load-induced NP dehydration, which locally affects oxygen and lactate levels, and reduced glucose concentration by 16% in the NP-annulus transition zone. Remarkably, CEP degeneration was a condition sine-qua-non to provoke cell starvation and death, while simulating the effect of extracellular matrix depletion in DD. This theoretical study cast doubts about the paradigm that CEP calcification is needed to provoke cell starvation, and suggests an alternative path for DD whereby the early degradation of the CEP plays a key role.

COMPUTATION OF THE DYNAMIC COMPRESSION EFFECTS IN SPINE DISCS USING INTEGRAL METHODS

The computational modeling using integral methods of dynamic loading and its effects on the nutrients transport in spine discs is addressed in this work. The numerical simulation and analysis was carried out using the Boundary Element Method (BEM) and a 3D model (axisymmetric) of the disc. The boundary model was discretized using linear interpolated elements and a multi-region approach. Concentration and production of three nutrients as lactate, oxygen and glucose were obtained. The maximum lactate concentration was observed very close to the interface between the nucleus and the inner annulus. A relatively simple model discretized with 130 boundary elements yielded very similar results to these coming from more complex FEM-based models. The numerical efforts in the domain and boundary discretizations were optimized using the BEM. Our results are in good agreement with those obtained using with finite element-based models. As expected, the dynamic loading increased the oxygen–glucose consumption and the lactate production, thus leading to a poor oxygen–glucose concentration at the nucleus of the disc. All of that is a favorable environment for a disc degeneration mechanism to be developed.

Mechanotransduction and cell biomechanics of the intervertebral disc

Mechanical loading of the intervertebral disc (IVD) initiates cell-mediated remodeling events that contribute to disc degeneration. Cells of the IVD, nucleus pulposus (NP) and anulus fibrosus (AF), will exhibit various responses to different mechanical stimuli which appear to be highly dependent on loading type, magnitude, duration, and anatomic zone of cell origin. Cells of the NP, the innermost region of the disc, exhibit an anabolic response to low-moderate magnitudes of static compression, osmotic pressure, or hydrostatic pressure, while higher magnitudes promote a catabolic response marked by increased protease expression and activity. Cells of the outer AF are responsive to physical forces in a manner that depends on frequency and magnitude, as are cells of the NP, though they experience different forces, deformations, pressure, and osmotic pressure in vivo. Much remains to be understood of the mechanotransduction pathways that regulate IVD cell responses to loading, including responses to specific stimuli and also differences among cell types. There is evidence that cytoskeletal remodeling and receptor-mediated signaling are important mechanotransduction events that can regulate downstream effects like gene expression and posttranslational biosynthesis, all of which may influence phenotype and bioactivity. These and other mechanotransduction events will be regulated by known and to-be-discovered cell-matrix and cell-cell interactions, and depend on composition of extracellular matrix ligands for cell interaction, matrix stiffness, and the phenotype of the cells themselves. Here, we present a review of the current knowledge of the role of mechanical stimuli and the impact upon the cellular response to loading and changes that occur with aging and degeneration of the IVD.

Intervertebral Disk Nutrients and Transport Mechanisms in Relation to Disk Degeneration: A Narrative Literature Review

Objective

The purpose of this paper was to review the literature regarding the mechanisms leading to degeneration in intervertebral disks and to discuss contributing mechanical and biological factors.

Methods

The inclusion criteria for the literature review were research studies conducted in the last 3 decades with free full-text available in English. Review articles and articles pertaining to temporomandibular joints and joints of the body other than the intervertebral disk were excluded. The following databases were searched: PubMed, EBSCOhost, and Google Scholar through September 9, 2016.

Results

A total of 57 articles were used in this review. Intervertebral disk cells require glucose for sustainability and oxygen to synthesize matrix components. Nutrients enter the disk via 2 vascular supply routes: capillary beds of end plates and the peripheral annulus fibrosus. Solute size, shape and charge, compression, and metabolic demand all influence the efficiency of nutrient transport, and alterations of any of these factors may have effects on nutrient transport and, potentially, disk degeneration.

Conclusions

Progressive nutrient transport disruptions may actively contribute in advancing the phases of degenerative disk disease. Such disruptions include dysfunctional loading and spinal position, lack of motion, high frequency loading, disk injury, aging, smoking, an acidic environment, and a lack of nutrient bioavailability.

Prediction of Necrotic Core and Hypoxic Zone of Multicellular Spheroids in a Microbioreactor with a U-Shaped Barrier

Microfluidic devices have been widely used for biological and cellular studies. Microbioreactors for three-dimensional (3D) multicellular spheroid culture are now considered as the next generation in in vitro diagnostic tools. The feasibility of using 3D cell aggregates to form multicellular spheroids in a microbioreactor with U-shaped barriers has been demonstrated experimentally. A barrier array is an alternative to commonly used microwell traps. The present study investigates oxygen and glucose concentration distributions as key parameters in a U-shaped array microbioreactor using finite element simulation. The effect of spheroid diameter, inlet concentration and flow rate of the medium are systematically studied. In all cases, the channel walls are considered to be permeable to oxygen. Necrotic and hypoxic or quiescent regions corresponding to both oxygen and glucose concentration distributions are identified for various conditions. The results show that the entire quiescent and necrotic regions become larger with increasing spheroid diameter and decreasing inlet and wall concentration. The shear stress (0.5⁻9 mPa) imposed on the spheroid surface by the fluid flow was compared with the critical values to predict possible damage to the cells. Finally, optimum range of medium inlet concentration (0.13⁻0.2 mM for oxygen and 3⁻11 mM for glucose) and flow rate (5⁻20 μL/min) are found to form the largest possible multicellular spheroid (500 μm), without any quiescent and necrotic regions with an acceptable shear stress. The effect of cell-trap types on the oxygen and glucose concentration inside the spheroid was also investigated. The levels of oxygen and glucose concentration for the microwell are much lower than those for the other two traps. The U-shaped barrier created with microposts allows for a continuous flow of culture medium, and so improves the glucose concentration compared to that in the integrated U-shaped barrier. Oxygen concentration for both types of U-shaped barriers is nearly the same. Due to the advantage of using U-shaped barriers to culture multicellular spheroids, the results of this paper can help to choose the experimental and design parameters of the microbioreactor.

A method for characterising human intervertebral disc glycosaminoglycan disaccharides using liquid chromatography-mass spectrometry with multiple reaction monitoring

Intervertebral disc (IVD) degeneration results in the depletion of proteoglycans and glycosaminoglycans (GAGs), which can lead to structural and mechanical loss of IVD function, ingrowth of nociceptive nerve fibres and eventually discogenic pain. Specific GAG types as well as their disaccharide patterns can be predictive of disease and degeneration in several tissues but have not been comprehensively studied within the IVD. A highly sensitive mass spectrometry based technique with multiple reaction monitoring (MRM) was used to provide characterisation of chondroitin sulphate (CS), hyaluronic acid (HA), heparan sulphate (HS) and their disaccharide sulphation patterns across different anatomical regions of human IVDs. Principal component analysis further distinguished important regional variations and proposed potential ageing variations in GAG profiles. CS was the GAG in greatest abundance in the IVD followed by HA and HS. Principal component analysis identified clear separation of GAG profiles between nucleus pulposus and annulus fibrosus in young and old specimens. Distinct patterns of predominantly expressed disaccharides of CS and HS between young and old IVD samples, provided preliminary evidence that important alterations in disaccharides occur within IVDs during ageing. This technique offered a novel approach to identify and quantify specific GAG disaccharides in human IVDs and the data presented were the first to offer insight into the spatial distribution as well as association with ageing of GAGs and GAG disaccharide sulphation patterns across the human IVD.

The nutrition of the human meniscus: A computational analysis investigating the effect of vascular recession on tissue homeostasis

The meniscus is essential to the functioning of the knee, offering load support, congruency, lubrication, and protection to the underlying cartilage. Meniscus degeneration affects ∼35% of the population, and potentially leads to knee osteoarthritis. The etiology of meniscal degeneration remains to be elucidated, although many factors have been considered. However, the role of nutritional supply to meniscus cells in the pathogenesis of meniscus degeneration has been so far overlooked. Nutrients are delivered to meniscal cells through the surrounding synovial fluid and the blood vessels present in the outer region of the meniscus. During maturation, vascularization progressively recedes up to the outer 10% of the tissue, leaving the majority avascular. It has been hypothesized that vascular recession might significantly reduce the nutrient supply to cells, thus contributing to meniscus degeneration. The objective of this study was to evaluate the effect of vascular recession on nutrient levels available to meniscus cells. This was done by developing a novel computational model for meniscus homeostasis based on mixture theory. It was found that transvascular transport of nutrients in the vascularized region of the meniscus contributes to more than 40% of the glucose content in the core of the tissue. However, vascular recession does not significantly alter nutrient levels in the meniscus, reducing at most 5% of the nutrient content in the central portion of the tissue. Therefore, our analysis suggests that reduced vascularity is not likely a primary initiating source in tissue degeneration. However, it does feasibly play a key role in inability for self-repair, as seen clinically.

Simulation of water content distributions in degenerated human intervertebral discs

The objective of this study was to investigate the spatial and temporal variations of water content in intervertebral discs during degeneration and repair processes. We hypothesized that the patterns of water content distribution in the discs are related to the intensity patterns observed in T2-weighted MRI images. Water content distributions in the mildly (e.g., 80% viable cells in the disc, 2.3% decrease in disc height) and moderately (e.g., 40% viable cells in the disc, 9.3% decrease in disc height) degenerated discs were predicted using a finite element model. The variation of water content in the degenerated discs treated with three biological therapies (i.e., increasing the cell density in the nucleus pulposus [Case I], increasing glycosaminoglycan synthesis rate in the nucleus pulposus [Case II], and decreasing glycosaminoglycan degradation rate in the nucleus pulposus [Case III]) were also predicted. It was found that two patterns of water content distributions, a horizontal region with lower water content at the mid-axial plane of nucleus pulposus and a spot with higher water content at the posterior region, were shown during the degeneration progress for the disc simulated in this study. These two patterns disappeared after treatment in Case I, but in Case II and Case III. The implication of these patterns for the horizontal gray band and high intensity zone in T2-weighted MRI images was discussed. This study provided new guidance to develop a novel method for diagnosing disc degeneration and assessing outcomes of biological therapies with MRI techniques. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:147-153, 2017.

A computational model for investigating the effects of changes in bioavailability of insulin-like growth factor-1 on the homeostasis of the intervertebral disc

Insulin-like growth factor-1 (IGF-1) is well-known for upregulating cell proliferation and biosynthesis of the extracellular matrix in the intervertebral disc (IVD). Pathological conditions, such as obesity or chronic kidney disease cause IGF-1 deficiency in plasma. How this deficiency impacts disc homeostasis remains unknown. Pro-anabolic approaches for the treatment of disc degeneration based on enhancing IGF-1 bioavailability to tissue-cells are considered, but knowledge of their effectiveness in enhancing cellular anabolism of a degenerated disc is limited. In this study, we developed a computational model for disc homeostasis specifically addressing the role of IGF-1 in modulating both extracellular matrix biosynthesis and cellularity in the IVD. This model was applied to investigate how changes in IGF-1 bioavailability, namely deficiency or enhancement of growth factor, affect disc health. In this study, it was found that IGF-1 deficiency mainly affects the biosynthesis of ECM components, especially in the most external regions of the IVD such as the cartilage endplates and the outer portion of annulus fibrosus. Also, a total of three approaches for increasing IGF-1 bioavailability as a therapy for degenerated IVDs were investigated. It was found that all these strategies are only beneficial to those disc regions receiving sufficient nutritional supply (i.e., the outmost IVD regions), while they exacerbate tissue degradation in malnourished regions (i.e., inner portion of the disc). This suggests that pro-anabolic growth factor-based therapies are limited in that their success strongly depends on an adequate nutritional supply to the IVD tissue, which is not guaranteed in degenerated discs.

MRI quantification of human spine cartilage endplate geometry: Comparison with age, degeneration, level, and disc geometry

Geometry is an important indicator of disc mechanical function and degeneration. While the geometry and associated degenerative changes in the nucleus pulposus and the annulus fibrosus are well-defined, the geometry of the cartilage endplate (CEP) and its relationship to disc degeneration are unknown. The objectives of this study were to quantify CEP geometry in three dimensions using an MRI FLASH imaging sequence and evaluate relationships between CEP geometry and age, degeneration, spinal level, and overall disc geometry. To do so, we assessed the MRI-based measurements for accuracy and repeatability. Next, we measured CEP geometry across a larger sample set and correlated CEP geometric parameters to age, disc degeneration, level, and disc geometry. The MRI-based measures resulted in thicknesses (0.3-1 mm) that are comparable to prior measurements of CEP thickness. CEP thickness was greatest at the anterior/posterior (A/P) margins and smallest in the center. The CEP A/P thickness, axial area, and lateral width decreased with age but were not related to disc degeneration. Age-related, but not degeneration-related, changes in geometry suggest that the CEP may not follow the progression of disc degeneration. Ultimately, if the CEP undergoes significant geometric changes with aging and if these can be related to low back pain, a clinically feasible translation of the FLASH MRI-based measurement of CEP geometry presented in this study may prove a useful diagnostic tool. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1410-1417, 2016.

Kinetics of charged antibiotic penetration into human intervertebral discs: A numerical study

Little quantitative information exists on the kinetics of charged antibiotic penetration into human intervertebral discs (IVD). This information is crucial for determining the dosage to use, timing of administration, and duration of treatment for infected IVDs. The objective of this study was to quantitatively analyze the transport of various charged antibiotics into human lumbar IVDs. Penetration of charged and uncharged antibiotics into a human lumbar disc was analyzed using a 3D finite element model. The valence (z) of the electrical charge of antibiotics varied from z=+2 (positively charged) to z=-2 (negatively charged). An uncharged antibiotic (z=0) was used as a control. Cases with intravenous (IV) administrations of different charged antibiotics were simulated. Our results showed that the electrical charge had great effects on kinetics of an antibiotic penetration into the IVD; with higher concentrations and uptakes for positively charged antibiotics than those for negatively charged ones. This study provides quantitative information on selecting antibiotics for treating intervertebral disc infections.

Simulation of biological therapies for degenerated intervertebral discs

The efficacy of biological therapies on intervertebral disc repair was quantitatively studied using a three-dimensional finite element model based on a cell-activity coupled multiphasic mixture theory. In this model, cell metabolism and matrix synthesis and degradation were considered. Three types of biological therapies-increasing the cell density (Case I), increasing the glycosaminoglycan (GAG) synthesis rate (Case II), and decreasing the GAG degradation rate (Case III)-to the nucleus pulposus (NP) of each of two degenerated discs [one mildly degenerated (e.g., 80% viable cells in the NP) and one severely degenerated (e.g., 30% viable cells in the NP)] were simulated. Degenerated discs without treatment were also simulated as a control. The cell number needed, nutrition level demanded, time required for the repair, and the long-term outcomes of these therapies were analyzed. For Case I, the repair process was predicted to be dependent on the cell density implanted and the nutrition level at disc boundaries. With sufficient nutrition supply, this method was predicted to be effective for treating both mildly and severely degenerated discs. For Case II, the therapy was predicted to be effective for repairing the mildly degenerated disc, but not for the severely degenerated disc. Similar results were predicted for Case III. No change in cell density for Cases II and III were predicted under normal nutrition level. This study provides a quantitative guide for choosing proper strategies of biological therapies for different degenerated discs.

Influences of Nutrition Supply and Pathways on the Degenerative Patterns in Human Intervertebral Disc

STUDY DESIGN

Investigation of the effects of the impairment of different nutritional pathways on the intervertebral disc degeneration patterns in terms of spatial distributions of cell density, glycosaminoglycan content, and water content.

OBJECTIVE

The aim of this study was to test the hypothesis that impairment of different nutritional pathways would result in different degenerative patterns in human discs.

SUMMARY OF BACKGROUND DATA

Impairment of nutritional pathways has been found to affect cell viability in the disc. However, details on how impairment of different nutritional pathways affects the disc degeneration patterns are unknown.

METHODS

A 3D finite element model was used for this study. This finite element method was based on the cell-activity coupled mechano-electrochemical theory for cartilaginous tissues. Impairment of the nutritional pathways was simulated by lowering the nutrition level at the disc boundaries. Effects of the impartment of cartilaginous endplate-nucleus pulposus (CEP-NP) pathway only (Case 1), annulus fibrosus (AF) pathway only (Case 2), and both pathways (Case 3) on disc degeneration patterns were studied.

RESULTS

The predicted critical levels of nutrition for Case 1, Case 2, and Case 3 were around 30%, 20%, and 50% of the reference values, respectively. Below this critical level, the disc degeneration would occur. Disc degeneration appeared mainly in the NP for Case 1, in the outer AF for Case 2, and in both the NP and inner to middle AF for Case 3. For Cases 1 and 3, the loss of water content was primarily located in the mid-axial plane, which is consistent with the horizontal gray band seen in some T2-weighted magnetic resonance imaging (MRI). For the disc geometry used in this study, it was predicted that there existed a high-intensity zone (for Case 3), as seen in some T2-weighted MRI images.

CONCLUSION

Impairment of different nutrition pathways results in different degenerative patterns.

LEVEL OF EVIDENCE

N/A.

Human cartilage endplate permeability varies with degeneration and intervertebral disc site

Despite the critical functions the human cartilage endplate (CEP) plays in the intervertebral disc, little is known about its structural and mechanical properties and their changes with degeneration. Quantifying these changes with degeneration is important for understanding how the CEP contributes to the function and pathology of the disc. Therefore the objectives of this study were to quantify the effect of disc degeneration on human CEP mechanical properties, determine the influence of superior and inferior disc site on mechanics and composition, and simulate the role of collagen fibers in CEP and disc mechanics using a validated finite element model. Confined compression data and biochemical composition data were used in a biphasic-swelling model to calculate compressive extrafibrillar elastic and permeability properties. Tensile properties were obtained by applying published tensile test data to an ellipsoidal fiber distribution. Results showed that with degeneration CEP permeability decreased 50-60% suggesting that transport is inhibited in the degenerate disc. CEP fibers are organized parallel to the vertebrae and nucleus pulposus and may contribute to large shear strains (0.1-0.2) and delamination failure of the CEP commonly seen in herniated disc tissue. Fiber-reinforcement also reduces CEP axial strains thereby enhancing fluid flux by a factor of 1.8. Collectively, these results suggest that the structure and mechanics of the CEP may play critical roles in the solute transport and disc mechanics.

The role of extracellular matrix elasticity and composition in regulating the nucleus pulposus cell phenotype in the intervertebral disc: a narrative review

Intervertebral disc (IVD) disorders are a major contributor to disability and societal health care costs. Nucleus pulposus (NP) cells of the IVD exhibit changes in both phenotype and morphology with aging-related IVD degeneration that may impact the onset and progression of IVD pathology. Studies have demonstrated that immature NP cell interactions with their extracellular matrix (ECM) may be key regulators of cellular phenotype, metabolism and morphology. The objective of this article is to review our recent experience with studies of NP cell-ECM interactions that reveal how ECM cues can be manipulated to promote an immature NP cell phenotype and morphology. Findings demonstrate the importance of a soft (&lt;700 Pa), laminin-containing ECM in regulating healthy, immature NP cells. Knowledge of NP cell-ECM interactions can be used for development of tissue engineering or cell delivery strategies to treat IVD-related disorders.

Finite element modeling of soft tissues: material models, tissue interaction and challenges

BACKGROUND

Musculoskeletal soft tissues, such as articular cartilage, ligaments, knee meniscus and intervertebral disk, have a complex structure, which provides elasticity and capability to support and distribute the body loads. Soft tissues describe an inhomogeneous and multiphasic structure, and exhibit a nonlinear, time-dependent behavior. Their mechanical response is governed by a substance composed of protein fiber-rich and proteoglycan-rich extracellular matrix and interstitial fluid. Protein fibers (e.g. collagen) give the tissue direction dependent stiffness and strength. To investigate these complex biological systems, the use of mathematical tools is well established, alone or in combination with experimental in vitro and in vivo tests. However, the development of these models poses many challenges due to the complex structure and mechanical response of soft tissues.

METHODS

Non-systematic literature review.

FINDINGS

This paper provides a summary of different modeling strategies with associated material properties, contact interactions between articulating tissues, validation and sensitivity of soft tissues with special focus on knee joint soft tissues and intervertebral disk. Furthermore, it reviews and discusses some salient clinical findings of reported finite element simulations.

INTERPRETATION

Model studies extensively contributed to the understanding of functional biomechanics of soft tissues. Models can be effectively used to elucidate clinically relevant questions. However, users should be aware of the complexity of such tissues and of the capabilities and limitations of these approaches to adequately simulate a specific in vivo or in vitro phenomenon.

Disc size markedly influences concentration profiles of intravenously administered solutes in the intervertebral disc: a computational study on glucosamine as a model solute

PURPOSE

Tests on animals of different species with large differences in intervertebral disc size are commonly used to investigate the therapeutic efficacy of intravenously injected solutes in the disc. We hypothesize that disc size markedly affects outcome.

METHODS

Here, using a small non-metabolized molecule, glucosamine (GL) as a model solute, we calculate the influence of disc size on transport of GL into rat, rabbit, dog and human discs for 10 h post intravenous-injection. We used transient finite element models and considered an identical GL supply for all animals.

RESULTS

Huge effects of disc size on GL concentration profiles were found. Post-injection GL concentration in the rat disc reached 70% blood concentration within 15 min but remained below 10% in the human disc nucleus throughout. The GL rapidly penetrated post-injection into smaller discs resulting in homogeneous concentrations. In contrast, GL concentration, albeit at much lower levels, increased with time in the human disc with a small outward flux at the annulus periphery at longer periods.

CONCLUSIONS

Changes in the disc size hugely influenced GL concentrations throughout the disc at all regions and times. Increases in administered dose can neither remedy the very low concentration levels in the disc center in larger human disc at early post-injection hours nor alter the substantial differences in concentration profiles estimated among various species. The size effect will only be exacerbated as molecular weight of the solute increases and as the endplate calcifies. Extrapolation of findings from animal to human discs on the efficacy of intravenously administered solutes must proceed with great caution.

Elastic, permeability and swelling properties of human intervertebral disc tissues: A benchmark for tissue engineering

The aim of functional tissue engineering is to repair and replace tissues that have a biomechanical function, i.e., connective orthopaedic tissues. To do this, it is necessary to have accurate benchmarks for the elastic, permeability, and swelling (i.e., biphasic-swelling) properties of native tissues. However, in the case of the intervertebral disc, the biphasic-swelling properties of individual tissues reported in the literature exhibit great variation and even span several orders of magnitude. This variation is probably caused by differences in the testing protocols and the constitutive models used to analyze the data. Therefore, the objective of this study was to measure the human lumbar disc annulus fibrosus (AF), nucleus pulposus (NP), and cartilaginous endplates (CEP) biphasic-swelling properties using a consistent experimental protocol and analyses. The testing protocol was composed of a swelling period followed by multiple confined compression ramps. To analyze the confined compression data, the tissues were modeled using a biphasic-swelling model, which augments the standard biphasic model through the addition of a deformation-dependent osmotic pressure term. This model allows considering the swelling deformations and the contribution of osmotic pressure in the analysis of the experimental data. The swelling stretch was not different between the disc regions (AF: 1.28±0.16; NP: 1.73±0.74; CEP: 1.29±0.26), with a total average of 1.42. The aggregate modulus (Ha) of the extra-fibrillar matrix was higher in the CEP (390kPa) compared to the NP (100kPa) or AF (30kPa). The permeability was very different across tissue regions, with the AF permeability (64 E(-16)m(4)/Ns) higher than the NP and CEP (~5.5 E(-16)m(4)/Ns). Additionally, a normalized time-constant (3000s) for the stress relaxation was similar for all the disc tissues. The properties measured in this study are important as benchmarks for tissue engineering and for modeling the disc's mechanical behavior and transport.

What have we learned from finite element model studies of lumbar intervertebral discs in the past four decades?

Finite element analysis is a powerful tool routinely used to study complex biological systems. For the last four decades, the lumbar intervertebral disc has been the focus of many such investigations. To understand the disc functional biomechanics, a precise knowledge of the disc mechanical, structural and biochemical environments at the microscopic and macroscopic levels is essential. In response to this need, finite element model studies have proven themselves as reliable and robust tools when combined with in vitro and in vivo measurements. This paper aims to review and discuss some salient findings of reported finite element simulations of lumbar intervertebral discs with special focus on their relevance and implications in disc functional biomechanics. Towards this goal, the earlier investigations are presented, discussed and summarized separately in three distinct groups of elastic, multi-phasic transient and transport model studies. The disc overall response as well as the relative role of its constituents are markedly influenced by loading rate, magnitude, combinations/preloads and posture. The nucleus fluid content and pressurizing capacity affect the disc compliance, annulus strains and failure sites/modes. Biodynamics of the disc is affected by not only the excitation characteristics but also preloads, existing mass and nucleus condition. The role of fluid pressurization and collagen fiber stiffening diminish with time during diurnal loading. The endplates permeability influences the time-dependent response of the disc in both loaded and unloaded recovery phases. The transport of solutes is substantially influenced by the disc size, tissue diffusivity and endplates permeability.

Effect of cartilage endplate on cell based disc regeneration: a finite element analysis

This study examines the effects of cartilage endplate (CEP) calcification and the injection of intervertebral disc (IVD) cells on the nutrition distributions inside the human IVD under physiological loading conditions using multiphasic finite element modeling. The human disc was modeled as an inhomogeneous mixture consisting of a charged elastic solid, water, ions (Na+ and Cl-), and nutrient solute (oxygen, glucose and lactate) phases. The effect of the endplate calcification was simulated by a reduction of the tissue porosity (i.e., water volume faction) from 0.60 to 0.48. The effect of cell injection was simulated by increasing the cell density in the nucleus pulposus (NP) region by 50%, 100%, and 150%. Strain-dependent transport properties (e.g., hydraulic permeability and solute diffusivities) were considered to couple the solute transport and the mechanical loading. The simulation results showed that nutrient solute distribution inside the disc is maintained at a stable state during the day and night. The physiological diurnal cyclic loading does not change the nutrient environment in the human IVD. The cartilage endplate plays a significant role in the nutrient supply to human IVD. Calcification of the cartilage endplate significantly reduces the nutrient levels in human IVD. Therefore, in cell based therapy for IVD regeneration, the increased nutrient demand as a result of cell injection needs to be addressed. Excessive numbers of injected cells may cause further deterioration of the nutrient environment in the degenerated disc. This study is important for understanding the pathology of IVD degeneration and providing new insights into cell based therapies for low back pain.

Osteoporosis delays intervertebral disc degeneration by increasing intradiscal diffusive transport of nutrients through both mechanical and vascular pathophysiological pathways

Several studies have demonstrated an inverse correlation between osteoporosis and degenerative disc disease, so that patients with lower bone mass index, despite presenting greater risks of vertebral fractures, would paradoxically present delayed intervertebral disc degeneration. However the exact pathophysiological mechanisms underlying such phenomenon are not yet completely elucidated. In this article the author provides a general scheme to explain the causal relation between osteoporosis and delayed intervertebral disc degeneration by two main pathophysiological pathways: a vascular and a mechanical one. According to such model, osteoporosis positively affects disc nutrient diffusion through several mechanisms such as: increased endplate vascularization, decreased endplate resistance and decreased intradiscal strain. In the sequence a comprehensive review of the current literature on the issue is performed in order to provide a general overview about the current degree of evidence about the role of each factor postulated to be involved in such pathophysiological scheme. Finally the author provides overall directions for future research on the issue with special attention to the causal links which are supported by weak scientific evidence or by evidence from single studies.

Factors regulating viable cell density in the intervertebral disc: blood supply in relation to disc height

The intervertebral disc is an avascular tissue, maintained by a small population of cells that obtain nutrients mainly by diffusion from capillaries at the disc-vertebral body interface. Loss of this nutrient supply is thought to lead to disc degeneration, but how nutrient supply influences viable cell density is unclear. We investigated two factors that influence nutrient delivery to disc cells and hence cell viability: disc height and blood supply. We used bovine caudal discs as our model as these show a gradation in disc height. We found that although disc height varied twofold from the largest to the smallest disc studied, it had no significant effect on cell density, unlike the situation found in articular cartilage. The density of blood vessels supplying the discs was markedly greater for the largest disc than the smallest disc, as was the density of pores allowing capillary penetration through the bony endplate. Results indicate that changes in blood vessels in the vertebral bodies supplying the disc, as well as changes in endplate architecture appear to influence density of cells in intervertebral discs.