Imbalance in Glucose Metabolism Regulates the Transition of Microglia from Homeostasis to Disease-Associated Microglia Stage 1.

Microglia undergo two-stage activation in neurodegenerative diseases, known as disease-associated microglia (DAM). TREM2 mediates the DAM2 stage transition, but what regulates the first DAM1 stage transition is unknown. We report that glucose dyshomeostasis inhibits DAM1 activation and PKM2 plays a role. As in tumors, PKM2 was aberrantly elevated in both male and female human AD brains, but unlike in tumors, it is expressed as active tetramers, as well as among TREM2+ microglia surrounding plaques in 5XFAD male and female mice. snRNAseq analyses of microglia without Pkm2 in 5XFAD mice revealed significant increases in DAM1 markers in a distinct metabolic cluster, which is enriched in genes for glucose metabolism, DAM1, and AD risk. 5XFAD mice incidentally exhibited a significant reduction in amyloid pathology without microglial Pkm2 Surprisingly, microglia in 5XFAD without Pkm2 exhibited increases in glycolysis and spare respiratory capacity, which correlated with restoration of mitochondrial cristae alterations. In addition, in situ spatial metabolomics of plaque-bearing microglia revealed an increase in respiratory activity. These results together suggest that it is not only glycolytic but also respiratory inputs that are critical to the development of DAM signatures in 5XFAD mice.

NMR 1H, 13C, 15N backbone resonance assignments of wild-type human K-Ras and its oncogenic mutants G12D and G12C bound to GTP.

Human K-Ras protein, which is a member of the GTPase Ras family, hydrolyzes GTP to GDP and concomitantly converts from its active to its inactive state. It is a key oncoprotein, because several mutations, particularly those at residue position 12, occur with a high frequency in a wide range of human cancers. The K-Ras protein is therefore an important target for developing therapeutic anti-cancer agents. In this work we report the almost complete sequence-specific resonance assignments of wild-type and the oncogenic G12C and G12D mutants in the GTP-complexed active forms, including the functionally important Switch I and Switch II regions. These assignments serve as the basis for a comprehensive functional dynamics study of wild-type K-Ras and its G12 mutants.

DEEP Picker1D and Voigt Fitter1D: a versatile tool set for the automated quantitative spectral deconvolution of complex 1D-NMR spectra.

The quantitative deconvolution of 1D-NMR spectra into individual resonances or peaks is a key step in many modern NMR workflows as it critically affects downstream analysis and interpretation. Depending on the complexity of the NMR spectrum, spectral deconvolution can be a notable challenge. Based on the recent deep neural network DEEP Picker and Voigt Fitter for 2D NMR spectral deconvolution, we present here an accurate, fully automated solution for 1D-NMR spectral analysis, including peak picking, fitting, and reconstruction. The method is demonstrated for complex 1D solution NMR spectra showing excellent performance also for spectral regions with multiple strong overlaps and a large dynamic range whose analysis is challenging for current computational methods. The new tool will help streamline 1D-NMR spectral analysis for a wide range of applications and expand their reach toward ever more complex molecular systems and their mixtures.

Excited-state observation of active K-Ras reveals differential structural dynamics of wild-type versus oncogenic G12D and G12C mutants.

Despite the prominent role of the K-Ras protein in many different types of human cancer, major gaps in atomic-level information severely limit our understanding of its functions in health and disease. Here, we report the quantitative backbone structural dynamics of K-Ras by solution nuclear magnetic resonance spectroscopy of the active state of wild-type K-Ras bound to guanosine triphosphate (GTP) nucleotide and two of its oncogenic P-loop mutants, G12D and G12C, using a new nanoparticle-assisted spin relaxation method, relaxation dispersion and chemical exchange saturation transfer experiments covering the entire range of timescales from picoseconds to milliseconds. Our combined experiments allow detection and analysis of the functionally critical Switch I and Switch II regions, which have previously remained largely unobservable by X-ray crystallography and nuclear magnetic resonance spectroscopy. Our data reveal cooperative transitions of K-Ras·GTP to a highly dynamic excited state that closely resembles the partially disordered K-Ras·GDP state. These results advance our understanding of differential GTPase activities and signaling properties of the wild type versus mutants and may thus guide new strategies for the development of therapeutics.

Differential metabolism between biofilm and suspended Pseudomonas aeruginosa cultures in bovine synovial fluid by 2D NMR-based metabolomics.

Total joint arthroplasty is a common surgical procedure resulting in improved quality of life; however, a leading cause of surgery failure is infection. Periprosthetic joint infections often involve biofilms, making treatment challenging. The metabolic state of pathogens in the joint space and mechanism of their tolerance to antibiotics and host defenses are not well understood. Thus, there is a critical need for increased understanding of the physiological state of pathogens in the joint space for development of improved treatment strategies toward better patient outcomes. Here, we present a quantitative, untargeted NMR-based metabolomics strategy for Pseudomonas aeruginosa suspended culture and biofilm phenotypes grown in bovine synovial fluid as a model system. Significant differences in metabolic pathways were found between the suspended culture and biofilm phenotypes including creatine, glutathione, alanine, and choline metabolism and the tricarboxylic acid cycle. We also identified 21 unique metabolites with the presence of P. aeruginosa in synovial fluid and one uniquely present with the biofilm phenotype in synovial fluid. If translatable in vivo, these unique metabolite and pathway differences have the potential for further development to serve as targets for P. aeruginosa and biofilm control in synovial fluid.

COLMARq: A Web Server for 2D NMR Peak Picking and Quantitative Comparative Analysis of Cohorts of Metabolomics Samples.

Highly quantitative metabolomics studies of complex biological mixtures are facilitated by the resolution enhancement afforded by 2D NMR spectra such as 2D 13C-1H HSQC spectra. Here, we describe a new public web server, COLMARq, for the semi-automated analysis of sets of 2D HSQC spectra of cohorts of samples. The workflow of COLMARq includes automated peak picking using the deep neural network DEEP Picker, quantitative cross-peak volume extraction by numerical fitting using Voigt Fitter, the matching of corresponding cross-peaks across cohorts of spectra, peak volume normalization between different spectra, database query for metabolite identification, and basic univariate and multivariate statistical analyses of the results. COLMARq allows the analysis of cross-peaks that belong to both known and unknown metabolites. After a user has uploaded cohorts of 2D 13C-1H HSQC and optionally 2D 1H-1H TOCSY spectra in their preferred format, all subsequent steps on the web server can be performed fully automatically, allowing manual editing if needed and the sessions can be saved for later use. The accuracy, versatility, and interactive nature of COLMARq enables quantitative metabolomics analysis, including biomarker identification, of a broad range of complex biological mixtures as is illustrated for cohorts of samples from bacterial cultures of Pseudomonas aeruginosa in both its biofilm and planktonic states.

Fundamental and practical aspects of machine learning for the peak picking of biomolecular NMR spectra.

Rapid progress in machine learning offers new opportunities for the automated analysis of multidimensional NMR spectra ranging from protein NMR to metabolomics applications. Most recently, it has been demonstrated how deep neural networks (DNN) designed for spectral peak picking are capable of deconvoluting highly crowded NMR spectra rivaling the facilities of human experts. Superior DNN-based peak picking is one of a series of critical steps during NMR spectral processing, analysis, and interpretation where machine learning is expected to have a major impact. In this perspective, we lay out some of the unique strengths as well as challenges of machine learning approaches in this new era of automated NMR spectral analysis. Such a discussion seems timely and should help define common goals for the NMR community, the sharing of software tools, standardization of protocols, and calibrate expectations. It will also help prepare for an NMR future where machine learning and artificial intelligence tools will be common place.

Cadaverine Is a Switch in the Lysine Degradation Pathway in Pseudomonas aeruginosa Biofilm Identified by Untargeted Metabolomics.

There is a critical need to accurately diagnose, prevent, and treat biofilms in humans. The biofilm forming P. aeruginosa bacteria can cause acute and chronic infections, which are difficult to treat due to their ability to evade host defenses along with an inherent antibiotic-tolerance. Using an untargeted NMR-based metabolomics approach, we identified statistically significant differences in 52 metabolites between P. aeruginosa grown in the planktonic and lawn biofilm states. Among them, the metabolites of the cadaverine branch of the lysine degradation pathway were systematically decreased in biofilm. Exogenous supplementation of cadaverine caused significantly increased planktonic growth, decreased biofilm accumulation by 49% and led to altered biofilm morphology, converting to a pellicle biofilm at the air-liquid interface. Our findings show how metabolic pathway differences directly affect the growth mode in P. aeruginosa and could support interventional strategies to control biofilm formation.

DEEP picker is a deep neural network for accurate deconvolution of complex two-dimensional NMR spectra.

The analysis of nuclear magnetic resonance (NMR) spectra for the comprehensive and unambiguous identification and characterization of peaks is a difficult, but critically important step in all NMR analyses of complex biological molecular systems. Here, we introduce DEEP Picker, a deep neural network (DNN)-based approach for peak picking and spectral deconvolution which semi-automates the analysis of two-dimensional NMR spectra. DEEP Picker includes 8 hidden convolutional layers and was trained on a large number of synthetic spectra of known composition with variable degrees of crowdedness. We show that our method is able to correctly identify overlapping peaks, including ones that are challenging for expert spectroscopists and existing computational methods alike. We demonstrate the utility of DEEP Picker on NMR spectra of folded and intrinsically disordered proteins as well as a complex metabolomics mixture, and show how it provides access to valuable NMR information. DEEP Picker should facilitate the semi-automation and standardization of protocols for better consistency and sharing of results within the scientific community.

Observation of Sub-Microsecond Protein Methyl-Side Chain Dynamics by Nanoparticle-Assisted NMR Spin Relaxation.

Amino-acid side-chain properties in proteins are key determinants of protein function. NMR spin relaxation of side chains is an important source of information about local protein dynamics and flexibility. However, traditional solution NMR relaxation methods are most sensitive to sub-nanosecond dynamics lacking information on slower ns-µs time-scale motions. Nanoparticle-assisted NMR spin relaxation (NASR) of methyl-side chains is introduced here as a window into these ns-µs dynamics. NASR utilizes the transient and nonspecific interactions between folded proteins and slowly tumbling spherical nanoparticles (NPs), whereby the increase of the relaxation rates reflects motions on time scales from ps all the way to the overall tumbling correlation time of the NPs ranging from hundreds of ns to µs. The observed motional amplitude of each methyl group can then be expressed by a model-free NASR S2 order parameter. The method is demonstrated for 2H-relaxation of CH2D methyl moieties and cross-correlated relaxation of CH3 groups for proteins Im7 and ubiquitin in the presence of anionic silica-nanoparticles. Both types of relaxation experiments, dominated by either quadrupolar or dipolar interactions, yield highly consistent results. Im7 shows additional dynamics on the intermediate time scales taking place in a functionally important loop, whereas ubiquitin visits the majority of its conformational substates on the sub-ns time scale. These experimental observations are in good agreement with 4-10 µs all-atom molecular dynamics trajectories. NASR probes side-chain dynamics on a much wider range of motional time scales than previously possible, thereby providing new insights into the interplay between protein structure, dynamics, and molecular interactions that govern protein function.

2D NMR-Based Metabolomics with HSQC/TOCSY NOAH Supersequences.

Sensitivity-improved versions of two-dimensional (2D) 13C-1H HSQC (heteronuclear single quantum coherence) and HSQC-TOCSY (HSQC-total correlation spectroscopy) NMR experiments optimized for small biological molecules and their complex mixtures encountered in metabolomics are presented that preserve the magnetization of 1H spins not directly attached to 13C spins. This allows (i) the application of rapid acquisition techniques to substantially shorten measurement time and (ii) their incorporation into supersequences (NOAH-NMR by ordered acquisition using 1H detection) for the compact acquisition of multiple 2D NMR data sets with significant gains in sensitivity, resolution, and/or time. The new pulse sequences, which are demonstrated for both metabolite model mixtures and mouse urine, offer an attractive approach for the efficient measurement of multiple 2D NMR spectra (HSQCsi and/or HSQCsi-TOCSY and TOCSY) of metabolomics samples in a single experiment for the accurate and comprehensive identification and quantitation of metabolites. These new methods bring to bear the advantages of 2D NMR to metabolomics studies with larger cohorts of samples.

Broadband Dynamics of Ubiquitin by Anionic and Cationic Nanoparticle Assisted NMR Spin Relaxation.

The quantitative and comprehensive description of the internal dynamics of proteins is critical for understanding their function. Nanoparticle-assisted 15 N NMR spin relaxation spectroscopy is a new method for the observation of picosecond to microsecond dynamics of proteins when transiently interacting with the surface of the nanoparticles (NPs). The method is applied here to the protein ubiquitin in the presence of anionic and cationic silica NPs (SNPs) of different sizes. The backbone dynamics profiles are reproducible and strikingly similar to each other, indicating that specific protein-SNP interactions are unimportant. The dynamics profiles closely match the sub-nanosecond dynamics S2 values observed by model-free analysis of standard 15 N relaxation of ubiquitin in free solution, indicating that the bulk of the ubiquitin backbone dynamics in solution is confined to sub-nanosecond timescales and, hence, it is dynamically more restrained than previous NMR studies have suggested.

COLMAR Lipids Web Server and Ultrahigh-Resolution Methods for Two-Dimensional Nuclear Magnetic Resonance- and Mass Spectrometry-Based Lipidomics.

Accurate identification of lipids in biological samples is a key step in lipidomics studies. Multidimensional nuclear magnetic resonance (NMR) spectroscopy is a powerful analytical tool for this purpose as it provides comprehensive structural information on lipid composition at atomic resolution. However, the interpretation of NMR spectra of complex lipid mixtures is currently hampered by limited spectral resolution and the absence of a customized lipid NMR database along with user-friendly spectral analysis tools. We introduce a new two-dimensional (2D) NMR metabolite database "COLMAR Lipids" that was specifically curated for hydrophobic metabolites presently containing 501 compounds with accurate experimental 2D 13C-1H heteronuclear single quantum coherence (HSQC) chemical shift data measured in CDCl3. A new module in the public COLMAR suite of NMR web servers was developed for the (semi)automated analysis of complex lipidomics mixtures (http://spin.ccic.osu.edu/index.php/colmarm/index2). To obtain 2D HSQC spectra with the necessary high spectral resolution along both 13C and 1H dimensions, nonuniform sampling in combination with pure shift spectroscopy was applied allowing the extraction of an abundance of unique cross-peaks belonging to hydrophobic compounds in complex lipidomics mixtures. As shown here, this information is critical for the unambiguous identification of underlying lipid molecules by means of the new COLMAR Lipids web server, also in combination with mass spectrometry, as is demonstrated for Caco-2 cell and lung tissue cell extracts.

Accurate and Efficient Determination of Unknown Metabolites in Metabolomics by NMR-Based Molecular Motif Identification.

Knowledge of the chemical identity of metabolite molecules is critical for the understanding of the complex biological systems to which they belong. Since metabolite identities and their concentrations are often directly linked to the phenotype, such information can be used to map biochemical pathways and understand their role in health and disease. A very large number of metabolites however are still unknown; i.e., their spectroscopic signatures do not match those in existing databases, suggesting unknown molecule identification is both imperative and challenging. Although metabolites are structurally highly diverse, the majority shares a rather limited number of structural motifs, which are defined by sets of 1H and 13C chemical shifts of the same spin system. This allows one to characterize unknown metabolites by a divide-and-conquer strategy that identifies their structural motifs first. Here, we present the structural motif-based approach "SUMMIT Motif" for the de novo identification of unknown molecular structures in complex mixtures, without the need for extensive purification, using NMR in tandem with two newly curated NMR molecular structural motif metabolomics databases (MSMMDBs). For the identification of structural motif(s), first, the 1H and 13C chemical shifts of all the individual spin systems are extracted from 2D and 3D NMR spectra of the complex mixture. Next, the molecular structural motifs are identified by querying these chemical shifts against the new MSMMDBs. One database, COLMAR MSMMDB, was derived from experimental NMR chemical shifts of known metabolites taken from the COLMAR metabolomics database, while the other MSMMDB, pNMR MSMMDB, is based on predicted chemical shifts of metabolites of several existing large metabolomics databases. For molecules consisting of multiple spin systems, spin systems are connected via long-range scalar J-couplings. When this motif-based identification method was applied to the hydrophilic extract of mouse bile fluid, two unknown metabolites could be successfully identified. This approach is both accurate and efficient for the identification of unknown metabolites and hence enables the discovery of new biochemical processes and potential biomarkers.

Extreme Nonuniform Sampling for Protein NMR Dynamics Studies in Minimal Time.

NMR spectroscopy is an extraordinarily rich source of quantitative dynamics of proteins in solution using spin relaxation or chemical exchange saturation transfer (CEST) experiments. However, 15N-CEST measurements require prolonged multidimensional, so-called pseudo-3D HSQC experiments where the pseudo dimension is a radio frequency offset Δω of a weak 15N saturation field. Nonuniform sampling (NUS) approaches have the potential to significantly speed up these measurements, but they also carry the risk of introducing serious artifacts and the systematic optimization of nonuniform sampling schedules has remained elusive. It is demonstrated here how this challenge can be addressed by using fitted cross-peaks of a reference 2D HSQC experiment as footprints, which are subsequently used to reconstruct cross-peak amplitudes of a pseudo-3D data set as a function of Δω by a linear least-squares fit. It is shown for protein Im7 how the approach can yield highly accurate CEST profiles based on an absolutely minimally sampled (AMS) data set allowing a speed-up of a factor 20-30. Spectrum-specific optimized nonuniform sampling (SONUS) schemes based on the Cramer-Rao lower bound metric were critical to achieve such a performance, revealing also more general properties of optimal sampling schedules. This is the first systematic exploration and optimization of NUS schedules for the dramatic speed-up of quantitative multidimensional NMR measurements that minimize unwanted errors.

Functional protein dynamics on uncharted time scales detected by nanoparticle-assisted NMR spin relaxation.

Protein function depends critically on intrinsic internal dynamics, which is manifested in distinct ways, such as loop motions that regulate protein recognition and catalysis. Under physiological conditions, dynamic processes occur on a wide range of time scales from subpicoseconds to seconds. Commonly used NMR spin relaxation in solution provides valuable information on very fast and slow motions but is insensitive to the intermediate nanosecond to microsecond range that exceeds the protein tumbling correlation time. Presently, very little is known about the nature and functional role of these motions. It is demonstrated here how transverse spin relaxation becomes exquisitely sensitive to these motions at atomic resolution when studying proteins in the presence of nanoparticles. Application of this novel cross-disciplinary approach reveals large-scale dynamics of loops involved in functionally critical protein-protein interactions and protein-calcium ion recognition that were previously unobservable.

Metabolic Regulation of Glycolysis and AMP Activated Protein Kinase Pathways during Black Raspberry-Mediated Oral Cancer Chemoprevention.

Oral cancer is a public health problem with an incidence of almost 50,000 and a mortality of 10,000 each year in the USA alone. Black raspberries (BRBs) have been shown to inhibit oral carcinogenesis in several preclinical models, but our understanding of how BRB phytochemicals affect the metabolic pathways during oral carcinogenesis remains incomplete. We used a well-established rat oral cancer model to determine potential metabolic pathways impacted by BRBs during oral carcinogenesis. F344 rats were exposed to the oral carcinogen 4-nitroquinoline-1-oxide in drinking water for 14 weeks, then regular drinking water for six weeks. Carcinogen exposed rats were fed a 5% or 10% BRB supplemented diet or control diet for six weeks after carcinogen exposure. RNA-Seq transcriptome analysis on rat tongue, and mass spectrometry and NMR metabolomics analysis on rat urine were performed. We tentatively identified 57 differentially or uniquely expressed metabolites and over 662 modulated genes in rats being fed with BRB. Glycolysis and AMPK pathways were modulated during BRB-mediated oral cancer chemoprevention. Glycolytic enzymes Aldoa, Hk2, Tpi1, Pgam2, Pfkl, and Pkm2 as well as the PKA-AMPK pathway genes Prkaa2, Pde4a, Pde10a, Ywhag, and Crebbp were downregulated by BRBs during oral cancer chemoprevention. Furthermore, the glycolysis metabolite glucose-6-phosphate decreased in BRB-administered rats. Our data reveal the novel metabolic pathways modulated by BRB phytochemicals that can be targeted during the chemoprevention of oral cancer.

Real-Time Pure Shift HSQC NMR for Untargeted Metabolomics.

Sensitivity and resolution are key considerations for NMR applications in general and for metabolomics in particular, where complex mixtures containing hundreds of metabolites over a large range of concentrations are commonly encountered. There is a strong demand for advanced methods that can provide maximal information in the shortest possible time frame. Here, we present the optimization and application of the recently introduced 2D real-time BIRD 1H-13C HSQC experiment for NMR-based metabolomics of aqueous samples at 13C natural abundance. For mouse urine samples, it is demonstrated how this real-time pure shift sensitivity-improved heteronuclear single quantum correlation method provides broadband homonuclear decoupling along the proton detection dimension and thereby significantly improves spectral resolution in regions that are affected by spectral overlap. Moreover, the collapse of the scalar multiplet structure of cross-peaks leads to a sensitivity gain of about 40-50% over a traditional 2D HSQC-SI experiment. The experiment works well over a range of magnetic field strengths and is particularly useful when resonance overlap in crowded regions of the HSQC spectra hampers accurate metabolite identification and quantitation.

Identification of Unknown Metabolomics Mixture Compounds by Combining NMR, MS, and Cheminformatics.

Metabolomics aims at the comprehensive identification of metabolites in complex mixtures to characterize the state of a biological system and elucidate their roles in biochemical pathways. For many biological samples, a large number of spectral features observed by NMR spectroscopy and mass spectrometry (MS) belong to unknowns, i.e., these features do not belong to metabolites that have been previously identified, and their spectral information is not available in databases. By combining NMR, MS, and combinatorial cheminformatics, the analysis of unknowns can be pursued in complex mixtures requiring minimal purification. This chapter describes the SUMMIT MS/NMR approach covering sample preparation, NMR and MS data collection and processing, and the identification of likely unknowns with the use of cheminformatics tools and the prediction of NMR spectral properties.

Time-Resolved Protein Side-Chain Motions Unraveled by High-Resolution Relaxometry and Molecular Dynamics Simulations.

Motions of proteins are essential for the performance of their functions. Aliphatic protein side chains and their motions play critical roles in protein interactions: for recognition and binding of partner molecules at the surface or serving as an entropy reservoir within the hydrophobic core. Here, we present a new NMR method based on high-resolution relaxometry and high-field relaxation to determine quantitatively both motional amplitudes and time scales of methyl-bearing side chains in the picosecond-to-nanosecond range. We detect a wide variety of motions in isoleucine side chains in the protein ubiquitin. We unambiguously identify slow motions in the low nanosecond range, which, in conjunction with molecular dynamics computer simulations, could be assigned to transitions between rotamers. Our approach provides unmatched detailed insight into the motions of aliphatic side chains in proteins and provides a better understanding of the nature and functional role of protein side-chain motions.