Fragment Assembly Approach Based on Graph/Network Theory with Quantum Chemistry Verifications for Assigning Multidimensional NMR Signals in Metabolite Mixtures

Metabolomic Analysis of Trehalose Alleviating Oxidative Stress in Myoblasts

Trehalose, a naturally occurring non-toxic disaccharide, has attracted considerable attention for its potential in alleviating oxidative stress in skeletal muscle. In this study, our aim was to elucidate the metabolic mechanisms underlying the protective effects of trehalose against hydrogen peroxide (H2O2)-induced oxidative stress in C2C12 myoblasts. Our results show that both trehalose treatment and pretreatment effectively alleviate the H2O2-induced decrease in cell viability, reduce intracellular reactive oxygen species (ROS), and attenuate lipid peroxidation. Furthermore, using NMR-based metabolomics analysis, we observed that trehalose treatment and pretreatment modulate the metabolic profile of myoblasts, specifically regulating oxidant metabolism and amino acid metabolism, contributing to their protective effects against oxidative stress. Importantly, our results reveal that trehalose treatment and pretreatment upregulate the expression levels of P62 and Nrf2 proteins, thereby activating the Nrf2-NQO1 axis and effectively reducing oxidative stress. These significant findings highlight the potential of trehalose supplementation as a promising and effective strategy for alleviating oxidative stress in skeletal muscle and provide valuable insights into its potential therapeutic applications.

Solubility Prediction from Molecular Properties and Analytical Data Using an In-phase Deep Neural Network (Ip-DNN)

Materials informatics is an emerging field that allows us to predict the properties of materials and has been applied in various research and development fields, such as materials science. In particular, solubility factors such as the Hansen and Hildebrand solubility parameters (HSPs and SP, respectively) and Log P are important values for understanding the physical properties of various substances. In this study, we succeeded at establishing a solubility prediction tool using a unique machine learning method called the in-phase deep neural network (ip-DNN), which starts exclusively from the analytical input data (e.g., NMR information, refractive index, and density) to predict solubility by predicting intermediate elements, such as molecular components and molecular descriptors, in the multiple-step method. For improving the level of accuracy of the prediction, intermediate regression models were employed when performing in-phase machine learning. In addition, we developed a website dedicated to the established solubility prediction method, which is freely available at "http://dmar.riken.jp/matsolca/".

Improved Prediction of Carbonless NMR Spectra by the Machine Learning of Theoretical and Fragment Descriptors for Environmental Mixture Analysis

As the first multidimensional NMR approach, 2D J-resolved (2DJ) spectroscopy is distinguished by signal resolution and detection sensitivity with remarkable advantages for the exhaustive evaluation of complex mixtures and environmental samples due to its carbonless feature without the requirement of ¹³C connectivity. Generally, the 2DJ signal assignment of metabolic mixtures is problematic in spite of references to experimental NMR databases, owing to the existence of metabolic "dark matter." In this study, a new method to predict 2DJ spectra was developed with a combination of quantum mechanical (QM) computation and machine learning (ML). The predictive accuracy of J-coupling constants was evaluated using validated data. The root-mean-square deviation (RMSD) for QM computation was 3.52 Hz, while the RMSD for QM + ML was 1.21 Hz, indicating a substantial increase in predictive accuracy. The proposed model was applied to predict the 2DJ spectra of 60 standard substances and 55 components of seawater. Furthermore, two practical environmental samples were used to evaluate the robustness of the constructed predictive model. A J-coupling tree and J-split spectra produced from QM + ML of aliphatic moieties had good consistency with the experimental data, as compared with the theoretical data produced by QM computation. The predicted J-coupling tree for the J-coupling multiplet analysis of freely rotating bonds in the complex mixture, which is traditionally difficult, was interpretable. In addition, in silico identification of the J-split ¹H NMR signals, which was independent of experimental databases, aided in the discovery of new components in a mixture.

NMR-TS: de novo molecule identification from NMR spectra

Nuclear magnetic resonance (NMR) spectroscopy is an effective tool for identifying molecules in a sample. Although many previously observed NMR spectra are accumulated in public databases, they cover only a tiny fraction of the chemical space, and molecule identification is typically accomplished manually based on expert knowledge. Herein, we propose NMR-TS, a machine-learning-based python library, to automatically identify a molecule from its NMR spectrum. NMR-TS discovers candidate molecules whose NMR spectra match the target spectrum by using deep learning and density functional theory (DFT)-computed spectra. As a proof-of-concept, we identify prototypical metabolites from their computed spectra. After an average 5451 DFT runs for each spectrum, six of the nine molecules are identified correctly, and proximal molecules are obtained in the other cases. This encouraging result implies that de novo molecule generation can contribute to the fully automated identification of chemical structures. NMR-TS is available at https://github.com/tsudalab/NMR-TS.

Spatial molecular-dynamically ordered NMR spectroscopy of intact bodies and heterogeneous systems

Noninvasive evaluation of the spatial distribution of chemical composition and diffusion behavior of materials is becoming possible by advanced nuclear magnetic resonance (NMR) pulse sequence editing. However, there is room for improvement in the spectral resolution and analytical method for application to heterogeneous samples. Here, we develop applications for comprehensively evaluating compounds and their dynamics in intact bodies and heterogeneous systems from NMR data, including spatial z-position, chemical shift, and diffusion or relaxation. This experiment is collectively named spatial molecular-dynamically ordered spectroscopy (SMOOSY). Pseudo-three-dimensional (3D) SMOOSY spectra of an intact shrimp and two heterogeneous systems are recorded to evaluate this methodology. Information about dynamics is mapped onto two-dimensional (2D) chemical shift imaging spectra using a pseudo-spectral imaging method with a processing tool named SMOOSY processor. Pseudo-2D SMOOSY spectral images can non-invasively assess the different dynamics of the compounds at each spatial z-position of the shrimp's body and two heterogeneous systems.

Exploratory machine-learned theoretical chemical shifts can closely predict metabolic mixture signals

Various chemical shift predictive methodologies have been studied and developed, but there remains the problem of prediction accuracy. Assigning the NMR signals of metabolic mixtures requires high predictive performance owing to the complexity of the signals. Here we propose a new predictive tool that combines quantum chemistry and machine learning. A scaling factor as the objective variable to correct the errors of 2355 theoretical chemical shifts was optimized by exploring 91 machine learning algorithms and using the partial structure of 150 compounds as explanatory variables. The optimal predictive model gave RMSDs between experimental and predicted chemical shifts of 0.2177 ppm for δ 1H and 3.3261 ppm for δ 13C in the test data; thus, better accuracy was achieved compared with existing empirical and quantum chemical methods. The utility of the predictive model was demonstrated by applying it to assignments of experimental NMR signals of a complex metabolic mixture.

Application of kernel principal component analysis and computational machine learning to exploration of metabolites strongly associated with diet

Computer-based technological innovation provides advancements in sophisticated and diverse analytical instruments, enabling massive amounts of data collection with relative ease. This is accompanied by a fast-growing demand for technological progress in data mining methods for analysis of big data derived from chemical and biological systems. From this perspective, use of a general “linear” multivariate analysis alone limits interpretations due to “non-linear” variations in metabolic data from living organisms. Here we describe a kernel principal component analysis (KPCA)-incorporated analytical approach for extracting useful information from metabolic profiling data. To overcome the limitation of important variable (metabolite) determinations, we incorporated a random forest conditional variable importance measure into our KPCA-based analytical approach to demonstrate the relative importance of metabolites. Using a market basket analysis, hippurate, the most important variable detected in the importance measure, was associated with high levels of some vitamins and minerals present in foods eaten the previous day, suggesting a relationship between increased hippurate and intake of a wide variety of vegetables and fruits. Therefore, the KPCA-incorporated analytical approach described herein enabled us to capture input–output responses, and should be useful not only for metabolic profiling but also for profiling in other areas of biological and environmental systems.

Environmental metabolomics with data science for investigating ecosystem homeostasis

A natural ecosystem can be viewed as the interconnections between complex metabolic reactions and environments. Humans, a part of these ecosystems, and their activities strongly affect the environments. To account for human effects within ecosystems, understanding what benefits humans receive by facilitating the maintenance of environmental homeostasis is important. This review describes recent applications of several NMR approaches to the evaluation of environmental homeostasis by metabolic profiling and data science. The basic NMR strategy used to evaluate homeostasis using big data collection is similar to that used in human health studies. Sophisticated metabolomic approaches (metabolic profiling) are widely reported in the literature. Further challenges include the analysis of complex macromolecular structures, and of the compositions and interactions of plant biomass, soil humic substances, and aqueous particulate organic matter. To support the study of these topics, we also discuss sample preparation techniques and solid-state NMR approaches. Because NMR approaches can produce a number of data with high reproducibility and inter-institution compatibility, further analysis of such data using machine learning approaches is often worthwhile. We also describe methods for data pretreatment in solid-state NMR and for environmental feature extraction from heterogeneously-measured spectroscopic data by machine learning approaches.

Knowns and unknowns in metabolomics identified by multidimensional NMR and hybrid MS/NMR methods

Metabolomics continues to make rapid progress through the development of new and better methods and their applications to gain insight into the metabolism of a wide range of different biological systems from a systems biology perspective. Customization of NMR databases and search tools allows the faster and more accurate identification of known metabolites, whereas the identification of unknowns, without a need for extensive purification, requires new strategies to integrate NMR with mass spectrometry, cheminformatics, and computational methods. For some applications, the use of covalent and non-covalent attachments in the form of labeled tags or nanoparticles can significantly reduce the complexity of these tasks.

Structure and Metabolic-Flow Analysis of Molecular Complexity in a (13) C-Labeled Tree by 2D and 3D NMR

Improved signal identification for biological small molecules (BSMs) in a mixture was demonstrated by using multidimensional NMR on samples from (13) C-enriched Rhododendron japonicum (59.5 atom%) cultivated in air containing (13) C-labeled carbon dioxide for 14 weeks. The resonance assignment of 386 carbon atoms and 380 hydrogen atoms in the mixture was achieved. 42 BSMs, including eight that were unlisted in the spectral databases, were identified. Comparisons between the experimental values and the (13) C chemical shift values calculated by density functional theory supported the identifications of unlisted BSMs. Tracing the (13) C/(12) C ratio by multidimensional NMR spectra revealed faster and slower turnover ratios of BSMs involved in central metabolism and those categorized as secondary metabolites, respectively. The identification of BSMs and subsequent flow analysis provided insight into the metabolic systems of the plant.

The Effect of Molecular Conformation on the Accuracy of Theoretical (1)H and (13)C Chemical Shifts Calculated by Ab Initio Methods for Metabolic Mixture Analysis

NMR spectroscopy is a powerful method for analyzing metabolic mixtures. The information obtained from an NMR spectrum is in the form of physical parameters, such as chemical shifts, and construction of databases for many metabolites will be useful for data interpretation. To increase the accuracy of theoretical chemical shifts for development of a database for a variety of metabolites, the effects of sets of conformations (structural ensembles) and the levels of theory on computations of theoretical chemical shifts were systematically investigated for a set of 29 small molecules in the present study. For each of the 29 compounds, 101 structures were generated by classical molecular dynamics at 298.15 K, and then theoretical chemical shifts for 164 (1)H and 123 (13)C atoms were calculated by ab initio quantum chemical methods. Six levels of theory were used by pairing Hartree-Fock, B3LYP (density functional theory), or second order Møller-Plesset perturbation with 6-31G or aug-cc-pVDZ basis set. The six average fluctuations in the (1)H chemical shift were ±0.63, ± 0.59, ± 0.70, ± 0.62, ± 0.75, and ±0.66 ppm for the structural ensembles, and the six average errors were ±0.34, ± 0.27, ± 0.32, ± 0.25, ± 0.32, and ±0.25 ppm. The results showed that chemical shift fluctuations with changes in the conformation because of molecular motion were larger than the differences between computed and experimental chemical shifts for all six levels of theory. In conclusion, selection of an appropriate structural ensemble should be performed before theoretical chemical shift calculations for development of an accurate database for a variety of metabolites.