Next Article in Journal
The Effect of Ferula communis Extract in Escherichia coli Lipopolysaccharide-Induced Neuroinflammation in Cultured Neurons and Oligodendrocytes
Next Article in Special Issue
Osteoblast Dysfunction in Non-Hereditary Sclerosing Bone Diseases
Previous Article in Journal
The Rab Geranylgeranyl Transferase Beta Subunit Is Essential for Embryo and Seed Development in Arabidopsis thaliana
Previous Article in Special Issue
Expression and Role of Ubiquitin-Specific Peptidases in Osteoblasts
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 15
10.3390/ijms22157908
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Proteomic and Metabolomic Analysis of Human Osteoblasts, Differentiated from Dental Pulp Stem Cells, Hinted Crucial Signaling Pathways Promoting Osteogenesis

by
Slavomíra Nováková
1,*,†,
Maksym Danchenko
2,†,
Terézia Okajčeková
1,
Eva Baranovičová
1,
Andrej Kováč
3,
Marián Grendár
1,
Gábor Beke
4,
Janka Pálešová
1,
Ján Strnádel
1,
Mária Janíčková
5,
Erika Halašová
1,6 and
Henrieta Škovierová
1,*
1
Biomedical Centre Martin, Jessenius Faculty of Medicine in Martin, Comenius University in Bratislava (JFM CU), Malá Hora 4C, 036 01 Martin, Slovakia
2
Plant Science and Biodiversity Center, Slovak Academy of Sciences, Dúbravská cesta 9, 845 23 Bratislava, Slovakia
3
Institute of Neuroimmunology, Slovak Academy of Sciences, Dúbravská cesta 9, 845 10 Bratislava, Slovakia
4
Institute of Molecular Biology, Slovak Academy of Sciences, Dúbravská cesta 21, 845 51 Bratislava, Slovakia
5
Department of Stomatology and Maxillofacial Surgery, University Hospital in Martin and JFM CU, Kollárova 2, 036 01 Martin, Slovakia
6
Department of Medical Biology, Jessenius Faculty of Medicine in Martin, Comenius University in Bratislava (JFM CU), Malá Hora 4C, 036 01 Martin, Slovakia
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2021, 22(15), 7908; https://doi.org/10.3390/ijms22157908
Submission received: 24 May 2021 / Revised: 18 July 2021 / Accepted: 20 July 2021 / Published: 24 July 2021
(This article belongs to the Special Issue Osteoblast Differentiation and Activity in Skeletal Diseases)
Browse Figures
Review Reports Versions Notes

Abstract

:
Population aging has been a global trend for the last decades, which increases the pressure to develop new cell-based or drug-based therapies, including those that may cure bone diseases. To understand molecular processes that underlie bone development and turnover, we followed osteogenic differentiation of human dental pulp stem cells (DPSCs) using a specific induction medium. The differentiation process imitating in vivo osteogenesis is triggered by various signaling pathways and is associated with massive proteome and metabolome changes. Proteome was profiled by ultrahigh-performance liquid chromatography and comprehensively quantified by ion mobility-enhanced mass spectrometry. From 2667 reproducibly quantified and identified proteins, 432 were differentially abundant by strict statistic criteria. Metabolome profiling was carried out by nuclear magnetic resonance. From 27 detected metabolites, 8 were differentially accumulated. KEGG and MetaboAnalyst hinted metabolic pathways that may be involved in the osteogenic process. Enrichment analysis of differentially abundant proteins highlighted PPAR, FoxO, JAK-STAT, IL-17 signaling pathways, biosynthesis of thyroid hormones and steroids, mineral absorption, and fatty acid metabolism as processes with prominent impact on osteoinduction. In parallel, metabolomic data showed that aminoacyl-tRNA biosynthesis, as well as specific amino acids, likely promote osteodifferentiation. Targeted immunoassays validated and complemented omic results. Our data underlined the complexity of the osteogenic mechanism. Finally, we proposed promising targets for future validation in patient samples, a step toward the treatment of bone defects.

Graphical Abstract

1. Introduction

The bone is a highly specific tissue that plays an essential role in the human body. It consists of specialized cells (osteoblasts, osteocytes, and osteoclasts), a mineralized and non-mineralized matrix, spaces including the vascular canals, bone marrow cavity, ducts, and gaps [1]. Bones form a skeleton that provides protection, mechanical support for soft tissues, and allows mobility. In addition, it is a mineral reservoir that supports calcium and phosphate homeostasis; therefore, the skeleton is closely linked to metabolic homeostasis as well as to the dynamic regulation by hormones and nutrients [2]. The main factor that affects the amount of bone mass is genetics. However, other variable factors such as diet, environmental involvement, exercises, current diseases, and drugs could influence bone mass [3]. Skeleton is an active organ: during growth and development, two bone-formation processes (modeling and remodeling) are involved in its sculpturing. Modeling [1] takes place during childhood and adolescence [2]. During this process, the bone mass is removed from one location and deposited to another location within the same bone for creating proper shape and size. After maturation, bone regeneration relies on a dynamic balance between bone formation (ossification) and resorption, mediated by osteoblasts and osteoclasts, respectively. This process is called remodeling [4]. The imbalance in such regulation leads to many bone diseases; for example, osteopetrosis due to excessive bone formation [5] or, in contrast, postmenopausal osteoporosis due to massive bone resorption [6].
One of the major recent challenges of regenerative medicine is bone regeneration/repair after various pathological problems such as fracture, metabolic diseases, osteoporosis, or tumors. In clinical practice, procedures for the treatment of these defects include drug application [7,8,9], using bone grafts, or the application of stem cells [10,11]. The use of medications for treating osteoporosis and fractures is limited by their efficacy and variable toxicities [3,12,13]. Therefore, clinicians recommend focusing on the prevention of fractures and defects by adopting a healthy lifestyle, using adequate nutritional supplements (calcium and vitamin D), avoiding smoking and alcohol consumption, and controlling weight [14].
The modern alternative to drug applications is tissue engineering using stem cells. Adult stem cells (in particular mesenchymal stem cells, MSCs) are very popular due to their self-renewal, differentiation capacity, multipotentiality, immune-modulatory and anti-inflammatory properties, and lack of ethical restrictions [15]. The MSCs can be directly injected into the injured bone, systemically infused, or loaded on three-dimensional scaffolds before implantation [16]. The latest option is preferable because it supports the growth of various cell types and this complex will trigger the healing capacity of bone tissue [17,18]. Human dental pulp stem cells (DPSCs), oral dental-derived stem cells, are part of the dental pulp, which ensures the formation, nutrition, protection, and restoration of teeth. Dental pulp produces dentin and maintains its physiological and biological vitality [19]. As DPSCs have the properties of MSCs, aside from their odontogenic potential, they can differentiate in vitro into osteoblasts, chondrocytes, myocytes, cardiomyocytes, active neurons, and hepatocyte-like cells, thus promoting healing of tissues [20].
To develop and practically implement the new advanced cell-based or drug-based strategies, it is necessary to understand molecular mechanisms responsible for natural osteogenesis and bone regeneration. Bone healing is a complex process including inflammation, repair, and remodeling. Various cell types (in particular MSCs) and multiple extracellular and intracellular signaling networks are involved in this well-organized physiological process [16,17]. Osteogenic microenvironment is often imitated in vitro by specific induction media [21]. The exact molecular mechanism by which cells differentiate and form a mineralized extracellular matrix is not yet fully understood.
Of note, the differences during in vitro differentiation of DPSCs into osteoblasts or odontoblasts are minor; consequently, inducing media are similar. While some researchers modulated in vitro odontogenesis by β-glycerophosphate [22], others used the medium with this supplement as osteogenic [23,24]. Specific successful differentiation of DPSCs with the proper structural arrangement into osteoblasts/odontoblasts is ensured only in vivo [25]. In the recent review, the authors argued that bone and dentin share similar characteristics in their composition, mineralization, biomarkers, and clinical manifestation of mineralization defects. It is difficult to discriminate bone/tooth differentiation markers because these proteins are expressed in both tissues [26]. Analogously, signaling pathways involved in osteogenic/odontogenic differentiation overlap [17,27,28]. Moreover, a large share of information about dentin mineralization was obtained from in vitro and in vivo bone studies [29]. Therefore, DPSCs are applicable not only in cranio-maxillofacial hard-tissue repair [30,31] but also have the potential for the treatment of various bone defects [32]. For example, dentin was applied as a graft to repair calvarial bone defect [33]. In this study, we focused on in vitro differentiation of DPSCs from the osteogenic perspective.
Omics technologies give the possibility of a complex view of the biological issue. Herein, we analyzed proteome and metabolome changes in DPSCs during osteogenic differentiation to understand the molecular background of the bone regeneration process. Comprehensive data were obtained using label-free mass spectrometry and NMR quantification. Next, we validated findings by the detection of specific proteins using Western blot analysis and immunocytochemistry. A complex understanding of signaling pathways that regulate the osteodifferentiation of stem cells advanced knowledge of the molecular mechanisms of bone development as well as offered ideas for more effective therapies for bone diseases.

2. Results and Discussion

The entire experimental workflow is shown in Figure 1 to illustrate the specific design of this work.
Bone remodeling is a complex process controlled by multiple local and systemic factors. In order to study the mechanisms that regulate osteodifferentiation, DPSCs were grown in OsteoMAX-XF differentiation medium selected in our previous study [21]. The efficiency of differentiation was demonstrated (i) microscopically, (ii) spectrophotometrically, and (iii) immunocytochemically.

2.1. DPSCs Possess In Vitro Osteogenic Potential

The differentiation process was monitored for 24 days based on cell morphology. Data in the literature point to cellular morphology correlations with the osteogenic differentiation potential as well as biochemical markers and provides strong predictive evidence in regenerative medicine [34]. During our experiment, cells showed the phenotype which was typical for a particular day of osteogenic development. This pattern was reproducible, as shown by three independent biological replicates (Figure 2 and Figure S1).
Osteogenic differentiation consists of several phases. Our work focused on the 24th day after the specific medium application. It falls into the mineralization phase, which is characterized by the production of osteocalcin (OSTCN) and osteopontin secreted by late osteoblasts as well as calcium and phosphate deposition [35]. On day 24, to confirm the proper differentiation process before performing omic analysis, OSTCN, bone sialoprotein 2 (BSP II), and bone morphogenetic protein 2 (BMP-2) as typical osteogenic markers [36] were imaged by immunocytochemistry. All three proteins showed high accumulation in osteodifferentiated cells compared to the weak signal in control (Figure 3).
Additionally, we followed the dynamics of calcium nodule production during the experimental time course. We revealed a remarkable production of mineral deposits of more than 100 times higher in differentiated cells compared to DPSCs on the 11th day after application of the osteogenic medium. Accumulation of deposits further increased almost three times from the 11th to the 24th day (from 139:1 to 396:1) (Figure 4).

2.2. Overview of Quantitative Proteomic and Metabolomic Profiling of DPSCs upon In Vitro Osteogenic Differentiation

To analyze the mechanisms of the differentiation process, a label-free quantitative approach was employed to acquire the peptide and fragment mass spectra. Dedicated software Progenesis QI with advanced chromatographic alignment, missing values free peak detection, and relative protein abundance determination, relying on the integration of peak areas of unique peptides, was used to process the mass spectrometry data. We identified and quantified 2667 proteins in the control and osteodifferentiated samples (Table S1). Among them, 432 proteins were differentially regulated based on criteria of statistical significance (p ≤ 0.01) and at least 1-fold of log2 transformed ratio (Table S2). Principal component analysis (PCA) revealed proper clustering of experimental replicates in two-dimensional space, increasing confidence in findings (Figure S2A). Then, proteins that met these stringent filtering criteria were categorized: (i) based on a biological process and (ii) based on signaling pathway. Proteins were clustered into twelve functional categories using GOTermFinder. The biological process categories affected during osteogenesis included cellular process (GO:0009987, 392 proteins), localization (GO:0051179, 223 proteins), response to stimulus (GO:0050896, 196 proteins), biological regulation (GO:0065007, 194 proteins), development (GO:0032502, 190 proteins), transport (GO:0006810, 182 proteins), cellular component organization (GO:0071840, 182 proteins), metabolism (GO:0008152, 166 proteins), cell death (GO:0008219, 89 proteins), secretion (GO:0046903, 77 proteins), immune system process (GO:0002376, 42 proteins), and cellular respiration (GO:0045333, 19 proteins). Fine granularities of ontologies are presented in Table S3. To enhance the biological context, differentially abundant proteins were also grouped into 265 metabolic and signaling pathways using KEGG (Table S4). The identified molecular pathways after filtering included several potentially crucial cascades promoting osteogenic differentiation. We described them in detail in the following section.
In parallel, we quantified 27 metabolites in the control and osteodifferentiated samples (Table S5). Among them, only eight were differentially regulated based on criteria of statistical significance (adjusted p ≤ 0.05) and at least 1-fold of log2 transformed ratio (Table S6). Then, metabolites satisfying these stringent filtering criteria were categorized into 14 metabolic pathways using MetaboAnalyst (Table S7). As in the proteins case, the PCA revealed proper clustering of three replicates in two-dimensional space (Figure S2B).

2.3. Major Signaling Pathways Involved in Osteogenesis Were Well Represented among Differentially Abundant Proteins

Initially, we focused on identifying and selecting several known signaling pathways highlighted in the literature (Table 1). We found several of the major pathways that were already described in connection with osteodifferentiation [17,37]. Namely, Ca2+, Wnt, TGF-β, PI3K-Akt, JAK-STAT, and MAPK signaling pathways were identified as affected in the mineralization phase after osteogenic medium induction on the 24th day (Figure 5, Table 1). Apart from the PI3K-Akt pathway (with a roughly equal number of proteins accumulated in control and osteodifferentiated cells), JAK-STAT and MAPK signaling cascades were both downregulated in differentiated cells. The remaining pathways were upregulated in differentiated cells, confirming their pivotal roles in osteogenesis. However, these pathways were not significantly enriched apart from JAK-STAT. Of note, the cascades crosstalk, hence form a signaling network. Differentially abundant proteins from these pathways may become either diagnostic markers in patients’ samples or even therapeutic targets.
JAK-STAT signaling pathway. MAPK and JAK-STAT signaling cascades were apparently downregulated in osteodifferentiated cells. JAK-STAT signaling is involved in many pathologies as well as bone homeostasis. STAT family consists of several members, particularly STAT1 and STAT3, that variously influence bone development [38]. Interestingly, it was shown that STAT1 absence led to the acceleration of bone healing [39]. Furthermore, STAT1/3 inhibition increased osteoblast activity, associated with OSTCN and alkaline phosphatase production [40]. STAT3 also controls osteoclast differentiation, bone degeneration, and homeostasis [41]. Our data are in line with those reports. We identified four members of JAK-STAT pathway (Table 1), signal transducer and activator of transcription 1 (STAT1), signal transducer and activator of transcription 3 (STAT3), four and half LIM domains 1 (FHL1), and glial fibrillary acidic protein (GFAP) as lower abundant in differentiated cells (−1.55, −1.11, −1.65, and −1.08-fold decrease, respectively). We hypothesize that in the mineralization phase, downregulation of STAT1 and STAT3 positively affects osteogenesis.
Ca2+ signaling pathway. Calcium and phosphate, major units of mineral metabolism, are closely linked to intermediary metabolism. Intracellular calcium also works as a secondary messenger that mediates the effect of membrane signals on the secretion of several mediators (e.g., insulin and epinephrine). Circulating calcium levels are tightly controlled through the mineral hormones PTH, 1,25-dihydroxyvitamin D, and calcitonin [42]. A higher level of cytoplasmic Ca2+ promotes phospholipase C activity and inositol-1,4,5-trisphosphate (IP3) signaling. This signaling leads to a release of Ca2+ from intracellular stores (mainly from the endoplasmic reticulum). Then, increased cytoplasmic Ca2+ further promotes regulation of several factors, such as IGF1 and IGF2, VEGF, TGF-β, BMP-2, and BMP-4, involved in other osteogenic signaling pathways [17]. We identified nine members of Ca2+ signaling pathway (Table 1). Specifically, voltage-dependent anion channels 1, 2, 3 (VDAC1, VDAC2, VDAC3; 1.18, 1.14, and 1.36-fold increase, respectively), solute carrier family 25 members 4 and 6 (SLC25-4, SLC25-6 both 1.05-fold increase), IP3 receptor type 2 (4.92-fold increase), aspartate beta-hydroxylase (ASPH, 1.07-fold increase), and G protein alpha subunits 11 and 14 (GNA11 and GNA14, −3.15 and −1.99-fold decrease, respectively). The identification of VDACs, ion channels located in the outer mitochondrial membrane, and SLC25s, solute carriers responsible for mitochondrial ATP transport, points to the role of mitochondria in the Ca2+ influx and efflux as more important than the endoplasmic reticulum during osteodifferentiation. This interpretation corresponds with the conclusions of Alves et al. [43]. Our results were strengthened by determining the IP3 receptor (validation of proteomic data) and BMP-2 (complementation) by Western blotting. Densitometry showed a statistically significant 4.32 and 7.65-fold increase, both p = 0.044 in osteodifferentiated cells compared to controls (Figure 6). A better understanding of the Ca2+-sensing mechanism may lead to novel targets for promoting osteogenesis because of close links with other cascades. For instance, it was reported that increased Ca2+ through its channels enhanced bone formation and prevented estrogen deficiency-induced bone loss [44].
Wnt signaling pathway. This pathway plays an important role during the growth, development, and homeostasis of numerous organs, including the bone system [45]. Wnts belong to a family of ligands for membrane-spanning frizzled (Fzd) receptors. Wnts can activate at least three signaling subcascades: Wnt/Ca2+ (noncanonical), Wnt/β-catenin (canonical), and Wnt/planar polarity (noncanonical) [17]. We detected a single Wnt protein (Table 1), Wnt family member 5A (1.24-fold increase in osteoblasts). Considering upregulated Ca2+ signaling, we suggest activation of noncanonical subcascade. It is initiated by Wnt ligands binding to Fzd receptors followed by the activation of G protein and downstream events releasing intracellular calcium from the endoplasmic reticulum [46]. Of note, it was reported that Wnt5A-induced calcium-dependent pathway promoted osteoblast differentiation through the upregulation of Lrp5 and Lrp6 and enhanced Wnt/β-catenin signaling [47]. Additionally, other studies pointed to the regulation of the Wnt signaling during osteogenesis; thus, it became an attractive therapeutic target for the treatment of bone defects [46].
TGF-β cascade. BMP/TGF-β signaling is a crucial network associated with maintaining bone development by regulating the balance between bone resorption and bone formation. Signaling takes place through the canonical (Smad-dependent) or noncanonical (Smad-independent) pathway. BMPs and TGF-βs are cytokines belonging to TGF-β superfamily [48]. It includes several subgroups, among others activins and inhibins. We identified only two members of TGF-β pathway (Table 1), decorin and inhibin subunit beta A, both accumulated in osteodifferentiated cells (1.75 and 5.20-fold increase, respectively). There is only scarce information in the literature about the modulation of osteogenesis by activins and inhibins, such as they act as modulators of TGF-β and BMP. Particularly, Gaddy-Kurten et al. reported that inhibins A and βA repress osteoblast differentiation and reduce alkaline phosphatase activity [49]. Decorin is a small leucine-rich proteoglycan associated with extracellular matrix deposition in the mineralization stage of bone and dentin development. It preferentially binds hydroxyapatite, regulating crystal size and morphology [50]. Covalent binding of decorin to TGF-β resulted in increased bone mass associated with a decreased level of TGF-β. Active TGF-β promotes osteoclast differentiation at a low dose but inhibits this process at a high dose [51]. Thus, we argue that a higher abundance of inhibin and decorin indicates the termination of active differentiation of osteoblasts, which reached the mineralization phase. Additionally, decorin is a potential marker or target to improve mineralization in bones. For instance, the authors showed low expression of decorin in the skin of patients with osteoporosis imperfecta [52]. A study on pigs demonstrated the importance of dentin sialophosphoprotein (DSPP) in maintaining TGF-β1 activity in dentin [53]. We detected (by complementary Western blotting) accumulation of its cleavage product (likely dentin sialoprotein, DSP) at 37 kDa in osteodifferentiated cells (14.91-fold increase, p = 0.0002) (Figure 6). The size of the observed DSP band is the same as reported by Yamakoshi et al. [54] as well as in DPSCs study of Teti et al. [55]. Our interpretation is that accumulation of decorin led to inactivation of TGF-β1, thus higher abundance of DSP in osteoblasts. Although the DSPP (also DSP) protein is primarily associated with odontogenesis, it was identified also in bones with function in hard tissue development. Verdelis et al. demonstrated DSPP expression at low level in the long bone; moreover, its disruption in mice resulted in modest but significant negative changes [56]. Jani et al. suggested that DSPP and DMP1 may function synergistically. Specifically, they showed that Dspp rescued the long bone defects in Dmp1-deficient mice [57].
MAPK signaling pathway. Mitogen-activated protein kinase (MAPK) cascade shares a few members with non-canonical TGF-β signaling. Seven from ten identified members were less abundant in osteodifferentiated cells (Table 1). Only insulin (INS), insulin-like growth factor 2 (IGF-2), and CD14 protein accumulated in osteoblasts (2.64, 1.12, and 2.42-fold increase, respectively). INS and IGF-2 are hormonal signaling molecules that affect a broad spectrum of metabolic and developmental processes. They regulate signaling from extracellular space until target expression in the cell. Several studies reported their impact on osteoblasts and bone development through PI3K-Akt and MAPK pathways [58,59]. Lower abundant proteins in osteoblasts included MAP3K20, MAP4K4, microtubule-associated protein tau, protein phosphatase 5, calcium voltage-gated channel, CRK proto-oncogene, and G protein subunit alpha 12 (−1.32, −3.12, −3.43, −1.17, −1.31, −1.08, and −1.73-fold decrease, respectively). The role of the MAPK pathway in osteoblasts remains controversial, as some studies have reported its stimulatory effect while others reported an inhibitory one [60]. Interestingly, the authors reported a direct connection of dentin matrix protein (DMP-1) to MAPK signaling activation and maintenance of osteocytes’ phenotype [61]. In contrast, another work suggested the involvement of DMP-1 as a negative regulator of bone formation [62]. Our data are consistent with the latter report. Complementary Western blotting showed a lower abundance of DMP-1 (−5.54-fold decrease, p = 0.18) in differentiated cells (Figure 6).
PI3K-Akt pathway. Phosphatidylinositol 3-kinase (PI3K) and its downstream targets, serine/threonine kinases (AKT), play an essential role in the differentiation of progenitor cells into mature osteoblasts [63]. We revealed nine members of this pathway (Table 1): four accumulated (insulin, insulin-like growth factor 2, heat shock protein 90, and laminin subunit beta 1) and five less abundant (collagen types 1/4/6 alpha 1 chain, laminin subunit beta 2, and reelin) in differentiated cells. Of note, several collagen types were reported as factors that promote osteodifferentiation [64,65]. We believe that a lower abundance of collagens in mature osteoblasts, together with accumulation of Hsp 90 (1.99-fold increase), a factor that induces osteoclast-associated genes and enhances osteoclastogenesis, indicates increased bone turnover. On the other hand, undifferentiated DPSCs (control), synthesize and secrete massive amounts of collagen, particularly types 1 and 6 [66]. Therefore, the lower level of collagen isoforms in differentiated cells is not surprising. Accumulation of Col6A1 in controls was validated by Western blotting (6.33-fold increase, p = 0.0002) (Figure 6).
Overall, our data on major signaling pathways connected with osteodifferentiation highlighted the following conclusions: (i) in the mineralization phase, downregulation of STAT1 and STAT3 positively affects osteogenesis; (ii) the accumulation of VDACs and SLC25s points to the primary role of mitochondria in the Ca2+ influx/efflux during osteodifferentiation; (iii) Wnt5a is an attractive therapeutic target for the treatment of bone defects; (iv) accumulation of decorin inactivated TGF-β, leading to accumulation of DSP in osteoblasts; (v) DMP-1 is a negative regulator of bone formation via MAPK cascade; and (vi) lower abundance of collagens belonging to PI3K-Akt cascade in mature osteoblasts indicates increased bone turnover.

2.4. Selected Pathways Potentially Associated with Bone Metabolism

Despite confirmation of several signaling pathways and molecules involved in osteogenesis, much remains unknown. We identified proteins in several cascades, which showed high enrichment scores in our analysis (Table 1 and Table S8, Figure 5). They are not considered major, since little is known about their connection to bone metabolism. Therefore, they need further exploration.

2.4.1. Prevention of Oxidation Stress by FoxO Cascade, Hormone Synthesis, and Mineral Absorption Are Factors That May Affect the Mineralization Phase of Osteodifferentiation

FoxO signaling pathway and pathway for the proteasome. Forkhead box O (FOXO) family transcription factors have an essential function in numerous biological processes, including cellular proliferation, differentiation, and apoptosis [67]. We identified four members of FoxO pathway (Table 1), protein arginine methyltransferase (PRMT1), superoxide dismutase (SOD2), INS more abundant (2.31, 3.17, 2.64-fold increase, respectively) and STAT3 less abundant (−1.11-fold decrease) in differentiated cells. FoxO signaling deficiency usually leads to bone mass reduction. Moreover, FOXOs maintain the lifespan of mature osteoblasts through antioxidants [68]. During differentiation ROS levels frequently increase due to enhanced cellular metabolism [69]. FOXO promotes the expression of antioxidant enzymes, which maintain low ROS levels, thus protecting the cells from osteoclastogenesis [70]. In line with these reports, we detected the accumulation of enzymes, antioxidant SOD2, and regulatory PRMT1, in differentiated cells. The authors reported that the latter methylated FOXO1 at conserved Arg248 and Arg250. This methylation directly blocked Akt-mediated phosphorylation of FOXO, which leads to translocation to the cytoplasm and subsequent degradation via the ubiquitin-proteasome system. Silencing of PRMT1 usually enhances proteasomal degradation of FOXO1 [71,72]. Methylation blocks the exclusion of FOXO1 from the nucleus, and consequently activates the expression of target genes, particularly, antioxidant enzymes [73]. This link is additionally supported by the identification of pathway for the proteasome downregulated in differentiated cells (Table 1). In addition, a study on the animal model showed that resveratrol prevented osteoporosis by increasing FOXO1 transcription [74].
IL-17 pathway. Interleukin 17 (IL-17), produced by Th17 cells, is a family of pro-inflammatory cytokines having a role in bone remodeling. They enhance MSCF and RANKL production, which leads to the formation of osteoclasts with bone-resorbing function. Furthermore, IL-17 induces chemokines with an impact on osteogenic differentiation [75]. We detected accumulation of six proteins from IL-17 pathway in differentiated cells (Table 1); particularly, chemokines (CXCL1, CXCL5, CXCL6 1.16, 2.04, 2.67-fold increase, respectively), matrix metallopeptidases (MMP1 and MMP3, 3.09 and 4.04-fold increase, respectively), and heat shock protein 90 (1.99-fold increase). Unique chemokines specifically affect bone remodeling, such as CXCL1 and CXCL5 function in osteoclastogenesis [76]. Additionally, MMPs are recognized osteoclastic resorption factors. The authors reported that MMPs inhibition can even block osteoclastogenesis [77,78]. Accumulation of proteins belonging to IL-17 pathway suggests the developmental switch on 24 days after triggering differentiation. Interestingly, it was reported that 1,25-dihydroxyvitamin D3 and thyroid hormones modulated MMP13 [79]—a link to a biosynthetic pathway described in the following section.
Biosynthesis of thyroid hormones and steroids. Among other factors, bone turnover is regulated by numerous hormones which maintain a fragile balance between bone formation and resorption. Hormones can act directly or indirectly as growth factors, receptor-binding molecules, and synthesis modulators. They include thyroid and parathyroid hormones, estrogens, growth hormones, and 1,25-dihydroxyvitamin D3. Thyroid-stimulating hormone (TSH) works via TSH receptors, which are activated by protein disulfide isomerase [80]. Additionally, TSH affects secondary messengers (for example, Ca2+ signaling cascade through IP3) in extrathyroidal tissues. It was reported that in embryonic stem cells, TSH stimulated osteoblastic differentiation via protein kinase C and noncanonical Wnt pathway. TSH is directly connected to bone homeostasis and development, for example, thyrotoxicosis in adults was accompanied by secondary osteoporosis and increased risk of fractures [81]. We detected five proteins involved in thyroid hormone synthesis pathway all upregulated in differentiated cells (Table 1). Namely, albumin, heat shock protein 70, inositol 1,4,5-trisphosphate receptor type 2, protein disulfide isomerase, and transthyretin (2.45, 1.76, 4.92, 1.50, and 2.51-fold increase, respectively). Interestingly, human transthyretin, TSH transporter, is involved in the bone mineralization process and its destabilization leads to changed crystal morphology [82]. All steroid hormones are synthesized from cholesterol. We identified three members of the steroid biosynthesis pathway accumulated in differentiated cells (Table 1): 7-dehydrocholesterol reductase (DHCR7), methylsterol monooxygenase 1, and NAD(P)-dependent steroid dehydrogenase-like (2.77, 2.38, 2.47-fold increase, respectively), which are involved in cholesterol metabolism. In general, LDL cholesterol promotes osteoclastogenesis and vice versa, HDL cholesterol protects osteoblasts from apoptosis [83]. Interestingly, cholesterol-mediated degradation of DHCR7 switched the cholesterol synthesis to vitamin D production [84] causing improved bone formation [85]. DHCR7 is the terminal enzyme of cholesterol synthesis converting 7-dehydrocholesterol (7DHC) to cholesterol. Since our data showed a high abundance of DHCR7 and apolipoprotein A2 (ApoA2) from upregulated PPAR pathway (discussed below in Section 2.4.2) in osteoblasts, we hypothesize that the former routes ApoA2 to HDL cholesterol for stabilizing structure.
Mineral absorption pathway. Bone remodeling is an active energy-consuming process, which requires a supply of macro- and micronutrients such as vitamins and minerals. It is well recognized that calcium is essential for preventing osteoporosis. The positive role of micronutrients in bone metabolism is dual: minerals can impact hydroxyapatite formation or can function as regulators and cofactors [86]. Three proteins belonging to the mineral absorption pathway were more abundant in osteodifferentiated cells (Table 1) (ferritin heavy chain, STEAP2 metalloreductase and metallothionein 1E, 1.77, 2.60, and 6.56-fold increase, respectively). Metallothionein 1E (MT1E), a metal-binding protein, has different roles. For example, it maintains the concentrations of free intracellular zinc/copper ions and acts as an antioxidant. Specifically, one of the probable functions is controlling bone differentiation through regulating the availability of zinc inside the cells [87]. Lin et al. reported that ZnCl2 promoted DPSCs differentiation through upregulation of MT1E [88]. We revealed the massive abundance of MT1E in osteoblasts confirming its essential role in osteogenesis. Iron overload can inhibit osteodifferentiation due to ferritin overexpression [89]. In contrast, yet in line with our results, Hu et al. reported more abundant ferritin heavy chain in osteodifferentiated cell [90]. Our data suggest a high potential of nutritional therapeutics for bone disorders.
Summing up, FoxO signaling may prevent apoptosis by inducing antioxidants, thus counterbalancing proinflammatory chemokines from the IL-17 pathway. Hence, differentially abundant chemokines and matrix metallopeptidases from IL-17 pathway are potential pharmacology targets for preventing osteoclastogenesis. Thyroid-stimulating hormone might promote differentiation of osteoblasts through activation of its receptors and accumulated transporter. Finally, the accumulation of proteins belonging to the mineral absorption pathway suggests a high potential for nutritional therapy of bone disorders.

2.4.2. Lipid Metabolism and PPAR Signaling Are Plausibly Connected to Energy Metabolism in the Context of Bone Development

PPAR signaling pathway. The nuclear peroxisome proliferator-activated receptors (PPARs) regulate many physiological processes, including adipogenesis, lipid metabolism, insulin sensitivity, inflammation, angiogenesis, and osteoclastogenesis. Additionally, they are key modulators of skeletal remodeling. PPARs belong to the superfamily of ligand-inducible nuclear hormone receptors, which includes three major members [91] with different expression profiles: (i) PPARα participates in fatty acid metabolism and lowers level of lipids, (ii) PPARβ/δ mainly influence fatty acid oxidation and transport, while (iii) PPARγ is primarily involved in the regulation of the adipogenesis and energy balance [91]. We identified differentially accumulated proteins from cascades of all three PPAR members (Table 1). PPARγ plays a negative role in osteoblast differentiation by operating in an adipogenesis-dependent manner due to the relationship between osteoblast and adipocyte differentiation. When PPARγ expression was ablated in mature osteoblasts and osteocytes, the mutant mice exhibited increased femoral bone mineral density and trabecular bone volume as well as reduced fat mass and higher energy expenditure [92]. In contrast, pharmacological induction of PPARα expression concomitantly enhanced expression of bone morphogenetic protein 2 (BMP-2) and calcium deposition [93]. We identified accumulation of ApoA2 and ApoC3 (5.14 and 2.18-fold increase, respectively) working downstream of PPARα pathway in osteoblasts, the former additionally validated by immunocytochemistry (Figure 3). Of note, it was reported that ApoA1 deficient mice developed osteoarthritis [94]. Information about ApoA2 and ApoC3 impact on osteogenesis is rare. ApoA2 constitutes roughly 20% of HDL cholesterol and stabilizes its structure, whereas ApoA1 is a major cholesterol-binding protein. Histological studies indicated that both mature osteoblasts and differentiating osteoblast progenitors contain stored lipids. Moreover, in vivo and in vitro studies showed that osteoblasts uptake circulating lipoproteins and free fatty acids [95].
Fatty acid biosynthesis and degradation. Besides peroxisomes, fatty acid biosynthesis and degradation take place in mitochondria. These pathways include acetyl-CoA acyltransferase 2, enoyl-CoA delta isomerase 2, enoyl-CoA hydratase short chain 1, acyl-CoA synthetase long chain family member 3, acyl-CoA synthetase long chain family member 4, hydroxyacyl-CoA dehydrogenase, acyl-CoA oxidase 1, acyl-CoA synthetase 2, and acyl-CoA oxidase 3 (1.06, 3.70, 1.45, 2.17, 2.03, 1.16, 2.41, 2.13, 3.11-fold increase, respectively). Overall, fatty acids represent one of the major energy sources involved in bone turnover. Furthermore, van Gastel et al. [96] proposed a role for fatty acid utilization during fracture healing and the determination of skeletal cell fate. In addition, the oxidation of fatty acids has an impact on bone metabolism and calcification [97]. Energy role of fatty acids on the mineralization phase of osteodifferentiation is further supported by the accumulation of amino acids valine, isoleucine, and threonine (discussed in Section 2.5), which could be precursors for lipogenesis through TCA cycle. Of note, our results are in line with the proteomic study of Wang et al. [98]. These authors reported a more prominent PPAR signaling and fatty acid biosynthesis pathway in in vitro osteoinduced DPSCs comparing to periodontal ligament stem cells.
Thermogenesis. Although thermogenesis pathway was not enriched, its identification underlined the context of molecular mechanism in osteoblast. The term thermogenesis primarily refers to the production of heat for maintaining body temperature (thermoregulation). It is also widely used to describe the regulation of the energy balance in the whole body [99]. We detected eleven members of thermogenesis pathway, accumulated in differentiated cells (Table 1). Particularly, proteins participating in energy production, ACSL3 and ACLS4 (2.17, 2.03-fold increase, respectively) as members of β-oxidation, ATP5MG (2.11-fold increase) responsible for ATP synthesis in mitochondria, NDUFS1, NDUFS3, NDUFV2, and NDUFAF2 (1.02, 1.01, 1.1, 2.97-fold increase, respectively) involved in electron transport chain in mitochondria, SDHA and SDHB (1.13, 1.47-fold increase, respectively) enzymes of TCA cycle. Bone remodeling is an energy-consuming process due to the synthesis of extracellular matrix proteins and the accumulation of mineral ions for hydroxyapatite crystal formation [95]. Interestingly, the authors reported that insulin signaling affected both bone metabolism and energy production [59]. Energy is harvested from different sources (glucose, glutamine, fatty acids, as well as autophagy) but the role of a particular substrate on the specific phase of the osteodifferentiation is not completely clear yet. Osteoclasts’ differentiation and their fusion are likely accompanied by oxidative phosphorylation [100], while another study reported increased oxidative phosphorylation during the differentiation of osteoblasts from bone marrow stem cells [101]. Our data pointed to fatty acids as major sources for energy production in the mineralization phase of osteodifferentiation.
Our data suggested the essential role of lipids on the mineralization phase of osteodifferentiation: (i) apolipoproteins A2 and C3 from PPARα cascade together with cholesterol biosynthesis enzyme DHCR7 stabilize HDL cholesterol and (ii) accumulated proteins involved in biosynthesis and degradation of fatty acids ensure preferential catabolism of lipids via thermogenesis pathway.

2.5. Pathways Potentially Essential for Osteodifferentiation Hinted by Metabolomic Data

Using metabolomics, we identified eight statistically significant differentially regulated metabolites, of which seven were amino acids and one acetate. All metabolites and/or metabolic pathways were upregulated in osteodifferentiated cells (Figure 7 and Figure S3). Only aminoacyl-tRNA biosynthesis and valine, leucine, isoleucine metabolism reported more than one metabolite incorporated in the pathway (Table S7, Figure S3).
Aminoacyl-tRNA synthesis. Genetic information is transformed into functional proteins through transcription and translation. Aminoacyl-tRNA synthesis plays a key role in determining how the genetic code is interpreted [102]. This metabolic pathway determined on the 24th day of osteodifferentiation is upregulated in differentiated cells. Affected protein synthesis was also reported by Klontzas et al. in a metabolomic study of osteodifferentiation [103]. We detected aminoacyl-tRNA biosynthesis at both levels, by incorporating proteins and metabolites into functional pathways. We identified two members on protein level (histidyl-tRNA synthetase and glutaminyl-tRNA synthetase, 4.20 and 2.0-fold increase in differentiated cells, respectively) and seven members by metabolomic analysis (histidine, phenylalanine, valine, isoleucine, threonine, tryptophan, and proline; 2.3, 2.96, 3.5, 4.0, 1.4, 4.8, and 2.9-fold increase, respectively). Both substrate and enzyme were determined only in the case of histidine and histidyl-tRNA synthetase, which may indicate their prominent role in the osteogenic process.
Histidine and phenylalanine metabolism. Histidyl-tRNA synthetase is necessary for the synthesis of histidyl-transfer RNA, which is essential for histidine incorporation into proteins [104]. HIStidin-rich proteins also include histatins, proteins with various function, among other with wound-healing properties. For example, histatins were used for the treatment of oral diseases [105]. It was reported that human salivary Hst1 protein promoted BMP-2 induced bone regeneration [106]. In our study, BMP-2 protein was detected (complementation) both by Western blotting and immunocytochemistry (Figure 3 and Figure 6). We also determined histidine as a solitary member of histidine and beta-alanine metabolism (2.3-fold increase in osteodifferentiated cells). Histidine metabolism contains several pathways. The most important is catabolism via urocanate to glutamate, relevant for HIS-rich proteins and HIS-containing dipeptides synthesis [105], especially L-carnosine (L-Car). L-Car is synthetized from L-histamine and beta-alanine and its zinc chelated form β-alanyl-L-histidinato zinc is a pharmacological tool for treatment of osteoporosis [107]. Chauhan and Singh reported that the presence of histidine led to hydroxyapatite mineralization with plate-like morphology in studies with mouse osteoblast precursor cells [108]. Hydroxyapatite is a basic building block of bones and teeth. These findings correspond to our observations. On the 24th day after application of osteogenic medium differentiated cells reached mineralization phase accompanied with production of calcium deposits (Figure 4). Aromatic amino acid phenylalanine (2.96-fold increase in differentiated cells) is a member of three metabolic pathways. The authors showed that supplementation of C57BL/6 mice model by aromatic amino acids (phenylalanine, tyrosine, and tryptophan) prevented bone loss [109].
Valine, leucine, and isoleucine biosynthesis and degradation. As we described in part thermogenesis, bone remodeling is energy consuming process and the involvement of certain amino acids in metabolism can bring another source of energy. Isoleucine, valine, and threonine (4.0, 3.5, and 1.4-fold increase, respectively) are three from nine essential amino acids. Their products can enter TCA cycle for energy production or they can be precursors for lipogenesis and/or ketone bodies production through acetyl-CoA and acetoacetate [110].
In summary, analysis of metabolites indicated their plausible support for protein synthesis and energy metabolism.

2.6. Proposed Proteins and Metabolites for Functional Validation in Patient Samples

Although Wang et al. [98] argued that different stem cell sources might employ different major pathways for the osteodifferentiation, we addressed this limitation by integrating and discussing our data in the context of published literature on different biological material and experimental approaches. Hence, we propose that the particular proteins or metabolites in pathways with high enrichment score are potential targets for validation in patient samples. Moreover, we suggest that several proteins/metabolites involved in osteogenesis and/or bone loss belonging to discussed pathways could be promising biomarkers. These are: (i) IP3R, histidine, decorin, and inhibin βA from known pathways; and (ii) apolipoproteins (ApoA2/C3), 7-dehydrocholesterol reductase, matrix metallopeptidases (MMP1/3), chemokines (CXCL1/5/6), metallothionein 1E, signal transducers and activators of transcription (STAT1/3), and PRMT1 from enriched pathways.

3. Materials and Methods

3.1. Cell Cultures and Osteogenic Treatment

The human commercial cell line, DPSCs (Lonza, Basel, Switzerland) was used. For all studies examining the molecular mechanism of osteogenic differentiation, cells were seeded in 25 cm2 flasks and grown under a humidified atmosphere of 5% CO2 at 37 °C. Cells were precultured in DPSCs growth medium consisted of α-MEM (Merck Sigma-Aldrich, Darmstadt, Germany) supplemented with 2 mM GlutaMAX (Thermo Fisher Scientific, Waltham, MA, USA), 0.2 mM ascorbic acid 2-phosphate (Sigma-Aldrich), 10% fetal bovine serum (Biosera, Nuaille, France), 100 U/mL penicillin, and 100 µg/mL streptomycin (Biosera). The medium was replaced every three days. During preculture, cells remained in an undifferentiated stage. After reaching 80–90% confluence, cells were harvested using TrypLE Express Enzyme (Thermo Fisher Scientific) and then seeded at 5000 cells/cm2. When DPSCs reached 50% confluence (day 0), the growth medium was replaced with osteoinductive medium OsteoMAX-XF (Merck Sigma-Aldrich) for osteogenic differentiation and mineralization. The cells were grown in this specific medium for 24 days. In parallel, undifferentiated DPSCs (control) were cultured in the DPSCs growth medium. Media were changed every 3 days. Cell morphology was monitored during the experiment by light microscopy (Optika, Via Rigla, Italy).

3.2. Alizarin Red S Staining and Quantification of Calcium Mineral Deposits

Control and differentiated cells were grown in 24 well plates. On a specific day (5th, 11th, 15th, 20th, 24th) after osteogenic medium application, cells were fixed by 4% paraformaldehyde (Cell Signaling Technology, Danvers, MA, USA) for 30 min at RT. After DPBS (Gibco, Paisley, UK) washing, formed Ca2+ deposits were stained by 2% Alizarin Red S (Merck Sigma-Aldrich, Darmstadt, Germany) for 3 min at RT. Next, the dye was aspired and cells were carefully washed three times by ultrapure water to remove the remaining stain. Finally, cells were kept in ultrapure water. Color intensity was monitored and images were captured by the inverted microscope XDS-2 (Optika). For quantification of total stained mineralized tissue, Alizarin Red S was eluted by destain solution (10% acetic acid, 20% methanol) for 15 min at RT [111]. The eluates were centrifuged (10,000× g) for 10 min at RT to remove mineral precipitates and cell debris. Then, the supernatant was transferred into 96 well plates in duplicates. The absorbance was measured at 405 nm using Synergy H1 Microplate Reader (BioTek, Winooski, VT, USA). The amount of mineral deposits is presented as optical density normalized to the total amount of proteins (mg), which were determined by BCA assay (Thermo Fisher Scientific).

3.3. Immunocytochemistry

Immunocytochemistry was performed according to the modified protocol [21]. Briefly, cells were seeded in lumox 96-multiwell plates (Sarstedt, Nümbrecht, Germany). On the 24th day after osteogenic induction, cells were washed with DPBS (Gibco) and fixed with 4% paraformaldehyde (Cell Signaling Technology) for 30 min at room temperature (RT). Then, fixed cells were washed with DPBS for 10 min, permeabilized with 0.2% Triton-X in DPBS (washing buffer, WB) and blocked with 5% goat serum in WB for 1 h at RT. Consequently, samples were incubated with appropriately diluted primary antibodies in WB at 4 °C overnight. Then, the cells were washed with WB three times for 10 min and further incubated with goat anti-rabbit Alexa Fluor 488-conjugated IgG (1:500, #ab150077, Abcam, Cambridge, UK) or goat anti-mouse Alexa Fluor 488-conjugated IgG (1:500, #ab150113, Abcam) secondary antibody in WB for 1 h at RT in the dark. After washing, nuclei were stained with DAPI (1 µg/mL) for 10 min at RT in the dark. Cells were again washed as described above, cellular proteins and nuclei were visualized by fluorescent microscope WiScan (Hermes IDEA Bio-Medical, Rehovot, Israel), magnification 100×.
The following primary antibodies were used: anti-bone morphogenetic protein 2, BMP-2 (1:100, #ab14933, Abcam), anti-apolipoprotein A2, ApoA2 (1:50, #ab92478, Abcam), anti-bone sialoprotein, BSP II (1:50, #5468, Cell Signaling Technology), and anti-osteocalcin, OSTCN (1:80, #ab13418, Abcam).

3.4. Protein Extraction

Control and differentiated cells were harvested in biological triplicates on the 24th day after application of the osteogenic medium. Cells were washed with ice-cold DPBS (Gibco) and scraped-off from flasks using SDS buffer (4% SDS, 0.1 M Tris pH 7.6, 100 mM DTT) supplemented with protease (Roche, Basel, Switzerland) and phosphatase inhibitors (Sigma-Aldrich). Then, lysates were transferred to the tubes, incubated at 10 °C for 15 min and sonicated (3 × 10 s, 20% pulse) followed by 20 min centrifugation (14,000× g) at 12 °C. The final supernatant was aliquoted and stored at −80 °C until use. Protein concentration was determined by the Bradford detergent compatible assay (Thermo Fisher Scientific). Control and differentiated cells were processed in parallel.

3.5. SDS-PAGE and Western Blotting

Equal amounts (40 µg) of protein aliquots in the Laemmli buffer (BioRad, Hercules, CA, USA) from both control (undifferentiated cells) and osteoblasts (differentiated cells) were loaded on Tris-glycine gradient (4–20%) or 12% polyacrylamide gels (BioRad) and electrophoretically separated in Mini-PROTEAN Tetra cell chamber (BioRad). Separation was performed at 60 V for 20 min, followed by 200 V until tracking dye reached the bottom of the gel. Then, proteins profiled in the gel were transferred to a nitrocellulose membrane (Thermo Fisher Scientific) using a semi-dry blotting unit (Biometra, Jena, Germany) and transferring condition 1 mA/cm2. The membrane was blocked (5% non-fat dry milk or 5% BSA in Tris-buffered saline with Tween 20, TBS-T) for 1 h at room temperature and incubated overnight with primary antibodies to target or reference proteins at 4 °C. Antibodies were diluted in 5% non-fat dry milk or 5% BSA. After incubation, the membrane was thoroughly washed with TBS-T three times and further incubated with horseradish peroxidase-conjugated goat anti-rabbit IgG (1:2000, #7074, Cell Signaling Technology) or horseradish peroxidase-conjugated goat anti-mouse IgG (1:2500, #ab6789, Abcam) for 2 h at RT. Proteins were detected by chemiluminescence after application of the developing solution according to manufacturer’s instructions (#34080, Thermo Fisher Scientific) and visualized using ChemiDoc XRS+ system (BioRad). Densitometric analysis of protein abundance was accomplished via ImageLab software (version 6.0) (BioRad). Signals were normalized to glyceraldehyde 3-phosphate dehydrogenase (GAPDH).
The following primary antibodies were used: anti-dentin sialophosphoprotein, DSPP (1:500, #PA5-76382, Invitrogen, Waltham, MA, USA), anti-dentin matrix protein 1, DMP-1 (1:500, #PA5-103323, Invitrogen), anti-bone morphogenetic protein 2, BMP-2 (1:250, #ab14933, Abcam), anti-collagen type VI alpha-1 chain, Col6A1 (1:250, #sc-377143, Santa Cruz Biotechnology, Dallas, TX, USA), anti-inositol 1,4,5-trisphosphate receptor, IP3R (1:250, #8568, Cell Signaling Technology), and anti-GAPDH (1:500, #sc-166574, Santa Cruz Biotechnology).

3.6. Protein Digestion by Filter-Aided Sample Preparation

Total protein extracts (50 µg) adjusted to 200 µL with urea buffer (8 M urea in 50 mM ammonium bicarbonate) were digested using filter-aided sample preparation [112] on Microcon Ultracel YM-10 centrifugal filters (Merck Sigma Aldrich, Darmstadt, Germany). Before digestion, the filters were pre-washed by 50 mM ammonium bicarbonate. Samples were incubated overnight at 37 °C with trypsin gold (Promega, Madison, WI, USA) in the proportion 1:75 to total proteins and consequently slightly acidified by 10% trifluoroacetic acid. The peptide mixtures were purified on Sep-Pak C18 light columns (Waters, Millford, MA, USA) and vacuum-concentrated. The concentration of peptides was measured spectrophotometrically by NanoDrop 2000 (Thermo Fisher Scientific) and adjusted to 100 ng/μL with 0.1% trifluoroacetic acid.

3.7. Relative Label-Free Quantification by Mass Spectrometry

Complex peptide mixtures (400 ng) were separated using Acquity M-Class UHPLC (Waters). Firstly, samples were desalted and concentrated on the nanoEase Symmetry C18 trap column: 20 mm length, 180 μm diameter, 5 μm particle size (Waters). Secondly, peptides were comprehensively profiled on the nanoEase HSS T3 C18 analytical column: 250 mm length, 75 μm diameter, 1.8 μm particle size (Waters) using a 118 min gradient of 5–35% acetonitrile with 0.1% formic acid at a flow rate of 300 nL/min. Next, samples were nanosprayed by 2.9 kV applied to PicoTip emitter (New Objective, Woburn, MA, USA) to the quadrupole time-of-flight mass spectrometer Synapt G2-Si (Waters). The instrument operated in resolution mode with enabled ion mobility, which provided additional gas-phase separation of ions, particularly useful for filtering non-peptide contaminants. Synapt G2-Si was tuned using parameters described in the literature [112] for maximum sensitivity and both mass and mobility resolutions. Spectra were recorded in a data-independent mode referred as high definition MSE. Ions with 50–2000 m/z were detected in both low energy and high energy channels. Voltage profile setting on the first quadrupole allowed efficient deflection of low mass precursor ions ˂ 400 m/z. The external standard Glu1-Fibrinopeptide B was infused for mass correction.
Compression and Archival Tool 1.0 (Waters) reduced noise, removing ion counts below 15. Subsequent processing was done in Progenesis QI 4.1 (Waters). Peaks were modeled with low energy threshold 400 counts and high energy threshold 40 counts. Chromatographic elution profiles in low and high energy traces were correlated to assign fragment ions to appropriate precursor ions. Retention times of peaks were aligned to automatically selected reference chromatogram. Next, peak intensities were normalized to the median distribution of all ions. The label-free quantification relied on the integration of peak areas of the three most intense precursor peptides. Arcsinh transformation of data was done before statistics. We used Ion Accounting 4.0 (Waters) search algorithm for protein identification. Spectra were searched against human proteome sequences downloaded from UniProt in September 2019 (74,823 entries, uniprot.org). Search parameters were: (i) one trypsin miscleavage, (ii) fixed carbamidomethyl cysteine, (iii) variable oxidized methionine and deamidated asparagine/glutamine, (iv) automatic mass tolerance, (v) 4% false discovery rate against the randomized database. Protein hits having unique peptides were accepted if at least two reliable peptides (score ≥ 5.93, mass accuracy ≤ 20 ppm) matched the sequence. Integrated principal component analysis visualized clustering of biological replicates according to experimental groups.

3.8. Metabolite Extraction

Control and differentiated cells were harvested in biological triplicates on the 24th day after application of osteogenic medium according to the modified protocol [113]. Briefly, cells were washed with cold DPBS (Gibco) and scraped-off from the flask using 90% methanol. Then, lysed cells were transferred to the tubes, incubated 15 min at −20 °C followed by 20 min centrifugation (16,000× g) at 4 °C. The final supernatant was stored at −80 °C until use. Control and differentiated cells were processed in parallel.

3.9. Identification and Quantification of Metabolites by Nuclear Magnetic Resonance

The supernatant was vacuum-dried and then the residue was dissolved in 550 µL of 200 mM phosphate buffer pH 7.4, containing 0.2 mM TMS-d4 (trimethylsilylpropionic acid-d4) as chemical shift reference, in deuterated water.
Data were acquired on 600 MHz NMR spectrometer Avance III (Bruker, Bremen, Germany) equipped with cryoprobe at 310 K. Samples were freshly prepared before acquisition and measured in random order. Before measurement, each sample was prewarmed to 310 K for 5 min. An exponential noise filter was used to introduce 0.3 Hz line broadening before Fourier transformation. The proton NMR chemical shifts are reported relative to TMS-d4 signal and data were zero filled once. For each sample, one-dimensional as well as two-dimensional NMR spectra were acquired. We modified standard profiling protocol (Bruker) as follows: (i) noesy with presaturation (FID size 64 k, dummy scans 4, number of scans 4, spectral width 20.4750 ppm); (ii) cosy with presaturation (FID size 4 k, dummy scans 8, number of scans 8, spectral width 16.0125 ppm); (iii) homonuclear j-resolved (FID size 8 k, dummy scans 16, number of scans 64, spectral width 16.6087); (iv) profiling cpmg (FID size 64 k, dummy scans 4, number of scans 2048, spectral width 20.0156 ppm). All experiments were conducted with a relaxation delay of 4 s.
All spectra were binned from 0.00 ppm to 10.00 ppm with bins 0.001 ppm. The multiplicity of peaks was confirmed in j-resolved spectra, and homonuclear cross peaks were confirmed in cosy spectra. Spectra were solved using Human Metabolome Database (www.hmda.ca, accessed on 15 September 2019) [114], internal metabolite database, Chenomx software (version 8.31), and metabolomic literature [115]. After identification of metabolites, we chose spectra subregions with unique assignments or minimally affected by interfering signals. In 0.001 ppm binned spectra, we summed integrals of selected metabolites. These data represent a relative concentration of a particular metabolite in the sample. Metabolites showing peaks with low intensity or strong peaks´ overlap were excluded from the evaluation. Next, metabolites´ relative abundances were normalized to the total amount of proteins [116].

3.10. Statistical Analysis and Bioinformatics

Differentially abundant proteins were filtered on effect size (at least 1-fold of log2 transformed ratio) and ANOVA p ≤ 0.01 (false discovery rate Q ˂ 1.05%). For statistical evaluation of metabolomic and densitometry data, we used Student’s t-test at adjusted (Benjamini-Hochberg correction) p ≤ 0.05 with the same effect size threshold. Differentially abundant proteins were sorted into biological processes using GOTermFinder (version 0.86) [117] and metabolic pathways by KEGG (release 97.0) [118]. Proteins with identification numbers absent in KEGG were queried for homology using BLASTP. Differentially abundant metabolites were uploaded in MetaboAnalyst (version 5.0) for functional analysis [119]. Pathway enrichment analysis was performed in R programming language (version 4.0.3) by ClusterProfiler package (version 3.18.1) [120] of R/Bioconductor (version 3.11) using a permutation test (p ≤ 0.05) without stringent false discovery rate filtering as discussed in specialized literature [121]. Dot plot was generated by enrichplot (version 1.10.2) R/Bioconductor package. A volcano plot was generated using VolcaNoseR (version 1.0.3) [122].

4. Conclusions

Understanding of molecular mechanisms which participate in various diseases and normal development is necessary both for expanding fundamental knowledge and for developing new treatment strategies. Complex proteomic and metabolomic data allowed us to identify important proteins and metabolites and to suggest their role in the osteodifferentiation, summarized in Figure 8. In essence, we explained the molecular mechanisms behind late osteoblasts’ differentiation (mineralization phase) by discussing several known (described in the literature) as well as unconventional signaling pathways in relation to osteogenesis. Our proteomic data indicated that the osteodifferentiation in the mineralization phase is energy consuming, whereas the cells harvested energy probably from fatty acids. Furthermore, this process likely required intensive antioxidant protection and was plausibly regulated by hormones. Metabolomic data enriched results by highlighting pathways related to energy metabolism and protein synthesis. Of note, the presented network of signaling pathways is largely hypothetical and should be validated in follow-up functional experiments. For instance, inhibitors of particular signaling pathways may be used for direct verification if specific molecular mechanisms indeed promote or inhibit osteodifferentiation of mesenchymal stem cells.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms22157908/s1.

Author Contributions

Conceptualization, S.N.; methodology, S.N., M.D., T.O., E.B., and H.Š.; formal analysis, M.D., M.G., and G.B.; investigation, S.N., M.D., T.O., E.B., J.P., and M.J.; validation, S.N., J.S., and H.Š.; resources, A.K., E.H., and H.Š.; data curation, A.K. and E.B.; writing—original draft, S.N.; writing—review and editing, M.D. and H.Š.; visualization, S.N.; supervision, S.N. and H.Š.; project administration, H.Š. and E.H.; funding acquisition, E.H. and H.Š. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by projects VEGA 1/0310/21 of the Scientific Grant Agency of the Ministry of Education, Science, Research and Sports of the Slovak Republic and the Slovak Academy of Sciences and APVV-15-0217 from the Slovak Research and Development Agency. This publication also has been produced with the support of the Integrated Infrastructure Operational Program for the project: Systemic Public Research Infrastructure—Biobank for Cancer and Rare diseases, ITMS: 313011AFG5, co-financed by the European Regional Development Fund.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analyzed during this study are included in the main text of the article and its electronic Supplementary Materials. The mass spectrometry proteomics data were deposited to the ProteomeXchange Consortium via PRIDE partner repository [123] with dataset identifier PXD025223. NMR spectra are available on request from Eva Baranovičová.

Acknowledgments

The authors are thankful to Alžbeta Kráľová Trančíková and Kristina Nystøl for excellent technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Manolagas, S.C. Birth and Death of Bone Cells: Basic Regulatory Mechanisms and Implications for the Pathogenesis and Treatment of Osteoporosis. Endocr. Rev. 2000, 21, 115–137. [Google Scholar] [CrossRef] [Green Version]
  2. Office of the Surgeon General (US). The basics of bone in health and disease. In Bone Health and Osteoporosis: A Report of the Surgeon General; Department of Health and Human Services, Office of the Surgeon General: Rockville, MD, USA, 2004; pp. 16–38. [Google Scholar]
  3. Tu, K.N.; Lie, J.D.; Wan, C.K.V.; Cameron, M.; Austel, A.G.; Nguyen, J.K.; Van, K.; Hyun, D. Osteoporosis: A Review of Treatment Options. Pharm. Ther. 2018, 43, 92–104. [Google Scholar]
  4. Abdel Meguid, E.; Ke, Y.; Ji, J.; El-Hashash, A.H.K. Stem Cells Applications in Bone and Tooth Repair and Regeneration: New Insights, Tools, and Hopes. J. Cell. Physiol. 2018, 233, 1825–1835. [Google Scholar] [CrossRef] [PubMed]
  5. Sobacchi, C.; Schulz, A.; Coxon, F.P.; Villa, A.; Helfrich, M.H. Osteopetrosis: Genetics, Treatment and New Insights into Osteoclast Function. Nat. Rev. Endocrinol. 2013, 9, 522–536. [Google Scholar] [CrossRef]
  6. Raggatt, L.J.; Partridge, N.C. Cellular and Molecular Mechanisms of Bone Remodeling. J. Biol. Chem. 2010, 285, 25103–25108. [Google Scholar] [CrossRef] [Green Version]
  7. Iñiguez-Ariza, N.M.; Clarke, B.L. Bone Biology, Signaling Pathways, and Therapeutic Targets for Osteoporosis. Maturitas 2015, 82, 245–255. [Google Scholar] [CrossRef] [PubMed]
  8. Hoang-Kim, A.; Gelsomini, L.; Luciani, D.; Moroni, A.; Giannini, S. Fracture Healing and Drug Therapies in Osteoporosis. Clin. Cases Miner. Bone Metab. 2009, 6, 136–143. [Google Scholar]
  9. Aspenberg, P. Drugs and Fracture Repair. Acta Orthop. 2005, 76, 741–748. [Google Scholar] [CrossRef] [Green Version]
  10. Iaquinta, M.R.; Mazzoni, E.; Bononi, I.; Rotondo, J.C.; Mazziotta, C.; Montesi, M.; Sprio, S.; Tampieri, A.; Tognon, M.; Martini, F. Adult Stem Cells for Bone Regeneration and Repair. Front. Cell Dev. Biol. 2019, 7. [Google Scholar] [CrossRef] [Green Version]
  11. Antebi, B.; Pelled, G.; Gazit, D. Stem Cell Therapy for Osteoporosis. Curr. Osteoporos. Rep. 2014, 12, 41–47. [Google Scholar] [CrossRef] [PubMed]
  12. Pavone, V.; Testa, G.; Giardina, S.M.C.; Vescio, A.; Restivo, D.A.; Sessa, G. Pharmacological Therapy of Osteoporosis: A Systematic Current Review of Literature. Front. Pharmacol. 2017, 8, 803. [Google Scholar] [CrossRef]
  13. Hegde, V.; Jo, J.E.; Andreopoulou, P.; Lane, J.M. Effect of Osteoporosis Medications on Fracture Healing. Osteoporos. Int. 2016, 27, 861–871. [Google Scholar] [CrossRef] [PubMed]
  14. Testa, G.; Pavone, V.; Mangano, S.; Riccioli, M.; Arancio, A.; Evola, F.R.; Avonda, S.; Sessa, G. Normal Nutritional Components and Effects on Bone Metabolism in Prevention of Osteoporosis. J. Biol. Regul. Homeost. Agents 2015, 29, 729–736. [Google Scholar] [PubMed]
  15. Caplan, A.I.; Dennis, J.E. Mesenchymal Stem Cells as Trophic Mediators. J. Cell. Biochem. 2006, 98, 1076–1084. [Google Scholar] [CrossRef] [PubMed]
  16. Oryan, A.; Kamali, A.; Moshiri, A.; Baghaban Eslaminejad, M. Role of Mesenchymal Stem Cells in Bone Regenerative Medicine: What Is the Evidence? Cells Tissues Organs 2017, 204, 59–83. [Google Scholar] [CrossRef]
  17. Hayrapetyan, A.; Jansen, J.A.; van den Beucken, J.J.J.P. Signaling Pathways Involved in Osteogenesis and Their Application for Bone Regenerative Medicine. Tissue Eng. Part B Rev. 2015, 21, 75–87. [Google Scholar] [CrossRef]
  18. Yang, B.; Qiu, Y.; Zhou, N.; Ouyang, H.; Ding, J.; Cheng, B.; Sun, J. Application of Stem Cells in Oral Disease Therapy: Progresses and Perspectives. Front. Physiol. 2017, 8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Bindal, P.; Kasim, N.H.A.; Ramasamy, T.S.; Dabbagh, A.; Moharamzadeh, K.; Chai, W.L. 18—Dental pulp tissue engineering and regenerative endodontic therapy. In Biomaterials for Oral and Dental Tissue Engineering; Tayebi, L., Moharamzadeh, K., Eds.; Woodhead Publishing: Sawston, UK, 2017; pp. 297–318. ISBN 978-0-08-100961-1. [Google Scholar]
  20. Liu, J.; Yu, F.; Sun, Y.; Jiang, B.; Zhang, W.; Yang, J.; Xu, G.-T.; Liang, A.; Liu, S. Concise Reviews: Characteristics and Potential Applications of Human Dental Tissue-Derived Mesenchymal Stem Cells. Stem Cells 2015, 33, 627–638. [Google Scholar] [CrossRef]
  21. Okajcekova, T.; Strnadel, J.; Pokusa, M.; Zahumenska, R.; Janickova, M.; Halasova, E.; Skovierova, H. A Comparative In Vitro Analysis of the Osteogenic Potential of Human Dental Pulp Stem Cells Using Various Differentiation Conditions. Int. J. Mol. Sci. 2020, 21, 2280. [Google Scholar] [CrossRef] [Green Version]
  22. Liu, M.; Sun, Y.; Liu, Y.; Yuan, M.; Zhang, Z.; Hu, W. Modulation of the Differentiation of Dental Pulp Stem Cells by Different Concentrations of β-Glycerophosphate. Molecules 2012, 17, 1219–1232. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Karbanová, J.; Soukup, T.; Suchánek, J.; Pytlík, R.; Corbeil, D.; Mokrý, J. Characterization of Dental Pulp Stem Cells from Impacted Third Molars Cultured in Low Serum-Containing Medium. Cells Tissues Organs 2011, 193, 344–365. [Google Scholar] [CrossRef] [Green Version]
  24. Kim, B.-C.; Bae, H.; Kwon, I.-K.; Lee, E.-J.; Park, J.-H.; Khademhosseini, A.; Hwang, Y.-S. Osteoblastic/Cementoblastic and Neural Differentiation of Dental Stem Cells and Their Applications to Tissue Engineering and Regenerative Medicine. Tissue Eng. Part B Rev. 2012, 18, 235–244. [Google Scholar] [CrossRef] [Green Version]
  25. Karbanová, J.; Soukup, T.; Suchánek, J.; Mokrý, J. Osteogenic Differentiation of Human Dental Pulp-Derived Stem Cells under Various Ex-Vivo Culture Conditions. Acta Med. 2010, 53, 79–84. [Google Scholar] [CrossRef]
  26. Kim, J.; Lee, G.; Chang, W.S.; Ki, S.H.; Park, J.-C. Comparison and Contrast of Bone and Dentin in Genetic Disorder, Morphology and Regeneration: A Review. J. Bone Metab. 2021, 28, 1–10. [Google Scholar] [CrossRef]
  27. MacDougall, M.J.; Javed, A. Dentin and Bone: Similar Collagenous Mineralized Tissues. In Bone and Development; Bronner, F., Farach-Carson, M.C., Roach, H.I., Eds.; Topics in Bone Biology; Springer: London, UK, 2010; pp. 183–200. ISBN 978-1-84882-822-3. [Google Scholar]
  28. Li, Z.; Yu, M.; Tian, W. An Inductive Signalling Network Regulates Mammalian Tooth Morphogenesis with Implications for Tooth Regeneration. Cell Prolif. 2013, 46, 501–508. [Google Scholar] [CrossRef] [PubMed]
  29. Sharma, V.; Srinivasan, A.; Nikolajeff, F.; Kumar, S. Biomineralization Process in Hard Tissues: The Interaction Complexity within Protein and Inorganic Counterparts. Acta Biomater. 2021, 120, 20–37. [Google Scholar] [CrossRef] [PubMed]
  30. Maraldi, T.; Riccio, M.; Pisciotta, A.; Zavatti, M.; Carnevale, G.; Beretti, F.; La Sala, G.B.; Motta, A.; De Pol, A. Human Amniotic Fluid-Derived and Dental Pulp-Derived Stem Cells Seeded into Collagen Scaffold Repair Critical-Size Bone Defects Promoting Vascularization. Stem Cell Res. Ther. 2013, 4, 53. [Google Scholar] [CrossRef] [Green Version]
  31. D’Aquino, R.; De Rosa, A.; Lanza, V.; Tirino, V.; Laino, L.; Graziano, A.; Desiderio, V.; Laino, G.; Papaccio, G. Human Mandible Bone Defect Repair by the Grafting of Dental Pulp Stem/Progenitor Cells and Collagen Sponge Biocomplexes. Eur. Cell Mater. 2009, 18, 75–83. [Google Scholar] [CrossRef]
  32. Campos, J.M.; Sousa, A.C.; Caseiro, A.R.; Pedrosa, S.S.; Pinto, P.O.; Branquinho, M.V.; Amorim, I.; Santos, J.D.; Pereira, T.; Mendonça, C.M.; et al. Dental Pulp Stem Cells and Bonelike® for Bone Regeneration in Ovine Model. Regen. Biomater. 2019, 6, 49–59. [Google Scholar] [CrossRef] [PubMed]
  33. Koga, T.; Minamizato, T.; Kawai, Y.; Miura, K.-I.; Takashi, I.; Nakatani, Y.; Sumita, Y.; Asahina, I. Bone Regeneration Using Dentin Matrix Depends on the Degree of Demineralization and Particle Size. PLoS ONE 2016, 11, e0147235. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Matsuoka, F.; Takeuchi, I.; Agata, H.; Kagami, H.; Shiono, H.; Kiyota, Y.; Honda, H.; Kato, R. Morphology-Based Prediction of Osteogenic Differentiation Potential of Human Mesenchymal Stem Cells. PLoS ONE 2013, 8, e55082. [Google Scholar] [CrossRef] [PubMed]
  35. Carluccio, M.; Ziberi, S.; Zuccarini, M.; Giuliani, P.; Caciagli, F.; Di Iorio, P.; Ciccarelli, R. Adult Mesenchymal Stem Cells: Is There a Role for Purine Receptors in Their Osteogenic Differentiation? Purinergic Signal. 2020, 16, 263–287. [Google Scholar] [CrossRef] [PubMed]
  36. Kirkham, G.R.; Cartmell, S. Genes and Proteins Involved in the Regulation of Osteogenesis. Available online: https://www.oulu.fi/spareparts/ebook_topics_in_t_e_vol3/index.html (accessed on 24 June 2021).
  37. Chau, J.F.L.; Leong, W.F.; Li, B. Signaling Pathways Governing Osteoblast Proliferation, Differentiation and Function. Histol. Histopathol. 2009, 24, 1593–1606. [Google Scholar] [CrossRef]
  38. Sanpaolo, E.R.; Rotondo, C.; Cici, D.; Corrado, A.; Cantatore, F.P. JAK/STAT Pathway and Molecular Mechanism in Bone Remodeling. Mol. Biol. Rep. 2020, 47, 9087–9096. [Google Scholar] [CrossRef]
  39. Tajima, K.; Takaishi, H.; Takito, J.; Tohmonda, T.; Yoda, M.; Ota, N.; Kosaki, N.; Matsumoto, M.; Ikegami, H.; Nakamura, T.; et al. Inhibition of STAT1 Accelerates Bone Fracture Healing. J. Orthop. Res. 2010, 28, 937–941. [Google Scholar] [CrossRef] [PubMed]
  40. Xu, L.; Zhang, L.; Zhang, H.; Yang, Z.; Qi, L.; Wang, Y.; Ren, S. The Participation of Fibroblast Growth Factor 23 (FGF23) in the Progression of Osteoporosis via JAK/STAT Pathway. J. Cell. Biochem. 2018, 119, 3819–3828. [Google Scholar] [CrossRef] [PubMed]
  41. Yang, Y.; Chung, M.R.; Zhou, S.; Gong, X.; Xu, H.; Hong, Y.; Jin, A.; Huang, X.; Zou, W.; Dai, Q.; et al. STAT3 Controls Osteoclast Differentiation and Bone Homeostasis by Regulating NFATc1 Transcription. J. Biol. Chem. 2019, 294, 15395–15407. [Google Scholar] [CrossRef] [PubMed]
  42. De Paula, F.J.A.; Rosen, C.J. Bone Remodeling and Energy Metabolism: New Perspectives. Bone Res. 2013, 1, 72–84. [Google Scholar] [CrossRef] [Green Version]
  43. Alves, R.D.A.M.; Eijken, M.; Swagemakers, S.; Chiba, H.; Titulaer, M.K.; Burgers, P.C.; Luider, T.M.; van Leeuwen, J.P.T.M. Proteomic Analysis of Human Osteoblastic Cells: Relevant Proteins and Functional Categories for Differentiation. J. Proteome Res. 2010, 9, 4688–4700. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Cao, C.; Ren, Y.; Barnett, A.S.; Mirando, A.J.; Rouse, D.; Mun, S.H.; Park-Min, K.-H.; McNulty, A.L.; Guilak, F.; Karner, C.M.; et al. Increased Ca2+ Signaling through CaV1.2 Promotes Bone Formation and Prevents Estrogen Deficiency–Induced Bone Loss. JCI Insight 2017, 2. [Google Scholar] [CrossRef] [Green Version]
  45. Krishnan, V.; Bryant, H.U.; Macdougald, O.A. Regulation of Bone Mass by Wnt Signaling. J. Clin. Investig. 2006, 116, 1202–1209. [Google Scholar] [CrossRef]
  46. Kim, J.H.; Liu, X.; Wang, J.; Chen, X.; Zhang, H.; Kim, S.H.; Cui, J.; Li, R.; Zhang, W.; Kong, Y.; et al. Wnt Signaling in Bone Formation and Its Therapeutic Potential for Bone Diseases. Ther. Adv. Musculoskelet. Dis. 2013, 5, 13–31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Okamoto, M.; Udagawa, N.; Uehara, S.; Maeda, K.; Yamashita, T.; Nakamichi, Y.; Kato, H.; Saito, N.; Minami, Y.; Takahashi, N.; et al. Noncanonical Wnt5a Enhances Wnt/β-Catenin Signaling during Osteoblastogenesis. Sci. Rep. 2014, 4, 4493. [Google Scholar] [CrossRef] [Green Version]
  48. Chen, G.; Deng, C.; Li, Y.-P. TGF-β and BMP Signaling in Osteoblast Differentiation and Bone Formation. Int. J. Biol. Sci. 2012, 8, 272–288. [Google Scholar] [CrossRef] [Green Version]
  49. Gaddy-Kurten, D.; Coker, J.K.; Abe, E.; Jilka, R.L.; Manolagas, S.C. Inhibin Suppresses and Activin Stimulates Osteoblastogenesis and Osteoclastogenesis in Murine Bone Marrow Cultures. Endocrinology 2002, 143, 74–83. [Google Scholar] [CrossRef]
  50. Waddington, R.J.; Roberts, H.C.; Sugars, R.V.; Schönherr, E. Differential Roles for Small Leucine-Rich Proteoglycans in Bone Formation. Eur. Cell. Mater. 2003, 6, 12–21, discussion 21. [Google Scholar] [CrossRef]
  51. Wu, M.; Chen, G.; Li, Y.-P. TGF-β and BMP Signaling in Osteoblast, Skeletal Development, and Bone Formation, Homeostasis and Disease. Bone Res. 2016, 4, 1–21. [Google Scholar] [CrossRef] [PubMed]
  52. Dyne, K.M.; Valli, M.; Forlino, A.; Mottes, M.; Kresse, H.; Cetta, G. Deficient Expression of the Small Proteoglycan Decorin in a Case of Severe/Lethal Osteogenesis Imperfecta. Am. J. Med. Genet. 1996, 63, 161–166. [Google Scholar] [CrossRef]
  53. Yamakoshi, Y.; Kinoshita, S.; Izuhara, L.; Karakida, T.; Fukae, M.; Oida, S. DPP and DSP Are Necessary for Maintaining TGF-Β1 Activity in Dentin. J. Dent. Res. 2014, 93, 671–677. [Google Scholar] [CrossRef] [Green Version]
  54. Yamakoshi, Y.; Simmer, J.P. Structural Features, Processing Mechanism and Gene Splice Variants of Dentin Sialophosphoprotein. Jpn. Dent. Sci. Rev. 2018, 54, 183–196. [Google Scholar] [CrossRef] [PubMed]
  55. Teti, G.; Salvatore, V.; Ruggeri, A.; Manzoli, L.; Gesi, M.; Orsini, G.; Falconi, M. In Vitro Reparative Dentin: A Biochemical and Morphological Study. Eur. J. Histochem. 2013, 57. [Google Scholar] [CrossRef]
  56. Verdelis, K.; Ling, Y.; Sreenath, T.; Haruyama, N.; MacDougall, M.; van der Meulen, M.C.H.; Lukashova, L.; Spevak, L.; Kulkarni, A.B.; Boskey, A.L. DSPP Effects on in Vivo Bone Mineralization. Bone 2008, 43, 983–990. [Google Scholar] [CrossRef] [Green Version]
  57. Jani, P.H.; Gibson, M.P.; Liu, C.; Zhang, H.; Wang, X.; Lu, Y.; Qin, C. Transgenic Expression of Dspp Partially Rescued the Long Bone Defects of Dmp1-Null Mice. Matrix Biol. 2016, 52–54, 95–112. [Google Scholar] [CrossRef] [Green Version]
  58. Yang, J.; Zhang, X.; Wang, W.; Liu, J. Insulin Stimulates Osteoblast Proliferation and Differentiation through ERK and PI3K in MG-63 Cells. Cell Biochem. Funct. 2010, 28, 334–341. [Google Scholar] [CrossRef] [PubMed]
  59. Pramojanee, S.N.; Phimphilai, M.; Chattipakorn, N.; Chattipakorn, S.C. Possible Roles of Insulin Signaling in Osteoblasts. Endocr. Res. 2014, 39, 144–151. [Google Scholar] [CrossRef] [PubMed]
  60. Ge, C.; Xiao, G.; Jiang, D.; Franceschi, R.T. Critical Role of the Extracellular Signal–Regulated Kinase–MAPK Pathway in Osteoblast Differentiation and Skeletal Development. J. Cell. Biol. 2007, 176, 709–718. [Google Scholar] [CrossRef] [PubMed]
  61. Wu, H.; Teng, P.-N.; Jayaraman, T.; Onishi, S.; Li, J.; Bannon, L.; Huang, H.; Close, J.; Sfeir, C. Dentin Matrix Protein 1 (DMP1) Signals via Cell Surface Integrin. J. Biol. Chem. 2011, 286, 29462–29469. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Zhang, S.; Wan, H.; Wang, P.; Liu, M.; Li, G.; Zhang, C.; Sun, Y. Extracellular Matrix Protein DMP1 Suppresses Osteogenic Differentiation of Mesenchymal Stem Cells. Biochem. Biophys. Res. Commun. 2018, 501, 968–973. [Google Scholar] [CrossRef]
  63. McGonnell, I.M.; Grigoriadis, A.E.; Lam, E.W.-F.; Price, J.S.; Sunters, A. A Specific Role for Phosphoinositide 3-Kinase and AKT in Osteoblasts? Front. Endocrinol. 2012, 3. [Google Scholar] [CrossRef] [Green Version]
  64. Tsai, K.-S.; Kao, S.-Y.; Wang, C.-Y.; Wang, Y.-J.; Wang, J.-P.; Hung, S.-C. Type I Collagen Promotes Proliferation and Osteogenesis of Human Mesenchymal Stem Cells via Activation of ERK and Akt Pathways. J. Biomed. Mater. Res. A 2010, 94, 673–682. [Google Scholar] [CrossRef] [PubMed]
  65. Izu, Y.; Ezura, Y.; Koch, M.; Birk, D.E.; Noda, M. Collagens VI and XII Form Complexes Mediating Osteoblast Interactions during Osteogenesis. Cell Tissue Res. 2016, 364, 623–635. [Google Scholar] [CrossRef] [Green Version]
  66. Chen, Y.-Y.; He, S.-T.; Yan, F.-H.; Zhou, P.-F.; Luo, K.; Zhang, Y.-D.; Xiao, Y.; Lin, M.-K. Dental Pulp Stem Cells Express Tendon Markers under Mechanical Loading and Are a Potential Cell Source for Tissue Engineering of Tendon-like Tissue. Int. J. Oral Sci. 2016, 8, 213–222. [Google Scholar] [CrossRef] [Green Version]
  67. Chen, D.; Gong, Y.; Xu, L.; Zhou, M.; Li, J.; Song, J. Bidirectional Regulation of Osteogenic Differentiation by the FOXO Subfamily of Forkhead Transcription Factors in Mammalian MSCs. Cell Prolif. 2019, 52, e12540. [Google Scholar] [CrossRef] [Green Version]
  68. Ambrogini, E.; Almeida, M.; Martin-Millan, M.; Paik, J.-H.; Depinho, R.A.; Han, L.; Goellner, J.; Weinstein, R.S.; Jilka, R.L.; O’Brien, C.A.; et al. FoxO-Mediated Defense against Oxidative Stress in Osteoblasts Is Indispensable for Skeletal Homeostasis in Mice. Cell Metab. 2010, 11, 136–146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Redza-Dutordoir, M.; Averill-Bates, D.A. Activation of Apoptosis Signalling Pathways by Reactive Oxygen Species. Biochim. Biophys. Acta 2016, 1863, 2977–2992. [Google Scholar] [CrossRef]
  70. Ma, X.; Su, P.; Yin, C.; Lin, X.; Wang, X.; Gao, Y.; Patil, S.; War, A.R.; Qadir, A.; Tian, Y.; et al. The Roles of FoxO Transcription Factors in Regulation of Bone Cells Function. Int. J. Mol. Sci. 2020, 21, 692. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Yamagata, K.; Daitoku, H.; Takahashi, Y.; Namiki, K.; Hisatake, K.; Kako, K.; Mukai, H.; Kasuya, Y.; Fukamizu, A. Arginine Methylation of FOXO Transcription Factors Inhibits Their Phosphorylation by Akt. Mol. Cell. 2008, 32, 221–231. [Google Scholar] [CrossRef]
  72. Wang, Z.; Yu, T.; Huang, P. Post-Translational Modifications of FOXO Family Proteins (Review). Mol. Med. Rep. 2016, 14, 4931–4941. [Google Scholar] [CrossRef] [PubMed]
  73. Sergi, C.; Shen, F.; Liu, S.-M. Insulin/IGF-1R, SIRT1, and FOXOs Pathways—An Intriguing Interaction Platform for Bone and Osteosarcoma. Front. Endocrinol. 2019, 10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Feng, Y.-L.; Jiang, X.-T.; Ma, F.-F.; Han, J.; Tang, X.-L. Resveratrol Prevents Osteoporosis by Upregulating FoxO1 Transcriptional Activity. Int. J. Mol. Med. 2018, 41, 202–212. [Google Scholar] [CrossRef] [Green Version]
  75. Lee, Y. The Role of Interleukin-17 in Bone Metabolism and Inflammatory Skeletal Diseases. BMB Rep. 2013, 46, 479–483. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Brylka, L.J.; Schinke, T. Chemokines in Physiological and Pathological Bone Remodeling. Front. Immunol. 2019, 10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Almalki, S.G.; Agrawal, D.K. Effects of Matrix Metalloproteinases on the Fate of Mesenchymal Stem Cells. Stem Cell Res. Ther. 2016, 7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Hayami, T.; Kapila, Y.L.; Kapila, S. MMP-1 (Collagenase-1) and MMP-13 (Collagenase-3) Differentially Regulate Markers of Osteoblastic Differentiation in Osteogenic Cells. Matrix Biol. 2008, 27, 682–692. [Google Scholar] [CrossRef] [Green Version]
  79. Fratzl-Zelman, N.; Glantschnig, H.; Rumpler, M.; Nader, A.; Ellinger, A.; Varga, F. The Expression of Matrix Metalloproteinase-13 and Osteocalcin in Mouse Osteoblasts Is Related to Osteoblastic Differentiation and Is Modulated by 1,25-Dihydroxyvitamin D3 and Thyroid Hormones. Cell Biol. Int. 2003, 27, 459–468. [Google Scholar] [CrossRef]
  80. Campos, J.L.O.; Doratioto, T.R.; Videira, N.B.; Ribeiro Filho, H.V.; Batista, F.A.H.; Fattori, J.; Indolfo, N.d.C.; Nakahira, M.; Bajgelman, M.C.; Cvoro, A.; et al. Protein Disulfide Isomerase Modulates the Activation of Thyroid Hormone Receptors. Front. Endocrinol. 2019, 9. [Google Scholar] [CrossRef] [Green Version]
  81. Bassett, J.H.D.; Williams, G.R. Role of Thyroid Hormones in Skeletal Development and Bone Maintenance. Endocr. Rev. 2016, 37, 135–187. [Google Scholar] [CrossRef] [Green Version]
  82. Wieczorek, E.; Chitruń, A.; Ożyhar, A. Destabilised Human Transthyretin Shapes the Morphology of Calcium Carbonate Crystals. Biochim. Biophys. Acta Gen. Subj. 2019, 1863, 313–324. [Google Scholar] [CrossRef]
  83. Yin, W.; Li, Z.; Zhang, W. Modulation of Bone and Marrow Niche by Cholesterol. Nutrients 2019, 11, 1394. [Google Scholar] [CrossRef] [Green Version]
  84. Prabhu, A.V.; Luu, W.; Sharpe, L.J.; Brown, A.J. Cholesterol-Mediated Degradation of 7-Dehydrocholesterol Reductase Switches the Balance from Cholesterol to Vitamin D Synthesis. J. Biol. Chem. 2016, 291, 8363–8373. [Google Scholar] [CrossRef] [Green Version]
  85. Laird, E.; Ward, M.; McSorley, E.; Strain, J.J.; Wallace, J. Vitamin D and Bone Health; Potential Mechanisms. Nutrients 2010, 2, 693–724. [Google Scholar] [CrossRef] [Green Version]
  86. Pepa, G.D.; Brandi, M.L. Microelements for Bone Boost: The Last but Not the Least. Clin. Cases Miner. Bone Metab. 2016, 13, 181–185. [Google Scholar] [CrossRef] [PubMed]
  87. Sankaramanivel, S.; Molgakar, M.Y.; Prakash, P.A.; Jaabir, M.S.M.; Gurunathan, S. Stem Cells and Metallothionein-A Review. Int. J. Curr. Res. Rev. 2017, 9, 54–61. [Google Scholar] [CrossRef]
  88. Lin, C.-Y.; Lin, H.-H.; Tsai, M.-H.; Lin, S.-P.; Chen, M.-H. Zinc Chloride for Odontogenesis of Dental Pulp Stem Cells via Metallothionein Up-Regulation. J. Endodont. 2011, 37, 211–216. [Google Scholar] [CrossRef]
  89. Zarjou, A.; Jeney, V.; Arosio, P.; Poli, M.; Zavaczki, E.; Balla, G.; Balla, J. Ferritin Ferroxidase Activity: A Potent Inhibitor of Osteogenesis. J. Bone Miner. Res. 2010, 25, 164–172. [Google Scholar] [CrossRef]
  90. Hu, J.-J.; Liu, Y.-W.; He, M.-Y.; Jin, D.; Zhao, H.; Yu, B. Proteomic Analysis on Effectors Involved in BMP-2-Induced Osteogenic Differentiation of Beagle Bone Marrow Mesenchymal Stem Cells. Proteome Sci. 2014, 12, 13. [Google Scholar] [CrossRef] [Green Version]
  91. Grygiel-Górniak, B. Peroxisome Proliferator-Activated Receptors and Their Ligands: Nutritional and Clinical Implications—A Review. Nutr. J. 2014, 13, 17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Brun, J.; Berthou, F.; Trajkovski, M.; Maechler, P.; Foti, M.; Bonnet, N. Bone Regulates Browning and Energy Metabolism through Mature Osteoblast/Osteocyte PPARγ Expression. Diabetes 2017, 66, 2541–2554. [Google Scholar] [CrossRef] [Green Version]
  93. Kim, Y.-H.; Jang, W.-G.; Oh, S.-H.; Kim, J.-W.; Lee, M.N.; Song, J.H.; Yang, J.-W.; Zang, Y.; Koh, J.-T. Fenofibrate Induces PPARα and BMP2 Expression to Stimulate Osteoblast Differentiation. Biochem. Biophys. Res. Commun. 2019, 520, 459–465. [Google Scholar] [CrossRef]
  94. Blair, H.C.; Kalyvioti, E.; Papachristou, N.I.; Tourkova, I.L.; Syggelos, S.A.; Deligianni, D.; Orkoula, M.G.; Kontoyannis, C.G.; Karavia, E.A.; Kypreos, K.E.; et al. Apolipoprotein A-1 Regulates Osteoblast and Lipoblast Precursor Cells in Mice. Lab. Investig. 2016, 96, 763–772. [Google Scholar] [CrossRef]
  95. Alekos, N.S.; Moorer, M.C.; Riddle, R.C. Dual Effects of Lipid Metabolism on Osteoblast Function. Front. Endocrinol. 2020, 11. [Google Scholar] [CrossRef]
  96. Van Gastel, N.; Stegen, S.; Eelen, G.; Schoors, S.; Carlier, A.; Daniëls, V.W.; Baryawno, N.; Przybylski, D.; Depypere, M.; Stiers, P.-J.; et al. Lipid Availability Determines Fate of Skeletal Progenitor Cells via SOX9. Nature 2020, 579, 111–117. [Google Scholar] [CrossRef] [PubMed]
  97. Kim, S.P.; Li, Z.; Zoch, M.L.; Frey, J.L.; Bowman, C.E.; Kushwaha, P.; Ryan, K.A.; Goh, B.C.; Scafidi, S.; Pickett, J.E.; et al. Fatty Acid Oxidation by the Osteoblast Is Required for Normal Bone Acquisition in a Sex- and Diet-Dependent Manner. JCI Insight 2017, 2. [Google Scholar] [CrossRef]
  98. Wang, H.; Ma, D.; Zhang, X.; Xu, S.; Ning, T.; Wu, B. Comparative Proteomic Profiling of Human Dental Pulp Stem Cells and Periodontal Ligament Stem Cells under in Vitro Osteogenic Induction. Arch. Oral Biol. 2018, 89, 9–19. [Google Scholar] [CrossRef] [PubMed]
  99. Trayhurn, P. Thermogenesis. In Encyclopedia of Food Sciences and Nutrition, 2nd ed.; Caballero, B., Ed.; Academic Press: Oxford, UK, 2003; pp. 5762–5767. ISBN 978-0-12-227055-0. [Google Scholar]
  100. Lee, W.-C.; Guntur, A.R.; Long, F.; Rosen, C.J. Energy Metabolism of the Osteoblast: Implications for Osteoporosis. Endocr. Rev. 2017, 38, 255–266. [Google Scholar] [CrossRef] [PubMed]
  101. Shum, L.C.; White, N.S.; Mills, B.N.; Bentley, K.L.d.M.; Eliseev, R.A. Energy Metabolism in Mesenchymal Stem Cells During Osteogenic Differentiation. Stem Cells Dev. 2016, 25, 114–122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  102. Ibba, M.; Soll, D. Aminoacyl-TRNA Synthesis. Annu. Rev. Biochem. 2000, 69, 617–650. [Google Scholar] [CrossRef]
  103. Klontzas, M.E.; Vernardis, S.I.; Heliotis, M.; Tsiridis, E.; Mantalaris, A. Metabolomics Analysis of the Osteogenic Differentiation of Umbilical Cord Blood Mesenchymal Stem Cells Reveals Differential Sensitivity to Osteogenic Agents. Stem Cells Dev. 2017, 26, 723–733. [Google Scholar] [CrossRef]
  104. Freist, W.; Verhey, J.F.; Rühlmann, A.; Gauss, D.H.; Arnez, J.G. Histidyl-TRNA Synthetase. Biol. Chem. 1999, 380, 623–646. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  105. Holeček, M. Histidine in Health and Disease: Metabolism, Physiological Importance, and Use as a Supplement. Nutrients 2020, 12, 848. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  106. Sun, P.; Shi, A.; Shen, C.; Liu, Y.; Wu, G.; Feng, J. Human Salivary Histatin-1 (Hst1) Promotes Bone Morphogenetic Protein 2 (BMP2)-Induced Osteogenesis and Angiogenesis. FEBS Open Bio 2020, 10, 1503–1515. [Google Scholar] [CrossRef] [PubMed]
  107. Yamaguchi, M. β-Alanyl-l-Histidinato Zinc and Bone Resorption. Gen. Pharmacol. 1995, 26, 1179–1183. [Google Scholar] [CrossRef]
  108. Chauhan, N.; Singh, Y. L-Histidine Controls the Hydroxyapatite Mineralization with Plate-like Morphology: Effect of Concentration and Media. Mater. Sci. Eng. 2021, 120, 111669. [Google Scholar] [CrossRef]
  109. Refaey, M.E.; Zhong, Q.; Ding, K.-H.; Shi, X.; Xu, J.; Bollag, W.B.; Hill, W.D.; Chutkan, N.; Robbins, R.; Nadeau, H.; et al. Impact of Dietary Aromatic Amino Acids on Osteoclastic Activity. Calcif. Tissue Int. 2014, 95, 174–182. [Google Scholar] [CrossRef] [Green Version]
  110. Manoli, I.; Venditti, C.P. Disorders of Branched Chain Amino Acid Metabolism. Transl. Sci. Rare Dis. 2016, 1, 91–110. [Google Scholar] [CrossRef] [Green Version]
  111. Kunimatsu, R.; Nakajima, K.; Awada, T.; Tsuka, Y.; Abe, T.; Ando, K.; Hiraki, T.; Kimura, A.; Tanimoto, K. Comparative Characterization of Stem Cells from Human Exfoliated Deciduous Teeth, Dental Pulp, and Bone Marrow-Derived Mesenchymal Stem Cells. Biochem. Biophys. Res. Commun. 2018, 501, 193–198. [Google Scholar] [CrossRef] [PubMed]
  112. Distler, U.; Kuharev, J.; Navarro, P.; Tenzer, S. Label-Free Quantification in Ion Mobility-Enhanced Data-Independent Acquisition Proteomics. Nat. Protoc. 2016, 11, 795–812. [Google Scholar] [CrossRef]
  113. Sapcariu, S.C.; Kanashova, T.; Weindl, D.; Ghelfi, J.; Dittmar, G.; Hiller, K. Simultaneous Extraction of Proteins and Metabolites from Cells in Culture. MethodsX 2014, 1, 74–80. [Google Scholar] [CrossRef] [PubMed]
  114. Wishart, D.S.; Feunang, Y.D.; Marcu, A.; Guo, A.C.; Liang, K.; Vázquez-Fresno, R.; Sajed, T.; Johnson, D.; Li, C.; Karu, N.; et al. HMDB 4.0: The Human Metabolome Database for 2018. Nucleic Acids Res. 2018, 46, D608–D617. [Google Scholar] [CrossRef] [PubMed]
  115. Nagana Gowda, G.A.; Gowda, Y.N.; Raftery, D. Expanding the Limits of Human Blood Metabolite Quantitation Using NMR Spectroscopy. Anal. Chem. 2015, 87, 706–715. [Google Scholar] [CrossRef] [Green Version]
  116. Fritsche-Guenther, R.; Bauer, A.; Gloaguen, Y.; Lorenz, M.; Kirwan, J.A. Modified Protocol of Harvesting, Extraction, and Normalization Approaches for Gas Chromatography Mass Spectrometry-Based Metabolomics Analysis of Adherent Cells Grown Under High Fetal Calf Serum Conditions. Metabolites 2019, 10, 2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  117. Boyle, E.I.; Weng, S.; Gollub, J.; Jin, H.; Botstein, D.; Cherry, J.M.; Sherlock, G. GO::TermFinder—Open Source Software for Accessing Gene Ontology Information and Finding Significantly Enriched Gene Ontology Terms Associated with a List of Genes. Bioinformatics 2004, 20, 3710–3715. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  118. Kanehisa, M.; Goto, S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  119. Chong, J.; Soufan, O.; Li, C.; Caraus, I.; Li, S.; Bourque, G.; Wishart, D.S.; Xia, J. MetaboAnalyst 4.0: Towards More Transparent and Integrative Metabolomics Analysis. Nucleic Acids Res. 2018, 46, W486–W494. [Google Scholar] [CrossRef] [Green Version]
  120. Yu, G.; Wang, L.-G.; Han, Y.; He, Q.-Y. ClusterProfiler: An R Package for Comparing Biological Themes among Gene Clusters. OMICS J. Integr. Biol. 2012, 16, 284–287. [Google Scholar] [CrossRef] [PubMed]
  121. Reimand, J.; Isserlin, R.; Voisin, V.; Kucera, M.; Tannus-Lopes, C.; Rostamianfar, A.; Wadi, L.; Meyer, M.; Wong, J.; Xu, C.; et al. Pathway Enrichment Analysis and Visualization of Omics Data Using g:Profiler, GSEA, Cytoscape and EnrichmentMap. Nat. Protoc. 2019, 14, 482–517. [Google Scholar] [CrossRef]
  122. Goedhart, J.; Luijsterburg, M.S. VolcaNoseR Is a Web App for Creating, Exploring, Labeling and Sharing Volcano Plots. Sci. Rep. 2020, 10, 20560. [Google Scholar] [CrossRef]
  123. Perez-Riverol, Y.; Csordas, A.; Bai, J.; Bernal-Llinares, M.; Hewapathirana, S.; Kundu, D.J.; Inuganti, A.; Griss, J.; Mayer, G.; Eisenacher, M.; et al. The PRIDE Database and Related Tools and Resources in 2019: Improving Support for Quantification Data. Nucleic Acids Res. 2019, 47, D442–D450. [Google Scholar] [CrossRef] [PubMed]
Ijms 22 07908 g001 550
Figure 1. Experimental design—analytical workflow. We followed undifferentiated and differentiated human dental pulp stem cells (DPSCs). Osteogenic differentiation was verified by calcium deposits stain, Western blotting, and immunocytochemistry. Proteins extracted from three biological replicates were digested using filter-aided sample preparation (FASP). Following that, we directly quantified proteins by mass spectrometry upon comprehensive two-dimensional separation of analytes. Metabolites were directly quantified by nuclear magnetic resonance (NMR). Dedicated KEGG, GOTermFinder, and MetaboAnalyst tools accompanied with enrichment analysis were used for revealing the biological meaning of the extensive dataset.
Figure 1. Experimental design—analytical workflow. We followed undifferentiated and differentiated human dental pulp stem cells (DPSCs). Osteogenic differentiation was verified by calcium deposits stain, Western blotting, and immunocytochemistry. Proteins extracted from three biological replicates were digested using filter-aided sample preparation (FASP). Following that, we directly quantified proteins by mass spectrometry upon comprehensive two-dimensional separation of analytes. Metabolites were directly quantified by nuclear magnetic resonance (NMR). Dedicated KEGG, GOTermFinder, and MetaboAnalyst tools accompanied with enrichment analysis were used for revealing the biological meaning of the extensive dataset.
Ijms 22 07908 g001
Ijms 22 07908 g002 550
Figure 2. The visual appearance of undifferentiated (control) and differentiated cells on the 24th day after osteogenic medium application. Controls were cultured in a DPSCs growth medium during the experiment. Cellular morphology correlates with the osteogenic differentiation potential. Representative micrographs were taken by phase contrast technique at magnification 100×.
Figure 2. The visual appearance of undifferentiated (control) and differentiated cells on the 24th day after osteogenic medium application. Controls were cultured in a DPSCs growth medium during the experiment. Cellular morphology correlates with the osteogenic differentiation potential. Representative micrographs were taken by phase contrast technique at magnification 100×.
Ijms 22 07908 g002
Ijms 22 07908 g003 550
Figure 3. Content of osteogenic markers and apolipoprotein A2. Cells were fixed with 4% paraformaldehyde. Specific proteins apolipoprotein A2 (ApoA2), osteocalcin (OSTCN), bone sialoprotein (BSP II), and bone morphogenetic protein (BMP-2) were detected by immunocytochemistry using WiScan microscope. Nuclei were visualized by DAPI stain. Bar = 50 µm. Areas with white rectangles are shown as magnified insets below each micrograph.
Figure 3. Content of osteogenic markers and apolipoprotein A2. Cells were fixed with 4% paraformaldehyde. Specific proteins apolipoprotein A2 (ApoA2), osteocalcin (OSTCN), bone sialoprotein (BSP II), and bone morphogenetic protein (BMP-2) were detected by immunocytochemistry using WiScan microscope. Nuclei were visualized by DAPI stain. Bar = 50 µm. Areas with white rectangles are shown as magnified insets below each micrograph.
Ijms 22 07908 g003
Ijms 22 07908 g004 550
Figure 4. Morphological changes of DPSCs during osteogenic differentiation on the 5th, 11th, 15th, 20th, and 24th day and quantitative evaluation of Ca2+ deposits. The control cells (upper line) and osteodifferentiated cells (lower line) were stained by 2% Alizarin Red S and monitored by light microscopy at magnification 100× (A). The mineral deposits production was quantified at each time point spectrophotometrically and represented as a relative ratio to control cells (B).
Figure 4. Morphological changes of DPSCs during osteogenic differentiation on the 5th, 11th, 15th, 20th, and 24th day and quantitative evaluation of Ca2+ deposits. The control cells (upper line) and osteodifferentiated cells (lower line) were stained by 2% Alizarin Red S and monitored by light microscopy at magnification 100× (A). The mineral deposits production was quantified at each time point spectrophotometrically and represented as a relative ratio to control cells (B).
Ijms 22 07908 g004
Ijms 22 07908 g005 550
Figure 5. KEGG pathway analysis. Dot plot (A) of enriched KEGG pathways (results from ClusterProfiler) produced by enrichplot Bioconductor package. The color of the dots represents the p-value, the size of the dots represents the number of differentially abundant proteins. Column chart (B) of KEGG pathway with a known role in osteogenesis. Specific proteins belonging to these functional groups are detailed in Table 1.
Figure 5. KEGG pathway analysis. Dot plot (A) of enriched KEGG pathways (results from ClusterProfiler) produced by enrichplot Bioconductor package. The color of the dots represents the p-value, the size of the dots represents the number of differentially abundant proteins. Column chart (B) of KEGG pathway with a known role in osteogenesis. Specific proteins belonging to these functional groups are detailed in Table 1.
Ijms 22 07908 g005
Ijms 22 07908 g006 550
Figure 6. Confirmation of several differentially abundant proteins by Western blotting. Densitometric analysis of protein abundance was performed using ImageLab software. Values represent the means of three biological replicates. Error bars denote the standard deviations, * indicate p ≤ 0.05. Arrows indicate densitometrically analyzed major protein bands: DSP a likely cleavage product of DSPP, BMP-2 mature protein, and functional IP3R tetramer. Abbreviations: DSP (dentin sialoprotein), DMP-1 (dentin matrix protein 1), BMP-2 (bone morphogenetic protein 2), IP3R (inositol 1,4,5-trisphosphate receptor), Col6A1 (collagen 6 alpha 1 chain), GAPDH (glyceraldehyde 3-phosphate dehydrogenase).
Figure 6. Confirmation of several differentially abundant proteins by Western blotting. Densitometric analysis of protein abundance was performed using ImageLab software. Values represent the means of three biological replicates. Error bars denote the standard deviations, * indicate p ≤ 0.05. Arrows indicate densitometrically analyzed major protein bands: DSP a likely cleavage product of DSPP, BMP-2 mature protein, and functional IP3R tetramer. Abbreviations: DSP (dentin sialoprotein), DMP-1 (dentin matrix protein 1), BMP-2 (bone morphogenetic protein 2), IP3R (inositol 1,4,5-trisphosphate receptor), Col6A1 (collagen 6 alpha 1 chain), GAPDH (glyceraldehyde 3-phosphate dehydrogenase).
Ijms 22 07908 g006
Ijms 22 07908 g007 550
Figure 7. Volcano plot of metabolomic data. Differentially abundant metabolites are highlighted.
Figure 7. Volcano plot of metabolomic data. Differentially abundant metabolites are highlighted.
Ijms 22 07908 g007
Ijms 22 07908 g008 550
Figure 8. Summary of the proposed molecular mechanisms of osteodifferentiation. Potential biomarkers are bold and italic denoted. Abbreviations: ALB (albumin), ACSBG2 (acyl-CoA synthetase bubblegum family member 2), ACSL1, ACSL3 (acyl-CoA synthetase long chain family members 1 and 3), ANGPTL4 (angiopoietin like 4), ACOX1, ACOX3 (acyl-CoA oxidases 1 and 3), APOA2, APOC3 (apolipoproteins A2 and C3), CXCL1, CXCL5, CXCL6 (C-X-C motif chemokine ligands 1, 5, and 6), FHL1 (four and a half LIM domains 1), FTH1 (ferritin heavy chain 1), GFAP (glial fibrillary acidic protein), GK (glycerol kinase), HSP90 (heat shock protein 90), HSPA5 (heat shock protein family A (Hsp70) member 5), INS (insulin), IP3R (inositol 1,4,5-trisphosphate receptor type 2), MMP1, MMP3 (matrix metallopeptidases 1 and 3), MSMO1 (methylsterol monooxygenase 1), MT1E (metallothionein 1E), NSDHL (NAD(P) dependent steroid dehydrogenase-like), PDIA4 (protein disulfide isomerase family A member 4), PLIN2 (perilipin 2), PRMT1 (protein arginine methyltransferase 1), SOD2 (superoxide dismutase 2), STAT1, STAT3 (signal transducers and activators of transcription 1 and 3), STEAP2 (STEAP2 metalloreductase), and TTR (transthyretin).
Figure 8. Summary of the proposed molecular mechanisms of osteodifferentiation. Potential biomarkers are bold and italic denoted. Abbreviations: ALB (albumin), ACSBG2 (acyl-CoA synthetase bubblegum family member 2), ACSL1, ACSL3 (acyl-CoA synthetase long chain family members 1 and 3), ANGPTL4 (angiopoietin like 4), ACOX1, ACOX3 (acyl-CoA oxidases 1 and 3), APOA2, APOC3 (apolipoproteins A2 and C3), CXCL1, CXCL5, CXCL6 (C-X-C motif chemokine ligands 1, 5, and 6), FHL1 (four and a half LIM domains 1), FTH1 (ferritin heavy chain 1), GFAP (glial fibrillary acidic protein), GK (glycerol kinase), HSP90 (heat shock protein 90), HSPA5 (heat shock protein family A (Hsp70) member 5), INS (insulin), IP3R (inositol 1,4,5-trisphosphate receptor type 2), MMP1, MMP3 (matrix metallopeptidases 1 and 3), MSMO1 (methylsterol monooxygenase 1), MT1E (metallothionein 1E), NSDHL (NAD(P) dependent steroid dehydrogenase-like), PDIA4 (protein disulfide isomerase family A member 4), PLIN2 (perilipin 2), PRMT1 (protein arginine methyltransferase 1), SOD2 (superoxide dismutase 2), STAT1, STAT3 (signal transducers and activators of transcription 1 and 3), STEAP2 (STEAP2 metalloreductase), and TTR (transthyretin).
Ijms 22 07908 g008
Table
Table 1. List of proteins belonging to KEGG pathways discussed in the text and visualized in Figure 5.
Table 1. List of proteins belonging to KEGG pathways discussed in the text and visualized in Figure 5.
KEGG
Identifier		SymbolProtein NameUniProt
Identifier		Log2 Ratio
Diff./Control
Ranked list of relevant enriched pathways
hsa03320 PPAR signaling pathway (11), p-value 0.00075
hsa:123PLIN2perilipin 2Q995411.50
hsa:2181ACSL3acyl-CoA synthetase long chain family member 3O955732.17
hsa:2182ACSL4acyl-CoA synthetase long chain family member 4O604882.03
hsa:2710GKglycerol kinaseP321891.53
hsa:336APOA2apolipoprotein A2P026525.14
hsa:345APOC3apolipoprotein C3P026562.18
hsa:4312MMP1matrix metallopeptidase 1P039563.09
hsa:51ACOX1acyl-CoA oxidase 1Q150672.41
hsa:51129ANGPTL4angiopoietin like 4Q9BY761.90
hsa:81616ACSBG2acyl-CoA synthetase bubblegum family member 2Q5FVE42.13
hsa:8310ACOX3acyl-CoA oxidase 3, pristanoylO152543.11
hsa00100 Steroid biosynthesis (3), p-value 0.00593
hsa:1717DHCR77-dehydrocholesterol reductaseQ9UBM72.77
hsa:50814NSDHLNAD(P)-dependent steroid dehydrogenase-likeQ157382.47
hsa:6307MSMO1methylsterol monooxygenase 1Q158002.38
hsa04657 IL-17 signaling pathway (6), p-value 0.01139
hsa:2919CXCL1C-X-C motif chemokine ligand 1P093411.16
hsa:3320HSP90AA1heat shock protein 90 alpha family class A member 1P079001.99
hsa:4312MMP1matrix metallopeptidase 1P039563.09
hsa:4314MMP3matrix metallopeptidase 3P082544.04
hsa:6372CXCL6C-X-C motif chemokine ligand 6P801622.67
hsa:6374CXCL5C-X-C motif chemokine ligand 5P428302.04
hsa00071 Fatty acid degradation (9), p-value 0.01192
hsa:10449ACAA2acetyl-CoA acyltransferase 2P427651.06
hsa:10455ECI2enoyl-CoA delta isomerase 2O755213.70
hsa:1892ECHS1enoyl-CoA hydratase, short chain 1P300841.23
hsa:2181ACSL3acyl-CoA synthetase long chain family member 3O955732.17
hsa:2182ACSL4acyl-CoA synthetase long chain family member 4O604882.03
hsa:3030HADHAhydroxyacyl-CoA dehydrogenase trifunctional multienzyme complex sb. αP409391.16
hsa:51ACOX1acyl-CoA oxidase 1Q150672.41
hsa:81616ACSBG2acyl-CoA synthetase bubblegum family member 2Q5FVE42.13
hsa:8310ACOX3acyl-CoA oxidase 3, pristanoylO152543.11
hsa04146 Peroxisome (9), p-value 0.01345
hsa:10455ECI2enoyl-CoA delta isomerase 2O755213.70
hsa:10901DHRS4dehydrogenase/reductase 4Q9BTZ21.14
hsa:1891ECH1enoyl-CoA hydratase 1Q130111.45
hsa:2181ACSL3acyl-CoA synthetase long chain family member 3O955732.17
hsa:2182ACSL4acyl-CoA synthetase long chain family member 4O604882.03
hsa:3418IDH2isocitrate dehydrogenase (NADP(+)) 2P487351.21
hsa:51ACOX1acyl-CoA oxidase 1Q150672.41
hsa:6648SOD2superoxide dismutase 2P041793.17
hsa:8310ACOX3acyl-CoA oxidase 3, pristanoylO152543.11
hsa04068 FoxO signaling pathway (4), p-value 0.02370
hsa:3276PRMT1protein arginine methyltransferase 1Q998732.31
hsa:3630INSinsulinP013082.64
hsa:6648SOD2superoxide dismutase 2P041793.17
hsa:6774STAT3signal transducer and activator of transcription 3P40763−1.11
hsa00061 Fatty acid biosynthesis (3), p-value 0.02800
hsa:2181ACSL3acyl-CoA synthetase long chain family member 3O955732.17
hsa:2182ACSL4acyl-CoA synthetase long chain family member 4O604882.03
hsa:81616ACSBG2acyl-CoA synthetase bubblegum family member 2Q5FVE42.13
hsa05323 Rheumatoid arthritis (8), p-value 0.03344
hsa:10312TCIRG1T cell immune regulator 1, ATPase H+ transporting V0 subunit a3Q134881.25
hsa:2919CXCL1C-X-C motif chemokine ligand 1P093411.16
hsa:4312MMP1matrix metallopeptidase 1P039563.09
hsa:4314MMP3matrix metallopeptidase 3P082544.04
hsa:6372CXCL6C-X-C motif chemokine ligand 6P801622.67
hsa:6374CXCL5C-X-C motif chemokine ligand 5P428302.04
hsa:6387CXCL12C-X-C motif chemokine ligand 12P480612.41
hsa:9114ATP6V0D1ATPase H+ transporting V0 subunit d1P614211.50
hsa04918 Thyroid hormone synthesis (5), p-value 0.04442
hsa:213ALBalbuminP027682.45
hsa:3309HSPA5heat shock protein family A (Hsp70) member 5P110211.76
hsa:3709ITPR2inositol 1,4,5-trisphosphate receptor type 2Q145714.92
hsa:7276TTRtransthyretinP027662.51
hsa:9601PDIA4protein disulfide isomerase family A member 4P136671.50
hsa04630 JAK-STAT signaling pathway (4), p-value 0.04770
hsa:2273FHL1four and a half LIM domains 1Q13642−1.65
hsa:2670GFAPglial fibrillary acidic proteinP14136−1.08
hsa:6772STAT1signal transducer and activator of transcription 1P42224−1.55
hsa:6774STAT3signal transducer and activator of transcription 3P40763−1.11
hsa04978 Mineral absorption (3), p-value 0.04882
hsa:2495FTH1ferritin heavy chain 1P027941.77
hsa:261729STEAP2STEAP2 metalloreductaseQ8NFT22.60
hsa:4493MT1Emetallothionein 1EP047326.56
List of not enriched pathways with recognized function in osteodifferentiation
hsa04020 Calcium signaling pathway (9)
hsa:2767GNA11G protein subunit alpha 11P29992−3.15
hsa:291SLC25A4solute carrier family 25 member 4P122351.05
hsa:293SLC25A6solute carrier family 25 member 6P122361.05
hsa:3709ITPR2inositol 1,4,5-trisphosphate receptor type 2Q145714.92
hsa:444ASPHaspartate beta-hydroxylaseQ127971.07
hsa:7416VDAC1voltage dependent anion channel 1P217961.18
hsa:7417VDAC2voltage dependent anion channel 2P458801.14
hsa:7419VDAC3voltage dependent anion channel 3Q9Y2771.36
hsa:9630GNA14G protein subunit alpha 14O95837−1.99
hsa04310 Wnt signaling pathway (1)
hsa:7474WNT5AWnt family member 5AP412211.24
hsa04350 TGF-beta signaling pathway (2)
hsa:1634DCNdecorinP075851.75
hsa:3624INHBAinhibin subunit beta AP084765.20
hsa04010 MAPK signaling pathway (10)
hsa:1398CRKCRK proto-oncogene, adaptor proteinP46108−1.08
hsa:2768GNA12G protein subunit alpha 12Q03113−1.73
hsa:3481IGF2insulin like growth factor 2P013441.12
hsa:3630INSinsulinP013082.64
hsa:4137MAPTmicrotubule associated protein tauP10636−3.43
hsa:51776MAP3K20mitogen-activated protein kinase kinase kinase 20Q9NYL2−1.32
hsa:5536PPP5Cprotein phosphatase 5 catalytic subunitP53041−1.17
hsa:781CACNA2D1calcium voltage-gated channel auxiliary subunit alpha2delta 1P54289−1.31
hsa:929CD14CD14 moleculeP085712.42
hsa:9448MAP4K4mitogen-activated protein kinase kinase kinase kinase 4O95819−3.12
hsa04151 PI3K-Akt signaling pathway (9)
hsa:1277COL1A1collagen type I alpha 1 chainP02452−1.19
hsa:1282COL4A1collagen type IV alpha 1 chainP02462−1.66
hsa:1291COL6A1collagen type VI alpha 1 chainP12109−1.43
hsa:3320HSP90AA1heat shock protein 90 alpha family class A member 1P079001.99
hsa:3481IGF2insulin like growth factor 2P013441.12
hsa:3630INSinsulinP013082.64
hsa:3912LAMB1laminin subunit beta 1P079421.17
hsa:3913LAMB2laminin subunit beta 2P55268−1.28
hsa:5649RELNreelinP78509−1.17
hsa04714 Thermogenesis (11)
hsa:10632ATP5MGATP synthase membrane subunit gO759642.11
hsa:2181ACSL3acyl-CoA synthetase long chain family member 3O955732.17
hsa:2182ACSL4acyl-CoA synthetase long chain family member 4O604882.03
hsa:4719NDUFS1NADH:ubiquinone oxidoreductase core subunit S1P283311.02
hsa:4722NDUFS3NADH:ubiquinone oxidoreductase core subunit S3O754891.01
hsa:4729NDUFV2NADH:ubiquinone oxidoreductase core subunit V2P194041.10
hsa:57492ARID1BAT-rich interaction domain 1BQ8NFD51.20
hsa:60ACTBactin betaP607091.55
hsa:6389SDHAsuccinate dehydrogenase complex flavoprotein subunit AP310401.13
hsa:6390SDHBsuccinate dehydrogenase complex iron sulfur subunit BP219121.47
hsa:91942NDUFAF2NADH:ubiquinone oxidoreductase complex assembly factor 2Q8N1832.97
hsa00970 Aminoacyl-tRNA biosynthesis (2)
hsa:23438HARS2histidyl-tRNA synthetase 2, mitochondrialP495902.07
hsa:5859QARS1glutaminyl-tRNA synthetase 1P478971.04
hsa03050 Proteasome (2)
hsa:143471PSMA8proteasome 20S subunit alpha 8Q8TAA3−1.13
hsa:5715PSMD9proteasome 26S subunit, non-ATPase 9O00233−2.00
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nováková, S.; Danchenko, M.; Okajčeková, T.; Baranovičová, E.; Kováč, A.; Grendár, M.; Beke, G.; Pálešová, J.; Strnádel, J.; Janíčková, M.; et al. Comparative Proteomic and Metabolomic Analysis of Human Osteoblasts, Differentiated from Dental Pulp Stem Cells, Hinted Crucial Signaling Pathways Promoting Osteogenesis. Int. J. Mol. Sci. 2021, 22, 7908. https://doi.org/10.3390/ijms22157908

AMA Style

Nováková S, Danchenko M, Okajčeková T, Baranovičová E, Kováč A, Grendár M, Beke G, Pálešová J, Strnádel J, Janíčková M, et al. Comparative Proteomic and Metabolomic Analysis of Human Osteoblasts, Differentiated from Dental Pulp Stem Cells, Hinted Crucial Signaling Pathways Promoting Osteogenesis. International Journal of Molecular Sciences. 2021; 22(15):7908. https://doi.org/10.3390/ijms22157908

Chicago/Turabian Style

Nováková, Slavomíra, Maksym Danchenko, Terézia Okajčeková, Eva Baranovičová, Andrej Kováč, Marián Grendár, Gábor Beke, Janka Pálešová, Ján Strnádel, Mária Janíčková, and et al. 2021. "Comparative Proteomic and Metabolomic Analysis of Human Osteoblasts, Differentiated from Dental Pulp Stem Cells, Hinted Crucial Signaling Pathways Promoting Osteogenesis" International Journal of Molecular Sciences 22, no. 15: 7908. https://doi.org/10.3390/ijms22157908

APA Style

Nováková, S., Danchenko, M., Okajčeková, T., Baranovičová, E., Kováč, A., Grendár, M., Beke, G., Pálešová, J., Strnádel, J., Janíčková, M., Halašová, E., & Škovierová, H. (2021). Comparative Proteomic and Metabolomic Analysis of Human Osteoblasts, Differentiated from Dental Pulp Stem Cells, Hinted Crucial Signaling Pathways Promoting Osteogenesis. International Journal of Molecular Sciences, 22(15), 7908. https://doi.org/10.3390/ijms22157908

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop