Skip to main content

Potential molecular targets and drugs for basement membranes-related intervertebral disk degeneration through bioinformatics analysis and molecular docking

Abstract

Background

Through bioinformatics analysis to identify the hub genes of Intervertebral disc degeneration (IVDD) associated with basement membranes (BMs) and find out the potential molecular targets and drugs for BMs-related annulus fibrosus (AF) degeneration based on bioinformatic analysis and molecular approach.

Methods

Intervertebral disc degeneration (IVDD) related targets were obtained from GeneCards, DisGenet and OMIM databases. BMs related genes were obtained from Basement membraneBASE database. The intersection targets were identified and subjected to protein-to-protein interaction (PPI) construction via STRING. Hub genes were identified and conducted Gene ontology (GO) and pathway enrichment analysis through MCODE and Clue GO in Cytospace respectively. DSigDB database was retrieved to predict therapeutic drugs and molecular docking was performed through PyMOL, AutoDock 1.5.6 to verify the binding energy between the drug and the different expressed hub genes. Finally, GSE70362 from GEO database was obtained to verify the different expression and correlation of each hub gene for AF degeneration.

Results

We identified 41 intersection genes between 3 disease targets databases and Basement membraneBASE database. PPI network revealed 25 hub genes and they were mainly enriched in GO terms relating to glycosaminoglycan catabolic process, the TGF-β signaling pathway. 4 core targets were found to be significant via comparison of microarray samples and they showed strong correlation. The molecular docking results showed that the core targets have strong binding energy with predicting drugs including chitosamine and retinoic acid.

Conclusions

In this study, we identified hub genes, pathways, potential targets, and drugs for treatment in BMs-related AF degeneration and IVDD.

Peer Review reports

Background

Intervertebral disc degeneration (IVDD) is a leading cause of low back pain, which significantly impacts patients’ quality of life and poses a significant burden on the economy [1, 2]. The annulus fibrosus (AF), a critical component of the intervertebral disc (IVD), encloses the nucleus pulposus (NP) and collaboratively accommodates mechanical loads, providing stability to the spinal column. Notably, IVDD arises due to an exploitation of the structural integrity of the AF before disc degeneration occurs [3]. Although the mechanisms underlying IVDD are not fully understood, an imbalance of anabolic processes in the AF may result in decreased tensile resistance, leading to IVDD [4]. Reduced abundance and quality of extracellular matrix (ECM) are hallmark features of disc degeneration [5]. Basement membranes (BMs) are a type of ECM that form sheet-like structures surrounding most tissues, including the IVD [6, 7]. Melrose et al. [8] used immunolocalization patterns to observe the developing human fetal spine and found that perlecan, the multidomin heparan sulfate proteoglycan of BMs, is a major component of the cartilage primordium in the developing human fetal spine. In their study, aggrecan staining was generally stronger in AF and heparan sulfate was diffusely distributed in the outer AF. Several studies [9, 10] have conducted bioinformatic analyses to investigate the relevant genes and pathways associated with IVDD. However, these studies did not explore the potential role of BM-related molecules in the pathogenesis of IVDD. In addition, AF may serve as a target tissue for the treatment of IVDD. Therapies based on the use of biological molecules, cell-based therapies, gene therapies and tissue engineering are applied in AF in recent decades [11]. Despite extensive research on the pathogenesis of IVDD, there is currently a lack of bioinformatics studies focused on drug prediction in the context of AF repair and regeneration. Therefore, it is imperative to gain a comprehensive understanding of the biological mechanisms underlying IVDD and the molecules associated with its BM, as this knowledge may facilitate the identification of novel therapeutic targets and drugs for AF repair and regeneration.

In this study, we aimed to identify genes associated with IVDD and BMs by conducting an extensive search of multiple disease databases and the latest BMs database [7]. We then constructed a protein-to-protein interaction (PPI) network and identified hub genes. Gene ontology (GO) and pathways enrichment analysis were performed to gain a deeper understanding of the functional roles of these genes. To further validate our findings, we utilized transcriptome sequencing gene data from the Gene Expression Omnibus (GEO) database to verify the hub genes and combined the Drug Signature database (DSigDB) to predict potential therapeutic targets for AF. Finally, we conducted molecular docking to validate the interactions between the predicting drugs and targets. A flowchart detailing the methodology of our study is presented in Fig. 1.

Fig. 1
figure 1

The flow chart of this study. IVDD = intervertebral disc degeneration; BM = basement membrane; GO = Gene Ontology; KEGG = Kyoto Encyclopedia of Genes and Genomes

Methods

Retrieving data

The keyword “intervertebral disc degeneration” was searched in GeneCards database (https://www.genecards.org), DisGeNet (http://www.disgenet.org), and the Online Mendelian Inheritance in Man (OMIM) (http://omim.org) to obtain the disease targets. BMs related genes were obtained from Basement membraneBASE (https://bmbase.manchester.ac.uk/). The R software (version 4.2.0) was used to identify intersection targets associated with IVDD and BMs through Venn diagram.

PPI analysis and hub genes screening

The PPI network was constructed using the Search Tool for the Retrieval of Interacting Genes (STRING) (http://string-db.org) with a high confidence score ≥ 0.7. The PPI network was built based on the interaction data between the diseases targets of IVDD and BMs. The node were scored by the innate algorithm of the STRING database. Cytoscape software (version 3.8.0) were used to visualize. Hub genes were selected by using the plug-in molecular complex detection technology (MCODE). The parameter of MCODE is as follow: K-score = 2, degree cutoff = 2, max depth = 100, and node score cutoff = 0.2.

Enrichment analysis

To further elucidate the biological functions if the hub genes, the Clue GO plug-in from the cytoscape software to conduct GO annotation and KEGG pathway analysis. The GO annotation focused on the biological processes relevant to our study. The results of the enrichment were filtered based on their P-value, with a P-value < 0.05 deemed statistically significant. The biological processes obtained were fused and clustered using ClueGo to provide a better understanding of the results. The KEGG was utilized to identify common pathways and maps, as well as to determine potential associations between key genes and pathways, and between pathways themselves. This approach can provide valuable insights for future research and aid in the development of novel therapeutic approaches for a range of diseases and conditions.

Drug prediction and molecular docking

DsigDB database (http://dsigdb.tanlab.org) was utilized to analyze the hub genes and identified potential drug candidates. To investigate the potential docking mode and binding energy between targets and drugs, molecular docking was performed by using AutoDock software (version 1.5.6). The structures of critical targets were retrieved from RCSB Protein Data Bank (http://www.rcsb.org/) in PDB format. The 3D structures of candidate drugs were obtained from PubChem (http://pubchem.ncbi.nlm.nih.gov), which were transformed by Open Babel Toolkit (version 2.4.1) into a MOL2 file format. The AutoDock software was used to optimize proteins’ structure for molecular docking. Finally, the PyMOL software (version 2.4.1) was conducted for the visualization of docking results.

Hub genes verification and correlation analysis

Microarray dataset GSE70362 was downloaded from the GEO database (http://www.ncbi.nlm.nih.gov/geo/). According to the Thompson grading of the fibrous ring tissue samples, those with Thompson grade 1–2 were defined as the control group, and those with Thompson grade 3–5 were defined as the degeneration group. There were 8 samples in the control group and 12 samples in the degeneration group. Furthermore, the expression of hub genes was verified in GSE70362. The comparison between the expression of hub genes from the dataset was conducted with the student t-test. After that, the correlation of the significant genes was analyzed by Spearman correlation analysis. P-value < 0.05 was considered significant.

Results

Intersection targets of IVDD and BMs

After eliminating duplicate genes, we obtained a total of 718 IVDD-related genes from reliable sources such as GeneCards, DisGeNet, and OMIM. Additionally, we retrieved 222 genes associated with BMs-related conditions from the Basement membraneBASE database. By intersecting these two gene sets, we identified 41 overlapping targets that were both IVDD-related and associated with BMs (Fig. 2).

Fig. 2
figure 2

The Venn diagram. IVDD: Intervertebral disc degeneration. BM: Basement membrane

PPI network and module analysis

The PPI network was constructed using STRING, consisting of 41 nodes and 381 edges, illustrating the connections between each protein (Fig. 3A). The average node degree was 18.6, indicating a significant level of interaction among these proteins, and the average local clustering coefficient was 0.741, reflecting a satisfactory level of clustering within the network. To identify key proteins within the network, we applied the MCODE plug-in and identified 25 hub genes. Subsequently, we constructed a separate PPI network focusing specifically on these 25 hub genes (Fig. 3B). This network comprised 25 nodes and 103 edges, with an average node degree of 7.63 and an average local clustering coefficient of 0.491. The increased number of edges in the network signifies a greater degree of interconnectedness and highlights the importance of the identified targets within the network.

Fig. 3
figure 3

(A) Protein and protein interaction (PPI) network according to intersection targets. (B) The 25 hub genes screened by MCODE plug-in

Enrichment analysis of hub genes

We utilized Clue GO for GO annotation and KEGG pathway analysis of the previously identified hub genes. By applying a significance threshold of P-value < 0.05, we identified the significantly enriched GO terms, as shown in Fig. 4. Our analysis revealed that these hub genes, associated with both IVDD and BMs-related conditions, were primarily involved in various biological processes. These processes included glycosaminoglycan (GAG) catabolic process, extracellular matrix (ECM) disassembly, endoderm development, collagen fibril organization, ECM organization, negative regulation of angiogenesis, and negative regulation of cell-matrix adhesion (Fig. 4A). The focus of these biological processes predominantly centered around the synthesis and metabolic functions of substances such as ECM and GAG. In terms of KEGG pathway analysis, the major pathways enriched by these hub genes included the TGF-β signaling pathway, ECM-receptor interaction, AGE-RAGE signaling pathway in diabetic complications, and bladder cancer (Fig. 4B). Notably, the TGF-β pathway, which is closely linked to ECM, and the receptor-related pathway of ECM, were prominently enriched, aligning with the results obtained from the GO analysis.

Fig. 4
figure 4

(A) Gene Ontology (GO) enrichment according to intersection targets of IVDD and BMs. The pie chart on the left represents the rate of pathways enrichment. (B) Pathway enrichment according to intersection targets of IVDD and BMs. The pie chart on the left represents the rate of enrichment. The network chart on the left represents the connection of biology processes and pathways. **P-value < 0.01

Prediction of potential drugs and molecular docking results

Using the DSigDB database for prediction, we identified the top 5 predicted drugs for preventing AF degeneration based on p-value: Chromium, Chitosamine, Healon, Retinoic acid, and CGS-27,023 A. We then considered the safety of these drugs in humans and their biological activity to select candidates for further study in the context of both MBs-related AF degeneration and IVDD. Among the selected drugs, Chitosamine and Retinoic acid emerged as promising candidates. We investigated their potential interaction with critical genes involved in both MBs-related AF degeneration and IVDD, including ACAN, DCN, FBN1, and HSPG2 (Fig. 5). Notably, our analysis revealed that Chitosamine displayed binding energy of -4.52 kcal/mol with ACAN. In comparison, Retinoic acid exhibited binding energies of -6.45 kcal/mol, -8.62 kcal/mol, and − 5.88 kcal/mol with FBN1, HSPG2, and DCN, respectively (Fig. 6). By incorporating the identification of unique drug selections associated with both MBs-related AF degeneration and IVDD, our study highlights the potential therapeutic relevance of Chitosamine and Retinoic acid. These findings warrant further investigation into their effectiveness in mitigating both conditions and provide valuable insights into potential therapeutic targets for the development of future treatments.

Fig. 5
figure 5

Sankey diagram presented the relationship with hub genes and predicted drugs. The right is the bubble diagram of gene enrichment. The gradual color represents the P-value and the increased bubble size represents the gene count

Fig. 6
figure 6

Molecular docking diagram of candidate drugs to target. (A) ACAN-chitosamine; (B) FBN1-retonoic acid; (C) HSPG2-retonoic acid; (D) DCN-retonoic acid

Hub genes verification and correlation analysis by GEO

We utilized samples from AF sources in the GSE70362 dataset to validate the hub genes. Among the 11 genes selected based on the PPI network and MCODE plug-in analysis, we found that ACAN, DCN, FBN1, and HSPG2 exhibited more significant differential expression (Fig. 7A). These hub genes demonstrated a more vital binding ability with the predicted drugs according to the previous molecular docking results. Furthermore, we conducted correlation analysis and observed that, except for FBN1 and ACAN, and DCN and ACAN, the expression of the remaining genes showed a positive correlation (Fig. 7B). This analysis provides additional evidence supporting the potential involvement of these hub genes in the context of AF degeneration and IVDD.

Fig. 7
figure 7

(A) Validation of hub genes in GSE70362. Con = healthy groups; IVDD = intervertebral disc degeneration group. ns represents no significant. *P-value < 0.05; **P-value < 0.01. (B) Correlation of hub genes with significantly different expression

Discussion

IVDD is a multifaceted disease related to the spine, and its progression cannot be reversed by current treatments. In healthy IVD tissue, the ECM provides a dynamic microenvironment that constantly undergoes remodeling, enabling it to withstand compression and torsion forces [12, 13]. BMs are a type of cell-adherent ECM composed of glycoproteins and proteoglycans, such as laminins and collagens, which are also abundant in IVD tissue, including the NP and AF [14, 15]. Despite limited research on BMs in IVDD, understanding the metabolic signature and relevant mechanisms of ECM is crucial for comprehending the cellular mechanisms of musculoskeletal degenerative diseases, including IVDD [12]. Therefore, it is reasonable to assume that the metabolic disorder of BMs and ECM may contribute to the pathogenesis of IVDD. However, the specific mechanisms of IVDD associated with BMs and ECM remain inconclusive. Consequently, using multiple databases, we identified a series of BMs-related genes to gain insight into the pathogenic mechanism of IVDD and explore potential drug candidates.

In this study, we conducted GO and pathway enrichment analyses in order to identify key biological processes and signaling pathways involved in IVDD. The essential biological processes are related to GAG catabolic, ECM disassembly, endoderm development, collagen fibril organization, and ECM organization. Of particular importance is the role of GAG catabolism, which has been shown to play a crucial role in human IVD. Proteoglycans, which contain GAGs, are degraded and depleted during IVDD, leading to the loss of ECM integrity [16]. Aberrant changes in GAG chain substitutions during degeneration can exacerbate this process [17]. To maintain normal IVD function, it is essential to maintain high concentrations and charges of proteoglycans, with functional proteoglycans possessing as many sulfated GAG chain substitutions as possible. However, enzymes like XT-1 and GlcAT-I for GAG biosynthesis are downregulated with age and degeneration, leading to diminished functional capacity of proteoglycans and abnormal GAG chain elongation [18, 19]. Additionally, our enrichment analysis revealed that the metabolism and regulation of ECM are also crucial in IVDD. As IVDD progresses, collagen type II and proteoglycans decrease notably, resulting in less resistance to loading of the ECM [20, 21]. Furthermore, the increase of fibrotic-like collagen type I contributes to the stiffness of IVD [22]. Therefore, a comprehensive understanding of the role of ECM in IVDD is essential for the development of novel treatments. Moreover, the disassembly and organization of the extracellular matrix (ECM) play a crucial role in the pathogenesis of IVDD with BMs, as revealed by the GO analysis. The ECM is a dynamic three-dimensional structure that appears in various tissues and can remodel and degenerate, contributing to the maturation and degeneration of the IVD [23, 24]. With the aging of the IVDD process, impaired ECM organization occurs simultaneously. In the worsened microenvironment, degenerative IVD tissue significantly affects the expression of ECM proteins. Razaq et al. [25] found that the synthesis of ECM proteins decreased in NP cells as the extracellular pH declined from 7.4 to 6.3. Furthermore, disassembled ECM has been detected in diseased IVDs. Patel et al. [26] demonstrated the presence of aggrecan fragmentation in human IVD tissue with disorders. The destruction of aggregating proteoglycans may be due to the failure of newly synthesized aggrecan to aggregate, the degradation of aggregating proteoglycans, and the dysfunction of link protein [27, 28].

The results of the pathway enrichment analysis indicate that the TGF-β signaling pathway is the primary pathway involved in IVDD. Previous studies have demonstrated a strong correlation between TGF-βs and TGF-β receptors [29]. Additionally, an animal study has suggested that the expression of TGF-βs and TGF-β receptors decreases in senescence-accelerated mouse models [30]. However, mRNA expression of TGF-βs and TGF-β receptors in both the NP and AF of older rabbits was found to be higher than that of younger rabbits [31]. The TGF-β signaling pathway is known to repair IVD degeneration by inhibiting ECM degradation and increasing ECM synthesis [32]. However, proinflammatory cytokines such as TNF-α and IL-1β can be secreted by immune cells and IVD cells, which promote ECM degradation by upregulating enzymes such as metalloprotease within thrombospondin motifs (ADAMTS), disintegrins, and MMPs [33]. Despite this, TGF-β has the ability to partially retard the upregulation of proinflammatory cytokines-induced matrix-degrading proteins by regulating the MAPK and NF-κB signaling pathways [34, 35]. Hu et al. [36] demonstrated that TGF-β can improve GAG biosynthesis through SMAD/3, RhoA/ROCK1, and MAPK pathways, which further highlights the importance of the TGF-β pathway in preventing IVD degeneration. In our study, we found that chitosamine and retinoic acid may have the potential to treat fibro-ring degeneration, and such an effect may be associated with the TGF-β pathway. This inspired us that future studies may need to focus on the role of drugs on this signaling pathway.

In this study, we identified four hub genes, namely ACAN, DCN, FBN1, and HSPG2, through the PPI network and validation of the GEO dataset. These hub genes exhibited significant differential expression between healthy and degenerated AF tissue. ACAN, also known as aggrecan, is a proteoglycan that creates an ion concentration imbalance between tissue and ambient fluid, resulting in an osmotic gradient for maintaining tissue hydration [37]. Previous studies have shown that the ECM of the inner AF is rich in aggrecan, which provides swelling pressure for necessary fiber elasticity and resistance against inner compressive, radical, and shear stresses [38,39,40]. Decorin (DCN) is an archetypal small leucine-rich proteoglycan (SLRP) that has been observed in lesioned AF tissue [41]. Previous studies have linked decorin to ECM damage [42] and remodeling and have detected it in damaged IVDs [34]. Zwambag et al. [43] demonstrated that decorin can bind to Toll-like receptor (TLR) 2/4 and improve pro-inflammatory chemokine and cytokine expression in AF cells. However, in the current study, decorin is differentially expressed in non-degenerative AF tissue. Therefore, future studies should investigate whether decorin expression changes with time course and mechanical properties.

Fibrillin-1 (FBN1) is a significant factor in the pathogenesis of Marfan syndrome, a genetic disorder characterized by connective tissue abnormalities [44]. In addition to its role in Marfan syndrome, Fibrillin-1 has been found to be involved in the regulation of the TGF-β signaling pathway, which controls the dispersal, storage, and release of ligands [45]. Multiple studies have demonstrated that Fibrillin-1 is widely distributed in human AF tissue and plays a crucial role in maintaining the homeostasis of the intervertebral disc bioenvironment [46, 47]. Heparin sulfate proteoglycan 2 (HSPG2) encodes for perlecan, a protein that is deposited in BMs and found in matrices in muscle, cartilage, and bone marrow. Hayes et al. [48] confirmed that Fibrillin-1 was differentially expressed in the HS-deficient HSPG mutant mouse, which highlights the role of perlecan in influencing the matrix organization.

According to predictions from the DSigDB database, chitosamine and retinoic acid are significant drugs that target hub genes, including ACAN, DCN, FBN1, and HSPG2. To investigate the interaction of these proteins with the predicted drugs, molecular docking was performed for further validation. The results of docking indicate that chitosamine can bind well with ACAN and retinoic acid can bind well with DCN, FBN1, and HSPG2. Chitosamine, also known as glucosamine, is a naturally occurring component of glycosaminoglycans found in cartilage. It has been demonstrated to have a therapeutic effect on osteoarthritis [49]. However, it is essential to conduct further research to investigate the effect of glucosamine on AF tissue in vivo, under conditions of inflammatory and mechanical stimulation, to ensure its efficacy in living organisms. Previous studies have shown that the treatment of glucosamine could decrease the expression of aggrecan in rabbit AF cells [50]. Retinoic acid, an active metabolite of vitamin A, plays a crucial role in the modulation and maintenance of the immune response [51]. Despite its potential therapeutic application in treating AF degeneration, there have been few studies investigating the pharmacological effect of retinoic acid. Therefore, additional research is necessary to explore the potential molecular mechanism of retinoic acid and its predicted targets. A mechanism diagram of the pharmacologic interaction is presented in Fig. 8.

Fig. 8
figure 8

Schematic diagram of drug-target interaction. AF: annulus fibrosus; NP: nucleus pulposus

Limitation

There are certain limitations inherent in our study: [1] Notably, not all of the hub genes identified in our analysis were present in the GSE70362 dataset, resulting in the inability to validate the expression of 14 genes [2]. While molecular docking provided valuable insights through computer simulations, it is crucial to acknowledge that this study did not experimentally verify the actual protein-drug molecule docking [3]. Further investigations are warranted in future studies to fully elucidate the specific molecular mechanisms underlying the action of drugs on AF.

Conclusion

In conclusion, our results show that the metabolism of ECM and the TGF-β signaling pathways are important in the degeneration of IVD. BMs-related genes including ACAN, DCN, FBN1, and HSPG2 are the main targets for impeding the degeneration of AF tissue. Chitosamine and retinoic acid may be the potential therapeutic drugs for treating IVDD through molecular interaction with the key targets above.

Data Availability

The GSE70362 dataset used in this study can be found in the Gene Expression Omnibus (http://www.nicb.nlm.nih.gov/geo/).

Abbreviations

IVDD:

Intervertebral disc degeneration

IVD:

Intervertebral disc

AF:

Annulus fibrosus

NP:

Nucleus pulposus

ECM:

Extracellular matrix

BMs:

Basement membranes

PPI:

Protein-to-protein interaction

GO:

Gene ontology

GEO:

Gene Expression Omnibus

DSigDB:

Drug Signature database

OMIM:

Online Mendelian Inheritance in Man

STRING:

Search Tool for the Retrieval of Interacting Genes

MCODE:

molecular complex detection technology

BP:

Biological processes

GAG:

Glycosaminoglycan

ADAMTS:

Metalloprotease within thrombospondin motifs

DCN:

decorin

FBN1:

Fibrillin-1

HSPG2:

Heparin sulfate proteoglycan 2

References

  1. Yan Q, Xiao Q, Ge J, Wu C, Wang Y, Yu H, et al. Bioinformatics-Based Research on Key genes and pathways of intervertebral disc degeneration. Cartilage. 2021;13(2suppl):582s–91s.

    PubMed  CAS  Google Scholar 

  2. Chu G, Shi C, Lin J, Wang S, Wang H, Liu T, et al. Biomechanics in Annulus Fibrosus Degeneration and Regeneration. Adv Exp Med Biol. 2018;1078:409–20.

    PubMed  CAS  Google Scholar 

  3. Jing D, Wu W, Deng X, Peng Y, Yang W, Huang D, et al. FoxO1a mediated cadmium-induced annulus fibrosus cells apoptosis contributes to intervertebral disc degeneration in smoking. J Cell Physiol. 2021;236(1):677–87.

    PubMed  CAS  Google Scholar 

  4. Chu G, Shi C, Wang H, Zhang W, Yang H, Li B. Strategies for Annulus Fibrosus Regeneration: from Biological Therapies to tissue Engineering. Front Bioeng Biotechnol. 2018;6:90.

    PubMed  PubMed Central  Google Scholar 

  5. Novais EJ, Tran VA, Johnston SN, Darris KR, Roupas AJ, Sessions GA, et al. Long-term treatment with senolytic drugs Dasatinib and Quercetin ameliorates age-dependent intervertebral disc degeneration in mice. Nat Commun. 2021;12(1):5213.

    PubMed  PubMed Central  CAS  Google Scholar 

  6. Melrose J, Smith S, Ghosh P, Whitelock J. Perlecan, the multidomain heparan sulfate proteoglycan of basement membranes, is also a prominent component of the cartilaginous primordia in the developing human fetal spine. J Histochem Cytochemistry: Official J Histochem Soc. 2003;51(10):1331–41.

    CAS  Google Scholar 

  7. Jayadev R, Morais M, Ellingford JM, Srinivasan S, Naylor RW, Lawless C, et al. A basement membrane discovery pipeline uncovers network complexity, regulators, and human disease associations. Sci Adv. 2022;8(20):eabn2265.

    PubMed  PubMed Central  CAS  Google Scholar 

  8. Melrose J, Smith S, Ghosh P, Whitelock J. Perlecan, the multidomain heparan sulfate proteoglycan of basement membranes, is also a prominent component of the cartilaginous primordia in the developing human fetal spine. J Histochem Cytochem. 2003;51(10):1331–41.

    PubMed  CAS  Google Scholar 

  9. Wang H, Liu W, Yu B, Yu X, Chen B. Identification of key modules and hub genes of Annulus Fibrosus in Intervertebral Disc Degeneration. Front Genet. 2020;11:596174.

    PubMed  CAS  Google Scholar 

  10. Xunlu Y, Minshan F, Liguo Z, Jiawen Z, Xu W, Jie Y, et al. Integrative Bioinformatics Analysis reveals potential gene biomarkers and analysis of function in human degenerative disc Annulus Fibrosus cells. Biomed Res Int. 2019;2019:9890279.

    PubMed  PubMed Central  Google Scholar 

  11. Sharifi S, Bulstra SK, Grijpma DW, Kuijer R. Treatment of the degenerated intervertebral disc; closure, repair and regeneration of the annulus fibrosus. J Tissue Eng Regen Med. 2015;9(10):1120–32.

    PubMed  Google Scholar 

  12. Sekiguchi R, Yamada KM. Basement membranes in Development and Disease. Curr Top Dev Biol. 2018;130:143–91.

    PubMed  PubMed Central  CAS  Google Scholar 

  13. Zhang S, Liu W, Chen S, Wang B, Wang P, Hu B, et al. Extracellular matrix in intervertebral disc: basic and translational implications. Cell Tissue Res. 2022;390(1):1–22.

    PubMed  Google Scholar 

  14. Kibble MJ, Domingos M, Hoyland JA, Richardson SM. Importance of Matrix Cues on intervertebral Disc Development, Degeneration, and regeneration. Int J Mol Sci. 2022;23(13).

  15. Pozzi A, Yurchenco PD, Iozzo RV. The nature and biology of basement membranes. Matrix Biology: Journal of the International Society for Matrix Biology. 2017;57–58:1–11.

    PubMed  Google Scholar 

  16. Liu X, Krishnamoorthy D, Lin L, Xue P, Zhang F, Chi L, et al. A method for characterising human intervertebral disc glycosaminoglycan disaccharides using liquid chromatography-mass spectrometry with multiple reaction monitoring. Eur Cells Mater. 2018;35:117–31.

    CAS  Google Scholar 

  17. Silagi ES, Shapiro IM, Risbud MV. Glycosaminoglycan synthesis in the nucleus pulposus: Dysregulation and the pathogenesis of disc degeneration. Matrix Biology: Journal of the International Society for Matrix Biology. 2018;71–72:368–79.

    PubMed  Google Scholar 

  18. Ye W, Zhou J, Markova DZ, Tian Y, Li J, Anderson DG, et al. Xylosyltransferase-1 expression is refractory to inhibition by the inflammatory cytokines tumor necrosis factor α and IL-1β in nucleus pulposus cells: novel regulation by AP-1, Sp1, and Sp3. Am J Pathol. 2015;185(2):485–95.

    PubMed  PubMed Central  CAS  Google Scholar 

  19. Collin EC, Carroll O, Kilcoyne M, Peroglio M, See E, Hendig D, et al. Ageing affects chondroitin sulfates and their synthetic enzymes in the intervertebral disc. Signal Transduct Target Therapy. 2017;2:17049.

    Google Scholar 

  20. Weiler C, Schietzsch M, Kirchner T, Nerlich AG, Boos N, Wuertz K. Age-related changes in human cervical, thoracal and lumbar intervertebral disc exhibit a strong intra-individual correlation. Eur Spine J. 2012;21(Suppl 6):810–8.

    Google Scholar 

  21. Antoniou J, Steffen T, Nelson F, Winterbottom N, Hollander AP, Poole RA, et al. The human lumbar intervertebral disc: evidence for changes in the biosynthesis and denaturation of the extracellular matrix with growth, maturation, ageing, and degeneration. J Clin Investig. 1996;98(4):996–1003.

    PubMed  PubMed Central  CAS  Google Scholar 

  22. Iatridis JC, Setton LA, Weidenbaum M, Mow VC. Alterations in the mechanical behavior of the human lumbar nucleus pulposus with degeneration and aging. J Orthop Res. 1997;15(2):318–22.

    PubMed  CAS  Google Scholar 

  23. Zhang S, Liu W, Chen S, et al. Extracellular matrix in intervertebral disc: basic and translational implications. Cell Tissue Res. 2022;390(1):1–22.

    PubMed  Google Scholar 

  24. Peng Z, Sun H, Bunpetch V, et al. The regulation of cartilage extracellular matrix homeostasis in joint cartilage degeneration and regeneration. Biomaterials. 2021;268:120555.

    PubMed  CAS  Google Scholar 

  25. Razaq S, Wilkins RJ, Urban JP. The effect of extracellular pH on matrix turnover by cells of the bovine nucleus pulposus. Eur Spine J. 2003;12(4):341–9.

    PubMed  PubMed Central  Google Scholar 

  26. Patel KP, Sandy JD, Akeda K, et al. Aggrecanases and aggrecanase-generated fragments in the human intervertebral disc at early and advanced stages of disc degeneration. Spine (Phila Pa 1976). 2007;32(23):2596–603.

    PubMed  Google Scholar 

  27. Pearce RH, Mathieson JM, Mort JS, Roughley PJ. Effect of age on the abundance and fragmentation of link protein of the human intervertebral disc. J Orthop Res. 1989;7(6):861–7.

    PubMed  CAS  Google Scholar 

  28. Sivan SS, Wachtel E, Roughley P. Structure, function, aging and turnover of aggrecan in the intervertebral disc. Biochim Biophys Acta. 2014;1840(10):3181–9.

    PubMed  CAS  Google Scholar 

  29. Chen S, Liu S, Ma K, Zhao L, Lin H, Shao Z. TGF-β signaling in intervertebral disc health and disease. Osteoarthr Cartil. 2019;27(8):1109–17.

    CAS  Google Scholar 

  30. Nagano S, Matsunaga S, Takae R, Morimoto N, Suzuki S, Yoshida H. Immunolocalization of transforming growth factor-betas and their receptors in the intervertebral disk of senescence-accelerated mouse. Int J Oncol. 2000;17(3):461–6.

    PubMed  CAS  Google Scholar 

  31. Murakami H, Yoon ST, Attallah-Wasif ES, Tsai KJ, Fei Q, Hutton WC. The expression of anabolic cytokines in intervertebral discs in age-related degeneration. Spine (Phila Pa 1976). 2006;31(16):1770–4.

    PubMed  Google Scholar 

  32. Risbud MV, Shapiro IM. Role of cytokines in intervertebral disc degeneration: pain and disc content. Nat Rev Rheumatol. 2014;10(1):44–56.

    PubMed  CAS  Google Scholar 

  33. Yang H, Liu H, Li X, Pan H, Li Z, Wang J, et al. TNF-α and TGF-β1 regulate Syndecan-4 expression in nucleus pulposus cells: role of the mitogen-activated protein kinase and NF-κB pathways. Connect Tissue Res. 2015;56(4):281–7.

    PubMed  CAS  Google Scholar 

  34. Yang H, Gao F, Li X, Wang J, Liu H, Zheng Z. TGF-β1 antagonizes TNF-α induced up-regulation of matrix metalloproteinase 3 in nucleus pulposus cells: role of the ERK1/2 pathway. Connect Tissue Res. 2015;56(6):461–8.

    PubMed  Google Scholar 

  35. Hu B, Xu C, Cao P, Tian Y, Zhang Y, Shi C, et al. TGF-β stimulates expression of Chondroitin Polymerizing factor in Nucleus Pulposus cells through the Smad3, RhoA/ROCK1, and MAPK signaling pathways. J Cell Biochem. 2018;119(1):566–79.

    PubMed  CAS  Google Scholar 

  36. Bhattacharya S, Dubey DK. Effect of aggrecan degradation on the nanomechanics of hyaluronan in extra-fibrillar matrix of annulus fibrosus: a molecular dynamics investigation. J Mech Behav Biomed Mater. 2020;107:103752.

    PubMed  CAS  Google Scholar 

  37. Iu J, Santerre JP, Kandel RA. Inner and outer annulus fibrosus cells exhibit differentiated phenotypes and yield changes in extracellular matrix protein composition in vitro on a polycarbonate urethane scaffold. Tissue Eng Part A. 2014;20(23–24):3261–9.

    PubMed  CAS  Google Scholar 

  38. Han WM, Nerurkar NL, Smith LJ, Jacobs NT, Mauck RL, Elliott DM. Multi-scale structural and tensile mechanical response of annulus fibrosus to osmotic loading. Ann Biomed Eng. 2012;40(7):1610–21.

    PubMed  PubMed Central  Google Scholar 

  39. Ayturk UM, Garcia JJ, Puttlitz CM. The micromechanical role of the annulus fibrosus components under physiological loading of the lumbar spine. J Biomech Eng. 2010;132(6):061007.

    PubMed  Google Scholar 

  40. Melrose J, Ghosh P, Taylor TK, Vernon-Roberts B, Latham J, Moore R. Elevated synthesis of biglycan and decorin in an ovine annular lesion model of experimental disc degeneration. Eur Spine J. 1997;6(6):376–84.

    PubMed  PubMed Central  CAS  Google Scholar 

  41. Melrose J, Smith SM, Fuller ES, Young AA, Roughley PJ, Dart A, et al. Biglycan and fibromodulin fragmentation correlates with temporal and spatial annular remodelling in experimentally injured ovine intervertebral discs. Eur Spine J. 2007;16(12):2193–205.

    PubMed  PubMed Central  Google Scholar 

  42. Brown S, Melrose J, Caterson B, Roughley P, Eisenstein SM, Roberts S. A comparative evaluation of the small leucine-rich proteoglycans of pathological human intervertebral discs. Eur Spine J. 2012;21(Suppl 2):154–9.

    PubMed Central  Google Scholar 

  43. Zwambag DP, Molladavoodi S, Guerreiro MJ, DeWitte-Orr SJ, Gregory DE. Immuno-stimulatory capacity of decorin in the rat tail intervertebral disc and the mechanical consequence of resultant inflammation. Eur Spine J. 2020;29(7):1641–8.

    PubMed  Google Scholar 

  44. Asano K, Cantalupo A, Sedes L, Ramirez F. The multiple functions of Fibrillin-1 Microfibrils in Organismal Physiology. Int J Mol Sci. 2022;23(3).

  45. Robertson IB, Rifkin DB. Regulation of the bioavailability of TGF-β and TGF-β-Related proteins. Cold Spring Harb Perspect Biol. 2016;8(6).

  46. Hayes AJ, Gibson MA, Shu C, Melrose J. Confocal microscopy demonstrates association of LTBP-2 in fibrillin-1 microfibrils and colocalisation with perlecan in the disc cell pericellular matrix. Tissue Cell. 2014;46(3):185–97.

    PubMed  CAS  Google Scholar 

  47. Hayes AJ, Smith SM, Gibson MA, Melrose J. Comparative immunolocalization of the elastin fiber-associated proteins fibrillin-1, LTBP-2, and MAGP-1 with components of the collagenous and proteoglycan matrix of the fetal human intervertebral disc. Spine (Phila Pa 1976). 2011;36(21):E1365–72.

    PubMed  Google Scholar 

  48. Hayes AJ, Smith SM, Melrose J. Comparative immunolocalisation of fibrillin-1 and perlecan in the human foetal, and HS-deficient hspg2 exon 3 null mutant mouse intervertebral disc. Histochem Cell Biol. 2013;139(1):1–11.

    PubMed  CAS  Google Scholar 

  49. Towheed TE, Maxwell L, Anastassiades TP, Shea B, Houpt J, Robinson V, et al. Glucosamine therapy for treating osteoarthritis. Cochrane Database Syst Rev. 2005;2005(2):Cd002946.

    PubMed  PubMed Central  Google Scholar 

  50. Sowa GA, Coelho JP, Jacobs LJ, Komperda K, Sherry N, Vo NV, et al. The effects of glucosamine sulfate on intervertebral disc annulus fibrosus cells in vitro. The Spine Journal: Official Journal of the North American Spine Society. 2015;15(6):1339–46.

    PubMed  Google Scholar 

  51. Abdelhamid L, Luo XM. Retinoic acid, leaky gut, and Autoimmune Diseases. Nutrients. 2018;10(8).

Download references

Acknowledgements

We appreciate GEO database for providing their platforms and thank for those contributors for uploading their meaningful datasets. Special thank for Ms. Zi Shen for drawing the schematic diagram.

Funding

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Author information

Authors and Affiliations

Authors

Contributions

Z.Z, X.J, H.R and Y.S designed the research. Z.Z wrote the original article; G.Z, W.Q and P.Z collected data. J.H and Z.C provided the illustrations. Y.G, D.L provided the analytic methods.

Corresponding author

Correspondence to Yuping Sun.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhou, Z., Qin, W., Zhang, P. et al. Potential molecular targets and drugs for basement membranes-related intervertebral disk degeneration through bioinformatics analysis and molecular docking. BMC Musculoskelet Disord 24, 772 (2023). https://doi.org/10.1186/s12891-023-06891-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12891-023-06891-z

Keywords