Next Article in Journal
Life-Threatening Cardiogenic Shock Related to Venlafaxine Poisoning—A Case Report with Metabolomic Approach
Previous Article in Journal
Enhanced Antioxidant and Anticancer Potential of Artemisia carvifolia Buch Transformed with rol A Gene
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Machine Learning in Automated Monitoring of Metabolic Changes Accompanying the Differentiation of Adipose-Tissue-Derived Human Mesenchymal Stem Cells Employing 1H-1H TOCSY NMR

1
Image Analysis Group, TU Dortmund, 44227 Dortmund, Germany
2
Leibniz-Institut für Analytische Wissenschaften—ISAS-e.V., 44139 Dortmund, Germany
3
Cell Therapy Center, University of Jordan, Amman 11942, Jordan
4
Department of Anatomy and Histology, College of Medicine, University of Jordan, Amman 11942, Jordan
5
Nanotechnology Center, The University of Jordan, Amman 11942, Jordan
*
Author to whom correspondence should be addressed.
Metabolites 2023, 13(3), 352; https://doi.org/10.3390/metabo13030352
Submission received: 20 January 2023 / Revised: 12 February 2023 / Accepted: 23 February 2023 / Published: 27 February 2023
(This article belongs to the Section Bioinformatics and Data Analysis)

Abstract

:
The ability to monitor the dynamics of stem cell differentiation is a major goal for understanding biochemical evolution pathways. Automating the process of metabolic profiling using 2D NMR helps us to understand the various differentiation behaviors of stem cells, and therefore sheds light on the cellular pathways of development, and enhances our understanding of best practices for in vitro differentiation to guide cellular therapies. In this work, the dynamic evolution of adipose-tissue-derived human Mesenchymal stem cells (AT-derived hMSCs) after fourteen days of cultivation, adipocyte and osteocyte differentiation, was inspected based on 1H-1H TOCSY using machine learning. Multi-class classification in addition to the novelty detection of metabolites was established based on a control hMSC sample after four days’ cultivation and we successively detected the changes of metabolites in differentiated MSCs following a set of 1H-1H TOCSY experiments. The classifiers Kernel Null Foley-Sammon Transform and Kernel Density Estimation achieved a total classification error between 0% and 3.6% and false positive and false negative rates of 0%. This approach was successfully able to automatically reveal metabolic changes that accompanied MSC cellular evolution starting from their undifferentiated status to their prolonged cultivation and differentiation into adipocytes and osteocytes using machine learning supporting the research in the field of metabolic pathways of stem cell differentiation.

Graphical Abstract

1. Introduction

Mesenchymal stem cells (MSCs) are multipotent stem cells with a high capacity to proliferate and differentiate, while exhibiting low immunogenicity and providing immunosuppressive properties [1]. These potentials put MSCs in the lead as a promising candidate for several innovative strategies of cellular therapy and tissue engineering. MSCs are obtained from several body tissue, and their potential in regeneration and differentiation is highly dependent on their [2,3]. Adipose tissue is considered a highly valued source to isolate MSCs being a byproduct that generates a good yield of primary cells, with high potential to proliferate and differentiate. Therefore, adipose tissue-derived MSCs are widely used in tissue engineering and regenerative medicine [4]. Metabolic adaptation of MSCs is highly dependent on their surrounding environment; MSCs cultivated under hypoxic conditions show limited proliferation rate and high production of glycolytic enzymes, while under normoxic conditions, they show high proliferation rate and an additional reliance on oxidation phosphorylation during glycolysis; described as Warburg effect [5]. Alternatively, studies have shown that the differentiation of MSCs into osteocytes is negatively affected by normoxic conditions [6]. The switch between the glycolytic and oxidative phosphorylation pathway shows the flexibility of MSCs in adapting a metabolism that enables them to fulfil regenerative and/or immunomodulatory roles at specific sites and environments. New approaches are required to reveal novel biomarkers and information in the metabolism of MSCs and to monitor the dynamics of their metabolism in response to stimuli, and metabolic adaptation associated with several biological processes, including differentiation [7,8]. This information may unveil their behavior which would enable researchers to control and guide these cells toward successful tailor-made therapies through providing the proper culture conditions and handling [9].
Nuclear magnetic resonance (NMR) spectroscopy is a powerful technique for the identification of the components of complex mixtures consist of small molecules, like metabolites in a mixture sample. NMR has proven its vital and powerful role as an analytical technique in metabolomics extracted from biofluid, tissue extract, or semisolid samples such as intact tissues or organs [10]. Chemical shifts are fingerprints that characterize the chemical composition of a biological compound [10,11]. The non-destructiveness and the reproducibility of NMR results enables high-throughput identification and quantitative accuracy of the metabolic concentration in biological mixtures [10,12,13]. However, due to the low sensitivity of NMR, spectral resolution and spectral overlapping, obtaining the metabolic profiling data from the NMR spectra is one of the main challenges in analyzing complex biological mixtures. Overlapping of the NMR signal and the shift in the NMR peaks are affected by the pH and ionic strength variations of the biological sample in which the metabolites are measured [13,14,15]. Two-dimensional NMR (2D NMR) has a significant resolving ability through adding a second frequency domain and dispersing the peaks into this added dimension [10]. Nevertheless, metabolic profiling of a 2D NMR spectrum with low concentration or overlapped peaks is an elaborate task [16]. Moreover, the analysis of biological sample is related to the complexity of biological mixtures, phase and baseline distortion and noise [15,17]. The 2D NMR TOCSY (Total Correlation Spectroscopy) experiment provides correlations between all the protons in the spin system allowing distinguish spin systems from different molecules [10]. Still, due to the dense content of 2D NMR spectrum of complex mixtures, manual analysis is a highly demanding task and it is dependent on the researcher’s experience [18].
In this work, machine learning has been applied to automate the monitoring of the MSCs differentiation and to resolve the convolution of the associated 1H-1H TOCSY NMR spectra. The analysis is based on observing the accompanying differentiation of AT-derived hMSCs cultivated in MSCs basal culture media in addition to their adipogenic or osteogenic differentiation. Identification of compounds related to MSCs differentiation based on non-targeted metabolic profiling is a significant task and has the potential to enable effective stem cell therapy [7]. Non-targeted metabolic profiling is an all-inclusive and comprehensive analysis of the whole NMR spectrum which requires the intensive analysis of 2D NMR TOCSY profiles to reveal novel occurrences of metabolites in response to various conditions and stimuli [19,20]. Methodologically, an unbiased classification approach is mandatory to overcome variations in the biological mixtures and the corresponding complexity of the NMR-generated data. Introducing machine learning as an analysis tool for 2D NMR appears to be a reasonable approach.
Several computer implementations have been proposed to enable NMR spectral processing and cross-peak identification of 2D NMR spectra. The COLMARm [21] web server is an online available platform that incorporates three types of 2D NMR spectra for the purpose of simultaneous analysis. COLMARm operates in two stages. First, an HSQC (Heteronuclear Single Quantum Coherence) spectrum is uploaded by the user, being compared against a unified database from the Biological Magnetic Resonance Data Bank (BMRB) [22] and the Human Metabolome Database (HMAB) [11] and a matched list of metabolites is created. In the next step, the matched list is validated against the corresponding TOCSY and/or HSQC-TOCSY spectrum. A Bayesian framework is used for the problem of the assignment of peaks in 2D NMR spectra. In [23], 2D NMR spectrum were modeled as a mixture of bivariate Gaussian densities. To estimate the positions of the peaks, the adaptive Markov chain Monte Carlo (MCMC) algorithm was used. A list of candidate peaks of the highest amplitude was created and the posterior probability of each candidate peak was calculated [23]. Another peak assignment approach which incorporates the shape of the peak on the 2D spectrum was introduced in [24], where images of the peak were translated into a matrix of features through shape mapping. These features are trained and tested using a SVM classifier [24]. Neural networks have utilized in NMR for the reconstruction and denoising of spectra, chemical shift prediction and automatic peak picking [25]. Mostly, these applications are implemented using mainstream libraries, such as TensorFlow [26,27,28] or the MATLAB Deep Learning Toolbox [29]. For the purpose of chemical shift prediction, multiple types of features have been used as a feature space for the training dataset. SMART and SMART 2.0 [30,31] are based on training a deep convolutional neural network (CNN) of Siamese architecture [32] to characterize new compounds, in addition to annotate known compounds in biological mixtures. Another tool that uses CNN to analyze 2D NMR is NMRNet [33]. NMRNet starts by sorting the cross peaks according to their intensity and eliminating peaks with low intensities. The images of the selected peaks are fed as cropped segments from the spectrum to a CNN. The output of the CNN is a sigmoid function that indicates the probability of the peak assignment [33]. [34]. Stem cell osteogenic differentiation of hMSCs for 21 days based on 1D NMR has recently been studied [7,35]. They mainly considered the lipidomic and amino acid characterization of osteogenic stem cells using Principal Component Analysis (PCA) and partial least squares discriminant analysis. Human embryonic stem cells were studied to monitor the intracellular and extracellular metabolic dynamics through directed and non-directed differentiation using 1D NMR. Similarly, PCA, least squares analysis and the ANOVA test were used to compare the differentiated and undifferentiated cells [36,37].

2. Materials and Methods

2.1. Machine Learning and Novelty Detection

In this work, the procedure creates an automatic metabolic profiling system of 2D NMR TOCSY spectra based on machine learning methodologies. The Kernel Null Foley–Sammon Transform classifier (KNFST) and Kernel density estimation (KDE) were tested to monitor the dynamic evolution of adipose-tissue-derived human MSCs. Novelty Detection (ND) or outlier detection is the task of distinguishing new samples that differ from the data on which the classifier has been trained. ND is involved in applications where new categories or classes are expected to appear in the future. Although the practicable training dataset is complete and contains all classification information at a given time, it is inapplicable to encompass all variations of the possible classes that might be encountered. Therefore, the classifier is supposed to detect new extreme conditions rather than classifying them into already-available classes. In the former situation, the training model developed during the training phase is not representative of the actual classification problem and the domain of expected categories. Extreme situations could be the emergence of unknown or unexpected metabolites during the dynamic biological evolution of samples, fault detection in industrial systems or the detection of new words in hand-writing applications [38,39,40]. A vital principle in ND is the novelty threshold which acts as a discrimination criterion between recognized and novel categories. The novelty threshold defines a decision function that is applied on the output scores of the classifiers [38,40,41].
The KNFST classifier is based on achieving the best separability between classes by maximizing the between-class scatter S b and minimizing the within-class scatter S w in a high-dimensional space φ . The objective of KNFST is to learn a discriminate projection direction ω which is calculated under the conditions ω T S w φ ω = 0 and   ω T S b φ ω > 0 . The previous conditions assure the best discrimination between classes [42,43,44,45]. Class identification can be achieved by calculating the Euclidean distance between the center of the projected training instances and the projection of the tested sample. KDE is a non-parametric probability-based classifier which measures the density of independent and identically distributed random points within a neighborhood range defined by the neighborhood width h. The Parzen window estimator [46] is a common KDE approach in which the class-conditional probability p ( x ) is calculated as a linear combination of the neighboring kernels k at each point in the dataset x . Parzen estimators can be defined as p ( x ) = 1 N i = 1 N 1 h d k ( x x i h ) , where N and d define the number of samples x in the training data and the dimension of the feature space, respectively. The bandwidth h defines the smoothness of the kernel function known as the Parzen window [47]. KDE and KNFST have previously been used in the metabolic profiling of 2D TOCSY spectra of breast cancer tissue samples [45,48].

2.2. Sample Preparation and Experimentation

AT-derived hMSCs were obtained from the Cell Therapy Center (CTC)/The University of Jordan. The samples belong to consenting healthy females in the age range of 35–43, and the donors’ recruitment and sample collection were approved by the Institutional Review Board, the University of Jordan (IRB: CTC/1-2020/04 and approved on 10 March 2020).

2.2.1. Cultivation of AT-Derived hMSCs

MSCs were maintained in basal MSCs culture media composed of alpha MEM medium with Earle’s Salts (Gibco) supplemented with 5% human platelet lysate (hPL), at a concentration of 3 I.U Heparin-Sodium 5000 I.U/mL, 1% penicillin streptomycin and 2 mM L-glutamine [49]. The cells were cultured in an adherent plate at a seeding density of 4000 cells/cm2, and subculture was performed every time the cells reached a confluence of 80% until reaching cell division in passage number 4 (P4). The passage number indicates the number of times that cells have been collected and recultured into new cell culture flasks [50].

2.2.2. Adipogenic and Osteogenic Differentiation of AT-Derived hMSCs

AT-MSCs were induced to differentiate into adipocytes or osteocytes using StemPro Adipogenesis and the Osteogenesis Differentiation Kit (Gibco), respectively, as described by the manufacturer. In brief, MSCs at P4 were cultivated in MSCs basal culture media (BCM) at a seeding density of 4000 cells/cm2. When cells reached 70% confluence, basal culture media was aspirated and the cells were washed twice with PBS, before the addition of complete adipogenic (ADM) or osteogenic differentiation media (ODM). The cells were maintained in standard culture conditions (37 °C, 5% CO2) in a humidified incubator for 14 days, while refeeding the cells every 3–4 days with completely fresh media. Through the duration of differentiation, morphological changes in MSCs were monitored using inverted microscopy. To confirm the differentiation of MSCs into adipocytes and osteocytes at the end of the differentiation duration, the generated monolayer of adipogenic- or osteogenic-induced MSCs went through a staining procedure using oil red O for adipocytes, or Alizarin red staining for osteocytes [51]. Oil red staining illustrates the internal neutral lipids generated in adipocytes [52,53], whereas alizarin red staining illustrates mineral deposits, such as calcium, generated by osteocytes [54]. BCM is supposed to maintain the stemness of MSCs without triggering their differentiation, and this was confirmed by the lack of coloration in AT-derived MSCs after 4 days of cultivation, as seen in Figure 1a. However, prolonged culture duration of AT-derived MSCs triggered their differentiation even in BMC which is detectable through the formation of lipid droplets and a faded oil red staining on Figure 1b. As seen in Figure 1c, MSCs cultivated in ADM for 14 days showed a clear alteration in their morphology due to the formation of large oil droplets in their cytoplasm as presented by the intense oil red staining. On the other hand, MSCs cultivated in ODM exhibited an intense deposition of minerals represented by the intense alizarin red staining, as shown on Figure 1d.

2.2.3. Intracellular Metabolites Extraction

Following differentiation, intracellular metabolites were extracted using the methanol extraction method, as previously described [55]. Briefly, differentiation media were aspirated, and the cultured cells were washed three times with phosphate-buffered saline (PBS). Immediately after washing, absolute methanol stored at −20 °C and water ice were added to the cells in a ratio of 2 parts:0.8 parts MeOH:H2O to quench metabolism. Culture plates were stored at −80 °C for 10 min, then, the cells were scraped off the cell culture plate, and the obtained cells/methanol mixture were centrifuged at a speed of 14,000 rpm for 10 min. To obtain the intracellular metabolite in powder form, the samples were lyophilized, and the obtained powder from each sample was stored at −80 °C until further use [51].

2.3. High Resolution 1D and 2D NMR Experiments

The NMR measurements were performed at Leibniz Institute for Analytical Sciences—ISAS, Dortmund, Germany. For 1H NMR profiling, 600 μL of deuterium oxide (D2O) (Sigma-Aldrich, Taufkirchen, Germany) was added to the lyophilized metabolite, in addition to an appropriate concentration of 3-(trimethylsilyl) propionate-2,2,3,3-d4 (TSP) as an internal reference and mixed thoroughly. Later, the samples were transferred into high resolution 5 mm borosilicate glass NMR tubes (Boro-600-5-8, Deutero GmbH, Kastellaun, Germany). The high resolution 1H NMR spectra of the intracellular extracted samples in addition to two reference samples were acquired using a broadband high resolution 600.13 MHz (B0 = 14.1 T) NMR Bruker spectrometer (Avance III 600, Bruker BioSpin GmbH, Rheinstetten, Germany) and a room-temperature NMR probe (BBO model-Bruker) at 279 K. Acquisition and processing of NMR spectra was achieved by using the Bruker TopSpin 3.6. The 1D NMR spectra were acquired using the 90° single-pulse experiment (Bruker pulse sequence zg) with embedded excitation sculpting for water suppression. 1H-1H TOCSY was acquired by employing the phase-sensitive TOCSY experiment, using z-axis decoupling in the presence of scalar interactions (DIPSI)-2 spin-lock implemented in the Bruker pulse sequence dipsi2esgpph. The spectral range was set to 7 kHz in both dimensions, 16K and 128 data points acquired in the horizontal and the vertical dimension (F2, F1), respectively. Before 2D Fourier Transform, zero filling was performed to 32K and 1K data points in the horizontal and the vertical dimension, respectively. The spectral widths in the two dimensions were 12.00 ppm.

2.4. Metabolic Profiling Assignment

The metabolites were initially assigned using the high resolution 1D 1H NMR spectra of the studied MSCs in this work. One-dimensional 1H NMR high resolution spectra processing and pre-analysis were achieved using the TopSpin3.6, and metabolic assignment was accomplished using BMRB [22], HMAB [11] and Chenomx NMR Analysis Software. The detected metabolites were identified and annotated in the 1D spectra as shown in Figure 2.
The initial processing of the 2D NMR spectra was conducted using TopSpin 3.6 as follows: the spectra were referenced to the 2D contour of TSP and base levels were equalized to eliminate background noise. Later, automated peak picking at a proper threshold was performed by applying the automatic method using the pp2 function in TopSpin 3.6, and then the obtained F2 and F1 frequencies were determined.
In agreement with the 1D spectra, a total of 32 metabolites were assigned from the 2D NMR spectra as shown in Table 1. Metabolites with only a single signal do not appear in the TOCSY spectrum. Taurine and asparagine were only detectable in the 2D spectra because they fully overlapped in the 1D spectra. Table 1 contains the metabolite name in the first column, and F2 and F1 measured the frequencies of the “Ct d4”, “Ct d14”, “AT d14” and “OS d14” samples. In the last column of Table 1, the standard F2 and F1 frequencies are listed. It can be observed that some metabolites appear and disappear during the cultivation and differentiation of the cells. In Table 1, the abbreviation ‘NP’ stands for ‘not present’ and it exposes the disappearance of metabolites during the dynamic evolution of the cells. Looking at the obtained metabolic 1D and 2D NMR spectra, metabolic changes occurring in the MSCs in response to prolonged cultivation or differentiation are noticeable and are mainly presented in their lipid profiles; this was shown by the different chemical groups corresponding to fatty acids. Multiple peaks, corresponding to the presence of chemical groups related to fatty acids that are normally produced by adipocytes, were predominant in the 1D and 2D NMR spectra of the differentiated and prolonged cultivation. MSCs differentiation is strongly related to remodeling in lipidomic metabolism directed by a variation in membrane demands depending on the differentiation characteristics and functional phenotypes [7,56,57,58]. Due to the variation of the level of intracellular metabolites, equalizing the signal intensities between all TOSCY NMR spectra leads to the disappearance of peaks with a signal to noise ratio (SNR) of less than three, as shown in Figure 2. A schematic diagram of the experimental results of this work is shown in Figure 3. AT-derived hMSCs are cultivated in a basal culture media and measured after four days using NMR. Non-targeted metabolic profiling of 2D NMR TOCSY is generated based on the four days’ cultivation where all collected peaks are manually assigned by the expert. AT-MSCs were subdivided into three experiments. In the first one, the MSCs were maintained in basal MSCs culture for prolonged cultivation. In the second and third experiments, AT-MSCs were induced to differentiate into adipocytes or osteocytes, respectively. After fourteen days, the adipogenic and osteogenic differentiation of the AT-derived hMSCs in addition to their control group were measured using 2D NMR TOCSY. Similarly, peak-picking was applied and the cross peaks were assigned by an expert. To evaluate the performance of our methodology, the manual assignments were compared to the automated method.

3. Datasets

In machine learning, creating a training model using diverse and large training dataset is crucial. Nevertheless, a reliable, large and labeled data set which considers the chemical shift and peak overlap does not exist. Using data augmentation [45,59], an extended data set of the peaks corresponding to the metabolites appearing in Table 1 is created. Multiple versions of the same metabolite are created by shifting the experimental chemical shift right and left up to 30 Hz to create the training dataset and adding random Gaussian noise to create the validation dataset [45,60]. Data augmentation is applied on the “control group at 4 days cultivation (Ct d4)” to create the training dataset.
The training dataset consists of 4000 independent data instances comprising all metabolites found on “Ct d4”. The horizontal and vertical frequencies of the TOCSY spectrum represent the features of the metabolites and the corresponding multiplet.
Due to the different number of multiples per metabolite, an uneven distribution of classes in the training dataset is observed and a class imbalance problem can arise. To overcome this issue, under-sampling of metabolites with more than two multiples has been applied during the data augmentation procedure. Figure 4 shows the feature space of the metabolites contained in the cross peaks of the metabolites contained in the samples Ct d4, Ct d14 (control group at 14 days of cultivation), AT d14 (after 14 days of adipogenic differentiation) and OS d14 (after 14 days of osteogenic differentiation). It can be observed that the peaks overlap on the horizontal and vertical axes and cannot be linearly separated.

4. Results and Discussion of the Metabolic Evolution of AT-Derived hMSCs

To observe the dynamic of the AT-derived hMSCs at 14 days of cultivation (Ct d14), adipocytes (AT d14) and osteocytes (OS d14) after 14 days of differentiation, the training dataset created from (Ct d4) is used to create the main training model θ Ct   d 4 using KNFST and KDE. Three independent testing datasets are constructed using Ct d14, AT d14 and OS d14 using the corresponding frequencies in Table 1, and are introduced to the classifiers and tested against θ Ct   d 4 .
The results are reported as multi-class confusion matrices that compare the human-based metabolic profiling with the predicted assignments of the frequencies of the TOCSY spectra. In addition, Figure 5 shows the novelty scores produced by the classifiers are plotted to show the separation ability of the classifier in terms of projection distance for KNFST and probability estimation for KDE. The scores are color-coded to distinguish the scores of the different representations of classifier outputs as follows: the scores of known instances in the training set in blue, the scores of known instances in the testing dataset in green, the scores of missed novel instances in pink, the scores of correctly classified novel classes in red and the scores of misclassified known instances in the testing dataset in black. In ideal cases, the scores of known classes in the training dataset and testing dataset are similar. On the other hand, the scores of novel instances must be relatively different to those known classes. Novelty thresholds are created based on the validating dataset choosing the thresholds with a minimum validation error.
Ct d14: Figure 6 shows the confusion matrices for the output of the classifier KNFST and KDE for the Ct d14 sample. Both classifiers were able to detect all the sixteen novel frequencies which belong the fatty acids, 1-methylnicotinamide, myo-inositol, and taurine in the sample. No misclassification was encountered in KDE. This can be observed in Figure 5b, where the output of the known testing data, training data and novel classes are clearly distinct. Nevertheless, KNFST had two misclassifications within known classes, where the two instances of valine were misclassified as proline. This can be seen in Figure 5a, where two instances were plotted in pink, indicating the misclassification within known classes.
AT d14: It can be seen in Figure 7, that both classifiers predicted all the sixteen novel metabolites which belong to the fatty acids, 1-methylnicotinamide, myo-inositol, and taurine in the sample. Nevertheless, both classifiers had misclassification within already known classes. KNFST and KDE misclassified methionine as glutamine. In addition, KNFST misclassified one of the instances of valine and proline as well as misclassified one instance of leucine as threonine. This can also be seen in Figure 5c,d, where misclassifications of known classes were plotted in pink.
OS d14: Figure 8 shows the confusion matrices for the output of the classifier KNFST and KDE for the OS d14 sample. Both classifiers were able to detect all six novel instances in the sample, such as myo-inositol, Fat2 and taurine. However, it can be observed that valine was misclassified as proline in KDE. This may be due to the overlap in the vertical and horizontal frequencies between these metabolites, which can be seen in Table 1 and Figure 4d. Except for this single misclassification, no misclassification was encountered in both classifiers. This can be also observed in Figure 5e,f.
Depending on the test sample, the number and type of novel metabolites differ. For instance, there are 16 identical novel (but shifted in frequency) metabolites in Ct d14 and AT d14 in comparison to Ct d4. Nevertheless, the disappearance of metabolites in both of these samples is also different. In sample OS d14, six metabolites were found in comparison to Ct d4, and more metabolites disappeared during the differentiation. For both classifiers and all samples, the disappearance of metabolites during the biological pathway did not affect the classification performance. For instance, though the main training model θ Ct   d 4 was created on specific metabolites that disappeared in the spectra of Ct d14, AT d14 and OS d14, both classifiers proved their classification flexibility in observing metabolites presence and absence. Hence, the classifiers were able to detect both the presence and the absence of individual metabolites in accordance with θ Ct   d 4 .
Following the novelty detection metrics used in [61], the assessment measures are:
False   negative   rate   = ( 100 * F n ) / ( N n ) False   positive   rate   = ( 100 * F p ) / ( N N n ) Total   error = ( 100 * E r r o r ) / ( N )
where F n is the number of novel metabolites classified as known, N n is the number of novel instances in the test dataset, N is the total number of instances in the test dataset, F p is the number of known metabolites misclassified as novel metabolites and F e is the misclassifications within known metabolites and E r r o r = F n + F e + F p . Table 2 shows the results following these assessment measures. It can be seen that no false positive or false negative error was encountered. However, the most recurrent error was related to the misclassification within known classes which can be associated with the overlap in the frequency between these metabolites.

5. Conclusions

This article demonstrates using machine learning to perform an automatic analysis of 1H-1H TOCSY spectra acquired on cultivated and differentiated adipose-tissue-derived human MSCs (AT-derived hMSCs). Multi-class classification in addition to the novelty detection of metabolites were established based on four different 2D NMR TOCSY spectra. The primary training model was built using TOCSY spectrum of AT-derived hMSCs at four days of cultivation. Subsequently, the metabolic changes of AT-derived hMSCs control sample was monitored under three different biological settings employing the classifiers KDE and KNFST. In spite of the severe overlapping in the frequencies in TOCSY spectra, the classification outputs proved the efficiency of the used method. KDE and KNFST achieved a total classification error between 0% and 3.6% and false positive and false negative rates of 0%. The investigation in this work confirms the common metabolic pathways associated to stem cell biology. In the future, further features can be added to the dataset to produce a higher discriminative ability. Furthermore, chemical structure information or integrating other 2D NMR spectra can be included in the classification process. This work provides methodological approaches to track information of MSCs metabolism and their biological pathways, including detecting novel metabolites related to diverse stimuli in terms of prolonged cultivation and varied differentiation. This work can be extended to monitor further kinds of MSCs proliferation and recognize spectral signatures of pathways and processes.

Author Contributions

Conceptualization, L.M. and A.T.; methodology, L.M., A.T. and N.S.; software, L.M.; validation, L.M., A.T. and N.S.; formal analysis, L.M., A.T. and N.S.; investigation, L.M.; resources, L.M., A.T., H.J. and N.S.; data curation, L.M., A.T. and N.S.; writing—original draft preparation, L.M., A.T. and N.S.; writing—review and editing, L.M., A.T., R.H., C.W. and N.S.; visualization, L.M. and A.T.; supervision, A.T., R.H., H.J. and C.W.; project administration, L.M., R.H. and C.W. All authors have read and agreed to the published version of the manuscript.

Funding

Financial support by the Ministerium für Innovation, Wissenschaft und Forschung des Landes Nordrhein-Westfalen, the Senatsverwaltung für Wirtschaft, Technologie und Forschung des Landes Berlin, the Bundesministerium für Bildung und Forschung and the German Academic Exchange Service (DAAD Project no. 57587918) is acknowledged.

Institutional Review Board Statement

This research was approved by the institutional review board by the Institutional Review Board from the Cell Therapy Center/the University of Jordan (IRB: CTC/1-2020/04 and approved on 10 March 2020).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data available in a publicly accessible repository. The data presented in this study are openly available on Zenodo at: https://doi.org/10.5281/zenodo.7276518. (Accessed on 24 February 2023).

Conflicts of Interest

The authors declare no competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Saeedi, P.; Halabian, R.; Fooladi, A.A.I. A revealing review of mesenchymal stem cells therapy, clinical perspectives and Modification strategies. Stem Cell Investig. 2019, 6, 34. [Google Scholar] [CrossRef]
  2. Andrzejewska, A.; Lukomska, B.; Janowski, M. Concise review: Mesenchymal stem cells: From roots to boost. Stem Cells 2019, 37, 855–864. [Google Scholar] [CrossRef] [Green Version]
  3. Alhattab, D.; Jamali, F.; Ali, D.; Hammad, H.; Adwan, S.; Rahmeh, R.; Samarah, O.; Salah, B.; Hamdan, M.; Awidi, A. An insight into the whole transcriptome profile of four tissue-specific human mesenchymal stem cells. Regen. Med. 2019, 14, 841–865. [Google Scholar] [CrossRef]
  4. Chu, D.-T.; Nguyen Thi Phuong, T.; Tien, N.L.B.; Tran, D.K.; Minh, L.B.; Van Thanh, V.; Gia Anh, P.; Pham, V.H.; Thi Nga, V. Adipose Tissue Stem Cells for Therapy: An Update on the Progress of Isolation, Culture, Storage, and Clinical Application. J. Clin. Med. 2019, 8, 917. [Google Scholar] [CrossRef] [Green Version]
  5. Funes, J.M.; Quintero, M.; Henderson, S.; Martinez, D.; Qureshi, U.; Westwood, C.; Clements, M.O.; Bourboulia, D.; Pedley, R.B.; Moncada, S.; et al. Transformation of human mesenchymal stem cells increases their dependency on oxidative phosphorylation for energy production. Proc. Natl. Acad. Sci. USA 2007, 104, 6223–6228. [Google Scholar] [CrossRef] [Green Version]
  6. Rocha, B.; Calamia, V.; Mateos, J.; Fernández-Puente, P.; Blanco, F.J.; Ruiz-Romero, C. Metabolic Labeling of Human Bone Marrow Mesenchymal Stem Cells for the Quantitative Analysis of their Chondrogenic Differentiation. J. Proteome Res. 2012, 11, 5350–5361. [Google Scholar] [CrossRef]
  7. Bispo, D.S.C.; Jesus, C.S.H.; Correia, M.; Ferreira, F.; Bonifazio, G.; Goodfellow, B.J.; Oliveira, M.B.; Mano, J.F.; Gil, A.M. NMR Metabolomics Assessment of Osteogenic Differentiation of Adipose-Tissue-Derived Mesenchymal Stem Cells. J. Proteome Res. 2022, 21, 654–670. [Google Scholar] [CrossRef]
  8. Salazar-Noratto, G.E.; Luo, G.; Denoeud, C.; Padrona, M.; Moya, A.; Bensidhoum, M.; Bizios, R.; Potier, E.; Logeart-Avramoglou, D.; Petite, H. Understanding and leveraging cell metabolism to enhance mesenchymal stem cell transplantation survival in tissue engineering and regenerative medicine applications. Stem Cells 2019, 38, 22–33. [Google Scholar] [CrossRef] [Green Version]
  9. Zhu, H.; Sun, A.; Zou, Y.; Ge, J. Inducible Metabolic Adaptation Promotes Mesenchymal Stem Cell Therapy for Ischemia. Arter. Thromb. Vasc. Biol. 2014, 34, 870–876. [Google Scholar] [CrossRef] [Green Version]
  10. Emwas, A.-H.; Roy, R.; McKay, R.T.; Tenori, L.; Saccenti, E.; Gowda, G.A.N.; Raftery, D.; Alahmari, F.; Jaremko, L.; Jaremko, M.; et al. NMR Spectroscopy for Metabolomics Research. Metabolites 2019, 9, 123. [Google Scholar] [CrossRef] [Green Version]
  11. Wishart, D.S.; Jewison, T.; Guo, A.C.; Wilson, M.; Knox, C.; Liu, Y.; Djoumbou, Y.; Mandal, R.; Aziat, F.; Dong, E.; et al. HMDB 3.0—The Human Metabolome Database in 2013. Nucleic Acids Res. 2013, 41, D801–D807. [Google Scholar] [CrossRef]
  12. Xu, Z.-F.; Pan, A.-Z.; Yong, F.; Shen, C.-Y.; Chen, Y.-W.; Wu, R.-H. Human umbilical mesenchymal stem cell and its adipogenic differentiation: Profiling by nuclear magnetic resonance spectroscopy. World J. Stem Cells 2012, 4, 21–27. [Google Scholar] [CrossRef]
  13. Gowda, G.A.N.; Zhang, S.; Gu, H.; Asiago, V.; Shanaiah, N.; Raftery, D. Metabolomics-based methods for early disease diagnostics. Expert Rev. Mol. Diagn. 2008, 8, 617–633. [Google Scholar] [CrossRef] [Green Version]
  14. Bingol, K.; Zhang, F.; Bruschweiler-Li, L.; Brüschweiler, R. Quantitative Analysis of Metabolic Mixtures by Two-Dimensional 13C Constant-Time TOCSY NMR Spectroscopy. Anal. Chem. 2013, 85, 6414–6420. [Google Scholar] [CrossRef] [Green Version]
  15. Hao, J.; Liebeke, M.; Astle, W.; De Iorio, M.; Bundy, J.G.; Ebbels, T. Bayesian deconvolution and quantification of metabolites in complex 1D NMR spectra using BATMAN. Nat. Protoc. 2014, 9, 1416–1427. [Google Scholar] [CrossRef]
  16. Dona, A.C.; Kyriakides, M.; Scott, F.; Shephard, E.A.; Varshavi, D.; Veselkov, K.; Everett, J.R. A guide to the identification of metabolites in NMR-based metabonomics/metabolomics experiments. Comput. Struct. Biotechnol. J. 2016, 14, 135–153. [Google Scholar] [CrossRef] [Green Version]
  17. Güntert, P. Automated structure determination from NMR spectra. Eur. Biophys. J. 2008, 38, 129–143. [Google Scholar] [CrossRef]
  18. Ross, A.; Schlotterbeck, G.; Dieterle, F.; Senn, H. Chapter 3—NMR Spectroscopy Techniques for Application to Metabonomics. In The Handbook of Metabonomics and Metabolomics; Lindon, J., Nicholson, J., Holmes, E., Eds.; Elsevier Science B.V.: Amsterdam, The Netherlands, 2007; pp. 55–112. [Google Scholar]
  19. Roberts, L.D.; Souza, A.L.; Gerszten, R.E.; Clish, C.B. Targeted Metabolomics. Curr. Protoc. Mol. Biol. 2012, 98, 30.2.1–30.2.24. [Google Scholar] [CrossRef]
  20. Ai, Z.; Zhang, Y.; Li, X.; Sun, W.; Liu, Y. Widely Targeted Metabolomics Analysis to Reveal Transformation Mechanism of Cistanche Deserticola Active Compounds During Steaming and Drying Processes. Front. Nutr. 2021, 8, 743. [Google Scholar] [CrossRef]
  21. Bingol, K.; Li, D.-W.; Zhang, B.; Brüschweiler, R. Comprehensive Metabolite Identification Strategy Using Multiple Two-Dimensional NMR Spectra of a Complex Mixture Implemented in the COLMARm Web Server. Anal. Chem. 2016, 88, 12411–12418. [Google Scholar] [CrossRef] [Green Version]
  22. Ulrich, E.L.; Akutsu, H.; Doreleijers, J.F.; Harano, Y.; Ioannidis, Y.E.; Lin, J.; Livny, M.; Mading, S.; Maziuk, D.; Miller, Z.; et al. BioMagResBank. Nucleic Acids Res. 2007, 36, D402–D408. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Cheng, Y.; Gao, X.; Liang, F. Bayesian Peak Picking for NMR Spectra. Genom. Proteom. Bioinform. 2014, 12, 39–47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Klukowski, P.; Walczak, M.J.; Gonczarek, A.; Boudet, J.; Wider, G. Computer vision-based automated peak picking applied to protein NMR spectra. Bioinformatics 2015, 31, 2981–2988. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Chen, D.; Wang, Z.; Guo, D.; Orekhov, V.; Qu, X. Review and Prospect: Deep Learning in Nuclear Magnetic Resonance Spectroscopy. Chem. A Eur. J. 2020, 26, 10391–10401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Abadi, M.; Agarwal, A.; Barham, P.; Brevdo, E.; Chen, Z.; Citro, C.; Corrado, G.S.; Davis, A.; Dean, J.; Devin, M.; et al. Tensorflow: Large-scale machine learning on heterogeneous distributed systems. In 12th USENIX Symposium on Operating Systems Design and Implementation (OSDI 16); USENIX Association: Savannah, GA, USA, 2016; Available online: https://www.usenix.org/conference/osdi16/technical-sessions/presentation/abadi (accessed on 30 December 2022).
  27. Hansen, D.F. Using Deep Neural Networks to Reconstruct Non-uniformly Sampled NMR Spectra. J. Biomol. NMR 2019, 73, 577–585. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Kwon, Y.; Lee, D.; Choi, Y.-S.; Kang, M.; Kang, S. Neural Message Passing for NMR Chemical Shift Prediction. J. Chem. Inf. Model. 2020, 60, 2024–2030. [Google Scholar] [CrossRef]
  29. Lee, H.H.; Kim, H. Intact metabolite spectrum mining by deep learning in proton magnetic resonance spectroscopy of the brain. Magn. Reson. Med. 2019, 82, 33–48. [Google Scholar] [CrossRef]
  30. Zhang, C.; Idelbayev, Y.; Roberts, N.; Tao, Y.; Nannapaneni, Y.; Duggan, B.M.; Min, J.; Lin, E.C.; Gerwick, E.C.; Cottrell, G.W.; et al. Small Molecule Accurate Recognition Technology (SMART) to Enhance Natural Products Research. Sci. Rep. 2017, 7, 14243. [Google Scholar] [CrossRef] [Green Version]
  31. Reher, R.; Kim, H.W.; Zhang, C.; Mao, H.H.; Wang, M.; Nothias, L.-F.; Caraballo-Rodriguez, A.M.; Glukhov, E.; Teke, B.; Leao, T.; et al. A Convolutional Neural Network-Based Approach for the Rapid Annotation of Molecularly Diverse Natural Products. J. Am. Chem. Soc. 2020, 142, 4114–4120. [Google Scholar] [CrossRef]
  32. Chopra, S.; Hadsell, R.; LeCun, Y. Learning a similarity metric discriminatively, with application to face verification. In Proceedings of the 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR’05), San Diego, CA, USA, 20–25 June 2005. [Google Scholar]
  33. Klukowski, P.; Augoff, M.; Zięba, M.; Drwal, M.; Gonczarek, A.; Walczak, M.J. NMRNet: A deep learning approach to automated peak picking of protein NMR spectra. Bioinformatics 2018, 34, 2590–2597. [Google Scholar] [CrossRef] [Green Version]
  34. Jang, M.-Y.; Chun, S.-I.; Mun, C.-W.; Hong, K.S.; Shin, J.-W. Evaluation of Metabolomic Changes as a Biomarker of Chondrogenic Differentiation in 3D-cultured Human Mesenchymal Stem Cells Using Proton (1H) Nuclear Magnetic Resonance Spectroscopy. PLoS ONE 2013, 8, e78325. [Google Scholar] [CrossRef] [Green Version]
  35. Bispo, D.S.C.; Michálková, L.; Correia, M.; Jesus, C.S.H.; Duarte, I.F.; Goodfellow, B.J.; Oliveira, M.B.; Mano, J.F.; Gil, A.M. Endo- and Exometabolome Crosstalk in Mesenchymal Stem Cells Undergoing Osteogenic Differentiation. Cells 2022, 11, 1257. [Google Scholar] [CrossRef] [PubMed]
  36. Castiglione, F.; Ferro, M.; Mavroudakis, E.; Pellitteri, R.; Bossolasco, P.; Zaccheo, D.; Morbidelli, M.; Silani, V.; Mele, A.; Moscatelli, D.; et al. NMR Metabolomics for Stem Cell type discrimination. Sci. Rep. 2017, 7, 15808. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Coope, A.; Ghanameh, Z.; Kingston, O.; Sheridan, C.M.; Barrett-Jolley, R.; Phelan, M.M.; Oldershaw, R.A. 1H NMR Metabolite Monitoring during the Differentiation of Human Induced Pluripotent Stem Cells Provides New Insights into the Molecular Events That Regulate Embryonic Chondrogenesis. Int. J. Mol. Sci. 2022, 23, 9266. [Google Scholar] [CrossRef] [PubMed]
  38. Roberts, S.J. Extreme value statistics for novelty detection in biomedical signal processing. In Proceedings of the 2000 First International Conference Advances in Medical Signal and Information Processing, Bristol, UK, 4–6 September 2000; pp. 168–172. [Google Scholar]
  39. Markou, M.; Singh, S. Novelty detection: A review—part 1: Statistical approaches. Signal Process. 2003, 83, 2481–2497. [Google Scholar] [CrossRef]
  40. Pimentel, M.A.; Clifton, D.A.; Clifton, L.; Tarassenko, L. A review of novelty detection. Signal Process. 2014, 99, 215–249. [Google Scholar] [CrossRef]
  41. Bishop, C. Novelty detection and neural network validation. IEE Proc.—Vision, Image, Signal Process. 1994, 141, 217–222. [Google Scholar] [CrossRef]
  42. Bodesheim, P.; Freytag, A.; Rodner, E.; Kemmler, M.; Denzler, J. Kernel Null Space Methods for Novelty Detection. In Proceedings of the IEEE Computer Society Conference on Computer Vision and Pattern Recognition, Portland, OR, USA, 23–28 June 2013. [Google Scholar]
  43. Huang, X.; Xu, J.; Guo, G. Incremental Kernel Null Foley-Sammon Transform for Person Re-identification. In Proceedings of the 2018 24th International Conference on Pattern Recognition (ICPR), Beijing, China, 20–24 August 2018; pp. 1683–1688. [Google Scholar] [CrossRef]
  44. Gu, G.; Liu, H.; Shen, J. Kernel null foley-sammon transform. In Proceedings of the 2008 International Conference on Computer Science and Software Engineering, Washington, DC, USA, 12–14 December 2008. [Google Scholar]
  45. Migdadi, L.; Lambert, J.; Telfah, A.; Hergenröder, R.; Wöhler, C. Automated metabolic assignment: Semi-supervised learning in metabolic analysis employing two dimensional Nuclear Magnetic Resonance (NMR). Comput. Struct. Biotechnol. J. 2021, 19, 5047–5058. [Google Scholar] [CrossRef]
  46. Parzen, E. On Estimation of a Probability Density Function and Mode. Ann. Math. Stat. 1962, 33, 1065–1076. [Google Scholar] [CrossRef]
  47. Bishop, C.M. Pattern Recognition and Machine Learning 2006; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  48. Migdadi, L.; Telfah, A.; Hergenröder, R.; Wöhler, C. Novelty detection for metabolic dynamics established on breast cancer tissue using 2D NMR TOCSY spectra. Comput. Struct. Biotechnol. J. 2022, 20, 2965–2977. [Google Scholar] [CrossRef]
  49. Abuarqoub, D.; Awidi, A.; Abuharfeil, N. Comparison of osteo/odontogenic differentiation of human adult dental pulp stem cells and stem cells from apical papilla in the presence of platelet lysate. Arch. Oral Biol. 2015, 60, 1545–1553. [Google Scholar] [CrossRef] [PubMed]
  50. Yoon, J.-Y. Cell Culture. In Tissue Engineering: A Primer with Laboratory Demonstrations; Springer: Berlin/Heidelberg, Germany, 2022. [Google Scholar]
  51. Ghorbani, A.; Jalali, S.A.; Varedi, M. Isolation of adipose tissue mesenchymal stem cells without tissue destruction: A non-enzymatic method. Tissue Cell 2014, 46, 54–58. [Google Scholar] [CrossRef] [PubMed]
  52. Sathishkumar, S.; Mohanashankar, P.; Boopalan, P. Cell surface protein expression of stem cells from human adipose tissue at early passage with reference to mesenchymal stem cell phenotype. Int. J. Med. Med. Sci. 2011, 3, 129–134. [Google Scholar]
  53. Zhang, A.-X.; Yu, W.-H.; Ma, B.-F.; Yu, X.-B.; Mao, F.F.; Liu, W.; Zhang, J.-Q.; Zhang, X.-M.; Li, S.-N.; Li, M.-T.; et al. Proteomic identification of differently expressed proteins responsible for osteoblast differentiation from human mesenchymal stem cells. Mol. Cell. Biochem. 2007, 304, 167–179. [Google Scholar] [CrossRef] [PubMed]
  54. Umrath, F.; Weber, M.; Reinert, S.; Wendel, H.-P.; Avci-Adali, M.; Alexander, D. iPSC-Derived MSCs versus Originating Jaw Periosteal Cells: Comparison of Resulting Phenotype and Stem Cell Potential. Int. J. Mol. Sci. 2020, 21, 587. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Martineau, E.; Tea, I.; Loaëc, G.; Giraudeau, P.; Akoka, S. Strategy for choosing extraction procedures for NMR-based metabolomic analysis of mammalian cells. Anal. Bioanal. Chem. 2011, 401, 2133–2142. [Google Scholar] [CrossRef] [PubMed]
  56. Levental, K.R.; Surma, M.A.; Skinkle, A.D.; Lorent, J.H.; Zhou, Y.; Klose, C.; Chang, J.T.; Hancock, J.F.; Levental, I. ω-3 polyunsaturated fatty acids direct differentiation of the membrane phenotype in mesenchymal stem cells to potentiate osteogenesis. Sci. Adv. 2017, 3, eaao1193. [Google Scholar] [CrossRef] [Green Version]
  57. Assis-Ribas, T.; Forni, M.F.; Winnischofer, S.M.B.; Sogayar, M.; Trombetta-Lima, M. Extracellular matrix dynamics during mesenchymal stem cells differentiation. Dev. Biol. 2018, 437, 63–74. [Google Scholar] [CrossRef] [PubMed]
  58. Shi, C.; Wang, X.; Wu, S.; Zhu, Y.; Chung, L.W.K.; Mao, H. HRMAS1H-NMR measured changes of the metabolite profile as mesenchymal stem cells differentiate to targeted fat cells in vitro: Implications for non-invasive monitoring of stem cell differentiation in vivo. J. Tissue Eng. Regen. Med. 2008, 2, 482–490. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Wong, S.C.; Gatt, A.; Stamatescu, V.; McDonnell, M.D. Understanding Data Augmentation for Classification: When to Warp? In Proceedings of the 2016 International Conference on Digital Image Computing: Techniques and Applications (DICTA), Gold Coast, Australia, 30 November–2 December 2016. [Google Scholar]
  60. Liu, J.; Osadchy, M.; Ashton, L.; Foster, M.; Solomon, C.J.; Gibson, S.J. Deep convolutional neural networks for Raman spectrum recognition: A unified solution. Analyst 2017, 142, 4067–4074. [Google Scholar] [CrossRef] [Green Version]
  61. Masud, M.; Gao, J.; Khan, L.; Han, J.; Thuraisingham, B.M. Classification and Novel Class Detection in Concept-Drifting Data Streams under Time Constraints. IEEE Trans. Knowl. Data Eng. 2010, 23, 859–874. [Google Scholar] [CrossRef]
Figure 1. Light microscopy images (10x magnification, using inverted microscope, Zeiss, Oberkochen, Germany) of (a) AT-derived hMSCs after 4 days, and (b) 14 days of cultivation in basal culture media (BCM). (c) Oil red staining illustrating adipogenic differentiation of AT-derived hMSCs after 14 days of cultivation in adipogenic differentiation media (ADM). (d) Alizarin red staining illustrating osteogenic differentiation of AT-derived hMSCs after 14 days of cultivation in osteogenic differentiation media (ODM).
Figure 1. Light microscopy images (10x magnification, using inverted microscope, Zeiss, Oberkochen, Germany) of (a) AT-derived hMSCs after 4 days, and (b) 14 days of cultivation in basal culture media (BCM). (c) Oil red staining illustrating adipogenic differentiation of AT-derived hMSCs after 14 days of cultivation in adipogenic differentiation media (ADM). (d) Alizarin red staining illustrating osteogenic differentiation of AT-derived hMSCs after 14 days of cultivation in osteogenic differentiation media (ODM).
Metabolites 13 00352 g001
Figure 2. Representative high resolution 1H NMR spectra of intracellular metabolite extracts obtained from AT-derived hMSC samples collected at day 14 of differentiation into adipocytes and osteocytes, and their control samples represented in AT-derived hMSC collected at day 4 and 14 of cultivation in BCM. (a) 0.4–5 ppm region; (b) 5–10 ppm region. Peak assignment: Ile: Isoleucine; Leu: Leucine; Val: Valine; Thr: Threonine; Lac: Lactate; Ala: Alanine; Glu: Glutamine; Gln: Glutamate; Pro: Proline; Met: Methionine; Lys: Lysine; Arg: Arginine; GPC: Glycerophosphorylcholine; α-Glc: Alfa-Glucose; β-Glc: Beta-Glucose; MI: myo-inositol; ChoP: O-Phosphocholine; pEtN: Phosphorylethanolamine; GroPEtn: Glycerophosphorylethanolamine; ATP: Adenosine triphosphate; ADP: Adenosine diphosphate; Tyr: Tyrosine; Phe: Phenylalanine; NAD+: Nicotinamide adenine dinucleotide; Tau: Taurine; Asp: Asparagine; 1-MNA: 1-methylnicotinamide; AcO-: Acetate; DMA; Dimethylamine. In addition to the fatty acids signals; namely FAT 1, FAT 2, FAT 3, FAT 4, and FAT 5, representing methyl group -CH3, Acyl chains-(CH2)n-, methylene group -CH2-CH=CH, vinyl hydrogen -CH=CH, and diallyl methylene group =CH-CH2-CH=, respectively. The presence of ETOH (ethanol) and MeOH (methanol) was observed to represent residues from the cleaning and extraction procedures.
Figure 2. Representative high resolution 1H NMR spectra of intracellular metabolite extracts obtained from AT-derived hMSC samples collected at day 14 of differentiation into adipocytes and osteocytes, and their control samples represented in AT-derived hMSC collected at day 4 and 14 of cultivation in BCM. (a) 0.4–5 ppm region; (b) 5–10 ppm region. Peak assignment: Ile: Isoleucine; Leu: Leucine; Val: Valine; Thr: Threonine; Lac: Lactate; Ala: Alanine; Glu: Glutamine; Gln: Glutamate; Pro: Proline; Met: Methionine; Lys: Lysine; Arg: Arginine; GPC: Glycerophosphorylcholine; α-Glc: Alfa-Glucose; β-Glc: Beta-Glucose; MI: myo-inositol; ChoP: O-Phosphocholine; pEtN: Phosphorylethanolamine; GroPEtn: Glycerophosphorylethanolamine; ATP: Adenosine triphosphate; ADP: Adenosine diphosphate; Tyr: Tyrosine; Phe: Phenylalanine; NAD+: Nicotinamide adenine dinucleotide; Tau: Taurine; Asp: Asparagine; 1-MNA: 1-methylnicotinamide; AcO-: Acetate; DMA; Dimethylamine. In addition to the fatty acids signals; namely FAT 1, FAT 2, FAT 3, FAT 4, and FAT 5, representing methyl group -CH3, Acyl chains-(CH2)n-, methylene group -CH2-CH=CH, vinyl hydrogen -CH=CH, and diallyl methylene group =CH-CH2-CH=, respectively. The presence of ETOH (ethanol) and MeOH (methanol) was observed to represent residues from the cleaning and extraction procedures.
Metabolites 13 00352 g002aMetabolites 13 00352 g002b
Figure 3. Schematic diagram of the experimental setting to observe the metabolic evolution of AT-derived hMSCs using 2D TOCSY of intracellular extracts of MSCs cultivated in basal culture media at 4 and 14 days and MSCs cultivated for a duration of 14 days in adipogenic and osteogenic differentiation media.
Figure 3. Schematic diagram of the experimental setting to observe the metabolic evolution of AT-derived hMSCs using 2D TOCSY of intracellular extracts of MSCs cultivated in basal culture media at 4 and 14 days and MSCs cultivated for a duration of 14 days in adipogenic and osteogenic differentiation media.
Metabolites 13 00352 g003
Figure 4. Feature space of the cross peaks of the metabolites contained in the samples (a) Ct d4, (b) Ct d14, (c) AT d14 and (d) OS d14.
Figure 4. Feature space of the cross peaks of the metabolites contained in the samples (a) Ct d4, (b) Ct d14, (c) AT d14 and (d) OS d14.
Metabolites 13 00352 g004
Figure 5. Novelty scores and threshold values of KDE and KNFST classifiers for Ct d14 (a,b), AT d14 (c,d) and OS d14 (e,f).
Figure 5. Novelty scores and threshold values of KDE and KNFST classifiers for Ct d14 (a,b), AT d14 (c,d) and OS d14 (e,f).
Metabolites 13 00352 g005
Figure 6. Confusion matrices of the output of classifiers KNFST and KDE for the spectra of 14 days of cultivation.
Figure 6. Confusion matrices of the output of classifiers KNFST and KDE for the spectra of 14 days of cultivation.
Metabolites 13 00352 g006
Figure 7. Confusion matrices of the output of classifiers KNFST and KDE for the spectra of 14 days of adipocyte differentiation.
Figure 7. Confusion matrices of the output of classifiers KNFST and KDE for the spectra of 14 days of adipocyte differentiation.
Metabolites 13 00352 g007
Figure 8. Confusion matrices of the output of classifiers KNFST and KDE for the spectra of 14 days of osteocyte differentiation.
Figure 8. Confusion matrices of the output of classifiers KNFST and KDE for the spectra of 14 days of osteocyte differentiation.
Metabolites 13 00352 g008
Table 1. Intracellular metabolites detected in AT-derived hMSCs. Frequencies: control group at 4 days of cultivation (Ct d4), 14 days of cultivation (Ct d14), 14 days of differentiation into adipocytes (AT d14) and osteocytes (OS d14). NP: Not present.
Table 1. Intracellular metabolites detected in AT-derived hMSCs. Frequencies: control group at 4 days of cultivation (Ct d4), 14 days of cultivation (Ct d14), 14 days of differentiation into adipocytes (AT d14) and osteocytes (OS d14). NP: Not present.
MetaboliteCt d4Ct d14 AT d14OS d14Standard
F2F1F2F1F2F1F2F1F2F1
Leu752494762484762489772494720540
Leu950590930584930589934560900600
Leu900720914720913750910760900720
Leu10406101073631106063810466081080600
Leu10406951073705106070910467491080720
Leu10408931073923106091210469221080900
Leu22296112210608221261021986202220600
Leu21707372171750217175221637492220720
Leu21639432168943215790821759212220900
Ile10205401020540102054010205401020540
Ile10234081009403104741710454171020600
Ile21805702170580217857822205402220540
Ile2160103221651042219510522210100322201020
Tyr1920177819201787192017821920179019201830
Tyr2396178822531784235517862358178023401830
Tyr2304187722531848235618362356184823621920
Tyr4073406740903946409439614095395243164139
Phe2340185323541906226117782360184823901868
Phe2340193422541960225419262260191323901970
Phe4362419343624193436841934370419744534422
Phe4362427343624275436842734370428644534394
Glu1373110213541106135411061349110614701260
Glu2337104323371057234410482333106222581278
Glu2341127823441269234112882333127822581468
Gln1295110012811100128410711288110012601200
Gln1378122013841200138912101370123013801200
Gln2190126921861288219112882194128822201200
Gln2225137022271367221013802208136922201380
Lys1740809NPNP1736783NPNP1800840
Lys1844893NPNP1836898NPNP1800900
Lys18361062NPNP18361058NPNP18061032
Lys2295962NPNP2290962NPNP2220900
Lys22821057NPNP22781044NPNP22501032
Lys22821118NPNP22861119NPNP22501137
FAT 1NPNP616405600420NPNP600420
FAT 2NPNP789545789531789545785535
FAT 3NPNP12306141245620NPNP1260600
FAT 3NPNP1240108012601080NPNP12601050
FAT 4NPNP17157721705778NPNP1792766
FAT 5NPNP31396073150607NPNP3180540
FAT 5NPNP3138105231501052NPNP31801080
FAT 5NPNP3140121731501219NPNP31801260
Lac24997152494709249471524947202463790
Thr21607892160790216072021607202160780
Thr25787892537790258272025737202580780
Pro187912381869123018731238NPNP19801200
Pro240812462408123424081238NPNP24721213
Pro240814382435144824051439NPNP24721402
Ala22707232295696229570522917012256876
Val13506321383619139461913836191380617
Val187512371890125918801244187012402160617
Met151811871523119715181177NPNP15601260
Met233812702342127423511277NPNP23401260
Met233813702342136723511380NPNP23401320
pEtN2317185223051848233118652307184924301950
GroPEtn2300179122911781229317912292179123001791
ChoP2454195024551953245819542454194725722187
GPC2127194321211943212419392123193921601980
GPC2552233325352338254423592533234825802340
Arg114410001140975114398411469801120920
Arg19109601920960194498819289881920960
Arg1974113419781134198611151969113019201140
MINPNP20391880204418692036186920401800
MINPNP20881882209318722087186621121959
MINPNP21541970215919772152197221671959
MINPNP24602156245221402452214924232113
Asp23901574NPNP23981578NPNP24001800
Asp18701750NPNP18761768NPNP18001740
TauNPNP20641809206518122063181220401980
α-Glc3135211931252139313721323140211231302112
α-Glc3135223831252254313722633140228031302224
α-Glc3135257331252558313225623140256531302568
β-Glc2760193727741928276519312759193627781938
β-Glc2717205527142065271720632717206827782084
β-Glc2717200827142000271220022714208927782081
ATP3620258736402581NPNPNPNP36202587
ATP3620268036402628NPNPNPNP36202680
ADP3569249635662503356625013570249835692496
ADP3569270035662706356927083569269035692700
ADP3569276235662759356927693569276535692760
ADP3569288235662885356928703569286835692880
NAD+53105110NPNP53025106NPNP52005110
1-MNANPNP5218471852184725NPNP53414921
1-MNANPNP5412471854124725NPNP55814921
1-MNANPNP5520532855125321NPNP55815341
Table 2. A summary of the performance of KDE and KNFST classifiers for Ct d14, AT d14 and OS d14.
Table 2. A summary of the performance of KDE and KNFST classifiers for Ct d14, AT d14 and OS d14.
Ct d14AT d14OS d14
KNFSTKDEKNFSTKDEKNFSTKDE
False negative rate0%0%0%0%0%0%
False positive rate0%0%0%0%0%0%
Total error2.6%0%3.6%1.2%0%1.7%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Migdadi, L.; Sharar, N.; Jafar, H.; Telfah, A.; Hergenröder, R.; Wöhler, C. Machine Learning in Automated Monitoring of Metabolic Changes Accompanying the Differentiation of Adipose-Tissue-Derived Human Mesenchymal Stem Cells Employing 1H-1H TOCSY NMR. Metabolites 2023, 13, 352. https://doi.org/10.3390/metabo13030352

AMA Style

Migdadi L, Sharar N, Jafar H, Telfah A, Hergenröder R, Wöhler C. Machine Learning in Automated Monitoring of Metabolic Changes Accompanying the Differentiation of Adipose-Tissue-Derived Human Mesenchymal Stem Cells Employing 1H-1H TOCSY NMR. Metabolites. 2023; 13(3):352. https://doi.org/10.3390/metabo13030352

Chicago/Turabian Style

Migdadi, Lubaba, Nour Sharar, Hanan Jafar, Ahmad Telfah, Roland Hergenröder, and Christian Wöhler. 2023. "Machine Learning in Automated Monitoring of Metabolic Changes Accompanying the Differentiation of Adipose-Tissue-Derived Human Mesenchymal Stem Cells Employing 1H-1H TOCSY NMR" Metabolites 13, no. 3: 352. https://doi.org/10.3390/metabo13030352

APA Style

Migdadi, L., Sharar, N., Jafar, H., Telfah, A., Hergenröder, R., & Wöhler, C. (2023). Machine Learning in Automated Monitoring of Metabolic Changes Accompanying the Differentiation of Adipose-Tissue-Derived Human Mesenchymal Stem Cells Employing 1H-1H TOCSY NMR. Metabolites, 13(3), 352. https://doi.org/10.3390/metabo13030352

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